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

BuffExDb: web-based tissue-specific gene expression resource for breeding and conservation programmes in Bubalus bubalis

Amidst the global challenge of extreme poverty, the livestock sector can significantly contribute to global sustainable development goals by enhancing resilience, smallholder productivity, and market participation. The Indian livestock sector is one of the largest in the world with a total livestock population of 535.82 million, ∼10.7% of the world's livestock population. Buffalo (Bubalus bubalis) holds significant importance in India and other Asian countries, notably contributing to their economies by surpassing cattle in milk production and providing various valuable products. The limited availability of genomic and transcriptomic resources for buffaloes hinders the efforts to enhance their traits for increased milk and meat production. To address this gap, this study adopted the state-of-the-art bioinformatics tools to analyse 2429 transcriptomes representing 438 BioSamples from 23 BioProjects obtained from a public domain database, representing 76 different types of tissues and cell types from all major organ systems in buffalo species (river and swamp). The outcome of this exhaustive genomic data led to the development of a relational buffalo expression database based on a three-tier architecture named as BuffExDb (http://46.202.167.198/buffex/). The user-friendliness and flexibilities in retrieval of tissue-specific genes (TSGs) and their functional annotation are the major characteristics of BuffExDb. This is the first of its kind that offers an effortlessly navigable and filterable database, enabling users to examine and visualize the expression levels of each tissue across multiple samples, simultaneously. It also provides the Tau score parameter for the identification of TSGs along with their essential roles in tissue development, maintenance, and function as observed through the enrichment test for gene ontologies. The exhaustive outcome of this work would pave the way for the biological, functional, and evolutionary studies for easy access. This prior information based on tissue-specific mechanisms can be used for future genomic research, especially in association studies in endeavour of enhanced buffalo breeding and conservation programmes. Database URL: http://46.202.167.198/buffex/.

Unraveling Ruminant Feed Efficiency Through Metabolomics: A Systematic Review

BACKGROUND

Advancements in metabolomic technologies have revolutionized our understanding of feed efficiency (FE) in livestock, offering new pathways to enhance both profitability and sustainability in ruminant production.

METHODS

This review offers a critical and systematic evaluation of the metabolomics methods used to measure and assess FE in ruminants. We conducted a comprehensive search of PubMed, Web of Science, and Scopus databases, covering publications from 1971 to 2023. This review synthesizes findings from 71 studies that applied metabolomic approaches to uncover the biological mechanisms driving interindividual variations in FE across cattle, sheep, goats, and buffaloes.

RESULTS

Most studies focused on cattle and employed targeted metabolomics to identify key biomarkers, including amino acids, fatty acids, and other metabolites linked to critical pathways such as energy metabolism, nitrogen utilization, and muscle development. Despite promising insights, challenges remain, including small sample sizes, methodological inconsistencies, and a lack of validation studies, particularly for non-cattle species.

CONCLUSIONS

By leveraging state-of-the-art metabolomic methods, this review highlights the potential of metabolomics to provide cost-effective, non-invasive molecular markers for FE evaluation, paving the way for more efficient and sustainable livestock management. Future research should prioritize larger, species-specific studies with standardized methods to validate identified biomarkers and enhance practical applications in livestock production systems.

Application of GWAS and mGWAS in Livestock and Poultry Breeding

In recent years, genome-wide association studies (GWAS) and metabolome genome-wide association studies (mGWAS) have emerged as crucial methods for investigating complex traits in animals and plants. These have played pivotal roles in research on livestock and poultry breeding, facilitating a deeper understanding of genetic diversity, the relationship between genes, and genetic bases in livestock and poultry. This article provides a review of the applications of GWAS and mGWAS in animal genetic breeding, aiming to offer reference and inspiration for relevant researchers, promote innovation in animal genetic improvement and breeding methods, and contribute to the sustainable development of animal husbandry.

Some plasma biomarkers of residual feed intake in beef cattle remain consistent regardless of intake level

This study investigated whether plasma biomarkers of residual feed intake (RFI), identified under ad libitum feeding conditions in beef cattle, remained consistent during feed restriction. Sixty Charolais crossbred young bulls were divided into two groups for a crossover study. Group A was initially fed ad libitum (first test) and then restricted (second test) on the same diet, while Group B experienced the opposite sequence. Blood samples were collected from the 12 most divergent RFI animals in each group at the end of the first test and again after the second test. 12 plasma variables consistently increased, while three consistently decreased during feed restriction (FDR < 0.05). Only two metabolites, α-aminoadipic acid for Group A and 5-aminovaleric acid for Group B, were associated with RFI independent of feed intake level (FDR < 0.05), demonstrating moderate-to-high repeatability across feeding levels (intraclass correlation coefficient ≥ 0.59). Notably, both metabolites belong to the same metabolic pathway: lysine degradation. These metabolites consistently correlated with RFI, irrespective of fluctuations in feed intake, indicating a connection to individual metabolic processes influencing feed efficiency. These findings suggest that a portion of RFI phenotypic variance is inherent to an individual's metabolic efficiency beyond variations in feed intake.

Effects of a liquid and dry Saccharomyces cerevisiae fermentation product feeding program on ruminal fermentation, total tract digestibility, and plasma metabolome of Holstein steers receiving a grain-based diet

The study aimed to determine the effects of a postbiotic feeding program consisting of liquid and dry Saccharomyces cerevisiae fermentation product (SCFP) on ruminal fermentation, digestibility, and plasma metabolome of Holstein steers receiving a grain-based diet. Eight Holstein steers (body weight, BW, 467 ± 13.9 kg) equipped with rumen cannulas were used in a crossover design study, with 21 d per period and a 7-d washout period in between periods. Steers were stratified by initial BW and assigned to 1 of 2 treatments. The treatments were 1) Control, basal finishing diet only (CON); 2) SCFP, 1-d feeding of liquid SCFP (infused into the rumen via the cannula at 11 mL/100 kg BW) followed by daily feeding of dry SCFP (12 g/d, top-dressed). Feed and spot fecal samples were collected during days 17 to 20 for determination of digestibility and fecal excretion of N, P, Cu, and Zn. Digestibility was measured using acid-insoluble ash as an internal marker. Blood samples were collected on day 21 before the morning feeding. Rumen fluid samples were collected on days 0, 1, 2, 3, 5, and 21 via rumen cannula. Results were analyzed with the GLIMMIX procedure of SAS 9.4 (SAS, 2023). Treatment did not affect dry matter intake (P = 0.15) and digestibility (P ≥ 0.62). The fecal output and absorption of Zn, Cu, P, and N were not affected (P > 0.22) by treatment. On day 1, the liquid SCFP supplementation tended to reduce (P = 0.07) ruminal VFA concentration and increased (P < 0.01) the molar proportion of valerate. Feeding SCFP tended to increase total ruminal VFA on day 5 (P = 0.08) and significantly increased total VFA on day 21 (P = 0.05). Ruminal NH3-N was reduced (P = 0.02) on day 21 by supplementing SCFP. Treatment did not affect the production of proinflammatory cytokines, interleukin (IL)-1β (P > 0.19), and IL-6 (P > 0.12) in the whole blood in response to various toll-like receptor stimulants in vitro. Feeding SCFP enriched (P ≤ 0.05) plasma metabolic pathways, including citric acid cycle, pyrimidine metabolism, glycolysis/gluconeogenesis, retinol metabolism, and inositol phosphate metabolism pathways. In summary, supplementing liquid SCFP with subsequent dry SCFP enhanced ruminal total VFA production and reduced NH3-N concentration in the rumen. Furthermore, feeding SCFP enriched several important pathways in lipid, protein, and glucose metabolism, which may improve feed efficiency of energy and protein in Holstein steers.

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.

Untargeted metabolomics confirms the association between plasma branched chain amino acids and residual feed intake in beef heifers

This study explored plasma biomarkers and metabolic pathways underlying feed efficiency measured as residual feed intake (RFI) in Charolais heifers. A total of 48 RFI extreme individuals (High-RFI, n = 24; Low-RFI, n = 24) were selected from a population of 142 heifers for classical plasma metabolite and hormone quantification and plasma metabolomic profiling through untargeted LC-MS. Most efficient heifers (Low-RFI) had greater (P = 0.03) plasma concentrations of IGF-1 and tended to have (P = 0.06) a lower back fat depth compared to least efficient heifers. However, no changes were noted (P ≥ 0.10) for plasma concentrations of glucose, insulin, non-esterified fatty acids, β-hydroxybutyrate and urea. The plasma metabolomic dataset comprised 3,457 ions with none significantly differing between RFI classes after false discovery rate correction (FDR > 0.10). Among the 101 ions having a raw P < 0.05 for the RFI effect, 13 were putatively annotated by using internal databases and 6 compounds were further confirmed with standards. Metabolic pathway analysis from these 6 confirmed compounds revealed that the branched chain amino acid metabolism was significantly (FDR < 0.05) impacted by the RFI classes. Our results confirmed for the first time in beef heifers previous findings obtained in male beef cattle and pointing to changes in branched-chain amino acids metabolism along with that of body composition as biological mechanisms related to RFI. Further studies are warranted to ascertain whether there is a cause-and-effect relationship between these mechanisms and RFI.

Identification of biomarkers associated with the feed efficiency by metabolomics profiling: results from the broiler lines divergent for high or low abdominal fat content

Background

Improving feed efficiency (FE) is one of the main objectives in broiler breeding. It is difficult to directly measure FE traits, and breeders hence have been trying to identify biomarkers for the indirect selection and improvement of FE traits. Metabolome is the "bridge" between genome and phenome. The metabolites may potentially account for more of the phenotypic variation and can suitably serve as biomarkers for selecting FE traits. This study aimed to identify plasma metabolite markers for selecting high-FE broilers. A total of 441 birds from Northeast Agricultural University broiler lines divergently selected for abdominal fat content were used to analyze plasma metabolome and estimate the genetic parameters of differentially expressed metabolites.

Results

The results identified 124 differentially expressed plasma metabolites (P < 0.05) between the lean line (high-FE birds) and the fat line (low-FE birds). Among these differentially expressed plasma metabolites, 44 were found to have higher positive or negative genetic correlations with FE traits (|r_g| ≥ 0.30). Of these 44 metabolites, 14 were found to display moderate to high heritability estimates (h² ≥ 0.20). However, among the 14 metabolites, 4 metabolites whose physiological functions have not been reported were excluded. Ultimately, 10 metabolites were suggested to serve as the potential biomarkers for breeding the high-FE broilers. Based on the physiological functions of these metabolites, reducing inflammatory and improving immunity were proposed to improve FE and increase production efficiency.

Conclusions

According to the pipeline for the selection of the metabolite markers established in this study, it was suggested that 10 metabolites including 7-ketocholesterol, dimethyl sulfone, epsilon-(gamma-glutamyl)-lysine, gamma-glutamyltyrosine, 2-oxoadipic acid, L-homoarginine, testosterone, adenosine 5'-monophosphate, adrenic acid, and calcitriol could be used as the potential biomarkers for breeding the "food-saving broilers".

Genomic prediction with whole-genome sequence data in intensely selected pig lines

Background

Early simulations indicated that whole-genome sequence data (WGS) could improve the accuracy of genomic predictions within and across breeds. However, empirical results have been ambiguous so far. Large datasets that capture most of the genomic diversity in a population must be assembled so that allele substitution effects are estimated with high accuracy. The objectives of this study were to use a large pig dataset from seven intensely selected lines to assess the benefits of using WGS for genomic prediction compared to using commercial marker arrays and to identify scenarios in which WGS provides the largest advantage.

Methods

We sequenced 6931 individuals from seven commercial pig lines with different numerical sizes. Genotypes of 32.8 million variants were imputed for 396,100 individuals (17,224 to 104,661 per line). We used BayesR to perform genomic prediction for eight complex traits. Genomic predictions were performed using either data from a standard marker array or variants preselected from WGS based on association tests.

Results

The accuracies of genomic predictions based on preselected WGS variants were not robust across traits and lines and the improvements in prediction accuracy that we achieved so far with WGS compared to standard marker arrays were generally small. The most favourable results for WGS were obtained when the largest training sets were available and standard marker arrays were augmented with preselected variants with statistically significant associations to the trait. With this method and training sets of around 80k individuals, the accuracy of within-line genomic predictions was on average improved by 0.025. With multi-line training sets, improvements of 0.04 compared to marker arrays could be expected.

Conclusions

Our results showed that WGS has limited potential to improve the accuracy of genomic predictions compared to marker arrays in intensely selected pig lines. Thus, although we expect that larger improvements in accuracy from the use of WGS are possible with a combination of larger training sets and optimised pipelines for generating and analysing such datasets, the use of WGS in the current implementations of genomic prediction should be carefully evaluated against the cost of large-scale WGS data on a case-by-case basis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00756-0.

Integrative analyses of genomic and metabolomic data reveal genetic mechanisms associated with carcass merit traits in beef cattle

Improvement of carcass merit traits is a priority for the beef industry. Discovering DNA variants and genes associated with variation in these traits and understanding biological functions/processes underlying their associations are of paramount importance for more effective genetic improvement of carcass merit traits in beef cattle. This study integrates 10,488,742 imputed whole genome DNA variants, 31 plasma metabolites, and animal phenotypes to identify genes and biological functions/processes that are associated with carcass merit traits including hot carcass weight (HCW), rib eye area (REA), average backfat thickness (AFAT), lean meat yield (LMY), and carcass marbling score (CMAR) in a population of 493 crossbred beef cattle. Regression analyses were performed to identify plasma metabolites associated with the carcass merit traits, and the results showed that 4 (3-hydroxybutyric acid, acetic acid, citric acid, and choline), 6 (creatinine, L-glutamine, succinic acid, pyruvic acid, L-lactic acid, and 3-hydroxybutyric acid), 4 (fumaric acid, methanol, D-glucose, and glycerol), 2 (L-lactic acid and creatinine), and 5 (succinic acid, fumaric acid, lysine, glycine, and choline) plasma metabolites were significantly associated with HCW, REA, AFAT, LMY, and CMAR (P-value < 0.1), respectively. Combining the results of metabolome-genome wide association studies using the 10,488,742 imputed SNPs, 103, 160, 83, 43, and 109 candidate genes were identified as significantly associated with HCW, REA, AFAT, LMY, and CMAR (P-value < 1 × 10^-5), respectively. By applying functional enrichment analyses for candidate genes of each trait, 26, 24, 26, 24, and 28 significant cellular and molecular functions were predicted for HCW, REA, AFAT, LMY, and CMAR, respectively. Among the five topmost significantly enriched biological functions for carcass merit traits, molecular transport and small molecule biochemistry were two top biological functions associated with all carcass merit traits. Lipid metabolism was the most significant biological function for LMY and CMAR and it was also the second and fourth highest biological function for REA and HCW, respectively. Candidate genes and enriched biological functions identified by the integrative analyses of metabolites with phenotypic traits and DNA variants could help interpret the results of previous genome-wide association studies for carcass merit traits. Our integrative study also revealed additional potential novel genes associated with these economically important traits. Therefore, our study improves understanding of the molecular and biological functions/processes that influence carcass merit traits, which could help develop strategies to enhance genomic prediction of carcass merit traits with incorporation of metabolomic data. Similarly, this information could guide management practices, such as nutritional interventions, with the purpose of boosting specific carcass merit traits.