- 1Department of Oncology, The First Hospital of Lanzhou University, Lanzhou, China
- 2The First School of Clinical Medicine, Lanzhou University, Lanzhou, China
- 3Department of Oncology, The First People’s Hospital of Lanzhou, Lanzhou, China
Objective: To establish a prediction model based on autophagy-related lncRNAs and investigate the functional enrichment of autophagy-related lncRNAs in colorectal cancer.
Methods: TCGA database was used to extract the transcriptome data and clinical features of colorectal cancer patients. HADb was used to obtain autophagy-related genes. Pearson correlation analysis was performed to identify autophagy-related lncRNAs. The autophagy-related lncRNAs with prognostic values were selected. Based on the selected lncRNAs, the risk score model and nomogram were constructed, respectively. Calibration curve, concordance index, and ROC curve were performed to evaluate the predictive efficacy of the prediction model. GSEA was performed to figure out the functional enrichment of autophagy-related lncRNAs.
Results: A total of 13413 lncRNAs and 938 autophagy-related genes were obtained. A total of 709 autophagy-related genes were identified in colon cancer tissues, and 11 autophagy-related lncRNAs (AL138756.1, LINC01063, CD27-AS1, LINC00957, EIF3J-DT, LINC02474, SNHG16, AC105219.1, AC068580.3, LINC02381, and LINC01011) were finally selected and set as prognosis-related lncRNAs. According to the risk score, patients were divided into the high-risk and low-risk groups, respectively. The survival K–M (Kaplan–Meier) curve showed the low-risk group exhibits better overall survival than the high-risk group. The AUCs under the ROC curves were 0.72, 0.814, and 0.83 at 1, 3, and 5 years, respectively. The C-index (concordance index) of the model was 0.814. The calibration curves at 1, 3, and 5 years showed the predicting values were consistent with the actual values. Functional enrichment analysis showed that autophagy-related lncRNAs were enriched in several pathways.
Conclusions: A total of 11 specific autophagy-related lncRNAs were identified to own prognostic value in colon cancer. The predicting model based on the lncRNAs and clinical features can effectively predict the OS. Furthermore, functional enrichment analysis showed that autophagy-related genes were enriched in various biological pathways.
Introduction
Colorectal cancer (CRC) is considered to be one of the leading causes of cancer death, and with the increasing incidence of colorectal cancer, its degree of harm has received increasing attention (Carr, 2018). The epidemiological investigation showed that CRC has been set as the third most incident cancer and the second most common in mortality among malignancies (Siegel et al., 2020). A previous study showed that the 5-year survival rate of CRC patients was merely 10–50%, which strongly affects the living quality (Miller et al., 2019). The histopathological types of CRC include adenocarcinoma, squamous cell carcinoma, and mucinous carcinoma. Among them, adenocarcinoma is the most common type, which occupied about 95% of CRC patients (Barresi et al., 2019).
The previous predicting analysis showed that CRC was associated with various clinical and environmental features including alcohol abuse, smoking, and different regions and countries (Øines et al., 2017; Nana and Edward, 2019). And nowadays, the tumor-node-metastasis (TNM) stage is set as the most common tool to predict the risk of patients suffering from malignancies (Weiser, 2018). However, the TNM stage predicting tool showed limited effects in prognostic prediction. With the development of transcriptome technology, a large number of sequencing results have been used to evaluate and predict the prognosis of cancer patients (Wang et al., 2020). Previous studies showed that autophagy-related genes and long non-coding RNA (lncRNAs), immune-related genes, and lncRNAs own preferable prognostic predicting effect on CRC (Wiener et al., 2014; Tokunaga et al., 2020). However, few studies investigate the combined predicting effects of sequencing data and clinical features.
Autophagy, a programmed cellular process, mainly contributes to the degradation and recycling of damaged cellular organelles and macromolecules (Kimmelman and White, 2017; Li et al., 2020; Ashrafizadeh et al., 2022). Dysfunction of the autophagic process has been associated with the onset and development of many human chronic pathologies, such as cardiovascular, metabolic, and neurodegenerative diseases as well as cancer (Francesca et al., 2018). Therefore, a large number of researchers tried to determine the regulation mechanism of autophagy (Lin et al., 2021; Liu et al., 2021; Zhou et al., 2021). Autophagy plays an important role in maintaining homeostasis. The process is controlled by autophagy-related genes (Zhou, Dong, 2021) and autophagy plays an important role in maintaining the homeostasis of the body. Previous studies showed that the autophagy process was already verified to be involved in the progress of cancer (Francesca, Sadia, 2018). The decreased expression of autophagy-related genes was associated with the acceleration of the cancer progress (Burada et al., 2015). And based on the specific function of autophagy, autophagy can be set as a potential prognostic predicting strategy in different types of tumors. Previous studies showed that autophagy-related genes can be an effective way to predict the prognosis of CRC (Thongchot et al., 2018).
Long non-coding RNAs (lncRNAs) were defined as RNA links with more than 200 nucleotides. The main function of lncRNAs was set to regulate the gene expression and protein synthesis (Liu et al., 2019). Also, lncRNAs are associated with cell proliferation, differentiation, and microRNA regulation (Jin et al., 2019). Previous studies showed that lncRNA participated in the regulation of cell autophagy by regulating the function of autophagy-related genes and proteins (Ulitsky and Bartel, 2013; Ashrafizaveh et al., 2021; Mirzaei et al., 2021; Mirzaei et al., 2022). HAGLROS as a lncRNA with 699bp can inhibit the autophagy process by regulating the mTOR signaling pathway (Frankel et al., 2017).
Therefore, we performed this study aiming to establish a prognostic model based on autophagy-related lncRNAs and clinical features in CRC patients.
Materials and Methods
Patient Datasets
Figure 1 illustrates the whole process of autophagy-related lncRNA CRC prediction model establishing. CRC datasets were downloaded from The Cancer Genome Atlas database (TCGA, https://cancergenome.nih.gov/), which includes clinical features, mRNA expression profiles, and lncRNA profiles. CRC patients whose overall survival (OS) time was less than 30 days or whose survival status was unknown were excluded from this study. Clinical features containing survival time, survival status, age, sex, grade, and TNM stage were included for further statistical analysis. All data were obtained from the open access TCGA database and therefore did not require the medical ethics committee approval.
Autophagy-related genes were downloaded from the Human Autophagy database (HADb, http://autophagy.lu/clustering/index. html). The “limma” package in R software was used to extract autophagy gene data from TCGA-CRC mRNA expression profiles. Pearson correlation analysis was performed to determine the autophagy-related lncRNAs. In this study, the correlation coefficient |R2| >0.3 and p < 0.001 (Wu et al., 2021) were defined as autophagy-related lncRNAs.
Univariate COX Regression and Least Absolute Shrinkage and Selection Operator Regression Analysis
As Figure 1 shows, prognostic autophagy-related lncRNAs were confirmed by univariate COX regression analysis and absolute shrinkage and selection operator (LASSO) regression. Firstly, autophagy-related lncRNAs in univariate COX regression analysis whose p values were less than 0.05 were selected. Then, the selected lncRNAs were set into LASSO regression analysis.
Multivariate COX Regression Analysis and Prognostic Model Establishment
All selected lncRNAs from LASSO regression were included in the multivariate COX regression model to generate their coefficient. Then, a prognosis prediction model was established with risk scores = ΣCoef * exp (genes). The risk score model was established based on coefficient values and expression levels. In this study, the median risk score was set as the cut-off value. And according to the cut-off value, patients were divided into high-risk groups and low-risk groups. The OS difference between groups was calculated by Kaplan-Meier (K-M) analysis.
Based on the clinical features and risk score of each colon cancer patient, we set univariate and multivariate COX regression analysis to construct the prognostic model. The nomogram was set to predict the 1-, 3-, and 5-years overall survival (OS) of each patient.
Nomogram and Prognostic Model Evaluation
The predictive efficacy of the prediction model was evaluated by calculating the area under the ROC curve (AUC). And what is more, we also calculated the concordance index (C-index), a calibration curve to evaluate the predicted efficacy of the mixed model.
Co-Expression Network and Gene Functional Enrichment Analysis
Based on the 11 prognosis-related lncRNAs identified, a co-expression network connecting these autophagy-related lncRNAs and autophagy-related genes was established. Protein-protein network (PPI) autophagy-related lncRNAs and their target genes were established based on Cytoscape 3.7.2 software (Cytoscape Consortium, San Diego, CA, United States). Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis between the high-risk group and low-risk group was performed by Gene Set Enrichment Analysis (GSEA, http://www.broadinstitute. org/gsea/index.jsp). In total,1000 genome permutations were performed per analysis. The enrichment pathway screening conditions were: FDR< 0.25 and NOM p value <0.05.
Statistical Analysis
All statistical analysis procedures were executed based on R (v4.0.3, New Zealand) software. Univariate cox analysis, LASSO regression, and multivariate cox analysis were used sequentially, and combined autophagy-related lncRNAs and clinical features to build prognostic models. Kaplan-Meier was performed to generate the survival curves. Log-rank analysis was utilized to compare the survival difference between the two groups. All statistical significances were set as p value less than 0.05.
Results
Differentially Expressed Autophagy-Related Gene and lncRNA Identification
A total of 417 CRC patients were screened from the TCGA database, the clinical baseline data of CRC patients were shown in Table1. Meanwhile, a total of 13413 lncRNAs were identified from the TCGA-COAD. A total of 938 autophagy-related genes were obtained from HADb. Finally, 709 autophagy-related genes were expressed in colon cancer tissues. Pearson correlation analysis showed that a total of 1342 autophagy-related lncRNAs were identified and selected.
Univariate COX Regression and Least Absolute Shrinkage and Selection Operator Regression Analysis
Univariate COX regression analysis selected 56 autophagy-related lncRNAs with p < 0.05. Eleven autophagy-related lncRNAs were finally selected by LASSO regression at the optimal values by using the 1 standard error (SE) of the minimum criteria (Figures 2A–C). K-M analysis based on 11 autophagy-related genes was shown in Figure 3.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/906900/fgene-13-906900-HTML/image_m/fgene-13-906900-g002.jpg)
FIGURE 2. Autophagy related lncRNAs selection by COX regression and LASSO regression analysis. (A) The lasso regression model and cross validation method were used to screen autophagy-related lncNRAs. (B) Cross validation results. The regression coefficient map of genes in the LASSO model. (C) Cox proportional hazards regression analysis of 11 prognostic genes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/906900/fgene-13-906900-HTML/image_m/fgene-13-906900-g003.jpg)
FIGURE 3. Survival curve of each autophagy related lncRNAs. Kaplan-Meier survival curves of 11 autophagy related lncRNAs were selected with p < 0.05 as the screening criteria. The red plots represent the high expression and the blue plots represent the median/low expression.
Risk Scoring Model Construction
The risk scoring model was established by multivariate COX regression analysis. As we set the median risk score as the cut-off value, colon cancer patients included in our study were divided into high-risk group and low-risk group. The K-M survival curve showed that the difference of survival rate in the two groups was statistically significant (p < 0.01, Figure 4A). Meanwhile, the risk curve and scatter plot confirmed that compared with the high-risk group, the low-risk group had a lower risk factor and mortality, and the difference was significant (Figures 4B,C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/906900/fgene-13-906900-HTML/image_m/fgene-13-906900-g004.jpg)
FIGURE 4. Risk score model evaluation. (A)The K-M curve c omparison between the high-risk group and low-risk group. (B) The Risk curve of high-risk and low-risk group. (C) Patients’ survival status distribution.
The univariate and multivariate COX regression analysis based on the clinical features and risk score showed that the risk score was an independent prognostic factor of CRC patients (Figures 5A,B and Table 2). A time-dependent ROC curve was drawn and the AUCs were 0.72, 0.814, and 0.83 at 1-, 3-, and 5-years, respectively (Figure 5C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/906900/fgene-13-906900-HTML/image_m/fgene-13-906900-g005.jpg)
FIGURE 5. Prognostics predicting model establishing and validation. (A) Univariate COX regression analysis based on the clinical features and risk score. (B) Multivariate COX regression analysis. (C) Time-dependent ROC curve. (D) Nomogram based on multivariate COX regression analysis. (E–G) Calibration curve at 1-, 3- and 5-year.
Prognostic Model Establishing and Evaluation
Nomogram was established based on the risk score, age, gender, and TNM stages (Figure 5D). The results showed that 11 kinds of lincRNA may be the prognostic factors of colorectal cancer patients (all p < 0.05). Calibration curves at 1-, 3-, and 5-years were shown in Figures 5E–G. The AUCs of each factor at 5-year showed that the risk score, stage, and N stage showed certain prediction ability (0.808, 0.730, and 0.716, Figure 6A). A summary of gene ontology involving eleven lncRNAs (Table 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/906900/fgene-13-906900-HTML/image_m/fgene-13-906900-g006.jpg)
FIGURE 6. Autophagy-related lncNRAs and their target genes in CRC. (A) The ROC curve of each factor. (B) Network of autophagy-related lncNRAs and their target genes. (C) Sankey diagram of lncRNAs and linked genes.
Co-Expression Network and Gene Functional Enrichment Analysis
A total of 256 autophagy-related genes were found to be associated with 11 prognosis-related lncRNAs. At the same time, we established a co-expression network based on the relationship between these genes and lncRNAs, with 267 nodes and 359 edges (Figure 6B). The Sankey diagram visually and clearly shows the association between each node and the patient’s prognosis (Figure 6C). Between the high-risk group and low-risk group, autophagy-related lncRNAs were enriched in GO terms: protein localization to the cilium, pseudouridine synthesis, tricarboxylic acid cycle, ciliary base, ciliary plasm, nucleoid, deacetylase activity, intramolecular transferase activity, NAD-dependent protein deacetylase activity, prenyltransferase activity, and KEGG terms: AXON guidance, citrate cycle TCA cycle, glycerophospholipid metabolism, glycosaminoglycan biosynthesis chondroitin sulfate, snare interactions in vesicular transport, terpenoid backbone biosynthesis, VEGF signaling pathway (Figures 7A–D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/906900/fgene-13-906900-HTML/image_m/fgene-13-906900-g007.jpg)
FIGURE 7. Gene functional enrichment analysis of autophagy-related lncNRAs. (A) The significant terms of MF (molecular function) enrichment. (B) The significant terms of CC (cellular component) enrichment. (C) The significant terms of BP (biological process) enrichment. (D) The top 10 significant terms of KEGG analysis.
Discussion
The enhancement of autophagy can make CRC cells survive in the status of nutrition and energy deficiency. The enhanced autophagy also strengthened radiotherapy and chemotherapy resistance of CRC cells. And during biological therapy procedures, physicians tended to inhibit the cell autophagy to mediate the CRC cell programming death (Liu et al., 2018). Furthermore, previous studies have shown that cell autophagy was associated with the CRC prognosis (F, 2018). Based on these backgrounds, the autophagy-related lncRNA signature can be an effective and reliable indicator to predict the prognosis of CRC.
Until now, the specific mechanisms of autophagy in cancer progression have remained unclear, but substantial evidence suggests that autophagy has a great potential as a prognostic predictor in CRC (Lv et al., 2022; Zhu et al., 2022). We obtained the open access data from the TCGA database (TCGA-CRC) and comprehensively analyzed the association between autophagy-related lncRNAs and CRC by bioinformatics analysis. We aimed to screen for features that could be used to predict the CRC prognosis and guide treatment strategies, as these features may be new prognostic markers.
We extract autophagy-related lncRNAs from the TCGA database and HADb database by Pearson correlation analysis. We performed univariate cox analysis, LASSO regression, and multivariate cox, respectively, a total of 11 autophagy-related lncRNAs were finally selected: AL138756.1, LINC01063, CD27-AS1, LINC00957, EIF3J-DT, LINC02474, SNHG16, AC105219.1, AC068580.3, LINC02381, and LINC01011. Among the 11 autophagy-related lncRNAs, 8 lncRNAs were already verified to be associated with cancer development: LINC01063, CD27-AS1, LINC00957, EIF3J-DT, LINC02474, SNHG16, LINC02381, and LINC01011. Soudeh et al. compared breast cancer tissue with normal breast tissue and find that LINC01063 was significantly decreased in breast cancer tissue (Ghafouri-Fard et al., 2021). Anirban et al. concluded that down-regulated expression of CD27-AS1 was associated with the prognosis of cervical carcinoma (Ar et al., 2020). Zhang et al. reported that LINC00957 can be set as a prognostic marker in human CRC (Zhang et al., 2007). Luo et al. concluded that EIF3J-DT regulated the autophagy procedure in gastric cancer by targeting ATG14 (Luo et al., 2021). Du et al. found that LINC02474 mediated the CRC cell’s apoptosis by inhibiting the expression of GZMB (Du et al., 2021). Wu et al. concluded that SNHG16 facilitated the nasopharyngeal carcinoma progression by sponging miR-520a-3p to upregulate MAPK1 expression (Wu et al., 2021). Through the RT-PCR technique, Sun et al reported that LINC02381 enhanced CBX5 expression by binding the CBX5 promoter and activating CBX5 transcription in glioma (Sun et al., 2021). Song et al. reported that LINC01011 controlled mitochondrial fission by inhibiting BRCA1 transcription in squamous cell carcinoma (Fan et al., 2020). The remaining 3 lncRNAs (al138756.1, AC105219.1, and AC068580.3) have not been proved to be correlated with the progress and prognosis in CRC. There still needs further research on these lncRNAs.
Wu et al. had already confirmed that autophagy-related lncRNAs were associated with the poor prognosis of CRC patients (Wu, Yin, 2021). However, the AUCs of the model were 0.70, 0.76, and 0.68 at 1-, 3-, and 5-years, which meant that the predictive value of the model was at a moderate degree and needed to be improved. In our study, according to the risk score formula, patients were divided into the low-risk group and high-risk group by the median risk score. Similar to the previous studies (Fan, Tian, 2020; Yang et al., 2021), the K-M curve showed that low-risk group survival was longer than the high-risk group with significance (p < 0.05). The AUCs at 1-, 3-, and 5-years were 0.72, 0.814, and 0.83, which indicate that the predictive ability of the risk score signatured by autophagy-related lncRNAs were stable and the risk score can be set as an independent prognostic indicator. Based on the previous analysis, the nomogram was built by multivariate COX regression analysis and shown in Figure 5D. As shown in Figures 5A,B, the risk score occupied the largest concentration. C-index, AUCs, and calibration curve showed that the nomogram can be used as an effective tool to predict the prognosis of CRC patients. We were able to conclude that our prognostic model was valid and independent of other clinical factors such as TNM classification, clinical stage, age, and gender.
To further investigate the functions of prognostic autophagy-related lncRNAs, GO enrichment analysis and KEGG pathway analysis were performed between the high-risk group and low-risk group. As for GO terms, the lncRNAs were enriched in the protein localization to the cilium, pseudouridine synthesis, tricarboxylic acid cycle of biological process, ciliary base, ciliary plasm, nucleoid of cellular components and deacetylase activity, intramolecular transferase activity, NAD-dependent protein deacetylase activity, and prenyltransferase activity of molecular function. As for KEGG pathway analysis, lncRNAs were enriched in AXON guidance, citrate cycle TCA cycle, glycerophospholipid metabolism, glycosaminoglycan biosynthesis chondroitin sulfate, snare interactions in vesicular transport, terpenoid backbone biosynthesis, and VEGF signaling pathway. A previous study revealed that the AXON guidance genes were involved in cancer development (Li et al., 2009; Wei et al., 2020). EM et al. revealed that in CRC patients, AXON-related genes such as ROBO1 and ROBO2 were critical genes in cancer pathogenesis (Je et al., 2013). VEGF signaling pathway is tightly associated with tumor angiogenesis. Chao decreased the VEGF expression in CRC tumor tissues and found that tumor growth was significantly suppressed (Fang et al., 2021).
There are several limitations to this study. First, the study type was a retrospective cohort study, which may contain statistical bias. Second, as a retrospective study, the data derived from the public database lacks information including the treatment methods and recurrence time. Third, the sample size was relatively small which may influence the accuracy of the results. Finally, this study was just an analysis based on a public database. Further analyses in vivo and in vitro are needed to verify the conclusion.
Conclusion
All in all, eleven specific autophagy-related lncRNAs were identified to own prognostic value in colon cancer. The predicting model based on lncRNAs and clinical features can effectively predict the OS.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found at: TCGA, https://cancergenome.nih.gov/.
Author Contributions
LD and XH conceived and directed the project and wrote the manuscript. YX and CL performed data bioinformatics analyses. NL helped with part of English writing and checking. All authors read and approved the final manuscript.
Funding
This work was supported by the Hospital Fund Project of the First Hospital of Lanzhou University (ldyyyn 2019-56), the Innovation Fund Project of Colleges and Universities in Gansu Province (2021B-003 and 2021B-007), and the Health Science and Technology Development Project in Lanzhou of Gansu Province (2021016), the Natural Science Foundation of Gansu Province China (21JR7RA377).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Ar, A., Ss, A., Pd, B., Dim, C., Ac, D., Sa, E., et al. (2020). Deregulation of H19 Is Associated with Cervical Carcinoma. Genomics 112 (1), 961–970. doi:10.1016/j.ygeno.2019.06.012
Ashrafizadeh, M., Paskeh, M. D. A., Mirzaei, S., Gholami, M. H., Zarrabi, A., Hashemi, F., et al. (2022). Targeting Autophagy in Prostate Cancer: Preclinical and Clinical Evidence for Therapeutic Response. J. Exp. Clin. Cancer Res. 41 (1), 105. doi:10.1186/s13046-022-02293-6
Ashrafizaveh, S., Ashrafizadeh, M., Zarrabi, A., Husmandi, K., Zabolian, A., Shahinozzaman, M., et al. (2021). Long Non-coding RNAs in the Doxorubicin Resistance of Cancer Cells. Cancer Lett. 508, 104–114. doi:10.1016/j.canlet.2021.03.018
Barresi, V., Bonetti, L. R., Leni, A., Caruso, R. A., and Tuccari, G. (2015) Histological Grading in Colorectal Cancer: New Insights and Perspectives Histology and Histopathology Cellular and Molecular Biology. Histol. Histopathol. 30(9):1059–1067.doi:10.14670/HH-11-633
Burada, F., Nicoli, E. R., Ciurea, M. E., Uscatu, D. C., Ioana, M., and Gheonea, D. I. (2015). Autophagy in Colorectal Cancer: An Important Switch from Physiology to Pathology. Wjgo 7 (11), 271–284. doi:10.4251/wjgo.v7.i11.271
Carr, P. R., Weigl, K., Jansen, L., Walter, V., Erben, V., Chang-Claude, J., et al. (2018). Healthy Lifestyle Factors Associated with Lower Risk of Colorectal Cancer Irrespective of Genetic Risk. Gastroenterology 155 (6), 1805–1815.e5. doi:10.1053/j.gastro.2018.08.044
Chen, F. J., WuChen, J. F., Xia, R., Yang, J., Huo, X. Y., Gu, D. Y., et al. (2018). STAT3-induced lncRNA HAGLROS Overexpression Contributes to the Malignant Progression of Gastric Cancer Cells via mTOR Signal-Mediated Inhibition of Autophagy. Mol. Cancer 17 (1), 6–16. doi:10.1186/s12943-017-0756-y
Du, T., Gao, Q., Zhao, Y., Gao, J., Li, J., Wang, L., et al. (2021). Long Non-coding RNA LINC02474 Affects Metastasis and Apoptosis of Colorectal Cancer by Inhibiting the Expression of GZMB. Front. Oncol. 11, 651796. doi:10.3389/fonc.2021.651796
Fan, S., Tian, T., Lv, X., Lei, X., Yang, Z., Liu, M., et al. (2020). lncRNA CISAL Inhibits BRCA1 Transcription by Forming a Tertiary Structure at its Promoter. iScience 23 (2), 100835. doi:10.1016/j.isci.2020.100835
Fang, C., Dai, L., Wang, C., Fan, C., Yu, Y., Yang, L., et al. (2021). Secretogranin II Impairs Tumor Growth and Angiogenesis by Promoting Degradation of Hypoxia‐inducible Factor‐1α in Colorectal Cancer. Mol. Oncol. 15 (12), 3513–3526. doi:10.1002/1878-0261.13044
Francesca, G., Sadia, A., Forbes-Hernandez, T. Y., Massimiliano, G., Danila, C., Patricia, R. R., et al. (20182017). Autophagy in Human Health and Disease: Novel Therapeutic Opportunities. Antioxidants Redox Signal. 30, 7234.
Frankel, L. B., Lubas, M., and Lund, A. H. (2017). Emerging Connections between RNA and Autophagy. Autophagy 13 (1), 3–23. doi:10.1080/15548627.2016.1222992
Ghafouri-Fard, S., Asadi, M., Sohrabi, B., Arsang-Jang, S., Mehravaran, E., Taheri, M., et al. (2021). Down-regulation of a Panel of Immune-Related lncRNAs in Breast Cancer. Pathology - Res. Pract. 224, 153534. doi:10.1016/j.prp.2021.153534
Jin, P., Jiang, J., Xie, N., Zhou, L., Huang, Z., Zhang, L., et al. (2019). MCT1 Relieves Osimertinib-Induced CRC Suppression by Promoting Autophagy through the LKB1/AMPK Signaling. Cell Death Dis. 10 (8), 615. doi:10.1038/s41419-019-1844-2
Kimmelman, A. C., and White, E. (2017). Autophagy and Tumor Metabolism. Cell Metab. 25 (5), 1037–1043. doi:10.1016/j.cmet.2017.04.004
Li, V. S. W., Yuen, S. T., Chan, T. L., Yan, H. H. N., Law, W. L., Yeung, B. H. Y., et al. (2009). Frequent Inactivation of Axon Guidance Molecule RGMA in Human Colon Cancer through Genetic and Epigenetic Mechanisms. Gastroenterology 137 (1), 176–187. doi:10.1053/j.gastro.2009.03.005
Li, Z.-L., Zhang, H.-L., Huang, Y., Huang, J.-H., Sun, P., Zhou, N.-N., et al. (2020). Autophagy Deficiency Promotes Triple-Negative Breast Cancer Resistance to T Cell-Mediated Cytotoxicity by Blocking Tenascin-C Degradation. Nat. Commun. 11 (1), 3806. doi:10.1038/s41467-020-17395-y
Lin, T., Cheng, H., Liu, D., Wen, L., Kang, J., Xu, L., et al. (2021). A Novel Six Autophagy-Related Genes Signature Associated with Outcomes and Immune Microenvironment in Lower-Grade Glioma. Front. Genet. 12, 698284. doi:10.3389/fgene.2021.698284
Liu, M., Han, T., Shi, S., and Chen, E. (2018). Long Noncoding RNA HAGLROS Regulates Cell Apoptosis and Autophagy in Lipopolysaccharides-Induced WI-38 Cells via Modulating miR-100/NF-κB axis. Biochem. biophysical Res. Commun. 500 (3), 589–596. doi:10.1016/j.bbrc.2018.04.109
Liu, M., Sun, T., Li, N., Peng, J., Fu, D., Li, W., et al. (2019). BRG1 Attenuates Colonic Inflammation and Tumorigenesis through Autophagy-dependent Oxidative Stress Sequestration. Nat. Commun. 10 (1), 4614. doi:10.1038/s41467-019-12573-z
Liu, X., Ma, B., Chen, M., Zhang, Y., Ma, Z., and Chen, H. (2021). Prognostic Autophagy-Related Genes of Gastric Cancer Patients on Chemotherapy. Front. Genet. 12, 720849. doi:10.3389/fgene.2021.720849
Luo, Y., Zheng, S., Wu, Q., Wu, J., Zhou, R., Wang, C., et al. (2021). Long Noncoding RNA (lncRNA) EIF3J-DT Induces Chemoresistance of Gastric Cancer via Autophagy Activation. Autophagy 17 (12), 4083–4101. doi:10.1080/15548627.2021.1901204
Lv, Y., Wang, Y., Zhang, Z., Bao, J., and Su, H. (2022). Potentials of Long Non-coding RNAs as Biomarkers of Colorectal Cancer. Clin. Transl. Oncol. doi:10.1007/s12094-022-02834-7Epub 2022/05/18Available: https://www.ncbi.nlm.nih.gov/pubmed/35581419
Mi Je, E., Gwak, M., Oh, H., Ryoung Choi, M., Jin Choi, Y., Lee, S. H., et al. (2013). Frameshift Mutations of Axon Guidance Genes ROBO1 and ROBO2 in Gastric and Colorectal Cancers with Microsatellite Instability. Pathology 45 (7), 645–650. doi:10.1097/pat.0000000000000007
Miller, K. D., Nogueira, L., Mariotto, A. B., Rowland, J. H., Yabroff, K. R., Alfano, C. M., et al. (20192019). Cancer Treatment and Survivorship Statistics, 2019. CA Cancer J. Clin. 69 (5), 363–385. doi:10.3322/caac.21565
Mirzaei, S., Gholami, M. H., Hushmandi, K., Hashemi, F., Zabolian, A., Canadas, I., et al. (2022). The Long and Short Non-coding RNAs Modulating EZH2 Signaling in Cancer. J. Hematol. Oncol. 15 (1), 18. doi:10.1186/s13045-022-01235-1
Mirzaei, S., Zarrabi, A., Hashemi, F., Zabolian, A., Saleki, H., Ranjbar, A., et al. (2021). Regulation of Nuclear Factor-KappaB (NF-κB) Signaling Pathway by Non-coding RNAs in Cancer: Inhibiting or Promoting Carcinogenesis? Cancer Lett. 509, 63–80. doi:10.1016/j.canlet.2021.03.025
Nana, K., and Edward, G. (2019). Global Burden of Colorectal Cancer: Emerging Trends, Risk Factors and Prevention Strategies. Nat. Rev. Gastroenterology Hepatology 16 (12), 713–732. doi:10.1038/s41575-019-0189-8
Øines, M., Helsingen, L. M., Bretthauer, M., and Emilsson, L. (2017). Epidemiology and Risk Factors of Colorectal Polyps. Best Pract. Res. Clin. Gastroenterology 31 (4), 419–424. doi:10.1016/j.bpg.2017.06.004
Siegel, R. L., Miller, K. D., and Jemal, A. (2020). Cancer Statistics, 2020. CA A Cancer J. Clin. 70 (1), 7–30. doi:10.3322/caac.21590
Sun, Y., Wang, X., and Bu, X. (2021). LINC02381 Contributes to Cell Proliferation and Hinders Cell Apoptosis in Glioma by Transcriptionally Enhancing CBX5. Brain Res. Bull. 176, 121–129. doi:10.1016/j.brainresbull.2021.07.009
Thongchot, S., Vidoni, C., Ferraresi, A., Loilome, W., Yongvanit, P., Namwat, N., et al. (2018). Dihydroartemisinin Induces Apoptosis and Autophagy-dependent Cell Death in Cholangiocarcinoma through a DAPK1-BECLIN1 Pathway. Mol. Carcinog. 57 (12), 1735–1750. doi:10.1002/mc.22893
Tokunaga, R., Nakagawa, S., Sakamoto, Y., Nakamura, K., Naseem, M., Izumi, D., et al. (2020). 12‐Chemokine Signature, a Predictor of Tumor Recurrence in Colorectal Cancer. Int. J. Cancer 147 (2), 532–541. doi:10.1002/ijc.32982
Ulitsky, I., and Bartel, D. P. (2013). lincRNAs: Genomics, Evolution, and Mechanisms. Cell 154 (1), 26–46. doi:10.1016/j.cell.2013.06.020
Wang, J., Yu, S., Chen, G., Kang, M., Jin, X., Huang, Y., et al. (2020). A Novel Prognostic Signature of Immune-Related Genes for Patients with Colorectal Cancer. J. Cell Mol. Med. 24 (1), 8491–8504. doi:10.1111/jcmm.15443
Wei, J., Ge, X., Tang, Y., Qian, Y., Lu, W., Jiang, K., et al. (2020). An Autophagy-Related Long Noncoding RNA Signature Contributes to Poor Prognosis in Colorectal Cancer. J. Oncol. 2020, 4728947. doi:10.1155/2020/4728947
Weiser, M. R. (2018). AJCC 8th Edition: Colorectal Cancer. Ann. Surg. Oncol. 25 (6), 1454–1455. doi:10.1245/s10434-018-6462-1
Wiener, Z., Högström, J., Hyvönen, V., Band, A. M., Kallio, P., Holopainen, T., et al. (2014). Prox1 Promotes Expansion of the Colorectal Cancer Stem Cell Population to Fuel Tumor Growth and Ischemia Resistance. Cell Rep. 8 (6), 1943–1956. doi:10.1016/j.celrep.2014.08.034
Wu, D., Yin, Z., Ji, Y., Li, L., Li, Y., Meng, F., et al. (2021). Identification of Novel Autophagy-Related lncRNAs Associated with a Poor Prognosis of Colon Adenocarcinoma through Bioinformatics Analysis. Sci. Rep. 11 (1), 8069. doi:10.1038/s41598-021-87540-0
Wu, Q., Zhao, Y., Shi, R., and Wang, T. (2021). LncRNA SNHG16 Facilitates Nasopharyngeal Carcinoma Progression by Acting as ceRNA to Sponge miR-520a-3p and Upregulate MAPK1 Expression. Cmar 13, 4103–4114. doi:10.2147/cmar.S305544
Yang, Y., Feng, M., Bai, L., Zhang, M., Zhou, K., Liao, W., et al. (2021). The Effects of Autophagy-Related Genes and lncRNAs in Therapy and Prognosis of Colorectal Cancer. Front. Oncol. 11, 582040. doi:10.3389/fonc.2021.582040
Zhang, L. H., Li, L. H., Zhang, P. F., Cai, Y. F., and Hua, D. (2007). LINC00957 Acted as Prognostic Marker Was Associated with Fluorouracil Resistance in Human Colorectal Cancer. Front. Oncol. 9, 776. doi:10.3389/fonc.2019.00776
Zhou, J., Dong, Y., Cai, X., Yang, H., and Guo, T. (2021). Identification and Validation of Autophagy-Related Genes as Potential Biomarkers and Therapeutic Targets in Atrial Fibrillation. Ijgm 14, 7783–7796. doi:10.2147/ijgm.S337855
Zhou, W., Zhang, S., Li, H.-b., Cai, Z., Tang, S., Chen, L.-x., et al. (2020). Development of Prognostic Indicator Based on Autophagy-Related lncRNA Analysis in Colon Adenocarcinoma. BioMed Res. Int. 2020, 9807918. doi:10.1155/2020/9807918Epub 2020/09/24Available: https://www.ncbi.nlm.nih.gov/pubmed/32964050
Zhu, G., Pei, L., Yang, F., and Zhang, C. (2022). A Robust Immune-Related lncRNA Signature for the Prognosis of Human Colorectal Cancer. Biosci. Rep. BSR20220078. doi:10.1042/BSR20220078Epub 2022/05/05Available: https://www.ncbi.nlm.nih.gov/pubmed/35506372
Keywords: colorectal adenocarcinoma, long non-coding RNA, autophagy, overall survival, prediction model
Citation: Duan L, Xia Y, Li C, Lan N and Hou X (2022) Identification of Autophagy-Related LncRNA to Predict the Prognosis of Colorectal Cancer. Front. Genet. 13:906900. doi: 10.3389/fgene.2022.906900
Received: 29 March 2022; Accepted: 17 June 2022;
Published: 11 August 2022.
Edited by:
Ángeles Carlos Reyes, Instituto Nacional de Enfermedades Respiratorias (INER), MexicoReviewed by:
Milad Ashrafizadeh, Sabancı University, TurkeyLeila Jahangiri, University of Cambridge, United Kingdom
Srikanth Talluri, Dana–Farber Cancer Institute, United States
Copyright © 2022 Duan, Xia, Li, Lan and Hou. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xiaoming Hou, ldyy_duanl@lzu.edu.cn