Mutations in "Ras-like in all tissues" (RIT1) occur in up to 2% of lung adenocarcinomas and are mutually exclusive with KRAS and EGFR mutations, suggesting that RIT1 may act as a non-canonical driver oncogene in lung cancer. However, the lack of a RIT1-mutant lung cancer model has hindered the development and testing of RIT1-targeted therapeutics. Here, we report a mouse model with conditional regulation of the cancer-associated RIT1M90I variant. We show that autochthonous expression of RIT1M90I and combined inactivation of Nf2 and p53 drives an aggressive lung cancer with 100% penetrance and short latency. Oncogenic cooperation between RIT1M90I and p53/Nf2 loss is driven by synergistic activation of AP-1 transcription factors and can be reversed by the combined inhibition of MEK and TEAD. These data identify YAP/TEAD as a mediator of RIT1's oncogenic capability and nominate TEAD as a potential drug target in RIT1-mutant lung cancer.

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.

Alternative mRNA splicing is dysregulated in many cancers including lung adenocarcinoma. These aberrant splicing events can sometimes be explained by mutations in splice sites or splicing factors. However, the majority of mis-splicing in cancers remains unexplained. We hypothesize that oncogenic signaling pathways play a role in regulating alternative splicing. In this study, we focus on mRNA splicing activity regulated by the Ras signaling cascade, which is frequently mutated in cancers. These pathways offer avenues for therapeutic modulation through small molecule inhibitors.We took a global proteomic and transcriptomic approach to study the effects of oncogenic Ras signaling on RNA splicing in vitro. Specifically, we performed LC-MS/MS and RNA-seq profiling on normal human lung airway epithelial (AALE) cells overexpressing wild-type or mutant KRAS or RIT1 (N = 3 replicates per allele). Phosphorylation of splicing factors was preferentially downregulated in KRASmut cells compared to KRASWT cells, suggesting that oncogenic KRAS signaling regulates splicing factor activity. Additionally, of the 2227 and 2452 skipped exon (SE) events in KRASG12V and KRASQ61H cells, respectively, compared to KRASWT overexpressing cells, 1013 events were shared between the two mutants. This statistically significant overlap (p < 0.001, hypergeometric test) indicates that KRAS has a role in splicing regulation that is disrupted when the protein is mutated. To determine how KRAS-regulated splicing activity compared to other signaling proteins and oncogenes, we also profiled the whole transcriptomes of isogenic A549 lung adenocarcinoma cell lines overexpressing 86 wild-type or variant alleles across 27 genes implicated in lung cancers (N = 4 to 8 replicates per allele). Of all pairs of mutant and wild-type alleles tested, KRASmut cells compared to KRASWT overexpressing cells exhibited the second highest levels of alternative splicing, behind only RNA binding protein RBM45 and its mutant allele RBM45M126I. In particular, ten SE events were observed to be differentially regulated in a KRAS dependent manner in both AALE and A549 cells. One of these 10 events is the differential splicing of the Myc Associated Zinc Finger (MAZ) protein to increase the expression of its splice variant that acts as a negative regulator of transcription. 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 KRAS mutants. Our proteomic and transcriptomic profiling in lung epithelial and lung adenocarcinoma cells uncovers splicing factor activity and mRNA splicing events regulated by oncogenic KRAS. With further studies of the alternative splicing events described, it will be possible to design or repurpose therapies in order to target aberrant splicing caused by disrupted Ras signaling and its downstream effectors. Citation Format: April Lo, Maria McSharry, Alice H. Berger. Dysregulation of alternative mRNA splicing by oncogenic KRAS in lung adenocarcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 778.

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 KRAS(G12V) and KRAS(Q61H) 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 RBM45(M126I). 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.

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.

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.

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.

Aberrant activation of RAS oncogenes is prevalent in lung adenocarcinoma, with somatic mutation of KRAS occurring in ∼30% of tumors. Recently, we identified somatic mutation of the RAS-family GTPase RIT1 in lung adenocarcinoma, but relatively little is known about the biological pathways regulated by RIT1 and how these relate to the oncogenic KRAS network. Here we present quantitative proteomic and transcriptomic profiles from KRAS -mutant and RIT1 -mutant isogenic lung epithelial cells and globally characterize the signaling networks regulated by each oncogene. We find that both mutant KRAS and mutant RIT1 promote S6 kinase, AKT, and RAF/MEK signaling, and promote epithelial-to-mesenchymal transition and immune evasion via HLA protein loss. However, KRAS and RIT1 diverge in regulation of phosphorylation sites on EGFR, USO1, and AHNAK proteins. The majority of the proteome changes are related to altered transcriptional regulation, but a small subset of proteins are differentially regulated by both oncoproteins at the post-transcriptional level, including intermediate filament proteins, metallothioneins, and MHC Class I proteins. These data provide the first global, unbiased characterization of oncogenic RIT1 network and identify the shared and divergent functions of oncogenic RIT1 and KRAS GTPases in lung cancer.

The KPC mouse model, driven by the Kras and Trp53 transgenes, is well regarded for faithful recapitulation of human pancreatic cancer biology. However, the extent that this model recapitulates the subclonal evolution of this tumor type is unknown. Here we report evidence of continuing subclonal evolution after tumor initiation that largely reflect copy number alterations that target cellular processes of established significance in human pancreatic cancer. The evolutionary trajectories of the mouse tumors show both linear and branching patterns as well as clonal mixing. We propose the KPC model and derivatives have unexplored utility as a functional system to model the mechanisms and modifiers of tumor evolution.

RNA splicing is dysregulated in a widespread manner in cancers including lung adenocarcinoma. In some cases, splicing changes can be attributed to cis-acting splice site mutations or trans-acting mutations in splicing factors. However, in most cases, the underlying causes of splicing changes are unknown. We hypothesized that upstream signaling inputs to alternative splicing regulation can explain some of these unknowns. Specifically, we studied how oncogenic lung cancer signaling pathways (EGFR/Ras, KEAP1/NRF2, MYC, and others) regulate alternative splicing and the expression and activity of splicing factors. We focus on signaling pathways because they can be readily therapeutically modulated with small molecule inhibitors (e.g. tyrosine kinase inhibitors), offering opportunities for therapeutic suppression of downstream splicing effects. To experimentally determine how signaling pathway perturbation affects alternative splicing, we perturbed A549 lung cancer cells with each of 82 alleles of 27 genes (n = 4 to 8 biological replicates per allele). Genes and variants were selected based on the occurrence of the variants in lung adenocarcinoma tumors. In total, 417 whole transcriptome profiles were generated using the Smart-Seq v4 method (Clontech) and Nextera XT library preparation (Illumina). We performed differential expression analysis using edgeR and differential splicing analysis using MISO. Using this approach, we identified a high-confidence set of 2430 alternative splicing events differentially spliced in the perturbed samples compared to controls. Of these alternative splicing events, 1219 are skipped exons, 469 are mutually exclusive exons, 235 are alternative 5’ splice sites, 255 are alternative 3’ splice sites, and 252 are retained introns. Among the perturbations tested, overexpression of the RBM45 wild-type allele and the RBM45 D434Y variant allele resulted in the greatest number of alternatively spliced events, with 367 and 323 events respectively. A closer look at these events reveals that when the RBM45 wild-type allele is overexpressed, exon 3 of the cyclin gene CCNG1 is skipped. Interestingly, when RBM45 variant alleles (D434Y, M126I) are overexpressed instead of the wild-type, the exon is not skipped, suggesting loss or change of function. These results suggest that, among other roles, RBM45 may regulate cell cycle patterns by regulating alternative splicing. In sum, our screen identifies splicing events which are regulated by oncogenic signaling pathways. With this information, it will be possible to propose therapeutic options that mitigate aberrant splicing driven by signaling pathway components. Such therapy may already exist in the form of drugs for other targeted uses but could be repurposed to address aberrant splicing. Importantly, therapy of this type has the potential to be more tolerable compared to treatments directly aimed at splicing factors and splice sites. Citation Format: April Lo, Maria McSharry, Alice Berger. A large-scale RNA-seq screen to identify regulators of alternative splicing in cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 4345.