MOTIVATION:The identification of minimal genetic interventions that modulate metabolic processes constitutes one of the most relevant applications of genome-scale metabolic models (GEMs). The concept of Minimal Cut Sets (MCSs) and its extension at the gene level, genetic Minimal Cut Sets (gMCSs), have attracted increasing interest in the field of Systems Biology to address this task. Different computational tools have been developed to calculate MCSs and gMCSs using both commercial and open-source software. RESULTS:Here, we present gMCSpy, an efficient Python package to calculate gMCSs in GEMs using both commercial and non-commercial optimization solvers. We show that gMCSpy substantially overperforms our previous computational tool GMCS, which exclusively relied on commercial software. Moreover, we compared gMCSpy with recently published competing algorithms in the literature, finding significant improvements in both accuracy and computation time. All these advances make gMCSpy an attractive tool for researchers in the field of Systems Biology for different applications in health and biotechnology. AVAILABILITY AND IMPLEMENTATION:The Python package gMCSpy and the data underlying this manuscript can be accessed at: https://github.com/PlanesLab/gMCSpy.

Motivation Alternative splicing plays a pivotal role in various biological processes. In the context of cancer, aberrant splicing patterns can lead to disease progression and treatment resistance. Understanding the regulatory mechanisms underlying alternative splicing is crucial for elucidating disease mechanisms and identifying potential therapeutic targets. Results We present DeepRBP, a deep learning (DL) based framework to identify potential RNA-binding proteins (RBP)-Gene regulation pairs for further in-vitro validation. DeepRBP is composed of a DL model that predicts transcript abundance given RBP and gene expression data coupled with an explainability module that computes informative RBP-Gene scores. We show that the proposed framework is able to identify known RBP-Gene regulations, demonstrating its applicability to identify new ones. Availability and Implementation DeepRBP is implemented in PyTorch, and all the code and material used in this work is available at <https://github.com/ML4BM-Lab/DeepRBP>. Contact iochoal{at}unav.es Supplementary information Supplementary data are available at Bioinformatics online. ### Competing Interest Statement The authors have declared no competing interest.

Cancer metabolism is a marvellously complex topic, in part, due to the reprogramming of its pathways to self-sustain the malignant phenotype in the disease, to the detriment of its healthy counterpart. Understanding these adjustments can provide novel targeted therapies that could disrupt and impair proliferation of cancerous cells. For this very purpose, genome-scale metabolic models (GEMs) have been developed, with Human1 being the most recent reconstruction of the human metabolism. Based on GEMs, we introduced the genetic Minimal Cut Set (gMCS) approach, an uncontextualized methodology that exploits the concepts of synthetic lethality to predict metabolic vulnerabilities in cancer. gMCSs define a set of genes whose knockout would render the cell unviable by disrupting an essential metabolic task in GEMs, thus, making cellular proliferation impossible. Here, we summarize the gMCS approach and review the current state of the methodology by performing a systematic meta-analysis based on two datasets of gene essentiality in cancer. First, we assess several thresholds and distinct methodologies for discerning highly and lowly expressed genes. Then, we address the premise that gMCSs of distinct length should have the same predictive power. Finally, we question the importance of a gene partaking in multiple gMCSs and analyze the importance of all the essential metabolic tasks defined in Human1. Our meta-analysis resulted in parameter evaluation to increase the predictive power for the gMCS approach, as well as a significant reduction of computation times by only selecting the crucial gMCS lengths, proposing the pertinency of particular parameters for the peak processing of gMCS.

ABSTRACTThe development of computational tools for the systematic prediction of metabolic vulnerabilities of cancer cells constitutes a central question in systems biology. Here, we presentgMCStool, a freely accessible and online tool that allows us to carry out this task in a simple, efficient and intuitive environment.gMCStoolexploits the concept of genetic Minimal Cut Sets (gMCSs), a theoretical approach to synthetic lethality based on genome-scale metabolic networks, including a unique database of thousands of synthetic lethals computed from Human1, the most recent metabolic reconstruction of human cells. Based on RNA-seq data,gMCStoolextends and improves our previously developed algorithms to predict, visualize and analyze metabolic essential genes in cancer, demonstrating a superior performance than competing algorithms in both accuracy and computational performance. A detailed illustration ofgMCStoolis presented for multiple myeloma (MM), an incurable hematological malignancy. gMCStool could identify a synthetic lethal that explains the dependency on CTP Synthase 1 (CTPS1) in a sub-group of MM patients. We providein vitroexperimental evidence that supports this hypothesis, which opens a new research area to treat MM.

Synthetic lethality (SL) is a promising concept in cancer research. A wide array of computational tools has been developed to predict and exploit synthetic lethality for the identification of tumour-specific vulnerabilities. Previously, we introduced the concept of genetic Minimal Cut Sets (gMCSs), a theoretical approach to SL developed for genome-scale metabolic networks. The major challenge in our gMCS framework is to go beyond metabolic networks and extend existing algorithms to more complex protein-protein interactions. In this article, we take a step further and incorporate linear regulatory pathways into our gMCS approach. Extensive algorithmic modifications to compute gMCSs in integrated metabolic and regulatory models are presented in detail. Our extended approach is applied to calculate gMCSs in integrated models of human cells. In particular, we integrate the most recent genome-scale metabolic network, Human1, with 3 different regulatory network databases: Omnipath, Dorothea and TRRUST. Based on the computed gMCSs and transcriptomic data, we discovered new essential genes and their associated synthetic lethal for different cancer cell lines. The performance of the different integrated models is assessed with available large-scale in-vitro gene silencing data. Finally, we discuss the most relevant gene essentiality predictions based on published literature in cancer research.

Alternative splicing (AS) plays a key role in cancer: all its hallmarks have been associated with different mechanisms of abnormal AS. The improvement of the human transcriptome annotation and the availability of fast and accurate software to estimate isoform concentrations has boosted the analysis of transcriptome profiling from RNA-seq. The statistical analysis of AS is a challenging problem not yet fully solved. We have included in EventPointer (EP), a Bioconductor package, a novel statistical method that can use the bootstrap of the pseudoaligners. We compared it with other state-of-the-art algorithms to analyze AS. Its performance is outstanding for shallow sequencing conditions. The statistical framework is very flexible since it is based on design and contrast matrices. EP now includes a convenient tool to find the primers to validate the discoveries using PCR. We also added a statistical module to study alteration in protein domain related to AS. Applying it to 9514 patients from TCGA and TARGET in 19 different tumor types resulted in two conclusions: i) aberrant alternative splicing alters the relative presence of Protein domains and, ii) the number of enriched domains is strongly correlated with the age of the patients.

Eradicating leukemia requires a deep understanding of the interaction between leukemic cells and their protective microenvironment. The CXCL12/CXCR4 axis has been postulated as a critical pathway dictating leukemia stem cell (LSC) chemoresistance in AML due to its role in controlling cellular egress from the marrow. Nevertheless, the cellular source of CXCL12 in the acute myeloid leukemia (AML) microenvironment and the mechanism by which CXCL12 exerts its protective role in vivo remain unresolved. Here, we show that CXCL12 produced by Prx1+ mesenchymal cells but not by mature osteolineage cells provide the necessary cues for the maintenance of LSCs in the marrow of an MLL::AF9-induced AML model. Prx1+ cells promote survival of LSCs by modulating energy metabolism and the REDOX balance in LSCs. Deletion of Cxcl12 leads to the accumulation of reactive oxygen species and DNA damage in LSCs, impairing their ability to perpetuate leukemia in transplantation experiments, a defect that can be attenuated by antioxidant therapy. Importantly, our data suggest that this phenomenon appears to be conserved in human patients. Hence, we have identified Prx1+ mesenchymal cells as an integral part of the complex niche-AML metabolic intertwining, pointing towards CXCL12/CXCR4 as a target to eradicate parenchymal LSCs in AML.

The presence of anti-myelin lipid-specific oligoclonal IgM bands (LS-OCMBs) has been defined as an accurate predictor of an aggressive evolution of multiple sclerosis. However, the detection of this biomarker is performed in cerebrospinal fluid, a quite invasive liquid biopsy. In the present study we aimed at studying the expression profile of miRNA, snoRNA, circRNA and linearRNA in peripheral blood mononuclear cells (PBMCs) from patients with lipid-specific oligoclonal IgM band characterization. We included a total of 89 MS patients, 47 with negative LS-OCMB status and 42 with positive status. Microarray (miRNA and snoRNA) and RNA-seq (circular and linear RNAs) were used to perform the profiling study in the discovery cohort and candidates were validated by RT-qPCR in the whole cohort. The biomarker potential of the candidates was evaluated by ROC curve analysis. RNA-seq and RT-qPCR validation revealed that two circular (hsa_circ_0000478 and hsa_circ_0116639) and two linear RNAs (IRF5 and MTRNR2L8) are downregulated in PBMCs from patients with positive LS-OCMBs. Finally, those RNAs show a performance of a 70% accuracy in some of the combinations. The expression of hsa_circ_0000478, hsa_circ_0116639, IRF5 and MTRNR2L8 might serve as minimally invasive biomarkers of highly active disease.