3.1. ERG11
By 2010, over 160 amino acid substitutions had been identified in the ERG11 gene, each with varying genetic consequences
[68][69][70][87,88,89]. Many single substitutions are synonymous and have no impact on gene function. In addition, non-synonymous single-nucleotide changes occurring in
Candida strains do not inherently contribute to antifungal resistance. For example, White et al. (2002) identified D116E and E266D, the two most frequent ERG11 substitutions in one sample set, using restriction fragment length polymorphism (RFLP) analysis, with no consistent correlation
[31][36]. These instances indicate that there is natural genetic variation within each fungal species and that allelic polymorphisms are relatively common. Mutations present in both resistant and susceptible samples are an indicator that the alteration is not directly implicated in conferring antifungal resistance
[67][86].
Candida spp. have acquired azole resistance via ERG11 point mutations that typically lower azole binding affinity to the Erg11 active site. Additionally, gain-of-function mutations in upstream transcriptional regulators can increase ERG11 expression and confer resistance.
[29][67][34,86]. Point mutations resulting in a defective Erg11 enzyme unable to bind to azoles have been clustered in three hotspot regions: 105–165, 266–287 and 405–488
[71][94].
Polymorphisms can have different biological impacts, depending on whether they are singly present or in combination with other relevant SNPs. Resistance levels can be enhanced when some mutations with moderately low impact alone are present simultaneously with other resistance mutations
[67][72][73][86,103,104].
The variety of currently available azole medications have different structural features (e.g., short and long chains); therefore, mutations affecting their efficacy may differ. For example, ERG11 point mutations K128T and Y132H may affect the ability of fluconazole or voriconazole molecules to enter or bind the target active site. Mutations in other gene sequences can also confer resistance, such as the G464S mutation, which affects haem coordination due to its location near a key cysteine residue
[74][93]. These mutations do not have the same binding inefficiency for posaconazole and itraconazole treatments, suggesting that these two antifungals have other key interaction sites within the Erg11 protein
[74][93].
3.2. Mutations in Transcriptional Regulators
Zn2-Cys6 transcription factor
uptake control 2 (Upc2) regulates most of the genes in the ergosterol biosynthesis pathway on some level. Gene overexpression of UPC2 can be induced upon azole exposure and can sufficiently compensate for the inhibition of target enzymes. Gain-of-function (GOF) mutations in UPC2 can drive this gene overexpression and fluconazole resistance
[75][101]. For example, a series of studies of azole-resistant clinical isolates identified the A643V substitution in the UPC2 gene, which was validated in vitro to confer resistance
[75][76][101,107]. This mutation and possibly others within this gene sequence region (G648D) may cause Upc2 to be released from a repressor, inducing hyperactivity
[76][107]. However, alternate models exist that involve the sterol regulator, SREBP
[76][107].
A genome-wide ChIP (chromatin immunoprecipitation) study used to identify Upc2-bound gene promoters identified up to 202 genes, including UPC2 itself
[77][111]. Other upregulated genes were found to be involved in ergosterol biosynthesis; oxidoreductase activity; and numerous drug efflux pumps, including MDR1 (MFS-transporter) and CDR1 (ABC-transporter)
[75][101].
3.3. Other ERG Genes and Toxic Diol Formation
The inhibition of Erg11 alters pathway products and induces the synthesis of a fungistatic toxic diol (14α-methylergosta 8–24 (28) dienol) by downstream enzymes (Erg3, Erg6, Erg25, Erg26 and Erg27)
[1]. Erg3 is a C5 sterol desaturase enzyme needed for the conversion of episterol to ergostatrienol
[78][112]. When its expression is inhibited by mutation or deletion, the reduction in toxicity due to the inhibition of toxic diol formation is sufficient to confer resistance in some
Candida spp.
[79][80][81][113,114,115]. However, ERG3 inactivation appeared to minimally contribute to azole resistance in a wide range of studied clinical
C. albicans [82][116]. Q139A substitution in ERG3 has been identified from
N. glabrata clinical isolates with azole resistance
[83][99].
Deletion of ERG pathway genes such as the ERG6 gene can impact resistance to other antifungals. This ∆24 sterol C-methyl transferase is non-essential for ergosterol biosynthesis but it is needed for toxic diol formation. Its disruption contributes to azole resistance in
C. albicans [84][85][117,118].
One differential gene expression analysis of a lab-generated
C. albicans strain with resistance to fluconazole and amphotericin B identified numerous upregulated ergosterol pathway genes, including ERG5, ERG6 and ERG25
[78][112]. Additionally, genes involved in cell stress responses were found to be upregulated, including DDR48 and RTA2
[78][112]. In
C. albicans, RTA2 is a key gene in calcium signaling pathways, and it has been shown to modulate azole resistance, including biofilms
[86][120].
4. Cell Membrane Proteins and Antifungal Resistance
Two types of membrane transporters have been implicated in azole resistance: ABC-Ts (ATP-binding cassette transporters) and MFS-Transporters (major facilitator superfamily transporters)
[31][87][88][89][36,110,121,122]. ABC-Ts facilitate the movement of molecules across membranes using energy derived from ATP hydrolysis, while MFS-Ts require a proton gradient across the plasma membrane to transport foreign molecules out of the cell
[1]. Both types of transporters can bind azoles as a substrate, and the drug can be exported out of the cell. This decreases the intracellular drug concentration and allows cells to circumvent the antifungal effects
[1]. An additional MLT1 (ABC-T) transporter has been implicated in
C. albicans resistance. This multidrug resistance protein (MRP) is localized to the vacuolar membrane and can import azole molecules into the vacuole for sequestration. Mutations in the MLT1 sequence can cause incorrect localization or the inability to bind and transport azoles
[90][123].
4.1. Drug Efflux Pump/Transporter Genes and Resistant Mutations
Despite the large number of transporters found in the
C. albicans genome, evidence of transporter overexpression in resistant clinical isolates is currently limited to ABC-Ts CDR1 and CDR2 and MFS-Ts MDR1 and PDR16
[1][77][91][92][1,111,131,133]. CDR1 and CDR2 overexpression is frequently observed in clinical isolates, and coregulation of these two pumps is evident
[31][36]. In fact, prolonged exposure to azoles can result in trisomy development of chromosome 3, which encodes CDR1 and CDR2, as well as subsequent overexpression of these genes
[93][124]. For azole-resistant
N. glabrata clinical isolates, overexpression of ABC-T genes CDR1, PDH1 (CDR2 in
C. albicans) and SNQ2 and MFS-T genes QDR2, FLR1 and PDR16 has been cited
[73][94][95][96][97][98][99][100][62,104,134,135,136,137,138,139].
4.2. Transcriptional Regulators of Transporter Genes
Efflux pump overexpression can be further induced by upstream GOF mutations, regulating transcription factor genes
[47][101][102][103][104][105][106][107][108][59,125,126,127,128,129,130,144,145]. In
C. albicans, the CDR1, CDR2 and PDR16 transporters, as well as MDR1, are regulated by zinc-cluster (Zn2-Cys6) transcription factors Tac1 and Mrr1, respectively
[77][88][91][92][107][109][111,121,131,133,144,146]. Other potential transcriptional regulators of CDR1 expression include Tup1 (thymidine uptake 1) and Ncb2 (β subunit of the NC2 complex)
[110][147]. For MDR1, additional transcriptional regulators are Cap1 (bZIP transcription factor) and Mcm1, but no antifungal-resistant mutations have been identified to date
[111][112][148,149].
4.3. Post-Translational Regulation of Transporter Genes
In addition to transcriptional and translational control of efflux pump genes, there is evidence of post-translational regulation. For instance, mitochondrial biogenesis gene FZO1 is important for directing Cdr1 to the correct membrane
[113][153]. In FZO1 deletion mutants, Cdr1 was found to be mis-sorted to the vacuole, which was correlated with increased azole susceptibility
[113][153].
5. The Cell Wall Biosynthesis Pathway and Antifungal Resistance
5.1. FKS1 and FKS2 Sequence Mutations
Echinocandin resistance is associated with modifications to the FKS1 or FKS2 gene sequences, which code for the β1–3 glucan synthase enzyme
[114][156]. FKS1 mutations have been identified in resistant
C. albicans,
C. tropicalis,
P. kudriavzevii and
N. glabrata, while FKS2 mutations have only been identified in
N. glabrata [115][116][117][118][157,158,159,160]. No definitive intrinsic resistance has been established in any
Candida species, but secondary resistance can be acquired in individual isolates through point mutations. Notably,
C. parapsilosis and
Meyerozyma guilliermondii have a higher rate of spontaneously occurring FKS1 point mutations and may be considered more intrinsically resistant
[58][59][77,78].
Significant GOF mutations in the FKS1 gene sequence in clinical isolates have been described, particularly in the hotspot regions of 637–654 and 1345–1365
[56][115][119][120][121][122][123][124][125][74,157,162,163,164,165,166,167,168]. Walker et al. (2010) noted that non-synonymous substitutions at aa position 645 have been commonly observed. Here, serine substitutions with phenylalanine, proline or tyrosine have been cited
[120][121][126][163,164,169]. Hotspot mutations are usually dominant, and
C. albicans fungal cells only require one mutant allele for resistance to be conferred across the three echinocandins
[56][74].
Evidence of FKS2-resistant mutations were previously limited to in vitro experiments, but recently, such mutations have been identified in clinical
N. glabrata isolates
[29][30][127][34,35,173]. This coincides with evidence suggesting that FKS2 has higher levels of expression in
N. glabrata than FKS1 and that mutations may have a greater influence on echinocandin resistance in this species
[116][158].
5.2. Transcriptional Regulators of fks Genes
Transcriptional regulators upstream of
fks genes can also affect echinocandin resistance. In particular, point mutations in transcription factor PDR1 have been found in numerous resistant
N. glabrata isolates (
Figure 3).
5.3. Protein Analysis Associated with Echinocandin Resistance
Cell wall remodeling enzymes upregulated in previous protein studies include glucanosyl transferases Phr1, Phr2 and Crh, as well as chitin-glucanosyl transferase family proteins
[128][129][178,179]. Proteomic analysis conducted using mass spectrometry (LC-MS/MS) revealed that levels of cell wall organization and maintenance proteins can differ between drug-resistant and susceptible strains in response to caspofungin treatment
[130][131][180,181]. Differentially expressed enzymes related to cell wall synthesis and remodeling include Sun41, Gsc1, Pmt1, Mnt1, Als3, Als4, Ecm33 and Pga31
[131][181]. Validated caspofungin tolerance regulators Cas5, Mkc1, Swi4, Gin4, Stt4, Ahr1 and Pkc1 were detected in an alternate screen
[132][182]. Furthermore, metabolic enzymes with immunogenic activity including Eno1, Fba1, Gpm1 and Pgk1 have also been observed to be released in
Candida cells exposed to caspofungin
[131][133][181,183]. Finally, the Hsp90 molecule may have a regulatory role with key resistance regulators like Mkc1 from the Pkc1 signaling pathway
[134][135][184,185]. These protein subsets may be good
Candidates for diagnostic markers that predict echinocandin resistance in
Candida species. Further validation across a wide range of antifungal-resistant clinical isolates is still needed
[131][181].
6. The Nucleic Acid Biosynthesis Pathway and Antifungal Resistance
The biosynthesis of nucleic acids (DNA and RNA) and subsequent protein synthesis in pathogenic fungi can be targeted with nucleoside analogue 5-fluorocytosine (5-FC). As a prodrug, it requires activation within the fungal cell via metabolism by the pyrimidine salvage pathway [65][83]. Then, it is incorporated as a toxic substrate, and the affected nucleotides have damaging effects on cell viability [65][83]. Membrane permeases encoded by FCY2 (cytosine permease) and other homologs (FCY21 and FCY22) are responsible for the active transport of 5-FC into the cell (Figure 3) [65][83]. 5-FC is then converted to toxic 5-fluoro-uridylate by enzymes encoded by fcy1 (cytosine deaminase) and FUR1 (uracil phosphoribosyltransferase (UPRT)) [65][83]. The FCY1 homologue in C. albicans and other Candida species is the FCA1 gene [136][137][189,190]. The lack of cytosine deaminase in mammalian cells prevents 5-FC conversion and subsequent toxic effects [138][191].
Resistance to 5-FC could arise with mutation or loss of any of the three key enzymes (FCY1, FCY2 or FUR1), as discovered in model organism yeast Saccharomyces cerevisiae [139][140][192,193]. Increased pyrimidine production in the fungal cell can also serve to circumvent toxic antifungal activity [64][65][82,83]. Kern et al. (1991) were among the first to identify the correlation between a point mutation (Arg134Ser) in the FUR1 gene and 5-FC resistance in S. cerevisiae yeast cells [141][194].
7. Biofilm Formation and Antifungal Resistance
7.1. Biofilm Formation during Antifungal Treatment
Over the course of antifungal treatment, a biofilm can develop and enable yeast cells to become more resistant
[142][143][209,215]. Numerous classes of antifungals, including polyenes and azoles, have been cited as less effective over time, even within 72 h of biofilm development/maturation
[142][209]. For fluconazole treatment, resistance in
C. albicans with biofilms can be increased by up to 1000-fold compared to planktonic cells
[143][144][215,217].
Multiple features of biofilms enable antifungal resistance aside from genetic mutations that typically drive resistance in planktonic cells, which adds to the complexity of treating this type of infection
[145][219]. These include physical protection, increased cell density of the microbe community, persister cells and extracellular vesicular secretion
[145][146][147][208,219,220].
7.2. The Roles of β-1,3 Glucan and Biofilm-Associated Antifungal Resistance
One of the main ECM components is β-1,3 glucan, which is synthesized by echinocandin target gene FKS1
[148][223]. The polysaccharide molecules are primarily responsible for the sequestration of antifungal molecules, making them a prime target for treatment of
Candida infections with biofilms
[149][224]. Supplemental treatment with the β-1,3 glucanase enzyme can break down this molecule and subsequently disrupt the biofilm architecture
[150][151][225,226].
7.3. Relevant Antifungal Resistance Genes in Biofilm-Associated Candida Infections
Candida strains associated with biofilm formation and resistance to antifungal treatment may have similar key genes implicated in planktonic resistant isolates. Changes in the expression of ergosterol biosynthesis pathway genes and in biofilm membrane composition appear to be linked to subsequent azole and polyene resistance
[152][153][154][155][230,231,232,233]. Differential gene expression of beta-glucan synthesis-associated genes SKN1 and KRE1 was observed in biofilm-associated-resistant
Candida exposed to amphotericin B, in agreement with previous studies
[152][156][218,230]. These relevant genes may be highlighted in the search for a suitable resistance biomarker.
Drug efflux pumps CDR1, CDR2 and MDR1, which are associated with azole resistance, may also be upregulated in biofilm-associated
Candida strains
[157][158][159][160][161][162][235,236,237,238,239,240]. Interestingly, neither polyenes nor echinocandins have been implicated as a substrate for drug efflux pumps in
Candida infections of this type. This suggests that echinocandins could be a preferable treatment choice over azoles
[158][159][163][236,237,241].
8. Summary
With antifungal resistance being a continued problem, there is a need for the development of quick and reliable molecular diagnostic tests that detect organisms with intrinsic and/or secondary resistance due to the genetic mechanisms [29][34]. Additionally, faster methods of species identification would be useful, given the differences in frequency and impact of antifungal resistance among Candida species. Currently, PCR-based methods using fungal cultures are still the first option for species identification and detection of antifungal resistance in individual strains [124][164][105,167]. The usefulness of real-time testing for resistance in Candida species to modify the treatment course has been well documented.
Resistance can be acquired through a dynamic combination of numerous point mutations and other genetic or transcriptional alterations. A large-scale comparison between matched fluconazole-resistant and -susceptible C. albicans clinical isolates using microarray analysis identified almost two hundred genes (n = 198) that were differentially expressed [165][248]. In resistant isolates, multidrug resistance and oxidative stress response genes, among others, were found to be upregulated compared to susceptible samples [165][248]. This highlights the fact that a dynamic response that leads to antifungal resistance and identification of reliable biomarkers or gene expression profiles that different strains have in common would be beneficial to improve future treatment decisions.
To develop a dependable point-of-care (POCT) diagnostic assay for antifungal resistance, reliable biomarkers in
Candida species need to be established. It is important to note that different non-synonymous substitutions and other genetic alterations may result in similar genetic effects.