The Applications of Single-Cell RNA Sequencing in Atherosclerotic Disease

Atherosclerosis still is the primary cause of death worldwide. Our characterization of the atherosclerotic lesion is mainly rooted in definitions based on pathological descriptions. We often speak in absolutes regarding plaque phenotypes: vulnerable vs. stable plaques or plaque rupture vs. plaque erosion. By focusing on these concepts, we may have oversimplified the atherosclerotic disease and its mechanisms. The widely used definitions of pathology-based plaque phenotypes can be fine-tuned with observations made with various -omics techniques. Recent advancements in single-cell transcriptomics provide the opportunity to characterize the cellular composition of the atherosclerotic plaque. This additional layer of information facilitates the in-depth characterization of the atherosclerotic plaque. In this review, we discuss the impact that single-cell transcriptomics may exert on our current understanding of atherosclerosis

B-cell precursor acute lymphoblastic leukemia elicits an interferon-α/b response in bone marrow-derived mesenchymal stroma

B-cell precursor acute lymphoblastic leukemia (BCP-ALL) can hijack the normal bone marrow microenvironment to create a leukemic niche which facilitates blast cell survival and promotes drug resistance. Bone marrow-derived mesenchymal stromal cells (MSC) mimic this protective environment in ex vivo co-cultures with leukemic cells obtained from children with newly diagnosed BCP-ALL. We examined the potential mechanisms of this protection by RNA sequencing of flow-sorted MSC after co-culture with BCP-ALL cells. Leukemic cells induced an interferon (IFN)-related gene signature in MSC, which was partially dependent on direct cell-cell signaling. The signature was selectively induced by BCP-ALL cells, most profoundly by ETV6-RUNX1-positive ALL cells, as co-culture of MSC with healthy immune cells did not provoke a similar IFN signature. Leukemic cells and MSC both secreted IFNα and IFNb, but not IFNγ. In line, the IFN gene signature was sensitive to blockade of IFNα/b signaling, but less to that of IFNγ. The viability of leukemic cells and level of resistance to three chemotherapeutic agents was not affected by interference with IFN signaling using selective IFNα/b inhibitors or silencing of IFN-related genes. Taken together, our data suggest that the leukemia-induced expression of IFNα/b-related genes by MSC does not support survival of BCP-ALL cells but may serve a different role in the pathobiology of BCP-ALL

Human Plaque Myofibroblasts to Study Mechanisms of Atherosclerosis

Background Plaque myofibroblasts are critical players in the initiation and advancement of atherosclerotic disease. They are involved in the production of extracellular matrix, the formation of the fibrous cap, and the underlying lipidic core via modulation processes in response to different environmental cues. Despite clear phenotypic differences between myofibroblast cells and healthy vascular smooth muscle cells, smooth muscle cells are still widely used as a cellular model in atherosclerotic research. Methods and Results Here, we present a conditioned outgrowth method to isolate and culture myofibroblast cells from plaques. We obtained these cells from 27 donors (24 carotid and 3 femoral endarterectomies). We show that they keep their proliferative capacity for 8 passages, are transcriptionally stable, retain donor-specific gene expression programs, and express extracellular matrix proteins (FN1, COL1A1, and DCN) and smooth muscle cell markers (ACTA2, MYH11, and CNN1). Single-cell transcriptomics reveals that the cells in culture closely resemble the plaque myofibroblasts. Chromatin immunoprecipitation sequencing shows the presence of histone H3 lysine 4 dimethylation at the MYH11 promoter, pointing to their smooth muscle cell origin. Finally, we demonstrated that plaque myofibroblasts can be efficiently transduced (>97%) and are capable of taking up oxidized low-density lipoprotein and undergoing calcification. Conclusions In conclusion, we present a method to isolate and culture cells that retain plaque myofibroblast phenotypical and functional capabilities, making them a suitable in vitro model for studying selected mechanisms of atherosclerosis

Transcriptomic-based clustering of human atherosclerotic plaques identifies subgroups with different underlying biology and clinical presentation

Histopathological studies have revealed key processes of atherosclerotic plaque thrombosis. However, the diversity and complexity of lesion types highlight the need for improved sub-phenotyping. Here we analyze the gene expression profiles of 654 advanced human carotid plaques. The unsupervised, transcriptome-driven clustering revealed five dominant plaque types. These plaque phenotypes were associated with clinical presentation and showed differences in cellular compositions. Validation in coronary segments showed that the molecular signature of these plaques was linked to coronary ischemia. One of the plaque types with the most severe clinical symptoms pointed to both inflammatory and fibrotic cell lineages. Further, we did a preliminary analysis of potential circulating biomarkers that mark the different plaques phenotypes. In conclusion, the definition of the plaque at risk for a thrombotic event can be fine-tuned by in-depth transcriptomic-based phenotyping. These differential plaque phenotypes prove clinically relevant for both carotid and coronary artery plaques and point to distinct underlying biology of symptomatic lesions

Microanatomy of the Human Atherosclerotic Plaque by Single-Cell Transcriptomics

Rationale:Atherosclerotic lesions are known for their cellular heterogeneity, yet the molecular complexity within the cells of human plaques has not been fully assessed.Objective:Using single-cell transcriptomics and chromatin accessibility, we gained a better understanding of the pathophysiology underlying human atherosclerosis.Methods and Results:We performed single-cell RNA and single-cell ATAC sequencing on human carotid atherosclerotic plaques to define the cells at play and determine their transcriptomic and epigenomic characteristics. We identified 14 distinct cell populations including endothelial cells, smooth muscle cells, mast cells, B cells, myeloid cells, and T cells and identified multiple cellular activation states and suggested cellular interconversions. Within the endothelial cell population, we defined subsets with angiogenic capacity plus clear signs of endothelial to mesenchymal transition. CD4(+) and CD8(+) T cells showed activation-based subclasses, each with a gradual decline from a cytotoxic to a more quiescent phenotype. Myeloid cells included 2 populations of proinflammatory macrophages showing IL (interleukin) 1B or TNF (tumor necrosis factor) expression as well as a foam cell-like population expressing TREM2 (triggering receptor expressed on myeloid cells 2) and displaying a fibrosis-promoting phenotype. ATACseq data identified specific transcription factors associated with the myeloid subpopulation and T cell cytokine profiles underlying mutual activation between both cell types. Finally, cardiovascular disease susceptibility genes identified using public genome-wide association studies data were particularly enriched in lesional macrophages, endothelial, and smooth muscle cells.Conclusions:This study provides a transcriptome-based cellular landscape of human atherosclerotic plaques and highlights cellular plasticity and intercellular communication at the site of disease. This detailed definition of cell communities at play in atherosclerosis will facilitate cell-based mapping of novel interventional targets with direct functional relevance for the treatment of human disease

Atherosclerotic Plaque Epigenetic Age Acceleration Predicts a Poor Prognosis and Is Associated With Endothelial-to-Mesenchymal Transition in Humans

BACKGROUND: Epigenetic age estimators (clocks) are predictive of human mortality risk. However, it is not yet known whether the epigenetic age of atherosclerotic plaques is predictive for the risk of cardiovascular events. METHODS: Whole-genome DNA methylation of human carotid atherosclerotic plaques (n=485) and of blood (n=93) from the Athero-Express endarterectomy cohort was used to calculate epigenetic age acceleration (EAA). EAA was linked to clinical characteristics, plaque histology, and future cardiovascular events (n=136). We studied whole-genome DNA methylation and bulk and single-cell transcriptomics to uncover molecular mechanisms of plaque EAA. We experimentally confirmed our in silico findings using in vitro experiments in primary human coronary endothelial cells. RESULTS: Male and female patients with severe atherosclerosis had a median chronological age of 69 years. The median epigenetic age was 65 years in females (median EAA, -2.2 [interquartile range, -4.3 to 2.2] years) and 68 years in males (median EAA, -0.3 [interquartile range, -2.9 to 3.8] years). Patients with diabetes and a high body mass index had higher plaque EAA. Increased EAA of plaque predicted future events in a 3-year follow-up in a Cox regression model (univariate hazard ratio, 1.7; P=0.0034) and adjusted multivariate model (hazard ratio, 1.56; P=0.02). Plaque EAA predicted outcome independent of blood EAA (hazard ratio, 1.3; P=0.018) and of plaque hemorrhage (hazard ratio, 1.7; P=0.02). Single-cell RNA sequencing in plaque samples from 46 patients in the same cohort revealed smooth muscle and endothelial cells as important cell types in plaque EAA. Endothelial-to-mesenchymal transition was associated with EAA, which was experimentally confirmed by TGFβ-triggered endothelial-to-mesenchymal transition inducing rapid epigenetic aging in coronary endothelial cells. CONCLUSIONS: Plaque EAA is a strong and independent marker of poor outcome in patients with severe atherosclerosis. Plaque EAA was linked to mesenchymal endothelial and smooth muscle cells. Endothelial-to-mesenchymal transition was associated with EAA, which was experimentally validated. Epigenetic aging mechanisms may provide new targets for treatments that reduce atherosclerosis complications

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

Female Gene Networks Are Expressed in Myofibroblast-Like Smooth Muscle Cells in Vulnerable Atherosclerotic Plaques

BACKGROUND: Women presenting with coronary artery disease more often present with fibrous atherosclerotic plaques, which are currently understudied. Phenotypically modulated smooth muscle cells (SMCs) contribute to atherosclerosis in women. How these phenotypically modulated SMCs shape female versus male plaques is unknown. METHODS: Gene regulatory networks were created using RNAseq gene expression data from human carotid atherosclerotic plaques. The networks were prioritized based on sex bias, relevance for smooth muscle biology, and coronary artery disease genetic enrichment. Network expression was linked to histologically determined plaque phenotypes. In addition, their expression in plaque cell types was studied at single-cell resolution using single-cell RNAseq. Finally, their relevance for disease progression was studied in female and male Apoe -/- mice fed a Western diet for 18 and 30 weeks. RESULTS: Here, we identify multiple sex-stratified gene regulatory networks from human carotid atherosclerotic plaques. Prioritization of the female networks identified 2 main SMC gene regulatory networks in late-stage atherosclerosis. Single-cell RNA sequencing mapped these female networks to 2 SMC phenotypes: a phenotypically modulated myofibroblast-like SMC network and a contractile SMC network. The myofibroblast-like network was mostly expressed in plaques that were vulnerable in women. Finally, the mice ortholog of key driver gene MFGE8 (milk fat globule EGF and factor V/VIII domain containing) showed retained expression in advanced plaques from female mice but was downregulated in male mice during atherosclerosis progression. CONCLUSIONS: Female atherosclerosis is characterized by gene regulatory networks that are active in fibrous vulnerable plaques rich in myofibroblast-like SMCs

Female gene networks are expressed in myofibroblast-like smooth muscle cells in vulnerable atherosclerotic plaques

Women presenting with coronary artery disease (CAD) more often present with fibrous atherosclerotic plaques, which are currently understudied. Phenotypically modulated smooth muscle cells (SMCs) contribute to atherosclerosis in women. How these phenotypically modulated SMCs shape female versus male plaques is unknown. Here, we show sex-stratified gene regulatory networks (GRNs) from human carotid atherosclerotic tissue. Prioritization of these networks identified two main SMC GRNs in late-stage atherosclerosis. Single-cell RNA-sequencing mapped these GRNs to two SMC phenotypes: a phenotypically modulated myofibroblast-like SMC network and a contractile SMC network. The myofibroblast-like GRN was mostly expressed in plaques that were vulnerable in females. Finally, mice orthologs of the female myofibroblast-like genes showed retained expression in advanced plaques from female mice but were downregulated in male mice during atherosclerosis progression. Female atherosclerosis is driven by GRNs that promote a fibrous vulnerable plaque rich in myofibroblast-like SMCs

The applications of single-cell RNA sequencing in atherosclerotic research

At the base of the work described here is the overview of all different cell populations present in the carotid atherosclerotic plaque, which were previously not characterized in such detail. We unraveled cell populations such as smooth muscle cells (SMCs), endothelial cells (ECs) and various immune cells via single cell RNA-sequencing (scRNA-seq) which were characterized by examining their specific marker genes. We examined potential ligand-receptor interactions to study cell-cell communications within the plaque and found that myeloid, EC and SMCs are most likely to be interacting. In addition, we compared the human myeloid populations with mouse populations and reported a decent concordance in cell type diversity between species. We leveraged this new knowledge of the cell populations to study the wealth of genome wide association studies (GWAS) loci associated with cardiovascular disease. Identifying the genes from GWAS that are potential candidate genes for clinical care is not straightforward. We describe a two-step approach how we can use scRNA-seq data to translate GWAS susceptibility loci into targets for clinical care by directly examining expression in disease relevant tissue. Loci associated with coronary artery disease are dominantly associated with plaque SMC and EC populations. After identifying (novel) cell- gene pairs, expression of SKI, KANK2 and EDNRA was correlated with migration, proliferation and calcification in vascular SMCs. The molecular mechanisms of endothelial to mesenchymal transition (EndoMT) in human atherosclerosis are poorly understood. We combined the power of in vitro experimental data and scRNA-seq of a lineage tracing mouse model to help study EndoMT in human lesions. We assessed the temporal gene expression under EndoMT promoting conditions and defined assorted patterns of gene expression over the course of time. Next we studied scRNA-seq data from a lineage tracing mouse model were we composed endoMT trajectories. This data was used to demonstrate that the temporal gene expression changes from the in vitro model could identify EndoMT trajectories. Finally, we constructed three candidate EndoMT lineages across multiple subpopulations of ECs and SMCs in human carotid scRNA-seq data, similar to the approach that was applied in mouse data. We examined gene expression over the course of these lineages and identified 73 potential markers for mid-stage EndoMT. From these genes, NRG1 and DEPP1 are associated with mid-stage EndoMT in human atherosclerotic cells. Lastly we studied sex differences in human atherosclerosis via (sc)RNA-seq data. We compared SMC and EC sex-stratified gene expression profiles from 20 female and 26 male atherosclerotic lesions. And used deconvolution to estimate cell population proportions of RNA seq from 169 female and 485 male samples. Differences in sex-related biological pathways in sex stratified scRNA-seq data revealed that the differences in genes found between sexes do not have a common denominator such as “apoptosis” or “endoMT”, but could individually point to different processes. Deconvolution results indicated that female plaques showed a trend towards a higher SMC content compared to males and that males showed a trend towards a higher macrophage and immune cell content. Consequently, females leaned towards a higher ratio of structural vs. immune cells