Citations

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.

Abstract Introduction The CYP2D6 enzyme metabolizes opioids commonly prescribed for cancer‐related pain, and CYP2D6 polymorphisms may contribute to variability in opioid response. We evaluated the feasibility of implementing CYP2D6 ‐guided opioid prescribing for patients with cancer and reported pilot outcome data. Methods Adult patients from two cancer centers were prospectively enrolled into a hybrid implementation‐effectiveness clinical trial and randomized to CYP2D6 ‐genotype‐guided opioid selection, with clinical recommendations, or usual care. Implementation metrics, including provider response, medication changes consistent with recommendations, and patient‐reported pain and symptom scores at baseline and up to 8 weeks, were assessed. Results Most (87/114, 76%) patients approached for the study agreed to participate. Of 85 patients randomized, 71% were prescribed oxycodone at baseline. The median (range) time to receive CYP2D6 test results was 10 (3–37) days; 24% of patients had physicians acknowledge genotype results in a clinic note. Among patients with CYP2D6 ‐genotype‐guided recommendations to change therapy ( n = 11), 18% had a change congruent with recommendations. Among patients who completed baseline and follow‐up questionnaires ( n = 48), there was no difference in change in mean composite pain score (−1.01 ± 2.1 vs. −0.41 ± 2.5; p = 0.19) or symptom severity at last follow‐up (3.96 ± 2.18 vs. 3.47 ± 1.78; p = 0.63) between the usual care arm ( n = 26) and genotype‐guided arm ( n = 22), respectively. Conclusion Our study revealed high acceptance of pharmacogenetic testing as part of a clinical trial among patients with cancer pain. However, provider response to genotype‐guided recommendations was low, impacting assessment of pain‐related outcomes. Addressing barriers to utility of pharmacogenetics results and clinical recommendations will be critical for implementation success.

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.

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.

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.

Background PCSK9 (proprotein convertase subtilisin/kexin type 9) inhibitors effectively lower LDL (low‐density lipoprotein) cholesterol and have been shown to reduce cardiovascular outcomes in high‐risk patients. We used real‐world electronic health record data to characterize use of PCSK9 inhibitors, in addition to standard therapies, according to cardiovascular risk status. Methods and Results Data were obtained from 18 health systems with data marts within the National Patient‐Centered Clinical Research Network (PCORnet) using a common data model. Participating sites identified >17.5 million adults, of whom 3.6 million met study criteria. Patients were categorized into 3 groups: (1) dyslipidemia, (2) untreated LDL ≥130 mg/dL, and (3) coronary artery disease or coronary heart disease. Demographics, comorbidities, estimated 10‐year atherosclerotic cardiovascular disease risk, and lipid‐lowering pharmacotherapies were summarized for each group. Participants’ average age was 62 years, 50% were female, and 11% were black. LDL cholesterol ranged from 85 to 151 mg/dL. Among patients in groups 1 and 3, 54% received standard lipid‐lowering therapies and a PCSK9 inhibitor was prescribed in <1%. PCSK9 inhibitor prescribing was greatest for patients with coronary artery disease or coronary heart disease and, although prescribing increased during the study period, overall PCSK9 inhibitor prescribing was low. Conclusions We successfully used electronic health record data from 18 PCORnet data marts to identify >3.6 million patients meeting criteria for 3 patient groups. Approximately half of patients had been prescribed lipid‐lowering medication, but <1% were prescribed PCSK9 inhibitors. PCSK9 inhibitor prescribing increased over time for patients with coronary artery disease or coronary heart disease but not for those with dyslipidemia.

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.