A Functional Taxonomy of Tumor Suppression in Oncogenic KRAS-Driven Lung Cancer

Integrative multiomic approaches reveal ZMAT3 and p21 as conserved hubs in the p53 tumor suppression network

TP53, the most frequently mutated gene in human cancer, encodes a transcriptional activator that induces myriad downstream target genes. Despite the importance of p53 in tumor suppression, the specific p53 target genes important for tumor suppression remain unclear. Recent studies have identified the p53-inducible gene Zmat3 as a critical effector of tumor suppression, but many questions remain regarding its p53-dependence, activity across contexts, and mechanism of tumor suppression alone and in cooperation with other p53-inducible genes. To address these questions, we used Tuba-seqUltra somatic genome editing and tumor barcoding in a mouse lung adenocarcinoma model, combinatorial in vivo CRISPR/Cas9 screens, meta-analyses of gene expression and Cancer Dependency Map data, and integrative RNA-sequencing and shotgun proteomic analyses. We established Zmat3 as a core component of p53-mediated tumor suppression and identified Cdkn1a as the most potent cooperating p53-induced gene in tumor suppression. We discovered that ZMAT3/CDKN1A serve as near-universal effectors of p53-mediated tumor suppression that regulate cell division, migration, and extracellular matrix organization. Accordingly, combined Zmat3-Cdkn1a inactivation dramatically enhanced cell proliferation and migration compared to controls, akin to p53 inactivation. Together, our findings place ZMAT3 and CDKN1A as hubs of a p53-induced gene program that opposes tumorigenesis across various cellular and genetic contexts.

In Vivo Multiplexed Modeling Reveals Diverse Roles of the TBX2 Subfamily and Egr1 in Ras -Driven Lung Adenocarcinoma

The TBX2 subfamily of T-box transcription factors (including Tbx2 , Tbx3 , Tbx4 , Tbx5 ) plays an essential role in lung development. Downregulation of these genes in human Lung adenocarcinoma (LUAD) suggests that these genes may be tumor suppressive, however because downregulation appears to occur primarily via epigenetic change, it remains unclear if these changes causally drive tumor progression or are merely the consequence of upstream events. Herein, we developed the first multiplexed mouse model to study the impact of TBX2 subfamily loss, alongside associated signaling genes Egr1 , Chd2 , Tnfaip3a , and Atf3 , in Ras -driven lung cancer. Using TuBa-seq, a high-throughput tumor-barcoding system, we quantified the growth effects of these knockouts during early and late tumorigenesis. Chd2 loss consistently suppressed tumor progression, while Tbx2 loss exhibited stage-dependent effects. Notably, Egr1 emerged as a potent tumor suppressor, with its knockout increasing tumor size (∼5x) at 20 weeks, surpassing Rb1 loss. Transcriptomic analyses of Egr1 -deficient tumors suggested immune dysregulation, including heightened inflammation and potential markers of T cell exhaustion in the tumor microenvironment. These findings indicate that Egr1 may play a role in suppressing tumor growth through modulating immune dynamics, offering new insights into the interplay between tumor progression and immune regulation in LUAD.

Combinatorial In Vivo Genome Editing Identifies Widespread Epistasis and an Accessible Fitness Landscape During Lung Tumorigenesis

Lung adenocarcinoma, the most common subtype of lung cancer, is genomically complex, with tumors containing tens to hundreds of non-synonymous mutations. However, little is understood about how genes interact with each other to enable the evolution of cancer in vivo, largely due to a lack of methods for investigating genetic interactions in a high-throughput and quantitative manner. Here, we employed a novel platform to generate tumors with inactivation of pairs of ten diverse tumor suppressor genes within an autochthonous mouse model of oncogenic KRAS-driven lung cancer. By quantifying the fitness of tumors with every single and double mutant genotype, we show that most tumor suppressor genetic interactions exhibited negative epistasis, with diminishing returns on tumor fitness. In contrast, Apc inactivation showed positive epistasis with the inactivation of several other genes, including synergistic effects on tumor fitness in combination with Lkb1 or Nf1 inactivation. Sign epistasis was extremely rare, suggesting a surprisingly accessible fitness landscape during lung tumorigenesis. These findings expand our understanding of the interactions that drive tumorigenesis in vivo.

Emerging roles of cohesin-STAG2 in cancer

Cohesin, a crucial regulator of genome organisation, plays a fundamental role in maintaining chromatin architecture as well as gene expression. Among its subunits, STAG2 stands out because of its frequent deleterious mutations in various cancer types, such as bladder cancer and melanoma. Loss of STAG2 function leads to significant alterations in chromatin structure, disrupts transcriptional regulation, and impairs DNA repair pathways. In this review, we explore the molecular mechanisms underlying cohesin-STAG2 function, highlighting its roles in healthy cells and its contributions to cancer biology, showing how STAG2 dysfunction promotes tumourigenesis and presents opportunities for targeted therapeutic interventions.

Tumor suppressor genotype influences the extent and mode of immunosurveillance in lung cancer

The impact of cancer driving mutations in regulating immunosurveillance throughout tumor development remains poorly understood. To better understand the contribution of tumor genotype to immunosurveillance, we generated and validated lentiviral vectors that create an epi-allelic series of increasingly immunogenic neoantigens. This vector system is compatible with autochthonous Cre-regulated cancer models, CRISPR/Cas9-mediated somatic genome editing, and tumor barcoding. Here, we show that in the context of KRAS-driven lung cancer and strong neoantigen expression, tumor suppressor genotype dictates the degree of immune cell recruitment, positive selection of tumors with neoantigen silencing, and tumor outgrowth. By quantifying the impact of 11 commonly inactivated tumor suppressor genes on tumor growth across neoantigenic contexts, we show that the growth promoting effects of tumor suppressor gene inactivation correlate with increasing sensitivity to immunosurveillance. Importantly, specific genotypes dramatically increase or decrease sensitivity to immunosurveillance independently of their growth promoting effects. We propose a model of immunoediting in which tumor suppressor gene inactivation works in tandem with neoantigen expression to shape tumor immunosurveillance and immunoediting such that the same neoantigens uniquely modulate tumor immunoediting depending on the genetic context.

One Sentence Summary

Here we uncover an under-appreciated role for tumor suppressor gene inactivation in shaping immunoediting upon neoantigen expression.

A STAG2-PAXIP1/PAGR1 axis suppresses lung tumorigenesis

The cohesin complex is a critical regulator of gene expression. STAG2 is the most frequently mutated cohesin subunit across several cancer types and is a key tumor suppressor in lung cancer. Here, we coupled somatic CRISPR-Cas9 genome editing and tumor barcoding with an autochthonous oncogenic KRAS-driven lung cancer model and showed that STAG2 is uniquely tumor-suppressive among all core and auxiliary cohesin components. The heterodimeric complex components PAXIP1 and PAGR1 have highly correlated effects with STAG2 in human lung cancer cell lines, are tumor suppressors in vivo, and are epistatic to STAG2 in oncogenic KRAS-driven lung tumorigenesis in vivo. STAG2 inactivation elicits changes in gene expression, chromatin accessibility, and 3D genome conformation that impact the cancer cell state. Gene expression and chromatin accessibility similarities between STAG2- and PAXIP1-deficient neoplastic cells further relate STAG2-cohesin to PAXIP1/PAGR1. These findings reveal a STAG2-PAXIP1/PAGR1 tumor-suppressive axis and uncover novel PAXIP1-dependent and PAXIP1-independent STAG2-cohesin-mediated mechanisms of lung tumor suppression.

Kelch-like ECH-associated Protein 1/Nuclear Factor Erythroid 2-related Factor 2 Pathway and Its Interplay with Oncogenes in Lung Tumorigenesis

Nuclear factor erythroid 2-related factor 2 (NRF2), a transcription factor regulating cellular redox homeostasis, exhibits a complex role in cancer biology. Genetic mutations in the Kelch-like ECH-associated protein 1 (KEAP1)/NRF2 system, which lead to NRF2 hyperactivation, are found in 20% to 30% of lung cancer cases. This review explores the intricate interplay between NRF2 and key oncogenic pathways in lung cancer, focusing on the interaction of KEAP1/NRF2 system with Kirsten rat sarcoma virus (KRAS), tumor protein P53 (TP53), epidermal growth factor receptor (EGFR), and phosphatidylinositol 3-kinases (PI3K)/AKT signaling. While NRF2 activation alone is insufficient to initiate tumorigenesis, it can significantly impact tumor initiation and progression when combined with oncogenic drivers such as KRAS. The review highlights the context-dependent effects of NRF2, from its protective role against chemical carcinogen-induced tumor initiation to its potential promotion of tumor growth in established cancers. These findings suggest the need for nuanced, stage-specific approaches to targeting the NRF2 pathway in cancer therapy.

An essential role for Cmtr2 in mammalian embryonic development

CMTR2 is an mRNA cap methyltransferase with poorly understood physiological functions. It catalyzes 2'-O-ribose methylation of the second transcribed nucleotide of mRNAs, potentially serving to mark RNAs as "self" to evade the cellular innate immune response. Here we analyze the consequences of Cmtr2 deficiency in mice. We discover that constitutive deletion of Cmtr2 results in mouse embryos that die during mid-gestation, exhibiting defects in embryo size, placental malformation and yolk sac vascularization. Endothelial cell deletion of Cmtr2 in mice results in vascular and hematopoietic defects, and perinatal lethality. Detailed characterization of the constitutive Cmtr2 KO phenotype shows an activation of the p53 pathway and decreased proliferation, but no evidence of interferon pathway activation. In summary, our study reveals the essential roles of Cmtr2 in mammalian cells beyond its immunoregulatory function.

Gene-based burden tests of rare germline variants identify six cancer susceptibility genes

Discovery of cancer risk variants in the sequence of the germline genome can shed light on carcinogenesis. Here we describe gene burden association analyses, aggregating rare missense and loss of function variants, at 22 cancer sites, including 130,991 cancer cases and 733,486 controls from Iceland, Norway and the United Kingdom. We identified four genes associated with increased cancer risk; the pro-apoptotic BIK for prostate cancer, the autophagy involved ATG12 for colorectal cancer, TG for thyroid cancer and CMTR2 for both lung cancer and cutaneous melanoma. Further, we found genes with rare variants that associate with decreased risk of cancer; AURKB for any cancer, irrespective of site, and PPP1R15A for breast cancer, suggesting that inhibition of PPP1R15A may be a preventive strategy for breast cancer. Our findings pinpoint several new cancer risk genes and emphasize autophagy, apoptosis and cell stress response as a focus point for developing new therapeutics.

Integrative multiomic approaches reveal ZMAT3 and p21 as conserved hubs in the p53 tumor suppression network

TP53 , the most frequently mutated gene in human cancer, encodes a transcriptional activator that induces myriad downstream target genes. Despite the importance of p53 in tumor suppression, the specific p53 target genes important for tumor suppression remain unclear. Recent studies have identified the p53-inducible gene Zmat3 as a critical effector of tumor suppression, but many questions remain regarding its p53-dependence, activity across contexts, and mechanism of tumor suppression alone and in cooperation with other p53-inducible genes. To address these questions, we used Tuba-seq Ultra somatic genome editing and tumor barcoding in a mouse lung adenocarcinoma model, combinatorial in vivo CRISPR/Cas9 screens, meta-analyses of gene expression and Cancer Dependency Map data, and integrative RNA-sequencing and shotgun proteomic analyses. We established Zmat3 as a core component of p53-mediated tumor suppression and identified Cdkn1a as the most potent cooperating p53-induced gene in tumor suppression. We discovered that ZMAT3/CDKN1A serve as near-universal effectors of p53-mediated tumor suppression that regulate cell division, migration, and extracellular matrix organization. Accordingly, combined Zmat3 - Cdkn1a inactivation dramatically enhanced cell proliferation and migration compared to controls, akin to p53 inactivation. Together, our findings place ZMAT3 and CDKN1A as hubs of a p53-induced gene program that opposes tumorigenesis across various cellular and genetic contexts.

Functional mapping of epigenetic regulators uncovers coordinated tumor suppression by the HBO1 and MLL1 complexes

Epigenetic dysregulation is widespread in cancer. However, the specific epigenetic regulators and the processes they control to drive cancer phenotypes are poorly understood. Here, we employed a novel, scalable and high-throughput in vivo method to perform iterative functional screens of over 250 epigenetic regulatory genes within autochthonous oncogenic KRAS-driven lung tumors. We identified multiple novel epigenetic tumor suppressor and tumor dependency genes. We show that a specific HBO1 complex and the MLL1 complex are among the most impactful tumor suppressive epigenetic regulators in lung. The histone modifications generated by the HBO1 complex are frequently absent or reduced in human lung adenocarcinomas. The HBO1 and MLL1 complexes regulate chromatin accessibility of shared genomic regions, lineage fidelity and the expression of canonical tumor suppressor genes. The HBO1 and MLL1 complexes are epistatic during lung tumorigenesis, and their functional correlation is conserved in human cancer cell lines. Together, these results demonstrate the value of quantitative methods to generate a phenotypic roadmap of epigenetic regulatory genes in tumorigenesis in vivo .

Aging represses lung tumorigenesis and alters tumor suppression

Most cancers are diagnosed in persons over the age of sixty, but little is known about how age impacts tumorigenesis. While aging is accompanied by mutation accumulation - widely understood to contribute to cancer risk - it is also associated with numerous other cellular and molecular changes likely to impact tumorigenesis. Moreover, cancer incidence decreases in the oldest part of the population, suggesting that very old age may reduce carcinogenesis. Here we show that aging represses tumor initiation and growth in genetically engineered mouse models of human lung cancer. Moreover, aging dampens the impact of inactivating many, but not all, tumor suppressor genes with the impact of inactivating PTEN, a negative regulator of the PI3K/AKT pathway, weakened to a disproportionate extent. Single-cell transcriptomic analysis revealed that neoplastic cells from tumors in old mice retain many age-related transcriptomic changes, showing that age has an enduring impact that persists through oncogenic transformation. Furthermore, the consequences of PTEN inactivation were strikingly age-dependent, with PTEN deficiency reducing signatures of aging in cancer cells and the tumor microenvironment. Our findings suggest that the relationship between age and lung cancer incidence may reflect an integration of the competing effects of driver mutation accumulation and tumor suppressive effects of aging.

Efficient and multiplexed somatic genome editing with Cas12a mice

Somatic genome editing in mouse models has increased our understanding of the in vivo effects of genetic alterations in areas ranging from neuroscience to cancer biology and beyond. However, existing models are limited in their ability to create multiple targeted edits. Thus, our understanding of the complex genetic interactions that underlie development, homeostasis, and disease remains incomplete. Cas12a is an RNA-guided endonuclease with unique attributes that enable simple targeting of multiple genes with crRNA arrays containing tandem guides. To accelerate and expand the generation of complex genotypes in somatic cells, we generated transgenic mice with Cre-regulated and constitutive expression of enhanced Acidaminococcus sp. Cas12a (enAsCas12a). In these mice, enAsCas12a-mediated somatic genome editing robustly generated compound genotypes, as exemplified by the initiation of diverse cancer types driven by homozygous inactivation of trios of tumor suppressor genes. We further integrated these modular crRNA arrays with clonal barcoding to quantify the size and number of tumors with each array, as well as the efficiency of each crRNA. These Cas12a alleles will enable the rapid generation of disease models and broadly facilitate the high-throughput investigation of coincident genomic alterations in somatic cells in vivo .

How CRISPR Is Revolutionizing the Generation of New Models for Cancer Research

Cancers arise through acquisition of mutations in genes that regulate core biological processes like cell proliferation and cell death. Decades of cancer research have led to the identification of genes and mutations causally involved in disease development and evolution, yet defining their precise function across different cancer types and how they influence therapy responses has been challenging. Mouse models have helped define the in vivo function of cancer-associated alterations, and genome-editing approaches using CRISPR have dramatically accelerated the pace at which these models are developed and studied. Here, we highlight how CRISPR technologies have impacted the development and use of mouse models for cancer research and discuss the many ways in which these rapidly evolving platforms will continue to transform our understanding of this disease.

Combinatorial in vivo genome editing identifies widespread epistasis during lung tumorigenesis

Lung adenocarcinoma, the most common subtype of lung cancer, is genomically complex, with tumors containing tens to hundreds of non-synonymous mutations. However, little is understood about how genes interact with each other to enable tumorigenesis in vivo , largely due to a lack of methods for investigating genetic interactions in a high-throughput and multiplexed manner. Here, we employed a novel platform to generate tumors with all pairwise inactivation of ten tumor suppressor genes within an autochthonous mouse model of oncogenic KRAS-driven lung cancer. By quantifying the fitness of tumors with every single and double mutant genotype, we show that most tumor suppressor genetic interactions exhibited negative epistasis, with diminishing returns on tumor fitness. In contrast, Apc inactivation showed positive epistasis with the inactivation of several other genes, including dramatically synergistic effects on tumor fitness in combination with Lkb1 or Nf1 inactivation. This approach has the potential to expand the scope of genetic interactions that may be functionally characterized in vivo , which could lead to a better understanding of how complex tumor genotypes impact each step of carcinogenesis.

Regulation of antioxidants in cancer

Scientists in this field often joke, "If you don't have a mechanism, say it's ROS." Seemingly connected to every biological process ever described, reactive oxygen species (ROS) have numerous pleiotropic roles in physiology and disease. In some contexts, ROS act as secondary messengers, controlling a variety of signaling cascades. In other scenarios, they initiate damage to macromolecules. Finally, in their worst form, ROS are deadly to cells and surrounding tissues. A set of molecules with detoxifying abilities, termed antioxidants, is the direct counterpart to ROS. Notably, antioxidants exist in the public domain, touted as a "cure-all" for diseases. Research has disproved many of these claims and, in some cases, shown the opposite. Of all the diseases, cancer stands out in its paradoxical relationship with antioxidants. Although the field has made numerous strides in understanding the roles of antioxidants in cancer, many questions remain.

Applications of CRISPR screening to lung cancer treatment

Lung cancer is an extremely aggressive and highly prevalent disease worldwide, and it is one of the leading causes of cancer death. Deciphering intrinsic genetic mechanism, finding new targets, and overcoming drug resistance are the key to lung cancer treatment. High-throughput CRISPR screening has been extensively used to obtain the genes related to cancers including lung cancer. This review describes CRISPR/Cas9 or CRISPR/dCas9-based technologies for high-throughput screening. We summarize the applications of CRISPR screening technology in exploring the mechanism of lung cancer development in vivo or in vitro, overcoming drug resistance, improving the effect of immunotherapy, and discovering new therapeutic targets. This review highlights the potential of CRISPR screening in combination with tumor barcoding and high-throughput sequencing (Tuba-seq) to precisely quantify the impact of alterations in many tumor suppressor genes on lung cancer.

Applications and advancements of CRISPR-Cas in the treatment of lung cancer

Lung cancer is one of the most malignant diseases and a major contributor to cancer-related deaths worldwide due to the deficiency of early diagnosis and effective therapy that are of great importance for patient prognosis and quality of life. Over the past decade, the advent of clustered regularly interspaced short palindromic repeats/CRISPR associated protein (CRISPR/Cas) system has significantly propelled the progress of both fundamental research and clinical trials of lung cancer. In this review, we review the current applications of the CRISPR/Cas system in diagnosis, target identification, and treatment resistance of lung cancer. Furthermore, we summarize the development of lung cancer animal models and delivery methods based on CRISPR system, providing novel insights into clinical diagnosis and treatment strategies of lung cancer.

Oncogenic context shapes the fitness landscape of tumor suppression

Tumors acquire alterations in oncogenes and tumor suppressor genes in an adaptive walk through the fitness landscape of tumorigenesis. However, the interactions between oncogenes and tumor suppressor genes that shape this landscape remain poorly resolved and cannot be revealed by human cancer genomics alone. Here, we use a multiplexed, autochthonous mouse platform to model and quantify the initiation and growth of more than one hundred genotypes of lung tumors across four oncogenic contexts: KRAS G12D, KRAS G12C, BRAF V600E, and EGFR L858R. We show that the fitness landscape is rugged-the effect of tumor suppressor inactivation often switches between beneficial and deleterious depending on the oncogenic context-and shows no evidence of diminishing-returns epistasis within variants of the same oncogene. These findings argue against a simple linear signaling relationship amongst these three oncogenes and imply a critical role for off-axis signaling in determining the fitness effects of inactivating tumor suppressors.

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.

Liquid biopsy in detecting early non-small cell lung cancer

Lung cancer screening programs, particularly in the UK, have shown a decrease in lung cancer-related deaths among individuals who underwent low-dose computed tomography (CT) screening. Researchers are now focusing on evaluating cell-free DNA through various methods to determine if pre-diagnostic mutations can be detected years before clinical diagnosis. This could help identify individuals at high risk of developing lung cancer. However, while this approach has successfully identified precursors of follicular lymphoma, the presence of occult lung preneoplasia in non-small-cell lung cancer still requires further investigation. The TRACERx consortium is conducting extensive research to comprehensively assess the detection and progression of non-small cell lung cancers (NSCLC). Liquid biopsy is being used in advanced stages of the disease to monitor disease progression, predict treatment response, and identify targetable driver oncogenic mutations and fusion genes. Intense research is also underway to identify numerous diagnostic gene signatures with high accuracy for early-stage lung cancer. However, a more focused clinical approach is needed, with a mechanistic focus on the key pathways of cancer development. Loss of liver kinase B1 (LKB1) function and deactivation due to 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanol (NNAL), a metabolite of tobacco-specific carcinogens, could potentially be traced and contribute to the development of new biomarkers. This testing could complement machine-learning approaches. The discovery of epidermal growth factor receptor (EGFR) and Kirsten rat sarcoma viral oncogene homolog (KRAS) mutations in healthy lung tissues by TRACERx investigators may also lead to the development of novel diagnostic tools. Tumor protein 53 (TP53) loss should also be considered as a marker that could contribute to malignant transformation. Intercepting aggressive non-small-cell lung cancer is a pressing priority. In this review, we discuss our experience and explore other research on exosomes and plasma circular RNA as potential biomarkers. Circular RNAs, formed through non-sequential back-splicing of pre-mRNA transcripts, play a role in epithelial-mesenchymal transition, with many of them regulated by the RNA-binding protein Quaking. Platelet RNA has shown promise in detecting early and late-stage cancer. The extensive exploration of liquid biopsy aims to provide affordable methods for tracing circulating precursors of non-small-cell lung cancer, highlighting the importance of its mission.

A functional sgRNA-CRISPR screening method for generating murine RET and NTRK1 rearranged oncogenes

CRISPR/Cas9 gene editing represents a powerful tool for investigating fusion oncogenes in cancer biology. Successful experiments require that sgRNAs correctly associate with their target sequence and initiate double stranded breaks which are subsequently repaired by endogenous DNA repair systems yielding fusion chromosomes. Simple tests to ensure sgRNAs are functional are not generally available and often require single cell cloning to identify successful CRISPR-editing events. Here, we describe a novel method relying on acquisition of IL3-independence in Ba/F3 cells to identify sgRNA pairs that generate oncogenic gene rearrangements of the Ret and Ntrk1 tyrosine kinases. The rearrangements were confirmed with PCR, RT-PCR and sequencing and Ba/F3 cells harboring Ret or Ntrk1 rearrangements acquired sensitivity to RET and TRK inhibitors, respectively. Adenoviruses encoding Cas9 and sgRNA pairs inducing the Kif5b-Ret and Trim24-Ret rearrangements were intratracheally instilled into mice and yielded lung adenocarcinomas. A cell line (TR.1) established from a Trim24-Ret positive tumor exhibited high in vitro sensitivity to the RET inhibitors LOXO-292 and BLU-667 and orthotopic TR.1 cell-derived tumors underwent marked shrinkage upon LOXO-292 treatment. Thus, the method offers an efficient means to validate sgRNAs that successfully target their intended loci for the generation of novel, syngeneic murine oncogene-driven tumor models.

High-Throughput Identification, Modeling, and Analysis of Cancer Driver Genes In Vivo

The vast number of genomic and molecular alterations in cancer pose a substantial challenge to uncovering the mechanisms of tumorigenesis and identifying therapeutic targets. High-throughput functional genomic methods in genetically engineered mouse models allow for rapid and systematic investigation of cancer driver genes. In this review, we discuss the basic concepts and tools for multiplexed investigation of functionally important cancer genes in vivo using autochthonous cancer models. Furthermore, we highlight emerging technical advances in the field, potential opportunities for future investigation, and outline a vision for integrating multiplexed genetic perturbations with detailed molecular analyses to advance our understanding of the genetic and molecular basis of cancer.

Distinct Nrf2 Signaling Thresholds Mediate Lung Tumor Initiation and Progression

Mutations in the KEAP1-NRF2 (Kelch-like ECH-associated protein 1-nuclear factor-erythroid 2 p45-related factor 2) pathway occur in up to a third of non-small cell lung cancer (NSCLC) cases and often confer resistance to therapy and poor outcomes. Here, we developed murine alleles of the KEAP1 and NRF2 mutations found in human NSCLC and comprehensively interrogated their impact on tumor initiation and progression. Chronic NRF2 stabilization by Keap1 or Nrf2 mutation was not sufficient to induce tumorigenesis, even in the absence of tumor suppressors, p53 or LKB1. When combined with KrasG12D/+, constitutive NRF2 activation promoted lung tumor initiation and early progression of hyperplasia to low-grade tumors but impaired their progression to advanced-grade tumors, which was reversed by NRF2 deletion. Finally, NRF2 overexpression in KEAP1 mutant human NSCLC cell lines was detrimental to cell proliferation, viability, and anchorage-independent colony formation. Collectively, these results establish the context-dependence and activity threshold for NRF2 during the lung tumorigenic process.

SIGNIFICANCE

Stabilization of the transcription factor NRF2 promotes oncogene-driven tumor initiation but blocks tumor progression, indicating distinct, threshold-dependent effects of the KEAP1/NRF2 pathway in different stages of lung tumorigenesis.

Dissecting metastasis using preclinical models and methods

Metastasis has long been understood to lead to the overwhelming majority of cancer-related deaths. However, our understanding of the metastatic process, and thus our ability to prevent or eliminate metastases, remains frustratingly limited. This is largely due to the complexity of metastasis, which is a multistep process that likely differs across cancer types and is greatly influenced by many aspects of the in vivo microenvironment. In this Review, we discuss the key variables to consider when designing assays to study metastasis: which source of metastatic cancer cells to use and where to introduce them into mice to address different questions of metastasis biology. We also examine methods that are being used to interrogate specific steps of the metastatic cascade in mouse models, as well as emerging techniques that may shed new light on previously inscrutable aspects of metastasis. Finally, we explore approaches for developing and using anti-metastatic therapies, and how mouse models can be used to test them.

A Rapid, Functional sgRNA Screening Method for Generating Murine RET and NTRK1 Fusion Oncogenes

CRISPR/Cas9 gene editing technology is an indispensable and powerful tool in the field of cancer biology. To conduct successful CRISPR-based experiments, it is crucial that sgRNAs generate their designed alterations. Here, we describe a simple and efficient sgRNA screening method for validating sgRNAs that generate oncogenic gene rearrangements. We used IL3-independence in Ba/F3 cells as an assay to identify sgRNA pairs that generate fusion oncogenes involving the Ret and Ntrk1 tyrosine kinases. We confirmed these rearrangements with PCR or RT-PCR as well as sequencing. Ba/F3 cells harboring Ret or Ntrk1 rearrangements acquired sensitivity to RET and TRK inhibitors, respectively. Adenoviruses encoding Cas9 and sgRNAs that catalyze the Kif5b-Ret and Trim24-Ret rearrangements were intratracheally instilled into mice and yielded lung adenocarcinomas. A cell line (TR.1) was established from a Trim24-Ret positive tumor that exhibited high in vitro sensitivity to RET-specific TKIs. Moreover, orthotopic transplantation of TR.1 cells into the left lung yielded well-defined tumors that shrank in response to LOXO-292 treatment. The method offers an efficient means to validate sgRNAs that successfully target their intended loci for the generation of novel murine oncogene-driven tumor models.

A multi-organoid platform identifies CIART as a key factor for SARS-CoV-2 infection

COVID-19 is a systemic disease involving multiple organs. We previously established a platform to derive organoids and cells from human pluripotent stem cells to model SARS-CoV-2 infection and perform drug screens^1,2. This provided insight into cellular tropism and the host response, yet the molecular mechanisms regulating SARS-CoV-2 infection remain poorly defined. Here we systematically examined changes in transcript profiles caused by SARS-CoV-2 infection at different multiplicities of infection for lung airway organoids, lung alveolar organoids and cardiomyocytes, and identified several genes that are generally implicated in controlling SARS-CoV-2 infection, including CIART, the circadian-associated repressor of transcription. Lung airway organoids, lung alveolar organoids and cardiomyocytes derived from isogenic CIART^-/- human pluripotent stem cells were significantly resistant to SARS-CoV-2 infection, independently of viral entry. Single-cell RNA-sequencing analysis further validated the decreased levels of SARS-CoV-2 infection in ciliated-like cells of lung airway organoids. CUT&RUN, ATAC-seq and RNA-sequencing analyses showed that CIART controls SARS-CoV-2 infection at least in part through the regulation of NR4A1, a gene also identified from the multi-organoid analysis. Finally, transcriptional profiling and pharmacological inhibition led to the discovery that the Retinoid X Receptor pathway regulates SARS-CoV-2 infection downstream of CIART and NR4A1. The multi-organoid platform identified the role of circadian-clock regulation in SARS-CoV-2 infection, which provides potential therapeutic targets for protection against COVID-19 across organ systems.

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 additional tumor suppressor mutations. These results suggest that while predicting cancer evolution will be challenging, acquiring the multiple alterations required for the growth of oncogene-negative tumors can be facilitated by the lack of constraints on mutational order.

mRNA ageing shapes the Cap2 methylome in mammalian mRNA

The mRNA cap structure is a major site of dynamic mRNA methylation. mRNA caps exist in either the Cap1 or Cap2 form, depending on the presence of 2'-O-methylation on the first transcribed nucleotide or both the first and second transcribed nucleotides, respectively1,2. However, the identity of Cap2-containing mRNAs and the function of Cap2 are unclear. Here we describe CLAM-Cap-seq, a method for transcriptome-wide mapping and quantification of Cap2. We find that unlike other epitranscriptomic modifications, Cap2 can occur on all mRNAs. Cap2 is formed through a slow continuous conversion of mRNAs from Cap1 to Cap2 as mRNAs age in the cytosol. As a result, Cap2 is enriched on long-lived mRNAs. Large increases in the abundance of Cap1 leads to activation of RIG-I, especially in conditions in which expression of RIG-I is increased. The methylation of Cap1 to Cap2 markedly reduces the ability of RNAs to bind to and activate RIG-I. The slow methylation rate of Cap2 allows Cap2 to accumulate on host mRNAs, yet ensures that low levels of Cap2 occur on newly expressed viral RNAs. Overall, these results reveal an immunostimulatory role for Cap1, and that Cap2 functions to reduce activation of the innate immune response.

A multiplexed in vivo approach to identify driver genes in small cell lung cancer

Small cell lung cancer (SCLC) is a lethal form of lung cancer. Here, we develop a quantitative multiplexed approach on the basis of lentiviral barcoding with somatic CRISPR-Cas9-mediated genome editing to functionally investigate candidate regulators of tumor initiation and growth in genetically engineered mouse models of SCLC. We found that naphthalene pre-treatment enhances lentiviral vector-mediated SCLC initiation, enabling high multiplicity of tumor clones for analysis through high-throughput sequencing methods. Candidate drivers of SCLC identified from a meta-analysis across multiple human SCLC genomic datasets were tested using this approach, which defines both positive and detrimental impacts of inactivating 40 genes across candidate pathways on SCLC development. This analysis and subsequent validation in human SCLC cells establish TSC1 in the PI3K-AKT-mTOR pathway as a robust tumor suppressor in SCLC. This approach should illuminate drivers of SCLC, facilitate the development of precision therapies for defined SCLC genotypes, and identify therapeutic targets.

Multiplexed screens identify RAS paralogues HRAS and NRAS as suppressors of KRAS-driven lung cancer growth

Oncogenic KRAS mutations occur in approximately 30% of lung adenocarcinoma. Despite several decades of effort, oncogenic KRAS-driven lung cancer remains difficult to treat, and our understanding of the regulators of RAS signalling is incomplete. Here to uncover the impact of diverse KRAS-interacting proteins on lung cancer growth, we combined multiplexed somatic CRISPR/Cas9-based genome editing in genetically engineered mouse models with tumour barcoding and high-throughput barcode sequencing. Through a series of CRISPR/Cas9 screens in autochthonous lung cancer models, we show that HRAS and NRAS are suppressors of KRASG12D-driven tumour growth in vivo and confirm these effects in oncogenic KRAS-driven human lung cancer cell lines. Mechanistically, RAS paralogues interact with oncogenic KRAS, suppress KRAS-KRAS interactions, and reduce downstream ERK signalling. Furthermore, HRAS and NRAS mutations identified in oncogenic KRAS-driven human tumours partially abolished this effect. By comparing the tumour-suppressive effects of HRAS and NRAS in oncogenic KRAS- and oncogenic BRAF-driven lung cancer models, we confirm that RAS paralogues are specific suppressors of KRAS-driven lung cancer in vivo. Our study outlines a technological avenue to uncover positive and negative regulators of oncogenic KRAS-driven cancer in a multiplexed manner in vivo and highlights the role RAS paralogue imbalance in oncogenic KRAS-driven lung cancer.

A new system for multiplexed mosaic analysis of gene function in the mouse

In a recent issue of Cell, Liu et al. present an innovative mouse model system in which Cre/lox stochastically turns on transgenic expression of one out of up to 100 sgRNAs in somatic cells, creating genetic mosaicism that enables the multiplexed assessment of gene function in vivo.

Tutorial: design and execution of CRISPR in vivo screens

Here we provide a detailed tutorial on CRISPR in vivo screening. Using the mouse as the model organism, we introduce a range of CRISPR tools and applications, delineate general considerations for 'transplantation-based' or 'direct in vivo' screening design, and provide details on technical execution, sequencing readouts, computational analyses and data interpretation. In vivo screens face unique pitfalls and limitations, such as delivery issues or library bottlenecking, which must be counteracted to avoid screening failure or flawed conclusions. A broad variety of in vivo phenotypes can be interrogated such as organ development, hematopoietic lineage decision and evolutionary licensing in oncogenesis. We describe experimental strategies to address various biological questions and provide an outlook on emerging CRISPR applications, such as genetic interaction screening. These technological advances create potent new opportunities to dissect the molecular underpinnings of complex organismal phenotypes.

Combinatorial Inactivation of Tumor Suppressors Efficiently Initiates Lung Adenocarcinoma with Therapeutic Vulnerabilities

Lung cancer is the leading cause of cancer death worldwide, with lung adenocarcinoma being the most common subtype. Many oncogenes and tumor suppressor genes are altered in this cancer type, and the discovery of oncogene mutations has led to the development of targeted therapies that have improved clinical outcomes. However, a large fraction of lung adenocarcinomas lacks mutations in known oncogenes, and the genesis and treatment of these oncogene-negative tumors remain enigmatic. Here, we perform iterative in vivo functional screens using quantitative autochthonous mouse model systems to uncover the genetic and biochemical changes that enable efficient lung tumor initiation in the absence of oncogene alterations. Generation of hundreds of diverse combinations of tumor suppressor alterations demonstrates that inactivation of suppressors of the RAS and PI3K pathways drives the development of oncogene-negative lung adenocarcinoma. Human genomic data and histology identified RAS/MAPK and PI3K pathway activation as a common feature of an event in oncogene-negative human lung adenocarcinomas. These Onc-negativeRAS/PI3K tumors and related cell lines are vulnerable to pharmacologic inhibition of these signaling axes. These results transform our understanding of this prevalent yet understudied subtype of lung adenocarcinoma.

SIGNIFICANCE

To address the large fraction of lung adenocarcinomas lacking mutations in proto-oncogenes for which targeted therapies are unavailable, this work uncovers driver pathways of oncogene-negative lung adenocarcinomas and demonstrates their therapeutic vulnerabilities.

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.

Tumor suppressor pathways shape EGFR-driven lung tumor progression and response to treatment

In vivo modeling combined with CRISPR/Cas9-mediated somatic genome editing has contributed to elucidating the functional importance of specific genetic alterations in human tumors. Our recent work uncovered tumor suppressor pathways that affect EGFR-driven lung tumor growth and sensitivity to tyrosine kinase inhibitors and reflect the mutational landscape and treatment outcomes in the human disease.

RNA-binding proteins of COSMIC importance in cancer

Herculean efforts by the Wellcome Sanger Institute, the National Cancer Institute, and the National Human Genome Research Institute to sequence thousands of tumors representing all major cancer types have yielded more than 700 genes that contribute to neoplastic growth when mutated, amplified, or deleted. While some of these genes (now included in the COSMIC Cancer Gene Census) encode proteins previously identified in hypothesis-driven experiments (oncogenic transcription factors, protein kinases, etc.), additional classes of cancer drivers have emerged, perhaps none more surprisingly than RNA-binding proteins (RBPs). Over 40 RBPs responsible for virtually all aspects of RNA metabolism, from synthesis to degradation, are recurrently mutated in cancer, and just over a dozen are considered major cancer drivers. This Review investigates whether and how their RNA-binding activities pertain to their oncogenic functions. Focusing on several well-characterized steps in RNA metabolism, we demonstrate that for virtually all cancer-driving RBPs, RNA processing activities are either abolished (the loss-of-function phenotype) or carried out with low fidelity (the LoFi phenotype). Conceptually, this suggests that in normal cells, RBPs act as gatekeepers maintaining proper RNA metabolism and the "balanced" proteome. From the practical standpoint, at least some LoFi phenotypes create therapeutic vulnerabilities, which are beginning to be exploited in the clinic.

Effect of ISM1 on the Immune Microenvironment and Epithelial-Mesenchymal Transition in Colorectal Cancer

Background: An increasing number of studies have shown that Isthmin 1 (ISM1), a secreted protein, is important in tumorigenesis and invasion, including in colorectal cancer (CRC). However, the mechanisms are still unclear. This study aims to explore the function and prognosis capacity of ISM1 in CRC.

Methods: We investigated the expression of ISM1 in 18 CRC tissues vs. adjacent normal tissues from GSE50760, 473 CRC tissues vs. 41 normal tissues from The Cancer Genome Atlas (TCGA), and across gastrointestinal cancer types. Differences were further confirmed in CRC tissues via quantitative real-time polymerase chain reaction (qRT-PCR). Then, we analyzed correlations between clinicopathologic features and ISM1 expression, including prognostic prediction value, using the Kaplan–Meier method and multivariate Cox regression. Gene set enrichment analysis (GSEA) was performed to identify ISM1-related pathways. In vitro experiments were performed to verify the role of ISM1 in epithelial-mesenchymal transition (EMT) and CRC progression.

Results: Multiple datasets showed that ISM1 is upregulated in CRC tissues, which was validated. Patients with higher ISM1 expression had shorter overall survival (OS), and ISM1 expression served as an independent prognostic factor. Enrichment analysis showed that ISM1 upregulation was positively correlated with cancer-related pathways, such as EMT, hypoxia, and the Notch and KRAS signaling pathways. We were exclusively interested in the connection between ISM1 and EMT because 71% of genes in this pathway were significantly positively co-expressed with ISM1, which may account for why patients with higher ISM1 expression are prone to regional lymph node involvement and progression to advanced stages. In addition, we found that ISM1 was positively correlated with multiple immunosuppressive pathways such as IL2/STAT5, TNF-α/NF-κB, and TGF-β, and immune checkpoints, including PD-L1, PD-1, CTLA-4, and LAG3, which may account for upregulation of ISM1 in immunotherapy-resistant patients. Notably, through in vitro experiments, we found that ISM1 promoted EMT and colon cancer cell migration and proliferation.

Conclusion: ISM1 is critical for CRC development and progression, which enhances our understanding of the low response rate of CRC to immunotherapy via immunosuppressive signaling pathways.

Quantitative In Vivo Analyses Reveal a Complex Pharmacogenomic Landscape in Lung Adenocarcinoma

The lack of knowledge about the relationship between tumor genotypes and therapeutic responses remains one of the most critical gaps in enabling the effective use of cancer therapies. Here, we couple a multiplexed and quantitative experimental platform with robust statistical methods to enable pharmacogenomic mapping of lung cancer treatment responses in vivo. The complex map of genotype-specific treatment responses uncovered that over 20% of possible interactions show significant resistance or sensitivity. Known and novel interactions were identified, and one of these interactions, the resistance of KEAP1-mutant lung tumors to platinum therapy, was validated using a large patient response data set. These results highlight the broad impact of tumor suppressor genotype on treatment responses and define a strategy to identify the determinants of precision therapies. SIGNIFICANCE: An experimental and analytical framework to generate in vivo pharmacogenomic maps that relate tumor genotypes to therapeutic responses reveals a surprisingly complex map of genotype-specific resistance and sensitivity.

Identification and Validation of Autophagy-Related Gene Nomograms to Predict the Prognostic Value of Patients with Cervical Cancer

Autophagy is a process of engulfing one's own cytoplasmic proteins or organelles and coating them into vesicles, fusing with lysosomes to form autophagic lysosomes, and degrading the contents it encapsulates. Increasing studies have shown that autophagy disorders are closely related to the occurrence of tumors. However, the prognostic role of autophagy genes in cervical cancer is still unclear. In this study, we constructed risk signatures of autophagy-related genes (ARGs) to predict the prognosis of cervical cancer. The expression profiles and clinical information of autophagy gene sets were downloaded from TCGA and GSE52903 queues as training and validation sets. The normal cervical tissue expression profile data from the UCSC XENA website (obtained from GTEx) were used as a supplement to the TCGA normal cervical tissue. Univariate COX regression analysis of 17 different autophagy genes was performed with the consensus approach. Tumor samples from TCGA were divided into six subtypes, and the clinical traits of the six subtypes had different distributions. Further absolute shrinkage and selection operator (LASSO) and multivariable COX regression yielded an autophagy genetic risk model consisting of eight genes. In the training set, the survival rate of the high-risk group was lower than that of the low-risk group (p < 0.0001). In the validation set, the AUC area of the receiver operating characteristic (ROC) curve was 0.772 for the training set and 0.889 for the verification set. We found that high and low risk scores were closely related to TNM stage (p < 0.05). The nomogram shows that the risk score combined with other indicators, such as G, T, M, and N, better predicts 1-, 3-, and 5-year survival rates. Decline curve analysis (DCA) shows that the risk model combined with other indicators produces better clinical efficacy. Immune cells with an enrichment score of 28 showed statistically significant differences related to high and low risk. GSEA enrichment analysis showed the main enrichment being in KRAS activation, genes defining epithelial and mesenchymal transition (EMT), raised in response to the low oxygen level (hypoxia) gene and NF-kB in response to TNF. These pathways are closely related to the occurrence of tumors. Our constructed autophagy risk signature may be a prognostic tool for cervical cancer.