Global liver gene expression differences in Nelore steers with divergent residual feed intake phenotypes

Differentially expressed messenger RNA isoforms in beef cattle skeletal muscle with different fatty acid profiles

This study aimed to identify mRNA isoforms that were expressed differently in the muscle tissue of Nellore cattle based on their intramuscular fatty acid profile. Forty-eight young bulls were used to quantify beef fatty acids (FA) and perform RNA sequencing analysis. The young bulls were divided into three different groups based on quantifying FA using k-means analysis. The Grp1 clustered animals with significantly elevated levels of PUFA, ω6, ω3, linoleic acid, α-linolenic acid, and PUFA/SFA ratios, indicating a more favorable fatty acid profile. Animals in Group 3 demonstrated significantly higher levels of palmitic acid, stearic acid, myristic acid, and SFA, which are less favorable fatty acid profiles. Grp2 included bulls with intermediate levels of fatty acids, positioned between the profiles of Grp1 and Grp3. Differential expression (DE) analyses were performed to compare these three distinct groups through the contrasts: Grp1 vs. Grp2, Grp1 vs. Grp3, and Grp2 vs. Grp3. The DE analyses identified a total of 62, 26, and 30 transcripts for the contrasts Grp1 vs. Grp2, Grp1 vs. Grp3, and Grp2 vs. Grp3, respectively. In the comparison between the Grp1 and Grp2 groups, we identified three mRNA isoforms, C7-203, ADD1-204, and OXT-201, which are involved in glycogen and lipid metabolism. These mRNA isoforms regulate the key genes responsible for fatty acid synthesis, leading to a higher PUFA content profile. On the other hand, in the comparison between the Grp1 and Grp3 groups, the mRNA isoforms RBM3-202 and TRAG1-202 were identified and play a crucial role in muscle development, adipogenesis, and concentration of PUFA and MUFA, respectively. The downregulation of THRSP-201 and FABP4-201, isoforms identified in both, Grp1 vs. Grp2 and Grp2 vs. Grp3, contrasts may contribute to lower levels of MUFA and SFA and higher levels of PUFA. In addition, these mRNA isoforms were associated with lipogenesis, fatty acid transport, and inhibition of lipolysis. Our findings suggest that the identified mRNA isoforms could serve as promising candidates to help develop strategies to select animals of higher nutritional meat quality.

Melatonin alleviates endometrial fibrosis in bovine endometritis by regulating TGF-β/Smad and MAPK signaling pathways via MT2

Bovine endometritis can lead to abnormal endometrial function and fibrosis, resulting in difficulties in successful embryo implantation and intrauterine adhesions. Melatonin is well known for its profitable effects against inflammation and pathological fibrosis in discrepant organs. Considering the potential therapeutic benefits of melatonin, this study aimed to investigate its effects on endometrial fibrosis in cows with endometritis. Firstly, we evaluated the expression patterns of various factors associated with fibrosis, such as transforming growth factor-β1 (TGF-β1), extracellular matrix (ECM)-related markers (COL1A1 and COL3A1), epithelial-mesenchymal transformation (EMT)-related proteins (α-SMA and Vimentin) in healthy and endometritis-affected bovine uterine tissues. The results showed that diseased tissues presented significantly higher TGF-β1 expression, ECM production, and EMT progression versus normal tissues. Moreover, we established an LPS-induced fibrosis model in endometrial stromal cells (ESCs), and found that melatonin inhibited the fibrosis process of ESCs in a dose-dependent manner. The MT2 inhibitor 4P-PDOT blocked the antifibrotic effects of melatonin and inhibited the phosphorylation of Smad2/3, ERK1/2, and JNK1/2 in LPS-induced fibrosis of ESCs, but not P38 MAPK. These data implied that melatonin supplementation attenuated LPS-induced fibrosis in ESCs by modulating the inhibition of TGF-β/Smad, ERK, and JNK signaling pathways via MT2.

Identification of biological pathways and putative candidate genes for residual feed intake in a tropically adapted beef cattle breed by plasma proteome analysis

This study identified potential biomarkers for feed efficiency by blood plasma proteome analysis of a tropically adapted beef cattle breed. Two experimental groups were selected based on residual feed intake (RFI). The proteome was investigated by LC-MS/MS in a data-dependent acquisition mode. After quality control, 123 differentially abundant proteins (DAPs) were identified between the two experimental groups. Among DAPs with the highest absolute log-fold change values, the PRDM2, KRT5, UGGT1, DENND5B, B2M, SLC44A2, SLC7A2, PTPRC, and FETUB were highlighted as potential biomarkers because of their functions that may contribute to RFI. Furthermore, functional enrichment analysis revealed several biological processes, molecular functions and pathways that contributes to RFI, such as cell signaling, cellular responses to stimuli, immune system, calcium, hormones, metabolism and functions of proteins, lipids and carbohydrates. Protein-protein interaction analysis identified 32 and 11 DAPs as important nodes based on their interactions in the high- and low-RFI groups, respectively. This study represents the first comprehensive profiling of the blood plasma proteome of a tropically adapted beef cattle breed and provides valuable insights into the potential roles of these DAPs in key biological processes and pathways, contributing to our understanding of the mechanisms underlying feed efficiency in tropically adapted beef cattle. SIGNIFICANCE: LC-MS/MS analysis was performed to investigate changes in the blood plasma proteome associated with residual feed intake (RFI) in a tropically adapted beef cattle breed (Bos taurus taurus). Some putative biomarkers were identified to distinguish the high-RFI to low-RFI animals, based on their log-fold change value or on their protein-protein interaction network, which provide helpful sources in developing novel selection strategies for breeding programs. Our findings also revealed valuable insights into the metabolic pathways and biological processes that contribute to RFI in beef cattle, such as those closely linked to cell signaling, cellular responses to stimuli, immune system, calcium, hormones, metabolism and functions of proteins, lipids and carbohydrates.

Differential gene expression in neonatal calf muscle tissues from Hanwoo cows overfed during mid to late pregnancy period

Maternal nutrition significantly influences fetal development and postnatal outcomes. This study investigates the impact of maternal overfeeding during mid to late pregnancy on gene expression in the round and sirloin muscles of Hanwoo neonatal calves. Eight cows were assigned to either a control group receiving standard nutrition (100%) or a treated group receiving overnutrition (150%). After birth, tissue samples from the round and sirloin muscles of neonatal calves were collected and subjected to RNA sequencing to assess differentially expressed genes (DEGs). RNA sequencing identified 43 DEGs in round muscle and 15 in sirloin muscle, involving genes related to myogenesis, adipogenesis, and energy regulation. Key genes, including PPARGC1A, THBS1, CD44, JUND, CNN1, ENAH, and RUNX1, were predominantly downregulated. Gene ontology (GO) enrichment analyses revealed terms associated with muscle development, such as "biological regulation," "cellular process," and "response to stimulus." Protein-protein interaction networks highlighted complex interactions among DEGs. Random Forest analysis identified ARC, SLC1A5, and GNPTAB as influential genes for distinguishing between control and treated groups. Overall, maternal overnutrition during mid-to-late pregnancy results in the downregulation of genes involved in muscle development and energy metabolism in neonatal Hanwoo calves. These findings provide insights into the molecular effects of maternal nutrition on muscle development.

Effect of breed and dietary composition on the miRNA profile of beef steers divergent for feed efficiency

Identifying and breeding cattle that are more feed efficient is of great benefit to beef production. Additionally, it is crucial that genes contributing to feed efficiency are robust across varying management settings including dietary source as well as being relevant across contrasting breeds of cattle. The aim of this study was to determine miRNAs that are contributing to the expression of residual feed intake (RFI) across two breeds and dietary sources. miRNA profiling was undertaken in Longissimus dorsi tissue of Charolais and Holstein-Friesian steers divergent for RFI phenotype following two contrasting consecutive diets (high-forage and high-concentrate). Ten miRNA were identified as differentially expressed (adj. P < 0.1) across the breed and diet contrasts examined. Of particular interest was the differential expression of miR-2419-5p and miR-2415-3p, both of which were up-regulated in the Low-RFI Charolais steers across each dietary phase. Pathway analysis of target mRNA genes of differentially expressed miRNA revealed enrichment (P < 0.05) for pathways including metabolic related pathways, insulin receptor signalling, adipogenesis as well as pathways related to skeletal muscle growth. These results provide insight into the skeletal muscle miRNAome of beef cattle and their potential molecular regulatory mechanisms relating to feed efficiency in beef cattle.

Unveiling Genetic Potential for Equine Meat Production: A Bioinformatics Approach

In view of the predicted significant increase in global meat production, alternative sources such as horsemeat are becoming increasingly important due to their lower environmental impact and high nutritional value. This study aimed to identify SNP markers on the GeneSeek® Genomic Profiler™ Equine (Neogen, Lansing, MI, USA) that are important for horsemeat production traits. First, orthologous genes related to meat yield in cattle and common genes between horses and cattle within QTLs for body size and weight were identified. Markers for these genes were then evaluated based on predicted variant consequences, GERP scores, and positions within constrained elements and orthologous regulatory regions in pigs. A total of 268 markers in 57 genes related to meat production were analyzed. This resulted in 27 prioritized SNP markers in 22 genes, including notable markers in LCORL, LASP1, IGF1R, and MSTN. These results will benefit smallholder farmers by providing genetic insights for selective breeding that could improve meat yield. This study also supports future large-scale genetic analyses such as GWAS and Genomic Best Linear Unbiased Prediction (GBLUP). The results of this study may be helpful in improving the accuracy of genomic breeding values. However, limitations include reliance on bioinformatics without experimental validation. Future research can validate these markers and consider a wider range of traits to ensure accuracy in equine breeding.

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.

A network-based approach to understanding gene-biological processes affecting economically important traits of Nelore cattle

This study aimed to build gene-biological process networks with differentially expressed genes associated with economically important traits of Nelore cattle from 17 previous studies. The genes were clustered into three groups by evaluated traits: group 1, production traits; group 2, carcass traits; and group 3, meat quality traits. For each group, a gene-biological process network analysis was performed with the differentially expressed genes in common. For production traits, 37 genes were found in common, of which 13 genes were enriched for six Gene Ontology (GO) terms; these terms were not functionally grouped. However, the enriched GO terms were related to homeostasis, the development of muscles and the immune system. For carcass traits, four genes were found in common. Thus, it was not possible to functionally group these genes into a network. For meat quality traits, the analysis revealed 222 genes in common. CSRP3 was the only gene differentially expressed in all three groups. Non-redundant biological terms for clusters of genes were functionally grouped networks, reflecting the cross-talk between all biological processes and genes involved. Many biological processes and pathways related to muscles, the immune system and lipid metabolism were enriched, such as striated muscle cell development and triglyceride metabolic processes. This study provides insights into the genetic mechanisms of production, carcass and meat quality traits of Nelore cattle. This information is fundamental for a better understanding of the complex traits and could help in planning strategies for the production and selection systems of Nelore cattle.

Genome-wide association study for reproduction traits in Colombian Creole Blanco Orejinegro cattle

The profitability of the beef cattle production system relies heavily on reproductive traits. Unfortunately, certain traits, such as age at first calving (AFC), calving interval (CI), and gestation length (GL), can pose challenges in traditional breeding programs because of their low heritability (0.01-0.12) and sex-limited characteristics. Another important aspect is the conservation of the genetic resources of animals adapted to the Colombian regions, which implies the preservation and rational use of the creole breeds in the country market. Therefore, this study aimed to identify genomic regions in the creole cattle breed Blanco Orejinegro (BON) that influence the reproductive traits in females. The dataset comprised 439 animals and 118,116 single-nucleotide polymorphisms' (SNPs) markers. The GS3 program was used to identify the SNP effects employing the BAYES Cπ methodology. The number of SNPs with effect for AFC was 25, 1527 for CI, and 23 for GL. Some of the genes found associated with reproductive and growth traits as well as immune response and environmental adaptation ECE1, EPH, EPHB2, SMARCAL1, IGFBP5, IGFBP2, FCGRT, EGFR, MUL1, PINK1, STPG1, CNGB1, TGFB1, OXTR, IL22RA1, MYOM3, OXTR, CNR2, HIVEP3, CTPS1, CXCL8, FCGRT, MREG, TMEM169, PECR, and MC1R. Our results evidenced a high contribution of the genetic architecture of the Colombian creole cattle breed Blanco Orejinegro that may impact should be included in implementing genetic improvement and conservation programs.

Transcriptomic analysis of yak longissimus dorsi muscle identifies genes associated with tenderness

Meat tenderness is an important sensory index when consumers choose meat products, which determines the value of meat products and consumers' buying intentions. Yak meat is rich in nutrition and unique in flavor, which is favored by consumers. However, its meat has the deficiencies of low tenderness and poor taste, which has a negative impact on the value of its meat products and customer satisfaction. To identify the genes affecting the yak meat tenderness, we used RNA-seq to analyze the longissimus dorsi muscle of yaks with different tenderness, screened a total of 1120 differentially expressed genes (DEGs). Meanwhile, 23 pathways were significantly enriched. By further analysis, we identified eight genes related to yak meat tenderness (WNT5A, ARID5B, SERPINE1 KLHL40, RUNX1, MAFF, RFX7 and ARID5A). Notably, SERPINE1 was involved in the significant enrichment pathways of 'complement and coagulation cascade pathway', 'HIF-1 signaling pathway' and 'AGE-RAGE signaling pathway in diabetic complications' which can regulate meat tenderness. This implies that SERPINE1 may play an important regulatory role among them. The DEGs associated with yak meat quality screened in this work will be helpful to identify potential biomarkers related to meat tenderness.

Analysis of liver miRNA in Hu sheep with different residual feed intake

Feed efficiency (FE), an important economic trait in sheep production, is indirectly assessed by residual feed intake (RFI). However, RFI in sheep is varied, and the molecular processes that regulate RFI are unclear. It is thus vital to investigate the molecular mechanism of RFI to developing a feed-efficient sheep. The miRNA-sequencing (RNA-Seq) was utilized to investigate miRNAs in liver tissue of 6 out of 137 sheep with extreme RFI phenotypic values. In these animals, as a typical metric of FE, RFI was used to distinguish differentially expressed miRNAs (DE_miRNAs) between animals with high (n = 3) and low (n = 3) phenotypic values. A total of 247 miRNAs were discovered in sheep, with four differentially expressed miRNAs (DE_miRNAs) detected. Among these DE_miRNAs, three were found to be upregulated and one was downregulated in animals with low residual feed intake (Low_RFI) compared to those with high residual feed intake (High_RFI). The target genes of DE_miRNAs were primarily associated with metabolic processes and biosynthetic process regulation. Furthermore, they were also considerably enriched in the FE related to glycolysis, protein synthesis and degradation, and amino acid biosynthesis pathways. Six genes were identified by co-expression analysis of DE_miRNAs target with DE_mRNAs. These results provide a theoretical basis for us to understand the sheep liver miRNAs in RFI molecular regulation.

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.

Combined transcriptome and metabolome analysis reveal key regulatory genes and pathways of feed conversion efficiency of oriental river prawn Macrobrachium nipponense

BACKGROUND

Oriental river prawn Macrobrachium nipponense is an economically important aquaculture species in China, Japan, and Vietnam. In commercial prawn farming, feed cost constitutes about 50 to 65% of the actual variable cost. Improving feed conversion efficiency in prawn culture will not only increase economic benefit, but also save food and protect the environment. The common indicators used for feed conversion efficiency include feed conversion ratio (FCR), feed efficiency ratio (FER), and residual feed intake (RFI). Among these, RFI is much more suitable than FCR and FER during the genetic improvement of feed conversion efficiency for aquaculture species.

RESULTS

In this study, the transcriptome and metabolome of hepatopancreas and muscle of M. nipponense from high RFI low RFI groups, which identified after culture for 75 days, were characterized using combined transcriptomic and metabolomic analysis. A total of 4540 differentially expressed genes (DEGs) in hepatopancreas, and 3894 DEGs in muscle were identified, respectively. The DEGs in hepatopancreas were mainly enriched in KEGG pathways including the metabolism of xenobiotics by cytochrome P450 (down-regulated), fat digestion and absorption (down-regulated) and aminoacyl-tRNA biosynthesis (up-regulated), etc. The DEGs in muscle were mainly enriched in KEGG pathways including the protein digestion and absorption (down-regulated), glycolysis/gluconeogenesis (down-regulated), and glutathione metabolism (up-regulated), etc. At the transcriptome level, the RFI of M. nipponense was mainly controlled in biological pathways such as the high immune expression and the reduction of nutrients absorption capacity. A total of 445 and 247 differently expressed metabolites (DEMs) were identified in the hepatopancreas and muscle, respectively. At the metabolome level, the RFI of M. nipponense was affected considerably by amino acid and lipid metabolism.

CONCLUSIONS

M. nipponense from higher and lower RFI groups have various physiological and metabolic capability processes. The down-regulated genes, such as carboxypeptidase A1, 6-phosphofructokinase, long-chain-acyl-CoA dehydrogenase, et. al., in digestion and absorption of nutrients, and the up-regulated metabolites, such as aspirin, lysine, et. al., in response to immunity could be potential candidate factors contributed to RFI variation for M. nipponense. Overall, these results would provide new insights into the molecular mechanism of feed conversion efficiency and assist in selective breeding to improve feed conversion efficiency in M. nipponense.

Association of residual feed intake with intestinal microbiome and metabolome in laying period of ducks

Introduction

Residual feed intake (RFI) is a indicator to evaluate animal feed. This experiment was explored to study the relationship between intestinal microbiome and metabolome of ducks with different residual feed intake during laying period.

Methods

A total of 300 Shaoxing ducks aged 42 weeks were randomly selected and fed a diet of 60 d. At the end of the trial, 20 samples were selected according to the phenotype of RFI and divided into two groups (HRFI and LRFI). The cecal microbiota composition was explored by 16S ribosomal RNA gene sequencing and rectal metabolomics uses liquid chromatography-mass spectrometry (LC-MS) to identify the composition of metabolites in a non-targeted manner.

Results

Results show feed intake and feed conversion ratio in the group HRFI were significantly higher than those in the group LRFI (p < 0.05). Chao1 indices were higher in the group LRFI than in the HRFI (p < 0.05), Shannon and Simpson indices were higher in the group LRFI than in the HRFI (p < 0.01). After linear discriminant analysis effect size (p < 0.05, LDA score > 3), Rikenellaceae, Rikenellaceae_RC9_gut_group, Lactobacillales and Ruminococcus_2, etc. were significantly enriched in the group LRFI at the genus level, while Prevotellaceae_NK3B31_group and Bacteria were significantly enriched in the group HRFI. After LC-MS analysis we found 338 metabolic difference products and 10 metabolic pathways, including the ABC transporter system, cysteine and methionine metabolism, arginine and proline metabolism, and vitamin B6 metabolism, were identified to be associated with the significantly differentially expressed between the groups LRFI and HRFI (p < 0.05). We hypothesize that the difference between ducks with different RFIs is mainly due to the fact that ducks with LRFI have more SCFAs-producing bacteria in their gut microorganisms, which regulate the RFI of animals. This process we found that Phascolarctobaterium and Anaerobiospirillum may provide energy for ABC transporter system by producing SCFAs, and regulate RFI to improve feed utilization efficiency.

Discussion

These results revealed the relationship between microbiome and metabonomics in laying ducks with different RFI, and provided theoretical basis for further study on the relationship between them.

Feed efficiency in dairy sheep: An insight from the milk transcriptome

Introduction

As higher feed efficiency in dairy ruminants means a higher capability to transform feed nutrients into milk and milk components, differences in feed efficiency are expected to be partly linked to changes in the physiology of the mammary glands. Therefore, this study aimed to determine the biological functions and key regulatory genes associated with feed efficiency in dairy sheep using the milk somatic cell transcriptome.

Material and methods

RNA-Seq data from high (H-FE, n = 8) and low (L-FE, n = 8) feed efficiency ewes were compared through differential expression analysis (DEA) and sparse Partial Least Square-Discriminant analysis (sPLS-DA).

Results

In the DEA, 79 genes were identified as differentially expressed between both conditions, while the sPLS-DA identified 261 predictive genes [variable importance in projection (VIP) &gt; 2] that discriminated H-FE and L-FE sheep.

Discussion

The DEA between sheep with divergent feed efficiency allowed the identification of genes associated with the immune system and stress in L-FE animals. In addition, the sPLS-DA approach revealed the importance of genes involved in cell division (e.g., KIF4A and PRC1) and cellular lipid metabolic process (e.g., LPL, SCD, GPAM, and ACOX3) for the H-FE sheep in the lactating mammary gland transcriptome. A set of discriminant genes, commonly identified by the two statistical approaches, was also detected, including some involved in cell proliferation (e.g., SESN2, KIF20A, or TOP2A) or encoding heat-shock proteins (HSPB1). These results provide novel insights into the biological basis of feed efficiency in dairy sheep, highlighting the informative potential of the mammary gland transcriptome as a target tissue and revealing the usefulness of combining univariate and multivariate analysis approaches to elucidate the molecular mechanisms controlling complex traits.

Transcriptome Profiling of the Liver in Nellore Cattle Phenotypically Divergent for RFI in Two Genetic Groups

The identification and selection of genetically superior animals for residual feed intake (RFI) could enhance productivity and minimize environmental impacts. The aim of this study was to use RNA-seq data to identify the differentially expressed genes (DEGs), known non-coding RNAs (ncRNAs), specific biomarkers and enriched biological processes associated with RFI of the liver in Nellore cattle in two genetic groups. In genetic group 1 (G1), 24 extreme RFI animals (12 low RFI (LRFI) versus 12 high RFI (HRFI)) were selected from a population of 60 Nellore bulls. The RNA-seq of the samples from their liver tissues was performed using an Illumina HiSeq 2000. In genetic group 2 (G2), 20 samples of liver tissue of Nellore bulls divergent for RFI (LRFI, n = 10 versus HRFI, n = 10) were selected from 83 animals. The raw data of the G2 were chosen from the ENA repository. A total of 1811 DEGs were found for the G1 and 2054 for the G2 (p-value ≤ 0.05). We detected 88 common genes in both genetic groups, of which 33 were involved in the immune response and in blocking oxidative stress. In addition, seven (B2M, ADSS, SNX2, TUBA4A, ARHGAP18, MECR, and ABCF3) possible gene biomarkers were identified through a receiver operating characteristic analysis (ROC) considering an AUC > 0.70. The B2M gene was overexpressed in the LRFI group. This gene regulates the lipid metabolism protein turnover and inhibits cell death. We also found non-coding RNAs in both groups. MIR25 was up-regulated and SNORD16 was down-regulated in the LRFI for G1. For G2, up-regulated RNase_MRP and SCARNA10 were found. We highlight MIR25 as being able to act by blocking cytotoxicity and oxidative stress and RMRP as a blocker of mitochondrial damage. The biological pathways associated with RFI of the liver in Nellore cattle in the two genetic groups were for energy metabolism, protein turnover, redox homeostasis and the immune response. The common transcripts, biomarkers and metabolic pathways found in the two genetic groups make this unprecedented work even more relevant, since the results are valid for different herds raised in different ways. The results reinforce the biological importance of these known processes but also reveal new insights into the complexity of the liver tissue transcriptome of Nellore cattle.

Transcript Characteristics on the Susceptibility Difference of Bovine Respiratory Disease

Bovine respiratory disease (BRD) is one of the major health issues in the cattle industry, resulting in significant financial crises globally. There is currently no good treatment, and cattle are made resistant to pneumonia through disease-resistant breeding. The serial blood samples from six Xinjiang brown (XJB) calves were collected for the RNA sequencing (RNA-seq). The obtained six samples were grouped into two groups, in each group as infected with BRD and healthy calves, respectively. In our study, the differential expression mRNAs were detected by using RNA-seq and constructed a protein-protein interaction (PPI) network related to the immunity in cattle. The key genes were identified by protein interaction network analysis, and the results from RNA-seq were verified using reverse transcription-quantitative polymerase chain reaction (RT-qPCR). A total of 488 differentially expressed (DE) mRNAs were identified. Importantly, the enrichment analysis of these identified DEGs classified them as mainly enriched in the regulation and immune response processes. The 16 hub genes were found to be related to immune pathways categorized by PPIs analysis. Results revealed that many hub genes were related to the immune response to respiratory disease. These results will provide the basis for a better understanding of the molecular mechanism of bovine resistance to BRD.

Multi-transcriptomics reveals RLMF axis-mediated signaling molecules associated with bovine feed efficiency

The regulatory axis plays a vital role in interpreting the information exchange and interactions among mammal organs. In this study on feed efficiency, it was hypothesized that a rumen-liver-muscle-fat (RLMF) regulatory axis exists and scrutinized the flow of energy along the RLMF axis employing consensus network analysis from a spatial transcriptomic standpoint. Based on enrichment analysis and protein-protein interaction analysis of the consensus network and tissue-specific genes, it was discovered that carbohydrate metabolism, energy metabolism, immune and inflammatory responses were likely to be the biological processes that contribute most to feed efficiency variation on the RLMF regulatory axis. In addition, clusters of genes related to the electron respiratory chain, including ND (2,3,4,4L,5,6), NDUF (A13, A7, S6, B3, B6), COX (1,3), CYTB, UQCR11, ATP (6,8), clusters of genes related to fatty acid metabolism including APO (A1, A2, A4, B, C3), ALB, FG (A, G), as well as clusters of the ribosomal-related gene including RPL (8,18A,18,15,13, P1), the RPS (23,27A,3A,4X), and the PSM (A1-A7, B6, C1, C3, D2-D4, D8 D9, E1) could be the primary effector genes responsible for feed efficiency variation. The findings demonstrate that high feed efficiency cattle, through the synergistic action of the regulatory axis RLMF, may improve the efficiency of biological processes (carbohydrate metabolism, protein ubiquitination, and energy metabolism). Meanwhile, high feed efficiency cattle might enhance the ability to respond to immunity and inflammation, allowing nutrients to be efficiently distributed across these organs associated with digestion and absorption, energy-producing, and energy-storing organs. Elucidating the distribution of nutrients on the RLMF regulatory axis could facilitate an understanding of feed efficiency variation and achieve the study on its molecular regulation.

Detection of potential functional variants based on systems-biology: the case of feed efficiency in beef cattle

BACKGROUND

Potential functional variants (PFVs) can be defined as genetic variants responsible for a given phenotype. Ultimately, these are the best DNA markers for animal breeding and selection, especially for polygenic and complex phenotypes. Herein, we described the identification of PFVs for complex phenotypes (in this case, Feed Efficiency in beef cattle) using a systems-biology driven approach based on RNA-seq data from physiologically relevant organs.

RESULTS

The systems-biology coupled with deep molecular phenotyping by RNA-seq of liver, muscle, hypothalamus, pituitary, and adrenal glands of animals with high and low feed efficiency (FE) measured by residual feed intake (RFI) identified 2,000,936 uniquely variants. Among them, 9986 variants were significantly associated with FE and only 78 had a high impact on protein expression and were considered as PFVs. A set of 169 significant uniquely variants were expressed in all five organs, however, only 27 variants had a moderate impact and none of them a had high impact on protein expression. These results provide evidence of tissue-specific effects of high-impact PFVs. The PFVs were enriched (FDR < 0.05) for processing and presentation of MHC Class I and II mediated antigens, which are an important part of the adaptive immune response. The experimental validation of these PFVs was demonstrated by the increased prediction accuracy for RFI using the weighted G matrix (ssGBLUP+wG; Acc = 0.10 and b = 0.48) obtained in the ssGWAS in comparison to the unweighted G matrix (ssGBLUP; Acc = 0.29 and b = 1.10).

CONCLUSION

Here we identified PFVs for FE in beef cattle using a strategy based on systems-biology and deep molecular phenotyping. This approach has great potential to be used in genetic prediction programs, especially for polygenic phenotypes.

Prune homolog 2 with BCH domain (PRUNE2) gene expression is associated with feed efficiency-related traits in Nelore steers

Animal feeding is a critical factor in increasing producer profitability. Improving feed efficiency can help reduce feeding costs and reduce the environmental impact of beef production. Candidate genes previously identified for this trait in differential gene expression studies (e.g., case-control studies) have not examined continuous gene-phenotype variation, which is a limitation. The aim of this study was to investigate the association between the expression of five candidate genes in the liver, measured by quantitative real-time PCR and feed-related traits. We adopted a linear mixed model to associate liver gene expression from 52 Nelore steers with the following production traits: average daily gain (ADG), body weight (BW), dry matter intake (DMI), feed conversion ratio (FCR), feed efficiency (FE), Kleiber index (KI), metabolic body weight (MBW), residual feed intake (RFI), and relative growth ratio (RGR). The total expression of the prune homolog 2 (PRUNE2) gene was significantly associated with DMI, FCR, FE, and RFI (P < 0.05). Furthermore, we have identified a new transcript of PRUNE2 (TCONS_00027692, GenBank MZ041267) that was inversely correlated with FCR and FE (P < 0.05), in contrast to the originally identified PRUNE2 transcript. The cytochrome P450 subfamily 2B (CYP2B6), early growth response protein 1 (EGR1), collagen type I alpha 1 chain (COL1A1), and connective tissue growth factor (CTGF) genes were not associated with any feed efficiency-related traits (P > 0.05). The findings reported herein suggest that PRUNE2 expression levels affects feed efficiency-related traits variation in Nelore steers.

EEF1A1 transcription cofactor gene polymorphism is associated with muscle gene expression and residual feed intake in Nelore cattle

Cis-acting effects of noncoding variants on gene expression and regulatory molecules constitute a significant factor for phenotypic variation in complex traits. To provide new insights into the impacts of single-nucleotide polymorphisms (SNPs) on transcription factors (TFs) and transcription cofactors (TcoF) coding genes, we carried out a multi-omic analysis to identify cis-regulatory effects of SNPs on these genes' expression in muscle and describe their association with feed efficiency-related traits in Nelore cattle. As a result, we identified one SNP, the rs137256008C > T, predicted to impact the EEF1A1 gene expression (β = 3.02; P-value = 3.51E-03) and the residual feed intake trait (β = - 3.47; P-value = 0.02). This SNP was predicted to modify transcription factor sites and overlaps with several QTL for feed efficiency traits. In addition, co-expression network analyses showed that animals containing the T allele of the rs137256008 SNP may be triggering changes in the gene network. Therefore, our analyses reinforce and contribute to a better understanding of the biological mechanisms underlying gene expression control of feed efficiency traits in bovines. The cis-regulatory SNP can be used as biomarker for feed efficiency in Nelore cattle.

Residual Feed Intake in Beef Cattle Is Associated With Differences in Hepatic mRNA Expression of Fatty Acid, Amino Acid, and Mitochondrial Energy Metabolism Genes

We evaluated the mRNA expression of genes involved in hepatic fatty acid, amino acid, and mitochondrial energy metabolism in crossbred beef steers with divergent low and high residual feed intake (RFI). Low-RFI beef steers (n = 8; RFI = - 1.93 kg/d) and high-RFI beef steers (n = 8; RFI = + 2.01kg/d) were selected from a group of 56 growing crossbred beef steers (average BW = 261 ± 18.5 kg) fed a high-forage total mixed ration after a 49-d performance testing period. At the end of the 49-d performance testing period, liver biopsies were collected from the low-RFI and high-RFI beef steers for RNA extraction and cDNA synthesis. The mRNA expression of 84 genes each related to fatty acid metabolism, amino acid metabolism, and mitochondrial energy metabolism were analyzed using pathway-focused PCR-based arrays. The mRNA expression of 8 genes (CRAT, SLC27A5, SLC27A2, ACSBG2, ACADL, ACADSB, ACAA1, and ACAA2) involved fatty acid transport and β-oxidation were upregulated (FC ≥ 2.0, FDR ≤ 0.05) in low-RFI, compared to high-RFI steers. Among those involved in amino acid metabolism, hepatic mRNA expression of a gene encoding for aminoadipate aminotransferase, an enzyme related to lysine degradation, was downregulated (FC = -5.45, FDR = 0.01) in low-RFI steers, whereas those of methionine adenosyltransferase I and aspartate aminotransferase 2, which both link amino acid and lipid metabolism, were upregulated (FC ≥ 2.0, FDR ≤ 0.05). Two mitochondrial energy metabolism genes (UQCRC1 and ATP5G1) involved in ATP synthesis via oxidative phosphorylation were upregulated (FC ≥ 2.0, FDR ≤ 0.05) in low-RFI beef steers, compared to high-RFI beef steers. The results of this study demonstrated that low-RFI beef steers exhibit upregulation of molecular mechanisms related to fatty acid transport, fatty acid β-oxidation, and mitochondrial ATP synthesis, which suggest that low-RFI beef steers have enhanced metabolic capacity to maximize capture of energy and nutrients from feeds consumed.

Proteomic analysis of adipose tissue revealing differentially abundant proteins in highly efficient mid-lactating dairy cows

The improvement of nutrient utilization efficiency in dairy cows represents an important task in view of the current rising demand for animal products and sustainable resource usage. In this perspective, the identification of appropriate markers to identify the most efficient animals for dairy production becomes a crucial factor. Residual feed intake (RFI), which represents the difference between predicted and actual intake, is used to define the efficiency of cows. In this study, subcutaneous adipose tissue (AT) was collected from five high efficient (HEF) and five low efficient (LEF) mid-lactation Holstein dairy cows, that represented subgroups of the 20% lowest RFI values (HEF) and highest 20% RFI values (LEF), out of a cohort of 155 cows that were examined for feed efficiency at the individual dairy barn at Volcani Institute, Israel. Adipose samples were examined for proteomic analysis by nano-LC/MS–MS and gene expression by RT-PCR. A total of 101 differential proteins (P ≤ 0.05 and fold change ± 1.5) and two protein networks related to feed efficiency were found between HEF and LEF cows. Among the enriched top canonical pathways, FAT10 signaling, EIF2 signaling, Sirtuin signaling, Acute phase response signaling, Protein ubiquitination and mTOR signaling pathways were related to feed efficiency in AT. Furthermore, abundance of transferrin (TF; FC = 78.35, P = 0.02) enriched pathways, including mTOR signaling, LXR/RXR and FXR/RXR activation was found in AT of HEF cows. Relative mRNA expression of RBM39, which is involved in energy metabolism, was decreased in AT of HEF versus LEF. The relationship found between the AT proteins and/or metabolic pathways and the feed efficiency demonstrates that AT may reflect metabolic adaptations to high efficiency, and suggests that these proteins together with their metabolic mechanisms are suitable candidates as biomarkers to identify efficient cows for dairy production.

Comparison of Growth Performance, Immunity, Antioxidant Capacity, and Liver Transcriptome of Calves between Whole Milk and Plant Protein-Based Milk Replacer under the Same Energy and Protein Levels

High-cost milk proteins necessitate cheaper, effective milk replacer alternatives, such as plant proteins. To examine plant protein-based milk replacer’s impact on growth performance, serum immune and antioxidant indicators, and liver transcriptome profiles in suckling calves. We assigned 28 newborn Holstein calves (41.60 ± 3.67 kg of body weight at birth) to milk (M) or milk replacer (MR) and starter diets pre-weaning (0–70 d of age) but with the same starter diet post-weaning (71–98 d of age). During the pre-weaning period, compared with the M group, MR group had significantly lower body weight, withers height, heart girth, average daily gain, feed efficiency, serum immunoglobulin (Ig) M concentration, superoxide dismutase concentration, and total antioxidant capacity; whereas they had significantly higher serum aspartate aminotransferase concentration. During the post-weaning period, MR group presented significantly higher average daily gain, alanine transaminase, aspartate aminotransferase, and malonaldehyde concentrations; whereas they had significantly lower serum IgA and IgM concentrations than the M group. Transcriptome analysis revealed 1, 120 and 293 differentially expressed genes (DEGs; MR vs. M group) in the calves from pre- and post-weaning periods, respectively. The DEGs related to xenobiotic and lipid metabolism and those related to energy metabolism, immune function, and mineral metabolism were up- and downregulated, respectively, during the pre-weaning period; during the post-weaning period, the DEGs related to osteoclast differentiation and metabolic pathways showed difference. In this study, compared with M group, MR group had the same growth performance during the overall experimental period; however, MR affected the hepatic metabolism, immune, and antioxidant function of calves. These observations can facilitate future studies on milk replacers.

Unraveling Admixture, Inbreeding, and Recent Selection Signatures in West African Indigenous Cattle Populations in Benin

The Dwarf Lagune and the Savannah Somba cattle in Benin are typical representatives of the endangered West African indigenous Shorthorn taurine. The Lagune was previously exported to African and European countries and bred as Dahomey cattle, whereas the Somba contributed to the formation of two indigenous hybrids known as Borgou and Pabli cattle. These breeds are affected by demographic, economic, and environmental pressures in local production systems. Considering current and historical genomic data, we applied a formal test of admixture, estimated admixture proportions, and computed genomic inbreeding coefficients to characterize the five breeds. Subsequently, we unraveled the most recent selection signatures using the cross-population extended haplotype homozygosity approach, based on the current and historical genotypes. Results from principal component analyses and high proportion of Lagune ancestry confirm the Lagune origin of the European Dahomey cattle. Moreover, the Dahomey cattle displayed neither indicine nor European taurine (EUT) background, but they shared on average 40% of autozygosity from common ancestors, dated approximately eight generations ago. The Lagune cattle presented inbreeding coefficients larger than 0.13; however, the Somba and the hybrids (Borgou and Pabli) were less inbred (≤0.08). We detected evidence of admixture in the Somba and Lagune cattle, but they exhibited a similar African taurine (AFT) ancestral proportion (≥96%) to historical populations, respectively. A moderate and stable AFT ancestral proportion (62%) was also inferred for less admixed hybrid cattle including the Pabli. In contrast, the current Borgou samples displayed a lower AFT ancestral proportion (47%) than historical samples (63%). Irrespective of the admixture proportions, the hybrid populations displayed more selection signatures related to economic traits (reproduction, growth, and milk) than the taurine. In contrast, the taurine, especially the Somba, presented several regions known to be associated with adaptive traits (immunity and feed efficiency). The identified subregion of bovine leukocyte antigen (BoLA) class IIb (including DSB and BOLA-DYA) in Somba cattle is interestingly uncommon in other African breeds, suggesting further investigations to understand its association with specific adaptation to endemic diseases in Benin. Overall, our study provides deeper insights into recent evolutionary processes in the Beninese indigenous cattle and their aptitude for conservation and genetic improvement.

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.

Genome-Wide Association Study Identifies Candidate Genes Associated with Feet and Leg Conformation Traits in Chinese Holstein Cattle

Simple Summary

Feet and leg problems are among the major reasons for dairy cows leaving the herd, as well as having direct association with production and reproduction efficiency, health (e.g., claw disorders and lameness) and welfare. Hence, understanding the genetic architecture underlying feet and conformation traits in dairy cattle offers new opportunities toward the genetic improvement and long-term selection. Through a genome-wide association study on Chinese Holstein cattle, we identified several candidate genes associated with feet and leg conformation traits. These results could provide useful information about the molecular breeding basis of feet and leg traits, thus improving the longevity and productivity of dairy cattle.

Abstract

Feet and leg conformation traits are considered one of the most important economical traits in dairy cattle and have a great impact on the profitability of milk production. Therefore, identifying the single nucleotide polymorphisms (SNPs), genes and pathways analysis associated with these traits might contribute to the genomic selection and long-term plan selection for dairy cattle. We conducted genome-wide association studies (GWASs) using the fixed and random model circulating probability unification (FarmCPU) method to identify SNPs associated with bone quality, heel depth, rear leg side view and rear leg rear view of Chinese Holstein cows. Phenotypic measurements were collected from 1000 individuals of Chinese Holstein cattle and the GeneSeek Genomic Profiler Bovine 100 K SNP chip was utilized for individual genotyping. After quality control, 984 individual cows and 84,906 SNPs remained for GWAS work; as a result, we identified 20 significant SNPs after Bonferroni correction. Several candidate genes were identified within distances of 200 kb upstream or downstream to the significant SNPs, including ADIPOR2, INPP4A, DNMT3A, ALDH1A2, PCDH7, XKR4 and CADPS. Further bioinformatics analyses showed 34 gene ontology terms and two signaling pathways were significantly enriched (p ≤ 0.05). Many terms and pathways are related to biological quality, metabolism and development processes; these identified SNPs and genes could provide useful information about the genetic architecture of feet and leg traits, thus improving the longevity and productivity of Chinese Holstein dairy cattle.

Differences in net fat oxidation, heat production, and liver mitochondrial DNA copy numbers between high and low feed-efficient dairy cows

Improving feed utilization efficiency in dairy cattle could have positive economic and environmental effects that would support the sustainability of the dairy industry. Identifying key differences in metabolism between high and low feed-efficient animals is vital to enhancing feed conversion efficiency. Therefore, our objectives were (1) to determine whether cows grouped by either high or low feed efficiency have measurable differences in net fat and carbohydrate metabolism that account for differences in heat production (HP), and if so, whether these differences also exists under conditions of feed withdrawal when the effect of feeding on HP is minimized, and (2) to determine whether the abundance of mitochondria in the liver can be related to the high or low feed-efficient groups. Ten dairy cows from a herd of 15 (parity = 2) were retrospectively grouped into either a high (H) or a low (L) feed-efficient group (n = 5 per group) based on weekly energy-corrected milk (ECM) divided by dry mater intake (DMI) from wk 4 through 30 of lactation. Livers were biopsied at wk -4, 2, and 12, and blood was sampled weekly from wk -3 to 12 relative to parturition. Blood was subset to be analyzed for the transition period (wk -3 to 3) and from wk 4 to 12. In wk 5.70 ± 0.82 (mean ± SD) postpartum (PP), cows spent 2 d in respiration chambers (RC), in which CO₂, O₂, and CH₄ gases were measured every 6 min for 24 h. Fatty acid oxidation (FOX), carbohydrate oxidation (COX), metabolic respiratory quotient (RQ), and HP were calculated from gas measurements for 23 h. Cows were fed ad libitum (AD-LIB) on d 1 and had feed withdrawn (RES, restricted diet) on d 2. Additional blood samples were taken at the end of the AD-LIB and RES feeding periods in the RC. During wk 4 to 30 PP, H had greater DMI/kg of metabolic body weight (BW^0.75), ECM per kilogram of BW^0.75 yield, and ECM/DMI ratio, compared with L, but a lower body condition score between wk 4 and 12 PP. In the RC period, we detected no differences in BW, DMI, or milk yield between groups. We also detected no significant group or group by feeding period interactions for plasma metabolites except for Revised Quantitative Insulin Sensitivity Check Index, which tended to have a group by feeding period interaction. The H group had lower HP and HP per kilogram of BW^0.75 compared with L. Additionally, H had lower FOX and FOX per kilogram of BW^0.75 compared with L during the AD-LIB period. Methane, CH₄ per kilogram of BW^0.75, and CH₄ per kilogram of milk yield were lower in H compared with L, but, when adjusted for DMI, CH₄/DMI did not differ between groups, nor did HP/DMI. Relative mitochondrial DNA copy numbers in the liver were lower in the L than in the H group. These results suggest that lower feed efficiency in dairy cows may result from fewer mitochondria per liver cell as well as a greater whole-body HP, which likely partially results from higher net fat oxidation.

An examination of skeletal muscle and hepatic tissue transcriptomes from beef cattle divergent for residual feed intake

The selection of cattle with enhanced feed efficiency is of importance with regard to reducing feed costs in the beef industry. Global transcriptome profiling was undertaken on liver and skeletal muscle biopsies from Simmental heifers and bulls divergent for residual feed intake (RFI), a widely acknowledged feed efficiency phenotype, in order to identify genes that may be associated with this trait. We identified 5 genes (adj. p < 0.1) to be differentially expressed in skeletal muscle between high and low RFI heifers with all transcripts involved in oxidative phosphorylation and mitochondrial homeostasis. A total of 11 genes (adj. p < 0. 1) were differentially expressed in liver tissue between high and low RFI bulls with differentially expressed genes related to amino and nucleotide metabolism as well as endoplasmic reticulum protein processing. No genes were identified as differentially expressed in either heifer liver or bull muscle analyses. Results from this study show that the molecular control of RFI in young cattle is modified according to gender, which may be attributable to differences in physiological maturity between heifers and bulls of the same age. Despite this we have highlighted a number of genes that may hold potential as molecular biomarkers for RFI cattle.

Characterization of Breed Specific Differences in Spermatozoal Transcriptomes of Sheep in Australia

Reduced reproductive efficiency results in economic losses to the Australian sheep industry. Reproductive success, particularly after artificial insemination, is dependent on a number of contributing factors on both ewe and ram sides. Despite considerable emphasis placed on characterising ewe side contributions, little emphasis has been placed on characterising ram side contributions to conception success. Over 14,000 transcripts are in spermatozoa of other species, which are transferred to the ova on fertilisation. These transcripts conceivably influence early embryonic development and whether conception is successful. Semen was collected (n = 45) across three breeds; Merino, Dohne, and Poll Dorset. Following collection, each ejaculate was split in two; an aliquot was assessed utilising Computer Assisted Semen Analysis (CASA) and the remaining was utilised for RNA extraction and subsequent next-generation sequencing. Overall, 754 differentially expressed genes were identified in breed contrasts and contrast between ejaculates of different quality. Downstream analysis indicated that these genes could play significant roles in a broad range of physiological functions, including maintenance of spermatogenesis, fertilisation, conception, embryonic development, and offspring production performance. Overall results provide evidence that the spermatozoal transcriptome could be a crucial contributing factor in improving reproductive performance as well as in the overall productivity and profitability of sheep industries.

Effect of dietary restriction and subsequent realimentation on hepatic oxidative phosphorylation in cattle

Compensatory growth (CG) is a naturally accelerated growth which occurs upon realimentation, following a prior period of dietary restriction. The process is harnessed worldwide as a management practice to reduce feed costs in beef cattle production. The objective of this study was to assess the potential contribution of hepatic cellular mitochondrial capacity to CG through global hepatic oxidative phosphorylation gene expression analyses as well as functional mitochondrial enzyme activity assays. Holstein-Friesian bulls were separated into two groups: (i) restricted feed allowance for 125 days (Period 1) (RES; n = 30) followed by ad-libitum feeding for 55 days (Period 2) or (ii) ad-libitum access to feed throughout (Periods 1 and 2) (ADLIB; n = 30). At the end of each period, 15 animals from each treatment group were slaughtered and hepatic tissue was collected. Tissue samples were subjected to RNAseq and spectrophotometric analysis for the functional assessment of mitochondria. RES and ADLIB groups grew at 0.6 kg/day and 1.9 kg/day, respectively, during Period 1. During Period 2, the RES group underwent CG growing at 2.5 kg/day, with ADLIB animals gaining 1.4 kg/day. Oxidative phosphorylation genes were differentially expressed in response to both dietary restriction and CG. Spectrophotometric assays indicated that mitochondrial abundance was greater in animals undergoing dietary restriction at the end of Period 1 and subsequently reduced during realimentation (P < 0.02). Results indicate that mitochondrial capacity may be enhanced during dietary restriction to more effectively utilize diet-derived nutrients. However, enhanced mitochondrial capacity does not appear to be directly contributing to CG in cattle.

Genome-wide association study and pathway analysis for fat deposition traits in nellore cattle raised in pasture-based systems

Genome-wide association study (GWAS) is a powerful tool to identify candidate genes and genomic regions underlying key biological mechanisms associated with economically important traits. In this context, the aim of this study was to identify genomic regions and metabolic pathways associated with backfat thickness (BFT) and rump fat thickness (RFT) in Nellore cattle, raised in pasture-based systems. Ultrasound-based measurements of BFT and RFT (adjusted to 18 months of age) were collected in 11,750 animals, with 39,903 animals in the pedigree file. Additionally, 1,440 animals were genotyped using the GGP-indicus 35K SNP chip, containing 33,623 SNPs after the quality control. The single-step GWAS analyses were performed using the BLUPF90 family programs. Candidate genes were identified through the Ensembl database incorporated in the BioMart tool, while PANTHER and REVIGO were used to identify the key metabolic pathways and gene networks. A total of 18 genomic regions located on 10 different chromosomes and harbouring 23 candidate genes were identified for BFT. For RFT, 22 genomic regions were found on 14 chromosomes, with a total of 29 candidate genes identified. The results of the pathway analyses showed important genes for BFT, including TBL1XR1, AHCYL2, SLC4A7, AADAT, VPS53, IDH2 and ETS1, which are involved in lipid metabolism, synthesis of cellular amino acids, transport of solutes, transport between Golgi Complex membranes, cell differentiation and cellular development. The main genes identified for RFT were GSK3β, LRP1B, EXT1, GRB2, SORCS1 and SLMAP, which are involved in metabolic pathways such as glycogen synthesis, lipid transport and homeostasis, polysaccharide and carbohydrate metabolism. Polymorphisms located in these candidate genes can be incorporated in commercial genotyping platforms to improve the accuracy of imputation and genomic evaluations for carcass fatness. In addition to uncovering biological mechanisms associated with carcass quality, the key gene pathways identified can also be incorporated in biology-driven genomic prediction methods.

Transcriptome assisted label free proteomics of hepatic tissue in response to both dietary restriction and compensatory growth in cattle

Compensatory growth (CG) is a naturally occurring phenomenon where, following a period of under nutrition, an animal exhibits accelerated growth upon re-alimentation. The objective was to identify and quantify hepatic proteins involved in the regulation of CG in cattle. Forty Holstein Friesian bulls were equally assigned to one of four groups. Groups; A1 and A2 had ad libitum access to feed for 125 days, groups R1 and R2 were feed restricted. Following this, R1 and A1 animals were slaughtered. Remaining animals (R2 and A2) were slaughtered following ad libitum feeding for a successive 55 days. At slaughter hepatic tissue samples were collected and label-free quantitative proteomics undertaken with spectra searched against a custom built transcriptome database specific to the animals in this study. 24 differentially abundant proteins were identified during CG (R2 vs. R1) including; PSPH, ASNS and GSTM1, which are involved in nutrient metabolism, immune response and cellular growth. Proteins involved in biochemical pathways related to nutrient metabolism were down-regulated during CG, indicating a possible adaptive response by the liver to a period of fluctuating nutrient availability. The livers ability to regulate its metabolic activity may have profound effects on the efficiency of whole body energy utilization during CG. SIGNIFICANCE: This study is the first to unravel the effect of compensatory growth on the hepatic proteome of cattle using transcriptome-assisted shot gun proteomics. Proteins identified as being affected by dietary restriction and subsequent expression of compensatory growth in this study may, following appropriate validation, contribute to the identification of functional genetic variants. Such information could be harnessed within the context of genomic selection in cattle breeding programs to identify animals with a greater genetic potential to undergo compensatory growth, thus increasing the profitability of the beef sector and accelerating genetic gain.

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.

Factors affecting expression and transcription of uncoupling protein 2 gene

Previous studies suggest a negative relationship between hepatic oxidative stress and productivity in beef cattle. Uncoupling protein 2 (UCP2) is involved in the disappearance of reactive oxygen species, suggesting the defensive role of UCP2 against oxidative stress. The present study examined the relationship between oxidative stress and expression levels of UCP2/Ucp2 in cultured human and mouse liver-derived cells. We also explored factors regulating bovine Ucp2 transcription. As oxidative stress inducers, hydrogen peroxide, ethanol, and cumene hydroperoxide (CmHP) were used. Expression levels of hemoxygenase 1 (HMOX1), a representative gene induced by oxidative stress, were not affected by any oxidative stress inducers in HepG2 human liver-derived cells. The levels of UCP2 mRNA were also unaffected by the oxidative stress inducers. Treatment with CmHP increased expression of Hmox1 in Hepa1-6 mouse liver-derived cells, but Ucp2 expression was not changed. Stimulus screening for regulator of transcription (SSRT) revealed that expression of p50 or p65, transcription factors conferring response to oxidative stress, did not stimulate bovine Ucp2 transcrition in HepG2 cells. SSRT also showed 11 molecules that induced Ucp2 transcription more than 4-fold; among them, endoplasmic reticulum (ER) stress-related transcription factors such as XBP1, c-JUN, JUNB, and C/EBPβ were identified. However, treatment with ER stress inducers did not increase Ucp2 expression in HepG2 and Hepa1-6 cells. The present results suggest that 1) neither oxidative stress nor ER stress induces Ucp2 expression in liver-derived cells, and 2) Ucp2 transcription is stimulated by several transcription factors.

Development and comparison of RNA-sequencing pipelines for more accurate SNP identification: practical example of functional SNP detection associated with feed efficiency in Nellore beef cattle

Background

Optimization of an RNA-Sequencing (RNA-Seq) pipeline is critical to maximize power and accuracy to identify genetic variants, including SNPs, which may serve as genetic markers to select for feed efficiency, leading to economic benefits for beef production. This study used RNA-Seq data (GEO Accession ID: PRJEB7696 and PRJEB15314) from muscle and liver tissue, respectively, from 12 Nellore beef steers selected from 585 steers with residual feed intake measures (RFI; n = 6 low-RFI, n = 6 high-RFI). Three RNA-Seq pipelines were compared including multi-sample calling from i) non-merged samples; ii) merged samples by RFI group, iii) merged samples by RFI and tissue group. The RNA-Seq reads were aligned against the UMD3.1 bovine reference genome (release 94) assembly using STAR aligner. Variants were called using BCFtools and variant effect prediction (VeP) and functional annotation (ToppGene) analyses were performed.

Results

On average, total reads detected for Approach i) non-merged samples for liver and muscle, were 18,362,086.3 and 35,645,898.7, respectively. For Approach ii), merging samples by RFI group, total reads detected for each merged group was 162,030,705, and for Approach iii), merging samples by RFI group and tissues, was 324,061,410, revealing the highest read depth for Approach iii). Additionally, Approach iii) merging samples by RFI group and tissues, revealed the highest read depth per variant coverage (572.59 ± 3993.11) and encompassed the majority of localized positional genes detected by each approach. This suggests Approach iii) had optimized detection power, read depth, and accuracy of SNP calling, therefore increasing confidence of variant detection and reducing false positive detection. Approach iii) was then used to detect unique SNPs fixed within low- (12,145) and high-RFI (14,663) groups. Functional annotation of SNPs revealed positional candidate genes, for each RFI group (2886 for low-RFI, 3075 for high-RFI), which were significantly (P < 0.05) associated with immune and metabolic pathways.

Conclusion

The most optimized RNA-Seq pipeline allowed for more accurate identification of SNPs, associated positional candidate genes, and significantly associated metabolic pathways in muscle and liver tissues, providing insight on the underlying genetic architecture of feed efficiency in beef cattle.

Exploring the Regulatory Potential of Long Non-Coding RNA in Feed Efficiency of Indicine Cattle

Long non-coding RNA (lncRNA) can regulate several aspects of gene expression, being associated with complex phenotypes in humans and livestock species. In taurine beef cattle, recent evidence points to the involvement of lncRNA in feed efficiency (FE), a proxy for increased productivity and sustainability. Here, we hypothesized specific regulatory roles of lncRNA in FE of indicine cattle. Using RNA-Seq data from the liver, muscle, hypothalamus, pituitary gland and adrenal gland from Nellore bulls with divergent FE, we submitted new transcripts to a series of filters to confidently predict lncRNA. Then, we identified lncRNA that were differentially expressed (DE) and/or key regulators of FE. Finally, we explored lncRNA genomic location and interactions with miRNA and mRNA to infer potential function. We were able to identify 126 relevant lncRNA for FE in Bos indicus, some with high homology to previously identified lncRNA in Bos taurus and some possible specific regulators of FE in indicine cattle. Moreover, lncRNA identified here were linked to previously described mechanisms related to FE in hypothalamus-pituitary-adrenal axis and are expected to help elucidate this complex phenotype. This study contributes to expanding the catalogue of lncRNA, particularly in indicine cattle, and identifies candidates for further studies in animal selection and management.

Genomic Tools for the Conservation and Genetic Improvement of a Highly Fragmented Breed-The Ramo Grande Cattle from the Azores

Simple Summary

Inbreeding control is a key concern in managing local endangered breeds, which often have developed unique adaptation features. Ramo Grande is a local cattle breed raised in the Azores archipelago under very harsh conditions, with a census of about 1300 cows dispersed by various islands. This fragmentation is a challenge when the goal is to keep inbreeding under control. Currently, panels of genetic markers are available which enable the assessment of inbreeding and the occurrence of previous bottlenecks in a population. These panels also allow the identification of genes associated with specific production traits, if reliable phenotypic information is available. We used a panel of genetic markers and estimated that the degree of inbreeding was approaching a level of concern, while some exotic gene inflow may have occurred in the past. We were able to identify genetic markers significantly associated with longevity, which reflects the ability of these cattle to remain productive under severe environmental conditions. Genetic markers were also identified as significantly associated with age at first calving and calf growth rate. The results indicate that genomic information can be used to control inbreeding and to implement genomic selection in Ramo Grande cattle to enhance adaptation and production traits.

Abstract

Ramo Grande is a local cattle breed raised in the archipelago of Azores, with a small and dispersed census, where inbreeding control is of utmost importance. A single nucleotide polymorphism (SNP) Beadchip array was used to assess inbreeding, by analysis of genomic regions harboring contiguous homozygous genotypes named runs of homozygosity (ROH), and to estimate past effective population size by analysis of linkage disequilibrium (LD). Genetic markers associated with production traits were also investigated, exploiting the unique genetic and adaptation features of this breed. A total of 639 ROH with length >4 Mb were identified, with mean length of 14.96 Mb. The mean genomic inbreeding was 0.09, and long segments of ROH were common, indicating recent inbred matings. The LD pattern indicates a large effective population size, suggesting the inflow of exotic germplasm in the past. The genome-wide association study identified novel markers significantly affecting longevity, age at first calving and direct genetic effects on calf weight. These results provide the first evidence of the association of longevity with genes related with DNA recognition and repair, and the association of age at first calving with aquaporin proteins, which are known to have a crucial role in reproduction.

Allele-specific expression is widespread in Bos indicus muscle and affects meat quality candidate genes

Differences between the expression of the two alleles of a gene are known as allele-specific expression (ASE), a common event in the transcriptome of mammals. Despite ASE being a source of phenotypic variation, its occurrence and effects on genetic prediction of economically relevant traits are still unexplored in bovines. Furthermore, as ASE events are likely driven by cis-regulatory mutations, scanning them throughout the bovine genome represents a significant step to elucidate the mechanisms underlying gene expression regulation. To address this question in a Bos indicus population, we built the ASE profile of the skeletal muscle tissue of 190 Nelore steers, using RNA sequencing data and SNPs genotypes from the Illumina BovineHD BeadChip (770 K bp). After quality control, 820 SNPs showed at least one sample with ASE. These SNPs were widespread among all autosomal chromosomes, being 32.01% found in 3'UTR and 31.41% in coding regions. We observed a considerable variation of ASE profile among individuals, which highlighted the need for biological replicates in ASE studies. Functional analysis revealed that ASE genes play critical biological functions in the development and maintenance of muscle tissue. Additionally, some of these genes were previously reported as associated with beef production and quality traits in livestock, thus indicating a possible source of bias on genomic predictions for these traits.

Landscape of multi-tissue global gene expression reveals the regulatory signatures of feed efficiency in beef cattle

MOTIVATION

Feed efficiency is an important trait for sustainable beef production that is regulated by the complex biological process, but the mode of action behinds it has not been clearly defined. Here, we aimed to elucidate the regulatory mechanisms of this trait through studying the landscape of the genome-wide gene expression of rumen, liver, muscle and backfat tissues, the key ones involved in the energy metabolism.

RESULTS

The transcriptome of 189 samples across four tissues from 48 beef steers with varied feed efficiency were generated using Illumina HiSeq4000. The analysis of global gene expression profiles of four tissues, functional analysis of tissue-shared and -unique genes, co-expressed network construction of tissue-shared genes, weighted correlations analysis between gene modules and feed efficiency-related traits in each tissue were performed. Among four tissues, the transcriptome of muscle tissue was distinctive from others, while those of rumen and backfat tissues were similar. The associations between co-expressed genes and feed efficiency related traits at single or all tissues level exhibited that the gene expression in the rumen, liver, muscle and backfat were the most correlated with feed conversion ratio, dry matter intake, average daily gain and residual feed intake, respectively. The 19 overlapped genes identified from the strongest module-trait relationships in four tissues are potential generic gene markers for feed efficiency.

AVAILABILITY AND IMPLEMENTATION

The distribution of gene expression data can be accessed at https://www.cattleomics.com/transcriptome.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identification and Annotation of Potential Function of Regulatory Antisense Long Non-Coding RNAs Related to Feed Efficiency in Bos taurus Bulls

Long non-coding RNAs (lncRNAs) can influence transcriptional and translational processes in mammalian cells and are associated with various developmental, physiological and phenotypic conditions. However, they remain poorly understood and annotated in livestock species. We combined phenotypic, metabolomics and liver transcriptomic data of bulls divergent for residual feed intake (RFI) and fat accretion. Based on a project-specific transcriptome annotation for the bovine reference genome ARS-UCD.1.2 and multiple-tissue total RNA sequencing data, we predicted 3590 loci to be lncRNAs. To identify lncRNAs with potential regulatory influence on phenotype and gene expression, we applied the regulatory impact factor algorithm on a functionally prioritized set of loci (n = 4666). Applying the algorithm of partial correlation and information theory, significant and independent pairwise correlations were calculated and co-expression networks were established, including plasma metabolites correlated with lncRNAs. The network hub lncRNAs were assessed for potential cis-actions and subjected to biological pathway enrichment analyses. Our results reveal a prevalence of antisense lncRNAs positively correlated with adjacent protein-coding genes and suggest their participation in mitochondrial function, acute phase response signalling, TCA-cycle, fatty acid β-oxidation and presumably gluconeogenesis. These antisense lncRNAs indicate a stabilizing function for their cis-correlated genes and a putative regulatory role in gene expression.

Genome-wide detection of signatures of selection in three Valdostana cattle populations

The Valdostana is a local dual purpose cattle breed developed in Italy. Three populations are recognized within this breed, based on coat colour, production level, morphology and temperament: Valdostana Red Pied (VPR), Valdostana Black Pied (VPN) and Valdostana Chestnut (VCA). Here, we investigated putative genomic regions under selection among these three populations using the Bovine 50K SNP array by combining three different statistical methods based either on allele frequencies (F_ST ) or extended haplotype homozygosity (iHS and Rsb). In total, 8, 5 and 8 chromosomes harbouring 13, 13 and 16 genomic regions potentially under selection were identified by at least two approaches in VPR, VPN and VCA, respectively. Most of these candidate regions were population-specific but we found one common genomic region spanning 2.38 Mb on BTA06 which either overlaps or is located close to runs of homozygosity islands detected in the three populations. This region included inter alia two well-known genes: KDR, a well-established coat colour gene, and CLOCK, which plays a central role in positive regulation of inflammatory response and in the regulation of the mammalian circadian rhythm. The other candidate regions identified harboured genes associated mainly with milk and meat traits as well as genes involved in immune response/inflammation or associated with behavioural traits. This last category of genes was mainly identified in VCA, which is selected for fighting ability. Overall, our results provide, for the first time, a glimpse into regions of the genome targeted by selection in Valdostana cattle. Finally, this study illustrates the relevance of using multiple complementary approaches to identify genomic regions putatively under selection in livestock.

Potential Biomarkers for Feed Efficiency-Related Traits in Nelore Cattle Identified by Co-expression Network and Integrative Genomics Analyses

Feed efficiency helps to reduce environmental impacts from livestock production, improving beef cattle profitability. We identified potential biomarkers (hub genes) for feed efficiency, by applying co-expression analysis in Longissimus thoracis RNA-Seq data from 180 Nelore steers. Six co-expression modules were associated with six feed efficiency-related traits (p-value ≤ 0.05). Within these modules, 391 hub genes were enriched for pathways as protein synthesis, muscle growth, and immune response. Trait-associated transcription factors (TFs) ELF1, ELK3, ETS1, FLI1, and TCF4, were identified with binding sites in at least one hub gene. Gene expression of CCDC80, FBLN5, SERPINF1, and OGN was associated with multiple feed efficiency-related traits (FDR ≤ 0.05) and were previously related to glucose homeostasis, oxidative stress, fat mass, and osteoblastogenesis, respectively. Potential regulatory elements were identified, integrating the hub genes with previous studies from our research group, such as the putative cis-regulatory elements (eQTLs) inferred as affecting the PCDH18 and SPARCL1 hub genes related to immune system and adipogenesis, respectively. Therefore, our analyses contribute to a better understanding of the biological mechanisms underlying feed efficiency in bovine and the hub genes disclosed can be used as biomarkers for feed efficiency-related traits in Nelore cattle.

Integrated metabolomic and transcriptomic analyses identify critical genes in eicosapentaenoic acid biosynthesis and metabolism in the sea urchin Strongylocentrotus intermedius

Gonads are the only edible part of the sea urchin and have great potential as a health-promoting food for human consumption. Polyunsaturated fatty acids (PUFAs) are important necessary nutrients that determine not only the nutritional value of sea urchins but guarantee their normal growth and reproduction. However, the information on the molecular mechanisms of PUFA biosynthesis and metabolism in this species remains elusive. In this study, we used Strongylocentrotus intermedius as our model species and conducted integrated metabolomic and transcriptomic analyses of potentially critical genes involved in PUFA biosynthesis and metabolism during gonad growth and development, mainly focusing on eicosapentaenoic acid (EPA). We found six differentially accumulated metabolites associated with PUFA in the metabolomic analysis. More differentially expressed genes (DEGs) were related to PUFA in testis than ovary (1823 DEGs in testis and 1499 DEGs in ovary). We verified 12 DEGs by RNA-Seq results and found that Aldh7a1, Ecm3, Fads2, and Hsd17b12 genes had similar expression patterns in EPA concentration during gonad growth and development. In contrast, the other DEGs were downregulated and we inferred that EPA or PUFA may be metabolized as energy during certain periods. Our metabolic and genetic data will facilitate a better understanding of PUFA regulation networks during gonad growth and development in S. intermedius.

Identification of genetic loci associated with growth traits at weaning in yak through a genome-wide association study

A multilocus GWAS was performed to explore the genetic architecture of four growth traits in yak. In total, 354 female yaks for which measurements of body weight (BW), withers height (WH), body length (BL) and chest girth (CG) at weaning were available underwent genotyping with the Illumina BovineHD BeadChip (770K). After quality control, we retained 98 688 SNPs and 354 animals for GWAS analysis. We identified seven, 18, seven and nine SNPs (corresponding to seven, 17, seven and eight candidate genes) associated with BW, WH, BL and CG at weaning respectively. Interestingly, most of these candidate genes were reported to be involved in growth-related processes such as muscle formation, lipid deposition, feed efficiency, carcass composition and development of the central and peripheral nervous system. Our results offer novel insight into the molecular architecture underpinning yak growth traits. Further functional analyses are thus warranted to explore the molecular mechanisms whereby these genes affect these traits of interest.

Combining information from genome-wide association and multi-tissue gene expression studies to elucidate factors underlying genetic variation for residual feed intake in Australian Angus cattle

BACKGROUND

Genome-wide association studies (GWAS) are extensively used to identify single nucleotide polymorphisms (SNP) underlying the genetic variation of complex traits. However, much uncertainly often still exists about the causal variants and genes at quantitative trait loci (QTL). The aim of this study was to identify QTL associated with residual feed intake (RFI) and genes in these regions whose expression is also associated with this trait. Angus cattle (2190 steers) with RFI records were genotyped and imputed to high density arrays (770 K) and used for a GWAS approach to identify QTL associated with RFI. RNA sequences from 126 Angus divergently selected for RFI were analyzed to identify the genes whose expression was significantly associated this trait with special attention to those genes residing in the QTL regions.

RESULTS

The heritability for RFI estimated for this Angus population was 0.3. In a GWAS, we identified 78 SNPs associated with RFI on six QTL (on BTA1, BTA6, BTA14, BTA17, BTA20 and BTA26). The most significant SNP was found on chromosome BTA20 (rs42662073) and explained 4% of the genetic variance. The minor allele frequencies of significant SNPs ranged from 0.05 to 0.49. All regions, except on BTA17, showed a significant dominance effect. In 1 Mb windows surrounding the six significant QTL, we found 149 genes from which OAS2, STC2, SHOX, XKR4, and SGMS1 were the closest to the most significant QTL on BTA17, BTA20, BTA1, BTA14, and BTA26, respectively. In a 2 Mb windows around the six significant QTL, we identified 15 genes whose expression was significantly associated with RFI: BTA20) NEURL1B and CPEB4; BTA17) RITA1, CCDC42B, OAS2, RPL6, and ERP29; BTA26) A1CF, SGMS1, PAPSS2, and PTEN; BTA1) MFSD1 and RARRES1; BTA14) ATP6V1H and MRPL15.

CONCLUSIONS

Our results showed six QTL regions associated with RFI in a beef Angus population where five of these QTL contained genes that have expression associated with this trait. Therefore, here we show that integrating information from gene expression and GWAS studies can help to better understand the genetic mechanisms that determine variation in complex traits.

Liver transcriptome profiling of beef steers with divergent growth rate, feed intake, or metabolic body weight phenotypes1

Average daily gain (ADG) and daily dry matter intake (DMI) are key determinants of beef industry profitability. These traits together with metabolic body weight (MWT) are combined as component traits to calculate residual feed intake (RFI), a common measure of feed efficiency in beef cattle. Recently, there have been significant efforts towards molecular genetic characterization of RFI through transcriptomic studies in different breeds and tissues. However, molecular mechanisms of RFI component traits still remain predominately unexplored. Therefore, in the current study, we investigated the hepatic transcriptomic profiles and their associations with ADG, DMI, and MWT in Angus, Charolais, and Kinsella Composite (KC) populations through global RNAseq analyses. In each population and for each trait, 12 steers with extreme phenotypes (n = 6 low and n = 6 high) were analyzed for differential gene expression. These animals were from 20 beef steers of each Angus, Charolais, and KC breed population that were initially selected for a transcriptome study of RFI. At a false discovery rate <0.05 and fold change >1.5, we identified 123, 102, and 78 differentially expressed (DE) genes between high- and low-ADG animals of Angus, Charolais, and KC populations, respectively. For DMI, 108, 180, and 156 DE genes were identified between high- and low-DMI from Angus, Charolais, and KC populations, respectively, while for MWT, 80, 82, and 84 genes were differentially expressed between high- and low-MWT animals in Angus, Charolais, and KC populations, respectively. The identified DE genes were largely breed specific (81.7% for ADG, 82.7% for DMI, and 83% for MWT), but were largely involved in the same biological functions across the breeds. Among the most enriched biological functions included metabolism of major nutrients (lipids, carbohydrates, amino acids, vitamins, and minerals), small molecule biochemistry, cellular movement, cell morphology, and cell-to-cell signaling and interaction. Notably, we identified multiple DE genes that are involved in cholesterol biosynthesis, and immune response pathways for the 3 studied traits. Thus, our findings present potential molecular genetic mechanisms and candidate genes that influence feed intake, growth, and MWT of beef cattle.