Matthew S. McCabe

Effect of Phenotypic Residual Feed Intake and Dietary Forage Content on the Rumen Microbial Community of Beef Cattle

ABSTRACT Feed-efficient animals have lower production costs and reduced environmental impact. Given that rumen microbial fermentation plays a pivotal role in host nutrition, the premise that rumen microbiota may contribute to host feed efficiency is gaining momentum. Since diet is a major factor in determining rumen community structure and fermentation patterns, we investigated the effect of divergence in phenotypic residual feed intake (RFI) on ruminal community structure of beef cattle across two contrasting diets. PCR-denaturing gradient gel electrophoresis (DGGE) and quantitative PCR (qPCR) were performed to profile the rumen bacterial population and to quantify the ruminal populations of Entodinium spp., protozoa, Fibrobacter succinogenes, Ruminococcus flavefaciens, Ruminococcus albus, Prevotella brevis, the genus Prevotella, and fungi in 14 low (efficient)- and 14 high (inefficient)-RFI animals offered a low-energy, high-forage diet, followed by a high-energy, low-forage diet. Canonical correspondence and Spearman correlation analyses were used to investigate associations between physiological variables and rumen microbial structure and specific microbial populations, respectively. The effect of RFI on bacterial profiles was influenced by diet, with the association between RFI group and PCR-DGGE profiles stronger for the higher forage diet. qPCR showed that Prevotella abundance was higher (P < 0.0001) in inefficient animals. A higher (P < 0.0001) abundance of Entodinium and Prevotella spp. and a lower (P < 0.0001) abundance of Fibrobacter succinogenes were observed when animals were offered the low-forage diet. Thus, differences in the ruminal microflora may contribute to host feed efficiency, although this effect may also be modulated by the diet offered.

Next Generation Sequencing Reveals the Expression of a Unique miRNA Profile in Response to a Gram-Positive Bacterial Infection

MicroRNAs (miRNAs) are short, non-coding RNAs, which post-transcriptionally regulate gene expression and are proposed to play a key role in the regulation of innate and adaptive immunity. Here, we report a next generation sequencing (NGS) approach profiling the expression of miRNAs in primary bovine mammary epithelial cells (BMEs) at 1, 2, 4 and 6 hours post-infection with Streptococcus uberis, a causative agent of bovine mastitis. Analysing over 450 million sequencing reads, we found that 20% of the approximately 1,300 currently known bovine miRNAs are expressed in unchallenged BMEs. We also identified the expression of more than 20 potentially novel bovine miRNAs. There is, however, a significant dynamic range in the expression of known miRNAs. The top 10 highly expressed miRNAs account for >80% of all aligned reads, with the remaining miRNAs showing much lower expression. Twenty-one miRNAs were identified as significantly differentially expressed post-infection with S. uberis. Several of these miRNAs have characterised roles in the immune systems of other species. This miRNA response to the Gram-positive S. uberis is markedly different, however, to lipopolysaccharide (LPS) induced miRNA expression. Of 145 miRNAs identified in the literature as being LPS responsive, only 9 were also differentially expressed in response to S. uberis. Computational analysis has also revealed that the predicted target genes of miRNAs, which are down-regulated in BMEs following S. uberis infection, are statistically enriched for roles in innate immunity. This suggests that miRNAs, which potentially act as central regulators of gene expression responses to a Gram-positive bacterial infection, may significantly regulate the sentinel capacity of mammary epithelial cells to mobilise the innate immune system.

Transcriptomic analysis of the stress response to weaning at housing in bovine leukocytes using RNA-seq technology

BackgroundWeaning of beef calves is a necessary husbandry practice and involves separating the calf from its mother, resulting in numerous stressful events including dietary change, social reorganisation and the cessation of the maternal-offspring bond and is often accompanied by housing. While much recent research has focused on the physiological response of the bovine immune system to stress in recent years, little is known about the molecular mechanisms modulating the immune response. Therefore, the objective of this study was to provide new insights into the molecular mechanisms underlying the physiological response to weaning at housing in beef calves using Illumina RNA-seq.ResultsThe leukocyte transcriptome was significantly altered for at least 7 days following either housing or weaning at housing. Analysis of differentially expressed genes revealed that four main pathways, cytokine signalling, transmembrane transport, haemostasis and G-protein-coupled receptor (GPRC) signalling were differentially regulated between control and weaned calves and underwent significant transcriptomic alterations in response to weaning stress on day 1, 2 and 7. Of particular note, chemokines, cytokines and integrins were consistently found to be up-regulated on each day following weaning. Evidence for alternative splicing of genes was also detected, indicating a number of genes involved in the innate and adaptive immune response may be alternatively transcribed, including those responsible for toll receptor cascades and T cell receptor signalling.ConclusionsThis study represents the first application of RNA-Seq technology for genomic studies in bovine leukocytes in response to weaning stress. Weaning stress induces the activation of a number of cytokine, chemokine and integrin transcripts and may alter the immune system whereby the ability of a number of cells of the innate and adaptive immune system to locate and destroy pathogens is transcriptionally enhanced. Stress alters the homeostasis of the transcriptomic environment of leukocytes for at least 7 days following weaning, indicating long term effects of stress exposure in the bovine. The identification of gene signature networks that are stress activated provides a mechanistic framework to characterise the multifaceted nature of weaning stress adaptation in beef calves. Thus, capturing subtle transcriptomic changes provides insight into the molecular mechanisms that underlie the physiological response to weaning stress.

Illumina MiSeq Phylogenetic Amplicon Sequencing Shows a Large Reduction of an Uncharacterised Succinivibrionaceae and an Increase of the Methanobrevibacter gottschalkii Clade in Feed Restricted Cattle.

Periodic feed restriction is used in cattle production to reduce feed costs. When normal feed levels are resumed, cattle catch up to a normal weight by an acceleration of normal growth rate, known as compensatory growth, which is not yet fully understood. Illumina Miseq Phylogenetic marker amplicon sequencing of DNA extracted from rumen contents of 55 bulls showed that restriction of feed (70% concentrate, 30% grass silage) for 125 days, to levels that caused a 60% reduction of growth rate, resulted in a large increase of relative abundance of Methanobrevibacter gottschalkii clade (designated as OTU-M7), and a large reduction of an uncharacterised Succinivibrionaceae species (designated as OTU-S3004). There was a strong negative Spearman correlation (ρ = -0.72, P = <1x10-20) between relative abundances of OTU-3004 and OTU-M7 in the liquid rumen fraction. There was also a significant increase in acetate:propionate ratio (A:P) in feed restricted animals that showed a negative Spearman correlation (ρ = -0.69, P = <1x10-20) with the relative abundance of OTU-S3004 in the rumen liquid fraction but not the solid fraction, and a strong positive Spearman correlation with OTU-M7 in the rumen liquid (ρ = 0.74, P = <1x10-20) and solid (ρ = 0.69, P = <1x10-20) fractions. Reduced A:P ratios in the rumen are associated with increased feed efficiency and reduced production of methane which has a global warming potential (GWP 100 years) of 28. Succinivibrionaceae growth in the rumen was previously suggested to reduce methane emissions as some members of this family utilise hydrogen, which is also utilised by methanogens for methanogenesis, to generate succinate which is converted to propionate. Relative abundance of OTU-S3004 showed a positive Spearman correlation with propionate (ρ = 0.41, P = <0.01) but not acetate in the liquid rumen fraction.

Process stability and microbial community composition in pig manure and food waste anaerobic co-digesters operated at low HRTs

This study assessed the effects of reducing hydraulic retention times (HRTs) from 21 days to 10.5 days when anaerobically co-digesting pig manure and food waste. Continuously stirred tank reactors of 3.75 L working volume were operated in triplicate at 42°C. Digester HRT was progressively decreased from 21 to 15 days to 10.5 days, with an associated increase in organic loading rate (OLR) from 3.1 kg volatile solids (VS)·m–3·day–1 to 5.1 kg VS·m–3·day–1 to 7.25 kg VS·m–3·day–1. Reducing HRT from 21 days to 15 days caused a decrease in specific methane yields and VS removal rates. Operation at a HRT of 10.5 days initially resulted in the accumulation of isobutyric acid in each reactor. High throughput 16S rRNA gene sequencing revealed that this increase coincided with a shift in acidogenic bacterial populations, which most likely resulted in the increased isobutyric acid concentrations. This may in turn have caused the increase in relative abundance of Clocamonaceae bacteria, which syntrophically degrade non-acetate volatile fatty acids (VFAs) into H2 and CO2. This, along with the increase in abundance of other syntrophic VFA oxidizers, such as Spiorchatetes, suggests that VFA oxidation plays a role in digester operation at low HRTs. Reducing the HRT to below 21 days compromised the ability of the anaerobic digestion system to reduce enteric indicator organism counts below regulatory limits.

Profiling microRNA expression in bovine alveolar macrophages using RNA-seq

MicroRNAs (miRNAs) are important regulators of gene expression and are known to play a key role in regulating both adaptive and innate immunity. Bovine alveolar macrophages (BAMs) help maintain lung homeostasis and constitute the front line of host defense against several infectious respiratory diseases, such as bovine tuberculosis. Little is known, however, about the role miRNAs play in these cells. In this study, we used a high-throughput sequencing approach, RNA-seq, to determine the expression levels of known and novel miRNAs in unchallenged BAMs isolated from lung lavages of eight different healthy Holstein-Friesian male calves. Approximately 80 million sequence reads were generated from eight BAM miRNA Illumina sequencing libraries, and 80 miRNAs were identified as being expressed in BAMs at a threshold of at least 100 reads per million (RPM). The expression levels of miRNAs varied over a large dynamic range, with a few miRNAs expressed at very high levels (up to 800,000RPM), and the majority lowly expressed. Notably, many of the most highly expressed miRNAs in BAMs have known roles in regulating immunity in other species (e.g. bta-let-7i, bta-miR-21, bta-miR-27, bta-miR-99b, bta-miR-146, bta-miR-147, bta-miR-155 and bta-miR-223). The most highly expressed miRNA in BAMs was miR-21, which has been shown to regulate the expression of antimicrobial peptides in Mycobacterium leprae-infected human monocytes. Furthermore, the predicted target genes of BAM-expressed miRNAs were found to be statistically enriched for roles in innate immunity. In addition to profiling the expression of known miRNAs, the RNA-seq data was also analysed to identify potentially novel bovine miRNAs. One putatively novel bovine miRNA was identified. To the best of our knowledge, this is the first RNA-seq study to profile miRNA expression in BAMs and provides an important reference dataset for investigating the regulatory roles miRNAs play in this important immune cell type.

Effect of Dietary Restriction and Subsequent Re-Alimentation on the Transcriptional Profile of Bovine Skeletal Muscle

Compensatory growth (CG), an accelerated growth phenomenon which occurs following a period of dietary restriction is exploited worldwide in animal production systems as a method to lower feed costs. However the molecular mechanisms regulated CG expression remain to be elucidated fully. This study aimed to uncover the underlying biology regulating CG in cattle, through an examination of skeletal muscle transcriptional profiles utilising next generation mRNA sequencing technology. Twenty Holstein Friesian bulls were fed either a restricted diet for 125 days, with a target growth rate of 0.6 kg/day (Period 1), following which they were allowed feed ad libitum for a further 55 days (Period 2) or fed ad libitum for the entirety of the trial. M. longissimus dorsi biopsies were harvested from all bulls on days 120 and 15 of periods 1 and 2 respectively and RNAseq analysis was performed. During re-alimentation in Period 2, previously restricted animals displayed CG, growing at 1.8 times the rate of the ad libitum control animals. Compensating animals were also more feed efficient during re-alimentation and compensated for 48% of their previous dietary restriction. 1,430 and 940 genes were identified as significantly differentially expressed (Benjamini Hochberg adjusted P < 0.1) in periods 1 and 2 respectively. Additionally, 2,237 genes were differentially expressed in animals undergoing CG relative to dietary restriction. Dietary restriction in Period 1 was associated with altered expression of genes involved in lipid metabolism and energy production. CG expression in Period 2 occurred in association with greater expression of genes involved in cellular function and organisation. This study highlights some of the molecular mechanisms regulating CG in cattle. Differentially expressed genes identified are potential candidate genes for the identification of biomarkers for CG and feed efficiency, which may be incorporated into future breeding programmes.

The Structural and Functional Capacity of Ruminal and Cecal Microbiota in Growing Cattle Was Unaffected by Dietary Supplementation of Linseed Oil and Nitrate

Microorganisms in the digestive tract of ruminants differ in their functionality and ability to use feed constituents. While cecal microbiota play an important role in post-rumen fermentation of residual substrates undigested in the rumen, limited knowledge exists regarding its structure and function. In this trial we investigated the effect of dietary supplementation with linseed oil and nitrate on methane emissions and on the structure of ruminal and cecal microbiota of growing bulls. Animals were allocated to either a CTL (control) or LINNIT (CTL supplemented with 1.9% linseed and 1.0% nitrates) diet. Methane emissions were measured using the GreenFeed system. Microbial diversity was assessed using amplicon sequencing of microbial genomic DNA. Additionally, total RNA was extracted from ruminal contents and functional mcrA and mtt genes were targeted in amplicon sequencing approach to explore the diversity of functional gene expression in methanogens. LINNIT had no effect on methane yield (g/kg DMI) even though it decreased methane production by 9% (g/day; P < 0.05). Methanobrevibacter- and Methanomassiliicoccaceae-related OTUs were more abundant in cecum (72 and 24%) compared to rumen (60 and 11%) irrespective of the diet (P < 0.05). Feeding LINNIT reduced the relative abundance of Methanomassiliicoccaceae mcrA cDNA reads in the rumen. Principal component analysis revealed significant differences in taxonomic composition and abundance of bacterial communities between rumen and cecum. Treatment decreased the relative abundance of a few Ruminococcaceae genera, without affecting global bacterial community structure. Our research confirms a high level of heterogeneity in species composition of microbial consortia in the main gastrointestinal compartments where feed is fermented in ruminants. There was a parallel between the lack of effect of LINNIT on ruminal and cecal microbial community structure and functions on one side and methane emission changes on the other. These results suggest that the sequencing strategy used here to study microbial diversity and function accurately reflected the absence of effect on methane phenotypes in bulls treated with linseed plus nitrate.

Divergent functional isoforms drive niche specialisation for nutrient acquisition and use in rumen microbiome

Many microbes in complex competitive environments share genes for acquiring and utilising nutrients, questioning whether niche specialisation exists and if so, how it is maintained. We investigated the genomic signatures of niche specialisation in the rumen microbiome, a highly competitive, anaerobic environment, with limited nutrient availability determined by the biomass consumed by the host. We generated individual metagenomic libraries from 14 cows fed an ad libitum diet of grass silage and calculated functional isoform diversity for each microbial gene identified. The animal replicates were used to calculate confidence intervals to test for differences in diversity of functional isoforms between microbes that may drive niche specialisation. We identified 153 genes with significant differences in functional isoform diversity between the two most abundant bacterial genera in the rumen (Prevotella and Clostridium). We found Prevotella possesses a more diverse range of isoforms capable of degrading hemicellulose, whereas Clostridium for cellulose. Furthermore, significant differences were observed in key metabolic processes indicating that isoform diversity plays an important role in maintaining their niche specialisation. The methods presented represent a novel approach for untangling complex interactions between microorganisms in natural environments and have resulted in an expanded catalogue of gene targets central to rumen cellulosic biomass degradation.

Differences in leukocyte profile, gene expression, and metabolite status of dairy cows with or without sole ulcers

Sole ulcers are one of the most severe pathologies causing lameness in dairy cows and are associated with abnormal behavior and impaired production performance. However, little is known about how or whether lameness caused by sole ulcers affects the cow systemically. This study compared hematology profile, leukocyte gene expression, and physiological responses [metabolite, cortisol, the endogenous steroid hormone dehydroepiandrosterone (DHEA), and haptoglobin concentrations] of cows with sole ulcers and healthy cows. Twelve clinically lame cows (lame) were identified as having at least one sole ulcer and no other disorder, and matched with a cow that had good locomotion and no disorders (sound), using days in milk, liveweight, body condition score, and diet. Blood samples were taken from all 24 cows within 24h of sole ulcer diagnosis. Leukocyte counts were obtained using an automated cell counter, cortisol and DHEA concentration by ELISA, and plasma haptoglobin, urea, total protein, creatine kinase, and glucose were analyzed on an Olympus analyzer. Expression of 16 genes associated with lameness or stress were estimated using reverse transcription-PCR. Data were analyzed using the MIXED procedure in SAS software (version 9.3; SAS Institute Inc., Cary, NC). Lame cows had a higher neutrophil percentage, a numerically lower lymphocyte percentage, and tended to have a higher neutrophil:lymphocyte ratio than sound cows. Serum cortisol and DHEA concentrations were higher in lame than in sound cows. Lame cows also tended to have higher haptoglobin and glucose levels than sound, as well as higher protein yet lower urea levels. Sound cows tended to have higher relative expression of the gene coding for colony-stimulating factor 2 than lame, but in all other cases where differences were detected in cytokine gene expression (IL-1α, IL-1β, CXCL8, and IL-10), relative gene expression in sound cows tended to be, or was, lower than in lame. Relative expression of MMP-13, GR-α, Fas, haptoglobin, and CD62L were, or tended to be, higher in lame than sound cows. A high neutrophil:lymphocyte ratio in combination with higher cortisol levels in cows with ulcers is indicative of physiological stress. Moreover, increased DHEA and a higher cortisol:DHEA ratio, as well as a tendency for higher haptoglobin levels and increased haptoglobin mRNA expression, are indicative of systemic inflammation. Increased cytokine mRNA expression indicates activation of the immune system compared with healthy cows. Increased expression of MMP-13 mRNA has been found in cows with impaired locomotion and thus could be implicated in development of claw horn disorders.