RNAi epimutations conferring antifungal drug resistance are inheritable.

Epimutations modify gene expression and lead to phenotypic variation while the encoding DNA sequence remains unchanged. Epimutations mediated by RNA interference (RNAi) and/or chromatin modifications can confer antifungal drug resistance and may impact virulence traits in fungi (1-5). However, whether these epigenetic modifications can be transmitted across generations following sexual reproduction was unclear. This study demonstrates that RNAi epimutations conferring antifungal drug resistance are transgenerationally inherited in the human fungal pathogen Mucor circinelloides . Our research revealed that RNAi-based antifungal resistance follows a DNA sequence-independent, non-Mendelian inheritance pattern. Small RNAs (sRNAs) are the exclusive determinants of inheritance, transmitting drug resistance independently of other known repressive epigenetic modifications. Unique sRNA signature patterns can be traced through inheritance from parent to progeny, further supporting RNA as an alternative molecule for transmitting information across generations. This study marks a significant advance in understanding epigenetic inheritance, highlighting RNAi-exclusive epimutations as a widespread phenomenon in the pathogenic M. circinelloides species complex. Understanding how epimutations occur, propagate, and confer resistance may enable their detection in other eukaryotic pathogens, provide solutions for challenges posed by rising antimicrobial drug resistance (AMR), and also advance research on phenotypic adaptability and its evolutionary implications.

Structure-guided design, synthesis, and characterization of C22- and C32-modified FK520 analogs with enhanced activity against human pathogenic fungi.

Invasive fungal infections are a leading cause of death worldwide. Translating molecular insights into clinical benefits is challenging because fungal pathogens and their hosts share similar eukaryotic physiology. Consequently, current antifungal treatments have limited efficacy, may be poorly fungicidal in the host, can exhibit toxicity, and are increasingly compromised by emerging resistance. We have established that the phosphatase calcineurin (CaN) is required for invasive fungal disease and an attractive target for antifungal drug development. CaN is a druggable target, and there is vast clinical experience with the CaN inhibitors FK506 and cyclosporin A (CsA). However, while FK506 and its natural analog FK520 exhibit antifungal activity, they are also immunosuppressive in the host and thus not fungal-selective. We leverage our pathogenic fungal CaN-FK506-FKBP12 complex X-ray structures and biophysical data to support structure-based ligand design as well as structure-activity relationship analyses of broad-spectrum FK506/FK520 derivatives with potent antifungal activity and reduced immunosuppressive activity. Here we apply molecular docking studies to develop antifungal C22- or C32-modified FK520 derivatives with improved therapeutic index scores. Among them, the C32-modified FK520 derivative JH-FK-44 ( 7 ) demonstrates a significantly improved therapeutic index compared to JH-FK-08, our lead compound to date. NMR binding studies with C32-derivatives are consistent with our hypothesis that C32 modifications disrupt the hydrogen bonding network in the human complex while introducing favorable electrostatic and cation-π interactions with the fungal FKBP12 R86 residue. These findings further reinforce calcineurin inhibition as a promising strategy for antifungal therapy. Significance: Invasive fungal infections cause significant mortality worldwide, and current antifungal treatments are often ineffective, toxic, or face growing resistance. This research identifies calcineurin (CaN), a critical protein for fungal survival, as a potential target for developing new antifungal drugs. Although existing CaN inhibitors such as FK506 (tacrolimus) and FK520 (ascomycin) possess antifungal properties, their immunosuppressive effects limit their clinical utility. By studying the structure of human and fungal FKBP12-FK506 or FK520 complexes with CaN, we have designed and synthesized modified FK520 derivatives with strong antifungal activity and reduced immunosuppressive effects. These new derivatives are expected to have significantly improved therapeutic profiles, offering hope for more effective and safer antifungal treatments.

Focus on fungi.

Fungi play critical roles in the homeostasis of ecosystems globally and have emerged as significant causes of an expanding repertoire of devastating diseases in plants, animals, and humans. In this Commentary, we highlight the importance of fungal pathogens and argue for concerted research efforts to enhance understanding of fungal virulence, antifungal immunity, novel drug targets, antifungal resistance, and the mycobiota to improve human health.

Malassezia responds to environmental pH signals through the conserved Rim/Pal pathway.

During mammalian colonization and infection, microorganisms must be able to rapidly sense and adapt to changing environmental conditions including alterations in extracellular pH. The fungus-specific Rim/Pal signaling pathway is one process that supports microbial adaptation to alkaline pH. This cascading series of interacting proteins terminates in the proteolytic activation of the highly conserved Rim101/PacC protein, a transcription factor that mediates microbial responses that favor survival in neutral/alkaline pH growth conditions, including many mammalian tissues. We identified the putative Rim pathway proteins Rim101 and Rra1 in the human skin colonizing fungus Malassezia sympodialis. Gene deletion by transconjugation and homologous recombination revealed that Rim101 and Rra1 are required for M. sympodialis growth at higher pH. In addition, comparative transcriptional analysis of the mutant strains compared to wild-type suggested mechanisms for fungal adaptation to alkaline conditions. These pH-sensing signaling proteins are required for optimal growth in a murine model of atopic dermatitis, a pathological condition associated with increased skin pH. Together, these data elucidate both conserved and phylum-specific features of microbial adaptation to extracellular stresses.IMPORTANCEThe ability to adapt to host pH has been previously associated with microbial virulence in several pathogenic fungal species. Here we demonstrate that a fungal-specific alkaline response pathway is conserved in the human skin commensal fungus Malassezia sympodialis (Ms). This pathway is characterized by the pH-dependent activation of the Rim101/PacC transcription factor that controls cell surface adaptations to changing environmental conditions. By disrupting genes encoding two predicted components of this pathway, we demonstrated that the Rim/Pal pathway is conserved in this fungal species as a facilitator of alkaline pH growth. Moreover, targeted gene mutation and comparative transcriptional analysis support the role of the Ms Rra1 protein as a cell surface pH sensor conserved within the basidiomycete fungi, a group including plant and human pathogens. Using an animal model of atopic dermatitis, we demonstrate the importance of Ms Rim/Pal signaling in this common inflammatory condition characterized by increased skin pH.

Distinct evolutionary trajectories following loss of RNA interference in Cryptococcus neoformans.

While increased mutation rates typically have negative consequences in multicellular organisms, hypermutation can be advantageous for microbes adapting to the environment. Previously, we identified two hypermutator Cryptococcus neoformans clinical isolates that rapidly develop drug resistance due to transposition of a retrotransposon, Cnl1. Cnl1-mediated hypermutation is caused by a nonsense mutation in the gene encoding a novel RNAi component, Znf3, combined with a tremendous transposon burden. To elucidate adaptative mechanisms following RNAi loss, two bioinformatic pipelines were developed to identify RNAi loss-of-function mutations in a collection of 387 sequenced C. neoformans isolates. Remarkably, several RNAi-loss isolates were identified that are not hypermutators and have not accumulated transposons. To test if these RNAi loss-of-function mutations can cause hypermutation, the mutations were introduced into a non-hypermutator strain with a high transposon burden, which resulted in a hypermutator phenotype. To further investigate if RNAi-loss isolates can become hypermutators, in vitro passaging was performed. Although no hypermutators were found in two C. neoformans RNAi-loss strains after short-term passage, hypermutation was observed in a passaged C. deneoformans strain with increased transposon burden. Consistent with a two-step evolution, when an RNAi-loss isolate was crossed with an isolate containing a high Cnl1 burden, F1 hypermutator progeny inheriting a high transposon burden were identified. In addition to Cnl1 transpositions, insertions of a novel gigantic DNA transposon KDZ1 (~11 kb), contributed to hypermutation in the progeny. Our results suggest that RNAi loss is relatively common (7/387, ~1.8%) and enables distinct evolutionary trajectories: hypermutation following transposon accumulation or survival without hypermutation.

The Cryptococcus neoformans STRIPAK complex controls genome stability, sexual development, and virulence.

The eukaryotic serine/threonine protein phosphatase PP2A is a heterotrimeric enzyme composed of a scaffold A subunit, a regulatory B subunit, and a catalytic C subunit. Of the four known B subunits, the B"' subunit (known as striatin) interacts with the multi-protein striatin-interacting phosphatase and kinase (STRIPAK) complex. Orthologs of STRIPAK components were identified in C. neoformans, namely PP2AA/Tpd3, PP2AC/Pph22, PP2AB"'/Far8, STRIP/Far11, SLMAP/Far9, and Mob3. Structural modeling, protein domain analysis, and detected protein-protein interactions suggest C. neoformans STRIPAK is assembled similarly to the human and fungal orthologs. Here, STRIPAK components Pph22, Far8, and Mob3 were functionally characterized. Whole-genome sequencing revealed that mutations in STRIPAK complex subunits lead to increased segmental and chromosomal aneuploidy, suggesting STRIPAK functions in maintaining genome stability. We demonstrate that PPH22 is a haploinsufficient gene: heterozygous PPH22/pph22Δ mutant diploid strains exhibit defects in hyphal growth and sporulation and have a significant fitness disadvantage when grown in competition against a wild-type diploid. Deletion mutants pph22Δ, far8Δ, and mob3Δ exhibit defects in mating and sexual differentiation, including impaired hyphae, basidia, and basidiospore production. Loss of either PPH22 or FAR8 leads to growth defects at 30°C, severely reduced growth at elevated temperature, abnormal cell morphology, and impaired virulence. The pph22Δ and far8Δ mutants are also unable to grow in the presence of the calcineurin inhibitors cyclosporine A or FK506, and thus these mutations are synthetically lethal with loss of calcineurin activity. Conversely, mob3Δ mutants display increased thermotolerance, capsule production, and melanization, and are hypervirulent in a murine infection model. Taken together, these findings reveal that the C. neoformans STRIPAK complex plays an important role in genome stability, vegetative growth, sexual development, and virulence in this prominent human fungal pathogen.

Alternative ergosterol biosynthetic pathways confer antifungal drug resistance in the human pathogens within the Mucor species complex.

Mucormycoses are emerging fungal infections caused by a variety of heterogeneous species within the Mucorales order. Among the Mucor species complex, Mucor circinelloides is the most frequently isolated pathogen in mucormycosis patients and despite its clinical significance, there is an absence of established genome manipulation techniques to conduct molecular pathogenesis studies. In this study, we generated a spontaneous uracil auxotrophic strain and developed a genetic transformation procedure to analyze molecular mechanisms conferring antifungal drug resistance. With this new model, phenotypic analyses of gene deletion mutants were conducted to define Erg3 and Erg6a as key biosynthetic enzymes in the M. circinelloides ergosterol pathway. Erg3 is a C-5 sterol desaturase involved in growth, sporulation, virulence, and azole susceptibility. In other fungal pathogens, erg3 mutations confer azole resistance because Erg3 catalyzes the production of a toxic diol upon azole exposure. Surprisingly, M. circinelloides produces only trace amounts of this toxic diol and yet, it is still susceptible to posaconazole and isavuconazole due to alterations in membrane sterol composition. These alterations are severely aggravated by erg3Δ mutations, resulting in ergosterol depletion and, consequently, hypersusceptibility to azoles. We also identified Erg6a as the main C-24 sterol methyltransferase, whose activity may be partially rescued by the paralogs Erg6b and Erg6c. Loss of Erg6a function diverts ergosterol synthesis to the production of cholesta-type sterols, resulting in resistance to amphotericin B. Our findings suggest that mutations or epimutations causing loss of Erg6 function may arise during human infections, resulting in antifungal drug resistance to first-line treatments against mucormycosis. IMPORTANCE: The Mucor species complex comprises a variety of opportunistic pathogens known to cause mucormycosis, a potentially lethal fungal infection with limited therapeutic options. The only effective first-line treatments against mucormycosis consist of liposomal formulations of amphotericin B and the triazoles posaconazole and isavuconazole, all of which target components within the ergosterol biosynthetic pathway. This study uncovered M. circinelloides Erg3 and Erg6a as key enzymes to produce ergosterol, a vital constituent of fungal membranes. Absence of any of those enzymes leads to decreased ergosterol and consequently, resistance to ergosterol-binding polyenes such as amphotericin B. Particularly, losing Erg6a function poses a higher threat as the ergosterol pathway is channeled into alternative sterols similar to cholesterol, which maintain membrane permeability. As a result, erg6a mutants survive within the host and disseminate the infection, indicating that Erg6a deficiency may arise during human infections and confer resistance to the most effective treatment against mucormycoses.

Malassezia responds to environmental pH signals through the conserved Rim/Pal pathway.

During mammalian colonization and infection, microorganisms must be able to rapidly sense and adapt to changing environmental conditions including alterations in extracellular pH. The fungus-specific Rim/Pal signaling pathway is one process that supports microbial adaptation to alkaline pH. This cascading series of interacting proteins terminates in the proteolytic activation of the highly conserved Rim101/PacC protein, a transcription factor that mediates microbial responses that favor survival in neutral/alkaline pH growth conditions, including many mammalian tissues. We identified the putative Rim pathway proteins Rim101 and Rra1 in the human skin colonizing fungus Malassezia sympodialis. Gene deletion by transconjugation and homologous recombination revealed that Rim101 and Rra1 are required for M. sympodialis growth at higher pH. Additionally, comparative transcriptional analysis of the mutant strains compared to wild-type suggested mechanisms for fungal adaptation to alkaline conditions. These pH-sensing signaling proteins are required for optimal growth in a murine model of atopic dermatitis, a pathological condition associated with increased skin pH. Together these data elucidate both conserved and phylum-specific features of microbial adaptation to extracellular stresses.

Duplication and neofunctionalization of a horizontally transferred xyloglucanase as a facet of the Red Queen coevolutionary dynamic.

Oomycete protists share phenotypic similarities with fungi, including the ability to cause plant diseases, but branch in a distant region of the tree of life. It has been suggested that multiple horizontal gene transfers (HGTs) from fungi-to-oomycetes contributed to the evolution of plant-pathogenic traits. These HGTs are predicted to include secreted proteins that degrade plant cell walls, a barrier to pathogen invasion and a rich source of carbohydrates. Using a combination of phylogenomics and functional assays, we investigate the diversification of a horizontally transferred xyloglucanase gene family in the model oomycete species Phytophthora sojae. Our analyses detect 11 xyloglucanase paralogs retained in P. sojae. Using heterologous expression in yeast, we show consistent evidence that eight of these paralogs have xyloglucanase function, including variants with distinct protein characteristics, such as a long-disordered C-terminal extension that can increase xyloglucanase activity. The functional variants analyzed subtend a phylogenetic node close to the fungi-to-oomycete transfer, suggesting the horizontally transferred gene was a bona fide xyloglucanase. Expression of three xyloglucanase paralogs in Nicotiana benthamiana triggers high-reactive oxygen species (ROS) generation, while others inhibit ROS responses to bacterial immunogens, demonstrating that the paralogs differentially stimulate pattern-triggered immunity. Mass spectrometry of detectable enzymatic products demonstrates that some paralogs catalyze the production of variant breakdown profiles, suggesting that secretion of variant xyloglucanases increases efficiency of xyloglucan breakdown as well as diversifying the damage-associated molecular patterns released. We suggest that this pattern of neofunctionalization and the variant host responses represent an aspect of the Red Queen host-pathogen coevolutionary dynamic.

Comparative genomics of the closely related fungal genera Cryptococcus and Kwoniella reveals karyotype dynamics and suggests evolutionary mechanisms of pathogenesis.

In exploring the evolutionary trajectories of both pathogenesis and karyotype dynamics in fungi, we conducted a large-scale comparative genomic analysis spanning the Cryptococcus genus, encompassing both global human fungal pathogens and nonpathogenic species, and related species from the sister genus Kwoniella. Chromosome-level genome assemblies were generated for multiple species, covering virtually all known diversity within these genera. Although Cryptococcus and Kwoniella have comparable genome sizes (about 19.2 and 22.9 Mb) and similar gene content, hinting at preadaptive pathogenic potential, our analysis found evidence of gene gain (via horizontal gene transfer) and gene loss in pathogenic Cryptococcus species, which might represent evolutionary signatures of pathogenic development. Genome analysis also revealed a significant variation in chromosome number and structure between the 2 genera. By combining synteny analysis and experimental centromere validation, we found that most Cryptococcus species have 14 chromosomes, whereas most Kwoniella species have fewer (11, 8, 5, or even as few as 3). Reduced chromosome number in Kwoniella is associated with formation of giant chromosomes (up to 18 Mb) through repeated chromosome fusion events, each marked by a pericentric inversion and centromere loss. While similar chromosome inversion-fusion patterns were observed in all Kwoniella species with fewer than 14 chromosomes, no such pattern was detected in Cryptococcus. Instead, Cryptococcus species with less than 14 chromosomes showed reductions primarily through rearrangements associated with the loss of repeat-rich centromeres. Additionally, Cryptococcus genomes exhibited frequent interchromosomal translocations, including intercentromeric recombination facilitated by transposons shared between centromeres. Overall, our findings advance our understanding of genetic changes possibly associated with pathogenicity in Cryptococcus and provide a foundation to elucidate mechanisms of centromere loss and chromosome fusion driving distinct karyotypes in closely related fungal species, including prominent global human pathogens.

Comparative genomics of Cryptococcus and Kwoniella reveals pathogenesis evolution and contrasting karyotype dynamics via intercentromeric recombination or chromosome fusion.

A large-scale comparative genomic analysis was conducted for the global human fungal pathogens within the Cryptococcus genus, compared to non-pathogenic Cryptococcus species, and related species from the sister genus Kwoniella. Chromosome-level genome assemblies were generated for multiple species of both genera, resulting in a dataset encompassing virtually all of their known diversity. Although Cryptococcus and Kwoniella have comparable genome sizes (about 19.2 and 22.9 Mb) and similar gene content, hinting at pre-adaptive pathogenic potential, our analysis found evidence in pathogenic Cryptococcus species of specific examples of gene gain (via horizontal gene transfer) and gene loss, which might represent evolutionary signatures of pathogenic development. Genome analysis also revealed a significant variation in chromosome number and structure between the two genera. By combining synteny analysis and experimental centromere validation, we found that most Cryptococcus species have 14 chromosomes, whereas most Kwoniella species have fewer (11, 8, 5 or even as few as 3). Reduced chromosome number in Kwoniella is associated with formation of giant chromosomes (up to 18 Mb) through repeated chromosome fusion events, each marked by a pericentric inversion and centromere loss. While similar chromosome inversion-fusion patterns were observed in all Kwoniella species with fewer than 14 chromosomes, no such pattern was detected in Cryptococcus. Instead, Cryptococcus species with less than 14 chromosomes, underwent chromosome reductions primarily through rearrangements associated with the loss of repeat-rich centromeres. Additionally, Cryptococcus genomes exhibited frequent interchromosomal translocations, including intercentromeric recombination facilitated by transposons shared between centromeres. Taken together, our findings advance our understanding of genomic changes possibly associated with pathogenicity in Cryptococcus and provide a foundation to elucidate mechanisms of centromere loss and chromosome fusion driving distinct karyotypes in closely related fungal species, including prominent global human pathogens.

Calcineurin contributes to RNAi-mediated transgene silencing and small interfering RNA production in the human fungal pathogen Cryptococcus neoformans.

Adaptation to external environmental challenges at the cellular level requires rapid responses and involves relay of information to the nucleus to drive key gene expression changes through downstream transcription factors. Here, we describe an alternative route of adaptation through a direct role for cellular signaling components in governing gene expression via RNA interference-mediated small RNA production. Calcium-calcineurin signaling is a highly conserved signaling cascade that plays central roles in stress adaptation and virulence of eukaryotic pathogens, including the human fungal pathogen Cryptococcus neoformans. Upon activation in C. neoformans, calcineurin localizes to P-bodies, membraneless organelles that are also the site for RNA processing. Here, we studied the role of calcineurin and its substrates in RNAi-mediated transgene silencing. Our results reveal that calcineurin regulates both the onset and the reversion of transgene silencing. We found that some calcineurin substrates that localize to P-bodies also regulate transgene silencing but in opposing directions. Small RNA sequencing in mutants lacking calcineurin or its targets revealed a role for calcineurin in small RNA production. Interestingly, the impact of calcineurin and its substrates was found to be different in genome-wide analysis, suggesting that calcineurin may regulate small RNA production in C. neoformans through additional pathways. Overall, these findings define a mechanism by which signaling machinery induced by external stimuli can directly alter gene expression to accelerate adaptative responses and contribute to genome defense.

Casein kinase 2 complex: a central regulator of multiple pathobiological signaling pathways in Cryptococcus neoformans.

The casein kinase 2 (CK2) complex has garnered extensive attention over the past decades as a potential therapeutic target for diverse human diseases, including cancer, diabetes, and obesity, due to its pivotal roles in eukaryotic growth, differentiation, and metabolic homeostasis. While CK2 is also considered a promising antifungal target, its role in fungal pathogens remains unexplored. In this study, we investigated the functions and regulatory mechanisms of the CK2 complex in Cryptococcus neoformans, a major cause of fungal meningitis. The cryptococcal CK2 complex consists of a single catalytic subunit, Cka1, and two regulatory subunits, Ckb1 and Ckb2. Our findings show that Cka1 plays a primary role as a protein kinase, while Ckb1 and Ckb2 have major and minor regulatory functions, respectively, in growth, cell cycle control, morphogenesis, stress response, antifungal drug resistance, and virulence factor production. Interestingly, triple mutants lacking all three subunits (cka1Δ ckb1Δ ckb2Δ) exhibited more severe phenotypic defects than the cka1Δ mutant alone, suggesting that Ckb1/2 may have Cka1-independent functions. In a murine model of systemic cryptococcosis, cka1Δ and cka1Δ ckb1Δ ckb2Δ mutants showed severely reduced virulence. Transcriptomic, proteomic, and phosphoproteomic analyses further revealed that the CK2 complex controls a wide array of effector proteins involved in transcriptional regulation, cell cycle control, nutrient metabolisms, and stress responses. Most notably, CK2 disruption led to dysregulation of key signaling cascades central to C. neoformans pathogenicity, including the Hog1, Mpk1 MAPKs, cAMP/PKA, and calcium/calcineurin signaling pathways. In summary, our study provides novel insights into the multifaceted roles of the fungal CK2 complex and presents a compelling case for targeting it in the development of new antifungal drugs.IMPORTANCEThe casein kinase 2 (CK2) complex, crucial for eukaryotic growth, differentiation, and metabolic regulation, presents a promising therapeutic target for various human diseases, including cancer, diabetes, and obesity. Its potential as an antifungal target is further highlighted in this study, which explores CK2's functions in C. neoformans, a key fungal meningitis pathogen. The CK2 complex in C. neoformans, comprising the Cka1 catalytic subunit and Ckb1/2 regulatory subunits, is integral to processes like growth, cell cycle, morphogenesis, stress response, drug resistance, and virulence. Our findings of CK2's role in regulating critical signaling pathways, including Hog1, Mpk1 MAPKs, cAMP/PKA, and calcium/calcineurin, underscore its importance in C. neoformans pathogenicity. This study provides valuable insights into the fungal CK2 complex, reinforcing its potential as a target for novel antifungal drug development and pointing out a promising direction for creating new antifungal agents.

Calcineurin contributes to RNAi-mediated transgene silencing and small interfering RNA production in the human fungal pathogen Cryptococcus neoformans.

Adaptation to external environmental challenges at the cellular level requires rapid responses and involves relay of information to the nucleus to drive key gene expression changes through downstream transcription factors. Here, we describe an alternative route of adaptation through a direct role for cellular signaling components in governing gene expression via RNA interference-mediated small RNA production. Calcium-calcineurin signaling is a highly conserved signaling cascade that plays central roles in stress adaptation and virulence of eukaryotic pathogens, including the human fungal pathogen Cryptococcus neoformans. Upon activation in C. neoformans, calcineurin localizes to P-bodies, membrane-less organelles that are also the site for RNA processing. Here, we studied the role of calcineurin and its substrates in RNAi-mediated transgene silencing. Our results reveal that calcineurin regulates both the onset and the reversion of transgene silencing. We found that some calcineurin substrates that localize to P-bodies also regulate transgene silencing but in opposing directions. Small RNA sequencing in mutants lacking calcineurin or its targets revealed a role for calcineurin in small RNA production. Interestingly, the impact of calcineurin and its substrates was found to be different in genome-wide analysis, suggesting that calcineurin may regulate small RNA production in C. neoformans through additional pathways. Overall, these findings define a mechanism by which signaling machinery induced by external stimuli can directly alter gene expression to accelerate adaptative responses and contribute to genome defense.

Alternative ergosterol biosynthetic pathways confer antifungal drug resistance in the human pathogens within the Mucor species complex.

Mucormycoses are emerging fungal infections caused by a variety of heterogeneous species within the order Mucorales. Among the Mucor species complex, Mucor circinelloides is the most frequently isolated pathogen in mucormycosis patients and despite its clinical significance, there is an absence of established genome manipulation techniques to conduct molecular pathogenesis studies. In this study, we generated a spontaneous uracil auxotrophic strain and developed a genetic transformation procedure to analyze molecular mechanisms conferring antifungal drug resistance. With this new model, phenotypic analyses of gene deletion mutants were conducted to define Erg3 and Erg6a as key biosynthetic enzymes in the M. circinelloides ergosterol pathway. Erg3 is a C-5 sterol desaturase involved in growth, sporulation, virulence, and azole susceptibility. In other fungal pathogens, erg3 mutations confer azole resistance because Erg3 catalyzes the production of a toxic diol upon azole exposure. Surprisingly, M. circinelloides produces only trace amounts of this toxic diol and yet, it is still susceptible to posaconazole and isavuconazole due to alterations in membrane sterol composition. These alterations are severely aggravated by erg3 Δ mutations, resulting in ergosterol depletion and consequently, hypersensitivity to azoles. We also identified Erg6a as the main C-24 sterol methyltransferase, whose activity may be partially rescued by the paralogs Erg6b and Erg6c. Loss of Erg6a function diverts ergosterol synthesis to the production of cholestatype sterols, resulting in resistance to amphotericin B. Our findings suggest that mutations or epimutations causing loss of Erg6 function may arise during human infections, resulting in antifungal drug resistance to first-line treatments against mucormycoses. Importance: The Mucor species complex includes a variety of opportunistic pathogens known to cause mucormycosis, a potentially lethal fungal infection with limited therapeutic options. The only effective first-line treatments against mucormycosis consist of liposomal formulations of amphotericin B and the triazoles posaconazole and isavuconazole, all of which target components within the ergosterol biosynthetic pathway. This study uncovered M. circinelloides Erg3 and Erg6a as key enzymes to produce ergosterol, a vital constituent of fungal membranes. Absence of any of those enzymes leads to decreased ergosterol and consequently, resistance to ergosterol-binding polyenes such as amphotericin B. Particularly, losing Erg6a function pose a higher threat as the ergosterol pathway is channeled into alternative sterols similar to cholesterol, which maintain membrane permeability. As a result, erg6a mutants survive within the host and disseminate the infection, indicating that Erg6a deficiency may arise during human infections and confer resistance to the most effective treatment against mucormycoses.

Amoeba predation of Cryptococcus: A quantitative and population genomic evaluation of the accidental pathogen hypothesis.

The "Amoeboid Predator-Fungal Animal Virulence Hypothesis" posits that interactions with environmental phagocytes shape the evolution of virulence traits in fungal pathogens. In this hypothesis, selection to avoid predation by amoeba inadvertently selects for traits that contribute to fungal escape from phagocytic immune cells. Here, we investigate this hypothesis in the human fungal pathogens Cryptococcus neoformans and Cryptococcus deneoformans. Applying quantitative trait locus (QTL) mapping and comparative genomics, we discovered a cross-species QTL region that is responsible for variation in resistance to amoeba predation. In C. neoformans, this same QTL was found to have pleiotropic effects on melanization, an established virulence factor. Through fine mapping and population genomic comparisons, we identified the gene encoding the transcription factor Bzp4 that underlies this pleiotropic QTL and we show that decreased expression of this gene reduces melanization and increases susceptibility to amoeba predation. Despite the joint effects of BZP4 on amoeba resistance and melanin production, we find no relationship between BZP4 genotype and escape from macrophages or virulence in murine models of disease. Our findings provide new perspectives on how microbial ecology shapes the genetic architecture of fungal virulence, and suggests the need for more nuanced models for the evolution of pathogenesis that account for the complexities of both microbe-microbe and microbe-host interactions.

Duplication and neofunctionalization of a horizontally-transferred xyloglucanase as a facet of the red queen co-evolutionary dynamic.

Oomycetes are heterotrophic protists that share phenotypic similarities with fungi, including the ability to cause plant diseases, but branch in a separate and distant region of the eukaryotic tree of life. It has been suggested that multiple horizontal gene transfers (HGTs) from fungi-to-oomycetes contributed to the evolution of plant-pathogenic traits. These HGTs are predicted to include secreted proteins that degrade plant cell walls. This is a key trait in the pathology of many oomycetes, as the plant cell wall represents a primary barrier to pathogen invasion and a rich source of carbohydrates. Many of the HGT gene families identified have undergone multiple rounds of duplication. Using a combination of phylogenomic analysis and functional assays, we investigate the diversification of a horizontally-transferred xyloglucanase gene family in the model oomycete species Phytophthora sojae. Our analyses detect 11 genes retained in P. sojae among a complex pattern of gene duplications and losses. Using a phenotype assay, based on heterologous expression in yeast, we show that eight of these paralogs have xyloglucanase function, including variants with distinct protein characteristics, such as a long-disordered C-terminal extension that can increase xyloglucanase activity. The functional xyloglucanase variants analysed subtend an ancestral node close to the fungi-oomycetes gene transfer, suggesting the horizontally-transferred gene was a bona fide xyloglucanase. Expression of xyloglucanase paralogs in Nicotiana benthamiana triggers distinct patterns of reactive oxygen species (ROS) generation, demonstrating that enzyme variants differentially stimulate pattern-triggered immunity in plants. Mass spectrometry of detectable enzymatic products demonstrates that some paralogs catalyze production of variant breakdown profiles, suggesting that secretion of multiple xyloglucanase variants increases efficiency of xyloglucan breakdown, as well as potentially diversifying the range of Damage-Associated Molecular Patterns (DAMPs) released during pathogen attack. We suggest that such patterns of protein neofunctionalization, and variant host responses, represent an aspect of the Red Queen host-pathogen co-evolutionary dynamic. Significance Statement: The oomycetes are a diverse group of eukaryotic microbes that include some of the most devastating pathogens of plants. Oomycetes perceive, invade, and colonize plants in similar ways to fungi, in part because they acquired the genes to attack and feed on plants from fungi. These genes are predicted to be useful to oomycete plant pathogens because they have undergone multiple rounds of gene duplication. One key enzyme for attacking plant cell wall structures is called xyloglucanase. Xyloglucanase in the oomycetes has undergone multiple rounds of gene duplication, leading to variants including an enzyme with a C-terminal extension that increases activity. Some xyloglucanase variants trigger unique patterns of reactive oxygen species (ROS) in planta, and generate different profiles of cell wall breakdown products - such outcomes could act to mystify and increase the workload of the plant immune system, allowing successful pathogens to proliferate.

Enhanced fungal specificity and in vivo therapeutic efficacy of a C-22-modified FK520 analog against C. neoformans.

IMPORTANCE: Fungal infections cause significant morbidity and mortality globally. The therapeutic armamentarium against these infections is limited, and the development of antifungal drugs has been hindered by the evolutionary conservation between fungi and the human host. With rising resistance to the current antifungal arsenal and an increasing at-risk population, there is an urgent need for the development of new antifungal compounds. The FK520 analogs described in this study display potent antifungal activity as a novel class of antifungals centered on modifying an existing orally active FDA-approved therapy. This research advances the development of much-needed newer antifungal treatment options with novel mechanisms of action.

Frequent transitions in mating-type locus chromosomal organization in Malassezia and early steps in sexual reproduction.

Fungi in the basidiomycete genus Malassezia are the most prevalent eukaryotic microbes resident on the skin of human and other warm-blooded animals and have been implicated in skin diseases and systemic disorders. Analysis of Malassezia genomes revealed that key adaptations to the skin microenvironment have a direct genomic basis, and the identification of mating/meiotic genes suggests a capacity to reproduce sexually, even though no sexual cycle has yet been observed. In contrast to other bipolar or tetrapolar basidiomycetes that have either two linked mating-type-determining (MAT) loci or two MAT loci on separate chromosomes, in Malassezia species studied thus far the two MAT loci are arranged in a pseudobipolar configuration (linked on the same chromosome but capable of recombining). By generating additional chromosome-level genome assemblies, and an improved Malassezia phylogeny, we infer that the pseudobipolar arrangement was the ancestral state of this group and revealed six independent transitions to tetrapolarity, seemingly driven by centromere fission or translocations in centromere-flanking regions. Additionally, in an approach to uncover a sexual cycle, Malassezia furfur strains were engineered to express different MAT alleles in the same cell. The resulting strains produce hyphae reminiscent of early steps in sexual development and display upregulation of genes associated with sexual development as well as others encoding lipases and a protease potentially relevant for pathogenesis of the fungus. Our study reveals a previously unseen genomic relocation of mating-type loci in fungi and provides insight toward the identification of a sexual cycle in Malassezia, with possible implications for pathogenicity.