RNA sequencing revealed novel actors of the acquisition of drug resistance in Candida albicans
- Sanjiveeni Dhamgaye1,
- Maria Bernard3, 4, 5,
- Gaelle Lelandais6,
- Odile Sismeiro7,
- Sophie Lemoine3, 4, 5,
- Jean-Yves Coppée7,
- Stéphane Le Crom3, 4, 5,
- Rajendra Prasad†1Email author and
- Frédéric Devaux†2, 8
© Dhamgaye et al.; licensee BioMed Central Ltd. 2012
Received: 25 May 2012
Accepted: 3 August 2012
Published: 16 August 2012
Drug susceptible clinical isolates of Candida albicans frequently become highly tolerant to drugs during chemotherapy, with dreadful consequences to patient health. We used RNA sequencing (RNA-seq) to analyze the transcriptomes of a CDR (Candida Drug Resistance) strain and its isogenic drug sensitive counterpart.
RNA-seq unveiled differential expression of 228 genes including a) genes previously identified as involved in CDR, b) genes not previously associated to the CDR phenotype, and c) novel transcripts whose function as a gene is uncharacterized. In particular, we show for the first time that CDR acquisition is correlated with an overexpression of the transcription factor encoding gene CZF1. CZF1 null mutants were susceptible to many drugs, independently of known multidrug resistance mechanisms. We show that CZF1 acts as a repressor of β-glucan synthesis, thus negatively regulating cell wall integrity. Finally, our RNA-seq data allowed us to identify a new transcribed region, upstream of the TAC1 gene, which encodes the major CDR transcriptional regulator.
Our results open new perspectives of the role of Czf1 and of our understanding of the transcriptional and post-transcriptional mechanisms that lead to the acquisition of drug resistance in C. albicans, with potential for future improvements of therapeutic strategies.
The yeast Candida albicans is the major cause of opportunistic fungal infections in humans. In case of systemic infections, the mortality rate can reach 50% . Azoles, which target the fungal P450 cytochrome 14alpha-lanosterol demethylase encoded by the ERG11 gene, are the most commonly used antifungal molecules for candidosis treatment . Unfortunately azoles only have a fungistatic effect and therefore have allowed the emergence of multidrug resistance strains in patients . There are two main groups of mutations which cause azole resistance in Candida albicans. The first group directly targets ERG11, either in cis by creating ERG11 alleles which encode a protein variant insensitive to azoles , or in trans by increasing the expression of ERG11 through gain of function mutations in the UPC2 gene, which encode a transcription factor regulating ERG11. The second type acts by increasing the expression of membrane transporters which export the drugs and therefore decrease their intracellular concentration . Two types of drug resistance, called MDR and CDR respectively, are distinguished, depending on the type of transporters which is involved. In most of the clinical isolates studied so far, the cells present either a drug resistance of the MDR or the CDR type [5, 6]. Only in rare occurrences, the two kind of resistance can be observed in the same cells [6–8]. The MDR resistant cells overexpress the Mdr1 permease encoding gene, which is under the control of the Mrr1 and Cap1 transcription factors . MDR1 is homologous to the Saccharomyces cerevisiae FLR1 transporter, which is involved in the detoxification of several drugs, including the antifungal benomyl [8, 10, 11]. The CDR resistant strains overexpress two ABC transporters encoding genes, CDR1 and CDR2, which are homologous to the pleiotropic drug resistance transporter Pdr5 from S. cerevisiae[12, 13]. These transporters are major determinants of the resistance to fluconazole, the most widely used azole. Their expression is controlled mainly by the Tac1 transcription factor . Most of the CDR resistant strains present a loss of heterozygocity and/or aneuploidy at the TAC1 locus, combined to gain of function mutations of TAC1[15, 16]. Microarray and proteomic analyses identified about 30 genes which are co-overexpressed with CDR1 and CDR2 in all CDR clinical isolates examined [17, 18]. Chromatine immunoprecipitation and tac1Δ strain analyses identified the TAC1 regulon, which is composed of about ten genes mostly involved in membrane properties, including CDR1 CDR2, the phospholipid transferase gene PDR16, the putative flippase gene RTA3 or the putative sphingosine kinase gene LCB4[18, 19].
The recently developed RNA sequencing (RNA seq) technology is revolutionizing our ability to analyse eukaryotic transcriptomes [20–22] and has given a chance to unravel the actual complexity of the CDR phenotype better. RNA seq has been proven to be more sensitive than microarrays . RNA seq also results in more accurate measurements of gene expression changes . Moreover, relative to standard microarray approaches, RNA seq and tiling arrays provide crude measurement and identification of transcripts, without any a priori on which region of the genome is transcribed [24, 25]. Hence, two recent RNA seq and one tiling array analyses of the Candida albicans transcriptome identified more than 1000 new transcripts, many of which are expressed in a condition-specific way [25–27]. Most of these transcripts do not have a coding potential and may be long non coding regulatory RNAs.
In this study, we have used RNA seq to analyse the transcriptomes of a CDR strain and its isogenic drug susceptible counterpart. In addition to the genes previously shown to be associated with CDR, we could identify about 50 genes which were overexpressed in the CDR strain. In particular, we show that the transcription factor encoding gene CZF1, which is involved in hyphal transition and white/opaque switching, is induced together with CDR1 and CDR2 in Gu5. Czf1 is likely to play an important role in CDR acquisition since its overexpression is a general feature of all the CDR strains that we have tested, but was not found in MDR strains. Moreover, its deletion caused susceptibility to several unrelated drugs. Additionally, the inactivation of CZF1 increased the resistance of the cells to cell wall perturbating agents, through the overexpression of beta glucan synthesis genes. We propose that Czf1 has a positive role on drug resistance and a negative role on cell wall integrity. Finally, we characterized a new transcribed region, previously undetected, just upstream of the TAC1 gene, which strongly suggests that TAC1 is subjected to complex post-transcriptional regulations, yet to be characterized. Taken together, our results open several new ways to our understanding of drug resistance acquisition in Candida albicans and may provide new targets for antifungal therapies.
Transcriptional landscape of Gu4 and Gu5 strains
In order to conduct a comparative analysis of the transcriptomes of CDR versus drug susceptible C. albicans cells, we performed high-throughput sequencing of cDNA made from poly(A) RNAs obtained from the Gu4 and Gu5 strains. Gu4 is a fluconazole susceptible clinical isolate obtained from an early infection episode. Gu5 is the corresponding fluconazole resistant clinical isolate obtained from later episode in the same patient treated with fluconazole. Gu4 and Gu5 are therefore supposed to be isogenic and mainly differ by their resistance to fluconazole. Gu5 has been characterised as a CDR strain, showing overexpression of CDR1 and CDR2, as compared with Gu4 . For each strain, two biologically independent samples of RNAs were analysed by RNA seq. We obtained 27 763 991 unique mappable reads for Gu4 (96% of total reads) and 20 570 932 mappable reads for Gu5 (97% of total reads). As noticed in previous RNA seq studies, the background was very low, with no or few reads mapped in intergenic regions. The correlation coefficient between the biological replicates was very high (0.991 for Gu5 and 0.999 for Gu4). We could detect expression for 5727 (92%) of the 6177 ORFs annotated in the CGD. The 450 ORFs which were undetected in our experiments are mostly dubious ORFs or gene without clear function (71%). We arbitrarily grouped the expressed ORFs in five categories: very high expression (more than 1000 reads/nucleotides), high expression (between 100 and 1000 reads/nucleotides), medium expression (between 10 and 100 reads/nucleotides), low expression (below 10 reads/nucleotides) and very low expression (below 3 reads/nucleotides). These categories contained 3%, 20%, 57%, 12% and 8% of the genes, respectively. As described previously, the “very high expression” category was enriched in genes encoding ribosomal proteins, glycolytic enzymes, translation regulators, histones and components of respiratory complexes. The “high expression” category was enriched in genes encoding enzymes involved in various metabolisms. The “medium” and “low” expression categories did not show any significant functional enrichment.
Finally, we used our data to detect new transcripts. To do so, the reads from the Gu4 and Gu5 strains were summed and used to identify regions of the genome with a minimum coverage of at least 3 reads/nucleotides (with at least one read in three independent samples) for a stretch of at least 50 nucleotides. With these criteria, we identified 1130 new transcription units, 800 of which partially or totally overlapped new transcribed regions recently identified by RNA seq [25, 26] or tiling arrays  ( Additional file 1: Table S2). We named these new transcripts NGuTs for “New Gu strains Transcripts”. The median size of NGuTs was 425 nucleotides. The majority of these NGuTs did not contain any obvious ORF and are therefore likely to be non-coding RNAs.
Transcripts differentially expressed in the drug resistant and the drug susceptible strains
The identification of genes differentially expressed in Gu4 and Gu5 was performed combining two different statistical methods (see methods section) and using the gene annotations from the CGD, the new transcribed regions defined by previous studies [25–27] and the new transcribed regions defined in this study (see above). We found 130 genes significantly overexpressed in Gu5 and 92 genes which were underexpressed in this strain compared to Gu4 (a full list and description of these genes can be found in Additional file 1: Table S3).
Comparison with previous microarray experiments
We first compared our results with the one obtained in a previous study, which used microarrays to analyse the transcriptomes of 4 different CDR isolates . Among the 130 genes induced in Gu5, 74 (56.5%) were previously shown to be overexpressed in one or several of these CDR clinical isolates ( Additional file 1: Table S3). Among these 74 genes, 53 (72%) were found overexpressed in at least 3 of the 4 CDR strains, suggesting that most of the genes identified in our RNA seq analyses are general features of the CDR isolates. These genes include the CDR1 and CDR2 ABC transporter encoding genes, together with the TAC1 gene encoding their transcriptional regulator. Remarkably, all the previously identified targets of Tac1p are found in our list, including PDR16 RTA3 and LCB4.
Of the 92 genes underexpressed in Gu5, 42 (46%) were previously shown to be repressed in CDR isolates. This includes genes involved in iron and metal homeostasis (FTR1, SIT1, FET34 for instance) and amino acid metabolism (MET3, OPT7, GAP4, SAM4 for instance).
ORFs not previously associated with CDR
We identified 53 ORFs which were overexpressed in Gu5 and were not previously associated with CDR phenotype. About half of them (45%) are uncharacterized and little information is available on their function.
Fifty genes were shown for the first time to be downregulated in Gu5. Of these, 24 (49%) are uncharacterized. Among the verified ORFs (Figure 2 and Additional file 1: Table S3), we could find the gene encoding the ABC transporter Cdr4, which is homologous to Cdr1 and Cdr2, and RTA4, which is homologous to the Tac1 target, RTA3. We also found in this group additional genes involved in copper and iron homeostasis (HEM13 CRP1 and SMF3) and genes encoding extracellular proteins (PLB1 ALS2 MAL2). Two genes encoding potential steroid targets, EBP1 and CBP1, were also underexpressed in Gu5, which is interesting considering the previously established link between steroid exposure and CDR response .
New transcripts differentially expressed in Gu5 and Gu4 strains
Among the new transcripts identified in this study, 5 were overexpressed and 8 underexpressed in Gu5 compared with Gu4. Four of the five induced transcripts were located immediately in 5’ of an induced ORF. Transcript NGuT-chr4-123 (previously identified as NOVEL-Ca21chr4-073 in ) was in 5’ of CZF1 (Figure 1), NGuT-chr5-037 was in 5’ of TAC1 (Figure 1), NGUT-chr2-038 was in 5’ of ADAEC (Figure 1) and NGuT-chr3-128 (previously identified as TF3W223 in ), was close to orf19.6713. The fifth transcript, NGuT-chr3-105 (NOVEL-Ca21chr3-056), is located just upstream of the AAF1 gene, which encodes an adhesin and which is not induced in the Gu5 strain.
All the new transcripts which were repressed in Gu5 are located at more than 1 kb from the adjacent ORFs. NGUT-chrR-157 (NOVEL-Ca21chrR-103) and 160 (TFRW269) belong to the same cluster of non coding RNAs. NGUT-chr4-069 (NOVEL-Ca21chr4-040) is located in a zone rich in tRNA genes. NGUT-chrR-041 (NOVEL-Ca21chrR-026) is close to a long terminal repeat.
Semi-quantitative RT-PCR analyses for 4 repressed and 3 induced NGuTs, using actin as a reference confirmed the RNA seq results, except for NGuT-chr3-105 and NGUT-chrR-157 (Figure 3).
New roles of the Czf1p transcription factor in multidrug resistance and cell wall maintenance
CZF1 deletion alters drug resistance
CZF1 deletion does not impair the expression of standard MDR and CDR actors
CZF1 negatively controls beta-glucan synthesis and resistance to cell wall perturbation agents
The sensitivity of the czf1 null mutant to various drugs prompted us to test the phenotype of this mutant in presence of drugs targeting the calcineurin pathway or cell wall integrity. Calcineurin has been shown to play multiple roles in stress responses including resistance to FLC and to membrane perturbation [30, 31]. It has also been shown to influence the dimorphic switch, a process which is controlled by Czf1 . We tested growth of wild type, czf1Δ/czf1Δ or czf1Δ/CZF1 strains in presence of the membrane perturbation agent SDS and the specific calcineurin inhibitor FK506. We observed that the double mutant was more sensitive to these drugs than the single mutant and the wild type strain (Figure 5C). We tested the expression of several direct or indirect actors of the calcineurin pathway (CMP1 CNB1 CRZ1 HSP90 and HOG1). We did not detect any effect of the partial or total inactivation of CZF1 on the expression of these genes (Figure 5B). Reciprocally, we tested the expression of CZF1 in various strains mutated for one or several actors of the calcineurin pathway. Again, we did not see any effects of the inactivation of the calcineurin pathway on the steady-state expression of CZF1 (data not shown).
In C. albicans, the cell wall plays important roles in virulence, defense against host and resistance to some drugs . Surprisingly, we observed that the czf1Δ/czf1Δ strain was more resistant than the wild type to congo red (a cell wall perturbation agent) (Figure 5C). The heterozygous mutant exhibited an intermediate phenotype (haploinsufficiency), suggesting a dosage effect for the role of CZF1 in cell wall maintenance. Congo red is considered to act on the cell wall glucan structure. Therefore, we measured the expression of the three beta 1,3-glucan synthase encoding genes (GSC1 GSL1 and GSL2) in presence and absence of CZF1 (Figure 5B). We observed no or few differences in GSC1 and GSL2 expression. In contrast, we observed a clear overexpression of GSL1 in the strains deleted for one or two copies of CZF1, compared to the wild type. This result suggests that Czf1 negatively controls the expression of this gene and that the resistance to congo red that we observed in the czf1 mutants could be linked to the partial deregulation of beta-glucan synthesis. Interestingly, a defect in beta-glucan activity has been shown to cause sensitivity to many different drugs, including azoles . However, the expression of GSL1 GSC1 and GSL2, tested by RNA-seq and semi-quantitative RT-PCR, was unchanged in the Gu5 strain as compared with Gu4 (data not shown).
CZF1 overexpression is commonly associated with CDR acquisition
Taken together these results strongly suggest that the overexpression of CZF1 is generally correlated to the acquisition of the CDR phenotype, with the notable exception of the DSY544/DSY775 pair.
In this work, we have used RNA-seq to analyze the transcriptome changes associated with CDR acquisition in a pair of clinical isolates in the human pathogen C. albicans. Our aim was to obtain an accuracy, sensitivity and completion unreached by previous microarray analyses of the same strains. Based on previous RNA-seq studies, we expected a better detection of mRNAs and of their expression changes, the identification of unannotated transcripts possibly up or down regulated in the CDR strain, and the detection of changes in the 5’ and 3’ UTR lengths compared with available annotations. Indeed, our data pointed out about one hundred genes, which were either induced or repressed in the Gu5 strain compared with Gu4, and which were not associated to CDR acquisition by previous studies. Although we identified more than one thousand new transcripts (some previously annotated by RNA seq or tiling arrays, some completely new), very few of them changed expression in the CDR strain. Four of these new transcripts were found to be induced in Gu5, but their position just upstream of annotated genes, which are similarly induced in Gu5, raises doubts on their actual existence as independent transcripts, or as long, previously undetected, 5’ UTR of the downstream gene.
The analyses of this large amount of data led to two important conclusions, which open new perspectives for our understanding of multidrug resistance in this pathogenic yeast. First, we revised the annotation of the TAC1 genomic region, identifying a new transcribed region in 5’ of the TAC1 mRNA. Second, we suggested new roles for the Czf1 transcription factor in CDR acquisition and cell wall maintenance.
A new transcribed region upstream of TAC1 and its potential impact on TAC1 expression regulation
RNA seq and tiling arrays experiments have allowed the systematic annotation of transcript boundaries in three yeast species: S. cerevisiae C. albicans and C. parapsilosis[21, 25, 27, 35]. In all three species, the 5’ and 3’UTR have been shown to be generally short, with an average size ranging between 50 and 80 nucleotides. However, in all these species, many exceptions have been described. In C. albicans, about 100 genes have 5’ UTR longer than 500 bp . Interestingly, this group is enriched in genes encoding transcriptional regulators and genes involved in filamentation. These include the CZF1 transcription factor encoding gene, which was induced in our experiments and which has a 5’UTR of about 2 kb ([27, 36]; this study). This observation is also true in mammals in which many genes with long 5’UTR encode regulators and proto-oncogenes . Long 5’UTR usually have important roles in the post-transcriptional and translational regulation of gene expression. In the yeast Saccharomyces cerevisiae, the expression of several transcription factors, including Gcn4 (amino acid starvation), Yap1 (oxidative stress response) and Yap2 (cadmium response), is regulated by the presence of short upstream ORFs (uORFs) in their long 5’UTR region, which exert negative or positive actions on the translation of the main ORF or on mRNA stability, depending on the cellular context . Recent studies have identified about 250 genes containing potential uORFs in this species .
New roles for Czf1 in CDR acquisition and cell wall maintenance
Czf1 is a zinc finger transcription factor mostly known for its role in the epigenetic white/opaque switch, which controls cell and colony morphogenesis and properties . Czf1 acts synergistically to Wor1 and Wor2 transcription factors, and antagonistically to Efg1, to promote the opaque state. Czf1 also has a role in hyphal growth, (also called filamentation). Czf1 mutants filament normally under standard filamentation induction conditions, but they are defective in hyphal development in anaerobic conditions . This phenotype is dependent on the composition of the media, for instance in terms of carbon source. The positive role of Czf1 on hypoxic filamentation is under the negative control of the TOR pathway and more specifically, of the Sch9 kinase . As for the white/opaque switch, Efg1 plays an important and complex role in filamentation. In aerobic conditions, it is a positive regulator of hyphal development, whereas in hypoxia, it counteracts the positive action of Czf1 on this process . Direct cross-talks exist between Efg1 and Czf1: Efg1 binds to the promoter of CZF1 and repress its expression, while Czf1 has a similar effect on the EFG1 promoter [36, 44]. More recently, Czf1 has been shown to be required for wild type adherence to silicone devices . Despite these roles in morphogenesis, the list of Czf1 targets has not been established to date. A partial transcriptome analysis of the czf1 mutant, conducted on 250 genes using the Nanochip technology, detected very few expression changes in this mutant .
Our data strongly suggest that Czf1 plays a role in the acquisition of multidrug resistance of the CDR type. First, CZF1 is overexpressed in all CDR mutants that we tested. Noteworthy, some other actors of white/opaque switch and hyphal growth were also overexpressed in the Gu5 strain, including SCH9, ADAEC, WH11 or SFL2. Second, the czf1 null mutant is sensitive to many unrelated drugs, including azoles and molecules targeting the calcineurin pathway. Remarkably, CZF1 does not seem to be involved in drug resistance phenotypes of the MDR kind, since its expression was unchanged in MDR mutant strains. The drug sensitivity associated with the inactivation of Czf1 does not seem to be due to standard drug resistance pathways, since the level of expression of the main actors of the CDR, MDR and calcineurin pathways were unchanged in this mutant. However, the fact that Czf1 does not control the steady-state expression levels of those genes does not mean that it is not involved in their induction following drug exposure. This hypothesis will be tested by inactivating CZF1 in the Gu5 and Gu4 backgrounds and testing the effect of this inactivation on drug resistance and gene expression.
Surprisingly, we observed that the czf1 mutant was resistant to congo red, a cell wall perturbation agent which targets beta-glucans, the most abundant component of Candida albicans cell wall. This phenotype was haploinsufficient, since the heterozygous mutant exhibited an intermediate resistance between the null mutant and the wild type. Congo red is thought to target the glucan architecture by directly interacting with glucan chains, therefore inhibiting synthesis of glucan polymers . We observed that the GSL1 gene encoding one of the three beta 1,3 glucan synthase subunits, was strongly overexpressed in the czf1 mutant. This suggests that Czf1 has a negative effect on the expression of some of the actors of the glucan synthesis pathway, which could explain the congo red resistance observed in the null mutant. Many functional connections have been established between cell wall components and hyphal growth. GSL1 itself has been shown to be down-regulated during filamentation , a process that is positively regulated by Czf1. Notably, Efg1 is also involved in cell wall architecture and efg1 mutants show a haploinsufficiency phenotype [47, 48]. The deletion of one or two copies of EFG1 leads to a decrease in the expression of many enzymes involved in cell wall biogenesis, including beta1-3 glucan synthases. This causes important changes in cell wall thickness and composition . Again, Czf1 and Efg1 seem to antagonize each other in the regulation of cell wall maintenance. The over-expression of some cell wall protein encoding genes has been extensively described as a compensatory mechanism in response to impairments in cell wall biosynthesis or to cell wall damages, both in S. cerevisiae and C. albicans (see for instance [49–51]). This phenomenon is largely controlled by the cell wall integrity (CWI) pathway  but part of it is independent of CWI . For instance, treatment with caspofungin (which directly inhibits the activity of the essential beta-glucanase subunit encoded by GSC1) triggers a compensatory process which leads to the overexpression of GSL1, similarly to what was observed in the czf1 null mutant . We suggest that Czf1 could play a role in this process.
All the media components were obtained from HiMedia (Mumbai, India). Amphotericin B (AMB), terbinafin, anisomycin, ketoconazole, sodium dodecyl sulphate, FK520, congo red were purchased from Sigma chemicals Co. (St. Louis, MO). Ranbaxy India generously provided FLC. Reverse transcriptase kit was purchased from fermentas, India. Caspofungin was generously provided by Dr. Neeraj Chauhan (Public Health Research Institute and Department of Microbiology and Molecular Genetics, New Jersey, USA). dATP, dGTP, dTTP and dCTP were bought from Sigma chemicals Co. (St. Louis, MO).
Yeast strains and growth conditions
Gu4 and Gu5 strains were described elsewhere . The strains of Candida albicans mutated in various genes (ALS1 MNN4 CZF1) that were used for phenotypic tests, were obtained from Dr. Ibrahim, A.S. (David Geffen School of Medicine at UCLA Los Angeles Biomedical Research Institute), Dr. Gow, N. (The university of Aberdeen, Scotland) and Dr. Kumamoto, C.A. (Tufts University, Boston, USA) respectively. The complete list of strains can be found in Additional file 1: Table S4. Cells were grown in YPD rich medium (2% glucose, 1% bactopeptone, 1% yeast extract) at 30°C for 48 hrs.
Twenty ml of exponentially growing (OD = 0.6-0.7) Gu4 or Gu5 cells were flash frozen in 30 ml of cold ethanol, centrifuged for 5 minutes at 3000 g, washed with sterile water and centrifuged again. The cell pellets were stored at −80°C. RNA were extracted using the RNeasy kit (Qiagen), with DNAse treatment. The RNA quality was checked on a bioanalyzer Nanochip (Agilent).
cDNA libraries and sequencing
Strand non-specific cDNA libraries (mRNA-Seq 8-Sample Prep Kit) were prepared according to manufacturer’s instructions (Illumina). The validation of libraries was made with Bioanalyser DNA 1000 chip (Agilent). All libraries were sequenced in single read mode, using a Genome Analyser instrument (Illumina).
Reads mapping and counting
Before read mapping, we trimmed polyN read tails, removed reads with a length lower or equals to 11 bases and discarded reads with quality mean lower or equals to 12. We aligned from the 4 samples a total of 48334923 bp single end reads produced from the Illumina GAII device using the SOAP2 software (version 2.20)  on the Candida albicans genome (assembly 21). The mapping tool was configured to return all alignment hits by setting the seed length to 28 bp. In addition we authorized with SOAP2 up to 5 mismatches with 2 in the seed.
Before expression estimation, alignments from reads that matches more than one time on the reference genome were removed. In addition, we searched for new transcripts by combining the Gu4 and Gu5 reads and identifying regions of the genome with a minimum coverage of at least 3 reads/nucleotides (with at least one read in three independent samples) for a stretch of at least 100 nucleotides. To compute gene expression abundance, we used Candida albicans GFF genome annotation from assembly 21 together with the annotation file for new transcripts and we counted all overlapping regions between alignments and listed exons. The complete results have been deposited at the GEO database (accession number: GSE38298).
Functional analyses of the genome-wide data were conducted using the CGD GO (gene ontology) term finder , with default parameters.
Genes differentially expressed between Gu5 and Gu4 strains were identified applying two different methodologies available in the R programming language, i.e. DESeq  and edgeR [59, 60]. Default parameters were used and genes with p-value lower than 0.01 with the two methods were finally selected for further analysis ( Additional file 1: Table S3).
Final point semi-quantitative RT-PCR analyses
RT-PCR was done using the RevertAidTM H Minus kit (MBI, Fermentas). Briefly, 1 μg isolated RNA was primed with oligo (dT)18 for cDNA synthesis at 42°C for 60 min. Reverse transcription reaction was terminated by heating at 70°C for 5 min. The synthesized cDNA product (2 μl) was directly used for PCR amplification reaction (50 μl) using gene specific forward and reverse primers. The amplified products were gel electrophoresed and quantitated by using Quantity One software by Bio Rad gel documentation system. The RNA samples used for these validation experiments were independent from the samples used for RNA seq. The sequences of the primers used for PCR reactions can be found in Additional file 1: Table S5.
Drug-susceptibility was tested by broth-microdilution assay according to CLSI (Clinical and Laboratory Standards Institute) and serial dilution assay essentially as described previously (CLSI, 2008). Following drugs stocks (in parenthesis) and working concentration were used 2.5 μg/ml of FLC (5 mg/ml), 4 μg/ml of TER (4 mg/ml), 4 μg/ml of AMB (4 mg/ml), 150 μg/ml of ANS (20 mg/ml), 0.5 μg/ml of KTC (1 mg/ml), 50 ng/ml of CAS (1 mg/ml), 0.08% SDS, 200 μg/ml of CR (10 mg/ml), and 200 μg/ml FK520 (20 mg/ml).
The work presented in this paper has been supported in part by grants to R.P. from Indo-French Center for the Promotion of Advanced Scientific Research (IFC-3403-2/2006) and DST INT/SWISS/P-31/2009. SD is thankful to Council of Scientific and Industrial Research, Government of India, for awarding Senior Research Fellowship. We thank Claire Torchet for useful discussions regarding long 5’UTR in eukaryotes.
- Pfaller MA, Diekema DJ, Gibbs DL, Newell VA, Meis JF, Gould IM, Fu W, Colombo AL, Rodriguez-Noriega E: Global Antifungal Surveillance Study: Results from the ARTEMIS DISK Global Antifungal Surveillance study, 1997 to 2005: an 8.5-year analysis of susceptibilities ofCandidaspecies and other yeast species to fluconazole and voriconazole determined by CLSI standardized disk diffusion testing. J Clin Microbiol. 2007, 45: 1735-1745. 10.1128/JCM.00409-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Odds FC, Brown AJ, Gow NA: Antifungal agents: mechanisms of action. Trends Microbiol. 2003, 11: 272-279. 10.1016/S0966-842X(03)00117-3.View ArticlePubMedGoogle Scholar
- Anderson JB: Evolution of antifungal-drug resistance: mechanisms and pathogen fitness. Nat Rev Microbiol. 2005, 3: 547-556. 10.1038/nrmicro1179.View ArticlePubMedGoogle Scholar
- Akins RA: An update on antifungal targets and mechanisms of resistance in Candida albicans. Med Mycol. 2005, 43: 285-318. 10.1080/13693780500138971.View ArticlePubMedGoogle Scholar
- Barker KS, Rogers PD: Recent insights into the mechanisms of antifungal resistance. Curr Infect Dis Rep. 2006, 8: 449-456. 10.1007/s11908-006-0019-3.View ArticlePubMedGoogle Scholar
- Lopez-Ribot JL, McAtee RK, Lee LN, Kirkpatrick WR, White TC, Sanglard D, Patterson TF: Distinct patterns of gene expression associated with development of fluconazole resistance in serial Candida albicans isolates from human immunodeficiency virus-infected patients with oropharyngeal candidiasis. Antimicrob Agents Chemother. 1998, 42: 2932-2937.PubMed CentralPubMedGoogle Scholar
- Saidane S, Weber S, De Deken X, St-Germain G, Raymond M: PDR16-mediated azole resistance in Candida albicans. Mol Microbiol. 2006, 60: 1546-1562. 10.1111/j.1365-2958.2006.05196.x.View ArticlePubMedGoogle Scholar
- Sanglard D, Kuchler K, Ischer F, Pagani JL, Monod M, Bille J: Mechanisms of resistance to azole antifungal agents in Candida albicans isolates from AIDS patients involve specific multidrug transporters. Antimicrob Agents Chemother. 1995, 39: 2378-2386. 10.1128/AAC.39.11.2378.PubMed CentralView ArticlePubMedGoogle Scholar
- Schubert S, Barker KS, Znaidi S, Schneider S, Dierolf F, Dunkel N, Aïd M, Boucher G, Rogers PD, Raymond M, Morschhäuser J: Regulation of efflux pump expression and drug resistance by the transcription factors Mrr1, Upc2, and Cap1 in Candida albicans. Antimicrob Agents Chemother. 2011, 55: 2212-2223. 10.1128/AAC.01343-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Ben-Yaacov R, Knoller S, Caldwell GA, Becker JM, Koltin Y: Candida albicans gene encoding resistance to benomyl and methotrexate is a multidrug resistance gene. Antimicrob Agents Chemother. 1994, 38: 648-652. 10.1128/AAC.38.4.648.PubMed CentralView ArticlePubMedGoogle Scholar
- Lucau-Danila A, Lelandais G, Kozovska Z, Tanty V, Delaveau T, Devaux F, Jacq C: Early expression of yeast genes affected by chemical stress. Mol Cell Biol. 2005, 25: 1860-1868. 10.1128/MCB.25.5.1860-1868.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Prasad R, De Wergifosse P, Goffeau A, Balzi E: Molecular cloning and characterization of a novel gene of Candida albicans, CDR1, conferring multiple resistance to drugs and antifungals. Curr Genet. 1995, 27: 320-329. 10.1007/BF00352101.View ArticlePubMedGoogle Scholar
- Sanglard D, Ischer F, Monod M, Bille J: Cloning of Candida albicans genes conferring resistance to azole antifungal agents: characterization of CDR2, a new multidrug ABC transporter gene. Microbiology. 1997, 143: 405-416. 10.1099/00221287-143-2-405.View ArticlePubMedGoogle Scholar
- Coste AT, Karababa M, Ischer F, Bille J, Sanglard D: TAC1, transcriptional activator of CDR genes, is a new transcription factor involved in the regulation of Candida albicans ABC transporters CDR1 and CDR2. Eukaryot Cell. 2004, 3: 1639-1652. 10.1128/EC.3.6.1639-1652.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Selmecki A, Gerami-Nejad M, Paulson C, Forche A, Berman J: An isochromosome confers drug resistance in vivo by amplification of two genes, ERG11 and TAC1. Mol Microbiol. 2008, 68: 624-641. 10.1111/j.1365-2958.2008.06176.x.View ArticlePubMedGoogle Scholar
- Coste A, Selmecki A, Forche A, Diogo D, Bougnoux ME, d'Enfert C, Berman J, Sanglard D: Genotypic evolution of azole resistance mechanisms in sequential Candida albicans isolates. Eukaryot Cell. 2007, 6: 1889-1904. 10.1128/EC.00151-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Hoehamer CF, Cummings ED, Hilliard GM, Morschhäuser J, Rogers PD: Proteomic analysis of Mrr1p- and Tac1p-associated differential protein expression in azole-resistant clinical isolates of Candida albicans. Proteomics Clin Appl. 2009, 3: 968-978. 10.1002/prca.200800252.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu TT, Znaidi S, Barker KS, Xu L, Homayouni R, Saidane S, Morschhäuser J, Nantel A, Raymond M, Rogers PD: Genome-wide expression and location analyses of the Candida albicans Tac1p regulon. Eukaryot Cell. 2007, 6: 2122-2138. 10.1128/EC.00327-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Znaidi S, De Deken X, Weber S, Rigby T, Nantel A, Raymond M: The zinc cluster transcription factor Tac1p regulates PDR16 expression in Candida albicans. Mol Microbiol. 2006, 66: 440-452.View ArticleGoogle Scholar
- Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009, 10: 57-63. 10.1038/nrg2484.PubMed CentralView ArticlePubMedGoogle Scholar
- Nagalakshmi U, Wang Z, Waern K, Shou C, Raha D, Gerstein M, Snyder M: The transcriptional landscape of the yeast genome defined by RNA sequencing. Science. 2008, 320: 1344-1349. 10.1126/science.1158441.PubMed CentralView ArticlePubMedGoogle Scholar
- Sultan M, Schulz MH, Richard H, Magen A, Klingenhoff A, Scherf M, Seifert M, Borodina T, Soldatov A, Parkhomchuk D, Schmidt D, O'Keeffe S, Haas S, Vingron M, Lehrach H, Yaspo ML: A global view of gene activity and alternative splicing by deep sequencing of the human transcriptome. Science. 2008, 321: 956-960. 10.1126/science.1160342.View ArticlePubMedGoogle Scholar
- Xiong Y, Chen X, Chen Z, Wang X, Shi S, Wang X, Zhang J, He X: RNA sequencing shows no dosage compensation of the active X-chromosome. Nat Genet. 2010, 42: 1043-1047. 10.1038/ng.711.View ArticlePubMedGoogle Scholar
- Guttman M, Garber M, Levin JZ, Donaghey J, Robinson J, Adiconis X, Fan L, Koziol MJ, Gnirke A, Nusbaum C, Rinn JL, Lander ES, Regev A: Ab initio reconstruction of cell type-specific transcriptomes in mouse reveals the conserved multi-exonic structure of lincRNAs. Nat Biotechnol. 2010, 28: 503-510. 10.1038/nbt.1633.PubMed CentralView ArticlePubMedGoogle Scholar
- Bruno VM, Wang Z, Marjani SL, Euskirchen GM, Martin J, Sherlock G, Snyder M: Comprehensive annotation of the transcriptome of the human fungal pathogen Candida albicans using RNA-seq. Genome Res. 2010, 20: 1451-1458. 10.1101/gr.109553.110.PubMed CentralView ArticlePubMedGoogle Scholar
- Tuch BB, Mitrovich QM, Homann OR, Hernday AD, Monighetti CK, De La Vega FM, Johnson AD: The transcriptomes of two heritable cell types illuminate the circuit governing their differentiation. PLoS Genet. 2010, 6: e1001070-10.1371/journal.pgen.1001070.PubMed CentralView ArticlePubMedGoogle Scholar
- Sellam A, Hogues H, Askew C, Tebbji F, van Het Hoog M, Lavoie H, Kumamoto CA, Whiteway M, Nantel A: Experimental annotation of the human pathogen Candida albicans coding and noncoding transcribed regions using high-resolution tiling arrays. Genome Biol. 2010, 11: R71-10.1186/gb-2010-11-7-r71.PubMed CentralView ArticlePubMedGoogle Scholar
- Franz R, Ruhnke M, Morschhäuser J: Molecular aspects of fluconazole resistance development in Candida albicans. Mycoses. 1999, 42: 453-458. 10.1046/j.1439-0507.1999.00498.x.View ArticlePubMedGoogle Scholar
- Banerjee D, Lelandais G, Shukla S, Mukhopadhyay G, Jacq C, Devaux F, Prasad R: Responses of pathogenic and nonpathogenic yeast species to steroids reveal the functioning and evolution of multidrug resistance transcriptional networks. Eukaryot Cell. 2008, 7: 68-77. 10.1128/EC.00256-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Cruz MC, Goldstein AL, Blankenship JR, Del Poeta M, Davis D, Cardenas ME, Perfect JR, McCusker JH, Heitman J: Calcineurin is essential for survival during membrane stress in Candida albicans. EMBO J. 2002, 21: 546-559. 10.1093/emboj/21.4.546.PubMed CentralView ArticlePubMedGoogle Scholar
- Sanglard D, Ischer F, Marchetti O, Entenza J, Bille J: Calcineurin A of Candida albicans: involvement in antifungal tolerance, cell morphogenesis and virulence. Mol Microbiol. 2003, 48: 959-976. 10.1046/j.1365-2958.2003.03495.x.View ArticlePubMedGoogle Scholar
- Klis FM, Sosinska GJ, de Groot PW, Brul S: Covalently linked cell wall proteins of Candida albicans and their role in fitness and virulence. FEMS Yeast Res. 2009, 9: 1013-1028. 10.1111/j.1567-1364.2009.00541.x.View ArticlePubMedGoogle Scholar
- Cantón E, Pemán J, Gobernado M, Alvarez E, Baquero F, Cisterna R, Gil J, Martín-Mazuelos E, Rubio C, Sánchez-Sousa A, Serrano C: Sensititre YeastOne caspofungin susceptibility testing of Candida clinical isolates: correlation with results of NCCLS M27-A2 multicenter study. Antimicrob Agents Chemother. 2005, 49: 1604-1607. 10.1128/AAC.49.4.1604-1607.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- White TC: Increased mRNA levels of ERG16, CDR, and MDR1 correlate with increases in azole resistance in Candida albicans isolates from a patient infected with human immunodeficiency virus. Antimicrob Agents Chemother. 1997, 41: 1482-1487.PubMed CentralPubMedGoogle Scholar
- Guida A, Lindstädt C, Maguire SL, Ding C, Higgins DG, Corton NJ, Berriman M, Butler G: Using RNA-seq to determine the transcriptional landscape and the hypoxic response of the pathogenic yeast Candida parapsilosis. BMC Genomics. 2011, 12: 628-10.1186/1471-2164-12-628.PubMed CentralView ArticlePubMedGoogle Scholar
- Vinces MD, Haas C, Kumamoto CA: Expression of the Candida albicans morphogenesis regulator gene CZF1 and its regulation by Efg1p and Czf1p. Eukaryot Cell. 2006, 5: 825-835. 10.1128/EC.5.5.825-835.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- Pickering BM, Willis AE: The implications of structured 5' untranslated regions on translation and disease. Semin Cell Dev Biol. 2005, 16: 39-47. 10.1016/j.semcdb.2004.11.006.View ArticlePubMedGoogle Scholar
- Vilela C, McCarthy JE: Regulation of fungal gene expression via short open reading frames in the mRNA 5'untranslated region. Mol Microbiol. 2003, 49: 859-867. 10.1046/j.1365-2958.2003.03622.x.View ArticlePubMedGoogle Scholar
- Lachke SA, Srikantha T, Soll DR: The regulation of EFG1 in white-opaque switching in Candida albicans involves overlapping promoters. Mol Microbiol. 2003, 48: 523-536. 10.1046/j.1365-2958.2003.t01-1-03448.x.View ArticlePubMedGoogle Scholar
- Lohse MB, Johnson AD: White-opaque switching in Candida albicans. Curr Opin Microbiol. 2009, 12: 650-654. 10.1016/j.mib.2009.09.010.PubMed CentralView ArticlePubMedGoogle Scholar
- Brown DH, Giusani AD, Chen X, Kumamoto CA: Filamentous growth of Candida albicans in response to physical environmental cues and its regulation by the unique CZF1 gene. Mol Microbiol. 1999, 34: 651-662. 10.1046/j.1365-2958.1999.01619.x.View ArticlePubMedGoogle Scholar
- Stichternoth C, Fraund A, Setiadi E, Giasson L, Vecchiarelli A, Ernst JF: Sch9 kinase integrates hypoxia and CO2 sensing to suppress hyphal morphogenesis in Candida albicans. Eukaryot Cell. 2011, 10: 502-511. 10.1128/EC.00289-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Giusani AD, Vinces M, Kumamoto CA: Invasive filamentous growth of Candida albicans is promoted by Czf1p-dependent relief of Efg1p-mediated repression. Genetics. 2002, 160: 1749-1753.PubMed CentralPubMedGoogle Scholar
- Lassak T, Schneider E, Bussmann M, Kurtz D, Manak JR, Srikantha T, Soll DR, Ernst JF: Target specificity of the Candida albicans Efg1 regulator. Mol Microbiol. 2011, 82: 602-618. 10.1111/j.1365-2958.2011.07837.x.View ArticlePubMedGoogle Scholar
- Finkel JS, Xu W, Huang D, Hill EM, Desai JV, Woolford CA, Nett JE, Taff H, Norice CT, Andes DR, Lanni F, Mitchell AP: Portrait of Candida albicans adherence regulators. PLoS Pathog. 2012, 8: e1002525-10.1371/journal.ppat.1002525.PubMed CentralView ArticlePubMedGoogle Scholar
- Kopecká M, Gabriel M: The influence of congo red on the cell wall and (1–3)-beta-D-glucan microfibril biogenesis in Saccharomyces cerevisiae. Arch Microbiol. 1992, 158: 115-126. 10.1007/BF00245214.View ArticlePubMedGoogle Scholar
- Sohn K, Urban C, Brunner H, Rupp S: EFG1 is a major regulator of cell wall dynamics in Candida albicans as revealed by DNA microarrays. Mol Microbiol. 2003, 47: 89-102.View ArticlePubMedGoogle Scholar
- Zavrel M, Majer O, Kuchler K, Rupp S: Transcription factor Efg1 shows a haploinsufficiency phenotype in modulating the cell wall architecture and immunogenicity of Candida albicans. Eukaryot Cell. 2012, 11: 129-140. 10.1128/EC.05206-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Agarwal AK, Rogers PD, Baerson SR, Jacob MR, Barker KS, Cleary JD, Walker LA, Nagle DG, Clark AM: Genome-wide expression profiling of the response to polyene, pyrimidine, azole, and echinocandin antifungal agents in Saccharomyces cerevisiae. J Biol Chem. 2003, 278: 34998-35015. 10.1074/jbc.M306291200.View ArticlePubMedGoogle Scholar
- Lagorce A, Hauser NC, Labourdette D, Rodriguez C, Martin-Yken H, Arroyo J, Hoheisel JD, François J: Genome-wide analysis of the response to cell wall mutations in the yeast Saccharomyces cerevisiae. J Biol Chem. 2003, 278: 20345-20357. 10.1074/jbc.M211604200.View ArticlePubMedGoogle Scholar
- Liu TT, Lee RE, Barker KS, Lee RE, Wei L, Homayouni R, Rogers PD: Genome-wide expression profiling of the response to azole, polyene, echinocandin, and pyrimidine antifungal agents in Candida albicans. Antimicrob Agents Chemother. 2005, 49: 2226-2236. 10.1128/AAC.49.6.2226-2236.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- García R, Rodríguez-Peña JM, Bermejo C, Nombela C, Arroyo J: The high osmotic response and cell wall integrity pathways cooperate to regulate transcriptional responses to zymolyase-induced cell wall stress in Saccharomyces cerevisiae. J Biol Chem. 2009, 284: 10901-10911. 10.1074/jbc.M808693200.PubMed CentralView ArticlePubMedGoogle Scholar
- Arias P, Díez-Muñiz S, García R, Nombela C, Rodríguez-Peña JM, Arroyo J: Genome-wide survey of yeast mutations leading to activation of the yeast cell integrity MAPK pathway: novel insights into diverse MAPK outcomes. BMC Genomics. 2011, 12: 390-10.1186/1471-2164-12-390.PubMed CentralView ArticlePubMedGoogle Scholar
- Douglas CM, D'Ippolito JA, Shei GJ, Meinz M, Onishi J, Marrinan JA, Li W, Abruzzo GK, Flattery A, Bartizal K, Mitchell A, Kurtz MB: Identification of the FKS1 gene of Candida albicans as the essential target of 1,3-beta-D-glucan synthase inhibitors. Antimicrob Agents Chemother. 1997, 41: 2471-2479.PubMed CentralPubMedGoogle Scholar
- Park S, Kelly R, Kahn JN, Robles J, Hsu MJ, Register E, Li W, Vyas V, Fan H, Abruzzo G, Flattery A, Gill C, Chrebet G, Parent SA, Kurtz M, Teppler H, Douglas CM, Perlin DS: Specific substitutions in the echinocandin target Fks1p account for reduced susceptibility of rare laboratory and clinical Candida sp. isolates. Antimicrob Agents Chemother. 2005, 49: 3264-3273. 10.1128/AAC.49.8.3264-3273.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Li R, Yu C, Li Y, Lam TW, Yiu SM, Kristiansen K, Wang J: SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009, 25: 1966-1967. 10.1093/bioinformatics/btp336.View ArticlePubMedGoogle Scholar
- Boyle EI, Weng S, Gollub J, Jin H, Botstein D, Cherry JM, Sherlock G: GO:: TermFinder–open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. Bioinformatics. 2004, 20: 3710-3715. 10.1093/bioinformatics/bth456.PubMed CentralView ArticlePubMedGoogle Scholar
- Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11: R106-10.1186/gb-2010-11-10-r106.PubMed CentralView ArticlePubMedGoogle Scholar
- Robinson MD, McCarthy DJ, Smyth GK: edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010, 26: 139-140. 10.1093/bioinformatics/btp616.PubMed CentralView ArticlePubMedGoogle Scholar
- Thorvaldsdóttir H, Robinson JT, Mesirov JP: Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2012Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.