FSHB -211 G>T is a major genetic modulator of reproductive physiology and health in childbearing age women.

STUDY QUESTION: Are the genetic variants FSHB -211 G>T (rs10835638), FSHR c.2039 A>G (Asn680Ser, rs6166) and FSHR -29 G>A (rs1394205) associated with serum FSH, LH and anti-Müllerian hormone (AMH) levels in reproductive age women, their menstrual cycle parameters and risk of infertility? SUMMARY ANSWER: Only the FSHB -211 G>T variant was a major genetic determinant of serum gonadotropin levels in both, eumenorrheic healthy women and female infertility patients, and the T-allele carrier status was enriched among idiopathic infertility cases. WHAT IS KNOWN ALREADY: There are accumulating data on common genetic variants modulating reproductive parameters and fertility potential. FSHB -211 G>T represents the strongest acknowledged genetic factor contributing to male circulating gonadotropins levels. Respective data in women are limited and the two previously published studies have reached conflicting results. In addition, previous studies have consistently associated FSHR c.2039 A>G (but not FSHR -29 G>A) with female serum FSH level. STUDY DESIGN, SIZE, DURATION: The study aimed to test robust and clinically meaningful genetic effects (if present) of the FSHB -211 G>T, FSHR c.2039 A>G and FSHR -29 G>A variants on female basal FSH, LH and AMH levels, and linked reproductive parameters. Genetic association testing was performed in two independent and clinically different study groups (i) eumenorrheic healthy women without known fertility problems (n = 169; 27.6 ± 6.1 years) and (ii) female partners of infertile couples (n = 186; 32.4 ± 4.7 years). The study groups were compared for allelic and genotypic distributions of the analysed variants. PARTICIPANTS/MATERIALS, SETTING, METHODS: All participants were recruited during the HAPPY PREGNANCY study (2013-2015) at the Women's Clinic, Tartu University Hospital, Estonia. Serum FSH, LH and AMH were measured in the follicular phase (Days 2-6) of the menstrual cycle. All three single nucleotide polymorphisms (SNPs) were genotyped by PCR and Taqman allelic discrimination assay. The effect of the analysed variants on hormonal measurements and menstrual cycle data was assessed using linear regression under additive and recessive models adjusted by age, BMI and smoking status. Results of the two subgroups were combined in a meta-analysis applying the fixed effects model. Restricted maximum likelihood analysis was applied to estimate the proportion of total phenotypic variance of analysed reproductive parameters, explainable by the tested genetic variants. In case-control analysis, genetic association with infertility status was tested using Fisher's exact test and logistic regression adjusted by age, BMI and smoking status. MAIN RESULTS AND THE ROLE OF CHANCE: In both study groups, T-allele of the FSHB -211 G>T was associated with significantly higher serum levels of FSH and LH. Results of the meta-analysis (additive genetic model) remained significant after Bonferroni correction for multiple testing: FSH, T-allele effect 0.80 IU/L, P = 1.2 × 10-3; LH, 1.58 IU/L, P = 1.8×10-8. A more pronounced effect of T-allele of the FSHB -211 G>T on circulating LH was identified as a driving factor to increased LH/FSH ratio (meta-analysis, P = 4.7 × 10-3). In healthy women, the FSHB -211 G>T variant was estimated to explain 3.5 and 7.1% of the total variance of the measured serum FSH and LH levels, respectively. The corresponding numbers for the infertility patients were 1.6 and 10.5%. Women with idiopathic infertility compared to controls exhibited a doubled T-allele frequency (23.6 versus 12.4%; P = 8.9 × 10-3) and a >3-fold excess of TT homozygotes (5.6 versus 1.8%; P = 3.5 × 10-2). The only association of the FSHR c.2039 A>G was detected with serum FSH levels in eumenorrheic healthy women, explaining 3.9% of the total parameter variance (G-allele effect 0.56 IU/L, P = 4.6 × 10-3). In the study group of healthy reproductive age women, the highest serum FSH levels were detected among the FSHB -211 T-allele carriers with the FSHR c.2039 GG-genotype (median 7.7 IU/L). In contrast, the lowest hormone concentrations were measured for the women carrying the combination of the FSHB -211 GG- and the FSHR c.2039 AA-homozygosity (median 5.8 IU/L, P = 9.6 × 10-3). None of the analysed reproductive parameters was associated with the FSHR -29 G>A variant. In our study groups, the tested polymorphisms did not reach significant associations with serum AMH measurements, menstrual cycle length or age at menarche. LIMITATIONS, REASONS FOR CAUTION: Small sample size and the design involving two clinical groups with different reproductive histories may have limited the capacity to replicate the associations with the age at menarche and length of menstrual cycle, initially reported in large genome-wide association studies. Small sample size may have also affected the accuracy in estimating the contribution of the tested variants to the total phenotypic variance of measured gonadotropin concentrations. The group of eumenorrheic healthy women had its limitations as a control to estimate the true effect of analysed genetic variants on individual's fertility potential as the recruitment strategy had been targeted mostly towards younger women, who may not yet have planned to conceive a child by this age. WIDER IMPLICATIONS OF THE FINDINGS: We propose that like in men, also in women the FSHB -211 G>T represents a key genetic modulator of circulating gonadotropin, leading to various possible downstream effects on reproductive physiology. This claim is strongly supported by the reports of genome-wide association studies on various female reproductive traits and diseases. In perspective, FSHB -211 G>T may have a diagnostic value in fertility clinics to detect female patients with genetically inherited elevated basal FSH and LH levels. STUDY FUNDING/COMPETING INTEREST(S): The study was supported by Estonian Science Foundation Grant (ETF9030 for M.L.); Institutional Research Grant (IUT34-12 for M.L.) and European Union through the European Regional Development Fund (project HAPPY PREGNANCY, 3.2.0701.12-0047; for M.L. and K.R.). The funders had no role in study design, data collection and analysis, decision to publish or preparation of the article. We have no competing interests to declare. TRAIL REGISTRATION NUMBER: Not applicable.

Blood pressure loci identified with a gene-centric array.

Raised blood pressure (BP) is a major risk factor for cardiovascular disease. Previous studies have identified 47 distinct genetic variants robustly associated with BP, but collectively these explain only a few percent of the heritability for BP phenotypes. To find additional BP loci, we used a bespoke gene-centric array to genotype an independent discovery sample of 25,118 individuals that combined hypertensive case-control and general population samples. We followed up four SNPs associated with BP at our p < 8.56 × 10(-7) study-specific significance threshold and six suggestively associated SNPs in a further 59,349 individuals. We identified and replicated a SNP at LSP1/TNNT3, a SNP at MTHFR-NPPB independent (r(2) = 0.33) of previous reports, and replicated SNPs at AGT and ATP2B1 reported previously. An analysis of combined discovery and follow-up data identified SNPs significantly associated with BP at p < 8.56 × 10(-7) at four further loci (NPR3, HFE, NOS3, and SOX6). The high number of discoveries made with modest genotyping effort can be attributed to using a large-scale yet targeted genotyping array and to the development of a weighting scheme that maximized power when meta-analyzing results from samples ascertained with extreme phenotypes, in combination with results from nonascertained or population samples. Chromatin immunoprecipitation and transcript expression data highlight potential gene regulatory mechanisms at the MTHFR and NOS3 loci. These results provide candidates for further study to help dissect mechanisms affecting BP and highlight the utility of studying SNPs and samples that are independent of those studied previously even when the sample size is smaller than that in previous studies.

Novel polymorphic AluYb8 insertion in the WNK1 gene is associated with blood pressure variation in Europeans.

Mutations in WNK1 and WNK4 cause familial hypertension, the Gordon syndrome. WNK1 and WNK4 conserved noncoding regions were targeted to polymorphism screening using DHPLC and DGGE. The scan identified an undescribed polymorphic AluYb8 insertion in WNK1 intron 10. Screening in primates revealed that this Alu-insertion has probably occurred in human lineage. Genotyping in 18 populations from Europe, Asia, and Africa (n = 854) indicated an expansion of the WNK1 AluYb8 bearing chromosomes out of Africa. The allele frequency in Sub-Saharan Africa was ~3.3 times lower than in other populations (4.8 vs. 15.8%; P = 9.7 × 10(-9) ). Meta-analysis across three European sample sets (n = 3,494; HYPEST, Estonians; BRIGHT, the British; CADCZ, Czech) detected significant association of the WNK1 AluYb8 insertion with blood pressure (BP; systolic BP, P = 4.03 × 10(-3) , effect 1.12; diastolic BP, P = 1.21 × 10(-2) , effect 0.67). Gender-stratified analysis revealed that this effect might be female-specific (n = 2,088; SBP, P = 1.99 × 10(-3) , effect 1.59; DBP P = 3.64 × 10(-4) , effect 1.23; resistant to Bonferroni correction), whereas no statistical support was identified for the association with male BP (n = 1,406). In leucocytes, the expressional proportions of the full-length WNK1 transcript and the splice-form skipping exon 11 were significantly shifted in AluYb8 carriers compared to noncarriers. The WNK1 AluYb8 insertion might affect human BP via altering the profile of alternatively spliced transcripts.

Genome-wide scan identifies CDH13 as a novel susceptibility locus contributing to blood pressure determination in two European populations.

Hypertension is a complex disease that affects a large proportion of adult population. Although approximately half of the inter-individual variance in blood pressure (BP) level is heritable, identification of genes responsible for its regulation has remained challenging. Genome-wide association study (GWAS) is a novel approach to search for genetic variants contributing to complex diseases. We conducted GWAS for three BP traits [systolic and diastolic blood pressure (SBP and DBP); hypertension (HYP)] in the Kooperative Gesundheitsforschung in der Region Augsburg (KORA) S3 cohort (n = 1644) recruited from general population in Southern Germany. GWAS with 395,912 single nucleotide polymorphisms (SNPs) identified an association between BP traits and a common variant rs11646213 (T/A) upstream of the CDH13 gene at 16q23.3. The initial associations with HYP and DBP were confirmed in two other European population-based cohorts: KORA S4 (Germans) and HYPEST (Estonians). The associations between rs11646213 and three BP traits were replicated in combined analyses (dominant model: DBP, P = 5.55 x 10(-5), effect -1.40 mmHg; SBP, P = 0.007, effect -1.56 mmHg; HYP, P = 5.30 x 10(-8), OR = 0.67). Carriers of the minor allele A had a decreased risk of hypertension. A non-significant trend for association was also detected with severe family based hypertension in the BRIGHT sample (British). The novel susceptibility locus, CDH13, encodes for an adhesion glycoprotein T-cadherin, a regulator of vascular wall remodeling and angiogenesis. Its function is compatible with the BP biology and may improve the understanding of the pathogenesis of hypertension.

HYPEST study: profile of hypertensive patients in Estonia.

BACKGROUND: More than one third of adult population in Estonia has problems with elevated blood pressure (BP). The Hypertension in Estonia (HYPEST) study represents the country's first hypertension-targeted sample collection aiming to examine the epidemiological and genetic determinants for hypertension (HTN) and related cardiovascular diseases (CVD) in Estonian population. The HYPEST subjects (n = 1,966) were recruited across Estonia between 2004-2007 including clinically diagnosed HTN cases and population-based controls. The present report is focused on the clinical and epidemiological profile of HYPEST cases, and gender-specific effects on the pathophysiology of hypertension. METHODS: Current analysis was performed on 1,007 clinically diagnosed HTN patients (617 women and 390 men) aged 18-85 years. The hypertensives were recruited to the study by BP specialists at the North Estonia Medical Center, Centre of Cardiology, Tallinn or at the Cardiology Clinic, Tartu University Hospital, Estonia. Longitudinal BP data was extracted retrospectively from clinical records. Current and retrospective data of patient's medical history, medication intake and lifestyle habits were derived from self-administrated questionnaire and each variable was examined separately for men and women. Eleven biochemical parameters were measured from fasting serum samples of 756 patients. RESULTS: The distribution of recruited men and women was 39% and 61% respectively. Majority of Estonian HTN patients (85%) were overweight (BMI ≥ 25 kg/m2) and a total of 79% of patients had additional complications with cardiovascular system. In men, the hypertension started almost 5 years earlier than in women (40.5 ± 14.5 vs 46.1 ± 12.7 years), which led to earlier age of first myocardial infarction (MI) and overall higher incidence rate of MI among male patients (men 21.2%, women 8.9%, P < 0.0001). Heart arrhythmia, thyroid diseases, renal tubulo-intestinal diseases and hyperlipidemia were more prevalent in hypertensive women compared to men (P < 0.0001). An earlier age of HTN onset was significantly associated with smoking (P = 0.00007), obesity (BMI ≥ 30 kg/m2; P = 0.0003), increased stress (P = 0.0003) and alcohol consumption (P = 0.004). CONCLUSION: Understanding the clinical profile of HTN patients contributes to CVD management. Estonian hypertension patients exhibited different disease and risk profiles of male and female patients. This well-characterized sample set provides a good resource for studying hypertension and other cardiovascular phenotypes.

Hypervariable intronic region in NCX1 is enriched in short insertion-deletion polymorphisms and showed association with cardiovascular traits.

BACKGROUND: Conserved non-coding regions (CNR) have been shown to harbor gene expression regulatory elements. Genetic variations in these regions may potentially contribute to complex disease susceptibility. METHODS: We targeted CNRs of cardiovascular disease (CVD) candidate gene, Na(+)-Ca(2+) exchanger (NCX1) with polymorphism screening among CVD patients (n = 46) using DHPLC technology. The flanking region (348 bp) of the 14 bp indel in intron 2 was further genotyped by DGGE assay in two Eastern-European CVD samples: essential hypertension (HYPEST; 470 cases, 652 controls) and coronary artery disease, CAD (CADCZ; 257 cases, controls 413). Genotype-phenotype associations were tested by regression analysis implemented in PLINK. Alignments of primate sequences were performed by ClustalW2. RESULTS: Nine of the identified NCX1 variants were either singletons or targeted by commercial platforms. The 14 bp intronic indel (rs11274804) was represented with substantial frequency in HYPEST (6.82%) and CADCZ (14.58%). Genotyping in Eastern-Europeans (n = 1792) revealed hypervariable nature of this locus, represented by seven alternative alleles. The alignments of human-chimpanzee-macaque sequences showed that the major human variant (allele frequency 90.45%) was actually a human-specific deletion compared to other primates. In humans, this deletion was surrounded by other short (5-43 bp) deletion variants and a duplication (40 bp) polymorphism possessing overlapping breakpoints. This indicates a potential indel hotspot, triggered by the initial deletion in human lineage. An association was detected between the carrier status of 14 bp indel ancestral allele and CAD (P = 0.0016, OR = 2.02; Bonferroni significance level alpha = 0.0045), but not with hypertension. The risk for the CAD development was even higher among the patients additionally diagnosed with metabolic syndrome (P = 0.0014, OR = 2.34). Consistent with the effect on metabolic processes, suggestive evidence for the association with heart rate, serum triglyceride and LDL levels was detected (P = 0.04). CONCLUSIONS: Compared to SNPs targeted by large number of locus-specific and genome-wide assays, considerably less attention has been paid to short indel variants in the human genome. The data of genome dynamics, mutation rate and population genetics of short indels, as well as their impact on gene expressional profile and human disease susceptibility is limited. The characterization of NCX1 intronic hypervariable non-coding region enriched in human-specific indel variants contributes to this gap of knowledge.

N-acetyltransferase 8, a positional candidate for blood pressure and renal regulation: resequencing, association and in silico study.

BACKGROUND: Kidneys have an important function in blood pressure (BP) regulation and elevated BP may lead to kidney failure. Chr2p12-p13 region linked to BP traits in multiple studies harbours a potential candidate for BP and renal function, N-acetyltransferase 8 (NAT8) expressed in embryonic and adult kidney and associated with nephrotoxicity response. METHODS/RESULTS: We report the first study exploring NAT8 as a potential candidate gene for blood pressure and kidney function. The resequencing (n = 42, random Estonian samples) identified 15 NAT8 polymorphisms, including 6 novel variants. The diversity of NAT8 5' upstream region (pi/bp = 0.00320) exceeded up to 10 times the variation in the NAT8 genic region (pi/bp = 0.00037) as well as the average variation (pi/bp = 0.00040) for the promoters of 29 reference genes associated with hypertension. We suggest that a potential source for such high variation could be an active gene conversion process from NAT8B duplicate gene to NAT8. Similarly to NAT8, several reference genes with the most variable upstream regions have also duplicate copies. The NAT8 promoter SNPs were targeted with pilot quantitative association studies for blood pressure (n = 137, healthy unrelated individuals) and for the index of kidney function - estimated glomerular filtration rate (eGFR; n = 157 hypertensives with and without nephropathy). Minor alleles of these polymorphisms revealed a significant protective effect against elevated systolic BP as well as kidney failure in hypertension patients (p < 0.05; linear regression model, addictive effect). CONCLUSION: The full resequencing and pilot association study of a novel positional candidate gene for blood pressure and renal function, human N-acetyltransferase 8, suggested a contribution of highly variable NAT8 promoter polymorphisms in determination of systolic blood pressure and eGFR. Based on in silico analysis, we raise the hypothesis that the alternative SNP alleles of the NAT8 upstream region may have differential effect on gene expression.

Resequencing PNMT in European hypertensive and normotensive individuals: no common susceptibilily variants for hypertension and purifying selection on intron 1.

BACKGROUND: Human linkage and animal QTL studies have indicated the contribution of genes on Chr17 into blood pressure regulation. One candidate gene is PNMT, coding for phenylethanolamine-N-methyltransferase, catalyzing the synthesis of epinephrine from norepinephrine. METHODS: Fine-scale variation of PNMT was screened by resequencing hypertensive (n = 50) and normotensive (n = 50) individuals from two European populations (Estonians and Czechs). The resulting polymorphism data were analyzed by statistical genetics methods using Genepop 3.4, PHASE 2.1 and DnaSP 4.0 software programs. In silico prediction of transcription factor binding sites for intron 1 was performed with MatInspector 2.2 software. RESULTS: PNMT was characterized by minimum variation and excess of rare SNPs in both normo- and hypertensive individuals. None of the SNPs showed significant differences in allelic frequencies among population samples, as well as between screened hypertensives and normotensives. In the joint case-control analysis of the Estonian and the Czech samples, hypertension patients had a significant excess of heterozygotes for two promoter region polymorphisms (SNP-184; SNP-390). The identified variation pattern of PNMT reflects the effect of purifying selection consistent with an important role of PNMT-synthesized epinephrine in the regulation of cardiovascular and metabolic functions, and as a CNS neurotransmitter. A striking feature is the lack of intronic variation. In silico analysis of PNMT intron 1 confirmed the presence of a human-specific putative Glucocorticoid Responsive Element (GRE), inserted by Alu-mediated transfer. Further analysis of intron 1 supported the possible existence of a full Glucocorticoid Responsive Unit (GRU) predicted to consist of multiple gene regulatory elements known to cooperate with GRE in driving transcription. The role of these elements in regulating PNMT expression patterns and thus determining the dynamics of the synthesis of epinephrine is still to be studied. CONCLUSION: We suggest that the differences in PNMT expression between normotensives and hypertensives are not determined by the polymorphisms in this gene, but rather by the interplay of gene expression regulators, which may vary among individuals. Understanding the determinants of PNMT expression may assist in developing PNMT inhibitors as potential novel therapeutics.

Stanniocalcin-1 Hormone in Nonpreeclamptic and Preeclamptic Pregnancy: Clinical, Life-Style, and Genetic Modulators.

CONTEXT AND OBJECTIVES: The study represents the first comprehensive analysis of Stanniocalcin-1 (STC1) hormone in human pregnancy, assessing clinical, lifestyle, and genetic determinants of circulating STC1 at term. DESIGN, SETTING, AND PARTICIPANTS: Participants included women with (n = 50) and without (n = 316) preeclampsia (PE) at delivery, recruited in the REPROgrammed fetal and/or maternal METAbolism (REPROMETA) study (2006-2011, Estonia). Genetic association analysis combined PE cases (n = 597) and controls (n = 623) from the REPROMETA and Finnish Genetics of Preeclampsia Consortium (2008-2011) studies. MAIN OUTCOME MEASURE(S): Maternal postpartum plasma STC1 was measured by ELISA (n = 366) and placental STC1 gene expression by TaqMan quantitative RT-PCR (n = 120). Genotyping was performed using Sequenom MassArray. RESULTS: Significantly higher STC1 plasma level was measured for the PE (median, 1952 pg/mL; 1030-4284 pg/mL) compared with non-PE group (median, 1562 pg/mL; 423-3781 pg/mL; P = 3.7 × 10-4, Mann-Whitney U test). Statistical significance was enhanced after adjustment for cofactors (linear regression, P = 1.8 × 10-6). STC1 measurements were negatively correlated with maternal smoking. Prepregnancy body mass index had a positive correlation with STC1 only among PE patients (r = 0.45; P = .001). The strongest genetic association with hormone concentrations was detected for STC1 single nucleotide polymorphisms rs3758089 (C allele: minor allele frequency, 5%; linear regression: ß = 249.2 pg/mL; P = .014) and rs12678447 (G allele: minor allele frequency, 7%; ß = 147.0 pg/mL; P = .082). rs12678447 placental genotypes were significantly associated with STC1 gene expression (P = .014). The REPROMETA/Finnish Genetics of Preeclampsia Consortium meta-analysis suggested an increased risk to develop late-onset PE for the rs12678447 G allele carriers (P = .05; odds ratio = 1.38 [0.98-1.93]). CONCLUSIONS: Increased STC1 hormone represents a hallmark of late-onset PE. STC1 gene variants modulate placental gene expression and maternal hormone levels.

Targeting 160 candidate genes for blood pressure regulation with a genome-wide genotyping array.

The outcome of Genome-Wide Association Studies (GWAS) has challenged the field of blood pressure (BP) genetics as previous candidate genes have not been among the top loci in these scans. We used Affymetrix 500K genotyping data of KORA S3 cohort (n = 1,644; Southern-Germany) to address (i) SNP coverage in 160 BP candidate genes; (ii) the evidence for associations with BP traits in genome-wide and replication data, and haplotype analysis. In total, 160 gene regions (genic region+/-10 kb) covered 2,411 SNPs across 11.4 Mb. Marker densities in genes varied from 0 (n = 11) to 0.6 SNPs/kb. On average 52.5% of the HAPMAP SNPs per gene were captured. No evidence for association with BP was obtained for 1,449 tested SNPs. Considerable associations (P<10(-3)) were detected for the genes, where >50% of HAPMAP SNPs were tagged. In general, genes with higher marker density (>0.2 SNPs/kb) revealed a better chance to reach close to significance associations. Although, none of the detected P-values remained significant after Bonferroni correction (P<0.05/2319, P<2.15 x 10(-5)), the strength of some detected associations was close to this level: rs10889553 (LEPR) and systolic BP (SBP) (P = 4.5 x 10(-5)) as well as rs10954174 (LEP) and diastolic BP (DBP) (P = 5.20 x 10(-5)). In total, 12 markers in 7 genes (ADRA2A, LEP, LEPR, PTGER3, SLC2A1, SLC4A2, SLC8A1) revealed considerable association (P<10(-3)) either with SBP, DBP, and/or hypertension (HYP). None of these were confirmed in replication samples (KORA S4, HYPEST, BRIGHT). However, supportive evidence for the association of rs10889553 (LEPR) and rs11195419 (ADRA2A) with BP was obtained in meta-analysis across samples stratified either by body mass index, smoking or alcohol consumption. Haplotype analysis highlighted LEPR and PTGER3. In conclusion, the lack of associations in BP candidate genes may be attributed to inadequate marker coverage on the genome-wide arrays, small phenotypic effects of the loci and/or complex interaction with life-style and metabolic parameters.

Polymorphisms in the WNK1 gene are associated with blood pressure variation and urinary potassium excretion.

WNK1--a serine/threonine kinase involved in electrolyte homeostasis and blood pressure (BP) control--is an excellent candidate gene for essential hypertension (EH). We and others have previously reported association between WNK1 and BP variation. Using tag SNPs (tSNPs) that capture 100% of common WNK1 variation in HapMap, we aimed to replicate our findings with BP and to test for association with phenotypes relating to WNK1 function in the British Genetics of Hypertension (BRIGHT) study case-control resource (1700 hypertensive cases and 1700 normotensive controls). We found multiple variants to be associated with systolic blood pressure, SBP (7/28 tSNPs min-p = 0.0005), diastolic blood pressure, DBP (7/28 tSNPs min-p = 0.002) and 24 hour urinary potassium excretion (10/28 tSNPs min-p = 0.0004). Associations with SBP and urine potassium remained significant after correction for multiple testing (p = 0.02 and p = 0.01 respectively). The major allele (A) of rs765250, located in intron 1, demonstrated the strongest evidence for association with SBP, effect size 3.14 mmHg (95%CI:1.23-4.9), DBP 1.9 mmHg (95%CI:0.7-3.2) and hypertension, odds ratio (OR: 1.3 [95%CI: 1.0-1.7]).We genotyped this variant in six independent populations (n = 14,451) and replicated the association between rs765250 and SBP in a meta-analysis (p = 7 x 10(-3), combined with BRIGHT data-set p = 2 x 10(-4), n = 17,851). The associations of WNK1 with DBP and EH were not confirmed. Haplotype analysis revealed striking associations with hypertension and BP variation (global permutation p<10(-7)). We identified several common haplotypes to be associated with increased BP and multiple low frequency haplotypes significantly associated with lower BP (>10 mmHg reduction) and risk for hypertension (OR<0.60). Our data indicates that multiple rare and common WNK1 variants contribute to BP variation and hypertension, and provide compelling evidence to initiate further genetic and functional studies to explore the role of WNK1 in BP regulation and EH.