Application of non-negative matrix factorization in oncology: one approach for establishing precision medicine

A disulfidptosis-associated long noncoding RNA signature to predict low-grade glioma classification, prognosis, tumor microenvironment, and therapy regimens: Observational study

This study aimed to investigate the function of disulfidptosis-associated long noncoding RNAs (DAlncRNAs) in low-grade gliomas (LGG) through bioinformatics analysis and construct a signature to predict the classification, prognosis, tumor microenvironment, and selection of immunotherapy and chemotherapy in LGG. Genomic, clinical, and mutational information of 526 patients with LGG was retrieved from The Cancer Genome Atlas repository. A nonnegative matrix factorization algorithm was applied to classify patients with LGG. Univariate, LASSO regression, and multivariate Cox regression analyses were performed to determine prognostic DAlncRNAs. Following the median risk score, we defined the sample as a high-risk (HR) or low-risk group. Finally, survival, receiver operating characteristic curve, risk curve, principal component, independent prognosis, risk difference, functional enrichment, tumor microenvironment, immune cell infiltration, mutation, and drug sensitivity analyses were performed. Patients were classified into C1 and C2 subtypes associated with disulfidptosis. Eight prognostic DAlncRNAs (AC003035.2, AC010157.2, AC010273.3, AC011444.3, AC092667.1, AL450270.1, AL645608.2, and LINC01571) were identified, and a prognostic signature of LGG was developed. The DAlncRNA-based signature was found to be an independent prognostic factor in patients with LGG, thereby constructing a nomogram. In addition, in the HR group, immune function was more active and the tumor mutation burden was higher. The patients were mainly composed of subtype C2, and their prognosis was worse. Immunotherapy and chemotherapy were predicted in the HR and low-risk groups, respectively. Our study, based on DAlncRNAs, highlights 2 disulfidptosis-associated LGG subtypes with different prognostic and immune characteristics and creates a novel disulfidptosis-associated prognostic signature, which may inform the classification, prognosis, molecular pathogenesis, and therapeutic strategies for patients with LGG.

A BERT-GNN Approach for Metastatic Breast Cancer Prediction Using Histopathology Reports

Metastatic breast cancer (MBC) continues to be a leading cause of cancer-related deaths among women. This work introduces an innovative non-invasive breast cancer classification model designed to improve the identification of cancer metastases. While this study marks the initial exploration into predicting MBC, additional investigations are essential to validate the occurrence of MBC. Our approach combines the strengths of large language models (LLMs), specifically the bidirectional encoder representations from transformers (BERT) model, with the powerful capabilities of graph neural networks (GNNs) to predict MBC patients based on their histopathology reports. This paper introduces a BERT-GNN approach for metastatic breast cancer prediction (BG-MBC) that integrates graph information derived from the BERT model. In this model, nodes are constructed from patient medical records, while BERT embeddings are employed to vectorise representations of the words in histopathology reports, thereby capturing semantic information crucial for classification by employing three distinct approaches (namely univariate selection, extra trees classifier for feature importance, and Shapley values to identify the features that have the most significant impact). Identifying the most crucial 30 features out of 676 generated as embeddings during model training, our model further enhances its predictive capabilities. The BG-MBC model achieves outstanding accuracy, with a detection rate of 0.98 and an area under curve (AUC) of 0.98, in identifying MBC patients. This remarkable performance is credited to the model's utilisation of attention scores generated by the LLM from histopathology reports, effectively capturing pertinent features for classification.

Glycometabolism and lipid metabolism related genes predict the prognosis of endometrial carcinoma and their effects on tumor cells

BACKGROUND

Glycometabolism and lipid metabolism are critical in cancer metabolic reprogramming. The primary aim of this study was to develop a prognostic model incorporating glycometabolism and lipid metabolism-related genes (GLRGs) for accurate prognosis assessment in patients with endometrial carcinoma (EC).

METHODS

Data on gene expression and clinical details were obtained from publicly accessible databases. GLRGs were obtained from the Genecards database. Through nonnegative matrix factorization (NMF) clustering, molecular groupings with various GLRG expression patterns were identified. LASSO Cox regression analysis was employed to create a prognostic model. Use rich algorithms such as GSEA, GSVA, xCELL ssGSEA, EPIC,CIBERSORT, MCPcounter, ESTIMATE, TIMER, TIDE, and Oncoppredict to analyze functional pathway characteristics of the forecast signal, immune status, anti-tumor therapy, etc. The expression was assessed using Western blot and quantitative real-time PCR techniques. A total of 113 algorithm combinations were combined to screen out the most significant GLRGs in the signature for in vitro experimental verification, such as colony formation, EdU cell proliferation, wound healing, apoptosis, and Transwell assays.

RESULTS

A total of 714 GLRGs were found, and 227 of them were identified as prognostic-related genes. And ten GLRGs (AUP1, ESR1, ERLIN2, ASS1, OGDH, BCKDHB, SLC16A1, HK2, LPCAT1 and PGR-AS1) were identified to construct the prognostic model of patients with EC. Based on GLRGs, the risk model's prognosis and independent prognostic value were established. The signature of GLRGs exhibited a robust correlation with the infiltration of immune cells and the sensitivity to drugs. In cytological experiments, we selected HK2 as candidate gene to verify its value in the occurrence and development of EC. Western blot and qRT-PCR revealed that HK2 was substantially expressed in EC cells. According to in vitro experiments, HK2 knockdown can increase EC cell apoptosis while suppressing EC cell migration, invasion, and proliferation.

CONCLUSION

The GLRGs signature constructed in this study demonstrated significant prognostic value for patients with endometrial carcinoma, thereby providing valuable guidance for treatment decisions.

q-Diffusion leverages the full dimensionality of gene coexpression in single-cell transcriptomics

Unlocking the full dimensionality of single-cell RNA sequencing data (scRNAseq) is the next frontier to a richer, fuller understanding of cell biology. We introduce q-diffusion, a framework for capturing the coexpression structure of an entire library of genes, improving on state-of-the-art analysis tools. The method is demonstrated via three case studies. In the first, q-diffusion helps gain statistical significance for differential effects on patient outcomes when analyzing the CALGB/SWOG 80405 randomized phase III clinical trial, suggesting precision guidance for the treatment of metastatic colorectal cancer. Secondly, q-diffusion is benchmarked against existing scRNAseq classification methods using an in vitro PBMC dataset, in which the proposed method discriminates IFN-γ stimulation more accurately. The same case study demonstrates improvements in unsupervised cell clustering with the recent Tabula Sapiens human atlas. Finally, a local distributional segmentation approach for spatial scRNAseq, driven by q-diffusion, yields interpretable structures of human cortical tissue.

Independent prognostic value of CLDN6 in bladder cancer based on M2 macrophages related signature

M2 macrophages are associated with the prognosis of bladder cancer. CLDN6 has been linked to immune infiltration and is crucial for predicting the prognosis in multi-tumor. The effect of CLDN6 on M2 macrophages in bladder cancer remains elusive. Here, we compared a total of 40 machine learning algorithms, then selected optimal algorithm to develop M2 macrophages-related signature (MMRS) based on the identified M2 macrophages related module. MMRS predicted the prognosis better than other models and associated to immunotherapy response. CLDN6, as an important variable in MMRS, was an independent factor for poor prognosis. We found that CLDN6 was highly expressed and affected immune infiltration, immunotherapy response, and M2 macrophages polarization. Meanwhile, CLDN6 promoted the growth of bladder cancer and enhanced the carcinogenic effect by inducing polarization of M2 macrophages. In total, CLDN6 is an independent risk factor in MMRS to predict the prognosis of bladder cancer.

Advances in cancer DNA methylation analysis with methPLIER: use of non-negative matrix factorization and knowledge-based constraints to enhance biological interpretability

DNA methylation is an epigenetic modification that results in dynamic changes during ontogenesis and cell differentiation. DNA methylation patterns regulate gene expression and have been widely researched. While tools for DNA methylation analysis have been developed, most of them have focused on intergroup comparative analysis within a dataset; therefore, it is difficult to conduct cross-dataset studies, such as rare disease studies or cross-institutional studies. This study describes a novel method for DNA methylation analysis, namely, methPLIER, which enables interdataset comparative analyses. methPLIER combines Pathway Level Information Extractor (PLIER), which is a non-negative matrix factorization (NMF) method, with regularization by a knowledge matrix and transfer learning. methPLIER can be used to perform intersample and interdataset comparative analysis based on latent feature matrices, which are obtained via matrix factorization of large-scale data, and factor-loading matrices, which are obtained through matrix factorization of the data to be analyzed. We used methPLIER to analyze a lung cancer dataset and confirmed that the data decomposition reflected sample characteristics for recurrence-free survival. Moreover, methPLIER can analyze data obtained via different preprocessing methods, thereby reducing distributional bias among datasets due to preprocessing. Furthermore, methPLIER can be employed for comparative analyses of methylation data obtained from different platforms, thereby reducing bias in data distribution due to platform differences. methPLIER is expected to facilitate cross-sectional DNA methylation data analysis and enhance DNA methylation data resources.

Mass Spectrometry-Based Proteogenomics: New Therapeutic Opportunities for Precision Medicine

Proteogenomics refers to the integration of comprehensive genomic, transcriptomic, and proteomic measurements from the same samples with the goal of fully understanding the regulatory processes converting genotypes to phenotypes, often with an emphasis on gaining a deeper understanding of disease processes. Although specific genetic mutations have long been known to drive the development of multiple cancers, gene mutations alone do not always predict prognosis or response to targeted therapy. The benefit of proteogenomics research is that information obtained from proteins and their corresponding pathways provides insight into therapeutic targets that can complement genomic information by providing an additional dimension regarding the underlying mechanisms and pathophysiology of tumors. This review describes the novel insights into tumor biology and drug resistance derived from proteogenomic analysis while highlighting the clinical potential of proteogenomic observations and advances in technique and analysis tools.

Longitudinal clustering of Life's Essential 8 health metrics: application of a novel unsupervised learning method in the CARDIA study

OBJECTIVE

Changes in cardiovascular health (CVH) during the life course are associated with future cardiovascular disease (CVD). Longitudinal clustering analysis using subgraph augmented non-negative matrix factorization (SANMF) could create phenotypic risk profiles of clustered CVH metrics.

MATERIALS AND METHODS

Life's Essential 8 (LE8) variables, demographics, and CVD events were queried over 15 ears in 5060 CARDIA participants with 18 years of subsequent follow-up. LE8 subgraphs were mined and a SANMF algorithm was applied to cluster frequently occurring subgraphs. K-fold cross-validation and diagnostics were performed to determine cluster assignment. Cox proportional hazard models were fit for future CV event risk and logistic regression was performed for cluster phenotyping.

RESULTS

The cohort (54.6% female, 48.7% White) produced 3 clusters of CVH metrics: Healthy & Late Obesity (HLO) (29.0%), Healthy & Intermediate Sleep (HIS) (43.2%), and Unhealthy (27.8%). HLO had 5 ideal LE8 metrics between ages 18 and 39 years, until BMI increased at 40. HIS had 7 ideal LE8 metrics, except sleep. Unhealthy had poor levels of sleep, smoking, and diet but ideal glucose. Race and employment were significantly different by cluster (P < .001) but not sex (P = .734). For 301 incident CV events, multivariable hazard ratios (HRs) for HIS and Unhealthy were 0.73 (0.53-1.00, P = .052) and 2.00 (1.50-2.68, P < .001), respectively versus HLO. A 15-year event survival was 97.0% (HIS), 96.3% (HLO), and 90.4% (Unhealthy, P < .001).

DISCUSSION AND CONCLUSION

SANMF of LE8 metrics identified 3 unique clusters of CVH behavior patterns. Clustering of longitudinal LE8 variables via SANMF is a robust tool for phenotypic risk assessment for future adverse cardiovascular events.

Identification of S1PR4 as an immune modulator for favorable prognosis in HNSCC through machine learning

G protein-coupled receptors (GPCRs) are the largest family of membrane proteins and play a critical role as pharmacological targets. An improved understanding of GPCRs' involvement in tumor microenvironment may provide new perspectives for cancer therapy. This study used machine learning to classify head and neck squamous cell carcinoma (HNSCC) patients into two GPCR-based subtypes. Notably, these subtypes showed significant differences in prognosis, gene expression, and immune microenvironment, particularly CD8+ T cell infiltration. S1PR4 emerged as a key regulator distinguishing the subtypes, positively correlated with CD8+ T cell proportion and cytotoxicity in HNSCC. It was predominantly expressed in CX3CR1+CD8+ T cells among T cells. Upregulation of S1PR4 enhanced T cell function during CAR-T cell therapy, suggesting its potential in cancer immunotherapy. These findings highlight S1PR4 as an immune modulator for favorable prognosis in HNSCC, and offer a potential GPCR-targeted therapeutic option for HNSCC treatment.

Exploration of prognosis and immunometabolism landscapes in ER+ breast cancer based on a novel lipid metabolism-related signature

Introduction

Lipid metabolic reprogramming is gaining attention as a hallmark of cancers. Recent mounting evidence indicates that the malignant behavior of breast cancer (BC) is closely related to lipid metabolism. Here, we focus on the estrogen receptor-positive (ER+) subtype, the most common subgroup of BC, to explore immunometabolism landscapes and prognostic significance according to lipid metabolism-related genes (LMRGs).

Methods

Samples from The Cancer Genome Atlas (TCGA) database were used as training cohort, and samples from the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC), Gene Expression Omnibus (GEO) datasets and our cohort were applied for external validation. The survival-related LMRG molecular pattern and signature were constructed by unsupervised consensus clustering and least absolute shrinkage and selection operator (LASSO) analysis. A lipid metabolism-related clinicopathologic nomogram was established. Gene enrichment and pathway analysis were performed to explore the underlying mechanism. Immune landscapes, immunotherapy and chemotherapy response were further explored. Moreover, the relationship between gene expression and clinicopathological features was assessed by immunohistochemistry.

Results

Two LMRG molecular patterns were identified and associated with distinct prognoses and immune cell infiltration. Next, a prognostic signature based on nine survival-related LMRGs was established and validated. The signature was confirmed to be an independent prognostic factor and an optimal nomogram incorporating age and T stage (AUC of 5-year overall survival: 0.778). Pathway enrichment analysis revealed differences in immune activities, lipid biosynthesis and drug metabolism by comparing groups with low- and high-risk scores. Further exploration verified different immune microenvironment profiles, immune checkpoint expression, and sensitivity to immunotherapy and chemotherapy between the two groups. Finally, arachidonate 15-lipoxygenase (ALOX15) was selected as the most prominent differentially expressed gene between the two groups. Its expression was positively related to larger tumor size, more advanced tumor stage and vascular invasion in our cohort (n = 149).

Discussion

This is the first lipid metabolism-based signature with value for prognosis prediction and immunotherapy or chemotherapy guidance for ER+ BC.

Identification and Characterization of Metabolic Subtypes of Endometrial Cancer Using a Systems-Level Approach

Endometrial cancer (EC) is the most common gynecological cancer worldwide. Understanding metabolic adaptation and its heterogeneity in tumor tissues may provide new insights and help in cancer diagnosis, prognosis, and treatment. In this study, we investigated metabolic alterations of EC to understand the variations in metabolism within tumor samples. Integration of transcriptomics data of EC (RNA-Seq) and the human genome-scale metabolic network was performed to identify the metabolic subtypes of EC and uncover the underlying dysregulated metabolic pathways and reporter metabolites in each subtype. The relationship between metabolic subtypes and clinical variables was explored. Further, we correlated the metabolic changes occurring at the transcriptome level with the genomic alterations. Based on metabolic profile, EC patients were stratified into two subtypes (metabolic subtype-1 and subtype-2) that significantly correlated to patient survival, tumor stages, mutation, and copy number variations. We observed the co-activation of the pentose phosphate pathway, one-carbon metabolism, and genes involved in controlling estrogen levels in metabolic subtype-2, which is linked to poor survival. PNMT and ERBB2 are also upregulated in metabolic subtype-2 samples and present on the same chromosome locus 17q12, which is amplified. PTEN and TP53 mutations show mutually exclusive behavior between subtypes and display a difference in survival. This work identifies metabolic subtypes with distinct characteristics at the transcriptome and genome levels, highlighting the metabolic heterogeneity within EC.

Deconvolution algorithms for inference of the cell-type composition of the spatial transcriptome

The spatial transcriptome has enabled researchers to resolve transcriptome expression profiles while preserving information about cell location to better understand the complex biological processes that occur in organisms. Due to technical limitations, the current high-throughput spatial transcriptome sequencing methods (known as next-generation sequencing with spatial barcoding methods or spot-based methods) cannot achieve single-cell resolution. A single measurement site, called a spot, in these technologies frequently contains multiple cells of various types. Computational tools for determining the cellular composition of a spot have emerged as a way to break through these limitations. These tools are known as deconvolution tools. Recently, a couple of deconvolution tools based on different strategies have been developed and have shown promise in different aspects. The resulting single-cell resolution expression profiles and/or single-cell composition of spots will significantly affect downstream data mining; thus, it is crucial to choose a suitable deconvolution tool. In this review, we present a list of currently available tools for spatial transcriptome deconvolution, categorize them based on the strategies they employ, and explain their advantages and limitations in detail in order to guide the selection of these tools in future studies.

Introducing AI to the molecular tumor board: one direction toward the establishment of precision medicine using large-scale cancer clinical and biological information

Since U.S. President Barack Obama announced the Precision Medicine Initiative in his New Year's State of the Union address in 2015, the establishment of a precision medicine system has been emphasized worldwide, particularly in the field of oncology. With the advent of next-generation sequencers specifically, genome analysis technology has made remarkable progress, and there are active efforts to apply genome information to diagnosis and treatment. Generally, in the process of feeding back the results of next-generation sequencing analysis to patients, a molecular tumor board (MTB), consisting of experts in clinical oncology, genetic medicine, etc., is established to discuss the results. On the other hand, an MTB currently involves a large amount of work, with humans searching through vast databases and literature, selecting the best drug candidates, and manually confirming the status of available clinical trials. In addition, as personalized medicine advances, the burden on MTB members is expected to increase in the future. Under these circumstances, introducing cutting-edge artificial intelligence (AI) technology and information and communication technology to MTBs while reducing the burden on MTB members and building a platform that enables more accurate and personalized medical care would be of great benefit to patients. In this review, we introduced the latest status of elemental technologies that have potential for AI utilization in MTB, and discussed issues that may arise in the future as we progress with AI implementation.

Multi-omics analysis: Paving the path toward achieving precision medicine in cancer treatment and immuno-oncology

The acceleration of large-scale sequencing and the progress in high-throughput computational analyses, defined as omics, was a hallmark for the comprehension of the biological processes in human health and diseases. In cancerology, the omics approach, initiated by genomics and transcriptomics studies, has revealed an incredible complexity with unsuspected molecular diversity within a same tumor type as well as spatial and temporal heterogeneity of tumors. The integration of multiple biological layers of omics studies brought oncology to a new paradigm, from tumor site classification to pan-cancer molecular classification, offering new therapeutic opportunities for precision medicine. In this review, we will provide a comprehensive overview of the latest innovations for multi-omics integration in oncology and summarize the largest multi-omics dataset available for adult and pediatric cancers. We will present multi-omics techniques for characterizing cancer biology and show how multi-omics data can be combined with clinical data for the identification of prognostic and treatment-specific biomarkers, opening the way to personalized therapy. To conclude, we will detail the newest strategies for dissecting the tumor immune environment and host–tumor interaction. We will explore the advances in immunomics and microbiomics for biomarker identification to guide therapeutic decision in immuno-oncology.