Disease Markers

Disease Markers / 2021 / Article
Special Issue

Biomarkers for Urological Malignancies

View this Special Issue

Research Article | Open Access

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

Yundong Xuan, Weihao Chen, Kan Liu, Yu Gao, Shidong Zuo, Baojun Wang, Xin Ma, Xu Zhang, "A Risk Signature with Autophagy-Related Long Noncoding RNAs for Predicting the Prognosis of Clear Cell Renal Cell Carcinoma: Based on the TCGA Database and Bioinformatics", Disease Markers, vol. 2021, Article ID 8849977, 18 pages, 2021. https://doi.org/10.1155/2021/8849977

A Risk Signature with Autophagy-Related Long Noncoding RNAs for Predicting the Prognosis of Clear Cell Renal Cell Carcinoma: Based on the TCGA Database and Bioinformatics

Academic Editor: Mingjun Shi
Received26 Aug 2020
Revised02 Apr 2021
Accepted23 Apr 2021
Published07 May 2021

Abstract

Background. Disorders of autophagic processes have been reported to affect the survival outcome of clear cell renal cell carcinoma (ccRCC) patients. The purpose of our study was to identify and validate the candidate prognostic long noncoding RNA signature of autophagy. Methods. Transcriptome profiles were obtained from The Cancer Genome Atlas. The autophagy gene list was obtained from the Human Autophagy Database. Based on coexpression analysis, we obtained a list of autophagy-related lncRNAs (ARlncRNAs). GO enrichment analysis and KEGG pathway analysis were conducted to explore the functional annotation of these ARlncRNAs. Univariate and multivariate Cox regression analyses were conducted to elucidate the correlation between overall survival and the expression level of each ARlncRNAs. We then established a prognostic signature that was a linear combination of the regression coefficients from the multivariate Cox regression model () multiplied by the expression levels of the respective ARlncRNAs in the training cohort. The predictive performance was tested in the validation cohort. Additionally, the independence of the risk signature was assessed, and the relationship between the risk signature and conventional clinicopathological features was explored. Results. Seven autophagy-related lncRNAs with prognostic value (SNHG3, SNHG17, MELTF-AS1, HOTAIRM1, EPB41L4A-DT, AP003352.1, and AC145423.2) were identified and integrated into a risk signature, dividing patients into low-risk and high-risk groups. The risk signature was independent of conventional clinical characteristics as a prognostic indicator of ccRCC (HR, 1.074, 95% confidence interval: 1.036-1.113, ) and was valuable in the prediction of ccRCC progression. Conclusion. Our risk signature has potential prognostic value in ccRCC, and these ARlncRNAs may play a significant role in ccRCC tumor biology.

1. Introduction

Renal cell carcinoma (RCC), a principal malignancy of the renal tubular epithelium, ranks third among urinary cancers1 [1]. Characterized by multifarious genetic features [2], RCC is composed of different histopathologic subtypes, of which clear cell renal cell carcinoma (ccRCC) represents the principle pathologic subtype, accounting for 70% to 80% of RCCs [3]. The prognosis of ccRCC varies between patients with complicated genetic mutations, such as von Hippel-Lindau (VHL), PBRM1, and BAP1 [4]. Surgery remains the principal treatment for ccRCC due to its therapeutic value, but 40% of ccRCC patients will eventually suffer distant metastasis [5], and few exhibit a positive response to radiotherapy or chemotherapy [6]. Neither immunologic nor targeted therapy can definitely benefit ccRCC patients [7]. Therefore, predicting the progression and prognosis of ccRCC by seeking novel effective biomarkers might facilitate the therapeutic schedule and the evaluation of survival status.

Autophagy, a process in which cells engulf proteins and organelles by forming double-membraned autophagic vesicles where degradation occurs, is regarded as a recycling of organelles and an adaptation of metabolism. Autophagy plays a tumor-suppressor role by inhibiting the selection and expansion of tumor-initiating cells early in tumor development [8]. However, increasing evidence in established tumors suggests that autophagy can help cope with environmental or intracellular stresses, such as hypoxia, nutrient shortage, or cancer therapy, thereby promoting tumor growth [911]. Although an increasing number of studies have sought to identify novel potential targets by probing the autophagy pathway, and autophagic drugs have been reported to induce renal cancer cell death, the mechanism is still to be clarified [12].

Long noncoding RNAs (lncRNAs), a variety of noncoding RNA, participate in many cellular processes with multitudinous functions by modulating gene expression at the epigenetic, transcriptional, and posttranscriptional levels [1315]. Accumulating evidence suggests that lncRNAs can target autophagy-related genes at both the transcriptional and posttranscriptional levels to regulate the autophagy pathway, and lncRNAs regulate various proteins that function in the autophagy process [16, 17]. Dysregulated lncRNAs are involved in ccRCC progression and dissemination [18]. It has been suggested that autophagy-related lncRNAs can exert their action in tumor regulation [19].

The potential value of autophagy-related lncRNAs in evaluating the prognosis of ccRCC patients and their role as potential therapeutic targets have yet to be fully explored. Here, we sought to identify an autophagy-related lncRNA signature in ccRCC and to advance more personalized treatment guidelines for ccRCC through bioinformatic analysis.

2. Materials and Methods

2.1. Data Acquisition Processing

In our study, RNA-seq transcriptome data of which the workflow type was HTSeq-FPKM and clinical information data were downloaded from the TCGA database (https://portal.gdc.cancer.gov/) [20]. The lncRNA matrix and mRNA matrix were extracted from transcriptome profiling, respectively, by gene annotation. A list consisting of 232 autophagy genes was obtained from the Human Autophagy Database (HADb, http://autophagy.lu/clustering/index.html) (Table S1). Then, the expression matrix of 210 autophagy genes was extracted combing the mRNA matrix with autophagy genes list (Table S2), and 21 genes were abandoned due to lacking expressing information. In addition, the clinicopathological features included survival status, survival time, age, sex, International Society of Urological Pathology (ISUP) grade, and American Joint Committee on Cancer (AJCC) stage.

2.2. Autophagy-Related lncRNAs Screening

The Pearson correlation test was performed to screen autophagy-related lncRNAs based on 210 autophagy genes and all lncRNAs as mentioned. A lncRNA with a correlation coefficient and was considered an autophagy-related lncRNA. The “limma” package in R software was then used to screen [21] differentially expressed autophagy-related lncRNAs (ARlncRNAs) (Table S3).

2.3. Construction of a Coexpression Network of the Autophagy-Related lncRNAs and Building a Sankey Diagram

Firstly, the clinicopathological data and ARlncRNAs were merged. Then, we constructed a mRNA-lncRNA interaction network of prognostic ARlncRNAs by applying Cytoscape software 3.8.0 [22]. A Sankey diagram was built to describe the relationship between the autophagy genes and lncRNAs by using the R software packages “ggalluvial” [23] and “ggplot2” [24].

2.4. Enrichment Analysis via DAVID Bioinformatics Resources

In this study, we conducted enrichment analyses via DAVID web server, a free server resource (https://david.ncifcrf.gov/summary.jsp) to elucidate the biological functions of those lncRNAs. 17 differentially expressed genes and 19 ARlncRNAs (Table S4) were taken to perform the analyses. Gene Ontology annotation and KEGG pathway were both involved.

2.5. Establishment of the ARlncRNA Signature and Performance Evaluation

First, the entire cohort obtained from the TCGA dataset was randomly divided into two cohorts: a training cohort (266 patients) and a validation cohort (264 patients). The training cohort was used to establish the Cox regression risk signature, and the validation cohort was then used to assess the performance of the signature. Then, univariate Cox regression and Kaplan-Meier (KM) analysis were used to identify ARlncRNAs with prognostic value. Seven lncRNAs considered significant with a value < 0.0001 by both analyses were included in the resulting model. Next, multivariate stepwise regression Cox analysis was performed to establish a prognostic risk model. The risk score was calculated as the sum of multivariate Cox regression coefficient (β)-weighted expression levels of lncRNAs: . The training cohort was then separated into high-risk and low-risk groups based on the median risk score. We generated a receiver operating characteristic (ROC) curve to evaluate the predictive power of the risk score. The same analysis was adopted in the testing cohort to assess the performance of the signature. Additionally, we built a nomogram using the R package “rms” based on the results of multivariate Cox regression analysis to evaluate the prognosis of ccRCC. Lastly, the C-index and calibration curve were used to assess the performance value.

2.6. Independence Verification of the Risk Signature and Performance Evaluation

Univariate and multivariate Cox regression analyses were performed to investigate the independence of the risk signature as a predictive factor from the traditional clinical features (including age, sex, ISUP grade, and AJCC stage) in both the training and validation cohorts.

2.7. Statistical Analysis

In this study, all analyses were performed using R 4.0.2. A value < 0.05 was considered statistically significant. The correlation matrix was constructed by R software based on Pearson correlation coefficients. The relationship between autophagy-related lncRNAs and overall survival was analyzed through the Kaplan-Meier curve, which was evaluated by the log-rank test. Time-dependent ROC curves were used to analyze the sensitivity and specificity of the prognostic prediction model. The nomogram was constructed with the regression coefficients based on the Cox analysis, and its performance was assessed by the c-index and calibration curve. Univariate and multivariate Cox regression analyses were performed to investigate the independence of the risk signature as a predictive factor, and Student’s -test was used to compare the clinicopathological features between different risk-score groups.

3. Results

3.1. Data Processing and Prognostic ARlncR Acquisition

We obtained transcriptome data and clinical data from the TCGA dataset. The basic clinical characteristics of the ccRCC patients in the TCGA database are shown in Table 1. Based on a total of 232 autophagy-related genes from the Human Autophagy Database (HADb, http://autophagy.lu/clustering/index.html), we screened 813 autophagy-related lncRNAs. Nineteen ARlncRNAs closely connected to prognosis were identified in the training cohort through univariate Cox regression analysis ( value < 0.01), and then, a network of prognostic 19 lncRNAs with coexpressed 17 autophagy genes in ccRCC and a Sankey diagram for visualization were built. The results are shown in Figures 1 and 2.


CharacteristicsTotal%

Age at diagnosis (y)58 (26~90)
GenderFemale19135.57
Male34664.43
StageI26950.37
II5710.67
III12523.41
IV8315.55
GradeG1142.65
G223043.48
G320739.13
G47814.74
T stageT127551.21
T26912.85
T318233.89
T4112.05
N stageN024093.39
N1176.61
M stageM042684.36
M17915.64

3.2. Enrichment Analyses

To further identify the Gene Ontology (GO) annotation and KEGG pathway in which the above lncRNAs were enriched, GO term and KEGG pathway enrichment analyses were performed via DAVID functional annotation tool. The visualization of results was achieved by “ggplot2” packages in R software. GO analysis showed that changes in the biological processes (BPs) of autophagy-related lncRNAs were enriched in mitophagy, peptidyl-threonine phosphorylation, positive regulation of translation, macroautophagy, and others (Table 2, Figure 3(a)). Changes in cell components (CCs) included the autophagosome, late endosome, and mitochondrion (Table 2, Figure 3(a)). Changes in molecular functions (MFs) were mainly enriched in kinase activity, protein kinase activity, protein serine or threonine activity, and others (Table 2, Figure 3(a)). KEGG pathway enrichment analysis indicated that those autophagy-related lncRNAs were involved in multiple tumor progressions and signaling pathways such as NOD-like receptor signaling pathway and mTOR signaling pathway (Table S5, Figure 3(b)).


CategoryIDDescription valueFDRGenes

Biological processGO:0000422Mitophagy4.97-040.038125PINK1, CDKN2A, ATG4B
Biological processGO:0018017Peptidyl-threonine phosphorylation5.87-040.038125MAPK8, BCL2, MAPK1
Biological processGO:0045727Positive regulation of translation1.14-030.063524PINK1, ERBB2, MAPK1
Biological processGO:0016236Macroautophagy4.78-050.006217PINK1, NBR1, ATG4B, PIK3C3
Biological processGO:0006950Response to stress1.50-030.073488PINK1, MAPK8, MAPK1
Biological processGO:0031647Regulation of protein stability1.98-030.002612CDKN2A, BCL2, MAPK1
Biological processGO:0006914Autophagy6.36-060.001240PINK1, NBR1, ATG4B, PIK3C3
Biological processGO:0018105Peptidyl-serine phosphorylation2.10-040.020483PINK1, MAPK8, BCL2, MAPK1
Biological processGO:0006468Protein phosphorylation2.47-069.62E-04IKBKB, PINK1, MAPK8, ERBB2, MAPK1, BICR5, PIK3C3
Cellular componentGO:0005776Autophagosome1.4-030.029687NBR1, RAB24, WDFY3
Cellular componentGO:0005770Late endosome5.0-030.069394NBR1, MAPK1, PIK3C3
Cellular componentGO:0005739Mitochondrion6.94-050.002878PINK1, MAPK8, CDKN2A, CAPN10, RAB24, BCL2, MAPK1
Cellular componentGO:0005622Intracellular6.37-040.017625IKBKB, MAPK8, HGS, CAPN10, RAB24, MAPK1, PIK3C3
Cellular componentGO:0005737Cytoplasm5.0-030.069395IKBKB, PINK1, CDKN2A, HGS, ERBB2, BCL2, ULK3, ATG4B, MAPK1, BIRC5, WDFY3
Molecular functionGO:0016301Kinase activity3.58-082.68E-06PINK1, MAPK8, CDKN2A, HGS, ERBB2, MAPK1, PIK3C3,
Molecular functionGO:0004672Protein kinase activity2.28-040.006804IKBKB, PINK1, ERBB2, ULK3, PIK3C3
Molecular functionGO:0004674Protein serine/threonine kinase activity0.00406432.72-040.020892
0.006803
BIRC5, PINK1
IKBKB, PINK1, MAPK8, ULK3, MAPK1
Molecular functionGO:0046982Protein heterodimerization activity0.0073820.092273IKBKB, ERBB2, BCL2, BIRC5
Molecular functionGO:0005524ATP binding1.2-030.017773IKBKB, PINK1, MAPK8, ERBB2, ULK3, MAPK1, PIK3C3
Molecular functionGO:0004861Protein binding8.21-040.015391CDKN2A, IKBKB, PINK1, MAPK8, HGS, NBR1, ERBB2, RAB24, ULK3, BCL2, BIRC5, ATG4B, MAPK1, PIK3C3, WDFY3

3.3. Establishment of the Autophagy-Related lncRNA Signature for ccRCC

To improve our prognostic ability further, we employed stepwise multivariate Cox regression analysis to further evaluate the aforementioned ARlncRNAs. As a result, 7 lncRNAs were eventually pulled out to construct the signature. Herein, the risk score was assigned using a linear combination of the expression levels of the 7 identified lncRNAs weighted by their regression coefficients () (Table 3, Figure 4). Then, the ccRCC patients were divided into high-risk and low-risk groups around the median risk score. As a result, the risk score distribution of the patients on the basis of the prognostic signature is shown in Figure 5(a). Survival status scatter plots for the patients based on the prognostic model are shown in Figure 5(b), indicating that patients in the high-risk group had a higher mortality than those in the low-risk group. A significant difference in overall survival between the two groups was seen (Figure 5(c), ), and the AUCs at one, three, and five years were 0.754, 0.791, and 0.808, respectively (Figures 5(d)).


IDCoefHRHR.95 LHR.95 H value

SNHG3-0.6477620.5232150.2819570.9709040.040017
SNHG170.9668782.6297241.4530594.7592320.001400
MELTF-AS10.5282761.6960061.0940152.6292470.018193
HOTAIRM10.4962321.6425201.0980482.4569720.015726
EPB41L4A-DT-1.6172470.1984440.1042380.3777898.51E-07
AP003352.1-0.9109120.4021570.1688550.9578020.039653
AC145423.20.5115531.6678800.9234323.0124840.089907

3.4. Evaluation of the Prediction Performance of the Signature

Next, we assessed the predictive ability of the prognostic signature in the validation cohort to further evaluate its performance. The risk scores of the above 7 lncRNAs were recalculated for each patient in the validation cohort. As shown in (Figure 6(a)), the risk score distribution of the patients based on the prognostic model presented similar results as in the training cohort. Survival status scatter plots for the patients on the basis prognostic model are presented in (Figure 6(b)). Survival differences were significant in the validation cohort (Figure 6(c)), ). The AUCs at one, three, and five years were 0.686, 0.673, and 0.711, respectively (Figures 6(d)). In addition, we constructed a nomogram to forecast the overall survival in the entire cohort with ccRCC according to the clinicopathological features and risk scores (Figure 7(a)), of which the C-index was 0.773. The calibration curve for the nomogram suggested good performance. Thus, the nomogram proved to be of value in the prediction of the prognosis of ccRCC patients (Figures 7(b) and 7(c)).

3.5. Independence Verification of the Signature as a Prognostic Predictor

To further investigate whether the risk score could be predictive for ccRCC independent of conventional clinicopathological features (age, sex, ISUP grade, and AJCC state), univariate and multivariate Cox regression analyses were run. The risk score proved to be independent of the aforementioned clinical features in the entire cohort (Figure 8) (HR, 1.074, 95% confidence interval: 1.036-1.113, ). The independence of the signature was further validated by stratified clinical features.

We conducted analyses to explore the value of the ARlncRNA signature in different clinicopathological subgroups, including age, sex, ISUP grade, AJCC stage, T stage, and M stage. Comparing the two results of each stratified feature revealed that the overall survival time of the high-risk group was shorter than that of the low-risk group in all clinicopathological subgroups (Figure 9). These results further indicated that the autophagy-related lncRNA signature can independently predict the prognosis of ccRCC patients.

3.6. Exploring the Predictive Value of the Signature for Tumor Progression

In order to get a deeper insight into the predictive value of the signature for ccRCC patients with regard to tumor progression, correlation analyses between the autophagy-related prognostic signature and clinicopathological features were performed. The risk score of stage III–IV ccRCC was higher than that of stage I–II (, Figure 10(a)). The risk score of G3–4 was higher than that of G1–2 (p = 4.247e−07, Figure 10(b)). The risk score of T3–4 was higher than that of T1–2 (, Figure 10(c)). The risk score of M1 was higher than that of M0 (, Figure 10(d)), and the risk score of N1 was higher than that of N0 (, Figure 10(e)). Taken together, these results indicated that as the risk score increased, the malignancy of ccRCC increased. Thus, the prognostic signature was of unique value for predicting the progression of ccRCC.

3.7. Relationships between the Prognostic ARlncRNAs and Clinicopathological Features

Subsequently, we investigated the connection between those ARlncRNAs and clinicopathological features. The purpose was to help develop a deeper understanding of the autophagy process, and the results shown in Table 4 revealed that each of those prognostic ARlncRNAs was significantly associated with clinicopathological features that are closely connected with tumor progression, including ISUP grade, AJCC stage, T stage, and N stage. SNHG17 and AP003352.1 were significantly associated with age. However, only EPB41L4A-DT was significantly associated with sex. Overall, these results indicated that the aforementioned ARlncRNAs can promote tumor progression.


lncRNAAgeSexGradeStageT stageM stage
≤65>65FemaleMaleG1-2G3-4I-IIIII-IVT1-2T3-4M0M1

16632322126828920030618341277
SNHG3 value1.810.12463.6604.5654.4052.791
value0.0710.9090.0003<0.0001<0.00010.0054
SNHG17 value2.241.173.1525.1494.9593.228
value0.02560.00170.0017<0.0001<0.00010.0013
MELTF-AS1 value1.4271.1354.3196.6896.1854.611
value0.15420.2571<0.0001<0.0001<0.0001<0.0001
HOTAIRM1 value1.0820.90083.4804.4074.3532.503
value0.280.36820.0005<0.0001<0.00010.0126
EPB41L4A-DT value0.69852.7015.9205.8267.7124.475
value0.48520.0072<0.0001<0.0001<0.0001<0.0001
AP003352.1 value2.2070.66702.6494.724.4242.814
value0.0280.50510.0083<0.0001<0.00010.0051
AC145423.2 value0.21750.21753.9213.7422.8362.302
value0.82790.82790.00010.00020.00480.0218

4. Discussion

ccRCC is well known for its heterogeneity, exhibiting molecular diversity, morphological variability, and metabolic reprogramming. Late diagnosis without early warning signs and limited response to chemotherapy or radiotherapy are the main culprits of poor prognosis [25]. Although novel biomarkers, especially autophagy-related genes and molecules, are emerging as predictive factors thanks to in-depth cancer genetics and molecular biology discoveries [26, 27], the value of autophagy-related lncRNAs as prognostic indicators has not been addressed. Unlike previous studies that focused on the role of autophagy-related genes in tumorigenesis and progression [2830], our study is aimed at improving prognostic prediction by finding autophagy-related lncRNAs associated with the poor prognosis of ccRCC through comprehensive bioinformatics analysis in TCGA databases.

We first identified 813 lncRNAs on the basis of the lncRNA-autophagy gene coexpression network. By using a univariate Cox regression model, we identified 19 ARlncRNAs associated with the prognosis of ccRCC patients. Seven ARlncRNAs were further screened using multivariate Cox regression analysis, including SNHG3, SNHG17, MELTF-AS1, HOTAIRM1, EPB41L4A-DT, AP003352.1, and AC145423.2. GO analysis was conducted to discover the main biological characteristics of these ARlncRNAs. Next, we constructed a risk score-based prognostic signature that separated ccRCC patients into low-risk and high-risk groups. The OS time in the high-risk group was shorter than that in the low-risk group. The prediction performance was validated in the validation cohorts.

Furthermore, through univariate and multivariate Cox regression analyses, the risk score based on the signature was shown to be a prognostic factor for ccRCC independent of conventional clinicopathological features (age, sex, ISUP grade, and AJCC stage). Further evaluation demonstrated that the ARlncRNA signature can independently predict the progression of ccRCC, which means that the higher the risk score was, the worse the prognosis and the greater the degree of malignancy were. Finally, we established a nomogram based on the risk score of the signature, and the C index and calibration curve indicated that the predictive performance of the nomogram was good. Overall, these results indicate that our ARlncRNA signature can play an important role in predicting the prognosis of ccRCC patients.

Autophagy is significantly connected to the prognosis of cancer; however, the complicated process and the numerous molecular interactions make autophagy play contradictory roles in cancer [31]. Long noncoding RNA small nucleolar RNA host gene 17 (SNHG17) was reported to be a critical regulator of tumorigenesis, and studies have reported its role in promoting tumor invasion and proliferation by activating the PI3K/AKT pathway [32, 33]. Jiang et al. [34] have reported that MELTF-AS1 can serve as a prognostic indicator and is associated with immunological processes. Our study further validates its prognostic value in association with the autophagy process, indicating that MELTF-AS can play distinct roles in multiple physio-pathological processes. HOTAIRM1 can promote glioblastoma progression [35], but in ccRCC, it is downregulated, serving as a suppressor of HIF1-dependent angiogenic pathways [36]. Thus, in vivo and in vitro experiments are required to illustrate the full effects of HOTAIRM1. Although we are the first to reveal the predictive value of EPB41L4A-DT, AP003352.1, and AC145423.2, their actual functions remain to be determined.

The common methods to detect lncRNA for clinical applications currently include microarray, lncRNA sequencing, and quantitative RT-PCR. Although the microarray can facilitate large-scale detections, its cost is relatively high [37]. LncRNA sequencing enables the detection to be more efficient and break the limit of traditional methods but leaving the library construction work tedious and expensive [38]. Quantitative RT-PCR is also a common method used to detect lncRNAs. It has been widely used because of its simplicity and low cost though it is still a low-throughput and low-specificity method [39]. Given those shortages, a comprehensive and personalized detection method of lncRNAs with further research achievements could be actually integrated into future precision medicine strategies.

Despite our novel findings, there are still limitations to this study. First, more basic experiments are required to clarify the mechanisms of action of ARlncRNAs in ccRCC tumor progression. Second, an external data set should be employed for validation, instead of our limited internal data set, to assess the consistency, reliability, and applicability of the autophagy-related signature.

5. Conclusion

In conclusion, we established an autophagy-related signature that independently predicted the prognosis of ccRCC patients. It can be used to guide individualized treatment regimens. Our study could help broaden the understanding of autophagy-related lncRNAs and narrow the gap between theoretical research and clinical practice, but the underlying mechanisms still urgently need to be understood to clarify the importance of our findings.

Data Availability

The data included in the current study were available in TCGA database (https://cancergenome.nih.gov/), the Human Autophagy Database (http://autophagy.lu/clustering/index.html), and DAVID (https://david.ncifcrf.gov/summary.jsp).

Ethical Approval

Because the dataset in the current study was downloaded from TCGA and data acquiring and application complied with TCGA publication guidelines and data access policies, additional approval by an ethics committee and consent to participate were not needed.

Conflicts of Interest

The authors declared that they have no conflict of interest.

Authors’ Contributions

Yundong Xuan and Shidong Zuo performed the data collection and the data analysis. Yundong Xuan and Kan Liu visualized and illustrated the results. Weihao Chen and Shidong Zuo conducted the statistic analyses. Yundong Xuan and Weihao Chen underwent the illustration of the results and carried out the initial manuscript writing. Baojun Wang and Xin Ma provided feedback on the manuscript and contributed to the study idea. All authors read and approved the final manuscript. Yundong Xuan and Weihao Chen contributed equally to this work.

Supplementary Materials

Table S1: 232 autophagy genes obtained from the Human Autophagy Database. Table S2: 210 autophagy genes expression matrix. Table S3: 813 autophagy-related lncRNAs. Table S4: 17 autophagy genes and 19 lncRNAs used to enrichment analyses. Table S5: results of enrichment via DAVID web serve. (Supplementary Materials)

References

  1. R. R. Kotecha, R. J. Motzer, and M. H. Voss, “Towards individualized therapy for metastatic renal cell carcinoma,” Nature Reviews. Clinical Oncology, vol. 16, no. 10, pp. 621–633, 2019. View at: Publisher Site | Google Scholar
  2. D. J. Clark, S. M. Dhanasekaran, F. Petralia et al., “Integrated proteogenomic characterization of clear cell renal cell carcinoma,” Cell, vol. 180, no. 1, p. 207, 2020. View at: Publisher Site | Google Scholar
  3. M. B. Atkins and N. M. Tannir, “Current and emerging therapies for first-line treatment of metastatic clear cell renal cell carcinoma,” Cancer Treatment Reviews, vol. 70, pp. 127–137, 2018. View at: Publisher Site | Google Scholar
  4. Y. Sato, T. Yoshizato, Y. Shiraishi et al., “Integrated molecular analysis of clear-cell renal cell carcinoma,” Nature Genetics, vol. 45, no. 8, pp. 860–867, 2013. View at: Publisher Site | Google Scholar
  5. A. Rao, C. Wiggins, and R. C. Lauer, “Survival outcomes for advanced kidney cancer patients in the era of targeted therapies,” Ann Transl Med, vol. 6, no. 9, p. 165, 2018. View at: Publisher Site | Google Scholar
  6. E. M. Posadas, S. Limvorasak, and R. A. Figlin, “Targeted therapies for renal cell carcinoma,” Nature Reviews. Nephrology, vol. 13, no. 8, pp. 496–511, 2017. View at: Publisher Site | Google Scholar
  7. A. K. Lalani, B. A. McGregor, L. Albiges et al., “Systemic treatment of metastatic clear cell renal cell carcinoma in 2018: current paradigms, use of immunotherapy, and future directions,” European Urology, vol. 75, no. 1, pp. 100–110, 2019. View at: Publisher Site | Google Scholar
  8. R. Amaravadi, A. C. Kimmelman, and E. White, “Recent insights into the function of autophagy in cancer,” Genes & Development, vol. 30, no. 17, pp. 1913–1930, 2016. View at: Publisher Site | Google Scholar
  9. G. Kroemer, G. Mariño, and B. Levine, “Autophagy and the integrated stress response,” Molecular Cell, vol. 40, no. 2, pp. 280–293, 2010. View at: Publisher Site | Google Scholar
  10. J. M. Farrow, J. C. Yang, and C. P. Evans, “Autophagy as a modulator and target in prostate cancer,” Nature Reviews. Urology, vol. 11, no. 9, pp. 508–516, 2014. View at: Publisher Site | Google Scholar
  11. E. White, J. M. Mehnert, and C. S. Chan, “Autophagy, metabolism, and cancer,” Clinical Cancer Research, vol. 21, no. 22, pp. 5037–5046, 2015. View at: Publisher Site | Google Scholar
  12. H. M. Kang, K. H. Noh, T. K. Chang et al., “Ubiquitination of MAP1LC3B by pVHL is associated with autophagy and cell death in renal cell carcinoma,” Cell Death & Disease, vol. 10, no. 4, p. 279, 2019. View at: Publisher Site | Google Scholar
  13. T. R. Mercer and J. S. Mattick, “Structure and function of long noncoding RNAs in epigenetic regulation,” Nature Structural & Molecular Biology, vol. 20, no. 3, pp. 300–307, 2013. View at: Publisher Site | Google Scholar
  14. Y. Long, X. Wang, D. T. Youmans, and T. R. Cech, “How do lncRNAs regulate transcription?” Science Advances, vol. 3, no. 9, article eaao2110, 2017. View at: Publisher Site | Google Scholar
  15. F. Kopp and J. T. Mendell, “Functional classification and experimental dissection of long noncoding RNAs,” Cell, vol. 172, no. 3, pp. 393–407, 2018. View at: Publisher Site | Google Scholar
  16. Z.-H. Chen, W.-T. Wang, W. Huang et al., “The lncRNA HOTAIRM1 regulates the degradation of PML-RARA oncoprotein and myeloid cell differentiation by enhancing the autophagy pathway,” Cell Death and Differentiation, vol. 24, no. 2, pp. 212–224, 2017. View at: Publisher Site | Google Scholar
  17. L.-L. Wang, L. Zhang, and X.-F. Cui, “Downregulation of long noncoding RNA LINC01419 inhibits cell migration, invasion, and tumor growth and promotes autophagyviainactivation of the PI3K/Akt1/mTOR pathway in gastric cancer,” Therapeutic advances in medical oncology, vol. 11, p. 175883591987465, 2019. View at: Publisher Site | Google Scholar
  18. 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
  19. L. B. Frankel, M. Lubas, and A. H. Lund, “Emerging connections between RNA and autophagy,” Autophagy, vol. 13, no. 1, pp. 3–23, 2017. View at: Publisher Site | Google Scholar
  20. Z. Wang, M. A. Jensen, and J. C. Zenklusen, “A practical guide to the cancer genome atlas (TCGA),” Methods in Molecular Biology, vol. 1418, pp. 111–141, 2016. View at: Publisher Site | Google Scholar
  21. M. E. Ritchie, B. Phipson, D. I. 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
  22. P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003. View at: Publisher Site | Google Scholar
  23. J. Brunson, “Ggalluvial: layered grammar for alluvial plots,” Journal of Open Source Software, vol. 5, no. 49, 2017. View at: Publisher Site | Google Scholar
  24. H. Wickham, ggplot2, Springer International Publishing, Cham, 2016.
  25. Cancer Genome Atlas Research Network, “Comprehensive molecular characterization of clear cell renal cell carcinoma,” Nature, vol. 499, no. 7456, pp. 43–49, 2013. View at: Publisher Site | Google Scholar
  26. A. C. Kimmelman and E. White, “Autophagy and tumor metabolism,” Cell Metabolism, vol. 25, no. 5, pp. 1037–1043, 2017. View at: Publisher Site | Google Scholar
  27. B. Levine and G. Kroemer, “Biological functions of autophagy genes: a disease perspective,” Cell, vol. 176, no. 1–2, pp. 11–42, 2019. View at: Publisher Site | Google Scholar
  28. P. Chen, M. Cescon, and P. Bonaldo, “Autophagy-mediated regulation of macrophages and its applications for cancer,” Autophagy, vol. 10, no. 2, pp. 192–200, 2014. View at: Publisher Site | Google Scholar
  29. G.-M. Jiang, Y. Tan, H. Wang et al., “The relationship between autophagy and the immune system and its applications for tumor immunotherapy,” Molecular Cancer, vol. 18, no. 1, p. 17, 2019. View at: Publisher Site | Google Scholar
  30. A. Yang, N. V. Rajeshkumar, X. Wang et al., “Autophagy is critical for pancreatic tumor growth and progression in tumors with p53 alterations,” Cancer Discovery, vol. 4, no. 8, pp. 905–913, 2014. View at: Publisher Site | Google Scholar
  31. B. Levine and G. Kroemer, “Autophagy in the pathogenesis of disease,” Cell, vol. 132, no. 1, pp. 27–42, 2008. View at: Publisher Site | Google Scholar
  32. T. Xu, S. Yan, L. Jiang et al., “Gene amplification-driven long noncoding RNA SNHG17 regulates cell proliferation and migration in human non-small-cell lung cancer,” Mol Ther Nucleic Acids, vol. 17, pp. 405–413, 2019. View at: Publisher Site | Google Scholar
  33. H. Gao, R. Liu, and X. Sun, “STAT3-induced upregulation of lncRNA SNHG17 predicts a poor prognosis of melanoma and promotes cell proliferation and metastasis through regulating PI3K-AKT pathway,” European Review for Medical and Pharmacological Sciences, vol. 23, no. 18, pp. 8000–8010, 2019. View at: Publisher Site | Google Scholar
  34. Y. Jiang, X. Gou, Z. Wei et al., “Bioinformatics profiling integrating a three immune-related long non-coding RNA signature as a prognostic model for clear cell renal cell carcinoma,” Cancer Cell International, vol. 20, no. 1, 2020. View at: Publisher Site | Google Scholar
  35. T. R. Spitzer, H. M. Lazarus, R. J. Creger, and N. A. Berger, “High-dose melphalan, misonidazole, and autologous bone marrow transplantation for the treatment of metastatic colorectal carcinoma. A phase I study,” American Journal of Clinical Oncology, vol. 12, no. 2, pp. 145–151, 1989. View at: Publisher Site | Google Scholar
  36. M. J. Hamilton, M. Young, K. Jang et al., “_HOTAIRM1_ lncRNA is downregulated in clear cell renal cell carcinoma and inhibits the hypoxia pathway,” Cancer Letters, vol. 472, pp. 50–58, 2020. View at: Publisher Site | Google Scholar
  37. C. Zhang, W. Gao, S. Wen et al., “Potential key molecular correlations in laryngeal squamous cell carcinoma revealed by integrated analysis of mRNA, miRNA and lncRNA microarray profiles,” Neoplasma, vol. 63, no. 6, pp. 888–900, 2016. View at: Publisher Site | Google Scholar
  38. X. Qian, J. Zhao, P. Y. Yeung, Q. C. Zhang, and C. K. Kwok, “Revealing lncRNA structures and interactions by sequencing-based approaches,” Trends in Biochemical Sciences, vol. 44, no. 1, pp. 33–52, 2019. View at: Publisher Site | Google Scholar
  39. M. Kazimierczyk, M. K. Kasprowicz, M. E. Kasprzyk, and J. Wrzesinski, “Human long noncoding RNA interactome: detection, characterization and function,” International journal of molecular sciences, vol. 21, no. 3, p. 1027, 2020. View at: Publisher Site | Google Scholar

Copyright © 2021 Yundong Xuan 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.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views272
Downloads403
Citations

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2021, as selected by our Chief Editors. Read the winning articles.