LKB1 drives stasis and C/EBP-mediated reprogramming to an alveolar type II fate in lung cancer.

LKB1 is among the most frequently altered tumor suppressors in lung adenocarcinoma. Inactivation of Lkb1 accelerates the growth and progression of oncogenic KRAS-driven lung tumors in mouse models. However, the molecular mechanisms by which LKB1 constrains lung tumorigenesis and whether the cancer state that stems from Lkb1 deficiency can be reverted remains unknown. To identify the processes governed by LKB1 in vivo, we generated an allele which enables Lkb1 inactivation at tumor initiation and subsequent Lkb1 restoration in established tumors. Restoration of Lkb1 in oncogenic KRAS-driven lung tumors suppressed proliferation and led to tumor stasis. Lkb1 restoration activated targets of C/EBP transcription factors and drove neoplastic cells from a progenitor-like state to a less proliferative alveolar type II cell-like state. We show that C/EBP transcription factors govern a subset of genes that are induced by LKB1 and depend upon NKX2-1. We also demonstrate that a defining factor of the alveolar type II lineage, C/EBPα, constrains oncogenic KRAS-driven lung tumor growth in vivo. Thus, this key tumor suppressor regulates lineage-specific transcription factors, thereby constraining lung tumor development through enforced differentiation.

LKB1 inactivation modulates chromatin accessibility to drive metastatic progression.

Metastasis is the leading cause of cancer-related deaths and enables cancer cells to compromise organ function by expanding in secondary sites. Since primary tumours and metastases often share the same constellation of driver mutations, the mechanisms that drive their distinct phenotypes are unclear. Here we show that inactivation of the frequently mutated tumour suppressor gene LKB1 (encoding liver kinase B1) has evolving effects throughout the progression of lung cancer, which leads to the differential epigenetic re-programming of early-stage primary tumours compared with late-stage metastases. By integrating genome-scale CRISPR-Cas9 screening with bulk and single-cell multi-omic analyses, we unexpectedly identify LKB1 as a master regulator of chromatin accessibility in lung adenocarcinoma primary tumours. Using an in vivo model of metastatic progression, we further show that loss of LKB1 activates the early endoderm transcription factor SOX17 in metastases and a metastatic-like sub-population of cancer cells within primary tumours. The expression of SOX17 is necessary and sufficient to drive a second wave of epigenetic changes in LKB1-deficient cells that enhances metastatic ability. Overall, our study demonstrates how the downstream effects of an individual driver mutation can change throughout cancer development, with implications for stage-specific therapeutic resistance mechanisms and the gene regulatory underpinnings of metastatic evolution.

High-throughput single-cell chromatin accessibility CRISPR screens enable unbiased identification of regulatory networks in cancer.

Chromatin accessibility profiling can identify putative regulatory regions genome wide; however, pooled single-cell methods for assessing the effects of regulatory perturbations on accessibility are limited. Here, we report a modified droplet-based single-cell ATAC-seq protocol for perturbing and evaluating dynamic single-cell epigenetic states. This method (Spear-ATAC) enables simultaneous read-out of chromatin accessibility profiles and integrated sgRNA spacer sequences from thousands of individual cells at once. Spear-ATAC profiling of 104,592 cells representing 414 sgRNA knock-down populations reveals the temporal dynamics of epigenetic responses to regulatory perturbations in cancer cells and the associations between transcription factor binding profiles.

Finding needles in a haystack: dissecting tumor heterogeneity with single-cell transcriptomic and chromatin accessibility profiling.

Tumor evolution often results in a wealth of heterogeneous cancer cell types within a single tumor - heterogeneity that can include epigenetic and gene expression changes that are impossible to identify from histological features alone. The invasion of cancer cells into nearby healthy tissue, accompanied by the infiltration of responding immune cells, results in an even more complex architecture of tumor and non-tumor cells. However, bulk genomics-based methods can only assay the aggregate transcriptomic and epigenetic profiles across all of this rich cellular diversity. Such bulk averaging hides small subpopulations of tumor cells with unique phenotypes that might result in therapeutic resistance or metastatic progression. The advent of single-cell-based genomics assays for measuring transcription and chromatin accessibility - particularly scRNA-seq and scATAC-seq - has enabled the dissection of cell-types within tumors at a scale and resolution capable of unraveling the epigenetic and gene expression programs of rare and unique cellular subpopulations. This Review focuses on recent advances in scRNA-seq and scATAC-seq technologies and their application to cancer biology in the context of furthering our understanding of tumor heterogeneity.

Zmat3 Is a Key Splicing Regulator in the p53 Tumor Suppression Program.

Although TP53 is the most commonly mutated gene in human cancers, the p53-dependent transcriptional programs mediating tumor suppression remain incompletely understood. Here, to uncover critical components downstream of p53 in tumor suppression, we perform unbiased RNAi and CRISPR-Cas9-based genetic screens in vivo. These screens converge upon the p53-inducible gene Zmat3, encoding an RNA-binding protein, and we demonstrate that ZMAT3 is an important tumor suppressor downstream of p53 in mouse KrasG12D-driven lung and liver cancers and human carcinomas. Integrative analysis of the ZMAT3 RNA-binding landscape and transcriptomic profiling reveals that ZMAT3 directly modulates exon inclusion in transcripts encoding proteins of diverse functions, including the p53 inhibitors MDM4 and MDM2, splicing regulators, and components of varied cellular processes. Interestingly, these exons are enriched in NMD signals, and, accordingly, ZMAT3 broadly affects target transcript stability. Collectively, these studies reveal ZMAT3 as a novel RNA-splicing and homeostasis regulator and a key component of p53-mediated tumor suppression.

CRISPR screens in cancer spheroids identify 3D growth-specific vulnerabilities.

Cancer genomics studies have identified thousands of putative cancer driver genes1. Development of high-throughput and accurate models to define the functions of these genes is a major challenge. Here we devised a scalable cancer-spheroid model and performed genome-wide CRISPR screens in 2D monolayers and 3D lung-cancer spheroids. CRISPR phenotypes in 3D more accurately recapitulated those of in vivo tumours, and genes with differential sensitivities between 2D and 3D conditions were highly enriched for genes that are mutated in lung cancers. These analyses also revealed drivers that are essential for cancer growth in 3D and in vivo, but not in 2D. Notably, we found that carboxypeptidase D is responsible for removal of a C-terminal RKRR motif2 from the α-chain of the insulin-like growth factor 1 receptor that is critical for receptor activity. Carboxypeptidase D expression correlates with patient outcomes in patients with lung cancer, and loss of carboxypeptidase D reduced tumour growth. Our results reveal key differences between 2D and 3D cancer models, and establish a generalizable strategy for performing CRISPR screens in spheroids to reveal cancer vulnerabilities.

Massively parallel single-cell chromatin landscapes of human immune cell development and intratumoral T cell exhaustion.

Understanding complex tissues requires single-cell deconstruction of gene regulation with precision and scale. Here, we assess the performance of a massively parallel droplet-based method for mapping transposase-accessible chromatin in single cells using sequencing (scATAC-seq). We apply scATAC-seq to obtain chromatin profiles of more than 200,000 single cells in human blood and basal cell carcinoma. In blood, application of scATAC-seq enables marker-free identification of cell type-specific cis- and trans-regulatory elements, mapping of disease-associated enhancer activity and reconstruction of trajectories of cellular differentiation. In basal cell carcinoma, application of scATAC-seq reveals regulatory networks in malignant, stromal and immune cells in the tumor microenvironment. Analysis of scATAC-seq profiles from serial tumor biopsies before and after programmed cell death protein 1 blockade identifies chromatin regulators of therapy-responsive T cell subsets and reveals a shared regulatory program that governs intratumoral CD8+ T cell exhaustion and CD4+ T follicular helper cell development. We anticipate that scATAC-seq will enable the unbiased discovery of gene regulatory factors across diverse biological systems.

Integrin-ß4 identifies cancer stem cell-enriched populations of partially mesenchymal carcinoma cells.

Neoplastic cells within individual carcinomas often exhibit considerable phenotypic heterogeneity in their epithelial versus mesenchymal-like cell states. Because carcinoma cells with mesenchymal features are often more resistant to therapy and may serve as a source of relapse, we sought to determine whether such cells could be further stratified into functionally distinct subtypes. Indeed, we find that a basal epithelial marker, integrin-ß4 (ITGB4), can be used to enable stratification of mesenchymal-like triple-negative breast cancer (TNBC) cells that differ from one another in their relative tumorigenic abilities. Notably, we demonstrate that ITGB4+ cancer stem cell (CSC)-enriched mesenchymal cells reside in an intermediate epithelial/mesenchymal phenotypic state. Among patients with TNBC who received chemotherapy, elevated ITGB4 expression was associated with a worse 5-year probability of relapse-free survival. Mechanistically, we find that the ZEB1 (zinc finger E-box binding homeobox 1) transcription factor activity in highly mesenchymal SUM159 TNBC cells can repress expression of the epithelial transcription factor TAp63α (tumor protein 63 isoform 1), a protein that promotes ITGB4 expression. In addition, we demonstrate that ZEB1 and ITGB4 are important in modulating the histopathological phenotypes of tumors derived from mesenchymal TNBC cells. Hence, mesenchymal carcinoma cell populations are internally heterogeneous, and ITGB4 is a mechanistically driven prognostic biomarker that can be used to identify the more aggressive subtypes of mesenchymal carcinoma cells in TNBC. The ability to rapidly isolate and mechanistically interrogate the CSC-enriched, partially mesenchymal carcinoma cells should further enable identification of novel therapeutic opportunities to improve the prognosis for high-risk patients with TNBC.

Characterization of aziridinylbenzoquinone DNA cross-links by liquid chromatography-infrared multiphoton dissociation-mass spectrometry.

DNA cross-linking was evaluated by liquid chromatography-tandem mass spectrometry to determine the relative cross-linking abilities of two aziridinylbenzoquinones. Reactivities of RH1 (2,5-diaziridinyl-3-[hydroxymethyl]-6-methyl-1,4-benzoquinone), a clinically studied antitumor cross-linking agent, and an analogue containing a phenyl group (2,5-diaziridinyl-3-[hydroxymethyl]-6-phenyl-1,4-benzoquinone, PhRH1) rather than a methyl group were compared. The bulky phenyl substituent was added to determine the impact of steric hindrance on the formation of cross-links within a double helical structure. Cross-links formed by RH1 and PhRH1 were observed at 5'-dGNC sites as well as 5'-dGAAC/dGTTC sites. RH1 was more effective at forming cross-links than PhRH1 for a variety of duplexes. Infrared multiphoton dissociation (IRMPD) and collision-induced dissociation results confirmed the presence and the location of the cross-links within the duplexes, and IRMPD was used to identify the dissociation pathways of the cross-linked duplexes.

Identification of the adduct between a 4-Aza-3-ene-1,6-diyne and DNA using electrospray ionization mass spectrometry.

The interactions between a novel enediyne [1-methyl-2-(phenylethynyl)-3-(3-phenylprop-2-ynyl)-3H-benzimidazolium] (1) and various cytosine-containing oligonucleotides were studied using electrospray ionization mass spectrometry (ESI-MS) in a flow injection analysis mode useful for small volumes. This enediyne ligand, developed as a potential alternative to the highly cytotoxic natural enediynes, some of which have been successfully used as anti-tumor agents, has previously been shown to interact with DNA through frank strand scission as well as via the formation of adducts that lead to 2'-deoxycytidine-specific cleavage. Through ESI-MS, the structures of these adducts were examined and a sequence dependence of the 2'-deoxycytidine-specific cleavage was noted. Collisionally activated dissociation of the observed adducts confirmed the strength of the interactions between the enediyne and DNA and supports a direct linkage between the enediyne and the cytosine nucleobase, likely the result of a nucleophilic attack of the phenylethynyl group by the cytosine amine.