Large-scale genome-wide association study of coronary artery disease in genetically diverse populations

Leveraging Random Effects in Cistrome-Wide Association Studies for Decoding the Genetic Determinants of Prostate Cancer

Cistrome-wide association studies (CWAS) are pivotal for identifying genetic determinants of diseases by correlating genetically regulated cistrome states with phenotypes. Traditional CWAS typically develops a model based on cistrome and genotype data to associate predicted cistrome states with phenotypes. The random effect cistrome-wide association study (RECWAS), reevaluates the necessity of cistrome state prediction in CWAS. RECWAS utilizes either a linear model or marginal effect for initial feature selection, followed by kernel-based feature aggregation for association testing is introduced. Through simulations and analysis of prostate cancer data, a thorough evaluation of CWAS and RECWAS is conducted. The results suggest that RECWAS offers improved power compared to traditional CWAS, identifying additional genomic regions associated with prostate cancer. CWAS identified 102 significant regions, while RECWAS found 50 additional significant regions compared to CWAS, many of which are validated. Validation encompassed a range of biological evidence, including risk signals from the GWAS catalog, susceptibility genes from the DisGeNET database, and enhancer-domain scores. RECWAS consistently demonstrated improved performance over traditional CWAS in identifying genomic regions associated with prostate cancer. These findings demonstrate the benefits of incorporating kernel methods into CWAS and provide new insights for genetic discovery in complex diseases.

Genome-wide discovery for biomarkers using quantile regression at biobank scale

Genome-wide association studies (GWAS) for biomarkers important for clinical phenotypes can lead to clinically relevant discoveries. Conventional GWAS for quantitative traits are based on simplified regression models modeling the conditional mean of a phenotype as a linear function of genotype. We draw attention here to an alternative, lesser known approach, namely quantile regression that naturally extends linear regression to the analysis of the entire conditional distribution of a phenotype of interest. Quantile regression can be applied efficiently at biobank scale, while having some unique advantages such as (1) identifying variants with heterogeneous effects across quantiles of the phenotype distribution; (2) accommodating a wide range of phenotype distributions including non-normal distributions, with invariance of results to trait transformations; and (3) providing more detailed information about genotype-phenotype associations even for those associations identified by conventional GWAS. We show in simulations that quantile regression is powerful across both homogeneous and various heterogeneous models. Applications to 39 quantitative traits in the UK Biobank demonstrate that quantile regression can be a helpful complement to linear regression in GWAS and can identify variants with larger effects on high-risk subgroups of individuals but with lower or no contribution overall.

Polygenic risk score informed clinical model for improving abdominal aortic aneurysm screening

INTRODUCTION

Abdominal aortic aneurysm (AAA) is a complex disease with environmental and genetic risk factors. Polygenic risk scores (PRS) based on disease-specific risk-associated single nucleotide variants (SNVs) have demonstrated effectiveness in stratifying individual-level disease risk for cardiovascular diseases. This prospective cohort study assessed associations of PRS of AAA (PRSAAA) with risk of incident AAA, analyzed the effectiveness of a combined clinical-genetic risk model, and explored the clinical utility of the model in identifying high-risk individuals for AAA screening.

METHODS

PRSAAA was calculated using 911,440 SNVs and PRSCAD was calculated using 2,324,683 SNVs derived from mixed ancestry genome wide association studies. The UK Biobank was used as the study cohort. All individuals with complete genetic data available and no diagnosis of AAA at time of recruitment were included in the analysis and followed prospectively to assess for incident AAA. A PRS informed clinical model, Prob-AAA, was developed using clinically significant variables and PRSAAA.

RESULTS

481,105 individuals were included in the analysis with 2,668 incident AAA cases. Incident AAA increased from 0.30% to 0.93% between the lowest and highest decile of PRSAAA; similarly, severe AAA, requiring surgery and/or presenting with rupture, increased from 23% to 39% of incident AAA cases across deciles. PRSAAA was a predictor of incident AAA diagnosis (HR 2.06 [1.70 - 2.48]) independent of other clinical risk factors including male sex, older age, and smoking history. Prob-AAA was an independent predictor of incident AAA (HR 1.92 [1.69 - 2.20]), and identified 9.6% of cases of incident AAA compared to only 4.2% by PRSAAA. Current screening guidelines captured 5.7% of the overall cohort, with an incident AAA rate of approximately 3.2%. Amongst males not included by current guidelines, Prob-AAA identified an additional cohort, approximately 2% of the overall cohort, with a similar rate of incident AAA.

CONCLUSIONS

Prob-AAA, a PRS informed clinical model for AAA, improved upon the predictive power of current, clinical risk factor informed, screening guidelines for AAA.

Application of polygenic scores to a deeply phenotyped sample enriched for substance use disorders reveals extensive pleiotropy with psychiatric and somatic traits

Co-occurring psychiatric, medical, and substance use disorders (SUDs) are common, but the complex pathways leading to such comorbidities are poorly understood. A greater understanding of genetic influences on this phenomenon could inform precision medicine efforts. We used the Yale-Penn dataset, a cross-sectional sample enriched for individuals with SUDs, to examine pleiotropic effects of genetic liability for psychiatric and somatic traits. Participants completed an in-depth interview that provides information on demographics, environment, medical illnesses, and psychiatric and SUDs. Polygenic scores (PGS) for psychiatric disorders and somatic traits were calculated in European-ancestry (EUR; n = 5691) participants and, when discovery datasets were available, for African-ancestry (AFR; n = 4918) participants. Phenome-wide association studies (PheWAS) were then conducted. In AFR participants, the only PGS with significant associations was bipolar disorder (BD), all of which were with substance use phenotypes. In EUR participants, PGS for major depressive disorder (MDD), generalized anxiety disorder (GAD), post-traumatic stress disorder (PTSD), schizophrenia (SCZ), body mass index (BMI), coronary artery disease (CAD), and type 2 diabetes (T2D) all showed significant associations, the majority of which were with phenotypes in the substance use categories. For instance, PGSMDD was associated with over 200 phenotypes, 15 of which were depression-related (e.g., depression criterion count), 55 of which were other psychiatric phenotypes, and 126 of which were substance use phenotypes; and PGSBMI was associated with 138 phenotypes, 105 of which were substance related. Genetic liability for psychiatric and somatic traits is associated with numerous phenotypes across multiple categories, indicative of the broad genetic liability of these traits.

Effect of disclosing a polygenic risk score for coronary heart disease on adverse cardiovascular events: 10-year follow-up of the MI-GENES randomized clinical trial

Background

The MI-GENES clinical trial (NCT01936675), in which participants at intermediate risk of coronary heart disease (CHD) were randomized to receive a Framingham risk score (FRSg, n=103), or an integrated risk score (IRSg, n=104) that additionally included a polygenic risk score (PRS), demonstrated that after 6 months, participants randomized to IRSg had higher statin initiation and lower low-density lipoprotein cholesterol (LDL-C).

Objectives

In a post hoc 10-year follow-up analysis of the MI-GENES trial, we investigated whether disclosure of a PRS for CHD was associated with a reduction in adverse cardiovascular events.

Methods

Participants were followed from randomization beginning in October 2013 until September 2023 to ascertain adverse cardiovascular events, testing for CHD, and changes in risk factors, by blinded review of electronic health records. The primary outcome was the time from randomization to the occurrence of the first major adverse cardiovascular event (MACE), defined as cardiovascular death, non-fatal myocardial infarction, coronary revascularization, and non-fatal stroke. Statistical analyses were conducted using Cox proportional hazards regression and linear mixed-effects models.

Results

We followed all 203 participants who completed the MI-GENES trial, 100 in FRSg and 103 in IRSg (mean age at the end of follow-up: 68.2±5.2, 48% male). During a median follow-up of 9.5 years, 9 MACEs occurred in FRSg and 2 in IRSg (hazard ratio (HR), 0.20; 95% confidence interval (CI), 0.04 to 0.94; P=0.042). In FRSg, 47 (47%) underwent at least one test for CHD, compared to 30 (29%) in IRSg (HR, 0.51; 95% CI, 0.32 to 0.81; P=0.004). IRSg participants had a longer duration of statin therapy during the first four years post-randomization and a greater reduction in LDL-C for up to 3 years post-randomization. No significant differences between the two groups were observed for hemoglobin A1C, systolic and diastolic blood pressures, weight, and smoking cessation rate during follow-up.

Conclusions

The disclosure of an IRS that included a PRS to individuals at intermediate risk for CHD was associated with a lower incidence of MACE after a decade of follow-up, likely due to a higher rate of initiation and longer duration of statin therapy, leading to lower LDL-C levels.

Identification and correction for collider bias in a genome-wide association study of diabetes-related heart failure

Type 2 diabetes (T2D) is a major risk factor for heart failure (HF) and has elevated incidence among individuals with HF. Since genetics and HF can independently influence T2D, collider bias may occur when T2D (i.e., collider) is controlled for by design or analysis. Thus, we conducted a genome-wide association study (GWAS) of diabetes-related HF with correction for collider bias. We first performed a GWAS of HF to identify genetic instrumental variables (GIVs) for HF and to enable bidirectional Mendelian randomization (MR) analysis between T2D and HF. We identified 61 genomic loci, significantly associated with all-cause HF in 114,275 individuals with HF and over 1.5 million controls of European ancestry. Using a two-sample bidirectional MR approach with 59 and 82 GIVs for HF and T2D, respectively, we estimated that T2D increased HF risk (odds ratio [OR] 1.07, 95% confidence interval [CI] 1.04-1.10), while HF also increased T2D risk (OR 1.60, 95% CI 1.36-1.88). Then we performed a GWAS of diabetes-related HF corrected for collider bias due to the study design of index cases. After removing the spurious association of TCF7L2 locus due to collider bias, we identified two genome-wide significant loci close to PITX2 (chromosome 4) and CDKN2B-AS1 (chromosome 9) associated with diabetes-related HF in the Million Veteran Program and replicated the associations in the UK Biobank. Our MR findings provide strong evidence that HF increases T2D risk. As a result, collider bias leads to spurious genetic associations of diabetes-related HF, which can be effectively corrected to identify true positive loci.

Smooth muscle expression of RNA editing enzyme ADAR1 controls vascular integrity and progression of atherosclerosis

Mapping the genomic architecture of complex disease has been predicated on the understanding that genetic variants influence disease risk through modifying gene expression. However, recent discoveries have revealed that a significant burden of disease heritability in common autoinflammatory disorders and coronary artery disease is mediated through genetic variation modifying post-transcriptional modification of RNA through adenosine-to-inosine (A-to-I) RNA editing. This common RNA modification is catalyzed by ADAR enzymes, where ADAR1 edits specific immunogenic double stranded RNA (dsRNA) to prevent activation of the double strand RNA (dsRNA) sensor MDA5 ( IFIH1 ) and stimulation of an interferon stimulated gene (ISG) response. Multiple lines of human genetic data indicate impaired RNA editing and increased dsRNA sensing to be an important mechanism of coronary artery disease (CAD) risk. Here, we provide a crucial link between observations in human genetics and mechanistic cell biology leading to progression of CAD. Through analysis of human atherosclerotic plaque, we implicate the vascular smooth muscle cell (SMC) to have a unique requirement for RNA editing, and that ISG induction occurs in SMC phenotypic modulation, implicating MDA5 activation. Through culture of human coronary artery SMCs, generation of a conditional SMC specific Adar1 deletion mouse model on a pro-atherosclerosis background, and with incorporation of single cell RNA sequencing cellular profiling, we further show that Adar1 controls SMC phenotypic state, is required to maintain vascular integrity, and controls progression of atherosclerosis and vascular calcification. Through this work, we describe a fundamental mechanism of CAD, where cell type and context specific RNA editing and sensing of dsRNA mediates disease progression, bridging our understanding of human genetics and disease causality.

Distinct genetic liability profiles define clinically relevant patient strata across common diseases

Stratified medicine holds great promise to tailor treatment to the needs of individual patients. While genetics holds great potential to aid patient stratification, it remains a major challenge to operationalize complex genetic risk factor profiles to deconstruct clinical heterogeneity. Contemporary approaches to this problem rely on polygenic risk scores (PRS), which provide only limited clinical utility and lack a clear biological foundation. To overcome these limitations, we develop the CASTom-iGEx approach to stratify individuals based on the aggregated impact of their genetic risk factor profiles on tissue specific gene expression levels. The paradigmatic application of this approach to coronary artery disease or schizophrenia patient cohorts identified diverse strata or biotypes. These biotypes are characterized by distinct endophenotype profiles as well as clinical parameters and are fundamentally distinct from PRS based groupings. In stark contrast to the latter, the CASTom-iGEx strategy discovers biologically meaningful and clinically actionable patient subgroups, where complex genetic liabilities are not randomly distributed across individuals but rather converge onto distinct disease relevant biological processes. These results support the notion of different patient biotypes characterized by partially distinct pathomechanisms. Thus, the universally applicable approach presented here has the potential to constitute an important component of future personalized medicine paradigms.

lncRNA CDKN2B-AS1 regulates collagen expression

The long noncoding RNA CDKN2B-AS1 harbors a major coronary artery disease risk haplotype, which is also associated with progressive forms of the oral inflammatory disease periodontitis as well as myocardial infarction (MI). Despite extensive research, there is currently no broad consensus on the function of CDKN2B-AS1 that would explain a common molecular role of this lncRNA in these diseases. Our aim was to investigate the role of CDKN2B-AS1 in gingival cells to better understand the molecular mechanisms underlying the increased risk of progressive periodontitis. We downregulated CDKN2B-AS1 transcript levels in primary gingival fibroblasts with LNA GapmeRs. Following RNA-sequencing, we performed differential expression, gene set enrichment analyses and Western Blotting. Putative causal alleles were searched by analyzing associated DNA sequence variants for changes of predicted transcription factor binding sites. We functionally characterized putative functional alleles using luciferase-reporter and antibody electrophoretic mobility shift assays in gingival fibroblasts and HeLa cells. Of all gene sets analysed, collagen biosynthesis was most significantly upregulated (Padj=9.7 × 10- 5 (AUC > 0.65) with the CAD and MI risk gene COL4A1 showing strongest upregulation of the enriched gene sets (Fold change = 12.13, Padj = 4.9 × 10- 25). The inflammatory "TNFA signaling via NFKB" gene set was downregulated the most (Padj=1 × 10- 5 (AUC = 0.60). On the single gene level, CAPNS2, involved in extracellular matrix organization, was the top upregulated protein coding gene (Fold change = 48.5, P < 9 × 10- 24). The risk variant rs10757278 altered a binding site of the pathogen responsive transcription factor STAT1 (P = 5.8 × 10- 6). rs10757278-G allele reduced STAT1 binding 14.4% and rs10757278-A decreased luciferase activity in gingival fibroblasts 41.2% (P = 0.0056), corresponding with GTEx data. CDKN2B-AS1 represses collagen gene expression in gingival fibroblasts. Dysregulated collagen biosynthesis through allele-specific CDKN2B-AS1 expression in response to inflammatory factors may affect collagen synthesis, and in consequence tissue barrier and atherosclerotic plaque stability.

Exome sequence analysis identifies rare coding variants associated with a machine learning-based marker for coronary artery disease

Coronary artery disease (CAD) exists on a spectrum of disease represented by a combination of risk factors and pathogenic processes. An in silico score for CAD built using machine learning and clinical data in electronic health records captures disease progression, severity and underdiagnosis on this spectrum and could enhance genetic discovery efforts for CAD. Here we tested associations of rare and ultrarare coding variants with the in silico score for CAD in the UK Biobank, All of Us Research Program and BioMe Biobank. We identified associations in 17 genes; of these, 14 show at least moderate levels of prior genetic, biological and/or clinical support for CAD. We also observed an excess of ultrarare coding variants in 321 aggregated CAD genes, suggesting more ultrarare variant associations await discovery. These results expand our understanding of the genetic etiology of CAD and illustrate how digital markers can enhance genetic association investigations for complex diseases.

Atherosclerosis Is a Smooth Muscle Cell-Driven Tumor-Like Disease

BACKGROUND

Atherosclerosis, a leading cause of cardiovascular disease, involves the pathological activation of various cell types, including immunocytes (eg, macrophages and T cells), smooth muscle cells (SMCs), and endothelial cells. Accumulating evidence suggests that transition of SMCs to other cell types, known as phenotypic switching, plays a central role in atherosclerosis development and complications. However, the characteristics of SMC-derived cells and the underlying mechanisms of SMC transition in disease pathogenesis remain poorly understood. Our objective is to characterize tumor cell-like behaviors of SMC-derived cells in atherosclerosis, with the ultimate goal of developing interventions targeting SMC transition for the prevention and treatment of atherosclerosis.

METHODS

We used SMC lineage tracing mice and human tissues and applied a range of methods, including molecular, cellular, histological, computational, human genetics, and pharmacological approaches, to investigate the features of SMC-derived cells in atherosclerosis.

RESULTS

SMC-derived cells in mouse and human atherosclerosis exhibit multiple tumor cell-like characteristics, including genomic instability, evasion of senescence, hyperproliferation, resistance to cell death, invasiveness, and activation of comprehensive cancer-associated gene regulatory networks. Specific expression of the oncogenic mutant KrasG12D in SMCs accelerates phenotypic switching and exacerbates atherosclerosis. Furthermore, we provide proof of concept that niraparib, an anticancer drug targeting DNA damage repair, attenuates atherosclerosis progression and induces regression of lesions in advanced disease in mouse models.

CONCLUSIONS

Our findings demonstrate that atherosclerosis is an SMC-driven tumor-like disease, advancing our understanding of its pathogenesis and opening prospects for innovative precision molecular strategies aimed at preventing and treating atherosclerotic cardiovascular disease.

Sex and sex steroids as determinants of cardiovascular risk

There are considerable sex differences regarding the risk of cardiovascular disease (CVD), including arterial hypertension, coronary artery disease (CAD) and stroke, as well as chronic renal disease. Women are largely protected from these conditions prior to menopause, and the risk increases following cessation of endogenous estrogen production or after surgical menopause. Cardiovascular diseases in women generally begin to occur at a later age than in men (on average with a delay of 10 years). Cessation of estrogen production also impacts metabolism, increasing the risk of developing obesity and diabetes. In middle-aged individuals, hypertension develops earlier and faster in women than in men, and smoking increases cardiovascular risk to a greater degree in women than it does in men. It is not only estrogen that affects female cardiovascular health and plays a protective role until menopause: other sex hormones such as progesterone and androgen hormones generate a complex balance that differentiates heart and blood vessel function in women compared to men. Estrogens improve vasodilation of epicardial coronary arteries and the coronary microvasculature by augmenting the release of vasodilating factors such as nitric oxide and prostacyclin, which are mechanisms of coronary vasodilatation that are more pronounced in women compared to men. Estrogens are also powerful inhibitors of inflammation, which in part explains their protective effects on CVD and chronic renal disease. Emerging evidence suggests that sex chromosomes also play a significant role in shaping cardiovascular risk. The cardiovascular protection conferred by endogenous estrogens may be extended by hormone therapy, especially using bioidentical hormones and starting treatment early after menopause.

Coronary Artery Disease Risk Variant Dampens the Expression of CALCRL by Reducing HSF Binding to Shear Stress Responsive Enhancer in Endothelial Cells In Vitro

BACKGROUND

CALCRL (calcitonin receptor-like) protein is an important mediator of the endothelial fluid shear stress response, which is associated with the genetic risk of coronary artery disease. In this study, we functionally characterized the noncoding regulatory elements carrying coronary artery disease that risks single-nucleotide polymorphisms and studied their role in the regulation of CALCRL expression in endothelial cells.

METHODS

To functionally characterize the coronary artery disease single-nucleotide polymorphisms harbored around the gene CALCRL, we applied an integrative approach encompassing statistical, transcriptional (RNA-seq), and epigenetic (ATAC-seq [transposase-accessible chromatin with sequencing], chromatin immunoprecipitation assay-quantitative polymerase chain reaction, and electromobility shift assay) analyses, alongside luciferase reporter assays, and targeted gene and enhancer perturbations (siRNA and clustered regularly interspaced short palindromic repeats/clustered regularly interspaced short palindromic repeat-associated 9) in human aortic endothelial cells.

RESULTS

We demonstrate that the regulatory element harboring rs880890 exhibits high enhancer activity and shows significant allelic bias. The A allele was favored over the G allele, particularly under shear stress conditions, mediated through alterations in the HSF1 (heat shock factor 1) motif and binding. CRISPR deletion of rs880890 enhancer resulted in downregulation of CALCRL expression, whereas HSF1 knockdown resulted in a significant decrease in rs880890-enhancer activity and CALCRL expression. A significant decrease in HSF1 binding to the enhancer region in endothelial cells was observed under disturbed flow compared with unidirectional flow. CALCRL knockdown and variant perturbation experiments indicated the role of CALCRL in mediating eNOS (endothelial nitric oxide synthase), APLN (apelin), angiopoietin, prostaglandins, and EDN1 (endothelin-1) signaling pathways leading to a decrease in cell proliferation, tube formation, and NO production.

CONCLUSIONS

Overall, our results demonstrate the existence of an endothelial-specific HSF (heat shock factor)-regulated transcriptional enhancer that mediates CALCRL expression. A better understanding of CALCRL gene regulation and the role of single-nucleotide polymorphisms in the modulation of CALCRL expression could provide important steps toward understanding the genetic regulation of shear stress signaling responses.

MRAS in coronary artery disease-Unchartered territory

Genome-wide association studies (GWAS) have identified coronary artery disease (CAD) susceptibility locus on chromosome 3q22.3. This locus contains a cluster of several genes that includes muscle rat sarcoma virus (MRAS). Common MRAS variants are also associated with CAD causing risk factors such as hypertension, dyslipidemia, obesity, and type II diabetes. The MRAS gene is an oncogene that encodes a membrane-bound small GTPase. It is involved in a variety of signaling pathways, regulating cell differentiation and cell survival (mitogen-activated protein kinase [MAPK]/extracellular signal-regulated kinase and phosphatidylinositol 3-kinase) as well as acute phase response signaling (tumor necrosis factor [TNF] and interleukin 6 [IL6] signaling). In this review, we will summarize the role of genetic MRAS variants in the etiology of CAD and its comorbidities with the focus on tissue distribution of MRAS isoforms, cell type/tissue specificity, and mode of action of single nucleotide variants in MRAS associated complex traits. Finally, we postulate that CAD risk variants in the MRAS locus are specific to smooth muscle cells and lead to higher levels of MRAS, particularly in arterial and cardiac tissue, resulting in MAPK-dependent tissue hypertrophy or hyperplasia.

Multi-Ancestry Polygenic Risk Score for Coronary Heart Disease Based on an Ancestrally Diverse Genome-Wide Association Study and Population-Specific Optimization

BACKGROUND

Predictive performance of polygenic risk scores (PRS) varies across populations. To facilitate equitable clinical use, we developed PRS for coronary heart disease (CHD; PRSCHD) for 5 genetic ancestry groups.

METHODS

We derived ancestry-specific and multi-ancestry PRSCHD based on pruning and thresholding (PRSPT) and ancestry-based continuous shrinkage priors (PRSCSx) applied to summary statistics from the largest multi-ancestry genome-wide association study meta-analysis for CHD to date, including 1.1 million participants from 5 major genetic ancestry groups. Following training and optimization in the Million Veteran Program, we evaluated the best-performing PRSCHD in 176,988 individuals across 9 diverse cohorts.

RESULTS

Multi-ancestry PRSPT and PRSCSx outperformed ancestry-specific PRSPT and PRSCSx across a range of tuning values. Two best-performing multi-ancestry PRSCHD (ie, PRSPTmult and PRSCSxmult) and 1 ancestry-specific (PRSCSxEUR) were taken forward for validation. PRSPTmult demonstrated the strongest association with CHD in individuals of South Asian ancestry and European ancestry (odds ratio per 1 SD [95% CI, 2.75 [2.41-3.14], 1.65 [1.59-1.72]), followed by East Asian ancestry (1.56 [1.50-1.61]), Hispanic/Latino ancestry (1.38 [1.24-1.54]), and African ancestry (1.16 [1.11-1.21]). PRSCSxmult showed the strongest associations in South Asian ancestry (2.67 [2.38-3.00]) and European ancestry (1.65 [1.59-1.71]), lower in East Asian ancestry (1.59 [1.54-1.64]), Hispanic/Latino ancestry (1.51 [1.35-1.69]), and the lowest in African ancestry (1.20 [1.15-1.26]).

CONCLUSIONS

The use of summary statistics from a large multi-ancestry genome-wide meta-analysis improved the performance of PRSCHD in most ancestry groups compared with single-ancestry methods. Despite the use of one of the largest and most diverse sets of training and validation cohorts to date, improvement of predictive performance was limited in African ancestry. This highlights the need for larger genome-wide association study datasets of underrepresented populations to enhance the performance of PRSCHD.

Bangladeshi Atherosclerosis Biobank and Hub: The BANGABANDHU Study

Background

Genetic factors contribute significantly to the risk of ischaemic heart disease (IHD), which is the leading cause of mortality in Bangladesh. The BANGABANDHU (Bangladeshi Atherosclerosis Biobank AND Hub) study will allow a hypothesis-free genome-wide association study (GWAS) to identify genetic risk factors associated with ischaemic heart disease patients undergoing coronary artery bypass graft (CABG) surgery in Bangladesh.

Methods

 This is a multi-centre population-based case-control study aimed to evaluate 1500 (Fifteen Hundred) adult (≥18 years of age) people divided into 2 study groups: Case/Proband (750 IHD patients undergoing CABG surgery) and Control (750 healthy people). Spouses or family members are preferred as healthy control subjects due to their shared geographic location and similar environmental exposure.

Results

This will be the first largest DNA repository of CABG patients in Bangladesh, and identifying novel gene loci among CABG patients might help to discover novel therapeutic targets for Bangladeshi IHD patients. Further, identifying and comparing novel gene loci among CABG patients with other ancestry might help devise national guidelines for treating coronary artery disease.

Conclusion

Promising current study results will encourage Bangladeshi researchers and pharmaceutical companies to conduct further studies into the genetic basis of Bangladeshi complex coronary artery disease, which might identify novel genes for therapeutic targets for Bangladeshi patients and strengthen the healthcare standards in Bangladesh.

Enhancing prediction accuracy of coronary artery disease through machine learning-driven genomic variant selection

Machine learning (ML) methods are increasingly becoming crucial in genome-wide association studies for identifying key genetic variants or SNPs that statistical methods might overlook. Statistical methods predominantly identify SNPs with notable effect sizes by conducting association tests on individual genetic variants, one at a time, to determine their relationship with the target phenotype. These genetic variants are then used to create polygenic risk scores (PRSs), estimating an individual's genetic risk for complex diseases like cancer or cardiovascular disorders. Unlike traditional methods, ML algorithms can identify groups of low-risk genetic variants that improve prediction accuracy when combined in a mathematical model. However, the application of ML strategies requires addressing the feature selection challenge to prevent overfitting. Moreover, ensuring the ML model depends on a concise set of genomic variants enhances its clinical applicability, where testing is feasible for only a limited number of SNPs. In this study, we introduce a robust pipeline that applies ML algorithms in combination with feature selection (ML-FS algorithms), aimed at identifying the most significant genomic variants associated with the coronary artery disease (CAD) phenotype. The proposed computational approach was tested on individuals from the UK Biobank, differentiating between CAD and non-CAD individuals within this extensive cohort, and benchmarked against standard PRS-based methodologies like LDpred2 and Lassosum. Our strategy incorporates cross-validation to ensure a more robust evaluation of genomic variant-based prediction models. This method is commonly applied in machine learning strategies but has often been neglected in previous studies assessing the predictive performance of polygenic risk scores. Our results demonstrate that the ML-FS algorithm can identify panels with as few as 50 genetic markers that can achieve approximately 80% accuracy when used in combination with known risk factors. The modest increase in accuracy over PRS performances is noteworthy, especially considering that PRS models incorporate a substantially larger number of genetic variants. This extensive variant selection can pose practical challenges in clinical settings. Additionally, the proposed approach revealed novel CAD-genetic variant associations.

Single-cell 'omic profiles of human aortic endothelial cells in vitro and human atherosclerotic lesions ex vivo reveal heterogeneity of endothelial subtype and response to activating perturbations

Heterogeneity in endothelial cell (EC) sub-phenotypes is becoming increasingly appreciated in atherosclerosis progression. Still, studies quantifying EC heterogeneity across whole transcriptomes and epigenomes in both in vitro and in vivo models are lacking. Multiomic profiling concurrently measuring transcriptomes and accessible chromatin in the same single cells was performed on six distinct primary cultures of human aortic ECs (HAECs) exposed to activating environments characteristic of the atherosclerotic microenvironment in vitro. Meta-analysis of single-cell transcriptomes across 17 human ex vivo arterial specimens was performed and two computational approaches quantitatively evaluated the similarity in molecular profiles between heterogeneous in vitro and ex vivo cell profiles. HAEC cultures were reproducibly populated by four major clusters with distinct pathway enrichment profiles and modest heterogeneous responses: EC1-angiogenic, EC2-proliferative, EC3-activated/mesenchymal-like, and EC4-mesenchymal. Quantitative comparisons between in vitro and ex vivo transcriptomes confirmed EC1 and EC2 as most canonically EC-like, and EC4 as most mesenchymal with minimal effects elicited by siERG and IL1B. Lastly, accessible chromatin regions unique to EC2 and EC4 were most enriched for coronary artery disease (CAD)-associated single-nucleotide polymorphisms from Genome Wide Association Studies (GWAS), suggesting that these cell phenotypes harbor CAD-modulating mechanisms. Primary EC cultures contain markedly heterogeneous cell subtypes defined by their molecular profiles. Surprisingly, the perturbations used here only modestly shifted cells between subpopulations, suggesting relatively stable molecular phenotypes in culture. Identifying consistently heterogeneous EC subpopulations between in vitro and ex vivo models should pave the way for improving in vitro systems while enabling the mechanisms governing heterogeneous cell state decisions.

Advancements in Omics and Breakthrough Gene Therapies: A Glimpse into the Future of Peripheral Artery Disease

Peripheral artery disease (PAD), a highly prevalent global disease, associates with significant morbidity and mortality in affected patients. Despite progress in endovascular and open revascularization techniques for advanced PAD, these interventions grapple with elevated rates of arterial restenosis and vein graft failure attributed to intimal hyperplasia (IH). Novel multiomics technologies, coupled with sophisticated analyses tools recently powered by advances in artificial intelligence, have enabled the study of atherosclerosis and IH with unprecedented single-cell and spatial precision. Numerous studies have pinpointed gene hubs regulating pivotal atherogenic and atheroprotective signaling pathways as potential therapeutic candidates. Leveraging advancements in viral and nonviral gene therapy (GT) platforms, gene editing technologies, and cutting-edge biomaterial reservoirs for delivery uniquely positions us to develop safe, efficient, and targeted GTs for PAD-related diseases. Gene therapies appear particularly fitting for ex vivo genetic engineering of IH-resistant vein grafts. This manuscript highlights currently available state-of-the-art multiomics approaches, explores promising GT-based candidates, and details GT delivery modalities employed by our laboratory and others to thwart mid-term vein graft failure caused by IH, as well as other PAD-related conditions. The potential clinical translation of these targeted GTs holds the promise to revolutionize PAD treatment, thereby enhancing patients' quality of life and life expectancy.

Transforming the cardiometabolic disease landscape: Multimodal AI-powered approaches in prevention and management

The rise of artificial intelligence (AI) has revolutionized various scientific fields, particularly in medicine, where it has enabled the modeling of complex relationships from massive datasets. Initially, AI algorithms focused on improved interpretation of diagnostic studies such as chest X-rays and electrocardiograms in addition to predicting patient outcomes and future disease onset. However, AI has evolved with the introduction of transformer models, allowing analysis of the diverse, multimodal data sources existing in medicine today. Multimodal AI holds great promise in more accurate disease risk assessment and stratification as well as optimizing the key driving factors in cardiometabolic disease: blood pressure, sleep, stress, glucose control, weight, nutrition, and physical activity. In this article we outline the current state of medical AI in cardiometabolic disease, highlighting the potential of multimodal AI to augment personalized prevention and treatment strategies in cardiometabolic disease.

The Role of Endothelial Cells in Atherosclerosis: Insights from Genetic Association Studies

Endothelial cells (ECs) mediate several biological functions that are relevant to atherosclerosis and coronary artery disease (CAD), regulating an array of vital processes including vascular tone, wound healing, reactive oxygen species, shear stress response, and inflammation. Although which of these functions is linked causally with CAD development and/or progression is not yet known, genome-wide association studies have implicated more than 400 loci associated with CAD risk, among which several have shown EC-relevant functions. Given the arduous process of mechanistically interrogating single loci to CAD, high-throughput variant characterization methods, including pooled Clustered Regularly Interspaced Short Palindromic Repeats screens, offer exciting potential to rapidly accelerate the discovery of bona fide EC-relevant genetic loci. These discoveries in turn will broaden the therapeutic avenues for CAD beyond lipid lowering and behavioral risk modification to include EC-centric modalities of risk prevention and treatment.

Nuclear Control of Vascular Smooth Muscle Cell Plasticity during Vascular Remodeling

Control of vascular smooth muscle cell (SMC) gene expression is an essential process for establishing and maintaining lineage identity, contractility, and plasticity. Most mechanisms (epigenetic, transcriptional, and post-transcriptional) implicated in gene regulation occur in the nucleus. Still, intranuclear pathways are directly impacted by modifications in the extracellular environment in conditions of adaptive or maladaptive remodeling. Integration of extracellular, cellular, and genomic information into the nucleus through epigenetic and transcriptional control of genome organization plays a major role in regulating SMC functions and phenotypic transitions during vascular remodeling and diseases. This review aims to provide a comprehensive update on nuclear mechanisms, their interactions, and their integration in controlling SMC homeostasis and dysfunction. It summarizes and discusses the main nuclear mechanisms preponderant in SMCs in the context of vascular disease, such as atherosclerosis, with an emphasis on studies employing in vivo cell-specific loss-of-function and single-cell omics approaches.

Omics Approaches Unveiling the Biology of Human Atherosclerotic Plaques

Atherosclerosis is a chronic inflammatory disease of the arterial wall, characterized by the buildup of plaques with the accumulation and transformation of lipids, immune cells, vascular smooth muscle cells, and necrotic cell debris. Plaques with collagen-poor thin fibrous caps infiltrated by macrophages and lymphocytes are considered unstable because they are at the greatest risk of rupture and clinical events. However, the current histologic definition of plaque types may not fully capture the complex molecular nature of atherosclerotic plaque biology and the underlying mechanisms contributing to plaque progression, rupture, and erosion. The advances in omics technologies have changed the understanding of atherosclerosis plaque biology, offering new possibilities to improve risk prediction and discover novel therapeutic targets. Genomic studies have shed light on the genetic predisposition to atherosclerosis, and integrative genomic analyses expedite the translation of genomic discoveries. Transcriptomic, proteomic, metabolomic, and lipidomic studies have refined the understanding of the molecular signature of atherosclerotic plaques, aiding in data-driven hypothesis generation for mechanistic studies and offering new prospects for biomarker discovery. Furthermore, advancements in single-cell technologies and emerging spatial analysis techniques have unveiled the heterogeneity and plasticity of plaque cells. This review discusses key omics-based discoveries that have advanced the understanding of human atherosclerotic plaque biology, focusing on insights derived from omics profiling of human atherosclerotic vascular specimens.

Somatic and Germline Variants and Coronary Heart Disease in a Chinese Population

Importance

The genetic basis of coronary heart disease (CHD) has expanded from a germline to somatic genome, including clonal hematopoiesis of indeterminate potential (CHIP). How CHIP confers CHD risk in East Asian individuals, especially those with small clones (variant allele fraction [VAF] 0.5%-2%) and different genetic backgrounds, was completely unknown.

Objective

To investigate the CHIP profile in a general Chinese cohort by deep sequencing and further explore the association between CHIP and incident CHD considering germline predisposition.

Design, Setting, and Participants

This cohort study used data from 3 prospective cohorts in the project Prediction for Atherosclerotic Cardiovascular Disease Risk in China. Participants without cardiovascular disease or cancer at baseline were enrolled in 2001 and 2008 and had a median follow-up of 12.17 years extending into 2021.

Exposures

CHIP mutations were detected by targeted sequencing (mean depth, 916×). A predefined CHD polygenic risk score (PRS) comprising 531 variants was used to evaluate germline predisposition.

Main Outcomes and Measures

The main outcome was first incident CHD.

Results

Among 6181 participants, the median (IQR) age was 53.83 years (45.35-62.39 years); 3082 participants (49.9%) were female, and 3099 (50.1%) were male. A total of 1100 individuals (17.80%) harbored 1372 CHIP mutations at baseline. CHIP was independently associated with incident CHD (hazard ratio [HR], 1.42; 95% CI, 1.18-1.72; P = 2.82 × 10-4) and presented a risk gradient with increasing VAF (P = 3.98 × 10-3 for trend). Notably, individuals with small clones, nearly half of CHIP carriers, also demonstrated a higher CHD risk compared with non-CHIP carriers (HR, 1.33; 95% CI, 1.02-1.74; P = .03) and were 4 years younger than those with VAF of 2% or greater (median age, 58.52 vs 62.70 years). Heightened CHD risk was not observed among CHIP carriers with low PRS (HR, 1.02; 95% CI, 0.64-1.64; P = .92), while high PRS and CHIP jointly contributed a 2.23-fold increase in risk (95% CI, 1.51-3.29; P = 6.29 × 10-5) compared with non-CHIP carriers with low PRS. Interestingly, the diversity in CHIP-related CHD risk within each PRS group was substantially diminished when removing variants in the inflammatory pathway from the PRS.

Conclusions

This study revealed that elevated CHD risk attributed to CHIP was nonnegligible even for small clones. Inflammation genes involved in CHD could aggravate or abrogate CHIP-related CHD risk.

The significance of PCSK-9's level and polymorphism in premature coronary artery disease: Relation to risk and severity

BACKGROUND

Proprotein convertase subtilisin/kexin type 9 (PCSK-9) is a circulating protein that plays an important role in lipid metabolism and is linked to inflammation, which has implications for atherosclerosis and its severe cardiac effects. We studied the potential association of the PCSK-9 gene single nucleotide polymorphism (SNP), Oxidized low-density lipoprotein receptor 1- (OLR-1), and caspase-3 serum levels with the risk and severity of premature coronary artery disease (PCAD). The potential contribution of PCSK-9 serum level to the severity of PCAD patients was also assessed.

METHOD

This case-control study included 120 PCAD patients (age < 45), and 60 age matched healthy controls. Serum PCSK-9 and caspase-3 levels and clinical characteristics were recorded. SYNTAX score was calculated to estimate the severity of the coronary artery lesions. The SNP rs2483205 of the PCSK-9 gene and the rs11053646 of the OLR-1gene were genotyped in all participants.

RESULTS

Serum PCSK-9 levels were higher in PCAD patients and were significantly different among the three SYNTAX score groups (SS ≤ 12, 12 < SS ≤ 21.5, and SS > 21.5). The diagnostic cutoff values of PCSK-9 and caspase-3 levels for PCAD were > 3.2 ng/mL for both, yielding an area under the curve (AUC) of 0.98 and 0.92, sensitivity of 85 %, 98 %, and specificity of 99.5 %, 93 % for PCSK-9 and caspase-3, respectively. The genotypes TT + CT vs. CC of PCSK-9's rs2483205 SNP presented a higher risk for PCAD and higher SYNTAX scores. Furthermore, the rs11053646 SNP of OLR-1 presented the CG genotype as more risky and having higher SYNTAX scores.

CONCLUSION

Circulating PCSK9 and caspase-3 concentrations were higher in PCAD patients and were associated with CAD severity. The SNPs of PCSK-9 (rs2483205) and OLR-1 (rs11053646) were associated with PCAD and its severity.

Association of Tenascin-C Gene Polymorphisms with Risk of Acute Coronary Syndrome in South Indian Population: A Case-Control Genetic Association Study

Background: The extracellular matrix (ECM) glycoprotein changes are associated with the pathogenesis and complications of atherosclerosis, leading to acute coronary syndrome (ACS). Tenascin-C (TNC), an ECM protein, has been implemented in the pathogenesis, diagnosis, and prognosis of patients with cardiovascular disease. Aim: The study aimed to compare the genetic variants of the TNC gene (rs13321, rs2104772, and rs12347433) between South Indians with ACS and healthy participants. Materials and Methods: This case-control study recruited 150 ACS patients as cases and 150 healthy participants as controls. TNC genotyping was performed using TaqMan 5'-exonuclease allele discrimination assay. Serum TNC levels were measured by enzyme-linked immunosorbent assay. Results: Serum TNC levels were significantly higher in cases compared with controls. No significant difference was observed in allele and genotype frequencies of rs13321, rs2104772, and rs12347433 between cases and controls, which was confirmed by dominant, recessive, codominant, and homozygotic genetic models. The patients with heterozygous genotypes of rs13321, rs2104772, and rs12347433 had significantly lower serum TNC levels than patients with respective homozygous genotypes. Haplotype analyses revealed that the C-T-A haplotype in the block of rs13321-rs12347433-rs2104772 was associated with lower ACS risk (OR = 0.33, 95% CI: 0.15 - 0.75; p = 0.005). Also, the C-T-T and G-T-A haplotypes of the TNC gene were associated with higher and lower serum TNC levels, respectively. Conclusion: Our study demonstrated no genetic association between single nucleotide polymorphisms of the TNC gene and ACS risk; however, the C-T-A haplotype of the TNC gene might be associated with reduced ACS risk in South Indians.

ADAMTS7 Promotes Smooth Muscle Cell Foam Cell Expansion in Atherosclerosis

Human genetic studies have repeatedly associated SNPs near the gene ADAMTS7 with atherosclerotic cardiovascular disease. Subsequent investigations in mice demonstrated that ADAMTS7 is proatherogenic, induced in response to vascular injury, and alters smooth muscle cell function. However, the mechanisms governing this function and its relationship to atherosclerosis remain unclear. Here, we report the first conditional Adamts7 transgenic mouse in which the gene can be conditionally overexpressed in smooth muscle cells, mimicking its induction in atherosclerosis. We observed that smooth muscle cell Adamts7 overexpression results in a 3.5-fold increase in peripheral atherosclerosis, coinciding with an expansion of smooth muscle foam cells. RNA sequencing of Adamts7 overexpressed primary smooth muscle cells revealed an upregulation in the expression of lipid uptake genes. Subsequent experiments in primary smooth muscle cells demonstrated that increased Spi1 and Cd36 expression leads to increased smooth muscle cell oxLDL uptake. To uncover ADAMTS7 expression in human disease, we have interrogated the largest scRNA-seq dataset of human carotid atherosclerosis. This analysis discovered that endothelial cells had the highest expression level of ADAMTS7 with lesser expression in smooth muscle cells, fibroblasts, and mast cells. Subsequent conditional knockout studies in smooth muscle cells surprisingly showed no change in atherosclerosis, suggesting redundant expression of this secreted factor in the vessel wall. Finally, mice overexpressing Adamts7 in endothelial cells also exhibit increased atherosclerosis, suggesting that multiple vascular cell types can contribute to ADAMTS7-mediated foam cell expansion. In summary, Adamts7 is expressed by multiple vascular cell types in atherosclerosis, and ADAMTS7 promotes oxLDL uptake in smooth muscle cells, increasing smooth muscle foam cell formation and peripheral atherosclerosis in mice.

EGR1 transcriptionally regulates SVEP1 to promote proliferation and migration in human coronary artery smooth muscle cells

A low-frequency variant of sushi, von Willebrand factor type A, EGF, and pentraxin domain-containing protein 1 (SVEP1) is associated with the risk of coronary artery disease, as determined by a genome-wide association study. SVEP1 induces vascular smooth muscle cell proliferation and an inflammatory phenotype to promote atherosclerosis. In the present study, qRT‒PCR demonstrated that the mRNA expression of SVEP1 was significantly increased in atherosclerotic plaques compared to normal tissues. Bioinformatics revealed that EGR1 was a transcription factor for SVEP1. The results of the luciferase reporter assay, siRNA interference or overexpression assay, mutational analysis and ChIP confirmed that EGR1 positively regulated the transcriptional activity of SVEP1 by directly binding to its promoter. EGR1 promoted human coronary artery smooth muscle cell (HCASMC) proliferation and migration via SVEP1 in response to oxidized low-density lipoprotein (ox-LDL) treatment. Moreover, the expression level of EGR1 was increased in atherosclerotic plaques and showed a strong linear correlation with the expression of SVEP1. Our findings indicated that EGR1 binding to the promoter region drive SVEP1 transcription to promote HCASMC proliferation and migration.

Polygenic risk score predicts all-cause death in East Asian patients with prior coronary artery disease

Introduction

Coronary artery disease (CAD) is a highly heritable and multifactorial disease. Numerous genome-wide association studies (GWAS) facilitated the construction of polygenic risk scores (PRS) for predicting future incidence of CAD, however, exclusively in European populations. Furthermore, identifying CAD patients with elevated risks of all-cause death presents a critical challenge in secondary prevention, which will contribute largely to reducing the burden for public healthcare.

Methods

We recruited a cohort of 1,776 Chinese CAD patients and performed medical follow-up for up to 11 years. A pruning and thresholding method was used to calculate PRS of CAD and its 14 risk factors. Their correlations with all-cause death were computed via Cox regression.

Results and discussion

We found that the PRS for CAD and its seven risk factors, namely myocardial infarction, ischemic stroke, angina, heart failure, low-density lipoprotein cholesterol, total cholesterol and C-reaction protein, were significantly associated with death (P ≤ 0.05), whereas the PRS of body mass index displayed moderate association (P < 0.1). Elastic-net Cox regression with 5-fold cross-validation was used to integrate these nine PRS models into a meta score, metaPRS, which performed well in stratifying patients at different risks for death (P < 0.0001). Combining metaPRS with clinical risk factors further increased the discerning power and a 4% increase in sensitivity. The metaPRS generated from the genetic susceptibility to CAD and its risk factors can well stratify CAD patients by their risks of death. Integrating metaPRS and clinical risk factors may contribute to identifying patients at higher risk of poor prognosis.

The Inclusion of Underrepresented Populations in Cardiovascular Genetics and Epidemiology

Novel genetic risk markers have helped us to advance the field of cardiovascular epidemiology and refine our current understanding and risk stratification paradigms. The discovery and analysis of variants can help us to tailor prognostication and management. However, populations underrepresented in cardiovascular epidemiology and cardiogenetics research may experience inequities in care if prediction tools are not applicable to them clinically. Therefore, the purpose of this article is to outline the barriers that underrepresented populations can face in participating in genetics research, to describe the current efforts to diversify cardiogenetics research, and to outline strategies that researchers in cardiovascular epidemiology can implement to include underrepresented populations. Mistrust, a lack of diverse research teams, the improper use of sensitive biodata, and the constraints of genetic analyses are all barriers for including diverse populations in genetics studies. The current work is beginning to address the paucity of ethnically diverse genetics research and has already begun to shed light on the potential benefits of including underrepresented and diverse populations. Reducing barriers for individuals, utilizing community-driven research processes, adopting novel recruitment strategies, and pushing for organizational support for diverse genetics research are key steps that clinicians and researchers can take to develop equitable risk stratification tools and improve patient care.

Effects of Ninjurin 2 polymorphisms on susceptibility to coronary heart disease

OBJECTIVE

The aim of this study was to explore the effects of Ninjurin 2 (NINJ2) polymorphisms on susceptibility to coronary heart disease (CHD).

METHODS

We conducted a case-control study with 499 CHD cases and 505 age and gender-matched controls. Five single nucleotide polymorphisms (SNPs) in NINJ2 (rs118050317, rs75750647, rs7307242, rs10849390, and rs11610368) were genotyped by the Agena MassARRAY platform. Odds ratios (ORs) and 95% confidence intervals (CIs) were calculated using logistic regression analysis to assess the association of NINJ2 polymorphisms and CHD risk-adjusted for age and gender. What's more, risk genes and molecular functions were screened via protein-protein interaction (PPI) network and functional enrichment analysis.

RESULTS

Rs118050317 in NINJ2 significantly increased CHD risk in people aged more than 60 years and women. Rs118050317 and rs7307242 had strong relationships with hypertension risk in CHD patients. Additionally, rs75750647 exceedingly raised diabetes risk in cases under multiple models, whereas rs10849390 could protect CHD patients from diabetes in allele, homozygote, and additive models. We also observed two blocks in NINJ2. Further interaction network and enrichment analysis showed that NINJ2 played a greater role in the pathogenesis and progression of CHD.

CONCLUSION

Our results suggest that NINJ2 polymorphisms are associated with CHD risk.

Convergence of coronary artery disease genes onto endothelial cell programs

Linking variants from genome-wide association studies (GWAS) to underlying mechanisms of disease remains a challenge1-3. For some diseases, a successful strategy has been to look for cases in which multiple GWAS loci contain genes that act in the same biological pathway1-6. However, our knowledge of which genes act in which pathways is incomplete, particularly for cell-type-specific pathways or understudied genes. Here we introduce a method to connect GWAS variants to functions. This method links variants to genes using epigenomics data, links genes to pathways de novo using Perturb-seq and integrates these data to identify convergence of GWAS loci onto pathways. We apply this approach to study the role of endothelial cells in genetic risk for coronary artery disease (CAD), and discover 43 CAD GWAS signals that converge on the cerebral cavernous malformation (CCM) signalling pathway. Two regulators of this pathway, CCM2 and TLNRD1, are each linked to a CAD risk variant, regulate other CAD risk genes and affect atheroprotective processes in endothelial cells. These results suggest a model whereby CAD risk is driven in part by the convergence of causal genes onto a particular transcriptional pathway in endothelial cells. They highlight shared genes between common and rare vascular diseases (CAD and CCM), and identify TLNRD1 as a new, previously uncharacterized member of the CCM signalling pathway. This approach will be widely useful for linking variants to functions for other common polygenic diseases.

Genome-Wide Genetic Associations Prioritize Evaluation of Causal Mechanisms of Atherosclerotic Disease Risk

OBJECTIVE

The goal of this review is to discuss the implementation of genome-wide association studies to identify causal mechanisms of vascular disease risk.

APPROACH AND RESULTS

The history of genome-wide association studies is described, the use of imputation and the creation of consortia to conduct meta-analyses with sufficient power to arrive at consistent associated loci for vascular disease. Genomic methods are described that allow the identification of causal variants and causal genes and how they impact the disease process. The power of single-cell analyses to promote genome-wide association studies of causal gene function is described.

CONCLUSIONS

Genome-wide association studies represent a paradigm shift in the study of cardiovascular disease, providing identification of genes, cellular phenotypes, and disease pathways that empower the future of targeted drug development.

Unlocking the potential of big data and AI in medicine: insights from biobanking

Big data and artificial intelligence are key elements in the medical field as they are expected to improve accuracy and efficiency in diagnosis and treatment, particularly in identifying biomedically relevant patterns, facilitating progress towards individually tailored preventative and therapeutic interventions. These applications belong to current research practice that is data-intensive. While the combination of imaging, pathological, genomic, and clinical data is needed to train algorithms to realize the full potential of these technologies, biobanks often serve as crucial infrastructures for data-sharing and data flows. In this paper, we argue that the 'data turn' in the life sciences has increasingly re-structured major infrastructures, which often were created for biological samples and associated data, as predominantly data infrastructures. These have evolved and diversified over time in terms of tackling relevant issues such as harmonization and standardization, but also consent practices and risk assessment. In line with the datafication, an increased use of AI-based technologies marks the current developments at the forefront of the big data research in life science and medicine that engender new issues and concerns along with opportunities. At a time when secure health data environments, such as European Health Data Space, are in the making, we argue that such meta-infrastructures can benefit both from the experience and evolution of biobanking, but also the current state of affairs in AI in medicine, regarding good governance, the social aspects and practices, as well as critical thinking about data practices, which can contribute to trustworthiness of such meta-infrastructures.

Genetic Predisposition of Atherosclerotic Cardiovascular Disease in Ancient Human Remains

Background

Several computed tomographic studies have shown the presence of atherosclerosis in ancient human remains. However, while it is important to understand the development of atherosclerotic cardiovascular disease (ASCVD), genetic data concerning the prevalence of the disease-associated single nucleotide polymorphisms (SNPs) in our ancestors are scarce.

Objective

For a better understanding of the role of genetics in the evolution of ASCVD, we applied an enrichment capture sequencing approach to mummified human remains from different geographic regions and time periods.

Methods

Twenty-two mummified individuals were analyzed for their genetic predisposition of ASCVD. Next-generation sequencing methods were applied to ancient DNA (aDNA) samples, including a novel enrichment approach specifically designed to capture SNPs associated with ASCVD in genome-wide association studies of modern humans.

Findings

Five out of 22 ancient individuals passed all filter steps for calculating a weighted polygenic risk score (PRS) based on 87 SNPs in 56 genes. PRSs were correlated to scores obtained from contemporary people from around the world and cover their complete range. The genetic results of the ancient individuals reflect their phenotypic results, given that the only two mummies showing calcified atherosclerotic arterial plaques on computed tomography scans are the ones exhibiting the highest calculated PRSs.

Conclusions

These data show that alleles associated with ASCVD have been widespread for at least 5,000 years. Despite some limitations due to the nature of aDNA, our approach has the potential to lead to a better understanding of the interaction between environmental and genetic influences on the development of ASCVD.

Application of polygenic scores to a deeply phenotyped sample enriched for substance use disorders reveals extensive pleiotropy with psychiatric and medical traits

Co-occurring psychiatric, medical, and substance use disorders (SUDs) are common, but the complex pathways leading to such comorbidities are poorly understood. A greater understanding of genetic influences on this phenomenon could inform precision medicine efforts. We used the Yale-Penn dataset, a cross-sectional sample enriched for individuals with SUDs, to examine pleiotropic effects of genetic liability for psychiatric and medical traits. Participants completed an in-depth interview that provides information on demographics, environment, medical illnesses, and psychiatric and SUDs. Polygenic scores (PGS) for psychiatric disorders and medical traits were calculated in European-ancestry (EUR; n=5,691) participants and, when discovery datasets were available, for African-ancestry (AFR; n=4,918) participants. Phenome-wide association studies (PheWAS) were then conducted. In AFR participants, the only PGS with significant associations was bipolar disorder (BD), all of which were with substance use phenotypes. In EUR participants, PGS for major depressive disorder (MDD), generalized anxiety disorder (GAD), post-traumatic stress disorder (PTSD), schizophrenia (SCZ), body mass index (BMI), coronary artery disease (CAD), and type 2 diabetes (T2D) all showed significant associations, the majority of which were with phenotypes in the substance use categories. For instance, PGS MDD was associated with over 200 phenotypes, 15 of which were depression-related (e.g., depression criterion count), 55 of which were other psychiatric phenotypes, and 126 of which were substance use phenotypes; and PGS BMI was associated with 138 phenotypes, 105 of which were substance related. Genetic liability for psychiatric and medical traits is associated with numerous phenotypes across multiple categories, indicative of the broad genetic liability of these traits.

Single cell 'omic profiles of human aortic endothelial cells in vitro and human atherosclerotic lesions ex vivo reveals heterogeneity of endothelial subtype and response to activating perturbations

Objective

Endothelial cells (ECs), macrophages, and vascular smooth muscle cells (VSMCs) are major cell types in atherosclerosis progression, and heterogeneity in EC sub-phenotypes are becoming increasingly appreciated. Still, studies quantifying EC heterogeneity across whole transcriptomes and epigenomes in both in vitro and in vivo models are lacking.

Approach and Results

To create an in vitro dataset to study human EC heterogeneity, multiomic profiling concurrently measuring transcriptomes and accessible chromatin in the same single cells was performed on six distinct primary cultures of human aortic ECs (HAECs). To model pro-inflammatory and activating environments characteristic of the atherosclerotic microenvironment in vitro, HAECs from at least three donors were exposed to three distinct perturbations with their respective controls: transforming growth factor beta-2 (TGFB2), interleukin-1 beta (IL1B), and siRNA-mediated knock-down of the endothelial transcription factor ERG (siERG). To form a comprehensive in vivo/ex vivo dataset of human atherosclerotic cell types, meta-analysis of single cell transcriptomes across 17 human arterial specimens was performed. Two computational approaches quantitatively evaluated the similarity in molecular profiles between heterogeneous in vitro and in vivo cell profiles. HAEC cultures were reproducibly populated by 4 major clusters with distinct pathway enrichment profiles: EC1-angiogenic, EC2-proliferative, EC3-activated/mesenchymal-like, and EC4-mesenchymal. Exposure to siERG, IL1B or TGFB2 elicited mostly distinct transcriptional and accessible chromatin responses. EC1 and EC2, the most canonically 'healthy' EC populations, were affected predominantly by siERG; the activated cluster EC3 was most responsive to IL1B; and the mesenchymal population EC4 was most affected by TGFB2. Quantitative comparisons between in vitro and in vivo transcriptomes confirmed EC1 and EC2 as most canonically EC-like, and EC4 as most mesenchymal with minimal effects elicited by siERG and IL1B. Lastly, accessible chromatin regions unique to EC2 and EC4 were most enriched for coronary artery disease (CAD)-associated SNPs from GWAS, suggesting these cell phenotypes harbor CAD-modulating mechanisms.

Conclusion

Primary EC cultures contain markedly heterogeneous cell subtypes defined by their molecular profiles. Surprisingly, the perturbations used here, which have been reported by others to be involved in the pathogenesis of atherosclerosis as well as induce endothelial-to-mesenchymal transition (EndMT), only modestly shifted cells between subpopulations, suggesting relatively stable molecular phenotypes in culture. Identifying consistently heterogeneous EC subpopulations between in vitro and in vivo models should pave the way for improving in vitro systems while enabling the mechanisms governing heterogeneous cell state decisions.

Multi-ancestry genetic analysis of gene regulation in coronary arteries prioritizes disease risk loci

Genome-wide association studies (GWASs) have identified hundreds of risk loci for coronary artery disease (CAD). However, non-European populations are underrepresented in GWASs, and the causal gene-regulatory mechanisms of these risk loci during atherosclerosis remain unclear. We incorporated local ancestry and haplotypes to identify quantitative trait loci for expression (eQTLs) and splicing (sQTLs) in coronary arteries from 138 ancestrally diverse Americans. Of 2,132 eQTL-associated genes (eGenes), 47% were previously unreported in coronary artery; 19% exhibited cell-type-specific expression. Colocalization revealed subgroups of eGenes unique to CAD and blood pressure GWAS. Fine-mapping highlighted additional eGenes, including TBX20 and IL5. We also identified sQTLs for 1,690 genes, among which TOR1AIP1 and ULK3 sQTLs demonstrated the importance of evaluating splicing to accurately identify disease-relevant isoform expression. Our work provides a patient-derived coronary artery eQTL resource and exemplifies the need for diverse study populations and multifaceted approaches to characterize gene regulation in disease processes.

Genetic ancestry and diagnostic yield of exome sequencing in a diverse population

It has been suggested that diagnostic yield (DY) from Exome Sequencing (ES) may be lower among patients with non-European ancestries than those with European ancestry. We examined the association of DY with estimated continental/subcontinental genetic ancestry in a racially/ethnically diverse pediatric and prenatal clinical cohort. Cases (N = 845) with suspected genetic disorders underwent ES for diagnosis. Continental/subcontinental genetic ancestry proportions were estimated from the ES data. We compared the distribution of genetic ancestries in positive, negative, and inconclusive cases by Kolmogorov-Smirnov tests and linear associations of ancestry with DY by Cochran-Armitage trend tests. We observed no reduction in overall DY associated with any genetic ancestry (African, Native American, East Asian, European, Middle Eastern, South Asian). However, we observed a relative increase in proportion of autosomal recessive homozygous inheritance versus other inheritance patterns associated with Middle Eastern and South Asian ancestry, due to consanguinity. In this empirical study of ES for undiagnosed pediatric and prenatal genetic conditions, genetic ancestry was not associated with the likelihood of a positive diagnosis, supporting the equitable use of ES in diagnosis of previously undiagnosed but potentially Mendelian disorders across all ancestral populations.

Potential Association of The Pathogenic Kruppel-like Factor 14 (KLF14) and Adiponectin (ADIPOQ) SNVs with Susceptibility to T2DM

AIM

To evaluate the associations of the pathogenic variants in Kruppel-like Factor 14 (KLF 14) and Adiponectin (ADIPOQ) with susceptibility to type 2 diabetes mellitus (T2DM).

BACKGROUND

Type 2 diabetes mellitus (T2DM) is a pandemic metabolic disease characterized by increased blood sugar and caused by resistance to insulin in peripheral tissues and damage to pancreatic beta cells. Kruppel-like Factor 14 (KLF-14) is proposed to be a regulator of metabolic diseases, such as diabetes mellitus (DM) and obesity. Adiponectin (ADIPOQ) is an adipocytokine produced by the adipocytes and other tissues and was reported to be involved in T2DM.

OBJECTIVES

To study the possible association of the KLF-14 rs972283 and ADIPOQ-rs266729 with the risk of T2DM in the Saudi population.

METHODS

We have evaluated the association of KLF-14 rs972283 C>T and ADIPOQ-rs266729 C>G SNV with the risk to T2D in the Saudi population using the Amplification Refractory Mutation System PCR (ARMS-PCR), and blood biochemistry analysis. For the KLF-14 rs972283 C>T SNV we included 115 cases and 116 healthy controls, and ADIPOQ-rs266729 C>G SNV, 103 cases and 104 healthy controls were included.

RESULTS

Results indicated that the KLF-14 rs972283 GA genotype and A allele were associated with T2D risk with OR=2.14, p-value= 0.014 and OR=1.99, p-value=0.0003, respectively. Results also ADIPOQ-rs266729 CG genotype and C allele were associated with an elevated T2D risk with an OR=2.53, p=0.003 and OR=1.66, p-value =0.012, respectively.

CONCLUSION

We conclude that SNVs in KLF-14 and ADIPOQ are potential loci for T2D risk. Future large-scale studies to verify these findings are recommended. These results need further verifications in protein functional and large-scale case control studies before being introduced for genetic testing.

Principles and methods for transferring polygenic risk scores across global populations

Polygenic risk scores (PRSs) summarize the genetic predisposition of a complex human trait or disease and may become a valuable tool for advancing precision medicine. However, PRSs that are developed in populations of predominantly European genetic ancestries can increase health disparities due to poor predictive performance in individuals of diverse and complex genetic ancestries. We describe genetic and modifiable risk factors that limit the transferability of PRSs across populations and review the strengths and weaknesses of existing PRS construction methods for diverse ancestries. Developing PRSs that benefit global populations in research and clinical settings provides an opportunity for innovation and is essential for health equity.

Polygenic Risk Score Associates With Atherosclerotic Plaque Characteristics at Autopsy

BACKGROUND

Polygenic risk scores (PRSs) for coronary artery disease (CAD) potentially improve cardiovascular risk prediction. However, their relationship with histopathologic features of CAD has never been examined systematically.

METHODS

From 4327 subjects referred to CVPath by the State of Maryland Office Chief Medical Examiner for sudden death between 1994 and 2015, 2455 cases were randomly selected for genotyping. We generated PRS from 291 known CAD risk loci. Detailed histopathologic examination of the coronary arteries was performed in all subjects. The primary study outcome measurements were histopathologic plaque features determining severity of atherosclerosis, including %stenosis, calcification, thin-cap fibroatheromas, and thrombotic CAD.

RESULTS

After exclusion of cases with insufficient DNA sample quality or with missing data, 954 cases (mean age, 48.8±14.7 years; 75.7% men) remained in the final study cohort. Subjects in the highest PRS quintile exhibited more severe atherosclerosis compared with subjects in the lowest quintile, with greater %stenosis (80.3%±27.0% versus 50.4%±38.7%; adjusted P<0.001) and a higher frequency of calcification (69.6% versus 35.8%; adjusted P=0.004) and thin-cap fibroatheroma (26.7% versus 9.5%; adjusted P=0.007). Even after adjustment for traditional CAD risk factors, subjects within the highest PRS quintile had higher odds of severe atherosclerosis (ie, ≥75% stenosis; adjusted odds ratio, 3.77 [95% CI, 2.10-6.78]; P<0.001) and plaque rupture (adjusted odds ratio, 4.05 [95% CI, 2.26-7.24]; P<0.001). Moreover, subjects within the highest quintile had higher odds of CAD-associated cause of death, especially among those aged ≤50 years (adjusted odds ratio, 4.08 [95% CI, 2.01-8.30]; P<0.001). No statistically significant associations were observed with plaque erosion after adjusting for covariates.

CONCLUSIONS

This is the first autopsy study investigating associations between PRS and atherosclerosis severity at the histopathologic level in subjects with sudden death. Our pathological analysis suggests PRS correlates with plaque burden and features of advanced atherosclerosis and may be useful as a method for CAD risk stratification, especially in younger subjects.

Integrative proteomic analyses across common cardiac diseases yield new mechanistic insights and enhanced prediction

Cardiac diseases represent common highly morbid conditions for which underlying molecular mechanisms remain incompletely understood. Here, we leveraged 1,459 protein measurements in 44,313 UK Biobank participants to characterize the circulating proteome associated with incident coronary artery disease, heart failure, atrial fibrillation, and aortic stenosis. Multivariable-adjusted Cox regression identified 820 protein-disease associations-including 441 proteins-at Bonferroni-adjusted P <8.6×10 -6 . Cis -Mendelian randomization suggested causal roles that aligned with epidemiological findings for 6% of proteins identified in primary analyses, prioritizing novel therapeutic targets for different cardiac diseases (e.g., interleukin-4 receptor for heart failure and spondin-1 for atrial fibrillation). Interaction analyses identified seven protein-disease associations that differed Bonferroni-significantly by sex. Models incorporating proteomic data (vs. clinical risk factors alone) improved prediction for coronary artery disease, heart failure, and atrial fibrillation. These results lay a foundation for future investigations to uncover novel disease mechanisms and assess the clinical utility of protein-based prevention strategies for cardiac diseases.

Shared and distinct pathways and networks genetically linked to coronary artery disease between human and mouse

Mouse models have been used extensively to study human coronary artery disease (CAD) or atherosclerosis and to test therapeutic targets. However, whether mouse and human share similar genetic factors and pathogenic mechanisms of atherosclerosis has not been thoroughly investigated in a data-driven manner. We conducted a cross-species comparison study to better understand atherosclerosis pathogenesis between species by leveraging multiomics data. Specifically, we compared genetically driven and thus CAD-causal gene networks and pathways, by using human GWAS of CAD from the CARDIoGRAMplusC4D consortium and mouse GWAS of atherosclerosis from the Hybrid Mouse Diversity Panel (HMDP) followed by integration with functional multiomics human (STARNET and GTEx) and mouse (HMDP) databases. We found that mouse and human shared >75% of CAD causal pathways. Based on network topology, we then predicted key regulatory genes for both the shared pathways and species-specific pathways, which were further validated through the use of single cell data and the latest CAD GWAS. In sum, our results should serve as a much-needed guidance for which human CAD-causal pathways can or cannot be further evaluated for novel CAD therapies using mouse models.

From mouse to human

A deep analysis of multiple genomic datasets reveals which genetic pathways associated with atherosclerosis and coronary artery disease are shared between mice and humans.

ExplaiNAble BioLogical Age (ENABL Age): an artificial intelligence framework for interpretable biological age

BACKGROUND

Biological age is a measure of health that offers insights into ageing. The existing age clocks, although valuable, often trade off accuracy and interpretability. We introduce ExplaiNAble BioLogical Age (ENABL Age), a computational framework that combines machine-learning models with explainable artificial intelligence (XAI) methods to accurately estimate biological age with individualised explanations.

METHODS

To construct the ENABL Age clock, we first predicted an age-related outcome (eg, all-cause or cause-specific mortality), and then rescaled these predictions to estimate biological age, using UK Biobank and National Health and Nutrition Examination Survey (NHANES) datasets. We adapted existing XAI methods to decompose individual ENABL Ages into contributing risk factors. For broad accessibility, we developed two versions: ENABL Age-L, based on blood tests, and ENABL Age-Q, based on questionnaire characteristics. Finally, we validated diverse ageing mechanisms captured by each ENABL Age clock through genome-wide association studies (GWAS) association analyses.

FINDINGS

Our ENABL Age clock was significantly correlated with chronological age (r=0·7867, p<0·0001 for UK Biobank; r=0·7126, p<0·0001 for NHANES). These clocks distinguish individuals who are healthy (ie, their ENABL Age is lower than their chronological age) from those who are unhealthy (ie, their ENABL Age is higher than their chronological age), predicting mortality more effectively than existing clocks. Groups of individuals who were unhealthy showed approximately three to 12 times higher log hazard ratio than healthy groups, as per ENABL Age. The clocks achieved high mortality prediction power with an area under the receiver operating characteristic curve of 0·8179 for 5-year mortality and 0·8115 for 10-year mortality on the UK Biobank dataset, and 0·8935 for 5-year mortality and 0·9107 for 10-year mortality on the NHANES dataset. The individualised explanations that revealed the contribution of specific characteristics to ENABL Age provided insights into the important characteristics for ageing. An association analysis with risk factors and ageing-related morbidities and GWAS results on ENABL Age clocks trained on different mortality causes showed that each clock captures distinct ageing mechanisms.

INTERPRETATION

ENABL Age brings an important leap forward in the application of XAI for interpreting biological age clocks. ENABL Age also carries substantial potential in practical settings, assisting medical professionals in untangling the complexity of ageing mechanisms, and potentially becoming a valuable tool in informed clinical decision-making processes.

FUNDING

National Science Foundation and National Institutes of Health.

Type 2 Diabetes Modifies the Association of CAD Genomic Risk Variants With Subclinical Atherosclerosis

BACKGROUND

Individuals with type 2 diabetes (T2D) have an increased risk of coronary artery disease (CAD), but questions remain about the underlying pathology. Identifying which CAD loci are modified by T2D in the development of subclinical atherosclerosis (coronary artery calcification [CAC], carotid intima-media thickness, or carotid plaque) may improve our understanding of the mechanisms leading to the increased CAD in T2D.

METHODS

We compared the common and rare variant associations of known CAD loci from the literature on CAC, carotid intima-media thickness, and carotid plaque in up to 29 670 participants, including up to 24 157 normoglycemic controls and 5513 T2D cases leveraging whole-genome sequencing data from the Trans-Omics for Precision Medicine program. We included first-order T2D interaction terms in each model to determine whether CAD loci were modified by T2D. The genetic main and interaction effects were assessed using a joint test to determine whether a CAD variant, or gene-based rare variant set, was associated with the respective subclinical atherosclerosis measures and then further determined whether these loci had a significant interaction test.

RESULTS

Using a Bonferroni-corrected significance threshold of P<1.6×10-4, we identified 3 genes (ATP1B1, ARVCF, and LIPG) associated with CAC and 2 genes (ABCG8 and EIF2B2) associated with carotid intima-media thickness and carotid plaque, respectively, through gene-based rare variant set analysis. Both ATP1B1 and ARVCF also had significantly different associations for CAC in T2D cases versus controls. No significant interaction tests were identified through the candidate single-variant analysis.

CONCLUSIONS

These results highlight T2D as an important modifier of rare variant associations in CAD loci with CAC.