- Research
- Open access
- Published:
Mechanism and bioinformatics analysis of the effect of berberine-enhanced fluconazole against drug-resistant Candida albicans
BMC Microbiology volume 24, Article number: 196 (2024)
Abstract
Biofilms produced by Candida albicans present a challenge in treatment with antifungal drug. Enhancing the sensitivity to fluconazole (FLC) is a reasonable method for treating FLC-resistant species. Moreover, several lines of evidence have demonstrated that berberine (BBR) can have antimicrobial effects. The aim of this study was to clarify the underlying mechanism of these effects. We conducted a comparative study of the inhibition of FLC-resistant strain growth by FLC treatment alone, BBR treatment alone, and the synergistic effect of combined FLC and BBR treatment. Twenty-four isolated strains showed distinct biofilm formation capabilities. The antifungal effect of combined FLC and BBR treatment in terms of the growth and biofilm formation of Candida albicans species was determined via checkerboard, time-kill, and fluorescence microscopy assays. The synergistic effect of BBR and FLC downregulated the expression of the efflux pump genes CDR1 and MDR, the hyphal gene HWP1, and the adhesion gene ALS3; however, the gene expression of the transcriptional repressor TUP1 was upregulated following treatment with this drug combination. Furthermore, the addition of BBR led to a marked reduction in cell surface hydrophobicity. To identify resistance-related genes and virulence factors through genome-wide sequencing analysis, we investigated the inhibition of related resistance gene expression by the combination of BBR and FLC, as well as the associated signaling pathways and metabolic pathways. The KEGG metabolic map showed that the metabolic genes in this strain are mainly involved in amino acid and carbon metabolism. The metabolic pathway map showed that several ergosterol (ERG) genes were involved in the synthesis of cell membrane sterols, which may be related to drug resistance. In this study, BBR + FLC combination treatment upregulated the expression of the ERG1, ERG3, ERG4, ERG5, ERG24, and ERG25 genes and downregulated the expression of the ERG6 and ERG9 genes compared with fluconazole treatment alone (p < 0.05).
Background
The infection rate of Candida albicans (C. albicans) has increased significantly, especially in immunocompromised patients and patients with organ intubation. According to statistics, systemic infection of C. albicans has a mortality rate of up to 31.8%, but the hospital intensive care unit is also the fourth highest risk factor for bloodstream infection [1, 2]. The rate of C. albicans infection was more than 50% in hospitals in the U.S., with a 40% mortality rate [3, 4]. The formation of biofilms can significantly reduce the sensitivity of fungi to azole drugs, resulting in drug resistance and sustained infection [5,6,7]. The formation of C. albicans biofilms can led to the production of extracellular matrix and a structure that encapsulates fungal cell populations, thereby hindering or delaying drug penetration; however, biofilms can also increase the expression of fungal efflux pump genes [8, 9]. FLC is a third-generation antifungal drug in the azole family that possesses a high degree of lipophilicity and low-toxicity and is a common first-line drug used to prevent and treat invasive C. albicans infections [10, 11]. In recent years, however, the formation of C. albicans biofilms has made the treatment of fungi ineffective, resulting in reduced of C. albicans susceptibility to fluconazole. This increase in drug resistance has hampered the penetration and absorption of azole drugs and reduced the number of antifungal drug targets [12,13,14,15]. In addition, mutations and the overexpression of ERG11 weaken the binding of the drug to the target enzyme. The upregulated expression of the coding efflux pump genes CDR1 and CDR2, which belong to the ATP family, and MDR can decrease the intracellular drug concentration in C. albicans, as well as the expression level of the gene HWP1, which encodes the mycelium cell wall-specific protein of C. albicans when the biofilm is formed [16,17,18]. At present, fluconazole can only inhibit the growth of fungal flora, not prevent it, and the emergence of fluconazole-resistant strains, as they can lead to repeated infections and hinder drug delivery, has presented a challenge. Accordingly, some researchers have aimed to identify drugs and treatment strategies that increase the sensitivity of C. albicans to fluconazole.
Berberine (BBR) is a common isoquinoline alkaloid found in many plants, particularly Berberis vulgaris. According to Ayurvedic and Chinese historical records dating back at least 3000 years, BBR was used to treat disease, with its medicinal properties having been reported to include significant antimicrobial activity against viral infections, fungi and inflammation. In terms of its clinical uses, BBR has been shown to help treat multiple diseases in vivo, such as tumors, diabetes, hyperlipidemia, and antiarrhythmic and immuno-suppressive diseases [19,20,21,22]. Recently, reports have also demonstrated that BBR, in combination with some antifungal agents, can be effective against Staphylococcus aureus, Staphylococcus epidermidis, C. tropicalis, and C. glabrata [23,24,25,26]. In this study, we found synergy between BBR and amphotericin B against candidiasis in vivo. To decrease the drug resistance and increase the sensitivity of C. albicans species, the ideal drug model was a combination of BBR and FLC, with FLC being one of the azole drugs used against candidiasis as a first-line treatment. Accordingly, the evaluation of such a drug model in vitro would contribute to the development of fungi-related therapeutics and targets.
In this study, we report that the molecular mechanism of BBR enhanced the susceptibility of C. albicans to FLC and reduced biofilm formation by altering the expression of efflux pump and adhesion genes, acting as a transcriptional repressor and altering cellular surface hydrophobicity (CSH), which are possible signaling pathways involved in drug resistance. In this study, we showed that the addition of BBR increased C. albicans susceptibility to FLC.
Methods
Strain cultivation and identification
The SC5314 reference strain was used as the control; this strain was procured from the American Type Culture Collection. Thirty-six C. albicans strains that were isolated from the medical laboratory center of Ning Xia Medical University Hospital between January 2021 and June 2022 were used in this study and were inspected by an automated microbiological VITEK-2 COMPACT system. C. albicans species were collected from patient intubating catheters, stored at -80 ℃, and generally maintained in YPD medium (peptone 2%, dextrose 2%, and yeast extract 1%) at 30℃ in a rotatory shaker for 24 h [27]. The cells were washed three times in sterile PBS after centrifugation (2,000 × g, 15 min) and resuspended in RPMI 1640 at the required concentration before the experiment [28].
Agents
Dimethyl sulfoxide (DMSO) solutions, including DMSO solutions of FLC and BBR (National Institutes Drug Control, Hang Zhou, China), were filtered to remove impurities, resulting in an initial concentration of 10,240 µg/ml. Sterile distilled water was then used to dilute the drugs to the desired concentration. The molecular structure of berberine is shown in Fig. 1.
Candida albicans biofilm formation assay
Biofilm formation assays of C. albicans species were performed using 96-well microtiter polystyrene plates as stated previously [29]. Overnight cultures of C. albicans species underwent inoculation followed by shaking in YPD at 37 ℃ for 24 h. Subsequently, the bacteria were washed and resuspended in RPMI 1640 at a concentration of 1 × 106 cells/ml for biofilm formation. Briefly, the experiment was implemented in 96-well microtiter polystyrene plates with 100 µL cell suspensions; the cells were subjected to a consecutive inoculation in a 5% CO2 atmosphere at 37 °C for 36 h or longer until mature biofilm formation was achieved without shaking. The nonadherent cells were washed in sterile PBS three times, then 200 µL methanol per well was added and the plates were left to dry for 15 min. The well plates then underwent crystal violet staining, with additional stain flushed under running water, and the biomass was quantified based on crystal violet OD570 nm absorbance.
Susceptibility testing
The reference M44-A method is described by the Clinical and Laboratory Standards Institute (CLSI) [30,31,32]. Mueller Hinton agar (MH-MB, 90 mm in diameter) was complemented with 2% glucose and 0.5 mg/ml methylene blue. Sterile cotton swabs were utilized to adjust the turbidity of organisms in sterile physiological saline to a 0.5 McFarlane standard and the organisms were inoculated on the agar surface. The agar plate surface was covered with 25 µg FLC antifungal drug, and the plates were examined after 24 h of incubation at 37 ℃. The zones of inhibition were detected based on the CLSI criteria for in vitro testing. The results were described as zone diameters; ≤ 14 mm was considered resistant, ≥ 19 mm was considered susceptible, and 15 ~ 18 mm was considered close dependent SDD. We used the SC5314 reference strain as the control strain.
Checkerboard assay
The minimal inhibitory concentrations of C. albicans species were determined using the Clinical and Laboratory Standards Institute guidelines for M27-A2 via the broth microdilution method, as previously reported [33, 34]. Briefly, the fungal liquid was adjusted to 0.5 CFU/ml using a turbidimeter, the concentration was diluted 20 times and then 50 times, resulting in a final concentration of 2 × 103 CFU/ml; an initial concentration of cell suspension (2 × 103 CFU/ml, 100 µl) was added to 96-well microtiter plates. The test concentrations of BBR and FLC were in the range of 1024 to 0.25 µg/ml. After aerobic incubation for 24 h at 37 ℃, the cells were observed. The minimum inhibitory concentration of MIC60 was determined based on the absorbance at 630 nm. The interaction between BBR and FLC was evaluated using the checkerboard method, and the results are represented by the inhibitory concentration index (FICI), with FICI = (MICFLC in combination/MICFLC alone) + (MICBBR in combination/MICBBR alone); synergism was FICI ≤ 0.5, indifference was 0.5 ≤ FICI ≤ 4.0, and antagonism was FICI> 4.0.
Time-kill curve test
In accordance with previous research [35, 36], a single larger colony of C. albicans cells in the exponential growth phase was picked from the medium. After washing with PBS three times, a standard turbidimeter was used to adjust the cell suspension to 15 × 104 CFU/ml, which was then immediately diluted with RPMI 1640 to 5 × 104 CFU/ml. The concentrations tested were 16 µg/ml FLC and 4 µg/ml BBR. The four experimental groups were the control group with no treatment, the FLC alone group, the BBR alone group, and the FLC + BBR combination treatment group. Samples from the experimental groups of C. albicans species were collected at predetermined time points (0, 4, 6, 8, 12, 24, 48 h), with aerobic agitation at 37 ℃ prior to counting. One milliliter samples of each solution were diluted 10 times, spotted on 1 ml YPD medium, and incubated at 37 ℃ for 24 h in a humid 5% CO2 atmosphere before colony counting. The experiment was repeated three times, and the CFU/ml was recorded. The combination treatment increased the killing by ≥ 2 log10 CFU/m compared to either drug alone, which was defined as synergy, while a killing effect within 2 log10 CFU/ml was defined as indifference.
Fluorescence microscopy analysis
Briefly, a stock solution of fluorescein diacetate (FDA, Sigma) was prepared according to previously described methods [37, 38]. A concentration of 10 mg/ml FDA was dissolved in acetone and stock solution at -20 ℃ prior to use. A 1:50 working solution was freshly diluted with sterile PBS before each assay. Subsequently, the concentration of the fungal broth was adjusted to 1 × 106 CFU/ml and incubated with 4 µg/ml BBR plus 16 µg/ml FLC for 4 h, 18 h, and 36 h at 37 ℃. In 96-well microtiter polystyrene plates, 100 µl of fungal liquid was added to 100 µl of FDA working solution. The microtiter polystyrene plates were placed on a shaking table for 30 min at 37 ℃ in the dark before being washed with unbound dye by sterile PBS and scanned under a fluorescence microscope at an emission wavelength of 518 nm and excitation wavelength of 494 nm.
Cell surface hydrophobicity assay
The CSH of C. albicans species was tested using a water‒hydrocarbon two-phase assay [39]. Briefly, C. albicans species were collected in EP tubes, centrifuged at 3000 rpm for 10 min and then washed twice with PBS. Subsequently, the washed solution was removed to obtain 1 × 106 cells/ml (OD600 = 1.0 in YPD medium). Then, 1.2 ml of cell suspension was moved into a clean glass tube and overlaid with 0.3 ml of octane. The cell suspension mixture was gently vortexed for 3 min and placed at 37 °C for 5-10 min for phase separation measurement at OD600 nm [40]. Each group was analyzed three times. The hydrophobicity was calculated using the following formula: [OD600 of control group - OD600 of test group)/OD600 of control group] × 100%.
qRT‒PCR gene expression analysis
In a 24-well cell culture plate, 10% RPMI 1640 cell culture medium was added and the plates were incubated at 37 ℃ overnight and washed three times with sterile PBS solution. The concentration of the broth was adjusted to 2 × 106 CFU/ml using a Mike turbidimeter, and 1 ml of the fungal suspension was added to each well of the 24-well plates; the plates were then placed in a constant-temperature shaker at 75 r/min and incubated at 37 ℃ for 2 h. The adhesion period of the yeast was followed by gently aspirating the suspension with a pipette and gently washing the wells two to three times with PBS. Then, 1 ml of RPMI 1640 cell culture medium was added and the plates were incubated at 37 ℃ 75 r/min for 24 h. During the maturation of the biofilms of the C. albicans species, the supernatant was gently removed with a pipette and the cells were washed twice to three times with PBS. After 36 h, the samples from four groups were collected: group A with 32 µg/ml FLC, group B with 8 µg/ml BBR, group C with 32 µg/ml FLC and 8 µg/ml BBR, and group D as the control group. For each group, RNA was extracted using an RNA purification kit (TaKaRa Ra Biotechnology). All the experimental steps took place on ice to prevent RNA or DNA degradation. First, cDNA was obtained using a reverse transcription kit (Takara Biotechnology). The thermal cycling conditions were as follows: an initial step at 94 ℃ for 5 min followed by 38 cycles and maintenance at 95 ℃ for 10 s, 58 ℃ for 30 s, and 72 ℃ for 30 s. After the reaction, the Light Cycler Real-Time PCR system was used to measure the CT of different genes. An endogenous reference control, 18 S RNA, was used. The primers of the ALS3, CDR1, MDR, TUP1, HWP1, and 18 S sequences are listed in Table 1. The gene levels were calculated using the following formula: 2 - (△△CT) [41]. The independent experiment was replicated three times.
Bioinformatics analysis
The biological significance of C. albicans species was investigated, and KEGG pathway enrichment analysis was used to analyze the pathways and functional annotations of different proteins. KEGG pathway and gene ontology (GO) analyses are common methods for functional systematic analyses of genes. These analyses were implemented for different proteins, with a p value less than 0.05 deemed significant for every comparison. The KEGG enrichment pathway analysis of the differently expressed genes between the groups and differentially methylated sites resulted in the identification of related genes, including ERG1, ERG3, ERG4, ERG5, ERG6, ERG9, and ERG24. The ERG25 gene participates in C. albicans signaling pathways. The utilized primer sequences are listed in Table 2 [42,43,44].
Statistical analysis
All experiments were performed in triplicate on independent occasions. The results are displayed as the mean ± standard deviation, and p < 0.05 was regarded as statistically significant. The statistical analyses were performed using SPSS 22.0.
Results
Evaluation of the biofilm-producing capacity of clinical Candida albicans species
To assess the biofilm formation capacity of the clinical species, 24 species of Candida albicans analyzed via crystal violet staining. The results showed that the C. albicans species had distinct biofilm-producing capacities, which could be categorized as high biofilm formation (HBF), intermediate biofilm formation (IBF), or low biofilm formation (LBF) based on the absorbance at OD570nm. Eight clinical species fell into the LBF group with absorbance values (Q1 OD570nm = 0.25) less than the first quartile, 10 clinical strains were in the HBF group with absorbance values greater than the second quartile (Q2 OD570nm = 0.5), and six strains were categorized in the IBF group with OD570nm values between the first and second quartiles Q1–Q2 (Fig. 2d). Based on macroscopic observation, the HBF strains had intact filamentous cells and a small number of yeast cells (Fig. 2c), while the LBF strains showed a lack of hyphal cells and a predominance of yeast cells (Fig. 2a). The species with capacities between the HBF and LBF strains were defined as IBF species, which had moderate filamentous phenotypes (Fig. 2b).
In vitro antimicrobial susceptibility test
The sensitivity of C. albicans to fluconazole was determined by the Kirby-Bauer(K-B) method. The capacity of FLC to inhibit strain growth was determined by measuring the size of the diameter. The results showed that the SC5314 reference strain and twenty C. albicans clinical strains were sensitive to FLC, with only two strains resistant to FLC. Among the library of pathogenic microorganisms, C. albicans had very little resistance to FLC; thus, the study of this fungus was considered to be of great significance. The susceptibility test results revealed that the diameter of the inhibition zone of the SC5314 reference strains was 23 mm (Fig. 3a), and the diameter of the inhibition zone of the C. albicans species was 13 mm. Therefore, C. albicans species were selected as experimental strains (Fig. 3b).
In vitro synergism of drugs against FLC-resistant C. Albicans
BBR is a natural isoquinoline alkaloid that is derived from coptis and phellodendron [45]. In this report, we evaluated the effect of the combination of FLC and BBR on the SC5314 reference strain and C. albicans species via the checkerboard method, and the results are summarized in Table 3. When combining BBR and FLC treatment for the SC5314 reference strain, FLC only needed to be given at 4 µg/ml to achieve the MIC60 antifungal rate, and if FLC is used alone, 32 µg/ml was required. For C. albicans species, to reach the MIC60 antifungal rate, fluconazole only needed to be given at 32 µg/ml when used in combination, which is significantly lower than 512 µg/ml required when used alone, and the FICI < 0.5, which means that the two agents were synergistic (Table 3). Furthermore, FLC inhibition of C. albicans was dose dependent; this dose dependent activity was confirmed by a concentration growth curve assay. As shown in Fig. 4, the bacteriostatic effect of FLC alone was poor, and the rate of inhibition was still less than 40% when the concentration of FLC reached 64 µg/ml, but the addition of BBR improved the inhibition of C. albicans growth. Meanwhile, the dose of FLC and BBR used was decreased in the combination.
Time-kill curves
To further investigate the dynamics of antifungal agents, we used time-kill curves, as they provided kinetic information on C. albicans growth and detailed information on the antifungal effects of the agents. The combination of FLC and BBR had a synergistic inhibitory effect on the growth of C. albicans species, and the growth was significantly different than that of those treated with single agents (p < 0.05) (Figs. 5 and 6). Time-kill curve analysis was performed according to a previous protocol and revealed the variation in C. albicans susceptibility to the different agents (Fig. 7). We found that 4 µg/ml BBR had a poor antifungal effect, while the addition of 16 µg/ml FLC had a stronger effect than BBR alone. The combination of the two agents led to a 2.1 log10 CFU/ml decrease in bacteria compared to that of FLC alone.
The sensitivity of C. Albicans to FLC increased after BBR treatment
We used a fluorescence microscope to observe the morphological changes in C. albicans species after exposure to the combination of BBR and FLC (16 µg/ml FLC and 4 µg/ml BBR). BBR has been shown to enhance the susceptibility of C. albicans species to FLC. The addition of 4 µg/ml BBR with FLC dramatically increased the susceptibility of biofilm-producing C. albicans species to FLC. Morphological analysis using fluorescence microscopy further showed that a combination of 16 µg/ml FLC and 4 µg/ml BBR inhibited C. albicans growth and hyphal formation (Fig. 8).
BBR and FLC synergistically decreased the hydrophobicity of C. Albicans
The effect of the different treatments on cell hydrophobicity was determined based on the principle of water‒hydrocarbon two-phase separation. The measurement results showed that when FLC alone was added at a concentration of 32 µg/ml or BBR alone was added at 8 µg/ml, the CSH was reduced when compared with that of the untreated group (p < 0.05). When the agents were used in combination, the decrease in CSH decrease was obvious and the was significant when compared with that of FLC or BBR alone (p < 0.05), as shown in Fig. 9. The results revealed that the combination of FLC and BBR synergistically decreased the cell surface hydrophobicity of C. albicans species.
The combination of FLC and BBR reduced the expression of drug resistance genes
To investigate the mechanism of the effects of the FLC and BBR combinations on C. albicans species, we further investigated the underlying molecular mechanism involved in increasing the sensitivity of biofilm-producing C. albicans species to combination treatment; the overexpression of genes is one of the usual mechanisms of drug resistance in C. albicans species, such as those encoding drug targets such as ERG11, the efflux pump genes CDR1 and MDR, the adhesion function gene ALS3, and the hypha-related gene HWP1. The transcription of the TUP1 virulence genes in the cells under 32 µg/ml FLC and 8 µg/ml BBR treatment was measured by qRT‒PCR assay. Consistent with previous reports, our results showed that FLC and BBR combination treatment, which was remarkably different when compared with treatment with the same agents alone (p < 0.05), inhibited the expression of CDR1, MDR, ALS3, and HWP1. Of note, both FLC, BBR, and the synergistic combination treatment failed to inhibit TUP1 gene expression in the clinical C. albicans species (Fig. 10).
DNA extraction, genome sequencing, assembly, and annotation
The genomes of the C. albicans species were sequenced using the Illumina HiSeq2000 platform. The final assembly included 15,445,198 bp and consisted of 72 scaffolds and 72 contigs. The overall GC content of the genome was 35.13% (Table 4). The data generated with the PACBIO platform were corrected to the generated assembly sequence by the Illumina platform to alleviate the sequencing error rate caused by the relatively high pressure. The average gene length was 1435 bp. The number of tRNA, SnRNA genes, and microsatellite sequences (144, 32, and 2384, respectively) identified in the assembly were > 2X that of other Candida sp.
In total, 6276 protein-coding genes were predicted with a total length of 9,006,661 bp in the genomes of the C. albicans species. Additionally, the VFDB database predicted that the C. albicans species had 763 virulence-related genes involved in the formation and regulation of toxins and invasive enzymes. The CARD database annotation showed 82 related resistance genes were associated with C. albicans. GO terms were assigned to the 14,972 protein-coding genes of the C. albicans species. Almost 14,972 genes were assigned to each of the following categories: molecular function, cellular components, and biological process (Fig. 11). A total of 8149 genes participating in 40 pathways were annotated using the KEGG database (Fig. 12). Six of these pathways were significantly enriched, including cellular processes, environmental information processing, genetic information processing, human diseases, and metabolism in organismal systems. KEGG analysis was implemented to categorize the functions of the identified genes, with the amino acid metabolic pathway revealed to be the most significantly enriched pathway; the metabolic genes of the C. albicans species were mainly involved in acid and carbon metabolism. These results indicate that these two types of genes account for a significant proportion of the biological processes involved in Candida growth, death, signal transduction, and translation. This result may explain how C. albicans acquires drug resistance as a conditional pathogen.
The phylogenetic tree of the 13 C. albicans species was constructed using homologous protein analysis (Fig. 13). In the phylogenetic tree, nodes represent the closest common ancestor of each branched strain, and the length of the line segment between the nodes represents the evolutionary distance of each strain. As seen from the figure, in this study, we showed that the C. albicans species were closely related to P34048, which is from the same family, but the SC5314, Ca6, and other strains were distant. The collinearity analysis showed that the collinearity of the nucleic acids between the C. albicans species and SC5314 was distinctive. Phylogenetic analysis and the GO database were used to further analyze the gene information to extract effective information from the large amount of gene data; this information could provide valuable clues for future research.
Gene duplication and collinearity analysis
The interspecific synteny among the WO-1, SC5314, and P34048 strains and the C. albicans species were analyzed and compared to further explore the evolution of the resistance genes. Interestingly, the C. albicans species and P3408 showed a highly similar collinear relationship (Fig. 14). However, the gene correspondence between the C. albicans species and the SC5314 and WO-1 strains was poor. These results are consistent with what was seen in the phylogenetic tree. The figures show the structural variation that occurs in the process of evolution, such as the insertion and deletion of sequences and the genetic changes caused by transposition, recombination, and other mechanisms.
Expression of genes in the steroid biosynthesis pathway
Figure 15 shows the steroid biosynthesis pathway, where genes marked in red had upregulated expression. Ergosterol is a key target gene in cell membrane synthesis and a therapeutic target for antifungal drugs. Most antifungal drugs counter the synthesis of steroids by inhibiting the activity of genes in the ERG family. Compared with the control group, ERG1, ERG3, ERG4, ERG5, ERG24, and ERG25 expression in the FLC alone or BBR alone group was upregulated; however, in the FLC and BBR groups, there was significantly downregulated expression of these genes compared with the other groups (P < 0.05) (Fig. 16). No significant difference was detected between the FLC or BBR alone treatment groups. ERG6 and ERG9 expression was significantly increased in the FLC and BBR groups compared to the FLC alone or BBR alone groups (P < 0.05). The figure shows that the combination of BBR and FLC increased drug sensitivity by regulating the expression of genes involved in intracellular mechanisms in Candida albicans.
Discussion
Candida albicans is an opportunistic pathogen that often infects immunocompromised patients [5]. Although FLC agents are used to prevent Candida albicans infections, the emergence of drug resistance remains a clinical challenge. With the increasing resistance of Candida albicans to FLC, an increasing number of experiments have focused on the development of new antifungal drugs. In this regard, we investigated possible risk factors, including biofilm formation. While the molecular mechanism of Candida albicans resistance to FLC has been extensively researched in the past few years, to the best of our knowledge, there is no effective treatment for Candida albicans infections in clinical settings [46, 47]. An increasing amount of evidence has suggested that Candida albicans biofilms could accelerate resistance to antifungal agents. Accordingly, studies with the aim of identifying potential antifungal drugs or targets to prevent further drug-resistance are urgently needed.
In recent years, several studies have reported that a combination of BBR and antibiotics could synergistically increase the susceptibility of Candida albicans to the antifungal FLC [26]. Previous studies have shown that BBR can induce DNA damage [48]. In this study, our results demonstrated that BBR can increase the sensitivity of Candida albicans by inhibiting biofilm and hyphal formation, drug transportation, and cell adhesion while simultaneously decreasing hydrophobicity and increasing the concentration of drugs entering the cell.
The biofilm of Candida albicans is composed of a mixture of yeast cells and mycelium, and the formation of biofilms is important for drug resistance [49, 50]. In our research, BBR and FLC treatment was accompanied by downregulation of the expression of the CDR1, MDR, ALS3, ERG11, and HWP1 genes. Among them, CDR1 and MDR are drug resistance genes that encode ATP-binding transporters. Researchers have shown that one of the main reasons behind FLC resistance in Candida albicans is the overexpression of efflux pump genes [51, 52]. Consistent with this discovery, we found downregulated expression of the CDR1 and MDR genes after BBR and FLC combination treatment, indicating that BBR increased the susceptibility of Candida albicans to FLC in part by inhibiting drug transporter gene expression. In contrast, Marchetti et al. [53] indicated that the efflux transporter genes CDR1, CDR2, and MDR had independent synergistic antifungal effects with FLC against Candida albicans. ALS is one of the main gene adhesion families that controls biofilm production. In particular, the ALS3 gene is related to biofilm formation, with its expression gradually increasing after biofilm formation [54, 55]. In our experiment, the same conclusion was obtained. Compared with FLC alone, BBR combined with FLC inhibited ALS3 expression, which disrupted biofilm production. These results demonstrated upregulated expression of TUP1 mRNA with filamentous growth and FLC treatment alone. TUP1 was associated with a transcriptional repressor in the Candida albicans-resistant strain. Increased expression of the TUP1 gene could inhibit the activity of hyphae after FLC and BBR combination treatment, thereby increasing drug sensitivity. HWP1 gene production is related to the ability of cell wall proteins to stimulate a host immune response, which then causes changes at the molecular level [56]. In a previous report, the overexpression of HWP1 significantly increased ALS3 expression and reduced hyphal development. Consistent with this research, our results suggested that FLC had weak antifungal effects on Candida albicans-resistant species. Therefore, the FLC and BBR combination will be effective in reducing side effects and protecting against biofilm formation.
Ergosterol is a part of the cell membrane and has a variety of cellular functions, such as maintaining membrane integrity and membrane enzyme binding [57]. One crucial mechanism of fungal resistance to drugs is changes in the transcription of genes related to ergosterol biosynthesis. KEGG metabolic pathway analysis following whole-genome sequencing verified the previously described results regarding the effects of drugs on ergosterol biosynthesis; the sterol genes were found to be related to FLC resistance. Based on the metabolic pathway analysis, upregulated genes and other ERG family genes were selected for experimental testing. Of note, FLC and BBR combined treatment decreased the gene expression of ERG1, ERG3, ERG4, ERG5, ERG24, and ERG25 and had a remarkably differential effect compared to single treatments. Interestingly, ERG6 and ERG9 expression was significantly increased compared to that in the other groups, indicating that ergosterol is closely related to the drug resistance mechanism. ERG6 is a rate-limiting enzyme in the ergosterol biosynthesis pathway, while ERG9 is a synthase of ergosterol and participates in feedback regulation. When combined with FLC, BBR inhibited the expression of other sterol genes, causing an increase in ERG6 and ERG9 gene expression via a feedback mechanism and accelerating the synthesis of ergosterol. In conclusion, BBR can enhance sensitivity to FLC by reducing the expression of mycelial adhesion genes, efflux-related proteins and ergosterol genes, thereby effectively reducing the hydrophobicity of the Candida albicans cell surface to increase sensitivity. Whole-genome sequencing can be used to identify as many resistance mechanisms as possible.
Data availability
All the data supporting our findings are contained within the manuscript.
References
Ho J, Camilli G, Griffiths JS, Richardson JP, Kichik N, Naglik JR. Candida albicans and candidalysin in inflammatory disorders and cancer. Immunology. 2021;162:11–6.
Duggan S, Leonhardt I, Hünniger K, Kurzai O. Host response to Candida albicans bloodstream infection and sepsis. Virulence. 2015;6:316–26.
Hoenigl M, Sprute R, Arastehfar A, Perfect JR, Lass-Flörl C, Bellmann R, et al. Invasive candidiasis: investigational drugs in the clinical development pipeline and mechanisms of action. Expert Opin Investig Drugs. 2022;31:795–812.
Chen H, Zhou X, Ren B, Cheng L. The regulation of hyphae growth in Candida albicans. Virulence. 2020;11:337–48.
Pereira R, Fontenelle ROS, Brito EHS, Morais SM. Biofilm of Candida albicans: formation, regulation and resistance. J Appl Microbiol. 2021;131:11–22.
Wall G, Montelongo-Jauregui D, Bonifacio BV, Lopez-Ribot JL, Uppuluri P. Candida albicans biofilm growth and dispersal: contributions to pathogenesis. Curr Opin Microbiol. 2019;52:1–6.
Huang X, Zheng M, Yi Y, Patel A, Song Z, Li Y. Inhibition of berberine hydrochloride on Candida albicans biofilm formation. Biotechnol Lett. 2020;42:2263–9.
Holmes AR, Tsao S, Ong SW, Lamping E, Niimi K, Monk BC, et al. Heterozygosity and functional allelic variation in the Candida albicans efflux pump genes CDR1 and CDR2. Mol Microbiol. 2006;62:170–86.
El-Ganiny AM, Kamel HA, Yossef NE, Mansour B, El-Baz AM. Repurposing pantoprazole and haloperidol as efflux pump inhibitors in azole resistant clinical Candida albicans and non-albicans isolates. Saudi Pharm J. 2022;30:245–55.
Castagnola E, Jacqz-Aigrain E, Kaguelidou F, Maragliano R, Stronati M, Rizzollo S, et al. Fluconazole use and safety in the nursery. Early Hum Dev. 2012;88(Suppl 2):S11–5.
Wang WY, Cai HQ, Qu SY, Lin WH, Liang CC, Liu H, et al. Genomic variation-mediating fluconazole resistance in yeast. Biomolecules. 2022;12:845.
Jacobs SE, Jacobs JL, Dennis EK, Taimur S, Rana M, Patel D, et al. Candida Auris pan-drug-resistant to four classes of antifungal agents. Antimicrob Agents Chemother. 2022;66:e0005322.
Dovo EE, Zohoncon TM, Tovo SF, Soubeiga ST, Kiendrebeogo IT, Yonli AT, et al. First detection of mutated ERG11 gene in vulvovaginal Candida albicans isolates at Ouagadougou/Burkina Faso. BMC Infect Dis. 2022;22:678.
Esfahani A, Omran AN, Salehi Z, Shams-Ghahfarokhi M, Ghane M, Eybpoosh S, et al. Molecular epidemiology, antifungal susceptibility, and ERG11 gene mutation of Candida species isolated from vulvovaginal candidiasis: comparison between recurrent and non-recurrent infections. Microb Pathog. 2022;170:105696.
Navarro-RodrÃguez P, López-Fernández L, Martin-Vicente A, Guarro J, Capilla J. ERG11 polymorphism in voriconazole-resistant Candida tropicalis: weak role of ERG11 expression, ergosterol content, and membrane permeability. Antimicrob Agents Chemother. 2020;65:e00325–20.
Liu J-Y, Wei B, Wang Y, Shi C, Li W-J, Zhao Y, et al. The H741D mutation in Tac1p contributes to the upregulation of CDR1 and CDR2 expression in Candida albicans. Braz J Microbiol. 2020;51:1553–61.
Nademi A, Shams-Ghahfarokhi M, Razzaghi-Abyaneh M. Effect of Allium cepa loaded polyacrylonitrile and polyvinyl pyrrolidone nanofibers on Candida albicans growth and the expression of CDR1 and CDR2 genes. Curr Med Mycol. 2021;7:28–33.
Nikmanesh B, Ahmadikia K, Getso MI, Gharehbolagh SA, Aboutalebian S, Mirhendi H, et al. Candida Africana and Candida Dubliniensis as causes of pediatric candiduria: a study using HWP1 gene size polymorphism. AIMS Microbiol. 2020;6:272–9.
Song D, Hao J, Fan D. Biological properties and clinical applications of berberine. Front Med. 2020;14:564–82.
Feng X, Sureda A, Jafari S, Memariani Z, Tewari D, Annunziata G, et al. Berberine in cardiovascular and metabolic diseases: from mechanisms to therapeutics. Theranostics. 2019;9:1923–51.
Zhu C, Li K, Peng X-X, Yao T-J, Wang Z-Y, Hu P, et al. Berberine a traditional Chinese drug repurposing: its actions in inflammation-associated ulcerative colitis and cancer therapy. Front Immunol. 2022;13:1083788.
Qin X, Jiang M, Zhao Y, Gong J, Su H, Yuan F, et al. Berberine protects against diabetic kidney disease via promoting PGC-1α-regulated mitochondrial energy homeostasis. Br J Pharmacol. 2020;177:3646–61.
Asila A, Liu J, Liu J, Li L, Liao J. Immunomodulatory effects of berberine on Staphylococcus aureus-induced septic arthritis through down-regulation of Th17 and Treg signaling pathways. Acta Biochim Pol. 2022;69:215–26.
Wang X, Qiu S, Yao X, Tang T, Dai K, Zhu Z. Berberine inhibits Staphylococcus Epidermidis adhesion and biofilm formation on the surface of titanium alloy. J Orthop Res. 2009;27:1487–92.
Shi G, Shao J, Wang T, Wu D, Wang C. Mechanism of berberine-mediated fluconazole-susceptibility enhancement in clinical fluconazole-resistant Candida tropicalis isolates. Biomed Pharmacother. 2017;93:709–12.
Xie Y, Liu X, Zhou P. In vitro antifungal effects of berberine against Candida spp. in planktonic and biofilm conditions. Drug Des Devel Ther. 2020;14:87–101.
Boone CHT, Gutzmann DJ, Kramer JJ, Atkin AL, Nickerson KW. Quantitative assay for farnesol and the aromatic fusel alcohols from the fungus Candida albicans. Appl Microbiol Biotechnol. 2022;106:6759–73.
Mat-Rani S, Chotprasert N, Srimaneekarn N, Choonharuangdej S. Fungicidal effect of lemongrass essential oil on Candida albicans biofilm pre-established on maxillofacial silicone specimens. J Int Soc Prev Community Dent. 2021;11:525–30.
Pierce CG, Uppuluri P, Tummala S, Lopez-Ribot JL. A 96 well microtiter plate-based method for monitoring formation and antifungal susceptibility testing of Candida albicans biofilms. J Vis Exp. 2010;21:2287.
de la Rosa-GarcÃa E, Olalde-Hernández MJ, Irigoyen-Camacho ME, Mondragón-Padilla A, Mendoza-Juache A, Sánchez-Vargas LO. Antifungal susceptibility of oral isolates of Candida species from chronic kidney disease patients on chronic dialysis. J Mycol Med. 2020;30:101009.
Tasneem U, Siddiqui MT, Faryal R, Shah AA. Prevalence and antifungal susceptibility of Candida species in a tertiary care hospital in Islamabad, Pakistan. J Pak Med Assoc. 2017;67:986–91.
Aydin F, Bayramoglu G, Guler NC, Kaklikkaya N, Tosun I. Bloodstream yeast infections in a university hospital in Northeast Turkey: a 4-year survey. Med Mycol. 2011;49:316–9.
Sharifzadeh A, Shokri H, Abbaszadeh S. Interaction of carvacroland voriconazole against drug – resistant Candida strains isolated from patients with candidiasis. J Mycol Med. 2019;29:44–8.
Goel N, Rohilla H, Singh G, Punia P. Antifungal activity of cinnamon oil and olive oil against Candida Spp. isolated from blood stream infections. J Clin Diagn Res. 2016;10:DC09–11.
Mijatovic S, Stankovic JA, Calovski IC, Dubljanin E, Pljevljakusic D, Bigovic D, et al. Antifungal activity of Lavandula angustifolia essential oil against Candida albicans: time-kill study on pediatric sputum isolates. Molecules. 2022;27:6300.
Partha ADSL, Widodo ADW, Endraswari PD. Evaluation of fluconazole, itraconazole, and voriconazole activity on Candida albicans: a case control study. Ann Med Surg (Lond). 2022;84:104882.
Peeters E, Nelis HJ, Coenye T. Comparison of multiple methods for quantification of microbial biofilms grown in microtiter plates. J Microbiol Methods. 2008;72:157–65.
Honraet K, Goetghebeur E, Nelis HJ. Comparison of three assays for the quantification of Candida biomass in suspension and CDC reactor grown biofilms. J Microbiol Methods. 2005;63:287–95.
Suchodolski J, Muraszko J, Korba A, Bernat P, Krasowska A. Lipid composition and cell surface hydrophobicity of Candida albicans influence the efficacy of fluconazole-gentamicin treatment. Yeast. 2020;37:117–29.
Haghighi F, Andriasian L, Tran NC, Lux R. Effect of cigarette and e-cigarette smoke condensates on Candida albicans biofilm formation and gene expression. Int J Environ Res Public Health. 2022;19:4626.
Lu L, Li Z, Shan C, Ma S, Nie W, Wang H, et al. Whole transcriptome analysis of schinifoline treatment in Caenorhabditis elegans infected with Candida albicans. Mol Immunol. 2021;135:312–9.
Silva MC, Carbone DCB, Diniz PF, Fernandes FF, Fuzo CA, Trindade CSPC, et al. Modulation of ERG genes expression in clinical isolates of Candida tropicalis susceptible and resistant to fluconazole and itraconazole. Mycopathologia. 2020;185:675–84.
Feng W, Yang J, Xi Z, Qiao Z, Lv Y, Wang Y, et al. Mutations and/or overexpressions of ERG4 and ERG11 genes in clinical azoles-resistant isolates of Candida albicans. Microb Drug Resist. 2017;23:563–70.
Yu L-H, Wei X, Ma M, Chen X-J, Xu S-B. Possible inhibitory molecular mechanism of farnesol on the development of fluconazole resistance in Candida albicans biofilm. Antimicrob Agents Chemother. 2012;56:770–5.
Almatroodi SA, Alsahli MA, Rahmani AH. Berberine: an important emphasis on its anticancer effects through modulation of various cell signaling pathways. Molecules. 2022;27:5889.
Fattouh N, Hdayed D, Geukgeuzian G, Tokajian S, Khalaf RA. Molecular mechanism of fluconazole resistance and pathogenicity attributes of Lebanese Candida albicans hospital isolates. Fungal Genet Biol. 2021;153:103575.
Lin W, Yuan D, Deng Z, Niu B, Chen Q. The cellular and molecular mechanism of glutaraldehyde-didecyldimethylammonium bromide as a disinfectant against Candida albicans. J Appl Microbiol. 2019;126:102–12.
Zhu Y, Ma N, Li H-X, Tian L, Ba Y-F, Hao B. Berberine induces apoptosis and DNA damage in MG–63 human osteosarcoma cells. Mol Med Rep. 2014;10:1734–8.
Li L, Zhang T, Xu J, Wu J, Wang Y, Qiu X, et al. The synergism of the small molecule ENOblock and fluconazole against fluconazole-resistant Candida albicans. Front Microbiol. 2019;10:2071.
Alves F, Mima EGDO, Passador RCP, Bagnato VS, Jorge JH, Pavarina AC. Virulence factors of fluconazole-susceptible and fluconazole-resistant Candida albicans after antimicrobial photodynamic therapy. Lasers Med Sci. 2017;32:815–26.
Maheronnaghsh M, Teimoori A, Dehghan P, Fatahinia M. The evaluation of the overexpression of the ERG-11, MDR-1, CDR-1, and CDR-2 genes in fluconazole-resistant Candida albicans isolated from Ahvazian cancer patients with oral candidiasis. J Clin Lab Anal. 2022;36:e24208.
Gil F, Laiolo J, Bayona-Pacheco B, Cannon RD, Ferreira-Pereira A, Carpinella MC. Extracts from Argentinian native plants reverse fluconazole resistance in Candida species by inhibiting the efflux transporters Mdr1 and Cdr1. BMC Complement Med Ther. 2022;22:264.
Marchetti O, Moreillon P, Entenza JM, Vouillamoz J, Glauser MP, Bille J, et al. Fungicidal synergism of fluconazole and cyclosporine in Candida albicans is not dependent on multidrug efflux transporters encoded by the CDR1, CDR2, CaMDR1, and FLU1 genes. Antimicrob Agents Chemother. 2003;47:1565–70.
Singh S, Nabeela S, Barbarino A, Ibrahim AS, Uppuluri P. Antibodies targeting Candida albicans Als3 and Hyr1 antigens protect neonatal mice from candidiasis. Front Immunol. 2022;13:925821.
Bruder-Nascimento A, Camargo CH, Mondelli AL, Sugizaki MF, Sadatsune T, Bagagli E. Candida species biofilm and Candida albicans ALS3 polymorphisms in clinical isolates. Braz J Microbiol. 2015;45:1371–7.
Orsi CF, Borghi E, Colombari B, Neglia RG, Quaglino D, Ardizzoni A et al. Impact of Candida albicans hyphal wall protein 1 (HWP1) genotype on biofilm production and fungal susceptibility to microglial cells. Microb Pathog. 2014;69–70:20–7.
Suchodolski J, Muraszko J, Bernat P, Krasowska A. Lactate like fluconazole reduces ergosterol content in the plasma membrane and synergistically kills Candida albicans. Int J Mol Sci. 2021;22:5219.
Acknowledgements
We thank the Ningxia Key Laboratory of Clinical and Pathogenic Microbiology for Wei Jia. We thank Gang Li, Xia Qiao, Caiyun Li, Le Du, and Yu Lu for their valuable discussion and assistance.
Funding
This work was supported by the Liuzhou Science and Technology Project (2021YB0104B026). The funders had a role in the study design.
Author information
Authors and Affiliations
Contributions
LD and WJ conceived and designed the experiments, SW analyzed the data and drafted the manuscript; YL and HJ performed experiments and acquired data; HC and Shifu Tang collected samples; All authors read and approved the final version of the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Wu, S., Jia, W., Lu, Y. et al. Mechanism and bioinformatics analysis of the effect of berberine-enhanced fluconazole against drug-resistant Candida albicans. BMC Microbiol 24, 196 (2024). https://doi.org/10.1186/s12866-024-03334-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12866-024-03334-0