Alex P. Reiner

A Genome-Wide Association Meta-Analysis of Circulating Sex Hormone–Binding Globulin Reveals Multiple Loci Implicated in Sex Steroid Hormone Regulation

Sex hormone-binding globulin (SHBG) is a glycoprotein responsible for the transport and biologic availability of sex steroid hormones, primarily testosterone and estradiol. SHBG has been associated with chronic diseases including type 2 diabetes (T2D) and with hormone-sensitive cancers such as breast and prostate cancer. We performed a genome-wide association study (GWAS) meta-analysis of 21,791 individuals from 10 epidemiologic studies and validated these findings in 7,046 individuals in an additional six studies. We identified twelve genomic regions (SNPs) associated with circulating SHBG concentrations. Loci near the identified SNPs included SHBG (rs12150660, 17p13.1, pu200a=u200a1.8×10−106), PRMT6 (rs17496332, 1p13.3, pu200a=u200a1.4×10−11), GCKR (rs780093, 2p23.3, pu200a=u200a2.2×10−16), ZBTB10 (rs440837, 8q21.13, pu200a=u200a3.4×10−09), JMJD1C (rs7910927, 10q21.3, pu200a=u200a6.1×10−35), SLCO1B1 (rs4149056, 12p12.1, pu200a=u200a1.9×10−08), NR2F2 (rs8023580, 15q26.2, pu200a=u200a8.3×10−12), ZNF652 (rs2411984, 17q21.32, pu200a=u200a3.5×10−14), TDGF3 (rs1573036, Xq22.3, pu200a=u200a4.1×10−14), LHCGR (rs10454142, 2p16.3, pu200a=u200a1.3×10−07), BAIAP2L1 (rs3779195, 7q21.3, pu200a=u200a2.7×10−08), and UGT2B15 (rs293428, 4q13.2, pu200a=u200a5.5×10−06). These genes encompass multiple biologic pathways, including hepatic function, lipid metabolism, carbohydrate metabolism and T2D, androgen and estrogen receptor function, epigenetic effects, and the biology of sex steroid hormone-responsive cancers including breast and prostate cancer. We found evidence of sex-differentiated genetic influences on SHBG. In a sex-specific GWAS, the loci 4q13.2-UGT2B15 was significant in men only (men pu200a=u200a2.5×10−08, women pu200a=u200a0.66, heterogeneity pu200a=u200a0.003). Additionally, three loci showed strong sex-differentiated effects: 17p13.1-SHBG and Xq22.3-TDGF3 were stronger in men, whereas 8q21.12-ZBTB10 was stronger in women. Conditional analyses identified additional signals at the SHBG gene that together almost double the proportion of variance explained at the locus. Using an independent study of 1,129 individuals, all SNPs identified in the overall or sex-differentiated or conditional analyses explained ∼15.6% and ∼8.4% of the genetic variation of SHBG concentrations in men and women, respectively. The evidence for sex-differentiated effects and allelic heterogeneity highlight the importance of considering these features when estimating complex trait variance.

Imputation of Exome Sequence Variants into Population- Based Samples and Blood-Cell-Trait-Associated Loci in African Americans: NHLBI GO Exome Sequencing Project

Researchers have successfully applied exome sequencing to discover causal variants in selected individuals with familial, highly penetrant disorders. We demonstrate the utility of exome sequencing followed by imputation for discovering low-frequency variants associated with complex quantitative traits. We performed exome sequencing in a reference panel of 761 African Americans and then imputed newly discovered variants into a larger sample of more than 13,000 African Americans for association testing with the blood cell traits hemoglobin, hematocrit, white blood count, and platelet count. First, we illustrate the feasibility of our approach by demonstrating genome-wide-significant associations for variants that are not covered by conventional genotyping arrays; for example, one such association is that between higher platelet count and an MPL c.117G>T (p.Lys39Asn) variant encoding a p.Lys39Asn amino acid substitution of the thrombopoietin receptor gene (p = 1.5 × 10(-11)). Second, we identified an association between missense variants of LCT and higher white blood count (p = 4 × 10(-13)). Third, we identified low-frequency coding variants that might account for allelic heterogeneity at several known blood cell-associated loci: MPL c.754T>C (p.Tyr252His) was associated with higher platelet count; CD36 c.975T>G (p.Tyr325(∗)) was associated with lower platelet count; and several missense variants at the α-globin gene locus were associated with lower hemoglobin. By identifying low-frequency missense variants associated with blood cell traits not previously reported by genome-wide association studies, we establish that exome sequencing followed by imputation is a powerful approach to dissecting complex, genetically heterogeneous traits in large population-based studies.

Combined admixture mapping and association analysis identifies a novel blood pressure genetic locus on 5p13: contributions from the CARe consortium

Admixture mapping based on recently admixed populations is a powerful method to detect disease variants with substantial allele frequency differences in ancestral populations. We performed admixture mapping analysis for systolic blood pressure (SBP) and diastolic blood pressure (DBP), followed by trait-marker association analysis, in 6303 unrelated African-American participants of the Candidate Gene Association Resource (CARe) consortium. We identified five genomic regions (P< 0.001) harboring genetic variants contributing to inter-individual BP variation. In follow-up association analyses, correcting for all tests performed in this study, three loci were significantly associated with SBP and one significantly associated with DBP (P< 10(-5)). Further analyses suggested that six independent single-nucleotide polymorphisms (SNPs) contributed to the phenotypic variation observed in the admixture mapping analysis. These six SNPs were examined for replication in multiple, large, independent studies of African-Americans [Womens Health Initiative (WHI), Maywood, Genetic Epidemiology Network of Arteriopathy (GENOA) and Howard University Family Study (HUFS)] as well as one native African sample (Nigerian study), with a total replication sample size of 11 882. Meta-analysis of the replication set identified a novel variant (rs7726475) on chromosome 5 between the SUB1 and NPR3 genes, as being associated with SBP and DBP (P< 0.0015 for both); in meta-analyses combining the CARe samples with the replication data, we observed P-values of 4.45 × 10(-7) for SBP and 7.52 × 10(-7) for DBP for rs7726475 that were significant after accounting for all the tests performed. Our study highlights that admixture mapping analysis can help identify genetic variants missed by genome-wide association studies because of drastically reduced number of tests in the whole genome.

Fifteen new risk loci for coronary artery disease highlight arterial-wall-specific mechanisms

Coronary artery disease (CAD) is a leading cause of morbidity and mortality worldwide. Although 58 genomic regions have been associated with CAD thus far, most of the heritability is unexplained, indicating that additional susceptibility loci await identification. An efficient discovery strategy may be larger-scale evaluation of promising associations suggested by genome-wide association studies (GWAS). Hence, we genotyped 56,309 participants using a targeted gene array derived from earlier GWAS results and performed meta-analysis of results with 194,427 participants previously genotyped, totaling 88,192 CAD cases and 162,544 controls. We identified 25 new SNP–CAD associations (P < 5 × 10−8, in fixed-effects meta-analysis) from 15 genomic regions, including SNPs in or near genes involved in cellular adhesion, leukocyte migration and atherosclerosis (PECAM1, rs1867624), coagulation and inflammation (PROCR, rs867186 (p.Ser219Gly)) and vascular smooth muscle cell differentiation (LMOD1, rs2820315). Correlation of these regions with cell-type-specific gene expression and plasma protein levels sheds light on potential disease mechanisms.

Collaborative pooled analysis of data on C-reactive protein gene variants and coronary disease: judging causality by Mendelian randomisation

Many prospective studies have reported associations between circulating C-reactive protein (CRP) levels and risk of coronary heart disease (CHD), but causality remains uncertain. Studies of CHD are being conducted that involve measurement of common polymorphisms of the CRP gene known to be associated with circulating concentrations, thereby utilising these variants as proxies for circulating CRP levels. By analysing data from several studies examining the association between relevant CRP polymorphisms and CHD risk, the present collaboration will undertake a Mendelian randomisation analysis to help assess the likelihood of any causal relevance of CRP levels to CHD risk. A central database is being established containing individual data on CRP polymorphisms, circulating CRP levels, and major coronary outcomes as well as age, sex and other relevant characteristics. Associations between CRP polymorphisms or haplotypes and CHD will be evaluated under different circumstances. This collaboration comprises, at present, about 37,000 CHD outcomes and about 120,000 controls, which should yield suitably precise findings to help judge causality. This work should advance understanding of the relevance of low-grade inflammation to CHD and indicate whether or not CRP itself is involved in long-term pathogenesis.

Rare Exome Sequence Variants in CLCN6 Reduce Blood Pressure Levels and Hypertension Risk

Background—Rare genetic variants influence blood pressure (BP). Methods and Results—Whole-exome sequencing was performed on DNA samples from 17 956 individuals of European ancestry and African ancestry (14 497, first-stage discovery and 3459, second-stage discovery) to examine the effect of rare variants on hypertension and 4 BP traits: systolic BP, diastolic BP, pulse pressure, and mean arterial pressure. Tests of ≈170 000 common variants (minor allele frequency, ≥1%; statistical significance, P⩽2.9×10−7) and gene-based tests of rare variants (minor allele frequency, <1%; ≈17 000 genes; statistical significance, P⩽1.5×10−6) were evaluated for each trait and ancestry, followed by multiethnic meta-analyses. In the first-stage discovery, rare coding variants (splicing, stop-gain, stop-loss, nonsynonymous variants, or indels) in CLCN6 were associated with lower diastolic BP (cumulative minor allele frequency, 1.3%; &bgr;=−3.20; P=4.1×10−6) and were independent of a nearby common variant (rs17367504) previously associated with BP. CLCN6 rare variants were also associated with lower systolic BP (&bgr;=−4.11; P=2.8×10−4), mean arterial pressure (&bgr;=−3.50; P=8.9×10−6), and reduced hypertension risk (odds ratio, 0.72; P=0.017). Meta-analysis of the 2-stage discovery samples showed that CLCN6 was associated with lower diastolic BP at exome-wide significance (cumulative minor allele frequency, 1.1%; &bgr;=−3.30; P=5.0×10−7). Conclusions—These findings implicate the effect of rare coding variants in CLCN6 in BP variation and offer new insights into BP regulation.

Associations between incident ischemic stroke events and stroke and cardiovascular disease-related genome-wide association studies single nucleotide polymorphisms in the population architecture using genomics and epidemiology study

Background— Genome-wide association studies (GWAS) have identified loci associated with ischemic stroke (IS) and cardiovascular disease (CVD) in European-descent individuals, but their replication in different populations has been largely unexplored. Methods and Results— Nine single nucleotide polymorphisms (SNPs) selected from GWAS and meta-analyses of stroke, and 86 SNPs previously associated with myocardial infarction and CVD risk factors, including blood lipids (high density lipoprotein [HDL], low density lipoprotein [LDL], and triglycerides), type 2 diabetes, and body mass index (BMI), were investigated for associations with incident IS in European Americans (EA) N=26 276, African-Americans (AA) N=8970, and American Indians (AI) N=3570 from the Population Architecture using Genomics and Epidemiology Study. Ancestry-specific fixed effects meta-analysis with inverse variance weighting was used to combine study-specific log hazard ratios from Cox proportional hazards models. Two of 9 stroke SNPs (rs783396 and rs1804689) were associated with increased IS hazard in AA; none were significant in this large EA cohort. Of 73 CVD risk factor SNPs tested in EA, 2 (HDL and triglycerides SNPs) were associated with IS. In AA, SNPs associated with LDL, HDL, and BMI were significantly associated with IS (3 of 86 SNPs tested). Out of 58 SNPs tested in AI, 1 LDL SNP was significantly associated with IS. Conclusions— Our analyses showing lack of replication in spite of reasonable power for many stroke SNPs and differing results by ancestry highlight the need to follow up on GWAS findings and conduct genetic association studies in diverse populations. We found modest IS associations with BMI and lipids SNPs, though these findings require confirmation.

Whole-exome imputation of sequence variants identified two novel alleles associated with adult body height in African Americans

Adult body height is a quantitative trait for which genome-wide association studies (GWAS) have identified numerous loci, primarily in European populations. These loci, comprising common variants, explain <10% of the phenotypic variance in height. We searched for novel associations between height and common (minor allele frequency, MAF ≥5%) or infrequent (0.5% < MAF < 5%) variants across the exome in African Americans. Using a reference panel of 1692 African Americans and 471 Europeans from the National Heart, Lung, and Blood Institutes (NHLBI) Exome Sequencing Project (ESP), we imputed whole-exome sequence data into 13 719 African Americans with existing array-based GWAS data (discovery). Variants achieving a height-association threshold of P < 5E-06 in the imputed dataset were followed up in an independent sample of 1989 African Americans with whole-exome sequence data (replication). We used P < 2.5E-07 (=0.05/196 779 variants) to define statistically significant associations in meta-analyses combining the discovery and replication sets (N = 15 708). We discovered and replicated three independent loci for association: 5p13.3/C5orf22/rs17410035 (MAF = 0.10, β = 0.64 cm, P = 8.3E-08), 13q14.2/SPRYD7/rs114089985 (MAF = 0.03, β = 1.46 cm, P = 4.8E-10) and 17q23.3/GH2/rs2006123 (MAF = 0.30; β = 0.47 cm; P = 4.7E-09). Conditional analyses suggested 5p13.3 (C5orf22/rs17410035) and 13q14.2 (SPRYD7/rs114089985) may harbor novel height alleles independent of previous GWAS-identified variants (r(2) with GWAS loci <0.01); whereas 17q23.3/GH2/rs2006123 was correlated with GWAS-identified variants in European and African populations. Notably, 13q14.2/rs114089985 is infrequent in African Americans (MAF = 3%), extremely rare in European Americans (MAF = 0.03%), and monomorphic in Asian populations, suggesting it may be an African-American-specific height allele. Our findings demonstrate that whole-exome imputation of sequence variants can identify low-frequency variants and discover novel variants in non-European populations.

Tissue Factor Pathway Inhibitor, Activated Protein C Resistance, and Risk of Coronary Heart Disease Due To Combined Estrogen Plus Progestin Therapy

Objective— To examine whether tissue factor pathway inhibitor or acquired activated protein C (APC) resistance influences the increased risk of coronary heart disease (CHD) due to estrogen plus progestin therapy. Approach and Results— Prospective nested case–control study of 205 cases of CHD and 481 matched controls in the Women’s Health Initiative randomized trial of estrogen plus progestin therapy. After multivariable covariate adjustment, both baseline tissue factor pathway activity (P=0.01) and APC resistance (P=0.004) were associated positively with CHD risk. Baseline tissue factor pathway activity and APC resistance singly or jointly did not significantly modify the effect of estrogen plus progestin on CHD risk. Compared with placebo, estrogen plus progestin decreased tissue factor pathway inhibitor activity and increased APC resistance but these changes did not seem to modify or mediate the effect of estrogen plus progestin on CHD risk. Conclusions— Tissue factor pathway inhibitor activity and APC resistance are related to CHD risk in women, but may not explain the increased CHD risk due to estrogen plus progestin therapy. The data from this study do not support the clinical use of measuring these hemostatic factors to help stratify risk before hormone therapy. Clinical Trial Registration— URL: http://www.clinicaltrials.gov. Unique identifier: NCT00000611.

Pharmacogenomics study of thiazide diuretics and QT interval in multi-ethnic populations: the cohorts for heart and aging research in genomic epidemiology

Thiazide diuretics, commonly used antihypertensives, may cause QT interval (QT) prolongation, a risk factor for highly fatal and difficult to predict ventricular arrhythmias. We examined whether common single-nucleotide polymorphisms (SNPs) modified the association between thiazide use and QT or its component parts (QRS interval, JT interval) by performing ancestry-specific, trans-ethnic and cross-phenotype genome-wide analyses of European (66%), African American (15%) and Hispanic (19%) populations (N=78u2009199), leveraging longitudinal data, incorporating corrected standard errors to account for underestimation of interaction estimate variances and evaluating evidence for pathway enrichment. Although no loci achieved genome-wide significance (P<5 × 10−8), we found suggestive evidence (P<5 × 10−6) for SNPs modifying the thiazide-QT association at 22 loci, including ion transport loci (for example, NELL1, KCNQ3). The biologic plausibility of our suggestive results and simulations demonstrating modest power to detect interaction effects at genome-wide significant levels indicate that larger studies and innovative statistical methods are warranted in future efforts evaluating thiazide–SNP interactions.