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

An across breed, diet and tissue analysis reveals the transcription factor NR1H3 as a key mediator of residual feed intake in beef cattle

BACKGROUND

Provision of feed is a major determinant of overall profitability in beef production systems, accounting for up to 75% of the variable costs. Thus, improving cattle feed efficiency, by way of determining the underlying genomic control and subsequently selecting for feed efficient cattle, provides a method through which feed input costs may be reduced. The objective of this study was to undertake gene co-expression network analysis using RNA-Sequence data generated from Longissimus dorsi and liver tissue samples collected from steers of two contrasting breeds (Charolais and Holstein-Friesian) divergent for residual feed intake (RFI), across two consecutive distinct dietary phases (zero-grazed grass and high-concentrate). Categories including differentially expressed genes (DEGs) based on the contrasts of RFI phenotype, breed and dietary source, as well as key transcription factors and proteins secreted in plasma were utilised as nodes of the gene co-expression network.

RESULTS

Of the 2,929 DEGs within the network analysis, 1,604 were reported to have statistically significant correlations (≥ 0.80), resulting in a total of 43,876 significant connections between genes. Pathway analysis of clusters of co-expressed genes revealed enrichment of processes related to lipid metabolism (fatty acid biosynthesis, fatty acid β-oxidation, cholesterol biosynthesis), immune function, (complement cascade, coagulation system, acute phase response signalling), and energy production (oxidative phosphorylation, mitochondrial L-carnitine shuttle pathway) based on genes related to RFI, breed and dietary source contrasts.

CONCLUSIONS

Although similar biological processes were evident across the three factors examined, no one gene node was evident across RFI, breed and diet contrasts in both liver and muscle tissues. However within the liver tissue, the IRX4, NR1H3, HOXA13 and ZNF648 gene nodes, which all encode transcription factors displayed significant connections across the RFI, diet and breed comparisons, indicating a role for these transcription factors towards the RFI phenotype irrespective of diet and breed. Moreover, the NR1H3 gene encodes a protein secreted into plasma from the hepatocytes of the liver, highlighting the potential for this gene to be explored as a robust biomarker for the RFI trait in beef cattle.

Copy Number Variation Discovery in South African Nguni-Sired and Bonsmara-Sired Crossbred Cattle

Crossbreeding forms part of Climate-Smart beef production and is one of the strategies to mitigate the effects of climate change. Two Nguni-sired and three Bonsmara-sired crossbred animals underwent whole genome sequencing. Following quality control and file preparation, the sequence data were investigated for genome-wide copy number variation (CNV) using the panelcn.MOPS tool. A total of 355 CNVs were identified in the crossbreds, of which 274 were unique in Bonsmara-sired crossbreds and 81 unique in the Nguni-sired crossbreds. Genes that differed in copy number in both crossbreds included genes related to growth (SCRN2, LOC109572916) and fertility-related factors (RPS28, LOC1098562432, LOC109570037). Genes that were present only in the Bonsmara-sired crossbreds included genes relating to lipid metabolism (MAF1), olfaction (LOC109569114), body size (HES7), immunity (LOC10957335, LOC109877039) and disease (DMBT1). Genes that were present only in the Nguni-sired crossbreds included genes relating to ketosis (HMBOX1) and amino acid transport (LOC109572916). Results of this study indicate that Nguni and Bonsmara cattle can be utilized in crossbreeding programs as they may enhance the presence of economically important traits associated with both breeds. This will produce crossbred animals that are good meat producers, grow faster, have high fertility, strong immunity and a better chance of producing in South Africa's harsh climate conditions. Ultimately, this study provides new genetic insights into the adaptability of Nguni and Bonsmara crossbred cattle.

Effect of breed and diet on the M. longissimus thoracis et lumborum transcriptome of steers divergent for residual feed intake

Improving cattle feed efficiency through selection of residual feed intake (RFI) is a widely accepted approach to sustainable beef production. A greater understanding of the molecular control of RFI in various breeds offered contrasting diets is necessary for the accurate identification of feed efficient animals and will underpin accelerated genetic improvement of the trait. The aim of this study was to determine genes and biological processes contributing to RFI across varying breed type and dietary sources in skeletal muscle tissue. Residual feed intake was calculated in Charolais and Holstein-Friesian steers across multiple dietary phases (phase-1: high concentrate (growing-phase); phase-2: zero-grazed grass (growing-phase); phase-3: high concentrate (finishing-phase). Steers divergent for RFI within each breed and dietary phase were selected for muscle biopsy collection, and muscle samples subsequently subjected to RNAseq analysis. No gene was consistently differentially expressed across the breed and diet types examined. However, pathway analysis revealed commonality across breeds and diets for biological processes including fatty acid metabolism, immune function, energy production and muscle growth. Overall, the lack of commonality of individual genes towards variation in RFI both within the current study and compared to the published literature, suggests other genomic features warrant further evaluation in relation to RFI.

Enhancing early life nutrition alters the hepatic transcriptome of Angus × Holstein-Friesian heifer calves

Early life nutrition has a major influence on subsequent lifetime performance in cattle. The aim of this experiment was to investigate the effect of plane of nutrition from 3 to 21 weeks of age on the liver transcriptome. Holstein-Friesian × Angus heifer calves with a mean (±SD) age and BW of 19 (±5) days and 51.2 (±7.8) kg, respectively, were assigned to either a high-energy diet to support a mean average daily gain (ADG) of 1.2 kg/day (HI; n = 15) or a moderate diet (MOD; n = 15) to support a mean ADG of 0.5 kg/day. At 145 ± 3 days of age, all calves were euthanised, liver tissue samples collected and flash-frozen in liquid nitrogen. Following RNA sequence analysis, the total number of differentially expressed genes (DEGs) (at false discovery rate (FDR) > 0.05) was 537; 308 upregulated and 229 downregulated in HI compared to MOD. The number of DEGs mapped to IPA (at FDR > 0.05) was 460; 264 upregulated and 196 downregulated. There was greater expression of genes associated with cellular development and metabolism in heifers on the HI compared to the MOD diet. The genes (fold change) of the somatotrophic axis; IGF1 (3.7), IGFALS (2.6) and GHR (1.5) were upregulated in the HI compared to MOD diet. The cytokine receptor genes, IL17RB (1.7) and IL20RA (3.3), were upregulated in the HI heifers, which were detected in a network interacting with metabolically regulated genes. The potential enhanced cell-to-cell communication evident from DEGs would increase the calves' ability to combat health challenges. The findings of this study indicate that enhancing the early life plane of nutrition in heifer calves results in the upregulation of genes that are associated with increased metabolic activity and thus metabolic capacity. Moreover, the interaction between metabolic and immune communication genes indicates that enhanced nutrition has the potential to improve the immune response in the liver which will play a central role in ensuring optimal lifetime performance.

Physiological responses and adaptations to high methane production in Japanese Black cattle

In this study, using enteric methane emissions, we investigated the metabolic characteristics of Japanese Black cattle. Their methane emissions were measured at early (age 13 months), middle (20 months), and late fattening phases (28 months). Cattle with the highest and lowest methane emissions were selected based on the residual methane emission values, and their liver transcriptome, blood metabolites, hormones, and rumen fermentation characteristics were analyzed. Blood β-hydroxybutyric acid and insulin levels were high, whereas blood amino acid levels were low in cattle with high methane emissions. Further, propionate and butyrate levels differed depending on the enteric methane emissions. Hepatic genes, such as SERPINI2, SLC7A5, ATP6, and RRAD, which were related to amino acid transport and glucose metabolism, were upregulated or downregulated during the late fattening phase. The above mentioned metabolites and liver transcriptomes could be used to evaluate enteric methanogenesis in Japanese Black cattle.

Bovine rumen epithelial miRNA-mRNA dynamics reveals post-transcriptional regulation of gene expression upon transition to high-grain feeding and phytogenic supplementation

The rumen epithelium has a pivotal role in nutrient uptake and host health. This study aimed to explore the role of microRNAs (miRNAs) in the epithelial transcriptome during diet transition from forage to high-grain feeding and the modulation through supplementation with a phytogenic feed additive. Rumen biopsies were collected from 9 ruminally-cannulated non-lactating Holstein cows fed a baseline forage diet (FD) and then transitioned to high-grain feeding (HG; 65% concentrate on a dry matter basis). Cows were randomly allocated into a control group (CON, n = 5) and a group supplemented with a phytogenic feed additive (PHY, n = 4). MiRNA and mRNA sequencing was performed in parallel and transcripts were analyzed for differential expression, pathway enrichment analysis, and miRNA-mRNA interaction networks. We identified 527 miRNAs shared by all samples of the rumen epithelium, from which, bta-miR-21-5p, bta-miR-143 and bta-miR-24-3p were the most expressed. Six miRNAs were differentially expressed between CON and PHY and 8 miRNAs between FD and HG feeding, which were mainly associated with fat metabolism. Transcriptome analysis identified 9481 differentially expressed genes (DEGs) between FD and HG, whereas PHY supplementation resulted in 5 DEGs. DEGs were mainly involved in epithelium development and morphogenesis. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways associated with tricarboxylic acid and short chain fatty acid (SCFA) metabolism were enriched in DEGs between diets. MiRNA target prediction and anti-correlation analysis was used to construct networks and identify DEGs targeted by DE miRNAs responsive to diet or PHY. This study allowed the identification of potential miRNA regulation mechanisms of gene expression during transition from FD to HG feeding and phytogenic supplementation, evidencing a direct role of miRNAs in host responses to nutrition.

Identification and characterization of circular RNAs in association with the feed efficiency in Hu lambs

Background

Circular RNA (circRNA), as a new members of noncoding RNA family, have vital functions in many biological processes by as microRNA sponges or competing endogenous RNAs (ceRNAs). However, little has been reported about the genetic mechanism of circRNAs regulation of feed efficiency in sheep.

Results

This study aimed to explore the expression of circRNAs in the liver of Hu sheep with High-RFI (High residual feed intake) and Low-RFI (Low residual feed intake) using transcriptome sequencing. A total of 20,729 circRNAs were identified in two groups, in which 219 circRNAs were found as significantly differentially expressed. Several circRNAs were validated by using RT-PCR, sanger sequencing and RT-qPCR methods. These results demonstrated that the RNA-seq result and expression level of circRNAs identified are reliable. Subsequently, GO and KEGG enrichment analysis of the parental genes of the differentially expressed (DE) circRNAs were mainly involved in immunity response and metabolic process. Finally, the ceRNA regulatory networks analysis showed that the target binding sites for miRNA such as novel_41, novel_115, novel_171 and oar-miR-485-3p in the identified DE cirRNAs. Importantly, two metabolic (SHISA3 and PLEKHH2) and four (RTP4, CD274, OAS1, and RFC3) immune-related target mRNAs were identified from 4 miRNAs. Association analysis showed that the polymorphism (RTP4 c.399 A > G) in the target gene RTP4 were significantly associated with RFI (P < 0.05).

Conclusions

Analysis of sequencing data showed some candidate ceRNAs that may play key roles in the feed efficiency in sheep by regulating animal immune and metabolic. These results provide the basis data for further study of the biological functions of circRNAs in regulating sheep feed efficiency.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08517-5.

Carcass Quality Profiles and Associated Genomic Regions of South African Goat Populations Investigated Using Goat SNP50K Genotypes

Simple Summary

South Africa is one of the major goat producing countries on the African continent. However, despite a large number of goats being produced, there is still a growing demand for chevon, which leads to producers being unable to reach demand, resulting in an absence of chevon in retail markets. Carcass quality is an important economic trait that plays a major role in influencing consumer preferences and high nutrient provision. Even though chevon is an easily accessible meat for smallholder farmers and has health benefits, it is still less preferred due to perceptions of low meat quality attributes such as toughness, off-odours and flavour, and unappealing colour. The majority of goat populations are village ecotypes whose genetic potential for meat and carcass quality is unknown.

Abstract

Carcass quality includes a battery of essential economic meat traits that play a significant role in influencing farmer breed preferences. A preliminary study was undertaken to investigate the carcass quality and the associated genomic regions in a small nucleus of animals that are representative of South African goat genetic resources. Samples of the South African Boer (n = 14), Northern Cape Speckled (n = 14), Eastern Cape Xhosa Lob ear (n = 12), Nguni/Mbuzi (n = 13), and Village (n = 20) were genotyped using the Illumina goat SNP50K and were phenotyped for carcass quality traits. SA Boer goats had heavier warm and cold carcass weights (17.2 ± 2.3 kg and 16.3 ± 2.3 kg). Pella village goats raised under an intensive system had significantly (p < 0.05) heavier warm and cold carcass weights (9.9 ± 1.1 kg and 9.2 ± 1.2 kg) compared to the village goats that are raised extensively (9.1 ± 2.0 kg and 8.4 ± 1.9). A total of 40 SNPs located on chromosomes 6, 10, 12, 13, 19, and 21 were significantly associated with carcass traits at (−log10 [p < 0.05]). Candidate genes that were associated with carcass characteristics (GADD45G, IGF2R, GAS1, VAV3, CAPN8, CAPN7, CAPN2, GHSR, COLQ, MRAS, and POU1F1) were also observed. Results from this study will inform larger future studies that will ultimately find use in breed improvement programs.

Chemical Group-Based Metabolome Analysis Identifies Candidate Plasma Biomarkers Associated With Residual Feed Intake in Beef Steers

We applied chemical group-based metabolomics to identify blood metabolic signatures associated with residual feed intake in beef cattle. A group of 56 crossbred growing beef steers (average BW = 261.3 ± 18.5 kg) were adapted to a high-forage total mixed ration in a confinement dry lot equipped with GrowSafe intake nodes for period of 49 d to determine their residual feed intake classification (RFI). After RFI determination, weekly blood samples were collected three times from beef steers with the lowest RFI [most efficient (HFE); n = 8] and highest RFI and least-efficient [least efficient (LFE); n = 8]. Plasma was prepared by centrifugation and composited for each steer. Metabolome analysis was conducted using a chemical isotope labeling (CIL)/liquid chromatography–mass spectrometry, which permitted the analysis of metabolites containing amine/phenol-, carboxylic acid-, and carbonyl-chemical groups, which are metabolites associated with metabolisms of amino acids, fatty acids, and carbohydrates, respectively. A total number of 495 amine/phenol-containing metabolites were detected and identified; pathway analysis of all these metabolites showed that arginine biosynthesis and histidine metabolism were enriched (P &lt; 0.10) in HFE, relative to LFE steers. Biomarker analyses of the amine/phenol-metabolites identified methionine, 5-aminopentanoic acid, 2-aminohexanedioic acid, and 4-chlorolysine as candidate biomarkers of RFI [false discovery rate ≤ 0.05; Area Under the Curve (AUC) &gt; 0.90]. A total of 118 and 330 metabolites containing carbonyl- and carboxylic acid-chemical groups, respectively were detected and identified; no metabolic pathways associated with these metabolites were altered and only one candidate biomarker (methionine sulfoxide) was identified. These results identified five candidate metabolite biomarkers of RFI in beef cattle which are mostly associated with amino acid metabolism. Further validation using a larger cohort of beef cattle of different genetic pedigree is required to confirm these findings.

Identifying the key genes and functional enrichment pathways associated with feed efficiency in cattle

Residual feed intake (RFI) is a measurement of feed efficiency, and is inversely correlated with feed efficiency. The differentially expressed genes (DEGs) associated with RFI vary substantially among studies, posing great challenges in finding the RFI-related marker genes. This study attempted to resolve this issue by integrating and comparing the multiple transcriptome sequencing data associated with RFI in the cattle liver, using differential, functional enrichment, protein-protein interaction (PPI) network, weighted co-expression network (WGCNA), and gene set enrichment analyses (GSEA) to identify the candidate genes and functional enrichment pathways that are closely associated with RFI. Four candidate genes namely SHC1, GPX4, ACADL, and IGF1 were identified and validated as the marker genes for RFI. Four functional enrichment pathways, namely the fatty acid metabolism, sugar metabolism, energy metabolism, and protein ubiquitination were also found to be closely related to RFI. This study identified several genes and signaling pathways with shared characteristics, which will provide new insights into the molecular mechanisms related to the regulation of feed efficiency, and provide basis for molecular markers related to feed efficiency in beef cattle.

Identification of candidate genes and enriched biological functions for feed efficiency traits by integrating plasma metabolites and imputed whole genome sequence variants in beef cattle

BACKGROUND

Feed efficiency is one of the key determinants of beef industry profitability and sustainability. However, the cellular and molecular background behind feed efficiency is largely unknown. This study combines imputed whole genome DNA variants and 31 plasma metabolites to dissect genes and biological functions/processes that are associated with residual feed intake (RFI) and its component traits including daily dry matter intake (DMI), average daily gain (ADG), and metabolic body weight (MWT) in beef cattle.

RESULTS

Regression analyses between feed efficiency traits and plasma metabolites in a population of 493 crossbred beef cattle identified 5 (L-valine, lysine, L-tyrosine, L-isoleucine, and L-leucine), 4 (lysine, L-lactic acid, L-tyrosine, and choline), 1 (citric acid), and 4 (L-glutamine, glycine, citric acid, and dimethyl sulfone) plasma metabolites associated with RFI, DMI, ADG, and MWT (P-value < 0.1), respectively. Combining the results of metabolome-genome wide association studies using 10,488,742 imputed SNPs, 40, 66, 15, and 40 unique candidate genes were identified as associated with RFI, DMI, ADG, and MWT (P-value < 1 × 10-5), respectively. These candidate genes were found to be involved in some key metabolic processes including metabolism of lipids, molecular transportation, cellular function and maintenance, cell morphology and biochemistry of small molecules.

CONCLUSIONS

This study identified metabolites, candidate genes and enriched biological functions/processes associated with RFI and its component traits through the integrative analyses of metabolites with phenotypic traits and DNA variants. Our findings could enhance the understanding of biochemical mechanisms of feed efficiency traits and could lead to improvement of genomic prediction accuracy via incorporating metabolite data.

Characterization and Duodenal Transcriptome Analysis of Chinese Beef Cattle With Divergent Feed Efficiency Using RNA-Seq

Residual feed intake (RFI) is an important measure of feed efficiency for agricultural animals. Factors associated with cattle RFI include physiology, dietary factors, and the environment. However, a precise genetic mechanism underlying cattle RFI variations in duodenal tissue is currently unavailable. The present study aimed to identify the key genes and functional pathways contributing to variance in cattle RFI phenotypes using RNA sequencing (RNA-seq). Six bulls with extremely high or low RFIs were selected for detecting differentially expressed genes (DEGs) by RNA-seq, followed by conducting GO, KEGG enrichment, protein-protein interaction (PPI), and co-expression network (WGCNA, n = 10) analysis. A total of 380 differentially expressed genes was obtained from high and low RFI groups, including genes related to energy metabolism (ALDOA, HADHB, INPPL1), mitochondrial function (NDUFS1, RFN4, CUL1), and feed intake behavior (CCK). Two key sub-networks and 26 key genes were detected using GO analysis of DEGs and PPI analysis, such as TPM1 and TPM2, which are involved in mitochondrial pathways and protein synthesis. Through WGCNA, a gene network was built, and genes were sorted into 27 modules, among which the blue (r = 0.72, p = 0.03) and salmon modules (r = -0.87, p = 0.002) were most closely related with RFI. DEGs and genes from the main sub-networks and closely related modules were largely involved in metabolism; oxidative phosphorylation; glucagon, ribosome, and N-glycan biosynthesis, and the MAPK and PI3K-Akt signaling pathways. Through WGCNA, five key genes, including FN1 and TPM2, associated with the biological regulation of oxidative processes and skeletal muscle development were identified. Taken together, our data suggest that the duodenum has specific biological functions in regulating feed intake. Our findings provide broad-scale perspectives for identifying potential pathways and key genes involved in the regulation of feed efficiency in beef cattle.

Identification of Predictor Genes for Feed Efficiency in Beef Cattle by Applying Machine Learning Methods to Multi-Tissue Transcriptome Data

Machine learning (ML) methods have shown promising results in identifying genes when applied to large transcriptome datasets. However, no attempt has been made to compare the performance of combining different ML methods together in the prediction of high feed efficiency (HFE) and low feed efficiency (LFE) animals. In this study, using RNA sequencing data of five tissues (adrenal gland, hypothalamus, liver, skeletal muscle, and pituitary) from nine HFE and nine LFE Nellore bulls, we evaluated the prediction accuracies of five analytical methods in classifying FE animals. These included two conventional methods for differential gene expression (DGE) analysis (t-test and edgeR) as benchmarks, and three ML methods: Random Forests (RFs), Extreme Gradient Boosting (XGBoost), and combination of both RF and XGBoost (RX). Utility of a subset of candidate genes selected from each method for classification of FE animals was assessed by support vector machine (SVM). Among all methods, the smallest subsets of genes (117) identified by RX outperformed those chosen by t-test, edgeR, RF, or XGBoost in classification accuracy of animals. Gene co-expression network analysis confirmed the interactivity existing among these genes and their relevance within the network related to their prediction ranking based on ML. The results demonstrate a great potential for applying a combination of ML methods to large transcriptome datasets to identify biologically important genes for accurately classifying FE animals.

Effects of Intensive Fattening With Total Mixed Rations on Carcass Characteristics, Meat Quality, and Meat Chemical Composition of Yak and Mechanism Based on Serum and Transcriptomic Profiles

The objective of this study was to investigate the effects of intensive fattening with total mixed rations (TMR) on carcass characteristics, meat quality, and chemical composition of the yak meat. Theoretical data has been provided for evaluating the quality of yak meat during natural grazing and short-term fattening. Based on the analysis, we found that in fattening yak, the carcass weight (CWT) was increased by 106.43%, whereas the cooking loss, tenderness, and drop loss were significantly improved due to higher intramuscular fat content and lower moisture (P < 0.05). Protein, fat, calcium, and amino acids were also much higher (P < 0.01) in fattening yak compared with the grazing yak. The levels of albumin (ALB), lactate dehydrogenase (LDH), triglyceride (TRIG), and amylase (AMYL) in serum indicated better nutritional status for fattening yaks. The transcriptomics analysis showed that the high expression of ACSL1 and ACACB genes improved the synthesis and deposition of fat in fattening yak, whereas the regulation of SLC7A8, ATP1A4, ATP1A1, SLC3A2, and CPA3 gene expression weakened the proteolysis. These results indicated that fattening with TMR improves the yield and quality of the yak meat.

Weighted Single-Step GWAS Identified Candidate Genes Associated with Growth Traits in a Duroc Pig Population

Growth traits are important economic traits of pigs that are controlled by several major genes and multiple minor genes. To better understand the genetic architecture of growth traits, we performed a weighted single-step genome-wide association study (wssGWAS) to identify genomic regions and candidate genes that are associated with days to 100 kg (AGE), average daily gain (ADG), backfat thickness (BF) and lean meat percentage (LMP) in a Duroc pig population. In this study, 3945 individuals with phenotypic and genealogical information, of which 2084 pigs were genotyped with a 50 K single-nucleotide polymorphism (SNP) array, were used for association analyses. We found that the most significant regions explained 2.56–3.07% of genetic variance for four traits, and the detected significant regions (>1%) explained 17.07%, 18.59%, 23.87% and 21.94% for four traits. Finally, 21 genes that have been reported to be associated with metabolism, bone growth, and fat deposition were treated as candidate genes for growth traits in pigs. Moreover, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses implied that the identified genes took part in bone formation, the immune system, and digestion. In conclusion, such full use of phenotypic, genotypic, and genealogical information will accelerate the genetic improvement of growth traits in pigs.

Identification of functional candidate variants and genes for feed efficiency in Holstein and Jersey cattle breeds using RNA-sequencing

The identification of functional genetic variants and associated candidate genes linked to feed efficiency may help improve selection for feed efficiency in dairy cattle, providing economic and environmental benefits for the dairy industry. This study used RNA-sequencing data obtained from liver tissue from 9 Holstein cows [n = 5 low residual feed intake (RFI), n = 4 high RFI] and 10 Jersey cows (n = 5 low RFI, n = 5 high RFI), which were selected from a single population of 200 animals. Using RNA-sequencing, 3 analyses were performed to identify: (1) variants within low or high RFI Holstein cattle; (2) variants within low or high RFI Jersey cattle; and (3) variants within low or high RFI groups, which are common across both Holstein and Jersey cattle breeds. From each analysis, all variants were filtered for moderate, modifier, or high functional effect, and co-localized quantitative trait loci (QTL) classes, enriched biological processes, and co-localized genes related to these variants, were identified. The overlapping of the resulting genes co-localized with functional SNP from each analysis in both breeds for low or high RFI groups were compared. For the first two analyses, the total number of candidate genes associated with moderate, modifier, or high functional effect variants fixed within low or high RFI groups were 2,810 and 3,390 for Holstein and Jersey breeds, respectively. The major QTL classes co-localized with these variants included milk and reproduction QTL for the Holstein breed, and milk, production, and reproduction QTL for the Jersey breed. For the third analysis, the common variants across both Holstein and Jersey breeds, uniquely fixed within low or high RFI groups were identified, revealing a total of 86,209 and 111,126 functional variants in low and high RFI groups, respectively. Across all 3 analyses for low and high RFI cattle, 12 and 31 co-localized genes were overlapping, respectively. Among the overlapping genes across breeds, 9 were commonly detected in both the low and high RFI groups (INSRR, CSK, DYNC1H1, GAB1, KAT2B, RXRA, SHC1, TRRAP, PIK3CB), which are known to play a key role in the regulation of biological processes that have high metabolic demand and are related to cell growth and regeneration, metabolism, and immune function. The genes identified and their associated functional variants may serve as candidate genetic markers and can be implemented into breeding programs to help improve the selection for feed efficiency in dairy cattle.

Bovine hepatic miRNAome profiling and differential miRNA expression analyses between beef steers with divergent feed efficiency phenotypes

MicroRNAs (miRNAs) are small RNA molecules involved in regulation of multiple biological processes through modulating expression of their target genes. Here we employed RNAseq to profile liver tissue miRNAome of 60 steers from Angus, Charolais, and Kinsella Composite (KC) populations. Of these animals, 36 animals (n = 12 for each breed) were utilized to identify differentially expressed (DE) miRNAs between animals with high (n = 6) or low (n = 6) phenotypic values of residual feed intake (RFI), a common measurement of feed efficiency. At a threshold of fold-change > 1.5 and P-value < 0.05, we detected 12 (7 up- and 5 downregulated in low-RFI animals), 18 (12 up- and 6 downregulated), and 13 (8 up- and 5 downregulated) DE miRNAs for Angus, Charolais, and KC steers, respectively. Most of the DE miRNAs were breed specific, with bta-miR-449a and bta-miR-AB-2 being differentially expressed in all three breeds. The predicted target genes of the identified DE miRNA are mainly involved in cell cycle, cell death and survival, cell signaling, cellular growth and proliferation, protein trafficking, cell morphology, cell-to-cell signaling and interaction, cellular development, molecular transport, post-translational modification, as well as nutrient metabolism (lipids, carbohydrates, protein and amino acid). Our results provide insights into the bovine hepatic miRNAome and their potential roles in molecular regulation of RFI in beef cattle.

Genes associated with body weight gain and feed intake identified by meta-analysis of the mesenteric fat from crossbred beef steers

Mesenteric fat is a visceral fat depot that increases with cattle maturity and can be influenced by diet. There may be a relationship between the accumulation of mesenteric fat and feed efficiency in beef cattle. The purpose of this study was to identify genes that may be differentially expressed in steers with high and low BW gain and feed intake. RNA-Seq was used to evaluate the transcript abundance of genes in the mesenteric fat from a total of 78 steers collected over 5 different cohorts. A meta-analysis was used to identify genes involved with gain, feed intake or the interaction of both phenotypes. The interaction analysis identified 11 genes as differentially expressed. For the main effect of gain, a total of 87 differentially expressed genes (DEG) were identified (P_ADJ<0.05), and 24 were identified in the analysis for feed intake. Genes identified for gain were involved in functions and pathways including lipid metabolism, stress response/protein folding, cell proliferation/growth, axon guidance and inflammation. The genes for feed intake did not cluster into pathways, but some of the DEG for intake had functions related to inflammation, immunity, and/or signal transduction (JCHAIN, RIPK1, LY86, SPP1, LYZ, CD5, CD53, SRPX, and NF2). At P_ADJ<0.1, only 4 genes (OLFML3, LOC100300716, MRPL15, and PUS10) were identified as differentially expressed in two or more cohorts, highlighting the importance of evaluating the transcriptome of more than one group of animals and incorporating a meta-analysis. This meta-analysis has produced many mesenteric fat DEG that may be contributing to gain and feed intake in cattle.

Transcriptomic analyses suggest a dominant role of insulin in the coordinated control of energy metabolism and ureagenesis in goat liver

BACKGROUND

The ureagenesis plays a central role in the homeostatic control of nitrogen metabolism. This process occurs in the liver, the key metabolic organ in the maintenance of energy homeostasis in the body. To date, the understanding of the influencing factors and regulators of ureagenesis in ruminants is still poor. The aim of this study was to investigate the relationship between energy metabolism and ureagenesis and detect the direct regulators of ureagenesis in the liver by using RNA-seq technology.

RESULTS

Eighteen four-month-old male goats were divided into two groups randomly and received a diet containing 10% (LNFC group, n = 9) or 30% non-fiber carbohydrate (MNFC group, n = 9), respectively, for four weeks. The global gene expression analysis of liver samples showed that, compared with a LNFC diet, the MNFC diet promoted the expression of genes required for synthesis of fatty acid and glycerol, whereas it suppressed those related to fatty acid oxidation, gluconeogenesis from amino acids and ureagenesis. Additionally, gene expression for rate-limiting enzymes of ureagenesis were highly correlated to the gene expression of key enzymes of both fatty acid synthesis and glycerol synthesis (Spearman correlation coefficient > 0.8 and p < 0.05). In the differentially expressed signaling pathways related to the endocrine system, the MNFC diet activated the insulin and PPAR signaling pathway, whereas it suppressed the leptin-JAK/STAT signaling pathway, compared with the LNFC diet. Reverse transcription quantitative PCR analyses of 40 differentially expressed genes confirmed the RNA-seq results (R² = 0.78).

CONCLUSION

Our study indicated that a dietary NFC-induced increase of energy supply promoted lipid anabolism and decreased ureagenesis in the caprine liver. By combining our results with previously published reports, insulin signaling can be suggested to play the dominant role in the coordinated control of hepatic energy metabolism and ureagenesis.