Functional gene profiling through metaRNAseq approach reveals diet-dependent variation in rumen microbiota of buffalo (Bubalus bubalis)

Integrated metabolome and microbiome analysis reveals the effect of rumen-protected sulfur-containing amino acids on the meat quality of Tibetan sheep meat

Introduction

This study investigated the effects of rumen-protected sulfur-containing amino acids (RPSAA) on the rumen and jejunal microbiota as well as on the metabolites and meat quality of the longissimus lumborum (LL) in Tibetan sheep.

Methods

By combining 16S rDNA sequencing with UHPLC-Q-TOF MS and Pearson correlation analysis, the relationship between gastrointestinal microbiota, muscle metabolites and meat quality was identified.

Results

The results showed that feeding RPSAA can increase the carcass weight, abdominal fat thickness (AP-2 group), and back fat thickness (AP-2 and AP-3 group) of Tibetan sheep. The water holding capacity (WHC), texture, and shear force (SF) of LL in the two groups also increased although the fatty acids content and brightness (L*) value significantly decreased in the AP-2 group. Metabolomics and correlation analysis further showed that RPSAA could significantly influence the metabolites in purine metabolism, thereby affecting L* and SF. In addition, RPSAA was beneficial for the fermentation of the rumen and jejunum. In both groups, the abundance of Prevotella 1, Lachnospiraceae NK3A20 group, Prevotella UCG-003, Lachnospiraceae ND3007 group in the rumen as well as the abundance of Eubacterium nodatum group and Mogibacterium group in the jejunum increased. In contrast, that of Turicibacter pathogens in the jejunum was reduced. The above microorganisms could regulate meat quality by regulating the metabolites (inosine, hypoxanthine, linoleic acid, palmitic acid, etc.) in purine and fatty acids metabolism.

Discussion

Overall, reducing the levels of crude proteins in the diet and feeding RPSAA is likely to improve the carcass quality of Tibetan sheep, with the addition of RPMET (AP-2) yielding the best edible quality, possibly due to its ability to influence the gastrointestinal microbiota to subsequently regulate muscle metabolites.

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.

Rumen microbial diversity, enteric methane emission and nutrient utilization of crossbred Karan-Fries cattle (Bos taurus) and Murrah buffalo (Bubalus bubalis) consuming varied roughage concentrate ratio

Dietary mix and host species have both been shown to have a significant impact on rumen microbial diversity, enteric methane emission and animal performance. The goal of this study was to see how the roughage concentrate ratio 70:30 (Low concentrate; LC) vs 40:60 (High concentrate; HC) and the host species crossbred cattle vs buffalo affected rumen microbial diversity, enteric methane emissions and nutrient utilization. Dry matter intake (kg/d) and dry matter percent digestibility were considerably (p < 0.05) higher in the HC ration and buffalo compared to LC ration and crossbred cattle, respectively. Both dietary mix and host species had a substantial (p < 0.05) impact on intake of various nutrients, including organic matter (OM), crude protein (CP), ether extract (EE), neutral detergent fiber (NDF), and acid detergent fiber (ADF). Increased concentrate proportion in the ration improved nitrogen balance, resulting in increased average daily gain and considerably reduced methane (g/d) output (p < 0.05). Furthermore, 16S rRNA genes were sequenced using Oxford Nanopore Technology (ONT) and subsequently annotated using the Centrifuge workflow to uncover ruminal bacterial diversity. Firmicutes was considerably (p < 0.01) greater in the LC diet, whereas, Bacteroidetes was higher in the HC ration. Genus Prevotella dominated all rumen samples, and buffalo fed LC ration had significantly (p < 0.01) higher Oscillospira abundance. At the species level, simple sugar-utilizing bacteria such as Prevotella spp. and Selenomonas ruminantium predominated in the crossbred cattle, but fibrolytic bacteria such as Oscillospira guilliermondii were statistically (p < 0.01) more abundant in the buffalo. Overall, dietary mix and host species have both been shown to have a significant impact on rumen microbial diversity, enteric methane emission and animal performance, however, host species remained a major driving force to change ruminal community composition as compared to roughage concentrate ratio under similar environmental conditions.

Characterizing rumen microbiota and CAZyme profile of Indian dromedary camel (Camelus dromedarius) in response to different roughages

In dromedary camels, which are pseudo-ruminants, rumen or C1 section of stomach is the main compartment involved in fiber degradation, as in true ruminants. However, as camels are adapted to the harsh and scarce grazing conditions of desert, their ruminal microbiota makes an interesting target of study. The present study was undertaken to generate the rumen microbial profile of Indian camel using 16S rRNA amplicon and shotgun metagenomics. The camels were fed three diets differing in the source of roughage. The comparative metagenomic analysis revealed greater proportions of significant differences between two fractions of rumen content followed by diet associated differences. Significant differences were also observed in the rumen microbiota collected at different time-points of the feeding trial. However, fraction related differences were more highlighted as compared to diet dependent changes in microbial profile from shotgun metagenomics data. Further, 16 genera were identified as part of the core rumen microbiome of Indian camels. Moreover, glycoside hydrolases were observed to be the most abundant among all Carbohydrate-Active enzymes and were dominated by GH2, GH3, GH13 and GH43. In all, this study describes the camel rumen microbiota under different dietary conditions with focus on taxonomic, functional, and Carbohydrate-Active enzymes profiles.

MetaRNAseq analysis of surti buffalo rumen content reveals that transcriptionally active microorganisms need not be abundant

The present study describes rumen microbiota composition and their functional profiles in Indian Surti buffaloes by metagenomic (MG) and metatranscriptomic (MT) approaches. The study compares samples from buffaloes fed three different proportion of roughages; green and dry type of roughage; and different rumen liquor fractions. Irrespective of sample, Bacteroidetes and Firmicutes were the most predominant bacterial phyla, followed by Proteobacteria, Fibrobacteres and Actinobacteria while, Prevotella, Bacteroides, Ruminococcus and Clostridium were the most abundant genera. Different proportions of taxa were observed in both MG and MT approaches indicating the differences in organisms present and organisms active in the rumen. Higher proportions of fungal taxa were observed in MT while important organisms like Fibrobacter and Butyrivibrio and abundant organisms like Bacteroides and Prevotella were underrepresented in MT data. Functionally, higher proportions of genes involved in Carbohydrate metabolism, Amino acid metabolism and Translation were observed in both data. Genes involved in Metabolism were observed to be underrepresented in MT data while, those involved in Genetic information processing were overrepresented in MT data. Further, genes involved in Carbohydrate metabolism were overexpressed compared to genes involved in Amino acid metabolism in MT data compared to MG data which had higher proportion of genes involved in Amino acid metabolism than Carbohydrate metabolism. In all significant differences were observed between both approaches, different fractions of rumen liquor (liquid and solid) and different proportions of roughage in diet.

Metatranscriptomics of the Hu sheep rumen microbiome reveals novel cellulases

BACKGROUND

Cellulosic biomass has great potential as a renewable biofuel resource. Robust, high-performance enzymes are needed to effectively utilize this valuable resource. In this study, metatranscriptomics was used to explore the carbohydrate-active enzymes (CAZymes), especially glycoside hydrolases (GHs), present in the rumen microbiome of Hu sheep. Select CAZymes were experimentally verified and characterized after cloning and expression in E. coli.

RESULTS

The metatranscriptomes of six Hu sheep rumen microbiomes yielded 42.3 Gbp of quality-checked sequence data that represented in total 2,380,783 unigenes after de novo assembling using Trinity and clustered with CD-HIT-EST. Annotation using the CAZy database revealed that 2.65% of the unigenes encoded GHs, which were assigned to 111 different CAZymes families. Firmicutes (18.7%) and Bacteroidetes (13.8%) were the major phyla to which the unigenes were taxonomically assigned. In total, 14,489 unigenes were annotated to 15 cellulase-containing GH families, with GH3, GH5 and GH9 being the predominant. From these putative cellulase-encoding unigenes, 4225 open reading frames (ORFs) were predicted to contain 2151 potential cellulase catalytic modules. Additionally, 147 ORFs were found to encode proteins that contain carbohydrate-binding modules (CBMs). Heterogeneous expression of 30 candidate cDNAs from the GH5 family in E. coli BL21 showed that 17 of the tested proteins had endoglucanase activity, while 7 exhibited exoglucanase activity. Interestingly, two of the GH5 proteins (Cel5A-h28 and Cel5A-h11) showed high specific activity against carboxymethylcellulose (CMC) and p-nitrophenyl-β-d-cellobioside (pNPC) (222.2 and 142.8 U/mg), respectively. The optimal pH value for activity of Cel5A-h11 and Cel5A-h28 was 6.0 for both enzymes, and optimal temperatures were 40 and 50 °C, respectively. Both enzymes retained over 70 and 60%, respectively, of their original activities after incubation at 40 °C for 60 min. However, their activities were rapidly diminished upon exposure to higher temperatures. Cel5A-h11 and Cel5A-h28 retained more than 80 and 60% of their maximal enzymatic activities after incubation for 16 h in buffered solutions in the pH range from 4.0 to 9.0.

CONCLUSION

The metatranscriptomic results revealed that the rumen microbiome of Hu sheep encoded a repertoire of new enzymes capable of cellulose degradation and metatranscriptomics was an effective method to discover novel cellulases for biotechnological applications.

STAble: a novel approach to de novo assembly of RNA-seq data and its application in a metabolic model network based metatranscriptomic workflow

Background

De novo assembly of RNA-seq data allows the study of transcriptome in absence of a reference genome either if data is obtained from a single organism or from a mixed sample as in metatranscriptomics studies. Given the high number of sequences obtained from NGS approaches, a critical step in any analysis workflow is the assembly of reads to reconstruct transcripts thus reducing the complexity of the analysis. Despite many available tools show a good sensitivity, there is a high percentage of false positives due to the high number of assemblies considered and it is likely that the high frequency of false positive is underestimated by currently used benchmarks. The reconstruction of not existing transcripts may false the biological interpretation of results as – for example – may overestimate the identification of “novel” transcripts. Moreover, benchmarks performed are usually based on RNA-seq data from annotated genomes and assembled transcripts are compared to annotations and genomes to identify putative good and wrong reconstructions, but these tests alone may lead to accept a particular type of false positive as true, as better described below.

Results

Here we present a novel methodology of de novo assembly, implemented in a software named STAble (Short-reads Transcriptome Assembler). The novel concept of this assembler is that the whole reads are used to determine possible alignments instead of using smaller k-mers, with the aim of reducing the number of chimeras produced. Furthermore, we applied a new set of benchmarks based on simulated data to better define the performance of assembly method and carefully identifying true reconstructions.

STAble was also used to build a prototype workflow to analyse metatranscriptomics data in connection to a steady state metabolic modelling algorithm. This algorithm was used to produce high quality metabolic interpretations of small gene expression sets obtained from already published RNA-seq data that we assembled with STAble.

Conclusions

The presented results, albeit preliminary, clearly suggest that with this approach is possible to identify informative reactions not directly revealed by raw transcriptomic data.

Microbiota composition, gene pool and its expression in Gir cattle (Bos indicus) rumen under different forage diets using metagenomic and metatranscriptomic approaches

Zebu (Bos indicus) is a domestic cattle species originating from the Indian subcontinent and now widely domesticated on several continents. In this study, we were particularly interested in understanding the functionally active rumen microbiota of an important Zebu breed, the Gir, under different dietary regimes. Metagenomic and metatranscriptomic data were compared at various taxonomic levels to elucidate the differential microbial population and its functional dynamics in Gir cattle rumen under different roughage dietary regimes. Different proportions of roughage rather than the type of roughage (dry or green) modulated microbiome composition and the expression of its gene pool. Fibre degrading bacteria (i.e. Clostridium, Ruminococcus, Eubacterium, Butyrivibrio, Bacillus and Roseburia) were higher in the solid fraction of rumen (P<0.01) compared to the liquid fraction, whereas bacteria considered to be utilizers of the degraded product (i.e. Prevotella, Bacteroides, Parabacteroides, Paludibacter and Victivallis) were dominant in the liquid fraction (P<0.05). Likewise, expression of fibre degrading enzymes and related carbohydrate binding modules (CBMs) occurred in the solid fraction. When metagenomic and metatranscriptomic data were compared, it was found that some genera and species were transcriptionally more active, although they were in low abundance, making an important contribution to fibre degradation and its further metabolism in the rumen. This study also identified some of the transcriptionally active genera, such as Caldicellulosiruptor and Paludibacter, whose potential has been less-explored in rumen. Overall, the comparison of metagenomic shotgun and metatranscriptomic sequencing appeared to be a much richer source of information compared to conventional metagenomic analysis.

PCR and Omics Based Techniques to Study the Diversity, Ecology and Biology of Anaerobic Fungi: Insights, Challenges and Opportunities

Anaerobic fungi (phylum Neocallimastigomycota) are common inhabitants of the digestive tract of mammalian herbivores, and in the rumen, can account for up to 20% of the microbial biomass. Anaerobic fungi play a primary role in the degradation of lignocellulosic plant material. They also have a syntrophic interaction with methanogenic archaea, which increases their fiber degradation activity. To date, nine anaerobic fungal genera have been described, with further novel taxonomic groupings known to exist based on culture-independent molecular surveys. However, the true extent of their diversity may be even more extensively underestimated as anaerobic fungi continue being discovered in yet unexplored gut and non-gut environments. Additionally many studies are now known to have used primers that provide incomplete coverage of the Neocallimastigomycota. For ecological studies the internal transcribed spacer 1 region (ITS1) has been the taxonomic marker of choice, but due to various limitations the large subunit rRNA (LSU) is now being increasingly used. How the continued expansion of our knowledge regarding anaerobic fungal diversity will impact on our understanding of their biology and ecological role remains unclear; particularly as it is becoming apparent that anaerobic fungi display niche differentiation. As a consequence, there is a need to move beyond the broad generalization of anaerobic fungi as fiber-degraders, and explore the fundamental differences that underpin their ability to exist in distinct ecological niches. Application of genomics, transcriptomics, proteomics and metabolomics to their study in pure/mixed cultures and environmental samples will be invaluable in this process. To date the genomes and transcriptomes of several characterized anaerobic fungal isolates have been successfully generated. In contrast, the application of proteomics and metabolomics to anaerobic fungal analysis is still in its infancy. A central problem for all analyses, however, is the limited functional annotation of anaerobic fungal sequence data. There is therefore an urgent need to expand information held within publicly available reference databases. Once this challenge is overcome, along with improved sample collection and extraction, the application of these techniques will be key in furthering our understanding of the ecological role and impact of anaerobic fungi in the wide range of environments they inhabit.

Unraveling the Microbial Interactions and Metabolic Potentials in Pre- and Post-treated Sludge from a Wastewater Treatment Plant Using Metagenomic Studies

Sewage waste represents an ecosystem of complex and interactive microbial consortia which proliferate with different kinetics according to their individual genetic as well as metabolic potential. We performed metagenomic shotgun sequencing on Ion-Torrent platform, to explore the microbial community structure, their biological interactions and associated functional capacity of pre-treated/raw sludge (RS) and post-treated/dried sludge (DS) of wastewater treatment plant. Bacterial phylotypes belonging to Epsilonproteobacteria (∼45.80%) dominated the RS with relatively few Archaea (∼1.94%) whereas DS has the dominance of beta- (30.23%) and delta- (13.38%) classes of Proteobacteria with relatively greater abundance of Archaea (∼7.18%). In particular, Epsilonproteobacteria appears as a primary energy source in RS and sulfur-reducing bacteria with methanogens seems to be in the potential syntrophic association in DS. These interactions could be ultimately responsible for carrying out amino-acid degradation, aromatic compound degradation and degradation of propionate and butyrate in DS. Our data also reveal the presence of key genes in the sludge microbial community responsible for degradation of polycyclic aromatic hydrocarbons. Potential pathogenic microbes and genes for the virulence factors were found to be relatively abundant in RS which clearly reflect the necessity of treatment of RS. After treatment, potential pathogens load was reduced, indicating the sludge hygienisation in DS. Additionally, the interactions found in this study would reveal the biological and environmental cooperation among microbial communities for domestic wastewater treatment.