Genome-Wide Meta-Analysis of Blood Pressure Response to β1-Blockers: Results From ICAPS (International Consortium of Antihypertensive Pharmacogenomics Studies)

Leveraging large-scale biobank EHRs to enhance pharmacogenetics of cardiometabolic disease medications

Electronic health records (EHRs) coupled with large-scale biobanks offer great promises to unravel the genetic underpinnings of treatment efficacy. However, medication-induced biomarker trajectories stemming from such records remain poorly studied. Here, we extract clinical and medication prescription data from EHRs and conduct GWAS and rare variant burden tests in the UK Biobank (discovery) and the All of Us program (replication) on ten cardiometabolic drug response outcomes including lipid response to statins, HbA1c response to metformin and blood pressure response to antihypertensives (N = 740-26,669). Our findings at genome-wide significance level recover previously reported pharmacogenetic signals and also include novel associations for lipid response to statins (N = 26,669) near LDLR and ZNF800. Importantly, these associations are treatment-specific and not associated with biomarker progression in medication-naive individuals. Furthermore, we demonstrate that individuals with higher genetically determined low-density and total cholesterol baseline levels experience increased absolute, albeit lower relative biomarker reduction following statin treatment. In summary, we systematically investigated the common and rare pharmacogenetic contribution to cardiometabolic drug response phenotypes in over 50,000 UK Biobank and All of Us participants with EHR and identified clinically relevant genetic predictors for improved personalized treatment strategies.

Drug classes affecting intracranial aneurysm risk: Genetic correlation and Mendelian randomization

INTRODUCTION

There is no non-invasive treatment to prevent aneurysmal subarachnoid hemorrhage (ASAH) caused by intracranial aneurysm (IA) rupture. We aimed to identify drug classes that may affect liability to IA using a genetic approach.

PATIENTS AND METHODS

Using genome-wide association summary statistics we calculated genetic correlation between unruptured IA (N = 2140 cases), ASAH (N = 5140) or the combined group, and liability to drug usage from 23 drug classes (N up to 320,000) independent of the risk factor high blood pressure. Next, we evaluated the causality and therapeutic potential of correlated drug classes using three different Mendelian randomization frameworks.

RESULTS

Correlations with IA were found for antidepressants, paracetamol, acetylsalicylic acid, opioids, beta-blockers, and peptic ulcer and gastro-esophageal reflux disease drugs. MR showed no evidence that genetically predicted usage of these drug classes caused IA. Genetically predicted high responders to antidepressant drugs were at higher risk of IA (odds ratio [OR] = 1.61, 95% confidence interval (CI) = 1.09-2.39, p = 0.018) and ASAH (OR = 1.68, 95% CI = 1.07-2.65, p = 0.024) if they used antidepressant drugs. This effect was absent in non-users. For beta-blockers, additional analyses showed that this effect was not independent of blood pressure after all. A complex and likely pleiotropic relationship was found between genetic liability to chronic multisite pain, pain medication usage (paracetamol, acetylsalicylic acid, and opioids), and IA.

CONCLUSIONS

We did not find drugs decreasing liability to IA and ASAH but found that antidepressant drugs may increase liability. We observed pleiotropic relationships between IA and other drug classes and indications. Our results improve understanding of pathogenic mechanisms underlying IA.

Genetic polymorphisms in ADRB1, ADRB2 and CYP2D6 genes and response to beta-blockers in patients with acute coronary syndrome

Betablockers (BBs) are prescribed for ischaemia in patients with acute coronary syndrome (ACS). In Spain, bisoprolol and carvedilol are the most prescribed BBs, but patients often had to discontinue them due to adverse effects. Single nucleotide polymorphisms (SNPs) in ADRB1, ADRB2 and CYP2D6 genes have strong evidence of pharmacogenetic association with BBs in heart failure or hypertension, but the evidence in ACS is limited. Therefore, our study focuses on investigating how these genes influence the response to BBs in ACS patients. We analysed the association between SNPs in ADRB1 Gly389Arg (rs1801253) and Ser49Gly (rs1801252), ADRB2 Gly16Arg (rs1042713) and Glu27Gln (rs1042714), and CYP2D* 6 (*2- rs1080985, *4- rs3892097, *10 - rs1065852) and the occurrence of bradycardia/hypotension events during one year of follow-up. We performed an observational study and included 285 ACS-PCI-stent patients. A first analysis including patients treated with bisoprolol and a second analysis including patients treated with other BBs were performed. We found that the presence of the G allele (Glu) of the ADRB2 gene (rs1042714; Glu27Gln) conferred a protective effect against hypotension-induced by BBs; OR (CI 95%) = 0,14 (0,03-0,60), p < 0.01. The ADRB2 (rs1042713; Gly16Arg) GG genotype could also prevent hypotensive events; OR (CI 95%) = 0.49 (0.28-0.88), p = 0015. SNPs in ADRB1 and CYP2D6 * 2, CYP2D6 * 4 weren´t associated with primary events. The effect of CYP2D6 * 10 does not seem to be relevant for the response to BBs. According to our findings, SNPs in ADRB2 (rs1042713, rs1042714) could potentially affect the response and tolerance to BBs in ACS-patients. Further studies are necessary to clarify the impact of ADRB2 polymorphisms.

Genome-Wide Association Study of Beta-Blocker Survival Benefit in Black and White Patients with Heart Failure with Reduced Ejection Fraction

In patients with heart failure with reduced ejection fraction (HFrEF), individual responses to beta-blockers vary. Candidate gene pharmacogenetic studies yielded significant but inconsistent results, and they may have missed important associations. Our objective was to use an unbiased genome-wide association study (GWAS) to identify loci influencing beta-blocker survival benefit in HFrEF patients. Genetic variant × beta-blocker exposure interactions were tested in Cox proportional hazards models for all-cause mortality stratified by self-identified race. The models were adjusted for clinical risk factors and propensity scores. A prospective HFrEF registry (469 black and 459 white patients) was used for discovery, and linkage disequilibrium (LD) clumped variants with a beta-blocker interaction of p < 5 × 10-5, were tested for Bonferroni-corrected validation in a multicenter HFrEF clinical trial (288 black and 579 white patients). A total of 229 and 18 variants in black and white HFrEF patients, respectively, had interactions with beta-blocker exposure at p < 5 × 10-5 upon discovery. After LD-clumping, 100 variants and 4 variants in the black and white patients, respectively, remained for validation but none reached statistical significance. In conclusion, genetic variants of potential interest were identified in a discovery-based GWAS of beta-blocker survival benefit in HFrEF patients, but none were validated in an independent dataset. Larger cohorts or alternative approaches, such as polygenic scores, are needed.

miR-4432 Targets FGFBP1 in Human Endothelial Cells

MicroRNAs (miRs) are small non-coding RNAs that modulate the expression of several target genes. Fibroblast growth factor binding protein 1 (FGFBP1) has been associated with endothelial dysfunction at the level of the blood-brain barrier (BBB). However, the underlying mechanisms are mostly unknown and there are no studies investigating the relationship between miRs and FGFBP1. Thus, the overarching aim of the present study was to identify and validate which miR can specifically target FGFBP1 in human brain microvascular endothelial cells, which represent the best in vitro model of the BBB. We were able to identify and validate miR-4432 as a fundamental modulator of FGFBP1 and we demonstrated that miR-4432 significantly reduces mitochondrial oxidative stress, a well-established pathophysiological hallmark of hypertension.

Genetically predicted on-statin LDL response is associated with higher intracerebral haemorrhage risk

Statins lower low-density lipoprotein cholesterol and are widely used for the prevention of atherosclerotic cardiovascular disease. Whether statin-induced low-density lipoprotein reduction increases risk of intracerebral haemorrhage has been debated for almost two decades. Here, we explored whether genetically predicted on-statin low-density lipoprotein response is associated with intracerebral haemorrhage risk using Mendelian randomization. Using genomic data from randomized trials, we derived a polygenic score from 35 single nucleotide polymorphisms of on-statin low-density lipoprotein response and tested it in the population-based UK Biobank. We extracted statin drug and dose information from primary care data on a subset of 225 195 UK Biobank participants covering a period of 29 years. We validated the effects of the genetic score on longitudinal low-density lipoprotein measurements with generalized mixed models and explored associations with incident intracerebral haemorrhage using Cox regression analysis. Statins were prescribed at least once to 75 973 (31%) of the study participants (mean 57 years, 55% females). Among statin users, mean low-density lipoprotein decreased by 3.45 mg/dl per year [95% confidence interval (CI): (-3.47, -3.42)] over follow-up. A higher genetic score of statin response [1 standard deviation (SD) increment] was associated with significant additional reductions in low-density lipoprotein levels [-0.05 mg/dl per year, (-0.07, -0.02)], showed concordant lipidomic effects on other lipid traits as statin use and was associated with a lower risk for incident myocardial infarction [hazard ratio per SD increment 0.98 95% CI (0.96, 0.99)] and peripheral artery disease [hazard ratio per SD increment 0.93 95% CI (0.87, 0.99)]. Over a 11-year follow-up period, a higher genetically predicted statin response among statin users was associated with higher intracerebral haemorrhage risk in a model adjusting for statin dose [hazard ratio per SD increment 1.16, 95% CI (1.05, 1.28)]. On the contrary, there was no association with intracerebral haemorrhage risk among statin non-users (P = 0.89). These results provide further support for the hypothesis that statin-induced low-density lipoprotein reduction may be causally associated with intracerebral haemorrhage risk. While the net benefit of statins for preventing vascular disease is well-established, these results provide insights about the personalized response to statin intake and the role of pharmacological low-density lipoprotein lowering in the pathogenesis of intracerebral haemorrhage.

Chromosomal Region 11p14.1 is Associated with Pharmacokinetics and Pharmacodynamics of Bisoprolol

Purpose

Bisoprolol is a widely used beta-blocker in patients with cardiovascular diseases. As with other beta-blockers, there is variability in response to bisoprolol, but the underlying reasons for this have not been clearly elucidated. Our aim was to investigate genetic factors that affect bisoprolol pharmacokinetics (PK) and pharmacodynamics (PD), and potentially the clinical outcomes.

Patients and Methods

Patients with non-ST elevation acute coronary syndrome were recruited prospectively on admission to hospital and followed up for up to 2 years. Patients from this cohort who were on treatment with bisoprolol, at any dose, had bisoprolol adherence data and a plasma sample, one month after discharge from index hospitalisation were included in the study. Individual bisoprolol clearance values were estimated using population pharmacokinetic modeling. Genome-wide association analysis after genotyping was undertaken using an Illumina HumanOmniExpressExome-8 v1.0 BeadChip array, while CYP2D6 copy number variations were determined by PCR techniques and phenotypes for CYP2D6 and CYP3A were inferred from the genotype. GWAS significant SNPs were analysed for heart rate response to bisoprolol in an independent cohort of hypertensive subjects.

Results

Six hundred twenty-two patients on bisoprolol underwent both PK and genome wide analysis. The mean (IQR) of the estimated clearance in this population was 13.6 (10.0-18.0) L/h. Bisoprolol clearance was associated with rs11029955 (p=7.17×10-9) mapped to the region of coiled-coil domain containing 34 region (CCDC34) on chromosome 11, and with rs116702638 (p=2.54×10-8). Each copy of the minor allele of rs11029955 was associated with 2.2 L/h increase in clearance. In an independent cohort of hypertensive subjects, rs11029955 was associated with 24-hour heart rate response to 4-week treatment with bisoprolol (p= 9.3×10-5), but not with rs116702638.

Conclusion

A novel locus on the chromosomal region 11p14.1 was associated with bisoprolol clearance in a real-world cohort of patients and was validated in independent cohort with a pharmacodynamic association.

Pharmacoepigenetics of hypertension: genome-wide methylation analysis of responsiveness to four classes of antihypertensive drugs using a double-blind crossover study design

Essential hypertension remains the leading risk factor of global disease burden, but its treatment goals are often not met. We investigated whether DNA methylation is associated with antihypertensive responses to a diuretic, a beta-blocker, a calcium channel blocker or an angiotensin receptor antagonist. In addition, since we previously showed an SNP at the transcription start site (TSS) of the catecholamine biosynthesis-related ACY3 gene to associate with blood pressure (BP) response to beta-blockers, we specifically analysed the association of methylation sites close to the ACY3 TSS with BP responses to beta-blockers. We conducted an epigenome-wide association study between leukocyte DNA methylation and BP responses to antihypertensive monotherapies in two hypertensive Finnish cohorts: the GENRES (https://clinicaltrials.gov/ct2/show/NCT03276598; amlodipine 5 mg, bisoprolol 5 mg, hydrochlorothiazide 25 mg, or losartan 50 mg daily) and the LIFE-Fin studies (https://clinicaltrials.gov/ct2/show/NCT00338260; atenolol 50 mg or losartan 50 mg daily). The monotherapy groups consisted of approximately 200 individuals each. We identified 64 methylation sites to suggestively associate (P < 1E-5) with either systolic or diastolic BP responses to a particular study drug in GENRES. These associations did not replicate in LIFE-Fin . Three methylation sites close to the ACY3 TSS were associated with systolic BP responses to bisoprolol in GENRES but not genome-wide significantly (P < 0.05). No robust associations between DNA methylation and BP responses to four different antihypertensive drugs were identified. However, the findings on the methylation sites close to the ACY3 TSS may support the role of ACY3 genetic and epigenetic variation in BP response to bisoprolol.

Pharmacogenomics Informs Cardiovascular Pharmacotherapy

Precision medicine exemplifies the emergence of personalized treatment options which may benefit specific patient populations based upon their genetic makeup. Application of pharmacogenomics requires an understanding of how genetic variations impact pharmacokinetic and pharmacodynamic properties. This particular approach in pharmacotherapy is helpful because it can assist in and improve clinical decisions. Application of pharmacogenomics to cardiovascular pharmacotherapy provides for the ability of the medical provider to gain critical knowledge on a patient's response to various treatment options and risk of side effects.

Adverse Cardiovascular Outcomes and Antihypertensive Treatment: A Genome-Wide Interaction Meta-Analysis in the International Consortium for Antihypertensive Pharmacogenomics Studies

We sought to identify genome-wide variants influencing antihypertensive drug response and adverse cardiovascular outcomes, utilizing data from four randomized controlled trials in the International Consortium for Antihypertensive Pharmacogenomics Studies (ICAPS). Genome-wide antihypertensive drug-single nucleotide polymorphism (SNP) interaction tests for four drug classes (β-blockers, n = 9,195; calcium channel blockers (CCBs), n = 10,511; thiazide/thiazide-like diuretics, n = 3,516; ACE-inhibitors/ARBs, n = 2,559) and cardiovascular outcomes (incident myocardial infarction, stroke, or death) were analyzed among patients with hypertension of European ancestry. Top SNPs from the meta-analyses were tested for replication of cardiovascular outcomes in an independent Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) study (n = 21,267), blood pressure (BP) response in independent ICAPS studies (n = 1,552), and ethnic validation in African Americans from the Genetics of Hypertension Associated Treatment study (GenHAT; n = 5,115). One signal reached genome-wide significance in the β-blocker-SNP interaction analysis (rs139945292, Interaction P = 1.56 × 10^-8 ). rs139945292 was validated through BP response to β-blockers, with the T-allele associated with less BP reduction (systolic BP response P = 6 × 10^-4 , Beta = 3.09, diastolic BP response P = 5 × 10^-3 , Beta = 1.53). The T-allele was also associated with increased adverse cardiovascular risk within the β-blocker treated patients' subgroup (P = 2.35 × 10^-4 , odds ratio = 1.57, 95% confidence interval = 1.23-1.99). The locus showed nominal replication in CHARGE, and consistent directional trends in β-blocker treated African Americans. rs139945292 is an expression quantitative trait locus for the 50 kb upstream gene NTM (neurotrimin). No SNPs attained genome-wide significance for any other drugs classes. Top SNPs were located near CALB1 (CCB), FLJ367777 (ACE-inhibitor), and CES5AP1 (thiazide). The NTM region is associated with increased risk for adverse cardiovascular outcomes and less BP reduction in β-blocker treated patients. Further investigation into this region is warranted.

Pharmacogenomic Effects of β-Blocker Use on Femoral Neck Bone Mineral Density

CONTEXT

Recent studies have shown that β-blocker (BB) users have a decreased risk of fracture and higher bone mineral density (BMD) compared to nonusers, likely due to the suppression of adrenergic signaling in osteoblasts, leading to increased BMD. There is also variability in the effect size of BB use on BMD in humans, which may be due to pharmacogenomic effects.

OBJECTIVE

To investigate potential single-nucleotide variations (SNVs) associated with the effect of BB use on femoral neck BMD, we performed a cross-sectional analysis using clinical data, dual-energy x-ray absorptiometry, and genetic data from the Framingham Heart Study's (FHS) Offspring Cohort. We then sought to validate our top 4 genetic findings using data from the Rotterdam Study, the BPROOF Study, the Malta Osteoporosis Fracture Study (MOFS), and the Hertfordshire Cohort Study.

METHODS

We used sex-stratified linear mixed models to determine SNVs that had a significant interaction effect with BB use on femoral neck (FN) BMD across 11 gene regions. We also evaluated the association of our top SNVs from the FHS with microRNA (miRNA) expression in blood and identified potential miRNA-mediated mechanisms by which these SNVs may affect FN BMD.

RESULTS

One variation (rs11124190 in HDAC4) was validated in females using data from the Rotterdam Study, while another (rs12414657 in ADRB1) was validated in females using data from the MOFS. We performed an exploratory meta-analysis of all 5 studies for these variations, which further validated our findings.

CONCLUSION

This analysis provides a starting point for investigating the pharmacogenomic effects of BB use on BMD measures.

Pharmacogenomics in Cardiovascular Diseases

Cardiovascular pharmacogenomics is the study and identification of genomic markers that are associated with variability in cardiovascular drug response, cardiovascular drug-related outcomes, or cardiovascular drug-related adverse events. This overview presents an introduction and historical background to cardiovascular pharmacogenomics, and a protocol for designing a cardiovascular pharmacogenomics study. Important considerations are also included for constructing a cardiovascular pharmacogenomics phenotype, designing the replication or validation strategy, common statistical approaches, and how to put the results in context with the cardiovascular drug or cardiovascular disease under investigation. © 2021 Wiley Periodicals LLC. Basic Protocol: Designing a cardiovascular pharmacogenomics study.

Therapeutic Properties of Highly Selective β-blockers With or Without Additional Vasodilator Properties: Focus on Bisoprolol and Nebivolol in Patients With Cardiovascular Disease

Bisoprolol and nebivolol are highly selective β₁-adrenoceptor antagonists, with clinical indications in many countries within the management of heart failure with reduced left ventricular ejection fraction (HFrEF), ischaemic heart disease (IHD), and hypertension. Nebivolol has additional vasodilator actions, related to enhanced release of NO in the vascular wall. In principle, this additional mechanism compared with bisoprolol might lead to more potent vasodilatation, which in turn might influence the effectiveness of nebivolol in the management of HFrEF, IHD and hypertension. In this article, we review the therapeutic properties of bisoprolol and nebivolol, as representatives of “second generation” and “third generation” β-blockers, respectively. Although head-to-head trials are largely lacking, there is no clear indication from published studies of an additional effect of nebivolol on clinical outcomes in patients with HFrEF or the magnitude of reductions of BP in patients with hypertension.

A Genome-Wide Association Study for Hypertensive Kidney Disease in Korean Men

Hypertension is one of the major risk factors for chronic kidney disease (CKD), and the coexistence of hypertension and CKD increases morbidity and mortality. Although many genetic factors have been identified separately for hypertension and kidney disease, studies specifically focused on hypertensive kidney disease (HKD) have been rare. Therefore, this study aimed to identify loci or genes associated with HKD. A genome-wide association study (GWAS) was conducted using two Korean cohorts, the Health Examinee (HEXA) and Korean Association REsource (KARE). Consequently, 19 single nucleotide polymorphisms (SNPs) were found to be significantly associated with HKD in the discovery and replication phases (p < 5 × 10⁻⁸, p < 0.05, respectively). We further analyzed HKD-related traits such as the estimated glomerular filtration rate (eGFR), creatinine, blood urea nitrogen (BUN), systolic blood pressure (SBP) and diastolic blood pressure (DBP) at the 14q21.2 locus, which showed a strong linkage disequilibrium (LD). Expression quantitative trait loci (eQTL) analysis was also performed to determine whether HKD-related SNPs affect gene expression changes in glomerular and arterial tissues. The results suggested that the FANCM gene may affect the development of HKD through an integrated analysis of eQTL and GWAS and was the most significantly associated candidate gene. Taken together, this study indicated that the FANCM gene is involved in the pathogenesis of HKD. Additionally, our results will be useful in prioritizing other genes for further experiments.

Human essential hypertension: no significant association of polygenic risk scores with antihypertensive drug responses

Polygenic risk scores (PRSs) for essential hypertension, calculated from > 900 genomic loci, were recently found to explain a significant fraction of hypertension heritability and complications. To investigate whether variation of hypertension PRS also captures variation of antihypertensive drug responsiveness, we calculated two different PRSs for both systolic and diastolic blood pressure: one based on the top 793 independent hypertension-associated single nucleotide polymorphisms and another based on over 1 million genome-wide variants. Using our pharmacogenomic GENRES study comprising four different antihypertensive monotherapies (n ~ 200 for all drugs), we identified a weak, but (after Bonferroni correction) statistically nonsignificant association of higher genome-wide PRSs with weaker response to a diuretic. In addition, we noticed a correlation between high genome-wide PRS and electrocardiographic left ventricular hypertrophy. Finally, using data of the Finnish arm of the LIFE study (n = 346), we found that PRSs for systolic blood pressure were slightly higher in patients with drug-resistant hypertension than in those with drug-controlled hypertension (p = 0.03, not significant after Bonferroni correction). In conclusion, our results indicate that patients with elevated hypertension PRSs may be predisposed to difficult-to-control hypertension and complications thereof. No general association between a high PRS and less efficient drug responsiveness was noticed.

Association of West African ancestry and blood pressure control among African Americans taking antihypertensive medication in the Jackson Heart Study

African Americans have a wide range of continental genetic ancestry. It is unclear whether racial differences in blood pressure (BP) control are related to ancestral background. The authors analyzed data from the Jackson Heart Study, a cohort exclusively comprised of self-identified African Americans, to assess the association between estimated West African ancestry (WAA) and BP control (systolic and diastolic BP < 140/90 mm Hg). Three nested modified Poisson regression models were used to calculate prevalence ratios for BP control associated with the three upper quartiles, separately, vs the lowest quartile of West African ancestry. The authors analyzed data from 1658 participants with hypertension who reported taking all of their antihypertensive medications in the previous 24 hours. WAA was estimated using 389 ancestry informative markers and categorized into quartiles (Q1: <73.7%, Q2: >73.7%-81.0%, Q3: >81.0%-86.3%, and Q4: >86.3%). The proportion of participants with controlled BP in the lowest-to-highest WAA quartile was 75.2%, 76.1%, 76.6%, and 74.4%. The prevalence ratios (95% CI) for controlled BP comparing Q2, Q3, and Q4 to Q1 of WAA were 1.00 (0.93-1.08), 1.02 (0.94-1.10), and 0.99 (0.91-1.07), respectively. Among African Americans in the Jackson Heart Study taking antihypertensive medication, BP control rates did not differ across quartiles of WAA.

Pharmacogenomics And Hypertension: Current Insights

Hypertension is a multifactorial disease that affects approximately one billion subjects worldwide and is a major risk factor associated with cardiovascular events, including coronary heart disease and cerebrovascular accidents. Therefore, adequate blood pressure control is important to prevent these events, reducing premature mortality and disability. However, only one third of patients have the effective control of blood pressure, despite several classes of antihypertensive drugs available. These disappointing outcomes may be at least in part explained by interpatient variability in drug response due to genetic polymorphisms. To address the effects of genetic polymorphisms on blood pressure responses to the antihypertensive drug classes, studies have applied candidate genes and genome wide approaches. More recently, a third approach that considers gene-gene interactions has also been applied in hypertension pharmacogenomics. In this article, we carried out a comprehensive review of recent findings on the pharmacogenomics of antihypertensive drugs, including diuretics, β-blockers, angiotensin-converting enzyme inhibitors and angiotensin II receptor blockers, and calcium channel blockers. We also discuss the limitations and inconsistences that have been found in hypertension pharmacogenomics and the challenges to implement this valuable approach in clinical practice.

Genome-Wide Meta-Analysis of Blood Pressure Response to β1-Blockers: Results From ICAPS (International Consortium of Antihypertensive Pharmacogenomics Studies)

BackgroundThere exists a wide interindividual variability in blood pressure (BP) response to β₁-blockers. To identify the genetic determinants of this variability, we performed a pharmacogenomic genome-wide meta-analysis of genetic variants influencing β₁-blocker BP response.Methods and ResultsGenome-wide association analysis for systolic BP and diastolic BP response to β₁-blockers from 5 randomized clinical trials consisting of 1254 patients with hypertension of European ancestry were combined in meta-analysis and single nucleotide polymorphisms (SNPs) with P<10^-4 were tested for replication in 2 independent randomized clinical trials of β₁-blocker-treated patients of European ancestry (n=1552). Regions harboring the replicated SNPs were validated in a β₁-blocker-treated black cohort from 2 randomized clinical trials (n=315). A missense SNP rs28404156 in BST1 was associated with systolic BP response to β₁-blockers in the discovery meta-analysis (P=9.33×10^-5, β=-3.21 mm Hg) and replicated at Bonferroni significance (P=1.85×10^-4, β=-4.86 mm Hg) in the replication meta-analysis with combined meta-analysis approaching genome-wide significance (P=2.18×10^-7). This SNP in BST1 is in linkage disequilibrium with several SNPs with putative regulatory functions in nearby genes, including CD38, FBXL5, and FGFBP1, all of which have been implicated in BP regulation. SNPs in this genetic region were also associated with BP response in the black cohort.ConclusionsData from randomized clinical trials of 8 European ancestry and 2 black cohorts support the assumption that BST1 containing locus on chromosome 4 is associated with β₁-blocker BP response. Given the previous associations of this region with BP, this is a strong candidate region for future functional studies and potential use in precision medicine approaches for BP management and risk prediction.