Skip to main content

An immune-related lncRNA pairs signature to identify the prognosis and predict the immune landscape of laryngeal squamous cell carcinoma

Abstract

Background

Laryngeal squamous cell carcinoma (LSCC) is the most common squamous cell carcinoma. Though significant effort has been focused on molecular pathogenesis, development, and recurrence of LSCC, little is known about its relationship with the immune-related long non-coding RNA (lncRNA) pairs.

Methods

After obtaining the transcriptome profiling data sets and the corresponding clinical characteristics of LSCC patients and normal samples from The Cancer Genome Atlas (TCGA) database, a series of bioinformatic analysis was conducted to select the differently expressed immune-related lncRNAs and build a signature of immune-related lncRNA pairs. Then, the effectiveness of the signature was validated.

Results

A total of 111 LSCC patients and 12 normal samples’ transcriptome profiling data sets were retrieved from TCGA. 301 differently expressed immune-related lncRNAs were identified and 35,225 lncRNA pairs were established. After univariate Cox analysis, LASSO regression and multivariate Cox analysis, 7 lncRNA pairs were eventually selected to construct a signature. The riskscore was computed using the following formula: Riskscore = 0.95 × (AL133330.1|AC132872.3) + (-1.23)  × (LINC01094|LINC02154) + 0.65 × (LINC02575|AC122685.1) + (-1.15)  × (MIR9-3HG|LINC01748) + 1.45 × (AC092687.3|SNHG12) + (-0.87)  × (AC090204.1|AL158166.1) + 0.64 × (LINC01063|Z82243.1). Patients were classified into the high-risk group (> 1.366) and the low-risk group (< 1.366) according to the cutoff value (1.366), which is based on the 5-year riskscore ROC curve. The survival analysis showed that the low-risk group had a better prognosis (P < 0.001). The riskscore was better than other clinical characteristics in prognostic prediction and the area under the curves (AUCs) for the 1-, 3-, and 5-year survivals were 0.796, 0.946, and 0.895, respectively. Combining age, gender, grade, stage, and riskscore, a nomograph was developed to predict survival probability in LSCC patients. Then, the riskscore was confirmed to be related with the content of tumor-infiltration immune cells and the model could serve as a potential predictor for chemosensitivity.

Conclusion

We successfully established a more stable signature of 7 immune-related lncRNA pairs, which has demonstrated a better prognostic ability for LSCC patients and may assist clinicians to precisely prescribe chemo drugs.

Peer Review reports

Introduction

Laryngeal squamous cell carcinoma (LSCC) is the most common squamous cell carcinoma of the head and neck and the second most prevalent malignant tumor of the respiratory tract [1]. Important risk factors for LSCC include smoking, drinking, low body mass index, laryngopharyngeal reflux, etc. [2]. LSCC treatment mainly integrates surgery, radiation, and chemotherapy [3]. However, the survival rates for most patients with advanced LSCC have not improved significantly in the last 30 years [4]. The larynx plays a central role in maintaining physiologic processes, including respiration, speech [5]. Therefore, understanding the molecular mechanisms underlying the occurrence and development of LSCC is essential for its diagnosis, treatment, and prognosis.

Long non-coding RNAs (lncRNAs) are a class of RNA molecules that are longer than 200 nucleotides in length and rarely encode for proteins [6, 7]. Recently, the accumulation of lncRNAs has been implicated in physiological and pathological processes, especially in cancers [8]. Specifically, lncRNA SATB2-AS1was reported to inhibit colorectal cancer (CRC) cell metastasis and regulate the immune response of CRC by cis-activating SATB2 [6]. LncRNAs have also been shown to promote malignant biological behaviors [3] and chemotherapy resistance in LSCC [1, 9]. Growing evidence suggests that lncRNAs may be potential biomarkers and therapeutic targets in various cancer types [10]. However, to the best of our knowledge, no prognostic signature of immune-related lncRNA pairs has been studied in LSCC.

Herein, we constructed a signature of immune-related lncRNA pairs and validated its effectiveness in LSCC.

Materials and methods

Data acquisition and preprocessing

Transcriptome profiling data sets and the corresponding clinical characteristics of LSCC patients and normal samples were obtained from The Cancer Genome Atlas (TCGA) database (https://portal.gdc.cancer.gov).

Differential expression analysis and construction of lncRNA pairs

Immune-related genes were screened according to the immune-related gene list obtained from the IMMPORT data source (https://www.immport.org/resources). Immune-related lncRNAs were selected by co-expression analysis based on the expression levels of lncRNAs and immune-related genes in the downloaded samples (|correlation coefficient|> 0.4, P < 0.001). Differently expressed immune-related lncRNAs were identified based on the previously selected expression level of lncRNAs between tumor and non-tumor specimens (|logFC|> 1.5, FDR < 0.05). The identified differently expressed immune-related lncRNAs were used to develop differential immune-related lncRNA pairs. The differential immune-related lncRNAs were cyclically singly paired. In the sample, when the level of expression of the first lncRNA was higher compared to the second, it was marked as 1, otherwise, it was marked as 0. Matches were considered effective when the expression level of 0 or 1 lncRNA pairs was between 20 and 80% of the total pairs.

Cell cultures and quantitative real-time PCR

Quantitative real-time PCR was carried out to verify the expression of these genes. The human HEK293T and LSCC cell line TU686 were purchased from BeNa Culture Collection (BNCC, Beijing, China). TU686 cells were cultured in RPMI 1640 medium (Procell Life Science &Technology, Wuhan, China) supplemented with 10% fetal bovine serum (FBS) (Procell) . HEK293T cells were cultured in DMEM (Procell) containing 10% fetal bovine serum, 100 U/mL penicillin, and 100 μg/mL streptomycin. All cells were cultured at 37˚C under 5% CO2. Total RNA was isolated from cell lines using TRIzol reagent (Invitrogen, Shanghai, China). The cDNA was synthesized using 2×NovoScript Plus 1 st Strand cDNA Synthesis SuperMix (Novoprotein, Suzhou,China). The 20 μL reaction system contained 10 μL SYBR (Vazyme Biotech, Nanjing, China), 0.5μL forward primer, 0.5 μL reverse primer, 1.0 μL cDNA, and 8.0 μL DEPC water. The amplification program was as follows: initial denaturation step at 95◦C for 2 mins, followed by 45 cycles at 95◦C for 10 s, and 60◦C for 45 s. The primers were designed and synthesized by(Sangon Biotech , Shanghai, China). The expression of each gene was normalized by expression of GAPDH. We selected 2 lncRNA pairs from 7 lncRNA pairs.The primer sequences used in this study were as follows: GAPDH forward: 5’ -ATGGCCTTCCGTGTCCC-3’; GAPDH reverse: 5’‐GAGGAGTGGGTGTCGCTGT‐3’. MIR9-3HG forward: 5’-GCTGGGCTTCTCCTGCAATTC-3’; MIR9-3HG reverse: 5’-AACCCTGTAAACCCCTCCTCTCAC-3’. LINC01748 forward: 5’-GCACTCCACAGAAGGAAGGTTCAG-3’; LINC01748 reverse: 5’-GGAAGTGATTGCCAGGAACAGGAC-3’; AC092687.3 forward: 5’-TGGAGTCCGTGCTCAGGTATGG-3’, ACO92687.3 reverse: 5’- ACGCTTGGTGGATGACTTCTTATGC-3’; SNHG12 forward: 5’-CTGAGGAGGTGAGCTTGTTTCTGG-3’, SNHG12 reverse: 5’-GTCCTTGCCTTCTGCTTCCCATAG-3’.

Construction of the prognostic signature and calculation of the risk score

Based on clinical data of LSCC patients in the TCGA, univariate Cox analysis was used to screen lncRNA pairs related to prognosis. LASSO regression was then performed to further select valuable lncRNA pairs. Then, multivariate Cox regression was performed to construct a prognostic signature of lncRNA pairs and calculate the coefficients. The riskscore of this signature was computed using the following formula: Riskscore = ∑ coef (pairi) × value (pairi) (i = 1-n), where coef indicates the contribution to the risk score and value represented the lncRNA pair results,0 or 1, in the sample. The 1-, 3-, and 5-year receiver operating characteristic (ROC) curves of the model were plotted.

Validation of the established risk model

The cutoff value was calculated based on the 5-year survival ROC curve and used to stratify patients into high- and low-risk groups. The survival state was compared between the two groups. Univariate and multivariate Cox regression analyses were performed to confirm whether the riskscore could independently predict the prognosis of LSCC patients. The prognostic relationship between the riskscore and other clinical parameters was explored. A nomograph was developed for predicting survival probability in LSCC patients.

Correlation analysis of tumor-infiltrating immune cells and riskscore values

The infiltration estimation results for all TCGA tumors were downloaded from the TIMER2.0 database (http://timer.cistrome.org/). These results were analyzed with TIMER [11], CIBERSORT, CIBERSORT-ABS [12], QUANTISEQ [13], xCell [14], MCP-counter [15], and EPIC [16] methods. Spearman correlation analysis was employed to analyze the relationship between tumor-infiltrating immune cells and riskscore values. Meanwhile, the differences in the expression level of immune-related genes were compared between the high- and low-risk groups.

Analysis of the significance of the model in the clinical treatment

To explore the clinical value of the model, the half-maximal inhibitory concentration (IC50) of common chemotherapeutic drugs, such as cisplatin, bleomycin, and paclitaxel, in the TCGA LSCC dataset was calculated. The Wilcoxon signed-rank test was used to analyze the difference in the IC50 between the high- and low-risk groups, using the pRRophetic package in R.

Statistical analysis

LncRNA expression in the tumor and normal tissues was compared by the Wilcoxon test. Univariate and multivariate Cox proportional hazard regression analyses were utilized to explore the correlation between riskscore values and overall survival (OS). The Kaplan–Meier method was used to evaluate the OS. A log-rank test was used to compare the differences between the two groups. Results from at least three independent experiments were presented as mean ± standard deviation (SD). Student’s t test was used to determine the significance of differences in groups.The statistical analysis software included R 4.0.3 and Graphpad Prism 6.0. A P value < 0.05 was considered significant.

Results

Data acquisition and preprocessing

LSCC transcriptome profiling data sets comprising of 111 patients and 12 normal samples were retrieved from TCGA(Table1) . Patients’ age ranged from 38-83 years. Ninety patients were males and 21 were females. Fifty patients were deceased. The flowchart of the study is shown in Figure 1.

Table 1 Clinicopathological features of 111 LSCC patients
Fig. 1
figure 1

Flow chart of the study

Identification of the immune-related lncRNAs and construction of lncRNA pairs

A total of 1,710 immune-related mRNAs (Supplementary File 1) were selected according to the gene list downloaded from the IMMPORT database, followed by co-expression analysis of known immune-related genes and lncRNAs. Overall, 1105 immune-related lncRNAs (Supplementary File 2) and 301differently expressed immune-related lncRNAs (Supplementary File 3) were identified (Figure 2A, B). Using iterative cycling and 0-or-1 matrix screening, 35,225 lncRNA pairs were established (Supplementary File 4).

Fig. 2
figure 2

Identification of differently expressed immune-related lncRNAs and construction of a lncRNA pairs signature. A The heatmap and (B) volcano plot of the differently expressed immune-related lncRNAs. C The generalized cross-validation curve of paired likelihood deviance. D The generalized cross-validation curve of coefficients. E A forest map showed 7 differently expressed immune-related lncRNA pairs identified by Cox proportional hazard regression in the stepwise method

Construction of the prognostic signature of immune-related lncRNA pairs

After univariate Cox analysis and 711 prognosis-related lncRNA pairs were identified. Further screening by modified LASSO regression retrieved 13 lncRNA pairs (Fig. 2C, D). Multivariate Cox analysis was then performed and 7 lncRNA pairs were eventually selected to construct a signature of immune-related lncRNA pairs (Table 2). The riskscore was computed using the following formula: Riskscore = 0.95 × (AL133330.1|AC132872.3) + (-1.23) × (LINC01094|LINC02154) + 0.65 × (LINC02575|AC122685.1) + (-1.15) × (MIR9-3HG|LINC01748) + 1.45 × (AC092687.3|SNHG12) + (-0.87) × (AC090204.1|AL158166.1) + 0.64 × (LINC01063|Z82243.1). A forest map of the 7 lncRNA pairs is shown in Fig. 2E. Each sample was scored using this formula. The ROC curves for 1-, 3-, and 5- years were drawn using the riskscore values (Fig. 3A). To verify the expression of these lncRNAs in LSCC, we selected two pairs of more common lncRNA pairs for qRT-PCR, and the results were shown in Supplementary Fig. 1A,B, which showed that the expression of these genes was elevated in LSCC and in accordance with our matching rules.

Table 2 Seven prognostic immune-related lncRNA pairs identified from univariate Cox regression analysis, LASSO regression, and multivariate Cox regression analysis
Fig. 3
figure 3

The cut-off value and survival curve between high- and low-risk groups. A The 1-, 3-, and 5-year ROC of the optimal model. B The maximum inflection point is the cut-off point obtained by the AIC in 5-year ROC curve. C The survival curve between high- and low-risk groups. D Riskscores and survival outcome (E) of each case

Validation of the established risk model

Based on the 5-year riskscore ROC curve, the cutoff value was set as 1.366 (Figure 3B). Patients were then classified into the high-risk group (>1.366) and the low-risk group (<1.366). The survival rate was compared between the two groups and the result showed that the low-risk group had a better prognosis (P<0.001) (Figure 3C, D, E). Univariate and multivariate Cox analyses confirmed that the riskscore was an independent prognostic factor of LSCC (Figure 4A, B). Moreover, the results from the ROC curves revealed that compared with other clinical characteristics, the riskscore was better in prognostic prediction in the 1, 3, and 5 years and the area under the curves (AUCs) were 0.796, 0.946, and 0.895, respectively (Figure 4C, D, E). We also compared the differences in the riskscore between other clinical characteristics (Figure 5A-H). Our results demonstrated that N3 patients had a higher riskscore than N0 patients. To accurately predict survival probability in LSCC patients, we developed a nomograph combining age, gender, grade, stage, and riskscore (Figure 6). In the nomograph, a higher point represented a poor prognosis.

Fig. 4
figure 4

Riskscore was an independent factor in predicting prognosis. (A)Univariate and (B) multivariate Cox regression analysis results were showed by forest maps. (C-E) Riskscore was better than other clinical characteristics in prognostic prediction in 1,3,5- year.

Fig. 5
figure 5

Correlations between the riskscore and other clinical characteristics. A strip chart (A) and the scatter diagrams showed the correlations between riskscore and age (B), gender (C), grade (D), stage (E), T(F), N(G), M(H)

Fig. 6
figure 6

A nomogram was developed to predict survival probability. A A nomogram was developed in which age, gender, grade, stage and riskscore were integrated. B The calibration curve of 3-year survival. C The calibration curve of 5-year survival

Correlation analysis of tumor-infiltrating immune cells and riskscore values

Combining with the downloaded results from the TIMER2.0 database, we studied the correlation between the riskscore and the content of tumor-infiltration immune cells (Fig. 7A-I). We found that B cell plasma, B cell, cytotoxicity score, T cell CD4+ effector memory, T cell CD8+, and T cell follicular helper were higher in the low-risk group. We also assessed the expression of immune-related genes in different groups. TNFRSF4 and TNFRSF18 were highly expressed in the low-risk group while versican (VCAN) was highly expressed in the high-risk group (Fig. 8A-C).

Fig. 7
figure 7

Exploration the relationship between tumor-infiltration immune cells and riskscore. (A) tumor-infiltration cells, such as B cell plasma, B cell, cytotoxicity score, T cell CD4 + effector memory, T cell CD8 + , T cell follicular helper, were negatively associated with the riskscore. The scatter diagrams showed that B cell plasma (B), B cell (C), cytotoxicity score (D), T cell CD4 + effector memory (E), T cell CD8 + (F, G), T cell follicular helper (H, I) were higher in low-risk group

Fig. 8
figure 8

Immune related genes were differently expressed between high-and low risk groups and riskscore was a potential predictor for chemotherapeutics.VCAN (A) was highly expressed in high-risk group while TNFRSF4 (B) and TNFRSF18 (C) were highly expressed in low-risk group. Bexarotene (D) and bicalutamide (E) had lower IC50 value in high-risk group

Association between the risk model and chemotherapeutics

We also investigated the correlation between risk and efficacy of common chemotherapeutics in the TCGA LSCC dataset. The results indicated that a high riskscore was associated with a lower IC50 of chemotherapeutics such as bexarotene and bicalutamide, suggesting that the model served as a potential predictor for chemosensitivity (Figure 8D, E).

Discussion

Based on LSCC datasets from TCGA, we selected differently expressed immune-related lncRNAs, matched lncRNA pairs, and then successfully established a signature of immune-related lncRNA pairs. Using this signature, we scored each sample and drew 1-, 3-, and 5- year riskscore ROC curves, respectively. Subsequently, we divided patients into high- and low-risk groups according to the cutoff value of 5- year riskscore ROC curve. Riskscore was confirmed to be an independent prognostic factor of LSCC. Moreover, our result provided an accurate nomograph for predicting survival probability in LSCC patients.

LSCC is one of the most common subtypes of head and neck squamous cell carcinoma (HNSCC). Over the past decades, the mortality of LSCC has significantly increased worldwide; therefore, early detection and precise treatment are crucial [17]. Significant effort has been focused on molecular pathogenesis, development, and recurrence of LSCC. Numerous studies have proven that lncRNAs could serve as potential oncogenes or tumor suppressors in tumorigenesis and tumor progression of several types of cancers [18]. For example, the expression of lncRNA retinoblastoma-associated transcript-1 (RBAT1) was significantly higher in both retinoblastoma (Rb) and bladder cancer (BCa) clinical tissues compared with normal tissues and RBAT promoted tumorigenesis by activating E2F3 transcription [19]. In lung cancer, lncRNA JPX was reported to upregulate Twist1 by competitively sponging miR-33a-5p, subsequently inducing EMT and lung cancer cell invasion [7]. Specifically, overexpression of lncRNA MIR31HG was correlated with aggressive clinicopathological features and could serve as a poor prognostic biomarker in LSCC [20]. While LINC‐PINT expression was downregulated in laryngeal carcinoma tissues, it acted as a suppressor through the LINC-PINT/miR-425-5p/PTCH1 axis [9]. Because tumor heterogeneity limits therapeutic efficacy, using a single prognostic biomarker in LSCC patients could be unreliable [21]. Therefore, the use of multiple biomarkers has been suggested in the predictive prognosis of laryngeal cancer and the development of new therapeutic strategies. Gong et al. successfully established a 6-lncRNA model, which had good performance for predicting survival in LSCC patients to aid in risk stratification and provide precise therapeutic recommendations [22]. Zhang et al. identified a 4-lncRNA signature predicting the prognosis of patients with laryngeal cancer [23]. They reported that the signature could influence the prognosis of laryngeal cancer by regulating tumor apoptosis, metastasis, invasion, immunity, and other characteristics through the Notch and Wnt signaling pathway and voltage-gated calcium channels. In the present study, we successfully constructed a signature of 7-lncRNA pairs. To reduce the differences in the expression levels of the same lncRNAs from different samples, we adopted an lncRNA-pairing method where only the expression levels of two lncRNAs in the same sample are compared instead of comparing these lncRNAs in different samples. In addition, two pairs of lncRNAs were selected for qRT-PCR verification of their expression in LSCC, and the results were consistent with our matching principle. Practically, our signature had better prognostic prediction performance compared with other clinical characteristics such as stage, age, gender, and grade.

Malignant tumors are mainly composed of cancerous cells and are influenced by a complex microenvironment. The contexture of immune cells in the tumor microenvironment (TME) profoundly influences tumor progression and the effectiveness of anti-cancer therapies [24]. Although HNSCC modalities such as surgery, radiotherapy, and chemotherapy have significantly improved over the past decades, recurrence rates are still high [25]. Given that HNSCCs are tumors with a rich immune infiltrate, HNSCC treatment is geared towards immunotherapy [25]. Association between TME and immune infiltration in LSCC has been previously investigated [26,27,28,29,30]. To explore the association between riskscores and tumor-infiltrating immune cells, we used common acceptable methods, including TIMER, CIBERSORT, QUANTISEQ, xCELL, MCP-counter, and EPIC, to estimate immune-infiltrating cells. After analysis, our results demonstrated that B cell plasma, B cell, cytotoxicity score, T cell CD4+ effector memory, T cell CD8+, and T cell follicular helper were higher in the low-risk group. This suggested that more immune-infiltrating cells were associated with a better prognosis, which was in accord with previous results. Zhou et al. [28] reported higher densities of CD8+ cell infiltration in early tumor stages than in late tumor stages in LSCC. They also found that a high density of CD8+ immune cells in both peritumoral and intratumoral regions was significantly associated with a favorable OS. In our study, we assessed the expression of immune-related genes in the high- and low-risk groups, and our results indicated that TNFRSF4 and TNFRSF18 were highly expressed in the low-risk group while VCAN was highly expressed in the high-risk group. In line with our result, it was reported that gastric cancer patients with high VCAN expression had worse prognoses than those with low VCAN expression and VCAN was an independent risk factor for OS [31]. Contrary to our result, Gu et al. [32] demonstrated that a higher TNFRSF4 expression level contributed to the poor prognosis of patients with non-M3 acute myeloid leukemia (AML). This discrepancy in results could be due to different disease types or sample sizes. Since chemotherapy is one of the important treatments in LSCC, we compared the IC50 of different chemotherapeutic drugs in high- and low-risk groups. We found that a high riskscore was associated with a lower IC50 of bexarotene and bicalutamide, suggesting that our model served as a potential predictor for chemosensitivity.

With the advancement of gene sequencing techniques coupled with traditional single biomarker analyses, signatures composed of multiple genes, miRNAs, or lncRNAs have attracted more attention and may offer practical solutions by providing treatment suggestions and predicting prognosis. Zhang et al. [33] established a five-gene signature that could effectively predict the prognosis of laryngeal cancer. Cui et al. [34] explored a genome-wide integrated analysis of methylation and the transcriptome to construct a methylation-driven gene-based prognostic model to precisely predict recurrence probability and optimize therapeutic strategies for LSCC.

Meanwhile, a novel 4-lncRNA signature was identified in patients with laryngeal cancer [35]. In our study, we combined the immune-related mRNA and lncRNA expression level to successfully construct a signature of immune-related lncRNA pairs and explored its clinical significance. To the best of our knowledge, this is the first report to investigate a signature of immune-related lncRNA pairs in LSCC. We selected 7 lncRNA pairs: AL133330.1|AC132872.3, LINC01094|LINC02154, LINC02575|AC122685.1, MIR9-3HG|LINC01748, AC092687.3|SNHG12, AC090204.1|AL158166.1, and LINC01063|Z82243.1. Among these lncRNAs, some have been reported to play essential roles in different types of tumors. For example, LINC01094 was found to upregulate SLC2A3 by targeting microRNA-184, subsequently promoting the development of clear cell renal cell carcinoma [36]. SNHG12 is a potential biomarker in various cancer types, including gastric cancer, non-small cell lung cancer, triple-negative breast cancer, hepatocellular carcinoma, etc. Therefore, targeting SNHG12 could lead to advances in the diagnoses, prognosis, and/or treatment of these cancers [37].

However, this study has some shortcomings. Firstly, our results were based on the TCGA LSCC data set, which was relatively insufficient. Secondly, although we confirmed our signature in the downloaded data, some lncRNAs in the signature have not been published before, limiting validation by other databases. Therefore, our signature should be verified with large sample sizes and basic experiments.

Conclusion

In conclusion, we successfully established a signature of immune-related lncRNA pairs that is more stable by neglecting the differences of samples. Our signature has demonstrated a better prognostic ability for LSCC patients and may assist clinicians to precisely prescribe chemo drugs. Certainly, there are also several limitations in our results. Our current datas are all downloaded from public databases, and we have not carried out comprehensive experimental verification.

Availability of data and materials

The datasets of this study are available on request to the corresponding author.

Abbreviations

LSCC:

Laryngeal Squamous Cell Carcinoma

lncRNA:

Long non-coding RNA

TCGA:

The Cancer Genome Atlas

FDR:

False Discovery Rate

OS:

Overall Survival

ROC:

Receiver Operating Characteristic

AUC:

Area Under Curve

HNSCC:

Head and Neck Squamous Cell Carcinoma

CRC:

Colorectal cancer

TIMER:

Tumor Immune Estimation Resource

VCAN:

Versican

TNFRSF:

Tumor Necrosis Factor Receptor Super Family

Rb:

Retinoblastoma

BCa:

Bladder cancer

E2F3:

E2F Transcription Factor 3

EMT:

Epithelial-Mesenchymal Transition

MIR31HG:

MIR31 Host Gene

PTCH1:

Patched 1

TME:

Tumor micro-environment

AML:

Acute Myeloid Leukemia

SLC2A3:

Solute Carrier Family 2 Member 3

SNHG:

Small Nucleolar RNA Host Gene

References

  1. Li R, Chen S, Zhan J, Li X, Liu W, Sheng X, Lu Z, Zhong R, Chen L, Luo X, et al. Long noncoding RNA FOXD2-AS1 enhances chemotherapeutic resistance of laryngeal squamous cell carcinoma via STAT3 activation. Cell Death Dis. 2020;11(1):41.

    Article  Google Scholar 

  2. Wei Q, Yu D, Liu M, Wang M, Zhao M, Liu M, Jia W, Ma H, Fang J, Xu W, et al. Genome-wide association study identifies three susceptibility loci for laryngeal squamous cell carcinoma in the Chinese population. Nat Genet. 2014;46(10):1110–4.

    CAS  Article  Google Scholar 

  3. Li X, Xu F, Meng Q, Gong N, Teng Z, Xu R, Zhao M, Xia M. Long noncoding RNA DLEU2 predicts a poor prognosis and enhances malignant properties in laryngeal squamous cell carcinoma through the miR-30c-5p/PIK3CD/Akt axis. Cell Death Dis. 2020;11(6):472.

    Article  Google Scholar 

  4. Gao W, Zhang C, Li W, Li H, Sang J, Zhao Q, Bo Y, Luo H, Zheng X, Lu Y, et al. Promoter Methylation-Regulated miR-145-5p Inhibits Laryngeal Squamous Cell Carcinoma Progression by Targeting FSCN1. Mol Ther. 2019;27(2):365–79.

    CAS  Article  Google Scholar 

  5. Zhang L, Wu Y, Zheng B, Su L, Chen Y, Ma S, Hu Q, Zou X, Yao L, Yang Y, et al. Rapid histology of laryngeal squamous cell carcinoma with deep-learning based stimulated Raman scattering microscopy. Theranostics. 2019;9(9):2541–54.

    CAS  Article  Google Scholar 

  6. Xu M, Xu X, Pan B, Chen X, Lin K, Zeng K, Liu X, Xu T, Sun L, Qin J, et al. LncRNA SATB2-AS1 inhibits tumor metastasis and affects the tumor immune cell microenvironment in colorectal cancer by regulating SATB2. Mol Cancer. 2019;18(1):135.

    Article  Google Scholar 

  7. Pan J, Fang S, Tian H, Zhou C, Zhao X, Tian H, He J, Shen W, Meng X, Jin X, et al. lncRNA JPX/miR-33a-5p/Twist1 axis regulates tumorigenesis and metastasis of lung cancer by activating Wnt/beta-catenin signaling. Mol Cancer. 2020;19(1):9.

    CAS  Article  Google Scholar 

  8. Liang Y, Song X, Li Y, Chen B, Zhao W, Wang L, Zhang H, Liu Y, Han D, Zhang N, et al. LncRNA BCRT1 promotes breast cancer progression by targeting miR-1303/PTBP3 axis. Mol Cancer. 2020;19(1):85.

    CAS  Article  Google Scholar 

  9. Yuan Z, Xiu C, Liu D, Zhou G, Yang H, Pei R, Ding C, Cui X, Sun J, Song K. Long noncoding RNA LINC-PINT regulates laryngeal carcinoma cell stemness and chemoresistance through miR-425-5p/PTCH1/SHH axis. J Cell Physiol. 2019;234(12):23111–22.

    CAS  Article  Google Scholar 

  10. Cao R, Yuan L, Ma B, Wang G, Tian Y. Immune-related long non-coding RNA signature identified prognosis and immunotherapeutic efficiency in bladder cancer (BLCA). Cancer Cell Int. 2020;20:276.

    CAS  Article  Google Scholar 

  11. 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.

    Article  Google Scholar 

  12. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, Hoang CD, Diehn M, Alizadeh AA. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.

    CAS  Article  Google Scholar 

  13. Finotello F, Mayer C, Plattner C, Laschober G, Rieder D, Hackl H, Krogsdam A, Loncova Z, Posch W, Wilflingseder D, et al. Molecular and pharmacological modulators of the tumor immune contexture revealed by deconvolution of RNA-seq data. Genome Med. 2019;11(1):34.

    Article  Google Scholar 

  14. Aran D, Hu Z, Butte AJ. xCell: digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. 2017;18(1):220.

    Article  Google Scholar 

  15. Becht E, Giraldo NA, Lacroix L, Buttard B, Elarouci N, Petitprez F, Selves J, Laurent-Puig P, Sautes-Fridman C, Fridman WH, et al. Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression. Genome Biol. 2016;17(1):218.

    Article  Google Scholar 

  16. Racle J, de Jonge K, Baumgaertner P, Speiser DE, Gfeller D. Simultaneous enumeration of cancer and immune cell types from bulk tumor gene expression data. Elife. 2017;6:e26476.

    Article  Google Scholar 

  17. Meng W, Cui W, Zhao L, Chi W, Cao H, Wang B. Aberrant methylation and downregulation of ZNF667-AS1 and ZNF667 promote the malignant progression of laryngeal squamous cell carcinoma. J Biomed Sci. 2019;26(1):13.

    Article  Google Scholar 

  18. He A, He S, Li X, Zhou L. ZFAS1: A novel vital oncogenic lncRNA in multiple human cancers. Cell Prolif. 2019;52(1):e12513.

    Article  Google Scholar 

  19. He X, Chai P, Li F, Zhang L, Zhou C, Yuan X, Li Y, Yang J, Luo Y, Ge S, et al. A novel LncRNA transcript, RBAT1, accelerates tumorigenesis through interacting with HNRNPL and cis-activating E2F3. Mol Cancer. 2020;19(1):115.

    CAS  Article  Google Scholar 

  20. Wang R, Ma Z, Feng L, Yang Y, Tan C, Shi Q, Lian M, He S, Ma H, Fang J. LncRNA MIR31HG targets HIF1A and P21 to facilitate head and neck cancer cell proliferation and tumorigenesis by promoting cell-cycle progression. Mol Cancer. 2018;17(1):162.

    CAS  Article  Google Scholar 

  21. Song L, Zhang S, Yu S, Ma F, Wang B, Zhang C, Sun J, Mao X, Wei L. Cellular heterogeneity landscape in laryngeal squamous cell carcinoma. Int J Cancer. 2020;147(10):2879–90.

    CAS  Article  Google Scholar 

  22. Gong S, Xu M, Zhang Y, Shan Y, Zhang H. The Prognostic Signature and Potential Target Genes of Six Long Non-coding RNA in Laryngeal Squamous Cell Carcinoma. Front Genet. 2020;11:413.

    CAS  Article  Google Scholar 

  23. Zhang G, Fan E, Zhong Q, Feng G, Shuai Y, Wu M, Chen Q, Gou X. Identification and potential mechanisms of a 4-lncRNA signature that predicts prognosis in patients with laryngeal cancer. Hum Genomics. 2019;13(1):36.

    CAS  Article  Google Scholar 

  24. Sturm G, Finotello F, Petitprez F, Zhang JD, Baumbach J, Fridman WH, List M, Aneichyk T. Comprehensive evaluation of transcriptome-based cell-type quantification methods for immuno-oncology. Bioinformatics. 2019;35(14):i436–45.

    CAS  Article  Google Scholar 

  25. Horton JD, Knochelmann HM, Day TA, Paulos CM, Neskey DM. Immune Evasion by Head and Neck Cancer: Foundations for Combination Therapy. Trends Cancer. 2019;5(4):208–32.

    CAS  Article  Google Scholar 

  26. Alessandrini L, Franz L, Ottaviano G, Ghi MG, Lanza C, Blandamura S, Marioni G. Prognostic role of programmed death ligand 1 (PD-L1) and the immune microenvironment in laryngeal carcinoma. Oral Oncol. 2020;108:104836.

    CAS  Article  Google Scholar 

  27. Zeng H, Huang Y, Chen L, Li H, Ma X. Exploration and validation of the effects of robust co-expressed immune-related genes on immune infiltration patterns and prognosis in laryngeal cancer. Int Immunopharmacol. 2020;85:106622.

    CAS  Article  Google Scholar 

  28. Zhou L, Li Y, Gao W, Huangfu H, Wen S, Zhang C, Zhao Q, Dong Z, Qu C, Li G, et al. Assessment of tumor-associated immune cells in laryngeal squamous cell carcinoma. J Cancer Res Clin Oncol. 2019;145(7):1761–72.

    CAS  Article  Google Scholar 

  29. Zhang D, Zhou J, Tang D, Zhou L, Chou L, Chou KY, Tao L, Lu LM. Neutrophil infiltration mediated by CXCL5 accumulation in the laryngeal squamous cell carcinoma microenvironment: A mechanism by which tumour cells escape immune surveillance. Clin Immunol. 2017;175:34–40.

    CAS  Article  Google Scholar 

  30. Karpathiou G, Casteillo F, Giroult JB, Forest F, Fournel P, Monaya A, Froudarakis M, Dumollard JM, Prades JM, Peoc’H M. Prognostic impact of immune microenvironment in laryngeal and pharyngeal squamous cell carcinoma: Immune cell subtypes, immuno-suppressive pathways and clinicopathologic characteristics. Oncotarget. 2017;8(12):19310–22.

  31. Li W, Han F, Fu M, Wang Z. High expression of VCAN is an independent predictor of poor prognosis in gastric cancer. J Int Med Res. 2020;48(1):1219690823.

    Google Scholar 

  32. Gu S, Zi J, Han Q, Song C, Ge Z. Elevated TNFRSF4 gene expression is a predictor of poor prognosis in non-M3 acute myeloid leukemia. Cancer Cell Int. 2020;20:146.

    CAS  Article  Google Scholar 

  33. Zhang G, Fan E, Yue G, Zhong Q, Shuai Y, Wu M, Feng G, Chen Q, Gou X. Five genes as a novel signature for predicting the prognosis of patients with laryngeal cancer. J CELL BIOCHEM. 2019;1–10. https://doi.org/10.1002/jcb.29535.

  34. Cui J, Wang L, Zhong W, Chen Z, Chen J, Yang H, Liu G. Identification and validation of methylation-driven genes prognostic signature for recurrence of laryngeal squamous cell carcinoma by integrated bioinformatics analysis. Cancer Cell Int. 2020;20:472.

    CAS  Article  Google Scholar 

  35. Zhang G, Fan E, Zhong Q, Feng G, Shuai Y, Wu M, Chen Q, Gou X. Identification and potential mechanisms of a 4-lncRNA signature that predicts prognosis in patients with laryngeal cancer. Hum Genomics. 2019;13(1):36.

    CAS  Article  Google Scholar 

  36. Xu H, Wang X, Wu J, Ji H, Chen Z, Guo H, Hou J. Long Non-coding RNA LINC01094 Promotes the Development of Clear Cell Renal Cell Carcinoma by Upregulating SLC2A3 via MicroRNA-184. Front Genet. 2020;11:562967.

    CAS  Article  Google Scholar 

  37. Tamang S, Acharya V, Roy D, Sharma R, Aryaa A, Sharma U, Khandelwal A, Prakash H, Vasquez KM, Jain A. SNHG12: An LncRNA as a Potential Therapeutic Target and Biomarker for Human Cancer. Front Oncol. 2019;9:901.

    Article  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This study was supported by Jiangsu Post-doctoral Foundation Research Project, China (No.2020Z136).

Author information

Authors and Affiliations

Authors

Contributions

HSN conceived this work and LSQ wrote the paper. TTN collected and preprocessed the data from TCGA. HS performed the analysis and prepared the figures and tables. HSN and TTN helped to interpret the results. LSQ revised the manuscript. BF was responsible for the revision and experimental work of the paper, and checked and corrected the paper again. HSN supervised the entire study. All authors have read and approved the manuscript. LSQ, TTN and BF contributed equally to this study.

Corresponding author

Correspondence to Haosheng Ni.

Ethics declarations

Ethics approval and consent to participate

Not applicable. All data in this study are publicly available and no permission was required to perform this study.

The study protocol was approved by the Clinical Research Ethics Committee of the Affiliated Hospital of Nantong University.

Consent for publication

Not Applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

(TIF 55 kb)

Additional file 2.

 A total of 1710 immune-related mRNAs were identified.

Additional file 3.

 A total of 1105 immune-related lncRNAs were identified by co-expression analysis.

Additional file 4.

 A total of 301 differently expressed immune-related lncRNAs were identified.

Additional file 5.

 A total of 35225 lncRNA pairs were constructed.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visithttp://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Qian, L., Ni, T., Fei, B. et al. An immune-related lncRNA pairs signature to identify the prognosis and predict the immune landscape of laryngeal squamous cell carcinoma. BMC Cancer 22, 545 (2022). https://doi.org/10.1186/s12885-022-09524-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-022-09524-1

Keywords

  • Laryngeal squamous cell carcinoma (LSCC)
  • The Cancer Genome Atlas (TCGA)
  • Long non-coding RNA (lncRNA)
  • Tumor-infiltrating immune cell
  • Prognostic signature