Assembly of 913 microbial genomes from metagenomic sequencing of the cow rumen

Single-cell transcriptomics across 2,534 microbial species reveals functional heterogeneity in the rumen microbiome

Deciphering the activity of individual microbes within complex communities and environments remains a challenge. Here we describe the development of microbiome single-cell transcriptomics using droplet-based single-cell RNA sequencing and pangenome-based computational analysis to characterize the functional heterogeneity of the rumen microbiome. We generated a microbial genome database (the Bovine Gastro Microbial Genome Map) as a functional reference map for the construction of a single-cell transcriptomic atlas of the rumen microbiome. The atlas includes 174,531 microbial cells and 2,534 species, of which 172 are core active species grouped into 12 functional clusters. We detected single-cell-level functional roles, including a key role for Basfia succiniciproducens in the carbohydrate metabolic niche of the rumen microbiome. Furthermore, we explored functional heterogeneity and reveal metabolic niche trajectories driven by biofilm formation pathway genes within B. succiniciproducens. Our results provide a resource for studying the rumen microbiome and illustrate the diverse functions of individual microbial cells that drive their ecological niche stability or adaptation within the ecosystem.

Genome reduction and horizontal gene transfer in the evolution of Endomicrobia-rise and fall of an intracellular symbiosis with termite gut flagellates

Bacterial endosymbionts of eukaryotic hosts typically experience massive genome reduction, but the underlying evolutionary processes are often obscured by the lack of free-living relatives. Endomicrobia, a family-level lineage of host-associated bacteria in the phylum Elusimicrobiota that comprises both free-living representatives and endosymbionts of termite gut flagellates, are an excellent model to study evolution of intracellular symbionts. We reconstructed 67 metagenome-assembled genomes (MAGs) of Endomicrobiaceae among more than 1,700 MAGs from the gut microbiota of a wide range of termites. Phylogenomic analysis confirmed a sister position of representatives from termites and ruminants, and allowed to propose eight new genera in the radiation of Endomicrobiaceae. Comparative genome analysis documented progressive genome erosion in the new genus Endomicrobiellum, which comprises all flagellate endosymbionts characterized to date. Massive gene losses were accompanied by the acquisition of new functions by horizontal gene transfer, which led to a shift from a glucose-based energy metabolism to one based on sugar phosphates. The breakdown of glycolysis and many anabolic pathways for amino acids and cofactors in several subgroups was compensated by the independent acquisition of new uptake systems, including an ATP/ADP antiporter, from other gut microbiota. The putative donors are mostly flagellate endosymbionts from other bacterial phyla, including several, hitherto unknown lineages of uncultured Alphaproteobacteria, documenting the importance of horizontal gene transfer in the convergent evolution of these intracellular symbioses. The loss of almost all biosynthetic capacities in some lineages of Endomicrobiellum suggests that their originally mutualistic relationship with flagellates is on its decline.IMPORTANCEUnicellular eukaryotes are frequently colonized by bacterial and archaeal symbionts. A prominent example are the cellulolytic gut flagellates of termites, which harbor diverse but host-specific bacterial symbionts that occur exclusively in termite guts. One of these lineages, the so-called Endomicrobia, comprises both free-living and endosymbiotic representatives, which offers the unique opportunity to study the evolutionary processes underpinning the transition from a free-living to an intracellular lifestyle. Our results revealed a progressive gene loss in energy metabolism and biosynthetic pathways, compensated by the acquisition of new functions via horizontal gene transfer from other gut bacteria, and suggest the eventual breakdown of an initially mutualistic symbiosis. Evidence for convergent evolution of unrelated endosymbionts reflects adaptations to the intracellular environment of termite gut flagellates.

Exploration of mobile genetic elements in the ruminal microbiome of Nellore cattle

Metagenomics has made it feasible to elucidate the intricacies of the ruminal microbiome and its role in the differentiation of animal production phenotypes of significance. The search for mobile genetic elements (MGEs) has taken on great importance, as they play a critical role in the transfer of genetic material between organisms. Furthermore, these elements serve a dual purpose by controlling populations through lytic bacteriophages, thereby maintaining ecological equilibrium and driving the evolutionary progress of host microorganisms. In this study, we aimed to identify the association between ruminal bacteria and their MGEs in Nellore cattle using physical chromosomal links through the Hi-C method. Shotgun metagenomic sequencing and the proximity ligation method ProxiMeta were used to analyze DNA, getting 1,713,111,307 bp, which gave rise to 107 metagenome-assembled genomes from rumen samples of four Nellore cows maintained on pasture. Taxonomic analysis revealed that most of the bacterial genomes belonged to the families Lachnospiraceae, Bacteroidaceae, Ruminococcaceae, Saccharofermentanaceae, and Treponemataceae and mostly encoded pathways for central carbon and other carbohydrate metabolisms. A total of 31 associations between host bacteria and MGE were identified, including 17 links to viruses and 14 links to plasmids. Additionally, we found 12 antibiotic resistance genes. To our knowledge, this is the first study in Brazilian cattle that connect MGEs with their microbial hosts. It identifies MGEs present in the rumen of pasture-raised Nellore cattle, offering insights that could advance biotechnology for food digestion and improve ruminant performance in production systems.

Multi-omics strategy reveals potential role of antimicrobial resistance and virulence factor genes responsible for Simmental diarrheic calves caused by Escherichia coli

Escherichia coli (E. coli) is reported to be an important pathogen associated with calf diarrhea. Antibiotic resistance genes (ARGs) and virulence factor genes (VFGs) pose a considerable threat to both animal and human health. However, little is known about the characterization of ARGs and VFGs presented in the gut microbiota of diarrheic calves caused by E. coli. In this study, we used multi-omics strategy to analyze the ARG and VFG profiles of Simmental calves with diarrhea caused by E. coli K99. We found that gut bacterial composition and their microbiome metabolic functions varied greatly in diarrheic calves compared to healthy calves. In total, 175 ARGs were identified, and diarrheal calves showed a significantly higher diversity and abundance of ARGs than healthy calves. Simmental calves with diarrhea showed higher association of VFGs with pili function, curli assembly, and ferrienterobactin transport of E. coli. Co-occurrence patterns based on Pearson correlation analysis revealed that E. coli had a highly significant (P < 0.0001) correlation coefficient (>0.8) with 16 ARGs and 7 VFGs. Metabolomics analysis showed that differentially expressed metabolites in Simmental calves with diarrhea displayed a high correlation with the aforementioned ARGs and VFGs. Phylotype analysis of E. coli genomes showed that the predominant phylogroup B1 in diarrheic Simmental calves was associated with 10 ARGs and 3 VFGs. These findings provide an overview of the diversity and abundance of the gut microbiota in diarrheic calves caused by E. coli and pave the way for further studies on the mechanisms of antibiotic resistance and virulence in the calves affected with diarrhea.IMPORTANCESimmental is a well-recognized beef cattle breed worldwide. They also suffer significant economic losses due to diarrhea. In this study, fecal metagenomic analysis was applied to characterize the antibiotic resistance gene (ARG) and virulence factor gene (VFG) profiles of diarrheic Simmental calves. We identified key ARGs and VFGs correlated with Escherichia coli isolated from Simmental calves. Additionally, metabolomics analysis showed that differentially expressed metabolites in Simmental calves with diarrhea displayed a high correlation with the aforementioned ARGs and VFGs. Our findings provide an insight into the diversity and abundance of the gut microbiota in diarrheic calves caused by Escherichia coli and pave the way for further studies on the mechanisms of antibiotic resistance and virulence in the diarrheal calves from cattle hosts.

A Review of the Rumen Microbiota and the Different Molecular Techniques Used to Identify Microorganisms Found in the Rumen Fluid of Ruminants

Variations in environments, including climate, diet, and agricultural practices, significantly impact the composition and microbial activity. A profound understanding of these adaptations allows for the improvement of nutrition and ruminant production. Therefore, this review aims to compile data from the literature on the rumen microbiota and molecular techniques for identifying the different types of microorganisms from the rumen fluid of ruminants. Analyzing the literature on rumen microbiology in different ruminants is complex due to microbial interactions, influenced by the environment and nutrition of these animals. In addition, it is worth noting that the genera of protozoa and fungi most evident in the studies used in this review on the microbiology of rumen fluid were Entodinium spp. and Aspergillus spp., respectively, and Fibrobacter spp. for bacteria. About the techniques used, it can be seen that DNA extraction, amplification, and sequencing were the most cited in the studies evaluated. Therefore, this review describes what is present in the literature and provides an overview of the main microbial agents in the rumen and the molecular techniques used.

Metagenomic and metabolomic analyses reveal differences in rumen microbiota between grass- and grain-fed Sanhe heifers

Introduction

The aim of this study was to investigate the effects of diets on the composition and function of rumen microbiome and metabolites in Sanhe heifers.

Methods

Metagenomic and metabolomic analyses were performed using rumen fluid samples collected from Sanhe heifers (n = 20) with similar body weights and ages from grass-fed and grain-fed systems.

Results

The grain-fed group exhibited more intensive rumen fermentation than the grass-fed group. However, the grass-fed group exhibited carbohydrate metabolism and methane production higher than that of the grain-fed group; these increases were observed as a higher abundance of various bacterial phyla (Firmicutes, Bacteroidetes, Actinobacteria, Lentisphaerae, and Verrucomicrobia), families (Lachnospiraceae, Eubacteriaceae, and Eggerthellaceae), and the archaeal family Methanobacteriaceae. A comparison of genes encoding carbohydrate-active enzymes, using Kyoto Encyclopedia of Genes and Genome profiles, revealed noteworthy differences in the functions of rumen microbiota; these differences were largely dependent on the feeding system.

Conclusion

These results could help manipulate and regulate feed efficiency in Sanhe cattle.

Viruses contribute to microbial diversification in the rumen ecosystem and are associated with certain animal production traits

BACKGROUND

The rumen microbiome enables ruminants to digest otherwise indigestible feedstuffs, thereby facilitating the production of high-quality protein, albeit with suboptimal efficiency and producing methane. Despite extensive research delineating associations between the rumen microbiome and ruminant production traits, the functional roles of the pervasive and diverse rumen virome remain to be determined.

RESULTS

Leveraging a recent comprehensive rumen virome database, this study analyzes virus-microbe linkages, at both species and strain levels, across 551 rumen metagenomes, elucidating patterns of microbial and viral diversity, co-occurrence, and virus-microbe interactions. Additionally, this study assesses the potential role of rumen viruses in microbial diversification by analyzing prophages found in rumen metagenome-assembled genomes. Employing CRISPR-Cas spacer-based matching and virus-microbe co-occurrence network analysis, this study suggests that the viruses in the rumen may regulate microbes at strain and community levels through both antagonistic and mutualistic interactions. Moreover, this study establishes that the rumen virome demonstrates responsiveness to dietary shifts and associations with key animal production traits, including feed efficiency, lactation performance, weight gain, and methane emissions.

CONCLUSIONS

These findings provide a substantive framework for further investigations to unravel the functional roles of the virome in the rumen in shaping the microbiome and influencing overall animal production performance. Video Abstract.

Molecular detection of ruminal micro-flora and micro-fauna in Saudi Arabian camels: Effects of season and region

This study investigated and explored the availability of micro-flora and micro-fauna in the ruminal contents of Arabian camel (Camelus dromedarius) from three different regions in Saudi Arabia along with two seasons. Samples were prepared and tested by conventional polymerase chain reaction (PCR). This study confirmed that the bacterial flora were dominating over other microbes. Different results of the availability of each microbe in each region and season were statistically analyzed and discussed. There was no significant effect of season on the micro-flora or micro-fauna however, the location revealed a positive effect with Ruminococcus flavefaciens (p < 0 0.03) in the eastern region. This study was the first to investigate the abundance of micro-flora and micro-fauna in the ruminal contents of camels of Saudi Arabia. This study underscores the significance of camel ruminal micro-flora and micro-fauna abundance, highlighting their correlation with both seasonality and geographic location. This exploration enhances our comprehension of camel rumination and digestion processes. The initial identification of these microbial communities serves as a foundational step, laying the groundwork for future in-depth investigations into camel digestibility and nutritional requirements.

Reference gene catalog and metagenome-assembled genomes from the gut microbiome reveal the microbial composition, antibiotic resistome, and adaptability of a lignocellulose diet in the giant panda

The giant panda, a strict herbivore that feeds on bamboo, still retains a typical carnivorous digestive system. Reference catalogs of microbial genes and genomes are lacking, largely limiting the antibiotic resistome and functional exploration of the giant panda gut microbiome. Here, we integrated 177 fecal metagenomes of captive and wild giant pandas to construct a giant panda integrated gene catalog (GPIGC) comprised of approximately 4.5 million non-redundant genes and reconstruct 393 metagenome-assembled genomes (MAGs). Taxonomic and functional characterization of genes revealed that the captivity of the giant panda significantly changed the core microbial composition and the distribution of microbial genes. Higher abundance and prevalence of antibiotic resistance genes (ARGs) were detected in the guts of captive giant pandas, and ARG distribution was influenced by geography, for both captive and wild individuals. Escherichia, as the prevalent genus in the guts of captive giant pandas, was the main carrier of ARGs, meaning there is a high risk of ARG transmission by Escherichia. We also found that multiple mcr gene variants, conferring plasmid-mediated mobile colistin resistance, were widespread in the guts of captive and wild giant pandas. There were low proportions of carbohydrate-active enzyme (CAZyme) genes in GPIGC and MAGs compared with several omnivorous and herbivorous mammals. Many members of Clostridium MAGs were significantly enriched in the guts of adult, old and wild giant pandas. The genomes of isolates and MAGs of Clostridiaceae harbored key genes or enzymes in complete pathways for degrading lignocellulose and producing short-chain fatty acids (SCFAs), indicating the potential of these bacteria to utilize the low-nutrient bamboo diet. Overall, our data presented an exhaustive reference gene catalog and MAGs in giant panda gut and provided a comprehensive understanding of the antibiotic resistome and microbial adaptability for a high-lignocellulose diet.

Ten common issues with reference sequence databases and how to mitigate them

Metagenomic sequencing has revolutionized our understanding of microbiology. While metagenomic tools and approaches have been extensively evaluated and benchmarked, far less attention has been given to the reference sequence database used in metagenomic classification. Issues with reference sequence databases are pervasive. Database contamination is the most recognized issue in the literature; however, it remains relatively unmitigated in most analyses. Other common issues with reference sequence databases include taxonomic errors, inappropriate inclusion and exclusion criteria, and sequence content errors. This review covers ten common issues with reference sequence databases and the potential downstream consequences of these issues. Mitigation measures are discussed for each issue, including bioinformatic tools and database curation strategies. Together, these strategies present a path towards more accurate, reproducible and translatable metagenomic sequencing.

Including microbiome information in a multi-trait genomic evaluation: a case study on longitudinal growth performance in beef cattle

BACKGROUND

Growth rate is an important component of feed conversion efficiency in cattle and varies across the different stages of the finishing period. The metabolic effect of the rumen microbiome is essential for cattle growth, and investigating the genomic and microbial factors that underlie this temporal variation can help maximize feed conversion efficiency at each growth stage.

RESULTS

By analysing longitudinal body weights during the finishing period and genomic and metagenomic data from 359 beef cattle, our study demonstrates that the influence of the host genome on the functional rumen microbiome contributes to the temporal variation in average daily gain (ADG) in different months (ADG1, ADG2, ADG3, ADG4). Five hundred and thirty-three additive log-ratio transformed microbial genes (alr-MG) had non-zero genomic correlations (rg) with at least one ADG-trait (ranging from |0.21| to |0.42|). Only a few alr-MG correlated with more than one ADG-trait, which suggests that a differential host-microbiome determinism underlies ADG at different stages. These alr-MG were involved in ribosomal biosynthesis, energy processes, sulphur and aminoacid metabolism and transport, or lipopolysaccharide signalling, among others. We selected two alternative subsets of 32 alr-MG that had a non-uniform or a uniform rg sign with all the ADG-traits, regardless of the rg magnitude, and used them to develop a microbiome-driven breeding strategy based on alr-MG only, or combined with ADG-traits, which was aimed at shaping the rumen microbiome towards increased ADG at all finishing stages. Combining alr-MG information with ADG records increased prediction accuracy of genomic estimated breeding values (GEBV) by 11 to 22% relative to the direct breeding strategy (using ADG-traits only), whereas using microbiome information, only, achieved lower accuracies (from 7 to 41%). Predicted selection responses varied consistently with accuracies. Restricting alr-MG based on their rg sign (uniform subset) did not yield a gain in the predicted response compared to the non-uniform subset, which is explained by the absence of alr-MG showing non-zero rg at least with more than one of the ADG-traits.

CONCLUSIONS

Our work sheds light on the role of the microbial metabolism in the growth trajectory of beef cattle at the genomic level and provides insights into the potential benefits of using microbiome information in future genomic breeding programs to accurately estimate GEBV and increase ADG at each finishing stage in beef cattle.

BASALT refines binning from metagenomic data and increases resolution of genome-resolved metagenomic analysis

Metagenomic binning is an essential technique for genome-resolved characterization of uncultured microorganisms in various ecosystems but hampered by the low efficiency of binning tools in adequately recovering metagenome-assembled genomes (MAGs). Here, we introduce BASALT (Binning Across a Series of Assemblies Toolkit) for binning and refinement of short- and long-read sequencing data. BASALT employs multiple binners with multiple thresholds to produce initial bins, then utilizes neural networks to identify core sequences to remove redundant bins and refine non-redundant bins. Using the same assemblies generated from Critical Assessment of Metagenome Interpretation (CAMI) datasets, BASALT produces up to twice as many MAGs as VAMB, DASTool, or metaWRAP. Processing assemblies from a lake sediment dataset, BASALT produces ~30% more MAGs than metaWRAP, including 21 unique class-level prokaryotic lineages. Functional annotations reveal that BASALT can retrieve 47.6% more non-redundant opening-reading frames than metaWRAP. These results highlight the robust handling of metagenomic sequencing data of BASALT.

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.

Red pandas with different diets and environments exhibit different gut microbial functional composition and capacity

The red panda (Ailurus fulgens) is a distinctive mammal known for its reliance on a diet primarily consisting of bamboo. The gut microbiota and overall health of animals are strongly influenced by diets and environments. Therefore, conducting research to explore the taxonomical and functional variances within the gut microbiota of red pandas exposed to various dietary and environmental conditions could shed light on the dynamic complexities of their microbial communities. In this study, normal fecal samples were obtained from red pandas residing in captive and semi-free environments under different dietary regimes and used for metabolomic, 16S rRNA, and metagenomic sequencing analysis, with the pandas classified into four distinct cohorts according to diet and environment. In addition, metagenomic sequencing was conducted on mucus fecal samples to elucidate potential etiological agents of disease. Results revealed an increased risk of gastrointestinal diseases in red pandas consuming bamboo shoots due to the heightened presence of pathogenic bacteria, although an increased presence of microbiota-derived tryptophan metabolites appeared to facilitate intestinal balance. The red pandas fed bamboo leaves also exhibited a decrease in gut microbial diversity, which may be attributed to the antibacterial flavonoids and lower protein levels in leaves. Notably, red pandas residing in semi-free environments demonstrated an enriched gut microbial diversity. Moreover, the occurrence of mucus secretion may be due to an increased presence of species associated with diarrhea and a reduced level of microbiota-derived tryptophan metabolites. In summary, our findings substantiate the influential role of diet and environment in modulating the gut microbiota of red pandas, offering potential implications for improved captive breeding practices.

Exploration of genes encoding KEGG pathway enzymes in rhizospheric microbiome of the wild plant Abutilon fruticosum

The operative mechanisms and advantageous synergies existing between the rhizobiome and the wild plant species Abutilon fruticosum were studied. Within the purview of this scientific study, the reservoir of genes in the rhizobiome, encoding the most highly enriched enzymes, was dominantly constituted by members of phylum Thaumarchaeota within the archaeal kingdom, phylum Proteobacteria within the bacterial kingdom, and the phylum Streptophyta within the eukaryotic kingdom. The ensemble of enzymes encoded through plant exudation exhibited affiliations with 15 crosstalking KEGG (Kyoto Encyclopaedia of Genes and Genomes) pathways. The ultimate goal underlying root exudation, as surmised from the present investigation, was the biosynthesis of saccharides, amino acids, and nucleic acids, which are imperative for the sustenance, propagation, or reproduction of microbial consortia. The symbiotic companionship existing between the wild plant and its associated rhizobiome amplifies the resilience of the microbial community against adverse abiotic stresses, achieved through the orchestration of ABA (abscisic acid) signaling and its cascading downstream effects. Emergent from the process of exudation are pivotal bioactive compounds including ATP, D-ribose, pyruvate, glucose, glutamine, and thiamine diphosphate. In conclusion, we hypothesize that future efforts to enhance the growth and productivity of commercially important crop plants under both favorable and unfavorable environmental conditions may focus on manipulating plant rhizobiomes.

- Invited Review - Understanding the functionality of the rumen microbiota: searching for better opportunities for rumen microbial manipulation

Rumen microbiota play a central role in the digestive process of ruminants. Their remarkable ability to break down complex plant fibers and proteins, converting them into essential organic compounds that provide animals with energy and nutrition. Research on rumen microbiota not only contributes to improving animal production performance and enhancing feed utilization efficiency but also holds the potential to reduce methane emissions and environmental impact. Nevertheless, studies on rumen microbiota face numerous challenges, including complexity, difficulties in cultivation, and obstacles in functional analysis. This review provides an overview of microbial species involved in the degradation of macromolecules, the fermentation processes, and methane production in the rumen, all based on cultivation methods. Additionally, the review introduces the applications, advantages, and limitations of emerging omics technologies such as metagenomics, metatranscriptomics, metaproteomics, and metabolomics, in investigating the functionality of rumen microbiota. Finally, the article offers a forward-looking perspective on the new horizons and technologies in the field of rumen microbiota functional research. These emerging technologies, with continuous refinement and mutual complementation, have deepened our understanding of rumen microbiota functionality, thereby enabling effective manipulation of the rumen microbial community.

Multi-omics revealed the effects of dietary energy levels on the rumen microbiota and metabolites in yaks under house-feeding conditions

Yak (Bos grunniens) is a unique large ruminant species in the Qinghai-Tibetan Plateau (QTP). Changing the energy levels of their rations can significantly improve their growth performance. Therefore, studying the effects of dietary energy levels on the rumen microflora and metabolites of yak is crucial for enhancing the development of the yak industry. Currently, there is a lack of understanding regarding the impact of feeding energy diets on rumen fermentation parameters, microbial functions, and metabolites. This study was designed to determine the appropriate energy level for feeding yak. Three test diets with metabolizable energy levels of 7.57 MJ/kg, 9.44 MJ/kg, and 11.9 MJ/kg were used and the concentration of volatile fatty acids (VFA) in rumen fluid was measured. The microbial communities, functions, and metabolites in yaks were studied by 16S rRNA sequencing, metagenome, and LC-MS non-targeted metabolomics to investigate the relationships among rumen fermentation parameters, microbial diversity, and metabolites. Ration energy levels significantly affect total VFA, acetate, propionate, butyrate, iso-valerate, valerate, and acetate/propionate (p < 0.05). At the phylum level, the dominant phyla in all three treatment groups were Bacteroidota, Firmicutes, and Actinobacteriota. At the genus level, the abundance of the unclassified_o__Bacteroidales, norank_f_Muribaculaceae, Lachnospiraceae_NK4A136_group, and Family _XIII_AD3011_group showed significant differences (p < 0.05) and were significantly correlated with differential metabolites screened for phosphatidylcholine [PC(16:0/0:0), PC(18:3/0:0)], uridine 3'-monophosphate, and adenosine monophosphate, etc. CAZymes family analysis showed that GHs and CEs differed significantly among the three groups. In addition, differential metabolites were mainly enriched in the pathways of lipid metabolism, nucleotide metabolism, and biosynthesis of other secondary metabolites, and the concentrations of differential metabolites were correlated with microbial abundance. In summary, this study analyzed the effects of ration energy levels on rumen microorganisms and metabolites of yaks and their relationships. The results provided a scientific basis for the selection of dietary energy for yaks in the house feeding period in the future.

Temporal colonization and metabolic regulation of the gut microbiome in neonatal oxen at single nucleotide resolution

The colonization of microbes in the gut is key to establishing a healthy host-microbiome symbiosis for newborns. We longitudinally profiled the gut microbiome in a model consisting of 36 neonatal oxen from birth up to 2 months postpartum and carried out microbial transplantation to reshape their gut microbiome. Genomic reconstruction of deeply sequenced fecal samples resulted in a total of 3931 metagenomic-assembled genomes from 472 representative species, of which 184 were identified as new species when compared with existing databases of oxen. Single nucleotide level metagenomic profiling shows a rapid influx of microbes after birth, followed by dynamic shifts during the first few weeks of life. Microbial transplantation was found to reshape the genetic makeup of 33 metagenomic-assembled genomes (FDR < 0.05), mainly from Prevotella and Bacteroides species. We further linked over 20 million microbial single nucleotide variations to 736 plasma metabolites, which enabled us to characterize 24 study-wide significant associations (P < 4.4 × 10-9) that identify the potential microbial genetic regulation of host immune and neuro-related metabolites, including glutathione and L-dopa. Our integration analyses further revealed that microbial genetic variations may influence the health status and growth performance by modulating metabolites via structural regulation of their encoded proteins. For instance, we found that the albumin levels and total antioxidant capacity were correlated with L-dopa, which was determined by single nucleotide variations via structural regulations of metabolic enzymes. The current results indicate that temporal colonization and transplantation-driven strain replacement are crucial for newborn gut development, offering insights for enhancing newborn health and growth.

Ecological niches and assembly dynamics of diverse microbial consortia in the gastrointestine of goat kids

Goats are globally invaluable ruminants that balance food security and environmental impacts, and their commensal microbiome residing in the gastrointestinal tract (GIT) is associated with animal health and productivity. However, the reference genomes and functional repertoires of GIT microbes in goat kids have not been fully elucidated. Herein, we performed a comprehensive landscape survey of the GIT microbiome of goat kids using metagenomic sequencing and binning, spanning a dense sampling regime covering three gastrointestinal compartments spatially and five developmental ages temporally. We recovered 1002 high-quality metagenome-assembled genomes (termed the goat kid GIT microbial catalog [GKGMC]), 618 of which were novel. They encode more than 2.3 million nonredundant proteins, and represent a variety of carbohydrate-degrading enzymes and metabolic gene clusters. The GKGMC-enriched microbial taxa, particularly Sodaliphilus, expanded the microbial tree of life in goat kids. Using this GKGMC, we first deciphered the prevalence of fiber-degrading bacteria for carbohydrate decomposition in the rumen and colon, while the ileal microbiota specialized in the uptake and conversion of simple sugars. Moreover, GIT microorganisms were rapidly assembled after birth, and their carbohydrate metabolic adaptation occurred in three phases of progression. Finally, phytobiotics modified the metabolic cascades of the ileal microbiome, underpinned by the enrichment of Sharpea azabuensis and Olsenella spp. implicated in lactate formation and utilization. This GKGMC reference provides novel insights into the early-life microbial developmental dynamics in distinct compartments, and offers expanded resources for GIT microbiota-related research in goat kids.

Microbial composition play the leading role in volatile fatty acid production in the fermentation of different scale of corn stover with rumen fluid

Rumen fluid is a natural and green biocatalyst that can efficiently degrade biomass into volatile fatty acid (VFA) used to produce value-added materials. But the essence of high degradation efficiency in the rumen has not been fully analyzed. This study investigated the contribution of substrate structure and microbial composition to volatile fatty acid production in the fermentation of corn stover. The ball milled corn stover were innovatively applied to ferment with the rumen fluid collected at different digestion times. Exogeneous cellulase was also added to the ruminal fermentation to further reveal the inner mechanism. With prolonged digestion time, the microbial community relative abundance levels of Bacteroidetes and Firmicutes increased from 29.98% to 72.74% and decreased from 51.76% to 22.11%, respectively. The highest VFA production of the corn stover was achieved via treatment with the rumen fluid collected at 24 h which was up to 9508 mg/L. The ball milled corn stover achieved high VFA production because of the more accessible substrate structure. The application of exogenous cellulase has no significant influence to the ruminal fermentation. The microbial community abundance contributed more to the VFA production compared with the substrate structures.

Comparative Rumen Metagenome and CAZyme Profiles in Cattle and Buffaloes: Implications for Methane Yield and Rumen Fermentation on a Common Diet

A study was undertaken to compare the rumen microbial community composition, methane yield, rumen fermentation, and CAZyme profiles between cattle and buffaloes. The primary aim of this study was to ascertain the impact of the host species on the above when diet and environmental factors are fixed. A total of 43 phyla, 200 orders, 458 families, and 1722 microbial genera were identified in the study. Bacteroidetes was the most prominent bacterial phylum and constituted >1/3rd of the ruminal microbiota; however, their abundances were comparable between cattle and buffaloes. Firmicutes were the second most abundant bacteria, found to be negatively correlated with the Bacteroidetes. The abundances of Firmicutes as well as the F/B ratio were not different between the two host species. In this study, archaea affiliated with the nine phyla were identified, with Euryarchaeota being the most prominent. Like bacterial phyla, the abundances of Euryarchaeota methanogens were also similar between the cattle and buffaloes. At the order level, Methanobacteriales dominated the archaea. Methanogens from the Methanosarcinales, Methanococcales, Methanomicrobiales, and Methanomassiliicoccales groups were also identified, but at a lower frequency. Methanobrevibacter was the most prevalent genus of methanogens, accounting for approximately three percent of the rumen metagenome. However, their distribution was not different between the two host species. CAZymes affiliated with five classes, namely CBM, CE, GH, GT, and PL, were identified in the metagenome, where the GH class was the most abundant and constituted ~70% of the total CAZymes. The protozoal numbers, including Entodiniomorphs and Holotrichs, were also comparable between the cattle and buffaloes. Results from the study did not reveal any significant difference in feed intake, nutrient digestibility, and rumen fermentation between cattle and buffaloes fed on the same diet. As methane yield due to the similar diet composition, feed ingredients, rumen fermentation, and microbiota composition did not vary, these results indicate that the microbiota community structure and methane emissions are under the direct influence of the diet and environment, and the host species may play only a minor role until the productivity does not vary. More studies are warranted to investigate the effect of different diets and environments on microbiota composition and methane yield. Further, the impact of variable productivity on both the cattle and buffaloes when the diet and environmental factors are fixed needs to be ascertained.

Comparative analysis of functional diversity of rumen microbiome in bison and beef heifers

IMPORTANCE

Ruminants play a key role in the conversion of cellulolytic plant material into high-quality meat and milk protein for humans. The rumen microbiome is the driver of this conversion, yet there is little information on how gene expression within the microbiome impacts the efficiency of this conversion process. The current study investigates gene expression in the rumen microbiome of beef heifers and bison and how transplantation of ruminal contents from bison to heifers alters gene expression. Understanding interactions between the host and the rumen microbiome is the key to developing informed approaches to rumen programming that will enhance production efficiency in ruminants.

Effect of Dietary and Age Changes on Ruminal Microbial Diversity in Holstein Calves

Ruminal microorganisms play a crucial role in the energy supply of ruminants and animal performance. We analyzed the variations in rumen bacteria and fungi at 45 d, 75 d, and 105 d by using 16SrRNA and ITS sequencing data and investigated their correlation with rumen fermentation. According to the results, rumen microflora tended to gradually mature with age, and bacterial and fungal establishment gradually stabilized. Upon comparing the three periods, the concentration of propionic acid increased significantly (p < 0.05) after weaning, and weaning accompanied by a transition in diet remarkably decreased (p < 0.05) rumen diversity in the short term and induced a corresponding change in the rumen microbiota composition. Bacteroidota, Actinobacteriota, and Firmicutes were the core bacterial phyla for all age periods. Ruminococcus, NK4A214_group, Sharpea, Rikenellaceae_RC9_gut_group, and norank_f__Butyricicoccaceae were the markedly abundant bacterial genera in pre-weaning. After weaning, the relative abundance of Erysipelotrichaceae_ UCG-002, Eubacterium_ruminantium_group, and Solobacterium significantly increased (p < 0.05). The relative abundance of Acetitomaculum increased with age with the greatest abundance noted at 105 d (37%). The dominant fungal phyla were Ascomycota and Basidiomycota, and Aspergillus and Xeromyces were the most abundant fungal genera after weaning. Trichomonascus, Phialosimplex, and Talaromyces were enriched at 105 d. However, the low abundance of Neocallimastigomycota was not detected throughout the study, which is worthy of further investigation. In addition, correlations were observed between age-related abundances of specific genera and microbiota functions and rumen fermentation-related parameters. This study revealed that rumen microbiota and rumen fermentation capacity are correlated, which contributed to a better understanding of the effects of age and diet on rumen microbiology and fermentation in calves.

Exploring the microbial diversity and characterization of cellulase and hemicellulase genes in goat rumen: a metagenomic approach

BACKGROUND

Goat rumen microbial communities are perceived as one of the most potential biochemical reservoirs of multi-functional enzymes, which are applicable to enhance wide array of bioprocesses such as the hydrolysis of cellulose and hemi-cellulose into fermentable sugar for biofuel and other value-added biochemical production. Even though, the limited understanding of rumen microbial genetic diversity and the absence of effective screening culture methods have impeded the full utilization of these potential enzymes. In this study, we applied culture independent metagenomics sequencing approach to isolate, and identify microbial communities in goat rumen, meanwhile, clone and functionally characterize novel cellulase and xylanase genes in goat rumen bacterial communities.

RESULTS

Bacterial DNA samples were extracted from goat rumen fluid. Three genomic libraries were sequenced using Illumina HiSeq 2000 for paired-end 100-bp (PE100) and Illumina HiSeq 2500 for paired-end 125-bp (PE125). A total of 435gb raw reads were generated. Taxonomic analysis using Graphlan revealed that Fibrobacter, Prevotella, and Ruminococcus are the most abundant genera of bacteria in goat rumen. SPAdes assembly and prodigal annotation were performed. The contigs were also annotated using the DOE-JGI pipeline. In total, 117,502 CAZymes, comprising endoglucanases, exoglucanases, beta-glucosidases, xylosidases, and xylanases, were detected in all three samples. Two genes with predicted cellulolytic/xylanolytic activities were cloned and expressed in E. coli BL21(DE3). The endoglucanases and xylanase enzymatic activities of the recombinant proteins were confirmed using substrate plate assay and dinitrosalicylic acid (DNS) analysis. The 3D structures of endoglucanase A and endo-1,4-beta xylanase was predicted using the Swiss Model. Based on the 3D structure analysis, the two enzymes isolated from goat's rumen metagenome are unique with only 56-59% similarities to those homologous proteins in protein data bank (PDB) meanwhile, the structures of the enzymes also displayed greater stability, and higher catalytic activity.

CONCLUSIONS

In summary, this study provided the database resources of bacterial metagenomes from goat's rumen fluid, including gene sequences with annotated functions and methods for gene isolation and over-expression of cellulolytic enzymes; and a wealth of genes in the metabolic pathways affecting food and nutrition of ruminant animals.

Rumen multi-omics addressing diet-host-microbiome interplay in farm animals: a review

Continuous improvement in the living standards of developing countries, calls for an urgent need of high quality meat and dairy products. The farm animals have a micro-ecosystem in gastro-intestinal tract, comprising of a wide variety of flora and fauna which converts roughages and agricultural byproducts as well as nutrient rich concentrate sources into the useful products such as volatile fatty acids and microbial crude proteins. The microbial diversity changes according to composition of the feed, host species/breed and host's individual genetic makeup. From culture methods to next-generation sequencing technologies, the knowledge has emerged a lot to know-how of microbial world viz. their identification, enzymatic activities and metabolites which are the keys of ruminant's successful existence. The structural composition of ruminal community revealed through metagenomics can be elaborated by metatranscriptomics and metabolomics through deciphering their functional role in metabolism and their responses to the external and internal stimuli. These highly sophisticated analytical tools have made possible to correlate the differences in the feed efficiency, nutrients utilization and methane emissions to their rumen microbiome. The comprehensively understood rumen microbiome will enhance the knowledge in the fields of animal nutrition, biotechnology and climatology through deciphering the significance of each and every domain of residing microbial entity. The present review undertakes the recent investigations regarding rumen multi-omics viz. taxonomic and functional potential of microbial populations, host-diet-microbiome interactions and correlation with metabolic dynamics.

Review: Towards the next-generation models of the rumen microbiome for enhancing predictive power and guiding sustainable production strategies

The rumen ecosystem harbours a galaxy of microbes working in syntrophy to carry out a metabolic cascade of hydrolytic and fermentative reactions. This fermentation process allows ruminants to harvest nutrients from a wide range of feedstuff otherwise inaccessible to the host. The interconnection between the ruminant and its rumen microbiota shapes key animal phenotypes such as feed efficiency and methane emissions and suggests the potential of reducing methane emissions and enhancing feed conversion into animal products by manipulating the rumen microbiota. Whilst significant technological progress in omics techniques has increased our knowledge of the rumen microbiota and its genome (microbiome), translating omics knowledge into effective microbial manipulation strategies remains a great challenge. This challenge can be addressed by modelling approaches integrating causality principles and thus going beyond current correlation-based approaches applied to analyse rumen microbial genomic data. However, existing rumen models are not yet adapted to capitalise on microbial genomic information. This gap between the rumen microbiota available omics data and the way microbial metabolism is represented in the existing rumen models needs to be filled to enhance rumen understanding and produce better predictive models with capabilities for guiding nutritional strategies. To fill this gap, the integration of computational biology tools and mathematical modelling frameworks is needed to translate the information of the metabolic potential of the rumen microbes (inferred from their genomes) into a mathematical object. In this paper, we aim to discuss the potential use of two modelling approaches for the integration of microbial genomic information into dynamic models. The first modelling approach explores the theory of state observers to integrate microbial time series data into rumen fermentation models. The second approach is based on the genome-scale network reconstructions of rumen microbes. For a given microorganism, the network reconstruction produces a stoichiometry matrix of the metabolism. This matrix is the core of the so-called genome-scale metabolic models which can be exploited by a plethora of methods comprised within the constraint-based reconstruction and analysis approaches. We will discuss how these methods can be used to produce the next-generation models of the rumen microbiome.

Effect of castration timing and weaning strategy on the taxonomic and functional profile of ruminal bacteria and archaea of beef calves

BACKGROUND

Beef cattle experience several management challenges across their lifecycle. Castration and weaning, two major interventions in the early life of beef cattle, can have a substantial impact on animal performance. Despite the key role of the rumen microbiome on productive traits of beef cattle, the effect of castration timing and weaning strategy on this microbial community has not been formally described. We assessed the effect of four castration time windows (at birth, turnout, pre-weaning and weaning) and two weaning strategies (fence-line and truck transportation) on the rumen microbiome in a randomized controlled study with 32 male calves across 3 collection days (i.e., time points). Ruminal fluid samples were submitted to shotgun metagenomic sequencing and changes in the taxonomic (microbiota) and functional profile (metagenome) of the rumen microbiome were described.

RESULTS

Using a comprehensive yet stringent taxonomic classification approach, we identified 10,238 unique taxa classified under 40 bacterial and 7 archaeal phyla across all samples. Castration timing had a limited long-term impact on the rumen microbiota and was not associated with changes in alpha and beta diversity. The interaction of collection day and weaning strategy was associated with changes in the rumen microbiota, which experienced a significant decrease in alpha diversity and shifts in beta diversity within 48 h post-weaning, especially in calves abruptly weaned by truck transportation. Calves weaned using a fence-line weaning strategy had lower relative abundance of Bacteroides, Lachnospira, Fibrobacter and Ruminococcus genera compared to calves weaned by truck transportation. Some genes involved in the hydrogenotrophic methanogenesis pathway (fwdB and fwdF) had higher relative abundance in fence-line-weaned calves post-weaning. The antimicrobial resistance gene tetW consistently represented more than 50% of the resistome across time, weaning and castration groups, without significant changes in relative abundance.

CONCLUSIONS

Within the context of this study, castration timing had limited long-term effects on the rumen microbiota, while weaning strategy had short-term effects on the rumen microbiota and methane-associated metagenome, but not on the rumen resistome.

Carbohydrate flow through agricultural ecosystems: Implications for synthesis and microbial conversion of carbohydrates

Carbohydrates are chemically and structurally diverse biomolecules, serving numerous and varied roles in agricultural ecosystems. Crops and horticulture products are inherent sources of carbohydrates that are consumed by humans and non-human animals alike; however carbohydrates are also present in other agricultural materials, such as soil and compost, human and animal tissues, milk and dairy products, and honey. The biosynthesis, modification, and flow of carbohydrates within and between agricultural ecosystems is intimately related with microbial communities that colonize and thrive within these environments. Recent advances in -omics techniques have ushered in a new era for microbial ecology by illuminating the functional potential for carbohydrate metabolism encoded within microbial genomes, while agricultural glycomics is providing fresh perspective on carbohydrate-microbe interactions and how they influence the flow of functionalized carbon. Indeed, carbohydrates and carbohydrate-active enzymes are interventions with unrealized potential for improving carbon sequestration, soil fertility and stability, developing alternatives to antimicrobials, and circular production systems. In this manner, glycomics represents a new frontier for carbohydrate-based biotechnological solutions for agricultural systems facing escalating challenges, such as the changing climate.

Hi-C metagenome sequencing reveals soil phage-host interactions

Bacteriophages are abundant in soils. However, the majority are uncharacterized, and their hosts are unknown. Here, we apply high-throughput chromosome conformation capture (Hi-C) to directly capture phage-host relationships. Some hosts have high centralities in bacterial community co-occurrence networks, suggesting phage infections have an important impact on the soil bacterial community interactions. We observe increased average viral copies per host (VPH) and decreased viral transcriptional activity following a two-week soil-drying incubation, indicating an increase in lysogenic infections. Soil drying also alters the observed phage host range. A significant negative correlation between VPH and host abundance prior to drying indicates more lytic infections result in more host death and inversely influence host abundance. This study provides empirical evidence of phage-mediated bacterial population dynamics in soil by directly capturing specific phage-host interactions.

Gastrointestinal microbiome, resistance genes, and risk assessment of heavy metals in wild giant pandas

The gastrointestinal microbiome (GM) of giant panda (GP) plays an important role in food utilization and health and is also an essential reservoir of resistance genes. Currently, little knowledge is available on the GM, acid resistance genes (AcRGs), antibiotic resistance genes (ARGs), metal resistance genes (MRGs), and mobile genetic elements (MGEs) in wild GPs. We sampled the gastrointestinal tract of a dead GP and explored the composition and function of GM and resistance genes through cryo-scanning electron microscopy, metagenomic sequencing, and genome-resolved metagenomics. The concentration of metals in the gastrointestinal lumen, feces, bamboo, and soil was measured by inductively coupled plasma mass spectrometry. Results showed that the composition of the microbiota varied in different gastrointestinal regions. Fecal microbiota was highly associated with small intestinal and colonic microbes. The lignocellulosic cross-linked structure of bamboo was destroyed in the stomach initially and destroying degree increased from stomach to anus. Reconstruction of metagenome-assembled-genomes confirmed that core GM, e.g., Streptococcus, Clostridium, Lactococcus, Leuconostoc, and Enterococcus, carried genes encoding the lignocellulose degradation enzyme. There were no significant differences of resistance genes between gastrointestinal and fecal samples, except MGEs. Multidrug and multi-metal resistance genes were predominant in all samples, while the transposase gene tnpA was the major type of MGE. Significant correlations were observed among the abundance of GM, resistance genes, and MGEs. Gastrointestinal and fecal mercury and chromium were the main metals influencing GM and resistance genes. The content of gastrointestinal and fecal metals was significantly associated with the presence of the same metals in bamboo, which could pose a threat to the health of wild GPs. This study characterized the gastrointestinal microbiome of wild GPs, providing new evidence for the role of the gastrointestinal microbiome in degrading lignocellulose from bamboo and highlighting the urgent need to monitor metal levels in soil and bamboo.

Evolving understanding of rumen methanogen ecophysiology

Production of methane by methanogenic archaea, or methanogens, in the rumen of ruminants is a thermodynamic necessity for microbial conversion of feed to volatile fatty acids, which are essential nutrients for the animals. On the other hand, methane is a greenhouse gas and its production causes energy loss for the animal. Accordingly, there are ongoing efforts toward developing effective strategies for mitigating methane emissions from ruminant livestock that require a detailed understanding of the diversity and ecophysiology of rumen methanogens. Rumen methanogens evolved from free-living autotrophic ancestors through genome streamlining involving gene loss and acquisition. The process yielded an oligotrophic lifestyle, and metabolically efficient and ecologically adapted descendants. This specialization poses serious challenges to the efforts of obtaining axenic cultures of rumen methanogens, and consequently, the information on their physiological properties remains in most part inferred from those of their non-rumen representatives. This review presents the current knowledge of rumen methanogens and their metabolic contributions to enteric methane production. It also identifies the respective critical gaps that need to be filled for aiding the efforts to mitigate methane emission from livestock operations and at the same time increasing the productivity in this critical agriculture sector.

Transcriptional control of carbohydrate catabolism by the CcpA protein in the ruminal bacterium Streptococcus bovis

As a potent, pleiotropic regulatory protein in Gram-positive bacteria, catabolite control protein A (CcpA) mediates the transcriptional control of carbohydrate metabolism in Streptococcus bovis, a lactate-producing bacterium that plays an essential role in rumen acidosis in dairy cows. Although the rumen uptake of carbohydrates is multi-substrate, the focus of S. bovis research thus far has been on the glucose. With the aid of gene deletion, whole-genome sequencing, and transcriptomics, we have unraveled the role of CcpA in carbohydrate metabolism, on the one hand, and acidosis, on the other, and we show that the S. bovis strain S1 encodes "Carbohydrate-Active Enzymes" and that ccpA deletion slows the organism's growth rate and modulates the organic acid fermentation pathways toward lower lactate, higher formate, and acetate in the maltose and cellobiose. Furthermore, this study revealed the different regulatory functions of the CcpA protein in rumen metabolism and acidosis.IMPORTANCEThis study is important as it illustrates the varying regulatory role of the Streptococcus bovis catabolite control protein A protein in carbohydrate metabolism and the onset of acidosis in dairy cattle.

Recovery of 52 bacterial genomes from the fecal microbiome of the domestic cat (Felis catus) using Hi-C proximity ligation and shotgun metagenomics

We used Hi-C proximity ligation with shotgun sequencing to retrieve metagenome-assembled genomes (MAGs) from the fecal microbiomes of two domestic cats (Felis catus). The genomes were assessed for completeness and contamination, classified taxonomically, and annotated for putative antimicrobial resistance (AMR) genes.

The hologenome of Daphnia magna reveals possible DNA methylation and microbiome-mediated evolution of the host genome

Properties that make organisms ideal laboratory models in developmental and medical research are often the ones that also make them less representative of wild relatives. The waterflea Daphnia magna is an exception, by both sharing many properties with established laboratory models and being a keystone species, a sentinel species for assessing water quality, an indicator of environmental change and an established ecotoxicology model. Yet, Daphnia's full potential has not been fully exploited because of the challenges associated with assembling and annotating its gene-rich genome. Here, we present the first hologenome of Daphnia magna, consisting of a chromosomal-level assembly of the D. magna genome and the draft assembly of its metagenome. By sequencing and mapping transcriptomes from exposures to environmental conditions and from developmental morphological landmarks, we expand the previously annotates gene set for this species. We also provide evidence for the potential role of gene-body DNA-methylation as a mutagen mediating genome evolution. For the first time, our study shows that the gut microbes provide resistance to commonly used antibiotics and virulence factors, potentially mediating Daphnia's environmental-driven rapid evolution. Key findings in this study improve our understanding of the contribution of DNA methylation and gut microbiota to genome evolution in response to rapidly changing environments.

PandaGUT provides new insights into bacterial diversity, function, and resistome landscapes with implications for conservation

BACKGROUND

The gut microbiota play important roles in host adaptation and evolution, but are understudied in natural population of wild mammals. To address host adaptive evolution and improve conservation efforts of threatened mammals from a metagenomic perspective, we established a high-quality gut microbiome catalog of the giant panda (pandaGUT) to resolve the microbiome diversity, functional, and resistome landscapes using approximately 7 Tbp of long- and short-read sequencing data from 439 stool samples.

RESULTS

The pandaGUT catalog comprises 820 metagenome-assembled genomes, including 40 complete closed genomes, and 64.5% of which belong to species that have not been previously reported, greatly expanding the coverage of most prokaryotic lineages. The catalog contains 2.37 million unique genes, with 74.8% possessing complete open read frames, facilitating future mining of microbial functional potential. We identified three microbial enterotypes across wild and captive panda populations characterized by Clostridium, Pseudomonas, and Escherichia, respectively. We found that wild pandas exhibited host genetic-specific microbial structures and functions, suggesting host-gut microbiota phylosymbiosis, while the captive cohorts encoded more multi-drug resistance genes.

CONCLUSIONS

Our study provides largely untapped resources for biochemical and biotechnological applications as well as potential intervention avenues via the rational manipulation of microbial diversity and reducing antibiotic usage for future conservation management of wildlife. Video Abstract.

MetaCC allows scalable and integrative analyses of both long-read and short-read metagenomic Hi-C data

Metagenomic Hi-C (metaHi-C) can identify contig-to-contig relationships with respect to their proximity within the same physical cell. Shotgun libraries in metaHi-C experiments can be constructed by next-generation sequencing (short-read metaHi-C) or more recent third-generation sequencing (long-read metaHi-C). However, all existing metaHi-C analysis methods are developed and benchmarked on short-read metaHi-C datasets and there exists much room for improvement in terms of more scalable and stable analyses, especially for long-read metaHi-C data. Here we report MetaCC, an efficient and integrative framework for analyzing both short-read and long-read metaHi-C datasets. MetaCC outperforms existing methods on normalization and binning. In particular, the MetaCC normalization module, named NormCC, is more than 3000 times faster than the current state-of-the-art method HiCzin on a complex wastewater dataset. When applied to one sheep gut long-read metaHi-C dataset, MetaCC binning module can retrieve 709 high-quality genomes with the largest species diversity using one single sample, including an expansion of five uncultured members from the order Erysipelotrichales, and is the only binner that can recover the genome of one important species Bacteroides vulgatus. Further plasmid analyses reveal that MetaCC binning is able to capture multi-copy plasmids.

The phenotype and genotype of fermentative prokaryotes

Fermentation is a type of metabolism pervasive in oxygen-deprived environments. Despite its importance, we know little about the range and traits of organisms that carry out this metabolism. Our study addresses this gap with a comprehensive analysis of the phenotype and genotype of fermentative prokaryotes. We assembled a dataset with phenotypic records of 8350 organisms plus 4355 genomes and 13.6 million genes. Our analysis reveals fermentation is both widespread (in ~30% of prokaryotes) and complex (forming ~300 combinations of metabolites). Furthermore, it points to previously uncharacterized proteins involved in this metabolism. Previous studies suggest that metabolic pathways for fermentation are well understood, but metabolic models built in our study show gaps in our knowledge. This study demonstrates the complexity of fermentation while showing that there is still much to learn about this metabolism. All resources in our study can be explored by the scientific community with an online, interactive tool.

Holistic view of the seascape dynamics and environment impact on macro-scale genetic connectivity of marine plankton populations

BACKGROUND

Plankton seascape genomics studies have revealed different trends from large-scale weak differentiation to microscale structures. Previous studies have underlined the influence of the environment and seascape on species differentiation and adaptation. However, these studies have generally focused on a few single species, sparse molecular markers, or local scales. Here, we investigated the genomic differentiation of plankton at the macro-scale in a holistic approach using Tara Oceans metagenomic data together with a reference-free computational method.

RESULTS

We reconstructed the FST-based genomic differentiation of 113 marine planktonic taxa occurring in the North and South Atlantic Oceans, Southern Ocean, and Mediterranean Sea. These taxa belong to various taxonomic clades spanning Metazoa, Chromista, Chlorophyta, Bacteria, and viruses. Globally, population genetic connectivity was significantly higher within oceanic basins and lower in bacteria and unicellular eukaryotes than in zooplankton. Using mixed linear models, we tested six abiotic factors influencing connectivity, including Lagrangian travel time, as proxies of oceanic current effects. We found that oceanic currents were the main population genetic connectivity drivers, together with temperature and salinity. Finally, we classified the 113 taxa into parameter-driven groups and showed that plankton taxa belonging to the same taxonomic rank such as phylum, class or order presented genomic differentiation driven by different environmental factors.

CONCLUSION

Our results validate the isolation-by-current hypothesis for a non-negligible proportion of taxa and highlight the role of other physicochemical parameters in large-scale plankton genetic connectivity. The reference-free approach used in this study offers a new systematic framework to analyse the population genomics of non-model and undocumented marine organisms from a large-scale and holistic point of view.

Interrogating the viral dark matter of the rumen ecosystem with a global virome database

The diverse rumen virome can modulate the rumen microbiome, but it remains largely unexplored. Here, we mine 975 published rumen metagenomes for viral sequences, create a global rumen virome database (RVD), and analyze the rumen virome for diversity, virus-host linkages, and potential roles in affecting rumen functions. Containing 397,180 species-level viral operational taxonomic units (vOTUs), RVD substantially increases the detection rate of rumen viruses from metagenomes compared with IMG/VR V3. Most of the classified vOTUs belong to Caudovirales, differing from those found in the human gut. The rumen virome is predicted to infect the core rumen microbiome, including fiber degraders and methanogens, carries diverse auxiliary metabolic genes, and thus likely impacts the rumen ecosystem in both a top-down and a bottom-up manner. RVD and the findings provide useful resources and a baseline framework for future research to investigate how viruses may impact the rumen ecosystem and digestive physiology.

MT-MAG: Accurate and interpretable machine learning for complete or partial taxonomic assignments of metagenomeassembled genomes

We propose MT-MAG, a novel machine learning-based software tool for the complete or partial hierarchically-structured taxonomic classification of metagenome-assembled genomes (MAGs). MT-MAG is alignment-free, with k-mer frequencies being the only feature used to distinguish a DNA sequence from another (herein k = 7). MT-MAG is capable of classifying large and diverse metagenomic datasets: a total of 245.68 Gbp in the training sets, and 9.6 Gbp in the test sets analyzed in this study. In addition to complete classifications, MT-MAG offers a "partial classification" option, whereby a classification at a higher taxonomic level is provided for MAGs that cannot be classified to the Species level. MT-MAG outputs complete or partial classification paths, and interpretable numerical classification confidences of its classifications, at all taxonomic ranks. To assess the performance of MT-MAG, we define a "weighted classification accuracy," with a weighting scheme reflecting the fact that partial classifications at different ranks are not equally informative. For the two benchmarking datasets analyzed (genomes from human gut microbiome species, and bacterial and archaeal genomes assembled from cow rumen metagenomic sequences), MT-MAG achieves an average of 87.32% in weighted classification accuracy. At the Species level, MT-MAG outperforms DeepMicrobes, the only other comparable software tool, by an average of 34.79% in weighted classification accuracy. In addition, MT-MAG is able to completely classify an average of 67.70% of the sequences at the Species level, compared with DeepMicrobes which only classifies 47.45%. Moreover, MT-MAG provides additional information for sequences that it could not classify at the Species level, resulting in the partial or complete classification of 95.13%, of the genomes in the datasets analyzed. Lastly, unlike other taxonomic assignment tools (e.g., GDTB-Tk), MT-MAG is an alignment-free and genetic marker-free tool, able to provide additional bioinformatics analysis to confirm existing or tentative taxonomic assignments.

The unique gut microbiome of giant pandas involved in protein metabolism contributes to the host's dietary adaption to bamboo

BACKGROUND

The gut microbiota of the giant panda (Ailuropoda melanoleuca), a global symbol of conservation, are believed to be involved in the host's dietary switch to a fibrous bamboo diet. However, their exact roles are still largely unknown.

RESULTS

In this study, we first comprehensively analyzed a large number of gut metagenomes giant pandas (n = 322), including 98 pandas sequenced in this study with deep sequencing (Illumina) and third-generation sequencing (nanopore). We reconstructed 408 metagenome-assembled genomes (MAGs), and 148 of which (36.27%) were near complete. The most abundant MAG was classified as Streptococcus alactolyticus. A pairwise comparison of the metagenomes and meta-transcriptomes in 14 feces revealed genes involved in carbohydrate metabolism were lower, but those involved in protein metabolism were greater in abundance and expression in giant pandas compared to those in herbivores and omnivores. Of note, S. alactolyticus was positively correlated to the KEGG modules of essential amino-acid biosynthesis. After being isolated from pandas and gavaged to mice, S. alactolyticus significantly increased the relative abundance of essential amino acids in mice jejunum.

CONCLUSIONS

The study highlights the unique protein metabolic profiles in the giant panda's gut microbiome. The findings suggest that S. alactolyticus is an important player in the gut microbiota that contributes to the giant panda's dietary adaptation by more involvement in protein rather than carbohydrate metabolism. Video Abstract.

A comparison of single-coverage and multi-coverage metagenomic binning reveals extensive hidden contamination

Metagenomic binning has revolutionized the study of uncultured microorganisms. Here we compare single- and multi-coverage binning on the same set of samples, and demonstrate that multi-coverage binning produces better results than single-coverage binning and identifies contaminant contigs and chimeric bins that other approaches miss. While resource expensive, multi-coverage binning is a superior approach and should always be performed over single-coverage binning.

Metagenomic approach to infer rumen microbiome derived traits of cattle

Ruminants enable the conversion of indigestible plant material into animal consumables, including dairy products, meat, and valuable fibers. Microbiome research is gaining popularity in livestock species because it aids in the knowledge of illnesses and efficiency processes in animals. In this study, we use WGS metagenomic data to thoroughly characterize the ruminal ecosystem of cows to infer positive and negative livestock traits determined by the microbiome. The rumen of cows from Argentina were described by combining different gene biomarkers, pathways composition and taxonomic information. Taxonomic characterization indicated that the two major phyla were Bacteroidetes and Firmicutes; in third place, Proteobacteria was highly represented followed by Actinobacteria; Prevotella, and Bacteroides were the most abundant genera. Functional profiling of carbohydrate-active enzymes indicated that members of the Glycoside Hydrolase (GH) class accounted for 52.2 to 55.6% of the total CAZymes detected, among them the most abundant were the oligosaccharide degrading enzymes. The diversity of GH families found suggested efficient hydrolysis of complex biomass. Genes of multidrug, macrolides, polymyxins, beta-lactams, rifamycins, tetracyclines, and bacitracin resistance were found below 0.12% of relative abundance. Furthermore, the clustering analysis of genera and genes that correlated to methane emissions or feed efficiency, suggested that the cows analysed could be regarded as low methane emitters and clustered with high feed efficiency reference animals. Finally, the combination of bioinformatic analyses used in this study can be applied to assess cattle traits difficult to measure and guide enhanced nutrition and breeding methods.

Expanding the rumen Prevotella collection: The description of Prevotella communis, sp. nov. of ovine origin

27 strains representing eight new Prevotella species were isolated from rumen of a single sheep in eight weeks interval. One of the putative species encompassing the highest number of isolated strains which also exhibited some genetic variability in preliminary data, was then selected for description of a novel species. We examined six strains in genomic and phenotypic detail, two of which may actually be the same strain isolated nearly three weeks apart. Other strains formed clearly diverged intraspecies lineages as evidenced by core genome phylogeny and phenotypic differences. Strains of the proposed new Prevotella species are strictly saccharolytic as is usual for rumen Prevotella, and use plant cell-wall xylans and pectins for growth. However, the range of cell-wall polysaccharides utilised for growth is rather limited compared to rumen generalists such as Prevotella bryantii or Prevotella ruminicola and this extends also to the inability to utilise starch, which is unexpected for the members of the genus Prevotella. Based on the data obtained, we propose Prevotella communis sp. nov. to accommodate strain E1-9T as well as other strains with the similar properties. The proposed species is widespread: two other strains were previously isolated from sheep in Japan and is also common in metagenomic data of cattle and sheep rumen samples from Scotland and New Zealand. It was also found in a collection of metagenome-assembled genomes originating from cattle in Scotland. Thus, it is a ubiquitous bacterium of domesticated ruminants specialising in degradation of a somewhat restricted set of plant cell wall components.

Age-related alterations in metabolome and microbiome provide insights in dietary transition in giant pandas

We conducted UPLC-MS-based metabolomics, 16S rRNA, and metagenome sequencing on the fecal samples of 44 captive giant pandas (Ailuropoda melanoleuca) from four age groups (i.e., Cub, Young, Adult, and Old) to comprehensively understand age-related changes in the metabolism and gut microbiota of giant pandas. We characterized the metabolite profiles of giant pandas based on 1,376 identified metabolites, with 152 significantly differential metabolites (SDMs) found across the age groups. We found that the metabolites and the composition/function of the gut microbiota changed in response to the transition from a milk-dominant diet in panda cubs to a bamboo-specific diet in young and adult pandas. Lipid metabolites such as choline and hippuric acid were enriched in the Cub group, and many plant secondary metabolites were significantly higher in the Young and Adult groups, while oxidative stress and inflammatory related metabolites were only found in the Old group. However, there was a decrease in the α-diversity of gut microbiota in adult and old pandas, who exclusively consume bamboo. The abundance of bacteria related to the digestion of cellulose-rich food, such as Firmicutes, Streptococcus, and Clostridium, significantly increased from the Cub to the Adult group, while the abundance of beneficial bacteria such as Faecalibacterium, Sarcina, and Blautia significantly decreased. Notably, several potential pathogenic bacteria had relatively high abundances, especially in the Young group. Metagenomic analysis identified 277 CAZyme genes including cellulose degrading genes, and seven of the CAZymes had abundances that significantly differed between age groups. We also identified 237 antibiotic resistance genes (ARGs) whose number and diversity increased with age. We also found a significant positive correlation between the abundance of bile acids and gut bacteria, especially Lactobacillus and Bifidobacterium. Our results from metabolome, 16S rRNA, and metagenome data highlight the important role of the gut microbiota-bile acid axis in the regulation of age-related metabolism and provide new insights into the lipid metabolism of giant pandas. IMPORTANCE The giant panda is a member of the order Carnivora but is entirely herbivorous. The giant panda's specialized diet and related metabolic mechanisms have not been fully understood. It is therefore crucial to investigate the dynamic changes in metabolites as giant pandas grow and physiologically adapt to their herbivorous diet. This study conducted UPLC-MS-based metabolomics 16S rRNA, and metagenome sequencing on the fecal samples of captive giant pandas from four age groups. We found that metabolites and the composition/function of gut microbiota changed in response to the transition from a milk-dominant diet in cubs to a bamboo-specific diet in young and adult pandas. The metabolome, 16S rRNA, and metagenome results highlight that the gut microbiota-bile acid axis has an important role in the regulation of age-related metabolism, and our study provides new insights into the lipid metabolism of giant pandas.

Metagenomic analysis reveals the efficient digestion mechanism of corn stover in Angus bull rumen: Microbial community succession, CAZyme composition and functional gene expression

Ruminant rumen is a biological fermentation system that can efficiently degrade lignocellulosic biomass. The knowledge about mechanisms of efficient lignocellulose degradation with rumen microorganisms is still limited. In this study, composition and succession of bacteria and fungi, carbohydrate-active enzymes (CAZymes), and functional genes involved in hydrolysis and acidogenesis were revealed during fermentation in Angus bull rumen via metagenomic sequencing. Results showed that degradation efficiency of hemicellulose and cellulose reached 61.2% and 50.4% at 72 h fermentation, respectively. Main bacterial genera were composed of Prevotella, Butyrivibrio, Ruminococcus, Eubacterium, and Fibrobacter, and main fungal genera were composed of Piromyces, Neocallimastix, Anaeromyces, Aspergillus, and Orpinomyces. Principal coordinates analysis indicated that community structure of bacteria and fungi dynamically changed during 72 h fermentation. Bacterial networks with higher complexity had stronger stability than fungal networks. Most CAZyme families showed a significant decrease trend after 48 h fermentation. Functional genes related to hydrolysis decreased at 72 h, while functional genes involved in acidogenesis did not change significantly. These findings provide a in-depth understanding of mechanisms of lignocellulose degradation in Angus bull rumen, and may guide the construction and enrichment of rumen microorganisms in anaerobic fermentation of waste biomass.

Metagenomics analysis reveals differences in rumen microbiota in cows with low and high milk protein percentage

Variation exists in milk protein concentration of dairy cows of the same breed that are fed and managed in the same environment, and little information was available on this variation which might be attributed to differences in rumen microbial composition as well as their fermentation metabolites. This study is aimed at investigating the difference in the composition and functions of rumen microbiota as well as fermentation metabolites in Holstein cows with high and low milk protein concentrations. In this study, 20 lactating Holstein cows on the same diet were divided into two groups (10 cows each), high degree of milk protein group (HD), and low degree of milk protein (LD) concentrations based on previous milk composition history. Rumen content samples were obtained to explore the rumen fermentation parameters and rumen microbial composition. Shotgun metagenomics sequencing was employed to investigate the rumen microbial composition and sequences were assembled via the metagenomics binning technique. Metagenomics revealed that 6 Archaea genera, 5 Bacteria genera, 7 Eukaryota genera, and 7 virus genera differed significantly between the HD and LD group. The analysis of metagenome-assembled genomes (MAGs) showed that 2 genera (g__Eubacterium_H and g__Dialister) were significantly enriched (P < 0.05, linear discriminant analysis (LDA) > 2) in the HD group. However, the LD group recorded an increased abundance (P < 0.05, LDA > 2) of 8 genera (g__CAG-603, g__UBA2922, g__Ga6A1, g__RUG13091, g__Bradyrhizobium, g__Sediminibacterium, g__UBA6382, and g__Succinivibrio) when compared to the HD group. Furthermore, investigation of the KEGG genes revealed an upregulation in a higher number of genes associated with nitrogen metabolism and lysine biosynthesis pathways in the HD group as compared to the LD group. Therefore, the high milk protein concentration in the HD group could be explained by an increased ammonia synthesis by ruminal microbes which were converted to microbial amino acids and microbial protein (MCP) in presence of an increased energy source made possible by higher activities of carbohydrate-active enzymes (CAZymes). This MCP gets absorbed in the small intestine as amino acids and might be utilized for the synthesis of milk protein. KEY POINTS: • Rumen microbiota and their functions differed between cows with high milk protein % and those with low milk protein %. • The rumen microbiome of cows with high milk protein recorded a higher number of enriched genes linked to the nitrogen metabolism pathway and lysine biosynthesis pathway. • The activities of carbohydrate-active enzymes were found to be higher in the rumen of cows with high milk protein %.

Noise reduction strategies in metagenomic chromosome confirmation capture to link antibiotic resistance genes to microbial hosts

The gut microbiota is a reservoir for antimicrobial resistance genes (ARGs). With current sequencing methods, it is difficult to assign ARGs to their microbial hosts, particularly if these ARGs are located on plasmids. Metagenomic chromosome conformation capture approaches (meta3C and Hi-C) have recently been developed to link bacterial genes to phylogenetic markers, thus potentially allowing the assignment of ARGs to their hosts on a microbiome-wide scale. Here, we generated a meta3C dataset of a human stool sample and used previously published meta3C and Hi-C datasets to investigate bacterial hosts of ARGs in the human gut microbiome. Sequence reads mapping to repetitive elements were found to cause problematic noise in, and may importantly skew interpretation of, meta3C and Hi-C data. We provide a strategy to improve the signal-to-noise ratio by discarding reads that map to insertion sequence elements and to the end of contigs. We also show the importance of using spike-in controls to quantify whether the cross-linking step in meta3C and Hi-C protocols has been successful. After filtering to remove artefactual links, 87 ARGs were assigned to their bacterial hosts across all datasets, including 27 ARGs in the meta3C dataset we generated. We show that commensal gut bacteria are an important reservoir for ARGs, with genes coding for aminoglycoside and tetracycline resistance being widespread in anaerobic commensals of the human gut.

Characterization of metagenome-assembled genomes from the International Space Station

BACKGROUND

Several investigations on the microbial diversity and functional properties of the International Space Station (ISS) environment were carried out to understand the influence of spaceflight conditions on the microbial population. However, metagenome-assembled genomes (MAGs) of ISS samples are yet to be generated and subjected to various genomic analyses, including phylogenetic affiliation, predicted functional pathways, antimicrobial resistance, and virulence characteristics.

RESULTS

In total, 46 MAGs were assembled from 21 ISS environmental metagenomes, in which metaSPAdes yielded 20 MAGs and metaWRAP generated 26 MAGs. Among 46 MAGs retrieved, 18 bacterial species were identified, including one novel genus/species combination (Kalamiella piersonii) and one novel bacterial species (Methylobacterium ajmalii). In addition, four bins exhibited fungal genomes; this is the first-time fungal genomes were assembled from ISS metagenomes. Phylogenetic analyses of five bacterial species showed ISS-specific evolution. The genes pertaining to cell membranes, such as transmembrane transport, cell wall organization, and regulation of cell shape, were enriched. Variations in the antimicrobial-resistant (AMR) and virulence genes of the selected 20 MAGs were characterized to predict the ecology and evolution of biosafety level (BSL) 2 microorganisms in space. Since microbial virulence increases in microgravity, AMR gene sequences of MAGs were compared with genomes of respective ISS isolates and corresponding type strains. Among these 20 MAGs characterized, AMR genes were more prevalent in the Enterobacter bugandensis MAG, which has been predominantly isolated from clinical samples. MAGs were further used to analyze if genes involved in AMR and biofilm formation of viable microbes in ISS have variation due to generational evolution in microgravity and radiation pressure.

CONCLUSIONS

Comparative analyses of MAGs and whole-genome sequences of related ISS isolates and their type strains were characterized to understand the variation related to the microbial evolution under microgravity. The Pantoea/Kalamiella strains have the maximum single-nucleotide polymorphisms found within the ISS strains examined. This may suggest that Pantoea/Kalamiella strains are much more subjective to microgravity changes. The reconstructed genomes will enable researchers to study the evolution of genomes under microgravity and low-dose irradiation compared to the evolution of microbes here on Earth. Video Abstract.

Recovery of metagenomic data from the Aedes aegypti microbiome using a reproducible snakemake pipeline: MINUUR

Background: Ongoing research of the mosquito microbiome aims to uncover novel strategies to reduce pathogen transmission. Sequencing costs, especially for metagenomics, are however still significant. A resource that is increasingly used to gain insights into host-associated microbiomes is the large amount of publicly available genomic data based on whole organisms like mosquitoes, which includes sequencing reads of the host-associated microbes and provides the opportunity to gain additional value from these initially host-focused sequencing projects. Methods: To analyse non-host reads from existing genomic data, we developed a snakemake workflow called MINUUR (Microbial INsights Using Unmapped Reads). Within MINUUR, reads derived from the host-associated microbiome were extracted and characterised using taxonomic classifications and metagenome assembly followed by binning and quality assessment. We applied this pipeline to five publicly available Aedes aegypti genomic datasets, consisting of 62 samples with a broad range of sequencing depths. Results: We demonstrate that MINUUR recovers previously identified phyla and genera and is able to extract bacterial metagenome assembled genomes (MAGs) associated to the microbiome. Of these MAGS, 42 are high-quality representatives with >90% completeness and <5% contamination. These MAGs improve the genomic representation of the mosquito microbiome and can be used to facilitate genomic investigation of key genes of interest. Furthermore, we show that samples with a high number of KRAKEN2 assigned reads produce more MAGs. Conclusions: Our metagenomics workflow, MINUUR, was applied to a range of Aedes aegypti genomic samples to characterise microbiome-associated reads. We confirm the presence of key mosquito-associated symbionts that have previously been identified in other studies and recovered high-quality bacterial MAGs. In addition, MINUUR and its associated documentation are freely available on GitHub and provide researchers with a convenient workflow to investigate microbiome data included in the sequencing data for any applicable host genome of interest.