Genetic Changes in Experimental Populations of a Hybrid in the Cryptococcus neoformans Species Complex

Uncovering the role of mitochondrial genome in pathogenicity and drug resistance in pathogenic fungi

Fungal infections are becoming more prevalent globally, particularly affecting immunocompromised populations, such as people living with HIV, organ transplant recipients and those on immunomodulatory therapy. Globally, approximately 6.55 million people are affected by invasive fungal infections annually, leading to serious health consequences and death. Mitochondria are membrane-bound organelles found in almost all eukaryotic cells and play an important role in cellular metabolism and energy production, including pathogenic fungi. These organelles possess their own genome, the mitochondrial genome, which is usually circular and encodes proteins essential for energy production. Variation and evolutionary adaptation within and between species' mitochondrial genomes can affect mitochondrial function, and consequently cellular energy production and metabolic activity, which may contribute to pathogenicity and drug resistance in certain fungal species. This review explores the link between the mitochondrial genome and mechanisms of fungal pathogenicity and drug resistance, with a particular focus on Cryptococcus neoformans and Candida albicans. These insights deepen our understanding of fungal biology and may provide new avenues for developing innovative therapeutic strategies.

Genomic analyses reveal high diversity and rapid evolution of Pichia kudriavzevii within a neonatal intensive care unit in Delhi, India

Pichia kudriavzevii causes life-threatening infections in immunocompromised hosts, including hospitalized neonates. This pathogen is intrinsically resistant to fluconazole, while uncommon P. kudriavzevii strains resistant to multiple antifungal drugs, including voriconazole, amphotericin B, and echinocandins, have also been reported from healthcare environments. Thus, understanding how P. kudriavzevii spread, persist, and adapt to healthcare settings could help us develop better infection management strategies. In this study, whole genome sequencing identifies multiple outbreaks of bloodstream infections in a single neonatal intensive care unit (NICU) over 5 years caused by genetically diverse strains of P. kudriavzevii. Interestingly, two genetically distinct clusters of P. kudriavzevii strains showed frequent loss of heterozygosity (LOH) events between two temporal samples. The first outbreak cluster (2015-2016) showed LOH at chromosomes 1, 4, and 5, and the other outbreak cluster (2020) exhibited LOH at chromosome 2. The circulation of two separate strain clusters of P. kudriavzevii suggests nosocomial transmission in the NICU in different time periods. Furthermore, we compared the transcriptomic profiles of three isolates of clusters I and II that exhibited distinct fluconazole and itraconazole MICs. While no significant difference in gene expression was found at the azole-target gene ERG11 or the ATP-binding cassette (ABC) transporter genes, such differences were found in genes involved in cell division and filamentation, such as SIR2 (sirtuin deacetylase) and RFA1 (replication factor A). Interestingly, increased filamentation was observed in clade I isolate exhibiting high fluconazole MICs. Together, our study indicates significant diversity, persistence, and rapid evolution of P. kudriavzevii within a single NICU.

Impact of inter-species hybridisation on antifungal drug response in the Saccharomyces genus

BACKGROUND

Antifungal drug resistance presents one of the major concerns for global public health, and hybridization allows the development of high fitness organisms that can better survive in restrictive conditions or in presence of antifungal agents. Hence, understanding how allelic variation can influence antifungal susceptibility in hybrid organisms is important for the development of targeted treatments. Here, we exploited recent advances in multigenerational breeding of hemiascomycete hybrids to study the impact of hybridisation on antifungal resistance and identify quantitative trait loci responsible for the phenotype.

RESULTS

The offspring of Saccharomyces cerevisiae x S. kudriavzevii hybrids were screened in the presence of six antifungal drugs and revealed a broad phenotypic diversity across the progeny. QTL analysis was carried out comparing alleles between pools of high and low fitness offspring, identifying hybrid-specific genetic regions involved in resistance to fluconazole, micafungin and flucytosine. We found both drug specific and pleiotropic regions, including 41 blocks containing genes not previously associated with resistance phenotypes. We identified linked genes that influence the same trait, namely a hybrid specific 'super' QTL, and validated, via reciprocal hemizygosity analysis, two causal genes, BCK2 and DNF1. The co-location of genes with similar phenotypic impact supports the notion of an adaption process that limits the segregation of advantageous alleles via recombination.

CONCLUSIONS

This study demonstrates the value of QTL studies to elucidate the hybrid-specific mechanisms of antifungal susceptibility. We also show that an inter-species hybrid model system in the Saccharomyces background, can help to decipher the trajectory of antifungal drug resistance in pathogenic hybrid lineages.

Study on the chemical composition and anti-fungi activities of anthraquinones and its glycosides from Rumex japonicus Houtt

Fungi, such as Trichophyton rubrum (T. rubrum) and Microsporum canis Bodin Anamorph (M. canis Bodin Anamorph) are the main pathogens of dermatophysis. According to ancient books records, Rumex japonicus Houtt. (RJH) has a miraculous effect on the treatment of dermatophysis. To reveal the anti-fungi (T. rubrum and M. canis Bodin Anamorph) components and its mechanism of the Rumex japonicus Houtt. The vinegar extraction and alcohol precipitation, HPLC and nuclear magnetic resonance spectroscopy (NMR) were employed for analyzing the chemical compositions of RJH; in vitro anti-fungal experiment was investigated including test the minimum inhibitory concentration (MIC) and the minimum fungicidal concentration (MFC), spore germination rate, nucleic acid, protein leakage rate, biofilm structure, and the mechanism of anti-fungal and anti-fungal biofilms in RJH. Seven anthraquinones and their glycoside compounds were obtained in this study respectively, such as chrysophanol, physcion, aloe-emodin, emodin, rhein, emodin-8-O-β-D-glucoside and chrysophanol-8-O-β-D-glucoside. In vitro anti-fungal experiment results showed that RJH extracts have good anti-fungal activity for dermatophytic fungi. Among them, the MIC of the rhein, emodin and aloe-emodin against T. rubrum are 1.9 µg/ml, 3.9 µg/ml and 15.6 µg/ml, respectively; the MIC of emodin and aloe-emodin against M. canis Bodin Anamorph are 7.8 µg/ml and 62.5 µg/ml, respectively. In addition, its active components can inhibit fungal spore germination and the formation of bud tube, change cell membrane permeability, prevent hyphal growth, destroy biofilm structure, and down-regulate the expression of agglutinin-like sequence family 1 of the adhesion phase of biofilm growth. The study shows that RJH play a fungicidal role.

Multi-locus sequence typing and phylogenetics of Cryptococcus neoformans AD hybrids

Hybrid AD strains of the human pathogenic Cryptococcus neoformans species complex have been reported from many parts of the world. However, their origin, diversity, and evolution are incompletely understood. In this study, we analyzed 102 AD hybrid strains representing 21 countries on five continents. For each strain, we obtained its mating type and its allelic sequences at each of the seven loci that have been used for genotyping haploid serotypes A and D strains of the species complex by the Cryptococcus research community. Our results showed that most AD hybrids exhibited loss of heterozygosity at one or more of the seven analyzed loci. Phylogenetic and population genetic analyses of the allelic sequences revealed multiple origins of the hybrids within each continent, dating back to one million years ago in Africa and up to the present in other continents. We found evidence for clonal reproduction and long-distance dispersal of these hybrids in nature. Comparisons with the global haploid serotypes A and D strains identified new alleles and new haploid multi-locus genotypes in AD hybrids, consistent with the presence of yet-to-be discovered genetic diversity in haploid populations of this species complex in nature. Together, our results indicate that AD hybrids can be effectively genotyped using the same multi-locus sequencing type approach as that established for serotypes A and D strains. Our comparisons of the AD hybrids among each other as well as with the global haploid serotypes A and D strains revealed novel genetic diversity as well as evidence for multiple origins and dynamic evolution of these hybrids in nature.

Genomic Analyses of a Fungemia Outbreak Caused by Lodderomyces elongisporus in a Neonatal Intensive Care Unit in Delhi, India

Lodderomyces elongisporus is a recently emerging yeast pathogen predominantly reported in adult patients who had immunosuppression and/or intravenous access devices. Here, we report a fungemia outbreak caused by L. elongisporus in a neonatal intensive care unit (NICU) in Delhi, India, from September 2021 to February 2022. All 10 neonates had low birthweight, and nine of the patients survived after amphotericin B treatment. Whole-genome sequence analyses of the patient isolates as well as those from other sources in India grouped them into two clusters: one cluster consists of isolates exclusively from stored apples and the other cluster includes isolates from patients, clinical environments, and stored apples. All outbreak strains from patients were closely related to each other and showed highly similar heterozygosity patterns across all 11 major scaffolds. While overall very similar, strains from the inanimate environment of the same neonatal intensive care unit showed loss of heterozygosity at scaffold 2 (NW_001813676) compared to the patient strains. Interestingly, evidence for recombination was found in all samples. All clinical strains were susceptible to 10 tested antifungal drugs, and comparisons with strains with high fluconazole MICs derived from the surface of stored apples revealed significant genome divergence between the clinical and apple surface strains, including 119 nonsynonymous single nucleotide polymorphisms (SNPs) in 24 triazole resistance-related genes previously found in other Candida spp. Together, our results indicate significant diversity, recombination, and persistence in the hospital setting and a high rate of evolution in this emerging yeast pathogen. IMPORTANCE Lodderomyces elongisporus was initially considered a teleomorph of Candida parapsilosis. However, DNA sequence analyses revealed it as a distinctive species. Invasive infections due to L. elongisporus have been reported globally. We report an outbreak of fungemia due to L. elongisporus in a NICU affecting 10 preterm, low-birthweight neonates during a period of 6 months. The outbreak investigation identified two environmental sites, the railing and the temperature panel of the neonate open care warmer, harboring L. elongisporus. Whole-genome sequencing confirmed that the neonate isolates were closely related to each other whereas strains from the inanimate clinical environment were related to clinical strains but showed a marked loss of heterozygosity. Further, L. elongisporus strains recovered previously from the surface of stored apples showed high fluconazole MICs and alterations in triazole resistance-related genes. Genome-wide SNP comparisons revealed recombination as an important source for genomic diversity during adaptation of L. elongisporus to different environments.

Loss of Heterozygosity and Its Importance in Evolution

Loss of heterozygosity (LOH) is a mitotic recombination event that converts heterozygous loci to homozygous loci. This mutation event is widespread in organisms that have asexual reproduction like budding yeasts, and is also an important and frequent mutation event in tumorigenesis. Mutation accumulation studies have demonstrated that LOH occurs at a rate higher than the point mutation rate, and can impact large portions of the genome. Laboratory evolution experiments of heterozygous yeasts have revealed that LOH often unmasks beneficial recessive alleles that can confer large fitness advantages. Here, I highlight advances in understanding dominance, fitness, and phenotypes in laboratory evolved heterozygous yeast strains. I discuss best practices for detecting LOH in intraspecific and interspecific evolved clones and populations. Utilizing heterozygous strain backgrounds in laboratory evolution experiments offers an opportunity to advance our understanding of this important mutation type in shaping adaptation and genome evolution in wild, domesticated, and clinical populations.

Analyses of the Global Multilocus Genotypes of the Human Pathogenic Yeast Cryptococcus neoformans Species Complex

Cryptococcus neoformans species complex (CNSC) is a globally distributed human opportunistic yeast pathogen consisting of five major molecular types (VNI, VNII, VNB, VNIII and VNIV) belonging to two species, C. neoformans (VNI, VNII and VNB, collectively called serotype A) and C. deneoformans (VNIV, commonly called serotype D), and their hybrids (VNIII, serotype AD). Over the years, many studies have analyzed the geographical distribution and genetic diversity of CNSC. However, the global population structure and mode of reproduction remain incompletely described. In this study, we analyze the published multilocus sequence data at seven loci for CNSC. The combined sequences at the seven loci identified a total of 657 multilocus sequence types (STs), including 296 STs with known geographic information, representing 4200 non-redundant isolates from 31 countries and four continents. Among the 296 STs, 78 and 52 were shared among countries and continents, respectively, representing 3643 of the 4200 isolates. Except for the clone-corrected serotype D sample among countries, our analysis of the molecular variance of the 4200 isolates revealed significant genetic differentiations among countries and continents in populations of CNSC, serotype A, and serotype D. Phylogenetic analyses of the concatenated sequences of all 657 STs revealed several large clusters corresponding to the major molecular types. However, several rare but distinct STs were also found, representing potentially novel molecular types and/or hybrids of existing molecular types. Phylogenetic incompatibility analyses revealed evidence for recombination within all four major molecular types-VNI, VNII, VNIV and VNB-as well as within two VNB subclades, VNBI and VNBII, and two ST clusters around the most common STs, ST5 and ST93. However, linkage disequilibrium analyses rejected the hypothesis of random recombination across most samples. Together, our results suggest evidence for historical differentiation, frequent recent gene flow, clonal expansion and recombination within and between lineages of the global CNSC population.

Synthesis and Evaluation of the Antifungal and Toxicological Activity of Nitrofuran Derivatives

Fungal diseases affect more than 1 billion people worldwide. The constant global changes, the advent of new pandemics, and chronic diseases favor the diffusion of fungal pathogens such as Candida, Cryptococcus, Aspergillus, Trichophyton, Histoplasma capsulatum, and Paracoccidioides brasiliensis. In this work, a series of nitrofuran derivatives were synthesized and tested against different fungal species; most of them showed inhibitory activity, fungicide, and fungistatic profile. The minimal inhibitory concentration (MIC90) values for the most potent compounds range from 0.48 µg/mL against H. capsulatum (compound 11) and P. brasiliensis (compounds 3 and 9) to 0.98 µg/mL against Trichophyton rubrum and T. mentagrophytes (compounds 8, 9, 12, 13 and 8, 12, 13, respectively), and 3.9 µg/mL against Candida and Cryptococcus neoformans strains (compounds 1 and 5, respectively). In addition, all compounds showed low toxicity when tested in vitro on lung cell lines (A549 and MRC-5) and in vivo in Caenorhabditis elegans larvae. Many of them showed high selectivity index values. Thus, these studied nitrofuran derivatives proved to be potent against different fungal species, characterized by low toxicity and high selectivity; for these reasons, they may become promising compounds for the treatment of mycoses.

Clonal evolution in serially passaged Cryptococcus neoformans × deneoformans hybrids reveals a heterogenous landscape of genomic change

Cryptococcus neoformans × deneoformans hybrids (also known as serotype AD hybrids) are basidiomycete yeasts that are common in a clinical setting. Like many hybrids, the AD hybrids are largely locked at the F1 stage and are mostly unable to undergo normal meiotic reproduction. However, these F1 hybrids, which display a high (∼10%) sequence divergence are known to genetically diversify through mitotic recombination and aneuploidy, and this diversification may be adaptive. In this study, we evolved a single AD hybrid genotype in six diverse environments by serial passaging and then used genome resequencing of evolved clones to determine evolutionary mechanisms of adaptation. The evolved clones generally increased fitness after passaging, accompanied by an average of 3.3 point mutations, 2.9 loss of heterozygosity (LOH) events, and 0.7 trisomic chromosomes per clone. LOH occurred through nondisjunction of chromosomes, crossing over consistent with break-induced replication, and gene conversion, in that order of prevalence. The breakpoints of these recombination events were significantly associated with regions of the genome with lower sequence divergence between the parents and clustered in sub-telomeric regions, notably in regions that had undergone introgression between the two parental species. Parallel evolution was observed, particularly through repeated homozygosity via nondisjunction, yet there was little evidence of environment-specific parallel change for either LOH, aneuploidy, or mutations. These data show that AD hybrids have both a remarkable genomic plasticity and yet are challenged in the ability to recombine through sequence divergence and chromosomal rearrangements, a scenario likely limiting the precision of adaptive evolution to novel environments.

Phenotypic Plasticity in the Productions of Virulence Factors Within and Among Serotypes in the Cryptococcus neoformans Species Complex

The Cryptococcus neoformans species complex (CNSC) is a common opportunistic human fungal pathogen and the most frequent cause of fungal meningitis. There are three major serotypes in CNSC: A, D, and their hybrids AD, and they have different geographic distributions and medical significance. Melanin pigment and a polysaccharide capsule are the two major virulence factors in CNSC. However, the relationships between serotype and virulence factor production and how environmental factors might impact their relationships are not known. This study investigated the expressions of melanin and capsular polysaccharide in a genetically diverse group of CNSC strains and how their phenotypic expressions were influenced by oxidative and nitrosative stress levels. We found significant differences in melanin and capsular polysaccharide productions among serotypes and across stress conditions. Under oxidative stress, the laboratory hybrids exhibited the highest phenotypic plasticity for melanin production while serotype A showed the highest for capsular polysaccharide production. In contrast, serotype D exhibited the highest phenotypic plasticity for capsular polysaccharide production and clinical serotype AD the highest phenotypic plasticity for melanin production under nitrosative stress. These results demonstrated that different serotypes have different environmental condition-specific mechanisms to modulate the expression of virulence factors.

Genetic and Phenotypic Diversities in Experimental Populations of Diploid Inter-Lineage Hybrids in the Human Pathogenic Cryptococcus

To better understand the potential factors contributing to genome instability and phenotypic diversity, we conducted mutation accumulation (MA) experiments for 120 days for 7 diploid cryptococcal hybrids under fluconazole (10 MA lines each) and non-fluconazole conditions (10 MA lines each). The genomic DNA content, loss of heterozygosity (LOH) rate, growth ability, and fluconazole susceptibility were determined for all 140 evolved cultures. Compared to that of their ancestral clones, the evolved clones showed: (i) genomic DNA content changes ranging from ~22% less to ~27% more, and (ii) reduced, similar, and increased phenotypic values for each tested trait, with most evolved clones displaying increased growth at 40 °C and increased fluconazole resistance. Aside from the ancestral multi-locus genotypes (MLGs) and heterozygosity patterns (MHPs), 77 unique MLGs and 70 unique MPHs were identified among the 140 evolved cultures at day 120. The average LOH rates of the MA lines in the absence and presence of fluconazole were similar at 1.27 × 10⁻⁴ and 1.38 × 10⁻⁴ LOH events per MA line per mitotic division, respectively. While LOH rates varied among MA lines from different ancestors, there was no apparent correlation between the genetic divergence of the parental haploid genomes within ancestral clones and LOH rates. Together, our results suggest that hybrids between diverse lineages of the human pathogenic Cryptococcus can generate significant genotypic and phenotypic diversities during asexual reproduction.

Pathogenesis of Fungal and Bacterial Microbes

The pathogenesis of fungal and bacterial microbes is a complex process involving distinct parameters, including virulence factors, nutrient sensing and availability, microbial signals, as well as host status and defense responses. Defining pathogenesis improves our understanding of how an organism causes diseases and provides insight into novel prospects to combat infection. The effects of pathogenic microbes have significant impact on diverse sectors, including health, agriculture, and economics, underscoring their immense importance in society. Articles in this Special Issue address unique aspects of microbial pathogenesis by exploring interactions between host and pathogen during infection, defining inflammatory immune responses, profiling the importance of essential microbial structures associated with virulence, and outlining critical considerations driving complex diseases.

Fungal species concepts in the genomics era

The 140 000 or so fungal species reported so far are heterogeneously defined based on varying criteria such as morphological, physiological, mating, and (or) molecular features. Incongruences are common among traits used to separating closely related species and it is often difficult to compare fungal taxonomic groups defined based on different species recognition criteria. Though DNA sequence-based classification and identification have been made, a consensus has not been reached, primarily due to intrinsic limitations in the proposed one or a few genes. Here, I argue that the fundamental reason for the observed inconsistencies is that speciation is a stochastic process with the emergence and fixation of different traits influenced differently by many non-deterministic factors such as population size, random mutation, mode(s) of reproduction, selection imposed by interacting biotic and abiotic factors, and chance events. Each species concept attempts to capture one or a few traits emerged in the continuous process of speciation. I propose that a genome sequence-based classification and identification system could unify and stabilize fungal taxonomy and help integrate taxonomy with other fields of fungal biology. The genomic species concept could be similarly argued for other groups of eukaryotic microbes as well as for plants and animals.

Mitochondrial Genome Polymorphisms in the Human Pathogenic Fungus Cryptococcus neoformans

The Cryptococcus complex consists of at least seven evolutionary divergent lineages and causes ∼200,000 fatal human infections each year worldwide. The dominant lineage is Cryptococcus neoformans which consists of three haploid clades VNI, VNII, and VNB, their haploid hybrids, and various diploids derived from intra- and inter-clade mating events. In this study, we analyzed the mitogenomes of 184 strains of C. neoformans. Our analyses revealed that all 184 mitogenomes contained the same 15 protein-coding genes in the same gene order. However, their mitogenome sizes varied between 24,740 and 31,327 bp, primarily due to differences in the number and size of mitochondrial introns. Twelve nucleotide sites within five mitochondrial genes were found to contain introns in at least one of the 184 strains, ranging from 2 to 7 introns within each mitogenome. The concatenated mitochondrial exon sequences of the 15 protein-coding genes and two rRNA genes showed that VNI, VNII, and VNB strains were separated into distinct clades or sub-clades, largely consistent with results based on nuclear genome SNPs. However, several novel findings were observed. First, one strain of the VNB clade contained mitogenome exon sequences identical to the main VNI mitogenome type but was distant to other VNB mitogenomes. Second, hybrids among clades VNI, VNII, and VNB identified based on their nuclear genome SNPs contained mitogenomes from different clades, with evidence of their mitogenomes inherited from either the MAT a or the MAT α parents. Third, the eight diploid VNB (C. neoformans) × VNIV (C. deneoformans) hybrids contained recombinant mitogenomes. Fourth, analyses of intron distribution and the paired exon-intron phylogenies for each of the 12 exon-intron pairs suggested frequent gains and losses of mitochondrial introns during the evolution of C. neoformans. The combined mitogenome exon-based phylogeny and intron distributions suggested that clades VNI, VNII and VNB could be further divided into sub-clades. Together, our results revealed a dynamic evolution of mitochondrial genomes in this important human fungal pathogen.

Hybridization Facilitates Adaptive Evolution in Two Major Fungal Pathogens

Hybridization is increasingly recognized as an important force impacting adaptation and evolution in many lineages of fungi. During hybridization, divergent genomes and alleles are brought together into the same cell, potentiating adaptation by increasing genomic plasticity. Here, we review hybridization in fungi by focusing on two fungal pathogens of animals. Hybridization is common between the basidiomycete yeast species Cryptococcus neoformans × Cryptococcus deneoformans, and hybrid genotypes are frequently found in both environmental and clinical settings. The two species show 10-15% nucleotide divergence at the genome level, and their hybrids are highly heterozygous. Though largely sterile and unable to mate, these hybrids can propagate asexually and generate diverse genotypes by nondisjunction, aberrant meiosis, mitotic recombination, and gene conversion. Under stress conditions, the rate of such genetic changes can increase, leading to rapid adaptation. Conversely, in hybrids formed between lineages of the chytridiomycete frog pathogen Batrachochytrium dendrobatidis (Bd), the parental genotypes are considerably less diverged (0.2% divergent). Bd hybrids are formed from crosses between lineages that rarely undergo sex. A common theme in both species is that hybrids show genome plasticity via aneuploidy or loss of heterozygosity and leverage these mechanisms as a rapid way to generate genotypic/phenotypic diversity. Some hybrids show greater fitness and survival in both virulence and virulence-associated phenotypes than parental lineages under certain conditions. These studies showcase how experimentation in model species such as Cryptococcus can be a powerful tool in elucidating the genotypic and phenotypic consequences of hybridization.