A general calculus of fitness landscapes finds genes under selection in cancers

Genetic variants drive the evolution of traits and diseases. We previously modeled these variants as small displacements in fitness landscapes and estimated their functional impact by differentiating the evolutionary relationship between genotype and phenotype. Conversely, here we integrate these derivatives to identify genes steering specific traits. Over cancer cohorts, integration identified 460 likely tumor-driving genes. Many have literature and experimental support but had eluded prior genomic searches for positive selection in tumors. Beyond providing cancer insights, these results introduce a general calculus of evolution to quantify the genotype-phenotype relationship and discover genes associated with complex traits and diseases.

EPIMUTESTR: a nearest neighbor machine learning approach to predict cancer driver genes from the evolutionary action of coding variants

Discovering rare cancer driver genes is difficult because their mutational frequency is too low for statistical detection by computational methods. EPIMUTESTR is an integrative nearest-neighbor machine learning algorithm that identifies such marginal genes by modeling the fitness of their mutations with the phylogenetic Evolutionary Action (EA) score. Over cohorts of sequenced patients from The Cancer Genome Atlas representing 33 tumor types, EPIMUTESTR detected 214 previously inferred cancer driver genes and 137 new candidates never identified computationally before of which seven genes are supported in the COSMIC Cancer Gene Census. EPIMUTESTR achieved better robustness and specificity than existing methods in a number of benchmark methods and datasets.

K-Ras and p53 mouse model with molecular characteristics of human rhabdomyosarcoma and translational applications

Rhabdomyosarcoma (RMS) is the most common soft tissue sarcoma in children, with overall long-term survival rates of ∼65-70%. Thus, additional molecular insights and representative models are critical for identifying and evaluating new treatment modalities. Using MyoD-Cre-mediated introduction of mutant K-RasG12D and perturbations in p53, we developed a novel genetically engineered mouse model (GEMM) for RMS. The anatomic sites of primary RMS development recapitulated human disease, including tumors in the head, neck, extremities and abdomen. We confirmed RMS histology and diagnosis through Hematoxylin and Eosin staining, and positive immunohistochemical staining for desmin, myogenin, and phosphotungstic acid-Hematoxylin. Cell lines from GEMM tumors were established with the ability to engraft in immunocompetent mice with comparable histological and staining features as the primary tumors. Tail vein injection of cell lines had high metastatic potential to the lungs. Transcriptomic analyses of p53R172H/K-RasG12D GEMM-derived tumors showed evidence of high molecular homology with human RMS. Finally, pre-clinical use of these murine RMS lines showed similar therapeutic responsiveness to chemotherapy and targeted therapies as human RMS cell lines.

A pediatric brain tumor atlas of genes deregulated by somatic genomic rearrangement

The global impact of somatic structural variants (SSVs) on gene expression in pediatric brain tumors has not been thoroughly characterised. Here, using whole-genome and RNA sequencing from 854 tumors of more than 30 different types from the Children's Brain Tumor Tissue Consortium, we report the altered expression of hundreds of genes in association with the presence of nearby SSV breakpoints. SSV-mediated expression changes involve gene fusions, altered cis-regulation, or gene disruption. SSVs considerably extend the numbers of patients with tumors somatically altered for critical pathways, including receptor tyrosine kinases (KRAS, MET, EGFR, NF1), Rb pathway (CDK4), TERT, MYC family (MYC, MYCN, MYB), and HIPPO (NF2). Compared to initial tumors, progressive or recurrent tumors involve a distinct set of SSV-gene associations. High overall SSV burden associates with TP53 mutations, histone H3.3 gene H3F3C mutations, and the transcription of DNA damage response genes. Compared to adult cancers, pediatric brain tumors would involve a different set of genes with SSV-altered cis-regulation. Our comprehensive and pan-histology genomic analyses reveal SSVs to play a major role in shaping the transcriptome of pediatric brain tumors.

Integrated Analysis of TP53 Gene and Pathway Alterations in The Cancer Genome Atlas

The TP53 tumor suppressor gene is frequently mutated in human cancers. An analysis of five data platforms in 10,225 patient samples from 32 cancers reported by The Cancer Genome Atlas (TCGA) enables comprehensive assessment of p53 pathway involvement in these cancers. More than 91% of TP53-mutant cancers exhibit second allele loss by mutation, chromosomal deletion, or copy-neutral loss of heterozygosity. TP53 mutations are associated with enhanced chromosomal instability, including increased amplification of oncogenes and deep deletion of tumor suppressor genes. Tumors with TP53 mutations differ from their non-mutated counterparts in RNA, miRNA, and protein expression patterns, with mutant TP53 tumors displaying enhanced expression of cell cycle progression genes and proteins. A mutant TP53 RNA expression signature shows significant correlation with reduced survival in 11 cancer types. Thus, TP53 mutation has profound effects on tumor cell genomic structure, expression, and clinical outlook.

WIP1 dephosphorylation of p27Kip1 Serine 140 destabilizes p27Kip1 and reverses anti-proliferative effects of ATM phosphorylation

The phosphoinositide-3-kinase like kinases (PIKK) such as ATM and ATR play a key role in initiating the cellular DNA damage response (DDR). One key ATM target is the cyclin-dependent kinase inhibitor p27^Kip1 that promotes G1 arrest. ATM activates p27^Kip1-induced arrest in part through phosphorylation of p27^Kip1 at Serine 140. Here we show that this site is dephosphorylated by the type 2C serine/threonine phosphatase, WIP1 (Wildtype p53-Induced Phosphatase-1), encoded by the PPM1D gene. WIP1 has been shown to dephosphorylate numerous ATM target sites in DDR proteins, and its overexpression and/or mutation has often been associated with oncogenesis. We demonstrate that wildtype, but not phosphatase-dead WIP1, efficiently dephosphorylates p27^Kip1 Ser140 both in vitro and in cells and that this dephosphorylation is sensitive to the WIP1-specific inhibitor GSK 2830371. Increased expression of wildtype WIP1 reduces stability of p27^Kip1 while increased expression of similar amounts of phosphatase-dead WIP1 has no effect on p27^Kip1 protein stability. Overexpression of wildtype p27^Kip1 reduces cell proliferation and colony forming capability relative to the S140A (constitutively non-phosphorylated) form of p27. Thus, WIP1 plays a significant role in homeostatic modulation of p27^Kip1 activity following activation by ATM.

PPM1D Mutations Drive Clonal Hematopoiesis in Response to Cytotoxic Chemotherapy

Clonal hematopoiesis (CH), in which stem cell clones dominate blood production, becomes increasingly common with age and can presage malignancy development. The conditions that promote ascendancy of particular clones are unclear. We found that mutations in PPM1D (protein phosphatase Mn2+/Mg2+-dependent 1D), a DNA damage response regulator that is frequently mutated in CH, were present in one-fifth of patients with therapy-related acute myeloid leukemia or myelodysplastic syndrome and strongly correlated with cisplatin exposure. Cell lines with hyperactive PPM1D mutations expand to outcompete normal cells after exposure to cytotoxic DNA damaging agents including cisplatin, and this effect was predominantly mediated by increased resistance to apoptosis. Moreover, heterozygous mutant Ppm1d hematopoietic cells outcompeted their wild-type counterparts in vivo after exposure to cisplatin and doxorubicin, but not during recovery from bone marrow transplantation. These findings establish the clinical relevance of PPM1D mutations in CH and the importance of studying mutation-treatment interactions. VIDEO ABSTRACT.

Transglutaminase-2 promotes metastatic and stem-like phenotypes in osteosarcoma

Osteosarcoma (OS) is a highly aggressive mesenchymal malignancy and the most common primary bone tumor in the pediatric population. OS frequently presents with or develops distal metastases. Patients with metastatic disease have extremely poor survival rates, thus necessitating improved molecular insights into OS metastatic biology. Utilizing our previously characterized genetically engineered mouse model (GEMM) of metastatic OS, we identified enhanced differential expression of Transglutaminase-2 (TGM2) in metastatic OS. However, the role of TGM2 in sarcoma development and metastatic progression remains largely undefined. To further investigate the role of TGM2 in OS metastasis, we performed both gain- and loss-of-function studies for TGM2 in human and mouse OS cell lines. Our data provide evidence that enhanced expression of TGM2 in metastatic OS contributes to migratory and invasive phenotypes. Besides the effects on metastatic phenotypes, we also observed that TGM2 contributes to OS stem-like properties. In addition, treatment with transglutaminase inhibitors had analogous effects on proliferation and migration to TGM2 knockdown. Finally, in vivo xenograft studies demonstrated that TGM2 functionally alters metastatic potential and survival outcome. Together, these data highlight TGM2 as a pro-metastatic factor in OS and a potential avenue for future therapeutic intervention to inhibit metastatic disease.

Genomic and Molecular Landscape of DNA Damage Repair Deficiency across The Cancer Genome Atlas

DNA damage repair (DDR) pathways modulate cancer risk, progression, and therapeutic response. We systematically analyzed somatic alterations to provide a comprehensive view of DDR deficiency across 33 cancer types. Mutations with accompanying loss of heterozygosity were observed in over 1/3 of DDR genes, including TP53 and BRCA1/2. Other prevalent alterations included epigenetic silencing of the direct repair genes EXO5, MGMT, and ALKBH3 in ∼20% of samples. Homologous recombination deficiency (HRD) was present at varying frequency in many cancer types, most notably ovarian cancer. However, in contrast to ovarian cancer, HRD was associated with worse outcomes in several other cancers. Protein structure-based analyses allowed us to predict functional consequences of rare, recurrent DDR mutations. A new machine-learning-based classifier developed from gene expression data allowed us to identify alterations that phenocopy deleterious TP53 mutations. These frequent DDR gene alterations in many human cancers have functional consequences that may determine cancer progression and guide therapy.

Enhanced inflammation and attenuated tumor suppressor pathways are associated with oncogene-induced lung tumors in aged mice

Aging is often accompanied by a dramatic increase in cancer susceptibility. To gain insights into how aging affects tumor susceptibility, we generated a conditional mouse model in which oncogenic Kras^G12D was activated specifically in lungs of young (3–5 months) and old (19–24 months) mice. Activation of Kras^G12D in old mice resulted in shorter survival and development of higher‐grade lung tumors. Six weeks after Kras^G12D activation, old lung tissues contained higher numbers of adenomas than their young tissue counterparts. Lung tumors in old mice displayed higher proliferation rates, as well as attenuated DNA damage and p53 tumor suppressor responses. Gene expression comparison of lung tumors from young and old mice revealed upregulation of extracellular matrix‐related genes in young tumors, indicative of a robust cancer‐associated fibroblast response. In old tumors, numerous inflammation‐related genes such as Ccl7,IL‐1β, Cxcr6, and IL‐15ra were consistently upregulated. Increased numbers of immune cells were localized around the periphery of lung adenomas from old mice. Our experiments indicate that more aggressive lung tumor formation in older Kras^G12D mice may be in part the result of subdued tumor suppressor and DNA damage responses, an enhanced inflammatory milieu, and a more accommodating tissue microenvironment.

Recommended Guidelines for Validation, Quality Control, and Reporting of TP53 Variants in Clinical Practice

Accurate assessment of TP53 gene status in sporadic tumors and in the germline of individuals at high risk of cancer due to Li-Fraumeni Syndrome (LFS) has important clinical implications for diagnosis, surveillance, and therapy. Genomic data from more than 20,000 cancer genomes provide a wealth of information on cancer gene alterations and have confirmed TP53 as the most commonly mutated gene in human cancer. Analysis of a database of 70,000 TP53 variants reveals that the two newly discovered exons of the gene, exons 9β and 9γ, generated by alternative splicing, are the targets of inactivating mutation events in breast, liver, and head and neck tumors. Furthermore, germline rearrange-ments in intron 1 of TP53 are associated with LFS and are frequently observed in sporadic osteosarcoma. In this context of constantly growing genomic data, we discuss how screening strategies must be improved when assessing TP53 status in clinical samples. Finally, we discuss how TP53 alterations should be described by using accurate nomenclature to avoid confusion in scientific and clinical reports. Cancer Res; 77(6); 1250-60. ©2017 AACR.

miR-130b directly targets ARHGAP1 to drive activation of a metastatic CDC42-PAK1-AP1 positive feedback loop in Ewing sarcoma

Ewing Sarcoma (ES) is a highly aggressive bone tumor with peak incidence in the adolescent population. It has a high propensity to metastasize, which is associated with dismal survival rates of approximately 25%. To further understand mechanisms of metastasis we investigated microRNA regulatory networks in ES. Our studies focused on miR-130b due to our analysis that enhanced expression of this microRNA has clinical relevance in multiple sarcomas, including ES. Our studies provide insights into a novel positive feedback network involving the direct regulation of miR-130b and activation of downstream signaling events contributing toward sarcoma metastasis. Specifically, we demonstrated miR-130b induces proliferation, invasion, and migration in vitro and increased metastatic potential in vivo. Using microarray analysis of ES cells with differential miR-130b expression we identified alterations in downstream signaling cascades including activation of the CDC42 pathway. We identified ARHGAP1, which is a negative regulator of CDC42, as a novel, direct target of miR-130b. In turn, downstream activation of PAK1 activated the JNK and AP-1 cascades and downstream transcriptional targets including IL-8, MMP1 and CCND1. Furthermore, chromatin immunoprecipitation of endogenous AP-1 in ES cells demonstrated direct binding to an upstream consensus binding site within the miR-130b promoter. Finally, small molecule inhibition of PAK1 blocked miR-130b activation of JNK and downstream AP-1 target genes, including primary miR-130b transcripts, and miR-130b oncogenic properties, thus identifying PAK1 as a novel therapeutic target for ES. Taken together, our findings identify and characterize a novel, targetable miR-130b regulatory network that promotes ES metastasis.

LZAP is a novel Wip1 binding partner and positive regulator of its phosphatase activity in vitro

The phosphatase Wip1 attenuates the DNA damage response (DDR) by removing phosphorylation marks from a number of DDR proteins (p53, MDM2, Chk1/2, p38). Wip1 also dephosphorylates and inactivates RelA. Notably, LZAP, a putative tumor suppressor, has been linked to dephosphorylation of several of these substrates, including RelA, p38, Chk1, and Chk2. LZAP has no known catalytic activity or functional motifs, suggesting that it exerts its effects through interaction with other proteins. Here we show that LZAP binds Wip1 and stimulates its phosphatase activity. LZAP had been previously shown to bind many Wip1 substrates (RelA, p38, Chk1/2), and our results show that LZAP also binds the previously identified Wip1 substrate, MDM2. This work identifies 2 novel Wip1 substrates, ERK1 and HuR, and demonstrates that HuR is a binding partner of LZAP. Pleasingly, LZAP potentiated Wip1 catalytic activity toward each substrate tested, regardless of whether full-length substrates or phosphopeptides were utilized. Since this effect was observed on ERK1, which does not bind LZAP, as well as for each of 7 peptides tested, we hypothesize that LZAP binding to the substrate is not required for this effect and that LZAP directly binds Wip1 to augment its phosphatase activity.

Secreted Frizzled-Related Protein 2 (sFRP2) promotes osteosarcoma invasion and metastatic potential

Background

Osteosarcoma (OS), which has a high potential for developing metastatic disease, is the most frequent malignant bone tumor in children and adolescents. Molecular analysis of a metastatic genetically engineered mouse model of osteosarcoma identified enhanced expression of Secreted Frizzled-Related Protein 2 (sFRP2), a putative regulator of Wnt signaling within metastatic tumors. Subsequent analysis correlated increased expression in the human disease, and within highly metastatic OS cells. However, the role of sFRP2 in osteosarcoma development and progression has not been well elucidated.

Methods

Studies using stable gain or loss-of-function alterations of sFRP2 within human and mouse OS cells were performed to assess changes in cell proliferation, migration, and invasive ability in vitro, via both transwell and 3D matrigel assays. In additional, xenograft studies using overexpression of sFRP2 were used to assess effects on in vivo metastatic potential.

Results

Functional studies revealed stable overexpression of sFRP2 within localized human and mouse OS cells significantly increased cell migration and invasive ability in vitro and enhanced metastatic potential in vivo. Additional studies exploiting knockdown of sFRP2 within metastatic human and mouse OS cells demonstrated decreased cell migration and invasion ability in vitro, thus corroborating a critical biological phenotype carried out by sFRP2. Interestingly, alterations in sFRP2 expression did not alter OS proliferation rates or primary tumor development.

Conclusions

While future studies further investigating the molecular mechanisms contributing towards this sFRP2-dependent phenotype are needed, our studies clearly provide evidence that aberrant expression of sFRP2 can contribute to the invasive and metastatic potential for osteosarcoma.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-016-2909-6) contains supplementary material, which is available to authorized users.

Ampullary Cancers Harbor ELF3 Tumor Suppressor Gene Mutations and Exhibit Frequent WNT Dysregulation

The ampulla of Vater is a complex cellular environment from which adenocarcinomas arise to form a group of histopathologically heterogenous tumors. To evaluate the molecular features of these tumors, 98 ampullary adenocarcinomas were evaluated and compared to 44 distal bile duct and 18 duodenal adenocarcinomas. Genomic analyses revealed mutations in the WNT signaling pathway among half of the patients and in all three adenocarcinomas irrespective of their origin and histological morphology. These tumors were characterized by a high frequency of inactivating mutations of ELF3, a high rate of microsatellite instability, and common focal deletions and amplifications, suggesting common attributes in the molecular pathogenesis are at play in these tumors. The high frequency of WNT pathway activating mutation, coupled with small-molecule inhibitors of β-catenin in clinical trials, suggests future treatment decisions for these patients may be guided by genomic analysis.

Genomic profiling of Sézary syndrome identifies alterations of key T cell signaling and differentiation genes

Sézary syndrome is a rare leukemic form of cutaneous T cell lymphoma characterized by generalized redness, scaling, itching and increased numbers of circulating atypical T lymphocytes. It is rarely curable, with poor prognosis. Here we present a multiplatform genomic analysis of 37 patients with Sézary syndrome that implicates dysregulation of cell cycle checkpoint and T cell signaling. Frequent somatic alterations were identified in TP53, CARD11, CCR4, PLCG1, CDKN2A, ARID1A, RPS6KA1 and ZEB1. Activating CCR4 and CARD11 mutations were detected in nearly one-third of patients. ZEB1, encoding a transcription repressor essential for T cell differentiation, was deleted in over one-half of patients. IL32 and IL2RG were overexpressed in nearly all cases. Our results demonstrate profound disruption of key signaling pathways in Sézary syndrome and suggest potential targets for new therapies.

Ghrelin Prevents Cisplatin-Induced Testicular Damage by Facilitating Repair of DNA Double Strand Breaks Through Activation of p53 in Mice

Cisplatin administration induces DNA damage resulting in germ cell apoptosis and subsequent testicular atrophy. Although 50 percent of male cancer patients receiving cisplatin-based chemotherapy develop long-term secondary infertility, medical treatment to prevent spermatogenic failure after chemotherapy is not available. Under normal conditions, testicular p53 promotes cell cycle arrest, which allows time for DNA repair and reshuffling during meiosis. However, its role in the setting of cisplatin-induced infertility has not been studied. Ghrelin administration ameliorates the spermatogenic failure that follows cisplatin administration in mice, but the mechanisms mediating these effects have not been well established. The aim of the current study was to characterize the mechanisms of ghrelin and p53 action in the testis after cisplatin-induced testicular damage. Here we show that cisplatin induces germ cell damage through inhibition of p53-dependent DNA repair mechanisms involving gamma-H2AX and ataxia telangiectasia mutated protein kinase. As a result, testicular weight and sperm count and motility were decreased with an associated increase in sperm DNA damage. Ghrelin administration prevented these sequelae by restoring the normal expression of gamma-H2AX, ataxia telangiectasia mutated, and p53, which in turn allows repair of DNA double stranded breaks. In conclusion, these findings indicate that ghrelin has the potential to prevent or diminish infertility caused by cisplatin and other chemotherapeutic agents by restoring p53-dependent DNA repair mechanisms.

Cross-species identification of a plasma microRNA signature for detection, therapeutic monitoring, and prognosis in osteosarcoma

Osteosarcoma (OS) is the primary bone tumor in children and young adults. Currently, there are no reliable, noninvasive biologic markers to detect the presence or progression of disease, assess therapy response, or provide upfront prognostic insights. MicroRNAs (miRNAs) are evolutionarily conserved, stable, small noncoding RNA molecules that are key posttranscriptional regulators and are ideal candidates for circulating biomarker development due to their stability in plasma, ease of isolation, and the unique expressions associated with specific disease states. Using a qPCR-based platform that analyzes more than 750 miRNAs, we analyzed control and diseased-associated plasma from a genetically engineered mouse model of OS to identify a profile of four plasma miRNAs. Subsequent analysis of 40 human patient samples corroborated these results. We also identified disease-specific endogenous reference plasma miRNAs for mouse and human studies. Specifically, we observed plasma miR-205-5p was decreased 2.68-fold in mice with OS compared to control mice, whereas, miR-214, and miR-335-5p were increased 2.37- and 2.69-fold, respectively. In human samples, the same profile was seen with miR-205-5p decreased 1.75-fold in patients with OS, whereas miR-574-3p, miR-214, and miR-335-5p were increased 3.16-, 8.31- and 2.52-fold, respectively, compared to healthy controls. Furthermore, low plasma levels of miR-214 in metastatic patients at time of diagnosis conveyed a significantly better overall survival. This is the first study to identify plasma miRNAs that could be used to prospectively identify disease, potentially monitor therapeutic efficacy and have prognostic implications for OS patients.

WIP1 phosphatase as a potential therapeutic target in neuroblastoma

The wild-type p53-induced phosphatase 1 (WIP1) is a serine/threonine phosphatase that negatively regulates multiple proteins involved in DNA damage response including p53, CHK2, Histone H2AX, and ATM, and it has been shown to be overexpressed or amplified in human cancers including breast and ovarian cancers. We examined WIP1 mRNA levels across multiple tumor types and found the highest levels in breast cancer, leukemia, medulloblastoma and neuroblastoma. Neuroblastoma is an exclusively TP53 wild type tumor at diagnosis and inhibition of p53 is required for tumorigenesis. Neuroblastomas in particular have previously been shown to have 17q amplification, harboring the WIP1 (PPM1D) gene and associated with poor clinical outcome. We therefore sought to determine whether inhibiting WIP1 with a selective antagonist, GSK2830371, can attenuate neuroblastoma cell growth through reactivation of p53 mediated tumor suppression. Neuroblastoma cell lines with wild-type TP53 alleles were highly sensitive to GSK2830371 treatment, while cell lines with mutant TP53 were resistant to GSK2830371. The majority of tested neuroblastoma cell lines with copy number gains of the PPM1D locus were also TP53 wild-type and sensitive to GSK2830371A; in contrast cell lines with no copy gain of PPM1D were mixed in their sensitivity to WIP1 inhibition, with the primary determinant being TP53 mutational status. Since WIP1 is involved in the cellular response to DNA damage and drugs used in neuroblastoma treatment induce apoptosis through DNA damage, we sought to determine whether GSK2830371 could act synergistically with standard of care chemotherapeutics. Treatment of wild-type TP53 neuroblastoma cell lines with both GSK2830371 and either doxorubicin or carboplatin resulted in enhanced cell death, mediated through caspase 3/7 induction, as compared to either agent alone. Our data suggests that WIP1 inhibition represents a novel therapeutic approach to neuroblastoma that could be integrated with current chemotherapeutic approaches.

Dynamic analyses of alternative polyadenylation from RNA-seq reveal a 3'-UTR landscape across seven tumour types

Alternative polyadenylation (APA) is a pervasive mechanism in the regulation of most human genes, and its implication in diseases including cancer is only beginning to be appreciated. Since conventional APA profiling has not been widely adopted, global cancer APA studies are very limited. Here we develop a novel bioinformatics algorithm (DaPars) for the de novo identification of dynamic APAs from standard RNA-seq. When applied to 358 TCGA Pan-Cancer tumour/normal pairs across seven tumour types, DaPars reveals 1,346 genes with recurrent and tumour-specific APAs. Most APA genes (91%) have shorter 3'-untranslated regions (3' UTRs) in tumours that can avoid microRNA-mediated repression, including glutaminase (GLS), a key metabolic enzyme for tumour proliferation. Interestingly, selected APA events add strong prognostic power beyond common clinical and molecular variables, suggesting their potential as novel prognostic biomarkers. Finally, our results implicate CstF64, an essential polyadenylation factor, as a master regulator of 3'-UTR shortening across multiple tumour types.

Homeostatic Regulation of Base Excision Repair by a p53-Induced Phosphatase: Linking Stress Response Pathways with DNA Repair Proteins

The p53 protein plays a central role in the integration of cellular stress responses. If the cell incurs DNA damage, p53-induced cell cycle arrest is accompanied by p53-facilitated DNA repair. In particular, p53 has been demonstrated to promote both nucleotide excision repair (NER) and base excision repair (BER). Once these repair processes are completed, p53 activity declines and the cell can reenter the cell cycle. A critical mediator of this p53 negative regulatory feedback loop is Mdm2, a p53 transcriptional target whose protein mediates p53 proteolytic degradation. Another such p53 transcriptional target that may function in a p53 negative regulation is the PPM1D phosphatase. PPM1D may inhibit p53 activity through inactivating dephosphorylation of the p38 MAP kinase. We have recently shown that PPM1D suppresses BER in part through dephosphorylation of a key BER effector, the nuclear isoform of uracil DNA glycosylase, or UNG2. This finding further links p53 signaling to DNA repair pathways and illustrates a mechanism by which activated DNA repair systems are returned to a deactivated, homeostatic state.

Abstract 2211: Viral subtypes in hepatocellular carcinoma are associated with different mechanisms of WNT/CTNNB1 alteration

Hepatocellular carcinoma (HCC) is one of the most prevalent human cancers, especially in Asia, and regions of endemic hepatitis infection. The primary risk factors for HCC are chronic hepatitis infection by hepatitis C (HCV) or hepatitis B (HBV). We performed whole exome sequencing on 549 matched tumor and normal samples from three centers (BCM HGSC, RCAST, NCC) to identify somatic changes in protein coding regions and changes in copy number. Here we describe variants associated with specific viral etiologies in HCC.

We found TP53, CTNNB1, ARID1A, ARID2, and AXIN1 to be among the most frequently mutated genes in the entire dataset, implicating DNA damage, WNT signaling, and chromatin remodeling as recurrently altered pathways in HCC. Our validation rate for mutations in significantly mutated genes was over 95%. We identified 18 genes which differed significantly (p < 0.05) in mutation frequency between viral subtypes. Of these, AXIN1 was mutated in over 15% of HBV-associated HCC but less than 5% of HCV and non-viral-associated HCC. Additionally, we performed GISTIC analysis based on read depth from our whole exome sequencing data. We found that WNT3A and MYC were significantly amplified and AXIN1 and TP53 were significantly deleted in HCV-associated HCC but not HBV or non-virus-associated HCC.

Our findings indicate that WNT/beta-catenin signaling is a frequently altered pathway in HCC. While aberrations in this pathway were present in all viral etiologies, HBV-associated HCC tended to mutate AXIN1 while HCV-associated HCC tended to have AXIN1 deletion. These differences may result from virus-specific oncogenic processes.

Citation Format: Kyle R. Covington, Lawrence A. Donehower, Chad Creighton, Betty L. Slagle, John A. Goss, Ronald T. Cotton, Marie-Claude Gingras, Eve Shinbrot, Jacfranz J. Guiteau, Thao N. Nguyen, Theresa R. Harring, Donna M. Muzny, Kimberly Walker, HarshaVardhan Doddapaneni, Richard A. Gibbs, Hiroyuki Aburatani, Tatsuhiro Shibata, David A. Wheeler, Japan ICGC HCC Project. Viral subtypes in hepatocellular carcinoma are associated with different mechanisms of WNT/CTNNB1 alteration. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 2211. doi:10.1158/1538-7445.AM2014-2211

Notch activation as a driver of osteogenic sarcoma

Osteogenic sarcoma (OS) is a deadly skeletal malignancy whose cause is unknown. We report here a mouse model of OS based on conditional expression of the intracellular domain of Notch1 (NICD). Expression of the NICD in immature osteoblasts was sufficient to drive the formation of bone tumors, including OS, with complete penetrance. These tumors display features of human OS; namely, histopathology, cytogenetic complexity, and metastatic potential. We show that Notch activation combined with loss of p53 synergistically accelerates OS development in mice, although p53-driven OS is not Rbpj dependent, which demonstrates a dual dominance of the Notch oncogene and p53 mutation in the development of OS. Using this model, we also reveal the osteoblasts as the potential sources of OS.

The somatic genomic landscape of chromophobe renal cell carcinoma

We describe the landscape of somatic genomic alterations of 66 chromophobe renal cell carcinomas (ChRCCs) on the basis of multidimensional and comprehensive characterization, including mtDNA and whole-genome sequencing. The result is consistent that ChRCC originates from the distal nephron compared with other kidney cancers with more proximal origins. Combined mtDNA and gene expression analysis implicates changes in mitochondrial function as a component of the disease biology, while suggesting alternative roles for mtDNA mutations in cancers relying on oxidative phosphorylation. Genomic rearrangements lead to recurrent structural breakpoints within TERT promoter region, which correlates with highly elevated TERT expression and manifestation of kataegis, representing a mechanism of TERT upregulation in cancer distinct from previously observed amplifications and point mutations.

Effects of TP53 mutational status on gene expression patterns across 10 human cancer types

Mutations in the TP53 tumour suppressor gene occur in half of all human cancers, indicating its critical importance in inhibiting cancer development. Despite extensive studies, the mechanisms by which mutant p53 enhances tumour progression remain only partially understood. Here, using data from the Cancer Genome Atlas (TCGA), genomic and transcriptomic analyses were performed on 2256 tumours from 10 human cancer types. We show that tumours with TP53 mutations have altered gene expression profiles compared to tumours retaining two wild-type TP53 alleles. Among 113 known p53-up-regulated target genes identified from cell culture assays, 10 were consistently up-regulated in at least eight of 10 cancer types that retain both copies of wild-type TP53. RPS27L, CDKN1A (p21(CIP1)) and ZMAT3 were significantly up-regulated in all 10 cancer types retaining wild-type TP53. Using this p53-based expression analysis as a discovery tool, we used cell-based assays to identify five novel p53 target genes from genes consistently up-regulated in wild-type p53 cancers. Global gene expression analyses revealed that cell cycle regulatory genes and transcription factors E2F1, MYBL2 and FOXM1 were disproportionately up-regulated in many TP53 mutant cancer types. Finally, > 93% of tumours with a TP53 mutation exhibited greatly reduced wild-type p53 messenger expression, due to loss of heterozygosity or copy neutral loss of heterozygosity, supporting the concept of p53 as a recessive tumour suppressor. The data indicate that tumours with wild-type TP53 retain some aspects of p53-mediated growth inhibitory signalling through activation of p53 target genes and suppression of cell cycle regulatory genes.

Combinatorial therapy discovery using mixed integer linear programming

MOTIVATION

Combinatorial therapies play increasingly important roles in combating complex diseases. Owing to the huge cost associated with experimental methods in identifying optimal drug combinations, computational approaches can provide a guide to limit the search space and reduce cost. However, few computational approaches have been developed for this purpose, and thus there is a great need of new algorithms for drug combination prediction.

RESULTS

Here we proposed to formulate the optimal combinatorial therapy problem into two complementary mathematical algorithms, Balanced Target Set Cover (BTSC) and Minimum Off-Target Set Cover (MOTSC). Given a disease gene set, BTSC seeks a balanced solution that maximizes the coverage on the disease genes and minimizes the off-target hits at the same time. MOTSC seeks a full coverage on the disease gene set while minimizing the off-target set. Through simulation, both BTSC and MOTSC demonstrated a much faster running time over exhaustive search with the same accuracy. When applied to real disease gene sets, our algorithms not only identified known drug combinations, but also predicted novel drug combinations that are worth further testing. In addition, we developed a web-based tool to allow users to iteratively search for optimal drug combinations given a user-defined gene set.

AVAILABILITY

Our tool is freely available for noncommercial use at http://www.drug.liuzlab.org/.

CONTACT

zhandong.liu@bcm.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Abstract 3867: Novel mouse models to investigate the molecular pathogenesis of metastatic osteosarcoma

Background: Osteosarcoma (OS) is the primary bone tumor in the pediatric population with the main determinant for patient prognosis being the presence of metastatic disease. Approximately 25-30% of patients will present with metastatic osteosarcoma. Patients with metastatic disease have long-term survival rates often <20%.

Objectives: We hypothesize that the molecular pathogenesis of metastatic osteosarcoma is different from localized disease. We intend this work to provide a relevant, endogenous model of metastatic OS that can be utilized to advance our understanding of the molecular pathogenesis of the disease, insights into novel therapeutic targets and as a pre-clinical model for investigating the efficacy of novel therapies.

Design: We have developed a tissue-specific alteration of the p53 status by using osteoblast specific Cre-recombinase expressing mice to generate progeny that spontaneously form endogenous osteosarcomas. Through the use of a mutated, gain of function form of p53, shown previously to be associated with metastatic disease, we have developed a novel immunocompetent model that significantly enhances the endogenous development of metastatic OS.

Results: Tumor analysis has revealed genetic insights in the metastatic progression of osteosarcoma. These include the significant downregulation of Wnt-signaling inhibitors, such as NKD2, APCDD1 and Wnt5a in the metastatic tumors. Functional studies have determined that overexpression of NKD2, also downregulated in several human OS metastatic cell lines, in metastatic mouse OS cell lines leads to a significant decrease in metastatic lung lesions upon transplantation into immunodeficient mice. Possible NKD2-dependent functions include regulation of not only the Wnt signaling pathway, but also blood vessel formation, regulation of cell migration and cell adhesion.

We have also noted the dysregulation of several critical microRNAs in metastatic OS, including the upregulation of mir-130b in the metastatic lesions. This particular microRNA has recently been shown to have clinical correlation in Ewing's sarcoma, with higher levels of tumor mir-130b having a significantly poorer outcome. We are actively pursuing the role of this (and other) microRNA(s) in the molecular pathogenesis of metastatic osteosarcoma.

Conclusion: This novel model has enabled valuable molecular insights into the development and progression of metastatic OS that can lead to the identification of novel therapeutic targets.

Citation Format: Shuying Zhao, Lyazat Kurenbekova, Lawrence A. Donehower, Jason T. Yustein. Novel mouse models to investigate the molecular pathogenesis of metastatic osteosarcoma. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 3867. doi:10.1158/1538-7445.AM2013-3867

Abstract 1114: Polymerase epsilon (POLE) mutations and mutator phenotypes in colorectal and endometrial tumors

Colorectal rectal and endometrial cancers are divided into microsatellite instable (MSI) and microsatellite stable (MSS) types. MSS tumors have chromosome instability with mutation rates of 1-10/Mb. MSI patients have a better prognosis than MSS patients and are treated less aggressively. These tumors exhibit microsatellite instability (MSI), and the CpG island methylator phenotype (CIMP), an increased mutation rate (10-100/Mb, hypermutated). Here we demonstrate a novel class of tumor in these two diseases exhibiting an ultramutator phenotype with mutation rates exceeding 100/ Mb, MSS and chromosome stable. The ultra-high mutation rates appear to be caused by recurrent mutations in the exonuclease domain of DNA polymerase epsilon (POLE).

Matched tumor and normal whole exomes for 512 colorectal and 248 endometrial cancers (TCGA data control center) were evaluated for mutation frequency and MSI status. Among the colorectal tumors, 70/512(14%) were hyper mutated, endometroid tumors had 65/248(26%) hypermutated tumors. Within the hypermutated tumors, we identified tumors with extremely high mutation rates (100 mutations/Mb) termed “ultramutated”. 14 (3%) colorectal and 17 (7%) endometrial were ultramutated. 100% of the colorectal and endometrial ultramutated contain an exonuclease domain mutation in POLE with three recurrently mutated positions: P286R/H/S and V411L and S459F, these mutations are found exclusively in ultramutated samples. The functional relevance of these mutations has been demonstrated by previous mutational analysis of the exonuclease domain in bacteria phage (T4) yeast,and mice. Disruption of residues in the active site leads to high error rates, a mutator phenotype and tumor formation. The defined exonuclease active site aligns closely with the recurrent mutation found in the ultramutated tumors.

The ultramuted tumors have a distinct mutation spectrum, occur in context and exhibit strand biases. POLE exonuclease domain mutation tumors have high rates of C to A and T to G mutations with low rates of C to G and T to A. Changes of C to A occur predominately in context of 3’ and 5’ flanking T bases. This mutation spectrum and context of flanking bases suggests replicative strand bias. We found that mutations on the leading strand, exhibited a preference for C to A over G to T (60:40). Suggesting transcription-coupled repair (TCR) plays a role in ultramutation.

Our results demonstrate that ultra mutated tumors are driven by recurrent mutations in the exonuclease domain of DNA polymerase epsilon. The distinct mutation spectrum, strand bias and mutation context all suggest TCR is involved in ultramutated tumors. We describe a previously unknown role of POLE in colorectal tumors as well as identification of a previously undescribed hot spot for cancer mutations.

Citation Format: Eve Shinbrot, Nils Weinhold, Nikolaus Schultz, Lawrence A. Donehower, Jennifer Drummond, Kyle Chang, Richard Gibbs, Chris Sander, David A. Wheeler. Polymerase epsilon (POLE) mutations and mutator phenotypes in colorectal and endometrial tumors. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 1114. doi:10.1158/1538-7445.AM2013-1114

MLH1-silenced and non-silenced subgroups of hypermutated colorectal carcinomas have distinct mutational landscapes

Approximately 15% of colorectal carcinomas (CRCs) exhibit a hypermutated genotype accompanied by high levels of microsatellite instability (MSI-H) and defects in DNA mismatch repair. These tumours, unlike the majority of colorectal carcinomas, are often diploid, exhibit frequent epigenetic silencing of the MLH1 DNA mismatch repair gene, and have a better clinical prognosis. As an adjunct study to The Cancer Genome Atlas consortium that recently analysed 224 colorectal cancers by whole exome sequencing, we compared the 35 CRCs (15.6%) with a hypermutated genotype to those with a non-hypermutated genotype. We found that 22 (63%) of the hypermutated CRCs exhibited transcriptional silencing of the MLH1 gene, a high frequency of BRAF V600E gene mutations, and infrequent APC and KRAS mutations, a mutational pattern significantly different from their non-hypermutated counterparts. However, the remaining 13 (37%) hypermutated CRCs lacked MLH1 silencing, contained tumours with the highest mutation rates ('ultramutated' CRCs), and exhibited higher incidences of APC and KRAS mutations, but infrequent BRAF mutations. These patterns were confirmed in an independent validation set of 250 exome-sequenced CRCs. Analysis of mRNA and microRNA expression signatures revealed that hypermutated CRCs with MLH1 silencing had greatly reduced levels of WNT signalling and increased BRAF signalling relative to non-hypermutated CRCs. Our findings suggest that hypermutated CRCs include one subgroup with fundamentally different pathways to malignancy than the majority of CRCs. Examination of MLH1 expression status and frequencies of APC, KRAS, and BRAF mutation in CRC may provide a useful diagnostic tool that could supplement the standard microsatellite instability assays and influence therapeutic decisions.

Activation of WIP1 phosphatase by HTLV-1 Tax mitigates the cellular response to DNA damage

Genomic instability stemming from dysregulation of cell cycle checkpoints and DNA damage response (DDR) is a common feature of many cancers. The cancer adult T cell leukemia (ATL) can occur in individuals infected with human T cell leukemia virus type 1 (HTLV-1), and ATL cells contain extensive chromosomal abnormalities, suggesting that they have defects in the recognition or repair of DNA damage. Since Tax is the transforming protein encoded by HTLV-1, we asked whether Tax can affect cell cycle checkpoints and the DDR. Using a combination of flow cytometry and DNA repair assays we showed that Tax-expressing cells exit G₁ phase and initiate DNA replication prematurely following damage. Reduced phosphorylation of H2AX (γH2AX) and RPA2, phosphoproteins that are essential to properly initiate the DDR, was also observed in Tax-expressing cells. To determine the cause of decreased DDR protein phosphorylation in Tax-expressing cells, we examined the cellular phosphatase, WIP1, which is known to dephosphorylate γH2AX. We found that Tax can interact with Wip1 in vivo and in vitro, and that Tax-expressing cells display elevated levels of Wip1 mRNA. In vitro phosphatase assays showed that Tax can enhance Wip1 activity on a γH2AX peptide target by 2-fold. Thus, loss of γH2AX in vivo could be due, in part, to increased expression and activity of WIP1 in the presence of Tax. siRNA knockdown of WIP1 in Tax-expressing cells rescued γH2AX in response to damage, confirming the role of WIP1 in the DDR. These studies demonstrate that Tax can disengage the G₁/S checkpoint by enhancing WIP1 activity, resulting in reduced DDR. Premature G₁ exit of Tax-expressing cells in the presence of DNA lesions creates an environment that tolerates incorporation of random mutations into the host genome.

Abstract A85: Novel mouse models to investigate the molecular pathogenesis of metastatic osteosarcoma

Background: Osteosarcoma (OS) is the primary bone tumor in the pediatric population with the main determinant for patient prognosis being the presence of metastatic disease. Approximately 25-30% of patients will present with metastatic osteosarcoma. Patients with metastatic disease have long-term survival rates often <20%.

Objectives: We hypothesize that the molecular pathogenesis of metastatic osteosarcoma is different from localized disease. We intend this work to provide a relevant, endogenous model of metastatic OS that can be utilized to advance our understanding of the molecular pathogenesis of the disease, insights into novel therapeutic targets and as a pre-clinical model for investigating the efficacy of novel therapies.

Design: We have developed a tissue-specific alteration of the p53 status by using osteoblast specific Cre-recombinase expressing mice to generate progeny that spontaneously form endogenous osteosarcomas. Through the use of a mutated, gain of function form of p53, shown previously to be associated with metastatic disease, we have developed a novel immunocompetent model that significantly enhances the endogenous development of metastatic OS.

Results: Tumor analysis has revealed genetic insights in the metastatic progression of osteosarcoma. These include the significant downregulation of Wnt-signaling inhibitors, such as NKD2, APCDD1 and Wnt5a in the metastatic tumors. Functional studies have determined that overexpression of NKD2, also downregulated in several human OS metastatic cell lines, in metastatic mouse OS cell lines leads to a significant decrease in metastatic lung lesions upon transplantation into immunodeficient mice. Possible NKD2-dependent functions include regulation of not only the Wnt signaling pathway, but also blood vessel formation, regulation of cell migration and cell adhesion.

We have also noted the dysregulation of several critical microRNAs in metastatic OS, including the upregulation of mir-130b in the metastatic lesions. This particular microRNA has recently been shown to have clinical correlation in Ewing's sarcoma, with higher levels of tumor mir-130b having a significantly poorer outcome. We are actively pursuing the role of this (and other) microRNA(s) in the molecular pathogenesis of metastatic osteosarcoma.

Conclusion: This novel model has enabled valuable molecular insights into the development and progression of metastatic OS that can lead to the identification of novel therapeutic targets.

Citation Format: Shuying Zhao, Lyazat Kurenbekova, Lawrence A. Donehower, Jason T. Yustein. Novel mouse models to investigate the molecular pathogenesis of metastatic osteosarcoma. [abstract]. In: Proceedings of the AACR Special Conference on Tumor Invasion and Metastasis; Jan 20-23, 2013; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2013;73(3 Suppl):Abstract nr A85.

Rapamycin as longevity enhancer and cancer preventative agent in the context of p53 deficiency

Comment on: 1. Komarova et al. Rapamycin extends lifespan and delays tumorigenesis in heterozygous p53+/- mice. Aging. 2012; 4:10 2. Comas et al. New nanoformulation of rapamycin Rapatar extends lifespan in homozygous p53-/- mice by delaying carcinogenesis. Aging.2012; 4;10.

Mouse tissues that undergo neoplastic progression after K-Ras activation are distinguished by nuclear translocation of phospho-Erk1/2 and robust tumor suppressor responses

Mutation of K-Ras is a frequent oncogenic event in human cancers, particularly cancers of lungs, pancreas, and colon. It remains unclear why some tissues are more susceptible to Ras-induced transformation than others. Here, we globally activated a mutant oncogenic K-Ras allele (K-Ras(G12D)) in mice and examined the tissue-specific effects of this activation on cancer pathobiology, Ras signaling, tumor suppressor, DNA damage, and inflammatory responses. Within 5 to 6 weeks of oncogenic Ras activation, mice develop oral and gastric papillomas, lung adenomas, and hematopoietic hyperproliferation and turn moribund. The oral, gastric, and lung premalignant lesions display activated extracellular signal-regulated kinases (Erk)1/2 and NF-κB signaling as well as activated tumor suppressor and DNA damage responses. Other organs such as pancreas, liver, and small intestine do not exhibit neoplastic progression within 6 weeks following K-Ras(G12D) activation and do not show a potent tumor suppressor response. Even though robust Erk1/2 signaling is activated in all the tissues examined, the pErk1/2 distribution remains largely cytoplasmic in K-Ras(G12D)-refractory tissues (pancreas, liver, and intestines) as opposed to a predominantly nuclear localization in K-Ras(G12D)-induced neoplasms of lung, oral, and gastric mucosa. The downstream targets of Ras signaling, pElk-1 and c-Myc, are elevated in K-Ras(G12D)-induced neoplastic lesions but not in K-Ras(G12D)-refractory tissues. We propose that oncogenic K-Ras-refractory tissues delay oncogenic progression by spatially limiting the efficacy of Ras/Raf/Erk1/2 signaling, whereas K-Ras-responsive tissues exhibit activated Ras/Raf/Erk1/2 signaling, rapidly form premalignant tumors, and activate potent antitumor responses that effectively prevent further malignant progression.

Absence of Wip1 partially rescues Atm deficiency phenotypes in mice

Wild-type p53-induced phosphatase 1 (WIP1) is a serine/threonine phosphatase that dephosphorylates proteins in the ataxia telangiectasia mutated (ATM)-initiated DNA damage response pathway. WIP1 may have a homeostatic role in ATM signaling by returning the cell to a normal pre-stress state following completion of DNA repair. To better understand the effects of WIP1 on ATM signaling, we crossed Atm-deficient mice to Wip1-deficient mice and characterized phenotypes of the double knockout progeny. We hypothesized that the absence of Wip1 might rescue Atm deficiency phenotypes. Atm null mice, like ATM-deficient humans with the inherited syndrome ataxia telangiectasia, exhibit radiation sensitivity, fertility defects, and are T-cell lymphoma prone. Most double knockout mice were largely protected from lymphoma development and had a greatly extended lifespan compared with Atm null mice. Double knockout mice had increased p53 and H2AX phosphorylation and p21 expression compared with their Atm null counterparts, indicating enhanced p53 and DNA damage responses. Additionally, double knockout splenocytes displayed reduced chromosomal instability compared with Atm null mice. Finally, doubly null mice were partially rescued from gametogenesis defects observed in Atm null mice. These results indicate that inhibition of WIP1 may represent a useful strategy for cancer treatment in general and A-T patients in particular.

Rpl27a mutation in the sooty foot ataxia mouse phenocopies high p53 mouse models

Ribosomal stress is an important, yet poorly understood, mechanism that results in activation of the p53 tumour suppressor. We present a mutation in the ribosomal protein Rpl27a gene (sooty foot ataxia mice), isolated through a sensitized N-ethyl-N-nitrosourea (ENU) mutagenesis screen for p53 pathway defects, that shares striking phenotypic similarities with high p53 mouse models, including cerebellar ataxia, pancytopenia and epidermal hyperpigmentation. This phenocopy is rescued in a haploinsufficient p53 background. A detailed examination of the bone marrow in these mice identified reduced numbers of haematopoietic stem cells and a p53-dependent c-Kit down-regulation. These studies suggest that reduced Rpl27a increases p53 activity in vivo, further evident with a delay in tumorigenesis in mutant mice. Taken together, these data demonstrate that Rpl27a plays a crucial role in multiple tissues and that disruption of this ribosomal protein affects both development and transformation.

Mdm2-p53 signaling regulates epidermal stem cell senescence and premature aging phenotypes in mouse skin

The p53 transcription factor is activated by various types of cell stress or DNA damage and induces the expression of genes that control cell growth and inhibit tumor formation. Analysis of mice that express mutant forms of p53 suggest that inappropriate p53 activation can alter tissue homeostasis and life span, connecting p53 tumor suppressor functions with accelerated aging. However, other mouse models that display increased levels of wildtype p53 in various tissues fail to corroborate a link between p53 and aging phenotypes, possibly due to the retention of signaling pathways that negatively regulate p53 activity in these models. In this present study, we have generated mice lacking Mdm2 in the epidermis. Deletion of Mdm2, the chief negative regulator of p53, induced an aging phenotype in the skin of mice, including thinning of the epidermis, reduced wound healing, and a progressive loss of fur. These phenotypes arise due to an induction of p53-mediated senescence in epidermal stem cells and a gradual loss of epidermal stem cell function. These results reveal that activation of endogenous p53 by ablation of Mdm2 can induce accelerated aging phenotypes in mice.

P53 genotype as a determinant of ER expression and tamoxifen response in the MMTV-Wnt-1 model of mammary carcinogenesis

Clinical studies show that estrogen receptor-α (ER) expressing tumors tend to have better prognosis, respond to antiestrogen therapy and have wild-type p53. Conversely, tumors with inactivating mutations in p53 tend to have worse outcomes and to be ER-negative and unresponsive to antihormone treatment. Previous studies from our laboratory have shown that p53 regulates ER expression transcriptionally, by binding the ER promoter and forming a complex with CARM1, CBP, c-Jun, RNA polymerase II and Sp1. In this study, the MMTV-Wnt-1 transgenic mouse model was used to demonstrate that p53 regulation of ER expression and function is not solely an in vitro phenomenon, but it is also operational in mammary tumorigenesis in vivo. The expression of ER and the ability to respond to tamoxifen were determined in mammary tumors arising in p53 wild type (WT) or p53 heterozygous (HT) animals carrying the Wnt-1 transgene. In p53 WT mice, development of ER-positive tumors was delayed by tamoxifen treatment, while tumors arising in p53 HT mice had significantly reduced levels of ER and were not affected by tamoxifen. P53 null tumors were also found in the p53 HT mice and these tumors were ER-negative. ER expression was upregulated in mouse mammary tumor cell lines following transfection with WT p53 or treatment with doxorubicin. These data demonstrate that p53 regulates ER expression in vivo, and affects response to tamoxifen. Results also provide an explanation for the concordant relationship between these prognostic proteins in human breast tumors.

Dual-specificity phosphatase 26 is a novel p53 phosphatase and inhibits p53 tumor suppressor functions in human neuroblastoma

Chemoresistance is a major cause of treatment failure and poor outcome in neuroblastoma. In this study, we investigated the expression and function of dual-specificity phosphatase 26 (DUSP26), also known as mitogen-activated protein kinase phophatase-8, in human neuroblastoma. We found that DUSP26 was expressed in a majority of neuroblastoma cell lines and tissue specimens. Importantly, we found that DUSP26 promotes the resistance of human neuroblastoma to doxorubicin-induced apoptosis by acting as a p53 phosphatase to downregulate p53 tumor suppressor function in neuroblastoma cells. Inhibiting DUSP26 expression in the IMR-32 neuroblastoma cell line enhanced doxorubicin-induced p53 phosphorylation at Ser20 and Ser37, p21, Puma, Bax expression as well as apoptosis. In contrast, DUSP26 overexpression in the SK-N-SH cell line inhibited doxorubicin-induced p53 phosphorylation at Ser20 and Ser37, p21, Puma, Bax expression and apoptosis. Using in vitro and in vivo assays, we found that DUSP26 binds to p53 and dephosphorylates p53 at Ser20 and Ser37. In this report, we show that DUSP26 functions as a p53 phosphatase, which suppresses downstream p53 activity in response to genotoxic stress. This suggests that inhibition of this phosphatase may increase neuroblastoma chemosensitivity and DUSP26 is a novel therapeutic target for this aggressive pediatric malignancy.

Disrupting circadian homeostasis of sympathetic signaling promotes tumor development in mice

Background

Cell proliferation in all rapidly renewing mammalian tissues follows a circadian rhythm that is often disrupted in advanced-stage tumors. Epidemiologic studies have revealed a clear link between disruption of circadian rhythms and cancer development in humans. Mice lacking the circadian genes Period1 and 2 (Per) or Cryptochrome1 and 2 (Cry) are deficient in cell cycle regulation and Per2 mutant mice are cancer-prone. However, it remains unclear how circadian rhythm in cell proliferation is generated in vivo and why disruption of circadian rhythm may lead to tumorigenesis.

Methodology/Principal Findings

Mice lacking Per1 and 2, Cry1 and 2, or one copy of Bmal1, all show increased spontaneous and radiation-induced tumor development. The neoplastic growth of Per-mutant somatic cells is not controlled cell-autonomously but is dependent upon extracellular mitogenic signals. Among the circadian output pathways, the rhythmic sympathetic signaling plays a key role in the central-peripheral timing mechanism that simultaneously activates the cell cycle clock via AP1-controlled Myc induction and p53 via peripheral clock-controlled ATM activation. Jet-lag promptly desynchronizes the central clock-SNS-peripheral clock axis, abolishes the peripheral clock-dependent ATM activation, and activates myc oncogenic potential, leading to tumor development in the same organ systems in wild-type and circadian gene-mutant mice.

Conclusions/Significance

Tumor suppression in vivo is a clock-controlled physiological function. The central circadian clock paces extracellular mitogenic signals that drive peripheral clock-controlled expression of key cell cycle and tumor suppressor genes to generate a circadian rhythm in cell proliferation. Frequent disruption of circadian rhythm is an important tumor promoting factor.

Dephosphorylation of γ-H2AX by WIP1: an important homeostatic regulatory event in DNA repair and cell cycle control

DNA double strand breaks are a particularly toxic form of DNA damage and the mammalian cell has evolved an intricate set of responses to repair this type of DNA lesion. A key early event in the DNA damage response (DDR) is ATM phosphorylation of the histone variant H2AX at serine 139 at the site of the DNA break. Phosphorylated S139 H2AX, or γ-H2AX, forms a docking site for binding of MDC1, leading to sustained recruitment of other DNA repair factors that mediate the repair of the DNA double strand break.  Moreover, recruitment of MDC1 to the break site activates cell cycle checkpoints, protecting the cell from replication of damaged DNA templates. While the molecular events leading to DNA double strand break repair have been well described, the deactivating or homeostatic mechanisms following completion of repair remain largely unexplored. Recent publications by our laboratories and the Medema laboratory shed new light on this issue. Both publications showed that the Wild-type p53-Induced Phosphatase 1 (WIP1) directly dephosphorylates γ-H2AX. WIP1 migrates to the sites of irradiation-induced foci (IRIF), though at a delayed rate relative to MDC1 and mediates γ-H2AX dephosphorylation, presumably after DNA repair is complete. This prevents recruitment of other repair factors such as MDC1 and 53BP1 to the DNA damage sites and promotes the dissolution of IRIF.  In addition, overexpression of WIP1 has a suppressive effect on DNA double strand break repair. Taken together, these reports further implicate WIP1 as a critical homeostatic regulator of the DDR.

55K isoform of CDK9 associates with Ku70 and is involved in DNA repair

Positive elongation factor b (P-TEFb) is a cellular protein kinase that is required for RNA polymerase II (RNAP II) transcriptional elongation of protein coding genes. P-TEFb is a set of different molecular complexes, each containing CDK9 as the catalytic subunit. There are two isoforms of the CDK9 protein - the major 42KDa CDK9 isoform and the minor 55KDa isoform that is translated from an in-frame mRNA that arises from an upstream transcriptional start site. We found that shRNA depletion of the 55K CDK9 protein in HeLa cells induces apoptosis and double-strand DNA breaks (DSBs). The levels of apoptosis and DSBs induced by the depletion were reduced by expression of a 55K CDK9 protein variant resistant to the shRNA, indicating that these phenotypes are the consequence of depletion of the 55K protein and not off-target effects. We also found that the 55K CDK9 protein, but not the 42K CDK9 protein, specifically associates with Ku70, a protein involved in DSB repair. Our findings suggest that the 55K CDK9 protein may function in repair of DNA through an association with Ku70.

Impact papers on aging in 2009

The editorial board of Aging reviews research papers published in 2009, which they believe have or will have a significant impact on aging research. Among many others, the topics include genes that accelerate aging or in contrast promote longevity in model organisms, DNA damage responses and telomeres, molecular mechanisms of life span extension by calorie restriction and pharmacologic interventions into aging. The emerging message in 2009 is that aging is not random but determined by a genetically-regulated longevity network and can be decelerated both genetically and pharmacologically.

Wild-type p53-induced phosphatase 1 dephosphorylates histone variant gamma-H2AX and suppresses DNA double strand break repair

In response to DNA double strand breaks, the histone variant H2AX at the break site is phosphorylated at serine 139 by DNA damage sensor kinases such as ataxia telangiectasia-mutated, forming gamma-H2AX. This phosphorylation event is critical for sustained recruitment of other proteins to repair the break. After repair, restoration of the cell to a prestress state is associated with gamma-H2AX dephosphorylation and dissolution of gamma-H2AX-associated damage foci. The phosphatases PP2A and PP4 have previously been shown to dephosphorylate gamma-H2AX. Here, we demonstrate that the wild-type p53-induced phosphatase 1 (WIP1) also dephosphorylates gamma-H2AX at serine 139 in vitro and in vivo. Overexpression of WIP1 reduces formation of gamma-H2AX foci in response to ionizing and ultraviolet radiation and blocks recruitment of MDC1 (mediator of DNA damage checkpoint 1) and 53BP1 (p53 binding protein 1) to DNA damage foci. Finally, these inhibitory effects of WIP1 on gamma-H2AX are accompanied by WIP1 suppression of DNA double strand break repair. Thus, WIP1 has a homeostatic role in reversing the effects of ataxia telangiectasia-mutated phosphorylation of H2AX.

Widespread and tissue specific age-related DNA methylation changes in mice

Aberrant methylation of promoter CpG islands in cancer is associated with silencing of tumor-suppressor genes, and age-dependent hypermethylation in normal appearing mucosa may be a risk factor for human colon cancer. It is not known whether this age-related DNA methylation phenomenon is specific to human tissues. We performed comprehensive DNA methylation profiling of promoter regions in aging mouse intestine using methylated CpG island amplification in combination with microarray analysis. By comparing C57BL/6 mice at 3-mo-old versus 35-mo-old for 3627 detectable autosomal genes, we found 774 (21%) that showed increased methylation and 466 (13%) that showed decreased methylation. We used pyrosequencing to quantitatively validate the microarray data and confirmed linear age-related methylation changes for all 12 genomic regions examined. We then examined 11 changed genomic loci for age-related methylation in other tissues. Of these, three of 11 showed similar changes in lung, seven of 11 changed in liver, and six of 11 changed in spleen, though to a lower degree than the changes seen in colon. There was partial conservation between age-related hypermethylation in human and mouse intestines, and Polycomb targets in embryonic stem cells were enriched among the hypermethylated genes. Our findings demonstrate a surprisingly high rate of hyper- and hypomethylation as a function of age in normal mouse small intestine tissues and a strong tissue-specificity to the process. We conclude that epigenetic deregulation is a common feature of aging in mammals.

Exercise effects on tumorigenesis in a p53-deficient mouse model of breast cancer

PURPOSE

Physically active women have a reduced risk of breast cancer, but the dose of activity necessary and the role of energy balance and other potential mechanisms have not been fully explored in animal models. We examined treadmill and wheel running effects on mammary tumorigenesis and biomarkers in p53-deficient (p53(+/-)):MMTV-Wnt-1 transgenic mice.

METHODS

Female mice (9 wk old) were randomly assigned to the following groups in experiment 1: treadmill exercise 5 d x wk(-1), 45 min x d(-1), 5% grade at 20 m x min(-1), approximately 0.90 km x d(-1) (TREX1, n = 20) or at 24 m x min(-1), approximately 1.08 km x d(-1) (TREX2, n = 21); or a nonexercise control (CON-TREX, n = 22). In experiment 2, mice were randomly assigned to voluntary wheel running (WHL, n = 21, 2.46 +/- 1.11 km x d(-1) (mean +/- SD)) or to a nonexercise control (CON-WHL, n = 22). Body composition was measured at approximately 9 wk and serum insulin-like growth factor 1 (IGF-1) at two to three monthly time points beginning at approximately 9 wk on study. Mice were sacrificed when tumors reached 1.5 cm, mice became moribund, or there was only one mouse per treatment group remaining.

RESULTS

TREX1 (24 wk) and TREX2 (21 wk) had shorter median survival times than CON-TREX (34 wk; P < 0.01), whereas those of WHL and CON-WHL were similar (23 vs 24 wk; P = 0.32). TREX2 had increased multiplicity of mammary gland carcinomas compared with CON-TREX; WHL had a higher tumor incidence than CON-WHL. All exercising animals were lighter than their respective controls, and WHL had lower body fat than CON-WHL (P < 0.01). There was no difference in IGF-1 between groups (P > 0.05).

CONCLUSIONS

Despite beneficial or no effects on body weight, body fat, or IGF-1, exercise had detrimental effects on tumorigenesis in this p53-deficient mouse model of spontaneous mammary cancer.

Using mice to examine p53 functions in cancer, aging, and longevity

The p53 tumor suppressor is a multifaceted transcription factor that responds to a diverse array of stresses that include DNA damage and aberrant oncogene signaling. On activation, p53 prevents the emergence of cancer cells by initiating cell cycle arrest, senescence (terminal cell cycle arrest), or apoptosis. Although its role in assuring longevity by suppressing cancer is well established, recent studies obtained largely from genetically engineered mouse models suggest that p53 may regulate longevity and aging. In some contexts, it appears that altered p53 activity may enhance longevity, and in others, it appears to suppress longevity and accelerate aging phenotypes. Here, we discuss how genetically engineered mouse models have been used to explore antiproliferative functions of p53 in cancer suppression and how mouse models with altered aging phenotypes have shed light on how p53 might influence the aging process.

20 years studying p53 functions in genetically engineered mice

Cell and molecular biological studies of p53 functions over the past 30 years have been complemented in the past 20 years by studies that use genetically engineered mice. As expected, mice that have mutant Trp53 alleles usually develop cancers of various types more rapidly than their counterparts that have wild-type Trp53 genes. These mouse studies have been instrumental in providing important new insights into p53 tumour suppressor function. Such studies have been facilitated by the development of increasingly sophisticated genetic engineering approaches, which allow the more precise manipulation of p53 structure and function in a mammalian model.