Cross-species genomics identifies DLG2 as a tumor suppressor in osteosarcoma

Oncogene, Aug 2018

Leveraging the conserved cancer genomes across mammals has the potential to transform driver gene discovery in orphan cancers. Here, we combine cross-species genomics with validation across human–dog–mouse systems to uncover a new bone tumor suppressor gene. Comparative genomics of spontaneous human and dog osteosarcomas (OS) expose Disks Large Homolog 2 (DLG2) as a tumor suppressor candidate. DLG2 copy number loss occurs in 42% of human and 56% of canine OS. Functional validation through pertinent human and canine OS DLG2-deficient cell lines identifies a regulatory role of DLG2 in cell division, migration and tumorigenesis. Moreover, osteoblast-specific deletion of Dlg2 in a clinically relevant genetically engineered mouse model leads to acceleration of OS development, establishing DLG2 as a critical determinant of OS. This widely applicable cross-species approach serves as a platform to expedite the search of cancer drivers in rare human malignancies, offering new targets for cancer therapy.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://www.nature.com/articles/s41388-018-0444-4.pdf

Cross-species genomics identifies DLG2 as a tumor suppressor in osteosarcoma

Abstract Leveraging the conserved cancer genomes across mammals has the potential to transform driver gene discovery in orphan cancers. Here, we combine cross-species genomics with validation across human–dog–mouse systems to uncover a new bone tumor suppressor gene. Comparative genomics of spontaneous human and dog osteosarcomas (OS) expose Disks Large Homolog 2 (DLG2) as a tumor suppressor candidate. DLG2 copy number loss occurs in 42% of human and 56% of canine OS. Functional validation through pertinent human and canine OS DLG2-deficient cell lines identifies a regulatory role of DLG2 in cell division, migration and tumorigenesis. Moreover, osteoblast-specific deletion of Dlg2 in a clinically relevant genetically engineered mouse model leads to acceleration of OS development, establishing DLG2 as a critical determinant of OS. This widely applicable cross-species approach serves as a platform to expedite the search of cancer drivers in rare human malignancies, offering new targets for cancer therapy. Introduction Extensive genomic alterations in orphan cancers make it difficult to identify driver mutations. Osteosarcoma (OS) presents a major challenge due to its chaotic genome and limited patient populations. Cytogenetic and genomic studies report a high degree of chromothripsis, chromosomal instability, copy number aberrations, and structural variations in OS [1,2,3,4,5]. Multiple patient cohorts have undergone aCGH amounting to ~100 OS genomes, and whole-genome sequencing was applied to 34 OS cases [6]. OS, the most common bone malignancy in adolescents, severely compromises the quality of life in children and adolescents. Further, despite improvements in orthopedic surgery and combination chemotherapy over the past 4 decades that have enhanced patient outcome, the 5-year survival rate of OS has plateaued at ~65% [7]. New strategies are urgently required to uncover the molecular determinants of this cancer in order to design novel targeted therapies. Naturally occurring canine cancers offer a cross-species approach to understanding the pathogenesis of human cancers due to the considerable homology between dog and human genomes. OS is one of the most frequent malignancies in dogs and shares remarkable clinical similarity with the human disease including radiological and histological features, and the pattern of metastatic spread. Compared to mice, dogs are relatively out-bred and share the same living environment as humans. Publication of high-resolution dog reference genomes now permit a detailed genomic comparison between these species [8,9,10,11]. We created a novel comparative oncogenomics platform in this study, compiling high-resolution human and newly generated canine OS aCGH data, to identify known and novel drivers in bone cancer. Discs large homolog 2 (DLG2) was highly mutated in both human and dog, and was validated as a tumor suppressor through a workflow that spanned human and canine cell lines, and a clinically relevant murine model of OS. Results and discussion To create a comparative oncogenomics platform (Fig. 1a), we selected a published dataset of 52 human OS samples which contain 40 tumor specimens and 12 OS cell lines (GSE12789). In parallel, 9 canine OS specimens representing a number of dog breeds were subjected to aCGH (Supplemental Table S1). Among the multiple tools used to segment the human and canine data, we found the greatest concordance between Circular Binary Segmentation (CBS) and Gain and Loss Analysis of DNA (GLAD); these were used to define the copy number profile [12, 13]. Segmented data were next analyzed by GISTIC (FDR Q-value < 0.1) to shortlist recurrent genomic regions for each species [14]. While candidate genes in human were directly produced by GISTIC, those in dog were pinpointed through cross species synteny matching of 38 canine chromosomes to the human reference genome, hg19. Human and canine OS genomes contained high levels of chromosomal instability marked by frequent whole chromosomal arm gains and losses (Fig. 1b). Mutational frequency of recurrent CNAs across samples showed that among the top eight highly significant genes, six were common between human and canine OS (Fig. 1c and Supplemental Fig. 1A), providing proof-of-principle for our approach. These included known oncogenes (MYC) and tumor suppressors (CDKN2A/B, RB1, and PTEN) indicating the high level of disease similarity and the overlapping mutational landscape across these species [15,16,17,18]. We recorded a novel focal deletion in Chr11q14 across multiple cases, harboring potential tumor suppressor gene(s) (Fig. 1c). Specifically, 42% (22/52; Chr11q14) of human OS and 55.6% (5/9; Chr21) of canine OS exhibited this focal loss (Supplemental Fig. 1B). This frequency is consistent with that reported for next-generation sequencing of 34 human OS cases (52.9%, 18/34) [6], although its biological significance was not pursued. We noted that DLG2 is the only gene residing in the focal deletion in human (chr11: 83835663–84011662: 180 kb) and in dog (chr21:17457665–17858164: 400 kb). The relative distribution of DLG2 deletion and other common alterations across human and dog OS is shown in Fig. 1d. Fig. 1 Whole-genome landscape of CNAs in human and canine OS identifies DLG2 locus. a A schematic of cross-species oncogenomics strategy begins with comparison of curated human and new dog aCGH data on OS. Recurrent events serve as a filter to limit the number of candidate genes and selected new gene candidate undergoes biological validation through gain-of-function and loss-of-function studies designed to assess its impact on tumorigenic properties in vitro and oncogenic or tumor suppressor function in vivo. b Segmented CNA profiles of 52 human OS genome (left panel) and 9 canine OS genome (right panel) were mapped. X-axis presents individual samples. Y-axis lists chromosomal locations from Chr1 to Chr22 for human, or Chr1 to Chr38 for dog. c Significant deletion regions calculated by GISTIC across the human and canine OS genome. Important genes in the peaks are labeled. X-axis: significant score, Y-axis: chromosomal locations. d GRID heat map showing copy number status of top mutated genes. Samples are ordered by DLG2 status. e Of 31 human OS samples, 9 exhibited DLG2 deletion in an independent cohort. f Pattern of point mutations in DLG2 genes based on data from COMSIC. Most mutations result in amino acid changes. g Mutation hot spots across DLG2 amino acid sequences Full size image We then acquired an independent cohort of 31 human OS samples to examine the DLG2 status using genomic qPCR (courtesy of Dr. David Malkin). Of 31 samples, 9 displayed DLG2 deletion (29%; Fig. 1e), providing a validation in this independent patient cohort. Next, interrogation of the Catalogue of Somatic Mutations in Cancer (COSMIC) database exposed DLG2 to be highly mutated in a variety of human cancers where most of the mutations resulted in drastic amino acid changes (Fig. 1f and Supplemental Table 2) [19]. DLG2 is a member of the Membrane-Associated Guanylate Kinase (MAGUK) protein family that is characterized by concatenations of PDZ, SH3 and GUK domains. These domains enable MAGUKs to function as scaffold proteins that orchestrate the assembly of multiple signal transduction networks [20]. The mammalian DLG family of proteins are primarily known for their role in epithelial polarity and polarity during cell division [21]. Several mutational hotspots were seen: V161 and V191 within the first PDZ domain; and S624 in the connector sequence adjacent to the SH3 domain (Fig. 1g). To examine the impact of restoring DLG2 expression on tumorigenic properties, we sought relevant human and canine cell lines that naturally harbored DLG2 deletions. Data from GSE12789 and GSE36003 revealed that DLG2 was deleted in SAOS-2, SJSA-1 and 2 primary OS cell lines (HuO-3N1 & HuO9), while DLG2 loss was also detected by aCGH analysis in a canine OS cell line that we generated in house (OVC-cOSA-31); Supplemental Fig. 1C). Consistently, endogenous DLG2 protein expression was very low in these cells (Supplemental Fig. 1D). A human DLG2 expression vector was then constructed for the dominant, longest isoform1 (Supplemental Fig. 1D). DLG2 expression did not affect proliferation in normal or low-serum conditions (Supplemental Fig. 2A). However, when tested in anchorage-independent soft agar and 3D-Matrigel assays, DLG2 overexpression in SJSA-1 or OVC-COSA-31 cells significantly reduced their colony formation capacity (Fig. 2a, b). HuO9 and SAOS-2 did not grow in anchorage-independent conditions. We then tested the effect of DLG2 restoration in vivo. Luciferase-tagged [22], parental SJSA-1 and OVC-COSA-31 cell lines formed xenografts when injected subcutaneously into NOD-SCID mice, whereas DLG2 restored cell lines showed slower xenograft growth, resulting in smaller tumors (Fig. 2c, d). Fig. 2 DLG2 restoration hinders cell and tumor growth in in vitro 3D cultures and xenograft experiments. a Representative pictures of colonies in 3D-matrigel. DLG2 re-expression induced a significant reduction in both colony number and size. Quantification of colony area is shown on the right. b Representative soft agar images and their respective quantifications. DLG2 restoration reduces colony forming capacities. c NOD-SCID mice injected with luciferase- expressing SJSA-1 parental (left) and DLG2 (right) cells imaged at 3 weeks post injection (n = 10 per group). Photon influx value was used as a surrogate marker for tumor size. Photon influx of each tumor was quantified and average values are shown on the right. d Tumors taken from NOD-SCID mice 4 weeks after innoculation with OVC-COSA-31 parental (left) and DLG2 (right) canine OS cells (n = 10 per group). Tumor size was calculated and is shown on the right. White bar represents 1 cm. Values shown are the mean ± SD of three separate determinations. *P < 0.05 by two-tailed Student’s t test Full size image We next investigated the gene expression profiles associated with DLG2 deletion in canine OS. Microarrays generated from 7 aCGH-matched canine primary tumors (4 DLG2 deleted, 3 DLG2 wild type) showed 167 differentially expressed probes/genes (DEGs) in the two groups (Fig. 3a, b, Supplemental Table 3). Interestingly, G protein-coupled receptor (GPCR) signaling was one of the top altered pathways (Fig. 3c). The functional interactions between PDZ domain-containing proteins such as DLG1 and DLG4, and GPCRs have been implicated in many diseases including cancer [23]. Indeed, biological process prediction of DLG2 using Integrative Multi-species Prediction (IMP) database (http://imp.princeton.edu) showed that DLG2 is also involved in GPCR signaling pathway (GO term ID: GO:0086103, GO:0007186, GO:0008277). Importantly, GTPase-activating proteins RGS9 (Regulator of G-protein signaling 9) [24] and GARNL3 (GTPase activating Rap/RanGAP domain-like 3) [25] called as GAPs were among the most downregulated genes between DLG2-deleted and DLG2-wild-type groups (Fig. 3b). GAPs are a family of regulatory proteins that bind to activated G proteins and result in “switching off” of GTPase signaling [26]. Therefore, RGS9 and GARNL3 downregulation in DLG2-deleted tumor may lead to increased GTPase signaling activity. Fig. 3 DLG2 is involved in small GTPase activity and its restoration reduces migratory capacities. a Differentially expressed probes/genes in Dlg2-deleted (n = 4, blue) vs. Dlg2-wild-type (n = 3, orange) groups. b Volcano plots show differential probes/genes identified using T-test with FDR correction. c Top altered pathways were analyzed using Enrichr online tool. G- protein-coupled receptor protein signaling pathway was among the top significant pathways. d Representative images of HuO9 and OVC-COSA-31 control and DLG2-re-expressing cells 24 h after a scratch-wound. Quantification of blank areas is shown on the right. e Representative images of Hu-O9 and OVC-COSA-31 cells 24 h after seeding in transwells. Transwell insert pore size of 8 and 12 μm were used. Quantification of cells that migrated through the transwell is shown on the right. f Cell cycle analysis showing DLG2 induced cell cycle changes in OVC-COSA-31 and Hu-O9 cells cultured in ultra-low attachment plate. DAPI (x-axis) was used to determine DNA content and BrdU (Y-axis) was used to determine nucleotide uptake. DLG2 restoration resulted in significant accumulation of cells in G2/M and S phases. Quantification of cells at specific stages of the cell cycle is shown on the right. Values shown are the mean ± SD of three separate determinations. *P < 0.05 by two-tailed Student’s t test Full size image Small GTPase signaling pathway, such as Rho and Ran, has been previously implicated in cytoskeletal re-organization and migration [27]. We therefore performed scratch-wound and transwell migration assays. DLG2-expressing cell lines showed significantly decreased motility in both assays compared to their respective parental control cells (Fig. 3d, e). This phenotype is in agreement with our gene expression data indicating a negative regulatory role for DLG2 in the GTPase signaling pathway, which in turn influences human and canine OS cell migration. Since GTPase signaling also orchestrates multiple stages of the cell cycle, we asked whether DLG2 status affects the cell cycle using flow cytometry [28]. We observed more cells in either G2/M or S phases when DLG2-re-expressing cells were grown under low-attachment conditions (Fig. 3f) but not when grown as a monolayer (data not shown). This observation together with slower proliferation under 3D, but not 2D conditions (Fig. 2a, b, and Supplemental Fig. 2A) suggests that DLG2 participates in regulating cell cycle progression. It is likely that the tumor suppressor function of DLG2 is mediated in part through modulating the small GTPase signaling pathway. The single Drosophila DLG homolog is a tumor suppressor [29]. To examine whether mammalian DLG2 functions as a tumor suppressor, we utilized a genetically engineered mouse model (GEMM) of OS. RB and TP53 mutations are the most common events in OS and germline mutations also predispose humans to this cancer. It has been shown that combined osteoblast-specific deletion of p53 and Rb, achieved through the Collagen Type 1 Alpha-1 (Col1a1) promoter-driven Cre recombination in mice leads to spontaneous OS development [30, 31]. We bred Dlg2-floxed mice into the above clinically relevant model to generate triple osteoblast-specific deletions (p53ΔOB/ΔOBRbΔOB/ΔOBDlg2ΔOB/ΔOB: Dlg2−/−, Fig. 4a). Dlg2 homozygous deletion significantly accelerated tumor development and shortened survival in mice (Fig. 4b–d). Specifically, the median lifespan for the Dlg2−/− cohort was 184 days vs. 239 days in the control group (p53ΔOB/ΔOBRbΔOB/ΔOBDlg2WT/WT: Dlg2+/+). The anatomical distribution and histological features of OS tumors developed in Dlg2−/− mice remained comparable to littermate controls. The abundant osteoid deposition at the invasive the front of thoracic vertebra osteosarcoma (Fig. 4e) illustrates the aggressive nature, osteoblastic nature of tumors developed in these mice. Histological features of osteosarcoma, such as plump tumorous osteoblasts and copious amount of osteoid were also seen in tumors that developed at different sites in the Dlg2-null mice (Fig. 4f). Moreover, tumor cells derived from Dlg2-deficient mice were more proliferative than those of wild-type mice, as indicated by Ki-67 staining (Fig. 4g). Fig. 4 Osteoblast-specific Dlg2 knockout reduces overall survival and accelerates tumor onset. a Breeding strategy for generating osteoblast-specific compound deletion of Dlg2, p53, and Rb1. b Representative X-ray and c MicroCT images. At 22 weeks of age, Dlg2 wild type mice show no visible signs of tumor while Dlg2 conditional knockout mice developed tumors on tibia and rib. d Kaplan–Meier curve of mice with different Dlg2 status. Dlg2 homozygous conditional knockout shortens overall survival compared to Dlg2 wild type. e H&E images of a thoracic vertebra osteosarcoma. At higher magnification, osteosarcoma cells (arrows) and osteoid deposition (asterisk) were seen invasive fronts. f H&E images of a femoral osteosarcoma located on the Ilium. Note the plump tumorous osteoblasts (arrows) and massive osteoid deposition (asterisk). g Ki-67 staining of indicated tumor tissues. Scale bar: 100 μm Full size image Altogether, through cross-species oncogenomics we have uncovered DLG2 deletion to be a highly frequent event in human (42%) and dog (56%) OS. We extended our cross-species approach to molecular and functional analyses of DLG2 by utilizing pertinent human and dog OS cell lines as well as p53/Rb-driven OS GEMM. These series of studies establish DLG2 as a bone tumor suppressor and provide insight into its involvement in GTPase signaling and regulation of cell cycle and proliferation in a 3D tissue culture model. This work presents a blue print for cross-species comparative genomics platform as a useful tool for driver gene discovery in rare human cancers. Methods aCGH analysis Published human OS genomic datasets on Gene Expression Omnibus (GEO) were selected (GSE12789). Canine OS specimens were obtained from University of Guelph veterinarian clinics. Genomic DNA was extracted using phenol chloroform and treated with RNaseA (Invitrogen) for 1 h at 37 °C. Canine aCGH assay was performed by SurePrint G3 Canine CGH Microarray 4 × 180 K (Agilent) according to manufacturers’ specifications at the University Health Network Microarray Centre. CBS and GLAD were used for segmentation of the genome and GISTIC was used to determine the common altered region and the significance of genetic events. Canine aCGH data were deposited in the GEO database (GSE111637). Gene expression microarray analysis Total RNA from canine OS tumors was prepared using the RNeasy mini-kit (QIAGEN), RNA concentration and integrity were assessed using Nanodrop 2000 and Bioanalyzer (Thermo Scientific and Agilent). Gene expression microarray on Canine (V2) Gene Expression Microarray, 4 × 44k (Agilent) were performed at the University Health Network Microarray Centre. Multi-array average (RMA) normalization was used across the sample sets and batch correction was performed. Differentially expressed genes were identified using T-test with FDR correction (P < 0.05 and FDR < 0.1) or Samtools. Analyses were performed using R (version 2.82), Agilent Genespring GX (version 11) or BRBArrayTools. Gene ontology and pathway analysis were performed using Enrichr. Total Ensemble genes were used as background. Microarray data were deposited in the GEO database (GSE111638). Mice Osteoblast-specific Dlg2 knockout OS mouse model (p53ΔOB/ΔOB RbΔOB/ΔOBDlg2ΔOB/ΔOB) was generated by breeding Dlgflox/+ (Dlg2tm1a(EUCOMM)Wtsi, obtained from Wellcome Trust Sanger Institute) with Col1a1-Cre-p53flox/floxRbflox/flox (n = 12). p53ΔOB/ΔOBRbΔOB/ΔOBDlg2wt/wt littermates (n = 10) were used as the control. X-ray and MicroCT scan were carried out as previously reported [28]. For xenograft experiments, 8-week-old female recipient non-obese diabetic/severe combined immunodeficiency (NOD-SCID) mice were obtained from Charles River Laboratories. Mice were assigned to experimental groups by stratified randomization according to body weight. OVC-COSA-31 and luciferase-expressing SJSA-1 cells were injected subcutaneously on both flanks (n = 10 per group) and monitored weekly for 4–8 weeks by electronic caliper or bioluminescence imaging using Xenogen IVIS Imaging System 100 (STTARR facility, UHN). Investigators who measured and analyzed the samples were blind to the group information. Statistics Two-tailed Student’s t-test was used to determine if two sets of normally distributed data are significantly different from each other. Biological replicates were performed independently at least three times. Log-rank test was performed for survival curves. For all studies, P < 0.05 was considered significant. Specialized analyses for bioinformatic data, including aCGH analysis and gene expression microarray analysis, were performed using dedicated tools mentioned above. Study approval All animal protocols were approved by the University Health Network Animal Care Committee and performed in accordance with the standards of the Canadian Council on Animal Care. Mice were housed in the Princess Margaret Hospital Animal Care Facility. User protocol AUP 849.1 was strictly followed. Written informed consent was obtained from the patients or their guardians before sample collection. References 1. dos Santos Aguiar S, de Jesus Girotto Zambaldi L, dos Santos AM, Pinto W Jr, Brandalise SR. Comparative genomic hybridization analysis of abnormalities in chromosome 21 in childhood osteosarcoma. Cancer Genet Cytogenet. 2007;175:35–40. ArticlePubMed Google Scholar2. Nathrath MH, Kuosaite V, Rosemann M, Kremer M, Poremba C, Wakana S, et al. Two novel tumor suppressor gene loci on chromosome 6q and 15q in human osteosarcoma identified through comparative study of allelic imbalances in mouse and man. Oncogene. 2002;21:5975–80. ArticlePubMed Google Scholar3. Kuijjer ML, Rydbeck H, Kresse SH, Buddingh EP, Lid AB, Roelofs H, et al. Identification of osteosarcoma driver genes by integrative analysis of copy number and gene expression data. Genes Chromosomes Cancer. 2012;51:696–706. ArticlePubMed Google Scholar4. Chen X, Bracht JR, Goldman AD, Dolzhenko E, Clay DM, Swart EC, et al. The architecture of a scrambled genome reveals massive levels of genomic rearrangement during development. Cell. 2014;158:1187–98. ArticlePubMedPubMed Central Google Scholar5. Stephens PJ, Greenman CD, Fu B, Yang F, Bignell GR, Mudie LJ, et al. Massive genomic rearrangement acquired in a single catastrophic event during cancer development. Cell. 2011;144:27–40. ArticlePubMedPubMed Central Google Scholar6. Chen X, Bahrami A, Pappo A, Easton J, Dalton J, Hedlund E, et al. Recurrent somatic structural variations contribute to tumorigenesis in pediatric osteosarcoma. Cell Rep. 2014;7:104–12. ArticlePubMedPubMed Central Google Scholar7. Aljubran AH, Griffin A, Pintilie M, Blackstein M. Osteosarcoma in adolescents and adults: survival analysis with and without lung metastases. Ann Oncol. 2009;20:1136–41. ArticlePubMed Google Scholar8. Tsai PC, Breen M. Array-based comparative genomic hybridization-guided identification of reference genes for normalization of real-time quantitative polymerase chain reaction assay data for lymphomas, histiocytic sarcomas, and osteosarcomas of dogs. Am J Vet Res. 2012;73:1335–43. ArticlePubMed Google Scholar9. Angstadt AY, Thayanithy V, Subramanian S, Modiano JF, Breen M. A genome-wide approach to comparative oncology: high-resolution oligonucleotide aCGH of canine and human osteosarcoma pinpoints shared microaberrations. Cancer Genet. 2012;205:572–87. ArticlePubMed Google Scholar10. Angstadt AY, Motsinger-Reif A, Thomas R, Kisseberth WC, Guillermo Couto C, Duval DL, et al. Characterization of canine osteosarcoma by array comparative genomic hybridization and RT-qPCR: signatures of genomic imbalance in canine osteosarcoma parallel the human counterpart. Genes Chromosomes Cancer. 2011;50:859–74. ArticlePubMed Google Scholar11. Scott MC, Sarver AL, Gavin KJ, Thayanithy V, Getzy DM, Newman RA, et al. Molecular subtypes of osteosarcoma identified by reducing tumor heterogeneity through an interspecies comparative approach. Bone. 2011;49:356–67. ArticlePubMedPubMed Central Google Scholar12. Venkatraman ES, Olshen AB. A faster circular binary segmentation algorithm for the analysis of array CGH data. Bioinformatics. 2007;23:657–63. ArticlePubMed Google Scholar13. Hupe P, Stransky N, Thiery JP, Radvanyi F, Barillot E. Analysis of array CGH data: from signal ratio to gain and loss of DNA regions. Bioinformatics. 2004;20:3413–22. ArticlePubMed Google Scholar14. Beroukhim R, Getz G, Nghiemphu L, Barretina J, Hsueh T, Linhart D, et al. Assessing the significance of chromosomal aberrations in cancer: methodology and application to glioma. Proc Natl Acad Sci USA. 2007;104:20007–12. ArticlePubMed Google Scholar15. Shimizu T, Ishikawa T, Sugihara E, Kuninaka S, Miyamoto T, Mabuchi Y, et al. c-MYC overexpression with loss of Ink4a/Arf transforms bone marrow stromal cells into osteosarcoma accompanied by loss of adipogenesis. Oncogene. 2010;29:5687–99. ArticlePubMed Google Scholar16. Toguchida J, Ishizaki K, Sasaki MS, Nakamura Y, Ikenaga M, Kato M, et al. Preferential mutation of paternally derived RB gene as the initial event in sporadic osteosarcoma. Nature. 1989;338:156–8. ArticlePubMed Google Scholar17. Freeman SS, Allen SW, Ganti R, Wu J, Ma J, Su X, et al. Copy number gains in EGFR and copy number losses in PTEN are common events in osteosarcoma tumors. Cancer. 2008;113:1453–61. ArticlePubMedPubMed Central Google Scholar18. Mohseny AB, Tieken C, van der Velden PA, Szuhai K, de Andrea C, Hogendoorn PC, et al. Small deletions but not methylation underlie CDKN2A/p16 loss of expression in conventional osteosarcoma. Genes Chromosomes Cancer. 2010;49:1095–103. ArticlePubMed Google Scholar19. Forbes SA, Beare D, Gunasekaran P, Leung K, Bindal N, Boutselakis H, et al. COSMIC: exploring the world’s knowledge of somatic mutations in human cancer. Nucleic Acids Res. 2015;43:D805–811. ArticlePubMed Google Scholar20. Craven SE, Bredt DS. PDZ proteins organize synaptic signaling pathways. Cell. 1998;93:495–8. ArticlePubMed Google Scholar21. Roberts S, Delury C, Marsh E. The PDZ protein discs-large (DLG): the ‘Jekyll and Hyde’ of the epithelial polarity proteins. FEBS J. 2012;279:3549–58. ArticlePubMed Google Scholar22. Ray P, De A, Min JJ, Tsien RY, Gambhir SS. Imaging tri-fusion multimodality reporter gene expression in living subjects. Cancer Res. 2004;64:1323–30. ArticlePubMedPubMed Central Google Scholar23. Dunn HA, Ferguson SS. PDZ protein regulation of G protein-coupled receptor trafficking and signaling pathways. Mol Pharmacol. 2015;88:624–39. ArticlePubMed Google Scholar24. Hooks SB, Waldo GL, Corbitt J, Bodor ET, Krumins AM, Harden TK. RGS6, RGS7, RGS9, and RGS11 stimulate GTPase activity of Gi family G-proteins with differential selectivity and maximal activity. J Biol Chem. 2003;278:10087–93. ArticlePubMed Google Scholar25. van Dam TJ, Bos JL, Snel B. Evolution of the Ras-like small GTPases and their regulators. Small GTPases. 2011;2:4–16. ArticlePubMedPubMed Central Google Scholar26. Ross EM, Wilkie TM. GTPase-activating proteins for heterotrimeric G proteins: regulators of G protein signaling (RGS) and RGS-like proteins. Annu Rev Biochem. 2000;69:795–827. ArticlePubMed Google Scholar27. Stramer B, Wood W, Galko MJ, Redd MJ, Jacinto A, Parkhurst SM, et al. Live imaging of wound inflammation in Drosophila embryos reveals key roles for small GTPases during in vivo cell migration. J Cell Biol. 2005;168:567–73. ArticlePubMedPubMed Central Google Scholar28. Coleman ML, Marshall CJ, Olson MF. RAS and RHO GTPases in G1-phase cell-cycle regulation. Nat Rev Mol Cell Biol. 2004;5:355–66. ArticlePubMed Google Scholar29. Woods DF, Hough C, Peel D, Callaini G, Bryant PJ. Dlg protein is required for junction structure, cell polarity, and proliferation control in Drosophila epithelia. J Cell Biol. 1996;134:1469–82. ArticlePubMed Google Scholar30. Walkley CR, Qudsi R, Sankaran VG, Perry JA, Gostissa M, Roth SI, et al. Conditional mouse osteosarcoma, dependent on p53 loss and potentiated by loss of Rb, mimics the human disease. Genes Dev. 2008;22:1662–76. ArticlePubMedPubMed Central Google Scholar31. Chen Y, Di Grappa MA, Molyneux SD, McKee TD, Waterhouse P, Penninger JM, et al. RANKL blockade prevents and treats aggressive osteosarcomas. Sci Transl Med. 2015;7:317ra197. ArticlePubMed Google Scholar Download references Acknowledgements We thank all members of our laboratory for helpful discussions. Funding This study is supported by funding from the Canadian Institutes of Health Research (CIHR), Ontario Institute for Cancer Research (OICR), and Global Leadership Round in Genomics and Life Sciences (GL2) to RK and the OVC Pet Trust to GW. J.Lu is also supported by National Natural Science Foundation of China (81702665), National Postdoctoral Program for Innovative Talents (BX201600196) and China Postdoctoral Science Foundation (2017M610571). Author information Author notes These authors contributed equally: Yang W. Shao, Geoffrey A. Wood, Jinchang Lu. AffiliationsPrincess Margaret Cancer Centre, University of Toronto, Toronto, Ontario, M5G 1L7, CanadaYang W. Shao, Jinchang Lu, Qing-Lian Tang, Sam Molyneux, Yan Chen, Hui Fang, Trevor McKee, Paul Waterhouse & Rama KhokhaDepartment of Pathobiology, Ontario Veterinary College, University of Guelph, Guelph, Ontario, N1G 2W1, CanadaGeoffrey A. Wood, Jonathan Liu & Hibret AdissuDepartment of Orthopedic Oncology, the First Affiliated Hospital, Sun Yat-Sen Univeristy, 510080, Guangzhou, ChinaJinchang Lu & Qing-Lian Tang AuthorsSearch for Yang W. Shao in:Nature Research journals • PubMed • Google Scholar Search for Geoffrey A. Wood in:Nature Research journals • PubMed • Google Scholar Search for Jinchang Lu in:Nature Research journals • PubMed • Google Scholar Search for Qing-Lian Tang in:Nature Research journals • PubMed • Google Scholar Search for Jonathan Liu in:Nature Research journals • PubMed • Google Scholar Search for Sam Molyneux in:Nature Research journals • PubMed • Google Scholar Search for Yan Chen in:Nature Research journals • PubMed • Google Scholar Search for Hui Fang in:Nature Research journals • PubMed • Google Scholar Search for Hibret Adissu in:Nature Research journals • PubMed • Google Scholar Search for Trevor McKee in:Nature Research journals • PubMed • Google Scholar Search for Paul Waterhouse in:Nature Research journals • PubMed • Google Scholar Search for Rama Khokha in:Nature Research journals • PubMed • Google Scholar Conflict of interest The authors declare that they have no conflict of interest. Corresponding author Correspondence to Rama Khokha. Electronic supplementary material Supplemental data Rights and permissions Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. About this article Publication history Received 19 March 2018 Revised 22 June 2018 Accepted 24 July 2018 Published 09 August 2018 DOI https://doi.org/10.1038/s41388-018-0444-4


This is a preview of a remote PDF: https://www.nature.com/articles/s41388-018-0444-4.pdf

Yang W. Shao, Geoffrey A. Wood, Jinchang Lu, Qing-Lian Tang, Jonathan Liu, Sam Molyneux, Yan Chen, Hui Fang, Hibret Adissu, Trevor McKee, Paul Waterhouse, Rama Khokha. Cross-species genomics identifies DLG2 as a tumor suppressor in osteosarcoma, Oncogene, 2018, DOI: 10.1038/s41388-018-0444-4