Prognosis of clear cell renal cell carcinoma (ccRCC) based on a six-lncRNA-based risk score: an investigation based on RNA-sequencing data

Boosting Clear Cell Renal Carcinoma-Specific Drug Discovery Using a Deep Learning Algorithm and Single-Cell Analysis

Clear cell renal carcinoma (ccRCC), the most common subtype of renal cell carcinoma, has the high heterogeneity of a highly complex tumor microenvironment. Existing clinical intervention strategies, such as target therapy and immunotherapy, have failed to achieve good therapeutic effects. In this article, single-cell transcriptome sequencing (scRNA-seq) data from six patients downloaded from the GEO database were adopted to describe the tumor microenvironment (TME) of ccRCC, including its T cells, tumor-associated macrophages (TAMs), endothelial cells (ECs), and cancer-associated fibroblasts (CAFs). Based on the differential typing of the TME, we identified tumor cell-specific regulatory programs that are mediated by three key transcription factors (TFs), whilst the TF EPAS1/HIF-2α was identified via drug virtual screening through our analysis of ccRCC's protein structure. Then, a combined deep graph neural network and machine learning algorithm were used to select anti-ccRCC compounds from bioactive compound libraries, including the FDA-approved drug library, natural product library, and human endogenous metabolite compound library. Finally, five compounds were obtained, including two FDA-approved drugs (flufenamic acid and fludarabine), one endogenous metabolite, one immunology/inflammation-related compound, and one inhibitor of DNA methyltransferase (N4-methylcytidine, a cytosine nucleoside analogue that, like zebularine, has the mechanism of inhibiting DNA methyltransferase). Based on the tumor microenvironment characteristics of ccRCC, five ccRCC-specific compounds were identified, which would give direction of the clinical treatment for ccRCC patients.

The lncRNA epigenetics: The significance of m6A and m5C lncRNA modifications in cancer

Most of our transcribed RNAs are represented by non-coding sequences. Long non-coding RNAs (lncRNAs) are transcripts with no or very limited protein coding ability and a length >200nt. They can be epigenetically modified. N6-methyladenosine (m6A), N1-methyladenosine (m1A), 5-methylcytosine (m5C), 7-methylguanosine (m7G) and 2’-O-methylation (Nm) are some of the lncRNAs epigenetic modifications. The epigenetic modifications of RNA are controlled by three classes of enzymes, each playing a role in a specific phase of the modification. These enzymes are defined as “writers”, “readers” and “erasers”. m6A and m5C are the most studied epigenetic modifications in RNA. These modifications alter the structure and properties, thus modulating the functions and interactions of lncRNAs. The aberrant expression of several lncRNAs is linked to the development of a variety of cancers and the epigenetic signatures of m6A- or m5C-related lncRNAs are increasingly recognized as potential biomarkers of prognosis, predictors of disease stage and overall survival. In the present manuscript, the most up to date literature is reviewed with the focus on m6A and m5C modifications of lncRNAs and their significance in cancer.

Molecular subtypes based on metabolic genes are potential biomarkers for predicting prognosis and immune responses of clear cell renal cell carcinoma

Due to the existence of tumor molecular heterogeneity, even patients having similar clinicopathological features could have vastly different survival rates. Hence, we aimed to explore novel metabolism-associated genes (MAGs) related molecular subtypes for clear cell renal cell carcinoma (ccRCC) and their immune landscapes for predicting prognosis and immune responses. Gene matrices and clinical information were downloaded from TCGA and ICGC datasets. Consensus clustering was conducted by the R "ConsensusClusterPlus" package. ccRCC patients were successfully divided into three clusters (MC1, MC2, and MC3) based on MAGs in both TCGA and ICGC datasets. Our established three MAGs were significantly associated with chemokine/chemokine receptor, IFN, CYT, angiogenesis, immune checkpoint molecules, tumor-infiltrating immune cells, oncogenic pathways, pan-cancer immune subtypes, and tumor microenvironment (TME) scores or expressions. Moreover, these three metabolic ccRCC subtypes could predict immunotherapeutic responses. We further constructed a characteristic index (LDAscore) in three metabolic ccRCC subtypes and identified LDAscore-related modules by WGCNA. After deep data mining, 10 hub genes were obtained and seven genes (ATRX, BPTF, DHX9, EP300, POLR2B, SIN3A, UBE3A) were finally validated by qRT-PCR. Our results successfully established a novel ccRCC subtype based on MAGs, providing novel insights into metabolism-related ccRCC tumor heterogeneity and facilitating individualized therapy for future work.

Development of genomic instability-associated long non-coding RNA signature: A prognostic risk model of clear cell renal cell carcinoma

Purpose

Renal clear cell carcinoma (ccRCC) is the most lethal of all pathological subtypes of renal cell carcinoma (RCC). Genomic instability was recently reported to be related to the occurrence and development of kidney cancer. The biological roles of long non-coding RNAs (lncRNAs) in tumorigenesis have been increasingly valued, and various lncRNAs were found to be oncogenes or cancer suppressors. Herein, we identified a novel genomic instability-associated lncRNA (GILncs) model for ccRCC patients to predict the overall survival (OS).

Methods

The Cancer Genome Atlas (TCGA) database was utilized to obtain full transcriptome data, somatic mutation profiles, and clinical characteristics. The differentially expressed lncRNAs between the genome-unstable-like group (GU) and the genome-stable-like group (GS) were defined as GILncs, with |logFC| > 1 and an adjusted p-value< 0.05 for a false discovery rate. All samples were allocated into GU-like or GS-like types based on the expression of GILncs observed using hierarchical cluster analyses. A genomic instability-associated lncRNA signature (GILncSig) was constructed using parameters of the included lncRNAs. Quantitative real-time PCR analysis was used to detect the in vitro expression of the included lncRNAs. Validation of the risk model was performed by the log-rank test, time-dependent receiver operating characteristic (ROC) curves analysis, and multivariate Cox regression analysis.

Results

Forty-six lncRNAs were identified as GILncs. LINC00460, AL139351.1, and AC156455.1 were employed for GILncSig calculation based on the results of Cox analysis. GILncSig was confirmed as an independent predictor for OS of ccRCC patients. Additionally, it presented a higher efficiency and accuracy than other RCC prognostic models reported before.

Conclusion

GILncSig score was qualified as a critical indicator, independent of other clinical factors, for prognostic prediction of ccRCC patients.

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.

An Integrative Analysis Framework for Identifying the Prognostic Markers from Multidimensional RNA Data of Clear Cell Renal Cell Carcinoma

The altered regulatory status of long noncoding RNA (lncRNA), miRNA, and mRNA and their interactions play critical roles in tumor proliferation, metastasis, and progression, which ultimately influence cancer prognosis. However, there are limited studies of comprehensive identification of prognostic biomarkers from combined data sets of the three RNA types in the highly metastatic clear cell renal cell carcinoma (ccRCC). The current study employed an integrative analysis framework of functional genomics approaches and machine learning methods to the lncRNA, miRNA, and mRNA data and identified 16 RNAs (3 lncRNAs, 6 miRNAs, and 7 mRNAs) of prognostic value, with 9 of them novel. A 16 RNA-based score was established for prognosis prediction of ccRCC with significance (P < 0.0001). The area under the curve for the score model was 0.868 to 0.870 in the training cohort and 0.714 to 0.778 in the validation cohort. Construction of the lncRNA-miRNA-mRNA interaction network showed that the downstream mRNAs and upstream lncRNAs in the network initiated from the miRNA or lncRNA markers exhibit significant enrichment in functional classifications associated with cancer metastasis, proliferation, progression, or prognosis. The functional analysis provided clear support for the role of the RNA biomarkers in predicting cancer prognosis. This study provides promising biomarkers for predicting prognosis of ccRCC using multidimensional RNA data, and these findings are expected to facilitate potential clinical applications of the biomarkers.

A computational prognostic model of lncRNA signature for clear cell renal cell carcinoma with genome instability

Long non-coding RNAs (lncRNAs) play a critical role in genomic instability and prognosis of cancer patients, but the methods to identify genomic instability-related lncRNAs have yet to be established. In the present study, to assess the prognostic value of lncRNAs associated with genomic instability in clear cell renal cell carcinoma (ccRCC).A computational framework was established based on the mutation hypothesis and combined lncRNA expression and somatic mutation profiles of the ccRCC genome. Furthermore, a prognostic model was developed using the genome instability-derived lncRNA signature GILncSig based on three lncRNA genes (LINC02471, LINC01234, and LINC00460) and verified using multiple independent patient cohorts.This study established an effective computational method to study the role of lncRNAs in genomic instability, with potential applications in identifying new genomic instability-related cancer biomarkers.

N6-Methyladenosine-Related Long Non-coding RNA Signature Associated With Prognosis and Immunotherapeutic Efficacy of Clear-Cell Renal Cell Carcinoma

Increasing evidence suggests that N6-methyladenosine (m6A) and long non-coding RNAs (lncRNAs) play important roles in cancer progression and immunotherapeutic efficacy in clear-cell renal cell carcinoma (ccRCC). In this study, we conducted a comprehensive ccRCC RNA-seq analysis using The Cancer Genome Atlas data to establish an m6A-related lncRNA prognostic signature (m6A-RLPS) for ccRCC. Forty-four prognostic m6A-related lncRNAs (m6A-RLs) were screened using Pearson correlation analysis (|R| > 0.7, p < 0.001) and univariable Cox regression analysis (p < 0.01). Using consensus clustering, the patients were divided into two clusters with different overall survival (OS) rates and immune status according to the differential expression of the lncRNAs. Gene set enrichment analysis corroborated that the clusters were enriched in immune-related activities. Twelve prognostic m6A-RLs were selected and used to construct the m6A-RLPS through least absolute shrinkage and selection operator Cox regression. We validated the differential expression of the 12 lncRNAs between tumor and non-cancerous samples, and the expression levels of four m6A-RLs were further validated using Gene Expression Omnibus data and Lnc2Cancer 3.0 database. The m6A-RLPS was verified to be an independent and robust predictor of ccRCC prognosis using univariable and multivariable Cox regression analyses. A nomogram based on age, tumor grade, clinical stage, and m6A-RLPS was generated and showed high accuracy and reliability at predicting the OS of patients with ccRCC. The prognostic signature was found to be strongly correlated to tumor-infiltrating immune cells and immune checkpoint expression. In conclusion, we established a novel m6A-RLPS with a favorable prognostic value for patients with ccRCC. The 12 m6A-RLs included in the signature may provide new insights into the tumorigenesis and allow the prediction of the treatment response of ccRCC.

Genomic Instability Promotes the Progression of Clear Cell Renal Cell Carcinoma Through Influencing the Immune Microenvironment

Background: Long non-coding RNAs (lncRNAs) are now under discussion as novel promising biomarkers for clear cell renal cell carcinoma (ccRCC). However, the role of genomic instability-associated lncRNA signatures in tumors has not been thoroughly uncovered. The purpose of our study is to probe the role of genomic instability-derived lncRNA signature (GILncSig) and to further investigate the mechanism of genomic instability-mediated ccRCC progression. Methods: The transcriptome data and somatic mutation profiles of ccRCC as well as clinical characteristics used in this study were obtained from The Cancer Genome Atlas database and Gene Expression Omnibus database. Lasso regression analysis was performed to construct the GILncSig. Gene set enrichment analysis (GSEA) was performed to elucidate the biological functions and relative pathways. CIBERSORT and EPIC algorithm were applied to calculate the proportion of immune cells in ccRCC. ESTIMATE algorithm was utilized to compute the immune microenvironment scores. Results: In total, 148 novel genomic instability-derived lncRNAs in ccRCC were identified. Immediately, on the basis of univariate cox analysis and lasso analysis, a GILncSig was appraised, through which the patients were allocated into High-Risk and Low-Risk groups with significantly different characteristics and prognoses. In addition, we confirmed that the somatic mutation count, tumor mutation burden, and the expression of UBQLN4, which were ascertainably associated with genomic instability, were significantly correlated with the GILncSig, indicating its reliability as a measurement of the genomic instability. Furthermore, the efficiency of GILncSig in prognostic aspects was better than the single mutation gene in ccRCC. In addition, MNX1-AS1 was defined to be a potential biomarker characterized by strong correlation with clinical features. Moreover, GSEA results indicated that the IL6/JAK/STAT3/SIGNALING pathway could be considered as a potential mechanism of genomic instability to influence tumor progression. Besides, the immune microenvironment showed significant differences between the GS-like group and the GU-like group, which was specifically manifested as high expression of CTLA4, GITR, TNFSF14, and regulatory T cells (Tregs) as well as low expression of endothelial cells (ECs) in the GU-like group. Finally, the prognostic value and clinical relevance of GILncSig were verified in GEO datasets and other urinary tumors in TCGA dataset. Conclusion: In conclusion, our study provided a new perspective for the role of lncRNAs in genomic instability and revealed that genomic instability may mediate tumor progression by affecting immunity. Besides, MNX1-AS1 played critical roles in promoting the progression of ccRCC, which may be a potential therapeutic target. What is more, the immune atlas of genomic instability was characterized by high expression of CTLA4, GITR, TNFSF14, and Tregs, and low expression of ECs.

Nucleolar protein NOP2 could serve as a potential prognostic predictor for clear cell renal cell carcinoma

As an indispensable part for cancer precision medicine, biomarkers and signatures for predicting cancer prognosis and therapeutic benefits were urgently required. The purpose of this study was to investigate the prognostic roles of NOP2 in renal clear cell carcinoma (ccRCC) for overall survival (OS) and its relationships with immunity. NOP2-related gene expression matrix associated with clinical information was obtained from the Cancer Genome Atlas (TCGA) ccRCC dataset and NOP2-related pathways were identified by gene set enrichment analysis (GSEA). Associations among the NOP2 expression and MSI, TMB, TNB, and immunity were also explored. Both the NOP2 mRNA and protein/phosphoprotein had a higher expression in ccRCC tumor tissues than in normal kidney tissues (both P < 0.001) and elevated NOP2 expression was associated with poor OS (P < 0.001). Logistic regression analysis revealed the NOP2 expression was significantly linked to stage, age, grade, N stage, T stage, and M stage (all P < 0.05). Univariate/multivariate Cox hazard regression analysis results indicated that NOP2 was an independent prognostic factor for OS in ccRCC and GSEA revealed five NOP2-related signaling pathways. Nomogram based on NOP2 and eight clinical characteristic parameters (grade, age, stage, gender, T stage, race, M stage, N stage) was constructed and carefully evaluated. Furthermore, NOP2 gene expression was also found to be significantly related to MSI, TMB, and immunity. Our findings revealed that NOP2 might be a potential prognostic factor for OS in ccRCC and it was significantly associated with immunity, MSI, and TMB.

Identification of a Four Hypoxia-Associated Long Non-Coding RNA Signature and Establishment of a Nomogram Predicting Prognosis of Clear Cell Renal Cell Carcinoma

To identify novel hypoxia-associated long non-coding RNAs (lncRNAs) as potential biomarkers, we developed a risk stratification signature and constructed a prognosis prediction nomogram of clear cell renal cell carcinoma (ccRCC). Hypoxia-related lncRNAs were identified through Pearson correlation analysis between the expression profiles of hypoxia-related differentially expressed genes and lncRNAs from The Cancer Genome Atlas Kidney Renal Clear Cell Carcinoma (TCGA-KIRC) dataset. Then, a signature of four key lncRNAs (COMETT, EMX2OS, AC026462.3, and HAGLR) was developed. The four lncRNAs were downregulated in high-grade, advanced stage, and high-risk ccRCC. The signature had an independent and long-standing prognosis prediction ability up to a 10-year follow-up. Notably, the risk score was significantly positively correlated with the infiltration abundances of six immune cells from the Tumor IMmune Estimation Resource (TIMER). The gene set enrichment analysis (GSEA) also suggested that the signature was involved in metabolism and tumorigenesis, which were closely related to the hypoxic tumor microenvironment. Ultimately, a nomogram of signature, age, stage, and grade, was built to predict the individual long-term survival possibility. Finally, the expressions of four lncRNAs were validated by quantitative real-time PCR (qRT-PCR). Our study identified a four-lncRNA signature and established a prognostic nomogram that reliably predicts survival in ccRCC. The findings may be beneficial to therapeutic customization and medical decision-making.

Broadening Drug Design and Targets to Tumor Microenvironment? Cancer-Associated Fibroblast Marker Expression in Cancers and Relevance for Survival Outcomes

Solid tumors have complex biology and structure comprising cancer cells, stromal cells, and the extracellular matrix. While most therapeutics target the cancer cells, recent data suggest that cancer cell behavior and response to treatment are markedly influenced by the tumor microenvironment (TME). In particular, the cancer-associated fibroblasts (CAFs) are the most abundant stromal cells, and play a significant contextual role in shaping tumor initiation, progression, and metastasis. CAFs have therefore emerged as part of the next-generation cancer drug design and discovery innovation strategy. We report here new findings on differential expression and prognostic significance of CAF markers in several cancers. We utilized two publicly available resources: The Cancer Genomic Atlas and Gene Expression Profiling Interactive Analysis. We examined the expression of CAF markers, ACTA2, S100A4, platelet-derived growth factor receptor-beta [PDGFR-β], CD10, and fibroblast activation protein-alpha (FAP-α), in tumor tissues versus the adjacent normal tissues. We found that CAF markers were differentially expressed in various different tumors such as colon, breast, and esophageal cancers and melanoma. No CAF marker is expressed in the same pattern in all cancers, however. Importantly, we report that patients with colon adenocarcinoma and esophageal carcinoma expressing high FAP-α and CD10, respectively, had significantly shorter overall survival, compared with those with low levels of these CAF markers (p < 0.05). We call for continued research on TME biology and clinical evaluation of the CAF markers ACTA2, S100A4, PDGFR-β, CD10, and FAP-α in relation to prognosis of solid cancers in large population samples. An effective cancer drug design and discovery roadmap in the 21st century ought to be broadly framed, and include molecular targets informed by both cancer cell and TME variations.

Process validation, current good manufacturing practice production, dosimetry, and toxicity studies of the carbonic anhydrase IX imaging agent [111 In]In-XYIMSR-01 for phase I regulatory approval

[111 In]In-XYIMSR-01 is a promising single-photon emission computed tomography (SPECT) imaging agent for identification of tumors that overexpress carbonic anhydrase IX. To translate [111 In]In-XYIMSR-01 to phase I trials, we performed animal toxicity and dosimetry studies, determined the maximum dose for human use, and completed the chemistry, manufacturing, and controls component of a standard regulatory application. The production process, quality control testing, stability studies, and specifications for sterile drug product release were based on United States Pharmacopeia chapters <823> and <825>, FDA 21 CFR Part 212. Toxicity was evaluated by using nonradioactive [113/115 In]In-XYIMSR-01 according to 21 CFR Part 58 guidelines. Organ Level INternal Dose Assessment/EXponential Modeling (OLINDA/EXM) was used to calculate the maximum single dose for human studies. Three process validation runs at starting radioactivities of ~800 MBq were completed with a minimum concentration of 407 MBq/ml and radiochemical purity of ≥99% at the end of synthesis. A single intravenous dose of 55 μg/ml of [113/115 In]In-XYIMSR-01 was well tolerated in male and female Sprague-Dawley rats. The calculated maximum single dose for human injection from dosimetry studies was 390.35 MBq of [111 In]In-XYIMSR-01. We have completed toxicity and dosimetry studies as well as validated a manufacturing process to test [111 In]In-XYIMSR-01 in a phase I clinical trial.

KLHDC7B-DT aggravates pancreatic ductal adenocarcinoma development via inducing cross-talk between cancer cells and macrophages

Tumor microenvironment (TME) exerts key roles in pancreatic ductal adenocarcinoma (PDAC) development. However, the factors regulating the cross-talk between PDAC cells and TME are largely unknown. In the present study, we identified a long noncoding RNA (lncRNA) KLHDC7B divergent transcript (KLHDC7B-DT), which was up-regulated in PDAC and correlated with poor survival of PDAC patients. Functional assays demonstrated that KLHDC7B-DT enhanced PDAC cell proliferation, migration, and invasion. Mechanistically, KLHDC7B-DT was found to directly bind IL-6 promoter, induce open chromatin structure at IL-6 promoter region, activate IL-6 transcription, and up-regulate IL-6 expression and secretion. The expression of KLHDC7B-DT was positively correlated with IL-6 in PDAC tissues. Via inducing IL-6 secretion, KLHDC7B-DT activated STAT3 signaling in PDAC cells in an autocrine manner. Furthermore, KLHDC7B-DT also activated STAT3 signaling in macrophages in a paracrine manner, which induced macrophage M2 polarization. KLHDC7B-DT overexpressed PDAC cells-primed macrophages promoted PDAC cell proliferation, migration, and invasion. Blocking IL-6/STAT3 signaling reversed the effects of KLHDC7B-DT on macrophage M2 polarization and PDAC cell proliferation, migration, and invasion. In conclusion, KLHDC7B-DT enhanced malignant behaviors of PDAC cells via IL-6-induced macrophage M2 polarization and IL-6-activated STAT3 signaling in PDAC cells. The cross-talk between PDAC cells and macrophages induced by KLHDC7B-DT represents potential therapeutic target for PDAC.

A risk score model with five long non-coding RNAs for predicting prognosis in gastric cancer: an integrated analysis combining TCGA and GEO datasets

Background

Gastric cancer (GC) is one of the most common carcinomas of the digestive tract, and the prognosis for these patients may be poor. There is evidence that some long non-coding RNAs(lncRNAs) can predict the prognosis of patients with GC. However, few lncRNA signatures have been used to predict prognosis. Herein, we aimed to construct a risk score model based on the expression of five lncRNAs to predict the prognosis of patients with GC and provide new potential therapeutic targets.

Methods

We performed differentially expressed and survival analyses to identify differentially expressed survival-ralated lncRNAs by using GC patient expression profile data from The Cancer Genome Atlas (TCGA) database. We then established a formula including five lncRNAs to predict the prognosis of patients with GC. In addition, to verify the prognostic value of this risk score model, two independent Gene Expression Omnibus (GEO) datasets, GSE62254 (N = 300) and GSE15459 (N = 200), were employed as validation groups.

Results

Based on the characteristics of five lncRNAs, patients with GC were divided into high or low risk subgroups. The prognostic value of the risk score model with five lncRNAs was confirmed in both TCGA and the two independent GEO datasets. Furthermore, stratification analysis results showed that this model had an independent prognostic value in patients with stage II-IV GC. We constructed a nomogram model combining clinical factors and the five lncRNAs to increase the accuracy of prognostic prediction. Enrichment analysis based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) suggested that the five lncRNAs are associated with multiple cancer occurrence and progression-related pathways.

Conclusion

The risk score model including five lncRNAs can predict the prognosis of patients with GC, especially those with stage II-IV, and may provide potential therapeutic targets in future.

Effects of LncRNA HCP5/miR-214-3p/MAPK1 Molecular Network on Renal Cell Carcinoma Cells

Background

Recent researches have shown that long non-coding RNA (LncRNA) is often disordered and acts in many carcinomas. Clear cell renal cell carcinoma (ccRCC) is the main reason for carcinoma-related deaths, which are mainly caused by the metastasis. HCP5 is a newly discovered LcnRNA. Early studies have found that HCP5 acts in neoplasm metastasis, but the mechanism of HCP5 in ccRCC is still unclear.

Methods

The expression of HCP5 in human renal cell carcinoma (RCC) was detected by real-time quantitative PCR. The biological effect of LncRNAs in proliferation, migration, invasion and metastasis of RCC cells was explored by gain-of-function and loss-of-function tests. The molecular mechanism of LncRNAs was explored by RNA immunoprecipitation and Western blot.

Results

qRT-PCR revealed that HCP5 was enhanced in neoplasm tissues of ccRCC patients and correlated with the metastatic characteristics of RCC. Over-expression of HCP5 promoted the proliferation, migration and invasion of renal carcinoma cells. The deletion of HCP5 inhibited the proliferation, migration and invasion of RCC in vitro and the metastasis of RCC in vivo. Mechanically, HCP5 inhibited the growth and metastasis of ccRCC cells by regulating miR-214-3p/MAPK1 axis.

Conclusion

HCP5, as a key LncRNA, can promote ccRCC metastasis by regulating miR-214-3p/MAPK1 axis and may be a biomarker and be helpful for judging the prognosis of ccRCC.

Animal model and bioinformatics analyses suggest the TIMP1/MMP9 axis as a potential biomarker in oral squamous cell carcinoma

Oral squamous cell carcinoma (OSCC) is a common malignant tumor of the head and neck. However, the molecular mechanism underlying its development and progression is yet unclear. Genes that are differentially expressed, that is, differentially expressed genes (DEGs), between normal and diseased tissues are believed to be involved in disease development and progression. To identify the DEGs in OSCC and explore their role in occurrence and progression, we established a Chinese hamster OSCC model, determined the DEG, screened the identified DEGs, and performed Gene Ontology (GO) and KEGG enrichment analyses. A protein-protein interaction (PPI) network was generated to screen potential candidate genes. We then analyzed the expression, tumor stage and prognosis of candidate genes using the Gene Expression Profiling Interactive Analysis (GEPIA) database. Finally, we verified the candidate DEGs by quantitative real-time PCR and Gene Expression Omnibus analysis. The results showed 194 significantly DEGs, 140 enriched GO terms, and 8 KEGG pathways, which suggested that OSCC was closely related to the immune system, cell migration, and extracellular matrix. GEPIA and PPI network analysis revealed that SPP1, TNC, and ACTA1 were significantly related to tumor staging; SPP1, tissue inhibitors of matrix metallopeptidases (MMPs) 1 (TIMP1), and ACTA1 were closely related to prognosis. The scores for the top five highest degree genes were close, and the TIMP1/MMP9 axis appeared to be at the center of the PPI network, indicating that expression changes in the TIMP1/MMP9 axis and related genes may be involved in tumor invasion and metastasis. These findings provide novel insights into the mechanism of oral cancer.

A novel prognostic model based on immunogenomics for clear cell renal cell carcinoma

BACKGROUND

Immune cell infiltration into tumor tissue is closely related to the clinical outcomes of patients with clear cell renal cell carcinoma (ccRCC). This study aimed to screen out potential immune genes associated with ccRCC, analyze their relationships with clinical outcomes, and construct a signature to predict ccRCC.

METHODS

The transcriptome RNA-sequencing data in 539 ccRCC and 72 adjacent normal tissues were obtained from TCGA database. Biomedical computational algorithms were conducted to identify immune-related differential expressed genes (IRDGs) and enriched pathways. Then, LASSO Cox and multivariate Cox analyses were performed to screen out genes that were then used to construct the prognostic model.

RESULTS

A total of 116 down-regulated and 565 up-regulated IRDGs were identified. Pathway enrichment analysis suggested that IRDGs was mainly enriched in the pathway of "cytokines and cytokine receptors". The entire data of ccRCC were randomly divided into the training set and the test set with a ratio of 1:1. A 4-gene signature was then constructed using LASSO Cox analysis and multivariate Cox analysis in the training set. This prognostic signature could stratify patients into high- and low-risk groups successfully, and serve as an independent predictor when adjusted with clinical factors by univariate and multivariate Cox regression analysis. These results were verified in the test set and the entire set. Besides, the abundance of CD4 + T cells and dendritic cells increased in the high-risk group. Finally, we built a nomogram incorporating risk score and clinical factors to predict the overall survival of ccRCC patients.

CONCLUSIONS

These findings may contribute to the research of ccRCC in immunization part.

Radiomics-based machine learning model to predict risk of death within 5-years in clear cell renal cell carcinoma patients

PURPOSE

The aim of this study was to develop radiomics-based machine learning models based on extracted radiomic features and clinical information to predict the risk of death within 5 years for prognosis of clear cell renal cell carcinoma (ccRCC) patients.

METHODS

According to image quality and clinical data availability, we eventually selected 70 ccRCC patients that underwent CT scans. Manual volume-of-interest (VOI) segmentation of each image was performed by an experienced radiologist using the 3D slicer software package. Prior to feature extraction, image pre-processing was performed on CT images to extract different image features, including wavelet, Laplacian of Gaussian, and resampling of the intensity values to 32, 64 and 128 bin levels. Overall, 2544 3D radiomics features were extracted from each VOI for each patient. Minimum Redundancy Maximum Relevance (MRMR) algorithm was used as feature selector. Four classification algorithms were used, including Generalized Linear Model (GLM), Support Vector Machine (SVM), K-nearest Neighbor (KNN) and XGBoost. We used the Bootstrap resampling method to create validation sets. Area under the receiver operating characteristic (ROC) curve (AUROC), accuracy, sensitivity, and specificity were used to assess the performance of the classification models.

RESULTS

The best single performance among 8 different models was achieved by the XGBoost model using a combination of radiomic features and clinical information (AUROC, accuracy, sensitivity, and specificity with 95% confidence interval were 0.95-0.98, 0.93-0.98, 0.93-0.96 and ~1.0, respectively).

CONCLUSIONS

We developed a robust radiomics-based classifier that is capable of accurately predicting overall survival of RCC patients for prognosis of ccRCC patients. This signature may help identifying high-risk patients who require additional treatment and follow up regimens.

Long non-coding RNA LINC01234 regulates proliferation, migration and invasion via HIF-2α pathways in clear cell renal cell carcinoma cells

Long non-coding RNAs (lncRNAs) have been proved to have an important role in different malignancies including clear cell renal cell carcinoma (ccRCC). However, their role in disease progression is still not clear. The objective of the study was to identify lncRNA-based prognostic biomarkers and further to investigate the role of one lncRNA LINC01234 in progression of ccRCC cells. We found that six adverse prognostic lncRNA biomarkers including LINC01234 were identified in ccRCC patients by bioinformatic analysis using The Cancer Genome Atlas database. LINC01234 knockdown impaired cell proliferation, migration and invasion in vitro as compared to negative control. Furthermore, the epithelial-mesenchymal transition was inhibited after LINC01234 knockdown. Additionally, LINC01234 knockdown impaired hypoxia-inducible factor-2a (HIF-2α) pathways, including a suppression of the expression of HIF-2α, vascular endothelial growth factor A, epidermal growth factor receptor, c-Myc, Cyclin D1 and MET. Together, these datas showed that LINC01234 was likely to regulate the progression of ccRCC by HIF-2α pathways, and LINC01234 was both a promising prognostic biomarker and a potential therapeutic target for ccRCC.

Prognostic Value and Potential Biological Functions of CLDN8 in Patients with Clear Cell Renal Cell Carcinoma

Purpose

Clear cell renal cell carcinoma (ccRCC) is among the most common malignant tumors worldwide, with a high incidence rate and poor prognosis. Currently, there are no biomarkers that can accurately guide prognostic evaluation and therapeutic strategy for ccRCC. The prognostic value and potential biological function of claudin-8 (CLDN8), a critical component of tight junctions in ccRCC, remain unclear.

Methods

Sequencing data were obtained from The Cancer Genome Atlas, International Cancer Genome Consortium, and Gene Expression Omnibus databases. R packages were used to explore CLDN8 mRNA expression levels and analyze differentially expressed genes. Results were validated in clinical specimens and cell lines, and bioinformatics analyses were conducted to explore the potential biological functions of CLDN8. Finally, functional analyses were carried out using 786–O ccRCC cell line.

Results

Both CLDN8 mRNA and protein expression levels were significantly lower in ccRCC compared with the normal control tissues. Kaplan–Meier analyses showed that low CLDN8 expression levels were associated with the poor overall survival, while univariate and multivariate Cox regression indicated that CLDN8 could serve as an independent prognostic factor in patient with ccRCC. Bioinformatic and Western blot analyses showed that CLDN8 suppressed proliferation, migration, and invasion of 786–O ccRCC cells through the epithelial–mesenchymal transition and AKT pathways.

Conclusion

CLDN8 could serve as an independent prognostic factor in ccRCC, in which it suppresses 786–O proliferation, migration, and invasion through EMT and AKT pathways.

Influence of gene expression on survival of clear cell renal cell carcinoma

Approximately 10%‐20% of patients with clinically localized clear cell renal cell carcinoma (ccRCC) at time of surgery will subsequently experience metastatic progression. Although considerable progression was seen in the systemic treatment of metastatic ccRCC in last 20 years, once ccRCC spreads beyond the confines of the kidney, 5‐year survival is less than 10%. Therefore, significant clinical advances are urgently needed to improve overall survival and patient care to manage the growing number of patients with localized ccRCC. We comprehensively evaluated expression of 388 candidate genes related with survival of ccRCC by using TCGA RNAseq (n = 515), Total Cancer Care (TCC) expression array data (n = 298), and a well characterized Moffitt RCC cohort (n = 248). We initially evaluated all 388 genes for association with overall survival using TCGA and TCC data. Eighty‐one genes were selected for further analysis and tested on Moffitt RCC cohort using NanoString expression analysis. Expression of nine genes (AURKA, AURKB, BIRC5, CCNE1, MK167, MMP9, PLOD2, SAA1, and TOP2A) was validated as being associated with poor survival. Survival prognostic models showed that expression of the nine genes and clinical factors predicted the survival in ccRCC patients with AUC value: 0.776, 0.821 and 0.873 for TCGA, TCC and Moffitt data set, respectively. Some of these genes have not been previously implicated in ccRCC survival and thus potentially offer insight into novel therapeutic targets. Future studies are warranted to validate these identified genes, determine their biological mechanisms and evaluate their therapeutic potential in preclinical studies.

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

Clear cell renal cell carcinoma (ccRCC) is the most common type of renal cell carcinoma (RCC). Despite the existing extensive research, the molecular and pathogenic mechanisms of ccRCC are elusive. We aimed to identify the immune-related lncRNA signature and molecular subtypes associated with ccRCC. By integrating 4 microarray datasets from Gene Expression Omnibus database, we identified 49 immune-related genes. The corresponding immune-related lncRNAs were further identified in the TCGA dataset. 12-lncRNAs prognostic and independent signature was identified through survival analysis and survival difference between risk groups was further identified based on the risk score. Besides, we identified 3 molecular subtypes and survival analysis result showed that cluster 2 has a better survival outcome. Further, ssGSEA enrichment analysis for the immune-associated gene sets revealed that cluster 1 corresponded to a high immune infiltration level. While cluster 2 and cluster 3 corresponded to low and medium immune infiltration level, respectively. In addition, we validated the 12-lncRNA prognostic signature and molecular subtypes in an external validation dataset from the ICGC database. In summary, we identified a 12-lncRNA prognostic signature which may provide new insights into the molecular mechanisms of ccRCC and the molecular subtypes provided a theoretical basis for personalized treatment by clinicians.

Clinical significance of CCNE2 protein and mRNA expression in thyroid cancer tissues

PURPOSE

Thyroid carcinoma (TC) is the most common endocrinal malignancy worldwide. Cyclin E2 (CCNE2), a member of the cyclin family, acts as a regulatory subunit of cyclin-dependent kinases (CDKs). It controls the transition of quiescent cells into the cell cycle, regulates the G1/S transition, promotes DNA replication, and activates CDK2. This study explored the role and potential molecular mechanisms of CCNE2 expression in TC tissues.

MATERIAL/METHODS

Immunohistochemistry was used to evaluate the CCNE2 protein expression levels in TC. High-throughput data on CCNE2 in TC were obtained from RNA sequencing (RNA-seq), microarray, and literature data. The CCNE2 expression levels in TC were comprehensively assessed through an integrated analysis. Analyses of Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and protein-protein interaction (PPIs) data facilitated the investigation of the relative molecular mechanisms of CCNE2 in TC.

RESULTS

The immunohistochemical experiment showed a significant increase in the expression of CCNE2 in the TC tissues. For 505 TC and 59 non-cancerous samples from RNA-seq data, the area under the curve (AUC) was 0.8016 (95% confidence interval [CI] 0.742-0.8612; p<0.001). With another 14 microarrays, the pool standard mean difference [SMD] was 1.01 (95% CI [0.82-1.19]). The pooled SMD of CCNE2 was 1.12 (95% CI [0.60-1.64]), and the AUC was 0.87 (95% CI [0.84-0.90]) for 1157 TC samples and 366 non-cancerous thyroid samples from all possible sources. Nine hub genes were upregulated in TC.

CONCLUSIONS

A high expression of CCNE2 may lead to carcinogenesis and the development of TC.

Identification of a three-long noncoding RNA prognostic model involved competitive endogenous RNA in kidney renal clear cell carcinoma

Background

Long noncoding RNA (lncRNA) is generally identified as competing endogenous RNA (ceRNA) that plays a vital role in the pathogenesis of kidney renal clear cell carcinoma (KIRC), the most common subtype of renal cell carcinoma with poor prognosis and unclear pathogenesis. This study established a novel ceRNA network and thus identified a three-lncRNA prognostic model in KIRC patients.

Methods

Differentially expressed genes (DEGs) were screened out from The Cancer Genome Atlas (TCGA) database. The lncATLAS was applied to determine the differentially expressed lncRNAs (DElncRNAs) of the cytoplasm. The miRcode, miRDB, miRTarBase, and TargetScan databases were utilized to predict the interactions of DElncRNAs, DEmiRNAs, and DEmRNAs. Cytoscape was used to construct the ceRNA network. Then, a lncRNA prognostic model (LPM) was constructed based on ceRNA-related lncRNA that was significantly related to overall survival (OS), and its predictive ability was evaluated. Moreover, an LPM-based nomogram model was constructed. The significantly different expression of genes in the LPM was validated in an independent clinical cohort (N = 21) by quantitative RT-PCR.

Results

A novel ceRNA regulatory network, including 73 lncRNAs, 8 miRNAs, and 21 mRNAs was constructed. Functional enrichment analysis indicated that integral components of membrane and PI3K-Akt signaling pathway represented the most significant GO terms and pathway, respectively. The LPM established based on three lncRNAs (MIAT, LINC00460, and LINC00443) of great prognostic value from the ceRNA network was proven to be independent of conventional clinical parameters to differentiate patients with low or high risk of poor survival, with the AUC of 1-, 5- and 10-year OS were 0.723, 0.714 and 0.826 respectively. Furthermore, the nomogram showed a better predictive value in KIRC patients than individual prognostic parameters. The expression of MIAT and LINC00460 was significantly upregulated in the KIRC samples, while the expression of LINC00443 was significantly downregulated compared with the adjacent normal samples in the clinical cohort, TCGA, and GTEx.

Conclusion

This LPM based on three-lncRNA could serve as an independent prognostic factor with a tremendous predictive ability for KIRC patients, and the identified novel ceRNA network may provide insight into the prognostic biomarkers and therapeutic targets of KIRC.

Identification of an Autophagy-Related Prognostic Signature for Clear Cell Renal Cell Carcinoma

Abnormal autophagy is closely related to the development of cancer. Many studies have demonstrated that autophagy plays an important role in biological function in clear cell renal cell carcinoma (ccRCC). This study aimed to construct a prognostic signature for ccRCC based on autophagy-related genes (ARGs) to predict the prognosis of ccRCC. Differentially expressed ARGs were obtained from ccRCC RNA-seq data in The Cancer Genome Atlas (TCGA) database. ARGs were enriched by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). The prognostic ARGs used to construct the risk score models for overall survival (OS) and disease-free survival (DFS) were identified by Cox regression analyses. According to the median value of the risk score, patients were divided into a high-risk group and a low-risk group. The OS and DFS were analyzed by the Kaplan-Meier method. The predictive accuracy was determined by a receiver operating characteristic (ROC) curve analysis. Additionally, we performed stratification analyses based on different clinical variables and evaluated the correlation between the risk score and the clinical variables. The differentially expressed ARGs were mainly enriched in the platinum drug resistance pathway. The prognostic signatures based on 11 ARGs for OS and 5 ARGs for DFS were constructed and showed that the survive time was significantly shorter in the high-risk group than in the low-risk group (P < 0.001). The ROC curve for OS exhibited good predictive accuracy, with an area under the curve value of 0.738. In the stratification analyses, the OS time of the high-risk group was shorter than that of the low-risk group stratified by different clinical variables. In conclusion, an autophagy-related signature for OS we constructed can independently predict the prognosis of ccRCC patient, and provide a deep understanding of the potential biological mechanisms of autophagy in ccRCC.

Comprehensive analysis of competitive endogenous RNAs network reveals potential prognostic lncRNAs in gastric cancer

Long non-coding RNAs (lncRNAs) are key regulators of a range of human diseases, including various cancers, with multiple previous studies having explored lncRNA dysregulation in the context of gastric cancer (GC). The present study sought to expand upon these previous results by downloading lncRNA, mRNA, and microRNA (miRNA) expression profiles derived from 180 GC tissues and 24 normal control tissues within the Cancer Genome Atlas (TCGA) database. These datasets were then interrogated to identify GC-related differentially expressed (DE) RNAs (|fold change| ≥ 2, FDR< 0.01), leading to the identification of 1946 DE lncRNAs, 123 DE miRNAs, and 3159 DE mRNAs. These results were then used to generate a putative GC-related competitive endogenous RNA (ceRNA) network composed of 131 lncRNAs, 9 miRNAs, and 78 mRNAs. Subsequent survival analyses based upon this network revealed 17 of these lncRNAs to be significantly associated with GC patient survival (P < 0.05). Further multivariable Cox regression and lasso analyses allowed for the construction of an 8-lncRNA risk score that was able to effectively predict GC patient survival with good discriminative ability. The Kaplan-Meier Plotter database further confirmed that network hub genes that were related to these 8 lncRNAs were associated with GC patient prognosis (P < 0.05). As the ceRNA network in the present study was constructed with a focus on both disease stage and differential gene expression, it represents a key resource that will offer valuable insights into the mechanistic roles of ceRNA pathways in GC development and progression.

Prognostic Inflammasome-Related Signature Construction in Kidney Renal Clear Cell Carcinoma Based on a Pan-Cancer Landscape

OBJECTIVE

To investigate the expression patterns and prognostic characteristics of inflammasome-related genes (IRGs) across cancer types and develop a robust biomarker for the prognosis of KIRC.

METHODS

The differentially expressed IRGs and prognostic genes among 10 cancers were analyzed based on The Cancer Genome Atlas (TCGA) dataset. Subsequently, an IRGs risk signature was developed in KIRC. Its prognostic accuracy was evaluated by receiver operating characteristic (ROC) analysis. The independent predictive capacity was identified by stratification survival and multivariate Cox analyses. The gene ontology (GO) analysis and principal component analysis (PCA) were performed to explore biological functions of the IRGs signature in KIRC.

RESULTS

The expression patterns and prognostic association of IRGs varied from different cancers, while KIRC showed the most abundant survival-related dysregulated IRGs. The IRG signature for KIRC was able to independently predict survival, and the signature genes were mainly involved inimmune-related processes.

CONCLUSIONS

The pan-cancer analysis provided a comprehensive landscape of IRGs across cancer types and identified a strong association between IRGs and the prognosis of KIRC. Further IRGs signature represented a reliable prognostic predictor for KIRC and verified the prognostic value of inflammasomes in KIRC, contributing to our understanding of therapies targeting inflammasomes for human cancers.

Deciphering Immune-Associated Genes to Predict Survival in Clear Cell Renal Cell Cancer

Background

To elucidate the correlations between tumor microenvironment and clinical characteristics as well as prognosis in clear cell renal cell cancer (ccRCC) and investigate the immune-associated genes by a comprehensive analysis of The Cancer Genome Atlas (TCGA) database.

Methods

We collected mRNA expression profiles of 537 ccRCC samples from the TCGA database. Immune scores and stromal scores were calculated by applying the ESTIMATE algorithm. We evaluated the correlation between immune/stromal scores and clinical characteristics as well as prognosis. The differentially expressed genes (DEGs) were screened between high immune/stromal score and low immune/stromal score groups by the cutoff of |log (fold change)| > 1, P value <0.05 by using package "limma" in R. Functional enrichment analysis was performed by DAVID, and the protein-protein interaction network of intersected DEGs between stromal score and immune score groups was conducted using the STRING database. The Kaplan-Meier method was used to explore DEGs with predictive values in overall survival, and the prognostic DEGs were further validated in a Gene Expression Omnibus (GEO) dataset GSE29609.

Results

A higher immune score was associated with T3/4 (vs. T1/2, P < 0.001), N1 (vs. N0, P=0.05), M1 (vs. M0, P=0.004), G3/4 (vs. G1/2, P < 0.001), advanced AJCC stage (P < 0.001), and shorter overall survival (P=0.04). Intersected DEGs between immune and stromal score groups were 48 upregulated and 47 downregulated genes, with 43 DEGs associated with overall survival in ccRCC. After validation by a cohort of 39 ccRCC cases with detailed follow-up information from GSE29609, six immune-associated DEGs including CASP5, HSD11B1, VSIG4, HMGCS2, HSD11B2, and OGDHL were demonstrated to be predictive of prognosis in ccRCC.

Conclusions

Our study elucidated tight associations between immune score and clinical characteristics as well as prognosis in ccRCC. Moreover, six DEGs were explored and validated to exert predictive values in overall survival of ccRCC.

Development and validation of a metastasis-associated prognostic signature based on single-cell RNA-seq in clear cell renal cell carcinoma

Single-cell RNA sequencing (scRNA-seq) was recently adopted for deciphering intratumoral heterogeneity across cell sub-populations, including clear cell renal cell carcinoma (ccRCC). Here, we characterized the single-cell expression profiling of 121 cell samples and found 44 metastasis-associated marker genes. Accordingly, we trained and validated 17 pivotal metastasis-associated genes (MAGs) in 626 patients incorporating internal and external cohorts to evaluate the model for predicting overall survival (OS) and progression-free survival (PFS). Correlation analysis revealed that the MAGs correlated significantly with several risk clinical characteristics. Moreover, we conducted Cox regression analysis integrating these independent clinical variables into a MAGs nomogram with superior accuracy in predicting progression events. We further revealed the differential landscape of somatic tumor mutation burden (TMB) between two nomogram-score groups and observed that TMB was also a prognostic biomarker; patients with high MAGs-nomogram scores suffered from a higher TMB, potentially leading to worse prognosis. Last, higher MAGs-nomogram scores correlated with the upregulation of oxidative phosphorylation, the Wnt signaling pathway, and MAPK signaling crosstalk in ccRCC. Overall, we constructed the robust MAGs through scRNA-seq and validated the model in a large patient population, which was valuable for prognostic stratification and providing potential targets against metastatic ccRCC.