The genomic landscape of lung cancer in never-smokers from the Women's Health Initiative.

Over 200,000 individuals are diagnosed with lung cancer in the U.S. every year, with a growing proportion of cases, especially lung adenocarcinoma, occurring in individuals who have never smoked. Women over the age of 50 comprise the largest affected demographic. To understand the genomic drivers of lung adenocarcinoma and therapeutic response in this population, we performed whole genome and/or whole exome sequencing on 73 matched lung tumor/normal pairs from post-menopausal women who participated in the Women's Health Initiative. Somatic copy number alterations showed little variation by smoking status, suggesting that aneuploidy may be a general characteristic of lung cancer regardless of smoke exposure. Similarly, clock-like and APOBEC mutation signatures were prevalent but did not differ in tumors from smokers and never-smokers. However, mutations in both EGFR and KRAS showed unique allelic differences determined by smoking status that are known to alter tumor response to targeted therapy. Mutations in the MYC-network member MGA were more prevalent in tumors from smokers. Fusion events in ALK, RET, and ROS1 were absent, likely due to age-related differences in fusion prevalence. Our work underscores the profound impact of smoking status, age, and sex on the tumor mutational landscape and identifies areas of unmet medical need.

The deubiquitinase USP9X regulates RIT1 protein abundance and oncogenic phenotypes.

RIT1 is a rare and understudied oncogene in lung cancer. Despite structural similarity to other RAS GTPase proteins such as KRAS, oncogenic RIT1 activity does not appear to be tightly regulated by nucleotide exchange or hydrolysis. Instead, there is a growing understanding that the protein abundance of RIT1 is important for its regulation and function. We previously identified the deubiquitinase USP9X as a RIT1 dependency in RIT1-mutant cells. Here, we demonstrate that both wild-type and mutant forms of RIT1 are substrates of USP9X. Depletion of USP9X leads to decreased RIT1 protein stability and abundance and resensitizes cells to EGFR tyrosine kinase inhibitors. Our work expands upon the current understanding of RIT1 protein regulation and presents USP9X as a key regulator of RIT1-driven oncogenic phenotypes.

Chromosome-level genome editing uncovers aneuploidy addiction in cancer.

Publishing in Science, Girish and colleagues achieve chromosome-level genome editing to reveal a requirement for aneuploidy in breast and melanoma cancers. Authors developed and leveraged ReDACT (restoring disomy in aneuploid cells using CRISPR targeting) to generate isogenic models of aneuploidy and demonstrate that some cancers are addicted to increased copy number of specific chromosome arms.

pgMAP: a pipeline to enable guide RNA read mapping from dual-targeting CRISPR screens.

We developed pgMAP, an analysis pipeline to map gRNA sequencing reads from dual-targeting CRISPR screens. pgMAP output includes a dual gRNA read counts table and quality control metrics including the proportion of correctly-paired reads and CRISPR library sequencing coverage across all time points and samples. pgMAP is implemented using Snakemake and is available open-source under the MIT license at https://github.com/fredhutch/pgmap_pipeline.

Somatic mutation but not aneuploidy differentiates lung cancer in never-smokers and smokers.

Lung cancer in never-smokers disproportionately affects older women. To understand the mutational landscape of this cohort, we performed detailed genome characterization of 73 lung adenocarcinomas from participants of the Women’s Health Initiative (WHI). We find enrichment of EGFR mutations in never-/light-smokers and KRAS mutations in heavy smokers as expected, but we also show that the specific variants of these genes differ by smoking status, with important therapeutic implications. Mutational signature analysis revealed signatures of clock, APOBEC, and DNA repair deficiency in never-/light-smokers; however, the mutational load of these signatures did not differ significantly from those found in smokers. Last, tumors from both smokers and never-/light-smokers shared copy number subtypes, with no significant differences in aneuploidy. Thus, the genomic landscape of lung cancer in never-/light-smokers and smokers is predominantly differentiated by somatic mutations and not copy number alterations.

Oncogenic KRAS alters splicing factor phosphorylation and alternative splicing in lung cancer.

BACKGROUND: Alternative RNA splicing is widely dysregulated in cancers including lung adenocarcinoma, where aberrant splicing events are frequently caused by somatic splice site mutations or somatic mutations of splicing factor genes. However, the majority of mis-splicing in cancers is unexplained by these known mechanisms. We hypothesize that the aberrant Ras signaling characteristic of lung cancers plays a role in promoting the alternative splicing observed in tumors. METHODS: We recently performed transcriptome and proteome profiling of human lung epithelial cells ectopically expressing oncogenic KRAS and another cancer-associated Ras GTPase, RIT1. Unbiased analysis of phosphoproteome data identified altered splicing factor phosphorylation in KRAS-mutant cells, so we performed differential alternative splicing analysis using rMATS to identify significantly altered isoforms in lung epithelial cells. To determine whether these isoforms were uniquely regulated by KRAS, we performed a large-scale splicing screen in which we generated over 300 unique RNA sequencing profiles of isogenic A549 lung adenocarcinoma cells ectopically expressing 75 different wild-type or variant alleles across 28 genes implicated in lung cancer. RESULTS: Mass spectrometry data showed widespread downregulation of splicing factor phosphorylation in lung epithelial cells expressing mutant KRAS compared to cells expressing wild-type KRAS. We observed alternative splicing in the same cells, with 2196 and 2416 skipped exon events in KRASG12V and KRASQ61H cells, respectively, 997 of which were shared (p < 0.001 by hypergeometric test). In the high-throughput splicing screen, mutant KRAS induced the greatest number of differential alternative splicing events, second only to the RNA binding protein RBM45 and its variant RBM45M126I. We identified ten high confidence cassette exon events across multiple KRAS variants and cell lines. These included differential splicing of the Myc Associated Zinc Finger (MAZ). As MAZ regulates expression of KRAS, this splice variant may be a mechanism for the cell to modulate wild-type KRAS levels in the presence of oncogenic KRAS. CONCLUSION: Proteomic and transcriptomic profiling of lung epithelial cells uncovered splicing factor phosphorylation and mRNA splicing events regulated by oncogenic KRAS. These data suggest that in addition to widespread transcriptional changes, the Ras signaling pathway can promote post-transcriptional splicing changes that may contribute to oncogenic processes.

All About That Ras: Novel Fusion Drives Ras Pathway Activation in Lung Cancer.

Lung cancers in never- and light-smokers often harbor targetable oncogenic mutations in Ras pathway genes. Here, a novel OCLN-RASGRF1 fusion is identified in an otherwise Ras wild-type lung tumor. Studying this and other RASGRF1 fusions, the authors show that these fusions lead to malignant phenotypes that can be reversed by MEK inhibition. See related article by Hunihan et al., p. 3091.

Cell Painting predicts impact of lung cancer variants.

Most variants in most genes across most organisms have an unknown impact on the function of the corresponding gene. This gap in knowledge is especially acute in cancer, where clinical sequencing of tumors now routinely reveals patient-specific variants whose functional impact on the corresponding genes is unknown, impeding clinical utility. Transcriptional profiling was able to systematically distinguish these variants of unknown significance as impactful vs. neutral in an approach called expression-based variant-impact phenotyping. We profiled a set of lung adenocarcinoma-associated somatic variants using Cell Painting, a morphological profiling assay that captures features of cells based on microscopy using six stains of cell and organelle components. Using deep-learning-extracted features from each cell's image, we found that cell morphological profiling (cmVIP) can predict variants' functional impact and, particularly at the single-cell level, reveals biological insights into variants that can be explored at our public online portal. Given its low cost, convenient implementation, and single-cell resolution, cmVIP profiling therefore seems promising as an avenue for using non-gene specific assays to systematically assess the impact of variants, including disease-associated alleles, on gene function.

Massively parallel phenotyping of coding variants in cancer with Perturb-seq.

Genome sequencing studies have identified millions of somatic variants in cancer, but it remains challenging to predict the phenotypic impact of most. Experimental approaches to distinguish impactful variants often use phenotypic assays that report on predefined gene-specific functional effects in bulk cell populations. Here, we develop an approach to functionally assess variant impact in single cells by pooled Perturb-seq. We measured the impact of 200 TP53 and KRAS variants on RNA profiles in over 300,000 single lung cancer cells, and used the profiles to categorize variants into phenotypic subsets to distinguish gain-of-function, loss-of-function and dominant negative variants, which we validated by comparison with orthogonal assays. We discovered that KRAS variants did not merely fit into discrete functional categories, but spanned a continuum of gain-of-function phenotypes, and that their functional impact could not have been predicted solely by their frequency in patient cohorts. Our work provides a scalable, gene-agnostic method for coding variant impact phenotyping, with potential applications in multiple disease settings.

Multiomic characterization of oncogenic signaling mediated by wild-type and mutant RIT1.

Aberrant activation of the RAS family of guanosine triphosphatases (GTPases) is prevalent in lung adenocarcinoma, with somatic mutation of KRAS occurring in ~30% of tumors. We previously identified somatic mutations and amplifications of the gene encoding RAS family GTPase RIT1 in lung adenocarcinomas. To explore the biological pathways regulated by RIT1 and how they relate to the oncogenic KRAS network, we performed quantitative proteomic, phosphoproteomic, and transcriptomic profiling of isogenic lung epithelial cells in which we ectopically expressed wild-type or cancer-associated variants of RIT1 and KRAS. We found that both mutant KRAS and mutant RIT1 promoted canonical RAS signaling and that overexpression of wild-type RIT1 partially phenocopied oncogenic RIT1 and KRAS, including induction of epithelial-to-mesenchymal transition. Our findings suggest that RIT1 protein abundance is a factor in its pathogenic function. Therefore, chromosomal amplification of wild-type RIT1 in lung and other cancers may be tumorigenic.

Discovery of synthetic lethal and tumor suppressor paralog pairs in the human genome.

CRISPR screens have accelerated the discovery of important cancer vulnerabilities. However, single-gene knockout phenotypes can be masked by redundancy among related genes. Paralogs constitute two-thirds of the human protein-coding genome, so existing methods are likely inadequate for assaying a large portion of gene function. Here, we develop paired guide RNAs for paralog genetic interaction mapping (pgPEN), a pooled CRISPR-Cas9 single- and double-knockout approach targeting more than 2,000 human paralogs. We apply pgPEN to two cell types and discover that 12% of human paralogs exhibit synthetic lethality in at least one context. We recover known synthetic lethal paralogs MEK1/MEK2, important drug targets CDK4/CDK6, and other synthetic lethal pairs including CCNL1/CCNL2. Additionally, we identify ten tumor suppressor paralog pairs whose compound loss promotes cell proliferation. These findings nominate drug targets and suggest that paralog genetic interactions could shape the landscape of positive and negative selection in cancer.

Integrative oncogene-dependency mapping identifies RIT1 vulnerabilities and synergies in lung cancer.

CRISPR-based cancer dependency maps are accelerating advances in cancer precision medicine, but adequate functional maps are limited to the most common oncogenes. To identify opportunities for therapeutic intervention in other rarer subsets of cancer, we investigate the oncogene-specific dependencies conferred by the lung cancer oncogene, RIT1. Here, genome-wide CRISPR screening in KRAS, EGFR, and RIT1-mutant isogenic lung cancer cells identifies shared and unique vulnerabilities of each oncogene. Combining this genetic data with small-molecule sensitivity profiling, we identify a unique vulnerability of RIT1-mutant cells to loss of spindle assembly checkpoint regulators. Oncogenic RIT1M90I weakens the spindle assembly checkpoint and perturbs mitotic timing, resulting in sensitivity to Aurora A inhibition. In addition, we observe synergy between mutant RIT1 and activation of YAP1 in multiple models and frequent nuclear overexpression of YAP1 in human primary RIT1-mutant lung tumors. These results provide a genome-wide atlas of oncogenic RIT1 functional interactions and identify components of the RAS pathway, spindle assembly checkpoint, and Hippo/YAP1 network as candidate therapeutic targets in RIT1-mutant lung cancer.

eVIP2: Expression-based variant impact phenotyping to predict the function of gene variants.

While advancements in genome sequencing have identified millions of somatic mutations in cancer, their functional impact is poorly understood. We previously developed the expression-based variant impact phenotyping (eVIP) method to use gene expression data to characterize the function of gene variants. The eVIP method uses a decision tree-based algorithm to predict the functional impact of somatic variants by comparing gene expression signatures induced by introduction of wild-type (WT) versus mutant cDNAs in cell lines. The method distinguishes between variants that are gain-of-function, loss-of-function, change-of-function, or neutral. We present eVIP2, software that allows for pathway analysis (eVIP Pathways) and usage with RNA-seq data. To demonstrate the eVIP2 software and approach, we characterized two recurrent frameshift variants in RNF43, a negative regulator of Wnt signaling, frequently mutated in colorectal, gastric, and endometrial cancer. RNF43 WT, RNF43 R117fs, RNF43 G659fs, or GFP control cDNA were overexpressed in HEK293T cells. Analysis with eVIP2 predicted that the frameshift at position 117 was a loss-of-function mutation, as expected. The second frameshift at position 659 has been previously described as a passenger mutation that maintains the RNF43 WT function as a negative regulator of Wnt. Surprisingly, eVIP2 predicted G659fs to be a change-of-function mutation. Additional eVIP Pathways analysis of RNF43 G659fs predicted 10 pathways to be significantly altered, including TNF-α via NFκB signaling, KRAS signaling, and hypoxia, highlighting the benefit of a more comprehensive approach when determining the impact of gene variant function. To validate these predictions, we performed reporter assays and found that each pathway activated by expression of RNF43 G659fs, but not expression of RNF43 WT, was identified as impacted by eVIP2, supporting that RNF43 G659fs is a change-of-function mutation and its effect on the identified pathways. Pathway activation was further validated by Western blot analysis. Lastly, we show primary colon adenocarcinoma patient samples with R117fs and G659fs variants have transcriptional profiles similar to BRAF missense mutations with activated RAS/MAPK signaling, consistent with KRAS signaling pathways being GOF in both variants. The eVIP2 method is an important step towards overcoming the current challenge of variant interpretation in the implementation of precision medicine. eVIP2 is available at https://github.com/BrooksLabUCSC/eVIP2.

Loss of MGA repression mediated by an atypical polycomb complex promotes tumor progression and invasiveness.

MGA, a transcription factor and member of the MYC network, is mutated or deleted in a broad spectrum of malignancies. As a critical test of a tumor suppressive role, we inactivated Mga in two mouse models of non-small cell lung cancer using a CRISPR-based approach. MGA loss significantly accelerated tumor growth in both models and led to de-repression of non-canonical Polycomb ncPRC1.6 targets, including genes involved in metastasis and meiosis. Moreover, MGA deletion in human lung adenocarcinoma lines augmented invasive capabilities. We further show that MGA-MAX, E2F6, and L3MBTL2 co-occupy thousands of promoters and that MGA stabilizes these ncPRC1.6 subunits. Lastly, we report that MGA loss also induces a pro-growth effect in human colon organoids. Our studies establish MGA as a bona fide tumor suppressor in vivo and suggest a tumor suppressive mechanism in adenocarcinomas resulting from widespread transcriptional attenuation of MYC and E2F target genes mediated by MGA-MAX associated with a non-canonical Polycomb complex.

Whole-genome characterization of lung adenocarcinomas lacking the RTK/RAS/RAF pathway.

RTK/RAS/RAF pathway alterations (RPAs) are a hallmark of lung adenocarcinoma (LUAD). In this study, we use whole-genome sequencing (WGS) of 85 cases found to be RPA(-) by previous studies from The Cancer Genome Atlas (TCGA) to characterize the minority of LUADs lacking apparent alterations in this pathway. We show that WGS analysis uncovers RPA(+) in 28 (33%) of the 85 samples. Among the remaining 57 cases, we observe focal deletions targeting the promoter or transcription start site of STK11 (n = 7) or KEAP1 (n = 3), and promoter mutations associated with the increased expression of ILF2 (n = 6). We also identify complex structural variations associated with high-level copy number amplifications. Moreover, an enrichment of focal deletions is found in TP53 mutant cases. Our results indicate that RPA(-) cases demonstrate tumor suppressor deletions and genome instability, but lack unique or recurrent genetic lesions compensating for the lack of RPAs. Larger WGS studies of RPA(-) cases are required to understand this important LUAD subset.

CRISPR base editor screens identify variant function at scale.

Cuella-Martin et al. (2021) and Hanna et al. (2021) showcase CRISPR base editing in large-scale pooled screens in human cells to discover both loss- and gain-of-function variants, enabling protein structure/function insights and clinical variant interpretation.

Genome-wide CRISPR screening reveals novel therapeutic targets in RIT1-driven lung cancer.

In recent work, we performed CRISPR/Cas9 screening in RIT1 (Ras-like in all tissues)-mutant cancer cells. We found that RIT1-mutant cells are vulnerable to loss of mitotic regulators, and mutant RIT1 synergizes with YAP1 (yes-associated protein 1) in oncogenesis. These findings can be leveraged to identify targeted therapies for RIT1-mutant cancer.

RNA isoform screens uncover the essentiality and tumor-suppressor activity of ultraconserved poison exons.

While RNA-seq has enabled comprehensive quantification of alternative splicing, no correspondingly high-throughput assay exists for functionally interrogating individual isoforms. We describe pgFARM (paired guide RNAs for alternative exon removal), a CRISPR-Cas9-based method to manipulate isoforms independent of gene inactivation. This approach enabled rapid suppression of exon recognition in polyclonal settings to identify functional roles for individual exons, such as an SMNDC1 cassette exon that regulates pan-cancer intron retention. We generalized this method to a pooled screen to measure the functional relevance of 'poison' cassette exons, which disrupt their host genes' reading frames yet are frequently ultraconserved. Many poison exons were essential for the growth of both cultured cells and lung adenocarcinoma xenografts, while a subset had clinically relevant tumor-suppressor activity. The essentiality and cancer relevance of poison exons are likely to contribute to their unusually high conservation and contrast with the dispensability of other ultraconserved elements for viability.