Ovarian cancer-derived IL-4 promotes immunotherapy resistance.

Ovarian cancer is resistant to immunotherapy, and this is influenced by the immunosuppressed tumor microenvironment (TME) dominated by macrophages. Resistance is also affected by intratumoral heterogeneity, whose development is poorly understood. To identify regulators of ovarian cancer immunity, we employed a spatial functional genomics screen (Perturb-map), focused on receptor/ligands hypothesized to be involved in tumor-macrophage communication. Perturb-map recapitulated tumor heterogeneity and revealed that interleukin-4 (IL-4) promotes resistance to anti-PD-1. We find ovarian cancer cells are the key source of IL-4, which directs the formation of an immunosuppressive TME via macrophage control. IL-4 loss was not compensated by nearby IL-4-expressing clones, revealing short-range regulation of TME composition dictating tumor evolution. Our studies show heterogeneous TMEs can emerge from localized altered expression of cancer-derived cytokines/chemokines that establish immune-rich and immune-excluded neighborhoods, which drive clone selection and immunotherapy resistance. They also demonstrate the potential of targeting IL-4 signaling to enhance ovarian cancer response to immunotherapy.

Multiplexed 3D atlas of state transitions and immune interaction in colorectal cancer.

Advanced solid cancers are complex assemblies of tumor, immune, and stromal cells characterized by high intratumoral variation. We use highly multiplexed tissue imaging, 3D reconstruction, spatial statistics, and machine learning to identify cell types and states underlying morphological features of known diagnostic and prognostic significance in colorectal cancer. Quantitation of these features in high-plex marker space reveals recurrent transitions from one tumor morphology to the next, some of which are coincident with long-range gradients in the expression of oncogenes and epigenetic regulators. At the tumor invasive margin, where tumor, normal, and immune cells compete, T cell suppression involves multiple cell types and 3D imaging shows that seemingly localized 2D features such as tertiary lymphoid structures are commonly interconnected and have graded molecular properties. Thus, while cancer genetics emphasizes the importance of discrete changes in tumor state, whole-specimen imaging reveals large-scale morphological and molecular gradients analogous to those in developing tissues.

Nociceptor-immune interactomes reveal insult-specific immune signatures of pain.

Inflammatory pain results from the heightened sensitivity and reduced threshold of nociceptor sensory neurons due to exposure to inflammatory mediators. However, the cellular and transcriptional diversity of immune cell and sensory neuron types makes it challenging to decipher the immune mechanisms underlying pain. Here we used single-cell transcriptomics to determine the immune gene signatures associated with pain development in three skin inflammatory pain models in mice: zymosan injection, skin incision and ultraviolet burn. We found that macrophage and neutrophil recruitment closely mirrored the kinetics of pain development and identified cell-type-specific transcriptional programs associated with pain and its resolution. Using a comprehensive list of potential interactions mediated by receptors, ligands, ion channels and metabolites to generate injury-specific neuroimmune interactomes, we also uncovered that thrombospondin-1 upregulated by immune cells upon injury inhibited nociceptor sensitization. This study lays the groundwork for identifying the neuroimmune axes that modulate pain in diverse disease contexts.

Obesity Shapes Metabolism in the Tumor Microenvironment to Suppress Anti-Tumor Immunity.

Obesity is a major cancer risk factor, but how differences in systemic metabolism change the tumor microenvironment (TME) and impact anti-tumor immunity is not understood. Here, we demonstrate that high-fat diet (HFD)-induced obesity impairs CD8+ T cell function in the murine TME, accelerating tumor growth. We generate a single-cell resolution atlas of cellular metabolism in the TME, detailing how it changes with diet-induced obesity. We find that tumor and CD8+ T cells display distinct metabolic adaptations to obesity. Tumor cells increase fat uptake with HFD, whereas tumor-infiltrating CD8+ T cells do not. These differential adaptations lead to altered fatty acid partitioning in HFD tumors, impairing CD8+ T cell infiltration and function. Blocking metabolic reprogramming by tumor cells in obese mice improves anti-tumor immunity. Analysis of human cancers reveals similar transcriptional changes in CD8+ T cell markers, suggesting interventions that exploit metabolism to improve cancer immunotherapy.

Vascular Disease and Thrombosis in SARS-CoV-2-Infected Rhesus Macaques.

The COVID-19 pandemic has led to extensive morbidity and mortality throughout the world. Clinical features that drive SARS-CoV-2 pathogenesis in humans include inflammation and thrombosis, but the mechanistic details underlying these processes remain to be determined. In this study, we demonstrate endothelial disruption and vascular thrombosis in histopathologic sections of lungs from both humans and rhesus macaques infected with SARS-CoV-2. To define key molecular pathways associated with SARS-CoV-2 pathogenesis in macaques, we performed transcriptomic analyses of bronchoalveolar lavage and peripheral blood and proteomic analyses of serum. We observed macrophage infiltrates in lung and upregulation of macrophage, complement, platelet activation, thrombosis, and proinflammatory markers, including C-reactive protein, MX1, IL-6, IL-1, IL-8, TNFα, and NF-κB. These results suggest a model in which critical interactions between inflammatory and thrombosis pathways lead to SARS-CoV-2-induced vascular disease. Our findings suggest potential therapeutic targets for COVID-19.

The Human Tumor Atlas Network: Charting Tumor Transitions across Space and Time at Single-Cell Resolution.

Crucial transitions in cancer-including tumor initiation, local expansion, metastasis, and therapeutic resistance-involve complex interactions between cells within the dynamic tumor ecosystem. Transformative single-cell genomics technologies and spatial multiplex in situ methods now provide an opportunity to interrogate this complexity at unprecedented resolution. The Human Tumor Atlas Network (HTAN), part of the National Cancer Institute (NCI) Cancer Moonshot Initiative, will establish a clinical, experimental, computational, and organizational framework to generate informative and accessible three-dimensional atlases of cancer transitions for a diverse set of tumor types. This effort complements both ongoing efforts to map healthy organs and previous large-scale cancer genomics approaches focused on bulk sequencing at a single point in time. Generating single-cell, multiparametric, longitudinal atlases and integrating them with clinical outcomes should help identify novel predictive biomarkers and features as well as therapeutically relevant cell types, cell states, and cellular interactions across transitions. The resulting tumor atlases should have a profound impact on our understanding of cancer biology and have the potential to improve cancer detection, prevention, and therapeutic discovery for better precision-medicine treatments of cancer patients and those at risk for cancer.

A Cancer Cell Program Promotes T Cell Exclusion and Resistance to Checkpoint Blockade.

Immune checkpoint inhibitors (ICIs) produce durable responses in some melanoma patients, but many patients derive no clinical benefit, and the molecular underpinnings of such resistance remain elusive. Here, we leveraged single-cell RNA sequencing (scRNA-seq) from 33 melanoma tumors and computational analyses to interrogate malignant cell states that promote immune evasion. We identified a resistance program expressed by malignant cells that is associated with T cell exclusion and immune evasion. The program is expressed prior to immunotherapy, characterizes cold niches in situ, and predicts clinical responses to anti-PD-1 therapy in an independent cohort of 112 melanoma patients. CDK4/6-inhibition represses this program in individual malignant cells, induces senescence, and reduces melanoma tumor outgrowth in mouse models in vivo when given in combination with immunotherapy. Our study provides a high-resolution landscape of ICI-resistant cell states, identifies clinically predictive signatures, and suggests new therapeutic strategies to overcome immunotherapy resistance.

Combination Cancer Therapy Can Confer Benefit via Patient-to-Patient Variability without Drug Additivity or Synergy.

Combination cancer therapies aim to improve the probability and magnitude of therapeutic responses and reduce the likelihood of acquired resistance in an individual patient. However, drugs are tested in clinical trials on genetically diverse patient populations. We show here that patient-to-patient variability and independent drug action are sufficient to explain the superiority of many FDA-approved drug combinations in the absence of drug synergy or additivity. This is also true for combinations tested in patient-derived tumor xenografts. In a combination exhibiting independent drug action, each patient benefits solely from the drug to which his or her tumor is most sensitive, with no added benefit from other drugs. Even when drug combinations exhibit additivity or synergy in pre-clinical models, patient-to-patient variability and low cross-resistance make independent action the dominant mechanism in clinical populations. This insight represents a different way to interpret trial data and a different way to design combination therapies.

AKT and EZH2 inhibitors kill TNBCs by hijacking mechanisms of involution.

Triple-negative breast cancer (TNBC) is the most aggressive breast cancer subtype and has the highest rate of recurrence1. The predominant standard of care for advanced TNBC is systemic chemotherapy with or without immunotherapy; however, responses are typically short lived1,2. Thus, there is an urgent need to develop more effective treatments. Components of the PI3K pathway represent plausible therapeutic targets; more than 70% of TNBCs have alterations in PIK3CA, AKT1 or PTEN3-6. However, in contrast to hormone-receptor-positive tumours, it is still unclear whether or how triple-negative disease will respond to PI3K pathway inhibitors7. Here we describe a promising AKT-inhibitor-based therapeutic combination for TNBC. Specifically, we show that AKT inhibitors synergize with agents that suppress the histone methyltransferase EZH2 and promote robust tumour regression in multiple TNBC models in vivo. AKT and EZH2 inhibitors exert these effects by first cooperatively driving basal-like TNBC cells into a more differentiated, luminal-like state, which cannot be effectively induced by either agent alone. Once TNBCs are differentiated, these agents kill them by hijacking signals that normally drive mammary gland involution. Using a machine learning approach, we developed a classifier that can be used to predict sensitivity. Together, these findings identify a promising therapeutic strategy for this highly aggressive tumour type and illustrate how deregulated epigenetic enzymes can insulate tumours from oncogenic vulnerabilities. These studies also reveal how developmental tissue-specific cell death pathways may be co-opted for therapeutic benefit.

A noncanonical Frizzled2 pathway regulates epithelial-mesenchymal transition and metastasis.

Wnt signaling plays a critical role in embryonic development, and genetic aberrations in this network have been broadly implicated in colorectal cancer. We find that the Wnt receptor Frizzled2 (Fzd2) and its ligands Wnt5a/b are elevated in metastatic liver, lung, colon, and breast cancer cell lines and in high-grade tumors and that their expression correlates with markers of epithelial-mesenchymal transition (EMT). Pharmacologic and genetic perturbations reveal that Fzd2 drives EMT and cell migration through a previously unrecognized, noncanonical pathway that includes Fyn and Stat3. A gene signature regulated by this pathway predicts metastasis and overall survival in patients. We have developed an antibody to Fzd2 that reduces cell migration and invasion and inhibits tumor growth and metastasis in xenografts. We propose that targeting this pathway could provide benefit for patients with tumors expressing high levels of Fzd2 and Wnt5a/b.

S. cerevisiae chromosomes biorient via gradual resolution of syntely between S phase and anaphase.

Following DNA replication, eukaryotic cells must biorient all sister chromatids prior to cohesion cleavage at anaphase. In animal cells, sister chromatids gradually biorient during prometaphase, but current models of mitosis in S. cerevisiae assume that biorientation is established shortly after S phase. This assumption is based on the observation of a bilobed distribution of yeast kinetochores early in mitosis and suggests fundamental differences between yeast mitosis and mitosis in animal cells. By applying super-resolution imaging methods, we show that yeast and animal cells share the key property of gradual and stochastic chromosome biorientation. The characteristic bilobed distribution of yeast kinetochores, hitherto considered synonymous for biorientation, arises from kinetochores in mixed attachment states to microtubules, the length of which discriminates bioriented from syntelic attachments. Our results offer a revised view of mitotic progression in S. cerevisiae that augments the relevance of mechanistic information obtained in this powerful genetic system for mammalian mitosis.

Quality control for single-cell analysis of high-plex tissue profiles using CyLinter.

Tumors are complex assemblies of cellular and acellular structures patterned on spatial scales from microns to centimeters. Study of these assemblies has advanced dramatically with the introduction of high-plex spatial profiling. Image-based profiling methods reveal the intensities and spatial distributions of 20-100 proteins at subcellular resolution in 103-107 cells per specimen. Despite extensive work on methods for extracting single-cell data from these images, all tissue images contain artifacts such as folds, debris, antibody aggregates, optical aberrations and image processing errors that arise from imperfections in specimen preparation, data acquisition, image assembly and feature extraction. Here we show that these artifacts dramatically impact single-cell data analysis, obscuring meaningful biological interpretation. We describe an interactive quality control software tool, CyLinter, that identifies and removes data associated with imaging artifacts. CyLinter greatly improves single-cell analysis, especially for archival specimens sectioned many years before data collection, such as those from clinical trials.

OpenFold: retraining AlphaFold2 yields new insights into its learning mechanisms and capacity for generalization.

AlphaFold2 revolutionized structural biology with the ability to predict protein structures with exceptionally high accuracy. Its implementation, however, lacks the code and data required to train new models. These are necessary to (1) tackle new tasks, like protein-ligand complex structure prediction, (2) investigate the process by which the model learns and (3) assess the model's capacity to generalize to unseen regions of fold space. Here we report OpenFold, a fast, memory efficient and trainable implementation of AlphaFold2. We train OpenFold from scratch, matching the accuracy of AlphaFold2. Having established parity, we find that OpenFold is remarkably robust at generalizing even when the size and diversity of its training set is deliberately limited, including near-complete elisions of classes of secondary structure elements. By analyzing intermediate structures produced during training, we also gain insights into the hierarchical manner in which OpenFold learns to fold. In sum, our studies demonstrate the power and utility of OpenFold, which we believe will prove to be a crucial resource for the protein modeling community.

Sequential application of anticancer drugs enhances cell death by rewiring apoptotic signaling networks.

Crosstalk and complexity within signaling pathways and their perturbation by oncogenes limit component-by-component approaches to understanding human disease. Network analysis of how normal and oncogenic signaling can be rewired by drugs may provide opportunities to target tumors with high specificity and efficacy. Using targeted inhibition of oncogenic signaling pathways, combined with DNA-damaging chemotherapy, we report that time-staggered EGFR inhibition, but not simultaneous coadministration, dramatically sensitizes a subset of triple-negative breast cancer cells to genotoxic drugs. Systems-level analysis-using high-density time-dependent measurements of signaling networks, gene expression profiles, and cell phenotypic responses in combination with mathematical modeling-revealed an approach for altering the intrinsic state of the cell through dynamic rewiring of oncogenic signaling pathways. This process converts these cells to a less tumorigenic state that is more susceptible to DNA damage-induced cell death by reactivation of an extrinsic apoptotic pathway whose function is suppressed in the oncogene-addicted state.

SpatialCells: automated profiling of tumor microenvironments with spatially resolved multiplexed single-cell data.

Cancer is a complex cellular ecosystem where malignant cells coexist and interact with immune, stromal and other cells within the tumor microenvironment (TME). Recent technological advancements in spatially resolved multiplexed imaging at single-cell resolution have led to the generation of large-scale and high-dimensional datasets from biological specimens. This underscores the necessity for automated methodologies that can effectively characterize molecular, cellular and spatial properties of TMEs for various malignancies. This study introduces SpatialCells, an open-source software package designed for region-based exploratory analysis and comprehensive characterization of TMEs using multiplexed single-cell data. The source code and tutorials are available at https://semenovlab.github.io/SpatialCells. SpatialCells efficiently streamlines the automated extraction of features from multiplexed single-cell data and can process samples containing millions of cells. Thus, SpatialCells facilitates subsequent association analyses and machine learning predictions, making it an essential tool in advancing our understanding of tumor growth, invasion and metastasis.

Measuring and modeling apoptosis in single cells.

Cell death plays an essential role in the development of tissues and organisms, the etiology of disease, and the responses of cells to therapeutic drugs. Here we review progress made over the last decade in using mathematical models and quantitative, often single-cell, data to study apoptosis. We discuss the delay that follows exposure of cells to prodeath stimuli, control of mitochondrial outer membrane permeabilization, switch-like activation of effector caspases, and variability in the timing and probability of death from one cell to the next. Finally, we discuss challenges facing the fields of biochemical modeling and systems pharmacology.

MCMICRO: a scalable, modular image-processing pipeline for multiplexed tissue imaging.

Highly multiplexed tissue imaging makes detailed molecular analysis of single cells possible in a preserved spatial context. However, reproducible analysis of large multichannel images poses a substantial computational challenge. Here, we describe a modular and open-source computational pipeline, MCMICRO, for performing the sequential steps needed to transform whole-slide images into single-cell data. We demonstrate the use of MCMICRO on tissue and tumor images acquired using multiple imaging platforms, thereby providing a solid foundation for the continued development of tissue imaging software.

Differentiable biology: using deep learning for biophysics-based and data-driven modeling of molecular mechanisms.

Deep learning using neural networks relies on a class of machine-learnable models constructed using 'differentiable programs'. These programs can combine mathematical equations specific to a particular domain of natural science with general-purpose, machine-learnable components trained on experimental data. Such programs are having a growing impact on molecular and cellular biology. In this Perspective, we describe an emerging 'differentiable biology' in which phenomena ranging from the small and specific (for example, one experimental assay) to the broad and complex (for example, protein folding) can be modeled effectively and efficiently, often by exploiting knowledge about basic natural phenomena to overcome the limitations of sparse, incomplete and noisy data. By distilling differentiable biology into a small set of conceptual primitives and illustrative vignettes, we show how it can help to address long-standing challenges in integrating multimodal data from diverse experiments across biological scales. This promises to benefit fields as diverse as biophysics and functional genomics.

Mechanistic model of MAPK signaling reveals how allostery and rewiring contribute to drug resistance.

BRAF is prototypical of oncogenes that can be targeted therapeutically and the treatment of BRAFV600E melanomas with RAF and MEK inhibitors results in rapid tumor regression. However, drug-induced rewiring generates a drug adapted state thought to be involved in acquired resistance and disease recurrence. In this article, we study mechanisms of adaptive rewiring in BRAFV600E melanoma cells using an energy-based implementation of ordinary differential equation (ODE) modeling in combination with proteomic, transcriptomic and imaging data. We develop a method for causal tracing of ODE models and identify two parallel MAPK reaction channels that are differentially sensitive to RAF and MEK inhibitors due to differences in protein oligomerization and drug binding. We describe how these channels, and timescale separation between immediate-early signaling and transcriptional feedback, create a state in which the RAS-regulated MAPK channel can be activated by growth factors under conditions in which the BRAFV600E -driven channel is fully inhibited. Further development of the approaches in this article is expected to yield a unified model of adaptive drug resistance in melanoma.

Automated assembly of molecular mechanisms at scale from text mining and curated databases.

The analysis of omic data depends on machine-readable information about protein interactions, modifications, and activities as found in protein interaction networks, databases of post-translational modifications, and curated models of gene and protein function. These resources typically depend heavily on human curation. Natural language processing systems that read the primary literature have the potential to substantially extend knowledge resources while reducing the burden on human curators. However, machine-reading systems are limited by high error rates and commonly generate fragmentary and redundant information. Here, we describe an approach to precisely assemble molecular mechanisms at scale using multiple natural language processing systems and the Integrated Network and Dynamical Reasoning Assembler (INDRA). INDRA identifies full and partial overlaps in information extracted from published papers and pathway databases, uses predictive models to improve the reliability of machine reading, and thereby assembles individual pieces of information into non-redundant and broadly usable mechanistic knowledge. Using INDRA to create high-quality corpora of causal knowledge we show it is possible to extend protein-protein interaction databases and explain co-dependencies in the Cancer Dependency Map.