Search

The epithelial NaWe explored the association of low frequency and rare variants in the genes encoding ENaC subunits, with systolic blood pressure, diastolic blood pressure, mean arterial pressure, and pulse pressure. Using whole-genome sequencing data from 14 studies participating in the Trans-Omics in Precision Medicine Whole-Genome Sequencing Program, and sequence kernel association tests.We found that variants inOur results suggest that variants in extrarenal ENaCs, in addition to ENaCs expressed in kidneys, influence blood pressure and kidney function.

Background: The availability of whole-genome sequencing data in large studies has enabled the assessment of coding and noncoding variants across the allele frequency spectrum for their associations with blood pressure. Methods: We conducted a multiancestry whole-genome sequencing analysis of blood pressure among 51 456 Trans-Omics for Precision Medicine and Centers for Common Disease Genomics program participants (stage-1). Stage-2 analyses leveraged array data from UK Biobank (N=383 145), Million Veteran Program (N=318 891), and Reasons for Geographic and Racial Differences in Stroke (N=10 643) participants, along with whole-exome sequencing data from UK Biobank (N=199 631) participants. Results: Two blood pressure signals achieved genome-wide significance in meta-analyses of stage-1 and stage-2 single variant findings ( P -8 ). Among them, a rare intergenic variant at novel locus, LOC100506274 , was associated with lower systolic blood pressure in stage-1 (beta [SE]=−32.6 [6.0]; P =4.99×10 -8 ) but not stage-2 analysis ( P =0.11). Furthermore, a novel common variant at the known INSR locus was suggestively associated with diastolic blood pressure in stage-1 (beta [SE]=−0.36 [0.07]; P =4.18×10 -7 ) and attained genome-wide significance in stage-2 (beta [SE]=−0.29 [0.03]; P =7.28×10 -23 ). Nineteen additional signals suggestively associated with blood pressure in meta-analysis of single and aggregate rare variant findings ( P -6 and P -4 , respectively). Discussion: We report one promising but unconfirmed rare variant for blood pressure and, more importantly, contribute insights for future blood pressure sequencing studies. Our findings suggest promise of aggregate analyses to complement single variant analysis strategies and the need for larger, diverse samples, and family studies to enable robust rare variant identification.

Genome-wide association studies (GWAS) have successfully identified thousands of single nucleotide polymorphisms (SNPs) associated with complex traits; however, the identified SNPs account for a fraction of trait heritability, and identifying the functional elements through which genetic variants exert their effects remains a challenge. Recent evidence suggests that SNPs associated with complex traits are more likely to be expression quantitative trait loci (eQTL). Thus, incorporating eQTL information can potentially improve power to detect causal variants missed by traditional GWAS approaches. Using genomic, transcriptomic, and platelet phenotype data from the Genetic Study of Atherosclerosis Risk family-based study, we investigated the potential to detect novel genomic risk loci by incorporating information from eQTL in the relevant target tissues (i.e., platelets and megakaryocytes) using established statistical principles in a novel way. Permutation analyses were performed to obtain family-wise error rates for eQTL associations, substantially lowering the genome-wide significance threshold for SNP-phenotype associations. In addition to confirming the well known association between PEAR1 and platelet aggregation, our eQTL-focused approach identified a novel locus (rs1354034) and gene (ARHGEF3) not previously identified in a GWAS of platelet aggregation phenotypes. A colocalization analysis showed strong evidence for a functional role of this eQTL.

Introduction: Markers of systemic and vascular inflammation have been associated with cerebral small vessel disease (cSVD). However, their association with cSVD lesions within individual white matter regions remains unclear. We hypothesized that there are different associations of serum cytokines with the periventricular white matter hyperintensity (PVWMH) and deep WMH (DWMH). Methods: Apparently healthy participants aged ≥50 years with brain MRI in the Genetic Study of Atherosclerosis Risk were studied. Circulating inflammatory markers measured included: hs-CRP, IL-6, MCP1, ICAM, TNF-α, VCAM and IL-4. The relative association of these markers with PVWMH and DWMH within individuals were determined using multilevel linear models adjusted for age and cardiovascular risk factors. Results: A total of 435 subjects (age 58.9±6.1 years, 58.4% female, 39.5% black) with moderate WMH burden were studied (Total WMH: 1703[607-4626] mm 3 , PVWMH :871[152-2593] mm 3 , DWMH: 388[157-1023] mm 3 ). Significant negative associations of total WMH volumes and hs-CRP (β=-0.2, P=0.04) as well as IL-6 (β=-0.15, p=0.04) were observed. The associations of total WMH with MCP-1, ICAM, VCAM, IL-4 and TNF-α were not significant. In the multilevel models, designed to compare difference between DWMH and PVWMH within individuals, there were significant differences in the relative association of PVWMH vs DWMH with hs-CRP and IL-6 (p=0.001 & p=0.04, respectively). DWMH was more strongly associated with hs-CRP and IL-6 (β=0.08, β=0.14, respectively) than PVWMH, which had negative association with these markers (β=-0.35, β=-0.05, respectively). There were no differences in the association of PVWMH vs DWMH and the other markers. Conclusions: DWMH is positively associated with increasing IL-6 and CRP when compared to PVWMH which is negatively associated with them. This suggests that systemic inflammation is more preferentially associated with DWMH lesions.

Many common and rare variants associated with hematologic traits have been discovered through imputation on large-scale reference panels. However, the majority of genome-wide association studies (GWASs) have been conducted in Europeans, and determining causal variants has proved challenging. We performed a GWAS of total leukocyte, neutrophil, lymphocyte, monocyte, eosinophil, and basophil counts generated from 109,563,748 variants in the autosomes and the X chromosome in the Trans-Omics for Precision Medicine (TOPMed) program, which included data from 61,802 individuals of diverse ancestry. We discovered and replicated 7 leukocyte trait associations, including (1) the association between a chromosome X, pseudo-autosomal region (PAR), noncoding variant located between cytokine receptor genes (CSF2RA and CLRF2) and lower eosinophil count; and (2) associations between single variants found predominantly among African Americans at the S1PR3 (9q22.1) and HBB (11p15.4) loci and monocyte and lymphocyte counts, respectively. We further provide evidence indicating that the newly discovered eosinophil-lowering chromosome X PAR variant might be associated with reduced susceptibility to common allergic diseases such as atopic dermatitis and asthma. Additionally, we found a burden of very rare FLT3 (13q12.2) variants associated with monocyte counts. Together, these results emphasize the utility of whole-genome sequencing in diverse samples in identifying associations missed by European-ancestry-driven GWASs.

Platelets play a key role in thrombosis and hemostasis. Platelet count (PLT) and mean platelet volume (MPV) are highly heritable quantitative traits, with hundreds of genetic signals previously identified, mostly in European ancestry populations. We here utilize whole genome sequencing (WGS) from NHLBI’s Trans-Omics for Precision Medicine initiative (TOPMed) in a large multi-ethnic sample to further explore common and rare variation contributing to PLT (n = 61 200) and MPV (n = 23 485). We identified and replicated secondary signals at MPL (rs532784633) and PECAM1 (rs73345162), both more common in African ancestry populations. We also observed rare variation in Mendelian platelet-related disorder genes influencing variation in platelet traits in TOPMed cohorts (not enriched for blood disorders). For example, association of GP9 with lower PLT and higher MPV was partly driven by a pathogenic Bernard-Soulier syndrome variant (rs5030764, p.Asn61Ser), and the signals at TUBB1 and CD36 were partly driven by loss of function variants not annotated as pathogenic in ClinVar (rs199948010 and rs571975065). However, residual signal remained for these gene-based signals after adjusting for lead variants, suggesting that additional variants in Mendelian genes with impacts in general population cohorts remain to be identified. Gene-based signals were also identified at several genome-wide association study identified loci for genes not annotated for Mendelian platelet disorders (PTPRH, TET2, CHEK2), with somatic variation driving the result at TET2. These results highlight the value of WGS in populations of diverse genetic ancestry to identify novel regulatory and coding signals, even for well-studied traits like platelet traits.

Background There is interest in deriving megakaryocytes (MKs) from pluripotent stem cells (iPSC) for biological studies. We previously found that genomic structural integrity and genotype concordance is maintained in iPSC-derived MKs. Objective To establish a comprehensive dataset of genes and proteins expressed in iPSC-derived MKs. Methods iPSCs were reprogrammed from peripheral blood mononuclear cells (MNCs) and MKs were derived from the iPSCs in 194 healthy European American and African American subjects. mRNA was isolated and gene expression measured by RNA sequencing. Protein expression was measured in 62 of the subjects using mass spectrometry. Results and conclusions MKs expressed genes and proteins known to be important in MK and platelet function and demonstrated good agreement with previous studies in human MKs derived from CD34+ progenitor cells. The percent of cells expressing the MK markers CD41 and CD42a was consistent in biological replicates, but variable across subjects, suggesting that unidentified subject-specific factors determine differentiation of MKs from iPSCs. Gene and protein sets important in platelet function were associated with increasing expression of CD41/42a, while those related to more basic cellular functions were associated with lower CD41/42a expression. There was differential gene expression by the sex and race (but not age) of the subject. Numerous genes and proteins were highly expressed in MKs but not known to play a role in MK or platelet function; these represent excellent candidates for future study of hematopoiesis, platelet formation, and/or platelet function.

Summary Whole-genome sequencing (WGS), a powerful tool for detecting novel coding and non-coding disease-causing variants, has largely been applied to clinical diagnosis of inherited disorders. Here we leveraged WGS data in up to 62,653 ethnically diverse participants from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program and assessed statistical association of variants with seven red blood cell (RBC) quantitative traits. We discovered 14 single variant-RBC trait associations at 12 genomic loci, which have not been reported previously. Several of the RBC trait-variant associations (RPN1, ELL2, MIDN, HBB, HBA1, PIEZO1, and G6PD) were replicated in independent GWAS datasets imputed to the TOPMed reference panel. Most of these discovered variants are rare/low frequency, and several are observed disproportionately among non-European Ancestry (African, Hispanic/Latino, or East Asian) populations. We identified a 3 bp indel p.Lys2169del (g.88717175_88717177TCT[4]) (common only in the Ashkenazi Jewish population) of PIEZO1, a gene responsible for the Mendelian red cell disorder hereditary xerocytosis (MIM: 194380 ), associated with higher mean corpuscular hemoglobin concentration (MCHC). In stepwise conditional analysis and in gene-based rare variant aggregated association analysis, we identified several of the variants in HBB, HBA1, TMPRSS6, and G6PD that represent the carrier state for known coding, promoter, or splice site loss-of-function variants that cause inherited RBC disorders. Finally, we applied base and nuclease editing to demonstrate that the sentinel variant rs112097551 (nearest gene RPN1) acts through a cis-regulatory element that exerts long-range control of the gene RUVBL1 which is essential for hematopoiesis. Together, these results demonstrate the utility of WGS in ethnically diverse population-based samples and gene editing for expanding knowledge of the genetic architecture of quantitative hematologic traits and suggest a continuum between complex trait and Mendelian red cell disorders.

Genome-wide association studies have identified common variants associated with platelet-related phenotypes, but because these variants are largely intronic or intergenic, their link to platelet biology is unclear. In 290 normal subjects from the GeneSTAR Research Study (110 African Americans [AAs] and 180 European Americans [EAs]), we generated whole-genome sequence data from whole blood and RNA sequence data from extracted nonribosomal RNA from 185 induced pluripotent stem cell-derived megakaryocyte (MK) cell lines (platelet precursor cells) and 290 blood platelet samples from these subjects. Using eigenMT software to select the peak single-nucleotide polymorphism (SNP) for each expressed gene, and meta-analyzing the results of AAs and EAs, we identify (q-value < 0.05) 946 cis-expression quantitative trait loci (eQTLs) in derived MKs and 1830 cis-eQTLs in blood platelets. Among the 57 eQTLs shared between the 2 tissues, the estimated directions of effect are very consistent (98.2% concordance). A high proportion of detected cis-eQTLs (74.9% in MKs and 84.3% in platelets) are unique to MKs and platelets compared with peak-associated SNP-expressed gene pairs of 48 other tissue types that are reported in version V7 of the Genotype-Tissue Expression Project. The locations of our identified eQTLs are significantly enriched for overlap with several annotation tracks highlighting genomic regions with specific functionality in MKs, including MK-specific DNAse hotspots, H3K27-acetylation marks, H3K4-methylation marks, enhancers, and superenhancers. These results offer insights into the regulatory signature of MKs and platelets, with significant overlap in genes expressed, eQTLs detected, and enrichment within known superenhancers relevant to platelet biology.

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%., The goals, resources and design of the NHLBI Trans-Omics for Precision Medicine (TOPMed) programme are described, and analyses of rare variants detected in the first 53,831 samples provide insights into mutational processes and recent human evolutionary history.

Introduction: We have demonstrated that asymptomatic cerebral small vessel disease (cSVD) measured by white matter hyperintensity volume is associated with reduced manipulative manual dexterity on the Grooved Peg Board Test (GPBT) in middle-aged healthy individuals with a family history of early coronary artery disease. In this current study, we aim to identify the association of subcortical white matter microstructural impairment measured by diffusion tensor imaging, manual dexterity measured by GPBT and circulating serums ceramide, another marker for white matter injury. We hypothesize that lower regional fractional anisotropy (rFA) is associated with worse performance on GPBT and elevated serum ceramides in the same study population. Methods: rFA of 48 regions representing the subcortical white matters were analyzed in GeneSTAR participants in addition to serum ceramides and GPBT scores. Unadjusted univariable analyses with Bonferroni correction for multiple comparisons were completed using Spearman correlation for testing the associations between ceramides, rFA of subcortical white matter, and GPBT performance. Subsequently, sensitivity analyses were performed after excluding the participants that had any physical limitation that may influence their performance on GPBT. Finally, in the adjusted analysis using generalized estimating equation, linear regression models were performed for the areas that met significance threshold in the unadjusted analyses. Results: 112 subjects (age [49 ± 11], 51% female, 39.3% African American) were included. Adjusted analyses for the significant correlations that met the Bonferroni correction threshold in the unadjusted univariable analyses identified significant negative associations between rFA of the right fornix (RF) and log-GPBT score (β = −0.497, p = 0.037). In addition, rFA of RF negatively correlated with log serum ceramide levels (C18: β = −0.03, p = 0.003, C20: β = −0.0002, p = 0.004) and rFA of left genu of corpus callosum negatively correlated with log C18 level (β = −0.0103, p = 0.027). Conclusions: These results demonstrate that subcortical microstructural white matter disruption is associated with elevated serum ceramides and reduced manual dexterity in a population with cSVD. These findings suggest that injury to white matter tracts undermines neural networks, with functional consequences in a middle-aged population with cardiovascular risk factors.

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.

Introduction: The periventricular white matter is more sensitive to alterations in the hemodynamics of the systemic circulation than the deep white matter due to differences in its vascular supply and relationship to systemic circulation. We hypothesize that periventricular white matter hyperintensity volume (PVWMH) shows greater association than deep white matter hyperintensity (DWMH) volume with vascular property measures of arterial stiffness (AS) and cardiovascular remodeling, important indicators of the systemic circulation. Methods: Participants in the Genetic Study of Atherosclerosis Risk ≥ 50 years of age with brain MRI were selected. Vascular property measures included pulse pressure (PP), hypertension response to exercise (HRE), diastolic brachial artery (BA) diameter, diastolic common carotid artery (CCA) diameter, CCA distensibility coefficient, and left ventricular function. PVWMH and DWMH were measured using semi-automated segmentation approach. The relative associations of these parameters with PVWMH and DWMH volumes as multiple measures within the same individual further nested within families were determined using multilevel linear models. We also determined if differences between PVWMH and DWMH associations were modified by age. Results: 426 participants (age 59.0±6.1 years, 57.5% females, and 39.7% African American) were included. Within the same subject, PVWMH volume showed greater association than DWMH volume with PP (P=0.002), HRE (p=0.04), diastolic BA diameter (p=0.012), and diastolic CCA diameter (p=0.04) independent of age and cardiovascular risk factors. This difference in PVWMH versus DWMH association with vascular properties did not differ at any age threshold. Conclusions: We show for the first time that PVWMH volume has greater association than DWMH volume, independent of age, with vascular property measurements of AS and cardiovascular remodeling. Our findings suggest that changes in the systemic circulation may affect PVWMH and DWMH differently adding to a growing body of literature that PVWMH and DWMH may represent different phenotypes of cerebral small vessel disease.

Introduction: White matter hyperintensity (WMH) has been associated with manual dexterity performance in apparently healthy middle-aged population. However, the relative effect of WMH topographic distribution on cognition remains unknown. We hypothesized that periventricular WMH (PVWMH) is more strongly associated with manual dexterity and cognitive measures compared to deep WMH (DWMH) given its proximity to the white matter tracts involved in motor and executive functions. Methods: Apparently healthy middle-aged participants in the Genetic Study of Atherosclerosis Risk with brain MRI and cognitive measurements were selected. Cognitive measures included delayed word recall short (DWRS) and long (DWRL), digit symbol substitution (DSS), word fluency (WF), digit span forward test (DSFT), backward (DSBT) and total (DSTS), and grooved peg board test (GPBT) including the dominant and non-dominant hand (GPBTD, GPBTnD). PVWMH and DWMH volumes were measured using semi-automated segmentation paradigm. The relative associations of these measurements with PVWMH and DWMH volumes as multiple measurements within the same individual further nested within families were determined using multilevel linear models adjusting for intracranial volume, age, sex and cardiovascular risk factors. Results: 782 subjects were included (age 51.1±10.6, 58.4% female and 39% African American, PVWMH volume 437.4 [0-1646] mm 3 , DWMH volume 273.5 [76-665] mm 3 ). Within the same subject, PVWMH showed greater association than DWMH with worse performance on GPBT, GPBTD, GPBTnD, DWRS, and DWRL (P=0.006, 0.01, 0.01, 0.01 and 0.03, respectively) after adjusting for age and other cardiovascular risk factors. There was no difference in the association of PVWMH versus DWMH volumes with DSS, WF, DSFT, DSBT, DSTS (P=0.19, 0.52, 0.13, 0.6, and 0.6, respectively). Conclusions: PVWMH shows greater association than DWMH with cognitive measures of manipulative manual dexterity and recall suggesting that PVWMH and DWMH may have different cognitive phenotypes in middle-aged population with cerebral small vessel disease.

Large-scale gene sequencing studies for complex traits have the potential to identify causal genes with therapeutic implications. We performed gene-based association testing of blood lipid levels with rare (minor allele frequency < 1%) predicted damaging coding variation by using sequence data from >170,000 individuals from multiple ancestries: 97,493 European, 30,025 South Asian, 16,507 African, 16,440 Hispanic/Latino, 10,420 East Asian, and 1,182 Samoan. We identified 35 genes associated with circulating lipid levels; some of these genes have not been previously associated with lipid levels when using rare coding variation from population-based samples. We prioritize 32 genes in array-based genome-wide association study (GWAS) loci based on aggregations of rare coding variants; three (EVI5, SH2B3, and PLIN1) had no prior association of rare coding variants with lipid levels. Most of our associated genes showed evidence of association among multiple ancestries. Finally, we observed an enrichment of gene-based associations for low-density lipoprotein cholesterol drug target genes and for genes closest to GWAS index single-nucleotide polymorphisms (SNPs). Our results demonstrate that gene-based associations can be beneficial for drug target development and provide evidence that the gene closest to the array-based GWAS index SNP is often the functional gene for blood lipid levels. This work was supported by a grant from the Swedish Research Council (2016-06830) and grants from the National Heart, Lung, and Blood Institute (NHLBI): R01HL142711 and R01HL127564. Please refer to the supplemental information for the full acknowledgements. Sí

A diverse set of driver genes, such as regulators of DNA methylation, RNA splicing, and chromatin remodeling, have been associated with pre-malignant clonal expansion of hematopoietic stem cells (HSCs). The factors mediating expansion of these mutant clones remain largely unknown, partially due to a paucity of large cohorts with longitudinal blood sampling. To circumvent this limitation, we developed and validated a method to infer clonal expansion rate from single timepoint data called PACER (passenger-approximated clonal expansion rate). Applying PACER to 5,071 persons with clonal hematopoiesis accurately recapitulated the known fitness effects due to different driver mutations. A genome-wide association study of PACER revealed that a common inherited polymorphism in the TCL1A promoter was associated with slower clonal expansion. Those carrying two copies of this protective allele had up to 80% reduced odds of having driver mutations in TET2, ASXL1, SF3B1, SRSF2, and JAK2, but not DNMT3A. TCL1A was not expressed in normal or DNMT3A-mutated HSCs, but the introduction of mutations in TET2 or ASXL1 by CRISPR editing led to aberrant expression of TCL1A and expansion of HSCs in vitro. These effects were abrogated in HSCs from donors carrying the protective TCL1A allele. Our results indicate that the fitness advantage of multiple common driver genes in clonal hematopoiesis is mediated through TCL1A activation. PACER is an approach that can be widely applied to uncover genetic and environmental determinants of pre-malignant clonal expansion in blood and other tissues.

Genome-wide association studies (GWAS) have successfully identified thousands of single nucleotide polymorphisms (SNPs) associated with complex traits; however, the identified SNPs account for a fraction of trait heritability, and identifying the functional elements through which genetic variants exert their effects remains a challenge. Recent evidence suggests that SNPs associated with complex traits are more likely to be expression quantitative trait loci (eQTL). Thus, incorporating eQTL information can potentially improve power to detect causal variants missed by traditional GWAS approaches. Using genomic, transcriptomic, and platelet phenotype data from the Genetic Study of Atherosclerosis Risk family-based study, we investigated the potential to detect novel genomic risk loci by incorporating information from eQTL in the relevant target tissues (i.e. platelets and megakaryocytes). Permutation analyses were performed to obtain family-wise error rates for eQTL associations, substantially lowering the genome-wide significance threshold for SNP-phenotype associations. In addition to confirming the well known association between PEAR1 and platelet aggregation, our eQTL focused approach identified a novel locus (rs1354034) and gene (ARHGEF3) not previously identified in a GWAS of platelet aggregation phenotypes. A colocalization analysis showed strong evidence for a functional role of this eQTL.

Background and Purpose: We aim to determine, in healthy high-risk adults, the association between subclinical coronary artery disease and white matter hyperintensity (WMH) volume and location, independent of atherosclerotic risk factors. Methods: Seven hundred eighty-two asymptomatic first-degree relatives of index cases with early-onset coronary artery disease ( Results: Mean age was 51 years ±10, with 58% women and 39% African American people. Participants with any coronary plaque had 52% larger WMH volumes than those without plaque (95% CI, 0.24–0.59). Per 1% greater coronary plaque volume, total WMH volumes were 0.07% larger (95% CI, 0.04–0.10). Every 1% higher total coronary plaque volume was associated with 5.03% larger deep WMH volume (95% CI, 4.67–5.38), 5.10% PVWMH larger volume (95% CI, 4.72–5.48), and 2.74% larger cuff volume (95% CI, 2.38–3.09) with differences in this association when comparing deep WMH to PVWMH ( P interaction, 0.001) or cuff ( P interaction, Conclusions: In healthy, high-risk individuals, the presence and volume of coronary artery plaque are associated with larger WMH volumes, appearing the strongest for PVWMH. These findings in high-risk families suggest a disease relationship in 2 different vascular beds, beyond traditional risk factors, possibly due to genetic predisposition.

BackgroundUse of targeted exome-arrays with common, rare variants and functionally enriched variation has led to discovery of new genes contributing to population variation in risk factors. Plasminogen activator-inhibitor 1 (PAI-1), tissue plasminogen activator (tPA), and the plasma product D-dimer are important components of the fibrinolytic system. There have been few large-scale genome-wide or exome-wide studies of PAI-1, tPA, and D-dimer.ObjectivesWe sought to discover new genetic loci contributing to variation in these traits using an exome-array approach.MethodsCohort-level analyses and fixed effects meta-analyses of PAI-1 (n=15603), tPA (n=6876,) and D-dimer (n=19306) from 12 cohorts of European ancestry with diverse study design were conducted, including single-variant analyses and gene-based burden testing.ResultsFive variants located in NME7, FGL1, and the fibrinogen locus, all associated with D-dimer levels, achieved genome-wide significance (P&nbsp

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 Trans-Omics for Precision Medicine (TOPMed) program (N=63,302) and UK Biobank (N=48,658). Bidirectional Mendelian randomization studies were consistent with LTL lengthening increasing propensity to develop CHIP, but CHIP then in turn hastening LTL shortening. We also demonstrated evidence of modest mediation between CHIP and CAD by LTL. Our data promote an understanding of potential causal relationships across CHIP and LTL toward prevention of CAD.

Large-scale gene sequencing studies for complex traits have the potential to identify causal genes with therapeutic implications. We performed gene-based association testing of blood lipid levels with rare (minor allele frequency170,000 individuals from multiple ancestries: 97,493 European, 30,025 South Asian, 16,507 African, 16,440 Hispanic/Latino, 10,420 East Asian, and 1,182 Samoan. We identified 35 genes associated with circulating lipid levels. Ten of these: ALB, SRSF2, JAK2, CREB3L3, TMEM136, VARS, NR1H3, PLA2G12A, PPARG and STAB1 have not been implicated for lipid levels using rare coding variation in population-based samples. We prioritize 32 genes identified in array-based genome-wide association study (GWAS) loci based on gene-based associations, of which three: EVI5, SH2B3, and PLIN1, had no prior evidence of rare coding variant associations. Most of the associated genes showed evidence of association in multiple ancestries. Also, we observed an enrichment of gene-based associations for low-density lipoprotein cholesterol drug target genes, and for genes closest to GWAS index single nucleotide polymorphisms (SNP). Our results demonstrate that gene-based associations can be beneficial for drug target development and provide evidence that the gene closest to the array-based GWAS index SNP is often the functional gene for blood lipid levels.

Pharmaceutical drugs targeting dyslipidemia and cardiovascular disease (CVD) may increase the risk of fatty liver disease and other metabolic disorders. To identify potential novel CVD drug targets without these adverse effects, we perform genome-wide analyses of participants in the HUNT Study in Norway (n = 69,479) to search for protein-altering variants with beneficial impact on quantitative blood traits related to cardiovascular disease, but without detrimental impact on liver function. We identify 76 (11 previously unreported) presumed causal protein-altering variants associated with one or more CVD- or liver-related blood traits. Nine of the variants are predicted to result in loss-of-function of the protein. This includes ZNF529:p.K405X, which is associated with decreased low-density-lipoprotein (LDL) cholesterol (P = 1.3 × 10−8) without being associated with liver enzymes or non-fasting blood glucose. Silencing of ZNF529 in human hepatoma cells results in upregulation of LDL receptor and increased LDL uptake in the cells. This suggests that inhibition of ZNF529 or its gene product should be prioritized as a novel candidate drug target for treating dyslipidemia and associated CVD., Drugs targeting cardiovascular disease (CVD) can have negative consequences for liver function. Here, the authors combine genome wide analyses on 69,479 individuals to identify loss-of-function variants with beneficial effects on CVD-related traits without negative impacts on liver function.

Whole genome sequencing (WGS), a powerful tool for detecting novel coding and non-coding disease-causing variants, has largely been applied to clinical diagnosis of inherited disorders. Here we leveraged WGS data in up to 62,653 ethnically diverse participants from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program and assessed statistical association of variants with seven red blood cell (RBC) quantitative traits. We discovered 14 single variant-RBC trait associations at 12 genomic loci. Several of the RBC trait-variant associations (RPN1, ELL2, MIDN, HBB, HBA1, PIEZO1, G6PD) were replicated in independent GWAS datasets imputed to the TOPMed reference panel. Most of these newly discovered variants are rare/low frequency, and several are observed disproportionately among non-European Ancestry (African, Hispanic/Latino, or East Asian) populations. We identified a 3bp indel p.Lys2169del (common only in the Ashkenazi Jewish population) of PIEZO1, a gene responsible for the Mendelian red cell disorder hereditary xerocytosis [OMIM 194380], associated with higher MCHC. In stepwise conditional analysis and in gene-based rare variant aggregated association analysis, we identified several of the variants in HBB, HBA1, TMPRSS6, and G6PD that represent the carrier state for known coding, promoter, or splice site loss-of-function variants that cause inherited RBC disorders. Finally, we applied base and nuclease editing to demonstrate that the sentinel variant rs112097551 (nearest gene RPN1) acts through a cis-regulatory element that exerts long-range control of the gene RUVBL1 which is essential for hematopoiesis. Together, these results demonstrate the utility of WGS in ethnically-diverse population-based samples and gene editing for expanding knowledge of the genetic architecture of quantitative hematologic traits and suggest a continuum between complex trait and Mendelian red cell disorders.

Autosomal genetic analyses of blood lipids have yielded key insights for coronary heart disease (CHD). However, X chromosome genetic variation is understudied for blood lipids in large sample sizes. We now analyze genetic and blood lipid data in a high-coverage whole X chromosome sequencing study of 65,322 multi-ancestry participants and perform replication among 456,893 European participants. Common alleles on chromosome Xq23 are strongly associated with reduced total cholesterol, LDL cholesterol, and triglycerides (min P = 8.5 × 10−72), with similar effects for males and females. Chromosome Xq23 lipid-lowering alleles are associated with reduced odds for CHD among 42,545 cases and 591,247 controls (P = 1.7 × 10−4), and reduced odds for diabetes mellitus type 2 among 54,095 cases and 573,885 controls (P = 1.4 × 10−5). Although we observe an association with increased BMI, waist-to-hip ratio adjusted for BMI is reduced, bioimpedance analyses indicate increased gluteofemoral fat, and abdominal MRI analyses indicate reduced visceral adiposity. Co-localization analyses strongly correlate increased CHRDL1 gene expression, particularly in adipose tissue, with reduced concentrations of blood lipids., The influence of X chromosome genetic variation on blood lipids and coronary heart disease (CHD) is not well understood. Here, the authors analyse X chromosome sequencing data across 65,322 multi-ancestry individuals, identifying associations of the Xq23 locus with lipid changes and reduced risk of CHD and diabetes mellitus.

Introduction: Induced pluripotent stem cells (iPSC) have a significant capacity for self-renewal and can differentiate into any somatic cell type, such as megakaryocytes (MK) and vascular endothelial cells (EC). We have previously shown that iPSC-derived MKs (IPSC-MK) upregulate the gene expression levels of irreversible glycolytic enzymes committing glucose to glycolysis in comparison to their parental iPSC lines (Vaidya et al., 2016). This pattern differed from past comparisons between somatic source cells and iPSCs (Chung et al., 2007; Varum et al., 2011; Zhang et al., 2011). To further corroborate our data, we extended our analysis to iPSC-derived ECs (iPSC-EC). Methods: We studied expression of glycolytic enzyme genes in 30 human iPSC lines and paired iPSC-EC and iPSC-MK lines using RNA-seq. iPSC-EC to iPSC expression ratios for 11 glycolytic enzymes were estimated using a multilevel mixed-effect model regression, and a significance threshold of p < 0.05/11 = 0.0045 was used. A bivariate statistical analysis was used to assess the correlation between iPSC-EC and iPSC-MK data. Results: Of the 11 glycolytic enzymes studied, 4 were highly expressed in iPSC-ECs (Log2(iPSC-EC/iPSC) ≥ 1) while 1 was highly expressed in iPSCs (Log2(iPSC/iPSC-EC) ≥ 1). Transcripts for functionally irreversible enzymes including ADPGK, PFKP, and PKLR displayed higher expression levels in iPSC-ECs compared to iPSCs. Similarly, 6 glycolytic genes demonstrated higher expression in iPSC-MKs, while 2 were more highly expressed in iPSCs. The functionally irreversible enzymes HK1, ADPGK, and PFK exhibited higher expression levels in iPSC-MKs. The ratios of gene expression (derived-lines:iPSCs) for all 37 genes in the iPSC-MK results correlated with those of the iPSC-EC results: Pearson’s r = 0.52, p = 0.001. Conclusion: Both iPSC-ECs and iPSC-MKs upregulate the expression of glycolytic enzymes compared to their parental iPSCs. In light of data showing that differentiated cells are more metabolically active and have more mitochondria than iPSCs (Chandra et al., 2020), these results suggest that differentiation of iPSCs is accompanied by an increase in both glycolytic and TCA cycle activities, with glycolytic generation of pyruvate used for mitochondrial production of ATP.

Human induced pluripotent stem cell (hiPSC) lines have previously been generated through the NHLBI sponsored NextGen program at nine individual study sites. Here, we examined the structural integrity of 506 hiPSC lines as determined by copy number variations (CNVs). We observed that 149 hiPSC lines acquired 258 CNVs relative to donor DNA. We identified six recurrent regions of CNVs on chromosomes 1, 2, 3, 16 and 20 that overlapped with cancer associated genes. Furthermore, the genes mapping to regions of acquired CNVs show an enrichment in cancer related biological processes (IL6 production) and signaling cascades (JNK cascade & NFκB cascade). The genomic region of instability on chr20 (chr20q11.2) includes transcriptomic signatures for cancer associated genes such as ID1, BCL2L1, TPX2, PDRG1 and HCK. Of these HCK shows statistically significant differential expression between carrier and non-carrier hiPSC lines. Overall, while a low level of genomic instability was observed in the NextGen generated hiPSC lines, the observation of structural instability in regions with known cancer associated genes substantiates the importance of systematic evaluation of genetic variations in hiPSCs before using them as disease/research models.

Background and Purpose: Periventricular white matter hyperintensities (WMH; PVWMH) and deep WMH (DWMH) are regional classifications of WMH and reflect proposed differences in cause. In the first study, to date, we undertook genome-wide association analyses of DWMH and PVWMH to show that these phenotypes have different genetic underpinnings. Methods: Participants were aged 45 years and older, free of stroke and dementia. We conducted genome-wide association analyses of PVWMH and DWMH in 26,654 participants from CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology), ENIGMA (Enhancing Neuro-Imaging Genetics Through Meta-Analysis), and the UKB (UK Biobank). Regional correlations were investigated using the genome-wide association analyses -pairwise method. Cross-trait genetic correlations between PVWMH, DWMH, stroke, and dementia were estimated using LDSC. Results: In the discovery and replication analysis, for PVWMH only, we found associations on chromosomes 2 ( NBEAL ), 10q23.1 ( TSPAN14/FAM231A ), and 10q24.33 ( SH3PXD2A). In the much larger combined meta-analysis of all cohorts, we identified ten significant regions for PVWMH: chromosomes 2 (3 regions), 6, 7, 10 (2 regions), 13, 16, and 17q23.1. New loci of interest include 7q36.1 ( NOS3 ) and 16q24.2. In both the discovery/replication and combined analysis, we found genome-wide significant associations for the 17q25.1 locus for both DWMH and PVWMH. Using gene-based association analysis, 19 genes across all regions were identified for PVWMH only, including the new genes: CALCRL (2q32.1), KLHL24 (3q27.1), VCAN (5q27.1), and POLR2F (22q13.1). Thirteen genes in the 17q25.1 locus were significant for both phenotypes. More extensive genetic correlations were observed for PVWMH with small vessel ischemic stroke. There were no associations with dementia for either phenotype. Conclusions: Our study confirms these phenotypes have distinct and also shared genetic architectures. Genetic analyses indicated PVWMH was more associated with ischemic stroke whilst DWMH loci were implicated in vascular, astrocyte, and neuronal function. Our study confirms these phenotypes are distinct neuroimaging classifications and identifies new candidate genes associated with PVWMH only.

Background: Whether obesity in the absence of cardiometabolic abnormalities increases the risk of subclinical atherosclerosis remains debatable. Methods: We compared coronary artery calcium (CAC) across obesity phenotypes in 505 Whites and 306 African Americans (AAs) from the Genetic Study of Aspirin Responsiveness (GeneSTAR); and 875 South Asians (SAs) from MASALA (Mediators of Atherosclerosis in South Asians Living in America). Using BMI and the International Harmonized Definition of metabolic syndrome, participants were categorized as Metabolically Healthy Non-Obese (MHNO), Metabolically Healthy Obese (MHO), Metabolically Unhealthy Non-Obese (MUNO), and Metabolically Unhealthy Obese (MUO). Results: The prevalence of any detectable CAC was similar in MHNO versus MHO among all ethnic groups (Whites 35 versus 30%, AAs 30 versus 23%, SAs 29 versus 26%; respectively, all p>0.05). In contrast, the prevalence of CAC was the highest in the MUNO and MUO groups in all ethnic groups. After adjusting for potential confounders; the odds of CAC in MUNO and MUO compared to MHNO were significantly higher in Whites and SAs; but not in AAs. Conclusion: Metabolically unhealthy Whites, AAs, and SAs had higher CAC scores regardless of obesity status suggesting that metabolic health is more closely associated with atherosclerosis than obesity alone especially in Whites and SAs. Disclosure M. Al-Sofiani: Speaker’s Bureau; Self; Novo Nordisk Inc. D. Vaidya: None. A.M. Kanaya: None. L. Yanek: None. M. Mongraw-Chaffin: None. N.R. Kandula: None. L. Becker: None. B.G. Kral: None. D. Becker: None. N.M. Maruthur: Other Relationship; Self; Johns Hopkins HealthCare Solutions. C. Lee: Employee; Self; Eli Lilly and Company. Funding National Institutes of Health (R01HL093009, K24HL112827); National Center for Research Resources (UL1RR024131, M01-RR000052); National Heart, Lung, and Blood Institute (U01HL72518)

Introduction: We previously reported impaired manipulative manual dexterity measured by grooved peg-board test (GPBT) in healthypeople with increased total brain ischemic white matter hyperintensities. We hypothesized that preclinical white matter tract damage (WMTD), measured as lower regional fractional anisotropy (rFA), is associated with functional decrements measured with GPBT, and with a biomarker of white matter disease, serum ceramides. Methods: Participants were recruited from the GeneSTAR study in families with a history of coronary artery disease. Diffusion tensor imaging was acquired on a Phillips 3T MRI and rFA was calculated using the 181 Eve segmented atlas. 48 susceptible regions for WMTD were analyzed.Liquid chromatography was used to measure ceramides (C). Spearman’s rank correlation (SRC) was used to test statistical associations and Bonferroni method was used to correct for multiple comparisons. Multivariate analysis (GEE) adjusting for age, sex, race and education was performed for associations that were significant in the univariate analyses. Result: Subjects (n=112) were 49±11 years; 49% male and 39% African American. rFA was correlated with decrements on GPBT in the left (L) anterior corona radiata, L genu of the corpus callosum (LGCC), L external capsule, and L & right fornix (RF) (P=0.0009, 0.0008, 0.0002, 0.0003 and 8.2E-07, respectively). Additionally, rFA was negatively correlated with C18:0 in the RF and LGCC (P=1.27E-05 and 2.58E-05, respectively) and C20:0 in RF (P=6.09E-06). C18:0 was correlated with decrements on GPBT (P=0.0025). GEE showed associations between decremental GPBT and rFA in the RF (B = -0.52, 95% CI -1 - 0.03, p= 0.037), between rFA in the RF and Ceramide (C18:0, B= -0.03, 95% CI -0.046--0.0073, p=0.007, C20:0, B=-0.0002, 95%CI -0.0003-0, P= 0.018), as well as rFA in LGCC and Ceramide (C18:0, B=-0.0097, 95% CI -0.019--0.006, P=0.037). GBPT and ceramides were not associated when adjusted for rFA (B = 0.0002 95% CI -0.052 - 0.052 P= 0.99). Conclusion: These analyses support our hypothesis that preclinical WMTD represented by reduced rFA, may undermine neural networks affecting manipulative manual dexterity, with concomitant elevation in serum ceramides in healthy high risk people.

Background Genome-wide association studies have identified multiple genomic loci associated with coronary artery disease, but most are common variants in non-coding regions that provide limited information on causal genes and etiology of the disease. To overcome the limited scope that common variants provide, we focused our investigation on low-frequency and rare sequence variations primarily residing in coding regions of the genome. Methods and results Using samples of individuals of European ancestry from ten cohorts within the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, both cross-sectional and prospective analyses were conducted to examine associations between genetic variants and myocardial infarction (MI), coronary heart disease (CHD), and all-cause mortality following these events. For prevalent events, a total of 27,349 participants of European ancestry, including 1831 prevalent MI cases and 2518 prevalent CHD cases were used. For incident cases, a total of 55,736 participants of European ancestry were included (3,031 incident MI cases and 5,425 incident CHD cases). There were 1,860 all-cause deaths among the 3,751 MI and CHD cases from six cohorts that contributed to the analysis of all-cause mortality. Single variant and gene-based analyses were performed separately in each cohort and then meta-analyzed for each outcome. A low-frequency intronic variant (rs988583) in PLCL1 was significantly associated with prevalent MI (OR = 1.80, 95% confidence interval: 1.43, 2.27; P = 7.12 × 10−7). We conducted gene-based burden tests for genes with a cumulative minor allele count (cMAC) ≥ 5 and variants with minor allele frequency (MAF) < 5%. TMPRSS5 and LDLRAD1 were significantly associated with prevalent MI and CHD, respectively, and RC3H2 and ANGPTL4 were significantly associated with incident MI and CHD, respectively. No loci were significantly associated with all-cause mortality following a MI or CHD event. Conclusion This study identified one known locus (ANGPTL4) and four new loci (PLCL1, RC3H2, TMPRSS5, and LDLRAD1) associated with cardiovascular disease risk that warrant further investigation.

BACKGROUND AND PURPOSE: Traditional cardiovascular risk factors have been associated with white matter disease. Because hypertension results in vascular stiffness and impaired cerebral perfusion, we hypothesized that it would be the most relevant risk factor for microstructural white matter disruption in apparently healthy middle-aged individuals with a family history of early-onset coronary artery disease. MATERIALS AND METHODS: This was a cross-sectional analysis of participants in the Genetic Study of Atherosclerosis Risk with DTI. Regional fractional anisotropy of 181 segmented brain regions was measured using Eve WM Atlas. Risk factors were examined using univariate analysis for 48 regions representing deep WM structures. Minimal multivariable linear regression models adjusting for age, sex, and race and maximal linear regression models adjusting for cardiovascular risk factors were performed for regions meeting the Bonferroni threshold in the initial analysis. RESULTS: Included were 116 subjects (mean age, 49 ± 11 years; 57% men) with a moderate load of cardiovascular risk factors. Subjects with hypertension had significantly lower regional fractional anisotropy in the right cingulum and left stria terminalis in the minimal and maximal regression models. Additionally, there was lower regional fractional anisotropy in the left fornix in the maximal model and right sagittal stratum in the minimal model. Systolic blood pressure values were significantly associated with regional fractional anisotropy in the left superior longitudinal fasciculus in the maximal model. There were no significant differences among regional fractional anisotropy values for other cardiovascular risk factors. CONCLUSIONS: In middle-aged apparently healthy individuals with susceptibility to vascular disease, among all known cardiovascular risk factors, hypertension was associated with microstructural WM disruption.

Lipoprotein(a), Lp(a), is a modified low-density lipoprotein particle that contains apolipoprotein(a), encoded by LPA, and is a highly heritable, causal risk factor for cardiovascular diseases that varies in concentrations across ancestries. Here, we use deep-coverage whole genome sequencing in 8392 individuals of European and African ancestry to discover and interpret both single-nucleotide variants and copy number (CN) variation associated with Lp(a). We observe that genetic determinants between Europeans and Africans have several unique determinants. The common variant rs12740374 associated with Lp(a) cholesterol is an eQTL for SORT1 and independent of LDL cholesterol. Observed associations of aggregates of rare non-coding variants are largely explained by LPA structural variation, namely the LPA kringle IV 2 (KIV2)-CN. Finally, we find that LPA risk genotypes confer greater relative risk for incident atherosclerotic cardiovascular diseases compared to directly measured Lp(a), and are significantly associated with measures of subclinical atherosclerosis in African Americans., Circulating lipoprotein(a) is an important risk factor for cardiovascular disease and shows variability between different ethnic groups. Here, Zekavat et al. perform whole-genome sequencing in individuals of European and African ancestries and find ancestry-specific genetic determinants for lipoprotein(a) levels.

Coronary artery disease (CAD) remains a major cause of mortality and morbidity worldwide. The aggregation of activated platelets on a ruptured atherosclerotic plaque is a critical step in most acute cardiovascular events like myocardial infarction. Platelet aggregation both at baseline and after aspirin is highly heritable. Genome-wide association studies (GWAS) have identified a common variant within the first intron of the platelet endothelial aggregation receptor1 (PEAR1), to be robustly associated with platelet aggregation. In this study, we used targeted deep sequencing to fine-map the prior GWAS peak and identify additional rare variants of PEAR1 that account for missing heritability in platelet aggregation within the GeneSTAR families. In this study, 1709 subjects (1043 European Americans, EA and 666 African Americans, AA) from families in the GeneSTAR study were included. In vitro platelet aggregation in response to collagen, ADP and epinephrine was measured at baseline and 14 days after aspirin therapy (81 mg/day). Targeted deep sequencing of PEAR1 in addition to 2kb of upstream and downstream of the gene was performed. Under an additive genetic model, the association of single variants of PEAR1 with platelet aggregation phenotypes were examined. Additionally, we examined the association between the burden of PEAR1 rare non-synonymous variants and platelet aggregation phenotypes. Of 532 variants identified through sequencing, the intron 1 variant, rs12041331, was significantly associated with all platelet aggregation phenotypes at baseline and after platelet inhibition with aspirin therapy. rs12566888, which is in linkage disequilibrium with rs12041331, was associated with platelet aggregation phenotypes but to a lesser extent. In the EA families, the burden of PEAR1 missense variants was associated with platelet aggregation after aspirin therapy when the platelets were stimulated with epinephrine (p = 0.0009) and collagen (p = 0.03). In AAs, the burden of PEAR1 missense variants was associated, to a lesser degree, with platelet aggregation in response to epinephrine (p = 0.02) and ADP (p = 0.04). Our study confirmed that the GWAS-identified variant, rs12041331, is the strongest variant associated with platelet aggregation both at baseline and after aspirin therapy in our GeneSTAR families in both races. We identified additional association of rare missense variants in PEAR1 with platelet aggregation following aspirin therapy. However, we observed a racial difference in the contribution of these rare variants to the platelet aggregation, most likely due to higher residual missing heritability of platelet aggregation after accounting for rs12041331 in the EAs compared to AAs.

SummaryInhibition of platelet reactivity is a common therapeutic strategy in secondary prevention of cardiovascular disease. Genetic and environmental factors influence inter-individual variation in platelet reactivity. Identifying genes that contribute to platelet reactivity can reveal new biological mechanisms and possible therapeutic targets. Here, we examined rare coding variation to identify genes associated with platelet reactivity in a population-based cohort. To do so, we performed whole exome sequencing in the Framingham Heart Study and conducted single variant and gene-based association tests against platelet reactivity to collagen, adenosine diphosphate (ADP), and epinephrine agonists in up to 1,211 individuals. Single variant tests revealed no significant associations (pSupplementary Material to this article is available online at www.thrombosis-online.com.

Age is the dominant risk factor for most chronic human diseases; yet the mechanisms by which aging confers this risk are largely unknown.1 Recently, the age-related acquisition of somatic mutations in regenerating hematopoietic stem cell populations was associated with both hematologic cancer incidence2–4 and coronary heart disease prevalence.5 Somatic mutations with leukemogenic potential may confer selective cellular advantages leading to clonal expansion, a phenomenon termed ‘Clonal Hematopoiesis of Indeterminate Potential’ (CHIP).6 Simultaneous germline and somatic whole genome sequence analysis now provides the opportunity to identify root causes of CHIP. Here, we analyze high-coverage whole genome sequences from 97,691 participants of diverse ancestries in the NHLBI TOPMed program and identify 4,229 individuals with CHIP. We identify associations with blood cell, lipid, and inflammatory traits specific to different CHIP genes. Association of a genome-wide set of germline genetic variants identified three genetic loci associated with CHIP status, including one locus at TET2 that was African ancestry specific. In silico-informed in vitro evaluation of the TET2 germline locus identified a causal variant that disrupts a TET2 distal enhancer. Aggregates of rare germline loss-of-function variants in CHEK2, a DNA damage repair gene, predisposed to CHIP acquisition. Overall, we observe that germline genetic variation altering hematopoietic stem cell function and the fidelity of DNA-damage repair increase the likelihood of somatic mutations leading to CHIP.

© 2020, The Author(s), under exclusive licence to Springer Nature Limited. Age is the dominant risk factor for most chronic human diseases, but the mechanisms through which ageing confers this risk are largely unknown1. The age-related acquisition of somatic mutations that lead to clonal expansion in regenerating haematopoietic stem cell populations has recently been associated with both haematological cancer2–4 and coronary heart disease5—this phenomenon is termed clonal haematopoiesis of indeterminate potential (CHIP)6. Simultaneous analyses of germline and somatic whole-genome sequences provide the opportunity to identify root causes of CHIP. Here we analyse high-coverage whole-genome sequences from 97,691 participants of diverse ancestries in the National Heart, Lung, and Blood Institute Trans-omics for Precision Medicine (TOPMed) programme, and identify 4,229 individuals with CHIP. We identify associations with blood cell, lipid and inflammatory traits that are specific to different CHIP driver genes. Association of a genome-wide set of germline genetic variants enabled the identification of three genetic loci associated with CHIP status, including one locus at TET2 that was specific to individuals of African ancestry. In silico-informed in vitro evaluation of the TET2 germline locus enabled the identification of a causal variant that disrupts a TET2 distal enhancer, resulting in increased self-renewal of haematopoietic stem cells. Overall, we observe that germline genetic variation shapes haematopoietic stem cell function, leading to CHIP through mechanisms that are specific to clonal haematopoiesis as well as shared mechanisms that lead to somatic mutations across tissues.

Most genome-wide association and fine-mapping studies to date have been conducted in individuals of European descent, and genetic studies of populations of Hispanic/Latino and African ancestry are limited. In addition, these populations have more complex linkage disequilibrium structure. In order to better define the genetic architecture of these understudied populations, we leveraged >100,000 phased sequences available from deep-coverage whole genome sequencing through the multi-ethnic NHLBI Trans-Omics for Precision Medicine (TOPMed) program to impute genotypes into admixed African and Hispanic/Latino samples with genome-wide genotyping array data. We demonstrated that using TOPMed sequencing data as the imputation reference panel improves genotype imputation quality in these populations, which subsequently enhanced gene-mapping power for complex traits. For rare variants with minor allele frequency (MAF) < 0.5%, we observed a 2.3- to 6.1-fold increase in the number of well-imputed variants, with 11–34% improvement in average imputation quality, compared to the state-of-the-art 1000 Genomes Project Phase 3 and Haplotype Reference Consortium reference panels. Impressively, even for extremely rare variants with minor allele count 86%. Subsequent association analyses of TOPMed reference panel-imputed genotype data with hematological traits (hemoglobin (HGB), hematocrit (HCT), and white blood cell count (WBC)) in ~21,600 African-ancestry and ~21,700 Hispanic/Latino individuals identified associations with two rare variants in the HBB gene (rs33930165 with higher WBC [p = 8.8x10-15] in African populations, rs11549407 with lower HGB [p = 1.5x10-12] and HCT [p = 8.8x10-10] in Hispanics/Latinos). By comparison, neither variant would have been genome-wide significant if either 1000 Genomes Project Phase 3 or Haplotype Reference Consortium reference panels had been used for imputation. Our findings highlight the utility of the TOPMed imputation reference panel for identification of novel rare variant associations not previously detected in similarly sized genome-wide studies of under-represented African and Hispanic/Latino populations., Author summary Admixed African and Hispanic/Latino populations remain understudied in genetic studies of complex diseases. These populations have more complex linkage disequilibrium (LD) structure that can impair mapping of variants. Genotype imputation represents an approach to improve genome coverage, especially for rare or ancestry-specific variation; however, these understudied populations also have smaller relevant imputation reference panels. In this study, we leveraged >100,000 phased sequences generated from the multi-ethnic NHLBI TOPMed project for imputation in ~21,600 individuals of African ancestry (AAs) and ~21,700 Hispanics/Latinos. We demonstrated substantially higher imputation quality for low frequency and rare variants in comparison to the 1000 Genomes Project and Haplotype Reference Consortium reference panels. Analysis of quantitative hematological traits led to the discovery of associations with two rare variants in the HBB gene; one of these variants was replicated in an independent sample, and the other is known to cause anemia in the homozygous state. By comparison, the same HBB variants would not have been genome-wide significant using current reference panels due to lower imputation quality. Our findings demonstrate the power of TOPMed whole genome sequencing data for imputation and subsequent association analysis in admixed African and Hispanic/Latino populations.