Johnson, Julie A.

Objective Genome‐wide association studies (GWAS) have been successful at identifying associations with stroke and stroke subtypes, but have not yet identified any associations solely with small vessel stroke (SVS). SVS comprises one quarter of all ischemic stroke and is a major manifestation of cerebral small vessel disease, the primary cause of vascular cognitive impairment. Studies across neurological traits have shown that younger‐onset cases have an increased genetic burden. We leveraged this increased genetic burden by performing an age‐at‐onset informed GWAS meta‐analysis, including a large younger‐onset SVS population, to identify novel associations with stroke. Methods We used a three‐stage age‐at‐onset informed GWAS to identify novel genetic variants associated with stroke. On identifying a novel locus associated with SVS, we assessed its influence on other small vessel disease phenotypes, as well as on messenger RNA (mRNA) expression of nearby genes, and on DNA methylation of nearby CpG sites in whole blood and in the fetal brain. Results We identified an association with SVS in 4,203 cases and 50,728 controls on chromosome 16q24.2 (odds ratio [OR; 95% confidence interval {CI}] = 1.16 [1.10–1.22]; p = 3.2 × 10 −9 ). The lead single‐nucleotide polymorphism (rs12445022) was also associated with cerebral white matter hyperintensities (OR [95% CI] = 1.10 [1.05–1.16]; p = 5.3 × 10 −5 ; N = 3,670), but not intracerebral hemorrhage (OR [95% CI] = 0.97 [0.84–1.12]; p = 0.71; 1,545 cases, 1,481 controls). rs12445022 is associated with mRNA expression of ZCCHC14 in arterial tissues ( p = 9.4 × 10 −7 ) and DNA methylation at probe cg16596957 in whole blood ( p = 5.3 × 10 −6 ). Interpretation 16q24.2 is associated with SVS. Associations of the locus with expression of ZCCHC14 and DNA methylation suggest the locus acts through changes to regulatory elements. Ann Neurol 2017;81:383–394

Lower melatonin level, melatonin receptor gene variations, and atenolol treatment are associated with glucose dysregulation. We investigated whether atenolol‐related glucose and melatonin changes are correlated, and whether single nucleotide polymorphisms (SNPs) in melatonin candidate genes contribute to interindividual variation in glucose change. Hypertensive Caucasians (n = 232) from the Pharmacogenomic Evaluation of Antihypertensive Responses (PEAR) study treated with atenolol for 9 weeks were studied. Urinary 6‐sulfatoxymelatonin (aMT6s) was measured pre‐ and posttreatment and normalized to urinary creatinine. Pharmacogenetic effects on glucose change of 160 SNPs in 16 melatonin candidate genes were assessed with multiple linear regression. Atenolol was associated with increased glucose (1.8 ± 10.1mg/dl, P = 0.02) and decreased aMT6s (–4.5 ± 10.1 ng/mg, P < 0.0001). However, the aMT6s change was not correlated with post‐atenolol glucose change. SNP rs11649514 in PRKCB was associated with glucose change (P = 1.0×10−4). PRKCB is involved in the melatonin‐insulin regulatory pathway, and may be important in mediating clinically meaningful atenolol‐related hyperglycemia.

Background Identification of genetic markers of antihypertensive drug responses could assist in individualization of hypertension treatment. Methods and Results We conducted a genome‐wide association study to identify gene loci influencing the responsiveness of 228 male patients to 4 classes of antihypertensive drugs. The Genetics of Drug Responsiveness in Essential Hypertension ( GENRES ) study is a double‐blind, placebo‐controlled cross‐over study where each subject received amlodipine, bisoprolol, hydrochlorothiazide, and losartan, each as a monotherapy, in a randomized order. Replication analyses were performed in 4 studies with patients of European ancestry ( PEAR Study, N=386; GERA I and II Studies, N=196 and N=198; SOPHIA Study, N=372). We identified 3 single‐nucleotide polymorphisms within the ACY 3 gene that showed associations with bisoprolol response reaching genome‐wide significance ( P <5×10 −8 ); however, this could not be replicated in the PEAR Study using atenolol. In addition, 39 single‐nucleotide polymorphisms showed P values of 10 −5 to 10 −7 . The 20 top‐associated single‐nucleotide polymorphisms were different for each antihypertensive drug. None of these top single‐nucleotide polymorphisms co‐localized with the panel of >40 genes identified in genome‐wide association studies of hypertension. Replication analyses of GENRES results provided suggestive evidence for a missense variant (rs3814995) in the NPHS 1 (nephrin) gene influencing losartan response, and for 2 variants influencing hydrochlorothiazide response, located within or close to the ALDH 1A3 (rs3825926) and CLIC 5 (rs321329) genes. Conclusions These data provide some evidence for a link between biology of the glomerular protein nephrin and antihypertensive action of angiotensin receptor antagonists and encourage additional studies on aldehyde dehydrogenase–mediated reactions in antihypertensive drug action.

Background Resistant hypertension (RHTN), defined by lack of blood pressure (BP) control despite treatment with at least 3 antihypertensive drugs, increases cardiovascular risk compared with controlled hypertension. Yet, there are few data on genetic variants associated with RHTN. Methods and Results We used a gene‐centric array containing ≈50 000 single‐nucleotide polymorphisms (SNPs) to identify polymorphisms associated with RHTN in hypertensive participants with coronary artery disease (CAD) from INVEST‐GENES (the INnternational VErapamil‐SR Trandolapril STudy—GENEtic Substudy). RHTN was defined as BP≥140/90 on 3 drugs, or any BP on 4 or more drugs. Logistic regression analysis was performed in European Americans (n=904) and Hispanics (n=837), using an additive model adjusted for age, gender, randomized treatment assignment, body mass index, principal components for ancestry, and other significant predictors of RHTN. Replication of the top SNP was conducted in 241 European American women from WISE (Women's Ischemia Syndrome Evaluation), where RHTN was defined similarly. To investigate the functional effect of rs12817819, mRNA expression was measured in whole blood. We found ATP2B1 rs12817819 associated with RHTN in both INVEST European Americans ( P ‐value=2.44×10 −3 , odds ratio=1.57 [1.17 to 2.01]) and INVEST Hispanics ( P =7.69×10 −4 , odds ratio=1.76 [1.27 to 2.44]). A consistent trend was observed at rs12817819 in WISE, and the INVEST‐WISE meta‐analysis result reached chip‐wide significance ( P =1.60×10 −6 , odds ratio=1.65 [1.36 to 1.95]). Expression analyses revealed significant differences in ATP2B1 expression by rs12817819 genotype. Conclusions The ATP2B1 rs12817819 A allele is associated with increased risk for RHTN in hypertensive participants with documented CAD or suspected ischemic heart disease. Clinical Trial Registration URL: www.clinicaltrials.gov ; Unique identifiers: NCT00133692 (INVEST), NCT00000554 (WISE).

Study ObjectiveTo develop and validate a predictive model for glucose change and risk for new‐onset impaired fasting glucose in hypertensive participants following treatment with atenolol or hydrochlorothiazide (HCTZ).DesignRandomized multicenter clinical trial.PatientsA total of 735 white or African‐American men and women with uncomplicated hypertension.Measurements and Main ResultsPharmacogenomic Evaluation of Antihypertensive Responses (PEAR) is a randomized clinical trial to assess the genetic and nongenetic predictors of blood pressure response and adverse metabolic effects following treatment with atenolol or HCTZ. To develop and validate predictive models for glucose change, PEAR participants were randomly divided into a derivation cohort of 367 and a validation cohort of 368. Linear and logistic regression modeling were used to build models of drug‐associated glucose change and impaired fasting glucose (IFG), respectively, in the derivation cohorts. These models were then evaluated in the validation cohorts. For glucose change after atenolol or HCTZ treatment, baseline glucose was a significant (p<0.0001) predictor, explaining 13% of the variability in glucose change after atenolol and 12% of the variability in glucose change after HCTZ. Baseline glucose was also the strongest and most consistent predictor (p<0.0001) for development of IFG after atenolol or HCTZ monotherapy. The area under the receiver operating curve was 0.77 for IFG after atenolol and 0.71 after HCTZ treatment, respectively.ConclusionBaseline glucose is the primary predictor of atenolol or HCTZ‐associated glucose increase and development of IFG after treatment with either drug.

AbstractObjectiveElevations in uric acid (UA) and the associated hyperuricaemia are commonly observed secondary to treatment with thiazide diuretics. We sought to identify novel single nucleotide polymorphisms (SNPs) associated with hydrochlorothiazide (HCTZ)‐induced elevations in UA and hyperuricaemia.MethodsA genome‐wide association study of HCTZ‐induced changes in UA was performed in Caucasian and African American participants from the pharmacogenomic evaluation of antihypertensive responses (PEAR) study who were treated with HCTZ monotherapy. Suggestive SNPs were replicated in Caucasians and African Americans from the PEAR study who were treated with HCTZ add‐on therapy. Replicated regions were followed up through expression and pathway analysis.ResultsFive unique gene regions were identified in African Americans (LUC7L2, ANKRD17/COX18, FTO, PADI4 and PARD3B), and one region was identified in Caucasians (GRIN3A). Increases in UA of up to 1.8 mg dL−1 were observed following HCTZ therapy in individuals homozygous for risk alleles, with heterozygotes displaying an intermediate phenotype. Several risk alleles were also associated with an increased risk of HCTZ‐induced clinical hyperuricaemia. A composite risk score, constructed in African Americans using the ‘top’ SNP from each gene region, was strongly associated with HCTZ‐induced UA elevations (P = 1.79 × 10−7) and explained 11% of the variability in UA response. Expression studies in RNA from whole blood revealed significant differences in expression of FTO by rs4784333 genotype. Pathway analysis showed putative connections between many of the genes identified through common microRNAs.ConclusionSeveral novel gene regions were associated with HCTZ‐induced UA elevations in African Americans (LUC7L2, COX18/ANKRD17, FTO, PADI4 and PARD3B), and one region was associated with these elevations in Caucasians (GRIN3A).

Study ObjectiveTo assess the relationship of the 33 single nucleotide polymorphisms (SNPs) previously associated with fasting glucose in Caucasians in genome‐wide association studies (GWAS) with glucose response to antihypertensive drugs shown to increase risk for hyperglycemia and diabetes.DesignRandomized, multicenter clinical trial.PatientsA total of 456 Caucasian men and women with uncomplicated hypertension.Measurements and Main ResultsThe Pharmacogenomic Evaluation of Antihypertensives Responses study evaluated blood pressure and glucose response in uncomplicated hypertensive patients randomized to either atenolol or hydrochlorothiazide (HCTZ) monotherapy, followed by combination therapy with both agents. Association of these SNPs with atenolol‐ or HCTZ‐induced glucose response was evaluated in 456 Caucasian patients using linear regression adjusting for age, sex, body mass index, baseline glucose, baseline insulin, and principal component for ancestry.The SNP rs340874 in the 5′ region of PROX1 gene was significantly associated with atenolol‐induced glucose change (p=0.0013). Participants harboring the C allele of this SNP had greater glucose elevation after approximately 9 weeks of atenolol monotherapy (β = +2.39 mg/dl per C allele), consistent with the direction of effect in fasting glucose GWAS, that showed the C allele is associated with higher fasting glucose.ConclusionThese data suggest that PROX1 SNP rs340874, discovered in fasting glucose GWAS, may also be a pharmacogenetic risk factor for antihypertensive‐induced hyperglycemia. β‐blockers and thiazides may interact with genetic risk factors to increase risk for dysglycemia and diabetes.

Background and Purpose— Meta-analyses of extant genome-wide data illustrate the need to focus on subtypes of ischemic stroke for gene discovery. The National Institute of Neurological Disorders and Stroke SiGN (Stroke Genetics Network) contributes substantially to meta-analyses that focus on specific subtypes of stroke. Methods— The National Institute of Neurological Disorders and Stroke SiGN includes ischemic stroke cases from 24 genetic research centers: 13 from the United States and 11 from Europe. Investigators harmonize ischemic stroke phenotyping using the Web-based causative classification of stroke system, with data entered by trained and certified adjudicators at participating genetic research centers. Through the Center for Inherited Diseases Research, the Network plans to genotype 10 296 carefully phenotyped stroke cases using genome-wide single nucleotide polymorphism arrays and adds to these another 4253 previously genotyped cases, for a total of 14 549 cases. To maximize power for subtype analyses, the study allocates genotyping resources almost exclusively to cases. Publicly available studies provide most of the control genotypes. Center for Inherited Diseases Research–generated genotypes and corresponding phenotypes will be shared with the scientific community through the US National Center for Biotechnology Information database of Genotypes and Phenotypes, and brain MRI studies will be centrally archived. Conclusions— The Stroke Genetics Network, with its emphasis on careful and standardized phenotyping of ischemic stroke and stroke subtypes, provides an unprecedented opportunity to uncover genetic determinants of ischemic stroke.

To identify novel genes influencing blood pressure response to thiazide diuretic therapy for hypertension, we conducted genome-wide association meta-analyses of ≈1.1 million single-nucleotide polymorphisms in a combined sample of 424 European Americans with primary hypertension treated with hydrochlorothiazide from the Pharmacogenomic Evaluation of Antihypertensive Responses study (n=228) and the Genetic Epidemiology of Responses to Antihypertensive study (n=196). Polymorphisms associated with blood pressure response at P <10 –5 were tested for replication of the associations in independent samples of hydrochlorothiazide-treated European hypertensives. The rs16960228 polymorphism in protein kinase C, α replicated for same-direction association with diastolic blood pressure response in the Nordic Diltiazem study (n=420) and the Genetics of Drug Responsiveness in Essential Hypertension study (n=206), and the combined 4-study meta-analysis P value achieved genome-wide significance ( P =3.3×10 −8 ). Systolic or diastolic blood pressure responses were consistently greater in carriers of the rs16960228 A allele than in GG homozygotes (>4/4 mm Hg) across study samples. The rs2273359 polymorphism in the GNAS-EDN3 region also replicated for same-direction association with systolic blood pressure response in the Nordic Diltiazem study, and the combined 3-study meta-analysis P value approached genome-wide significance ( P =5.5×10 −8 ). The findings document clinically important effects of genetic variation at novel loci on blood pressure response to a thiazide diuretic, which may be a basis for individualization of antihypertensive drug therapy and identification of new drug targets.