Combined genetic, genomic and transcriptomic methods in the analysis of animal traits

Multi-omic data integration and analysis using systems genomics approaches: methods and applications in animal production, health and welfare

In the past years, there has been a remarkable development of high-throughput omics (HTO) technologies such as genomics, epigenomics, transcriptomics, proteomics and metabolomics across all facets of biology. This has spearheaded the progress of the systems biology era, including applications on animal production and health traits. However, notwithstanding these new HTO technologies, there remains an emerging challenge in data analysis. On the one hand, different HTO technologies judged on their own merit are appropriate for the identification of disease-causing genes, biomarkers for prevention and drug targets for the treatment of diseases and for individualized genomic predictions of performance or disease risks. On the other hand, integration of multi-omic data and joint modelling and analyses are very powerful and accurate to understand the systems biology of healthy and sustainable production of animals. We present an overview of current and emerging HTO technologies each with a focus on their applications in animal and veterinary sciences before introducing an integrative systems genomics framework for analysing and integrating multi-omic data towards improved animal production, health and welfare. We conclude that there are big challenges in multi-omic data integration, modelling and systems-level analyses, particularly with the fast emerging HTO technologies. We highlight existing and emerging systems genomics approaches and discuss how they contribute to our understanding of the biology of complex traits or diseases and holistic improvement of production performance, disease resistance and welfare.

An integrative systems genetics approach reveals potential causal genes and pathways related to obesity

Background

Obesity is a multi-factorial health problem in which genetic factors play an important role. Limited results have been obtained in single-gene studies using either genomic or transcriptomic data. RNA sequencing technology has shown its potential in gaining accurate knowledge about the transcriptome, and may reveal novel genes affecting complex diseases. Integration of genomic and transcriptomic variation (expression quantitative trait loci [eQTL] mapping) has identified causal variants that affect complex diseases. We integrated transcriptomic data from adipose tissue and genomic data from a porcine model to investigate the mechanisms involved in obesity using a systems genetics approach.

Methods

Using a selective gene expression profiling approach, we selected 36 animals based on a previously created genomic Obesity Index for RNA sequencing of subcutaneous adipose tissue. Differential expression analysis was performed using the Obesity Index as a continuous variable in a linear model. eQTL mapping was then performed to integrate 60 K porcine SNP chip data with the RNA sequencing data. Results were restricted based on genome-wide significant single nucleotide polymorphisms, detected differentially expressed genes, and previously detected co-expressed gene modules. Further data integration was performed by detecting co-expression patterns among eQTLs and integration with protein data.

Results

Differential expression analysis of RNA sequencing data revealed 458 differentially expressed genes. The eQTL mapping resulted in 987 cis-eQTLs and 73 trans-eQTLs (false discovery rate < 0.05), of which the cis-eQTLs were associated with metabolic pathways. We reduced the eQTL search space by focusing on differentially expressed and co-expressed genes and disease-associated single nucleotide polymorphisms to detect obesity-related genes and pathways. Building a co-expression network using eQTLs resulted in the detection of a module strongly associated with lipid pathways. Furthermore, we detected several obesity candidate genes, for example, ENPP1, CTSL, and ABHD12B.

Conclusions

To our knowledge, this is the first study to perform an integrated genomics and transcriptomics (eQTL) study using, and modeling, genomic and subcutaneous adipose tissue RNA sequencing data on obesity in a porcine model. We detected several pathways and potential causal genes for obesity. Further validation and investigation may reveal their exact function and association with obesity.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-015-0229-0) contains supplementary material, which is available to authorized users.

Systems biology of ovine intestinal parasite resistance: disease gene modules and biomarkers

This study reports on the molecular systems biology of gastrointestinal nematode (GIN) infection and potential biomarkers for GIN resistance in sheep. Microarray gene expression data were obtained for 3 different tissues at 4 time points from sheep artificially challenged with two types of nematodes, Haemonchus contortus (HC) and Trichostrongylus colubriformis (TC). We employed an integrated systems biology approach, integrating 3 main methods: standard differential gene expression analyses, weighted gene co-expression network analyses (WGCNA) and quantitative genetic analyses of gene expression traits of key biomarkers. Using standard differential gene expression analyses we identified differentially expressed genes (DE) which responded differently in sheep challenged with HC compared to those challenged with TC. These interaction genes (e.g. MRPL51, SMEK2, CAT, MAPK1IP1 and SLC25A20A) were enriched in Wnt receptor signalling pathway (p = 0.0132) and positive regulation of NFκβ transcription factor activity (p = 0.00208). We report FCER1A, a gene encoding a high-affinity receptor for the Fc region of immunoglobulin E, which is linked to innate immunity to GIN in sheep. Using weighted gene co-expression network analysis (WGCNA) methods, we identified gene modules that were correlated with the length of infection (disease modules). Hub genes (with high intramodular connectivity) were filtered further to identify biomarkers that are related to the length of infection (e.g. CAT, FBX033, COL15A1, IGFBP7, FBLN1 and IgCgamma). The biomarkers we found in HC networks were significantly associated with functions such as T-cell and B-cell regulations, TNF-alpha, interleukin and cytokine production. In TC networks, biomarkers were significantly associated with functions such as protein catabolic process, heat shock protein binding, protein targeting and localization, cytokine receptor binding, TNF receptor binding, apoptosis and IGF binding. These results provide specific gene targets for therapeutic interventions and provide insights into GIN infections in sheep which may be used to infer the same in related host species. This is also the first study to apply the concept of estimating breeding values of animals to expression traits and reveals 11 heritable candidate biomarkers (0.05 to 0.92) that could be used in selection of animals for GIN resistance.

Genetical systems biology in livestock: application to gonadotrophin releasing hormone and reproduction

Genetical systems biology or systems genetics treats the genome as the central reference point for all omics variations and is an emerging new branch of systems biology. Quantitative genetic principles were developed for high-throughput genomic, transcriptomic and metabolomic data observed in large populations. New statistical genetic models were developed for expression quantitative trait loci (eQTL), namely, marker regression eQTL mapping and marker-expression co-factor mapping. Evaluations of power to detect eQTL showed that sample size requirements are higher for detecting trans-acting genes than cis-acting genes. Power is higher for eQTL with high heritability than for eQTL with low heritability. These results will be valuable for systems genetic investigations. Gonadotrophin releasing hormone (GnRH) and its receptor gene (GnRH-R) are crucial for mammalian reproduction. Whole genome scan of eQTLs for GnRH-R gene expression in mouse showed three possible trans-eQTL regions on chr 13 and 19, harbouring regulatory genes. Applications of genetical genomics in systems biology were identified as: (1) detection and validation of causal gene for complex traits; (2) development of genetic interaction networks; (3) prediction of transcription factor binding sites and (4) in data-driven systems biology. These applications were illustrated using data on eQTL, protein network and signalling pathways for GnRH. Gpr54 (G protein-coupled receptor kinase 54), Prl (prolactin), Ins1 (insulin) and Fos (viral oncogenes) were found to be major regulators of GnRH and GnRH-R; thus validating their important role in reproduction, mammary gland development and sexual (im)maturity. These results will be useful for further study of mammalian reproductive biology.