Horizontal Gene Transfer as an Indispensable Driver for Evolution of Neocallimastigomycota into a Distinct Gut-Dwelling Fungal Lineage

INTERNATIONAL SYMPOSIUM ON RUMINANT PHYSIOLOGY: Rumen fungi, archaea and their interactions

Anaerobic gut fungi (AGF) were the last phylum to be identified within the rumen microbiome and account for 7-9% of microbial biomass. They produce potent lignocellulases that degrade recalcitrant plant cell walls, and rhizoids that can penetrate the cuticle of plant cells, exposing internal components to other microbiota. Interspecies H2 transfer between AGF and rumen methanogenic archaea is an essential metabolic process in the rumen that occurs during the reduction of CO2 to CH4 by methanogens. This symbiotic relationship is bolstered by hydrogensomes, fungal organelles that generate H2 and formate. Interspecies H2 transfer prevents the accumulation of reducing equivalents that would otherwise impede fermentation. The extent to which hydrogenosomes serve as a conduit for H2 flow to methanogens is unknown, but it is likely greater with low quality forages. Strategies that alter the production of CH4 could also have implications for H2 transfer by anaerobic fungi. Understanding the factors that drive these interactions and H2 flow could provide insight into the effect of reducing CH4 production on the activity of ruminal fungi and the digestion of low-quality feeds.

Recent advances in methane and hydrogen production from lignocellulosic degradation with anaerobic fungi

Anaerobic fungi (AF) efficiently degrade lignocellulosic biomass with unique pseudoroot system and enzymatic properties that can remove polysaccharides and some lignified components from plant cell walls, further releasing acetate, lactate, ethanol, hydrogen (H2), etc. As research on AF for bioengineering has become a hot topic, a review of lignocellulosic conversion with AF for methane (CH4) and H2 production is needed. Efficient degradation of lignocellulose with AF mainly relies on multiple free carbohydrate-active enzymes and cellulosomes in the free and bound state. Meanwhile, co-cultivation of AF and methanogens significantly improves the lignocellulose degradation and CH4 production, and the maximum CH4 yield reached 315 mL/g. Bioaugmentation of AF in anaerobic digestion increases the maximum CH4 yield by 330 %. Also, AF show H2 production potential, however, H2 yield from anaerobic fungal fermentation of lignocellulose remains low. Therefore, anaerobic fungi have great potential in the conversion of lignocellulosic biomass to CH4 and H2.

Genetic conflicts in budding yeast: The 2μ plasmid as a model selfish element

Antagonistic coevolution, arising from genetic conflict, can drive rapid evolution and biological innovation. Conflict can arise both between organisms and within genomes. This review focuses on budding yeasts as a model system for exploring intra- and inter-genomic genetic conflict, highlighting in particular the 2-micron (2μ) plasmid as a model selfish element. The 2μ is found widely in laboratory strains and industrial isolates of Saccharomyces cerevisiae and has long been known to cause host fitness defects. Nevertheless, the plasmid is frequently ignored in the context of genetic, fitness, and evolution studies. Here, I make a case for further exploring the evolutionary impact of the 2μ plasmid as well as other selfish elements of budding yeasts, discuss recent advances, and, finally, future directions for the field.

Symbioses between fungi and bacteria: from mechanisms to impacts on biodiversity

Symbiotic interactions between fungi and bacteria range from positive to negative. They are ubiquitous in free-living as well as host-associated microbial communities worldwide. Yet, the impact of fungal-bacterial symbioses on the organization and dynamics of microbial communities is uncertain. There are two reasons for this uncertainty: (1) knowledge gaps in the understanding of the genetic mechanisms underpinning fungal-bacterial symbioses and (2) prevailing interpretations of ecological theory that favor antagonistic interactions as drivers stabilizing biological communities despite the existence of models emphasizing contributions of positive interactions. This review synthesizes information on fungal-bacterial symbioses common in the free-living microbial communities of the soil as well as in host-associated polymicrobial biofilms. The interdomain partnerships are considered in the context of the relevant community ecology models, which are discussed critically.

Structural, Biochemical, and Phylogenetic Analysis of Bacterial and Fungal Carbohydrate Esterase Family 15 Glucuronoyl Esterases in the Rumen

Glucuronoyl esterases (GEs) are carbohydrate active enzymes in carbohydrate esterase family 15 which are involved in the hydrolysis of lignin-carbohydrate complexes. They are encoded by a wide range of aerobic and anaerobic fungi and bacteria inhabiting diverse environments. The rumen microbiome is a complex microbial community with a wide array of enzymes that specialize in deconstructing plant cell wall carbohydrates. Enzymes from the rumen tend to show low similarity to homologues found in other environments, making the rumen microbiome a promising source for the discovery of novel enzymes. Using a combination of phylogenetic and structural analysis, we investigated the structure-function relationship of GEs from the rumen bacteria Fibrobacter succinogenes and Ruminococcus flavefaciens, and from the rumen fungus, Piromyces rhizinflata. All adopt a canonical α/β hydrolase fold and possess a structurally conserved Ser-His-Glu/Asp catalytic triad. Structural variations in the enzymes are localized to loops surrounding the active site. Analysis of the active site structures in these enzymes emphasized the importance of structural plasticity in GEs with non-canonical active site conformations. We hypothesize that interkingdom HGT events may have contributed to the diversity of GEs in the rumen, and this is demonstrated by the phylogenetic and structural similarity observed between rumen bacterial and fungal GEs. This study advances our understanding of the structure-function relationship in glucuronoyl esterases and illuminates the evolutionary dynamics that contribute to enzyme diversity in the rumen microbiome.

Anaerobic fungi contain abundant, diverse, and transcriptionally active Long Terminal Repeat retrotransposons

Long Terminal Repeat (LTR) retrotransposons are a class of repetitive elements that are widespread in the genomes of plants and many fungi. LTR retrotransposons have been associated with rapidly evolving gene clusters in plants and virulence factor transfer in fungal-plant parasite-host interactions. We report here the abundance and transcriptional activity of LTR retrotransposons across several species of the early-branching Neocallimastigomycota, otherwise known as the anaerobic gut fungi (AGF). The ubiquity of LTR retrotransposons in these genomes suggests key evolutionary roles in these rumen-dwelling biomass degraders, whose genomes also contain many enzymes that are horizontally transferred from other rumen-dwelling prokaryotes. Up to 10% of anaerobic fungal genomes consist of LTR retrotransposons, and the mapping of sequences from LTR retrotransposons to transcriptomes shows that the majority of clusters are transcribed, with some exhibiting expression greater than 104 reads per kilobase million mapped reads (rpkm). Many LTR retrotransposons are strongly differentially expressed upon heat stress during fungal cultivation, with several exhibiting a nearly three-log10 fold increase in expression, whereas growth substrate variation modulated transcription to a lesser extent. We show that some LTR retrotransposons contain carbohydrate-active enzymes (CAZymes), and the expansion of CAZymes within genomes and among anaerobic fungal species may be linked to retrotransposon activity. We further discuss how these widespread sequences may be a source of promoters and other parts towards the bioengineering of anaerobic fungi.

Horizontal gene transfer in eukaryotes: aligning theory with data

Horizontal gene transfer (HGT), or lateral gene transfer, is the non-sexual movement of genetic information between genomes. It has played a pronounced part in bacterial and archaeal evolution, but its role in eukaryotes is less clear. Behaviours unique to eukaryotic cells - phagocytosis and endosymbiosis - have been proposed to increase the frequency of HGT, but nuclear genomes encode fewer HGTs than bacteria and archaea. Here, I review the existing theory in the context of the growing body of data on HGT in eukaryotes, which suggests that any increased chance of acquiring new genes through phagocytosis and endosymbiosis is offset by a reduced need for these genes in eukaryotes, because selection in most eukaryotes operates on variation not readily generated by HGT.

Anaerobic fungi in the tortoise alimentary tract illuminate early stages of host-fungal symbiosis and Neocallimastigomycota evolution

Anaerobic gut fungi (AGF, Neocallimastigomycota) reside in the alimentary tract of herbivores. While their presence in mammals is well documented, evidence for their occurrence in non-mammalian hosts is currently sparse. Culture-independent surveys of AGF in tortoises identified a unique community, with three novel deep-branching genera representing >90% of sequences in most samples. Representatives of all genera were successfully isolated under strict anaerobic conditions. Transcriptomics-enabled phylogenomic and molecular dating analyses indicated an ancient, deep-branching position in the AGF tree for these genera, with an evolutionary divergence time estimate of 104-112 million years ago (Mya). Such estimates push the establishment of animal-Neocallimastigomycota symbiosis from the late to the early Cretaceous. Further, tortoise-associated isolates (T-AGF) exhibited limited capacity for plant polysaccharides metabolism and lacked genes encoding several carbohydrate-active enzyme (CAZyme) families. Finally, we demonstrate that the observed curtailed degradation capacities and reduced CAZyme repertoire is driven by the paucity of horizontal gene transfer (HGT) in T-AGF genomes, compared to their mammalian counterparts. This reduced capacity was reflected in an altered cellulosomal production capacity in T-AGF. Our findings provide insights into the phylogenetic diversity, ecological distribution, evolutionary history, evolution of fungal-host nutritional symbiosis, and dynamics of genes acquisition in Neocallimastigomycota.

Herptile gut microbiomes: a natural system to study multi-kingdom interactions between filamentous fungi and bacteria

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.

Rumen microbes, enzymes, metabolisms, and application in lignocellulosic waste conversion - A comprehensive review

The rumen of ruminants is a natural anaerobic fermentation system that efficiently degrades lignocellulosic biomass and mainly depends on synergistic interactions between multiple microbes and their secreted enzymes. Ruminal microbes have been employed as biomass waste converters and are receiving increasing attention because of their degradation performance. To explore the application of ruminal microbes and their secreted enzymes in biomass waste, a comprehensive understanding of these processes is required. Based on the degradation capacity and mechanism of ruminal microbes and their secreted lignocellulose enzymes, this review concentrates on elucidating the main enzymatic strategies that ruminal microbes use for lignocellulose degradation, focusing mainly on polysaccharide metabolism-related gene loci and cellulosomes. Hydrolysis, acidification, methanogenesis, interspecific H2 transfer, and urea cycling in ruminal metabolism are also discussed. Finally, we review the research progress on the conversion of biomass waste into biofuels (bioethanol, biohydrogen, and biomethane) and value-added chemicals (organic acids) by ruminal microbes. This review aims to provide new ideas and methods for ruminal microbe and enzyme applications, biomass waste conversion, and global energy shortage alleviation.

Sulfur assimilation using gaseous carbonyl sulfide by the soil fungus Trichoderma harzianum

Fungi have the capacity to assimilate a diverse range of both inorganic and organic sulfur compounds. It has been recognized that all sulfur sources taken up by fungi are in soluble forms. In this study, we present evidence that fungi can utilize gaseous carbonyl sulfide (COS) for the assimilation of a sulfur compound. We found that the filamentous fungus Trichoderma harzianum strain THIF08, which has constitutively high COS-degrading activity, was able to grow with COS as the sole sulfur source. Cultivation with 34S-labeled COS revealed that sulfur atom from COS was incorporated into intracellular metabolites such as glutathione and ergothioneine. COS degradation by strain THIF08, in which as much of the moisture derived from the agar medium as possible was removed, indicated that gaseous COS was taken up directly into the cell. Escherichia coli transformed with a COS hydrolase (COSase) gene, which is clade D of the β-class carbonic anhydrase subfamily enzyme with high specificity for COS but low activity for CO2 hydration, showed that the COSase is involved in COS assimilation. Comparison of sulfur metabolites of strain THIF08 revealed a higher relative abundance of reduced sulfur compounds under the COS-supplemented condition than the sulfate-supplemented condition, suggesting that sulfur assimilation is more energetically efficient with COS than with sulfate because there is no redox change of sulfur. Phylogenetic analysis of the genes encoding COSase, which are distributed in a wide range of fungal taxa, suggests that the common ancestor of Ascomycota, Basidiomycota, and Mucoromycota acquired COSase at about 790-670 Ma.IMPORTANCEThe biological assimilation of gaseous CO2 and N2 involves essential processes known as carbon fixation and nitrogen fixation, respectively. In this study, we found that the fungus Trichoderma harzianum strain THIF08 can grow with gaseous carbonyl sulfide (COS), the most abundant and ubiquitous gaseous sulfur compound, as a sulfur source. When the fungus grew in these conditions, COS was assimilated into sulfur metabolites, and the key enzyme of this assimilation process is COS hydrolase (COSase), which specifically degrades COS. Moreover, the pathway was more energy efficient than the typical sulfate assimilation pathway. COSase genes are widely distributed in Ascomycota, Basidiomycota, and Mucoromycota and also occur in some Chytridiomycota, indicating that COS assimilation is widespread in fungi. Phylogenetic analysis of these genes revealed that the acquisition of COSase in filamentous fungi was estimated to have occurred at about 790-670 Ma, around the time that filamentous fungi transitioned to a terrestrial environment.

ContScout: sensitive detection and removal of contamination from annotated genomes

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.

Comparative genomics reveals probable adaptations for xylose use in Thermoanaerobacterium saccharolyticum

Second-generation ethanol, a promising biofuel for reducing greenhouse gas emissions, faces challenges due to the inefficient metabolism of xylose, a pentose sugar. Overcoming this hurdle requires exploration of genes, pathways, and organisms capable of fermenting xylose. Thermoanaerobacterium saccharolyticum is an organism capable of naturally fermenting compounds of industrial interest, such as xylose, and understanding evolutionary adaptations may help to bring novel genes and information that can be used for industrial yeast, increasing production of current bio-platforms. This study presents a deep evolutionary study of members of the firmicutes clade, focusing on adaptations in Thermoanaerobacterium saccharolyticum that may be related to overall fermentation metabolism, especially for xylose fermentation. One highlight is the finding of positive selection on a xylose-binding protein of the xylFGH operon, close to the annotated sugar binding site, with this protein already being found to be expressed in xylose fermenting conditions in a previous study. Results from this study can serve as basis for searching for candidate genes to use in industrial strains or to improve Thermoanaerobacterium saccharolyticum as a new microbial cell factory, which may help to solve current problems found in the biofuels' industry.

Genome-wide analysis of transmembrane 9 superfamily genes in wheat (Triticum aestivum) and their expression in the roots under nitrogen limitation and Bacillus amyloliquefaciens PDR1 treatment conditions

Introduction

Transmembrane 9 superfamily (TM9SF) proteins play significant roles in plant physiology. However, these proteins are poorly characterized in wheat (Triticum aestivum). The present study aimed at the genome-wide analysis of putative wheat TM9SF (TraesTM9SF) proteins and their potential involvement in response to nitrogen limitation and Bacillus amyloliquefaciens PDR1 treatments.

Methods

TraesTM9SF genes were retrieved from the wheat genome, and their physiochemical properties, alignment, phylogenetic, motif structure, cis-regulatory element, synteny, protein-protein interaction (PPI), and transcription factor (TF) prediction analyses were performed. Transcriptome sequencing and quantitative real-time polymerase reaction (qRT-PCR) were performed to detect gene expression in roots under single or combined treatments with nitrogen limitation and B. amyloliquefaciens PDR1.

Results and discussion

Forty-seven TraesTM9SF genes were identified in the wheat genome, highlighting the significance of these genes in wheat. TraesTM9SF genes were absent on some wheat chromosomes and were unevenly distributed on the other chromosomes, indicating that potential regulatory functions and evolutionary events may have shaped the TraesTM9SF gene family. Fifty-four cis-regulatory elements, including light-response, hormone response, biotic/abiotic stress, and development cis-regulatory elements, were present in the TraesTM9SF promoter regions. No duplication of TraesTM9SF genes in the wheat genome was recorded, and 177 TFs were predicted to target the 47 TraesTM9SF genes in a complex regulatory network. These findings offer valued data for predicting the putative functions of uncharacterized TM9SF genes. Moreover, transcriptome analysis and validation by qRT-PCR indicated that the TraesTM9SF genes are expressed in the root system of wheat and are potentially involved in the response of this plant to single or combined treatments with nitrogen limitation and B. amyloliquefaciens PDR1, suggesting their functional roles in plant growth, development, and stress responses.

Conclusion

These findings may be vital in further investigation of the function and biological applications of TM9SF genes in wheat.

The evolution of anaerobic growth in Saccharomycotina yeasts

Humans rely on the ability of budding yeasts to grow without oxygen in industrial scale fermentations that produce beverages, foods, and biofuels. Oxygen is deeply woven into the energy metabolism and biosynthetic capabilities of budding yeasts. While diverse ecological habitats may provide wide varieties of different carbon and nitrogen sources for yeasts to utilize, there is no direct substitute for molecular oxygen, only a range of availability. Understanding how a small subset of budding yeasts evolved the ability to grow without oxygen could expand the set of useful species in industrial scale fermentations as well as provide insight into the cryptic field of yeast ecology. However, we still do not yet appreciate the full breadth of species that can growth without oxygen, what genes underlie this adaptation, and how these genes have evolved.

Identification and expression of small multidrug resistance transporters in early-branching anaerobic fungi

Membrane-embedded transporters impart essential functions to cells as they mediate sensing and the uptake and extrusion of nutrients, waste products, and effector molecules. Promiscuous multidrug exporters are implicated in resistance to drugs and antibiotics and are highly relevant for microbial engineers who seek to enhance the tolerance of cell factory strains to hydrophobic bioproducts. Here, we report on the identification of small multidrug resistance (SMR) transporters in early-branching anaerobic fungi (Neocallimastigomycetes). The SMR class of transporters is commonly found in bacteria but has not previously been reported in eukaryotes. In this study, we show that SMR transporters from anaerobic fungi can be produced heterologously in the model yeast Saccharomyces cerevisiae, demonstrating the potential of these proteins as targets for further characterization. The discovery of these novel anaerobic fungal SMR transporters offers a promising path forward to enhance bioproduction from engineered microbial strains.

Genomes of fungi and relatives reveal delayed loss of ancestral gene families and evolution of key fungal traits

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.

Patterns and determinants of the global herbivorous mycobiome

Despite their role in host nutrition, the anaerobic gut fungal (AGF) component of the herbivorous gut microbiome remains poorly characterized. Here, to examine global patterns and determinants of AGF diversity, we generate and analyze an amplicon dataset from 661 fecal samples from 34 mammalian species, 9 families, and 6 continents. We identify 56 novel genera, greatly expanding AGF diversity beyond current estimates (31 genera and candidate genera). Community structure analysis indicates that host phylogenetic affiliation, not domestication status and biogeography, shapes the community rather than. Fungal-host associations are stronger and more specific in hindgut fermenters than in foregut fermenters. Transcriptomics-enabled phylogenomic and molecular clock analyses of 52 strains from 14 genera indicate that most genera with preferences for hindgut hosts evolved earlier (44-58 Mya) than those with preferences for foregut hosts (22-32 Mya). Our results greatly expand the documented scope of AGF diversity and provide an ecologically and evolutionary-grounded model to explain the observed patterns of AGF diversity in extant animal hosts.

Phylogenomic analysis of the Neocallimastigomycota: proposal of Caecomycetaceae fam. nov., Piromycetaceae fam. nov., and emended description of the families Neocallimastigaceae and Anaeromycetaceae

The anaerobic gut fungi (AGF) represent a coherent phylogenetic clade within the Mycota. Twenty genera have been described so far. Currently, the phylogenetic and evolutionary relationships between AGF genera remain poorly understood. Here, we utilized 52 transcriptomic datasets from 14 genera to resolve AGF inter-genus relationships using phylogenomics, and to provide a quantitative estimate (amino acid identity, AAI) for intermediate rank assignments. We identify four distinct supra-genus clades, encompassing all genera producing polyflagellated zoospores, bulbous rhizoids, the broadly circumscribed genus Piromyces, and the Anaeromyces and affiliated genera. We also identify the genus Khoyollomyces as the earliest evolving AGF genus. Concordance between phylogenomic outputs and RPB1 and D1/D2 LSU, but not RPB2, MCM7, EF1α or ITS1, phylogenies was observed. We combine phylogenomic analysis and AAI outputs with informative phenotypic traits to propose accommodating 14/20 AGF genera into four families: Caecomycetaceae fam. nov. (encompassing the genera Caecomyces and Cyllamyces), Piromycetaceae fam. nov. (encompassing the genus Piromyces), emend the description of the family Neocallimastigaceae to encompass the genera Neocallimastix, Orpinomyces, Pecoramyces, Feramyces, Ghazallomyces, Aestipascuomyces and Paucimyces, as well as the family Anaeromycetaceae to include the genera Oontomyces, Liebetanzomyces and Capellomyces in addition to Anaeromyces. We refrain from proposing families for the deeply branching genus Khoyollomyces and for genera with uncertain position (Buwchfawromyces, Joblinomyces, Tahromyces, Agriosomyces and Aklioshbomyces) pending availability of additional isolates and sequence data; and these genera are designated as 'genera incertae sedis' in the order Neocallimastigales. Our results establish an evolutionary-grounded Linnaean taxonomic framework for the AGF, provide quantitative estimates for rank assignments, and demonstrate the utility of RPB1 as an additional informative marker in Neocallimastigomycota taxonomy.

Comparative genomic study of the Penicillium genus elucidates a diverse pangenome and 15 lateral gene transfer events

The Penicillia are known to produce a wide range natural products-some with devastating outcome for the agricultural industry and others with unexploited potential in different applications. However, a large-scale overview of the biosynthetic potential of different species has been lacking. In this study, we sequenced 93 Penicillium isolates and, together with eleven published genomes that hold similar assembly characteristics, we established a species phylogeny as well as defining a Penicillium pangenome. A total of 5612 genes were shared between ≥ 98 isolates corresponding to approximately half of the average number of genes a Penicillium genome holds. We further identified 15 lateral gene transfer events that have occurred in this collection of Penicillium isolates, which might have played an important role, such as niche adaption, in the evolution of these fungi. The comprehensive characterization of the genomic diversity in the Penicillium genus supersedes single-reference genomes, which do not necessarily capture the entire genetic variation.

Definition of the Acceptor Substrate Binding Specificity in Plant Xyloglucan Endotransglycosylases Using Computational Chemistry

Xyloglucan endotransglycosylases (XETs) play key roles in the remodelling and reconstruction of plant cell walls. These enzymes catalyse homo-transglycosylation reactions with xyloglucan-derived donor and acceptor substrates and hetero-transglycosylation reactions with a variety of structurally diverse polysaccharides. In this work, we describe the basis of acceptor substrate binding specificity in non-specific Tropaeolum majus (TmXET6.3) and specific Populus tremula x tremuloides (PttXET16A) XETs, using molecular docking and molecular dynamics (MD) simulations combined with binding free energy calculations. The data indicate that the enzyme-donor (xyloglucan heptaoligosaccharide or XG-OS7)/acceptor complexes with the linear acceptors, where a backbone consisted of glucose (Glc) moieties linked via (1,4)- or (1,3)-β-glycosidic linkages, were bound stably in the active sites of TmXET6.3 and PttXET16A. Conversely, the acceptors with the (1,6)-β-linked Glc moieties were bound stably in TmXET6.3 but not in PttXET16A. When in the (1,4)-β-linked Glc containing acceptors, the saccharide moieties were replaced with mannose or xylose, they bound stably in TmXET6.3 but lacked stability in PttXET16A. MD simulations of the XET-donor/acceptor complexes with acceptors derived from (1,4;1,3)-β-glucans highlighted the importance of (1,3)-β-glycosidic linkages and side chain positions in the acceptor substrates. Our findings explain the differences in acceptor binding specificity between non-specific and specific XETs and associate theoretical to experimental data.

Simultaneous Metabarcoding and Quantification of Neocallimastigomycetes from Environmental Samples: Insights into Community Composition and Novel Lineages

Anaerobic fungi from the herbivore digestive tract (Neocallimastigomycetes) are primary lignocellulose modifiers and hold promise for biotechnological applications. Their molecular detection is currently difficult due to the non-specificity of published primer pairs, which impairs evolutionary and ecological research with environmental samples. We developed and validated a Neocallimastigomycetes-specific PCR primer pair targeting the D2 region of the ribosomal large subunit suitable for screening, quantifying, and sequencing. We evaluated this primer pair in silico on sequences from all known genera, in vitro with pure cultures covering 16 of the 20 known genera, and on environmental samples with highly diverse microbiomes. The amplified region allowed phylogenetic differentiation of all known genera and most species. The amplicon is about 350 bp long, suitable for short-read high-throughput sequencing as well as qPCR assays. Sequencing of herbivore fecal samples verified the specificity of the primer pair and recovered highly diverse and so far unknown anaerobic gut fungal taxa. As the chosen barcoding region can be easily aligned and is taxonomically informative, the sequences can be used for classification and phylogenetic inferences. Several new Neocallimastigomycetes clades were obtained, some of which represent putative novel lineages such as a clade from feces of the rodent Dolichotis patagonum (mara).

Taxonomy of the anaerobic gut fungi (Neocallimastigomycota): a review of classification criteria and description of current taxa

Members of the anaerobic gut fungi (Neocallimastigomycota) reside in the rumen and alimentary tract of larger mammalian and some reptilian, marsupial and avian herbivores. The recent decade has witnessed a significant expansion in the number of described Neocallimastigomycota genera and species. However, the difficulties associated with the isolation and maintenance of Neocallimastigomycota strains has greatly complicated comparative studies to resolve inter- and intra-genus relationships. Here, we provide an updated outline of Neocallimastigomycota taxonomy. We critically evaluate various morphological, microscopic and phylogenetic traits previously and currently utilized in Neocallimastigomycota taxonomy, and provide an updated key for quick characterization of all genera. We then synthesize data from taxa description manuscripts, prior comparative efforts and molecular sequence data to present an updated list of Neocallimastigomycota genera and species, with an emphasis on resolving relationships and identifying synonymy between recent and historic strains. We supplement data from published manuscripts with information and illustrations from strains in the authors' collections. Twenty genera and 36 species are recognized, but the status of 10 species in the genera Caecomyces, Piromyces, Anaeromyces and Cyllamyces remains uncertain due to the unavailability of culture and conferre (cf.) strains, lack of sequence data, and/or inadequacy of available microscopic and phenotypic data. Six cases of synonymy are identified in the genera Neocallimastix and Caecomyces, and two names in the genus Piromyces are rejected based on apparent misclassification.

An Integrative View of the Phyllosphere Mycobiome of Native Rubber Trees in the Brazilian Amazon

The rubber tree, Hevea brasiliensis, is a neotropical Amazonian species. Despite its high economic value and fungi associated with native individuals, in its original area in Brazil, it has been scarcely investigated and only using culture-dependent methods. Herein, we integrated in silico approaches with novel field/experimental approaches and a case study of shotgun metagenomics and small RNA metatranscriptomics of an adult individual. Scientific literature, host fungus, and DNA databases are biased to fungal taxa, and are mainly related to rubber tree diseases and in non-native ecosystems. Metabarcoding retrieved specific phyllospheric core fungal communities of all individuals, adults, plantlets, and leaves of the same plant, unravelling hierarchical structured core mycobiomes. Basidiomycotan yeast-like fungi that display the potential to produce antifungal compounds and a complex of non-invasive ectophytic parasites (Sooty Blotch and Flyspeck fungi) co-occurred in all samples, encompassing the strictest core mycobiome. The case study of the same adult tree (previously studied using culture-dependent approach) analyzed by amplicon, shotgun metagenomics, and small RNA transcriptomics revealed a high relative abundance of insect parasite-pathogens, anaerobic fungi and a high expression of Trichoderma (a fungal genus long reported as dominant in healthy wild rubber trees), respectively. Altogether, our study unravels new and intriguing information/hypotheses of the foliar mycobiome of native H. brasiliensis, which may also occur in other native Amazonian trees.

Disentangling the Evolutionary History of Feo, the Major Ferrous Iron Transport System in Bacteria

Iron acquisition is essential for almost all living organisms. In certain environments, ferrous iron is the most prevalent form of this element. Feo is the most widespread system for ferrous iron uptake in bacteria and is critical for virulence in some species. The canonical architecture of Feo consists of a large transmembrane nucleoside triphosphatase (NTPase) protein, FeoB, and two accessory cytoplasmic proteins, FeoA and FeoC. The role of the latter components and the mechanism by which Feo orchestrates iron transport are unclear. In this study, we conducted a comparative analysis of Feo protein sequences to gain insight into the evolutionary history of this transporter. We identified instances of how horizontal gene transfer contributed to the evolution of Feo. Also, we found that FeoC, while absent in most lineages, is largely present in the Gammaproteobacteria group, although its sequence is poorly conserved. We propose that FeoC, which may couple FeoB NTPase activity with pore opening, was an ancestral element that has been dispensed with through mutations in FeoA and FeoB in some lineages. We provide experimental evidence supporting this hypothesis by isolating and characterizing FeoC-independent mutants of the Vibrio cholerae Feo system. Also, we confirmed that the closely related species Shewanella oneidensis does not require FeoC; thus, Vibrio FeoC sequences may resemble transitional forms on an evolutionary pathway toward FeoC-independent transporters. Finally, by combining data from our bioinformatic analyses with this experimental evidence, we propose an evolutionary model for the Feo system in bacteria. IMPORTANCE Feo, a ferrous iron transport system composed of three proteins (FeoA, -B, and -C), is the most prevalent bacterial iron transporter. It plays an important role in iron acquisition in low-oxygen environments and some host-pathogen interactions. The large transmembrane protein FeoB provides the channel for the transport of iron into the bacterial cell, but the functions of the two small, required accessory proteins FeoA and FeoC are not well understood. Analysis of the evolution of this transporter shows that FeoC is poorly conserved and has been lost from many bacterial lineages. Experimental evidence indicates that FeoC may have different functions in different species that retain this protein, and the loss of FeoC is promoted by mutations in FeoA or by the fusion of FeoA and FeoB.

In Silico Predictions of Ecological Plasticity Mediated by Protein Family Expansions in Early-Diverging Fungi

Early-diverging fungi (EDF) are ubiquitous and versatile. Their diversity is reflected in their genome sizes and complexity. For instance, multiple protein families have been reported to expand or disappear either in particular genomes or even whole lineages. The most commonly mentioned are CAZymes (carbohydrate-active enzymes), peptidases and transporters that serve multiple biological roles connected to, e.g., metabolism and nutrients intake. In order to study the link between ecology and its genomic underpinnings in a more comprehensive manner, we carried out a systematic in silico survey of protein family expansions and losses among EDF with diverse lifestyles. We found that 86 protein families are represented differently according to EDF ecological features (assessed by median count differences). Among these there are 19 families of proteases, 43 CAZymes and 24 transporters. Some of these protein families have been recognized before as serine and metallopeptidases, cellulases and other nutrition-related enzymes. Other clearly pronounced differences refer to cell wall remodelling and glycosylation. We hypothesize that these protein families altogether define the preliminary fungal adaptasome. However, our findings need experimental validation. Many of the protein families have never been characterized in fungi and are discussed in the light of fungal ecology for the first time.

Cryogenian Origin and Subsequent Diversification of the Plant Cell-Wall Enzyme XTH Family

All land plants encode large multigene families of xyloglucan endotransglucosylase/hydrolases (XTHs), plant-specific enzymes that cleave and reconnect plant cell-wall polysaccharides. Despite the ubiquity of these enzymes, considerable uncertainty remains regarding the evolutionary history of the XTH family. Phylogenomic and comparative analyses in this study traced the non-plant origins of the XTH family to Alphaproteobacteria ExoKs, bacterial enzymes involved in loosening biofilms, rather than Firmicutes licheninases, plant biomass digesting enzymes, as previously supposed. The relevant horizontal gene transfer (HGT) event was mapped to the divergence of non-swimming charophycean algae in the Cryogenian geological period. This HGT event was the likely origin of charophycean EG16-2s, which are putative intermediates between ExoKs and XTHs. Another HGT event in the Cryogenian may have led from EG16-2s or ExoKs to fungal Congo Red Hypersensitive proteins (CRHs) to fungal CRHs, enzymes that cleave and reconnect chitin and glucans in fungal cell walls. This successive transfer of enzyme-encoding genes may have supported the adaptation of plants and fungi to the ancient icy environment by facilitating their sessile lifestyles. Furthermore, several protein evolutionary steps, including coevolution of substrate-interacting residues and putative intra-family gene fusion, occurred in the land plant lineage and drove diversification of the XTH family. At least some of those events correlated with the evolutionary gain of broader substrate specificities, which may have underpinned the expansion of the XTH family by enhancing duplicated gene survival. Together, this study highlights the Precambrian evolution of life and the mode of multigene family expansion in the evolutionary history of the XTH family.

Class-II dihydroorotate dehydrogenases from three phylogenetically distant fungi support anaerobic pyrimidine biosynthesis

Background

In most fungi, quinone-dependent Class-II dihydroorotate dehydrogenases (DHODs) are essential for pyrimidine biosynthesis. Coupling of these Class-II DHODHs to mitochondrial respiration makes their in vivo activity dependent on oxygen availability. Saccharomyces cerevisiae and closely related yeast species harbor a cytosolic Class-I DHOD (Ura1) that uses fumarate as electron acceptor and thereby enables anaerobic pyrimidine synthesis. Here, we investigate DHODs from three fungi (the Neocallimastigomycete Anaeromyces robustus and the yeasts Schizosaccharomyces japonicus and Dekkera bruxellensis) that can grow anaerobically but, based on genome analysis, only harbor a Class-II DHOD.

Results

Heterologous expression of putative Class-II DHOD-encoding genes from fungi capable of anaerobic, pyrimidine-prototrophic growth (Arura9, SjURA9, DbURA9) in an S. cerevisiae ura1Δ strain supported aerobic as well as anaerobic pyrimidine prototrophy. A strain expressing DbURA9 showed delayed anaerobic growth without pyrimidine supplementation. Adapted faster growing DbURA9-expressing strains showed mutations in FUM1, which encodes fumarase. GFP-tagged SjUra9 and DbUra9 were localized to S. cerevisiae mitochondria, while ArUra9, whose sequence lacked a mitochondrial targeting sequence, was localized to the yeast cytosol. Experiments with cell extracts showed that ArUra9 used free FAD and FMN as electron acceptors. Expression of SjURA9 in S. cerevisiae reproducibly led to loss of respiratory competence and mitochondrial DNA. A cysteine residue (C265 in SjUra9) in the active sites of all three anaerobically active Ura9 orthologs was shown to be essential for anaerobic activity of SjUra9 but not of ArUra9.

Conclusions

Activity of fungal Class-II DHODs was long thought to be dependent on an active respiratory chain, which in most fungi requires the presence of oxygen. By heterologous expression experiments in S. cerevisiae, this study shows that phylogenetically distant fungi independently evolved Class-II dihydroorotate dehydrogenases that enable anaerobic pyrimidine biosynthesis. Further structure–function studies are required to understand the mechanistic basis for the anaerobic activity of Class-II DHODs and an observed loss of respiratory competence in S. cerevisiae strains expressing an anaerobically active DHOD from Sch. japonicus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40694-021-00117-4.

Diversity of the Tryptophanase Gene and Its Evolutionary Implications in Living Organisms

Tryptophanase encoded by the gene tnaA is a pyridoxal phosphate-dependent enzyme that catalyses the conversion of tryptophan to indole, which is commonly used as an intra- and interspecies signalling molecule, particularly by microbes. However, the production of indole is rare in eukaryotic organisms. A nucleotide and protein database search revealed tnaA is commonly reported in various Gram-negative bacteria, but that only a few Gram-positive bacteria and archaea possess the gene. The presence of tnaA in eukaryotes, particularly protozoans and marine organisms, demonstrates the importance of this gene in the animal kingdom. Here, we document the distribution of tnaA and its acquisition and expansion among different taxonomic groups, many of which are usually categorized as non-indole producers. This study provides an opportunity to understand the intriguing role played by tnaA, and its distribution among various types of organisms.

Early-diverging fungal phyla: taxonomy, species concept, ecology, distribution, anthropogenic impact, and novel phylogenetic proposals

The increasing number of new fungal species described from all over the world along with the use of genetics to define taxa, has dramatically changed the classification system of early-diverging fungi over the past several decades. The number of phyla established for non-Dikarya fungi has increased from 2 to 17. However, to date, both the classification and phylogeny of the basal fungi are still unresolved. In this article, we review the recent taxonomy of the basal fungi and re-evaluate the relationships among early-diverging lineages of fungal phyla. We also provide information on the ecology and distribution in Mucoromycota and highlight the impact of chytrids on amphibian populations. Species concepts in Chytridiomycota, Aphelidiomycota, Rozellomycota, Neocallimastigomycota are discussed in this paper. To preserve the current application of the genus Nephridiophaga (Chytridiomycota: Nephridiophagales), a new type species, Nephridiophaga blattellae, is proposed.

A squalene-hopene cyclase in Schizosaccharomyces japonicus represents a eukaryotic adaptation to sterol-limited anaerobic environments

Biosynthesis of sterols, which are key constituents of canonical eukaryotic membranes, requires molecular oxygen. Anaerobic protists and deep-branching anaerobic fungi are the only eukaryotes in which a mechanism for sterol-independent growth has been elucidated. In these organisms, tetrahymanol, formed through oxygen-independent cyclization of squalene by a squalene-tetrahymanol cyclase, acts as a sterol surrogate. This study confirms an early report [C. J. E. A. Bulder, Antonie Van Leeuwenhoek, 37, 353-358 (1971)] that Schizosaccharomyces japonicus is exceptional among yeasts in growing anaerobically on synthetic media lacking sterols and unsaturated fatty acids. Mass spectrometry of lipid fractions of anaerobically grown Sch. japonicus showed the presence of hopanoids, a class of cyclic triterpenoids not previously detected in yeasts, including hop-22(29)-ene, hop-17(21)-ene, hop-21(22)-ene, and hopan-22-ol. A putative gene in Sch. japonicus showed high similarity to bacterial squalene-hopene cyclase (SHC) genes and in particular to those of Acetobacter species. No orthologs of the putative Sch. japonicus SHC were found in other yeast species. Expression of the Sch. japonicus SHC gene (Sjshc1) in Saccharomyces cerevisiae enabled hopanoid synthesis and stimulated anaerobic growth in sterol-free media, thus indicating that one or more of the hopanoids produced by SjShc1 could at least partially replace sterols. Use of hopanoids as sterol surrogates represents a previously unknown adaptation of eukaryotic cells to anaerobic growth. The fast anaerobic growth of Sch. japonicus in sterol-free media is an interesting trait for developing robust fungal cell factories for application in anaerobic industrial processes.

Genomic driven factors enhance biocatalyst-related cellulolysis potential in anaerobic digestion

Anaerobic digestion (AD) is a promising technology to recover bioenergy from biodegradable biomass, including cellulosic wastes. Through a few fractionation/separation techniques, cellulose has demonstrated its potential in AD, but the performance of the process is rather substrate-specific, as cellulolysis bacteria are sensitive to the enzyme-substrate interactions. Cellulosome is a self-assembled enzyme complex with many functionalized modules in the bacteria which has been gradually studied, however the genomic fingerprints of the culture-specific cellulosome in AD are relatively unclear especially under processing conditions. To clarify the key factors affecting the cellulosome induced cellulolysis, this review summarized the most recent publications of AD regarding the fates of cellulose, sources and functional genes of cellulosome, and omics methods for functional analyses. Different processes for organic treatment including applying food grinds in sewer, biomass valorization, cellulose fractionation, microaeration, and enzymatic hydrolysis enhanced fermentation, were highlighted to support the sustainable development of AD technology.

Identification of Oxygen-Independent Pathways for Pyridine Nucleotide and Coenzyme A Synthesis in Anaerobic Fungi by Expression of Candidate Genes in Yeast

Neocallimastigomycetes are unique examples of strictly anaerobic eukaryotes. This study investigates how these anaerobic fungi bypass reactions involved in synthesis of pyridine nucleotide cofactors and coenzyme A that, in canonical fungal pathways, require molecular oxygen. Analysis of Neocallimastigomycetes proteomes identified a candidate l-aspartate-decarboxylase (AdcA) and l-aspartate oxidase (NadB) and quinolinate synthase (NadA), constituting putative oxygen-independent bypasses for coenzyme A synthesis and pyridine nucleotide cofactor synthesis. The corresponding gene sequences indicated acquisition by ancient horizontal gene transfer (HGT) events involving bacterial donors. To test whether these enzymes suffice to bypass corresponding oxygen-requiring reactions, they were introduced into fms1Δ and bna2Δ Saccharomyces cerevisiae strains. Expression of nadA and nadB from Piromyces finnis and adcA from Neocallimastix californiae conferred cofactor prototrophy under aerobic and anaerobic conditions. This study simulates how HGT can drive eukaryotic adaptation to anaerobiosis and provides a basis for elimination of auxotrophic requirements in anaerobic industrial applications of yeasts and fungi. IMPORTANCE NAD (NAD⁺) and coenzyme A (CoA) are central metabolic cofactors whose canonical biosynthesis pathways in fungi require oxygen. Anaerobic gut fungi of the Neocallimastigomycota phylum are unique eukaryotic organisms that adapted to anoxic environments. Analysis of Neocallimastigomycota genomes revealed that these fungi might have developed oxygen-independent biosynthetic pathways for NAD⁺ and CoA biosynthesis, likely acquired through horizontal gene transfer (HGT) from prokaryotic donors. We confirmed functionality of these putative pathways under anaerobic conditions by heterologous expression in the yeast Saccharomyces cerevisiae. This approach, combined with sequence comparison, offers experimental insight on whether HGT events were required and/or sufficient for acquiring new traits. Moreover, our results demonstrate an engineering strategy for enabling S. cerevisiae to grow anaerobically in the absence of the precursor molecules pantothenate and nicotinate, thereby contributing to alleviate oxygen requirements and to move closer to prototrophic anaerobic growth of this industrially relevant yeast.

Distribution of methionine sulfoxide reductases in fungi and conservation of the free-methionine-R-sulfoxide reductase in multicellular eukaryotes

Methionine, either as a free amino acid or included in proteins, can be oxidized into methionine sulfoxide (MetO), which exists as R and S diastereomers. Almost all characterized organisms possess thiol-oxidoreductases named methionine sulfoxide reductase (Msr) enzymes to reduce MetO back to Met. MsrA and MsrB reduce the S and R diastereomers of MetO, respectively, with strict stereospecificity and are found in almost all organisms. Another type of thiol-oxidoreductase, the free-methionine-R-sulfoxide reductase (fRMsr), identified so far in prokaryotes and a few unicellular eukaryotes, reduces the R MetO diastereomer of the free amino acid. Moreover, some bacteria possess molybdenum-containing enzymes that reduce MetO, either in the free or protein-bound forms. All these Msrs play important roles in the protection of organisms against oxidative stress. Fungi are heterotrophic eukaryotes that colonize all niches on Earth and play fundamental functions, in organic matter recycling, as symbionts, or as pathogens of numerous organisms. However, our knowledge on fungal Msrs is still limited. Here, we performed a survey of msr genes in almost 700 genomes across the fungal kingdom. We show that most fungi possess one gene coding for each type of methionine sulfoxide reductase: MsrA, MsrB, and fRMsr. However, several fungi living in anaerobic environments or as obligate intracellular parasites were devoid of msr genes. Sequence inspection and phylogenetic analyses allowed us to identify non-canonical sequences with potentially novel enzymatic properties. Finaly, we identified several ocurences of msr horizontal gene transfer from bacteria to fungi.

Comprehensive and Functional Analysis of Horizontal Gene Transfer Events in Diatoms

Diatoms are a diverse group of mainly photosynthetic algae, responsible for 20% of worldwide oxygen production, which can rapidly respond to favorable conditions and often outcompete other phytoplankton. We investigated the contribution of horizontal gene transfer (HGT) to its ecological success. A large-scale phylogeny-based prokaryotic HGT detection procedure across nine sequenced diatoms showed that 3-5% of their proteome has a horizontal origin and a large influx occurred at the ancestor of diatoms. More than 90% of HGT genes are expressed, and species-specific HGT genes in Phaeodactylum tricornutum undergo strong purifying selection. Genes derived from HGT are implicated in several processes including environmental sensing and expand the metabolic toolbox. Cobalamin (vitamin B12) is an essential cofactor for roughly half of the diatoms and is only produced by bacteria. Five consecutive genes involved in the final synthesis of the cobalamin biosynthetic pathway, which could function as scavenging and repair genes, were detected as HGT. The full suite of these genes was detected in the cold-adapted diatom Fragilariopsis cylindrus. This might give diatoms originating from the Southern Ocean, a region typically depleted in cobalamin, a competitive advantage. Overall, we show that HGT is a prevalent mechanism that is actively used in diatoms to expand its adaptive capabilities.

A novel D-xylose isomerase from the gut of the wood feeding beetle Odontotaenius disjunctus efficiently expressed in Saccharomyces cerevisiae

Carbohydrate rich substrates such as lignocellulosic hydrolysates remain one of the primary sources of potentially renewable fuel and bulk chemicals. The pentose sugar d-xylose is often present in significant amounts along with hexoses. Saccharomyces cerevisiae can acquire the ability to metabolize d-xylose through expression of heterologous d-xylose isomerase (XI). This enzyme is notoriously difficult to express in S. cerevisiae and only fourteen XIs have been reported to be active so far. We cloned a new d-xylose isomerase derived from microorganisms in the gut of the wood-feeding beetle Odontotaenius disjunctus. Although somewhat homologous to the XI from Piromyces sp. E2, the new gene was identified as bacterial in origin and the host as a Parabacteroides sp. Expression of the new XI in S. cerevisiae resulted in faster aerobic growth than the XI from Piromyces on d-xylose media. The d-xylose isomerization rate conferred by the new XI was also 72% higher, while absolute xylitol production was identical in both strains. Interestingly, increasing concentrations of xylitol (up to 8 g L⁻¹) appeared not to inhibit d-xylose consumption. The newly described XI displayed 2.6 times higher specific activity, 37% lower K_M for d-xylose, and exhibited higher activity over a broader temperature range, retaining 51% of maximal activity at 30 °C compared with only 29% activity for the Piromyces XI.

Anaerobic Fungi: Past, Present, and Future

Anaerobic fungi (AF) play an essential role in feed conversion due to their potent fiber degrading enzymes and invasive growth. Much has been learned about this unusual fungal phylum since the paradigm shifting work of Colin Orpin in the 1970s, when he characterized the first AF. Molecular approaches targeting specific phylogenetic marker genes have facilitated taxonomic classification of AF, which had been previously been complicated by the complex life cycles and associated morphologies. Although we now have a much better understanding of their diversity, it is believed that there are still numerous genera of AF that remain to be described in gut ecosystems. Recent marker-gene based studies have shown that fungal diversity in the herbivore gut is much like the bacterial population, driven by host phylogeny, host genetics and diet. Since AF are major contributors to the degradation of plant material ingested by the host animal, it is understandable that there has been great interest in exploring the enzymatic repertoire of these microorganisms in order to establish a better understanding of how AF, and their enzymes, can be used to improve host health and performance, while simultaneously reducing the ecological footprint of the livestock industry. A detailed understanding of AF and their interaction with other gut microbes as well as the host animal is essential, especially when production of affordable high-quality protein and other animal-based products needs to meet the demands of an increasing human population. Such a mechanistic understanding, leading to more sustainable livestock practices, will be possible with recently developed -omics technologies that have already provided first insights into the different contributions of the fungal and bacterial population in the rumen during plant cell wall hydrolysis.

Fungal evolution: cellular, genomic and metabolic complexity

The question of how phenotypic and genomic complexity are inter-related and how they are shaped through evolution is a central question in biology that historically has been approached from the perspective of animals and plants. In recent years, however, fungi have emerged as a promising alternative system to address such questions. Key to their ecological success, fungi present a broad and diverse range of phenotypic traits. Fungal cells can adopt many different shapes, often within a single species, providing them with great adaptive potential. Fungal cellular organizations span from unicellular forms to complex, macroscopic multicellularity, with multiple transitions to higher or lower levels of cellular complexity occurring throughout the evolutionary history of fungi. Similarly, fungal genomes are very diverse in their architecture. Deep changes in genome organization can occur very quickly, and these phenomena are known to mediate rapid adaptations to environmental changes. Finally, the biochemical complexity of fungi is huge, particularly with regard to their secondary metabolites, chemical products that mediate many aspects of fungal biology, including ecological interactions. Herein, we explore how the interplay of these cellular, genomic and metabolic traits mediates the emergence of complex phenotypes, and how this complexity is shaped throughout the evolutionary history of Fungi.

Proteome specialization of anaerobic fungi during ruminal degradation of recalcitrant plant fiber

The rumen harbors a complex microbial mixture of archaea, bacteria, protozoa, and fungi that efficiently breakdown plant biomass and its complex dietary carbohydrates into soluble sugars that can be fermented and subsequently converted into metabolites and nutrients utilized by the host animal. While rumen bacterial populations have been well documented, only a fraction of the rumen eukarya are taxonomically and functionally characterized, despite the recognition that they contribute to the cellulolytic phenotype of the rumen microbiota. To investigate how anaerobic fungi actively engage in digestion of recalcitrant fiber that is resistant to degradation, we resolved genome-centric metaproteome and metatranscriptome datasets generated from switchgrass samples incubated for 48 h in nylon bags within the rumen of cannulated dairy cows. Across a gene catalog covering anaerobic rumen bacteria, fungi and viruses, a significant portion of the detected proteins originated from fungal populations. Intriguingly, the carbohydrate-active enzyme (CAZyme) profile suggested a domain-specific functional specialization, with bacterial populations primarily engaged in the degradation of hemicelluloses, whereas fungi were inferred to target recalcitrant cellulose structures via the detection of a number of endo- and exo-acting enzymes belonging to the glycoside hydrolase (GH) family 5, 6, 8, and 48. Notably, members of the GH48 family were amongst the highest abundant CAZymes and detected representatives from this family also included dockerin domains that are associated with fungal cellulosomes. A eukaryote-selected metatranscriptome further reinforced the contribution of uncultured fungi in the ruminal degradation of recalcitrant fibers. These findings elucidate the intricate networks of in situ recalcitrant fiber deconstruction, and importantly, suggest that the anaerobic rumen fungi contribute a specific set of CAZymes that complement the enzyme repertoire provided by the specialized plant cell wall degrading rumen bacteria.

Genomic and fossil windows into the secret lives of the most ancient fungi

Fungi have crucial roles in modern ecosystems as decomposers and pathogens, and they engage in various mutualistic associations with other organisms, especially plants. They have a lengthy geological history, and there is an emerging understanding of their impact on the evolution of Earth systems on a large scale. In this Review, we focus on the roles of fungi in the establishment and early evolution of land and freshwater ecosystems. Today, questions of evolution over deep time are informed by discoveries of new fossils and evolutionary analysis of new genomes. Inferences can be drawn from evolutionary analysis by comparing the genes and genomes of fungi with the biochemistry and development of their plant and algal hosts. We then contrast this emerging picture against evidence from the fossil record to develop a new, integrated perspective on the origin and early evolution of fungi.

Phylogenomic Analyses of Non-Dikarya Fungi Supports Horizontal Gene Transfer Driving Diversification of Secondary Metabolism in the Amphibian Gastrointestinal Symbiont, Basidiobolus

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.

Squalene-Tetrahymanol Cyclase Expression Enables Sterol-Independent Growth of Saccharomyces cerevisiae

The laboratory experiments described in this report simulate a proposed horizontal gene transfer event during the evolution of strictly anaerobic fungi. The demonstration that expression of a single heterologous gene sufficed to eliminate anaerobic sterol requirements in the model eukaryote Saccharomyces cerevisiae therefore contributes to our understanding of how sterol-independent eukaryotes evolved in anoxic environments. This report provides a proof of principle for a metabolic engineering strategy to eliminate sterol requirements in yeast strains that are applied in large-scale anaerobic industrial processes. The sterol-independent yeast strains described in this report provide a valuable platform for further studies on the physiological roles and impacts of sterols and sterol surrogates in eukaryotic cells.

Biosynthesis of sterols, which are considered essential components of virtually all eukaryotic membranes, requires molecular oxygen. Anaerobic growth of the yeast Saccharomyces cerevisiae therefore strictly depends on sterol supplementation of synthetic growth media. Neocallimastigomycota are a group of strictly anaerobic fungi which, instead of containing sterols, contain the pentacyclic triterpenoid “sterol surrogate” tetrahymanol, which is formed by cyclization of squalene. Here, we demonstrate that expression of the squalene-tetrahymanol cyclase gene TtTHC1 from the ciliate Tetrahymena thermophila enables synthesis of tetrahymanol by S. cerevisiae. Moreover, expression of TtTHC1 enabled exponential growth of anaerobic S. cerevisiae cultures in sterol-free synthetic media. After deletion of the ERG1 gene from a TtTHC1-expressing S. cerevisiae strain, native sterol synthesis was abolished and sustained sterol-free growth was demonstrated under anaerobic as well as aerobic conditions. Anaerobic cultures of TtTHC1-expressing S. cerevisiae on sterol-free medium showed lower specific growth rates and biomass yields than ergosterol-supplemented cultures, while their ethanol yield was higher. This study demonstrated that acquisition of a functional squalene-tetrahymanol cyclase gene offers an immediate growth advantage to S. cerevisiae under anaerobic, sterol-limited conditions and provides the basis for a metabolic engineering strategy to eliminate the oxygen requirements associated with sterol synthesis in yeasts.

IMPORTANCE The laboratory experiments described in this report simulate a proposed horizontal gene transfer event during the evolution of strictly anaerobic fungi. The demonstration that expression of a single heterologous gene sufficed to eliminate anaerobic sterol requirements in the model eukaryote Saccharomyces cerevisiae therefore contributes to our understanding of how sterol-independent eukaryotes evolved in anoxic environments. This report provides a proof of principle for a metabolic engineering strategy to eliminate sterol requirements in yeast strains that are applied in large-scale anaerobic industrial processes. The sterol-independent yeast strains described in this report provide a valuable platform for further studies on the physiological roles and impacts of sterols and sterol surrogates in eukaryotic cells.

Phylogenetic Analyses of Glycosyl Hydrolase Family 6 Genes in Tunicates: Possible Horizontal Transfer

Horizontal gene transfer (HGT) is the movement of genetic material between different species. Although HGT is less frequent in eukaryotes than in bacteria, several instances of HGT have apparently shaped animal evolution. One well-known example is the tunicate cellulose synthase gene, CesA, in which a gene, probably transferred from bacteria, greatly impacted tunicate evolution. A Glycosyl Hydrolase Family 6 (GH6) hydrolase-like domain exists at the C-terminus of tunicate CesA, but not in cellulose synthases of other organisms. The recent discovery of another GH6 hydrolase-like gene (GH6-1) in tunicate genomes further raises the question of how tunicates acquired GH6. To examine the probable origin of these genes, we analyzed the phylogenetic relationship of GH6 proteins in tunicates and other organisms. Our analyses show that tunicate GH6s, the GH6-1 gene, and the GH6 part of the CesA gene, form two independent, monophyletic gene groups. We also compared their sequence signatures and exon splice sites. All tunicate species examined have shared splice sites in GH6-containing genes, implying ancient intron acquisitions. It is likely that the tunicate CesA and GH6-1 genes existed in the common ancestor of all extant tunicates.

Assessing anaerobic gut fungal diversity in herbivores using D1/D2 large ribosomal subunit sequencing and multi-year isolation

The anaerobic gut fungi (AGF, Neocallimastigomycota) reside in the alimentary tracts of herbivores where they play a central role in the breakdown of plant material. Here, we report on the development of the hypervariable domains D1/D2 of the large ribosomal subunit (D1/D2 LSU) as a barcoding marker for the AGF. We generated a reference D1/D2 LSU database for all cultured AGF genera, as well as the majority of candidate genera encountered in prior internal transcribed spacer 1 (ITS1)-based surveys. Subsequently, a D1/D2 LSU-based diversity survey using long read PacBio SMRT sequencing was conducted on faecal samples from 21 wild and domesticated herbivores. Twenty-eight genera and candidate genera were identified, including multiple novel lineages that were predominantly, but not exclusively, identified in wild herbivores. Association between certain AGF genera and animal lifestyles, or animal host family was observed. Finally, to address the current paucity of AGF isolates, concurrent isolation efforts utilizing multiple approaches to maximize recovery yielded 216 isolates belonging to 12 different genera, several of which have no prior cultured-representatives. Our results establish the utility of D1/D2 LSU and PacBio sequencing for AGF diversity surveys, the culturability of multiple AGF taxa, and demonstrate that wild herbivores represent a yet-untapped reservoir of AGF diversity.

Toward a Fully Resolved Fungal Tree of Life

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

HGT in the human and skin commensal Malassezia: A bacterially derived flavohemoglobin is required for NO resistance and host interaction

The skin of humans and animals is colonized by commensal and pathogenic fungi and bacteria that share this ecological niche and have established microbial interactions. Malassezia are the most abundant fungal skin inhabitant of warm-blooded animals and have been implicated in skin diseases and systemic disorders, including Crohn's disease and pancreatic cancer. Flavohemoglobin is a key enzyme involved in microbial nitrosative stress resistance and nitric oxide degradation. Comparative genomics and phylogenetic analyses within the Malassezia genus revealed that flavohemoglobin-encoding genes were acquired through independent horizontal gene transfer events from different donor bacteria that are part of the mammalian microbiome. Through targeted gene deletion and functional complementation in Malassezia sympodialis, we demonstrated that bacterially derived flavohemoglobins are cytoplasmic proteins required for nitric oxide detoxification and nitrosative stress resistance under aerobic conditions. RNA-sequencing analysis revealed that endogenous accumulation of nitric oxide resulted in up-regulation of genes involved in stress response and down-regulation of the MalaS7 allergen-encoding genes. Solution of the high-resolution X-ray crystal structure of Malassezia flavohemoglobin revealed features conserved with both bacterial and fungal flavohemoglobins. In vivo pathogenesis is independent of Malassezia flavohemoglobin. Lastly, we identified an additional 30 genus- and species-specific horizontal gene transfer candidates that might have contributed to the evolution of this genus as the most common inhabitants of animal skin.

Genomic and proteomic biases inform metabolic engineering strategies for anaerobic fungi

Anaerobic fungi (Neocallimastigomycota) are emerging non-model hosts for biotechnology due to their wealth of biomass-degrading enzymes, yet tools to engineer these fungi have not yet been established. Here, we show that the anaerobic gut fungi have the most GC depleted genomes among 443 sequenced organisms in the fungal kingdom, which has ramifications for heterologous expression of genes as well as for emerging CRISPR-based genome engineering approaches. Comparative genomic analyses suggest that anaerobic fungi may contain cellular machinery to aid in sexual reproduction, yet a complete mating pathway was not identified. Predicted proteomes of the anaerobic fungi also contain an unusually large fraction of proteins with homopolymeric amino acid runs consisting of five or more identical consecutive amino acids. In particular, threonine runs are especially enriched in anaerobic fungal carbohydrate active enzymes (CAZymes) and this, together with a high abundance of predicted N-glycosylation motifs, suggests that gut fungal CAZymes are heavily glycosylated, which may impact heterologous production of these biotechnologically useful enzymes. Finally, we present a codon optimization strategy to aid in the development of genetic engineering tools tailored to these early-branching anaerobic fungi.

Rumen Protozoa Play a Significant Role in Fungal Predation and Plant Carbohydrate Breakdown

The rumen protozoa, alongside fungi, comprise the eukaryotic portion of the rumen microbiome. Rumen protozoa may account for up to 50% of biomass, yet their role in this ecosystem remains unclear. Early experiments inferred a role in carbohydrate and protein metabolism, but due to their close association with bacteria, definitively attributing these functions to the protozoa was challenging. The advent of ‘omic technologies has created opportunities to broaden our understanding of the rumen protozoa. This study aimed to utilize these methods to further our understanding of the role that protozoa play in the rumen in terms of their metabolic capacities, and in doing so, contribute valuable sequence data to reduce the chance of mis or under-representation of the rumen protozoa in meta’omic datasets. Rumen protozoa were isolated and purified using glucose-based sedimentation and differential centrifugation, extracted RNA was Poly(A) fraction enriched and DNase treated before use in a phage-based, cDNA metatranscriptomic library. Biochemical activity testing of the phage library showed 6 putatively positive plaques in response to carboxymethyl cellulose agar (indicative of cellulose activity), and no positive results for tributyrin (indicative of esterase/lipase activity) or egg yolk agar (indicative of proteolysis). Direct sequencing of the cDNA was also conducted using the Illumina HiSeq 2500. The metatranscriptome identified a wealth of carbohydrate-active enzymes which accounted for 8% of total reads. The most highly expressed carbohydrate-active enzymes were glycosyl hydrolases 5 and 11, polysaccharide lyases and deacetylases, xylanases and enzymes active against pectin, mannan and chitin; the latter likely used to digest rumen fungi which contain a chitin-rich cell membrane. Codon usage analysis of expressed genes also showed evidence of horizontal gene transfer, suggesting that many of these enzymes were acquired from the rumen bacteria in an evolutionary response to the carbohydrate-rich environment of the rumen. This study provides evidence of the significant contribution that the protozoa make to carbohydrate breakdown in the rumen, potentially using horizontally acquired genes, and highlights their predatory capacity.

Horizontal transfer of a natterin-like toxin encoding gene within the holobiont of the reef building coral Acropora digitifera (Cnidaria: Anthozoa: Scleractinia) and across multiple animal linages

Phylogenetic evidence is provided for horizontal transfer of a natterin-like toxin encoding gene from fungi into the genome of the coral Acropora digitifera. Sequencing analysis of the coral tissues supported that a fungal taxon predicted to be the most likely gene donor was represented in the coral microbiome. Further bioinformatics data suggested widespread recruitment of the natterin-like gene into venomous terrestrial invertebrates, and repositioning of this gene to non-toxic functions in non-venomous teleost fish.