Lineage tracing reveals the phylodynamics, plasticity, and paths of tumor evolution.

Tumor evolution is driven by the progressive acquisition of genetic and epigenetic alterations that enable uncontrolled growth and expansion to neighboring and distal tissues. The study of phylogenetic relationships between cancer cells provides key insights into these processes. Here, we introduced an evolving lineage-tracing system with a single-cell RNA-seq readout into a mouse model of Kras;Trp53(KP)-driven lung adenocarcinoma and tracked tumor evolution from single-transformed cells to metastatic tumors at unprecedented resolution. We found that the loss of the initial, stable alveolar-type2-like state was accompanied by a transient increase in plasticity. This was followed by the adoption of distinct transcriptional programs that enable rapid expansion and, ultimately, clonal sweep of stable subclones capable of metastasizing. Finally, tumors develop through stereotypical evolutionary trajectories, and perturbing additional tumor suppressors accelerates progression by creating novel trajectories. Our study elucidates the hierarchical nature of tumor evolution and, more broadly, enables in-depth studies of tumor progression.

Antigen dominance hierarchies shape TCF1+ progenitor CD8 T cell phenotypes in tumors.

CD8 T cell responses against different tumor neoantigens occur simultaneously, yet little is known about the interplay between responses and its impact on T cell function and tumor control. In mouse lung adenocarcinoma, we found that immunodominance is established in tumors, wherein CD8 T cell expansion is predominantly driven by the antigen that most stably binds MHC. T cells responding to subdominant antigens were enriched for a TCF1+ progenitor phenotype correlated with response to immune checkpoint blockade (ICB) therapy. However, the subdominant T cell response did not preferentially benefit from ICB due to a dysfunctional subset of TCF1+ cells marked by CCR6 and Tc17 differentiation. Analysis of human samples and sequencing datasets revealed that CCR6+ TCF1+ cells exist across human cancers and are not correlated with ICB response. Vaccination eliminated CCR6+ TCF1+ cells and dramatically improved the subdominant response, highlighting a strategy to optimally engage concurrent neoantigen responses against tumors.

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.

Endocrine-Exocrine Signaling Drives Obesity-Associated Pancreatic Ductal Adenocarcinoma.

Obesity is a major modifiable risk factor for pancreatic ductal adenocarcinoma (PDAC), yet how and when obesity contributes to PDAC progression is not well understood. Leveraging an autochthonous mouse model, we demonstrate a causal and reversible role for obesity in early PDAC progression, showing that obesity markedly enhances tumorigenesis, while genetic or dietary induction of weight loss intercepts cancer development. Molecular analyses of human and murine samples define microenvironmental consequences of obesity that foster tumorigenesis rather than new driver gene mutations, including significant pancreatic islet cell adaptation in obesity-associated tumors. Specifically, we identify aberrant beta cell expression of the peptide hormone cholecystokinin (Cck) in response to obesity and show that islet Cck promotes oncogenic Kras-driven pancreatic ductal tumorigenesis. Our studies argue that PDAC progression is driven by local obesity-associated changes in the tumor microenvironment and implicate endocrine-exocrine signaling beyond insulin in PDAC development.

Commensal Microbiota Promote Lung Cancer Development via γδ T Cells.

Lung cancer is closely associated with chronic inflammation, but the causes of inflammation and the specific immune mediators have not been fully elucidated. The lung is a mucosal tissue colonized by a diverse bacterial community, and pulmonary infections commonly present in lung cancer patients are linked to clinical outcomes. Here, we provide evidence that local microbiota provoke inflammation associated with lung adenocarcinoma by activating lung-resident γδ T cells. Germ-free or antibiotic-treated mice were significantly protected from lung cancer development induced by Kras mutation and p53 loss. Mechanistically, commensal bacteria stimulated Myd88-dependent IL-1ß and IL-23 production from myeloid cells, inducing proliferation and activation of Vγ6+Vδ1+ γδ T cells that produced IL-17 and other effector molecules to promote inflammation and tumor cell proliferation. Our findings clearly link local microbiota-immune crosstalk to lung tumor development and thereby define key cellular and molecular mediators that may serve as effective targets in lung cancer intervention.

Anatomically and Functionally Distinct Lung Mesenchymal Populations Marked by Lgr5 and Lgr6.

The diversity of mesenchymal cell types in the lung that influence epithelial homeostasis and regeneration is poorly defined. We used genetic lineage tracing, single-cell RNA sequencing, and organoid culture approaches to show that Lgr5 and Lgr6, well-known markers of stem cells in epithelial tissues, are markers of mesenchymal cells in the adult lung. Lgr6+ cells comprise a subpopulation of smooth muscle cells surrounding airway epithelia and promote airway differentiation of epithelial progenitors via Wnt-Fgf10 cooperation. Genetic ablation of Lgr6+ cells impairs airway injury repair in vivo. Distinct Lgr5+ cells are located in alveolar compartments and are sufficient to promote alveolar differentiation of epithelial progenitors through Wnt activation. Modulating Wnt activity altered differentiation outcomes specified by mesenchymal cells. This identification of region- and lineage-specific crosstalk between epithelium and their neighboring mesenchymal partners provides new understanding of how different cell types are maintained in the adult lung.

Lung Adenocarcinoma Distally Rewires Hepatic Circadian Homeostasis.

The circadian clock controls metabolic and physiological processes through finely tuned molecular mechanisms. The clock is remarkably plastic and adapts to exogenous "zeitgebers," such as light and nutrition. How a pathological condition in a given tissue influences systemic circadian homeostasis in other tissues remains an unanswered question of conceptual and biomedical importance. Here, we show that lung adenocarcinoma operates as an endogenous reorganizer of circadian metabolism. High-throughput transcriptomics and metabolomics revealed unique signatures of transcripts and metabolites cycling exclusively in livers of tumor-bearing mice. Remarkably, lung cancer has no effect on the core clock but rather reprograms hepatic metabolism through altered pro-inflammatory response via the STAT3-Socs3 pathway. This results in disruption of AKT, AMPK, and SREBP signaling, leading to altered insulin, glucose, and lipid metabolism. Thus, lung adenocarcinoma functions as a potent endogenous circadian organizer (ECO), which rewires the pathophysiological dimension of a distal tissue such as the liver. PAPERCLIP.

Conventional type I dendritic cells maintain a reservoir of proliferative tumor-antigen specific TCF-1+ CD8+ T cells in tumor-draining lymph nodes.

In tumors, a subset of CD8+ T cells expressing the transcription factor TCF-1 drives the response to immune checkpoint blockade. We examined the mechanisms that maintain these cells in an autochthonous model of lung adenocarcinoma. Longitudinal sampling and single-cell sequencing of tumor-antigen specific TCF-1+ CD8+ T cells revealed that while intratumoral TCF-1+ CD8+ T cells acquired dysfunctional features and decreased in number as tumors progressed, TCF-1+ CD8+ T cell frequency in the tumor draining LN (dLN) remained stable. Two discrete intratumoral TCF-1+ CD8+ T cell subsets developed over time-a proliferative SlamF6+ subset and a non-cycling SlamF6- subset. Blocking dLN egress decreased the frequency of intratumoral SlamF6+ TCF-1+ CD8+ T cells. Conventional type I dendritic cell (cDC1) in dLN decreased in number with tumor progression, and Flt3L+anti-CD40 treatment recovered SlamF6+ T cell frequencies and decreased tumor burden. Thus, cDC1s in tumor dLN maintain a reservoir of TCF-1+ CD8+ T cells and their decrease contributes to failed anti-tumor immunity.

SOX17 enables immune evasion of early colorectal adenomas and cancers.

A hallmark of cancer is the avoidance of immune destruction. This process has been primarily investigated in locally advanced or metastatic cancer1-3; however, much less is known about how pre-malignant or early invasive tumours evade immune detection. Here, to understand this process in early colorectal cancers (CRCs), we investigated how naive colon cancer organoids that were engineered in vitro to harbour Apc-null, KrasG12D and Trp53-null (AKP) mutations adapted to the in vivo native colonic environment. Comprehensive transcriptomic and chromatin analyses revealed that the endoderm-specifying transcription factor SOX17 became strongly upregulated in vivo. Notably, whereas SOX17 loss did not affect AKP organoid propagation in vitro, its loss markedly reduced the ability of AKP tumours to persist in vivo. The small fraction of SOX17-null tumours that grew displayed notable interferon-γ (IFNγ)-producing effector-like CD8+ T cell infiltrates in contrast to the immune-suppressive microenvironment in wild-type counterparts. Mechanistically, in both endogenous Apc-null pre-malignant adenomas and transplanted organoid-derived AKP CRCs, SOX17 suppresses the ability of tumour cells to sense and respond to IFNγ, preventing anti-tumour T cell responses. Finally, SOX17 engages a fetal intestinal programme that drives differentiation away from LGR5+ tumour cells to produce immune-evasive LGR5- tumour cells with lower expression of major histocompatibility complex class I (MHC-I). We propose that SOX17 is a transcription factor that is engaged during the early steps of colon cancer to orchestrate an immune-evasive programme that permits CRC initiation and progression.

Deciphering cell states and genealogies of human haematopoiesis.

The human blood system is maintained through the differentiation and massive amplification of a limited number of long-lived haematopoietic stem cells (HSCs)1. Perturbations to this process underlie diverse diseases, but the clonal contributions to human haematopoiesis and how this changes with age remain incompletely understood. Although recent insights have emerged from barcoding studies in model systems2-5, simultaneous detection of cell states and phylogenies from natural barcodes in humans remains challenging. Here we introduce an improved, single-cell lineage-tracing system based on deep detection of naturally occurring mitochondrial DNA mutations with simultaneous readout of transcriptional states and chromatin accessibility. We use this system to define the clonal architecture of HSCs and map the physiological state and output of clones. We uncover functional heterogeneity in HSC clones, which is stable over months and manifests as both differences in total HSC output and biases towards the production of different mature cell types. We also find that the diversity of HSC clones decreases markedly with age, leading to an oligoclonal structure with multiple distinct clonal expansions. Our study thus provides a clonally resolved and cell-state-aware atlas of human haematopoiesis at single-cell resolution, showing an unappreciated functional diversity of human HSC clones and, more broadly, paving the way for refined studies of clonal dynamics across a range of tissues in human health and disease.

Modeling diverse genetic subtypes of lung adenocarcinoma with a next-generation alveolar type 2 organoid platform.

Lung cancer is the leading cause of cancer-related death worldwide. Lung adenocarcinoma (LUAD), the most common histological subtype, accounts for 40% of all cases. While existing genetically engineered mouse models (GEMMs) recapitulate the histological progression and transcriptional evolution of human LUAD, they are time-consuming and technically demanding. In contrast, cell line transplant models are fast and flexible, but these models fail to capture the full spectrum of disease progression. Organoid technologies provide a means to create next-generation cancer models that integrate the most advantageous features of autochthonous and transplant-based systems. However, robust and faithful LUAD organoid platforms are currently lacking. Here, we describe optimized conditions to continuously expand murine alveolar type 2 (AT2) cells, a prominent cell of origin for LUAD, in organoid culture. These organoids display canonical features of AT2 cells, including marker gene expression, the presence of lamellar bodies, and an ability to differentiate into the AT1 lineage. We used this system to develop flexible and versatile immunocompetent organoid-based models of KRAS, BRAF, and ALK mutant LUAD. Notably, organoid-based tumors display extensive burden and complete penetrance and are histopathologically indistinguishable from their autochthonous counterparts. Altogether, this organoid platform is a powerful, versatile new model system to study LUAD.

KRAS and YAP1 converge to regulate EMT and tumor survival.

Cancer cells that express oncogenic alleles of RAS typically require sustained expression of the mutant allele for survival, but the molecular basis of this oncogene dependency remains incompletely understood. To identify genes that can functionally substitute for oncogenic RAS, we systematically expressed 15,294 open reading frames in a human KRAS-dependent colon cancer cell line engineered to express an inducible KRAS-specific shRNA. We found 147 genes that promoted survival upon KRAS suppression. In particular, the transcriptional coactivator YAP1 rescued cell viability in KRAS-dependent cells upon suppression of KRAS and was required for KRAS-induced cell transformation. Acquired resistance to Kras suppression in a Kras-driven murine lung cancer model also involved increased YAP1 signaling. KRAS and YAP1 converge on the transcription factor FOS and activate a transcriptional program involved in regulating the epithelial-mesenchymal transition (EMT). Together, these findings implicate transcriptional regulation of EMT by YAP1 as a significant component of oncogenic RAS signaling.

Genetic and clonal dissection of murine small cell lung carcinoma progression by genome sequencing.

Small cell lung carcinoma (SCLC) is a highly lethal, smoking-associated cancer with few known targetable genetic alterations. Using genome sequencing, we characterized the somatic evolution of a genetically engineered mouse model (GEMM) of SCLC initiated by loss of Trp53 and Rb1. We identified alterations in DNA copy number and complex genomic rearrangements and demonstrated a low somatic point mutation frequency in the absence of tobacco mutagens. Alterations targeting the tumor suppressor Pten occurred in the majority of murine SCLC studied, and engineered Pten deletion accelerated murine SCLC and abrogated loss of Chr19 in Trp53; Rb1; Pten compound mutant tumors. Finally, we found evidence for polyclonal and sequential metastatic spread of murine SCLC by comparative sequencing of families of related primary tumors and metastases. We propose a temporal model of SCLC tumorigenesis with implications for human SCLC therapeutics and the nature of cancer-genome evolution in GEMMs.

PKM2 isoform-specific deletion reveals a differential requirement for pyruvate kinase in tumor cells.

The pyruvate kinase M2 isoform (PKM2) is expressed in cancer and plays a role in regulating anabolic metabolism. To determine whether PKM2 is required for tumor formation or growth, we generated mice with a conditional allele that abolishes PKM2 expression without disrupting PKM1 expression. PKM2 deletion accelerated mammary tumor formation in a Brca1-loss-driven model of breast cancer. PKM2 null tumors displayed heterogeneous PKM1 expression, with PKM1 found in nonproliferating tumor cells and no detectable pyruvate kinase expression in proliferating cells. This suggests that PKM2 is not necessary for tumor cell proliferation and implies that the inactive state of PKM2 is associated with the proliferating cell population within tumors, whereas nonproliferating tumor cells require active pyruvate kinase. Consistent with these findings, variable PKM2 expression and heterozygous PKM2 mutations are found in human tumors. These data suggest that regulation of PKM2 activity supports the different metabolic requirements of proliferating and nonproliferating tumor cells.

Spatial genomics enables multi-modal study of clonal heterogeneity in tissues.

The state and behaviour of a cell can be influenced by both genetic and environmental factors. In particular, tumour progression is determined by underlying genetic aberrations1-4 as well as the makeup of the tumour microenvironment5,6. Quantifying the contributions of these factors requires new technologies that can accurately measure the spatial location of genomic sequence together with phenotypic readouts. Here we developed slide-DNA-seq, a method for capturing spatially resolved DNA sequences from intact tissue sections. We demonstrate that this method accurately preserves local tumour architecture and enables the de novo discovery of distinct tumour clones and their copy number alterations. We then apply slide-DNA-seq to a mouse model of metastasis and a primary human cancer, revealing that clonal populations are confined to distinct spatial regions. Moreover, through integration with spatial transcriptomics, we uncover distinct sets of genes that are associated with clone-specific genetic aberrations, the local tumour microenvironment, or both. Together, this multi-modal spatial genomics approach provides a versatile platform for quantifying how cell-intrinsic and cell-extrinsic factors contribute to gene expression, protein abundance and other cellular phenotypes.

Deciphering the immunopeptidome in vivo reveals new tumour antigens.

Immunosurveillance of cancer requires the presentation of peptide antigens on major histocompatibility complex class I (MHC-I) molecules1-5. Current approaches to profiling of MHC-I-associated peptides, collectively known as the immunopeptidome, are limited to in vitro investigation or bulk tumour lysates, which limits our understanding of cancer-specific patterns of antigen presentation in vivo6. To overcome these limitations, we engineered an inducible affinity tag into the mouse MHC-I gene (H2-K1) and targeted this allele to the KrasLSL-G12D/+Trp53fl/fl mouse model (KP/KbStrep)7. This approach enabled us to precisely isolate MHC-I peptides from autochthonous pancreatic ductal adenocarcinoma and from lung adenocarcinoma (LUAD) in vivo. In addition, we profiled the LUAD immunopeptidome from the alveolar type 2 cell of origin up to late-stage disease. Differential peptide presentation in LUAD was not predictable by mRNA expression or translation efficiency and is probably driven by post-translational mechanisms. Vaccination with peptides presented by LUAD in vivo induced CD8+ T cell responses in naive mice and tumour-bearing mice. Many peptides specific to LUAD, including immunogenic peptides, exhibited minimal expression of the cognate mRNA, which prompts the reconsideration of antigen prediction pipelines that triage peptides according to transcript abundance8. Beyond cancer, the KbStrep allele is compatible with other Cre-driver lines to explore antigen presentation in vivo in the pursuit of understanding basic immunology, infectious disease and autoimmunity.

Optimality in the development of intestinal crypts.

Intestinal crypts in mammals are comprised of long-lived stem cells and shorter-lived progenies. These two populations are maintained in specific proportions during adult life. Here, we investigate the design principles governing the dynamics of these proportions during crypt morphogenesis. Using optimal control theory, we show that a proliferation strategy known as a "bang-bang" control minimizes the time to obtain a mature crypt. This strategy consists of a surge of symmetric stem cell divisions, establishing the entire stem cell pool first, followed by a sharp transition to strictly asymmetric stem cell divisions, producing nonstem cells with a delay. We validate these predictions using lineage tracing and single-molecule fluorescence in situ hybridization of intestinal crypts in infant mice, uncovering small crypts that are entirely composed of Lgr5-labeled stem cells, which become a minority as crypts continue to grow. Our approach can be used to uncover similar design principles in other developmental systems.

Driving Rel-iant Tregs toward an Identity Crisis.

Inhibiting Treg cell function in tumors is an attractive strategy to improve anti-cancer immunity. In a pair of papers in Immunity and Cell, Ghosh and colleagues show that the canonical NF-κB subunits p65 and c-Rel have non-redundant, critical roles in promoting Treg cell development and function (Oh et al., 2017). Targeting c-Rel blunts Treg cell immunosuppressive activity in the tumor microenvironment and enhances anti-tumor T cell responses (Grinberg-Bleyer et al., 2017).

A large intergenic noncoding RNA induced by p53 mediates global gene repression in the p53 response.

Recently, more than 1000 large intergenic noncoding RNAs (lincRNAs) have been reported. These RNAs are evolutionarily conserved in mammalian genomes and thus presumably function in diverse biological processes. Here, we report the identification of lincRNAs that are regulated by p53. One of these lincRNAs (lincRNA-p21) serves as a repressor in p53-dependent transcriptional responses. Inhibition of lincRNA-p21 affects the expression of hundreds of gene targets enriched for genes normally repressed by p53. The observed transcriptional repression by lincRNA-p21 is mediated through the physical association with hnRNP-K. This interaction is required for proper genomic localization of hnRNP-K at repressed genes and regulation of p53 mediates apoptosis. We propose a model whereby transcription factors activate lincRNAs that serve as key repressors by physically associating with repressive complexes and modulate their localization to sets of previously active genes.