Bacteria-to-Human Protein Networks Reveal Origins of Endogenous DNA Damage.

DNA damage provokes mutations and cancer and results from external carcinogens or endogenous cellular processes. However, the intrinsic instigators of endogenous DNA damage are poorly understood. Here, we identify proteins that promote endogenous DNA damage when overproduced: the DNA "damage-up" proteins (DDPs). We discover a large network of DDPs in Escherichia coli and deconvolute them into six function clusters, demonstrating DDP mechanisms in three: reactive oxygen increase by transmembrane transporters, chromosome loss by replisome binding, and replication stalling by transcription factors. Their 284 human homologs are over-represented among known cancer drivers, and their RNAs in tumors predict heavy mutagenesis and a poor prognosis. Half of the tested human homologs promote DNA damage and mutation when overproduced in human cells, with DNA damage-elevating mechanisms like those in E. coli. Our work identifies networks of DDPs that provoke endogenous DNA damage and may reveal DNA damage-associated functions of many human known and newly implicated cancer-promoting proteins.

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.

Literature-based automated discovery of tumor suppressor p53 phosphorylation and inhibition by NEK2.

Scientific progress depends on formulating testable hypotheses informed by the literature. In many domains, however, this model is strained because the number of research papers exceeds human readability. Here, we developed computational assistance to analyze the biomedical literature by reading PubMed abstracts to suggest new hypotheses. The approach was tested experimentally on the tumor suppressor p53 by ranking its most likely kinases, based on all available abstracts. Many of the best-ranked kinases were found to bind and phosphorylate p53 (P value = 0.005), suggesting six likely p53 kinases so far. One of these, NEK2, was studied in detail. A known mitosis promoter, NEK2 was shown to phosphorylate p53 at Ser315 in vitro and in vivo and to functionally inhibit p53. These bona fide validations of text-based predictions of p53 phosphorylation, and the discovery of an inhibitory p53 kinase of pharmaceutical interest, suggest that automated reasoning using a large body of literature can generate valuable molecular hypotheses and has the potential to accelerate scientific discovery.

Systematic characterization of A-to-I RNA editing hotspots in microRNAs across human cancers.

RNA editing, a widespread post-transcriptional mechanism, has emerged as a new player in cancer biology. Recent studies have reported key roles for individual miRNA editing events, but a comprehensive picture of miRNA editing in human cancers remains largely unexplored. Here, we systematically characterized the miRNA editing profiles of 8595 samples across 20 cancer types from miRNA sequencing data of The Cancer Genome Atlas and identified 19 adenosine-to-inosine (A-to-I) RNA editing hotspots. We independently validated 15 of them by perturbation experiments in several cancer cell lines. These miRNA editing events show extensive correlations with key clinical variables (e.g., tumor subtype, disease stage, and patient survival time) and other molecular drivers. Focusing on the RNA editing hotspot in miR-200b, a key tumor metastasis suppressor, we found that the miR-200b editing level correlates with patient prognosis opposite to the pattern observed for the wild-type miR-200b expression. We further experimentally showed that, in contrast to wild-type miRNA, the edited miR-200b can promote cell invasion and migration through its impaired ability to inhibit ZEB1/ZEB2 and acquired concomitant ability to repress new targets, including LIFR, a well-characterized metastasis suppressor. Our study highlights the importance of miRNA editing in gene regulation and suggests its potential as a biomarker for cancer prognosis and therapy.

Analyzing Somatic Genome Rearrangements in Human Cancers by Using Whole-Exome Sequencing.

Although exome sequencing data are generated primarily to detect single-nucleotide variants and indels, they can also be used to identify a subset of genomic rearrangements whose breakpoints are located in or near exons. Using >4,600 tumor and normal pairs across 15 cancer types, we identified over 9,000 high confidence somatic rearrangements, including a large number of gene fusions. We find that the 5' fusion partners of functional fusions are often housekeeping genes, whereas the 3' fusion partners are enriched in tyrosine kinases. We establish the oncogenic potential of ROR1-DNAJC6 and CEP85L-ROS1 fusions by showing that they can promote cell proliferation in vitro and tumor formation in vivo. Furthermore, we found that ∼4% of the samples have massively rearranged chromosomes, many of which are associated with upregulation of oncogenes such as ERBB2 and TERT. Although the sensitivity of detecting structural alterations from exomes is considerably lower than that from whole genomes, this approach will be fruitful for the multitude of exomes that have been and will be generated, both in cancer and in other diseases.

Intronic miR-211 assumes the tumor suppressive function of its host gene in melanoma.

When it escapes early detection, malignant melanoma becomes a highly lethal and treatment-refractory cancer. Melastatin is greatly downregulated in metastatic melanomas and is widely believed to function as a melanoma tumor suppressor. Here we report that tumor suppressive activity is not mediated by melastatin but instead by a microRNA (miR-211) hosted within an intron of melastatin. Increasing expression of miR-211 but not melastatin reduced migration and invasion of malignant and highly invasive human melanomas characterized by low levels of melastatin and miR-211. An unbiased network analysis of melanoma-expressed genes filtered for their roles in metastasis identified three central node genes: IGF2R, TGFBR2, and NFAT5. Expression of these genes was reduced by miR-211, and knockdown of each gene phenocopied the effects of increased miR-211 on melanoma invasiveness. These data implicate miR-211 as a suppressor of melanoma invasion whose expression is silenced or selected against via suppression of the entire melastatin locus during human melanoma progression.

Differential UBE2C and HOXA1 expression in melanocytic nevi and melanoma.

BACKGROUND: Recent molecular advances suggest that Spitz nevi and other spitzoid neoplasms are biologically distinct from melanoma and conventional nevi. The ubiquitin ligase UBE2C and the homeobox transcription factor HOXA1 are candidate oncogenes in melanoma. METHODS: Using RNA expression analysis and immunohistochemistry, we evaluated these biomarkers in Spitz nevi (n = 20), melanoma (n = 20), and by immunohistochemistry in conventional nevi (n = 20). RESULTS: RNA analysis with branched DNA multiplex assay identified upregulation of UBE2C in melanomas vs Spitz nevi (P = .003), whereas HOXA1 was downregulated in melanoma (P < .0001). Immunohistochemical analysis confirmed increased nuclear expression of UBE2C in melanoma (mean = 18% of cells; range 3%-44%) when compared with Spitz nevi (mean = 9%; range 2%-28%; P = .001) and conventional nevi (mean = 1.5%; range 0-9%; P < .0001). Strong UBE2C staining was identified in cells undergoing mitosis. UBE2C RNA and protein detection correlated with mitotic rate (P < .0001). On the other hand, HOXA1 nuclear staining was low in melanoma (mean = 69%; range 5%-100%) when compared with Spitz nevi (mean = 94%; range 66%-100%; P = .0024) and conventional nevi (mean = 94%; range 83%-99%; P = .009). CONCLUSIONS: UBE2C and HOXA1 RNA and protein are differentially expressed in conventional and Spitz nevi and melanoma.

SMAD4-dependent barrier constrains prostate cancer growth and metastatic progression.

Effective clinical management of prostate cancer (PCA) has been challenged by significant intratumoural heterogeneity on the genomic and pathological levels and limited understanding of the genetic elements governing disease progression. Here, we exploited the experimental merits of the mouse to test the hypothesis that pathways constraining progression might be activated in indolent Pten-null mouse prostate tumours and that inactivation of such progression barriers in mice would engender a metastasis-prone condition. Comparative transcriptomic and canonical pathway analyses, followed by biochemical confirmation, of normal prostate epithelium versus poorly progressive Pten-null prostate cancers revealed robust activation of the TGFß/BMP-SMAD4 signalling axis. The functional relevance of SMAD4 was further supported by emergence of invasive, metastatic and lethal prostate cancers with 100% penetrance upon genetic deletion of Smad4 in the Pten-null mouse prostate. Pathological and molecular analysis as well as transcriptomic knowledge-based pathway profiling of emerging tumours identified cell proliferation and invasion as two cardinal tumour biological features in the metastatic Smad4/Pten-null PCA model. Follow-on pathological and functional assessment confirmed cyclin D1 and SPP1 as key mediators of these biological processes, which together with PTEN and SMAD4, form a four-gene signature that is prognostic of prostate-specific antigen (PSA) biochemical recurrence and lethal metastasis in human PCA. This model-informed progression analysis, together with genetic, functional and translational studies, establishes SMAD4 as a key regulator of PCA progression in mice and humans.

Mutations in FBXL4 cause mitochondrial encephalopathy and a disorder of mitochondrial DNA maintenance.

Nuclear genetic disorders causing mitochondrial DNA (mtDNA) depletion are clinically and genetically heterogeneous, and the molecular etiology remains undiagnosed in the majority of cases. Through whole-exome sequencing, we identified recessive nonsense and splicing mutations in FBXL4 segregating in three unrelated consanguineous kindreds in which affected children present with a fatal encephalopathy, lactic acidosis, and severe mtDNA depletion in muscle. We show that FBXL4 is an F-box protein that colocalizes with mitochondria and that loss-of-function and splice mutations in this protein result in a severe respiratory chain deficiency, loss of mitochondrial membrane potential, and a disturbance of the dynamic mitochondrial network and nucleoid distribution in fibroblasts from affected individuals. Expression of the wild-type FBXL4 transcript in cell lines from two subjects fully rescued the levels of mtDNA copy number, leading to a correction of the mitochondrial biochemical deficit. Together our data demonstrate that mutations in FBXL4 are disease causing and establish FBXL4 as a mitochondrial protein with a possible role in maintaining mtDNA integrity and stability.

Whole-exome sequencing combined with functional genomics reveals novel candidate driver cancer genes in endometrial cancer.

Endometrial cancer is the most common gynecological malignancy, with more than 280,000 cases occurring annually worldwide. Although previous studies have identified important common somatic mutations in endometrial cancer, they have primarily focused on a small set of known cancer genes and have thus provided a limited view of the molecular basis underlying this disease. Here we have developed an integrated systems-biology approach to identifying novel cancer genes contributing to endometrial tumorigenesis. We first performed whole-exome sequencing on 13 endometrial cancers and matched normal samples, systematically identifying somatic alterations with high precision and sensitivity. We then combined bioinformatics prioritization with high-throughput screening (including both shRNA-mediated knockdown and expression of wild-type and mutant constructs) in a highly sensitive cell viability assay. Our results revealed 12 potential driver cancer genes including 10 tumor-suppressor candidates (ARID1A, INHBA, KMO, TTLL5, GRM8, IGFBP3, AKTIP, PHKA2, TRPS1, and WNT11) and two oncogene candidates (ERBB3 and RPS6KC1). The results in the "sensor" cell line were recapitulated by siRNA-mediated knockdown in endometrial cancer cell lines. Focusing on ARID1A, we integrated mutation profiles with functional proteomics in 222 endometrial cancer samples, demonstrating that ARID1A mutations frequently co-occur with mutations in the phosphatidylinositol 3-kinase (PI3K) pathway and are associated with PI3K pathway activation. siRNA knockdown in endometrial cancer cell lines increased AKT phosphorylation supporting ARID1A as a novel regulator of PI3K pathway activity. Our study presents the first unbiased view of somatic coding mutations in endometrial cancer and provides functional evidence for diverse driver genes and mutations in this disease.

Circulating tumor cells from 4D model have less integrin beta 4 expression.

BACKGROUND: Currently, there is no in vitro or ex vivo model that can isolate circulating tumor cells (CTCs). Recently, we developed a four-dimensional (4D) lung cancer model that allows for the isolation of CTCs. We postulated that these cells have different properties than parental (2D) cells. MATERIALS AND METHODS: We obtained CTCs by growing A549, H1299, 393P, and 344SQ cell lines on the 4D lung model. The CTCs were functionally characterized in vitro and gene expression of the cell adhesion molecules was compared with respective 2D cells. Integrin beta 4 (ITGB4) was further investigated by stably transfecting the A549 and H1299 cells. RESULTS: We found that all cell lines produced CTCs, and that CTCs from the 4D model were less adherent to the plastic and have a slower growth rate than respective 2D cells (P < 0.01). Most of the cell adhesion molecules were downregulated (P < 0.05) in CTCs, and ITGB4 was the common molecule, significantly more underexpressed in CTCs from all cell lines than their respective 2D cells. The modulation of ITGB4 led to a differential function of 2D cells. CONCLUSIONS: CTCs from the 4D model have different transcriptional, translational, and in vitro characteristics than the same cells grown on a petri dish, and these CTCs from the 4D model have the properties of CTCs that are responsible for metastasis.

GOLPH3 modulates mTOR signalling and rapamycin sensitivity in cancer.

Genome-wide copy number analyses of human cancers identified a frequent 5p13 amplification in several solid tumour types, including lung (56%), ovarian (38%), breast (32%), prostate (37%) and melanoma (32%). Here, using integrative analysis of a genomic profile of the region, we identify a Golgi protein, GOLPH3, as a candidate targeted for amplification. Gain- and loss-of-function studies in vitro and in vivo validated GOLPH3 as a potent oncogene. Physically, GOLPH3 localizes to the trans-Golgi network and interacts with components of the retromer complex, which in yeast has been linked to target of rapamycin (TOR) signalling. Mechanistically, GOLPH3 regulates cell size, enhances growth-factor-induced mTOR (also known as FRAP1) signalling in human cancer cells, and alters the response to an mTOR inhibitor in vivo. Thus, genomic and genetic, biological, functional and biochemical data in yeast and humans establishes GOLPH3 as a new oncogene that is commonly targeted for amplification in human cancer, and is capable of modulating the response to rapamycin, a cancer drug in clinical use.

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.

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.

A functional genomic approach to actionable gene fusions for precision oncology.

Fusion genes represent a class of attractive therapeutic targets. Thousands of fusion genes have been identified in patients with cancer, but the functional consequences and therapeutic implications of most of these remain largely unknown. Here, we develop a functional genomic approach that consists of efficient fusion reconstruction and sensitive cell viability and drug response assays. Applying this approach, we characterize ~100 fusion genes detected in patient samples of The Cancer Genome Atlas, revealing a notable fraction of low-frequency fusions with activating effects on tumor growth. Focusing on those in the RTK-RAS pathway, we identify a number of activating fusions that can markedly affect sensitivity to relevant drugs. Last, we propose an integrated, level-of-evidence classification system to prioritize gene fusions systematically. Our study reiterates the urgent clinical need to incorporate similar functional genomic approaches to characterize gene fusions, thereby maximizing the utility of gene fusions for precision oncology.

mTORC1-dependent and -independent regulation of stem cell renewal, differentiation, and mobilization.

The Tuberous Sclerosis Complex component, TSC1, functions as a tumor suppressor via its regulation of diverse cellular processes, particularly cell growth. TSC1 exists in a complex with TSC2 and functions primarily as a key negative regulator of mammalian target of rapamycin complex 1 (mTORC1) signaling and protein synthesis, although the TSC1/TSC2 complex also shows mTORC1-independent outputs to other pathways. Here, we explored the role of TSC1 in various aspects of stem cell biology and dissected the extent to which TSC1 functions are executed via mTORC1-dependent versus mTORC1-independent pathways. Using hematopoietic stem cells (HSCs) as a model system, we demonstrate that somatic deletion of TSC1 produces striking stem cell and derivative effector cell phenotypes characterized by increased HSC cell cycling, mobilization, marked progressive depletion, defective long-term repopulating potential, and hematopoietic lineage developmental aberrations. On the mechanistic level, we further establish that TSC1 regulation of HSC quiescence and long-term repopulating potential and hematopoietic lineage development is mediated through mTORC1 signaling. In contrast, TSC1 regulation of HSC mobilization is effected in an mTORC1-independent manner, and gene profiling and functional analyses reveals the actin-bundling protein FSCN1 as a key TSC1/TSC2 target in the regulation of HSC mobilization. Thus, TSC1 is a critical regulator of HSC self-renewal, mobilization, and multilineage development and executes these actions via both mTORC1-dependent and -independent pathways.

POT1 Regulates Proliferation and Confers Sexual Dimorphism in Glioma.

Germline POT1 mutations are found in a spectrum of cancers and confer increased risk. Recently, we identified a series of novel germline POT1 mutations that predispose carrier families to the development of glioma. Despite these strong associations, how these glioma-associated POT1 mutations contribute to glioma tumorigenesis remains undefined. Here we show that POT1-G95C increases proliferation in glioma-initiating cells in vitro and in progenitor populations in the developing brain. In a native mouse model of glioma, loss of Pot1a/b resulted in decreased survival in females compared with males. These findings were corroborated in human glioma, where low POT1 expression correlated with decreased survival in females. Transcriptomic and IHC profiling of Pot1a/b-deficient glioma revealed that tumors in females exhibited decreased expression of immune markers and increased expression of cell-cycle signatures. Similar sex-dependent trends were observed in human gliomas that had low expression of POT1. Together, our studies demonstrate context-dependent functions for POT1 mutation or loss in driving progenitor proliferation in the developing brain and sexual dimorphism in glioma. SIGNIFICANCE: This study shows that manipulation of POT1 expression in glioma has sex-specific effects on tumorigenesis and associated immune signatures.

IMPAD1 and KDELR2 drive invasion and metastasis by enhancing Golgi-mediated secretion.

Non-small cell lung cancer (NSCLC) is the deadliest form of cancer worldwide, due in part to its proclivity to metastasize. Identifying novel drivers of invasion and metastasis holds therapeutic potential for the disease. We conducted a gain-of-function invasion screen, which identified two separate hits, IMPAD1 and KDELR2, as robust, independent drivers of lung cancer invasion and metastasis. Given that IMPAD1 and KDELR2 are known to be localized to the ER-Golgi pathway, we studied their common mechanism of driving in vitro invasion and in vivo metastasis and demonstrated that they enhance Golgi-mediated function and secretion. Therapeutically inhibiting matrix metalloproteases (MMPs) suppressed both IMPAD1- and KDELR2-mediated invasion. The hits from this unbiased screen and the mechanistic validation highlight Golgi function as one of the key cellular features altered during invasion and metastasis.

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.

The RNA-binding protein AKAP8 suppresses tumor metastasis by antagonizing EMT-associated alternative splicing.

Alternative splicing has been shown to causally contribute to the epithelial-mesenchymal transition (EMT) and tumor metastasis. However, the scope of splicing factors that govern alternative splicing in these processes remains largely unexplored. Here we report the identification of A-Kinase Anchor Protein (AKAP8) as a splicing regulatory factor that impedes EMT and breast cancer metastasis. AKAP8 not only is capable of inhibiting splicing activity of the EMT-promoting splicing regulator hnRNPM through protein-protein interaction, it also directly binds to RNA and alters splicing outcomes. Genome-wide analysis shows that AKAP8 promotes an epithelial cell state splicing program. Experimental manipulation of an AKAP8 splicing target CLSTN1 revealed that splice isoform switching of CLSTN1 is crucial for EMT. Moreover, AKAP8 expression and the alternative splicing of CLSTN1 predict breast cancer patient survival. Together, our work demonstrates the essentiality of RNA metabolism that impinges on metastatic breast cancer.