Johnson, Julie A.

Introduction: Hydrochlorothiazide (HCTZ) is among the most commonly prescribed antihypertensives in the US, yet <50% of HCTZ treated patients achieve blood pressure (BP) control. Hypothesis: Integrating metabolomic and genomic profiles of HCTZ treated patients could identify novel pathways/biomarkers for optimizing the use of HCTZ. Methods: Primary analysis included 228 white hypertensives treated with HCTZ from the Pharmacogenomic Evaluation of Antihypertensive Response (PEAR) study. Genotyping was determined via Illumina Omni 1M-Quad Chip and untargeted metabolomics was performed on baseline fasting plasma samples using a GC-TOF MS platform. Pathway analysis was used to integrate metabolites significantly associated with HCTZ BP response (FDR<.05) and genetic signals at p<5x10-5 from HCTZ BP genome-wide analysis. Replication was conducted in another 212 white HCTZ treated patients. Replicated genetic signals were used to create a BP response score by summing their BP lowering alleles. We validated this score by testing it in 196 white HCTZ treated patients in the Genetic Epidemiology of Responses to Antihypertensives (GERA) study. Results: Thirteen metabolites were significantly associated with HCTZ systolic BP (SBP) and diastolic BP (DBP) responses (FDR<.05). Metabolomics-genomics integration revealed significant metabolites along with rs2727563 PRKAG2, rs12604940 DCC, and rs13262930 EPHX2, in the netrin signalling pathway (p=1x10-5), as potential markers significantly influencing HCTZ BP response. We successfully replicated those 3 genetic signals and added their alleles to create a response score which explained 11.3% and 11.9% of the variability to HCTZ SBP and DBP responses, respectively. We found that patients carrying one “response” allele had a significantly worse response than carriers of six alleles (ΔSBP/ΔDBP: -1.5/1.2 vs -16.3/-10.4 mmHg, respectively, SBP score p=1x10-8 and DBP score p=3x10-9). We further validated this response score by testing it in the GERA (DBP score p=0.03, SBP score p=0.07). Conclusion: This study highlights the power of using different omics to identify novel pathways/biomarkers of drug response, and suggests that PRKAG2, DCC and EPHX2 might be important determinants of HCTZ BP response.

Introduction: β-blockers, including atenolol, are among the most widely prescribed antihypertensive drug classes. However, their use has been associated with hyperglycemia and incident diabetes, which poses concern for prescribers and patients. The Pharmacogenomic Evaluation of Antihypertensive Responses (PEAR) study is a multicenter clinical trial with one of its objectives to identify genetic influences on adverse metabolic events in patients treated with atenolol. Here we report a genomics-metabolomics integrative approach to identify biomarkers associated with atenolol-induced glucose change. Methods: PEAR participants were genotyped using the Illumina Omni 1M-Quad Chip. Metabolite profiles were acquired through GC-TOF-MS. Association between baseline levels of metabolites and glucose change was conducted using a multivariate linear regression model adjusted for age, gender and baseline fasting glucose. Genomics-metabolomics integration included selecting SNPs within genes involved in the synthesis and degradation of top metabolites and testing their effect on glucose change in 234 atenolol-treated Caucasians in PEAR, adjusting for age, gender, baseline glucose and principal components 1 and 2. A false discovery rate (FDR) of .05 was used to adjust for multiple comparisons. Replication of significant SNPs was conducted in an additional 228 participants treated with atenolol in PEAR. Results: Metabolomics analysis revealed 3 metabolites associated with glucose change that passed FDR adjustment. The top metabolite signal was beta-alanine. Increasing levels of baseline beta-alanine was significantly correlated with greater glucose increases after atenolol treatment (p=.0003). A total of 79 SNPs in 12 genes in the beta-alanine pathway were explored. SNP rs2669429 in DPYS was significantly associated with glucose change, after FDR adjustment. We found that G allele carriers had significantly greater glucose change compared with non-carriers (p=.0006; β=+2.76 mg/dL). This finding was successfully replicated in the other 228 patients who received atenolol (p=.04; β=+1.86 mg/dL). Conclusions: These results suggest that beta alanine and DPYS rs2669429 SNP might be an important predictor of atenolol-induced dysglycemia.

Introduction: Thiazide diuretics (TZDs) are the second most frequently prescribed class of antihypertensives. But <50% patients achieve goal blood pressure (BP) on TZD monotherapy, mostly because therapy is chosen overlooking the heterogeneity in patient responses. Circulating microRNAs (miRNAs) remain attractive as non-invasive biomarkers to predict drug response and tailor therapy for optimal outcomes. Hypothesis: Plasma miRNAs may play a role in TZD response; we aim to identify those associated with hydrochlorothiazide (HCTZ) response. Methods: We profiled 754 miRNAs in baseline (untreated) plasma samples of 36 uncomplicated hypertensive adults (18 Responders (R) and 18 Non-responders (NR) based on their diastolic BP (DBP) response to HCTZ therapy) from the Pharmacogenomic Evaluation of Antihypertensive Responses clinical trial, using the real-time qPCR based TaqMan OpenArray Human microRNA panel. We filtered out miRNAs with low amplification scores and Cq>38 and samples with low miRNA expression. Combination of miR-223 and miR-19b (identified as most stably expressed miRNAs by NormFinder Software) was used to normalize expression across samples. We used MannWhitney U test to identify microRNAs differentially expressed between R and NR; linear and logistic regressions for associations of miRNA expression with BP response. FDA corrected p<0.05 was used as significance threshold, unadjusted p<0.05 as suggestive threshold. Results: Mean systolic BP (SBP)/DBP change in R was 14/10 mmHg and in NR 0/1 mmHg. We found 77 miRNAs to be consistently expressed across samples, of which 11 were upregulated and 14 downregulated with >2-fold difference between R and NR. Though no miRNAs reached FDR p<0.05, miR-376a (fold change(fc)= 10.1, p= 0.035) and miR-17 (fc= -0.45, p= 0.017) were differentially expressed between R and NR and their expression associated with DBP and SBP change, along with miR-10b, miR-24, miR-134 associated with SBP change at p<0.05. Conclusions: In this first ever genome-wide miRNA study of antihypertensive drug response, we discovered circulating miRNAs associated with BP response to HCTZ. Future studies are planned to validate these findings in a larger cohort, across different race group and across different TZD.

Abstract Acute renal depletion of sorting nexin 1 (SNX1) in mice results in blunted natriuretic response and hypertension due to impaired dopamine D 5 receptor (D 5 R) activity. We elucidated the molecular mechanisms for these phenotypes in Snx1 −/− mice. These mice had increased renal expressions of angiotensin II type 1 receptor (AT 1 R), NADPH oxidase (NOX) subunits, D 5 R, and NaCl cotransporter. Basal reactive oxygen species (ROS), NOX activity, and blood pressure (BP) were also higher in Snx1 ‐/‐ mice, which were normalized by apocynin, a drug that prevents NOX assembly. Renal proximal tubule (RPT) cells from hypertensive (HT) Euro‐American males had deficient SNX1 activity, impaired D 5 R endocytosis, and increased ROS compared with cells from normotensive (NT) Euro‐American males. siRNA‐mediated depletion of SNX1 in RPT cells from NT subjects led to a blunting of D 5 R agonist‐induced increase in cAMP production and decrease in Na + transport, effects that were normalized by over‐expression of SNX1. Among HT African‐Americans, three of the 12 single nucleotide polymorphisms interrogated for the SNX1 gene were associated with a decrease in systolic BP in response to hydrochlorothiazide (HCTZ). The results illustrate a new paradigm for the development of hypertension and imply that the trafficking protein SNX1 may be a crucial determinant for hypertension and response to antihypertensive therapy.

Background There exists a wide interindividual variability in blood pressure ( BP ) response to β 1 ‐blockers. To identify the genetic determinants of this variability, we performed a pharmacogenomic genome‐wide meta‐analysis of genetic variants influencing β 1 ‐blocker BP response. Methods and Results Genome‐wide association analysis for systolic BP and diastolic BP response to β 1 ‐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 β 1 ‐blocker–treated patients of European ancestry (n=1552). Regions harboring the replicated SNP s were validated in a β 1 ‐blocker–treated black cohort from 2 randomized clinical trials (n=315). A missense SNP rs28404156 in BST 1 was associated with systolic BP response to β 1 ‐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 BST 1 is in linkage disequilibrium with several SNP s with putative regulatory functions in nearby genes, including CD 38 , FBXL 5 , and FGFBP 1 , all of which have been implicated in BP regulation. SNP s in this genetic region were also associated with BP response in the black cohort. Conclusions Data from randomized clinical trials of 8 European ancestry and 2 black cohorts support the assumption that BST 1 containing locus on chromosome 4 is associated with β 1 ‐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.

Selection of antihypertensive treatment according to self-defined ethnicity is recommended by some guidelines but might be better guided by individual genotype rather than ethnicity or race. We compared the extent to which variation in blood pressure response across different ethnicities may be explained by genetic factors: genetically defined ancestry and gene variants at loci known to be associated with blood pressure. We analyzed data from 5 trials in which genotyping had been performed (n=4696) and in which treatment responses to β-blockers, angiotensin-converting enzyme inhibitors, angiotensin receptor blocker, thiazide or thiazide-like diuretic and calcium channel blocker were available. Genetically defined ancestry for proportion of African ancestry was computed using the 1000 genomes population database as a reference. Differences in response to the thiazide diuretic hydrochlorothiazide, the β-blockers atenolol and metoprolol, the angiotensin-converting enzyme inhibitor lisinopril, and the angiotensin receptor blocker candesartan were more closely associated to genetically defined ancestry than self-defined ethnicity in admixed subjects. A relatively small number of gene variants related to loci associated with drug-signaling pathways ( KCNK3, SULT1C3, AMH, PDE3A, PLCE1, PRKAG2 ) with large effect size (−3.5 to +3.5 mm Hg difference in response per allele) and differing allele frequencies in black versus white individuals explained a large proportion of the difference in response to candesartan and hydrochlorothiazide between these groups. These findings suggest that a genomic precision medicine approach can be used to individualize antihypertensive treatment within and across populations without recourse to surrogates of genetic structure such as self-defined ethnicity.

Abstractβ‐Blockers’ heart rate (HR)‐lowering effect is an important determinant of the effectiveness for this class of drugs, yet it is variable among β‐blocker–treated patients. To date, genetic studies have revealed several genetic signals associated with HR response to β‐blockers. However, these genetic signals have not been consistently replicated across multiple independent cohorts. Here we sought to use data from 3 hypertension clinical trials to validate single‐nucleotide polymorphisms (SNPs) previously associated with the HR response to β‐blockers. Using linear regression analysis, we investigated the effects of 6 SNPs in 3 genes, including ADRB1, ADRB2, and GNB3, relative to the HR response following β‐blocker used in the PEAR (n = 757), PEAR‐2 (n = 368), and INVEST (n = 1401) trials, adjusting for baseline HR, age, sex, and ancestry. Atenolol was used in PEAR and INVEST, and metoprolol was used in PEAR‐2. We found that rs1042714 and rs1042713 in ADRB2 were significantly associated with HR response to both β‐blockers in whites (rs1042714 C‐allele carriers, meta‐analysis β = –0.95 beats per minute [bpm], meta‐analysis P = 3×10–4; rs1042713 A‐allele carriers, meta‐analysis β = –1.15 bpm, meta‐analysis P = 2×10–3). In conclusion, the results of our analyses provide strong evidence to support the hypothesis that rs1042714 and rs1042713 in the ADRB2 gene are important predictors of HR response to cardioselective β‐blockade in hypertensive patient cohorts.

European Americans (EA) have a better antihypertensive response to β‐blockers when compared with African Americans, albeit with some variability. We undertook a genomewide association study to elucidate the underlying genetic determinants in EA contributing to this variability in blood pressure (BP) response. A discovery genomewide association study of change in BP post–metoprolol treatment was performed in EA participants (n = 201) from the Pharmacogenomic Evaluation of Antihypertensive Responses‐2 (PEAR‐2) study and tested for replication in the atenolol‐treated EA from the PEAR study (n = 233). Rs294610 in the FGD5, which encodes for FYVE, RhoGEF and PH Domain Containing 5, (expression quantitative trait loci for FGD5 in the small intestine) was significantly associated with increased diastolic BP response to β‐blockers in the PEAR‐2 study (P = 3.41 × 10−6, β = −2.70) and replicated (P = 0.01, β = −1.17) in the PEAR study. Post–meta‐analysis of these studies, an additional single nucleotide polymorphism rs45545233 in the SLC4A1, encoding for Solute Carrier Family 4 Member 1, (expression quantitative trait loci for dual specificity phosphatase 3 in the artery tibial) was identified that was significantly associated with a poor response to β‐blockers (P = 3.43 × 10−6, β = 4.57) and was replicated in the atenolol add‐on cohort (P = 0.007, β = 4.97). We identified variants in FGD5 and SLC4A1, which have been previously cited as candidate genes for hypertension, to be associated with a β‐blocker BP response in EA. Further elucidation is warranted of the underlying mechanisms of these variants and genes by which they influence the BP response to β‐blockers.

Abstract BACKGROUND Interindividual variability in blood pressure (BP) response to antihypertensives has been reported. Although plasma renin activity (PRA) is a potential biomarker for personalizing antihypertensive therapy in European American (EA) and African American (AA) hypertensives, clinical utility of PRA-guided prescribing is incompletely understood. METHODS Using systematic-phased approach, PRA’s clinical utility was assessed. After categorizing by baseline PRA, clinic systolic BP (SBP) responses to metoprolol and chlorthalidone were compared in 134 EAs and 102 AAs enrolled in the Pharmacogenomics Evaluation of Antihypertensive Responses-2 (PEAR-2) trial. Receiver operating characteristic (ROC) analysis was conducted in EAs. Data from PEAR-2 AAs were used to estimate an optimal PRA cut point using multivariable linear regression models. The derived cut point in AAs was tested in a meta-analysis of 2 independent AA cohorts, and its sensitivity and specificity were assessed. RESULTS EAs with PRA &lt; 0.65 ng/ml/hour had a greater decrease in SBP to chlorthalidone than metoprolol (by –15.9 mm Hg, adjusted P &lt; 0.0001), whereas those with PRA ≥ 0.65 ng/ml/hour had a greater decrease in SBP to metoprolol than chlorthalidone (by 3.3 mm Hg, adjusted P = 0.04). Area under ROC curve (0.69, P = 0.0001) showed that PRA can predict SBP response among EAs. However, we observed no association between PRA and SBP response in PEAR-2 AAs. Among independent AA cohorts, those with PRA ≥ 1.3 ng/ml/hour (PEAR-2-derived cut point) responded better to atenolol/candesartan than hydrochlorothiazide (meta-analysis P = 0.01). However, sensitivity of the derived cut point was 10%. CONCLUSIONS PRA at the previously established 0.60–0.65 ng/ml/hour cut point is an effective predictive biomarker of BP response in EAs. However, we were unable to identify PRA cut point that could be used to guide antihypertensive selection in AAs. TRIAL REGISTRATION NCT01203852, NCT00246519, NCT00005520.