Long Wang1,2
Mingyue Wu
Jinyu Huang- 1Key Laboratory of Clinical Cancer Pharmacology and Toxicology Research of Zhejiang Province, Department of Cardiology, Affiliated Hangzhou First People's Hospital, Zhejiang University School of Medicine, Hangzhou, China
- 2Key Laboratory of Clinical Cancer Pharmacology and Toxicology Research of Zhejiang Province, Center for Translational Medicine, Affiliated Hangzhou First People's Hospital, Zhejiang University School of Medicine, Hangzhou, China
- 3Department of Cardiology, Nanjing University Medical School, Nanjing, China
- 4Department of Cardiology, Jiangsu Taizhou People's Hospital, Taizhou, China
- 5Key Laboratory of Combined Multi-Organ Transplantation, Ministry of Public Health, The First Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China
- 6Institute of Translational Medicine, Zhejiang University, Hangzhou, China
Since immune infiltration is closely associated with the progression and prognosis of atherosclerosis, we aimed to describe the abundance of 24 immune cell types within atherosclerotic tissues. In the current study, we used the Immune Cell Abundance Identifier (ImmuCellAI), a web-based tool, to estimate the abundance of 24 immune cells based on the microarray profiles of atherosclerotic carotid artery samples to analyze the proportions and the dysregulation of immune cell types within carotid atherosclerosis. We found that atherosclerotic immune cells had a diverse landscape dominated by T cells and myeloid cells and that macrophages and dendritic cells (DCs) showed different abundance in normal and atherosclerotic tissues. Moreover, the expression of macrophages was closely related to the level of the expression of DCs and of exhausted T cells, while the expression of T-helper type 1 (Th1) cells was strongly correlated with the expression of T-helper type 2 (Th2) cells and effector memory cells. Our data confirm a distinct profile of atherosclerosis-infiltrating immune cell subpopulations, which may inspire an immunological direction for research on atherosclerosis.
Introduction
The immune system, comprising various immune mediators and cells, is carefully orchestrated and is critical for the host defense (1). Innate immune cells [e.g., macrophages, neutrophils, dendritic cells (DCs), and natural killer (NK) cells] and adaptive immune cells (e.g., B and T cells) play pivotal roles in cancers and chronic inflammatory diseases (2, 3). Immune cell dysfunction, including abnormal distributions of abundance and type, has also been recognized to be associated with atherosclerosis (4). Consequently, investigation of the abundance of immune cells and of alterations in atherosclerotic (AS) tissues could provide insights into the pathogenesis, development, regression, and treatment of atherosclerosis.
Atherosclerosis is an inflammatory disease, with various stages of the disease involving immune cells (5). Subendothelial deposition of modified lipoproteins acts as a damage-associated molecular pattern to stimulate and recruit monocytes and elicit vascular inflammation. Then, monocytes are infiltrated to differentiate locally into macrophages, whose dysfunctional lipid metabolism and reduced efferocytosis led to unresolvable inflammation (6). Previous studies confirmed T cells as critical drivers and modifiers of the pathogenesis of atherosclerosis. CD4+ T-helper type 1 (Th1) cells promote the development of inflammatory atheroma while regulatory T (Tregs) cells have antiatherogenic functions (7). Little is understood regarding the role of other Th-cell subsets, such as T-helper type 2 (Th2), Th9, Th17, and Tfh cells, and the role of CD8+ T cells in atherosclerosis, which seem to have both proatherogenic and antiatherogenic functions (8). Several B lymphocyte subsets contribute to the inflammatory process through cellular and humoral responses (9). DCs are important modulators of atherosclerosis, which connect the innate and the adaptive immune responses (10). Compared with the Cardiovascular Inflammation Reduction Trial (CIRT) (11), the Canakinumab Anti-inflammatory Thrombosis Outcome Study (CANTOS) (12) resulted in fewer cardiovascular events in patients with stable atherosclerosis, which may be attributed to lessening the level of interleukin (IL)-1β, IL-6, or C-reactive protein. Therefore, it provided proof of the principle that anti-inflammatory treatments targeting special immune cells, pathways, or cytokines will have the potential for presenting atherosclerotic events in the future.
High-throughput sequencing technologies, such as microarray, RNA-sequencing, and single-cell sequencing, provide numerous gene expression profiles suitable for the study of the distribution of immune cells. Some methods, including cell-type identification by estimating relative subsets of RNA transcript (CIBERSORT) (13), Estimating the Proportions of Immune and Cancer cells (EPIC) (14), and Tumor Immune Estimation Resource (TIMER) (15), have advanced the notion of immune infiltration by enumerating immune cells from the bulk transcriptome data of tumor samples, and considering the critical role of T-cell subsets in atherosclerosis (8), we aimed to explore the proportion and alteration of immune cells in atherosclerosis. In the present study, we introduced the Immune Cell Abundance Identifier (ImmuCellAI), a web-based tool, to estimate the abundance of 24 immune cells from a gene expression data set, including 18 T-cell subsets in atherosclerosis (16).
In this study, we downloaded the microarray profiles of normal and AS tissues of the carotid artery from a gene expression omnibus (GEO) data set and applied the ImmuCellAI method to investigate the quantity of 24 immune cell types within these tissues. Furthermore, we explored the differences in immune cells and their functions in the initiation and development of atherosclerosis.
Materials and Methods
Data Acquisition
The data used in this study were acquired from the GEO database (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE100927). The GSE100927 microarray profile was generated on the GPL17077 platform (Agilent-039494 SurePrint G3 Human GE v2 8x60K Microarray) from the samples of human peripheral arteries consisting of atherosclerotic plaques and control tissues. Atheromatous plaques were harvested and collected from patients who underwent carotid, femoral, or infrapopliteal endarterectomy, while atherosclerotic lesion-free healthy arteries were obtained from organ donors. Demographic and clinical details about this biocollection have been published in another study (17). We selected all 25 carotid tissues, including 11 normal and 14 stable atherosclerotic plaque samples, for a carotid atherosclerosis analysis.
Evaluation of Infiltrating Immune Cells in Atherosclerosis
The ImmuCellAI was used to calculate the fractions of infiltrating immune cells with the default parameters. Recent studies using single-cell RNA-sequencing, ATAC-sequencing, and mass cytometry have shown a dominance of T cells in human carotid plaques (18–21), therefore, the focus of the study is on T cells. ImmuCellAI is an online tool that estimates the abundance of 24 immune cell types using the gene expression data based on a gene set signature method, including 18 T-cell subtypes and 6 other immune cells: B cells, NK cells, monocytes, macrophages, neutrophils, and DCs (16). The flow cytometric definition of immune cells was referenced at http://bioinfo.life.hust.edu.cn/ImmuCellAI. The infiltration score was defined as the sum of all percentages of 24 infiltrating immune cells.
Visualization of the Filtered Data
The filtered data were used to visualize the immune cell percentage in tissues by using R version 3.6.1. The ggplot2 package was used to perform principal component analysis (PCA) with the percent values of immune cells within each sample as the input. The bar plot and the heatmap were plotted to show the proportions of 24 immune cell types. The correlation heatmap was drawn by using the corrplot package to visualize the correlations between 22 types of infiltrated immune cells. The violin plots showed that a differential expression of 24 types of immune cells was plotted by using the vioplot package.
Statistical Analysis
The Wilcoxon test was conducted to compare the differences in cell composition between the two groups. The value of p < 0.05 was considered statistically significant.
Results
Before performing an ImmuCellAI analysis, the probe names from the microarrays were replaced with the corresponding gene symbols, and the data were normalized. The ImmuCellAI method allowed a novel speculation about the immune cell infiltration in atherosclerotic carotid tissues. About 22 subpopulations of immune cells in 25 carotid artery samples (11 normal tissues and 14 atherosclerotic tissues) are shown in Figure 1 and Table 1. We found that T cells, especially Th17, γδT, Th2, and Th1 cells, dominated in the carotid artery. Because of a significant variation of immune cells in atherosclerosis, we inferred that a variation in immune cell proportions may be used to describe the difference between individuals. The abundance of immune cells from 25 samples showed a significant group-bias clustering in the PCA plot (Figure 2). The first two principal components could explain most of the data variation. In addition, we further explored the correlation between the subpopulations of immune cells and found that some subgroups had obvious correlations (Figure 3). The presence of macrophages was closely related to the expression of DCs and exhausted T cells, while Th1 cells were associated with Th2 cells and effector memory cells. Compared with normal samples, atherosclerosis samples had increased infiltration scores, which suggest a greater abundance of immune cells and altered proportions of immune cell subpopulations (Table 1, Figure 4). Atherosclerotic arteries have a higher abundance of exhausted T cells, DC cells, mucosal-associated invariant T (MAIT) cells, B cells, and macrophages. Conversely, the proportion of Th cells decreased in atherosclerosis, which could be due to the overall increase in immune cell abundance. Finally, based on a hierarchical clustering of 24 immune cell subpopulations (Figure 5), we found some modules where specific types of immune cells had a similar expression, which suggest a synergetic role in atherosclerosis.
Figure 1. A bar plot of immune infiltration in normal samples and atherosclerosis samples. About 22 subpopulations of immune cells in 25 carotid artery samples [11 normal tissues (CT) and 14 atherosclerotic tissues (AS)] are shown.
Table 1. Immune cell abundance identifier (ImmuCellAI) fraction of infiltration of immune cells in atherosclerotic and normal carotid tissues.
Figure 2. A principal component analysis (PCA) is performed on normal and atherosclerotic arteries. The first two principal components explaining the most data variation are shown.
Figure 3. A correlation matrix of immune cell proportions between normal and atherosclerotic tissues.
Figure 4. A violin plot of the immune cell proportions in the two groups. The normal control group is marked in blue, and the AS group is marked in red. The values of p < 0.05 were considered to be statistically significant.
Figure 5. Heatmap of the proportions of 20 immune cell types in the normal group and the AS group.
Discussion
Inherent immune diversity is a pivotal component of atherosclerosis that contributes to lesion progression or regression and predicts the clinical outcomes (20). Therefore, the identification of a specific immune dysregulation at the lesion site could provide insights into the progression of atherosclerosis. In this study, we provided a description of infiltrating immune cells in atherogenic plaques and analyzed the differences in the abundance of immune cells between the groups. These immune cells can be considered as either proatherogenic or antiatherogenic until functional experimental validation in the future.
By using the ImmuCellAI method, we can directly compare the differences in the profiles of 24 immune cells between control and atherosclerosis samples. We found that DCs, Th17 cells, γδT cells, Th2 cells, and CD4+ T cells were the main immune cell subpopulations in carotid artery atherosclerosis, which were consistent with the fact that T cells represent the largest population of leukocytes in atherosclerotic plaques (22). Furthermore, AS tissues have higher percentages of MAIT cells, DCs, B cells, and macrophages and lower proportions of Th1 cells, Th2 cells, Th17 cells, effector memory cells, and neutrophils. These results appear to be in contrast with the current knowledge that proatherogenic cells, including Th1 (23) and Th17 cells (24), are found at higher levels in normal tissues. The immune cell number is a relative percentage among each group. AS tissues have more immune cell infiltration, so there is a relatively lower content of Th1 and Th17 cells in these tissues than in the control tissues, which mainly contain macrophages. Due to an increased infiltration of immune cells, the decrease in percentage may not reflect a reduced number of Thlls. In case of adjusting to the infiltration score, the absolute percentage of Th cells in the AS group may not be lower than that in the normal group.
Macrophages participate in the initiation of atherosclerosis when subjected to the stimulation with oxidized lipids, cholesterol crystals, and inflammatory cytokines (25). B lymphocyte subsets are also activated and participate in the regulation of the inflammatory process (25). B1 cells have been demonstrated to inhibit lesion formation, whereas B2 cells have been shown to promote lesion formation (9). The abundance of DCs is dramatically increased during the atherosclerosis development, and DCs participate in all stages of atherosclerosis (26). DCs link innate and adaptive immune responses by presenting antigens to T lymphocytes (27) and may promote or limit atherogenesis by modulating these responses (28). CD11b+CD11c+ cells directly participate in foam cell formation (29), while conventional DCs (cDCs) interact with T cells and natural killer T (NKT) cells, resulting in an increased secretion of interferon γ (IFNγ) and IL-17 by T cells. Mice deficient in IFNγ had a lower atherosclerotic lesion burden, which suggested that IFNγ had pathogenic effects on atherosclerosis (30). The effects of Th17 cells on atherosclerosis are controversial. Some studies have shown a proatherogenic role for cytokines secreted by Th17 cells (31), while other studies have suggested that Th17 cells are associated with reduced disease development (32). Plasmacytoid DCs (pDCs) produce IFN-α and IFN-β, which play a proatherogenic role (33). Moreover, the presence of macrophages was closely related to the presence of DCs and exhausted T cells, while Th1 cells were associated with Th2 cells and effector memory cells. Th1 cells promote the progression of inflammatory atheroma, while Tregs inhibit both the innate and adaptive inflammatory responses (5). CD4+ T and CD8+ T cells were associated with increased atherosclerosis in ApoE−/− mice (34). These results suggest that they may play a synergetic role in atherosclerosis. Single-cell sequencing reported that, among aortic CD45+ cells, the abundance of Th2 cells (as antiatherogenic cells) was more in the normal tissue while the abundance of macrophages was less in the normal tissue. Moreover, we used another deconvolution method (CIBERSORT) and found that the abundance of CD4 T cells is more in normal tissues than in AS tissues. However, CIBERSORT does not include Th1, Th2, or Th17 cells.
There were several limitations to our study. First, the included samples were unpaired, and we lacked the information about relevant plaque phenotypes, such as stability, vulnerability, or calcification. Second, the ImmuCellAI method focused on T-cell subsets and underestimated some macrophage or DC subsets. Finally, this was a bioinformatic analysis, and the predictions needed to be confirmed by immunofluorescent staining and flow cytometry in the future.
Conclusion
In summary, with the help of the ImmuCellAI analysis, we revealed the immune cell composition of human carotid artery atherosclerosis and identified dysregulated immune cells. Although the study revealed some differences in respect to the literature studies, the different immunological composition of the atherosclerotic lesion suggests further investigation in the atherosclerosis studies.
Data Availability Statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.
Author Contributions
WY, PL, and JH designed the study. LW performed a bioinformatic analysis. MW and BG contributed to manuscript revision. All authors approved the final submitted manuscript.
Funding
This research was supported by the National Natural Science Foundation of China (Grant No. 81800284), the Ministry of Science and Technology of China (Grant No. 2019YFC0120705), the Science and Technology Foundation of Zhejiang Province of China (Grant Nos. BQ20H020001 and 2020C03018), the Science and Technology Foundation of Hangzhou of China (Grant No. 20191203B82), and the Natural Science Foundation of Jiangsu Province of China (Grant No. BK20180144).
Conflict 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.
References
1. Chaplin DD. Overview of the immune response. J Allergy Clin Immunol. (2010) 125(2 Suppl. 2):S3–23. doi: 10.1016/j.jaci.2009.12.980
2. Brodin P, Davis MM. Human immune system variation. Nat Rev Immunol. (2017) 17:21–9. doi: 10.1038/nri.2016.125
3. Gajewski TF, Schreiber H, Fu YX. Innate and adaptive immune cells in the tumor microenvironment. Nat Immunol. (2013) 14:1014–22. doi: 10.1038/ni.2703
4. Schaftenaar F, Frodermann V, Kuiper J, Lutgens E. Atherosclerosis: the interplay between lipids and immune cells. Curr Opin Lipidol. (2016) 27:209–15. doi: 10.1097/MOL.0000000000000302
5. Libby P, Lichtman AH, Hansson GK. Immune effector mechanisms implicated in atherosclerosis: from mice to humans. Immunity. (2013) 38:1092–104. doi: 10.1016/j.immuni.2013.06.009
6. Tabas I, Bornfeldt KE. Macrophage phenotype and function in different stages of atherosclerosis. Circ Res. (2016) 118:653–67. doi: 10.1161/CIRCRESAHA.115.306256
7. Tabas I, Lichtman AH. Monocyte-Macrophages and T cells in atherosclerosis. Immunity. (2017) 47:621–34. doi: 10.1016/j.immuni.2017.09.008
8. Saigusa R, Winkels H, Ley K. T cell subsets and functions in atherosclerosis. Nat Rev Cardiol. (2020) 17:387–401. doi: 10.1038/s41569-020-0352-5
9. Tsiantoulas D, Diehl CJ, Witztum JL, Binder CJ. B cells and humoral immunity in atherosclerosis. Circ Res. (2014) 114:1743–56. doi: 10.1161/CIRCRESAHA.113.301145
10. Cybulsky MI, Cheong C, Robbins CS. Macrophages and dendritic cells: partners in atherogenesis. Circ Res. (2016) 118:637–52. doi: 10.1161/CIRCRESAHA.115.306542
11. Ridker PM, Everett BM, Pradhan A, MacFadyen JG, Solomon DH, Zaharris E, et al. Low-Dose methotrexate for the prevention of atherosclerotic events. N Engl J Med. (2019) 380:752–62. doi: 10.1056/NEJMoa1809798
12. Ridker PM, Everett BM, Thuren T, MacFadyen JG, Chang WH, Ballantyne C, et al. Antiinflammatory therapy with canakinumab for atherosclerotic disease. N Engl J Med. (2017) 377:1119–31. doi: 10.1056/NEJMoa1707914
13. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. (2015) 12:453–7. doi: 10.1038/nmeth.3337
14. Racle J, de Jonge K, Baumgaertner P, Speiser DE, Gfeller D. Simultaneous enumeration of cancer and immune cell types from bulk tumor gene expression data. Elife. (2017) 6:e26476. doi: 10.7554/eLife.26476.049
15. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. (2017) 77:e108–10. doi: 10.1158/0008-5472.CAN-17-0307
16. Miao YR, Zhang Q, Lei Q, Luo M, Xie GY, Wang H, et al. ImmuCellAI: a unique method for comprehensive T-cell subsets abundance prediction and its application in cancer immunotherapy. Adv Sci. (2020) 7:1902880. doi: 10.1002/advs.201902880
17. Steenman M, Espitia O, Maurel B, Guyomarch B, Heymann MF, Pistorius MA, et al. Identification of genomic differences among peripheral arterial beds in atherosclerotic and healthy arteries. Sci Rep. (2018) 8:3940. doi: 10.1038/s41598-018-22292-y
18. Depuydt MAC, Prange KHM, Slenders L, Ord T, Elbersen D, Boltjes A, et al. Microanatomy of the human atherosclerotic plaque by single-cell transcriptomics. Circ Res. (2020) 127:1437–55. doi: 10.1161/CIRCRESAHA.120.316770
19. Fernandez DM, Rahman AH, Fernandez NF, Chudnovskiy A, Amir ED, Amadori L, et al. Single-cell immune landscape of human atherosclerotic plaques. Nat Med. (2019) 25:1576–88. doi: 10.1038/s41591-019-0590-4
20. Libby P. Superficial erosion and the precision management of acute coronary syndromes: not one-size-fits-all. Eur Heart J. (2017) 38:801–3. doi: 10.1093/eurheartj/ehw599
21. Winkels H, Ehinger E, Vassallo M, Buscher K, Dinh HQ, Kobiyama K, et al. Atlas of the immune cell repertoire in mouse atherosclerosis defined by single-cell RNA-sequencing and mass cytometry. Circ Res. (2018) 122:1675–88. doi: 10.1161/CIRCRESAHA.117.312513
22. Winkels H, Wolf D. Heterogeneity of T cells in atherosclerosis defined by single-cell RNA-sequencing and cytometry by time of flight. Arterioscler Thromb Vasc Biol. (2020) 41:ATVBAHA120312137. doi: 10.1161/ATVBAHA.120.312137
23. Buono C, Binder CJ, Stavrakis G, Witztum JL, Glimcher LH, Lichtman AH. T-bet deficiency reduces atherosclerosis and alters plaque antigen-specific immune responses. Proc Natl Acad Sci USA. (2005) 102:1596–601. doi: 10.1073/pnas.0409015102
24. Zhang W, Tang T, Nie D, Wen S, Jia C, Zhu Z, et al. IL-9 aggravates the development of atherosclerosis in ApoE-/- mice. Cardiovasc Res. (2015) 106:453–64. doi: 10.1093/cvr/cvv110
25. Kuznetsova T, Prange KHM, Glass CK, de Winther MPJ. Transcriptional and epigenetic regulation of macrophages in atherosclerosis. Nat Rev Cardiol. (2020) 17:216–28. doi: 10.1038/s41569-019-0265-3
26. Koltsova EK, Ley K. How dendritic cells shape atherosclerosis. Trends Immunol. (2011) 32:540–7. doi: 10.1016/j.it.2011.07.001
27. Zhao TX, Mallat Z. Targeting the immune system in atherosclerosis: JACC state-of-the-art review. J Am Coll Cardiol. (2019) 73:1691–706. doi: 10.1016/j.jacc.2018.12.083
28. Ait-Oufella H, Sage AP, Mallat Z, Tedgui A. Adaptive (T and B cells) immunity and control by dendritic cells in atherosclerosis. Circ Res. (2014) 114:1640–60. doi: 10.1161/CIRCRESAHA.114.302761
29. Miller YI, Choi SH, Wiesner P, Fang L, Harkewicz R, Hartvigsen K, et al. Oxidation-specific epitopes are danger-associated molecular patterns recognized by pattern recognition receptors of innate immunity. Circ Res. (2011) 108:235–48. doi: 10.1161/CIRCRESAHA.110.223875
30. Whitman SC, Ravisankar P, Daugherty A. IFN-gamma deficiency exerts gender-specific effects on atherogenesis in apolipoprotein E-/- mice. J Interferon Cytokine Res. (2002) 22:661–70. doi: 10.1089/10799900260100141
31. Smith E, Prasad KM, Butcher M, Dobrian A, Kolls JK, Ley K, et al. Blockade of interleukin-17A results in reduced atherosclerosis in apolipoprotein E-deficient mice. Circulation. (2010) 121:1746–55. doi: 10.1161/CIRCULATIONAHA.109.924886
32. Taleb S, Romain M, Ramkhelawon B, Uyttenhove C, Pasterkamp G, Herbin O, et al. Loss of SOCS3 expression in T cells reveals a regulatory role for interleukin-17 in atherosclerosis. J Exp Med. (2009) 206:2067–77. doi: 10.1084/jem.20090545
33. Goossens P, Gijbels MJ, Zernecke A, Eijgelaar W, Vergouwe MN, van der Made I, et al. Myeloid type I interferon signaling promotes atherosclerosis by stimulating macrophage recruitment to lesions. Cell Metab. (2010) 12:142–53. doi: 10.1016/j.cmet.2010.06.008
Keywords: immune infiltration, immune cells, atherosclerosis, gene expression, T cells
Citation: Wang L, Gao B, Wu M, Yuan W, Liang P and Huang J (2021) Profiles of Immune Cell Infiltration in Carotid Artery Atherosclerosis Based on Gene Expression Data. Front. Immunol. 12:599512. doi: 10.3389/fimmu.2021.599512
Received: 16 October 2020; Accepted: 09 February 2021;
Published: 25 March 2021.
Edited by:
Silvia Piconese, Sapienza University of Rome, ItalyReviewed by:
Harlan Pierre Jones, University of North Texas Health Science Center, United StatesDavide Flego, Sapienza University of Rome, Italy
Copyright © 2021 Wang, Gao, Wu, Yuan, Liang and Huang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wei Yuan, eXVhbndlaWh0JiN4MDAwNDA7MTI2LmNvbQ==; Ping Liang, cGluZ2xpYW5nJiN4MDAwNDA7emp1LmVkdS5jbg==; Jinyu Huang, emRzeWhqeTA5MDImI3gwMDA0MDt6anUuZWR1LmNu