INTRODUCTION
Superficial and life-threatening systemic fungal infections have been increasing over the last 2 decades. The majority of systemic fungal infections are caused by
Candida,
Aspergillus, and
Cryptococcus species, but
Candida albicans and non-
albicans Candida species still account for most of the infections. A few treatment options exist in medical practice, including the use of at least four antifungal chemical classes (azoles, candins, pyrimidine analogues, and polyenes). Emergence of antifungal resistance is a consequence of long-term use of these agents, which is occurring in most immunocompromised patients with HIV or undergoing organ transplants or cancer chemotherapy (
1). Clinical criteria can define antifungal resistance, and this has been achieved by the setting of
Clinical
Break
Points (CPB) which indicate a drug concentration for a given fungal pathogen above/under which failure/success of a therapy can be expected (
2). For example and according to these criteria, antifungal resistance for azoles is currently the highest for
C. glabrata among other
Candida spp. and accounts for 10 to 20% of the
C. glabrata population (
3,
4). This yeast species is ranked as second after
C. albicans among bloodstream isolates. Recent studies report in several institutions an epidemiological shift of
C. glabrata at the expense of
C. albicans, but the reasons behind this phenomenon are still unexplained (
4).
Antifungal resistance involves different mechanisms, including principally enhanced drug efflux and target alterations by active-site mutations or overexpression (
5). In
C. glabrata, one of the prominent resistance mechanisms invokes the participation of multidrug transporters of the ABC transporter family (
CgCDR1,
CgCDR2, and
CgSNQ2) (
6–8). We and others have found that these transporters are upregulated in azole-resistant isolates with an MIC higher than 16 μg/ml for fluconazole. The upregulation is associated with mutations, so-called gain-of-function (GOF) mutations, in a transcription factor of the Zn2-Cys6 family,
CgPDR1 (
1,
9–13). In
C. albicans, antifungal resistance is multifactorial and involves the participation of efflux transporters, target mutations, compensatory mutations, and genome rearrangements (
5). As in the case of
C. glabrata, mutations in the transcription factors
TAC1,
MRR1, and
UPC2 result in the upregulation of target genes participating to the development of azole resistance (
14–18). The resistance levels achieved by
C. glabrata and
C. albicans address the need to overcome and avoid this phenomenon. Several concepts have been proposed in the past and utilize as basic principle the combination of one antifungal with another compound in order to increase antifungal activity (
19,
20). Given the importance of ABC-transporters for the development of azole resistance both in
C. glabrata and
C. albicans, one possible option to tackle resistance could be the use of transporter inhibitors. In cancer research, human P-gp, which are functional homologs of fungal ABC transporters, are important mediators of resistance to many anticancer drugs. A wide range of P-gp-inhibitory compounds identified include natural and synthetic polymers, P-gp substrates (such as FK506), calcium channel modulator (verapamil), calmodulin inhibitors, and quinine analogs (
21). Consequently, implementing the concept of combination therapy using an efflux inhibitor in order to increase the activity of fluconazole, a drug still widely used to treat fungal infections, seems an adequate strategy for overcoming resistance development in
C. glabrata.
Drug resistance acquisition can also be associated to fitness costs in several microbial systems (
22). In
C. glabrata, however, we found that azole resistance was on the contrary resulting in enhanced virulence and fitness in animal models. This feature was dependent on the presence of a GOF mutation in
CgPDR1 (
9). Recently, we found that this effect could be mediated partially by the ABC transporter
CgCDR1, which is upregulated to variable levels in azole-resistant isolates in
C. glabrata (
23).
Since CgCDR1 plays an important role in the development of azole resistance and that it can also contribute to increase virulence and fitness of C. glabrata, we reasoned that ABC transporters inhibitors such as milbemycins could be of potential therapeutic interest. We demonstrate here that specific milbemycin derivatives not only exhibit expected activities as efflux inhibitors but also that they possess intrinsic antifungal and fungicidal activities. Here we show that experimental treatments of C. glabrata infections by combination therapy are feasible and expanded this idea to C. albicans infections. Lastly, we perform transcriptional profiling of both Candida species exposed to milbemycins in order to understand the basis for their unexpected antifungal activity.
MATERIALS AND METHODS
Strains, media, and drugs.
The
C. albicans strains used in the present study are listed in
Table 1. Yeast strains were grown in liquid YEPD complete medium (1% Bacto peptone [Difco], 0.5% yeast extract [Difco], 2% glucose [Fluka]). To grow the strains on solid media, 2% agar (Difco) was added.
Escherichia coli DH5α was used as a host for plasmid construction and propagation. DH5α cells were grown in Luria-Bertani (LB) broth or on LB plates, which were supplemented with ampicillin (0.1 mg/ml) when required. Fluconazole was obtained from Sigma. Milbemycins were obtained from Novartis Animal Health (Basel, Switzerland).
Drug susceptibility testing.
Susceptibility assays were performed according to the standard broth microdilution protocols Edef. 7.1 (Subcommittee on Antifungal Susceptibility Testing of the ESCMID European Committee for Antimicrobial Susceptibility Testing [AFST-EUCAST]) (
28). Briefly, serial 2-fold dilutions of fluconazole in RPMI 1640 broth (with
l-glutamine, without bicarbonate and with phenol red as the pH indicator; Sigma), supplemented with 2%, (wt/vol) of
d-glucose for Edef. 7.1, were distributed in 50-μl volumes at four times the final desired concentration into the wells of flat-bottom microtiter plates. Fluconazole final concentrations ranged from 128 to 0.25 μg/ml. Cell suspensions were prepared in sterile saline solution from overnight cultures of yeast strains at 35°C in Sabouraud dextrose agar plates. The suspensions were diluted in the test medium and added in 150-μl volumes to the drug solutions in the microtiter plates to yield final inoculum sizes of 0.5 to 2.5 × 10
5 cells/ml. Drug-free cultures and sterility controls were included in each microtiter plate. The plates were incubated at 35°C, and the MIC values were read after 24 h by measuring the absorbance using a spectrophotometric microdilution plate reader set at 450 nm. Checkerboard tests were performed with a similar procedure but including serial dilutions of milbemycins in a 50-μl volume at four times the final desired concentration for each fluconazole-containing well and with final concentrations ranging from 0.8 to 25.6 μg/ml. Milbemycins were first diluted in dimethyl sulfoxide (DMSO) and diluted accordingly. DMSO concentration in final suspensions was 0.8%. The cell suspensions were added in 100-μl volumes to the drug solutions to obtain final inoculum sizes of 0.5 to 2.5 × 10
5 cells/ml.
The fractional inhibitory concentration (FIC) index for a drug combination in the checkerboard method is the minimum ΣFIC obtained with the following equation: ΣFIC = FIC
A + FIC
B = (C
FLC/MIC
FLC) + (C
Mil/MIC
Mil), where MIC
FLC and MIC
Mil are the MICs of fluconazole and the tested milbemycin when tested alone, and C
FLC and C
Mil are the concentrations of fluconazole and the tested milbemycin resulting in at least 50% growth inhibition in a given well in the 96-well microtiter plate. Interactions were categorized according to the method of Te Dorsthorst et al. (
29) by the following rules: synergism (FIC, ≤0.5), additivity (FIC, >0.5 to ≤1), and indifference (FIC, >1 to <4).
Efflux of rhodamine 6G.
A whole cells rhodamine 6G (R6G) efflux assay, adapted from a previously developed protocol, was used to measure the drug efflux capacity of
C. albicans and
C. glabrata isolates (
30). Each fungal species required a specific procedure for cell preparation. For
C. albicans, cell cultures grown overnight in YEPD were diluted in 5 ml of fresh medium and allowed to grow at 30°C under constant agitation to a density of 2 × 10
7 cells/ml. Cells were centrifuged, washed in 5 ml of phosphate-buffered saline (PBS; pH 7), and resuspended in 2 ml of PBS. These suspensions were incubated for 1 h at 30°C under constant agitation to energy-deprived cells, and R6G was next added at a concentration of 10 μg/ml. The incubation was continued for 1 h to allow R6G accumulation. After this incubation time, cells were sedimented by centrifugation, washed with PBS at 4°C, and resuspended in a final volume of 300 μl of PBS. For
C. glabrata, overnight cultures were diluted to 2 × 10
7 cells/ml in 50 ml of YEPD and agitated for 4 h at 30°C. Cells were next washed twice with PBS and resuspended in PBS at a 2% (wet weight) concentration. R6G was added to a 10 μM end concentration together with 2 mM deoxyglucose. After 1 h of incubation at 30°C, the cells were washed twice with PBS and resuspended to a concentration of 10
8 cells/ml.
Fifty-microliter portions of individual suspensions were diluted in 150 μl of PBS and divided into aliquots into a 96-well microtiter plate, which was then placed in a SpectraMax Gemini fluorimeter (Molecular Devices, Sunnyvale, CA, USA) with a temperature control set at 30°C. Baseline emission of fluorescence (excitation wavelength, 344 nm; emission wavelength, 555 nm) was recorded for 5 min as relative fluorescence units (RFU), and d-glucose was next added to each strain at a final concentration of 1% (wt/vol) to initiate R6G efflux. As negative controls, no glucose was added to a series of separate aliquots of each strain. The data points were recorded in duplicate for 60 min in 1-min intervals.
Efflux inhibition with milbemycins was carried out by adding the different compounds to each yeast sample. The milbemycin-yeast cell mixture was incubated for 10 min at room temperature before initiation of efflux experiments.
MitoSOX Red staining.
Intracellular reactive oxygen species (ROS) production was examined using MitoSOX Red (Molecular Probes). MitoSOX Red is a lipid soluble cation that accumulates in the mitochondrial matrix, where it can be oxidized to a fluorescent product by superoxide. Yeast strains from initial concentration of 2 × 106 cells/ml were grown in RPMI medium containing 10 μg of A3Ox/ml. After 1 h of incubation, 500-μl aliquots were washed twice with PBS and incubated in the dark for 20 min in 2.5 μM MitoSOX Red. The cells were washed three times with PBS and resuspended in PBS, and the percentage of cells positively stained with MitoSOX Red was determined by fluorescence microscopy using a Zeiss Axioplan 2 fluorescence microscope. Images were recorded using a Visitron Systems HistoScope Camera and VisiView Imaging Software. Fluorescence of cells was also determined using a FLUOstar OMEGA microplate reader (BMG) with excitation and emission wavelengths of 510 and 580 nm, respectively.
Construction of Candida microarrays.
C. glabrata microarrays were designed according to the Agilent eArray Design guidelines as previously described (
23), and custom arrays were manufactured in the 8×15 K format by Agilent Technologies.
C. albicans microarrays were also manufactured by Agilent Technologies; however, the array design (design ID 037331) was obtained from Synnott et al. (
31). A total of 6,101 genes (including 12 mitochondrial genes) are represented by two sets of probes, both spotted in duplicate. Four copies of each array were printed on a 4×44 K format.
cRNA synthesis, one-color labeling, and array hybridization.
Sample preparation was performed on biological quadruplicate cultures of the
C. glabrata strain DSY562 and biological triplicate for the
C. albicans strain SC5314. Log-phase cultures at 35°C with agitation in RPMI 1640 medium with 2%
d-glucose were treated for 1 h in the same conditions with either 10 μg of milbemycin A3 oxim/ml or the DMSO solvent. Total RNA was extracted after mechanical disruption of the cells with glass beads by a phenol-chloroform-lithium chloride procedure, as previously described (
32). The integrity of the input RNA template was determined prior to labeling or amplification using an Agilent RNA 6000 Nano LabChip kit and 2100 BioAnalyzer (Agilent Technologies). Agilent's One-Color Quick Amp labeling kit (Agilent Technologies) was used to generate fluorescent cRNA as previously described (
32). Briefly, a spike mix and T7 promoter primers were added to 400 ng of total RNA from each sample. cDNA synthesis was promoted by Moloney murine leukemia virus reverse transcriptase in the presence of deoxynucleoside triphosphates and RNaseOUT. Next, cRNA was produced from this first reaction with T7 RNA polymerase, which simultaneously amplifies target material and incorporates cyanine 3-labeled CTP. The labeled cRNAs were purified with an RNeasy minikit (Qiagen) and quantified using a NanoDrop ND-1000 UV/VIS spectrophotometer. A total of 600 ng of Cy3-labeled cRNAs from each sample were fragmented and hybridized for 17 h at 65°C to specific subarrays of the 8×15 K or 4×44 K format using a gene expression hybridization kit (Agilent Technologies) and a gasket slide.
Microarray data analysis.
Slides were washed and processed according to the Agilent 60-mer Oligo microarray processing protocol and scanned on an Agilent microarray scanner G2565BA (Agilent Technologies). The data were extracted from the images with feature extraction (FE) software (Agilent Technologies). Each slide was processed with spike quality controls to establish the dynamic range of signals which fitted a minimum R2 value of 0.99. The FE software flags outlier features (the percentages ranged from 0.01 to 0.15% of signals overall for cRNA hybridizations) and detects and removes spatial gradients and local backgrounds. The data were normalized using a combined rank consistency filtering with Lowess intensity normalization. The gene expression values obtained with FE software were imported into GeneSpring 11 software (Agilent Technologies) for preprocessing and data analysis. For inter-array comparisons, a linear scaling of the data was performed using the quantile normalization of one-color signal values of noncontrol probes on the microarray. The expression of each gene was normalized by its median expression across all samples. Statistical analyses were performed using unpaired t tests and a corrected P value of 0.05 was chosen as the cutoff for significance. Changes in expression for each gene of at least 2-fold between treated and nontreated samples were considered significant. Microarray data can be retrieved at the Gene Expression Omnibus NCBI site under the accession number GSE40232. Validation of gene expression results from microarray analysis was performed by quantitative real-time PCR (qPCR) as detailed below.
qPCR.
RNA sample preparation was performed as described above for microarray experiments. RNA samples were treated with RNase-free DNase (DNA-free; Ambion) to remove any contaminating genomic DNA and reverse transcribed using random hexamers as the priming method (Transcriptor First-Strand cDNA synthesis kit; Roche). Expression levels of
C. glabrata target genes (
ERG4,
MET3,
YPS1, and
YPS3) and the
C. albicans genes (
HXT6,
WOR1, and
OPT5) were determined by a SYBR green-based quantitative real-time PCR (MesaBlue qPCR kit for SYBR assay; Eurogentec) in a StepOne Plus real-time thermal cycler (Applied Biosystems). Primers used for
YPS1 and
YPS3 quantification were previously described (
23). Primers used for quantification of the remaining genes are listed in
Table 2. Relative gene expression in milbemycin A3 oxim-treated samples in comparison to nontreated controls was determined from
CgTEF3-normalized expression levels for
C. glabrata and
ACT1 for
C. albicans.
Animal experiments.
Female BALB/c mice (20 to 25 g) were purchased from Harlan Italy S.r.l (San Pietro al Natisone, Udine, Italy) and inbred in-house. The mice were housed in filter-top cages with free access to food and water. To establish C. glabrata infection, mice were injected into their lateral vein with saline suspensions of the C. glabrata strains (each in a volume of 200 μl).
Groups of 10 mice were established for each yeast strain. For tissue burden experiments, immunocompetent mice were inoculated with 4 × 107 CFU of each C. glabrata strains and 7 × 104 CFU of each C. albicans strains. After 7 days, mice were sacrificed by use of CO2 inhalation, and target organs (spleen and kidney) were excised aseptically, weighted individually, and homogenized in sterile saline by using a Stomacher 80 device (Pbi International, Milan, Italy) for 120 s at high speed. Organ homogenates were diluted and plated onto YPD. Colonies were counted after 2 days of incubation at 30°C, and the numbers of CFU g of organ−1 were calculated. CFU counts were analyzed with nonparametric Wilcoxon rank-sum tests. A P value of <0.05 was considered to be significant.
The milbemycin injectable solution was prepared by adding 1 volume of PEG 400 to 2 volumes of a milbemycin DMSO solution (0.2 mg/ml) in order to obtain a 0.1-mg/ml solution with DMSO/PEG 400 ratio of 2:1. The DMSO/PEG solution was injected in a 100-μl volume intraperitoneally each day. Fluconazole was injected intraperitoneally at dosages of 100 mg/kg/day dosage for C. glabrata and of 3 mg/kg/day for C. albicans.
DISCUSSION
In this work, we demonstrated that milbemycins are effective agents to treat
C. albicans and
C. glabrata infections in combination with fluconazole in animal models and even reverse azole resistance to levels of wild-type isolates. We showed that milbemycins A3 and A4 and especially their oxim derivatives have potent drug efflux inhibition activities. Several studies already established that this family of substances inhibit fungal ABC transporters (
20,
21,
50); however, it is the first time, to our knowledge, that this is demonstrated for oxim derivatives. These substances are active components of commercial antiparasitic drugs and thus are easy to obtain in large quantities.
One of the most remarkable properties of milbemycin oxims was their antifungal and fungicidal activities both in
C. albicans and
C. glabrata, which to our knowledge has not been reported yet for other efflux inhibitors. This activity seems to be mediated by the oxim prosthetic group, since A4 and A3 were not exhibiting the same activity. Furthermore, A3Ox was the most potent compound in both
Candida species. Thus, both the difference in chemical structure between A3 and A4 (
Fig. 8, one ethyl and methyl group at position 25, respectively) and the addition of an oxim group (
Fig. 8, at position 5) (
51) conferred a strong antifungal activity. Milbemycin oxims carry two activities, one including the inhibition of drug efflux at low concentration (efflux inhibition was detectable starting from 0.1 μg/ml) and the other including antifungal activity at higher drug concentrations (3.2 to 6.4 μg/ml). The antifungal activity is independent of the presence or absence of ABC transporters, since mutants lacking major ABC transporters involved in azole resistance (CgCdr1/CgCdr2, Cdr1, and Cdr2) exhibit similar A3Ox MIC values to wild-type isolates (data not shown). Interestingly, efflux inhibition is more potent against drug-resistant
C. albicans than
C. glabrata, as reflected by the respective inhibition of R6G efflux (
Vmax IC
50s are >20-fold higher for
C. glabrata,
Fig. 2). This is possibly due to structural differences of ABC transporters between the two species which may result in different inhibitor affinities. The intrinsic antifungal activity of A3Ox, on the other hand, shows no difference between the two types of isolates, therefore highlighting that it is mediated by an ABC-transporter-independent mechanism of action.
With the idea to better understand how the fungicidal effect of A3Ox was exerted on
C. albicans and
C. glabrata, we performed whole-genome transcript profiling with the same drug concentration and experiment duration. The results obtained revealed, as expected, several features common to drug stress that include the mobilization of oxidoreductive processes, of chaperoning of misfolded proteins and of protein degradation pathways. The upregulation of the caspase-like gene
MCA1 in
C. albicans could be related to the fungicidal effect of A3Ox. It has been shown that a striking caspase activity pattern can occur in
C. albicans during oxidative stress-induced programmed cell death (PCD) (
52). Together with the ROS-producing activity of A3Ox,
MCA1 upregulation suggests that A3Ox has the capacity to induce PCD in
C. albicans, although it can be expected that cell necrosis would be an important readout of A3Ox exposure.
One surprising observation was the apparent divergent transcriptional response of both
Candida species upon A3Ox exposure. We found that 204 gene orthologs between
C. glabrata (of 887 genes) and
C. albicans (of 1,809 genes) were commonly regulated by A3Ox, which reflects a common response mechanism, but also demonstrated quite different effects of this drug on the transcriptome of the two yeast species (see File S3 in the supplemental material). The cluster analysis of these 204 regulated genes showed that 62 were regulated in the opposite way between the two yeast species (
Fig. 9). The remaining genes that were commonly up- and downregulated could constitute a common core of genes responding to A3Ox in a species-independent way. The GO term analysis of the 144 coregulated genes revealed enrichments of genes involved in oxidoreductive processes and in oxidative stress response, which is consistent with a stress situation (see File S3 in the supplemental material). From the annotated genes involved in oxidoreductive processes, 75% belonged to genes commonly upregulated by A3Ox between both species (
Fig. 9). Oxidative stress is a signature common in response to xenobiotics in yeast (
53–55), and our transcript profile analyses therefore underline a similar effect for A3Ox. In addition to this oxidative stress response, the presence of chaperone genes (HSP gene family), of genes involved in protein ubiquitination process, and of genes involved in ESCRT vesicle trafficking was reminiscent of the response signature involving protein damage (
55). Interestingly, MAP kinases from both
C. albicans (
MKC1) and
C. glabrata (
SLT2) were also upregulated upon A3Ox treatment. Besides the role of these kinases in cell wall integrity as described above, it has been suggested that at least
MKC1 participates to oxidative stress response (
56), which is again in agreement with the idea that A3Ox stress responses involve oxidative processes.
At this stage, it is not possible to predict direct cellular targets of this substance in addition to efflux transporters in the tested
Candida isolates. A few studies comparing the effect of the same drug on different yeast species using microarrays are available. These studies have investigated the effect of antifungal drugs on fungal metabolism and generally could identify common response patterns between the different investigated species (
1). Here, we identified a limited overlap of gene expression between
C. glabrata and
C. albicans, thus suggesting species-specific response signatures. This could be due to the evolutionary distance between both yeast species, which determines how a specific pathogen behaves in contact with exogenous substances. However, even if species are more closely related, drug response can be still different as illustrated by transcriptome comparisons between
C. glabrata and
S. cerevisiae in the presence of benomyl (
57). The same conclusion was drawn for the response of diverse yeast species to azoles (
58).
Testing efflux inhibitors in animal models has been achieved in a few occasions. For example, Hayama et al. (
59) have used the
d-octapeptide derivative RC21v3, a Cdr1 inhibitor, in the treatment of murine oral candidiasis caused by azole-susceptible and azole-resistant
C. albicans clinical isolates. RC21v3 potentiated the therapeutic efficacy of fluconazole for mice infected with either strain. Sorensen et al. (
60) reported that milbemycin α
9 could potentiate fluconazole in an experimental systemic pyelonephritis with
C. albicans. We showed here that milbemycin oxims could potentiate fluconazole efficacy both in
C. albicans and
C. glabrata using a systemic model of infection, and thus we expand the concept of drug combination to both of these important fungal pathogens. We have used different milbemycin regimens ranging from 0.5 to 2.5 mg/kg/day in mice without signs of toxicity. These dosages correspond to those given in animal care (
61). A higher dosage (5 mg/kg/day) exhibited slight toxicity effects in treated mice without increasing fluconazole efficacy (data not shown). The serum levels obtained after 7-day treatments reached concentrations near to the MIC measured in
C. glabrata and
C. albicans. This could explain why milbemycin oxims on their own could decrease fungal burden in tissues infected with
C. glabrata and
C. albicans. On the other hand, we have shown in previous studies that the ABC transporter
CgCDR1 could participate to
C. glabrata virulence (
23). Thus, inhibition of this transporter by milbemycins could also impact on
C. glabrata virulence and therefore could contribute to the decrease of tissue burden upon inhibitor treatment.
Whether or not milbemycin oxims with fluconazole can be used in human is an open question. The primary targets of the substance are glutamate sensitive chloride channels in neurons and myocytes of invertebrates, leading to hyperpolarization of these cells and blocking of signal transfer. These glutamate channels are specific for invertebrates and are not expressed in mammalian hosts (
62). Some milbemycin derivatives have been used in humans already. For example, Cotreau et al. (
63) reported the use of moxidectin to treat onchocerciasis (river blindness). Onchocerciasis is a parasitic disease caused by the helminth
Onchocerca volvulus and is transmitted to humans through the bite of a black fly of the genus
Simulium. This study showed that such substances have low toxicity in humans when administered in the 3- to 36-mg range. Our preliminary data show that moxidectin acts synergistically with fluconazole in
C. glabrata and
C. albicans as in the case of milbemycin A3Ox. Therefore, the use of milbemycin for treating
Candida infections is potentially possible in humans.
Candida infections, and especially
C. glabrata infections, are rising in diverse countries and, together with azole resistance levels reached nowadays by
C. glabrata, alternative therapeutic approaches should be proposed in the future (
4).