Identification of an m6A-Related lncRNA Signature for Predicting the Prognosis in Patients With Kidney Renal Clear Cell Carcinoma

Correlation study of LINC02609 and SNHG17 as prognostic biomarkers of kidney renal clear cell carcinoma and therapeutic sensitivity based on public data and In Vitro analysis

Background

Cuprotosis, a newly identified form of regulated cell death, has emerged as a potential therapeutic target for cancers. Kidney renal clear cell carcinoma (KIRC) is frequently metastatic at diagnosis, resulting in poor prognosis. This study aimed to identify prognostic biomarkers and construct a risk model to improve survival prediction and guide therapeutic strategies for KIRC patients.

Methods

Differential expression analysis, Cox regression, and risk modeling were performed using transcriptomic and clinical data. The response to immunotherapy and the sensitivity to chemotherapy drugs were analyzed through the Tumor Immune Dysfunction and Exclusion (TIDE) database and the Genomics of Drug Sensitivity in Cancer2 (GDSC2) database. Functional validation of LINC02609 was conducted in renal carcinoma A498 cells using siRNA-mediated knockdown.

Results

LINC02609 and SNHG17 were significantly upregulated in KIRC tissues and independently associated with poor overall survival. The risk model constructed using those two candidate biomarkers (LINC02609 and SNHG17) exhibited high predictive accuracy as measured by the value of area under the curve (AUC). Immune status analysis showed that high- and low-risk KIRC patients exhibited abnormalities immune landscapes. TIDE analysis suggested that the risk model was significantly correlated with multiple immunotherapy-related signatures. RNA-sequencing (RNA-seq) analysis indicated that inhibition of LINC02609 would lead to abnormal activation of the mitogen-activated protein kinases (MAPK) signaling pathway. In vitro experiments confirmed that LINC02609 knockout inhibits the proliferation, migration, and invasion of A498 cells by suppressing the MAPK signaling pathway.

Conclusion

The candidate biomarker LINC02609 regulates the progression of renal cell carcinoma through the MAPK signaling pathway. The risk model constructed using LINC02609 and SNHG17 was significantly correlated with multiple immunotherapy-related signatures, suggesting that it might be used for the determination of immunotherapy options in KIRC.

Identification of fatty acid metabolism-related genes in the tumor microenvironment of breast cancer by a development and validation of prognostic index signature

BACKGROUND

Breast cancer (BRCA) is a malignancy originating in the breast cells, characterized by a poor overall survival rate. Post-resection, chemotherapy is commonly recommended as a primary therapeutic approach; however, its efficacy remains limited. Recent advancements in lipidomics and metabolomics have provided new insights into the intricate landscape of fatty acid metabolism (FAM) and the fatty acid lipidome in both health and disease. A growing body of evidence suggests that dysregulations in FAM and fatty acid levels play a significant role in cancer initiation and progression. Despite these advances, the precise mechanisms through which FAM mediates the anti-cancer effects of lobaplatin in BRCA remain poorly understood and warrant further investigation.

METHODS

GEO and TCGA data were classified into two types. We aimed to show how FAMGs influence immune function, immune checkpoints, and m6a in BRCA. A co-expression analysis discovered that gene expression is strongly connected to pyroptosis. The TCGA gathered information about mRNAsi, gene mutations, CNV, and clinical features.

RESULTS

In the low-risk group, overall survival (OS) is longer. GSEA was utilized to identify immune and tumor-related pathways. Most of the FAMG-derived prognostic signatures predominantly modulate immunological and oncogenic signaling pathways, including the Wnt, neurotrophin, chemokine, and calcium signaling cascades. Among the genes involved are CEL, WT1, and ULBP2. Expression levels varied as well. The prognostic model, CNVs, single nucleotide polymorphism (SNP), and drug sensitivity all pointed to the gene.

CONCLUSIONS

The primary objective of this study is to identify and validate BRCA-associated FAMGs that can serve as prognostic indicators and provide insights into immune system function, while also offering evidence to support the development of fatty acid metabolism-related molecularly targeted therapeutics. Consequently, FAMGs and their interactions with the immune system, as well as their role in BRCA, may emerge as promising therapeutic targets.

Artificial intelligence based on imaging data to predict rectal cancer recurrence: A meta-analysis

PURPOSE

The purpose of this study was to evaluate the diagnostic performance of artificial intelligence based on imaging data to predict rectal cancer recurrence using a meta-analysis system.

MATERIALS AND METHODS

Medline, Embase, Cochrane Library, Web of Science, and other databases were searched for all articles on artificial intelligence prediction of rectal cancer recurrence based on imaging data published publicly from the establishment of the library to December 31, 2023. The quality of the articles was assessed using Quality Assessment of Diagnostic Accuracy Studies 2 (QUADAS-2). Meta-analysis was performed by the software Revman 5.4 and Statistics data (Stata), and sensitivity analysis was used to detect potential sources of heterogeneity and test to assess the presence of publication bias. We evaluated how well imaging-based data can predict recurrence in patients with rectal cancer by analysing the pooled sensitivity, specificity, and area under the curve.

RESULTS

Ten studies were included. The pooled sensitivity, specificity, and area under the curve of imaging-based data for recurrence in patients with rectal cancer were respectively 0.84 (95 % confidence interval [CI]: 0.74-0.91), 0.87 (95 % CI: 0.82-0.91) and 0.92 (95 % CI: 0.89-0.94). Based on QUADAS-2, the quality of the article is acceptable. We found the causes of heterogeneity through meta-regression: recurrence time predesign Lasso. Based on Deeks' funnel plot, no publication bias was detected.

CONCLUSION

Artificial intelligence based on imaging data has a high predictive ability for rectal cancer recurrence.

Study of an N6-methyladenosine- and ferroptosis-related prognostic model and the mechanisms underlying the molecular network in neuroblastoma based on multiple datasets

Recent research highlights the pivotal role of N6-methyladenosine (m6A) modification and ferroptosis in the evolution of various cancers. This study aimed to establish a prognostic framework centered on genes associated with m6A and ferroptosis to enhance the accuracy of prognosis predictions for neuroblastoma (NB) patients, thereby improving targeted therapeutic strategies. Patient data, including expression profiles and clinical information from NB cases, were acquired from The Cancer Genome Atlas. Genes related to m6A modification and ferroptosis were identified, and those significant for prognosis were pinpointed using a combination of Cox regression analysis and the least absolute shrinkage and selection operator (LASSO) regression. For further validation, the study utilized external datasets GSE62564 and GSE85047. A prognostic index was computed for each NB patient, followed by analyses of immune cell infiltration and potential drug responsiveness based on the prognostic model. Additionally, enrichment analysis was conducted on the prognostic scores. These scores showed a strong association with the tumor immune environment and the efficacy of prevalent cancer therapies. Moreover, the model's prognostic score emerged as an independent predictive marker for NB. This research succeeded in creating and confirming a prognostic model rooted in m6A and ferroptosis-linked genes, promising to enrich the prognostic understanding and treatment approaches for NB.

Development and validation of a prognostic model based on m6A-related lncRNAs to predict prognosis for papillary renal cell cancer patients

To evaluate the predictive utility of N6-methyladenosine (m6A)-associated long non-coding RNAs (lncRNAs) for the prognosis and immunotherapy response in papillary renal cell carcinoma (pRCC). Transcriptomic data of pRCC samples were extracted from the TCGA database. The m6A-related lncRNAs were identified by Pearson correlation analysis. Univariate and LASSO regression analyses were used to develop a risk model. The discrimination and predictive ability were evaluated through survival analysis, ROC analysis and consensus clustering. Tumor mutation burden (TMB) and immune infiltration of the risk groups were compared. A prognostic nomogram was constructed using six m6A-related lncRNAs, and validated through calibration and decision curve analysis (DCA). The lncRNAs HCG25 and NOP14-AS1 were knocked down in a human pRCC cell line using specific siRNA constructs, and the proliferation and migration rates were assessed by the CCK-8 and transwell assays. We identified a total of 153 m6A-related lncRNAs in pRCC datasets, of which six were selected for constructing a m6A-related lncRNA pRCC prognostic model. Mutations in the SETD2 gene correlated with worse prognosis. Significant differences were observed in immune cell infiltration between the two risk groups. A clinical prognostic nomogram for pRCC was further established based on clinical variables. In vitro assays further showed that HCG25 and NOP14-AS1 regulate the proliferation and migration of pRCC cells. The results validated the discrimination ability of both the m6A-related lncRNA pRCC prognostic model and the pRCC clinical prognostic nomogram. We developed a clinical prognostic nomogram for pRCC using pRCC prognostic-associated m6A-related lncRNAs, which can be utilized for predicting the prognosis and immune landscape of pRCC patients.

Disulfidptosis-associated LncRNA signature predicts prognosis and immune response in kidney renal clear cell carcinoma

BACKGROUND

Kidney renal clear cell carcinoma (KIRC) represents a significant proportion of renal cell carcinomas and is characterized by high aggressiveness and poor prognosis despite advancements in immunotherapy. Disulfidptosis, a novel cell death pathway, has emerged as a critical mechanism in various cellular processes, including cancer. This study leverages machine learning to identify disulfidptosis-related long noncoding RNAs (DRlncRNAs) as potential prognostic biomarkers in KIRC, offering new insights into tumor pathogenesis and treatment avenues.

RESULTS

Our analysis of data from The Cancer Genome Atlas (TCGA) led to the identification of 431 DRlncRNAs correlated with disulfidptosis-related genes. Five key DRlncRNAs (SPINT1-AS1, AL161782.1, OVCH1-AS1, AC131009.3, and AC108673.3) were used to develop a prognostic model that effectively distinguished between low- and high-risk patients with significant differences in overall survival and progression-free survival. The low-risk group had a favorable prognosis associated with a protective immune microenvironment and a better response to targeted drugs. Conversely, the high-risk group displayed aggressive tumor features and poor immunotherapy outcomes. Validation through qRT‒PCR confirmed the differential expression of these DRlncRNAs in KIRC cells compared to normal kidney cells, underscoring their potential functional significance in tumor biology.

CONCLUSIONS

This study established a robust link between disulfidptosis-related lncRNAs and patient prognosis in KIRC, underscoring their potential as prognostic biomarkers and therapeutic targets. The differential expression of these lncRNAs in tumor versus normal tissue further highlights their relevance in KIRC pathogenesis. The predictive model not only enhances our understanding of KIRC biology but also provides a novel stratification tool for precision medicine approaches, improving treatment personalization and outcomes in KIRC patients.

Comprehensive Analysis for Predicting Prognoses and Immune Responses of m6A-Related lncRNAs in Women with Lung Adenocarcinoma

During the past decade, the average 5-year survival rate of patients with Lung adenocarcinoma (LUAD) has remained < 20%, although the targeted therapies and novel immunotherapy approaches have held promise. Epigenetic modifications could provide prognostic value as molecular biomarkers, and we aimed to identify the independent risk of m6A-related lncRNAs to establish a risk model for the clinical prediction of prognoses in women with LUAD. In this study, we first assessed 31 N6-methyladenosine (m6A)-related lncRNAs associated with overall survival. Moreover, we evaluated the expression of the oncogenic driver and the tumor immune microenvironment (TIME) in two female LUAD subtypes (clusters 1 and 2) using consensus clustering. We also found 16 m6A-related lncRNAs as the independent prognostic indicator of women with LUAD and established a risk model developed from these lncRNAs. We comprehensively investigated the correlation between the TIME and m6A-related lncRNA and found that m6A-related lncRNA may significantly affect the immune cell infiltration level in LUAD. In conclusion, our study provides evidence on the prognostic prediction in women with LUAD and may help elucidate the processes and mechanisms of m6A-regulated lncRNAs.

Comprehensive landscape of m6A regulator-related gene patterns and tumor microenvironment infiltration characterization in gastric cancer

The epigenetic regulation of N6-methyladenosine (m6A) has attracted considerable interest in tumor research, but the potential roles of m6A regulator-related genes, remain largely unknown within the context of gastric cancer (GC) and tumor microenvironment (TME). Here, a comprehensive strategy of data mining and computational biology utilizing multiple datasets based on 28 m6A regulators (including novel anti-readers) was employed to identify m6A regulator-related genes and patterns and elucidate their underlying mechanisms in GC. Subsequently, a scoring system was constructed to evaluate individual prognosis and immunotherapy response. Three distinct m6A regulator-related patterns were identified through the unsupervised clustering of 56 m6A regulator-related genes (all significantly associated with GC prognosis). TME characterization revealed that these patterns highly corresponded to immune-inflamed, immune-excluded, and immune-desert phenotypes, and their TME characteristics were highly consistent with different clinical outcomes and biological processes. Additionally, an m6A-related scoring system was developed to quantify the m6A modification pattern of individual samples. Low scores indicated high survival rates and high levels of immune activation, whereas high scores indicated stromal activation and tumor malignancy. Furthermore, the m6A-related scores were correlated with tumor mutation loads and various clinical traits, including molecular or histological subtypes and clinical stage or grade, and the score had predictive values across all digestive system tumors and even in all tumor types. Notably, a low score was linked to improved responses to anti-PD-1/L1 and anti-CTLA4 immunotherapy in three independent cohorts. This study has expanded the important role of m6A regulator-related genes in shaping TME diversity and clinical/biological traits of GC. The developed scoring system could help develop more effective immunotherapy strategies and personalized treatment guidance.

Machine Learning Gene Signature to Metastatic ccRCC Based on ceRNA Network

Clear-cell renal-cell carcinoma (ccRCC) is a silent-development pathology with a high rate of metastasis in patients. The activity of coding genes in metastatic progression is well known. New studies evaluate the association with non-coding genes, such as competitive endogenous RNA (ceRNA). This study aims to build a ceRNA network and a gene signature for ccRCC associated with metastatic development and analyze their biological functions. Using data from The Cancer Genome Atlas (TCGA), we constructed the ceRNA network with differentially expressed genes, assembled nine preliminary gene signatures from eight feature selection techniques, and evaluated the classification metrics to choose a final signature. After that, we performed a genomic analysis, a risk analysis, and a functional annotation analysis. We present an 11-gene signature: SNHG15, AF117829.1, hsa-miR-130a-3p, hsa-mir-381-3p, BTBD11, INSR, HECW2, RFLNB, PTTG1, HMMR, and RASD1. It was possible to assess the generalization of the signature using an external dataset from the International Cancer Genome Consortium (ICGC-RECA), which showed an Area Under the Curve of 81.5%. The genomic analysis identified the signature participants on chromosomes with highly mutated regions. The hsa-miR-130a-3p, AF117829.1, hsa-miR-381-3p, and PTTG1 were significantly related to the patient's survival and metastatic development. Additionally, functional annotation resulted in relevant pathways for tumor development and cell cycle control, such as RNA polymerase II transcription regulation and cell control. The gene signature analysis within the ceRNA network, with literature evidence, suggests that the lncRNAs act as "sponges" upon the microRNAs (miRNAs). Therefore, this gene signature presents coding and non-coding genes and could act as potential biomarkers for a better understanding of ccRCC.

Comprehensive Analysis of N6-Methylandenosine-Related lncRNAs in Clear Cell Renal Cell Carcinoma: A Correlation With Prognosis, Tumor Progression, and Therapeutic Response

This study aims to develop a prognostic signature based on m6A-related lncRNAs for clear cell renal cell carcinoma (ccRCC). Differential expression analysis and Pearson correlation analysis were used to identify m6A-related lncRNAs associated with patient outcomes in The Cancer Genome Atlas (TCGA) database. Our approach led to the development of an m6A-related lncRNA risk score (MRLrisk), formulated using six identified lncRNAs: NFE4, AL008729.2, AL139123.1, LINC02154, AC124854.1 and ARHGAP31-AS1. Higher MRLrisk was identified as a risk factor for patients' prognosis in ccRCC. Furthermore, an MRLrisk-based nomogram was developed and demonstrated as a reliable tool for prognosis prediction in ccRCC. Enrichment analysis and tumor mutation signature studies were conducted to investigate MRLrisk-related biological phenotypes. The tumor immune dysfunction and exclusion (TIDE) score was employed to infer patients' response to immunotherapy, indicating a negative correlation between high MRLrisk and immunotherapy response. Our focus then shifted to LINC02154 for deeper exploration. We assessed LINC02154 expression in 28 ccRCC/normal tissue pairs and 3 ccRCC cell lines through quantitative real-time polymerase chain reaction (qRT-PCR). Functional experiments, including EdU incorporation, flow cytometry and transwell assays, were performed to assess the role of LINC02154 in ccRCC cell functions, discovering that its downregulation hinders cancer cell proliferation and migration. Furthermore, the influence of LINC02154 on ccRCC cells' sensitivity to Sunitinib was explored using CCK-8 assays, demonstrating that decreased LINC02154 expression increases Sunitinib sensitivity. In summary, this study successfully developed an MRLrisk model with significant prognostic value for ccRCC and established LINC02154 as a critical biomarker and prospective therapeutic target in ccRCC management.

Cuproptosis gene-related, neural network-based prognosis prediction and drug-target prediction for KIRC

BACKGROUND

Kidney renal clear cell carcinoma (KIRC), as a common case in renal cell carcinoma (RCC), has the risk of postoperative recurrence, thus its prognosis is poor and its prognostic markers are usually based on imaging methods, which have the problem of low specificity. In addition, cuproptosis, as a novel mode of cell death, has been used as a biomarker to predict disease in many cancers in recent years, which also provides an important basis for prognostic prediction in KIRC. For postoperative patients with KIRC, an important means of preventing disease recurrence is pharmacological treatment, and thus matching the appropriate drug to the specific patient's target is also particularly important. With the development of neural networks, their predictive performance in the field of medical big data has surpassed that of traditional methods, and this also applies to the field of prognosis prediction and drug-target prediction.

OBJECTIVE

The purpose of this study is to screen for cuproptosis genes related to the prognosis of KIRC and to establish a deep neural network (DNN) model for patient risk prediction, while also developing a personalized nomogram model for predicting patient survival. In addition, sensitivity drugs for KIRC were screened, and a graph neural network (GNN) model was established to predict the targets of the drugs, in order to discover potential drug action sites and provide new treatment ideas for KIRC.

METHODS

We used the Cancer Genome Atlas (TCGA) database, International Cancer Genome Consortium (ICGC) database, and DrugBank database for our study. Differentially expressed genes (DEGs) were screened using TCGA data, and then a DNN-based risk prediction model was built and validated using ICGC data. Subsequently, the differences between high- and low-risk groups were analyzed and KIRC-sensitive drugs were screened, and finally a GNN model was trained using DrugBank data to predict the relevant targets of these drugs.

RESULTS

A prognostic model was built by screening 10 significantly different cuproptosis-related genes, the model had an AUC of 0.739 on the training set (TCGA data) and an AUC of 0.707 on the validation set (ICGC data), which demonstrated a good predictive performance. Based on the prognostic model in this paper, patients were also classified into high- and low-risk groups, and functional analyses were performed. In addition, 251 drugs were screened for sensitivity, and four drugs were ultimately found to have high sensitivity, with 5-Fluorouracil having the best inhibitory effect, and subsequently their corresponding targets were also predicted by GraphSAGE, with the most prominent targets including Cytochrome P450 2D6, UDP-glucuronosyltransferase 1A, and Proto-oncogene tyrosine-protein kinase receptor Ret. Notably, the average accuracy of GraphSAGE was 0.817 ± 0.013, which was higher than that of GAT and GTN.

CONCLUSION

Our KIRC risk prediction model, constructed using 10 cuproptosis-related genes, had good independent prognostic ability. In addition, we screened four highly sensitive drugs and predicted relevant targets for these four drugs that might treat KIRC. Finally, literature research revealed that four drug-target interactions have been demonstrated in previous studies and the remaining targets are potential sites of drug action for future research.

LncRNA AGAP2 antisense RNA 1 stabilized by insulin-like growth factor 2 mRNA binding protein 3 promotes macrophage M2 polarization in clear cell renal cell carcinoma through regulation of the microRNA-9-5p/THBS2/PI3K-Akt pathway

BACKGROUND

Increasing evidence highlights the potential role of long non-coding RNAs (lncRNAs) in the biological behaviors of renal cell carcinoma (RCC). Here, we explored the mechanism of AGAP2-AS1 in the occurrence and development of clear cell RCC (ccRCC) involving IGF2BP3/miR-9-5p/THBS2.

METHODS

The expressions of AGAP2-AS1, IGF2BP3, miR-9-5p, and THBS2 and their relationship were analyzed by bioinformatics. The targeting relationship between AGAP2-AS1 and miR-9-5p and between miR-9-5p and THBS2 was evaluated with their effect on cell biological behaviors and macrophage polarization assayed. Finally, we tested the effect of AGAP2-AS1 on ccRCC tumor formation in xenograft tumors.

RESULTS

IGF2BP3 could stabilize AGAP2-AS1 through m6A modification. AGAP2-AS1 was highly expressed in ccRCC tissues and cells. The lentivirus-mediated intervention of AGAP2-AS1 induced malignant behaviors of ccRCC cells and led to M2 polarization of macrophages. In addition, THBS2 promoted M2 polarization of macrophages by activating the PI3K/AKT signaling pathway. AGAP2-AS1 could directly bind with miR-9-5p and promote the expression of THBS2 downstream of miR-9-5p. These results were further verified by in vivo experiments.

CONCLUSION

AGAP2-AS1 stabilized by IGF2BP3 competitively binds to miR-9-5p to up-regulate THBS2, activating the PI3K/AKT signaling pathway and inducing macrophage M2 polarization, thus facilitating the development of RCC.

N6-methyladenosine methylation in kidney injury

Multiple mechanisms are involved in kidney damage, among which the role of epigenetic modifications in the occurrence and development of kidney diseases is constantly being revealed. However, N6-methyladenosine (M6A), a well-known post-transcriptional modification, has been regarded as the most prevalent epigenetic modifications in higher eukaryotic, which is involved in various biological processes of cells such as maintaining the stability of mRNA. The role of M6A modification in the mechanism of kidney damage has attracted widespread attention. In this review, we mainly summarize the role of M6A modification in the progression of kidney diseases from the following aspects: the regulatory pattern of N6-methyladenosine, the critical roles of N6-methyladenosine in chronic kidney disease, acute kidney injury and renal cell carcinoma, and then reveal its potential significance in the diagnosis and treatment of various kidney diseases. A better understanding of this field will be helpful for future research and clinical treatment of kidney diseases.

Prognostic model of kidney renal clear cell carcinoma using aging-related long noncoding RNA signatures identifies THBS1-IT1 as a potential prognostic biomarker for multiple cancers

Aging is responsible for the main intrinsic triggers of cancers; however, the studies of aging risk factors in cancer animal models and cancer patients are rare and insufficient to be represented in cancer clinical trials. For a better understanding of the complex regulatory networks of aging and cancers, 8 candidate aging related long noncoding RNAs (CarLncs) identified from the healthy aging models, centenarians and their offsprings, were selected and their association with kidney renal clear cell carcinoma (KIRC) was explored by series of cutting edge analyses such as support vector machine (SVM) and random forest (RF) algorithms. Using data downloaded from TCGA and GTEx databases, a regulatory network of CarLncs-miRNA-mRNA was constructed and five genes within the network were screened out as aging related feature genes for developing KIRC prognostic models. After a strict filtering pipeline for modeling, a formula using the transcript per million (TPM) values of feature genes "LncAging_score = 0.008* MMP11 + 0.066* THBS1-IT1 + (-0.014)* DYNLL2 + (-0.030)* RMND5A+ 0.008* PEG10" was developed. ROC analysis and nomogram suggest our model achieves a great performance in KIRC prognosis. Among the 8 CarLncs, we found that THBS1-IT1 was significantly dysregulated in 12 cancer types. A comprehensive pan-cancer analysis demonstrated that THBS1-IT1 is a potential prognostic biomarker in not only KIRC but also multiple cancers, such as LUSC, BLCA, GBM, LGG, MESO, PAAD, STAD and THCA, it was correlated with tumor microenvironment (TME) and tumor immune cell infiltration (TICI) and its high expression was related with poor survival.

A novel oxidative-stress related lncRNA signature predicts the prognosis of clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is a primary malignant tumour of tubular epithelial origin and is most common in the urinary tract. Growing evidence suggests that oxidative stress (OS), generates high levels of reactive oxygen species (ROS) and free radicals, and plays a critical role in cancer in humans. However, the predictive value of OS-related long non-coding RNAs (lncRNAs) in ccRCC remains unclear. We constructed a predictive signature of survival based on OS-related lncRNAs that were obtained from The Cancer Genome Atlas (TCGA-KIRC), to predict the prognosis of patients with ccRCC. The signature comprised seven lncRNAs: SPART-AS1, AL162586.1, LINC00944, LINC01550, HOXB-AS4, LINC02027, and DOCK9-DT. OS-related signature of lncRNAs had diagnostic efficiency higher than that of clinicopathological variables, with an area of 0.794 under the receiver operating characteristic curve. Additionally, the nomogram based on risk scores and clinicopathological variables (age, gender, grade, stage, M-stage, and N-stage) showed strong predictive performance. Patients with high-risk were found to be more sensitive to the therapeutic drugs ABT.888, AICAR, MS.275, sunitinib, AZD.2281, and GDC.0449. Our constructed the predictive signature can independently predict the prognosis of patients with ccRCC; however, the underlying mechanism needs further investigation.

SERPINE1 and its co-expressed genes are associated with the progression of clear cell renal cell carcinoma

BACKGROUND

Clear cell renal cell carcinoma(ccRCC) is a frequently occurring malignant tumor of the urinary system. Despite extensive research, the regulatory mechanisms underlying the pathogenesis and progression of ccRCC remain largely unknown.

METHODS

We downloaded 5 ccRCC expression profiles from the Gene Expression Omnibus (GEO) database and obtained the list of differentially expressed genes (DEGs). Using String and Cytoscape tools, we determined the hub genes of ccRCC, and then analyzed their relationship with ccRCC patient survival. Ultimately, we identified SERPINE1 as a prognostic factor in ccRCC. Meanwhile, we confirmed the role of SERPINE1 in 786-O cells by cell transfection and in vitro experiments.

RESULTS

Our analysis yielded a total of 258 differentially expressed genes, comprising 105 down-regulated genes and 153 up-regulated genes. Survival analysis of SERPINE1 expression in The Cancer Genome Atlas (TCGA) confirmed its association with the increase of tumor grade, lymph node metastasis, and tumor stage, as well as with shorter survival. Furthermore, we found that SERPINE1 expression levels were associated with CD8 + T cells, CD4 + T cells, B cells, macrophages, neutrophils, and dendritic cells. Cell experiments showed that knockdown SERPINE1 expression could inhibit the proliferation, migration and invasion of ccRCC cells. Among the co-expressed genes with the highest correlation, ITGA5, SLC2A3, SLC2A14, SHC1, CEBPB, and ADA were overexpressed and associated with shorter overall survival (OS) in ccRCC.

CONCLUSIONS

In this study, we identified hub genes that are strongly related to ccRCC, and highlights the potential utility of overexpressed SERPINE1 and its co-expressed genes could be used as prognostic and diagnostic biomarkers in ccRCC.

Constructing and validating of m6a-related genes prognostic signature for stomach adenocarcinoma and immune infiltration: Potential biomarkers for predicting the overall survival

Background

Stomach adenocarcinoma (STAD) arises from the mutations of stomach cells and has poor overall survival. Chemotherapy is commonly indicated for patients with stomach cancer following surgical resection. The most prevalent alteration that affects cancer growth is N6-methyladenosine methylation (m6A), although the possible function of m6A in STAD prognosis is not recognized.

Method

The research measured predictive FRGs in BLCA samples from the TCGA and GEO datasets. Data on the stemness indices (mRNAsi), gene mutations, copy number variations (CNV), tumor mutation burden (TMB), and corresponding clinical characteristics were obtained from TCGA and GEO. STAD from TCGA and GEO at 24 m6A was investigated. Lasso regression was used to construct the prediction model to assess the m6A prognostic signals in STAD. In addition, the correlation between m6a and immune infiltration in STAD patients was discussed using GSVA and ssGSEA analysis. Based on these genes, GO and KEGG analyses were performed to identify key biological functions and key pathways.

Result

A significant relationship was discovered between numerous m6A clusters and the tumor immune microenvironment, as well as three m6A alteration patterns with different clinical outcomes. Furthermore, GSVA and ssGSEA showed that m6A clusters were significantly associated with immune infiltration in the STAD. The low-m6Ascore group had a lower immunotherapeutic response than the high-m6Ascore group. ICIs therapy was more effective in the group with a higher m6Ascore. Three writers (VIRMA, ZC3H13, and METTL3) showed significantly lower expression, whereas five authors (METTL14, METTL16, WTAP, RBM15, and RBM15B) showed considerably higher expression. Three readers (YTHDC2, YTHDF2, and LRPPRC) had higher levels of expression, whereas eleven readers (YTHDC1, YTHDF1, YTHDF3, HNRNPC, FMR1, HNRNPA2B1, IGFBP1, IGFBP2, IGFBP3, and RBMX) had lower levels. As can be observed, the various types of m6 encoders have varied ramifications for STAD control.

Conclusion

STAD occurrence and progression are linked to m6A-genes. Corresponding prognostic models help forecast the prognosis of STAD patients. m6A-genes and associated immune cell infiltration in the tumor microenvironment (TME) may serve as potential therapeutic targets in STAD, which requires further trials. In addition, the m6a-related gene signature offers a viable alternative to predict bladder cancer, and these m6A-genes show a prospective research area for STAD targeted treatment in the future.

Prognostic value and immunological characteristics of a novel cuproptosis-related long noncoding RNAs risk signature in kidney renal clear cell carcinoma

Background: Cuproptosis has been found as a novel cell death mode significantly associated with mitochondrial metabolism, which may be significantly associated with the occurrence and growth of tumors. LncRNAs take on critical significance in regulating the development of kidney renal clear cell carcinoma (KIRC), whereas the correlation between cuproptosis-related LncRNAs (CRLs) and KIRC is not clear at present. Therefore, this study built a prognosis signature based on CRLs, which can achieve accurate prediction of the outcome of KIRC patients. Methods: The TCGA database provided the expression profile information and relevant clinical information of KIRC patients. Univariate Cox, Lasso, and multivariate Cox were employed for building a risk signature based on CRLs. Kaplan-Meier (K-M) survival analysis and time-dependent receiver operating characteristic (ROC) curve were employed for the verification and evaluation of the reliability and accuracy of risk signature. Then, qRT-PCR analysis of risk LncRNAs was conducted. Finally, the possible effect of the developed risk signature on the microenvironment for tumor immunization was speculated in accordance with ssGSEA and ESTIMATE algorithms. Results: A prognosis signature composed of APCDD1L-DT, MINCR, AL161782.1, and AC026401.3 was built based on CRLs. As revealed by the results of the K-M survival study, the OS rate and progression-free survival rate of highrisk KIRC patients were lower than those of lowrisk KIRC patients, and the areas under ROC curves of 1, 3, and 5 years were 0.828, 0.780, and 0.794, separately. The results of the immune analysis showed that there were significant differences in the status of immunization and the microenvironment of tumor between groups at low-risk and at high-risk. The qRT-PCR results showed that the relative expression level of MINCR and APCDD1L-DT were higher in 786-O and 769-P tumor cells than in HK-2 cells, which were normal renal tubular epithelial cells. Conclusion: The developed risk signature takes on critical significance in the prediction of the prognosis of patients with KIRC, and it can bring a novel direction for immunotherapy and clinical drug treatment of KIRC. In addition, 4 identified risk LncRNAs (especially APCDD1L-DT and MINCR) can be novel targets for immunotherapy of KIRC patients.

Clinical significance and immune landscape of cuproptosis-related lncRNAs in kidney renal clear cell carcinoma: a bioinformatical analysis

Background

Kidney renal clear cell carcinoma (KIRC) is considered an immunogenic tumor. Cuproptosis is a newly identified copper-induced regulated cell death that relies on mitochondria respiration. Long noncoding RNAs (lncRNAs) have emerged as significant players in tumorigenesis and metastasis. However, there is a huge knowledge gap on the prognostic role of cuproptosis-related lncRNAs in KIRC. And, the clinical value of them is still unknown. Here, we aimed to develop a cuproptosis-related lncRNA prognostic signature in KIRC.

Methods

The messenger RNA (mRNA)/lncRNA expression profiles and the clinical information including age, gender, tumor stage, grade, and overall survival (OS) were acquired from The Cancer Genome Atlas (TCGA) database. The included KIRC samples were further randomly assigned into training (n=258) or testing (n=257) data sets. We performed Pearson correlation analysis to identify the cuproptosis-related lncRNAs and then constructed the prognostic signature using Cox regression analysis and LASSO algorithm. Subsequently, Kaplan-Meier survival analysis, a nomogram, and receiver operating characteristic (ROC) curve were performed to assess the predictive performance of the signature. Moreover, the immune characteristics and drug sensitivity related to the signature were also explored.

Results

The signature comprised 7 cuproptosis-related lncRNAs. The patients with a low-risk score had superior OS compared with those with a high-risk score. The survival rates of the high- and low-risk groups were 44.96% and 83.72% (P<0.001). The area under the curve (AUC) value for 1-, 3-, 5-year survival rate reached 0.814, 0.762 and 0.825, respectively. In addition, a nomogram was also generated; the AUC was 0.785 for risk score, higher than that for age (0.593), gender (0.489), grade (0.679), and stage (0.721). The high-risk group had more enriched immune- and tumor-related genes. Patients with low-risk scores were more sensitive to immunotherapy and the small molecular drugs GSK1904529A, tipifarnib, BX-912, FR-180204, and GSK1070916. Meanwhile, the high-risk group tended to be more sensitive to pyrimethamine, MS-275, and CGP-60474.

Conclusions

Collectively, we constructed a cuproptosis-related lncRNA prognostic signature with a higher predictive accuracy compared to multiple clinicopathological parameters, which may provide vital guidance for therapeutic strategies in KIRC. Combination of more prognostic biomarkers may further improve the accuracy.

Comprehensive Multiomic Analysis Identified TUBA1C as a Potential Prognostic Biological Marker of Immune-Related Therapy in Pan-Cancer

TUBA1C is correlated with an unfavourable prognosis and the infiltration of immune cells in several cancers. However, its function as a significant biomarker for the prognosis of immunotherapy in pan-cancer remains unclear. This study aims at assessing the role of TUBA1C in pan-cancer at multiple levels, including mutations, gene expression, methylation, m6A methylation, and immune cell infiltration levels. Data retrieved from major public databases, such as TCGA, GEO, GTEx, GSCA, CancerSEA, HPA, and RNAactDrugs, revealed that TUBA1C expression was high in 33 cancer types. Survival analysis revealed that TUBA1C was a poor prognostic factor for 12 tumour types, and mutations, CNVs, and methylation affected the prognosis of some cancer types. Furthermore, TUBA1C was found to be related to immune-related genes, immune cell infiltration, and the immune microenvironment. In addition, the sensitivity of 10 anticancer drugs was associated with high TUBA1C expression. Therefore, TUBA1C may serve as a viable prognostic biomarker for immunotherapy of pan-cancer.

Calreticulin as a prognostic biomarker and correlated with immune infiltrate in kidney renal clear cell carcinoma

Background: Calreticulin (CALR) has been investigated in several malignant diseases and is associated with immune-cell infiltration. However, the prognostic value of CALR in kidney renal clear cell carcinoma (KIRC) is still unknown. Methods: Based on the computational analysis, data from 530 KIRC cases and 72 normal kidney samples from The Cancer Genome Atlas (TGCA-KIRC) database were analyzed in this study. The expression of CALR mRNA in pan-cancer and immune infiltrates was analyzed using the Tumor Immune Estimation Resource (TIMER) database. The CALR protein expression was obtained from the UALCAN and Human Protein Atlas (HPA) databases. Survival, functional, and statistical analyses were conducted using R software. Results: The CALR expression was higher in KIRC cases than in normal kidneys. A high CALR expression was correlated with TNM stage, pathological stage, and histological grade. Kaplan-Meier survival analysis showed that a high CALR expression was associated with poor overall survival, disease-specific survival, and progression-free interval. Gene set enrichment analysis (GSEA) indicated that CALR was enriched in IL-6 and IL-2 signaling, interferon signaling, TNF signaling, inflammatory response, apoptosis, and the p53 pathway. CALR is correlated with immune-infiltrating cells. A significant correlation was observed between CALR expression and immunomodulators. Conclusion: We identified CALR as a prognostic biomarker of KIRC. Meanwhile, the CALR expression associated with immune infiltration indicated that CALR might be a potential immunotherapy target for patients with KIRC.

Chromatin regulators-related lncRNA signature predicting the prognosis of kidney renal clear cell carcinoma and its relationship with immune microenvironment: A study based on bioinformatics and experimental validation

Background: Kidney Renal Clear cell carcinoma (KIRC) is a major concern in the urinary system. A lot of researches were focused on Chromatin Regulators (CRs) in tumors. In this study, CRs-related lncRNAs (CRlncRNAs) were investigated for their potential impact on the prognosis of KIRC and the immune microenvironment.

Methods: The TCGA database was used to obtain transcriptome and related clinical information. CRs were obtained from previous studies, whereas CRlncRNAs were obtained by differential and correlation analysis. We screened the lncRNAs for the signature construction using regression analysis and LASSO regression analysis. The effectiveness of the signature was evaluated using the Kaplan-Meier (K-M) curve and Receiver Operating Characteristic curve (ROC). Additionally, we examined the associations between the signature and Tumor Microenvironment (TME), and the efficacy of drug therapy. Finally, we further verified whether these lncRNAs could affect the biological function of KIRC cells by functional experiments such as CCK8 and transwell assay.

Results: A signature consisting of 8 CRlncRNAs was constructed to predict the prognosis of KIRC. Quantitative Real-Time PCR verified the expression of 8 lncRNAs at the cell line and tissue level. The signature was found to be an independent prognostic indicator for KIRC in regression analysis. This signature was found to predict Overall Survival (OS) better for patients in the subgroups of age, gender, grade, stage, M, N0, and T. Furthermore, a significant correlation was found between riskScore and immune cell infiltration and immune checkpoint. Finally, we discovered several drugs with different IC50 values in different risk groups using drug sensitivity analysis. And functional experiments showed that Z97200.1 could affect the proliferation, migration and invasion of KIRC cells.

Conclusion: Overall, the signature comprised of these 8 lncRNAs were reliable prognostic biomarkers for KIRC. Moreover, the signature had significant potential for assessing the immunological landscape of tumors and providing individualized treatment.

Identification of N6-methylandenosine related lncRNA signatures for predicting the prognosis and therapy response in colorectal cancer patients

Despite recent advances in surgical and multimodal therapies, the overall survival (OS) of advanced colorectal cancer (CRC) patients remains low. Thus, discerning sensitive prognostic biomarkers to give the optimistic treatment for CRC patients is extremely critical. N6-methyladenosine (m6A) and long noncoding RNAs (lncRNAs) play an important role in CRC progression. Nonetheless, few studies have focused on the impact of m6A-related lncRNAs on the prognosis, tumor microenvironment (TME) and treatment of CRC. In this study, 1707 m6A-related lncRNAs were identified through Pearson correlation analysis and Weighted co-expression network analysis (WGCNA) using The Cancer Genome Atlas (TCGA) cohort. Then, 28 m6A-related prognostic lncRNAs were screened by univariate Cox regression analysis, followed by identifying two clusters by consensus clustering analysis. A prognostic model consisted of 8 lncRNA signatures was constructed by the least absolute shrinkage and selection operator (LASSO). Kaplan–Meier curve analysis and a nomogram were performed to investigate the prognostic ability of this model. The risk score of prognostic model act as an independent risk factor for OS rate. Functional enrichment analysis indicated that lncRNA signatures related tumor immunity. The low-risk group characterized by increased microsatellite instability-high (MSI-H), mutation burden, and immunity activation, indicated favorable odds of OS. Moreover, the lncRNA signatures were significantly associated with the cancer stem cell (CSC) index and drug sensitivity. In addition, 3 common immune genes shared by the lncRNA signatures were screened out. We found that these immune genes were widely distributed in 2 cell types of TME. Finally, a ceRNA network was constructed to identify ZEB1-AS1 regulatory axis in CRC. We found that ZEB1-AS1 was significantly overexpressed in tumor tissues, and was related to the metastasis of EMT and the chemoresistance of 5-Fu in CRC. Therefore, our study demonstrated the important role of m6A-related lncRNAs in TME remodeling. Moreover, these results illustrated the levels of ZEB1-AS1 might be valuable for predicting the progression and prognosis of CRC, and further provided a new target for the diagnosis and treatment of CRC patients.

Pan-Cancer Gene Analysis of m6A Modification and Immune Infiltration in Uterine Corpus Endometrial Carcinoma

Objective

This investigation was to test the potential role of m6A-related long non-coding RNAs (lncRNAs) and immune infiltration as crucial factors in the diagnosis and treatment of uterine corpus endometrial cancer (UCEC).

Method

The UCEC RNA-seq data were downloaded in the Cancer Genome Atlas (TCGA, https://portal.gdc.cancer.gov/). There were 587 samples totally, containing 543 UCEC cases and 35 healthy cases. The clinical information of UCEC cases included survival time, survival status, gender, age, stage, and TMN stage. Twenty-three m6A-related genes were found in published journals. The RNA-seq documents of UCEC were downloaded in the Cancer Genome Atlas (TCGA). The hub gene data of UCEC were downloaded from GEPIA2 database. The different packages of R language were applied to calculate and analyze in this research.

Results

Among 587 cases in our study, we discovered 3039 lncRNAs in the TCGA-UCEC database. After the differential analysis, 23 m6A-associated genetics were screened and twenty-one m6A-associated differential genetics were found. In the end, we obtained 20 m6A-related lncRNAs. LNCTAM34A was considered as a predictive gene through univariate and multivariate Cox regression analysis. In addition to the above, patients with high LNCTAM34A expression had better outcomes than those with low LNCTAM34A expression. The high-risk cohort had greater scores of activated dendritic cells (aDCs), B cells, and T cell regulatory (Tregs) than low-risk cohort; in the meanwhile, high-risk cohort had lower scores of DCs and iDCs. Then, the high-risk cohort displayed greater scores in the immune functions of MHC class I, para-inflammation, and type I IFN response than those of low-risk cohort. Among 27 immune-inducible genes, the level of CD244, KIR3DLI, NRP1, PDCD1LG2, and TNFRSF8 was reduced in UCEC samples and the level of CD27, CD28, CD70, CD80, CD86, HAVCR2, ICOS, IDO1, LAIR1, PDCD1, TIGIT, TNFRSF18, -25, -9, -14, and VTCN1 was increased in UCEC samples.

Conclusion

The key role of M6A-related lncRNAs in immune microenvironment in high-risk patients of UCEC. The patients with strong expression of LNCTAM34A have a good prognosis, and LNCTAM34A can be used as a prognostic gene for UCEC. m6A-related lncRNAs can be used as a potential treatment for UCEC. Our observations can be used as a hypothetical basis for future in vitro and animal experiments.

A novel marker based on necroptosis-related long non-coding RNA for forecasting prognostic in patients with clear cell renal cell carcinoma

Background: The incidence of clear cell renal cell carcinoma (ccRCC) is high and has increased gradually in recent years. At present, due to the lack of effective prognostic indicators, the prognosis of ccRCC patients is greatly affected.Necroptosis is a type of cell death, and along with cell necrosis is considered a new cancer treatment strategy. The aim of this study was to construct a new marker for predicting the prognosis of ccRCC patients based on long non-coding RNA (nrlncRNAs) associated with necroptosis.

Methods: RNA sequence data and clinical information of ccRCC patients from the Cancer Genome Atlas database (TCGA) were downloaded. NrlncRNA was identified by Pearson correlation study. The differentially expressed nrlncRNA and nrlncRNA pairs were identified by univariate Cox regression and Lasso-Cox regression. Finally, a Kaplan-Meier survival study, Cox regression, clinicopathological features correlation study, and receiver operating characteristic (ROC) spectrum were used to evaluate the prediction ability of 25-nrlncrnas for markers. In addition, correlations between the risk values and sensitivity to tumor-infiltrating immune cells, immune checkpoint inhibitors, and targeted drugs were also investigated.

Results: In the current research, a novel marker of 25-nrlncRNAs pairs was developed to improve prognostic prediction in patients with ccRCC. Compared with clinicopathological features, nrlncRNAs had a higher diagnostic validity for markers, with the 1-year, 3-years, and 5-years operating characteristic regions being 0.902, 0.835, and 0.856, respectively, and compared with the stage of 0.868, an increase of 0.034. Cox regression and stratified survival studies showed that this marker could be an independent predictor of ccRCC patients. In addition, patients with different risk scores had significant differences in tumor-infiltrating immune cells, immune checkpoint, and semi-inhibitory concentration of targeted drugs. The feature could be used to evaluate the clinical efficacy of immunotherapy and targeted drug therapy.

Conclusion: 25-nrlncRNAs pair markers may help to evaluate the prognosis and molecular characteristics of ccRCC patients, which improve treatment methods and can be more used in clinical practice.

Novel amino acid metabolism‐related gene signature to predict prognosis in clear cell renal cell carcinoma

Background: Amino acid metabolism (AAM) deregulation, an emerging metabolic hallmark of malignancy, plays an essential role in tumour proliferation, invasion, and metastasis. However, the expression of AAM-related genes and their correlation with prognosis in clear cell renal cell carcinoma (ccRCC) remain elusive. This study aims to develop a novel consensus signature based on the AAM-related genes.

Methods: The RNA-seq expression data and clinical information for ccRCC were downloaded from the TCGA (KIRC as training dataset) and ArrayExpress (E-MTAB-1980 as validation dataset) databases. The AAM‐related differentially expressed genes were screened via the “limma” package in TCGA cohorts for further analysis. The machine learning algorithms (Lasso and stepwise Cox (direction = both)) were then utilised to establish a novel consensus signature in TCGA cohorts, which was validated by the E-MTAB-1980 cohorts. The optimal cutoff value determined by the “survminer” package was used to categorise patients into two risk categories. The Kaplan-Meier curve, the receiver operating characteristic (ROC) curve, and multivariate Cox regression were utilised to evaluate the prognostic value. The nomogram based on the gene signature was constructed, and its performance was analysed using ROC and calibration curves. Gene Set Enrichment Analysis (GSEA) and immune cell infiltration analysis were conducted on its potential mechanisms. The relationship between the gene signature and key immune checkpoint, N6-methyladenosine (m⁶A)-related genes, and sensitivity to chemotherapy was assessed.

Results: A novel consensus AMM‐related gene signature consisting of IYD, NNMT, ACADSB, GLDC, and PSAT1 is developed to predict prognosis in TCGA cohorts. Kaplan-Meier survival shows that overall survival in the high-risk group was more dismal than in the low-risk group in the TCGA cohort, validated by the E-MTAB-1980 cohort. Multivariate regression analysis also demonstrates that the gene signature is an independent predictor of ccRCC. Immune infiltration analysis highlighted that the high-risk group indicates an immunosuppressive microenvironment. It is also closely related to the level of key immune checkpoints, m⁶A modification, and sensitivity to chemotherapy drugs.

Conclusion: In this study, a novel consensus AAM-related gene signature is developed and validated as an independent predictor to robustly predict the overall survival from ccRCC, which would further improve the clinical outcomes.

Establishment of a prognostic ferroptosis- and immune-related long noncoding RNAs profile in kidney renal clear cell carcinoma

Background: Ferroptosis and immunity are novel treatments that target several cancers, including kidney renal clear cell carcinoma (KIRC). Long noncoding RNAs (lncRNAs) are an important class of gene expression regulators that play fundamental roles in the regulation of ferroptosis and immunity. We aimed to identify ferroptosis- and immune-related lncRNAs as biomarkers in patients with KIRC.

Methods: Corresponding data for each patient with KIRC were obtained from The Cancer Genome Atlas (TCGA) database. Univariate and multivariate Cox regression analyses were used to identify candidate biomarkers followed by least absolute shrinkage and selection operator (LASSO) regression analyses, weighted gene coexpression network analysis (WGCANA), and gene set enrichment analysis (GSEA).

Results: Three ferroptosis- and immune-related differentially expressed lncRNAs (FI-DELs) (AC124854.1, LINC02609, and ZNF503-AS2) were markedly and independently correlated with the overall survival (OS) of patients with KIRC. The area under the curve (AUC) value of the prognostic model in the entire group using the three FI-DELs was > 0.70. The sensitivity and specificity of the diagnostic model using the three FI-DELs were 0.8586 and 0.9583, respectively.

Conclusion: The present study found that AC124854.1, LINC02609, and ZNF503-AS2 were considerably and independently correlated with the OS of patients with KIRC, suggesting that the three FI-DELs could be used as prognostic and diagnostic biomarkers for patients with KIRC.

The integrated comprehension of lncRNA HOXA-AS3 implication on human diseases

Long non-coding RNA (lncRNA) is a non-protein-coding RNA with a length of more than 200 nucleotides. Studies have shown that lncRNAs have vital impacts on various pathological processes and participate in the development of human diseases, usually through acting as competing endogenous RNAs to modulate miRNA expression and biological functions. lncRNA HOXA Cluster Antisense RNA 3 (HOXA-AS3) was a newly discovered lncRNA and has been demonstrated to be abnormally expressed in many diseases. Moreover, HOXA-AS3 expression was closely correlated with the clinicopathologic characteristics in cancer patients. In addition, HOXA-AS3 exhibited significant properties in regulating several biological processes, including cell proliferation, invasion, and migration. Furthermore, HOXA-AS3 has provided promising values in the diagnosis, prognosis, and therapeutic strategies of several diseases such as liver cancer, glioma, lung cancer, oral cancer, gastric cancer, and even atherosclerosis. In this review, we discuss the abnormal expression of HOXA-AS3 in several human disorders and some pathobiological processes and its clinical characteristics, followed by a summary of HOXA-AS3 functions, regulatory mechanisms, and clinical application potential.

Identification and Validation of an m6A-Related LncRNA Signature to Predict Progression-Free Survival in Colorectal Cancer

The RNA methylation of N6 adenosine (m6A) plays a crucial role in various biological processes. Strong evidence reveals that the dysregulation of long non-coding RNAs (lncRNA) brings about the abnormality of downstream signaling in multiple ways, thus influencing tumor initiation and progression. Currently, it is essential to discover effective and succinct molecular biomarkers for predicting colorectal cancer (CRC) prognosis. However, the prognostic value of m6A-related lncRNAs for CRC remains unclear, especially for progression-free survival (PFS). Here, we screened 24 m6A-related lncRNAs in 622 CRC patients and identified five lncRNAs (SLCO4A1-AS1, MELTF-AS1, SH3PXD2A-AS1, H19 and PCAT6) associated with patient PFS. Compared to normal samples, their expression was up-regulated in CRC tumors from TCGA dataset, which was validated in 55 CRC patients from our in-house cohort. We established an m6A-Lnc signature for predicting patient PFS, which was an independent prognostic factor by classification analysis of clinicopathologic features. Moreover, the signature was validated in 1,077 patients from six independent datasets (GSE17538, GSE39582, GSE33113, GSE31595, GSE29621, and GSE17536), and it showed better performance than three known lncRNA signatures for predicting PFS. In summary, our study demonstrates that the m6A-Lnc signature is a promising biomarker for forecasting patient PFS in CRC.

LTBP2 inhibits prostate cancer progression and metastasis via the PI3K/AKT signaling pathway

Biochemical recurrence (BCR) is a cause of concern in advanced prostate cancer (PCa). Thus, novel diagnostic biomarkers are required to improve clinical care. However, research on PCa immunotherapy is also scarce. Hence, the present study aimed to explore promising BCR-related diagnostic biomarkers, and their expression pattern, prognostic value, immune response effects, biological functions, and possible molecular mechanisms were evaluated. GEO datasets (GSE46602, GSE70768, and GSE116918) were downloaded and merged as the training cohort, and differential expression analysis was performed. Lasso regression and SVM-RFE algorithm, as well as PPI analysis and MCODE algorithm, were then applied to filter BCR-related biomarker genes. The CIBERSORT and estimation of stromal and immune cells in malignant tumor tissues using expression data (ESTIMATE) methods were used to calculate the fractions of tumor-infiltrating immune cells. GO/DO enrichment analyses were used to identify the biological functions. The expression of latent transforming growth factor β-binding protein 2 (LTBP2) was determined by RT-qPCR and western blotting. The role of LTBP2 in PCa was determined by CCK-8, Transwell, and the potential mechanism was investigated by KEGG and GSEA and confirmed by western blotting. In total, 44 BCR-related differentially expressed genes (DEGs) in the training cohort were screened. LTBP2 was found to be a diagnostic biomarker of BCR in PCa and was associated with CD4⁺ T-cell infiltration and response to anti-PD-1/PD-L1 immunotherapy. Subsequently, using the ESTIMATE algorithm, it was identified that LTBP2 was associated with the tumor microenvironment and could be a predictor of the clinical benefit of immune checkpoint blockade. Finally, the expression and biological function of LTBP2 were evaluated via cellular experiments. The results showed that LTBP2 was downregulated in PCa cells and inhibited PCa proliferation and metastasis via the PI3K/AKT signaling pathway in vitro. In conclusion, LTBP2 was a promising diagnostic biomarker of BCR of PCa and had an important role in CD4⁺ T-cell recruitment. Moreover, it was associated with immunotherapy in patients with PCa who developed BCR, and it inhibited PCa proliferation and metastasis via the PI3K/AKT signaling pathway in vitro.

Bioinformatics analysis of KLF2 as a potential prognostic factor in ccRCC and association with epithelial‑mesenchymal transition

Clear cell renal cell carcinoma (ccRCC) is a primary pathological subtype of RCC and has poor clinical outcome. Krüppel-like factors (KLFs), which are zinc-finger proteins, may be involved in ccRCC development and progression. KLFs belong to the zinc-finger family of DNA-binding transcription factors and regulate transcription of downstream target genes. KLFs are involved in cancer development. The present study aimed to investigate the role of KLFs in ccRCC prognosis. The Cancer Genome Atlas database and multifactorial analysis showed that KLFs were widely expressed in pan-cancers and KLF2 was an independent protective factor for ccRCC prognosis. Patients with low KLF2 expression had a low survival probability and expression of KLF2 was downregulated in patients with ccRCC with high pathological grade (II + III vs. I). In addition, western blot and reverse transcription-quantitative PCR revealed that KLF2 was expressed at low levels in ccRCC cell lines and overexpression of KLF2 inhibited cell migration. In addition, KLF2 expression was negatively correlated with methylation. KLF2 expression was elevated following treatment of ccRCC cells with DNA methyltransferase inhibitor. A prognostic risk index prediction model was constructed based on multiple Cox regression. The receiver operating characteristic curve was 0.780 (area under curve >0.5). Furthermore, Gene Ontology enrichment analysis showed that ‘cell adhesion’ and ‘junction’ were negatively correlated with KLF2 and that high-risk group exhibited significantly activated ‘epithelial-mesenchymal transition’. Western blot analysis showed that overexpression of KLF2 increased expression of E-cadherin, while decreasing levels of N-cadherin and vimentin. The present study highlighted the role of KLFs in ccRCC prognosis prediction and provides a research base for the search of validated prognostic biological markers for ccRCC.

Construction of m6A-Related lncRNA Prognostic Signature Model and Immunomodulatory Effect in Glioblastoma Multiforme

Background

Glioblastoma multiforme (GBM), the most prevalent and aggressive of primary malignant central nervous system tumors (grade IV), has a poor clinical prognosis. This study aimed to assess and predict the survival of GBM patients by establishing an m6A-related lncRNA signaling model and to validate its validity, accuracy and applicability.

Methods

RNA sequencing data and clinical data of GBM patients were obtained from TCGA data. First, m6A-associated lncRNAs were screened and lncRNAs associated with overall survival in GBM patients were obtained. Subsequently, the signal model was established using LASSO regression analysis, and its accuracy and validity are further verified. Finally, GO enrichment analysis was performed, and the influence of this signature on the immune regulation response and anticancer drug sensitivity of GBM patients was discussed.

Results

The signature constructed by four lncRNAs AC005229.3, SOX21-AS1, AL133523.1, and AC004847.1 is obtained. Furthermore, the signature proved to be effective and accurate in predicting and assessing the survival of GBM patients and could function independently of other clinical characteristics (Age, Gender and IDH1 mutation). Finally, Immunosuppression-related factors, including APC co-inhibition, T-cell co-inhibition, CCR and Check-point, were found to be significantly up-regulated in GBM patients in the high-risk group. Some chemotherapeutic drugs (Doxorubicin and Methotrexate) and targeted drugs (AZD8055, BI.2536, GW843682X and Vorinostat) were shown to have higher IC50 values in patients in the high-risk group.

Conclusion

We constructed an m6A-associated lncRNA risk model to predict the prognosis of GBM patients and provide new ideas for the treatment of GBM. Further biological experiments can be conducted on this basis to validate the clinical value of the model.

N6-Methyladenosine-Related Long Non-Coding RNAs Are Identified as a Potential Prognostic Biomarker for Lung Squamous Cell Carcinoma and Validated by Real-Time PCR

Currently, the precise mechanism by which N6-methyladenosine (m6A) modification of long non-coding RNAs (lncRNAs) promotes the occurrence and development of lung squamous cell carcinoma (LUSC) and influences tumor microenvironment (TME) remains unclear. Therefore, we studied the prognostic value of m6A-related lncRNAs and their relationship with TME in 495 LUSC samples from The Cancer Genome Atlas (TCGA) database. Pearson’s correlation and univariate Cox regression analysis identified 6 m6A-related lncRNAs with prognostic values for LUSC patients. LUSC patients were divided into two subgroups (clusters 1 and 2) using principal component analysis. The expression of PD-L1 was lower in tumor tissues and cluster 2 of LUSC patients. Cluster 2 of LUSC patients had a high immune score, stromal score, and unique immune cell infiltration. The focal adhesion kinase (FAK) pathway and cytokine receptor pathways are enriched in cluster 1. The m6A-related lncRNA prognostic markers (m6A-LPMs) were established using the least absolute shrinkage and selection operator (LASSO) Cox regression analysis. The risk score was calculated by 4 m6A-LPMs and associated with OS, TME, clinicopathological characteristics of LUSC patients. After adjusting for age, gender, and stage, the risk score was also an independent prognostic factor for LUSC patients. Real-time PCR results showed that the expression of 4 m6A-LPMs was consistent with our prediction results. Our study found that 4 m6A-LPMs (AC138035.1, AC243919.2, HORMAD2-AS1, and AL122125.1) are closely associated with LUSC prognosis, in future, they may as novel diagnostic biomarkers for LUSC and provide new immunotherapy targets for LUSC patients.

Comprehensive Analysis of the Correlation Between Pyroptosis-Related LncRNAs and Tumor Microenvironment, Prognosis, and Immune Infiltration in Hepatocellular Carcinoma

Background: Accumulating evidence shows that pyroptosis plays a crucial role in hepatocellular carcinoma (HCC). However, the relationship between pyroptosis-related long non-coding RNAs (lncRNAs) and HCC tumor characteristics remains enigmatic. We aimed to explore the predictive effect of pyroptosis-related lncRNAs (PRLs) in the prognosis of HCC.

Methods: We comprehensively analyzed the role of the PRLs in the tumor microenvironment and HCC prognosis by integrating genomic data from patients of HCC. Consensus clustering analysis of PRLs was applied to identify HCC subtypes. A prognostic model was then established with a training cohort from The Cancer Genome Atlas (TCGA) using univariate and least absolute shrinkage and selection operator (LASSO) Cox regression analysis. Further, we evaluated the accuracy of this predictive model using a validation set. We predicted IC50s of commonly used chemotherapeutic and targeted drugs through the R package pRRophetic.

Results: Based on pyroptosis-related lncRNAs, a prognostic risk signature composed of seven PRLs (MKLN1AS, AL031985.3, SNHG4, GHRLOS, AC005479.2, AC099850.4, and AC026412.3) was established. For long-term prognosis of HCC patients, our model shows excellent accuracy to forecast overall survival of HCC individuals both in training set and testing set. We found a significant correlation between clinical features and the risk score. Patients in the high-risk group had tumor characteristics associated with progression such as aggressive pathological grade and stage. Besides that, gene set enrichment analysis (GSEA) showed that cell cycle and focal adhesion were significantly enriched in the high-risk group.

Conclusion: The association of the risk model constituted by these seven pyroptosis-related lncRNAs with clinical prognosis, tumor microenvironment, chemotherapy and small molecule drugs was evaluated. Our study provides strong evidence for individualized prediction of prognosis, shedding light on immunotherapy in HCC patients.

Transcriptome-Wide Dynamics of m6A Methylation in Tumor Livers Induced by ALV-J Infection in Chickens

Avian Leukosis Virus Subgroup J (ALV-J) is a tumorigenic virus with high morbidity and rapid transmission. N6-methyladenosine (m⁶A) is a common epigenetic modification that may be closely related to the pathogenicity of ALV-J. Currently, there are no reports on whether m⁶A modification is related to ALV-J induced tumor formation. In this study, we used methylated RNA immunoprecipitation sequencing (MeRIP-seq) and RNA sequencing (RNA-seq) to examine the differences in m⁶A methylation and gene expression in normal livers and ALV-J-induced tumor livers systematically, with functional enrichment and co-expression analysis. The results identified 6,541 m⁶A methylated peaks, mainly enriched in CDS, and more than 83% of the transcripts contained 1-2 m⁶A peaks. For RNA-seq, 1,896 and 1,757 differentially expressed mRNAs and lncRNAs were identified, respectively. Gene enrichment analysis indicated that they may be involved in biological processes and pathways such as immunology-related and apoptosis. Moreover, we identified 17 lncRNAs, commonly existing in differently expressed methylome and transcriptome. Through co-expression analysis, 126 differentially expressed lncRNAs, and 18 potentially m⁶A-related methyltransferases were finally identified and connected, suggesting that m⁶A modifications might affect gene expression of lncRNAs and play a role in ALV-J induced tumor formation. This study provides the first comprehensive description of the m⁶A expression profile in tumor livers induced by ALV-J infection in chickens, which provides a basis for studying the role of m⁶A modification in ALV-J induced tumorigenesis. This study provides clues for studying the epigenetic etiology and pathogenesis of ALV-J.

The Diagnostic and Prognostic Values of HOXA Gene Family in Kidney Clear Cell Renal Cell Carcinoma

Kidney renal clear cell carcinoma (KIRC) is one of the most common cancers with high mortality worldwide. As members of the homeobox (HOX) family, homeobox-A (HOXA) genes have been reported to play an increasingly important role in tumorigenesis and the progression of multiple cancers. However, limited studies have investigated the potential diagnostic and prognostic roles of HOXA genes in KIRC. In this research, we explored the expression pattern of the HOXA gene family in KIRC progression by differential analysis of expression profiles from The Cancer Genome Atlas (TCGA). By using univariate Cox analysis and lasso regression analysis, we comprehensively evaluated the prognostic value of HOXA genes and eventually identified a prognostic risk model consisting of five HOXA genes (HOXA2, HOXA3, HOXA7, HOXA11, and HOXA13). The risk model was further validated as a novel independent prognostic factor for KIRC patients based on the calculated risk score by Kaplan-Meier analysis, univariate and multivariate Cox regression analyses, and time-dependent receiver operating characteristic (ROC) curve analysis. Moreover, to explore the potential mechanism of tumorigenesis and clinical application of KIRC, we also developed the HOXA-based competing endogenous RNA (ceRNA) regulatory network and machine learning classification model. Valproic acid and tretinoin were predicted to be the most promising small molecules to adjuvant treatment of KIRC by mining the CMAP and DGIdb drug database. Subsequently, pathway and functional enrichment analyses provided us with new ways to search for a possible mechanism of action of drugs. Taken together, our study demonstrated the nonnegligible role of HOXA genes in KIRC and constructed an effective prognostic and diagnostic model, which offers novel insights into KIRC prognosis.

Pyroptosis-Related lncRNAs Predict the Prognosis and Immune Response in Patients With Breast Cancer

Background: Breast cancer (BC) is the most common malignant tumor and the leading cause of cancer-related death in women worldwide. Pyroptosis and long noncoding RNAs (lncRNAs) have been demonstrated to play vital roles in the tumorigenesis and development of BC. However, the clinical significance of pyroptosis-related lncRNAs in BC remains unclear.

Methods: Using the mRNA and lncRNA profiles of BC obtained from TCGA dataset, a risk model based on the pyroptosis-related lncRNAs for prognosis was constructed using univariate and multivariate Cox regression model, and least absolute shrinkage and selection operator. Patients were divided into high- and low-risk groups based on the risk model, and the prognosis value and immune response in different risk groups were analyzed. Furthermore, functional enrichment annotation, therapeutic signature, and tumor mutation burden were performed to evaluate the risk model we established. Moreover, the expression level and clinical significance of the selected pyroptosis-related lncRNAs were further validated in BC samples.

Results: 3,364 pyroptosis-related lncRNAs were identified using Pearson’s correlation analysis. The risk model we constructed comprised 10 pyroptosis-related lncRNAs, which was identified as an independent predictor of overall survival (OS) in BC. The nomogram we constructed based on the clinicopathologic features and risk model yielded favorable performance for prognosis prediction in BC. In terms of immune response and mutation status, patients in the low-risk group had a higher expression of immune checkpoint markers and exhibited higher fractions of activated immune cells, while the high-risk group had a highly percentage of TMB. Further analyses in our cohort BC samples found that RP11-459E5.1 was significantly upregulated, while RP11-1070N10.3 and RP11-817J15.3 were downregulated and significantly associated with worse OS.

Conclusion: The risk model based on the pyroptosis-related lncRNAs we established may be a promising tool for predicting the prognosis and personalized therapeutic response in BC patients.

Identification of a Hypoxia-Related Gene Model for Predicting the Prognosis and Formulating the Treatment Strategies in Kidney Renal Clear Cell Carcinoma

Purpose

The present study aimed to establish a hypoxia related genes model to predict the prognosis of kidney clear cell carcinoma (KIRC) patients using data accessed from The Cancer Genome Atlas (TCGA) database and International Cancer Genome Consortium (ICGC) database.

Methods

Patients’ data were downloaded from the TCGA and ICGC databases, and hypoxia related genes were accessed from the Molecular Signatures Database. The differentially expressed genes were evaluated and then the differential expressions hypoxia genes were screened. The TCGA cohort was randomly divided into a discovery TCGA cohort and a validation TCGA cohort. The discovery TCGA cohort was used for constructing the hypoxia genes risk model through Lasso regression, univariate and multivariate Cox regression analysis. Receiver operating characteristic (ROC) curves were used to assess the reliability and sensitivity of our model. Then, we established a nomogram to predict the probable one-, three-, and five-year overall survival rates. Lastly, the Tumor Immune Dysfunction and Exclusion (TIDE) score of patients was calculated.

Results

We established a six hypoxia-related gene prognostic model of KIRC patients in the TCGA database and validated in the ICGC database. The patients with high riskscore present poorer prognosis than those with low riskscore in the three TCGA cohorts and ICGC cohort. ROC curves show our six-gene model with a robust predictive capability in these four cohorts. In addition, we constructed a nomogram for KIRC patients in the TCGA database. Finally, the high risk-group had a high TIDE score than the patients with low riskscore.

Conclusions

We established a six hypoxia-related gene risk model for independent prediction of the prognosis of KIRC patients was established and constructed a robust nomogram. The different riskscores might be a biomarker for immunotherapy strategy.

Characterization of an Autophagy-Immune Related Genes Score Signature and Prognostic Model and its Correlation with Immune Response for Bladder Cancer

Purpose

The study aimed to identify an autophagy-related molecular subtype and characterize a novel defined autophagy-immune related genes score (AI-score) signature and prognosis model in bladder cancer (BLCA) patients using public databases.

Methods

The transcriptome cohorts downloaded from TCGA and GEO database were carried out with genomic analysis and unsupervised methods to obtain autophagy-related molecular subtypes. The single-sample gene-set enrichment analysis (ssGSEA) was utilized to perform immune subtype clustering. We defined a novel autophagy subtype and evaluated the role in TME cell infiltration. Then, the principal-component analysis (PCA) was applied to construct an AI-score signature. Subsequently, two immunotherapeutic cohorts were used to evaluate the predictive value in immunotherapeutic benefits and immune response. Finally, univariate, Lasso and multivariate Cox regression algorithm were used to construct and evaluate an autophagy-immune-related genes prognosis model. Also, qRT-PCR and IHC was applied to validate the expression of the 6 genes in the model.

Results

Three distinct autophagy clusters and immune-related clusters were identified, and a novel autophagy-related molecular subtypes were defined. Furthermore, the roles in TME cell infiltration and clinical traits for the autophagy subtypes were characterized. Meanwhile, we constructed an AI-score signature and demonstrated it could predict genetic mutation, clinicopathological traits, prognosis, and TME stromal activity. We found that it could accurately predict the clinicopathological characteristics and immune response of individual BLCA patients and provide guidance for selecting immunotherapy. Ultimately, we constructed and verified an autophagy-immune-related prognostic model of BLCA patients and established a prognostic nomogram with a good prediction accuracy.

Conclusion

We constructed AI-score signatures and prognosis risk model to characterize their role in clinical features and TME immune cell infiltration. It revealed that the AI-score signature and prognosis model could be a valid predictive tool, which could accurately predict the prognosis of BLCA patients and contribute to choosing effective personalized immunotherapy strategies.

Potential Prognostic Value of a Seven m6A-Related LncRNAs Signature and the Correlative Immune Infiltration in Colon Adenocarcinoma

Long non-coding RNAs (lncRNAs) and their N6-methyladenosine (m6A) modifications play an essential role in tumorigenesis and cancer progression. This study was designed to explore the value of m6A-related lncRNAs in prognosis and therapeutic applications of immune infiltration of colon adenocarcinoma (COAD). We downloaded the COAD gene expression and clinical data from The Cancer Genome Atlas project. By co-expression analysis, Lasso Cox regression analysis, and univariate and multivariate Cox regression, we constructed an independent prognostic signature of seven m6A-related lncRNAs. The prognostic lncRNAs were divided into two clusters by consistent clustering analysis, as well as into two groups of low–high risk based on the signature. Then we identified the relationship between the different groups with clinical features and immune cell infiltration. Cluster 2 had a higher risk score with a lower survival rate. The risk score was higher in groups with advanced clinical features, such as stage III–IV, N1-3, and M1. The expression of AC156455.1 was increased in tumor tissues and cluster 2, and the lncRNA ZEB1−AS1 was notably higher in the high-risk group. Five types of immune cells showed differences in two clusters, and most were upregulated in type 2. The expression of memory B cells was positively correlated with the risk score. The prognostic model was verified by the Gene Expression Omnibus (GEO) dataset. Besides, we found that the expression of these seven lncRNAs in tumor tissues was significantly higher than that in normal tissues, which verified the feasibility of the model. Thus, the signature of seven m6A-related lncRNAs can independently predict the prognosis of COAD. This signature is also closely associated with immune cell infiltration, and new therapeutic targets can be explored from this field.

Non-Coding RNA m6A Modification in Cancer: Mechanisms and Therapeutic Targets

Recently, N6-methyl-adenosine (m6A) ribonucleic acid (RNA) modification, a critical and common internal RNA modification in higher eukaryotes, has generated considerable research interests. Extensive studies have revealed that non-coding RNA m6A modifications (e.g. microRNAs, long non-coding RNAs, and circular RNAs) are associated with tumorigenesis, metastasis, and other tumour characteristics; in addition, they are crucial molecular regulators of cancer progression. In this review, we discuss the relationship between non-coding RNA m6A modification and cancer progression from the perspective of various cancers. In particular, we focus on important mechanisms in tumour progression such as proliferation, apoptosis, invasion and metastasis, tumour angiogenesis. In addition, we introduce clinical applications to illustrate more vividly that non-coding RNA m6A modification has broad research prospects. With this review, we aim to summarize the latest insights and ideas into non-coding RNA m6A modification in cancer progression and targeted therapy, facilitating further research.

Identification of a N6-Methyladenosine (m6A)-Related lncRNA Signature for Predicting the Prognosis and Immune Landscape of Lung Squamous Cell Carcinoma

Background

m6A-related lncRNAs emerged as potential targets for tumor diagnosis and treatment. This study aimed to identify m6A-regulated lncRNAs in lung squamous cell carcinoma (LUSC) patients.

Materials and Methods

RNA sequencing and the clinical data of LUSC patients were downloaded from The Cancer Genome Atlas (TCGA) database. The m6A-related lncRNAs were identified by using Pearson correlation assay. Univariate and multivariate Cox regression analyses were utilized to construct a risk model. The performance of the risk model was validated using Kaplan–Meier survival analysis and receiver operating characteristics (ROC). Immune estimation of LUSC was downloaded from TIMER, and the correlations between the risk score and various immune cells infiltration were analyzed using various methods. Differences in immune functions and expression of immune checkpoint inhibitors and m6A regulators between high-risk and low-risk groups were further explored. Finally, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were utilized to explore the biological functions of AL122125.1.

Results

A total of 351 m6A-related lncRNAs were obtained from TCGA. Seven lncRNAs demonstrated prognostic values. A further multivariate Cox regression assay constructed a risk model consisting of two lncRNAs (AL122125.1 and HORMAD2-AS1). The Kaplan–Meier analysis and area under the curve indicated that this risk model could be used to predict the prognosis of LUSC patients. The m6A-related lncRNAs were immune-associated. There were significant correlations between risk score and immune cell infiltration, immune functions, and expression of immune checkpoint inhibitors. Meanwhile, there were significant differences in the expression of m6A regulators between the high- and low-risk groups. Moreover, GO and KEGG analyses revealed that the upregulated expression of AL122125.1 was tumor-related.

Conclusion

In this study, we constructed an m6A-related lncRNA risk model to predict the survival of LUSC patients. This study could provide a novel insight to the prognosis and treatment of LUSC patients.

Angiogenesis Pathway in Kidney Renal Clear Cell Carcinoma and Its Prognostic Value for Cancer Risk Prediction

Angiogenesis, a process highly regulated by pro-angiogenic and anti-angiogenic factors, is disrupted and dysregulated in cancer. Despite the increased clinical use of angiogenesis inhibitors in cancer therapy, most molecularly targeted drugs have been less effective than expected. Therefore, an in-depth exploration of the angiogenesis pathway is warranted. In this study, the expression of angiogenesis-related genes in various cancers was explored using The Cancer Genome Atlas datasets, whereupon it was found that most of them were protective genes in the patients with kidney renal clear cell carcinoma (KIRC). We divided the samples from the KIRC dataset into three clusters according to the mRNA expression levels of these genes, with the enrichment scores being in the order of Cluster 2 (upregulated expression) > Cluster 3 (normal expression) > Cluster 1 (downregulated expression). The survival curves plotted for the three clusters revealed that the patients in Cluster 2 had the highest overall survival rates. Via a sensitivity analysis of the drugs listed on the Genomics of Drug Sensitivity in Cancer database, we generated IC50 estimates for 12 commonly used molecularly targeted drugs for KIRC in the three clusters, which can provide a more personalized treatment plan for the patients according to angiogenesis-related gene expression. Subsequently, we investigated the correlation between the angiogenesis pathway and classical cancer-related genes as well as that between the angiogenesis score and immune cell infiltration. Finally, we used the least absolute shrinkage and selection operator (LASSO)-Cox regression analysis to construct a risk score model for predicting the survival of patients with KIRC. According to the areas under the receiver operating characteristic (ROC) curves, this new survival model based on the angiogenesis-related genes had high prognostic prediction value. Our results should provide new avenues for the clinical diagnosis and treatment of patients with KIRC.

Identification of Hub Genes Associated With Clear Cell Renal Cell Carcinoma by Integrated Bioinformatics Analysis

BACKGROUND

Clear cell renal cell carcinoma (ccRCC) is a common genitourinary cancer type with a high mortality rate. Due to a diverse range of biochemical alterations and a high level of tumor heterogeneity, it is crucial to select highly validated prognostic biomarkers to be able to identify subtypes of ccRCC early and apply precision medicine approaches.

METHODS

Transcriptome data of ccRCC and clinical traits of patients were obtained from the GSE126964 dataset of Gene Expression Omnibus and The Cancer Genome Atlas Kidney Renal Clear Cell Carcinoma (TCGA-KIRC) database. Weighted gene co-expression network analysis (WGCNA) and differentially expressed gene (DEG) screening were applied to detect common differentially co-expressed genes. Gene Ontology, Kyoto Encyclopedia of Genes and Genomes analysis, survival analysis, prognostic model establishment, and gene set enrichment analysis were also performed. Immunohistochemical analysis results of the expression levels of prognostic genes were obtained from The Human Protein Atlas. Single-gene RNA sequencing data were obtained from the GSE131685 and GSE171306 datasets.

RESULTS

In the present study, a total of 2,492 DEGs identified between ccRCC and healthy controls were filtered, revealing 1,300 upregulated genes and 1,192 downregulated genes. Using WGCNA, the turquoise module was identified to be closely associated with ccRCC. Hub genes were identified using the maximal clique centrality algorithm. After having intersected the hub genes and the DEGs in GSE126964 and TCGA-KIRC dataset, and after performing univariate, least absolute shrinkage and selection operator, and multivariate Cox regression analyses, ALDOB, EFHD1, and ESRRG were identified as significant prognostic factors in patients diagnosed with ccRCC. Single-gene RNA sequencing analysis revealed the expression profile of ALDOB, EFHD1, and ESRRG in different cell types of ccRCC.

CONCLUSIONS

The present results demonstrated that ALDOB, EFHD1, and ESRRG may act as potential targets for medical therapy and could serve as diagnostic biomarkers for ccRCC.

Gene Characteristics and Prognostic Values of m6A RNA Methylation Regulators in Nonsmall Cell Lung Cancer

Background

N⁶-methyladenosine (m⁶A) is the most common internal modification present in mRNAs and long noncoding RNAs (lncRNAs), associated with tumorigenesis and cancer progression. However, little is known about the roles of m⁶A and its regulatory genes in nonsmall cell lung cancer (NSCLC). Here, we systematically explored the roles and prognostic significance of m⁶A-associated regulatory genes in NSCLC.

Methods

The copy number variation (CNV), mutation, mRNA expression data, and corresponding clinical pathology information of 1057 NSCLC patients were downloaded from the cancer genome atlas (TCGA) database. The gain and loss levels of CNVs were determined by utilizing segmentation analysis and GISTIC algorithm. The GSEA was conducted to explore the functions related to different levels of m⁶A regulatory genes. Logrank test was utilized to assess the prognostic significance of m⁶A-related gene's CNV.

Results

The genetic alterations of ten m⁶A-associated regulators were identified in 102 independent NSCLC samples and significantly related to advanced tumor stage. Deletions or shallow deletions corresponded to lower mRNA expression while copy number gains or amplifications were related to increased mRNA expression of m⁶A regulatory genes. Survival analysis showed the patients with copy number loss of FTO with worse disease-free survival (DFS) or overall survival (OS). Besides, copy number loss of YTHDC2 was also with poor OS for NSCLC patients. Moreover, high FTO expression was significantly associated with oxidative phosphorylation, translation, and metabolism of mRNA.

Conclusion

Our findings provide novel insight for better understanding of the roles of m⁶A regulators and RNA epigenetic modification in the pathogenesis of NSCLC.