Recent CYP2D6 phenotype standardization efforts by CYP2D6 activity score (AS) are based on limited pharmacokinetic (PK) and pharmacodynamic (PD) data. Using data from two independent clinical trials of metoprolol, we compared metoprolol PK and PD across CYP2D6 AS with the goal of determining whether the PK and PD data support the new phenotype classification. S-metoprolol apparent oral clearance (CLo), adjusted for clinical factors, was correlated with CYP2D6 AS (P < 0.001). The natural log of CLo was lower with an AS of 1 (7.6 ± 0.4 mL/minute) vs. 2-2.25 (8.3 ± 0.6 mL/minute; P = 0.012), similar between an AS of 1 and 1.25-1.5 (7.8 ± 0.5 mL/minute; P = 0.702), and lower with an AS of 1.25-1.5 vs. 2-2.25 (P = 0.03). There was also a greater reduction in heart rate with metoprolol among study participants with AS of 1 (-10.8 ± 5.5) vs. 2-2.25 (-7.1 ± 5.6; P < 0.001) and no significant difference between those with an AS of 1 and 1.25-1.5 (-9.2 ± 4.7; P = 0.095). These data highlight linear trends among CYP2D6 AS and metoprolol PK and PD, but inconsistencies with the phenotypes assigned by AS based on the current standards. Overall, this case study with metoprolol suggests that utilizing CYP2D6 AS, instead of collapsing AS into phenotype categories, may be the most precise approach for utilizing CYP2D6 pharmacogenomics in clinical practice.

Plasma renin activity (PRA) is a predictive biomarker of blood pressure (BP) response to antihypertensives in European-American hypertensive patients. We aimed to identify the metabolic signatures of baseline PRA and the linkages with BP response to β-blockers and thiazides. Using data from the Pharmacogenomic Evaluation of Antihypertensive Responses-2 (PEAR-2) trial, multivariable linear regression adjusting for age, sex and baseline systolic-BP (SBP) was performed on European-American individuals treated with metoprolol (n = 198) and chlorthalidone (n = 181), to test associations between 856 metabolites and baseline PRA. Metabolites with a false discovery rate (FDR) < 0.05 or p < 0.01 were tested for replication in 463 European-American individuals treated with atenolol or hydrochlorothiazide. Replicated metabolites were then tested for validation based on the directionality of association with BP response. Sixty-three metabolites were associated with baseline PRA, of which nine, including six lipids, were replicated. Of those replicated, two metabolites associated with higher baseline PRA were validated: caprate was associated with greater metoprolol SBP response (β = -1.7 ± 0.6, p = 0.006) and sphingosine-1-phosphate was associated with reduced hydrochlorothiazide SBP response (β = 7.6 ± 2.8, p = 0.007). These metabolites are clustered with metabolites involved in sphingolipid, phospholipid, and purine metabolic pathways. The identified metabolic signatures provide insights into the mechanisms underlying BP response.

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.

Chlorthalidone (CTD) is more potent than hydrochlorothiazide (HCTZ) in reducing blood pressure (BP) in hypertensive patients, though both are plagued with BP response variability. However, there is a void in the literature regarding the genetic determinants contributing to the variability observed in BP response to CTD. We performed a discovery genome wide association analysis of BP response post CTD treatment in African Americans (AA) and European Americans (EA) from the Pharmacogenomic Evaluation of Antihypertensive Responses-2 (PEAR-2) study and replication in an independent cohort of AA and EA treated with HCTZ from the PEAR study, followed by a race specific meta-analysis of the two studies. Successfully replicated SNPs were further validated in beta-blocker treated participants from PEAR-2 and PEAR for opposite direction of association. The replicated and validated signals were further evaluated by protein-protein interaction network analysis. An intronic SNP rs79237970 in the WDR92 (eQTL for PPP3R1) was significantly associated with better DBP response to CTD (p = 5.76 × 10-6, β = -15.75) in the AA cohort. This SNP further replicated in PEAR (p = 0.00046, β = -9.815) with a genome wide significant meta-analysis p-value of 8.49 × 10-9. This variant was further validated for opposite association in two β-blockers treated cohorts from PEAR-2 metoprolol (p = 9.9 × 10-3, β = 7.47) and PEAR atenolol (p = 0.04, β = 4.36) for association with DBP. Studies have implicated WDR92 in coronary artery damage. PPP3R1 is the regulatory subunit of the calcineurin complex. Use of calcineurin inhibitors is associated with HTN. Studies have also shown polymorphisms in PPP3R1 to be associated with ventricular hypertrophy in AA hypertensive patients. Protein-protein interaction analysis further identified important hypertension related pathways such as inositol phosphate-mediated signaling and calcineurin-NFAT signaling cascade as important biological process associated with PPP3R1 which further strengthen the potential importance of this signal. These data collectively suggest that WDR92 and PPP3R1 are novel candidates that may help explain the genetic underpinnings of BP response of thiazide and thiazide-like diuretics and help identify the patients better suited for thiazide and thiazide-like diuretics compared to β-blockers for improved BP management. This may further help advance personalized approaches to antihypertensive therapy.

Hypertension (HTN) disproportionately affects African Americans (AAs), who respond better to thiazide diuretics than other antihypertensives. Variants of the APOL1 gene found in AAs are associated with a higher rate of kidney disease and play a complex role in cardiovascular disease.