A Radioresponse-Related lncRNA Biomarker Signature for Risk Classification and Prognosis Prediction in Non-Small-Cell Lung Cancer

Song, Jiahang; Zhang, Shuming; Sun, Yuanyuan; Gu, Junjie; Ye, Ziqi; Sun, Xinchen; Tang, Qiyun

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

Journal of Oncology

On this page

Abstract Introduction Methods Results Discussion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Special Issue

Bioinformatics in Cancer and Immune Microenvironment

View this Special Issue

Research Article | Open Access

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

A Radioresponse-Related lncRNA Biomarker Signature for Risk Classification and Prognosis Prediction in Non-Small-Cell Lung Cancer

Jiahang Song,^1,2Shuming Zhang,¹Yuanyuan Sun,²Junjie Gu,²Ziqi Ye,²Xinchen Sun,²and Qiyun Tang¹

Academic Editor: Jimei Wang

Received17 Aug 2021

Accepted06 Sept 2021

Published21 Sept 2021

Abstract

Purpose. Radiotherapy resistance is now recognized as the major obstacle to the effective therapeutic management of non-small-cell lung cancer (NSCLC). As a single biomarker has limited effect in stratifying NSCLC patients, this research aimed to identify long non-coding RNAs (lncRNAs) correlated with radiotherapy response to ameliorate forecast of NSCLC prognosis. Methods. In a cohort of NSCLC patients with radiotherapy history (n = 96) from TCGA, genetic data of lncRNA expression profiling were performed. To identify radioresponse-related lncRNA sets which dysregulated significantly between radiosensitive (RS) and radioresistant (RR) groups, differential expression analysis was carried out. Cox relative regression was implemented to set up a radioresponse-related risk model. Moreover, we adopted survival analysis to measure the predictive potentiality of the prognosis model. Results. Four radioresponse-related lncRNAs (CASC19, LINC01977, LINC02471, and MAGI2-AS3) were screened to create a prognostic signature. Then, we described a lncRNA signature-based regulatory network and explored the correlation of the immune microenvironment and the signature. Additionally, in vitro assays uncovered inhibition of LINC01977 weakened radioresistance of NSCLC cells. Conclusion. We provided a novel radioresponse-related lncRNAs signature with excellent clinical potency for an effective prognostic forecast of patients.

1. Introduction

Lung cancer is the most prevalent thoracic tumor with the second incidence of malignancy in the world [1]. Non-small-cell lung cancer (NSCLC) makes up nearly 85% of lung cancers, whereas 50% to 75% of patients are already diagnosed at the late stages, losing the opportunity for early surgical resection [2]. Therefore, radiotherapy (RT) has increasingly become the main clinical treatment choice. In fact, radiotherapy, as an active and effective therapeutic modality, has played a central part in the whole process management of lung cancer [3]. Nonetheless, the persistence of local control of NSCLC by radiotherapy still has the problems of inevitable recurrence and tolerance [4]. It has been reported that the incidence of local failure of radiotherapy in patients with lung cancer is 24%–40%. Even with concurrent radical radiotherapy and chemotherapy, there is a local failure rate of up to 30%–50% [5]. In conclusion, radiation resistance will still limit the long-term control of lung cancer and eventually lead to local control failure and disease progression of lung cancer [6]. Consequently, we need to further explore the molecular mechanism and develop novel biomarkers, so as to improve the response rate of radiotherapy and overcome resistance to radiotherapy in NSCLC.

lncRNAs belong to the noncoding RNA family with a length of over 200 nucleotides, which usually are not capable of encoding proteins [7]. With the gradual in-depth study of lncRNA, they have been shown to participate in a variety of tumor initiation and development [8–10]. Currently, a total of 7,942 cancer-related lncRNAs were identified to be biomarkers for specific particular tumor types by Iyer and colleagues [11]. Accordingly, better comprehension of the effects of lncRNAs in tumors can favor the development of new diagnostic markers and the exploitation of prospective therapeutic targets. Increasing evidence has explicitly pointed out that radioresistance-related lncRNAs also play crucial roles in the regulation of tumor radiation sensitivity. For instance, lncRNA KCNQ1OT1 can interact with miR-372-3p to enhance stereotactic body radioresistance by inducing autophagy in lung adenocarcinoma [12]. HOX transcript antisense RNA (HOTAIR), best known for regulating transcription, contributes to tumor radioresistance via miR-93/ATG12 [13]. However, prognostic biomarkers based on expression profiles of radioresistance-related lncRNA have not been elaborated in NSCLC.

This article aimed to set up a radioresistance-related lncRNA prognostic model with an integrated bioinformatics analysis, which predicts the prognosis of patients with NSCLC undergoing RT. First, the lncRNAs which correlated with RT response were collected from 96 TCGA-NSCLC cases with lncRNA expression profiles. Then, we built up a radioresistance-related signature consisting of four lncRNAs and verified its accuracy and reliability in different cohorts. In addition, the relationship between immunity and signature was compared by combining the ssGSEA method and immunity-associated algorithm. Finally, we confirmed that these four lncRNAs were markedly dysregulated in NSCLC cell lines by a qRT-PCR assay. Moreover, we observed that the downregulation of LINC01977 significantly blocked the clonogenic survival of NSCLC cells after irradiation.

2. Methods

2.1. Data Source and Processing

The TCGA-NSCLC lncRNA-seq expression profiles and related clinical data were collected from the TCGA project (https:/portal.gdc.cancer.gov/). A total of 96 NSCLC patients who received RT were selected for this study. Our inclusion criteria for patients were as follows: (1) histologically diagnosed with NSCLC; (2) available expression profiles; and (3) overall survival time greater than 30 days. We divided 96 cases into the radiosensitive (RS) group (n = 53) and radioresistant (RR) group (n = 43) based on the radiotherapy response. Patients who presented complete remission after RT were recognized as radiosensitive whereas those showing stable disease and progressive disease after RT were considered radioresistant. All patients’ clinical characteristics are shown in Table 1.

In order to obtain radioresponse-related lncRNAs (RRlncRNAs) in NSCLC, we analyzed the differential expression of lncRNAs between the RS group and RR group by limma in R software (|fold change (FC)| = 1.0 and ) [14].

2.2. Identification of RRlncRNAs Risk Score Model

To get an optimal RRlncRNAs model, we randomly divided 96 patients into the training set and test set at 1 : 1 ratio. In the training set, candidate prognostic lncRNAs were firstly obtained by the univariate hazard Cox method based on RRlncRNAs. Next, least absolute shrinkage and selection operator (LASSO) penalized analysis with multivariate regression method were carried to set up a novel RRlncRNAs signature through the package “glmnet.” We utilized the following formula to generate the risk score of NSCLC patients: risk score = (exp RRlncRNA 1 × coef) + (exp RRlncRNA 2 × coef) + … + (exp RRlncRNA n × coef). The exp means the expression value of each RRlncRNA and the coef is the coefficient of each RRlncRNA generated by Cox relative analysis. According to the median value of risk score, all cases with RT history were split into high- and low-risk score groups, and the differences in survival outcome were compared using K-M analysis.

2.3. Establishment and Validation of Nomogram

The data about 96 patients with complete clinical information was integrated with the risk signature, and the independent prognostic analysis was conducted by the “survival” package [15]. Furthermore, the RT response-based nomogram was established according to risk scores and other clinicopathological factors to forecast the prognosis of patients. Subsequently, the accuracy of the prognostic nomogram was verified according to the calibration curves which were generated by comparing the nomogram prediction ability and the observation for the 1-year, 3-year, and 5-year OS rates.

2.4. Gene Set Enrichment Analysis (GSEA)

We employed GSEA to distinguish the important functional phenotypes between the high-risk group and the low-risk group. After 1000 substitutions, an enriched gene set was obtained based on a nominal value <0.05. The Hallmark (v7.4) gene sets were downloaded from the Molecular Signatures Database v7.4 download page (https://www.gsea-msigdb.org/gsea/downloads.jsp). GSEA was conducted based on the downloaded gene sets by GSEA software (v4.1.0, https://www.gsea-msigdb.org/).

2.5. Immunity Analysis of the Signature

To uncover the immune activity of the RRlncRNAs signature, six immunity algorithms, including TIMER, CIBERSORT, QUANTISEQ, MCPCOUNTER, XCELL, and EPIC, were applied to evaluate immune responses between two risk subgroups by single-sample gene set enrichment analysis (ssGSEA) [16].

2.6. Competing Endogenous RNA (ceRNA) Network Construction

The DIANA prediction tools were used to detect the target miRNAs of each lncRNA and the threshold for the relationship between lncRNAs and miRNAs was set to 0.95. Then, we explored the potential target mRNAs binding to miRNAs using three miRNA databases (miRDB, TargetScan, and miRTarBase). Based on the predicted correlation of lncRNA-miRNA and miRNA-mRNA, the ceRNA network was created by Cytoscape.

2.7. Functional Enrichment Analysis of the Target mRNAs

To determine the underlying function of mRNAs in the ceRNA regulatory network, we conducted the Gene Ontology (GO) function and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis by “clusterProfiler” package. was found statistically important.

2.8. Cell Culture and Cell Transfection

Two human NSCLC cell lines (A549 and NCI-H520) and one human lung epithelial cell line (BEAS-2B) were purchase from the Chinese Academy of Sciences (Shanghai, China). Above NSCLC cells and BEAS-2B were cultured in RPMI-1640 medium (KeyGEN, Nanjing, China) with 10% fetal bovine serum (FBS, PAN-Seratech, Germany) and 1% penicillin/streptomycin with 5% CO₂ at 37°C. We purchased siRNA negative control and si-LINC01977 from Ribobio Biotechnology and transfected them into NSCLC cell lines using the Lipotransfectamine 3000 (Thermo). The target sequence of si-NC and si-LINC01977 are presented in Supplementary Table S1. After 48 h, transfected NSCLC cells were collected for the next in vitro experiments.

2.9. Quantitative Real-Time PCR (qRT-PCR)

Total RNA was extracted from cell lines via RNA-easy isolation reagent (Vazyme biotech) and was reversed into synthesizing complementary DNA (cDNA) using PrimeScript Mix reagent (Takara). The PCR reaction system was prepared to utilize SYBR Green® Premix Ex Taq™ (Vazyme biotech). Results of individual lncRNAs were normalized to the expression of GAPDH. The primer sequences for each gene are shown in Supplementary Table S2.

2.10. Clonogenic Survival Assay

NSCLC cells were seeded in 6-well plates at 200, 400, 800, 1600, and 3200 cells/well. 24 hours later, we irradiated them with a single dose of 0, 2, 4, 6, and 8 Gy, respectively. Then, cultured for 2 weeks later, 1% crystal violet was applied in staining clonogenic cells. The clonogenic survival curves were generated according to the following formula model: SF = 1 − (1 − e^−kd)ⁿ as previously described [17].

2.11. Immunofluorescence Assay

The immunofluorescence of γH2AX was detected by the anti-γH2AX primary antibody (1 : 500, Abcam, UK) and the detailed process was described in the previous protocol [17]. Cells were snapped under a confocal fluorescence microscope (Leica).

2.12. Statistical Analysis

R software (4.0.1) was utilized for all statistical analyses. Differences in OS of both risk groups were assessed through the Kaplan–Meier analysis. By the combination of univariate and multivariate Cox regression methods, the independence of our risk model was detected. In addition, receiver operating characteristic (ROC) analysis was utilized to examine the reliability of our risk model. was considered statistically significant for each analysis.

3. Results

3.1. Determination of Radioresponse-Related lncRNAs

The infographic flowchart of the present study is shown in Figure 1. We selected 96 NSCLC cases for the entire analysis. The lncRNAs expressions were detected between 43 RR cases and 53 RS cases. Next, 60 highly expressed lncRNAs and 72 lowly expressed lncRNAs were obtained. Figures 2(a) and 2(b) revealed that radioresponse-related lncRNAs (RRlncRNAs) differed dramatically between RS and RR groups.

(a)

(b)

3.2. Construction of the RRlncRNAs Signature

Firstly, we randomly divided all samples into training cohort (n = 48) and test cohort (n = 48). After applying the univariate Cox method on the training set, we screened 26 representative prognostic lncRNAs significantly correlated with survival (Figure 3(a)). Then, LASSO-penalized regression was applied to diminish the overfitting of signature (Figures 3(b) and 3(c)). Moreover, a signature of four radioresponse-related lncRNAs (CASC19, LINC01977, LINC02471, and MAGI2-AS3) was built up by the multivariate Cox analysis (Table 2). The risk value of each NSCLC case was generated using the following formula: [exp CASC19 × (0.3999)] + [exp LINC01977 × (0.3693)] + [exp LINC02471 × (−0.0701)] + [exp MAGI2 − AS3 × (−0.0238)]. Subsequently, all patients were separated into the high-risk score group and low-risk score group by the median risk value. The layout of risk score and OS outcomes of NCSLS patients in two risk subgroups are presented in Figure 3(d).

(a)

(b)

(c)

(d)

(e)

(f)

Figure 3

Identification of the radioresponse-associated lncRNAs biomarker model in the training cohort. (a) Univariate Cox analysis for determination of prognostic radioresponse-associated lncRNAs. (b) LASSO regression method for construction of radioresponse-related lncRNAs signature. (c) LASSO coefficient of 26 radioresponse-related lncRNAs in NSCLC. (d) Distribution of increasing risk scores based on radioresponse-related lncRNAs (upper) and clinical outcome of NSCLC cases with increasing risk scores (below). (e) Kaplan–Meier curves of NSCLC cases were classified by a median value of risk. (f) ROC methods reveal the potentiality of employing the risk score to forecast the prognosis of NSCLC patients.

In the training cohort, the KM survival curves indicated a notable difference in survival outcome between the low-risk group and high-risk group (Figure 3(e)). ROC plots were applied to measure the reliability of our constructed model. The area under the ROC (AUC) values of 1-, 3-, and 5-year OS rates were 0.805, 0.904, and 0.903, pointing out the powerful predictive ability of the signature (Figure 3(f)). Similarly, the survival outcome of cases could be distinguished well by the model through KM survival curves in both the test set and entire cohort (Figures 4(a) and 4(b)). Meanwhile, the value of AUC for the test cohort and entire set is shown in Figures 4(c) and 4(d). In addition, risk score layout and OS outcome in the test set and entire cohort are shown in Figures 4(e) and 4(f).

(a)

(b)

(c)

(d)

(e)

(f)

3.3. Development of the Radioresponse-Based Nomogram

To determine the independence of the four-RRlncRNA signature, we carried out univariate and multivariate Cox methods to analyze the clinical characteristics and risk scores generated by the model. As shown in Figure 5(a), stage (), M stage (), N stage (), and risk score () were meaningful for predicting survival outcome. Multivariate Cox analysis disclosed that only risk score () was remarkably correlated with survival time, suggesting our signature could be independent of the clinical factors (Figure 5(b)). Afterward, a survival prognosis nomogram was established with clinical characteristics and a radioresponse-based model in the entire cohort (Figure 5(c)). At the same time, the calibration curves for the predictive ability of 1, 3, and 5 years presented no declination between the ideal reference line and prediction by prognostic nomogram (Figures 5(d)–5(f)).

(a)

(b)

(c)

(d)

(e)

(f)

3.4. GSEA Enrichment of the Signature

As displayed in Figure 6, the top six hallmark pathways were “p53 pathway,” “IL6/JAK/STAT3 pathway,” “Notch pathway,” “DNA repair,” “glycolysis,” and “apoptosis.” The results indicated that these hallmark pathways related to radioresistance regulation were activated in the patients with the high-risk group.

3.5. The Association between the Clinical Traits and the Signature in NSCLC

To evaluate the clinical potency of our signature, we employed clinical relative analysis. The heat map demonstrated that the high-risk score was remarkably related to the status, stage, gender, age, histological type, and radiotherapy response (Figure 7(a)).

(a)

(b)

(c)

3.6. Immune Activity Analysis of the Signature

Immune microenvironment has been proved to be closely bound up with the regulation of tumor progression and radioresistance [18, 19]. To investigate the power of our model for mirroring the condition of the immune microenvironment in NSCLC, we performed immunity analysis by the ssGSEA method. Figure 7(b) shows the difference in immune responses between both risk score groups. We then estimated the association between risk score and immunocyte subpopulations infiltration by ssGSEA of TCGA-NSCLC. As uncovered by (Figure 7(c)), except for B cells, CD8 T cells, and pCDs, the remaining immunocyte populations were highly enriched in the high-risk group.

3.7. Construction and Functional Analysis of the Radioresponse-Related ceRNA Network

Currently, a growing body of evidence indicates that lncRNAs could bind with microRNAs (miRNAs) to regulate the expression s of target mRNAs, which may modulate the biological process of malignant tumors [20–22]. We set up a ceRNA regulatory network consisting of four lncRNAs, 15 miRNAs, and 257 mRNAs based on interaction relationships predicted by DIAND, miRDB, TargetScan, and miRTarBase (Figure 8). Moreover, GO and KEGG analyses were implemented to better unearth the underlying function of the network. Figure 9(a) illustrates that the most distinctively enriched biological process contained “cell growth,” “positive regulation of secretion,” “ERK1 and ERK2 cascade,” “response to oxidative stress,” “cell junction assembly,” and “DNA damage stimulus.” In terms of the KEGG pathway, the main important pathways included “MAPK pathway,” “Hippo pathway,” and “Wnt pathway” and “p53 pathway” (Figure 9(b)).

(a)

(b)

3.8. Determination of the Expression Patterns and the Prognostic Performance of the Four-lncRNA Signature

As revealed by Figure 10(a), CASC19 and LINC01977 were greatly upregulated in NSCLC carcinoma than in normal tissues, whereas LINC02471 and MAGI2-AS3 were downregulated in normal tissues. In addition, CASC19 and LINC01977 were highly expressed in the RR group, while LINC02471 and MAGI2-AS3 had lower expressions in the RR group (Figure 10(b)). Based on the dataset of TCGA-NSCLC patients with radiotherapy history, KM survival analysis indicated that LINC01977 was accompanied with dismal OS in NSCLC patients with RT (Figure 10(c)).

(a)

(b)

(c)

3.9. Inhibition of LINC01977 Enhanced Radiosensitivity in NSCLC

Firstly, we applied the qRT-PCR assay to confirm the expression level of four lncRNAs in BEAS-2B, H520, and A549 (Figure 11(a)). Subsequently, we selected LINC01977 for further in vitro experiment. As observed from Figure 11(b), LINC01977 expression evidently downregulated within H520 and A549 cells by si-LINC01977 transfection. By performing a clonogenic survival assay, the roles of TRPM2-AS on the radiosensitivity of NSCLC cells were detected. After being exposed to an increasing dose of irradiation (IR), colony survival fractions were dramatically reduced by LINC01977 downregulation (Figure 11(c)). The same results were acquired by DNA damage assay. We found that NSCLC cells with si-LINC01977 greatly promoted the expression of γH2AX (DNA damage marker) after receiving an 8 Gy dose of RT (Figure 11(d)).

(a)

(b)

(c)

(d)

4. Discussion

Non-small-cell lung cancer (NSCLC) is a primary thoracic tumor with a terrifically malignant behavior, which is the leading cause of cancerous death. Despite the considerable success of radiotherapy in the treatment of NSCLC, radioresistance is still considered a major obstacle to effective radiotherapy. Therefore, it is highly essential to unearth novel biomarkers predicting radiation response, to obtain new therapeutic targets and develop new treatment options. Accumulated evidence has reported that lncRNAs act as regulators in NSCLC radioresistance [23]. Recently, researchers have attached more attention to lncRNA-based signatures because of their powerful predictive ability compared to TNM staging systems or standard benchmarks [24–26]. However, prognostic markers based on radioresponse associated lncRNA expression profiling have not been analyzed in NSCLC cases with RT history.

Here, according to the TCGA-NSCLC project, we created a newly radioresponse-related lncRNAs signature that could accurately determine high-risk patients and exploit the clinical utility in NSCLC patients. We first collected 132 differentially expressed radioresponse-related lncRNAs based on 96 NSCLC patients who received RT. Then, four radioresponse-related lncRNAs were screened by Cox relative regression methods in the training set. These four lncRNAs were used to build up a risk model for forecasting the survival outcome of NSCLC samples. KM survival method demonstrated a notable difference in OS between the two risk subgroups. Both the test cohort and the entire cohort were employed to confirm the above results. To simplify our proposed model, we created a radiobiological nomogram by combining the RRlncRNAs signature and other clinical factors. The calibration plots indicate that our model exhibits a favorable fit and better clinical effectiveness.

In order to exploit the clinical potency of our radiobiological model, we analyzed the relationship between clinical traits and the signature and observed that low-risk score negatively associated with stage, age, and radiotherapy response. We further explored the lncRNA-related ceRNA network to elucidate the underlying biological function and potential pathway of the selected four lncRNAs. There were many radioresistance associated signal pathways about ceRNA network regulation, such as MAPK pathway,” “Hippo pathway,” and “Wnt pathway” and “p53 pathway”.

To explore the possible biological function of the signature, the GSEA method was employed. GSEA revealed that the gene sets in the high-risk score group were enriched in “DNA repair” “glycolysis,” “apoptosis,” and pathways related to cancer, including p53 pathway, JAK-STAT3 signal pathway, and Notch signal pathway, providing strong evidence that these Hallmark pathways play a central part in in the resistance to radiotherapy in NSCLC. The main mechanism of irradiation-induced tumor cell death is DNA damage. However, cancer cells could correct DNA double-strand breaks by activating DNA damage repair, subsequently promoting radioresistance and tumor cell survival [27]. Glucose uptake supplies energy and biosynthetic materials for tumor cell proliferation [28]. Overexpression of glucose transporter proteins and glycolytic enzymes are common oncogenic signals [29, 30]. Currently, several studies have shown that the radiosensitivity of NSCLC can be enhanced by modulating glycolysis [31]. It has been shown that there is a large amount of P53-dependent apoptosis in irradiation-sensitive tissues, but p53-deficient proved notable radioresistance in the mouse model [32, 33]. Targeting the P53 signaling pathway could enhance radiosensitivity, which has been demonstrated in a variety of cancers [34–36]. Notch signaling mainly regulates the cell cycle, inhibits apoptosis by suppressing PTEN expression, and triggers tumor progression together with PI3K-AKT signaling [37]. Activation of the Notch pathway in NSCLC patients is closely bound up with a worse prognosis, and cancer cells with high Notch expression are more resistant to irradiation [38].

Radiotherapy-induced immunosuppressive changes in the TME are conducive to the immune escape of cancer cells, which in turn confers radioresistance in tumors [18]. Thus, we evaluate the condition of the NSCLC immune microenvironment by the ssGSEA algorithm. As revealed by immunocyte infiltration analysis, patients with high risk showed dramatically higher proportions of M2-like macrophages. M2-like macrophage, as a protumor subtype of tumor-associated macrophage (TAM), shows radioresistance activity in TME. Apoptotic cells, appearing after radiotherapy, activate macrophages with the M2 phenotype to secrete a range of cytokines such as TGF-B which could induce radioresistance and block the development of ”tumor vaccine” [39–41]. Zhang et al. revealed that M2 macrophage-derived exosome could heighten resistance to RT through transferring AGAP2-AS1 into lung cancer cells [42]. Other immunosuppressive immune Tregs were also highly infiltrated in a high-risk group. The existence of Tregs might curb the effectiveness of RT. With the increased recruitment level of Tregs after irradiation, the immunomodulatory effects induced by RT were greatly reduced [43]. In addition, Treg is more resistant to irradiation than other T cell populations due to its increased expression of Akt [44].

In the present study, our constructed risk signature consisted of four lncRNAs which were closely related to NSCLC prognosis. Among these five lncRNAs, CASC19 and LINC01977 are potential risky indicators, but LINC02471 and MAGI2-AS3 are potential protective indicators. CASC19 (Cancer Susceptibility 19) has been reported to mediate the progression and radioresistance of several tumors. In pancreatic cancer (PC), CASC19 could facilitate malignant growth and metastasis of PC via miR-148b/E2F7 ceRNA axis with sponge activity [45]. As suggested by Liu et al., CASC19 was highly expressed in radioresistant nasopharyngeal carcinoma (NPC) cells and could confer radioresistance in NPC through autophagy-related AMPK-mTOR pathway [46]. In addition, CASC19 also has an oncogenic effect on NSCLC. Wang et al. found that CASC19 could bind with miR-301b-3p to regulate the expression of low-density lipoprotein receptor (LDLR), which promotes NSCLC cell proliferation and metastasis [47]. The LINC01977 is a novel carcinogenic promoter reported in breast cancer (BC). Silencing LINC01977 could greatly repress BC development and chemoresistance to doxorubicin. In terms of molecular mechanism, LINC01977 might sponge miR-212-3p to attenuate the expression of the GOLM1 gene [48]. Notably, LINC01977 has not been studied in lung cancer, which suggests that we can do further investigation. LINC02471, a member of the intergenic lncRNA family, is recognized as a crucial enhancer in papillary thyroid carcinoma (PTC). Chen et al. observed that downregulation of LINC02471 could hinder PTC cell metastasis and trigger l apoptosis by binding with miR-375, indicating that it can be considered as a potent predictor for metastasis of patients with PTC [49]. MAGI2-AS3 was originally identified as an antisense RNA of Membrane Associated Guanylate Kinase 2 (MAGI2) involved in regulating the malignant behaviors of various types of cancers. As discovered by Cheng et al., MAGI2-AS3 could enhance sensitivity to radiotherapy in esophageal cancer by inhibiting the HOXB7 gene, offering a valuable molecular marker for radioresistance [50]. Another research studied that bladder cancer (BC) patients with upregulation of MAGI2-AS3 display a lower incidence of cancer metastasis. MAGI2-AS3 may reinforce the stability of MAGI2 to mediate epithelial-mesenchymal transition (EMT) of BC [51]. Our results are in line with these studies, suggesting the MAGI2-AS3 is a protective factor (HR < 1) in NSCLC.

Finally, we explored the relationship between LINC01977 expression and radioresistance in NSCLC cell lines. We found that silencing LINC01977 significantly suppressed clonogenic survival and boosted DNA damage repair, implying that LINC01977 could be a possible therapeutic target for promoting radiotherapy sensitivity in NSCLC patients.

There are some limitations to our study. First, all research populations processed in this subject were merely collected from the TCGA database. The large-scale clinical cohort data or other external datasets need to be warranted to confirm our signature. Moreover, considering our in vitro assays mainly focus on the cell phenotype of tumor, molecular mechanism and in vivo experiments are need to verify the results in further exploration.

In summary, we firstly created a radioresponse-related lncRNAs signature based on four prognostic lncRNAs associated with radiation therapy which offer a reliable reference for the prognostic forecast. Our subject brings new sights into the clinical strategy of NSCLC patients who received radiotherapy.

Data Availability

The public datasets to support the results of this subject can be gained from TCGA (https://portal.gdc.cancer.gov/).

Conflicts of Interest

All the authors declare that they have no conflicts of interest.

Authors’ Contributions

J. S., S. Z., and Y. S. contributed equally to this work. Q. T. and J. S. designed the original study. J. S., S. Z., and Z. Y. collected and analyzed the data. J. S., S. Z., Y. S., and J. G. explained the data. J. S. and S. Z. drafted the manuscript. Q. T. revised the manuscript. All authors read and approved the final manuscript.

Acknowledgments

The authors would like to thank the Core Facility of Jiangsu Provincial People’s Hospital for their help with the detection of experimental samples. This work was supported by the National Natural Science Foundation of China (82073344) and the Jiangsu Province Special Program for Young Medical Talent (Grant no. QNRC2016571).

Supplementary Materials

Supplementary Table 1: target sequence of si-NC and si-LINC01977. Supplementary Table 2: primers utilized for qRT-PCR. (Supplementary Materials)

References

M. Evison and U. K. L. AstraZeneca, “The current treatment landscape in the UK for stage III NSCLC,” British Journal of Cancer, vol. 123, pp. 3–9, 2020.
View at: Publisher Site | Google Scholar
E. N. Imyanitov, A. G. Iyevleva, and E. V. Levchenko, “Molecular testing and targeted therapy for non-small cell lung cancer: current status and perspectives,” Critical Reviews in Oncology, vol. 157, Article ID 103194, 2021.
View at: Publisher Site | Google Scholar
N. Giaj-Levra, P. Borghetti, A. Bruni et al., “Current radiotherapy techniques in NSCLC: challenges and potential solutions,” Expert Review of Anticancer Therapy, vol. 20, no. 5, pp. 387–402, 2020.
View at: Publisher Site | Google Scholar
E. C. Ko, D. Raben, and S. C. Formenti, “The integration of radiotherapy with immunotherapy for the treatment of non-small cell lung cancer,” Clinical Cancer Research, vol. 24, no. 23, pp. 5792–5806, 2018.
View at: Publisher Site | Google Scholar
Z. Liao, J. J. Lee, R. Komaki et al., “Bayesian adaptive randomization trial of passive scattering proton therapy and intensity-modulated photon radiotherapy for locally advanced non-small-cell lung cancer,” Journal of Clinical Oncology, vol. 36, no. 18, pp. 1813–1822, 2018.
View at: Publisher Site | Google Scholar
A. Merlotti, A. Bruni, P. Borghetti et al., “Sequential chemo-hypofractionated RT versus concurrent standard CRT for locally advanced NSCLC: GRADE recommendation by the Italian Association of Radiotherapy and Clinical Oncology (AIRO),” La radiologia medica, vol. 126, no. 8, pp. 1117–1128, 2021.
View at: Publisher Site | Google Scholar
C. P. Ponting, P. L. Oliver, and W. Reik, “Evolution and functions of long noncoding RNAs,” Cell, vol. 136, no. 4, pp. 629–641, 2009.
View at: Publisher Site | Google Scholar
M. C. Jiang, J. J. Ni, W. Y. Cui, B. Y. Wang, and W. Zhuo, “Emerging roles of lncRNA in cancer and therapeutic opportunities,” American journal of cancer research, vol. 9, no. 7, pp. 1354–1366, 2019.
View at: Google Scholar
E. S. Martens-Uzunova, R. Böttcher, C. M. Croce, G. Jenster, T. Visakorpi, and G. A. Calin, “Long noncoding RNA in prostate, bladder, and kidney cancer,” European Urology, vol. 65, no. 6, pp. 1140–1151, 2014.
View at: Publisher Site | Google Scholar
A. Bhan, M. Soleimani, and S. S. Mandal, “Long noncoding RNA and cancer: a new paradigm,” Cancer Research, vol. 77, no. 15, pp. 3965–3981, 2017.
View at: Publisher Site | Google Scholar
M. K. Iyer, Y. S. Niknafs, R. Malik et al., “The landscape of long noncoding RNAs in the human transcriptome,” Nature Genetics, vol. 47, no. 3, pp. 199–208, 2015.
View at: Publisher Site | Google Scholar
H. He, X. Song, Z. Yang et al., “Upregulation of KCNQ1OT1 promotes resistance to stereotactic body radiotherapy in lung adenocarcinoma by inducing ATG5/ATG12-mediated autophagy via miR-372-3p,” Cell Death & Disease, vol. 11, no. 10, p. 883, 2020.
View at: Publisher Site | Google Scholar
Y. Liu, X. Chen, X. Chen et al., “Long non-coding RNA HOTAIR knockdown enhances radiosensitivity through regulating microRNA-93/ATG12 axis in colorectal cancer,” Cell Death & Disease, vol. 11, no. 3, p. 175, 2020.
View at: Publisher Site | Google Scholar
M. E. Ritchie, B. Phipson, D. Wu et al., “Limma powers differential expression analyses for RNA-sequencing and microarray studies,” Nucleic Acids Research, vol. 43, no. 7, p. e47, 2015.
View at: Publisher Site | Google Scholar
A. Iasonos, D. Schrag, G. V. Raj, and K. S. Panageas, “How to build and interpret a nomogram for cancer prognosis,” Journal of Clinical Oncology, vol. 26, no. 8, pp. 1364–1370, 2008.
View at: Publisher Site | Google Scholar
M. Yi, D. V. Nissley, F. McCormick, and R. M. Stephens, “ssGSEA score-based Ras dependency indexes derived from gene expression data reveal potential Ras addiction mechanisms with possible clinical implications,” Scientific Reports, vol. 10, no. 1, p. 10258, 2020.
View at: Publisher Site | Google Scholar
F. Chen, B. Xu, J. Li et al., “Hypoxic tumour cell-derived exosomal miR-340-5p promotes radioresistance of oesophageal squamous cell carcinoma via KLF10,” Journal of Experimental & Clinical Cancer Research, vol. 40, no. 1, p. 38, 2021.
View at: Publisher Site | Google Scholar
H. E. Barker, J. T. E. Paget, A. A. Khan, and K. J. Harrington, “The tumour microenvironment after radiotherapy: mechanisms of resistance and recurrence,” Nature Reviews Cancer, vol. 15, no. 7, pp. 409–425, 2015.
View at: Publisher Site | Google Scholar
S. Gandhi and S. Chandna, “Radiation-induced inflammatory cascade and its reverberating crosstalks as potential cause of post-radiotherapy second malignancies,” Cancer and Metastasis Reviews, vol. 36, no. 2, pp. 375–393, 2017.
View at: Publisher Site | Google Scholar
X. Lin, S. Zhuang, X. Chen et al., “LncRNA ITGB8-AS1 functions as a ceRNA to promote colorectal cancer growth and migration through integrin-mediated focal adhesion signaling,” Molecular Therapy, vol. 454, 2021.
View at: Publisher Site | Google Scholar
D. Li, J. She, X. Hu, M. Zhang, R. Sun, and S. Qin, “The ELF3-regulated lncRNA UBE2CP3 is over-stabilized by RNA-RNA interactions and drives gastric cancer metastasis via miR-138-5p/ITGA2 axis,” Oncogene, vol. 40, 2021.
View at: Publisher Site | Google Scholar
H. Li, Q. Mu, G. Zhang et al., “Linc00426 accelerates lung adenocarcinoma progression by regulating miR-455-5p as a molecular sponge,” Cell Death & Disease, vol. 11, no. 12, p. 1051, 2020.
View at: Publisher Site | Google Scholar
X. Zhang, K. Xie, H. Zhou et al., “Role of non-coding RNAs and RNA modifiers in cancer therapy resistance,” Molecular Cancer, vol. 19, no. 1, p. 47, 2020.
View at: Publisher Site | Google Scholar
Y. Lu, W. Wang, Z. Liu, J. Ma, X. Zhou, and W. Fu, “Long non-coding RNA profile study identifies a metabolism-related signature for colorectal cancer,” Molecular Medicine, vol. 27, no. 1, p. 83, 2021.
View at: Publisher Site | Google Scholar
W. Ma, Y. Yao, G. Xu et al., “Identification of a seven-long non-coding RNA signature associated with Jab1/CSN5 in predicting hepatocellular carcinoma,” Cell Death Discovery, vol. 7, no. 1, p. 178, 2021.
View at: Publisher Site | Google Scholar
W. Hong, L. Liang, Y. Gu et al., “Immune-related lncRNA to construct novel signature and predict the immune landscape of human hepatocellular carcinoma,” Molecular Therapy - Nucleic Acids, vol. 22, pp. 937–947, 2020.
View at: Publisher Site | Google Scholar
M. A. Morgan and T. S. Lawrence, “Molecular pathways: overcoming radiation resistance by targeting DNA damage response pathways,” Clinical Cancer Research, vol. 21, no. 13, pp. 2898–2904, 2015.
View at: Publisher Site | Google Scholar
D. Hanahan and R. A. Weinberg, “Hallmarks of cancer: the next generation,” Cell, vol. 144, no. 5, pp. 646–674, 2011.
View at: Publisher Site | Google Scholar
M. G. Vander Heiden and R. J. DeBerardinis, “Understanding the intersections between metabolism and cancer biology,” Cell, vol. 168, no. 4, pp. 657–669, 2017.
View at: Publisher Site | Google Scholar
S. Ganapathy-Kanniappan and J.-F. H. Geschwind, “Tumor glycolysis as a target for cancer therapy: progress and prospects,” Molecular Cancer, vol. 12, no. 1, p. 152, 2013.
View at: Publisher Site | Google Scholar
T. Gong, L. Cui, H. Wang, H. Wang, and N. Han, “Knockdown of KLF5 suppresses hypoxia-induced resistance to cisplatin in NSCLC cells by regulating HIF-1α-dependent glycolysis through inactivation of the PI3K/Akt/mTOR pathway,” Journal of Translational Medicine, vol. 16, no. 1, p. 164, 2018.
View at: Publisher Site | Google Scholar
J. Zhang, L. Shen, and L.-Q. Sun, “The regulation of radiosensitivity by p53 and its acetylation,” Cancer Letters, vol. 363, no. 2, pp. 108–118, 2015.
View at: Publisher Site | Google Scholar
S. W. Lowe, E. M. Schmitt, S. W. Smith, B. A. Osborne, and T. Jacks, “p53 is required for radiation-induced apoptosis in mouse thymocytes,” Nature, vol. 362, no. 6423, pp. 847–849, 1993.
View at: Publisher Site | Google Scholar
K. B. Leszczynska, I. P. Foskolou, A. G. Abraham et al., “Hypoxia-induced p53 modulates both apoptosis and radiosensitivity via AKT,” Journal of Clinical Investigation, vol. 125, no. 6, pp. 2385–2398, 2015.
View at: Publisher Site | Google Scholar
Z. Wang, J.-W. Mao, G.-Y. Liu et al., “MicroRNA-372 enhances radiosensitivity while inhibiting cell invasion and metastasis in nasopharyngeal carcinoma through activating the PBK-dependent p53 signaling pathway,” Cancer Medicine, vol. 8, no. 2, pp. 712–728, 2019.
View at: Publisher Site | Google Scholar
T. He, J. Guo, H. Song et al., “Nutlin-3, an antagonist of MDM2, enhances the radiosensitivity of esophageal squamous cancer with wild-type p53,” Pathology and Oncology Research, vol. 24, no. 1, pp. 75–81, 2018.
View at: Publisher Site | Google Scholar
K. Thippu Jayaprakash and A. Michael, “Notch inhibition: a promising strategy to improve radiosensitivity and curability of radiotherapy,” Clinical Oncology, vol. 33, no. 1, pp. e44–e49, 2021.
View at: Publisher Site | Google Scholar
J. Theys, S. Yahyanejad, R. Habets et al., “High NOTCH activity induces radiation resistance in non small cell lung cancer,” Radiotherapy & Oncology, vol. 108, no. 3, pp. 440–445, 2013.
View at: Publisher Site | Google Scholar
C. Vanpouille-Box, J. M. Diamond, K. A. Pilones et al., “TGFβ is a master regulator of radiation therapy-induced antitumor immunity,” Cancer Research, vol. 75, no. 11, pp. 2232–2242, 2015.
View at: Publisher Site | Google Scholar
M. Shevtsov, H. Sato, G. Multhoff, and A. Shibata, “Novel approaches to improve the efficacy of immuno-radiotherapy,” Frontiers in Oncology, vol. 9, p. 156, 2019.
View at: Publisher Site | Google Scholar
M. Jarosz-Biej, R. Smolarczyk, T. Cichon, and N. Kulach, “Tumor microenvironment as A “game changer” in cancer radiotherapy,” International Journal of Molecular Sciences, vol. 20, no. 13, 2019.
View at: Publisher Site | Google Scholar
F. Zhang, Y. Sang, D. Chen et al., “M2 macrophage-derived exosomal long non-coding RNA AGAP2-AS1 enhances radiotherapy immunity in lung cancer by reducing microRNA-296 and elevating NOTCH2,” Cell Death & Disease, vol. 12, no. 5, p. 467, 2021.
View at: Publisher Site | Google Scholar
F. G. Herrera, J. Bourhis, and G. Coukos, “Radiotherapy combination opportunities leveraging immunity for the next oncology practice,” CA: A Cancer Journal for Clinicians, vol. 67, no. 1, pp. 65–85, 2017.
View at: Publisher Site | Google Scholar
E. Persa, A. Balogh, G. Sáfrány, and K. Lumniczky, “The effect of ionizing radiation on regulatory T cells in health and disease,” Cancer Letters, vol. 368, no. 2, pp. 252–261, 2015.
View at: Publisher Site | Google Scholar
T. Lu, G. H. Wei, J. Wang, and J. Shen, “LncRNA CASC19 contributed to the progression of pancreatic cancer through modulating miR-148b/E2F7 axis,” European Review for Medical and Pharmacological Sciences, vol. 24, no. 20, pp. 10462–10471, 2020.
View at: Publisher Site | Google Scholar
H. Liu, W. Zheng, Q. Chen et al., “lncRNA CASC19 contributes to radioresistance of nasopharyngeal carcinoma by promoting autophagy via AMPK-mTOR pathway,” International Journal of Molecular Sciences, vol. 22, no. 3, 2021.
View at: Publisher Site | Google Scholar
L. Wang, C. Lin, N. Sun, Q. Wang, X. Ding, and Y. Sun, “Long non-coding RNA CASC19 facilitates non-small cell lung cancer cell proliferation and metastasis by targeting the miR-301b-3p/LDLR axis,” The Journal of Gene Medicine, vol. 22, no. 12, Article ID e3254, 2020.
View at: Publisher Site | Google Scholar
Z. Li, Y. Li, X. Wang et al., “LINC01977 Promotes breast cancer progression and chemoresistance to doxorubicin by targeting miR-212-3p/GOLM1 Axis,” Frontiers in Oncology, vol. 11, Article ID 657094, 2021.
View at: Google Scholar
D. Chen, Z. Huang, Y. Ning, and C. Lou, “Knockdown of LINC02471 inhibits papillary thyroid carcinoma cell invasion and metastasis by targeting miR-375,” Cancer Management and Research, vol. 12, pp. 8757–8771, 2020.
View at: Publisher Site | Google Scholar
W. Cheng, X. Shi, M. Lin, Q. Yao, J. Ma, and J. Li, “LncRNA MAGI2-AS3 overexpression sensitizes esophageal cancer cells to irradiation through down-regulation of HOXB7 via EZH2,” Frontiers in Cell and Developmental Biology, vol. 8, Article ID 552822, 2020.
View at: Publisher Site | Google Scholar
D. Shen, J. Xu, X. Cao, X. Cao, H. Tan, and H. Deng, “Long noncoding RNA MAGI2-AS3 inhibits bladder cancer progression through MAGI2/PTEN/epithelial-mesenchymal transition (EMT) axis,” Cancer Biomarkers, vol. 30, no. 2, pp. 155–165, 2021.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Jiahang Song 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

410

Downloads

664

Citations