Nfib Promotes Metastasis through a Widespread Increase in Chromatin Accessibility.

Metastases are the main cause of cancer deaths, but the mechanisms underlying metastatic progression remain poorly understood. We isolated pure populations of cancer cells from primary tumors and metastases from a genetically engineered mouse model of human small cell lung cancer (SCLC) to investigate the mechanisms that drive the metastatic spread of this lethal cancer. Genome-wide characterization of chromatin accessibility revealed the opening of large numbers of distal regulatory elements across the genome during metastatic progression. These changes correlate with copy number amplification of the Nfib locus, and differentially accessible sites were highly enriched for Nfib transcription factor binding sites. Nfib is necessary and sufficient to increase chromatin accessibility at a large subset of the intergenic regions. Nfib promotes pro-metastatic neuronal gene expression programs and drives the metastatic ability of SCLC cells. The identification of widespread chromatin changes during SCLC progression reveals an unexpected global reprogramming during metastatic progression.

Machine-learning-optimized Cas12a barcoding enables the recovery of single-cell lineages and transcriptional profiles.

The development of CRISPR-based barcoding methods creates an exciting opportunity to understand cellular phylogenies. We present a compact, tunable, high-capacity Cas12a barcoding system called dual acting inverted site array (DAISY). We combined high-throughput screening and machine learning to predict and optimize the 60-bp DAISY barcode sequences. After optimization, top-performing barcodes had ∼10-fold increased capacity relative to the best random-screened designs and performed reliably across diverse cell types. DAISY barcode arrays generated ∼12 bits of entropy and ∼66,000 unique barcodes. Thus, DAISY barcodes-at a fraction of the size of Cas9 barcodes-achieved high-capacity barcoding. We coupled DAISY barcoding with single-cell RNA-seq to recover lineages and gene expression profiles from ∼47,000 human melanoma cells. A single DAISY barcode recovered up to ∼700 lineages from one parental cell. This analysis revealed heritable single-cell gene expression and potential epigenetic modulation of memory gene transcription. Overall, Cas12a DAISY barcoding is an efficient tool for investigating cell-state dynamics.

miR-200 deficiency promotes lung cancer metastasis by activating Notch signaling in cancer-associated fibroblasts.

Lung adenocarcinoma, the most prevalent lung cancer subtype, is characterized by its high propensity to metastasize. Despite the importance of metastasis in lung cancer mortality, its underlying cellular and molecular mechanisms remain largely elusive. Here, we identified miR-200 miRNAs as potent suppressors for lung adenocarcinoma metastasis. miR-200 expression is specifically repressed in mouse metastatic lung adenocarcinomas, and miR-200 decrease strongly correlates with poor patient survival. Consistently, deletion of mir-200c/141 in the KrasLSL-G12D/+ ; Trp53flox/flox lung adenocarcinoma mouse model significantly promoted metastasis, generating a desmoplastic tumor stroma highly reminiscent of metastatic human lung cancer. miR-200 deficiency in lung cancer cells promotes the proliferation and activation of adjacent cancer-associated fibroblasts (CAFs), which in turn elevates the metastatic potential of cancer cells. miR-200 regulates the functional interaction between cancer cells and CAFs, at least in part, by targeting Notch ligand Jagged1 and Jagged2 in cancer cells and inducing Notch activation in adjacent CAFs. Hence, the interaction between cancer cells and CAFs constitutes an essential mechanism to promote metastatic potential.

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.

The AMBRA1 E3 ligase adaptor regulates the stability of cyclin D.

The initiation of cell division integrates a large number of intra- and extracellular inputs. D-type cyclins (hereafter, cyclin D) couple these inputs to the initiation of DNA replication1. Increased levels of cyclin D promote cell division by activating cyclin-dependent kinases 4 and 6 (hereafter, CDK4/6), which in turn phosphorylate and inactivate the retinoblastoma tumour suppressor. Accordingly, increased levels and activity of cyclin D-CDK4/6 complexes are strongly linked to unchecked cell proliferation and cancer2,3. However, the mechanisms that regulate levels of cyclin D are incompletely understood4,5. Here we show that autophagy and beclin 1 regulator 1 (AMBRA1) is the main regulator of the degradation of cyclin D. We identified AMBRA1 in a genome-wide screen to investigate the genetic basis of  the response to CDK4/6 inhibition. Loss of AMBRA1 results in high levels of cyclin D in cells and in mice, which promotes proliferation and decreases sensitivity to CDK4/6 inhibition. Mechanistically, AMBRA1 mediates ubiquitylation and proteasomal degradation of cyclin D as a substrate receptor for the cullin 4 E3 ligase complex. Loss of AMBRA1 enhances the growth of lung adenocarcinoma in a mouse model, and low levels of AMBRA1 correlate with worse survival in patients with lung adenocarcinoma. Thus, AMBRA1 regulates cellular levels of cyclin D, and contributes to cancer development and the response of cancer cells to CDK4/6 inhibitors.

Hardwiring tissue-specific AAV transduction in mice through engineered receptor expression.

The development of transgenic mouse models that express genes of interest in specific cell types has transformed our understanding of basic biology and disease. However, generating these models is time- and resource-intensive. Here we describe a model system, SELective Expression and Controlled Transduction In Vivo (SELECTIV), that enables efficient and specific expression of transgenes by coupling adeno-associated virus (AAV) vectors with Cre-inducible overexpression of the multi-serotype AAV receptor, AAVR. We demonstrate that transgenic AAVR overexpression greatly increases the efficiency of transduction of many diverse cell types, including muscle stem cells, which are normally refractory to AAV transduction. Superior specificity is achieved by combining Cre-mediated AAVR overexpression with whole-body knockout of endogenous Aavr, which is demonstrated in heart cardiomyocytes, liver hepatocytes and cholinergic neurons. The enhanced efficacy and exquisite specificity of SELECTIV has broad utility in development of new mouse model systems and expands the use of AAV for gene delivery in vivo.

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.

Fully accessible fitness landscape of oncogene-negative lung adenocarcinoma.

Cancer genomes are almost invariably complex with genomic alterations cooperating during each step of carcinogenesis. In cancers that lack a single dominant oncogene mutation, cooperation between the inactivation of multiple tumor suppressor genes can drive tumor initiation and growth. Here, we shed light on how the sequential acquisition of genomic alterations generates oncogene-negative lung tumors. We couple tumor barcoding with combinatorial and multiplexed somatic genome editing to characterize the fitness landscapes of three tumor suppressor genes NF1, RASA1, and PTEN, the inactivation of which jointly drives oncogene-negative lung adenocarcinoma initiation and growth. The fitness landscape was surprisingly accessible, with each additional mutation leading to growth advantage. Furthermore, the fitness landscapes remained fully accessible across backgrounds with the inactivation of additional tumor suppressor genes. These results suggest that while predicting cancer evolution will be challenging, acquiring the multiple alterations that drive the growth of oncogene-negative tumors can be facilitated by the lack of constraints on mutational order.

Publisher Correction: Towards quantitative and multiplexed in vivo functional cancer genomics.

The originally published article failed to acknowledge the equal first authorship contribution of I. P. Winters and C. W. Murray. The article has now been corrected online. The editors apologize for this error.

Towards quantitative and multiplexed in vivo functional cancer genomics.

Large-scale sequencing of human tumours has uncovered a vast array of genomic alterations. Genetically engineered mouse models recapitulate many features of human cancer and have been instrumental in assigning biological meaning to specific cancer-associated alterations. However, their time, cost and labour-intensive nature limits their broad utility; thus, the functional importance of the majority of genomic aberrations in cancer remains unknown. Recent advances have accelerated the functional interrogation of cancer-associated alterations within in vivo models. Specifically, the past few years have seen the emergence of CRISPR-Cas9-based strategies to rapidly generate increasingly complex somatic alterations and the development of multiplexed and quantitative approaches to ascertain gene function in vivo.

Targeted deletion reveals essential and overlapping functions of the miR-17 through 92 family of miRNA clusters.

miR-17 approximately 92, miR-106b approximately 25, and miR-106a approximately 363 belong to a family of highly conserved miRNA clusters. Amplification and overexpression of miR-1792 is observed in human cancers, and its oncogenic properties have been confirmed in a mouse model of B cell lymphoma. Here we show that mice deficient for miR-17 approximately 92 die shortly after birth with lung hypoplasia and a ventricular septal defect. The miR-17 approximately 92 cluster is also essential for B cell development. Absence of miR-17 approximately 92 leads to increased levels of the proapoptotic protein Bim and inhibits B cell development at the pro-B to pre-B transition. Furthermore, while ablation of miR-106b approximately 25 or miR-106a approximately 363 has no obvious phenotypic consequences, compound mutant embryos lacking both miR-106b approximately 25 and miR-17 approximately 92 die at midgestation. These results provide key insights into the physiologic functions of this family of microRNAs and suggest a link between the oncogenic properties of miR-17 approximately 92 and its functions during B lymphopoiesis and lung development.

Microbial single-strand annealing proteins enable CRISPR gene-editing tools with improved knock-in efficiencies and reduced off-target effects.

Several existing technologies enable short genomic alterations including generating indels and short nucleotide variants, however, engineering more significant genomic changes is more challenging due to reduced efficiency and precision. Here, we developed RecT Editor via Designer-Cas9-Initiated Targeting (REDIT), which leverages phage single-stranded DNA-annealing proteins (SSAP) RecT for mammalian genome engineering. Relative to Cas9-mediated homology-directed repair (HDR), REDIT yielded up to a 5-fold increase of efficiency to insert kilobase-scale exogenous sequences at defined genomic regions. We validated our REDIT approach using different formats and lengths of knock-in templates. We further demonstrated that REDIT tools using Cas9 nickase have efficient gene-editing activities and reduced off-target errors, measured using a combination of targeted sequencing, genome-wide indel, and insertion mapping assays. Our experiments inhibiting repair enzyme activities suggested that REDIT has the potential to overcome limitations of endogenous DNA repair steps. Finally, our REDIT method is applicable across cell types including human stem cells, and is generalizable to different Cas9 enzymes.

Pancreatic cancer modeling using retrograde viral vector delivery and in vivo CRISPR/Cas9-mediated somatic genome editing.

Pancreatic ductal adenocarcinoma (PDAC) is a genomically diverse, prevalent, and almost invariably fatal malignancy. Although conventional genetically engineered mouse models of human PDAC have been instrumental in understanding pancreatic cancer development, these models are much too labor-intensive, expensive, and slow to perform the extensive molecular analyses needed to adequately understand this disease. Here we demonstrate that retrograde pancreatic ductal injection of either adenoviral-Cre or lentiviral-Cre vectors allows titratable initiation of pancreatic neoplasias that progress into invasive and metastatic PDAC. To enable in vivo CRISPR/Cas9-mediated gene inactivation in the pancreas, we generated a Cre-regulated Cas9 allele and lentiviral vectors that express Cre and a single-guide RNA. CRISPR-mediated targeting of Lkb1 in combination with oncogenic Kras expression led to selection for inactivating genomic alterations, absence of Lkb1 protein, and rapid tumor growth that phenocopied Cre-mediated genetic deletion of Lkb1. This method will transform our ability to rapidly interrogate gene function during the development of this recalcitrant cancer.

An AMPK-independent signaling pathway downstream of the LKB1 tumor suppressor controls Snail1 and metastatic potential.

The serine/threonine kinase LKB1 is a tumor suppressor whose loss is associated with increased metastatic potential. In an effort to define biochemical signatures of metastasis associated with LKB1 loss, we discovered that the epithelial-to-mesenchymal transition transcription factor Snail1 was uniquely upregulated upon LKB1 deficiency across cell types. The ability of LKB1 to suppress Snail1 levels was independent of AMPK but required the related kinases MARK1 and MARK4. In a screen for substrates of these kinases involved in Snail regulation, we identified the scaffolding protein DIXDC1. Similar to loss of LKB1, DIXDC1 depletion results in upregulation of Snail1 in a FAK-dependent manner, leading to increased cell invasion. MARK1 phosphorylation of DIXDC1 is required for its localization to focal adhesions and ability to suppress metastasis in mice. DIXDC1 is frequently downregulated in human cancers, which correlates with poor survival. This study defines an AMPK-independent phosphorylation cascade essential for LKB1-dependent control of metastatic behavior.

Differential Tks5 isoform expression contributes to metastatic invasion of lung adenocarcinoma.

Metastasis accounts for the vast majority of cancer-related deaths, yet the molecular mechanisms that drive metastatic spread remain poorly understood. Here we report that Tks5, which has been linked to the formation of proteolytic cellular protrusions known as invadopodia, undergoes an isoform switch during metastatic progression in a genetically engineered mouse model of lung adenocarcinoma. Nonmetastatic primary tumor-derived cells predominantly expressed a short isoform, Tks5short, while metastatic primary tumor- and metastasis-derived cells acquired increased expression of the full-length isoform Tks5long. This elevation of Tks5long to Tks5short ratio correlated with a commensurate increase in invadopodia activity in metastatic cells compared with nonmetastatic cells. Further characterization of these isoforms by knockdown and overexpression experiments demonstrated that Tks5long promoted invadopodia in vitro and increased metastasis in transplant models and an autochthonous model of lung adenocarcinoma. Conversely, Tks5short decreased invadopodia stability and proteolysis, acting as a natural dominant-negative inhibitor to Tks5long. Importantly, high Tks5long and low Tks5short expressions in human lung adenocarcinomas correlated with metastatic disease and predicted worse survival of early stage patients. These data indicate that tipping the Tks5 isoform balance to a high Tks5long to Tks5short ratio promotes invadopodia-mediated invasion and metastasis.

A quantitative and multiplexed approach to uncover the fitness landscape of tumor suppression in vivo.

Cancer growth is a multistage, stochastic evolutionary process. While cancer genome sequencing has been instrumental in identifying the genomic alterations that occur in human tumors, the consequences of these alterations on tumor growth remain largely unexplored. Conventional genetically engineered mouse models enable the study of tumor growth in vivo, but they are neither readily scalable nor sufficiently quantitative to unravel the magnitude and mode of action of many tumor-suppressor genes. Here, we present a method that integrates tumor barcoding with ultradeep barcode sequencing (Tuba-seq) to interrogate tumor-suppressor function in mouse models of human cancer. Tuba-seq uncovers genotype-dependent distributions of tumor sizes. By combining Tuba-seq with multiplexed CRISPR-Cas9-mediated genome editing, we quantified the effects of 11 tumor-suppressor pathways that are frequently altered in human lung adenocarcinoma. Tuba-seq enables the broad quantification of the function of tumor-suppressor genes with unprecedented resolution, parallelization, and precision.

HMGA2 functions as a competing endogenous RNA to promote lung cancer progression.

Non-small-cell lung cancer (NSCLC) is the most prevalent histological cancer subtype worldwide. As the majority of patients present with invasive, metastatic disease, it is vital to understand the basis for lung cancer progression. Hmga2 is highly expressed in metastatic lung adenocarcinoma, in which it contributes to cancer progression and metastasis. Here we show that Hmga2 promotes lung cancer progression in mouse and human cells by operating as a competing endogenous RNA (ceRNA) for the let-7 microRNA (miRNA) family. Hmga2 can promote the transformation of lung cancer cells independent of protein-coding function but dependent upon the presence of let-7 sites; this occurs without changes in the levels of let-7 isoforms, suggesting that Hmga2 affects let-7 activity by altering miRNA targeting. These effects are also observed in vivo, where Hmga2 ceRNA activity drives lung cancer growth, invasion and dissemination. Integrated analysis of miRNA target prediction algorithms and metastatic lung cancer gene expression data reveals the TGF-ß co-receptor Tgfbr3 (ref. 12) as a putative target of Hmga2 ceRNA function. Tgfbr3 expression is regulated by the Hmga2 ceRNA through differential recruitment to Argonaute 2 (Ago2), and TGF-ß signalling driven by Tgfbr3 is important for Hmga2 to promote lung cancer progression. Finally, analysis of NSCLC-patient gene-expression data reveals that HMGA2 and TGFBR3 are coordinately regulated in NSCLC-patient material, a vital corollary to ceRNA function. Taken together, these results suggest that Hmga2 promotes lung carcinogenesis both as a protein-coding gene and as a non-coding RNA; such dual-function regulation of gene-expression networks reflects a novel means by which oncogenes promote disease progression.

Quantitative proteomics identify Tenascin-C as a promoter of lung cancer progression and contributor to a signature prognostic of patient survival.

The extracellular microenvironment is an integral component of normal and diseased tissues that is poorly understood owing to its complexity. To investigate the contribution of the microenvironment to lung fibrosis and adenocarcinoma progression, two pathologies characterized by excessive stromal expansion, we used mouse models to characterize the extracellular matrix (ECM) composition of normal lung, fibrotic lung, lung tumors, and metastases. Using quantitative proteomics, we identified and assayed the abundance of 113 ECM proteins, which revealed robust ECM protein signatures unique to fibrosis, primary tumors, or metastases. These analyses indicated significantly increased abundance of several S100 proteins, including Fibronectin and Tenascin-C (Tnc), in primary lung tumors and associated lymph node metastases compared with normal tissue. We further showed that Tnc expression is repressed by the transcription factor Nkx2-1, a well-established suppressor of metastatic progression. We found that increasing the levels of Tnc, via CRISPR-mediated transcriptional activation of the endogenous gene, enhanced the metastatic dissemination of lung adenocarcinoma cells. Interrogation of human cancer gene expression data revealed that high TNC expression correlates with worse prognosis for lung adenocarcinoma, and that a three-gene expression signature comprising TNC, S100A10, and S100A11 is a robust predictor of patient survival independent of age, sex, smoking history, and mutational load. Our findings suggest that the poorly understood ECM composition of the fibrotic and tumor microenvironment is an underexplored source of diagnostic markers and potential therapeutic targets for cancer patients.

An in vivo multiplexed small-molecule screening platform.

Phenotype-based small-molecule screening is a powerful method to identify molecules that regulate cellular functions. However, such screens are generally performed in vitro under conditions that do not necessarily model complex physiological conditions or disease states. Here, we use molecular cell barcoding to enable direct in vivo phenotypic screening of small-molecule libraries. The multiplexed nature of this approach allows rapid in vivo analysis of hundreds to thousands of compounds. Using this platform, we screened >700 covalent inhibitors directed toward hydrolases for their effect on pancreatic cancer metastatic seeding. We identified multiple hits and confirmed the relevant target of one compound as the lipase ABHD6. Pharmacological and genetic studies confirmed the role of this enzyme as a regulator of metastatic fitness. Our results highlight the applicability of this multiplexed screening platform for investigating complex processes in vivo.

Nuclear factor I/B is an oncogene in small cell lung cancer.

Small cell lung cancer (SCLC) is an aggressive cancer often diagnosed after it has metastasized. Despite the need to better understand this disease, SCLC remains poorly characterized at the molecular and genomic levels. Using a genetically engineered mouse model of SCLC driven by conditional deletion of Trp53 and Rb1 in the lung, we identified several frequent, high-magnitude focal DNA copy number alterations in SCLC. We uncovered amplification of a novel, oncogenic transcription factor, Nuclear factor I/B (Nfib), in the mouse SCLC model and in human SCLC. Functional studies indicate that NFIB regulates cell viability and proliferation during transformation.