Previous Article | Next Article ![]()
Eukaryotic Cell, August 2009, p. 1174-1183, Vol. 8, No. 8
1535-9778/09/$08.00+0 doi:10.1128/EC.00074-09
Copyright © 2009, American Society for Microbiology. All Rights Reserved.
,
Biotechnology Research Institute, National Research Council of Canada, Montréal, Québec H4P 2R2, Canada,1 Department of Anatomy and Cell Biology, McGill University, Montréal, Québec H3A 1B1, Canada,2 Institut de Recherche en Biologie Végétale, Université de Montréal, Montréal, Québec H1X 2B2, Canada3
Received 4 March 2009/ Accepted 12 June 2009
|
|
|---|
|
|
|---|
The arsenal of clinically active antifungal compounds is limited. The eukaryotic nature of C. albicans makes it similar to its human host, thus reducing the number of potential drug targets (11). Most Candida infections can be treated with inhibitors that target either the biosynthesis of ergosterol, the main sterol of fungal membranes, or the biosynthesis of the key component of the fungal cell wall, (1,3)-β-D-glucan. However, the emergence of drug resistance in this pathogen sometimes reduces the effectiveness of antifungal drugs. C. albicans has evolved different resistance mechanisms to bypass the inhibitory effects of azole drugs (11), mainly through the overexpression of the ATP-binding cassette (ABC) transporter Cdr1p or Cdr2p or the major facilitator superfamily (MFS) transporter Mdr1p (36). Recently, Morschhauser et al. (32) have elucidated a mechanism that results in the overexpression of MDR1 in clinical isolates. Drug resistance was caused by a gain-of-function mutation in a zinc cluster transcription factor (TF) called Mrr1p. Inactivation of MRR1 abolishes the resistance of Mdr1p-overexpressing strains. Gain-of-function mutations in the zinc cluster TF Tac1p have also been shown to result in Cdr1p or Cdr2p overexpression (9, 10). Genome-wide location analysis using chromatin immunoprecipitation coupled to microarrays (ChIP-chip) showed that, in addition to Cdr1p and Cdr2p, Tac1p also targets the promoters of other resistance genes, such as the integral membrane flippase Rta3p and stress-related genes (29). Furthermore, Dunkel et al. (14) showed that azole resistance is linked to a gain-of-function mutation in the sterol metabolism regulator Upc2p, leading to constitutive expression of ergosterol biosynthesis genes. Genome-wide occupancy of Upc2p confirmed the direct binding of this TF to ergosterol genes, including the azole target Erg11p (45).
The TF Ndt80p has been identified as a positive regulator of CDR1 by screening a C. albicans expression library in a Saccharomyces cerevisiae strain carrying a CDR1-lacZ reporter construct (4). Furthermore, it has been shown that NDT80 is crucial for azole drug tolerance and activation of CDR1 expression in response to miconazole (4). The homologue of Ndt80p in S. cerevisiae regulates gene expression of middle sporulation genes and is required for exit from pachytene and for full meiotic recombination (23, 35, 37, 44). Thus far, the accurate role of this TF in regulating drug resistance in C. albicans remains to be assessed. To gain a more complete understanding of the role of Ndt80p in mediating azole resistance in C. albicans, we set out to investigate its genomic occupancy using genome-wide location analysis (ChIP-chip). By combining genome-wide location and expression analyses, we have characterized the Ndt80p fluconazole-dependent regulon for the first time and elucidated the key role of this TF as an activator of ergosterol metabolism genes.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Candida albicans strains used in this study
|
Gene disruption was accomplished as follows. Plasmid pFA-HIS1 was used as the template to prepare the HIS1-PCR cassette, and plasmid pFA-URA3 was used to prepare the URA3 PCR cassette, as previously described by Gola et al. (21). Ndt80 was deleted in two steps. In the first step, one allele was replaced by homologous recombination with a PCR cassette containing the HIS1 gene. The BWP17 strain was transformed with the PCR fragments, and the cells were plated onto YPD plates minus histidine. Genomic DNA of the positive colonies was analyzed by PCR for the proper integration site of the cassette. The strain produced by deleting one allele of Ndt80 was named AS30. The second allele disruption was accomplished as done for the first allele using the URA3 PCR cassette, and the cells were plated onto YPD plates minus uridine. The colonies were again analyzed by PCR to identify the proper integration site. The complete deletion of candidate genes was confirmed by PCR, using primers internal to the recombination sites amplifying open reading frame (ORF) regions. The strain produced by deleting both alleles from Ndt80 (named AS31) was treated with 5-fluoroorotic acid to recycle the URA3 marker, and the resulting uridine-negative strain was named AS32.
The absence of aneuploidy was confirmed for mutant and revertant strains using comparative genomic hybridization. C. albicans genomic DNA was isolated from a saturated overnight culture with the Qiagen genomic DNA extraction kit and labeled with either Cy3 or Cy5 dye with the BioPrime CGH labeling kit (Invitrogen). Unincorporated nucleotides were removed with Qiagen PCR columns, and the labeled probes were then hybridized to DNA microarrays, as described by Nantel et al. (34).
For reintegration experiments, the NDT80 gene was reintegrated into the null mutant ndt80 strain. The wild-type NDT80 gene was amplified from genomic DNA using oligonucleotides REVF1 and REVR1 (Table 2) and Expand high-fidelity polymerase (Roche). The PCR fragment was digested with restriction enzymes ClaI and XhoI and cloned in the same sites of the CIp10 vector (33). The plasmid was sequenced to confirm the integrity of the NDT80 gene. Plasmid CIp10-NDT80 was digested with the StuI restriction enzyme and used to transform the ndt80 mutant strain (AS32). The uridine-positive colonies were analyzed by PCR, and the obtained wild-type fragment confirmed the reintegration of the NDT80 gene.
|
View this table: [in a new window] |
TABLE 2. Primers used in this study
|
Whole-genome location profiling by ChIP-chip and data analysis. Cells were grown to an optical density at 600 nm (OD600) of 2 in 40 ml of YPD. The subsequent steps of DNA cross-linking, DNA shearing, ChIP, and DNA labeling with Cy dyes were conducted exactly as described by Lavoie et al. (26). Tiling arrays were cohybridized with tagged immunoprecipitated (Cy5-labeled) and mock immunoprecipitated (untagged BWP17 strain; Cy3-labeled) DNA samples. Microarray hybridizing, washing, and scanning were performed, as described by Nantel et al. (34). Prehybridization and hybridization solutions consisted of DIG Easy Hyb solution (Roche Diagnostics, Mannheim, Germany), with 0.45% salmon sperm DNA and 0.45% yeast tRNA. The hybridization was carried out at 42°C for 20 h in a SlideBooster Hyb chamber SB800 (Advalytix, Brunnthal, Germany), with regular microagitation of the sample. Slides were washed once in 1.0% SSC (1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate), once in 0.2% sodium dodecyl sulfate at 42°C for 5 min, twice in 0.1% SSC, twice in 0.2% sodium dodecyl sulfate at 42°C for 5 min, and once in 0.1% SSC at 24°C for 5 min, followed by four rinses in 0.1% SSC. Chips were air dried before being scanned using a ScanArray Lite microarray scanner (PerkinElmer). Fluorescence intensities were quantified using ImaGene software (BioDiscovery, Inc.), background corrected, and normalized for signal intensity (using Lowess normalization). The significance cutoff was determined using the distribution of log ratios for each factor. It was set at two standard deviations from the mean of log-transformed enrichments. Values shown are an average of two biological replicates derived from independently isolated transformants of tagged and mock constructs. Peak detection was performed using Gaussian edge detection applied to the smoothed probe signal curve, as described by Tuch et al. (39).
Motif discovery. For de novo identification of the consensus-binding site, the sequences were analyzed with the Multiple EM for Motif Elucidation (MEME) program (http://meme.sdsc.edu/meme/website/meme.html). For sequences of 150 bp upstream and downstream, the detected signal peak was used. The motifs were allowed to have any sequence and length (between 6 and 50 bp) and could be present anywhere in the sequence.
Gene expression profiling. Overnight cultures of designated C. albicans strains were diluted to a starting OD600 of 0.1 in 100 ml fresh YPD-uridine media. The cultures were incubated with shaking at 30°C to an OD600 of 0.8 and split into two 50-ml cultures. Fluconazole was added to the experimental culture to a final concentration of 100 µg per ml, while an equal volume of dimethyl sulfoxide was added to the control culture. The cultures were incubated for 60 min. Cells were harvested by centrifugation and stored at –80°C.
To extract RNA from cells, samples stored at –80°C were placed on ice, and RNeasy buffer RLT was added to pellets at a ratio of 10:1 (vol/vol) buffer/pellet. The pellet was allowed to thaw in the buffer while being vortexed briefly at high speed. The resuspended pellet was placed back on ice and divided into 1-ml aliquots in 2-ml screw-cap microcentrifuge tubes containing 0.6 ml of 3-mm-diameter acid-washed glass beads. Samples were homogenized five times, for 1 min each, at 4,200 rpm using BeadBeater. Samples were placed on ice for 1 min after each homogenization step. After homogenization, the Qiagen RNeasy protocol was followed, as recommended. Total RNA samples were eluted in RNase-free H2O. RNA quality and integrity were assessed using an Agilent 2100 bioanalyzer.
cDNA labeling and microarray production were performed, as previously described (34). Briefly, 20 µg of total RNA was reverse transcribed using oligo(dT)21 in the presence of Cy3 or Cy5-dCTP (Invitrogen) and SuperScript III reverse transcriptase (Invitrogen). Thereafter, template RNA was degraded by adding 2.5 units RNase H (USB) and 1 µg RNase A (Pharmacia), followed by incubation for 15 min at 37°C. The labeled cDNAs were purified with the QIAquick PCR purification kit (Qiagen). Prior to hybridization, Cy3/Cy5-labeled cDNA was quantified using a NanoDrop ND-1000 UV-Vis spectrophotometer (NanoDrop) to confirm dye incorporation. Microarray hybridization, washing, and scanning were performed as described for ChIP-chip experiments. Genes with statistically significant changes in transcript abundance in each experiment were identified using a 1.5 cutoff point and Welch's t test with a false discovery rate of less than 5%.
|
|
|---|
The vast majority of Ndt80p-bound sites were found at gene promoter regions. Intragenic occupancy was exclusively detected in ORFs that are considered spurious or dubious (2). Among the 90 ORFs bound by Ndt80p, 25 are annotated as spurious, 49 as dubious, and 16 as uncharacterized ORFs not experimentally validated. Dubious and spurious ORFs have no orthologs in other eukaryotes and had expression profiles that were not significantly correlated with those of other genes in the genome. Consequently, those regions are more likely to be intergenic rather than intragenic sequences. We thus conclude that Ndt80p binding occurs specifically in promoter regions.
Ndt80p targets are enriched in genes related to metabolism, stress, development, and drug resistance. To reveal general functional features of the biological processes modulated by Ndt80p, we conducted a gene ontology (GO) analysis by analyzing the 1,446 genes whose promoters are associated with Ndt80p. For these analyses, all GO categories of genes having an enrichment P value of <0.0001 were selected. As shown in Fig. 1A, we found that Ndt80p binds to a wide range of genes with diverse biological functions. This analysis revealed a significant enrichment in genes related to the cell wall (P = 8.31e–14), carbohydrate metabolism (P = 6.84e–09), response to heat (P = 4.43e–06) and osmotic stress (P = 5.56e–04), alcohol metabolism (P = 5.45e–06), cell adhesion (P = 2e–05), glycolysis (P = 8.39e–05), the G1/S cell cycle (P = 7.76e–04), hyphal growth (P = 5.1e–04), and amine transmembrane transport (P = 2.36e–04). Examination of Ndt80p target genes revealed that they were remarkably enriched in promoter genes coding for transcription factors and transcriptional regulators (P = 2.65e–04) associated with each of these different functions. Transcriptional regulators that were found among Ndt80p targets include regulators of hyphal growth (Efg1p, Nrg1p, Ume6p, Tec1p, Cph2p, Flo8p, Czf1p, Ssn6p, Rfg1p), carbohydrate metabolism (Rgt1p, Tye7p, Gal4p, Mig1p), the cell cycle (Swi4p, Ash1p), lipid metabolism (Ino2p, Opi1p, Ctf1p), translation and amino acid metabolism (Cbf1p, Gln3p, Gcn4p), stress (Cat8p, Hac1p, Cas5p), and general transcriptional regulators (Sua71p, Tbp1p, Stp1p, Stp2p, Stp3p, Stp4p).
![]() View larger version (22K): [in a new window] |
FIG. 1. (A) GO biological process annotation of Ndt80p-bound promoters. The P value was calculated using hypergeometric distribution, as described on the GO Term Finder website (www.candidagenome.org/cgi-bin/GO/goTermFinder). (B) CDR1 promoter occupancy by Ndt80p.
|
De novo motif analysis of Ndt80p-bound promoters. Structural studies demonstrated that the Ndt80p homologue of S. cerevisiae binds to the middle sporulation element (MSE) (18, 25, 31, 37). Based on both statistical and mutational analyses, the MSE has been well defined and corresponds to the 9-bp consensus, as follows: 5'-gNCRCAAAY-3' (where lowercase letters indicate semiconserved residues, R indicates a purine, N indicates any nucleotide, and Y indicates either a thymine or a cytosine) (25, 37, 43). Using 100 experimentally determined Ndt80-bound loci ranking at the top of highly enriched peaks, we sought to characterize properties of Ndt80p-DNA interaction by assessing motif enrichment to detect de novo consensus of Ndt80p. As shown in Fig. 2, our result illustrated that the top-scoring minimotifs correspond to the consensus, as follows: 5'-NaCacAAAa-3' (P = 10e–27). This motif is highly similar to that of the S. cerevisiae MSE, suggesting a conservation of the Ndt80p binding site in C. albicans. To assess the occurrence of the core MSE consensus, 5'-ACACAAA-3', in all Ndt80p targets, we scanned the 1,446 promoters and found that 1,326 of these contain the Ndt80p binding consensus.
![]() View larger version (24K): [in a new window] |
FIG. 2. Logo of the top-scoring motif discovered by MEME using 100 experimentally determined Ndt80-bound loci, ranking at the top of the highly enriched peaks. Saccharomyces cerevisiae MSE consensus is shown on the bottom of the logo.
|
![]() View larger version (38K): [in a new window] |
FIG. 3. The absence of Ndt80p causes sensitivity to ergosterol biosynthesis inhibitors. Serial dilutions (10-fold) of wild-type (wt), revertant, and ndt80 mutant strains were grown on YPD supplemented with fluconazole at 10 µg/ml (A) or terbinafine at 2 µg/ml and lovastatin at 20 µg/ml (B) and grown at 30°C for 48 h.
|
|
View this table: [in a new window] |
TABLE 3. Effect of fluconazole on gene expression of the C. albicans wild-type strain and the ndt80 mutanta
|
![]() View larger version (17K): [in a new window] |
FIG. 4. Real-time quantitative PCR validation of microarray data. Expression levels of ERG3, ERG4, ERG6, ERG7, ERG11, ERG24, and ERG25 were determined in the wild type (wt), the ndt80 mutant, and the revertant strain, using ACT1 as a reference. The reported values are the means ± standard deviations from two independent experiments.
|
![]() View larger version (25K): [in a new window] |
FIG. 5. Effect of fluconazole on expression of C. albicans post-squalene ergosterol biosynthesis genes carried by the wild-type strain and ndt80 mutant. Ndt80p promoter binding of ergosterol genes is also shown. The positions of the MSE motifs are shown.
|
|
|
|---|
In this study, we used ChIP-chip technology to investigate the roles of Ndt80p in C. albicans. Several lines of evidence, summarized below, demonstrate the usefulness of this strategy. Indeed, out of 100 promoters that were randomly selected for the confirmation of the ChIP tiling array results, none were shown to be false positive. Furthermore, our results confirmed the observation that Ndt80p binds to the promoter region of CDR1, as suggested by Chen et al. (4), who demonstrated that this TF alleviates the LacZ activity of a CDR1-lacZ construct. Additionally, de novo motif analysis revealed a high enrichment of the MSE motif in Ndt80p promoter targets, thus demonstrating a high degree of conservation for this binding site in C. albicans. This was expected since the Ndt80p DNA-binding domain is conserved between C. albicans and S. cerevisiae. Thus, our results demonstrate that the ChIP-chip technique used here is a compelling tool for genome-wide occupancy analysis and for devising the function of transcription factors.
The finding that Ndt80p occupies the promoter regions of 23% of C. albicans genes suggests that this TF plays a global role in transcriptional regulation. This hypothesis is buttressed by the fact that the list of Ndt80p targets shows remarkable enrichment for genes whose products are transcriptional regulators acting in different biological process. This suggests that Ndt80p is a high-level hierarchical regulator acting in the transcriptional cascade of diverse biological process. Intriguingly, Ndt80p was found to bind promoters of hyphal growth regulators, such as the TF Cph2p, Efg1p, Ume6p and Flo8p, which activate filamentation-specific genes when C. albicans undergoes morphological switching. Based on this observation, Ndt80p might play a critical role in regulating morphogenesis in this fungal pathogen. Indeed, the ndt80 mutant is not able to form hyphae under different filamentation-promoting conditions (A. Sellam and A. Nantel, unpublished data), and the precise mechanisms involving Ndt80p in hyphal growth are under investigation.
To help with the identification of genes that require Ndt80p for fluconazole-dependent expression, we first produced our own list of fluconazole-responsive genes carried by the wild-type strain. Only 10 of our fluconazole-responsive genes were shown to be activated during the exposure to the same drug by Lepak et al. (27), with 8 of these common genes being involved in ergosterol biosynthesis (ERG1, ERG3, ERG4, ERG5, ERG6, ERG11, ERG24, ERG26). Copping et al. (8) have also investigated the transcriptional response of C. albicans to fluconazole. Likewise, 10 upregulated genes were common with those of our study, including 8 ERG genes (ERG1, ERG3, ERG4, ERG5, ERG6, ERG7, ERG11, ERG251) and also the fatty acyl coenzyme A synthase FAS2. These findings suggests that, even when using different strains and growth and treatment conditions on different microarray platforms, C. albicans responds to fluconazole by activating mainly the ergosterol pathway, as was indeed demonstrated by several fungi following the exposure to other azole drugs (1, 12, 7, 27, 28). In the wild-type strain, upregulation of ERG genes is considered to be a compensatory response used to overcome the inhibition of the lanosterol demethylase Erg11p by fluconazole. Based on our transcriptional profiling data, ndt80 was not able to activate as many ERG genes as in the wild-type strain, suggesting that these compensatory mechanisms are altered in this mutant. Thus, our study suggests that Ndt80p affects drug resistance by controlling the expression of genes involved in the biosynthesis of ergosterol. Direct analysis using ChIP-chip revealed that many promoters of ERG genes, namely, ERG3, ERG4, ERG6, ERG7, ERG11, ERG13, ERG24, ERG25, and ERG251, were bound by Ndt80p in vivo, thus suggesting that this TF directly regulates their expression.
Ndt80p is a transcriptional activator that plays a key role in the progression of the meiotic divisions in the yeast S. cerevisiae (23). In C. albicans, Ndt80p has been identified as a positive regulator of CDR1p (4). C. albicans Ndt80p appears to have functionally diverged from the role of its S. cerevisiae homologue, since ScNdt80p fails to activate the expression of the CDR1p-lacZ heterologous reporter gene (42). In our study, we demonstrate that Ndt80p binds the promoter of ERG biosynthesis genes and is required for their transcriptional activation in response to fluconazole. Such function has not been demonstrated in S. cerevisiae. Indeed, genome-wide mapping of Ndt80p in S. cerevisiae was performed; however, no significant binding was observed in the promoters of ergosterol biosynthesis genes (22). Furthermore, even with the absence of a complete sexual cycle in C. albicans, meiosis genes are quite conserved (40). However, no significant enrichment for genes involved in meiosis was obtained in this study. Taken together, these findings are a further confirmation that the Ndt80p regulon has been rewired in C. albicans.
This work was supported by grant MOP-42516 from the Canadian Institutes of Health Research (CIHR) to A.N.
This is NRC manuscript 50651.
Published ahead of print on 19 June 2009. ![]()
Supplemental material for this article may be found at http://ec.asm.org/. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2010 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»