Intersecting single-cell transcriptomics and genome-wide association studies identifies crucial cell populations and candidate genes for atherosclerosis

High-Dimensional Single-Cell Multimodal Landscape of Human Carotid Atherosclerosis

BACKGROUND

Atherosclerotic plaques are complex tissues composed of a heterogeneous mixture of cells. However, our understanding of the comprehensive transcriptional and phenotypic landscape of the cells within these lesions is limited.

METHODS

To characterize the landscape of human carotid atherosclerosis in greater detail, we combined cellular indexing of transcriptomes and epitopes by sequencing and single-cell RNA sequencing to classify all cell types within lesions (n=21; 13 symptomatic) to achieve a comprehensive multimodal understanding of the cellular identities of atherosclerosis and their association with clinical pathophysiology.

RESULTS

We identified 25 cell populations, each with a unique multiomic signature, including macrophages, T cells, NK (natural killer) cells, mast cells, B cells, plasma cells, neutrophils, dendritic cells, endothelial cells, fibroblasts, and smooth muscle cells (SMCs). Among the macrophages, we identified 2 proinflammatory subsets enriched in IL-1B (interleukin-1B) or C1Q expression, 2 TREM2-positive foam cells (1 expressing inflammatory genes), and subpopulations with a proliferative gene signature and SMC-specific gene signature with fibrotic pathways upregulated. Further characterization revealed various subsets of SMCs and fibroblasts, including SMC-derived foam cells. These foamy SMCs were localized in the deep intima of coronary atherosclerotic lesions. Utilizing cellular indexing of transcriptomes and epitopes by sequencing data, we developed a flow cytometry panel, using cell surface proteins CD29, CD142, and CD90, to isolate SMC-derived cells from lesions. Lastly, we observed reduced proportions of efferocytotic macrophages, classically activated endothelial cells, and contractile and modulated SMC-derived cells, while inflammatory SMCs were enriched in plaques of clinically symptomatic versus asymptomatic patients.

CONCLUSIONS

Our multimodal atlas of cell populations within atherosclerosis provides novel insights into the diversity, phenotype, location, isolation, and clinical relevance of the unique cellular composition of human carotid atherosclerosis. These findings facilitate both the mapping of cardiovascular disease susceptibility loci to specific cell types and the identification of novel molecular and cellular therapeutic targets for the treatment of the disease.

Vascular smooth muscle cells in response to cholesterol crystals modulates inflammatory cytokines release and promotes neutrophil extracellular trap formation

BACKGROUND

The formation and accumulation of cholesterol crystals (CC) at the lesion site is a hallmark of atherosclerosis. Although studies have shown the importance of vascular smooth muscle cells (VSMCs) in the disease atherosclerosis, little is known about the molecular mechanism behind the uptake of CC in VSMCs and their role in modulating immune response.

METHODS

Human aortic smooth muscle cells were cultured and treated with CC. CC uptake and CC mediated signaling pathway and protein induction were studied using flow cytometry, confocal microscopy, western blot and Olink proteomics. Conditioned medium from CC treated VSMCs was used to study neutrophil adhesion, ROS production and phagocytosis. Neutrophil extracellular traps (NETs) formations were visualized using confocal microscopy.

RESULTS

VSMCs and macrophages were found around CC clefts in human carotid plaques. CC uptake in VSMCs are largely through micropinocytosis and phagocytosis via PI3K-AkT dependent pathway. The uptake of CC in VSMCs induce the release inflammatory proteins, including IL-33, an alarming cytokine. Conditioned medium from CC treated VSMCs can induce neutrophil adhesion, neutrophil reactive oxygen species (ROS) and neutrophil extracellular traps (NETs) formation. IL-33 neutralization in conditioned medium from CC treated VSMCs inhibited neutrophil ROS production and NETs formation.

CONCLUSION

We demonstrate that VSMCs due to its vicinity to CC clefts in human atherosclerotic lesion can modulate local immune response and we further reveal that the interaction between CC and VSMCs impart an inflammatory milieu in the atherosclerotic microenvironment by promoting IL-33 dependent neutrophil influx and NETs formation.

Translatome profiling reveals Itih4 as a novel smooth muscle cell-specific gene in atherosclerosis

AIMS

Vascular smooth muscle cells (SMCs) and their derivatives are key contributors to the development of atherosclerosis. However, studying changes in SMC gene expression in heterogeneous vascular tissues is challenging due to the technical limitations and high cost associated with current approaches. In this paper, we apply Translating Ribosome Affinity Purification sequencing (TRAP-Seq) to profile SMC-specific gene expression directly from tissue.

METHODS AND RESULTS

To facilitate SMC-specific translatome analysis, we generated SMCTRAP mice, a transgenic mouse line expressing EGFP-tagged ribosomal protein L10a (EGFP-L10a) under the control of the SMC-specific αSMA promoter. These mice were further crossed with the atherosclerosis model Ldlr-/-, ApoB100/100 to generate SMCTRAP-AS mice and used to profile atherosclerosis-associated SMCs in thoracic aorta samples of 15-month-old SMCTRAP and SMCTRAP-AS mice. Our analysis of SMCTRAP-AS mice showed that EGFP-L10a expression was localized to SMCs in various tissues, including the aortic wall and plaque. The TRAP fraction demonstrated high enrichment of known SMC-specific genes, confirming the specificity of our approach. We identified several genes, including Cemip, Lum, Mfge8, Spp1, and Serpina3, that are known to be involved in atherosclerosis-induced gene expression. Moreover, we identified several novel genes not previously linked to SMCs in atherosclerosis, such as Anxa4, Cd276, Itih4, Myof, Pcdh11x, Rab31, Serpinb6b, Slc35e4, Slc8a3, and Spink5. Among them, we confirmed the SMC-specific expression of Itih4 in atherosclerotic lesions using immunofluorescence staining of mouse aortic roots and spatial transcriptomics of human carotid arteries. Furthermore, our more detailed analysis of Itih4 showed its link to coronary artery disease (CAD) through the colocalization of GWAS, splice-QTL, and protein-QTL signals.

CONCLUSIONS

We generated a SMC-specific TRAP mouse line to study atherosclerosis and identified Itih4 as a novel SMC-expressed gene in atherosclerotic plaques, warranting further investigation of its putative function in extracellular matrix stability and genetic evidence of causality.

Emerging applications of single-cell profiling in precision medicine of atherosclerosis

Atherosclerosis is a chronic, progressive, inflammatory disease that occurs in the arterial wall. Despite recent advancements in treatment aimed at improving efficacy and prolonging survival, atherosclerosis remains largely incurable. In this review, we discuss emerging single-cell sequencing techniques and their novel insights into atherosclerosis. We provide examples of single-cell profiling studies that reveal phenotypic characteristics of atherosclerosis plaques, blood, liver, and the intestinal tract. Additionally, we highlight the potential clinical applications of single-cell analysis and propose that combining this approach with other techniques can facilitate early diagnosis and treatment, leading to more accurate medical interventions.

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.

Clonal Proliferation Within Smooth Muscle Cells in Unstable Human Atherosclerotic Lesions

BACKGROUND

Studies in humans and mice using the expression of an X-linked gene or lineage tracing, respectively, have suggested that clones of smooth muscle cells (SMCs) exist in human atherosclerotic lesions but are limited by either spatial resolution or translatability of the model.

METHODS

Phenotypic clonality can be detected by X-chromosome inactivation patterns. We investigated whether clones of SMCs exist in unstable human atheroma using RNA in situ hybridization (BaseScope) to identify a naturally occurring 24-nucleotide deletion in the 3'UTR of the X-linked BGN (biglycan) gene, a proteoglycan highly expressed by SMCs. BGN-specific BaseScope probes were designed to target the wild-type or deletion mRNA. Three different coronary artery plaque types (erosion, rupture, and adaptive intimal thickening) were selected from heterozygous females for the deletion BGN. Hybridization of target RNA-specific probes was used to visualize the spatial distribution of mutants. A clonality index was calculated from the percentage of each probe in each region of interest. Spatial transcriptomics were used to identify differentially expressed transcripts within clonal and nonclonal regions.

RESULTS

Less than one-half of regions of interest in the intimal plaque were considered clonal with the mean percent regions of interest with clonality higher in the intimal plaque than in the media. This was consistent for all plaque types. The relationship of the dominant clone in the intimal plaque and media showed significant concordance. In comparison with the nonclonal lesions, the regions with SMC clonality had lower expression of genes encoding cell growth suppressors such as CD74, SERF-2 (small EDRK-rich factor 2), CTSB (cathepsin B), and HLA-DPA1 (major histocompatibility complex, class II, DP alpha 1), among others.

CONCLUSIONS

Our novel approach to examine clonality suggests atherosclerosis is primarily a disease of polyclonally and to a lesser extent clonally expanded SMCs and may have implications for the development of antiatherosclerotic therapies.

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.

The Greatly Under-Represented Role of Smooth Muscle Cells in Atherosclerosis

PURPOSE OF REVIEW

This article summarizes previous and recent research on the fundamental role of arterial smooth muscle cells (SMCs) as drivers of initial and, along with macrophages, later stages of human atherosclerosis.

RECENT FINDINGS

Studies using human tissues and SMC lineage-tracing mice have reinforced earlier observations that SMCs drive initial atherogenesis in humans and contribute a multitude of phenotypes including foam cell formation hitherto attributed primarily to macrophages in atherosclerosis. Arterial smooth muscle cells (SMCs) are the primary cell type in human pre-atherosclerotic intima and are responsible for the retention of lipoproteins that drive the development of atherosclerosis. Despite this, images of atherogenesis still depict the process as initially devoid of SMCs, primarily macrophage driven, and indicate only relatively minor roles such as fibrous cap formation to intimal SMCs. This review summarizes historical and recent observations regarding the importance of SMCs in the formation of a pre-atherosclerotic intima, initial and later foam cell formation, and the phenotypic changes that give rise to multiple different roles for SMCs in human and mouse lesions. Potential SMC-specific therapies in atherosclerosis are presented.

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.

Multi-ancestry genome-wide study identifies effector genes and druggable pathways for coronary artery calcification

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

Tissue Inhibitor of Metalloproteinases-1 Interacts with CD74 to Promote AKT Signaling, Monocyte Recruitment Responses, and Vascular Smooth Muscle Cell Proliferation

Tissue inhibitor of metalloproteinases-1 (TIMP-1), an important regulator of matrix metalloproteinases (MMPs), has recently been shown to interact with CD74, a receptor for macrophage migration inhibitory factor (MIF). However, the biological effects mediated by TIMP-1 through CD74 remain largely unexplored. Using sequence alignment and in silico protein-protein docking analysis, we demonstrated that TIMP-1 shares residues with both MIF and MIF-2, crucial for CD74 binding, but not for CXCR4. Subcellular colocalization, immunoprecipitation, and internalization experiments supported these findings, demonstrating that TIMP-1 interacts with surface-expressed CD74, resulting in its internalization in a dose-dependent manner, as well as with a soluble CD74 ectodomain fragment (sCD74). This prompted us to study the effects of the TIMP-1-CD74 axis on monocytes and vascular smooth muscle cells (VSCMs) to assess its impact on vascular inflammation. A phospho-kinase array revealed the activation of serine/threonine kinases by TIMP-1 in THP-1 pre-monocytes, in particular AKT. Similarly, TIMP-1 dose-dependently triggered the phosphorylation of AKT and ERK1/2 in primary human monocytes. Importantly, Transwell migration, 3D-based Chemotaxis, and flow adhesion assays demonstrated that TIMP-1 engagement of CD74 strongly promotes the recruitment response of primary human monocytes, while live cell imaging studies revealed a profound activating effect on VSMC proliferation. Finally, re-analysis of scRNA-seq data highlighted the expression patterns of TIMP-1 and CD74 in human atherosclerotic lesions, thus, together with our experimental data, indicating a role for the TIMP-1-CD74 axis in vascular inflammation and atherosclerosis.

High-Dimensional Single-Cell Multimodal Landscape of Human Carotid Atherosclerosis

Background

Atherosclerotic plaques are complex tissues composed of a heterogeneous mixture of cells. However, we have limited understanding of the comprehensive transcriptional and phenotypical landscape of the cells within these lesions.

Methods

To characterize the landscape of human carotid atherosclerosis in greater detail, we combined cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) and single-cell RNA sequencing (scRNA-seq) to classify all cell types within lesions (n=21; 13 symptomatic) to achieve a comprehensive multimodal understanding of the cellular identities of atherosclerosis and their association with clinical pathophysiology.

Results

We identified 25 distinct cell populations each having a unique multi-omic signature, including macrophages, T cells, NK cells, mast cells, B cells, plasma cells, neutrophils, dendritic cells, endothelial cells, fibroblasts, and smooth muscle cells (SMCs). Within the macrophage populations, we identified 2 proinflammatory subsets that were enriched in IL1B or C1Q expression, 2 distinct TREM2 positive foam cell subsets, one of which also expressed inflammatory genes, as well as subpopulations displaying a proliferative gene expression signature and one expressing SMC-specific genes and upregulation of fibrotic pathways. An in-depth characterization uncovered several subsets of SMCs and fibroblasts, including a SMC-derived foam cell. We localized this foamy SMC to the deep intima of coronary atherosclerotic lesions. Using CITE-seq data, we also developed the first flow cytometry panel, using cell surface proteins CD29, CD142, and CD90, to isolate SMC-derived cells from lesions. Last, we found that the proportion of efferocytotic macrophages, classically activated endothelial cells, contractile and modulated SMC-derived cell types were reduced, and inflammatory SMCs were enriched in plaques of clinically symptomatic vs. asymptomatic patients.

Conclusions

Our multimodal atlas of cell populations within atherosclerosis provides novel insights into the diversity, phenotype, location, isolation, and clinical relevance of the unique cellular composition of human carotid atherosclerosis. This facilitates both the mapping of cardiovascular disease susceptibility loci to specific cell types as well as the identification of novel molecular and cellular therapeutic targets for treatment of the disease.

Pyruvate dehydrogenase kinase regulates vascular inflammation in atherosclerosis and increases cardiovascular risk

AIMS

Recent studies have revealed a close connection between cellular metabolism and the chronic inflammatory process of atherosclerosis. While the link between systemic metabolism and atherosclerosis is well established, the implications of altered metabolism in the artery wall are less understood. Pyruvate dehydrogenase kinase (PDK)-dependent inhibition of pyruvate dehydrogenase (PDH) has been identified as a major metabolic step regulating inflammation. Whether the PDK/PDH axis plays a role in vascular inflammation and atherosclerotic cardiovascular disease remains unclear.

METHODS AND RESULTS

Gene profiling of human atherosclerotic plaques revealed a strong correlation between PDK1 and PDK4 transcript levels and the expression of pro-inflammatory and destabilizing genes. Remarkably, the PDK1 and PDK4 expression correlated with a more vulnerable plaque phenotype, and PDK1 expression was found to predict future major adverse cardiovascular events. Using the small-molecule PDK inhibitor dichloroacetate (DCA) that restores arterial PDH activity, we demonstrated that the PDK/PDH axis is a major immunometabolic pathway, regulating immune cell polarization, plaque development, and fibrous cap formation in Apoe-/- mice. Surprisingly, we discovered that DCA regulates succinate release and mitigates its GPR91-dependent signals promoting NLRP3 inflammasome activation and IL-1β secretion by macrophages in the plaque.

CONCLUSIONS

We have demonstrated for the first time that the PDK/PDH axis is associated with vascular inflammation in humans and particularly that the PDK1 isozyme is associated with more severe disease and could predict secondary cardiovascular events. Moreover, we demonstrate that targeting the PDK/PDH axis with DCA skews the immune system, inhibits vascular inflammation and atherogenesis, and promotes plaque stability features in Apoe-/- mice. These results point toward a promising treatment to combat atherosclerosis.

Exploring the genetic basis of coronary artery disease using functional genomics

Genome-wide Association Studies (GWAS) have identified more than 300 loci associated with coronary artery disease (CAD), defining the genetic risk map of the disease. However, the translation of the association signals into biological-pathophysiological mechanisms constitute a major challenge. Through a group of examples of studies focused on CAD, we discuss the rationale, basic principles and outcomes of the main methodologies implemented to prioritize and characterize causal variants and their target genes. Additionally, we highlight the strategies as well as the current methods that integrate association and functional genomics data to dissect the cellular specificity underlying the complexity of disease mechanisms. Despite the limitations of existing approaches, the increasing knowledge generated through functional studies helps interpret GWAS maps and opens novel avenues for the clinical usability of association data.

RNA-binding proteins in vascular inflammation and atherosclerosis

Atherosclerotic cardiovascular disease (ASCVD) remains the major cause of premature death and disability worldwide, even when patients with an established manifestation of atherosclerotic heart disease are optimally treated according to the clinical guidelines. Apart from the epigenetic control of transcription of the genetic information to messenger RNAs (mRNAs), gene expression is tightly controlled at the post-transcriptional level before the initiation of translation. Although mRNAs are traditionally perceived as the messenger molecules that bring genetic information from the nuclear DNA to the cytoplasmic ribosomes for protein synthesis, emerging evidence suggests that processes controlling RNA metabolism, driven by RNA-binding proteins (RBPs), affect cellular function in health and disease. Over the recent years, vascular endothelial cell, smooth muscle cell and immune cell RBPs have emerged as key co- or post-transcriptional regulators of several genes related to vascular inflammation and atherosclerosis. In this review, we provide an overview of cell-specific function of RNA-binding proteins involved in all stages of ASCVD and how this knowledge may be used for the development of novel precision medicine therapeutics.

Dissecting the polygenic basis of atherosclerosis via disease-associated cell state signatures

Coronary artery disease (CAD) is a pandemic disease where up to half of the risk is explained by genetic factors. Advanced insights into the genetic basis of CAD require deeper understanding of the contributions of different cell types, molecular pathways, and genes to disease heritability. Here, we investigate the biological diversity of atherosclerosis-associated cell states and interrogate their contribution to the genetic risk of CAD by using single-cell and bulk RNA sequencing (RNA-seq) of mouse and human lesions. We identified 12 disease-associated cell states that we characterized further by gene set functional profiling, ligand-receptor prediction, and transcription factor inference. Importantly, Vcam1+ smooth muscle cell state genes contributed most to SNP-based heritability of CAD. In line with this, genetic variants near smooth muscle cell state genes and regulatory elements explained the largest fraction of CAD-risk variance between individuals. Using this information for variant prioritization, we derived a hybrid polygenic risk score (PRS) that demonstrated improved performance over a classical PRS. Our results provide insights into the biological mechanisms associated with CAD risk, which could make a promising contribution to precision medicine and tailored therapeutic interventions in the future.

The Genetic Determinants of Aortic Distention

BACKGROUND

As the largest conduit vessel, the aorta is responsible for the conversion of phasic systolic inflow from ventricular ejection into more continuous peripheral blood delivery. Systolic distention and diastolic recoil conserve energy and are enabled by the specialized composition of the aortic extracellular matrix. Aortic distensibility decreases with age and vascular disease.

OBJECTIVES

In this study, we sought to discover epidemiologic correlates and genetic determinants of aortic distensibility and strain.

METHODS

We trained a deep learning model to quantify thoracic aortic area throughout the cardiac cycle from cardiac magnetic resonance images and calculated aortic distensibility and strain in 42,342 UK Biobank participants.

RESULTS

Descending aortic distensibility was inversely associated with future incidence of cardiovascular diseases, such as stroke (HR: 0.59 per SD; P = 0.00031). The heritabilities of aortic distensibility and strain were 22% to 25% and 30% to 33%, respectively. Common variant analyses identified 12 and 26 loci for ascending and 11 and 21 loci for descending aortic distensibility and strain, respectively. Of the newly identified loci, 22 were not significantly associated with thoracic aortic diameter. Nearby genes were involved in elastogenesis and atherosclerosis. Aortic strain and distensibility polygenic scores had modest effect sizes for predicting cardiovascular outcomes (delaying or accelerating disease onset by 2%-18% per SD change in scores) and remained statistically significant predictors after accounting for aortic diameter polygenic scores.

CONCLUSIONS

Genetic determinants of aortic function influence risk for stroke and coronary artery disease and may lead to novel targets for medical intervention.

An automatic entropy method to efficiently mask histology whole-slide images

Tissue segmentation of histology whole-slide images (WSI) remains a critical task in automated digital pathology workflows for both accurate disease diagnosis and deep phenotyping for research purposes. This is especially challenging when the tissue structure of biospecimens is relatively porous and heterogeneous, such as for atherosclerotic plaques. In this study, we developed a unique approach called 'EntropyMasker' based on image entropy to tackle the fore- and background segmentation (masking) task in histology WSI. We evaluated our method on 97 high-resolution WSI of human carotid atherosclerotic plaques in the Athero-Express Biobank Study, constituting hematoxylin and eosin and 8 other staining types. Using multiple benchmarking metrics, we compared our method with four widely used segmentation methods: Otsu's method, Adaptive mean, Adaptive Gaussian and slideMask and observed that our method had the highest sensitivity and Jaccard similarity index. We envision EntropyMasker to fill an important gap in WSI preprocessing, machine learning image analysis pipelines, and enable disease phenotyping beyond the field of atherosclerosis.

Dipeptidyl Peptidase 4/Midline-1 Axis Promotes T Lymphocyte Motility in Atherosclerosis

T cells play a crucial role in atherosclerosis, with its infiltration preceding the formation of atheroma. However, how T-cell infiltration is regulated in atherosclerosis remains largely unknown. Here, this work demonstrates that dipeptidyl peptidase-4 (DPP4) is a novel regulator of T-cell motility in atherosclerosis. Single-cell ribonucleic acid (RNA) sequencing and flow cytometry show that CD4+ T cells in atherosclerotic patients display a marked increase of DPP4. Lack of DPP4 in hematopoietic cells or T cells reduces T-cell infiltration and atherosclerotic plaque volume in atherosclerosis mouse models. Mechanistically, DPP4 deficiency reduces T-cell motility by suppressing the expression of microtubule associated protein midline-1 (Mid1) in T cells. Deletion of either DPP4 or Mid1 inhibits chemokine-induced shape change and motility, while restitution of Mid1 in Dpp4-/- T cell largely restores its migratory ability. Thus, DPP4/Mid1, as a novel regulator of T-cell motility, may be a potential inflammatory target in atherosclerosis.

Diversity of arterial cell and phenotypic heterogeneity induced by high-fat and high-cholesterol diet

Lipid metabolism disorder is the basis of atherosclerotic lesions, in which cholesterol and low-density lipoprotein (LDL) is the main factor involved with the atherosclerotic development. A high-fat and high-cholesterol diet can lead to this disorder in the human body, thus accelerating the process of disease. The development of single-cell RNA sequencing in recent years has opened the possibility to unbiasedly map cellular heterogeneity with high throughput and high resolution; alterations mediated by a high-fat and high-cholesterol diet at the single-cell transcriptomic level can be explored with this mean afterward. We assessed the aortic arch of 16-week old Apoe-/- mice of two control groups (12 weeks of chow diet) and two HFD groups (12 weeks of high fat, high cholesterol diet) to process single-cell suspension and use single-cell RNA sequencing to anatomize the transcripts of 5,416 cells from the control group and 2,739 from the HFD group. Through unsupervised clustering, 14 cell types were divided and defined. Among these cells, the cellular heterogeneity exhibited in endothelial cells and immune cells is the most prominent. Subsequent screening delineated ten endothelial cell subsets with various function based on gene expression profiling. The distribution of endothelial cells and immune cells differs significantly between the control group versus the HFD one. The existence of pathways that inhibit atherosclerosis was found in both dysfunctional endothelial cells and foam cells. Our data provide a comprehensive transcriptional landscape of aortic arch cells and unravel the cellular heterogeneity brought by a high-fat and high-cholesterol diet. All these findings open new perspectives at the transcriptomic level to studying the pathology of atherosclerosis.

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.

Study Design Overview

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.

TXNIP Suppresses the Osteochondrogenic Switch of Vascular Smooth Muscle Cells in Atherosclerosis

BACKGROUND

The osteochondrogenic switch of vascular smooth muscle cells (VSMCs) is a pivotal cellular process in atherosclerotic calcification. However, the exact molecular mechanism of the osteochondrogenic transition of VSMCs remains to be elucidated. Here, we explore the regulatory role of TXNIP (thioredoxin-interacting protein) in the phenotypical transitioning of VSMCs toward osteochondrogenic cells responsible for atherosclerotic calcification.

METHODS

The atherosclerotic phenotypes of Txnip^-/- mice were analyzed in combination with single-cell RNA-sequencing. The atherosclerotic phenotypes of Tagln-Cre; Txnip^flox/flox mice (smooth muscle cell-specific Txnip ablation model), and the mice transplanted with the bone marrow of Txnip^-/- mice were analyzed. Public single-cell RNA-sequencing dataset (GSE159677) was reanalyzed to define the gene expression of TXNIP in human calcified atherosclerotic plaques. The effect of TXNIP suppression on the osteochondrogenic phenotypic changes in primary aortic VSMCs was analyzed.

RESULTS

Atherosclerotic lesions of Txnip^-/- mice presented significantly increased calcification and deposition of collagen content. Subsequent single-cell RNA-sequencing analysis identified the modulated VSMC and osteochondrogenic clusters, which were VSMC-derived populations. The osteochondrogenic cluster was markedly expanded in Txnip^-/- mice. The pathway analysis of the VSMC-derived cells revealed enrichment of bone- and cartilage-formation-related pathways and bone morphogenetic protein signaling in Txnip^-/- mice. Reanalyzing public single-cell RNA-sequencing dataset revealed that TXNIP was downregulated in the modulated VSMC and osteochondrogenic clusters of human calcified atherosclerotic lesions. Tagln-Cre; Txnip^flox/flox mice recapitulated the calcification and collagen-rich atherosclerotic phenotypes of Txnip^-/- mice, whereas the hematopoietic deficiency of TXNIP did not affect the lesion phenotype. Suppression of TXNIP in cultured VSMCs accelerates osteodifferentiation and upregulates bone morphogenetic protein signaling. Treatment with the bone morphogenetic protein signaling inhibitor K02288 abrogated the effect of TXNIP suppression on osteodifferentiation.

CONCLUSIONS

Our results suggest that TXNIP is a novel regulator of atherosclerotic calcification by suppressing bone morphogenetic protein signaling to inhibit the transition of VSMCs toward an osteochondrogenic phenotype.

Translational opportunities of single-cell biology in atherosclerosis

The advent of single-cell biology opens a new chapter for understanding human biological processes and for diagnosing, monitoring, and treating disease. This revolution now reaches the field of cardiovascular disease (CVD). New technologies to interrogate CVD samples at single-cell resolution are allowing the identification of novel cell communities that are important in shaping disease development and direct towards new therapeutic strategies. These approaches have begun to revolutionize atherosclerosis pathology and redraw our understanding of disease development. This review discusses the state-of-the-art of single-cell analysis of atherosclerotic plaques, with a particular focus on human lesions, and presents the current resolution of cellular subpopulations and their heterogeneity and plasticity in relation to clinically relevant features. Opportunities and pitfalls of current technologies as well as the clinical impact of single-cell technologies in CVD patient care are highlighted, advocating for multidisciplinary and international collaborative efforts to join the cellular dots of CVD.

New Insights into Cerebral Vessel Disease Landscapes at Single-Cell Resolution: Pathogenetic and Therapeutic Perspectives

Cerebrovascular diseases are a leading cause of death and disability globally. The development of new therapeutic targets for cerebrovascular diseases (e.g., ischemic, and hemorrhagic stroke, vascular dementia) is limited by a lack of knowledge of the cellular and molecular biology of health and disease conditions and the factors that cause injury to cerebrovascular structures. Here, we describe the role of advances in omics technology, particularly RNA sequencing, in studying high-dimensional, multifaceted profiles of thousands of individual blood and vessel cells at single-cell resolution. This analysis enables the dissection of the heterogeneity of diseased cerebral vessels and their atherosclerotic plaques, including the microenvironment, cell evolutionary trajectory, and immune response pathway. In animal models, RNA sequencing permits the tracking of individual cells (including immunological, endothelial, and vascular smooth muscle cells) that compose atherosclerotic plaques and their alteration under experimental settings such as phenotypic transition. We describe how single-cell RNA transcriptomics in humans allows mapping to the molecular and cellular levels of atherosclerotic plaques in cerebral arteries, tracking individual lymphocytes and macrophages, and how these data can aid in identifying novel immune mechanisms that could be exploited as therapeutic targets for cerebrovascular diseases. Single-cell multi-omics approaches will likely provide the unprecedented resolution and depth of data needed to generate clinically relevant cellular and molecular signatures for the precise treatment of cerebrovascular diseases.

Opportunities and Challenges in Understanding Atherosclerosis by Human Biospecimen Studies

Over the last few years, new high-throughput biotechnologies and bioinformatic methods are revolutionizing our way of deep profiling tissue specimens at the molecular levels. These recent innovations provide opportunities to advance our understanding of atherosclerosis using human lesions aborted during autopsies and cardiac surgeries. Studies on human lesions have been focusing on understanding the relationship between molecules in the lesions with tissue morphology, genetic risk of atherosclerosis, and future adverse cardiovascular events. This review will highlight ways to utilize human atherosclerotic lesions in translational research by work from large cardiovascular biobanks to tissue registries. We will also discuss the opportunities and challenges of working with human atherosclerotic lesions in the era of next-generation sequencing.

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.