Molecular subtyping of cancer: current status and moving toward clinical applications

Comprehensive analysis of senescence-related genes identifies prognostic clusters with distinct characteristics in glioma

Cellular senescence, defined as a state of permanent arrest in cell growth, is regarded as a crucial tumor suppression mechanism. However, accumulating scientific evidence suggests that senescent cells play a detrimental role in the progression of cancer. Unfortunately, the current lack of reliable markers that specifically reflect the level of senescence in cancer greatly hinders our in-depth understanding of this important biological foundation. Therefore, the search for more specific and reliable markers to reveal the specific role of senescent cells in cancer progression is particularly urgent and important. To uncover the role of senescence in gliomas, we collected senescence-related genes for integrated analysis. Consensus clustering was used to subtype gliomas based on the senescence gene set, and we identified two robust prognostic clusters of gliomas with distinct survival outcomes, multi-omics landscapes, immune characteristics, and differential drug responses. Multiple external datasets were used to validate the stability of our subtypes. Various computational and experimental methods, including WGCNA (Weighted Gene Co-expression Network Analysis), ssGSEA (single-sample Gene Set Enrichment Analysis), and machine learning algorithms (lasso regression, support vector machines, random forests), were employed for analysis. We found that CEBPB and LMNA are associated with poor prognosis in gliomas and may mediate immunosuppression and tumor proliferation. Drug prediction indicated that dasatinib is a potential therapeutic agent. Our findings provide insights into the role of the senescence gene set in patient stratification and precision medicine.

NetSDR: Drug repurposing for cancers based on subtype-specific network modularization and perturbation analysis

Cancer, a heterogeneous disease, presents significant challenges for drug development due to its complex etiology. Drug repurposing, particularly through network medicine approaches, offers a promising avenue for cancer treatment by analyzing how drugs influence cellular networks on a systemic scale. The advent of large-scale proteomics data provides new opportunities to elucidate regulatory mechanisms specific to cancer subtypes. Herein, we present NetSDR, a Network-based Subtype-specific Drug Repurposing framework for prioritizing repurposed drugs specific to certain cancer subtypes, guided by subtype-specific proteomic signatures and network perturbations. First, by integrating cancer subtype information into a network-based method, we developed a pipeline to recognize subtype-specific functional modules. Next, we conducted drug response analysis for each module to identify the "therapeutic module" and then used deep learning to construct weighted drug response network for the particular subtype. Finally, we employed a perturbation response scanning-based drug repurposing method, which incorporates dynamic information, to facilitate the prioritization of candidate drugs. Applying the framework to gastric cancer, we attested the significance of the extracellular matrix module in treatment strategies and discovered a promising potential drug target, LAMB2, as well as a series of possible repurposed drugs. This study demonstrates a systems biology framework for precise drug repurposing in cancer and other complex diseases.

Integrative analysis of senescence-related genes identifies robust prognostic clusters with distinct features in hepatocellular carcinoma

INTRODUCTION

Senescence refers to a state of permanent cell growth arrest and is regarded as a tumor suppressive mechanism, whereas accumulative evidence demonstrate that senescent cells play an adverse role during cancer progression. The scarcity of specific and reliable markers reflecting senescence level in cancer impede our understanding of this biological basis.

OBJECTIVES

Senescence-related genes (SRGs) were collected for integrative analysis to reveal the role of senescence in hepatocellular carcinoma (HCC).

METHODS

Consensus clustering was used to subtype HCC based on SRGs. Several computational methods, including single sample gene set enrichment analysis (ssGSEA), fuzzy c-means algorithm, were performed. Data of drug sensitivities were utilized to screen potential therapeutic agents for different senescence patients. Additionally, we developed a method called signature-related gene analysis (SRGA) for identification of markers relevant to phenotype of interest. Experimental strategies consisting quantitative real-time PCR (qRT-PCR), β-galactosidase assay, western blot, and tumor-T cell co-culture system were used to validate the findings in vitro.

RESULTS

We identified three robust prognostic clusters of HCC patients with distinct survival outcome, mutational landscape, and immune features. We further extracted signature genes of senescence clusters to construct the senescence scoring system and profile senescence level in HCC at bulk and single-cell resolution. Senescence-induced stemness reprogramming was confirmed both in silico and in vitro. HCC patients with high senescence were immune suppressed and sensitive to Tozasertib and other drugs. We suggested that MAFG, PLIN3, and 4 other genes were pertinent to HCC senescence, and MAFG potentially mediated immune suppression, senescence, and stemness.

CONCLUSION

Our findings provide insights into the role of SRGs in patients stratification and precision medicine.

Association of epigenetic landscapes with heterogeneity and plasticity in pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) has a poor prognosis. Due to a lack of clear symptoms, patients often present with advanced disease, with limited clinical intervention options. The high mortality rate of PDAC is, however, also a result of several other factors that include a high degree of heterogeneity and treatment resistant cellular phenotypes. Molecular subtypes of PDAC have been identified that are thought to represent cellular phenotypes at the tissue level. The epigenetic landscape is an important factor that dictates these subtypes. Permissive epigenetic landscapes serve as drivers of molecular heterogeneity and cellular plasticity in developing crypts as well as metaplastic lesions. Drawing parallels with other cancers, we hypothesize that epigenetic permissiveness is a potential driver of cellular plasticity in PDAC. In this review will explore the epigenetic alterations that underlie PDAC cell states and relate them to cellular plasticity from other contexts. In doing so, we aim to highlight epigenomic drivers of PDAC heterogeneity and plasticity and, with that, offer some insight to guide pre-clinical research.

Robust Cluster Prediction Across Data Types Validates Association of Sex and Therapy Response in GBM

BACKGROUND

Previous studies have described sex-specific patient subtyping in glioblastoma. The cluster labels associated with these "legacy data" were used to train a predictive model capable of recapitulating this clustering in contemporary contexts.

METHODS

We used robust ensemble machine learning to train a model using gene microarray data to perform multi-platform predictions including RNA-seq and potentially scRNA-seq.

RESULTS

The engineered feature set was composed of many previously reported genes that are associated with patient prognosis. Interestingly, these well-known genes formed a predictive signature only for female patients, and the application of the predictive signature to male patients produced unexpected results.

CONCLUSIONS

This work demonstrates how annotated "legacy data" can be used to build robust predictive models capable of multi-target predictions across multiple platforms.

Comprehensive Evaluation of Multi-Omics Clustering Algorithms for Cancer Molecular Subtyping

As a highly heterogeneous and complex disease, the identification of cancer's molecular subtypes is crucial for accurate diagnosis and personalized treatment. The integration of multi-omics data enables a comprehensive interpretation of the molecular characteristics of cancer at various biological levels. In recent years, an increasing number of multi-omics clustering algorithms for cancer molecular subtyping have been proposed. However, the absence of a definitive gold standard makes it challenging to evaluate and compare these methods effectively. In this study, we developed a general framework for the comprehensive evaluation of multi-omics clustering algorithms and introduced an innovative metric, the accuracy-weighted average index, which simultaneously considers both clustering performance and clinical relevance. Using this framework, we performed a thorough evaluation and comparison of 11 state-of-the-art multi-omics clustering algorithms, including deep learning-based methods. By integrating the accuracy-weighted average index with computational efficiency, our analysis reveals that PIntMF demonstrates the best overall performance, making it a promising tool for molecular subtyping across a wide range of cancers.

Utility of Multiparametric Breast MRI Radiomics to Predict Cyclin D1 and TGF-β1 Expression

OBJECTIVE

To develop a machine learning model that integrates clinical features and multisequence MRI radiomics for noninvasively predicting the expression status of prognostic-related factors cyclin D1 and TGF-β1 in breast cancer, providing additional information for the clinical development of personalized treatment plans.

METHODS

A total of 123 breast cancer patients confirmed by surgical pathology were retrospectively enrolled in our Hospital from January 2016 to July 2022. The patients were randomly divided into a training group (87 cases) and a validation group (36 cases). Preoperative routine and dynamic contrast-enhanced magnetic resonance imaging scans of the breast were performed for treatment subjects. The region of interest was manually outlined, and texture features were extracted using AK software. Subsequently, the LASSO algorithm was employed for dimensionality reduction and feature selection to establish the MRI radiomics labels. The diagnostic efficacy and clinical value were assessed through receiver operating characteristic (ROC) curve analysis and decision curve analysis (DCA).

RESULTS

In the cyclin D1 cohort, the area under the receiver operating characteristic (ROC) curve in the clinical prediction model training and validation groups was 0.738 and 0.656, respectively. The multisequence MRI radiomics prediction model achieved an AUC of 0.874 and 0.753 in these respective groups, while the combined prediction model yielded an AUC of 0.892 and 0.785. In the TGF-β1 cohort, the ROC AUC for the clinical prediction model was found to be 0.693 and 0.645 in the training and validation groups, respectively. For the multiseries MRI radiomics prediction model, it achieved an AUC of 0.875 and 0.760 in these respective groups; whereas for the combined prediction model, it reached an AUC of 0.904 and 0.833. Decision curve analysis (DCA) demonstrated that both cohorts indicated a higher clinical application value for the combined prediction model compared with both individual models-clinical prediction model alone or radiomics model.

CONCLUSION

The integration of clinical features and multisequence MRI radiomics in a combined modeling approach holds significant predictive value for the expression status of cyclin D1 and TGF-β1. The model provides a noninvasive, dynamic evaluation method that provides effective guidance for clinical treatment.

Survival prediction and molecular subtyping of squamous cell lung cancer based on network embedding

Squamous cell lung cancer (SQCLC), which is fatal to humans, is heterogeneous with different genetic and histological features. We used SBMOI, a multi-omics data integration method from previous study, to integrate clinical, gene expression, and somatic mutation data of SQCLC to construct new patient features. Next, random survival forest (RSF) model and SimpleMKL model were constructed to predict the survival of SQCLC patients, and K-means model was constructed to perform molecular subtyping. The results of the RSF model showed that when the dimension of the patient features were 11 × 364 and the hard threshold was 0.2, we obtained the best results, and the AUC value of the 1-year time-dependent ROC curve was 0.706. The SimpleMKL model, constructed using the same patient features, performed exceptionally well, with 1-year, 5-year, and 10-year survival prediction AUC values of 0.944, 0.947 and 0.950, respectively. We used K-means analysis to identify three SQCLC molecular subtypes with significant survival differences. The patient features constructed by SBMOI were used to effectively predict the survival and molecular subtyping of SQCLC patients. In addition, our study further confirmed the effectiveness in multi-omics data integration task and broad applicability of SBMOI.

Deciphering the molecular heterogeneity of intermediate- and (very-)high-risk non-muscle-invasive bladder cancer using multi-layered -omics studies

Bladder cancer (BC) is the most common malignancy of the urinary tract. About 75% of all BC patients present with non-muscle-invasive BC (NMIBC), of which up to 70% will recur, and 15% will progress in stage and grade. As the recurrence and progression rates of NMIBC are strongly associated with some clinical and pathological factors, several risk stratification models have been developed to individually predict the short- and long-term risks of disease recurrence and progression. The NMIBC patients are stratified into four risk groups as low-, intermediate-, high-risk, and very high-risk by the European Association of Urology (EAU). Significant heterogeneity in terms of oncological outcomes and prognosis has been observed among NMIBC patients within the same EAU risk group, which has been partly attributed to the intrinsic heterogeneity of BC at the molecular level. Currently, we have a poor understanding of how to distinguish intermediate- and (very-)high-risk NMIBC with poor outcomes from those with a more benign disease course and lack predictive/prognostic tools that can specifically stratify them according to their pathologic and molecular properties. There is an unmet need for developing a more accurate scoring system that considers the treatment they receive after TURBT to enable their better stratification for further follow-up regimens and treatment selection, based also on a better response prediction to the treatment. Based on these facts, by employing a multi-layered -omics (namely, genomics, epigenetics, transcriptomics, proteomics, lipidomics, metabolomics) and immunohistopathology approach, we hypothesize to decipher molecular heterogeneity of intermediate- and (very-)high-risk NMIBC and to better stratify the patients with this disease. A combination of different -omics will provide a more detailed and multi-dimensional characterization of the tumor and represent the broad spectrum of NMIBC phenotypes, which will help to decipher the molecular heterogeneity of intermediate- and (very-)high-risk NMIBC. We think that this combinatorial multi-omics approach has the potential to improve the prediction of recurrence and progression with higher precision and to develop a molecular feature-based algorithm for stratifying the patients properly and guiding their therapeutic interventions in a personalized manner.

Spatial oncology: Translating contextual biology to the clinic

Microscopic examination of cells in their tissue context has been the driving force behind diagnostic histopathology over the past two centuries. Recently, the rise of advanced molecular biomarkers identified through single cell profiling has increased our understanding of cellular heterogeneity in cancer but have yet to significantly impact clinical care. Spatial technologies integrating molecular profiling with microenvironmental features are poised to bridge this translational gap by providing critical in situ context for understanding cellular interactions and organization. Here, we review how spatial tools have been used to study tumor ecosystems and their clinical applications. We detail findings in cell-cell interactions, microenvironment composition, and tissue remodeling for immune evasion and therapeutic resistance. Additionally, we highlight the emerging role of multi-omic spatial profiling for characterizing clinically relevant features including perineural invasion, tertiary lymphoid structures, and the tumor-stroma interface. Finally, we explore strategies for clinical integration and their augmentation of therapeutic and diagnostic approaches.

A Multi-Omics Prognostic Model Capturing Tumor Stemness and the Immune Microenvironment in Clear Cell Renal Cell Carcinoma

Background: Cancer stem-like cells (CSCs), a distinct subset recognized for their stem cell-like abilities, are intimately linked to the resistance to radiotherapy, metastatic behaviors, and self-renewal capacities in tumors. Despite their relevance, the definitive traits and importance of CSCs in the realm of oncology are still not fully comprehended, particularly in the context of clear cell renal cell carcinoma (ccRCC). A comprehensive understanding of these CSCs' properties in relation to stemness, and their impact on the efficacy of treatment and resistance to medication, is of paramount importance. Methods: In a meticulous research effort, we have identified new molecular categories designated as CRCS1 and CRCS2 through the application of an unsupervised clustering algorithm. The analysis of these subtypes included a comprehensive examination of the tumor immune environment, patterns of metabolic activity, progression of the disease, and its response to immunotherapy. In addition, we have delved into understanding these subtypes' distinctive clinical presentations, the landscape of their genomic alterations, and the likelihood of their response to various pharmacological interventions. Proceeding from these insights, prognostic models were developed that could potentially forecast the outcomes for patients with ccRCC, as well as inform strategies for the surveillance of recurrence after treatment and the handling of drug-resistant scenarios. Results: Compared with CRCS1, CRCS2 patients had a lower clinical stage/grading and a better prognosis. The CRCS2 subtype was in a hypoxic state and was characterized by suppression and exclusion of immune function, which was sensitive to gefitinib, erlotinib, and saracatinib. The constructed prognostic risk model performed well in both training and validation cohorts, helping to identify patients who may benefit from specific treatments or who are at risk of recurrence and drug resistance. A novel therapeutic target, SAA2, regulating neutrophil and fibroblast infiltration, and, thus promoting ccRCC progression, was identified. Conclusions: Our findings highlight the key role of CSCs in shaping the ccRCC tumor microenvironment, crucial for therapy research and clinical guidance. Recognizing tumor stemness helps to predict treatment efficacy, recurrence, and drug resistance, informing treatment strategies and enhancing ccRCC patient outcomes.

Pattern-centric transformation of omics data grounded on discriminative gene associations aids predictive tasks in TCGA while ensuring interpretability

The increasing prevalence of omics data sources is pushing the study of regulatory mechanisms underlying complex diseases such as cancer. However, the vast quantities of molecular features produced and the inherent interplay between them lead to a level of complexity that hampers both descriptive and predictive tasks, requiring custom-built algorithms that can extract relevant information from these sources of data. We propose a transformation that moves data centered on molecules (e.g., transcripts and proteins) to a new data space focused on putative regulatory modules given by statistically relevant co-expression patterns. To this end, the proposed transformation extracts patterns from the data through biclustering and uses them to create new variables with guarantees of interpretability and discriminative power. The transformation is shown to achieve dimensionality reductions of up to 99% and increase predictive performance of various classifiers across multiple omics layers. Results suggest that omics data transformations from gene-centric to pattern-centric data supports both prediction tasks and human interpretation, notably contributing to precision medicine applications.

Multiomics2Targets identifies targets from cancer cohorts profiled with transcriptomics, proteomics, and phosphoproteomics

The availability of data from profiling of cancer patients with multiomics is rapidly increasing. However, integrative analysis of such data for personalized target identification is not trivial. Multiomics2Targets is a platform that enables users to upload transcriptomics, proteomics, and phosphoproteomics data matrices collected from the same cohort of cancer patients. After uploading the data, Multiomics2Targets produces a report that resembles a research publication. The uploaded matrices are processed, analyzed, and visualized using the tools Enrichr, KEA3, ChEA3, Expression2Kinases, and TargetRanger to identify and prioritize proteins, genes, and transcripts as potential targets. Figures and tables, as well as descriptions of the methods and results, are automatically generated. Reports include an abstract, introduction, methods, results, discussion, conclusions, and references and are exportable as citable PDFs and Jupyter Notebooks. Multiomics2Targets is applied to analyze version 3 of the Clinical Proteomic Tumor Analysis Consortium (CPTAC3) pan-cancer cohort, identifying potential targets for each CPTAC3 cancer subtype. Multiomics2Targets is available from https://multiomics2targets.maayanlab.cloud/.

Comprehensive Transcriptome Analysis Expands lncRNA Functional Profiles in Breast Cancer

Breast cancer is a heterogeneous disease that arises as a multi-stage process involving multiple cell types. Patients diagnosed with the same clinical stage and pathological classification may have different prognoses and therapeutic responses due to alterations in molecular genetics. As an essential marker for the molecular subtyping of breast cancer, long non-coding RNAs (lncRNAs) play a crucial role in gene expression regulation, cell differentiation, and the maintenance of genomic stability. Here, we developed a modular framework for lncRNA identification and applied it to a breast cancer cohort to identify novel lncRNAs not previously annotated. To investigate the potential biological function, regulatory mechanisms, and clinical relevance of the novel lncRNAs, we elucidated the genomic and chromatin features of these lncRNAs, along with the associated protein-coding genes and putative enhancers involved in the breast cancer regulatory networks. Furthermore, we uncovered that the expression patterns of novel and annotated lncRNAs identified in breast cancer were related to the hormone response in the PAM50 subtyping criterion, as well as the immune response and progression states of breast cancer across different immune cells and immune checkpoint genes. Collectively, the comprehensive identification and functional analysis of lncRNAs revealed that these lncRNAs play an essential role in breast cancer by altering gene expression and participating in the regulatory networks, contributing to a better insight into breast cancer heterogeneity and potential avenues for therapeutic intervention.

Tumor Subtype Classification Tool for HPV-associated Head and Neck Cancers

Importance

Molecular subtypes of HPV-associated Head and Neck Squamous Cell Carcinoma (HNSCC), named IMU (immune strong) and KRT (highly keratinized), are well-recognized and have been shown to have distinct mechanisms of carcinogenesis, clinical outcomes, and potentially differing optimal treatment strategies. Currently, no standardized method exists to subtype a new HPV+ HNSCC tumor. Our paper introduces a machine learning-based classifier and webtool to reliably subtype HPV+ HNSCC tumors using the IMU/KRT paradigm and highlights the importance of subtype in HPV+ HNSCC.

Objective

To develop a robust, accurate machine learning-based classification tool that standardizes the process of subtyping HPV+ HNSCC, and to investigate the clinical, demographic, and molecular features associated with subtype in a meta-analysis of four patient cohorts.

Data Sources

We conducted RNA-seq on 67 HNSCC FFPE blocks from University of Michigan hospital. Combining this with three publicly available datasets, we utilized a total of 229 HPV+ HNSCC RNA-seq samples. All participants were HPV+ according to RNA expression. An ensemble machine learning approach with five algorithms and three different input training gene sets were developed, with final subtype determined by majority vote. Several additional steps were taken to ensure rigor and reproducibility throughout.

Study Selection

The classifier was trained and tested using 84 subtype-labeled HPV+ RNA-seq samples from two cohorts: University of Michigan (UM; n=18) and TCGA-HNC (n=66). The classifier robustness was validated with two independent cohorts: 83 samples from the HPV Virome Consortium and 62 additional samples from UM. We revealed 24 of 39 tested clinicodemographic and molecular variables significantly associated with subtype.

Results

The classifier achieved 100% accuracy in the test set. Validation on two additional cohorts demonstrated successful separation by known features of the subtypes. Investigating the relationship between subtype and 39 molecular and clinicodemographic variables revealed IMU is associated with epithelial-mesenchymal transition (p=2.25×10-4), various immune cell types, and lower radiation resistance (p=0.0050), while KRT is more highly keratinized (p=2.53×10-8), and more likely female than IMU (p=0.0082).

Conclusions and Relevance

This study provides a reliable classifier for subtyping HPV+ HNSCC tumors as either IMU or KRT based on bulk RNA-seq data, and additionally, improves our understanding of the HPV+ HNSCC subtypes.

The "Road" to Malignant Transformation from Endometriosis to Endometriosis-Associated Ovarian Cancers (EAOCs): An mTOR-Centred Review

Ovarian cancer is an umbrella term covering a number of distinct subtypes. Endometrioid and clear-cell ovarian carcinoma are endometriosis-associated ovarian cancers (EAOCs) frequently arising from ectopic endometrium in the ovary. The mechanistic target of rapamycin (mTOR) is a crucial regulator of cellular homeostasis and is dysregulated in both endometriosis and endometriosis-associated ovarian cancer, potentially favouring carcinogenesis across a spectrum from benign disease with cancer-like characteristics, through an atypical phase, to frank malignancy. In this review, we focus on mTOR dysregulation in endometriosis and EAOCs, investigating cancer driver gene mutations and their potential interaction with the mTOR pathway. Additionally, we explore the complex pathogenesis of transformation, considering environmental, hormonal, and epigenetic factors. We then discuss postmenopausal endometriosis pathogenesis and propensity for malignant transformation. Finally, we summarize the current advancements in mTOR-targeted therapeutics for endometriosis and EAOCs.

Subtype-MGTP: a cancer subtype identification framework based on multi-omics translation

MOTIVATION

The identification of cancer subtypes plays a crucial role in cancer research and treatment. With the rapid development of high-throughput sequencing technologies, there has been an exponential accumulation of cancer multi-omics data. Integrating multi-omics data has emerged as a cost-effective and efficient strategy for cancer subtyping. While current methods primarily rely on genomics data, protein expression data offers a closer representation of phenotype. Therefore, integrating protein expression data holds promise for enhancing subtyping accuracy. However, the scarcity of protein expression data compared to genomics data presents a challenge in its direct incorporation into existing methods. Moreover, striking a balance between omics-specific learning and cross-omics learning remains a prevalent challenge in current multi-omics integration methods.

RESULTS

We introduce Subtype-MGTP, a novel cancer subtyping framework based on the translation of Multiple Genomics To Proteomics. Subtype-MGTP comprises two modules: a translation module, which leverages available protein data to translate multi-type genomics data into predicted protein expression data, and an improved deep subspace clustering module, which integrates contrastive learning to cluster the predicted protein data, yielding refined subtyping results. Extensive experiments conducted on benchmark datasets demonstrate that Subtype-MGTP outperforms nine state-of-the-art cancer subtyping methods. The interpretability of clustering results is further supported by the clinical and survival analysis. Subtype-MGTP also exhibits strong robustness against varying rates of missing protein data and demonstrates distinct advantages in integrating multi-omics data with imbalanced multi-omics data.

AVAILABILITY AND IMPLEMENTATION

The code and results are available at https://github.com/kybinn/Subtype-MGTP.

Omics-based molecular classifications empowering in precision oncology

BACKGROUND

In the past decades, cancer enigmatical heterogeneity at distinct expression levels could interpret disparities in therapeutic response and prognosis. It built hindrances to precision medicine, a tactic to tailor customized treatment informed by the tumors' molecular profile. Single-omics analysis dissected the biological features associated with carcinogenesis to some extent but still failed to revolutionize cancer treatment as expected. Integrated omics analysis incorporated tumor biological networks from diverse layers and deciphered a holistic overview of cancer behaviors, yielding precise molecular classification to facilitate the evolution and refinement of precision medicine.

CONCLUSION

This review outlined the biomarkers at multiple expression layers to tutor molecular classification and pinpoint tumor diagnosis, and explored the paradigm shift in precision therapy: from single- to multi-omics-based subtyping to optimize therapeutic regimens. Ultimately, we firmly believe that by parsing molecular characteristics, omics-based typing will be a powerful assistant for precision oncology.

Themis: advancing precision oncology through comprehensive molecular subtyping and optimization

Recent advances in tumor molecular subtyping have revolutionized precision oncology, offering novel avenues for patient-specific treatment strategies. However, a comprehensive and independent comparison of these subtyping methodologies remains unexplored. This study introduces 'Themis' (Tumor HEterogeneity analysis on Molecular subtypIng System), an evaluation platform that encapsulates a few representative tumor molecular subtyping methods, including Stemness, Anoikis, Metabolism, and pathway-based classifications, utilizing 38 test datasets curated from The Cancer Genome Atlas (TCGA) and significant studies. Our self-designed quantitative analysis uncovers the relative strengths, limitations, and applicability of each method in different clinical contexts. Crucially, Themis serves as a vital tool in identifying the most appropriate subtyping methods for specific clinical scenarios. It also guides fine-tuning existing subtyping methods to achieve more accurate phenotype-associated results. To demonstrate the practical utility, we apply Themis to a breast cancer dataset, showcasing its efficacy in selecting the most suitable subtyping methods for personalized medicine in various clinical scenarios. This study bridges a crucial gap in cancer research and lays a foundation for future advancements in individualized cancer therapy and patient management.

Identification and classification of glioma subtypes based on RNA-binding proteins

BACKGROUND

Glioma is a common and aggressive primary malignant cancer known for its high morbidity, mortality, and recurrence rates. Despite this, treatment options for glioma are currently restricted. The dysregulation of RBPs has been linked to the advancement of several types of cancer, but their precise role in glioma evolution is still not fully understood. This study sought to investigate how RBPs may impact the development and prognosis of glioma, with potential implications for prognosis and therapy.

METHODS

RNA-seq profiles of glioma and corresponding clinical data from the CGGA database were initially collected for analysis. Unsupervised clustering was utilized to identify crucial tumor subtypes in glioma development. Subsequent time-series analysis and MS model were employed to track the progression of these identified subtypes. RBPs playing a significant role in glioma progression were then pinpointed using WGCNA and Lasso Cox regression models. Functional analysis of these key RBP-related genes was conducted through GSEA. Additionally, the CIBERSORT algorithm was utilized to estimate immune infiltrating cells, while the STRING database was consulted to uncover potential mechanisms of the identified biomarkers.

RESULTS

Six tumor subgroups were identified and found to be highly homogeneous within each subgroup. The progression stages of these tumor subgroups were determined using time-series analysis and a MS model. Through WGCNA, Lasso Cox, and multivariate Cox regression analysis, it was confirmed that BCLAF1 is correlated with survival in glioma patients and is closely linked to glioma progression. Functional annotation suggests that BCLAF1 may impact glioma progression by influencing RNA splicing, which in turn affects the cell cycle, Wnt signaling pathway, and other cancer development pathways.

CONCLUSIONS

The study initially identified six subtypes of glioma progression and assessed their malignancy ranking. Furthermore, it was determined that BCLAF1 could serve as an RBP-related prognostic marker, offering significant implications for the clinical diagnosis and personalized treatment of glioma.

A pathway-based computational framework for identification of a new modal of multi-omics biomarkers and its application in esophageal cancer

BACKGROUND

The pathway-based strategy has been recently proposed for identifying biomarkers with the advantages of higher biological interpretability and cross-data robustness than the conventional gene-based strategy. However, its utility in clinical applications has been limited due to the high computational complexity and ill-defined performance.

OBJECTIVE

The current study presents a machine learning-based computational framework using multi-omics data for identifying a new modal of biomarkers, called pathway-derived core biomarkers, which have the advantages of both gene-based and pathway-based biomarkers.

METHODS

Machine-learning methods and gene-pathway network were integrated to select the pathway-derived core biomarkers. Multiple machine-learning algorithms were used to construct and validate the diagnostic models of the biomarkers based on more than 1400 multi-omics clinical samples of esophageal squamous cell carcinoma (ESCC).

RESULTS

The results showed that the classifier models based on the new modal biomarkers achieved superior performance in the training datasets with an average AUC/accuracy of 0.98/0.95 and 0.89/0.81 for mRNAs and miRNA, respectively, higher than the currently known classifier models based on the conventional gene-based strategy and pathway-based strategy. In the testing cohorts, the AUC/accuracy increased by 6.1 %/7.3 % than the models based on the native gene-based biomarkers. The improved performance was further confirmed in independent validation cohorts. Specifically, the sensitivity/specificity increased by ∼3 % and the variance significantly decreased by ∼69 % compared with that of the native gene-based biomarkers. Importantly, the pathway-derived core biomarkers also recovered 45 % more previously reported biomarkers than the gene-based biomarkers and are more functionally relevant to the ESCC etiology (involved in 14 versus 7 pathways related with ESCC or other cancer), highlighting the cross-data robustness of this new modal of biomarkers via enhanced functional relevance.

CONCLUSIONS

The results demonstrated that the new modal of biomarkers not only have improved predicting performance and robustness, but also exhibit higher functional interpretability thus leading to the potential application in cancer diagnosis.

GraphPath: a graph attention model for molecular stratification with interpretability based on the pathway-pathway interaction network

MOTIVATION

Studying the molecular heterogeneity of cancer is essential for achieving personalized therapy. At the same time, understanding the biological processes that drive cancer development can lead to the identification of valuable therapeutic targets. Therefore, achieving accurate and interpretable clinical predictions requires paramount attention to thoroughly characterizing patients at both the molecular and biological pathway levels.

RESULTS

Here, we present GraphPath, a biological knowledge-driven graph neural network with multi-head self-attention mechanism that implements the pathway-pathway interaction network. We train GraphPath to classify the cancer status of patients with prostate cancer based on their multi-omics profiling. Experiment results show that our method outperforms P-NET and other baseline methods. Besides, two external cohorts are used to validate that the model can be generalized to unseen samples with adequate predictive performance. We reduce the dimensionality of latent pathway embeddings and visualize corresponding classes to further demonstrate the optimal performance of the model. Additionally, since GraphPath's predictions are interpretable, we identify target cancer-associated pathways that significantly contribute to the model's predictions. Such a robust and interpretable model has the potential to greatly enhance our understanding of cancer's biological mechanisms and accelerate the development of targeted therapies.

AVAILABILITY AND IMPLEMENTATION

https://github.com/amazingma/GraphPath.

Deep learning in cancer genomics and histopathology

Histopathology and genomic profiling are cornerstones of precision oncology and are routinely obtained for patients with cancer. Traditionally, histopathology slides are manually reviewed by highly trained pathologists. Genomic data, on the other hand, is evaluated by engineered computational pipelines. In both applications, the advent of modern artificial intelligence methods, specifically machine learning (ML) and deep learning (DL), have opened up a fundamentally new way of extracting actionable insights from raw data, which could augment and potentially replace some aspects of traditional evaluation workflows. In this review, we summarize current and emerging applications of DL in histopathology and genomics, including basic diagnostic as well as advanced prognostic tasks. Based on a growing body of evidence, we suggest that DL could be the groundwork for a new kind of workflow in oncology and cancer research. However, we also point out that DL models can have biases and other flaws that users in healthcare and research need to know about, and we propose ways to address them.

USP3 promotes osteosarcoma progression via deubiquitinating EPHA2 and activating the PI3K/AKT signaling pathway

Ubiquitin-specific protease 3 (USP3) plays an important role in the progression of various tumors. However, the role of USP3 in osteosarcoma (OS) remains poorly understood. The aim of this study was to explore the biological function of USP3 in OS and the underlying molecular mechanism. We found that OS had higher USP3 expression compared with that of normal bone tissue, and high expression of USP3 was associated with poor prognosis in patients with OS. Overexpression of USP3 significantly increased OS cell proliferation, migration, and invasion. Mechanistically, USP3 led to the activation of the PI3K/AKT signaling pathway in OS by binding to EPHA2 and then reducing its protein degradation. Notably, the truncation mutant USP3-F2 (159-520) interacted with EPHA2, and amino acid 203 was found to play an important role in this process. And knockdown of EPHA2 expression reversed the pro-tumour effects of USP3-upregulating. Thus, our study indicates the USP3/EPHA2 axis may be a novel potential target for OS treatment.

A signature of enhanced proliferation associated with response and survival to anti-PD-L1 therapy in early-stage non-small cell lung cancer

In early-stage non-small cell lung cancer, the combination of neoadjuvant anti-PD-L1 and subablative stereotactic body radiation therapy (SBRT) is associated with higher rates of major pathologic response compared to anti-PD-L1 alone. Here, we identify a 140-gene set, enriched in genes characteristic of highly proliferating cells, associated with response to the dual therapy. Analysis of on-treatment transcriptome data indicate roles for T and B cells in response. The 140-gene set is associated with disease-free survival when applied to the combined trial arms. This 140-gene set identifies a subclass of tumors in all 7 of The Cancer Genome Atlas tumor types examined. Worse survival is associated with the 140-gene signature in 5 of these tumor types. Collectively, our data support that this 140-gene set, discovered in association with response to combined anti-PD-L1 and SBRT, identifies a clinically aggressive subclass of solid tumors that may be more likely to respond to immunotherapies.

Tumor molecular landscape of Epstein-Barr virus (EBV) related nasopharyngeal carcinoma in EBV-endemic and non-endemic areas: Implications for improving treatment modalities

Epstein-Barr virus (EBV) related- nasopharyngeal carcinoma (NPC) is a squamous carcinoma of the nasopharyngeal mucosal lining. Endemic areas (EA) are east and Southeast Asia, were NPC was recorded with higher incidence and longer estimated survival than in non-endemic area (NEA) such as Europe, We analyzed the gene expression and microenvironment properties of NPC in both areas to identify molecular subtypes and assess biological and clinical correlates that might explain the differences in incidence and outcome between EA- and NEA-NPCs. Six EA-NPC transcriptomic datasets, including tumor and normal samples, were integrated in a meta-analysis to identify molecular subtypes using a ConsensusClusterPlus bioinformatic approach. Based on the biological/functional characterization of four identified clusters were identified: Cl1, Immune-active; Cl2, defense-response; Cl3, proliferation; Cl4, perineural-interaction/EBV-exhaustion. Kaplan-Meier survival analysis, applied to the single dataset with available disease-free survival indicated Cl3 as the cluster with the worst prognosis (P = 0.0476), confirmed when applying four previously published prognostic signatures. A Cl3 classifier signature was generated and its prognostic performance was confirmed (P = 0.0368) on a validation dataset. Prediction of treatment response suggested better responses to: radiotherapy and immune checkpoint inhibitors immune-active and defense-response clusters; chemotherapy proliferation cluster; cisplatin perineural-interaction/EBV-exhaustion cluster. RNA sequencing for gene expression profiling was performed on 50 NEA-NPC Italian samples. In the NEA cohort, Cl1, Cl2 and Cl3 were represented, while perineural-interaction/EBV-exhaustion was almost absent. The immune/biological characterization and treatment-response prediction analyses of NEA-NPC partially replicated the EA-NPC results. Well characterized EA- and NEA-NPC retrospective and prospective cohorts are needed to validate the obtained results and can help designing future clinical studies.

A deep-learning workflow to predict upper tract urothelial carcinoma protein-based subtypes from H&E slides supporting the prioritization of patients for molecular testing

Upper tract urothelial carcinoma (UTUC) is a rare and aggressive, yet understudied, urothelial carcinoma (UC). The more frequent UC of the bladder comprises several molecular subtypes, associated with different targeted therapies and overlapping with protein-based subtypes. However, if and how these findings extend to UTUC remains unclear. Artificial intelligence-based approaches could help elucidate UTUC's biology and extend access to targeted treatments to a wider patient audience. Here, UTUC protein-based subtypes were identified, and a deep-learning (DL) workflow was developed to predict them directly from routine histopathological H&E slides. Protein-based subtypes in a retrospective cohort of 163 invasive tumors were assigned by hierarchical clustering of the immunohistochemical expression of three luminal (FOXA1, GATA3, and CK20) and three basal (CD44, CK5, and CK14) markers. Cluster analysis identified distinctive luminal (N = 80) and basal (N = 42) subtypes. The luminal subtype mostly included pushing, papillary tumors, whereas the basal subtype diffusely infiltrating, non-papillary tumors. DL model building relied on a transfer-learning approach by fine-tuning a pre-trained ResNet50. Classification performance was measured via three-fold repeated cross-validation. A mean area under the receiver operating characteristic curve of 0.83 (95% CI: 0.67-0.99), 0.8 (95% CI: 0.62-0.99), and 0.81 (95% CI: 0.65-0.96) was reached in the three repetitions. High-confidence DL-based predicted subtypes showed significant associations (p < 0.001) with morphological features, i.e. tumor type, histological subtypes, and infiltration type. Furthermore, a significant association was found with programmed cell death ligand 1 (PD-L1) combined positive score (p < 0.001) and FGFR3 mutational status (p = 0.002), with high-confidence basal predictions containing a higher proportion of PD-L1 positive samples and high-confidence luminal predictions a higher proportion of FGFR3-mutated samples. Testing of the DL model on an independent cohort highlighted the importance to accommodate histological subtypes. Taken together, our DL workflow can predict protein-based UTUC subtypes, associated with the presence of targetable alterations, directly from H&E slides.

Trends in orphan medicinal products approvals in the European Union between 2010-2022

BACKGROUND

Over the last twenty years of orphan drug regulation in Europe, the regulatory framework has increased its complexity, with different regulatory paths and tools engineered to facilitate the innovation and accelerate approvals. Recently, the proposal of the new Pharmaceutical Legislation for the European Union, which will replace at least three Regulations and one Directive, was released and its new framework is raising many questions. The aim of this study was to present a characterisation of the Orphan Medicinal Products (OMPs) authorised by the European Commission (EC), between 2010 and 2022, looking into eighteen variables, contributing to the ongoing discussion on the proposal and implementation of the new Pharmaceutical Legislation proposed.

METHODS

Data of the OMPs identified and approved between 2010 and 2022 were extracted from the European Public Assessment Reports (EPARs) produced by the European Medicines Agency. Information regarding legal basis of the application, applicant, protocol assistance received, type of authorization, registration status, type of molecule, ATC code, therapeutic area, target age, disease prevalence, number of pivotal clinical trials supporting the application, clinical trial designs, respective efficacy endpoints and number of patients enrolled in the pivotal clinical trials were extracted. A descriptive statistical analysis was applied.

RESULTS

We identified 192 OMPs approved in the period between 2010 and 2022. 89% of the OMPs have legal basis of "full application". 86% of the sponsors received protocol assistance whereas 64% of the MAA benefited from the accelerated assessment. 53% of the active substances are small molecules; about 1 in 5 molecules are repurposed. 40% of the OMPs have oncological therapeutic indications and 56% of the OMPs are intended to treat only adults. 71% of the products were approved based on a single pivotal trial.

CONCLUSIONS

This analysis of OMPs approved between 2010 and 2022 shows that a shift has occurred in the rare disease medicine development space. Through the period studied we observe an increase of non-small molecules approved, accelerated assessment received and non-standard MA's granted.

Autoimmune inflammatory myopathy biomarkers

The autoimmune inflammatory myopathy disease spectrum, commonly known as myositis, is a group of systemic diseases that mainly affect the muscles, skin and lungs. Biomarker assessment helps in understanding disease mechanisms, allowing for the implementation of precise strategies in the classification, diagnosis, and management of these diseases. This review examines the pathogenic mechanisms and highlights current data on blood and tissue biomarkers of autoimmune inflammatory myopathies.

YRNA and tRNA fragments can differentiate benign from malignant canine mammary gland tumors

Mammary gland tumors (MGT) are the most common tumors in sexually intact female dogs. The functional regulation of miRNAs, a type of noncoding RNAs (ncRNAs), in canine MGT has been extensively investigated. However, the expression of other ncRNAs, such as YRNAs and transfer RNA-derived fragments (tRFs) in canine MGT is unknown. We investigated ncRNAs other than miRNAs from our small RNA project (PRJNA716131) in different canine MGT histologic subtypes. This study included benign tumors (benign mixed tumor, complex adenoma) and malignant tumors (carcinoma in benign tumor and carcinoma with metastasis) samples. Aberrantly expressed ncRNAs were examined by comparisons among MGT subtypes. The relative expression trends were validated in canine MGT tissues, plasma, extracellular vesicles, and MGT cell lines using quantitative reverse transcription PCR. Three aberrantly expressed ncRNAs were identified by comparisons among MGT subtypes. YRNA and tRNA-Gly-GCC distinguished benign mixed tumor from other MGT histologic subtypes, while tRNA-Val differentiated complex adenoma, carcinoma in benign tumors, and carcinoma with metastasis. The ROC curve of the three ncRNAs showed they might be potential biomarkers to discriminate malignant from benign MGT. YRNA and tRFs expression levels were decreased in metastatic compared with primary canine MGT cell lines. To the best of our knowledge, this is the first investigation of YRNA and tRFs in canine MGT. The three identified ncRNAs may be biomarkers for differentiating MGT histologic subtypes. Suggested Reviewers: Powered by Editorial Manager® and ProduXion Manager® from Aries Systems Corporatio.

The use of SP/Neurokinin-1 as a Therapeutic Target in Colon and Rectal Cancer

Different studies have highlighted the role of Substance P / Neurokinin 1 Receptor (SP/NK-1R) axis in multiple hallmarks of cancer including cell transformation, proliferation, and migration as well as angiogenesis and metastasis of a wide range of solid tumors including colorectal cancer. Until now, the selective high-affinity antagonist of human SP/NK1-R aprepitant (Emend) has been authorized by the Food and Drug Administration as a low dosage medication to manage and treat chemotherapy-induced nausea. However, increasing evidence in recent years support the potential utility of high doses of aprepitant as an antitumor agent and thus, opening the possibility to the pharmacological repositioning of SP/NK1-R antagonists as an adjuvant therapy to conventional cancer treatments. In this review, we summarize current knowledge on the molecular basis of colorectal cancer as well as the pathophysiological importance of SP/NK1-R and the potential utility of SP/NK-1R axis as a therapeutic target in this malignancy.

A comparative study of clustering methods on gene expression data for lung cancer prognosis

Lung cancer subtyping based on gene expression data is important for identifying patient subgroups with differing survival prognosis to facilitate customized treatment strategies for each subtype of patients. Unsupervised clustering methods are the traditional approach for clustering patients into subtypes. However, since those methods cluster patients based only on gene expression data, the resulting clusters may not always be relevant to the survival outcome of interest. In recent years, semi-supervised and supervised methods have been proposed, which leverage the survival outcome data to identify clusters more relevant to survival prognosis. This paper aims to compare the performance of different clustering methods for identifying clinically prognostic lung cancer subtypes based on two lung adenocarcinoma datasets. For each method, we clustered patients into two clusters and assessed the difference in patient survival time between clusters. Unsupervised methods were found to have large logrank p-values and no significant results in most cases. Semi-supervised and supervised methods had improved performance over unsupervised methods and very significant p-values. These results indicate that unsupervised methods are not capable of identifying clusters with significant differences in survival prognosis in most cases, while supervised and semi-supervised methods can better cluster patients into clinically useful subtypes.

Single-cell characterization of macrophages in uveal melanoma uncovers transcriptionally heterogeneous subsets conferring poor prognosis and aggressive behavior

Uveal melanoma (UM) is the most frequent primary intraocular malignancy with high metastatic potential and poor prognosis. Macrophages represent one of the most abundant infiltrating immune cells with diverse functions in cancers. However, the cellular heterogeneity and functional diversity of macrophages in UM remain largely unexplored. In this study, we analyzed 63,264 single-cell transcriptomes from 11 UM patients and identified four transcriptionally distinct macrophage subsets (termed MΦ-C1 to MΦ-C4). Among them, we found that MΦ-C4 exhibited relatively low expression of both M1 and M2 signature genes, loss of inflammatory pathways and antigen presentation, instead demonstrating enhanced signaling for proliferation, mitochondrial functions and metabolism. We quantified the infiltration abundance of MΦ-C4 from single-cell and bulk transcriptomes across five cohorts and found that increased MΦ-C4 infiltration was relevant to aggressive behaviors and may serve as an independent prognostic indicator for poor outcomes. We propose a novel subtyping scheme based on macrophages by integrating the transcriptional signatures of MΦ-C4 and machine learning to stratify patients into MΦ-C4-enriched or MΦ-C4-depleted subtypes. These two subtypes showed significantly different clinical outcomes and were validated through bulk RNA sequencing and immunofluorescence assays in both public multicenter cohorts and our in-house cohort. Following further translational investigation, our findings highlight a potential therapeutic strategy of targeting macrophage subsets to control metastatic disease and consistently improve the outcome of patients with UM.

Identification of a prognostic evaluator from glutamine metabolic heterogeneity studies within and between tissues in hepatocellular carcinoma

Background: The liver is the major metabolic organ of the human body, and abnormal metabolism is the main factor influencing hepatocellular carcinoma (HCC). This study was designed to determine the effect of glutamine metabolism on HCC heterogeneity and to develop a prognostic evaluator based on the heterogeneity study of glutamine metabolism within HCC tumors and between tissues. Methods: Single-cell transcriptome data were extracted from the GSE149614 dataset and processed using the Seurat package in R for quality control of these data. HCC subtypes in the Cancer Genome Atlas and the GSE14520 dataset were identified via consensus clustering based on glutamine family amino acid metabolism (GFAAM) process genes. The machine learning algorithms gradient boosting machine, support vector machine, random forest, eXtreme gradient boosting, decision trees, and least absolute shrinkage and selection operator were utilized to develop the prognosis model of differentially expressed genes among the molecular gene subtypes. Results: The samples in the GSE149614 dataset included 10 cell types, and there was no significant difference in the GFAAM pathway. HCC was classified into three molecular subtypes according to GFAAM process genes, showing molecular heterogeneity in prognosis, clinicopathological features, and immune cell infiltration. C1 showed the worst survival rate and the highest immune score and immune cell infiltration. A six-gene model for prognostic and immunotherapy responses was constructed among subtypes, and the calculated high-risk score was significantly correlated with poor prognosis, high immune abundance, and a low response rate of immunotherapy in HCC. Conclusion: Our discovery of GFAAM-associated marker genes may help to further decipher the role in HCC occurrence and progression. In particular, this six-gene prognostic model may serve as a predictor of treatment and prognosis in HCC patients.

Artificial intelligence (AI) and machine learning (ML) in precision oncology: a review on enhancing discoverability through multiomics integration

Multiomics data including imaging radiomics and various types of molecular biomarkers have been increasingly investigated for better diagnosis and therapy in the era of precision oncology. Artificial intelligence (AI) including machine learning (ML) and deep learning (DL) techniques combined with the exponential growth of multiomics data may have great potential to revolutionize cancer subtyping, risk stratification, prognostication, prediction and clinical decision-making. In this article, we first present different categories of multiomics data and their roles in diagnosis and therapy. Second, AI-based data fusion methods and modeling methods as well as different validation schemes are illustrated. Third, the applications and examples of multiomics research in oncology are demonstrated. Finally, the challenges regarding the heterogeneity data set, availability of omics data, and validation of the research are discussed. The transition of multiomics research to real clinics still requires consistent efforts in standardizing omics data collection and analysis, building computational infrastructure for data sharing and storing, developing advanced methods to improve data fusion and interpretability, and ultimately, conducting large-scale prospective clinical trials to fill the gap between study findings and clinical benefits.

Artificial intelligence & clinical nutrition: What the future might have in store

Artificial Intelligence (AI) is a rapidly emerging technology in healthcare that has the potential to revolutionise clinical nutrition. AI can assist in analysing complex data, interpreting medical images, and providing personalised nutrition interventions for patients. Clinical nutrition is a critical aspect of patient care, and AI can help clinicians make more informed decisions regarding patients' nutritional requirements, disease prevention, and management. AI algorithms can analyse large datasets to identify novel associations between diet and disease outcomes, enabling clinicians to make evidence-based nutritional recommendations. AI-powered devices and applications can also assist in tracking dietary intake, providing feedback, and motivating patients towards healthier food choices. However, the adoption of AI in clinical nutrition raises several ethical and regulatory concerns, such as data privacy and bias. Further research is needed to assess the clinical effectiveness and safety of AI-powered nutrition interventions. In conclusion, AI has the potential to transform clinical nutrition, but its integration into clinical practice should be carefully monitored to ensure patient safety and benefit. This article discusses the current and future applications of AI in clinical nutrition and highlights its potential benefits.

Pan-cancer characterization of ncRNA synergistic competition uncovers potential carcinogenic biomarkers

Non-coding RNAs (ncRNAs) act as important modulators of gene expression and they have been confirmed to play critical roles in the physiology and development of malignant tumors. Understanding the synergism of multiple ncRNAs in competing endogenous RNA (ceRNA) regulation can provide important insights into the mechanisms of malignant tumors caused by ncRNA regulation. In this work, we present a framework, SCOM, for identifying ncRNA synergistic competition. We systematically construct the landscape of ncRNA synergistic competition across 31 malignant tumors, and reveal that malignant tumors tend to share hub ncRNAs rather than the ncRNA interactions involved in the synergistic competition. In addition, the synergistic competition ncRNAs (i.e. ncRNAs involved in the synergistic competition) are likely to be involved in drug resistance, contribute to distinguishing molecular subtypes of malignant tumors, and participate in immune regulation. Furthermore, SCOM can help to infer ncRNA synergistic competition across malignant tumors and uncover potential diagnostic and prognostic biomarkers of malignant tumors. Altogether, the SCOM framework (https://github.com/zhangjunpeng411/SCOM/) and the resulting web-based database SCOMdb (https://comblab.cn/SCOMdb/) serve as a useful resource for exploring ncRNA regulation and to accelerate the identification of carcinogenic biomarkers.

A workflow to study mechanistic indicators for driver gene prediction with Moonlight

Prediction of driver genes (tumor suppressors and oncogenes) is an essential step in understanding cancer development and discovering potential novel treatments. We recently proposed Moonlight as a bioinformatics framework to predict driver genes and analyze them in a system-biology-oriented manner based on -omics integration. Moonlight uses gene expression as a primary data source and combines it with patterns related to cancer hallmarks and regulatory networks to identify oncogenic mediators. Once the oncogenic mediators are identified, it is important to include extra levels of evidence, called mechanistic indicators, to identify driver genes and to link the observed gene expression changes to the underlying alteration that promotes them. Such a mechanistic indicator could be for example a mutation in the regulatory regions for the candidate gene. Here, we developed new functionalities and released Moonlight2 to provide the user with a mutation-based mechanistic indicator as a second layer of evidence. These functionalities analyze mutations in a cancer cohort to classify them into driver and passenger mutations. Those oncogenic mediators with at least one driver mutation are retained as the final set of driver genes. We applied Moonlight2 to the basal-like breast cancer subtype, lung adenocarcinoma and thyroid carcinoma using data from The Cancer Genome Atlas. For example, in basal-like breast cancer, we found four oncogenes (COPZ2, SF3B4, KRTCAP2 and POLR2J) and nine tumor suppressor genes (KIR2DL4, KIF26B, ARL15, ARHGAP25, EMCN, GMFG, TPK1, NR5A2 and TEK) containing a driver mutation in their promoter region, possibly explaining their deregulation. Moonlight2R is available at https://github.com/ELELAB/Moonlight2R.

Deep multi-view contrastive learning for cancer subtype identification

Cancer heterogeneity has posed great challenges in exploring precise therapeutic strategies for cancer treatment. The identification of cancer subtypes aims to detect patients with distinct molecular profiles and thus could provide new clues on effective clinical therapies. While great efforts have been made, it remains challenging to develop powerful computational methods that can efficiently integrate multi-omics datasets for the task. In this paper, we propose a novel self-supervised learning model called Deep Multi-view Contrastive Learning (DMCL) for cancer subtype identification. Specifically, by incorporating the reconstruction loss, contrastive loss and clustering loss into a unified framework, our model simultaneously encodes the sample discriminative information into the extracted feature representations and well preserves the sample cluster structures in the embedded space. Moreover, DMCL is an end-to-end framework where the cancer subtypes could be directly obtained from the model outputs. We compare DMCL with eight alternatives ranging from classic cancer subtype identification methods to recently developed state-of-the-art systems on 10 widely used cancer multi-omics datasets as well as an integrated dataset, and the experimental results validate the superior performance of our method. We further conduct a case study on liver cancer and the analysis results indicate that different subtypes might have different responses to the selected chemotherapeutic drugs.

Deep multi-omics integration by learning correlation-maximizing representation identifies prognostically stratified cancer subtypes

Motivation

Molecular subtyping by integrative modeling of multi-omics and clinical data can help the identification of robust and clinically actionable disease subgroups; an essential step in developing precision medicine approaches.

Results

We developed a novel outcome-guided molecular subgrouping framework, called Deep Multi-Omics Integrative Subtyping by Maximizing Correlation (DeepMOIS-MC), for integrative learning from multi-omics data by maximizing correlation between all input -omics views. DeepMOIS-MC consists of two parts: clustering and classification. In the clustering part, the preprocessed high-dimensional multi-omics views are input into two-layer fully connected neural networks. The outputs of individual networks are subjected to Generalized Canonical Correlation Analysis loss to learn the shared representation. Next, the learned representation is filtered by a regression model to select features that are related to a covariate clinical variable, for example, a survival/outcome. The filtered features are used for clustering to determine the optimal cluster assignments. In the classification stage, the original feature matrix of one of the -omics view is scaled and discretized based on equal frequency binning, and then subjected to feature selection using RandomForest. Using these selected features, classification models (for example, XGBoost model) are built to predict the molecular subgroups that were identified at clustering stage. We applied DeepMOIS-MC on lung and liver cancers, using TCGA datasets. In comparative analysis, we found that DeepMOIS-MC outperformed traditional approaches in patient stratification. Finally, we validated the robustness and generalizability of the classification models on independent datasets. We anticipate that the DeepMOIS-MC can be adopted to many multi-omics integrative analyses tasks.

Availability and implementation

Source codes for PyTorch implementation of DGCCA and other DeepMOIS-MC modules are available at GitHub (https://github.com/duttaprat/DeepMOIS-MC).

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Network diffusion-based approach for survival prediction and identification of biomarkers using multi-omics data of papillary renal cell carcinoma

Identification of cancer subtypes based on molecular knowledge is crucial for improving the patient diagnosis, prognosis, and treatment. In this work, we integrated copy number variations (CNVs) and transcriptomic data of Kidney Papillary Renal Cell Carcinoma (KIRP) using a network diffusion strategy to stratify cancers into clinically and biologically relevant subtypes. We constructed GeneNet, a KIRP specific gene expression network from RNA-seq data. The copy number variation data was projected onto GeneNet and propagated on the network for clustering. We identified robust subtypes that are biologically informative and significantly associated with patient survival, tumor stage and clinical subtypes of KIRP. We performed a Singular Value Decomposition (SVD) analysis of KIRP subtypes, which revealed the genes/silent players related to poor survival. A differential gene expression analysis between subtypes showed that genes related to immune, extracellular matrix organization, and genomic instability are upregulated in the poor survival group. Overall, the network-based approach revealed the molecular subtypes of KIRP and captured the relationship between gene expression and CNVs. This framework can be further expanded to integrate other omics data.

Multi-omic integration via similarity network fusion to detect molecular subtypes of ageing

Molecular subtyping of brain tissue provides insights into the heterogeneity of common neurodegenerative conditions, such as Alzheimer's disease. However, existing subtyping studies have mostly focused on single data modalities and only those individuals with severe cognitive impairment. To address these gaps, we applied similarity network fusion, a method capable of integrating multiple high-dimensional multi-omic data modalities simultaneously, to an elderly sample spanning the full spectrum of cognitive ageing trajectories. We analyzed human frontal cortex brain samples characterized by five omic modalities: bulk RNA sequencing (18 629 genes), DNA methylation (53 932 CpG sites), histone acetylation (26 384 peaks), proteomics (7737 proteins) and metabolomics (654 metabolites). Similarity network fusion followed by spectral clustering was used for subtype detection, and subtype numbers were determined by Eigen-gap and rotation cost statistics. Normalized mutual information determined the relative contribution of each modality to the fused network. Subtypes were characterized by associations with 13 age-related neuropathologies and cognitive decline. Fusion of all five data modalities (n = 111) yielded two subtypes (n S1 = 53, n S2 = 58), which were nominally associated with diffuse amyloid plaques; however, this effect was not significant after correction for multiple testing. Histone acetylation (normalized mutual information = 0.38), DNA methylation (normalized mutual information = 0.18) and RNA abundance (normalized mutual information = 0.15) contributed most strongly to this network. Secondary analysis integrating only these three modalities in a larger subsample (n = 513) indicated support for both three- and five-subtype solutions, which had significant overlap, but showed varying degrees of internal stability and external validity. One subtype showed marked cognitive decline, which remained significant even after correcting for tests across both three- and five-subtype solutions (p Bonf = 5.9 × 10-3). Comparison to single-modality subtypes demonstrated that the three-modal subtypes were able to uniquely capture cognitive variability. Comprehensive sensitivity analyses explored influences of sample size and cluster number parameters. We identified highly integrative molecular subtypes of ageing derived from multiple high dimensional, multi-omic data modalities simultaneously. Fusing RNA abundance, DNA methylation, and histone acetylation measures generated subtypes that were associated with cognitive decline. This work highlights the potential value and challenges of multi-omic integration in unsupervised subtyping of post-mortem brain.

Robust correlation estimation and UMAP assisted topological analysis of omics data for disease subtyping

Deciphering information hidden in the gene expression assays for identifying disease subtypes has significant importance in precision medicine. However, computational limitations thwart this process due to the intricacy of the biological networks and the curse of dimensionality of gene expression data. Therefore, clustering in such scenarios often becomes the first choice of exploratory data analysis to identify natural structures and intrinsic patterns in the data. However, sparse and high dimensional nature of omics data prevents conventional clustering algorithms to discover subtypes that are clinically relevant and statistically significant. Hence, non-linear dimensionality reduction techniques coupled with clustering in such scenarios often becomes imperative to improve the clustering results. In this study, we present a robust pipeline to discover disease subtypes with clinical relevance. Specifically, we focus on discovering patient sub-groups that have a residual life patterns remarkably different from other sub-groups. This is significant because by refining prognosis, subtyping can reduce uncertainty in approximating patients expected outcome. The methodology present is based on robust correlation estimation, UMAP- a non-linear dimensionality reduction method and mapper- a tool from topology. Notably, we suggest a method for improving the robustness of the correlation matrix of gene expression data for improving the clustering results. The performance of the model is evaluated by applying to five cancer datasets obtained through TCGA and comparisons are performed with some state of the art methods of NEMO, RSC-OTRI and SNF with regard to log-rank test and Restricted Life Expectancy Difference. For example in GBM dataset, the minimum separation for any two discovered subtypes is 221 days which is significantly higher than the other methodologies. We also compared the results without using the robust correlation based estimate and observed that robust correlation improves separability between survival curves significantly. From the results we infer that our methodology performs better compared to other methodologies with regard to separating survival curves of patient sub-groups despite using single omics profiles of patients compared to multiple omics profiles of SNF and NEMO. Pathway over-representation analysis is performed on the final clustering results to investigate the biological underpinnings characterizing each subtype.

Classification of the mitochondrial ribosomal protein-associated molecular subtypes and identified a serological diagnostic biomarker in hepatocellular carcinoma

Purpose

The objective of this study was to sort out innovative molecular subtypes associated with mitochondrial ribosomal proteins (MRPs) to predict clinical therapy response and determine the presence of circulating markers in hepatocellular carcinoma (HCC) patients.

Methods

Using an unsupervised clustering method, we categorized the relative molecular subtypes of MRPs in HCC patients. The prognosis, biological properties, immune checkpoint inhibitor and chemotherapy response of the patients were clarified. A signature and nomogram were developed to evaluate the prognosis. Enzyme-linked immunosorbent assay (ELISA) measured serum mitochondrial ribosomal protein L9 (MRPL9) levels in liver disease patients and normal individuals. Receiver operating characteristic (ROC) curves were conducted to calculate the diagnostic effect. The Cell Counting Kit 8 was carried out to examine cell proliferation, and flow cytometry was used to investigate the cell cycle. Transwell assay was applied to investigate the potential of cell migration and invasion. Western blot detected corresponding changes of biological markers.

Results

Participants were classified into two subtypes according to MRPs expression levels, which were characterized by different prognoses, biological features, and marked differences in response to chemotherapy and immune checkpoint inhibitors. Serum MRPL9 was significantly higher in HCC patients than in normal individuals and the benign liver disease group. ROC curve analysis showed that MRPL9 was superior to AFP and Ferritin in differentiating HCC from healthy and benign patients, or alone. Overexpressed MRPL9 could enhance aggressiveness and facilitate the G1/S progression in HCC cells.

Conclusion

We constructed novel molecular subtypes based on MRPs expression in HCC patients, which provided valuable strategies for the prediction of prognosis and clinical personalized treatment. MRPL9 might act as a reliable circulating diagnostic biomarker and therapeutic target for HCC patients.

Multiview Robust Graph-Based Clustering for Cancer Subtype Identification

Cancer subtype identification is to classify cancer into groups according to their molecular characteristics and clinical manifestations and is the basis for more personalized diagnosis and therapy. Public datasets such as The Cancer Genome Atlas (TCGA) have collected a massive number of multi-omics data. The accumulation of these datasets provides unprecedented opportunities to study the mechanism of cancers and further identify cancer subtypes at a comprehensive level. In this paper, we propose a multi-view robust graph-based clustering (MRGC) method to effectively identify cancer subtypes. Our method first learns robust latent representations from the raw omics data to alleviate the influences of the noise, where a set of similarity matrices are then adaptively learned based on these new representations. Finally, a global similarity graph is obtained by exploiting the consensus structure from the graphs. As a result, the three parts in our method can reinforce each other in a mutual iterative manner. We conduct extensive experiments on both generic machine learning datasets and cancer datasets. The experimental results confirm that our model can achieve satisfactory clustering performance compared to several state-of-the-art approaches. Moreover, we convey the practicability of MRGC by carrying out a case study on hepatocellular carcinoma.

A Bayesian model for identifying cancer subtypes from paired methylation profiles

Aberrant DNA methylation is the most common molecular lesion that is crucial for the occurrence and development of cancer, but has thus far been underappreciated as a clinical tool for cancer classification, diagnosis or as a guide for therapeutic decisions. Partly, this has been due to a lack of proven algorithms that can use methylation data to stratify patients into clinically relevant risk groups and subtypes that are of prognostic importance. Here, we proposed a novel Bayesian model to capture the methylation signatures of different subtypes from paired normal and tumor methylation array data. Application of our model to synthetic and empirical data showed high clustering accuracy, and was able to identify the possible epigenetic cause of a cancer subtype.

The path towards consensus genome classification of diffuse large B-cell lymphoma for use in clinical practice

Diffuse large B-cell lymphoma (DLBCL) is a widely heterogeneous disease in presentation, treatment response and outcome that results from a broad biological heterogeneity. Various stratification approaches have been proposed over time but failed to sufficiently capture the heterogeneous biology and behavior of the disease in a clinically relevant manner. The most recent DNA-based genomic subtyping studies are a major step forward by offering a level of refinement that could serve as a basis for exploration of personalized and targeted treatment for the years to come. To enable consistent trial designs and allow meaningful comparisons between studies, harmonization of the currently available knowledge into a single genomic classification widely applicable in daily practice is pivotal. In this review, we investigate potential avenues for harmonization of the presently available genomic subtypes of DLBCL inspired by consensus molecular classifications achieved for other malignancies. Finally, suggestions for laboratory techniques and infrastructure required for successful clinical implementation are described.