The effect of breed and diet type on the global transcriptome of hepatic tissue in beef cattle divergent for feed efficiency

Background

Feed efficiency is an important economic and environmental trait in beef production, which can be measured in terms of residual feed intake (RFI). Cattle selected for low-RFI (feed efficient) have similar production levels but decreased feed intake, while also emitting less methane. RFI is difficult and expensive to measure and is not widely adopted in beef production systems. However, development of DNA-based biomarkers for RFI may facilitate its adoption in genomic-assisted breeding programmes. Cattle have been shown to re-rank in terms of RFI across diets and age, while also RFI varies by breed. Therefore, we used RNA-Seq technology to investigate the hepatic transcriptome of RFI-divergent Charolais (CH) and Holstein-Friesian (HF) steers across three dietary phases to identify genes and biological pathways associated with RFI regardless of diet or breed.

Results

Residual feed intake was measured during a high-concentrate phase, a zero-grazed grass phase and a final high-concentrate phase. In total, 322 and 33 differentially expressed genes (DEGs) were identified across all diets for CH and HF steers, respectively. Three genes, GADD45G, HP and MID1IP1, were differentially expressed in CH when both the high-concentrate zero-grazed grass diet were offered. Two canonical pathways were enriched across all diets for CH steers. These canonical pathways were related to immune function.

Conclusions

The absence of common differentially expressed genes across all dietary phases and breeds in this study supports previous reports of the re-ranking of animals in terms of RFI when offered differing diets over their lifetime. However, we have identified biological processes such as the immune response and lipid metabolism as potentially associated with RFI divergence emphasising the previously reported roles of these biological processes with respect to RFI.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5906-8) contains supplementary material, which is available to authorized users.

Effect of Early Calf-Hood Nutrition on the Transcriptional Regulation of the Hypothalamic-Pituitary-Testicular axis in Holstein-Friesian Bull Calves

The aim of this study was to investigate the effect of early calf-hood nutrition on the transcriptomic profile of the arcuate nucleus of the hypothalamus, anterior pituitary and testes in Holstein-Friesian bulls. Holstein-Friesian bull calves with a mean (±S.D.) age and bodyweight of 19 (±8.2) days and 47.5 (±5.3) kg, respectively, were offered a high (n = 10) or low (n = 10) plane of nutrition in order to achieve an overall growth rate of 1.2 and 0.5 kg/day. At 126 (±3) days of age, calves were euthanized, hypothalamus (arcuate region), anterior pituitary and testicular parenchyma samples were harvested and RNAseq analysis was performed. There were 0, 49 and 1,346 genes differentially expressed in the arcuate nucleus, anterior pituitary and testicular tissue of bull calves on the low relative to the high plane of nutrition, respectively (P < 0.05; False Discovery Rate <0.05). Cell cycle processes in the anterior pituitary were down regulated in the low relative to the high plane of nutrition; there was no differential expression of genes related to reproductive processes. Gene expression involved in cholesterol and androgen biosynthesis in the testes were down regulated in animals on the low plane of nutrition. This study provides insight into the effect of early life plane of nutrition on the regulation of the HPT axis.

Blood immune transcriptome analysis of artificially fed dairy calves and naturally suckled beef calves from birth to 7 days of age

Neonatal calves possess a very immature and naïve immune system and are reliant on the intake of maternal colostrum for passive transfer of immunoglobulins. Variation in colostrum management of beef and dairy calves is thought to affect early immune development. Therefore, the objective of this study was to examine changes in gene expression and investigate molecular pathways involved in the immune-competence development of neonatal Holstein dairy calves and naturally suckled beef calves using next generation RNA-sequencing during the first week of life. Jugular whole blood samples were collected from Holstein (H) dairy calves (n = 8) artificially fed 5% B.W. colostrum, and from beef calves which were the progenies of Charolais-Limousin (CL; n = 7) and Limousin-Friesian beef suckler cows (LF; n = 7), for subsequent RNA isolation. In dairy calves, there was a surge in pro-inflammatory cytokine gene expression possibly due to the stress of separation from the dam. LF calves exhibited early signs of humoral immune development with observed increases in the expression genes coding for Ig receptors, which was not evident in the other breeds by 7 days of age. Immune and health related DEGs identified as upregulated in beef calves are prospective contender genes for the classification of biomarkers for immune-competence development, and will contribute towards a greater understanding of the development of an immune response in neonatal calves.

Effect of a butyrate-fortified milk replacer on gastrointestinal microbiota and products of fermentation in artificially reared dairy calves at weaning

Enrichment of calf diets with exogenous butyrate has shown promise as a promotor of calf growth and intestinal development. However, the impact of dietary derived butyrate on the gut microbiota and their potential role, in turn, as mediators of its effect on calf growth and development is not known. Here, the effects of butyrate supplementation on rumen and hindgut microbiota and fermentation profiles were assessed in 16 Holstein-Friesian bull calves randomly assigned to one of two groups: Control (CON) fed conventional milk replacer or Sodium-Butyrate (SB – added to milk replacer) from days 7 to 56 of life. In the colon, total short chain fatty acid (SCFA), propionate and acetate concentrations were increased by SB (P < 0.05). 16S rRNA gene amplicon sequencing showed cecal abundance of butyrate producers Butyrivibrio and Shuttleworthia were decreased by SB (P < 0.05), while that of the propionate producer Phascolarctobacterium was higher (P < 0.05). Mogibacterium is associated with impaired gut health and was reduced in the cecum of SB calves (P < 0.05). These data show that the beneficial effects of SB on growth and performance occur in tandem with changes in the abundance of important SCFA producing and health-associated bacteria in the hindgut in milk-fed calves.

Role of early life nutrition on regulating the hypothalamic–anterior pituitary–testicular axis of the bull

The objective of this study was to examine the effect of nutrition during the first 18 weeks of life on the physiological and transcriptional functionality of the hypothalamic (arcuate nucleus region), anterior pituitary and testes in Holstein–Friesian bull calves. Holstein–Friesian bull calves with a mean (±S.D.) age and bodyweight of 19 (±8.2) days and 47.5 (±5.3) kg, respectively, were assigned to either a HIGH (n = 10) or LOW (n = 10) plane of nutrition, to achieve an overall target growth rate of 1.2 or 0.5 kg/day, respectively. At 126 ± 1.1 days of age, all calves were euthanised. Animal performance (weekly) and systemic concentrations of metabolic (monthly) and reproductive hormones (fortnightly) were assessed. Testicular histology, targeted gene and protein expression of the arcuate nucleus region, anterior pituitary and testes were also assessed using qPCR and immunohistochemistry, respectively. The expression of candidate genes in testicular tissue from post pubertal 19-month-old Holstein–Friesian bulls (n = 10) was compared to that of the 18-week-old calves. Metabolite and metabolic hormone profiles generally reflected the improved metabolic status of the calves on the HIGH (P< 0.001). Calves offered a HIGH plane of nutrition were heavier at slaughter (P < 0.001), had larger testes (P < 0.001), larger seminiferous tubule diameter (P < 0.001), more mature spermatogenic cells (P < 0.001) and more Sertoli cells (P < 0.05) in accordance with both morphological and transcriptional data. Overall, testicular gene expression profiles suggested a more mature stage of development in HIGH compared with LOW and were more closely aligned to that of mature bulls. Ghrelin receptor was the only differentially expressed gene between LOW and HIGH calves in either the anterior pituitary (P < 0.05) or arcuate nucleus region of the hypothalamus (P < 0.10) and was upregulated in LOW for both tissues. This study indicates that an enhanced plane of nutrition during early calfhood favourably alters the biochemical regulation of the hypothalamus–anterior pituitary–testicular axis, advancing testicular development and hastening spermatogenesis.

GWAS and eQTL analysis identifies a SNP associated with both residual feed intake and GFRA2 expression in beef cattle

Residual feed intake (RFI), a measure of feed efficiency, is an important economic and environmental trait in beef production. Selection of low RFI (feed efficient) cattle could maintain levels of production, while decreasing feed costs and methane emissions. However, RFI is a difficult and expensive trait to measure. Identification of single nucleotide polymorphisms (SNPs) associated with RFI may enable rapid, cost effective genomic selection of feed efficient cattle. Genome-wide association studies (GWAS) were conducted in multiple breeds followed by meta-analysis to identify genetic variants associated with RFI and component traits (average daily gain (ADG) and feed intake (FI)) in Irish beef cattle (n = 1492). Expression quantitative trait loci (eQTL) analysis was conducted to identify functional effects of GWAS-identified variants. Twenty-four SNPs were associated (P < 5 × 10^-5) with RFI, ADG or FI. The variant rs43555985 exhibited strongest association for RFI (P = 8.28E-06). An eQTL was identified between this variant and GFRA2 (P = 0.0038) where the allele negatively correlated with RFI was associated with increased GFRA2 expression in liver. GFRA2 influences basal metabolic rates, suggesting a mechanism by which genetic variation may contribute to RFI. This study identified SNPs that may be useful both for genomic selection of RFI and for understanding the biology of feed efficiency.

RNA-seq analysis of bovine adipose tissue in heifers fed diets differing in energy and protein content

Adipose tissue is no longer considered a mere energy reserve, but a metabolically and hormonally active organ strongly associated with the regulation of whole-body metabolism. Knowledge of adipose metabolic regulatory function is of great importance in cattle management, as it affects the efficiency and manner with which an animal converts feedstuff to milk, meat and fat. However, the molecular mechanisms regulating metabolism in bovine adipose tissue are still not fully elucidated. The emergence of next-generation sequencing technologies has facilitated the analysis of metabolic function and regulation at the global gene expression level. The aim of this study was to investigate the effect of diets differing in protein and energy density level on gene expression in adipose tissue of growing replacement dairy heifers using next-generation RNA sequencing (RNAseq). Norwegian Red heifers were fed either a high- or low-protein concentrate (HP/LP) and a high- or low-energy roughage (HE/LE) diet from 3 months of age until confirmed pregnancy to give four treatments (viz, HPHE, HPLE, LPHE, LPLE) with different growth profiles. Subcutaneous adipose tissue sampled at 12 months of age was analyzed for gene expression differences using RNAseq. The largest difference in gene expression was found between LPHE and LPLE heifers, for which 1092 genes were significantly differentially expressed, representing an up-regulation of mitochondrial function, lipid, carbohydrate and amino acid metabolism as well as changes in the antioxidant system in adipose tissue of LPHE heifers. Differences between HPHE and HPLE heifers were much smaller, and dominated by genes representing NAD biosynthesis, as was the significantly differentially expressed genes (DEG) common to both HE-LE contrasts. Differences between HP and LP groups within each energy treatment were minimal. This study emphasizes the importance of transcriptional regulation of adipose tissue energy metabolism, and identifies candidate genes for further studies on early-stage obesity and glucose load in dairy cattle.

16S rRNA Sequencing Reveals Relationship Between Potent Cellulolytic Genera and Feed Efficiency in the Rumen of Bulls

The rumen microbial population dictates the host’s feed degradation capacity and subsequent nutrient supply. The rising global human population and intensifying demand for animal protein is creating environmental challenges. As a consequence, there is an increasing requirement for livestock with enhanced nutrient utilization capacity in order to more efficiently convert plant material to high quality edible muscle. In the current study, residual feed intake (RFI), a widely used and a highly accepted measure of feed efficiency in cattle, was calculated for a combination of three cohorts of Simmental bulls. All animals were managed similarly from birth and offered concentrate ad libitum in addition to 3 kg of grass silage daily during the finishing period. Solid and liquid rumen digesta samples collected at slaughter and were analyzed using amplicon sequencing targeting the 16S rRNA gene utilizing the Illumina MiSeq platform. Volatile fatty acid analysis was also conducted on the liquid digesta samples. Spearman’s correlation coefficient was utilized to determine the association between RFI and bacterial and archaeal taxa and inter-taxonomic relationships. The data indicate a tendency toward an increase in butyrate (P = 0.06), which corresponds with an increase in plasma β-hydroxybutyrate concentration in low RFI (LRFI) bulls in comparison to their high RFI (HRFI) contemporaries (P < 0.05). A decrease in propionate (P < 0.05) was also recorded in the rumen of LRFI in comparison to HRFI bulls. These results indicate alternate fermentation patterns in the rumen of LRFI bulls. The data also identified that OTUs within the phyla Tenericutes, Fibrobacteres, and Cyanobacteria may potentially influence RFI phenotype. In particular, a negative association between F. succinogenes and RFI was evident. The unique cellulolytic metabolism of F. succinogenes suggests it could contribute to host efficiency by providing substrate to the host ruminant and other microbial populations (e.g., Selenomonas ruminantium, Methanobrevibacter, and Methanomassiliicoccaceae) in the rumen. This study provides evidence that bacterial OTUs within common phyla could influence ruminant feed efficiency phenotype through their role in ruminal degradation of complex plant polysaccharides or increased capability to harvest nutrients from ingested feed.

Evaluating Established Methods for Rumen 16S rRNA Amplicon Sequencing With Mock Microbial Populations

The rumen microbiome scientific community has utilized amplicon sequencing as an aid in identifying potential community compositional trends that could be used as an estimation of various production and performance traits including methane emission, animal protein production efficiency, and ruminant health status. In order to translate rumen microbiome studies into executable application, there is a need for experimental and analytical concordance within the community. The objective of this study was to assess these factors in relation to selected currently established methods for 16S phylogenetic community analysis on a microbial community standard (MC) and a DNA standard (DS; ZymoBIOMICS^TM). DNA was extracted from MC using the RBBC method commonly used for microbial DNA extraction from rumen digesta samples. 16S rRNA amplicon libraries were generated for the MC and DS using primers routinely used for rumen bacterial and archaeal community analysis. The primers targeted the V4 and V3–V4 region of the 16S rRNA gene and samples were subjected to both 20 and 28 polymerase chain reaction (PCR) cycles under identical cycle conditions. Sequencing was conducted using the Illumina MiSeq platform. As the bacteria contained in the microbial mock community were well-classified species, and for ease of explanation, we used the results of the Basic Local Alignment Search Tool classification to assess the DNA, PCR cycle number, and primer type. Sequence classification methodology was assessed independently. Spearman’s correlation analysis indicated that utilizing the repeated bead beating and column method for DNA extraction in combination with primers targeting the 16S rRNA gene using 20 first-round PCR cycles was sufficient for amplicon sequencing to generate a relatively accurate depiction of the bacterial communities present in rumen samples. These results also emphasize the requirement to develop and utilize positive mock community controls for all rumen microbiomic studies in order to discern errors which may arise at any step during a next-generation sequencing protocol.

Transcriptome analyses reveal reduced hepatic lipid synthesis and accumulation in more feed efficient beef cattle

The genetic mechanisms controlling residual feed intake (RFI) in beef cattle are still largely unknown. Here we performed whole transcriptome analyses to identify differentially expressed (DE) genes and their functional roles in liver tissues between six extreme high and six extreme low RFI steers from three beef breed populations including Angus, Charolais, and Kinsella Composite (KC). On average, the next generation sequencing yielded 34 million single-end reads per sample, of which 87% were uniquely mapped to the bovine reference genome. At false discovery rate (FDR) < 0.05 and fold change (FC) > 2, 72, 41, and 175 DE genes were identified in Angus, Charolais, and KC, respectively. Most of the DE genes were breed-specific, while five genes including TP53INP1, LURAP1L, SCD, LPIN1, and ENSBTAG00000047029 were common across the three breeds, with TP53INP1, LURAP1L, SCD, and LPIN1 being downregulated in low RFI steers of all three breeds. The DE genes are mainly involved in lipid, amino acid and carbohydrate metabolism, energy production, molecular transport, small molecule biochemistry, cellular development, and cell death and survival. Furthermore, our differential gene expression results suggest reduced hepatic lipid synthesis and accumulation processes in more feed efficient beef cattle of all three studied breeds.

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.

Pregnancy losses in cattle: potential for improvement

For heifers, beef and moderate-yielding dairy cows, it appears that the fertilisation rate generally lies between 90% and 100%. For high-producing dairy cows, there is a less substantive body of literature, but it would appear that the fertilisation rate is somewhat lower and possibly more variable. In cattle, the major component of embryo loss occurs in the first 16 days following breeding (Day 0), with emerging evidence of greater losses before Day 8 in high-producing dairy cows. In cattle, late embryo mortality causes serious economic losses because it is often recognised too late to rebreed females. Systemic concentrations of progesterone during both the cycle preceding and following insemination affect embryo survival, with evidence of either excessive or insufficient concentrations being negatively associated with survival rate. The application of direct progesterone supplementation or treatments to increase endogenous output of progesterone to increase embryo survival cannot be recommended at this time. Energy balance and dry matter intake during the first 4 weeks after calving are critically important in determining pregnancies per AI when cows are inseminated at 70-100 days after calving. Level of concentrate supplementation of cows at pasture during the breeding period has minimal effects on conception rates, although sudden reductions in dietary intake should be avoided. For all systems of milk production, more balanced breeding strategies with greater emphasis on fertility and feed intake and/or energy must be developed. There is genetic variability within the Holstein breed for fertility traits, which can be exploited. Genomic technology will not only provide scientists with an improved understanding of the underlying biological processes involved in fertilisation and the establishment of pregnancy, but also, in the future, could identify genes responsible for improved embryo survival. Such information could be incorporated into breeding objectives in order to increase the rate of genetic progress for embryo survival. In addition, there is a range of easily adoptable management factors, under producer control, that can either directly increase embryo survival or ameliorate the consequences of low embryo survival rates. The correction of minor deficits in several areas can have a substantial cumulative positive effect on herd reproductive performance.

Examination of the molecular control of ruminal epithelial function in response to dietary restriction and subsequent compensatory growth in cattle

BACKGROUND

The objective of this study was to investigate the effect of dietary restriction and subsequent compensatory growth on the relative expression of genes involved in volatile fatty acid transport, metabolism and cell proliferation in ruminal epithelial tissue of beef cattle. Sixty Holstein Friesian bulls (mean liveweight 370 ± 35 kg; mean age 479 ± 15 d) were assigned to one of two groups: (i) restricted feed allowance (RES; n = 30) for 125 d (Period 1) followed by ad libitum access to feed for 55 d (Period 2) or (ii) ad libitum access to feed throughout (ADLIB; n = 30). Target growth rate for RES was 0.6 kg/d during Period 1. At the end of each dietary period, 15 animals from each treatment group were slaughtered and ruminal epithelial tissue and liquid digesta harvested from the ventral sac of the rumen. Real-time qPCR was used to quantify mRNA transcripts of 26 genes associated with ruminal epithelial function. Volatile fatty acid analysis of rumen fluid from individual animals was conducted using gas chromatography.

RESULTS

Diet × period interactions were evident for genes involved in ketogenesis (BDH2, P = 0.017), pyruvate metabolism (LDHa, P = 0.048; PDHA1, P = 0.015) and cellular transport and structure (DSG1, P = 0.019; CACT, P = 0.027). Ruminal concentrations of propionic acid (P = 0.018) and n-valeric acid (P = 0.029) were lower in RES animals, compared with ADLIB, throughout the experiment. There was also a strong tendency (P = 0.064) toward a diet × period interaction for n-butyric with higher concentrations in RES animals, compared with ADLIB, during Period 1.

CONCLUSIONS

These data suggest that following nutrient restriction, the structural integrity of the rumen wall is compromised and there is upregulation of genes involved in the production of ketone bodies and breakdown of pyruvate for cellular energy. These results provide an insight into the potential molecular mechanisms regulating ruminal epithelial absorptive metabolism and growth following nutrient restriction and subsequent compensatory growth.

Effect of manipulating progesterone before timed artificial insemination on reproductive and endocrine parameters in seasonal-calving, pasture-based Holstein-Friesian cows

Fertility to timed AI (TAI) is profoundly affected by progesterone (P4) levels during hormonal synchronization protocols. Holstein-Friesian dairy cows managed in a seasonal-calving, pasture-based production system were randomly assigned to 2 treatments to manipulate P4 before TAI during growth of the preovulatory follicle. Cows in the first treatment (High P4; n=30) were submitted to a Double-Ovsynch protocol {Pre-Ovsynch [GnRH; 7 d, PGF2α; 3 d, GnRH] followed 7 d later by Breeding-Ovsynch [GnRH (G1); 7 d PGF2α; 24 h, PGF2α; 32 h, GnRH (G2); 16 h, TAI]}. Cows in the second treatment (n=30; Low P4) received the same Double-Ovsynch protocol but with an additional PGF2α treatment 24 h after G1. Overall, synchronization rate did not differ between treatments and was 92% (55/60). Unexpectedly, 37% of Low P4 cows were detected in estrus ~24 h before scheduled TAI and were inseminated ~16 h before scheduled TAI. Overall, P4 did not differ between treatments at G1, whereas High P4 cows had greater P4 concentrations at PGF2α and G2 than Low P4 cows. High P4 cows had the smallest mean follicle diameter at G2, whereas Low P4 cows with no estrus before TAI had intermediate mean follicle diameter at G2, and Low P4 cows with estrus before TAI had the largest mean follicle diameter. Low P4 cows with estrus before TAI had larger corpora lutea 15 d after TAI than Low P4 cows without estrus before TAI or High P4 cows. In accordance with corpus luteum size on d 15, High P4 cows and Low P4 cows without estrus before TAI had lower P4 from 4 to 46 d after TAI than Low P4 cows with estrus before TAI. Relative mRNA levels of the interferon-stimulated genes ISG15, MX1, MX2, and OAS1 were greater for Low P4 than for High P4 cows, whereas relative mRNA levels of RTP4 were greater for High P4 than for Low P4 cows 18 d after TAI. Treatment did not affect plasma pregnancy-associated glycoprotein concentrations after TAI; however, pregnancy-associated glycoprotein concentrations were affected by pregnancy status and parity. Treatment did not affect pregnancy per artificial insemination at 29, 39, or 60 d after TAI, and no pregnancy losses were observed from 39 to 60 d after TAI. We concluded that (1) Low P4 cows were more likely to express estrus than High P4 cows; (2) the subpopulation of Low P4 cows that expressed estrus had larger preovulatory follicles and greater P4 concentrations after TAI; and (3) regardless of estrus before TAI, all Low P4 cows had greater mRNA expression for 5 of 6 interferon-stimulated genes than High P4 cows 18 d after TAI.

Endometrial gene expression in high- and low-fertility heifers in the late luteal phase of the estrous cycle and a comparison with midluteal gene expression

Embryonic mortality is a major constraint to improving reproductive efficiency and profitability in livestock enterprises. We previously reported differential expression of genes with identified roles in cellular growth and proliferation, lipid metabolism, endometrial remodeling, inflammation, angiogenesis, and metabolic exchange in endometrial tissue on day 7 of the estrous cycle (D7), between heifers ranked as either high (HF) or low (LF) for fertility. The aim of the current study was to further elucidate the underlying molecular mechanisms contributing to early embryo loss by examining differential endometrial gene expression in HF or LF heifers at a later stage of the estrous cycle;day 14(D14). A second objective was to compare these expression profiles with those from midluteal HF and LF endometrium. Using the same animal model as employed in the previous study, we slaughtered HF and LF animals on D14, harvested endometrial tissue, and carried out global gene expression analysis using the Affymetrix Bovine GeneChip. Microarray analysis detected 430 differentially expressed genes (DEG) between HF and LF animals. Ingenuity Pathway Analysis revealed enrichment for a host of biological pathways including lipid metabolism, molecular transport, immune response, cell morphology and development, and cell growth and proliferation. Important DEG includedALB, BMPR2, CCL28, COL4A3/4, FADS1, ITGA6, LDLR, PLCB3, PPARG, PTGS2, and SLC27A4 Furthermore, DEG expressed on both D7 and D14 included:PCCB,SLC25A24,DAP, and COL4A4 This study highlights some of the pathways and mechanisms underpinning late luteal bovine endometrial physiology and endometrial-related conception rate variance.

Insulin secretion and signaling in response to dietary restriction and subsequent re-alimentation in cattle

The objectives of this study were to examine systemic insulin response to a glucose tolerance test (GTT) and transcript abundance of genes of the insulin signaling pathway in skeletal muscle, during both dietary restriction and re-alimentation-induced compensatory growth. Holstein Friesian bulls were blocked to one of two groups: 1) restricted feed allowance for 125 days (period 1) (RES, n = 15) followed by ad libitum feeding for 55 days (period 2) or 2) ad libitum access to feed throughout (periods 1 and 2) (ADLIB, n = 15). On days 90 and 36 of periods 1 and 2, respectively, a GTT was performed. M. longissimus dorsi biopsies were harvested from all bulls on days 120 and 15 of periods 1 and 2, respectively, and RNA-Seq analysis was performed. RES displayed a lower growth rate during period 1 (RES: 0.6 kg/day, ADLIB: 1.9 kg/day; P < 0.001), subsequently gaining more during re-alimentation (RES: 2.5 kg/day, ADLIB: 1.4 kg/day; P < 0.001). Systemic insulin response to glucose administration was lower in RES in period 1 (P < 0.001) with no difference observed during period 2. The insulin signaling pathway in M. longissimus dorsi was enriched (P < 0.05) in response to dietary restriction but not during re-alimentation (P > 0.05). Genes differentially expressed in the insulin signaling pathway suggested a greater sensitivity to insulin in skeletal muscle, with pleiotropic effects of insulin signaling interrupted during dietary restriction. Collectively, these results indicate increased sensitivity to glucose clearance and skeletal muscle insulin signaling during dietary restriction; however, no overall role for insulin was apparent in expressing compensatory growth.

Effect of feed restriction and subsequent re-alimentation on hormones and genes of the somatotropic axis in cattle

The objective of this study was to characterize the effect of feed restriction and compensatory growth during re-alimentation on the functionality of the somatotropic axis. We blocked 60 bulls into one of two groups: 1) restricted feed allowance for 125 days (period 1) (RES, n = 30) followed by ad libitum feeding for 55 days (period 2) or 2) ad libitum access to feed throughout (ADLIB, n = 30). A growth hormone releasing hormone (GHRH) challenge was performed during each period. At the end of each period, 15 animals from each treatment were slaughtered and hepatic tissue collected. Hepatic expression of 13 genes of the somatotropic axis was measured by qRT-PCR. RES displayed a lower growth rate during period 1 (0.6 vs. 1.9 kg/day; P < 0.001), subsequently gaining more than ADLIB animals during period 2 (2.5 vs. 1.4 kg/day; P < 0.001). Growth hormone response to GHRH was not different between treatments at either time-point (P > 0.05); however, resultant plasma IGF-1 was lower in period 1 and greater in period 2 in RES animals (P < 0.05). Expression of IGFBP2 was higher (P < 0.01) and IGF1 (P < 0.001) and GHRIA (P < 0.05) lower in RES compared with ADLIB during period 1, with no difference evident in period 2 (P > 0.05). Collectively, the results of this study are consistent with uncoupling of the somatotropic axis following feed restriction. However, there is no evidence from this study that the somatotropic axis per se is a significant contributor to compensatory growth.

Diversity and community composition of methanogenic archaea in the rumen of Scottish upland sheep assessed by different methods

Ruminal archaeomes of two mature sheep grazing in the Scottish uplands were analysed by different sequencing and analysis methods in order to compare the apparent archaeal communities. All methods revealed that the majority of methanogens belonged to the Methanobacteriales order containing the Methanobrevibacter, Methanosphaera and Methanobacteria genera. Sanger sequenced 1.3 kb 16S rRNA gene amplicons identified the main species of Methanobrevibacter present to be a SGMT Clade member Mbb. millerae (≥91% of OTUs); Methanosphaera comprised the remainder of the OTUs. The primers did not amplify ruminal Thermoplasmatales-related 16S rRNA genes. Illumina sequenced V6–V8 16S rRNA gene amplicons identified similar Methanobrevibacter spp. and Methanosphaera clades and also identified the Thermoplasmatales-related order as 13% of total archaea. Unusually, both methods concluded that Mbb. ruminantium and relatives from the same clade (RO) were almost absent. Sequences mapping to rumen 16S rRNA and mcrA gene references were extracted from Illumina metagenome data. Mapping of the metagenome data to16S rRNA gene references produced taxonomic identification to Order level including 2–3% Thermoplasmatales, but was unable to discriminate to species level. Mapping of the metagenome data to mcrA gene references resolved 69% to unclassified Methanobacteriales. Only 30% of sequences were assigned to species level clades: of the sequences assigned to Methanobrevibacter, most mapped to SGMT (16%) and RO (10%) clades. The Sanger 16S amplicon and Illumina metagenome mcrA analyses showed similar species richness (Chao1 Index 19–35), while Illumina metagenome and amplicon 16S rRNA analysis gave lower richness estimates (10–18). The values of the Shannon Index were low in all methods, indicating low richness and uneven species distribution. Thus, although much information may be extracted from the other methods, Illumina amplicon sequencing of the V6–V8 16S rRNA gene would be the method of choice for studying rumen archaeal communities.

Quantitative analysis of ruminal methanogenic microbial populations in beef cattle divergent in phenotypic residual feed intake (RFI) offered contrasting diets

Background

Methane (CH₄) emissions in cattle are an undesirable end product of rumen methanogenic fermentative activity as they are associated not only with negative environmental impacts but also with reduced host feed efficiency. The aim of this study was to quantify total and specific rumen microbial methanogenic populations in beef cattle divergently selected for residual feed intake (RFI) while offered (i) a low energy high forage (HF) diet followed by (ii) a high energy low forage (LF) diet. Ruminal fluid was collected from 14 high (H) and 14 low (L) RFI animals across both dietary periods. Quantitative real time PCR (qRT-PCR) analysis was conducted to quantify the abundance of total and specific rumen methanogenic microbes. Spearman correlation analysis was used to investigate the association between the relative abundance of methanogens and animal performance, rumen fermentation variables and diet digestibility.

Results

Abundance of methanogens, did not differ between RFI phenotypes. However, relative abundance of total and specific methanogen species was affected (P < 0.05) by diet type, with greater abundance observed while animals were offered the LF compared to the HF diet.

Conclusions

These findings suggest that differences in abundance of specific rumen methanogen species may not contribute to variation in CH₄ emissions between efficient and inefficient animals, however dietary manipulation can influence the abundance of total and specific methanogen species.

Effects of phenotypic residual feed intake on response to a glucose tolerance test and gene expression in the insulin signaling pathway in longissimus dorsi in beef cattle

The objectives of this study were to determine the insulinogenic response to an intravenous glucose tolerance test (GTT) and examine gene expression profiles in the insulin signaling pathway (ISP) in beef animals of differing phenotypic residual feed intake (RFI). Two experiments were conducted. In Exp. 1, a total of 39 Simmental heifers, over 2 yr (yr 1, n = 22, and yr 2, n = 17; mean initial BW = 472 kg [SD = 52.4 kg]), were offered grass silage ad libitum for 104 d. Heifers were subjected to a GTT on d 8 and 65 of the RFI measurement period in yr 1 and 2, respectively. Concentrations of plasma glucose and insulin were measured at -45, -30, -15, 0, 5, 10, 15, 20, 30, 45, 60, 90, 120, 150, and 180 min relative to glucose infusion (0 min). In Exp. 2, a total of 67 Simmental bulls, over 3 yr (yr 1, n = 20; yr 2, n = 33; and yr 3, n = 14; mean initial BW = 431 kg [SD = 63.7 kg]), were offered concentrates ad libitum for 105 d. Biopsies of LM were harvested during the RFI measurement period (yr 1, d 49 and 91; yr 2, d 52 and 92; and yr 3, d 50 and 92). The residuals of the regression of DMI on ADG, midtest metabolic BW (BW(0.75)), and the fixed effect of year, using all animals, were used to compute individual RFI coefficients. Animals were ranked on RFI and assigned to high (inefficient), medium, or low groupings by dividing them into terciles, resulting in 13 heifers and 22, 23, and 22 bulls in their respective RFI groups. In Exp. 1, data from 13 heifers from each of the high- and low-RFI groups, and in Exp. 2, data from the 15 highest and 15 lowest ranking bulls on RFI are reported. In Exp. 1, glucose and insulin response and area under the response curve for glucose and insulin were similar (P > 0.05) between high- and low-RFI heifers. In Exp. 2, no differences (P > 0.05) were found for mRNA expression of 22 genes of the ISP in muscle tissue; however, expression of the transcription factor SREBP1c tended to be positively correlated (r = 0.25, P = 0.07) with RFI. Expression of GLUT4, INPPL1, and SHC increased (P < 0.05) over time, while there was no effect of sample time for any other genes measured. Collectively, these results suggest that insulin response, sensitivity, and associated expression of genes in the ISP within muscle tissue are not contributory factors to variation in RFI. However, further examination of target genes of SREBP1c, which is involved in lipogenesis, may explain some of the biochemical processes underlying variation in phenotypic RFI.

Gastrointestinal tract size, total-tract digestibility, and rumen microflora in different dairy cow genotypes

The superior milk production efficiency of Jersey (JE) and Jersey × Holstein-Friesian (JE × HF) cows compared with Holstein-Friesian (HF) has been widely published. The biological differences among dairy cow genotypes, which could contribute to the milk production efficiency differences, have not been as widely studied however. A series of component studies were conducted using cows sourced from a longer-term genotype comparison study (JE, JE × HF, and HF). The objectives were to (1) determine if differences exist among genotypes regarding gastrointestinal tract (GIT) weight, (2) assess and quantify whether the genotypes tested differ in their ability to digest perennial ryegrass, and (3) examine the relative abundance of specific rumen microbial populations potentially relating to feed digestibility. Over 3 yr, the GIT weight was obtained from 33 HF, 35 JE, and 27 JE × HF nonlactating cows postslaughter. During the dry period the cows were offered a perennial ryegrass silage diet at maintenance level. The unadjusted GIT weight was heavier for the HF than for JE and JE × HF. When expressed as a proportion of body weight (BW), JE and JE × HF had a heavier GIT weight than HF. In vivo digestibility was evaluated on 16 each of JE, JE × HF, and HF lactating dairy cows. Cows were individually stalled, allowing for the total collection of feces and were offered freshly cut grass twice daily. During this time, daily milk yield, BW, and dry matter intake (DMI) were greater for HF and JE × HF than for JE; milk fat and protein concentration ranked oppositely. Daily milk solids yield did not differ among the 3 genotypes. Intake capacity, expressed as DMI per BW, tended to be different among treatments, with JE having the greatest DMI per BW, HF the lowest, and JE × HF being intermediate. Production efficiency, expressed as milk solids per DMI, was higher for JE than HF and JE × HF. Digestive efficiency, expressed as digestibility of dry matter, organic matter, N, neutral detergent fiber, and acid detergent fiber, was higher for JE than HF. In grazing cows (n=15 per genotype) samples of rumen fluid, collected using a transesophageal sampling device, were analyzed to determine the relative abundance of rumen microbial populations of cellulolytic bacteria, protozoa, and fungi. These are critically important for fermentation of feed into short-chain fatty acids. A decrease was observed in the relative abundance of Ruminococcus flavefaciens in the JE rumen compared with HF and JE × HF. We can deduce from this study that the JE genotype has greater digestibility and a different rumen microbial population than HF. Jersey and JE × HF cows had a proportionally greater GIT weight than HF. These differences are likely to contribute to the production efficiency differences among genotypes previously reported.

Global gene expression in endometrium of high and low fertility heifers during the mid-luteal phase of the estrous cycle

BACKGROUND

In both beef and dairy cattle, the majority of early embryo loss occurs within the first 14 days following insemination. During this time-period, embryos are completely dependent on their maternal uterine environment for development, growth and ultimately survival, therefore an optimum uterine environment is critical to their survival. The objective of this study was to investigate whether differences in endometrial gene expression during the mid-luteal phase of the estrous cycle exist between crossbred beef heifers ranked as either high (HF) or low fertility (LF) (following four rounds of artificial insemination (AI)) using the Affymetrix® 23 K Bovine Gene Chip.

RESULTS

Conception rates for each of the four rounds of AI were within a normal range: 70-73.3%. Microarray analysis of endometrial tissue collected on day 7 of the estrous cycle detected 419 differentially expressed genes (DEG) between HF (n = 6) and LF (n = 6) animals. The main gene pathways affected were, cellular growth and proliferation, angiogenesis, lipid metabolism, cellular and tissue morphology and development, inflammation and metabolic exchange. DEG included, FST, SLC45A2, MMP19, FADS1 and GALNT6.

CONCLUSIONS

This study highlights, some of the molecular mechanisms potentially controlling uterine endometrial function during the mid-luteal phase of the estrous cycle, which may contribute to uterine endometrial mediated impaired fertility in cattle. Differentially expressed genes are potential candidate genes for the identification of genetic variation influencing cow fertility, which may be incorporated into future breeding programmes.

Effect of dietary n-3 polyunsaturated fatty acids on transcription factor regulation in the bovine endometrium

Dietary n-3 polyunsaturated fatty acid (n-3 PUFA) supplementation is postulated to have positive effects on fertility. The impact of dietary n-3 PUFA supplementation on physiological and biochemical processes involved in reproduction is likely to be associated with significant alterations in gene expression in key reproductive tissues which is in turn regulated by transcription factors. Beef heifers were supplemented with a rumen protected source of either a saturated fatty acid or high n-3 PUFA diet per animal per day for 45 days and uterine endometrial tissue was harvested post slaughter. A microarray analysis was conducted and bioinformatic tools were employed to evaluate the effect of n-3 PUFA supplementation on gene expression in the bovine endometrium. Clustering of microarray gene expression data was performed to identify co-expressed genes. Functional annotation of each cluster of genes was carried out using Ingenuity Pathway Analysis. Furthermore, oPOSSUM was employed to identify transcription factors involved in gene expression changes due to supplementary PUFA. Gene functions which showed a significant response to n-3 PUFA supplementation included tissue development, immune function and reproductive function. Numerous transcription factors such as FOXD1, FOXD3, NFKB1, ESR1, PGR, FOXA2, NKX3-1 and PPARα were identified as potential regulators of gene expression in the endometrium of cattle supplemented with n-3 PUFA. This study demonstrates the complex nature of the alterations in the transcriptional regulation process in the uterine endometrium of cattle following dietary supplementation which may positively influence the uterine environment.