1,990
Views
4
CrossRef citations to date
0
Altmetric
Research Paper

Set1-mediated H3K4 methylation is required for Candida albicans virulence by regulating intracellular level of reactive oxygen species

, , , ORCID Icon & ORCID Icon
Pages 2648-2658 | Received 14 May 2021, Accepted 10 Sep 2021, Published online: 26 Oct 2021

ABSTRACT

Candida albicans is an opportunistic human fungal pathogen that exists in normal flora but can cause infection in immunocompromised individuals. The transition to pathogenic C. albicans requires a change of various gene expressions. Because histone-modifying enzymes can regulate gene expression, they are thought to control the virulence of C. albicans. Indeed, the absence of H3 lysine 4 (H3K4) methyltransferase Set1 has been shown to reduce the virulence of C. albicans; however, Set1-regulated genes responsible for this attenuated virulence phenotype remain unknown. Here, we demonstrated that Set1 positively regulates the expression of mitochondrial protein genes by methylating H3K4. In particular, levels of cellular mitochondrial reactive oxygen species (ROS) were higher in Δset1 than in the wild-type due to the defect of those genes’ expression. Set1 deletion also increases H2O2 sensitivity and prevents proper colony formation when interacting with macrophage in vitro, consistent with its attenuated virulence in vivo. Together, these findings suggest that Set1 is required to regulate proper cellular ROS production by positively regulating the expression of mitochondrial protein genes and subsequently sustaining mitochondrial membrane integrity. Consequently, C. albicans maintains proper ROS levels via Set1-mediated transcriptional regulation, thus establishing a rapid defense against external ROS generated by the host.

Introduction

Candida albicans is the most common fungal pathogen in humans. Although it is a member of the normal flora of the skin and mucosal surfaces in most healthy people, C. albicans can cause opportunistic infection in response to environmental changes in the host [Citation1]. Indeed, C. albicans overgrowth can cause local or systemic infections, from candidiasis to life-threatening candidemia; therefore, it is necessary to study its pathogenesis to more effectively prevent and treat C. albicans infections [Citation2].

When transitioning from a non-virulent commensal to a virulent pathogen, C. albicans undergoes various phenotypic switching events, including hyphal formation [Citation1,Citation3], white-opaque switching [Citation4], biofilm formation [Citation5], and secreted aspartyl proteinase expression [Citation2,Citation6–8]. In addition, C. albicans experiences a considerable change in the expression of genes encoding virulence factors, which are involved in recognizing the host environment [Citation9]. Consequently, this leads to phenotypic switching and survival within the host, allowing C. albicans to proliferate rapidly by neutralizing and resisting host attacks [Citation2,Citation6,Citation10]. Therefore, it is important to understand the differentially expressed genes and transcriptional regulation mechanisms that alter virulence-related genes’ expression.

One of the mechanisms that regulate transcription is the post-translational modification of histone proteins, which is a well-conserved phenomenon in all eukaryotes. Histone protein residues can be chemically modified in various ways, including acetylation, phosphorylation, methylation, and ubiquitination [Citation11,Citation12]. The methylation of histone H3 at lysine 4 (H3K4) is a well-conserved and widely studied histone modification due to its positive role in transcription [Citation13,Citation14]. In general, histone lysine residues can be mono-, di-, or trimethylated, and H3K4 trimethylation (H3K4me3) is a marker of active transcription given that it is enriched in the 5′ regions of actively transcribed genes [Citation13,Citation14].

A previous study revealed that the deletion of Set1, the only H3K4 methyltransferase in C. albicans, attenuated its virulence in mice [Citation15]. H3K4 methyltransferase-deletion in other pathogenic fungi has also been shown to attenuate the virulence [Citation16–21]. H3K4 methyltransferases positively regulate the expression of secondary metabolite genes in some fungal pathogens, thereby contributing to their virulence [Citation17,Citation20]; however, the mechanisms linking the effect of Set1 on virulence to transcriptional regulation by Set1-mediated H3K4 methylation in C. albicans remain unclear.

In this study, we investigated the effect of Set1 on the pathogenesis of C. albicans based on the finding that Set1-regulated genes are responsible for the pathogenicity of C. albicans. Whole-transcriptome sequencing (RNA-seq) revealed that the decreased expression of mitochondrial genes in Set1-deleted mutants increased cellular reactive oxygen species (ROS). Therefore, we suggest that C. albicans maintains proper ROS levels via Set1-mediated transcriptional regulation, thus establishing a rapid defense against external ROS generated by the host. Therefore, we suggest that Set1-mediated gene expression enables C. albicans to respond more rapidly to ROS generated by the host and protect against it.

Materials and methods

Strains and media

The SET1 deleted strain (Δset1) was given from Dr Clancy group [Citation15]. C. albicans strains were grown in YPD media (1% yeast extracts, 2% peptone, 2% glucose) at 30°C. The Raw264.7 mouse macrophage cell line was cultured in DMEM (Cellgro, 10–013-CV) with 10% FBS (Gibco, 26,140–079) and 1% penicillin/streptomycin (Gibco 15,140–122), in an atmosphere of 5% CO2 and 95% humidity at 37°C. For oxidative stress assay, overnight cultured cells were 10-fold serially diluted from starting OD600 of 0.5, and each 3 μl was spotted to YPD plate, including 5 mM H2O2 (Sigma, 216,763). The colony was observed after 2 days of culture at 37°C.

RNA-seq analysis

Total RNA was extracted using NucleoSpin® RNA (MN, MN740955) according to the protocol of the manufacturer using each duplicated sample. C. albicans were grown in YPD and harvested at exponential phase (OD600 = 1.0). For sequencing, mRNA was captured using NEBNext® Poly(A) mRNA Magnetic Isolation Module (NEB, E7490), and a strand-specific sequencing library was synthesized using NEBNext® UltraTM Directional RNA Library Prep Kit for Illumina (NEB, E7420) according to the instruction manual. To compare the differential expression between wild-type (WT) and Δset1, we used DEseq2 normalization. Expression data was visualized using Heatmap generated by Pheatmap R package and Integrative Genomics Viewer (IGV) genome browser track.

Chromatin immunoprecipitation (ChIP)

The ChIP assay was performed, as previously described [Citation22]. Antibody used in the ChIP assay was anti-H3K4me3. ChIP DNAs were analyzed by quantitative PCR (qPCR) using the SYBR Green PCR mix (Toyobo, TOQPS-201) and the Applied Biosystems 7500 Real-Time PCR System. The sequences of primers used in this study are listed in Table S1.

Cellular ROS observation

Overnight cultured cells were diluted to OD600 of 0.5 into fresh 5 ml YPD. After 1 h incubation at 30°C with shaking, cells were harvested and washed with PBS. For making spheroplast, cells were resuspended with zymolyase buffer (1 M sorbitol, 50 mM Tris-Cl, pH 7.4), added 10 μg zymolyase, and incubated at 30°C for 10 min. To observe cellular ROS, cells were incubated with 5 mM CellROX® Green (Thermo, C10444) at 37°C for 30 min, and sequentially DAPI (Sigma, D9542) was added and incubated for 10 min. After staining, cells were washed with PBS three times and observed under a fluorescence microscope.

Macrophage interaction assay

For the survivability against macrophage attacks, harvested C. albicans cells were resuspended at 107 cells ml−1 in cold PBS containing 10% FBS. RAW264.7 cells were seeded in each well of 96-well culture dish at 2.5 × 104 cells 150 μl−1 per well or 5 × 104 cells 150 μl−1 per well. Prepared C. albicans cells were serially diluted, and 50 μl cells were co-incubated with macrophages. The plate was incubated on ice for 30 min and followed by cultured for 24 h at 37°C in 5% CO2. After 24 h incubation, C. albicans colony was observed and counted to calculate survived cells.

Mouse survival test

Animal experiments were performed at the Kangwon National University Animal Laboratory Center with approval through the Institutional Animal Care and Use Committee (IACUC) of Kangwon University (Approval Number KW-170302-5). Five-week-old female BALB/c mice were acclimated for 1 week, and tail vein injected with 106 CFU of C. albicans. In 25 days after injection, survived Δset1-injected mice were euthanized with CO2 inhalation. Kidney, spleen, and liver sections were stained with H&E (Hematoxylin and Eosin) or PAS (Periodic Acid Schiff) and observed under a microscope.

Results

Set1 is required for full virulence but has a marginal effect on overall gene expression

Set1 is the only H3K4 methyltransferase in C. albicans and is necessary for its full virulence in ICR mice [Citation15]; however, the relationship between Set1-mediated H3K4 methylation and virulence remains unknown. Therefore, we investigated the mechanism by which Set1 regulates the pathogenicity of C. albicans through H3K4 methylation. First, we examined whether the virulence of Δset1 was similar in BALB/c mice, which are inbred mice compared to outbred ICR mice. Briefly, two groups of mice (n = 5 per group) were infected with 1 × 106 CFU of wild-type (WT) and Δset1 C. albicans strains via tail vein injection, and their survival was observed for 21 days. We observed that BALB/c mice infected with Δset1 survived longer than those infected with the WT (), indicating that the Δset1 mutant displays attenuated virulence in this mouse model. Because the Δset1-infected mice that survived after 21 days appeared to be healthy, we examined the tissue condition of the infected but survived mice. No significant damage was observed in any C. albicans-infected tissues, including the kidney, spleen, and liver, compared to uninfected tissues. Moreover, C. albicans was not detected in any mouse tissues by staining with Periodic acid-Schiff (PAS), indicating that C. albicans had been cleared (Fig. S1).

Figure 1. Set1 has little effect on overall gene expression although it is required for virulence of C. albicans. A, mouse survival analysis. the 6-week-old female BALB/c mice are injected intravenously with 0.1 ml of 107 ml−1 cells. five mice were used for each strain. The Δset1 strain has attenuated virulence in mice as reported by Dr Clancy group [Citation15]. black closed circles, WT (SC5314, clinical isolated); Gray closed squares, Δset1. B, heatmap of total mRNA expression in WT and Δset1. Raw RNA-seq data were normalized by library size using DESeq2 and transformed into a log10 scale. the values were clustered and visualized using Pheatmap R Package. C, Venn diagram of top 50% genes with the highest expression levels. Top 50% expressed genes (3,097 genes in WT) were classified into three groups. Only 30 genes downregulated significantly in Δset1 (adjusted p value < 0.05)

Figure 1. Set1 has little effect on overall gene expression although it is required for virulence of C. albicans. A, mouse survival analysis. the 6-week-old female BALB/c mice are injected intravenously with 0.1 ml of 107 ml−1 cells. five mice were used for each strain. The Δset1 strain has attenuated virulence in mice as reported by Dr Clancy group [Citation15]. black closed circles, WT (SC5314, clinical isolated); Gray closed squares, Δset1. B, heatmap of total mRNA expression in WT and Δset1. Raw RNA-seq data were normalized by library size using DESeq2 and transformed into a log10 scale. the values were clustered and visualized using Pheatmap R Package. C, Venn diagram of top 50% genes with the highest expression levels. Top 50% expressed genes (3,097 genes in WT) were classified into three groups. Only 30 genes downregulated significantly in Δset1 (adjusted p value < 0.05)

For nonpathogenic C. albicans to become pathogenic in the host, large-scale alterations in gene expression are required. Previous studies have reported a positive correlation between H3K4 methylation mediated by Set1 and transcription [Citation13,Citation14,Citation23]; therefore, we performed RNA-seq to determine which genes were differentially transcribed in Δset1. Surprisingly, we found that most genes’ overall expression patterns were unaltered, even without H3K4 methylation (). Moreover, several pathogenesis-related genes exhibited an increased expression in Δset1 mutant (adjusted p-value <0.05) (Fig. S2A), contrary to our expectations from the defective virulence phenotype of Δset1 mutant. Most of these genes with the increased expression in Δset1 mutant were classified as cell wall proteins (Table S2).

We observed that only 3% of all annotated genes were differentially expressed in the absence of Set1 (adjusted p-value <0.05). In other words, the expression of 97% of genes was unchanged despite the lack of Set1-mediated H3K4 methylation. Remarkably, only 41 genes were downregulated in Δset1 (Table S3 and ), adjusted p-value <0.05). To allow the meaningful statistical analysis of expression levels, we selected the top 50% more highly expressed genes for further analysis (3,097 of 6,194 genes). Of these genes, only 30 genes showed a decrease in expression in the Δset1 mutant (, adjusted p-value <0.05). Because a previous study reported that Set1 can regulate pathogenicity by regulating adherence to host epithelial cells [Citation15], we selected 100 genes with GO slim for adhesion among all C. albicans genes using the Candida Genome Database (http://www.candidagenome.org/) and checked the expression levels of selected 100 genes. As a result, we confirmed that the expression of most genes (93 of 100 genes) was not significantly different in Δset1 mutant (p-value ≥0.05) (Table S4). In addition, adhesion of C. albicans is the first step in biofilm formation. When we checked the amount of biofilm formation through XTT assay, we did not observe significant differences in the amount of biofilm between WT and Δset1 mutant (Fig. S2B).

Figure 2. Set1 is required for proper expression of mitochondria-related genes and oxidative stress response genes. A, Heatmap of 41 Set1-regulated genes. The 41 genes are significantly downregulated in Δset1 versus WT (adjusted p value < 0.05). Raw RNA-seq data were normalized by library size using DESeq2 and transformed into a log10 scale. The values were clustered and visualized as a heatmap using Pheatmap R Package. The mitochondrial genes are indicated in bold. B, RNA-seq data was visualized in Integrative Genomics Viewer (IGV) for 10 Set1-regulated genes related to mitochondria. ADH1 gene is used as control. The y-axis range varies depending on the gene expression levels for each gene

Figure 2. Set1 is required for proper expression of mitochondria-related genes and oxidative stress response genes. A, Heatmap of 41 Set1-regulated genes. The 41 genes are significantly downregulated in Δset1 versus WT (adjusted p value < 0.05). Raw RNA-seq data were normalized by library size using DESeq2 and transformed into a log10 scale. The values were clustered and visualized as a heatmap using Pheatmap R Package. The mitochondrial genes are indicated in bold. B, RNA-seq data was visualized in Integrative Genomics Viewer (IGV) for 10 Set1-regulated genes related to mitochondria. ADH1 gene is used as control. The y-axis range varies depending on the gene expression levels for each gene

It was previously reported that the expression of most genes does not change in the Saccharomyces cerevisiae strain lacking Set1 [Citation24,Citation25]; therefore, the authors of this genome-wide study suggested that Set1-mediated H3K4 methylation does not regulate gene expression [Citation24]. Because the virulence phenotype was clearly reduced in C. albicans Δset1, we focused on the 41 genes whose expression decreased depending on H3K4 methylation. To examine the functional role of these downregulated genes in Δset1, we analyzed the functional annotation of each orthologous gene in S. cerevisiae, which is well studied in terms of functional genomics (Table S3). Although no gene was found to have virulence-related functions, functional annotation clustering (https://david.ncifcrf.gov/tools.jsp) analysis revealed that approximately 25% of the 41 Δset1-downregulated gene-encoded proteins were localized in the mitochondria, while most of the other gene products were located in the cytoplasm ( and Table S3). Approximately 13% of the 6,100 total proteins in yeast cells are generally localized in mitochondria [Citation26]; therefore, there were considerably more mitochondrial proteins among the Set1-dependent genes than average.

Set1 regulates the expression of mitochondrial genes involved in biogenesis and protection against oxidative damage

The mitochondrion is the primary site for cellular ROS production since its inner membrane houses the electron transfer system (ETS), a series of four large complexes (I to IV) that transfer electrons from donors to acceptors via redox reactions and pump protons out of the mitochondrial matrix, forming a proton gradient [Citation27–29]. Because this proton gradient drives proton-motive force, accumulated protons enter the matrix via ATP synthase to lower the proton gradient, leading to ATP synthesis. Electrons must pass through complexes to finally be absorbed as H2O, yet when electrons leak out, they react with O2 and become unstable superoxide anions (O2•). Generally, O2• is converted to hydrogen peroxide (H2O2) by superoxide dismutase (SOD), thereby reducing its reactivity; however, hydroxyl radicals (•OH) can be generated from H2O2 via the Fenton reaction and cause cellular damage through lipid oxidation, protein denaturation, and DNA mutations [Citation27–29].

We hypothesized that the virulence attenuation observed in Δset1 may be due to mitochondrial dysregulation. Thus, we focused on the 10 mitochondrial protein-coding genes, which were downregulated in the absence of Set1 in C. albicans (). The S. cerevisiae orthologs of the Set1-regulated genes, SOM1 (orf19.6359) and TOM5 (orf19.6247.1), play roles in the assembly and translocation of mitochondrial proteins, including the respiratory chain complex ( and Table S3) [Citation30–33]. A reduction in the expression of these genes could therefore reduce overall mitochondrial membrane integrity. Subsequently, the localization defect of the respiratory chains involved in oxidative phosphorylation could generate more mitochondrial ROS and eventually cause cellular damage. We also detected the downregulation of genes known to protect against oxidative stress damage in mitochondria, including NFU1 (orf19.2067) and AIM25 (orf19.3929; and Table S3) [Citation34–36]. We found that the expression of mitochondrial protein genes was downregulated in the absence of SET1, which may result in the production of many mitochondria-driven ROS. In addition, we found that the expression of some oxidative stress-responsive genes was downregulated (Table S3), suggesting that the SET1-deleted strain produces too much cellular ROS to remove, making Candida cells more sensitive to external attacks.

H3K4 methylation is highly enriched in set1-regulated genes

Because H3K4me3 is an important histone modification for transcriptional initiation, we confirmed that Set1 mediated all three H3K4 methylation states using antibodies that recognize H3K4me1, -me2, and -me3, respectively (). We already selected 10 mitochondrial protein-coding genes whose expression decreased in Δset1 strain (). We then carried out chromatin immunoprecipitation (ChIP) followed by qPCR (ChIP-qPCR) to determine whether these 10 mitochondrial protein-coding genes were regulated by H3K4 methylation. ChIP-qPCR with the H3K4me3 antibody revealed that H3K4me3 levels were high in the 5′ ORF of the mitochondrial protein-coding genes among the Set1-regulated genes in the WT (). Together, these results indicate that Set1 positively regulates target gene expression via H3K4me3, and that SET1 deletion reduces target gene expression.

Figure 3. H3K4me3 is enriched in 5ʹ ORF of Set1-regulated genes. A, Western blot analysis of H3K4 methylation in S. cerevisiae (Sc) and C. albicans (Ca). FM391 is used as a WT control for S. cerevisiae Δset1. Set1 is the sole methyltransferase in C. albicans. Histone H3 is used as a loading control. B, H3K4me3 ChIP followed by qPCR was performed in WT and Δset1. The 5ʹ end sequences of Set1-regulated genes related to mitochondria were used as amplicon. Intergenic region (IGR) is used as a negative control. All ChIP analyses were performed in two independent biological replicates and qPCR was performed in triplicated. *p < 0.05 and **P < 0.01

Figure 3. H3K4me3 is enriched in 5ʹ ORF of Set1-regulated genes. A, Western blot analysis of H3K4 methylation in S. cerevisiae (Sc) and C. albicans (Ca). FM391 is used as a WT control for S. cerevisiae Δset1. Set1 is the sole methyltransferase in C. albicans. Histone H3 is used as a loading control. B, H3K4me3 ChIP followed by qPCR was performed in WT and Δset1. The 5ʹ end sequences of Set1-regulated genes related to mitochondria were used as amplicon. Intergenic region (IGR) is used as a negative control. All ChIP analyses were performed in two independent biological replicates and qPCR was performed in triplicated. *p < 0.05 and **P < 0.01

Set1 deletion causes cellular ROS accumulation and renders Δset1 cells hypersensitive to external oxidative stress

We found that the expression of mitochondrial protein genes and oxidative stress response-related genes were downregulated in Δset1 cells ( and Table S4); therefore, we observed the cells under a fluorescence microscope to measure the levels of ROS generated using CellROX Green staining, which detects oxidative stress by binding to DNA when oxidized. Briefly, cells were treated mildly with zymolyase to produce spheroplasts for increased permeability and then stained with CellROX Green and DAPI. The CellROX Green signal was more potent in Δset1 than in WT cells (), suggesting that the absence of Set1 increases cellular ROS levels.

Figure 4. Cellular ROS is more generated in Δset1. A, Cellular oxidative status analysis in WT and Δset1. Cells were harvested in the early exponential phase and permeabilized with zymolyase followed by staining with 5 mM CellROX® Green. Cellular ROS was detected by fluorescence microscopy. C, Comparative growth of WT and Δset1 on H2O2-containing media. Overnight cultured cells were diluted into 107 CFU ml−1. 3 μl of 5-fold serial diluted cells were spotting on YPD or YPD added 5 mM H2O2. Each experiment was repeated at the same condition

Figure 4. Cellular ROS is more generated in Δset1. A, Cellular oxidative status analysis in WT and Δset1. Cells were harvested in the early exponential phase and permeabilized with zymolyase followed by staining with 5 mM CellROX® Green. Cellular ROS was detected by fluorescence microscopy. C, Comparative growth of WT and Δset1 on H2O2-containing media. Overnight cultured cells were diluted into 107 CFU ml−1. 3 μl of 5-fold serial diluted cells were spotting on YPD or YPD added 5 mM H2O2. Each experiment was repeated at the same condition

We hypothesized that SET1 deletion causes the mitochondrial membrane leakage, thereby accumulating cellular ROS and making cells susceptible to external oxidative stress. To test this hypothesis, we determined the sensitivity of the Δset1 strain to oxidative stress by spotting them on H2O2-containing media. We observed that the Δset1 strain was more sensitive to H2O2 than the WT strain (), consistent with our hypothesis.

Set1-dependent gene expression is required for C. albicans survival in the host

Macrophages play a vital role in the innate immune system by degrading pathogens and sterilizing tissues [Citation37]. When phagocytosis is induced by the host recognizing C. albicans, macrophages form inflammasomes and secrete ROS to attack the pathogen [Citation38]. In response to the host, C. albicans expresses hypha-specific genes, resulting in morphogenesis and establishing a defense against oxidative stress [Citation39]. Depending on the outcome of these processes, C. albicans is either cleared by macrophages or escapes from the phagosome and kills the macrophages [Citation38]. To determine whether Δset1 C. albicans was more vulnerable to attack by ROS-releasing macrophages in the host, we performed an in vitro interaction analysis between mouse macrophages and C. albicans.

To determine colony formation of C. albicans surviving macrophage attack, properly diluted C. albicans strains were cultured with the mouse macrophage cell line Raw264.7 (). We observed that the colonies formed by the Δset1 strain were too small to be characterized as colonies, unlike those formed by the WT strain (). Moreover, the Δset1 cells formed 50% fewer colonies than the WT cells (). These data suggest that Set1 contributes to C. albicans virulence by affecting resistance to and survival from macrophage attack. Because the virulence of Δset1 is undoubtedly attenuated both in vivo and in vitro, we concluded that Set1 is required for the full virulence phenotype of C. albicans when infecting host tissues.

Figure 5. Δset1 has attenuated virulence because of mitochondrial ROS leakage. A, The colony formation of C. albicans in the presence of RAW264.7 murine macrophage cells after 24 h incubation. The used amount of C. albicans and macrophage cells was represented in box. The representative microscopic images are shown. CA, C. albicans; MP, macrophage. B, The number of colonies was counted in the indicated ratio between C. albicans and macrophages and represented the average form. C, Proposed model. In WT, naturally occurring cellular ROS was cleared by ROS scavengers, represented as packman. On the other hand, because the mitochondrial protein genes are not properly expressed in Δset1, the mitochondrial membrane becomes loose and generates more cellular mitochondrial ROS than WT. Since the amount of ROS scavenger is the same as WT, the ROS accumulated in the cell is continuously neutralized by the antioxidant enzymes. However, if they are attacked by macrophage or treated with external ROS, Δset1 is more susceptible to oxidative stress because of the already generated cellular ROS

Figure 5. Δset1 has attenuated virulence because of mitochondrial ROS leakage. A, The colony formation of C. albicans in the presence of RAW264.7 murine macrophage cells after 24 h incubation. The used amount of C. albicans and macrophage cells was represented in box. The representative microscopic images are shown. CA, C. albicans; MP, macrophage. B, The number of colonies was counted in the indicated ratio between C. albicans and macrophages and represented the average form. C, Proposed model. In WT, naturally occurring cellular ROS was cleared by ROS scavengers, represented as packman. On the other hand, because the mitochondrial protein genes are not properly expressed in Δset1, the mitochondrial membrane becomes loose and generates more cellular mitochondrial ROS than WT. Since the amount of ROS scavenger is the same as WT, the ROS accumulated in the cell is continuously neutralized by the antioxidant enzymes. However, if they are attacked by macrophage or treated with external ROS, Δset1 is more susceptible to oxidative stress because of the already generated cellular ROS

The antioxidant enzymes that protect cells from ROS are conserved in C. albicans including SOD, glutathione peroxidase (GPX), and catalase (CAT) [Citation40–43] (Candida Genome Database; http://candidagenome.org/). The expression of antioxidant-encoded genes generally increases under oxidative stress conditions to neutralize ROS and prevent cell damage [Citation44]; however, we found that the expression of these genes did not change in the Δset1 strain, which contains high cellular ROS levels (Table S5), suggesting that Set1 deletion does not affect the ability to neutralize ROS. Because the Δset1 strain contains high cellular ROS levels, its antioxidants are saturated faster than in the WT, making Δset1 more vulnerable to ROS attack.

Overall, we propose that Set1 contributes to the virulence of C. albicans as shown in (). In WT cells, mitochondrial protein genes and oxidative stress response-related genes are generally expressed and localized correctly; therefore, mitochondrial ROS levels are low because WT C. albicans can immediately detoxify ROS as they are produced. Conversely, the mitochondrial proteins in Δset1 are not correctly localized, causing electrons to escape from the respiratory chain and increasing mitochondrial ROS levels. Although the ROS scavenging system in Δset1 eliminates ROS generated from the defective mitochondrial membrane, Δset1 results in high ROS levels even in the absence of external oxidative stress. When C. albicans is attacked by an external ROS derived from host macrophages or is grown in a medium supplemented with H2O2, the WT strain can remove ROS under the oxidative stress, but Δset1 cells cannot due to their saturated ROS scavenging system. Consequently, ROS quickly accumulates to a level that cannot be neutralized and decreases the survival rate of C. albicans Δset1 as the cells are unable to resist oxidative stress (). Therefore, we propose that Set1-mediated gene expression is required for the survival of C. albicans against oxidative stress and its pathogenicity.

Discussion

Opportunistic pathogens can become pathogenic via direct or indirect mechanisms: direct mechanisms involve a change in the ability of pathogens to attack hosts by altering their virulence factor expression, whereas indirect mechanisms involve pathogens passively enduring and surviving against attacks by the host immune system [Citation45]. To understand why the virulence of C. albicans decreases in the absence of Set1, we identified genes regulated by Set1 using RNA-seq. The expression of genes directly involved in virulence did not change or was even upregulated in the Δset1strain; however, the expression of many genes required to resist oxidative stress, including mitochondrial protein genes, decreased significantly. The reduced expression of these genes increased the generation of mitochondrial ROS but did not reduce the growth rate compared to the WT under normal conditions (data not shown), since the generated cellular ROS could be neutralized continuously by antioxidant enzymes, such as SOD. However, the Δset1 strain was more sensitive than the WT strain when treated with H2O2 or co-incubated with macrophages. This hyper-sensitivity was likely due to the inability of the Δset1 strain to neutralize the excessive amounts of ROS generated during H2O2 treatment or co-incubation with macrophages. Indeed, the failure of the Δset1 strain to induce oxidative stress-related genes in response to excessive ROS levels led to C. albicans cell death. These findings suggest that Set1 did not directly affect virulence by promoting the expression of virulence-related genes but instead lowered the production of mitochondrial ROS by regulating the expression of mitochondrial protein genes. Therefore, antioxidant enzymes remained available for external ROS attack to reduce ROS levels to normal physiological conditions. Together, these findings indicate that Set1 allows C. albicans to survive longer in the host and eventually switch to its pathogenic form ().

Histone H3K4 methylation and the H3K4 methyltransferase, Set1, are well conserved in many eukaryotes, and several studies have reported that Set1 and Set1-mediated H3K4 methylation are important for the pathogenicity of some fungal pathogens, in which Set1 regulates virulence and the stress response by controlling the expression of specific genes. For instance, Set1 has been shown to be involved in the virulence of the plant pathogen Fusarium verticillioides by regulating the expression of fumonisin B1 toxin-encoded genes [Citation17], whereas Set1-mediated H3K4 methylation activates the expression of TR1, which encodes the toxin deoxynivalenol (DON), in the plant fungal pathogen Fusarium graminearum [Citation20]. In addition, Set1 is involved in Magnaporthe oryzae fungal virulence by activating virulence-related genes via H3K4 methylation [Citation21], while an H3K4 methyltransferase is required to induce genes in the entomopathogenic fungus Metarhizium robertsii under host conditions and induce virulence to mosquito infection [Citation19]. Importantly, C. albicans is the only example that the virulent effects of Set1 have been reported in a human pathogen [Citation15].

H3K4me3 is a marker of active transcription that is enriched in actively transcribed genes; however, we found that the absence of SET1 did not change 97% of the total gene expression (). Similarly, recent genome-wide studies revealed that the absence of H3K4 methyltransferase did not dramatically change the overall gene expression in other fungal organisms, such as S. cerevisiae [Citation46,Citation47], M. oryzae [Citation21], and Fusarium fujikuroi [Citation18]. Because we observed that the expression of some genes increased in the absence of H3K4 methyltransferase, it is difficult to assume that the role of H3K4 methylation simply correlates with active transcription. Although the effect of Set1 was not significant in terms of the genome-wide expression profile, H3K4 methyltransferase has been implicated in the pathogenicity phenotype of various pathogenic fungi, including C. albicans [Citation48]. In this study, we described the mechanism via which the absence of Set1 in C. albicans attenuates virulence, with a particular focus on the genes regulated by H3K4 methylation. Transcriptome analysis revealed that RNA levels did not change significantly in C. albicans in the presence or absence of Set1 under normal conditions; however, there may be a set of genes whose expression differs depending on the presence or absence of Set1 under oxidative stress conditions or during interactions with macrophage. Therefore, the role of H3K4 methylation in transcriptional regulation cannot simply be defined as active or negative regulation, but must be described as a function that regulates each gene differently in specific environments and should be investigated in future studies.

Author contributions

J.-S. L. designed the research; J.K. and S.P. performed the research; S.K revised the paper; J.K, J.-S.L, and E.-J.L. wrote the paper.

Supplemental material

Supplemental Material

Download (6 MB)

Acknowledgments

We thank Ali Shilatifard for providing the antibodies and Cornelius Clancy for kindly giving the Set1 deletion mutant. This work was supported by the National Research Foundation of Korea under Grant [No. NRF-2015R1D1A1A02061743; NRF-2018R1D1A1A02048280; and NRF-2020R1I1A3072234].

Disclosure statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Availability of data and materials

Sequencing data sets generated and/or analyzed during this study have been deposited to the Gene Expression Omnibus Database (GEO) under the accession number: GSE158481 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE158481). All data generated or analyzed during this study are included in this published article and its Additional files.

Supplementary material

Supplemental data for this article can be accessed here.

Additional information

Funding

This work was supported by the National Research Foundation of Korea under Grant [No. NRF-2015R1D1A1A02061743; NRF-2018R1D1A1A02048280; NRF-2020R1I1A3072234; NRF-2020M3A9H5104235; 2021R1I1A1A01057137].

References

  • Gow NAR, Van DV, Brown AJP, et al. Candida albicans morphogenesis and host defence: discriminating invasion from colonization. Nat Rev Microbiol. 2011;10(2):112.
  • Calderone R, Fonzi W. Virulence factors of Candida albicans. Trends Microbiol. 2001;9(7):327–335.
  • Kumamoto CA, Vinces MD. Contributions of hyphae and hypha-co-regulated genes to Candida albicans virulence. Cell Microbiol. 2005;7(11):1546–1554.
  • Soll DR. Why does Candida albicans switch? FEMS Yeast Res. 2009;9(7):973–989.
  • Tsui C, Kong EF, Jabra-Rizk MA. Pathogenesis of Candida albicans biofilm. Pathog Dis. 2016;74(4):ftw018.
  • Mayer FL, Wilson D, Hube B. Candida albicans pathogenicity mechanisms. Virulence. 2013;4(2):119–128.
  • Naglik JR, Challacombe SJ, Hube B. Candida albicans Secreted Aspartyl Proteinases in Virulence and Pathogenesis. Microbiol Mol Biol Rev. 2003;67(3):400–428.
  • Noble SM, Gianetti BA, Witchley JN. Candida albicans cell-type switching and functional plasticity in the mammalian host. Nat Rev Microbiol. 2017;15(2):96–108.
  • Hube B. From commensal to pathogen: stage- and tissue-specific gene expression of Candida albicans. Curr Opin Microbiol. 2004;7(4):336–341.
  • Biswas S, Van Dijck P, Datta A. Environmental Sensing and Signal Transduction Pathways Regulating Morphopathogenic Determinants of Candida albicans. Microbiol Mol Biol Rev. 2007;71:348.
  • Berger SL. An embarrassment of niches: the many covalent modifications of histones in transcriptional regulation. Oncogene. 2001;20(24):3007.
  • Li B, Carey M, Workman JL. The Role of Chromatin during Transcription. Cell. 2007;128(4):707–719.
  • Bernstein BE, Humphrey EL, Erlich RL, et al. Methylation of histone H3 Lys 4 in coding regions of active genes. Proc Natl Acad Sci U S A. 2002;99(13):8695.
  • Santos-Rosa H, Schneider R, Bannister AJ, et al. Active genes are tri-methylated at K4 of histone H3. Nature. 2002;419(6905):407–411.
  • Raman SB, Hong Nguyen M, Zhang Z, et al. Candida albicans SET1 encodes a histone 3 lysine 4 methyltransferase that contributes to the pathogenesis of invasive candidiasis. Mol Microbiol. 2006;60(3):697–709.
  • Dallery J, Adelin É, Le Goff G, et al. H3K4 trimethylation by CclA regulates pathogenicity and the production of three families of terpenoid secondary metabolites in Colletotrichum higginsianum. Mol Plant Pathol. 2019;20(6):831–842.
  • Gu Q, Tahir AH, Zhang H, et al. Involvement of FvSet1 in Fumonisin B1 Biosynthesis, Vegetative Growth, Fungal Virulence, and Environmental Stress Responses in Fusarium verticillioides. Toxins (Basel). 2017;9(2):43.
  • Janevska S, Güldener U, Sulyok M, et al. Set1 and Kdm5 are antagonists for H3K4 methylation and regulators of the major conidiation-specific transcription factor gene ABA1 in Fusarium fujikuroi. Environ Microbiol. 2018;20(9):3343–3362.
  • Lai Y, Cao X, Chen J, et al. Coordinated regulation of infection-related morphogenesis by the KMT2-Cre1-Hyd4 regulatory pathway to facilitate fungal infection. Sci Adv. 2020;6(13):eaaz1659–eaaz1659.
  • Liu Y, Zhang M, Xie R, et al. The Methyltransferase AflSet1 Is Involved in Fungal Morphogenesis, AFB1 Biosynthesis, and Virulence of Aspergillus flavus. Front Microbiol. 2020;11:234.
  • Liu Y, Liu N, Yin Y, et al. Histone H3K4 methylation regulates hyphal growth, secondary metabolism and multiple stress responses in Fusarium graminearum. Environ Microbiol. 2015;17(11):4615–4630.
  • Pham KTM, Inoue Y, Vu BV, et al. MoSET1 (Histone H3K4 Methyltransferase in Magnaporthe oryzae) Regulates Global Gene Expression during Infection-Related Morphogenesis. PLoS Genet. 2015;11(7):e1005385.
  • Kim J, Lee J. Rapid method for chromatin immunoprecipitation (ChIP) assay in a dimorphic fungus, Candida albicans. J Microbiol. 2020;58(1):11–16.
  • Ng HH, Robert F, Young RA, et al. Targeted Recruitment of Set1 Histone Methylase by Elongating Pol II Provides a Localized Mark and Memory of Recent Transcriptional Activity. Mol Cell. 2003;11(3):709–719.
  • Howe FS, Fischl H, Murray SC, et al. Is H3K4me3 instructive for transcription activation? Bioessays. 2017;39(1):e201600095.
  • Murray SC, Lorenz P, Howe FS, et al. H3K4me3 is neither instructive for, nor informed by, transcription. bioRxiv. 2019;709014. 10.1101/709014
  • Kumar A, Agarwal S, Heyman JA, et al. Subcellular localization of the yeast proteome. Genes Dev. 2002;16(6):707–719.
  • Cadenas E. Mitochondrial free radical production and cell signaling. Mol Aspects Med. 2004;25(1–2):17–26.
  • Calderone R, Li D, Traven A. System-level impact of mitochondria on fungal virulence: to metabolism and beyond. FEMS Yeast Res. 2015;15(4). 10.1093/femsyr/fov027
  • Turrens JF. Mitochondrial formation of reactive oxygen species. J Physiol. 2003;552(2):335–344.
  • Dietmeier K, Hönlinger A, Bömer U, et al. Tom5 functionally links mitochondrial preprotein receptors to the general import pore. Nature. 1997;388(6638):195–200.
  • Jan PS, Esser K, Pratje E, et al. Som1, a third component of the yeast mitochondrial inner membrane peptidase complex that contains Imp1 and Imp2. Mol. Gen. Genet. 2000;263(3):483–491
  • Liang H, Luo W, Green N, et al. Cargo Sequences Are Important for Som1p-dependent Signal Peptide Cleavage in Yeast Mitochondria. J Biol Chem. 2004;279(38):39396–39400.
  • Mokranjac D, Neupert W. Thirty years of protein translocation into mitochondria: unexpectedly complex and still puzzling. Biochim. Biophys. Acta Mol. Cell Res. 2009;1793(1):33–41
  • Aguilar-Lopez J, Laboy R, Jaimes-Miranda F, et al. Slm35 links mitochondrial stress response and longevity through TOR signaling pathway. Aging (Albany NY). 2016;8(12):3255–3271.
  • Blankenship JR, Fanning S, Hamaker JJ, et al. An Extensive Circuitry for Cell Wall Regulation in Candida albicans. PLOS Pathog. 2010;6(2):e1000752.
  • Melber A, Na U, Vashisht A, et al. Role of Nfu1 and Bol3 in iron-sulfur cluster transfer to mitochondrial clients. eLife. 2016;5:e15991.
  • Martin CJ, Peters KN, Behar SM. Macrophages clean up: efferocytosis and microbial control. Curr Opin Microbiol. 2014;17:17–23.
  • Dantas AD, Lee KK, Raziunaite I, et al. Cell biology of Candida albicans–host interactions. Curr Opin Microbiol. 2016;34:111–118.
  • Dantas AD, Day A, Ikeh M, et al. Oxidative Stress Responses in the Human Fungal Pathogen, Candida albicans. Biomolecules. 2015;5:
  • Hwang C, Rhie G, Oh J, et al. Copper- and zinc-containing superoxide dismutase (Cu/ZnSOD) is required for the protection of Candida albicans against oxidative stresses and the expression of its full virulence. Microbiology. 2002;148(11):3705–3713.
  • Miramón P, Dunker C, Kasper L, et al. A family of glutathione peroxidases contributes to oxidative stress resistance in Candida albicans. Med Mycol. 2014;52(3):223–239.
  • Nakagawa Y, Watanabe KK, Kiyomi, Morita Y, et al. The Expression of the Pathogenic Yeast Candida albicans Catalase Gene in Response to Hydrogen Peroxide. Microbiol Immunol. 1999;43(7):645–651.
  • Rhie G, Hwang C, Brady M, et al. Manganese-containing superoxide dismutase and its gene from Candida albicans. Biochim. Biophys. Acta, Gen. Subj. 1999;1426(3):409–419
  • Scandalios JG. Oxidative stress: molecular perception and transduction of signals triggering antioxidant gene defenses. Braz J Med Biol Res. 2005;38(7):995–1014.
  • Casadevall A, Pirofski L. Accidental Virulence, Cryptic Pathogenesis, Martians, Lost Hosts, and the Pathogenicity of Environmental Microbes. Eukaryot Cell. 2007;6(12):2169.
  • Martín GM, King DA, Green EM, et al. Set5 and Set1 cooperate to repress gene expression at telomeres and retrotransposons. Epigenetics. 2014;9(4):513–522.
  • Ramakrishnan S, Pokhrel S, Palani S, et al. Counteracting H3K4 methylation modulators Set1 and Jhd2 co-regulate chromatin dynamics and gene transcription. Nat Commun. 2016;7(1):11949.