Most genetic variants associated with complex traits are hypothesized to regulate gene expression. To understand the genetics underlying gene expression variability, we characterized 14,324 RNA-sequencing samples from the Trans-Omics for Precision Medicine program and performed expression and splicing quantitative trait locus (e/sQTL) analyses in six tissues and cell types, including whole blood ( n = 6454) and lung ( n = 1291). We detected tens of thousands of secondary cis-e/sQTLs, showing that secondary cis-e/sQTL discovery remains unsaturated. We fine-mapped UK Biobank–derived genome-wide association study (GWAS) signals from 164 traits and identified e/sQTL colocalizations for 10,611 GWAS signals, including 7096 that colocalize with secondary e/sQTLs. Our results suggest that even larger e/sQTL analyses will uncover additional secondary e/sQTLs, further benefiting GWAS interpretation.
Association Between Whole Blood–Derived Mitochondrial DNA Copy Number, Low‐Density Lipoprotein Cholesterol, and Cardiovascular Disease Risk
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.
Myocardial Fibrosis and Cardiomyopathy Risk: A Genetic Link in the MESA
BACKGROUND: Common genetic variants are associated with risk for hypertrophic cardiomyopathy and dilated cardiomyopathy and with left ventricular (LV) traits. Whether these variants are associated with myocardial fibrosis, an important pathophysiological mediator of cardiomyopathy, is unknown. METHODS: Multi-Ethnic Study of Atherosclerosis participants with T1-mapping cardiac magnetic resonance imaging in-whom extracellular volume was assessed, and genotyping information was available were included (N=1255). Log extracellular volume (%) was regressed on 50 candidate single nucleotide polymorphisms (previously identified to be associated with hypertrophic cardiomyopathy, dilated cardiomyopathy, and LV traits) adjusting for age, sex, diabetes, blood pressure, and principal components of ancestry. Ancestry-specific results were pooled by fixed-effect meta-analyses. Gene knockdown experiments were performed in human cardiac fibroblasts. RESULTS: The SMARCB1 rs2186370 intronic variant (minor allele frequency: 0.18 in White and 0.50 in Black participants), previously identified as a risk variant for dilated cardiomyopathy and hypertrophic cardiomyopathy, was significantly associated with increased extracellular volume ( P= 0.0002) after adjusting for confounding clinical variables. The SMARCB1 rs2070458 locus previously associated with increased LV wall thickness and mass was similarly significantly associated with increased extracellular volume ( P= 0.0002). The direction of effect was similar in all 4 ancestry groups, but the effect was strongest in Black participants. The variants are strong expression quantitative loci in human LV tissue and associated with genotype-dependent decreased expression of SMARCB1 ( P= 7.3×10 -22 ). SMARCB1 knockdown in human cardiac fibroblasts resulted in increased TGF (transforming growth factor)-β1–mediated α-smooth muscle actin and collagen expression. CONCLUSIONS: Common genetic variation in SMARCB1 previously associated with risk for cardiomyopathies and increased LV wall thickness is associated with increased cardiac magnetic resonance imaging–based myocardial fibrosis and increased TGF-β1 mediated myocardial fibrosis in vitro. Whether these findings suggest a pathophysiologic link between myocardial fibrosis and cardiomyopathy risk remains to be proven.
Abstract 11869: Myocardial Fibrosis and Cardiomyopathy Risk: A Genetic Link in the Multi-Ethnic Study of Atherosclerosis
Introduction: Genome-wide association studies have identified common genetic variants associated with hypertrophic (HCM) and dilated cardiomyopathy (DCM) and left ventricular (LV) characteristics. However, the association of these variants with myocardial fibrosis (MF), an important pathophysiological mediator of cardiomyopathy, is unknown. Hypothesis: We used the Multi-Ethnic Study of Atherosclerosis (MESA), to evaluate the association of these loci with interstitial MF measured by cardiac MRI (CMR). Methods: Unrelated MESA participants with T1-mapping CMR and no history of MI or heart failure were included (n=1,235). We used the extracellular volume (ECV) fraction as a well-established surrogate of interstitial MF. Genotyping was performed by Affymetrix Human SNP array 6.0, and imputation was performed using the TOPMED reference panel. Log ECV(%) was regressed on 50 candidate SNPs, in each ethnicity/race group, adjusting for age, sex, and population stratification principal component. The race-specific results were pooled by fixed-effect meta-analyses using METASOFT. Results: The SMARCB1 rs2186370 intronic variant, previously identified as a risk variant for DCM and HCM, was associated with ECV (p=0.0004). The rs2186370 A allele (frequency: 0.18 in Whites, and 0.50 in African Americans) was associated with a 1.6% increase in ECV. The rs2070458 A allele (another SMARCB1 intronic variant), in complete linkage disequilibrium with rs2186370, was known previously to be associated with increased LV wall thickness, mass, and concentricity. The direction of effect was similar in all four ancestry groups, but the effect was strongest in African Americans. These variants are strong expression quantitative loci (eQTLs) in human LV tissue and are associated with genotype-dependent expression of SMARCB1 ( p =7.3х10 -22 ) and MMP11 ( p =2.5х10 -24 ). Conclusions: In conclusion, variants in two linked SMARCB1 loci previously associated with HCM, DCM and LV traits, are associated with CMR-based interstitial MF. These loci are eQTL not only for SMARCB1 but also for MMP11 an important mediator of MF . These findings underscore the important link between MF and CM, and potentially provide novel mechanistic targets in understanding the development of CM.
The genetic determinants of recurrent somatic mutations in 43,693 blood genomes
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.
Cross-Ancestry Investigation of Venous Thromboembolism Genomic Predictors
Background: Venous thromboembolism (VTE) is a life-threatening vascular event with environmental and genetic determinants. Recent VTE genome-wide association studies (GWAS) meta-analyses involved nearly 30 000 VTE cases and identified up to 40 genetic loci associated with VTE risk, including loci not previously suspected to play a role in hemostasis. The aim of our research was to expand discovery of new genetic loci associated with VTE by using cross-ancestry genomic resources. Methods: We present new cross-ancestry meta-analyzed GWAS results involving up to 81 669 VTE cases from 30 studies, with replication of novel loci in independent populations and loci characterization through in silico genomic interrogations. Results: In our genetic discovery effort that included 55 330 participants with VTE (47 822 European, 6320 African, and 1188 Hispanic ancestry), we identified 48 novel associations, of which 34 were replicated after correction for multiple testing. In our combined discovery-replication analysis (81 669 VTE participants) and ancestry-stratified meta-analyses (European, African, and Hispanic), we identified another 44 novel associations, which are new candidate VTE-associated loci requiring replication. In total, across all GWAS meta-analyses, we identified 135 independent genomic loci significantly associated with VTE risk. A genetic risk score of the significantly associated loci in Europeans identified a 6-fold increase in risk for those in the top 1% of scores compared with those with average scores. We also identified 31 novel transcript associations in transcriptome-wide association studies and 8 novel candidate genes with protein quantitative-trait locus Mendelian randomization analyses. In silico interrogations of hemostasis and hematology traits and a large phenome-wide association analysis of the 135 GWAS loci provided insights to biological pathways contributing to VTE, with some loci contributing to VTE through well-characterized coagulation pathways and others providing new data on the role of hematology traits, particularly platelet function. Many of the replicated loci are outside of known or currently hypothesized pathways to thrombosis. Conclusions: Our cross-ancestry GWAS meta-analyses identified new loci associated with VTE. These findings highlight new pathways to thrombosis and provide novel molecules that may be useful in the development of improved antithrombosis treatments.
Genomic and phenomic landscape of clonal hematopoiesis in over a million ancestrally diverse participants
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.
Abstract P256: Association of Genetic Variants on the Metabochip With Urine Albumin-Creatinine Ratio in African Americans: The Population Architecture using Genomics and Epidemiology (PAGE) Study
Several loci have been discovered and replicated for urine albumin-creatinine ratio (ACR), a measure of kidney function, in European populations; however, other ethnic populations such as African-Americans have been less studied for this trait. We examined the association of genetic variants on the Metabochip with natural log transformed ACR (ln(ACR)) in 4,629 African-Americans from the Population Architecture using Genomics and Epidemiology (PAGE) Study which includes participants from the Atherosclerosis Risk in Communities (ARIC), Cardiovascular Health Study (CHS), Coronary Artery Risk Development in Young Adults (CARDIA), and Epidemiologic Architecture for Genes Linked to the Environment (EAGLE-BioVU) studies. We then replicated our findings in 1,533 African-American participants from the Multi-Ethnic Study of Atherosclerosis (MESA). Linear regression models were used to examine associations, assuming additive effects of the SNPs. Within each study, models were adjusted for age, field center (where applicable), and principal components of ancestry, and stratified by sex. Fixed effects inverse variance weighted meta-analysis was used to combine results from contributing studies. Array-wide statistical significance was concluded for SNP associations of p-value less than 2.8 x 10-7. Mean ln(ACR) ranged from 1.4 to 2.9 mg/g, average age ranged from 44 to 77 years, and percent female participants ranged from 59-65% in the contributing studies. Four SNPs were significantly associated array-wide with ACR in the discovery sample: rs17399841 in GRHL2 (β=-0.67, MAF=0.01, p=2.7x10-9), rs153914 near SPATA9/RHOBTB3 (β=-0.59, MAF=0.01, p=1.0x10-7), rs17439845 in LEF1 (β=-0.49, MAF=0.02, p=1.3x10-7), and rs11836243 near ZNF641/ANP32D (β=0.23, MAF=0.12, p=1.8x10-7). Of these four loci, only the association of rs17439845 with ACR replicated in MESA (β=-0.36, p=0.038), with a combined discovery-replication beta of -0.46 and p value of 1.8x10-8. Combined discovery-replication for rs17399841, rs153914, and rs11836243 no longer achieved array-wide significance (p-values 9.84x10-7, 1.98x10-5, and 1.32x10-4, respectively). LEF1 is a mediator in the Wnt signaling pathway, and has been previously implicated in systemic lupus erythematosus and osteoporosis, as well as previously associated with waist circumference, echocardiographic traits, bone mineral density, and blood count traits. Continuing studies in diverse populations are warranted to identify a set of loci important in contributing to the variation in ACR.
Genetic architecture of smoking: Evaluating rare variant contribution from deep whole-genome sequencing of up to 26,000 individuals
Abstract Background Across complex traits, common variants explain only a modest amount of variance, with SNP-heritability consistently below heritability estimates from close relatives. Here, we examined the contribution of rare variant to tobacco use risk in up to 26,000 individuals of European ancestry in the Trans-Omics for Precision Medicine (TOPMed) program with whole genome sequence (WGS;~30X coverage).Method We grouped about 35 million genetic variants by their minor allele frequencies (MAF) and linkage disequilibrium (LD) and estimated SNP-heritability for age of smoking initiation (N = 14,747), cigarettes smoked per day (N = 15,425), smoking cessation (N = 17,871) and initiation (N = 26,340) using linear mixed model. Rare variant population structure is detected and adjusted for by permutation procedure. We estimated an upper bound for narrow-sense heritability for tobacco use using available pedigrees consisting of close relatives in TOPMed.Results Rare variants with MAF 0.1–0.01%, mostly from non-protein altering region, accounted for 26% of variation in age of initiation and 15% for cessation. Follow-up analysis indicated that about one-third of these rare variants contribtion is potentially confounded with rare variants structure even after adjusting for principal components. After further conservative adjustment of population structure, we estimated SNP-based heritability to be 0.21 (SE = 0.08) for age of initiation, 0.15 (0.06) for cigarettes per day, 0.21 (0.09) for cessation, and 0.24 (0.07) for initiation, 1.8–4.5 times higher than previous SNP-based estimates. Our pedigree-based upper-bound for SNP-based heritability ranged from 0.18–0.35.Conclusion The substantial contribution of rare variants for several smoking phenotypes sheds light on the missing heritability and genetic etiology of tobacco use. It also informs fine-mapping strategies since the majority of the rare variant contribution was located in non-coding regulatory regions.
Genomic and transcriptomic association studies identify 16 novel susceptibility loci for venous thromboembolism
In this work related to familial aggregation of familial venous thromboembolism, the investigators report genomic and transcriptomic association of 16 novel susceptibility loci for venous thromboembolism.