Cardiovascular Therapeutics

Cardiovascular Therapeutics / 2020 / Article

Research Article | Open Access

Volume 2020 |Article ID 1230513 |

Bin Zhao, Dan Wang, Yanling Liu, Xiaohong Zhang, Zheng Wan, Jinling Wang, Ting Su, Linshan Duan, Yan Wang, Yuehua Zhang, Yilin Zhao, "Six-Gene Signature Associated with Immune Cells in the Progression of Atherosclerosis Discovered by Comprehensive Bioinformatics Analyses", Cardiovascular Therapeutics, vol. 2020, Article ID 1230513, 13 pages, 2020.

Six-Gene Signature Associated with Immune Cells in the Progression of Atherosclerosis Discovered by Comprehensive Bioinformatics Analyses

Academic Editor: Victor Garcia
Received20 Apr 2020
Revised15 Jun 2020
Accepted29 Jun 2020
Published01 Aug 2020


Background. As a multifaceted disease, atherosclerosis is often characterized by the formation and accumulation of plaque anchored to the inner wall of the arteries and causes some cardiovascular diseases and vascular embolism. Numerous studies have reported on the pathogenesis of atherosclerosis. However, fewer studies focused on both genes and immune cells, and the correlation of genes and immune cells was evaluated via comprehensive bioinformatics analyses. Methods. 29 samples of atherosclerosis-related gene expression profiling, including 16 human advanced atherosclerosis plaque (AA) and 13 human early atherosclerosis plaque (EA) samples from the Gene Expression Omnibus (GEO) database, were analyzed to get differentially expressed genes (DEGs) and the construction of protein and protein interaction (PPI) networks. Besides, we detected the relative fraction of 22 immune cell types in atherosclerosis by using the deconvolution algorithm of “cell type identification by estimating relative subsets of RNA transcripts (CIBERSORT).” Ultimately, based on the significantly changed types of immune cells, we executed the correlation analysis between DEGs and immune cells to discover the potential genes and pathways associated with immune cells. Results. We identified 17 module genes and 6 types of significantly changed immune cells. Correlation analysis showed that the relative percentage of T cell CD8 has negative correlation with the C1QB expression (, ), and the relative percentage of macrophage M2 has positive correlation with the CD86 expression (, ) in EA. Meanwhile, four gene expressions (CD53, C1QC, NCF2, and ITGAM) have a high correlation with the percentages of T cell CD8 and macrophages (M0 and M2) in AA samples. Conclusions. In this study, we suggested that the progression of atherosclerosis might be related to CD86, C1QB, CD53, C1QC, NCF2, and ITGAM and that it plays a role in regulating immune-competent cells such as T cell CD8 and macrophages M0 and M2. These results will enable studies of the potential genes associated with immune cells in the progression of atherosclerosis, as well as provide insight for discovering new treatments and drugs.

1. Background

Atherosclerosis is a multifaceted, progressive, and chronic inflammatory arterial disease that is recognized to be the leading cause of morbidity and mortality around the world [1]). It is characterized by the formation and build-up of atherosclerosis plaque inside the damaged arteries [2, 3]. Plaque is composed of low-density lipoprotein (LDL) cholesterol, fat, calcium, and other substances existing in the blood, which can harden and narrow the arteries [47]. Many studies have shown that atherosclerosis can affect any arterial blood vessels in the body and can lead to the atherosclerosis-related diseases, including ischemic heart, carotid artery, peripheral artery, and chronic kidney diseases [813]. Risk factors include high blood pressure, abnormal cholesterol levels, diabetes, obesity, family genetic history, smoking, age, and an unhealthy lifestyle.

Data mining has been used in various applications, including sequencing [14], microarray gene expression analysis [1517], single-nucleotide polymorphism detection [18, 19], and genomic loss and amplification (copy number variation) analysis [20, 21]. Using microarrays, integrated bioinformatics enables researchers to quickly identify differentially expressed target genes between atherosclerosis samples in a single experiment [22, 23]. CIBERSORT is a deconvolution computational method for quantifying immune cell fractions from bulk tissue gene expression profiles. This method can accurately calculate the relative proportion of 22 types of immune cell compositions in lesion samples [24, 25].

The detailed mechanism of the pathogenesis of atherosclerosis is unclear. Although studies have revealed that chronic inflammation can drive atherosclerosis, which is the leading cause of cardiovascular disease as confirmed by molecular and cellular experiments, fewer studies have been conducted to analyze the correlation of genes and immune cells in atherosclerosis-related big data.

In this study, we reanalyzed the GSE28829 dataset reported previously in Doring et al.’s team research [2] and detected potential target genes for atherosclerosis treatment from the perspective of big data analysis. We firstly identified candidate DEGs and significant immune cells. Then, we explored the correlation between gene expressions and the relative percentages of immune cells to identify potential gene signatures useful for the diagnosis and therapeutic treatment of atherosclerosis.

2. Materials and Methods

2.1. Data Acquisition

The robust multiarray averaging normalized microarray expression profile GSE28829 [2] and affiliated annotation file were downloaded from the National Center Biotechnology Information Gene Expression Omnibus ( website [26], which was tested on the GPL570 platform based on the Affymetrix Human Genome U133 Plus 2.0 array. GSE28829 contains 13 early atherosclerotic plaque samples (EA group) and 16 advanced atherosclerotic plaque samples (AA group) from the human carotid artery. Figure 1(a) provides an overview of the analysis workflow.

2.2. Data Preprocessing

After the GSE28829 expression matrix was downloaded, probe identification was matched to the corresponding gene symbol. For multiprobes to one gene, we retained the probe showing a significant gene expression value after deleting the non-mRNA probe. Based on this gene expression matrix information, we identified the significant differentially expressed genes and immune cells.

2.3. Identification of DEGs

The limma package was utilized to identify differentially expressed genes (DEGs) between advanced atherosclerotic plaques and early atherosclerotic plaque samples in RStudio [2729]. The criteria were as follows: (1) the adjusted , a moderate -test corrected by Benjamini and Hochberg’s method [30]; (2) log fold change (FC) of or log fold change (FC) of .

2.4. Gene Ontology and Pathway Analysis

Gene Ontology (GO) is used to describe the roles of genes and gene products in any organism based on existing biological knowledge and is divided into three independent branches: biological process (BP), cellular component (CC), and molecular function (MF) (Harris et al., 2004; Harris et al., 2006). Metabolic pathways and gene signaling networks based on available databases such as KEGG [31] and Reactome [32] were used to describe the pathway enrichment analyses. We used the DAVID website (https://david. for gene annotation and visualization to perform the GO and pathway analysis. , calculated via Fisher’s exact test [33], was used as the threshold for statistical significance [32, 34].

2.5. PPI Network Construction and Module Analysis

First, the identified DEGs were uploaded to the STRING [35] (version 11.0) website which includes 2 billion interactions associated with 24.6 million proteins referred to 5090 organs. STRING was used to determine the PPIs between DEG-encoded proteins. Second, the minimum interaction score was set to 0.4. The PPI networks were constructed using Cytoscape software [36]. The built-in Molecular Complex Detection (MCODE), a well-known automated method for detecting highly interconnected subgraphs as molecular complexes or clusters in large PPI networks was utilized to screen the modules in the PPI network. The correlated parameter criteria were set by default, except . Moreover, functional enrichment analysis was performed for DEGs in the significant module with , calculated via Fisher’s exact test [34], as the cutoff criterion.

2.6. Immune Cell Infiltration Analysis

Normalized gene expression data were utilized to evaluate the relative proportions of 22 types of infiltrating immune cells via using the CIBERSORT algorithm [25]. The gene expression matrix was uploaded to the CIBERSORT online website ( by setting the default signature matrix at 1000 permutations. CIBERSORT is a deconvolution algorithm that depends on a set of reference gene expression values (a “signature matrix” of 547 genes in 22 types of immune cells). Next, significant immune cells between EA and AA samples were identified with the threshold Wilcoxon test at .

2.7. Correlation Analysis of Genes and Immune Cells

Pearson correlation test analysis was carried out to illustrate the relationship between gene expressions and the relative percentages of immune cells in EA and AA samples, respectively [37]. The value of the correlation coefficient between gene expressions and the relative proportion of immune cells could indicate the strong, weak, or no correlation. Based on the paired -test, was considered statistically significant.

2.8. Statistical Analysis

The moderate -test was used to identify differentially expressed genes. Fisher’s exact test was applied to perform GO and KEGG analysis. The Wilcoxon test was applied to immune cell analysis. Paired -test was applied to correlation analysis between genes and cells. All statistical analyses were carried out in R version 3.5.2 software.

3. Results

3.1. Identification of DEGs

In the study, we identified 91 differentially expressed genes (DEGs) in the AA group compared to the EA group (Figure 1(b) and Table 1). Among them, 59 DEGs were upregulated (adjusted and ), and the remaining 32 DEGs were downregulated (adjusted and ).

DEGsGene name


3.2. GO and Pathway Analysis

DEGs were uploaded to the DAVID website to identify the GO and pathway terms. As shown in Figure 1(c), significantly enriched GO and pathway DEGs were involved in defense response (BP), extracellular space (CC), immunoglobulin receptor binding (MF), and Staphylococcus aureus infection (KEGG pathway). As shown in Table 2, the significant GO terms of upregulated DEGs were mainly enriched in defense response (BP), immune response (BP), and regulation of immune response (BP), while significant GO terms of downregulated DEGs mainly were enriched in muscle contraction (BP), muscle system process (BP), and movement of cell or subcellular component (BP). The pathway terms of upregulated DEGs were enriched in Staphylococcus aureus infection (KEGG), phagosome (KEGG), and leukocyte transendothelial migration (KEGG), while the pathways of downregulated DEGs were unavailable (Table 3).

CategoryTermCount value

Upregulated genes
GOTERM_BP_FATGO:0006952~defense response341.78-22
GOTERM_BP_FATGO:0006955~immune response329.60-20
GOTERM_BP_FATGO:0050776~regulation of immune response234.74-15
GOTERM_BP_FATGO:0045087~innate immune response211.47-13
GOTERM_BP_FATGO:0002684~positive regulation of immune system process221.68-13
GOTERM_BP_FATGO:0002682~regulation of immune system process252.27-13
GOTERM_BP_FATGO:0050778~positive regulation of immune response188.88-12
GOTERM_BP_FATGO:0048584~positive regulation of response to stimulus261.26-10
GOTERM_BP_FATGO:0009605~response to external stimulus261.99-10
GOTERM_BP_FATGO:0002250~adaptive immune response142.78-10
GOTERM_BP_FATGO:0006954~inflammatory response164.15-10
GOTERM_BP_FATGO:0002253~activation of immune response156.56-10
GOTERM_BP_FATGO:0002764~immune response-regulating signaling pathway144.79-09
GOTERM_BP_FATGO:0007166~cell surface receptor signaling pathway277.57-09
GOTERM_BP_FATGO:0050900~leukocyte migration121.43-08
GOTERM_BP_FATGO:0002757~immune response-activating signal transduction132.47-08
GOTERM_BP_FATGO:0032101~regulation of response to external stimulus141.75-07
GOTERM_BP_FATGO:0002252~immune effector process142.44-07
GOTERM_CC_FATGO:0044421~extracellular region part334.63-07
GOTERM_CC_FATGO:0005615~extracellular space204.84-07
Downregulated genes
GOTERM_BP_FATGO:0006936~muscle contraction50.00187531
GOTERM_BP_FATGO:0003012~muscle system process50.003885779
GOTERM_BP_FATGO:0006928~movement of cell or subcellular component90.006091447
GOTERM_BP_FATGO:0015672~monovalent inorganic cation transport50.007385488
GOTERM_BP_FATGO:0002028~regulation of sodium ion transport30.007948344
GOTERM_BP_FATGO:0043269~regulation of ion transport50.011584817
GOTERM_BP_FATGO:0010959~regulation of metal ion transport40.014867703
GOTERM_MF_FATGO:0003779~actin binding40.017959842
GOTERM_BP_FATGO:0006812~cation transport60.018744883
GOTERM_BP_FATGO:0030007~cellular potassium ion homeostasis20.019293154
GOTERM_BP_FATGO:0042391~regulation of membrane potential40.019968055
GOTERM_BP_FATGO:0034765~regulation of ion transmembrane transport40.024929627
GOTERM_BP_FATGO:0034762~regulation of transmembrane transport40.02729303
GOTERM_BP_FATGO:0055075~potassium ion homeostasis20.028804816
GOTERM_MF_FATGO:0008092~cytoskeletal protein binding50.029503547
GOTERM_BP_FATGO:0071805~potassium ion transmembrane transport30.034876856
GOTERM_BP_FATGO:0071804~cellular potassium ion transport30.034876856
GOTERM_BP_FATGO:0032412~regulation of ion transmembrane transporter activity30.03593287
GOTERM_BP_FATGO:0048738~cardiac muscle tissue development30.036287643

PathwayTermCount value

Upregulated genes
KEGG_PATHWAYhsa05150: Staphylococcus aureus infection72.40-07
KEGG_PATHWAYhsa04145: phagosome79.36-05
KEGG_PATHWAYhsa04670: leukocyte transendothelial migration62.60-04
KEGG_PATHWAYhsa05133: pertussis55.28-04
KEGG_PATHWAYhsa04060: cytokine-cytokine receptor interaction70.001252687
KEGG_PATHWAYhsa04062: chemokine signaling pathway60.002299717
KEGG_PATHWAYhsa05134: legionellosis40.002506044
BioCartah_blymphocytePathway: B lymphocyte cell surface molecules30.003628217
KEGG_PATHWAYhsa04610: complement and coagulation cascades40.005025012
KEGG_PATHWAYhsa05323: rheumatoid arthritis40.009854418
KEGG_PATHWAYhsa05152: tuberculosis50.011906527
KEGG_PATHWAYhsa04672: intestinal immune network for IgA production30.023560079
KEGG_PATHWAYhsa04380: osteoclast differentiation40.028350582
KEGG_PATHWAYhsa05416: viral myocarditis30.033703679
KEGG_PATHWAYhsa04514: cell adhesion molecules (CAMs)40.034833291
Downregulated genes

No significant signal pathway ( value < 0.05) available.
3.3. PPI Network Construction and Module Analysis

After uploading the 91 DEGs into the STRING online database and downloading the TSV format file of the interaction of multiple genes to Cytoscape software for PPI network construction, 59 DEGs (48 upregulated and 11 downregulated genes) were filtered from the 91 DEGs to construct the PPI networks, which contained 59 nodes/genes and 306 edges (Figure 2(a)); 32 genes did not participate in the PPI networks. Among these 59 nodes/genes, 17 central nodes/genes in module 1 were identified by the MCODE app and significantly associated with immune system function (Figure 2(b), Figure S1).

3.4. Immune Cell Infiltration Analysis

We first used the CIBERSORT algorithm to investigate the relative proportion of the 22 subpopulations of immune cells among EA and AA samples (Figure 3(a)). The relative proportions of 6 types of immune cells were significantly different between the EA and AA groups (Figure 3(b)). The cell types were T cell CD8 (), T cell gamma delta (), monocytes (), macrophage M0 (), macrophage M2 (), and dendritic cells (activated) (). Among these 6 types of immune cells, we found that T cell CD8, monocytes, and dendritic cells (activated) in the EA group were present at higher fractions than in the AA group, while the other three types of immune cells showed the opposite results (Figure 3(b)).

3.5. Correlation Analysis of Genes and Immune Cells

Correlation analysis (Pearson test) was carried out to illustrate and display the relationship between candidate genes and immune cells in EA and AA samples (Figures 4 and 5). As shown, CD86 and macrophage M2 (, ) and C1QB and T cell CD8 (, ) represented good correlation (Figures 4(a) and 4(b)) in early atherosclerosis plaque samples. Moreover, most genes have a close correlation with immune cells in advanced atherosclerosis plaque samples. It is worth noting that four common genes (CD53, C1QC, NCF2, and ITGAM) from module 1 have a close correlation with T cell CD8 and macrophages M0 and M2 (Figures 5(a) and 5(b)). However, we have discarded that situation: although the value was less than 0.05 and had a high absolute value of , the scatter plot shows that the point distribution was aggregated to zero (supplement Figures S2-3).

4. Discussion

Atherosclerosis is a disease caused by plaque accumulation within the arteries. Current studies have confirmed that immune cells including dendritic cells, several T cells, monocyte/macrophage subsets, and neutrophils are associated with atherosclerosis [2, 3840]. It has been shown that specific therapies targeting the pro/anti-inflammatory cytokines such as CCL2, TNFα, and IL-6 have suggested slowing in the progression of atherosclerosis in animal models and might improve cardiovascular outcomes in human subjects in large-scale phase III trials. [41, 42]. Notably, the monoclonal antibody canakinumab, targeting to IL-1β, has reduced the risk of adverse cardiovascular events [41, 42].

In the current study, we aimed to identify the potential molecular gene signatures associated with the immune system during the progression of atherosclerotic disease. We first figured out the genes in module 1 (17 genes) and significantly changed types of immune cells (6 types of immune cells) between advanced atherosclerosis and early atherosclerosis. Afterward, according to the correlation analysis between genes and immune cells, we inferred that CD86 and C1QB have a good correlation with macrophage M2 and T cell CD8 in EA, respectively. What is more, most of the genes have associated with T cell CD8, macrophages (M0 and M2), and four common genes (CD53, C1QC, NCF2, and ITGAM), and all have correlation with the three types of immune cells in AA.

CD86 (cluster of differentiation 86) is a protein encoded by CD86, which is expressed on antigen-presenting cells (APCs) and provides costimulatory signals to T cells [43]. Meletta et al. used CD86/CD80 as a probe for atherosclerosis imaging [43]. Transfer of native Foxp3+ T cells showed a protective effect against experimental atherosclerosis (Ait-Oufella et al.; [44, 45]). CD53 (leukocyte surface antigen) is a member of the “tetraspan family” of membrane proteins and is expressed on various immune cells [46, 47]. CD53 can contribute to improving the transduction of CD2-generated signals in T cells and natural killer cells [48]. C1QB (complement C1q B chain) or C1QC (complement C1q C chain) encodes the C-chain or B-chain polypeptide of serum complement subcomponent C1q, respectively, and deficiency of C1q is associated with glomerulonephritis and lupus erythematosus. The Bos et al. team has suggested that C1QB might be associated with atherosclerosis and coronary artery disease [49]. What is more, the Khoonsari et al. group has revealed that the lower levels of C1QB and C1QC were involved in cell adhesion, migration, regulation of the synapse, and the immune system [50]. NCF2 (neutrophil cytosol factor 2) encodes a subunit of NADPH oxidase, and mutation in this gene can result in chronic granulomatous disease [51]. However, no research has revealed that NCF2 was involved in atherosclerosis. ITGAM (integrin alpha M) is known as complement receptor 3A (CR3A) or cluster of differentiation molecule 11B (CD11B) [52] and primarily expressed on the surface of innate immune cells [53]. Recent reports revealed that NCF2 and ITGAM play significant immune-regulatory roles in autoimmune disease [54, 55]. To date, fewer studies were focused on how genes (C1QB, C1QC, NCF2, and ITGAM) regulate immune cells (T cell CD8 and macrophages M0 or M2) and their relationship between gene expressions and the relative percentages of immune cells.

Based on this fact, (1) atherosclerotic diseases are related to immune cells and (2) there are examples of researchers using modified genes (such as IL-1β, CCL2, and IL-6) to treat atherosclerotic diseases [41, 42]. The innovation of this study is to screen out gene signatures associated with immune cells in the progression of atherosclerosis. According to the relationship between gene expressions and the relative proportions of immune cells, these genes may interact with immune cells through some unknown cell membrane receptors or ligands. It may be possible to modify the interactions of these genes with membrane receptors or ligands to identify new therapeutic treatments for atherosclerosis, as well as mechanisms of atherosclerosis regulation.

Nevertheless, there exist some inevitable difficulties in this study that should be taken into consideration. For instance, the atherosclerosis-related datasets were fewer and not as easily acquired and collected from the open public database as cancer datasets, which leads to the lack of comprehensiveness of the study which cannot be used to verify our results. Although the limited sample size of atherosclerosis may reduce the confidence, the approaches and ideas in this study are helpful in enlightening the inspiration of other researchers. Of course, additional molecular and cellular experiments should be performed to assess their characteristics.

5. Conclusions

The identification of 6 specific gene signatures and correlated immune cells in the progression of atherosclerosis may give us a clue to explore the mechanism of cardiovascular disease and its therapeutic treatments.


NCBI:National Center for Biotechnology Information
GEO:Gene Expression Omnibus
DEGs:Differentially expressed genes
PPI:Protein and protein interaction
CIBERSORT:Cell type identification by estimating relative subsets of RNA transcripts
LDL:Low-density lipoprotein
EA:Early atherosclerosis
AA:Advanced atherosclerosis
GO:Gene Ontology
KEGG:Kyoto Encyclopedia of Genes and Genomes
BP:Biological process
CC:Cellular component
MF:Molecular function
MCODE:Molecular complex detection.

Data Availability

The data associated with this article has been deposited in the NCBI-GEO website (

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Authors’ Contributions

BZ performed comparative analysis using bioinformatics tools. DW, YLL, XHZ, ZW, JLW, TS, LSD, YW, and YHZ participated in data analysis and discussion. BZ interpreted the data and wrote the manuscript. YLZ organized and offered funding for the project. All authors read and approved the final manuscript.


We also would like to thank Huabo Yang from the College of Foreign Languages and Cultures, Xiamen University, for correcting the grammar in the manuscript. We also thank Yingying Zhang (the General Hospital of PLA Rocket Force) for the statistics analysis. This study was supported by the National Natural Science Foundation of China (No. 81770294).

Supplementary Materials

Figure S1: the chord plot for functional enrichments of module 1 genes. Figure S2: the scatterplot of correlation between CTSS expression and the relative proportion of T cell gamma delta. Gray-shaded areas in scatterplots represent the standard errors of the blue regression lines. : correlation coefficient. Figure S3: the scatterplot of correlation between 17 gene expressions and the proportion of dendritic cells activated. Gray-shaded areas in scatterplots represent the standard errors of the blue regression lines. : correlation coefficient. (Supplementary Materials)


  1. I. Hilgendorf, F. K. Swirski, and C. S. Robbins, “Monocyte fate in atherosclerosis,” Arteriosclerosis, Thrombosis, and Vascular Biology, vol. 35, no. 2, pp. 272–279, 2015. View at: Publisher Site | Google Scholar
  2. Y. Döring, H. D. Manthey, M. Drechsler et al., “Auto-antigenic protein-DNA complexes stimulate plasmacytoid dendritic cells to promote atherosclerosis,” Circulation, vol. 125, no. 13, pp. 1673–1683, 2012. View at: Publisher Site | Google Scholar
  3. K. J. Woollard and F. Geissmann, “Monocytes in atherosclerosis: subsets and functions,” Nature Reviews Cardiology, vol. 7, no. 2, pp. 77–86, 2010. View at: Publisher Site | Google Scholar
  4. N. Arnold and W. Koenig, “Atherosklerose als inflammatorische Erkrankung – Pathophysiologie, klinische Relevanz und therapeutische Implikationen,” Deutsche Medizinische Wochenschrift, vol. 144, no. 5, pp. 315–321, 2019. View at: Publisher Site | Google Scholar
  5. W. Jin, W. Wu, K. Yang et al., “The single nucleotide polymorphisms of chromosome 9p21 and CD147 were relevant with the carotid plaque risk in acute cerebral infarction patients among Chinese Han population,” Journal of Molecular Neuroscience, vol. 70, pp. 1–11, 2020. View at: Publisher Site | Google Scholar
  6. J. Nilsson and G. K. Hansson, “Vaccination strategies and immune modulation of atherosclerosis,” Circulation Research, vol. 126, no. 9, pp. 1281–1296, 2020. View at: Publisher Site | Google Scholar
  7. T. N. Tulenko, L. Laury-Kleintop, M. F. Walter, and R. P. Mason, “Cholesterol, calcium and atherosclerosis: is there a role for calcium channel blockers in atheroprotection?” International Journal of Cardiology, vol. 62, pp. S55–S66, 1997. View at: Publisher Site | Google Scholar
  8. V. Aboyans, J.-B. Ricco, M.-L. E. L. Bartelink et al., “2017 ESC guidelines on the diagnosis and treatment of peripheral arterial diseases, in collaboration with the European Society for Vascular Surgery (ESVS),” European Heart Journal, vol. 39, no. 9, pp. 763–816, 2018. View at: Publisher Site | Google Scholar
  9. S. S. Anand, S. Yusuf, V. Vuksan et al., “Differences in risk factors, atherosclerosis, and cardiovascular disease between ethnic groups in Canada: the Study of Health Assessment and Risk in Ethnic groups (SHARE),” The Lancet, vol. 356, no. 9226, pp. 279–284, 2000. View at: Publisher Site | Google Scholar
  10. M. H. Criqui and V. Aboyans, “Correction,” Circulation Research, vol. 117, no. 1, pp. e12–e26, 2015. View at: Publisher Site | Google Scholar
  11. R. A. DeFronzo and E. Ferrannini, “Insulin resistance: a multifaceted syndrome responsible for NIDDM, obesity, hypertension, dyslipidemia, and atherosclerotic cardiovascular disease,” Diabetes Care, vol. 14, no. 3, pp. 173–194, 1991. View at: Publisher Site | Google Scholar
  12. I. J. Kullo and T. W. Rooke, “Peripheral artery disease,” New England Journal of Medicine, vol. 374, no. 9, pp. 861–871, 2016. View at: Publisher Site | Google Scholar
  13. A. S. Levey, K. U. Eckardt, Y. Tsukamoto et al., “Definition and classification of chronic kidney disease: a position statement from Kidney Disease: Improving Global Outcomes (KDIGO),” Kidney International, vol. 67, no. 6, pp. 2089–2100, 2005. View at: Publisher Site | Google Scholar
  14. A. A. Brown, A. Viñuela, O. Delaneau, T. D. Spector, K. S. Small, and E. T. Dermitzakis, “Predicting causal variants affecting expression by using whole-genome sequencing and RNA-seq from multiple human tissues,” Nature Genetics, vol. 49, no. 12, pp. 1747–1751, 2017. View at: Publisher Site | Google Scholar
  15. J. Martorell-Marugan, D. Toro-Dominguez, M. E. Alarcon-Riquelme, and P. Carmona-Saez, “MetaGenyo: a web tool for meta-analysis of genetic association studies,” BMC Bioinformatics, vol. 18, no. 1, p. 563, 2017. View at: Publisher Site | Google Scholar
  16. R. L. Stears, T. Martinsky, and M. Schena, “Trends in microarray analysis,” Nature Methods, vol. 9, no. 1, pp. 140–145, 2003. View at: Publisher Site | Google Scholar
  17. D. Toro-Domínguez, J. Martorell-Marugán, R. López-Domínguez et al., “ImaGEO: integrative gene expression meta-analysis from GEO database,” Bioinformatics, vol. 35, no. 5, pp. 880–882, 2019. View at: Publisher Site | Google Scholar
  18. J. B. Fan, X. Chen, M. K. Halushka et al., “Parallel genotyping of human SNPs using generic high-density oligonucleotide tag arrays,” Genome Research, vol. 10, no. 6, pp. 853–860, 2000. View at: Publisher Site | Google Scholar
  19. D. G. Wang, J. B. Fan, C. J. Siao et al., “Large-scale identification, mapping, and genotyping of single-nucleotide polymorphisms in the human genome,” Science, vol. 280, no. 5366, pp. 1077–1082, 1998. View at: Publisher Site | Google Scholar
  20. T. Walsh, J. M. McClellan, S. E. McCarthy et al., “Rare structural variants disrupt multiple genes in neurodevelopmental pathways in schizophrenia,” Science, vol. 320, no. 5875, pp. 539–543, 2008. View at: Publisher Site | Google Scholar
  21. D. E. Watkins-Chow and W. J. Pavan, “Genomic copy number and expression variation within the C57BL/6J inbred mouse strain,” Genome Research, vol. 18, no. 1, pp. 60–66, 2007. View at: Publisher Site | Google Scholar
  22. L. Perisic, S. Aldi, Y. Sun et al., “Gene expression signatures, pathways and networks in carotid atherosclerosis,” Journal of Internal Medicine, vol. 279, no. 3, pp. 293–308, 2016. View at: Publisher Site | Google Scholar
  23. X. W. Tan, X. T. Zhang, L. L. Pan, X. X. Tian, and P. Z. Dong, “Identification of key pathways and genes in advanced coronary atherosclerosis using bioinformatics analysis,” BioMed Research International, vol. 2017, 12 pages, 2017. View at: Publisher Site | Google Scholar
  24. A. J. Gentles, A. M. Newman, C. L. Liu et al., “The prognostic landscape of genes and infiltrating immune cells across human cancers,” Nature Medicine, vol. 21, no. 8, pp. 938–945, 2015. View at: Publisher Site | Google Scholar
  25. A. M. Newman, C. L. Liu, M. R. Green et al., “Robust enumeration of cell subsets from tissue expression profiles,” Nature Methods, vol. 12, no. 5, pp. 453–457, 2015. View at: Publisher Site | Google Scholar
  26. T. Barrett, D. B. Troup, S. E. Wilhite et al., “NCBI GEO: mining tens of millions of expression profiles--database and tools update,” Nucleic Acids Research, vol. 35, pp. D760–D765, 2007. View at: Publisher Site | Google Scholar
  27. B. Phipson, S. Lee, I. J. Majewski, W. S. Alexander, and G. K. Smyth, “Robust hyperparameter estimation protects against hypervariable genes and improves power to detect differential expression,” The Annals of Applied Statistics, vol. 10, no. 2, pp. 946–963, 2016. View at: Publisher Site | Google Scholar
  28. K. C. Ruegg, E. C. Anderson, R. J. Harrigan et al., “Genetic assignment with isotopes and habitat suitability (GAIAH), a migratory bird case study,” Methods in Ecology and Evolution, vol. 8, no. 10, pp. 1241–1252, 2017. View at: Publisher Site | Google Scholar
  29. G. K. Smyth, “Linear models and empirical Bayes methods for assessing differential expression in microarray experiments,” Statistical Applications in Genetics and Molecular Biology, vol. 3, no. 1, pp. 1–25, 2004. View at: Publisher Site | Google Scholar
  30. Y. Benjamini and Y. Hochberg, “Controlling the false discovery rate - a practical and powerful approach to multiple testing,” Journal of the Royal Statistical Society, vol. 57, no. 1, pp. 289–300, 1995. View at: Publisher Site | Google Scholar
  31. M. Kanehisa, S. Goto, S. Kawashima, Y. Okuno, and M. Hattori, “The KEGG resource for deciphering the genome,” Nucleic Acids Research, vol. 32, no. 90001, pp. 277D–2280, 2004. View at: Publisher Site | Google Scholar
  32. G. Joshi-Tope, M. Gillespie, I. Vastrik et al., “Reactome: a knowledgebase of biological pathways,” Nucleic Acids Research, vol. 33, pp. D428–D432, 2004. View at: Publisher Site | Google Scholar
  33. R. A. Fisher, “On the interpretation of χ2 from contingency tables, and the calculation of P,” Journal of the Royal Statistical Society, vol. 85, no. 1, pp. 87–94, 1922. View at: Publisher Site | Google Scholar
  34. G. Dennis, B. T. Sherman, D. A. Hosack et al., “DAVID: Database for Annotation, Visualization, and Integrated Discovery,” Genome Biology, vol. 4, no. 9, p. R60, 2003. View at: Publisher Site | Google Scholar
  35. A. Franceschini, D. Szklarczyk, S. Frankild et al., “STRING v9.1: protein-protein interaction networks, with increased coverage and integration,” Nucleic Acids Research, vol. 41, no. D1, pp. D808–D815, 2012. View at: Publisher Site | Google Scholar
  36. M. Kohl, S. Wiese, and B. Warscheid, “Cytoscape: software for visualization and analysis of biological networks,” Methods in Molecular Biology, vol. 696, pp. 291–303, 2011. View at: Publisher Site | Google Scholar
  37. P. Sedgwick, “Pearson's correlation coefficient,” British Medical Journal, vol. 345, 2012. View at: Publisher Site | Google Scholar
  38. M. M. Vazquez, M. V. Gutierrez, S. R. Salvatore et al., “Nitro-oleic acid, a ligand of CD36, reduces cholesterol accumulation by modulating oxidized-LDL uptake and cholesterol efflux in RAW264.7 macrophages,” Redox Biology, vol. 36, article 101591, 2020. View at: Publisher Site | Google Scholar
  39. C. Weber, A. Zernecke, and P. Libby, “The multifaceted contributions of leukocyte subsets to atherosclerosis: lessons from mouse models,” Nature Reviews Immunology, vol. 8, no. 10, pp. 802–815, 2008. View at: Publisher Site | Google Scholar
  40. A. Zernecke, I. Bot, Y. Djalali-Talab et al., “Protective role of CXC receptor 4/CXC ligand 12 unveils the importance of neutrophils in atherosclerosis,” Circulation Research, vol. 102, no. 2, pp. 209–217, 2008. View at: Publisher Site | Google Scholar
  41. J. Khambhati, M. Engels, M. Allard-Ratick, P. B. Sandesara, A. A. Quyyumi, and L. Sperling, “Immunotherapy for the prevention of atherosclerotic cardiovascular disease: promise and possibilities,” Atherosclerosis, vol. 276, pp. 1–9, 2018. View at: Publisher Site | Google Scholar
  42. P. M. Ridker and T. F. Luscher, “Anti-inflammatory therapies for cardiovascular disease,” European Heart Journal, vol. 35, no. 27, pp. 1782–1791, 2014. View at: Publisher Site | Google Scholar
  43. R. Meletta, A. M. Herde, P. Dennler, E. Fischer, R. Schibli, and S. D. Krämer, “Preclinical imaging of the co-stimulatory molecules CD80 and CD86 with indium-111-labeled belatacept in atherosclerosis,” EJNMMI Research, vol. 6, no. 1, p. 1, 2016. View at: Publisher Site | Google Scholar
  44. A. Mor, D. Planer, G. Luboshits et al., “Role of naturally occurring CD4+ CD25+ regulatory T cells in experimental atherosclerosis,” Arteriosclerosis Thrombosis and Vascular Biology, vol. 27, no. 4, pp. 893–900, 2007. View at: Publisher Site | Google Scholar
  45. H. Ait-Oufella, B. L. Salomon, S. Potteaux et al., “Natural regulatory T cells control the development of atherosclerosis in mice,” Nature Medicine, vol. 12, no. 2, pp. 178–180, 2006. View at: Publisher Site | Google Scholar
  46. H. S. Jin, J. E. Cho, and S. Park, “Association between CD53 genetic polymorphisms and tuberculosis cases,” Genes & Genomics, vol. 41, no. 4, pp. 389–395, 2019. View at: Publisher Site | Google Scholar
  47. V. Korinek and V. Horejsi, “Genomic structure of the human CD53 gene,” Immunogenetics, vol. 38, no. 4, pp. 272–279, 1993. View at: Publisher Site | Google Scholar
  48. C. J. Comerci, E. M. Mace, P. P. Banerjee, and J. S. Orange, “CD2 promotes human natural killer cell membrane nanotube formation,” PLoS One, vol. 7, no. 10, article e47664, 2012. View at: Publisher Site | Google Scholar
  49. S. Bos, M. Phillips, G. F. Watts, A. J. M. Verhoeven, E. J. G. Sijbrands, and N. C. Ward, “Novel protein biomarkers associated with coronary artery disease in statin- treated patients with familial hypercholesterolemia,” Journal of Clinical Lipidology, vol. 11, no. 3, pp. 682–693, 2017. View at: Publisher Site | Google Scholar
  50. P. E. Khoonsari, A. Häggmark, M. Lönnberg et al., “Analysis of the cerebrospinal fluid proteome in Alzheimer’s disease,” PLoS One, vol. 11, no. 3, article e0150672, 2016. View at: Publisher Site | Google Scholar
  51. R. T. Kenney and T. L. Leto, “A HindIII polymorphism in the human NCF2 gene,” Nucleic Acids Research, vol. 18, no. 23, p. 7193, 1990. View at: Publisher Site | Google Scholar
  52. D. A. Solovjov, E. Pluskota, and E. F. Plow, “Distinct roles for the α and β subunits in the functions of integrin αMβ2,” Journal of Biological Chemistry, vol. 280, no. 2, pp. 1336–1345, 2005. View at: Publisher Site | Google Scholar
  53. S. Park, H. Sim, H. I. Kim et al., “CD11b regulates antibody class switching via induction of AID,” Molecular Immunology, vol. 87, pp. 47–59, 2017. View at: Publisher Site | Google Scholar
  54. J. Chou, J. T. Hsu, W. Bainter, R.'a. Al-Attiyah, W. Al-Herz, and R. S. Geha, “A novel mutation in _NCF2_ associated with autoimmune disease and a solitary late-onset infection,” Clinical Immunology, vol. 161, no. 2, pp. 128–130, 2015. View at: Publisher Site | Google Scholar
  55. F. Rosetti and T. N. Mayadas, “The many faces of Mac-1 in autoimmune disease,” Immunological Reviews, vol. 269, no. 1, pp. 175–193, 2016. View at: Publisher Site | Google Scholar

Copyright © 2020 Bin Zhao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.