TY - JOUR

T1 - A network-based conditional genetic association analysis of the human metabolome

AU - Tsepilov, Y. A.

AU - Sharapov, S. Z.

AU - Zaytseva, O. O.

AU - Krumsek, J.

AU - Prehn, C.

AU - Adamski, J.

AU - Kastenmüller, G.

AU - Wang-Sattler, R.

AU - Strauch, K.

AU - Gieger, C.

AU - Aulchenko, Y. S.

N1 - Publisher Copyright: © The Author(s) 2018.

PY - 2018/11/29

Y1 - 2018/11/29

N2 - Background: Genome-wide association studies have identified hundreds of loci that influence a wide variety of complex human traits; however, little is known regarding the biological mechanism of action of these loci. The recent accumulation of functional genomics ("omics"), including metabolomics data, has created new opportunities for studying the functional role of specific changes in the genome. Functional genomic data are characterized by their high dimensionality, the presence of (strong) statistical dependency between traits, and, potentially, complex genetic control. Therefore, the analysis of such data requires specific statistical genetics methods. Results: To facilitate our understanding of the genetic control of omics phenotypes, we propose a trait-centered, network-based conditional genetic association (cGAS) approach for identifying the direct effects of genetic variants on omics-based traits. For each trait of interest, we selected from a biological network a set of other traits to be used as covariates in the cGAS. The network can be reconstructed either from biological pathway databases (a mechanistic approach) or directly from the data, using a Gaussian graphical model applied to the metabolome (a data-driven approach). We derived mathematical expressions that allow comparison of the power of univariate analyses with conditional genetic association analyses. We then tested our approach using data from a population-based Cooperative Health Research in the region of Augsburg (KORA) study (n = 1,784 subjects, 1.7 million single-nucleotide polymorphisms) with measured data for 151 metabolites. Conclusions: We found that compared to single-trait analysis, performing a genetic association analysis that includes biologically relevant covariates can either gain or lose power, depending on specific pleiotropic scenarios, for which we provide empirical examples. In the context of analyzed metabolomics data, the mechanistic network approach had more power compared to the data-driven approach. Nevertheless, we believe that our analysis shows that neither a prior-knowledge-only approach nor a phenotypic-data-only approach is optimal, and we discuss possibilities for improvement.

AB - Background: Genome-wide association studies have identified hundreds of loci that influence a wide variety of complex human traits; however, little is known regarding the biological mechanism of action of these loci. The recent accumulation of functional genomics ("omics"), including metabolomics data, has created new opportunities for studying the functional role of specific changes in the genome. Functional genomic data are characterized by their high dimensionality, the presence of (strong) statistical dependency between traits, and, potentially, complex genetic control. Therefore, the analysis of such data requires specific statistical genetics methods. Results: To facilitate our understanding of the genetic control of omics phenotypes, we propose a trait-centered, network-based conditional genetic association (cGAS) approach for identifying the direct effects of genetic variants on omics-based traits. For each trait of interest, we selected from a biological network a set of other traits to be used as covariates in the cGAS. The network can be reconstructed either from biological pathway databases (a mechanistic approach) or directly from the data, using a Gaussian graphical model applied to the metabolome (a data-driven approach). We derived mathematical expressions that allow comparison of the power of univariate analyses with conditional genetic association analyses. We then tested our approach using data from a population-based Cooperative Health Research in the region of Augsburg (KORA) study (n = 1,784 subjects, 1.7 million single-nucleotide polymorphisms) with measured data for 151 metabolites. Conclusions: We found that compared to single-trait analysis, performing a genetic association analysis that includes biologically relevant covariates can either gain or lose power, depending on specific pleiotropic scenarios, for which we provide empirical examples. In the context of analyzed metabolomics data, the mechanistic network approach had more power compared to the data-driven approach. Nevertheless, we believe that our analysis shows that neither a prior-knowledge-only approach nor a phenotypic-data-only approach is optimal, and we discuss possibilities for improvement.

KW - genome-wide association study

KW - multivariate model

KW - metabolomics

KW - conditional analysis

KW - pleiotropy

KW - GENOME-WIDE ASSOCIATION

KW - GWAS

KW - ATLAS

KW - Genome-Wide Association Study

KW - Humans

KW - Genotype

KW - Metabolome/genetics

KW - Genetic Loci

KW - Metabolic Networks and Pathways/genetics

KW - Phenotype

KW - Algorithms

KW - Polymorphism, Single Nucleotide

KW - Metabolomics/methods

KW - Metabolomics

KW - Pleiotropy

KW - Genome-wide association study

KW - Conditional analysis

KW - Multivariate model

UR - http://www.scopus.com/inward/record.url?scp=85058610423&partnerID=8YFLogxK

U2 - 10.1093/gigascience/giy137

DO - 10.1093/gigascience/giy137

M3 - Article

C2 - 30496450

AN - SCOPUS:85058610423

VL - 7

JO - GigaScience

JF - GigaScience

SN - 2047-217X

IS - 12

M1 - 137

ER -