Results and lessons learned from the sbv IMPROVER metagenomics diagnostics for inflammatory bowel disease challenge.

A growing body of evidence links gut microbiota changes with inflammatory bowel disease (IBD), raising the potential benefit of exploiting metagenomics data for non-invasive IBD diagnostics. The sbv IMPROVER metagenomics diagnosis for inflammatory bowel disease challenge investigated computational metagenomics methods for discriminating IBD and nonIBD subjects. Participants in this challenge were given independent training and test metagenomics data from IBD and nonIBD subjects, which could be wither either raw read data (sub-challenge 1, SC1) or processed Taxonomy- and Function-based profiles (sub-challenge 2, SC2). A total of 81 anonymized submissions were received between September 2019 and March 2020. Most participants' predictions performed better than random predictions in classifying IBD versus nonIBD, Ulcerative Colitis (UC) versus nonIBD, and Crohn's Disease (CD) versus nonIBD. However, discrimination between UC and CD remains challenging, with the classification quality similar to the set of random predictions. We analyzed the class prediction accuracy, the metagenomics features by the teams, and computational methods used. These results will be openly shared with the scientific community to help advance IBD research and illustrate the application of a range of computational methodologies for effective metagenomic classification.

Identify, quantify and characterize cellular communication from single-cell RNA sequencing data with scSeqComm.

MOTIVATION: Recently, single-cell RNA-seq (scRNA-seq) data have been used to study cellular communication. Most bioinformatics methods infer only the intercellular signaling between groups of cells, mainly exploiting ligand-receptor expression levels. Only few methods consider the entire intercellular + intracellular signaling, mainly inferring lists/networks of signaling involved genes. RESULTS: Here, we present scSeqComm, a computational method to identify and quantify the evidence of ongoing intercellular and intracellular signaling from scRNA-seq data, and at the same time providing a functional characterization of the inferred cellular communication. The possibility to quantify the evidence of ongoing communication assists the prioritization of the results, while the combined evidence of both intercellular and intracellular signaling increase the reliability of inferred communication. The application to a scRNA-seq dataset of tumor microenvironment, the agreement with independent bioinformatics analysis, the validation using spatial transcriptomics data and the comparison with state-of-the-art intercellular scoring schemes confirmed the robustness and reliability of the proposed method. AVAILABILITY AND IMPLEMENTATION: scSeqComm R package is freely available at https://gitlab.com/sysbiobig/scseqcomm and https://sysbiobig.dei.unipd.it/software/#scSeqComm. Submitted software version and test data are available in Zenodo, at https://dx.doi.org/10.5281/zenodo.5833298. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MAST: a hybrid Multi-Agent Spatio-Temporal model of tumor microenvironment informed using a data-driven approach.

Motivation: Recently, several computational modeling approaches, such as agent-based models, have been applied to study the interaction dynamics between immune and tumor cells in human cancer. However, each tumor is characterized by a specific and unique tumor microenvironment, emphasizing the need for specialized and personalized studies of each cancer scenario. Results: We present MAST, a hybrid Multi-Agent Spatio-Temporal model which can be informed using a data-driven approach to simulate unique tumor subtypes and tumor-immune dynamics starting from high-throughput sequencing data. It captures essential components of the tumor microenvironment by coupling a discrete agent-based model with a continuous partial differential equations-based model.The application to real data of human colorectal cancer tissue investigating the spatio-temporal evolution and emergent properties of four simulated human colorectal cancer subtypes, along with their agreement with current biological knowledge of tumors and clinical outcome endpoints in a patient cohort, endorse the validity of our approach. Availability and implementation: MAST, implemented in Python language, is freely available with an open-source license through GitLab (https://gitlab.com/sysbiobig/mast), and a Docker image is provided to ease its deployment. The submitted software version and test data are available in Zenodo at https://dx.doi.org/10.5281/zenodo.7267745. Supplementary information: Supplementary data are available at Bioinformatics Advances online.

Nilotinib-induced bone marrow CD34+/lin-Ph+ cells early clearance in newly diagnosed CP-Chronic Myeloid Leukemia: Final report of the PhilosoPhi34 study.

Chronic Myeloid Leukemia is a clonal disorder characterized by the presence of the Ph-chromosome and the BCR-ABL tyrosine-kinase (TK). Target-therapy with Imatinib has greatly improved its outcome. Deeper and faster responses are reported with the second-generation TKI Nilotinib. Sustained responses may enable TKI discontinuation. However, even in a complete molecular response, some patients experience disease recurrence possibly due to persistence of quiescent leukemic CD34+/lin-Ph+ stem cells (LSCs). Degree and mechanisms of LSCs clearance during TKI treatment are not clearly established. The PhilosoPhi34 study was designed to verify the in-vivo activity and timecourse of first-line Nilotinib therapy on BM CD34+/lin-Ph+ cells clearance. Eighty-seven CP-CML patients were enrolled. BM cells were collected and tested for Ph+ residual cells, at diagnosis, 3, 6 and 12 months of treatment. FISH analysis of unstimulated CD34+/lin- cells in CCyR patients were positive in 8/65 (12.3%), 5/71 (7%), 0/69 (0%) evaluable tests, respectively. Per-Protocol analysis response rates were as follows: CCyR 95% at 12 months, MR4.5 31% and 46% at 12 and 36 months, respectively. An exploratory Gene Expression Profiling (GEP) study of CD34+/lin- cells was performed on 30 patients at diagnosis and after, on 79 patients at diagnosis vs 12 months of nilotinib treatment vs 10 healthy subjects. Data demonstrated some genes significantly different expressed: NFKBIA, many cell cycle genes, ABC transporters, JAK-STAT signaling pathway (JAK2). In addition, a correlation between different expression of some genes (JAK2, OLFM4, ICAM1, NFKBIA) among patients at diagnosis and their achievement of an early and deeper MR was observed.

Nilotinib interferes with cell cycle, ABC transporters and JAK-STAT signaling pathway in CD34+/lin- cells of patients with chronic phase chronic myeloid leukemia after 12 months of treatment.

Chronic myeloid leukemia (CML) is characterized by the constitutive tyrosine kinase activity of the oncoprotein BCR-ABL1 in myeloid progenitor cells that activates multiple signal transduction pathways leading to the leukemic phenotype. The tyrosine-kinase inhibitor (TKI) nilotinib inhibits the tyrosine kinase activity of BCR-ABL1 in CML patients. Despite the success of nilotinib treatment in patients with chronic-phase (CP) CML, a population of Philadelphia-positive (Ph+) quiescent stem cells escapes the drug activity and can lead to drug resistance. The molecular mechanism by which these quiescent cells remain insensitive is poorly understood. The aim of this study was to compare the gene expression profiling (GEP) of bone marrow (BM) CD34+/lin- cells from CP-CML patients at diagnosis and after 12 months of nilotinib treatment by microarray, in order to identify gene expression changes and the dysregulation of pathways due to nilotinib action. We selected BM CD34+/lin- cells from 78 CP-CML patients at diagnosis and after 12 months of first-line nilotinib therapy and microarray analysis was performed. GEP bioinformatic analyses identified 2,959 differently expressed probes and functional clustering determined some significantly enriched pathways between diagnosis and 12 months of nilotinib treatment. Among these pathways, we observed the under expression of 26 genes encoding proteins belonging to the cell cycle after 12 months of nilotinib treatment which led to the up-regulation of chromosome replication, cell proliferation, DNA replication, and DNA damage checkpoint at diagnosis. We demonstrated the under expression of the ATP-binding cassette (ABC) transporters ABCC4, ABCC5, and ABCD3 encoding proteins which pumped drugs out of the cells after 12 months of nilotinib. Moreover, GEP data demonstrated the deregulation of genes involved in the JAK-STAT signaling pathway. The down-regulation of JAK2, IL7, STAM, PIK3CA, PTPN11, RAF1, and SOS1 key genes after 12 months of nilotinib could demonstrate the up-regulation of cell cycle, proliferation and differentiation via MAPK and PI3K-AKT signaling pathways at diagnosis.