Original Research ARTICLE
Interactions between copy number and expression level of genes involved in fluconazole resistance in Candida glabrata
- 1Laboratoire de Biologie Moléculaire Parasitaire et Fongique, Faculté de médecine, University of Sfax, Sfax, Tunisie
- 2Laboratoire de Parasitologie-Mycologie, Aix-Marseille Université, AP-HM, Marseille, France
Objectives: This study aimed to elucidate the relative involvement of drug resistance gene copy number and overexpression in fluconazole resistance in clinical C. glabrata isolates using a population-based approach.
Methods: Fluconazole resistance levels were quantified using the minimal inhibitory concentration (MIC) via Etest method. Both gene expression levels and gene copy number of CgCDR1, CgPDH1, CgERG11, and CgSNQ2 were assessed via quantitative real-time PCR. The influence of the main effects and first-level interactions of both the expression level and copy number of these genes on fluconazole resistance levels were analyzed using a multivariate statistical model.
Results: Forty-three C. glabrata isolates were collected from 30 patients during in a hospital survey. In the multivariate analysis, C. glabrata fluconazole MICs were independently increased by CgSNQ2 overexpression (p < 10−4) and the interaction between CgPDH1 gene copy number and CgPDH1 expression level (p = 0.038). In contrast, both CgPDH1 overexpression (p = 0.049) and the interaction between CgSNQ2 and CgERG11 expression (p = 0.003) led to a significant decrease in fluconazole MICs.
Conclusion: Fluconazole resistance in C. glabrata involves complex interactions between drug resistance gene expression and/or copy number. The population-based multivariate analysis highlighted the involvement of the CgSNQ2 gene in fluconazole resistance and the complex effect of the other genes such as PDH1 for which overexpression was associated with reduced fluconazole resistance levels, while the interaction between PDH1 overexpression and copy number was associated with increased resistance levels.
Candida glabrata is responsible for the increase in cases of systemic and superficial candidiasis in many countries (Diekema et al., 2002; Li et al., 2007; Pfaller and Diekema, 2007). Several reports have revealed that a significant percentage of C. glabrata clinical isolates are resistant to fluconazole. Pfaller et al. have found that between 2001 and 2007, isolation frequencies and fluconazole resistance among C. glabrata blood stream infection isolates was increased compared with the 1992–2001 period (Pfaller and Diekema, 2004; Pfaller et al., 2009). Similarly, antifungal susceptibility surveys have detected 17 and 18% resistance levels in Europe and North America, respectively (Pfaller et al., 2009). Many studies have shown that drug efflux mediated by the overexpression of ATP-binding cassette (ABC) transporter family membrane proteins was the main molecular mechanism associated with C. glabrata fluconazole resistance (Sanglard et al., 1999, 2001; Bouchara et al., 2000; Niimi et al., 2002; Wada et al., 2002; Cernicka and Subik, 2006; Coleman and Mylonakis, 2009). Studies have also shown that among ABC transporters, CgCDR1, and CgPDH1, but not CgSNQ2, play a major role in fluconazole resistance (Sanguinetti et al., 2005; Torelli et al., 2008). Azole resistance in C. albicans was primarily associated with mutation or overexpression of ERG11, which encodes the key azole target protein lanosterol 14-alpha demethylase (Chau et al., 2004; Cernicka and Subik, 2006). In contrast, CgERG11 mutations were marginally involved in C. glabrata fluconazole resistance (Sanguinetti et al., 2005). However, CgERG11 overexpression has been associated with chromosome duplication and aneuploidy due to genome plasticity as well as fluconazole resistance in Cryptococcus neoformans, C. albicans, and C. glabrata (Marichal et al., 1997; Kwon-Chung and Chang, 2012). This study aimed to dissect the relative contribution of drug resistance gene copy number and overexpression on fluconazole resistance in clinical C. glabrata isolates.
Patients and Isolates
The population of Candida glabrata clinical isolates included in the present study were recovered during an epidemiological survey of antifungal resistance conducted at the Habib Bourguiba University Hospital in Sfax from January 2005 to December 2007 as previously described (Abbes et al., 2011, 2012; Amouri et al., 2011; Sellami et al., 2011). The specimens were isolated from 30 patients with urinary tract infection, vaginal infection or systemic infections. Candida glabrata was identified by assessing the isolates using Candiselect medium (bioMérieux, Craponne, France) and the ID32C profile (bioMérieux). Isolate fluconazole minimum inhibitory concentrations (MICs) were then measured using the E-test method (AB Biodisk, Sweden) on RPMI 1640 medium (AES, France) as recommended by the manufacturer. Isolates with MICs >64 mg/l were defined as fluconazole resistant (Clinical and Laboratory Standards Institute, 2008).
RNA and DNA Isolation
Total RNA was extracted from yeast cells isolated from RPMI plates during the logarithmic growth phase at 48 h of incubation using an RNAeasy Protect™ mini kit (Qiagen, Courtaboeuf, France) according to the manufacturer's instructions. RNA extracts were treated with RNase-free DNase (Qiagen) to avoid DNA contamination. Genomic DNA was extracted using the NucleoSpin™ Tissue method (Macherey-Nagel, Dueren, Germany) according to the manufacturer's instructions. Genomic DNA (gDNA) extracts were treated with DNase-free RNase (Qiagen) to eliminate RNA contamination.
Quantitative Real-Time PCR
Quantitative real-time PCR analysis was performed to measure the expression levels and copy number of each of the four target genes: CgCDR1, CgPDH1, CgSNQ2, and CgERG11 (Table 1). Gene transcript levels and copy number values were normalized to URA3, a single copy gene encoding orotidine 5-phosphate decarboxylase.
Table 1. Primers and probes used to amplify CgCDR1, CgPDH1 (formerly denoted CDR2), CgSNQ2, CgERG11, and CgURA3 in separate quantitative Real-Time PCR reactions as described in Wada et al. (2002).
Reverse transcriptase PCR (RT-PCR)
First, the RNA extract of each isolate was transcribed into cDNA using reverse transcriptase (RT), which was performed with a 25-μl reaction volume containing 2.5 μl MuLva buffer (Biology, France), 2 U of MuLva RT (Biology), 5 μl RNA sample and 0.5 μM random primer and polyA mixture (Eurogentec, France). The samples were subjected to reverse transcription at 37°C for 60 min.
Quantitative real-time PCR (qPCR)
Two μl of either cDNA or gDNA template was used to separately assess expression levels and gene copy number of the four genes: CgCDR1, CgPDH1, CgSNQ2, and CgERG11, as previously described (Sanguinetti et al., 2005). The primers utilized are detailed in Table 1. Briefly, 200 pM of each primer, 200 pM hydrolysis probe and 25 μl TaqMan mix (Applied Biosystems, Courtaboeuf, France) were included in the PCR assay. Copy number PCR cycling conditions consisted of an initial step at 94°C for 10 min and then 40 cycles of 15 s at 94°C, 30 s at 55°C and 20 s at 72°C. Fluorescence data were collected during the extension step and analyzed with a MX 4000 real-time PCR system (Stratagene, France). Each reaction was performed in duplicate. Copy number efficiency was determined for each gene. The normalized relative quantities of the gene-specific product or copy number for each sample were then calculated, i.e., the quantity of the samples of interest was compared with a calibrator (TU10, a fluconazole susceptible control, MIC = 0.125 mg/l). The quantity of both the calibrator and the samples of interest were normalized to URA3.
C. glabrata Genotyping
The C. glabrata isolates were genotyped using six markers, five microsatellites (RPM2, MTI, GLM4, GLM4, and GLM6) and irregular patterns in the ERG3 gene as previously described (Abbes et al., 2011, 2012).
Normalized gene copy numbers and expression levels are indicated using range, median, and interquartile range (IQR). The correlation between fluconazole MICs and the relative expression level and copy number of each gene was analyzed using the Spearman rank correlation coefficient. For qualitative analysis, we selected a threshold value of 2.5 normalized gene expression or normalized gene copy number to define gene overexpression or an increase in gene copy number, respectively. Regarding quantitative analysis, fluconazole MIC values were log-transformed, standardized and centered. The effect on fluconazole MICs due to the gene copy number and expression levels of each of the four genes studied was modeled using the Genmod procedure. All first-level interactions (the product of the two main effect terms) between each covariate were tested. When an interaction was statistically significant, each of the two main effects was retained in the model, regardless of the level of statistical significance. The most parsimonious best-fitting model was selected using the log-likelihood ratio test. After selection of the model, the generalized estimating equation option was used to account for the non-independence of isolates sampled from a same patient. Statistical analyses were performed with SAS 9.2 statistical software (Cary, NC, USA). All statistical tests were two-sided, and p < 0.05 was considered significant.
The characteristics of the 43 C. glabrata clinical isolates collected during a hospital survey from 30 patients (9 patients had ≥2 episodes) at the Habib Bourguiba University hospital in Sfax (Tunisia) are detailed in Table 2. Briefly, most isolates were grown from urine samples, and 15 isolates were involved in invasive infections. Fluconazole MICs ranged from 0.016 to 256 mg/l; MIC50 and MIC90 were 4 mg/l and 256 mg/l, respectively. Nine isolates were classified as resistant to fluconazole with a MIC >256 mg/l. Microsatellite analysis displayed 17 distinct genotypes. Twenty-two isolates were sequentially collected (7–220 days lag in the same patient). In patients 1, 2, 3, and 4 treated with fluconazole, MICs were significantly increased in sequential C. glabrata isolates, whereas the microsatellite-based genotypes remained stable (Tables 3, 4). These patients were exposed to 150, 30, 220, or 120 days of fluconazole therapy with a cumulative dose of 4, 5.6, 8, and 12 g fluconazole, respectively (Table 3). In patients 6 and 11, the fluconazole MICs were increased in sequential isolates with a distinct microsatellite genotype, thereby indicating that they probably acquired resistant C. glabrata strains upon exposure to fluconazole (Table 2). Finally, both microsatellite genotypes and fluconazole MICs remained stable in the sequential isolates derived from patients 19, 21, and 26 (Table 2).
Table 2. Relative gene expression levels and copy numbers of CgCDR1, CgPDH1, CgERG11, and CgSNQ2 in fluconazole sensitive and resistant Candida glabrata isolates.
Table 3. Characteristics and fluconazole treatment features in six patients infected with C. glabrata that acquired in vitro resistance under fluconazole treatment.
Table 4. Details of the relative expression levels and copy numbers of CDR1, PDH1, ERG11, and SNQ2 genes in four patients' genetically stable Candida glabrata sequential isolates that displayed distinct fluconazole MIC levels.
CgCDR1, CgPDH1, CgSNQ2, and CgERG11 expression
To investigate the mechanism of fluconazole resistance in C. glabrata we analyzed the gene expression levels of CgCDR1, CgPDH1, CgERG11, and CgSNQ2. In the 43 C. glabrata isolates, the normalized gene expression levels were as follows: CDR1 levels ranged from 0.010 to 92.61, median 0.415 IQR [0.120–1.825]; CDR2 levels ranged from 0.04 to 375.57, median 0.91 IQR [0.46–3.59]; SNQ2 levels ranged from 0.010 to 73.61, median 0.575 IQR [0.34–1.72] and ERG11 levels ranged from 0 to 4.34, median 0.54 IQR [0.27–1.00] (Table 2).
When considering the nine fluconazole resistant isolates, CgCDR1 was overexpressed in five isolates, CgPDH1 was overexpressed in six isolates, CgSNQ2 was overexpressed in three isolates and CgERG11 was overexpressed in one isolate (Table 2). Overall, the expression of the efflux pump genes CgCDR1, CgPDH1, and CgSNQ2 was increased in two resistant isolates compared with the susceptible initial isolates. One resistant isolate displayed simultaneous overexpression of CgCDR1, CgPDH1 and CgERG11, while another isolate displayed overexpression of both CgPDH1 and CgSNQ2. CgCDR1 and CgPDH1 were each overexpressed in two resistant isolates. In contrast, at least one gene was overexpressed in nine fluconazole susceptible isolates. Although overexpression of CgERG11 was never observed, CgPDH1, CgCDR1 and CgSNQ2 were overexpressed in seven, three and two fluconazole-susceptible isolates, respectively (Table 2). In four of these isolates, CgPDH1 overexpression was observed in parallel with either CgSNQ2 or CgCDR1 overexpression. When considering sequential isolates derived from the same patients, the acquisition of fluconazole resistance was associated with increased expression of CgCDR1 and CgPDH1 in patient 1 and increased expression of CgPDH1 in patient 4 (Table 4).
CgCDR1, CgPDH1, CgSNQ2, and CgERG11 Copy Number
To investigate one possible mechanism of C. glabrata fluconazole resistance we analyzed gene copy number of CgCDR1, CgPDH1, CgERG11, and CgSNQ2. The normalized gene copy number values were as follows: CDR1 copy numbers ranged from 1.02 to 4.70, median 1.71 IQR [1.2–2.18]; CDR2 copy numbers ranged from 0.00 to 5.70, median 1.43 IQR [1.14–1.98]; SNQ2 copy numbers ranged from 0.18 to 364.12, median 1.22 IQR [0.77–1.61] and ERG11 copy numbers ranged from 0.55 to 1549.94, median 1.51 IQR [1.13–2.20]. In some isolates, either overexpression or increased gene copy number was observed for CgCDR1 and CgPDH1 (Table 2). Nevertheless, three out of five resistant isolates that displayed CgCDR1 overexpression also displayed an increase in CgCDR1 gene copy number (4.7, 4.04 and 3.64 copies for isolates 1.2, 6.2, and 8, respectively) (Table 2). Three of the six resistant isolates that displayed CgPDH1 overexpression also exhibited an increase in CgPDH1 gene copy number (2.71, 4.13, and 5.7 copies in isolates 7, 1.2, and 6.2, respectively). Similarly, resistant isolates displaying CgERG11 overexpression also exhibited an increase in CgERG11 gene copy number (12.7 times more than TU10). For CgSNQ2, gene overexpression did not correlate with gene copy number (Table 2). Among the fluconazole susceptible isolates, the copy number of at least one gene was increased in seven isolates. More precisely, copy number of the four genes was increased in isolate 21.2, while only one gene copy number was increased in the remaining isolates as follows: CgERG11 copy number was increased in four isolates and the gene copy number of CgSNQ2 and CgCDR1 was increased, each in one isolate.
Population-Based Gene Interaction Study
Univariate correlation analysis indicated a complex interaction between gene copy number and expression levels (Table 5). A correlation between CgPDH1, CgCDR1, and CgSNQ2 gene copy numbers was observed, while a correlation between the expression levels of CgPDH1, CgERG11, and CgSNQ2 was also found. Fluconazole MIC values positively correlated with both CgCDR1 and CgSNQ2 expression levels. A positive correlation between gene copy number and expression level was observed only for the CgPDH1 gene. Figure 1 summarizes the significant combined main effects and first-level interactions between gene copy number and gene expression levels influencing fluconazole MIC levels in C. glabrata.
Table 5. Spearman rank correlation coefficients between the fluconazole MICs and the normalized expression level and copy number of CDR1, PDH1, ERG11, and SNQ2 genes.
Figure 1. Illustration of the effect of the complex interactions between CgPDH1, CgSNQ2, and CgERG11 drug resistance gene copy numbers and expression levels on fluconazole resistance in Candida glabrata clinical isolates. CgPDH1 gene copy number and expression were correlated, and their interaction led to an increase in fluconazole MIC levels. In contrast, the main effect of CgPDH1 overexpression independently led to a decrease in MIC levels. The main effect of CgSNQ2 overexpression independently promoted an increase in fluconazole MIC levels, while its interaction with CgERG11 expression led to a decrease in MIC levels.
The best-fitting multivariate model of the factors influencing fluconazole MIC values included (i) the main effects of the expression levels of CgSNQ2, CgERG11, and CgPDH1 and the gene copy number of CgPDH1 and (ii) the interaction between CgSNQ2 and CgERG11 gene expression and the interaction between CgPDH1 copy number and expression. In this model, C. glabrata fluconazole MICs were significantly increased due to both CgSNQ2 overexpression (p < 10−4) and the interaction between CgPDH1 copy number and expression (p = 0.038). In contrast, both CgPDH1 overexpression (p = 0.049) and the interaction between CgSNQ2 and CgERG11 expression (p = 0.003) led to a statistically significant decreased in fluconazole MICs. Neither CgERG11 expression (p = 0.823) nor CgPDH1 gene copy number (p = 0.985) independently displayed a statistically significant influence on fluconazole MICs. These multivariate analysis findings are illustrated in Figure 1.
This population-based approach allowed us to confirm the finding that ABC transporter family membrane gene overexpression is a pivotal fluconazole resistance mechanism in clinical C. glabrata isolates, which had been previously evidenced primarily using laboratory strains. Moreover, this study enabled us to further elucidate the main effects and complex interactions between gene expression and copy number associated with in vitro fluconazole resistance. The population-based analysis revealed that CgPDH1 and CgERG11 overexpression indirectly influence fluconazole MIC levels by increasing CDR1 and SNQ2 expression levels. Although CgERG11, CgSNQ2, and CgPDH1 expression and, to a lesser extent, copy number levels were significantly associated with fluconazole resistance, our findings highlight the complex interaction between these genes and suggest the involvement of unmeasured confounding effects. Indeed, mutations in CgPDR1, CgERG11 or others genes that might impact C. glabrata fluconazole resistance were not analyzed in this study. Univariate analysis (Table 5) suggested that CgCDR1 overexpression plays a major role in the regulation of MIC levels. However, its effect did not remain significant in the multivariate analysis, which highlighted the major independent effect of CgSNQ2 overexpression. Both univariate and multivariate analyses showed that copy number and overexpression were markedly correlated for CgPDH1 but not for the other genes. In patient 1, for whom sequential C. glabrata isolates displayed a stable genotype, switching to fluconazole resistance was associated with both overexpression and increased gene copy number of CgPDH1 (Table 4). A high correlation between the expression levels of all membrane transporter genes was detected, which is probably because they are controlled by the same transcription factor CgPDR1 (Ferrari et al., 2009; Paul et al., 2011). This transcription factor has been associated with both resistance to azole antifungals via the overexpression of ABC transporter genes and enhanced virulence in C. glabrata due to gain-of-function mutations (Ferrari et al., 2009).
CgCDR1 and CgPDH1 overexpression-driven drug efflux is the most frequent fluconazole resistance mechanism observed in C. glabrata (Miyazaki et al., 1998; Sanglard et al., 1999, 2001; Wada et al., 2002; Coleman and Mylonakis, 2009). Numerous reports have documented acquired resistance upon exposure to long-term fluconazole therapy (Coste et al., 2007; Shin et al., 2007), as observed in the four patients who acquired resistant isolates (Table 3). Although CgCDR1 and CgPDH1 expression could not be determined in the fluconazole susceptible isolate, overexpression of both genes was associated with fluconazole resistance in the isolate derived from patient 2. Acquisition of both CgPDH1 overexpression and increased gene copy number was observed in isolates derived from patient 1 (Table 4).
An overproduction of lanosterol-demethylase in C. glabrata has initially been described by Marichal et al. who found that duplication of the entire chromosome E (Marichal et al., 1997), the chromosome on which ERG11 is located, was associated with a four-fold increase in ERG11 copy number in fluconazole resistant compared with susceptible isolates. A recent study has also shown that lanosterol-demethylase overproduction, for which the mechanism is unknown, was involved in fluconazole resistance (Redding et al., 2003). In our study, an increase in CgERG11 gene copy number was associated with CgERG11 overexpression in only one isolate. Furthermore, CgCDR1, CgPDH1, and/or CgERG11 gene copy number was inconsistently associated with overexpression, thereby suggesting that an increase in gene copy number might sometime result in gene silencing via either heterochromatin formation or DNA methylation (Merrick and Duraisingh, 2006; Mishra et al., 2011).
A strong correlation between CgCDR1, CgPDH1, and CgERG11 gene copy number was observed (Table 4), although these genes are located on chromosomes M, F, and E, respectively. Further genome hybridization experiments may shed light on the underlying mechanisms and enable the localization of these amplified genes. Many gene duplication mechanisms have been described in various regions of the C. glabrata genome. Poláková et al. have suggested that duplication occurs in chromosomal segments that encode for some ABC transporter genes (Poláková et al., 2009). Similarly, Mûller et al. have shown that the mechanisms of microevolution of this asexual species include recombination within tandem arrays of repeated genes, which often encode cell wall proteins, thereby suggesting a possible adaptive function (Mûller et al., 2009). The fact that some genotypically identical isolates (genotype 3 and 7, Table 2) displayed different fluconazole resistance mechanisms support the hypothesis that fluconazole resistance is chiefly associated with the overexpression of normally present genes within the C. glabrata genome rather than gene mutation.
In conclusion, this population-based approach provides further evidence that overexpression of membrane transporter genes is a pivotal fluconazole resistance mechanism in C. glabrata clinical isolates. A high correlation between the expressions of these genes, which are controlled by the same transcription factor, was observed. However, the multivariate analysis revealed, for the first time, that CgSNQ2 exerts a major independent effect on fluconazole resistance. Regarding the other genes assessed, interactions played a greater role than the main effects in fluconazole resistance. Notably, CgPDH1 overexpression was significantly associated with increased gene copy number, which was not observed with the other genes analyzed. Furthermore, a correlation between gene copy number and fluconazole resistance was only observed for CgERG11, thereby suggesting that an increase in copy number of C. glabrata membrane transporter genes might result in gene silencing. This study provides only partial insight into C. glabrata fluconazole resistance mechanisms, as many other mechanisms have not been analyzed, including those targeting other membrane components. Nevertheless, our findings highlight the important role that complex gene copy number and overexpression interactions play in fluconazole resistance, which should been taken into account in further studies.
Salma Abbes carried out the collection of the clinical isolates, the molecular genetic studies and drafted the manuscript. Charles Mary participated in the design of the molecular genetic studies and drafted the manuscript. Hayet Sellami participated in the study design and isolates and clinical data collection. Annie Michel-Nguyen participated in C. glabrata isolates' identification and in vitro susceptibility assays analysis. Stéphane Ranque performed the statistical analysis. Ali Ayadi and Stéphane Ranque conceived the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Salma Abbes received travel grants from Aix Marseille University. We gratefully acknowledge Sandra Moore for English editing.
Abbes, S., Sellami, H., Sellami, A., Hadrich, I., Amouri, I., Mahfoudh, N., et al. (2012). Candida glabrata strain relatedness by new microsatellite markers. Eur. J. Clin. Microbiol. Infect. Dis. 31, 83–91. doi: 10.1007/s10096-011-1280-4
Abbes, S., Sellami, H., Sellami, A., Makni, F., Mahfoudh, N., Makni, H., et al. (2011). Microsatellite analysis and susceptibility to FCZ of Candida glabrata invasive isolates in Sfax Hospital, Tunisia. Med. Mycol. 49, 10–15. doi: 10.3109/13693786.2010.493561
Amouri, I., Sellami, H., Borji, N., Abbes, S., Sellami, A., Cheikhrouhou, F., et al. (2011). Epidemiological survey of vulvovaginal candidosis in Sfax, Tunisia. Mycoses 54, e499–e505. doi: 10.1111/j.1439-0507.2010.01965.x
Bouchara, J. P., Zouhair, R., Le Boudouil, S., Renier, G., Filmon, R., Chabasse, D., et al. (2000). In-vivo selection of an azole-resistant petite mutant of Candida glabrata. J. Med. Microbiol. 49, 977–984.
Cernicka, J., and Subik, J. (2006). Resistance mechanisms in fluconazole-resistant Candida albicans isolates from vaginal candidiasis. Int. J. Antimicrob. Agents 27, 403–408. doi: 10.1016/j.ijantimicag.2005.12.005
Chau, A. S., Mendrick, C. A., Sabatelli, F. J., Loebenberg, D., and McNicholas, P. M. (2004). Application of real-time quantitative PCR to molecular analysis of Candida albicans strains exhibiting reduced susceptibility to azoles. Antimicrob. Agents Chemother. 48, 2124–2131. doi: 10.1128/AAC.48.6.2124-2131.2004
Clinical and Laboratory Standards Institute. (2008). Reference Method for Broth Dilution Antifungal Susceptibility Testing of Yeasts. Approved standard CLSI document M27-A3. Wayne, PA: Clinical and Laboratory Standards Institute.
Coste, A., Selmecki, A., Forche, A., Diogo, D., Bougnoux, M.-E., d' Enfert, C., et al. (2007). Genotypic evolution of azole resistance mechanisms in sequential Candida albicans isolates. Eukaryot. Cell. 6, 1889–1904. doi: 10.1128/EC.00151-07
Diekema, D. J., Messer, S. A., Brueggemann, A. B., Coffman, S. L., Doern, G. V., Herwaldt, L. A., et al. (2002). Epidemiology of candidemia: 3-year results from the emerging infections and the epidemiology of Iowa organisms study. J. Clin. Microbiol. 40, 1298–1302. doi: 10.1128/JCM.40.4.1298-1302.2002
Ferrari, S., Ischer, F., Calabrese, D., Posteraro, B., Sanguinetti, M., Fadda, G., et al. (2009). Gain of function mutations in CgPDR1 of Candida glabrata not only mediate antifungal resistance but also enhance virulence. PLoS. Pathog. 5:e1000268. doi: 10.1371/journal.ppat.1000268
Marichal, P., Vanden Bossche, H., Odds, F. C., Nobels, G., Warnock, D. W., Timmerman, V., et al. (1997). Molecular biological characterization of an azole-resistant Candida glabrata isolate. Antimicrob. Agents Chemother. 41, 2229–2237.
Mishra, P. K., Baum, M., and Carbon, J. (2011). DNA methylation regulates phenotype-dependent transcriptional activity in Candida albicans. Proc. Natl. Acad. Sci. U.S.A. 108, 11965–11970. doi: 10.1073/pnas.1109631108
Miyazaki, H., Miyazaki, Y., Geber, A., Parkinson, T., Hitchcock, C., Falconer, D. J., et al. (1998). Fluconazole resistance associated with drug efflux and increased transcription of a drug transporter gene, PDH1, in Candida glabrata. Antimicrob. Agents. Chemother. 42, 1695–1701.
Mûller, H., Thierry, A., Coppée, J.-Y., Gouyette, C., Hennequin, C., Sismeiro, O., et al. (2009). Genomic polymorphism in the population of Candida glabrata: gene copy-number variation and chromosomal translocations. Fungal. Genet. Biol. 46, 264–276. doi: 10.1016/j.fgb.2008.11.006
Niimi, M., Nagai, Y., Niimi, K., Wada, S. I., Cannon, R. D., Uehara, Y., et al. (2002). Identification of two proteins induced by exposure of the pathogenic fungus Candida glabrata to fluconazole. J. Chromatogr. B. Analyt. Technol. Biomed. Life. Sci. 782, 245–252. doi: 10.1016/S1570-0232(02)00668-2
Pfaller, M. A., and Diekema, D. J. (2004). Twelve years of fluconazole in clinical practice: global trends in species distribution and fluconazole susceptibility of bloodstream isolates of Candida. Clin. Microbiol. Infect. 10(Suppl. 1), 11–23. doi: 10.1111/j.1470-9465.2004.t01-1-00844.x
Pfaller, M. A., Messer, S. A., Hollis, R. J., Boyken, L., Tendolkar, S., Kroeger, J., et al. (2009). Variation in susceptibility of bloodstream isolates of Candida glabrata to fluconazole according to patient age and geographic location in the United States in 2001 to 2007. J. Clin. Microbiol. 47, 3185–3190. doi: 10.1128/JCM.00946-09
Poláková, S., Blume, C., Zárate, J. A., Mentel, M., Jørck-Ramberg, D., Stenderup, J., et al. (2009). Formation of new chromosomes as a virulence mechanism in yeast Candida glabrata. Proc. Natl. Acad. Sci. U.S.A. 106, 2688–2693. doi: 10.1073/pnas.0809793106
Redding, S. W., Kirkpatrick, W. R., Saville, S., Coco, B. J., White, W., Fothergill, A., et al. (2003). Multiple patterns of resistance to fluconazole in Candida glabrata isolates from a patient with oropharyngeal candidiasis receiving head and neck radiation. J. Clin. Microbiol. 41, 619–622. doi: 10.1128/JCM.41.2.619-622.2003
Sanglard, D., Ischer, F., and Bille, J. (2001). Role of ATP-binding-cassette transporter genes in high-frequency acquisition of resistance to azole antifungals in Candida glabrata. Antimicrob. Agents Chemother. 45, 1174–1183. doi: 10.1128/AAC.45.4.1174-1183.2001
Sanglard, D., Ischer, F., Calabrese, D., Majcherczyk, P. A., and Bille, J. (1999). The ATP binding cassette transporter gene CgCDR1 from Candida glabrata is involved in the resistance of clinical isolates to azole antifungal agents. Antimicrob. Agents Chemother. 43, 2753–2765.
Sanguinetti, M., Posteraro, B., Fiori, B., Ranno, S., Torelli, R., and Fadda, G. (2005). Mechanisms of azole resistance in clinical isolates of Candida glabrata collected during a hospital survey of antifungal resistance. Antimicrob. Agents Chemother. 49, 668–679. doi: 10.1128/AAC.49.2.668-679.2005
Sellami, A., Sellami, H., Néji, S., Makni, F., Abbes, S., Cheikhrouhou, F., et al. (2011). Antifungal susceptibility of bloodstream Candida isolates in Sfax hospital: Tunisia. Mycopathologia 171, 417–422. doi: 10.1007/s11046-010-9388-0
Shin, J. H., Chae, M. J., Song, J. W., Jung, S.-I., Cho, D., Kee, S. J., et al. (2007). Changes in karyotype and azole susceptibility of sequential bloodstream isolates from patients with Candida glabrata candidemia. J. Clin. Microbiol. 45, 2385–2391. doi: 10.1128/JCM.00381-07
Torelli, R., Posteraro, B., Ferrari, S., La Sorda, M., Fadda, G., Sanglard, D., et al. (2008). The ATP-binding cassette transporter-encoding gene CgSNQ2 is contributing to the CgPDR1-dependent azole resistance of Candida glabrata. Mol. Microbiol. 68, 186–201. doi: 10.1111/j.1365-2958.2008.06143.x
Wada, S. I., Niimi, M., Niimi, K., Holmes, A. R., Monk, B. C., Cannon, R. D., et al. (2002). Candida glabrata ATP-binding cassette transporters CDR1p and PDH1p expressed in a Saccharomyces cerevisiae strain deficient in membrane transporters show phosphorylation-dependent pumping properties. J. Biol. Chem. 277, 46809–46821. doi: 10.1074/jbc.M207817200
Keywords: Candida glabrata, fluconazole, resistance mechanisms, gene expression, gene copy number, quantitative real-time PCR, human infection
Citation: Abbes S, Mary C, Sellami H, Michel-Nguyen A, Ayadi A and Ranque S (2013) Interactions between copy number and expression level of genes involved in fluconazole resistance in Candida glabrata. Front. Cell. Infect. Microbiol. 3:74. doi: 10.3389/fcimb.2013.00074
Received: 06 September 2013; Accepted: 18 October 2013;
Published online: 11 November 2013.
Edited by:Kunihiko Nishino, Institute of Scientific and Industrial Research, Japan
Reviewed by:Margaret E. Bauer, Indiana University School of Medicine, USA
Janio M. Santurio, Universidade Federal de Santa Maria, Brazil
Copyright © 2013 Abbes, Mary, Sellami, Michel-Nguyen, Ayadi and Ranque. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Stéphane Ranque, Laboratoire de Parasitologie-Mycologie, AP-HM Timone, 264 rue Saint Pierre, 13005 Marseille, France e-mail: firstname.lastname@example.org