The 9p21.3 coronary artery disease risk locus drives vascular smooth muscle cells to an osteochondrogenic state.

Genome-wide association studies have identified common genetic variants at ~400 human genomic loci linked to coronary artery disease (CAD) susceptibility. Among these genomic regions, the most impactful is the 9p21.3 CAD risk locus, which spans a 60 kb gene desert and encompasses ~80 SNPs in high linkage disequilibrium. Despite nearly two decades since its discovery, the functional mechanism of this genomic region remains incompletely resolved. To investigate the transcriptional gene programs mediated by 9p21.3 risk locus, we applied a model of induced pluripotent stem cells (iPSCs) from risk and non-risk donors at 9p21.3, as well as isogenic lines with a full haplotype deletion. Upon differentiation to vascular smooth muscle cells (VSMC), single-cell transcriptomic profiling demonstrated iPSC-VSMC phenotypes resembling those from native human coronary arteries, confirming the robustness of this model. Remarkably, our analyses revealed that VSMCs harboring the 9p21.3 risk haplotype preferentially adopt an osteochondrogenic state. Importantly, we identified LIMCH1 and CRABP1 as signature genes critical for defining this transcriptional program. Our study provides new insights into the mechanism at the 9p21.3 risk locus and defines its role in driving a disease-prone transcriptional state in VSMCs.

CD163+ Macrophages Induce Endothelial-to-Mesenchymal Transition in Atheroma.

BACKGROUND: Cell phenotype switching is increasingly being recognized in atherosclerosis. However, our understanding of the exact stimuli for such cellular transformations and their significance for human atherosclerosis is still evolving. Intraplaque hemorrhage is thought to be a major contributor to plaque progression in part by stimulating the influx of CD163+ macrophages. Here, we explored the hypothesis that CD163+ macrophages cause plaque progression through the induction of proapoptotic endothelial-to-mesenchymal transition (EndMT) within the fibrous cap. METHODS: Human coronary artery sections from CVPath's autopsy registry were selected for pathological analysis. Athero-prone ApoE-/- and ApoE-/-/CD163-/- mice were used for in vivo studies. Human peripheral blood mononuclear cell-induced macrophages and human aortic endothelial cells were used for in vitro experiments. RESULTS: In 107 lesions with acute coronary plaque rupture, 55% had pathological evidence of intraplaque hemorrhage in nonculprit vessels/lesions. Thinner fibrous cap, greater CD163+ macrophage accumulation, and a larger number of CD31/FSP-1 (fibroblast specific protein-1) double-positive cells and TUNEL (terminal deoxynucleotidyl transferase-dUTP nick end labeling) positive cells in the fibrous cap were observed in nonculprit intraplaque hemorrhage lesions, as well as in culprit rupture sections versus nonculprit fibroatheroma sections. Human aortic endothelial cells cultured with supernatants from hemoglobin/haptoglobin-exposed macrophages showed that increased mesenchymal marker proteins (transgelin and FSP-1) while endothelial markers (VE-cadherin and CD31) were reduced, suggesting EndMT induction. Activation of NF-κB (nuclear factor kappa ß) signaling by proinflammatory cytokines released from CD163+ macrophages directly regulated the expression of Snail, a critical transcription factor during EndMT induction. Western blot analysis for cleaved caspase-3 and microarray analysis of human aortic endothelial cells indicated that apoptosis was stimulated during CD163+ macrophage-induced EndMT. Additionally, CD163 deletion in athero-prone mice suggested that CD163 is required for EndMT and plaque progression. Using single-cell RNA sequencing from human carotid endarterectomy lesions, a population of EndMT was detected, which demonstrated significant upregulation of apoptosis-related genes. CONCLUSIONS: CD163+ macrophages provoke EndMT, which may promote plaque progression through fibrous cap thinning.

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.

Integrative single-cell meta-analysis reveals disease-relevant vascular cell states and markers in human atherosclerosis.

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.

FHL5 Controls Vascular Disease-Associated Gene Programs in Smooth Muscle Cells.

BACKGROUND: Genome-wide association studies have identified hundreds of loci associated with common vascular diseases, such as coronary artery disease, myocardial infarction, and hypertension. However, the lack of mechanistic insights for many GWAS loci limits their translation into the clinic. Among these loci with unknown functions is UFL1-four-and-a-half LIM (LIN-11, Isl-1, MEC-3) domain 5 (FHL5; chr6q16.1), which reached genome-wide significance in a recent coronary artery disease/ myocardial infarction GWAS meta-analysis. UFL1-FHL5 is also associated with several vascular diseases, consistent with the widespread pleiotropy observed for GWAS loci. METHODS: We apply a multimodal approach leveraging statistical fine-mapping, epigenomic profiling, and ex vivo analysis of human coronary artery tissues to implicate FHL5 as the top candidate causal gene. We unravel the molecular mechanisms of the cross-phenotype genetic associations through in vitro functional analyses and epigenomic profiling experiments in coronary artery smooth muscle cells. RESULTS: We prioritized FHL5 as the top candidate causal gene at the UFL1-FHL5 locus through expression quantitative trait locus colocalization methods. FHL5 gene expression was enriched in the smooth muscle cells and pericyte population in human artery tissues with coexpression network analyses supporting a functional role in regulating smooth muscle cell contraction. Unexpectedly, under procalcifying conditions, FHL5 overexpression promoted vascular calcification and dysregulated processes related to extracellular matrix organization and calcium handling. Lastly, by mapping FHL5 binding sites and inferring FHL5 target gene function using artery tissue gene regulatory network analyses, we highlight regulatory interactions between FHL5 and downstream coronary artery disease/myocardial infarction loci, such as FOXL1 and FN1 that have roles in vascular remodeling. CONCLUSIONS: Taken together, these studies provide mechanistic insights into the pleiotropic genetic associations of UFL1-FHL5. We show that FHL5 mediates vascular disease risk through transcriptional regulation of downstream vascular remodeling gene programs. These transacting mechanisms may explain a portion of the heritable risk for complex vascular diseases.

Gene-repressing epigenetic reader EED unexpectedly enhances cyclinD1 gene activation.

Epigenetically switched, proliferative vascular smooth muscle cells (SMCs) form neointima, engendering stenotic diseases. Histone-3 lysine-27 trimethylation (H3K27me3) and acetylation (H3K27ac) marks are associated with gene repression and activation, respectively. The polycomb protein embryonic ectoderm development (EED) reads H3K27me3 and also enhances its deposition, hence is a canonical gene repressor. However, herein we found an unexpected role for EED in activating the bona fide pro-proliferative gene Ccnd1 (cyclinD1). EED overexpression in SMCs increased Ccnd1 mRNA, seemingly contradicting its gene-repressing function. However, consistently, EED co-immunoprecipitated with gene-activating H3K27ac reader BRD4, and they co-occupied at both mitogen-activated Ccnd1 and mitogen-repressed P57 (bona fide anti-proliferative gene), as indicated by chromatin immunoprecipitation qPCR. These results were abolished by an inhibitor of either the EED/H3K27me3 or BRD4/H3K27ac reader function. In accordance, elevating BRD4 increased H3K27me3. In vivo, while EED was upregulated in rat and human neointimal lesions, selective EED inhibition abated angioplasty-induced neointima and reduced cyclinD1 in rat carotid arteries. Thus, results uncover a previously unknown role for EED in Ccnd1 activation, likely via its cooperativity with BRD4 that enhances each other's reader function; i.e., activating pro-proliferative Ccnd1 while repressing anti-proliferative P57. As such, this study confers mechanistic implications for the epigenetic intervention of neointimal pathology.

Integrative multi-ancestry genetic analysis of gene regulation in coronary arteries prioritizes disease risk loci.

Genome-wide association studies (GWAS) have identified hundreds of genetic risk loci for coronary artery disease (CAD). However, non-European populations are underrepresented in GWAS and the causal gene-regulatory mechanisms of these risk loci during atherosclerosis remain unclear. We incorporated local ancestry and haplotype information to identify quantitative trait loci (QTL) for gene expression and splicing in coronary arteries obtained from 138 ancestrally diverse Americans. Of 2,132 eQTL-associated genes (eGenes), 47% were previously unreported in coronary arteries and 19% exhibited cell-type-specific expression. Colocalization analysis with GWAS identified subgroups of eGenes unique to CAD and blood pressure. Fine-mapping highlighted additional eGenes of interest, including TBX20 and IL5 . Splicing (s)QTLs for 1,690 genes were also identified, among which TOR1AIP1 and ULK3 sQTLs demonstrated the importance of evaluating splicing events to accurately identify disease-relevant gene expression. Our work provides the first human coronary artery eQTL resource from a patient sample and exemplifies the necessity of diverse study populations and multi-omic approaches to characterize gene regulation in critical disease processes.

CD163+ macrophages restrain vascular calcification, promoting the development of high-risk plaque.

Vascular calcification (VC) is concomitant with atherosclerosis, yet it remains uncertain why rupture-prone high-risk plaques do not typically show extensive calcification. Intraplaque hemorrhage (IPH) deposits erythrocyte-derived cholesterol, enlarging the necrotic core and promoting high-risk plaque development. Pro-atherogenic CD163+ alternative macrophages engulf hemoglobin:haptoglobin (HH) complexes at IPH sites. However, their role in VC has never been examined to our knowledge. Here we show, in human arteries, the distribution of CD163+ macrophages correlated inversely with VC. In vitro experiments using vascular smooth muscle cells (VSMCs) cultured with HH-exposed human macrophage - M(Hb) - supernatant reduced calcification, while arteries from ApoE-/- CD163-/- mice showed greater VC. M(Hb) supernatant-exposed VSMCs showed activated NF-κB, while blocking NF-κB attenuated the anticalcific effect of M(Hb) on VSMCs. CD163+ macrophages altered VC through NF-κB-induced transcription of hyaluronan synthase (HAS), an enzyme that catalyzes the formation of the extracellular matrix glycosaminoglycan, hyaluronan, within VSMCs. M(Hb) supernatants enhanced HAS production in VSMCs, while knocking down HAS attenuated its anticalcific effect. NF-κB blockade in ApoE-/- mice reduced hyaluronan and increased VC. In human arteries, hyaluronan and HAS were increased in areas of CD163+ macrophage presence. Our findings highlight an important mechanism by which CD163+ macrophages inhibit VC through NF-κB-induced HAS augmentation and thus promote the high-risk plaque development.

PlaqView 2.0: A comprehensive web portal for cardiovascular single-cell genomics.

Single-cell RNA-seq (scRNA-seq) is a powerful genomics technology to interrogate the cellular composition and behaviors of complex systems. While the number of scRNA-seq datasets and available computational analysis tools have grown exponentially, there are limited systematic data sharing strategies to allow rapid exploration and re-analysis of single-cell datasets, particularly in the cardiovascular field. We previously introduced PlaqView, an open-source web portal for the exploration and analysis of published atherosclerosis single-cell datasets. Now, we introduce PlaqView 2.0 (www.plaqview.com), which provides expanded features and functionalities as well as additional cardiovascular single-cell datasets. We showcase improved PlaqView functionality, backend data processing, user-interface, and capacity. PlaqView brings new or improved tools to explore scRNA-seq data, including gene query, metadata browser, cell identity prediction, ad hoc RNA-trajectory analysis, and drug-gene interaction prediction. PlaqView serves as one of the largest central repositories for cardiovascular single-cell datasets, which now includes data from human aortic aneurysm, gene-specific mouse knockouts, and healthy references. PlaqView 2.0 brings advanced tools and high-performance computing directly to users without the need for any programming knowledge. Lastly, we outline steps to generalize and repurpose PlaqView's framework for single-cell datasets from other fields.

Single-nucleus chromatin accessibility profiling highlights regulatory mechanisms of coronary artery disease risk.

Coronary artery disease (CAD) is a complex inflammatory disease involving genetic influences across cell types. Genome-wide association studies have identified over 200 loci associated with CAD, where the majority of risk variants reside in noncoding DNA sequences impacting cis-regulatory elements. Here, we applied single-nucleus assay for transposase-accessible chromatin with sequencing to profile 28,316 nuclei across coronary artery segments from 41 patients with varying stages of CAD, which revealed 14 distinct cellular clusters. We mapped ~320,000 accessible sites across all cells, identified cell-type-specific elements and transcription factors, and prioritized functional CAD risk variants. We identified elements in smooth muscle cell transition states (for example, fibromyocytes) and functional variants predicted to alter smooth muscle cell- and macrophage-specific regulation of MRAS (3q22) and LIPA (10q23), respectively. We further nominated key driver transcription factors such as PRDM16 and TBX2. Together, this single-nucleus atlas provides a critical step towards interpreting regulatory mechanisms across the continuum of CAD risk.

GenomicDistributions: fast analysis of genomic intervals with Bioconductor.

BACKGROUND: Epigenome analysis relies on defined sets of genomic regions output by widely used assays such as ChIP-seq and ATAC-seq. Statistical analysis and visualization of genomic region sets is essential to answer biological questions in gene regulation. As the epigenomics community continues generating data, there will be an increasing need for software tools that can efficiently deal with more abundant and larger genomic region sets. Here, we introduce GenomicDistributions, an R package for fast and easy summarization and visualization of genomic region data. RESULTS: GenomicDistributions offers a broad selection of functions to calculate properties of genomic region sets, such as feature distances, genomic partition overlaps, and more. GenomicDistributions functions are meticulously optimized for best-in-class speed and generally outperform comparable functions in existing R packages. GenomicDistributions also offers plotting functions that produce editable ggplot objects. All GenomicDistributions functions follow a uniform naming scheme and can handle either single or multiple region set inputs. CONCLUSIONS: GenomicDistributions offers a fast and scalable tool for exploratory genomic region set analysis and visualization. GenomicDistributions excels in user-friendliness, flexibility of outputs, breadth of functions, and computational performance. GenomicDistributions is available from Bioconductor ( https://bioconductor.org/packages/release/bioc/html/GenomicDistributions.html ).

Enhanced single-cell RNA-seq workflow reveals coronary artery disease cellular cross-talk and candidate drug targets.

BACKGROUND AND AIMS: The atherosclerotic plaque microenvironment is highly complex, and selective agents that modulate plaque stability are not yet available. We sought to develop a scRNA-seq analysis workflow to investigate this environment and uncover potential therapeutic approaches. We designed a user-friendly, reproducible workflow that will be applicable to other disease-specific scRNA-seq datasets. METHODS: Here we incorporated automated cell labeling, pseudotemporal ordering, ligand-receptor evaluation, and drug-gene interaction analysis into a ready-to-deploy workflow. We applied this pipeline to further investigate a previously published human coronary single-cell dataset by Wirka et al. Notably, we developed an interactive web application to enable further exploration and analysis of this and other cardiovascular single-cell datasets. RESULTS: We revealed distinct derivations of fibroblast-like cells from smooth muscle cells (SMCs), and showed the key changes in gene expression along their de-differentiation path. We highlighted several key ligand-receptor interactions within the atherosclerotic environment through functional expression profiling and revealed several avenues for future pharmacological development for precision medicine. Further, our interactive web application, PlaqView (www.plaqview.com), allows lay scientists to explore this and other datasets and compare scRNA-seq tools without prior coding knowledge. CONCLUSIONS: This publicly available workflow and application will allow for more systematic and user-friendly analysis of scRNA datasets in other disease and developmental systems. Our analysis pipeline provides many hypothesis-generating tools to unravel the etiology of coronary artery disease. We also highlight potential mechanisms for several drugs in the atherosclerotic cellular environment. Future releases of PlaqView will feature more scRNA-seq and scATAC-seq atherosclerosis-related datasets to provide a critical resource for the field, and to promote data harmonization and biological interpretation.

Nuclear lipid droplets and nuclear damage in Caenorhabditis elegans.

Fat stored in the form of lipid droplets has long been considered a defining characteristic of cytoplasm. However, recent studies have shown that nuclear lipid droplets occur in multiple cells and tissues, including in human patients with fatty liver disease. The function(s) of stored fat in the nucleus has not been determined, and it is possible that nuclear fat is beneficial in some situations. Conversely, nuclear lipid droplets might instead be deleterious by disrupting nuclear organization or triggering aggregation of hydrophobic proteins. We show here that nuclear lipid droplets occur normally in C. elegans intestinal cells and germ cells, but appear to be associated with damage only in the intestine. Lipid droplets in intestinal nuclei can be associated with novel bundles of microfilaments (nuclear actin) and membrane tubules that might have roles in damage repair. To increase the normal, low frequency of nuclear lipid droplets in wild-type animals, we used a forward genetic screen to isolate mutants with abnormally large or abundant nuclear lipid droplets. Genetic analysis and cloning of three such mutants showed that the genes encode the lipid regulator SEIP-1/seipin, the inner nuclear membrane protein NEMP-1/Nemp1/TMEM194A, and a component of COPI vesicles called COPA-1/α-COP. We present several lines of evidence that the nuclear lipid droplet phenotype of copa-1 mutants results from a defect in retrieving mislocalized membrane proteins that normally reside in the endoplasmic reticulum. The seip-1 mutant causes most germ cells to have nuclear lipid droplets, the largest of which occupy more than a third of the nuclear volume. Nevertheless, the nuclear lipid droplets do not trigger apoptosis, and the germ cells differentiate into gametes that produce viable, healthy progeny. Thus, our results suggest that nuclear lipid droplets are detrimental to intestinal nuclei, but have no obvious deleterious effect on germ nuclei.