Genomic Heritability and Genome-Wide Association Studies of Plasma Metabolites in Crossbred Beef Cattle

Merging metabolomics and genomics provides a catalog of genetic factors that influence molecular phenotypes in pigs linking relevant metabolic pathways

BACKGROUND

Metabolomics opens novel avenues to study the basic biological mechanisms underlying complex traits, starting from characterization of metabolites. Metabolites and their levels in a biofluid represent simple molecular phenotypes (metabotypes) that are direct products of enzyme activities and relate to all metabolic pathways, including catabolism and anabolism of nutrients. In this study, we demonstrated the utility of merging metabolomics and genomics in pigs to uncover a large list of genetic factors that influence mammalian metabolism.

RESULTS

We obtained targeted characterization of the plasma metabolome of more than 1300 pigs from two populations of Large White and Duroc pig breeds. The metabolomic profiles of these pigs were used to identify genetically influenced metabolites by estimating the heritability of the level of 188 metabolites. Then, combining breed-specific genome-wide association studies of single metabolites and their ratios and across breed meta-analyses, we identified a total of 97 metabolite quantitative trait loci (mQTL), associated with 126 metabolites. Using these results, we constructed a human-pig comparative catalog of genetic factors influencing the metabolomic profile. Whole genome resequencing data identified several putative causative mutations for these mQTL. Additionally, based on a major mQTL for kynurenine level, we designed a nutrigenetic study feeding piglets that carried different genotypes at the candidate gene kynurenine 3-monooxygenase (KMO) varying levels of tryptophan and demonstrated the effect of this genetic factor on the kynurenine pathway. Furthermore, we used metabolomic profiles of Large White and Duroc pigs to reconstruct metabolic pathways using Gaussian Graphical Models, which included perturbation of the identified mQTL.

CONCLUSIONS

This study has provided the first catalog of genetic factors affecting molecular phenotypes that describe the pig blood metabolome, with links to important metabolic pathways, opening novel avenues to merge genetics and nutrition in this livestock species. The obtained results are relevant for basic and applied biology and to evaluate the pig as a biomedical model. Genetically influenced metabolites can be further exploited in nutrigenetic approaches in pigs. The described molecular phenotypes can be useful to dissect complex traits and design novel feeding, breeding and selection programs in pigs.

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.

Genome-wide association analysis identify candidate genes for feed efficiency and growth traits in Wenchang chickens

BACKGROUND

Wenchang chickens are one of the most popular local chicken breeds in the Chinese chicken industry. However, the low feed efficiency is the main shortcoming of this breed. Therefore, there is a need to find a more precise breeding method to improve the feed efficiency of Wenchang chickens. In this study, we explored important candidate genes and variants for feed efficiency and growth traits through genome-wide association study (GWAS) analysis.

RESULTS

Estimates of genomic heritability for growth and feed efficiency traits, including residual feed intake (RFI) of 0.05, average daily food intake (ADFI) of 0.21, average daily weight gain (ADG) of 0.24, body weight (BW) at 87, 95, 104, 113 days of age (BW87, BW95, BW104 and BW113) ranged from 0.30 to 0.44. Important candidate genes related to feed efficiency and growth traits were identified, such as PLCE1, LAP3, MED28, QDPR, LDB2 and SEL1L3 genes.

CONCLUSION

The results identified important candidate genes for feed efficiency and growth traits in Wenchang chickens and provide a theoretical basis for the development of new molecular breeding technology.

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.

Serum metabolic profile and metabolome genome-wide association study in chicken

BACKGROUND

Chickens provide globally important livestock products. Understanding the genetic and molecular mechanisms underpinning chicken economic traits is crucial for improving their selective breeding. Influenced by a combination of genetic and environmental factors, metabolites are the ultimate expression of physiological processes and can provide key insights into livestock economic traits. However, the serum metabolite profile and genetic architecture of the metabolome in chickens have not been well studied.

RESULTS

Here, comprehensive metabolome detection was performed using non-targeted LC-MS/MS on serum from a chicken advanced intercross line (AIL). In total, 7,191 metabolites were used to construct a chicken serum metabolomics dataset and to comprehensively characterize the serum metabolism of the chicken AIL population. Regulatory loci affecting metabolites were identified in a metabolome genome-wide association study (mGWAS). There were 10,061 significant SNPs associated with 253 metabolites that were widely distributed across the entire chicken genome. Many functional genes affect metabolite synthesis, metabolism, and regulation. We highlight the key roles of TDH and AASS in amino acids, and ABCB1 and CD36 in lipids.

CONCLUSIONS

We constructed a chicken serum metabolite dataset containing 7,191 metabolites to provide a reference for future chicken metabolome characterization work. Meanwhile, we used mGWAS to analyze the genetic basis of chicken metabolic traits and metabolites and to improve chicken breeding.

GWAS and genetic and phenotypic correlations of plasma metabolites with complete blood count traits in healthy young pigs reveal implications for pig immune response

Introduction: In this study estimated genetic and phenotypic correlations between fifteen complete blood count (CBC) traits and thirty-three heritable plasma metabolites in young healthy nursery pigs. In addition, it provided an opportunity to identify candidate genes associated with variation in metabolite concentration and their potential association with immune response, disease resilience, and production traits.

Methods: The blood samples were collected from healthy young pigs and Nuclear Magnetic Resonance (NMR) was used to quantify plasma metabolites. CBC was determined using the ADVIA® 2120i Hematology System. Genetic correlations of metabolite with CBC traits and single step genome-wide association study (ssGWAS) were estimated using the BLUPF90 programs.

Results: Results showed low phenotypic correlation estimates between plasma metabolites and CBC traits. The highest phenotypic correlation was observed between lactic acid and plasma basophil concentration (0.36 ± 0.04; p &lt; 0.05). Several significant genetic correlations were found between metabolites and CBC traits. The plasma concentration of proline was genetically positively correlated with hemoglobin concentration (0.94 ± 0.03; p &lt; 0.05) and L-tyrosine was negatively correlated with mean corpuscular hemoglobin (MCH; −0.92 ± 0.74; p &lt; 0.05). The genomic regions identified in this study only explained a small percentage of the genetic variance of metabolites levels that were genetically correlated with CBC, resilience, and production traits.

Discussion: The results of this systems approach suggest that several plasma metabolite phenotypes are phenotypically and genetically correlated with CBC traits, suggesting that they may be potential genetic indicators of immune response following disease challenge. Genomic analysis revealed genes and pathways that might interact to modulate CBC, resilience, and production traits.

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.

Genome-wide association study and genomic heritabilities for blood protein levels in Lori-Bakhtiari sheep

Serum protein levels are related to physiological and pathological status of animals and could be affected by both genetic and environmental factors. This study aimed to evaluate genetic variation of serum protein profile in sheep. Blood samples were randomly collected from 96 Lori-Bakhtiari ewes, a heavy meat-type breed. Total protein, albumin, globulin, α1, α2, β and γ globulins and IgG levels were measured in blood serum. The samples were genotyped using the Illumina OvineSNP50 BeadChip. The studied traits adjusted for age, birth type, birth season and estimate of breeding value for body weight were considered as pseudo-phenotypes in genome-wide association analysis. In the GWAS model, the first five principal components were fitted as covariates to correct the biases due to possible population stratification. The Plink, R and GCTA software were used for genome-wide association analysis, construction of Q-Q and Manhattan plots and estimation of genetic variances, respectively. Noticeable genomic heritabilities ± SE were estimated for total and γ globulins (0.868 ± 0.262 and 0.831 ± 0.364, respectively), but other protein fractions had zero or close to zero estimates. Based on the Bonferroni adjusted p values, four QTLs located on 181.7 Mbp of OAR3, 107.7 Mbp of OAR4, 86.3 Mbp of OAR7 and 83.0 Mbp of OAR8 were significantly associated with α1, β, β and γ globulins, respectively. The results showed that the PKP2, IGF2R, SLC22A1 and SLC22A2 genes could be considered as candidate genes for blood serum proteins. The present study showed significant genetic variations of some blood protein fractions.

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.

Overview of Metabolomic Analysis and the Integration with Multi-Omics for Economic Traits in Cattle

Metabolomics has been applied to measure the dynamic metabolic responses, to understand the systematic biological networks, to reveal the potential genetic architecture, etc., for human diseases and livestock traits. For example, the current published results include the detected relevant candidate metabolites, identified metabolic pathways, potential systematic networks, etc., for different cattle traits that can be applied for further metabolomic and integrated omics studies. Therefore, summarizing the applications of metabolomics for economic traits is required in cattle. We here provide a comprehensive review about metabolomic analysis and its integration with other omics in five aspects: (1) characterization of the metabolomic profile of cattle; (2) metabolomic applications in cattle; (3) integrated metabolomic analysis with other omics; (4) methods and tools in metabolomic analysis; and (5) further potentialities. The review aims to investigate the existing metabolomic studies by highlighting the results in cattle, integrated with other omics studies, to understand the metabolic mechanisms underlying the economic traits and to provide useful information for further research and practical breeding programs in cattle.

Heritability and genetic correlations of plasma metabolites of pigs with production, resilience and carcass traits under natural polymicrobial disease challenge

Metabolites in plasma of healthy nursery pigs were quantified using nuclear magnetic resonance. Heritabilities of metabolite concentration were estimated along with their phenotypic and genetic correlations with performance, resilience, and carcass traits in growing pigs exposed to a natural polymicrobial disease challenge. Variance components were estimated by GBLUP. Heritability estimates were low to moderate (0.11 ± 0.08 to 0.19 ± 0.08) for 14 metabolites, moderate to high (0.22 ± 0.09 to 0.39 ± 0.08) for 17 metabolites, and highest for l-glutamic acid (0.41 ± 0.09) and hypoxanthine (0.42 ± 0.08). Phenotypic correlation estimates of plasma metabolites with performance and carcass traits were generally very low. Significant genetic correlation estimates with performance and carcass traits were found for several measures of growth and feed intake. Interestingly the plasma concentration of oxoglutarate was genetically negatively correlated with treatments received across the challenge nursery and finisher (− 0.49 ± 0.28; P < 0.05) and creatinine was positively correlated with mortality in the challenge nursery (0.85 ± 0.76; P < 0.05). These results suggest that some plasma metabolite phenotypes collected from healthy nursery pigs are moderately heritable and genetic correlations with measures of performance and resilience after disease challenge suggest they may be potential genetic indicators of disease resilience.