Skip to main content
Download PDF

Genome-wide expression profiling of the response to short-term exposure to fluconazole in Cryptococcus neoformans serotype A

BMC Microbiology volume 11, Article number: 97 (2011) Cite this article

Abstract

Background

Fluconazole (FLC), a triazole antifungal drug, is widely used for the maintenance therapy of cryptococcal meningoencephalitis, the most common opportunistic infection in AIDS patients. In this study, we examined changes in the gene expression profile of the C. neoformans reference strain H99 (serotype A) following FLC treatment in order to investigate the adaptive cellular responses to drug stress.

Results

Simultaneous analysis of over 6823 transcripts revealed that 476 genes were responsive to FLC. As expected up-regulation of genes involved in ergosterol biosynthesis was observed, including the azole target gene ERG11 and ERG13, ERG1, ERG7, ERG25, ERG2, ERG3 and ERG5. In addition, SRE1 which is a gene encoding a well-known regulator of sterol homeostasis in C. neoformans was up-regulated. Several other genes such as those involved in a variety of important cellular processes (i.e. lipid and fatty acid metabolism, cell wall maintenance, stress and virulence) were found to be up-regulated in response to FLC treatment. Conversely, expression of AFR1, the major transporter of azoles in C. neoformans, was not regulated by FLC.

Conclusions

Short-term exposure of C. neoformans to FLC resulted in a complex altered gene expression profile. Some of the observed changes could represent specific adaptive responses to the antifungal agent in this pathogenic yeast.

Background

Cryptococcus neoformans is a basidiomycetous fungal pathogen that causes meningoencephalitis in predominantly immunocompromised hosts [1, 2], that is the most devastating manifestation of cryptococcal disease and is fatal unless treated [3]. Cryptococcosis appears to be a significant opportunistic infection in solid-organ transplant recipients, with a prevalence rate ranging from 0.26% to 5% and overall mortality of 42% [4]. Notably, cryptococcal meningitis was reported to occur in 46% of patients from an Indian HIV-positive cohort [5]. Although the introduction of highly active antiretroviral therapy has led to a decrease in the number of cryptococcal infections in AIDS patients in most developed countries, this is not the case in developing countries where the incidence of HIV/AIDS and cryptococcal meningitis continue to rise [6]. As fluconazole (FLC) became increasingly used due to the need for life-long maintenance therapy in HIV/AIDS patients, FLC resistance was hence detected at relatively high frequency in C. neoformans clinical isolates from India, Africa and Cambodia [7–9].

Increased FLC resistance in vitro was shown to be predictive of treatment failures and infection relapses [10]. Recently, the mechanism underlying the heteroresistance to FLC was elucidated [11], that is an adaptive mode of azole resistance previously associated with FLC therapy failure cases [12]. This mechanism is based on duplications of multiple chromosomes in response to drug pressure [13]. Interestingly, Sionov et al. [13] observed that the number of disomic chromosomes positively correlated with the duration of exposure to FLC, whereas the duplication of chromosome 1 was closely associated with two genes, ERG11, the target of FLC [14], and AFR1, the major transporter of azoles in C. neoformans [11, 15]. Such genomic plasticity enables cells to cope with drug stress and was observed in C. neoformans strains of both serotypes, A (C. neoformans var. grubii) and D (C. neoformans var. neoformans) [13].

The recent sequencing of the C. neoformans genome [16] has stimulated the development of C. neoformans-specific microarrays that made possible to address hypotheses about global responses to overcome stresses during growth in the human host [17, 18]. Regardless of the source (i.e. host-derived or antifungal drugs), toxic compounds exert constant selective pressure on the fungus that responds by developing mechanisms necessary for survival [19].

With the aim to identify genes required for adaptive growth in the presence of sub-inhibitory concentrations of FLC, we investigated here the transient response of C. neoformans to FLC by analyzing differences in gene expression prior and after FLC exposure of strain H99, a reference strain of serotype A. Thus, genome-wide transcriptional profiling of over 6823 C. neoformans genes identified 476 genes with significant expression changes. Apart from genes involved in ergosterol biosynthesis (e.g. ERG11), genes involved in other important cellular functions, such as those encoding the sterol homeostasis regulator Sre1 [20] or phospholipase B1 (Plb1) [21], were shown to be induced by FLC treatment. In addition, AFR1 was not found FLC-responsive, suggesting indirectly that this gene is responsible for long-term FLC adaptation in C. neoformans.

Methods

Strain, growth conditions and RNA isolation

C. neoformans var. grubii serotype A strain (H99) was obtained from David S. Perlin [22], kept as 20% glycerol stock at -80°C and sub-cultured, as required, on YEPD (1% yeast extract, 2% peptone, 2% glucose) agar plates at 30°C. For RNA isolation independent overnight cultures were diluted 1:100 in liquid YEPD and grown at 30°C or 37°C with agitation for 3 h to reach a density of 3 × 107 CFU/ml. At this point cultures were equally divided into two aliquots to which either FLC at a concentration of 10 mg/l or distilled water was added, followed by incubation at 30°C or 37°C for 90 min. After this treatment, cultures were centrifuged at 4°C and 5500 × g and total RNA was extracted as previously described [23].

Microarray design and preparation

C. neoformans H99 microarrays were designed following the Agilent Array Design guidelines (Earray platform) by first creating two separate sets of 60-base nucleotide probes for each of 6967 open reading frame (ORF) sequences as downloaded from the Broad Institute website http://www.broadinstitute.org/annotation/genome/cryptococcusneoformans/MultiHome.html. The probe selection was performed using the GE Probe Design Tool; probes were filtered following their base composition and distribution, cross-hybridization potential, and melting temperature, to yield final duplicate probes representing 6823 ORFs to cover 97.9% of the whole C. neoformans H99 genome. C. neoformans custom arrays were manufactured in the 8 × 15k format by Agilent Technologies (Santa Clara, CA, USA). For quality control and normalization purposes, 157 probes were selected randomly and spotted 10 times throughout each array. Standard controls (Agilent Technologies) were also included.

cRNA synthesis, labeling and hybridization

RNA sample preparation was performed on three biological triplicates of H99 cells grown at 30°C, as described above. Prior to the labeling/amplification step, purity and integrity of the RNA samples were determined using Agilent RNA 6000 Nano LabChip kit on the Agilent 2100 bioanalyzer (Agilent Technologies). Agilent's One-Color Quick Amp Labeling kit (Agilent Technologies) was used to generate fluorescently labeled cRNA probes according to the manufacturer's instructions. The method uses T7 RNA polymerase, which simultaneously amplifies target material and incorporates cyanine 3-labeled CTP. The labeled cRNAs were purified with the RNeasy Mini kit (Qiagen, Hilden, Germany) and quantified using NanoDrop ND-1000 UV-VIS spectrophotometer. Aliquots (600 ng) of Cy3-labeled cRNAs were fragmented and hybridized for 17 h at 65°C to each array using the Gene Expression Hybridization kit (Agilent Technologies) and according to the manufacturer's instructions.

Microarray imaging and data analysis

Slides were washed and processed according to the Agilent 60-mer Oligo Microarray Processing protocol and scanned on a Agilent microarray scanner G2565BA (Agilent Technologies). Data were extracted from the images with Feature Extraction (FE) software (Agilent Technologies). FE software flags outlier features, and detects and removes spatial gradients and local backgrounds. Data were normalized using a combined rank consistency filtering with LOWESS intensity normalization. The gene expression values obtained from FE software were imported into GeneSpring 10.0.2 software (Agilent Technologies) for pre-processing and data analysis. For inter-array comparisons, a linear scaling of the data was performed using the 75th percentile signal value of all of non-control probes on the microarray to normalize one-colour signal values. Probe sets with a signal intensity value below the 20th percentile were considered as absent and discarded from subsequent analysis. The expression of each gene was normalized by its median expression across all samples. Genes were included in the final data set if their expression changed by at least twofold between strain H99 FLC-exposed or -not exposed (control sample) in at least two independent experiments, together with a P-value cut-off of < 0.05 (by one-way analysis of variance [ANOVA] corrected). Genes listed in Table 1 were categorized by reported or putative functions by the BROAD Institute database with NCBI blastP http://www.ncbi.nlm.nih.gov/BLAST/ editing, and also by the Uniprot http://www.uniprot.org/ and Saccharomyces genome http://www.yeastgenome.org/cgi-bin/blast-sgd.pl databases. As indicated in Table 1, each S. cerevisiae gene name was assigned by blastP search with the C. neoformans H99 gene sequence (e-value cutoff: e-6) according to Kim et al. [24]. Gene Ontology (GO) term analysis was carried to help categorize a list of genes into functional groups. The whole microarray data have been deposited in National Center for Biotechnology Information's Gene Expression Omnibus [25] and are accessible through GEO Series accession number GSE24927.

Table 1 Changes in the gene expression of C. neoformans H99 cells exposed to FLC
Full size table

Quantitative RT-PCR (qRT-PCR) validation of gene expression

Expression of selected differentially regulated genes as identified by the microarray analysis was quantitatively assessed with qRT-PCR in an i-Cycler iQ system (Bio-Rad Laboratories, Hercules, CA, USA). All primers and probes (see Additional file 1) were designed with Beacon Designer 2 (version 2.06) software (Premier Biosoft International, Palo Alto, CA, USA) and synthesized by MWG Biotech (Florence, Italy). qRT-PCRs were carried out as previously described [23]. The annealing temperature used for all primers was 65°C. Each reaction was run in triplicate on three separate occasions. For relative quantification of target gene expression, ACT1 was used as a normalizer gene [23]. Changes (n-fold) in gene expression relative to that of the control were determined from mean ACT1-normalized expression levels.

Oxidative stress and cell wall inhibitor assays

Susceptibilities to hydrogen peroxide (H2O2) and cell wall inhibitors were measured with exponentially growing cells in liquid YEPD at 30°C or 37°C pre-treated or not with FLC (10 mg/l) for 90 min as described elsewhere with modifications [26, 27]. The cells were next washed with sterile PBS and diluted to an OD650 of 1.0 in PBS. For the oxidative stress assays, aliquots of the cell suspensions were transferred to Eppendorf tubes where H2O2 (Sigma, Milan, Italy) was added to 20 mM and incubated at 30°C or 37°C for 2 h. Viability was determined after appropriate dilution of the samples with PBS by plating 100 μl in triplicate on solid YEPD. The CFU were counted after incubation for 72 h at 30°C or 37°C. For the cell wall inhibitor assays, dilutions of the cell suspensions were made in PBS and 5 μl of these were grown on YEPD plates containing 0.5% Congo red (Sigma, C-6767), 0.5, 1.0 and 1.5 mg ml-1 calcofluor white (Sigma, F-3543), 0.01%, 0.03% and 0.06% SDS (Sigma) and 0.2, 0.5 and 1.0 mg ml-1 caffeine (Sigma, C-0750). Plates were incubated for 48 h at 30°C or 37°C and photographed.

Results and Discussion

Experimental design and global gene expression results

The transcript profiles of C. neoformans H99 cells exposed to 10 mg/l of FLC (1/2 × MIC) for one doubling time (90 min) at 30°C were compared with profiles of untreated cells. A total of 476 genes were found responsive to FLC treatment under the test conditions, consisting of a single concentration and a single time point as described elsewhere [28–30]. The threshold value used in the present analysis was at least a twofold difference of gene expression between the experimental conditions, which is a value generally accepted in fungal genome-wide expression profiling [31]. Given that approximately 95% of the genes (6434/6823) spotted on the microarrays gave validated data, the above mentioned number indicate that 7.4% of the total number of genes in the C. neoformans H99 genome exhibited transcriptional changes, with 231 genes being upregulated and 245 downregulated upon FLC treatment.

In order to verify the changes in gene expression identified by our microarray analysis, we randomly selected 10 target genes (CNAG_00747, CNAG_01858, CNAG_02048, CNAG_02226, CNAG_03007, CNAG_03204, CNAG_04632, CNAG_03433, CNAG_05264, CNAG_05602) including those regulated and not regulated by FLC for validation of microarray data. A strong correlation (r = 0.94) was found between relative expression levels obtained by microarray or qRT-PCR analysis (Figure 1). In addition, qRT-PCR experiments performed with RNA extracted from H99 cells FLC-treated at 37°C demonstrated that expression of the target genes also including AFR1 was comparable to that obtained when H99 cells were pre-treated with FLC at 30°C (Figure 2).

Figure 1
figure 1

Scatter plot of the results by microarray and quantitative RT-PCR analyses for ten selected differentially regulated genes in H99 cells FLC-treated (H99F) compared to untreated control cells.

Full size image
Figure 2
figure 2

Results of qRT-PCR analysis performed with RNAs extracted from H99 cells FLC-treated (H99F) at 30°C and 37°C. The values, which are means of three separated experiments, represent the increase in gene expression relative to untreated control cells (set at 1.00). Error bars show standard deviations

Full size image

The genes listed in Table 1 were categorized in 10 main groups by functional profiles as described in Methods. The category with the largest number of genes was "transport" with 31 genes, followed by categories that include genes (n = 18) involved in carbohydrate metabolism or protein processes (i.e. biosynthesis, modification, transport and degradation). While up- or down-regulated genes were distributed homogenously within almost all the function groups, some categories included more up-regulated genes (ergosterol biosynthesis) or down-regulated genes (TCA cycle). As it will be discussed below, the finding of a large number of genes differentially regulated adds support to the concept that azole activity is beyond the inhibition of the lanosterol demethylase target encoded by ERG11 [32], whose overexpression has been associated with fungal resistance [33]. To further classify the genes regulated by FLC exposure, we performed GO term analysis. As expected, GO analysis of genes induced by FLC revealed a significant enrichment of genes involved in sterol metabolism, particularly ergosterol biosynthetic process (Table 2). Enrichment of genes repressed by FLC was observed in processes involving metabolism of amino acids and derivatives (Table 2).

Table 2 Gene Ontology (GO) term analysis for the C. neoformans FLC response
Full size table

Effect of FLC on genes involved in ergosterol biosynthesis and related pathways

Earlier efforts to profile the response of yeast cells (S. cerevisiae or C. albicans) to the short-term exposure to azole drugs implicated genes in the ergosterol biosynthetic pathway as major players [28, 29], thus indicating that this pathway is the target of azoles and is responsive to modulations in ergosterol levels. As shown in Table 1, we found that eight ERG genes (ERG1, ERG2, ERG3, ERG5, ERG7, ERG11, ERG13 and ERG25) exhibited increases in expression (2.09- to 3.95-fold) upon FLC treatment. This was a predictable result from the inhibition of Erg11 function by FLC, which is the rate-limiting step of the ergosterol biosynthetic pathway. Indeed, the idea of a compensatory response to re-establish the plasma membrane ergosterol levels [30] may account for the observed upregulation of either early (ERG13, ERG7 and ERG1) or late (ERG25, ERG2, ERG3 and ERG5) genes of the ergosterol pathway, in addition to upregulation of ERG11 itself (Table 1, ergosterol biosynthesis).

ERG13 encodes the enzyme hydroxymethylglutaryl-CoA synthase that catalyzes the production of hydroxymethylglutaryl-CoA from acetyl-CoA and acetoacetyl-CoA, and acts in the mevalonate biosynthesis, a precursor required for the biosynthesis of ergosterol. Acetyl-CoA is converted to carbon dioxide and water by enzymes (e.g. isocitrate dehydrogenase) that function in the TCA cycle, a central metabolic process in the mitochondria leading to produce, after oxidative phosphorylation, chemical energy in the form of ATP and NADH. Presumably, as a result of feedback control, we observed that several TCA cycle enzymes were downregulated in response to FLC (Table 1, TCA cycle), suggesting that C. neoformans may direct the cellular acetyl-CoA content to lipid (sterol) biosynthesis and metabolism to counterbalance ergosterol alteration.

Our particular interest was the up-regulation (4.04-fold) of SRE1, that belongs to a group of sterol regulatory element-binding proteins (SREBPs), first characterized in mammalian cells as regulator of lipid homeostasis [34]. While C. neoformans Sre1 regulates genes encoding ergosterol biosynthetic enzymes, SRE1 was shown to be required for growth and survival in the presence of azoles and also for virulence in a mouse model of cryptococcosis [18, 20, 35]. In addition, C. neoformans Sre1 stimulates ergosterol production in response to sterol depletion when the oxygen-dependent ergosterol synthesis is limited by hypoxia [36]. Consistently, C. neoformans mutants in the SREBP pathway showed reduction in ergosterol levels, increased sensitivity not only to low oxygen but also to several chemical agents, including azole antifungals, CoCl2 and reactive oxygen species (ROS)-generating compounds. Most importantly, these mutants showed reduced virulence in mice [37].

Effect of FLC on genes involved in cell structure and maintenance

Consequent to depletion of ergosterol and the concomitant accumulation of 14-methylated sterols, several plausible hypotheses on the mode of action of azoles were suggested by Vanden Bossche [32] two decades ago including alterations in membrane functions, synthesis and activity of membrane-bound enzymes, mitochondrial activities and uncoordinated activation of chitin synthesis. Transcript levels of several genes involving lipid and fatty acid metabolism decreased in the current study (Table 1), possibly in agreement with a remodelling of the cell membrane in response to reduced ergosterol levels. Conversely, expression of PLB1, that encodes Plb1, a known virulence factor in C. neoformans, was increased 2.18-fold. Phospholipases cleave fatty acid moieties from larger lipid molecules, releasing arachidonic acid for the production of eicosanoids that are utilized by the pathogenic yeasts C. neoformans and C. albicans to produce immunomodulatory prostaglandins [38]. In addition, cell wall-linked cryptococcal Plb1 contributes to cell wall integrity and is a source of secreted enzyme [39].

It was also expected that exposure to FLC would affect genes responsible for cell wall integrity. Two chitin synthase genes were found to be significantly up-regulated (2.20-fold for CHS2 and 3.62-fold for CHS7), concomitantly with down-regulated expression (4.35-fold) of the chitin deacetylase CDA3 (homolog to S. cerevisiae CDA2) (Table 1, cell wall maintenance). In C. albicans, activation of chitin synthesis, which is mediated by the PKC-, Ca2+/calcineurin-, and HOG- cell wall signalling pathways, appears to be an adaptive response to caspofungin treatment. Hence, subculturing caspofungin-resistant cells in the absence of caspofungin resulted in wild-type levels of chitin content [40]. While this form of drug tolerance is rationally accepted for a drug damaging the cell wall integrity (caspofungin is known to reduce β-glucan synthesis), it is also possible that exposure to azoles induces a salvage mechanism involving the up-regulation of chitin synthesis. Although known as a relatively minor cell wall component, chitin is thought to contribute significantly to cryptococcal wall strength and integrity [3]. Chitosan, the enzymatically deacetytaled form of chitin, helps to maintain cell integrity and is necessary for maintaining normal capsule width and retention of cell wall melanin [41]. Consistently, up-regulation was observed for BGL2 (2.61-fold) that encodes the glucantransferase (also termed glucosyltransferase) Bgl2, a major cell wall constituent described in a wide range of yeast species.

Effect of FLC on genes involved in cell stress and virulence

We found that FLC induced the expression of several genes involved in oxidative-stress response (Table 1, cell stress). One of these genes, GRE2, was induced 3.54-fold, consistent with the previous observation that transcripts from GRE2 and other stress-induced genes (YDR453C and SOD2) were increased in S. cerevisiae exposed to azoles [28]. Interestingly, loss of Gre2 is impairing tolerance to ergosterol biosynthesis disrupting agents (i.e. clotrimazole and ketoconazole), further supporting an association between GRE2 and ergosterol metabolism [42]. YHB1 that encodes a flavo-haemoglobin able to detoxify nitric oxide in C. albicans and C. neoformans was down-regulated 2.32-fold in our study, which is opposed to its established relevance in vivo [43]. A strong reduction in the expression of FHB1 (the C. neoformans ortholog of YHB1) was also observed during growth of C. neoformans at 37°C compared to 25°C, indicating that regulation of this gene or its product at the posttranslational level may occur in response to environmental changes [44]. In contrast, CTA1 encoding catalase in S. cerevisiae was induced (2.81-fold) by FLC exposure. Together with TSA3 (2.09-fold) encoding thiol-specific antioxidant protein 3 (Table 1, cell stress) and other responsive genes with oxidoreductase activity (Table 1, oxidoreduction), these genes may function in response to oxidative stress. Accordingly, the stress-related gene encoding Ssa1 was also up-regulated (2.48-fold). This C. neoformans protein (Hsp70 family member) acts in vivo as transcriptional co-activator of laccase [45] and is important for the production of melanin, which is a free-radical scavenger playing a protective role in stress resistance [17].

The C. neoformans polysaccharide capsule is a complex structure that is required for virulence [46, 47]. Interestingly, the capsule-associated gene CAS3 [48] was found to be up-regulated (12.16-fold) upon exposure to the drug (Table 1, capsule synthesis). This gene encodes a protein belonging to a seven-member protein family that includes Cap64. Treatment with FLC did not significantly change expression of the essential capsule-producing genes, CAP10, CAP59, CAP60 and CAP64. Since the cryptococcal cell wall is needed for the localization or attachment of known or putative virulence factors other than capsule (i.e. melanin, Plb1 and Bgl2), it could be hypothesized that FLC induces alterations in the cell wall which in turns affects the expression of these factors. An alternative hypothesis would be that FLC acts as a stress-generating molecule and triggers enhanced expression of virulence determinant(s) that enable to survive in hostile environments.

Effect of FLC on genes involved in cellular transport

Several genes involved in small molecule transport and vesicular transport were either up- or down-regulated in response to FLC (Table 1, transport). These include DUR3 (plasma membrane transporter for urea, up-regulated by 4.78-fold), MEP2/AMT2 (ammonium permease, up-regulated by 3.78-fold) and AQY1 (aquaporin water channel, up-regulated by 2.73-fold), which all belong to the group of C. neoformans genes regulated by osmotic stress [49]. It is possible that defects in the plasma membrane resulting from inhibition of ergosterol biosynthesis by FLC affects transport of small molecules through the membrane. Analysis of the H99 genome sequence [16] predicted 54 ATP-Binding Cassette (ABC) transporters and 159 major facilitator superfamily (MFS) transporters, suggesting wide transport capabilities of this environmental yeast [50]. However, we found only two S. cerevisiae transporter homologues with significant increased expression. One is PDR15 that is a member of the ABC transporter subfamily exporting antifungals and other xenobiotics in fungi [51]. The other gene is ATR1 that encodes a multidrug resistance transport protein belonging to the MFS class of transporters. ATR1 expression was recently shown to be upregulated by boron and several stress conditions [52]. To date, Afr1 (encoded by AFR1; also termed CneAfr1) and CneMdr1 are the only two efflux pumps associated with antifungal drug resistance in C. neoformans [50]. Since Afr1 is the major efflux pump mediating azole resistance in C. neoformans [11, 15], the absence of altered AFR1 expression could be expected. Not surprisingly, we noticed downregulated expression (2.35-fold) of FLR1 (for fluconazole resistance) encoding a known MFS multidrug transporter in yeast, that is able to confer resistance to a wide range of dissimilar drugs and other chemicals [53]. This may suggest that both AFR1 and FLR1 do not participate to the short-term stress induced by FLC in C. neoformans.

Effect of FLC on the susceptibility to cell wall inhibitors

It was demonstrated that compounds interfering with normal cell wall formation (Congo red, calcofluor white, SDS and caffeine) affect growth of C. neoformans strains with altered cell wall integrity [27]. For instance, several deletion strains for genes involved in the PKC1 signal transduction pathway were found to be sensitive to SDS and Congo red and to a lesser extent caffeine. To test the hypothesis that FLC treatment might induce cell wall stress, we analyzed H99 cells for susceptibility to the cell wall perturbing agents, before and after the cells were exposed for 90 min to FLC at sub-MIC concentration (10 mg/l) at 30°C. Phenotypes of H99 cells on cell wall inhibitor plates are shown in Figure 3. The FLC pre-treated H99 cells were slightly more resistant to all four cell wall inhibitors as compared to untreated cells. These findings are consistent with expression changes of cell wall associated genes identified in our microarray analysis. Particularly, since calcofluor white (which binds to chitin) disrupts the cell wall and Congo red (which binds to β-glucans) interferes with the cell wall biogenesis [27], the altered regulation of genes involved in the chitin (CHS2 and CHS7) and glucan (BGL2) synthesis may explain the phenotype of decreased susceptibility to cell wall stress exhibited by FLC-exposed cells. Similar results were obtained when H99 cells were pre-treated with FLC at 37°C (see Additional file 2).

Figure 3
figure 3

Cell wall integrity assays with H99 C. neoformans cells left untreated (H99) or exposed to FLC (H99F) at a sub-MIC concentration of 10 mg/l for 90 min at 30°C. Cells were grown at the same temperature for 48 h on YEPD supplemented with calcofluor white (CFW), Congo red, sodium dodecyl sulphate (SDS) and caffeine. Aliquots of cells were applied onto the agar surface with 10-fold serial dilutions.

Full size image

Effect of FLC on the susceptibility to H2O2

Because a number of FLC-responsive transcriptional changes was found to affect genes involved in the oxidative stress response (i.e. CTA1, GRE2), it seemed reasonable to examine whether FLC at sub-inhibitory concentrations could induce oxidative stress resistance in vitro. For this purpose, exponentially growing H99 cells that were treated with 10 mg/l FLC for 90 min were subjected to an additional challenge with 20 mM H2O2. The viable cells were next quantified on YEPD plates after 0.5, 1, 1.5 and 2 h of additional growth. As shown in Figure 4, while untreated cells showed a high degree of cell death, cells treated with FLC exhibited gained more viability upon oxidative exposure at the endpoints of 1, 1.5 and 2 h. Similar results were obtained when H99 cells were pre-treated with FLC at 37°C (see Additional file 3). These findings indicate that FLC exposure is able to generate protection against oxidative stress in vitro, possibly as a result of a transcriptional adaptive response.

Figure 4
figure 4

Survival of C. neoformans after oxidative treatment. Exponentially growing cells were left untreated (H99) or exposed to 10 mg/l FLC (H99F) for 90 min at 30°C and then challenged with 20 mM H2O2 for 2 h. Aliquots were harvested at given time points and cell viability performed as described in Methods. Plotted values are means of three experiments

Full size image

Conclusions

Although exposure to azoles has been already investigated in several other fungal species and the transcriptional profile of differentially expressed genes was obtained using a single FLC concentration and time point, our study reveals several interesting findings. First, we demonstrated that short-term exposure of C. neoformans to FLC resulted in a complex altered gene expression profile. These genes included not only genes commonly responding to diverse environmental stresses, such as oxidative and drug stresses, but also genes encoding virulence factors (i.e. Plb1, Sre1 and capsule). Second, we corroborated the potential of genome-wide transcriptional analyses to envisage alternative therapeutic strategies for cryptococcosis. Apart from ergosterol and its biosynthesis, there are yet few other targets to be exploited in anticryptococcal therapy. Therefore, elucidation of molecular processes underlying the physiological responses of cryptococcal cells to FLC could serve not only to identify novel treatment approaches but also to potentiate the inhibitory effects of existing azole drugs. Our findings show that the phenomena described can apply to the in vivo situation, i.e. during azole maintenance therapy in the host, but transcriptional analyses using different growth conditions of H99 cells, mimicking stress conditions encountered during a human meningeal infection, may reveal new fields to pursue for anticryptococcal therapy.

References

  1. Perfect JR, Casadevall A: Cryptococcosis. Infect Dis Clin North Am. 2002, 16: 837-874. 10.1016/S0891-5520(02)00036-3.

    PubMed  Article  Google Scholar 

  2. Bicanic T, Harrison TS: Cryptococcal meningitis. Br Med Bull. 2005, 72: 99-118. 10.1093/bmb/ldh043.

    PubMed  Article  Google Scholar 

  3. Doering TL: How sweet it is! Cell wall biogenesis and polysaccharide capsule formation in Cryptococcus neoformans. Annu Rev Microbiol. 2009, 63: 223-247. 10.1146/annurev.micro.62.081307.162753.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  4. Silveira FP, Husain S: Fungal infections in solid organ transplantation. Med Mycol. 2007, 45: 305-320. 10.1080/13693780701200372.

    PubMed  Article  Google Scholar 

  5. Thakur R, Sarma S, Kushwaha S: Prevalence of HIV-associated cryptococcal meningitis and utility of microbiological determinants for its diagnosis in a tertiary care center. Indian J Pathol Microbiol. 2008, 51: 212-214. 10.4103/0377-4929.41689.

    PubMed  Article  Google Scholar 

  6. Park BJ, Wannemuehler KA, Marston BJ, Govender N, Pappas PG, Chiller TM: Estimation of the current global burden of cryptococcal meningitis among persons living with HIV/AIDS. AIDS. 2009, 23: 525-530. 10.1097/QAD.0b013e328322ffac.

    PubMed  Article  Google Scholar 

  7. Datta K, Jain N, Sethi S, Rattan A, Casadevall A, Banerjee U: Fluconazole and itraconazole susceptibilities of clinical isolates of Cryptococcus neoformans at a tertiary care centre in India; a need for care. J Antimicrob Chemother. 2003, 52: 683-686. 10.1093/jac/dkg399.

    PubMed  CAS  Article  Google Scholar 

  8. Sar B, Monchy D, Vann M, Keo C, Sarthou JL, Buisson Y: Increasing in vitro resistance to fluconazole in Cryptococcus neoformans Cambodian isolates: April 2000 to March 2002. J Antimicrob Chemother. 2004, 54: 563-565. 10.1093/jac/dkh361.

    PubMed  CAS  Article  Google Scholar 

  9. Bii CC, Makimura K, Abe S, Taguchi H, Mugasia OM, Revathi G, Wamae NC, Kamiya S: Antifungal drug susceptibility of Cryptococcus neoformans from clinical sources in Nairobi, Kenya. Mycoses. 2007, 50: 25-30. 10.1111/j.1439-0507.2006.01293.x.

    PubMed  CAS  Article  Google Scholar 

  10. Perfect JR, Cox GM: Drug resistance in Cryptococcus neoformans. Drug Resist Updat. 1999, 2: 259-269. 10.1054/drup.1999.0090.

    PubMed  CAS  Article  Google Scholar 

  11. Sionov E, Chang YC, Garraffo HM, Kwon-Chung KJ: Heteroresistance to fluconazole in Cryptococcus neoformans is intrinsic and associated with virulence. Antimicrob Agents Chemother. 2009, 53: 2804-2815. 10.1128/AAC.00295-09.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  12. Mondon P, Petter R, Amalfitano G, Luzzati R, Concia E, Polacheck I, Kwon-Chung KJ: Heteroresistance to fluconazole and voriconazole in Cryptococcus neoformans. Antimicrob Agents Chemother. 1999, 43: 1856-1861.

    PubMed  CAS  PubMed Central  Google Scholar 

  13. Sionov E, Lee H, Chang YC, Kwon-Chung KJ: Cryptococcus neoformans overcomes stress of azole drugs by formation of disomy in specific multiple chromosomes. PLoS Pathog. 2010, 6: e1000848-10.1371/journal.ppat.1000848.

    PubMed  PubMed Central  Article  Google Scholar 

  14. Lupetti A, Danesi R, Campa M, Del Tacca M, Kelly S: Molecular basis of resistance to azole antifungals. Trends Mol Med. 2002, 8: 76-81. 10.1016/S1471-4914(02)02280-3.

    PubMed  CAS  Article  Google Scholar 

  15. Posteraro B, Sanguinetti M, Sanglard D, La Sorda M, Boccia S, Romano L, Morace G, Fadda G: Identification and characterization of a Cryptococcus neoformans ATP binding cassette (ABC) transporter-encoding gene, CnAFR1, involved in the resistance to fluconazole. Mol Microbiol. 2003, 47: 357-371. 10.1046/j.1365-2958.2003.03281.x.

    PubMed  CAS  Article  Google Scholar 

  16. Loftus BJ, Fung E, Roncaglia P, Rowley D, Amedeo P, Bruno D, Vamathevan J, Miranda M, Anderson IJ, Fraser JA, Allen JE, Bosdet IE, Brent MR, Chiu R, Doering TL, Donlin MJ, D'Souza CA, Fox DS, Grinberg V, Fu J, Fukushima M, Haas BJ, Huang JC, Janbon G, Jones SJ, Koo HL, Krzywinski MI, Kwon-Chung JK, Lengeler KB, Maiti R, Marra MA, Marra RE, Mathewson CA, Mitchell TG, Pertea M, Riggs FR, Salzberg SL, Schein JE, Shvartsbeyn A, Shin H, Shumway M, Specht CA, Suh BB, Tenney A, Utterback TR, Wickes BL, Wortman JR, Wye NH, Kronstad JW, Lodge JK, Heitman J, Davis RW, Fraser CM, Hyman RW: The genome of the basidiomycetous yeast and human pathogen Cryptococcus neoformans. Science. 2005, 307: 1321-1324. 10.1126/science.1103773.

    PubMed  PubMed Central  Article  Google Scholar 

  17. Brown SM, Campbell LT, Lodge JK: Cryptococcus neoformans, a fungus under stress. Curr Opin Microbiol. 2007, 10: 320-325. 10.1016/j.mib.2007.05.014.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  18. Chun CD, Liu OW, Madhani HD: A link between virulence and homeostatic responses to hypoxia during infection by the human fungal pathogen Cryptococcus neoformans. PLoS Pathog. 2007, 3: e22-10.1371/journal.ppat.0030022.

    PubMed  PubMed Central  Article  Google Scholar 

  19. Coleman JJ, Mylonakis E: Efflux in fungi: la pièce de résistance. PLoS Pathog. 2009, 5: e1000486-10.1371/journal.ppat.1000486.

    PubMed  PubMed Central  Article  Google Scholar 

  20. Chang YC, Bien CM, Lee H, Espenshade PJ, Kwon-Chung KJ: Sre1p, a regulator of oxygen sensing and sterol homeostasis, is required for virulence in Cryptococcus neoformans. Mol Microbiol. 2007, 64: 614-629. 10.1111/j.1365-2958.2007.05676.x.

    PubMed  CAS  Article  Google Scholar 

  21. Eisenman HC, Casadevall A, McClelland EE: New insights on the pathogenesis of invasive Cryptococcus neoformans infection. Curr Infect Dis Rep. 2007, 9: 457-464. 10.1007/s11908-007-0070-8.

    PubMed  Article  Google Scholar 

  22. Soteropoulos P, Vaz T, Santangelo R, Paderu P, Huang DY, Tamás MJ, Perlin DS: Molecular characterization of the plasma membrane H+-ATPase, an antifungal target in Cryptococcus neoformans. Antimicrob Agents Chemother. 2000, 44: 2349-2355. 10.1128/AAC.44.9.2349-2355.2000.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  23. Sanguinetti M, Posteraro B, La Sorda M, Torelli R, Fiori B, Santangelo R, Delogu G, Fadda G: Role of AFR1, an ABC transporter-encoding gene, in the in vivo response to fluconazole and virulence of Cryptococcus neoformans. Infect Immun. 2006, 74: 1352-1359. 10.1128/IAI.74.2.1352-1359.2006.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  24. Kim MS, Ko YJ, Maeng S, Floyd A, Heitman J, Bahn YS: Comparative transcriptome analysis of the CO2 sensing pathway via differential expression of carbonic anhydrase in Cryptococcus neoformans. Genetics. 2010, 185: 1207-1219. 10.1534/genetics.110.118315.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  25. Barrett ER: Gene Expression Omnibus (GEO): Microarray data storage, submission, retrieval, and analysis. Methods Enzymol. 2006, 411: 352-369.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  26. Arana DM, Nombela C, Pla J: Fluconazole at subinhibitory concentrations induces the oxidative- and nitrosative-responsive genes TRR1, GRE2 and YHB1, and enhances the resistance of Candida albicans to phagocytes. J Antimicrob Chemother. 2010, 65: 54-62. 10.1093/jac/dkp407.

    PubMed  CAS  Article  Google Scholar 

  27. Gerik KJ, Donlin MJ, Soto CE, Banks AM, Banks IR, Maligie MA, Selitrennikoff CP, Lodge JK: Cell wall integrity is dependent on the PKC1 signal transduction pathway in Cryptococcus neoformans. Mol Microbiol. 2005, 58: 393-408. 10.1111/j.1365-2958.2005.04843.x.

    PubMed  CAS  Article  Google Scholar 

  28. Bammert GF, Fostel JM: Genome-wide expression patterns in Saccharomyces cerevisiae: comparison of drug treatments and genetic alterations affecting biosynthesis of ergosterol. Antimicrob Agents Chemother. 2000, 44: 1255-1265. 10.1128/AAC.44.5.1255-1265.2000.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  29. De Backer MD, Ilyina T, Ma XJ, Vandoninck S, Luyten WH, Vanden Bossche H: Genomic profiling of the response of Candida albicans to itraconazole treatment using a DNA microarray. Antimicrob Agents Chemother. 2001, 45: 1660-1670. 10.1128/AAC.45.6.1660-1670.2001.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  30. Gamarra S, Rocha EM, Zhang YQ, Park S, Rao R, Perlin DS: Mechanism of the synergistic effect of amiodarone and fluconazole in Candida albicans. Antimicrob Agents Chemother. 2010, 54: 1753-1761. 10.1128/AAC.01728-09.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  31. Karababa M, Coste AT, Rognon B, Bille J, Sanglard D: Comparison of gene expression profiles of Candida albicans azole-resistant clinical isolates and laboratory strains exposed to drugs inducing multidrug transporters. Antimicrob Agents Chemother. 2004, 48: 3064-3079. 10.1128/AAC.48.8.3064-3079.2004.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  32. Vanden Bossche H: Biochemical targets for antifungal azole derivatives: hypothesis on the mode of action. Curr Top Med Mycol. 1985, 1: 313-351.

    PubMed  CAS  Article  Google Scholar 

  33. Henry KW, Nickels JT, Edlind TD: Upregulation of ERG genes in Candida species by azoles and other sterol biosynthesis inhibitors. Antimicrob Agents Chemother. 2000, 44: 2693-2700. 10.1128/AAC.44.10.2693-2700.2000.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  34. Goldstein JL, Debose-Boyd RA, Brown MS: Protein sensors for membrane sterols. Cell. 2006, 124: 35-46. 10.1016/j.cell.2005.12.022.

    PubMed  CAS  Article  Google Scholar 

  35. Bien CM, Chang YC, Nes WD, Kwon-Chung KJ, Espenshade PJ: Cryptococcus neoformans Site-2 protease is required for virulence and survival in the presence of azole drugs. Mol Microbiol. 2009, 74: 672-690. 10.1111/j.1365-2958.2009.06895.x.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  36. Bien CM, Espenshade PJ: Sterol regulatory element binding proteins in fungi: hypoxic transcription factors linked to pathogenesis. Eukaryot Cell. 2010, 9: 352-359. 10.1128/EC.00358-09.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  37. Chang YC, Ingavale SS, Bien C, Espenshade P, Kwon-Chung KJ: Conservation of the sterol regulatory element-binding protein pathway and its pathobiological importance in Cryptococcus neoformans. Eukaryot Cell. 2009, 8: 1770-1779. 10.1128/EC.00207-09.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  38. Rhome R, Del Poeta M: Lipid signaling in pathogenic fungi. Annu Rev Microbiol. 2009, 63: 119-131. 10.1146/annurev.micro.091208.073431.

    PubMed  CAS  Article  Google Scholar 

  39. Siafakas AR, Sorrell TC, Wright LC, Wilson C, Larsen M, Boadle R, Williamson PR, Djordjevic JT: Cell wall-linked cryptococcal phospholipase B1 is a source of secreted enzyme and a determinant of cell wall integrity. J Biol Chem. 2007, 282: 37508-37514. 10.1074/jbc.M707913200.

    PubMed  CAS  Article  Google Scholar 

  40. Walker LA, Gow NA, Munro CA: Fungal echinocandin resistance. Fungal Genet Biol. 2010, 47: 117-126. 10.1016/j.fgb.2009.09.003.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  41. Baker LG, Specht CA, Donlin MJ, Lodge JK: Chitosan, the deacetylated form of chitin, is necessary for cell wall integrity in Cryptococcus neoformans. Eukaryot Cell. 2007, 6: 855-867. 10.1128/EC.00399-06.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  42. Warringer J, Blomberg A: Involvement of yeast YOL151W/GRE2 in ergosterol metabolism. Yeast. 2006, 23: 389-398. 10.1002/yea.1363.

    PubMed  CAS  Article  Google Scholar 

  43. de Jesús-Berríos M, Liu L, Nussbaum JC, Cox GM, Stamler JS, Heitman J: Enzymes that counteract nitrosative stress promote fungal virulence. Curr Biol. 2003, 13: 1963-1968. 10.1016/j.cub.2003.10.029.

    PubMed  Article  Google Scholar 

  44. Kraus PR, Boily MJ, Giles SS, Stajich JE, Allen A, Cox GM, Dietrich FS, Perfect JR, Heitman J: Identification of Cryptococcus neoformans temperature-regulated genes with a genomic-DNA microarray. Eukaryot Cell. 2004, 3: 1249-1260. 10.1128/EC.3.5.1249-1260.2004.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  45. Zhang S, Hacham M, Panepinto J, Hu G, Shin S, Zhu X, Williamson PR: The Hsp70 member, Ssa1, acts as a DNA-binding transcriptional co-activator of laccase in Cryptococcus neoformans. Mol Microbiol. 2006, 62: 1090-1101. 10.1111/j.1365-2958.2006.05422.x.

    PubMed  CAS  Article  Google Scholar 

  46. Fromtling RA, Shadomy HJ, Jacobson ES: Decreased virulence in stable, acapsular mutants of Cryptococcus neoformans. Mycopathologia. 1982, 79: 23-29. 10.1007/BF00636177.

    PubMed  CAS  Article  Google Scholar 

  47. Chang YC, Kwon-Chung KJ: Complementation of a capsule-deficient mutation of Cryptococcus neoformans restores its virulence. Mol Cell Biol. 1994, 14: 4912-4919.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  48. Moyrand F, Chang YC, Himmelreich U, Kwon-Chung KJ, Janbon G: Cas3p belongs to a seven-member family of capsule structure designer proteins. Eukaryot Cell. 2004, 3: 1513-1524. 10.1128/EC.3.6.1513-1524.2004.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  49. Ko YJ, Yu YM, Kim GB, Lee GW, Maeng PJ, Kim S, Floyd A, Heitman J, Bahn YS: Remodeling of global transcription patterns of Cryptococcus neoformans genes mediated by the stress-activated HOG signaling pathways. Eukaryot Cell. 2009, 8: 1197-1217. 10.1128/EC.00120-09.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  50. Cannon RD, Lamping E, Holmes AR, Niimi K, Baret PV, Keniya MV, Tanabe K, Niimi M, Goffeau A, Monk BC: Efflux-mediated antifungal drug resistance. Clin Microbiol Rev. 2009, 22: 291-321. 10.1128/CMR.00051-08.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  51. Seret ML, Diffels JF, Goffeau A, Baret PV: Combined phylogeny and neighborhood analysis of the evolution of the ABC transporters conferring multiple drug resistance in hemiascomycete yeasts. BMC Genomics. 2009, 10: 459-10.1186/1471-2164-10-459.

    PubMed  PubMed Central  Article  Google Scholar 

  52. Kaya A, Karakaya HC, Fomenko DE, Gladyshev VN, Koc A: Identification of a novel system for boron transport: Atr1 is a main boron exporter in yeast. Mol Cell Biol. 2009, 29: 3665-3674. 10.1128/MCB.01646-08.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  53. Sá-Correia I, dos Santos SC, Teixeira MC, Cabrito TR, Mira NP: Drug:H+ antiporters in chemical stress response in yeast. Trends Microbiol. 2009, 17: 22-31. 10.1016/j.tim.2008.09.007.

    PubMed  Article  Google Scholar 

Download references

Acknowledgements

This work was supported by grants from the Istituto di Ricovero e Cura a Carattere Scientifico (IRCCS) Lazzaro Spallanzani (Strategic Research Program 2006) to GF, from the Università Cattolica del S. Cuore (Fondi Ateneo Linea D1-2009) to MS, and from the Swiss Research National Foundation 31003A_127378 to DS.

Author information

Affiliations

  1. Istituto di Microbiologia, Università Cattolica del Sacro Cuore, Largo F. Vito 1, 00168, Rome, Italy

    Ada Rita Florio, Elena De Carolis, Riccardo Torelli, Giovanni Fadda, Maurizio Sanguinetti & Brunella Posteraro

  2. Institute of Microbiology, University of Lausanne and University Hospital Center, Lausanne, Rue du Bugnon 48, CH-1011, Lausanne, Switzerland

    Selene Ferrari & Dominique Sanglard

Authors
  1. Ada Rita Florio
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Selene Ferrari
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Elena De Carolis
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Riccardo Torelli
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Giovanni Fadda
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Maurizio Sanguinetti
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Dominique Sanglard
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Brunella Posteraro
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Maurizio Sanguinetti.

Additional information

Authors' contributions

MS, DS and BP designed the study; ARF and SF carried out the experimental work; ARF, EDC and RT analysed the data; ARF and BP wrote the manuscript. GF and DS corrected the manuscript. All the authors read and approved the final manuscript.

Electronic supplementary material

12866_2010_1395_MOESM1_ESM.DOC

Additional file 1: Table A1 Primers and fluorescent probes used in qRT-PCR. Contains Table A1 showing the qRT-PCR primers and probes. (DOC 58 KB)

12866_2010_1395_MOESM2_ESM.DOC

Additional file 2: Figure A1 Cell wall integrity assays with H99 C. neoformanscells left untreated (H99) or exposed to FLC (H99F) at a sub-MIC concentration of 10 mg/l for 90 min at 37°C. Cells were grown at the same temperature for 48 h on YEPD supplemented with calcofluor white (CFW), Congo red, sodium dodecyl sulphate (SDS) and caffeine. Aliquots of cells were applied onto the agar surface with 10-fold serial dilutions. Contains Figure A1 showing the results of cell wall inhibitors susceptibility assays for H99 cells pre-treated with FLC at 37°C. (DOC 122 KB)

12866_2010_1395_MOESM3_ESM.DOC

Additional file 3: Figure A2 Survival of C. neoformans after oxidative treatment. Exponentially growing cells were left untreated (H99) or exposed to 10 mg/l FLC (H99F) for 90 min at 37°C and then challenged with 20 mM H2O2 for 2 h. Aliquots were harvested at given time points and cell viability performed as described in Methods. Plotted values are means of three experiments. Contains Figure A2 showing the results of H2O2 susceptibility assays for H99 cells pre-treated with FLC at 37°C. (DOC 262 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Florio, A.R., Ferrari, S., De Carolis, E. et al. Genome-wide expression profiling of the response to short-term exposure to fluconazole in Cryptococcus neoformans serotype A. BMC Microbiol 11, 97 (2011). https://doi.org/10.1186/1471-2180-11-97

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2180-11-97

Keywords

  • Azole
  • Ergosterol
  • Caspofungin
  • Cryptococcosis
  • Cryptococcal Meningitis

BMC Microbiology

ISSN: 1471-2180

Contact us