Profiles of Immune Cell Infiltration and Immune-Related Genes in the Tumor Microenvironment of Esophageal Squamous Cell Carcinoma

doi:10.21203/rs.3.rs-122000/v1

Download PDF

Research article

Profiles of Immune Cell Infiltration and Immune-Related Genes in the Tumor Microenvironment of Esophageal Squamous Cell Carcinoma

https://doi.org/10.21203/rs.3.rs-122000/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 10 Mar, 2021

Read the published version in BMC Medical Genomics →

You are reading this latest preprint version

Backgrounds: As a complex system participating in tumor development and progression, the tumor microenvironment was poorly understood in esophageal cancer especially squamous cell carcinoma (ESCC).

Methods: ESTIMATE algorithm is used to investigate tumor-infiltrating immune cells and prognostic genes which were associated with the tumor microenvironment in ESCC.

Results: Based on the immune and stromal scores, ESCC samples were divided into high and low score groups and 299 overlapping differentially expressed genes were identified. Functional enrichment analysis showed that these genes were mainly involved in muscle-related function. Prognostic genes including COL9A3, GFRA2, and VSIG4 were used to establish a risk prediction model using Cox regression analyses. Then multivariate analysis showed that COL9A3 was an independent discriminator of a better prognosis. Kaplan-Meier survival analysis showed that the expression of COL9A3 was significantly correlated with the overall survival of ESCC patients. The area under the curve for the risk model in predicting 1- and 3- year survival rates were 0.660 and 0.942, respectively. The risk score was negatively correlated with plasma cells, while positively correlated with the proportions of activated CD4 memory T cells, M1 Macrophages and M2 Macrophages (P<0.001 for each comparison). Gene set enrichment analysis suggested that both immune response and immune system process gene sets were significantly enriched in high-risk group.

Conclusions: Our study provided a comprehensive understanding of the TME in ESCC patients. The establishment of the risk model is valuable for the early identification of high-risk patients to facilitate individualized treatment and improve the possibility of immunotherapy response.

Epigenetics & Genomics

Esophageal squamous cell carcinoma

tumor microenvironment

tumor-infiltrating immune cells

prognosis

overall survival

Esophageal cancer is a gastrointestinal malignancy with extremely aggressive nature and poor prognosis [1]. It is the eighth most common cancer and the sixth most common cause of cancer death globally [2]. In Iran, esophageal cancer is more popular than any other countries or regions in the world [3]. Classified by histology, esophageal cancer is divisible into adenocarcinoma and squamous cell carcinoma (ESCC) [4]. The effective methods for treatment of ESCC include chemotherapy or chemoradiotherapy followed by extensive surgery, which will obviously reduce health-related quality of life. Although recent developments have improved prognosis and survivorship, the molecular mechanism behind ESCC is not clear till now. Therefore, it is still important to identify potential biomarkers to increase the effectiveness of therapy and survival rate of ESCC patients. As an effective therapeutic option, immunotherapy, especially immune checkpoint inhibitors, shows obviously clinical benefits in various cancers [5, 6]. Instead of two well-known immune checkpoint molecules PD-1 (programmed cell death protein 1, also named CD279) and PD-L1 (programmed cell death-ligand 1), other molecules such as CD155, CD226, and LAG3 are also recognized as new immune-related molecules which contribute to tumor-mediated immune suppression and promote tumor immunity escape in ESCC [7]. Tumor-infiltrating immune cells (TIICs), as the main components of the tumor microenvironment (TME) which composed of stromal cells, endothelial cells, and TIICs[8], have a significant impact on tumor progression, treatment, and outcomes of patients. Recently, more researches have paid attention to antitumor immunity regulated by immune microenvironment [9–11], but the mechanism regulating the infiltration of immunocytes in ESCC is poorly understood. Thus, identifying the TME related therapeutic targets may improve immunotherapy efficacy and give a new clue for clinical strategy. Based on DNA copy number, ESTIMATE (Estimation of Stromal and Immune cells in Malignant Tumor tissues using Expression data) uses gene expression signatures to infer the fraction of stromal and immune cells in different tumor samples [12]. The researches in pancreatic adenocarcinoma, lung cancer, glioblastoma, osteosarcoma, and other types of cancers[13–16] showed that this newly developed method is reliable for molecule screen.

In this current study, we first calculated immune and stromal scores of 81 ESCC tissues in the TCGA database using ESTIMATE algorithm and then retrieved immune-associated differentially expressed genes (DEGs). Then, the correlation between immune/stromal scores and clinical characteristics, prognosis of ESCC patients including age(years), gender, pathologic TNM tumor stage, and tumor grade were analyzed respectively. A predictive risk model to estimate patient outcome was established and the associations of the TME-related risk score with the levels of TIICs and immune pathways were analyzed.

Dataset and estimation of stromal and immune scores

81 ESCC samples with clinical data including age, gender, tumor grades, tumor stages, survival and outcome were obtained from the TCGA data portal (https://tcga-data.nci.nih.gov/tcga/). Both the gene expression and clinical data used in this research are publicly available and classified as open-access. To evaluate the infiltrating levels of the immune and stromal cells in the ESCC tissues, “estimate” R package was used as a tool of ESTIAMTE algorithm for calculation of immune and stromal scores. Based on the median value of immune/stromal scores, the ESCC patients were divided into high and low score groups to identify a possible association of these scores with overall survival.

Identification of DEGs based on immune and stromal scores

We used the “limma” R package to identify the DEGs between high and low score groups. Log2FC > 1.5 or log2FC < (-1.5) and p < 0.05 were set as the threshold for genes screening.

Functional enrichment analysis

Then, the functional enrichment analyses of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway were analyzed using “clusterProfiler” package, and p < 0.05 was considered to indicate statistical significance.

Establishment of TME-related risk model and survival analysis

Using the “survival” R package, the univariate and multivariate Cox regression analyses were performed for analyzing associations of the levels of DEGs with overall survival. Consequently, powerful prognostic genes were screened out and then risk scores were calculated for all ESCC patients. Based on the median of risk scores, all patients were divided into high- (n = 41) and low- risk (n = 40) group. We used Kaplan-Meier curves to exhibit the association between risk score and overall survival. The log-rank test was employed to test the statistic difference with the significance level p < 0.05. Using the “timeROC” package, the survival receiver operating characteristic curve (ROC) with the area under the curve (AUC) value was visualized.

Gene set enrichment analysis (GSEA) was carried out to evaluate associations of immune pathways with the TME-related risk score. P < 0.05 indicated statistical significance.

Association between risk scores and the levels of TIICs

The online analytical platform CIBERSORT with an arrangement of 1000 default statistical parameter was used to quantify the relative proportions of 22 TIICs in the TME of ESCC. Wilcoxon rank-sum test was used to determine the association between risk scores and the differential proportions of 22 TIICs in ESCC tissues. P < 0.05 was considered as statistical significance. Spearman rank analysis was performed using GraphPad Prism version 8.0.0 for Windows, (GraphPad Software, San Diego, California USA, www.graphpad.com).

The immune and stromal scores were tightly associated with tumor grade.

Using the ESTIMATE algorithm, we first determined immune and stromal scores of 81 ESCC samples based on the gene expression data obtained from the TCGA database (Supplementary Table 1). Statistical analysis showed that age, gender, and TNM stage had no correlation with both immune score (Fig. 1A-1E) and stromal score (Fig. 1G-1K). But high-grade tumors (G2-G3) had higher immune and stromal scores than G1 tumors (p = 0.039 and 0.003, respectively), while G2 tumor had the highest immune score and G3 tumor had the highest stromal score (Fig. 1F, 1L).

Based on the median of immune and stromal scores of every sample, we divided 81 samples into two groups, that is, 41 samples were in the high-score group while 40 in the low-score group. Kaplan-Meier survival curves revealed that no significant results were found neither in immune score groups or stromal score groups. (p = 0.264 and 0.276, Fig. 2A-2B).

Table 1

Univariate and multivariate Cox regression analysis of prognostic factors for overall survival
Variables	Univariate analysis		Multivariate analysis
Variables	HR (95%CI)	p value	HR (95%CI)	p value
COL9A3	0.689(0.538–0.881)	0.003	0.675(0.519–0.877)	0.003
VSIG4	1.321(1.008–1.732)	0.044	1.155(0.782–1.705)	0.468
GFRA2	1.331(1.007–1.759)	0.045	1.210(0.798–1.832)	0.370
Notes: Abbreviations: HR, Hazard ratio; 95%CI, 95% confidence interval

Identification of differentially expressed genes between high and low immune/stromal scores

To identify the immune-related and stromal-related genes, we performed differential analysis using “limma” packages in R based on genes expression level. As shown in volcano plots, there were a total of 1048 and 1219 TME-related DEGs between high and low immune/stromal score groups, respectively (Fig. 3A-3B). Compared with the low score group, there were 756 up-regulated and 292 down-regulated genes in the high immune score group, while 382 up-regulated and 837 down-regulated genes in the high stromal score group (Fig. 3C-3D). Subsequently, 299 overlapping genes in Venn diagrams were selected for further analysis.

Functional enrichment analysis

Functional enrichment analysis was performed to comprehend the functional properties of the 299 TME-related DEGs using the “clusterProfiler” package in R. We selected the top 10 GO terms in each biological process (Fig. 4A), cellular component (Fig. 4B), and molecular function (Fig. 4C). The results showed that these DEGs mainly enriched in the biological process including muscle system process, muscle contraction, and striated muscle cell differentiation. Similarly, contractile fiber part, contractile fiber, and myofibril which mainly related to muscle structure were most enriched by selected DEGs. As for molecular function, substrate-specific channel activity, channel activity, and passive transmembrane transporter activity were mainly annotated. In the KEGG pathway annotation and enrichment analysis, we found that enriched pathways were associated with neuroactive ligand-receptor interaction, vascular smooth muscle contraction, and pancreatic secretion. (Fig. 4D).

Establishment of the risk prediction model

To further explore the prognostic significance of the TME-related DEGs, we revealed 3 genes including COL9A3, VSIG4 and GFRA2, which were significantly associated with overall survival (p < 0.05 for each comparison) using univariate Cox regression analysis (Table 1). And multivariate Cox regression results showed that COL9A3 was an independent prognosis factor for ESCC patients (p = 0.003). Constantly, the risk score was obtained based on relative coefficients in Cox regression according to the formula: (-0.3937 * COL9A3 expression level) + (0.1902 * GFRA2 expression level) + (0.1441 * VSIG4 expression level) (Supplementary Table 2). The median value of risk scores was used as a cutoff to divide samples into high- and low-risk groups (Fig. 5A). As shown in the scatter plot, patients with lower risk scores showed higher survival probability than those with higher risk scores (Fig. 5B). Moreover, compared with the low-risk group, higher expression levels of GFRA2 and VSIG4 were observed in the high-risk group. Conversely, the cases with the higher expression level of COL9A3 were in the low-risk group, confirming that COL9A3 is a good prognostic factor (Fig. 5C). The Kaplan-Meier survival analysis revealed that high-risk score was significantly associated with a poor outcome (Fig. 5D). As for ROC curve, the AUC for the risk model in predicting 1- and 3- year survival rates were 0.660 and 0.942, respectively (Fig. 5E), which indicated a high sensitivity and specificity of the risk model. Moreover, the survival curve showed that a high COL9A3 expression level was related to low survival probabilities (Fig. 5F).

The relationship between immune scores and TIICs

Using spearman rank analysis in GraphPad Prism 8.0, we identify the association between risk score and the infiltration levels of 22 TIICs in ESCC tissues (Table 2). The results showed that the risk score was negatively correlated with the proportions of plasma cells (r = -0.265, p = 0.017) (Fig. 6A) while positively correlated with the proportions of activated CD4 memory T cells (r = 0.247, p = 0.026), M1 Macrophages (r = 0.233, p = 0.036) and M2 Macrophages (r = 0.391, p < 0.001) (Fig. 6B-6D).

Table 2

Spearman rank analysis to determine the association between risk score and the levels of 22 TIICs in ESCC tissues
Tumor-infiltrating immune cell	Risk score
Tumor-infiltrating immune cell	Spearman r	p-value
B cells memory	0.091	0.417
B cells naive	-0.164	0.143
Dendritic cells activated	-0.047	0.674
Dendritic cells resting	0.127	0.259
Eosinophils	-0.039	0.730
Macrophages M0	-0.159	0.157
Macrophages M1	0.233	0.036
Macrophages M2	0.391	< 0.001
Mast cells activated	-0.172	0.124
Mast cells resting	0.188	0.093
Monocytes	-0.021	0.852
Neutrophils	0.088	0.436
NK cells activated	0.057	0.616
NK cells resting	-0.046	0.686
Plasma cells	-0.265	0.017
T cells CD4 memory activated	0.247	0.026
T cells CD4 memory resting	-0.178	0.112
T cells CD4 naive	-0.098	0.384
T cells CD8	0.105	0.349
T cells follicular helper	-0.114	0.311
T cells gamma delta	-0.014	0.899
T cells regulatory (Tregs)	-0.033	0.771

The involvement of immune pathways predicted by risk model

With online database Molecular Signatures Databases v4.0 (http://software.broadinstitute.org/gsea/downloads.jsp), we retrieved two immune gene sets, involving immune response (M19817) and immune system process (M13664). GSEA results showed that both immune response and immune system process gene sets were significantly enriched in the high-risk group (p = 0.010 and 0.012, respectively) (Fig. 7A-7B).

TME is one of the main hallmarks of cancer, so it is important to identify the key druggable factors and pathways in the TME. The response to cancer immunotherapy especially immune checkpoint inhibitors were impacted by tumor immune microenvironment. Wang et al. found that an “infiltrated-excluded” or “cole” tumor immune microenvironment is predictive of poor response and low-dose metformin reprograms the TME in ESCC [17]. Research by Strizova et al. showed that FasR + NK cells, CD4+, and CD8 + T cells infiltrated lymph nodes at the lowest levels and that the FasR + DR3 + CD4 + T cells were enhanced in esophageal cancer [18]. These compartmental proportions correlated with tumor stage and tumor grade suggested new possibilities for personalized immunotherapy for patients.

Therefore, in the current study, we conducted the bioinformatics analysis of TIICs and TME-related genes in ESCC and established contacts with the clinical outcome and prognosis of ESCC patients for potential prognostic biomarker selection. ESTIMATE, as a common method for calculation of immune and stromal scores, can bring scientific evidence to further analysis. We first obtained 81 ESCC samples with clinical data from the TCGA database. Then, using the ESTIMATE algorithm, we calculated immune and stromal scores of these patients. The correlation between these scores and clinical characteristics were also analyzed. As a result, there was a significance between immune and stromal scores and tumor grades. It is suggested that the tumor immune microenvironment in ESCC had a potential influence on tumor differentiation. Besides, identifying the biomarkers related to TME may predict and even improve the prognosis of ESCC patients.

Continuously, 299 TME-related DEGs were obtained between the high and low immune/stromal score groups. GO annotation results showed that these genes were not enriched in immune-related signaling but almost be relative to muscle function and structure, which indicated a new clue for ESCC development. Skeletal muscles contain resident immune cells and there is a cross-talk between muscle and innate immune cells in physiological and pathogenic conditions, including inflammatory myopathies, endotoxemia, or different types of muscle injury/insult [19]. Paracrine/autocrine and contact interactions have been proven to be involved in these pathological events [19]. In addition, innate immune receptors such as toll-like receptors and NOD-like receptors have influences on skeletal muscle metabolism and the muscle cells have the ability to secrete factors affecting the immune system [20]. These findings showed the correlation between immune response and muscle physiological effect, but there was little research followed with interest of tumor genesis. In our study, 299 TME-related DEGs were mainly involved in muscle system process, muscle contraction, and striated muscle cell differentiation in biological process enrichment analysis. And contractile fiber part, contractile fiber, and myofibril which mainly related to muscle structure were most enriched in cellular component. In molecular function enrichment, substrate-specific channel activity, channel activity and passive transmembrane transporter activity were 3 most significant signal in selected DEGs. In the enrichment analysis of KEGG pathways, we found that only six pathways had significantly statistics including neuroactive ligand-receptor interaction, vascular smooth muscle contraction, pancreatic secretion, cGMP-PKG signaling pathway, insulin secretion and staphylococcus aureus infection. These results indicated that the TME-related genes were involved in not only tumor immune microenvironment but also other undiscovered relative signal pathways.

Based on univariate and multivariate Cox regression analyses, three prognostic genes were identified and used to establish a risk model for predicting the prognosis of ESCC patients. The AUC value of 3-year survival was infinitely close to 1, which indicated a strong capability for predicting survival in ESCC patients. Among these three genes, COL9A3 was identified as an independent prognosis factor in ESCC. And its expression was positively correlated with the clinical outcome, that is, the patients with high expression level of COL9A3 has longer survival time than low expression group (p = 0.003).

COL9A3 encodes the major collagen component of hyaline cartilage, which is one of the three alpha chains of type IX collagen. Type IX collagen, a heterotrimeric molecule, was usually found in tissues containing type II collagen, a fibrillar collagen [21]. Mutations in this gene were usually found in the patients with multiple epiphyseal dysplasia type 3. Previous study has proven that the allelic variants in the collagen IX gene-COL9A3 was a genetic risk factor for intervertebral disc disease [22], and two single nucleotide polymorphisms introducing in COL9A3 were linked to an increased risk of lumbar disc disease[23]. In addition, in X-linked adrenoleukodystrophy patients, the combination of methylation levels of SPG20, UNC45A, and COL9A3 and also the expression levels of ID4 and MYRF would be a good marker for distinguishing the discriminating childhood from adult inflammatory phenotypes [24]. As for tumor-related research, COL9A3 was identified as tumor suppressor gene in rectal cancer [25], and it was also significantly associated with the prognosis of triple-negative breast cancer as an independent prognostic signature [26]. Seldom research about the relationship between COL9A3 and ESCC was taken by now.

GFRA2 named Glial cell line-derived neurotrophic factor family receptor alpha 2. In human neuroblastoma cells and tissues, GFRA2 was upregulated. It can promote cell proliferation by interacting with the tumor suppressor PTEN in neuroblastoma [27]. Similarly, a high expression level of GFRA2 leads to PTEN inactivation and then promotes tumor cell growth and chemoresistance in pancreatic cancer [28, 29]. It is suggested that GFRA2 may have the same effects on ESCC, as another tumor type in the digestive system. But the specific relationship between GFRA2 expression and ESCC development deserves further test and verification.

VSIG4 encodes a protein that may be a negative regulator of T-cell response. It broadly expressed in placenta, lung, and 19 other tissues. Byun et al. showed that high VSIG4 expression of cancer tissue was associated with a longer disease-free interval in benign ovarian tumors [30] and hepatitis B virus-related hepatocellular carcinoma [31]. Both Xu et al. and Hu et al. found that VISG4 could be used as a prognostic factor and a potential immunotherapeutic target for glioma [32] and clear cell renal cell carcinoma[33]. Similarly, Waldera-Lupa et al. found that together with other 2 genes, VSIG4 could be a novel biomarker for supporting the diagnosis of primary central nervous system lymphomas [34]. Tumor-associated macrophage is the prominent component of lung cancer stroma and VSIG4 may play a cancer-promoting effect in lung carcinoma development [35]. In summary, specific targeting of VSIG4 may prove to be an efficacious strategy for the treatment of ESCC, but more research should be taken for further investigation.

Plasma cell was derived from small B lymphocytes after their activation and related with some important process in tumor progression. It showed that tumor-associated plasma cell signatures emerged as a significant signal of survival for diverse solid tumors, but its infiltrated levels was associated with poor prognosis of patients both in breast and lung adenocarcinomas [36]. In triple-negative breast cancers, the infiltration level of plasma cells was highly connected with the disease recurrence [37]. With context-dependent immune responses influenced by oncogenic drivers and the presence of inflammation, CD4 + T cells carried complex and important roles within tumor microenvironments[38]. Tumor-associated macrophages are heterogeneous with diverse functions. For example, M1 macrophages inhibit tumor growth, as M2 macrophages promote tumor growth. And their phenotype and functions are regulated by the surrounding micro-environment especially TME. Due to the key roles in tumor progression, cell invasion, and metastasis [39], direct targeting tumor-related macrophages may be a potential therapy strategy for patients. These results indicated that the infiltration level may have potential significance in ESCC. Using spearman rank analysis, we found that the risk score calculated by risk model was negatively correlated with the proportions of plasma cells and positively correlated with the proportions of activated CD4 memory T cells, M1 Macrophages and M2 Macrophages. It is suggested that ESCC patients with high infiltration level in activated CD4 memory T cells, M1 Macrophages and M2 Macrophages need more attention in clinical therapy. In contrast, patients with high infiltration level in plasma cells may have better prognosis and more survival time than other patients.

Finally, GSEA was performed and confirmed the close relationship between the risk scores and immune pathways. As it shown that the pathways significantly enriched in the high-risk group involved immune response and immune system process, suggesting that immunosuppression exists in high-risk ESCC patients and these high-risk patients may have a poor outcome due to un-worked immune response.

In conclusion, our study provided a comprehensive understanding of the TME and identified a list of TME-related prognostic genes for ESCC patients. The establishment of the risk model is valuable for the early identification of high-risk patients to facilitate individualized treatment and improve the possibility of immunotherapy response.

ESCC: esophageal squamous cell carcinoma;

PD-1: programmed cell death protein 1;

PD-L1: programmed cell death-ligand 1;

TIICs: tumor-infiltrating immune cells;

TME: tumor microenvironment;

ESTIMATE: Estimation of Stromal and Immune cells in Malignant Tumor tissues using Expression data;

DEGs: differentially expressed genes;

GO: Gene Ontology;

KEGG: Kyoto Encyclopedia of Genes and Genomes;

ROC: receiver operating characteristic curve;

AUC: area under the curve;

GSEA: Gene set enrichment analysis.

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

Availability of data and materials

The datasets generated and/or analysed during the current study are available in the TCGA data portal repository (https://tcga-data.nci.nih.gov/tcga/).

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

ML conceptualized and designed the experiments. FL performed the experiments. BL, YJ and DZ analyzed the results and prepared the manuscript. HZ, YW and ZX supervised this study. All authors read and approved the final manuscript.

Acknowledgements

Not applicable.

Alsop BR, Sharma P. Esophageal Cancer. Gastroenterol Clin North Am. 2016;45(3):399–412.
Ai D, Zhu H, Ren W, Chen Y, Liu Q, Deng J, Ye J, Fan J, Zhao K. Patterns of distant organ metastases in esophageal cancer: a population-based study. J Thorac Dis. 2017;9(9):3023–30.
Hosseinzadeh S, Alipanah-Moghadam R, Isapanah Amlashi F, Nemati A. Evaluation of Haptoglobin Genotype and Some Risk Factors of Cancer in Patients with Early Stage Esophageal Cancer. Asian Pac J Cancer Prev. 2019;20(10):2897–901.
Network CGAR, University AWGA, Agency BC, Hospital BaWs, Institute B, University B, University CWR, Institute D-FC, University D, Centre GPC, et al. Integrated genomic characterization of oesophageal carcinoma. Nature. 2017;541(7636):169–75.
Brahmer JR, Tykodi SS, Chow LQ, Hwu WJ, Topalian SL, Hwu P, Drake CG, Camacho LH, Kauh J, Odunsi K, et al. Safety and activity of anti-PD-L1 antibody in patients with advanced cancer. N Engl J Med. 2012;366(26):2455–65.
Ghahremanloo A, Soltani A, Modaresi SMS, Hashemy SI. Recent advances in the clinical development of immune checkpoint blockade therapy. Cell Oncol (Dordr). 2019;42(5):609–26.
Yoshida J, Ishikawa T, Doi T, Ota T, Yasuda T, Okayama T, Sakamoto N, Inoue K, Dohi O, Yoshida N, et al. Clinical significance of soluble forms of immune checkpoint molecules in advanced esophageal cancer. Med Oncol. 2019;36(7):60.
Keren L, Bosse M, Marquez D, Angoshtari R, Jain S, Varma S, Yang SR, Kurian A, Van Valen D, West R, et al. A Structured Tumor-Immune Microenvironment in Triple Negative Breast Cancer Revealed by Multiplexed Ion Beam Imaging. Cell. 2018;174(6):1373–87.e1319.
Kelly RJ. The emerging role of immunotherapy for esophageal cancer. Curr Opin Gastroenterol 2019.
Zhao Q, Yu J, Meng X. A good start of immunotherapy in esophageal cancer. Cancer Med. 2019;8(10):4519–26.
Huang TX, Fu L. The immune landscape of esophageal cancer. Cancer Commun (Lond). 2019;39(1):79.
Yoshihara K, Shahmoradgoli M, Martínez E, Vegesna R, Kim H, Torres-Garcia W, Treviño V, Shen H, Laird PW, Levine DA, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612.
Xu M, Li Y, Li W, Zhao Q, Zhang Q, Le K, Huang Z, Yi P. Immune and Stroma Related Genes in Breast Cancer: A Comprehensive Analysis of Tumor Microenvironment Based on the Cancer Genome Atlas (TCGA) Database. Front Med (Lausanne). 2020;7:64.
Qu Y, Cheng B, Shao N, Jia Y, Song Q, Tan B, Wang J. Prognostic value of immune-related genes in the tumor microenvironment of lung adenocarcinoma and lung squamous cell carcinoma. Aging. 2020;12(6):4757–77.
Jia D, Li S, Li D, Xue H, Yang D, Liu Y. Mining TCGA database for genes of prognostic value in glioblastoma microenvironment. Aging. 2018;10(4):592–605.
Färkkilä A, Gulhan DC, Casado J, Jacobson CA, Nguyen H, Kochupurakkal B, Maliga Z, Yapp C, Chen YA, Schapiro D, et al. Immunogenomic profiling determines responses to combined PARP and PD-1 inhibition in ovarian cancer. Nat Commun. 2020;11(1):1459.
Wang S, Lin Y, Xiong X, Wang L, Guo Y, Chen Y, Chen S, Wang G, Lin P, Chen H, et al: Low-dose metformin reprograms the tumor immune microenvironment in human esophageal cancer: Results of a phase II clinical trial. Clin Cancer Res 2020.
Strizova Z, Snajdauf M, Stakheev D, Taborska P, Vachtenheim J Jr, Biskup J, Lischke R, Bartunkova J, Smrz D. The paratumoral immune cell signature reveals the potential for the implementation of immunotherapy in esophageal carcinoma patients. J Cancer Res Clin Oncol. 2020;146(8):1979–92.
Pillon NJ, Bilan PJ, Fink LN, Klip A. Cross-talk between skeletal muscle and immune cells: muscle-derived mediators and metabolic implications. Am J Physiol Endocrinol Metab. 2013;304(5):E453–65.
Pillon NJ, Krook A. Innate immune receptors in skeletal muscle metabolism. Exp Cell Res. 2017;360(1):47–54.
Landis WJ, Hodgens KJ. Visualization of sulfur-containing components associated with proliferating chondrocytes from rat epiphyseal growth plate cartilage: possible proteoglycan and collagen co-migration. Anat Rec. 1990;226(2):153–67.
Janeczko Ł, Janeczko M, Chrzanowski R, Zieliński G. The role of polymorphisms of genes encoding collagen IX and XI in lumbar disc disease. Neurol Neurochir Pol. 2014;48(1):60–2.
Zhu Y, Wu JJ, Weis MA, Mirza SK, Eyre DR. Type IX collagen neo-deposition in degenerative discs of surgical patients whether genotyped plus or minus for COL9 risk alleles. Spine. 2011;36(24):2031–8.
Decramer S, Bascands JL, Schanstra JP. Non-invasive markers of ureteropelvic junction obstruction. World J Urol. 2007;25(5):457–65.
Ge C, Wu M, Chen G, Yu G, Ji D, Wang S. Identification of molecular characteristics induced by radiotherapy in rectal cancer based on microarray data. Oncol Lett. 2017;13(4):2777–83.
Lv X, He M, Zhao Y, Zhang L, Zhu W, Jiang L, Yan Y, Fan Y, Zhao H, Zhou S, et al. Identification of potential key genes and pathways predicting pathogenesis and prognosis for triple-negative breast cancer. Cancer Cell Int. 2019;19:172.
Li Z, Xie J, Fei Y, Gao P, Xie Q, Gao W, Xu Z. GDNF family receptor alpha 2 promotes neuroblastoma cell proliferation by interacting with PTEN. Biochem Biophys Res Commun, 510(3):339–344.
Gu J, Wang D, Zhang J, Zhu Y, Li Y, Chen H, Shi M, Wang X, Shen B, Deng X, et al: GFRα2 prompts cell growth and chemoresistance through down-regulating tumor suppressor gene PTEN via Mir-17-5p in pancreatic cancer. Cancer Lett, 380(2):434–441.
Mathioudakis N, Sundaresh R, Larsen A, Ruff W, Schiller J, Guerrero-Cázares H, Burger P, Salvatori R, Quiñones-Hinojosa A. Expression of the pituitary stem/progenitor marker GFRα2 in human pituitary adenomas and normal pituitary. Pituitary. 2015;18(1):31–41.
Byun JM, Jeong DH, Choi IH, Lee DS, Kang MS, Jung KO, Jeon YK, Kim YN, Jung EJ, Lee KB, et al. The Significance of VSIG4 Expression in Ovarian Cancer. Int J Gynecol Cancer. 2017;27(5):872–8.
Zhu S, Tan W, Li W, Zhou R, Wu X, Chen X, Li W, Shang C, Chen Y. Low expression of VSIG4 is associated with poor prognosis in hepatocellular carcinoma patients with hepatitis B infection. Cancer Manag Res. 2018;10:3697–705.
Xu T, Jiang Y, Yan Y, Wang H, Lu C, Xu H, Li W, Fu D, Lu Y, Chen J. VSIG4 is highly expressed and correlated with poor prognosis of high-grade glioma patients. Am J Transl Res. 2015;7(6):1172–80.
Hu D, Zhou M, Zhu X. Deciphering Immune-Associated Genes to Predict Survival in Clear Cell Renal Cell Cancer. Biomed Res Int. 2019;2019:2506843.
Waldera-Lupa DM, Poschmann G, Kirchgaessler N, Etemad-Parishanzadeh O, Baberg F, Brocksieper M, Seidel S, Kowalski T, Brunn A, Haghikia A, et al: A Multiplex Assay for the Stratification of Patients with Primary Central Nervous System Lymphoma Using Targeted Mass Spectrometry. Cancers (Basel) 2020, 12(7).
Liao Y, Guo S, Chen Y, Cao D, Xu H, Yang C, Fei L, Ni B, Ruan Z. VSIG4 expression on macrophages facilitates lung cancer development. Lab Invest. 2014;94(7):706–15.
Gentles AJ, Newman AM, Liu CL, Bratman SV, Feng W, Kim D, Nair VS, Xu Y, Khuong A, Hoang CD, et al. The prognostic landscape of genes and infiltrating immune cells across human cancers. Nat Med. 2015;21(8):938–45.
Yeong J, Lim JCT, Lee B, Li H, Chia N, Ong CCH, Lye WK, Putti TC, Dent R, Lim E, et al. High Densities of Tumor-Associated Plasma Cells Predict Improved Prognosis in Triple Negative Breast Cancer. Front Immunol. 2018;9:1209.
Caza T, Landas S. Functional and Phenotypic Plasticity of CD4(+) T Cell Subsets. Biomed Res Int. 2015;2015:521957.
Shapouri-Moghaddam A, Mohammadian S, Vazini H, Taghadosi M, Esmaeili SA, Mardani F, Seifi B, Mohammadi A, Afshari JT, Sahebkar A. Macrophage plasticity, polarization, and function in health and disease. J Cell Physiol. 2018;233(9):6425–40.

Download PDF

Journal Publication

published 10 Mar, 2021

Read the published version in BMC Medical Genomics →

Editorial decision: Major revision
30 Dec, 2020
Review #1 received at journal
20 Dec, 2020
Review #2 received at journal
20 Dec, 2020
Reviewer #2 agreed at journal
06 Dec, 2020
Reviewers invited by journal
06 Dec, 2020
Reviewer #1 agreed at journal
06 Dec, 2020
Editor assigned by journal
02 Dec, 2020
Submission checks completed at journal
02 Dec, 2020
Editor invited by journal
02 Dec, 2020
First submitted to journal
29 Nov, 2020

You are reading this latest preprint version

Profiles of Immune Cell Infiltration and Immune-Related Genes in the Tumor Microenvironment of Esophageal Squamous Cell Carcinoma

Status:

Journal Publication

Version 1

Abstract

Figures

Background

Methods

Dataset and estimation of stromal and immune scores

Identification of DEGs based on immune and stromal scores

Functional enrichment analysis

Establishment of TME-related risk model and survival analysis

Association between risk scores and the levels of TIICs

Results

Identification of differentially expressed genes between high and low immune/stromal scores

Functional enrichment analysis

Establishment of the risk prediction model

The relationship between immune scores and TIICs

The involvement of immune pathways predicted by risk model

Discussion

Conclusion

Abbreviations

Declarations

References

Supplementary Files

Status:

Journal Publication

Version 1