Malignancy is accompanied by changes in the metabolism of both cells and the organism1,2. Pancreatic ductal adenocarcinoma (PDAC) is associated with wasting of peripheral tissues, a metabolic syndrome that lowers quality of life and has been proposed to decrease survival of patients with cancer3,4. Tissue wasting is a multifactorial disease and targeting specific circulating factors to reverse this syndrome has been mostly ineffective in the clinic5,6. Here we show that loss of both adipose and muscle tissue occurs early in the development of pancreatic cancer. Using mouse models of PDAC, we show that tumour growth in the pancreas but not in other sites leads to adipose tissue wasting, suggesting that tumour growth within the pancreatic environment contributes to this wasting phenotype. We find that decreased exocrine pancreatic function is a driver of adipose tissue loss and that replacement of pancreatic enzymes attenuates PDAC-associated wasting of peripheral tissues. Paradoxically, reversal of adipose tissue loss impairs survival in mice with PDAC. When analysing patients with PDAC, we find that depletion of adipose and skeletal muscle tissues at the time of diagnosis is common, but is not associated with worse survival. Taken together, these results provide an explanation for wasting of adipose tissue in early PDAC and suggest that early loss of peripheral tissue associated with pancreatic cancer may not impair survival.

Immune checkpoint blockade (ICB) results in durable disease control in a subset of patients with advanced renal cell carcinoma (RCC), but mechanisms driving resistance are poorly understood. We characterize the single-cell transcriptomes of cancer and immune cells from metastatic RCC patients before or after ICB exposure. In responders, subsets of cytotoxic T cells express higher levels of co-inhibitory receptors and effector molecules. Macrophages from treated biopsies shift toward pro-inflammatory states in response to an interferon-rich microenvironment but also upregulate immunosuppressive markers. In cancer cells, we identify bifurcation into two subpopulations differing in angiogenic signaling and upregulation of immunosuppressive programs after ICB. Expression signatures for cancer cell subpopulations and immune evasion are associated with PBRM1 mutation and survival in primary and ICB-treated advanced RCC. Our findings demonstrate that ICB remodels the RCC microenvironment and modifies the interplay between cancer and immune cell populations critical for understanding response and resistance to ICB.

Recombinant granulocyte colony-stimulating factor (G-CSF) is routinely administered for prophylaxis or treatment of chemotherapy-induced neutropenia. Chronic myelopoiesis and granulopoiesis in patients with cancer has been shown to induce immature monocytes and neutrophils that contribute to both systemic and local immunosuppression in the tumor microenvironment. The effect of recombinant G-CSF (pegfilgrastim or filgrastim) on the production of myeloid-derived suppressive cells is unknown. Here we examined patients with pancreatic cancer, a disease known to induce myeloid-derived suppressor cells (MDSCs), and for which pegfilgrastim is routinely administered concurrently with FOLFIRINOX but not with gemcitabine-based chemotherapy regimens.

Approximately 8% to 10% of pancreatic ductal adenocarcinomas (PDAC) do not harbor mutations in KRAS. Understanding the unique molecular and clinical features of this subset of pancreatic cancer is important to guide patient stratification for clinical trials of molecularly targeted agents.

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.

MicroRNAs (miRNAs) are post-transcriptional regulators of gene expression critical for organismal viability. Changes in miRNA activity are common in cancer, but how these changes relate to subsequent alterations in transcription and the process of tumorigenesis is not well understood. Here, we report a deep transcriptional, oncogenic network regulated by miRNAs. We present analysis of the gene expression and phenotypic changes associated with global miRNA restoration in miRNA-deficient fibroblasts. This analysis uncovers a miRNA-repressed network containing oncofetal genes Imp1, Imp2, and Imp3 (Imp1-3) that is up-regulated primarily transcriptionally >100-fold upon Dicer loss and is resistant to resilencing by complete restoration of miRNA activity. This Dicer-resistant epigenetic switch confers tumorigenicity to these cells. Let-7 targets Imp1-3 are required for this tumorigenicity and feed back to reinforce and sustain expression of the oncogenic network. Together, these Dicer-resistant genes constitute an mRNA expression signature that is present in numerous human cancers and is associated with poor survival.

Immune checkpoint inhibitors have revolutionized cancer treatment, but the benefits in refractory patients with esophageal cancer have been modest. Predictors of response as well as new targets for novel therapeutic combinations are needed. In this phase 2 clinical trial, we tested single-agent pembrolizumab in patients with advanced esophageal cancer, who received at least one prior line of therapy.

With the growing number of available targeted therapeutics and molecular biomarkers, the optimal care of patients with cancer now depends on a comprehensive understanding of the rapidly evolving landscape of precision oncology, which can be challenging for oncologists to navigate alone.

Transient activation of Src oncoprotein in non-transformed, breast epithelial cells can initiate an epigenetic switch to the stably transformed state via a positive feedback loop that involves the inflammatory transcription factors STAT3 and NF-κB. Here, we develop an experimental and computational pipeline that includes 1) a Bayesian network model (AccessTF) that accurately predicts protein-bound DNA sequence motifs based on chromatin accessibility, and 2) a scoring system (TFScore) that rank-orders transcription factors as candidates for being important for a biological process. Genetic experiments validate TFScore and suggest that more than 40 transcription factors contribute to the oncogenic state in this model. Interestingly, individual depletion of several of these factors results in similar transcriptional profiles, indicating that a complex and interconnected transcriptional network promotes a stable oncogenic state. The combined experimental and computational pipeline represents a general approach to comprehensively identify transcriptional regulators important for a biological process.

Immune responses to cancer are highly variable, with mismatch repair-deficient (MMRd) tumors exhibiting more anti-tumor immunity than mismatch repair-proficient (MMRp) tumors. To understand the rules governing these varied responses, we transcriptionally profiled 371,223 cells from colorectal tumors and adjacent normal tissues of 28 MMRp and 34 MMRd individuals. Analysis of 88 cell subsets and their 204 associated gene expression programs revealed extensive transcriptional and spatial remodeling across tumors. To discover hubs of interacting malignant and immune cells, we identified expression programs in different cell types that co-varied across tumors from affected individuals and used spatial profiling to localize coordinated programs. We discovered a myeloid cell-attracting hub at the tumor-luminal interface associated with tissue damage and an MMRd-enriched immune hub within the tumor, with activated T cells together with malignant and myeloid cells expressing T cell-attracting chemokines. By identifying interacting cellular programs, we reveal the logic underlying spatially organized immune-malignant cell networks.

Metastatic castration-resistant prostate cancer is typically lethal, exhibiting intrinsic or acquired resistance to second-generation androgen-targeting therapies and minimal response to immune checkpoint inhibitors1. Cellular programs driving resistance in both cancer and immune cells remain poorly understood. We present single-cell transcriptomes from 14 patients with advanced prostate cancer, spanning all common metastatic sites. Irrespective of treatment exposure, adenocarcinoma cells pervasively coexpressed multiple androgen receptor isoforms, including truncated isoforms hypothesized to mediate resistance to androgen-targeting therapies2,3. Resistance to enzalutamide was associated with cancer cell-intrinsic epithelial-mesenchymal transition and transforming growth factor-β signaling. Small cell carcinoma cells exhibited divergent expression programs driven by transcriptional regulators promoting lineage plasticity and HOXB5, HOXB6 and NR1D2 (refs. 4-6). Additionally, a subset of patients had high expression of dysfunction markers on cytotoxic CD8+ T cells undergoing clonal expansion following enzalutamide treatment. Collectively, the transcriptional characterization of cancer and immune cells from human metastatic castration-resistant prostate cancer provides a basis for the development of therapeutic approaches complementing androgen signaling inhibition.

Treatment of BRAF-mutant melanomas with MAP kinase pathway inhibitors is paradigmatic of the promise of precision cancer therapy but also highlights problems with drug resistance that limit patient benefit. We use live-cell imaging, single-cell analysis, and molecular profiling to show that exposure of tumor cells to RAF/MEK inhibitors elicits a heterogeneous response in which some cells die, some arrest, and the remainder adapt to drug. Drug-adapted cells up-regulate markers of the neural crest (e.g., NGFR), a melanocyte precursor, and grow slowly. This phenotype is transiently stable, reverting to the drug-naïve state within 9 days of drug withdrawal. Transcriptional profiling of cell lines and human tumors implicates a c-Jun/ECM/FAK/Src cascade in de-differentiation in about one-third of cell lines studied; drug-induced changes in c-Jun and NGFR levels are also observed in xenograft and human tumors. Drugs targeting the c-Jun/ECM/FAK/Src cascade as well as BET bromodomain inhibitors increase the maximum effect (Emax) of RAF/MEK kinase inhibitors by promoting cell killing. Thus, analysis of reversible drug resistance at a single-cell level identifies signaling pathways and inhibitory drugs missed by assays that focus on cell populations.

Neuroblastoma is a pediatric cancer arising from the developing sympathoadrenal lineage with complex inter- and intra-tumoral heterogeneity. To chart this complexity, we generated a comprehensive cell atlas of 55 neuroblastoma patient tumors, collected from two pediatric cancer institutions, spanning a range of clinical, genetic, and histologic features. Our atlas combines single-cell/nucleus RNA-seq (sc/scRNA-seq), bulk RNA-seq, whole exome sequencing, DNA methylation profiling, spatial transcriptomics, and two spatial proteomic methods. Sc/snRNA-seq revealed three malignant cell states with features of sympathoadrenal lineage development. All of the neuroblastomas had malignant cells that resembled sympathoblasts and the more differentiated adrenergic cells. A subset of tumors had malignant cells in a mesenchymal cell state with molecular features of Schwann cell precursors. DNA methylation profiles defined four groupings of patients, which differ in the degree of malignant cell heterogeneity and clinical outcomes. Using spatial proteomics, we found that neuroblastomas are spatially compartmentalized, with malignant tumor cells sequestered away from immune cells. Finally, we identify spatially restricted signaling patterns in immune cells from spatial transcriptomics. To facilitate the visualization and analysis of our atlas as a resource for further research in neuroblastoma, single cell, and spatial-omics, all data are shared through the Human Tumor Atlas Network Data Commons at www.humantumoratlas.org.

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.

Serological assays are important diagnostic tools for surveying exposure to the pathogen, monitoring immune response post vaccination, and managing spread of the infectious agent among the population. Current serological laboratory assays are often limited because they require the use of specialized laboratory technology and/or work with a limited number of sample types. Here, we evaluate an alternative by developing time-resolved Förster resonance energy transfer (TR-FRET) homogeneous assays that exhibited exceptional versatility, scalability, and sensitivity and outperformed or matched currently used strategies in terms of sensitivity, specificity, and precision. We validated the performance of the assays measuring total immunoglobulin G (IgG) levels; antibodies against severe acute respiratory syndrome coronavirus (SARS-CoV) or Middle Eastern respiratory syndrome (MERS)-CoV spike (S) protein; and SARS-CoV-2 S and nucleocapsid (N) proteins and applied it to several large sample sets and real-world applications. We further established a TR-FRET-based ACE2-S competition assay to assess the neutralization propensity of the antibodies. Overall, these TR-FRET-based serological assays can be rapidly extended to other antigens and are compatible with commonly used plate readers.

Single-cell genomics is essential to chart tumor ecosystems. Although single-cell RNA-Seq (scRNA-Seq) profiles RNA from cells dissociated from fresh tumors, single-nucleus RNA-Seq (snRNA-Seq) is needed to profile frozen or hard-to-dissociate tumors. Each requires customization to different tissue and tumor types, posing a barrier to adoption. Here, we have developed a systematic toolbox for profiling fresh and frozen clinical tumor samples using scRNA-Seq and snRNA-Seq, respectively. We analyzed 216,490 cells and nuclei from 40 samples across 23 specimens spanning eight tumor types of varying tissue and sample characteristics. We evaluated protocols by cell and nucleus quality, recovery rate and cellular composition. scRNA-Seq and snRNA-Seq from matched samples recovered the same cell types, but at different proportions. Our work provides guidance for studies in a broad range of tumors, including criteria for testing and selecting methods from the toolbox for other tumors, thus paving the way for charting tumor atlases.