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

Scans for Signatures of Selection in Genomes of Wagyu and Buryat Cattle Breeds Reveal Candidate Genes and Genetic Variants for Adaptive Phenotypes and Production Traits

Past and ongoing selection shapes the genomes of livestock breeds. Identifying such signatures of selection allows for uncovering the genetic bases of affected phenotypes, including economically important traits and environmental adaptations, for the further improvement of breed genetics to respond to climate and economic challenges. Turano-Mongolian cattle are a group of taurine breeds known for their adaptation to extreme environmental conditions and outstanding production performance. Buryat Turano-Mongolian cattle are among the few breeds adapted to cold climates and poor forage. Wagyu, on the other hand, is famous for high productivity and unique top-quality marbled meat. We used hapFLK, the de-correlated composite of multiple signals (DCMS), PBS, and FST methods to search for signatures of selection in their genomes. The scans revealed signals in genes related to cold adaptation (e.g., STAT3, DOCK5, GSTM3, and CXCL8) and food digestibility (SI) in the Buryat breed, and growth and development traits (e.g., RBFOX2 and SHOX2) and marbling (e.g., DGAT1, IQGAP2, RSRC1, and DIP2B) in Wagyu. Several putatively selected genes associated with reproduction, immunity, and resistance to pathogens were found in both breed genomes. The results of our work could be used for creating new productive adapted breeds or improving the extant breeds.

Genes involved in the cholecystokinin receptor signaling map were differentially expressed in the jejunum of steers with variation in residual feed intake

The jejunum is a critical site for nutrient digestion and absorption, and variation in its ability to take up nutrients within the jejunum is likely to affect feed efficiency. The purpose of this study was to determine differences in gene expression in the jejunum of beef steers divergent for residual feed intake (RFI) in one cohort of steers (Year 1), and to validate those genes in animals from a second study (Year 2). Steers from Year 1 (n = 16) were selected for high and low RFI. Jejunum mucosal tissue was obtained for RNA-seq. Thirty-two genes were differentially expressed (PFDR≤0.15), and five were over-represented in pathways including inflammatory mediator, cholecystokinin receptor (CCKR) signaling, and p38 MAPK pathways. Several differentially expressed genes (ALOX12, ALPI, FABP6, FABP7, FLT1, GSTA2, MEF2B, PDK4, SPP1, and TTF2) have been previously associated with RFI in other studies. Real-time qPCR was used to validate nine differentially expressed genes in the Year 1 steers used for RNA-seq, and in the Year 2 validation cohort. Six genes were validated as differentially expressed (P < 0.1) using RT-qPCR in the Year 1 population. In the Year 2 population, five genes displayed the same direction of expression as the Year 1 population and 3 were differentially expressed (P < 0.1). The CCKR pathway is involved in digestion, appetite control, and regulation of body weight making it a compelling candidate for feed efficiency in cattle, and the validation of these genes in a second population of cattle is suggestive of a role in feed efficiency.

Analyses of plasma metabolites using a high performance four-channel CIL LC-MS method and identification of metabolites associated with enteric methane emissions in beef cattle

Reducing enteric methane (one greenhouse gas) emissions from beef cattle not only can be beneficial in reducing global warming, but also improve efficiency of nutrient utilization in the production system. However, direct measurement of enteric methane emissions on individual cattle is difficult and expensive. The objective of this study was to detect plasma metabolites that are associated with enteric methane emissions in beef cattle. Average enteric methane emissions (CH4) per day (AVG_DAILYCH4) for each individual cattle were measured using the GreenFeed emission monitoring (GEM) unit system, and beef cattle with divergent AVG_DAILYCH4 from Angus (n = 10 for the low CH4 group and 9 for the high CH4 group), Charolais (n = 10 for low and 10 for = high), and Kinsella Composite (n = 10 for low and 10 for high) populations were used for plasma metabolite quantification and metabolite-CH4 association analyses. Blood samples of these cattle were collected near the end of the GEM system tests and a high performance four-channel chemical isotope labeling (CIL) liquid chromatography (LC) mass spectrometer (MS) method was applied to identify and quantify concentrations of metabolites. The four-channel CIL LC-MS method detected 4235 metabolites, of which 1105 were found to be significantly associated with AVG_DAILYCH4 by a t-test, while 1305 were significantly associated with AVG_DAILYCH4 by a regression analysis at p<0.05. Both the results of the t-test and regression analysis revealed that metabolites that were associated with enteric methane emissions in beef cattle were largely breed-specific whereas 4.29% to 6.39% CH4 associated metabolites were common across the three breed populations and 11.07% to 19.08% were common between two breed populations. Pathway analyses of the CH4 associated metabolites identified top enriched molecular processes for each breed population, including arginine and proline metabolism, arginine biosynthesis, butanoate metabolism, and glutathione metabolism for Angus; beta-alanine metabolism, pyruvate metabolism, glycolysis / gluconeogenesis, and citrate cycle (TCA cycle) for Charolais; phenylalanine, tyrosine and tryptophan biosynthesis, phenylalanine metabolism, arginine biosynthesis, and arginine and proline metabolism for Kinsella Composite. The detected CH4 associated metabolites and enriched molecular processes will help understand biological mechanisms of enteric methane emissions in beef cattle. The detected CH4 associated plasma metabolites will also provide valuable resources to further characterize the metabolites and verify their utility as biomarkers for selection of cattle with reduced methane emissions.

Identification of Candidate Genes and Functional Pathways Associated with Body Size Traits in Chinese Holstein Cattle Based on GWAS Analysis

Body size is one of the most economically important traits of dairy cattle, as it is significantly associated with cow longevity, production, health, fertility, and environmental adaptation. The identification and application of genetic variants using a novel genetic approach, such as genome-wide association studies (GWASs), may give more insights into the genetic architecture of complex traits. The identification of genes, single nucleotide polymorphisms (SNPs), and pathways associated with the body size traits may offer a contribution to genomic selection and long-term planning for selection in dairy cows. In this study, we performed GWAS analysis to identify the genetic markers and genes associated with four body size traits (body height, body depth, chest width, and angularity) in 1000 Chinese Holstein cows. We performed SNPs genotyping in 1000 individuals, based on the GeneSeek Genomic Profiler Bovine 100 K. In total, we identified 11 significant SNPs in association with body size traits at the threshold of Bonferroni correction (5.90 × 10−7) using the fixed and random model circulating probability unification (FarmCPU) model. Several genes within 200 kb distances (upstream or downstream) of the significant SNPs were identified as candidate genes, including MYH15, KHDRBS3, AIP, DCC, SQOR, and UBAP1L. Moreover, genes within 200 kb of the identified SNPs were significantly enriched (p ≤ 0.05) in 25 Gene Ontology terms and five Kyoto Encyclopedia of Genes and Genomes pathways. We anticipate that these results provide a foundation for understanding the genetic architecture of body size traits. They will also contribute to breeding programs and genomic selection work on Chinese Holstein cattle.

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

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

High-Throughput Sequencing Reveals Transcriptome Signature of Early Liver Development in Goat Kids

As a vital metabolic and immune organ in animals, the liver plays an important role in protein synthesis, detoxification, metabolism, and immune defense. The primary research purpose of this study was to reveal the effect of breast-feeding, weaning transition, and weaning on the gene expression profile in the goat kid liver and to elucidate the transcriptome-level signatures associated with liver metabolic adaptation. Therefore, transcriptome sequencing was performed on liver tissues, which was collected at 1 day (D1), 2 weeks (W2), 4 weeks (W4), 8 weeks (W8), and 12 weeks (W12) after birth in Laiwu black goats at five different time-points, with five goats at each time point. From 25 libraries, a total of 37497 mRNAs were found to be expressed in goat kid livers, and 1271 genes were differentially expressed between at least two of the five time points. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses revealed that these genes were annotated as an extracellular region fraction, exhibiting monooxygenase activity, positive regulation of T cell activation, mitotic spindle mid-region assembly, cytokinesis, cytoskeleton-dependent cytokinesis, regulation of cytokinesis, regulation of lymphocyte proliferation, and so on. In addition, it mainly deals with metabolism, endocrine, cell proliferation and apoptosis, and immune processes. Finally, a gene regulatory network was constructed, and a total of 14 key genes were screened, which were mainly enriched for cell growth and development, endocrine, immune, and signal transduction-related pathways. Our results provide new information on the molecular mechanisms and pathways involved in liver development, metabolism, and immunity of goats.

Cross-population selection signatures in Canchim composite beef cattle

Analyses of livestock genomes have been used to detect selection signatures, which are genomic regions associated with traits under selection leading to a change in allele frequency. The objective of the present study was to characterize selection signatures in Canchim composite beef cattle using cross-population analyses with the founder Nelore and Charolais breeds. High-density single nucleotide polymorphism genotypes were available on 395 Canchim representing the target population, along with genotypes from 809 Nelore and 897 Charolais animals representing the reference populations. Most of the selection signatures were co-located with genes whose functions agree with the expectations of the breeding programs; these genes have previously been reported to associate with meat quality, as well as reproductive traits. Identified genes were related to immunity, adaptation, morphology, as well as behavior, could give new perspectives for understanding the genetic architecture of Canchim. Some selection signatures identified genes that were recently introduced in Canchim, such as the loci related to the polled trait.

Genome-wide association study for meat tenderness in beef cattle identifies patterns of the genetic contribution in different post-mortem stages

The beef tenderization process during the post-mortem period is one of the most important sensorial attributes and it is well-established. The aim of this study was to identify the genetic contribution pattern to meat tenderness at 7-(LMD7), 14-(LMD14), and 21-(LMD21) days post-mortem. The heritabilities for LMD7 (0.194), LMD14 (0.142) and LMD21 (0.048) are well established in the population evaluated here. However, its genetic contribution in terms of genomic candidate regions is still poorly understood. Tenderness was measured in the Longissiums thoracis using Warner-Bratzler shear force in the three post-mortem periods. A total of 4323 crossbred beef cattle were phenotyped and genotyped using the Illumina BovineSNP50K. The percentage of the total genetic variance was estimated using the weighted single-step genomic best linear unbiased prediction method. The main candidate windows for LMD7 were associated with proteolysis of myofibrillar structures and the weakening endomysium and perimysium. Candidate windows for LMD14 and LMD21 were mapped in bovine QTLs for body composition, height and growth. Results presented herein highlight, the largest contribution of proteolysis related processes before 14-days post-mortem and body composition characteristics in later stages for meat tenderness.

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.

Integrated Environmental and Genomic Analysis Reveals the Drivers of Local Adaptation in African Indigenous Chickens

Breeding for climate resilience is currently an important goal for sustainable livestock production. Local adaptations exhibited by indigenous livestock allow investigating the genetic control of this resilience. Ecological niche modeling (ENM) provides a powerful avenue to identify the main environmental drivers of selection. Here, we applied an integrative approach combining ENM with genome-wide selection signature analyses (XPEHH and Fst) and genotype-environment association (redundancy analysis), with the aim of identifying the genomic signatures of adaptation in African village chickens. By dissecting 34 agro-climatic variables from the ecosystems of 25 Ethiopian village chicken populations, ENM identified six key drivers of environmental challenges: One temperature variable-strongly correlated with elevation, three precipitation variables as proxies for water availability, and two soil/land cover variables as proxies of food availability for foraging chickens. Genome analyses based on whole-genome sequencing (n = 245), identified a few strongly supported genomic regions under selection for environmental challenges related to altitude, temperature, water scarcity, and food availability. These regions harbor several gene clusters including regulatory genes, suggesting a predominantly oligogenic control of environmental adaptation. Few candidate genes detected in relation to heat-stress, indicates likely epigenetic regulation of thermo-tolerance for a domestic species originating from a tropical Asian wild ancestor. These results provide possible explanations for the rapid past adaptation of chickens to diverse African agro-ecologies, while also representing new landmarks for sustainable breeding improvement for climate resilience. We show that the pre-identification of key environmental drivers, followed by genomic investigation, provides a powerful new approach for elucidating adaptation in domestic animals.

An integrated approach in gene-expression landscape profiling to identify housekeeping and tissue-specific genes in cattle

Context High-throughput transcriptome sequencing (RNA-Seq) has been widely applied in cattle studies. Public databases such as the National Center for Biotechnology Information (NCBI) contain large collections of gene expression data from various cattle tissues that can be used in gene expression analysis research Aims This study was conducted to investigate patterns of transcriptome variation across tissues of cattle through large-scale identification of housekeeping genes (i.e. those crucial to maintaining basic cellular activity) and tissue-specific genes in cattle tissues. Methods Using data available in the NCBI Sequence Read Archive database, we analysed 1377 transcriptome data sequences from 60 bovine tissue types, identified tissue-specific and housekeeping genes, and set up a web-based bovine gene expression analysis tool. Key results We found 101 genes widely expressed in almost all tissue and screened out five housekeeping genes: RPL35A, eIF4A2, GAPDH, IPO5 and PAK2. Focusing on 12 major organs, we found 861 genes specifically expressing in these tissues. Furthermore, 187 significantly differentially expressed genes were found among six types of muscle tissues. All expression data were made available at our new website http://cattleExp.org, which can be freely accessed for future gene expression analyses. Conclusions The housekeeping genes and tissue-specific genes identified will provide more information for researchers studying gene expression in cattle. Implications The web-based cattle gene expression analysis tool will make it easy for researchers to access large public datasets. Users can easily access all publicly available RNA data and upload their own RNA-Seq data.

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.

Weighted single-step genome-wide association study and pathway analyses for feed efficiency traits in Nellore cattle

The aim was to conduct a weighted single-step genome-wide association study to detect genomic regions and putative candidate genes related to residual feed intake, dry matter intake, feed efficiency (FE), feed conversion ratio, residual body weight gain, residual intake and weight gain in Nellore cattle. Several protein-coding genes were identified within the genomic regions that explain more than 0.5% of the additive genetic variance for these traits. These genes were associated with insulin, leptin, glucose, protein and lipid metabolisms; energy balance; heat and oxidative stress; bile secretion; satiety; feed behaviour; salivation; digestion; and nutrient absorption. Enrichment analysis revealed functional pathways (p-value < .05) such as neuropeptide signalling (GO:0007218), negative regulation of canonical Wingless/Int-1 (Wnt) signalling (GO:0090090), bitter taste receptor activity (GO:0033038), neuropeptide hormone activity (GO:0005184), bile secretion (bta04976), taste transduction (bta0742) and glucagon signalling pathway (bta04922). The identification of these genes, pathways and their respective functions should contribute to a better understanding of the genetic and physiological mechanisms regulating Nellore FE-related traits.

Integrative Analysis of Metabolomic and Transcriptomic Profiles Uncovers Biological Pathways of Feed Efficiency in Pigs

Feed efficiency (FE) is an economically important trait. Thus, reliable predictors would help to reduce the production cost and provide sustainability to the pig industry. We carried out metabolome-transcriptome integration analysis on 40 purebred Duroc and Landrace uncastrated male pigs to identify potential gene-metabolite interactions and explore the molecular mechanisms underlying FE. To this end, we applied untargeted metabolomics and RNA-seq approaches to the same animals. After data quality control, we used a linear model approach to integrate the data and find significant differently correlated gene-metabolite pairs separately for the breeds (Duroc and Landrace) and FE groups (low and high FE) followed by a pathway over-representation analysis. We identified 21 and 12 significant gene-metabolite pairs for each group. The valine-leucine-isoleucine biosynthesis/degradation and arginine-proline metabolism pathways were associated with unique metabolites. The unique genes obtained from significant metabolite-gene pairs were associated with sphingolipid catabolism, multicellular organismal process, cGMP, and purine metabolic processes. While some of the genes and metabolites identified were known for their association with FE, others are novel and provide new avenues for further research. Further validation of genes, metabolites, and gene-metabolite interactions in larger cohorts will elucidate the regulatory mechanisms and pathways underlying FE.

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.