We previously identified two distinct molecular subtypes of osteosarcoma through gene expression profiling. These subtypes are associated with distinct tumor behavior and clinical outcomes. Here, we ...describe mechanisms that give rise to these molecular subtypes. Using bioinformatic analyses, we identified a significant association between deregulation of the retinoblastoma (RB)-E2F pathway and the molecular subtype with worse clinical outcomes. Xenotransplantation models recapitulated the corresponding behavior for each osteosarcoma subtype; thus, we used cell lines to validate the role of the RB-E2F pathway in regulating the prognostic gene signature. Ectopic RB resets the patterns of E2F regulated gene expression in cells derived from tumors with worse clinical outcomes (molecular phenotype 2) to those comparable with those observed in cells derived from tumors with less aggressive outcomes (molecular phenotype 1), providing a functional association between RB-E2F dysfunction and altered gene expression in osteosarcoma. DNA methyltransferase and histone deacetylase inhibitors similarly reset the transcriptional state of the molecular phenotype 2 cells from a state associated with RB deficiency to one seen with RB sufficiency. Our data indicate that deregulation of RB-E2F pathway alters the epigenetic landscape and biological behavior of osteosarcoma.
Background: Gene expression signatures define prognostically significant osteosarcoma phenotypes.
Results: Deregulation of the RB-E2F pathway establishes more aggressive phenotype. Inhibitors of DNA and chromatin remodeling promote comparable transcriptional changes as genetic restoration of RB.
Conclusion: Aberrant RB-E2F pathway alters epigenetic landscape and biological behavior of osteosarcoma.
Significance: Epigenetic remodeling regulated by RB-E2F gives rise to patterns of gene expression that are associated with different biological behavior and progression of osteosarcoma.
Osteosarcoma has a guarded prognosis. A major hurdle in developing more effective osteosarcoma therapies is the lack of disease-specific biomarkers to predict risk, prognosis, or therapeutic ...response. Exosomes are secreted extracellular microvesicles emerging as powerful diagnostic tools. However, their clinical application is precluded by challenges in identifying disease-associated cargo from the vastly larger background of normal exosome cargo. We developed a method using canine osteosarcoma in mouse xenografts to distinguish tumor-derived from host-response exosomal messenger RNAs (mRNAs). The model allows for the identification of canine osteosarcoma-specific gene signatures by RNA sequencing and a species-differentiating bioinformatics pipeline. An osteosarcoma-associated signature consisting of five gene transcripts (SKA2, NEU1, PAF1, PSMG2, and NOB1) was validated in dogs with spontaneous osteosarcoma by real-time quantitative reverse transcription PCR (qRT-PCR), while a machine learning model assigned dogs into healthy or disease groups. Serum/plasma exosomes were isolated from 53 dogs in distinct clinical groups (“healthy”, “osteosarcoma”, “other bone tumor”, or “non-neoplastic disease”). Pre-treatment samples from osteosarcoma cases were used as the training set, and a validation set from post-treatment samples was used for testing, classifying as “osteosarcoma detected” or “osteosarcoma-NOT detected”. Dogs in a validation set whose post-treatment samples were classified as “osteosarcoma-NOT detected” had longer remissions, up to 15 months after treatment. In conclusion, we identified a gene signature predictive of molecular remissions with potential applications in the early detection and minimal residual disease settings. These results provide proof of concept for our discovery platform and its utilization in future studies to inform cancer risk, diagnosis, prognosis, and therapeutic response.
Intrinsic biological fluctuation and/or measurement error can obscure the association of gene expression patterns between RNA and protein levels. Appropriate normalization of reverse-transcription ...quantitative PCR (RT-qPCR) data can reduce technical noise in transcript measurement, thus uncovering such relationships. The accuracy of gene expression measurement is often challenged in the context of cancer due to the genetic instability and “splicing weakness” involved. Here, we sequenced the poly(A) cancer transcriptome of canine osteosarcoma using mRNA-Seq. Expressed sequences were resolved at the level of two consecutive exons to enable the design of exon-border spanning RT-qPCR assays and ranked for stability based on the coefficient of variation (
CV
). Using the same template type for RT-qPCR validation, i.e. poly(A) RNA, avoided skewing of stability assessment by circular RNAs (circRNAs) and/or rRNA deregulation. The strength of the relationship between mRNA expression of the tumour marker S100A4 and its proportion score of quantitative immunohistochemistry (qIHC) was introduced as an experimental readout to fine-tune the normalization choice. Together with the essential logit transformation of qIHC scores, this approach reduced the noise of measurement as demonstrated by uncovering a highly significant, strong association between mRNA and protein expressions of
S100A4
(Spearman’s coefficient
ρ
= 0.72 (
p
= 0.006)).
Key messages
• RNA-seq identifies stable pairs of consecutive exons in a heterogeneous tumour.
• Poly(A) RNA templates for RT-qPCR avoid bias from circRNA and rRNA deregulation.
• HNRNPL
is stably expressed across various cancer tissues and osteosarcoma.
• Logit transformed qIHC score better associates with mRNA amount.
• Quantification of minor
S100A4
mRNA species requires poly(A) RNA templates and dPCR.
Osteosarcoma is an aggressive tumor of the bone that primarily affects young adults and adolescents. Osteosarcoma is characterized by genomic chaos and heterogeneity. While inactivation of tumor ...protein p53 (TP53) is nearly universal other high frequency mutations or structural variations have not been identified. Despite this genomic heterogeneity, key conserved transcriptional programs associated with survival have been identified across human, canine and induced murine osteosarcoma. The epigenomic landscape, including DNA methylation, plays a key role in establishing transcriptional programs in all cell types. The role of epigenetic dysregulation has been studied in a variety of cancers but has yet to be explored at scale in osteosarcoma. Here we examined genome-wide DNA methylation patterns in 24 human and 44 canine osteosarcoma samples identifying groups of highly correlated DNA methylation marks in human and canine osteosarcoma samples. We also link specific DNA methylation patterns to key transcriptional programs in both human and canine osteosarcoma. Building on previous work, we built a DNA methylation-based measure for the presence and abundance of various immune cell types in osteosarcoma. Finally, we determined that the underlying state of the tumor, and not changes in cell composition, were the main driver of differences in DNA methylation across the human and canine samples.
Genome wide comparison of DNA methylation patterns in osteosarcoma across two species lays the ground work for the exploration of DNA methylation programs that help establish conserved transcriptional programs in the context of varied mutational landscapes.
•Several shared DNA methylation profiles between human (n = 24) and canine (n = 44) primary osteosarcomas•Shared DNA methylation profiles align with shared transcriptional profiles•New method to quantify immune cell infiltration using genome wide DNA methylation profiles•Epigentic role in the conserfed transriptional profiles across the varied mutational landscapes seen in osteosarcoma.
The etiology of hemangiosarcoma remains incompletely understood. Its common occurrence in dogs suggests predisposing factors favor its development in this species. These factors could represent a ...constellation of heritable characteristics that promote transformation events and/or facilitate the establishment of a microenvironment that is conducive for survival of malignant blood vessel-forming cells. The hypothesis for this study was that characteristic molecular features distinguish hemangiosarcoma from non-malignant endothelial cells, and that such features are informative for the etiology of this disease.
We first investigated mutations of VHL and Ras family genes that might drive hemangiosarcoma by sequencing tumor DNA and mRNA (cDNA). Protein expression was examined using immunostaining. Next, we evaluated genome-wide gene expression profiling using the Affymetrix Canine 2.0 platform as a global approach to test the hypothesis. Data were evaluated using routine bioinformatics and validation was done using quantitative real time RT-PCR.
Each of 10 tumor and four non-tumor samples analyzed had wild type sequences for these genes. At the genome wide level, hemangiosarcoma cells clustered separately from non-malignant endothelial cells based on a robust signature that included genes involved in inflammation, angiogenesis, adhesion, invasion, metabolism, cell cycle, signaling, and patterning. This signature did not simply reflect a cancer-associated angiogenic phenotype, as it also distinguished hemangiosarcoma from non-endothelial, moderately to highly angiogenic bone marrow-derived tumors (lymphoma, leukemia, osteosarcoma).
The data show that inflammation and angiogenesis are important processes in the pathogenesis of vascular tumors, but a definitive ontogeny of the cells that give rise to these tumors remains to be established. The data do not yet distinguish whether functional or ontogenetic plasticity creates this phenotype, although they suggest that cells which give rise to hemangiosarcoma modulate their microenvironment to promote tumor growth and survival. We propose that the frequent occurrence of canine hemangiosarcoma in defined dog breeds, as well as its similarity to homologous tumors in humans, offers unique models to solve the dilemma of stem cell plasticity and whether angiogenic endothelial cells and hematopoietic cells originate from a single cell or from distinct progenitor cells.
Lymphocytic choriomeningitis virus (LCMV) is the prototypic arenavirus and a natural mouse pathogen. LCMV-Armstrong, an acutely resolved strain, and LCMV-clone 13, a mutant that establishes chronic ...infection, have provided contrasting infection models that continue to inform the fundamental biology of T cell differentiation, regulation of exhaustion, and response to checkpoint blockade. In this study, we report the isolation and characterization of LCMV-Minnesota (LCMV-MN), which was naturally transmitted to laboratory mice upon cohousing with pet shop mice and shares 80-95% amino acid homology with previously characterized LCMV strains. Infection of laboratory mice with purified LCMV-MN resulted in viral persistence that was intermediate between LCMV-Armstrong and -clone 13, with widely disseminated viral replication and viremia that was controlled within 15-30 d, unless CD4 T cells were depleted prior to infection. LCMV-MN-responding CD8
T cells biased differentiation toward the recently described programmed death-1 (PD-1)
CXCR5
Tim-3
stemlike CD8
T cell population (also referred to as progenitor exhausted T cells) that effectuates responses to PD-1 blockade checkpoint inhibition, a therapy that rejuvenates responses against chronic infections and cancer. This subset resembled previously characterized PD-1
TCF1
stemlike CD8
T cells by transcriptional, phenotypic, and functional assays, yet was atypically abundant. LCMV-MN may provide a tool to better understand the breadth of immune responses in different settings of chronic Ag stimulation as well as the ontogeny of progenitor exhausted T cells and the regulation of responsiveness to PD-1 blockade.
Osteosarcoma is an aggressive primary bone tumor in humans and is among the most common cancer afflicting dogs. Despite surgical advancements and intensification of chemo- and targeted therapies, the ...survival outcome for osteosarcoma patients is, as of yet, suboptimal. The presence of metastatic disease at diagnosis or its recurrence after initial therapy is a major factor for the poor outcomes. It is thought that most human and canine patients have at least microscopic metastatic lesions at diagnosis. Osteosarcoma in dogs occurs naturally with greater frequency and shares many biological and clinical similarities with osteosarcoma in humans. From a genetic perspective, osteosarcoma in both humans and dogs is characterized by complex karyotypes with highly variable structural and numerical chromosomal aberrations. Similar molecular abnormalities have been observed in human and canine osteosarcoma. For instance, loss of
and
regulated pathways are common. While there are several oncogenes that are commonly amplified in both humans and dogs, such as
and
, no commonly activated proto-oncogene has been identified that could form the basis for targeted therapies. It remains possible that recurrent aberrant gene expression changes due to gene amplification or epigenetic alterations could be uncovered and these could be used for developing new, targeted therapies. However, the remarkably high genomic complexity of osteosarcoma has precluded their definitive identification. Several advantageous murine models of osteosarcoma have been generated. These include spontaneous and genetically engineered mouse models, including a model based on forward genetics and transposon mutagenesis allowing new genes and genetic pathways to be implicated in osteosarcoma development. The proposition of this review is that careful comparative genomic studies between human, canine and mouse models of osteosarcoma may help identify commonly affected and targetable pathways for alternative therapies for osteosarcoma patients. Translational research may be found through a path that begins in mouse models, and then moves through canine patients, and then human patients.
Abstract The heterogeneous and chaotic nature of osteosarcoma has confounded accurate molecular classification, prognosis, and prediction for this tumor. The occurrence of spontaneous osteosarcoma is ...largely confined to humans and dogs. While the clinical features are remarkably similar in both species, the organization of dogs into defined breeds provides a more homogeneous genetic background that may increase the likelihood to uncover molecular subtypes for this complex disease. We thus hypothesized that molecular profiles derived from canine osteosarcoma would aid in molecular subclassification of this disease when applied to humans. To test the hypothesis, we performed genome wide gene expression profiling in a cohort of dogs with osteosarcoma, primarily from high-risk breeds. To further reduce inter-sample heterogeneity, we assessed tumor-intrinsic properties through use of an extensive panel of osteosarcoma-derived cell lines. We observed strong differential gene expression that segregated samples into two groups with differential survival probabilities. Groupings were characterized by the inversely correlated expression of genes associated with ‘G2/M transition and DNA damage checkpoint’ and ‘microenvironment-interaction’ categories. This signature was preserved in data from whole tumor samples of three independent dog osteosarcoma cohorts, with stratification into the two expected groups. Significantly, this restricted signature partially overlapped a previously defined, predictive signature for soft tissue sarcomas, and it unmasked orthologous molecular subtypes and their corresponding natural histories in five independent data sets from human patients with osteosarcoma. Our results indicate that the narrower genetic diversity of dogs can be utilized to group complex human osteosarcoma into biologically and clinically relevant molecular subtypes. This in turn may enhance prognosis and prediction, and identify relevant therapeutic targets.
Interleukin-8 (IL-8) gene expression is highly up-regulated in canine hemangiosarcoma (HSA); however, its role in the pathogenesis of this disease is unknown. We investigated the expression of IL-8 ...in canine HSA tissues and cell lines, as well and the effects of IL-8 on canine HSA in vitro, and in vivo using a mouse xenograft model for the latter. Constitutive expression of IL-8 mRNA, IL-8 protein, and IL-8 receptor were variable among different tumor samples and cell lines, but they showed stable steady states in each cell line. Upon the addition of IL-8, HSA cells showed transient intracellular calcium fluxes, suggesting that their IL-8 receptors are functional and that IL-8 binding activates relevant signaling pathways. Yet, neither addition of exogenous IL-8 nor blockade of endogenous IL-8 by neutralizing anti-IL-8 antibody (α-IL-8 Ab) affected HSA cell proliferation or survival in vitro. To assess potential effects of IL-8 in other tumor constituents, we stratified HSA cell lines and whole tumor samples into “IL-8 high” and “IL-8 low” groups. Genome-wide gene expression profiling showed that samples in the “IL-8 high” tumor group were enriched for genes associated with a “reactive microenvironment,” including activation of coagulation, inflammation, and fibrosis networks. Based on these findings, we hypothesized that the effects of IL-8 on these tumors were mostly indirect, regulating interactions with the microenvironment. This hypothesis was supported by in vivo xenograft experiments where survival and engraftment of tumor cells was inhibited by administration of neutralizing α-IL-8 Ab. Together, our results suggest that IL-8 contributes to establishing a permissive microenvironment during the early stages of tumorigenesis in HSA.
•IL-8 is expressed in canine hemangiosarcoma tumor samples and cell lines.•IL-8 transduces a relevant biological signal in canine hemangiosarcoma cells.•IL-8 gene signature is associated with reactive tumor microenvironments.•IL-8 potentiates tumor cell survival and engraftment into host tissues.•Canine hemangiosarcoma provides a unique comparative model for IL-8 studies.
Osteosarcoma (OS) is a heterogeneous and rare disease with a disproportionate impact because it mainly affects children and adolescents. Lamentably, more than half of patients with OS succumb to ...metastatic disease. Clarification of the etiology of the disease, development of better strategies to manage progression, and methods to guide personalized treatments are among the unmet health needs for OS patients. Progress in managing the disease has been hindered by the extreme heterogeneity of OS; thus, better models that accurately recapitulate the natural heterogeneity of the disease are needed. For this study, we used cell lines derived from two spontaneous canine OS tumors with distinctly different biological behavior (OS-1 and OS-2) for heterotypic in vivo modeling that recapitulates the heterogeneous biology and behavior of this disease. Both cell lines demonstrated stability of the transcriptome when grown as orthotopic xenografts in athymic nude mice. Consistent with the behavior of the original tumors, OS-2 xenografts grew more rapidly at the primary site and had greater propensity to disseminate to lung and establish microscopic metastasis. Moreover, OS-2 promoted formation of a different tumor-associated stromal environment than OS-1 xenografts. OS-2-derived tumors comprised a larger percentage of the xenograft tumors than OS-1-derived tumors. In addition, a robust pro-inflammatory population dominated the stromal cell infiltrates in OS-2 xenografts, whereas a mesenchymal population with a gene signature reflecting myogenic signaling dominated those in the OS-1 xenografts. Our studies show that canine OS cell lines maintain intrinsic features of the tumors from which they were derived and recapitulate the heterogeneous biology and behavior of bone cancer in mouse models. This system provides a resource to understand essential interactions between tumor cells and the stromal environment that drive the progression and metastatic propensity of OS.