Boerwinkle, Eric

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.

Background The relationship between mitochondrial DNA copy number (mtDNA CN) and cardiovascular disease remains elusive. Methods and Results We performed cross‐sectional and prospective association analyses of blood‐derived mtDNA CN and cardiovascular disease outcomes in 27 316 participants in 8 cohorts of multiple racial and ethnic groups with whole‐genome sequencing. We also performed Mendelian randomization to explore causal relationships of mtDNA CN with coronary heart disease (CHD) and cardiometabolic risk factors (obesity, diabetes, hypertension, and hyperlipidemia). P <0.01 was used for significance. We validated most of the previously reported associations between mtDNA CN and cardiovascular disease outcomes. For example, 1‐SD unit lower level of mtDNA CN was associated with 1.08 (95% CI, 1.04–1.12; P <0.001) times the hazard for developing incident CHD, adjusting for covariates. Mendelian randomization analyses showed no causal effect from a lower level of mtDNA CN to a higher CHD risk (β=0.091; P =0.11) or in the reverse direction (β=−0.012; P =0.076). Additional bidirectional Mendelian randomization analyses revealed that low‐density lipoprotein cholesterol had a causal effect on mtDNA CN (β=−0.084; P <0.001), but the reverse direction was not significant ( P =0.059). No causal associations were observed between mtDNA CN and obesity, diabetes, and hypertension, in either direction. Multivariable Mendelian randomization analyses showed no causal effect of CHD on mtDNA CN, controlling for low‐density lipoprotein cholesterol level ( P =0.52), whereas there was a strong direct causal effect of higher low‐density lipoprotein cholesterol on lower mtDNA CN, adjusting for CHD status (β=−0.092; P <0.001). Conclusions Our findings indicate that high low‐density lipoprotein cholesterol may underlie the complex relationships between mtDNA CN and vascular atherosclerosis.

Abstract Exonic variants present some of the strongest links between genotype and phenotype. However, these variants can have significant inter-individual pathogenicity differences, known as variable penetrance. In this study, we propose a model where genetically controlled mRNA splicing modulates the pathogenicity of exonic variants. By first cataloging exonic inclusion from RNA-sequencing data in GTEx V8, we find that pathogenic alleles are depleted on highly included exons. Using a large-scale phased whole genome sequencing data from the TOPMed consortium, we observe that this effect may be driven by common splice-regulatory genetic variants, and that natural selection acts on haplotype configurations that reduce the transcript inclusion of putatively pathogenic variants, especially when limiting to haploinsufficient genes. Finally, we test if this effect may be relevant for autism risk using families from the Simons Simplex Collection, but find that splicing of pathogenic alleles has a penetrance reducing effect here as well. Overall, our results indicate that common splice-regulatory variants may play a role in reducing the damaging effects of rare exonic variants.

Nononcogenic somatic mutations are thought to be uncommon and inconsequential. To test this, we analyzed 43,693 National Heart, Lung and Blood Institute Trans-Omics for Precision Medicine blood whole genomes from 37 cohorts and identified 7131 non-missense somatic mutations that are recurrently mutated in at least 50 individuals. These recurrent non-missense somatic mutations (RNMSMs) are not clearly explained by other clonal phenomena such as clonal hematopoiesis. RNMSM prevalence increased with age, with an average 50-year-old having 27 RNMSMs. Inherited germline variation associated with RNMSM acquisition. These variants were found in genes involved in adaptive immune function, proinflammatory cytokine production, and lymphoid lineage commitment. In addition, the presence of eight specific RNMSMs associated with blood cell traits at effect sizes comparable to Mendelian genetic mutations. Overall, we found that somatic mutations in blood are an unexpectedly common phenomenon with ancestry-specific determinants and human health consequences.

Background: Risk for venous thromboembolism has a strong genetic component. Whole genome sequencing from the TOPMed program (Trans-Omics for Precision Medicine) allowed us to look for new associations, particularly rare variants missed by standard genome-wide association studies. Methods: The 3793 cases and 7834 controls (11.6% of cases were individuals of African, Hispanic/Latino, or Asian ancestry) were analyzed using a single variant approach and an aggregate gene-based approach using our primary filter (included only loss-of-function and missense variants predicted to be deleterious) and our secondary filter (included all missense variants). Results: Single variant analyses identified associations at 5 known loci. Aggregate gene-based analyses identified only PROC (odds ratio, 6.2 for carriers of rare variants; P =7.4×10 −14 ) when using our primary filter. Employing our secondary variant filter led to a smaller effect size at PROC (odds ratio, 3.8; P =1.6×10 −14 ), while excluding variants found only in rare isoforms led to a larger one (odds ratio, 7.5). Different filtering strategies improved the signal for 2 other known genes: PROS1 became significant (minimum P =1.8×10 −6 with the secondary filter), while SERPINC1 did not (minimum P =4.4×10 −5 with minor allele frequency <0.0005). Results were largely the same when restricting the analyses to include only unprovoked cases; however, one novel gene, MS4A1 , became significant ( P =4.4×10 −7 using all missense variants with minor allele frequency <0.0005). Conclusions: Here, we have demonstrated the importance of using multiple variant filtering strategies, as we detected additional genes when filtering variants based on their predicted deleteriousness, frequency, and presence on the most expressed isoforms. Our primary analyses did not identify new candidate loci; thus larger follow-up studies are needed to replicate the novel MS4A1 locus and to identify additional rare variation associated with venous thromboembolism.

Background: Underrepresented persons are often not included in biomedical research. It is unknown if the general Asian American population is being represented in All of Us . The purpose of this study was to compare the Asian demographic data in the All of Us cohort with the Asian nationally representative data from the American Community Survey. Method: Demographic characteristics and health literacy of Asians in All of Us were examined. Findings were qualitatively compared with the Asian data in the 2019 American Community Survey 1-year estimate. Results: Compared with the national composition of Asians, less All of Us participants were born outside the United States (64% vs 79%), were younger, and had higher levels of education (76% vs 52%). Over 60% of All of Us participants reported high levels of health literacy. Conclusion: This study had implications for the development of strategies that ensure diverse populations are represented in biomedical research.

Abstract With aging, somatic mutations in hematopoietic stem and progenitor cells (HSPC) can give rise to clonal hematopoiesis of indeterminate potential (CHIP), a premalignant state associated with diverse age-related diseases. Here we report the largest multi-ancestry genome-wide analysis of CHIP to date (N = 1,018,305), including individuals of African (N = 85,978), Admixed American (N = 191,371), East Asian (N = 13,532), European (N = 694,015), and South Asian (N = 13,193) ancestry. Multi-ancestry meta-analyses identified 72 genome-wide significant loci, including 44 novel associations implicating genes such as AFF1 , ATF7IP , ATP8B4 , BCL2 , CEBPA , CYRIA , DNM2 , ELF1 , NKX2-3 , PIK3CB , PRDM16 , RPN1 , TERC , and TRIM4 . Notably, variants at MECOM and PHF20L1 showed opposite allelic effects between DNMT3A - and non- DNMT3A -driven CHIP, highlighting driver-specific germline influences. Implicated loci converge on pathways regulating telomere maintenance, cell-cycle control, hematopoietic transcription, DNA damage response and immune signaling. These findings support a model in which germline variation both expands the HSPC pool and biases clonal selection in a driver-dependent manner, providing a mechanistic basis for inter-individual heterogeneity in CHIP. Phenome-wide association analyses further linked CHIP to hematologic, neoplastic, and circulatory traits, with enrichment across hematopoietic and non-hematopoietic cell types. Together, this work expands the genomic and phenomic landscape of CHIP and reveals germline–somatic interactions that shape clonal evolution during aging.

Human genetic studies support an inverse causal relationship between leukocyte telomere length (LTL) and coronary artery disease (CAD), but directionally mixed effects for LTL and diverse malignancies. Clonal hematopoiesis of indeterminate potential (CHIP), characterized by expansion of hematopoietic cells bearing leukemogenic mutations, predisposes both hematologic malignancy and CAD. TERT (which encodes telomerase reverse transcriptase) is the most significantly associated germline locus for CHIP in genome-wide association studies. Here, we investigated the relationship between CHIP, LTL, and CAD in the Trans-Omics for Precision Medicine (TOPMed) program ( n = 63,302) and UK Biobank ( n = 47,080). Bidirectional Mendelian randomization studies were consistent with longer genetically imputed LTL increasing propensity to develop CHIP, but CHIP then, in turn, hastens to shorten measured LTL (mLTL). We also demonstrated evidence of modest mediation between CHIP and CAD by mLTL. Our data promote an understanding of potential causal relationships across CHIP and LTL toward prevention of CAD.

37 Knowledge of genetic factors that influence pharmacodynamics of responses to antihypertensive drugs may provide insights into molecular mechanisms contributing to elevation of BP and response to treatment. The T-allele of the C825T polymorphism of the gene coding for the G-protein β 3 subunit ( GNB3 ) has been associated with increased Na + -H + exchange and low plasma renin in individuals with EHT. Moreover, diminished counterregulatory activity of the renin-angiotensin system has been associated with greater BP declines in response to negative sodium balance induced by diuretic therapy. To assess the possible influence of GNB3 genotype on BP response to diuretic therapy, we measured the C825T polymorphism in community-based samples of 204 blacks (139 women, 65 men) and 192 whites (77 women, 115 men) with EHT (mean age=48±7 years) who underwent monotherapy with hydrochlorothiazide (HCTZ) for 4 weeks. In univariate analyses, the TT genotype was significantly associated with lower plasma renin activity ( p <0.001) and higher plasma aldosterone concentration ( p <0.019). In addition, mean declines in systolic blood pressure (SBP) and diastolic blood pressure (DBP) were 6±2 mmHg ( p <0.001) and 5±1 mmHg ( p <0.001) greater, respectively, in TT than in CC homozygotes. Responses in heterozygotes were intermediate between homozygote groups. Other univariate predictors of greater declines in both SBP and DBP were higher pretreatment BP ( p <0.001), black race ( p <0.001), female gender ( p <0.005), older age ( p ≤0.064), lower waist/hip ratio (p≤0.087), and lower plasma renin activity ( p <0.001). Even after considering effects of these other predictors, the TT genotype remained associated with greater declines in SBP ( p <0.058) and DBP ( p <0.001). Thus, the C825T polymorphism of GNB3 may help identify individuals with EHT who are more responsive to diuretic therapy and may have Na + -dependent BP elevation.

P10 Coronary artery calcification (CAC), quantified by electron beam computed tomography, provides a measure of extent of coronary atherosclerosis in asymptomatic and symptomatic individuals. It is well recognized that coronary artery disease and CAC are complex processes that are influenced by many genetic and environmental factors. Several aspects of CAC quantity make genetic studies of it particularly challenging. These include the non-normality of CAC frequency distributions and the strong age- and sex-dependence of CAC levels. To search for regions that may contain genes influencing CAC quantity, a genome-wide scan was conducted on a sample enriched for hypertension from the Epidemiology of Coronary Artery Calcification Study. The genome-wide linkage scan was based on 370 anonymous markers spanning the 22 autosomal chromosomes in 29 families. The analysis sample consisted of 20 affected sib-pairs, 8 affected sib-trios and 1 family with 6 affected siblings. Here, affection status was defined as being above the 70th sex- and age-specific percentile based on a reference population. Mode-of-inheritance-free qualitative linkage analysis was used to analyze sib-pair data. Maximum multipoint lod scores were calculated using the computer software Genehunter-2. Weighted analyses were conducted using Holmans “possible triangle” constraint. The overall maximum lod scores were observed on chromosomes 6 (lod=2.22) and 10 (lod=3.24). Additional lod scores greater than 1.0 were observed on chromosomes 2, 14, 15 and 17. These results identify regions that may harbor genes involved in high sex- and age-specific levels of CAC (supported by HL46292).