Framingham Heart Study 100K Project: genome-wide associations for blood pressure and arterial stiffness

Background About one quarter of adults are hypertensive and high blood pressure carries increased risk for heart disease, stroke, kidney disease and death. Increased arterial stiffness is a key factor in the pathogenesis of systolic hypertension and cardiovascular disease. Substantial heritability of blood-pressure (BP) and arterial-stiffness suggests important genetic contributions. Methods In Framingham Heart Study families, we analyzed genome-wide SNP (Affymetrix 100K GeneChip) associations with systolic (SBP) and diastolic (DBP) BP at a single examination in 1971–1975 (n = 1260), at a recent examination in 1998–2001 (n = 1233), and long-term averaged SBP and DBP from 1971–2001 (n = 1327, mean age 52 years, 54% women) and with arterial stiffness measured by arterial tonometry (carotid-femoral and carotid-brachial pulse wave velocity, forward and reflected pressure wave amplitude, and mean arterial pressure; 1998–2001, n = 644). In primary analyses we used generalized estimating equations in models for an additive genetic effect to test associations between SNPs and phenotypes of interest using multivariable-adjusted residuals. A total of 70,987 autosomal SNPs with minor allele frequency ≥ 0.10, genotype call rate ≥ 0.80, and Hardy-Weinberg equilibrium p ≥ 0.001 were analyzed. We also tested for association of 69 SNPs in six renin-angiotensin-aldosterone pathway genes with BP and arterial stiffness phenotypes as part of a candidate gene search. Results In the primary analyses, none of the associations attained genome-wide significance. For the six BP phenotypes, seven SNPs yielded p values < 10-5. The lowest p-values for SBP and DBP respectively were rs10493340 (p = 1.7 × 10-6) and rs1963982 (p = 3.3 × 10-6). For the five tonometry phenotypes, five SNPs had p values < 10-5; lowest p-values were for reflected wave (rs6063312, p = 2.1 × 10-6) and carotid-brachial pulse wave velocity (rs770189, p = 2.5 × 10-6) in MEF2C, a regulator of cardiac morphogenesis. We found only weak association of SNPs in the renin-angiotensin-aldosterone pathway with BP or arterial stiffness. Conclusion These results of genome-wide association testing for blood pressure and arterial stiffness phenotypes in an unselected community-based sample of adults may aid in the identification of the genetic basis of hypertension and arterial disease, help identify high risk individuals, and guide novel therapies for hypertension. Additional studies are needed to replicate any associations identified in these analyses.

arterial disease, help identify high risk individuals, and guide novel therapies for hypertension. Additional studies are needed to replicate any associations identified in these analyses.

Background
Hypertension affects about one quarter of adults in industrialized countries [1] and carries a substantial burden of risk for cardiovascular disease (CVD), kidney disease, and death [2]. Increased arterial stiffness is a key factor in the pathogenesis of hypertension in older people and it contributes to the development of hypertensive target organ damage, CVD, and death [3][4][5]. Substantial heritability of blood pressure [6] and arterial stiffness [7]), as measured by arterial tonometry, points to genetic contributions to these cardiovascular phenotypes.
The search for genetic variants contributing to hypertension and arterial stiffness has focused on complementary approaches: linkage applied to rare Mendelian blood pressure disorders and to large family-based studies to identify positional candidate genes, and the study of biologically plausible candidate genes selected by virtue of their role in blood pressure regulation or vascular properties. A great deal is known about mutations responsible for Mendelian blood pressure disorders [8], but neither these rare variants nor more common variants in these genes account for substantial blood pressure variation in the general population. Similarly, although numerous linkage [9] and candidate gene association studies [10] have been conducted, there is a paucity of evidence that common genetic variation contributes to alterations in blood pressure or arterial stiffness in the general population.
Genome-wide association offers the opportunity to conduct analysis of common genetic variants unconstrained by prior knowledge of biological pathways in relation to phenotypes of interest. This approach succeeded in identifying the association of complement factor H with agerelated macular degeneration [11]. The Framingham Heart Study, which enrolled participants without regard to phenotype status, provides a setting for a genome-wide association study in a community-based sample in which selection bias is inherently low. In addition, because of the familial structure of the study, it also provides an opportunity to use genome-wide SNP data for family based association testing (FBAT) and linkage analyses.
In this report we provide results of a genome-wide association study of blood pressure and arterial stiffness, including results of generalized estimating equation (GEE) association testing, FBAT, and linkage, as well as a summary of associations of these phenotypes with candidate genes in the renin-angiotensin-aldosterone pathways.

Study sample
The Framingham Heart Study began in 1948 when 5209 men and women from Framingham, Mass, who were between 28 and 62 years of age were recruited to participate in an observational study [12]. Subjects underwent a medical history, physician-administered physical examination including blood pressure measurement, laboratory tests, and electrocardiography. Examinations have been repeated every 2 years. In 1971, 5124 offspring and spouses of offspring of original participants were recruited into the Framingham Offspring Cohort [13]. The offspring cohort was reexamined approximately every 4 years, except for an 8 year interval between their initial and second visit. All subjects gave written informed consent before each clinic visit, and the examination protocol was approved by the Institutional Review Board at Boston Medical Center (Boston, Mass).

Blood pressure phenotypes
At each clinic visit, the examining physician measured the systolic and diastolic BP in the left arm using a mercury column sphygmomanometer. BP was measured twice at each visit, with the exception of the first Offspring Cohort clinic visit, when it was measured once in about half the participants. Systolic and diastolic pressures were determined by the first and fifth Korotkoff sounds, respectively, and the two BP measurements were averaged to derive the systolic and diastolic pressures for that examination.
Examination cycles for the two cohorts were overlaid temporally as follows [offspring cohort/original cohort (earliest -latest year)]: examination 1/examination 12 (1971-1975), examination 2/examination 16 (1979-1983), examination 3/examination 18 (1983-1987), examination 4/examination 20 (1986-1991), examination 5/ examination 22 (1990-1995), examination 6/examination 24 (1995-1998) and examination 7/examination 26 (1998)(1999)(2000)(2001). Referring to offspring cycle numbers, the six BP phenotypes analyzed for this investigation were residuals for SBP and DBP at Examination 1, at Examination 7, and average of residuals from available Examinations 1 to 7. BP was imputed for treated observations as previously described [6]. No adjustment was made for untreated observations, which constituted the vast majority of BP values. Systolic and diastolic BP phenotypes were analyzed independently. Residuals were obtained from cohort-and examination-specific regression models accounting for sex, age and BMI; for DBP, age-squared was added. For inclusion in long-term BP analyses, each par-ticipant had to have BP measured on at least three examinations over a period of 12 years or more.

Arterial stiffness phenotypes
Arterial tonometry for assessment of arterial stiffness was conducted on Offspring Cohort participants attending their 7th clinic examination. Five primary tonometry phenotypes were analyzed: carotid-femoral and carotid-brachial pulse wave velocity, forward and reflected pressure wave amplitude, and mean arterial pressure. Tonometry was performed in the supine position after 5 minutes of rest. Arterial tonometry with simultaneous ECG recording was obtained from brachial, radial, femoral and carotid arteries using a commercially available tonometer (SPT-301, Millar Instruments, Houston, TX). Carotid-brachial, carotid-radial and carotid-femoral PWV were calculated as previously described [14]. Mean arterial pressure was calculated from the planimetered brachial arterial tracing after calibration to the brachial blood pressure, which was obtained by an oscillometric device. Forward pressure wave amplitude was defined as the difference between pressure at the waveform foot and pressure at the first systolic inflection point or peak of the carotid pressure waveform; reflected pressure wave amplitude was defined as the difference between the central systolic pressure and the pressure at the forward wave peak. Sex-specific regressions were conducted for each tonometry phenotype with the following covariates: age, age 2 , height, weight, to generate sex-specific residuals.

Genotyping methods
Details of the genotyping methods are available in the Executive Summary [15]. Briefly, 112990 autosomal SNPs on the Affymetrix 100K chip were genotyped in the Boston University School of Medicine Genetics Laboratory on the Framingham Heart Study family plate set. SNPs were excluded for the following reasons: minor allele frequency <10% (n = 38062); call rate <80% (n = 2346); Hardy Weinberg equilibrium p value < 0.001 (n = 1595), leaving 70,987 SNPs available for analysis.

Statistical methods
Standardized multivariable-adjusted blood pressure and tonometry residuals were generated as described above. Table 1 lists the covariates used for each phenotype. As described in the Executive Summary [15], we conducted association testing using family based association testing (FBAT), and generalized estimating equations (GEE) applied to the additive genetic effects model. In secondary analyses that used the GEE general genetic effects model, which is more sensitive to recessive genetic effects, to be more conservative, we limited analyses to two phenotypes: long-term SBP and long-term DBP, and we limited eligible SNPs to those with a minor allele frequency >= 0.20 and Hardy-Weinberg equilibrium p value >= 0.05.
The software package Merlin [16] was used to compute exact identity by descent linkage probabilities for allele sharing, and linkage analysis by variance component method was carried out SOLAR using 11,200 SNPs and STRs. Heritability was estimated using variance-components methods (SOLAR). For BP, 2155 study participants were used for examination 1 SBP and DBP, 1479 for examination 7, and 2009 for long-term average; 770 individuals were used in heritability analysis of arterial stiffness phenotypes.
Candidate gene analyses GEE and FBAT additive genetic effect models were run for SNPs in or near 6 genes in the renin-angiotensin-aldosterone pathways. These genes were selected a priori because of a substantial body of literature implicating them in hypertension and altered vascular properties. All SNPs from 200 Kb proximal to the start and extending to 200 kb of the terminus of each gene were included in analysis providing the minor allele frequency was >= 0.1, the genotype call rate was 0.8, and the Hardy-Weinberg equilibrium p value was >= 0.001.

Results
The six primary BP phenotypes were examination 1 SBP and DBP (n = 1260), examination 7 SBP and DBP (n = 1233), and long-term averaged SBP and DBP (n = 1327). The five primary arterial stiffness phenotypes were carotid-femoral and carotid-brachial pulse wave velocity, forward and reflected pressure wave amplitude, and mean arterial pressure (n = 644). The study sample available for BP phenotypes included up to 1327 individuals (mean age 52 years, 54% women for the long-term SBP and DBP phenotypes). The complete list of blood pressure and arterial stiffness phenotypes analyzed and the covariates used in generating sex-specific standardized residuals for each phenotype are listed in Table 1. Full disclosure of all GEE and FBAT associations for the traits listed in Table 1 can be found at the National Center for Biotechnology Information dbGaP website: http://web.ncbi.nlm.nih.gov/ projects/gap/framingham/cgi-bin/ study.cgi?id=phs000007.
Linkage analyses (Table 2c) yielded a LOD score of 3 for long-term SBP on chromosome 15 at 100 Mb. Several tonometry linkage peaks exceeded a LOD score of 3, including a LOD of 5.0 for reflected wave (chromosome 8 at 19 Mb). Heritability estimates (Table 2d) were high for long-term average DBP (h 2 = 0.55) and SBP (h 2 = 0.57), and intermediate for the other BP phenotypes (h 2 = 0.28-0.45). Among the arterial stiffness phenotypes, heritability was high for the reflected arterial waveform (h 2 = 0.66), low for carotid-brachial PWV (h 2 = 0.02), and intermediate for the other phenotypes (h 2 = 0.22-0.43). These heritability results are consistent with our prior findings [6,7].
Secondary analyses using the GEE general genetic effects model (2 degrees of freedom; more sensitive in detecting recessive effects) are presented in Table 3. The lowest p value for long-term DBP was in CCL20 (rs7591163, p = 2.3 × 10 -7 ) and for SBP was in CDH13 (rs3096277, p = 9.9 × 10 -8 ). Of note, SNPs in CDH13, CCL20, and WDR69 were associated with DBP and SBP. GEE general effects models for the tonometry phenotypes identified association of mean arterial pressure with TGFBR2 (rs3773643, p = 2 × 10 -7 ).
Geometric means of GEE association results (additive genetic effect model) for SBP and DBP considered jointly are summarized in Table 4. The lowest p values were noted for Examination 1 BP values (rs10493340, p = 1.5 × 10 -5 ). Geometric means of association results for the 5 tonometry phenotypes considered concurrently yielded its lowest p value (rs10518082, p = 0.002) for DCK.
SNPs in 6 renin-angiotensin-aldosterone pathway genes were analyzed for association with the BP and tonometry phenotypes (Table 5). A total of 69 SNPs qualified for analysis (minor allele frequency >= 0.1, Hardy Weinberg equilibrium p >= 0.001, call rate >= 0.8). For the primary traits there were few associations from GEE models for an additive genetic effect with p values < 0.05 and none with p < 0.001.

Discussion and conclusion
We provide results of genome-wide association study for 6 blood pressure and 5 arterial stiffness phenotypes in a carefully characterized study sample. Association analyses and linkage reveal a number of intriguing results. For the GEE model of additive genetic effects (Table 2a) there were 7 SNPs with p values < 10 -5 for blood pressure and 5 for arterial stiffness phenotypes. Among the GEE additive effect model results the most likely candidate genes were MEF2C, SYNE1, and TNFSF11, which were associated with arterial stiffness. We have not yet attempted replication of our results. Follow-up genotyping of the top SBP and DBP SNPs reported in our study sample in additional Framingham participants is planned; additional replication attempts will be needed in independent samples to confirm any of the association results we report.
FBAT (Table 2b) identified association of COL8A1 with arterial stiffness (p value 6 × 10 -6 for rs792833). This gene codes for type VIII collagen, which is produced by aortic endothelial cells [17], suggesting a biologically plausible association.
Linkage yielded a LOD score of 3, approaching genomewide significance, for long-term SBP on chromosome 15. A meta-analysis of blood pressure and hypertension linkage studies did not identify this as a region of interest [9]. The lower LOD scores for long-term SBP on chromosome 17 (~67 cM) in this investigation compared with our prior findings [6] appears to be largely due to differences in phenotype definition of long-term SBP with the exclusion of early examination BP values in the original cohort participants and the inclusion of offspring cohort examination 7 blood pressures in this analysis. When linkage analyses were repeated with the inclusion of the early original cohort exams using the prior phenotype definitions, the same linkage peak on chromosome 17 emerged (LOD > 4).
For tonometry phenotypes, we found LOD scores for reflected wave amplitude of 5.0 (chromosome 8 at 19 Mb) and 3.2 (chromosome 4, 169 Mb) near peaks for this phenotype that we previously reported in a largely overlapping study sample [7]. Similarly, we once again identified a linkage peak for carotid-femoral pulse wave velocity (LOD 3.0; chromosome 2 at 74 mb).
Compared with the primary GEE model for additive genetic effects (Table 2a), a different set of SNPs was identified in secondary GEE general effects models (Table 3) for long-term DBP and SBP, including 2 SNPs with p values < 10 -6 . Differences in model results may be due to the greater sensitivity of the general model to detect recessive genotype effects. SNPs in CCL20, CDH13, and LPP were associated with both long-term SBP and DBP. GEE general genetic effects models for arterial stiffness phenotypes yielded the lowest p value (p = 1.99 × 10 -7 ) for rs3773643 in TGFBR2, which has been implicated in aortic aneurysm and Loeys-Dietz syndrome [18,19]. Disruption of the aortic wall would be expected to affect arterial stiffness.  Due to high correlations of SBP and DBP (within examination r = 0.77; long-term r = 0.82), joint analyses of SBP and DBP added little to what was identified in individual phenotype analyses. In contrast, joint analyses of the five tonometry phenotypes, which are less highly correlated, identified LOXL2, SYNE1, and MEF2C as attractive candidates. LOXL2 is a member of the lysyl oxidase family of enzymes that initiate cross-linking of collagens and elastin, and alter arterial elasticity [20]. Collagen and elastin cross-links are critical to tensile strength of the extracellular matrix. Mice null for lysyl oxidase (LOX) die perinatally from aortic aneurysm [21]. MEF2C is involved in cardiac morphogenesis and extracellular matrix remode-ling [22]. SYNE1 is involved in aortic vascular smooth muscle differentiation [23]. To our knowledge, genetic variation in these genes has not previously been shown to be associated with alterations in arterial properties in humans. Whether our results provide nominal evidence of such association or merely chance findings remains to be determined.
Since none of the primary associations attained genomewide significance, this investigation should be viewed as hypothesis generating. Association analyses for SNPs in six renin-angiotensin-aldosterone pathway genes showed weak evidence of association. Negative results for these candidate genes may be due in part to incomplete linkage disequilibrium coverage of these genes by the SNPs in this genome-wide scan. It is likely that the vast majority of low p values from association analyses are due to chance. Replication studies in other populations, using a genomewide approach or selective genotyping is needed to establish if any of our results are indicative of true positive associations.
We provide results of genome-wide association testing for blood pressure and arterial stiffness phenotypes obtained in a carefully described community-based sample of adults who were recruited without regard to disease status. Additional studies are needed to validate these results. Finding genetic variants associated with hypertension or altered arterial properties may aid in the identification of high risk individuals and in the development of new targeted therapies for hypertension. Our report is one of the earlier genome-wide association studies of blood pressure. Several additional studies, some with larger sample size and others with more dense genome-wide coverage of common variation will follow. In that regard, a 550 k SNP genome-wide association study in approximately 9400 Framingham Heart Study participants across three generations is underway and results from that study will help in the interpretation of the findings we report in this manuscript.
results, and drafted the manuscript. GFM conceived of the FHS tonometry project and assisted in securing funding, planned the analyses, and critically revised the manuscript. RSV provided critical input in conceiving the project, securing the funding, planning the analyses and critically revising the manuscript. TJW contributed to design, analysis, and critical review of the manuscript.
We also wish to acknowledge the contributions of Dr. Christopher J. O'Donnell to this project.
This article has been published as part of BMC Medical Genetics Volume 8 Supplement 1, 2007: The Framingham Heart Study 100,000 single nucleotide polymorphisms resource. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2350/8?issue=S1.