Genome-wide study of resistant hypertension identified from electronic health records

Whole genome sequence analysis of apparent treatment resistant hypertension status in participants from the Trans-Omics for Precision Medicine program

Introduction: Apparent treatment-resistant hypertension (aTRH) is characterized by the use of four or more antihypertensive (AHT) classes to achieve blood pressure (BP) control. In the current study, we conducted single-variant and gene-based analyses of aTRH among individuals from 12 Trans-Omics for Precision Medicine cohorts with whole-genome sequencing data. Methods: Cases were defined as individuals treated for hypertension (HTN) taking three different AHT classes, with average systolic BP ≥ 140 or diastolic BP ≥ 90 mmHg, or four or more medications regardless of BP (n = 1,705). A normotensive control group was defined as individuals with BP < 140/90 mmHg (n = 22,079), not on AHT medication. A second control group comprised individuals who were treatment responsive on one AHT medication with BP < 140/ 90 mmHg (n = 5,424). Logistic regression with kinship adjustment using the Scalable and Accurate Implementation of Generalized mixed models (SAIGE) was performed, adjusting for age, sex, and genetic ancestry. We assessed variants using SKAT-O in rare-variant analyses. Single-variant and gene-based tests were conducted in a pooled multi-ethnicity stratum, as well as self-reported ethnic/racial strata (European and African American). Results: One variant in the known HTN locus, KCNK3, was a top finding in the multi-ethnic analysis (p = 8.23E-07) for the normotensive control group [rs12476527, odds ratio (95% confidence interval) = 0.80 (0.74-0.88)]. This variant was replicated in the Vanderbilt University Medical Center's DNA repository data. Aggregate gene-based signals included the genes AGTPBP, MYL4, PDCD4, BBS9, ERG, and IER3. Discussion: Additional work validating these loci in larger, more diverse populations, is warranted to determine whether these regions influence the pathobiology of aTRH.

Prospective, multi-site study of healthcare utilization after actionable monogenic findings from clinical sequencing

As large-scale genomic screening becomes increasingly prevalent, understanding the influence of actionable results on healthcare utilization is key to estimating the potential long-term clinical impact. The eMERGE network sequenced individuals for actionable genes in multiple genetic conditions and returned results to individuals, providers, and the electronic health record. Differences in recommended health services (laboratory, imaging, and procedural testing) delivered within 12 months of return were compared among individuals with pathogenic or likely pathogenic (P/LP) findings to matched individuals with negative findings before and after return of results. Of 16,218 adults, 477 unselected individuals were found to have a monogenic risk for arrhythmia (n = 95), breast cancer (n = 96), cardiomyopathy (n = 95), colorectal cancer (n = 105), or familial hypercholesterolemia (n = 86). Individuals with P/LP results more frequently received services after return (43.8%) compared to before return (25.6%) of results and compared to individuals with negative findings (24.9%; p < 0.0001). The annual cost of qualifying healthcare services increased from an average of $162 before return to $343 after return of results among the P/LP group (p < 0.0001); differences in the negative group were non-significant. The mean difference-in-differences was $149 (p < 0.0001), which describes the increased cost within the P/LP group corrected for cost changes in the negative group. When stratified by individual conditions, significant cost differences were observed for arrhythmia, breast cancer, and cardiomyopathy. In conclusion, less than half of individuals received billed health services after monogenic return, which modestly increased healthcare costs for payors in the year following return.

Opportunities and challenges for biomarker discovery using electronic health record data

Electronic health records (EHRs) have become increasingly relied upon as a source for biomedical research. One important research application of EHRs is the identification of biomarkers associated with specific patient states, especially within complex conditions. However, using EHRs for biomarker identification can be challenging because the EHR was not designed with research as the primary focus. Despite this challenge, the EHR offers huge potential for biomarker discovery research to transform our understanding of disease etiology and treatment and generate biological insights informing precision medicine initiatives. This review paper provides an in-depth analysis of how EHR data is currently used for phenotyping and identifying molecular biomarkers, current challenges and limitations, and strategies we can take to mitigate challenges going forward.

Importance of Diversity in Precision Medicine: Generalizability of Genetic Associations Across Ancestry Groups Toward Better Identification of Disease Susceptibility Variants

Genome-wide association studies (GWAS) revolutionized our understanding of common genetic variation and its impact on common human disease and traits. Developed and adopted in the mid-2000s, GWAS led to searchable genotype-phenotype catalogs and genome-wide datasets available for further data mining and analysis for the eventual development of translational applications. The GWAS revolution was swift and specific, including almost exclusively populations of European descent, to the neglect of the majority of the world's genetic diversity. In this narrative review, we recount the GWAS landscape of the early years that established a genotype-phenotype catalog that is now universally understood to be inadequate for a complete understanding of complex human genetics. We then describe approaches taken to augment the genotype-phenotype catalog, including the study populations, collaborative consortia, and study design approaches aimed to generalize and then ultimately discover genome-wide associations in non-European descent populations. The collaborations and data resources established in the efforts to diversify genomic findings undoubtedly provide the foundations of the next chapters of genetic association studies with the advent of budget-friendly whole-genome sequencing.

Genome-wide Association Studies Categorized by Class of Antihypertensive Drugs Reveal Complex Pathogenesis of Hypertension with Drug Resistance

Resistant hypertension is defined as uncontrolled blood pressure (BP) despite the use of three or more antihypertensive drugs of different classes. Although genetic factors may greatly contribute to hypertension with resistance to multiple drug classes, more than for general hypertension, its pathogenesis remains unknown. To reveal the genetic background of resistant hypertension, we categorized 32,239 patients whose data were obtained from the BioBank Japan Project, by prescription of 7 classes of antihypertensive drugs and performed genome-wide association studies (GWAS). Our GWAS identified four loci with significant association (P < 5 × 10-8 ): rs6445583 in CACNA1D and rs12308051 in the intergenic region on chromosome 12 for angiotensin II receptor blockers, rs35497065 in FOXA3 for calcium channel blockers, and rs11066280 in HECTD4 for αβ-blockers. Because these loci are known to be susceptibility loci for hypertension and/or BP, our results indicate that resistant hypertension is caused by a combination of excessive BP and drug resistance to each antihypertensive pharmacological class. Furthermore, to investigate the genetic difference between BP traits and the treatment effectiveness of antihypertensive drugs, we performed gene-set analysis and calculated the genetic correlation continuously. Most of the genetic factors were in common between BP traits and antihypertensive effectiveness, but it seems that the genetic architecture of the drug response to antihypertensive treatment is more complicated than BP traits. This corresponds to the well-known mosaic theory of hypertension. Our findings reveal the complex pathogenesis of hypertension with resistance to multiple classes of antihypertensive drugs.

A flexible symbolic regression method for constructing interpretable clinical prediction models

Machine learning (ML) models trained for triggering clinical decision support (CDS) are typically either accurate or interpretable but not both. Scaling CDS to the panoply of clinical use cases while mitigating risks to patients will require many ML models be intuitively interpretable for clinicians. To this end, we adapted a symbolic regression method, coined the feature engineering automation tool (FEAT), to train concise and accurate models from high-dimensional electronic health record (EHR) data. We first present an in-depth application of FEAT to classify hypertension, hypertension with unexplained hypokalemia, and apparent treatment-resistant hypertension (aTRH) using EHR data for 1200 subjects receiving longitudinal care in a large healthcare system. FEAT models trained to predict phenotypes adjudicated by chart review had equivalent or higher discriminative performance (p < 0.001) and were at least three times smaller (p < 1 × 10^-6) than other potentially interpretable models. For aTRH, FEAT generated a six-feature, highly discriminative (positive predictive value = 0.70, sensitivity = 0.62), and clinically intuitive model. To assess the generalizability of the approach, we tested FEAT on 25 benchmark clinical phenotyping tasks using the MIMIC-III critical care database. Under comparable dimensionality constraints, FEAT's models exhibited higher area under the receiver-operating curve scores than penalized linear models across tasks (p < 6 × 10^-6). In summary, FEAT can train EHR prediction models that are both intuitively interpretable and accurate, which should facilitate safe and effective scaling of ML-triggered CDS to the panoply of potential clinical use cases and healthcare practices.

A genome-wide association study identifies a novel association between SDC3 and apparent treatment-resistant hypertension

BACKGROUND

Compared with patients who require fewer antihypertensive agents, those with apparent treatment-resistant hypertension (aTRH) are at increased risk for cardiovascular and all-cause mortality, independent of blood pressure control. However, the etiopathogenesis of aTRH is still poorly elucidated.

METHODS

We performed a genome-wide association study (GWAS) in first cohort including 586 aTRHs and 871 healthy controls. Next, expression quantitative trait locus (eQTL) analysis was used to identify genes that are regulated by single nucleotide polymorphisms (SNPs) derived from the GWAS. Then, we verified the genes obtained from the eQTL analysis in the validation cohort including 65 aTRHs, 96 hypertensives, and 100 healthy controls through gene expression profiling analysis and real-time quantitative polymerase chain reaction (RT-qPCR) assay.

RESULTS

The GWAS in first cohort revealed four suggestive loci (1p35, 4q13.2-21.1, 5q22-23.2, and 15q11.1-q12) represented by 23 SNPs. The 23 significant SNPs were in or near LAPTM5, SDC3, UGT2A1, FTMT, and NIPA1. eQTL analysis uncovered 14 SNPs in 1p35 locus all had same regulation directions for SDC3 and LAPTM5. The disease susceptible alleles of SNPs in 1p35 locus were associated with lower gene expression for SDC3 and higher gene expression for LAPTM5. The disease susceptible alleles of SNPs in 4q13.2-21.1 were associated with higher gene expression for UGT2B4. GTEx database did not show any statistically significant eQTLs between the SNPs in 5q22-23.2 and 15q11.1-q12 loci and their influenced genes. Then, gene expression profiling analysis in the validation cohort confirmed lower expression of SDC3 in aTRH but no significant differences on LAPTM5 and UGT2B4, when compared with controls and hypertensives, respectively. RT-qPCR assay further verified the lower expression of SDC3 in aTRH.

CONCLUSIONS

Our study identified a novel association of SDC3 with aTRH, which contributes to the elucidation of its etiopathogenesis and provides a promising therapeutic target.

Application of omics in hypertension and resistant hypertension

Hypertension is a major modifiable risk factor that affects the global health burden. Despite the availability of multiple antihypertensive drugs, blood pressure is often not optimally controlled. The prevalence of true resistant hypertension in treated hypertensive patients is ~2-20%, and these patients are at higher risk for adverse events and poor clinical outcomes. Therefore, an in-depth dissection of the pathophysiological mechanisms of hypertension and resistant hypertension is needed to identify more effective targets for regulating blood pressure. Omics technologies, such as genomics, transcriptomics, proteomics, metabolomics, and microbiomics, can accurately present the characteristics of organisms at varying molecular levels. Integrative omics can further reveal the network of interactions between molecular levels and provide a complete dynamic view of the organism. In this review, we describe the applications, progress, and challenges of omics technologies in hypertension. Specifically, we discuss the application of omics in resistant hypertension. We believe that omics approaches will produce a better understanding of the pathogenesis of hypertension and resistant hypertension and improve diagnostic and therapeutic strategies, thus increasing rates of blood pressure control and reducing the public health burden of hypertension.

Replication of European hypertension associations in a case-control study of 9,534 African Americans

Objective

Hypertension is more prevalent in African Americans (AA) than other ethnic groups. Genome-wide association studies (GWAS) have identified loci associated with hypertension and other cardio-metabolic traits like type 2 diabetes, coronary artery disease, and body mass index (BMI), however the AA population is underrepresented in these studies. In this study, we examined a large AA cohort for the generalizability of 14 Metabochip array SNPs with previously reported European hypertension associations.

Methods

To evaluate associations, we analyzed genotype data of 14 SNPs for their associations with a diagnosis of hypertension, systolic blood pressure (SBP), and diastolic blood pressure (DBP) in a case-control study of an AA population (N = 9,534). We also performed an age-stratified analysis (>30, 30≥59 and ≥60 years) following the hypertension definition described by the 8th Joint National Committee (JNC). Associations were adjusted for BMI, age, age², sex, clinical confounders, and genetic ancestry using multivariable regression models to estimate odds ratios (ORs) and beta-coefficients. Analyses stratified by sex were also conducted. Meta-analyses (including both BioVU and COGENT-BP cohorts) were performed using a random-effects model.

Results

We found rs880315 to be associated with systolic hypertension (SBP≥140 mmHg) in the entire cohort (OR = 1.14, p = 0.003) and within women only (OR = 1.16, p = 0.012). Variant rs17080093 associated with lower SBP and DBP (β = -2.99, p = 0.0352 and - β = 1.69, p = 0.0184) among younger individuals, particularly in younger women (β = -3.92, p = 0.0025 and β = -1.87, p = 0.0241 for SBP and DBP respectively). SNP rs1530440 associated with higher SBP and DBP measurements (younger individuals β = 4.1, p = 0.039 and β = 2.5, p = 0.043 for SBP and DBP; (younger women β = 4.5, p = 0.025 and β = 2.9, p = 0.028 for SBP and DBP), and hypertension risk in older women (OR = 1.4, p = 0.050). rs16948048 increases hypertension risk in younger individuals (OR = 1.31, p = 0.011). Among mid-age women rs880315 associated with higher risk of hypertension (OR = 1.20, p = 0.027). rs1361831 associated with DBP (β = -1.96, p = 0.02) among individuals older than 60 years. rs3096277 increases hypertension risk among older individuals (OR = 1.26 p = 0.0015), however, this variant also reduces SBP among younger women (β = -2.63, p = 0.0102).

Conclusion

These findings suggest that European-descent and AA populations share genetic loci that contribute to blood pressure traits and hypertension. However, the OR and beta-coefficient estimates differ, and some are age-dependent. Additional genetic studies of hypertension in AA are warranted to identify new loci associated with hypertension and blood pressure traits in this population.

A genome-wide association study identifies a novel candidate locus at the DLGAP1 gene with susceptibility to resistant hypertension in the Japanese population

Numerous genetic variants associated with hypertension and blood pressure are known, but there is a paucity of evidence from genetic studies of resistant hypertension, especially in Asian populations. To identify novel genetic loci associated with resistant hypertension in the Japanese population, we conducted a genome-wide association study with 2705 resistant hypertension cases and 21,296 mild hypertension controls, all from BioBank Japan. We identified one novel susceptibility candidate locus, rs1442386 on chromosome 18p11.3 (DLGAP1), achieving genome-wide significance (odds ratio (95% CI) = 0.85 (0.81–0.90), P = 3.75 × 10⁻⁸) and 18 loci showing suggestive association, including rs62525059 of 8q24.3 (CYP11B2) and rs3774427 of 3p21.1 (CACNA1D). We further detected biological processes associated with resistant hypertension, including chemical synaptic transmission, regulation of transmembrane transport, neuron development and neurological system processes, highlighting the importance of the nervous system. This study provides insights into the etiology of resistant hypertension in the Japanese population.

The Role of Electronic Health Records in Advancing Genomic Medicine

Recent advances in genomic technology and widespread adoption of electronic health records (EHRs) have accelerated the development of genomic medicine, bringing promising research findings from genome science into clinical practice. Genomic and phenomic data, accrued across large populations through biobanks linked to EHRs, have enabled the study of genetic variation at a phenome-wide scale. Through new quantitative techniques, pleiotropy can be explored with phenome-wide association studies, the occurrence of common complex diseases can be predicted using the cumulative influence of many genetic variants (polygenic risk scores), and undiagnosed Mendelian syndromes can be identified using EHR-based phenotypic signatures (phenotype risk scores). In this review, we trace the role of EHRs from the development of genome-wide analytic techniques to translational efforts to test these new interventions to the clinic. Throughout, we describe the challenges that remain when combining EHRs with genetics to improve clinical care.

Hypertension prevalence in the All of Us Research Program among groups traditionally underrepresented in medical research

The All of Us Research Program was designed to enable broad-based precision medicine research in a cohort of unprecedented scale and diversity. Hypertension (HTN) is a major public health concern. The validity of HTN data and definition of hypertension cases in the All of Us (AoU) Research Program for use in rule-based algorithms is unknown. In this cross-sectional, population-based study, we compare HTN prevalence in the AoU Research Program to HTN prevalence in the 2015-2016 National Health and Nutrition Examination Survey (NHANES). We used AoU baseline data from patient (age ≥ 18) measurements (PM), surveys, and electronic health record (EHR) blood pressure measurements. We retrospectively examined the prevalence of HTN in the EHR cohort using Systemized Nomenclature of Medicine (SNOMED) codes and blood pressure medications recorded in the EHR. We defined HTN as the participant having at least 2 HTN diagnosis/billing codes on separate dates in the EHR data AND at least one HTN medication. We calculated an age-standardized HTN prevalence according to the age distribution of the U.S. Census, using 3 groups (18-39, 40-59, and ≥ 60). Among the 185,770 participants enrolled in the AoU Cohort (mean age at enrollment = 51.2 years) available in a Researcher Workbench as of October 2019, EHR data was available for at least one SNOMED code from 112,805 participants, medications for 104,230 participants, and 103,490 participants had both medication and SNOMED data. The total number of persons with SNOMED codes on at least two distinct dates and at least one antihypertensive medication was 33,310 for a crude prevalence of HTN of 32.2%. AoU age-adjusted HTN prevalence was 27.9% using 3 groups compared to 29.6% in NHANES. The AoU cohort is a growing source of diverse longitudinal data to study hypertension nationwide and develop precision rule-based algorithms for use in hypertension treatment and prevention research. The prevalence of hypertension in this cohort is similar to that in prior population-based surveys.

Rapid response to the alpha-1 adrenergic agent phenylephrine in the perioperative period is impacted by genomics and ancestry

The emergence of genomic data in biobanks and health systems offers new ways to derive medically important phenotypes, including acute phenotypes occurring during inpatient clinical care. Here we study the genetic underpinnings of the rapid response to phenylephrine, an α1-adrenergic receptor agonist commonly used to treat hypotension during anesthesia and surgery. We quantified this response by extracting blood pressure (BP) measurements 5 min before and after the administration of phenylephrine. Based on this derived phenotype, we show that systematic differences exist between self-reported ancestry groups: European-Americans (EA; n = 1387) have a significantly higher systolic response to phenylephrine than African-Americans (AA; n = 1217) and Hispanic/Latinos (HA; n = 1713) (31.3% increase, p value < 6e-08 and 22.9% increase, p value < 5e-05 respectively), after adjusting for genetic ancestry, demographics, and relevant clinical covariates. We performed a genome-wide association study to investigate genetic factors underlying individual differences in this derived phenotype. We discovered genome-wide significant association signals in loci and genes previously associated with BP measured in ambulatory settings, and a general enrichment of association in these genes. Finally, we discovered two low frequency variants, present at ~1% in EAs and AAs, respectively, where patients carrying one copy of these variants show no phenylephrine response. This work demonstrates our ability to derive a quantitative phenotype suited for comparative statistics and genome-wide association studies from dense clinical and physiological measures captured for managing patients during surgery. We identify genetic variants underlying non response to phenylephrine, with implications for preemptive pharmacogenomic screening to improve safety during surgery.

Optimizing identification of resistant hypertension: Computable phenotype development and validation

PURPOSE

Computable phenotypes are constructed to utilize data within the electronic health record (EHR) to identify patients with specific characteristics; a necessary step for researching a complex disease state. We developed computable phenotypes for resistant hypertension (RHTN) and stable controlled hypertension (HTN) based on the National Patient-Centered Clinical Research Network (PCORnet) common data model (CDM). The computable phenotypes were validated through manual chart review.

METHODS

We adapted and refined existing computable phenotype algorithms for RHTN and stable controlled HTN to the PCORnet CDM in an adult HTN population from the OneFlorida Clinical Research Consortium (2015-2017). Two independent reviewers validated the computable phenotypes through manual chart review of 425 patient records. We assessed precision of our computable phenotypes through positive predictive value (PPV) and test validity through interrater reliability (IRR).

RESULTS

Among the 156 730 HTN patients in our final dataset, the final computable phenotype algorithms identified 24 926 patients with RHTN and 19 100 with stable controlled HTN. The PPV for RHTN in patients randomly selected for validation of the final algorithm was 99.1% (n = 113, CI: 95.2%-99.9%). The PPV for stable controlled HTN in patients randomly selected for validation of the final algorithm was 96.5% (n = 113, CI: 91.2%-99.0%). IRR analysis revealed a raw percent agreement of 91% (152/167) with Cohen's kappa statistic = 0.87.

CONCLUSIONS

We constructed and validated a RHTN computable phenotype algorithm and a stable controlled HTN computable phenotype algorithm. Both algorithms are based on the PCORnet CDM, allowing for future application to epidemiological and drug utilization based research.

Predictive Utility of Polygenic Risk Scores for Coronary Heart Disease in Three Major Racial and Ethnic Groups

Because polygenic risk scores (PRSs) for coronary heart disease (CHD) are derived from mainly European ancestry (EA) cohorts, their validity in African ancestry (AA) and Hispanic ethnicity (HE) individuals is unclear. We investigated associations of "restricted" and genome-wide PRSs with CHD in three major racial and ethnic groups in the U.S. The eMERGE cohort (mean age 48 ± 14 years, 58% female) included 45,645 EA, 7,597 AA, and 2,493 HE individuals. We assessed two restricted PRSs (PRS_Tikkanen and PRS_Tada; 28 and 50 variants, respectively) and two genome-wide PRSs (PRS_metaGRS and PRS_LDPred; 1.7 M and 6.6 M variants, respectively) derived from EA cohorts. Over a median follow-up of 11.1 years, 2,652 incident CHD events occurred. Hazard and odds ratios for the association of PRSs with CHD were similar in EA and HE cohorts but lower in AA cohorts. Genome-wide PRSs were more strongly associated with CHD than restricted PRSs were. PRS_metaGRS, the best performing PRS, was associated with CHD in all three cohorts; hazard ratios (95% CI) per 1 SD increase were 1.53 (1.46-1.60), 1.53 (1.23-1.90), and 1.27 (1.13-1.43) for incident CHD in EA, HE, and AA individuals, respectively. The hazard ratios were comparable in the EA and HE cohorts (p_interaction = 0.77) but were significantly attenuated in AA individuals (p_interaction= 2.9 × 10^-3). These results highlight the potential clinical utility of PRSs for CHD as well as the need to assemble diverse cohorts to generate ancestry- and ethnicity PRSs.

Electronic health records and polygenic risk scores for predicting disease risk

Accurate prediction of disease risk based on the genetic make-up of an individual is essential for effective prevention and personalized treatment. Nevertheless, to date, individual genetic variants from genome-wide association studies have achieved only moderate prediction of disease risk. The aggregation of genetic variants under a polygenic model shows promising improvements in prediction accuracies. Increasingly, electronic health records (EHRs) are being linked to patient genetic data in biobanks, which provides new opportunities for developing and applying polygenic risk scores in the clinic, to systematically examine and evaluate patient susceptibilities to disease. However, the heterogeneous nature of EHR data brings forth many practical challenges along every step of designing and implementing risk prediction strategies. In this Review, we present the unique considerations for using genotype and phenotype data from biobank-linked EHRs for polygenic risk prediction.

Optimizing Antihypertensive Medication Classification in Electronic Health Record-Based Data: Classification System Development and Methodological Comparison

Background

Computable phenotypes have the ability to utilize data within the electronic health record (EHR) to identify patients with certain characteristics. Many computable phenotypes rely on multiple types of data within the EHR including prescription drug information. Hypertension (HTN)-related computable phenotypes are particularly dependent on the correct classification of antihypertensive prescription drug information, as well as corresponding diagnoses and blood pressure information.

Objective

This study aimed to create an antihypertensive drug classification system to be utilized with EHR-based data as part of HTN-related computable phenotypes.

Methods

We compared 4 different antihypertensive drug classification systems based off of 4 different methodologies and terminologies, including 3 RxNorm Concept Unique Identifier (RxCUI)–based classifications and 1 medication name–based classification. The RxCUI-based classifications utilized data from (1) the Drug Ontology, (2) the new Medication Reference Terminology, and (3) the Anatomical Therapeutic Chemical Classification System and DrugBank, whereas the medication name–based classification relied on antihypertensive drug names. Each classification system was applied to EHR-based prescription drug data from hypertensive patients in the OneFlorida Data Trust.

Results

There were 13,627 unique RxCUIs and 8025 unique medication names from the 13,879,046 prescriptions. We observed a broad overlap between the 4 methods, with 84.1% (691/822) to 95.3% (695/729) of terms overlapping pairwise between the different classification methods. Key differences arose from drug products with multiple dosage forms, drug products with an indication of benign prostatic hyperplasia, drug products that contain more than 1 ingredient (combination products), and terms within the classification systems corresponding to retired or obsolete RxCUIs.

Conclusions

In total, 2 antihypertensive drug classifications were constructed, one based on RxCUIs and one based on medication name, that can be used in future computable phenotypes that require antihypertensive drug classifications.

Computational and Pharmacogenomic Insights on Hypertension Treatment: Rational Drug Design and Optimization Strategies

Background::

Hypertension is a prevalent cardiovascular complication caused by genetic and nongenetic factors. Blood pressure (BP) management is difficult because most patients become resistant to monotherapy soon after treatment initiation. Although many antihypertensive drugs are available, some patients do not respond to multiple drugs. Identification of personalized antihypertensive treatments is a key for better BP management.

Objective::

This review aimed to elucidate aspects of rational drug design and other methods to develop better hypertension management.

Results::

Among hypertension-related signaling mechanisms, the renin-angiotensin-aldosterone system is the leading genetic target for hypertension treatment. Identifying a single drug that acts on multiple targets is an emerging strategy for hypertension treatment, and could be achieved by discovering new drug targets with less mutated and highly conserved regions. Extending pharmacogenomics research to include patients with hypertension receiving multiple antihypertensive drugs could help identify the genetic markers of hypertension. However, available evidence on the role of pharmacogenomics in hypertension is limited and primarily focused on candidate genes. Studies on hypertension pharmacogenomics aim to identify the genetic causes of response variations to antihypertensive drugs. Genetic association studies have identified single nucleotide polymorphisms affecting drug responses. To understand how genetic traits alter drug responses, computational screening of mutagenesis can be utilized to observe drug response variations at the protein level, which can help identify new inhibitors and drug targets to manage hypertension.

Conclusions::

Rational drug design facilitates the discovery and design of potent inhibitors. However, further research and clinical validation are required before novel inhibitors can be clinically used as antihypertensive therapies.

Genome-Wide Association Study of Apparent Treatment-Resistant Hypertension in the CHARGE Consortium: The CHARGE Pharmacogenetics Working Group

BACKGROUND

Only a handful of genetic discovery efforts in apparent treatment-resistant hypertension (aTRH) have been described.

METHODS

We conducted a case-control genome-wide association study of aTRH among persons treated for hypertension, using data from 10 cohorts of European ancestry (EA) and 5 cohorts of African ancestry (AA). Cases were treated with 3 different antihypertensive medication classes and had blood pressure (BP) above goal (systolic BP ≥ 140 mm Hg and/or diastolic BP ≥ 90 mm Hg) or 4 or more medication classes regardless of BP control (nEA = 931, nAA = 228). Both a normotensive control group and a treatment-responsive control group were considered in separate analyses. Normotensive controls were untreated (nEA = 14,210, nAA = 2,480) and had systolic BP/diastolic BP < 140/90 mm Hg. Treatment-responsive controls (nEA = 5,266, nAA = 1,817) had BP at goal (<140/90 mm Hg), while treated with one antihypertensive medication class. Individual cohorts used logistic regression with adjustment for age, sex, study site, and principal components for ancestry to examine the association of single-nucleotide polymorphisms with case-control status. Inverse variance-weighted fixed-effects meta-analyses were carried out using METAL.

RESULTS

The known hypertension locus, CASZ1, was a top finding among EAs (P = 1.1 × 10-8) and in the race-combined analysis (P = 1.5 × 10-9) using the normotensive control group (rs12046278, odds ratio = 0.71 (95% confidence interval: 0.6-0.8)). Single-nucleotide polymorphisms in this locus were robustly replicated in the Million Veterans Program (MVP) study in consideration of a treatment-responsive control group. There were no statistically significant findings for the discovery analyses including treatment-responsive controls.

CONCLUSION

This genomic discovery effort for aTRH identified CASZ1 as an aTRH risk locus.

Renin-Angiotensin-Aldosterone System Modulates Blood Pressure Response During Vascular Endothelial Growth Factor Receptor Inhibition

Objectives

This study postulated that antihypertensive therapy with renin-angiotensin-aldosterone system (RAAS) inhibition may mitigate vascular endothelial growth factor inhibitor (VEGFi)–mediated increases in blood pressure more effectively than other antihypertensive medications in patients receiving VEGFi therapy.

Background

VEGFi therapy is commonly used in the treatment of cancer. One common side effect of VEGFi therapy is elevated blood pressure. Evidence suggests that the RAAS may be involved in VEGFi-mediated increases in blood pressure.

Methods

This retrospective cohort analysis was performed using a de-identified version of the electronic health record at Vanderbilt University Medical Center in Nashville, Tennessee. Subjects with cancer who were exposed to VEGFi therapy were identified, and blood pressure and medication data were extracted. Changes in mean systolic and diastolic blood pressure in response to VEGFi therapy in patients receiving RAAS inhibitor (RAASi) therapy before VEGFi initiation were compared with changes in mean systolic and diastolic blood pressure in patients not receiving RAASi therapy before VEGFi initiation.

Results

Mean systolic and diastolic blood pressure rose in both groups after VEGFi use; however, patients who had RAASi therapy before VEGFi initiation had a significantly lower increase in systolic blood pressure as compared with patients with no RAASi therapy (2.46 mm Hg [95% confidence interval: 0.7 to 4.2] compared with 4.56 mm Hg [95% confidence interval: 3.5 to 5.6], respectively; p = 0.034).

Conclusions

In a real-world clinical population, RAASi therapy before VEGFi initiation may ameliorate VEGFi-mediated increases in blood pressure. Randomized clinical trials are needed to further our understanding of the role of RAASi therapy in VEGFi-mediated increases in blood pressure.

Generality of genomic findings on blood pressure traits and its usefulness in precision medicine in diverse populations: A systematic review

Remarkable findings from genome-wide association studies (GWAS) on blood pressure (BP) traits have made new insights for developing precision medicine toward more effective screening measures. However, generality of GWAS findings in diverse populations is hampered by some technical limitations. There is no comprehensive study to evaluate source(s) of the non-generality of GWAS results on BP traits, so to fill the gap, this systematic review study was carried out. Using MeSH terms, 1545 records were detected through searching in five databases and 49 relevant full-text articles were included in our review. Overall, 749 unique variants were reported, of those, majority of variants have been detected in Europeans and were associated to systolic and diastolic BP traits. Frequency of genetic variants with same position was low in European and non-European populations (n = 38). However, more than 200 (>25%) single nucleotide polymorphisms were found on same loci or linkage disequilibrium blocks (r² ≥ 80%). Investigating for locus position and linkage disequilibrium of infrequent unique variants showed modest to high reproducibility of findings in Europeans that in some extent was generalizable in other populations. Beyond theoretical limitations, our study addressed other possible sources of non-generality of GWAS findings for BP traits in the same and different origins.

Using Electronic Health Records To Generate Phenotypes For Research

Electronic health records contain patient-level data collected during and for clinical care. Data within the electronic health record include diagnostic billing codes, procedure codes, vital signs, laboratory test results, clinical imaging, and physician notes. With repeated clinic visits, these data are longitudinal, providing important information on disease development, progression, and response to treatment or intervention strategies. The near universal adoption of electronic health records nationally has the potential to provide population-scale real-world clinical data accessible for biomedical research, including genetic association studies. For this research potential to be realized, high-quality research-grade variables must be extracted from these clinical data warehouses. We describe here common and emerging electronic phenotyping approaches applied to electronic health records, as well as current limitations of both the approaches and the biases associated with these clinically collected data that impact their use in research. © 2018 by John Wiley & Sons, Inc.

A case study evaluating the portability of an executable computable phenotype algorithm across multiple institutions and electronic health record environments

Electronic health record (EHR) algorithms for defining patient cohorts are commonly shared as free-text descriptions that require human intervention both to interpret and implement. We developed the Phenotype Execution and Modeling Architecture (PhEMA, http://projectphema.org) to author and execute standardized computable phenotype algorithms. With PhEMA, we converted an algorithm for benign prostatic hyperplasia, developed for the electronic Medical Records and Genomics network (eMERGE), into a standards-based computable format. Eight sites (7 within eMERGE) received the computable algorithm, and 6 successfully executed it against local data warehouses and/or i2b2 instances. Blinded random chart review of cases selected by the computable algorithm shows PPV ≥90%, and 3 out of 5 sites had >90% overlap of selected cases when comparing the computable algorithm to their original eMERGE implementation. This case study demonstrates potential use of PhEMA computable representations to automate phenotyping across different EHR systems, but also highlights some ongoing challenges.

The eMERGE genotype set of 83,717 subjects imputed to ~40 million variants genome wide and association with the herpes zoster medical record phenotype

The Electronic Medical Records and Genomics (eMERGE) network is a network of medical centers with electronic medical records linked to existing biorepository samples for genomic discovery and genomic medicine research. The network sought to unify the genetic results from 78 Illumina and Affymetrix genotype array batches from 12 contributing medical centers for joint association analysis of 83,717 human participants. In this report, we describe the imputation of eMERGE results and methods to create the unified imputed merged set of genome‐wide variant genotype data. We imputed the data using the Michigan Imputation Server, which provides a missing single‐nucleotide variant genotype imputation service using the minimac3 imputation algorithm with the Haplotype Reference Consortium genotype reference set. We describe the quality control and filtering steps used in the generation of this data set and suggest generalizable quality thresholds for imputation and phenotype association studies. To test the merged imputed genotype set, we replicated a previously reported chromosome 6 HLA‐B herpes zoster (shingles) association and discovered a novel zoster‐associated loci in an epigenetic binding site near the terminus of chromosome 3 (3p29).

Genome-wide association analysis of common genetic variants of resistant hypertension

Resistant hypertension (RHTN), defined as uncontrolled blood pressure (BP) ≥ 140/90 using three or more drugs or controlled BP (<140/90) using four or more drugs, is associated with adverse outcomes, including decline in kidney function. We conducted a genome-wide association analysis in 1194 White and Hispanic participants with hypertension and coronary artery disease from the INternational VErapamil-SR Trandolapril STudy-GENEtic Substudy (INVEST-GENES). Top variants associated with RHTN at p < 10^-4 were tested for replication in 585 White and Hispanic participants with hypertension and subcortical strokes from the Secondary Prevention of Subcortical Strokes GENEtic Substudy (SPS3-GENES). A genetic risk score for RHTN was created by summing the risk alleles of replicated RHTN signals. rs11749255 in MSX2 was associated with RHTN in INVEST (odds ratio (OR) (95% CI) = 1.50 (1.2-1.8), p = 7.3 × 10^-5) and replicated in SPS3 (OR = 2.0 (1.4-2.8), p = 4.3 × 10^-5), with genome-wide significance in meta-analysis (OR = 1.60 (1.3-1.9), p = 3.8 × 10^-8). Other replicated signals were in IFLTD1 and PTPRD. IFLTD1 rs6487504 was associated with RHTN in INVEST (OR = 1.90 (1.4-2.5), p = 1.1 × 10^-5) and SPS3 (OR = 1.70 (1.2-2.5), p = 4 × 10^-3). PTPRD rs324498, a previously reported RHTN signal, was among the top signals in INVEST (OR = 1.60 (1.3-2.0), p = 3.4 × 10^-5) and replicated in SPS3 (OR = 1.60 (1.1-2.4), one-sided p = 0.005). Participants with the highest number of risk alleles were at increased risk of RHTN compared to participants with a lower number (p-trend = 1.8 × 10^-15). Overall, we identified and replicated associations with RHTN in the MSX2, IFLTD1, and PTPRD regions, and combined these associations to create a genetic risk score.

Large-Scale Analysis of Genetic and Clinical Patient Data

Biomedical data science has experienced an explosion of new data over the past decade. Abundant genetic and genomic data are increasingly available in large, diverse data sets due to the maturation of modern molecular technologies. Along with these molecular data, dense, rich phenotypic data are also available on comprehensive clinical data sets from health care provider organizations, clinical trials, population health registries, and epidemiologic studies. The methods and approaches for interrogating these large genetic/genomic and clinical data sets continue to evolve rapidly, as our understanding of the questions and challenges continue to emerge. In this review, the state-of-the-art methodologies for genetic/genomic analysis along with complex phenomics will be discussed. This field is changing and adapting to the novel data types made available, as well as technological advances in computation and machine learning. Thus, I will also discuss the future challenges in this exciting and innovative space. The promises of precision medicine rely heavily on the ability to marry complex genetic/genomic data with clinical phenotypes in meaningful ways.

Characteristics and treatment of African-American and European-American patients with resistant hypertension identified using the electronic health record in an academic health centre: a case-control study

OBJECTIVE

To identify patients with hypertension with resistant and controlled blood pressure (BP) using electronic health records (EHRs) in order to elucidate practices in the real-world clinical treatment of hypertension and to enable future genetic studies.

DESIGN

Using EHRs, we developed and validated algorithms to identify patients with resistant and controlled hypertension.

SETTING

An academic medical centre in Nashville, Tennessee.

POPULATION

European-American (EA) and African-American (AA) patients with hypertension.

MAIN OUTCOME MEASURES

Demographic characteristics: race, age, gender, body mass index, outpatient BPs and the history of diabetes mellitus, chronic kidney disease stage 3, ischaemic heart disease, transient ischaemic attack, atrial fibrillation and sleep apnoea.

MEDICATION TREATMENT

All antihypertensive medication classes prescribed to a patient at the time of classification and ever prescribed following classification.

RESULTS

The algorithms had performance metrics exceeding 92%. The prevalence of resistant hypertension in the total hypertensive population was 7.3% in EA and 10.5% in AA. At diagnosis, AA were younger, heavier, more often female and had a higher incidence of type 2 diabetes and higher BPs than EA. AA with resistant hypertension were more likely to be treated with vasodilators, dihydropyridine calcium channel blockers and alpha-2 agonists while EA were more likely to be treated with angiotensin receptor blockers, renin inhibitors and beta blockers. Mineralocorticoid receptor antagonists use was increased in patients treated with more than four antihypertensive medications compared with patients treated with three (12.4% vs 2.6% in EA, p<0.001; 12.3% vs 2.8% in AA, p<0.001). The number of patients treated with a mineralocorticoid receptor antagonist increased to 37.4% in EA and 41.2% in AA over a mean follow-up period of 7.4 and 8.7 years, respectively.

CONCLUSIONS

Clinical treatment of resistant hypertension differs in EA and AA patients. These results demonstrate the feasibility of identifying resistant hypertension using an EHR.

Genetics of Resistant Hypertension: the Missing Heritability and Opportunities

PURPOSE OF THE REVIEW

Blood pressure regulation in humans has long been known to be a genetically determined trait. The identification of causal genetic modulators for this trait has been unfulfilling at the least. Despite the recent advances of genome-wide genetic studies, loci associated with hypertension or blood pressure still explain a very low percentage of the overall variation of blood pressure in the general population. This has precluded the translation of discoveries in the genetics of human hypertension to clinical use. Here, we propose the combined use of resistant hypertension as a trait for mapping genetic determinants in humans and the integration of new large-scale technologies to approach in model systems the multidimensional nature of the problem.

RECENT FINDINGS

New large-scale efforts in the genetic and genomic arenas are paving the way for an increased and granular understanding of genetic determinants of hypertension. New technologies for whole genome sequence and large-scale forward genetic screens can help prioritize gene and gene-pathways for downstream characterization and large-scale population studies, and guided pharmacological design can be used to drive discoveries to the translational application through better risk stratification and new therapeutic approaches. Although significant challenges remain in the mapping and identification of genetic determinants of hypertension, new large-scale technological approaches have been proposed to surpass some of the shortcomings that have limited progress in the area for the last three decades. The incorporation of these technologies to hypertension research may significantly help in the understanding of inter-individual blood pressure variation and the deployment of new phenotyping and treatment approaches for the condition.

Genomic medicine for kidney disease

Technologies such as next-generation sequencing and chromosomal microarray have advanced the understanding of the molecular pathogenesis of a variety of renal disorders. Genetic findings are increasingly used to inform the clinical management of many nephropathies, enabling targeted disease surveillance, choice of therapy, and family counselling. Genetic analysis has excellent diagnostic utility in paediatric nephrology, as illustrated by sequencing studies of patients with congenital anomalies of the kidney and urinary tract and steroid-resistant nephrotic syndrome. Although additional investigation is needed, pilot studies suggest that genetic testing can also provide similar diagnostic insight among adult patients. Reaching a genetic diagnosis first involves choosing the appropriate testing modality, as guided by the clinical presentation of the patient and the number of potential genes associated with the suspected nephropathy. Genome-wide sequencing increases diagnostic sensitivity relative to targeted panels, but holds the challenges of identifying causal variants in the vast amount of data generated and interpreting secondary findings. In order to realize the promise of genomic medicine for kidney disease, many technical, logistical, and ethical questions that accompany the implementation of genetic testing in nephrology must be addressed. The creation of evidence-based guidelines for the utilization and implementation of genetic testing in nephrology will help to translate genetic knowledge into improved clinical outcomes for patients with kidney disease.

Genetic Variants Associated With Uncontrolled Blood Pressure on Thiazide Diuretic/β-Blocker Combination Therapy in the PEAR (Pharmacogenomic Evaluation of Antihypertensive Responses) and INVEST (International Verapamil-SR Trandolapril Study) Trials

BACKGROUND

The majority of hypertensive individuals require combination antihypertensive therapy to achieve adequate blood pressure (BP) control. This study aimed to identify genetic variants associated with uncontrolled BP on combination therapy with a thiazide diuretic and a β-blocker.

METHODS AND RESULTS

A genome-wide association study of uncontrolled BP on combination therapy was conducted among 314 white participants of the PEAR (Pharmacogenomic Evaluation of Antihypertensive Responses) trial. Multivariable logistic regression analysis was used. Genetic variants meeting a suggestive level of significance (P<1.0E-05) were tested for replication in an external cohort, INVEST (International Verapamil-SR Trandolapril study). We also examined genome-wide variant associations with systolic and diastolic BP response on combination therapy and tested for replication. We discovered a single nucleotide polymorphism, the rs261316 major allele, at chromosome 15 in the gene ALDH1A2 associated with an increased odds of having uncontrolled BP on combination therapy (odds ratio: 2.56, 95% confidence interval, 1.69-3.88, P=8.64E-06). This single nucleotide polymorphism replicated (odds ratio: 1.86, 95% confidence interval, 1.35-2.57, P=0.001) and approached genome-wide significance in the meta-analysis between discovery and replication cohorts (odds ratio: 2.16, 95% confidence interval, 1.63-2.86, P=8.60E-08). Other genes in the region surrounding rs261316 (ALDH1A2) include AQP9 and LIPC.

CONCLUSIONS

A single nucleotide polymorphism in the gene ALDH1A2 may be associated with uncontrolled BP following treatment with a thiazide diuretic/β-blocker combination.

CLINICAL TRIAL REGISTRATION

URL: https://www.clinicaltrials.gov. Unique identifier: NCT00246519.