Metabolomics uses advanced analytical chemistry techniques to comprehensively measure large numbers of small molecule metabolites in cells, tissues and biofluids. The ability to rapidly detect and ...quantify hundreds or even thousands of metabolites within a single sample is helping scientists paint a far more complete picture of system-wide metabolism and biology. Metabolomics is also allowing researchers to focus on measuring the end-products of complex, hard-to-decipher genetic, epigenetic and environmental interactions. As a result, metabolomics has become an increasingly popular "omics" approach to assist with the robust phenotypic characterization of humans, crop plants and model organisms. Indeed, metabolomics is now routinely used in biomedical, nutritional and crop research. It is also being increasingly used in livestock research and livestock monitoring. The purpose of this systematic review is to quantitatively and objectively summarize the current status of livestock metabolomics and to identify emerging trends, preferred technologies and important gaps in the field. In conducting this review we also critically assessed the applications of livestock metabolomics in key areas such as animal health assessment, disease diagnosis, bioproduct characterization and biomarker discovery for highly desirable economic traits (i.e., feed efficiency, growth potential and milk production). A secondary goal of this critical review was to compile data on the known composition of the livestock metabolome (for 5 of the most common livestock species namely cattle, sheep, goats, horses and pigs). These data have been made available through an open access, comprehensive livestock metabolome database (LMDB, available at http://www.lmdb.ca). The LMDB should enable livestock researchers and producers to conduct more targeted metabolomic studies and to identify where further metabolome coverage is needed.
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.
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
), 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.
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.
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
), 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.
Abstract
The objective was to estimate the genetic parameters of performance and resilience of growing pigs under disease. Data were from 3,139 Yorkshire × Landrace wean-to-finish pigs that were ...exposed to a natural polymicrobial disease challenge that was established by entering naturally infected animals into a nursery barn, targeting various viral and bacterial diseases. The challenge was maintained by entering batches of 60 or 75 healthy nursery pigs every 3 wk in a continuous flow system. Traits analyzed included average daily gain (ADG), feed intake (ADFI) and duration (ADFD); feed conversion ratio (FCR); residual feed intake (RFI); mortality (MOR); number of health treatments (TRT); health scores (HScore); carcass weight (CWT), back fat (CBF) and loin depth (CLD); dressing percentage (DRS); lean yield (LYLD); day-to-day variation in feed intake and duration (VARFI and VARDUR); and the proportion of off-feed days (OFFFI and OFFDUR). Analyses were performed by mixed linear models with genomic relationships. The resilience traits, such as TRT, MOR, and HScore, were lowly heritable (up to 0.15) but had high genetic correlations with each other. Performance traits, such as ADG, ADFI, ADFD, FCR, RFI, and carcass traits, were moderate to highly heritable (0.17 to 0.49). Heritabilities of resilience indicator traits such as OFF and VAR had low to moderate heritabilities (0.08 to 0.23) but were higher when based on duration vs. amount. ADFI had a low genetic correlation with ADFD (0.13). ADG in the challenge nursery had stronger negative genetic correlations with both TRT and MOR than ADG in the finisher (−0.37 to −0.74 vs. −0.15 to −0.56). ADFI and FCR had moderate negative (−0.21 to −0.39) and positive (0.34 to 0.49) genetic correlations, respectively, with TRT and MOR. ADFD and RFI had very low genetic correlations with TRT and MOR. CWT and DRS were moderately negatively correlated with TRT and MOR (−0.33 to −0.59). Resilience indicator traits based on feed intake or duration had moderate to high positive genetic correlations with TRT (0.18 to 0.81) and MOR (0.33 to 0.87). In conclusion, performance and resilience traits under a polymicrobial disease challenge are heritable and can be changed by selection. Phenotypes extracted from feed intake patterns can be used as genetic indicator traits for disease resilience. Most promising is day-to-day variation in intake duration, which had a sizeable heritability (0.23) and favorable genetic correlations with MOR (0.79) and treatment rate (0.20).
Disease resilience, defined as an animal's ability to maintain productive performance in the face of infection, provides opportunities to manage the polymicrobial challenge common in pig production. ...Disease resilience can deliver a number of benefits, including more sustainable production as well as improved animal health and the potential for reduced antimicrobial use. However, little progress has been made to date in the application of disease resilience in breeding programs due to a number of factors, including (1) confusion around definitions of disease resilience and its component traits disease resistance and tolerance, and (2) the difficulty in characterizing such a complex trait consisting of multiple biological functions and dynamic elements of rates of response and recovery from infection. Accordingly, this review refines the definitions of disease resistance, tolerance, and resilience based on previous studies to help improve the understanding and application of these breeding goals and traits under different scenarios. We also describe and summarize results from a "natural disease challenge model" designed to provide inputs for selection of disease resilience. The next steps for managing polymicrobial challenges faced by the pig industry will include the development of large-scale multi-omics data, new phenotyping technologies, and mathematical and statistical methods adapted to these data. Genome editing to produce pigs resistant to major diseases may complement selection for disease resilience along with continued efforts in the more traditional areas of biosecurity, vaccination and treatment. Altogether genomic approaches provide exciting opportunities for the pig industry to overcome the challenges provided by hard-to-manage diseases as well as new environmental challenges associated with climate change.
We report the sequencing at 131× coverage, de novo assembly and analyses of the genome of a female Tibetan wild boar. We also resequenced the whole genomes of 30 Tibetan wild boars from six major ...distributed locations and 18 geographically related pigs in China. We characterized genetic diversity, population structure and patterns of evolution. We searched for genomic regions under selection, which includes genes that are involved in hypoxia, olfaction, energy metabolism and drug response. Comparing the genome of Tibetan wild boar with those of neighboring Chinese domestic pigs further showed the impact of thousands of years of artificial selection and different signatures of selection in wild boar and domestic pig. We also report genetic adaptations in Tibetan wild boar that are associated with high altitudes and characterize the genetic basis of increased salivation in domestic pig.
Porcine circovirus 2 (PCV2) is a circular single-stranded DNA virus responsible for a group of diseases collectively known as PCV2 Associated Diseases (PCVAD). Variation in the incidence and severity ...of PCVAD exists between pigs suggesting a host genetic component involved in pathogenesis. A large-scale genome-wide association study of experimentally infected pigs (n = 974), provided evidence of a host genetic role in PCV2 viremia, immune response and growth during challenge. Host genotype explained 64% of the phenotypic variation for overall viral load, with two major Quantitative Trait Loci (QTL) identified on chromosome 7 (SSC7) near the swine leukocyte antigen complex class II locus and on the proximal end of chromosome 12 (SSC12). The SNP having the strongest association, ALGA0110477 (SSC12), explained 9.3% of the genetic and 6.2% of the phenotypic variance for viral load. Dissection of the SSC12 QTL based on gene annotation, genomic and RNA-sequencing, suggested that a missense mutation in the SYNGR2 (SYNGR2 p.Arg63Cys) gene is potentially responsible for the variation in viremia. This polymorphism, located within a protein domain conserved across mammals, results in an amino acid variant SYNGR2 p.63Cys only observed in swine. PCV2 titer in PK15 cells decreased when the expression of SYNGR2 was silenced by specific-siRNA, indicating a role of SYNGR2 in viral replication. Additionally, a PK15 edited clone generated by CRISPR-Cas9, carrying a partial deletion of the second exon that harbors a key domain and the SYNGR2 p.Arg63Cys, was associated with a lower viral titer compared to wildtype PK15 cells (>24 hpi) and supernatant (>48hpi)(P < 0.05). Identification of a non-conservative substitution in this key domain of SYNGR2 suggests that the SYNGR2 p.Arg63Cys variant may underlie the observed genetic effect on viral load.
Campylobacter jejuni
is a common cause of diarrheal disease worldwide. Human infection typically occurs through the ingestion of contaminated poultry products. We previously demonstrated that an ...attenuated
Escherichia coli
live vaccine strain expressing the
C. jejuni
N-glycan on its surface reduced the
Campylobacter
load in more than 50% of vaccinated leghorn and broiler birds to undetectable levels (responder birds), whereas the remainder of the animals was still colonized (non-responders). To understand the underlying mechanism, we conducted three vaccination and challenge studies using 135 broiler birds and found a similar responder/non-responder effect. Subsequent genome-wide association studies (GWAS), analyses of bird sex and levels of vaccine-induced IgY responses did not correlate with the responder versus non-responder phenotype. In contrast, antibodies isolated from responder birds displayed a higher
Campylobacter
-opsonophagocytic activity when compared to antisera from non-responder birds. No differences in the N-glycome of the sera could be detected, although minor changes in IgY glycosylation warrant further investigation. As reported before, the composition of the microbiota, particularly levels of OTU classified as
Clostridium
spp.,
Ruminococcaceae
and
Lachnospiraceae
are associated with the response. Transplantation of the cecal microbiota of responder birds into new birds in combination with vaccination resulted in further increases in vaccine-induced antigen-specific IgY responses when compared to birds that did not receive microbiota transplants. Our work suggests that the IgY effector function and microbiota contribute to the efficacy of the
E. coli
live vaccine, information that could form the basis for the development of improved vaccines targeted at the elimination of
C. jejuni
from poultry.
Bovine respiratory disease (BRD) is the most common and costly infectious disease affecting the wellbeing and productivity of beef cattle in North America. BRD is a complex disease whose development ...is dependent on environmental factors and host genetics. Due to the polymicrobial nature of BRD, our understanding of the genetic and molecular mechanisms underlying the disease is still limited. This knowledge would augment the development of better genetic/genomic selection strategies and more accurate diagnostic tools to reduce BRD prevalence. Therefore, this study aimed to utilize multi-omics data (genomics, transcriptomics, and metabolomics) analyses to study the genetic and molecular mechanisms of BRD infection. Blood samples of 143 cattle (80 BRD; 63 non-BRD animals) were collected for genotyping, RNA sequencing, and metabolite profiling. Firstly, a genome-wide association study (GWAS) was performed for BRD susceptibility using 207,038 SNPs. Two SNPs (Chr5:25858264 and BovineHD1800016801) were identified as associated (
-value <1 × 10
) with BRD susceptibility. Secondly, differential gene expression between BRD and non-BRD animals was studied. At the significance threshold used (log
FC>2, logCPM>2, and FDR<0.01), 101 differentially expressed (DE) genes were identified. These DE genes significantly (
-value <0.05) enriched several immune responses related functions such as inflammatory response. Additionally, we performed expression quantitative trait loci (eQTL) analysis and identified 420 cis-eQTLs and 144 trans-eQTLs significantly (FDR <0.05) associated with the expression of DE genes. Interestingly, eQTL results indicated the most significant SNP (Chr5:25858264) identified
GWAS was a cis-eQTL for DE gene
This analysis also demonstrated that an important SNP (rs209419196) located in the promoter region of the DE gene
significantly influenced the expression of this gene. Finally, the abundance of 31 metabolites was significantly (FDR <0.05) different between BRD and non-BRD animals, and 17 of them showed correlations with multiple DE genes, which shed light on the interactions between immune response and metabolism. This study identified associations between genome, transcriptome, metabolome, and BRD phenotype of feedlot crossbred cattle. The findings may be useful for the development of genomic selection strategies for BRD susceptibility, and for the development of new diagnostic and therapeutic tools.
It is possible to identify sub-populations of sows in every pig herd that consistently give birth to low birth weight (BW) piglets, irrespective of the litter size. A previous study from our group ...demonstrated that placental development is a main factor affecting the litter birth weight phenotype (LBWP) in sows, thereby impacting the BW of entire litters, but the biological and molecular pathways behind this phenomenon are largely unknown. The aim of this study was to investigate the differential gene expression in placental tissues at day 30 of gestation between low LBWP (LLBWP) vs. high LBWP (HLBWP) sows from a purebred Large White maternal line. Using mRNA sequencing, we found 45 differentially expressed genes (DEGs) in placental tissues of LLBWP and HLBWP sows. Furthermore, (GO) enrichment of upregulated DEGs predicted that there were two biological processes significantly related to cornification and regulation of cell population proliferation. To better understand the molecular interaction between cell proliferation and cornification, we conducted transcriptional factor binding site (TFBS) prediction analysis. The results indicated that a highly significant TFBS was located at the 5′ upstream of all four upregulated genes (CDSN, DSG3, KLK14, KRT17), recognized by transcription factors EGR4 and FOSL1. Our findings provide novel insight into how transcriptional regulation of two different biological processes interact in placental tissues of LLBWP sows.