Metabolic plasticity, the ability of cells to dynamically adapt their metabolic pathways in response to changing environments, is a hallmark of rheumatoid arthritis (RA) pathogenesis and plays a critical role in immune dysfunction. However, scalable methods to quantify inter-pathway crosstalk in RA remain lacking. To address this gap, we present MMP3C v2, an updated network-based framework that integrates gene expression with protein-protein interaction network topology to compute directed pairwise metabolic plasticity (PMP) scores. We applied MMP3C v2 to ~3400 bulk transcriptomes (RA, osteoarthritis, systemic lupus erythematosus, and healthy controls) and ~228 000 single-cell transcriptomics from blood and synovium to profile RA-associated PMP alterations and develop diagnostic classifiers. We found that a single PMP-derived signature demonstrated strong predictive capability for diagnosis. Then, we developed a feature selection pipeline and combined it with 110 machine learning model combinations, by which we established the optimal ensemble classifier (stepwise forward selection + ridge regression), achieving robust and generalized performance (mean area under the curve (AUC) = 0.935; mean F1 score = 0.915) across 12 independent validation cohorts, outperforming seven previously published models. Single-cell analysis revealed cell-type-specific PMP remodeling: a Warburg-like shift in synovial macrophages (↑glycolysis, ↑pentose phosphate pathway, ↓oxidative phosphorylation). Cell-cell communication analysis highlighted FN1-centered signaling linked to glucose metabolic remodeling in myofibroblasts. Collectively, MMP3C v2 establishes metabolic pathway crosstalk as a core diagnostic feature of RA, enabling interpretable and cross-platform diagnostic modeling and the identification of cell-type-specific PMP patterns. The open-source R package mmp3c supports reproducible analysis and broad application.

Interactions between Tumor microenvironment (TME) cells shape the unique growth environment, sustaining tumor growth and causing the immune escape of tumor cells. Nonetheless, no studies have reported a systematic analysis of cellular interactions in the identification of cancer-related TME cells. Here, we proposed a novel network-based computational method, named as iATMEcell, to identify the abnormal TME cells associated with the biological outcome of interest based on a cell-cell crosstalk network. In the method, iATMEcell first manually collected TME cell types from multiple published studies and obtained their corresponding gene signatures. Then, a weighted cell-cell crosstalk network was constructed in the context of a specific cancer bulk tissue transcriptome data, where the weight between cells reflects both their biological function similarity and the transcriptional dysregulated activities of gene signatures shared by them. Finally, it used a network propagation algorithm to identify significantly dysregulated TME cells. Using the cancer genome atlas (TCGA) Bladder Urothelial Carcinoma training set and two independent validation sets, we illustrated that iATMEcell could identify significant abnormal cells associated with patient survival and immunotherapy response. iATMEcell was further applied to a pan-cancer analysis, which revealed that four common abnormal immune cells play important roles in the patient prognosis across multiple cancer types. Collectively, we demonstrated that iATMEcell could identify potentially abnormal TME cells based on a cell-cell crosstalk network, which provided a new insight into understanding the effect of TME cells in cancer. iATMEcell is developed as an R package, which is freely available on GitHub (https://github.com/hanjunwei-lab/iATMEcell).

The link between tumor genetic variations and immunotherapy benefits has been widely recognized. Recent studies suggested that the key biological pathways activated by accumulated genetic mutations may act as an effective biomarker for predicting the efficacy of immune checkpoint inhibitor (ICI) therapy. Here, we developed a novel individual Pathway Mutation Perturbation (iPMP) method that measures the pathway mutation perturbation level by combining evidence of the cumulative effect of mutated genes with the position of mutated genes in the pathways. In iPMP, somatic mutations on a single sample were first mapped to genes in a single pathway to infer the pathway mutation perturbation score (PMPscore), and then, an integrated PMPscore profile was produced, which can be used in place of the original mutation dataset to identify associations with clinical outcomes. To illustrate the effect of iPMP, we applied it to a melanoma cohort treated with ICIs and identified seven significant perturbation pathways, which jointly constructed a pathway-based signature. With the signature, patients were classified into two subgroups with significant distinctive overall survival and objective response rate to immunotherapy. Moreover, the pathway-based signature was consistently validated in two independent melanoma cohorts. We further applied iPMP to two non-small cell lung cancer cohorts and also obtained good performance. Altogether, the iPMP method could be used to identify the significant mutation perturbation pathways for constructing the pathway-based biomarker to predict the clinical outcomes of immunotherapy. The iPMP method has been implemented as a freely available R-based package (https://CRAN.R-project.org/package=PMAPscore).

Biological pathways reflect the key cellular mechanisms that dictate disease states, drug response and altered cellular function. The local areas of pathways are defined as subpathways (SPs), whose dysfunction has been reported to be associated with the occurrence and development of cancer. With the development of high-throughput sequencing technology, identifying dysfunctional SPs by using multi-omics data has become possible. Moreover, the SPs are not isolated in the biological system but interact with each other. Here, we propose a network-based calculated method, CNA2Subpathway, to identify dysfunctional SPs is driven by somatic copy number alterations (CNAs) in cancer through integrating pathway topology information, multi-omics data and SP crosstalk. This provides a novel way of SP analysis by using the SP interactions in the system biological level. Using data sets from breast cancer and head and neck cancer, we validate the effectiveness of CNA2Subpathway in identifying cancer-relevant SPs driven by the somatic CNAs, which are also shown to be associated with cancer immune and prognosis of patients. We further compare our results with five pathway or SP analysis methods based on CNA and gene expression data without considering SP crosstalk. With these analyses, we show that CNA2Subpathway could help to uncover dysfunctional SPs underlying cancer via the use of SP crosstalk. CNA2Subpathway is developed as an R-based tool, which is freely available on GitHub (https://github.com/hanjunwei-lab/CNA2Subpathway).