Vasan, Ramachandran S

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.

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.

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: The study of hypertrophic cardiomyopathy (HCM) can yield insight into the mechanisms underlying the complex trait of cardiac hypertrophy. To date, most genetic variants associated with HCM have been found in sarcomeric genes. Here, we describe a novel HCM-associated variant in the noncanonical Wnt signaling interactor WTIP (Wilms tumor interacting protein) and provide evidence of a role for WTIP in complex disease. Methods: In a family affected by HCM, we used exome sequencing and identity-by-descent analysis to identify a novel variant in WTIP (p.Y233F). We knocked down WTIP in isolated neonatal rat ventricular myocytes with lentivirally delivered short hairpin ribonucleic acids and in Danio rerio via morpholino injection. We performed weighted gene coexpression network analysis for WTIP in human cardiac tissue, as well as association analysis for WTIP variation and left ventricular hypertrophy. Finally, we generated induced pluripotent stem cell–derived cardiomyocytes from patient tissue, characterized size and calcium cycling, and determined the effect of verapamil treatment on calcium dynamics. Results: WTIP knockdown caused hypertrophy in neonatal rat ventricular myocytes and increased cardiac hypertrophy, peak calcium, and resting calcium in D rerio . Network analysis of human cardiac tissue indicated WTIP as a central coordinator of prohypertrophic networks, while common variation at the WTIP locus was associated with human left ventricular hypertrophy. Patient-derived WTIP p.Y233F-induced pluripotent stem cell–derived cardiomyocytes recapitulated cellular hypertrophy and increased resting calcium, which was ameliorated by verapamil. Conclusions: We demonstrate that a novel genetic variant found in a family with HCM disrupts binding to a known Wnt signaling protein, misregulating cardiomyocyte calcium dynamics. Further, in orthogonal model systems, we show that expression of the gene WTIP is important in complex cardiac hypertrophy phenotypes. These findings, derived from the observation of a rare Mendelian disease variant, uncover a novel disease mechanism with implications across diverse forms of cardiac hypertrophy.

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.

Abstract Aims Mitral valve prolapse (MVP) is a common valvular heart disease with a prevalence of &gt;2% in the general adult population. Despite this high incidence, there is a limited understanding of the molecular mechanism of this disease, and no medical therapy is available for this disease. We aimed to elucidate the genetic basis of MVP in order to better understand this complex disorder. Methods and results We performed a meta-analysis of six genome-wide association studies that included 4884 cases and 434 649 controls. We identified 14 loci associated with MVP in our primary analysis and 2 additional loci associated with a subset of the samples that additionally underwent mitral valve surgery. Integration of epigenetic, transcriptional, and proteomic data identified candidate MVP genes including LMCD1, SPTBN1, LTBP2, TGFB2, NMB, and ALPK3. We created a polygenic risk score (PRS) for MVP and showed an improved MVP risk prediction beyond age, sex, and clinical risk factors. Conclusion We identified 14 genetic loci that are associated with MVP. Multiple analyses identified candidate genes including two transforming growth factor-β signalling molecules and spectrin β. We present the first PRS for MVP that could eventually aid risk stratification of patients for MVP screening in a clinical setting. These findings advance our understanding of this common valvular heart disease and may reveal novel therapeutic targets for intervention.

Genome-wide association studies have identified >50 common variants associated with kidney function, but these variants do not fully explain the variation in eGFR. We performed a two-stage meta-analysis of associations between genotypes from the Illumina exome array and eGFR on the basis of serum creatinine (eGFRcrea) among participants of European ancestry from the CKDGen Consortium ( n Stage1 : 111,666; n Stage2 : 48,343). In single-variant analyses, we identified single nucleotide polymorphisms at seven new loci associated with eGFRcrea ( PPM1J , EDEM3, ACP1, SPEG, EYA4, CYP1A1 , and ATXN2L ; P Stage1 <3.7×10 −7 ), of which most were common and annotated as nonsynonymous variants. Gene-based analysis identified associations of functional rare variants in three genes with eGFRcrea, including a novel association with the SOS Ras/Rho guanine nucleotide exchange factor 2 gene, SOS2 ( P =5.4×10 −8 by sequence kernel association test). Experimental follow-up in zebrafish embryos revealed changes in glomerular gene expression and renal tubule morphology in the embryonic kidney of acp1- and sos2 -knockdowns. These developmental abnormalities associated with altered blood clearance rate and heightened prevalence of edema. This study expands the number of loci associated with kidney function and identifies novel genes with potential roles in kidney formation.

Background— Data regarding the familial aggregation of left ventricular (LV) geometry and its relations to parental heart failure (HF) are limited. Methods and Results— We evaluated concordance of LV geometry within 1093 nuclear families in 5758 participants of the original (parents) (n=2351) and offspring (n=3407) cohorts of the Framingham Heart Study undergoing routine echocardiography in mid- to late adulthood. LV geometry was categorized based on cohort- and sex-specific 80th percentile cutoffs of LV mass and relative wall thickness (RWT) into normal (both <80th percentile), concentric remodeling (LV mass <80th percentile; RWT >80th percentile), concentric hypertrophy (both >80th percentile), and eccentric hypertrophy (LV mass >80th percentile; RWT <80th percentile). Within nuclear families, LV geometry was concordant among related pairs (parent-child, sibling-sibling) ( P =0.0015) but not among unrelated spousal pairs ( P =0.60), a finding that remained unchanged after adjusting for clinical covariates known to influence LV remodeling (age, systolic blood pressure, body mass index), excluding individuals with prevalent HF and myocardial infarction, and varying the thresholds for defining LV geometry. The prevalence of abnormal LV geometry was higher in family members of affected individuals, with recurrence risks of 1.4 for concentric remodeling (95% CI, 1.2 to 1.7) and eccentric hypertrophy (95% CI, 1.1 to 1.8) and 3.9 (95% CI, 3.2 to 4.6) for concentric hypertrophy. In a subset of 1497 offspring, we observed an association between parental HF (n=458) and eccentric hypertrophy in offspring ( P <0.0001). Conclusions— Our investigation of a 2-generational community-based sample demonstrates familial aggregation of LV geometry, with the greatest recurrence risk for concentric LV geometry, and establishes an association between eccentric LV geometry and parental HF.