Advertisement
Regular article Tumorigenesis and neoplastic progression| Volume 188, ISSUE 12, P2902-2911, December 2018

Characterization of a Clival Chordoma Xenograft Model Reveals Tumor Genomic Instability

  • Roberto J. Diaz
    Correspondence
    Address correspondence to Roberto Jose Diaz, M.D., Ph.D., F.R.C.S.C., Department of Neurology and Neurosurgery, Montreal Neurological Institute and Hospital, 3801 Rue University, Room 109, Montreal, Quebec H3A 2B4, Canada.
    Affiliations
    Department of Neurology and Neurosurgery, Montreal Neurological Institute and Hospital, McGill University, Montreal, Quebec, Canada
    Search for articles by this author
  • Amanda Luck
    Affiliations
    Division of Neurosurgery, The Hospital for Sick Children, Toronto, Ontario, Canada
    Search for articles by this author
  • Andrew Bondoc
    Affiliations
    Division of Neurosurgery, The Hospital for Sick Children, Toronto, Ontario, Canada
    Search for articles by this author
  • Brian Golbourn
    Affiliations
    Division of Neurosurgery, The Hospital for Sick Children, Toronto, Ontario, Canada

    Department of Laboratory Medicine and Pathobiology, University of Toronto, Toronto, Ontario, Canada
    Search for articles by this author
  • Daniel Picard
    Affiliations
    Department of Pediatric Neuro-Oncogenomics, German Cancer Consortium (DKTK) and German Cancer Research Center (DKFZ), Heidelberg, Germany

    Department of Pediatric Oncology, Hematology, and Clinical Immunology, Medical Faculty, University Hospital Düsseldorf, Düsseldorf, Germany
    Search for articles by this author
  • Marc Remke
    Affiliations
    Department of Pediatric Neuro-Oncogenomics, German Cancer Consortium (DKTK) and German Cancer Research Center (DKFZ), Heidelberg, Germany

    Department of Pediatric Oncology, Hematology, and Clinical Immunology, Medical Faculty, University Hospital Düsseldorf, Düsseldorf, Germany

    Department of Neuropathology, Medical Faculty, Heinrich-Heine University Düsseldorf, Düsseldorf, Germany
    Search for articles by this author
  • James Loukides
    Affiliations
    Division of Neurosurgery, The Hospital for Sick Children, Toronto, Ontario, Canada
    Search for articles by this author
  • Nesrin Sabha
    Affiliations
    Division of Neurosurgery, The Hospital for Sick Children, Toronto, Ontario, Canada
    Search for articles by this author
  • Christian Smith
    Affiliations
    Division of Neurosurgery, The Hospital for Sick Children, Toronto, Ontario, Canada
    Search for articles by this author
  • Michael D. Cusimano
    Affiliations
    Division of Neurosurgery, St. Michael's Hospital, Toronto, Ontario, Canada

    Division of Neurosurgery, Department of Surgery, University of Toronto, Toronto, Ontario, Canada
    Search for articles by this author
  • James T. Rutka
    Affiliations
    Division of Neurosurgery, The Hospital for Sick Children, Toronto, Ontario, Canada

    Department of Laboratory Medicine and Pathobiology, University of Toronto, Toronto, Ontario, Canada

    Division of Neurosurgery, Department of Surgery, University of Toronto, Toronto, Ontario, Canada
    Search for articles by this author
Open ArchivePublished:September 21, 2018DOI:https://doi.org/10.1016/j.ajpath.2018.08.004
      Patient-derived xenografts retain the genotype of the parent tumors more readily than tumor cells maintained in culture. The two previously reported clival chordoma xenografts were derived from recurrent tumors after radiation. To study the genetics of clival chordoma in the absence of prior radiation exposure we established a patient-derived xenograft at primary resection of a clival chordoma. Epicranial grafting of clival chordoma collected during surgery was performed. Tumor growth was established in a nonobese diabetic/severe combined immunodeficiency mouse and tumors have been passaged serially for seven generations. Physaliferous cell architecture was shown in the regenerated tumors, which stained positive for Brachyury, cytokeratin, and S100 protein. The tumors showed bone invasion. Single-nucleotide polymorphism analysis of the tumor xenograft was compared with the parental tumor. Copy number gain of the T gene (brachyury) and heterozygous loss of cyclin dependent kinase inhibitor 2A (CDKN2A) was observed. Heterozygous loss of the tumor-suppressor fragile histidine triad (FHIT) gene also was observed, although protein expression was preserved. Accumulation of copy number losses and gains as well as increased growth rate was observed over three generations. The patient-derived xenograft reproduces the phenotype of clival chordoma. This model can be used in the future to study chordoma biology and to assess novel treatments.
      Chordoma is a very rare and slowly growing tumor arising in the ventral axial skeleton. Tumor cells express Brachyury, a transcription factor that is restricted topographically to nascent and migrating mesoderm and the notochord during development.
      • Herrmann B.G.
      Expression pattern of the Brachyury gene in whole-mount TWis/TWis mutant embryos.
      • Wilkinson D.G.
      • Bhatt S.
      • Herrmann B.G.
      Expression pattern of the mouse T gene and its role in mesoderm formation.
      Notochordal cell rests presently are hypothesized as precursors to chordoma.
      • Salisbury J.R.
      • Deverell M.H.
      • Cookson M.J.
      • Whimster W.F.
      Three-dimensional reconstruction of human embryonic notochords: clue to the pathogenesis of chordoma.
      • Shen J.
      • Li C.D.
      • Yang H.L.
      • Lu J.
      • Zou T.M.
      • Wang D.L.
      • Deng M.
      Classic chordoma coexisting with benign notochordal cell rest demonstrating different immunohistological expression patterns of Brachyury and galectin-3.
      • Yamaguchi T.
      • Suzuki S.
      • Ishiiwa H.
      • Ueda Y.
      Intraosseous benign notochordal cell tumours: overlooked precursors of classic chordomas?.
      Approximately a third of all patients with a diagnosis of chordoma have tumor located at the skull base.
      • McMaster M.L.
      • Goldstein A.M.
      • Bromley C.M.
      • Ishibe N.
      • Parry D.M.
      Chordoma: incidence and survival patterns in the United States, 1973-1995.
      These patients more commonly are younger in age and present most frequently with cranial nerve deficits.
      • Diaz R.J.
      • Maggacis N.
      • Zhang S.
      • Cusimano M.D.
      Determinants of quality of life in patients with skull base chordoma.
      The median survival for patients with clival chordoma is 9.2 years.
      • Chambers K.J.
      • Lin D.T.
      • Meier J.
      • Remenschneider A.
      • Herr M.
      • Gray S.T.
      Incidence and survival patterns of cranial chordoma in the United States.
      These tumors are exceptionally resistant to standard DNA-altering chemotherapy, limiting the effectiveness of adjuvant therapies subsequent to surgical resection.
      • Diaz R.J.
      • Cusimano M.D.
      The biological basis for modern treatment of chordoma.
      To better understand the tumor response to specific adjuvant therapy, experimental models that allow the assessment of tumor cells in a living microenvironment are needed. Presently, four chordoma xenografts (two clival and two sacral) derived from primary tumor, without cell culture, have been reported in the literature (Table 1).
      • Davies J.M.
      • Robinson A.E.
      • Cowdrey C.
      • Mummaneni P.V.
      • Ducker G.S.
      • Shokat K.M.
      • Bollen A.
      • Hann B.
      • Phillips J.J.
      Generation of a patient-derived chordoma xenograft and characterization of the phosphoproteome in a recurrent chordoma.
      • Siu I.M.
      • Salmasi V.
      • Orr B.A.
      • Zhao Q.
      • Binder Z.A.
      • Tran C.
      • Ishii M.
      • Riggins G.J.
      • Hann C.L.
      • Gallia G.L.
      Establishment and characterization of a primary human chordoma xenograft model.
      • Owen J.H.
      • Komarck C.M.
      • Wang A.C.
      • Abuzeid W.M.
      • Keep R.F.
      • McKean E.L.
      • Sullivan S.
      • Fan X.
      • Prince M.E.P.
      UM-Chor1: establishment and characterization of the first validated clival chordoma cell line.
      • Bozzi F.
      • Manenti G.
      • Conca E.
      • Stacchiotti S.
      • Messina A.
      • Dagrada G.
      • Gronchi A.
      • Panizza P.
      • Pierotti M.A.
      • Tamborini E.
      • Pilotti S.
      Development of transplantable human chordoma xenograft for preclinical assessment of novel therapeutic strategies.
      • Trucco M.M.
      • Awad O.
      • Wilky B.A.
      • Goldstein S.D.
      • Huang R.
      • Walker R.L.
      • Shah P.
      • Katuri V.
      • Gul N.
      • Zhu Y.J.
      • McCarthy E.F.
      • Paz-Priel I.
      • Meltzer P.S.
      • Austin C.P.
      • Xia M.
      • Loeb D.M.
      A novel chordoma xenograft allows in vivo drug testing and reveals the importance of NF-kappaB signaling in chordoma biology.
      • Karikari I.O.
      • Gilchrist C.L.
      • Jing L.
      • Alcorta D.A.
      • Chen J.
      • Richardson W.J.
      • Gabr M.A.
      • Bell R.D.
      • Kelley M.J.
      • Bagley C.A.
      • Setton L.A.
      Molecular characterization of chordoma xenografts generated from a novel primary chordoma cell source and two chordoma cell lines.
      • Hsu W.
      • Mohyeldin A.
      • Shah S.R.
      • ap Rhys C.M.
      • Johnson L.F.
      • Sedora-Roman N.I.
      • Kosztowski T.A.
      • Awad O.A.
      • McCarthy E.F.
      • Loeb D.M.
      • Wolinsky J.P.
      • Gokaslan Z.L.
      • Quinones-Hinojosa A.
      Generation of chordoma cell line JHC7 and the identification of Brachyury as a novel molecular target.
      • Presneau N.
      • Shalaby A.
      • Ye H.
      • Pillay N.
      • Halai D.
      • Idowu B.
      • Tirabosco R.
      • Whitwell D.
      • Jacques T.S.
      • Kindblom L.G.
      • Bruderlein S.
      • Moller P.
      • Leithner A.
      • Liegl B.
      • Amary F.M.
      • Athanasou N.N.
      • Hogendoorn P.C.
      • Mertens F.
      • Szuhai K.
      • Flanagan A.M.
      Role of the transcription factor T (brachyury) in the pathogenesis of sporadic chordoma: a genetic and functional-based study.
      • Liu X.
      • Nielsen G.P.
      • Rosenberg A.E.
      • Waterman P.R.
      • Yang W.
      • Choy E.
      • Sassi S.
      • Yang S.
      • Harmon D.C.
      • Yang C.
      • Schwab J.H.
      • Kobayashi E.
      • Mankin H.J.
      • Xavier R.
      • Weissleder R.
      • Duan Z.
      • Hornicek F.J.
      Establishment and characterization of a novel chordoma cell line: CH22.
      • DeComas A.M.
      • Penfornis P.
      • Harris M.R.
      • Meyer M.S.
      • Pochampally R.R.
      Derivation and characterization of an extra-axial chordoma cell line (EACH-1) from a scapular tumor.
      These tumors were established in the flanks of immunodeficient mice. The clival chordoma xenografts previously described have been derived from tumors with prior radiation exposure and therefore may harbor genetic aberrations that are not representative of the originating oncogenic mutations in the tumor. Therefore, we present the first clival chordoma primary xenograft derived from previously untreated tumor.
      Table 1Summary of Chordoma Xenograft Models
      Study (year)NameSourceBrachyuryS100CKGenomic
      Present studySMH5Clival PDX+ (nuc)NT+SNP
      SNP Affymetrix 6.0.
      Davies et al
      • Davies J.M.
      • Robinson A.E.
      • Cowdrey C.
      • Mummaneni P.V.
      • Ducker G.S.
      • Shokat K.M.
      • Bollen A.
      • Hann B.
      • Phillips J.J.
      Generation of a patient-derived chordoma xenograft and characterization of the phosphoproteome in a recurrent chordoma.
      (2014)
      SF8894Clival PDX+ (nuc)NT+None
      Siu et al
      • Siu I.M.
      • Salmasi V.
      • Orr B.A.
      • Zhao Q.
      • Binder Z.A.
      • Tran C.
      • Ishii M.
      • Riggins G.J.
      • Hann C.L.
      • Gallia G.L.
      Establishment and characterization of a primary human chordoma xenograft model.
      (2012)
      JHH-2009-011Clival PDX+ (nuc)++SNP
      SNP Illumina 600W (Illumina, San Diego, CA).
      Owen et al
      • Owen J.H.
      • Komarck C.M.
      • Wang A.C.
      • Abuzeid W.M.
      • Keep R.F.
      • McKean E.L.
      • Sullivan S.
      • Fan X.
      • Prince M.E.P.
      UM-Chor1: establishment and characterization of the first validated clival chordoma cell line.
      (2017)
      UM-Chor1Clival CLD+ (cyto)NTNTNone
      + (nuc)
      Bozzi et al
      • Bozzi F.
      • Manenti G.
      • Conca E.
      • Stacchiotti S.
      • Messina A.
      • Dagrada G.
      • Gronchi A.
      • Panizza P.
      • Pierotti M.A.
      • Tamborini E.
      • Pilotti S.
      Development of transplantable human chordoma xenograft for preclinical assessment of novel therapeutic strategies.
      (2014)
      4Sacral PDX+ (nuc)NTNTFISH
      Trucco et al
      • Trucco M.M.
      • Awad O.
      • Wilky B.A.
      • Goldstein S.D.
      • Huang R.
      • Walker R.L.
      • Shah P.
      • Katuri V.
      • Gul N.
      • Zhu Y.J.
      • McCarthy E.F.
      • Paz-Priel I.
      • Meltzer P.S.
      • Austin C.P.
      • Xia M.
      • Loeb D.M.
      A novel chordoma xenograft allows in vivo drug testing and reveals the importance of NF-kappaB signaling in chordoma biology.
      (2013)
      Sacral PDX+ (cyto)NTNTSNP
      HumanCNV370-Duo version 1.0 (Illumina).
      Karikari et al
      • Karikari I.O.
      • Gilchrist C.L.
      • Jing L.
      • Alcorta D.A.
      • Chen J.
      • Richardson W.J.
      • Gabr M.A.
      • Bell R.D.
      • Kelley M.J.
      • Bagley C.A.
      • Setton L.A.
      Molecular characterization of chordoma xenografts generated from a novel primary chordoma cell source and two chordoma cell lines.
      (2014)
      DVC-4Sacral CLD+ (nuc)NTNTPCR at brachy locus
      Hsu et al
      • Hsu W.
      • Mohyeldin A.
      • Shah S.R.
      • ap Rhys C.M.
      • Johnson L.F.
      • Sedora-Roman N.I.
      • Kosztowski T.A.
      • Awad O.A.
      • McCarthy E.F.
      • Loeb D.M.
      • Wolinsky J.P.
      • Gokaslan Z.L.
      • Quinones-Hinojosa A.
      Generation of chordoma cell line JHC7 and the identification of Brachyury as a novel molecular target.
      (2011)
      JHC7Sacral CLD+ (nuc)++Karyotype
      G-banding
      Presneau et al
      • Presneau N.
      • Shalaby A.
      • Ye H.
      • Pillay N.
      • Halai D.
      • Idowu B.
      • Tirabosco R.
      • Whitwell D.
      • Jacques T.S.
      • Kindblom L.G.
      • Bruderlein S.
      • Moller P.
      • Leithner A.
      • Liegl B.
      • Amary F.M.
      • Athanasou N.N.
      • Hogendoorn P.C.
      • Mertens F.
      • Szuhai K.
      • Flanagan A.M.
      Role of the transcription factor T (brachyury) in the pathogenesis of sporadic chordoma: a genetic and functional-based study.
      (2011)
      U-CH1Sacral CLD+ (nuc)NTNTFISH
      Liu et al
      • Liu X.
      • Nielsen G.P.
      • Rosenberg A.E.
      • Waterman P.R.
      • Yang W.
      • Choy E.
      • Sassi S.
      • Yang S.
      • Harmon D.C.
      • Yang C.
      • Schwab J.H.
      • Kobayashi E.
      • Mankin H.J.
      • Xavier R.
      • Weissleder R.
      • Duan Z.
      • Hornicek F.J.
      Establishment and characterization of a novel chordoma cell line: CH22.
      (2012)
      CH22Sacral CLD+ (nuc)++None
      DeComas et al
      • DeComas A.M.
      • Penfornis P.
      • Harris M.R.
      • Meyer M.S.
      • Pochampally R.R.
      Derivation and characterization of an extra-axial chordoma cell line (EACH-1) from a scapular tumor.
      (2010)
      EACH-1Scapular CLD+ (nuc)++Karyotype
      CK, cytokeratin; CLD, cell-line–derived xenograft; cyto, cytoplasmic; FISH, fluorescence in situ hybridization; Genomic, type of genomic analysis conducted; NT, not tested; nuc, nuclear; PDX, patient-derived xenograft; SNP, single-nucleotide polymorphism; +, positive immunostaining.
      SNP Affymetrix 6.0.
      SNP Illumina 600W (Illumina, San Diego, CA).
      HumanCNV370-Duo version 1.0 (Illumina).
      The patient-derived chordoma xenograft has been characterized histopathologically and by molecular analysis to determine if phenotype and genotype are maintained in serial generations. We previously identified frequent loss of fragile histidine triad (FHIT) protein expression in chordoma and postulated a potential role in chordoma pathogenesis.
      • Diaz R.J.
      • Guduk M.
      • Romagnuolo R.
      • Smith C.A.
      • Northcott P.
      • Shih D.
      • Berisha F.
      • Flanagan A.
      • Munoz D.G.
      • Cusimano M.D.
      • Pamir M.N.
      • Rutka J.T.
      High-resolution whole-genome analysis of skull base chordomas implicates FHIT loss in chordoma pathogenesis.
      Therefore, FHIT expression and genetic aberration at the FHIT locus were investigated in this patient-derived chordoma xenograft. FHIT is a well-characterized tumor suppressor and a potential target for therapy.
      • Diaz R.J.
      • Guduk M.
      • Romagnuolo R.
      • Smith C.A.
      • Northcott P.
      • Shih D.
      • Berisha F.
      • Flanagan A.
      • Munoz D.G.
      • Cusimano M.D.
      • Pamir M.N.
      • Rutka J.T.
      High-resolution whole-genome analysis of skull base chordomas implicates FHIT loss in chordoma pathogenesis.

      Materials and Methods

      Tumor Harvesting and Initial Implantation

      Patient consent for use of tumor material for research purposes was obtained and tissue collection was conducted as per approval by the St. Michael's Hospital Institutional Review Board. The patient was a 68-year-old man who presented with a 1-year history of neck pain and tongue weakness. Preoperative magnetic resonance imaging and computed tomography scan of the skull base showed a central mass in the lower clivus, extending to the level of the C2 vertebra. No evidence of metastatic disease was found. Heterogenous contrast enhancement, T1-weighted hypointensity, T2-weighted hyperintensity, and bony erosion were observed (Figure 1). An extended image-guided endoscopic transnasal transphenoidal surgical resection was performed. After ensuring sufficient tumor tissue had been obtained for pathologic diagnosis, several pieces of tumor were wrapped in sterile gauze and placed in a partially filled specimen container containing saline to keep the tissue moist. The consulting pathologist confirmed the diagnosis of classic chordoma, not otherwise specified. The tumor showed integrase interactor 1 and p53 expression by immunohistochemistry. The patient had rapid progression of disease after surgical treatment.
      Figure thumbnail gr1
      Figure 1Imaging features of originating clival chordoma. A: Axial T1-weighted magnetic resonance image (MRI) showing a midline clival mass (T). B: T1-weighted gadolinium-enhanced MRI showing a tumor (T) with heterogeneous enhancement. C: Axial T2-weighted MRI showing a hyperintense microcystic epidural tumor (T) in the region of the lower clivus. D: Axial noncontrast computed tomography showing bony erosion of the clivus (arrowhead).

      Initial Implantation Procedure

      The specimen container was transported on ice. Under a biological containment hood the tumor sample was divided and a portion was frozen in liquid nitrogen and stored at −80°C. The remainder of the tissue was minced (≤1-mm fragments) and suspended in 750 μL of phosphate-buffered saline (PBS) and kept cold on ice. Two nonobese diabetic/severe combined immunodeficiency γ mice (The Jackson Laboratory, Bar Harbor, ME) were anesthetized by isoflurane gas and a small transverse incision (approximately 4 mm) was made over the posterior parietal skull and suboccipital musculature (n = 2) or over the flank (n = 1). The parietal bone was scraped with the scalpel to thin the outer cortex. The 250 μL of the minced tissue suspension was mixed with 250 μL of Matrigel (Becton Dickinson, Franklin Lakes, NJ) and then implanted by wide-bore pipette into the subcutaneous epicranial space above the posterior parietal bone and suboccipital musculature or into the subcutaneous space over the flank. The incisions were closed with absorbable suture and the mice recovered. Serial inspection of the wound and implantation site was performed on a daily basis for 1 week and then on a weekly basis. Mice received 5 mg/kg ketoprofen analgesic and a bolus of 0.5 mL 0.9% saline subcutaneously in the immediate period after surgery.

      Serial Passage of Tumors and Growth Assessment

      Animal experiments were approved by the Hospital for Sick Children Animal Care Committee (AUP0204-H) and conducted in accordance with the Ontario Animals for Research Act and the Canadian Council for Animal Care guidelines. Palpable tumors were measured with a digital caliper to estimate the tumor volume (in mm3) using the following formula: tumor volume = length (mm) × width2 (mm2)/2.
      Serial transplantation of tumors was performed when the tumor exceeded 1.5 cm in maximum diameter. Mice were anesthetized using isoflurane gas and the tumor was excised and immediately placed into cold PBS for mincing. The tumor bulk was divided into four equal portions: one portion was frozen in liquid nitrogen, one portion was fixed in 4% formaldehyde or cryopreserved in liquid nitrogen isopentane for histologic analysis, and the remaining two portions were resuspended in 1 mL PBS with 1 mL Matrigel. From the PBS:Matrigel mix, 400 μL was implanted into the epicranial space above the posterior parietal bone and suboccipital musculature of anesthetized nonobese diabetic/severe combined immunodeficiency γ mice (n = 4) as per the procedure described in Initial Implantation Procedure. This process was continued for subsequent generations of tumor growth with two modifications. After P1, the implantation ratio was changed to one tumor to three mice for P2 to P4. From P5 and onward, the implantation ratio was one tumor to two mice in order to achieve faster tumor regeneration. The implantation volume was changed to a minced tissue suspension of 75 μL PBS and 75 μL Matrigel to allow for a smaller epicranial pocket. Tumor growth was monitored on a weekly basis and tumor volume measurements were initiated every 15 days when tumors reached a minimum of 0.2 mm in maximum diameter. Mice were euthanized for collection of tumor if any pain, severe disability, or tumor maximum diameter greater than 1.5 cm was reached in accordance with institutional humane animal use protocols. Upon death, the dorsal subcutaneous space was opened and inspected for invasion of the skull by tumor.

      Immunohistochemical Analysis

      Paraffin tissue sections (5 μm) were deparaffinized, rehydrated, and pretreated in citrate buffer, pH 6.0, for 15 minutes. Sections were blocked in 10% goat serum and incubated in primary antibodies for 1 hour at room temperature: Brachyury (1:100, 04-135; Millipore, Burlington, MA), FHIT (1:100, HPA018909; Sigma, St. Louis, MO), pan-cytokeratin (1:100, NBP2-29429; Novus Bio, Littleton, CO), Ki67 (1:50, CRM 325 A; Biocare Medical, Pacheco, CA), INI1 (1:50, 612110; BD Biosciences), and p53 (1:50, 554294; BD Biosciences). Flash-frozen OCT-embedded tumor tissue sections (8 μm) were blocked in 10% goat serum and incubated with anti-nuclei antibody clone 235-1 (1:100, MAB1812; Millipore)
      • Bissig-Choisat B.
      • Wang L.
      • Legras X.
      • Saha P.K.
      • Chen L.
      • Bell P.
      • Pankowicz F.P.
      • Hill M.C.
      • Barzi M.
      • Kettlun Leyton C.
      • Leung H.C.
      • Kruse R.L.
      • Himes R.W.
      • Goss J.A.
      • Wilson J.M.
      • Chan L.
      • Lagor W.R.
      • Bissig K.D.
      Development and rescue of human familial hypercholesterolaemia in a xenograft mouse model.
      for 1 hour at room temperature. This antibody recognizes human nuclei, but not mouse or rat nuclei. After incubation with primary antibody and washing three times, biotinylated secondary antibody 1:100 (ABC kit; Vector Labs, Burlington, Ontario, Canada) was applied for 30 minutes at room temperature. Sections then were exposed to the avidin-biotin detection system for 40 minutes at room temperature. After three washes, signal was detected using diaminobenzidine, counterstained in hematoxylin, and mounted in Permount (Thermo Fisher Scientific, Waltham, MA). Micrographs were captured with an Infinity1 camera (Lumenera Corporation, Ottawa, Ontario, Canada) with Infinity1 software visualized through an Olympus BX43 light microscope (Olympus, Tokyo, Japan). The Ki67 index was calculated by counting the number of positive nuclei in 10 high-power fields (20× objective).

      Whole-Genome Single-Nucleotide Polymorphism Profiling

      Genomic DNA was isolated from tumor samples and snap-frozen in liquid nitrogen at the time of passaging. A proteinase K/SDS buffer was used to dissociate the tissue (10 mmol/L Tris-Cl [pH 8.0], 0.1 mol/L EDTA [pH 8.0], 0.5% (w/v) SDS, 20 mg/mL proteinase K) overnight at 55°C. Protein extraction was achieved with phenol-chloroform-isoamyl alcohol using phase-lock centrifuge gels (5-Prime, Gaithersburg, MD) and DNA was precipitated with 100% ethanol and 7.5 mol/L ammonium acetate. The pellet was washed with 75% ethanol and then dried. It was resuspended in Tris-ethanolamine (pH 8.0) buffer and quantified by NanoDrop (NanoDrop, Willmington, DE). Microarray processing was performed in The Centre for Applied Genomics, The Hospital for Sick Children. Briefly, genomic DNA was cleaved with Nsp and Sty nucleases, biotinylated, and hybridized to the Genome-Wide Human Single-Nucleotide Polymorphism (SNP) Array 6.0 (Affymetrix, Santa Clara, CA). Scanning of the microarray was performed on a Affymetrix GeneChip Scanner 3000. Copy number data were analyzed using Partek Genomic Suite 6.6 v6.15.1207 (Partek, St. Louis, MO). Samples were imported, background was normalized using Partek Defaults, and the copy number was calculated from allele-specific intensity. Genomic segmentation was performed using the following parameters: i) 150 minimum genomic markers were required; ii) P value threshold was 0.001; iii) signal-to-noise ratio was 0.3; and iv) diploid status was in the range of 1.7 and 2.3, losses had a copy number <1.7, and gains had a copy number >2.3. Copy number alteration counts were determined by summing the segmentation results <1.7 and >2.3.
      Data sets GSE25720 and GSE9023 were downloaded for comparative analysis from GEO at the NCBI website (https://www.ncbi.nlm.nih.gov/geo). GSE25720 samples were analyzed on the Genome-Wide Human SNP Array 6.0 and data analysis was performed as described in the previous paragraph. For GSE9023 samples, which were processed using the klingenkhh BAC 32K Full array, normalized log2 ratios were imported into Partek Genomic Suite to be visualized. Genomic segmentation was performed using the following parameters: i) 10 minimum genomic markers were required; ii) P value threshold was 0.001; iii) signal-to-noise ratio was 0.3; and iv) diploid status was in the range of −0.15 and 0.15, losses had a copy number <−0.15, and gains had a copy number >0.15.

      RNAscope in Situ Hybridization

      The RNAscope in situ hybridization (ACDBio, Newark, CA) assay was performed on formalin-fixed, paraffin-embedded sections from the originating patient-derived xenograft (PDX0) and third-generation (PDX3) tumors. Briefly, tissue sections were deparaffinized and treated with hydrogen peroxide, followed by RNA target retrieval and protease treatment. After the sections were prepared, a commercial probe against human T gene mRNA was used to hybridize to the target, followed by staining of the targets with RNAscope 2.5 HD Red Detection reagent (ACDBio). Slides were counterstained with hematoxylin and mounted using EcoMount (BioCare Medical, Pacheco, CA). Whole slide scans were acquired using an Aperio AT2 Scanner (Leica Biosystems, Wetzlar, Germany) at ×40 magnification. Images then were analyzed using Definiens Developer XD software version 2.6 (Definiens, Munich, Germany), using a custom-programed algorithm to count the number of spots per nucleus (indicating single RNA molecules) in the entire tumor region.

      Statistical Analysis

      Prism 7 (GraphPad Software, Inc., La Jolla, CA) was used for statistical analysis. All measures are reported as means ± SEM. Means for survival time, tumor diameter, and Ki67 index were compared by an unpaired Welch t-test with a P value <0.05 set as significant. RNAscope data were compared using negative binomial regression.

      Results

      Human Clival Chordoma Can Be Regenerated in Nonobese Diabetic/Severe Combined Immunodeficiency Mouse Epicranial Space

      Clival chordoma tissue obtained at the time of surgical resection (Figure 1) was implanted in a Matrigel (Corning, Corning, NY) slurry within the epicranial space over the posterior parietal bone and in the subcutaneous space of the flank of nonobese diabetic/severe combined immunodeficiency mice. After 4.5 months a palpable and visible tumor formed in the epicranial space in 2 of 2 mice (Figure 2). No tumor formed in the subcutaneous space at the flank of the third host mouse by 4.5 months. The tumor from the epicranial space from one of the mice was passaged serially and is now in the 7th generation (Table 2). The patient-derived xenograft has been made available to the research community through the Chordoma Foundation (https://www.chordomafoundation.org, last accessed September 4, 2017). On gross pathology, bony invasion was observed in each generation from PDX3 to PDX7; however, it was not apparent for every tumor within each generation (Table 2). Microscopic bony invasion consisted of nests of tumor cells surrounded by mesenchymal cells (Figure 2). A pattern of successive growth rate increase over successive generations was observed as the tumor to host transplantation ratio was changed from 1:4 to 1:2. The median time to the appearance of the first tumor measuring >0.2 mm in maximum diameter was 104 ± 12 days (n = 15) in PDX4 (1:3 transplant ratio) compared with 44 ± 11 days (n = 10) in PDX7 (1:2 transplant ratio) (P = 0.001). The estimated volumetric tumor growth rates showed a wide range of variability in the PDX4 compared with the PDX7 generation (Figure 3A). The Ki67 labeling index was significantly higher in the PDX7 generation (10% ± 2%) compared with parental tumor (3% ± 1%; P = 0.002) (Figure 3B). Along with the faster tumor growth a decrease in mouse survival was observed over generations, with mean survival times of 35 ± 3 weeks (n = 15) in PDX4 compared with 27 ± 1 week (n = 10) in PDX7 (P = 0.026).
      Figure thumbnail gr2
      Figure 2A: Vivisection of xenograft tumor (arrowheads) overlying the calvarium and extension to subcutaneous tissue overlying the thoracic spine. B: Hematoxylin and eosin–stained chordoma xenograft tumor (arrowheads) invading into mouse skull bone. Scale bar = 50 μm.
      Table 2Tally of Tumors Formed in Each Generation after Serial Host-To-Host Implantation
      GenerationMice implanted, nTumors formed, nTime to first measurable tumor >0.2 mm length, daysTumors with gross bony invasion, n
      PDX0221260
      PDX144890
      PDX2552240
      PDX399514
      PDX41715754
      PDX599305
      PDX61010306
      PDX71010157
      Figure thumbnail gr3
      Figure 3A: Tumor growth curves for individual mice over time since implantation of tumor xenograft. Growth in the mice in the fourth generation (PDX4) and seventh generation (PDX7) is compared showing earlier tumor formation in PDX7. B: Histologic sections of original tumor, PDX4, and PDX7 showing Ki67 immunohistochemistry. Graph quantifying the Ki67 index in original tumor, PDX4, and PDX7. The Ki67 index is higher in PDX7 tumor compared with original tumor. **P < 0.01. Scale bars = 50 μm.

      Expression of Chordoma Markers over Xenograft Generations

      The xenograft tumor showed similar cytology as the original tumor with lobules of physaliferous cells (Figure 4). Nuclear Brachyury expression and cytokeratin expression was conserved over generations from original tumor to PDX3 (Figure 4). The original and PDX3 tumor did not express epithelial membrane antigen, as expected for a chordoma (Figure 4). The identity of the chordoma xenograft was confirmed by immunohistochemistry independently by Valasciences, Inc. (San Diego, CA). The physaliferous cells of the patient-derived xenograft expressed human FHIT in a cytoplasmic and perinuclear location, although cytoplasmic expression was prominent (Figure 5). FHIT expression was conserved over the first three generations. Absent anti-human FHIT immunolabeling of connective tissue components in the xenografts was observed, as would be expected owing to the antibody specificity (Figure 5). Concordantly, anti-human nucleus antibody staining of physaliferous cells, but not the connective tissue component of the xenograft tumors, was observed (Figure 5).
      Figure thumbnail gr4
      Figure 4Histologic sections of original tumor and third-generation xenograft (PDX3) showing immunohistochemistry positive for Pan-cytokeratin (Pan-CK) and nuclear Brachyury, and absence of epithelial membrane antigen (EMA). Tumor cytology showing lobules of physalipherous cells surrounded by a thin layer of spindle-like cells is preserved in the xenograft on hematoxylin and eosin (H&E) staining. Scale bars = 40 μm.
      Figure thumbnail gr5
      Figure 5A: Histologic section of parental chordoma immunolabeled with anti-human FHIT antibody. B: Histologic section of third-generation xenograft immunolabeled with anti-human FHIT antibody. Prominent immunolabeling with anti-human FHIT antibody is seen in physaliferous tumor cells (arrow), but not in adjacent mouse connective tissue (arrowheads). C: Histologic section of U87 human glioblastoma xenograft (negative control) immunolabeled with anti-human FHIT. D: Histologic section of human kidney (positive control) immunolabeled with anti-human FHIT. E: Histologic section of eighth-generation xenograft immunolabeled with anti-human nuclear antibody identifying physaliferous tumor cells as human and associated connective tissue (arrowheads) as originating from mouse. F: Histologic section of mouse muscle (negative control) immunolabeled with anti-human nuclear antibody. Scale bars = 20 μm.

      Chordoma Genome Accumulates Copy Number Gains and Losses over Generations

      To show the resemblance of the original tumor to other clival chordomas a karyogram was generated from the SNP data collected (Supplemental Figure S1). Copy number analysis shows gains and losses throughout the tumor genome (Figure 6A). The tumor established in the immunocompromised mouse at PDX0 and PDX3 carries the same key copy number aberrations as the original patient tumor, namely heterozygous loss of chromosome 9 (which includes the CDKN2A locus) and copy number gain of a region of chromosome 6 encompassing the T gene (brachyury) (Supplemental Figure S2 and Figure 6B). The PDX3 tumor showed significantly higher levels of T gene labeling per nucleus as compared with PDX0 (P < 0.001) (Figure 6C). Heterozygous loss of FHIT was observed over the first three generations (Supplemental Figure S3). The major genomic difference between the original tumor and the xenograft tumor are chromosome 19 losses (Supplemental Figure S4). Furthermore, copy number gains in chromosome 21 are lost in the transition from PDX0 to PDX3 (Supplemental Figure S5). The presence of gains or losses in 150 probe segments across the genome increased as follows: original tumor (n = 323), PDX0 (n = 937), and PDX3 (n = 1407). Unlike the MUG-Chor1 cell line
      • Owen J.H.
      • Komarck C.M.
      • Wang A.C.
      • Abuzeid W.M.
      • Keep R.F.
      • McKean E.L.
      • Sullivan S.
      • Fan X.
      • Prince M.E.P.
      UM-Chor1: establishment and characterization of the first validated clival chordoma cell line.
      parental tumor, which does not show T gene copy number gain, the parental tumor in this study showed T gene copy number gain and this persisted in PDX1 and PDX3 (Supplemental Figure S6). MUG-Chor1 cells in passages 5 and 13 showed acquisition of T gene copy number gain. Among 11 chordoma tumor samples studied in GSE9023, only two samples showed T gene copy number gain. No highly recurrent copy number change was observed across tumor, cell line, and patient-derived xenograft samples. There was significant variation in copy number gains and losses across samples; however, chromosomal alterations in 1p, 3, 4, 7, 9p, 10, 12, 14, 19, and 22 were frequent. SNP data from the present study are available at (https://www.ncbi.nlm.nih.gov/geo; accession number GSE103342).
      Figure thumbnail gr6
      Figure 6A: Copy number losses (blue) and gains (red) across chromosomes in parental tumor (O), originating xenograft (PDX0), and third-generation xenograft (PDX3). B: Chromosomal map showing copy number gain of a segment of chromosome 6, which includes the T gene. C: RNAscope analysis using a probe targeting the T gene. Each spot indicates a single RNA molecule. A custom algorithm was used to count the number of spots per nucleus. Nuclei and spots were segmented by Definiens XD version 2.6 and delineated. PDX3 showed significantly higher (P < 0.001) levels of T gene labeling compared with PDX0, thus indicating copy gain of the gene. Scale bar = 10 μm. Chr6, chromosome 6.

      Discussion

      There is growing interest in establishing validated models for investigation of novel therapeutics for chordoma, a tumor that is highly resistant to standard chemotherapy.
      • Diaz R.J.
      • Cusimano M.D.
      The biological basis for modern treatment of chordoma.
      Here, we report the growth and serial passage in a nonobese diabetic/severe combined immunodeficiency mouse of a patient-derived clival chordoma that had not been irradiated previously. This is the third patient-derived clival chordoma model reported in the literature; however, it is the first derived from a previously nonirradiated tumor.
      • Davies J.M.
      • Robinson A.E.
      • Cowdrey C.
      • Mummaneni P.V.
      • Ducker G.S.
      • Shokat K.M.
      • Bollen A.
      • Hann B.
      • Phillips J.J.
      Generation of a patient-derived chordoma xenograft and characterization of the phosphoproteome in a recurrent chordoma.
      • Siu I.M.
      • Salmasi V.
      • Orr B.A.
      • Zhao Q.
      • Binder Z.A.
      • Tran C.
      • Ishii M.
      • Riggins G.J.
      • Hann C.L.
      • Gallia G.L.
      Establishment and characterization of a primary human chordoma xenograft model.
      The advantage of this model is that it harbors the initiating mutations required for tumor formation. Tumor formation in the mouse was initially slow, taking longer than 3 months. However, with higher starting cell numbers, tumors now are palpable as early as 15 days after implantation. The slow growth of the tumor followed by a period of rapid expansion and subsequent death of the mice in all generations follows the same clinical course as observed in humans. The xenografts have the same cytology as the original tumor, namely physaliferous cells with round nuclei in lobules surrounded by connective tissue stroma (Figure 5E). The connective tissue within the xenografts is derived from the host mouse cells. Thus, physaliferous cells are able to direct the formation of supportive stroma. This is an important finding because this model would allow investigation of changes in the tumor microenvironment as a way of suppressing tumor growth. The xenograft tumor cells invade bone as is seen with human chordoma. The preoperative imaging of the parental tumor showed destruction of the clivus in an endophytic (type 1) pattern (Figure 1). Clival chordomas with this pattern of growth recently have been associated with differential expression of cell motility and invasion proteins compared with exophytic (type 2) clival chordomas.
      • Wu Z.
      • Wang L.
      • Guo Z.
      • Wang K.
      • Zhang Y.
      • Tian K.
      • Zhang J.
      • Sun W.
      • Yu C.
      Experimental study on differences in clivus chordoma bone invasion: an iTRAQ-based quantitative proteomic analysis.
      The present xenograft clival chordoma model could be used to further study the cell interactions that result in the observed tumor cytology, signaling processes in the tumor microenvironment, and mechanisms of bony invasion. Studying tumor in its native environment avoids alterations in matrix protein expression and cellular behavior that may be seen in 2-dimensional cell culture systems.
      • Reig G.
      • Pulgar E.
      • Concha M.L.
      Cell migration: from tissue culture to embryos.
      The established tumor xenograft, which had not been radiated previously, also will be useful to assess genomic and proteomic changes that occur upon radiation treatment and to study the effect of radiation sensitizers.
      At a genetic level the initial xenograft (PDX0) and its third generation (PDX3) show copy number changes compared with the parental tumor, indicating instability of the genome in this clival chordoma model. An associated increase in tumor growth was noticed over the PDX generations with an increased Ki67 index in PDX4 compared with the parental tumor. The preservation of copy number gains of a region of chromosome 6 involving the T gene locus in the clival chordoma PDX model reported here and in a sacral chordoma cell line (MUG-Chor1)
      • Owen J.H.
      • Komarck C.M.
      • Wang A.C.
      • Abuzeid W.M.
      • Keep R.F.
      • McKean E.L.
      • Sullivan S.
      • Fan X.
      • Prince M.E.P.
      UM-Chor1: establishment and characterization of the first validated clival chordoma cell line.
      is striking because it was preserved through three generations in our model and through multiple passages in the MUG-Chor1 cell line. By using RNAscope analysis a propensity to acquire increased T gene copy number from the PDX0 to the PDX3 generation was identified. These observations support the importance of T gene expression in chordoma tumor maintenance. Brachyury expression has been shown to be important in the maintenance of chordoma cell lines.
      • Hsu W.
      • Mohyeldin A.
      • Shah S.R.
      • ap Rhys C.M.
      • Johnson L.F.
      • Sedora-Roman N.I.
      • Kosztowski T.A.
      • Awad O.A.
      • McCarthy E.F.
      • Loeb D.M.
      • Wolinsky J.P.
      • Gokaslan Z.L.
      • Quinones-Hinojosa A.
      Generation of chordoma cell line JHC7 and the identification of Brachyury as a novel molecular target.
      • Shah S.R.
      • David J.M.
      • Tippens N.D.
      • Mohyeldin A.
      • Martinez-Gutierrez J.C.
      • Ganaha S.
      • Schiapparelli P.
      • Hamilton D.H.
      • Palena C.
      • Levchenko A.
      • Quinones-Hinojosa A.
      Brachyury-YAP regulatory axis drives stemness and growth in cancer.
      In addition, the increased T gene copy number may contribute to the enhanced proliferation of later generations in the xenograft model because increased copy number has been associated with phosphatidylinositol 3-kinase/Akt pathway activation.
      • Otani R.
      • Mukasa A.
      • Shin M.
      • Omata M.
      • Takayanagi S.
      • Tanaka S.
      • Ueki K.
      • Saito N.
      Brachyury gene copy number gain and activation of the PI3K/Akt pathway: association with upregulation of oncogenic Brachyury expression in skull base chordoma.
      The SNP analysis performed on the patient-derived chordoma xenograft described here showed similar patterns of copy number gains and losses as that observed in SNP data from both chordoma cell lines and tumors. Previously, heterozygous loss of CDKN2A was reported in 58% and homozygous loss was reported in 12% of a set of spinal and sacral chordomas.
      • Hallor K.H.
      • Staaf J.
      • Jonsson G.
      • Heidenblad M.
      • Vult von Steyern F.
      • Bauer H.C.
      • Ijszenga M.
      • Hogendoorn P.C.
      • Mandahl N.
      • Szuhai K.
      • Mertens F.
      Frequent deletion of the CDKN2A locus in chordoma: analysis of chromosomal imbalances using array comparative genomic hybridisation.
      In this model, it was shown that CDKN2A is lost heterozygously, but as part of widespread heterozygous chromosome 9 loss. The 9p deletions, including the CDKN2A locus, occurred in 4 of 18 (22%) clival chordomas.
      • Diaz R.J.
      • Guduk M.
      • Romagnuolo R.
      • Smith C.A.
      • Northcott P.
      • Shih D.
      • Berisha F.
      • Flanagan A.
      • Munoz D.G.
      • Cusimano M.D.
      • Pamir M.N.
      • Rutka J.T.
      High-resolution whole-genome analysis of skull base chordomas implicates FHIT loss in chordoma pathogenesis.
      Similarly, heterozygous loss of FHIT is seen as part of the loss of the short arm of chromosome 3. FHIT encodes a tumor-suppressor protein that is lost in a number of cancers.
      • Virgilio L.
      • Shuster M.
      • Gollin S.M.
      • Veronese M.L.
      • Ohta M.
      • Huebner K.
      • Croce C.M.
      FHIT gene alterations in head and neck squamous cell carcinomas.
      • Baffa R.
      • Veronese M.L.
      • Santoro R.
      • Mandes B.
      • Palazzo J.P.
      • Rugge M.
      • Santoro E.
      • Croce C.M.
      • Huebner K.
      Loss of FHIT expression in gastric carcinoma.
      • Hendricks D.T.
      • Taylor R.
      • Reed M.
      • Birrer M.J.
      FHIT gene expression in human ovarian, endometrial, and cervical cancer cell lines.
      • Sozzi G.
      • Veronese M.L.
      • Negrini M.
      • Baffa R.
      • Cotticelli M.G.
      • Inoue H.
      • Tornielli S.
      • Pilotti S.
      • De Gregorio L.
      • Pastorino U.
      • Pierotti M.A.
      • Ohta M.
      • Huebner K.
      • Croce C.M.
      The FHIT gene 3p14.2 is abnormal in lung cancer.
      • Ishii H.
      • Vecchione A.
      • Furukawa Y.
      • Sutheesophon K.
      • Han S.Y.
      • Druck T.
      • Kuroki T.
      • Trapasso F.
      • Nishimura M.
      • Saito Y.
      • Ozawa K.
      • Croce C.M.
      • Huebner K.
      • Furukawa Y.
      Expression of FRA16D/WWOX and FRA3B/FHIT genes in hematopoietic malignancies.
      FHIT protein expression was reduced or absent in 67% of clival chordomas and in 98% of sacral chordomas.
      • Diaz R.J.
      • Guduk M.
      • Romagnuolo R.
      • Smith C.A.
      • Northcott P.
      • Shih D.
      • Berisha F.
      • Flanagan A.
      • Munoz D.G.
      • Cusimano M.D.
      • Pamir M.N.
      • Rutka J.T.
      High-resolution whole-genome analysis of skull base chordomas implicates FHIT loss in chordoma pathogenesis.
      Single allelic loss of FHIT was observed in lung precancerous lesions resulting from smoke exposure.
      • D'Agostini F.
      • Izzotti A.
      • Balansky R.
      • Zanesi N.
      • Croce C.M.
      • De Flora S.
      Early loss of Fhit in the respiratory tract of rodents exposed to environmental cigarette smoke.
      In mice lacking one allele of the FHIT gene, 78% developed forestomach tumors after a single dose of the carcinogen N-nitrosylbenzylamine and the spontaneous tumor formation occurred at an average of 0.94 tumors per mouse after 21 months.
      • Zanesi N.
      • Fidanza V.
      • Fong L.Y.
      • Mancini R.
      • Druck T.
      • Valtieri M.
      • Rudiger T.
      • McCue P.A.
      • Croce C.M.
      • Huebner K.
      The tumor spectrum in FHIT-deficient mice.
      Thus, a single-hit event such as mono-allelic FHIT loss may be sufficient to trigger chordomagenesis. Despite expression of FHIT protein as observed in the established xenograft, protein function may be altered because altered transcripts have been reported in the FHIT gene in cervix cancer.
      • Hendricks D.T.
      • Taylor R.
      • Reed M.
      • Birrer M.J.
      FHIT gene expression in human ovarian, endometrial, and cervical cancer cell lines.
      The marked expression of FHIT in tumor cells suggests an important role of this protein in the cell biology of these cells, similar to expression seen in renal collecting tubule cells. It is interesting to note that both notochord cells, as well as the embryonic cells that give rise to intermediate mesoderm and subsequently renal collecting tubules, both showed expression of Brachyury.
      • Chapman D.L.
      • Agulnik I.
      • Hancock S.
      • Silver L.M.
      • Papaioannou V.E.
      Tbx6, a mouse T-Box gene implicated in paraxial mesoderm formation at gastrulation.
      Alterations of the mRNA transcript from the FHIT allele were not explored in this study; thus, the question of whether functional FHIT protein is present in this model remains unanswered at this time.
      Large- and small-scale genomic alterations including chromosomal imbalance and increasing copy number gains/losses through the whole genome across the tumor generations from parental to the third generation were observed. This study is the first to show chromosomal genomic instability in chordoma as the xenograft tumor grows through multiple generations. It is uncertain if the accumulation of copy number changes is an inherent characteristic of all chordoma cells or the result of selection bias toward tumor cells that are more likely to adapt to a nonhuman extracellular environment owing to loss of DNA repair mechanisms. However, the finding that at least a subpopulation of chordoma cells show accumulation of copy number changes may explain the resistance to radiation and chemotherapy shown by this tumor. It is known that genomic instability imparts resistance to chemotherapy in other tumors such as melanoma and colorectal cancer.
      • Brkic G.
      • Gopas J.
      • Tanic N.
      • Dedovic-Tanic N.
      • Benharroch D.
      • Finkelstein-Jaworowsky E.
      • Kedar I.
      • Dimitrijevic B.
      Genomic instability in drug-resistant human melanoma cell lines detected by Alu-I-arbitrary-primed PCR.
      • Lee A.J.
      • Endesfelder D.
      • Rowan A.J.
      • Walther A.
      • Birkbak N.J.
      • Futreal P.A.
      • Downward J.
      • Szallasi Z.
      • Tomlinson I.P.
      • Howell M.
      • Kschischo M.
      • Swanton C.
      Chromosomal instability confers intrinsic multidrug resistance.

      Conclusions

      We have established a patient-derived xenograft in immunocompromised mice from a previously nonirradiated clival chordoma. The tumor can be passaged and reproducibly generates a mass with cytology and immunohistology consistent with chordoma. Genomic gains/losses accumulate over generations, suggesting that chromosomal instability may play an important biological role in chordoma progression and chemotherapeutic resistance.

      Acknowledgment

      We thank the Chordoma Foundation for providing liaison support for independent immunohistochemical characterization of the xenograft.

      Supplemental Data

      • Supplemental Figure S2

        Heterozygous loss of CDKN2A is maintained from parental to PDX0 and PDX3 tumor. Chr9, chromosome 9; hg, human genome; MTAP, methylthioadenosine phosphorylase.

      • Supplemental Figure S3

        Heterozygous loss of FHIT as part of a large 3p isochromosome deletion is maintained from parental to PDX0 and PDX3 tumor. Chr3, chromosome 3; hg, human genome; PTPRG, protein tyrosine phosphatase receptor type G.

      • Supplemental Figure S4

        Chromosome 19 (Chr19) single-nucleotide polymorphism map shows acquisition of large isochromosome deletions (blue) and copy number gain (red) in small regions near the centromere in the transition from parental tumor to the first generation.

      • Supplemental Figure S5

        Chromosome 21 (Chr21) single-nucleotide polymorphism map shows acquisition of isochromosome deletions (blue) in the q arm and loss of large copy gain regions (red) in the transition from PDX0 to PDX3.

      • Supplemental Figure S6

        Comparison of chordoma samples across multiple data sets. Top: Data from this study showing segmentation across the entire genome of clival tumor and two PDX generations of 0 and 3. Middle: Data from a study looking at copy number of the original sacral tumor and the genomic stability in the derived cell line over several passages. Bottom: Patient samples analyzed using a chromosome tiling array showing recurrent copy number alterations in chordoma tumor samples. Cases 1, 2, 3, 4, 8, 10, and 11 are sacral and cases 5, 6, 7, and 9 are clival.

      References

        • Herrmann B.G.
        Expression pattern of the Brachyury gene in whole-mount TWis/TWis mutant embryos.
        Development. 1991; 113: 913-917
        • Wilkinson D.G.
        • Bhatt S.
        • Herrmann B.G.
        Expression pattern of the mouse T gene and its role in mesoderm formation.
        Nature. 1990; 343: 657-659
        • Salisbury J.R.
        • Deverell M.H.
        • Cookson M.J.
        • Whimster W.F.
        Three-dimensional reconstruction of human embryonic notochords: clue to the pathogenesis of chordoma.
        J Pathol. 1993; 171: 59-62
        • Shen J.
        • Li C.D.
        • Yang H.L.
        • Lu J.
        • Zou T.M.
        • Wang D.L.
        • Deng M.
        Classic chordoma coexisting with benign notochordal cell rest demonstrating different immunohistological expression patterns of Brachyury and galectin-3.
        J Clin Neurosci. 2011; 18: 96-99
        • Yamaguchi T.
        • Suzuki S.
        • Ishiiwa H.
        • Ueda Y.
        Intraosseous benign notochordal cell tumours: overlooked precursors of classic chordomas?.
        Histopathology. 2004; 44: 597-602
        • McMaster M.L.
        • Goldstein A.M.
        • Bromley C.M.
        • Ishibe N.
        • Parry D.M.
        Chordoma: incidence and survival patterns in the United States, 1973-1995.
        Cancer Causes Control. 2001; 12: 1-11
        • Diaz R.J.
        • Maggacis N.
        • Zhang S.
        • Cusimano M.D.
        Determinants of quality of life in patients with skull base chordoma.
        J Neurosurg. 2014; 120: 528-537
        • Chambers K.J.
        • Lin D.T.
        • Meier J.
        • Remenschneider A.
        • Herr M.
        • Gray S.T.
        Incidence and survival patterns of cranial chordoma in the United States.
        Laryngoscope. 2014; 124: 1097-1102
        • Diaz R.J.
        • Cusimano M.D.
        The biological basis for modern treatment of chordoma.
        J Neurooncol. 2011; 104: 411-422
        • Davies J.M.
        • Robinson A.E.
        • Cowdrey C.
        • Mummaneni P.V.
        • Ducker G.S.
        • Shokat K.M.
        • Bollen A.
        • Hann B.
        • Phillips J.J.
        Generation of a patient-derived chordoma xenograft and characterization of the phosphoproteome in a recurrent chordoma.
        J Neurosurg. 2014; 120: 331-336
        • Siu I.M.
        • Salmasi V.
        • Orr B.A.
        • Zhao Q.
        • Binder Z.A.
        • Tran C.
        • Ishii M.
        • Riggins G.J.
        • Hann C.L.
        • Gallia G.L.
        Establishment and characterization of a primary human chordoma xenograft model.
        J Neurosurg. 2012; 116: 801-809
        • Owen J.H.
        • Komarck C.M.
        • Wang A.C.
        • Abuzeid W.M.
        • Keep R.F.
        • McKean E.L.
        • Sullivan S.
        • Fan X.
        • Prince M.E.P.
        UM-Chor1: establishment and characterization of the first validated clival chordoma cell line.
        J Neurosurg. 2018; 128: 701-709
        • Bozzi F.
        • Manenti G.
        • Conca E.
        • Stacchiotti S.
        • Messina A.
        • Dagrada G.
        • Gronchi A.
        • Panizza P.
        • Pierotti M.A.
        • Tamborini E.
        • Pilotti S.
        Development of transplantable human chordoma xenograft for preclinical assessment of novel therapeutic strategies.
        Neuro Oncol. 2014; 16: 72-80
        • Trucco M.M.
        • Awad O.
        • Wilky B.A.
        • Goldstein S.D.
        • Huang R.
        • Walker R.L.
        • Shah P.
        • Katuri V.
        • Gul N.
        • Zhu Y.J.
        • McCarthy E.F.
        • Paz-Priel I.
        • Meltzer P.S.
        • Austin C.P.
        • Xia M.
        • Loeb D.M.
        A novel chordoma xenograft allows in vivo drug testing and reveals the importance of NF-kappaB signaling in chordoma biology.
        PLoS One. 2013; 8: e79950
        • Karikari I.O.
        • Gilchrist C.L.
        • Jing L.
        • Alcorta D.A.
        • Chen J.
        • Richardson W.J.
        • Gabr M.A.
        • Bell R.D.
        • Kelley M.J.
        • Bagley C.A.
        • Setton L.A.
        Molecular characterization of chordoma xenografts generated from a novel primary chordoma cell source and two chordoma cell lines.
        J Neurosurg Spine. 2014; 21: 386-393
        • Hsu W.
        • Mohyeldin A.
        • Shah S.R.
        • ap Rhys C.M.
        • Johnson L.F.
        • Sedora-Roman N.I.
        • Kosztowski T.A.
        • Awad O.A.
        • McCarthy E.F.
        • Loeb D.M.
        • Wolinsky J.P.
        • Gokaslan Z.L.
        • Quinones-Hinojosa A.
        Generation of chordoma cell line JHC7 and the identification of Brachyury as a novel molecular target.
        J Neurosurg. 2011; 115: 760-769
        • Presneau N.
        • Shalaby A.
        • Ye H.
        • Pillay N.
        • Halai D.
        • Idowu B.
        • Tirabosco R.
        • Whitwell D.
        • Jacques T.S.
        • Kindblom L.G.
        • Bruderlein S.
        • Moller P.
        • Leithner A.
        • Liegl B.
        • Amary F.M.
        • Athanasou N.N.
        • Hogendoorn P.C.
        • Mertens F.
        • Szuhai K.
        • Flanagan A.M.
        Role of the transcription factor T (brachyury) in the pathogenesis of sporadic chordoma: a genetic and functional-based study.
        J Pathol. 2011; 223: 327-335
        • Liu X.
        • Nielsen G.P.
        • Rosenberg A.E.
        • Waterman P.R.
        • Yang W.
        • Choy E.
        • Sassi S.
        • Yang S.
        • Harmon D.C.
        • Yang C.
        • Schwab J.H.
        • Kobayashi E.
        • Mankin H.J.
        • Xavier R.
        • Weissleder R.
        • Duan Z.
        • Hornicek F.J.
        Establishment and characterization of a novel chordoma cell line: CH22.
        J Orthop Res. 2012; 30: 1666-1673
        • DeComas A.M.
        • Penfornis P.
        • Harris M.R.
        • Meyer M.S.
        • Pochampally R.R.
        Derivation and characterization of an extra-axial chordoma cell line (EACH-1) from a scapular tumor.
        J Bone Joint Surg Am. 2010; 92: 1231-1240
        • Diaz R.J.
        • Guduk M.
        • Romagnuolo R.
        • Smith C.A.
        • Northcott P.
        • Shih D.
        • Berisha F.
        • Flanagan A.
        • Munoz D.G.
        • Cusimano M.D.
        • Pamir M.N.
        • Rutka J.T.
        High-resolution whole-genome analysis of skull base chordomas implicates FHIT loss in chordoma pathogenesis.
        Neoplasia. 2012; 14: 788-798
        • Bissig-Choisat B.
        • Wang L.
        • Legras X.
        • Saha P.K.
        • Chen L.
        • Bell P.
        • Pankowicz F.P.
        • Hill M.C.
        • Barzi M.
        • Kettlun Leyton C.
        • Leung H.C.
        • Kruse R.L.
        • Himes R.W.
        • Goss J.A.
        • Wilson J.M.
        • Chan L.
        • Lagor W.R.
        • Bissig K.D.
        Development and rescue of human familial hypercholesterolaemia in a xenograft mouse model.
        Nat Commun. 2015; 6: 7339
        • Wu Z.
        • Wang L.
        • Guo Z.
        • Wang K.
        • Zhang Y.
        • Tian K.
        • Zhang J.
        • Sun W.
        • Yu C.
        Experimental study on differences in clivus chordoma bone invasion: an iTRAQ-based quantitative proteomic analysis.
        PLoS One. 2015; 10: e0119523
        • Reig G.
        • Pulgar E.
        • Concha M.L.
        Cell migration: from tissue culture to embryos.
        Development. 2014; 141: 1999-2013
        • Shah S.R.
        • David J.M.
        • Tippens N.D.
        • Mohyeldin A.
        • Martinez-Gutierrez J.C.
        • Ganaha S.
        • Schiapparelli P.
        • Hamilton D.H.
        • Palena C.
        • Levchenko A.
        • Quinones-Hinojosa A.
        Brachyury-YAP regulatory axis drives stemness and growth in cancer.
        Cell Rep. 2017; 21: 495-507
        • Otani R.
        • Mukasa A.
        • Shin M.
        • Omata M.
        • Takayanagi S.
        • Tanaka S.
        • Ueki K.
        • Saito N.
        Brachyury gene copy number gain and activation of the PI3K/Akt pathway: association with upregulation of oncogenic Brachyury expression in skull base chordoma.
        J Neurosurg. 2018; 128: 1428-1437
        • Hallor K.H.
        • Staaf J.
        • Jonsson G.
        • Heidenblad M.
        • Vult von Steyern F.
        • Bauer H.C.
        • Ijszenga M.
        • Hogendoorn P.C.
        • Mandahl N.
        • Szuhai K.
        • Mertens F.
        Frequent deletion of the CDKN2A locus in chordoma: analysis of chromosomal imbalances using array comparative genomic hybridisation.
        Br J Cancer. 2008; 98: 434-442
        • Virgilio L.
        • Shuster M.
        • Gollin S.M.
        • Veronese M.L.
        • Ohta M.
        • Huebner K.
        • Croce C.M.
        FHIT gene alterations in head and neck squamous cell carcinomas.
        Proc Natl Acad Sci U S A. 1996; 93: 9770-9775
        • Baffa R.
        • Veronese M.L.
        • Santoro R.
        • Mandes B.
        • Palazzo J.P.
        • Rugge M.
        • Santoro E.
        • Croce C.M.
        • Huebner K.
        Loss of FHIT expression in gastric carcinoma.
        Cancer Res. 1998; 58: 4708-4714
        • Hendricks D.T.
        • Taylor R.
        • Reed M.
        • Birrer M.J.
        FHIT gene expression in human ovarian, endometrial, and cervical cancer cell lines.
        Cancer Res. 1997; 57: 2112-2115
        • Sozzi G.
        • Veronese M.L.
        • Negrini M.
        • Baffa R.
        • Cotticelli M.G.
        • Inoue H.
        • Tornielli S.
        • Pilotti S.
        • De Gregorio L.
        • Pastorino U.
        • Pierotti M.A.
        • Ohta M.
        • Huebner K.
        • Croce C.M.
        The FHIT gene 3p14.2 is abnormal in lung cancer.
        Cell. 1996; 85: 17-26
        • Ishii H.
        • Vecchione A.
        • Furukawa Y.
        • Sutheesophon K.
        • Han S.Y.
        • Druck T.
        • Kuroki T.
        • Trapasso F.
        • Nishimura M.
        • Saito Y.
        • Ozawa K.
        • Croce C.M.
        • Huebner K.
        • Furukawa Y.
        Expression of FRA16D/WWOX and FRA3B/FHIT genes in hematopoietic malignancies.
        Mol Cancer Res. 2003; 1: 940-947
        • D'Agostini F.
        • Izzotti A.
        • Balansky R.
        • Zanesi N.
        • Croce C.M.
        • De Flora S.
        Early loss of Fhit in the respiratory tract of rodents exposed to environmental cigarette smoke.
        Cancer Res. 2006; 66: 3936-3941
        • Zanesi N.
        • Fidanza V.
        • Fong L.Y.
        • Mancini R.
        • Druck T.
        • Valtieri M.
        • Rudiger T.
        • McCue P.A.
        • Croce C.M.
        • Huebner K.
        The tumor spectrum in FHIT-deficient mice.
        Proc Natl Acad Sci U S A. 2001; 98: 10250-10255
        • Chapman D.L.
        • Agulnik I.
        • Hancock S.
        • Silver L.M.
        • Papaioannou V.E.
        Tbx6, a mouse T-Box gene implicated in paraxial mesoderm formation at gastrulation.
        Dev Biol. 1996; 180: 534-542
        • Brkic G.
        • Gopas J.
        • Tanic N.
        • Dedovic-Tanic N.
        • Benharroch D.
        • Finkelstein-Jaworowsky E.
        • Kedar I.
        • Dimitrijevic B.
        Genomic instability in drug-resistant human melanoma cell lines detected by Alu-I-arbitrary-primed PCR.
        Anticancer Res. 2003; 23: 2601-2608
        • Lee A.J.
        • Endesfelder D.
        • Rowan A.J.
        • Walther A.
        • Birkbak N.J.
        • Futreal P.A.
        • Downward J.
        • Szallasi Z.
        • Tomlinson I.P.
        • Howell M.
        • Kschischo M.
        • Swanton C.
        Chromosomal instability confers intrinsic multidrug resistance.
        Cancer Res. 2011; 71: 1858-1870