Microenvironment drives cell state, plasticity, and drug response in pancreatic cancer.

Prognostically relevant RNA expression states exist in pancreatic ductal adenocarcinoma (PDAC), but our understanding of their drivers, stability, and relationship to therapeutic response is limited. To examine these attributes systematically, we profiled metastatic biopsies and matched organoid models at single-cell resolution. In vivo, we identify a new intermediate PDAC transcriptional cell state and uncover distinct site- and state-specific tumor microenvironments (TMEs). Benchmarking models against this reference map, we reveal strong culture-specific biases in cancer cell transcriptional state representation driven by altered TME signals. We restore expression state heterogeneity by adding back in vivo-relevant factors and show plasticity in culture models. Further, we prove that non-genetic modulation of cell state can strongly influence drug responses, uncovering state-specific vulnerabilities. This work provides a broadly applicable framework for aligning cell states across in vivo and ex vivo settings, identifying drivers of transcriptional plasticity and manipulating cell state to target associated vulnerabilities.

Aneuploidy renders cancer cells vulnerable to mitotic checkpoint inhibition.

Selective targeting of aneuploid cells is an attractive strategy for cancer treatment1. However, it is unclear whether aneuploidy generates any clinically relevant vulnerabilities in cancer cells. Here we mapped the aneuploidy landscapes of about 1,000 human cancer cell lines, and analysed genetic and chemical perturbation screens2-9 to identify cellular vulnerabilities associated with aneuploidy. We found that aneuploid cancer cells show increased sensitivity to genetic perturbation of core components of the spindle assembly checkpoint (SAC), which ensures the proper segregation of chromosomes during mitosis10. Unexpectedly, we also found that aneuploid cancer cells were less sensitive than diploid cells to short-term exposure to multiple SAC inhibitors. Indeed, aneuploid cancer cells became increasingly sensitive to inhibition of SAC over time. Aneuploid cells exhibited aberrant spindle geometry and dynamics, and kept dividing when the SAC was inhibited, resulting in the accumulation of mitotic defects, and in unstable and less-fit karyotypes. Therefore, although aneuploid cancer cells could overcome inhibition of SAC more readily than diploid cells, their long-term proliferation was jeopardized. We identified a specific mitotic kinesin, KIF18A, whose activity was perturbed in aneuploid cancer cells. Aneuploid cancer cells were particularly vulnerable to depletion of KIF18A, and KIF18A overexpression restored their response to SAC inhibition. Our results identify a therapeutically relevant, synthetic lethal interaction between aneuploidy and the SAC.

WRN helicase is a synthetic lethal target in microsatellite unstable cancers.

Synthetic lethality-an interaction between two genetic events through which the co-occurrence of these two genetic events leads to cell death, but each event alone does not-can be exploited for cancer therapeutics1. DNA repair processes represent attractive synthetic lethal targets, because many cancers exhibit an impairment of a DNA repair pathway, which can lead to dependence on specific repair proteins2. The success of poly(ADP-ribose) polymerase 1 (PARP-1) inhibitors in cancers with deficiencies in homologous recombination highlights the potential of this approach3. Hypothesizing that other DNA repair defects would give rise to synthetic lethal relationships, we queried dependencies in cancers with microsatellite instability (MSI), which results from deficient DNA mismatch repair. Here we analysed data from large-scale silencing screens using CRISPR-Cas9-mediated knockout and RNA interference, and found that the RecQ DNA helicase WRN was selectively essential in MSI models in vitro and in vivo, yet dispensable in models of cancers that are microsatellite stable. Depletion of WRN induced double-stranded DNA breaks and promoted apoptosis and cell cycle arrest selectively in MSI models. MSI cancer models required the helicase activity of WRN, but not its exonuclease activity. These findings show that WRN is a synthetic lethal vulnerability and promising drug target for MSI cancers.

Next-generation characterization of the Cancer Cell Line Encyclopedia.

Large panels of comprehensively characterized human cancer models, including the Cancer Cell Line Encyclopedia (CCLE), have provided a rigorous framework with which to study genetic variants, candidate targets, and small-molecule and biological therapeutics and to identify new marker-driven cancer dependencies. To improve our understanding of the molecular features that contribute to cancer phenotypes, including drug responses, here we have expanded the characterizations of cancer cell lines to include genetic, RNA splicing, DNA methylation, histone H3 modification, microRNA expression and reverse-phase protein array data for 1,072 cell lines from individuals of various lineages and ethnicities. Integration of these data with functional characterizations such as drug-sensitivity, short hairpin RNA knockdown and CRISPR-Cas9 knockout data reveals potential targets for cancer drugs and associated biomarkers. Together, this dataset and an accompanying public data portal provide a resource for the acceleration of cancer research using model cancer cell lines.

Aneuploidy increases resistance to chemotherapeutics by antagonizing cell division.

Aneuploidy, defined as whole chromosome gains and losses, is associated with poor patient prognosis in many cancer types. However, the condition causes cellular stress and cell cycle delays, foremost in G1 and S phase. Here, we investigate how aneuploidy causes both slow proliferation and poor disease outcome. We test the hypothesis that aneuploidy brings about resistance to chemotherapies because of a general feature of the aneuploid condition-G1 delays. We show that single chromosome gains lead to increased resistance to the frontline chemotherapeutics cisplatin and paclitaxel. Furthermore, G1 cell cycle delays are sufficient to increase chemotherapeutic resistance in euploid cells. Mechanistically, G1 delays increase drug resistance to cisplatin and paclitaxel by reducing their ability to damage DNA and microtubules, respectively. Finally, we show that our findings are clinically relevant. Aneuploidy correlates with slowed proliferation and drug resistance in the Cancer Cell Line Encyclopedia (CCLE) dataset. We conclude that a general and seemingly detrimental effect of aneuploidy, slowed proliferation, provides a selective benefit to cancer cells during chemotherapy treatment.

Bridging the gap between cancer cell line models and tumours using gene expression data.

Cancer cell line models are a cornerstone of cancer research, yet our understanding of how well they represent the molecular features of patient tumours remains limited. Our recent work provides a computational approach to systematically compare large gene expression datasets to better understand which cell lines most closely resemble each tumour type, as well as identify potential gaps in our current cancer models.

Variability and Correlations in Primary Visual Cortical Neurons Driven by Fixational Eye Movements.

UNLABELLED: The ability to distinguish between elements of a sensory neuron's activity that are stimulus independent versus driven by the stimulus is critical for addressing many questions in systems neuroscience. This is typically accomplished by measuring neural responses to repeated presentations of identical stimuli and identifying the trial-variable components of the response as noise. In awake primates, however, small "fixational" eye movements (FEMs) introduce uncontrolled trial-to-trial differences in the visual stimulus itself, potentially confounding this distinction. Here, we describe novel analytical methods that directly quantify the stimulus-driven and stimulus-independent components of visual neuron responses in the presence of FEMs. We apply this approach, combined with precise model-based eye tracking, to recordings from primary visual cortex (V1), finding that standard approaches that ignore FEMs typically miss more than half of the stimulus-driven neural response variance, creating substantial biases in measures of response reliability. We show that these effects are likely not isolated to the particular experimental conditions used here, such as the choice of visual stimulus or spike measurement time window, and thus will be a more general problem for V1 recordings in awake primates. We also demonstrate that measurements of the stimulus-driven and stimulus-independent correlations among pairs of V1 neurons can be greatly biased by FEMs. These results thus illustrate the potentially dramatic impact of FEMs on measures of signal and noise in visual neuron activity and also demonstrate a novel approach for controlling for these eye-movement-induced effects. SIGNIFICANCE STATEMENT: Distinguishing between the signal and noise in a sensory neuron's activity is typically accomplished by measuring neural responses to repeated presentations of an identical stimulus. For recordings from the visual cortex of awake animals, small "fixational" eye movements (FEMs) inevitably introduce trial-to-trial variability in the visual stimulus, potentially confounding such measures. Here, we show that FEMs often have a dramatic impact on several important measures of response variability for neurons in primary visual cortex. We also present an analytical approach for quantifying signal and noise in visual neuron activity in the presence of FEMs. These results thus highlight the importance of controlling for FEMs in studies of visual neuron function, and demonstrate novel methods for doing so.

Inferring Cortical Variability from Local Field Potentials.

The responses of sensory neurons can be quite different to repeated presentations of the same stimulus. Here, we demonstrate a direct link between the trial-to-trial variability of cortical neuron responses and network activity that is reflected in local field potentials (LFPs). Spikes and LFPs were recorded with a multielectrode array from the middle temporal (MT) area of the visual cortex of macaques during the presentation of continuous optic flow stimuli. A maximum likelihood-based modeling framework was used to predict single-neuron spiking responses using the stimulus, the LFPs, and the activity of other recorded neurons. MT neuron responses were strongly linked to gamma oscillations (maximum at 40 Hz) as well as to lower-frequency delta oscillations (1-4 Hz), with consistent phase preferences across neurons. The predicted modulation associated with the LFP was largely complementary to that driven by visual stimulation, as well as the activity of other neurons, and accounted for nearly half of the trial-to-trial variability in the spiking responses. Moreover, the LFP model predictions accurately captured the temporal structure of noise correlations between pairs of simultaneously recorded neurons, and explained the variation in correlation magnitudes observed across the population. These results therefore identify signatures of network activity related to the variability of cortical neuron responses, and suggest their central role in sensory cortical function. SIGNIFICANCE STATEMENT: The function of sensory neurons is nearly always cast in terms of representing sensory stimuli. However, recordings from visual cortex in awake animals show that a large fraction of neural activity is not predictable from the stimulus. We show that this variability is predictable given the simultaneously recorded measures of network activity, local field potentials. A model that combines elements of these signals with the stimulus processing of the neuron can predict neural responses dramatically better than current models, and can predict the structure of correlations across the cortical population. In identifying ways to understand stimulus processing in the context of ongoing network activity, this work thus provides a foundation to understand the role of sensory cortex in combining sensory and cognitive variables.

Inferring nonlinear neuronal computation based on physiologically plausible inputs.

The computation represented by a sensory neuron's response to stimuli is constructed from an array of physiological processes both belonging to that neuron and inherited from its inputs. Although many of these physiological processes are known to be nonlinear, linear approximations are commonly used to describe the stimulus selectivity of sensory neurons (i.e., linear receptive fields). Here we present an approach for modeling sensory processing, termed the Nonlinear Input Model (NIM), which is based on the hypothesis that the dominant nonlinearities imposed by physiological mechanisms arise from rectification of a neuron's inputs. Incorporating such 'upstream nonlinearities' within the standard linear-nonlinear (LN) cascade modeling structure implicitly allows for the identification of multiple stimulus features driving a neuron's response, which become directly interpretable as either excitatory or inhibitory. Because its form is analogous to an integrate-and-fire neuron receiving excitatory and inhibitory inputs, model fitting can be guided by prior knowledge about the inputs to a given neuron, and elements of the resulting model can often result in specific physiological predictions. Furthermore, by providing an explicit probabilistic model with a relatively simple nonlinear structure, its parameters can be efficiently optimized and appropriately regularized. Parameter estimation is robust and efficient even with large numbers of model components and in the context of high-dimensional stimuli with complex statistical structure (e.g. natural stimuli). We describe detailed methods for estimating the model parameters, and illustrate the advantages of the NIM using a range of example sensory neurons in the visual and auditory systems. We thus present a modeling framework that can capture a broad range of nonlinear response functions while providing physiologically interpretable descriptions of neural computation.

Spontaneous persistent activity and inactivity in vivo reveals differential cortico-entorhinal functional connectivity.

Understanding the functional connectivity between brain regions and its emergent dynamics is a central challenge. Here we present a theory-experiment hybrid approach involving iteration between a minimal computational model and in vivo electrophysiological measurements. Our model not only predicted spontaneous persistent activity (SPA) during Up-Down-State oscillations, but also inactivity (SPI), which has never been reported. These were confirmed in vivo in the membrane potential of neurons, especially from layer 3 of the medial and lateral entorhinal cortices. The data was then used to constrain two free parameters, yielding a unique, experimentally determined model for each neuron. Analytic and computational analysis of the model generated a dozen quantitative predictions about network dynamics, which were all confirmed in vivo to high accuracy. Our technique predicted functional connectivity; e. g. the recurrent excitation is stronger in the medial than lateral entorhinal cortex. This too was confirmed with connectomics data. This technique uncovers how differential cortico-entorhinal dialogue generates SPA and SPI, which could form an energetically efficient working-memory substrate and influence the consolidation of memories during sleep. More broadly, our procedure can reveal the functional connectivity of large networks and a theory of their emergent dynamics.

Germline variation contributes to false negatives in CRISPR-based experiments with varying burden across ancestries.

Reducing disparities is vital for equitable access to precision treatments in cancer. Socioenvironmental factors are a major driver of disparities, but differences in genetic variation likely also contribute. The impact of genetic ancestry on prioritization of cancer targets in drug discovery pipelines has not been systematically explored due to the absence of pre-clinical data at the appropriate scale. Here, we analyze data from 611 genome-scale CRISPR/Cas9 viability experiments in human cell line models to identify ancestry-associated genetic dependencies essential for cell survival. Surprisingly, we find that most putative associations between ancestry and dependency arise from artifacts related to germline variants. Our analysis suggests that for 1.2-2.5% of guides, germline variants in sgRNA targeting sequences reduce cutting by the CRISPR/Cas9 nuclease, disproportionately affecting cell models derived from individuals of recent African descent. We propose three approaches to mitigate this experimental bias, enabling the scientific community to address these disparities.

The effects of GluA1 deletion on the hippocampal population code for position.

The GluA1 subunit of AMPA receptors (AMPARs) is critical for hippocampal synaptic transmission and plasticity. Here, we measured the activity of single units from the CA1 region of the hippocampus while GluA1 knock-out (GluA1⁻/⁻) and wild-type (WT) mice traversed a linear track. Although overall firing rates were similar, GluA1⁻/⁻ neurons were more likely to spike in bursts, but at lower burst frequencies, compared with WT neurons. GluA1⁻/⁻ neurons showed large reductions in all measures of spatial and directional selectivity compared with WT neurons. Consistent with these alterations of single-neuron properties, the accuracy of the population code for position was substantially reduced in GluA1⁻/⁻, yet it is predicted to approach the accuracy of WT with increasing population size. The absolute representation of space, independent of movement direction, was greatly diminished in GluA1⁻/⁻ mice and is predicted to remain reduced even for larger populations. Finally, we found that the rate maps of GluA1⁻/⁻ neurons showed increased trial-by-trial variability but reduced experiential plasticity compared with the WT. These results reveal the critical contribution of GluA1-containing AMPARs to individual place cells and the hippocampal population code for space, which could explain the selective behavioral impairments observed in these mice.

Partial gene suppression improves identification of cancer vulnerabilities when CRISPR-Cas9 knockout is pan-lethal.

BACKGROUND: Hundreds of functional genomic screens have been performed across a diverse set of cancer contexts, as part of efforts such as the Cancer Dependency Map, to identify gene dependencies-genes whose loss of function reduces cell viability or fitness. Recently, large-scale screening efforts have shifted from RNAi to CRISPR-Cas9, due to superior efficacy and specificity. However, many effective oncology drugs only partially inhibit their protein targets, leading us to question whether partial suppression of genes using RNAi could reveal cancer vulnerabilities that are missed by complete knockout using CRISPR-Cas9. Here, we compare CRISPR-Cas9 and RNAi dependency profiles of genes across approximately 400 matched cancer cell lines. RESULTS: We find that CRISPR screens accurately identify more gene dependencies per cell line, but the majority of each cell line's dependencies are part of a set of 1867 genes that are shared dependencies across the entire collection (pan-lethals). While RNAi knockdown of about 30% of these genes is also pan-lethal, approximately 50% have selective dependency patterns across cell lines, suggesting they could still be cancer vulnerabilities. The accuracy of the unique RNAi selectivity is supported by associations to multi-omics profiles, drug sensitivity, and other expected co-dependencies. CONCLUSIONS: Incorporating RNAi data for genes that are pan-lethal knockouts facilitates the discovery of a wider range of gene targets than could be detected using the CRISPR dataset alone. This can aid in the interpretation of contrasting results obtained from CRISPR and RNAi screens and reinforce the importance of partial gene suppression methods in building a cancer dependency map.

A Ubiquitination Cascade Regulating the Integrated Stress Response and Survival in Carcinomas.

Systematic identification of signaling pathways required for the fitness of cancer cells will facilitate the development of new cancer therapies. We used gene essentiality measurements in 1,086 cancer cell lines to identify selective coessentiality modules and found that a ubiquitin ligase complex composed of UBA6, BIRC6, KCMF1, and UBR4 is required for the survival of a subset of epithelial tumors that exhibit a high degree of aneuploidy. Suppressing BIRC6 in cell lines that are dependent on this complex led to a substantial reduction in cell fitness in vitro and potent tumor regression in vivo. Mechanistically, BIRC6 suppression resulted in selective activation of the integrated stress response (ISR) by stabilization of the heme-regulated inhibitor, a direct ubiquitination target of the UBA6/BIRC6/KCMF1/UBR4 complex. These observations uncover a novel ubiquitination cascade that regulates ISR and highlight the potential of ISR activation as a new therapeutic strategy. SIGNIFICANCE: We describe the identification of a heretofore unrecognized ubiquitin ligase complex that prevents the aberrant activation of the ISR in a subset of cancer cells. This provides a novel insight on the regulation of ISR and exposes a therapeutic opportunity to selectively eliminate these cancer cells. See related commentary Leli and Koumenis, p. 535. This article is highlighted in the In This Issue feature, p. 517.

Genome-scale functional genomics identify genes preferentially essential for multiple myeloma cells compared to other neoplasias.

Clinical progress in multiple myeloma (MM), an incurable plasma cell (PC) neoplasia, has been driven by therapies that have limited applications beyond MM/PC neoplasias and do not target specific oncogenic mutations in MM. Instead, these agents target pathways critical for PC biology yet largely dispensable for malignant or normal cells of most other lineages. Here we systematically characterized the lineage-preferential molecular dependencies of MM through genome-scale clustered regularly interspaced short palindromic repeats (CRISPR) studies in 19 MM versus hundreds of non-MM lines and identified 116 genes whose disruption more significantly affects MM cell fitness compared with other malignancies. These genes, some known, others not previously linked to MM, encode transcription factors, chromatin modifiers, endoplasmic reticulum components, metabolic regulators or signaling molecules. Most of these genes are not among the top amplified, overexpressed or mutated in MM. Functional genomics approaches thus define new therapeutic targets in MM not readily identifiable by standard genomic, transcriptional or epigenetic profiling analyses.

Sparse dictionary learning recovers pleiotropy from human cell fitness screens.

In high-throughput functional genomic screens, each gene product is commonly assumed to exhibit a singular biological function within a defined protein complex or pathway. In practice, a single gene perturbation may induce multiple cascading functional outcomes, a genetic principle known as pleiotropy. Here, we model pleiotropy in fitness screen collections by representing each gene perturbation as the sum of multiple perturbations of biological functions, each harboring independent fitness effects inferred empirically from the data. Our approach (Webster) recovered pleiotropic functions for DNA damage proteins from genotoxic fitness screens, untangled distinct signaling pathways upstream of shared effector proteins from cancer cell fitness screens, and predicted the stoichiometry of an unknown protein complex subunit from fitness data alone. Modeling compound sensitivity profiles in terms of genetic functions recovered compound mechanisms of action. Our approach establishes a sparse approximation mechanism for unraveling complex genetic architectures underlying high-dimensional gene perturbation readouts.

Phosphate dysregulation via the XPR1-KIDINS220 protein complex is a therapeutic vulnerability in ovarian cancer.

Despite advances in precision medicine, the clinical prospects for patients with ovarian and uterine cancers have not substantially improved. Here, we analyzed genome-scale CRISPR-Cas9 loss-of-function screens across 851 human cancer cell lines and found that frequent overexpression of SLC34A2-encoding a phosphate importer-is correlated with sensitivity to loss of the phosphate exporter XPR1, both in vitro and in vivo. In patient-derived tumor samples, we observed frequent PAX8-dependent overexpression of SLC34A2, XPR1 copy number amplifications and XPR1 messenger RNA overexpression. Mechanistically, in SLC34A2-high cancer cell lines, genetic or pharmacologic inhibition of XPR1-dependent phosphate efflux leads to the toxic accumulation of intracellular phosphate. Finally, we show that XPR1 requires the novel partner protein KIDINS220 for proper cellular localization and activity, and that disruption of this protein complex results in acidic "vacuolar" structures preceding cell death. These data point to the XPR1-KIDINS220 complex and phosphate dysregulation as a therapeutic vulnerability in ovarian cancer.

FOXR2 Is an Epigenetically Regulated Pan-Cancer Oncogene That Activates ETS Transcriptional Circuits.

Forkhead box R2 (FOXR2) is a forkhead transcription factor located on the X chromosome whose expression is normally restricted to the testis. In this study, we performed a pan-cancer analysis of FOXR2 activation across more than 10,000 adult and pediatric cancer samples and found FOXR2 to be aberrantly upregulated in 70% of all cancer types and 8% of all individual tumors. The majority of tumors (78%) aberrantly expressed FOXR2 through a previously undescribed epigenetic mechanism that involves hypomethylation of a novel promoter, which was functionally validated as necessary for FOXR2 expression and proliferation in FOXR2-expressing cancer cells. FOXR2 promoted tumor growth across multiple cancer lineages and co-opted ETS family transcription circuits across cancers. Taken together, this study identifies FOXR2 as a potent and ubiquitous oncogene that is epigenetically activated across the majority of human cancers. The identification of hijacking of ETS transcription circuits by FOXR2 extends the mechanisms known to active ETS transcription factors and highlights how transcription factor families cooperate to enhance tumorigenesis. SIGNIFICANCE: This work identifies a novel promoter that drives aberrant FOXR2 expression and delineates FOXR2 as a pan-cancer oncogene that specifically activates ETS transcriptional circuits across human cancers. See related commentary by Liu and Northcott, p. 2977.

Global computational alignment of tumor and cell line transcriptional profiles.

Cell lines are key tools for preclinical cancer research, but it remains unclear how well they represent patient tumor samples. Direct comparisons of tumor and cell line transcriptional profiles are complicated by several factors, including the variable presence of normal cells in tumor samples. We thus develop an unsupervised alignment method (Celligner) and apply it to integrate several large-scale cell line and tumor RNA-Seq datasets. Although our method aligns the majority of cell lines with tumor samples of the same cancer type, it also reveals large differences in tumor similarity across cell lines. Using this approach, we identify several hundred cell lines from diverse lineages that present a more mesenchymal and undifferentiated transcriptional state and that exhibit distinct chemical and genetic dependencies. Celligner could be used to guide the selection of cell lines that more closely resemble patient tumors and improve the clinical translation of insights gained from cell lines.