Differential gene expression and pathway induction in reconstituted human vaginal epithelium following C. albicans challenge.
The reconstituted human vaginal epithelium (RVE) model is an excellent
in vitro surrogate to study epithelium-specific responses of vaginal candidiasis, as the tissue layer is sufficiently differentiated and supports robust hyphal invasion, and infected RVE tissue largely resembles
in vivo infection dynamics (
13,
14). In order to elucidate global host transcriptomic changes in the vaginal epithelium in response to challenge with
C. albicans (compared to phosphate-buffered saline [PBS] mock control), total epithelial RNA was selectively isolated from three independent RVE at 6 and 24 h postchallenge and subjected to microarray analysis. As with oral epithelium, the intermediate (6-h) time point is associated with initial fungal adherence and microbial recognition, while the late (24-h) time point is associated with fungal invasion and cellular damage (
10,
15). Approximately 800 and nearly 4,000 genes were differentially expressed (
P < 0.001) at 6 and 24 h, respectively, in response to
C. albicans (
Fig. 1). Comparatively few genes were regulated in response to PBS mock treatment at the same time points (see Fig. S1 in the supplemental material). At the intermediate stage of infection (6 h postinfection), the majority of differentially expressed genes were upregulated (
Fig. 1A), with only 65 genes strongly upregulated (>4-fold) and none showing strong downregulation (>4-fold). However, by late stages of infection (24 h), an increase in the proportion of genes showing downregulated expression was observed (
Fig. 1B). Approximately 320 genes were strongly upregulated at 24 h (>4-fold), and over half of the genes showing upregulation at 6 h were also strongly upregulated at the later time point (
Fig. 1C). Surprisingly, relatively few genes were strongly (>4-fold) downregulated in response to fungal challenge at either time point.
Gene Ontology, pathway, and network mapping revealed profiles from
C. albicans-infected cells to be consistent with MAPK, NF-κB, phosphatidylinositol 3-kinase (PI3K), ErbB receptor, and tumor necrosis factor (TNF) signaling pathways (
Fig. 2). Pathways involving extracellular matrix remodeling, including proteoglycans in cancer, focal adhesion, adherens junctions, and tight junctions, were also significantly enriched during
C. albicans infection. Pathways involved in responses to infection by other microbes, including Epstein-Barr virus,
Shigella, hepatitis B virus, influenza A virus, herpesvirus,
Salmonella, and trypanosomes, were also predicted to be activated, suggesting conservation of epithelial responses to a broad array of pathogens. Pathways predicted to be activated were generally conserved at the 6-h and 24-h time points. A list of individually expressed genes may be found in Table S1 in the supplemental material.
Genes involved in innate inflammatory signaling were strongly induced by C. albicans, including the cytokine genes IL-8 (100-fold), IL-1A (18-fold), IL-1B (3.8-fold), CXCL1 (19-fold), CXCL2 (26-fold), GM-CSF (10-fold), and PTGS2 (prostaglandin synthase) (7.3-fold), many of which play critical roles in recruiting inflammatory cells (particularly neutrophils) to the site of infection. Similar to previous findings, there was clear induction of genes associated with MAPK activity: MAP3K2 (6.8-fold), MAP2K3 (4-fold), MAP3K9 (4-fold), and MAP4K4 (2.7-fold). Additionally, C. albicans infection led to epithelial induction of c-FOS (32-fold) and c-JUN (17.7-fold), which encode members of two families that form the heterodimeric transcription factor AP1, a major effector of MAPK activation. The dual-specificity phosphatase 1 (DUSP1) gene, encoding a regulator of MAPK signaling, was also elevated (6.7-fold) in response to C. albicans.
A number of genes involved in tissue repair, wound healing, or dampening of active inflammation were also upregulated during
C. albicans infection, including the genes coding for IL-24 (2.3-fold) and the IL-1 receptor antagonist (IL-1RN) (4-fold) (
16,
17). Interestingly, a number of other related genes were also induced, including genes coding for heparin-binding EGF-like growth factor (HBEGF) (39.5-fold) and epiregulin (EREG) (6-fold), which are members of the epidermal growth factors (EGFs). They exert their functions by binding to their cognate receptors (EGFR) or the v-
erb-b2 oncogene homolog (ERBB) to induce cellular proliferation and healing of skin and epidermal tissues (
18,
19).
Candidalysin damages and activates vaginal epithelial cells.
As we observed upregulated expression of genes encoding several proinflammatory cytokines (e.g.,
IL-1A,
IL-1B,
IL-8, and
GM-CSF) and chemokines during RVE challenge with
C. albicans at time points when hyphae invaded the vaginal tissue, we sought to determine whether the hypha-associated peptide toxin Candidalysin similarly elicited these effector and damage responses. Indeed, there was a dose-dependent release of LDH when Candidalysin was applied to A431 cells (
Fig. 3A). Significant levels of cellular damage were observed with doses above 15 μM compared to treatment with the vehicle control.
Vaginal epithelial cells respond to
C. albicans hyphae by activating the p38-MAPK and ERK1/2-MAPK signaling pathways, resulting in the regulated secretion of proinflammatory cytokines (
8). To assess whether Candidalysin is capable of activating these pathways, epithelial cells were exposed to Candidalysin
in vitro, and c-Fos production/MKP1 phosphorylation was assessed by Western blotting (
Fig. 3B). The c-Fos/p-MKP1 response was induced strongly by 15 and 70 μM Candidalysin, whereas the vehicle was unable to activate signaling. Concomitant with damage, treatment with Candidalysin caused a dose-dependent increase in the release of IL-1α, IL-1β, granulocyte colony-stimulating factor (G-CSF), granulocyte-macrophage colony-stimulating factor (GM-CSF), and IL-8 in spent culture supernatants (
Fig. 3C to
H). The lone exception was IL-6, which was significantly elevated only at the highest Candidalysin concentration (70 μM). With the exception of IL-6, all the cytokines assayed were significantly induced at Candidalysin doses above 3 μM; however, this dose was insufficient to cause significant damage (
Fig. 3A), suggesting that Candidalysin exhibits dual functionality, playing both immunostimulatory and lytic roles against vaginal epithelial cells, similar to what is observed in oral epithelia (
12).
Candidalysin is required for vaginitis immunopathology.
We next questioned whether Ece1p and/or Candidalysin contributes to immunopathology in an established estrogen-dependent mouse model of VVC. Therefore, we utilized strains of
C. albicans that had had both copies of
ECE1 deleted (
ece1Δ/Δ) and restored with one full-length allele (
ece1Δ/Δ +
ECE1) or one mutant allele lacking the Candidalysin-encoding region of
ECE1 (
ece1Δ/Δ + ECE1Δ184–279), along with the appropriate parental isogenic control (BWP17 + CIp30, referred to as here the wild type [WT]). Somewhat surprisingly, recovered fungal burdens from the vaginal lavage fluid were not significantly different between strains at either day 3 (
Fig. 4A) or day 7 (
Fig. 4B) postinoculation (p.i.). However, there was a significant reduction in the number of neutrophils recruited into the vaginal lumen during challenge with either the
ece1Δ/Δ or
ece1Δ/Δ + ECE1Δ184–279 strain, which was restored to WT levels during infection with the
ece1Δ/Δ +
ECE1 reintegrant strain (
Fig. 4C,
D, and
G, yellow arrows). Consistent with this phenotype, levels of the damage biomarker LDH were significantly reduced with the same mutants compared to infection with the WT or the
ece1Δ/Δ +
ECE1 reintegrant (
Fig. 4E and
F). Given our previous data using hypha-deficient strains, a morphogenesis defect may account for this phenotype (
6). However,
ece1Δ/Δ and
ece1Δ/Δ + ECE1Δ184–279 strains robustly formed hyphae at the vaginal mucosa, as did WT and
ece1Δ/Δ +
ECE1 strains (
Fig. 4G, green arrows). Thus, these results demonstrate that Candidalysin is required for vaginal immunopathogenesis
in vivo and that hypha formation alone is insufficient to elicit the hallmark immunopathology.
Candidalysin-dependent innate cytokine expression is conserved between mouse and human.
We also wanted to determine whether the Candidalysin-induced innate immune response observed in human vaginal epithelial cells paralleled cytokine expression in the murine vaginal mucosa
in vivo. RNA was isolated from whole vaginas of mice challenged with WT,
ece1Δ/Δ,
ece1Δ/Δ +
ECE1,
ece1Δ/Δ + ECE1Δ184–279, and PBS, and gene expression was assessed by quantitative PCR (qPCR). Overall, cytokine gene expression patterns were similar between
in vitro and
in vivo samples, including Candidalysin-induced expression of the genes
Il-6,
Cxcl2,
Il-1a, and
Il-1b (
Fig. 5A and
C to
E). There was a similar trend for expression of the genes
Cxcl1 and
Gm-csf, although only the
ECE1-null mutant (
ece1Δ/Δ) demonstrated a statistically significant reduction in cytokine gene induction (
Fig. 5B and
F). Unexpectedly,
G-csf gene expression was not increased during challenge with any of the fungal strains, unlike that observed with Candidalysin treatment (
Fig. 5G). In the oral cavity,
C. albicans induces expression of the antimicrobial peptide (AMP) cathelicidin, the murine equivalent of which is the cathelicidin-related AMP (CAMP) (
20). Interestingly, the gene encoding CAMP was not induced in the vagina by Candidalysin and in fact was downregulated similarly by all the strains compared to mock treatment (
Fig. 5H). However, induction of the antimicrobial peptide murine β-defensin 3 (mBD3) gene was Candidalysin dependent (
Fig. 5I).
We also sought to determine if two inflammatory markers previously identified as associated with VVC immunopathology were regulated in a Candidalysin-dependent manner. Expression of the gene coding for S100A8 (
S100a8), a calcium-binding protein with important functions in antifungal defense and danger responses and strongly induced during
C. albicans infection, was almost completely absent during infection with Candidalysin deletion strains (
Fig. 5J) (
21,
22). Similarly, the gene encoding serum amyloid A3 (
Saa3), an inducible acute-phase apolipoprotein capable of recruiting immune cells to inflammatory sites, was similarly increased in a Candidalysin-dependent fashion (
Fig. 5K) (
23,
24).
Finally, we investigated whether production of cytokines at the protein level (at both day 3 and day 7 p.i.) was Candidalysin dependent. Indeed,
C. albicans-mediated secretion of IL-1α, IL-1β, CXCL2, and S100A8 into the vaginal lavage fluid required expression of functional Candidalysin (
Fig. 6A to
H). Despite increased expression of the genes encoding IL-6, CXCL1, and GM-CSF in vaginal tissue, we were unable to demonstrably quantify these cytokines at the protein level in the lavage fluid of mice inoculated with any of the
C. albicans strains tested (data not shown).