Search

Blood lipid traits are treatable and heritable risk factors for heart disease, a leading cause of mortality worldwide. Although genome-wide association studies (GWASs) have discovered hundreds of variants associated with lipids in humans, most of the causal mechanisms of lipids remain unknown. To better understand the biological processes underlying lipid metabolism, we investigated the associations of plasma protein levels with total cholesterol (TC), triglycerides (TG), high-density lipoprotein (HDL) cholesterol, and low-density lipoprotein (LDL) cholesterol in blood. We trained protein prediction models based on samples in the Multi-Ethnic Study of Atherosclerosis (MESA) and applied them to conduct proteome-wide association studies (PWASs) for lipids using the Global Lipids Genetics Consortium (GLGC) data. Of the 749 proteins tested, 42 were significantly associated with at least one lipid trait. Furthermore, we performed transcriptome-wide association studies (TWASs) for lipids using 9,714 gene expression prediction models trained on samples from peripheral blood mononuclear cells (PBMCs) in MESA and 49 tissues in the Genotype-Tissue Expression (GTEx) project. We found that although PWASs and TWASs can show different directions of associations in an individual gene, 40 out of 49 tissues showed a positive correlation between PWAS and TWAS signed p values across all the genes, which suggests high-level consistency between proteome-lipid associations and transcriptome-lipid associations.

BackgroundInflammation is a feature of coronary heart disease (CHD), but the role of proinflammatory microbial infection in CHD remains understudied.Methods and resultsCHD was defined in the MESA (Multi-Ethnic Study of Atherosclerosis) as myocardial infarction (251 participants), resuscitated arrest (2 participants), and CHD death (80 participants). We analyzed sequencing reads from 4421 MESA participants in the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine program using the PathSeq workflow of the Genome Analysis Tool Kit and a 65-gigabase microbial reference. Paired reads aligning to 840 microbes were detected in >1% of participants. The association of the presence of microbe reads with incident CHD (follow-up, ~18 years) was examined. First, important variables were ascertained using a single regularized Cox proportional hazard model, examining change of risk as a function of presence of microbe with age, sex, education level, Life's Simple 7, and inflammation. For variables of importance, the hazard ratio (HR) was estimated in separate (unregularized) Cox proportional hazard models including the same covariates (significance threshold Bonferroni corrected P

Clonal hematopoiesis of indeterminate potential (CHIP), whereby somatic mutations in hematopoietic stem cells confer a selective advantage and drive clonal expansion, not only correlates with age but also confers increased risk of morbidity and mortality. Here, we leverage genetically predicted traits to identify factors that determine CHIP clonal expansion rate. We used the passenger-approximated clonal expansion rate method to quantify the clonal expansion rate for 4,370 individuals in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) cohort and calculated polygenic risk scores for DNA methylation aging, inflammation-related measures and circulating protein levels. Clonal expansion rate was significantly associated with both genetically predicted and measured epigenetic clocks. No associations were identified with inflammation-related lab values or diseases and CHIP expansion rate overall. A proteome-wide search identified predicted circulating levels of myeloid zinc finger 1 and anti-Müllerian hormone as associated with an increased CHIP clonal expansion rate and tissue inhibitor of metalloproteinase 1 and glycine N-methyltransferase as associated with decreased CHIP clonal expansion rate. Together, our findings identify epigenetic and proteomic patterns associated with the rate of hematopoietic clonal expansion.

Abstract Aims Intercellular adhesion molecule‐1 (ICAM‐1) facilitates inflammation via leucocyte recruitment and has been implicated in heart failure with preserved ejection fraction (HFpEF). Approximately 35% of African American individuals carry a copy of an ICAM1 missense variant (rs5491; p.K56M), which is associated with an increased risk of HFpEF. The pathways by which rs5491 increases HFpEF risk are not well defined. We evaluated the circulating immune cell profile of rs5491. Methods Among African American individuals in the Multi‐Ethnic Study of Atherosclerosis, we evaluated the associations of rs5491 with 29 circulating peripheral blood mononuclear cell subsets. The top immune cells were then related to echocardiographic measures of structure and function. Results Among 502 individuals with immune cell profiling (mean age 63 years, 51% female), 191 individuals (38%) had at least one copy of rs5491. Each additional rs5491 allele was significantly associated with higher proportions of Tc17 CD8+ cytotoxic T cells (β = 1.34, SE = 0.45, P = 9.5 × 10−5) and Tc2 CD8+ cytotoxic T cells (β = 1.19, SE = 0.44, P = 0.00012). There were no other associations noted between rs5491 and the remaining immune cells. A higher proportion of Tc17 cells was significantly associated with a higher left ventricular ejection fraction, E/e′ average and right ventricular systolic pressure (RVSP), while a higher proportion of Tc2 cells was significantly associated with a higher RVSP. Conclusions The ICAM1 p.K56M variant (rs5491) carries a distinct and inflammatory T‐cell subset profile. These cytotoxic T cells are in turn associated with alterations in cardiac function and adverse haemodynamics later in life, thus providing insight into pathways by which rs5491 may increase the risk of HFpEF.

Background: Black adults have higher incidence of all-cause mortality and worse cardiovascular disease (CVD) outcomes when compared to other U.S. populations. The Duffy chemokine receptor is not expressed on erythrocytes in a large majority of Black adults, but the clinical implications of this are unclear. Methods: Here, we investigated the relationship of Duffy receptor status, high-sensitivity C-reactive protein (hs-CRP), and mortality and incident CVD events (coronary heart disease, stroke, and heart failure) in self-identified Black members of three contemporary, longitudinal cohort studies (the Women's Health Initiative, Jackson Heart Study, and Multi-Ethnic Study of Atherosclerosis). Data on 14,358 Black participants (9023 Duffy-null and 5335 Duffy-receptor-positive, as defined using single-nucleotide polymorphism (SNP) rs2814778) were included in this analysis. Results: Duffy null was strongly associated with higher hs-CRP (meta-analysis p = 2.62 × 10-9), but the association was largely attenuated, though still marginally significant (p = 0.005), after conditioning on known CRP locus alleles in linkage disequilibrium with the Duffy gene. In our discovery cohorts, Duffy-null status appeared to be associated with a higher risk of all-cause mortality and incident stroke, though these associations were attenuated and non-significant following adjustment for traditional risk factors including hs-CRP. Moreover, the association of Duffy-null status with mortality could not be replicated in an independent sample of Black adults from the UK Biobank. Conclusions: These findings suggest that the higher levels of hs-CRP found in Duffy-null individuals may be in part independent of CRP alleles known to influence circulating levels of hs-CRP. During the follow-up of this community-based sample of Black participants, Duffy-null status was not associated with mortality or incident CVD events independently of traditional risk factors including hs-CRP.

Chronic kidney disease is a leading cause of death and disability globally and impacts individuals of African ancestry (AFR) or with ancestry in the Americas (AMS) who are under-represented in genome-wide association studies (GWASs) of kidney function. To address this bias, we conducted a large meta-analysis of GWASs of estimated glomerular filtration rate (eGFR) in 145,732 AFR and AMS individuals. We identified 41 loci at genome-wide significance (p

ObjectiveClonal hematopoiesis of indeterminate potential (CHIP) is an aging-related accumulation of somatic mutations in hematopoietic stem cells, leading to clonal expansion. CHIP presence has been implicated in atherosclerotic coronary heart disease (CHD) and all-cause mortality, but its association with incident type 2 diabetes (T2D) is unknown. We hypothesized that CHIP is associated with elevated risk of T2D.Research design and methodsCHIP was derived from whole-genome sequencing of blood DNA in the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine (TOPMed) prospective cohorts. We performed analysis for 17,637 participants from six cohorts, without prior T2D, cardiovascular disease, or cancer. We evaluated baseline CHIP versus no CHIP prevalence with incident T2D, including associations with DNMT3A, TET2, ASXL1, JAK2, and TP53 variants. We estimated multivariable-adjusted hazard ratios (HRs) and 95% CIs with adjustment for age, sex, BMI, smoking, alcohol, education, self-reported race/ethnicity, and combined cohorts' estimates via fixed-effects meta-analysis.ResultsMean (SD) age was 63.4 (11.5) years, 76% were female, and CHIP prevalence was 6.0% (n = 1,055) at baseline. T2D was diagnosed in n = 2,467 over mean follow-up of 9.8 years. Participants with CHIP had 23% (CI 1.04, 1.45) higher risk of T2D than those with no CHIP. Specifically, higher risk was for TET2 (HR 1.48; CI 1.05, 2.08) and ASXL1 (HR 1.76; CI 1.03, 2.99) mutations; DNMT3A was nonsignificant (HR 1.15; CI 0.93, 1.43). Statistical power was limited for JAK2 and TP53 analyses.ConclusionsCHIP was associated with higher incidence of T2D. CHIP mutations located on genes implicated in CHD and mortality were also related to T2D, suggesting shared aging-related pathology.

BACKGROUND. Most GWAS of plasma proteomics have focused on White individuals of European ancestry, limiting biological insight from other ancestry-enriched protein quantitative loci (pQTL). METHODS. We conducted a discovery GWAS of approximately 3,000 plasma proteins measured by the antibody-based Olink platform in 1,054 Black adults from the Jackson Heart Study (JHS) and validated our findings in the Multi-Ethnic Study of Atherosclerosis (MESA). The genetic architecture of identified pQTLs was further explored through fine mapping and admixture association analysis. Finally, using our pQTL findings, we performed a phenome-wide association study (PheWAS) across 2 large multiethnic electronic health record (EHR) systems in All of Us and BioMe. RESULTS. We identified 1,002 pQTLs for 925 protein assays. Fine mapping and admixture analyses suggested allelic heterogeneity of the plasma proteome across diverse populations. We identified associations for variants enriched in African ancestry, many in diseases that lack precise biomarkers, including c/s-pQTLs for cathepsin L (CTSL) and Siglec-9, which were linked with sarcoidosis and non-Hodgkin's lymphoma, respectively. We found concordant associations across clinical diagnoses and laboratory measurements, elucidating disease pathways, including a cis-pQTL associated with circulating CD58, WBC count, and multiple sclerosis. CONCLUSIONS. Our findings emphasize the value of leveraging diverse populations to enhance biological insights from proteomics GWAS, and we have made this resource readily available as an interactive web portal. FUNDING. NIH K08 HL161445-01A1; 5T32HL160522-03; HHSN268201600034I; HL133870., Introduction Large-scale GWAS have uncovered thousands of loci implicated in human disease. However, genotype is often far removed from phenotype, limiting insight into the pathological processes that result in disease. [...]

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

BackgroundClonal hematopoiesis of indeterminate potential (CHIP) is an age-related somatic mutation associated with incident hematologic cancer. Environmental stressors which, like air pollution, generate oxidative stress at the cellular level, may induce somatic mutations and some mutations may provide a selection advantage for persistence and expansion of specific clones.MethodsWe used data from the Multi-Ethnic Study of Atherosclerosis (MESA) N = 4,379 and the Women's Health Initiative (WHI) N = 7,701 to estimate cross-sectional associations between annual average air pollution concentrations at participant address the year before blood draw using validated spatiotemporal models. We used covariate-adjusted logistic regression to estimate risk of CHIP per interquartile range increases in particulate matter (PM2.5; 4 μg/m3) and nitrogen dioxide (NO2; 10 ppb) as ORs (95% confidence intervals).ResultsPrevalence of CHIP at blood draw (variant allele fraction > 2%) was 4.4% and 8.7% in MESA and WHI, respectively. The most common CHIP driver mutation was in DNMT3A. Neither pollutant was associated with CHIP: ORMESA PM2.5 = 1.00 (0.68-1.45), ORMESA NO2 = 1.05 (0.69-1.61), ORWHI PM2.5 = 0.97 (0.86-1.09), ORWHI NO2 = 0.98 (0.88-1.10); or with DNMT3A-driven CHIP.ConclusionsWe did not find evidence that air pollution contributes to CHIP prevalence in two large observational cohorts.ImpactThis is the first study to estimate associations between air pollution and CHIP.

Nonalcoholic fatty liver disease (NAFLD) is common and partially heritable and has no effective treatments. We carried out a genome-wide association study (GWAS) meta-analysis of imaging (n = 66,814) and diagnostic code (3,584 cases versus 621,081 controls) measured NAFLD across diverse ancestries. We identified NAFLD-associated variants at torsin family 1 member B (TOR1B), fat mass and obesity associated (FTO), cordon-bleu WH2 repeat protein like 1 (COBLL1)/growth factor receptor-bound protein 14 (GRB14), insulin receptor (INSR), sterol regulatory element-binding transcription factor 1 (SREBF1) and patatin-like phospholipase domain-containing protein 2 (PNPLA2), as well as validated NAFLD-associated variants at patatin-like phospholipase domain-containing protein 3 (PNPLA3), transmembrane 6 superfamily 2 (TM6SF2), apolipoprotein E (APOE), glucokinase regulator (GCKR), tribbles homolog 1 (TRIB1), glycerol-3-phosphate acyltransferase (GPAM), mitochondrial amidoxime-reducing component 1 (MARC1), microsomal triglyceride transfer protein large subunit (MTTP), alcohol dehydrogenase 1B (ADH1B), transmembrane channel like 4 (TMC4)/membrane-bound O-acyltransferase domain containing 7 (MBOAT7) and receptor-type tyrosine-protein phosphatase δ (PTPRD). Implicated genes highlight mitochondrial, cholesterol and de novo lipogenesis as causally contributing to NAFLD predisposition. Phenome-wide association study (PheWAS) analyses suggest at least seven subtypes of NAFLD. Individuals in the top 10% and 1% of genetic risk have a 2.5-fold to 6-fold increased risk of NAFLD, cirrhosis and hepatocellular carcinoma. These genetic variants identify subtypes of NAFLD, improve estimates of disease risk and can guide the development of targeted therapeutics.

Transcriptome prediction models built with data from European-descent individuals are less accurate when applied to different populations because of differences in linkage disequilibrium patterns and allele frequencies. We hypothesized that methods that leverage shared regulatory effects across different conditions, in this case, across different populations, may improve cross-population transcriptome prediction. To test this hypothesis, we made transcriptome prediction models for use in transcriptome-wide association studies (TWASs) using different methods (elastic net, joint-tissue imputation [JTI], matrix expression quantitative trait loci [Matrix eQTL], multivariate adaptive shrinkage in R [MASHR], and transcriptome-integrated genetic association resource [TIGAR]) and tested their out-of-sample transcriptome prediction accuracy in population-matched and cross-population scenarios. Additionally, to evaluate model applicability in TWASs, we integrated publicly available multiethnic genome-wide association study (GWAS) summary statistics from the Population Architecture using Genomics and Epidemiology (PAGE) study and Pan-ancestry genetic analysis of the UK Biobank (PanUKBB) with our developed transcriptome prediction models. In regard to transcriptome prediction accuracy, MASHR models performed better or the same as other methods in both population-matched and cross-population transcriptome predictions. Furthermore, in multiethnic TWASs, MASHR models yielded more discoveries that replicate in both PAGE and PanUKBB across all methods analyzed, including loci previously mapped in GWASs and loci previously not found in GWASs. Overall, our study demonstrates the importance of using methods that benefit from different populations' effect size estimates in order to improve TWASs for multiethnic or underrepresented populations.

Coronary artery calcification (CAC), a measure of subclinical atherosclerosis, predicts future symptomatic coronary artery disease (CAD). Identifying genetic risk factors for CAC may point to new therapeutic avenues for prevention. Currently, there are only four known risk loci for CAC identified from genome-wide association studies (GWAS) in the general population. Here we conducted the largest multi-ancestry GWAS meta-analysis of CAC to date, which comprised 26,909 individuals of European ancestry and 8,867 individuals of African ancestry. We identified 11 independent risk loci, of which eight were new for CAC and five had not been reported for CAD. These new CAC loci are related to bone mineralization, phosphate catabolism and hormone metabolic pathways. Several new loci harbor candidate causal genes supported by multiple lines of functional evidence and are regulators of smooth muscle cell-mediated calcification ex vivo and in vitro. Together, these findings help refine the genetic architecture of CAC and extend our understanding of the biological and potential druggable pathways underlying CAC.

Each human genome has tens of thousands of rare genetic variants; however, identifying impactful rare variants remains a major challenge. We demonstrate how use of personal multi-omics can enable identification of impactful rare variants by using the Multi-Ethnic Study of Atherosclerosis, which included several hundred individuals, with whole-genome sequencing, transcriptomes, methylomes, and proteomes collected across two time points, 10 years apart. We evaluated each multi-omics phenotype's ability to separately and jointly inform functional rare variation. By combining expression and protein data, we observed rare stop variants 62 times and rare frameshift variants 216 times as frequently as controls, compared to 13-27 times as frequently for expression or protein effects alone. We extended a Bayesian hierarchical model, "Watershed," to prioritize specific rare variants underlying multi-omics signals across the regulatory cascade. With this approach, we identified rare variants that exhibited large effect sizes on multiple complex traits including height, schizophrenia, and Alzheimer's disease.

BackgroundCommon 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.MethodsMulti-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.ResultsThe 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.ConclusionsCommon 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.

Although many novel gene-metabolite and gene-protein associations have been identified using high-throughput biochemical profiling, systematic studies that leverage human genetics to illuminate causal relationships between circulating proteins and metabolites are lacking. Here, we performed protein-metabolite association studies in 3,626 plasma samples from three human cohorts. We detected 171,800 significant protein-metabolite pairwise correlations between 1,265 proteins and 365 metabolites, including established relationships in metabolic and signaling pathways such as the protein thyroxine-binding globulin and the metabolite thyroxine, as well as thousands of new findings. In Mendelian randomization (MR) analyses, we identified putative causal protein-to-metabolite associations. We experimentally validated top MR associations in proof-of-concept plasma metabolomics studies in three murine knockout strains of key protein regulators. These analyses identified previously unrecognized associations between bioactive proteins and metabolites in human plasma. We provide publicly available data to be leveraged for studies in human metabolism and disease.

Multi-omics datasets are becoming more common, necessitating better integration methods to realize their revolutionary potential. Here, we introduce multi-set correlation and factor analysis (MCFA), an unsupervised integration method tailored to the unique challenges of high-dimensional genomics data that enables fast inference of shared and private factors. We used MCFA to integrate methylation markers, protein expression, RNA expression, and metabolite levels in 614 diverse samples from the Trans-Omics for Precision Medicine/Multi-Ethnic Study of Atherosclerosis multi-omics pilot. Samples cluster strongly by ancestry in the shared space, even in the absence of genetic information, while private spaces frequently capture dataset-specific technical variation. Finally, we integrated genetic data by conducting a genome-wide association study (GWAS) of our inferred factors, observing that several factors are enriched for GWAS hits and trans-expression quantitative trait loci. Two of these factors appear to be related to metabolic disease. Our study provides a foundation and framework for further integrative analysis of ever larger multi-modal genomic datasets.

Summary: Blood lipid traits are treatable and heritable risk factors for heart disease, a leading cause of mortality worldwide. Although genome-wide association studies (GWASs) have discovered hundreds of variants associated with lipids in humans, most of the causal mechanisms of lipids remain unknown. To better understand the biological processes underlying lipid metabolism, we investigated the associations of plasma protein levels with total cholesterol (TC), triglycerides (TG), high-density lipoprotein (HDL) cholesterol, and low-density lipoprotein (LDL) cholesterol in blood. We trained protein prediction models based on samples in the Multi-Ethnic Study of Atherosclerosis (MESA) and applied them to conduct proteome-wide association studies (PWASs) for lipids using the Global Lipids Genetics Consortium (GLGC) data. Of the 749 proteins tested, 42 were significantly associated with at least one lipid trait. Furthermore, we performed transcriptome-wide association studies (TWASs) for lipids using 9,714 gene expression prediction models trained on samples from peripheral blood mononuclear cells (PBMCs) in MESA and 49 tissues in the Genotype-Tissue Expression (GTEx) project. We found that although PWASs and TWASs can show different directions of associations in an individual gene, 40 out of 49 tissues showed a positive correlation between PWAS and TWAS signed p values across all the genes, which suggests high-level consistency between proteome-lipid associations and transcriptome-lipid associations.

Distinct tissue-specific mechanisms mediate insulin action in fasting and postprandial states. Previous genetic studies have largely focused on insulin resistance in the fasting state, where hepatic insulin action dominates. Here we studied genetic variants influencing insulin levels measured 2 h after a glucose challenge in >55,000 participants from three ancestry groups. We identified ten new loci (P