Functional annotation of rare gene aberration drivers of pancreatic cancer

Simultaneous screening of overexpressed genes in breast cancer for oncogenic drivers and tumor dependencies

There are hundreds of genes typically overexpressed in breast cancer cells and it's often assumed that their overexpression contributes to cancer progression. However, the precise proportion of these overexpressed genes contributing to tumorigenicity remains unclear. To address this gap, we undertook a comprehensive screening of a diverse set of seventy-two genes overexpressed in breast cancer. This systematic screening evaluated their potential for inducing malignant transformation and, concurrently, assessed their impact on breast cancer cell proliferation and viability. Select genes including ALDH3B1, CEACAM5, IL8, PYGO2, and WWTR1, exhibited pronounced activity in promoting tumor formation and establishing gene dependencies critical for tumorigenicity. Subsequent investigations revealed that CEACAM5 overexpression triggered the activation of signaling pathways involving β-catenin, Cdk4, and mTOR. Additionally, it conferred a growth advantage independent of exogenous insulin in defined medium and facilitated spheroid expansion by inducing multiple layers of epithelial cells while preserving a hollow lumen. Furthermore, the silencing of CEACAM5 expression synergized with tamoxifen-induced growth inhibition in breast cancer cells. These findings underscore the potential of screening overexpressed genes for both oncogenic drivers and tumor dependencies to expand the repertoire of therapeutic targets for breast cancer treatment.

A de novo missense variant in EZH1 associated with developmental delay exhibits functional deficits in Drosophila melanogaster

EZH1, a polycomb repressive complex-2 component, is involved in a myriad of cellular processes. EZH1 represses transcription of downstream target genes through histone 3 lysine27 (H3K27) trimethylation (H3K27me3). Genetic variants in histone modifiers have been associated with developmental disorders, while EZH1 has not yet been linked to any human disease. However, the paralog EZH2 is associated with Weaver syndrome. Here we report a previously undiagnosed individual with a novel neurodevelopmental phenotype identified to have a de novo missense variant in EZH1 through exome sequencing. The individual presented in infancy with neurodevelopmental delay and hypotonia and was later noted to have proximal muscle weakness. The variant, p.A678G, is in the SET domain, known for its methyltransferase activity, and an analogous somatic or germline mutation in EZH2 has been reported in patients with B-cell lymphoma or Weaver syndrome, respectively. Human EZH1/2 are homologous to fly Enhancer of zeste (E(z)), an essential gene in Drosophila, and the affected residue (p.A678 in humans, p.A691 in flies) is conserved. To further study this variant, we obtained null alleles and generated transgenic flies expressing wildtype [E(z)WT] and the variant [E(z)A691G]. When expressed ubiquitously the variant rescues null-lethality similar to the wildtype. Overexpression of E(z)WT induces homeotic patterning defects but notably the E(z)A691G variant leads to dramatically stronger morphological phenotypes. We also note a dramatic loss of H3K27me2 and a corresponding increase in H3K27me3 in flies expressing E(z)A691G, suggesting this acts as a gain-of-function allele. In conclusion, here we present a novel EZH1 de novo variant associated with a neurodevelopmental disorder. Furthermore, we found that this variant has a functional impact in Drosophila.

Nicotinamide adenine dinucleotide kinase promotes lymph node metastasis of NSCLC via activating ID1 expression through BMP pathway

Metastasis is a significant cause of high mortality in lung cancer. Lymph node (LN) metastasis is the most common metastatic pathway in non-small cell lung cancer and the most crucial factor affecting the prognosis of NSCLC. Nevertheless, the molecular mechanism underlying metastasis is unknown. We demonstrated that higher NADK expression suggests worsened survival prognosis, and NADK expression positively correlates with the lymph node metastasis rate and TNM and AJCC stages in NSCLC patients. Moreover, patients with LN metastasis show higher NADK expression than those without LN metastasis. NADK can promote NSCLC progression by enhancing the migration, invasion, lymph node metastasis and growth of NSCLC cells. Mechanistically, NADK inhibits the ubiquitination and degradation of BMPR1A by interacting with Smurf1, further activating the BMPs signalling pathway and promoting ID1 transcription. In conclusion, NADK may be a potential diagnostic indicator and a novel therapeutic target for metastatic NSCLC.

The genomic landscape of familial glioma

Glioma is a rare brain tumor with a poor prognosis. Familial glioma is a subset of glioma with a strong genetic predisposition that accounts for approximately 5% of glioma cases. We performed whole-genome sequencing on an exploratory cohort of 203 individuals from 189 families with a history of familial glioma and an additional validation cohort of 122 individuals from 115 families. We found significant enrichment of rare deleterious variants of seven genes in both cohorts, and the most significantly enriched gene was HERC2 (P = 0.0006). Furthermore, we identified rare noncoding variants in both cohorts that were predicted to affect transcription factor binding sites or cause cryptic splicing. Last, we selected a subset of discovered genes for validation by CRISPR knockdown screening and found that DMBT1, HP1BP3, and ZCH7B3 have profound impacts on proliferation. This study performs comprehensive surveillance of the genomic landscape of familial glioma.

In vivo functional characterization of EGFR variants identifies novel drivers of glioblastoma

BACKGROUND

Glioblastoma is the most common and aggressive primary brain tumor. Large-scale sequencing initiatives have cataloged its mutational landscape in hopes of elucidating mechanisms driving this deadly disease. However, a major bottleneck in harnessing this data for new therapies is deciphering "driver" and "passenger" events amongst the vast volume of information.

METHODS

We utilized an autochthonous, in vivo screening approach to identify driver, EGFR variants. RNA-Seq identified unique molecular signatures of mouse gliomas across these variants, which only differ by a single amino acid change. In particular, we identified alterations to lipid metabolism, which we further validated through an unbiased lipidomics screen.

RESULTS

Our screen identified A289I as the most potent EGFR variant, which has previously not been characterized. One of the mechanisms through which A289I promotes gliomagenesis is to alter cellular triacylglycerides through MTTP. Knockout of Mttp in mouse gliomas, reduces gliomagenesis in multiple models.

CONCLUSIONS

EGFR variants that differ by a single amino acid residue differentially promote gliomagenesis. Among the identified mechanism that drives glioma growth include lipid metabolism through MTTP. Understanding triacylglyceride accumulation may present a prospective therapeutic pathway for this deadly disease.

Glioma epileptiform activity and progression are driven by IGSF3-mediated potassium dysregulation

Seizures are a frequent pathophysiological feature of malignant glioma. Recent studies implicate peritumoral synaptic dysregulation as a driver of brain hyperactivity and tumor progression; however, the molecular mechanisms that govern these phenomena remain elusive. Using scRNA-seq and intraoperative patient ECoG recordings, we show that tumors from seizure patients are enriched for gene signatures regulating synapse formation. Employing a human-to-mouse in vivo functionalization pipeline to screen these genes, we identify IGSF3 as a mediator of glioma progression and dysregulated neural circuitry that manifests as spreading depolarization (SD). Mechanistically, we discover that IGSF3 interacts with Kir4.1 to suppress potassium buffering and found that seizure patients exhibit reduced expression of potassium handlers in proliferating tumor cells. In vivo imaging reveals that dysregulated synaptic activity emanates from the tumor-neuron interface, which we confirm in patients. Our studies reveal that tumor progression and seizures are enabled by ion dyshomeostasis and identify SD as a driver of disease.

A de novo missense variant in EZH1 associated with developmental delay exhibits functional deficits in Drosophila melanogaster

EZH1 ( Enhancer of Zeste, homolog 1) , a Polycomb Repressive Complex-2 (PRC2) component, is involved in a myriad of cellular processes through modifying histone 3 lysine27 (H3K27) residues. EZH1 represses transcription of downstream target genes through H3K27 trimethylation (H3K27me3). Genetic mutations in histone modifiers have been associated with developmental disorders, while EZH1 has not yet been linked to any human disease. However, the paralog EZH2 is associated with Weaver syndrome. Here we report a previously undiagnosed individual with a novel neurodevelopmental phenotype identified to have a de novo variant in EZH1 , p.Ala678Gly, through exome sequencing. The individual presented in infancy with neurodevelopmental delay and hypotonia and was later noted to have proximal muscle weakness. The variant, p.A678G, is in the SET domain, known for its methyltransferase activity, and was the best candidate variant found in the exome. Human EZH1 / 2 are homologous to fly Enhancer of zeste E(z) , an essential gene in flies, and the residue (A678 in humans, A691 in Drosophila ) is conserved. To further study this variant, we obtained Drosophila null alleles and generated transgenic flies expressing wild-type (E(z) WT ) and the variant (E(z) A691G ) . The E(z) A691G variant led to hyper H3K27me3 while the E(z) WT did not, suggesting this is as a gain-of-function allele. When expressed under the tubulin promotor in vivo the variant rescued null-lethality similar to wild-type but the E(z) A691G flies exhibit bang sensitivity and shortened lifespan. In conclusion, here we present a novel EZH1 de novo variant associated with a neurodevelopmental disorder. Furthermore, we found that this variant has a functional impact in Drosophila . Biochemically this allele leads to increased H3K27me3 suggesting gain-of-function, but when expressed in adult flies the E(z) A691G has some characteristics of partial loss-of-function which may suggest it is a more complex allele in vivo .

Molecular properties and regulation of NAD+ kinase (NADK)

Nicotinamide adenine dinucleotide (NAD⁺) kinase (NADK) phosphorylates NAD⁺, thereby producing nicotinamide adenine dinucleotide phosphate (NADP). Both NADK genes and the NADP(H)-producing mechanism are evolutionarily conserved among archaea, bacteria, plants and mammals. In mammals, NADK is activated by phosphorylation and protein-protein interaction. Recent studies conducted using genetically altered models validate the essential role of NADK in cellular redox homeostasis and metabolism in multicellular organisms. Here, we describe the evolutionary conservation, molecular properties, and signaling mechanisms and discuss the pathophysiological significance of NADK.

(-)-Epigallocatechin Gallate is a Noncompetitive Inhibitor of NAD Kinase

Nicotinamide adenine dinucleotide kinase (NADK) controls the intracellular NADPH content and provides reducing power for the synthesis of macromolecules and anti-ROS. Moreover, NADK is considered to be a synthetic lethal gene for KRAS mutations. To discover NADK-targeted probes, a high-throughput screening assay was established and optimized with a Z factor of 0.71. The natural product (-)-epigallocatechin gallate (EGCG) was found to be a noncompetitive inhibitor of NADK with K i = 3.28 ± 0.32 μΜ. The direct binding of EGCG to NADK was determined by several biophysical methods, including NMR spectroscopy, surface plasmon resonance (SPR) assay, and hydrogen-deuterium exchange mass spectrometry (HDX-MS). The SPR assay showed a K d of 1.78 ± 1.15 μΜ. The HDX-MS experiment showed that EGCG was bound at the non-substrate-binding sites of NADK. Besides, binding mode prediction and derivative activity analysis revealed a potential structure-activity relationship between EGCG and NADK. Furthermore, EGCG can specifically inhibit the proliferation of KRAS-mutated lung cancer cell lines without affecting KRAS wild-type lung cancer cell lines.

G6PD-mediated increase in de novo NADP+ biosynthesis promotes antioxidant defense and tumor metastasis

Metastasizing cancer cells are able to withstand high levels of oxidative stress through mechanisms that are poorly understood. Here, we show that under various oxidative stress conditions, pancreatic cancer cells markedly expand NADPH and NADP+ pools. This expansion is due to up-regulation of glucose-6-phosphate dehydrogenase (G6PD), which stimulates the cytoplasmic nicotinamide adenine dinucleotide kinase (NADK1) to produce NADP+ while converting NADP+ to NADPH. G6PD is activated by the transcription factor TAp73, which is, in turn, regulated by two pathways. Nuclear factor-erythroid 2 p45-related factor-2 suppresses expression of the ubiquitin ligase PIRH2, stabilizing the TAp73 protein. Checkpoint kinases 1/2 and E2F1 induce expression of the TAp73 gene. Levels of G6PD and its upstream activators are elevated in metastatic pancreatic cancer. Knocking down G6PD impedes pancreatic cancer metastasis, whereas forced G6PD expression promotes it. These findings reveal an intracellular network that maintains redox homeostasis through G6PD-mediated increase in de novo NADP+ biosynthesis, which may be co-opted by tumor cells to enable metastasis.

The role of ROS in tumour development and progression

Eukaryotic cells have developed complex systems to regulate the production and response to reactive oxygen species (ROS). Different ROS control diverse aspects of cell behaviour from signalling to death, and deregulation of ROS production and ROS limitation pathways are common features of cancer cells. ROS also function to modulate the tumour environment, affecting the various stromal cells that provide metabolic support, a blood supply and immune responses to the tumour. Although it is clear that ROS play important roles during tumorigenesis, it has been difficult to reliably predict the effect of ROS modulating therapies. We now understand that the responses to ROS are highly complex and dependent on multiple factors, including the types, levels, localization and persistence of ROS, as well as the origin, environment and stage of the tumours themselves. This increasing understanding of the complexity of ROS in malignancies will be key to unlocking the potential of ROS-targeting therapies for cancer treatment.

Drosophila functional screening of de novo variants in autism uncovers damaging variants and facilitates discovery of rare neurodevelopmental diseases

Individuals with autism spectrum disorder (ASD) exhibit an increased burden of de novo mutations (DNMs) in a broadening range of genes. While these studies have implicated hundreds of genes in ASD pathogenesis, which DNMs cause functional consequences in vivo remains unclear. We functionally test the effects of ASD missense DNMs using Drosophila through "humanization" rescue and overexpression-based strategies. We examine 79 ASD variants in 74 genes identified in the Simons Simplex Collection and find 38% of them to cause functional alterations. Moreover, we identify GLRA2 as the cause of a spectrum of neurodevelopmental phenotypes beyond ASD in 13 previously undiagnosed subjects. Functional characterization of variants in ASD candidate genes points to conserved neurobiological mechanisms and facilitates gene discovery for rare neurodevelopmental diseases.

Restriction of extracellular lipids renders pancreatic cancer dependent on autophagy

BACKGROUND

KRAS is the predominant oncogene mutated in pancreatic ductal adenocarcinoma (PDAC), the fourth cause of cancer-related deaths worldwide. Mutant KRAS-driven tumors are metabolically programmed to support their growth and survival, which can be used to identify metabolic vulnerabilities. In the present study, we aimed to understand the role of extracellularly derived fatty acids in KRAS-driven pancreatic cancer.

METHODS

To assess the dependence of PDAC cells on extracellular fatty acids we employed delipidated serum or RNAi-mediated suppression of ACSL3 (to inhibit the activation and cellular retention of extracellular fatty acids) followed by cell proliferation assays, qPCR, apoptosis assays, immunoblots and fluorescence microscopy experiments. To assess autophagy in vivo, we employed the KrasG12D/+;p53flox/flox;Pdx1-CreERT2 (KPC) mice crossed with Acsl3 knockout mice, and to assess the efficacy of the combination therapy of ACSL3 and autophagy inhibition we used xenografted human cancer cell-derived tumors in immunocompromised mice.

RESULTS

Here we show that depletion of extracellularly derived lipids either by serum lipid restriction or suppression of ACSL3, triggers autophagy, a process that protects PDAC cells from the reduction of bioenergetic intermediates. Combined extracellular lipid deprivation and autophagy inhibition exhibits anti-proliferative and pro-apoptotic effects against PDAC cell lines in vitro and promotes suppression of xenografted human pancreatic cancer cell-derived tumors in mice. Therefore, we propose lipid deprivation and autophagy blockade as a potential co-targeting strategy for PDAC treatment.

CONCLUSIONS

Our work unravels a central role of extracellular lipid supply in ensuring fatty acid provision in cancer cells, unmasking a previously unappreciated metabolic vulnerability of PDAC cells.

Connections between metabolism and epigenetic modifications in cancer

How cells sense and respond to environmental changes is still a key question. It has been identified that cellular metabolism is an important modifier of various epigenetic modifications, such as DNA methylation, histone methylation and acetylation and RNA N6-methyladenosine (m6A) methylation. This closely links the environmental nutrient availability to the maintenance of chromatin structure and gene expression, and is crucial to regulate cellular homeostasis, cell growth and differentiation. Cancer metabolic reprogramming and epigenetic alterations are widely observed, and facilitate cancer development and progression. In cancer cells, oncogenic signaling-driven metabolic reprogramming modifies the epigenetic landscape via changes in the key metabolite levels. In this review, we briefly summarized the current evidence that the abundance of key metabolites, such as S-adenosyl methionine (SAM), acetyl-CoA, α-ketoglutarate (α-KG), 2-hydroxyglutarate (2-HG), uridine diphospho-N-acetylglucosamine (UDP-GlcNAc) and lactate, affected by metabolic reprogramming plays an important role in dynamically regulating epigenetic modifications in cancer. An improved understanding of the roles of metabolic reprogramming in epigenetic regulation can contribute to uncover the underlying mechanisms of metabolic reprogramming in cancer development and identify the potential targets for cancer therapies.

NADK is activated by oncogenic signaling to sustain pancreatic ductal adenocarcinoma

Metabolic adaptations and the signaling events that control them promote the survival of pancreatic ductal adenocarcinoma (PDAC) at the fibrotic tumor site, overcoming stresses associated with nutrient and oxygen deprivation. Recently, rewiring of NADPH production has been shown to play a key role in this process. NADPH is recycled through reduction of NADP+ by several enzymatic systems in cells. However, de novo NADP+ is synthesized only through one known enzymatic reaction, catalyzed by NAD+ kinase (NADK). In this study, we show that oncogenic KRAS promotes protein kinase C (PKC)-mediated NADK phosphorylation, leading to its hyperactivation, thus sustaining both NADP+ and NADPH levels in PDAC cells. Together, our data show that increased NADK activity is an important adaptation driven by oncogenic signaling. Our findings indicate that NADK could serve as a much-needed therapeutic target for PDAC.

Next-generation sequencing-based identification of EGFR and NOTCH2 complementary mutations in non-small cell lung cancer

Although targeted therapy has emerged as an effective treatment strategy for non-small cell lung cancer (NSCLC), some patients cannot benefit from such therapy due to the limited number of therapeutic targets. The present study aimed to identify mutated genes associated with clinicopathological characteristics and prognosis and to screen for mutations that are not concurrent with applicable drug target sites in patients with NSCLC. Tumor tissue and blood samples were obtained from 97 patients with NSCLC. A lung cancer-specific panel of 55 genes was established and analyzed using next-generation sequencing (NGS). The results obtained from the clinical cohort were compared with the NSCLC dataset from The Cancer Genome Atlas (TCGA). Subsequently, 25 driver genes were identified by taking the intersection of the 55 lung-cancer-specific genes with three databases, namely, the Catalog of Somatic Mutations in Cancer database, the Network of Cancer Genes database and Vogelstein's list. Functional annotation and protein-protein interaction analysis were conducted on these 25 driver genes. The χ² test and logistic regression were used to evaluate the association between mutations in the 25 driver genes and the clinicopathological characteristics of 97 patients, and phosphatase and tensin homolog (PTEN) and kirsten rat sarcoma viral oncogene homolog (KRAS) were associated with stage at diagnosis and sex, respectively, while epidermal growth factor receptor (EGFR) was associated with sex, stage at diagnosis, metastasis, CEA and CYFRA21-1. Moreover, the association between the 25 driver gene mutations and overall survival were examined using Cox regression analysis. Age and Notch homolog 2 (NOTCH2) mutations were independent prognostic factors in TCGA dataset. The correlations between statistically significant mutations in EGFR, KRAS, PTEN and NOTCH2 were further examined, both in the clinical data and TCGA dataset. There was a negative correlation between EGFR and NOTCH2 mutations (correlation coefficient, −0.078; P=0.027). Thus, the present study highlights the importance of NOTCH2 mutations and might provide novel therapeutic options for patients with NSCLC who do not harbor EGFR mutations.

Precision oncology for breast cancer through clinical trials

Precision treatment for breast cancers has made several notable advances in recent decades, but challenges of tumor heterogeneity, drug resistance, and aggressive recurrence and metastases remain. To meet and overcome these challenges, we must refine our understanding of breast subtypes and treatment biomarkers according to the knowledge afforded across the spectrum of 'omics assays. A critical aspect of harnessing this knowledge into actionable biomarkers for treatment decision relies on our ability to integrate knowledge across data types and leverage our insight in evidence-based clinical trials. We review recent advances in cutting-edge clinical trials for precision treatment of breast cancer, including chemotherapies, targeted therapies, immunotherapies, and combination therapies. We comment on promising future areas of development for this exciting point in precision breast cancer research.

In cancer, all roads lead to NADPH

Cancer cells require increased levels of NADPH for increased nucleotide synthesis and for protection from ROS. Recent studies show that increased NADPH is generated in several ways. Activated AKT phosphorylates NAD kinase (NADK), increasing its activity. NADP formed, is rapidly converted to NADPH by glucose 6-phosphate dehydrogenase and malic enzymes, overexpressed in tumor cells with mutant p53. Calmodulin, overexpressed in some cancers, also increases NADK activity. Also, in IDH1/2 mutant cancer, NADPH serves as the cofactor to generate D-2 hydroxyglutarate, an oncometabolite. The requirement of cancer cells for elevated levels of NADPH provides an opportunity to target its synthesis for cancer treatment.

Biological significance of KRAS mutant allele expression in ovarian endometriosis

KRAS is the most frequently mutated in ovarian endometriosis. However, it is unclear whether the KRAS mutant allele's mRNA is expressed and plays a biological role in ovarian endometriosis. Here, we performed mutation-specific RNA in situ hybridization to evaluate mutant allele expression of KRAS p.G12V, the most frequently detected mutation in ovarian endometriosis in our previous study, in formalin-fixed paraffin-embedded tissue (FFPE) samples of ovarian endometriosis, cancer cell lines, and ovarian cancers. First, we verified that mutant or wild-type allele of KRAS were expressed in all 5 cancer cell lines and 9 ovarian cancer cases corresponding to the mutation status. Next, we applied this assay to 26 ovarian endometriosis cases, and observed mutant allele expression of KRAS p.G12V in 10 cases. Mutant or wild-type allele of KRAS were expressed in line with mutation status in 12 available endometriosis cases for which KRAS gene sequence was determined. Comparison of clinical features between ovarian endometriosis with KRAS p.G12V mutant allele expression and with KRAS wild-type showed that KRAS p.G12V mutant allele expression was significantly associated with inflammation in ovarian endometriosis. Finally, we assessed the spatial distribution of KRAS mutant allele expression in 5 endometriosis cases by performing multiregional sampling. Intratumor heterogeneity of KRAS mutant allele expression was observed in two endometriosis cases, whereas the spatial distribution of KRAS p.G12V mutation signals were diffuse and homogenous in ovarian cancer. In conclusion, evaluation of oncogene mutant expression will be useful for clarifying the biological significance of oncogene mutations in benign tumors.

Protein kinase C activates NAD kinase in human neutrophils

NAD kinase (NADK) is required for the de novo synthesis of NADP⁺ from NAD⁺. In neutrophils, NADK plays an essential role by providing sufficient levels of NADPH to support a robust oxidative burst. Activation of NADPH oxidase-2 (NOX-2) in neutrophils by stimulators of protein kinase C (PKC), such as phorbol myristate acetate (PMA), results in the rapid generation of superoxide at the expense of oxidation of NADPH to NADP⁺. In this study, we measured the levels of pyridine nucleotides following the addition of PMA to neutrophils. PMA elicited a rapid increase in NADP⁺ in neutrophils, which was not due to oxidation of NADPH, the levels of which also rose. This was mirrored by a rapid reduction in NAD⁺ levels, suggesting that NADK had been activated. PMA-induced depletion of NAD⁺ in neutrophils was blocked by PKC inhibitors, but was not dependent on NOX-2, as it was not blocked by the NOX inhibitor, diphenyleneiodonium. PMA also increased NADK activity in neutrophil lysates as well as NADK phosphorylation, as revealed by a monoclonal antibody selective for phospho-NADK. Human recombinant NADK was phosphorylated by PKCδ, resulting in increased immunoreactivity, but unchanged enzyme activity, suggesting that PKC-induced phosphorylation alone is insufficient to increase NADK activity in neutrophils. This leads us to speculate that phosphorylation of NADK promotes the dissociation of an inhibitory molecule from a complex, thereby increasing enzyme activity. Activation of NADK by PKC in phagocytic cells could be critical for the rapid provision of sufficient levels of superoxide for host defence against invading microorganisms.

NADPH homeostasis in cancer: functions, mechanisms and therapeutic implications

Nicotinamide adenine dinucleotide phosphate (NADPH) is an essential electron donor in all organisms, and provides the reducing power for anabolic reactions and redox balance. NADPH homeostasis is regulated by varied signaling pathways and several metabolic enzymes that undergo adaptive alteration in cancer cells. The metabolic reprogramming of NADPH renders cancer cells both highly dependent on this metabolic network for antioxidant capacity and more susceptible to oxidative stress. Modulating the unique NADPH homeostasis of cancer cells might be an effective strategy to eliminate these cells. In this review, we summarize the current existing literatures on NADPH homeostasis, including its biological functions, regulatory mechanisms and the corresponding therapeutic interventions in human cancers, providing insights into therapeutic implications of targeting NADPH metabolism and the associated mechanism for cancer therapy.

In vivo functional screening for systems-level integrative cancer genomics

With the genetic portraits of all major human malignancies now available, we next face the challenge of characterizing the function of mutated genes, their downstream targets, interactions and molecular networks. Moreover, poorly understood at the functional level are also non-mutated but dysregulated genomes, epigenomes or transcriptomes. Breakthroughs in manipulative mouse genetics offer new opportunities to probe the interplay of molecules, cells and systemic signals underlying disease pathogenesis in higher organisms. Herein, we review functional screening strategies in mice using genetic perturbation and chemical mutagenesis. We outline the spectrum of genetic tools that exist, such as transposons, CRISPR and RNAi and describe discoveries emerging from their use. Genome-wide or targeted screens are being used to uncover genomic and regulatory landscapes in oncogenesis, metastasis or drug resistance. Versatile screening systems support experimentation in diverse genetic and spatio-temporal settings to integrate molecular, cellular or environmental context-dependencies. We also review the combination of in vivo screening and barcoding strategies to study genetic interactions and quantitative cancer dynamics during tumour evolution. These scalable functional genomics approaches are transforming our ability to interrogate complex biological systems.

ACSL3-PAI-1 signaling axis mediates tumor-stroma cross-talk promoting pancreatic cancer progression

Pancreatic ductal adenocarcinoma (PDAC) is characterized by marked fibrosis and low immunogenicity, features that are linked to treatment resistance and poor clinical outcomes. Therefore, understanding how PDAC regulates the desmoplastic and immune stromal components is of great clinical importance. We found that acyl-CoA synthetase long-chain 3 (ACSL3) is up-regulated in PDAC and correlates with increased fibrosis. Our in vivo results show that Acsl3 knockout hinders PDAC progression, markedly reduces tumor fibrosis and tumor-infiltrating immunosuppressive cells, and increases cytotoxic T cell infiltration. This effect is, at least in part, due to decreased plasminogen activator inhibitor-1 (PAI-1) secretion from tumor cells. Accordingly, PAI-1 expression in PDAC positively correlates with markers of fibrosis and immunosuppression and predicts poor patient survival. We found that PAI-1 pharmacological inhibition strongly enhances chemo- and immunotherapeutic response against PDAC, increasing survival of mice. Thus, our results unveil ACSL3-PAI-1 signaling as a requirement for PDAC progression with druggable attributes.

Bioinformatics Data Mining Repurposes the JAK2 (Janus Kinase 2) Inhibitor Fedratinib for Treating Pancreatic Ductal Adenocarcinoma by Reversing the KRAS (Kirsten Rat Sarcoma 2 Viral Oncogene Homolog)-Driven Gene Signature

Pancreatic ductal adenocarcinoma (PDAC) is still one of the most aggressive and lethal cancer types due to the late diagnosis, high metastatic potential, and drug resistance. The development of novel therapeutic strategies is urgently needed. KRAS (Kirsten rat sarcoma 2 viral oncogene homolog) is the major driver mutation gene for PDAC tumorigenesis. In this study, we mined cancer genomics data and identified a common KRAS-driven gene signature in PDAC, which is related to cell–cell and cell–extracellular matrix (ECM) interactions. Higher expression of this gene signature was associated with poorer overall survival of PDAC patients. Connectivity Map (CMap) analysis and drug sensitivity profiling predicted that a clinically approved JAK2 (Janus kinase 2)-selective inhibitor, fedratinib (also known as TG-101348), could reverse the KRAS-driven gene signature and exhibit KRAS-dependent anticancer activity in PDAC cells. As an approved treatment for myelofibrosis, the pharmacological and toxicological profiles of fedratinib have been well characterized. It may be repurposed for treating KRAS-driven PDAC in the future.

Beyond the Genomic Mutation: Rethinking the Molecular Biomarkers of K-RAS Dependency in Pancreatic Cancers

Oncogenic v-Ki-ras2 Kirsten rat sarcoma viral oncogene homolog (K-RAS) plays a key role in the development and maintenance of pancreatic ductal adenocarcinoma (PDAC). The targeting of K-RAS would be beneficial to treat tumors whose growth depends on active K-RAS. The analysis of K-RAS genomic mutations is a clinical routine; however, an emerging question is whether the mutational status is able to identify tumors effectively dependent on K-RAS for tailoring targeted therapies. With the emergence of novel K-RAS inhibitors in clinical settings, this question is relevant. Several studies support the notion that the K-RAS mutation is not a sufficient biomarker deciphering the effective dependency of the tumor. Transcriptomic and metabolomic profiles of tumors, while revealing K-RAS signaling complexity and K-RAS-driven molecular pathways crucial for PDAC growth, are opening the opportunity to specifically identify K-RAS-dependent- or K-RAS-independent tumor subtypes by using novel molecular biomarkers. This would help tumor selection aimed at tailoring therapies against K-RAS. In this review, we will present studies about how the K-RAS mutation can also be interpreted in a state of K-RAS dependency, for which it is possible to identify specific K-RAS-driven molecular biomarkers in certain PDAC subtypes, beyond the genomic K-RAS mutational status.

Molecular Landscape of BRAF-Mutant NSCLC Reveals an Association Between Clonality and Driver Mutations and Identifies Targetable Non-V600 Driver Mutations

INTRODUCTION

Approximately 4% of NSCLC harbor BRAF mutations, and approximately 50% of these are non-V600 mutations. Treatment of tumors harboring non-V600 mutations is challenging because of functional heterogeneity and lack of knowledge regarding their clinical significance and response to targeted agents.

METHODS

We conducted an integrative analysis of BRAF non-V600 mutations using genomic profiles of BRAF-mutant NSCLC from the Guardant360 database. BRAF mutations were categorized by clonality and class (1 and 2: RAS-independent; 3: RAS-dependent). Cell viability assays were performed in Ba/F3 models. Drug screens were performed in NSCLC cell lines.

RESULTS

A total of 305 unique BRAF mutations were identified. Missense mutations were most common (276, 90%), and 45% were variants of unknown significance. F468S and N581Y were identified as novel activating mutations. Class 1 to 3 mutations had higher clonality than mutations of unknown class (p < 0.01). Three patients were treated with MEK with or without BRAF inhibitors. Patients harboring G469V and D594G mutations did not respond, whereas a patient with the L597R mutation had a durable response. Trametinib with or without dabrafenib, LXH254, and lifirafenib had more potent inhibition of BRAF non-V600-mutant NSCLC cell lines than other MEK, BRAF, and ERK inhibitors, comparable with the inhibition of BRAF V600E cell line.

CONCLUSIONS

In BRAF-mutant NSCLC, clonality is higher in known functional mutations and may allow identification of variants of unknown significance that are more likely to be oncogenic drivers. Our data indicate that certain non-V600 mutations are responsive to MEK and BRAF inhibitors. This integration of genomic profiling and drug sensitivity may guide the treatment for BRAF-mutant NSCLC.

Spatially resolved analyses link genomic and immune diversity and reveal unfavorable neutrophil activation in melanoma

Complex tumor microenvironmental (TME) features influence the outcome of cancer immunotherapy (IO). Here we perform immunogenomic analyses on 67 intratumor sub-regions of a PD-1 inhibitor-resistant melanoma tumor and 2 additional metastases arising over 8 years, to characterize TME interactions. We identify spatially distinct evolution of copy number alterations influencing local immune composition. Sub-regions with chromosome 7 gain display a relative lack of leukocyte infiltrate but evidence of neutrophil activation, recapitulated in The Cancer Genome Atlas (TCGA) samples, and associated with lack of response to IO across three clinical cohorts. Whether neutrophil activation represents cause or consequence of local tumor necrosis requires further study. Analyses of T-cell clonotypes reveal the presence of recurrent priming events manifesting in a dominant T-cell clonotype over many years. Our findings highlight the links between marked levels of genomic and immune heterogeneity within the physical space of a tumor, with implications for biomarker evaluation and immunotherapy response.

Differential expression of MAGEA6 toggles autophagy to promote pancreatic cancer progression

The melanoma-associated antigen family A (MAGEA) antigens are expressed in a wide variety of malignant tumors but not in adult somatic cells, rendering them attractive targets for cancer immunotherapy. Here we show that a number of cancer-associated MAGEA mutants that undergo proteasome-dependent degradation in vitro could negatively impact their utility as immunotherapeutic targets. Importantly, in pancreatic ductal adenocarcinoma cell models, MAGEA6 suppresses macroautophagy (autophagy). The inhibition of autophagy is released upon MAGEA6 degradation, which can be induced by nutrient deficiency or by acquisition of cancer-associated mutations. Using xenograft mouse models, we demonstrated that inhibition of autophagy is critical for tumor initiation whereas reinstitution of autophagy as a consequence of MAGEA6 degradation contributes to tumor progression. These findings could inform cancer immunotherapeutic strategies for targeting MAGEA antigens and provide mechanistic insight into the divergent roles of MAGEA6 during pancreatic cancer initiation and progression.

PBMC gene expression profiles of female Bangladeshi adults chronically exposed to arsenic-contaminated drinking water

Arsenic, a class I human carcinogen, is ubiquitously found throughout the environment and around the globe, posing a great public health concern. Notably, Bangladesh and regions of West Bengal have been found to have high levels (0.5-4600 μg/L) of arsenic drinking water contamination, and approximately 50 million of the world's 200 million people chronically exposed to arsenic in Bangladesh alone. This study was carried out to examine genome-wide gene expression changes in individuals chronically exposed to arsenic-contaminated drinking water. Our study population includes twenty-nine Bangladeshi female participants with urinary arsenic levels ranging from 22.32 to 1828.12 μg/g creatinine. RNA extracted from peripheral blood mononuclear cells (PBMCs) were evaluated using RNA-Sequencing analysis. Our results indicate that a total of 1,054 genes were significantly associated with increasing urinary arsenic levels (FDR p < 0.05), which include 418 down-regulated and 636 up-regulated genes. Further Ingenuity Pathway Analysis revealed potential target genes (DAPK1, EGR2, APP), microRNAs (miR-155, -338, -210) and pathways (NOTCH signaling pathway) related to arsenic carcinogenesis. The selection of female-only participants provides a homogenous study population since arsenic has significant sex dependent effects, and the wide exposure range provides new insight for key gene expression changes that correlate with increasing urinary arsenic levels.

Comprehensive assessment of computational algorithms in predicting cancer driver mutations

BACKGROUND

The initiation and subsequent evolution of cancer are largely driven by a relatively small number of somatic mutations with critical functional impacts, so-called driver mutations. Identifying driver mutations in a patient's tumor cells is a central task in the era of precision cancer medicine. Over the decade, many computational algorithms have been developed to predict the effects of missense single-nucleotide variants, and they are frequently employed to prioritize mutation candidates. These algorithms employ diverse molecular features to build predictive models, and while some algorithms are cancer-specific, others are not. However, the relative performance of these algorithms has not been rigorously assessed.

RESULTS

We construct five complementary benchmark datasets: mutation clustering patterns in the protein 3D structures, literature annotation based on OncoKB, TP53 mutations based on their effects on target-gene transactivation, effects of cancer mutations on tumor formation in xenograft experiments, and functional annotation based on in vitro cell viability assays we developed including a new dataset of ~ 200 mutations. We evaluate the performance of 33 algorithms and found that CHASM, CTAT-cancer, DEOGEN2, and PrimateAI show consistently better performance than the other algorithms. Moreover, cancer-specific algorithms show much better performance than those designed for a general purpose.

CONCLUSIONS

Our study is a comprehensive assessment of the performance of different algorithms in predicting cancer driver mutations and provides deep insights into the best practice of computationally prioritizing cancer mutation candidates for end-users and for the future development of new algorithms.

PIK3CA variants selectively initiate brain hyperactivity during gliomagenesis

Glioblastoma is a universally lethal form of brain cancer that exhibits an array of pathophysiological phenotypes, many of which are mediated by interactions with the neuronal microenvironment1,2. Recent studies have shown that increases in neuronal activity have an important role in the proliferation and progression of glioblastoma3,4. Whether there is reciprocal crosstalk between glioblastoma and neurons remains poorly defined, as the mechanisms that underlie how these tumours remodel the neuronal milieu towards increased activity are unknown. Here, using a native mouse model of glioblastoma, we develop a high-throughput in vivo screening platform and discover several driver variants of PIK3CA. We show that tumours driven by these variants have divergent molecular properties that manifest in selective initiation of brain hyperexcitability and remodelling of the synaptic constituency. Furthermore, secreted members of the glypican (GPC) family are selectively expressed in these tumours, and GPC3 drives gliomagenesis and hyperexcitability. Together, our studies illustrate the importance of functionally interrogating diverse tumour phenotypes driven by individual, yet related, variants and reveal how glioblastoma alters the neuronal microenvironment.

The PI3K-AKT network at the interface of oncogenic signalling and cancer metabolism

The altered metabolic programme of cancer cells facilitates their cell-autonomous proliferation and survival. In normal cells, signal transduction pathways control core cellular functions, including metabolism, to couple the signals from exogenous growth factors, cytokines or hormones to adaptive changes in cell physiology. The ubiquitous, growth factor-regulated phosphoinositide 3-kinase (PI3K)-AKT signalling network has diverse downstream effects on cellular metabolism, through either direct regulation of nutrient transporters and metabolic enzymes or the control of transcription factors that regulate the expression of key components of metabolic pathways. Aberrant activation of this signalling network is one of the most frequent events in human cancer and serves to disconnect the control of cell growth, survival and metabolism from exogenous growth stimuli. Here we discuss our current understanding of the molecular events controlling cellular metabolism downstream of PI3K and AKT and of how these events couple two major hallmarks of cancer: growth factor independence through oncogenic signalling and metabolic reprogramming to support cell survival and proliferation.

Computer-aided drug repurposing for cancer therapy: Approaches and opportunities to challenge anticancer targets

Despite huge efforts made in academic and pharmaceutical worldwide research, current anticancer therapies achieve effective treatment in a limited number of neoplasia cases only. Oncology terms such as big killers - to identify tumours with yet a high mortality rate - or undruggable cancer targets, and chemoresistance, represent the current therapeutic debacle of cancer treatments. In addition, metastases, tumour microenvironments, tumour heterogeneity, metabolic adaptations, and immunotherapy resistance are essential features controlling tumour response to therapies, but still, lack effective therapeutics or modulators. In this scenario, where the pharmaceutical productivity and drug efficacy in oncology seem to have reached a plateau, the so-called drug repurposing - i.e. the use of old drugs, already in clinical use, for a different therapeutic indication - is an appealing strategy to improve cancer therapy. Opportunities for drug repurposing are often based on occasional observations or on time-consuming pre-clinical drug screenings that are often not hypothesis-driven. In contrast, in-silico drug repurposing is an emerging, hypothesis-driven approach that takes advantage of the use of big-data. Indeed, the extensive use of -omics technologies, improved data storage, data meaning, machine learning algorithms, and computational modeling all offer unprecedented knowledge of the biological mechanisms of cancers and drugs' modes of action, providing extensive availability for both disease-related data and drugs-related data. This offers the opportunity to generate, with time and cost-effective approaches, computational drug networks to predict, in-silico, the efficacy of approved drugs against relevant cancer targets, as well as to select better responder patients or disease' biomarkers. Here, we will review selected disease-related data together with computational tools to be exploited for the in-silico repurposing of drugs against validated targets in cancer therapies, focusing on the oncogenic signaling pathways activation in cancer. We will discuss how in-silico drug repurposing has the promise to shortly improve our arsenal of anticancer drugs and, likely, overcome certain limitations of modern cancer therapies against old and new therapeutic targets in oncology.

Predictive Signatures Inform the Effective Repurposing of Decitabine to Treat KRAS-Dependent Pancreatic Ductal Adenocarcinoma

Mutated KRAS protein is a pivotal tumor driver in pancreatic cancer. However, despite comprehensive efforts, effective therapeutics that can target oncogenic KRAS are still under investigation or awaiting clinical approval. Using a specific KRAS-dependent gene signature, we implemented a computer-assisted inspection of a drug-gene network to in silico repurpose drugs that work like inhibitors of oncogenic KRAS. We identified and validated decitabine, an FDA-approved drug, as a potent inhibitor of growth in pancreatic cancer cells and patient-derived xenograft models that showed KRAS dependency. Mechanistically, decitabine efficacy was linked to KRAS-driven dependency on nucleotide metabolism and its ability to specifically impair pyrimidine biosynthesis in KRAS-dependent tumors cells. These findings also showed that gene signatures related to KRAS dependency might be prospectively used to inform on decitabine sensitivity in a selected subset of patients with KRAS-mutated pancreatic cancer. Overall, the repurposing of decitabine emerged as an intriguing option for treating pancreatic tumors that are addicted to mutant KRAS, thus offering opportunities for improving the arsenal of therapeutics for this extremely deadly disease. SIGNIFICANCE: Decitabine is a promising drug for cancer cells dependent on RAS signaling.

Subcellular compartmentalization of NAD+ and its role in cancer: A sereNADe of metabolic melodies

Nicotinamide adenine dinucleotide (NAD+) is an essential biomolecule involved in many critical processes. Its role as both a driver of energy production and a signaling molecule underscores its importance in health and disease. NAD+ signaling impacts multiple processes that are dysregulated in cancer, including DNA repair, cell proliferation, differentiation, redox regulation, and oxidative stress. Distribution of NAD+ is highly compartmentalized, with each subcellular NAD+ pool differentially regulated and preferentially involved in distinct NAD+-dependent signaling or metabolic events. Emerging evidence suggests that targeting NAD+ metabolism is likely to repress many specific mechanisms underlying tumor development and progression, including proliferation, survival, metabolic adaptations, invasive capabilities, heterotypic interactions with the tumor microenvironment, and stress response including notably DNA maintenance and repair. Here we provide a comprehensive overview of how compartmentalized NAD+ metabolism in mitochondria, nucleus, cytosol, and extracellular space impacts cancer formation and progression, along with a discussion of the therapeutic potential of NAD+-targeting drugs in cancer.

Metabolic Regulation of Redox Balance in Cancer

Reactive oxygen species (ROS) are chemically active free radicals produced by partial reduction of oxygen that can activate discrete signaling pathways or disrupt redox homeostasis depending on their concentration. ROS interacts with biomolecules, including DNA, and can cause mutations that can transform normal cells into cancer cells. Furthermore, certain cancer-causing mutations trigger alterations in cellular metabolism that can increase ROS production, resulting in genomic instability, additional DNA mutations, and tumor evolution. To prevent excess ROS-mediated toxicity, cancer-causing mutations concurrently activate pathways that manage this oxidative burden. Hence, an understanding of the metabolic pathways that regulate ROS levels is imperative for devising therapies that target tumor cells. In this review, we summarize the dual role of metabolism as a generator and inhibitor of ROS in cancer and discuss current strategies to target the ROS axis.

Functional and mechanistic studies reveal MAGEA3 as a pro-survival factor in pancreatic cancer cells

Background

In the era of personalized therapy, functional annotation of less frequent genetic aberrations will be instrumental in adapting effective therapeutic in clinic. Overexpression of Melanoma associated antigen A3 (MAGEA3) is reported in certain pancreatic cancer (PCA) patients. The major objective of the current study was to investigate the functional role of MAGEA3 in pancreatic cancer cells (PCCs) growth and survival.

Methods

Using overexpression (tet-on regulated system and constitutive expression system) and knockdown (by siRNA and shRNA) approach, we dissected the mechanistic role of MAGEA3 in pancreatic cancer pathogenesis. We generated MAGEA3 expressing stable PCA cell lines and mouse primary pancreatic epithelial cells. MAGEA3 was also depleted in certain MAGEA3 positive PCCs by siRNA or shRNA. The stable cells were subjected to in vitro assays like proliferation and survival assays under growth factor deprivation or in the presence of cytotoxic drugs. The MAGEA3 overexpressing or depleted stable PCCs were evaluated in vivo using xenograft model to check the role of MAGEA3 in tumor progression. We also dissected the mechanism behind the MAGEA3 role in tumor progression using western blot analysis and CCL2 neutralization.

Results

MAGEA3 overexpression in PCA cells did not alter the cell proliferation but protected the cells during growth factor deprivation and also in the presence of cytotoxic drugs. However, depletion of MAGEA3 in MAGEA3 positive cells resulted in reduced cell proliferation and increased apoptosis upon growth factor deprivation and also in response to cytotoxic drugs. The in vivo xenograft study revealed that overexpression of MAGEA3 promoted tumor growth however depleting the same hindered the tumor progression. Mechanistically, our in vitro and in vivo study revealed that MAGEA3 has tumor-promoting role by reducing macro-autophagy and overexpressing pro-survival molecules like CCL2 and survivin.

Conclusion

Our data proves tumor-promoting role of MAGEA3 and provides the rationale to target MAGEA3 and/or its functional mediators like CCL2 for PCA, which may have a better impact in PCA therapy.

Electronic supplementary material

The online version of this article (10.1186/s13046-019-1272-2) contains supplementary material, which is available to authorized users.

Exploiting rare driver mutations for precision cancer medicine

Catalyzed by the ability to develop precision therapies targeting the unique genetic changes that drive individual tumors, sequencing patients' tumor genomes is an increasingly common practice in oncology. In most cancer types, however, a limited number of common mutations are accompanied by a plethora of low-frequency variants whose functional consequences and clinical actionability are often unknown. We here illustrate that this 'long tail' of infrequent molecular alterations includes oncogenic drivers of biological significance that can be the genetic basis of extraordinary responses to systemic cancer therapies. Furthermore, we review current strategies to identify, prioritize, and experimentally validate novel long-tail driver mutations, efforts that will likely provide new insights into the clinically actionable genome and improve outcomes for patients.

NADP+ is an endogenous PARP inhibitor in DNA damage response and tumor suppression

ADP-ribosylation is a unique posttranslational modification catalyzed by poly(ADP-ribose) polymerases (PARPs) using NAD⁺ as ADP-ribose donor. PARPs play an indispensable role in DNA damage repair and small molecule PARP inhibitors have emerged as potent anticancer drugs. However, to date, PARP inhibitor treatment has been restricted to patients with BRCA1/2 mutation-associated breast and ovarian cancer. One of the major challenges to extend the therapeutic potential of PARP inhibitors to other cancer types is the absence of predictive biomarkers. Here, we show that ovarian cancer cells with higher level of NADP⁺, an NAD⁺ derivative, are more sensitive to PARP inhibitors. We demonstrate that NADP⁺ acts as a negative regulator and suppresses ADP-ribosylation both in vitro and in vivo. NADP⁺ impairs ADP-ribosylation-dependent DNA damage repair and sensitizes tumor cell to chemically synthesized PARP inhibitors. Taken together, our study identifies NADP⁺ as an endogenous PARP inhibitor that may have implications in cancer treatment.

Cancer cells respond differently to inhibitors of Poly (ADP-ribose) polymerase. Here the authors reveal that ovarian cancer cells with higher cellular NADP+ levels are more sensitive to clinically relevant PARP1 inhibitors and show that NADP+ act as an endogenous inhibitor of PARP enzymes.

Metabolic Dependencies in Pancreatic Cancer

Pancreatic ductal adenocarcinoma (PDA) is a highly lethal cancer with a long-term survival rate under 10%. Available cytotoxic chemotherapies have significant side effects, and only marginal therapeutic efficacy. FDA approved drugs currently used against PDA target DNA metabolism and DNA integrity. However, alternative metabolic targets beyond DNA may prove to be much more effective. PDA cells are forced to live within a particularly severe microenvironment characterized by relative hypovascularity, hypoxia, and nutrient deprivation. Thus, PDA cells must possess biochemical flexibility in order to adapt to austere conditions. A better understanding of the metabolic dependencies required by PDA to survive and thrive within a harsh metabolic milieu could reveal specific metabolic vulnerabilities. These molecular requirements can then be targeted therapeutically, and would likely be associated with a clinically significant therapeutic window since the normal tissue is so well-perfused with an abundant nutrient supply. Recent work has uncovered a number of promising therapeutic targets in the metabolic domain, and clinicians are already translating some of these discoveries to the clinic. In this review, we highlight mitochondria metabolism, non-canonical nutrient acquisition pathways (macropinocytosis and use of pancreatic stellate cell-derived alanine), and redox homeostasis as compelling therapeutic opportunities in the metabolic domain.

Comprehensive Molecular Characterization of the Hippo Signaling Pathway in Cancer

Hippo signaling has been recognized as a key tumor suppressor pathway. Here, we perform a comprehensive molecular characterization of 19 Hippo core genes in 9,125 tumor samples across 33 cancer types using multidimensional "omic" data from The Cancer Genome Atlas. We identify somatic drivers among Hippo genes and the related microRNA (miRNA) regulators, and using functional genomic approaches, we experimentally characterize YAP and TAZ mutation effects and miR-590 and miR-200a regulation for TAZ. Hippo pathway activity is best characterized by a YAP/TAZ transcriptional target signature of 22 genes, which shows robust prognostic power across cancer types. Our elastic-net integrated modeling further reveals cancer-type-specific pathway regulators and associated cancer drivers. Our results highlight the importance of Hippo signaling in squamous cell cancers, characterized by frequent amplification of YAP/TAZ, high expression heterogeneity, and significant prognostic patterns. This study represents a systems-biology approach to characterizing key cancer signaling pathways in the post-genomic era.

TMEM106B drives lung cancer metastasis by inducing TFEB-dependent lysosome synthesis and secretion of cathepsins

Metastatic lung cancer is the leading cause of cancer-associated mortality worldwide, therefore necessitating novel approaches to identify specific genetic drivers for lung cancer progression and metastasis. We recently performed an in vivo gain-of-function genetic screen to identify driver genes of lung cancer metastasis. In the study reported here, we identify TMEM106B as a primary robust driver of lung cancer metastasis. Ectopic expression of TMEM106B could significantly promote the synthesis of enlarged vesicular lysosomes that are laden with elevated levels of active cathepsins. In a TFEB-dependent manner, TMEM106B could modulate the expression of lysosomal genes of the coordinated lysosomal expression and regulation (CLEAR) pathway in lung cancer cells and patient samples. We also demonstrate that TMEM106B-induced lysosomes undergo calcium-dependent exocytosis, thereby releasing active lysosomal cathepsins necessary for TMEM106B-mediated cancer cell invasion and metastasis in vivo, which could be therapeutically prevented by pharmacological inhibition of cathepsins. Further, in TCGA LUAD data sets, 19% of patients show elevated expression of TMEM106B, which predicts for poor disease-free and overall-survival.

In vivo screening identifies GATAD2B as a metastasis driver in KRAS-driven lung cancer

Genetic aberrations driving pro-oncogenic and pro-metastatic activity remain an elusive target in the quest of precision oncology. To identify such drivers, we use an animal model of KRAS-mutant lung adenocarcinoma to perform an in vivo functional screen of 217 genetic aberrations selected from lung cancer genomics datasets. We identify 28 genes whose expression promoted tumor metastasis to the lung in mice. We employ two tools for examining the KRAS-dependence of genes identified from our screen: 1) a human lung cell model containing a regulatable mutant KRAS allele and 2) a lentiviral system permitting co-expression of DNA-barcoded cDNAs with Cre recombinase to activate a mutant KRAS allele in the lungs of mice. Mechanistic evaluation of one gene, GATAD2B, illuminates its role as a dual activity gene, promoting both pro-tumorigenic and pro-metastatic activities in KRAS-mutant lung cancer through interaction with c-MYC and hyperactivation of the c-MYC pathway.

Systematic Functional Annotation of Somatic Mutations in Cancer

The functional impact of the vast majority of cancer somatic mutations remains unknown, representing a critical knowledge gap for implementing precision oncology. Here, we report the development of a moderate-throughput functional genomic platform consisting of efficient mutant generation, sensitive viability assays using two growth factor-dependent cell models, and functional proteomic profiling of signaling effects for select aberrations. We apply the platform to annotate >1,000 genomic aberrations, including gene amplifications, point mutations, indels, and gene fusions, potentially doubling the number of driver mutations characterized in clinically actionable genes. Further, the platform is sufficiently sensitive to identify weak drivers. Our data are accessible through a user-friendly, public data portal. Our study will facilitate biomarker discovery, prediction algorithm improvement, and drug development.

Genome-Wide CRISPR Screen for Essential Cell Growth Mediators in Mutant KRAS Colorectal Cancers

Targeting mutant KRAS signaling pathways continues to attract attention as a therapeutic strategy for KRAS-driven tumors. In this study, we exploited the power of the CRISPR-Cas9 system to identify genes affecting the tumor xenograft growth of human mutant KRAS (KRAS^MUT) colorectal cancers. Using pooled lentiviral single-guide RNA libraries, we conducted a genome-wide loss-of-function genetic screen in an isogenic pair of human colorectal cancer cell lines harboring mutant or wild-type KRAS. The screen identified novel and established synthetic enhancers or synthetic lethals for KRAS^MUT colorectal cancer, including targetable metabolic genes. Notably, genetic disruption or pharmacologic inhibition of the metabolic enzymes NAD kinase or ketohexokinase was growth inhibitory in vivo In addition, the chromatin remodeling protein INO80C was identified as a novel tumor suppressor in KRAS^MUT colorectal and pancreatic tumor xenografts. Our findings define a novel targetable set of therapeutic targets for KRAS^MUT tumors. Cancer Res; 77(22); 6330-9. ©2017 AACR.

Rational combination therapy with PARP and MEK inhibitors capitalizes on therapeutic liabilities in RAS mutant cancers

Mutant RAS has remained recalcitrant to targeted therapy efforts. We demonstrate that combined treatment with poly(adenosine diphosphate-ribose) polymerase (PARP) inhibitors and mitogen-activated protein kinase (MAPK) kinase (MEK) inhibitors evokes unanticipated, synergistic cytotoxic effects in vitro and in vivo in multiple RAS mutant tumor models across tumor lineages where RAS mutations are prevalent. The effects of PARP and MEK inhibitor combinations are independent of BRCA1/2 and p53 mutation status, suggesting that the synergistic activity is likely to be generalizable. Synergistic activity of PARP and MEK inhibitor combinations in RAS mutant tumors is associated with (i) induction of BIM-mediated apoptosis, (ii) decrease in expression of components of the homologous recombination DNA repair pathway, (iii) decrease in homologous recombination DNA damage repair capacity, (iv) decrease in DNA damage checkpoint activity, (v) increase in PARP inhibitor-induced DNA damage, (vi) decrease in vascularity that could increase PARP inhibitor efficacy by inducing hypoxia, and (vii) elevated PARP1 protein, which increases trapping activity of PARP inhibitors. Mechanistically, enforced expression of FOXO3a, which is a target of the RAS/MAPK pathway, was sufficient to recapitulate the functional consequences of MEK inhibitors including synergy with PARP inhibitors. Thus, the ability of mutant RAS to suppress FOXO3a and its reversal by MEK inhibitors accounts, at least in part, for the synergy of PARP and MEK inhibitors in RAS mutant tumors. The rational combination of PARP and MEK inhibitors warrants clinical investigation in patients with RAS mutant tumors where there are few effective therapeutic options.

Engineering and Functional Characterization of Fusion Genes Identifies Novel Oncogenic Drivers of Cancer

Oncogenic gene fusions drive many human cancers, but tools to more quickly unravel their functional contributions are needed. Here we describe methodology permitting fusion gene construction for functional evaluation. Using this strategy, we engineered the known fusion oncogenes, BCR-ABL1, EML4-ALK, and ETV6-NTRK3, as well as 20 previously uncharacterized fusion genes identified in The Cancer Genome Atlas datasets. In addition to confirming oncogenic activity of the known fusion oncogenes engineered by our construction strategy, we validated five novel fusion genes involving MET, NTRK2, and BRAF kinases that exhibited potent transforming activity and conferred sensitivity to FDA-approved kinase inhibitors. Our fusion construction strategy also enabled domain-function studies of BRAF fusion genes. Our results confirmed other reports that the transforming activity of BRAF fusions results from truncation-mediated loss of inhibitory domains within the N-terminus of the BRAF protein. BRAF mutations residing within this inhibitory region may provide a means for BRAF activation in cancer, therefore we leveraged the modular design of our fusion gene construction methodology to screen N-terminal domain mutations discovered in tumors that are wild-type at the BRAF mutation hotspot, V600. We identified an oncogenic mutation, F247L, whose expression robustly activated the MAPK pathway and sensitized cells to BRAF and MEK inhibitors. When applied broadly, these tools will facilitate rapid fusion gene construction for subsequent functional characterization and translation into personalized treatment strategies. Cancer Res; 77(13); 3502-12. ©2017 AACR.

Inhibition of DNA2 nuclease as a therapeutic strategy targeting replication stress in cancer cells

Replication stress is a characteristic feature of cancer cells, which is resulted from sustained proliferative signaling induced by activation of oncogenes or loss of tumor suppressors. In cancer cells, oncogene-induced replication stress manifests as replication-associated lesions, predominantly double-strand DNA breaks (DSBs). An essential mechanism utilized by cells to repair replication-associated DSBs is homologous recombination (HR). In order to overcome replication stress and survive, cancer cells often require enhanced HR repair capacity. Therefore, the key link between HR repair and cellular tolerance to replication-associated DSBs provides us with a mechanistic rationale for exploiting synthetic lethality between HR repair inhibition and replication stress. DNA2 nuclease is an evolutionarily conserved essential enzyme in replication and HR repair. Here we demonstrate that DNA2 is overexpressed in pancreatic cancers, one of the deadliest and more aggressive forms of human cancers, where mutations in the KRAS are present in 90–95% of cases. In addition, depletion of DNA2 significantly reduces pancreatic cancer cell survival and xenograft tumor growth, suggesting the therapeutic potential of DNA2 inhibition. Finally, we develop a robust high-throughput biochemistry assay to screen for inhibitors of the DNA2 nuclease activity. The top inhibitors were shown to be efficacious against both yeast Dna2 and human DNA2. Treatment of cancer cells with DNA2 inhibitors recapitulates phenotypes observed upon DNA2 depletion, including decreased DNA double strand break end resection and attenuation of HR repair. Similar to genetic ablation of DNA2, chemical inhibition of DNA2 selectively attenuates the growth of various cancer cells with oncogene-induced replication stress. Taken together, our findings open a new avenue to develop a new class of anticancer drugs by targeting druggable nuclease DNA2. We propose DNA2 inhibition as new strategy in cancer therapy by targeting replication stress, a molecular property of cancer cells that is acquired as a result of oncogene activation instead of targeting currently undruggable oncoprotein itself such as KRAS.

Gene Essentiality Profiling Reveals Gene Networks and Synthetic Lethal Interactions with Oncogenic Ras

The genetic dependencies of human cancers widely vary. Here, we catalog this heterogeneity and use it to identify functional gene interactions and genotype-dependent liabilities in cancer. By using genome-wide CRISPR-based screens, we generate a gene essentiality dataset across 14 human acute myeloid leukemia (AML) cell lines. Sets of genes with correlated patterns of essentiality across the lines reveal new gene relationships, the essential substrates of enzymes, and the molecular functions of uncharacterized proteins. Comparisons of differentially essential genes between Ras-dependent and -independent lines uncover synthetic lethal partners of oncogenic Ras. Screens in both human AML and engineered mouse pro-B cells converge on a surprisingly small number of genes in the Ras processing and MAPK pathways and pinpoint PREX1 as an AML-specific activator of MAPK signaling. Our findings suggest general strategies for defining mammalian gene networks and synthetic lethal interactions by exploiting the natural genetic and epigenetic diversity of human cancer cells.