Human cancer cell lines are the workhorse of cancer research. Although cell lines are known to evolve in culture, the extent of the resultant genetic and transcriptional heterogeneity and its functional consequences remain understudied. Here we use genomic analyses of 106 human cell lines grown in two laboratories to show extensive clonal diversity. Further comprehensive genomic characterization of 27 strains of the common breast cancer cell line MCF7 uncovered rapid genetic diversification. Similar results were obtained with multiple strains of 13 additional cell lines. Notably, genetic changes were associated with differential activation of gene expression programs and marked differences in cell morphology and proliferation. Barcoding experiments showed that cell line evolution occurs as a result of positive clonal selection that is highly sensitive to culture conditions. Analyses of single-cell-derived clones demonstrated that continuous instability quickly translates into heterogeneity of the cell line. When the 27 MCF7 strains were tested against 321 anti-cancer compounds, we uncovered considerably different drug responses: at least 75% of compounds that strongly inhibited some strains were completely inactive in others. This study documents the extent, origins and consequences of genetic variation within cell lines, and provides a framework for researchers to measure such variation in efforts to support maximally reproducible cancer research.

The Library of Integrated Network-Based Cellular Signatures (LINCS) is an NIH Common Fund program that catalogs how human cells globally respond to chemical, genetic, and disease perturbations. Resources generated by LINCS include experimental and computational methods, visualization tools, molecular and imaging data, and signatures. By assembling an integrated picture of the range of responses of human cells exposed to many perturbations, the LINCS program aims to better understand human disease and to advance the development of new therapies. Perturbations under study include drugs, genetic perturbations, tissue micro-environments, antibodies, and disease-causing mutations. Responses to perturbations are measured by transcript profiling, mass spectrometry, cell imaging, and biochemical methods, among other assays. The LINCS program focuses on cellular physiology shared among tissues and cell types relevant to an array of diseases, including cancer, heart disease, and neurodegenerative disorders. This Perspective describes LINCS technologies, datasets, tools, and approaches to data accessibility and reusability.

Diffuse large B cell lymphoma (DLBCL), the most common lymphoid malignancy in adults, is a clinically and genetically heterogeneous disease that is further classified into transcriptionally defined activated B cell (ABC) and germinal center B cell (GCB) subtypes. We carried out a comprehensive genetic analysis of 304 primary DLBCLs and identified low-frequency alterations, captured recurrent mutations, somatic copy number alterations, and structural variants, and defined coordinate signatures in patients with available outcome data. We integrated these genetic drivers using consensus clustering and identified five robust DLBCL subsets, including a previously unrecognized group of low-risk ABC-DLBCLs of extrafollicular/marginal zone origin; two distinct subsets of GCB-DLBCLs with different outcomes and targetable alterations; and an ABC/GCB-independent group with biallelic inactivation of TP53, CDKN2A loss, and associated genomic instability. The genetic features of the newly characterized subsets, their mutational signatures, and the temporal ordering of identified alterations provide new insights into DLBCL pathogenesis. The coordinate genetic signatures also predict outcome independent of the clinical International Prognostic Index and suggest new combination treatment strategies. More broadly, our results provide a roadmap for an actionable DLBCL classification.

In the version of this article originally published, some text above the "Tri-nucleotide sequence motifs" label in Fig. 2a appeared incorrectly. The text was garbled and should have appeared as nucleotide codes.Additionally, the labels on the bars in Fig. 2c were not italicized in the original publication. These are gene symbols, and they should have been italicized.The colored labels above the graphs in Fig. 4b were also erroneously not italicized. These labels represent gene names and loci, and they should have been italicized.

Pharmacologically difficult targets, such as MYC transcription factors, represent a major challenge in cancer therapy. For the childhood cancer neuroblastoma, amplification of the oncogene MYCN is associated with high-risk disease and poor prognosis. Here, we deployed genome-scale CRISPR-Cas9 screening of MYCN-amplified neuroblastoma and found a preferential dependency on genes encoding the polycomb repressive complex 2 (PRC2) components EZH2, EED, and SUZ12. Genetic and pharmacological suppression of EZH2 inhibited neuroblastoma growth in vitro and in vivo. Moreover, compared with neuroblastomas without MYCN amplification, MYCN-amplified neuroblastomas expressed higher levels of EZH2. ChIP analysis showed that MYCN binds at the EZH2 promoter, thereby directly driving expression. Transcriptomic and epigenetic analysis, as well as genetic rescue experiments, revealed that EZH2 represses neuronal differentiation in neuroblastoma in a PRC2-dependent manner. Moreover, MYCN-amplified and high-risk primary tumors from patients with neuroblastoma exhibited strong repression of EZH2-regulated genes. Additionally, overexpression of IGFBP3, a direct EZH2 target, suppressed neuroblastoma growth in vitro and in vivo. We further observed strong synergy between histone deacetylase inhibitors and EZH2 inhibitors. Together, these observations demonstrate that MYCN upregulates EZH2, leading to inactivation of a tumor suppressor program in neuroblastoma, and support testing EZH2 inhibitors in patients with MYCN-amplified neuroblastoma.

Glioblastoma (GBM) is the most common primary adult brain tumor. Despite extensive efforts, the median survival for GBM patients is approximately 14 months. GBM therapy could benefit greatly from patient-specific targeted therapies that maximize treatment efficacy. Here we report a platform termed SynergySeq to identify drug combinations for the treatment of GBM by integrating information from The Cancer Genome Atlas (TCGA) and the Library of Integrated Network-Based Cellular Signatures (LINCS). We identify differentially expressed genes in GBM samples and devise a consensus gene expression signature for each compound using LINCS L1000 transcriptional profiling data. The SynergySeq platform computes disease discordance and drug concordance to identify combinations of FDA-approved drugs that induce a synergistic response in GBM. Collectively, our studies demonstrate that combining disease-specific gene expression signatures with LINCS small molecule perturbagen-response signatures can identify preclinical combinations for GBM, which can potentially be tested in humans.

Drug resistance to approved systemic therapies in estrogen receptor-positive (ER+) breast cancer remains common. We hypothesized that factors present in the human tumor microenvironment (TME) drive drug resistance. Screening of a library of recombinant secreted microenvironmental proteins revealed fibroblast growth factor 2 (FGF2) as a potent mediator of resistance to anti-estrogens, mTORC1 inhibition, and phosphatidylinositol 3-kinase inhibition in ER+ breast cancer. Phosphoproteomic analyses identified ERK1/2 as a major output of FGF2 signaling via FGF receptors (FGFRs), with consequent up-regulation of Cyclin D1 and down-regulation of Bim as mediators of drug resistance. FGF2-driven drug resistance in anti-estrogen-sensitive and -resistant models, including patient-derived xenografts, was reverted by neutralizing FGF2 or FGFRs. A transcriptomic signature of FGF2 signaling in primary tumors predicted shorter recurrence-free survival independently of age, grade, stage, and FGFR amplification status. These findings delineate FGF2 signaling as a ligand-based drug resistance mechanism and highlights an underdeveloped aspect of precision oncology: characterizing and treating patients according to their TME constitution.

This Letter is being retracted owing to issues with Fig. 1d and Supplementary Fig. 31b, and the unavailability of original data for these figures that raise concerns regarding the integrity of the figures. Nature published two previous corrections related to this Letter(1,2). These issues in aggregate undermine the confidence in the integrity of this study. Authors Michael Foley, Monica Schenone, Nicola J. Tolliday, Todd R. Golub, Steven A. Carr, Alykhan F. Shamji, Andrew M. Stern and Stuart L. Schreiber agree with the Retraction. Authors Lakshmi Raj, Takao Ide, Aditi U. Gurkar, Anna Mandinova and Sam W. Lee disagree with the Retraction. Author Xiaoyu Li did not respond.

Understanding complex biological systems requires the system-wide characterization of both molecular and cellular features. Existing methods for spatial mapping of biomolecules in intact tissues suffer from information loss caused by degradation and tissue damage. We report a tissue transformation strategy named stabilization under harsh conditions via intramolecular epoxide linkages to prevent degradation (SHIELD), which uses a flexible polyepoxide to form controlled intra- and intermolecular cross-link with biomolecules. SHIELD preserves protein fluorescence and antigenicity, transcripts and tissue architecture under a wide range of harsh conditions. We applied SHIELD to interrogate system-level wiring, synaptic architecture, and molecular features of virally labeled neurons and their targets in mouse at single-cell resolution. We also demonstrated rapid three-dimensional phenotyping of core needle biopsies and human brain cells. SHIELD enables rapid, multiscale, integrated molecular phenotyping of both animal and clinical tissues.

In the version of this article originally published, an asterisk was omitted from Fig. 1a. The asterisk has been added to the figure. Additionally, a "NOTCH2" label was erroneously included in Fig. 4a. The label has been removed. The errors have been corrected in the PDF and HTML versions of this article.