Highlights d Phosphoproteomics analysis of SARS-CoV-2-infected cells uncovers signaling rewiring d Infection promotes host p38 MAPK cascade activity and shutdown of mitotic kinases d Infection stimulates CK2-containing filopodial protrusions with budding virus d Kinase activity analysis identifies potent antiviral drugs and compounds

Bone marrow fibrosis (BMF) develops in various hematological and non-hematological conditions and is a central pathological feature of myelofibrosis. Effective cell-targeted therapeutics are needed, but the cellular origin of BMF remains elusive. Here, we show using genetic fate tracing in two murine models of BMF that Gli1 mesenchymal stromal cells (MSCs) are recruited from the endosteal and perivascular niche to become fibrosis-driving myofibroblasts in the bone marrow. Genetic ablation of Gli1 cells abolished BMF and rescued bone marrow failure. Pharmacological targeting of Gli proteins with GANT61 inhibited Gli1 cell expansion and myofibroblast differentiation and attenuated fibrosis severity. The same pathway is also active in human BMF, and Gli1 expression in BMF significantly correlates with the severity of the disease. In addition, GANT61 treatment reduced the myofibroblastic phenotype of human MSCs isolated from patients with BMF, suggesting that targeting of Gli proteins could be a relevant therapeutic strategy.

SummaryThe bioenergetics and molecular determinants of the metabolic response to mitochondrial dysfunction are incompletely understood, in part due to a lack of appropriate isogenic cellular models of primary mitochondrial defects. Here, we capitalize on a recently developed cell model with defined levels of m.8993T>G mutation heteroplasmy, mTUNE, to investigate the metabolic underpinnings of mitochondrial dysfunction. We found that impaired utilization of reduced nicotinamide adenine dinucleotide (NADH) by the mitochondrial respiratory chain leads to cytosolic reductive carboxylation of glutamine as a new mechanism for cytosol-confined NADH recycling supported by malate dehydrogenase 1 (MDH1). We also observed that increased glycolysis in cells with mitochondrial dysfunction is associated with increased cell migration in an MDH1-dependent fashion. Our results describe a novel link between glycolysis and mitochondrial dysfunction mediated by reductive carboxylation of glutamine.

Supplementary data are available at Bioinformatics online.

The growing availability of single-cell data, especially transcriptomics, has sparked an increased interest in the inference of cell-cell communication. Many computational tools were developed for this purpose. Each of them consists of a resource of intercellular interactions prior knowledge and a method to predict potential cell-cell communication events. Yet the impact of the choice of resource and method on the resulting predictions is largely unknown. To shed light on this, we systematically compare 16 cell-cell communication inference resources and 7 methods, plus the consensus between the methods’ predictions. Among the resources, we find few unique interactions, a varying degree of overlap, and an uneven coverage of specific pathways and tissue-enriched proteins. We then examine all possible combinations of methods and resources and show that both strongly influence the predicted intercellular interactions. Finally, we assess the agreement of cell-cell communication methods with spatial colocalisation, cytokine activities, and receptor protein abundance and find that predictions are generally coherent with those data modalities. To facilitate the use of the methods and resources described in this work, we provide LIANA, a LIgand-receptor ANalysis frAmework as an open-source interface to all the resources and methods.

Summary Many methods allow us to extract biological activities from omics data using information from prior knowledge resources, reducing the dimensionality for increased statistical power and better interpretability. Here, we present decoupleR, a Bioconductor and Python package containing computational methods to extract these activities within a unified framework. decoupleR allows us to flexibly run any method with a given resource, including methods that leverage mode of regulation and weights of interactions, which are not present in other frameworks. Moreover, it leverages OmniPath, a meta-resource comprising over 100 databases of prior knowledge. Using decoupleR, we evaluated the performance of methods on transcriptomic and phospho-proteomic perturbation experiments. Our findings suggest that simple linear models and the consensus score across top methods performs better than other methods at predicting perturbed regulators. Availability and Implementation decoupleR’s open source code is available in Bioconductor (https://www.bioconductor.org/packages/release/bioc/html/decoupleR.html) for R and in GitHub (https://github.com/saezlab/decoupler-py) for Python. The code to reproduce the results is in GitHub (https://github.com/saezlab/decoupleR_manuscript) and the data in Zenodo (https://zenodo.org/record/5645208).

Multi-omics datasets can provide molecular insights beyond the sum of individual omics. Various tools have been recently developed to integrate such datasets, but there are limited strategies to systematically extract mechanistic hypotheses from them. Here, we present COSMOS (Causal Oriented Search of Multi-Omics Space), a method that integrates phosphoproteomics, transcriptomics, and metabolomics datasets. COSMOS combines extensive prior knowledge of signaling, metabolic, and gene regulatory networks with computational methods to estimate activities of transcription factors and kinases as well as network-level causal reasoning. COSMOS provides mechanistic hypotheses for experimental observations across multi-omics datasets. We applied COSMOS to a dataset comprising transcriptomics, phosphoproteomics, and metabolomics data from healthy and cancerous tissue from eleven clear cell renal cell carcinoma (ccRCC) patients. COSMOS was able to capture relevant crosstalks within and between multiple omics layers, such as known ccRCC drug targets. We expect that our freely available method will be broadly useful to extract mechanistic insights from multi-omics studies.

No abstract