Gums, John G.

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.

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.

Background Thiazide and thiazide‐like diuretics are first‐line medications for treating uncomplicated hypertension. However, their use has been associated with adverse metabolic events, including hyperglycemia and incident diabetes mellitus, with incompletely understood mechanisms. Our goal was to identify genomic variants associated with thiazide‐like diuretic/chlorthalidone‐induced glucose change. Methods and Results Genome‐wide analysis of glucose change after treatment with chlorthalidone was performed by race among the white (n=175) and black (n=135) participants from the PEAR ‐2 (Pharmacogenomic Evaluation of Antihypertensive Responses‐2). Single‐nucleotide polymorphisms with P <5×10 −8 were further prioritized using in silico analysis based on their expression quantitative trait loci function. Among blacks, an intronic single‐nucleotide polymorphism (rs9943291) in the HMGCS 2 was associated with increase in glucose levels following chlorthalidone treatment (ß=12.5; P =4.17×10 −8 ). G‐allele carriers of HMGCS 2 had higher glucose levels (glucose change=+16.29 mg/dL) post chlorthalidone treatment compared with noncarriers of G allele (glucose change=+2.80 mg/dL). This association was successfully replicated in an independent replication cohort of hydrochlorothiazide‐treated participants from the PEAR study (ß=5.54; P =0.023). A meta‐analysis of the 2 studies was performed by race in Meta‐Analysis Helper, where this single‐nucleotide polymorphism, rs9943291, was genome‐wide significant with a meta‐analysis P value of 3.71×10 −8 . HMGCS 2 , a part of the HMG ‐CoA synthase family, is important for ketogenesis and cholesterol synthesis pathways that are essential in glucose homeostasis. Conclusions These results suggest that HMGCS 2 is a promising candidate gene involved in chlorthalidone and Hydrochlorothiazide (HCTZ)‐induced glucose change. This may provide insights into the mechanisms involved in thiazide‐induced hyperglycemia that may ultimately facilitate personalized approaches to antihypertensive selection for hypertension treatment. Clinical Trial Registration URL : http://www.clinicaltrials.gov . Unique identifiers: NCT 00246519 and NCT 01203852.

Background For many indications, the negative chronotropic effect of β‐blockers is important to their efficacy, yet the heart rate ( HR ) response to β‐blockers varies. Herein, we sought to use a genome‐wide association approach to identify novel single nucleotide polymorphisms ( SNP s) associated with HR response to β‐blockers. Methods and Results We first performed 4 genome‐wide association analyses for HR response to atenolol (a β1‐adrenergic receptor blocker) as: (1) monotherapy or (2) add‐on therapy, in 426 whites and 273 blacks separately from the PEAR (Pharmacogenomic Evaluation of Antihypertensive Responses) study. A meta‐analysis was then performed between the genome‐wide association analysis performed in PEAR atenolol monotherapy and add‐on therapy, in each race separately, using the inverse variance method assuming fixed effects. From this analysis, SNP s associated with HR response to atenolol at a P <1E‐05 were tested for replication in whites (n=200) and blacks (n=168) treated with metoprolol (a β1‐adrenergic receptor blocker). From the genome‐wide association meta‐analyses, SNP rs17117817 near olfactory receptor family10 subfamily‐p‐member1 ( OR 10P1 ), and SNP rs2364349 in sorting nexin‐9 ( SNX 9 ) replicated in blacks. The combined studies meta‐analysis P values for the rs17117817 and rs2364349 reached genome‐wide significance (rs17117817G‐allele; Meta‐β=5.53 beats per minute, Meta‐ P =2E‐09 and rs2364349 A‐allele; Meta‐β=3.5 beats per minute, Meta‐ P =1E‐08). Additionally, SNP s in the OR 10P1 and SNX 9 gene regions were also associated with HR response in whites. Conclusions This study highlights OR 10P1 and SNX 9 as novel genes associated with changes in HR in response to β‐blockers. Clinical Trial Registration URL : http://www.clinicaltrials.gov . Unique identifier: NCT 00246519.