Ploidy dynamics and evolvability in fungi

Chromosome-level subgenome-aware de novo assembly provides insight into Saccharomyces bayanus genome divergence after hybridization

Interspecies hybridization is prevalent in various eukaryotic lineages and plays important roles in phenotypic diversification, adaptation, and speciation. To better understand the changes that occurred in the different subgenomes of a hybrid species and how they facilitate adaptation, we have completed chromosome-level de novo assemblies of all chromosomes for a recently formed hybrid yeast, Saccharomyces bayanus strain CBS380, using Oxford Nanopore Technologies' MinION long-read sequencing. We characterize the S. bayanus genome and compare it with its parent species, Saccharomyces uvarum and Saccharomyces eubayanus, and other S. bayanus genomes to better understand genome evolution after a relatively recent hybridization event. We observe multiple recombination events between the subgenomes in each chromosome, followed by loss of heterozygosity (LOH) in nine chromosome pairs. In addition to maintaining nearly all gene content and synteny from its parental genomes, S. bayanus has acquired many genes from other yeast species, primarily through the introgression of Saccharomyces cerevisiae, such as those involved in the maltose metabolism. Finally, the patterns of recombination and LOH suggest an allotetraploid origin of S. bayanus The gene acquisition and rapid LOH in the hybrid genome probably facilitated its adaptation to maltose brewing environments and mitigated the maladaptive effect of hybridization. This paper describes the first in-depth study using long-read sequencing technology of an S. bayanus hybrid genome, which may serve as an excellent reference for future studies of this important yeast and other yeast strains.

Conserved signalling functions for Mps1, Mad1 and Mad2 in the Cryptococcus neoformans spindle checkpoint

Cryptococcus neoformans is an opportunistic, human fungal pathogen which undergoes fascinating switches in cell cycle control and ploidy when it encounters stressful environments such as the human lung. Here we carry out a mechanistic analysis of the spindle checkpoint which regulates the metaphase to anaphase transition, focusing on Mps1 kinase and the downstream checkpoint components Mad1 and Mad2. We demonstrate that Cryptococcus mad1Δ or mad2Δ strains are unable to respond to microtubule perturbations, continuing to re-bud and divide, and die as a consequence. Fluorescent tagging of Chromosome 3, using a lacO array and mNeonGreen-lacI fusion protein, demonstrates that mad mutants are unable to maintain sister-chromatid cohesion in the absence of microtubule polymers. Thus, the classic checkpoint functions of the SAC are conserved in Cryptococcus. In interphase, GFP-Mad1 is enriched at the nuclear periphery, and it is recruited to unattached kinetochores in mitosis. Purification of GFP-Mad1 followed by mass spectrometric analysis of associated proteins show that it forms a complex with Mad2 and that it interacts with other checkpoint signalling components (Bub1) and effectors (Cdc20 and APC/C sub-units) in mitosis. We also demonstrate that overexpression of Mps1 kinase is sufficient to arrest Cryptococcus cells in mitosis, and show that this arrest is dependent on both Mad1 and Mad2. We find that a C-terminal fragment of Mad1 is an effective in vitro substrate for Mps1 kinase and map several Mad1 phosphorylation sites. Some sites are highly conserved within the C-terminal Mad1 structure and we demonstrate that mutation of threonine 667 (T667A) leads to loss of checkpoint signalling and abrogation of the GAL-MPS1 arrest. Thus Mps1-dependent phosphorylation of C-terminal Mad1 residues is a critical step in Cryptococcus spindle checkpoint signalling. We conclude that CnMps1 protein kinase, Mad1 and Mad2 proteins have all conserved their important, spindle checkpoint signalling roles helping ensure high fidelity chromosome segregation.

Chromosome-level Subgenome-aware de novo Assembly of Saccharomyces bayanus Provides Insight into Genome Divergence after Hybridization

Interspecies hybridization is prevalent in various eukaryotic lineages and plays important roles in phenotypic diversification, adaption, and speciation. To better understand the changes that occurred in the different subgenomes of a hybrid species and how they facilitated adaptation, we completed chromosome-level de novo assemblies of all 16 pairs chromosomes for a recently formed hybrid yeast, Saccharomyces bayanus strain CBS380 (IFO11022), using Nanopore MinION long-read sequencing. Characterization of S. bayanus subgenomes and comparative analysis with the genomes of its parent species, S. uvarum and S. eubayanus, provide several new insights into understanding genome evolution after a relatively recent hybridization. For instance, multiple recombination events between the two subgenomes have been observed in each chromosome, followed by loss of heterozygosity (LOH) in most chromosomes in nine chromosome pairs. In addition to maintaining nearly all gene content and synteny from its parental genomes, S. bayanus has acquired many genes from other yeast species, primarily through the introgression of S. cerevisiae, such as those involved in the maltose metabolism. In addition, the patterns of recombination and LOH suggest an allotetraploid origin of S. bayanus. The gene acquisition and rapid LOH in the hybrid genome probably facilitated its adaption to maltose brewing environments and mitigated the maladaptive effect of hybridization.

Origin of fungal hybrids with pathogenic potential from warm seawater environments

Hybridisation is a common event in yeasts often leading to genomic variability and adaptation. The yeast Candida orthopsilosis is a human-associated opportunistic pathogen belonging to the Candida parapsilosis species complex. Most C. orthopsilosis clinical isolates are hybrids resulting from at least four independent crosses between two parental lineages, of which only one has been identified. The rare presence or total absence of parentals amongst clinical isolates is hypothesised to be a consequence of a reduced pathogenicity with respect to their hybrids. Here, we sequence and analyse the genomes of environmental C. orthopsilosis strains isolated from warm marine ecosystems. We find that a majority of environmental isolates are hybrids, phylogenetically closely related to hybrid clinical isolates. Furthermore, we identify the missing parental lineage, thus providing a more complete overview of the genomic evolution of this species. Additionally, we discover phenotypic differences between the two parental lineages, as well as between parents and hybrids, under conditions relevant for pathogenesis. Our results suggest a marine origin of C. orthopsilosis hybrids, with intrinsic pathogenic potential, and pave the way to identify pre-existing environmental adaptations that rendered hybrids more prone than parental lineages to colonise and infect the mammalian host.

Candida albicans exhibits heterogeneous and adaptive cytoprotective responses to antifungal compounds

Candida albicans, an opportunistic human pathogen, poses a significant threat to human health and is associated with significant socio-economic burden. Current antifungal treatments fail, at least in part, because C. albicans can initiate a strong drug tolerance response that allows some cells to grow at drug concentrations above their minimal inhibitory concentration. To better characterize this cytoprotective tolerance program at the molecular single-cell level, we used a nanoliter droplet-based transcriptomics platform to profile thousands of individual fungal cells and establish their subpopulation characteristics in the absence and presence of antifungal drugs. Profiles of untreated cells exhibit heterogeneous expression that correlates with cell cycle stage with distinct metabolic and stress responses. At 2 days post-fluconazole exposure (a time when tolerance is measurable), surviving cells bifurcate into two major subpopulations: one characterized by the upregulation of genes encoding ribosomal proteins, rRNA processing machinery, and mitochondrial cellular respiration capacity, termed the Ribo-dominant (Rd) state; and the other enriched for genes encoding stress responses and related processes, termed the Stress-dominant (Sd) state. This bifurcation persists at 3 and 6 days post-treatment. We provide evidence that the ribosome assembly stress response (RASTR) is activated in these subpopulations and may facilitate cell survival.

Evolution of loss of heterozygosity patterns in hybrid genomes of Candida yeast pathogens

BACKGROUND

Hybrids are chimeric organisms with highly plastic heterozygous genomes that may confer unique traits enabling the adaptation to new environments. However, most evolutionary theory frameworks predict that the high levels of genetic heterozygosity present in hybrids from divergent parents are likely to result in numerous deleterious epistatic interactions. Under this scenario, selection is expected to favor recombination events resulting in loss of heterozygosity (LOH) affecting genes involved in such negative interactions. Nevertheless, it is so far unknown whether this phenomenon actually drives genomic evolution in natural populations of hybrids. To determine the balance between selection and drift in the evolution of LOH patterns in natural yeast hybrids, we analyzed the genomic sequences from fifty-five hybrid strains of the pathogenic yeasts Candida orthopsilosis and Candida metapsilosis, which derived from at least six distinct natural hybridization events.

RESULTS

We found that, although LOH patterns in independent hybrid clades share some level of convergence that would not be expected from random occurrence, there is an apparent lack of strong functional selection. Moreover, while mitosis is associated with a limited number of inter-homeologous chromosome recombinations in these genomes, induced DNA breaks seem to increase the LOH rate. We also found that LOH does not accumulate linearly with time in these hybrids. Furthermore, some C. orthopsilosis hybrids present LOH patterns compatible with footprints of meiotic recombination. These meiotic-like patterns are at odds with a lack of evidence of sexual recombination and with our inability to experimentally induce sporulation in these hybrids.

CONCLUSIONS

Our results suggest that genetic drift is the prevailing force shaping LOH patterns in these hybrid genomes. Moreover, the observed LOH patterns suggest that these are likely not the result of continuous accumulation of sporadic events-as expected by mitotic repair of rare chromosomal breaks-but rather of acute episodes involving many LOH events in a short period of time.

Aneuploidy and gene dosage regulate filamentation and host colonization by Candida albicans

Aneuploidy is a frequent occurrence in fungal species where it can alter gene expression and promote adaptation to a variety of environmental cues. Multiple forms of aneuploidy have been observed in the opportunistic fungal pathogen Candida albicans, which is a common component of the human gut mycobiome but can escape this niche and cause life-threatening systemic disease. Using a barcode sequencing (Bar-seq) approach, we evaluated a set of diploid C. albicans strains and found that a strain carrying a third copy of chromosome (Chr) 7 was associated with increased fitness during both gastrointestinal (GI) colonization and systemic infection. Our analysis revealed that the presence of a Chr 7 trisomy resulted in decreased filamentation, both in vitro and during GI colonization, relative to isogenic euploid controls. A target gene approach demonstrated that NRG1, encoding a negative regulator of filamentation located on Chr 7, contributes to increased fitness of the aneuploid strain due to inhibition of filamentation in a gene dosage-dependent fashion. Together, these experiments establish how aneuploidy enables the reversible adaptation of C. albicans to its host via gene dosage-dependent regulation of morphology.

Description and Genome Characterization of Three Novel Fungal Strains Isolated from Mars 2020 Mission-Associated Spacecraft Assembly Facility Surfaces-Recommendations for Two New Genera and One Species

National Aeronautics and Space Administration’s (NASA) spacecraft assembly facilities are monitored for the presence of any bacteria or fungi that might conceivably survive a transfer to an extraterrestrial environment. Fungi present a broad and diverse range of phenotypic and functional traits to adapt to extreme conditions, hence the detection of fungi and subsequent eradication of them are needed to prevent forward contamination for future NASA missions. During the construction and assembly for the Mars 2020 mission, three fungal strains with unique morphological and phylogenetic properties were isolated from spacecraft assembly facilities. The reconstruction of phylogenetic trees based on several gene loci (ITS, LSU, SSU, RPB, TUB, TEF1) using multi-locus sequence typing (MLST) and whole genome sequencing (WGS) analyses supported the hypothesis that these were novel species. Here we report the genus or species-level classification of these three novel strains via a polyphasic approach using phylogenetic analysis, colony and cell morphology, and comparative analysis of WGS. The strain FJI-L9-BK-P1 isolated from the Jet Propulsion Laboratory Spacecraft Assembly Facility (JPL-SAF) exhibited a putative phylogenetic relationship with the strain Aaosphaeria arxii CBS175.79 but showed distinct morphology and microscopic features. Another JPL-SAF strain, FJII-L3-CM-DR1, was phylogenetically distinct from members of the family Trichomeriaceae and exhibited morphologically different features from the genera Lithohypha and Strelitziana. The strain FKI-L1-BK-DR1 isolated from the Kennedy Space Center facility was identified as a member of Dothideomycetes incertae sedis and is closely related to the family Kirschsteiniotheliaceae according to a phylogenetic analysis. The polyphasic taxonomic approach supported the recommendation for establishing two novel genera and one novel species. The names Aaosphaeria pasadenensis (FJI-L9-BK-P1 = NRRL 64424 = DSM 114621), Pasadenomyces melaninifex (FJII-L3-CM-DR1 = NRRL 64433 = DSM 114623), and Floridaphiala radiotolerans (FKI-L1-BK-DR1 = NRRL 64434 = DSM 114624) are proposed as type species. Furthermore, resistance to ultraviolet-C and presence of specific biosynthetic gene cluster(s) coding for metabolically active compounds are unique to these strains.

Tetraploidy accelerates adaptation under drug selection in a fungal pathogen

Baseline ploidy significantly impacts evolutionary trajectories and, specifically, tetraploidy is associated with higher rates of adaptation relative to haploidy and diploidy. While the majority of experimental evolution studies investigating ploidy use the budding yeast Saccharomyces cerivisiae, the fungal pathogen Candida albicans is a powerful system to investigate ploidy dynamics, particularly in the context of acquiring antifungal drug resistance. C. albicans laboratory and clinical strains are predominantly diploid, but have been isolated as haploid and polyploid. Here, we evolved diploid and tetraploid C. albicans for ~60 days in the antifungal drug caspofungin. Tetraploid-evolved lines adapted faster than diploid-evolved lines and reached higher levels of caspofungin resistance. While diploid-evolved lines generally maintained their initial genome size, tetraploid-evolved lines rapidly underwent genome-size reductions and did so prior to caspofungin adaptation. While clinical resistance was largely due to mutations in FKS1, these mutations were caused by substitutions in diploid, and indels in tetraploid isolates. Furthermore, fitness costs in the absence of drug selection were significantly less in tetraploid-evolved lines compared to the diploid-evolved lines. Taken together, this work supports a model of adaptation in which the tetraploid state is transient but its ability to rapidly transition ploidy states improves adaptive outcomes and may drive drug resistance in fungal pathogens.

Karyon: a computational framework for the diagnosis of hybrids, aneuploids, and other nonstandard architectures in genome assemblies

BACKGROUND

Recent technological developments have made genome sequencing and assembly highly accessible and widely used. However, the presence in sequenced organisms of certain genomic features such as high heterozygosity, polyploidy, aneuploidy, heterokaryosis, or extreme compositional biases can challenge current standard assembly procedures and result in highly fragmented assemblies. Hence, we hypothesized that genome databases must contain a nonnegligible fraction of low-quality assemblies that result from such type of intrinsic genomic factors.

FINDINGS

Here we present Karyon, a Python-based toolkit that uses raw sequencing data and de novo genome assembly to assess several parameters and generate informative plots to assist in the identification of nonchanonical genomic traits. Karyon includes automated de novo genome assembly and variant calling pipelines. We tested Karyon by diagnosing 35 highly fragmented publicly available assemblies from 19 different Mucorales (Fungi) species.

CONCLUSIONS

Our results show that 10 (28.57%) of the assemblies presented signs of unusual genomic configurations, suggesting that these are common, at least for some lineages within the Fungi.

Genome plasticity driven by aneuploidy and loss of heterozygosity in Trypanosoma cruzi

Trypanosoma cruzi the causative agent of Chagas disease shows a marked genetic diversity and divided into at least six Discrete Typing Units (DTUs). High intra genetic variability has been observed in the TcI DTU, the most widely distributed DTU, where patterns of genomic diversity can provide information on ecological and evolutionary processes driving parasite population structure and genome organization. Chromosomal aneuploidies and rearrangements across multigene families represent an evidence of T. cruzi genome plasticity. We explored genomic diversity among 18 Colombian T. cruzi I clones and 15 T. cruzi I South American strains. Our results confirm high genomic variability, heterozygosity and presence of a clade compatible with the TcI_dom genotype, described for strains from humans in Colombia and Venezuela. TcI showed high structural plasticity across the geographical region studied. Differential events of whole and segmental aneuploidy (SA) along chromosomes even between clones from the same strain were found and corroborated by the depth and allelic frequency. We detected loss of heterozygosity (LOH) events in different chromosomes, however, the size and location of segments under LOH varied between clones. Genes adjacent to breakpoints were evaluated, and retrotransposon hot spot genes flanked the beginning of segmental aneuploidies. Our results suggest that T. cruzi genomes, like those of Leishmania, may have a highly unstable structure and there is now an urgent need to design experiments to explore any potential adaptive role for the plasticity observed.

Genome plasticity in Candida albicans: A cutting-edge strategy for evolution, adaptation, and survival

Candida albicans is the most implicated fungal species that grows as a commensal or opportunistic pathogen in the human host. It is associated with many life-threatening infections, especially in immunocompromised persons. The genome of Candida albicans is very flexible and can withstand a wide assortment of variations in a continuously changing environment. Thus, genome plasticity is central to its adaptation and has long been of considerable interest. C. albicans has a diploid heterozygous genome that is highly dynamic and can display variation from small to large scale chromosomal rearrangement and aneuploidy, which have implications in drug resistance, virulence, and pathogenicity. This review presents an up-to-date overview of recent genomic studies involving C. albicans. It discusses the accumulating evidence that shows how mitotic recombination events, ploidy dynamics, aneuploidy, and loss of heterozygosity (LOH) influence evolution, adaptation, and survival in C. albicans. Understanding the factors that affect the genome is crucial for a proper understanding of species and rapid development and adjustment of therapeutic strategies to mitigate their spread.

Brettanomyces bruxellensis: Overview of the genetic and phenotypic diversity of an anthropized yeast

Human-associated microorganisms are ideal models to study the impact of environmental changes on species evolution and adaptation because of their small genome, short generation time, and their colonization of contrasting and ever-changing ecological niches. The yeast Brettanomyces bruxellensis is a good example of organism facing anthropogenic-driven selective pressures. It is associated with fermentation processes in which it can be considered either as a spoiler (e.g. winemaking, bioethanol production) or as a beneficial microorganism (e.g. production of specific beers, kombucha). Besides its industrial interests, noteworthy parallels and dichotomies with Saccharomyces cerevisiae propelled B. bruxellensis as a valuable complementary yeast model. In this review, we emphasize that the broad genetic and phenotypic diversity of this species is only beginning to be uncovered. Population genomic studies have revealed the co-existence of auto- and allotriploidization events with different evolutionary outcomes. The different diploid, autotriploid and allotriploid subpopulations are associated with specific fermented processes, suggesting independent adaptation events to anthropized environments. Phenotypically, B. bruxellensis is renowned for its ability to metabolize a wide variety of carbon and nitrogen sources, which may explain its ability to colonize already fermented environments showing low-nutrient contents. Several traits of interest could be related to adaptation to human activities (e.g. nitrate metabolization in bioethanol production, resistance to sulphite treatments in winemaking). However, phenotypic traits are insufficiently studied in view of the great genomic diversity of the species. Future work will have to take into account strains of varied substrates, geographical origins as well as displaying different ploidy levels to improve our understanding of an anthropized yeast's phenotypic landscape.

Evaluation of Bacteria and Fungi DNA Abundance in Human Tissues

Whereas targeted and shotgun sequencing approaches are both powerful in allowing the study of tissue-associated microbiota, the human: microorganism abundance ratios in tissues of interest will ultimately determine the most suitable sequencing approach. In addition, it is possible that the knowledge of the relative abundance of bacteria and fungi during a treatment course or in pathological conditions can be relevant in many medical conditions. Here, we present a qPCR-targeted approach to determine the absolute and relative amounts of bacteria and fungi and demonstrate their relative DNA abundance in nine different human tissue types for a total of 87 samples. In these tissues, fungi genomes are more abundant in stool and skin samples but have much lower levels in other tissues. Bacteria genomes prevail in stool, skin, oral swabs, saliva, and gastric fluids. These findings were confirmed by shotgun sequencing for stool and gastric fluids. This approach may contribute to a more comprehensive view of the human microbiota in targeted studies for assessing the abundance levels of microorganisms during disease treatment/progression and to indicate the most informative methods for studying microbial composition (shotgun versus targeted sequencing) for various samples types.

Loss of heterozygosity results in rapid but variable genome homogenization across yeast genetic backgrounds

The dynamics and diversity of the appearance of genetic variants play an essential role in the evolution of the genome and the shaping of biodiversity. Recent population-wide genome sequencing surveys have highlighted the importance of loss of heterozygosity (LOH) events and have shown that they are a neglected part of the genetic diversity landscape. To assess the extent, variability, and spectrum, we explored the accumulation of LOH events in 169 heterozygous diploid Saccharomyces cerevisiae mutation accumulation lines across nine genetic backgrounds. In total, we detected a large set of 22,828 LOH events across distinct genetic backgrounds with a heterozygous level ranging from 0.1% to 1%. LOH events are very frequent with a rate consistently much higher than the mutation rate, showing their importance for genome evolution. We observed that the interstitial LOH (I-LOH) events, resulting in internal short LOH tracts, were much frequent (n = 19,660) than the terminal LOH (T-LOH) events, that is, tracts extending to the end of the chromosome (n = 3168). However, the spectrum, the rate, and the fraction of the genome under LOH vary across genetic backgrounds. Interestingly, we observed that the more the ancestors were heterozygous, the more they accumulated T-LOH events. In addition, frequent short I-LOH tracts are a signature of the lines derived from hybrids with low spore fertility. Finally, we found lines showing almost complete homozygotization during vegetative progression. Overall, our results highlight that the variable dynamics of the LOH accumulation across distinct genetic backgrounds might lead to rapid differential genome evolution during vegetative growth.

Adaptation to the dietary sugar D-tagatose via genome instability in polyploid Candida albicans cells

The opportunistic fungal pathogen Candida albicans undergoes an unusual parasexual cycle wherein diploid cells mate to form tetraploid cells that can generate genetically diverse progeny via a non-meiotic program of chromosome loss. The genetic diversity afforded by parasex impacts clinically relevant features including drug resistance and virulence, and yet the factors influencing genome instability in C. albicans are not well defined. To understand how environmental cues impact genome instability, we monitored ploidy change following tetraploid cell growth in a panel of different carbon sources. We found that growth in one carbon source, D-tagatose, led to high levels of genomic instability and chromosome loss in tetraploid cells. This sugar is a stereoisomer of L-sorbose which was previously shown to promote karyotypic changes in C. albicans. However, while expression of the SOU1 gene enabled utilization of L-sorbose, overexpression of this gene did not promote growth in D-tagatose, indicating differences in assimilation of the two sugars. In addition, genome sequencing of multiple progeny recovered from D-tagatose cultures revealed increased relative copy numbers of chromosome 4, suggestive of chromosome-level regulation of D-tagatose metabolism. Together, these studies identify a novel environmental cue that induces genome instability in C. albicans, and further implicate chromosomal changes in supporting metabolic adaptation in this species.

Genetic Diversity of the Flavohemoprotein Gene of Giardia lamblia: Evidence for High Allelic Heterozygosity and Copy Number Variation

Purpose

The flavohemoprotein (gFlHb) in Giardia plays an important role in managing nitrosative and oxidative stress, and potentially also in virulence and nitroimidazole drug tolerance. The aim of this study was to analyze the genetic diversity of gFlHb in Giardia assemblages A and B clinical isolates.

Methods

gFlHb genes from 20 cultured clinical Giardia isolates were subjected to PCR amplification and cloning, followed by Sanger sequencing. Sequences of all cloned PCR fragments from each isolate were analyzed for single nucleotide variants (SNVs) and compared to genomic Illumina sequence data. Identical clone sequences were sorted into alleles, and diversity was further analyzed. The number of gFlHb gene copies was assessed by mining PacBio de novo assembled genomes in eight isolates. Homology models for assessment of SNV's potential impact on protein function were created using Phyre2.

Results

A variable copy number of the gFlHb gene, between two and six copies, depending on isolate, was found. A total of 37 distinct sequences, representing different alleles of the gFlHb gene, were identified in AII isolates, and 41 were identified in B isolates. In some isolates, up to 12 different alleles were found. The total allelic diversity was high for both assemblages (>0.9) and was coupled with a nucleotide diversity of <0.01. The genetic variation (SNVs per CDS length) was 4.8% in sub-assemblage AII and 5.4% in assemblage B. The number of non-synonymous (ns) SNVs was high in gFIHb of both assemblages, 1.6% in A and 3.0% in B, respectively. Some of the identified nsSNV are predicted to alter protein structure and possibly function.

Conclusion

In this study, we present evidence that gFlHb, a putative protective enzyme against oxidative and nitrosative stress in Giardia, is a variable copy number gene with high allelic diversity. The genetic variability of gFlHb may contribute metabolic adaptability against metronidazole toxicity.

Genomic Diversification, Structural Plasticity, and Hybridization in Leishmania (Viannia) braziliensis

Leishmania (Viannia) braziliensis is an important Leishmania species circulating in several Central and South American countries. Among Leishmania species circulating in Brazil, Argentina and Colombia, L. braziliensis has the highest genomic variability. However, genomic variability at the whole genome level has been only studied in Brazilian and Peruvian isolates; to date, no Colombian isolates have been studied. Considering that in Colombia, L. braziliensis is a species with great clinical and therapeutic relevance, as well as the role of genetic variability in the epidemiology of leishmaniasis, we analyzed and evaluated intraspecific genomic variability of L. braziliensis from Colombian and Bolivian isolates and compared them with Brazilian isolates. Twenty-one genomes were analyzed, six from Colombian patients, one from a Bolivian patient, and 14 Brazilian isolates downloaded from public databases. The results obtained of Phylogenomic analysis showed the existence of four well-supported clades, which evidenced intraspecific variability. The whole-genome analysis revealed structural variations in the somy, mainly in the Brazilian genomes (clade 1 and clade 3), low copy number variations, and a moderate number of single-nucleotide polymorphisms (SNPs) in all genomes analyzed. Interestingly, the genomes belonging to clades 2 and 3 from Colombia and Brazil, respectively, were characterized by low heterozygosity (~90% of SNP loci were homozygous) and regions suggestive of loss of heterozygosity (LOH). Additionally, we observed the drastic whole genome loss of heterozygosity and possible hybridization events in one genome belonging to clade 4. Unique/shared SNPs between and within the four clades were identified, revealing the importance of some of them in biological processes of L. braziliensis. Our analyses demonstrate high genomic variability of L. braziliensis in different regions of South America, mainly in Colombia and suggest that this species exhibits striking genomic diversity and a capacity of genomic hybridization; additionally, this is the first study to report whole-genome sequences of Colombian L. braziliensis isolates.

What do we know about the biology of the emerging fungal pathogen of humans Candida auris?

Candida auris is a worrisome fungal pathogen of humans which emerged merely about a decade ago. Ever since then the scientific community worked hard to understand clinically relevant traits, such as virulence factors, antifungal resistance mechanisms, and its ability to adhere to human skin and medical devices. Whole-genome sequencing of clinical isolates and epidemiological studies outlining the path of nosocomial outbreaks have been the focus of research into this pathogenic and multidrug-resistant yeast since its first description in 2009. More recently, work was started by several laboratories to explore the biology of C. auris. Here, we review the insights of studies characterizing the mechanisms underpinning antifungal drug resistance, biofilm formation, morphogenetic switching, cell aggregation, virulence, and pathogenicity of C. auris. We conclude that, although some progress has been made, there is still a long journey ahead of us, before we fully understand this novel pathogen. Critically important is the development of molecular tools for C. auris to make this fungus genetically tractable and traceable. This will allow an in-depth molecular dissection of the life cycle of C. auris, of its characteristics while interacting with the human host, and the mechanisms it employs to avoid being killed by antifungals and the immune system.

The genetic basis of aneuploidy tolerance in wild yeast

Aneuploidy is highly detrimental during development yet common in cancers and pathogenic fungi – what gives rise to differences in aneuploidy tolerance remains unclear. We previously showed that wild isolates of Saccharomyces cerevisiae tolerate chromosome amplification while laboratory strains used as a model for aneuploid syndromes do not. Here, we mapped the genetic basis to Ssd1, an RNA-binding translational regulator that is functional in wild aneuploids but defective in laboratory strain W303. Loss of SSD1 recapitulates myriad aneuploidy signatures previously taken as eukaryotic responses. We show that aneuploidy tolerance is enabled via a role for Ssd1 in mitochondrial physiology, including binding and regulating nuclear-encoded mitochondrial mRNAs, coupled with a role in mitigating proteostasis stress. Recapitulating ssd1Δ defects with combinatorial drug treatment selectively blocked proliferation of wild-type aneuploids compared to euploids. Our work adds to elegant studies in the sensitized laboratory strain to present a mechanistic understanding of eukaryotic aneuploidy tolerance.

A Double-Edged Sword: Aneuploidy is a Prevalent Strategy in Fungal Adaptation

Aneuploidy, a deviation from a balanced genome by either gain or loss of chromosomes, is generally associated with impaired fitness and developmental defects in eukaryotic organisms. While the general physiological impact of aneuploidy remains largely elusive, many phenotypes associated with aneuploidy link to a common theme of stress adaptation. Here, we review previously identified mechanisms and observations related to aneuploidy, focusing on the highly diverse eukaryotes, fungi. Fungi, which have conquered virtually all environments, including several hostile ecological niches, exhibit widespread aneuploidy and employ it as an adaptive strategy under severe stress. Gambling with the balance between genome plasticity and stability has its cost and in fact, most aneuploidies have fitness defects. How can this fitness defect be reconciled with the prevalence of aneuploidy in fungi? It is likely that the fitness cost of the extra chromosomes is outweighed by the advantage they confer under life-threatening stresses. In fact, once the selective pressures are withdrawn, aneuploidy is often lost and replaced by less drastic mutations that possibly incur a lower fitness cost. We discuss representative examples across hostile environments, including medically and industrially relevant cases, to highlight potential adaptive mechanisms in aneuploid yeast.

Rapid and extensive karyotype diversification in haploid clinical Candida auris isolates

Candida auris is a newly emerged pathogenic microbe, having been identified as a medically relevant fungus as recently as 2009. It is one of the most drug-resistant yeast species known to date and its emergence and population structure are unusual. Because of its recent emergence, we are largely ignorant about fundamental aspects of its general biology, life cycle, and population dynamics. Here, we report the karyotype variability of 26 C. auris strains representing the four main clades. We demonstrate that all strains are haploid and have a highly plastic karyotype containing five to seven chromosomes, which can undergo marked alterations within a short time frame when the fungus is put under genotoxic, heat, or osmotic stress. No simple correlation was found between karyotype pattern, drug resistance, and clade affiliation indicating that karyotype heterogeneity is rapidly evolving. As with other Candida species, these marked karyotype differences between isolates are likely to have an important impact on pathogenic traits of C. auris.

Adaptation by Loss of Heterozygosity in Saccharomyces cerevisiae Clones Under Divergent Selection

Loss of heterozygosity (LOH) is observed during vegetative growth and reproduction of diploid genotypes through mitotic crossovers, aneuploidy caused by nondisjunction, and gene conversion. We aimed to test the role that LOH plays during adaptation of two highly heterozygous Saccharomyces cerevisiae genotypes to multiple environments over a short time span in the laboratory. We hypothesized that adaptation would be observed through parallel LOH events across replicate populations. Using genome resequencing of 70 clones, we found that LOH was widespread with 5.2 LOH events per clone after ∼500 generations. The most common mode of LOH was gene conversion (51%) followed by crossing over consistent with either break-induced replication or double Holliday junction resolution. There was no evidence that LOH involved nondisjunction of whole chromosomes. We observed parallel LOH in both an environment-specific and environment-independent manner. LOH largely involved recombining existing variation between the parental genotypes, but also was observed after de novo, presumably beneficial, mutations occurred in the presence of canavanine, a toxic analog of arginine. One highly parallel LOH event involved the ENA salt efflux pump locus on chromosome IV, which showed repeated LOH to the allele from the European parent, an allele originally derived by introgression from S. paradoxus Using CRISPR-engineered LOH we showed that the fitness advantage provided by this single LOH event was 27%. Overall, we found extensive evidence that LOH could be adaptive and is likely to be a greater source of initial variation than de novo mutation for rapid evolution of diploid genotypes.

Seasons of change: Mechanisms of genome evolution in human fungal pathogens

Fungi are a diverse kingdom of organisms capable of thriving in various niches across the world including those in close association with multicellular eukaryotes. Fungal pathogens that contribute to human disease reside both within the host as commensal organisms of the microbiota and the environment. Their niche of origin dictates how infection initiates but also places specific selective pressures on the fungal pathogen that contributes to its genome organization and genetic repertoire. Recent efforts to catalogue genomic variation among major human fungal pathogens have unveiled evolutionary themes that shape the fungal genome. Mechanisms ranging from large scale changes such as aneuploidy and ploidy cycling as well as more targeted mutations like base substitutions and gene copy number variations contribute to the evolution of these species, which are often under multiple competing selective pressures with their host, environment, and other microbes. Here, we provide an overview of the major selective pressures and mechanisms acting to evolve the genome of clinically important fungal pathogens of humans.

Proteins that physically interact with the phosphatase Cdc14 in Candida albicans have diverse roles in the cell cycle

The chromosome complement of the human fungal pathogen Candida albicans is unusually unstable, suggesting that the process of nuclear division is error prone. The Cdc14 phosphatase plays a key role in organising the intricate choreography of mitosis and cell division. In order to understand the role of Cdc14 in C. albicans we used quantitative proteomics to identify proteins that physically interact with Cdc14. To distinguish genuine Cdc14-interactors from proteins that bound non-specifically to the affinity matrix, we used a substrate trapping mutant combined with mass spectrometry analysis using Stable Isotope Labelling with Amino Acids in Cell Culture (SILAC). The results identified 126 proteins that interact with Cdc14 of which 80% have not previously been identified as Cdc14 interactors in C. albicans or S. cerevisiae. In this set, 55 proteins are known from previous research in S. cerevisiae and S. pombe to play roles in the cell cycle, regulating the attachment of the mitotic spindle to kinetochores, mitotic exit, cytokinesis, licensing of DNA replication by re-activating pre-replication complexes, and DNA repair. Five Cdc14-interacting proteins with previously unknown functions localised to the Spindle Pole Bodies (SPBs). Thus, we have greatly increased the number of proteins that physically interact with Cdc14 in C. albicans.

Dynamic ploidy changes drive fluconazole resistance in human cryptococcal meningitis

BACKGROUND

Cryptococcal meningitis (CM) causes an estimated 180,000 deaths annually, predominantly in sub-Saharan Africa, where most patients receive fluconazole (FLC) monotherapy. While relapse after FLC monotherapy with resistant strains is frequently observed, the mechanisms and impact of emergence of FLC resistance in human CM are poorly understood. Heteroresistance (HetR) - a resistant subpopulation within a susceptible strain - is a recently described phenomenon in Cryptococcus neoformans (Cn) and Cryptococcus gattii (Cg), the significance of which has not previously been studied in humans.

METHODS

A cohort of 20 patients with HIV-associated CM in Tanzania was prospectively observed during therapy with either FLC monotherapy or in combination with flucytosine (5FC). Total and resistant subpopulations of Cryptococcus spp. were quantified directly from patient cerebrospinal fluid (CSF). Stored isolates underwent whole genome sequencing and phenotypic characterization.

RESULTS

Heteroresistance was detectable in Cryptococcus spp. in the CSF of all patients at baseline (i.e., prior to initiation of therapy). During FLC monotherapy, the proportion of resistant colonies in the CSF increased during the first 2 weeks of treatment. In contrast, no resistant subpopulation was detectable in CSF by day 14 in those receiving a combination of FLC and 5FC. Genomic analysis revealed high rates of aneuploidy in heteroresistant colonies as well as in relapse isolates, with chromosome 1 (Chr1) disomy predominating. This is apparently due to the presence on Chr1 of ERG11, which is the FLC drug target, and AFR1, which encodes a drug efflux pump. In vitro efflux levels positively correlated with the level of heteroresistance.

CONCLUSION

Our findings demonstrate for what we believe is the first time the presence and emergence of aneuploidy-driven FLC heteroresistance in human CM, association of efflux levels with heteroresistance, and the successful suppression of heteroresistance with 5FC/FLC combination therapy.

FUNDING

This work was supported by the Wellcome Trust Strategic Award for Medical Mycology and Fungal Immunology 097377/Z/11/Z and the Daniel Turnberg Travel Fellowship.

FTIR and Raman Spectroscopy-Based Biochemical Profiling Reflects Genomic Diversity of Clinical Candida Isolates That May Be Useful for Diagnosis and Targeted Therapy of Candidiasis

Despite the fact that Candida albicans is documented to be the main cause of human candidiasis, non-C. albicans Candida (NCAC) species, such as Candida glabrata and Candida tropicalis, are also suggested to be implicated in the etiopathogenesis of opportunistic fungal infections. As biology, epidemiology, pathogenicity, and antifungal resistance of NCAC species may be affected as a result of genomic diversity and plasticity, rapid and unambiguous identification of Candida species in clinical samples is essential for proper diagnosis and therapy. In the present study, 25 clinical isolates of C. albicans, C. glabrata, and C. tropicalis species were characterized in terms of their karyotype patterns, DNA content, and biochemical features. Fourier transform infrared (FTIR) spectra- and Raman spectra-based molecular fingerprints corresponded to the diversity of chromosomal traits and DNA levels that provided correct species identification. Moreover, Raman spectroscopy was documented to be useful for the evaluation of ergosterol content that may be associated with azole resistance. Taken together, we found that vibrational spectroscopy-based biochemical profiling reflects the variability of chromosome patterns and DNA content of clinical Candida species isolates and may facilitate the diagnosis and targeted therapy of candidiasis.

Competition experiments between Brettanomyces bruxellensis strains reveal specific adaptation to sulfur dioxide and complex interactions at intraspecies level

Recent studies have suggested a strong niche adaptation for Brettanomyces bruxellensis strains according to human-related fermentation environments, including beer, wine and bioethanol. This is further supported by a correlation between B. bruxellensis genetic grouping and tolerance to SO2, the main antimicrobial used in wine. The allotriploid AWRI1499-like cluster, in particular, shows high SO2 tolerance suggesting that the genetic configuration observed for these strains may confer a selective advantage in winemaking conditions. To test this hypothesis, we evaluated the relative selective advantage of representatives of the three main B. bruxellensis genetic groups in presence of SO2. As a proof-of-concept and using recently developed transformation cassettes, we compared strains under different SO2 concentrations using pairwise competitive fitness experiments. Our results showed that AWRI1499 is specifically adapted to environments with high SO2 concentrations compared to other B. bruxellensis wine strains, indicating a potential correlation between allotriploidisation origin and environmental adaptation in this species. Additionally, our findings suggest different types of competition between strains, such as coexistence and exclusion, revealing new insights on B. bruxellensis interactions at intraspecies level.

Methodologies for in vitro and in vivo evaluation of efficacy of antifungal and antibiofilm agents and surface coatings against fungal biofilms

Unlike superficial fungal infections of the skin and nails, which are the most common fungal diseases in humans, invasive fungal infections carry high morbidity and mortality, particularly those associated with biofilm formation on indwelling medical devices. Therapeutic management of these complex diseases is often complicated by the rise in resistance to the commonly used antifungal agents. Therefore, the availability of accurate susceptibility testing methods for determining antifungal resistance, as well as discovery of novel antifungal and antibiofilm agents, are key priorities in medical mycology research. To direct advancements in this field, here we present an overview of the methods currently available for determining (i) the susceptibility or resistance of fungal isolates or biofilms to antifungal or antibiofilm compounds and compound combinations; (ii) the in vivo efficacy of antifungal and antibiofilm compounds and compound combinations; and (iii) the in vitro and in vivo performance of anti-infective coatings and materials to prevent fungal biofilm-based infections.

Characterization of phenotypic variation and genome aberrations observed among Phytophthora ramorum isolates from diverse hosts

BACKGROUND

Accumulating evidence suggests that genome plasticity allows filamentous plant pathogens to adapt to changing environments. Recently, the generalist plant pathogen Phytophthora ramorum has been documented to undergo irreversible phenotypic alterations accompanied by chromosomal aberrations when infecting trunks of mature oak trees (genus Quercus). In contrast, genomes and phenotypes of the pathogen derived from the foliage of California bay (Umbellularia californica) are usually stable. We define this phenomenon as host-induced phenotypic diversification (HIPD). P. ramorum also causes a severe foliar blight in some ornamental plants such as Rhododendron spp. and Viburnum spp., and isolates from these hosts occasionally show phenotypes resembling those from oak trunks that carry chromosomal aberrations. The aim of this study was to investigate variations in phenotypes and genomes of P. ramorum isolates from non-oak hosts and substrates to determine whether HIPD changes may be equivalent to those among isolates from oaks.

RESULTS

We analyzed genomes of diverse non-oak isolates including those taken from foliage of Rhododendron and other ornamental plants, as well as from natural host species, soil, and water. Isolates recovered from artificially inoculated oak logs were also examined. We identified diverse chromosomal aberrations including copy neutral loss of heterozygosity (cnLOH) and aneuploidy in isolates from non-oak hosts. Most identified aberrations in non-oak hosts were also common among oak isolates; however, trisomy, a frequent type of chromosomal aberration in oak isolates was not observed in isolates from Rhododendron.

CONCLUSION

This work cross-examined phenotypic variation and chromosomal aberrations in P. ramorum isolates from oak and non-oak hosts and substrates. The results suggest that HIPD comparable to that occurring in oak hosts occurs in non-oak environments such as in Rhododendron leaves. Rhododendron leaves are more easily available than mature oak stems and thus can potentially serve as a model host for the investigation of HIPD, the newly described plant-pathogen interaction.

Brettanomyces bruxellensis population survey reveals a diploid-triploid complex structured according to substrate of isolation and geographical distribution

Brettanomyces bruxellensis is a unicellular fungus of increasing industrial and scientific interest over the past 15 years. Previous studies revealed high genotypic diversity amongst B. bruxellensis strains as well as strain-dependent phenotypic characteristics. Genomic assemblies revealed that some strains harbour triploid genomes and based upon prior genotyping it was inferred that a triploid population was widely dispersed across Australian wine regions. We performed an intraspecific diversity genotypic survey of 1488 B. bruxellensis isolates from 29 countries, 5 continents and 9 different fermentation niches. Using microsatellite analysis in combination with different statistical approaches, we demonstrate that the studied population is structured according to ploidy level, substrate of isolation and geographical origin of the strains, underlying the relative importance of each factor. We found that geographical origin has a different contribution to the population structure according to the substrate of origin, suggesting an anthropic influence on the spatial biodiversity of this microorganism of industrial interest. The observed clustering was correlated to variable stress response, as strains from different groups displayed variation in tolerance to the wine preservative sulfur dioxide (SO₂). The potential contribution of the triploid state for adaptation to industrial fermentations and dissemination of the species B. bruxellensis is discussed.

Tackling emerging fungal threats to animal health, food security and ecosystem resilience

Emerging infections caused by fungi have become a widely recognized global phenomenon. Their notoriety stems from their causing plagues and famines, driving species extinctions, and the difficulty in treating human mycoses alongside the increase of their resistance to antifungal drugs. This special issue comprises a collection of articles resulting from a Royal Society discussion meeting examining why pathogenic fungi are causing more disease now than they did in the past, and how we can tackle this rapidly emerging threat to the health of plants and animals worldwide.This article is part of the themed issue 'Tackling emerging fungal threats to animal health, food security and ecosystem resilience'.

Hybridization and emergence of virulence in opportunistic human yeast pathogens

Hybridization between different species can result in the emergence of new lineages and adaptive phenotypes. Occasionally, hybridization in fungal organisms can drive the appearance of opportunistic lifestyles or shifts to new hosts, resulting in the emergence of novel pathogens. In recent years, an increasing number of studies have documented the existence of hybrids in diverse yeast clades, including some comprising human pathogens. Comparative and population genomics studies performed on these clades are enabling us to understand what roles hybridization may play in the evolution and emergence of a virulence potential towards humans. Here we survey recent genomic studies on several yeast pathogenic clades where hybrids have been identified, and discuss the broader implications of hybridization in the evolution and emergence of pathogenic lineages. © 2017 The Authors. Yeast published by John Wiley & Sons, Ltd.

Parasex Generates Phenotypic Diversity de Novo and Impacts Drug Resistance and Virulence in Candida albicans

Candida albicans is a diploid fungus that is a frequent cause of mucosal and systemic infections in humans. This species exhibits an unusual parasexual cycle in which mating produces tetraploid cells that undergo a nonmeiotic program of concerted chromosome loss to return to a diploid or aneuploid state. In this work, we used a multipronged approach to examine the capacity of parasex to generate diversity in C. albicans First, we compared the phenotypic properties of 32 genotyped progeny and observed wide-ranging differences in fitness, filamentation, biofilm formation, and virulence. Strikingly, one parasexual isolate displayed increased virulence relative to parental strains using a Galleria mellonella model of infection, establishing that parasex has the potential to enhance pathogenic traits. Next, we examined parasexual progeny derived from homothallic, same-sex mating events, and reveal that parasex can generate diversity de novo from identical parental strains. Finally, we generated pools of parasexual progeny and examined resistance of these pools to environmental stresses. Parasexual progeny were generally less fit than control strains across most test conditions, but showed an increased ability to grow in the presence of the antifungal drug fluconazole (FL). FL-resistant progeny were aneuploid isolates, often being diploid strains trisomic for both Chr3 and Chr6. Passaging of these aneuploid strains frequently led to loss of the supernumerary chromosomes and a concomitant decrease in drug resistance. These experiments establish that parasex generates extensive phenotypic diversity de novo, and that this process has important consequences for both virulence and drug resistance in C. albicans populations.

Microbe Profile: Candida albicans: a shape-changing, opportunistic pathogenic fungus of humans

Candida albicans is normally a harmless commensal of human beings, but it can cause superficial infections of the mucosa (oral/vaginal thrush) in healthy individuals and (rarely) infections of the skin or nails. It can also become invasive, causing life-threatening systemic and bloodstream infections in immunocompromised hosts, where the mortality rate can be as high as 50 %. It is the most common cause of serious fungal infection and is a common cause of nosocomial infections in hospitals. Some strains have been recognized that are resistant to azoles or echinocandins, which are the first-line antifungals for treatment of C. albicans infections.

Tolerance to Caspofungin in Candida albicans Is Associated with at Least Three Distinctive Mechanisms That Govern Expression of FKS Genes and Cell Wall Remodeling

Expanding echinocandin use to prevent or treat invasive fungal infections has led to an increase in the number of breakthrough infections due to resistant Candida species. Although it is uncommon, echinocandin resistance is well documented for Candida albicans, which is among the most prevalent bloodstream organisms. A better understanding is needed to assess the cellular factors that promote tolerance and predispose infecting cells to clinical breakthrough. We previously showed that some mutants that were adapted to growth in the presence of toxic sorbose due to loss of one chromosome 5 (Ch5) also became more tolerant to caspofungin. We found here, following direct selection of mutants on caspofungin, that tolerance can be conferred by at least three mechanisms: (i) monosomy of Ch5, (ii) combined monosomy of the left arm and trisomy of the right arm of Ch5, and (iii) an aneuploidy-independent mechanism. Tolerant mutants possessed cell walls with elevated chitin and showed downregulation of genes involved in cell wall biosynthesis, namely, FKS, located outside Ch5, and CHT2, located on Ch5, irrespective of Ch5 ploidy. Also irrespective of Ch5 ploidy, the CNB1 and MID1 genes on Ch5, which are involved in the calcineurin signaling pathway, were expressed at the diploid level. Thus, multiple mechanisms can affect the relative expression of the aforementioned genes, controlling them in similar ways. Although breakthrough mutations in two specific regions of FKS1 have previously been associated with caspofungin resistance, we found mechanisms of caspofungin tolerance that are independent of FKS1 and thus represent an earlier event in resistance development.