1
|
Liang J, Cade BE, Wang H, Chen H, Gleason KJ, Larkin EK, Saxena R, Lin X, Redline S, Zhu X. Comparison of Heritability Estimation and Linkage Analysis for Multiple Traits Using Principal Component Analyses. Genet Epidemiol 2016; 40:222-32. [PMID: 27027516 DOI: 10.1002/gepi.21957] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2015] [Revised: 11/30/2015] [Accepted: 12/14/2015] [Indexed: 12/16/2022]
Abstract
A disease trait often can be characterized by multiple phenotypic measurements that can provide complementary information on disease etiology, physiology, or clinical manifestations. Given that multiple phenotypes may be correlated and reflect common underlying genetic mechanisms, the use of multivariate analysis of multiple traits may improve statistical power to detect genes and variants underlying complex traits. The literature, however, has been unclear as to the optimal approach for analyzing multiple correlated traits. In this study, heritability and linkage analysis was performed for six obstructive sleep apnea hypopnea syndrome (OSAHS) related phenotypes, as well as principal components of the phenotypes and principal components of the heritability (PCHs) using the data from Cleveland Family Study, which include both African and European American families. Our study demonstrates that principal components generally result in higher heritability and linkage evidence than individual traits. Furthermore, the PCHs can be transferred across populations, strongly suggesting that these PCHs reflect traits with common underlying genetic mechanisms for OSAHS across populations. Thus, PCHs can provide useful traits for using data on multiple phenotypes and for genetic studies of trans-ethnic populations.
Collapse
Affiliation(s)
- Jingjing Liang
- Department of Epidemiology and Biostatistics, School of Medicine, Case Western Reserve University, Cleveland, Ohio, United States of America
| | - Brian E Cade
- Division of Sleep and Circadian Disorders, Brigham and Women's Hospital, Boston, Massachusetts, United States of America.,Division of Sleep Medicine, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Heming Wang
- Department of Epidemiology and Biostatistics, School of Medicine, Case Western Reserve University, Cleveland, Ohio, United States of America
| | - Han Chen
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States of America
| | - Kevin J Gleason
- Division of Sleep and Circadian Disorders, Brigham and Women's Hospital, Boston, Massachusetts, United States of America.,Division of Sleep Medicine, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Emma K Larkin
- Department of Medicine, Division of Allergy, Pulmonary and Critical Care, Vanderbilt University Medical Center, Nashville, Tennessee, United States of America
| | - Richa Saxena
- Division of Sleep and Circadian Disorders, Brigham and Women's Hospital, Boston, Massachusetts, United States of America.,Division of Sleep Medicine, Harvard Medical School, Boston, Massachusetts, United States of America.,Center for Human Genetic Research and Department of Anesthesia, Pain, and Critical Care Medicine, Massachusetts General Hospital, Boston, Massachusetts, United States of America.,Program in Medical and Population Genetics, Broad Institute, Cambridge, Massachusetts, United States of America
| | - Xihong Lin
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States of America
| | - Susan Redline
- Division of Sleep and Circadian Disorders, Brigham and Women's Hospital, Boston, Massachusetts, United States of America.,Division of Sleep Medicine, Harvard Medical School, Boston, Massachusetts, United States of America.,Division of Pulmonary, Critical Care, and Sleep Medicine, Beth Israel Deaconess Medical Center, Boston, Massachusetts, United States of America
| | - Xiaofeng Zhu
- Department of Epidemiology and Biostatistics, School of Medicine, Case Western Reserve University, Cleveland, Ohio, United States of America
| |
Collapse
|
2
|
Chen G, de las Fuentes L, Gu CC, He J, Gu D, Kelly T, Hixson J, Jacquish C, Rao DC, Rice TK. Aggregate blood pressure responses to serial dietary sodium and potassium intervention: defining responses using independent component analysis. BMC Genet 2015; 16:64. [PMID: 26088064 PMCID: PMC4474450 DOI: 10.1186/s12863-015-0226-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2014] [Accepted: 06/10/2015] [Indexed: 11/30/2022] Open
Abstract
BACKGROUND Hypertension is a complex trait that often co-occurs with other conditions such as obesity and is affected by genetic and environmental factors. Aggregate indices such as principal components among these variables and their responses to environmental interventions may represent novel information that is potentially useful for genetic studies. RESULTS In this study of families participating in the Genetic Epidemiology Network of Salt Sensitivity (GenSalt) Study, blood pressure (BP) responses to dietary sodium interventions are explored. Independent component analysis (ICA) was applied to 20 variables indexing obesity and BP measured at baseline and during low sodium, high sodium and high sodium plus potassium dietary intervention periods. A "heat map" protocol that classifies subjects based on risk for hypertension is used to interpret the extracted components. ICA and heat map suggest four components best describe the data: (1) systolic hypertension, (2) general hypertension, (3) response to sodium intervention and (4) obesity. The largest heritabilities are for the systolic (64%) and general hypertension (56%) components. There is a pattern of higher heritability for the component response to intervention (40-42%) as compared to those for the traditional intervention responses computed as delta scores (24%-40%). CONCLUSIONS In summary, the present study provides intermediate phenotypes that are heritable. Using these derived components may prove useful in gene discovery applications.
Collapse
Affiliation(s)
- Gengsheng Chen
- Washington University in St. Louis School of Medicine, St. Louis, MO, USA.
| | | | - Chi C Gu
- Washington University in St. Louis School of Medicine, St. Louis, MO, USA.
| | - Jiang He
- Tulane University Health Sciences Center, New Orleans, LA, USA.
| | - Dongfeng Gu
- Chinese Academy of Medical Sciences, Beijing, China.
| | - Tanika Kelly
- Tulane University Health Sciences Center, New Orleans, LA, USA.
| | - James Hixson
- University of Texas Health Sciences Center at Houston, Houston, TX, USA.
| | | | - D C Rao
- Washington University in St. Louis School of Medicine, St. Louis, MO, USA.
| | - Treva K Rice
- Washington University in St. Louis School of Medicine, St. Louis, MO, USA.
| |
Collapse
|
3
|
Climer S, Yang W, de las Fuentes L, Dávila-Román VG, Gu CC. A custom correlation coefficient (CCC) approach for fast identification of multi-SNP association patterns in genome-wide SNPs data. Genet Epidemiol 2014; 38:610-21. [PMID: 25168954 DOI: 10.1002/gepi.21833] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2014] [Revised: 05/07/2014] [Accepted: 05/19/2014] [Indexed: 01/27/2023]
Abstract
Complex diseases are often associated with sets of multiple interacting genetic factors and possibly with unique sets of the genetic factors in different groups of individuals (genetic heterogeneity). We introduce a novel concept of custom correlation coefficient (CCC) between single nucleotide polymorphisms (SNPs) that address genetic heterogeneity by measuring subset correlations autonomously. It is used to develop a 3-step process to identify candidate multi-SNP patterns: (1) pairwise (SNP-SNP) correlations are computed using CCC; (2) clusters of so-correlated SNPs identified; and (3) frequencies of these clusters in disease cases and controls compared to identify disease-associated multi-SNP patterns. This method identified 42 candidate multi-SNP associations with hypertensive heart disease (HHD), among which one cluster of 22 SNPs (six genes) included 13 in SLC8A1 (aka NCX1, an essential component of cardiac excitation-contraction coupling) and another of 32 SNPs had 29 from a different segment of SLC8A1. While allele frequencies show little difference between cases and controls, the cluster of 22 associated alleles were found in 20% of controls but no cases and the other in 3% of controls but 20% of cases. These suggest that both protective and risk effects on HHD could be exerted by combinations of variants in different regions of SLC8A1, modified by variants from other genes. The results demonstrate that this new correlation metric identifies disease-associated multi-SNP patterns overlooked by commonly used correlation measures. Furthermore, computation time using CCC is a small fraction of that required by other methods, thereby enabling the analyses of large GWAS datasets.
Collapse
Affiliation(s)
- Sharlee Climer
- Division of Biostatistics, Washington University School of Medicine, Missouri, United States of America; Cardiovascular Imaging and Clinical Research Core Laboratory, Cardiovascular Division, Department of Medicine, Washington University School of Medicine, Missouri, United States of America; Current address: Sharlee Climer, Department of Computer Science and Engineering, Washington University School of Engineering, Missouri, United States of America
| | | | | | | | | |
Collapse
|
4
|
|
5
|
Bruchas RR, de las Fuentes L, Carney RM, Reagan JL, Bernal-Mizrachi C, Riek AE, Gu CC, Bierhals A, Schootman M, Malmstrom TK, Burroughs TE, Stein PK, Miller DK, Dávila-Román VG. The St. Louis African American health-heart study: methodology for the study of cardiovascular disease and depression in young-old African Americans. BMC Cardiovasc Disord 2013; 13:66. [PMID: 24011389 PMCID: PMC3847628 DOI: 10.1186/1471-2261-13-66] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2013] [Accepted: 08/13/2013] [Indexed: 01/20/2023] Open
Abstract
BACKGROUND Coronary artery disease (CAD) is a major cause of death and disability worldwide. Depression has complex bidirectional adverse associations with CAD, although the mechanisms mediating these relationships remain unclear. Compared to European Americans, African Americans (AAs) have higher rates of morbidity and mortality from CAD. Although depression is common in AAs, its role in the development and features of CAD in this group has not been well examined. This project hypothesizes that the relationships between depression and CAD can be explained by common physiological pathways and gene-environment interactions. Thus, the primary aims of this ongoing project are to: a) determine the prevalence of CAD and depression phenotypes in a population-based sample of community-dwelling older AAs; b) examine the relationships between CAD and depression phenotypes in this population; and c) evaluate genetic variants from serotoninP and inflammatory pathways to discover potential gene-depression interactions that contribute significantly to the presence of CAD in AAs. METHODS/DESIGN The St. Louis African American Health (AAH) cohort is a population-based panel study of community-dwelling AAs born in 1936-1950 (inclusive) who have been followed from 2000/2001 through 2010. The AAH-Heart study group is a subset of AAH participants recruited in 2009-11 to examine the inter-relationships between depression and CAD in this population. State-of-the-art CAD phenotyping is based on cardiovascular characterizations (coronary artery calcium, carotid intima-media thickness, cardiac structure and function, and autonomic function). Depression phenotyping is based on standardized questionnaires and detailed interviews. Single nucleotide polymorphisms of selected genes in inflammatory and serotonin-signaling pathways are being examined to provide information for investigating potential gene-depression interactions as modifiers of CAD traits. Information from the parent AAH study is being used to provide population-based prevalence estimates. Inflammatory and other biomarkers provide information about potential pathways. DISCUSSION This population-based investigation will provide valuable information on the prevalence of both depression and CAD phenotypes in this population. The study will examine interactions between depression and genetic variants as modulators of CAD, with the intent of detecting mechanistic pathways linking these diseases to identify potential therapeutic targets. Analytic results will be reported as they become available.
Collapse
Affiliation(s)
- Robin R Bruchas
- Cardiovascular Imaging and Clinical Research Core Laboratory, Cardiovascular Division, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8086, St. Louis, MO 63110, USA
- Endocrinology, Metabolism and Lipid Research Division, Washington University School of Medicine, 660 South Euclid Avenue, St. Louis, MO 63110, USA
| | - Lisa de las Fuentes
- Cardiovascular Imaging and Clinical Research Core Laboratory, Cardiovascular Division, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8086, St. Louis, MO 63110, USA
- Division of Biostatistics, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8067, St. Louis, MO 63110, USA
| | - Robert M Carney
- Department of Psychiatry, Washington University School of Medicine, 4320 Forest Park Avenue Suite 301, St. Louis, MO 63108, USA
| | - Joann L Reagan
- Cardiovascular Imaging and Clinical Research Core Laboratory, Cardiovascular Division, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8086, St. Louis, MO 63110, USA
| | - Carlos Bernal-Mizrachi
- Endocrinology, Metabolism and Lipid Research Division, Washington University School of Medicine, 660 South Euclid Avenue, St. Louis, MO 63110, USA
| | - Amy E Riek
- Endocrinology, Metabolism and Lipid Research Division, Washington University School of Medicine, 660 South Euclid Avenue, St. Louis, MO 63110, USA
| | - Chi Charles Gu
- Division of Biostatistics, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8067, St. Louis, MO 63110, USA
| | - Andrew Bierhals
- Mallinckrodt Institute of Radiology, Washington University School of Medicine, 660 South Euclid Avenue, St. Louis, MO 63110, USA
| | - Mario Schootman
- Division of Health Behavior Research, Washington University School of Medicine, 660 south Euclid Avenue, St. Louis, MO 63110, USA
| | - Theodore K Malmstrom
- Department of Neurology & Psychiatry, School of Medicine, Saint Louis University, St. Louis, MO, USA
| | - Thomas E Burroughs
- Center for Outcomes Research, Saint Louis University, 3545 Lafayette Avenue, St. Louis, MO 63104, USA
| | - Phyllis K Stein
- Cardiovascular Imaging and Clinical Research Core Laboratory, Cardiovascular Division, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8086, St. Louis, MO 63110, USA
| | - Douglas K Miller
- Regenstrief Institute, Inc., and Indiana University Center for Aging Research, School of Medicine, Indiana University, 410 West 10th Street, Indianapolis, IN 46202, USA
| | - Victor G Dávila-Román
- Cardiovascular Imaging and Clinical Research Core Laboratory, Cardiovascular Division, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8086, St. Louis, MO 63110, USA
| |
Collapse
|
6
|
Juang JMJ, de las Fuentes L, Waggoner AD, Gu CC, Dávila-Román VG. Association and interaction of PPAR-complex gene variants with latent traits of left ventricular diastolic function. BMC MEDICAL GENETICS 2010; 11:65. [PMID: 20426853 PMCID: PMC2874543 DOI: 10.1186/1471-2350-11-65] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/05/2009] [Accepted: 04/28/2010] [Indexed: 01/04/2023]
Abstract
BACKGROUND Abnormalities in myocardial metabolism and/or regulatory genes have been implicated in left ventricular systolic dysfunction. However, the extent to which these modulate left ventricular diastolic function (LVDF) is uncertain. METHODS Independent component analysis was applied to extract latent LVDF traits from 14 measured echocardiography-derived endophenotypes of LVDF in 403 Caucasians. Genetic association was assessed between measured and latent LVDF traits and 64 single nucleotide polymorphisms (SNPs) in three peroxisome proliferator-activated receptor (PPAR)-complex genes involved in the transcriptional regulation of fatty acid metabolism. RESULTS By linear regression analysis, 7 SNPs (4 in PPARA, 2 in PPARGC1A, 1 in PPARG) were significantly associated with the latent LVDF trait, whereas a range of 0-4 SNPs were associated with each of the 14 measured echocardiography-derived endophenotypes. Frequency distribution of P values showed a greater proportion of significant associations with the latent LVDF trait than for the measured endophenotypes, suggesting that analyses of the latent trait improved detection of the genetic underpinnings of LVDF. Ridge regression was applied to investigate within-gene and gene-gene interactions. In the within-gene analysis, there were five significant pair-wise interactions in PPARGC1A and none in PPARA or PPARG. In the gene-gene analysis, significant interactions were found between rs4253655 in PPARA and rs1873532 (p = 0.02) and rs7672915 (p = 0.02), both in PPARGC1A, and between rs1151996 in PPARG and rs4697046 in PPARGC1A (p = 0.01). CONCLUSIONS Myocardial metabolism PPAR-complex genes, including within and between genes interactions, may play an important role modulating left ventricular diastolic function.
Collapse
Affiliation(s)
- Jyh-Ming Jimmy Juang
- Cardiovascular Imaging and Clinical Research Core Laboratory, Cardiovascular Division, Department of Medicine, Washington University School of Medicine, St. Louis, Missouri, USA
- Division of Biostatistics, Washington University School of Medicine, St. Louis, Missouri, USA
- Division of Cardiology, Department of Internal Medicine, National Taiwan University Hospital and National Taiwan University College of Medicine, Taipei, Taiwan
| | - Lisa de las Fuentes
- Cardiovascular Imaging and Clinical Research Core Laboratory, Cardiovascular Division, Department of Medicine, Washington University School of Medicine, St. Louis, Missouri, USA
| | - Alan D Waggoner
- Cardiovascular Imaging and Clinical Research Core Laboratory, Cardiovascular Division, Department of Medicine, Washington University School of Medicine, St. Louis, Missouri, USA
| | - C Charles Gu
- Division of Biostatistics, Washington University School of Medicine, St. Louis, Missouri, USA
- Department of Genetics, Washington University School of Medicine, St. Louis, Missouri, USA
| | - Víctor G Dávila-Román
- Cardiovascular Imaging and Clinical Research Core Laboratory, Cardiovascular Division, Department of Medicine, Washington University School of Medicine, St. Louis, Missouri, USA
| |
Collapse
|
7
|
Gu CC, Yang W(W, Kraja AT, de las Fuentes L, Dávila-Román VG. Genetic association analysis of coronary heart disease by profiling gene-environment interaction based on latent components in longitudinal endophenotypes. BMC Proc 2009; 3 Suppl 7:S86. [PMID: 20018082 PMCID: PMC2795989 DOI: 10.1186/1753-6561-3-s7-s86] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Studies of complex diseases collect panels of disease-related traits, also known as secondary phenotypes or endophenotypes. They reflect intermediate responses to environment exposures, and as such, are likely to contain hidden information of gene-environment (G x E) interactions. The information can be extracted and used in genetic association studies via latent-components analysis. We present such a method that extracts G x E information in longitudinal data of endophenotypes, and apply the method to repeated measures of multiple phenotypes related to coronary heart disease in Genetic Analysis Workshop 16 Problem 2. The new method identified many genes, including SCNN1B (sodium channel nonvoltage-gated 1 beta) and PKP2 (plakophilin 2), with potential time-dependent G x E interactions; and several others including a novel cardiac-specific kinase gene (TNNI3K), with potential G x E interactions independent of time and marginal effects.
Collapse
Affiliation(s)
- C Charles Gu
- Division of Biostatistics, Washington University School of Medicine, 660 South Euclid Avenue, Box 8067, St. Louis, Missouri 63110, USA
- Department of Genetics, Washington University School of Medicine, 660 South Euclid Avenue, Box 8067, St. Louis, Missouri 63110-1093, USA
| | - Wei (Will) Yang
- Division of Biostatistics, Washington University School of Medicine, 660 South Euclid Avenue, Box 8067, St. Louis, Missouri 63110, USA
| | - Aldi T Kraja
- Division of Statistical Genomics, Department of Genetics, 4444 Forest Park Boulevard, Box 8506, St. Louis, MO 63108, USA
| | - Lisa de las Fuentes
- Cardiovascular Imaging and Clinical Research Core Laboratory, Cardiovascular Division, Department of Medicine, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8086, St. Louis, Missouri 63110, USA
| | - Victor G Dávila-Román
- Cardiovascular Imaging and Clinical Research Core Laboratory, Cardiovascular Division, Department of Medicine, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8086, St. Louis, Missouri 63110, USA
| |
Collapse
|