Search

Ongoing large-scale genome sequencing projects are forecasting a data deluge that will almost certainly overwhelm current analytical capabilities of evolutionary genomics. In contrast to population genomics, there are no standardized methods in evolutionary genomics for extracting evolutionary and functional (e.g. gene-trait association) signal from genomic data. Here, we examine how current practices of multi-species comparative genomics perform in this aspect and point out that many genomic datasets are under-utilized due to the lack of powerful methodologies. As a result, many current analyses emphasize gene families for which some functional data is already available, resulting in a growing gap between functionally well-characterized genes/organisms and the universe of unknowns. This leaves unknown genes on the 'dark side' of genomes, a problem that will not be mitigated by sequencing more and more genomes, unless we develop tools to infer functional hypotheses for unknown genes in a systematic manner. We provide an inventory of recently developed methods capable of predicting gene-gene and gene-trait associations based on comparative data, then argue that realizing the full potential of whole genome datasets requires the integration of phylogenetic comparative methods into genomics, a rich but underutilized toolbox for looking into the past., (© The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Parasola represents an enigmatic lineage of veil-less, coprinoid fungi in Psathyrellaceae (Agaricales). The species-level taxonomy of the genus has been in a flux recently, resulting in the elimination of some long-established names and the description of new taxa. Here, we reconstruct the phylogeny of Parasola using two nuc rDNA loci, the internal transcribed spacer region (ITS1-5.8S-ITS2) and 28S and identify several putatively undescribed species, of which three are formally described here (Parasola crataegi, P. ochracea, and P. plicatilis-similis) based on molecular and morphological data. Morphological descriptions for the new species and an identification key to accepted Parasola species are given. We revise and discuss our current understanding of the phylogeny of Parasola.

The genomic era has been transformative for many fields, including our understanding of the phylogenetic relationships between organisms. The wide availability of whole-genome sequences practically eliminated data availability as a limiting factor for inferring phylogenetic trees, providing hundreds to thousands of loci for analyses, leading to molecular phylogenetics gradually being replaced by phylogenomics. The new era has also brought new challenges: systematic errors (resulting from, e.g., model violation) can be more pronounced in phylogenomic datasets and can lead to strongly supported incorrect relationships, creating significant incongruence among studies. Here, we review common practices, technical and biological challenges of phylogenomic analyses, with examples illustrated from fungi. We compare major approaches of phylogenetic inference, and illustrate the advantages conferred and challenges presented in phylogenomic case studies across the fungal tree of life, including cases where genome-scale data could conclusively resolve contentious relationships, and others that remain challenging despite the flood of genomic data., (© 2017 Elsevier Inc. All rights reserved.)

Reconciling traditional classifications, morphology, and the phylogenetic relationships of brown-spored agaric mushrooms has proven difficult in many groups, due to extensive convergence in morphological features. Here, we address the monophyly of the Bolbitiaceae, a family with over 700 described species and examine the higher-level relationships within the family using a newly constructed multilocus dataset (ITS, nrLSU rDNA and EF1-alpha). We tested whether the fast-evolving Internal Transcribed Spacer (ITS) sequences can be accurately aligned across the family, by comparing the outcome of two iterative alignment refining approaches (an automated and a manual) and various indel-treatment strategies. We used PRANK to align sequences in both cases. Our results suggest that--although PRANK successfully evades overmatching of gapped sites, referred previously to as alignment overmatching--it infers an unrealistically high number of indel events with natively generated guide-trees. This 'alignment undermatching' could be avoided by using more rigorous (e.g. ML) guide trees. The trees inferred in this study support the monophyly of the core Bolbitiaceae, with the exclusion of Panaeolus, Agrocybe, and some of the genera formerly placed in the family. Bolbitius and Conocybe were found monophyletic, however, Pholiotina and Galerella require redefinition. The phylogeny revealed that stipe coverage type is a poor predictor of phylogenetic relationships, indicating the need for a revision of the intrageneric relationships within Conocybe.

Sections Lanatuli and Atramentarii of the genus Coprinopsis contain some of the best known and most important agaric species, including C. cinerea and C. lagopus, yet a critical, phylogeny-based assessment of the species limits has not been carried out. Taxa have been characterized chiefly on the basis of morphological characters, which however show little discriminatory power and/or considerable overlap between several species pairs. We used ITS and LSU sequence data of 29 described taxa in Coprinopsis sections Lanatuli and Atramentarii to infer species limits and the correspondence between morphological characters and species lineages, as well as to examine the phylogenetic affinities of sections Lanatuli and Atramentarii. Our analyses recovered three large clades, implying a paraphyly for section Lanatuli. Based on morphology and clade structure, we estimate ca. 38 species in the two sections, including several potentially new taxa, three of which are described herein. Coprinopsis pachyderma, C. lagopus var. vacillans, C. acuminata, C. spelaiophila, Coprinus citrinovelatus and Cop. brunneistrangulatus were found to be synonymous with other, earlier described species. Congruent with previous mating studies, our analyses recovered multiple, morphologically indistinguishable lineages within C. lagopus, which included C. lagopus var. vacillans, an ephemeral, developmental variant. Morphological traits supporting the inferred clade structure are discussed. Three new taxa (C. fusispora, C. babosiae, C. villosa), and one new combination (C. mitraespora) are proposed.

Bursts of diversification are known to have contributed significantly to the extant morphological and species diversity, but evidence for many of the theoretical predictions about adaptive radiations have remained contentious. Despite their tremendous diversity, patterns of evolutionary diversification and the contribution of explosive episodes in fungi are largely unknown. Here, using the genus Coprinellus (Psathyrellaceae, Agaricales) as a model, we report the first explosive fungal radiation and infer that the onset of the radiation correlates with a change from a multilayered to a much simpler defense structure on the fruiting bodies. We hypothesize that this change constitutes a key innovation, probably relaxing constraints on diversification imposed by nutritional investment into the development of protective tissues of fruiting bodies. Fossil calibration suggests that Coprinellus mushrooms radiated during the Miocene coinciding with global radiation of large grazing mammals following expansion of dry open grasslands. In addition to diversification rate-based methods, we test the hard polytomy hypothesis, by analyzing the resolvability of internal nodes of the backbone of the putative radiation using Reversible-Jump MCMC. We discuss potential applications and pitfalls of this approach as well as how biologically meaningful polytomies can be distinguished from alignment shortcomings. Our data provide insights into the nature of adaptive radiations in general by revealing a deceleration of morphological diversification through time. The dynamics of morphological diversification was approximated by obtaining the temporal distribution of state changes in discrete traits along the trees and comparing it with the tempo of lineage accumulation. We found that the number of state changes correlate with the number of lineages, even in parts of the tree with short internal branches, and peaks around the onset of the explosive radiation followed by a slowdown, most likely because of the decrease in available niches.

Fruiting body evolution is one of the central topics in fungal evolutionary biology. A number of hypotheses have been developed to explain the contemporary diversity of fruiting body forms, but their evaluation has been hampered by the lack of well-sampled data sets and suitable statistical methods. Phylogenetic evidence of the physiological changes that accompany switches in fruiting body type is lacking, and very little is known about the age of major events of fruiting body evolution. Based on a new multigene phylogeny, by using Bayesian methods, we demonstrate the existence of correlation between a number of morphological features and switches from nondeliquescent to deliquescent (autodigesting) fruiting bodies in the mushroom family Psathyrellaceae. Our results show that switches in the anatomy of two types of spacer cells (cystidia and pseudoparaphyses) and basidia (bimorphic or monomorphic) as well as the structure of the mushroom cap follow the evolution of deliquescent fruiting bodies, which suggests strong functional linkage between these traits. We performed Bayes factor-based tests, referred hereafter to as evolutionary pathway test (EPT), to decide which of the correlated characters were gained first during evolution. The EPTs strongly suggest that deliquescence was gained first, followed after short waiting times by the other morphological features. Bayesian relaxed molecular clock analyses suggest that the various events of switching between fruiting body types occurred independently at various ages during the history of the family. The utility of two mushroom fossils (Archaemarasmius and Protomycena), the only ones with unambiguous taxonomic positions, for the calibration of agaric trees were also examined. Based on our results, we suggest that the evolutionary benefit of deliquescence may be prevention against desiccation via accelerated ontogeny of the fruiting body. Hypotheses regarding the functional significance of the correlated evolution are presented and discussed. Further, we argue that the changes in fruiting body types in mushrooms in general can be attributed to independent events (e.g., dispersal and adaptation) and not to particular geologic ages.

Although the fungal order Mortierellales constitutes one of the largest classical groups of Zygomycota, its phylogeny is poorly understood and no modern taxonomic revision is currently available. In the present study, 90 type and reference strains were used to infer a comprehensive phylogeny of Mortierellales from the sequence data of the complete ITS region and the LSU and SSU genes with a special attention to the monophyly of the genus Mortierella. Out of 15 alternative partitioning strategies compared on the basis of Bayes factors, the one with the highest number of partitions was found optimal (with mixture models yielding the best likelihood and tree length values), implying a higher complexity of evolutionary patterns in the ribosomal genes than generally recognized. Modeling the ITS1, 5.8S, and ITS2, loci separately improved model fit significantly as compared to treating all as one and the same partition. Further, within-partition mixture models suggests that not only the SSU, LSU and ITS regions evolve under qualitatively and/or quantitatively different constraints, but that significant heterogeneity can be found within these loci also. The phylogenetic analysis indicated that the genus Mortierella is paraphyletic with respect to the genera Dissophora, Gamsiella and Lobosporangium and the resulting phylogeny contradict previous, morphology-based sectional classification of Mortierella. Based on tree structure and phenotypic traits, we recognize 12 major clades, for which we attempt to summarize phenotypic similarities. M. longicollis is closely related to the outgroup taxon Rhizopus oryzae, suggesting that it belongs to the Mucorales. Our results demonstrate that traits used in previous classifications of the Mortierellales are highly homoplastic and that the Mortierellales is in a need of a reclassification, where new, phylogenetically informative phenotypic traits should be identified, with molecular phylogenies playing a decisive role.

Recently developed comparative phylogenetic methods offer a wide spectrum of applications in evolutionary biology, although it is generally accepted that their statistical properties are incompletely known. Here, we examine and compare the statistical power of the ML and Bayesian methods with regard to selection of best-fit models of fruiting-body evolution and hypothesis testing of ancestral states on a real-life data set of a physiological trait (autodigestion) in the family Psathyrellaceae. Our phylogenies are based on the first multigene data set generated for the family. Two different coding regimes (binary and multistate) and two data sets differing in taxon sampling density are examined. The Bayesian method outperformed Maximum Likelihood with regard to statistical power in all analyses. This is particularly evident if the signal in the data is weak, i.e. in cases when the ML approach does not provide support to choose among competing hypotheses. Results based on binary and multistate coding differed only modestly, although it was evident that multistate analyses were less conclusive in all cases. It seems that increased taxon sampling density has favourable effects on inference of ancestral states, while model parameters are influenced to a smaller extent. The model best fitting our data implies that the rate of losses of deliquescence equals zero, although model selection in ML does not provide proper support to reject three of the four candidate models. The results also support the hypothesis that non-deliquescence (lack of autodigestion) has been ancestral in Psathyrellaceae, and that deliquescent fruiting bodies represent the preferred state, having evolved independently several times during evolution., (Copyright © 2010 Elsevier Inc. All rights reserved.)

Mycena s.s. is a ubiquitous mushroom genus whose members degrade multiple dead plant substrates and opportunistically invade living plant roots. Having sequenced the nuclear genomes of 24 Mycena species, we find them to defy the expected patterns for fungi based on both their traditionally perceived saprotrophic ecology and substrate specializations. Mycena displayed massive genome expansions overall affecting all gene families, driven by novel gene family emergence, gene duplications, enlarged secretomes encoding polysaccharide degradation enzymes, transposable element (TE) proliferation, and horizontal gene transfers. Mainly due to TE proliferation, Arctic Mycena species display genomes of up to 502 Mbp (2-8× the temperate Mycena), the largest among mushroom-forming Agaricomycetes, indicating a possible evolutionary convergence to genomic expansions sometimes seen in Arctic plants. Overall, Mycena show highly unusual, varied mosaic-like genomic structures adaptable to multiple lifestyles, providing genomic illustration for the growing realization that fungal niche adaptations can be far more fluid than traditionally believed.

Contamination of genomes is an increasingly recognized problem affecting several downstream applications, from comparative evolutionary genomics to metagenomics. Here we introduce ContScout, a precise tool for eliminating foreign sequences from annotated genomes. It achieves high specificity and sensitivity on synthetic benchmark data even when the contaminant is a closely related species, outperforms competing tools, and can distinguish horizontal gene transfer from contamination. A screen of 844 eukaryotic genomes for contamination identified bacteria as the most common source, followed by fungi and plants. Furthermore, we show that contaminants in ancestral genome reconstructions lead to erroneous early origins of genes and inflate gene loss rates, leading to a false notion of complex ancestral genomes. Taken together, we offer here a tool for sensitive removal of foreign proteins, identify and remove contaminants from diverse eukaryotic genomes and evaluate their impact on phylogenomic analyses.

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.

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.

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.

Agaricomycetes are fruiting body-forming fungi that produce some of the most efficient enzyme systems to degrade wood. Despite decades-long interest in their biology, the evolution and functional diversity of both wood-decay and fruiting body formation are incompletely known. We performed comparative genomic and transcriptomic analyses of wood-decay and fruiting body development in Auriculariopsis ampla and Schizophyllum commune (Schizophyllaceae), species with secondarily simplified morphologies, an enigmatic wood-decay strategy and weak pathogenicity to woody plants. The plant cell wall-degrading enzyme repertoires of Schizophyllaceae are transitional between those of white rot species and less efficient wood-degraders such as brown rot or mycorrhizal fungi. Rich repertoires of suberinase and tannase genes were found in both species, with tannases restricted to Agaricomycetes that preferentially colonize bark-covered wood, suggesting potential complementation of their weaker wood-decaying abilities and adaptations to wood colonization through the bark. Fruiting body transcriptomes revealed a high rate of divergence in developmental gene expression, but also several genes with conserved expression patterns, including novel transcription factors and small-secreted proteins, some of the latter which might represent fruiting body effectors. Taken together, our analyses highlighted novel aspects of wood-decay and fruiting body development in an important family of mushroom-forming fungi.

Mushroom-forming fungi (Agaricomycetes) have the greatest morphological diversity and complexity of any group of fungi. They have radiated into most niches and fulfil diverse roles in the ecosystem, including wood decomposers, pathogens or mycorrhizal mutualists. Despite the importance of mushroom-forming fungi, large-scale patterns of their evolutionary history are poorly known, in part due to the lack of a comprehensive and dated molecular phylogeny. Here, using multigene and genome-based data, we assemble a 5,284-species phylogenetic tree and infer ages and broad patterns of speciation/extinction and morphological innovation in mushroom-forming fungi. Agaricomycetes started a rapid class-wide radiation in the Jurassic, coinciding with the spread of (sub)tropical coniferous forests and a warming climate. A possible mass extinction, several clade-specific adaptive radiations and morphological diversification of fruiting bodies followed during the Cretaceous and the Paleogene, convergently giving rise to the classic toadstool morphology, with a cap, stalk and gills (pileate-stipitate morphology). This morphology is associated with increased rates of lineage diversification, suggesting it represents a key innovation in the evolution of mushroom-forming fungi. The increase in mushroom diversity started during the Mesozoic-Cenozoic radiation event, an era of humid climate when terrestrial communities dominated by gymnosperms and reptiles were also expanding.

Fungal decomposition of plant cell walls (PCW) is a complex process that has diverse industrial applications and huge impacts on the carbon cycle. White rot (WR) is a powerful mode of PCW decay in which lignin and carbohydrates are both degraded. Mechanistic studies of decay coupled with comparative genomic analyses have provided clues to the enzymatic components of WR systems and their evolutionary origins, but the complete suite of genes necessary for WR remains undetermined. Here, we use phylogenomic comparative methods, which we validate through simulations, to identify shifts in gene family diversification rates that are correlated with evolution of WR, using data from 62 fungal genomes. We detected 409 gene families that appear to be evolutionarily correlated with WR. The identified gene families encode well-characterized decay enzymes, e.g., fungal class II peroxidases and cellobiohydrolases, and enzymes involved in import and detoxification pathways, as well as 73 gene families that have no functional annotation. About 310 of the 409 identified gene families are present in the genome of the model WR fungus Phanerochaete chrysosporium and 192 of these (62%) have been shown to be upregulated under ligninolytic culture conditions, which corroborates the phylogeny-based functional inferences. These results illuminate the complexity of WR and suggest that its evolution has involved a general elaboration of the decay apparatus, including numerous gene families with as-yet unknown exact functions.

Evolution of lignocellulose decomposition was one of the most ecologically important innovations in fungi. White-rot fungi in the Agaricomycetes (mushrooms and relatives) are the most effective microorganisms in degrading both cellulose and lignin components of woody plant cell walls (PCW). However, the precise evolutionary origins of lignocellulose decomposition are poorly understood, largely because certain early-diverging clades of Agaricomycetes and its sister group, the Dacrymycetes, have yet to be sampled, or have been undersampled, in comparative genomic studies. Here, we present new genome sequences of ten saprotrophic fungi, including members of the Dacrymycetes and early-diverging clades of Agaricomycetes (Cantharellales, Sebacinales, Auriculariales, and Trechisporales), which we use to refine the origins and evolutionary history of the enzymatic toolkit of lignocellulose decomposition. We reconstructed the origin of ligninolytic enzymes, focusing on class II peroxidases (AA2), as well as enzymes that attack crystalline cellulose. Despite previous reports of white rot appearing as early as the Dacrymycetes, our results suggest that white-rot fungi evolved later in the Agaricomycetes, with the first class II peroxidases reconstructed in the ancestor of the Auriculariales and residual Agaricomycetes. The exemplars of the most ancient clades of Agaricomycetes that we sampled all lack class II peroxidases, and are thus concluded to use a combination of plesiomorphic and derived PCW degrading enzymes that predate the evolution of white rot.

Abstract: Aims: A homology modeling methodology was developed and used to obtain the 3D structures for four putative catalases of Rhizopus oryzae in order to assess their functionality. Main methods: Homology models were built using different modeling strategies using non-protein compounds as steric constraints, a symmetry constraint to force identical chains and an additional loop modeling algorithm. Percent structural overlap values (SO) were calculated for each model–template pair to qualify the homology models. Key findings: Comparing the different modeling strategies by the SO values revealed that the quality of the models, i.e. the similarity to the template was greatly increased in the presence of the prosthetic groups, modeling multiple protein chains together, enforcing symmetrical chains and applying additional loop modeling. For the best homology models achieved this way, the SO values express similar evolutionary relationships between the proteins modeled and the templates that were previously established by phylogenetic analysis. In three out of the four catalases of R. oryzae the highest quality models, the active center, i.e. the heme molecule and the surrounding amino acids showed a spatial arrangement identical to that observed experimentally in other catalases. The remaining protein is missing an 11 residue long fragment and has mutated residues within the active center. Significance: Better homology models can be obtained with templates chosen by phylogenetic relationship, alth-ough building an accurate model needs structural constraints too. Calculating the structural overlap between the models and the templates may also help to find the appropriate templates. [Copyright &y& Elsevier]