Systematic discovery of mutation-directed neo-protein-protein interactions in cancer.

Comprehensive sequencing of patient tumors reveals genomic mutations across tumor types that enable tumorigenesis and progression. A subset of oncogenic driver mutations results in neomorphic activity where the mutant protein mediates functions not engaged by the parental molecule. Here, we identify prevalent variant-enabled neomorph-protein-protein interactions (neoPPI) with a quantitative high-throughput differential screening (qHT-dS) platform. The coupling of highly sensitive BRET biosensors with miniaturized coexpression in an ultra-HTS format allows large-scale monitoring of the interactions of wild-type and mutant variant counterparts with a library of cancer-associated proteins in live cells. The screening of 17,792 interactions with 2,172,864 data points revealed a landscape of gain of interactions encompassing both oncogenic and tumor suppressor mutations. For example, the recurrent BRAF V600E lesion mediates KEAP1 neoPPI, rewiring a BRAFV600E/KEAP1 signaling axis and creating collateral vulnerability to NQO1 substrates, offering a combination therapeutic strategy. Thus, cancer genomic alterations can create neo-interactions, informing variant-directed therapeutic approaches for precision medicine.

Structural basis for RNA surveillance by the human nuclear exosome targeting (NEXT) complex.

RNA quality control relies on co-factors and adaptors to identify and prepare substrates for degradation by ribonucleases such as the 3' to 5' ribonucleolytic RNA exosome. Here, we determined cryogenic electron microscopy structures of human nuclear exosome targeting (NEXT) complexes bound to RNA that reveal mechanistic insights to substrate recognition and early steps that precede RNA handover to the exosome. The structures illuminate ZCCHC8 as a scaffold, mediating homodimerization while embracing the MTR4 helicase and flexibly anchoring RBM7 to the helicase core. All three subunits collaborate to bind the RNA, with RBM7 and ZCCHC8 surveying sequences upstream of the 3' end to facilitate RNA capture by MTR4. ZCCHC8 obscures MTR4 surfaces important for RNA binding and extrusion as well as MPP6-dependent recruitment and docking onto the RNA exosome core, interactions that contribute to RNA surveillance by coordinating RNA capture, translocation, and extrusion from the helicase to the exosome for decay.

Humans with inherited T cell CD28 deficiency are susceptible to skin papillomaviruses but are otherwise healthy.

We study a patient with the human papilloma virus (HPV)-2-driven "tree-man" phenotype and two relatives with unusually severe HPV4-driven warts. The giant horns form an HPV-2-driven multifocal benign epithelial tumor overexpressing viral oncogenes in the epidermis basal layer. The patients are unexpectedly homozygous for a private CD28 variant. They have no detectable CD28 on their T cells, with the exception of a small contingent of revertant memory CD4+ T cells. T cell development is barely affected, and T cells respond to CD3 and CD2, but not CD28, costimulation. Although the patients do not display HPV-2- and HPV-4-reactive CD4+ T cells in vitro, they make antibodies specific for both viruses in vivo. CD28-deficient mice are susceptible to cutaneous infections with the mouse papillomavirus MmuPV1. The control of HPV-2 and HPV-4 in keratinocytes is dependent on the T cell CD28 co-activation pathway. Surprisingly, human CD28-dependent T cell responses are largely redundant for protective immunity.

Revealing the Activity of Trimeric G-proteins in Live Cells with a Versatile Biosensor Design.

Heterotrimeric G-proteins (Gαßγ) are the main transducers of signals from GPCRs, mediating the action of countless natural stimuli and therapeutic agents. However, there are currently no robust approaches to directly measure the activity of endogenous G-proteins in cells. Here, we describe a suite of optical biosensors that detect endogenous active G-proteins with sub-second resolution in live cells. Using a modular design principle, we developed genetically encoded, unimolecular biosensors for endogenous Gα-GTP and free Gßγ: the two active species of heterotrimeric G-proteins. This design was leveraged to generate biosensors with specificity for different heterotrimeric G-proteins or for other G-proteins, such as Rho GTPases. Versatility was further validated by implementing the biosensors in multiple contexts, from characterizing cancer-associated G-protein mutants to neurotransmitter signaling in primary neurons. Overall, the versatile biosensor design introduced here enables studying the activity of endogenous G-proteins in live cells with high fidelity, temporal resolution, and convenience.

Functional Enhancers Shape Extrachromosomal Oncogene Amplifications.

Non-coding regions amplified beyond oncogene borders have largely been ignored. Using a computational approach, we find signatures of significant co-amplification of non-coding DNA beyond the boundaries of amplified oncogenes across five cancer types. In glioblastoma, EGFR is preferentially co-amplified with its two endogenous enhancer elements active in the cell type of origin. These regulatory elements, their contacts, and their contribution to cell fitness are preserved on high-level circular extrachromosomal DNA amplifications. Interrogating the locus with a CRISPR interference screening approach reveals a diversity of additional elements that impact cell fitness. The pattern of fitness dependencies mirrors the rearrangement of regulatory elements and accompanying rewiring of the chromatin topology on the extrachromosomal amplicon. Our studies indicate that oncogene amplifications are shaped by regulatory dependencies in the non-coding genome.

Tracing Oncogene Rearrangements in the Mutational History of Lung Adenocarcinoma.

Mutational processes giving rise to lung adenocarcinomas (LADCs) in non-smokers remain elusive. We analyzed 138 LADC whole genomes, including 83 cases with minimal contribution of smoking-associated mutational signature. Genomic rearrangements were not correlated with smoking-associated mutations and frequently served as driver events of smoking-signature-low LADCs. Complex genomic rearrangements, including chromothripsis and chromoplexy, generated 74% of known fusion oncogenes, including EML4-ALK, CD74-ROS1, and KIF5B-RET. Unlike other collateral rearrangements, these fusion-oncogene-associated rearrangements were frequently copy-number-balanced, representing a genomic signature of early oncogenesis. Analysis of mutation timing revealed that fusions and point mutations of canonical oncogenes were often acquired in the early decades of life. During a long latency, cancer-related genes were disrupted or amplified by complex rearrangements. The genomic landscape was different between subgroups-EGFR-mutant LADCs had frequent whole-genome duplications with p53 mutations, whereas fusion-oncogene-driven LADCs had frequent SETD2 mutations. Our study highlights LADC oncogenesis driven by endogenous mutational processes.

Structural Basis of Smoothened Activation in Hedgehog Signaling.

The seven-transmembrane-spanning protein Smoothened is the central transducer in Hedgehog signaling, a pathway fundamental in development and in cancer. Smoothened is activated by cholesterol binding to its extracellular cysteine-rich domain (CRD). How this interaction leads to changes in the transmembrane domain and Smoothened activation is unknown. Here, we report crystal structures of sterol-activated Smoothened. The CRD undergoes a dramatic reorientation, allosterically causing the transmembrane domain to adopt a conformation similar to active G-protein-coupled receptors. We show that Smoothened contains a unique inhibitory π-cation lock, which is broken on activation and is disrupted in constitutively active oncogenic mutants. Smoothened activation opens a hydrophobic tunnel, suggesting a pathway for cholesterol movement from the inner membrane leaflet to the CRD. All Smoothened antagonists bind the transmembrane domain and block tunnel opening, but cyclopamine also binds the CRD, inducing the active transmembrane conformation. Together, these results define the mechanisms of Smoothened activation and inhibition.

KRAS Dimerization Impacts MEK Inhibitor Sensitivity and Oncogenic Activity of Mutant KRAS.

The mechanism by which the wild-type KRAS allele imparts a growth inhibitory effect to oncogenic KRAS in various cancers, including lung adenocarcinoma (LUAD), is poorly understood. Here, using a genetically inducible model of KRAS loss of heterozygosity (LOH), we show that KRAS dimerization mediates wild-type KRAS-dependent fitness of human and murine KRAS mutant LUAD tumor cells and underlies resistance to MEK inhibition. These effects are abrogated when wild-type KRAS is replaced by KRASD154Q, a mutant that disrupts dimerization at the α4-α5 KRAS dimer interface without changing other fundamental biochemical properties of KRAS, both in vitro and in vivo. Moreover, dimerization has a critical role in the oncogenic activity of mutant KRAS. Our studies provide mechanistic and biological insights into the role of KRAS dimerization and highlight a role for disruption of dimerization as a therapeutic strategy for KRAS mutant cancers.

Targeting KRAS Mutant Cancers with a Covalent G12C-Specific Inhibitor.

KRASG12C was recently identified to be potentially druggable by allele-specific covalent targeting of Cys-12 in vicinity to an inducible allosteric switch II pocket (S-IIP). Success of this approach requires active cycling of KRASG12C between its active-GTP and inactive-GDP conformations as accessibility of the S-IIP is restricted only to the GDP-bound state. This strategy proved feasible for inhibiting mutant KRAS in vitro; however, it is uncertain whether this approach would translate to in vivo. Here, we describe structure-based design and identification of ARS-1620, a covalent compound with high potency and selectivity for KRASG12C. ARS-1620 achieves rapid and sustained in vivo target occupancy to induce tumor regression. We use ARS-1620 to dissect oncogenic KRAS dependency and demonstrate that monolayer culture formats significantly underestimate KRAS dependency in vivo. This study provides in vivo evidence that mutant KRAS can be selectively targeted and reveals ARS-1620 as representing a new generation of KRASG12C-specific inhibitors with promising therapeutic potential.

Myc Cooperates with Ras by Programming Inflammation and Immune Suppression.

The two oncogenes KRas and Myc cooperate to drive tumorigenesis, but the mechanism underlying this remains unclear. In a mouse lung model of KRasG12D-driven adenomas, we find that co-activation of Myc drives the immediate transition to highly proliferative and invasive adenocarcinomas marked by highly inflammatory, angiogenic, and immune-suppressed stroma. We identify epithelial-derived signaling molecules CCL9 and IL-23 as the principal instructing signals for stromal reprogramming. CCL9 mediates recruitment of macrophages, angiogenesis, and PD-L1-dependent expulsion of T and B cells. IL-23 orchestrates exclusion of adaptive T and B cells and innate immune NK cells. Co-blockade of both CCL9 and IL-23 abrogates Myc-induced tumor progression. Subsequent deactivation of Myc in established adenocarcinomas triggers immediate reversal of all stromal changes and tumor regression, which are independent of CD4+CD8+ T cells but substantially dependent on returning NK cells. We show that Myc extensively programs an immune suppressive stroma that is obligatory for tumor progression.

Kinome-wide siRNA screen identifies a DCLK2-TBK1 oncogenic signaling axis in clear cell renal cell carcinoma.

TANK-binding kinase 1 (TBK1) is a potential therapeutic target in multiple cancers, including clear cell renal cell carcinoma (ccRCC). However, targeting TBK1 in clinical practice is challenging. One approach to overcome this challenge would be to identify an upstream TBK1 regulator that could be targeted therapeutically in cancer specifically. In this study, we perform a kinome-wide small interfering RNA (siRNA) screen and identify doublecortin-like kinase 2 (DCLK2) as a TBK1 regulator in ccRCC. DCLK2 binds to and directly phosphorylates TBK1 on Ser172. Depletion of DCLK2 inhibits anchorage-independent colony growth and kidney tumorigenesis in orthotopic xenograft models. Conversely, overexpression of DCLK2203, a short isoform that predominates in ccRCC, promotes ccRCC cell growth and tumorigenesis in vivo. Mechanistically, DCLK2203 elicits its oncogenic signaling via TBK1 phosphorylation and activation. Taken together, these results suggest that DCLK2 is a TBK1 activator and potential therapeutic target for ccRCC.

Higher-order SPOP assembly reveals a basis for cancer mutant dysregulation.

The speckle-type POZ protein (SPOP) functions in the Cullin3-RING ubiquitin ligase (CRL3) as a receptor for the recognition of substrates involved in cell growth, survival, and signaling. SPOP mutations have been attributed to the development of many types of cancers, including prostate and endometrial cancers. Prostate cancer mutations localize in the substrate-binding site of the substrate recognition (MATH) domain and reduce or prevent binding. However, most endometrial cancer mutations are dispersed in seemingly inconspicuous solvent-exposed regions of SPOP, offering no clear basis for their cancer-causing and peculiar gain-of-function properties. Herein, we present the first structure of SPOP in its oligomeric form, uncovering several new interfaces important for SPOP self-assembly and normal function. Given that many previously unaccounted-for cancer mutations are localized in these newly identified interfaces, we uncover molecular mechanisms underlying dysregulation of SPOP function, with effects ranging from gross structural changes to enhanced self-association, and heightened stability and activity.

RAD51 protects human cells from transcription-replication conflicts.

Oncogene activation during tumorigenesis promotes DNA replication stress (RS), which subsequently drives the formation of cancer-associated chromosomal rearrangements. Many episodes of physiological RS likely arise due to conflicts between the DNA replication and transcription machineries operating simultaneously at the same loci. One role of the RAD51 recombinase in human cells is to protect replication forks undergoing RS. Here, we have identified a key role for RAD51 in preventing transcription-replication conflicts (TRCs) from triggering replication fork breakage. The genomic regions most affected by RAD51 deficiency are characterized by being replicated and transcribed in early S-phase and show significant overlap with loci prone to cancer-associated amplification. Consistent with a role for RAD51 in protecting against transcription-replication conflicts, many of the adverse effects of RAD51 depletion are ameliorated by inhibiting early S-phase transcription. We propose a model whereby RAD51 suppresses fork breakage and subsequent inadvertent amplification of genomic loci prone to experiencing TRCs.

METTL1-mediated m7G modification of Arg-TCT tRNA drives oncogenic transformation.

The emerging "epitranscriptomics" field is providing insights into the biological and pathological roles of different RNA modifications. The RNA methyltransferase METTL1 catalyzes N7-methylguanosine (m7G) modification of tRNAs. Here we find METTL1 is frequently amplified and overexpressed in cancers and is associated with poor patient survival. METTL1 depletion causes decreased abundance of m7G-modified tRNAs and altered cell cycle and inhibits oncogenicity. Conversely, METTL1 overexpression induces oncogenic cell transformation and cancer. Mechanistically, we find increased abundance of m7G-modified tRNAs, in particular Arg-TCT-4-1, and increased translation of mRNAs, including cell cycle regulators that are enriched in the corresponding AGA codon. Accordingly, Arg-TCT expression is elevated in many tumor types and is associated with patient survival, and strikingly, overexpression of this individual tRNA induces oncogenic transformation. Thus, METTL1-mediated tRNA modification drives oncogenic transformation through a remodeling of the mRNA "translatome" to increase expression of growth-promoting proteins and represents a promising anti-cancer target.

Clonal selection of stable aneuploidies in progenitor cells drives high-prevalence tumorigenesis.

Chromosome gains and losses are a frequent feature of human cancers. However, how these aberrations can outweigh the detrimental effects of aneuploidy remains unclear. An initial comparison of existing chromosomal instability (CIN) mouse models suggests that aneuploidy accumulates to low levels in these animals. We therefore developed a novel mouse model that enables unprecedented levels of chromosome missegregation in the adult animal. At the earliest stages of T-cell development, cells with random chromosome gains and/or losses are selected against, but CIN eventually results in the expansion of progenitors with clonal chromosomal imbalances. Clonal selection leads to the development of T-cell lymphomas with stereotypic karyotypes in which chromosome 15, containing the Myc oncogene, is gained with high prevalence. Expressing human MYC from chromosome 6 (MYCChr6) is sufficient to change the karyotype of these lymphomas to include universal chromosome 6 gains. Interestingly, while chromosome 15 is still gained in MYCChr6 tumors after genetic ablation of the endogenous Myc locus, this chromosome is not efficiently gained after deletion of one copy of Rad21, suggesting a synergistic effect of both MYC and RAD21 in driving chromosome 15 gains. Our results show that the initial detrimental effects of random missegregation are outbalanced by clonal selection, which is dictated by the chromosomal location and nature of certain genes and is sufficient to drive cancer with high prevalence.

Solving the chromosome puzzle of aneuploidy in cancer.

Chromosome instability (CIN) and aneuploidy are hallmarks of cancer cells, typically associated with aggressiveness and poor outcomes. Historically, the causative link between aneuploidy and cancer has been difficult to study due to its intrinsic complexity and the poor fitness of aneuploid cells. In this issue of Genes & Development, two companion papers (Trakala and colleagues [pp. 1079-1092] and Shoshani and colleagues [pp. 1093-1108]) exploited sophisticated mouse models to study the progression of aneuploidy from early phases to established tumors. Both groups observed that, while in the early nontumoral cells aneuploidy is characterized by random chromosomal gains, established tumors display a stereotypic karyotype with recurrent gains of only a few chromosomes. Thus, aneuploidy in tumors is not random but shows reproducible patterns of chromosomal changes induced by mechanisms that these two studies are beginning to unveil.

Poison Exon Splicing Regulates a Coordinated Network of SR Protein Expression during Differentiation and Tumorigenesis.

The RNA isoform repertoire is regulated by splicing factor (SF) expression, and alterations in SF levels are associated with disease. SFs contain ultraconserved poison exon (PE) sequences that exhibit greater identity across species than nearby coding exons, but their physiological role and molecular regulation is incompletely understood. We show that PEs in serine-arginine-rich (SR) proteins, a family of 14 essential SFs, are differentially spliced during induced pluripotent stem cell (iPSC) differentiation and in tumors versus normal tissues. We uncover an extensive cross-regulatory network of SR proteins controlling their expression via alternative splicing coupled to nonsense-mediated decay. We define sequences that regulate PE inclusion and protein expression of the oncogenic SF TRA2ß using an RNA-targeting CRISPR screen. We demonstrate location dependency of RS domain activity on regulation of TRA2ß-PE using CRISPR artificial SFs. Finally, we develop splice-switching antisense oligonucleotides to reverse the increased skipping of TRA2ß-PE detected in breast tumors, altering breast cancer cell viability, proliferation, and migration.

4D Genome Rewiring during Oncogene-Induced and Replicative Senescence.

To understand the role of the extensive senescence-associated 3D genome reorganization, we generated genome-wide chromatin interaction maps, epigenome, replication-timing, whole-genome bisulfite sequencing, and gene expression profiles from cells entering replicative senescence (RS) or upon oncogene-induced senescence (OIS). We identify senescence-associated heterochromatin domains (SAHDs). Differential intra- versus inter-SAHD interactions lead to the formation of senescence-associated heterochromatin foci (SAHFs) in OIS but not in RS. This OIS-specific configuration brings active genes located in genomic regions adjacent to SAHDs in close spatial proximity and favors their expression. We also identify DNMT1 as a factor that induces SAHFs by promoting HMGA2 expression. Upon DNMT1 depletion, OIS cells transition to a 3D genome conformation akin to that of cells in replicative senescence. These data show how multi-omics and imaging can identify critical features of RS and OIS and discover determinants of acute senescence and SAHF formation.

Integrated mutational landscape analysis of poorly differentiated high-grade neuroendocrine carcinoma of the uterine cervix.

High-grade neuroendocrine cervical cancers (NETc) are exceedingly rare, highly aggressive tumors. We analyzed 64 NETc tumor samples by whole-exome sequencing (WES). Human papillomavirus DNA was detected in 65.6% (42/64) of the tumors. Recurrent mutations were identified in PIK3CA, KMT2D/MLL2, K-RAS, ARID1A, NOTCH2, and RPL10. The top mutated genes included RB1, ARID1A, PTEN, KMT2D/MLL2, and WDFY3, a gene not yet implicated in NETc. Somatic CNV analysis identified two copy number gains (3q27.1 and 19q13.12) and five copy number losses (1p36.21/5q31.3/6p22.2/9q21.11/11p15.5). Also, gene fusions affecting the ACLY-CRHR1 and PVT1-MYC genes were identified in one of the eight samples subjected to RNA sequencing. To resolve evolutionary history, multiregion WES in NETc admixed with adenocarcinoma cells was performed (i.e., mixed-NETc). Phylogenetic analysis of mixed-NETc demonstrated that adenocarcinoma and neuroendocrine elements derive from a common precursor with mutations typical of adenocarcinomas. Over one-third (22/64) of NETc demonstrated a mutator phenotype of C > T at CpG consistent with deficiencies in MBD4, a member of the base excision repair (BER) pathway. Mutations in the PI3K/AMPK pathways were identified in 49/64 samples. We used two patient-derived-xenografts (PDX) (i.e., NET19 and NET21) to evaluate the activity of pan-HER (afatinib), PIK3CA (copanlisib), and ATR (elimusertib) inhibitors, alone and in combination. PDXs harboring alterations in the ERBB2/PI3K/AKT/mTOR/ATR pathway were sensitive to afatinib, copanlisib, and elimusertib (P < 0.001 vs. controls). However, combinations of copanlisib/afatinib and copanlisib/elimusertib were significantly more effective in controlling NETc tumor growth. These findings define the genetic landscape of NETc and suggest that a large subset of these highly lethal malignancies might benefit from existing targeted therapies.

PML isoforms: a molecular basis for PML pleiotropic functions.

PML is a stress-responsive protein that coordinates assembly of phase-separated nuclear aggregates, known as PML nuclear bodies (PML-NBs), where a large number of protein interactors and chromatin processes are finely regulated. Tampering with the PML gene produces a variety of phenotypic consequences that include promoting or interfering with tumor progression but the molecular underpinnings of PML pleiotropy are still elusive. In this review, we explore the contribution of PML splicing isoforms to PML-NB assorted activities. We describe recent literature indicating that distinct PML isoforms drive formation of specialized PML-NBs and perform unique functions and we suggest that future research efforts should delve into the contribution of isoform specificity to help elucidate the complex functionality of the PML gene.