Thomas Lumley

Genome-wide association study of blood pressure and hypertension

Blood pressure is a major cardiovascular disease risk factor. To date, few variants associated with interindividual blood pressure variation have been identified and replicated. Here we report results of a genome-wide association study of systolic (SBP) and diastolic (DBP) blood pressure and hypertension in the CHARGE Consortium (n = 29,136), identifying 13 SNPs for SBP, 20 for DBP and 10 for hypertension at P < 4 × 10−7. The top ten loci for SBP and DBP were incorporated into a risk score; mean BP and prevalence of hypertension increased in relation to the number of risk alleles carried. When ten CHARGE SNPs for each trait were included in a joint meta-analysis with the Global BPgen Consortium (n = 34,433), four CHARGE loci attained genome-wide significance (P < 5 × 10−8) for SBP (ATP2B1, CYP17A1, PLEKHA7, SH2B3), six for DBP (ATP2B1, CACNB2, CSK-ULK3, SH2B3, TBX3-TBX5, ULK4) and one for hypertension (ATP2B1). Identifying genes associated with blood pressure advances our understanding of blood pressure regulation and highlights potential drug targets for the prevention or treatment of hypertension.

Genome-wide Analysis of Genetic Loci Associated With Alzheimer Disease

CONTEXT Genome-wide association studies (GWAS) have recently identified CLU, PICALM, and CR1 as novel genes for late-onset Alzheimer disease (AD). OBJECTIVES To identify and strengthen additional loci associated with AD and confirm these in an independent sample and to examine the contribution of recently identified genes to AD risk prediction in a 3-stage analysis of new and previously published GWAS on more than 35,000 persons (8371 AD cases). DESIGN, SETTING, AND PARTICIPANTS In stage 1, we identified strong genetic associations (P < 10(-3)) in a sample of 3006 AD cases and 14,642 controls by combining new data from the population-based Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (1367 AD cases [973 incident]) with previously reported results from the Translational Genomics Research Institute and the Mayo AD GWAS. We identified 2708 single-nucleotide polymorphisms (SNPs) with P < 10(-3). In stage 2, we pooled results for these SNPs with the European AD Initiative (2032 cases and 5328 controls) to identify 38 SNPs (10 loci) with P < 10(-5). In stage 3, we combined data for these 10 loci with data from the Genetic and Environmental Risk in AD consortium (3333 cases and 6995 controls) to identify 4 SNPs with P < 1.7x10(-8). These 4 SNPs were replicated in an independent Spanish sample (1140 AD cases and 1209 controls). Genome-wide association analyses were completed in 2007-2008 and the meta-analyses and replication in 2009. MAIN OUTCOME MEASURE Presence of Alzheimer disease. RESULTS Two loci were identified to have genome-wide significance for the first time: rs744373 near BIN1 (odds ratio [OR],1.13; 95% confidence interval [CI],1.06-1.21 per copy of the minor allele; P = 1.59x10(-11)) and rs597668 near EXOC3L2/BLOC1S3/MARK4 (OR, 1.18; 95% CI, 1.07-1.29; P = 6.45x10(-9)). Associations of these 2 loci plus the previously identified loci CLU and PICALM with AD were confirmed in the Spanish sample (P < .05). However, although CLU and PICALM were confirmed to be associated with AD in this independent sample, they did not improve the ability of a model that included age, sex, and APOE to predict incident AD (improvement in area under the receiver operating characteristic curve from 0.847 to 0.849 in the Rotterdam Study and 0.702 to 0.705 in the Cardiovascular Health Study). CONCLUSIONS Two genetic loci for AD were found for the first time to reach genome-wide statistical significance. These findings were replicated in an independent population. Two recently reported associations were also confirmed. These loci did not improve AD risk prediction. While not clinically useful, they may implicate biological pathways useful for future research.

Systematic identification of trans eQTLs as putative drivers of known disease associations

Identifying the downstream effects of disease-associated SNPs is challenging. To help overcome this problem, we performed expression quantitative trait locus (eQTL) meta-analysis in non-transformed peripheral blood samples from 5,311 individuals with replication in 2,775 individuals. We identified and replicated trans eQTLs for 233 SNPs (reflecting 103 independent loci) that were previously associated with complex traits at genome-wide significance. Some of these SNPs affect multiple genes in trans that are known to be altered in individuals with disease: rs4917014, previously associated with systemic lupus erythematosus (SLE), altered gene expression of C1QB and five type I interferon response genes, both hallmarks of SLE. DeepSAGE RNA sequencing showed that rs4917014 strongly alters the 3′ UTR levels of IKZF1 in cis, and chromatin immunoprecipitation and sequencing analysis of the trans-regulated genes implicated IKZF1 as the causal gene. Variants associated with cholesterol metabolism and type 1 diabetes showed similar phenomena, indicating that large-scale eQTL mapping provides insight into the downstream effects of many trait-associated variants.

Genetic loci for retinal arteriolar microcirculation

Narrow arterioles in the retina have been shown to predict hypertension as well as other vascular diseases, likely through an increase in the peripheral resistance of the microcirculatory flow. In this study, we performed a genome-wide association study in 18,722 unrelated individuals of European ancestry from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium and the Blue Mountain Eye Study, to identify genetic determinants associated with variations in retinal arteriolar caliber. Retinal vascular calibers were measured on digitized retinal photographs using a standardized protocol. One variant (rs2194025 on chromosome 5q14 near the myocyte enhancer factor 2C MEF2C gene) was associated with retinal arteriolar caliber in the meta-analysis of the discovery cohorts at genome-wide significance of P-value <5×10−8. This variant was replicated in an additional 3,939 individuals of European ancestry from the Australian Twins Study and Multi-Ethnic Study of Atherosclerosis (rs2194025, P-value = 2.11×10−12 in combined meta-analysis of discovery and replication cohorts). In independent studies of modest sample sizes, no significant association was found between this variant and clinical outcomes including coronary artery disease, stroke, myocardial infarction or hypertension. In conclusion, we found one novel loci which underlie genetic variation in microvasculature which may be relevant to vascular disease. The relevance of these findings to clinical outcomes remains to be determined.

Genomewide Association Studies of Stroke

BACKGROUND The genes underlying the risk of stroke in the general population remain undetermined. METHODS We carried out an analysis of genomewide association data generated from four large cohorts composing the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium, including 19,602 white persons (mean [+/-SD] age, 63+/-8 years) in whom 1544 incident strokes (1164 ischemic strokes) developed over an average follow-up of 11 years. We tested the markers most strongly associated with stroke in a replication cohort of 2430 black persons with 215 incident strokes (191 ischemic strokes), another cohort of 574 black persons with 85 incident strokes (68 ischemic strokes), and 652 Dutch persons with ischemic stroke and 3613 unaffected persons. RESULTS Two intergenic single-nucleotide polymorphisms on chromosome 12p13 and within 11 kb of the gene NINJ2 were associated with stroke (P<5x10(-8)). NINJ2 encodes an adhesion molecule expressed in glia and shows increased expression after nerve injury. Direct genotyping showed that rs12425791 was associated with an increased risk of total (i.e., all types) and ischemic stroke, with hazard ratios of 1.30 (95% confidence interval [CI], 1.19 to 1.42) and 1.33 (95% CI, 1.21 to 1.47), respectively, yielding population attributable risks of 11% and 12% in the discovery cohorts. Corresponding hazard ratios were 1.35 (95% CI, 1.01 to 1.79; P=0.04) and 1.42 (95% CI, 1.06 to 1.91; P=0.02) in the large cohort of black persons and 1.17 (95% CI, 1.01 to 1.37; P=0.03) and 1.19 (95% CI, 1.01 to 1.41; P=0.04) in the Dutch sample; the results of an underpowered analysis of the smaller black cohort were nonsignificant. CONCLUSIONS A genetic locus on chromosome 12p13 is associated with an increased risk of stroke.

Meta-analyses of genome-wide association studies identify multiple loci associated with pulmonary function

Spirometric measures of lung function are heritable traits that reflect respiratory health and predict morbidity and mortality. We meta-analyzed genome-wide association studies for two clinically important lung-function measures: forced expiratory volume in the first second (FEV1) and its ratio to forced vital capacity (FEV1/FVC), an indicator of airflow obstruction. This meta-analysis included 20,890 participants of European ancestry from four CHARGE Consortium studies: Atherosclerosis Risk in Communities, Cardiovascular Health Study, Framingham Heart Study and Rotterdam Study. We identified eight loci associated with FEV1/FVC (HHIP, GPR126, ADAM19, AGER-PPT2, FAM13A, PTCH1, PID1 and HTR4) and one locus associated with FEV1 (INTS12-GSTCD-NPNT) at or near genome-wide significance (P < 5 × 10−8) in the CHARGE Consortium dataset. Our findings may offer insights into pulmonary function and pathogenesis of chronic lung disease.

Multiple loci associated with indices of renal function and chronic kidney disease

Chronic kidney disease (CKD) has a heritable component and is an important global public health problem because of its high prevalence and morbidity. We conducted genome-wide association studies (GWAS) to identify susceptibility loci for glomerular filtration rate, estimated by serum creatinine (eGFRcrea) and cystatin C (eGFRcys), and CKD (eGFRcrea < 60 ml/min/1.73 m2) in European-ancestry participants of four population-based cohorts (ARIC, CHS, FHS, RS; n = 19,877; 2,388 CKD cases), and tested for replication in 21,466 participants (1,932 CKD cases). We identified significant SNP associations (P < 5 × 10−8) with CKD at the UMOD locus, with eGFRcrea at UMOD, SHROOM3 and GATM-SPATA5L1, and with eGFRcys at CST and STC1. UMOD encodes the most common protein in human urine, Tamm-Horsfall protein, and rare mutations in UMOD cause mendelian forms of kidney disease. Our findings provide new insights into CKD pathogenesis and underscore the importance of common genetic variants influencing renal function and disease.

Common variants at ten loci influence QT interval duration in the QTGEN Study.

QT interval duration reflecting myocardial repolarization on the electrocardiogram is a heritable risk factor for sudden cardiac death and drug-induced arrhythmias. We conducted a meta-analysis of 3 genome-wide association studies in 13,685 individuals of European ancestry from the Framingham Heart Study, the Rotterdam Study and the Cardiovascular Health Study. We observed associations at P < 5×10−8 with variants in NOS1AP, KCNQ1, KCNE1, KCNH2 and SCN5A, known to be involved in myocardial repolarization and Mendelian Long QT Syndromes. Associations at five novel loci included 16q21 near NDRG4 and GINS3, 6q22 near PLN, 1p36 near RNF207, 16p13 near LITAF and 17q12 near LIG3 and RIFFL. Collectively, the 14 independent variants at these 10 loci explain 5.4–6.5% of variation in QT interval. Identifying the causal variants and defining their impact on myocardial repolarization may add incrementally to the prevention of sudden cardiac death and drug-induced arrhythmias.QT interval duration, reflecting myocardial repolarization on the electrocardiogram, is a heritable risk factor for sudden cardiac death and drug-induced arrhythmias. We conducted a meta-analysis of three genome-wide association studies in 13,685 individuals of European ancestry from the Framingham Heart Study, the Rotterdam Study and the Cardiovascular Health Study, as part of the QTGEN consortium. We observed associations at P < 5 × 10−8 with variants in NOS1AP, KCNQ1, KCNE1, KCNH2 and SCN5A, known to be involved in myocardial repolarization and mendelian long-QT syndromes. Associations were found at five newly identified loci, including 16q21 near NDRG4 and GINS3, 6q22 near PLN, 1p36 near RNF207, 16p13 near LITAF and 17q12 near LIG3 and RFFL. Collectively, the 14 independent variants at these 10 loci explain 5.4–6.5% of the variation in QT interval. These results, together with an accompanying paper, offer insights into myocardial repolarization and suggest candidate genes that could predispose to sudden cardiac death and drug-induced arrhythmias.

Case-crossover analyses of air pollution exposure data: referent selection strategies and their implications for bias.

The case–crossover design has been widely used to study the association between short-term air pollution exposure and the risk of an acute adverse health event. The design uses cases only; for each individual case, exposure just before the event is compared with exposure at other control (or “referent”) times. Time-invariant confounders are controlled by making within-subject comparisons. Even more important in the air pollution setting is that time-varying confounders can also be controlled by design by matching referents to the index time. The referent selection strategy is important for reasons in addition to control of confounding. The case–crossover design makes the implicit assumption that there is no trend in exposure across the referent times. In addition, the statistical method that is used—conditional logistic regression—is unbiased only with certain referent strategies. We review here the case–crossover literature in the air pollution context, focusing on key issues regarding referent selection. We conclude with a set of recommendations for choosing a referent strategy with air pollution exposure data. Specifically, we advocate the time-stratified approach to referent selection because it ensures unbiased conditional logistic regression estimates, avoids bias resulting from time trend in the exposure series, and can be tailored to match on specific time-varying confounders.

Referent selection in case-crossover analyses of acute health effects of air pollution.

The case-crossover design was proposed for the study of a transient effect of an intermittent exposure on the subsequent occurrence of a rare acute-onset disease. This design can be an alternative to Poisson time series regression for studying the health effects of fine particulate matter air pollution. Characteristics of time-series of particulate matter, including long-term time trends, seasonal trends, and short-term autocorrelations, require that referent selection in the case-crossover design be considered carefully and adapted to minimize bias. We performed simulations to evaluate the bias associated with various referent selection strategies for a proposed case-crossover study of associations between particulate matter and primary cardiac arrest. Some a priori reasonable strategies were associated with a relative bias as large as 10%, but for most strategies the relative bias was less than 2% with confidence interval coverage within 1% of the nominal level. We show that referent selection for case–crossover designs raises the same issues as selection of smoothing method for time series analyses. In addition, conditional logistic regression analysis is not strictly valid for some case–crossover designs, introducing further bias.