scPerturb: harmonized single-cell perturbation data.

Analysis across a growing number of single-cell perturbation datasets is hampered by poor data interoperability. To facilitate development and benchmarking of computational methods, we collect a set of 44 publicly available single-cell perturbation-response datasets with molecular readouts, including transcriptomics, proteomics and epigenomics. We apply uniform quality control pipelines and harmonize feature annotations. The resulting information resource, scPerturb, enables development and testing of computational methods, and facilitates comparison and integration across datasets. We describe energy statistics (E-statistics) for quantification of perturbation effects and significance testing, and demonstrate E-distance as a general distance measure between sets of single-cell expression profiles. We illustrate the application of E-statistics for quantifying similarity and efficacy of perturbations. The perturbation-response datasets and E-statistics computation software are publicly available at scperturb.org. This work provides an information resource for researchers working with single-cell perturbation data and recommendations for experimental design, including optimal cell counts and read depth.

Disulfiram use is associated with lower risk of COVID-19: A retrospective cohort study.

Effective, low-cost therapeutics are needed to prevent and treat COVID-19. Severe COVID-19 disease is linked to excessive inflammation. Disulfiram is an approved oral drug used to treat alcohol use disorder that is a potent anti-inflammatory agent and an inhibitor of the viral proteases. We investigated the potential effects of disulfiram on SARS-CoV-2 infection and disease severity in an observational study using a large database of clinical records from the national US Veterans Affairs healthcare system. A multivariable Cox regression adjusted for demographic information and diagnosis of alcohol use disorder revealed a reduced risk of SARS-CoV-2 infection with disulfiram use at a hazard ratio of 0.66 (34% lower risk, 95% confidence interval 24-43%). There were no COVID-19 related deaths among the 188 SARS-CoV-2 positive patients treated with disulfiram, in contrast to 5-6 statistically expected deaths based on the untreated population (P = 0.03). Our epidemiological results suggest that disulfiram may contribute to the reduced incidence and severity of COVID-19. These results support carefully planned clinical trials to assess the potential therapeutic effects of disulfiram in COVID-19.

Molecular response to PARP1 inhibition in ovarian cancer cells as determined by mass spectrometry based proteomics.

BACKGROUND: Poly (ADP)-ribose polymerase (PARP) inhibitors have entered routine clinical practice for the treatment of high-grade serous ovarian cancer (HGSOC), yet the molecular mechanisms underlying treatment response to PARP1 inhibition (PARP1i) are not fully understood. METHODS: Here, we used unbiased mass spectrometry based proteomics with data-driven protein network analysis to systematically characterize how HGSOC cells respond to PARP1i treatment. RESULTS: We found that PARP1i leads to pronounced proteomic changes in a diverse set of cellular processes in HGSOC cancer cells, consistent with transcript changes in an independent perturbation dataset. We interpret decreases in the levels of the pro-proliferative transcription factors SP1 and ß-catenin and in growth factor signaling as reflecting the anti-proliferative effect of PARP1i; and the strong activation of pro-survival processes NF-κB signaling and lipid metabolism as PARPi-induced adaptive resistance mechanisms. Based on these observations, we nominate several protein targets for therapeutic inhibition in combination with PARP1i. When tested experimentally, the combination of PARPi with an inhibitor of fatty acid synthase (TVB-2640) has a 3-fold synergistic effect and is therefore of particular pre-clinical interest. CONCLUSION: Our study improves the current understanding of PARP1 function, highlights the potential that the anti-tumor efficacy of PARP1i may not only rely on DNA damage repair mechanisms and informs on the rational design of PARP1i combination therapies in ovarian cancer.

CellBox: Interpretable Machine Learning for Perturbation Biology with Application to the Design of Cancer Combination Therapy.

Systematic perturbation of cells followed by comprehensive measurements of molecular and phenotypic responses provides informative data resources for constructing computational models of cell biology. Models that generalize well beyond training data can be used to identify combinatorial perturbations of potential therapeutic interest. Major challenges for machine learning on large biological datasets are to find global optima in a complex multidimensional space and mechanistically interpret the solutions. To address these challenges, we introduce a hybrid approach that combines explicit mathematical models of cell dynamics with a machine-learning framework, implemented in TensorFlow. We tested the modeling framework on a perturbation-response dataset of a melanoma cell line after drug treatments. The models can be efficiently trained to describe cellular behavior accurately. Even though completely data driven and independent of prior knowledge, the resulting de novo network models recapitulate some known interactions. The approach is readily applicable to various kinetic models of cell biology. A record of this paper's Transparent Peer Review process is included in the Supplemental Information.

Analyses of non-coding somatic drivers in 2,658 cancer whole genomes.

The discovery of drivers of cancer has traditionally focused on protein-coding genes1-4. Here we present analyses of driver point mutations and structural variants in non-coding regions across 2,658 genomes from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium5 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). For point mutations, we developed a statistically rigorous strategy for combining significance levels from multiple methods of driver discovery that overcomes the limitations of individual methods. For structural variants, we present two methods of driver discovery, and identify regions that are significantly affected by recurrent breakpoints and recurrent somatic juxtapositions. Our analyses confirm previously reported drivers6,7, raise doubts about others and identify novel candidates, including point mutations in the 5' region of TP53, in the 3' untranslated regions of NFKBIZ and TOB1, focal deletions in BRD4 and rearrangements in the loci of AKR1C genes. We show that although point mutations and structural variants that drive cancer are less frequent in non-coding genes and regulatory sequences than in protein-coding genes, additional examples of these drivers will be found as more cancer genomes become available.

Discovering in vivo cytokine-eQTL interactions from a lupus clinical trial.

BACKGROUND: Cytokines are critical to human disease and are attractive therapeutic targets given their widespread influence on gene regulation and transcription. Defining the downstream regulatory mechanisms influenced by cytokines is central to defining drug and disease mechanisms. One promising strategy is to use interactions between expression quantitative trait loci (eQTLs) and cytokine levels to define target genes and mechanisms. RESULTS: In a clinical trial for anti-IL-6 in patients with systemic lupus erythematosus, we measure interferon (IFN) status, anti-IL-6 drug exposure, and whole blood genome-wide gene expression at three time points. We show that repeat transcriptomic measurements increases the number of cis eQTLs identified compared to using a single time point. We observe a statistically significant enrichment of in vivo eQTL interactions with IFN status and anti-IL-6 drug exposure and find many novel interactions that have not been previously described. Finally, we find transcription factor binding motifs interrupted by eQTL interaction SNPs, which point to key regulatory mediators of these environmental stimuli and therefore potential therapeutic targets for autoimmune diseases. In particular, genes with IFN interactions are enriched for ISRE binding site motifs, while those with anti-IL-6 interactions are enriched for IRF4 motifs. CONCLUSIONS: This study highlights the potential to exploit clinical trial data to discover in vivo eQTL interactions with therapeutically relevant environmental variables.