- Research
- Open access
- Published:
Identification of hub genes and regulatory networks in histologically unstable carotid atherosclerotic plaque by bioinformatics analysis
BMC Medical Genomics volume 15, Article number: 145 (2022)
Abstract
Objective
This study identified underlying genetic molecules associated with histologically unstable carotid atherosclerotic plaques through bioinformatics analysis that may be potential biomarkers and therapeutic targets.
Methods
Three transcriptome datasets (GSE41571, GSE120521 and E-MTAB-2055) and one non-coding RNA dataset (GSE111794) that met histological grouping criteria of unstable plaque were downloaded. The common differentially expressed genes (co-DEGs) of unstable plaques identified from three mRNA datasets were annotated by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomics (KEGG). A protein–protein interaction (PPI) network was constructed to present the interaction between co-DEGs and screen out hub genes. MiRNet database and GSE111794 dataset were used to identify the miRNAs targeting hub genes. Associated transcription factors (TFs) and drugs were also predicted. These predicted results were used to construct miRNA/TFs-hub gene and drug-hub gene regulatory networks.
Results
A total of 105 co-DEGs were identified, including 42 up-regulated genes and 63 down-regulated genes, which were mainly enriched in collagen-containing extracellular matrix, focal adhesion, actin filament bundle, chemokine signaling pathway and regulates of actin cytoskeleton. Ten hub genes (up-regulated: HCK, C1QC, CD14, FCER1G, LCP1 and RAC2;Â down-regulated: TPM1, MYH10, PLS3 and FMOD) were screened. HCK and RAC2 were involved in chemokine signaling pathway, MYH10 and RAC2 were involved in regulation of actin cytoskeleton. We also predicted 12 miRNAs, top5 TFs and 25 drugs targeting hub genes. In the miRNA/TF-hub gene regulatory network, PLS3 was the most connected hub genes and was targeted by six miRNAs and all five screened TFs. In the drug-hub gene regulatory network, HCK was targeted by 20 drugs including 10 inhibitors.
Conclusions
We screened 10 hub genes and predicted miRNAs and TFs targeting them. These molecules may play a crucial role in the progression of histologically unstable carotid plaques and serve as potential biomarkers and therapeutic targets.
Introduction
Carotid artery stenosis, most commonly caused by atherosclerosis, is a main reason of ischemic stroke [1]. However, there is increasing evidence that in addition to atherosclerotic stenosis, unstable plaques also play an important role in promoting symptomatic stroke [2, 3]. In advanced atherosclerosis, plaques may have a large necrotic lipid core, a weak fibrous cap, intraplaque neovascularization and hemorrhage. These histological changes make the plaque unstable and prone to rupture, which is more likely to result in thrombosis and embolization of plaque material, thereby leading to vascular occlusion and subsequent ischemic stroke [3, 4]. Moreover, unstable plaques cause more severe neurological damage in patients with acute cerebral infarction than stable plaques [5]. Therefore, clarifying the pathogenesis associated with unstable plaque formation is significant and may reveal potential biomarkers and therapeutic targets.
With the further research of the pathogenesis, bioinformatics analysis has become an essential technique for discovering genetic alteration in various diseases. By analyzing microarray or sequencing data, the differentially expressed genes (DEGs) of unstable carotid plaques can be found. The functions and pathways involved in DEGs were then predicted by functional enrichment analysis. Further construction of protein–protein interaction (PPI) network can understand their interactions and screen out the hub genes related to unstable carotid plaques. Finally, prediction of miRNAs, transcription factors (TFs) and drugs targeting hub genes can make the pathogenesis network more comprehensive.
Several unstable carotid plaque datasets based on clinical and/or histological criteria have been published. In fact, unstable carotid plaques identified by clinical criteria alone, also known as symptomatic carotid stenosis, do not necessarily meet the histological criteria [6]. Some studies have pooled the datasets of unstable plaques that meet symptomatic or histological criteria, which may produce inaccurate results. In this study, three transcriptome datasets and one non-coding RNA dataset met histological criteria of unstable carotid plaques were selected from the Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo/) database and the ArrayExpress (https://www.ebi.ac.uk/arrayexpress/) database for bioinformatics analysis, hoping to provide more insights into unstable carotid plaques.
Methods
The study flow was shown in Fig. 1.
Data source
Three gene expression profiles were downloaded. The matrix data of GSE41571 (microarray) and the Fragments Per Kilobase Million (FPKM) data of GSE120521 (RNA-sequencing) were obtained from the GEO database, and the processed data of E-MTAB-2055 (microarray) was downloaded from the ArrayExpress database. GSE41571, including 5 ruptured plaques and 6 stable plaques, was assayed at the platform of GPL570. GSE120521, containing 4 unstable plaques and 4 stable plaques, was measured at the platform of GPL16791. The platform of A-MEXP-931 was used by E-MTAB-2055 which included 25 ruptured plaques and 22 stable plaques. GSE111794 is a microRNA expression profile downloaded from the GEO database. It used GPL16384 platform to assay 9 unstable and 9 stable plaques. The local research ethics committee authorized construction of these three datasets.
Identification of DEGs
NetworkAnalyst [7] (https://www.networkanalyst.ca/), a comprehensive gene-centric platform supporting gene expression profiling, was used to identify the DEGs of GSE41571 and GSE120521. The DEGs of E-MTAB-2055 were screened directly from its processed data. DEGs were defined as genes with p-value < 0.05 and |log fold change (FC)|> 1. The screening results were visualized using volcano maps drawn by GraphPad Prism 9.0.0. An online drawing tool (http://bioinformatics.psb.ugent.be/webtools/Venn/) was used to select the common DEGs (co-DEGs) of three datasets and draw a Venn diagram.
Functional enrichment analysis
To further investigate the potential functions of these co-DEGs, the list of co-DEGs was uploaded to Metascape [8] (https://metascape.org/) for Gene Ontology (GO) [9] and Kyoto Encyclopedia of Genes and Genomes (KEGG) [10] pathway enrichment analysis. Metascape is a web-based portal leveraging over 40 independent knowledgebases to provide a comprehensive gene list annotation and analysis resource for users [8]. Cellular compartments (CC), biological processes (BP), molecular function (MF) and KEGG pathway were selected for analysis. The threshold for P-value was set at 0.05 and the minimum enrichment score was 1.5. The visualization results were generated by Metascape and GraphPad Prism 9.0.0.
PPI network and hub gene identification
The Search Tool for the Retrieval of Interacting Genes [11] (STRING, v11.5, https://cn.string-db.org/) is a database of known and predicted PPIs. After uploading the co-DEGs list to the STRING database, PPI results were obtained with required interaction score > 0.4. Cytoscape (v3.9.0) [12] was then used to construct a visual PPI network and identify the hub genes. In order to avoid inaccurate screening results caused by using a single rank method, maximal clique centrality (MCC), maximum neighborhood component (MNC), node connect degree (Degree), edge percolated component (EPC), closeness and radiality rank methods were used for screening top 10 genes, respectively. All of these rank methods are included in Cytoscape’s application named Cytohubba (v0.1) [13]. Their results were integrated, then the genes with the top 10 frequencies were identified as the hub genes.
Prediction of miRNA/TF‑hub gene regulatory network
The miRNAs targeting hub genes were predicted using miRNet [14] (v2.0, https://www.mirnet.ca/), which is a miRNA-centric network visual analytics platform. The miRNA-hub gene data were collected from the well-annotated database miRTarBase v8.0 [15], DIANA-TarBase v8.0 [16], and miRecords [17]. The differentially expressed miRNAs of GSE111794 were analyzed using GEO2R, and then intersected with the predicted results of miRNet 2.0 to identify meaningful miRNA. TFs were also predicted in miRNet, using JASPAR database[18] resources. The top5 TFs were selected to construct the regulatory network. Finally, we visualized the miRNA/TF-hub gene regulatory network using Cytoscape 3.9.0.
Prediction of drug-hub gene regulatory network
The drugs targeting all hub genes were predicted using the Drug Gene Interaction Database [19] (DGIdb, v4.2.0, https://dgidb.genome.wustl.edu/). DGIdb used a combination of expert curation and text-mining to mine drug-gene interactions from DrugBank [20], PharmGKB [21], Drug Target Commons [22] and others. The Cytoscape3.9.0 was used to construct the drug–gene regulatory network.
Results
Identification of DEGs
DEGs between unstable and stable carotid plaques were identified. Seven hundred and ninety-six DEGs were found in GSE41571, including 266 up-regulated and 530 down-regulated DEGs (Fig. 2A). There were 796 DEGs in GSE41571, including 266 up-regulated and 530 down-regulated DEGs (Fig. 2B). And we obtained 418 DEGs in E-MTAB-2055, including 222 up-regulated and 196 down-regulated DEGs (Fig. 2C). After integration, a total of 105 co-DEGs including 42 up-regulated and 63 down-regulated DEGs overlapped in three datasets (Fig. 2D, E).
GO and KEGG enrichment analysis
All 105 co-DEGs were uploaded to Metascape database for functional enrichment analysis. The first 20 representative enriched terms (one for each cluster) are shown in Fig. 3A and Additional file 1: Table S1. To further capture the relationships between these terms, Metascape presented a network where terms with a similarity > 0.3 are connected by edges (Fig. 3B, C). The top 3 enriched terms were collagen-containing extracellular matrix, focal adhesion and actin filament bundle, all of which belong to CC (Figs. 3A, 4A). As summarized in Table 1, co-DEGs were mostly enriched in BP of extracellular matrix organization, leukocyte migration and circulatory system (Fig. 4B). For MF, these co-DEGs were particularly enriched in extracellular matrix structural constituent, actin binding and calcium ion binding (Fig. 4C). KEGG pathway analysis revealed co-DEGs were primarily enriched in chemokine signaling pathway and regulation of actin cytoskeleton (Fig. 4D).
PPI network and hub gene identification
A PPI network with 105 nodes and 116 interacting pairs was constructed (Fig. 5A). The top10 DEGs were identified as hub genes by six rank methods of Cytoscape Cytohubba (Table 2). Among them, four genes (HCK, C1QC, TPM1 and CD14) were selected in all six methods, and FMOD and PLS3 with the lowest frequency appeared in the results of four methods. The identified 10 hub genes included 6 up-regulated genes (HCK C1QC, CD14, FCER1G, LCP1 and RAC2) and 4 down-regulated genes (TPM1, MYH10, FMOD and PLS3) (Fig. 5B). All of the hub genes were present in the top enriched functional terms and may play critical roles in atherosclerotic carotid unstable plaques. Furthermore, HCK and RAC2 were enriched in chemokine signaling pathway, MYH10 and RAC2 were enriched in regulation of actin cytoskeleton.
miRNA/TF‑hub gene regulatory network
Two hundred and twenty-four miRNAs targeting hub genes were predicted using miRNet. There were 736 differentially expressed miRNAs in the GSE111794 dataset, among which 12 miRNAs overlapped with the miRNet results (Fig. 6A), including 1 up-regulated and 11 down-regulated miRNAs targeting five hub genes (TPM1, MYH10, PLS3, LCP1 and FMOD). The Top5 TFs targeting hub genes were predicted to be GATA2, TP53, FOXC1, FOXL1 and JUN. Then a miRNA/TF-hub gene regulatory network was constructed (Fig. 6B). PLS3 was the most targeted hub genes. PLS3 was targeted by six miRNAs and all five screened TFs.
Drug-hub gene regulatory network
As shown in Fig. 7, 25 potential drugs were found for 4 up-regulated hub genes using the DGIdb database. Among them, HCK was targeted by 20 drugs including 10 inhibitors. Drugs targeted the other 6 hub genes were not predicted.
Discussion
We screened 105 co-DEGs from three data sets, including 42 up-regulated and 63 down-regulated genes. Through enrichment analysis of co-DEGs, we identified biological processes, cellular compartments and molecular functions associated with unstable carotid plaques. Atherosclerosis is associated with uncontrolled extracellular matrix (ECM) remodeling [23]. Reduced collagen content in ECM is a typical feature of unstable plaques [24]. Apolipoprotein E-deficient mice expressing collagenase-resistant collagen-I or lacking matrix metalloproteinase-13 (MMP-13/collagenase-3) can obtain collagen accumulation resulting in more stable plaques [25, 26]. Focal adhesions (FAs) are key junctions between cells and ECM and points of termination of actin filaments. FAs play an important role in tissue remodeling, integrity and homeostasis. The changes of FA constituents, such as the low expression of Talin and Vinculin, can affect the tissue remodeling and healing capabilities, which promotes the development of vulnerable plaques [27]. Actin filament bundles are the assemblies of actin filaments that participate in the regulation of endothelial cell (EC) adhesion to adjacent cells and matrix. During the development of arteriosclerosis, actin bundles change dynamically, and eventually the decrease or even disappearance of central actin microfilaments leads to dysfunction of cell–matrix adhesion [28]. In addition, calcification has been found to affect plaque stability. Common in unstable plaques are microcalcifications that originate from matrix vesicles rich in calcium-binding protein [29, 30]. The infiltration and activation of leukocytes within lesions also contribute to plaque rupture. For example, invariant natural killer T (iNKT) cells can promote rupture by activating inflammatory cells and up-regulating MMP-2 in vascular tissue [31]. In our study, the top2 enriched pathways are chemokine signaling pathway and regulation of actin cytoskeleton, each of which has been mentioned in previous studies to be involved in the progression of atherosclerotic plaque [32,33,34].
After constructing the PPI network, we identified 10 hub genes of unstable carotid plaque, including 6 up-regulated genes (HCK, C1QC, CD14, FCER1G, LCP1 and RAC2) and 4 down-regulated genes (TPM1, MYH10, PLS3 and FMOD). MiRNAs, TFs and drugs targeting these hub genes were also predicted.
HCK, a member of the SRC family of cytoplasmic tyrosine kinases (SFKs), is confined to the hematopoietic system, such as cells of myeloid and B lymphocyte lineages [35]. HCK promotes monocyte/macrophage intravasation by participating in a broad spectrum of processes including monocyte/macrophage proliferation, migration and endothelial adhesion, which is an essential mechanism of atherosclerosis involving a series of signaling pathways initiated by integrin, immune and growth factors, Fcy and chemokine receptors [36]. In addition, HCK can facilitate ECM degradation by phosphorylating the Wiskott-Aldrich syndrome protein (WASP) [37], which may contribute to plaque instability. However, Medina et al. indicated that atherosclerotic plaques showed reduced size in HCK/FGR double knockout mice but presented vulnerable phenotype characterized by necrotic core expansion, and significant reductions in collagen and fibrous cap thickness [36]. This may be due to a more complex response caused by the complete absence of HCK or the simultaneous deletion of FGR which is another member of SFKs, but it deserves further study. RAC2, a member of the Rho GTPases, is only found in cells of myeloid origin and is a crucial factor in leukocyte chemotaxis, as well as involved in the regulation of actin and microtubule cytoskeletal dynamics and adhesion [38]. RAC2 has been demonstrated to participate in reactive oxygen species (ROS) production [39], which induces the release of MMPs leading to plaque vulnerability by degrading the fibrous wall of atheromatous plaques and the basal membrane of endothelial cells [40]. MYH10 encodes a constituent protein of non-muscle myosin II, which is an actin-dependent motor protein and plays fundamental roles in cell adhesion and migration [41]. Kim et al. suggested that structural changes in the actomyosin network and defective ECM remodeling due to MYH10 deficiency contribute to the pathogenesis of emphysema [42], which may also contribute to unstable atherosclerotic plaques.
The other 7 hub genes were not enriched in the top 2 KEGG pathways. C1QC encodes the C-chain polypeptide of serum complement subcomponent C1q, a deficiency of which is associated with lupus erythematosus and glomerulonephritis. Lubbers et al. demonstrated that C1q can induce changes in ECM collagen expression [43]. C1QC was thought to regulate immune-competent cells involved in the progression of atherosclerosis [44]. CD14 is a differentiated antigen preferentially expressed on monocytes/macrophages and associated with inflammation and immune response. Inflammation can cause plaque disruption through processes such as endothelial cell death and activation of MMPs [45], but the role of CD14 needs to be further elucidated. It was found that the increase of CD14 + monocytes in coronary atherosclerosis patients was significantly correlated with the severity [46]. FCER1G is involved in encoding high affinity immunoglobulin epsilon receptors and has repeatedly been found to be overexpressed in the progression of atherosclerosis and unstable plaques. Although LCP1 (also named PLS2) and PLS3 are both actin binding proteins, LCP1 is overexpressed in unstable plaques, while PLS3 is underexpressed. LCP1 is essential for the degradation of ECM by macrophages. When LCP1 was inhibited by nanobodies, actin turnover was hampered and matrix degradation was significantly decreased [47]. PLS3 deficiency has been reported to cause osteoporosis and neurodegeneration [48]. TPM1 is a member of the tropomyosin family. Simoneau et al. proposed that TPM1 is required to maintain the endothelial barrier integrity and demonstrated that phosphorylation at Ser283 of TPM1 protects against oxidative stress-related endothelial barrier dysfunction [49]. FMOD encodes fibromodulin. In arteriosclerosis lesions, FMOD has stimulatory or stabilizing effects on collagen, smooth muscle cell proliferation, plaque lipids, inflammatory and proinflammatory cytokines [50].
In present study, we predicted 12 miRNAs that might be associated with unstable carotid plaques, several of which have been conducted the related studies. For example, mir-1910-5p was found to be highly expressed under oxidative stress [51], so it may be involved in ROS induced plaque destruction. Ye et al. demonstrated that lncRNA myocardial infarction associated transcript (MIAT) regulates the size of atherosclerotic necrotic core through the sponging mir-149-5p, while large necrotic core is a marker of unstable plaques [52]. A circulating microRNAs study observed that plasma mir-30e-5p was positively correlated with the volume of necrotic core in coronary plaques [53], but we found that mir-30e-5p was low expressed in unstable plaque samples. Whether this conflict is caused by different sample types or other reasons still needs further exploration. In addition, another study suggested that circulating mir-330-3p can be used to distinguish the plaque phenotype in patients with ST-segment elevation myocardial infarction [54]. TFs can affect gene expression and may play a role in unstable plaques. Therefore, we predicted the top five suspected TFs to improve the regulatory network of unstable plaques.
Several related drugs were predicted, including ten and one inhibitors of HCK and CD14. These drugs may be able to prevent or delay the progression of atherosclerotic plaque rupture. The eleven predicted inhibitors are mainly kinase inhibitors, which have been widely used in the treatment of malignant tumors. Several of them, such as nintedanib and bosutinib, have also been shown to ameliorate the development of atherosclerosis [55, 56].
Some limitations exist in our study. First of all, our analysis was based on four datasets from public data, whose grouping criteria may differ slightly, even though they all meet the histological criteria. In addition, only one miRNA dataset was used without cross-validation of other datasets, which may cause bias. Finally, we predicted the relevant molecules only by bioinformatic analysis, the true and complete regulatory mechanisms need to be verified by further in-vivo and in-vitro experiments.
Conclusions
In conclusion, the rupture of carotid atherosclerotic plaque involves extremely complex molecular mechanisms and regulatory networks. Bioinformatics analysis was used to identify 105 DEGs that are mainly enriched in the chemokine signaling pathway and the regulation of actin cytoskeleton pathway to participate in the formation of histologically unstable carotid atherosclerotic plaques. We screened 10 hub genes and predicted miRNAs and TFs targeting them to construct the regulatory network. These molecules may play a crucial role in the progression of unstable carotid plaques and serve as potential biomarkers and therapeutic targets.
Availability of data and materials
Publicly available datasets were analyzed in this study. This data can be found here: https://www.ncbi.nlm.nih.gov/geo/; https://www.ebi.ac.uk/arrayexpress/.
Abbreviations
- DEG:
-
Differentially expressed gene
- PPI:
-
Protein–protein interaction
- TF:
-
Transcription factor
- GEO:
-
Gene Expression Omnibus
- FPKM:
-
Fragments per kilobase million
- FC:
-
Fold change
- GO:
-
Gene ontology
- KEGG:
-
Kyoto Encyclopedia of Genes and Genomes
- CC:
-
Cellular compartments
- BP:
-
Biological processes
- MF:
-
Molecular function
- STRING:
-
The Search Tool for the Retrieval of Interacting Genes
- MCC:
-
Maximal clique centrality
- MNC:
-
Maximum neighborhood component
- EPC:
-
Edge percolated component
- DGIdb:
-
The Drug Gene Interaction Database
- ECM:
-
Extracellular matrix
- MMP:
-
Matrix metalloproteinase
- FA:
-
Focal adhesion
- EC:
-
Endothelial cell
- SFKs:
-
The SRC family of cytoplasmic tyrosine kinases
- WASP:
-
Wiskott–Aldrich syndrome protein
- ROS:
-
Reactive oxygen species
- MIAT:
-
Myocardial infarction associated transcript
References
Boehme AK, Esenwa C, Elkind MSV. Stroke risk factors, genetics, and prevention. Circ Res. 2017;120:472–95.
Liapis CD, Kakisis JD, Kostakis AG. Carotid stenosis: factors affecting symptomatology. Stroke. 2001;32:2782–6.
Badimon L, Vilahur G. Thrombosis formation on atherosclerotic lesions and plaque rupture. J Intern Med. 2014;276:618–32.
Stary HC, Chandler AB, Dinsmore RE, Fuster V, Glagov S, Insull W, et al. A definition of advanced types of atherosclerotic lesions and a histological classification of atherosclerosis. A report from the Committee on Vascular Lesions of the Council on Arteriosclerosis, American Heart Association. Circulation. 1995;92:1355–74.
Chen L, Yang Q, Ding R, Liu D, Chen Z. Carotid thickness and atherosclerotic plaque stability, serum inflammation, serum MMP-2 and MMP-9 were associated with acute cerebral infarction. Exp Ther Med. 2018;16:5253–7.
Salem MK, Butt HZ, Choke E, Moore D, West K, Robinson TG, et al. Gene and protein expression of chemokine (C-C-Motif) ligand 19 is upregulated in unstable carotid atherosclerotic plaques. Eur J Vasc Endovasc Surg. 2016;52:427–36.
Zhou G, Soufan O, Ewald J, Hancock REW, Basu N, Xia J. NetworkAnalyst 3.0: a visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucl Acids Res. 2019;47:W234–41.
Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10:1523.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The gene ontology consortium. Nat Genet. 2000;25:25–9.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucl Acids Res. 2000;28:27–30.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucl Acids Res. 2019;47:D607–13.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
Chin C, Chen S, Wu H, Ho C, Ko M, Lin C. cytoHubba: identifying hub objects and sub-networks from complex interactome. Bmc Syst Biol. 2014;8(Suppl 4):S11.
Chang Le, Zhou G, Soufan O, Xia J. miRNet 2.0: network-based visual analytics for miRNA functional analysis and systems biology. Nucl Acids Res. 2020;48:W244–51.
Huang H, Lin Y, Li J, Huang K, Shrestha S, Hong H, et al. miRTarBase 2020: updates to the experimentally validated microRNA-target interaction database. Nucl Acids Res. 2020;48:D148–54.
Karagkouni D, Paraskevopoulou MD, Chatzopoulos S, Vlachos IS, Tastsoglou S, Kanellos I, et al. DIANA-TarBase v8: a decade-long collection of experimentally supported miRNA-gene interactions. Nucl Acids Res. 2018;46:D239–45.
Xiao F, Zuo Z, Cai G, Kang S, Gao X, Li T. miRecords: an integrated resource for microRNA-target interactions. Nucl Acids Res. 2009;37:D105–10.
Fornes O, Castro-Mondragon JA, Khan A, van der Lee R, Zhang X, Richmond PA, et al. JASPAR 2020: update of the open-access database of transcription factor binding profiles. Nucl Acids Res. 2020;48:D87-92.
Freshour SL, Kiwala S, Cotto KC, Coffman AC, McMichael JF, Song JJ, et al. Integration of the drug-gene interaction database (DGIdb 4.0) with open crowdsource efforts. Nucl Acids Res. 2021;49:D1144–51.
Wishart DS, Feunang YD, Guo AC, Lo EJ, Marcu A, Grant JR, et al. DrugBank 5.0: a major update to the DrugBank database for 2018. Nucl Acids Res. 2018;46:D1074–82.
Whirl-Carrillo M, Huddart R, Gong L, Sangkuhl K, Thorn CF, Whaley R, et al. An evidence-based framework for evaluating pharmacogenomics knowledge for personalized medicine. Clin Pharmacol Ther. 2021;110:563–72.
Tang J, Tanoli Z, Ravikumar B, Alam Z, Rebane A, Vähä-Koskela M, et al. Drug target commons: a community effort to build a consensus knowledge base for drug-target interactions. Cell Chem Biol. 2018;25:224–9.
Buchler A, Munch M, Farber G, Zhao X, Al-Haddad R, Farber E, et al. Selective imaging of matrix metalloproteinase-13 to detect extracellular matrix remodeling in atherosclerotic lesions. Mol Imag Biol. 2022;24:93–103.
Newby AC. Metalloproteinases and vulnerable atherosclerotic plaques. Trends Cardiovas Med. 2007;17:253–8.
Fukumoto Y, Deguchi J, Libby P, Rabkin-Aikawa E, Sakata Y, Chin MT, et al. Genetically determined resistance to collagenase action augments interstitial collagen accumulation in atherosclerotic plaques. Circulation. 2004;110:1953–9.
Deguchi J, Aikawa E, Libby P, Vachon JR, Inada M, Krane SM, et al. Matrix metalloproteinase-13/collagenase-3 deletion promotes collagen accumulation and organization in mouse atherosclerotic plaques. Circulation. 2005;112:2708–15.
von Essen M, Rahikainen R, Oksala N, Raitoharju E, Seppälä I, Mennander A, et al. Talin and vinculin are downregulated in atherosclerotic plaque. Tampere Vas Study Atherosclerosis. 2016;255:43–53.
Colangelo S, Langille BL, Steiner G, Gotlieb AI. Alterations in endothelial F-actin microfilaments in rabbit aorta in hypercholesterolemia. Arterioscler Thromb Vasc Biol. 1998;18:52–6.
Jinnouchi H, Sato Y, Sakamoto A, Cornelissen A, Mori M, Kawakami R, et al. Calcium deposition within coronary atherosclerotic lesion: Implications for plaque stability. Atherosclerosis. 2020;306:85–95.
Hutcheson JD, Maldonado N, Aikawa E. Small entities with large impact: microcalcifications and atherosclerotic plaque vulnerability. Curr Opin Lipidol. 2014;25:327–32.
Ohmura Y, Ishimori N, Saito A, Yokota T, Horii S, Tokuhara S, et al. Natural killer T cells are involved in atherosclerotic plaque instability in apolipoprotein-E knockout mice. Int J Mol Sci. 2021;22:12451.
Tan X, Zhang X, Pan L, Tian X, Dong P. Identification of key pathways and genes in advanced coronary atherosclerosis using bioinformatics analysis. Biomed Res Int. 2017;2017:4323496.
Qian C, Xia M, Yang X, Chen P, Ye Q. Long noncoding RNAs in the progression of atherosclerosis: an integrated analysis based on competing endogenous RNA theory. Dna Cell Biol. 2021;40:283–92.
Yang R, Yao L, Du C, Wu Y. Identification of key pathways and core genes involved in atherosclerotic plaque progression. Ann Transl Med. 2021;9:267.
Poh AR, O’Donoghue RJJ, Ernst M. Hematopoietic cell kinase (HCK) as a therapeutic target in immune and cancer cells. Oncotarget. 2015;6:15752–71.
Medina I, Cougoule C, Drechsler M, Bermudez B, Koenen RR, Sluimer J, et al. Hck/Fgr kinase deficiency reduces plaque growth and stability by blunting monocyte recruitment and intraplaque motility. Circulation. 2015;132:490–501.
Park H, Dovas A, Hanna S, Lastrucci C, Cougoule C, Guiet R, et al. Tyrosine phosphorylation of Wiskott-Aldrich syndrome protein (WASP) by Hck regulates macrophage function. J Biol Chem. 2014;289:7897–906.
Bokoch GM. Regulation of innate immunity by Rho GTPases. Trends Cell Biol. 2005;15:163–71.
Diebold BA, Bokoch GM. Rho GTPases and the control of the oxidative burst in polymorphonuclear leukocytes. Curr Top Microbiol. 2005;291:91–111.
Kattoor AJ, Pothineni NVK, Palagiri D, Mehta JL. Oxidative stress in atherosclerosis. Curr Atheroscler Rep. 2017;19:42.
Heissler SM, Manstein DJ. Nonmuscle myosin-2: mix and match. Cell Mol Life Sci. 2013;70:1–21.
Kim H, Yin W, Jin Y, Panza P, Gunawan F, Grohmann B, et al. Myh10 deficiency leads to defective extracellular matrix remodeling and pulmonary disease. Nat Commun. 2018;9:4600.
Lubbers R, van Schaarenburg RA, Kwekkeboom JC, Levarht EWN, Bakker AM, Mahdad R, et al. Complement component C1q is produced by isolated articular chondrocytes. Osteoarthr Cartilage. 2020;28:675–84.
Zhao B, Wang D, Liu Y, Zhang X, Wan Z, Wang J, et al. Six-gene signature associated with immune cells in the progression of atherosclerosis discovered by comprehensive bioinformatics analyses. Cardiovasc Ther. 2020;2020:1230513.
Libby P. Inflammation in atherosclerosis. Nature. 2002;420:868–74.
Lee M, Huang C, Kuo C, Lin CS, Lai W, Chiou S. Clinical proteomics identifies urinary CD14 as a potential biomarker for diagnosis of stable coronary artery disease. PLoS ONE. 2015;10:e117169.
De Clercq S, Boucherie C, Vandekerckhove J, Gettemans J, Guillabert A. L-plastin nanobodies perturb matrix degradation, podosome formation, stability and lifetime in THP-1 macrophages. PLoS ONE. 2013;8: e78108.
Wolff L, Strathmann EA, Müller I, Mählich D, Veltman C, Niehoff A, et al. Plastin 3 in health and disease: a matter of balance. Cell Mol Life Sci. 2021;78:5275–301.
Simoneau B, Houle F, Huot J. Regulation of endothelial permeability and transendothelial migration of cancer cells by tropomyosin-1 phosphorylation. Vasc cell. 2012;4:18.
Hultgårdh-Nilsson A, Borén J, Chakravarti S. The small leucine-rich repeat proteoglycans in tissue repair and atherosclerosis. J Intern Med. 2015;278:447–61.
Ayaz L, Dinç E. Evaluation of microRNA responses in ARPE-19 cells against the oxidative stress. Cutan Ocul Toxicol. 2018;37:121–6.
Ye Z, Yang S, Xia Y, Hu R, Chen S, Li B, et al. LncRNA MIAT sponges miR-149-5p to inhibit efferocytosis in advanced atherosclerosis through CD47 upregulation. Cell Death Dis. 2019;10:138.
Taraldsen MD, Wiseth R, Videm V, Bye A, Madssen E. Associations between circulating microRNAs and coronary plaque characteristics: potential impact from physical exercise. Physiol Genom. 2022;54:129–40.
Li J, Tan Y, Sheng Z, Zhou P, Liu C, Zhao H, et al. RNA-seq identifies circulating miRNAs as potential biomarkers for plaque rupture in patients with ST-segment elevation myocardial infarction. Genomics. 2021;113:1–10.
Li L, Chen Y, Shi C. Nintedanib ameliorates oxidized low-density lipoprotein -induced inflammation and cellular senescence in vascular endothelial cells. Bioengineered. 2022;13:6196–207.
Li B, He J, Lv H, Liu Y, Lv X, Zhang C, et al. c-Abl regulates YAPY357 phosphorylation to activate endothelial atherogenic responses to disturbed flow. J Clin Invest. 2019;129:1167–79.
Acknowledgements
We would like to thank all the research staff who made it possible to perform this study.
Funding
This study was supported by the National Key R&D Program of China (No.2021YFC2500500).
Author information
Authors and Affiliations
Contributions
YG and LG contributed to conception and design of the study. JG, YN and ZS organized the database and performed the statistical analysis. JG wrote the manuscript and prepared figures. All authors have read and agreed to the published version of the manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Table S1
. Top 20 terms (one per cluster) of functional enrichment analysis.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Guo, J., Ning, Y., Su, Z. et al. Identification of hub genes and regulatory networks in histologically unstable carotid atherosclerotic plaque by bioinformatics analysis. BMC Med Genomics 15, 145 (2022). https://doi.org/10.1186/s12920-022-01257-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12920-022-01257-1