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

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.

The multi-kingdom microbiome catalog of the chicken gastrointestinal tract

Chicken is an important food animal worldwide and plays an important role in human life by providing meat and eggs. Despite recent significant advances in gut microbiome studies, a comprehensive study of chicken gut bacterial, archaeal, and viral genomes remains unavailable. In this study, we constructed a chicken multi-kingdom microbiome catalog (CMKMC), including 18,201 bacterial, 225 archaeal, and 33,411 viral genomes, and annotated over 6,076,006 protein-coding genes by integrating 135 chicken gut metagenomes and publicly available metagenome-assembled genomes (MAGs) from ten countries. We found that 812 and 240 MAGs in our dataset were putative novel species and genera, respectively, far beyond what was previously reported. The newly unclassified MAGs were predominant in Phyla Firmicutes_A (n = 263), followed by Firmicutes (n = 126), Bacteroidota (n = 121), and Proteobacteria (n = 87). Most of the classified species-level viral operational taxonomic units belong to Caudovirales. Approximately, 63.24 % of chicken gut viromes are predicted to infect two or more hosts, including complete circular viruses. Moreover, we found that diverse auxiliary metabolic genes and antibiotic resistance genes were carried by viruses. Together, our CMKMC provides the largest integrated MAGs and viral genomes from the chicken gut to date, functional insights into the chicken gastrointestinal tract microbiota, and paves the way for microbial interventions for better chicken health and productivity.

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.

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.

Ensiled diet improved the growth performance of Tibetan sheep by regulating the rumen microbial community and rumen epithelial morphology

The aim of this study was to investigate the effects of ensiled agricultural byproducts from Qinghai-Tibet plateau on growth performance, rumen microbiota, ruminal epithelium morphology, and nutrient transport-related gene expression in Tibetan sheep. Fourteen male Tibetan sheep were randomly assigned to one of two diets: an untreated diet (without silage inoculum, CON, n = 7) or an ensiled diet (with silage inoculum, ESD, n = 7). The total experimental period lasted for 84 d, including early 14 d as adaption period and remaining 70 d for data collection. The ESD increased average daily gain (P = 0.046), dry matter intake (P < 0.001), ammonia nitrogen (P = 0.045), microbial crude protein (P = 0.034), and total volatile fatty acids concentration (P < 0.001), and decreased ruminal pH value (P = 0.014). The proportion of propionate (P = 0.006) and the copy numbers of bacteria (P = 0.01) and protozoa (P = 0.002) were higher, while the proportion of acetate (P = 0.028) was lower in the sheep fed ESD compared to CON. Pyrosequencing of the 16S ribosomal RNA gene revealed that ESD increased the relative abundance of Firmicutes, Ruminococcus, Lachnospiraceae_AC2044_group, Lachnospiraceae_XPB1014_group, and Christensenellaceae_R-7_group in the rumen (P < 0.05), while decreased the relative abundance of Bacteroidota, Prevotellaceae_UCG-003, and Veillonellaceae_UCG-001 (P < 0.05). Analyses with PICRUSt2 and STAMP indicated that the propionate metabolism pathway was enriched in the sheep fed ESD (P = 0.026). The ESD increased the rumen papillae height (P = 0.012), density (P = 0.036), and surface area (P = 0.001), and improved the thickness of the total epithelia (P = 0.018), stratum corneum (P = 0.040), stratum granulosum (P = 0.042), and stratum spinosum and basale (P = 0.004). The relative mRNA expression of cyclin-dependent Kinase 2, CyclinA2, CyclinD2, zonula occludens-1, Occludin, monocarboxylate transporter isoform 1 (MCT1), MCT4, sodium/potassium pump, and sodium/hydrogen antiporter 3 were higher in the rumen epithelial of sheep fed ESD than CON (P < 0.05). Conversely, the relative mRNA expressions of Caspase 3 and B-cell lymphoma-2 were lower in the sheep fed ESD than CON (P < 0.05). In conclusion, compared with an untreated diet, feeding an ensiled diet altered the rumen microbial community, enhanced nutrient transport through rumen epithelium, and improved the growth performance of Tibetan sheep.

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.

Gut microbiome signatures of extreme environment adaption in Tibetan pig

Tibetan pigs (TPs) can adapt to the extreme environments in the Tibetan plateau implicated by their self-genome signals, but little is known about roles of the gut microbiota in the host adaption. Here, we reconstructed 8210 metagenome-assembled genomes from TPs (n = 65) living in high-altitude and low-altitude captive pigs (87 from China-CPs and 200 from Europe-EPs) that were clustered into 1050 species-level genome bins (SGBs) at the threshold of 95% average nucleotide identity. 73.47% of SGBs represented new species. The gut microbial community structure analysis based on 1,048 SGBs showed that TPs was significantly different from low-altitude captive pigs. TP-associated SGBs enabled to digest multiple complex polysaccharides, including cellulose, hemicellulose, chitin and pectin. Especially, we found TPs showed the most common enrichment of phyla Fibrobacterota and Elusimicrobia, which were involved in the productions of short- and medium-chain fatty acids (acetic acid, butanoate and propanoate; octanomic, decanoic and dodecanoic acids), as well as in the biosynthesis of lactate, 20 essential amino acids, multiple B vitamins (B1, B2, B3, B5, B7 and B9) and cofactors. Unexpectedly, Fibrobacterota solely showed powerful metabolic capacity, including the synthesis of acetic acid, alanine, histidine, arginine, tryptophan, serine, threonine, valine, B2, B5, B9, heme and tetrahydrofolate. These metabolites might contribute to host adaptation to high-altitude, such as energy harvesting and resistance against hypoxia and ultraviolet radiation. This study provides insights into understanding the role of gut microbiome played in mammalian high-altitude adaptation and discovers some potential microbes as probiotics for improving animal health.

Catch me if you can: capturing microbial community transformation by extracellular DNA using Hi-C sequencing

The transformation of environmental microorganisms by extracellular DNA is an overlooked mechanism of horizontal gene transfer and evolution. It initiates the acquisition of exogenous genes and propagates antimicrobial resistance alongside vertical and conjugative transfers. We combined mixed-culture biotechnology and Hi-C sequencing to elucidate the transformation of wastewater microorganisms with a synthetic plasmid encoding GFP and kanamycin resistance genes, in the mixed culture of chemostats exposed to kanamycin at concentrations representing wastewater, gut and polluted environments (0.01-2.5-50-100 mg L^-1). We found that the phylogenetically distant Gram-negative Runella (102 Hi-C links), Bosea (35), Gemmobacter (33) and Zoogloea (24) spp., and Gram-positive Microbacterium sp. (90) were transformed by the foreign plasmid, under high antibiotic exposure (50 mg L^-1). In addition, the antibiotic pressure shifted the origin of aminoglycoside resistance genes from genomic DNA to mobile genetic elements on plasmids accumulating in microorganisms. These results reveal the power of Hi-C sequencing to catch and surveil the transfer of xenogenetic elements inside microbiomes.

Lignocellulose degradation by rumen bacterial communities: New insights from metagenome analyses

Ruminant animals house a dense and diverse community of microorganisms in their rumen, an enlarged compartment in their stomach, which provides a supportive environment for the storage and microbial fermentation of ingested feeds dominated by plant materials. The rumen microbiota has acquired diverse and functionally overlapped enzymes for the degradation of plant cell wall polysaccharides. In rumen Bacteroidetes, enzymes involved in degradation are clustered into polysaccharide utilization loci to facilitate coordinated expression when target polysaccharides are available. Firmicutes use free enzymes and cellulosomes to degrade the polysaccharides. Fibrobacters either aggregate lignocellulose-degrading enzymes on their cell surface or release them into the extracellular medium in membrane vesicles, a mechanism that has proven extremely effective in the breakdown of recalcitrant cellulose. Based on current metagenomic analyses, rumen Bacteroidetes and Firmicutes are categorized as generalist microbes that can degrade a wide range of polysaccharides, while other members adapted toward specific polysaccharides. Particularly, there is ample evidence that Verrucomicrobia and Spirochaetes have evolved enzyme systems for the breakdown of complex polysaccharides such as xyloglucans, peptidoglycans, and pectin. It is concluded that diversity in degradation mechanisms is required to ensure that every component in feeds is efficiently degraded, which is key to harvesting maximum energy by host animals.

Strain-Level Dynamics Reveal Regulatory Roles in Atopic Eczema by Gut Bacterial Phages

The vast population of bacterial phages or viruses (virome) plays pivotal roles in the ecology of human microbial flora and health conditions. Obstacles, including poor viral sequence inference, strain-sensitive virus-host relationship, and the high diversity among individuals, hinder the in-depth understanding of the human virome. We conducted longitudinal studies of the virome based on constructing a high-quality personal reference metagenome (PRM). By applying long-read sequencing for representative samples, we could build a PRM of high continuity that allows accurate annotation and abundance estimation of viruses and bacterial species in all samples of the same individual by aligning short sequencing reads to the PRM. We applied this approach to a series of fecal samples collected for 6 months from a 2-year-old boy who had experienced a 2-month flare-up of atopic eczema (dermatitis) in this period. We identified 31 viral strains in the patient's gut microbiota and deciphered their strain-level relationship to their bacterial hosts. Among them, a lytic crAssphage developed into a dozen substrains and coordinated downregulation in the catabolism of aromatic amino acids (AAAs) in their host bacteria which govern the production of immune-active AAA derivates. The metabolic alterations confirmed based on metabolomic assays cooccurred with symptom remission. Our PRM-based analysis provides an easy approach for deciphering the dynamics of the strain-level human gut virome in the context of entire microbiota. Close temporal correlations among virome alteration, microbial metabolism, and disease remission suggest a potential mechanism for how bacterial phages in microbiota are intimately related to human health. IMPORTANCE The vast populations of viruses or bacteriophages in human gut flora remain mysterious. However, poor annotation and abundance estimation remain obstacles to strain-level analysis and clarification of their roles in microbiome ecology and metabolism associated with human health and diseases. We demonstrate that a personal reference metagenome (PRM)-based approach provides strain-level resolution for analyzing the gut microbiota-associated virome. When applying such an approach to longitudinal samples collected from a 2-year-old boy who has experienced a 2-month flare-up of atopic eczema, we observed thriving substrains of a lytic crAssphage, showing temporal correlation with downregulated catabolism of aromatic amino acids, lower production of immune-active metabolites, and remission of the disease. The PRM-based approach is practical and powerful for strain-centric analysis of the human gut virome, and the underlying mechanism of how strain-level virome dynamics affect disease deserves further investigation.

Nutrient availability of roughages in isocaloric and isonitrogenous diets alters the bacterial networks in the whole gastrointestinal tract of Hu sheep

BACKGROUND

The nutrient availability of roughages could affect the dietary utilization efficiency of ruminants even in isocaloric and isonitrogenous diets. Here, we analyzed the bacterial composition and their metabolic pathways in the gastrointestinal tracts (GITs) of Hu sheep fed with wheat straw (WS) instead of alfalfa (AL) in isocaloric and isonitrogenous diets, trying to explore the reasons from the perspective of GITs bacterial network structure changes.

RESULTS

We employed 16S rRNA gene sequencing in combination with the Kruskal-Wallis test, Spearman correlation analysis, and other statistical methods to describe the microbiota composition in the GITs of Hu sheep. The results showed after the roughage was replaced from AL to WS, the most positive response occurred in the rumen microbiota, resulting in a more obvious microbiological and functional redundancy phenomenon. Whereas extended biogeographic studies of the GITs bacterial community found opposite results for the hindgut microbiota and metabolism networks compared to the forestomach. The abundance of fiber-degrading bacteria such as Prevotella, Oscillospiraceae NK4A214 group, and Treponema was significantly increased in GITs, but low-efficiency crude fiber degradation inhibited energy use efficiency, the pentose phosphate pathway, gluconeogenesis, and volatile acid synthesis. In addition, dietary shifting from AL to WS decreased the abundance of beneficial bacteria such as the Lachnospiraceae NK3A20 group and Alistipes, thereby enhancing the underlying inflammatory response.

CONCLUSIONS

These findings suggest that feeding untreated WS affected the structure and function of the bacterial network in the GITs due to limited total digestible nutrients, and in particular increases the complexity of the rumen bacterial network, and limit the abundance of bacteria involved in the crude fiber degradation in the hindgut.