Nonnegative matrix factorization analysis and multiple machine learning methods identified IL17C and ACOXL as novel diagnostic biomarkers for atherosclerosis
BMC Bioinformatics volume 24, Article number: 196 (2023)
Atherosclerosis is the common pathological basis for many cardiovascular and cerebrovascular diseases. The purpose of this study is to identify the diagnostic biomarkers related to atherosclerosis through machine learning algorithm.
Clinicopathological parameters and transcriptomics data were obtained from 4 datasets (GSE21545, GSE20129, GSE43292, GSE100927). A nonnegative matrix factorization algorithm was used to classify arteriosclerosis patients in GSE21545 dataset. Then, we identified prognosis-related differentially expressed genes (DEGs) between the subtypes. Multiple machine learning methods to detect pivotal markers. Discrimination, calibration and clinical usefulness of the predicting model were assessed using area under curve, calibration plot and decision curve analysis respectively. The expression level of the feature genes was validated in GSE20129, GSE43292, GSE100927.
2 molecular subtypes of atherosclerosis was identified, and 223 prognosis-related DEGs between the 2 subtypes were identified. These genes are not only related to epithelial cell proliferation, mitochondrial dysfunction, but also to immune related pathways. Least absolute shrinkage and selection operator, random forest, support vector machine- recursive feature elimination show that IL17C and ACOXL were identified as diagnostic markers of atherosclerosis. The prediction model displayed good discrimination and good calibration. Decision curve analysis showed that this model was clinically useful. Moreover, IL17C and ACOXL were verified in other 3 GEO datasets, and also have good predictive performance.
IL17C and ACOXL were diagnostic genes of atherosclerosis and associated with higher incidence of ischemic events.
Atherosclerosis (AS) is a chronic progressive inflammatory disease of blood vessels, which involves physiological processes such as disorders of lipid metabolism, vascular endothelial cell injury, inflammatory cell infiltration, and neo-capillary formation [1,2,3]. AS has no obvious symptoms in the early stage of the disease, and patients are mostly aware of it because of other cardiovascular and cerebrovascular complications . Although drugs are available to treat atherosclerosis, there are still many patients who do not benefit from current drug therapy without significant effects or who are intolerant to adverse effects . The search for highly sensitive and specific biomarkers can help reduce the morbidity and mortality of AS. However, current timely diagnostic biomarkers for AS are very limited and not well suited for the early and accurate diagnosis of AS. Therefore, it is important to find new diagnostic markers of AS for accurate diagnosis of AS.
Currently, the development of microarrays has led to extensive and in-depth analysis of genome-wide mRNA expression profiles. With the rapid development of gene chips, high-throughput sequencing, multi-omics analysis and other technologies, gene expression public databases are rapidly increasing. With the development and maturation of bioinformatics, bioinformatics techniques are widely used to analyze large number of expression profiling microarrays to find biomarkers related to disease diagnosis, treatment and prognosis [6,7,8]. Machine learning has been widely used in finding markers for disease diagnosis base on multi-omics analysis. Support vector machine (SVM), least absolute shrinkage and selection operator (LASSO) regression and random forest (RF) methods are 3 important techniques in machine learning [9,10,11]. Due to the three methods can identify the best classification feature factor and build a prediction model with generalizability and high prediction accuracy. Xiong et al. screened 2 mRNAs as potential diagnostic biomarkers for abdominal aortic aneurysm using machine learning .
In this study, we used multiple machine learning methods analysis to analyze the expression profile microarrays of AS, aiming to screen genes closely related to AS diagnosis and provide new genetic diagnostic markers for AS. We firstly identified 2 subtypes of atherosclerosis by using nonnegative matrix factorization (NMF) algorithm in GSE21545 dataset, and identified prognosis-related DEGs between the subtypes. After that, we identified 2 potential diagnosis genes by using LASSO, RF, SVM-RFE methods and developed a novel prediction model for AS diagnosis. We validated the model and found that the novel prediction model achieved a high AUC in 3 validation AS cohort.
Materials and methods
Data collection and pre-processing
The gene expression profiles of human atherosclerosis samples and healthy control samples were obtained from array-based data available in the Gene Expression Omnibus (GEO) database. The search strategy aimed to find published dataset which included a three-step search strategy that was carried out from inception to July, 2022. An initial limited search using the keywords: “Atherosclerosis”, “Atheromatosis”, “Homo sapiens”, “Expression profiling by array”. Dataset needs to meet the following points: (1) Homo sapiens; (2) Atherosclerosis; (3) Sample size greater than 20 cases. Four microarray data sets (GSE21545 , GSE20129 , GSE43292 , GSE100927 ) were utilized in our analysis (Additional file 1: Fig. S1). From the GSE21545 dataset, 126 AS samples were included. From GSE20129, 71 control samples and 48 AS samples were included. From GSE43292, the 32 control samples and 32 AS samples were included. From GSE100927, 35 control samples and 69 AS samples were included. The training set was obtained from GSE21545, and the validation set was obtained from GSE20129, GSE43292, GSE100927. The raw files from the four datasets were pre-processed and normalized using limma or RMA-affylmGUI in R Bioconductor.
Nonnegative matrix factorization (NMF) analysis in GSE21545
The R package “NMF” was performed to identify molecular subtypes based on the gene expression profiles, and patients were classified for follow up studies. We used a NMF algorithm to determine the number of clusters and their stability according to parameters such as cophenetic, dispersion, silhouette, and sparseness .
Identification of differentially expressed genes (DEGs) and prognosis genes between subtypes
DEGs between subtypes were identified using the R package limma with screening criteria of adjusted P value < 0.05 . Then, we performed univariate COX analysis to determine the prognostic value of each DEGs.
Functional and pathway enrichment analysis
To explore the biological functions mainly performed by DEGs and prognosis genes between subtypes, we performed functional enrichment analysis using the "clusterProfiler" package, including GO and KEGG analysis. The screening criteria were P < 0.05 and FDR < 0.05 [19, 20].
Evaluation of immune infiltrating cells in AS
Based on the normalized gene expression data from the disease and control samples, the web tool CIBERSORT (http://CIBERSORT.stanford.edu/) was used to calculate immune cell infiltration and explore the disease immune microenvironment. The 22 immune cell genes (LM22) were used as the reference set. The number of permutations set was 1000. A P value < 0.05 in the CIBERSORT results was retained .
Machine learning methods
LASSO method, which is suitable for the reduction in high-dimensional data, was used to select the optimal predictive features in risk factors from the patients with AS. Support vector machine-recursive feature elimination (SVM-RFE) approach is based on the VC dimensional theory of statistical learning theory and the structural risk minimization principle. Based on limited sample information, SVM-RFE seeks to find the best compromise between the complexity of the model (the learning accuracy) and the learning ability. Random forest (RF) refers to a classifier that uses multiple trees to train and predict samples. The three classifiable models’ overlapping genes were then figured out.
Construction of AS diagnosis nomogram
We used the expression level of predictors by the R-package “rms” to construct the nomogram and predict the risk of AS. Calibration curves were used to estimate the consistency between predicted and actual diagnosis, and the performance of the model in predicting diagnosis was evaluated by area under curve (AUC) .
R software (4.1.2) was employed to carry out all statistical analysis and graph plotting. Wilcoxon test was applied to analyze the differences between two groups. Kruskal–Wallis test was used for comparison among more than two groups of samples. The Kaplan–Meier method was used to plot survival curves for prognostic analysis, and the log-rank test was used to determine the significance of differences. The correlation test was performed using Spearman correlation analysis and distance correlation analysis. Comparisons of composition ratios among groups were performed by chi-square test. All statistical P values were two-tailed, and P < 0.05 was used as the truncated value.
Identification of molecular subtypes in AS
To explore the expression characteristics in AS, we qualitatively classified patients based on the expression profiles. By NMF algorithm, a cluster number of 2 was the best choice to classify the whole sample into C1 (n = 46) and C2 (n = 80) in GSE21545 dataset (Additional file 2: Fig. S2, Fig. 1A). Kaplan–Meier survival analysis indicated that patients with RPMRs.cluster.A had a worse ischemic events ((HR 4.08, 95% CI 1.22–13.63, P = 0.023, Fig. 1B). To explore the potential biological change between distinct cluster, firstly, the PCA demonstrated there is significant DEGs between the two clusters, and 223 prognostic related DEGs were identified (Fig. 1C, Additional file 3: Table S1). Then, we applied GO and KEGG enrichment analysis, which showed that C2 was significantly enriched in immune-related pathways (mast cell activation, regulation of interferon-gamma production, and IL-17 signaling pathway), epithelial cell proliferation (epithelial cell proliferation, regulation of epithelial cell proliferation), mitochondrial dysfunction (mitochondrial inner membrane), suggesting that C2 may play an important role in AS development and immune regulation (Fig. 1D, E).
Immune infiltrating cell analysis between the two molecular subtypes
To further investigate the role of two cluster in immune infiltrating of AS, we used CIBERSORT to explore the infiltration of various types of immune cells in AS samples. Figure 2A indicated the immune cell infiltration landscape and immune cell score correlation results in different samples of the GSE21545 dataset, respectively. Moreover, univariate Cox regression analysis was performed based on GSE21545 dataset, high immune cell score of eosinophils, mast cells activated, B cells memory corelated with worse ischemic events in AS (Fig. 2B). In addition, we also evaluate the association between molecular subtypes and immune cell subpopulations. The results showed T cells CD8, T cells regulatory (Tregs), M0 macrophages, and M1 macrophages were more abundant in C2 cluster, while M2 macrophages was significantly higher in C1 cluster (Fig. 2C). The above results suggest that C2 cluster had a higher inflammatory environment, which leads to the progression of the disease.
Machine learning methods to detect diagnostic markers in AS
To further explore the risk gene features, we carried out the LASSO regression model to screen out 33 potential predictors from 223 prognostic related DEGs (Fig. 3A, B, Additional file 4: Table S2). To evaluate the discrimination of the prediction model, the AUC of ROC was estimated. As shown in Fig. 3C, the prediction model achieved a AUC of 0.930 (95% CI 0.887–0.973), which indicated good discrimination of the model. Then, we carried out the SVM-RFE approach (k = 10, halve.above = 50) to screen out 17 potential predictors from 223 prognostic related DEGs (Fig. 3D, E, Additional file 5: Table S3). To evaluate the discrimination of the prediction model, the AUC of ROC was estimated. As shown in Fig. 3F, the prediction model achieved a AUC of 0.981 (95% CI 0.964–0.998), which indicated good discrimination of the model. In addition, we carried out the RF approach (ntree = 500) to screen out 3 potential predictors from 223 prognostic related DEGs (Fig. 3G, H, Additional file 6: Table S4). To evaluate the discrimination of the prediction model, the AUC of ROC was estimated. As shown in Fig. 3I, the prediction model achieved a AUC of 0.997 (95% CI 0.992–1.000), which indicated good discrimination of the model. The above results show that the three machine learning methods have good performance in identifying diagnostic markers.
IL17C and ACOXL were identified as diagnostic biomarkers for AS
To further explore the diagnostic biomarkers for AS, the three classifiable models’ overlapping genes were figured out (Fig. 4A). These genes included IL17C and ACOXL. Next, we compared the expression of IL17C and ACOXL in two molecular subtypes, of which the expression level of IL17C and ACOXL were significantly higher in C2 cluster than in C1 cluster (Fig. 4B, E). To evaluate the discrimination of the diagnostic biomarkers, the AUC of ROC was estimated. The IL17C achieved a AUC of 0.917 (95% CI 0.865–0.968) (Fig. 4C) and the ACOXL achieved a AUC of 0.899 (95% CI 0.842–0.955) (Fig. 4F), which indicated good discrimination of the diagnostic biomarkers. Next, Kaplan–Meier survival analysis was performed based on GSE21545 dataset, high expression of IL17C coorelated with worse ischemic events ((HR 2.43, 95% CI 1.05–5.64, P = 0.039, Fig. 4D), high expression of ACOXL coorelated with worse ischemic events (HR 2.68, 95% CI 1.12–6.43, P = 0.027, Fig. 4G). Then, we applied KEGG enrichment analysis, which showed that high expression IL17C was significantly enriched in immune-related pathways (B cell receptor signaling pathway, T cell receptor signaling pathway), metabolic pathways (citrate cycle TCA cycle, gylcosylphosphatidylinositol GPI anchor biosynthsis, sphingolipid metabolism), high expression ACOXL was significantly enriched in metabolic pathways (linoleic acid metabolism, alpha linolenic acid metabolism) (Fig. 4H, I).
Relationship between diagnostic biomarkers and immune cells
Next, we investigate the role of diagnostic biomarkers in immune infiltrating of AS (Fig. 5A, G). The results showed the expression of ACOXL was positively correlated with macrophages M1 (r = 0.27, P = 0.002), monocytes (r = 0.20, P = 0.027), T cells CD8 (r = 0.40, P < 0.001), T cells regulatory (Tregs) (r = 0.22, P = 0.012), and negatively correlated with T cells CD4 naive (r = -0.25, P = 0.004) (Fig. 5B–F). The expression of IL17C was positively correlated with dendritic cells resting (r = 0.19, P = 0.031), macrophages M1 (r = 0.24, P = 0.008), T cells CD8 (r = 0.57, P < 0.001), and T cells regulatory (Tregs) (r = 0.31, P < 0.001) (Fig. 5H–K).
Construction of AS diagnosis nomogram
Next, we developed a diagnosis model of AS. The model that incorporated the above independent predictors was developed and presented as the nomogram (Fig. 6A). The calibration curve of the AS diagnosis nomogram for the prediction of AS risk demonstrated good agreement in this cohort (Fig. 6B). The decision curve analysis for the nomogram was presented in Fig. 6C, D. The decision curve showed that if the threshold probability of a patient and a doctor is > 1 and < 96%, respectively, using this nomogram to predict AS risk adds more benefit than the scheme.
IL17C and ACOXL were a robust diagnostic biomarkers for AS in GSE20129, GSE43292, and GSE100927 datasets
Consistent with this finding, increased mRNA expression of IL17C and ACOXL were observed in atherosclerosis compared with that in normal in GSE20129 datasets (Fig. 7A, C). The IL17C achieved a AUC of 0.892 (95% CI 0.833–0.950) (Fig. 7B) and the ACOXL achieved a AUC of 0.750 (95% CI 0.661–0.838) (Fig. 7D), which indicated good discrimination of the diagnostic biomarkers. Moreover, we compared the expression of IL17C and ACOXL in GSE43292 dataset, of which the expression level of IL17C and ACOXL were significantly higher in AS than in the normal tissues (Fig. 7E, G). ROC curves with AUC values was 0.852 (95% CI 0.759–0.944) for IL17C, 0.925 (95% CI 0.865–0.984) for ACOXL (Fig. 7F, H). In addition, we compared the expression of IL17C and ACOXL in GSE100927 dataset, of which the expression level of IL17C and ACOXL were significantly higher in AS than in the normal tissues (Fig. 7I, K). ROC curves with AUC values was 0.818 (95% CI 0.736–0.899) for IL17C, 0.803 (95% CI 0.719–0.887) for ACOXL (Fig. 7J, L). The above results suggest that IL17C and ACOXL were a robust diagnostic biomarkers for AS.
In this study, we identified two molecular subtypes based on the expression profiles, with C2 cluster showing a worse ischemic events. In addition, mRNA transcriptome differential expression genes between distinct cluster were closely related to biological processes such as immune-related pathways, epithelial cell proliferation pathways, and mitochondrial dysfunction pathways. In addition, C2 cluster had a higher T cells CD8, T cells regulatory (Tregs), M0 macrophages, and M1 macrophages, and lower M2 macrophages.
The immune system is one of the major regulatory systems in the development and progression of atherosclerosis . In the early stages of atherosclerosis, low-density lipoprotein (LDL) is retained in the intima and is modified by oxidases, lipolytic enzymes, protein hydrolases and reactive oxygen species to form a variety of risk-related molecular patterns, thereby acquiring immunogenicity . Immunogenic LDL activates vascular endothelial cells, which chemotacticize various immune cells into the vessel wall, mainly monocytes and T cells . Histological analysis of human atherosclerotic plaques showed that M1 macrophages were mainly distributed in the lipid core, whereas M2 macrophages were mainly distributed in the plaque region away from the lipid core . In vitro mouse experiments showed that M1 macrophages promote plaque inflammation, while M2 macrophages promote plaque inflammation regression . We found C2 cluster had a higher M0 macrophages, and M1 macrophages, and loewer M2 macrophages, suggesting C2 cluster had a higher inflammatory environment, which leads to the progression of the disease.
CD4 + T cells receive antigens presented by antigen-presenting cells and differentiate into different Th cells (Th1, Th2, Th9, Th17, Th22, Tfh) and Treg cells through immune responses, whose role in atherosclerosis is multifaceted. Secretion of IFN-γ by Th1 cells affects macrophage polarization by inhibiting VSMC proliferation, thereby inhibiting plaque stability . In addition to IFN-γ, Th1 cells secrete IL-2, IL-3, tumor necrosis factor, and lymphotoxin, all of which activate macrophages, T cells, and other cells within the plaque, thereby accelerating the inflammatory response . At the same time, CD8+ T cells act on VSMC and release some inflammatory factors that make the atherosclerotic plaque unstable as well as aggravate the inflammatory response .
Wang et al. showed that CD68 (AUC = 0.80), PAM (AUC = 0.79), and IGFBP6 (AUC = 0.81) could be used as diagnostic markers to identify unstable plaques effectively by using LASSO and RF . Xu et al. showed that C1QA (AUC = 0.83) and ITGB2 (AUC = 0.83) could be used as diagnostic markers to identify unstable plaques effectively by using LASSO . In this study, multiple machine learning methods (LASSO, RF, SVM-RFE) identified IL17C (AUC = 0.92) and ACOXL (AUC = 0.90) as novel diagnostic biomarkers for atherosclerosis, and verified in other datasets.
Interleukin-17C (IL-17C) is one of the important members of the IL-17 cytokine family, which can be secreted by many types of cells or produced by the stimulation of pathogenic factors. IL-17C is mainly expressed in the mucosal surface of the gastrointestinal and respiratory tracts as well as the skin barrier. In the gastrointestinal tract, IL-17C is secreted by enteroendocrine cells and cupped cells. In the skin, it is mainly expressed by keratinocytes, monocytes and endothelial cells [33, 34]. In this present study, the expression level of IL17C was significantly higher in AS than in the normal tissues, which was consistent with previous findings. IL-17C exerts a proatherogenic effect by recruiting Th17 cells to atherosclerotic plaques .
Acyl coenzyme A oxidase like gene is a member of the acyl coenzyme A oxidase family. Paul et al. found that in mammals ACOXL is actively expressed at the transcriptome level, and that ACOXL is specifically expressed in the lung. ACOXL has a dehydrogenase activity of acyl coenzyme A and also catalyzes an important step in the β-oxidation pathway involving the oxidation of long-chain fatty acids . Gillian et al. identified ACOXL as a biomarker for the diagnosis of prostate cancer through transcriptomics and antibody analysis of the human prostate-specific proteome . In this study, we found ACOXL can be used a diagnostic biomarkers for AS, and metabolic pathways play an important role in AS disease progression. Metabolic intermediates or oxidation products produced during metabolism, such as oxidized LDL, ceramide, TMAO, and cholesterol crystals, can also be recognized by macrophages and cause activation of inflammatory pathways in the body, thus further aggravating the inflammatory response of the vasculature .
There are some limitations of our study. Although our analysis was based on a large sample, these cases were obtained retrospectively, and selection bias in the dataset may also affect the accuracy of the results. Large-scale prospective studies and in vivo, in vitro mechanistic studies are still needed to further confirm our results. In addition, some important clinical variables such as age, gender, and therapy information are missing in most of the datasets, we also need to combine more clinical characteristics to improve the prediction accuracy.
In conclusion, we identified IL17C and ACOXL were diagnostic genes of atherosclerosis and associated with higher incidence of ischemic events. These findings may provide a new strong scientific basis for the diagnosis and treatment of atherosclerotic.
Availability of data and materials
All data used in the study can be downloaded from the GEO database (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE21545; https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE20129;https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE43292;https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE100927; accessed 26 May 2022).
Bäck M, Yurdagul A Jr, Tabas I, Öörni K, Kovanen PT. Inflammation and its resolution in atherosclerosis: mediators and therapeutic opportunities. Nat Rev Cardiol. 2019;16(7):389–406.
Ahmadi A, Jamialahmadi T, Sahebkar A. Polyphenols and atherosclerosis: a critical review of clinical effects on LDL oxidation. Pharmacol Res. 2022;184:106414.
Ross R. Atherosclerosis–an inflammatory disease. N Engl J Med. 1999;340(2):115–26.
Gallino A, Aboyans V, Diehm C, Cosentino F, Stricker H, Falk E, Schouten O, Lekakis J, Amann-Vesti B, Siclari F, et al. Non-coronary atherosclerosis. Eur Heart J. 2014;35(17):1112–9.
Yang YW, Deng NH, Tian KJ, Liu LS, Wang Z, Wei DH, Liu HT, Jiang ZS. Development of hydrogen sulfide donors for anti-atherosclerosis therapeutics research: Challenges and future priorities. Front Cardiovasc Med. 2022;9:909178.
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.
Herrington DM, Mao C, Parker SJ, Fu Z, Yu G, Chen L, Venkatraman V, Fu Y, Wang Y, Howard TD, et al. Proteomic architecture of human coronary and aortic atherosclerosis. Circulation. 2018;137(25):2741–56.
Huang HM, Jiang X, Hao ML, Shan MJ, Qiu Y, Hu GF, Wang Q, Yu ZQ, Meng LB, Zou YY. Identification of biomarkers in macrophages of atherosclerosis by microarray analysis. Lipids Health Dis. 2019;18(1):107.
Wang H, Yang F, Luo Z. An experimental study of the intrinsic stability of random forest variable importance measures. BMC Bioinform. 2016;17:60.
Albaradei S, Thafar M, Alsaedi A, Van Neste C, Gojobori T, Essack M, Gao X. Machine learning and deep learning methods that use omics data for metastasis prediction. Comput Struct Biotechnol J. 2021;19:5008–18.
He Y, Ma J, Wang A, Wang W, Luo S, Liu Y, Ye X. A support vector machine and a random forest classifier indicates a 15-miRNA set related to osteosarcoma recurrence. Onco Targets Ther. 2018;11:253–69.
Xiong T, Lv XS, Wu GJ, Guo YX, Liu C, Hou FX, Wang JK, Fu YF, Liu FQ. Single-cell sequencing analysis and multiple machine learning methods identified G0S2 and HPSE as novel biomarkers for abdominal aortic aneurysm. Front Immunol. 2022;13:907309.
Folkersen L, Persson J, Ekstrand J, Agardh HE, Hansson GK, Gabrielsen A, Hedin U, Paulsson-Berne G. Prediction of ischemic events on the basis of transcriptomic and genomic profiling in patients undergoing carotid endarterectomy. Mol Med (Cambridge, Mass). 2012;18(1):669–75.
Huang CC, Lloyd-Jones DM, Guo X, Rajamannan NM, Lin S, Du P, Huang Q, Hou L, Liu K. Gene expression variation between African Americans and whites is associated with coronary artery calcification: the multiethnic study of atherosclerosis. Physiol Genomics. 2011;43(13):836–43.
Ayari H, Bricca G. Identification of two genes potentially associated in iron-heme homeostasis in human carotid plaque using microarray analysis. J Biosci. 2013;38(2):311–5.
Steenman M, Espitia O, Maurel B, Guyomarch B, Heymann MF, Pistorius MA, Ory B, Heymann D, Houlgatte R, Gouëffic Y, et al. Identification of genomic differences among peripheral arterial beds in atherosclerotic and healthy arteries. Sci Rep. 2018;8(1):3940.
Kim H, Park H. Sparse non-negative matrix factorizations via alternating non-negativity-constrained least squares for microarray data analysis. Bioinformatics (Oxford, England). 2007;23(12):1495–502.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
Altermann E, Klaenhammer TR. PathwayVoyager: pathway mapping using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. BMC Genomics. 2005;6:60.
Tweedie S, Ashburner M, Falls K, Leyland P, McQuilton P, Marygold S, Millburn G, Osumi-Sutherland D, Schroeder A, Seal R et al. FlyBase: enhancing Drosophila Gene Ontology annotations. Nucleic Acids Res 2009;37(Database issue):D555–559.
Newman AM, Steen CB, Liu CL, Gentles AJ, Chaudhuri AA, Scherer F, Khodadoust MS, Esfahani MS, Luca BA, Steiner D, et al. Determining cell type abundance and expression from bulk tissues with digital cytometry. Nat Biotechnol. 2019;37(7):773–82.
Park SY. Nomogram: An analogue tool to deliver digital knowledge. J Thorac Cardiovasc Surg. 2018;155(4):1793.
Wolf D, Ley K. Immunity and inflammation in atherosclerosis. Circ Res. 2019;124(2):315–27.
Gisterå A, Hansson GK. The immunology of atherosclerosis. Nat Rev Nephrol. 2017;13(6):368–80.
Abdolmaleki F, Gheibi Hayat SM, Bianconi V, Johnston TP, Sahebkar A. Atherosclerosis and immunity: a perspective. Trends Cardiovasc Med. 2019;29(6):363–71.
Chinetti-Gbaguidi G, Baron M, Bouhlel MA, Vanhoutte J, Copin C, Sebti Y, Derudas B, Mayi T, Bories G, Tailleux A, et al. Human atherosclerotic plaque alternative macrophages display low cholesterol handling but high phagocytosis because of distinct activities of the PPARγ and LXRα pathways. Circ Res. 2011;108(8):985–95.
Colin S, Chinetti-Gbaguidi G, Staels B. Macrophage phenotypes in atherosclerosis. Immunol Rev. 2014;262(1):153–66.
Orecchioni M, Ghosheh Y, Pramod AB, Ley K. Macrophage polarization: different gene signatures in M1(LPS+) vs. classically and M2(LPS-) vs. alternatively activated macrophages. Front Immunol 2019;10:1084.
Wolf D, Ley K. Immunity and inflammation in atherosclerosis. Herz. 2019;44(2):107–20.
Chan WL, Pejnovic N, Hamilton H, Liew TV, Popadic D, Poggi A, Khan SM. Atherosclerotic abdominal aortic aneurysm and the interaction between autologous human plaque-derived vascular smooth muscle cells, type 1 NKT, and helper T cells. Circ Res. 2005;96(6):675–83.
Wang J, Kang Z, Liu Y, Li Z, Liu Y, Liu J. Identification of immune cell infiltration and diagnostic biomarkers in unstable atherosclerotic plaques by integrated bioinformatics analysis and machine learning. Front Immunol. 2022;13:956078.
Xu J, Chen C, Yang Y. Identification and validation of candidate gene module along with immune cells infiltration patterns in atherosclerosis progression to plaque rupture via transcriptome analysis. Front Cardiovasc Med. 2022;9:894879.
Nies JF, Panzer U. IL-17C/IL-17RE: emergence of a unique axis in T(H)17 biology. Front Immunol. 2020;11:341.
McGeachy MJ, Cua DJ, Gaffen SL. The IL-17 family of cytokines in health and disease. Immunity. 2019;50(4):892–906.
Butcher MJ, Waseem TC, Galkina EV. Smooth muscle cell-derived interleukin-17C plays an atherogenic role via the recruitment of proinflammatory interleukin-17A+ T cells to the aorta. Arterioscler Thromb Vasc Biol. 2016;36(8):1496–506.
Hong MG, Karlsson R, Magnusson PK, Lewis MR, Isaacs W, Zheng LS, Xu J, Grönberg H, Ingelsson E, Pawitan Y, et al. A genome-wide assessment of variability in human serum metabolism. Hum Mutat. 2013;34(3):515–24.
O’Hurley G, Busch C, Fagerberg L, Hallström BM, Stadler C, Tolf A, Lundberg E, Schwenk JM, Jirström K, Bjartell A, et al. Analysis of the human prostate-specific proteome defined by transcriptomics and antibody-based profiling identifies TMEM79 and ACOXL as two putative, diagnostic markers in prostate cancer. PLoS ONE. 2015;10(8):e0133449.
Poznyak A, Grechko AV, Poggio P, Myasoedova VA, Alfieri V, Orekhov AN. The diabetes mellitus-atherosclerosis connection: the role of lipid and glucose metabolism and chronic inflammation. Int J Mol Sci 2020, 21(5).
We thank Jichang Hu help us do biological analysis, and we also thank Jichang Hu help us improved the English expression and corrected the grammatical errors.
This work was received no funding.
Ethics approval and consent to participate
Ethical approval was not required for this study because we used a public database.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Flowchart of dataset selecting.
. Nonnegative matrix factorization (NMF) clustering was conducted and two subgroups were identified the optimal value for consensus clustering.
. List of 223 prognostic related DEGs used in this study.
. Identification of hub genes by using LASSO.
. Identification of hub genes by using SVM-RFE.
. Identification of hub genes by using RF.
About this article
Cite this article
Rao, L., Peng, B. & Li, T. Nonnegative matrix factorization analysis and multiple machine learning methods identified IL17C and ACOXL as novel diagnostic biomarkers for atherosclerosis. BMC Bioinformatics 24, 196 (2023). https://doi.org/10.1186/s12859-023-05244-w