On representing the prognostic value of continuous gene expression biomarkers with the restricted mean survival curve

Multimodal radiopathological integration for prognosis and prediction of adjuvant chemotherapy benefit in resectable lung adenocarcinoma: A multicentre study

Lung adenocarcinoma (LUAD) has a heterogeneous prognosis and controversial postoperative treatment protocols. We aim to develop and validate a multimodal analysis framework that integrates CT images with H&E-stained whole-slide images (WSIs) to enhance risk stratification and predict adjuvant chemotherapy benefit in LUAD patients. We retrospectively collected data from 1039 resectable LUAD patients (stage I-III) across four centres, forming a training dataset (n = 303), two testing datasets (n = 197 and n = 228) for survival analysis, and a feature testing dataset (n = 311) for interpretability analysis. We extracted 487 tumour/peritumour radiomics features from CT images and 783 multiscale pathomics features from WSIs, characterising the shape of tumour (CT) and cancer nuclei (WSIs), as well as the intensity and texture of tumour/peritumour regions (CT) and tumour regions/epithelium/stroma (WSIs). A survival support vector machine (SVM) was employed to establish a radiopathomics signature using the optimal set of multimodal features, including 2 tumour radiomics features, 3 peritumour radiomics features, and 4 nuclei heterogeneity pathomics features. The radiopathomics signature outperformed both radiomics and pathomics signatures in predicting disease-free survival (DFS) (C-index: training dataset, 0.744 vs. 0.734 and 0.692; testing dataset 1, 0.719 vs. 0.701 and 0.638; testing dataset 2, 0.711 vs. 0.689 and 0.684), demonstrating greater robustness compared to the state-of-the-art deep learning integration approaches. It provided additional prognostic information beyond clinical risk factors (C-index of clinical plus radiopathomics vs. clinical models: training dataset, 0.763 vs. 0.676; testing dataset 1, 0.739 vs. 0.676; testing dataset 2, 0.711 vs. 0.699, p < 0.001). Compared to low-risk patients categorised by the radiopathomics signature, high-risk patients achieved comparable DFS when receiving adjuvant chemotherapy (training dataset, HR = 1.53, 95 % CI 0.85-2.73, p = 0.153; testing dataset 1 and 2, HR = 1.62, 95 % CI 0.92-2.85, p = 0.096), but had significantly worse DFS when only observed after surgery (training dataset, HR = 4.46, 95 % CI 2.82-7.05, p < 0.001; testing datasets 1 and 2, HR = 3.52, 95 % CI 2.26-5.49, p < 0.001), indicating the predictive value of the radiopathomics signature for adjuvant chemotherapy benefit (interaction p < 0.05). Further interpretability analysis revealed that the radiopathomics signature was associated with various prognostic/treatment-related biomarkers, including differentiation, immune phenotypes, and EGFR status. The multimodal integration framework offered a cost-effective approach for LUAD characterisation by leveraging complementary information from radiological and histopathological imaging. The radiopathomics signature demonstrated robust prognostic capabilities, providing valuable insights for postoperative treatment decisions.

STEAP4 facilitates growth, migration, and invasion of prostate carcinoma through upregulation of NOTCH4

STEAP4 manifested differential expression and aberrant methylation in prostate cancer (PCa). Therefore, this study proposed to explore the effect of STEAP4 on the PCa malignant phenotype in vivo and in vitro and the possible molecular mechanisms using RNA-seq. The expression of STEAP4 in PCa and its prognostic and diagnostic value was identified using bioinformatics. After exogenous modulation of STEAP4, the effect of STEAP4 on the malignant phenotype of PCa cells was examined using functional assays and nude mouse tumor models. The STEAP4-related differentially expressed genes (DEGs) and the hub genes were characterized using RNA-seq in conjunction with bioinformatics. STEAP4 exhibited high expression in PCa tissues from TCGA-PRAD and GEO datasets (GSE179321, GSE229904, and GSE237995), which predicted lower survival of patients. The STEAP4-associated nomogram model and diagnostic ROC curve had excellent predictive performance (AUC = 0.814). STEAP4 was overexpressed in PCa tissues and cells. Knockdown of STEAP4 effectively decreased the viability, number of invading cells, and wound healing of PCa cells and increased apoptosis. Overexpression of STEAP4 showed the opposite pattern. RNA-seq revealed that knockdown of STEAP4 resulted in 234 DEGs in PCa cells. FGF17, KCNQ2, PDGFRB, and NOTCH4 are hub genes in DEGs. Notably, NOTCH4 was likewise overexpressed in PCa tissues and cells and was regulated by STEAP4. In in vitro experiments, overexpression of NOTCH4 facilitated PCa cell proliferation, migration, and invasion, which was limited by knockdown of STEAP4. In in vivo experiments, overexpression of STEAP4 exacerbated PCa tumor burden, which was rescued by knockdown of NOTCH4. STEAP4 is a valid biomarker for predicting prognosis and diagnosis of PCa patients. STEAP4 contributes to PCa growth, migration, and invasion by upregulating NOTCH4.

Putative function and prognostic molecular marker of mast cells in colorectal cancer

BACKGROUND

The increased demand for markers for colorectal cancer (CRC) highlights the importance of investigating immune cells involved in CRC progression. This study aims to dissect the mast cells in CRC, characterize the role of mast cells in CRC development, coordinate molecular communication between mast cells and malignant cells, and construct and validate a prognostic classification model based on mast cell markers.

METHODS

Single-cell transcriptome data of CRC patients were extracted from GSE146771 for cell classification and annotation. The malignant cells were identified by copykat and the communication between mast cells and malignant cells was analyzed by CellChat. Least absolute shrinkage and selection operator (LASSO) regression analysis and Cox regression analysis of mast cell markers were performed in the TCGA-COAD cohort to construct a prognostic classification model. qRT-PCR was performed to detect the mRNA expression of the molecules in the classification model in P815 and MC-9 cells. The co-culture experiment of MC38 and P815 cells were performed in 12-well transwell dish. Wound healing assay and Transwell assay were performed to detect cell migration and invasion.

RESULTS

10,186 high-quality cells in GSE146771 were annotated to 9 cell types. Six markers in mast cells (HDC, GATA2, ASAH1, BTBD19, TIMP1, FAM110A) were selected to construct a classification model. The high-risk score defined showed high infiltration of immunosuppressive cells, including endothelial cells, CAFs, Tregs and high angiogenesis and epithelial-mesenchymal transition (EMT) activities. In the model, HDC were abnormally low expressed in P815 cells, while BTBD19, FAM110A, GATA2, ASAH1 and TIMP1 showed excessive expression in P815 cells. Knockdown of GATA2 in the co-culture system of P815 and MC38 cells blocked cell migration and invasion.

CONCLUSION

This study identified the cell types within CRC, elaborated the cellular functions of mast cells in CRC development and their molecular communication to coordinate malignant cells, and highlighted the molecular components and biological features that constitute promising prognostic classification model.

Construction and evaluation of a prognostic model based on the expression of the metabolism-related signatures in patients with osteosarcoma

BACKGROUND

The aim of this study was to screen three major substance metabolism-related genes and establish a prognostic model for osteosarcoma.

METHODS

RNA-seq expression data for osteosarcoma were downloaded from The Cancer Genome Atlas (TCGA) and GEO databases. Differentially expressed (DE) RNAs were selected, followed by the selection of metabolic-related DE mRNAs. Using Cox regression analysis, prognostic DE RNAs were identified to construct a prognostic model. Subsequently, independent prognostic clinical factors were screened, and the functions of the long non-coding RNAs (lncRNAs) were analyzed. Finally, the expression of signature genes was further tested in osteosarcoma cells using quantitative reverse transcription quantitative real-time polymerase chain reaction (qRT-PCR) and western blotting.

RESULTS

A total of 432 DE RNAs, comprising 79 DE lncRNAs and 353 DE mRNAs were obtained, and then 107 metabolic-related DE mRNAs. Afterwards signature genes (LINC00545, LINC01537, FOXC2-AS1, CYP27B1, PFKFB4, PHKG1, PHYKPL, PXMP2, and XYLB) served as optimal combinations, and a prognostic score model was successfully proposed. Three verification datasets (GSE16091, GSE21257, and GSE39055) showed that the model had high specificity and sensitivity. In addition, two independent prognostic clinical factors (age and tumor metastasis) were identified. Finally, the concordance rate between the in silico analysis, qRT-PCR, and western blotting analysis was 88.89% (8/9), suggesting the robustness of our analysis.

CONCLUSIONS

The prognostic model based on the nine signature genes accurately predicted the prognosis of patients with osteosarcoma; CYP27B1, PFKFB4, PHKG1, PHYKPL, PXMP2, and XYLB may serve as metabolism-related biomarkers in osteosarcoma.

NDRG1 alleviates Erastin-induced ferroptosis of hepatocellular carcinoma

BACKGROUND

NDRG1, a cell differentiation-associated factor, has recently emerged as a regulator ferroptosis. Nevertheless, its role in modulating ferroptosis within hepatocellular carcinoma (HCC) remains uncharacterized.

METHODS

The differential expression of NDRG1 and its prognostic value were analyzed in HCC using data from TCGA and GEO. Ferroptosis in HepG2 and Huh7 cells was assessed using flow cytometry, transmission electron microscopy, and propidium iodide staining following NDRG1 knockdown using shRNA. RNA-seq was performed to characterize the mRNA expression profiles in HepG2 cells, identifying differentially expressed mRNAs (DE-mRNAs) and NDRG1-related hub genes.

RESULTS

NDRG1 was overexpressed in multiple malignant tumors, including HCC, and was associated with a significantly poor prognosis in HCC patients. A nomogram model integrating NDRG1 expression and clinical parameters demonstrated robust prognostic accuracy. NDRG1 knockdown potentiated erastin-induced alterations in Fe2+, total ROS, lipid ROS, and ferroptosis markers (PTGS2, ACSL4, GPX4, SLC7A11, GSH, GSSG), while exacerbating mitochondrial ultrastructural damage in HepG2 and Huh7 cells. Erastin induction elicited 1,056 DE-mRNAs, while subsequent NDRG1 knockdown revealed 1,323 DE-mRNAs in HepG2 cells. These DE-mRNAs are mainly involved in metastasis, immunity, growth, ferroptosis, and are associated with AMPK, MAPK, and PI3K/AKT pathways. Moreover, NDRG1 potentially interacted with HSPA8, CDH1, ALDOC, ANGPTL4, ANKRD37, CA9, ERBB3, FOS. qRT-PCR confirmed their expression changes consistent with RNA-seq.

CONCLUSION

NDRG1 exhibits strong predictive value for HCC, and accelerates tumor progression by suppressing ferroptosis.

ATRX mutation modifies the DNA damage response in glioblastoma multiforme tumor cells and enhances patient prognosis

The presence of specific genetic mutations in patients with glioblastoma multiforme (GBM) is associated with improved survival outcomes. Disruption of the DNA damage response (DDR) pathway in tumor cells enhances the effectiveness of radiotherapy drugs, while increased mutational burden following tumor cell damage also facilitates the efficacy of immunotherapy. The ATRX gene, located on chromosome X, plays a crucial role in DDR. The aim of this research is to elucidate the correlation between ATRX mutations and GBM. Dataset obtained from TCGA-GBM were conducted an analysis on the genomic features, biological characteristics, immunopathological markers, and clinical prognosis of patients carrying ATRX mutations. Our findings revealed a significantly elevated level of microsatellite instability in individuals with ATRX mutants, along with significant alterations in the receptor-tyrosine kinase (RTK)-ras pathway among patients exhibiting combined ATRX mutations. TCGA-GBM patients with concurrent ATRX mutations exhibited sensitivity to 26 chemotherapeutic and anticancer drugs, which exerted their effects by modulating the DDR of tumor cells through highly correlated mechanisms involving the RTK-ras pathway. Additionally, we observed an enrichment of ATRX mutations in specific pathways associated with DDR among TCGA-GBM patients. Our model also demonstrated prolonged overall survival in patients carrying ATRX mutations, particularly showing strong predictive value for 3- and 5-year survival rates. Furthermore, additional protective factors such as younger age, female gender, combined IDH mutations, and TP53 mutations were identified. The results underscore the protective role and prognostic significance of ATRX mutations in GBM as a potential therapeutic target and biomarker for patient survival.

LOCC: a novel visualization and scoring of cutoffs for continuous variables with hepatocellular carcinoma prognosis as an example

BACKGROUND

The interpretation of large datasets, such as The Cancer Genome Atlas (TCGA), for scientific and research purposes, remains challenging despite their public availability. In this study, we focused on identifying gene expression profiles most relevant to patient prognosis and aimed to develop a method and database to address this issue. To achieve this, we introduced Luo's Optimization Categorization Curve (LOCC), an innovative tool for visualizing and scoring continuous variables against dichotomous outcomes. To demonstrate the efficacy of LOCC using real-world data, we analyzed gene expression profiles and patient data from TCGA hepatocellular carcinoma samples.

RESULTS

To showcase LOCC, we demonstrate an optimal cutoff for E2F1 expression in hepatocellular carcinoma, which was subsequently validated in an independent cohort. Compared to ROC curves and their AUC, LOCC offered a superior description of the predictive value of E2F1 expression across various cancer types. The LOCC score, comprised of factors representing significance, range, and impact of the biomarker, facilitated the ranking of all gene expression profiles in hepatocellular carcinoma, aiding in the evaluation and understanding of previously published prognostic gene signatures. We also demonstrate that LOCC does not have the same assumptions required of Cox proportional hazards modeling for accurate analysis. Repeated sampling demonstrated that LOCC scores outperformed ROC's AUC in discriminating predictors from non-predictors. Additionally, gene set enrichment analysis revealed significant associations between certain genes and prognosis, such as E2F target genes and G2M checkpoint with poor prognosis, and bile acid metabolism and oxidative phosphorylation with good prognosis.

CONCLUSION

In summary, we present LOCC as a novel visualization tool for the analysis of gene expression in cancer, particularly for understanding and selecting cutoffs. Our findings suggest that LOCC scores, which effectively rank genes based on their prognostic potential, represent a more suitable approach than ROC curves and Cox proportional hazard for prognostic modeling and understanding in cancer gene expression analysis. LOCC holds promise as an invaluable tool for advancing precision medicine and furthering biomarker research. Further research regarding multivariable integration and validation will help LOCC reach its full potential and establish its utility across diverse cancer types and clinical settings.

Characterization of RNA Processing Genes in Colon Cancer for Predicting Clinical Outcomes

Objective

Colon cancer is associated with multiple levels of molecular heterogeneity. RNA processing converts primary transcriptional RNA to mature RNA, which drives tumourigenesis and its maintenance. The characterisation of RNA processing genes in colon cancer urgently needs to be elucidated.

Methods

In this study, we obtained 1033 relevant samples from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases to explore the heterogeneity of RNA processing phenotypes in colon cancer. Firstly, Unsupervised hierarchical cluster analysis detected 4 subtypes with specific clinical outcomes and biological features via analysis of 485 RNA processing genes. Next, we adopted the least absolute shrinkage and selection operator (LASSO) as well as Cox regression model with penalty to characterise RNA processing-related prognostic features.

Results

An RNA processing-related prognostic risk model based on 10 genes including FXR1, MFAP1, RBM17, SAGE1, SNRPA1, SRRM4, ADAD1, DDX52, ERI1, and EXOSC7 was identified finally. A composite prognostic nomogram was constructed by combining this feature with the remaining clinical variables including TNM, age, sex, and stage. Genetic variation, pathway activation, and immune heterogeneity with risk signatures were also analysed via bioinformatics methods. The outcomes indicated that the high-risk subgroup was associated with higher genomic instability, increased proliferative and cycle characteristics, decreased tumour killer CD8+ T cells and poorer clinical prognosis than the low-risk group.

Conclusion

This prognostic classifier based on RNA-edited genes facilitates stratification of colon cancer into specific subgroups according to TNM and clinical outcomes, genetic variation, pathway activation, and immune heterogeneity. It can be used for diagnosis, classification and targeted treatment strategies comparable to current standards in precision medicine. It provides a rationale for elucidation of the role of RNA editing genes and their clinical significance in colon cancer as prognostic markers.

A ferroptosis-related signature predicts the clinical diagnosis and prognosis, and associates with the immune microenvironment of lung cancer

Targeting ferroptosis-related pathway is a potential strategy for treatment of lung cancer (LC). Consequently, exploration of ferroptosis-related markers is important for treating LC. We collected LC clinical data and mRNA expression profiles from TCGA and GEO database. Ferroptosis-related genes (FRGs) were obtained through FerrDB database. Expression analysis was performed to obtain differentially expressed FRGs. Diagnostic and prognostic models were constructed based on FRGs by LASSO regression, univariate, and multivariate Cox regression analysis, respectively. External verification cohorts GSE72094 and GSE157011 were used for validation. The interrelationship between prognostic risk scores based on FRGs and the tumor immune microenvironment was analyzed. Immunocytochemistry, Western blotting, and RT-qPCR detected the FRGs level. Eighteen FRGs were used for diagnostic models, 8 FRGs were used for prognostic models. The diagnostic model distinguished well between LC and normal samples in training and validation cohorts of TCGA. The prognostic models for TCGA, GSE72094, and GSE157011 cohorts significantly confirmed lower overall survival (OS) in high-risk group, which demonstrated excellent predictive properties of the survival model. Multivariate Cox regression analysis further confirmed risk score was an independent risk factor related with OS. Immunoassays revealed that in high-risk group, a significantly higher proportion of Macrophages_M0, Neutrophils, resting Natural killer cells and activated Mast cells and the level of B7H3, CD112, CD155, B7H5, and ICOSL were increased. In conclusion, diagnostic and prognostic models provided superior diagnostic and predictive power for LC and revealed a potential link between ferroptosis and TIME.

Oncogene-induced matrix reorganization controls CD8+ T cell function in the soft-tissue sarcoma microenvironment

CD8+ T cell dysfunction impedes antitumor immunity in solid cancers, but the underlying mechanisms are diverse and poorly understood. Extracellular matrix (ECM) composition has been linked to impaired T cell migration and enhanced tumor progression; however, impacts of individual ECM molecules on T cell function in the tumor microenvironment (TME) are only beginning to be elucidated. Upstream regulators of aberrant ECM deposition and organization in solid tumors are equally ill-defined. Therefore, we investigated how ECM composition modulates CD8+ T cell function in undifferentiated pleomorphic sarcoma (UPS), an immunologically active desmoplastic tumor. Using an autochthonous murine model of UPS and data from multiple human patient cohorts, we discovered a multifaceted mechanism wherein the transcriptional coactivator YAP1 promotes collagen VI (COLVI) deposition in the UPS TME. In turn, COLVI induces CD8+ T cell dysfunction and immune evasion by remodeling fibrillar collagen and inhibiting T cell autophagic flux. Unexpectedly, collagen I (COLI) opposed COLVI in this setting, promoting CD8+ T cell function and acting as a tumor suppressor. Thus, CD8+ T cell responses in sarcoma depend on oncogene-mediated ECM composition and remodeling.

RUVBL1 accelerates tongue squamous cell carcinoma by mediating CRaf/MEK/ERK pathway

RAF/MEK/ERK pathway is frequently activated in tumor. Therefore, this study will investigate the function of RUVBL1 (RAF-binding protein) in tongue squamous cell carcinoma (TSCC). Bioinformatics was performed to identify differentially expressed mRNAs (DE-mRNAs) in TCGA-oral squamous cell carcinoma, GSE13601, and GSE34105 datasets. A total of 672 shared DE-mRNAs were identified in three datasets, and they are regulating metastasis and angiogenesis. Patients with RUVBL1 low expression had high overall survival. Overexpressing RUVBL1 enhanced the viability, wound healing percentage, invasion, sphere formation, angiogenesis, and resistance to cisplatin and 5-fluorouracil in CAL-27 and SCC-4 cells, and the opposite results were obtained by knocking down RUVBL1. Moreover, overexpression of RUVBL1 bolstered tumor growth in vivo. Strikingly, RUVBL1 diminished the phosphorylation of CRAF Ser259, which led to activation of the MEK/ERK pathway. In conclusion, RUVBL1 contributes to the malignant biological behavior of TSCC via activating the CRAF/MEK/ERK pathway. This provides molecular mechanisms and perspectives for targeted therapy of the CRAF/MEK/ERK pathway.

Blunted Neurobehavioral Loss Anticipation Predicts Relapse to Stimulant Drug Use

BACKGROUND

Patients with stimulant use disorder experience high rates of relapse. While neurobehavioral mechanisms involved in initiating drug use have been studied extensively, less research has focused on relapse.

METHODS

To assess motivational processes involved in relapse and diagnosis, we acquired functional magnetic resonance imaging responses to nondrug (monetary) gains and losses in detoxified patients with stimulant use disorder (n = 68) and community control participants (n = 42). In a prospective multimodal design, we combined imaging of brain function, brain structure, and behavior to longitudinally track subsequent risk for relapse.

RESULTS

At the 6-month follow-up assessment, 27 patients remained abstinent, but 33 had relapsed. Patients with blunted anterior insula (AIns) activity during loss anticipation were more likely to relapse, an association that remained robust after controlling for potential confounds (i.e., craving, negative mood, years of use, age, and gender). Lower AIns activity during loss anticipation was associated with lower self-reported negative arousal to loss cues and slower behavioral responses to avoid losses, which also independently predicted relapse. Furthermore, AIns activity during loss anticipation was associated with the structural coherence of a tract connecting the AIns and the nucleus accumbens, as was functional connectivity between the AIns and nucleus accumbens during loss processing. However, these neurobehavioral responses did not differ between patients and control participants.

CONCLUSIONS

Taken together, the results of the current study show that neurobehavioral markers predicted relapse above and beyond conventional self-report measures, with a cross-validated accuracy of 72.7%. These findings offer convergent multimodal evidence that implicates blunted avoidance motivation in relapse to stimulant use and may therefore guide interventions targeting individuals who are most vulnerable to relapse.

Immune Cell Infiltration Types as Biomarkers for the Recurrence Diagnosis and Prognosis of Bladder Cancer

This study aimed to investigate the role of infiltrating immune cell types in diagnosing and predicting bladder cancer recurrence. This study mainly applied some algorithms, including Estimate the Proportion of Immune and Cancer Cells (EPIC), support vector machine-recursive feature elimination (SVM-RFE), random forest out-of-bag (RF-OOB) and least absolute shrinkage and selection operator (LASSO)-Cox regression analysis. We found six immune infiltrating cell types significantly associated with recurrence prognosis and two independent clinical prognostic factors. Infiltrating immune cell types (IICTs) based on the prognostic immune risk score (pIRS) models may provide significant biomarkers for the diagnosis and prognostic prediction of bladder cancer recurrence.

Novel molecular insights into pyroptosis in triple-negative breast cancer prognosis and immunotherapy

BACKGROUND

Patients with triple-negative breast cancer (TNBC) often have a poor prognostic outcome. Current treatment strategies cannot benefit all TNBC patients. Previous findings suggested pyroptosis as a novel target for suppressing cancer development, although the relationship between TNBC and pyroptosis-related genes (PRGs) was still unclear.

METHODS

Gene expression data and clinical follow-up of TNBC patients were collected from the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) and Gene Expression Omnibus (GEO). PRGs were screened using weighted gene co-expression network analysis. Cox regression analysis and the least absolute shrinkage and selection operator (i.e. LASSO) technique were applied to construct a pyroptosis-related prognostic risk score (PPRS) model, which was further combined with the clinicopathological characteristics of TNBC patients to develop a survival decision tree and a nomogram. The model was used to calculate the PPRS, and then the overall survival, immune infiltration, immunotherapy response and drug sensitivity of TNBC patients were analyzed based on the PPRS.

RESULTS

The PPRS model was closely related to clinicopathological features and can independently and accurately predict the prognosis of TNBC. According to normalized PPRS, patients in different cohorts were divided into two groups. Compared with the high-PPRS group, the low-PPRS group had significantly higher ESTIMATE (i.e. Estimation of STromal and Immune cells in MAlignantTumours using Expression data) score, immune score and stromal score, and it also had overexpressed immune checkpoints and significantly reduced Tumor Immune Dysfunction and Exclusion (TIDE) score, as well as higher sensitivity to paclitaxel, veliparib, olaparib and talazoparib. A decision tree and nomogram based on PPRS and clinical characteristics can improve the prognosis stratification and survival prediction for TNBC patients.

CONCLUSIONS

A PPRS model was developed to predict TNBC patients' immune characteristics and response to immunotherapy, chemotherapy and targeted therapy, as well as their survival outcomes.

Fatty acid metabolism-related genes as a novel module biomarker for kidney renal clear cell carcinoma: Bioinformatics modeling with experimental verification

BACKGROUNDS

Lipid metabolism reprogramming is a hallmark of cancer, however, the associations between fatty acid metabolism (FAM) and kidney renal clear cell carcinoma (KIRC) prognosis are still less investigated.

METHODS

The gene expression and clinical data of KIRC were obtained from TCGA. Using Cox regression and LASSO regression, a novel prognostic risk score model based on FAM-related genes was constructed, and a nomogram for prediction of overall survival rate of patients with KIRC was proposed. The correlation between risk score and the immune cell infiltration, immune-related function and tumor mutation burden (TMB) were explored. Finally, a hub gene was extracted from the model, and RT-qPCR, Western blot, Immunohistochemical, EdU, Scratch assay and Transwell experiments were conducted to validate and decipher the biomarker role of the hub gene in KIRC theranostics.

RESULTS

In this study, a novel risk score model and a nomogram were constructed based on 20 FAM-related genes to predict the prognosis of KIRC patients with AUC>0.7 at 1-, 3-, and 5-years. Patients in different subgroups showed different phenotypes in immune cell infiltration, immune-related function, TMB, and sensitivity to immunotherapy. In particular, the hub gene in the model, i.e., ACADM, was significantly down-expressed in human KIRC samples, and the knockdown of OCLN promoted proliferation, migration and invasion of KIRC cells in vitro.

CONCLUSIONS

In this study, a novel risk score model and a module biomarker based on FAM-related genes were screened for KIRC prognosis. More clinical carcinogenic validations will be performed for future translational applications of the findings.

Statistical considerations in long-term efficacy evaluation of anti-cancer therapies

Anti-cancer therapy has been a significant focus of research. Developing and marketing various types and mechanisms of anti-cancer therapies benefit a variety of patients significantly. The long-term benefit to patients in evaluating the risk-benefit ratio of anti-cancer therapy has become a significant concern. This paper discusses the evaluation of long-term efficacy within the estimand framework and summarizes the various strategies for addressing potential intercurrent events. Non-proportional hazards of survival data may arise with novel anti-cancer therapies, leading to potential bias in conventional evaluation methods. This paper reviews statistical methods for addressing this issue, including novel endpoints, hypothesis testing, and efficacy estimation methods. We also discuss the influences of treatment switching. Although advanced methods have been developed to address the non-proportional hazard, they still have limitations that require continued collaborative efforts to resolve issues.

OCLN as a novel biomarker for prognosis and immune infiltrates in kidney renal clear cell carcinoma: an integrative computational and experimental characterization

Background

Occludin (OCLN) is an important tight junction protein and has been reported to be abnormally expressed in the development of malignant tumors. However, its biomarker and carcinogenic roles in kidney renal clear cell carcinoma (KIRC) are less investigated.

Methods

The Cancer Genome Atlas database and Human Protein Atlas database were used to analyze the expression of OCLN in KIRC. UALCAN database and methylation-specific PCR assay were used to evaluate the methylation level of OCLN in KIRC. Univariate and multivariate Cox regression analyses were performed to model the prognostic significance of OCLN in KIRC patient cohorts. The correlation between OCLN expression and the immune cell infiltration, immune-related function and immune checkpoints were explored. Finally, EdU, scratch assay and transwell experiments were conducted to validate the role of OCLN in KIRC development.

Results

The expression of OCLN was significantly downregulated in KIRC, compared with normal renal tissues (p<0.001). Patients with low OCLN expression showed a worse prognosis and poorer clinicopathological characteristics. Functional enrichment analysis revealed that OCLN was mainly involved in biological processes such as immune response, immunoglobulin complex circulating and cytokine and chemokine receptor to mediate KIRC development. Immune-related analysis indicated that OCLN could potentially serve as a candidate target for KIRC immunotherapy. OCLN overexpression inhibited proliferation, migration and invasion of KIRC cells in vitro.

Conclusion

OCLN was validated as a candidate prognostic biomarker and therapeutic target of KIRC based both on computational and experimental approaches. More in vivo experiments will be conducted to decode its molecular mechanism in KIRC carcinogenesis in the future work.

Visualizing the (Causal) Effect of a Continuous Variable on a Time-To-Event Outcome

Visualization is a key aspect of communicating the results of any study aiming to estimate causal effects. In studies with time-to-event outcomes, the most popular visualization approach is depicting survival curves stratified by the variable of interest. This approach cannot be used when the variable of interest is continuous. Simple workarounds, such as categorizing the continuous covariate and plotting survival curves for each category, can result in misleading depictions of the main effects. Instead, we propose a new graphic, the survival area plot, to directly depict the survival probability over time and as a function of a continuous covariate simultaneously. This plot utilizes g-computation based on a suitable time-to-event model to obtain the relevant estimates. Through the use of g-computation, those estimates can be adjusted for confounding without additional effort, allowing a causal interpretation under the standard causal identifiability assumptions. If those assumptions are not met, the proposed plot may still be used to depict noncausal associations. We illustrate and compare the proposed graphics to simpler alternatives using data from a large German observational study investigating the effect of the Ankle-Brachial Index on survival. To facilitate the usage of these plots, we additionally developed the contsurvplot R-package, which includes all methods discussed in this paper.

Diagnostic value of preoperative examination for evaluating margin status in breast cancer

BACKGROUND

A positive resection margin is a major risk factor for local breast cancer recurrence after breast-conserving surgery (BCS). Preoperative imaging examinations are frequently employed to assess the surgical margin.

AIM

To investigate the role and value of preoperative imaging examinations [magnetic resonance imaging (MRI), molybdenum target, and ultrasound] in evaluating margins for BCS.

METHODS

A retrospective study was conducted on 323 breast cancer patients who met the criteria for BCS and consented to the procedure from January 2014 to July 2021. The study gathered preoperative imaging data (MRI, ultrasound, and molybdenum target examination) and intraoperative and postoperative pathological information. Based on their BCS outcomes, patients were categorized into positive and negative margin groups. Subsequently, the patients were randomly split into a training set (226 patients, approximately 70%) and a validation set (97 patients, approximately 30%). The imaging and pathological information was analyzed and summarized using R software. Non-conditional logistic regression and LASSO regression were conducted in the validation set to identify factors that might influence the failure of BCS. A column chart was generated and applied to the validation set to examine the relationship between pathological margin range and prognosis. This study aims to identify the risk factors associated with failure in BCS.

RESULTS

The multivariate non-conditional logistic regression analysis demonstrated that various factors raise the risk of positive margins following BCS. These factors comprise non-mass enhancement (NME) on dynamic contrast-enhanced MRI, multiple focal vascular signs around the lesion on MRI, tumor size exceeding 2 cm, type III time-signal intensity curve, indistinct margins on molybdenum target examination, unclear margins on ultrasound examination, and estrogen receptor (ER) positivity in immunohistochemistry. LASSO regression was additionally employed in this study to identify four predictive factors for the model: ER, molybdenum target tumor type (MT Xmd Shape), maximum intensity projection imaging feature, and lesion type on MRI. The model constructed with these predictive factors exhibited strong consistency with the real-world scenario in both the training set and validation set. Particularly, the outcomes of the column chart model accurately predicted the likelihood of positive margins in BCS.

CONCLUSION

The proposed column chart model effectively predicts the success of BCS for breast cancer. The model utilizes preoperative ultrasound, molybdenum target, MRI, and core needle biopsy pathology evaluation results, all of which align with the real-world scenario. Hence, our model can offer dependable guidance for clinical decision-making concerning BCS.

Diagnostic value of preoperative examination for evaluating margin status in breast cancer

BACKGROUND

A positive resection margin is a major risk factor for local breast cancer recurrence after breast-conserving surgery (BCS). Preoperative imaging examinations are frequently employed to assess the surgical margin.

AIM

To investigate the role and value of preoperative imaging examinations [magnetic resonance imaging (MRI), molybdenum target, and ultrasound] in evaluating margins for BCS.

METHODS

A retrospective study was conducted on 323 breast cancer patients who met the criteria for BCS and consented to the procedure from January 2014 to July 2021. The study gathered preoperative imaging data (MRI, ultrasound, and molybdenum target examination) and intraoperative and postoperative pathological information. Based on their BCS outcomes, patients were categorized into positive and negative margin groups. Subsequently, the patients were randomly split into a training set (226 patients, approximately 70%) and a validation set (97 patients, approximately 30%). The imaging and pathological information was analyzed and summarized using R software. Non-conditional logistic regression and LASSO regression were conducted in the validation set to identify factors that might influence the failure of BCS. A column chart was generated and applied to the validation set to examine the relationship between pathological margin range and prognosis. This study aims to identify the risk factors associated with failure in BCS.

RESULTS

The multivariate non-conditional logistic regression analysis demonstrated that various factors raise the risk of positive margins following BCS. These factors comprise non-mass enhancement (NME) on dynamic contrast-enhanced MRI, multiple focal vascular signs around the lesion on MRI, tumor size exceeding 2 cm, type III time-signal intensity curve, indistinct margins on molybdenum target examination, unclear margins on ultrasound examination, and estrogen receptor (ER) positivity in immunohistochemistry. LASSO regression was additionally employed in this study to identify four predictive factors for the model: ER, molybdenum target tumor type (MT Xmd Shape), maximum intensity projection imaging feature, and lesion type on MRI. The model constructed with these predictive factors exhibited strong consistency with the real-world scenario in both the training set and validation set. Particularly, the outcomes of the column chart model accurately predicted the likelihood of positive margins in BCS.

CONCLUSION

The proposed column chart model effectively predicts the success of BCS for breast cancer. The model utilizes preoperative ultrasound, molybdenum target, MRI, and core needle biopsy pathology evaluation results, all of which align with the real-world scenario. Hence, our model can offer dependable guidance for clinical decision-making concerning BCS.

Identification and validation of a DNA methylation-driven gene-based prognostic model for clear cell renal cell carcinoma

BACKGROUND

Clear cell renal cell carcinoma (ccRCC) is a malignant tumor with heterogeneous morphology and poor prognosis. This study aimed to establish a DNA methylation (DNAm)-driven gene-based prognostic model for ccRCC.

METHODS

Reduced representation bisulfite sequencing (RRBS) was performed on the DNA extracts from ccRCC patients. We analyzed the RRBS data from 10 pairs of patient samples to screen the candidate CpG sites, then trained and validated an 18-CpG site model, and integrated the clinical characters to establish a Nomogram model for the prognosis or risk evaluation of ccRCC.

RESULTS

We identified 2261 DMRs in the promoter region. After DMR selection, 578 candidates were screened, and was correspondence with 408 CpG dinucleotides in the 450 K array. We collected the DNAm profiles of 478 ccRCC samples from TCGA dataset. Using the training set with 319 samples, a prognostic panel of 18 CpGs was determined by univariate Cox regression, LASSO regression, and multivariate Cox proportional hazards regression analyses. We constructed a prognostic model by combining the clinical signatures. In the test set (159 samples) and whole set (478 samples), the Kaplan-Meier plot showed significant differences; and the ROC curve and survival analyses showed AUC greater than 0.7. The Nomogram integrated with clinicopathological characters and methylation risk score had better performance, and the decision curve analyses also showed a beneficial effect.

CONCLUSIONS

This work provides insight into the role of hypermethylation in ccRCC. The targets identified might serve as biomarkers for early ccRCC diagnosis and prognosis biomarkers for ccRCC. We believe our findings have implications for better risk stratification and personalized management of this disease.

Identification of fatty acid-related subtypes, the establishment of a prognostic signature, and immune infiltration characteristics in lung adenocarcinoma

Abnormal fatty acid (FA) metabolism can change the inflammatory microenvironment and promote tumor progression and metastasis, however, the potential association between FA-related genes (FARGs) and lung adenocarcinoma (LUAD) is still unclear. In this study, we described the genetic and transcriptomic changes of FARGs in LUAD patients and identified two different FA subtypes, which were significantly correlated with overall survival and tumor microenvironment infiltrating cells in LUAD patients. In addition, the FA score was also constructed through the LASSO Cox to evaluate the FA dysfunction of each patient. Multivariate Cox analysis proved that the FA score was an independent predictor and created the FA score integrated nomogram, which offered a quantitative tool for clinical practice. The performance of the FA score has been substantiated in numerous datasets for its commendable accuracy in estimating overall survival in LUAD patients. The groups with high and low FA scores exhibited different mutation spectrums, copy number variations, enrichment pathways, and immune status. Noteworthy differences between the two groups in terms of immunophenoscore and Tumor Immune Dysfunction and Exclusion were observed, suggesting that the group with a low FA score was more responsive to immunotherapy, and similar results were also confirmed in the immunotherapy cohort. In addition, seven potential chemotherapeutic drugs related to FA score targeting were predicted. Ultimately, we ascertained that the attenuation of KRT6A expression impeded the proliferation, migration, and invasion of LUAD cell lines. In summary, this research offers novel biomarkers to facilitate prognostic forecasting and clinical supervision for individuals afflicted with LUAD.

LOCC: a novel visualization and scoring of cutoffs for continuous variables

Objective

There is a need for new methods to select and analyze cutoffs employed to define genes that are most prognostic significant and impactful. We designed LOCC (Luo's Optimization Categorization Curve), a novel tool to visualize and score continuous variables for a dichotomous outcome.

Methods

To demonstrate LOCC with real world data, we analyzed TCGA hepatocellular carcinoma gene expression and patient data using LOCC. We compared LOCC visualization to receiver operating characteristic (ROC) curve for prognostic modeling to showcase its utility in understanding predictors in various TCGA datasets.

Results

Analysis of E2F1 expression in hepatocellular carcinoma using LOCC demonstrated appropriate cutoff selection and validation. In addition, we compared LOCC visualization and scoring to ROC curves and c-statistics, demonstrating that LOCC better described predictors. Analysis of a previously published gene signature showed large differences in LOCC scoring, and removing the lowest scoring genes did not affect prognostic modeling of the gene signature demonstrating LOCC scoring could distinguish which predictors were most critical.

Conclusion

Overall, LOCC is a novel visualization tool for understanding and selecting cutoffs, particularly for gene expression analysis in cancer. The LOCC score can be used to rank genes for prognostic potential and is more suitable than ROC curves for prognostic modeling.

Investigation of Hippo pathway-related prognostic lncRNAs and molecular subtypes in liver hepatocellular carcinoma

This study aimed to investigate Hippo pathway-related prognostic long noncoding RNAs (lncRNAs) and their prognostic value in liver hepatocellular carcinoma (LIHC). Expression and clinical data regarding LIHC were acquired from The Cancer Genome Atlas and European Bioinformatics Institute array databases. Hippo pathway-related lncRNAs and their prognostic value were revealed, followed by molecular subtype investigations. Differences in survival, clinical characteristics, immune cell infiltration, and checkpoint expression between the subtypes were explored. LASSO regression was used to determine the most valuable prognostic lncRNAs, followed by the establishment of a prognostic model. Survival and differential expression analyses were conducted between two groups (high- and low-risk). A total of 313 Hippo pathway-related lncRNAs were identified from LIHC, of which 88 were associated with prognosis, and two molecular subtypes were identified based on their expression patterns. These two subtypes showed significant differences in overall survival, pathological stage and grade, vascular invasion, infiltration abundance of seven immune cells, and expression of several checkpoints, such as CTLA-4 and PD-1/L1 (P < 0.05). LASSO regression identified the six most valuable independent prognostic lncRNAs for establishing a prognosis risk model. Risk scores calculated by the risk model assigned patients into two risk groups with an AUC of 0.913 and 0.731, respectively, indicating that the high-risk group had poor survival. The risk score had an independent prognostic value with an HR of 2.198. In total, 3007 genes were dysregulated between the two risk groups, and the expression of most genes was elevated in the high-risk group, involving the cell cycle and pathways in cancers. Hippo pathway-related lncRNAs could stratify patients for personalized treatment and predict the prognosis of patients with LIHC.

Landscape of the intratumroal microenvironment in bladder cancer: Implications for prognosis and immunotherapy

Introduction

This study aims to present the landscape of the intratumoral microenvironment and by which establish a classification system that can be used to predict the prognosis of bladder cancer patients and their response to anti-PD-L1 immunotherapy.

Methods

The expression profiles of 1554 bladder cancer cases were downloaded from seven public datasets. Single-sample gene set enrichment analysis (ssGSEA), univariate Cox regression analysis, and meta-analysis were employed to establish the bladder cancer immune prognostic index (BCIPI). Extensive analyses were executed to investigate the association between BCIPI and overall survival, tumor-infiltrated immunocytes, immunotherapeutic response, mutation load, etc.

Results

The results obtained from seven independent cohorts and meta-analyses suggested that the BCIPI is an effective classification system for estimating bladder cancer patients' overall survival. Patients in the BCIPI-High subgroup revealed different immunophenotypic outcomes from those in the BCIPI-Low subgroup regarding tumor-infiltrated immunocytes and mutated genes. Subsequent analysis suggested that patients in the BCIPI-High subgroup were more sensitive to anti-PD-L1 immunotherapy than those in the BCIPI-Low subgroup.

Conclusions

The newly established BCIPI is a valuable tool for predicting overall survival outcomes and immunotherapeutic responses in patients with bladder cancer.

Construction of a 10-gene prognostic score model of predicting recurrence for laryngeal cancer

We constructed a prognostic score (PS) model to predict the recurrence risk in patients previously diagnosed with laryngeal cancer (LC). Here the training dataset, consisting of 82 LC samples, was downloaded from The Cancer Genome Atlas (TCGA). The PS model then divided the LC samples into high- and low-risk groups, which predicted well the survival time of LC in three datasets (TCGA dataset: AUC = 0.899; GSE27020: AUC = 0.719; and GSE25727: AUC = 0.662). Therefore, the PS model based on the 10 genes and its nomogram is proposed to help predict the recurrence risk in patients with LC.

The types of tumor infiltrating lymphocytes are valuable for the diagnosis and prognosis of breast cancer

This study aimed at constructing a diagnostic immune risk score (dIRS) system and a prognostic immune risk score (pIRS) system for diagnose and prognosis of breast cancer (BC). The gene expression data of BC were downloaded from TCGA dataset (training set), and from GSE65194, GSE29044, GSE42568, and GSE20685 (validation sets). Then, the immune cell type proportions in each dataset were assessed using EPIC tool, and the dIRS system was built based on the SVM-RFE and RF-VIMP algorithms. Subsequently, the pIRS system and the nomogram survival model were established separately using penalized and rms packages. Finally, the differential expressed genes (DEGs) between low and high pIRS groups were screened, and submitted for functional analysis. The dIRS system consisted of B cells, CD8 + T cells, endothelial cells, NK cells, and other cells had high accuracy in distinguishing BC patients from the healthy controls (AUROC &gt;0.7). Subsequently, the pIRS system with the five prognosis-associated immune-infiltrating cell was constructed, and Kaplan-Meier analysis demonstrated that the survival rate of low pIRS group was significantly higher than that of high pIRS group (p &lt; 0.05). Based on age, pathologic stage and the pIRS values, the nomogram survival model was built. The AUROC value, Specificity value, Sensitivity value and C-index of the nomogram survival model were higher than 0.7000, and had a good predictive ability for BC. Finally, a total of 539 DEGs were identified, and significantly enriched in six pathways. The dIRS system and the pIRS system composed of immune cells might be critical for the diagnosis and prognosis of BC patients.

Increased expression of SYCP2 predicts poor prognosis in patients suffering from breast carcinoma

Overexpression of synaptonemal complex protein-2 (SYCP2) has been identified in various human papillomavirus (HPV)–related carcinomas, whereas its significant role in breast carcinoma remains unclear. The aim of this study was to elucidate the prognostic value and potential function of SYCP2 in breast carcinoma. Herein, data for breast carcinoma patients from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas database (TCGA) were analyzed. The enrichment analysis of SYCP2 including Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Friends, and GSEA was performed. Kaplan–Meier analysis, Cox regression, and receiver operating characteristic (ROC) curves were employed for determining the predictive value of SYCP2 on clinical outcomes in patients suffering from breast carcinoma. A nomogram was generated to predict the effect arising from SYCP2 on prognosis. The association analysis of SYCP2 gene expression and diverse immune infiltration levels was conducted through ssGSEA and ESTIMATE analysis, which consisted of dendritic cell (DC), neutrophil, eosinophil, macrophage, mast cell, NK cell, and other 18 cell subtypes. The results showed that SYCP2 expression was significantly elevated in breast carcinoma tissues as compared with that of normal tissues (p < 0.001). SYCP2 plays a certain role in pathways related to DNA methylation, keratinocyte differentiation, steroid hormone biosynthesis, and immune infiltration. The high expression of SYCP2 had a significant relationship to age, pathological type, ER expression, and PR expression (p < 0.001). Kaplan–Meier survival analysis showed that patients suffering from breast carcinoma characterized by high-SYCP2 expression had a poorer prognosis than patients with low-SYCP2 expression (p = 0.005). Univariate and multivariate Cox regression analyses revealed that SYCP2 had an independent relationship to overall survival (p = 0.049). Moreover, ROC curves suggested the significant diagnostic ability of SYCP2 for breast carcinoma, and as time went on, SYCP2 had more accurate prognostic efficacy. Furthermore, a high level of SYCP2 expression was found to have a relationship to poor prognosis of breast carcinoma in the subgroups of T3, N0, and M0, and infiltrating ductal carcinoma (HR > 1, p < 0.05). The calibration plot of the nomogram indicated that the SYCP2 model has an effective predictive performance for breast carcinoma patients. Conclusively, SYCP2 plays a vital role in the pathogenesis and progression of human breast carcinoma, so it may serve as a promising prognostic molecular marker of poor survival.

ATP6AP1 is a potential prognostic biomarker and is associated with iron metabolism in breast cancer

Cancer occurrence and progression may be facilitated by aberrant expression of ATPase H+ transporting accessory protein 1 (ATP6AP1). However, the clinical relevance of ATP6AP1 in breast cancer remains unclear. In this study, we investigated the association between ATP6AP1 and breast cancer. Data collected from patients with breast cancer from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) were used in this study. To determine the relationship between ATP6AP1 and breast cancer survival rates, Kaplan-Meier analysis was used. To determine the prognostic value of ATP6AP1, a receiver operating characteristic (ROC) curve was constructed. To identify the major pathways involving ATP6AP1, we performed functional enrichment analysis using gene set enrichment analysis (GSEA). We analyzed the association between ATP6AP1 expression and tumor immunity using the ESTIMATE algorithm and single-sample GSEA (ssGSEA). A nomogram based on a Cox regression analysis was constructed to predict the impact of ATP6AP1 on prognosis. ATP6AP1 expression was significantly upregulated in breast cancer tissues. Moreover, patients with elevated ATP6AP1 expression had shorter total survival rates than those with lower expression levels (p = 0.032). The area under the receiver operating characteristic curve for ATP6AP1 was 0.939. Gene set enrichment analysis revealed that reaction iron uptake and transport, proteasome degradation, glutathione metabolism, and pyruvate metabolism were enriched in the ATP6AP1 high expression phenotype. The relationship between immune infiltration cells and ATP6AP1 expression, including macrophages, B cells, dendritic cells, cytotoxic cells, NK cells, and T cells, was found to be negative, suggesting that ATP6AP1 overexpression results in immunosuppression. Based on the Cox regression analyses, the calibration plot of the nomogram demonstrated effective performance in predicting breast cancer patients. ATP6AP1 may facilitate breast cancer progression by inhibiting antitumor immunity and promoting iron metabolism and may be a biomarker for breast cancer prognosis.

Integrative Proteo-Genomic Analysis for Recurrent Survival Prognosis in Colon Adenocarcinoma

Background

The survival prognosis is the hallmark of cancer progression. Here, we aimed to develop a recurrence-related gene signature to predict the prognosis of colon adenocarcinoma (COAD).

Methods

The proteomic data from the Clinical Proteomic Tumor Analysis Consortium (CPTAC) and genomic data from the cancer genomic maps [The Cancer Genome Atlas (TCGA)] dataset were analyzed to identify co-differentially expressed genes (cDEGs) between recurrence samples and non-recurrence samples in COAD using limma package. Functional enrichment analysis, including Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway was conducted. Univariate and multivariate Cox regressions were applied to identify the independent prognostic feature cDEGs and establish the signature whose performance was evaluated by Kaplan–Meier curve, receiver operating characteristic (ROC), Harrell’s concordance index (C-index), and calibration curve. The area under the receiver operating characteristic (ROC) curve (AUROC) and a nomogram were calculated to assess the predictive accuracy. GSE17538 and GSE39582 were used for external validation. Quantitative real-time PCR and Western blot analysis were carried out to validate our findings.

Results

We identified 86 cDEGs in recurrence samples compared with non-recurrence samples. These genes were primarily enriched in the regulation of carbon metabolic process, fructose and mannose metabolism, and extracellular exosome. Then, an eight-gene-based signature (CA12, HBB, NCF1, KBTBD11, MMAA, DMBT1, AHNAK2, and FBLN2) was developed to separate patients into high- and low-risk groups. Patients in the low-risk group had significantly better prognosis than those in the high-risk group. Four prognostic clinical features, including pathological M, N, T, and RS model status, were screened for building the nomogram survival model. The PCR and Western blot analysis results suggested that CA12 and AHNAK2 were significantly upregulated, while MMAA and DMBT1 were downregulated in the tumor sample compared with adjacent tissues, and in non-recurrent samples compared with non-recurrent samples in COAD.

Conclusion

These identified recurrence-related gene signatures might provide an effective prognostic predictor and promising therapeutic targets for COAD patients.

Identification of Novel Characteristics in TP53-Mutant Hepatocellular Carcinoma Using Bioinformatics

Background: TP53 mutations are the most frequent mutations in hepatocellular carcinoma (HCC) and affect the occurrence and development of this cancer type. Therefore, it is essential to clarify the function and mechanism of TP53 mutations in HCC.

Methods: We performed a sequence of bioinformatic analyses to elucidate the characteristics of TP53 mutations in HCC. We downloaded the data of hepatocellular carcinoma from The Cancer Genome Atlas database and used different R packages for serial analyses, including gene mutation analysis, copy number variation analysis, analysis of the tumor mutational burden and microsatellite instability, differential gene expression analysis, and functional enrichment analysis of TP53 mutations, and performed gene set enrichment analysis. We established a protein-protein interaction network using the STRING online database and used the Cytoscape software for network visualization, and hub gene screening. In addition, we performed anticancer drug sensitivity analysis using data from the Genomics of Drug Sensitivity in Cancer. Immune infiltration and prognosis analyses were also performed.

Results: Missense mutations accounted for a great proportion of HCC mutations, the frequency of single nucleotide polymorphisms was high, and C > T was the most common form of single nucleotide variations. TP53 had a mutation rate of 30% and was the most commonly mutated gene in HCC. In the TP53 mutant group, the tumor mutational burden (p < 0.001), drug sensitivity (p < 0.05), ESTIMATE score (p = 0.038), and stromal score (p < 0.001) dramatically decreased. The Cytoscape software screened ten hub genes, including CT45A1, XAGE1B, CT55, GAGE2A, PASD1, MAGEA4, CTAG2, MAGEA10, MAGEC1, and SAGE1. The prognostic model showed a poor prognosis in the TP53 mutation group compared with that in the wild-type group (overall survival, p = 0.023). Univariate and multivariate cox regression analyses revealed that TP53 mutation was an independent risk factor for the prognosis of HCC patients (p <0.05). The constructed prognostic model had a favorable forecast value for the prognosis of HCC patients at 1 and 3 years (1-year AUC = 0.752, 3-years AUC = 0.702).

Conclusion: This study further deepened our understanding of TP53-mutated HCC, provided new insights into a precise individualized therapy for HCC, and has particular significance for prognosis prediction.

Identification and validation of an eight-lncRNA signature that predicts prognosis in patients with esophageal squamous cell carcinoma

Background

Esophageal squamous cell carcinoma (ESCC) is correlated with worse clinical prognosis and lacks available targeted therapy. Thus, identification of reliable biomarkers is required for the diagnosis and treatment of ESCC.

Methods

We downloaded the GSE53625 dataset as a training dataset to screen differentially expressed RNAs (DERs) with the criterion of false discovery rate (FDR) < 0.05 and |log₂fold change (FC)| > 1. A support vector machine classifier was used to find the optimal feature gene set that could conclusively distinguish different samples. An eight-lncRNA signature was identified by random survival forest algorithm and multivariate Cox regression analysis. The RNA sequencing data from The Cancer Genome Atlas (TCGA) database were used for external validation. The predictive value of the signature was assessed using Kaplan–Meier test, time-dependent receiver operating characteristic (ROC) curves, and dynamic area under the curve (AUC). Furthermore, a nomogram to predict patients’ 3-year and 5-year prognosis was constructed. CCK-8 assay, flow cytometry, and transwell assay were conducted in ESCC cells.

Results

A total of 1136 DERs, including 689 downregulated mRNAs, 318 upregulated mRNAs, 74 downregulated lncRNAs and 55 upregulated lncRNAs, were obtained in the GES53625 dataset. From the training dataset, we identified an eight-lncRNA signature, (ADAMTS9-AS1, DLX6-AS1, LINC00470, LINC00520, LINC01497, LINC01749, MAMDC2-AS1, and SSTR5-AS1). A nomogram based on the eight-lncRNA signature, age, and pathologic stage was developed and showed good accuracy for predicting 3-year and 5-year survival probability of patients with ESCC. Functionally, knockdown of LINC00470 significantly suppressed cell proliferation, G1/S transition, and migration in two ESCC cell lines (EC9706 and TE-9). Moreover, knockdown of LINC00470 downregulated the protein levels of PCNA, CDK4, and N-cadherin, while upregulating E-cadherin protein level in EC9706 and TE-9 cells.

Conclusion

Our eight-lncRNA signature and nomogram can provide theoretical guidance for further research on the molecular mechanism of ESCC and the screening of molecular markers.

Supplementary Information

The online version contains supplementary material available at 10.1186/s11658-022-00331-x.

Identification of Novel Molecular Therapeutic Targets and Their Potential Prognostic Biomarkers Based on Cytolytic Activity in Skin Cutaneous Melanoma

Skin cutaneous melanoma (SKCM) attracts attention worldwide for its extremely high malignancy. A novel term cytolytic activity (CYT) has been introduced as a potential immunotherapy biomarker associated with counter-regulatory immune responses and enhanced prognosis in tumors. In this study, we extracted all datasets of SKCM patients, namely, RNA sequencing data and clinical information from The Cancer Genome Atlas (TCGA) database and the Gene Expression Omnibus (GEO) database, conducted differential expression analysis to yield 864 differentially expressed genes (DEGs) characteristic of CYT and used non-negative matrix factorization (NMF) method to classify molecular subtypes of SKCM patients. Among all genes, 14 hub genes closely related to prognosis for SKCM were finally screen out. Based on these genes, we constructed a 14-gene prognostic risk model and its robustness and strong predictive performance were further validated. Subsequently, the underlying mechanisms in tumor pathogenesis and prognosis have been defined from a number of perspectives, namely, tumor mutation burden (TMB), copy number variation (CNV), tumor microenvironment (TME), infiltrating immune cells, gene set enrichment analysis (GSEA) and immune checkpoint inhibitors (ICIs). Furthermore, combined with GTEx database and HPA database, the expression of genes in the model was verified at the transcriptional level and protein level, and the relative importance of genes in the model was described by random forest algorithm. In addition, the model was used to predict the difference in sensitivity of SKCM patients to chemotherapy and immunotherapy. Finally, a nomogram was constructed to better aid clinical diagnosis.

A Ferroptosis-Related Genes Model Allows for Prognosis and Treatment Stratification of Clear Cell Renal Cell Carcinoma: A Bioinformatics Analysis and Experimental Verification

Introduction

Clear cell renal cell carcinoma (ccRCC) is a malignant tumor characterized by poor prognosis and difficult treatment. Ferroptosis is a relatively new form of programmed cell death that involved in cancer development and therapy resistance. Studies have shown that targeted ferroptosis may be a novel option for the treatment of ccRCC, but key genes and their roles between ferroptosis and ccRCC are limited so far. This study aims to develop a ccRCC stratified model based on ferroptosis-related genes to provide a reference for the prognosis prediction and the individualized treatment of ccRCC.

Materials and Methods

The mRNAs expression data of ccRCC and FRGs were obtained from TCGA and FerrDb database, respectively. Through multiple analysis, a 4-FRG based prognostic stratified model was constructed and its predictive performance was validated through various methods. Then, a nomogram based on the model was constructed and ccRCC patients stratified by the model were analyzed for tumor microenvironment, immune infiltration, sensitivity for immune checkpoint inhibitors (ICIs)/traditional anti-tumor therapy and tumor mutation burden (TMB). Functional enrichment analysis was performed to explore potential biological pathways. Finally, we verified our model by RT-qPCR, siRNA transfection, scratch assay and CCK-8 assay.

Results

In this study, the stratified model and a model-based nomogram can accurately predict the prognosis of ccRCC patients in TCGA database. The patients stratified by the model showed different tumor microenvironments, immune infiltration, TMB, resistance to traditional and ICIs therapy, and sensitivity to ferroptosis. Functional enrichment analysis suggested several biological pathways related to the process and prognosis of ccRCC. RT-qPCR confirmed the differential expression of ferroptosis-related genes. Scratch assay and CCK-8 assay indicated the promotion effects of CD44 on the proliferation and migration of ccRCC.

Conclusion

In this study, we established a novel ccRCC stratified model based on FRGs, which can accurately predict the prognosis of ccRCC patients and provide a reference for clinical individualized treatment.

Combined mRNAs and clinical factors model on predicting prognosis in patients with triple-negative breast cancer

OBJECTIVE

Triple-negative breast cancer (TNBC) is aggressive cancer usually diagnosed in young women with no effective prognosis prediction model to use. The present study was performed to develop a useful prognostic model for predicting overall survival (OS) for TNBC patients.

METHODS

The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) databases were used as training and validation data sets, respectively, in which the gene expression levels and clinical prognostic information of TNBC were collected. Differentially expressed genes (DEGs) between TNBC and non-TNBC (NTNBC) were identified with the thresholds of false discovery rate < 0.05 and |log2 Fold Change| > 1. DEGs in AmiGO2 and the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases were retained for further study. Univariate, multivariate Cox, and logistic regression analysis were conducted for detecting DEG signature with the threshold of log-rank P < 0.05. The prognosis models of mRNA signature, clinical factors were constructed and compared.

RESULTS

One five-DEG signature, including CHST4, COCH, CST9, SOX11, and TDGF1 was identified in DEG prognosis model. Stratified analysis showed that the patients aged over 60, with higher pathologic stage (III-IV) and recurrence induced a significantly lower survival rate than those aged below 60, lower pathologic stage and without recurrence. Compared with patients with low-risk scores, those presented high-risk scores demonstrated significantly lower survival rate in the subgroup aged over 60 [HR = 3.780 (1.801-7.933), P < 0.0001]. For patients who obtained a higher pathologic stage and recurrence, high-risk scores were correlated with a significantly lower survival rate than patients with low-risk scores. The five-mRNA signature combined with clinical model (AUC = 0.950) predicted better than single clinical model (AUC = 0.795) or five-mRNA signature model (AUC = 0.823).

CONCLUSION

Our present study identified a prognostic prediction model (combined with five-mRNA signature and clinical factors) for TNBC patients receiving immunotherapy, which will benefit future research and clinical therapies.

Kernel hazard estimation for visualisation of the effect of a continuous covariates on time-to-event endpoints

The problem of associating a continuous covariate, or biomarker, against a time-to-event outcome, is that it often requires categorisation of the covariate. This can lead to bias, loss of information and a poor representation of any underlying relationship. Here, two methods are proposed for estimating the effects of a continuous covariate on a time-to-event endpoint using weighted kernel estimators. The first method aims to estimate a density function for a time-to-event endpoint conditional on some covariate value whilst the second uses a joint density estimator. The results are visualisations in the form of surface plots that show the effects of a covariate without any need for categorisation. Both methods can aid interpretation and analysis of covariates against a time-to-event endpoint.

Development and validation of a four-lipid metabolism gene signature for diagnosis of pancreatic cancer

Abnormal lipid metabolism is closely related to the malignant biological behavior of tumor cells. Such abnormal lipid metabolism provides energy for rapid proliferation, and certain genes related to lipid metabolism encode important components of tumor signaling pathways. In this study, we analyzed pancreatic cancer datasets from The Cancer Genome Atlas and searched for prognostic genes related to lipid metabolism in the Molecular Signature Database. A risk score model was built and verified using the GSE57495 dataset and International Cancer Genome Consortium dataset. Four molecular subtypes and 4249 differentially expressed genes (DEGs) were identified. The DEGs obtained by Weighted Gene Coexpression Network Construction analysis were intersected with 4249 DEGs to obtain a total of 1340 DEGs. The final prognosis model included CA8, CEP55, GNB3 and SGSM2, and these had a significant effect on overall survival. The area under the curve at 1, 3 and 5 years was 0.72, 0.79 and 0.87, respectively. These same results were obtained using the validation cohort. Survival analysis data showed that the model could stratify the prognosis of patients with different clinical characteristics, and the model has clinical independence. Functional analysis indicated that the model is associated with multiple cancer-related pathways. Compared with published models, our model has a higher C-index and greater risk value. In summary, this four-gene signature is an independent risk factor for pancreatic cancer survival and may be an effective prognostic indicator.

Screening of Significant Biomarkers Related to Prognosis of Cervical Cancer and Functional Study Based on lncRNA-associated ceRNA Regulatory Network

BACKGROUND

Cervical cancer (CESC), which threatens the health of women, has a very high recurrence rate.

PURPOSES

This study aimed to identify the signature long non-coding RNAs (lncRNAs) associated with the prognosis of CESC and predict the prognostic survival rate with the clinical risk factors.

METHODS

The CESC gene expression profiling data were downloaded from TCGA database and NCBI Gene Expression Omnibus. Afterwards, the differentially expressed RNAs (DERs) were screened using limma package of R software. R package "survival" was then used to screen the signature lncRNAs associated with independently recurrence prognosis, and a nomogram recurrence rate model based on these signature lncRNAs was constructed to predict the 3-year and 5-year survival probability of CESC. Finally, a competing endogenous RNAs (ceRNA) regulatory network was proposed to study the functions of these genes.

RESULTS

We obtained 305 DERs significantly associated with prognosis. Afterwards, a risk score (RS) prediction model was established using the screened 5 signature lncRNAs associated with independently recurrence prognosis (DLEU1, LINC01119, RBPMS-AS1, RAD21-AS1 and LINC00323). Subsequently, a nomogram recurrence rate model, proposed with Pathologic N and RS model status, was found to have a good prediction ability for CESC. In ceRNA regulatory network, LINC00323 and DLEU1 were hub nodes which targeted more miRNAs and mRNAs. After that, 15 GO terms and 3 KEGG pathways were associated with recurrence prognosis and showed that the targeted genes PTK2, NRP1, PRKAA1 and HMGCS1 might influence the prognosis of CESC.

CONCLUSION

The signature lncRNAs can help improve our understanding of the development and recurrence of CESC and the nomogram recurrence rate model can be applied to predict the survival rate of CESC patients in clinical practice.

Comprehensive Transcriptomic Analysis of Critical RNA Regulation Associated With Metabolism and Prognosis in Clear Cell Renal Carcinoma

This study focuses on investigating the metabolism-related gene profile and prognosis of clear cell renal cell carcinoma (ccRCC) patients. The research data from the Gene Expression Omnibus database, including GSE40435, GSE53757, and GSE53000, were used to analyze the consistently differentially expressed RNAs (cDERs) by the MetaDE limma package. Gene expression profiling associated with metabolism was downloaded from the GSEA database. The cancer genome atlas (TCGA) dataset of ccRCC (the training set) and RNA sequencing data of E-MTAB-3267 from EBI ArrayExpress database (the validation set) were obtained to construct a prognostic model. A series of bioinformatics analysis, including functional enrichment analysis, Cox regression analysis, and constructing a prognostic score (PS) model, was performed. Further in vitro experiments including cell proliferation assay and flow cytometry were performed to validate our results. We constructed a metabolism-related prognostic model based on 27 DElncRNAs and 126 DEGs. Gene Set Enrichment Analysis revealed that 19 GO terms and 9 KEGG signaling pathways were significantly associated with lipid metabolic pathways. Furthermore, we generated a nomogram illustrating the association between the identified DERs and the tumor recurrence risk in ccRCC. The results from experimental validation showed that lncRNA SNHG20 was significantly upregulated in tumor tissues compared with adjacent tissues. Knockdown of SNHG20 suppressed the proliferation and induced cell cycle G0/G1 arrest, and apoptosis in ccRCC cells. Our study might contribute to a better understanding of metabolic pathways and to the further development of novel therapeutic approaches for ccRCC.

Tumor Mutation Burden-Associated LINC00638/miR-4732-3p/ULBP1 Axis Promotes Immune Escape via PD-L1 in Hepatocellular Carcinoma

Tumor mutation burden (TMB) is associated with immune infiltration, while its underlying mechanism in hepatocellular carcinoma (HCC) remains unclear. A long noncoding RNA (lncRNA)-related competitive endogenous RNA (ceRNA) network can regulate various tumor behaviors, and research about its correlation with TMB and immune infiltration is warranted. Data were downloaded from TCGA and ArrayExpress databases. Cox analysis and machine learning algorithms were employed to establish a lncRNA-based prognostic model for HCC. We then developed a nomogram model to predict overall survival and odds of death for HCC patients. The association of this prognostic model with TMB and immune infiltration was also analyzed. In addition, a ceRNA network was constructed by using DIANA-LncBasev2 and the starBase database and verified by luciferase reporter and colocalization analysis. Multiplex immunofluorescence was applied to determine the correlation between ULBP1 and PD-L1. An eight-lncRNA (SLC25A30-AS1, HPN-AS1, LINC00607, USP2-AS1, HCG20, LINC00638, MKLN1-AS and LINC00652) prognostic score model was constructed for HCC, which was highly associated with TMB and immune infiltration. Next, we constructed a ceRNA network, LINC00638/miR-4732-3p/ULBP1, that may be responsible for NK cell infiltration in HCC with high TMB. However, patients with high ULBP1 possessed a poorer prognosis. Using multiplex immunofluorescence, we found a significant correlation between ULBP1 and PD-L1 in HCC, and patients with high ULBP1 and PD-L1 had the worst prognosis. In brief, the eight-lncRNA model is a reliable tool to predict the prognosis of HCC patients. The LINC00638/miR-4732-3p/ULBP1 axis may regulate immune escape via PD-L1 in HCC with high TMB.

Comprehensive Characterization of RNA Processing Factors in Gastric Cancer Identifies a Prognostic Signature for Predicting Clinical Outcomes and Therapeutic Responses

RNA processing converts primary transcript RNA into mature RNA. Altered RNA processing drives tumor initiation and maintenance, and may generate novel therapeutic opportunities. However, the role of RNA processing factors in gastric cancer (GC) has not been clearly elucidated. This study presents a comprehensive analysis exploring the clinical, molecular, immune, and drug response features underlying the RNA processing factors in GC. This study included 1079 GC cases from The Cancer Genome Atlas (TCGA, training set), our hospital cohort, and two other external validation sets (GSE15459, GSE62254). We developed an RNA processing-related prognostic signature using Cox regression with the least absolute shrinkage and selection operator (LASSO) penalty. The prognostic value of the signature was evaluated using a multiple-method approach. The genetic variants, pathway activation, immune heterogeneity, drug response, and splicing features associated with the risk signature were explored using bioinformatics methods. Among the tested 819 RNA processing genes, we identified five distinct RNA processing patterns with specific clinical outcomes and biological features. A 10-gene RNA processing-related prognostic signature, involving ZBTB7A, METTL2B, CACTIN, TRUB2, POLDIP3, TSEN54, SUGP1, RBMS1, TGFB1, and PWP2, was further identified. The signature was a powerful and robust prognosis factor in both the training and validation datasets. Notably, it could stratify the survival of patients with GC in specific tumor-node-metastasis (TNM) classification subgroups. We constructed a composite prognostic nomogram to facilitate clinical practice by integrating this signature with other clinical variables (TNM stage, age). Patients with low-risk scores were characterized with good clinical outcomes, proliferation, and metabolism hallmarks. Conversely, poor clinical outcome, invasion, and metastasis hallmarks were enriched in the high-risk group. The RNA processing signature was also involved in tumor microenvironment reprogramming and regulating alternative splicing, causing different drug response features between the two risk groups. The low-risk subgroup was characterized by high genomic instability, high alternative splicing and might benefit from the immunotherapy. Our findings highlight the prognostic value of RNA processing factors for patients with GC and provide insights into the specific clinical and molecular features underlying the RNA processing-related signature, which may be important for patient management and targeting treatment.

Identification of the Prognostic Signatures of Glioma With Different PTEN Status

The high-grade glioma is characterized by cell heterogeneity, gene mutations, and poor prognosis. The deletions and mutations of the tumor suppressor gene PTEN (5%-40%) in glioma patients are associated with worse survival and therapeutic resistance. Characterization of unique prognosis molecular signatures by PTEN status in glioma is still unclear. This study established a novel risk model, screened optimal prognostic signatures, and calculated the risk score for the individual glioma patients with different PTEN status. Screening results revealed fourteen independent prognostic gene signatures in PTEN-wt and three in the -50PTEN-mut subgroup. Moreover, we verified risk score as an independent prognostic factor significantly correlated with tumor malignancy. Due to the higher malignancy of the PTEN-mut gliomas, we explored the independent prognostic signatures (CLCF1, AEBP1, and OS9) for a potential therapeutic target in PTEN-mut glioma. We further separated IDH wild-type glioma patients into GBM and LGG to verify the therapeutic target along with PTEN status, notably, the above screened therapeutic targets are also significant prognostic genes in both IDH-wt/PTEN-mut GBM and LGG patients. We further identified the small molecule compound (+)-JQ1 binds to all three targets, indicating a potential therapy for PTEN-mut glioma. In sum, gene signatures and risk scores in the novel risk model facilitate glioma diagnosis, prognosis prediction, and treatment.

Adjusting for covariates in analysis based on restricted mean survival times

We summarize extensions to the analysis of restricted mean survival time (RMST) in the context of time-to-event outcomes. The RMST estimate and its inference are based on the classical Kaplan-Meier curves. When covariate effects are considered, an adjusted RMST (ARMST) estimate can be derived analogously based on adjusted Kaplan-Meier curves. The adjusted Kaplan-Meier Estimator (AKME) was developed to reduce confounding by the method of inverse probability of treatment weighting. We will show how the ARMST method combines the concepts of the RMST and AKME to make inferences. Two regression based methods to adjust for potential covariate effect on the RMST estimates will be compared with the ARMST approach. Simulation studies are performed to compare the different methods with and without covariate adjustments. In addition, we will summarize the extension of RMST and ARMST to the setting with competing risks. The restricted mean time lost (RMTL) and adjusted RMTL (ARMTL) are estimates of interest from cumulative incidence curves. A phase 3 oncology clinical trial example is provided to demonstrate the applications of these methods.

Integrative analysis of ceRNA network reveals functional lncRNAs associated with independent recurrent prognosis in colon adenocarcinoma

Background

Long non-coding RNAs (lncRNAs), acting as competing endogenous RNA (ceRNA) have been reported to regulate the expression of targeted genes by sponging miRNA in colon adenocarcinoma (COAD).

Methods

However, their potential implications for recurrence free survival prognosis and functional roles remains largely unclear in COAD. In this study, we downloaded the TCGA dataset (training dataset) and GSE39582 (validation dataset) of COAD patients with prognostic information.

Results

A total of 411 differentially expressed genes (DElncRNAs: 12 downregulated and 43 upregulated), 18 DE miRNAs (9 downregulated and 9 upregulated) and 338 DEmRNAs (113 downregulated and 225 upregulated) were identified in recurrence samples compared with non-recurrence samples with the thresholds of FDR < 0.05 and |log₂FC|> 0.263. Based on six signature lncRNAs (LINC00899, LINC01503, PRKAG2-AS1, RAD21-AS1, SRRM2-AS1 and USP30-AS1), the risk score (RS) system was constructed. Two prognostic clinical features, including pathologic stage and RS model status were screened for building the nomogram survival model. Moreover, a recurrent-specific ceRNA network was successfully constructed with 2 signature lncRNAs, 4 miRNAs and 113 mRNAs. Furthermore, we further manifested that SRRM2-AS1 predicted a poor prognosis in COAD patients. Furthermore, knockdown of SRRM2-AS1 significantly suppressed cell proliferation, migration, invasion and EMT markers in HT-29 and SW1116 cells.

Conclusion

These identified novel lncRNA signature and ceRNA network associated with recurrence prognosis might provide promising therapeutic targets for COAD patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02069-6.

A Novel miRNA-Based Model Can Predict the Prognosis of Clear Cell Renal Cell Carcinoma

Clear cell renal cell carcinoma (ccRCC) is one of the most prevalent renal malignant cancer, whose survival rate and quality of life of patients are still not satisfactory. Nevertheless, the TNM staging system currently used in clinical cannot make accurate survival predictions and precise treatment decisions for ccRCC patients. Therefore, there is an urgent need for more reliable biomarkers to identify high-risk subgroups of ccRCC patients to guide timely intervention and treatment. Recently, MiRNAs have been shown to be closely related to the procession of a variety of tumors, and they have high stability in various tissues, which makes them suggested to have the potential as a prognostic biomarker of ccRCC. In this study, by analyzing and processing the miRNAs expression profile of ccRCC patients from the TCGA database, we finally constructed an excellent miRNAs signature and verified it through a variety of methods. In order to build a more accurate and reliable clinical predictive model, we integrated the miRNAs signature with other prognostic-related clinical parameters to construct a nomogram. Functional enrichment analysis showed that miRNAs in the signature may regulate the genes involved in the Hippo signaling pathway, Tight junction, and Wnt signaling pathway to cause different prognoses of ccRCC patients, which may provide a reference for subsequent basic research and targeted therapy. To conclude, our study constructed a useful miRNAs signature, which allows the prognosis stratification for ccRCC patients and thereby guides the timely and effective interventions on high-risk patients. At the same time, this study also found the potential biological pathways involved in the procession of ccRCC.

Prognostic score model-based signature genes for predicting the prognosis of metastatic skin cutaneous melanoma

PURPOSE

Cutaneous melanoma (SKCM) is the most invasive malignancy of skin cancer. Metastasis to distant lymph nodes or other system is an indicator of poor prognosis in melanoma patients. The aim of this study was to identify reliable prognostic biomarkers for SKCMs.

METHODS

Four RNA-sequencing datasets associated with SKCMs were downloaded from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) database as well as corresponding clinical information. Differentially expressed genes (DEGs) were screened between primary and metastatic samples by using MetaDE tool. Weighted gene co-expression network analysis (WGCNA) was conducted to screen functional modules. A prognostic score (PS)-based predictive model and nomogram model were constructed to identify signature genes and independent clinicopathologic factors.

RESULTS

Based on MetaDE analysis and WGCNA, a total of 456 overlapped genes were identified as hub genes related to SKCMs progression. Functional enrichment analysis revealed these genes were mainly involved in the hippo signaling pathway, signaling pathways regulating pluripotency of stem cells, pathways in cancer. In addition, eight optimal DEGs (RFPL1S, CTSV, EGLN3, etc.) were identified as signature genes by using PS model. Cox regression analysis revealed that pathologic stage T, N and recurrence were independent prognostic factors. Three clinical factors and PS status were incorporated to construct a nomogram predictive model for estimating the three years and five-year survival probability of individuals.

CONCLUSIONS

The prognosis prediction model of this study may provide a promising method for decision making in clinic and prognosis predicting of SKCM patients.

A genomic-clinicopathologic Nomogram for the preoperative prediction of lymph node metastasis in gastric cancer

BACKGROUND

Preoperative evaluation of lymph node (LN) state is of pivotal significance for informing therapeutic decisions in gastric cancer (GC) patients. However, there are no non-invasive methods that can be used to preoperatively identify such status. We aimed at developing a genomic biosignature based model to predict the possibility of LN metastasis in GC patients.

METHODS

We used the RNA profile retrieving strategy and performed RNA expression profiling in a large GC cohort (GSE62254, n = 300) from Gene Expression Ominus (GEO). In the exploratory stage, 300 GC patients from GSE62254 were involved and the differentially expressed RNAs (DERs) for LN-status were determined using the R software. GC samples in GSE62254 were randomly allocated into a learning set (n = 210) and a verification set (n = 90). By using the Least absolute shrinkage and selection operator (LASSO) regression approach, a set of 23-RNA signatures were established and the signature based nomogram was subsequently built for distinguishing LN condition. The diagnostic efficiency, as well as the clinical performance of this model were assessed using the decision curve analysis (DCA). Metascape was used for bioinformatic analysis of the DERs.

RESULTS

Based on the genomic signature, we established a nomogram that robustly distinguished LN status in the learning (AUC = 0.916, 95% CI 0.833-0.999) and verification sets (AUC = 0.775, 95% CI 0.647-0.903). DCA demonstrated the clinical value of this nomogram. Functional enrichment analysis of the DERs was performed using bioinformatics methods which revealed that these DERs were involved in several lymphangiogenesis-correlated cascades.

CONCLUSIONS

In this study, we present a genomic signature based nomogram that integrates the 23-RNA biosignature based scores and Lauren classification. This model can be utilized to estimate the probability of LN metastasis with good performance in GC. The functional analysis of the DERs reveals the prospective biogenesis of LN metastasis in GC.

Development and verification of an immune-related gene pairs prognostic signature in ovarian cancer

Ovarian cancer (OV) is the most common gynaecological cancer worldwide. Immunotherapy has recently been proven to be an effective treatment strategy. The work here attempts to produce a prognostic immune-related gene pair (IRGP) signature to estimate OV patient survival. The Gene Expression Omnibus (GEO) and Cancer Genome Atlas (TCGA) databases provided the genetic expression profiles and clinical data of OV patients. Based on the InnateDB database and the least absolute shrinkage and selection operator (LASSO) regression model, we first identified a 17-IRGP signature associated with survival. The average area under the curve (AUC) values of the training, validation, and all TCGA sets were 0.869, 0.712, and 0.778, respectively. The 17-IRGP signature noticeably split patients into high- and low-risk groups with different prognostic outcomes. As suggested by a functional study, some biological pathways, including the Toll-like receptor and chemokine signalling pathways, were significantly negatively correlated with risk scores; however, pathways such as the p53 and apoptosis signalling pathways had a positive correlation. Moreover, tumour stage III, IV, grade G1/G2, and G3/G4 samples had significant differences in risk scores. In conclusion, an effective 17-IRGP signature was produced to predict prognostic outcomes in OV, providing new insights into immunological biomarkers.

A Risk Score System Based on the Methylation Levels of 15 RNAs in Breast Cancer

Background: Breast cancer (BC) occurs in the epithelial tissues of the breast gland, which is the most common cancer in women. This study is implemented to construct a risk score system for BC. Methods: The methylation data of BC from The Cancer Genome Atlas database (the training set) and GSE37754 from Gene Expression Omnibus database (the validation set) were downloaded. The differentially methylated RNAs (DMRs) between BC and normal samples were screened by limma package, and the correlations between the expression levels and methylation levels of the DMRs were analyzed to calculate their Pearson correlation coefficients (PCCs) using the cor.test function. To build the risk score system, the optimal RNAs were identified by penalized package. Subsequently, the nomogram survival model was established using the rms package. The lncRNA-mRNA comethylation network was constructed by Cytoscape software, and then enrichment analysis was performed using DAVID tool. Results: From the 1170 DMRs between BC and normal samples, 800 DMRs with significant negative PCCs were screened. For building the risk score system, the 15 optimal RNAs were selected. Afterward, the nomogram survival model based on four independent clinical prognostic factors (including age, radiation therapy, tumor recurrence, and RS model status) was constructed. In the comethylation network, the long noncoding RNA (lncRNA) PRNT was comethylated with FAM19A5 and RBM24. For the mRNAs in the comethylation network, angiogenesis and pathways in cancer were enriched. Conclusion: The risk score system and the nomogram survival model might be of great importance for the prognosis prediction of BC patients.

Identifying biomolecules and constructing a prognostic risk prediction model for recurrence in osteosarcoma

INTRODUCTION

Osteosarcoma is a high-morbidity bone cancer with an unsatisfactory prognosis. The aim of this study is to develop novel potential prognostic biomarkers and construct a prognostic risk prediction model for recurrence in osteosarcoma.

METHODS

By analyzing microarray data, univariate and multivariate Cox regression analyses were performed to screen prognostic RNA signatures and to build a prognostic model. The RNA signatures were validated using Kaplan-Meier curves. Then, we developed and validated a nomogram combining age, recurrence, metastatic, and Prognostic score (PS) models to predict the individual's overall survival at the 3- and 5-year points. Pathway enrichment of RNA was conducted based on the significant co-expressed RNAs.

RESULTS

A total of 319 mRNAs and 14 lncRNAs were identified in the microarray data. One lncRNA (LINC00957) and six mRNAs (METL1, CA9, B3GALT4, ALDH1A1, LAMB3, and ITGB4) were identified as RNA signatures and showed good performances in survival prediction for both the training and validation cohorts. Cox regression analysis showed that the seven RNA signatures could independently predict overall survival. Furthermore, age, recurrence, metastatic, and PS models were identified as independent prognostic factors via univariate and multivariate Cox analyses (P < 0.05) and included in the prognostic nomogram. The C-index values for the 3- and 5-year overall survival predictions of the nomogram were 0.809 and 0.740, respectively.

CONCLUSIONS

The current study provides the novel potential of seven RNA candidates as prognostic biomarkers. Nomograms were constructed to provide accurate and individualized survival prediction for recurrence in osteosarcoma patients.