While atenolol is an effective antihypertensive agent, its use is also associated with adverse events including hyperglycemia and incident diabetes that may offset the benefits of blood pressure lowering. By combining metabolomic and genomic data acquired from hypertensive individuals treated with atenolol, it may be possible to better understand the pathways that most impact the development of an adverse glycemic state. To identify biomarkers that can help predict susceptibility to blood glucose excursions during exposure to atenolol. Plasma samples acquired from 234 Caucasian participants treated with atenolol in the Pharmacogenomic Evaluation of Antihypertensive Responses trial were analyzed by gas chromatography Time-Of-Flight Mass Spectroscopy. Metabolomics and genomics data were integrated by first correlating participant's metabolomic profiles to change in glucose after treatment with atenolol, and then incorporating genotype information from genes involved in metabolite pathways associated with glucose response. Our findings indicate that the baseline level of β-alanine was associated with glucose change after treatment with atenolol (Q = 0.007, β = 2.97 mg/dL). Analysis of genomic data revealed that carriers of the G allele for SNP rs2669429 in gene DPYS, which codes for dihydropyrimidinase, an enzyme involved in β-alanine formation, had significantly higher glucose levels after treatment with atenolol when compared with non-carriers (Q = 0.05, β = 2.76 mg/dL). This finding was replicated in participants who received atenolol as an add-on therapy (P = 0.04, β = 1.86 mg/dL). These results suggest that β-alanine and rs2669429 may be predictors of atenolol-induced hyperglycemia in Caucasian individuals and further investigation is warranted.

Hydrochlorothiazide is among the most commonly prescribed antihypertensives; yet, <50% of hydrochlorothiazide-treated patients achieve blood pressure (BP) control. Herein, we integrated metabolomic and genomic profiles of hydrochlorothiazide-treated patients to identify novel genetic markers associated with hydrochlorothiazide BP response. The primary analysis included 228 white hypertensives treated with hydrochlorothiazide from the Pharmacogenomic Evaluation of Antihypertensive Responses (PEAR) study. Genome-wide analysis was conducted using Illumina Omni 1 mol/L-Quad Chip, and untargeted metabolomics was performed on baseline fasting plasma samples using a gas chromatography-time-of-flight mass spectrometry platform. We found 13 metabolites significantly associated with hydrochlorothiazide systolic BP (SBP) and diastolic BP (DBP) responses (false discovery rate, <0.05). In addition, integrating genomic and metabolomic data revealed 3 polymorphisms (rs2727563 PRKAG2, rs12604940 DCC, and rs13262930 EPHX2) along with arachidonic acid, converging in the netrin signaling pathway (P=1×10(-5)), as potential markers, significantly influencing hydrochlorothiazide BP response. We successfully replicated the 3 genetic signals in 212 white hypertensives treated with hydrochlorothiazide and created a response score by summing their BP-lowering alleles. We found patients carrying 1 response allele had a significantly lower response than carriers of 6 alleles (∆SBP/∆DBP: -1.5/1.2 versus -16.3/-10.4 mm Hg, respectively, SBP score, P=1×10(-8) and DBP score, P=3×10(-9)). This score explained 11.3% and 11.9% of the variability in hydrochlorothiazide SBP and DBP responses, respectively, and was further validated in another independent study of 196 whites treated with hydrochlorothiazide (DBP score, P=0.03; SBP score, P=0.07). This study suggests that PRKAG2, DCC, and EPHX2 might be important determinants of hydrochlorothiazide BP response.

Although hydrochlorothiazide (HCTZ) is a well-established first-line antihypertensive in the United States, <50% of HCTZ treated patients achieve blood pressure (BP) control. Thus, identifying biomarkers that could predict the BP response to HCTZ is critically important. In this study, we utilized metabolomics, genomics, and lipidomics to identify novel pathways and biomarkers associated with HCTZ BP response. First, we conducted a pathway analysis for 13 metabolites we recently identified to be significantly associated with HCTZ BP response. From this analysis, we found the sphingolipid metabolic pathway as the most significant pathway (P=5.8E-05). Testing 78 variants, within 14 genes involved in the sphingolipid metabolic canonical pathway, with the BP response to HCTZ identified variant rs6078905, within the SPTLC3 gene, as a novel biomarker significantly associated with the BP response to HCTZ in whites (n=228). We found that rs6078905 C-allele carriers had a better BP response to HCTZ versus noncarriers (∆SBP/∆DBP: -11.4/-6.9 versus -6.8/-3.5 mm Hg; ∆SBP P=6.7E-04; ∆DBP P=4.8E-04). Additionally, in blacks (n=148), we found genetic signals in the SPTLC3 genomic region significantly associated with the BP response to HCTZ (P<0.05). Last, we observed that rs6078905 significantly affects the baseline level of 4 sphingomyelins (N24:2, N24:3, N16:1, and N22:1; false discovery rate <0.05), from which N24:2 sphingomyelin has a significant correlation with both HCTZ DBP-response (r=-0.42; P=7E-03) and SBP-response (r=-0.36; P=2E-02). This study provides insight into potential pharmacometabolomic and genetic mechanisms underlying HCTZ BP response and suggests that SPTLC3 is a potential determinant of the BP response to HCTZ. URL: http://www.clinicaltrials.gov. Unique identifier: NCT00246519.