Comparative metagenomic and metatranscriptomic analyses reveal the breed effect on the rumen microbiome and its associations with feed efficiency in beef cattle

Abstract Background Microorganisms are responsible for fermentation within the rumen and have been reported to contribute to the variation in feed efficiency of cattle. However, to what extent the breed affects the rumen microbiome and its association with host feed efficiency is unknown. Here, rumen microbiomes of beef cattle (n = 48) from three breeds (Angus, Charolais, Kinsella composite hybrid) with high and low feed efficiency were explored using metagenomics and metatranscriptomics, aiming to identify differences between functional potentials and activities of same rumen microbiomes and to evaluate the effects of host breed and feed efficiency on the rumen microbiome. Results Rumen metagenomes were more closely clustered together and thus more conserved among individuals than metatranscriptomes, suggesting that inter-individual functional variations at the RNA level were higher than those at the DNA level. However, while mRNA enrichment significantly increased the sequencing depth of mRNA and generated similar functional profiles to total RNA-based metatranscriptomics, it led to biased abundance estimation of several transcripts. We observed divergent rumen microbial composition (metatranscriptomic level) and functional potentials (metagenomic level) among three breeds, but differences in functional activity (metatranscriptomic level) were less apparent. Differential rumen microbial features (e.g., taxa, diversity indices, functional categories, and genes) were detected between cattle with high and low feed efficiency, and most of them were breed-specific. Conclusions Metatranscriptomes represent real-time functional activities of microbiomes and have the potential to better associate rumen microorganisms with host performances compared to metagenomics. As total RNA-based metatranscriptomics seem to avoid potential biases caused by mRNA enrichment and allow simultaneous use of rRNA for generation of compositional profiles, we suggest their use for linking the rumen microbiome with host phenotypes in future studies. However, if exploration of specific lowly expressed genes is desired, mRNA enrichment is recommended as it will enhance the resolution of mRNA. Finally, the differential microbial features observed between efficient and inefficient steers tended to be specific to breeds, suggesting that interactions between host breed genotype and the rumen microbiome contribute to the variations in feed efficiency observed. These breed-associated differences represent an opportunity to engineer specific rumen microbiomes through selective breeding of the hosts

Characterization of bovine miRNAs by sequencing and bioinformatics analysis

<p>Abstract</p> <p>Background</p> <p>MicroRNAs (miRNAs) are a family of ~22 nucleotide small RNA molecules which regulate gene expression by fully or partially binding to their complementary sequences in mRNAs or promoters. A large number of miRNAs and their expression patterns have been reported in human, mouse and rat. However, miRNAs and their expression patterns in live stock species such as beef cattle are not well studied.</p> <p>Results</p> <p>We constructed and sequenced small-RNA libraries to yield a total of 13,541 small-RNA sequences from 11 bovine tissues including brain, subcutaneous fat, muscle, liver, kidney, spleen and thymus. In total, 228 miRNAs including 29 novel miRNA candidates were identified. Of the 199 miRNAs, 101 have been previously reported as bovine miRNAs and the other 98 are bovine orthologs of known miRNAs that have been identified in at least one other mammalian species. Of the 29 novel miRNA candidates, 17 appeared at this point in time to be bovine specific, while the remaining 12 had evidence of evolutionary conservation in other mammalian species. Five miRNAs (miR-23a, -23b, -99a, -125b and -126-5p) were very abundant across the 11 tissues, accounting for 44.3% of all small RNA sequences. The expression analysis of selected miRNAs using qRT-PCR also showed that miR-26a and -99a were highly expressed in all tissues, while miR-122 and miR-133a were predominantly expressed in liver and muscle, respectively.</p> <p>Conclusion</p> <p>The miRNA expression patterns among 11 tissues from beef cattle revealed that most miRNAs were ubiquitously expressed in all tissues, while only a few miRNAs were tissue specific. Only 60% miRNAs in this study were found to display strand bias, suggesting that there are some key factors for mature miRNA selection other than internal stability. Most bovine miRNAs are highly conserved in other three mammalian species, indicating that these miRNAs may have a role in different species that are potential molecular markers for evolution.</p

Characterization of microRNA expression in bovine adipose tissues: a potential regulatory mechanism of subcutaneous adipose tissue development

<p>Abstract</p> <p>Background</p> <p>MicroRNAs (miRNAs), a family of small non-coding RNA molecules, appear to regulate animal lipid metabolism and preadipocyte conversion to form lipid-assimilating adipocytes (<it>i.e. </it>adipogenesis). However, no miRNA to date has been reported to modulate adipogenesis and lipid deposition in beef cattle.</p> <p>Results</p> <p>The expression patterns of 89 miRNAs including four bovine specific miRNAs in subcutaneous adipose tissues from three groups of crossbred steers differing in backfat thickness were compared using qRT-PCR analysis. Eighty-six miRNAs were detectable in all samples, with 42 miRNAs differing among crossbreds (P < 0.05) and 15 miRNAs differentially expressed between tissues with high and low backfat thickness (P < 0.05). The expression levels of 18 miRNAs were correlated with backfat thickness (P < 0.05). The miRNA most differentially expressed and the most strongly associated with backfat thickness was miR-378, with a 1.99-fold increase in high backfat thickness tissues (r = 0.72).</p> <p>Conclusions</p> <p>MiRNA expression patterns differed significantly in response to host genetic components. Approximately 20% of the miRNAs in this study were identified as being correlated with backfat thickness. This result suggests that miRNAs may play a regulatory role in white adipose tissue development in beef animals.</p

Proteome differences associated with fat accumulation in bovine subcutaneous adipose tissues

<p>Abstract</p> <p>Background</p> <p>The fat components of red meat products have been of interest to researchers due to the health aspects of excess fat consumption by humans. We hypothesized that differences in protein expression have an impact on adipose tissue formation during beef cattle development and growth. Therefore, in this study we evaluated the differences in the discernable proteome of subcutaneous adipose tissues of 35 beef crossbred steers [Charolais × Red Angus (CHAR) (n = 13) and Hereford × Angus (HEAN) (n = 22)] with different back fat (BF) thicknesses. The goal was to identify specific protein markers that could be associated with adipose tissue formation in beef cows.</p> <p>Results</p> <p>Approximately 541-580 protein spots were detected and compared in each crossbred group, and 33 and 36 protein spots showed expression differences between tissues with high and low BF thicknesses from HEAN and CHAR crossbed, respectively. The annexin 1 protein was highly expressed in both crossbred steers that had a higher BF thickness (p < 0.05) and this was further validated by a western blot analysis. In 13 tissues of CHAR animals and 22 tissues of HEAN animals, the relative expression of annexin 1 was significantly different (p < 0.05) between tissues with high and low BF thicknesses.</p> <p>Conclusion</p> <p>The increased expression of annexin 1 protein has been found to be associated with higher BF thickness in both crossbred steers. This result lays the foundation for future studies to develop the protein marker for assessing animals with different BF thickness.</p

Comparative transcriptomic analysis of rectal tissue from beef steers revealed reduced host immunity in Escherichia coli 0157:H7 super-shedders

Sherpa Romeo green journal. Open access, distributed under the terms of the Creative Commons Attribution LicenseSuper-shedder cattle are a major disseminator of E . coli O157:H7 into the environment, and the terminal rectum has been proposed as the primary E . coli O157:H7 colonization site. This study aimed to identify host factors that are associated with the super-shedding pro- cess by comparing transcriptomic profiles in rectal tissue collected from 5 super-shedder cattle and 4 non-shedder cattle using RNA-Seq. In total, 17,859 ± 354 genes and 399 ± 16 miRNAs were detected, and 11,773 genes were expressed in all animals. Fifty-eight differ- entially expressed (DE) genes (false discovery rate < 0.05) including 11 up-regulated and 47 down-regulated (log 2 (fold change) ranged from -5.5 to 4.2), and 2 up-regulated DE miRNAs (log 2 (fold change) = 2.1 and 2.5, respectively) were identified in super-shedders compared to non-shedders. Functional analysis of DE genes revealed that 31 down-regu- lated genes were potentially associated with reduced innate and adaptive immune functions in super-shedders, including 13 lymphocytes membrane receptors, 3 transcription factors and 5 cytokines, suggesting the decreased key host immune functions in the rectal tissue of super-shedders, including decreased quantity and migration of immune cells such as lym- phocytes, neutrophils and dendritic cells. The up-regulation of bta-miR-29d-3p and the down regulation of its predicted target gene, regulator of G-protein signaling 13 , suggested a potential regulatory role of this miRNA in decreased migration of lymphocytes in super- shedders. Based on these findings, the rectal tissue of super-shedders may inherently exhibit less effective innate and adaptive immune protection. Further study is required to confirm if such effect on host immunity is due to the nature of the host itself or due to actions mediated by E . coli O157:H7.Ye

Effect of Medicinal Plant By-products Supplementation to Total Mixed Ration on Growth Performance, Carcass Characteristics and Economic Efficacy in the Late Fattening Period of Hanwoo Steers

This study was conducted to evaluate the effect of medicinal plant by-products (MPB) supplementation to a total mixed ration (TMR) on growth, carcass characteristics and economic efficacy in the late fattening period of Hanwoo steers. Twenty seven steers (body weight [BW], 573±57 kg) were assigned to 3 treatment groups so that each treatment based on BW contained 9 animals. All groups received ad libitum TMR throughout the feeding trial until slaughter (from 24 to 30 months of age) and treatments were as follows: control, 1,000 g/kg TMR; treatment 1 (T1), 970 g/kg TMR and 30 g/kg MPB; treatment 2 (T2), 950 g/kg TMR and 50 g/kg MPB. Initial and final BW were not different among treatments. Resultant data were analyzed using general linear models of SAS. Average daily gain and feed efficiency were higher (p<0.05) for T1 than control, but there was no difference between control and T2. Plasma albumin showed low-, intermediate- and high-level (p<0.05) for control, T1 and T2, whereas non-esterified fatty acid was high-, intermediate- and high-level (p<0.05) for control, T1 and T2, respectively. Carcass weight, carcass rate, backfat thickness and rib eye muscle area were not affected by MPB supplementation, whereas quality and yield grades were highest (p<0.05) for T1 and T2, respectively. Daily feed costs were decreased by 0.5% and 0.8% and carcass prices were increased by 18.1% and 7.6% for T1 and T2 compared to control, resulting from substituting TMR with 30 and 50 g/kg MPB, respectively. In conclusion, the substituting TMR by 30 g/kg MPB may be a potential feed supplement approach to improve economic efficacy in the late fattening period of Hanwoo steers

Saccharomyces cerevisiae boulardii accelerates intestinal microbiota maturation and is correlated with increased secretory IgA production in neonatal dairy calves

Neonatal calves have a limited capacity to initiate immune responses due to a relatively immature adaptive immune system, which renders them susceptible to many on-farm diseases. At birth, the mucosal surfaces of the intestine are rapidly colonized by microbes in a process that promotes mucosal immunity and primes the development of the adaptive immune system. In a companion study, our group demonstrated that supplementation of a live yeast probiotic, Saccharomyces cerevisiae boulardii (SCB) CNCM I-1079, to calves from birth to 1 week of age stimulates secretory IgA (sIgA) production in the intestine. The objective of the study was to evaluate how SCB supplementation impacts the intestinal microbiota of one-week-old male calves, and how changes in the bacterial community in the intestine relate to the increase in secretory IgA. A total of 20 calves were randomly allocated to one of two treatments at birth: Control (CON, n = 10) fed at 5 g/d of carrier with no live yeast; and SCB (n = 10) fed at 5 g of live SCB per day (10 × 109 CFU/d). Our study revealed that supplementing calves with SCB from birth to 1 week of age had its most marked effects in the ileum, increasing species richness and phylogenetic diversity in addition to expediting the transition to a more interconnected bacterial community. Furthermore, LEfSe analysis revealed that there were several differentially abundant taxa between treatments and that SCB increased the relative abundance the family Eubacteriaceae, Corynebacteriaceae, Eggerthellaceae, Bacillaceae, and Ruminococcaceae. Furthermore, network analysis suggests that SCB promoted a more stable bacterial community and appears to reduce colonization with Shigella. Lastly, we observed that the probiotic-driven increase in microbial diversity was highly correlated with the enhanced secretory IgA capacity of the ileum, suggesting that the calf’s gut mucosal immune system relies on the development of a stable and highly diverse microbial community to provide the necessary cues to train and promote its proper function. In summary, this data shows that supplementation of SCB promoted establishment of a diverse and interconnected microbiota, prevented colonization of Escherichia Shigella and indicates a possible role in stimulating humoral mucosal immunity

Genetic architecture of gene expression underlying variation in host response to porcine reproductive and respiratory syndrome virus infection

It has been shown that inter-individual variation in host response to porcine reproductive and respiratory syndrome (PRRS) has a heritable component, yet little is known about the underlying genetic architecture of gene expression in response to PRRS virus (PRRSV) infection. Here, we integrated genome-wide genotype, gene expression, viremia level, and weight gain data to identify genetic polymorphisms that are associated with variation in inter-individual gene expression and response to PRRSV infection in pigs. RNA-seq analysis of peripheral blood samples collected just prior to experimental challenge (day 0) and at 4, 7, 11 and 14 days post infection from 44 pigs revealed 6,430 differentially expressed genes at one or more time points post infection compared to the day 0 baseline. We mapped genetic polymorphisms that were associated with inter-individual differences in expression at each day and found evidence of cis-acting expression quantitative trait loci (cis-eQTL) for 869 expressed genes (qval \u3c 0.05). Associations between cis-eQTL markers and host response phenotypes using 383 pigs suggest that host genotype-dependent differences in expression of GBP5, GBP6, CCHCR1 and CMPK2 affect viremia levels or weight gain in response to PRRSV infection

Improved taxonomic assignment of rumen bacterial 16S rRNA sequences using a revised SILVA taxonomic framework

The taxonomy and associated nomenclature of many taxa of rumen bacteria are poorly defined within databases of 16S rRNA genes. This lack of resolution results in inadequate definition of microbial community structures, with large parts of the community designated as incertae sedis, unclassified, or uncultured within families, orders, or even classes. We have begun resolving these poorly-defined groups of rumen bacteria, based on our desire to name these for use in microbial community profiling. We used the previously-reported global rumen census (GRC) dataset consisting of >4.5 million partial bacterial 16S rRNA gene sequences amplified from 684 rumen samples and representing a wide range of animal hosts and diets. Representative sequences from the 8,985 largest operational units (groups of sequence sharing >97% sequence similarity, and covering 97.8% of all sequences in the GRC dataset) were used to identify 241 pre-defined clusters (mainly at genus or family level) of abundant rumen bacteria in the ARB SILVA 119 framework. A total of 99 of these clusters (containing 63.8% of all GRC sequences) had no unique or had inadequate taxonomic identifiers, and each was given a unique nomenclature. We assessed this improved framework by comparing taxonomic assignments of bacterial 16S rRNA gene sequence data in the GRC dataset with those made using the original SILVA 119 framework, and three other frameworks. The two SILVA frameworks performed best at assigning sequences to genus-level taxa. The SILVA 119 framework allowed 55.4% of the sequence data to be assigned to 751 uniquely identifiable genus-level groups. The improved framework increased this to 87.1% of all sequences being assigned to one of 871 uniquely identifiable genus-level groups. The new designations were included in the SILVA 123 release (https://www.arb-silva.de/documentation/release-123/) and will be perpetuated in future releases