Search

The impacts of microsporidia on host individuals are frequently subtle and can be context dependent. A key example of the latter comes from a recently discovered microsporidian symbiont of Daphnia, the net impact of which was found to shift from negative to positive based on environmental context. Given this, we hypothesized low baseline virulence of the microsporidian; here, we investigated the impact of infection on hosts in controlled conditions and the absence of other stressors. We also investigated its phylogenetic position, ecology, and host range. The genetic data indicate that the symbiont is Ordospora pajunii, a newly described microsporidian parasite of Daphnia. We show that O. pajunii infection damages the gut, causing infected epithelial cells to lose microvilli and then rupture. The prevalence of this microsporidian could be high (up to 100% in the lab and 77% of adults in the field). Its overall virulence was low in most cases, but some genotypes suffered reduced survival and/or reproduction. Susceptibility and virulence were strongly host-genotype dependent. We found that North American O. pajunii were able to infect multiple Daphnia species, including the European species Daphnia longispina, as well as Ceriodaphnia spp. Given the low, often undetectable virulence of this microsporidian and potentially far-reaching consequences of infections for the host when interacting with other pathogens or food, this Daphnia-O. pajunii symbiosis emerges as a valuable system for studying the mechanisms of context-dependent shifts between mutualism and parasitism, as well as for understanding how symbionts might alter host interactions with resources.ImportanceThe net outcome of symbiosis depends on the costs and benefits to each partner. Those can be context dependent, driving the potential for an interaction to change between parasitism and mutualism. Understanding the baseline fitness impact in an interaction can help us understand those shifts; for an organism that is generally parasitic, it should be easier for it to become a mutualist if its baseline virulence is relatively low. Recently, a microsporidian was found to become beneficial to its Daphnia hosts in certain ecological contexts, but little was known about the symbiont (including its species identity). Here, we identify it as the microsporidium Ordospora pajunii. Despite the parasitic nature of microsporidia, we found O. pajunii to be, at most, mildly virulent; this helps explain why it can shift toward mutualism in certain ecological contexts and helps establish O. pajunii is a valuable model for investigating shifts along the mutualism-parasitism continuum.

The global panzootic lineage (GPL) of the pathogenic fungus Batrachochytrium dendrobatidis (Bd) has caused severe amphibian population declines, yet the drivers underlying the high frequency of GPL in regions of amphibian decline are unclear. Using publicly available Bd genome sequences, we identified multiple non-GPL Bd isolates that contain a circular Rep-encoding single-stranded (CRESS)-like DNA virus, which we named Bd DNA virus 1 (BdDV-1). We further sequenced and constructed genome assemblies with long read sequences to find that the virus is integrated into the nuclear genome in some strains. Attempts to cure virus-positive isolates were unsuccessful; however, phenotypic differences between naturally virus-positive and virus-negative Bd isolates suggested that BdDV-1 decreases the growth of its host in vitro but increases the virulence of its host in vivo. BdDV-1 is the first-described CRESS DNA mycovirus of zoosporic true fungi, with a distribution inversely associated with the emergence of the panzootic lineage.

Giant viruses of the phylum Nucleocytoviricota have captured researchers’ attention due to their increasingly recognized impacts on eukaryotic genome evolution. Their origins are hypothesized to predate or coincide with the diversification of eukaryotes and their natural hosts span the eukaryotic tree of life. But surprisingly, giant viruses have not been found in the fungal kingdom, though recent metagenomic work suggests a putative association. Here we show that early-diverging fungal lineages maintain both “viral fossils” and host active infections by giant viruses, which form a monophyletic clade that we name Mycodnaviridae . Viral genomes described here span up to 350kb and encode over 300 genes, including many genes characteristic of giant viruses. Interestingly, we observed variation in infection status among the isolates including apparent active infection and transcriptionally-suppressed states, suggestive that viral activation may be constrained to certain life stages of the host. Our experimental findings add to the scant few natural virus-host systems available in culture for the study of giant viruses. These viruses have likely shaped the early evolution of these fungal lineages and should prove a useful model for their study.

The Global Panzootic Lineage (GPL) of the pathogenic fungus Batrachochytrium dendrobatidis ( Bd ) has caused severe amphibian population declines, yet the drivers underlying the high frequency of GPL in regions of amphibian decline are unclear. Using publicly available Bd genome sequences, we identified multiple non-GPL Bd isolates that contain a circular Rep-encoding single stranded DNA-like virus which we named BdDV-1. We further sequenced and constructed genome assemblies with long read sequences to find that the virus is likely integrated into the nuclear genome. Attempts to cure virus positive isolates were unsuccessful, however, phenotypic differences between naturally virus positive and virus negative Bd isolates suggested that BdDV-1 decreases the growth of its host in vitro but increases the virulence of its host in vivo . BdDV-1 is the first described CRESS DNA mycovirus to increase host virulence with a distribution inversely associated with the emergence of the panzootic lineage.

The first genome sequenced of a eukaryotic organism was for Saccharomyces cerevisiae, as reported in 1996, but it was more than 10 years before any of the zygomycete fungi, which are the early-diverging terrestrial fungi currently placed in the phyla Mucoromycota and Zoopagomycota, were sequenced. The genome for Rhizopus delemar was completed in 2008; currently, more than 1000 zygomycete genomes have been sequenced. Genomic data from these early-diverging terrestrial fungi revealed deep phylogenetic separation of the two major clades-primarily plant-associated saprotrophic and mycorrhizal Mucoromycota versus the primarily mycoparasitic or animal-associated parasites and commensals in the Zoopagomycota. Genomic studies provide many valuable insights into how these fungi evolved in response to the challenges of living on land, including adaptations to sensing light and gravity, development of hyphal growth, and co-existence with the first terrestrial plants. Genome sequence data have facilitated studies of genome architecture, including a history of genome duplications and horizontal gene transfer events, distribution and organization of mating type loci, rDNA genes and transposable elements, methylation processes, and genes useful for various industrial applications. Pathogenicity genes and specialized secondary metabolites have also been detected in soil saprobes and pathogenic fungi. Novel endosymbiotic bacteria and viruses have been discovered during several zygomycete genome projects. Overall, genomic information has helped to resolve a plethora of research questions, from the placement of zygomycetes on the evolutionary tree of life and in natural ecosystems, to the applied biotechnological and medical questions.

The fungal kingdom represents an extraordinary diversity of organisms with profound impacts across animal, plant, and ecosystem health. Fungi simultaneously support life, by forming beneficial symbioses with plants and producing life-saving medicines, and bring death, by causing devastating diseases in humans, plants, and animals. With climate change, increased antimicrobial resistance, global trade, environmental degradation, and novel viruses altering the impact of fungi on health and disease, developing new approaches is now more crucial than ever to combat the threats posed by fungi and to harness their extraordinary potential for applications in human health, food supply, and environmental remediation. To address this aim, the Canadian Institute for Advanced Research (CIFAR) and the Burroughs Wellcome Fund convened a workshop to unite leading experts on fungal biology from academia and industry to strategize innovative solutions to global challenges and fungal threats. This report provides recommendations to accelerate fungal research and highlights the major research advances and ideas discussed at the meeting pertaining to 5 major topics: (1) Connections between fungi and climate change and ways to avert climate catastrophe; (2) Fungal threats to humans and ways to mitigate them; (3) Fungal threats to agriculture and food security and approaches to ensure a robust global food supply; (4) Fungal threats to animals and approaches to avoid species collapse and extinction; and (5) Opportunities presented by the fungal kingdom, including novel medicines and enzymes.

Most of the described species in kingdom Fungi are contained in two phyla, the Ascomycota and the Basidiomycota (subkingdom Dikarya). As a result, our understanding of the biology of the kingdom is heavily influenced by traits observed in Dikarya, such as aerial spore dispersal and life cycles dominated by mitosis of haploid nuclei. We now appreciate that Fungi comprises numerous phylum-level lineages in addition to those of Dikarya, but the phylogeny and genetic characteristics of most of these lineages are poorly understood due to limited genome sampling. Here, we addressed major evolutionary trends in the non-Dikarya fungi by phylogenomic analysis of 69 newly generated draft genome sequences of the zoosporic (flagellated) lineages of true fungi. Our phylogeny indicated five lineages of zoosporic fungi and placed Blastocladiomycota, which has an alternation of haploid and diploid generations, as branching closer to the Dikarya than to the Chytridiomyceta. Our estimates of heterozygosity based on genome sequence data indicate that the zoosporic lineages plus the Zoopagomycota are frequently characterized by diploid-dominant life cycles. We mapped additional traits, such as ancestral cell-cycle regulators, cell-membrane- and cell-wall-associated genes, and the use of the amino acid selenocysteine on the phylogeny and found that these ancestral traits that are shared with Metazoa have been subject to extensive parallel loss across zoosporic lineages. Together, our results indicate a gradual transition in the genetics and cell biology of fungi from their ancestor and caution against assuming that traits measured in Dikarya are typical of other fungal lineages.

We describe environmental microbial eukaryotes (EMEs) sample collection, single-cell isolation, lysis, and genome amplification, followed by the rDNA amplification and OTU screening for recovery of high-quality species-specific genomes for de novo assembly. These protocols are part of our pipeline that also includes bioinformatic methods. The pipeline and its application on a wide range of phyla of different sample complexity are described in our complementary paper. In addition, this protocol describes optimized lysis, genome amplification, and OTU screening steps of the pipeline. For complete details on the use and execution of this protocol, please refer to Ciobanu et al. (2021).

Plant diversity varies immensely over large-scale gradients in temperature, precipitation, and seasonality at global and regional scales. This relationship may be driven in part by climatic variation in the relative importance of abiotic and biotic interactions to the diversity and composition of plant communities. In particular, biotic interactions may become stronger and more host specific with increasing precipitation and temperature, resulting in greater plant species richness in wetter and warmer environments. This hypothesis predicts that the many defensive compounds found in plants’ metabolomes should increase in richness and decrease in interspecific similarity with precipitation, temperature, and plant diversity. To test this prediction, we compared patterns of chemical and morphological trait diversity of 140 woody plant species among seven temperate forests in North America representing 16.2°C variation in mean annual temperature (MAT), 2,115 mm variation in mean annual precipitation (MAP), and from 10 to 68 co-occurring species. We used untargeted metabolomics methods based on data generated with liquid chromatography-tandem mass spectrometry to identify, classify, and compare 13,480 unique foliar metabolites and to quantify the metabolomic similarity of species in each community with respect to the whole metabolome and each of five broad classes of metabolites. In addition, we compiled morphological trait data from existing databases and field surveys for three commonly measured traits (specific leaf area [SLA], wood density, and seed mass) for comparison with foliar metabolomes. We found that chemical defense strategies and growth and allocation strategies reflected by these traits largely represented orthogonal axes of variation. In addition, functional dispersion of SLA increased with MAP, whereas functional richness of wood density and seed mass increased with MAT. In contrast, chemical similarity of co-occurring species decreased with both MAT and MAP, and metabolite richness increased with MAT. Variation in metabolite richness among communities was positively correlated with species richness, but variation in mean chemical similarity was not. Our results are consistent with the hypothesis that plant metabolomes play a more important role in community assembly in wetter and warmer climates, even at temperate latitudes, and suggest that metabolomic traits can provide unique insight to studies of trait-based community assembly.

Single-cell sequencing of environmental microorganisms is an essential component of the microbial ecology toolkit. However, large-scale targeted single-cell sequencing for the whole-genome recovery of uncultivated eukaryotes is lagging. The key challenges are low abundance in environmental communities, large complex genomes, and cell walls that are difficult to break. We describe a pipeline composed of state-of-the art single-cell genomics tools and protocols optimized for poorly studied and uncultivated eukaryotic microorganisms that are found at low abundance. This pipeline consists of seven distinct steps, beginning with sample collection and ending with genome annotation, each equipped with quality review steps to ensure high genome quality at low cost. We tested and evaluated each step on environmental samples and cultures of early-diverging lineages of fungi and Chromista/SAR. We show that genomes produced using this pipeline are almost as good as complete reference genomes for functional and comparative genomics for environmental microbial eukaryotes.

Phylogenomic studies using genome-scale amounts of data have greatly improved understanding of the tree of life. Despite the diversity, ecological significance, and biomedical and industrial importance of fungi, evolutionary relationships among several major lineages remain poorly resolved, especially those near the base of the fungal phylogeny. To examine poorly resolved relationships and assess progress toward a genome-scale phylogeny of the fungal kingdom, we compiled a phylogenomic data matrix of 290 genes from the genomes of 1,644 species that includes representatives from most major fungal lineages. We also compiled 11 data matrices by subsampling genes or taxa from the full data matrix based on filtering criteria previously shown to improve phylogenomic inference. Analyses of these 12 data matrices using concatenation- and coalescent-based approaches yielded a robust phylogeny of the fungal kingdom, in which ∼85% of internal branches were congruent across data matrices and approaches used. We found support for several historically poorly resolved relationships as well as evidence for polytomies likely stemming from episodes of ancient diversification. By examining the relative evolutionary divergence of taxonomic groups of equivalent rank, we found that fungal taxonomy is broadly aligned with both genome sequence divergence and divergence time but also identified lineages where current taxonomic circumscription does not reflect their levels of evolutionary divergence. Our results provide a robust phylogenomic framework to explore the tempo and mode of fungal evolution and offer directions for future fungal phylogenetic and taxonomic studies.

In this review, we discuss the current status and future challenges for fully elucidating the fungal tree of life. In the last 15 years, advances in genomic technologies have revolutionized fungal systematics, ushering the field into the phylogenomic era. This has made the unthinkable possible, namely access to the entire genetic record of all known extant taxa. We first review the current status of the fungal tree and highlight areas where additional effort will be required. We then review the analytical challenges imposed by the volume of data and discuss methods to recover the most accurate species tree given the sea of gene trees. Highly resolved and deeply sampled trees are being leveraged in novel ways to study fungal radiations, species delimitation, and metabolic evolution. Finally, we discuss the critical issue of incorporating the unnamed and uncultured dark matter taxa that represent the vast majority of fungal diversity.

Phylogenomic studies based on genome-scale amounts of data have greatly improved understanding of the tree of life. Despite their diversity, ecological significance, and biomedical and industrial importance, large-scale phylogenomic studies of Fungi are lacking. Furthermore, several evolutionary relationships among major fungal lineages remain controversial, especially those at the base of the fungal phylogeny. To begin filling these gaps and assess progress toward a genome-scale phylogeny of the entire fungal kingdom, we compiled a phylogenomic data matrix of 290 genes from the genomes of 1,644 fungal species that includes representatives from most major fungal lineages; we also compiled 11 additional data matrices by subsampling genes or taxa based on filtering criteria previously shown to improve phylogenomic inference. Analyses of these 12 data matrices using concatenation- and coalescent-based approaches yielded a robust phylogeny of the kingdom in which ∼85% of internal branches were congruent across data matrices and approaches used. We found support for several relationships that have been historically contentious (e.g., for the placement of Wallemiomycotina (Basidiomycota), as sister to Agaricomycotina), as well as evidence for polytomies likely stemming from episodes of ancient diversification (e.g., at the base of Basidiomycota). By examining the relative evolutionary divergence of taxonomic groups of equivalent rank, we found that fungal taxonomy is broadly aligned with genome sequence divergence, but also identified lineages, such as the subphylum Saccharomycotina, where current taxonomic circumscription does not fully account for their high levels of evolutionary divergence. Our results provide a robust phylogenomic framework to explore the tempo and mode of fungal evolution and directions for future fungal phylogenetic and taxonomic studies.

Previous genome-scale phylogenetic analyses of Fungi have under sampled taxa from Zoopagales; this order contains many predacious or parasitic genera, and most have never been grown in pure culture. We sequenced the genomes of 4 zoopagalean taxa that are predators of amoebae, nematodes, or rotifers and the genome of one taxon that is a parasite of amoebae using single cell sequencing methods with whole genome amplification. Each genome was a metagenome, which was assembled and binned using multiple techniques to identify the target genomes. We inferred phylogenies with both super matrix and coalescent approaches using 192 conserved proteins mined from the target genomes and performed ancestral state reconstructions to determine the ancestral trophic lifestyle of the clade. Our results indicate that Zoopagales is monophyletic. Ancestral state reconstructions provide moderate support for mycoparasitism being the ancestral state of the clade.

The generation of thousands of fungal genomes is leading to a better understanding of genes and genomic organization within the kingdom. However, the epigenome, which includes DNA and chromatin modifications, remains poorly investigated in fungi. Large comparative studies in animals and plants have deepened our understanding of epigenomic variation, particularly of the modified base 5-methylcytosine (5mC), but taxonomic sampling of disparate groups is needed to develop unifying explanations for 5mC variation. Here, we utilize the largest phylogenetic resolution of 5mC methyltransferases (5mC MTases) and genome evolution to better understand levels and patterns of 5mC across fungi. We show that extant 5mC MTase genotypes are descendent from ancestral maintenance and de novo genotypes, whereas the 5mC MTases DIM-2 and RID are more recently derived, and that 5mC levels are correlated with 5mC MTase genotype and transposon content. Our survey also revealed that fungi lack canonical gene-body methylation, which distinguishes fungal epigenomes from certain insect and plant species. However, some fungal species possess independently derived clusters of contiguous 5mC encompassing many genes. In some cases, DNA repair pathways and the N6-methyladenine DNA modification negatively coevolved with 5mC pathways, which additionally contributed to interspecific epigenomic variation across fungi.

Interactions between fungi and plants, including parasitism, mutualism, and saprotrophy, have been invoked as key to their respective macroevolutionary success. Here we evaluate the origins of plant-fungal symbioses and saprotrophy using a time-calibrated phylogenetic framework that reveals linked and drastic shifts in diversification rates of each kingdom. Fungal colonization of land was associated with at least two origins of terrestrial green algae and preceded embryophytes (as evidenced by losses of fungal flagellum, ca. 720 Ma), likely facilitating terrestriality through endomycorrhizal and possibly endophytic symbioses. The largest radiation of fungi (Leotiomyceta), the origin of arbuscular mycorrhizae, and the diversification of extant embryophytes occurred ca. 480 Ma. This was followed by the origin of extant lichens. Saprotrophic mushrooms diversified in the Late Paleozoic as forests of seed plants started to dominate the landscape. The subsequent diversification and explosive radiation of Agaricomycetes, and eventually of ectomycorrhizal mushrooms, were associated with the evolution of Pinaceae in the Mesozoic, and establishment of angiosperm-dominated biomes in the Cretaceous.

Phylogenetic taxon definitions (PTDs) are explicit, phylogeny-based statements that specify clades. PTDs are central to the system of rank-free classification that is governed by the PhyloCode, but they can also be used to clarify the meanings of ranked names. We present PTDs for four major groups: Fungi, Dikarya, Ascomycota, and Basidiomycota.

Environmental DNA surveys reveal that most fungal diversity represents uncultured species. We sequenced the genomes of eight uncultured species across the fungal tree of life using a new single-cell genomics pipeline. We show that, despite a large variation in genome and gene space recovery from each single amplified genome (SAG), ≥90% can be recovered by combining multiple SAGs. SAGs provide robust placement for early-diverging lineages and infer a diploid ancestor of fungi. Early-diverging fungi share metabolic deficiencies and show unique gene expansions correlated with parasitism and unculturability. Single-cell genomics holds great promise in exploring fungal diversity, life cycles and metabolic potential.

N6-methyldeoxyadenine (6mA) is a noncanonical DNA base modification present at low levels in plant and animal genomes, but its prevalence and association with genome function in other eukaryotic lineages remains poorly understood. Here we report that abundant 6mA is associated with transcriptionally active genes in early-diverging fungal lineages. Using single-molecule long-read sequencing of 16 diverse fungal genomes, we observed that up to 2.8% of all adenines were methylated in early-diverging fungi, far exceeding levels observed in other eukaryotes and more derived fungi. 6mA occurred symmetrically at ApT dinucleotides and was concentrated in dense methylated adenine clusters surrounding the transcriptional start sites of expressed genes; its distribution was inversely correlated with that of 5-methylcytosine. Our results show a striking contrast in the genomic distributions of 6mA and 5-methylcytosine and reinforce a distinct role for 6mA as a gene-expression-associated epigenomic mark in eukaryotes.

Zygomycete fungi were classified as a single phylum, Zygomycota, based on sexual reproduction by zygospores, frequent asexual reproduction by sporangia, absence of multicellular sporocarps, and production of coenocytic hyphae, all with some exceptions. Molecular phylogenies based on one or a few genes did not support the monophyly of the phylum, however, and the phylum was subsequently abandoned. Here we present phylogenetic analyses of a genome-scale data set for 46 taxa, including 25 zygomycetes and 192 proteins, and we demonstrate that zygomycetes comprise two major clades that form a paraphyletic grade. A formal phylogenetic classification is proposed herein and includes two phyla, six subphyla, four classes and 16 orders. On the basis of these results, the phyla Mucoromycota and Zoopagomycota are circumscribed. Zoopagomycota comprises Entomophtoromycotina, Kickxellomycotina and Zoopagomycotina; it constitutes the earliest diverging lineage of zygomycetes and contains species that are primarily parasites and pathogens of small animals (e.g. amoeba, insects, etc.) and other fungi, i.e. mycoparasites. Mucoromycota comprises Glomeromycotina, Mortierellomycotina, and Mucoromycotina and is sister to Dikarya. It is the more derived clade of zygomycetes and mainly consists of mycorrhizal fungi, root endophytes, and decomposers of plant material. Evolution of trophic modes, morphology, and analysis of genome-scale data are discussed.

Fungi grow within their food, externally digesting it and absorbing nutrients across a semirigid chitinous cell wall. Members of the new phylum Cryptomycota were proposed to represent intermediate fungal forms, lacking a chitinous cell wall during feeding and known almost exclusively from ubiquitous environmental ribosomal RNA sequences that cluster at the base of the fungal tree [1, 2]. Here, we sequence the first Cryptomycotan genome (the water mold endoparasite Rozella allomycis) and unite the Cryptomycota with another group of endoparasites, the microsporidia, based on phylogenomics and shared genomic traits. We propose that Cryptomycota and microsporidia share a common endoparasitic ancestor, with the clade unified by a chitinous cell wall used to develop turgor pressure in the infection process [3, 4]. Shared genomic elements include a nucleotide transporter that is used by microsporidia for stealing energy in the form of ATP from their hosts [5]. Rozella harbors a mitochondrion that contains a very rapidly evolving genome and lacks complex I of the respiratory chain. These degenerate features are offset by the presence of nuclear genes for alternative respiratory pathways. The Rozella proteome has not undergone major contraction like microsporidia; instead, several classes have undergone expansion, such as host-effector, signal-transduction, and folding proteins.

Understanding the evolutionary history of microbial pathogens is critical for mitigating the impacts of emerging infectious diseases on economically and ecologically important host species. We used a genome resequencing approach to resolve the evolutionary history of an important microbial pathogen, the chytrid Batrachochytrium dendrobatidis (Bd), which has been implicated in amphibian declines worldwide. We sequenced the genomes of 29 isolates of Bd from around the world, with an emphasis on North, Central, and South America because of the devastating effect that Bd has had on amphibian populations in the New World. We found a substantial amount of evolutionary complexity in Bd with deep phylogenetic diversity that predates observed global amphibian declines. By investigating the entire genome, we found that even the most recently evolved Bd clade (termed the global panzootic lineage) contained more genetic variation than previously reported. We also found dramatic differences among isolates and among genomic regions in chromosomal copy number and patterns of heterozygosity, suggesting complex and heterogeneous genome dynamics. Finally, we report evidence for selection acting on the Bd genome, supporting the hypothesis that protease genes are important in evolutionary transitions in this group. Bd is considered an emerging pathogen because of its recent effects on amphibians, but our data indicate that it has a complex evolutionary history that predates recent disease outbreaks. Therefore, it is important to consider the contemporary effects of Bd in a broader evolutionary context and identify specific mechanisms that may have led to shifts in virulence in this system.

Parasitism and saprotrophic wood decay are two fungal strategies fundamental for succession and nutrient cycling in forest ecosystems. An opportunity to assess the trade-off between these strategies is provided by the forest pathogen and wood decayer Heterobasidion annosum sensu lato.We report the annotated genome sequence and transcript profiling, as well as the quantitative trait loci mapping, of one member of the species complex: H. irregulare. Quantitative trait loci critical for pathogenicity, and rich in transposable elements, orphan and secreted genes, were identified.A wide range of cellulose-degrading enzymes are expressed during wood decay. By contrast, pathogenic interaction between H. irregulare and pine engages fewer carbohydrate-active enzymes, but involves an increase in pectinolytic enzymes, transcription modules for oxidative stress and secondary metabolite production.Our results show a trade-off in terms of constrained carbohydrate decomposition and membrane transport capacity during interaction with living hosts. Our findings establish that saprotrophic wood decay and necrotrophic parasitism involve two distinct, yet overlapping, processes.