- Research article
- Open Access
Chromatin run-on sequencing analysis finds that ECM remodeling plays an important role in canine hemangiosarcoma pathogenesis
BMC Veterinary Research volume 16, Article number: 206 (2020)
Canine visceral hemangiosarcoma (HSA) is a highly aggressive cancer of endothelial origin that closely resembles visceral angiosarcoma in humans, both clinically and histopathologically. Currently there is an unmet need for new diagnostics and therapies for both forms of this disease. The goal of this study was to utilize Chromatin run-on sequencing (ChRO-seq) and immunohistochemistry (IHC) to identify gene and protein expression signatures that may be important drivers of HSA progression.
ChRO-seq was performed on tissue isolated from 17 HSA samples and 4 normal splenic samples. Computational analysis was then used to identify differentially expressed genes and these factors were subjected to gene ontology analysis. ChRO-seq analysis revealed over a thousand differentially expressed genes in HSA tissue compared with normal splenic tissue (FDR < 0.005). Interestingly, the majority of genes overexpressed in HSA tumor tissue were associated with extracellular matrix (ECM) remodeling. This observation correlated well with our histological analysis, which found that HSA tumors contain a rich and complex collagen network. Additionally, we characterized the protein expression patterns of two highly overexpressed molecules identified in ChRO-seq analysis, podoplanin (PDPN) and laminin alpha 4 (LAMA4). We found that the expression of these two ECM-associated factors appeared to be largely limited to transformed endothelial cells within the HSA lesions.
Outcomes from this study suggest that ECM remodeling plays an important role in HSA progression. Additionally, our study identified two potential novel biomarkers of HSA, PDPN and LAMA4. Interestingly, given that function-blocking anti-PDPN antibodies have shown anti-tumor effects in mouse models of canine melanoma, our studies raise the possibility that these types of therapeutic strategies could potentially be developed for treating canine HSA.
Angiosarcomas (AS) are highly aggressive malignant tumors originating from endothelial cells. They account for approximately 2% of all soft tissue sarcomas in humans with the number of cases increasing significantly over the past 30 years [1,2,3,4,5]. For patients presenting with non-metastatic AS, the reported five-year survival rate is ~ 35%. These patients also have a 75% chance of local recurrence within 24 months and a 50% likelihood of metastases developing despite local treatment [4, 6]. When metastases are already present at the initial presentation, the five-year survival rate is poor, with a median survival time of just 3 months . Although published research on this rare tumor is increasing, we still know very little about the pathogenesis of this disease.
Canine hemangiosarcoma (HSA) is histopathologically similar to angiosarcoma, with both forms of this disease following a similar clinical course . Canine HSA is most commonly observed in the spleen, but also occurs in other organs such as the heart, liver, and dermis with the later form often being associated with ultraviolet radiation-associated oncogenesis. Prognosis for the visceral forms of HSA is poor, with most dogs dying from this disease within months of their diagnosis [8, 9]. For dogs diagnosed with splenic HSA, surgical removal of the spleen can increase the life expectancy up to 6 months, and when combined with chemotherapy, can prolong their life for up to 12 months [8, 9].
While all breeds are susceptible to hemangiosarcoma, German Shepherds, Golden Retrievers, and Labrador Retrievers, are particularly predisposed to developing this disease , with the estimated lifetime risk of Golden Retrievers developing HSA being 20%. This breed predisposition suggests a genetic component for HSA . Genome-wide association studies have been carried out to identify risk loci [12, 13]. SNP array analysis of genomic DNA from Golden Retrievers with HSA identified a risk locus on chromosome 5 that was shared by ~ 20% of the cases . In another study, microarray analysis of genomic DNA from dogs with HSA revealed gained copy number aberrations on several genes such as PDGFRA, KDR and VEGFA . A recent whole exome sequencing study of dogs with HSA across a range of breeds revealed somatic mutations in tumor suppressor genes, including TP53, and two genes (PIK3CA and PIK3R1) in the Phosphoinositide 3-kinase (PI3K) pathway . Comparison of somatic copy number aberration profiles in human angiosarcoma and canine hemangiosarcoma identified recurrent copy number gains in KDR (31% in human, 22% in canine) and KIT (17% in both) . These genome-wide studies reveal that, while specific genetic aberrations are associated with HSA in some populations, these alterations are not sufficient to explain the majority of HSA cases. These observations support the hypothesis that HSA pathogenesis is heterogeneous in nature.
Transcriptome analysis has been previously performed on cell lines and tumor tissues to identify molecular features that define canine HSA. Gene expression profiling of HSA cell lines and non-malignant proliferating endothelial cells revealed that HSA cell lines appeared to overexpress genes associated with inflammation, angiogenesis, adhesion, invasion, metabolism, cell cycle progression, and patterning . Microarray and RNA-seq analysis of visceral HSA tumors identified three distinct molecular subtypes; angiogenesis, inflammation and lipogenesis. These molecular subtypes did not appear to be associated with a specific breed or tumor morphologic subtype .
A variety of analytic tools exist for assaying the transcriptome and molecular alterations in tissue. One of these tools, chromatin run-on and sequencing (ChRO-seq), uses RNA polymerase activity to measure transcription and, as such, provides for a highly-sensitive, base-pair level readout of gene expression [18, 19]. Given the lack of clarity regarding the molecular underpinnings of HSA, the major goal of this study was to utilize ChRO-seq to document changes in gene expression between normal splenic tissue and HSA tumor tissue. Our results show that the majority of genes that are upregulated in HSA appear to be associated with extracellular matrix (ECM) remodeling. Additionally, we further characterize two ECM-associated molecules that were highly overexpressed in HSA tumor tissue, podoplanin (PDPN) and laminin alpha 4 (LAMA4). We show by immunohistochemistry (IHC) that the expression of these two cancer-associated factors appears to be primarily limited to HSA tumor cells.
ChRO-Seq analysis of transcription in HSA and normal splenic tissue
For our genome-wide study, we performed ChRO-seq analysis on tumor tissue from dogs histopathologically diagnosed with HSA (20) and on normal splenic tissue (4) (Table 1). We quantified the similarity of transcription between these tissue samples using Pol II abundance in annotated gene bodies. Three HSA samples (B297, B675, B788) were excluded due to a low number of mappable reads (< 2 million reads). Sequencing data from the remaining samples was then analyzed to create a Spearman’s rank correlation matrix (Fig. 1a). Seven thousand seven genes (> 20 rpkm) were used to calculate the correlation coefficients by GENE-E. Four normal tissues and 3 HSA tissues (B307, B829 and C349) were found to form one cluster, with 2 HSA samples (C001 and C034) also being similar to the Normal/HSA cluster. DESeq2 analysis was then performed on the samples to identify differentially expressed genes. When compared to normal splenic tissue, a total of 906 genes were upregulated in the HSA samples (FDR < 0.005) while 358 genes were downregulated (FDR < 0.005). An MA plot showing the differentially expressed gene set is shown in Fig. 1b. In this plot, the X axis represents the average expression over all samples and the Y axis represents the log2 fold change between the tumor samples and normal spleen. Highly upregulated extracellular matrix (ECM)-associated genes and downregulated heme synthesis-associated genes are labelled with blue dots (Fig. 1b). Differentially expressed genes (FDR < 0.001, 659 genes) were then assessed to evaluate expression levels in individual samples (Fig. 1c). Two hundred and ninety-one genes (rpkm > 20) are shown in the resulting heatmap which also contains an associated hierarchical clustering dendrogram (one minus spearman rank correlation with complete linkage). Results show that the expression profiles of HSA samples appear to be distinct from normal splenic tissue.
Extracellular matrix remodeling appears to be a major feature of HSA
When evaluating the top 50 most highly upregulated genes in the HSA group, we observed that many of these molecules are associated with the ECM, examples of which include podoplanin (PDPN), laminin alpha 4 (LAMA4), ADAM12, and fibronectin 1 (Fig. 1b). To test the prevalence of ECM-associated molecules in our dataset, we carried out gene ontology (GO) analysis using PANTHER (protein analysis through evolutionary relationships) (Fig. 2). We utilized a total of 369 of the most highly upregulated genes (FDR < 0.0005) as an input list, with 282 of these genes then being identified as uniquely mapped genes in a Canis lupus familiaris reference database. A total of 20 annotation categories were detected (FDR < 0.05) from the reactome pathway database, with the top 10 annotation categories from this dataset being found to be related to ECM regulation. A more detailed list of the ECM-associated factors that were upregulated in HSA tumor tissue is shown in Table 2, along with their fold change (log2) and FDR. A full list of Reactome pathways terms and the list of genes is shown in the supplemental table.
Interestingly, when we queried genes in the downregulated dataset (124 genes with an FDR cut-off of 0.0005), we did not observe statistically significant functional associations using GO analysis. Therefore, we adjusted the FDR cut-off to 0.005 and this adjustment identified 358 genes as being downregulated in the HSA samples with 301 of these genes being uniquely mapped in the reference database. As shown in Fig. 2, five annotation categories were detected in this dataset, with the most enriched reactome pathway term being heme biosynthesis. Examples of specific genes that are associated with heme biosynthesis include, UROS, FECH, UROD, JCHAIN, and ALAD (Fig. 1b). This result suggests that the expression of genes involved in normal splenic function is suppressed in HSA tumor tissue.
Characterization of ECM-associated proteins in HSA tumor tissue
To more directly test the hypothesis that ECM-associated factors are overexpressed in HSA tumor tissue, we first stained paraffin sections of tumor and normal tissue with Masson’s trichrome to label collagen fibers (Fig. 3). Staining of the more vascular regions of the HSA tumor tissue show that collagen fibrils are abundant throughout these regions. This observation fits well with our finding that multiple collagen isoforms were highly overrepresented in our HSA-overexpression ChRO-Seq dataset. Interestingly, in these vascular regions, the malignant tumor cells appear to surround the collagen bundles, forming what resembles inverted blood vessels. To further support these findings, images and quantitative data from six additional HSA and two normal splenic samples are shown in supplemental figure, S1.
We next further characterized two ECM-associated factors that were highly overexpressed in HSA tumor tissue. For this analysis we focused on PDPN (FDR 6.43E-05, fold change 3.62) and LAMA4 (FDR 1.58E-09, fold change 3.77) because these molecules have been previously associated with cancer progression in other tumor types [20,21,22,23,24,25]. For further characterization, we first investigated PDPN and LAMA4 expression levels in our ChRO-seq dataset. A genome browser view of transcription levels from one normal and three HSA samples is shown in Fig. 4a. Expression levels from all ChRO-seq samples are shown in supplemental figure S2. Results show that, at the PDPN locus, two of the HSA cases (C085 and C356) showed increased RNA polymerase activity compared to normal tissue, while expression levels of the third HSA sample (C340) was comparable with normal tissue. LAMA4 was transcribed from the minus strand and all HSA samples showed higher transcriptional activity than normal tissue. We next investigated the expression of mature PDPN and LAMA4 transcripts in tumor tissue and normal spleen using RT-PCR. Results show that, while PDPN or LAMA4 transcripts were either absent or weakly expressed in normal spleen, LAMA4 was significantly expressed in all three tested HSA tumor samples while PDPN was robustly expressed in one (C356) of the three HSA samples (Fig. 4b and supplemental Figure S7).
Next, we performed IHC analysis of PDPN and LAMA4 on HSA tumor sections to confirm that these molecules are expressed in tumor tissue at the protein level and to investigate their localization. HSA samples B297, B783 and B848 are shown in Figs. 5 and 6 with images and quantitative data from additional samples being shown in supplemental Figure S3 and S4. IHC analysis found that anti-PDPN staining in normal splenic tissue was low to absent. In HSA tumor tissue, however, anti-PDPN staining was robust in the cytoplasm of tumor cells encircling the collagen fibrils (Fig. 5c and d, arrow). In addition to tumor cells, anti-PDPN staining was also evident in the cytoplasm of endothelial cells in tumor-adjacent blood vessels (Fig. 5b, arrowhead). The localization pattern of anti-LAMA4 staining was similar to that of PDPN, with little staining being seen in normal splenic tissue and strong anti-LAMA4 staining being observed in the cytoplasm of HSA tumor cells encircling collagen fibrils. Again, as with PDPN, anti-LAMA4 staining was also observed in the cytoplasm in the endothelial cells from tumor-adjacent blood vessels (Fig. 6c and d).
Serial sections of additional HSA samples were then stained with Masson’s trichrome and immunohistochemistry was performed on these additional tissues using anti-PDPN and LAMA4 antibodies (Supplemental figure S5 and S6). Sections from one sample (B176, Figure S5) showed solid/cavernous features with blood vessels (Figure S5 a). Masson’s trichrome staining highlighted the presence of collagen fibrils filling spaces between regions of neoplastic cells (Figure S5 a and b). Anti-PDPN staining appeared to be primarily cytosolic in these neoplastic cells and in endothelial cells lining adjacent blood vessels (Figure S5 c). Anti-LAMA4 staining was strong in extracellular regions of the tumor and in endothelial cells lining blood vessels (Figure S5 d). Sections from sample B554 primarily showed capillary/cavernous features (Figure S6). Masson’s trichrome staining of this sample revealed a core collagen region surrounded by malignant endothelial cells with spindle shaped nuclei (Figure S6 a and b). Anti-PDPN and LAMA4 staining appeared to be primarily cytosolic in tumor cells from this sample.
In this study, we first utilized ChRO-seq to identify RNA polymerase activity in hemangiosarcoma tumor tissue and normal splenic tissue. Analysis of the correlation matrix for our ChRO-seq dataset finds that 14 of the HSA tumor samples clustered together, while 3 of the HSA samples appeared to be more similar to normal splenic tissue (Fig. 1). While speculative, it is possible that these 3 samples represent a different subtype of HSA and, as such, displayed a gene expression profile that is more similar to normal tissue. Alternatively, it is also possible that these 3 presumptive HSA tissues were actually tumor-adjacent normal tissue or a mix of normal and tumor tissue. While grading of canine hemangiosarcoma is not often utilized due to the aggressive nature of the neoplasm, a grading scheme does exist and differences in gene expression between the samples may correlate with progression and eventual outcome .
Previous HSA transcriptomic profiling studies have implicated a number of signaling pathways in HSA pathogenesis. Tamburini et al., for example, found that, when compared to cell lines derived from splenic hematomas, HSA cell lines exhibited several different distinct gene expression profiles, including signatures for angiogenesis, inflammation, adhesion, invasion metabolism, cell cycle and signaling. In another study, Gorden et al. performed microarray and RNA-seq analyses on visceral HSA tumor tissues (12 spleen, 7 heart, 4 liver, 1 lung) and identified three distinct tumor subtypes that were associated with either angiogenesis, inflammation, and adipogenesis. Our finding that ECM remodeling appears to be a major gene expression signature in HSA differs from these previous studies. One possible reason for the difference in outcomes between our study and the Tamburini study is that our report compared gene expression patterns in HSA splenic tissue with normal splenic tissue while the Tamburini study was comparing the gene expression profiles of HSA cell lines with splenic hematoma cell lines . Thus, the differences may be due to the fact that the expression profile of normal splenic samples differs significantly from hematoma samples. Additionally, these differences may also be because the expression profile of cell lines may differ significantly from primary tissue samples due to decreased cellular complexity and prolonged cell culture. One potential reason why outcomes from our study differed from the Gorden study is that we evaluated differences in gene expression patterns between HSA and normal splenic tissue while the other study analyzed gene expression patterns within HSA tissues. Lastly, another reason why outcomes from our study may have differed from both of these previous studies is that ChRO-seq analysis detects nascently synthesized RNA as opposed to mature transcripts, whose levels can be affected by a variety of factors including transcript stability.
Regarding the specific types of ECM factors that were upregulated in our dataset (Table 2), we identified genes encoding fibrous components of the ECM, including 12 collagen genes, 3 laminin genes, and fibronectin. Several ECM binding proteins were also identified in this dataset including, Lumican and Biglycan, which are small leucine-rich proteoglycans that regulate collagen fibril and matrix assembly [27,28,29,30]. Integrin alpha 2 is directly associated with fibril-forming collagens (1,2,3,5,6,14, 18)  while DDR2 functions as a cell surface receptor for fibrillar collagen and regulates cell differentiation by remodeling ECM [32, 33]. Additionally, integrin alpha 5 binds directly to fibronectin while integrin alpha 6 binds to laminin [34, 35]. Several ECM-related enzymes were also found in this dataset including P4HA2 and PLOD1 which catalyze collagen biosynthesis [36, 37] along with ADAMTS14, ADAMTS4, ADAM12 and TLL1 which process procollagen to collagen by cleaving N-propeptide and C- propeptide [38, 39]. Lastly, two molecules directly involved in ECM turnover, tissue inhibitor of metalloproteinase 1 (TIMP1) and matrix metalloproteinase 14 (MMP14), were also found in this dataset [40,41,42].
The extent to which ECM genes were upregulated in HSA tumors in our study is highlighted by our GO analysis which found that 9 of the top 10 biological process categories were ECM-related (Fig. 2). Interestingly, 5 of these 9 categories related to collagen function. We more directly tested for the abundance of collagen fibers in the tumor samples using Masson’s trichrome stain and found extensive collagen deposition throughout the tumor tissue (Fig. 3). Collagen deposition was observed both in the more solid areas of the tumor and, in particular, in the tumor regions filled with malformed vascular channels. In these vascular regions, neoplastic endothelial cells are often found to encircle the collagen bundles forming “gumball”-shaped structures that look like inverted blood vessels. While stromal fibroblasts are presumably primarily responsible for the synthesis of these collagen fibers, it is also possible that the tumor cells may also be partly responsible for synthesis of these collagen fibers.
In addition to collagen, we also found that PDPN and LAMA4 were highly overexpressed in HSA tumor tissue when compared to normal tissue. These molecules were of particular interest to us given their close association with cancer progression in other types of cancer. Podoplanin is a mucin-type protein that contains an extracellular region, transmembrane domain, and intracellular tail. It is widely known as a marker for lymphatic endothelial cells and also to play a critical role in heart and lung development and in development of the lymphatic endothelial system [43,44,45]. PDPN appears to play several roles in cancer progression. A number of studies have shown that PDPN expression in cancer cells promotes tumor cell proliferation and invasion [20,21,22]. In addition to cell-intrinsic effects, PDPN is also believed to promote tumor metastasis by interacting with its receptor, CLEC2, on the platelets. This interaction then promotes the coating of tumor cells by platelets, thereby protecting tumor cells from the immune system . In addition to its role in human cancers, PDPN is also overexpressed in canine squamous cell carcinomas and melanomas  and PDPN mAbs were recently found to have potent anti-tumor activity in mouse xenograft models of canine melanoma . Interestingly, overexpression of PDPN in mice leads to disseminated intravascular coagulation , a condition that is strikingly similar to that seen in many dogs with HSA . Our ChRO-seq data demonstrated that PDPN expression varies significantly between HSA samples (Figs. 1 and 4). IHC analysis found that PDPN was robustly expressed in transformed endothelial cells in certain HSA tumor samples. However, similar to our ChRO-seq finding, we did find that PDPN appears to only be expressed in a subset of HSA tumors. In future studies, we plan to test whether PDPN expression in HSA tumors correlates with disease severity.
Laminins form large heterotrimeric αβγ protein complexes and are a prominent component of basement membranes. LAMA4 is distinct from other laminin isoforms in that it lacks a polymerization domain with the loss of this domain potentially facilitating tumor cell migration . Interestingly, LAMA4 was recently described as an “oncolaminin” due to its strong association with cancer cell migration and tumor progression in a range of cancers . These links to cancer include recent studies which found that LAMA4 and MCAM (melanoma cell adhesion molecule) are highly enriched in tumor blood vessels in renal cell carcinoma and colorectal carcinoma. Additionally another study found that expression of these molecules is enhanced in locally invasive and metastatic clear cell renal cell carcinoma . Further, antisense oligonucleotides against laminin-8 (LAMA4 and LAMB1) were found to block the invasion of glioma cells and neovascularization in vitro . Taken together, these published studies indicate that LAMA4 plays a key role for vascular development, tumor progression and metastasis. In our study, deseq2 analysis found that LAMA4 is highly transcribed in HSA tumor tissue. This observation was supported by PCR analysis of mRNA isolated from HSA tumor tissue (Fig. 4 and supplemental Figure S7). Additionally, our IHC analysis found that LAMA4 protein expression appears to be primarily limited to malignant endothelial cells (Fig. 5). Given the previously defined roles for LAMA4 in cancer progression and tumor metastasis, our findings raise the possibility that LAMA4 is an important mediator of canine HSA pathobiology.
Outcomes from our studies found that the majority of upregulated genes in HSA tumor tissue appear to be associated with ECM remodeling. This finding was supported by IHC analysis which found a robust collagen network throughout HSA tumor tissue. Additionally, we further characterized two ECM-associated factors, PDPN and LAMA4, and found that their expression was largely limited to tumor cells in HSA tissue. Both of these molecules have been previously identified as potential biomarkers in other types of cancer. Given the recent previous preclinical mouse studies showing that anti-canine PDPN antibodies can block the growth of melanomas, our findings raise the possibility that similar types of therapies may have utility for treating canine HSA.
Canine tissue samples were obtained from the Canine Comparative Oncology & Genomics Consortium and Cornell University Hospital for Animals. Pathology was independently confirmed, and patient demographics are described in the supplemental material (Table 1). ChRO-seq libraries were made from splenic hemangiosarcoma (n = 20) and normal splenic tissues from HSA dogs (n = 4). Normal splenic samples from three beagles were either snap-frozen or fixed in PFA and Paraffin-embedded (FFPE) for sectioning. Additional FFPE sections of HSA cases were obtained from Cornell University Hospital for Animals.
Chromatin- run-on and sequencing
Chromatin was extracted from each tissue sample and chromatin run-on was performed as described previously [18, 19, 52]. ChRO-seq library preparations were executed according to the Illumina protocol and were sequenced using Illumina NextSeq500 sequencing. Raw sequence FASTQ files were assessed by FastQC for quality control . Single-read sequencing data were preprocessed and mapped to the canine genome assembly using the proseq 2.0 pipeline (https://github.com/Danko-Lab/proseq2.0) . Briefly, Single-read sequencing data were preprocessed to remove the adapter sequences and trimmed based on base quality, and deduplicating the reads based on unique molecular identifiers in RNA adapters. Sequencing reads were mapped into a canine genome assembly (CanFam 3.1) using the Burrows-Wheeler Aligner software package . The number of mappable reads are listed in Table 1. Aligned BAM files were converted into bigWig format, which were used to create the matrix of read counts mapping to each annotated gene (CanFam3 ensGene). Seven thousand seven genes (> 20 rpkm) were used to calculate the Spearman’s rank correlation coefficients to create the matrix and dendrogram by GENE-E (https://software.broadinstitute.org/GENE-E/). The list of differentially expressed genes and an MA-plot were produced by running DESeq2 (false discovery rate < 0.01) . Upregulated genes (369 genes, FDR < 0.0005) and down regulated genes (358 genes, FDR < 0.005) in HSA were subjected to PANTHER (Protein Analysis Through Evolutionary Relationships) overrepresentation gene ontogeny analysis (GO) on Reactome pathway dataset (version 65, released 2019-12-22) . Differentially expressed genes (FDR < 0.001, 659 genes) were assessed to show the expression level in individual samples. Heatmaps and dendrograms were created by GENE-E.
Immunohistochemistry and histology staining
Prior to all staining protocols, paraffin section slides were rehydrated in xylene followed by sequential washes in 100, 90, 80, and 70% EtOH. H&E and Masson’s trichrome staining were performed on 10 HSA samples and 4 normal spleen sections. Gill’s hematoxylin and Eosin Y solutions were used for H&E. Masson’s trichrome staining was performed following the manufacture’s protocol (Masson’s Trichrome Stain Kit, Polysciences, Inc. #25088–100). Immunohistochemistry (IHC) experiments using anti-canine PDPN (clone PMab-38 , now commercially available at FUJIFILM Wako Pure Chemical Corporation, Osaka, Japan; 017–27,091, 1:200) and anti-LAMA4 (Sigma-Aldrich Corp. St. Louis, MO, HPA015693, Human Atlas antibodies for IHC, Rabbit polyclonal, epitope homology 85%, 1:300), were performed on 8 HSA samples and 2 normal spleens following a standard protocol . Briefly, slides were boiled for 40 min in 0.01 M sodium citrate to retrieve antigens. Slides were cooled and washed in running tap water for 20 min, then incubated for 15 min in 3% hydrogen peroxide to quench endogenous peroxidases. Slides were blocked in 2.5% normal goat serum for 30 min at room temperature. Primary antibodies were diluted in 1% BSA/TBST, and incubated overnight at 4C. Slides were washed three times with TBS-T, then incubated with a secondary antibody for 1 h at room temperature (Goat anti-mouse-HRP, Vector laboratories, Inc., Burlingame, CA; #MP7452 or Goat anti-rabbit-HRP, Vector laboratories, Inc., #MP7451). ImmPACT NovaRED peroxidase substrate (Vector laboratories, Inc., SK-4805) with hematoxylin counterstain was used for signal detection. Slides were dehydrated through 80, 95 and 100% EtOH to Xylene and coverslipped. Blue staining and NovaRED signals were quantified using image J (see supplemental method and figures S1-S7).
Snap frozen tissues (3 normal spleens and 3 HSA tissues, see Table 1 for demographics) were pulverized in a tissue crusher and total RNA was extracted with Trizol reagent. 1 μg of RNA was reverse-transcribed to cDNA using High-capacity cDNA reverse transcription kit (Applied Biosystems) following the manufacture’s protocol (25 °C for 10 min, 37 °C for 120 min, 85 °C for 5 min). All cDNA reactions were diluted 5-fold with water prior to use in PCR. 1 μl of diluted cDNA was used for PCR using Go-Taq master mix (Promega Corporation) with canine specific primers (0.4 μM each, sequences are shown in supplemental material), products were amplified using the following parameters: 95 °C for 5 min, 30 cycles of 3 steps (95 °C for 30 s, 58 °C for 30 s, and 72 °C for 20 s) and 72 °C for 5 min. The amplified products were analyzed using agarose gel electrophoresis (1.2%).
Availability of data and materials
The raw files (FASTQ) and processed files (BigWig) analysed during the current study are available in the Gene Expression Omnibus GSE150705 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE150705).
Chromatin run-on sequencing
Laminin alpha 4
Protein analysis through evolutionary relationships
Fixed in PFA and Paraffin-embedded
Reverse transcription – polymerase chain reaction
Shustef E, Kazlouskaya V, Prieto VG, Ivan D, Aung PP. Cutaneous angiosarcoma: a current update. J Clin Pathol. 2017;70:917–25. https://doi.org/10.1136/jclinpath-2017-204601.
Young RJ, Brown NJ, Reed MW, Hughes D, Woll PJ. Angiosarcoma. Lancet Oncol. 2010;11:983–91. https://doi.org/10.1016/S1470-2045(10)70023-1.
Ye J, Li X-F, Wang Y-D, Yuan Y. Long-term survival of a patient with scalp angiosarcoma and multiple metastases treated using combination therapy: a case report. Oncol Lett. 2015;9:1725–8. https://doi.org/10.3892/ol.2015.2919.
Buehler D, Rice SR, Moody JS, Rush P, Hafez G-R, Attia S, et al. Angiosarcoma outcomes and prognostic factors: a 25-year single institution experience. Am J Clin Oncol. 2014;37:473–9. https://doi.org/10.1097/COC.0b013e31827e4e7b.
Mery CM, George S, Bertagnolli MM, Raut CP. Secondary sarcomas after radiotherapy for breast cancer. Cancer. 2009;115:4055–63. https://doi.org/10.1002/cncr.24462.
Mark RJ, Poen JC, Tran LM, Fu YS, Juillard GF. Angiosarcoma: a report of 67 patients and a review of the literature. Cancer. 1996;77:2400–6. https://doi.org/10.1002/(SICI)1097-0142(19960601)77:11<2400::AID-CNCR32>3.0.CO;2-Z.
Vail DM, MacEwen EG. Spontaneously occurring tumors of companion animals as models for human cancer. Cancer Investig. 2000;18:781–92 http://www.ncbi.nlm.nih.gov/pubmed/11107448. Accessed 17 Apr 2019.
Wendelburg KM, Price LL, Burgess KE, Lyons JA, Lew FH, Berg J. Survival time of dogs with splenic hemangiosarcoma treated by splenectomy with or without adjuvant chemotherapy: 208 cases (2001–2012). J Am Vet Med Assoc. 2015;247:393–403. https://doi.org/10.2460/javma.247.4.393.
Sorenmo KU, Jeglum KA, Helfand SC. Chemotherapy of canine hemangiosarcoma with doxorubicin and cyclophosphamide. J Vet Intern Med. 1993;7:370–6. https://doi.org/10.1111/j.1939-1676.1993.tb01033.x.
Cole PA. Association of canine splenic hemangiosarcomas and hematomas with nodular lymphoid hyperplasia or siderotic nodules. J Vet Diagnostic Investig. 2012;24:759–62. https://doi.org/10.1177/1040638712447580.
Kent MS, Burton JH, Dank G, Bannasch DL, Rebhun RB. Association of cancer-related mortality, age and gonadectomy in golden retriever dogs at a veterinary academic center (1989-2016). PLoS One. 2018;13:e0192578. https://doi.org/10.1371/journal.pone.0192578.
Tonomura N, Elvers I, Thomas R, Megquier K, Turner-Maier J, Howald C, et al. Genome-wide association study identifies shared risk loci common to two malignancies in golden retrievers. PLoS Genet. 2015;11:1–24.
Thomas R, Borst L, Rotroff D, Motsinger-Reif A, Lindblad-Toh K, Modiano JF, et al. Genomic profiling reveals extensive heterogeneity in somatic DNA copy number aberrations of canine hemangiosarcoma. Chromosom Res. 2014;22:305–19. https://doi.org/10.1007/s10577-014-9406-z.
Wang G, Wu M, Maloneyhuss MA, Wojcik J, Durham AC, Mason NJ, et al. Actionable mutations in canine hemangiosarcoma. PLoS One. 2017;12:e0188667. https://doi.org/10.1371/journal.pone.0188667.
Megquier K, Turner-Maier J, Swofford R, Kim J-H, Sarver AL, Wang C, et al. Comparative genomics reveals shared mutational landscape in canine hemangiosarcoma and human angiosarcoma. Mol Cancer Res. 2019. https://doi.org/10.1158/1541-7786.MCR-19-0221.
Tamburini BA, Phang TL, Fosmire SP, Scott MC, Trapp SC, Duckett MM, et al. Gene expression profiling identifies inflammation and angiogenesis as distinguishing features of canine hemangiosarcoma. BMC Cancer. 2010;10:619. https://doi.org/10.1186/1471-2407-10-619.
Gorden BH, Kim J-H, Sarver AL, Frantz AM, Breen M, Lindblad-Toh K, et al. Identification of three molecular and functional subtypes in canine hemangiosarcoma through gene expression profiling and progenitor cell characterization. Am J Pathol. 2014;184:985–95. https://doi.org/10.1016/j.ajpath.2013.12.025.
Mahat DB, Kwak H, Booth GT, Jonkers IH, Danko CG, Patel RK, et al. Base-pair-resolution genome-wide mapping of active RNA polymerases using precision nuclear run-on (PRO-seq). Nat Protoc. 2016;11:1455.
Chu T, Rice EJ, Booth GT, Salamanca HH, Wang Z, Core LJ, et al. Chromatin run-on and sequencing maps the transcriptional regulatory landscape of glioblastoma multiforme. Nat Genet. 2018;50:1553–64. https://doi.org/10.1038/s41588-018-0244-3.
Neri S, Ishii G, Hashimoto H, Kuwata T, Nagai K, Date H, et al. Podoplanin-expressing cancer-associated fibroblasts lead and enhance the local invasion of cancer cells in lung adenocarcinoma. Int J Cancer. 2015;137:784–96. https://doi.org/10.1002/ijc.29464.
Astarita JL, Acton SE, Turley SJ. Podoplanin: emerging functions in development, the immune system, and cancer. Front Immunol. 2012;3:283. https://doi.org/10.3389/fimmu.2012.00283.
Miyata K, Takemoto A, Okumura S, Nishio M, Fujita N. Podoplanin enhances lung cancer cell growth in vivo by inducing platelet aggregation. Sci Rep. 2017;7:4059. https://doi.org/10.1038/s41598-017-04324-1.
Pozzi A, Yurchenco PD, Iozzo RV. The nature and biology of basement membranes. Matrix Biol. 2017;57–58:1–11. https://doi.org/10.1016/j.matbio.2016.12.009.
Ishikawa T, Wondimu Z, Oikawa Y, Gentilcore G, Kiessling R, Egyhazi Brage S, et al. Laminins 411 and 421 differentially promote tumor cell migration via α6β1 integrin and MCAM (CD146). Matrix Biol. 2014;38:69–83. https://doi.org/10.1016/j.matbio.2014.06.002.
Wwragg J, Finnity JP, Anderson JA, Ferguson HJM, Porfiri E, Bhatt RI, et al. MCAM and LAMA4 are highly enriched in tumor blood vessels of renal cell carcinoma and predict patient outcome. Cancer Res. 2016;76:2314.
Moore AS, Rassnick KM, Frimberger AE. Evaluation of clinical and histologic factors associated with survival time in dogs with stage II splenic hemangiosarcoma treated by splenectomy and adjuvant chemotherapy: 30 cases (2011–2014). J Am Vet Med Assoc. 2017;251:559–65. https://doi.org/10.2460/javma.251.5.559.
Zhang G, Chen S, Goldoni S, Calder BW, Simpson HC, Owens RT, et al. Genetic evidence for the coordinated regulation of collagen fibrillogenesis in the cornea by decorin and biglycan. J Biol Chem. 2009;284:8888–97. https://doi.org/10.1074/jbc.M806590200.
Wiberg C, Heinegård D, Wenglén C, Timpl R, Mörgelin M. Biglycan organizes collagen VI into hexagonal-like networks resembling tissue structures. J Biol Chem. 2002;277:49120–6. https://doi.org/10.1074/jbc.M206891200.
Chakravarti S, Magnuson T, Lass JH, Jepsen KJ, LaMantia C, Carroll H. Lumican regulates collagen fibril assembly: skin fragility and corneal opacity in the absence of lumican. J Cell Biol. 1998;141:1277–86 http://www.ncbi.nlm.nih.gov/pubmed/9606218. Accessed 29 Apr 2019.
Appunni S, Anand V, Khandelwal M, Gupta N, Rubens M, Sharma A. Small leucine rich proteoglycans (decorin, biglycan and lumican) in cancer. Clin Chim Acta. 2019;491:1–7. https://doi.org/10.1016/j.cca.2019.01.003.
Heino J. Cellular signaling by collagen-binding Integrins. Dordrecht: Springer; 2014. p. 143–55. https://doi.org/10.1007/978-94-017-9153-3_10.
Wall SJ, Werner E, Werb Z, DeClerck YA. Discoidin domain receptor 2 mediates tumor cell cycle arrest induced by fibrillar collagen. J Biol Chem. 2005;280:40187–94. https://doi.org/10.1074/jbc.M508226200.
Majkowska I, Shitomi Y, Ito N, Gray NS, Itoh Y. Discoidin domain receptor 2 mediates collagen-induced activation of membrane-type 1 matrix metalloproteinase in human fibroblasts. J Biol Chem. 2017;292:6633–43. https://doi.org/10.1074/jbc.M116.770057.
Huveneers S, Truong H, Fässler R, Sonnenberg A, Danen EHJ. Binding of soluble fibronectin to integrin alpha5 beta1 - link to focal adhesion redistribution and contractile shape. J Cell Sci. 2008;121(Pt 15):2452–62. https://doi.org/10.1242/jcs.033001.
Tözeren A, Kleinman HK, Wu S, Mercurio AM, Byers SW. Integrin alpha 6 beta 4 mediates dynamic interactions with laminin. J Cell Sci. 1994;107(Pt 11):3153–63 http://www.ncbi.nlm.nih.gov/pubmed/7699012. Accessed 29 Apr 2019.
Qi Y, Xu R. Roles of PLODs in collagen synthesis and cancer progression. Front Cell Dev Biol. 2018;6:66. https://doi.org/10.3389/fcell.2018.00066.
Fang M, Yuan J, Peng C, Li Y. Collagen as a double-edged sword in tumor progression. Tumour Biol. 2014;35:2871–82. https://doi.org/10.1007/s13277-013-1511-7.
Canty EG, Kadler KE. Procollagen trafficking, processing and fibrillogenesis. J Cell Sci. 2005;118(Pt 7):1341–53. https://doi.org/10.1242/jcs.01731.
Kelwick R, Desanlis I, Wheeler GN, Edwards DR. The ADAMTS (a disintegrin and metalloproteinase with thrombospondin motifs) family. Genome Biol. 2015;16:113. https://doi.org/10.1186/s13059-015-0676-3.
Albrechtsen R, Kveiborg M, Stautz D, Vikeså J, Noer JB, Kotzsh A, et al. ADAM12 redistributes and activates MMP-14, resulting in gelatin degradation, reduced apoptosis and increased tumor growth. J Cell Sci. 2013;126:4707–20. https://doi.org/10.1242/jcs.129510.
Arpino V, Brock M, Gill SE. The role of TIMPs in regulation of extracellular matrix proteolysis. Matrix Biol. 2015;44–46:247–54. https://doi.org/10.1016/J.MATBIO.2015.03.005.
Brew K, Nagase H. The tissue inhibitors of metalloproteinases (TIMPs): an ancient family with structural and functional diversity. Biochim Biophys Acta. 1803;2010:55–71. https://doi.org/10.1016/j.bbamcr.2010.01.003.
Schacht V, Ramirez MI, Hong Y-K, Hirakawa S, Feng D, Harvey N, et al. T1 /podoplanin deficiency disrupts normal lymphatic vasculature formation and causes lymphedema. EMBO J. 2003;22:3546–56. https://doi.org/10.1093/emboj/cdg342.
Ramirez MI, Millien G, Hinds A, Cao Y, Seldin DC, Williams MC. T1alpha, a lung type I cell differentiation gene, is required for normal lung cell proliferation and alveolus formation at birth. Dev Biol. 2003;256:61–72 http://www.ncbi.nlm.nih.gov/pubmed/12654292. Accessed 30 Apr 2019.
Mahtab EAF, Wijffels MCEF, Van Den Akker NMS, Hahurij ND, Lie-Venema H, Wisse LJ, et al. Cardiac malformations and myocardial abnormalities inpodoplanin knockout mouse embryos: correlation with abnormal epicardial development. Dev Dyn. 2008;237:847–57. https://doi.org/10.1002/dvdy.21463.
Fujita N, Takagi S. The impact of Aggrus/podoplanin on platelet aggregation and tumour metastasis. J Biochem. 2012;152:407.
Ogasawara S, Honma R, Kaneko MK, Fujii Y, Kagawa Y, Konnai S, et al. Podoplanin expression in canine melanoma. Monoclon Antib Immunodiagn Immunother. 2016;35:304–6.
Kato Y, Ohishi T, Kawada M, Maekawa N, Konnai S, Itai S, et al. The mouse–canine chimeric anti-dog podoplanin antibody P38B exerts antitumor activity in mouse xenograft models. Biochem Biophys Reports. 2019;17:23–6.
Cueni LN, Chen L, Zhang H, Marino D, Huggenberger R, Alitalo A, et al. Podoplanin-Fc reduces lymphatic vessel formation in vitro and in vivo and causes disseminated intravascular coagulation when transgenically expressed in the skin. Blood. 2010;116:4376.
Hammer AS, Couto CG, Swardson C, Getzy D. Hemostatic abnormalities in dogs with hemangiosarcoma. J Vet Intern Med. 1991;5:11.
Khazenzon NM, Ljubimov AV, Lakhter AJ, Fujita M, Fujiwara H, Sekiguchi K, et al. Antisense inhibition of laminin-8 expression reduces invasion of human gliomas in vitro. Mol Cancer Ther. 2003;2:985–94 http://www.ncbi.nlm.nih.gov/pubmed/14578463. Accessed 4 Apr 2019.
Horibata S, Rice EJ, Mukai C, Marks BA, Sams K, Zheng H, et al. ER-positive breast cancer cells are poised for RET-mediated endocrine resistance. PLoS One. 2018;13:e0194023.
Andrews S. FastQC A Quality control tool for high throughput sequence data. 2010. http://www.bioinformatics.babraham.ac.uk/projects/fastqc/. Accessed 29 Mar 2019.
Chu T, Wang Z, Chou S-P, Danko CG. Discovering transcriptional regulatory elements from run-on and sequencing data using the web-based dREG gateway. Curr Protoc Bioinforma. 2018:e70. https://doi.org/10.1002/cpbi.70.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–60. https://doi.org/10.1093/bioinformatics/btp324.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550. https://doi.org/10.1186/s13059-014-0550-8.
Mi H, Muruganujan A, Ebert D, Huang X, Thomas PD. PANTHER version 14: more genomes, a new PANTHER GO-slim and improvements in enrichment analysis tools. Nucleic Acids Res. 2019;47(Database issue):D419. https://doi.org/10.1093/NAR/GKY1038.
Honma R, Kaneko MK, Ogasawara S, Fujii Y, Konnai S, Takagi M, et al. Specific detection of dog podoplanin expressed in renal glomerulus by a novel monoclonal antibody PMab-38 in immunohistochemistry. Monoclon Antib Immunodiagn Immunother. 2016;35:212–6. https://doi.org/10.1089/mab.2016.0022.
Cherrington BD, Morency E, Struble AM, Coonrod SA, Wakshlag JJ. Potential role for peptidylarginine deiminase 2 (PAD2) in citrullination of canine mammary epithelial cell histones. PLoS One. 2010;5:e11768. https://doi.org/10.1371/journal.pone.0011768.
We would like to acknowledge Dr. Marta G. Castelhano, Ms. Susan J. Garrison (Cornell Veterinary BioBank), Dr. Daniel Fletcher (Emergency and Critical Care, Cornell University Companion Animal Hospital), and Ms. Luce E. Guanzini (Center for Animal Resources and Education, Cornell University) for collecting samples that were used for this manuscript. The authors also acknowledge generous resources provided by Canine Comparative Oncology and Genomics Consortium (CCOGC).
This research was supported by the Judy and Fred Wilpon Family Foundation, Baker Institute internal grant (dog project #3), Cornell Research Animal Health grant (SAC) in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript. We would also like to thank Kyle’s Legacy Inc., Mr. Jake Holshuh, Mrs. Alice C. Sinclair, and ClancysCure. The production of podoplanin antibody was supported by the following grants, JP19am0401013, JP19am0101078, JP19ae0101028 from AMED (YK).
Ethics approval and consent to participate
Not applicable. All samples in this study were purchased from the Canine Comparative Oncology & Genomics Consortium and Cornell Veterinary Biobank.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplemental methods and figure legends.
Supplemental tables for gene ontology analysis.
FS1. Trichrome staining quantification.
FS2. ChRO-seq counts on PDPN and LAMA4 gene bodies.
FS3. PDPN immunohistochemistry quantification.
FS4. LAMA4 immunohistochemistry quantification.
FS5. Masson’s trichrome staining and IHC on serial sections from case B176 (HSA).
FS6. Masson’s trichrome staining and IHC on serial sections from case B554 (HSA).
About this article
Cite this article
Mukai, C., Choi, E., Sams, K.L. et al. Chromatin run-on sequencing analysis finds that ECM remodeling plays an important role in canine hemangiosarcoma pathogenesis. BMC Vet Res 16, 206 (2020). https://doi.org/10.1186/s12917-020-02395-3