Search

Reptiles and amphibians (herptiles) are some of the most endangered and threatened species on the planet and numerous conservation strategies are being implemented with the goal of ensuring species recovery. Little is known, however, about the gut microbiome of wild herptiles and how it relates to the health of these populations. Here, we report results from the gut microbiome characterization of both a broad survey of herptiles, and the correlation between the fungus Basidiobolus, and the bacterial community supported by a deeper, more intensive sampling of Plethodon glutinosus, known as slimy salamanders. We demonstrate that bacterial communities sampled from frogs, lizards, and salamanders are structured by the host taxonomy and that Basidiobolus is a common and natural component of these wild gut microbiomes. Intensive sampling of multiple hosts across the ecoregions of Tennessee revealed that geography and host:geography interactions are strong predictors of distinct Basidiobolus operational taxonomic units present within a given host. Co-occurrence analyses of Basidiobolus and bacterial community diversity support a correlation and interaction between Basidiobolus and bacteria, suggesting that Basidiobolus may play a role in structuring the bacterial community. We further the hypothesis that this interaction is advanced by unique specialized metabolism originating from horizontal gene transfer from bacteria to Basidiobolus and demonstrate that Basidiobolus is capable of producing a diversity of specialized metabolites including small cyclic peptides.IMPORTANCEThis work significantly advances our understanding of biodiversity and microbial interactions in herptile microbiomes, the role that fungi play as a structural and functional members of herptile gut microbiomes, and the chemical functions that structure microbiome phenotypes. We also provide an important observational system of how the gut microbiome represents a unique environment that selects for novel metabolic functions through horizontal gene transfer between fungi and bacteria. Such studies are needed to better understand the complexity of gut microbiomes in nature and will inform conservation strategies for threatened species of herpetofauna.

The Lipomyces clade contains oleaginous yeast species with advantageous metabolic features for biochemical and biofuel production. Limited knowledge about the metabolic networks of the species and limited tools for genetic engineering have led to a relatively small amount of research on the microbes. Here, a genome-scale metabolic model (GSM) of Lipomyces starkeyi NRRL Y-11557 was built using orthologous protein mappings to model yeast species. Phenotypic growth assays were used to validate the GSM (66% accuracy) and indicated that NRRL Y-11557 utilized diverse carbohydrates but had more limited catabolism of organic acids. The final GSM contained 2,193 reactions, 1,909 metabolites, and 996 genes and was thus named iLst996. The model contained 96 of the annotated carbohydrate-active enzymes. iLst996 predicted a flux distribution in line with oleaginous yeast measurements and was utilized to predict theoretical lipid yields. Twenty-five other yeasts in the Lipomyces clade were then genome sequenced and annotated. Sixteen of the Lipomyces species had orthologs for more than 97% of the iLst996 genes, demonstrating the usefulness of iLst996 as a broad GSM for Lipomyces metabolism. Pathways that diverged from iLst996 mainly revolved around alternate carbon metabolism, with ortholog groups excluding NRRL Y-11557 annotated to be involved in transport, glycerolipid, and starch metabolism, among others. Overall, this study provides a useful modeling tool and data for analyzing and understanding Lipomyces species metabolism and will assist further engineering efforts in Lipomyces.

Fungi are ecologically important heterotrophs that have radiated into most niches on Earth and fulfil key ecological services. Despite intense interest in their origins, major genomic trends of their evolutionary route from a unicellular opisthokont ancestor to derived multicellular fungi remain poorly known. Here we provide a highly resolved genome-wide catalogue of gene family changes across fungal evolution inferred from the genomes of 123 fungi and relatives. We show that a dominant trend in early fungal evolution has been the gradual shedding of protist genes and the punctuated emergence of innovation by two main gene duplication events. We find that the gene content of non-Dikarya fungi resembles that of unicellular opisthokonts in many respects, owing to the conservation of protist genes in their genomes. The most rapidly duplicating gene groups included extracellular proteins and transcription factors, as well as ones linked to the coordination of nutrient uptake with growth, highlighting the transition to a sessile osmotrophic feeding strategy and subsequent lifestyle evolution as important elements of early fungal history. These results suggest that the genomes of pre-fungal ancestors evolved into the typical filamentous fungal genome by a combination of gradual gene loss, turnover and several large duplication events rather than by abrupt changes. Consequently, the taxonomically defined Fungi represents a genomically non-uniform assemblage of species.

Fungi have evolved over millions of years and their species diversity is predicted to be the second largest on the earth. Fungi have cross-kingdom interactions with many organisms that have mutually shaped their evolutionary trajectories. Zygomycete fungi hold a pivotal position in the fungal tree of life and provide important perspectives on the early evolution of fungi from aquatic to terrestrial environments. Phylogenomic analyses have found that zygomycete fungi diversified into two separate clades, the Mucoromycota which are frequently associated with plants and Zoopagomycota that are commonly animal-associated fungi. Genetic elements that contributed to the fitness and divergence of these lineages may have been shaped by the varied interactions these fungi have had with plants, animals, bacteria, and other microbes. To investigate this, we performed comparative genomic analyses of the two clades of zygomycetes in the context of Kingdom Fungi, benefiting from our generation of a new collection of zygomycete genomes, including nine produced for this study. We identified lineage-specific genomic content that may contribute to the disparate biology observed in these zygomycetes. Our findings include the discovery of undescribed diversity in CotH, a Mucormycosis pathogenicity factor, which was found in a broad set of zygomycetes. Reconciliation analysis identified multiple duplication events and an expansion of CotH copies throughout the Mucoromycotina, Mortierellomycotina, Neocallimastigomycota, and Basidiobolus lineages. A kingdom-level phylogenomic analysis also identified new evolutionary relationships within the subphyla of Mucoromycota and Zoopagomycota, including supporting the sister-clade relationship between Glomeromycotina and Mortierellomycotina and the placement of Basidiobolus as sister to other Zoopagomycota lineages.

Improved sequencing technologies have profoundly altered global views of fungal diversity and evolution. High-throughput sequencing methods are critical for studying fungi due to the cryptic, symbiotic nature of many species, particularly those that are difficult to culture. However, the low coverage genome sequencing (LCGS) approach to phylogenomic inference has not been widely applied to fungi. Here we analyzed 171 Kickxellomycotina fungi using LCGS methods to obtain hundreds of marker genes for robust phylogenomic reconstruction. Additionally, we mined our LCGS data for a set of nine rDNA and protein coding genes to enable analyses across species for which no LCGS data were obtained. The main goals of this study were to: 1) evaluate the quality and utility of LCGS data for both phylogenetic reconstruction and functional annotation, 2) test relationships among clades of Kickxellomycotina, and 3) perform comparative functional analyses between clades to gain insight into putative trophic modes. In opposition to previous studies, our nine-gene analyses support two clades of arthropod gut dwelling species and suggest a possible single evolutionary event leading to this symbiotic lifestyle. Furthermore, we resolve the mycoparasitic Dimargaritales as the earliest diverging clade in the subphylum and find four major clades of Coemansia species. Finally, functional analyses illustrate clear variation in predicted carbohydrate active enzymes and secondary metabolites (SM) based on ecology, that is biotroph versus saprotroph. Saprotrophic Kickxellales broadly lack many known pectinase families compared with saprotrophic Mucoromycota and are depauperate for SM but have similar numbers of predicted chitinases as mycoparasitic.

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.

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.

Fungi survive in diverse ecological niches by secreting proteins and other molecules into the environment to acquire food and interact with various biotic and abiotic stressors. Fungal secretome content is, therefore, believed to be tightly linked to fungal ecologies. We sampled 132 genomes from the early-diverging terrestrial fungal lineage zygomycetes (Mucoromycota and Zoopagomycota) and characterized their secretome composition. Our analyses revealed that phylogeny played an important role in shaping the secretome composition of zygomycete fungi with trophic mode contributing a smaller amount. Reconstruction of the evolution of secreted digestive enzymes revealed lineage-specific expansions, indicating that Mucoromycota and Zoopagomycota followed different trajectories early in their evolutionary history. We identified the presence of multiple pathogenicity-related proteins in the lineages known as saprotrophs, suggesting that either the ecologies of these fungi are incompletely known, and/or that these pathogenicity-related proteins have important functions associated with saprotrophic ecologies, both of which invite further investigation.

Marine fungi remain poorly covered in global genome sequencing campaigns; the 1000 fungal genomes (1KFG) project attempts to shed light on the diversity, ecology and potential industrial use of overlooked and poorly resolved fungal taxa. This study characterizes the genomes of three marine fungi: Emericellopsis sp. TS7, wood-associated Amylocarpus encephaloides and algae-associated Calycina marina. These species were genome sequenced to study their genomic features, biosynthetic potential and phylogenetic placement using multilocus data. Amylocarpus encephaloides and C. marina were placed in the Helotiaceae and Pezizellaceae (Helotiales), respectively, based on a 15-gene phylogenetic analysis. These two genomes had fewer biosynthetic gene clusters (BGCs) and carbohydrate active enzymes (CAZymes) than Emericellopsis sp. TS7 isolate. Emericellopsis sp. TS7 (Hypocreales, Ascomycota) was isolated from the sponge Stelletta normani. A six-gene phylogenetic analysis placed the isolate in the marine Emericellopsis clade and morphological examination confirmed that the isolate represents a new species, which is described here as E. atlantica. Analysis of its CAZyme repertoire and a culturing experiment on three marine and one terrestrial substrates indicated that E. atlantica is a psychrotrophic generalist fungus that is able to degrade several types of marine biomass. FungiSMASH analysis revealed the presence of 35 BGCs including, eight non-ribosomal peptide synthases (NRPSs), six NRPS-like, six polyketide synthases, nine terpenes and six hybrid, mixed or other clusters. Of these BGCs, only five were homologous with characterized BGCs. The presence of unknown BGCs sets and large CAZyme repertoire set stage for further investigations of E. atlantica. The Pezizellaceae genome and the genome of the monotypic Amylocarpus genus represent the first published genomes of filamentous fungi that are restricted in their occurrence to the marine habitat and form thus a valuable resource for the community that can be used in studying ecological adaptions of fungi using comparative genomics.

Because they comprise some of the most efficient wood-decayers, Polyporales fungi impact carbon cycling in forest environment. Despite continuous discoveries on the enzymatic machinery involved in wood decomposition, the vision on their evolutionary adaptation to wood decay and genome diversity remains incomplete. We combined the genome sequence information from 50 Polyporales species, including 26 newly sequenced genomes and sought for genomic and functional adaptations to wood decay through the analysis of genome composition and transcriptome responses to different carbon sources. The genomes of Polyporales from different phylogenetic clades showed poor conservation in macrosynteny, indicative of genome rearrangements. We observed different gene family expansion/contraction histories for plant cell wall degrading enzymes in core polyporoids and phlebioids and captured expansions for genes involved in signalling and regulation in the lineages of white rotters. Furthermore, we identified conserved cupredoxins, thaumatin-like proteins and lytic polysaccharide monooxygenases with a yet uncharacterized appended module as new candidate players in wood decomposition. Given the current need for enzymatic toolkits dedicated to the transformation of renewable carbon sources, the observed genomic diversity among Polyporales strengthens the relevance of mining Polyporales biodiversity to understand the molecular mechanisms of wood decay.

Because they comprise some of the most efficient wood-decayers, Polyporales fungi impact carbon cycling in forest environment. Despite continuous discoveries on the enzymatic machinery involved in wood decomposition, the vision on their evolutionary adaptation to wood decay and genome diversity remains incomplete. We combined the genome sequence information from 50 Polyporales species, including 26 newly sequenced genomes and sought for genomic and functional adaptations to wood decay through the analysis of genome composition and transcriptome responses to different carbon sources. The genomes of Polyporales from different phylogenetic clades showed poor conservation in macrosynteny, indicative of genome rearrangements. We observed different gene family expansion/contraction histories for plant cell wall degrading enzymes in core polyporoids and phlebioids and captured expansions for genes involved in signalling and regulation in the lineages of white rotters. Furthermore, we identified conserved cupredoxins, thaumatin-like proteins and lytic polysaccharide monooxygenases with a yet uncharacterized appended module as new candidate players in wood decomposition. Given the current need for enzymatic toolkits dedicated to the transformation of renewable carbon sources, the observed genomic diversity among Polyporales strengthens the relevance of mining Polyporales biodiversity to understand the molecular mechanisms of wood decay.

Marine fungi remain poorly covered in global genome sequencing campaigns; the 1000 fungal genomes (1KFG) project attempts to shed light on the diversity, ecology and potential industrial use of overlooked and poorly resolved fungal taxa. This study characterizes the genomes of three marine fungi: Emericellopsis sp. TS7, wood-associated Amylocarpus encephaloides and algae-associated Calycina marina. These species were genome sequenced to study their genomic features, biosynthetic potential and phylogenetic placement using multilocus data. Amylocarpus encephaloides and C. marina were placed in the Helotiaceae and Pezizellaceae (Helotiales), respectively, based on a 15-gene phylogenetic analysis. These two genomes had fewer biosynthetic gene clusters (BGCs) and carbohydrate active enzymes (CAZymes) than Emericellopsis sp. TS7 isolate. Emericellopsis sp. TS7 (Hypocreales, Ascomycota) was isolated from the sponge Stelletta normani. A six-gene phylogenetic analysis placed the isolate in the marine Emericellopsis clade and morphological examination confirmed that the isolate represents a new species, which is described here as E. atlantica. Analysis of its CAZyme repertoire and a culturing experiment on three marine and one terrestrial substrates indicated that E. atlantica is a psychrotrophic generalist fungus that is able to degrade several types of marine biomass. FungiSMASH analysis revealed the presence of 35 BGCs including, eight non-ribosomal peptide synthases (NRPSs), six NRPS-like, six polyketide synthases, nine terpenes and six hybrid, mixed or other clusters. Of these BGCs, only five were homologous with characterized BGCs. The presence of unknown BGCs sets and large CAZyme repertoire set stage for further investigations of E. atlantica. The Pezizellaceae genome and the genome of the monotypic Amylocarpus genus represent the first published genomes of filamentous fungi that are restricted in their occurrence to the marine habitat and form thus a valuable resource for the community that can be used in studying ecological adaptions of fungi using comparative genomics.

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.

The zoosporic obligate endoparasites, Olpidium, hold a pivotal position to the reconstruction of the flagellum loss in fungi, one of the key morphological transitions associated with the colonization of land by the early fungi. We generated genome and transcriptome data from non-axenic zoospores of Olpidium bornovanus and used a metagenome approach to extract phylogenetically informative fungal markers. Our phylogenetic reconstruction strongly supported Olpidium as the closest zoosporic relative of the non-flagellated terrestrial fungi. Super-alignment analyses resolved Olpidium as sister to the non-flagellated terrestrial fungi, whereas a super-tree approach recovered different placements of Olpidium, but without strong support. Further investigations detected little conflicting signal among the sampled markers but revealed a potential polytomy in early fungal evolution associated with the branching order among Olpidium, Zoopagomycota and Mucoromycota. The branches defining the evolutionary relationships of these lineages were characterized by short branch lengths and low phylogenetic content and received equivocal support for alternative phylogenetic hypotheses from individual markers. These nodes were marked by important morphological innovations, including the transition to hyphal growth and the loss of flagellum, which enabled early fungi to explore new niches and resulted in rapid and temporally concurrent Precambrian diversifications of the ancestors of several phyla of fungi.

The zoosporic obligate endoparasites, Olpidium, hold a pivotal position to the reconstruction of the flagellum loss in fungi, one of the key morphological transitions associated with the colonization of land by the early fungi. We generated genome and transcriptome data from non-axenic zoospores of Olpidium bornovanus and used a metagenome approach to extract phylogenetically informative fungal markers. Our phylogenetic reconstruction strongly supported Olpidium as the closest zoosporic relative of the non-flagellated terrestrial fungi. Super-alignment analyses resolved Olpidium as sister to the non-flagellated terrestrial fungi, whereas a super-tree approach recovered different placements of Olpidium, but without strong support. Further investigations detected little conflicting signal among the sampled markers but revealed a potential polytomy in early fungal evolution associated with the branching order among Olpidium, Zoopagomycota and Mucoromycota. The branches defining the evolutionary relationships of these lineages were characterized by short branch lengths and low phylogenetic content and received equivocal support for alternative phylogenetic hypotheses from individual markers. These nodes were marked by important morphological innovations, including the transition to hyphal growth and the loss of flagellum, which enabled early fungi to explore new niches and resulted in rapid and temporally concurrent Precambrian diversifications of the ancestors of several phyla of fungi.

Mycoviruses are widespread and purportedly common throughout the fungal kingdom, although most are known from hosts in the two most recently diverged phyla, Ascomycota and Basidiomycota, together called Dikarya. To augment our knowledge of mycovirus prevalence and diversity in underexplored fungi, we conducted a large-scale survey of fungi in the earlier-diverging lineages, using both culture-based and transcriptome-mining approaches to search for RNA viruses. In total, 21.6% of 333 isolates were positive for RNA mycoviruses. This is a greater proportion than expected based on previous taxonomically broad mycovirus surveys and is suggestive of a strong phylogenetic component to mycoviral infection. Our newly found viral sequences are diverse, composed of double-stranded RNA, positive-sense single-stranded RNA (ssRNA), and negative-sense ssRNA genomes and include novel lineages lacking representation in the public databases. These identified viruses could be classified into 2 orders, 5 families, and 5 genera; however, half of the viruses remain taxonomically unassigned. Further, we identified a lineage of virus-like sequences in the genomes of members of Phycomycetaceae and Mortierellales that appear to be novel genes derived from integration of a viral RNA-dependent RNA polymerase gene. The two screening methods largely agreed in their detection of viruses; thus, we suggest that the culture-based assay is a cost-effective means to quickly assess whether a laboratory culture is virally infected. This study used culture collections and publicly available transcriptomes to demonstrate that mycoviruses are abundant in laboratory cultures of early-diverging fungal lineages. The function and diversity of mycoviruses found here will help guide future studies into mycovirus origins and ecological functions.IMPORTANCE Viruses are key drivers of evolution and ecosystem function and are increasingly recognized as symbionts of fungi. Fungi in early-diverging lineages are widespread, ecologically important, and comprise the majority of the phylogenetic diversity of the kingdom. Viruses infecting early-diverging lineages of fungi have been almost entirely unstudied. In this study, we screened fungi for viruses by two alternative approaches: a classic culture-based method and by transcriptome-mining. The results of our large-scale survey demonstrate that early-diverging lineages have higher infection rates than have been previously reported in other fungal taxa and that laboratory strains worldwide are host to infections, the implications of which are unknown. The function and diversity of mycoviruses found in these basal fungal lineages will help guide future studies into mycovirus origins and their evolutionary ramifications and ecological impacts.

Mycorrhizal fungi are mutualists that play crucial roles in nutrient acquisition in terrestrial ecosystems. Mycorrhizal symbioses arose repeatedly across multiple lineages of Mucoromycotina, Ascomycota, and Basidiomycota. Considerable variation exists in the capacity of mycorrhizal fungi to acquire carbon from soil organic matter. Here, we present a combined analysis of 135 fungal genomes from 73 saprotrophic, endophytic and pathogenic species, and 62 mycorrhizal species, including 29 new mycorrhizal genomes. This study samples ecologically dominant fungal guilds for which there were previously no symbiotic genomes available, including ectomycorrhizal Russulales, Thelephorales and Cantharellales. Our analyses show that transitions from saprotrophy to symbiosis involve (1) widespread losses of degrading enzymes acting on lignin and cellulose, (2) co-option of genes present in saprotrophic ancestors to fulfill new symbiotic functions, (3) diversification of novel, lineage-specific symbiosis-induced genes, (4) proliferation of transposable elements and (5) divergent genetic innovations underlying the convergent origins of the ectomycorrhizal guild.

Ecological diversity in fungi is largely defined by metabolic traits, including the ability to produce secondary or "specialized" metabolites (SMs) that mediate interactions with other organisms. Fungal SM pathways are frequently encoded in biosynthetic gene clusters (BGCs), which facilitate the identification and characterization of metabolic pathways. Variation in BGC composition reflects the diversity of their SM products. Recent studies have documented surprising diversity of BGC repertoires among isolates of the same fungal species, yet little is known about how this population-level variation is inherited across macroevolutionary timescales. Here, we applied a novel linkage-based algorithm to reveal previously unexplored dimensions of diversity in BGC composition, distribution, and repertoire across 101 species of Dothideomycetes, which are considered the most phylogenetically diverse class of fungi and known to produce many SMs. We predicted both complementary and overlapping sets of clustered genes compared with existing methods and identified novel gene pairs that associate with known secondary metabolite genes. We found that variation among sets of BGCs in individual genomes is due to nonoverlapping BGC combinations and that several BGCs have biased ecological distributions, consistent with niche-specific selection. We observed that total BGC diversity scales linearly with increasing repertoire size, suggesting that secondary metabolites have little structural redundancy in individual fungi. We project that there is substantial unsampled BGC diversity across specific families of Dothideomycetes, which will provide a roadmap for future sampling efforts. Our approach and findings lend new insight into how BGC diversity is generated and maintained across an entire fungal taxonomic class.

Ecological diversity in fungi is largely defined by metabolic traits, including the ability to produce secondary or "specialized" metabolites (SMs) that mediate interactions with other organisms. Fungal SM pathways are frequently encoded in biosynthetic gene clusters (BGCs), which facilitate the identification and characterization of metabolic pathways. Variation in BGC composition reflects the diversity of their SM products. Recent studies have documented surprising diversity of BGC repertoires among isolates of the same fungal species, yet little is known about how this population-level variation is inherited across macroevolutionary timescales. Here, we applied a novel linkage-based algorithm to reveal previously unexplored dimensions of diversity in BGC composition, distribution, and repertoire across 101 species of Dothideomycetes, which are considered the most phylogenetically diverse class of fungi and known to produce many SMs. We predicted both complementary and overlapping sets of clustered genes compared with existing methods and identified novel gene pairs that associate with known secondary metabolite genes. We found that variation among sets of BGCs in individual genomes is due to nonoverlapping BGC combinations and that several BGCs have biased ecological distributions, consistent with niche-specific selection. We observed that total BGC diversity scales linearly with increasing repertoire size, suggesting that secondary metabolites have little structural redundancy in individual fungi. We project that there is substantial unsampled BGC diversity across specific families of Dothideomycetes, which will provide a roadmap for future sampling efforts. Our approach and findings lend new insight into how BGC diversity is generated and maintained across an entire fungal taxonomic class.

Ecological diversity in fungi is largely defined by metabolic traits, including the ability to produce secondary or "specialized" metabolites (SMs) that mediate interactions with other organisms. Fungal SM pathways are frequently encoded in biosynthetic gene clusters (BGCs), which facilitate the identification and characterization of metabolic pathways. Variation in BGC composition reflects the diversity of their SM products. Recent studies have documented surprising diversity of BGC repertoires among isolates of the same fungal species, yet little is known about how this population-level variation is inherited across macroevolutionary timescales. Here, we applied a novel linkage-based algorithm to reveal previously unexplored dimensions of diversity in BGC composition, distribution, and repertoire across 101 species of Dothideomycetes, which are considered the most phylogenetically diverse class of fungi and known to produce many SMs. We predicted both complementary and overlapping sets of clustered genes compared with existing methods and identified novel gene pairs that associate with known secondary metabolite genes. We found that variation among sets of BGCs in individual genomes is due to nonoverlapping BGC combinations and that several BGCs have biased ecological distributions, consistent with niche-specific selection. We observed that total BGC diversity scales linearly with increasing repertoire size, suggesting that secondary metabolites have little structural redundancy in individual fungi. We project that there is substantial unsampled BGC diversity across specific families of Dothideomycetes, which will provide a roadmap for future sampling efforts. Our approach and findings lend new insight into how BGC diversity is generated and maintained across an entire fungal taxonomic class.

Mycoviruses are widespread and purportedly common throughout the fungal kingdom, although most are known from hosts in the two most recently diverged phyla, Ascomycota and Basidiomycota, together called Dikarya. To augment our knowledge of mycovirus prevalence and diversity in underexplored fungi, we conducted a large-scale survey of fungi in the earlier-diverging lineages, using both culture-based and transcriptome-mining approaches to search for RNA viruses. In total, 21.6% of 333 isolates were positive for RNA mycoviruses. This is a greater proportion than expected based on previous taxonomically broad mycovirus surveys and is suggestive of a strong phylogenetic component to mycoviral infection. Our newly found viral sequences are diverse, composed of double-stranded RNA, positive-sense single-stranded RNA (ssRNA), and negative-sense ssRNA genomes and include novel lineages lacking representation in the public databases. These identified viruses could be classified into 2 orders, 5 families, and 5 genera; however, half of the viruses remain taxonomically unassigned. Further, we identified a lineage of virus-like sequences in the genomes of members of Phycomycetaceae and Mortierellales that appear to be novel genes derived from integration of a viral RNA-dependent RNA polymerase gene. The two screening methods largely agreed in their detection of viruses; thus, we suggest that the culture-based assay is a cost-effective means to quickly assess whether a laboratory culture is virally infected. This study used culture collections and publicly available transcriptomes to demonstrate that mycoviruses are abundant in laboratory cultures of early-diverging fungal lineages. The function and diversity of mycoviruses found here will help guide future studies into mycovirus origins and ecological functions.IMPORTANCE Viruses are key drivers of evolution and ecosystem function and are increasingly recognized as symbionts of fungi. Fungi in early-diverging lineages are widespread, ecologically important, and comprise the majority of the phylogenetic diversity of the kingdom. Viruses infecting early-diverging lineages of fungi have been almost entirely unstudied. In this study, we screened fungi for viruses by two alternative approaches: a classic culture-based method and by transcriptome-mining. The results of our large-scale survey demonstrate that early-diverging lineages have higher infection rates than have been previously reported in other fungal taxa and that laboratory strains worldwide are host to infections, the implications of which are unknown. The function and diversity of mycoviruses found in these basal fungal lineages will help guide future studies into mycovirus origins and their evolutionary ramifications and ecological impacts.

Research into secondary metabolism (SM) production by fungi has resulted in the discovery of diverse, biologically active compounds with significant medicinal applications. The fungi rich in SM production are taxonomically concentrated in the subkingdom Dikarya, which comprises the phyla Ascomycota and Basidiomycota. Here, we explore the potential for SM production in Mucoromycota and Zoopagomycota, two phyla of nonflagellated fungi that are not members of Dikarya, by predicting and identifying core genes and gene clusters involved in SM. The majority of non-Dikarya have few genes and gene clusters involved in SM production except for the amphibian gut symbionts in the genus Basidiobolus Basidiobolus genomes exhibit an enrichment of SM genes involved in siderophore, surfactin-like, and terpene cyclase production, all these with evidence of constitutive gene expression. Gene expression and chemical assays also confirm that Basidiobolus has significant siderophore activity. The expansion of SMs in Basidiobolus are partially due to horizontal gene transfer from bacteria, likely as a consequence of its ecology as an amphibian gut endosymbiont.

Research into secondary metabolism (SM) production by fungi has resulted in the discovery of diverse, biologically active compounds with significant medicinal applications. The fungi rich in SM production are taxonomically concentrated in the subkingdom Dikarya, which comprises the phyla Ascomycota and Basidiomycota. Here, we explore the potential for SM production in Mucoromycota and Zoopagomycota, two phyla of nonflagellated fungi that are not members of Dikarya, by predicting and identifying core genes and gene clusters involved in SM. The majority of non-Dikarya have few genes and gene clusters involved in SM production except for the amphibian gut symbionts in the genus BasidiobolusBasidiobolus genomes exhibit an enrichment of SM genes involved in siderophore, surfactin-like, and terpene cyclase production, all these with evidence of constitutive gene expression. Gene expression and chemical assays also confirm that Basidiobolus has significant siderophore activity. The expansion of SMs in Basidiobolus are partially due to horizontal gene transfer from bacteria, likely as a consequence of its ecology as an amphibian gut endosymbiont.

Dothideomycetes is the largest class of kingdom Fungi and comprises an incredible diversity of lifestyles, many of which have evolved multiple times. Plant pathogens represent a major ecological niche of the class Dothideomycetes and they are known to infect most major food crops and feedstocks for biomass and biofuel production. Studying the ecology and evolution of Dothideomycetes has significant implications for our fundamental understanding of fungal evolution, their adaptation to stress and host specificity, and practical implications with regard to the effects of climate change and on the food, feed, and livestock elements of the agro-economy. In this study, we present the first large-scale, whole-genome comparison of 101 Dothideomycetes introducing 55 newly sequenced species. The availability of whole-genome data produced a high-confidence phylogeny leading to reclassification of 25 organisms, provided a clearer picture of the relationships among the various families, and indicated that pathogenicity evolved multiple times within this class. We also identified gene family expansions and contractions across the Dothideomycetes phylogeny linked to ecological niches providing insights into genome evolution and adaptation across this group. Using machine-learning methods we classified fungi into lifestyle classes with >95 % accuracy and identified a small number of gene families that positively correlated with these distinctions. This can become a valuable tool for genome-based prediction of species lifestyle, especially for rarely seen and poorly studied species.

Dothideomycetes is the largest class of kingdom Fungi and comprises an incredible diversity of lifestyles, many of which have evolved multiple times. Plant pathogens represent a major ecological niche of the class Dothideomycetes and they are known to infect most major food crops and feedstocks for biomass and biofuel production. Studying the ecology and evolution of Dothideomycetes has significant implications for our fundamental understanding of fungal evolution, their adaptation to stress and host specificity, and practical implications with regard to the effects of climate change and on the food, feed, and livestock elements of the agro-economy. In this study, we present the first large-scale, whole-genome comparison of 101 Dothideomycetes introducing 55 newly sequenced species. The availability of whole-genome data produced a high-confidence phylogeny leading to reclassification of 25 organisms, provided a clearer picture of the relationships among the various families, and indicated that pathogenicity evolved multiple times within this class. We also identified gene family expansions and contractions across the Dothideomycetes phylogeny linked to ecological niches providing insights into genome evolution and adaptation across this group. Using machine-learning methods we classified fungi into lifestyle classes with >95 % accuracy and identified a small number of gene families that positively correlated with these distinctions. This can become a valuable tool for genome-based prediction of species lifestyle, especially for rarely seen and poorly studied species.

Research into secondary metabolism (SM) production by fungi has resulted in the discovery of diverse, biologically active compounds with significant medicinal applications. However, the fungi rich in SM production are taxonomically restricted to Dikarya, two phyla of Kingdom Fungi, Ascomycota and Basidiomycota. Here, we explore the potential for SM production in Mucoromycota and Zoopagomycota, two phyla of nonflagellated fungi that are not members of Dikarya, by predicting and identifying core genes and gene clusters involved in SM. The majority of non-Dikarya have few genes and gene clusters involved in SM production except for the amphibian gut symbionts in the genus Basidiobolus . Basidiobolus genomes exhibit an enrichment of SM genes involved in siderophore, surfactin-like, and terpene cyclase production, all these with evidence of constitutive gene expression. Gene expression and chemical assays confirm that Basidiobolus has significant siderophore activity. The expansion of SMs in Basidiobolus are partially due to horizontal gene transfer from bacteria, likely as a consequence of its ecology as an amphibian gut endosymbiont.

Abstract: Background Ecological diversity in fungi is largely defined by metabolic traits, including the ability to produce secondary or “specialized” metabolites (SMs) that mediate interactions with other organisms. Fungal SM pathways are frequently encoded in biosynthetic gene clusters (BGCs), which facilitate the identification and characterization of metabolic pathways. Variation in BGC composition reflects the diversity of their SM products. Recent studies have documented surprising diversity of BGC repertoires among isolates of the same fungal species, yet little is known about how this population-level variation is inherited across macroevolutionary timescales. Results Here, we applied a novel linkage-based algorithm to reveal previously unexplored dimensions of diversity in BGC composition, distribution, and repertoire across 101 species of Dothideomycetes, which are considered to be the most phylogenetically diverse class of fungi and are known to produce many SMs. We predicted both complementary and overlapping sets of clustered genes compared with existing methods and identified novel gene pairs that associate with known secondary metabolite genes. We found that variation in BGC repertoires is due to non-overlapping BGC combinations and that several BGCs have biased ecological distributions, consistent with niche-specific selection. We observed that total BGC diversity scales linearly with increasing repertoire size, suggesting that secondary metabolites have little structural redundancy in individual fungi. Conclusions We project that there is substantial unsampled BGC diversity across specific families of Dothideomycetes, which will provide a roadmap for future sampling efforts. Our approach and findings lend new insight into how BGC diversity is generated and maintained across an entire fungal taxonomic class.

Mycorrhizal fungi are mutualists that play crucial roles in nutrient acquisition in terrestrial ecosystems. Mycorrhizal symbioses arose repeatedly across multiple lineages of Mucoromycotina, Ascomycota, and Basidiomycota. Considerable variation exists in the capacity of mycorrhizal fungi to acquire carbon from soil organic matter. Here, we present a combined analysis of 135 fungal genomes from 73 saprotrophic, endophytic and pathogenic species, and 62 mycorrhizal species, including 29 new mycorrhizal genomes. This study samples ecologically dominant fungal guilds for which there were previously no symbiotic genomes available, including ectomycorrhizal Russulales, Thelephorales and Cantharellales. Our analyses show that transitions from saprotrophy to symbiosis involve (1) widespread losses of degrading enzymes acting on lignin and cellulose, (2) co-option of genes present in saprotrophic ancestors to fulfill new symbiotic functions, (3) diversification of novel, lineage-specific symbiosis-induced genes, (4) proliferation of transposable elements and (5) divergent genetic innovations underlying the convergent origins of the ectomycorrhizal guild.

Abstract: Background Ecological diversity in fungi is largely defined by metabolic traits, including the ability to produce secondary or “specialized” metabolites (SMs) that mediate interactions with other organisms. Fungal SM pathways are frequently encoded in biosynthetic gene clusters (BGCs), which facilitate the identification and characterization of metabolic pathways. Variation in BGC composition reflects the diversity of their SM products. Recent studies have documented surprising diversity of BGC repertoires among isolates of the same fungal species, yet little is known about how this population-level variation is inherited across macroevolutionary timescales. Results Here, we applied a novel linkage-based algorithm to reveal previously unexplored dimensions of diversity in BGC composition, distribution, and repertoire across 101 species of Dothideomycetes, which are considered to be the most phylogenetically diverse class of fungi and are known to produce many SMs. We predicted both complementary and overlapping sets of clustered genes compared with existing methods and identified novel gene pairs that associate with known secondary metabolite genes. We found that variation in BGC repertoires is due to non-overlapping BGC combinations and that several BGCs have biased ecological distributions, consistent with niche-specific selection. We observed that total BGC diversity scales linearly with increasing repertoire size, suggesting that secondary metabolites have little structural redundancy in individual fungi. Conclusions We project that there is substantial unsampled BGC diversity across specific families of Dothideomycetes, which will provide a roadmap for future sampling efforts. Our approach and findings lend new insight into how BGC diversity is generated and maintained across an entire fungal taxonomic class.

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.

The zoosporic obligate endoparasites, Olpidium , hold a pivotal position to the reconstruction of the flagellum loss in fungi, one of the key morphological transitions associated with the colonization of land by the early fungi. We generated genome and transcriptome data from non-axenic zoospores of Olpidium bornovanus and used a metagenome approach to extract phylogenetically informative fungal markers. Our phylogenetic reconstruction strongly supported Olpidium as the closest zoosporic relative of the non-flagellated terrestrial fungi. Super-alignment analyses resolved Olpidium as sister to the non-flagellated terrestrial fungi, whereas a super-tree approach recovered different placements of Olpidium , but without strong support. Further investigations detected little conflicting signal among the sampled markers but revealed a potential polytomy in early fungal evolution associated with the branching order among Olpidium , Zoopagomycota and Mucoromycota. The branches defining the evolutionary relationships of these lineages were characterized by short branch lengths and low phylogenetic content and received equivocal support for alternative phylogenetic hypotheses from individual markers. These nodes were marked by important morphological innovations, including the transition to hyphal growth and the loss of flagellum, which enabled early fungi to explore new niches and resulted in rapid and temporally concurrent Precambrian diversifications of the ancestors of several phyla of fungi.

Research into secondary metabolism (SM) production by fungi has resulted in the discovery of diverse, biologically active compounds with significant medicinal applications. However, the fungi rich in SM production are taxonomically restricted to Dikarya, two phyla of Kingdom Fungi, Ascomycota and Basidiomycota. Here, we explore the potential for SM production in Mucoromycota and Zoopagomycota, two phyla of nonflagellated fungi that are not members of Dikarya, by predicting and identifying core genes and gene clusters involved in SM. The majority of non-Dikarya have few genes and gene clusters involved in SM production except for the amphibian gut symbionts in the genus Basidiobolus . Basidiobolus genomes exhibit an enrichment of SM genes involved in siderophore, surfactin-like, and terpene cyclase production, all these with evidence of constitutive gene expression. Gene expression and chemical assays confirm that Basidiobolus has significant siderophore activity. The expansion of SMs in Basidiobolus are partially due to horizontal gene transfer from bacteria, likely as a consequence of its ecology as an amphibian gut endosymbiont.

Several Phyllosticta species are known as pathogens of Citrus spp., and are responsible for various disease symptoms including leaf and fruit spots. One of the most important species is P. citricarpa, which causes a foliar and fruit disease called citrus black spot. The Phyllosticta species occurring on citrus can most effectively be distinguished from P. citricarpa by means of multilocus DNA sequence data. Recent studies also demonstrated P. citricarpa to be heterothallic, and reported successful mating in the laboratory. Since the domestication of citrus, different clones of P. citricarpa have escaped Asia to other continents via trade routes, with obvious disease management consequences. This pathogen profile represents a comprehensive literature review of this pathogen and allied taxa associated with citrus, focusing on identification, distribution, genomics, epidemiology and disease management. This review also considers the knowledge emerging from seven genomes of Phyllosticta spp., demonstrating unknown aspects of these species, including their mating behaviour.TaxonomyPhyllosticta citricarpa (McAlpine) Aa, 1973. Kingdom Fungi, Phylum Ascomycota, Class Dothideomycetes, Order Botryosphaeriales, Family Phyllostictaceae, Genus Phyllosticta, Species citricarpa.Host rangeConfirmed on more than 12 Citrus species, Phyllosticta citricarpa has only been found on plant species in the Rutaceae.Disease symptomsP. citricarpa causes diverse symptoms such as hard spot, virulent spot, false melanose and freckle spot on fruit, and necrotic lesions on leaves and twigs.Useful websitesDOE Joint Genome Institute MycoCosm portals for the Phyllosticta capitalensis (https://genome.jgi.doe.gov/Phycap1), P. citriasiana (https://genome.jgi.doe.gov/Phycit1), P. citribraziliensis (https://genome.jgi.doe.gov/Phcit1), P. citrichinaensis (https://genome.jgi.doe.gov/Phcitr1), P. citricarpa (https://genome.jgi.doe.gov/Phycitr1, https://genome.jgi.doe.gov/Phycpc1), P. paracitricarpa (https://genome.jgi.doe.gov/Phy27169) genomes. All available Phyllosticta genomes on MycoCosm can be viewed at https://genome.jgi.doe.gov/Phyllosticta.

Several Phyllosticta species are known as pathogens of Citrus spp., and are responsible for various disease symptoms including leaf and fruit spots. One of the most important species is P. citricarpa, which causes a foliar and fruit disease called citrus black spot. The Phyllosticta species occurring on citrus can most effectively be distinguished from P. citricarpa by means of multilocus DNA sequence data. Recent studies also demonstrated P. citricarpa to be heterothallic, and reported successful mating in the laboratory. Since the domestication of citrus, different clones of P. citricarpa have escaped Asia to other continents via trade routes, with obvious disease management consequences. This pathogen profile represents a comprehensive literature review of this pathogen and allied taxa associated with citrus, focusing on identification, distribution, genomics, epidemiology and disease management. This review also considers the knowledge emerging from seven genomes of Phyllosticta spp., demonstrating unknown aspects of these species, including their mating behaviour.TaxonomyPhyllosticta citricarpa (McAlpine) Aa, 1973. Kingdom Fungi, Phylum Ascomycota, Class Dothideomycetes, Order Botryosphaeriales, Family Phyllostictaceae, Genus Phyllosticta, Species citricarpa.Host rangeConfirmed on more than 12 Citrus species, Phyllosticta citricarpa has only been found on plant species in the Rutaceae.Disease symptomsP. citricarpa causes diverse symptoms such as hard spot, virulent spot, false melanose and freckle spot on fruit, and necrotic lesions on leaves and twigs.Useful websitesDOE Joint Genome Institute MycoCosm portals for the Phyllosticta capitalensis (https://genome.jgi.doe.gov/Phycap1), P. citriasiana (https://genome.jgi.doe.gov/Phycit1), P. citribraziliensis (https://genome.jgi.doe.gov/Phcit1), P. citrichinaensis (https://genome.jgi.doe.gov/Phcitr1), P. citricarpa (https://genome.jgi.doe.gov/Phycitr1, https://genome.jgi.doe.gov/Phycpc1), P. paracitricarpa (https://genome.jgi.doe.gov/Phy27169) genomes. All available Phyllosticta genomes on MycoCosm can be viewed at https://genome.jgi.doe.gov/Phyllosticta.

The environmental accumulation of polycyclic aromatic hydrocarbons (PAHs) is of great concern due to potential carcinogenic and mutagenic risks, as well as their resistance to remediation. While many fungi have been reported to break down PAHs in environments, the details of gene-based metabolic pathways are not yet comprehensively understood. Specifically, the genome-scale transcriptional responses of fungal PAH degradation have rarely been reported. In this study, we report the genomic and transcriptomic basis of PAH bioremediation by a potent fungal degrader, Dentipellis sp. KUC8613. The genome size of this fungus was 36.71 Mbp long encoding 14,320 putative protein-coding genes. The strain efficiently removed more than 90% of 100 mg/l concentration of PAHs within 10 days. The genomic and transcriptomic analysis of this white rot fungus highlights that the strain primarily utilized non-ligninolytic enzymes to remove various PAHs, rather than typical ligninolytic enzymes known for playing important roles in PAH degradation. PAH removal by non-ligninolytic enzymes was initiated by both different PAH-specific and common upregulation of P450s, followed by downstream PAH-transforming enzymes such as epoxide hydrolases, dehydrogenases, FAD-dependent monooxygenases, dioxygenases, and glycosyl- or glutathione transferases. Among the various PAHs, phenanthrene induced a more dynamic transcriptomic response possibly due to its greater cytotoxicity, leading to highly upregulated genes involved in the translocation of PAHs, a defense system against reactive oxygen species, and ATP synthesis. Our genomic and transcriptomic data provide a foundation of understanding regarding the mycoremediation of PAHs and the application of this strain for polluted environments.

The environmental accumulation of polycyclic aromatic hydrocarbons (PAHs) is of great concern due to potential carcinogenic and mutagenic risks, as well as their resistance to remediation. While many fungi have been reported to break down PAHs in environments, the details of gene-based metabolic pathways are not yet comprehensively understood. Specifically, the genome-scale transcriptional responses of fungal PAH degradation have rarely been reported. In this study, we report the genomic and transcriptomic basis of PAH bioremediation by a potent fungal degrader, Dentipellis sp. KUC8613. The genome size of this fungus was 36.71 Mbp long encoding 14,320 putative protein-coding genes. The strain efficiently removed more than 90% of 100 mg/l concentration of PAHs within 10 days. The genomic and transcriptomic analysis of this white rot fungus highlights that the strain primarily utilized non-ligninolytic enzymes to remove various PAHs, rather than typical ligninolytic enzymes known for playing important roles in PAH degradation. PAH removal by non-ligninolytic enzymes was initiated by both different PAH-specific and common upregulation of P450s, followed by downstream PAH-transforming enzymes such as epoxide hydrolases, dehydrogenases, FAD-dependent monooxygenases, dioxygenases, and glycosyl- or glutathione transferases. Among the various PAHs, phenanthrene induced a more dynamic transcriptomic response possibly due to its greater cytotoxicity, leading to highly upregulated genes involved in the translocation of PAHs, a defense system against reactive oxygen species, and ATP synthesis. Our genomic and transcriptomic data provide a foundation of understanding regarding the mycoremediation of PAHs and the application of this strain for polluted environments.

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.

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.

Endogonales (Mucoromycotina), composed of Endogonaceae and Densosporaceae, is the only known non-Dikarya order with ectomycorrhizal members. They also form mycorrhizal-like association with some nonspermatophyte plants. It has been recently proposed that Endogonales were among the earliest mycorrhizal partners with land plants. It remains unknown whether Endogonales possess genomes with mycorrhizal-lifestyle signatures and whether Endogonales originated around the same time as land plants did. We sampled sporocarp tissue from four Endogonaceae collections and performed shotgun genome sequencing. After binning the metagenome data, we assembled and annotated the Endogonaceae genomes. We performed comparative analysis on plant-cell-wall-degrading enzymes (PCWDEs) and small secreted proteins (SSPs). We inferred phylogenetic placement of Endogonaceae and estimated the ages of Endogonaceae and Endogonales with expanded taxon sampling. Endogonaceae have large genomes with high repeat content, low diversity of PCWDEs, but without elevated SSP/secretome ratios. Dating analysis estimated that Endogonaceae originated in the Permian-Triassic boundary and Endogonales originated in the mid-late Silurian. Mycoplasma-related endobacterium sequences were identified in three Endogonaceae genomes. Endogonaceae genomes possess typical signatures of mycorrhizal lifestyle. The early origin of Endogonales suggests that the mycorrhizal association between Endogonales and plants might have played an important role during the colonization of land by plants.

Endogonales (Mucoromycotina), composed of Endogonaceae and Densosporaceae, is the only known non-Dikarya order with ectomycorrhizal members. They also form mycorrhizal-like association with some nonspermatophyte plants. It has been recently proposed that Endogonales were among the earliest mycorrhizal partners with land plants. It remains unknown whether Endogonales possess genomes with mycorrhizal-lifestyle signatures and whether Endogonales originated around the same time as land plants did. We sampled sporocarp tissue from four Endogonaceae collections and performed shotgun genome sequencing. After binning the metagenome data, we assembled and annotated the Endogonaceae genomes. We performed comparative analysis on plant-cell-wall-degrading enzymes (PCWDEs) and small secreted proteins (SSPs). We inferred phylogenetic placement of Endogonaceae and estimated the ages of Endogonaceae and Endogonales with expanded taxon sampling. Endogonaceae have large genomes with high repeat content, low diversity of PCWDEs, but without elevated SSP/secretome ratios. Dating analysis estimated that Endogonaceae originated in the Permian-Triassic boundary and Endogonales originated in the mid-late Silurian. Mycoplasma-related endobacterium sequences were identified in three Endogonaceae genomes. Endogonaceae genomes possess typical signatures of mycorrhizal lifestyle. The early origin of Endogonales suggests that the mycorrhizal association between Endogonales and plants might have played an important role during the colonization of land by plants.

We identified two poplar (Populus sp.)-associated microbes, the fungus, Mortierella elongata strain AG77, and the bacterium, Burkholderia strain BT03, that mutually promote each other's growth. Using culture assays in concert with a novel microfluidic device to generate time-lapse videos, we found growth specific media differing in pH and pre-conditioned by microbial growth led to increased fungal and bacterial growth rates. Coupling microfluidics and comparative metabolomics data results indicated that observed microbial growth stimulation involves metabolic exchange during two ordered events. The first is an emission of fungal metabolites, including organic acids used or modified by bacteria. A second signal of unknown nature is produced by bacteria which increases fungal growth rates. We find this symbiosis is initiated in part by metabolic exchange involving fungal organic acids.

We identified two poplar (Populus sp.)-associated microbes, the fungus, Mortierella elongata strain AG77, and the bacterium, Burkholderia strain BT03, that mutually promote each other's growth. Using culture assays in concert with a novel microfluidic device to generate time-lapse videos, we found growth specific media differing in pH and pre-conditioned by microbial growth led to increased fungal and bacterial growth rates. Coupling microfluidics and comparative metabolomics data results indicated that observed microbial growth stimulation involves metabolic exchange during two ordered events. The first is an emission of fungal metabolites, including organic acids used or modified by bacteria. A second signal of unknown nature is produced by bacteria which increases fungal growth rates. We find this symbiosis is initiated in part by metabolic exchange involving fungal organic acids.

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.

Tuberaceae is one of the most diverse lineages of symbiotic truffle-forming fungi. To understand the molecular underpinning of the ectomycorrhizal truffle lifestyle, we compared the genomes of Piedmont white truffle (Tuber magnatum), Périgord black truffle (Tuber melanosporum), Burgundy truffle (Tuber aestivum), pig truffle (Choiromyces venosus) and desert truffle (Terfezia boudieri) to saprotrophic Pezizomycetes. Reconstructed gene duplication/loss histories along a time-calibrated phylogeny of Ascomycetes revealed that Tuberaceae-specific traits may be related to a higher gene diversification rate. Genomic features in Tuber species appear to be very similar, with high transposon content, few genes coding lignocellulose-degrading enzymes, a substantial set of lineage-specific fruiting-body-upregulated genes and high expression of genes involved in volatile organic compound metabolism. Developmental and metabolic pathways expressed in ectomycorrhizae and fruiting bodies of T. magnatum and T. melanosporum are unexpectedly very similar, owing to the fact that they diverged ~100 Ma. Volatile organic compounds from pungent truffle odours are not the products of Tuber-specific gene innovations, but rely on the differential expression of an existing gene repertoire. These genomic resources will help to address fundamental questions in the evolution of the truffle lifestyle and the ecology of fungi that have been praised as food delicacies for centuries.

Tuberaceae is one of the most diverse lineages of symbiotic truffle-forming fungi. To understand the molecular underpinning of the ectomycorrhizal truffle lifestyle, we compared the genomes of Piedmont white truffle (Tuber magnatum), Périgord black truffle (Tuber melanosporum), Burgundy truffle (Tuber aestivum), pig truffle (Choiromyces venosus) and desert truffle (Terfezia boudieri) to saprotrophic Pezizomycetes. Reconstructed gene duplication/loss histories along a time-calibrated phylogeny of Ascomycetes revealed that Tuberaceae-specific traits may be related to a higher gene diversification rate. Genomic features in Tuber species appear to be very similar, with high transposon content, few genes coding lignocellulose-degrading enzymes, a substantial set of lineage-specific fruiting-body-upregulated genes and high expression of genes involved in volatile organic compound metabolism. Developmental and metabolic pathways expressed in ectomycorrhizae and fruiting bodies of T. magnatum and T. melanosporum are unexpectedly very similar, owing to the fact that they diverged ~100 Ma. Volatile organic compounds from pungent truffle odours are not the products of Tuber-specific gene innovations, but rely on the differential expression of an existing gene repertoire. These genomic resources will help to address fundamental questions in the evolution of the truffle lifestyle and the ecology of fungi that have been praised as food delicacies for centuries.

Ustilaginomycotina is home to a broad array of fungi including important plant pathogens collectively called smut fungi. Smuts are biotrophs that produce characteristic perennating propagules called teliospores, one of which, Ustilago maydis, is a model genetic organism. Broad exploration of smut biology has been hampered by limited phylogenetic resolution of Ustilaginiomycotina as well as an overall lack of genomic data for members of this subphylum. In this study, we sequenced eight Ustilaginomycotina genomes from previously unrepresented lineages, deciphered ordinal-level phylogenetic relationships for the subphylum, and performed comparative analyses. Unlike other Basidiomycota subphyla, all sampled Ustilaginomycotina genomes are relatively small and compact. Ancestral state reconstruction analyses indicate that teliospore formation was present at the origin of the subphylum. Divergence time estimation dates the divergence of most extant smut fungi after that of grasses (Poaceae). However, we found limited conservation of well-characterized genes related to smut pathogenesis from U. maydis, indicating dissimilar pathogenic mechanisms exist across other smut lineages. The genomes of Malasseziomycetes are highly diverged from the other sampled Ustilaginomycotina, likely due to their unique history as mammal-associated lipophilic yeasts. Despite extensive genomic data, the phylogenetic placement of this class remains ambiguous. Although the sampled Ustilaginomycotina members lack many core enzymes for plant cell wall decomposition and starch catabolism, we identified several novel carbohydrate active enzymes potentially related to pectin breakdown. Finally, ∼50% of Ustilaginomycotina species-specific genes are present in previously undersampled and rare lineages, highlighting the importance of exploring fungal diversity as a resource for novel gene discovery.

Ustilaginomycotina is home to a broad array of fungi including important plant pathogens collectively called smut fungi. Smuts are biotrophs that produce characteristic perennating propagules called teliospores, one of which, Ustilago maydis, is a model genetic organism. Broad exploration of smut biology has been hampered by limited phylogenetic resolution of Ustilaginiomycotina as well as an overall lack of genomic data for members of this subphylum. In this study, we sequenced eight Ustilaginomycotina genomes from previously unrepresented lineages, deciphered ordinal-level phylogenetic relationships for the subphylum, and performed comparative analyses. Unlike other Basidiomycota subphyla, all sampled Ustilaginomycotina genomes are relatively small and compact. Ancestral state reconstruction analyses indicate that teliospore formation was present at the origin of the subphylum. Divergence time estimation dates the divergence of most extant smut fungi after that of grasses (Poaceae). However, we found limited conservation of well-characterized genes related to smut pathogenesis from U. maydis, indicating dissimilar pathogenic mechanisms exist across other smut lineages. The genomes of Malasseziomycetes are highly diverged from the other sampled Ustilaginomycotina, likely due to their unique history as mammal-associated lipophilic yeasts. Despite extensive genomic data, the phylogenetic placement of this class remains ambiguous. Although the sampled Ustilaginomycotina members lack many core enzymes for plant cell wall decomposition and starch catabolism, we identified several novel carbohydrate active enzymes potentially related to pectin breakdown. Finally, ∼50% of Ustilaginomycotina species-specific genes are present in previously undersampled and rare lineages, highlighting the importance of exploring fungal diversity as a resource for novel gene discovery.