Abstract
Aberrant activation of the Hedgehog signaling pathway has been linked to the formation of numerous cancer types, including the myogenic soft tissue sarcoma, embryonal rhabdomyosarcoma (eRMS). Here, we report PCG2, a novel mouse model in which human GLI2A, a constitutive activator of Hedgehog signaling, induced undifferentiated sarcomas that were phenotypically divergent from eRMS. Rather, sarcomas arising in PCG2 mice featured some characteristics that were reminiscent of Ewing sarcoma. Even though it is widely understood that Ewing sarcoma formation is driven by EWS-ETS gene fusions, a genetically defined mouse model is not well-established. While EWS-ETS gene fusions were not present in PCG2 sarcomas, precluding their designation as Ewing sarcoma, we did find that GLI2A induced expression of known EWS-ETS gene targets essential to Ewing pathogenesis, most notably, Nkx2.2. Moreover, we found that naïve mesenchymal progenitors originate tumors in PCG2 mice. Altogether, our work provides a novel genetic mouse model, which directly connects oncogenic Hedgehog activity to the etiology of undifferentiated soft tissue sarcomas for the first time.
Implications: The finding that activation of Gli2 transcription factor is sufficient to induce Ewing-like sarcomas provides a direct transformative role of the Hedgehog signaling pathway in undifferentiated soft tissue sarcoma.
Introduction
Small round-cell sarcomas are a category of phenotypically diverse, often rare, but clinically aggressive cancers. The relative rareness of these soft tissue or bone-localized tumors versus other solid tumor types (lung, GI, skin, etc) and limited access to prospectively isolated clinical specimens needed for laboratory investigation has caused them to be less well studied by comparison (1). This deficit in accessibility may have hindered progress in testing potential diagnostic markers and development of new therapeutic interventions for cancers that can be difficult to manage clinically especially when not accurately diagnosed. Fortunately, for several sarcoma entities genetic mouse models exist as a tool that can mitigate the problem of specimen access. Models for alveolar rhabdomyosarcoma (aRMS), synovial sarcoma, clear cell sarcoma, and numerous embryonal rhabdomyosarcoma (eRMS) models have been well described, and often have succeeded owing to incorporation of a defining genetic lesion into their design.
For example, several fusion-positive sarcoma models have utilized clinically relevant gene fusions, such as PAX43-FKHR in aRMS, SYT-SSX in synovial sarcoma, and EWS-ATF1 in clear cell sarcoma (2–4). In addition, an array of models has been created for the fusion-negative sarcoma, eRMS, utilizing numerous genetic mechanisms, but haven often incorporated induction of the developmental signaling pathway (5, 6), Hedgehog, which alone can be sufficient to induce formation of eRMS in mice (7, 8). As more data accumulates related to gene fusions identified in other sarcoma variants, such as CIC-DUX4, FUS-NFATc2, and BCOR-CCNB3 for example, it will be important to test their respective activities in mouse model design. Although, one test case for this modeling approach that has consistently underwhelmed, despite novel and clever design approaches from many groups worldwide, is that of Ewing sarcomas (9), which are clinically defined by the presence of EWS and ETS gene fusions, most notably EWS-FLI1 (10).
In this article, we endeavor to describe a novel genetic mouse model called PCG2, in which induction of Hedgehog signaling was achieved using an activated human GLI2 (GLI2A). Much to our surprise, tumors formed by these mice diverged phenotypically from eRMS, and instead displayed pathologic, ultrastructural, immunophenotypic, and gene expression patterns reported by others in human Ewing sarcoma. Although key differences were noted as well, such as the absence of a relevant EWS-ETS gene fusion—the most critical diagnostic criterion for Ewing sarcoma. We determined that PCG2 sarcomas emerged from Pcp2-cre–expressing multipotent mesenchymal progenitor cells, which were not transformed by ectopic activation of EWS-FLI1 in separate experiments. We also observed a dramatic sensitivity to Gli1 gene dosage, as partial loss diminished the PCG2 phenotype. Altogether, we found that GLI2A is sufficient to induce undifferentiated soft tissue sarcomas in a new genetically defined mouse model, adding to the overall understanding of the transformative role for Hedgehog in sarcoma pathogenesis.
Materials and Methods
Animals
All experiments were performed using embryos from timed pregnancies or young neonatal and adult animals (ages P3–1.5 years) according to the NIH and VUMC Division of Animal Care. Mice with the following alleles, of either sex, were used for this study and were maintained on a mixed genetic background. Animals were either obtained from The Jackson Laboratory or the originating laboratory. Ai9 [Gt(ROSA)26Sortm9(CAG-tdTomato)Hze], CLEG2, EWSEWS-FLI1, Gli1nlacZ [Gli1tmAlj], Myf5Cre [Myf5tm1(cre)Mrc], Pax7-creER [Pax7tm2.1(cre/ERT2)Fan], Pcp2-cre [Tg(Pcp2-cre)1Amc], PtchlacZ [Ptch1tm1Mps], R26ReYFP [Gt(ROSA)26Sortm1(EYFP)Cos], SmoM2 [Gt(ROSA)26Sortm1(Smo/EYFP)Amc], and SpiB−/−.
Tamoxifen (Sigma) was administered to Pax7-creER; CLEG2 neonates at P3 or P7, and young adults at P21. A stock solution of 2 mg/mL was prepared in corn oil (Sigma), and a dose of 100 μg or 600 μg was administered to neonates or young adults, respectively.
Tumor volume was determined using digital calipers (Electron Microscopy Sciences) to measure the largest and smallest diameter, from the earliest signs of tumor palpability until reaching a largest diameter measuring no larger than 3 cm.
Human tumor specimens
Deidentified human sarcoma blocks and slides were retrieved from pathology archives at Vanderbilt University Medical Center (Nashville, TN).
Cell lines
The PCG2M226 tumor cell line was generated in our laboratory in 2011 and authenticated using cell grafting, IHC, and RNA sequencing (RNAseq)-based comparisons with parent tumor tissue. PCG2M226 tumor cells were isolated from freshly resected tumor tissue from PCG2 animals under sterile conditions. Small samples (1–2 mm2) were transferred to 60-mm culture dishes and minced with microdissection scissors. Cells were allowed to adhere to plastic prior to media change (24–48 hours), and were in culture for 1.5–2 weeks until large colonies were observed. Colonies were dissociated with 0.05% trypsin, spun down, washed, and plated for passaging. Following two passages (7–10 days) near homogeneous cultures were obtained, which could then be passaged every 3 days. Cell lines were maintained in standard culture conditions (37°C, 5% CO2) in DMEM (Gibco) supplemented with 10% FBS and penicillin–streptomycin. PCG2M226 tumor cells have been kept as “low passage” (<6) or “high passage” (>6).
Grafting experiments
Sarcoma cells (either mouse or human) were injected unilaterally into gastrocnemius muscle of NOD-SCID mice under anesthesia. By pressing the muscle of the lower hind limb together with thumb and index finger below the knee joint, the gastrocnemius muscle becomes readily conspicuous after removal of fur. Approximately 2 × 106 cells resuspended in HBSS were injected using 1-mL syringes capped with 30 gauge needles.
Tissue processing, protein isolation, IHC, and Western blotting
Tissues were dissected and fixed in 4% paraformaldehyde for either 4–6 hours or O/N at 4°C, and were either processed for frozen embedding in OCT compound or processed for paraffin embedding. Frozen tissues were sectioned on a Leica cryostat at 10 μm, paraffin-embedded tissues were cut at 5 μm. Specimens containing portions of bone were decalcified in 0.4 mol/L EDTA at 4°C for 2 weeks prior to dehydration. Protein was isolated from fresh or snap-frozen tumor and tibialis anterior muscle using standard lysis buffer, and a BCA Assay (Thermo Fisher Scientific) was used for measuring concentration. IHC and immunocytochemistry were performed with standard protocols, 1 mmol/L EDTA pH 9.0 was used for antigen retrieval with tumor sections. Mouse on mouse blocking reagent (Vector Laboratories) was used with mouse primary antibodies.
Antibodies
The following primary antibodies were used to perform IHC on frozen and/or paraffin tissue sections: mouse α-β-Catenin (Vanderbilt Antibody and Protein Resource), rabbit α-Cd99 (Dietmar Vestweber, Max Planck Institute for Molecular Biomedicine, Röntgenstr, Münster, Germany), rabbit α-cleaved Caspase 3 (Cell Signaling Technology), mouse α-CyclinD1 (DSHB), rabbit α-Desmin (Thermo Fisher Scientific), mouse α-EZH2 (Cell Signaling Technology), rabbit α-Foxd3 (Dr. Patricia Labosky, NIH, Bethesda, MD), chicken α-GFP (Aaves), rabbit α-Gli1 (Cell Signaling Technology), rabbit α-γH2AX (BETHYL), rabbit α-Keratin (Sigma), rabbit α-Ki67 (NeoMarkers), mouse α-Laminin (Thermo Fisher Scientific), mouse α-Myc (DSHB), mouse α-Myc Tag 9B11 (Cell Signaling Technology), mouse α-Myogenin (DSHB), rabbit α-Myogenin (EPITOMICS), mouse α-Myf5 (DSHB), mouse α-MyoD (DAKO), mouse α-Nkx2.2 (DSHB), mouse α-Pax7 (DSHB), mouse α-Pax3 (DSHB), rabbit α-p-histone H3 (Millipore), rabbit α-PPARγ (Cell Signaling Technology), rabbit α-SpiB (Cell Signaling Technology), rabbit α-Tnc (Dr. Herald Erickson, Duke University Medical Center, Durham, NC), and mouse α-tubulin (DSHB). Species-specific HRP-conjugated secondary antibodies (Invitrogen) were used followed by incubation in DAB reaction (Invitrogen). Double-labeling fluorescence IHC was performed using species-specific, AlexaFluor secondary antibodies (Invitrogen) followed by counterstaining with To-pro3 iodide (Invitrogen).
Chromatin immunoprecipitation and quantitative PCR
Chromatin immunoprecipitation (ChIP) was carried out on PCG2M226 cells of approximately 3 × 106 cells per epitope. Chromatin from 1% formaldehyde fixed was fragmented to a range between 200 and 700 bases using a BioRuptor (Diagenode). Solubilized chromatin (5 μg) was immunoprecipitated with the following antibodies; mouse α-Myc Tag 9B11 (Cell Signaling Technology catalog no. 2267) and rabbit α-H3K4me3 (Cell Signaling Technology catalog no. 9751). The latter antibody was used for trouble shooting purposes. Antibody–chromatin complexes were pulled down with either protein G-Dynabeads or protein A-Dynabeads (Life Technologies), for mouse and rabbit primaries, respectively, then washed and eluted. Following crosslink reversal with proteinase K treatment overnight at 65°C, the immunoprecipitated DNA was extracted using the PCR Purification Kit from Qiagen.
qRT-PCR was performed using specific primers (for Ptch1, Gli1, and Nkx2.2) and SYBR Green PCR Master Mix in a CFX96 Thermocycler (Bio-Rad). Percent input was used to calculate enrichment, and species-specific IgG antibodies were used as controls. Primers were designed on the basis of regulatory regions described previously (11), sequences are as follows:
Ptch1 forward TTCTTTGGAGCTCAATTTCC; Ptch1 reverse GTTTTCCCCGATTTTAAGGT
Gli1 forward AAGCCAGATGTGATGTAGCA; Gli1 reverse TTATGGAAGGACAGACAGCA
Nkx2.2 forward ATCCAAGCTCTGAGTTGAGG; Nkx2.2 reverse AAGTGCCAAACAGGTCTCTC
Gene expression analysis
RNAseq data analysis.
We first performed quality control on raw data to identify potential outliers before doing any advanced analysis by applying tool FastQC. RNA data alignment was performed by TopHat and gene quantification was done by Cufflinks. Unsupervised cluster analysis based on Pearson correlation coefficient on all genes shows clear separation between muscle and tumor samples. RPKM (reads per kilobase per million reads) based approaches with Cuffdiff were used to detect differentially expressed genes. FDR < 0.05 was used for multiple test correction.
Gene set enrichment analysis.
Gene set enrichment analysis (GSEA) was performed using software and datasets made available by the Broad Institute (http://software.broadinstitute.org/gsea/index.jsp), as well as several published modules as indicated. A rank was generated of the PCG2 transcriptome (average of three tumors) for GSEA, which included 20,461 gene values. Gene sets that were significant at FDR < 25% and significantly enriched at nominal P < 1% or <5% were considered significantly enriched.
Principal components analysis and hierarchal clustering analysis.
RNAseq data for PCG2 and PCG2M226 cells were generated as described at VANTAGE. We obtained normalized RNAseq data for Ewing sarcoma (n = 66) from the NCI (12), and data for other human cancers, including synovial sarcoma (n = 10), myxofibrosarcoma (n = 25), undifferentiated pleomorphic sarcoma (n = 21), and diffuse large B-cell lymphoma (n = 48), were obtained from TCGA. For each dataset, we removed genes whose median expression is zero and log transformed the gene expression after filtering. We manually inspected the distribution of the expression data in all datasets and confirmed that all have similar distribution. After centering and quantile normalization of the data we ran principal component analysis (PCA) and plotted all samples based on their first three principal components.
RNA isolation, reverse transcription, and quantitative PCR.
Total RNA was isolated from paraffin-imbedded tissue using the RecoverAll Total Nucleic Acid Isolation Kit (ambion by Life Technologies). cDNAs were synthesized with 300 ng total RNA input for all samples using a High-Capacity cDNA Reverse Transcription Kit (Applied Biosciences). qRT-PCR was performed using specific primers for Gli1 and Gapdh and SYBR Green PCR Master Mix in a CFX96 Thermocycler (Bio-Rad), and fold change was calculated using the ΔCt method. Sequences for Gli1 primers are as follows:
Gli1 forward CCTGAGCCTGAGTCTGTGTATG; Gli1 reverse CGTGGATATGCTCACTGTTGAT.
Microscopy
Bright-field and confocal imaging.
Bright-field images were collected on an Olympus BX51 upright microscope and a Leica M165 FC stereoscope. Fluorescent images were acquired using a Leica TCS SP5 laser-scanning confocal.
Transmission electron microscopy.
Specimens were processed for transmission electron microscopy (TEM) and imaged in the Vanderbilt Cell Imaging Shared Resource-Research Electron Microscopy facility. TEM was performed on a Philips/FEI T-12.
Embedding.
Samples were fixed in 2.5% glutaraldehyde in 0.1 mol/L cacodylate buffer, pH 7.4 at room temperature for 1 hour then transferred to 4°C, overnight. The samples were washed in 0.1 mol/L cacodylate buffer, incubated for 1 hour in 1% osmium tetraoxide at room temperature, and then washed again with 0.1 mol/L cacodylate buffer. Subsequently, the samples were dehydrated through a graded ethanol series and then three exchanges of 100% ethanol, followed by two exchanges of pure propylene oxide (PO). Samples were then infiltrated with 25% Epon 812 resin and 75% PO for 30 minutes at room temperature. Next, they were infiltrated with Epon 812 resin and PO [1:1] for 1 hour at room temperature then overnight at room temperature. The samples are subsequently infiltrated with resin for 48 hours then allowed to polymerize at 60°C for 48 hours.
Sectioning and imaging.
500-nm to 1-μm–thick sections were collected using a Leica Ultracut microtome. Thick sections were contrast stained with 1% toluidine blue and imaged with a Nikon AZ100 microscope. Ultra-thin sections (70–80 nm) were cut and collected on 300-mesh copper grids and poststained with 2% uranyl acetate and then with Reynold lead citrate. Samples were subsequently imaged on the Philips/FEI Tecnai T12 electron microscope at various magnifications.
Cell counting
For counting, entire tumor sections were imaged on a Leica SCN400 slide scanner. Cell counts were obtained from all intact tissue for Ki67, cleaved caspase 3, and phospho-Histone H3, and were quantified as positive cells per mm2. Slides from at least 5 animals per genotype were used for counting.
Transcript detection
For RNA in situ hybridization, Gli1 cDNA was used as a template for synthesizing digoxygenin-labeled riboprobes.
Spectral karyotyping
Spectral karyotyping (SKY) was performed by Roswell Park Cancer Institute's SKY laboratory (Buffalo, NY) on metaphase PCG2M226 and PCG2M656 cells (https://www.roswellpark.edu/shared-resources/pathology-resource-network/services-and-fees/sky/fish).
RNA isolation and sequencing
Total RNA from PCG2 tumor, tibialis anterior muscle, and PCG2M226 cell lines, was isolated using a two-step protocol in which tissue/cells were homogenized in TRizol and RNA was collected by phenol–chloroform extraction and column elution with RNeasy Mini Kit (Qiagen). RNA was reverse transcribed to create a cDNA library by the VANTAGE core at VUMC prior to paired-end sequencing on Illumina Hiseq 2500.
Statistical analysis
Bar charts and survival curves were generated using Prism software (GraphPad). Clustering and PCA were performed with R.
Accession numbers
The data accompanying this article were deposited into the GEO database under accession number GSE89419.
Results
PCG2 mice form undifferentiated soft tissue sarcomas
Initially, we sought to activate Hedgehog signaling in mouse cerebellar Purkinje neurons, which produce, but do not respond to, the ligand, Sonic hedgehog (Shh). We bred Pcp2-cre (13) and CLEG2 transgenic mice (14) to conditionally express a constitutively activated human GLI2 (GLI2A; Fig. 1A and C). After cre-mediated excision of an upstream GFP-coding sequence and stop signal, Purkinje neurons expressing GLI2A appeared qualitatively normal (Supplementary Fig. S1A). However, to our surprise, these Pcp2-cre; CLEG2 (PCG2) animals developed soft tissue tumors (Fig. 1A) similar to human small round-cell sarcomas, displaying cohesive sheets of small round cells with high nuclear to cytoplasmic (n/c) ratio and geographic necrosis (Fig. 1A; Supplementary Fig. S1B). Ninety-eight percent of male and female mice evaluated (n > 50) formed fast growing, encapsulated tumors by 8 weeks of life on average that targeted multiple anatomic sites (P < 0.0001; Fig. 1D and E). Invasion of tumor cells into adjacent muscle and bone marrow was found in long bones of ipsilateral limbs, but was not apparent in other tissues, likely owing to rapid growth kinetics that necessitated timely sacrifice (Supplementary Fig. S1B and S1C). Intriguingly, mitotic figures were rare in PCG2 tumors (Supplementary Fig. S1D).
GLI2A induced undifferentiated sarcomas distinct from PtchlacZ eRMS. A, Hematoxylin and eosins (H&E) stains of PCG2 and PtchlaZ mouse tumors, and human Ewing sarcoma and eRMS. B, Immunostaining as indicated, PCG2 tumors lacked myogenic markers Myogenin MyoD, Desmin, and Pax7, which were present in PtchlaZ RMS. Both tumors broadly expressed proliferation marker Ki67. C, Schematic showing CLEG allele and breeding strategy with Pcp2-cre mice. D, Frequency of tumor foci at various anatomic sites and external view of tumor masses located on the right hip/thigh and ipsilateral shoulder/chest (white arrows). E, Kaplan–Meier survival plot for PCG2 and CLEG2 mice, tumor margin on right thigh muscle of a separate mouse. Scale bars, 50 μm.
Because oncogenic Hedgehog signaling has previously been linked to the formation of eRMS, we compared PCG2 tumors with those from PtchlacZ mice. Distinctive features seen in PtchlacZ eRMS (i.e., rhabdomyoblasts and entrapped myofibers) were absent from PCG2 tumors (Fig. 1A). Comparison of immunophenotypes consistently indicated that PCG2 tumors lacked myogenic differentiation, although Myogenin, MyoD, Desmin, and Pax7 were readily detected in PtchlacZ eRMS (Fig. 1B). Both tumors displayed a large fraction of cycling cells indicated by Ki67+ status (Fig. 1B). We measured the mitotic and apoptotic indices of PCG2 tumors using phospho-histone H3 and cleaved caspase 3, respectively, and found on average approximately 479 cells/mm2 were mitotic and approximately 337 cells/mm2 apoptotic.
We evaluated GLI2A expression using the Myc epitope tag encoded by CLEG2, which revealed widespread nuclear staining (Supplementary Fig. S1E). Persistent Hedgehog signaling was evident in PCG2 tumors, indicated by expression of target gene Gli1, which was not detected in adjacent myofibers (Supplementary Fig. S1F). Gli1 and Cyclin D1 protein levels were also high in PCG2 tumor lysates compared with tibialis anterior muscle, further confirming persistent Hedgehog pathway activity (Supplementary Fig. S1G). We crossed Pcp2-cre and SmoM2 mice, to evaluate whether the upstream acting oncogene could differentially influence tumor phenotype. Unfortunately, these mice succumbed at weaning with severe brain malformations not seen in PCG2 animals (Supplementary Fig. S1H). These findings indicate that GLI2A is sufficient to drive the formation of soft tissue sarcomas distinct from eRMS.
Evaluating the gene expression profile of PCG2 sarcomas
Next, RNAseq was conducted for comparative gene expression profiling of PCG2 sarcomas. Hierarchal clustering of the top 100 most differentially expressed genes between tumor and unaffected tibialis anterior muscle revealed distinct gene profiles (Fig. 2A; Supplementary Table S1). As expected, Hedgehog target genes (Gli1, Ptch1, and Hhip) were highly upregulated in tumors. Unexpectedly, we also detected induction of the proneural transcription factor and known Hedgehog target Nkx2.2 in PCG2 tumors (Fig. 2A’). Immunostaining confirmed positivity for Nkx2.2 protein in all tumors analyzed (n > 25), while PtchlacZ eRMS were consistently negative (Fig. 2B). Because Nkx2.2 is well documented as an essential target gene of EWS-FLI1 in Ewing sarcoma pathogenesis (15, 16), but is not frequently seen in other sarcomas (15), we chose to incorporate additional characteristic elements of Ewing sarcoma described by other groups elsewhere into our analyses.
Exploring the PCG2 tumor transcriptome. A, Hierarchal clustering analysis of RNAseq data shows top 100 most differentially expressed genes from tumor versus unaffected tibialis anterior muscle in PCG2 mice (n = 3). A’, RPKM values show induction of Hedgehog targets, Ewing sarcoma–related genes, and select developmental genes in PCG2 tumors versus tibialis anterior muscle. B, Immunophenotype indicates PCG2 tumors, but not PtchlaZ eRMS, robustly express Ewing sarcoma markers Nkx2.2 and Cd99, as well as methyltransferase EZH2. C, TEM on ultra-thin PCG2 tumor sections. D, GSEA shows enrichment of indicated gene sets. Scale bars, 50 μm.
In line with previous animal model studies (17, 18), we attempted immunostaining against Cd99, which revealed striking membranous localization in all PCG2 sarcomas tested, while staining in PtchlacZ eRMS was clearly diffuse (Fig. 2B). This result was intriguing at first, given findings by others and that conservation of the latter is limited at the sequence level in mice, but is reconciled in that functionality of the protein may be somewhat well conserved (19–21). The histone methyltransferase, enhancer of zest homolog 2, another known EWS-FLI1 target (22), was detected in PCG2 tumors, but not in PtchlacZ eRMS (Fig. 2B). TEM revealed neurosecretory granule-like particles, abundant cytoplasmic glycogen, a high n/c ratio, and close cellular contacts in PCG2 sarcomas, all of which have been observed previously in Ewing sarcoma, but not in other sarcomas to our knowledge (ref. 23; Fig. 2C; Supplementary Fig. S2A, and S2A’).
We utilized GSEA (24) to explore the activity of previously defined gene modules in the PCG2 transcriptome (containing 20,461 ranked genes). Significantly enriched gene sets included epithelial-to-mesenchymal transition (NES = 1.68; Nom P = 0.0; FDR q-val < 0.03; FWER P = 0.0) and G2–M checkpoint modules (NES = 2.5; Nom P = 0.0; FDR q-val = 0.0; FWER P = 0.0; not shown), consistent with the fast-growing nature of PCG2 tumors. Meanwhile, genes upregulated in other cancers, such as mesenchymal glioblastoma multiforme, alveolar RMS, and malignant fibrous histiocytoma, did not show significant enrichment (NES = −0.53, Nom P = 0.97, FDR q-val = 0.97, FWER P = 0.92; NES = −0.58, Nom P = 0.33, FDR q-val = 0.33, FWER P = 0.28; and NES = 1.34, Nom P = 0.09, FDR q-val = 0.09, FWER P = 0.03, respectively; refs. 25–27). However, 53 genes upregulated in mouse Hedgehog-driven medulloblastoma were significantly enriched in PCG2 tumors (NES = 1.96; Nom P = 0.0; FDR q-val = 0.0; FWER P = 0.00; ref. 28).
Other tumor gene modules with significant enrichment included 500 genes upregulated by EWS-FLI1 in Ewing sarcoma cells (NES = 1.97; Nom P = 0.0; FDR q-val = 0.0; FWER P = 0.0; ref. 29), and the EWS-FLI1–core module (NES = 1.38; Nom P = 0.0; FDR q-val = 0.0; FWER P = 0.0; Fig. 2D; Supplementary Fig. S2B; ref. 30). Of the 293 genes downregulated by EWS-FLI1 in the latter module, 219 were strongly downregulated in PCG2 tumors. The Ewing sarcoma GLI1 gene module, consisting of 86 genes upregulated by GLI1 (31), also showed significant enrichment (NES = 1.46; Nom P = 0.023, FDR q-val = 0.023; FWER P = 0.004; Fig. 2D). Keratin 17, whose function is necessary for transformation and cellular adhesion in Ewing sarcoma (31), was the leading gene. E2F targets (NES = 2.7; Nom P = 0.0, FDR q-val = 0.0; FWER P = 0.0; Supplementary Fig. S2B), which ETS fusions are also known to deregulate, were upregulated (32). The PCG2 transcriptome showed upregulation of additional Ewing sarcoma–associated genes, including Hoxd13 (33), and proinvasion gene Trip6 (34), while downregulation of Stag2, whose loss confers poor prognosis in Ewing sarcoma, was observed (Fig. 2A’; ref. 12). Many genes commonly activated by EWS-FLI1 were not upregulated in PCG2 tumors, like Prkcb, while others reversed trend, such as Gstm4, Parp1, and Cav1 (Supplementary Fig. S2C). Nr0b1 did not appear in the PCG2 dataset, and therefore could not be evaluated.
To test whether EWS-FLI1 could induce tumors in our system, and to have a true point of comparison for PCG2 mice, we crossed Pcp2-cre mice with an allele expressing EWS-FLI1 from the endogenous EWS gene locus (35). However, EWS-FLI1 did not induce tumors in Pcp2-cre; EWSEWS-FLI1 mice (n = 7), indicating that this oncogene alone is not sufficient to transform Pcp2-lineage cells. However, activity of the oncogene was apparent, as all Pcp2-cre; EWSEWS-FLI1 mice displayed ataxia and exopthalmia (not shown), consistent with cerebellar and retinal expression of Pcp2-cre, and apoptosis induced by EWS-FLI1 in certain tissues (35, 36). We also checked for expression of Fli1 and other ETS genes in the PCG2 transcriptome, but found that only SpiB was differentially upregulated (Supplementary Fig. S2D). Immunostaining confirmed SpiB expression, but ultimately tumors indistinguishable from PCG2 formed in PCG2; SpiB−/− animals (not shown).
GLI2A-mediated induction of Nkx2.2 expression, genomic instability, and a distinctive gene profile
We isolated primary cells with serial passaging of a dissociated PCG2 tumor (Supplementary Fig. S3A). Expression of mitotic markers and Nkx2.2 was robust in the resulting PCG2M226 cell line, which, when injected into gastrocnemius muscle of NOD/SCID recipient mice, propagated aggressive tumors within 10–14 days (Supplementary Fig. S3A and S3B). These tumors showed a robust Nkx2.2+, Cd99+ signature and Myc expression, similar to parent PCG2 sarcomas (Fig. 3A; Supplementary Fig. S3B). In an effort to understand Nkx2.2 induction, ChIP was performed with chromatin from PCG2M226 cells (Supplementary Fig. S3C). By targeting the GLI2A-MYC epitope tag, we found enrichment at previously described regulatory elements in Hedgehog target genes Gli1 and Ptch1, and importantly, in Nkx2.2 (Supplementary Fig. S3C). GLI2A occupancy in this Nkx2.2 regulatory element is consistent with studies of Hedgehog-dependent patterning in the neural tube (11), and provides a mechanism for Nkx2.2 induction, which appears specific to PCG2, given that Nkx2.2 is absent in all other reported Hedgehog-driven mouse tumors.
Induction of Nkx2.2 expression and a distinctive gene signature. A, Immunostaining of allografted PCG2M226 cells shows Cd99 and Nkx2.2 localization in resultant tumor. B, PCA of transcriptomes of indicated human tumors and PCG2 tumors/cells. C, A heatmap and dendrogram of the same human tumors and PCG2 tumors/cells (x-axis) is shown by hierarchal clustering analysis of 13,345 genes (y-axis). Scale bars, 50 μm.
Because many sarcomas are defined by chromosomal translocation, we tested whether such changes accompany tumorigenesis in PCG2 mice. PCG2 sarcomas showed greater positivity for double-stranded break marker γH2AX than PtchlacZ eRMS (Supplementary Fig. S3D), indicative of heightened genomic instability (37). SKY of PCG2M226 cells revealed amplifications, deletions, and a t(1;7) translocation (Supplementary Fig. S3E). These findings indicate that Hedgehog activation through GLI2A, but not Ptch1 loss-of-function, induces persistent genomic instability, which might explain the more aggressive features of PCG2 tumors and is consistent with previous studies of GLI2A (38).
We used RNAseq to profile the transcriptome of PCG2M226 cells, which demonstrated a 0.98 correlation coefficient with a resequenced parent tumor. PCA was used to compare transcriptomes of PCG2 tumors/PCG2M226 cells (average of 3; red/n = 2; royal blue), and five different small round-cell human tumor types with data from The Cancer Genome Atlas (TCGA). The latter group included synovial sarcoma (n = 10; purple), myxofibrosarcoma (n = 25; orange), undifferentiated pleomorphic sarcoma (UPS; n = 21; brown), diffuse large B-cell lymphoma (DLBCL; n = 48; green), and clinically and experimentally defined Ewing sarcoma (n = 61; light blue; ref. 12). The results showed dissimilarity between PCG2 tumors/PCG2M226 cells and DLBCL, myxofibrosarcoma, and UPS, but similarity with certain human Ewing sarcoma samples (Fig. 3C). Separation of PCG2 into three individual tumors (red), and Ewing sarcoma into extraosseous (n = 6; light blue) and osseous (n = 55; dark green) did not alter these trends (Supplementary Fig. S3F).
We chose to focus on the three Ewing sarcoma samples that fell nearest to the PCG2 domain. EWS112 was a soft tissue tumor of the neck harboring an EWS-FLI1 type I fusion and expression loss of STAG2. EWS102 was a tumor of the femur with a novel FUS-NFATc2 fusion, and EWS140 was a soft tissue tumor of the thigh with no detected gene fusion (Fig. 3B), but both were clinically diagnosed as Ewing sarcoma (12). In hierarchal clustering analysis with the above tumor groups, PCG2 clustered with the bulk of primary Ewing sarcoma samples when 13,345 gene expression values per sample were compared, indicated by the accompanying dendrogram (Fig. 3C).
Gli1 gene dosage influenced sarcoma phenotypic outcome
To localize ectopic Hedgehog signals in PCG2 tissues, we utilized heterozygous Gli1nlacZ knock-in reporter mice, which carry only one functional Gli1 allele. The resulting PCG2; Gli1nlacZ animals formed soft tissue tumors with histologic features that were divergent from PCG2 tumors, showing instead nested cells with peripherally arranged nuclei, collagen-rich matrix (blue), and eosinophilic cysts (Fig. 4A; Supplementary Fig. S4A). Fewer PCG2; Gli1nlacZ mice developed tumors than PCG2 (27% vs. 98%), with an average onset of 6 months, while the remaining animals succumbed to old age (n = 10; P < 0.0001; Fig. 4B). Along with tumor latency, penetrance, and pathology, the immunophenotype of PCG2; Gli1nlacZ tumors was also markedly different. Although myogenic marker expression was still absent in PCG2; Gli1nlacZ tumors, which were still highly proliferative, the Nkx2.2+, Cd99+ signature was not maintained (Fig. 4A and C; Supplementary Fig. S4B). To determine whether phenotypic outcome can be associated with level of Hedgehog pathway activation, we compared Gli1 abundance between tumors from PCG2, PCG2; Gli1nlacZ, and PtchlacZ mice. We found that Gli1 induction in PCG2; Gli1nlacZ tumors was less than half the level detected in PCG2. Surprisingly, induction in PCG2 tumors did not differ considerably from PtchlacZ eRMS using this method (Supplementary Fig. S4C). These findings indicate that a threshold of Hedgehog induction must be reached for PCG2 features to be expressed, and that the level of pathway induction alone does not likely distinguish undifferentiated and myogenic sarcomas.
Gli1 gene dosage influenced tumor phenotype. A, Hematoxylin and eosin (H&E) staining, trichrome staining, and immunostaining on PCG2; Gli1nlaZ tumors for indicated markers. B, Kaplan–Meier survival curves for PCG2 and PCG2; Gli1nlaZ mice. C, Mitotic and apoptotic indices for PtchlacZ, PCG2; Gli1nlacZ, and PCG2 tumors. Bars, average ± SEM. D, Immunostaining on tumors from PtchlacZ, PCG2; Gli1nlacZ, and PCG2 mice for indicated markers. Scale bars, 50 μm.
We chose to further characterize the epithelial-like histology of PCG2; Gli1nlacZ tumors. Although absent in PtchlacZ eRMS, keratins had strong membranous localization in PCG2; Gli1nlacZ tumors, consistent with epithelial organization. Cytoplasmic keratin was seen occasionally in PCG2 tumor cells (Fig. 4D). For corroboration we evaluated β-Catenin localization. Similar to keratins, expression of β-Catenin was below the detection limit in PtchlacZ tumors, but showed strong enrichment at cell borders in PCG2; Gli1nlacZ tumors (Fig. 4D). In contrast, PCG2 tumors showed intermittent nuclear enrichment of β-Catenin (Fig. 4D). Interestingly, aggressiveness in a subset of Ewing sarcoma with more frequent nuclear β-Catenin appears to be driven through activation of Wnt/β-Catenin signaling concomitant with upregulation of the extracellular matrix component and known EWS-FLI1 target gene, Tenascin-C (TNC; refs. 39, 40). Scattered enrichment of Tnc was detected in PCG2 tumors, but not in those from PtchlacZ or PCG2; Gli1nlacZ mice (Fig. 4D). Together these findings suggest a possible mechanism supporting the aggressive growth kinetics of PCG2 sarcomas, but not a likely parallel connecting them with the broader collection of Ewing sarcoma.
PCG2 sarcomas emerge from a naïve mesenchymal ancestor
Knowledge of the ancestral cell type for a given tumor provides key insight into the specific molecular landscape initially permissive to oncogenic transformation, details that could be powerful in the design of more effective treatment options (41). Therefore, we sought to elucidate a prospective cellular origin for PCG2 sarcomas. To achieve this end, we performed genetically inducible fate mapping to mark Pcp2-lineage cells with indelible reporters that express either cytoplasmic eYFP (R26ReYFP) or tdTomato (Ai9). At embryonic day 11.5 (e11.5), eYFP reporter activity (red) was largely detected in forelimb level dermamyotome cells that expressed Pax7/3 (green; Fig. 5A–D). Lineage cells were less frequently seen in Myf5 or Myogenin expression domains (green; Fig. 5A–D). In addition, we evaluated neural crest marker Foxd3, but found no overlap with lineage cells (Supplementary Fig. S5A and S5A’). These findings suggest that Pcp2-lineage cells are multipotent mesenchymal progenitors that predominantly generate skeletal muscle, adipocytes, and satellite cells.
Prospective cell-of-origin for PCG2 sarcomas. A–D, Immunostaining on transverse sections from Pcp2-cre; R26ReYFP embryos (e11.5) shows Pcp2-lineage cells in the transient dermamyotome (white outline, A and C) at the forelimb level express Pax7, Pax3 more frequently than Myogenin or Myf5, as indicated. E and E’, Direct fluorescence for tdTomato indicates Pcp2-lineage cells (red) contribute extensively to musculature (blue outline) in adult Pcp2-cre; Ai9 mice. Immunostaining of cross sections from proximal hind limb (black dotted line) revealed that Pcp2-lineage cells become both mature myofibers (F) and PPARγ+ adipocytes (F’). G, Hematoxylin and eosin (H&E), immunophenotype, and survival for Myf5Cre; CLEG2 (M5G2) mice. White arrowheads indicate double labeled cells. Scale bars represent distances indicated.
Whole-mount imaging of skeletal muscle in young adult Pcp2-cre; Ai9 mice showed that myofibers of the head, neck, and proximal limbs were often tdTomato+ (Fig. 5E and E’). In tissue sections, myofibers and PPARγ+ adipocytes were clearly labeled, consistent with Pax7 expression in Pcp2-lineage cells at e11.5 (Fig. 5F and F’). In postnatal tibialis anterior muscle, which was used in RNAseq studies, only few interstitial cells were tdTomato+, indicating limited cre activity in this tissue (Supplementary Fig. S5B and S5B’). We also targeted GLI2A to precursors of muscle and fat using Myf5Cre mice. Approximately 50% of the resulting M5G2 animals formed aggressive soft tissue tumors by 4 weeks of age that were indistinguishable from PCG2 tumors in terms of pathology and molecular marker profile (n = 11; P = 0.0016; Fig. 5G). However, when activated in postnatal satellite cells at P3 or P7 using Pax7-creER; Ai9 mice, GLIA did not induce tumor formation by 18 months, even though recombination was detected with the Ai9 reporter allele (n = 10, not shown). Altogether, these findings argue that the undifferentiated soft tissue sarcomas seen in PCG2 mice likely originate in common progenitors of myogenic and adipogenic cells.
Discussion
Collectively, our description of PCG2 mice demonstrates that GLI2A is sufficient to induce formation of sarcomas that are distinctly different from Hedgehog-driven eRMS. Instead, PCG2 tumors displayed an unexpected assortment of pathologic, immunophenotypic, transcriptomic, and ultrastructural features typical of human Ewing sarcoma and Ewing-like tumors, the latter of which receive the same standard-of-care treatment as bona fide Ewing sarcoma. Our work argues that PCG2 sarcomas emerge from a multipotent mesenchymal origin cell, which, in our hands, could not be transformed by EWS-FLI1, and demonstrates that Gli1 dosage can determine sarcoma phenotypic outcome. Altogether, we have provided a novel genetically defined mouse model with potential preclinical utility. Despite the perceived similarity shared by Ewing sarcoma and Ewing-like tumors with those formed by PCG2 mice, it is essential to emphasize that the latter are decidedly not a mouse equivalent to Ewing sarcoma, as they are not defined by an EWS-ETS gene fusion. However, PCG2 mice do have utility to inform more generally about the etiology of undifferentiated soft tissue sarcomas, especially those that display Ewing-like features but lack expected gene fusions and representative animal models.
A clear impediment to furthering knowledge of sarcoma biology, particularly in the case of Ewing sarcoma, has been limited availability of faithful mouse models. Difficulty with producing a true, genetically defined Ewing sarcoma mouse may stem from nonconservation of microsatellites, which determine functional output of EWS-FLI1, and/or lethality induced by the oncogene (9, 30, 42). Promising work with flank grafting of isolated progenitor cells ectopically expressing EWS-FLI1 has yielded mouse “Ewing-like” tumors (17, 18), which leaves hope that the right design may ultimately yield a true “Ewing” mouse. Future approaches should consider harnessing direct activation of known downstream effectors of EWS-FLI1 as an alternative or in combination with the oncogene. GLI1 (43–45), one of three zinc finger transcription factors that regulate vertebrate Hedgehog target genes (46), is an established target and effector of EWS-FLI1, hence the sensitivity of Ewing sarcoma cell lines to GLI1/2 inhibitor arsenic trioxide (47).
It is likely that the subset of features shared by PCG2 sarcomas and Ewing sarcoma is attributable to this core connection, utilization of GLI activator function, which has not been observed in any other sarcoma variant beyond eRMS. Even though PCG2 mice leverage activity of a truncated human GLI2, GLI1 and GLI2 are known to operate through a consensus binding sequence to regulate common target genes such as PTCH1 (46). Because ubiquitous activation of Hedgehog signaling in mice, due to either widespread deactivation of receptor Patched1 (Ptch1; refs. 5, 6) or constitutively activated Smoothened (SmoM2), is so far only known to cause eRMS and not other sarcoma types (7, 8), suggesting these different modes of activating the pathway induce different target gene modules and ultimately very different tumor phenotypes.
However, it is critical to emphasize that in addition to the noted similarities, differences between genuine Ewing tumors and those formed by PCG2 mice were also observed. The most significant of which is that the latter do not appear to harbor any gene fusions, definitely not a fusion relevant in the context of the human sarcoma variant. Given this important difference, it is understandable that the resulting gene expression profiles different in a meaningful way; several notable “Ewing” genes were strongly downregulated or not detected in in our model, such as Prkcb and Nr0b1. Whether the epigenetic landscape in PCG2 tumors would be reflective of what has been documented in Ewing is to be determined, but given that EWS-FLI1 is thought to be a major driver of establishing that signature (42), there may not be much similarity in that regulatory dimension. When considering these distinguishing differences, it becomes clear why PCG2 specimens would group most closely with a subset of Ewing and “Ewing-like” tumors in PCA of their respective transcriptomic profiles. Furthermore, the most telling comparison will be testing the efficacy of therapeutic agents in PCG2 mice that one day may be considered for Ewing treatment.
Our work also demonstrated that Pcp2-lineage tumor origin cells are naïve mesenchymal progenitors of muscle and fat, not neural crest. Interestingly, the prospective cellular origin for Ewing sarcoma is thought to be either mesenchymal stem/progenitor cells, possibly from bone marrow, or cells in the neural crest lineage. This thinking is driven by Ewing sarcoma tumor pathology that can exhibit mesenchymal and/or neuroectodermal elements. The latter being more prevalent in peripheral neuroectodermal tumors (PNET) that share the hallmark t(11; 22) rearrangement and therefore the Ewing sarcoma classification. Whether originating cell type is an influential factor that helps shape tumor phenotype (osseous, extraosseous, or PNETs in the case of Ewing sarcoma) is unclear. Why EWS-FLI1 failed to induce tumors in Pcp2-cre; EWSEWS-FLI1 mice is also not apparent, but may suggest that an additional genetic alteration is needed to drive transformation over lethality in Pcp2-lineage mesenchymal progenitor cells.
Altogether, PCG2 represents the first genetically inducible, Hedgehog-driven mouse model of undifferentiated sarcoma and is therefore a unique tool for exploring tumor initiation/progression, and perhaps even for preclinical studies. Our data argue that future studies should evaluate the efficacy of targeting Hedgehog signaling at the transcriptional level. For example, use of inhibitors of BET bromodomain proteins (48) in the clinical management of Ewing sarcoma warrants consideration. Indeed, work from our laboratory and others has already demonstrated efficacy for these compounds in suppressing EWS-FLI1–mediated transcriptional regulation and Ewing sarcoma tumor growth (49–52).
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Authors' Contributions
Conception and design: J.T. Fleming, C. Chiang
Development of methodology: J.T. Fleming, A.A. Dlugosz, C. Chiang
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): J.T. Fleming, E. Brignola, L. Chen, A.N. Ermilov, A.A. Dlugosz, C. Chiang
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): J.T. Fleming, L. Chen, Y. Guo, S. Zhao, Q. Wang, B. Li, H. Correa, C. Chiang
Writing, review, and/or revision of the manuscript: J.T. Fleming, B. Li, C. Chiang
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): J.T. Fleming
Study supervision: J.T. Fleming, C. Chiang
Acknowledgments
The authors wish to extend many thanks to those who shared their reagents or assistance during the course of this study, including Javed Khan (NCI), Patricia Labosky (NIH), Dietmar Vestweber (Max Plank), Herald Erickson (Duke), and Leighton Stein (RPCI), and Vanderbilt colleagues: Patrick Grohar (VARI), Scott Borinstein, and Cheryl Coffin. We thank several VUMC core resources, TPSR, DHSR, CISR, and VANTAGE. We are also grateful to the DSHB, established in part through generosity from the NIH. This study was supported by grants to A. Dlugosz from the NIH RO1 CA087837 and University of Michigan Rogel Cancer Center Support Grant P30 CA046592, and to C. Chiang from the NIH RO1 NS097898, the Vanderbilt-Ingram Cancer Center Support Grant P30 CA068485 and Alex's Lemonade Stand Foundation.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Footnotes
Note: Supplementary data for this article are available at Molecular Cancer Research Online (http://mcr.aacrjournals.org/).
- Received February 20, 2018.
- Revision received June 6, 2018.
- Accepted January 22, 2019.
- Published first January 25, 2019.
- ©2019 American Association for Cancer Research.