Identification of 12 immune-related lncRNAs and molecular subtypes for the clear cell renal cell carcinoma based on RNA sequencing data

Screening of differential gene expression patterns through survival analysis for diagnosis, prognosis and therapies of clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is the most prevalent subtype of kidney cancer. Although there is increasing evidence linking ccRCC to genetic alterations, the exact molecular mechanism behind this relationship is not yet completely known to the researchers. Though drug therapies are the best choice after the metastasis, unfortunately, the majority of the patients progressively develop resistance against the therapeutic drugs after receiving it for almost 2 years. In this case, multi-targeted different variants of therapeutic drugs are essential for effective treatment against ccRCC. To understand molecular mechanisms of ccRCC development and progression, and explore multi-targeted different variants of therapeutic drugs, it is essential to identify ccRCC-causing key genes (KGs). In order to obtain ccRCC-causing KGs, at first, we detected 133 common differentially expressed genes (cDEGs) between ccRCC and control samples based on nine (9) microarray gene-expression datasets with NCBI accession IDs GSE16441, GSE53757, GSE66270, GSE66272, GSE16449, GSE76351, GSE66271, GSE71963, and GSE36895. Then, we filtered these cDEGs through survival analysis with the independent TCGA and GTEx database and obtained 54 scDEGs having significant prognostic power. Next, we used protein-protein interaction (PPI) network analysis with the reduced set of 54 scDEGs to identify ccRCC-causing top-ranked eight KGs (PLG, ENO2, ALDOB, UMOD, ALDH6A1, SLC12A3, SLC12A1, SERPINA5). The pan-cancer analysis with KGs based on TCGA database showed the significant association with different subtypes of kidney cancers including ccRCC. The gene regulatory network (GRN) analysis revealed some crucial transcriptional and post-transcriptional regulators of KGs. The scDEGs-set enrichment analysis significantly identified some crucial ccRCC-causing molecular functions, biological processes, cellular components, and pathways that are linked to the KGs. The results of DNA methylation study indicated the hypomethylation and hyper-methylation of KGs, which may lead the development of ccRCC. The immune infiltrating cell types (CD8+ T and CD4+ T cell, B cell, neutrophil, dendritic cell and macrophage) analysis with KGs indicated their significant association in ccRCC, where KGs are positively correlated with CD4+ T cells, but negatively correlated with the majority of other immune cells, which is supported by the literature review also. Then we detected 10 repurposable drug molecules (Irinotecan, Imatinib, Telaglenastat, Olaparib, RG-4733, Sorafenib, Sitravatinib, Cabozantinib, Abemaciclib, and Dovitinib.) by molecular docking with KGs-mediated receptor proteins. Their ADME/T analysis and cross-validation with the independent receptors, also supported their potent against ccRCC. Therefore, these outputs might be useful inputs/resources to the wet-lab researchers and clinicians for considering an effective treatment strategy against ccRCC.

A novel autophagy-related long non-coding RNAs prognostic risk score for clear cell renal cell carcinoma

BACKGROUND

As the main histological subtype of renal cell carcinoma, clear cell renal cell carcinoma (ccRCC) places a heavy burden on health worldwide. Autophagy-related long non-coding RNAs (ARlncRs) have shown tremendous potential as prognostic signatures in several studies, but the relationship between them and ccRCC still has to be demonstrated.

METHODS

The RNA-sequencing and clinical characteristics of 483 ccRCC patients were downloaded download from the Cancer Genome Atlas and International Cancer Genome Consortium. ARlncRs were determined by Pearson correlation analysis. Univariate and multivariate Cox regression analyses were applied to establish a risk score model. A nomogram was constructed considering independent prognostic factors. The Harrell concordance index calibration curve and the receiver operating characteristic analysis were utilized to evaluate the nomogram. Furthermore, functional enrichment analysis was used for differentially expressed genes between the two groups of high- and low-risk scores.

RESULTS

A total of 9 SARlncRs were established as a risk score model. The Kaplan-Meier survival curve, principal component analysis, and subgroup analysis showed that low overall survival of patients was associated with high-risk scores. Age, M stage, and risk score were identified as independent prognostic factors to establish a nomogram, whose concordance index in the training cohort, internal validation, and external ICGC cohort was 0.793, 0.671, and 0.668 respectively. The area under the curve for 5-year OS prediction in the training cohort, internal validation, and external ICGC cohort was 0.840, 0.706, and 0.708, respectively. GO analysis and KEGG analysis of DEGs demonstrated that immune- and inflammatory-related pathways are likely to be critically involved in the progress of ccRCC.

CONCLUSIONS

We established and validated a novel ARlncRs prognostic risk model which is valuable as a potential therapeutic target and prognosis indicator for ccRCC. A nomogram including the risk model is a promising clinical tool for outcomes prediction of ccRCC patients and further formulation of individualized strategy.

Novel Prognosis and Therapeutic Response Model of Immune-Related lncRNA Pairs in Clear Cell Renal Cell Carcinoma

Clear cell renal cell carcinoma (ccRCC) is the most common type of renal carcinoma. It is particularly important to accurately judge the prognosis of patients. Since most tumor prediction models depend on the specific expression level of related genes, a better model therefore needs to be constructed. To provide an immune-related lncRNA (irlncRNAs) tumor prognosis model that is independent of the specific gene expression levels, we first downloaded and sorted out the data on ccRCC in the TCGA database and screened irlncRNAs using co-expression analysis and then obtained the differently expressed irlncRNA (DEirlncRNA) pairs by means of univariate analysis. In addition, we modified LASSO penalized regression. Subsequently, the ROC curve was drawn, and we compared the area under the curve, calculated the Akaike information standard value of the 5-year receiver operating characteristic curve, and determined the cut-off point to establish the best model to distinguish the high- or low-disease-risk group of ccRCC. Subsequently, we reassessed the model from the perspectives of survival, clinic-pathological characteristics, tumor-infiltrating immune cells, chemotherapeutics efficacy, and immunosuppressed biomarkers. A total of 17 DEirlncRNAs pairs (AL031710.1|AC104984.5, AC020907.4|AC127-24.4,AC091185.1|AC005104.1, AL513218.1|AC079015.1, AC104564.3|HOXB-AS3, AC003070.1|LINC01355, SEMA6A-AS1|CR936218.1, AL513327.1|AS005785.1, AC084876.1|AC009704.2, IGFL2-AS1|PRDM16-DT, AC011462.4|MMP25-AS1, AL662844.3I|TGB2-AS1, ARHGAP27P1|AC116914.2, AC093788.1|AC007098.1, MCF2L-AS1|AC093001.1, SMIM25|AC008870.2, and AC027796.4|LINC00893) were identified, all of which were included in the Cox regression model. Using the cut-off point, we can better distinguish patients according to different factors, such as survival status, invasive clinic-pathological features, tumor immune infiltration, whether they are sensitive to chemotherapy or not, and expression of immunosuppressive biomarkers. We constructed the irlncRNA model by means of pairing, which can better eliminate the dependence on the expression level of the target genes. In other words, the signature established by pairing irlncRNA regardless of expression levels showed promising clinical prediction value.

Construction of immune-related lncRNA signature to predict aggressiveness, immune landscape, and drug resistance of colon cancer

BACKGROUND

Colon cancer remains one of the most common malignancies across the world. Thus far, a biomarker, which can comprehensively predict the survival outcomes, clinical characteristics, and therapeutic sensitivity, is still lacking.

METHODS

We leveraged transcriptomic data of colon cancer from the existing datasets and constructed immune-related lncRNA (irlncRNA) pairs. After integrating with clinical survival data, we performed differential analysis and identified 11 irlncRNAs signature using Lasso regression analysis. We next plotted the 1-, 5-, and 10-year curve lines of receiver operating characteristics, calculated the areas under the curve, and recognized the optimal cutoff point. Then, we validated the pair-risk model in terms of the survival outcomes of the patients involved. Moreover, we tested the reliability of the model for predicting tumor aggressiveness and therapeutic susceptibility of colon cancer. Additionally, we reemployed the 11 of irlncRNAs involved in the pair-risk model to construct an expression-risk model to predict the prognostic outcomes of the patients involved.

RESULTS

We recognized a total of 377 differentially expressed irlncRNAs (DEirlcRNAs), including 28 low-expressed and 349 high-expressed irlncRNAs in colon cancer patients. After performing a univariant Cox analysis, we identified 115 risk irlncRNAs that were significantly correlated with survival outcomes of patients involved. By taking the overlap of the DEirlcRNAs and the risk irlncRNAs, we ultimately recognized 55 irlncRNAs as core irlncRNAs. Then, we established a Cox HR model (pair-risk model) as well as an expression HR model (exp-risk model) based on 11 of the 55 core irlncRNAs. We found that both of the two models significantly outperformed the commonly used clinical characteristics, including age, T, N, and M stages when predicting survival outcomes. Moreover, we validated the pair-risk model as a potential tool for studying the tumor microenvironment of colon cancer and drug susceptibility. Additionally, we noticed that combinational use of the pair-risk model and the exp-risk model yielded a more robust approach for predicting the survival outcomes of patients with colon cancer.

CONCLUSIONS

We recognized 11 irlncRNAs and created a pair-risk model and an exp-risk model, which have the potential to predict clinical characteristics of colon cancer, either solely or conjointly.

ISPRF: a machine learning model to predict the immune subtype of kidney cancer samples by four genes

Background

Clear cell renal cell carcinoma (ccRCC) is the most common type of renal cell carcinoma (RCC). Immunotherapy, especially anti-PD-1, is becoming a pillar of ccRCC treatment. However, precise biomarkers and robust models are needed to select the proper patients for immunotherapy.

Methods

A total of 831 ccRCC transcriptomic profiles were obtained from 6 datasets. Unsupervised clustering was performed to identify the immune subtypes among ccRCC samples based on immune cell enrichment scores. Weighted correlation network analysis (WGCNA) was used to identify hub genes distinguishing subtypes and related to prognosis. A machine learning model was established by a random forest (RF) algorithm and used on an open and free online website to predict the immune subtype.

Results

In the identified immune subtypes, subtype2 was enriched in immune cell enrichment scores and immunotherapy biomarkers. WGCNA analysis identified four hub genes related to immune subtypes, CTLA4, FOXP3, IFNG, and CD19. The RF model was constructed by mRNA expression of these four hub genes, and the value of area under the receiver operating characteristic curve (AUC) was 0.78. Subtype2 patients in the independent validation cohort had a better drug response and prognosis for immunotherapy treatment. Moreover, an open and free website was developed by the RF model (https://immunotype.shinyapps.io/ISPRF/).

Conclusions

The current study constructs a model and provides a free online website that could identify suitable ccRCC patients for immunotherapy, and it is an important step forward to personalized treatment.

Identification of N6-Methyladenosine-Related LncRNAs for Predicting Overall Survival and Clustering of a Potentially Novel Molecular Subtype of Breast Cancer

We aimed to identify a signature comprising N6-methyladenosine (m6A)-related long non-coding RNAs (lncRNAs) and molecular subtypes associated with breast cancer (BRCA). We obtained data of BRCA samples from The Cancer Genome Atlas database. The m6A-related lncRNA prognostic signature (m6A-LPS) included 10 lncRNAs previously identified as prognostic m6A-related lncRNAs and was constructed using integrated bioinformatics analysis and validated. Accordingly, a risk score based on the m6A-LPS signature was established and shown to confirm differences in survival between high-risk and low-risk groups. Three distinct genotypes were identified, whose characteristics included features of the tumor immune microenvironment in each subtype. Our results indicated that patients in Cluster 2 might have a worse prognostic outcome than those in other clusters. The three genotypes and risk subgroups were enriched in different biological processes and pathways, respectively. We then constructed a competing endogenous RNA network based on the prognostic m6A-related lncRNAs. Finally, we validated the expression levels of target lncRNAs in 72 clinical samples. In summary, the m6A-LPS and the potentially novel genotype may provide a theoretical basis for further study of the molecular mechanism of BRCA and may provide novel insights into precision medicine.

Circulating RNA in Kidney Cancer: What We Know and What We Still Suppose

Renal cancer represents the 7th most common tumor worldwide, affecting 400,000 people annually. This malignancy, which is the third most frequent cancer among urological diseases, displays a completely different prognosis if the tumor is detected in the early stages or advance phases. Unfortunately, more than 50% of renal cancers are discovered incidentally, with a consistent percentage of cases where the tumor remains clinically silent till the metastatic process is established. In day-to-day clinical practice, no available predictive biomarkers exist, and the existent imaging diagnostic techniques harbor several gaps in terms of diagnosis and prognosis. In the last decade, many efforts have been reported to detect new predictive molecular biomarkers using liquid biopsies, which are less invasive in comparison to renal biopsy. However, until now, there has been no clear evidence that a liquid biopsy biomarker could be relevant to the creation of a precise and tailored medical management in these oncological patients, even though circulating RNA biomarkers remain among the most promising. Given the idea that liquid biopsies will play a future key role in the management of these patients, in the present review, we summarize the current state of circulating RNA (miRNA, lncRNAs, and circRNAs) as possible biomarkers of renal cancer presence and aggressiveness in patients.