Construction of Competitive Endogenous RNA Network and Verification of 3-Key LncRNA Signature Associated With Distant Metastasis and Poor Prognosis in Patients With Clear Cell Renal Cell Carcinoma

Established and emerging biomarkers of immunotherapy in renal cell carcinoma

Immunotherapies, such as immune checkpoint inhibitors, have heralded impressive progress for patient care in renal cell carcinoma (RCC). Despite this success, some patients' disease fails to respond, and other patients experience significant side effects. Thus, development of biomarkers is needed to ensure that patients can be selected to maximize benefit from immunotherapies. Improving clinicians' ability to predict which patients will respond to immunotherapy and which are most at risk of adverse events - namely through clinical biomarkers - is indispensable for patient safety and therapeutic efficacy. Accordingly, an evolving suite of therapeutic biomarkers continues to be investigated. This review discusses biomarkers for immunotherapy in RCC, highlighting current practices and emerging innovations, aiming to contribute to improved outcomes for patients with RCC.

Identification of AP002498.1 and LINC01871 as prognostic biomarkers and therapeutic targets for distant metastasis of colorectal adenocarcinoma

BACKGROUND

Increasing evidence suggests that lncRNA (Long non-coding RNA, lncRNA)-mediated ceRNA (competing endogenous RNA, ceRNA) networks are involved in the occurrence and progression of colorectal cancer (CRC). However, the roles of the lncRNA-miRNA-mRNA ceRNA network in distant metastasis of CRC are still unclear.

METHODS

In this study, we constructed a specific ceRNA network to identify potential biomarkers and therapeutic targets for distant metastasis of CRC. Specifically, RNA-Seq data from The Cancer Genome Atlas (TCGA) were used to screen for differentially expressed lncRNAs (DElncRNAs) and mRNAs (DEmRNAs) related to metastasis. After validation and selection by qRT-PCR and univariate and multivariate analysis of the metastasis- and prognosis-related lncRNAs, the regulated microRNAs (miRNAs) and coexpressed mRNAs were used to construct a ceRNA network for distant metastasis of CRC.

RESULTS

Two key distant metastasis-related DElncRNAs, AP002498.1 and LINC01871, were identified by univariate and multivariate analysis in combination with analyses of clinical data and expression levels. Furthermore, lncRNA-associated ceRNA subnetworks were constructed from the predicted miRNAs and 13 coexpressed DEmRNAs (SERPINA1, ITLN1, REG4, L1TD1, IGFALS, MUC5B, CIITA, CXCL9, CXCL10, GBP4, GNLY, IDO1, and NOS2). The AP002498.1- and LINC01871-associated ceRNA subnetworks regulated the expression of the target genes SERPINA1 and MUC5B and GNLY, respectively, through the associated miRNAs.

CONCLUSION

The DElncRNA AP002498.1 and the LINC01871/miR-4644 and miR-185-5p/GNLY axes were identified as being closely associated with distant metastasis and could represent independent prognostic biomarkers or therapeutic targets in colorectal adenocarcinoma.

Comprehensive analysis of basement membrane and immune checkpoint related lncRNA and its prognostic value in hepatocellular carcinoma via machine learning

Background

Hepatocellular carcinoma (HCC), which is characterized by its high malignancy, generally exhibits poor response to immunotherapy. As part of the tumor microenvironment, basement membranes (BMs) are involved in tumor development and immune activities. Presently, there is no integrated analysis linking the basement membrane with immune checkpoints, especially from the perspective of lncRNA.

Methods

Based on transcriptome data from The Cancer Genome Atlas, BMs-related and immune checkpoint-related lncRNAs were identified. By applying univariable Cox regression and Machine learning (LASSO and SVM-RFE algorithm), a 10-lncRNA prognosis signature was constructed. The prognostic significance of this signature was assessed by survival analysis. GSEA, ssGSEA, and drug sensitivity analysis were conducted to investigate potential functional pathways, immune status, and clinical implications of guiding individual treatments in HCC. Finally, the promoting migration effect of LINC01224 was validated via in vitro experiments.

Results

The multiple Cox regression, receiver operating characteristic curves, and stratified survival analysis of clinical subgroups exhibited the robust prognostic ability of the lncRNA signature. Results of the GSEA and drug sensitivity analysis revealed significant differences in potential functional pathways and response to drugs between the two risk groups. In addition, the risk level of HCC patients was distinctly correlated with immune cell infiltration status. More importantly, LINC01224 was independently associated with the OS of HCC patients (P < 0.05), suppressing the expression of LINC01224 inhibited the migration of HCC cells.

Conclusion

This study developed a reliable signature for the prognosis of HCC based on BM and immune checkpoint related lncRNA, revealing that LINC01224 might be a prognostic biomarker for HCC associated with the progression of HCC.

Epigenome-wide methylation study identified two novel CpGs associated with T2DM risk and a network of co-methylated CpGs capable of patient's classifications

Prevention of Type 2 diabetes mellitus (T2DM) pandemic needs markers that can precisely predict the disease risk in an individual. Alterations in DNA methylations due to exposure towards environmental risk factors are widely sought markers for T2DM risk prediction. To identify such individual DNA methylation signatures and their effect on disease risk, we performed an epigenome-wide association study (EWAS) in 844 Indian individuals of Indo-European origin. We identified and validated methylation alterations at two novel CpG sites in MIR1287 (cg01178710) and EDN2-SCMH1 (cg04673737) genes associated with T2DM risk at the epigenome-wide-significance-level (P < 1.2 × 10-7). Further, we also replicated the association of two known CpG sites in TXNIP, and CPT1A in the Indian population. With 535 EWAS significant CpGs (P < 1.2 × 10-7) identified in the discovery phase samples, we created a co-methylation network using weighted correlation network analysis and identified four modules among the CpGs. We observed that methylation of one of the module associates with T2DM risk factors (e.g. BMI, insulin and C-peptide) and can be used as markers to segregate T2DM patients with good glycemic control (e.g. low HbA1c) and dyslipidemia (low HDL and high TG) from the other patients. Additionally, an intronic SNP (rs6503650) in the JUP gene, a member of the same module, associated with methylation at all the 14 hub CpG sites of that module as methQTL. Our network-assisted EWAS is the first to systematically explore DNA methylation variations conferring risks to T2DM in Indians and use the identified risk CpG sites for patient segregation with different clinical outcomes. These findings can be useful for better stratification of patients to improve the clinical management and treatment effects.

Predicting prognosis, immunotherapy and distinguishing cold and hot tumors in clear cell renal cell carcinoma based on anoikis-related lncRNAs

Background

Clear cell renal cell carcinoma (ccRCC) is the most frequently occurring malignant tumor within the kidney cancer subtype. It has low sensitivity to traditional radiotherapy and chemotherapy, the optimal treatment for localized ccRCC has been surgical resection, but even with complete resection the tumor will be eventually developed into metastatic disease in up to 40% of localized ccRCC. For this reason, it is crucial to find early diagnostic and treatment markers for ccRCC.

Methods

We obtained anoikis-related genes (ANRGs) integrated from Genecards and Harmonizome dataset. The anoikis-related risk model was constructed based on 12 anoikis-related lncRNAs (ARlncRNAs) and verified by principal component analysis (PCA), Receiver operating characteristic (ROC) curves, and T-distributed stochastic neighbor embedding (t-SNE), and the role of the risk score in ccRCC immune cell infiltration, immune checkpoint expression levels, and drug sensitivity was evaluated by various algorithms. Additionally, we divided patients based on ARlncRNAs into cold and hot tumor clusters using the ConsensusClusterPlus (CC) package.

Results

The AUC of risk score was the highest among various factors, including age, gender, and stage, indicating that the model we built to predict survival was more accurate than the other clinical features. There was greater sensitivity to targeted drugs like Axitinib, Pazopanib, and Sunitinib in the high-risk group, as well as immunotherapy drugs. This shows that the risk-scoring model can accurately identify candidates for ccRCC immunotherapy and targeted therapy. Furthermore, our results suggest that cluster 1 is equivalent to hot tumors with enhanced sensitivity to immunotherapy drugs.

Conclusion

Collectively, we developed a risk score model based on 12 prognostic lncRNAs, expected to become a new tool for evaluating the prognosis of patients with ccRCC, providing different immunotherapy strategies by screening for hot and cold tumors.

Transcriptome profiling of skeletal muscles from Korean patients with Bethlem myopathy

Bethlem myopathy is one of the collagens VI-related muscular dystrophies caused by mutations in the collagen VI genes. The study was designed to analyze the gene expression profiles in the skeletal muscle of patients with Bethlem myopathy. Six skeletal muscle samples from 3 patients with Bethlem myopathy and 3 control subjects were analyzed by RNA-sequencing. 187 transcripts were significantly differentially expressed, with 157 upregulated and 30 downregulated transcripts in the Bethlem group. Particularly, 1 (microRNA-133b) was considerably upregulated, and 4 long intergenic non-protein coding RNAs, LINC01854, MBNL1-AS1, LINC02609, and LOC728975, were significantly downregulated. We categorized differentially expressed gene using Gene Ontology and showed that Bethlem myopathy is strongly associated with the organization of extracellular matrix (ECM). Kyoto Encyclopedia of Genes and Genomes pathway enrichment reflected themes with significant enrichment of the ECM-receptor interaction (hsa04512), complement and coagulation cascades (hsa04610), and focal adhesion (hsa04510). We confirmed that Bethlem myopathy is strongly associated with the organization of ECM and the wound healing process. Our results demonstrate transcriptome profiling of Bethlem myopathy, and provide new insights into the path mechanism of Bethlem myopathy associated with non-protein coding RNAs.

Development and validation of a novel ferroptosis‑related lncRNA prognostic signature for pancreatic adenocarcinoma

Long non‑coding RNAs (lncRNAs) serve a pivotal role in the regulation of cancer cell ferroptosis. However, the prognostic value of ferroptosis‑related lncRNAs in pancreatic adenocarcinoma (PAAD) largely remains unclear. We aimed at constructing a lncRNA‑based signature to improve the prognosis prediction of PAAD. In the present study, the transcriptome profiling data and clinical information of patients with PAAD were obtained from The Cancer Genome Atlas (TCGA) and International Cancer Gene Consortium (ICGC) databases. Univariate Cox regression analysis of the TCGA cohort demonstrated that 26 ferroptosis‑related lncRNAs had significant prognostic value for PAAD (all P<0.01). Least absolute shrinkage and selection operator regression and multivariate Cox proportional hazards regression analyses were performed to construct a prognostic ferroptosis‑related lncRNA signature (FRLS) comprising nine ferroptosis‑related lncRNAs. The efficacy of this FRLS was verified in the training (TCGA) and validation (ICGC) cohorts. Based on the risk model, high risk scores were significantly correlated with poor overall survival (OS) (hazard ratio, 1.314; 95% confidence interval, 1.218‑1.418; P<0.001). The receiver operating characteristic curves and principal component analysis further demonstrated the robust prognostic ability of the FRLS. Furthermore, a nomogram with favorable predictive efficacy for the prediction of OS was constructed based on the FRLS and clinical features. Gene set enrichment analysis demonstrated that the genes in the FRLS participated in a number of cancer‑associated immunoregulatory pathways. Importantly, it was demonstrated that immune infiltration and response to cancer immunotherapy differed significantly between the high and low‑risk groups according to the FRLS. In conclusion, the risk signature based on the FRLS has potential for the clinical prediction of prognosis and immunotherapy response in patients with PAAD.

Prediction of the Prognosis of Clear Cell Renal Cell Carcinoma by Cuproptosis-Related lncRNA Signals Based on Machine Learning and Construction of ceRNA Network

Background

Clear cell renal cell carcinoma's (ccRCC) occurrence and development are strongly linked to the metabolic reprogramming of tumors, and thus far, neither its prognosis nor treatment has achieved satisfying clinical outcomes.

Methods

The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, respectively, provided us with information on the RNA expression of ccRCC patients and their clinical data. Cuproptosis-related genes (CRGS) were discovered in recent massive research. With the help of log-rank testing and univariate Cox analysis, the prognostic significance of CRGS was examined. Different cuproptosis subtypes were identified using consensus clustering analysis, and GSVA was used to further investigate the likely signaling pathways between various subtypes. Univariate Cox, least absolute shrinkage and selection operator (Lasso), random forest (RF), and multivariate stepwise Cox regression analysis were used to build prognostic models. After that, the models were verified by means of the C index, Kaplan-Meier (K-M) survival curves, and time-dependent receiver operating characteristic (ROC) curves. The association between prognostic models and the tumor immune microenvironment as well as the relationship between prognostic models and immunotherapy were next examined using ssGSEA and TIDE analysis. Four online prediction websites-Mircode, MiRDB, MiRTarBase, and TargetScan-were used to build a lncRNA-miRNA-mRNA ceRNA network.

Results

By consensus clustering, two subgroups of cuproptosis were identified that represented distinct prognostic and immunological microenvironments.

Conclusion

A prognostic risk model with 13 CR-lncRNAs was developed. The immune microenvironment and responsiveness to immunotherapy are substantially connected with the model, which may reliably predict the prognosis of patients with ccRCC.

The Role of Long Noncoding RNA (lncRNAs) Biomarkers in Renal Cell Carcinoma

Renal cell carcinoma is one of the common cancers whose incidence and mortality are continuously growing worldwide. Initially, this type of tumour is usually asymptomatic. Due to the lack of reliable diagnostic markers, one-third of ccRCC patients already have distant metastases at the time of diagnosis. This underlines the importance of establishing biomarkers that would enable the prediction of the disease's course and the risk of metastasis. LncRNA, which modulates genes at the epigenetic, transcriptional, and post-transcriptional levels, appears promising. The actions of lncRNA involve sponging and sequestering target miRNAs, thus affecting numerous biological processes. Studies have confirmed the involvement of RNAs in various diseases, including RCC. In this review, we focused on MALAT1 (a marker of serious pathological changes and a factor in the promotion of tumorigenesis), RCAT1 (tumour promoter in RCC), DUXAP9 (a plausible marker of localized ccRCC), TCL6 (exerting tumour-suppressive effects in renal cancer), LINC00342 (acting as an oncogene), AGAP2 Antisense1 (plausible predictor of RCC progression), DLEU2 (factor promoting tumours growth via the regulation of epithelial-mesenchymal transition), NNT-AS1 (sponge of miR-22 contributing to tumour progression), LINC00460 (favouring ccRCC development and progression) and Lnc-LSG1 (a factor that may stimulate ccRCC metastasis).

Identification and validation of a ferroptosis-related lncRNA signature to robustly predict the prognosis, immune microenvironment, and immunotherapy efficiency in patients with clear cell renal cell carcinoma

Background

Ferroptosis is a new type of iron- and reactive oxygen species-dependent cell death, studies on ferroptosis-related long noncoding RNAs (FerLncRNAs) in clear cell renal cell carcinoma (ccRCC) are limited. The purpose of this study was to investigate the potential prognostic value of FerLncRNAs and their relationship with the immune microenvironment and immunotherapy response of ccRCC.

Methods

RNA sequencing data of 526 patients with ccRCC were downloaded from The Cancer Genome Atlas (TCGA) database. The patients with ccRCC in TCGA were randomly divided (1:1) into a training and testing cohort. ICGC and GEO databases were used for validation. Screening for FerLncRNAs was performed using Pearson's correlation analysis with the reported ferroptosis-related genes. A FerLncRNA signature was constructed using univariate, LASSO, and multivariate Cox regression analyses in the training cohort. Internal and external datasets were performed to verify the FRlncRNA signature. Four major FRlncRNAs were verified through in vitro experiment.

Results

We identified seven FerLncRNAs (LINC00894, DUXAP8, LINC01426, PVT1, PELATON, LINC02609, and MYG1-AS1), and established a risk signature and nomogram for predicting the prognosis of ccRCC. Four major FRlncRNAs were verified with the prognosis of ccRCC in the GEPIA and K-M Plotter databases, and their expressions were validated by realtime PCR. The risk signature can also effectively reflect the immune environment, immunotherapy response and drug sensitivity of ccRCC. These FRlncRNAs have great significance to the implementation of individualized treatment and disease monitoring of ccRCC patients.

Qualification of Necroptosis-Related lncRNA to Forecast the Treatment Outcome, Immune Response, and Therapeutic Effect of Kidney Renal Clear Cell Carcinoma

Background

Kidney renal clear cell carcinoma (KIRC) is considered as a highly immune infiltrative tumor. Necroptosis is an inflammatory programmed cell death associated with a wide range of diseases. Long noncoding RNAs (lncRNAs) play important roles in gene regulation and immune function. lncRNA associated with necroptosis could systematically explore the prognostic value, regulate tumor microenvironment (TME), etc.

Method

The patients' data was collected from TCGA datasets. We used the univariate Cox regression (UCR) to select prediction lncRNAs that are related to necroptosis. Meanwhile, risk models were constructed using LASSO Cox regression (LCR). Kaplan–Meier (KM) analysis, accompanied with receiver operating characteristic (ROC) curves, was performed to assess the independent risk factors of different clinical characteristics. The evaluated factors are age, gender, disease staging, grade, and their related risk score. Databases such as Gene Ontology (GO), Kyoto encyclopedia of genes and genomes (KEGG), and Gene set enrichment analysis (GSEA) were used to search the probable biological characteristics that could influence the risk groups, containing signaling pathway and immue-related pathways. The single-sample gene set enrichment analysis (ssGSEA) was chosen to perform gene set variation analysis (GSVA), and the GSEABase package was selected to detect the immune and inflammatory infiltration profiles. The TIDE and IC₅₀ evaluation were used to estimate the effectiveness of clinical treatment on KIRC.

Results

Based on the above analysis, we have got a conclusion that patients who show high risk had higher immune infiltration, immune checkpoint expression, and poorer prognosis. We identified 19 novel prognostic necroptosis-related lncRNAs, which could offer opinions for a deeper study of KIRC.

Conclusion

The risk model we constructed makes it possible to predict the prognosis of KIRC patients and offers directions for further research on the prognostication and treatment strategies for KIRC.

Pyroptosis-Related lncRNA Prognostic Model for Renal Cancer Contributes to Immunodiagnosis and Immunotherapy

Background

Renal clear cell cancer (ccRCC) is one of the most common cancers in humans. Thus, we aimed to construct a risk model to predict the prognosis of ccRCC effectively.

Methods

We downloaded RNA sequencing (RNA-seq) data and clinical information of 539 kidney renal clear cell carcinoma (KIRC) patients and 72 normal humans from The Cancer Genome Atlas (TCGA) database and divided the data into training and testing groups randomly. Pyroptosis-related lncRNAs (PRLs) were obtained through Pearson correlation between pyroptosis genes and all lncRNAs (p &lt; 0.05, coeff &gt; 0.3). Univariate and multivariate Cox regression analyses were then performed to select suitable lncRNAs. Next, a novel signature was constructed and evaluated by survival analysis and ROC analysis. The same observation applies to the testing group to validate the value of the signature. By gene set enrichment analysis (GSEA), we predicted the underlying signaling pathway. Furthermore, we calculated immune cell infiltration, immune checkpoint, the T-cell receptor/B-cell receptor (TCR/BCR), SNV, and Tumor Immune Dysfunction and Exclusion (TIDE) scores in TCGA database. We also validated our model with an immunotherapy cohort. Finally, the expression of PRLs was validated by quantitative PCR (qPCR).

Results

We constructed a prognostic signature composed of six key lncRNAs (U62317.1, MIR193BHG, LINC02027, AC121338.2, AC005785.1, AC156455.1), which significantly predict different overall survival (OS) rates. The efficiency was demonstrated using the receiver operating characteristic (ROC) curve. The signature was observed to be an independent prognostic factor in cohorts. In addition, we found the PRLs promote the tumor progression via immune-related pathways revealed in GSEA. Furthermore, the TCR, BCR, and SNV data were retrieved to screen immune features, and immune cell scores were calculated to measure the effect of the immune microenvironment on the risk model, indicating that high- and low-risk scores have different immune statuses. The TIDE algorithm was then used to predict the immune checkpoint blockade (ICB) response of our model, and subclass mapping was used to verify our model in another immunotherapy cohort data. Finally, qPCR validates the PRLs in cell lines.

Conclusion

This study provided a new risk model to evaluate ccRCC and may be pyroptosis-related therapeutic targets in the clinic.

Decoding the Immune Microenvironment of Clear Cell Renal Cell Carcinoma by Single-Cell Profiling to Aid Immunotherapy

Clear cell renal cell carcinoma (ccRCC) is the most common subtype of kidney cancer, and it is the major cause of kidney cancer death. Understanding tumor immune microenvironments (TMEs) is critical in cancer immunotherapies. Here, we studied the immune characterization at single-cell resolution by integrating public data of ccRCC across different tissue types, and comparing the transcriptome features and tumor TME differences in tumors, normal adjacent tissue, and peripheral blood. A total of 16 different types of cell components of ccRCC were identified. We revealed that there is an overall increase in T-cell and myeloid populations in tumor-infiltrated immune cells compared to normal renal tissue, and the B-cell population in the tumor showed a sharp decrease, which indicates that the cells in tumor tissue undergo strong immune stress. In addition, the cell-cell communication analysis revealed specific or conserved signals in different tissue types, which may aid to uncover the distinct immune response. By combining and analyzing publicly available ccRCC bulk RNA-seq datasets, 10 genes were identified as marker genes in specific cell types, which were significantly associated with poor prognosis. Of note, UBE2C, which may be a good indicator of tumor proliferation, is positively associated with reductions in overall survival and highly associated with tumor grade. Our integrated analysis provides single-cell transcriptomic profiling of ccRCC and their TME, and it unmasked new correlations between gene expression, survival outcomes, and immune cell-type components, enabling us to dissect the dynamic variables in the tumor development process. This resource provides deeper insight into the transcriptome features and immune response of ccRCC and will be helpful in kidney cancer immunotherapy.

Identification of Five Ferroptosis-Related LncRNAs as Novel Prognosis and Diagnosis Signatures for Renal Cancer

Background: Ferroptosis is a novel regulated cell death that is characterized by iron-dependent oxidative damage. Renal cancer is the second most common cancer of the urinary system, which is highly correlated with iron metabolism. The aim of our present study was to identify suitable ferroptosis-related prognosis signatures for renal cancer.

Methods: We downloaded the RNA-seq count data of renal cancer from The Cancer Genome Atlas database and used the DESeq2, Survival, and Cox regression analyses to screen the prognosis signatures.

Results: We identified 5 ferroptosis-related differentially expressed lncRNAs (FR-DELs) (DOCK8-AS1, SNHG17, RUSC1-AS1, LINC02609, and LUCAT1) to be independently correlated with the overall survival (OS) of patients with renal cancer. The risk assessment model and diagnosis model constructed by those 5 FR-DELs could well predict the outcome and the diagnosis of renal cancer.

Conclusion: Our present study not only suggested those 5 FR-DELs could be used as prognosis and diagnosis signatures for renal cancer but also provided strategies for other cancers in the screening of ferroptosis-related biomarkers.

Identification of a competing endogenous RNA network related to immune signature in clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is a fatal cancer of the urinary system. Long non-coding RNAs (lncRNAs) act as competitive endogenous RNAs (ceRNAs) involving the ccRCC progression. However, the relationship between the ceRNA network and immune signature is largely unknown. In this study, the ccRCC-related gene expression profiles retrieved from the TCGA database were used first to identify the differentially expressed genes through differential gene expression analysis and weighted gene co-expression network analysis. The interaction among differentially expressed lncRNAs, miRNAs, and mRNAs were matched using public databases. As a result, a ceRNA network was developed that contained 144 lncRNAs, 23 miRNAs, as well as 62 mRNAs. Four of 144 lncRNAs including LINC00943, SRD5A3-AS1, LINC02345, and U62317.3 were identified through LASSO regression and Cox regression analyses, and were used to create a prognostic risk model. Then, the ccRCC samples were divided into the high- and low-risk groups depending on their risk scores. ROC curves, Kaplan-Meier survival analysis, and the survival risk plots indicated that the predictive performance of our developed risk model was accurate. Moreover, the CIBERSORT algorithm was used to measure the infiltration levels of immune cells in the ccRCC samples. The further genomic analysis illustrated a positive correlation between most immune checkpoint blockade-related genes and the risk score. In conclusion, the present findings effectually contribute to the comprehensive understanding of the ccRCC pathogenesis, and may offer a reference for developing novel therapeutic and prognostic biomarkers.

An immune-related lncRNA risk coefficient model to predict the outcomes in clear cell renal cell carcinoma

Objective: Using model algorithms, we constructed an immune-related long non-coding RNAs (lncRNAs) risk coefficient model to predict outcomes for patients with clear cell renal cell carcinoma (ccRCC) to understand the infiltration of tumor immune cells and the sensitivity to immune-targeted drugs.

Methods: Open genes data were downloaded from The Cancer Genome Atlas and The Immunology Database and Analysis Portal, and immune-related lncRNAs were obtained through Pearson correlation analysis. R language software was used to obtain differentially expressed immune-related lncRNAs and immune-related lncRNA pairs. The model was constructed using least absolute shrinkage and selector operation regression analysis, and receiver operator characteristic curves were drawn. The Akaike information criterion was used to distinguish the high-risk from the low-risk group. We also conducted correlation analysis for the high- and low-risk subgroups.

Results: We identified 27 immune-related lncRNAs pairs, 16 of which were included in the model construction. After merging clinical data, the areas under the curve of 1 -year, 3-year, and 5-year survival times of ccRCC patients were 0.867, 0.832, and 0.838, respectively. Subgroup analyses were conducted according to the cut-off value. We found that the high-risk group was associated with poor outcomes. The risk score and tumor stage were independent predictors of the outcome of ccRCC. The risk model predicted specific immune cell infiltration, immune checkpoint gene expression levels, and high-risk groups more sensitive to sunitinib targeted therapy.

Conclusion: We obtained prognostic-related novel ccRCC markers and risk model that predicts the outcome of patients with ccRCC and helps identify those who can benefit from sunitinib.

Development and validation of ferroptosis-related lncRNAs prognosis signatures in kidney renal clear cell carcinoma

Background

Ferroptosis is a recently recognised new type of cell death which may be a potential target for cancer therapy. In the present study, we aimed to screen ferroptosis-related differentially expressed long non-coding RNAs as biomarkers to predict the outcome of kidney renal clear cell carcinoma.

Methods

RNAseq count data and corresponding clinical information were obtained from the Cancer Genome Atlas database. Lists of ferroptosis-related genes and long non-coding RNAs were obtained from the FerrDb and GENCODE databases, respectively. The candidate prognostic signatures were screened by Cox regression analyses and least absolute shrinkage and selection operator analyses.

Results

Three ferroptosis-related long non-coding RNAs (DUXAP8, LINC02609, and LUCAT1) were significantly correlated with the overall survival of kidney renal clear cell carcinoma independently. Kidney renal clear cell carcinoma patients with high-risk values displayed worse OS. Meanwhile, the expression of these three ferroptosis-related long non-coding RNAs and their risk scores were significantly correlated with clinicopathological features. Principal component analyses showed that patients with kidney renal clear cell carcinoma have differential risk values were well distinguished by the three ferroptosis-related long non-coding RNAs.

Conclusions

The present study suggests that the risk assessment model constructed by these three ferroptosis-related long non-coding RNAs could accurately predict the outcome of kidney renal clear cell carcinoma. We also provide a novel perspective for cancer prognosis screening.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02284-1.

Angiogenesis-Related Molecular Subtypes and a Novel Prognostic Signature in Clear Cell Renal Cell Carcinoma Patients

Background

This study aimed to develop and validate a novel angiogenesis-related gene (ARG) signature and molecular subtypes by bioinformatics analysis.

Materials and Methods

The transcriptome data and clinical data were obtained from TCGA and ICGC database. We performed consensus clustering analysis to identify angiogenesis molecular subtypes for ccRCC. Univariate and multivariate Cox regression analyses were used to develop a novel ARG-related signature as a prognostic biomarker for ccRCC. Internal and external validation were then performed in TCGA and ICGC cohort, respectively.

Results

We identified a total of two angiogenesis molecular subtypes of ccRCC. The overall survival (OS) of subtype 1 ccRCC was significantly decreased compared with that of subtype 2 ccRCC (P=0.001). These two molecular subtypes have significantly different tumor microenvironment and immune checkpoint inhibitor sensitivities (P<0.05). Besides, we developed a novel signature based on three ARGs (including MSX1, TIMP1 and JAG2) for subtype 1 ccRCC. The difference in OS between high- and low-risk group was statistically significant in training cohort (P=0.009), test cohort (P=0.024), the whole type 1 cohort (P<0.001), and validation cohort (P=0.041). The AUC for one-year OS prediction was 0.732, 0.710, 0.725, and 0.645 in training cohort, test cohort, the whole type 1 cohort, and validation cohort, respectively. Independent prognostic analysis showed that this signature was an independent predictor for OS of subtype 1 ccRCC (P=0.028914). The power of this prognostic signature was superior to other signatures reported in previous studies.

Conclusion

We developed and successfully validated a novel ARG signature for predicting prognosis of subtype 1 ccRCC, which was superior to several previous signatures.