Genome-wide association studies (GWAS) have been extensively used to study common complex diseases such as coronary artery disease (CAD), revealing 153 suggestive CAD loci, of which at least 46 have been validated as having genome-wide significance. However, these loci collectively explain <10% of the genetic variance in CAD. Thus, we must address the key question of what factors constitute the remaining 90% of CAD heritability. We review possible limitations of GWAS, and contextually consider some candidate CAD loci identified by this method. Looking ahead, we propose systems genetics as a complementary approach to unlocking the CAD heritability and etiology. Systems genetics builds network models of relevant molecular processes by combining genetic and genomic datasets to ultimately identify key "drivers" of disease. By leveraging systems-based genetic approaches, we can help reveal the full genetic basis of common complex disorders, enabling novel diagnostic and therapeutic opportunities.

Adeno-associated virus (AAV)-based vectors are promising vehicles for therapeutic gene delivery, including for the treatment for heart failure. It has been demonstrated for each of the AAV serotypes 1 through 8 that inhibition of the proteasome results in increased transduction efficiencies. For AAV9, however, the effect of proteasome inhibitors on in vivo transduction has until now not been evaluated. Here we demonstrate, in a well-established rodent heart failure model, that concurrent treatment with the proteasome inhibitor bortezomib does not enhance the efficacy of AAV9.SERCA2a to improve cardiac function as examined by echocardiography and pressure volume analysis. Western blot analysis of SERCA2a protein and reverse transcription-PCR of SERCA2a mRNA demonstrated that bortezomib had no effect on either endogenous rat SERCA2a levels nor on expression levels of human SERCA2a delivered by AAV9.SERCA2a. Similarly, the number of AAV9 genomes in heart samples was unaffected by bortezomib treatment. Interestingly, whereas transduction of HeLa cells and neonatal rat cardiomyocytes by AAV9 was stimulated by bortezomib, transduction of adult rat cardiomyocytes was inhibited. These results indicate an organ/cell-type-specific effect of proteasome inhibition on AAV9 transduction. A future detailed analysis of the underlying molecular mechanisms promises to facilitate the development of improved AAV vectors.

Endothelial to mesenchymal transition (EndMT) plays a major role during development, and also contributes to several adult cardiovascular diseases. Importantly, mesenchymal cells including fibroblasts are prominent in atherosclerosis, with key functions including regulation of: inflammation, matrix and collagen production, and plaque structural integrity. However, little is known about the origins of atherosclerosis-associated fibroblasts. Here we show using endothelial-specific lineage-tracking that EndMT-derived fibroblast-like cells are common in atherosclerotic lesions, with EndMT-derived cells expressing a range of fibroblast-specific markers. In vitro modelling confirms that EndMT is driven by TGF-β signalling, oxidative stress and hypoxia; all hallmarks of atherosclerosis. 'Transitioning' cells are readily detected in human plaques co-expressing endothelial and fibroblast/mesenchymal proteins, indicative of EndMT. The extent of EndMT correlates with an unstable plaque phenotype, which appears driven by altered collagen-MMP production in EndMT-derived cells. We conclude that EndMT contributes to atherosclerotic patho-biology and is associated with complex plaques that may be related to clinical events.

Atherosclerotic plaque rupture is the culprit event which underpins most acute vascular syndromes such as acute myocardial infarction. Novel biomarkers of plaque rupture could improve biological understanding and clinical management of patients presenting with possible acute vascular syndromes but such biomarker(s) remain elusive. Investigation of biomarkers in the context of de novo plaque rupture in humans is confounded by the inability to attribute the plaque rupture as the source of biomarker release, as plaque ruptures are typically associated with prompt down-stream events of myocardial necrosis and systemic inflammation.

PDGFRα+ mesenchymal progenitor cells are associated with pathological fibro-adipogenic processes. Conversely, a beneficial role for these cells during homeostasis or in response to revascularization and regeneration stimuli is suggested, but remains to be defined. We studied the molecular profile and function of PDGFRα+ cells in order to understand the mechanisms underlying their role in fibrosis versus regeneration. We show that PDGFRα+ cells are essential for tissue revascularization and restructuring through injury-stimulated remodeling of stromal and vascular components, context-dependent clonal expansion, and ultimate removal of pro-fibrotic PDGFRα+-derived cells. Tissue ischemia modulates the PDGFRα+ phenotype toward cells capable of remodeling the extracellular matrix and inducing cell-cell and cell-matrix adhesion, likely favoring tissue repair. Conversely, pathological healing occurs if PDGFRα+-derived cells persist as terminally differentiated mesenchymal cells. These studies support a context-dependent "yin-yang" biology of tissue-resident mesenchymal progenitor cells, which possess an innate ability to limit injury expansion while also promoting fibrosis in an unfavorable environment.

Reprogramming non-cardiomyocytes (non-CMs) into cardiomyocyte (CM)-like cells is a promising strategy for cardiac regeneration in conditions such as ischemic heart disease. Here, we used a modified mRNA (modRNA) gene delivery platform to deliver a cocktail, termed 7G-modRNA, of four cardiac-reprogramming genes-Gata4 (G), Mef2c (M), Tbx5 (T), and Hand2 (H)-together with three reprogramming-helper genes-dominant-negative (DN)-TGFβ, DN-Wnt8a, and acid ceramidase (AC)-to induce CM-like cells. We showed that 7G-modRNA reprogrammed 57% of CM-like cells in vitro. Through a lineage-tracing model, we determined that delivering the 7G-modRNA cocktail at the time of myocardial infarction reprogrammed ∼25% of CM-like cells in the scar area and significantly improved cardiac function, scar size, long-term survival, and capillary density. Mechanistically, we determined that while 7G-modRNA cannot create de novo beating CMs in vitro or in vivo, it can significantly upregulate pro-angiogenic mesenchymal stromal cells markers and transcription factors. We also demonstrated that our 7G-modRNA cocktail leads to neovascularization in ischemic-limb injury, indicating CM-like cells importance in other organs besides the heart. modRNA is currently being used around the globe for vaccination against COVID-19, and this study proves this is a safe, highly efficient gene delivery approach with therapeutic potential to treat ischemic diseases.

The majority of risk loci identified by genome-wide association studies (GWAS) are in non-coding regions, hampering their functional interpretation. Instead, transcriptome-wide association studies (TWAS) identify gene-trait associations, which can be used to prioritize candidate genes in disease-relevant tissue(s). Here, we aimed to systematically identify susceptibility genes for coronary artery disease (CAD) by TWAS. We trained prediction models of nine CAD-relevant tissues using EpiXcan based on two genetics-of-gene-expression panels, the Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task (STARNET) and the Genotype-Tissue Expression (GTEx). Based on these prediction models, we imputed gene expression of respective tissues from individual-level genotype data on 37,997 CAD cases and 42,854 controls for the subsequent gene-trait association analysis. Transcriptome-wide significant association (i.e. P < 3.85e-6) was observed for 114 genes. Of these, 96 resided within previously identified GWAS risk loci and 18 were novel. Stepwise analyses were performed to study their plausibility, biological function, and pathogenicity in CAD, including analyses for colocalization, damaging mutations, pathway enrichment, phenome-wide associations with human data and expression-traits correlations using mouse data. Finally, CRISPR/Cas9-based gene knockdown of two newly identified TWAS genes, RGS19 and KPTN, in a human hepatocyte cell line resulted in reduced secretion of APOB100 and lipids in the cell culture medium. Our CAD TWAS work (i) prioritized candidate causal genes at known GWAS loci, (ii) identified 18 novel genes to be associated with CAD, and iii) suggested potential tissues and pathways of action for these TWAS CAD genes.

Coronary artery disease (CAD) is characterized by atherosclerotic plaque formation in the arterial wall. CAD progression involves complex interactions and phenotypic plasticity among vascular and immune cell lineages. Single-cell RNA-seq (scRNA-seq) studies have highlighted lineage-specific transcriptomic signatures, but human cell phenotypes remain controversial. Here, we perform an integrated meta-analysis of 22 scRNA-seq libraries to generate a comprehensive map of human atherosclerosis with 118,578 cells. Besides characterizing granular cell-type diversity and communication, we leverage this atlas to provide insights into smooth muscle cell (SMC) modulation. We integrate genome-wide association study data and uncover a critical role for modulated SMC phenotypes in CAD, myocardial infarction, and coronary calcification. Finally, we identify fibromyocyte/fibrochondrogenic SMC markers (LTBP1 and CRTAC1) as proxies of atherosclerosis progression and validate these through omics and spatial imaging analyses. Altogether, we create a unified atlas of human atherosclerosis informing cell state-specific mechanistic and translational studies of cardiovascular diseases.

Non-atherosclerotic abnormalities of vessel calibre, aneurysm and ectasia, are challenging to quantify and are often overlooked in qualitative reporting. Utilising a novel 3-dimensional (3D) quantitative coronary angiography (QCA) application, we have evaluated the characteristics of normal, diabetic and aneurysmal or ectatic coronary arteries. We selected 131 individuals under 50 years-of-age, who had undergone coronary angiography for suspected myocardial ischaemia between 1st January 2011 and 31st December 2015, at the Bristol Heart Institute, Bristol, UK. This included 42 patients with angiographically normal coronary arteries, 36 diabetic patients with unobstructed coronaries, and 53 patients with abnormal coronary dilatation (aneurysm and ectasia). A total of 1105 coronary segments were analysed using QAngio XA 3D (Research Edition, Medis medical imaging systems, Leiden, The Netherlands). The combined volume of the major coronary arteries was significantly different between each group (1240 ± 476 mm3 diabetic group, 1646 ± 391 mm3 normal group, and 2072 ± 687 mm3 abnormal group). Moreover, the combined coronary artery volumes correlated with patient body surface area (r = 0.483, p < 0.01). Inter-observer variability was assessed and intraclass correlation coefficient of the total coronary artery volume demonstrated a low variability of 3D QCA (r = 0.996, p < 0.001). Dedicated 3D QCA facilitates reproducible coronary artery volume estimation and allows discrimination of normal and diseased vessels.

A library of well-characterized human induced pluripotent stem cell (hiPSC) lines from clinically healthy human subjects could serve as a useful resource of normal controls for in vitro human development, disease modeling, genotype-phenotype association studies, and drug response evaluation. We report generation and extensive characterization of a gender-balanced, racially/ethnically diverse library of hiPSC lines from 40 clinically healthy human individuals who range in age from 22 to 61 years. The hiPSCs match the karyotype and short tandem repeat identities of their parental fibroblasts, and have a transcription profile characteristic of pluripotent stem cells. We provide whole-genome sequencing data for one hiPSC clone from each individual, genomic ancestry determination, and analysis of mendelian disease genes and risks. We document similar transcriptomic profiles, single-cell RNA-sequencing-derived cell clusters, and physiology of cardiomyocytes differentiated from multiple independent hiPSC lines. This extensive characterization makes this hiPSC library a valuable resource for many studies on human biology.

Coronary artery disease (CAD) is a complex, multifactorial disease caused, in particular, by inflammation and cholesterol metabolism. At the molecular level, the role of tissue-specific signaling pathways leading to CAD is still largely unexplored. This study relied on two main resources: (1) genes with impact on atherosclerosis/CAD, and (2) liver-specific transcriptome analyses from human and mouse studies. The transcription factor activating transcription factor 3 (ATF3) was identified as a key regulator of a liver network relevant to atherosclerosis and linked to inflammation and cholesterol metabolism. ATF3 was predicted to be a direct and indirect (via MAF BZIP Transcription Factor F (MAFF)) regulator of low-density lipoprotein receptor (LDLR). Chromatin immunoprecipitation DNA sequencing (ChIP-seq) data from human liver cells revealed an ATF3 binding motif in the promoter regions of MAFF and LDLR. siRNA knockdown of ATF3 in human Hep3B liver cells significantly upregulated LDLR expression (p < 0.01). Inflammation induced by lipopolysaccharide (LPS) stimulation resulted in significant upregulation of ATF3 (p < 0.01) and subsequent downregulation of LDLR (p < 0.001). Liver-specific expression data from human CAD patients undergoing coronary artery bypass grafting (CABG) surgery (STARNET) and mouse models (HMDP) confirmed the regulatory role of ATF3 in the homeostasis of cholesterol metabolism. This study suggests that ATF3 might be a promising treatment candidate for lowering LDL cholesterol and reducing cardiovascular risk.

Elevated plasma cholesterol and type 2 diabetes (T2D) are associated with coronary artery disease (CAD). Individuals treated with cholesterol-lowering statins have increased T2D risk, while individuals with hypercholesterolemia have reduced T2D risk. We explore the relationship between lipid and glucose control by constructing network models from the STARNET study with sequencing data from seven cardiometabolic tissues obtained from CAD patients during coronary artery by-pass grafting surgery. By integrating gene expression, genotype, metabolomic, and clinical data, we identify a glucose and lipid determining (GLD) regulatory network showing inverse relationships with lipid and glucose traits. Master regulators of the GLD network also impact lipid and glucose levels in inverse directions. Experimental inhibition of one of the GLD network master regulators, lanosterol synthase (LSS), in mice confirms the inverse relationships to glucose and lipid levels as predicted by our model and provides mechanistic insights.

Definitive hematopoiesis emerges during embryogenesis via an endothelial-to-hematopoietic transition. We attempted to induce this process in mouse fibroblasts by screening a panel of factors for hemogenic activity. We identified a combination of four transcription factors, Gata2, Gfi1b, cFos, and Etv6, that efficiently induces endothelial-like precursor cells, with the subsequent appearance of hematopoietic cells. The precursor cells express a human CD34 reporter, Sca1, and Prominin1 within a global endothelial transcription program. Emergent hematopoietic cells possess nascent hematopoietic stem cell gene-expression profiles and cell-surface phenotypes. After transgene silencing and reaggregation culture, the specified cells generate hematopoietic colonies in vitro. Thus, we show that a simple combination of transcription factors is sufficient to induce a complex, dynamic, and multistep developmental program in vitro. These findings provide insights into the specification of definitive hemogenesis and a platform for future development of patient-specific stem and progenitor cells, as well as more-differentiated blood products.

Inferring molecular networks can reveal how genetic perturbations interact with environmental factors to cause common complex diseases. We analyzed genetic and gene expression data from seven tissues relevant to coronary artery disease (CAD) and identified regulatory gene networks (RGNs) and their key drivers. By integrating data from genome-wide association studies, we identified 30 CAD-causal RGNs interconnected in vascular and metabolic tissues, and we validated them with corresponding data from the Hybrid Mouse Diversity Panel. As proof of concept, by targeting the key drivers AIP, DRAP1, POLR2I, and PQBP1 in a cross-species-validated, arterial-wall RGN involving RNA-processing genes, we re-identified this RGN in THP-1 foam cells and independent data from CAD macrophages and carotid lesions. This characterization of the molecular landscape in CAD will help better define the regulation of CAD candidate genes identified by genome-wide association studies and is a first step toward achieving the goals of precision medicine.

Endothelial-mesenchymal transition (EndMT) is associated with various cardiovascular diseases and in particular with atherosclerosis and plaque instability. However, the molecular pathways that govern EndMT are poorly defined. Specifically, the role of epigenetic factors and histone deacetylases (HDACs) in controlling EndMT and the atherosclerotic plaque phenotype remains unclear. Here, we identified histone deacetylation, specifically that mediated by HDAC9 (a class IIa HDAC), as playing an important role in both EndMT and atherosclerosis. Using in vitro models, we found class IIa HDAC inhibition sustained the expression of endothelial proteins and mitigated the increase in mesenchymal proteins, effectively blocking EndMT. Similarly, ex vivo genetic knockout of Hdac9 in endothelial cells prevented EndMT and preserved a more endothelial-like phenotype. In vivo, atherosclerosis-prone mice with endothelial-specific Hdac9 knockout showed reduced EndMT and significantly reduced plaque area. Furthermore, these mice displayed a more favorable plaque phenotype, with reduced plaque lipid content and increased fibrous cap thickness. Together, these findings indicate that HDAC9 contributes to vascular pathology by promoting EndMT. Our study provides evidence for a pathological link among EndMT, HDAC9, and atherosclerosis and suggests that targeting of HDAC9 may be beneficial for plaque stabilization or slowing the progression of atherosclerotic disease.

Spontaneous Coronary Artery Dissection (SCAD) is an important cause of acute coronary syndromes, particularly in young to middle-aged women. Differentiating acute SCAD from coronary atherothrombosis remains a major clinical challenge.

To investigate percutaneous coronary intervention (PCI) practice in an international cohort of patients with spontaneous coronary artery dissection (SCAD). To explore factors associated with complications and study angiographic and longer term outcomes.

Direct evaluation of vascular inflammation in patients with COVID-19 would facilitate more efficient trials of new treatments and identify patients at risk of long-term complications who might respond to treatment. We aimed to develop a novel artificial intelligence (AI)-assisted image analysis platform that quantifies cytokine-driven vascular inflammation from routine CT angiograms, and sought to validate its prognostic value in COVID-19.

The associations between diseases/traits and copy number variants (CNVs) have not been systematically investigated in genome-wide association studies (GWASs), primarily due to a lack of robust and accurate tools for CNV genotyping. Herein, we propose a novel ensemble learning framework, ensembleCNV, to detect and genotype CNVs using single nucleotide polymorphism (SNP) array data. EnsembleCNV (a) identifies and eliminates batch effects at raw data level; (b) assembles individual CNV calls into CNV regions (CNVRs) from multiple existing callers with complementary strengths by a heuristic algorithm; (c) re-genotypes each CNVR with local likelihood model adjusted by global information across multiple CNVRs; (d) refines CNVR boundaries by local correlation structure in copy number intensities; (e) provides direct CNV genotyping accompanied with confidence score, directly accessible for downstream quality control and association analysis. Benchmarked on two large datasets, ensembleCNV outperformed competing methods and achieved a high call rate (93.3%) and reproducibility (98.6%), while concurrently achieving high sensitivity by capturing 85% of common CNVs documented in the 1000 Genomes Project. Given this CNV call rate and accuracy, which are comparable to SNP genotyping, we suggest ensembleCNV holds significant promise for performing genome-wide CNV association studies and investigating how CNVs predispose to human diseases.

Fibromuscular dysplasia (FMD) and spontaneous coronary artery dissection (SCAD) are related, non-atherosclerotic arterial diseases mainly affecting middle-aged women. Little is known about their physiopathological mechanisms. We aimed to identify rare genetic causes to elucidate molecular mechanisms implicated in FMD and SCAD.