Identification of Epithelial-Mesenchymal Transition- (EMT-) Related LncRNA for Prognostic Prediction and Risk Stratification in Esophageal Squamous Cell Carcinoma

Wang, Peipei; Chen, Yueyun; Zheng, Yue; Fu, Yang; Ding, Zhenyu

doi:https://doi.org/10.1155/2021/5340240

Disease Markers

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Abbreviations Data Availability Conflicts of Interest Acknowledgments Supplementary Materials References Copyright Related Articles

Special Issue

The Discovery of Nucleotide and Nucleic Acid-Associated Biomarkers and Cancer Personalized Therapies

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 5340240 | https://doi.org/10.1155/2021/5340240

Identification of Epithelial-Mesenchymal Transition- (EMT-) Related LncRNA for Prognostic Prediction and Risk Stratification in Esophageal Squamous Cell Carcinoma

Peipei Wang,¹Yueyun Chen,¹Yue Zheng,¹Yang Fu,¹and Zhenyu Ding¹

Academic Editor: Zhongqiu Xie

Received18 Jul 2021

Revised25 Sept 2021

Accepted29 Sept 2021

Published19 Oct 2021

Abstract

Background. Epithelial-mesenchymal transition (EMT) is significantly associated with the invasion and development of esophageal squamous cell carcinoma (ESCC). However, the importance of EMT-related long noncoding RNA (lncRNA) is little known in ESCC. Methods. GSE53624 () and GSE53622 () datasets retrieved from the Gene Expression Omnibus (GEO) database were used as training and external validation cohorts, respectively. GSE53624 and GSE53622 datasets were all sampled from China. Then, the prognostic value of EMT-related lncRNA was comprehensively investigated by weighted coexpression network analysis (WGCNA) and COX regression model. Results. High expression of PLA2G4E-AS1, AC063976.1, and LINC01592 significantly correlated with the favorable overall survival (OS) of ESCC patients, and LINC01592 had the greatest contribution to OS. Importantly, ESCC patients were divided into low- and high-risk groups based on the optimal cut-off value of risk score estimated by the multivariate COX regression model of these three lncRNA. Patients with high risk had a shorter OS rate and restricted mean survival time (RMST) than those with low risk. Moreover, univariate and multivariate COX regression revealed that risk stratification, age, and TNM were independent prognostic predictors, which were used to construct a nomogram model for individualized and visualized prognosis prediction of ESCC patients. The calibration curves and time-dependent ROC curves in the training and validation cohorts suggested that the nomogram model had a good performance. Interestingly, clear trends indicated that risk score positively correlated with tumor microenvironment (TME) scores and immune checkpoints TIGIT, CTLA4, and BTLA. In addition, the Kyoto Encyclopedia of Genes and Genomes (KEGG) showed that PLA2G4E-AS1, AC063976.1, and LINC01592 were primarily associated with TNF signaling pathway, NF-kappa B signaling pathway, and ECM-receptor interaction. Conclusion. We developed EMT-related lncRNA PLA2G4E-AS1, AC063976.1, and LINC01592 for prognostic prediction and risk stratification of Chinese ESCC patients, which might provide deep insight for personalized prognosis prediction in Chinese ESCC patients and be potential biomarkers for designing novel therapy.

1. Introduction

Esophageal cancer is one of the top ten malignant tumors globally and the sixth leading cause of cancer-related deaths. Esophageal cancer is mainly prevalent in eastern Asia and eastern and southern Africa [1, 2]. Esophageal squamous cell carcinoma (ESCC) is the most common histological subtype of esophageal cancer, accounting for about 90%. ESCC has a high tendency to be aggressive and metastatic, as well as a high chance of recurrence. Even with multimodality treatments (surgery, radiotherapy, chemotherapy, and targeted therapy), the prognosis of ESCC patients remains poor [3]. Insight into the molecules and mechanisms behind ESCC invasion and metastasis helps deepen the understanding of the disease. It is also urgent to discover novel biomarkers to develop new therapeutic strategies and improve the prognosis of ESCC patients.

Long noncoding RNA (LncRNA) is a group of evolutionarily conserved RNA molecules with more than 200 nucleotides in length, lacking protein-coding ability [4]. Abnormal expression of LncRNA plays a vital role in the occurrence and development of various tumors [5–9], including ESCC [10, 11]. Many studies have shown that lncRNA plays multiple roles in malignant behaviors such as tumor formation, invasion, migration, and immunogenicity. LncRNA BRCAT54 inhibits the tumorigenesis of non-small -ell lung cancer by binding to RPS9 to regulate JAK-STAT and calcium pathways [12]. Moreover, lncRNA LINC00472 regulates cell stiffness and inhibits the migration and invasion of lung adenocarcinoma by binding to YBX1 [13]. In addition, a study found that LIMIT may be a target for cancer immunotherapy [14]. LncRNA can participate in the occurrence and development of tumors by affecting chromatin remodeling, histone modification, DNA methylation, gene transcription, translation, etc.

Notably, the effect of lncRNA on the epithelial-mesenchymal transition (EMT) of tumor cells has also been confirmed in recent studies [15, 16]. EMT is an essential step in the metastasis of malignant tumors, which can transform epithelial-like cells into a mesenchymal-like cell state. By modifying the adhesion molecules expressed in cells, EMT reduces the adhesion ability of epithelial-derived tumor cells, thereby causing epithelial cells to separate from each other, increasing the metastatic potential of tumor cells, and further resisting antitumor treatments [17, 18]. Accordingly, EMT plays a crucial role in the metastasis of epithelial tumors. Few studies have explored the prognostic value of EMT-related lncRNA in cancer patients [19]. The role of EMT-related lncRNA in ESCC and related mechanisms is little known.

This study comprehensively investigated and validated the prognostic value of EMT-related lncRNA in Chinese ESCC patients from the Gene Expression Omnibus (GEO) database by weighted gene coexpression network analysis (WGCNA), COX proportional hazard regression, and Kaplan-Meier survival analysis. Furthermore, risk stratification and a nomogram model were constructed to personalize and visualize the overall survival (OS) rates of ESCC patients. Additionally, the correlation between weighted EMT-related lncRNA with tumor microenvironment (TME) scores, immune checkpoints (ICs), and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways was further explored.

2. Materials and Methods

2.1. ESCC Patients

The transcriptome data of 119 and 60 newly diagnosed ESCC patients in the GSE53624 and GSE53622 datasetswere downloaded from the GEO database (https://www.ncbi.nlm.nih.gov/geo/), assigned as training and validation cohorts, respectively. GSE53624 and GSE53622 datasets were all sampled from China. The clinical characteristics, including age, gender, tumor invasion depth (T), lymph node metastasis (N), tumor node metastasis (TNM) staging, tumor grade, overall survival (OS) time, and events, were obtained and listed in Table S1. The workflow of data analysis in this study was performed according to Figure 1. Since the GEO database is publicly available, no approval from the local ethics committee was required.

Figure 1

Schematic diagram of the study. The transcriptome sequencing data and clinical information of patients with esophageal squamous cell carcinoma (ESCC) were downloaded from the Gene Expression Omnibus (GEO) database. Different datasets were assigned as training and validation cohorts for weighted coexpression network analysis (WGCNA) and univariate COX proportional hazards regression analysis. Then, the epithelial-mesenchymal transition- (EMT-) related lncRNA determined by multivariate COX regression analysis was used for risk stratification and constructing a nomogram model. Finally, the pathways of EMT-related lncRNA and the correlation between lncRNA and tumor microenvironment (TME) or immune checkpoints (ICs) were explored.

2.2. Acquisition of lncRNA and EMT-Related mRNA

A total of 17,936 lncRNA (version 35) were downloaded from the GENCODE database (https://www.gencodegenes.org/human/) [20]. Furthermore, the 50 hallmark gene sets and the “HALLMARK_EPITHELIAL_MESENCHYMAL_ TRANSITION” gene list, including 200 EMT-related mRNA, were obtained from the Gene Set Enrichment Analysis (GSEA) database (http://www.gsea-msigdb.org/gsea/msigdb/collections.jsp#H) [21, 22].

2.3. Weighted Gene Coexpression Network Analysis (WGCNA)

As an unsupervised machine learning, WGCNA was applied for investigating the correlation between genes. The package “WGCNA” in software (version 4.0.2, https://www.r-project.org/) was used to construct a weighted coexpression network between the lncRNA and EMT-related mRNA [23]. In the network, the pairwise Pearson coefficient was used to evaluate the coexpression weight among all genes. The power β of the soft threshold was used to confirm a scale-free network. Notably, the genes with similar expression patterns were clustered into the same color module in the unsupervised coexpression network.

2.4. Nomogram Model

The packages “foreign” and “rms” in software (version 4.0.2, https://www.r-project.org/) were used to construct a nomogram model to personalize and visualize the OS rate of ESCC patients [24, 25]. Each variable was assigned a point according to the nomogram model. Then, the total points were obtained by summing the points of all variables for determining the OS rate of an ESCC patient. Finally, time-dependent receiver operating characteristic (ROC) and calibration curves of the training and validation cohort were used to evaluate the accuracy of the nomogram model that predicted the OS rate.

2.5. Estimation of Tumor Microenvironment (TME) Score

The ESTIMATE algorithm was used to calculate the fraction of immune and stromal cells in ESCC tissues based on gene expression levels [26]. The ESTIMATE algorithm was performed using package “estimate” to calculate TME, immune, and stromal scores in each ESCC patient.

2.6. Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathways

“DOSE,” “http://org.Hs.eg.db,” “topGO,” and “clusterProfiler” packages in software (version 4.0.2, https://www.r-project.org/) were used to obtain the KEGG pathways of EMT-related lncRNA in ESCC patients.

2.7. Statistical Analysis

All statistical analysis was performed using software (version 4.0.2, https://www.r-project.org/). A chi-square and Fisher tests were used to compare differences between two groups of categorical variables, as appropriate. Univariate and multivariate Cox proportional hazard regression analysis was performed using package “survival.” The “surv_cutpoint” function in the package “survminer” was used to determine the optimal cut-point of a gene (Figure S1). Kaplan-Meier curves were compared by the log-rank test. The package “survRM2” was used to obtain the restricted mean survival time (RMST). Correlation coefficients between two quantitative variables were obtained by Pearson’s method. The “survivalROC” package determined the area under the curve (AUC) in the time-dependent ROC curve. A two-tailed value <0.05 and a value <0.1 were considered statistically significant and a clear trend, respectively.

3. Results

3.1. WGCNA for lncRNA and EMT-Related mRNA in ESCC Patients

As shown in Table S1, the clinical characteristics were balanced between GSE53624 and GSE53622 datasets (). To identify EMT-related lncRNA in ESCC patients, 188 EMT-related mRNA and 5,506 lncRNA were included in the construction of WGCNA, and the workflow for data analysis was shown in Figure 1. The soft thresholds for building a scale-free network of training and validation cohorts were set to 4 and 5, respectively (Figures 2(a) and 2(b)). Then, a total of 3,900 lncRNA were coexpressed with EMT-related mRNA in the training cohort, which was distributed in 6 modules, including black, blue, brown, green, grey, and turquoise modules (Figure 2(a)). Moreover, 3,742 lncRNA and EMT-related mRNA were showed in 11 coexpression modules, including black, blue, green, greenyellow, grey, magenta, pink, purple, red, turquoise, and yellow, in the validation cohort (Figure 2(b)). Therefore, 3,900 and 3,742 EMT-related lncRNA in the training and validation cohorts, respectively, were used for the following univariate and multivariate COX regression analysis. Notably, the percent of overlapping EMT-related lncRNA between training and validation cohorts was 78.1% (3,045/3,900) and 81.4% (3,045/3,742), respectively.

(a)

(b)

3.2. Univariate and Multivariate COX Regression Analysis

After univariate COX regression analysis, 338 and 169 EMT-related lncRNA were significantly associated with the OS of ESCC patients in the training and validation cohorts, respectively (, Figure 3(a)). To further confirm that EMT-related lncRNA had prognostic value in both cohorts, lncRNA with a or HR<1 in the univariate regression model was overlapped between the training and validation cohorts, and the results showed that 3 EMT-related lncRNA expression, including PLA2G4E-AS1, AC063976.1, and LINC01592, significantly correlated with the favorable OS of ESCC patients (, Figure 3(b)). Then, multivariate COX regression analysis was used for weighted combination of AC063976.1, LINC01592, and PLA2G4E-AS1, which indicated that LINC01592 contributed the greatest to the OS of ESCC patients (). Notably, the time-dependent ROC curve results demonstrated that the multivariate COX regression model performs well in the training cohort (, Figure 3(c)). This finding was confirmed in the validation cohort (, Figure 3(c)).

(a)

(b)

(c)

Figure 3

Univariate and multivariate COX regression analysis of EMT-related lncRNA in the training and validation cohorts. (a) The heatmaps displayed the hazard ratio (HR) and 95% confidence interval (CI) of EMT-related lncRNA with by univariate COX regression analysis in the training (a) and validation (c) cohorts. Blue to red in the color scale represented HR from low to high. HR.95 L: the low value in the 95% confidence interval of HR; HR.95H: the high value in the 95% confidence interval of HR. (b) The overlapping number of EMT-related lncRNA with (a) or (c) in the training and validation cohorts. (c) Multivariate COX regression analysis of EMT-related lncRNA. The contribution of EMT-related lncRNA to the overall survival (OS) of ESCC patients in the training cohort (upper panel). The time-dependent receiver operator characteristic (ROC) curves of the training (b) and validation (d) cohorts were used to evaluate the performance of the multivariate COX regression model.

3.3. Establishment of Risk Stratification for ESCC Patients

To establish a risk stratification for ESCC patients, we first obtain the risk score based on the coefficients of multivariate COX regression, and the formula for calculating the risk score was as follows: (Figure 3(c)). Based on the optimal prognostic cut-point of risk score -6.57, ESCC patients were divided into high- and low-risk groups. ESCC patients with high risk were significantly associated with poor OS in the training cohort (, 95% confidence interval (CI): 2.05 to 6.66, , Figure 4(a)). This result was confirmed in the validation cohort (, 95% CI: 1.27 to 5.07, , Figure 4(c)). Furthermore, high-risk ESCC patients had a shorter RMST than low-risk patients in the training cohort (4-year RMST: 25 (95% CI: 22 to 29) vs. 40 (95% CI: 36 to 44) months) (Figure 4(a)). This result was again confirmed in the validation cohort (4-year RMST: 24 (95% CI: 17 to 31) vs. 37 (95% CI: 32 to 42) months) (Figure 4(c)). Interestingly, The Kaplan-Meier curves indicated that high expression of AC063976.1, LINC01592, and PLA2G4E-AS1 correlated with the favorable OS of ESCC patients in both the training and validation cohorts (, Figures 4(b) and 4(d)). Importantly, risk stratification was an independent prognostic predictor for ESCC patients by univariate and multivariate COX regression analysis in the training cohort (, 95% CI: 2.13 to 7.11, , Table 1). This was again confirmed in the validation cohort (, 95% CI: 1.29 to 5.78, , Table 1). In order to determine whether risk stratification has prognostic significance in a random population, Microsoft Excel 2016 was further used to randomly select a portion of samples in each dataset as a training cohort and then treat the other samples as a validation cohort. Patients with a high-risk score were associated with a poor OS in the training cohort of the GSE53624 dataset (, 95% CI: 1.21 to 5.91, ). This result was confirmed in the validation cohort of GSE53624 dataset (, 95% CI: 2.13 to 12.47, ) (Figure 5(a)). Interestingly, high-risk patients had a shorter OS than low-risk patients in the training cohort of the GSE53622 dataset; although, it is not statistically significant at that point (, 95% CI: 0.91 to 6.66, ). This finding was again confirmed in the validation cohort of GSE53622 dataset (, 95% CI: 0.96 to 6.62, ) (Figure 5(b)). This might be due to the small sample size in the GSE53622 dataset. To confirm the risk stratification that could better predict prognosis in which part of the population of ESCC patients, we conducted a subgroup analysis. There is a clear trend in the training and validation cohorts that risk stratification can better predict the prognosis in patients with TNM III/IV stage, N1-3, T3-4, male, and >60 years old. The prognosis can be predicted well through risk stratification regardless of the tumor grade (Figure 6).

(a)

(b)

3.4. Construction of a Nomogram Visualizing and Personalizing the OS Rate of ESCC Patients

Univariate and multivariate COX regression analysis was used to identify independent prognosis factors for constructing a nomogram model. In addition to risk stratification, age and TNM stage were independent prognostic factors for ESCC patients in both the training and validation cohorts (, , Table 1). Accordingly, a nomogram model constructed by the risk stratification, age, and TNM stage could visualize and personalize the 1-, 2-, 3-, and 4-year OS rates of ESCC patients (Figure 7(a)). Details of the points for the variables and OS rates in the nomogram model were listed in Table S2. The time-dependent ROC and calibration curves were further used to evaluate the predicted performance of the nomogram model. Notably, the time-dependent ROC curves illustrated that all the AUCs were ≥0.70 in both the training and validation cohorts (Figure 7(b)). Moreover, calibration curves indicated that the 1-, 2-, 3-, and 4-year OS rates predicted in the nomogram model were highly in line with actual observations in both the training and validation cohorts (Figures 7(c) and 7(d)). These results suggested that the nomogram model had good performance in predicting the OS rate of ESCC patients.

(a)

(b)

(c)

(d)

Figure 7

Construction of the nomogram model. (a) The combination of risk stratification, age, and TNM stage visualized and personalized the OS rate of ESCC patients. After a point was assigned to risk stratification, age, and TNM stage of each patient according to the nomogram model, the total points were obtained to predict the OS rate of ESCC patients. (b) The time-dependent receive operating characteristic curve (ROC) curve was used to evaluate the performance of the nomogram model in the training (a, b) and validation (c, d) cohorts. (c d) The calibration curves were used to evaluate the performance of the nomogram model in predicting the OS rate in training (c) and validation (d) cohorts. The criterion for good performance of the nomogram model was that the predicted OS rate was relatively consistent with the actual OS rate. AUC: the area under the curve.

3.5. KEGG Pathways for EMT-Related lncRNA

Based on WGCNA, 57 and 25 EMT-related genes were coexpressed with AC063976.1, LINC01592, or PLA2G4E-AS1 in the training and validation cohorts, respectively (Figure 8(a)). Then, the KEGG was applied to identify the significant pathway associated with the AC063976.1, LINC01592, or PLA2G4E-AS1 in ESCC patients. The results showed that in the training and validation cohorts, a total of 7 and 8 pathways were enriched, respectively. And three overlapped pathways, including TNF signaling pathway, NF-kappa B signaling pathway, and ECM-receptor interaction, were enriched in both cohorts (Figures 8(c) and 8(d)).

(a)

(b)

(c)

(d)

3.6. The Risk Score Was Positively Correlated with the TME Score and Immune Checkpoints (ICs)

Ample reports found that EMT-related genes were enriched in TME, whereas little was known about these genes in ESCC. This prompted us to investigate further the relationship between the risk score based on EMT-related lncRNA and TME score and the expression levels of ICs. TME score was composed of stromal score and immune score. The results suggested that risk score was positively correlated with TME score (, ). Further analysis found that the risk score had a positive correlation with immune score (, ); although, statistical significance was not reached at this point. Because ICs is closely related to immune score, we further analyzed the correlation between risk score and ICs. Notably, risk score had a significantly positive correlation with BTLA (, ) and CTLA4 (, ). What is more, there was a clear trend suggested that the risk score was positively correlated with TIGIT (, ); although, statistical significance was not reached at this point. However, there was no significant correlation between risk score and PD-1, PD-L1, PD-L2, LAG3, and CD276 (). Interestingly, risk score was also positively correlated with stromal scores (, ) (Figures 9(a) and 9(b)). We further analyzed the correlation between risk score and cancer-associated fibroblast- (CAFs-) related genes. The results suggested that the risk score had a significant positive correlation with MGP (, ), MFAP5 (, ), ITGA11 (, ), DCN (, ), or ACTA2 (, ). Moreover, there was a clear trend showing that risk score was positively correlated with COL11A1 (, ) and BMP4 (, ); although. the data were not yet significant enough at this point. However, there is no significant correlation between risk score and SPHK1, CSPG4, TGFBI, and TNNC1 () (Figures 9(a) and 9(b)).

(a)

(b)

Figure 9

The correlation between risk score estimated by nomogram and tumor microenvironment (TME) and immune checkpoints (ICs) in the training cohort. (a) Distribution of risk score, TME score, ICs, and cancer-associated fibroblast (CAF) related gene expression level. package “estimate” was used to calculate TME, immune, and stromal scores in each ESCC patient based on gene expression levels. (b) The relationship between risk score and TME, ICs (a) and CAFs (b). The size of the dot represents the size of the correlation coefficient, and the color scale from blue to yellow represents the value from large to small.

4. Discussion

ESCC is the prevalent histological subtype of esophageal cancer, with a poor prognosis and prone to distant metastasis [3, 27]. Therefore, exploring potential biomarkers is essential for managing and predicting the prognosis of ESCC patients. Increasing evidence suggests that EMT is highly correlated with cancer progression and metastasis [28]. In recent years, a few studies have investigated the role of EMT-related lncRNA in the prognosis and progression of ESCC [29, 30]. However, based on next-generation transcriptome sequencing, a comprehensive assessment of the prognostic importance, risk stratification, and visualization of OS rates by EMT-related lncRNA in ESCC was little known.

In this study, based on the analysis of two large datasets in the GEO database, the results suggested that the high expression of AC063976.1, LINC01592, or PLA2G4E-AS1 was significantly associated with favorable OS in ESCC patients. In addition, KEGG results indicated that AC063976.1, LINC01592, or PLA2G4E-AS1 were mainly enriched in the TNF signaling pathway, NF-kappa B signaling pathway, and ECM-receptor interaction. AC063976.1, as a novel EMT-related lncRNA, has not yet been explored in cancer. For the first time, we revealed that upregulation of AC063976.1 corrected with a favorable OS of ESCC patients. Li et al. reported that LINC01592 was a protective factor for ESCC patients [31], which was consistent with this study. Furthermore, LINC01592 contributed the greatest to the OS of ESCC patients. Although PLA2G4E-AS1 was downregulated in thyroid carcinoma, its prognostic importance in cancer patients has not been elucidated [32]. The results of this study demonstrated that the high expression of PLA2G4E-AS1 could predict favorable OS of ESCC patients. These findings will provide prognostic information for further exploring the functions and mechanisms of AC063976.1, LINC01592, or PLA2G4E-AS1 in ESCC in the future.

Risk stratification plays a vital role in guiding clinical treatment and management of cancer patients [33]. Notably, the risk stratification constructed by a weighted combination of AC063976.1, LINC01592, and PLA2G4E-AS1 divided ESCC patients into low- and high-risk groups, implying it was also an independent prognostic predictor for ESCC patients. Interestingly, a subgroup analysis found that risk stratification was mainly performed in ESCC patients with TNM III/IV stage, N1-3, T3-4, male, or >60 y. Furthermore, the risk stratification could be used regardless of the tumor grade. Notably, a nomogram model established by the risk stratification, age, and TNM stage could display and visualize the 1-, 2-, 3-, and 4-year OS rates of ESCC patients, which might contribute to the management of individualized treatment.

Previous studies reported that TME was associated with EMT in cancers [34]. Moreover, the stromal microenvironment and CAFs play an important role in EMT [35–38]. Hence, the relationship between risk scores calculated by AC063976.1, LINC01592, and PLA2G4E-AS1 and TME was further investigated. In this study, the risk score was positively correlated with TME, immune, and stromal scores. Furthermore, the risk score had a significant positive correlation with CAF genes, including MGP, ITGA11, DCN, ACTA2, COL11A1, or BMP4. However, high expression levels of ICs usually lead to T cell exhaustion in cancers [39, 40]. Interestingly, there was also a positive correlation between risk score and ICs, including BTLA, TIGIT, and CTLA4. These results can be interpreted that the antitumor effect of the high level of immune cell infiltration is offset by the strong immunosuppressive pathway activated by upregulated IC proteins [34, 41] and might provide the likelihood of immunotherapy for high-risk ESCC patients.

This study had several limitations: first, the results of the analysis and validation of transcriptome sequencing data in this study were based on publicly available datasets. Therefore, some important clinical information was incomplete, such as treatment options, which might produce potential biases in conclusions. Secondly, this study did not provide additional validation by original ESCC samples from our clinical center. Finally, EMT-related lncRNA should be validated through in vivo and in vitro experiments in the future.

5. Conclusions

We demonstrated that based on risk stratification constructed by PLA2G4E-AS1, AC063976.1, and LINC01592, ESCC patients were divided into low- and high-risk groups. Moreover, a nomogram model established by the risk stratification, age, and TNM stage could display and visualize the 1-, 2-, 3-, and 4-year OS rates of Chinese ESCC patients. In addition, the risk score was positively correlated with the TME score, ICs, and CAFs. These findings might provide deep insight for personalized prognosis prediction by EMT-related lncRNA in Chinese ESCC patients, and the three lncRNAs might be potential biomarkers for designing novel therapy.

Abbreviations

AUC:	Area under the curve
CI:	Confidence interval
EMT:	Epithelial-mesenchymal transition
ESCC:	Esophageal squamous cell carcinoma
GEO:	Gene expression omnibus
GSEA:	Gene set enrichment analysis
HR:	Hazard ratio
IC:	Immune checkpoint
KEGG:	Kyoto Encyclopedia of Genes and Genomes
lncRNA:	Long noncoding RNA
OS:	Overall survival
RMST:	Restricted mean survival time
ROC:	Receiver operating characteristic
TME:	Tumor microenvironment
TNM:	Tumor invasion depth, lymph node metastasis, tumor node metastasis
WGCNA:	Weighted coexpression network analysis.

Data Availability

The original contributions presented in the study are included in the article/Supplementary Material, and further inquiries can be directed to the corresponding authors.

Conflicts of Interest

The authors have declared that no competing interest exists.

Acknowledgments

This work was supported by the Clinical Research Incubation Project, West China Hospital, Sichuan University (2020HXFH056).

Supplementary Materials

Figure S1: the optimal cut-points for risk score, AC063976.1, LINC01592, and PLA2G4E-AS1 in the training (a) and validation (b) cohorts. Table S1: clinical characteristics of ESCC patients. Table S2: the points for the nomogram model. (Supplementary Materials)

References

H. Sung, J. Ferlay, R. L. Siegel et al., “Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries,” CA: a Cancer Journal for Clinicians, vol. 71, no. 3, pp. 209–249, 2021.
View at: Publisher Site | Google Scholar
J. Lagergren, E. Smyth, D. Cunningham, and P. Lagergren, “Oesophageal cancer,” Lancet, vol. 390, no. 10110, pp. 2383–2396, 2017.
View at: Publisher Site | Google Scholar
S. B. Chen, H. R. Weng, G. Wang et al., “Prognostic factors and outcome for patients with esophageal squamous cell carcinoma underwent surgical resection alone: evaluation of the seventh edition of the American joint committee on cancer staging system for esophageal squamous cell carcinoma,” Journal of Thoracic Oncology, vol. 8, no. 4, pp. 495–501, 2013.
View at: Publisher Site | Google Scholar
L. Ma, V. B. Bajic, and Z. Zhang, “On the classification of long non-coding RNAs,” RNA Biology, vol. 10, no. 6, pp. 924–933, 2013.
View at: Publisher Site | Google Scholar
C. T. Chen, P. P. Wang, W. J. Mo et al., “Expression profile analysis of prognostic long non-coding RNA in adult acute myeloid leukemia by weighted gene co-expression network analysis (WGCNA),” Journal of Cancer, vol. 10, no. 19, pp. 4707–4718, 2019.
View at: Publisher Site | Google Scholar
W. Xu, K. Li, C. Song et al., “Knockdown of lncRNA LINC01234 suppresses the tumorigenesis of liver cancer via sponging miR-513a-5p,” Frontiers in Oncology, vol. 10, article 571565, 2020.
View at: Publisher Site | Google Scholar
Y. Lv, Y. Wang, Y. Song et al., “LncRNA PINK1-AS promotes Gαi1-driven gastric cancer tumorigenesis by sponging microRNA-200a,” Oncogene, vol. 40, no. 22, pp. 3826–3844, 2021.
View at: Publisher Site | Google Scholar
K. C. Chang, S. D. Diermeier, A. T. Yu et al., “_MaTAR25_ lncRNA regulates the _Tensin1_ gene to impact breast cancer progression,” Nature Communications, vol. 11, no. 1, p. 6438, 2020.
View at: Publisher Site | Google Scholar
Q.-N. Wu, X.-J. Luo, J. Liu et al., “MYC-activated LncRNAMNX1-AS1Promotes the progression of colorectal cancer by stabilizing YB1,” Cancer Research, vol. 81, no. 10, pp. 2636–2650, 2021.
View at: Publisher Site | Google Scholar
C. Zhang, L. Xie, Y. Fu, J. Yang, and Y. Cui, “lncRNA MIAT promotes esophageal squamous cell carcinoma progression by regulating miR-1301-3p/INCENP axis and interacting with SOX2,” Journal of Cellular Physiology, vol. 235, no. 11, pp. 7933–7944, 2020.
View at: Publisher Site | Google Scholar
X. Li, Z. Wu, Q. Mei et al., “Long non-coding RNA HOTAIR, a driver of malignancy, predicts negative prognosis and exhibits oncogenic activity in oesophageal squamous cell carcinoma,” British Journal of Cancer, vol. 109, no. 8, pp. 2266–2278, 2013.
View at: Publisher Site | Google Scholar
W. Yang, Y. Qian, K. Gao et al., “LncRNA BRCAT54 inhibits the tumorigenesis of non-small cell lung cancer by binding to RPS9 to transcriptionally regulate JAK-STAT and calcium pathway genes,” Carcinogenesis, vol. 42, no. 1, pp. 80–92, 2021.
View at: Publisher Site | Google Scholar
X. Deng, W. Xiong, X. Jiang et al., “LncRNA LINC00472 regulates cell stiffness and inhibits the migration and invasion of lung adenocarcinoma by binding to YBX1,” Cell Death & Disease, vol. 11, no. 11, p. 945, 2020.
View at: Publisher Site | Google Scholar
G. Li, I. Kryczek, J. Nam et al., “LIMIT is an immunogenic lncRNA in cancer immunity and immunotherapy,” Nature Cell Biology, vol. 23, no. 5, pp. 526–537, 2021.
View at: Publisher Site | Google Scholar
J. Jarroux, D. Foretek, C. Bertrand et al., “HOTAIR lncRNA promotes epithelial-mesenchymal transition by redistributing LSD1 at regulatory chromatin regions,” EMBO Reports, vol. 22, no. 7, article e50193, 2021.
View at: Publisher Site | Google Scholar
P. Wang, W. Chen, T. Ma et al., “lncRNA lnc-TSI inhibits metastasis of clear cell renal cell carcinoma by suppressing TGF-β-induced epithelial-mesenchymal transition,” Mol Ther Nucleic Acids., vol. 22, pp. 1–16, 2020.
View at: Publisher Site | Google Scholar
M. Taki, K. Abiko, M. Ukita et al., “Tumor immune microenvironment during epithelial-Mesenchymal transition,” Clinical Cancer Research, vol. 27, no. 17, pp. 4669–4679, 2021.
View at: Publisher Site | Google Scholar
N. Erin, J. Grahovac, A. Brozovic, and T. Efferth, “Tumor microenvironment and epithelial mesenchymal transition as targets to overcome tumor multidrug resistance,” Drug Resistance Updates, vol. 53, article 100715, 2020.
View at: Publisher Site | Google Scholar
C. Liu, C. Hu, J. Li, L. Jiang, and C. Zhao, “Identification of epithelial-mesenchymal transition-related lncRNAs that associated with the prognosis and immune microenvironment in colorectal cancer,” Frontiers in Molecular Biosciences, vol. 8, article 633951, 2021.
View at: Publisher Site | Google Scholar
A. Frankish, M. Diekhans, A. M. Ferreira et al., “GENCODE reference annotation for the human and mouse genomes,” Nucleic Acids Research, vol. 47, no. D1, pp. D766–d773, 2019.
View at: Publisher Site | Google Scholar
A. Subramanian, P. Tamayo, V. K. Mootha et al., “Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 43, pp. 15545–15550, 2005.
View at: Publisher Site | Google Scholar
V. K. Mootha, C. M. Lindgren, K. F. Eriksson et al., “PGC-1α-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes,” Nature Genetics, vol. 34, no. 3, pp. 267–273, 2003.
View at: Publisher Site | Google Scholar
P. Langfelder and S. Horvath, “WGCNA: an R package for weighted correlation network analysis,” BMC Bioinformatics, vol. 9, no. 1, p. 559, 2008.
View at: Publisher Site | Google Scholar
D. Dong, M. J. Fang, L. Tang et al., “Deep learning radiomic nomogram can predict the number of lymph node metastasis in locally advanced gastric cancer: an international multicenter study,” Annals of Oncology, vol. 31, no. 7, pp. 912–920, 2020.
View at: Publisher Site | Google Scholar
P. Wang, Y. Fu, Y. Chen et al., “Nomogram personalizes and visualizes the overall survival of patients with triple-negative breast cancer based on the immune genome,” BioMed Research International, vol. 2020, Article ID 4029062, 2020.
View at: Publisher Site | Google Scholar
K. Yoshihara, M. Shahmoradgoli, E. Martínez et al., “Inferring tumour purity and stromal and immune cell admixture from expression data,” Nature Communications, vol. 4, no. 1, p. 2612, 2013.
View at: Publisher Site | Google Scholar
D. Zhang, C. Qian, H. Wei, and X. Qian, “Identification of the prognostic value of tumor microenvironment-related genes in esophageal squamous cell carcinoma,” Frontiers in Molecular Biosciences, vol. 7, article 599475, 2020.
View at: Publisher Site | Google Scholar
Y. Tian, P. Qi, Q. Niu, and X. Hu, “Combined snail and E-cadherin predicts overall survival of cervical carcinoma patients: comparison among various epithelial-mesenchymal transition proteins,” Frontiers in Molecular Biosciences, vol. 7, p. 22, 2020.
View at: Publisher Site | Google Scholar
G. D. Gao, X. Y. Liu, Y. Lin, H. F. Liu, and G. J. Zhang, “LncRNA CASC9 promotes tumorigenesis by affecting EMT and predicts poor prognosis in esophageal squamous cell cancer,” European Review for Medical and Pharmacological Sciences, vol. 22, no. 2, pp. 422–429, 2018.
View at: Publisher Site | Google Scholar
P. Ren, H. Zhang, L. Chang, X. D. Hong, and L. Xing, “LncRNA NR2F1-AS1 promotes proliferation and metastasis of ESCC cells via regulating EMT,” European Review for Medical and Pharmacological Sciences, vol. 24, no. 7, pp. 3686–3693, 2020.
View at: Publisher Site | Google Scholar
W. Li, J. Liu, and H. Zhao, “Identification of a nomogram based on long non-coding RNA to improve prognosis prediction of esophageal squamous cell carcinoma,” Aging (Albany NY), vol. 12, no. 2, pp. 1512–1526, 2020.
View at: Publisher Site | Google Scholar
Y. Zhang, T. Jin, H. Shen, J. Yan, M. Guan, and X. Jin, “Identification of long non-coding RNA expression profiles and co-expression genes in thyroid carcinoma based on the Cancer Genome Atlas (TCGA) Database,” Medical Science Monitor, vol. 25, pp. 9752–9769, 2019.
View at: Publisher Site | Google Scholar
M. S. Tallman, E. S. Wang, J. K. Altman et al., “Acute myeloid leukemia, version 3.2019, NCCN clinical practice guidelines in oncology,” Journal of the National Comprehensive Cancer Network, vol. 17, no. 6, pp. 721–749, 2019.
View at: Publisher Site | Google Scholar
W. Zhong, F. Zhang, C. Huang, Y. Lin, and J. Huang, “Identification of epithelial-mesenchymal transition-related lncRNA with prognosis and molecular subtypes in clear cell renal cell carcinoma,” Frontiers in Oncology, vol. 10, article 591254, 2020.
View at: Publisher Site | Google Scholar
X. Wu, L. Ruan, Y. Yang, and Q. Mei, “Analysis of gene expression changes associated with human carcinoma-associated fibroblasts in non-small cell lung carcinoma,” Biological Research, vol. 50, no. 1, p. 6, 2017.
View at: Publisher Site | Google Scholar
S. Wen, N. Qu, B. Ma et al., “Cancer-associated fibroblasts positively correlate with dedifferentiation and aggressiveness of thyroid cancer,” Oncotargets and Therapy, vol. Volume 14, pp. 1205–1217, 2021.
View at: Publisher Site | Google Scholar
M. Nurmik, P. Ullmann, F. Rodriguez, S. Haan, and E. Letellier, “In search of definitions: cancer-associated fibroblasts and their markers,” International Journal of Cancer, vol. 146, no. 4, pp. 895–905, 2020.
View at: Publisher Site | Google Scholar
A. Calon, E. Lonardo, A. Berenguer-Llergo et al., “Stromal gene expression defines poor-prognosis subtypes in colorectal cancer,” Nature Genetics, vol. 47, no. 4, pp. 320–329, 2015.
View at: Publisher Site | Google Scholar
C. Chen, C. Liang, S. Wang et al., “Expression patterns of immune checkpoints in acute myeloid leukemia,” Journal of Hematology & Oncology, vol. 13, no. 1, p. 28, 2020.
View at: Publisher Site | Google Scholar
C. Chen, L. Xu, R. Gao et al., “Transcriptome-based co-expression of BRD4 and PD-1/PD-L1 predicts poor overall survival in patients with acute myeloid leukemia,” Frontiers in Pharmacology, vol. 11, article 582955, 2021.
View at: Publisher Site | Google Scholar
Y. Yao, Z. C. Wang, D. Yu, and Z. Liu, “Role of allergen-specific T-follicular helper cells in immunotherapy,” Current Opinion in Allergy and Clinical Immunology, vol. 18, no. 6, pp. 495–501, 2018.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Peipei Wang 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.

PDF Download Citation

Download other formats

Order printed copies

Views

375

Downloads

565

Citations