Immune infiltrating cells in cholangiocarcinoma may become clinical diagnostic markers: based on bioinformatics analysis.

doi:10.21203/rs.3.rs-54963/v2

Download PDF

Research

Immune infiltrating cells in cholangiocarcinoma may become clinical diagnostic markers: based on bioinformatics analysis.

https://doi.org/10.21203/rs.3.rs-54963/v2

This work is licensed under a CC BY 4.0 License

Journal Publication

published 22 Feb, 2021

Read the published version in World Journal of Surgical Oncology →

You are reading this older preprint version

Read the latest preprint version →

Background: intrahepatic cholangiocarcinoma (ICC) is a malignant tumor originating from the secondary bile duct and its branch epithelium. Among primary liver tumors, the incidence of ICC is second only to hepatocellular carcinoma. tumor microenvironment can regulate the occurrence and development of tumors. This study is dedicated to finding more markers that can diagnose ICC by finding the differential tumor microenvironment cells between ICC and normal tissues.

Methods: We wanted to study the infiltration of immune cells between the cholangiocarcinoma of the same patient and its paired non-tumor tissues, to explore the difference of immune cells in the tumor microenvironment and adjacent non-tumor tissues

in the same organism. So, we searched the relevant data of patients with ICC from the GEO database and found that the GSE45001 data set meets our research needs. CIBERSORT database is used to calculate immune cell composition. Finally, perform visual analysis and data statistics to find out the differentially expressed immune cells.

Results: we found that the expression levels of Dendritic cells activated, Macrophages M2 and T cells regulatory (Tregs) in ICC were higher than normal tissues and the expression levels of Macrophages M1, Monocytes and T cells follicular helper in ICC were lower than normal tissues.

Conclusion: These 6 types of immune cells are expected to become molecular markers for clinical diagnosis of ICC.

General Surgery

Oncology

Bioinformatics analysis

prognosis

intrahepatic cholangiocarcinoma

immune cells.

Intrahepatic cholangiocarcinoma (ICC) is a malignant tumor originating from the secondary bile duct and its branch epithelium [1]. Among primary liver tumors, the incidence of ICC is second only to hepatocellular carcinoma, and the incidence of ICC in China is on the rise. Early diagnosis can give ICC patients better therapeutic effects, so it is very important to find more markers to diagnose ICC [2]. Recent studies have shown that the tumor microenvironment can regulate the occurrence and development of tumors. This study is dedicated to finding more markers that can diagnose ICC by finding the differential tumor microenvironment cells between tumors and normal tissues [3].

In this study, by searching the ICC related data set in the GEO database ，We established relevant selection and exclusion criteria 1. must be paired intrahepatic cholangiocarcinoma and para-cancer tissue; 2. The quantity must be at least 10 pairs; 3. These data sets have not previously been used to study immune cell infiltration in ICC [4]. it was found that GSE45001 [5] .(https://www.ncbi.nlm.nih.gov/gds/?term=GSE45001) met the research conditions (with 10 tumor samples and 10 normal tissue samples). We put it into the CIBERSORT (https://cibersort.stanford.edu/) [6] database to analyze their immune cell composition. Then, take histograms, principal component maps, heat maps [7] and other visualizations to show the differences in immune cells between tumor tissues and normal tissues. We found that there are indeed some differences in the content of immune cells between ICC tissues and normal tissues, and through co-expression heat maps, violin diagrams, etc., more accurately extract the difference between tumor tissues and normal tissues.

This study is dedicated to finding new markers that can diagnose ICC. The results show that the expression differences between Dendritic cells activated, Macrophage M1, Macrophage M2, Monocytes, T cells follicular helper and T cells regulatory (Tregs) between ICC and normal tissue are statistically significant (P <0.05), are expected to become a new diagnostic ICC marker.

2.1 Data collection and download

We searched the relevant data of patients with ICC from the GEO database and found that the GSE45001 data set meets our research needs. It contains data for 10 pairs of tumor tissues and normal tissues.

2.2 Evaluation of the content of immune cells in the sample.

CIBERSORT is a database that can calculate the immune cell content of each sample, from which we obtain the information of the immune cell content of the sample we want. We will download 20 samples of the data set, calculate the relevant immune cell content in CIBERSORT.

2.3 Drawing of heat maps related to samples and immune cells.

Use the R language pheatmap package to analyze and draw a heat map of the relationship between immune cells and sample expression.

2.4 Drawing and analysis of the correlation heat map of immune cells.

The corrplot package of the R language was used to analyze the correlation of the immune cells expressed in the samples in this study, and draw a correlation heat map.

2.5 The violin diagram of the expression of immune cells in adjacent tissues and ICC tissues.

Use the vioplot package of R language to analyze the sample's immune cells and draw a violin chart (Wilcox Test was used for statistical analysis, P value has been adjusted).. Therefore, yellow represents the tumor tissue group, and green represents the normal tissue group.

2.6 The paired difference graphs explain the expression of immune cells in the sample.

Use R language related codes to draw paired difference maps. The two shapes connected by black lines are two types of tissues of the same patient.

2.7 Survival analysis and Kaplan-Meier curve drawing.

Analyze the direct relationship between immune cells and survival prognosis in Tumor Immune Estimation Resource (http://timer.cistrome.org/)(TIMER 2.0) [8-9] database, and draw K-M curve.

2.8 Sample principal component analysis chart.

Using R language code, the main component analysis of the sample. Use two circles to distinguish between the tumor group and the normal group, and one dot represents one sample.

2.9 Statistics description

P<0.05 is considered statistically significant. P<0.01 is considered statistically significant.

3.1 Download and analysis of data.

We searched the GEO database for samples that meet the research requirements and found that GSE45001 meets our needs. It contains 10 ICC samples and 10 normal samples (Table1).

Table 1. Data from GSE45001; 20 samples from 10 paired CCA patients

GSM1095633	Normal Stroma of CCA	GSM1095634	Tumoral Stroma of CCA
GSM1095635	Normal Stroma of CCA	GSM1095636	Tumoral Stroma of CCA
GSM1095637	Normal Stroma of CCA	GSM1095638	Tumoral Stroma of CCA
GSM1095639	Normal Stroma of CCA	GSM1095640	Tumoral Stroma of CCA
GSM1095641	Normal Stroma of CCA	GSM1095642	Tumoral Stroma of CCA
GSM1095643	Normal Stroma of CCA	GSM1095644	Tumoral Stroma of CCA
GSM1095645	Normal Stroma of CCA	GSM1095646	Tumoral Stroma of CCA
GSM1095647	Normal Stroma of CCA	GSM1095648	Tumoral Stroma of CCA
GSM1095649	Normal Stroma of CCA	GSM1095650	Tumoral Stroma of CCA
GSM1095651	Normal Stroma of CCA	GSM1095652	Tumoral Stroma of CCA

3.2 Analysis of the content of immune cells in the sample.

After quality inspection, we found that 7 pairs of samples were available, we downloaded them and included them in the study. The 7 normal samples and 7 ICC samples obtained by the correction were put into the CIBERSORT database to calculate the immune cell content, and the obtained results were plotted as a histogram. Then, we further draw the results into a principal component diagram, from which we can find that the tumor tissue group and the normal tissue group have significantly different immune cell content (Figure 1). Finally, we draw a heat map and visually analyze the expression of immune cells in the sample (Figure 2). The results show that the expression of some immune cells shows differences between tumor tissues and normal tissues.

3.3 Co-expression analysis between immune cells.

Co-expression analysis of the expression levels of immune cells obtained in this study showed that the correlation between T cells regulatory (Tregs) [10-13]and Neutrophils [14-17] was the strongest, with a coefficient of 0.78, and the two cells with the weakest correlation were Macrophages [18-19] and T cells follicular helper [20-21], the correlation coefficient is -0.59 (Figure 3).

3.4 Immune cells with significantly different expressions in ICC and normal tissues are screened.

By drawing a violin chart (Figure 4), We found that the expression levels of Dendritic cells activated (Figure5A), Macrophages M2 (Figure5B) and T cells regulatory (Tregs) (Figure5C) in ICC were higher than normal tissues and the expression levels of Monocytes (Figure5D), T cells follicular helper (Figure5E) and Macrophages M1 (Figure5F) in ICC were lower than normal tissues

3.5 Principal Component Analysis (PCA)

After obtaining the matrix of immune cells, we wondered whether these immune cells could distinguish between the normal group and the tumor group. Then we did principal component analysis. The dimensionality was reduced to PCA1 and PCA2 by PCA, and the X-axis was labeled as PCA1 and Y-axis as PCA2. Then an ellipse was simulated for the normal group and the tumor group, respectively. If the two ellipses did not cross, it suggested that the 22 immune cells could distinguish the normal group from the tumor group well. We used the "GGplot2" package for analysis, and the results showed that the two ellipses did not cross (Figure 6), indicating that the 22 kinds of immune cells in this study could well distinguish the tumor group from the normal group.

3.6 Immune cells related to ICC survival prognosis.

We further explore the relationship between each immune cell in Figure 5 and ICC survival in the Timer2.0 database. It was found that the survival rate of the group with high Monocyte cell content was significantly better than that of the group with low Monocyte cell content.

The survival rate of the group with high T cell regulatory content was better than that of the group with low regulatory content. However, the data of other cells is not sufficient at present, and further research and supplement are needed.

Intrahepatic cholangiocarcinoma (ICC) is a type of primary malignant tumor of the liver that originates from the bile duct and its branches. Because the clinical symptoms of ICC are not obvious early, it often causes patients to lose the opportunity for timely surgery. Therefore, it is extremely critical to find a better method for early diagnosis of ICC. The tumor microenvironment is closely related to the occurrence and development of tumors. Many studies have shown that the tumor microenvironment plays a role in promoting tumorigenesis and development, such as promoting angiogenesis, producing growth factors that promote tumors, causing hypoxic environments and inflammatory reactions. However, ICC-related research has not explored the relationship between the tumor microenvironment and ICC. With the development of biochip technology, more and more bioinformatics databases and bioinformatics tools are available for our research and use. Therefore, we can analyze the cells that constitute the immune microenvironment of ICC through bioinformatics technology to find tumor markers that can better diagnose ICC.

In this study, we retrieved the ICC data set from the GEO database and found that GSE45001 met our research needs. Then we downloaded and grouped them reasonably, 10 samples were classified into tumor group, and 10 samples were classified into normal tissue group. Put these 20 samples into the CIBERSORT database to analyze the content of immune cells and draw the principal component map. The results show that the composition of immune cells between tumor tissues and normal tissues is not the same. Through the heat map we draw, we can find that some immune cells are expressed in tumors higher than normal tissues, but some immune cells are lower than normal tissues. From the co-expression heat maps, we draw, we can better discover which immune cells are more or less relevant in this study, and drawing a violin chart can more intuitively find the expression of each immune cell in different tissue. Finally, we found that the expression of Dendritic cells activated, Macrophages M2, T cells regulatory (Tregs) in tumor tissue was higher than that in normal tissue, while the expression of Macrophages M1, Monocytes, T cells follicular helper in tumor tissue was lower than that in normal tissue. This difference provides us with more basis for clinical diagnosis of ICC.

Dendritic cells have long been considered to be anti-tumor immune cells, but studies have found that Dendritic cells can promote tumor suppression of immunity. Some people have achieved the idea of inhibiting tumor growth by suppressing the number of Dendritic cells, and believe that Dendritic cells-related pathways may become new therapeutic targets [22]. Macrophages M1 cells can accelerate the apoptosis process of tumor cells, in line with our research results, its expression in normal tissues is higher, and it has a certain effect on inhibiting tumor cells. However, Macrophages M2 and Macrophages M1 have opposite functions. They can delay the apoptosis of tumor cells and promote the progression of tumors [23]. Monocytes can not only promote tumor growth by increasing angiogenesis, but also inhibit tumor growth by identifying and clearing tumor escape signal pathways. In this study, although the content of Monocytes was higher in normal tissues than in tumor tissues in general, there were obviously some of the opposite. According to our statistical analysis, Monocytes inhibit tumor growth in ICC [24]. By querying in the TIMER 2.0 database, two immune cells (Monocytes, T cells regulatory) were found to be related to the survival prognosis of ICC patients, and there was statistical significance. And these two cells also fit the logic of this study that both of these cells are protective factors for ICC.

In this study, we analyzed the components of the immune microenvironment of ICC and found that 6 types of immune cells are differentially expressed between ICC and normal tissues. These 6 types of immune cells are expected to become molecular markers for clinical diagnosis of ICC.

Intrahepatic cholangiocarcinoma (ICC); T cells regulatory (Tregs)；Tumor Immune Estimation Resource (TIMER)；Kaplan-Meier curve (K-M curve)；Principal Component Analysis (PCA).

Consent for publication: All authors agree to publish the paper.

Availability of data and material: The data used to support the findings of this study are included within the article.

Competing interests: The authors declare that there is no conflict of interests.

Funding: Natural Science Research Foundation of Anhui Province and Natural Science Foundation of China

Authors' contributions: Yongwei Zhang is responsible for writing and submitting the papers; Jun Li and Wei Dai are responsible for data analysis and collation; Sihan Chen are responsible for the production of pictures; Yeben Qian is responsible for the manuscript fees and ideas guidance.

Acknowledgements: Thanks to the Natural Science Research Foundation of Anhui Province and Natural Science Foundation of China for supporting this research.
Compliance with ethics requirements: This article does not contain any studies with human or animal subjects. All data comes from public database, So the study didn't require patients to sign a consent form.

Gao Y-X, Yang T-W, Yin J-M, Yang P-X, Kou B-X, Chai M-Y, Liu X-N, Chen D-X: Progress and prospects of biomarkers in primary liver cancer (Review). Int J Oncol 2020, 57(1):54-66.
Jiang M, Deng H, Liu S, Su S: Letter regarding Risk factors for intrahepatic and extrahepatic cholangiocarcinoma: A systematic review and meta-analysis. J Hepatol 2020, 72(6):1216-1216.
Wang G, Wang Q, Liang N, Xue H, Yang T, Chen X, Qiu Z, Zeng C, Sun T, Yuan W et al: Oncogenic driver genes and tumor microenvironment determine the type of liver cancer. Cell Death Dis 2020, 11(5).
Clough E, Barrett T: The Gene Expression Omnibus Database. Methods Mol Biol 2016, 1418:93-110.
Sulpice L, Rayar M, Desille M, Turlin B, Fautrel A, Boucher E, Llamas-Gutierrez F, Meunier B, Boudjema K, Clément B et al: Molecular profiling of stroma identifies osteopontin as an independent predictor of poor prognosis in intrahepatic cholangiocarcinoma. Hepatology 2013, 58(6):1992-2000.
Zheng R, Wan C, Mei S, Qin Q, Wu Q, Sun H, Chen CH, Brown M, Zhang X, Meyer CA et al: Cistrome Data Browser: expanded datasets and new tools for gene regulatory analysis. Nucleic Acids Res 2019, 47(D1):D729-d735.
Kellom M, Raymond J: Using Dendritic Heat Maps to Simultaneously Display Genotype Divergence with Phenotype Divergence. PLoS One 2016, 11(8):e0161292.
Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, Li B, Liu XS: TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res 2017, 77(21):e108-e110.
Li B, Severson E, Pignon JC, Zhao H, Li T, Novak J, Jiang P, Shen H, Aster JC, Rodig S et al: Comprehensive analyses of tumor immunity: implications for cancer immunotherapy. Genome Biol 2016, 17(1):174.
Duffy SS, Keating BA, Perera CJ, Moalem-Taylor G: The role of regulatory T cells in nervous system pathologies. J Neurosci Res 2018, 96(6):951-968.
Martin-Moreno PL, Tripathi S, Chandraker A: Regulatory T Cells and Kidney Transplantation. Clin J Am Soc Nephrol 2018, 13(11):1760-1764.
Peters C, Kabelitz D, Wesch D: Regulatory functions of γδ T cells. Cell Mol Life Sci 2018, 75(12):2125-2135.
Roncarolo MG, Gregori S, Bacchetta R, Battaglia M, Gagliani N: The Biology of T Regulatory Type 1 Cells and Their Therapeutic Application in Immune-Mediated Diseases. Immunity 2018, 49(6):1004-1019.
van Rees DJ, Szilagyi K, Kuijpers TW, Matlung HL, van den Berg TK: Immunoreceptors on neutrophils. Semin Immunol 2016, 28(2):94-108.
Panettieri RA, Jr.: The Role of Neutrophils in Asthma. Immunol Allergy Clin North Am 2018, 38(4):629-638.
Papayannopoulos V: Neutrophils Stepping Through (to the Other Side). Immunity 2018, 49(6):992-994.
Galkina SI, Golenkina EA, Viryasova GM, Romanova YM, Sud'ina GF: Nitric Oxide in Life and Death of Neutrophils. Curr Med Chem 2019, 26(31):5764-5780.
Varol C, Mildner A, Jung S: Macrophages: development and tissue specialization. Annu Rev Immunol 2015, 33:643-675.
Wynn TA, Vannella KM: Macrophages in Tissue Repair, Regeneration, and Fibrosis. Immunity 2016, 44(3):450-462.
Vinuesa CG, Linterman MA, Yu D, MacLennan IC: Follicular Helper T Cells. Annu Rev Immunol 2016, 34:335-368.
Kim SJ, Lee K, Diamond B: Follicular Helper T Cells in Systemic Lupus Erythematosus. Front Immunol 2018, 9:1793.
Veglia F, Gabrilovich DI: Dendritic cells in cancer: the role revisited. Curr Opin Immunol 2017, 45:43-51.
Genin M, Clement F, Fattaccioli A, Raes M, Michiels C: M1 and M2 macrophages derived from THP-1 cells differentially modulate the response of cancer cells to etoposide. BMC Cancer 2015, 15:577.
Olingy CE, Dinh HQ, Hedrick CC: Monocyte heterogeneity and functions in cancer. J Leukoc Biol 2019, 106(2):309-322.

Download PDF

Journal Publication

published 22 Feb, 2021

Read the published version in World Journal of Surgical Oncology →

Reviewer #2 agreed at journal
20 Jan, 2021
Reviewer #1 agreed at journal
20 Jan, 2021
Review #1 received at journal
20 Jan, 2021
Editor assigned by journal
19 Jan, 2021
Reviewers invited by journal
19 Jan, 2021
Submission checks completed at journal
19 Jan, 2021
Editor invited by journal
19 Jan, 2021

You are reading this older preprint version

Read the latest preprint version →

Immune infiltrating cells in cholangiocarcinoma may become clinical diagnostic markers: based on bioinformatics analysis.

Status:

Journal Publication

Version 2

Abstract

Figures

1. Background

2. Methods

3. Result

4. Discussion

5. Conclusion

Abbreviations

Declarations

References

Status:

Journal Publication

Version 2