==== Front Circ Cardiovasc GenetCirc Cardiovasc GenetHCGCirculation. Cardiovascular Genetics1942-325X1942-3268Lippincott Williams & Wilkins 0000910.1161/CIRCGENETICS.117.00171310015100581008210111Original ArticlesGenome-Wide Association Study Implicates Atrial Natriuretic Peptide Rather Than B-Type Natriuretic Peptide in the Regulation of Blood Pressure in the General Population Salo Perttu P. MScHavulinna Aki S. PhDTukiainen Taru PhDRaitakari Olli MD, PhDLehtimäki Terho MD, PhDKähönen Mika MD, PhDKettunen Johannes PhDMännikkö Minna PhDEriksson Johan G. MD, PhDJula Antti MD, PhDBlankenberg Stefan MD, PhDZeller Tanja PhDSalomaa Veikko MD, PhDKristiansson Kati PhD*Perola Markus MD, PhD*Correspondence to Perttu Salo, MSc, Institute for Molecular Medicine Finland, Biomedicum 1, Haartmaninkatu 8 rm. A331b, 00290 Helsinki, Finland. E-mail perttu.salo@thl.fi12 2017 13 12 2017 10 6 e001713 e001713 23 1 2017 3 10 2017 © 2017 The Authors.2017Circulation: Cardiovascular Genetics is published on behalf of the American Heart Association, Inc., by Wolters Kluwer Health, Inc. This is an open access article under the terms of the Creative Commons Attribution Non-Commercial-NoDerivs License, which permits use, distribution, and reproduction in any medium, provided that the original work is properly cited, the use is noncommercial, and no modifications or adaptations are made.Supplemental Digital Content is available in the text. Background— Cardiomyocytes secrete atrial natriuretic peptide (ANP) and B-type natriuretic peptide (BNP) in response to mechanical stretching, making them useful clinical biomarkers of cardiac stress. Both human and animal studies indicate a role for ANP as a regulator of blood pressure with conflicting results for BNP. Methods and Results— We used genome-wide association analysis (n=6296) to study the effects of genetic variants on circulating natriuretic peptide concentrations and compared the impact of natriuretic peptide–associated genetic variants on blood pressure (n=27 059). Eight independent genetic variants in 2 known (NPPA-NPPB and POC1B-GALNT4) and 1 novel locus (PPP3CC) associated with midregional proANP (MR-proANP), BNP, aminoterminal proBNP (NT-proBNP), or BNP:NT-proBNP ratio. The NPPA-NPPB locus containing the adjacent genes encoding ANP and BNP harbored 4 independent cis variants with effects specific to either midregional proANP or BNP and a rare missense single nucleotide polymorphism in NT-proBNP seriously altering its measurement. Variants near the calcineurin catalytic subunit gamma gene PPP3CC and the polypeptide N-acetylgalactosaminyltransferase 4 gene GALNT4 associated with BNP:NT-proBNP ratio but not with BNP or midregional proANP, suggesting effects on the post-translational regulation of proBNP. Out of the 8 individual variants, only those correlated with midregional proANP had a statistically significant albeit weak impact on blood pressure. The combined effect of these 3 single nucleotide polymorphisms also associated with hypertension risk (P=8.2×10−4). Conclusions— Common genetic differences affecting the circulating concentration of ANP associated with blood pressure, whereas those affecting BNP did not, highlighting the blood pressure–lowering effect of ANP in the general population. blood pressuregenesgenome-wide association studyhypertensionnatriuretic peptide, brainOPEN-ACCESSTRUE ==== Body The heart secretes atrial natriuretic peptide (ANP) and B-type natriuretic peptide (BNP) into the circulation in response to myocardial stretching. Atrial cardiomyocytes mainly secrete ANP, whereas ventricular cardiomyocytes predominantly produce BNP. Coded by the adjacent genes NPPA and NPPB in humans, the proANP and proBNP prohormones are cleaved to produce an inactive N-terminal fragment and the active hormone. ANP and BNP reduce cardiac load via increased natriuresis, vasorelaxation, and other physiological effects mediated by the natriuretic peptide receptor A.1 Both the active hormones and the N-terminal fragments may be used as biomarkers of cardiac stress. A particularly valuable clinical application is the use of low measured BNP or NT-proBNP (N-terminal proBNP) concentration to rule out suspected heart failure.2 ANP and BNP are, thus, regulators of cardiovascular function and useful clinical biomarkers. See Editorial by Armando See Clinical Perspective Natriuretic peptides are attractive therapeutic targets. Overexpression of either NPPA or NPPB in mice leads to pronounced hypotension.3,4 Deleting NPPA in mice predisposes them to hypertension, but knocking out NPPB triggers cardiac fibrosis instead of inducing hypertension.5–7 In contrast to mice, the deletion of NPPB in a hypertensive rat model has been reported to decrease survival and increase both systolic and diastolic blood pressure (BP).8 In humans, the effects of ANP or BNP infusions depend on baseline status. Infusions of ANP or BNP in patients having heart failure trigger various hemodynamic changes, including a decrease in arterial pressure, but in healthy males only induce natriuresis without affecting arterial pressure.9–12 Both ANP and BNP have a BP-lowering effect in those having essential hypertension, with BNP surprisingly showing a 2- to 3-fold greater potency than ANP despite similar receptor affinity.13,14 A lack of association or even a paradoxical negative association of ANP with BP has been reported in obese men.15 How heart failure, hypertension, or obesity may modify ANP and BNP function is incompletely understood. Recombinant BNP has also failed to show a clear clinical benefit in treating acute decompensated heart failure when used in addition to standard care.9 A more detailed understanding of ANP, BNP, and their physiological role may aid in successfully exploiting their potential. Genetic studies of ANP and BNP in humans are of particular interest as a large part of the literature regard knockout animal models and relatively high doses of intravenous infusions. Data on variation in their concentration within the normal physiological range are more scarce but necessary to understand the function of these peptides under nondiseased conditions. The association of human genetic variation with circulating ANP and BNP has been studied for selected single nucleotide polymorphisms (SNPs).16,17 Four genome-wide association studies (GWAS) have studied circulating BNP or NT-proBNP levels.18–21 The prior studies have associated variants near NPPA-NPPB with proANP, BNP, and NT-proBNP, and the GWAS have associated trans loci near LOXL2, SLC39A8, KLKB1, and GALNT4 with NT-proBNP. No genome-wide studies have been published on ANP. The prior studies, thus, either did not have genome-wide coverage of genetic variation or did not assay ANP, limiting the interpretation of their results. We performed genome-wide association tests of BNP, NT-proBNP, and midregional proANP (MR-proANP) and studied the impact of the natriuretic peptide–associated genetic variants on BP. Because proBNP is processed peripherally into BNP and NT-proBNP that have different circulating half-lives, we also studied the ratio of BNP to NT-proBNP concentrations (BNP:NT-proBNP ratio) as a potential proxy for the processing and degradation of BNP, proBNP, and NT-proBNP.1,22 Materials and Methods MR-proANP, NT-proBNP, and BNP were measured in the GWAS discovery (n=4932) and replication samples (n=1373), originally recruited for the FINRISK 1997 study. The National FINRISK Study cohorts are collected every 5 years as representative age- and sex-stratified samples of the populations of 5 geographical areas of Finland, described in more detail elsewhere.23,24 We tested the association of genetic variants with natriuretic peptide traits in the GWAS discovery and replication samples excluding participants who had prevalent diabetes mellitus, heart failure, stroke, or coronary heart disease. We then studied the BP associations of the genetic variants detected in the GWAS in an independent study population, comprising the FINRISK 1992 (n=4920), FINRISK 2002 (n=5,21), FINRISK 2007 (n=4996), the Northern Finland Birth Cohort 1966 (NFBC66, n=5363), the HBCS (Helsinki Birth Cohort Study, n=1619), the YFS (Young Finns Study, n=2443), and the Health2000 (n=1997) cohorts.23–28 All study cohorts were population-based samples of Finns, approved by their respective institutional review committees, and participants gave their informed consent. Natriuretic Peptide and BP Measurements Natriuretic peptide concentrations were measured in the MORGAM Biomarker Laboratory, University of Mainz, Germany, using the Abbott Architect i2000 BNP (BNP, UniProt acc. P16860, residues 103–134), Roche Elecsys 2010 proBNP (NT-proBNP, acc. P16860 residues 27–134), and B.R.A.H.M.S. MR-proANP KRYPTOR (MR-proANP, acc. P01160) assays, described in more detail previously.29 The inter/intra-assay coefficients of variation were 2.11%/4.28% (BNP), 2.58%/1.38% (NT-proBNP), and 3.65%/2.33% (MR-proANP). BP was measured from the participants’ right arm, and hypertension was defined as diastolic BP >90 mm Hg or systolic BP >140 mm Hg or known use of antihypertensive medication. Genotyping and Imputation The GWAS discovery sample and replication samples were genotyped using the Illumina HumanCoreExome beadchip at the Wellcome Trust Sanger Institute (Cambridge, UK) and at the Broad Institute of Harvard and MIT (MA, USA), respectively. The data were prephased and imputed using the 1000 Genomes project phase 1 and 3 haplotypes and a custom haplotype set of 2000 Finnish individuals. After quality control (Hardy–Weinberg equilibrium P value <0.01, minor allele frequency <1%, imputation quality <0.9, genotyping success rate <95%) and removal of rare SNPs (minor allele frequency <1%), the discovery phase GWAS data set contained a total of 7 358 451 SNPs and 4932 samples. Cohorts comprising the BP study population were genotyped on various genome-wide genotyping arrays and imputed using the same methods as used for the GWAS discovery sample (Data Supplement). All genomic coordinates are given using the GRCh37 human reference genome. Association Tests We used multiple imputation to account for any missing values for MR-proANP (Nmissing=0), BNP (Nmissing =131), and NT-proBNP (Nmissing=133) and random-effects meta-analysis to combine results from the different cohorts.30 We inverse-normal transformed the natriuretic peptide measurements and used linear regression with an additive genetic model adjusted for geographical sampling region, age2 sex, body mass index (BMI), current smoking (yes/no), systolic BP, glomerular filtration rate estimated using cystatin C and creatinine as proxies, and genotyping batch. We used least absolute shrinkage and selection operator regression implemented in the LLARMA package for fine-mapping the natriuretic peptide–associated loci to identify possible secondary independent variants.31 The genetic association tests are described in more detail in the Data Supplement. We used linear regression implemented in the glm function for R to study the association of genetic variants with systolic and diastolic BP. We log-transformed systolic (but not diastolic) BP and set the first 2 genomic principal components, age, sex, BMI, current BP medication use (yes/no, only for systolic and diastolic BP), study year, and genotyping batch as covariates (the latter 2 only for the FINRISK samples). For hypertension, we used logistic regression and the same covariates excluding BP medication. Phenotypic Variance Explained by SNPs Genome-Wide We used autosomal SNPs from the imputed data set to estimate the fraction of phenotypic variance explained by the SNPs genome-wide in the participants of the GWAS discovery sample using PLINK v1.90 and GCTA v1.25.3.32,33 As a quality control measure, we derived 4 genomic scores corresponding to each of the 4 estimates (for MR-proANP, BNP, NT-proBNP, and BNP:NT-proBNP ratio) and tested the association of the genomic scores with their respective phenotypes in the replication sample (Data Supplement). Coassociation With Gene Expression We investigated the coassociation of SNPs with both natriuretic peptides and gene expression in data from 190 left ventricular tissue samples and 159 atrial appendage samples from the Genotype-Tissue Expression (GTEx) consortium (release V6, October 6, 2016).34 We used 3 metrics to confirm that the same genetic variants correlated with both gene expression and natriuretic peptide concentrations in a consistent way (Data Supplement): we required that the most statistically significant natriuretic peptide–associated SNPs (lead SNPs) associated with the genes’ expression levels and that both the association P values and the effect estimates (βs) were correlated across the SNPs in the natriuretic peptide–associated regions. Because P values depend on allele frequencies, we used both Spearman rank (for P values) and Pearson product moment (for βs) correlation coefficients as measures of the correlation between the natriuretic peptide and gene expression associations and derived the P values empirically. Results Baseline Characteristics The baseline characteristics of the GWAS discovery sample, the replication sample, and the BP study population are described in Table I in the Data Supplement. The strata were broadly similar, and the main difference was that participants with prevalent cardiovascular disease were not excluded from the BP study population. GWAS and Variance Explained by All SNPs To quantify the total amount of genetic signal present in the data, we first estimated the proportion of variance in the natriuretic peptide traits jointly explained by all SNPs genome-wide. The point estimates were 13.9% for MR-proANP, 13.5% for BNP, 23.0% for NT-proBNP, and 17.9% for BNP:NT-proBNP ratio, but the coarse precision of the estimates prevents ranking the 4 phenotypes in any particular order in terms of variance explained (Table II in the Data Supplement). The magnitude of the 4 estimates nonetheless indicates that the SNPs together explained a moderate proportion of the phenotypic variance. Having estimated the proportion of variance explained by all SNPs genome-wide, we tested the SNPs individually for association with the phenotypes. Variants in 4 loci near NPPA-NPPB, PPP3CC, GALNT4, and NCOR12 met the prespecified threshold for genome-wide significance P<5×10−8 for association (Figure 1; Table 1; Figure I in the Data Supplement; Table III in the Data Supplement). We selected the SNP with the smallest P value (lead SNP) at each locus for replication. Only the association of rs701041 with MR-proANP near NCOR12 did not replicate (P=0.94). Associations near NPPA-NPPB and GALNT4 have been reported previously, whereas the association of rs7000551 with BNP:NT-proBNP ratio on chromosome 8 near PPP3CC is a novel finding.16–20 Fine-mapping the loci using least absolute shrinkage and selection operator regression identified independent secondary signals near NPPA-NPPB and GALNT4. Three independent SNPs near NPPA associated with MR-proANP levels, whereas 2 independent SNPs near GALNT4 associated with BNP:NT-proBNP ratio. Table 1. Association of Genetic Variants With Natriuretic Peptides in the Genome-Wide Significant Loci Figure 1. Genome-wide association study P values. P values of the genome-wide association tests and their genomic locations. Y axis cut at Y=18, the peak on chromosome 1 extends to Y=80. Previously detected associations replicated successfully in the present data in terms of the direction of association (Table IV in the Data Supplement). Of these, all but 1 of the cis associations near NPPA-NPPB also reached statistical significance. Two of the 3 previously published trans associations, rs13107325 in SLC39A8 and rs3733402 in KLKB1, associated with BNP:NT-proBNP ratio in the meta-analysis of the discovery and replication samples (rs13107325 P=2.19×10−9; rs3733402 P=0.00277) and the meta-analysis P value of rs13107325 with NT-proBNP (P=0.00496) was also nominally significant. The third, rs6557662 in LOXL2, did not reach statistical significance. None of the trans loci associated with BNP or MR-proANP. Most common variants are thought to affect phenotypes by altering gene expression.35,36 We, thus, studied data from 190 left ventricular tissue samples and 159 atrial appendage samples from the GTEx consortium to identify coassociation of SNPs with both natriuretic peptide traits and gene expression.34 The results of these tests, together with those of the fine-mapping tests with least absolute shrinkage and selection operator regression, are presented in detail below for the loci meeting genome-wide significance in the present study. NPPA-NPPB on Chromosome 1 SNPs associating with the natriuretic peptide on chromosome 1 were located near the NPPA and NPPB genes (Figure 2; Figure II in the Data Supplement). Previously, associations in this locus have been reported using a GWAS strategy for NT-proBNP and a candidate SNP approach for ANP.16–20 To extend the previously reported results, we focus here on the extensive panel of SNPs and the more detailed phenotyping, which were not available in the prior studies. Figure 2. SNPs associated with natriuretic peptides on chromosome 1 near NPPA and NPPB. Linkage disequilibrium on chromosome 1 near NPPA and NPPB. R2 and D′ calculated with Haploview v4.2 from 22 374 unrelated Finnish samples. Genes are depicted as annotated in GENCODE v19, potential regulatory regions identified by digital genomic footprinting in human cardiac myocytes (ENCODE: Encyclopedia of DNA Elements, experiment number ENCSR000ENJ) indicated by the black graph. SNPs independently associated with midregional proatrial natriuretic peptide (MR-proANP) colored with blue, SNPs independently associated with B-type natriuretic peptide (BNP) or NT-proBNP (N-terminal pro-B-type natriuretic peptide) colored in pink. SNPs associated with MR-proANP or NT-proBNP in previous studies colored in black. The NPPA-NPPB locus contained 3 initial association signals for BNP, NT-proBNP, and BNP:NT-proBNP ratio, depending on which of the phenotypes was tested (Table 1; Table III in the Data Supplement). Rs198379, situated 2055 base pairs downstream from the last exon of NPPB, associated with BNP (P=4.42×10−52). For NT-proBNP and BNP:NT-proBNP ratio, rs61761991 was the most statistically significant SNP (P=8.76×10−68 and P=4.81×10−103, respectively), and least absolute shrinkage and selection operator regression detected rs12406089 as a secondary signal for NT-proBNP (P=8.31×10−48). However, neither of these 2 SNPs associated with BNP when rs198379 was included in the model. The NPPA-NPPB locus, therefore, harbored only 1 variant, rs198379, independently associated with both BNP and NT-proBNP, with every C allele increasing BNP concentration by ≈4.5 pg/mL and NT-proBNP concentration by 9.6 pg/mL. Three SNPs associated independently with MR-proANP at the NPPA-NPPB locus (Table 1). The most statistically significant was rs3753584 (P=3.85×10−13), but the effect sizes of the 3 SNPs were broadly similar. Each allele of the SNPs correlated with a 2.5 to 5.0 pmol/L difference in MR-proANP concentration. The SNPs are found ≈40 kb downstream from NPPA within an area bound by regulatory proteins in human cardiomyocytes (ENCODE: Encyclopedia of DNA Elements, https://www.encodeproject.org, experiment ENCSR000ENJ).37 Because obesity disturbs the association of MR-proANP with BP, we studied the effect of body mass on the SNP associations by introducing body mass*SNP interaction terms to the regression models.15 The interaction terms were statistically nonsignificant (P>0.05) for both BMI as a continuous variable and obesity (BMI>30) as a categorical variable. Furthermore, because NPPA and NPPB are separated by <10 kb, any variant in this region might affect either both genes or only 1 of the 2. We explored this by fitting models containing all of the previously mentioned SNPs of the NPPA-NPPB locus and found that SNPs associated with MR-proANP did not associate with BNP or NT-proBNP and vice versa (Table V in the Data Supplement), indicating that their effects were specific to either MR-proANP or NT-proBNP (and BNP). Gene expression profiling in human cardiac tissue samples confirmed that the associations of SNPs with BNP or MR-proANP concentration and with NPPB and NPPA gene expression were positively correlated. In left ventricular tissue samples, SNPs associated with circulating MR-proANP concentration also associated with NPPA expression level (Figure III in the Data Supplement; Spearman rank correlation of P values, P=0.014), and SNPs associated with BNP concentration also associated with NPPB expression (P=0.004). Furthermore, the effect estimates for circulating MR-proANP concentration and NPPA expression in the left ventricle correlated (Pearson r=0.611; P=0.024) as did those for BNP and NPPB (r=0.735; P=0.001). A somewhat attenuated trend was also present in the atrial appendage samples, where the correlations between the effect estimates were statistically significant (MR-proANP versus NPPA r=0.508; P=0.021 and BNP versus NPPB r=0.481; P=0.033), but the correlations between association P values were not. In addition to NPPA and NPPB, the correlations were also significant for EXOSC10 and ENSG00000272482 (with MR-proANP) and EXOSC10 and MTHFR (with BNP). The regulatory effects underlying the MR-proANP and BNP associations near NPPA-NPPB may, therefore, be stronger in the left ventricle compared with the atrium and also selectively affect the expression of other nearby genes. PPP3CC and GALNT4 on Chromosomes 8 and 12 Rs7000551 on chromosome 8 near PPP3CC associated with BNP:NT-proBNP ratio (P=2.27×10−9). This correlation was driven by an effect on the NT-proBNP concentration as rs7000551 associated with NT-proBNP (P=3.72×10−5) but not with BNP (P=0.87) in the discovery GWAS sample. However, only the association of rs7000551 with BNP:NT-proBNP ratio met genome-wide significance and replicated. The genotype-specific mean BNP:NT-proBNP ratios for rs7000551 (AA=0.386; AG=0.412; GG=0.463) suggest an additive or multiplicative genetic effect with each G-allele raising the atio by ≈0.04 U or 10%. The association peak on chromosome 8 extends from the 3′ end of SLC39A14 into the promoter region and 5′ end of PPP3CC, with rs7000551 itself located in an intron of SLC39A14 (Figure 3). When we studied the coassociation of SNPs with BNP:NT-proBNP ratio and gene expression, PPP3CC and 2 antisense RNA genes ENSG00000245025 and ENSG00000248738 matched the prespecified criteria. SNPs associated with increased BNP:NT-proBNP ratio also associated with reduced expression of PPP3CC in both atrial and left ventricular tissue samples (atrial appendage, Pearson r=−0.70; P=0.006 and left ventricle, r=−0.81; P=0.021; Figure III in the Data Supplement). The coassociation with the 2 RNA genes was significant only in the left ventricular tissue samples. Both the physical location near the promoter of PPP3CC and the coassociation with its expression, therefore, suggest that the BNP:NT-proBNP ratio–associated SNPs tag a regulatory variant that alters the expression of PPP3CC in the heart. Figure 3. B-type natriuretic peptide (BNP):NT-proBNP (N-terminal pro-B-type natriuretic peptide) associated SNPs on chromosomes 8 and 12. Association of SNPs with BNP:NT-proBNP ratio on chromosomes 8 and 12 after meta-analyzing the results from the genome-wide association study (GWAS) and replication samples with lead SNPs from the GWAS indicated with purple diamonds. SNPs near POC1B and GALNT4 on chromosome 12 associated with NT-proBNP and BNP:NT-proBNP ratio. Rs11105298 and rs61378614, located in different introns of the POC1B gene (Figure 3), independently associated with the ratio (P=1.52×10−26 and P=3.98×10−9, respectively). The genes’ expression on chromosome 12 did not show a clear coassociation with BNP:NT-proBNP ratio because none of them was significant for all 3 predefined criteria. Association With BP Having identified the set of SNPs associated with the natriuretic peptide traits in the genome-wide significant loci, we next studied their correlation with systolic BP, diastolic BP, and hypertension in an independent sample. We fitted all SNPs simultaneously in each locus, excluding rs61761991 and rs12406089 on chromosome 1, which did not independently associate with BNP. After genotyping quality control, the study sample contained 27 059 participants with both BP measurements and SNP genotypes available. The 3 SNPs associated with MR-proANP also associated weakly with BP (Table 2; Figure IV in the Data Supplement). The point estimates of the MR-proANP increasing alleles’ effects were ≈0.25 mm Hg (diastolic BP) and 0.50 mm Hg (systolic BP). Only 1 of these SNPs was independently associated with hypertension as a binary end point (rs3753584; P=6.8×10−4). To assess the combined effect of the genetic differences in MR-proANP concentration on BP, we formed an allele-counting score of the 3 SNPs. The score explained 2.36% of the variance in MR-proANP concentration and a unit increase in the score associated with a 9% decrease in the odds ratio for hypertension (odds ratio=0.91; SE=0.0283; P=8.2×10−4). Table 2. Independent Effects of Genetic Variants on Natriuretic Peptides and Blood Pressure In contrast to MR-proANP, none of the SNPs correlated with BNP, NT-proBNP, or BNP:NT-proBNP ratio associated with BP. Rs198379, associated with NPPB expression and circulating BNP levels, did not associate with systolic or diastolic BP when adjusted for the nearby MR-proANP–correlated SNPs. SNPs near PPP3CC and GALNT4, correlated with NT-proBNP and BNP:NT-proBNP ratio, similarly did not associate with BP or hypertension. Discussion We performed a GWAS of circulating MR-proANP, BNP, and NT-proBNP concentration and BNP:NT-proBNP concentration ratio in 4932 samples with replication in 1373 samples. We then studied the effect of the natriuretic peptide–associated loci on systolic BP, diastolic BP, and hypertension in 27 059 additional samples. We detected a novel locus for BNP:NT-proBNP ratio on chromosome 8 near PPP3CC and fine-mapped 2 published loci on chromosomes 1 and 12 for their association with ANP and BNP and BP. The entire genome-wide SNP data explained from 14% to 23% of the variation in the natriuretic peptide traits in our population-based sample. These estimates are similar to those, for example, BMI (14%) or systolic BP (24%) published elsewhere, showing that the natriuretic peptide traits considered here have an additive genetic component comparable to traditional cardiovascular risk factors.38 The present study is the first to assess the NPPA-NPPB locus with a dense SNP panel simultaneously for MR-proANP, BNP, and NT-proBNP, extending the results of previous investigations.16–21 We identified 3 statistically independent cis variants associated with MR-proANP, and 1 variant associated with BNP and NT-proBNP. Analysis of gene expression data suggests that the protein-level cis associations stem from effects on NPPA and NPPB gene expression, affecting both atrial and ventricular tissue. Furthermore, even if the 2 genes are separated by <10 000 bp, their transcriptional regulation is decoupled to the extent that the ANP-associated SNPs had no observable effect on BNP and vice versa. Each of these SNPs, however, correlates with both MR-proANP and BNP concentrations, if the analysis is not adjusted for the other SNPs. This is crucial for the interpretation of results from Mendelian randomization studies using SNPs in this locus as instruments, such as those performed in relation with type 2 diabetes mellitus.39 SNPs on chromosome 8 near SLC39A14 and PPP3CC associate with BNP:NT-proBNP ratio. SLC39A14 belongs to the same large family of solute carrier proteins as SLC39A8 in the previously detected NT-proBNP associated locus on chromosome 4, but it is difficult to assess whether this is only coincidental.19,40 SNPs associated with increased BNP:NT-proBNP ratio correlated with decreased expression of PPP3CC in both left ventricular and atrial tissue samples, whereas no such correlation was present for SLC39A14. PPP3CC codes for 1 of the 3 alternative catalytic subunits of calcineurin, a phosphatase with a wide range of functions including the regulation of cardiac hypertrophic signaling.41 Originally characterized as a testis-specific calcineurin subunit, PPP3CC has been later detected in multiple tissues.34 Because of its central role in spermatogenesis, drugs inhibiting PPP3CC have been suggested as a potential male contraceptive.42 Our results indicate that systemic inhibition of calcineurin containing the subunit coded by PPP3CC may have unintended cardiovascular side effects. Two independent SNPs near POC1B and GALNT4 associated with NT-proBNP levels and BNP:NT-proBNP ratio in our study. An association of a SNP with NT-proBNP in this locus has been previously reported in whites.19 Analysis of gene expression in cardiac tissue failed to highlight any of the nearby genes but, as previously noted, GALNT4 is an attractive candidate.19 It codes for an aminoacyltranferase that initiates O-linked glycosylation, and proBNP is known to be O-glycosylated.43,44 According to data presented here, the association near GALNT4 is specific to NT-proBNP, supporting the hypothesis that proBNP may be a target of GALNT4. Because BNP and NT-proBNP are produced as a single polypeptide, deviations in their circulating concentration ratio should reflect their differential secretion or removal, the processing of proBNP, or factors disturbing the detection of the peptides. The latter is probably the case with rs61761991, located within the region of the NT-proBNP prohormone (NP_002512.1:p.Arg72His) used as the antigen to prepare the assay’s primary antibody.45,46 The variant, which effectively blocked the signal of the NT-proBNP assay, is rare or absent in other populations but significantly enriched in Finns, where the frequency of the T allele is ≈3%.47 One in 20 Finns will, therefore, have a measured concentration of NT-proBNP, which is ≈50% lower than the corresponding C-terminal BNP value, potentially causing false rule-out of suspected heart failure. The associations of SNPs near GALNT4 with BNP:NT-proBNP ratio may also relate to the detection of NT-proBNP rather than changes in its concentration, if they are indeed linked to the possible glycosylation of proBNP by GALNT4. How PPP3CC may affect BNP:NT-proBNP ratio is unclear. We adjusted the analysis for the estimated glomerular filtration rate, but confounding by kidney function cannot be ruled out. Experimental data has pointed to either similar or different cardiovascular effects of ANP and BNP, depending on the experimental setting.6–8,48 The results of this study are in line with some of the previous studies that identified ANP rather than BNP as an important regulator of BP. Genetically determined increases in ANP concentration decreased systolic and diastolic BP, but a smaller genetic decrease in BNP did not. Obesity did not modify the associations of SNPs with MR-proANP, showing that the transcriptional regulation of NPPA is at least partially unaffected by the reported ANP-decreasing effect of high body mass.15,49 According to data presented here, earlier genetic associations of the NPPA-NPPB locus with BP were driven by ANP-associated variants and should not be taken as evidence of any BP lowering effect of BNP.16 We conclude that there are interesting differences between ANP and BNP in humans that are yet to be fully elucidated and that genetics provides unique insights into the effects of lifelong alterations of these hormones. Appendix From the National Institute for Health and Welfare, Helsinki, Finland (P.P.S., A.S.H., J.K., J.G.E., A.J., V.S., K.K., M.P.); Institute for Molecular Medicine Finland, Helsinki (P.P.S., A.S.H., T.T., K.K., M.P.); Diabetes and Obesity Research Program (K.K., M.P.) and Department of General Practice and Primary Health Care, Helsinki University Hospital (J.G.E.), University of Helsinki, Finland; The Research Centre of Applied and Preventive Cardiovascular Medicine (O.R.) and Department of Clinical Physiology, Turku University Hospital (O.R.), University of Turku, Finland; Department of Clinical Chemistry, Fimlab Laboratories and Finnish Cardiovascular Research Center Tampere, Faculty of Medicine and Life Sciences, University of Tampere, Finland (T.L.); Department of Clinical Physiology, Tampere University Hospital, Finland (M.K.); Department of Clinical Physiology, University of Tampere School of Medicine, Finland (M.K.); Institute for Computational Medicine, Center for Life Course Health Research, Faculty of Medicine (J.K.), Biocenter Oulu (J.K.), and Center for Life Course Health Research, Faculty of Medicine (M.M.), University of Oulu, Finland; Folkhälsan Research Center, Helsinki, Finland (J.G.E.); Department of General and Interventional Cardiology, University Heart Center Hamburg, Germany (S.B., T.Z.); German Center for Cardiovascular research, partner site Hamburg/Lübeck/Kiel, Hamburg, Germany (S.B., T.Z.); and Estonian Genome Center, University of Tartu, Estonia (M.P.). Sources of Funding This study was supported by Aarne Koskelo Foundation; the Academy of Finland grants 269517, 250207, 269517, 283045, 297338, 286284 (Dr Lehtimäki), 134309(Eye), 126925, 121584, 124282, 129378(Salve), 117787(Gendi), and 41071(Skidi); Biomedicum Helsinki Foundation; the Competitive State Research Financing of the Expert Responsibility area of Tampere, Turku; Kuopio University Hospital (grant X51001); the Diabetes Research Foundation of the Finnish Diabetes Association; Emil Aaltonen Foundation; the EU FP7 grants 313010 (BBMRI-LPC), 305280 (MIMOmics), and HZ2020 633589 (Ageing with Elegans); Finnish Cultural Foundation; Finnish Foundation for Cardiovascular Research; Ida Montin Foundation; Integrative Life Science Doctoral Program of the University of Helsinki; Juho Vainio Foundation; Paavo Nurmi Foundation; Signe and Ane Gyllenberg Foundation; the Social Insurance Institution of Finland, Tampere Tuberculosis Foundation; and Yrjö Jahnsson Foundation. Disclosures None. Supplementary Material * Drs Kristiansson and Perola contributed equally to this work. The Data Supplement is available at http://circgenetics.ahajournals.org/lookup/suppl/doi:10.1161/CIRCGENETICS.117.001713/-/DC1. CLINICAL PERSPECTIVE Atrial natriuretic peptide and B-type natriuretic peptide are unique hormones secreted by cardiomyocytes, often used in the diagnostics of heart failure. They bind to the same receptor, but unexpected differences in their effects have been reported in both human and animal models. We used genome-wide association analysis to study genetic variation affecting their circulating concentration, identifying 8 variants near the genes NPPA, NPPB, PPP3CC, and GALNT4. Subsequently, we investigated the correlation between the natriuretic peptide–associated genetic variants and blood pressure. Genetic variants lowering the concentration of midregional proatrial natriuretic peptide associated with higher blood pressure, but we did not observe a similar blood pressure correlation with genetic variants affecting B-type natriuretic peptide or NT-proBNP (N-terminal pro-B-type natriuretic peptide). The effect sizes of the midregional proatrial natriuretic peptide correlated genetic variants on blood pressure were small, from 0.25 to 0.50 mm Hg per allele. Their combined effect, however, associated with a 9% difference in the odds ratio for hypertension, contributing significantly to the burden of high blood pressure in the general population. ==== Refs References 1. Kerkelä R Ulvila J Magga J Natriuretic peptides in the regulation of cardiovascular physiology and metabolic events. J Am Heart Assoc 2015 4 e002423 doi: 10.1161/JAHA.115.002423. 26508744 2. Ponikowski P Voors AA Anker SD Bueno H Cleland JG Coats AJ Authors/Task Force Members 2016 ESC Guidelines for the diagnosis and treatment of acute and chronic heart failure: The Task Force for the diagnosis and treatment of acute and chronic heart failure of the European Society of Cardiology (ESC)Developed with the special contribution of the Heart Failure Association (HFA) of the ESC. Eur Heart J 2016 37 2129 2200 doi: 10.1093/eurheartj/ehw128. 27206819 3. Steinhelper ME Cochrane KL Field LJ Hypotension in transgenic mice expressing atrial natriuretic factor fusion genes. Hypertension 1990 16 301 307 2144261 4. Ogawa Y Itoh H Tamura N Suga S Yoshimasa T Uehira M Molecular cloning of the complementary DNA and gene that encode mouse brain natriuretic peptide and generation of transgenic mice that overexpress the brain natriuretic peptide gene. J Clin Invest 1994 93 1911 1921 doi: 10.1172/JCI117182. 8182124 5. John SW Krege JH Oliver PM Hagaman JR Hodgin JB Pang SC Genetic decreases in atrial natriuretic peptide and salt-sensitive hypertension. Science 1995 267 679 681 7839143 6. Tamura N Ogawa Y Chusho H Nakamura K Nakao K Suda M Cardiac fibrosis in mice lacking brain natriuretic peptide. Proc Natl Acad Sci USA 2000 97 4239 4244 doi: 10.1073/pnas.070371497. 10737768 7. John SW Veress AT Honrath U Chong CK Peng L Smithies O Blood pressure and fluid-electrolyte balance in mice with reduced or absent ANP. Am J Physiol 1996 271 1 pt 2 R109 R114 8760210 8. Holditch SJ Schreiber CA Nini R Tonne JM Peng KW Geurts A B-type natriuretic peptide deletion leads to progressive hypertension, associated organ damage, and reduced survival: novel model for human hypertension. Hypertension 2015 66 199 210 doi: 10.1161/HYPERTENSIONAHA.115.05610. 26063669 9. O’Connor CM Starling RC Hernandez AF Armstrong PW Dickstein K Hasselblad V Effect of nesiritide in patients with acute decompensated heart failure. N Engl J Med 2011 365 32 43 doi: 10.1056/NEJMoa1100171. 21732835 10. Crozier IG Nicholls MG Ikram H Espiner EA Gomez HJ Warner NJ Haemodynamic effects of atrial peptide infusion in heart failure. Lancet 1986 2 1242 1245 2878131 11. Jensen KT Carstens J Pedersen EB Effect of BNP on renal hemodynamics, tubular function and vasoactive hormones in humans. Am J Physiol 1998 274 1 pt 2 F63 F72 9458824 12. Hynynen M Kupari M Salmenperä M Tikkanen I Heinonen J Fyhrquist F Hemodynamic effects of alpha-human atrial natriuretic peptide in healthy volunteers. J Cardiovasc Pharmacol 1988 11 711 715 2457767 13. Pidgeon GB Richards AM Nicholls MG Espiner EA Yandle TG Frampton C Differing metabolism and bioactivity of atrial and brain natriuretic peptides in essential hypertension. Hypertension 1996 27 906 913 8613267 14. Suga S Nakao K Hosoda K Mukoyama M Ogawa Y Shirakami G Receptor selectivity of natriuretic peptide family, atrial natriuretic peptide, brain natriuretic peptide, and C-type natriuretic peptide. Endocrinology 1992 130 229 239 doi: 10.1210/endo.130.1.1309330. 1309330 15. Asferg CL Andersen UB Linneberg A Hedley PL Christiansen M Goetze JP Serum proatrial natriuretic peptide does not increase with higher systolic blood pressure in obese men. Heart 2017 103 154 158 doi: 10.1136/heartjnl-2016-309462. 27496822 16. Newton-Cheh C Larson MG Vasan RS Levy D Bloch KD Surti A Association of common variants in NPPA and NPPB with circulating natriuretic peptides and blood pressure. Nat Genet 2009 41 348 353 doi: 10.1038/ng.328. 19219041 17. Pereira NL Tosakulwong N Scott CG Jenkins GD Prodduturi N Chai Y Circulating atrial natriuretic peptide genetic association study identifies a novel gene cluster associated with reduced NT-proANP, increased stroke and higher diastolic blood pressure. BMC Pharmacol Toxicol 2015 16 A37 18. Musani SK Fox ER Kraja A Bidulescu A Lieb W Lin H Genome-wide association analysis of plasma B-type natriuretic peptide in blacks: the Jackson Heart Study. Circ Cardiovasc Genet 2015 8 122 130 doi: 10.1161/CIRCGENETICS.114.000900. 25561047 19. Johansson Å Eriksson N Lindholm D Varenhorst C James S Syvänen AC PLATO Investigators Genome-wide association and Mendelian randomization study of NT-proBNP in patients with acute coronary syndrome. Hum Mol Genet 2016 25 1447 1456 doi: 10.1093/hmg/ddw012. 26908625 20. Del Greco M F Pattaro C Luchner A Pichler I Winkler T Hicks AA Genome-wide association analysis and fine mapping of NT-proBNP level provide novel insight into the role of the MTHFR-CLCN6-NPPA-NPPB gene cluster. Hum Mol Genet 2011 20 1660 1671 doi: 10.1093/hmg/ddr035. 21273288 21. Folkersen L Fauman E Sabater-Lleal M Strawbridge RJ Frånberg M Sennblad B IMPROVE study group Mapping of 79 loci for 83 plasma protein biomarkers in cardiovascular disease. PLoS Genet 2017 13 e1006706 doi: 10.1371/journal.pgen.1006706. 28369058 22. Semenov AG Tamm NN Seferian KR Postnikov AB Karpova NS Serebryanaya DV Processing of pro-B-type natriuretic peptide: furin and corin as candidate convertases. Clin Chem 2010 56 1166 1176 doi: 10.1373/clinchem.2010.143883. 20489134 23. Vartiainen E Jousilahti P Alfthan G Sundvall J Pietinen P Puska P Cardiovascular risk factor changes in Finland, 1972-1997. Int J Epidemiol 2000 29 49 56 10750603 24. Vartiainen E Laatikainen T Peltonen M Juolevi A Männistö S Sundvall J Thirty-five-year trends in cardiovascular risk factors in Finland. Int J Epidemiol 2010 39 504 518 doi: 10.1093/ije/dyp330. 19959603 25. Eriksson JG Forsén T Tuomilehto J Winter PD Osmond C Barker DJ Catch-up growth in childhood and death from coronary heart disease: longitudinal study. BMJ 1999 318 427 431 9974455 26. Järvelin MR Sovio U King V Lauren L Xu B McCarthy MI Early life factors and blood pressure at age 31 years in the 1966 northern Finland birth cohort. Hypertension 2004 44 838 846 doi: 10.1161/01.HYP.0000148304.33869.ee. 15520301 27. Johansson JK Puukka PJ Niiranen TJ Varis J Peltonen M Salomaa V Health 2000 score–development and validation of a novel cardiovascular risk score. Ann Med 2016 48 403 409 27187608 28. Viikari J Rönnemaa T Seppänen A Marniemi J Porkka K Räsänen L Serum lipids and lipoproteins in children, adolescents and young adults in 1980-1986. Ann Med 1991 23 53 59 2036206 29. Blankenberg S Zeller T Saarela O Havulinna AS Kee F Tunstall-Pedoe H MORGAM Project Contribution of 30 biomarkers to 10-year cardiovascular risk estimation in 2 population cohorts: the MONICA, risk, genetics, archiving, and monograph (MORGAM) biomarker project. Circulation 2010 121 2388 2397 doi: 10.1161/CIRCULATIONAHA.109.901413. 20497981 30. Viechtbauer W Conducting meta-analyses in R with the metafor package. J Stat Softw 2010 36 1 48 31. Sabourin J Nobel AB Valdar W Fine-mapping additive and dominant SNP effects using group-LASSO and fractional resample model averaging. Genet Epidemiol 2015 39 77 88 doi: 10.1002/gepi.21869. 25417853 32. Purcell S Neale B Todd-Brown K Thomas L Ferreira MA Bender D PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 2007 81 559 575 doi: 10.1086/519795. 17701901 33. Yang J Lee SH Goddard ME Visscher PM GCTA: a tool for genome-wide complex trait analysis. Am J Hum Genet 2011 88 76 82 doi: 10.1016/j.ajhg.2010.11.011. 21167468 34. GTEx Consortium Human genomics. the genotype-tissue expression (GTEx) pilot analysis: multitissue gene regulation in humans. Science 2015 348 648 660 25954001 35. Nicolae DL Gamazon E Zhang W Duan S Dolan ME Cox NJ Trait-associated SNPs are more likely to be eQTLs: annotation to enhance discovery from GWAS. PLoS Genet 2010 6 e1000888 doi: 10.1371/journal.pgen.1000888. 20369019 36. Maurano MT Humbert R Rynes E Thurman RE Haugen E Wang H Systematic localization of common disease-associated variation in regulatory DNA. Science 2012 337 1190 1195 doi: 10.1126/science.1222794. 22955828 37. ENCODE Project Consortium An integrated encyclopedia of DNA elements in the human genome. Nature 2012 489 57 74 22955616 38. Vattikuti S Guo J Chow CC Heritability and genetic correlations explained by common SNPs for metabolic syndrome traits. PLoS Genet 2012 8 e1002637 doi: 10.1371/journal.pgen.1002637. 22479213 39. Pfister R Sharp S Luben R Welsh P Barroso I Salomaa V Mendelian randomization study of B-type natriuretic peptide and type 2 diabetes: evidence of causal association from population studies. PLoS Med 2011 8 e1001112 doi: 10.1371/journal.pmed.1001112. 22039354 40. He L Vasiliou K Nebert DW Analysis and update of the human solute carrier (SLC) gene superfamily. Hum Genom 2009 3 1 41. Molkentin JD Lu JR Antos CL Markham B Richardson J Robbins J A calcineurin-dependent transcriptional pathway for cardiac hypertrophy. Cell 1998 93 215 228 9568714 42. Miyata H Satouh Y Mashiko D Muto M Nozawa K Shiba K Sperm calcineurin inhibition prevents mouse fertility with implications for male contraceptive. Science 2015 350 442 445 doi: 10.1126/science.aad0836. 26429887 43. Bennett EP Hassan H Mandel U Mirgorodskaya E Roepstorff P Burchell J Cloning of a human UDP-N-acetyl-alpha-D-galactosamine:polypeptide N-acetylgalactosaminyltransferase that complements other GalNAc-transferases in complete O-glycosylation of the MUC1 tandem repeat. J Biol Chem 1998 273 30472 30481 9804815 44. Schellenberger U O’Rear J Guzzetta A Jue RA Protter AA Pollitt NS The precursor to B-type natriuretic peptide is an O-linked glycoprotein. Arch Biochem Biophys 2006 451 160 166 doi: 10.1016/j.abb.2006.03.028. 16750161 45. Analytical characteristics of commercial MR-proANP, BNP and NT-proBNP assays as per the manufacturer http://www.ifcc.org/media/276711/IFCC%20NP%20Assay%20Table%20November%202014.pdf 46. Aken BL Ayling S Barrell D Clarke L Curwen V Fairley S The ensembl gene annotation system. Database (Oxford) 2016 2016 doi: 10.1093/database/baw093. 47. Lek M Karczewski K Minikel E Samocha K Banks E Fennell T Analysis of protein-coding genetic variation in 60,706 humans. BioRxiv 2016 030338 48. Hunt PJ Espiner EA Nicholls MG Richards AM Yandle TG Differing biological effects of equimolar atrial and brain natriuretic peptide infusions in normal man. J Clin Endocrinol Metab 1996 81 3871 3876 doi: 10.1210/jcem.81.11.8923831. 8923831 49. Wang TJ Larson MG Levy D Benjamin EJ Leip EP Wilson PW Impact of obesity on plasma natriuretic peptide levels. Circulation 2004 109 594 600 doi: 10.1161/01.CIR.0000112582.16683.EA. 14769680