Novel Loci for Metabolic Networks and Multi-Tissue Expression Studies Reveal Genes for Atherosclerosis

Association testing of multiple correlated phenotypes offers better power than univariate analysis of single traits. We analyzed 6,600 individuals from two population-based cohorts with both genome-wide SNP data and serum metabolomic profiles. From the observed correlation structure of 130 metabolites measured by nuclear magnetic resonance, we identified 11 metabolic networks and performed a multivariate genome-wide association analysis. We identified 34 genomic loci at genome-wide significance, of which 7 are novel. In comparison to univariate tests, multivariate association analysis identified nearly twice as many significant associations in total. Multi-tissue gene expression studies identified variants in our top loci, SERPINA1 and AQP9, as eQTLs and showed that SERPINA1 and AQP9 expression in human blood was associated with metabolites from their corresponding metabolic networks. Finally, liver expression of AQP9 was associated with atherosclerotic lesion area in mice, and in human arterial tissue both SERPINA1 and AQP9 were shown to be upregulated (6.3-fold and 4.6-fold, respectively) in atherosclerotic plaques. Our study illustrates the power of multi-phenotype GWAS and highlights candidate genes for atherosclerosis.

[1]  G. Abecasis,et al.  Genotype imputation. , 2009, Annual review of genomics and human genetics.

[2]  L. Peltonen,et al.  Genome scans provide evidence for low-HDL-C loci on chromosomes 8q23, 16q24.1-24.2, and 20q13.11 in Finnish families. , 2002, American journal of human genetics.

[3]  E. Xing,et al.  Statistical Estimation of Correlated Genome Associations to a Quantitative Trait Network , 2009, PLoS genetics.

[4]  P. Bork,et al.  A method and server for predicting damaging missense mutations , 2010, Nature Methods.

[5]  Y. Shibamoto,et al.  The ZFHX3 (ATBF1) transcription factor induces PDGFRB, which activates ATM in the cytoplasm to protect cerebellar neurons from oxidative stress , 2010, Disease Models & Mechanisms.

[6]  P. Donnelly,et al.  A Flexible and Accurate Genotype Imputation Method for the Next Generation of Genome-Wide Association Studies , 2009, PLoS genetics.

[7]  Jean Tichet,et al.  A Polymorphism Within the G6PC2 Gene Is Associated with Fasting Plasma Glucose Levels , 2008, Science.

[8]  T. Joh,et al.  ATBF1 enhances the suppression of STAT3 signaling by interaction with PIAS3. , 2004, Biochemical and biophysical research communications.

[9]  Nicholas J. Schork,et al.  Coexpression network analysis of neural tissue reveals perturbations in developmental processes in schizophrenia. , 2010, Genome research.

[10]  R. Flavell,et al.  Disruption of Myosin 1e promotes podocyte injury. , 2009, Journal of the American Society of Nephrology : JASN.

[11]  Bin Zhang,et al.  Defining clusters from a hierarchical cluster tree: the Dynamic Tree Cut package for R , 2008, Bioinform..

[12]  Tero Aittokallio,et al.  Genetic Variants and Their Interactions in the Prediction of Increased Pre-Clinical Carotid Atherosclerosis: The Cardiovascular Risk in Young Finns Study , 2010, PLoS genetics.

[13]  V. Mootha,et al.  Metabolite profiles and the risk of developing diabetes , 2011, Nature Medicine.

[14]  Kasper Lage,et al.  Pervasive Sharing of Genetic Effects in Autoimmune Disease , 2011, PLoS genetics.

[15]  T. Funahashi,et al.  Metabolic impact of adipose and hepatic glycerol channels aquaporin 7 and aquaporin 9 , 2008, Nature Clinical Practice Endocrinology &Metabolism.

[16]  Inês Barroso,et al.  Genetic Variants Influencing Circulating Lipid Levels and Risk of Coronary Artery Disease , 2010, Arteriosclerosis, thrombosis, and vascular biology.

[17]  P. Libby,et al.  Progress and challenges in translating the biology of atherosclerosis , 2011, Nature.

[18]  M. Daly,et al.  Estimation of the multiple testing burden for genomewide association studies of nearly all common variants , 2008, Genetic epidemiology.

[19]  宁北芳,et al.  疟原虫var基因转换速率变化导致抗原变异[英]／Paul H, Robert P, Christodoulou Z, et al//Proc Natl Acad Sci U S A , 2005 .

[20]  H. Hotelling.  Relations Between Two Sets of Variates , 1936 .

[21]  Risto Telama,et al.  Cohort profile: the cardiovascular risk in Young Finns Study. , 2008, International journal of epidemiology.

[22]  B. L. de Groot,et al.  Aquaporin-9 Protein Is the Primary Route of Hepatocyte Glycerol Uptake for Glycerol Gluconeogenesis in Mice* , 2011, The Journal of Biological Chemistry.

[23]  Markus Perola,et al.  Metabonomic, transcriptomic, and genomic variation of a population cohort , 2010, Molecular systems biology.

[24]  Yurii S. Aulchenko,et al.  Genome-Wide Association Study Identifies Novel Loci Associated with Circulating Phospho- and Sphingolipid Concentrations , 2012, PLoS genetics.

[25]  Anbupalam Thalamuthu,et al.  A Genome-Wide Association Study Identifies Novel and Functionally Related Susceptibility Loci for Kawasaki Disease , 2009, PLoS genetics.

[26]  C. Gieger,et al.  Human metabolic individuality in biomedical and pharmaceutical research , 2011, Nature.

[27]  P. Rantakallio,et al.  Groups at risk in low birth weight infants and perinatal mortality. , 1969, Acta paediatrica Scandinavica.

[28]  A. Arnold,et al.  Tissue-specific expression and regulation of sexually dimorphic genes in mice. , 2006, Genome research.

[29]  A. Burlingame,et al.  HDL antielastase activity prevents smooth muscle cell anoikis, a potential new antiatherogenic property , 2009, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[30]  Michael Karin,et al.  Ultraviolet Light and Osmotic Stress: Activation of the JNK Cascade Through Multiple Growth Factor and Cytokine Receptors , 1996, Science.

[31]  E. Ostap,et al.  Myo1e binds anionic phospholipids with high affinity. , 2010, Biochemistry.

[32]  Christian Gieger,et al.  Genetics Meets Metabolomics: A Genome-Wide Association Study of Metabolite Profiles in Human Serum , 2008, PLoS genetics.

[33]  Eric E Schadt,et al.  Identification of Pathways for Atherosclerosis in Mice: Integration of Quantitative Trait Locus Analysis and Global Gene Expression Data , 2007, Circulation research.

[34]  E. Renard,et al.  CXC ligand 5 is an adipose-tissue derived factor that links obesity to insulin resistance. , 2009, Cell metabolism.

[35]  L. Peltonen,et al.  Familial combined hyperlipidemia is associated with upstream transcription factor 1 (USF1) , 2004, Nature Genetics.

[36]  F. Agakov,et al.  Abundant pleiotropy in human complex diseases and traits. , 2011, American journal of human genetics.

[37]  C. Hoggart,et al.  Genome-wide association analysis of metabolic traits in a birth cohort from a founder population , 2008, Nature Genetics.

[38]  Manuel A. R. Ferreira,et al.  A multivariate test of association , 2009, Bioinform..

[39]  G. Remuzzi,et al.  MYO1E mutations and childhood familial focal segmental glomerulosclerosis. , 2011, The New England journal of medicine.

[40]  D B Allison,et al.  Multiple phenotype modeling in gene-mapping studies of quantitative traits: power advantages. , 1998, American journal of human genetics.

[41]  Reino Laatikainen,et al.  High-throughput serum NMR metabonomics for cost-effective holistic studies on systemic metabolism. , 2009, The Analyst.

[42]  F. Collins,et al.  Potential etiologic and functional implications of genome-wide association loci for human diseases and traits , 2009, Proceedings of the National Academy of Sciences.

[43]  B. Yandell,et al.  Bayesian Quantitative Trait Loci Mapping for Multiple Traits , 2008, Genetics.

[44]  Marcia M. Nizzari,et al.  Genome-Wide Association Analysis Identifies Loci for Type 2 Diabetes and Triglyceride Levels , 2007, Science.

[45]  John D. Storey,et al.  Mapping the Genetic Architecture of Gene Expression in Human Liver , 2008, PLoS biology.

[46]  Tanya M. Teslovich,et al.  Biological, Clinical, and Population Relevance of 95 Loci for Blood Lipids , 2010, Nature.

[47]  Niku Oksala,et al.  ADAM-9, ADAM-15, and ADAM-17 are upregulated in macrophages in advanced human atherosclerotic plaques in aorta and carotid and femoral arteries—Tampere vascular study , 2009, Annals of medicine.

[48]  C. Gieger,et al.  A sequence variant in ZFHX3 on 16q22 associates with atrial fibrillation and ischemic stroke , 2009, Nature Genetics.

[49]  Markus Perola,et al.  An Immune Response Network Associated with Blood Lipid Levels , 2010, PLoS genetics.

[50]  Zhaohui S. Qin,et al.  A second generation human haplotype map of over 3.1 million SNPs , 2007, Nature.

[51]  Christian Gieger,et al.  A genome-wide perspective of genetic variation in human metabolism , 2010, Nature Genetics.

[52]  Manuel A. R. Ferreira,et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. , 2007, American journal of human genetics.

[53]  E. Niki,et al.  In vivo complex formation of oxidized alpha(1)-antitrypsin and LDL. , 2001, Arteriosclerosis, thrombosis, and vascular biology.

[54]  Xia Yang,et al.  Systematic genetic and genomic analysis of cytochrome P450 enzyme activities in human liver. , 2010, Genome research.

[55]  W D Wagner,et al.  A definition of advanced types of atherosclerotic lesions and a histological classification of atherosclerosis. A report from the Committee on Vascular Lesions of the Council on Arteriosclerosis, American Heart Association. , 1995, Arteriosclerosis, thrombosis, and vascular biology.

[56]  A. Dobson,et al.  Estimation of contribution of changes in classic risk factors to trends in coronary-event rates across the WHO MONICA Project populations , 2000, The Lancet.

[57]  Markus Perola,et al.  Genome-wide association study identifies multiple loci influencing human serum metabolite levels , 2012, Nature Genetics.

[58]  G. Calamita,et al.  Insulin- and leptin-mediated control of aquaglyceroporins in human adipocytes and hepatocytes is mediated via the PI3K/Akt/mTOR signaling cascade. , 2011, The Journal of clinical endocrinology and metabolism.

[59]  E. Füchtbauer,et al.  Defective glycerol metabolism in aquaporin 9 (AQP9) knockout mice , 2007, Proceedings of the National Academy of Sciences.

[60]  T. Kodama,et al.  In Vivo Complex Formation of Oxidized α1-Antitrypsin and LDL , 2001 .

[61]  R. Eils,et al.  Comparison of normalization methods for Illumina BeadChip HumanHT-12 v3 , 2010, BMC Genomics.