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.

Computational estimation of quality and clinical relevance of cancer cell lines.

Immortal cancer cell lines (CCLs) are the most widely used system for investigating cancer biology and for the preclinical development of oncology therapies. Pharmacogenomic and genome-wide editing screenings have facilitated the discovery of clinically relevant gene-drug interactions and novel therapeutic targets via large panels of extensively characterised CCLs. However, tailoring pharmacological strategies in a precision medicine context requires bridging the existing gaps between tumours and in vitro models. Indeed, intrinsic limitations of CCLs such as misidentification, the absence of tumour microenvironment and genetic drift have highlighted the need to identify the most faithful CCLs for each primary tumour while addressing their heterogeneity, with the development of new models where necessary. Here, we discuss the most significant limitations of CCLs in representing patient features, and we review computational methods aiming at systematically evaluating the suitability of CCLs as tumour proxies and identifying the best patient representative in vitro models. Additionally, we provide an overview of the applications of these methods to more complex models and discuss future machine-learning-based directions that could resolve some of the arising discrepancies.

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.

Pan-ERBB kinase inhibition augments CDK4/6 inhibitor efficacy in oesophageal squamous cell carcinoma.

OBJECTIVE: Oesophageal squamous cell carcinoma (OSCC), like other squamous carcinomas, harbour highly recurrent cell cycle pathway alterations, especially hyperactivation of the CCND1/CDK4/6 axis, raising the potential for use of existing CDK4/6 inhibitors in these cancers. Although CDK4/6 inhibition has shown striking success when combined with endocrine therapy in oestrogen receptor positive breast cancer, CDK4/6 inhibitor palbociclib monotherapy has not revealed evidence of efficacy to date in OSCC clinical studies. Herein, we sought to elucidate the identification of key dependencies in OSCC as a foundation for the selection of targets whose blockade could be combined with CDK4/6 inhibition. DESIGN: We combined large-scale genomic dependency and pharmaceutical screening datasets with preclinical cell line models, to identified potential combination therapies in squamous cell cancer. RESULTS: We identified sensitivity to inhibitors to the ERBB family of receptor kinases, results clearly extending beyond the previously described minority of tumours with EGFR amplification/dependence, specifically finding a subset of OSCCs with dual dependence on ERBB3 and ERBB2. Subsequently. we demonstrated marked efficacy of combined pan-ERBB and CDK4/6 inhibition in vitro and in vivo. Furthermore, we demonstrated that squamous lineage transcription factor KLF5 facilitated activation of ERBBs in OSCC. CONCLUSION: These results provide clear rationale for development of combined ERBB and CDK4/6 inhibition in these cancers and raises the potential for KLF5 expression as a candidate biomarker to guide the use of these agents. These data suggested that by combining existing Food and Drug Administration (FDA)-approved agents, we have the capacity to improve therapy for OSCC and other squamous cancer.

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.

BRD-810 is a highly selective MCL1 inhibitor with optimized in vivo clearance and robust efficacy in solid and hematological tumor models.

The MCL1 gene is frequently amplified in cancer and codes for the antiapoptotic protein myeloid cell leukemia 1 (MCL1), which confers resistance to the current standard of care. Therefore, MCL1 is an attractive anticancer target. Here we describe BRD-810 as a potent and selective MCL1 inhibitor and its key design principle of rapid systemic clearance to potentially minimize area under the curve-driven toxicities associated with MCL1 inhibition. BRD-810 induced rapid cell killing within 4 h in vitro but, in the same 4-h window, had no impact on cell viability or troponin I release in human induced pluripotent stem cell-derived cardiomyocytes, even at suprapharmacologic concentrations. In vivo BRD-810 induced efficacy in xenograft hematological and solid tumor models despite the short residence time of BRD-810 in plasma. In totality, our data support the hypothesis that short-term inhibition of MCL1 with BRD-810 can induce apoptosis in tumor cells while maintaining an acceptable safety profile. We, therefore, intend to advance BRD-810 to clinical trials.

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.

Net-Zero CO2 by 2050 Scenarios for the United States in the Energy Modeling Forum 37 Study.

The Energy Modeling Forum (EMF) 37 study on deep decarbonization and high electrification analyzed a set of scenarios that achieve economy-wide net-zero carbon dioxide (CO2) emissions in North America by mid-century, exploring the implications of different technology evolutions, policies, and behavioral assumptions affecting energy supply and demand. For this paper, 16 modeling teams reported resulting emissions projections, energy system evolution, and economic activity. This paper provides an overview of the study, documents the scenario design, provides a roadmap for complementary forthcoming papers from this study, and offers an initial summary and comparison of results for net-zero CO2 by 2050 scenarios in the United States. We compare various outcomes across models and scenarios, such as emissions, energy use, fuel mix evolution, and technology adoption. Despite disparate model structure and sources for input assumptions, there is broad agreement in energy system trends across models towards deep decarbonization of the electricity sector coupled with increased end-use electrification of buildings, transportation, and to a lesser extent industry. All models deploy negative emissions technologies (e.g., direct air capture and bioenergy with carbon capture and storage) in addition to land sinks to achieve net-zero CO2 emissions. Important differences emerged in the results, showing divergent pathways among end-use sectors with deep electrification and grid decarbonization as necessary but not sufficient conditions to achieve net zero. These differences will be explored in the papers complementing this study to inform efforts to reach net-zero emissions and future research needs.

Emissions and energy impacts of the Inflation Reduction Act.

Economy-wide emissions drop 43 to 48% below 2005 levels by 2035 with accelerated clean energy deployment.

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.