Diversity and functional analysis of rumen and fecal microbial communities associated with dietary changes in crossbreed dairy cattle

The objective of this study was to investigate the effect of varying roughage and concentrate proportions, in diet of crossbreed dairy cattle, on the composition and associated functional genes of rumen and fecal microbiota. We also explored fecal samples as a proxy for rumen liquor samples. Six crossbred dairy cattle were reared on three diets with an increasing concentrate and reducing roughage amount in three consecutive 10-day periods. After each period, individual rumen liquor and fecal samples were collected and analyzed through shotgun metagenomic sequencing. Average relative abundance of identified Operational Taxonomic Units (OTU) and microbial functional roles from all animals were compared between diets and sample types (fecal and rumen liquor). Results indicated that dietary modifications significantly affected several rumen and fecal microbial OTUs. In the rumen, an increase in dietary concentrate resulted in an upsurge in the abundance of Proteobacteria, while reducing the proportions of Bacteroidetes and Firmicutes. Conversely, changes in microbial composition in fecal samples were not consistent with dietary modification patterns. Microbial functional pathway classification identified that carbohydrate metabolism and protein metabolism pathways dominated microbial roles. Assessment of dietary effects on the predicted functional roles of these microbiota revealed that a high amount of dietary concentrate resulted in an increase in central carbohydrate metabolism and a corresponding reduction in protein synthesis. Moreover, we identified several microbial stress-related responses linked to dietary changes. Bacteroides and Clostridium genera were the principal hosts of these microbial functions. Therefore, the roughage to concentrate proportion has more influence on the microbial composition and microbial functional genes in rumen samples than fecal samples. As such, we did not establish a significant relationship between the rumen and fecal metagenome profiles, and the rumen and fecal microbiota from one animal did not correlate more than those from different animals.

A Comparison of Two DNA Metagenomic Bioinformatic Pipelines While Evaluating the Microbial Diversity in Feces of Tanzanian Small Holder Dairy Cattle

Analysis of shotgun metagenomic data generated from next generation sequencing platforms can be done through a variety of bioinformatic pipelines. These pipelines employ different sets of sophisticated bioinformatics algorithms which may affect the results of this analysis. In this study, we compared two commonly used pipelines for shotgun metagenomic analysis: MG-RAST and Kraken 2, in terms of taxonomic classification, diversity analysis, and usability using their primarily default parameters. Overall, the two pipelines detected similar abundance distributions in the three most abundant taxa Proteobacteria, Firmicutes, and Bacteroidetes. Within bacterial domain, 497 genera were identified by both pipelines, while an additional 694 and 98 genera were solely identified by Kraken 2 and MG-RAST, respectively. 933 species were detected by the two algorithms. Kraken 2 solely detected 3550 species, while MG-RAST identified 557 species uniquely. For archaea, Kraken 2 generated 105 and 236 genera and species, respectively, while MG-RAST detected 60 genera and 88 species. 54 genera and 72 species were commonly detected by the two methods. Kraken 2 had a quicker analysis time (~4 hours) while MG-RAST took approximately 2 days per sample. This study revealed that Kraken 2 and MG-RAST generate comparable results and that a reliable high-level overview of sample is generated irrespective of the pipeline selected. However, Kraken 2 generated a more accurate taxonomic identification given the higher number of "Unclassified" reads in MG-RAST. The observed variations at the genus level show that a main restriction is using different databases for classification of the metagenomic data. The results of this research indicate that a more inclusive and representative classification of microbiomes may be achieved through creation of the combined pipelines.

RNA-seq profiling of skin in temperate and tropical cattle

Skin is a major thermoregulatory organ in the body controlling homeothermy, a critical function for climate adaptation. We compared genes expressed between tropical- and temperate-adapted cattle to better understand genes involved in climate adaptation and hence thermoregulation. We profiled the skin of representative tropical and temperate cattle using RNA-seq. A total of 214,754,759 reads were generated and assembled into 72,993,478 reads and were mapped to unique regions in the bovine genome. Gene coverage of unique regions of the reference genome showed that of 24,616 genes, only 13,130 genes (53.34%) displayed more than one count per million reads for at least two libraries and were considered suitable for downstream analyses. Our results revealed that of 255 genes expressed differentially, 98 genes were upregulated in tropically-adapted White Fulani (WF; Bos indicus) and 157 genes were down regulated in WF compared to Angus, AG (Bos taurus). Fifteen pathways were identified from the differential gene sets through gene ontology and pathway analyses. These include the significantly enriched melanin metabolic process, proteinaceous extracellular matrix, inflammatory response, defense response, calcium ion binding and response to wounding. Quantitative PCR was used to validate six representative genes which are associated with skin thermoregulation and epithelia dysfunction (mean correlation 0.92; p < 0.001). Our results contribute to identifying genes and understanding molecular mechanisms of skin thermoregulation that may influence strategic genomic selection in cattle to withstand climate adaptation, microbial invasion and mechanical damage.