Co‑cultivation of the anaerobic fungus Caecomyces churrovis with Methanobacterium bryantii enhances transcription of carbohydrate binding modules, dockerins, and pyruvate formate lyases on specific substrates

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.

Structural and functional analysis of the active cow rumen's microbial community provides a catalogue of genes and microbes participating in the deconstruction of cardoon biomass

BACKGROUND

Ruminal microbial communities enriched on lignocellulosic biomass have shown considerable promise for the discovery of microorganisms and enzymes involved in digesting cell wall compounds, a key bottleneck in the development of second-generation biofuels and bioproducts, enabling a circular bioeconomy. Cardoon (Cynara cardunculus) is a promising inedible energy crop for current and future cellulosic biorefineries and the emerging bioenergy and bioproducts industries. The rumen microbiome can be considered an anaerobic "bioreactor", where the resident microbiota carry out the depolymerization and hydrolysis of plant cell wall polysaccharides (PCWPs) through the catalytic action of fibrolytic enzymes. In this context, the rumen microbiota represents a potential source of microbes and fibrolytic enzymes suitable for biofuel production from feedstocks. In this study, metatranscriptomic and 16S rRNA sequencing were used to profile the microbiome and to investigate the genetic features within the microbial community adherent to the fiber fractions of the rumen content and to the residue of cardoon biomass incubated in the rumen of cannulated cows.

RESULTS

The metatranscriptome of the cardoon and rumen fibre-adherent microbial communities were dissected in their functional and taxonomic components. From a functional point of view, transcripts involved in the methanogenesis from CO2 and H2, and from methanol were over-represented in the cardoon-adherent microbial community and were affiliated with the Methanobrevibacter and Methanosphaera of the Euryarchaeota phylum. Transcripts encoding glycoside hydrolases (GHs), carbohydrate-binding modules (CBMs), carbohydrate esterases (CEs), polysaccharide lyases (PLs), and glycoside transferases (GTs) accounted for 1.5% (6,957) of the total RNA coding transcripts and were taxonomically affiliated to major rumen fibrolytic microbes, such as Oscillospiraceae, Fibrobacteraceae, Neocallimastigaceae, Prevotellaceae, Lachnospiraceae, and Treponemataceae. The comparison of the expression profile between cardoon and rumen fiber-adherent microbial communities highlighted that specific fibrolytic enzymes were potentially responsible for the breakdown of cardoon PCWPs, which was driven by specific taxa, mainly Ruminococcus, Treponema, and Neocallimastigaceae.

CONCLUSIONS

Analysis of 16S rRNA and metatranscriptomic sequencing data revealed that the cow rumen microbiome harbors a repertoire of new enzymes capable of degrading PCWPs. Our results demonstrate the feasibility of using metatranscriptomics of enriched microbial RNA as a potential approach for accelerating the discovery of novel cellulolytic enzymes that could be harnessed for biotechnology. This research contributes a relevant perspective towards degrading cellulosic biomass and providing an economical route to the production of advanced biofuels and high-value bioproducts.

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.

Targeted rRNA depletion enables efficient mRNA sequencing in diverse bacterial species and complex co-cultures

IMPORTANCE

Microbes present one of the most diverse sources of biochemistry in nature, and mRNA sequencing provides a comprehensive view of this biological activity by quantitatively measuring microbial transcriptomes. However, efficient mRNA capture for sequencing presents significant challenges in prokaryotes as mRNAs are not poly-adenylated and typically make up less than 5% of total RNA compared with rRNAs that exceed 80%. Recently developed methods for sequencing bacterial mRNA typically rely on depleting rRNA by tiling large probe sets against rRNAs; however, such approaches are expensive, time-consuming, and challenging to scale to varied bacterial species and complex microbial communities. Therefore, we developed EMBR-seq+, a method that requires fewer than 10 short oligonucleotides per rRNA to achieve up to 99% rRNA depletion in diverse bacterial species. Finally, EMBR-seq+ resulted in a deeper view of the transcriptome, enabling systematic quantification of how microbial interactions result in altering the transcriptional state of bacteria within co-cultures.

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.

High-quality RNA extraction and the regulation of genes encoding cellulosomes are correlated with growth stage in anaerobic fungi

Anaerobic fungi produce biomass-degrading enzymes and natural products that are important to harness for several biotechnology applications. Although progress has been made in the development of methods for extracting nucleic acids for genomic and transcriptomic sequencing of these fungi, most studies are limited in that they do not sample multiple fungal growth phases in batch culture. In this study, we establish a method to harvest RNA from fungal monocultures and fungal-methanogen co-cultures, and also determine an optimal time frame for high-quality RNA extraction from anaerobic fungi. Based on RNA quality and quantity targets, the optimal time frame in which to harvest anaerobic fungal monocultures and fungal-methanogen co-cultures for RNA extraction was 2-5 days of growth post-inoculation. When grown on cellulose, the fungal strain Anaeromyces robustus cocultivated with the methanogen Methanobacterium bryantii upregulated genes encoding fungal carbohydrate-active enzymes and other cellulosome components relative to fungal monocultures during this time frame, but expression patterns changed at 24-hour intervals throughout the fungal growth phase. These results demonstrate the importance of establishing methods to extract high-quality RNA from anaerobic fungi at multiple time points during batch cultivation.

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.

Expression and characterization of spore coat CotH kinases from the cellulosomes of anaerobic fungi (Neocallimastigomycetes)

Anaerobic fungi (Neocallimastigomycetes) found in the guts of herbivores are biomass deconstruction specialists with a remarkable ability to extract sugars from recalcitrant plant material. Anaerobic fungi, as well as many species of anaerobic bacteria, deploy multi-enzyme complexes called cellulosomes, which modularly tether together hydrolytic enzymes, to accelerate biomass hydrolysis. While the majority of genomically encoded cellulosomal genes in Neocallimastigomycetes are biomass degrading enzymes, the second largest family of cellulosomal genes encode spore coat CotH domains, whose contribution to fungal cellulosome and/or cellular function is unknown. Structural bioinformatics of CotH proteins from the anaerobic fungus Piromyces finnis shows anaerobic fungal CotH domains conserve key ATP and Mg²⁺ binding motifs from bacterial Bacillus CotH proteins known to act as protein kinases. Experimental characterization further demonstrates ATP hydrolysis activity in the presence and absence of substrate from two cellulosomal P. finnis CotH proteins when recombinantly produced in E. coli. These results present foundational evidence for CotH activity in anaerobic fungi and provide a path towards elucidating the functional contribution of this protein family to fungal cellulosome assembly and activity.

Metabolic influence of core ciliates within the rumen microbiome

Protozoa comprise a major fraction of the microbial biomass in the rumen microbiome, of which the entodiniomorphs (order: Entodiniomorphida) and holotrichs (order: Vestibuliferida) are consistently observed to be dominant across a diverse genetic and geographical range of ruminant hosts. Despite the apparent core role that protozoal species exert, their major biological and metabolic contributions to rumen function remain largely undescribed in vivo. Here, we have leveraged (meta)genome-centric metaproteomes from rumen fluid samples originating from both cattle and goats fed diets with varying inclusion levels of lipids and starch, to detail the specific metabolic niches that protozoa occupy in the context of their microbial co-habitants. Initial proteome estimations via total protein counts and label-free quantification highlight that entodiniomorph species Entodinium and Epidinium as well as the holotrichs Dasytricha and Isotricha comprise an extensive fraction of the total rumen metaproteome. Proteomic detection of protozoal metabolism such as hydrogenases (Dasytricha, Isotricha, Epidinium, Enoploplastron), carbohydrate-active enzymes (Epidinium, Diplodinium, Enoploplastron, Polyplastron), microbial predation (Entodinium) and volatile fatty acid production (Entodinium and Epidinium) was observed at increased levels in high methane-emitting animals. Despite certain protozoal species having well-established reputations for digesting starch, they were unexpectedly less detectable in low methane emitting-animals fed high starch diets, which were instead dominated by propionate/succinate-producing bacterial populations suspected of being resistant to predation irrespective of host. Finally, we reaffirmed our abovementioned observations in geographically independent datasets, thus illuminating the substantial metabolic influence that under-explored eukaryotic populations have in the rumen, with greater implications for both digestion and methane metabolism.

Lignin deconstruction by anaerobic fungi

Lignocellulose forms plant cell walls, and its three constituent polymers, cellulose, hemicellulose and lignin, represent the largest renewable organic carbon pool in the terrestrial biosphere. Insights into biological lignocellulose deconstruction inform understandings of global carbon sequestration dynamics and provide inspiration for biotechnologies seeking to address the current climate crisis by producing renewable chemicals from plant biomass. Organisms in diverse environments disassemble lignocellulose, and carbohydrate degradation processes are well defined, but biological lignin deconstruction is described only in aerobic systems. It is currently unclear whether anaerobic lignin deconstruction is impossible because of biochemical constraints or, alternatively, has not yet been measured. We applied whole cell-wall nuclear magnetic resonance, gel-permeation chromatography and transcriptome sequencing to interrogate the apparent paradox that anaerobic fungi (Neocallimastigomycetes), well-documented lignocellulose degradation specialists, are unable to modify lignin. We find that Neocallimastigomycetes anaerobically break chemical bonds in grass and hardwood lignins, and we further associate upregulated gene products with the observed lignocellulose deconstruction. These findings alter perceptions of lignin deconstruction by anaerobes and provide opportunities to advance decarbonization biotechnologies that depend on depolymerizing lignocellulose.

A Genetic Engineering Toolbox for the Lignocellulolytic Anaerobic Gut Fungus Neocallimastix frontalis

Anaerobic fungi are powerful platforms for biotechnology that remain unexploited due to a lack of genetic tools. These gut fungi encode the largest number of lignocellulolytic carbohydrate active enzymes (CAZymes) in the fungal kingdom, making them attractive for applications in renewable energy and sustainability. However, efforts to genetically modify anaerobic fungi have remained limited due to inefficient methods for DNA uptake and a lack of characterized genetic parts. We demonstrate that anaerobic fungi are naturally competent for DNA and leverage this to develop a nascent genetic toolbox informed by recently acquired genomes for transient transformation of anaerobic fungi. We validate multiple selectable markers (HygR and Neo), an anaerobic reporter protein (iRFP702), enolase and TEF1A promoters, TEF1A terminator, and a nuclear localization tag for protein compartmentalization. This work establishes novel methods to reliably transform the anaerobic fungus Neocallimastix frontalis, thereby paving the way for strain development and various synthetic biology applications.

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.

Understanding the microbial fibre degrading communities & processes in the equine gut

The equine gastrointestinal tract is a self-sufficient fermentation system, housing a complex microbial consortium that acts synergistically and independently to break down complex lignocellulolytic material that enters the equine gut. Despite being strict herbivores, equids such as horses and zebras lack the diversity of enzymes needed to completely break down plant tissue, instead relying on their resident microbes to carry out fibrolysis to yield vital energy sources such as short chain fatty acids. The bulk of equine digestion occurs in the large intestine, where digesta is fermented for 36-48 h through the synergistic activities of bacteria, fungi, and methanogenic archaea. Anaerobic gut dwelling bacteria and fungi break down complex plant polysaccharides through combined mechanical and enzymatic strategies, and notably possess some of the greatest diversity and repertoire of carbohydrate active enzymes among characterized microbes. In addition to the production of enzymes, some equid-isolated anaerobic fungi and bacteria have been shown to possess cellulosomes, powerful multi-enzyme complexes that further enhance break down. The activities of both anaerobic fungi and bacteria are further facilitated by facultatively aerobic yeasts and methanogenic archaea, who maintain an optimal environment for fibrolytic organisms, ultimately leading to increased fibrolytic microbial counts and heightened enzymatic activity. The unique interactions within the equine gut as well as the novel species and powerful mechanisms employed by these microbes makes the equine gut a valuable ecosystem to study fibrolytic functions within complex communities. This review outlines the primary taxa involved in fibre break down within the equine gut and further illuminates the enzymatic strategies and metabolic pathways used by these microbes. We discuss current methods used in analysing fibrolytic functions in complex microbial communities and propose a shift towards the development of functional assays to deepen our understanding of this unique ecosystem.

Co‑cultivation of anaerobic fungi with Clostridium acetobutylicum bolsters butyrate and butanol production from cellulose and lignocellulose

A system for co-cultivation of anaerobic fungi with anaerobic bacteria was established based on lactate cross-feeding to produce butyrate and butanol from plant biomass. Several co-culture formulations were assembled that consisted of anaerobic fungi (Anaeromyces robustus, Neocallimastix californiae, or Caecomyces churrovis) with the bacterium Clostridium acetobutylicum. Co-cultures were grown simultaneously (e.g., 'one pot'), and compared to cultures where bacteria were cultured in fungal hydrolysate sequentially. Fungal hydrolysis of lignocellulose resulted in 7-11 mM amounts of glucose and xylose, as well as acetate, formate, ethanol, and lactate to support clostridial growth. Under these conditions, one-stage simultaneous co-culture of anaerobic fungi with C. acetobutylicum promoted the production of butyrate up to 30 mM. Alternatively, two-stage growth slightly promoted solventogenesis and elevated butanol levels (∼4-9 mM). Transcriptional regulation in the two-stage growth condition indicated that this cultivation method may decrease the time required to reach solventogenesis and induce the expression of cellulose-degrading genes in C. acetobutylicum due to relieved carbon-catabolite repression. Overall, this study demonstrates a proof of concept for biobutanol and bio-butyrate production from lignocellulose using an anaerobic fungal-bacterial co-culture system.

Carbohydrate active enzyme system in rumen fungi: a review

Hydrolysis and dehydration reactions of carbohydrates, which are used as energy raw materials by all living things in nature, are controlled by Carbohydrate Active Enzyme (CAZy) systems. These enzymes are also used in different industrial areas today. There are different types of microorganisms that have the CAZy system and are used in the industrial sector. Apart from current organisms, there are also rumen fungi within the group of candidate microorganisms with the CAZy system. It has been reported that xylanase (EC3.2.1.8 and EC3.2.1.37) enzyme, a member of the glycoside hydrolase enzyme family obtained from Trichoderma sp. and used especially in areas such as bread, paper, and feed industry, is more synthesized in rumen fungi such as Orpinomyces sp. and Neocallimastix sp. Therefore, this study reviews Neocallimastixsp., Orpinomyces sp., Caecomyces sp., Piromyces sp., and Anaeromyces sp., registered in the CAZy and Mycocosm database for rumen fungi to have both CAZy enzyme activity and to be an alternative microorganism in the industry. Furthermore the CAZy enzyme activities of the strains are investigated. The review shows thatNeocallimax sp. and Orpinomyces sp. areconsidered as candidate microorganisms.

Biofilm disruption enhances growth rate and carbohydrate-active enzyme production in anaerobic fungi

Anaerobic gut fungi (AGF) are lignocellulose degraders that naturally form biofilms in the rumen of large herbivores and in standard culture techniques. While biofilm formation enhances biomass degradation and carbohydrate-active enzyme (CAZyme) production in some bacteria and aerobic fungi, gene expression and metabolism in AGF biofilms have not been compared to non-biofilm cultures. Here, using the tunable morphology of the non-rhizoidal AGF, Caecomyces churrovis, the impacts of biofilm formation on AGF gene expression, metabolic flux, growth rate, and xylan degradation rate are quantified to inform future industrial scale-up efforts. Contrary to previous findings, C. churrovis upregulated catabolic CAZymes in stirred culture relative to biofilm culture. Using a de novo transcriptome, 197 new transcripts with predicted CAZyme function were identified. Stirred cultures grew and degraded xylan significantly faster than biofilm-forming cultures with negligible differences in primary metabolic flux, offering a way to accelerate AGF biomass valorization without altering the fermentation product profile. The rhizoidal AGF, Neocallimastix lanati, also grew faster with stirring on a solid plant substrate, suggesting that the advantages of stirred C. churrovis cultures may apply broadly to other AGF.