Introduction

Atherosclerotic plaque rupture precipitates a thromboembolic event leading to ischemia. Clinical examples include acute neurological events, such as transient ischemic attacks and stroke, in the case of carotid arteries, and acute coronary syndromes in the case of coronary arteries1,2. The development of an atherosclerotic plaque is characterized by chronic inflammation, with the recruitment of monocyte-derived macrophages and T cells. These produce inflammatory mediators that enhance inflammation and promote vascular smooth muscle cells (VSMC) to proliferate and secrete collagen to form the fibrous cap3. In the late stage of atherosclerosis, thinning of the fibrous cap due to increased apoptosis of VSMCs and expression of matrix metalloproteinases create a rupture-prone thin cap fibroatheroma4.

Atherosclerotic plaques consist of a soft lipid-rich atheromatous core covered by a firm, collagen-rich fibrous cap which contains sclerotic tissue, VSMCs, and an intact endothelial lining5. Ruptured plaques have thinner overlying fibrous caps infiltrated by macrophages, T-lymphocytes, and B-lymphocytes2. The molecular mechanisms driving asymptomatic plaques to transition to those prone to rupture and cause atheroemboli remain unclear and require further study.

The features of the ruptured plaque have been primarily determined through histological analysis of tissue obtained from patients exhibiting symptoms within the past 180 days. While these data are quite valuable, our previous works suggest notable transient changes in RNA expression that occur with plaque rupture6. Here, we compared the transcriptomes of asymptomatic and recently (within 5 days) ruptured carotid plaques to determine molecular mechanisms active at the time of rupture.

Results

Patient population

Carotid plaques were obtained from 11 patients undergoing CEA, 5 were asymptomatic and 6 had an ischemic cerebrovascular event within the past 5 days (“recently ruptured”, Table 1). There was no significant difference in age, body mass index, lipid profile, serum creatinine, or estimated glomerular filtration rates between the two patient groups. Additionally, there similar rates of smoking, alcohol use and chronic kidney disease in both groups. All subjects were hypertensive while none were diabetic. Each group had similar levels of aspirin, statin, and clopidogrel use. The percent stenosis was similar between the two groups. The recently ruptured group consisted of 2 patients with crescendo transient ischemic attacks and 4 with NIH Stroke Severity Scores of 1, 3, 4, and 7. The time between the cerebrovascular event and the CEA of 3.5 ± 0.5 days.

Table 1 Patient characteristics.

Principal component analysis of highly variable transcripts identifies a principal component associated with rupture

PCA of the 1750 most variable transcripts found that the first principal component, which accounts for 26% of the total variability, represents the differences between the samples from asymptomatic versus recently ruptured groups (Fig. 1A). Gene ontology analysis of the transcripts with loadings associated with the recently symptomatic group, showed that this group was enriched for transcripts related to phagocytosis, B cell function, complement activation, inflammation and leukocyte chemotaxis (Fig. 1B). The loadings associated with the asymptomatic group were characterized by increased plasma lipoprotein particle clearance and remodeling and transcripts associated with regulation of muscle contraction (Fig. 1C).

Figure 1
figure 1

Principal component analysis (PCA) indicates plaque rupture is associated with alterations in the transcriptome. (A) PCA plot of the top 1750 most variable transcripts in asymptomatic and recently ruptured plaques. (B) Gene ontology of the transcripts associated with recently ruptured plaques. (C) Gene ontology of the transcripts associated with asymptomatic plaques.

Hierarchical clustering of highly variable transcripts separates asymptomatic and recently symptomatic

To more directly examine the changes in gene expression that occur with plaque rupture, we performed hierarchical clustering of the transcripts in the carotid plaques (Fig. 2). Clustering a minimum of the 400 most variable transcripts samples achieved separation between the asymptomatic and recently ruptured plaque samples. Interestingly, while the asymptomatic samples clustered tightly together, two of the recently ruptured samples clustered away from all the other samples, suggesting greater variability in the recently ruptured plaque samples. We compared the individual demographic data for each patient with the patient clustering to determine if the patient characteristics might explain the groupings, but no clinical factor correlated with the sample dendrogram.

Figure 2
figure 2

Hierarchical clustering identifies groups of co-expressed transcripts that distinguish recently ruptured plaques from asymptomatic plaques. Clustering of the top 400 transcripts in carotid plaques differentiates between asymptomatic and recently ruptured plaques. Bars adjacent to the cluster indicate locations of transcripts associated with B cell function (green), metalloproteinases (blue), and interferon responses (red). Supplementary Fig. 2 presents this data with the transcript names included.

Overall, the heatmap does not present distinctive clusters of transcripts unique to asymptomatic or recently ruptured. Close inspection of the heatmap reveals a cluster of immunoglobulins and B cell function transcripts where expression is increased in the recently ruptured group. This was also true for clusters that include matrix metalloproteinases (e.g., MMP1, MMP7, MMP8, MMP9, MMP12, MMP13, and MMP25) and responses to interferon (e.g., IFI6, IFI27, IFIT1, IFIT3, IFITM1, IFI44L, IRF7, ISG15, OAS2, OAS3, and OASL) though the magnitude of the difference is less pronounced. Transcripts encoding inflammatory genes were scattered throughout the clusters, including numerous cytokines (e.g., IL6, CXCL1, CXCL5, CXCL6, CXCL8, CXCL9, and CCL18) exhibiting increases in the recently ruptured group.

Differential gene expression underlying increased atherosclerotic plaque vulnerability and rupture

A total of 422 DEGs were identified from atherosclerotic plaque samples of recently ruptured compared to asymptomatic patients (p-adj < 0.05). Two hundred sixty-three transcripts were significantly increased, and 159 transcripts were significantly decreased (Supplementary Table 1). To better understand the biological processes behind the transition of stable to vulnerable atherosclerotic plaques, we explored significant transcripts of interest in more detail (Table 2).

Table 2 Differentially expressed genes of interest.

Recently ruptured atherosclerotic plaques demonstrated increased expression of pro-inflammatory genes involved in the mobilization and recruitment of leukocytes to the vessel wall. XCR1, which encodes a chemokine receptor, and CD177, which mediates neutrophil activation, were increased in the recently ruptured group. The recently ruptured group also exhibited elevated levels of several transcripts associated with B cell function, including numerous immunoglobins and transcripts related to B cell proliferation (MZB1) and activation (CD79A, SH2D3C, and ZAP70).

Recently ruptured atherosclerotic plaque samples were also found to have a change in transcript expression that may relate to thinning of the fibrous cap through decreased VSMC proliferation and migration and increased apoptosis. SDC1 and SIK1 are associated with decreased VSMC proliferation and migration. Additionally, the recently ruptured plaques had increased expression of transcripts related to inhibition of cell proliferation (NPDC1, DOK3, and ZBTB17). Reduced expression of a pro-migration transcript (TMSB15B) was also observed in the recently ruptured samples. Several transcripts encode Ras GTPase-activating proteins (RASA4, RASA4B, ARHGAP4, HMHA1, and SH3BP1) that promote inactivation of Rac1 and cdc42, inhibiting proliferation and migration.

In addition to its anti-proliferative functions, DOK3 is associated with increased apoptosis. Recently ruptured samples exhibited an increase in a pro-apoptotic transcript (TP63) and decreases in a transcript encoding an inhibitor of apoptosis (MTRNR2L13). Additionally, ADAMTS4, which encodes a proteinase able to degrade the extracellular matrix and has been identified as a potential marker of carotid plaque vulnerability, was upregulated in the recently ruptured samples7. The long non-coding RNA, MIAT, was significantly increased in the recently ruptured samples. This agrees with a previous report that MIAT is increased in ruptured compared to stable carotid plaques; however, MIAT promotes proliferation and inhibits apoptosis8.

Discussion

The RNA-Seq dataset presented here demonstrates that recently ruptured plaques have a unique transcriptome that includes inflammatory and anti-proliferative features. The PCA analysis revealed that rupture status is the most significant contributor to the variability between samples. The hierarchical clustering demonstrated that the transcriptomes of the plaques share similarities according to plaque status as well. While expected, this common expression profile supports that the development of plaque vulnerability is an active process mediated by significant changes in gene expression.

Papaspyridonos et al. performed a similar analysis where CEA samples were collected and assessed macroscopically to stratify the samples into stable and unstable and a supervised analysis was used to identify DEGs between these groups. This study found 27 DEGs of interest that were confirmed by PCR9. Our study found Fold Changes in the same direction as in their study for 18 of these mRNAs. These similarities suggest that our dataset does include similar changes as seen in unstable plaques that have not necessarily ruptured.

Examination of the transcripts underlying the clustering of the samples revealed that many transcripts associated with the recently ruptured plaques play a role in inflammation, especially B cell function and type I interferon responses. B cells have been found in the fibrous caps of ruptured carotid plaques2. Analysis of the DEGs further supported an essential role for B cells in plaque rupture, with significant increases in transcripts associated with B cell function observed in the recently ruptured plaques. The role of B cells in atherosclerosis is dependent on subtype. B1a lymphocytes reduce the necrotic core size and thus are atheroprotective10. In contrast, the B2 subset promotes plaque development and vulnerability, potentially through secretion of tumor necrosis factor-α and increased apoptosis11,12. While many immunoglobulin transcripts were increased in the recently ruptured plaques, the lack of a significant increase in the transcript encoding IgM suggests that the B2 subset may predominate. Overall, these findings play an essential role for B2 cells in developing plaque vulnerability and rupture.

In addition to the B cells, our data also support enhanced inflammation. Type I interferons can be produced by multiple cell types in atherosclerotic plaques, including B cells, and promote foam cell formation and amplify inflammation. The hierarchical clustering indicated a set of co-expressed interferon response genes. IFIT1 and IFIT3 have been implicated in pro-inflammatory polarization of macrophages and decreased collagen deposition leading to plaque vulnerability13,14. Loss of OAS2 and OAS3 expression has been linked to reduced atherosclerotic plaque development in a mouse model of atherosclerosis15. Both the PCA and hierarchical clustering pointed to increased cytokine signaling. Overall, the data suggest an increase in inflammation associated with plaque rupture.

Our data also points to an increase in mechanisms that promote thinning of the fibrous cap. An enrichment of transcripts involved with matrix disassembly was associated with the recently ruptured group, and several matrix metalloproteinases were co-expressed in the hierarchical clustering. Similar to previous reports7, ADAMTS4 was reduced in the recently ruptured group. Loss of ADAMTS4 is associated with increased collagen content and overall plaque stability in a mouse model of atherosclerosis. The matrix metalloproteinases co-expressed in the hierarchical clustering target a broad spectrum of substrates, including collagens (MMP1, MMP8, and MMP13), gelatins (MMP9), elastins (MMP12), and fibronectin (MMP7)16,17. In a mouse model of plaque vulnerability, MMP9 was associated with greater plaque stability, while MMP7 attenuated VSMC content and MMP12 promoted destabilization18. In contrast, overexpression of MMP9 has been proposed to promote plaque destabilization4.

Several transcripts associated with decreased proliferation and migration were also upregulated in the recently ruptured plaques. Increased SDC1 is associated with non-proliferating differentiated VSMCs, and SIK1 is associated with decreased vascular remodeling and was increased in the ruptured samples19,20. There was also upregulation of several Ras GTPase-activating proteins. Members of this protein family have been implicated in the development of plaque vulnerability21,22. Additional transcripts suggest decreased proliferation and increased apoptosis occur at the time of rupture. Thinning of the fibrous cap occurs through decreased VSMC proliferation and increased VSMC apoptosis23.

It is important to note that this study measured RNA from whole artery homogenates. Thus, we cannot attribute an individual change in RNA expression to a specific cell type. Likewise, the observed changes in RNA expression may represent changes in the cellular composition of the plaque. Indeed, our data suggest an increase in B cell infiltration occurred in the recently ruptured plaques. Additionally, the samples were collected 2 to 5 days post-rupture and therefore include changes associated with the response to the rupture. We cannot dissect those changes present at the time of rupture from those that occur after rupture. Finally, additional studies are required to confirm that changes in protein expression accompany the changes in transcript abundance we report here. To our knowledge, however, this is the first study to perform RNA-Seq on carotid samples collected this close to the plaque rupture. Most studies comparing CEA samples include symptomatic plaques collected greater than 100 days post-rupture. Previously, we have reported that miR-221 and -222 are down-regulated in recently ruptured plaques but increase to levels similar to that found in asymptomatic plaques in 7 days6. Thus, these data likely include changes in RNA expression that have previously been missed. Previous histologic analyses of carotid plaque rupture found "scarce" B cell infiltration22. Our data support a robust role for B cells in plaque rupture.

Conclusion

Overall, these data present a picture of significant changes in the inflammation within the plaque associated with rupture, including a vital role for B cells. Additionally, an increase in a broad spectrum of matrix metalloproteinases coupled with the loss of VSMCs underlies thinning of the fibrous cap. Future studies on changes in RNA expression that can be localized to specific regions or cells within the plaque will further enhance our understanding of the molecular mechanisms driving plaque rupture.

Materials and methods

Sample collection

Carotid plaque specimens were obtained from a total of 11 patients who underwent carotid endarterectomy (CEA) in the Section of Vascular/Endovascular Surgery in the Department of Surgery at the Ochsner Clinic in New Orleans, LA. Samples obtained from patients undergoing CEA were stratified into two distinct clinical phenotypes: asymptomatic and recently ruptured. The asymptomatic group included patients undergoing CEA without a previous neurologic event but high-grade carotid stenosis (n = 5). The recently ruptured group (n = 6) consists of patients undergoing CEA within 5 days of an acute neurologic event (transient ischemic attack or acute stroke). The samples were immediately snap-frozen and stored in liquid nitrogen until processing for RNA isolation.

All procedures in this study were performed in accordance with the ethical standards of the 1964 Helsinki Declaration. This study was approved by the Ochsner Health institutional review board (IRB protocol # 2010.130B). All of the participants gave written informed consent for their participation in the study.

RNA isolation

Carotid plaques were sectioned into 6 transverse segments of ~ 1 cm in length and macroscopic examination was performed to determine the sections with greatest plaque content. Sections containing a complicated lesion were homogenized in TRIzol Reagent (ThermoFisher) using a rotor stator homogenizer and total RNA was isolated according to the manufacturer's instructions. Total RNA was further purified using the miRNeasy Mini Kit (Qiagen).

RNA-sequencing

The University Wisconsin Biotechnology Center performed library preparation and sequencing. Ribosomal RNA was reduced using the RiboZero assay (EpiCentre-Illumina), and an RNA-Seq library was constructed with the stranded TruSeq RNA kit (Illumina). Acceptable quality of each RNA sample and DNA library was confirmed using the 2100 Bioanalyzer system. Using a HiSeq2000 (Illumina), 1 × 100 bp sequencing was performed at average depth of 42 M reads per sample. The Tulane Next Generation Sequence Analysis Core aligned RNA-Seq reads to the human genome (hg19 assembly) using Novoalign (Novocraft), and splice junctions were identified using TopHat24. Coverage files and FPKM values were generated using SAMMate, and coverage was visualized on the Integrated Genome Viewer25,26.

Transcript abundance estimates were imported using tximport and package and principal component analysis (PCA), hierarchical clustering, and identification of differentially expressed genes (DEGs) were performed in an unsupervised manner using the DESeq2 package27,28. DEGs were defined as those transcripts with a p-adjusted value < 0.05 with a fold change > 2. Gene ontology was performed using the PANTHER classification system with a minimum enrichment value of 1 and false discovery rate < 0.01229. The datasets generated and/or analysed during the current study are available in the [NCBI GEO (Gene Expression Omnibus)] repository, [PERSISTENT WEB INK OR GSE198600 TO DATASETS].

Validation of RNA-Seq dataset with Droplet Digital PCR

The results of the RNA sequencing were validated by comparing the log2(Fold Change) obtained from the RNA-Seq analysis to that measured by polymerase chain reaction (PCR) of mRNAs selected randomly from a preliminary list of DEGs in an additional set of asymptomatic and recently ruptured samples (Supplementary Fig. 1), n = 5–16/group). The appropriate Quantitect Primer Assays (Qiagen) was paired with QX200 ddPCR EvaGreen PCR supermix (BioRad) to perform a one-step RT-PCR. QX200 droplet reader (BioRad) paired with QuantaSoft Analysis Pro v1.0 software (BioRad) for direct quantification of individual target molecules from each sample.