Gene expression profile comparison of primary and pulmonary metastatic lesions in a dog with appendicular osteosarcoma and hypertrophic osteopathy

Gene expression profile comparison of primary and pulmonary metastatic lesions in a dog with appendicular osteosarcoma and hypertrophic osteopathy | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Gene expression profile comparison of primary and pulmonary metastatic lesions in a dog with appendicular osteosarcoma and hypertrophic osteopathy Keita Kitagawa, Jessica Dryfhout, Alexander I. Engleberg, Ya-Ting Yang, and 2 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4438777/v1 This work is licensed under a CC BY 4.0 License Status: Under Review Version 1 posted 11 You are reading this latest preprint version Abstract Hypertrophic osteopathy (HO) is a paraneoplastic syndrome, and the most notable cause in dogs is pulmonary metastatic osteosarcoma (OSA). Although many molecular factors in canine OSA have been shown in metastasis, little is known about the gene expression profile of HO secondary to metastatic OSA. Therefore, the purpose of this study was to compare the gene expression profiles between primary and metastatic OSA lesions from the same dog and to look for gene expression changes that can elucidate the molecular mechanism of metastases and HO. Tumoral samples were obtained from a 2-year-old, intact male, Labrador retriever. At the first visit, the patient presented with an appendicular OSA as the primary lesion. About 10 months later, the dog developed HO due to a single pulmonary metastasis. Using these primary and metastatic samples from the same dog, as well as normal canine osteoblasts, we investigated the gene expression profiling using the NanoString nCounter® Canine IO panel. A total of 180 differentially expressed genes were identified between malignant OSA cells and non-malignant canine osteoblasts. Furthermore, 5 genes (CCL17, VEGFC, C3, C4BPA, and FOS) were differentially expressed in comparison between primary and metastatic OSA samples. CCL17 and VEGFC were upregulated in the primary lesion compared to the metastatic lesion, while C3, C4BPA, and FOS were downregulated in the primary lesion relative to the metastatic lesion. Given that the metastatic lesion was relevant to the development of HO, the different gene expression profiles may be relevant to understanding the pathophysiology of HO. Hypertrophic osteopathy NanoString osteosarcoma paraneoplastic syndrome gene expression profile Figures Figure 1 Figure 2 BACKGROUND Canine hypertrophic osteopathy (HO) is a paraneoplastic syndrome associated with pulmonary tumoral lesions and is characterized by the symmetrical formation of periosteal new bone along the diaphysis of the long bones of the appendicular skeleton( 1 – 4 ). Clinical signs include lameness, lethargy, localized pain, and soft tissue swelling surrounding the affected bones( 1 – 4 ). The diagnosis of HO is made based on the clinical history, clinical signs, and radiographical changes( 2 – 4 ). Metastasectomy, a surgical intervention to remove one or more metastatic lesions, can result in the control of HO by achieving the resolution of the clinical signs and the radiographical changes.( 2 , 5 , 6 ). The pathogenesis of HO is not currently fully understood. The veterinary literature in this regard is null, and the scarce references in human literature propose two major hypotheses explaining the development of HO: neurological and humoral pathways ( 7 – 9 ). Resolution of symptoms following vagotomy empirically suggests neurologic involvement in HO, whereas humoral pathway involvement is supported by increased release of serum growth factors, including growth hormone, growth hormone-releasing hormone, platelet-derived growth factor (PDGF), and vascular endothelial growth factor (VEGF) by the metastatic lesion ( 8 – 10 ). Canine appendicular osteosarcoma (OSA) is the most common form of malignant bone cancer in dogs, which is locally aggressive and highly metastatic; OSA is also the most notable cause of paraneoplastic HO in dogs ( 2 , 11 ). Regarding the pathophysiology of canine OSA, many molecular and genetic factors have been investigated to explain the development of metastasis ( 12 , 13 ). The overexpression of delta Np63, a transcription factor belonging to p53 family, has been reported to be correlated with pulmonary metastases via increased secretion of VEGFA through activation of STAT3( 12 ). In addition, the downregulation of STAT3 decreases the expression of VEGF, and conversely, activation of STAT3 increases VEGF expression, and may thus facilitate metastasis in OSA ( 13 ). Also, in addition to STAT3 and delta Np63, TGF beta also regulates VEGF secretion. The inhibition of TGF beta1 decreases VEGF secretion, suggesting TGF beta1 is associated with angiogenesis via VEGF( 14 ). TGF beta signaling is suggested to promote distant metastasis by activating cell migration and neoangiogenesis in OSA ( 14 ). Furthermore, the downregulation of CXCR4 by zoledronate treatment alters the pattern of metastasis in OSA, curtailing metastasis to the lungs, but, not to other organs, such as the nerves and skin, suggesting CXCR4 has some contribution to pulmonary metastasis in OSA ( 15 ). Molecular investigations in metastasis in canine OSA until now have focused on the pathophysiology of the primary OSA lesion. Lately, however, proteomics for OSA has shown differences in protein expression between primary and metastatic lesions ( 16 ). Identifying the differences in gene expression profiling between primary and metastatic lesions would provide us with critical information for further understanding the pathophysiology of HO, as well as the metastatic process in OSA. Hence, the purpose of this study was to compare the gene expression profiles between primary and metastatic OSA lesions from the same dog who developed HO due to metastatic OSA. MATERIALS AND METHODS PRIMARY AND METASTATIC OSA SAMPLES Primary and metastatic OSA samples were collected from a 2-year-old, intact male, Labrador retriever that presented to the Michigan State University (MSU) Veterinary Medical Center. Both appendicular and solitary pulmonary metastatic lesions were surgically removed, and histopathologic diagnoses were confirmed by board-certified pathologists at MSU. CELL CULTURE Three normal osteoblast strains were used in the current study as controls. Among them, one canine osteoblast strain (CnOb) was purchased from Cell Application Inc. (San Diego, CA, USA). Additional two canine osteoblast cell strains were derived from two different healthy beagles in our lab as described ( 17 ). All three CnOb cell strains were maintained with CnOb Growth Medium (Cell Application Inc., San Diego, CA, USA) and supplemented with 10% fetal bovine serum (Thermo Fisher Scientific, Waltham, MA, USA), antibiotics (0.1% gentamycin, Life Technologies, Carlsbad, CA, USA) and incubated in a humidified incubator at 37˚C with 5% CO 2 . These three canine osteoblasts strains were shown to express common osteoblast markers, including osteocalcin and collagen type I alpha chain (COL1A1) via quantitative real-time PCR previously by our team ( 17 ). RNA ISOLATION Formalin-fixed paraffin-embedded (FFPE) canine OSA tissues from primary and metastatic lesions were obtained from the MSU Veterinary Diagnostic Library. Total RNA was isolated from FFPE tissue blocks using the RNeasy FFPE Kit (Qiagen, Hilden, Germany) according to manufacturer’s instructions. RNA was also isolated from canine osteoblasts using mirVana miRNA Isolation Kit with phenol (Thermo Fisher Scientific, Waltham, MA, USA). Each sample was processed in triplicate. Isolated RNA was immediately stored at -80˚C for further studies. DATA ANALYSIS BY NANOSTRING TECHNOLOGY Gene expression was evaluated using the NanoString nCounter® Canine IO (immuno-oncology) gene-expression panel comprising of 780 genes, developed using the CanFam3.1 reference genome. Targeted expression profiling was performed using isolated RNA from OSA samples and canine osteoblast cell lines. Differential gene expression analysis was performed comparing affected OSA samples to non-tumoral control osteoblast cells. Further analysis compared expression levels between primary and metastatic OSA samples. STATISTICAL ANALYSIS Data normalization, quality control, and functional and pathway analysis were performed by importing nCounter data onto the ROSALIND® analysis platform ( https://www.rosalind.bio ). Normalization, fold changes, and p-values were calculated in-platform following the nCounter® advanced analysis protocol comparing the geometric mean of housekeeping probes to sample counts from the same lane of the array. Housekeeping genes used for normalization were chosen based on the geNorm algorithm as used in the NormqPCR R library.( 18 ) In addition to differential gene expression, gene set analysis (GSA) and cell type profiling were also performed using the ROSALIND platform( 19 ). Differential expression was determined as significant under the parameters of a fold-change >|1.5|and and an adjusted P-value < 0.05, using the Benjamini-Hochberg correction method. The ROSALIND platform performs enrichment analysis by referencing several database sources, such as Interpro, NCBI, MSigDB, REACTOME, and WikiPathways ( 20 – 25 ). Volcano plots and heatmaps were generated on ROSALIND comparing OSA cells to osteoblast cells and primary OSA to metastatic OSA. RESULTS CASE HISTORY A 2-year-old, intact male, Labrador retriever was presented for a progressive right hindlimb lameness of 2-months duration at MSU Veterinary Medical Center. The dog was non-weight bearing and moderate muscle atrophy of the right hindlimb was observed. There was a large firm mass of approximately 10 cm in diameter in the right proximal tibia, causing a decreased range of motion of the right stifle, and pitting edema distal to the stifle joint. For staging purposes, we obtained bloodwork, thoracic radiographs, pre- and post-contrast Computed tomography (CT) scan. Staging results did not reveal any remarkable abnormalities except for an irregular shaped, large mass associated with the caudal aspect of the proximal tibia by CT scan. The histopathology report of surgical biopsy sample obtained from the proximal tibia mass confirmed the diagnosis of chondroblastic OSA with a mitotic index of 77 in 10 high powered fields. The dog then underwent a right hindlimb amputation; however, adjuvant chemotherapy was not pursued at that time per the owner’s decision. The dog was examined 314 days postoperatively because of lethargy and decreased appetite. The dog had progressive lameness, pyrexia of 40.4˚C, soft tissue swelling and pain in the extremities. Orthopedic examination revealed a decreased range of motion in both elbow and carpal joints and pain upon extension of elbows bilaterally. Blood work revealed inflammatory leukogram, mild-moderate elevated ALP [154 U/L (normal range: 10–92 U/L)], and low urea nitrogen [8 mg/dL (normal range: 12–27 mg/dL)], otherwise within reference intervals. For the diagnostic imaging, thoracic radiographs revealed the presence of a large intrathoracic mass and a smoothly marginated periosteal proliferation along the right humerus. Also, in CT scan, periosteal proliferations were noted along the fore and remaining hindlimb including the right humerus, the left tibia, the left metatarsal bones, and the left femur. The CT scan further indicated a solitary mass in the cranial mediastinum or right cranial lung lobe. Fine needle aspiration of the mediastinal mass was done; the result of cytology report was consistent with OSA. Integrating these findings, the dog was diagnosed as HO secondary to pulmonary metastatic OSA. For the treatment of HO, the cranial mediastinal mass was surgically removed, and submitted for histopathology. The histopathology evaluation was consistent with chondroblastic OSA, likely metastatic. As adjuvant therapy following metastasectomy, carboplatin (250-300mg/m 2 IV every 3 weeks) was given a total of 6 doses and without remarkable adverse effects throughout the chemotherapy. Postoperative follow-up was done by physical examination with either radiographs or CT scan. Clinical improvement of HO was achieved 9 days after metastasectomy; the dog regained mobility without associated pain. The first postoperative CT scan, performed 58 days after metastasectomy, showed that the dog retained radiologic evidence of HO. The previously noted HO lesions persisted and most of them were slightly more severe in right thoracic limb involving the right scapula, right humerus, right metacarpal bones, and the bilateral radii, ulnae, and second metacarpal bones. Interestingly, the mild improvement of HO in the right scapula was found radiologically based on the CT scan at that time. No lameness nor pain on the legs were appreciated at the second postoperative CT scan (i.e., 128 days after metastasectomy). Palpation revealed no abnormal findings and the previously noted periosteal proliferation (i.e., HO lesion) confirmed no longer present in the right humerus by CT scan. GENE EXPRESSION PROFILE ANALYSIS We used the NanoString nCounter® Canine IO panel which has great sensitivity and specificity to elucidate the differentially expressed genes in OSA cells compared to non-tumoral osteoblasts( 26 , 27 ). The overview of the gene expression data is presented in the heatmap in Fig. 1 . Distinct expression profiles for both up- and downregulated genes across all samples can be appreciated in Fig. 1 a. The volcano plot shows all samples plotted as a functional fold change vs. p-value (Fig. 1 b and 2 b). Genes that exhibited a significant (p < 0.05) and 1.5-fold change in expression compared to the control group were selected. Overall, 180 genes were differentially regulated; 121 of these were upregulated, whereas 59 were downregulated in OSA compared to control (Fig. 1 b). Some of the top upregulated include genes involved in the immune response and cytokine production, such as SPP1, CXCL10, CXCR4, whereas other genes were associated with the cell signaling pathway, such as TNFRSF11A and KIT (Table S1 ). On the other hand, of the 59 downregulated genes, the most down-regulated genes included those involved in cell migration and the cell cycle, such as PTHR1, THBS1, CCND1, CCDN2, and CDKN1A among others (Table S2). Of special interest to us were the differences in gene expression between the primary and metastatic lesions. Although the number of differential genes expressed was small, the heatmap showed distinct expression profiles for both up- and downregulated genes. (Fig. 2 a). Overall, 5 genes were differently expressed. The CCL17 and VEGFC were significantly upregulated in primary lesion when compared to metastatic lesion, while FOS, C3, and C4BPA were significantly downregulated in primary lesion when compared to metastatic lesion (Fig. 2 b). Fold changes of these differential gene expressions are summarized in Table. 1. Through ROSALIND®, several differentially expressed signaling pathways were identified. Pathway significance was determined by a Global Significance Score (GSS) rating based on gene set analysis (Table 2 ). The GSS measures the extent of the difference of overall differential gene expression between primary and metastasis OSAs within a distinct signaling pathway. GSSs are positive regardless of upregulation or downregulation, with a large value indicating greater significance. The pathway with the highest GSS rating between primary and metastatic OSA was the Complement pathway with a score of 1.60. The other significant pathway was the Complement System pathway (1.41). The pathways with GSS ratings less than 1 included the NK cell functions pathway (0.92), the Angiogenesis pathway (0.85), the Chemokines pathway (0.77), the Wnt Signaling pathway (0.77), the Interleukins pathway (0.76), the Epigenetic Regulation pathway (0.76), the Regulation pathway (0.74) and the TGF-beta signaling pathway (0.74). The detailed information of these pathways was summarized in the Table S3. Table 1 Genes with significant differential expression between the Primary vs Metastatic OS in this patient Gene Fold Change Log 2 fold change CCL17 + 17.5 4.1 VEGFC + 3.7 1.9 FOS − 11 -3.5 C3 − 34 -5.1 C4BPA -48 -5.6 Table 2 Global significance ratings comparing overall differential expression of selected pathways between primary vs metastatic OSA. Term Name Global Significance Score Complement 1.6066 Complement System 1.415 NK Cell Functions 0.9216 Angiogenesis 0.8553 Chemokines 0.774 Wnt Signaling 0.7688 Interleukins 0.7613 Epigenetic Regulation 0.7612 Regulation 0.7441 TGF-beta Signaling 0.7385 DISCUSSION The clinical history of our OSA-HO case, including blood work results, radiographic findings, and the radiological resolution of HO after metastasectomy, are consistent with the previous studies ( 2 – 5 ). Therefore, it is an archetypical subject of analysis for further understanding the pathophysiology of HO and OSA in dogs. When comparing the gene expression profile of OSA vs. control, 180 genes were identified show significant differential expression. The most upregulated gene was SPP1 (Secreted Phosphoprotein 1), the protein-coding gene for Osteopontin (OPN), a structural protein of bone with extracellular matrix binding and cytokine activities. OPN exerts a regulatory effect in the differentiation of osteoblasts and enhances the metastatic potential of OSA tumoral cells at high levels by activating metalloproteinases ( 28 ). OPN also promotes oncogenesis and metastasis by inhibiting apoptosis and stimulating the endurance of tumoral cells and neovascularization ( 29 – 31 ). OPN modulates the immune system and possesses multiple roles in the inflammatory process ( 29 , 32 ). It serves as a chemotactic molecule to promote the migration of inflammatory cells and regulates dendritic cell responses ( 29 , 33 ). OPN is currently being studied as a potential cancer biomarker and target for the treatment of cancer ( 28 ). CXCL10 encodes for a pro-inflammatory cytokine that binds to the CXCR3 receptor predominantly expressed on immune cells (e.g., lymphocytes, macrophages, dendritic cells). This protein binding results in pleiotropic effects, including alterations of the tumor microenvironment and cytokine storms, as documented in SARS Co V-2 infections ( 34 , 35 ). The overexpression of CXCL10 has been associated with both antitumoral attributes and disease severity (hematologic and solid) depending on the type of cancer ( 36 , 37 ). CXCR4 oncogene encodes a chemokine receptor, which activation transduces a signal resulting in modulation of the AKT signaling cascade and MAPK1/MAPK3 activation. High expression of CXCR4 correlated with the evidence of pulmonary metastasis in human OSA; additionally, the CXCR4 ligand, CXCL12, promotes metastasis, angiogenesis, and growth of OSA cells ( 38 – 40 ). Fan et al. has shown that the inhibition of CXCR4 by zoledronate alters the metastatic behavior in canine OSA ( 15 ). Some of the activities of the proteins encoded by the genes mentioned above (i.e., SPP1, CXCL10, CXCR4) could enable or partially explain some of the clinical findings (i.e., pulmonary metastasis) observed in our case. The most downregulated gene was PTHLH, which encodes a protein involved in the regulation of endochondral bone development, as well as epithelial-mesenchymal interactions during the formation of mammary glands and teeth ( 41 – 43 ). The overexpression of its receptor PTHR1 in OSA has been associated with increased invasion and proliferation and conversely, decreased mRNA expression of PTHR1 has been associated with inhibition of proliferation, migration, and invasion in human OSA cell lines ( 44 , 45 ). Patients with strong staining for PTHR1 in canine OSA tumors had reduced survival times compared to those with weak immunostaining intensity( 45 ). Another downregulated gene was THBS1, a gene that encodes for the TSP1 protein, a glycoprotein that mediates cell-to-cell and cell-to-matrix interactions. This protein promotes cell migration and metastasis via FAK pathway and is highly upregulated in lung metastatic OSA ( 46 ). Similar associations with tumorigenesis, metastasis, and endothelial to mesenchymal transition have been described for the overexpression of LGALS3 protein and monocyte chemotactic protein 3( also known as CCL7) that encodes and is widely expressed in multiple cell types ( 47 , 48 ). We suspect that the downregulation of those genes (i.e, PTHLH, THSB1, LGALS3, CCL7) could be a cause or an effect on the OSA & HO pathogenesis. For example, Mainetti et al, showed that the KIT upregulation and activation stimulated cell migration by suppressing BRCA2 in human prostate carcinoma cells ( 49 ). In our case, we observed the upregulation of KIT but not the downregulation of BRCA2. The KIT – BRCA2 regulatory interaction serves as an example of the complexity of cellular pathway interactions that could explain differences in gene expression profiles noted above. When comparing the gene expression profile between primary and metastatic OSA lesions, we found the upregulation of CCL17 and VEGFC in primary lesion compared to metastasis. CCL17 and VEGFC proteins function as a cytokine and growth factor, respectively ( 50 , 51 ). Regarding the upregulation of CCL17 in the primary lesion, mature regulatory dendritic cells (mreg DCs) are tumor-specific and express CCL17 in human OSA ( 52 ). Liu et al. showed mreg DCs were enriched in OSA but absent in PBMCs, suggesting that mreg DCs are tumor-associated and may be immunosuppressive in the tumor microenvironment ( 52 ). Hence, the upregulation of CCL17 in the primary OSA may reflect a similar tumor microenvironment to human OSA, where it has enriched mreg DCs expressing CCL17. Regarding the VEGFC upregulation found in the primary lesion, the study by Park et al. demonstrated that OSA cell lines expressed 2–5-fold higher amplification of VEGFC mRNA than the control, and the higher VEGFC expression was associated with high-grade histologic type.( 53 ) In addition, a study by Yang et al. found that a higher VEGFC expression correlated with a higher tumor grade in chondrosarcoma ( 54 ). Canine OSA is generally deemed a high-grade tumor, and the presence of a high mitotic index (77 in 10 high-powered fields) in our case is concordant with those findings. FOS, C3, and C4BPA genes had significantly lower expression in primary lesion relative to metastasis. FOS is a proto-oncogene in humans; the overexpression of endogenous c-fos results in the development of OSA in transgenic mice ( 55 ). Also, C4BPA is a protein-coding gene that promotes cancer cell proliferation in CD40-expressing pancreas cancer ( 56 ). Liu et al. reported that C4BPA and FOS were significantly upregulated in metastatic human OSA compared to non-metastatic OSA ( 57 ). Similar to the findings of Liu et al., in our data, downregulation of FOS and C4BPA in primary lesion was found compared to metastasis. Thus, the upregulation of FOS and C4BPA could be a gene expression profile associated with metastatic lesions in canine OSA. The comparison of GSA between primary and metastatic lesions revealed that the upregulation of C4BPA and C3 genes in the metastatic lesions mainly contributed to the predominance of the Complement and Complement System. C4BPA is an isoform of acute phase protein C4BP, and during inflammation, circulating C4BP levels may increase up to 400% mainly driven by the increase in isoform C4BPA level ( 58 , 59 ). C3 is also associated with inflammatory responses; it was upregulated at mRNA level in human endometriosis cells compared to non-inflammatory healthy control ( 60 ). In addition, patients with HO had elevated complement levels, including C3 and C4 ( 61 ). Hence, activating the complement pathways through the upregulation of C3 and C4BPA may explain the inflammatory process observed in HO, such as pyrexia or pain in the extremities. The five differentially expressed genes identified here provide clues for unraveling the pathophysiology of HO and pulmonary metastasis in canine OSA. Specifically, we found the upregulation of CCL17, VEGFC, and the downregulation of FOS, C3, and C4BPA in the primary lesion compared to the metastatic lesion in the same dog with HO caused by solitary metastasis of OSA. Further studies with more cases are needed to systematically explore the significance of these changes and explore gene expression more extensively in HO cases. The IO panel consists of 780 genes most of which are involved in tumorigenesis and as well as some genes reflecting immune pathways and immune cells and is thus not as comprehensive a look into the transcriptome as bulk RNA sequencing. However, it has great sensitivity and specificity, even when formalin fixed, and paraffin embedded samples are used and is therefore very useful for analysis of clinical samples. Our data presented here may form the basis of further comparisons if additional cases are acquired. Abbreviations CnOb Canine osteoblast COL1A1 Collagen type I alpha chain CT Computed tomography FFPE Formalin-fixed paraffin-embedded GSS Global Significance Score HO Hypertrophic osteopathy MSU Michigan State University OPN Osteopontin OSA Osteosarcoma PDGF Platelet-derived growth factor SPP1 Secreted Phosphoprotein 1 VEGF Vascular endothelial growth factor Declarations Ethics approval and consent to participate This study followed the guidelines approved by the VTH -MSU Ethics Committee. The dog in this study was examined and treated with the owner's written consent. Consent for publication The owner signed an informed consent form including the authorization to publish the data arising from this study. Availability of data and materials Data is available from the authors upon reasonable request and with permission of NanoString. Competing interests The authors declare no potential conflicts of interest with respect to the research, authorship, and/or publication of this article. Funding The authors received no external financial support for the research, authorship, and/or publication of this article. Authors' contributions P.V.S. and V.Y.G conceived the original idea and design. K.K., J.D., and P.V.S. were involved in clinical data acquisitions and interpretations. K.K., A.I.E., Y.-T. Y, P.V.S, and V.Y.-G. were involved in the production, acquisition, and interpretation of transcriptomic data. K.K., J. D, Y.-T. Y, A.I.E., V.Y.-G., and P.V.S. were involved in writing, reviewing, and editing the manuscript. P.V.S. and V.Y.-G. contributed to the final critical revision of the manuscript on their experience on the topic. All authors read and approved the final manuscript. Acknowledgements The authors thank the Genomics Core of the Research and Technology Support Facility of Michigan State University for the NanoString assay. Authors' information 1 Department of Small Animal Clinical Science, College of Veterinary Medicine, University of Florida, Gainesville, FL 32608 2 Veterinary Emergency Group, Anaheim, CA 92808 3 Department of Small Animal Clinical Sciences, College of Veterinary Medicine, Michigan State University, East Lansing, MI 48824, USA. 4 Department of Microbiology, Genetics and Immunology, College of Veterinary Medicine, Michigan State University, East Lansing, MI 48824 References Brodey RS. Hypertrophic osteoarthropathy in the dog: a clinicopathologic survey of 60 cases. J Am Vet Med Assoc. 1971;159(10):1242-56. Withers SS, Johnson EG, Culp WT, Rodriguez CO, Jr., Skorupski KA, Rebhun RB. Paraneoplastic hypertrophic osteopathy in 30 dogs. Vet Comp Oncol. 2015;13(3):157-65. Pachamé AV, Beltrán MJ, Recchiuti NE, Portiansky EL, Massone AR, Gimeno EJ. Lesions of hypertrophic osteopathy in the forelimbs of a dog associated with pulmonary metastasis from a periosteal osteosarcoma. 2017. Ramoo S. Hypertrophic osteopathy associated with two pulmonary tumours and myocardial metastases in a dog: a case report. New Zealand Veterinary Journal. 2013;61(1):45-8. Liptak JM, Monnet E, Dernell WS, Withrow SJ. Pulmonary metastatectomy in the management of four dogs with hypertrophic osteopathy. Vet Comp Oncol. 2004;2(1):1-12. Turner H, Séguin B, Worley DR, Ehrhart NP, Lafferty MH, Withrow SJ, et al. Prognosis for dogs with stage III osteosarcoma following treatment with amputation and chemotherapy with and without metastasectomy. J Am Vet Med Assoc. 2017;251(11):1293-305. Treasure T. Hypertrophic pulmonary osteoarthropathy and the vagus nerve: an historical note. J R Soc Med. 2006;99(8):388-90. Nomori H, Kobayashi R, Kubo A, Morinaga S, Shintani Y, Sano T. Lung cancer containing growth hormone-releasing hormone associated with hypertrophic osteoarthropathy. Scandinavian journal of thoracic and cardiovascular surgery. 1994;28(3-4):149-52. Mito K, Maruyama R, Uenishi Y, Arita K, Kawano H, Kashima K, et al. Hypertrophic Pulmonary Osteoarthropathy Associated with Non-small Cell Lung Cancer Demonstrated Growth Hormone-Releasing Hormone by Immunohistochemical Analysis. Internal Medicine. 2001;40(6):532-5. Ooi A, Saad RA, Moorjani N, Amer KM. Effective symptomatic relief of hypertrophic pulmonary osteoarthropathy by video-assisted thoracic surgery truncal vagotomy. Ann Thorac Surg. 2007;83(2):684-5. Fenger JM, London CA, Kisseberth WC. Canine osteosarcoma: a naturally occurring disease to inform pediatric oncology. ILAR J. 2014;55(1):69-85. Cam M, Gardner HL, Roberts RD, Fenger JM, Guttridge DC, London CA, et al. ΔNp63 mediates cellular survival and metastasis in canine osteosarcoma. Oncotarget. 2016;7(30):48533-46. Fossey SL, Liao AT, McCleese JK, Bear MD, Lin J, Li PK, et al. Characterization of STAT3 activation and expression in canine and human osteosarcoma. BMC Cancer. 2009;9:81. Portela RF, Fadl-Alla BA, Pondenis HC, Byrum ML, Garrett LD, Wycislo KL, et al. Pro-tumorigenic effects of transforming growth factor beta 1 in canine osteosarcoma. J Vet Intern Med. 2014;28(3):894-904. Byrum ML, Pondenis HC, Fredrickson RL, Wycislo KL, Fan TM. Downregulation of CXCR4 Expression and Functionality After Zoledronate Exposure in Canine Osteosarcoma. J Vet Intern Med. 2016;30(4):1187-96. Roy J, Wycislo KL, Pondenis H, Fan TM, Das A. Comparative proteomic investigation of metastatic and non-metastatic osteosarcoma cells of human and canine origin. PLoS One. 2017;12(9):e0183930. Yang Y-T, Engleberg AI, Yuzbasiyan-Gurkan V. Establishment and Characterization of Cell Lines from Canine Metastatic Osteosarcoma. Cells. 2024;13(1):25. Danaher P, Warren S, Dennis L, D'Amico L, White A, Disis ML, et al. Gene expression markers of Tumor Infiltrating Leukocytes. J Immunother Cancer. 2017;5:18. Alexa A, Rahnenfuhrer J. topGO: enrichment analysis for gene ontology. R package version. 2010;2(0):2010. Mitchell AL, Attwood TK, Babbitt PC, Blum M, Bork P, Bridge A, et al. InterPro in 2019: improving coverage, classification and access to protein sequence annotations. Nucleic Acids Res. 2019;47(D1):D351-d60. Geer LY, Marchler-Bauer A, Geer RC, Han L, He J, He S, et al. The NCBI BioSystems database. Nucleic Acids Res. 2010;38(Database issue):D492-6. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545-50. Liberzon A, Subramanian A, Pinchback R, Thorvaldsdóttir H, Tamayo P, Mesirov JP. Molecular signatures database (MSigDB) 3.0. Bioinformatics. 2011;27(12):1739-40. Fabregat A, Jupe S, Matthews L, Sidiropoulos K, Gillespie M, Garapati P, et al. The Reactome Pathway Knowledgebase. Nucleic Acids Res. 2018;46(D1):D649-d55. Slenter DN, Kutmon M, Hanspers K, Riutta A, Windsor J, Nunes N, et al. WikiPathways: a multifaceted pathway database bridging metabolomics to other omics research. Nucleic Acids Res. 2018;46(D1):D661-d7. Malkov VA, Serikawa KA, Balantac N, Watters J, Geiss G, Mashadi-Hossein A, et al. Multiplexed measurements of gene signatures in different analytes using the Nanostring nCounter Assay System. BMC Res Notes. 2009;2:80. Geiss GK, Bumgarner RE, Birditt B, Dahl T, Dowidar N, Dunaway DL, et al. Direct multiplexed measurement of gene expression with color-coded probe pairs. Nat Biotechnol. 2008;26(3):317-25. Han X, Wang W, He J, Jiang L, Li X. Osteopontin as a biomarker for osteosarcoma therapy and prognosis. Oncol Lett. 2019;17(3):2592-8. Lund SA, Giachelli CM, Scatena M. The role of osteopontin in inflammatory processes. J Cell Commun Signal. 2009;3(3-4):311-22. Wai PY, Kuo PC. The role of Osteopontin in tumor metastasis. J Surg Res. 2004;121(2):228-41. Zhao J, Dong L, Lu B, Wu G, Xu D, Chen J, et al. Down-regulation of osteopontin suppresses growth and metastasis of hepatocellular carcinoma via induction of apoptosis. Gastroenterology. 2008;135(3):956-68. Mori R, Shaw TJ, Martin P. Molecular mechanisms linking wound inflammation and fibrosis: knockdown of osteopontin leads to rapid repair and reduced scarring. J Exp Med. 2008;205(1):43-51. Ashkar S, Weber GF, Panoutsakopoulou V, Sanchirico ME, Jansson M, Zawaideh S, et al. Eta-1 (osteopontin): an early component of type-1 (cell-mediated) immunity. Science. 2000;287(5454):860-4. Gudowska-Sawczuk M, Mroczko B. What Is Currently Known about the Role of CXCL10 in SARS-CoV-2 Infection? Int J Mol Sci. 2022;23(7). Huang H, Zhou W, Chen R, Xiang B, Zhou S, Lan L. CXCL10 is a Tumor Microenvironment and Immune Infiltration Related Prognostic Biomarker in Pancreatic Adenocarcinoma. Frontiers in Molecular Biosciences. 2021;8. Ariyarathna H, Thomson N, Aberdein D, Munday JS. Chemokine gene expression influences metastasis and survival time of female dogs with mammary carcinoma. Vet Immunol Immunopathol. 2020;227:110075. Russo E, Santoni A, Bernardini G. Tumor inhibition or tumor promotion? The duplicity of CXCR3 in cancer. Journal of Leukocyte Biology. 2020;108(2):673-85. Bianchi ME, Mezzapelle R. The Chemokine Receptor CXCR4 in Cell Proliferation and Tissue Regeneration. Front Immunol. 2020;11:2109. Cao Y, Hunter ZR, Liu X, Xu L, Yang G, Chen J, et al. The WHIM-like CXCR4(S338X) somatic mutation activates AKT and ERK, and promotes resistance to ibrutinib and other agents used in the treatment of Waldenstrom's Macroglobulinemia. Leukemia. 2015;29(1):169-76. Laverdiere C, Hoang BH, Yang R, Sowers R, Qin J, Meyers PA, et al. Messenger RNA Expression Levels of CXCR4 Correlate with Metastatic Behavior and Outcome in Patients with Osteosarcoma. Clinical Cancer Research. 2005;11(7):2561-7. Wysolmerski JJ, Philbrick WM, Dunbar ME, Lanske B, Kronenberg H, Broadus AE. Rescue of the parathyroid hormone-related protein knockout mouse demonstrates that parathyroid hormone-related protein is essential for mammary gland development. Development. 1998;125(7):1285-94. Philbrick WM, Dreyer BE, Nakchbandi IA, Karaplis AC. Parathyroid hormone-related protein is required for tooth eruption. Proc Natl Acad Sci U S A. 1998;95(20):11846-51. Foley J, Longely BJ, Wysolmerski JJ, Dreyer BE, Broadus AE, Philbrick WM. PTHrP regulates epidermal differentiation in adult mice. J Invest Dermatol. 1998;111(6):1122-8. Yang R, Hoang BH, Kubo T, Kawano H, Chou A, Sowers R, et al. Over-expression of parathyroid hormone Type 1 receptor confers an aggressive phenotype in osteosarcoma. International Journal of Cancer. 2007;121(5):943-54. Al-Khan AA, Al Balushi NR, Richardson SJ, Danks JA. Roles of Parathyroid Hormone-Related Protein (PTHrP) and Its Receptor (PTHR1) in Normal and Tumor Tissues: Focus on Their Roles in Osteosarcoma. Front Vet Sci. 2021;8:637614. Hu C, Wen J, Gong L, Chen X, Wang J, Hu F, et al. Thrombospondin-1 promotes cell migration, invasion and lung metastasis of osteosarcoma through FAK dependent pathway. Oncotarget. 2017;8(44):75881-92. Liu Y, Cai Y, Liu L, Wu Y, Xiong X. Crucial biological functions of CCL7 in cancer. PeerJ. 2018;6:e4928. Ahmed H, AlSadek DM. Galectin-3 as a Potential Target to Prevent Cancer Metastasis. Clin Med Insights Oncol. 2015;9:113-21. Mainetti LE, Zhe X, Diedrich J, Saliganan AD, Cho WJ, Cher ML, et al. Bone-induced c-kit expression in prostate cancer: a driver of intraosseous tumor growth. Int J Cancer. 2015;136(1):11-20. Scheu S, Ali S, Ruland C, Arolt V, Alferink J. The C-C Chemokines CCL17 and CCL22 and Their Receptor CCR4 in CNS Autoimmunity. Int J Mol Sci. 2017;18(11). Rauniyar K, Jha SK, Jeltsch M. Biology of Vascular Endothelial Growth Factor C in the Morphogenesis of Lymphatic Vessels. Front Bioeng Biotechnol. 2018;6:7. Liu W, Hu H, Shao Z, Lv X, Zhang Z, Deng X, et al. Characterizing the tumor microenvironment at the single-cell level reveals a novel immune evasion mechanism in osteosarcoma. Bone Research. 2023;11(1):4. Park HR, Min K, Kim HS, Jung WW, Park YK. Expression of vascular endothelial growth factor-C and its receptor in osteosarcomas. Pathol Res Pract. 2008;204(8):575-82. Yang W-H, Chang A-C, Wang S-W, Wang S-J, Chang Y-S, Chang T-M, et al. Leptin promotes VEGF-C production and induces lymphangiogenesis by suppressing miR-27b in human chondrosarcoma cells. Scientific Reports. 2016;6(1):28647. Wang Z-Q, Liang J, Schellander K, Wagner EF, Grigoriadis AE. c-fos-induced Osteosarcoma Formation in Transgenic Mice: Cooperativity with c-jun and the Role of Endogenous c-fos1. Cancer Research. 1995;55(24):6244-51. Sasaki K, Takano S, Tomizawa S, Miyahara Y, Furukawa K, Takayashiki T, et al. C4b-binding protein α-chain enhances antitumor immunity by facilitating the accumulation of tumor-infiltrating lymphocytes in the tumor microenvironment in pancreatic cancer. Journal of Experimental & Clinical Cancer Research. 2021;40(1):212. Liu K, He Q, Liao G, Han J. Identification of critical genes and gene interaction networks that mediate osteosarcoma metastasis to the lungs. Exp Ther Med. 2015;10(5):1796-806. Criado García O, Sánchez-Corral P, Rodríguez de Córdoba S. Isoforms of human C4b-binding protein. II. Differential modulation of the C4BPA and C4BPB genes by acute phase cytokines. J Immunol. 1995;155(8):4037-43. García de Frutos P, Alim RI, Härdig Y, Zöller B, Dahlbäck B. Differential regulation of alpha and beta chains of C4b-binding protein during acute-phase response resulting in stable plasma levels of free anticoagulant protein S. Blood. 1994;84(3):815-22. Agostinis C, Zorzet S, Balduit A, Zito G, Mangogna A, Macor P, et al. The Inflammatory Feed-Forward Loop Triggered by the Complement Component C3 as a Potential Target in Endometriosis. Frontiers in Immunology. 2021;12. Pisarcik D, Kousourou G, Vinod N. A Case of Paraneoplastic Arthralgias. American Journal of Medical Case Reports. 2021;9(8):397-401. Additional Declarations No competing interests reported. Supplementary Files FINALsupplementalmaterialsOSAKK05172024.docx Cite Share Download PDF Status: Under Review Version 1 posted Editorial decision: Revision requested 17 Jul, 2024 Reviews received at journal 17 Jul, 2024 Reviewers agreed at journal 12 Jul, 2024 Reviewers agreed at journal 09 Jul, 2024 Reviews received at journal 11 Jun, 2024 Reviewers agreed at journal 22 May, 2024 Reviewers agreed at journal 21 May, 2024 Reviewers invited by journal 20 May, 2024 Submission checks completed at journal 19 May, 2024 Editor assigned by journal 19 May, 2024 First submitted to journal 17 May, 2024 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4438777","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":305595932,"identity":"9d891192-c6fe-412b-a615-8952e21c8dfa","order_by":0,"name":"Keita Kitagawa","email":"","orcid":"","institution":"University of Florida","correspondingAuthor":false,"prefix":"","firstName":"Keita","middleName":"","lastName":"Kitagawa","suffix":""},{"id":305595935,"identity":"693a880d-4dbe-4053-83a8-d22f0ce723bf","order_by":1,"name":"Jessica Dryfhout","email":"","orcid":"","institution":"Veterinary Emergency Group","correspondingAuthor":false,"prefix":"","firstName":"Jessica","middleName":"","lastName":"Dryfhout","suffix":""},{"id":305595937,"identity":"f9c35e61-a534-488f-9201-dea12b8b3a39","order_by":2,"name":"Alexander I. Engleberg","email":"","orcid":"","institution":"Michigan State University","correspondingAuthor":false,"prefix":"","firstName":"Alexander","middleName":"I.","lastName":"Engleberg","suffix":""},{"id":305595940,"identity":"6ac99ede-7c93-4683-bcc6-7c87c993e3f1","order_by":3,"name":"Ya-Ting Yang","email":"","orcid":"","institution":"Michigan State University","correspondingAuthor":false,"prefix":"","firstName":"Ya-Ting","middleName":"","lastName":"Yang","suffix":""},{"id":305595942,"identity":"7825e792-25bb-4d10-b04a-e7a65191f663","order_by":4,"name":"Vilma Yuzbasiyan-Gurkan","email":"","orcid":"","institution":"Michigan State University","correspondingAuthor":false,"prefix":"","firstName":"Vilma","middleName":"","lastName":"Yuzbasiyan-Gurkan","suffix":""},{"id":305595946,"identity":"a2d4f0aa-1470-4865-8187-521c4b22fbb7","order_by":5,"name":"Paulo Vilar Saavedra","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA+ElEQVRIiWNgGAWjYFACxsYDMKYEA4MNA3s7kMWDX0sDspY0Bp7DBLUwMCBrOUxYi7lEcsOBn3vq5PlnJB+8XVFzPrGHmYHxwds23FosZyQ2HOx5dthwxo20ZMszx26DtDAbzsWjxeDMwYYDPAcOMDacOWMm2cB2O3E/MwObNC8BLQf/HKizn3/m/DfJhn/nQLaw/8ar5Xhjw2GeA8yJG473sEk2th0AaWFjxqfFsh2oRebA4eSNx9uMLRv7ko17mBmbJeecw63FnJn94cM3B+ps5x1mfniz4ZudbA9788EPb8rwOAyLGGMDbvU4tIyCUTAKRsEoQAUAX3VZ2Zo/woMAAAAASUVORK5CYII=","orcid":"","institution":"University of Florida","correspondingAuthor":true,"prefix":"","firstName":"Paulo","middleName":"Vilar","lastName":"Saavedra","suffix":""}],"badges":[],"createdAt":"2024-05-17 22:53:21","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-4438777/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4438777/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":57728011,"identity":"b35c4eca-aa57-4fd2-aba2-bee24b4cc0df","added_by":"auto","created_at":"2024-06-04 21:41:25","extension":"jpeg","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":356453,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eDifferential gene expression between primary and metastatic OSA and osteoblasts\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003ea) Heatmap of the gene expressions in canine osteoblast cells (in orange bar on top) and OSA cells (in blue bar on top).\u003c/p\u003e\n\u003cp\u003eb) Volcano plot of differential gene expressions in OSA cells relative to canine osteoblast cells. Green dots mean the gene upregulated and purple dots mean the gene downregulated.\u003c/p\u003e","description":"","filename":"floatimage1.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-4438777/v1/cf48aaee85a8b3eb1a9a4ed5.jpeg"},{"id":57728013,"identity":"90e7855e-8607-47b8-957c-0af7f85eb830","added_by":"auto","created_at":"2024-06-04 21:41:26","extension":"jpeg","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":165756,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eDifferential gene expression between primary and metastatic OSA.\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003ea) The heatmap comparison of the gene expressions between metastatic OSA (in orange bar on top) and primary OSA (in blue bar on top)\u003c/p\u003e\n\u003cp\u003eb) Volcano plot of differential gene expressions in primary OSA relative to metastatic OSA. Green dots indicate upregulated and purple dots indicate downregulated genes in primary OSA.\u003c/p\u003e","description":"","filename":"floatimage2.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-4438777/v1/92a730cd279fa1c62eaaba90.jpeg"},{"id":57728803,"identity":"cc4a8692-2d1d-4fa5-a719-4135bc4f6790","added_by":"auto","created_at":"2024-06-04 21:49:30","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1009489,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4438777/v1/b2cfd589-dbb2-48df-9179-83458fc93842.pdf"},{"id":57728799,"identity":"2b2f773b-59f1-480b-85cc-35cccb7b4001","added_by":"auto","created_at":"2024-06-04 21:49:25","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":49708,"visible":true,"origin":"","legend":"","description":"","filename":"FINALsupplementalmaterialsOSAKK05172024.docx","url":"https://assets-eu.researchsquare.com/files/rs-4438777/v1/c42d353a045a105b2e824d8e.docx"}],"financialInterests":"No competing interests reported.","formattedTitle":"Gene expression profile comparison of primary and pulmonary metastatic lesions in a dog with appendicular osteosarcoma and hypertrophic osteopathy","fulltext":[{"header":"BACKGROUND","content":"\u003cp\u003eCanine hypertrophic osteopathy (HO) is a paraneoplastic syndrome associated with pulmonary tumoral lesions and is characterized by the symmetrical formation of periosteal new bone along the diaphysis of the long bones of the appendicular skeleton(\u003cspan additionalcitationids=\"CR2 CR3\" citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e). Clinical signs include lameness, lethargy, localized pain, and soft tissue swelling surrounding the affected bones(\u003cspan additionalcitationids=\"CR2 CR3\" citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e). The diagnosis of HO is made based on the clinical history, clinical signs, and radiographical changes(\u003cspan additionalcitationids=\"CR3\" citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e). Metastasectomy, a surgical intervention to remove one or more metastatic lesions, can result in the control of HO by achieving the resolution of the clinical signs and the radiographical changes.(\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e, \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e, \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e). The pathogenesis of HO is not currently fully understood. The veterinary literature in this regard is null, and the scarce references in human literature propose two major hypotheses explaining the development of HO: neurological and humoral pathways (\u003cspan additionalcitationids=\"CR8\" citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e). Resolution of symptoms following vagotomy empirically suggests neurologic involvement in HO, whereas humoral pathway involvement is supported by increased release of serum growth factors, including growth hormone, growth hormone-releasing hormone, platelet-derived growth factor (PDGF), and vascular endothelial growth factor (VEGF) by the metastatic lesion (\u003cspan additionalcitationids=\"CR9\" citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eCanine appendicular osteosarcoma (OSA) is the most common form of malignant bone cancer in dogs, which is locally aggressive and highly metastatic; OSA is also the most notable cause of paraneoplastic HO in dogs (\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e, \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e). Regarding the pathophysiology of canine OSA, many molecular and genetic factors have been investigated to explain the development of metastasis (\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e, \u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e). The overexpression of delta Np63, a transcription factor belonging to p53 family, has been reported to be correlated with pulmonary metastases via increased secretion of VEGFA through activation of STAT3(\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e). In addition, the downregulation of STAT3 decreases the expression of VEGF, and conversely, activation of STAT3 increases VEGF expression, and may thus facilitate metastasis in OSA (\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e). Also, in addition to STAT3 and delta Np63, TGF beta also regulates VEGF secretion. The inhibition of TGF beta1 decreases VEGF secretion, suggesting TGF beta1 is associated with angiogenesis via VEGF(\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e). TGF beta signaling is suggested to promote distant metastasis by activating cell migration and neoangiogenesis in OSA (\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e). Furthermore, the downregulation of CXCR4 by zoledronate treatment alters the pattern of metastasis in OSA, curtailing metastasis to the lungs, but, not to other organs, such as the nerves and skin, suggesting CXCR4 has some contribution to pulmonary metastasis in OSA (\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eMolecular investigations in metastasis in canine OSA until now have focused on the pathophysiology of the primary OSA lesion. Lately, however, proteomics for OSA has shown differences in protein expression between primary and metastatic lesions (\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e). Identifying the differences in gene expression profiling between primary and metastatic lesions would provide us with critical information for further understanding the pathophysiology of HO, as well as the metastatic process in OSA. Hence, the purpose of this study was to compare the gene expression profiles between primary and metastatic OSA lesions from the same dog who developed HO due to metastatic OSA.\u003c/p\u003e"},{"header":"MATERIALS AND METHODS","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003ePRIMARY AND METASTATIC OSA SAMPLES\u003c/h2\u003e \u003cp\u003ePrimary and metastatic OSA samples were collected from a 2-year-old, intact male, Labrador retriever that presented to the Michigan State University (MSU) Veterinary Medical Center. Both appendicular and solitary pulmonary metastatic lesions were surgically removed, and histopathologic diagnoses were confirmed by board-certified pathologists at MSU.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003eCELL CULTURE\u003c/h2\u003e \u003cp\u003eThree normal osteoblast strains were used in the current study as controls. Among them, one canine osteoblast strain (CnOb) was purchased from Cell Application Inc. (San Diego, CA, USA). Additional two canine osteoblast cell strains were derived from two different healthy beagles in our lab as described (\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e). All three CnOb cell strains were maintained with CnOb Growth Medium (Cell Application Inc., San Diego, CA, USA) and supplemented with 10% fetal bovine serum (Thermo Fisher Scientific, Waltham, MA, USA), antibiotics (0.1% gentamycin, Life Technologies, Carlsbad, CA, USA) and incubated in a humidified incubator at 37˚C with 5% CO\u003csub\u003e2\u003c/sub\u003e. These three canine osteoblasts strains were shown to express common osteoblast markers, including osteocalcin and collagen type I alpha chain (COL1A1) via quantitative real-time PCR previously by our team (\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec5\" class=\"Section2\"\u003e \u003ch2\u003eRNA ISOLATION\u003c/h2\u003e \u003cp\u003eFormalin-fixed paraffin-embedded (FFPE) canine OSA tissues from primary and metastatic lesions were obtained from the MSU Veterinary Diagnostic Library. Total RNA was isolated from FFPE tissue blocks using the RNeasy FFPE Kit (Qiagen, Hilden, Germany) according to manufacturer\u0026rsquo;s instructions. RNA was also isolated from canine osteoblasts using mirVana miRNA Isolation Kit with phenol (Thermo Fisher Scientific, Waltham, MA, USA). Each sample was processed in triplicate. Isolated RNA was immediately stored at -80˚C for further studies.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec6\" class=\"Section2\"\u003e \u003ch2\u003eDATA ANALYSIS BY NANOSTRING TECHNOLOGY\u003c/h2\u003e \u003cp\u003eGene expression was evaluated using the NanoString nCounter\u0026reg; Canine IO (immuno-oncology) gene-expression panel comprising of 780 genes, developed using the CanFam3.1 reference genome. Targeted expression profiling was performed using isolated RNA from OSA samples and canine osteoblast cell lines. Differential gene expression analysis was performed comparing affected OSA samples to non-tumoral control osteoblast cells. Further analysis compared expression levels between primary and metastatic OSA samples.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec7\" class=\"Section2\"\u003e \u003ch2\u003eSTATISTICAL ANALYSIS\u003c/h2\u003e \u003cp\u003eData normalization, quality control, and functional and pathway analysis were performed by importing nCounter data onto the ROSALIND\u0026reg; analysis platform (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.rosalind.bio\u003c/span\u003e\u003cspan address=\"https://www.rosalind.bio\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e). Normalization, fold changes, and p-values were calculated in-platform following the nCounter\u0026reg; advanced analysis protocol comparing the geometric mean of housekeeping probes to sample counts from the same lane of the array. Housekeeping genes used for normalization were chosen based on the geNorm algorithm as used in the NormqPCR R library.(\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e) In addition to differential gene expression, gene set analysis (GSA) and cell type profiling were also performed using the ROSALIND platform(\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e). Differential expression was determined as significant under the parameters of a fold-change \u0026gt;|1.5|and and an adjusted P-value\u0026thinsp;\u0026lt;\u0026thinsp;0.05, using the Benjamini-Hochberg correction method. The ROSALIND platform performs enrichment analysis by referencing several database sources, such as Interpro, NCBI, MSigDB, REACTOME, and WikiPathways (\u003cspan additionalcitationids=\"CR21 CR22 CR23 CR24\" citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e). Volcano plots and heatmaps were generated on ROSALIND comparing OSA cells to osteoblast cells and primary OSA to metastatic OSA.\u003c/p\u003e \u003c/div\u003e"},{"header":"RESULTS","content":"\u003cdiv id=\"Sec9\" class=\"Section2\"\u003e \u003ch2\u003eCASE HISTORY\u003c/h2\u003e \u003cp\u003eA 2-year-old, intact male, Labrador retriever was presented for a progressive right hindlimb lameness of 2-months duration at MSU Veterinary Medical Center. The dog was non-weight bearing and moderate muscle atrophy of the right hindlimb was observed. There was a large firm mass of approximately 10 cm in diameter in the right proximal tibia, causing a decreased range of motion of the right stifle, and pitting edema distal to the stifle joint. For staging purposes, we obtained bloodwork, thoracic radiographs, pre- and post-contrast Computed tomography (CT) scan. Staging results did not reveal any remarkable abnormalities except for an irregular shaped, large mass associated with the caudal aspect of the proximal tibia by CT scan. The histopathology report of surgical biopsy sample obtained from the proximal tibia mass confirmed the diagnosis of chondroblastic OSA with a mitotic index of 77 in 10 high powered fields. The dog then underwent a right hindlimb amputation; however, adjuvant chemotherapy was not pursued at that time per the owner\u0026rsquo;s decision.\u003c/p\u003e \u003cp\u003eThe dog was examined 314 days postoperatively because of lethargy and decreased appetite. The dog had progressive lameness, pyrexia of 40.4˚C, soft tissue swelling and pain in the extremities. Orthopedic examination revealed a decreased range of motion in both elbow and carpal joints and pain upon extension of elbows bilaterally. Blood work revealed inflammatory leukogram, mild-moderate elevated ALP [154 U/L (normal range: 10\u0026ndash;92 U/L)], and low urea nitrogen [8 mg/dL (normal range: 12\u0026ndash;27 mg/dL)], otherwise within reference intervals.\u003c/p\u003e \u003cp\u003eFor the diagnostic imaging, thoracic radiographs revealed the presence of a large intrathoracic mass and a smoothly marginated periosteal proliferation along the right humerus. Also, in CT scan, periosteal proliferations were noted along the fore and remaining hindlimb including the right humerus, the left tibia, the left metatarsal bones, and the left femur. The CT scan further indicated a solitary mass in the cranial mediastinum or right cranial lung lobe. Fine needle aspiration of the mediastinal mass was done; the result of cytology report was consistent with OSA. Integrating these findings, the dog was diagnosed as HO secondary to pulmonary metastatic OSA. For the treatment of HO, the cranial mediastinal mass was surgically removed, and submitted for histopathology. The histopathology evaluation was consistent with chondroblastic OSA, likely metastatic. As adjuvant therapy following metastasectomy, carboplatin (250-300mg/m\u003csup\u003e2\u003c/sup\u003e IV every 3 weeks) was given a total of 6 doses and without remarkable adverse effects throughout the chemotherapy.\u003c/p\u003e \u003cp\u003ePostoperative follow-up was done by physical examination with either radiographs or CT scan. Clinical improvement of HO was achieved 9 days after metastasectomy; the dog regained mobility without associated pain. The first postoperative CT scan, performed 58 days after metastasectomy, showed that the dog retained radiologic evidence of HO. The previously noted HO lesions persisted and most of them were slightly more severe in right thoracic limb involving the right scapula, right humerus, right metacarpal bones, and the bilateral radii, ulnae, and second metacarpal bones. Interestingly, the mild improvement of HO in the right scapula was found radiologically based on the CT scan at that time. No lameness nor pain on the legs were appreciated at the second postoperative CT scan (i.e., 128 days after metastasectomy). Palpation revealed no abnormal findings and the previously noted periosteal proliferation (i.e., HO lesion) confirmed no longer present in the right humerus by CT scan.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec10\" class=\"Section2\"\u003e \u003ch2\u003eGENE EXPRESSION PROFILE ANALYSIS\u003c/h2\u003e \u003cp\u003eWe used the NanoString nCounter\u0026reg; Canine IO panel which has great sensitivity and specificity to elucidate the differentially expressed genes in OSA cells compared to non-tumoral osteoblasts(\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e, \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e). The overview of the gene expression data is presented in the heatmap in Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e. Distinct expression profiles for both up- and downregulated genes across all samples can be appreciated in Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ea. The volcano plot shows all samples plotted as a functional fold change vs. p-value (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eb and \u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eb). Genes that exhibited a significant (p\u0026thinsp;\u0026lt;\u0026thinsp;0.05) and 1.5-fold change in expression compared to the control group were selected. Overall, 180 genes were differentially regulated; 121 of these were upregulated, whereas 59 were downregulated in OSA compared to control (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eb). Some of the top upregulated include genes involved in the immune response and cytokine production, such as SPP1, CXCL10, CXCR4, whereas other genes were associated with the cell signaling pathway, such as TNFRSF11A and KIT (Table \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e). On the other hand, of the 59 downregulated genes, the most down-regulated genes included those involved in cell migration and the cell cycle, such as PTHR1, THBS1, CCND1, CCDN2, and CDKN1A among others (Table S2).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eOf special interest to us were the differences in gene expression between the primary and metastatic lesions. Although the number of differential genes expressed was small, the heatmap showed distinct expression profiles for both up- and downregulated genes. (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ea). Overall, 5 genes were differently expressed. The CCL17 and VEGFC were significantly upregulated in primary lesion when compared to metastatic lesion, while FOS, C3, and C4BPA were significantly downregulated in primary lesion when compared to metastatic lesion (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eb). Fold changes of these differential gene expressions are summarized in Table. 1.\u003c/p\u003e \u003cp\u003eThrough ROSALIND\u0026reg;, several differentially expressed signaling pathways were identified. Pathway significance was determined by a Global Significance Score (GSS) rating based on gene set analysis (Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). The GSS measures the extent of the difference of overall differential gene expression between primary and metastasis OSAs within a distinct signaling pathway. GSSs are positive regardless of upregulation or downregulation, with a large value indicating greater significance. The pathway with the highest GSS rating between primary and metastatic OSA was the Complement pathway with a score of 1.60. The other significant pathway was the Complement System pathway (1.41). The pathways with GSS ratings less than 1 included the NK cell functions pathway (0.92), the Angiogenesis pathway (0.85), the Chemokines pathway (0.77), the Wnt Signaling pathway (0.77), the Interleukins pathway (0.76), the Epigenetic Regulation pathway (0.76), the Regulation pathway (0.74) and the TGF-beta signaling pathway (0.74). The detailed information of these pathways was summarized in the Table S3.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eGenes with significant differential expression between the Primary vs Metastatic OS in this patient\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"3\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003e Gene\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eFold Change\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eLog\u003csub\u003e2\u003c/sub\u003e fold change\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eCCL17\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e+\u0026thinsp;17.5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e4.1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eVEGFC\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e+\u0026thinsp;3.7\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e1.9\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eFOS\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u0026minus;\u0026thinsp;11\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-3.5\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eC3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u0026minus;\u0026thinsp;34\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-5.1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eC4BPA\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e-48\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-5.6\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab2\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 2\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eGlobal significance ratings comparing overall differential expression of selected pathways between primary vs metastatic OSA.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"2\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eTerm Name\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eGlobal Significance Score\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eComplement\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e1.6066\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eComplement System\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e1.415\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eNK Cell Functions\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.9216\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eAngiogenesis\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.8553\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eChemokines\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.774\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eWnt Signaling\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.7688\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eInterleukins\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.7613\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eEpigenetic Regulation\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.7612\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eRegulation\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.7441\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eTGF-beta Signaling\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.7385\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003c/div\u003e"},{"header":"DISCUSSION","content":"\u003cp\u003eThe clinical history of our OSA-HO case, including blood work results, radiographic findings, and the radiological resolution of HO after metastasectomy, are consistent with the previous studies (\u003cspan additionalcitationids=\"CR3 CR4\" citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e). Therefore, it is an archetypical subject of analysis for further understanding the pathophysiology of HO and OSA in dogs. When comparing the gene expression profile of OSA vs. control, 180 genes were identified show significant differential expression. The most upregulated gene was SPP1 (Secreted Phosphoprotein 1), the protein-coding gene for Osteopontin (OPN), a structural protein of bone with extracellular matrix binding and cytokine activities. OPN exerts a regulatory effect in the differentiation of osteoblasts and enhances the metastatic potential of OSA tumoral cells at high levels by activating metalloproteinases (\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e). OPN also promotes oncogenesis and metastasis by inhibiting apoptosis and stimulating the endurance of tumoral cells and neovascularization (\u003cspan additionalcitationids=\"CR30\" citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e). OPN modulates the immune system and possesses multiple roles in the inflammatory process (\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e, \u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e). It serves as a chemotactic molecule to promote the migration of inflammatory cells and regulates dendritic cell responses (\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e, \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e). OPN is currently being studied as a potential cancer biomarker and target for the treatment of cancer (\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e). CXCL10 encodes for a pro-inflammatory cytokine that binds to the CXCR3 receptor predominantly expressed on immune cells (e.g., lymphocytes, macrophages, dendritic cells). This protein binding results in pleiotropic effects, including alterations of the tumor microenvironment and cytokine storms, as documented in SARS Co V-2 infections (\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e, \u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e). The overexpression of CXCL10 has been associated with both antitumoral attributes and disease severity (hematologic and solid) depending on the type of cancer (\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e, \u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e). CXCR4 oncogene encodes a chemokine receptor, which activation transduces a signal resulting in modulation of the AKT signaling cascade and MAPK1/MAPK3 activation. High expression of CXCR4 correlated with the evidence of pulmonary metastasis in human OSA; additionally, the CXCR4 ligand, CXCL12, promotes metastasis, angiogenesis, and growth of OSA cells (\u003cspan additionalcitationids=\"CR39\" citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e). Fan et al. has shown that the inhibition of CXCR4 by zoledronate alters the metastatic behavior in canine OSA (\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e). Some of the activities of the proteins encoded by the genes mentioned above (i.e., SPP1, CXCL10, CXCR4) could enable or partially explain some of the clinical findings (i.e., pulmonary metastasis) observed in our case.\u003c/p\u003e \u003cp\u003eThe most downregulated gene was PTHLH, which encodes a protein involved in the regulation of endochondral bone development, as well as epithelial-mesenchymal interactions during the formation of mammary glands and teeth (\u003cspan additionalcitationids=\"CR42\" citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e43\u003c/span\u003e). The overexpression of its receptor PTHR1 in OSA has been associated with increased invasion and proliferation and conversely, decreased mRNA expression of PTHR1 has been associated with inhibition of proliferation, migration, and invasion in human OSA cell lines (\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e, \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e). Patients with strong staining for PTHR1 in canine OSA tumors had reduced survival times compared to those with weak immunostaining intensity(\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e). Another downregulated gene was THBS1, a gene that encodes for the TSP1 protein, a glycoprotein that mediates cell-to-cell and cell-to-matrix interactions. This protein promotes cell migration and metastasis via FAK pathway and is highly upregulated in lung metastatic OSA (\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e). Similar associations with tumorigenesis, metastasis, and endothelial to mesenchymal transition have been described for the overexpression of LGALS3 protein and monocyte chemotactic protein 3( also known as CCL7) that encodes and is widely expressed in multiple cell types (\u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e, \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e). We suspect that the downregulation of those genes (i.e, PTHLH, THSB1, LGALS3, CCL7) could be a cause or an effect on the OSA \u0026amp; HO pathogenesis. For example, Mainetti et al, showed that the KIT upregulation and activation stimulated cell migration by suppressing BRCA2 in human prostate carcinoma cells (\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e). In our case, we observed the upregulation of KIT but not the downregulation of BRCA2. The KIT \u0026ndash; BRCA2 regulatory interaction serves as an example of the complexity of cellular pathway interactions that could explain differences in gene expression profiles noted above.\u003c/p\u003e \u003cp\u003eWhen comparing the gene expression profile between primary and metastatic OSA lesions, we found the upregulation of CCL17 and VEGFC in primary lesion compared to metastasis. CCL17 and VEGFC proteins function as a cytokine and growth factor, respectively (\u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e50\u003c/span\u003e, \u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e51\u003c/span\u003e). Regarding the upregulation of CCL17 in the primary lesion, mature regulatory dendritic cells (mreg DCs) are tumor-specific and express CCL17 in human OSA (\u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e52\u003c/span\u003e). Liu et al. showed mreg DCs were enriched in OSA but absent in PBMCs, suggesting that mreg DCs are tumor-associated and may be immunosuppressive in the tumor microenvironment (\u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e52\u003c/span\u003e). Hence, the upregulation of CCL17 in the primary OSA may reflect a similar tumor microenvironment to human OSA, where it has enriched mreg DCs expressing CCL17.\u003c/p\u003e \u003cp\u003eRegarding the VEGFC upregulation found in the primary lesion, the study by Park et al. demonstrated that OSA cell lines expressed 2\u0026ndash;5-fold higher amplification of VEGFC mRNA than the control, and the higher VEGFC expression was associated with high-grade histologic type.(\u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e53\u003c/span\u003e) In addition, a study by Yang et al. found that a higher VEGFC expression correlated with a higher tumor grade in chondrosarcoma (\u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e). Canine OSA is generally deemed a high-grade tumor, and the presence of a high mitotic index (77 in 10 high-powered fields) in our case is concordant with those findings.\u003c/p\u003e \u003cp\u003eFOS, C3, and C4BPA genes had significantly lower expression in primary lesion relative to metastasis. FOS is a proto-oncogene in humans; the overexpression of endogenous c-fos results in the development of OSA in transgenic mice (\u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e55\u003c/span\u003e). Also, C4BPA is a protein-coding gene that promotes cancer cell proliferation in CD40-expressing pancreas cancer (\u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e56\u003c/span\u003e). Liu et al. reported that C4BPA and FOS were significantly upregulated in metastatic human OSA compared to non-metastatic OSA (\u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e57\u003c/span\u003e). Similar to the findings of Liu et al., in our data, downregulation of FOS and C4BPA in primary lesion was found compared to metastasis. Thus, the upregulation of FOS and C4BPA could be a gene expression profile associated with metastatic lesions in canine OSA.\u003c/p\u003e \u003cp\u003eThe comparison of GSA between primary and metastatic lesions revealed that the upregulation of C4BPA and C3 genes in the metastatic lesions mainly contributed to the predominance of the Complement and Complement System. C4BPA is an isoform of acute phase protein C4BP, and during inflammation, circulating C4BP levels may increase up to 400% mainly driven by the increase in isoform C4BPA level (\u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e58\u003c/span\u003e, \u003cspan citationid=\"CR59\" class=\"CitationRef\"\u003e59\u003c/span\u003e). C3 is also associated with inflammatory responses; it was upregulated at mRNA level in human endometriosis cells compared to non-inflammatory healthy control (\u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003e60\u003c/span\u003e). In addition, patients with HO had elevated complement levels, including C3 and C4 (\u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e61\u003c/span\u003e). Hence, activating the complement pathways through the upregulation of C3 and C4BPA may explain the inflammatory process observed in HO, such as pyrexia or pain in the extremities.\u003c/p\u003e \u003cp\u003eThe five differentially expressed genes identified here provide clues for unraveling the pathophysiology of HO and pulmonary metastasis in canine OSA. Specifically, we found the upregulation of CCL17, VEGFC, and the downregulation of FOS, C3, and C4BPA in the primary lesion compared to the metastatic lesion in the same dog with HO caused by solitary metastasis of OSA. Further studies with more cases are needed to systematically explore the significance of these changes and explore gene expression more extensively in HO cases. The IO panel consists of 780 genes most of which are involved in tumorigenesis and as well as some genes reflecting immune pathways and immune cells and is thus not as comprehensive a look into the transcriptome as bulk RNA sequencing. However, it has great sensitivity and specificity, even when formalin fixed, and paraffin embedded samples are used and is therefore very useful for analysis of clinical samples. Our data presented here may form the basis of further comparisons if additional cases are acquired.\u003c/p\u003e"},{"header":"Abbreviations","content":"\u003cdiv class=\"DefinitionList\"\u003e \u003cdiv class=\"DefinitionListEntry\"\u003e \u003cdiv class=\"Term\"\u003eCnOb\u003c/div\u003e \u003cdiv class=\"Description\"\u003e \u003cp\u003eCanine osteoblast\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003cdiv class=\"DefinitionListEntry\"\u003e \u003cdiv class=\"Term\"\u003eCOL1A1\u003c/div\u003e \u003cdiv class=\"Description\"\u003e \u003cp\u003eCollagen type I alpha chain\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003cdiv class=\"DefinitionListEntry\"\u003e \u003cdiv class=\"Term\"\u003eCT\u003c/div\u003e \u003cdiv class=\"Description\"\u003e \u003cp\u003eComputed tomography\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003cdiv class=\"DefinitionListEntry\"\u003e \u003cdiv class=\"Term\"\u003eFFPE\u003c/div\u003e \u003cdiv class=\"Description\"\u003e \u003cp\u003eFormalin-fixed paraffin-embedded\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003cdiv class=\"DefinitionListEntry\"\u003e \u003cdiv class=\"Term\"\u003eGSS\u003c/div\u003e \u003cdiv class=\"Description\"\u003e \u003cp\u003eGlobal Significance Score\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003cdiv class=\"DefinitionListEntry\"\u003e \u003cdiv class=\"Term\"\u003eHO\u003c/div\u003e \u003cdiv class=\"Description\"\u003e \u003cp\u003eHypertrophic osteopathy\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003cdiv class=\"DefinitionListEntry\"\u003e \u003cdiv class=\"Term\"\u003eMSU\u003c/div\u003e \u003cdiv class=\"Description\"\u003e \u003cp\u003eMichigan State University\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003cdiv class=\"DefinitionListEntry\"\u003e \u003cdiv class=\"Term\"\u003eOPN\u003c/div\u003e \u003cdiv class=\"Description\"\u003e \u003cp\u003eOsteopontin\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003cdiv class=\"DefinitionListEntry\"\u003e \u003cdiv class=\"Term\"\u003eOSA\u003c/div\u003e \u003cdiv class=\"Description\"\u003e \u003cp\u003eOsteosarcoma\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003cdiv class=\"DefinitionListEntry\"\u003e \u003cdiv class=\"Term\"\u003ePDGF\u003c/div\u003e \u003cdiv class=\"Description\"\u003e \u003cp\u003ePlatelet-derived growth factor\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003cdiv class=\"DefinitionListEntry\"\u003e \u003cdiv class=\"Term\"\u003eSPP1\u003c/div\u003e \u003cdiv class=\"Description\"\u003e \u003cp\u003eSecreted Phosphoprotein 1\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003cdiv class=\"DefinitionListEntry\"\u003e \u003cdiv class=\"Term\"\u003eVEGF\u003c/div\u003e \u003cdiv class=\"Description\"\u003e \u003cp\u003eVascular endothelial growth factor\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003c/div\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eEthics approval and consent to participate\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis study followed the guidelines approved by the VTH -MSU Ethics Committee. The dog in this study was examined and treated with the owner\u0026apos;s written consent.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConsent for publication\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe owner signed an informed consent form including the authorization to publish the data arising from this study.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of data and materials\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eData is available from the authors upon reasonable request and with permission of NanoString.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors received no external financial support for the research, authorship, and/or publication of this article.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthors\u0026apos; contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;P.V.S. and V.Y.G conceived the original idea and design. K.K., J.D., and P.V.S. were involved in clinical data acquisitions and interpretations. K.K., A.I.E., Y.-T. Y, P.V.S, and V.Y.-G. were involved in the production, acquisition, and interpretation of transcriptomic data.\u003c/p\u003e\n\u003cp\u003eK.K., J. D, Y.-T. Y, A.I.E., V.Y.-G., and P.V.S. were involved in writing, reviewing, and editing the manuscript. P.V.S. and V.Y.-G. contributed to the final critical revision of the manuscript on their experience on the topic. All authors read and approved the final manuscript.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgements\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors thank the Genomics Core of the Research and Technology Support Facility of Michigan State University for the NanoString assay.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthors\u0026apos; information\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003csup\u003e1\u003c/sup\u003e Department of Small Animal Clinical Science, College of Veterinary Medicine, University of Florida, Gainesville, FL 32608\u003c/p\u003e\n\u003cp\u003e\u003csup\u003e2\u003c/sup\u003e Veterinary Emergency Group, Anaheim, CA 92808\u003c/p\u003e\n\u003cp\u003e\u003csup\u003e3\u003c/sup\u003e Department of Small Animal Clinical Sciences, College of Veterinary Medicine, Michigan State University, East Lansing, MI 48824, USA.\u003c/p\u003e\n\u003cp\u003e\u003csup\u003e4\u0026nbsp;\u003c/sup\u003eDepartment of Microbiology, Genetics and Immunology, College of Veterinary Medicine, Michigan State University, East Lansing, MI 48824\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eBrodey RS. Hypertrophic osteoarthropathy in the dog: a clinicopathologic survey of 60 cases. J Am Vet Med Assoc. 1971;159(10):1242-56.\u003c/li\u003e\n\u003cli\u003eWithers SS, Johnson EG, Culp WT, Rodriguez CO, Jr., Skorupski KA, Rebhun RB. Paraneoplastic hypertrophic osteopathy in 30 dogs. Vet Comp Oncol. 2015;13(3):157-65.\u003c/li\u003e\n\u003cli\u003ePacham\u0026eacute; AV, Beltr\u0026aacute;n MJ, Recchiuti NE, Portiansky EL, Massone AR, Gimeno EJ. Lesions of hypertrophic osteopathy in the forelimbs of a dog associated with pulmonary metastasis from a periosteal osteosarcoma. 2017.\u003c/li\u003e\n\u003cli\u003eRamoo S. Hypertrophic osteopathy associated with two pulmonary tumours and myocardial metastases in a dog: a case report. New Zealand Veterinary Journal. 2013;61(1):45-8.\u003c/li\u003e\n\u003cli\u003eLiptak JM, Monnet E, Dernell WS, Withrow SJ. Pulmonary metastatectomy in the management of four dogs with hypertrophic osteopathy. Vet Comp Oncol. 2004;2(1):1-12.\u003c/li\u003e\n\u003cli\u003eTurner H, S\u0026eacute;guin B, Worley DR, Ehrhart NP, Lafferty MH, Withrow SJ, et al. Prognosis for dogs with stage III osteosarcoma following treatment with amputation and chemotherapy with and without metastasectomy. J Am Vet Med Assoc. 2017;251(11):1293-305.\u003c/li\u003e\n\u003cli\u003eTreasure T. Hypertrophic pulmonary osteoarthropathy and the vagus nerve: an historical note. J R Soc Med. 2006;99(8):388-90.\u003c/li\u003e\n\u003cli\u003eNomori H, Kobayashi R, Kubo A, Morinaga S, Shintani Y, Sano T. Lung cancer containing growth hormone-releasing hormone associated with hypertrophic osteoarthropathy. Scandinavian journal of thoracic and cardiovascular surgery. 1994;28(3-4):149-52.\u003c/li\u003e\n\u003cli\u003eMito K, Maruyama R, Uenishi Y, Arita K, Kawano H, Kashima K, et al. Hypertrophic Pulmonary Osteoarthropathy Associated with Non-small Cell Lung Cancer Demonstrated Growth Hormone-Releasing Hormone by Immunohistochemical Analysis. Internal Medicine. 2001;40(6):532-5.\u003c/li\u003e\n\u003cli\u003eOoi A, Saad RA, Moorjani N, Amer KM. Effective symptomatic relief of hypertrophic pulmonary osteoarthropathy by video-assisted thoracic surgery truncal vagotomy. Ann Thorac Surg. 2007;83(2):684-5.\u003c/li\u003e\n\u003cli\u003eFenger JM, London CA, Kisseberth WC. Canine osteosarcoma: a naturally occurring disease to inform pediatric oncology. ILAR J. 2014;55(1):69-85.\u003c/li\u003e\n\u003cli\u003eCam M, Gardner HL, Roberts RD, Fenger JM, Guttridge DC, London CA, et al. \u0026Delta;Np63 mediates cellular survival and metastasis in canine osteosarcoma. Oncotarget. 2016;7(30):48533-46.\u003c/li\u003e\n\u003cli\u003eFossey SL, Liao AT, McCleese JK, Bear MD, Lin J, Li PK, et al. Characterization of STAT3 activation and expression in canine and human osteosarcoma. BMC Cancer. 2009;9:81.\u003c/li\u003e\n\u003cli\u003ePortela RF, Fadl-Alla BA, Pondenis HC, Byrum ML, Garrett LD, Wycislo KL, et al. Pro-tumorigenic effects of transforming growth factor beta 1 in canine osteosarcoma. J Vet Intern Med. 2014;28(3):894-904.\u003c/li\u003e\n\u003cli\u003eByrum ML, Pondenis HC, Fredrickson RL, Wycislo KL, Fan TM. Downregulation of CXCR4 Expression and Functionality After Zoledronate Exposure in Canine Osteosarcoma. J Vet Intern Med. 2016;30(4):1187-96.\u003c/li\u003e\n\u003cli\u003eRoy J, Wycislo KL, Pondenis H, Fan TM, Das A. Comparative proteomic investigation of metastatic and non-metastatic osteosarcoma cells of human and canine origin. PLoS One. 2017;12(9):e0183930.\u003c/li\u003e\n\u003cli\u003eYang Y-T, Engleberg AI, Yuzbasiyan-Gurkan V. Establishment and Characterization of Cell Lines from Canine Metastatic Osteosarcoma. Cells. 2024;13(1):25.\u003c/li\u003e\n\u003cli\u003eDanaher P, Warren S, Dennis L, D\u0026apos;Amico L, White A, Disis ML, et al. Gene expression markers of Tumor Infiltrating Leukocytes. J Immunother Cancer. 2017;5:18.\u003c/li\u003e\n\u003cli\u003eAlexa A, Rahnenfuhrer J. topGO: enrichment analysis for gene ontology. R package version. 2010;2(0):2010.\u003c/li\u003e\n\u003cli\u003eMitchell AL, Attwood TK, Babbitt PC, Blum M, Bork P, Bridge A, et al. InterPro in 2019: improving coverage, classification and access to protein sequence annotations. Nucleic Acids Res. 2019;47(D1):D351-d60.\u003c/li\u003e\n\u003cli\u003eGeer LY, Marchler-Bauer A, Geer RC, Han L, He J, He S, et al. The NCBI BioSystems database. Nucleic Acids Res. 2010;38(Database issue):D492-6.\u003c/li\u003e\n\u003cli\u003eSubramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545-50.\u003c/li\u003e\n\u003cli\u003eLiberzon A, Subramanian A, Pinchback R, Thorvaldsd\u0026oacute;ttir H, Tamayo P, Mesirov JP. Molecular signatures database (MSigDB) 3.0. Bioinformatics. 2011;27(12):1739-40.\u003c/li\u003e\n\u003cli\u003eFabregat A, Jupe S, Matthews L, Sidiropoulos K, Gillespie M, Garapati P, et al. The Reactome Pathway Knowledgebase. Nucleic Acids Res. 2018;46(D1):D649-d55.\u003c/li\u003e\n\u003cli\u003eSlenter DN, Kutmon M, Hanspers K, Riutta A, Windsor J, Nunes N, et al. WikiPathways: a multifaceted pathway database bridging metabolomics to other omics research. Nucleic Acids Res. 2018;46(D1):D661-d7.\u003c/li\u003e\n\u003cli\u003eMalkov VA, Serikawa KA, Balantac N, Watters J, Geiss G, Mashadi-Hossein A, et al. Multiplexed measurements of gene signatures in different analytes using the Nanostring nCounter Assay System. BMC Res Notes. 2009;2:80.\u003c/li\u003e\n\u003cli\u003eGeiss GK, Bumgarner RE, Birditt B, Dahl T, Dowidar N, Dunaway DL, et al. Direct multiplexed measurement of gene expression with color-coded probe pairs. Nat Biotechnol. 2008;26(3):317-25.\u003c/li\u003e\n\u003cli\u003eHan X, Wang W, He J, Jiang L, Li X. Osteopontin as a biomarker for osteosarcoma therapy and prognosis. Oncol Lett. 2019;17(3):2592-8.\u003c/li\u003e\n\u003cli\u003eLund SA, Giachelli CM, Scatena M. The role of osteopontin in inflammatory processes. J Cell Commun Signal. 2009;3(3-4):311-22.\u003c/li\u003e\n\u003cli\u003eWai PY, Kuo PC. The role of Osteopontin in tumor metastasis. J Surg Res. 2004;121(2):228-41.\u003c/li\u003e\n\u003cli\u003eZhao J, Dong L, Lu B, Wu G, Xu D, Chen J, et al. Down-regulation of osteopontin suppresses growth and metastasis of hepatocellular carcinoma via induction of apoptosis. Gastroenterology. 2008;135(3):956-68.\u003c/li\u003e\n\u003cli\u003eMori R, Shaw TJ, Martin P. Molecular mechanisms linking wound inflammation and fibrosis: knockdown of osteopontin leads to rapid repair and reduced scarring. J Exp Med. 2008;205(1):43-51.\u003c/li\u003e\n\u003cli\u003eAshkar S, Weber GF, Panoutsakopoulou V, Sanchirico ME, Jansson M, Zawaideh S, et al. Eta-1 (osteopontin): an early component of type-1 (cell-mediated) immunity. Science. 2000;287(5454):860-4.\u003c/li\u003e\n\u003cli\u003eGudowska-Sawczuk M, Mroczko B. What Is Currently Known about the Role of CXCL10 in SARS-CoV-2 Infection? Int J Mol Sci. 2022;23(7).\u003c/li\u003e\n\u003cli\u003eHuang H, Zhou W, Chen R, Xiang B, Zhou S, Lan L. CXCL10 is a Tumor Microenvironment and Immune Infiltration Related Prognostic Biomarker in Pancreatic Adenocarcinoma. Frontiers in Molecular Biosciences. 2021;8.\u003c/li\u003e\n\u003cli\u003eAriyarathna H, Thomson N, Aberdein D, Munday JS. Chemokine gene expression influences metastasis and survival time of female dogs with mammary carcinoma. Vet Immunol Immunopathol. 2020;227:110075.\u003c/li\u003e\n\u003cli\u003eRusso E, Santoni A, Bernardini G. Tumor inhibition or tumor promotion? The duplicity of CXCR3 in cancer. Journal of Leukocyte Biology. 2020;108(2):673-85.\u003c/li\u003e\n\u003cli\u003eBianchi ME, Mezzapelle R. The Chemokine Receptor CXCR4 in Cell Proliferation and Tissue Regeneration. Front Immunol. 2020;11:2109.\u003c/li\u003e\n\u003cli\u003eCao Y, Hunter ZR, Liu X, Xu L, Yang G, Chen J, et al. The WHIM-like CXCR4(S338X) somatic mutation activates AKT and ERK, and promotes resistance to ibrutinib and other agents used in the treatment of Waldenstrom\u0026apos;s Macroglobulinemia. Leukemia. 2015;29(1):169-76.\u003c/li\u003e\n\u003cli\u003eLaverdiere C, Hoang BH, Yang R, Sowers R, Qin J, Meyers PA, et al. Messenger RNA Expression Levels of CXCR4 Correlate with Metastatic Behavior and Outcome in Patients with Osteosarcoma. Clinical Cancer Research. 2005;11(7):2561-7.\u003c/li\u003e\n\u003cli\u003eWysolmerski JJ, Philbrick WM, Dunbar ME, Lanske B, Kronenberg H, Broadus AE. Rescue of the parathyroid hormone-related protein knockout mouse demonstrates that parathyroid hormone-related protein is essential for mammary gland development. Development. 1998;125(7):1285-94.\u003c/li\u003e\n\u003cli\u003ePhilbrick WM, Dreyer BE, Nakchbandi IA, Karaplis AC. Parathyroid hormone-related protein is required for tooth eruption. Proc Natl Acad Sci U S A. 1998;95(20):11846-51.\u003c/li\u003e\n\u003cli\u003eFoley J, Longely BJ, Wysolmerski JJ, Dreyer BE, Broadus AE, Philbrick WM. PTHrP regulates epidermal differentiation in adult mice. J Invest Dermatol. 1998;111(6):1122-8.\u003c/li\u003e\n\u003cli\u003eYang R, Hoang BH, Kubo T, Kawano H, Chou A, Sowers R, et al. Over-expression of parathyroid hormone Type 1 receptor confers an aggressive phenotype in osteosarcoma. International Journal of Cancer. 2007;121(5):943-54.\u003c/li\u003e\n\u003cli\u003eAl-Khan AA, Al Balushi NR, Richardson SJ, Danks JA. Roles of Parathyroid Hormone-Related Protein (PTHrP) and Its Receptor (PTHR1) in Normal and Tumor Tissues: Focus on Their Roles in Osteosarcoma. Front Vet Sci. 2021;8:637614.\u003c/li\u003e\n\u003cli\u003eHu C, Wen J, Gong L, Chen X, Wang J, Hu F, et al. Thrombospondin-1 promotes cell migration, invasion and lung metastasis of osteosarcoma through FAK dependent pathway. Oncotarget. 2017;8(44):75881-92.\u003c/li\u003e\n\u003cli\u003eLiu Y, Cai Y, Liu L, Wu Y, Xiong X. Crucial biological functions of CCL7 in cancer. PeerJ. 2018;6:e4928.\u003c/li\u003e\n\u003cli\u003eAhmed H, AlSadek DM. Galectin-3 as a Potential Target to Prevent Cancer Metastasis. Clin Med Insights Oncol. 2015;9:113-21.\u003c/li\u003e\n\u003cli\u003eMainetti LE, Zhe X, Diedrich J, Saliganan AD, Cho WJ, Cher ML, et al. Bone-induced c-kit expression in prostate cancer: a driver of intraosseous tumor growth. Int J Cancer. 2015;136(1):11-20.\u003c/li\u003e\n\u003cli\u003eScheu S, Ali S, Ruland C, Arolt V, Alferink J. The C-C Chemokines CCL17 and CCL22 and Their Receptor CCR4 in CNS Autoimmunity. Int J Mol Sci. 2017;18(11).\u003c/li\u003e\n\u003cli\u003eRauniyar K, Jha SK, Jeltsch M. Biology of Vascular Endothelial Growth Factor C in the Morphogenesis of Lymphatic Vessels. Front Bioeng Biotechnol. 2018;6:7.\u003c/li\u003e\n\u003cli\u003eLiu W, Hu H, Shao Z, Lv X, Zhang Z, Deng X, et al. Characterizing the tumor microenvironment at the single-cell level reveals a novel immune evasion mechanism in osteosarcoma. Bone Research. 2023;11(1):4.\u003c/li\u003e\n\u003cli\u003ePark HR, Min K, Kim HS, Jung WW, Park YK. Expression of vascular endothelial growth factor-C and its receptor in osteosarcomas. Pathol Res Pract. 2008;204(8):575-82.\u003c/li\u003e\n\u003cli\u003eYang W-H, Chang A-C, Wang S-W, Wang S-J, Chang Y-S, Chang T-M, et al. Leptin promotes VEGF-C production and induces lymphangiogenesis by suppressing miR-27b in human chondrosarcoma cells. Scientific Reports. 2016;6(1):28647.\u003c/li\u003e\n\u003cli\u003eWang Z-Q, Liang J, Schellander K, Wagner EF, Grigoriadis AE. c-fos-induced Osteosarcoma Formation in Transgenic Mice: Cooperativity with c-jun and the Role of Endogenous c-fos1. Cancer Research. 1995;55(24):6244-51.\u003c/li\u003e\n\u003cli\u003eSasaki K, Takano S, Tomizawa S, Miyahara Y, Furukawa K, Takayashiki T, et al. C4b-binding protein \u0026alpha;-chain enhances antitumor immunity by facilitating the accumulation of tumor-infiltrating lymphocytes in the tumor microenvironment in pancreatic cancer. Journal of Experimental \u0026amp; Clinical Cancer Research. 2021;40(1):212.\u003c/li\u003e\n\u003cli\u003eLiu K, He Q, Liao G, Han J. Identification of critical genes and gene interaction networks that mediate osteosarcoma metastasis to the lungs. Exp Ther Med. 2015;10(5):1796-806.\u003c/li\u003e\n\u003cli\u003eCriado Garc\u0026iacute;a O, S\u0026aacute;nchez-Corral P, Rodr\u0026iacute;guez de C\u0026oacute;rdoba S. Isoforms of human C4b-binding protein. II. Differential modulation of the C4BPA and C4BPB genes by acute phase cytokines. J Immunol. 1995;155(8):4037-43.\u003c/li\u003e\n\u003cli\u003eGarc\u0026iacute;a de Frutos P, Alim RI, H\u0026auml;rdig Y, Z\u0026ouml;ller B, Dahlb\u0026auml;ck B. Differential regulation of alpha and beta chains of C4b-binding protein during acute-phase response resulting in stable plasma levels of free anticoagulant protein S. Blood. 1994;84(3):815-22.\u003c/li\u003e\n\u003cli\u003eAgostinis C, Zorzet S, Balduit A, Zito G, Mangogna A, Macor P, et al. The Inflammatory Feed-Forward Loop Triggered by the Complement Component C3 as a Potential Target in Endometriosis. Frontiers in Immunology. 2021;12.\u003c/li\u003e\n\u003cli\u003ePisarcik D, Kousourou G, Vinod N. A Case of Paraneoplastic Arthralgias. American Journal of Medical Case Reports. 2021;9(8):397-401.\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"veterinary-oncology","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"","sideBox":"Learn more about [Veterinary Oncology](https://veterinaryoncology.biomedcentral.com/)","snPcode":"44356","submissionUrl":"https://submission.springernature.com/new-submission/44356/3","title":"Veterinary Oncology","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"BMC/SO AJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"Hypertrophic osteopathy, NanoString, osteosarcoma, paraneoplastic syndrome, gene expression profile","lastPublishedDoi":"10.21203/rs.3.rs-4438777/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4438777/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eHypertrophic osteopathy (HO) is a paraneoplastic syndrome, and the most notable cause in dogs is pulmonary metastatic osteosarcoma (OSA). Although many molecular factors in canine OSA have been shown in metastasis, little is known about the gene expression profile of HO secondary to metastatic OSA. Therefore, the purpose of this study was to compare the gene expression profiles between primary and metastatic OSA lesions from the same dog and to look for gene expression changes that can elucidate the molecular mechanism of metastases and HO.\u003c/p\u003e \u003cp\u003eTumoral samples were obtained from a 2-year-old, intact male, Labrador retriever. At the first visit, the patient presented with an appendicular OSA as the primary lesion. About 10 months later, the dog developed HO due to a single pulmonary metastasis. Using these primary and metastatic samples from the same dog, as well as normal canine osteoblasts, we investigated the gene expression profiling using the NanoString nCounter\u0026reg; Canine IO panel.\u003c/p\u003e \u003cp\u003eA total of 180 differentially expressed genes were identified between malignant OSA cells and non-malignant canine osteoblasts. Furthermore, 5 genes (CCL17, VEGFC, C3, C4BPA, and FOS) were differentially expressed in comparison between primary and metastatic OSA samples. CCL17 and VEGFC were upregulated in the primary lesion compared to the metastatic lesion, while C3, C4BPA, and FOS were downregulated in the primary lesion relative to the metastatic lesion. Given that the metastatic lesion was relevant to the development of HO, the different gene expression profiles may be relevant to understanding the pathophysiology of HO.\u003c/p\u003e","manuscriptTitle":"Gene expression profile comparison of primary and pulmonary metastatic lesions in a dog with appendicular osteosarcoma and hypertrophic osteopathy","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-06-04 21:41:21","doi":"10.21203/rs.3.rs-4438777/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2024-07-17T18:28:43+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2024-07-17T16:29:10+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"308962329533221702661683315283924171966","date":"2024-07-12T07:39:06+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"313815889560627408940200398941175159906","date":"2024-07-09T16:56:45+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2024-06-12T03:25:49+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"242681070665775870807205996483035531198","date":"2024-05-22T16:19:03+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"14286148632401249030577937448071116354","date":"2024-05-21T21:53:36+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2024-05-20T21:23:13+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2024-05-19T23:49:40+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2024-05-19T23:49:40+00:00","index":"","fulltext":""},{"type":"submitted","content":"Veterinary Oncology","date":"2024-05-17T22:43:16+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"veterinary-oncology","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"","sideBox":"Learn more about [Veterinary Oncology](https://veterinaryoncology.biomedcentral.com/)","snPcode":"44356","submissionUrl":"https://submission.springernature.com/new-submission/44356/3","title":"Veterinary Oncology","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"BMC/SO AJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"9bbe7d4f-68a6-4eee-b420-a42350a880af","owner":[],"postedDate":"June 4th, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"under-review","subjectAreas":[],"tags":[],"updatedAt":"2024-09-03T13:13:07+00:00","versionOfRecord":[],"versionCreatedAt":"2024-06-04 21:41:21","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-4438777","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4438777","identity":"rs-4438777","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}