Hypertension and Genetic Variation in Endothelial-Specific Genes

Abstract

Genome-wide association (GWA) studies usually detect common genetic variants with low-to-medium effect sizes. Many contributing variants are not revealed, since they fail to reach significance after strong correction for multiple comparisons. The WTCCC study for hypertension, for example, failed to identify genome-wide significant associations. We hypothesized that genetic variation in genes expressed specifically in the endothelium may be important for hypertension development. Results from the WTCCC study were combined with previously published gene expression data from mice to specifically investigate SNPs located within endothelial-specific genes, bypassing the requirement for genome-wide significance. Six SNPs from the WTCCC study were selected for independent replication in 5205 hypertensive patients and 5320 population-based controls, and successively in a cohort of 16537 individuals. A common variant (rs10860812) in the DRAM (damage-regulated autophagy modulator) locus showed association with hypertension (P = 0.008) in the replication study. The minor allele (A) had a protective effect (OR = 0.93; 95% CI 0.88–0.98 per A-allele), which replicates the association in the WTCCC GWA study. However, a second follow-up, in the larger cohort, failed to reveal an association with blood pressure. We further tested the endothelial-specific genes for co-localization with a panel of newly discovered SNPs from large meta-GWAS on hypertension or blood pressure. There was no significant overlap between those genes and hypertension or blood pressure loci. The result does not support the hypothesis that genetic variation in genes expressed in endothelium plays an important role for hypertension development. Moreover, the discordant association of rs10860812 with blood pressure in the case control study versus the larger Malmö Preventive Project–study highlights the importance of rigorous replication in multiple large independent studies.

About the Authors:

Erik Larsson

* E-mail: [email protected] (EL); [email protected] (PL)

Affiliation: Department of Medical Biochemistry and Cell Biology, Institute of Biomedicine, The Sahlgrenska Academy at University of Gothenburg, Gothenburg, Sweden

Björn Wahlstrand

Affiliation: Institute of Medicine, The Sahlgrenska Academy at University of Gothenburg, Gothenburg, Sweden

Bo Hedblad

Affiliation: Clinical research Center (CRC), Malmö University Hospital, Malmö, Sweden

Thomas Hedner

Affiliation: Institute of Medicine, The Sahlgrenska Academy at University of Gothenburg, Gothenburg, Sweden

Sverre E. Kjeldsen

Affiliation: Department of Cardiology/Cardiovascular and Renal Research Center, Ullevaal University Hospital, Oslo, Norway

Olle Melander

Contributed equally to this work with: Olle Melander, Per Lindahl

Affiliation: Clinical research Center (CRC), Malmö University Hospital, Malmö, Sweden

Per Lindahl

Contributed equally to this work with: Olle Melander, Per Lindahl

* E-mail: [email protected] (EL); [email protected] (PL)

Affiliation: Wallenberg Laboratory, Institute of Medicine, The Sahlgrenska Academy at University of Gothenburg, Gothenburg, Sweden

Introduction

Hypertension is the major global risk factor for coronary heart disease and stroke. The pathogenesis is poorly understood and the primary cause is unknown in 90–95% of cases. Heritability has been estimated to between 30 and 50% and ambitious efforts have been made to elucidate the genetic basis. While genome-wide association (GWA) studies have been effective at identifying novel genetic risk loci for a wide range of diseases, hypertension has proved more challenging. The Welcome Trust Case Control Consortium (WTCCC) study [1], encompassing 2000 subjects for each of 7 diseases and 3000 common controls, was successful overall but failed in the case of hypertension. Robust associations of genetic variants with hypertension were eventually identified based on large meta-analyses and multi-center studies of 30,000–200,000 individuals [2]–[6].

One possible reason for this is that hypertension may have relatively few common risk alleles of large effect sizes [7]. Thus, discovery of novel variants for hypertension on the whole genome scale is challenging, and many true variants are likely fail to reach significance after correction for multiple testing. Although no single nucleotide polymorphisms (SNPs) reached the P<5×10−7 threshold required for significance in the WTCCC hypertension study, the number and distribution of signals in the 10−7 to 10−4 range was similar to that observed for other diseases [1]. It is therefore possible that several truly hypertension-associated variations are hidden among the long list of moderately significant SNPs.

Although results have been contradictive, previous association studies indicate that genetic variation in endothelial genes such as endothelin-1 and endothelial nitric oxide synthase, both related to blood pressure (BP) regulation, may contribute to essential hypertension [8]–[11]. Based on the hypothesis that genetic variation in endothelium-specific genes may influence BP, we combined data from the WTCCC GWA study with previously published microarray gene expression data from mice [12], to select a subset of six SNPs, all located within endothelial marker genes, for independent replication in a case control study comprising more than 10,000 individuals. One SNP was further evaluated in a population-based study with 16537 participants (Figure 1).

[Figure omitted. See PDF.]

Figure 1. Study design and workflow.

https://doi.org/10.1371/journal.pone.0062035.g001

Methods

Ethics Statement

The Nordic Diltiazem study (NORDIL) and the Malmö Diet and Cancer Cardiovascular Cohort study (MDC-CC) were approved by the human ethics committees at the University of Gothenburg and Lund University, respectively. The Malmö Preventive Project (MPP) study was approved by the local ethics committee of Southern Sweden. All study participants had given written consent.

The Nordic Diltiazem Study

The Nordic Diltiazem study (NORDIL) is an intervention trial that between 1992 and 1999 included 10881 Swedish and Norwegian patients diagnosed with severe (grade 2) hypertension based on repeated diastolic BP≥100 mmHg on different occasions, and prospectively compared cardiovascular outcome in patients randomized to diltiazem-based antihypertensive treatment as compared to patients randomized to diuretic and/or β-blocker based antihypertensive treatment. The design and main results have previously been described in detail [13]. The primary endpoint was fatal and non-fatal stroke, fatal and non-fatal myocardial infarction and other cardiovascular deaths, and there was no significant difference in cardiovascular outcome between the two treatment regimens. The Swedish subcohort of NORDIL participated in a genetic study and we obtained whole blood samples for DNA isolation from 5262 patients, of whom 5205 were successfully genotyped for the main SNP rs10860812 for the present investigation. BP was measured in the supine position after 10 minutes rest. These patients with diagnosed grade 2 hypertension formed our case group.

The Malmö Diet and Cancer Cardiovascular Cohort study

As a population control we used a Swedish cohort study, the population based Malmö Diet and Cancer Cardiovascular Cohort (MDC-CC) [14], which was designed to investigate the epidemiology of carotid artery disease. From a community-based prospective epidemiologic cohort of 28,449 persons enrolled between 1991 and 1996, 6,103 persons were randomly selected to participate in the MDC-CC. Whole blood samples for DNA extraction was obtained from 5445 subjects and genotypes for the main SNP rs10860812 were obtained from 5320 of these subjects. In the MDC-CC, BP was measured using a mercury-column sphygmomanometer after 10 minutes of rest in the supine position. Cardiovascular events (fatal and non-fatal stroke, fatal and non-fatal myocardial infarction and other cardiovascular deaths) were recorded during follow-up using national and local registers. Follow-up extended until December 31st 2005. Intima-media thickness (IMT) of the carotid artery was measured using 2D B-mode ultrasound as described previously [15]. The mean IMT of a 10 mm section of the common carotid artery (IMTmean) and the maximum IMT of the carotid bulb (IMTmax) was recorded.

Characteristics of the grade 2 hypertension patients (NORDIL) and the population based control sample (MDC-CC) are shown in Table 1.

[Figure omitted. See PDF.]

Table 1. Characteristics of study subjects.

https://doi.org/10.1371/journal.pone.0062035.t001

The Malmö Preventive Project

33,346 individuals from the city of Malmö (22,400 men and 10902 women) with a mean age of 49 years participated in a health study, the Malmö Preventive Project (MPP), during 1974–1992 (attendance rate 71%) [16]. The participants underwent a physical examination at the onset of the study with measurement of BP. Of those 33,346 participants, 5486 individuals were lost from follow-up (4931 died and 551 were lost from other reasons). 18,240 individuals participated in a re-screening visit during 2002–2006 (73% of the invited individuals) for renewed physical examination and BP determination. Of the 18,240 attending individuals, 1703 were excluded from the present investigation since clinical information or DNA was lacking. Thus, 16537 individuals were included for sequential replication. BP was monitored in different ways at the two occasions. At the onset of the study, a first BP reading was taken after 1 minute of rest in the supine position. A second BP reading was taken in the upright position after 1 minute of standing. The procedure was repeated after a 10 minute rest in the supine position. The average BP of subjects with at least three valid measurements was included in the present study. At the re-screening, BP was determined twice in the supine position and the average BP of individuals with two valid measurements was included in the present study.

Characteristics of the MPP study subjects are shown in Table 1.

Genotyping

DNA was extracted from granulocyte or buffy coat suspensions, maintained at –80°C from the time of enrolment. Samples were thawed rapidly at 37°C and 200 µL aliquots were subjected to QiaAmp mini-preps in 96-well format (Qiagen) according to the manufacturer’s instructions. SNPs rs893881, rs10060812, rs2269772, rs4684243, rs4981504 and rs6891143 were genotyped using 2.5 ng of DNA on the 7900HT instrument using TaqMan SNP Genotyping Assays (Applied Biosystems) in a total reaction volume of 6 µL in 384-well microtiter plates, according to the manufacturer’s instructions.

Animals

Adult C57Bl/6 mice were kept in groups at the Laboratory for Experimental Biomedicine at University of Gothenburg in a 12 h day/12 h night light cycle with food and water administered ad libitum in a temperature- and humidity-controlled room. All animal experiments were approved by the animal research ethics committee in Gothenburg, Sweden.

Isolation of Microvascular Fragments

Microvascular fragments were isolated from adult C57/Bl6 mouse brain and kidney. Brains or kidneys were dissected out, minced into pieces and digested with 5 mg Collagenase A (Roche Diagnostics GmbH, Mannheim, Germany) dissolved in Hanks’ balanced salt solution (HBSS, Invitrogen AB, Lidingö, Sweden) including 1% BSA and 100U DNase at 37°C for 15 min with gentle agitation. The tissue was then gently pressed through a 100 µm cell strainer (Falcon, BD Biosciences, Stockholm, Sweden). Cells were washed out from the strainer in 2 ml of HBSS/1%BSA/100U DNase, pelleted at 200 g for 5 min, suspended in 1.5 ml HBSS/1% BSA/100U DNase, and again pelleted and resuspended. Rat anti-PECAM (BD PharMingen, San Diego, CA, USA) antibody (Ab) -coated magnetic beads (Dynabeads M-450, sheep anti-Rat IgG, Dynal A.S., Oslo, Norway) were added, and after incubation at 4°C for 30 min with gentle agitation, microvascular fragments were isolated with a magnetic particle concentrator (MPC, Dynal) and washed three times with HBSS/1% BSA.

qPCR

mRNA expression levels were determined using SYBR Green quantitative qPCR (95°, 55°, 72°, 40 cycles) on a 7900HT instrument (Applied Biosystems) using the following primers: Adcy4, 5′-CTTTGGGTGGCTTCTCTCTG-3′ and 5′-ATGGCGTACACGGTGAAGAT-3′; Gpr116, 5′-AAGACATGAGATCGCCAAGG-3′ and 5′-TTGGCCTCAGTAGCTCTTCC-3′; Fgd5, 5′-GCTCAGGAGCTGCTGTCTTC-3′ and 5′-CCCTGGTGAAGGTCACAGATA-3′; Arap3, 5′-GACTGAGCCCAATCTTCTGG-3′ and 5′-TCGCCTGAGAACTATCTGGA-3′; Itga3, 5′-TTGAGGATATGTGGCTTGGA-3′ and 5′-ATGCCGGTCTGCAAGTAGTC-3′; Dram, 5′-GACACAGGAACAACTCCTCCA-3′ and 5′-AACGGGAGTGCTGAAGTAGC-3′; Nebl, 5′-ATGTTTCCACTGCGAGGTCT-3′ and 5′-TGCAGTTCACTTTGCTGCTT-3′. Expression of Gapdh and Tie2 were measured using TaqMan assays under standard cycling conditions (Mm99999915_g1 and Mm00443242_m1, Applied Biosystems). Relative expression levels were determined using the standard curve method [17].

Bioinformatical Identification of Candidate Hypertension SNPs

Tab-delimited text files with hypertension association statistics for 469,557 SNPs were obtained from the WTCCC [2]. A list of 71 genes predicted to be specifically expressed in the microvascular endothelium were obtained from a recently published study [12]. Reciprocal human orthologs and genomic coordinates for these genes were determined using ENSEMBL and the UCSC browser [18], [19]. For each gene region, here defined as the start of the first exon until the end of the last exon, the SNP with the lowest trend P-value in the WTCCC study was identified (Table 2).

[Figure omitted. See PDF.]

Table 2. SNPs in the WTCCC hypertension study mapped onto endothelial-specific genes.

https://doi.org/10.1371/journal.pone.0062035.t002

Co-localization of SNPs Associated with Hypertension or Blood Pressure in Meta-GWAS with EC-specific Genes

549 SNPs with P<10−5 based on data from 18 GWAS or meta-GWAS on hypertension or blood pressure were extracted from the PheGenI database (http://www.ncbi.nlm.nih.gov/gap/PheGenI) [2], [17]–[33]. 670 putative associated genes were within 100 kb these SNPs, and these were evaluated for overrepresentation of EC-specific genes (71 gene list) using Fisher’s exact test.

Statistical Analysis

We assumed an additive model of inheritance and calculated allelic odds ratios (OR) and 95% confidence intervals (95% CI) for each SNP in relation to the main dependent variable (belonging to the grade 2 hypertension case group or belonging to the population controls) and to the secondary outcome variable (presence of cardiovascular events during follow-up in NORDIL or MDC-CC or no cardiovascular event during follow-up in NORDIL or MDC-CC) using crude and multivariate adjusted logistic regression. Continuous variables were related to rs10860812 using linear regression assuming an additive model of inheritance. Deviation from Hardy-Weinberg equilibrium was evaluated with a χ2 test using Levene’s method (1949).

Results

Putative Hypertension-associated SNPs in Endothelial-specific Genes

Since the endothelium has a functional role in BP regulation, we tested the hypothesis that genetic variation in genes expressed specifically in the endothelium may influence hypertension development. SNPs from the published WTCCC GWA analysis for human hypertension [2] were mapped onto a list of 71 endothelial genes obtained from a published study, where putative endothelial-specific genes were identified through analysis of a large mouse microarray compendium [12]. For each gene region, the SNP with the strongest association to hypertension was identified (additive model). In consistency with our hypothesis, the number of moderately significant SNPs (P<0.005) was significantly larger among the set of EC-specific genes compared to remaining genes (10% v.s. 4%, P = 0.03, Fisher’s exact test). Genes were ranked according to best P-value and SNPs in the top six candidates were selected for replication in an independent material. Association P-values in the WTCCC study for these SNPs were in the range of 1.8×10−3 to 4.1×10−4 (Table 2, italic). In addition, a SNP in the NEBL gene was included for validation due to its relatively strong association in the WTCCC (rs893881, Padd = 1.5×10−5). NEBL is selectively expressed in endothelial cells and was highly ranked in the original analysis of EC-specific gene expression [12], but did not qualify for the 71-gene list due to additional strong expression in the heart.

To confirm expression of these genes in the vasculature, microvascular fragments were isolated from mouse brain and kidney using anti-CD31 (PECAM)-coated magnetic beads. Expression in vascular fragments vs. surrounding tissue was subsequently determined using real-time quantitative PCR. All genes were found to be significantly enriched in CD31+ fractions from both tissues (P<0.005) and all except one had >70-fold enrichment in at least one of the tissues (Figure 2). Gapdh was not differentially expressed, while Tie2, included as a positive control, was strongly enriched in CD31+ fractions from both tissues.

[Figure omitted. See PDF.]

Figure 2. Differential mRNA expression in CD31+ microvascular fragments.

Microvascular fragments were isolated from mouse tissues using anti-CD31 (PECAM)-coated magnetic beads. Reverse transcription was performed using equal amounts of RNA and expression levels were assessed using real-time quantitative PCR. The graph shows the ratio between expression in microvascular fragments and the remaining tissue fraction for all genes included in the study. Gapdh and an established endothelial marker, Tie2, were included for comparison. Error bars represent SEM (n = 4, P<0.005 in all cases using Student’s t-test).

https://doi.org/10.1371/journal.pone.0062035.g002

The DRAM Locus and Hypertension

An independent evaluation of the above selected SNPs was performed in 5205 hypertensive patients and 5320 population based controls (NORDIL vs. MDC-CC). Genotyping was successful (average genotyping success rate 98%) for six of the seven SNPs selected for validation, while one SNP (rs2021916, GPR116 locus) did not pass quality control. For all tested SNPs, allele frequencies were fairly similar in the WTCCC study and the replication studies, and Hardy-Weinberg expectations were satisfied in the population control group (P>0.15). The power to show an association to hypertension in the case-control material (NORDIL+MCD-CC), assuming additive model odds ratios of 1.05, 1.10 and 1.15 was 21%, 60% and 90%, respectively.

One common variant (rs10860812) in the DRAM locus showed significant association with hypertension (Table 3), whereas none of the other SNPs were significantly associated. After applying a Bonferroni correction of 6 to account for all tested loci, this variant was still significant at the P<0.05 level. The minor allele (A) showed a protective effect (OR = 0.93; 95% CI 0.88–0.98 per A-allele; P = 0.008), the directionality of which is in concordance with the association that was observed in the WTCCC GWA study. The protective effect remained significant (OR = 0.93; 95% CI 0.87–0.99 per A-allele; P = 0.01) after adjustment for age, sex, body mass index, diabetes, smoking, plasma cholesterol, and previous cardiovascular disease. Exclusion of subjects on antihypertensive therapy in the control group (n = 888) had a small positive effect on the strength of the association (OR = 0.93, 95% CI 0.87–0.98 per A-allele; P = 0.007). We performed interaction tests between the DRAM rs10860812 and each of the 5 other SNPs with hypertension as the outcome. There was no significant interaction between the rs10860812 and any of the 5 other SNPs (P = 0.159–0.887).

[Figure omitted. See PDF.]

Table 3. Association between hypertension and SNPs in the study in WTCCC and NORDIL/MDC-CC.

https://doi.org/10.1371/journal.pone.0062035.t003

Secondary Analyses of DRAM Locus in the MDC-CC Study

Secondary analyses of rs10860812 (additive models) in relation to the continuous BP variable in the MDC-CC population sample did not reach significance, although results pointed toward a protective effect for the A-allele (β-coefficient ± SEM) (−0.425±0.374 mmHg per A-allele, P = 0.26 for systolic BP and −0.200±0.187 mmHg per A-allele, P = 0.29 for diastolic BP). The risk of incident cardiovascular disease during follow-up (n = 894 cardiovascular events in MDC-CC+NORDIL) was reduced in carriers of the rs10860812 A-allele (OR = 0.90, 95% CI 0.82–1.0 per A-allele; P = 0.04). Finally, in patients from the MDC-CC study, where IMT was measured, the rs10860812 A-allele tended to be associated with lower (β-coefficient ± SEM) IMTmax (−0.028±0.015 mm per A-allele; P = 0.06) and IMTmean (−0.006±0.003 mm per A-allele; P = 0.07).

Secondary Analyses of DRAM Locus in the MPP Study

The rs10860812 variant was further evaluated (additive models) in relation to the continuous BP variable in the MPP study comprising 16537 individuals. There was no association between the rs10860812 variant and systolic or diastolic blood pressure neither at the MPP study baseline exam (β-coefficient ± SEM) (0.064±0.158 mmHg per A-allele, P = 0.688 for systolic BP and −0.068±0.097 mmHg per A-allele, P = 0.486 for diastolic BP) or at re-screening 15–25 years later (0.032±0.219 mmHg per A-allele, P = 0.883 for systolic BP and −0.121±0.116 mmHg per A-allele, P = 0.297 for diastolic BP). One possible confounding factor is that pharmacological treatment of hypertensive subjects in the study lowers the observed BP thus reducing the statistical power [20]. However, similar results were obtained when diastolic (+10 mmHg) and systolic BP (+15 mmHg) were corrected in subjects on antihypertensive treatment (data not shown), as described in [6]. Moreover, the risk of incident cardiovascular disease during follow-up was not different in carriers of the rs10860812 A-allele, compared to the control population (data not shown).

Co-localization of EC-specific Genes and SNPs from Meta-GWAS

The lack of significant associations in our study suggests that variants in EC-specific genes may not be associated with hypertension or blood pressure more often than random genes. However, the publication of large meta-GWAS has expanded the repertoire of SNPs associated with those phenotypes, raising the possibility that newly discovered SNPs co-localize with EC-specific genes. To assess this possibility, we tested the 71 EC-specific genes against a recent catalogue of SNPs associated with hypertension or blood pressure [21], comprising a total of 18 GWAS or meta-GWAS [1], [3], [5], [6], [22]–[35]. A set of 529 SNPs with association P<10−5 were downloaded from the genotype phenotype integrator (PheGenI) at NCBI (http://www.ncbi.nlm.nih.gov/gap/PheGenI), and mapped against the genomic locations of the 71 EC-specific genes. EC-specific genes were not significantly over-represented among 670 genes located within 100 kb of those SNPs (5 EC-specific genes located within 100 kb of associated SNPs, P = 0.10, Fisher’s exact test).

Discussion

Rare high-penetrance or common weakly penetrant genetic variants are not readily detected in GWA studies, opening up for alternative approaches that maintain the genome-wide perspective but reduce multiple testing. Here, we assess if genetic variants in endothelial-specific genes are associated with hypertension. The SNPs were selected for validation based on 1) endothelium-specific expression of the gene locus and 2) an association to hypertension that was observed in the WTCCC GWA study, although at a P-value below genome-wide significance. The selection process can thus be described as an intermediate between a candidate gene approach and a genome-wide approach.

Our case control study was substantially larger and our cases had a strictly defined and more severe clinical hypertension (grade 2) compared to the WTCCC, suggesting a greater power as the genetic contribution is likely to be larger in this case. Similar to the WTCCC, our control group was population based and thus included patients with hypertension. Our estimate of the genetic effect on risk therefore has relevance at the population level, but this may also have reduced the power in our study. To further validate effects on BP, SNPs were re-evaluated in relation to the continuous BP variable in the population based MPP study.

Our first replication suggested an association between a SNP (rs10860812) in the DRAM locus on chromosome 12q23 and hypertension. DRAM is a regulator of autophagy that plays a critical role in apoptosis [36]. Although open to question, the involvement of apoptosis in hypertension-related vascular remodeling has been suggested [37]. The minor allele (A) of this SNP was concluded to have a protective effect and this is consistent with results from the WTCCC study. Secondary analyses revealed an association with incident cardiovascular events and borderline significant associations with IMTmax and IMTmean. However, evaluation of rs10860812 in relation to the continuous blood pressure variable in the MDC-CC population study and sequentially in the larger MPP study did not show association with BP. Moreover, there was no association of rs10860812 with cardiovascular disease in the MPP study.

Recent meta-analyses have to some part overcome the problems of the early hypertension GWAS, and 29 loci are now robustly associated with blood pressure and hypertension [3]. The EC-specific gene NRP3 is located near one of the 29 validated loci, rs1173771, and encodes a receptor for natriuretic peptides that are implicated in the maintenance of blood pressure [38], supporting a causal role for the endothelium. A large proportion of blood pressure heritability remains unknown [3], and alternative methods such as the one we describe could help identify additional loci. However, our replication studies of moderately significant SNPs from the WTCCC study did not support a connection between genetic variation in the endothelium and hypertension. Similarly, the expanded investigation of 529 potential SNPs from the PheGenI catalogue failed to show significant overlap with EC-specific genes.

In conclusion, our study does not support that genetic variation in genes expressed in endothelium plays a major role in the development of hypertension. The result further underscores the importance of rigorous validation of genetic associations in large and independent populations.

Acknowledgments

This study makes use of data generated by the Wellcome Trust Case-Control Consortium. A full list of the investigators who contributed to the generation of the data is available from www.wtccc.org.uk.

Author Contributions

Conceived and designed the experiments: EL BW BH TH SEK OM PL. Performed the experiments: EL BH OM. Analyzed the data: EL OM PL. Contributed reagents/materials/analysis tools: EL TH SE BH OM PL. Wrote the paper: EL PL.

Citation: Larsson E, Wahlstrand B, Hedblad B, Hedner T, Kjeldsen SE, Melander O, et al. (2013) Hypertension and Genetic Variation in Endothelial-Specific Genes. PLoS ONE 8(4): e62035. https://doi.org/10.1371/journal.pone.0062035

References

1. Wellcome Trust Case Control Consortium (2007) Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature 447: 661–678.

2. Wellcome Trust Case Control Consortium (2007) Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature 447: 661–678.

3. International Consortium for Blood Pressure Genome-Wide Association Studies, Ehret GB, Munroe PB, Rice KM, Bochud M, et al. (2011) Genetic variants in novel pathways influence blood pressure and cardiovascular disease risk. Nature 478: 103–109.

4. Saxena R, Voight BF, Lyssenko V, Burtt NP, de Bakker PI, et al. (2007) Genome-wide association analysis identifies loci for type 2 diabetes and triglyceride levels. Science 316: 1331–1336.

5. Levy D, Ehret GB, Rice K, Verwoert GC, Launer LJ, et al. (2009) Genome-wide association study of blood pressure and hypertension. Nat Genet 41: 677–687.

6. Newton-Cheh C, Johnson T, Gateva V, Tobin MD, Bochud M, et al. (2009) Genome-wide association study identifies eight loci associated with blood pressure. Nat Genet 41: 666–676.

7. McCarthy MI, Abecasis GR, Cardon LR, Goldstein DB, Little J, et al. (2008) Genome-wide association studies for complex traits: consensus, uncertainty and challenges. Nat Rev Genet 9: 356–369.

8. Asai T, Ohkubo T, Katsuya T, Higaki J, Fu Y, et al. (2001) Endothelin-1 gene variant associates with blood pressure in obese Japanese subjects: the Ohasama Study. Hypertension 38: 1321–1324.

9. Glenn CL, Wang WY, Morris BJ (1999) Different frequencies of inducible nitric oxide synthase genotypes in older hypertensives. Hypertension 33: 927–932.

10. Rutherford S, Johnson MP, Curtain RP, Griffiths LR (2001) Chromosome 17 and the inducible nitric oxide synthase gene in human essential hypertension. Hum Genet 109: 408–415.

11. Tiret L, Poirier O, Hallet V, McDonagh TA, Morrison C, et al. (1999) The Lys198Asn polymorphism in the endothelin-1 gene is associated with blood pressure in overweight people. Hypertension 33: 1169–1174.

12. Wallgard E, Larsson E, He L, Hellstrom M, Armulik A, et al. (2008) Identification of a core set of 58 gene transcripts with broad and specific expression in the microvasculature. Arterioscler Thromb Vasc Biol 28: 1469–1476.

13. Hansson L, Hedner T, Lund-Johansen P, Kjeldsen SE, Lindholm LH, et al. (2000) Randomised trial of effects of calcium antagonists compared with diuretics and beta-blockers on cardiovascular morbidity and mortality in hypertension: the Nordic Diltiazem (NORDIL) study. Lancet 356: 359–365.

14. Kathiresan S, Melander O, Anevski D, Guiducci C, Burtt NP, et al. (2008) Polymorphisms associated with cholesterol and risk of cardiovascular events. N Engl J Med 358: 1240–1249.

15. Hedblad B, Wikstrand J, Janzon L, Wedel H, Berglund G (2001) Low-dose metoprolol CR/XL and fluvastatin slow progression of carotid intima-media thickness: Main results from the Beta-Blocker Cholesterol-Lowering Asymptomatic Plaque Study (BCAPS). Circulation 103: 1721–1726.

16. Berglund G, Eriksson KF, Israelsson B, Kjellstrom T, Lindgarde F, et al. (1996) Cardiovascular risk groups and mortality in an urban swedish male population: the Malmo Preventive Project. J Intern Med 239: 489–497.

17. Rutledge RG, Cote C (2003) Mathematics of quantitative kinetic PCR and the application of standard curves. Nucleic Acids Res 31: e93.

18. Hubbard T, Barker D, Birney E, Cameron G, Chen Y, et al. (2002) The Ensembl genome database project. Nucleic Acids Res 30: 38–41.

19. Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, et al. (2002) The human genome browser at UCSC. Genome Res 12: 996–1006.

20. Tobin MD, Sheehan NA, Scurrah KJ, Burton PR (2005) Adjusting for treatment effects in studies of quantitative traits: antihypertensive therapy and systolic blood pressure. Stat Med 24: 2911–2935.

21. Mailman MD, Feolo M, Jin Y, Kimura M, Tryka K, et al. (2007) The NCBI dbGaP database of genotypes and phenotypes. Nat Genet 39: 1181–1186.

22. Levy D, Larson MG, Benjamin EJ, Newton-Cheh C, Wang TJ, et al. (2007) Framingham Heart Study 100K Project: genome-wide associations for blood pressure and arterial stiffness. BMC Med Genet 8 Suppl 1S3.

23. O’Donnell CJ, Cupples LA, D’Agostino RB, Fox CS, Hoffmann U, et al. (2007) Genome-wide association study for subclinical atherosclerosis in major arterial territories in the NHLBI’s Framingham Heart Study. BMC Med Genet 8 Suppl 1S4.

24. Wang Y, O’Connell JR, McArdle PF, Wade JB, Dorff SE, et al. (2009) From the Cover: Whole-genome association study identifies STK39 as a hypertension susceptibility gene. Proc Natl Acad Sci U S A 106: 226–231.

25. Lowe JK, Maller JB, Pe’er I, Neale BM, Salit J, et al. (2009) Genome-wide association studies in an isolated founder population from the Pacific Island of Kosrae. PLoS Genet 5: e1000365.

26. Org E, Eyheramendy S, Juhanson P, Gieger C, Lichtner P, et al. (2009) Genome-wide scan identifies CDH13 as a novel susceptibility locus contributing to blood pressure determination in two European populations. Hum Mol Genet 18: 2288–2296.

27. Cho YS, Go MJ, Kim YJ, Heo JY, Oh JH, et al. (2009) A large-scale genome-wide association study of Asian populations uncovers genetic factors influencing eight quantitative traits. Nat Genet 41: 527–534.

28. Arnett DK, Li N, Tang W, Rao DC, Devereux RB, et al. (2009) Genome-wide association study identifies single-nucleotide polymorphism in KCNB1 associated with left ventricular mass in humans: the HyperGEN Study. BMC Med Genet 10: 43.

29. Adeyemo A, Gerry N, Chen G, Herbert A, Doumatey A, et al. (2009) A genome-wide association study of hypertension and blood pressure in African Americans. PLoS Genet 5: e1000564.

30. Zabaneh D, Balding DJ (2010) A genome-wide association study of the metabolic syndrome in Indian Asian men. PLoS One 5: e11961.

31. Padmanabhan S, Melander O, Johnson T, Di Blasio AM, Lee WK, et al. (2010) Genome-wide association study of blood pressure extremes identifies variant near UMOD associated with hypertension. PLoS Genet 6: e1001177.

32. Kato N, Takeuchi F, Tabara Y, Kelly TN, Go MJ, et al. (2011) Meta-analysis of genome-wide association studies identifies common variants associated with blood pressure variation in east Asians. Nat Genet 43: 531–538.

33. Wain LV, Verwoert GC, O’Reilly PF, Shi G, Johnson T, et al. (2011) Genome-wide association study identifies six new loci influencing pulse pressure and mean arterial pressure. Nat Genet 43: 1005–1011.

34. Kim SJ, Lee SK, Kim SH, Yun CH, Kim JH, et al. (2012) Genetic association of short sleep duration with hypertension incidence–a 6-year follow-up in the Korean genome and epidemiology study. Circ J 76: 907–913.

35. Guo Y, Tomlinson B, Chu T, Fang YJ, Gui H, et al. (2012) A genome-wide linkage and association scan reveals novel loci for hypertension and blood pressure traits. PLoS One 7: e31489.

36. Crighton D, Wilkinson S, O’Prey J, Syed N, Smith P, et al. (2006) DRAM, a p53-induced modulator of autophagy, is critical for apoptosis. Cell 126: 121–134.

37. Intengan HD, Schiffrin EL (2001) Vascular remodeling in hypertension: roles of apoptosis, inflammation, and fibrosis. Hypertension 38: 581–587.

38. Potter LR, Abbey-Hosch S, Dickey DM (2006) Natriuretic peptides, their receptors, and cyclic guanosine monophosphate-dependent signaling functions. Endocr Rev 27: 47–72.