Truncation- and motif-based pan-cancer analysis reveals tumor-suppressing kinases

HER2-driven breast cancer suppression by the JNK signaling pathway

The HER2⁺ subtype of human breast cancer is associated with the malignant transformation of luminal ductal cells of the mammary epithelium. The sequence analysis of tumor DNA identifies loss of function mutations and deletions of the MAP2K4 and MAP2K7 genes that encode direct activators of the JUN NH₂-terminal kinase (JNK). We report that in vitro studies of human mammary epithelial cells with CRISPR-induced mutations in the MAPK and MAP2K components of the JNK pathway caused no change in growth in 2D culture, but these mutations promoted epithelial cell proliferation in 3D culture. Analysis of gene expression signatures in 3D culture demonstrated similar changes caused by HER2 activation and JNK pathway loss. The mechanism of signal transduction cross-talk may be mediated, in part, by JNK-suppressed expression of integrin α6β4 that binds HER2 and amplifies HER2 signaling. These data suggest that HER2 activation and JNK pathway loss may synergize to promote breast cancer. To test this hypothesis, we performed in vivo studies using a mouse model of HER2⁺ breast cancer with Cre/loxP-mediated ablation of genes encoding JNK (Mapk8 and Mapk9) and the MAP2K (Map2k4 and Map2k7) that activate JNK in mammary epithelial cells. Kaplan-Meier analysis of tumor development demonstrated that JNK pathway deficiency promotes HER2⁺-driven breast cancer. Collectively, these data identify JNK pathway genes as potential suppressors for HER2⁺ breast cancer.

Characterization of genomic alterations and neoantigens and analysis of immune infiltration identified therapeutic and prognostic biomarkers in adenocarcinoma at the gastroesophageal junction

Objectives

Adenocarcinoma at the gastroesophageal junction (ACGEJ) refers to a malignant tumor that occurs at the esophagogastric junction. Despite some progress in targeted therapies for HER2, FGFR2, EGFR, MET, Claudin 18.2 and immune checkpoints in ACGEJ tumors, the 5-year survival rate of patients remains poor. Thus, it is urgent to explore genomic alterations and neoantigen characteristics of tumors and identify CD8+ T-cell infiltration-associated genes to find potential therapeutic targets and develop a risk model to predict ACGEJ patients’ overall survival (OS).

Methods

Whole-exome sequencing (WES) was performed on 55 paired samples from Chinese ACGEJ patients. Somatic mutations and copy number variations were detected by Strelka2 and FACETS, respectively. SigProfiler and SciClone were employed to decipher the mutation signature and clonal structure of each sample, respectively. Neoantigens were predicted using the MuPeXI pipeline. RNA sequencing (RNA-seq) data of ACGEJ samples from our previous studies and The Cancer Genome Atlas (TCGA) were used to identify genes significantly associated with CD8+ T-cell infiltration by weighted gene coexpression network analysis (WGCNA). To construct a risk model, we conducted LASSO and univariate and multivariate Cox regression analyses.

Results

Recurrent MAP2K7, RNF43 and RHOA mutations were found in ACGEJ tumors. The COSMIC signature SBS17 was associated with ACGEJ progression. CCNE1 and VEGFA were identified as putative CNV driver genes. PI3KCA and TP53 mutations conferred selective advantages to cancer cells. The Chinese ACGEJ patient neoantigen landscape was revealed for the first time, and 58 potential neoantigens common to TSNAdb and IEDB were identified. Compared with Siewert type II samples, Siewert type III samples had significant enrichment of the SBS17 signature, a lower TNFRSF14 copy number, a higher proportion of samples with complex clonal architecture and a higher neoantigen load. We identified 10 important CD8+ T-cell infiltration-related Hub genes (CCL5, CD2, CST7, GVINP1, GZMK, IL2RB, IKZF3, PLA2G2D, P2RY10 and ZAP70) as potential therapeutic targets from the RNA-seq data. Seven CD8+ T-cell infiltration-related genes (ADAM28, ASPH, CAMK2N1, F2R, STAP1, TP53INP2, ZC3H3) were selected to construct a prognostic model. Patients classified as high risk based on this model had significantly worse OS than low-risk patients, which was replicated in the TCGA-ACGEJ cohort.

Conclusions

This study provides new neoantigen-based immunotherapeutic targets for ACGEJ treatment and effective disease prognosis biomarkers.

Running title: EEF1AKMT3/MAP2K7/TP53 axis in gastric cancerThe EEF1AKMT3/MAP2K7/TP53 axis suppresses tumor invasiveness and metastasis in gastric cancer

The importance of methylation in the tumorigenic responses of nonhistone proteins, such as TP53, PTEN, RB1, AKT, and STAT3, has been emphasized in numerous studies. In parallel, the corresponding nonhistone protein methyltransferases have been acknowledged in the pathophysiology of cancer. Thus, this study aimed to explore the pathological role of a nonhistone methyltransferase in gastric cancer (GC), identify nonhistone substrate protein, and understand the underlying mechanism. Interestingly, among the 24 methyltransferases and methyltransferase family 16 (MTF16) proteins, EEF1AKMT3 (METTL21B) expression was prominently lower in GC tissues than in normal adjacent tissues and was associated with a worse prognosis. In addition, EEF1AKMT3-knockdown induced gastric tumor invasiveness and migration. Through gain and loss-of-function studies, mass spectrometry analysis, RNA-seq, and phospho-antibody array, we identified EEF1AKMT3 as a novel tumor-suppressive methyltransferase that catalyzes the monomethylation of MAP2K7 (MKK7) at K296, thereby decreasing the phosphorylation, ubiquitination, and degradation of TP53. Furthermore, EEF1AKMT3, p-MAP2K7, and TP53 protein levels were positively correlated in GC tissues. Collectively, our results delineate the tumor-suppressive function of the EEF1AKMT3/MAP2K7/TP53 signaling axis and suggest the dysregulation of the signaling axis as potential targeted therapy in GC.

A Signature of N6-methyladenosine Regulator-Related Genes Predicts Prognoses and Immune Responses for Head and Neck Squamous Cell Carcinoma

This study aimed to construct a signature of N⁶-methyladenosine (m6A) regulator-related genes that could be used for the prognosis of head and neck squamous cell carcinoma (HNSCC) and to clarify the molecular and immune characteristics and benefits of immune checkpoint inhibitor (ICI) therapy using the prognostic signature to define the subgroups of HNSCC. This study showed that eighteen m6A regulators were abnormally expressed in the Cancer Genome Atlas (TCGA) HNSCC tissues compared with those in normal tissues. We constructed a signature of 12 m6A regulator-related genes using the Cox risk model, combined with the least absolute shrinkage and selection operator (Lasso) variable screening algorithm. Based on the median of the signature risk score, the patients were divided into high- and low-risk groups. The Kaplan-Meier survival analyses showed that patients with high-risk scores demonstrated poorer overall survival (OS) than those with low-risk scores based on TCGA-HNSCC data (p <0.001). The OS of high-risk patients was significantly worse than that of low-risk patients in the GSE65858 (p <0.001) and International Cancer Genome Consortium (ICGC) oral cancer cohorts (p = 0.0089). Furthermore, immune infiltration analyses showed that 8 types of immune cell infiltration showed highly significant differences between the two risk groups (p <0.001). In the Imvigor210CoreBiologies dataset of patients who received ICIs, the objective response rate (ORR) of the low-risk group (32%) was significantly higher than that of the high-risk group (13%). Additionally, patients in the high-risk group presented with a more significant adverse OS than that of the low-risk group (p = 0.00032). GSE78220 also showed that the ORR of the low-risk group (64%) was higher than that of the high-risk group (43%) and the OS of low-risk patients was better than that of high-risk patients (p = 0.0064). The constructed prognostic signature, based on m6A regulator-related genes, could be used to effectively distinguish between prognoses for HNSCC patients. The prognostic signature was found to be related to the immune cell infiltration of HNSCC; it might help predict the responses and prognoses of ICIs during treatment.

DYRK3 contributes to differentiation and hypoxic control in neuroblastoma

Neuroblastoma (NB), a pediatric cancer of the peripheral sympathetic nervous system, represents the most frequent solid malignancy in infants. Treatment of high-risk patients is still challenging and, depending on the genetic make-up and involved risk factors, the 5-year survival rate can drop to only 30%. Here, we found that the expression of the Dual Specificity Tyrosine Phosphorylation Regulated Kinase 3 (DYRK3) is increased in NB and is associated with decreased survival in NB patients. We further identified DYRK3 as a cytoplasmic kinase in NB cells and found that its levels are increased by hypoxic conditions. Further mechanistic studies revealed that DYRK3 acts as a negative regulator of HIF-driven transcriptional responses, suggesting that it functions in a negative feedback loop controlling the hypoxic response. Moreover, DYRK3 negatively impacted on NB cell differentiation, proposing an oncogenic role of this kinase in the etiology of NB. In summary, we describe novel functions of the DYRK3 kinase in NB, which will help to further improve the understanding of this disease eventually leading to the design of improved therapeutic concepts.

The autoinhibited state of MKK4: Phosphorylation, putative dimerization and R134W mutant studied by molecular dynamics simulations

Protein kinases are crucial components of the cell-signalling machinery that orchestrate and convey messages to their downstream targets. Most often, kinases are activated upon a phosphorylation to their activation loop, which will shift the kinase into the active conformation. The Dual specificity mitogen-activated protein kinase kinase 4 (MKK4) exists in a unique conformation in its inactive unphosphorylated state, where its activation segment appears in a stable α-helical conformation. However, the precise role of this unique conformational state of MKK4 is unknown. Here, by all-atom molecular dynamics simulations (MD simulations), we show that this inactive state is unstable as monomer even when unphosphorylated and that the phosphorylation of the activation segment further destabilizes the autoinhibited α-helix. The specific phosphorylation pattern of the activation segment has also a unique influence on MKK4 dynamics. Furthermore, we observed that this specific inactive state is stable as a dimer, which becomes destabilized upon phosphorylation. Finally, we noticed that the most frequent MKK4 mutation observed in cancer, R134W, which role has not been disclosed to date, contributes to the dimer stability. Based on these data we postulate that MKK4 occurs as a dimer in its inactive autoinhibited state, providing an additional layer for its activity regulation.

A tumor-associated splice-isoform of MAP2K7 drives dedifferentiation in MBNL1-low cancers via JNK activation

Master splicing regulator MBNL1 shapes large transcriptomic changes that drive cellular differentiation during development. Here we demonstrate that MBNL1 is a suppressor of tumor dedifferentiation. We surveyed MBNL1 expression in matched tumor/normal pairs across The Cancer Genome Atlas and found that MBNL1 was down-regulated in several common cancers. Down-regulation of MBNL1 predicted poor overall survival in breast, lung, and stomach adenocarcinomas and increased relapse and distant metastasis in triple-negative breast cancer. Down-regulation of MBNL1 led to increased tumorigenic and stem/progenitor-like properties in vitro and in vivo. A discrete set of alternative splicing events (ASEs) are shared between MBNL1-low cancers and embryonic stem cells including a MAP2K7∆exon2 splice variant that leads to increased stem/progenitor-like properties via JNK activation. Accordingly, JNK inhibition is capable of reversing MAP2K7∆exon2-driven tumor dedifferentiation in MBNL1-low cancer cells. Our work elucidates an alternative-splicing mechanism that drives tumor dedifferentiation and identifies biomarkers that predict enhanced susceptibility to JNK inhibition.

Orange is the new black: Kinases are the new master regulators of tumor suppression

For many decades, kinases have predominantly been characterized as oncogenes and drivers of tumorigenesis, because activating mutations in kinases occur in cancer with high frequency. The oncogenic functions of kinases relate to their roles as growth factor receptors and as critical mediators of mitogen-activated pathways. Indeed, some of the most promising cancer therapeutic agents are kinase inhibitors. However, cancer genomics studies, especially screens that utilize high-throughput identification of loss-of-function somatic mutations, are beginning to shed light on a widespread role for kinases as tumor suppressors. The initial characterization of tumor-suppressing kinases- in particular members of the protein kinase C (PKC) family, MKK4 of the mitogen-activated protein kinase kinase family, and DAPK3 of the death-associated protein kinase family- laid the foundation for bioinformatic approaches that enable the identification of other tumor-suppressing kinases. In this review, we discuss the important role that kinases play as tumor suppressors, using several examples to illustrate the history of their discovery and highlight the modern approaches that presently aid in the identification of tumor-suppressing kinases. © 2018 IUBMB Life, 71(6):738-748, 2019.