Comprehensive structure-function analysis reveals gain- and loss-of-function mechanisms impacting oncogenic KRAS activity.

To dissect variant-function relationships in the KRAS oncoprotein, we performed deep mutational scanning (DMS) screens for both wild-type and KRAS G12D mutant alleles. We defined the spectrum of oncogenic potential for nearly all possible KRAS variants, identifying several novel transforming alleles and elucidating a model to describe the frequency of KRAS mutations in human cancer as a function of transforming potential, mutational probability, and tissue-specific mutational signatures. Biochemical and structural analyses of variants identified in a KRAS G12D second-site suppressor DMS screen revealed that attenuation of oncogenic KRAS can be mediated by protein instability and conformational rigidity, resulting in reduced binding affinity to effector proteins, such as RAF and PI3-kinases, or reduced SOS-mediated nucleotide exchange activity. These studies define the landscape of single amino acid alterations that modulate the function of KRAS, providing a resource for the clinical interpretation of KRAS variants and elucidating mechanisms of oncogenic KRAS inactivation for therapeutic exploitation.

YAP1 and PRDM14 converge to promote cell survival and tumorigenesis.

The transcriptional co-activator YAP1 oncogene is the downstream effector of the Hippo pathway, which regulates tissue homeostasis, organ size, regeneration, and tumorigenesis. Multiple cancers are dependent on sustained expression of YAP1 for cell proliferation, survival, and tumorigenesis, but the molecular basis of this oncogene dependency is not well understood. To identify genes that can functionally substitute for YAP1, we performed a genome-scale genetic rescue screen in YAP1-dependent colon cancer cells expressing an inducible YAP1-specific shRNA. We found that the transcription factor PRDM14 rescued cell proliferation and tumorigenesis upon YAP1 suppression in YAP1-dependent cells, xenografts, and colon cancer organoids. YAP1 and PRDM14 individually activated the transcription of calmodulin 2 (CALM2) and a glucose transporter SLC2A1 upon YAP1 suppression, and CALM2 or SLC2A1 expression was required for the rescue of YAP1 suppression. Together, these findings implicate PRDM14-mediated transcriptional upregulation of CALM2 and SLC2A1 as key components of oncogenic YAP1 signaling and dependency.

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.

CREB5 Promotes Resistance to Androgen-Receptor Antagonists and Androgen Deprivation in Prostate Cancer.

Androgen-receptor (AR) inhibitors, including enzalutamide, are used for treatment of all metastatic castration-resistant prostate cancers (mCRPCs). However, some patients develop resistance or never respond. We find that the transcription factor CREB5 confers enzalutamide resistance in an open reading frame (ORF) expression screen and in tumor xenografts. CREB5 overexpression is essential for an enzalutamide-resistant patient-derived organoid. In AR-expressing prostate cancer cells, CREB5 interactions enhance AR activity at a subset of promoters and enhancers upon enzalutamide treatment, including MYC and genes involved in the cell cycle. In mCRPC, we found recurrent amplification and overexpression of CREB5. Our observations identify CREB5 as one mechanism that drives resistance to AR antagonists in prostate cancers.

An alternative splicing switch in FLNB promotes the mesenchymal cell state in human breast cancer.

Alternative splicing of mRNA precursors represents a key gene expression regulatory step and permits the generation of distinct protein products with diverse functions. In a genome-scale expression screen for inducers of the epithelial-to-mesenchymal transition (EMT), we found a striking enrichment of RNA-binding proteins. We validated that QKI and RBFOX1 were necessary and sufficient to induce an intermediate mesenchymal cell state and increased tumorigenicity. Using RNA-seq and eCLIP analysis, we found that QKI and RBFOX1 coordinately regulated the splicing and function of the actin-binding protein FLNB, which plays a causal role in the regulation of EMT. Specifically, the skipping of FLNB exon 30 induced EMT by releasing the FOXC1 transcription factor. Moreover, skipping of FLNB exon 30 is strongly associated with EMT gene signatures in basal-like breast cancer patient samples. These observations identify a specific dysregulation of splicing, which regulates tumor cell plasticity and is frequently observed in human cancer.

Mutational processes shape the landscape of TP53 mutations in human cancer.

Unlike most tumor suppressor genes, the most common genetic alterations in tumor protein p53 (TP53) are missense mutations1,2. Mutant p53 protein is often abundantly expressed in cancers and specific allelic variants exhibit dominant-negative or gain-of-function activities in experimental models3-8. To gain a systematic view of p53 function, we interrogated loss-of-function screens conducted in hundreds of human cancer cell lines and performed TP53 saturation mutagenesis screens in an isogenic pair of TP53 wild-type and null cell lines. We found that loss or dominant-negative inhibition of wild-type p53 function reliably enhanced cellular fitness. By integrating these data with the Catalog of Somatic Mutations in Cancer (COSMIC) mutational signatures database9,10, we developed a statistical model that describes the TP53 mutational spectrum as a function of the baseline probability of acquiring each mutation and the fitness advantage conferred by attenuation of p53 activity. Collectively, these observations show that widely-acting and tissue-specific mutational processes combine with phenotypic selection to dictate the frequencies of recurrent TP53 mutations.