Skip to main content

RHAMM regulates MMTV-PyMT-induced lung metastasis by connecting STING-dependent DNA damage sensing to interferon/STAT1 pro-apoptosis signaling

Abstract

Background

RHAMM is a multifunctional protein that is upregulated in breast tumors, and the presence of strongly RHAMM+ve cancer cell subsets associates with elevated risk of peripheral metastasis. Experimentally, RHAMM impacts cell cycle progression and cell migration. However, the RHAMM functions that contribute to breast cancer metastasis are poorly understood.

Methods

We interrogated the metastatic functions of RHAMM using a loss-of-function approach by crossing the MMTV-PyMT mouse model of breast cancer susceptibility with Rhamm−/− mice. In vitro analyses of known RHAMM functions were performed using primary tumor cell cultures and MMTV-PyMT cell lines. Somatic mutations were identified using a mouse genotyping array. RNA-seq was performed to identify transcriptome changes resulting from Rhamm-loss, and SiRNA and CRISPR/Cas9 gene editing was used to establish cause and effect of survival mechanisms in vitro.

Results

Rhamm-loss does not alter initiation or growth of MMTV-PyMT-induced primary tumors but unexpectedly increases lung metastasis. Increased metastatic propensity with Rhamm-loss is not associated with obvious alterations in proliferation, epithelial plasticity, migration, invasion or genomic stability. SNV analyses identify positive selection of Rhamm−/− primary tumor clones that are enriched in lung metastases. Rhamm−/− tumor clones are characterized by an increased ability to survive with ROS-mediated DNA damage, which associates with blunted expression of interferon pathway and target genes, particularly those implicated in DNA damage-resistance. Mechanistic analyses show that ablating RHAMM expression in breast tumor cells by siRNA knockdown or CRISPR-Cas9 gene editing blunts interferon signaling activation by STING agonists and reduces STING agonist-induced apoptosis. The metastasis-specific effect of RHAMM expression-loss is linked to microenvironmental factors unique to tumor-bearing lung tissue, notably high ROS and TGFB levels. These factors promote STING-induced apoptosis of RHAMM+ve tumor cells to a significantly greater extent than RHAMM−ve comparators. As predicted by these results, colony size of Wildtype lung metastases is inversely related to RHAMM expression.

Conclusion

RHAMM expression-loss blunts STING-IFN signaling, which offers growth advantages under specific microenvironmental conditions of lung tissue. These results provide mechanistic insight into factors controlling clonal survival/expansion of metastatic colonies and has translational potential for RHAMM expression as a marker of sensitivity to interferon therapy.

Introduction

Recent advances in early detection, and the development of targeted therapies have significantly improved clinical outcome of breast cancer patients [1,2,3], but metastasis and tumor recurrence remain as major obstacles to successful disease management [4]. Metastasis requires an ability of tumor cells to survive and expand in alien microenvironments and occurs as a result of a complex set of tumor and host traits that are distinct from those facilitating primary tumor initiation and growth. The metastatic process is still poorly understood at a mechanistic level [5, 6] although the discovery of promoter and suppressor genes that selectively affect metastases but not primary tumors has facilitated identifying tumor-intrinsic and microenvironmental properties necessary for secondary site colonization [6, 7].

One oncogenic signaling hub that has been implicated in the metastases of breast and other cancers is the multifunctional, intracellular/extracellular protein, RHAMM (gene name HMMR) [8,9,10,11]. RHAMM expression is heterogeneous in breast cancer, and the presence of strongly RHAMM-positive tumor cell subsets is linked to increased peripheral metastasis and poor clinical outcome [12]. Experimental models of breast cancer confirm a role for RHAMM in promoting functions associated with breast tumor initiation and metastasis [13,14,15] predicting that targeting this protein may improve clinical management of breast cancer metastases. Mechanistically, RHAMM performs multiple extracellular and intracellular functions relevant to metastasis [12, 16] including regulation of mitosis [17], genomic stability [18], cell motility, cellular plasticity [19,20,21], pluripotency of progenitor cells [20, 22] and oncogenic driver pathway activation [13]. In addition to tumor-intrinsic functions, RHAMM also regulates host cell responses that can impact tumor cell survival [8].

To better define the mechanisms by which RHAMM functions affect metastasis, we assessed the consequences of Rhamm-loss to mammary tumor progression using the MMTV-PyMT transgenic mouse model of breast cancer susceptibility [23]. This model was chosen for its rapid progression to metastatic disease, and its molecular similarity to both luminal B breast cancer and basal-like breast tumors, which typically express high levels of RHAMM clinically [24, 25]. Rhamm-loss in this model has no detectable effect on primary tumor initiation or growth but unexpectedly increases, rather than decreases, lung metastasis. This effect is traced to clonal selection of Rhamm−/− tumor cells with an intrinsic resistance to DNA-damage-induced apoptosis that is sensed by STING/interferon signaling. This mechanism provides a survival advantage in lung but not the mammary microenvironment that is linked to the higher ROS and TGFB levels in lung tissue, which enhance STING-dependent apoptosis of RHAMM+ve tumor cells but spare RHAMM−ve comparators. These results identify RHAMM as a novel tissue-specific metastasis regulator and document tumor intrinsic and microenvironmental contexts that trigger its apparent metastasis suppressor functions.

Material and methods

Mouse breeding, genotyping, tumor measurements and whole mount preparations

C57Bl/6 Rhamm−/− and Wildtype mice were crossed to MMTV-PyMT mice on an FVB background (purchased from Jackson Labs) as described by Lopez et al. [26] to obtain Rhamm−/−:MMTV-PyMT, Rhamm±:MMTV-PyMT and Rhamm+/+:MMTV-PyMT genotypes. The degree of SNP homozygosity, measured using the mouse diversity genotyping array (MDGA), was similar between Rhamm−/−:MMTV-PyMT and Rhamm+/+:MMTV-PyMT genotypes. Littermate heterozygotes (Rhamm±) from the Rhamm±:MMTV-PyMT x Rhamm−/−:MMTV-PyMT cross produced a similar tumor profile as the Rhamm+/+:MMTV-PyMT x C57Bl/6 Wildtype cross mice. The preparation and breeding of Rhamm−/− mice, tumor measurements and whole mount preparations are described in Additional file 1: Methods.

Methods for IHC, IF, and immunoblot are described in [27] and in Additional file 1: Methods.

CRISPR cell line generation

The MDA-MB-231 RHAMM CRISPR cell line was generated by transfection with paired guide RNA’s (5′–3′) GTATTGTATTTGATTAGAAT (within exon 3 of the RHAMM gene) and GAATTTGAGAATTCTAAGCT (within exon 6) in plasmid pCR4-TOPO-U6-HPRT-gRNA. Guide RNA’s were co-transfected with plasmid expressing the CAS9 enzyme (pT3.5 Caggs-FLAG-hCas9) as well as plasmids for puromycin and GFP selection, pcDNA-PB7 and pPB SB-CG-LUC-GFP (Puro)(+CRE), using Lipofectamine 2000 reagent (Invitrogen, cat#11668-019) following the manufacturer’s suggested protocol. Mock cell lines were generated by transfection of parent cells with selection plasmids only and selected as a pool by culture in puromycin containing medium (0.6 µg/ml). RHAMM-CRISPR knockout cell lines were selected by clonal plating in puromycin containing media (0.6 µg/ml). Single cell derived colonies were expanded and screened by genomic PCR for the corresponding deletion within the RHAMM gene (primers 5′–3′ AGATACTACCTTGCCTGCTTCA and ACCTGCAGCTTCATCTCCAT), and by immunoblot for loss of RHAMM protein.

Primary cultures of tumor cells

Tumor cell isolation is described in Additional file 1:  Methods.

Cultured cell treatment

For quantification of H2O2-induced apoptosis and STAT1 activation, RHAMM-CRISPR knockout and mock transfected cells were plated on cover slips in DMEM medium containing 10% FBS, resulting in sub-confluent cultures after 24 h incubation. STAT1 activation and apoptosis was induced by incubation in culture medium containing 50–200 µM H2O2 for either 4 (STAT1 activation) or 48–72 h (apoptosis). Cells were stained for either STAT1 or cleaved CASPASE 3 as described in Additional file 1: Methods. Staining quantification by ImageJ used confocal images.

To induce DNA double strand breaks or apoptosis in tumor cells that were isolated from primary tumors, cells were plated at low cell density on fibronectin coated coverslips and cultured in growth medium for 24–48 h. Cultures were treated with either Cis-platin or 300 µM H2O2 at the indicated concentrations and durations. Cells were stained for either γH2AX or ApopTag as described in Additional file 1:  Methods. Staining quantification by ImageJ used confocal images.

Py8119 cells were obtained from the ATCC (ATCC CRL-3278) and cultured in F12K medium (Wisent) supplemented with 5% FBS and Mito + Serum Extender. These cells were originally isolated and cloned from tumors that arose in C57Bl/6 MMTV-PyMT mice and therefore do not contain genomic sequences from FVB mice [28].

siRNA transfection

Sub-confluent Py8119 tumor cells cultures were transfected using Lipofectamine RNAiMAX (Invitrogen) following the manufacturer’s instructions. Culture medium was changed to CTS Opti-MEM (Gibco). Lipofectamine RNAiMAX reagent was diluted 1:50 with CTS Opti-MEM. siRNA (ID: 151008, 159287, s279, negative control siRNA#1, Ambion) was diluted to a concentration of 500 nM in CTS Opti-MEM. Diluted transfection reagent and siRNA were mixed 1:1, incubated at RT for 20 min and then, added to the cell cultures. After 4–5 h incubation at 37 °C, medium was changed to culture medium containing STING agonists (Vadimezan (DMXAA), 33 µM, G10 for MDA-MB-231 tumor cells, 40 µM), H2O2 (300 µM), DMSO and/or TGFB1 (5 ng/ml). STAT1, CASPASE 3 staining or cell survival were analyzed 20–72 h later.

AlamarBlue assay

Py8119 cells were plated at a density of 3000 cells/well of a 96 well plate using complete culture medium. After ON incubation at 37 °C, cells were transfected with Stat1 siRNA, Rhamm siRNA or negative control siRNA. After transfection, cells were treated with F12K culture medium containing STING Agonist Vadimezan (33 µM), DMSO, H2O2 (300 µM), and/or TGFB1 (5 ng/ml). After 72 h incubation at 37 °C, the number of surviving cells was quantified by adding AlamarBlue reagent (1/10 Vol.) followed by 1–2 h incubation at 37 °C. Fluorescence was measured using a plate reader.

DNA and RNA isolation are described in Additional file 1:  Methods.

Analysis of de novo mutation genotypes

The somatic (de novo) mutation burden in Wildtype and Rhamm−/− MMTV-PyMT tumors were compared by identifying germline and mammary tumor mutations using a mouse genomic diversity array (MGDA) unbiased platform [29, 30]. The MDGA detects large, de novo postzygotic deletions and duplications as CNVs using close to 900,000 markers, and also de novo postzygotic base substitutions at single nucleotide polymorphic loci as SNVs, assayed at close to 500,000 SNV loci distributed across the mouse genome. Mouse diversity genotyping array (MDGA) hybridization was performed at the London Regional Genomics Centre (Robarts Research Institute, London, ON) according to instructions in the Affymetrix® Genome-Wide Human single nucleotide polymorphism (SNP) Nsp/Sty 6.0 manual (Affymetrix 2007; https://assests.thermofisher.com/TFS-Assets/LSG/manuals/snp6_atp_userguide.pdf). The resulting CEL files were then used for single nucleotide variant (SNV) genotyping and copy number variant (CNV) identification as described (Additional file 1: Methods). Ilk/Rrp8, Taf10 genes (Mm00232271_cn), which overlap the same CNV region, were used for confirmation by ddPCR. The transferrin receptor gene (Tfrc) was used as the diploid copy number reference for the Ilk, Taf10 and Rhamm assays. No-template controls and two technical replicates were used in all assays. DdPCR procedures are described in Additional file 1: Methods. Recurrent and unique CNVs were determined using HD-CNV (Additional file 1: Methods) [31]. To visualize CNV occurrence across the genome for individual mouse samples, a timeline-style plot was generated in R using ggplot2 (v3.2.1). Genic annotation used to identify the genic content of CNVs was obtained from Ensembl’s BioMart (Ensembl genes 67, NCBIM37). Protein-coding genes, non-coding genes, and pseudogenes that completely overlapped CNVs of the same state, in all three samples of a group (shared Rhamm genotype and tumor type), were considered recurrent within that group.

Phenogram construction was done as previously described [29] and is described in Additional file 1: Methods. Filtering procedures are described in Additional file 1: methods. Identification of candidate de novo mutations is described in Additional file 1: Methods. The spatial distribution of candidate de novo mutations across the genome was visualized using rainfall plots (Additional file 1: Methods) [32]. Chromosomes were treated linearly, and the genomic position of de novo mutations was assigned in an additive manner based on the position of the locus in relation to the whole GRCm38.p4 genome.

RNA sequencing is described in Additional file 1: Methods.

Pathway analysis

Gene expression differences between Rhamm−/− and Wildtype tumors of 1.5-fold with a p-value of less than 0.05 were analyzed for enriched Gene Ontology (GO) and KEGG pathways. The functions of differentially expressed genes were further probed using the hallmark gene set from the Molecular Signatures Database v 7.1 (https://www.gsea-msigdb.org/gsea/msigdb/index.jsp), used in the GSEA analysis, and Metascape (https://metascape.org/gp/index.html#/main/step1). Significantly down-regulated genes that grouped into the top hallmark gene sets in the GSEA analysis were assessed for mRNA co-expression (p < 0.05) with RHAMM (HMMR) using breast cancer data sets in cBioPortal (cbioportal.org). Invasive breast cancer subtypes that express the highest RHAMM mRNA levels were identified using the molecular subtypes in the TCGA PanCancer breast invasive carcinoma data set in cBioPortal (cbioportal.org).

Apoptag® staining, ROS/NOS and 8-oxodG ELISA are described in Additional file 1: Methods.

Results

Rhamm-loss increases metastasis in MMTV-PyMT mice but does not affect primary tumor initiation or growth

Data set analyses (cBioportal.org) confirm a previous report that RHAMM expression is highest in basal and luminal B breast cancer subtypes [33] (Fig. 1A, TCGA PanCancer Atlas, cBioportal). We utilized the MMTV-PyMT transgenic mouse, which models aggressive luminal B breast cancer [24, 34] to uncover oncogenic functions of RHAMM in an immune intact microenvironment, using loss-of-function, germline ablation of Rhamm [35].

Fig. 1
figure 1

Rhamm-loss in MMTV-PyMT transgenic mice increases lung metastasis without affecting primary mammary tumor initiation, incidence or size. A Comparison of RHAMM mRNA expression in breast cancer subtypes (Breast Invasive Carcinoma, TCGA, PanCancer Atlas). RHAMM expression is significantly higher in basal and luminal B breast cancer subtypes than in luminal A comparators. B Rhamm-loss detected by the Mouse Diversity Genotyping Array (MDGA) genotyping was confirmed by ddPCR, which is the standard ultrasensitive method for verification of genotyping array-based CNV detection. Two additional ddPCR assays confirm another CNV deletion overlapping Ilk and Taf10 genes that is detected by MDGA CNV genotyping, confirming accuracy of the MDGA and consistency of the two detection methods. WT = Wildtype; Rh−/− = Rhamm−/−. C Western blot assays were performed using primary tumor lysates and RHAMM antibodies to the N-terminal sequence as described in Methods. Full-length RHAMM is expressed in Wildtype tumors (arrow) and is not detected in Rhamm−/− tumors. D Histology sections of primary tumors from 16-week-old mice. Tissue was paraffin-processed, stained for RHAMM protein and imaged with a Nikon confocal microscope as described in Methods. Results show heterogeneous RHAMM staining in primary tumors and also in the host microenvironment (arrows). The lack of staining in histology sections from Rhamm−/− tissue confirms the specificity of the anti-RHAMM antibody. E Whole mounts of mammary fat pads were prepared as described in Methods and show primary tumor initiation at 5 weeks in both Wildtype and Rhamm−/− mice. Detectable primary tumor masses of Wildtype and Rhamm−/− mice were measured weekly with calipers. Differences between genotypes are not statistically significant. Values are the Mean and S.E.M. n = 15 mice. F The lungs of Wildtype and Rhamm−/− mice were harvested and metastatic nodules identified in hematoxylin/eosin-stained histology sections. 100% of Rhamm−/− mice contain metastatic nodules in lungs in contrast to 58% of Wildtype mice, n = 11 mice/genotype. G Metastatic nodules were quantified from serial histology sections as described in Methods. Rhamm-loss significantly increases the number of metastatic colonies/lung compared to Wildtype. Values are the Mean and S.E.M. n = 15 mice. *p < 0.05

For these experiments, we followed the breeding strategy of Lopez et al. [26] that resulted in Cd44−/−:MMTV-PyMT mice on a mixed C57Bl/6:FVB background. This study showed that Cd44-loss increased lung metastases, which was proposed to result from functional compensation by Rhamm [36]. We confirmed that Cd44-loss on this mixed background increases the incidence of lung metastasis (Additional file 2: Fig. S1). Successful deletion of Rhamm in Rhamm−/−:MMTV-PyMT mice was confirmed using the Mouse Diversity Genotyping Array (MDGA) and ddPCR (Fig. 1B), while absence of RHAMM protein expression was verified in Rhamm−/− tumors by Western blot and immunohistochemistry (Fig. 1C, D). RHAMM full-length protein was confirmed to be expressed in Wildtype but not Rhamm−/− MMTV-PyMT tumors (Fig. 1C), and immunofluorescent analyses further show that RHAMM protein expression is heterogeneous in Wildtype primary tumors, similar to human breast cancers [37], and also occurs in host cells of the tumor microenvironment (Fig. 1D).

As previously documented [38], Wildtype MMTV-PyMT primary tumors initiate in the mammary glands by 5 weeks (Fig. 1E) and reach a palpable size by 8 weeks (Fig. 1E, graph), while lung metastases are detected visibly as nodules and by histology at 14–16 weeks (Fig. 1F). Rhamm-loss does not detectably alter primary tumor initiation, incidence or growth (Fig. 1E) but results in increased lung colonies (Fig. 1F, G). Thus, 100% of Rhamm−/− MMTV-PyMT mice develop lung tumors compared to 56% of Wildtype comparators (Fig. 1F), and metastatic colony number/lung is also significantly increased with Rhamm-loss (Fig. 1G). Dual loss of Rhamm and Cd44 does not modify primary tumor initiation or growth, and these mice also exhibit elevated metastases similar to the single knockout of Cd44 or Rhamm (Additional file 2: Fig. S1).

Known RHAMM functions do not account for increased lung metastasis

RHAMM is well-documented to regulate a number of cellular processes in vitro that potentially affect metastasis including regulating tumor cell motility, invasion, cellular plasticity, proliferation and genomic instability [8, 15]. Analyses of primary Rhamm−/− MMTV-PyMT tumor cell cultures show that Rhamm-loss significantly reduces both directed migration (scratch wound assays) and random motility speed (Timelapse assays) (Additional file 2: Fig. S2A, B). Rhamm-loss has no detectable effect on invasion, using standard Boyden chamber assays (Additional file 2: Fig. S2C), tumor cell proliferation (tumor size, Ki67 IHC, Additional file 2: Fig. S2D, E) or epithelial-to-mesenchymal (EMT) plasticity. Evidence for EMT was assessed by RNA-seq and immunofluorescence. RNA-seq analyses reveal no significant change in the expression of commonly used markers [39] for EMT (CDH1, CDH2, SNAI2, TWIST1, and FN1), while immunofluorescence analyses show no difference in protein staining for VIM and ZEB1 (Additional file 2: Fig. S2F).

As expected, and independent of genotype, the total de novo mutation burden in MMTV-PyMT tumors is increased and CNVs are longer, have more copy number losses than gains, and occur more frequently within genes than the inherited germline CNV burden (Fig. 2A, Additional file 2: Tables S1, S2). The total mutation burden (Additional file 2: Table S1, S2, S3) and its genome-wide distribution in primary and metastatic tumors (Additional file 2: Fig. S3A) are not significantly (e.g., p > 0.05) altered by Rhamm-loss. There is also no obvious difference between genotypes in double- and single-strand DNA breaks as detected by Comet assays in Wildtype and Rhamm−/− primary mammary tumor cells maintained in standard culture conditions (Fig. 3A) although CNV’s are uniquely shorter in Rhamm−/− primary tumors (Additional file 2: Fig. S3B). Collectively, these results rule out obvious changes in motility, proliferation, plasticity and increased genomic instability associated with unrepaired DNA breaks and total mutation burden as mechanisms for the enhanced lung colonization observed in Rhamm−/−:MMTV-PyMT mice.

Fig. 2
figure 2figure 2

Wildtype and Rhamm−/− tumor mutation burden is similar, but Rhamm−/− tumors uniquely exhibit a strong inter-animal homogeneity in mutation characteristics. A The number of CNV gains and losses detected in primary tumor and metastases from Wildtype and Rhamm−/− mice (n = 3/group) were compared with two inbred mouse stocks, C57Bl/6J (n = 8) and FVB/NJ (n = 1) and a reference set of inbred mice used for comparison with the inherited CNV burden (n = 114) as described in Methods. CNV gains are shown in red, while losses are shown in blue. Statistically significant differences between genotypes are not detected. Error bars represent standard error. B. Variation in the properties of de novo CNVs and SNVs in Wildtype metastases (WM), Wildtype primary (WP) tumors Rhamm−/− metastases (RM) and Rhamm−/− primary tumors (RP), (n = 3/group) is displayed as a principal component analysis (PCA). Results show clustering of these properties in Rhamm−/− but not Wildtype tumors, indicating that the two genotypes differ in the CNV and SNV properties. C, D Phenograms representing SNV (C) and CNV (D) genetic differences. The numeric values in the sample labels represent the individual mouse identifier number. Scale bars represent genetic differences between samples. WP = Wildtype primary tumor; WM = Wildtype lung metastases; RP = Rhamm−/− primary tumor; RM = Rhamm−/− lung metastases. SNV and CNV differences were calculated as described in Methods. SNVs of Rhamm−/− tumors clearly segregate from Wildtype comparators (C). Clear segregation of CNV differences between genotypes is not observed (C). E Different colors represent different mutation types: C:G > A:T (red), C:G > G:C (blue), C:G > T:A (green), T:A > A:T (purple), T:A > C:G (orange), T:A > G:C (yellow). Rhamm−/− rainfall plots for both primary and metastatic tumors are similar between all three mice, while Wildtype rainfall plots vary, showing that a strong similarity in the pattern of SNV is acquired with Rhamm-loss. Results predict positive selection of a limited number of Rhamm−/− tumor clones. F. SNV burden that is unique to Wildtype primary tumors, Wildtype lung metastases, Rhamm−/− primary tumors and Rhamm−/− lung metastases. There is a greater number of lung-specific SNVs in Rhamm−/− metastasis vs. primary tumors, but these are unchanged in Wildtype tumors. Values are the Mean and S.E.M. n = 3 mice/group. **p < 0.01

Fig. 3
figure 3

Rhamm−/− tumor cells are resistant to DNA damage. A DNA fragmentation was quantified by comet assays as a measure of DNA damage using frozen samples of isolated, unstimulated Wildtype and Rhamm−/− primary tumor cells. Quantification of % DNA in tail and tail moment is not affected by Rhamm-loss consistent with similar mutation burdens in the two genotypes identified by the genotyping array. Values are the Mean and S.E.M of n = 3 cell samples/genotype and n = 21 wells/genotype. B The amount of ROS-induced DNA damage in Wildtype and Rhamm−/− tumors was quantified using ELISA as described in Methods. Results show that DNA from lung metastases of both genotypes carry more ROS-induced DNA damage that can cause point mutations and CNVs than primary tumors, which is consistent with the higher levels of ROS in lungs versus mammary fat pads. Values are the Mean and S.E.M. N = 4 biological replicates and 3 technical replicates. *p < 0.05, ***p < 0.001. C, D Cultures of Wildtype and Rhamm−/− primary tumor cells were exposed to H2O2 (C) and cisplatin (D) to induce DNA damage. DNA damage was detected by gamma-H2AX staining, and cell death was quantified by Apoptag reactivity as described in Methods. Rhamm−/− tumor cells survive with more DNA damage than Wildtype comparators. Values are the Mean and S.E.M. n = 3 biological replicates *p < 0.05, **p < 0.01

Rhamm-loss alters the nature of the SNV mutation burden

A principal component analysis (PCA) of the collective de novo CNV and SNVs in Wildtype and Rhamm−/− MMTV-PyMT primary and metastatic tumors was constructed to identify additional mutation patterns that might differ between the two genotypes (Fig. 2B) to provide clues for the increase in lung colonies resulting from Rhamm-loss. For CNVs, the number of gains, losses, average length of losses, average length of gains, number of losses overlapping or encompassing a gene and number of gains overlapping or encompassing a gene were assessed. SNV metrics included in the PCA are total mutations and number of de novo SNV genotype changes unique to each group (Additional file 2: Table S3). The PCA predicts that genetic variation in these collective properties is reduced by Rhamm-loss in both primary and metastatic tumors. A phenogram comparing genetic properties between tumors of the two genotypes was calculated using pairwise mutations at CNV and SNV loci. This shows that SNV properties clearly segregate Rhamm−/− from Wildtype tumors (Fig. 2C) while, in contrast, CNV properties do not (Fig. 2D). These results predict that recurrent SNVs but not CNVs are associated with the metastatic progression as has previously been demonstrated in other pre-clinical models of breast cancer [40]. These results prompted us to further analyze the nature of SNVs in Rhamm−/− vs. Wildtype tumors and their relationship to the increase in lung metastasis resulting from Rhamm-loss.

Graphic visualization of individual SNVs relative to the genome position (rainfall plots) provides an overview of mutation distribution and allows for pattern recognition (random vs. rainfall). Rainfall plots of SNVs show that the genome-wide distribution is similar in both genotypes (Fig. 2E), and that inter-tumor SNV heterogeneity is high in Wildtype tumors. However, SNV pattern heterogeneity is remarkably decreased in both the primary and metastatic tumors of Rhamm−/− mice (Fig. 2E). Furthermore, while the total SNV burden and number of shared somatic SNV are similar in primary and lung Wildtype tumors, they are strongly increased by 100-fold in Rhamm−/− lung vs. primary tumors (Fig. 2F, Additional file 2: Tables S1, S3). The majority of these shared de novo SNVs are non-genic (Additional file 2: Table S3).

These collective genetic analyses demonstrate that primary Rhamm−/− and Wildtype primary tumors are significantly different from each other despite their similar growth and initiation properties in mammary tissue. The homogeneity of somatic SNV patterns predicts that total animal Rhamm-loss exerts positive or purifying selection for tumor cell clones within the mammary fat pad. We hypothesized that while these clones do not have a growth advantage in the mammary gland, their enrichment in lung metastases provide a selective advantage for survival in a lung microenvironment. Since the microenvironment of the lung is subject to greater oxidative stress than other tissues, the possibility that Rhamm−/− tumor cells might be able to survive with more ROS-induced DNA damage than Wildtype comparators was next assessed.

Rhamm-loss de-sensitizes tumor cells to DNA damage

Differences in oxidative stress of mammary vs. lung tissues were assessed by quantifying 8-Oxo-DG, used as a marker for oxidative damage (Fig. 3B). Rhamm−/− primary tumors and lung metastases have significantly higher levels of 8-Oxo-DG than Wildtype comparators. While the comet assay did not reveal differences in single or double strand DNA breaks (Fig. 3A), the primary tumor cells used in this assay were not exposed to ROS-inducing agents. We therefore next compared the tolerance of primary Rhamm−/− and Wildtype tumor cells to DNA damage in vitro by exposing primary cultures of tumor cells to ROS (H2O2, Fig. 3C) or Cisplatin (Fig. 3D), both of which result in DNA damage and SNVs as well as CNVs. Rhamm−/− tumor cells exposed to H2O2 or Cisplatin contain more DNA damage than Wildtype tumor cells, as detected by gamma H2AX staining. Rhamm−/− tumor cells are also significantly more resistant to DNA-damage-induced apoptosis than Wildtype counterparts (Cisplatin shown, Fig. 3E), indicating that Rhamm-loss provides an intrinsic ability to survive DNA damage. The difference in this ability does not associate with altered expression of DNA repair enzymes, as detected by PCR arrays and RNA-seq, or with functions of shared de novo genic SNVs of Rhamm−/− lung metastases predicted by IPA pathway analyses (Additional file 2: Fig. S4). Further analysis using RNA-seq was therefore performed to uncover mechanisms responsible for the increased ability of Rhamm−/− tumor cells to survive with DNA damage.

Rhamm-loss blunts interferon signaling

RNA-seq shows that Rhamm-loss decreases expression of 459 genes and increases expression of 168 genes in primary mammary tumors (p < 0.05). Unbiased pathway analyses of significantly (p < 0.05) upregulated genes using the Molecular Signatures Database (v7.1) and Metascape reveal a weak association with aberrant G2M checkpoint (3.03e−4), E2F (1.69e−3), MYC signaling (1.06e−2) and cell cycle (Additional file 2: Fig. S5B), respectively. Metascape analyses of down-regulated genes predict strongly suppressed interferon signaling (e.g., Fig. 4A), and hallmark gene set analyses (GSEA) reveal significant down-regulation of gene sets in Rhamm−/− tumor cells that respond to inflammatory stimuli (Additional file 2: Table S4). The top gene sets are regulated by IFNG (p = 2.43e−61) and IFNA (p = 4.54e−39) (Additional file 2: Table S4). Expression of IFNG and IFNA cytokines are not downregulated, but the expression of 53 (IFNG) and 31 (IFNA) pathway and target genes are significantly reduced (approximately twofold) by Rhamm-loss (p < 0.05). Twenty-two of these genes are coordinately expressed with RHAMM mRNA in 2 breast cancer data sets (METABRIC, TCGA Cancer Atlas cBioportal, Fig. 4B). Notably the breast cancer tumor suppressor and interferon-regulated transcription factor STAT1 [41] is most strongly co-expressed with RHAMM in these data sets (Fig. 4B, C, TCGA pan-cancer atlas shown). STAT1 protein levels are significantly reduced in Rhamm−/− MMTV-PyMT primary tumors as detected by Western blots (Fig. 4D–F). Immunofluorescence assays show that STAT1 protein is reduced in both tumor and host cells of Rhamm−/− primary tumors and lung metastases (Fig. 4E). Silencing RHAMM expression by siRNA knockdown in Py8119 MMTV-PyMT tumor cells also reduces STAT1 protein expression (Western blot, Fig. 4F).

Fig. 4
figure 4

Rhamm-loss results in reduced expression of interferon alpha and gamma pathway components. A The transcriptome of primary MMTV-PyMT tumors was analyzed using RNA-seq as described in the Methods. Analysis of genes that are differentially expressed in Wildtype versus Rhamm−/− primary tumors using Metascape and GSEA, identified interferon signaling as the top down-regulated pathway in Rhamm−/− primary tumors. B GSEA Hallmark gene set analyses identified 51 IFNG and IFNA regulated genes as top hallmark gene sets that are down-regulated by Rhamm-loss. Twenty-six of these genes are significantly (p < 0.001) co-expressed with RHAMM in 2 datasets of breast cancer patients (METABRIC and TCGA Pan-Cancer Atlas, c-Bioportal.org) and are shown with Spearman correlation coefficients from METABRIC data sets in the histogram. C STAT1 is an example of an interferon pathway gene that is linearly co-expressed with RHAMM in luminal B breast cancer subtype. Graph is derived from the TCGA Pan-Cancer Atlas data set. D STAT1 protein levels in Wildtype and Rhamm−/− primary mammary tumors were analyzed by Western blot. Wildtype tumors express significantly more STAT1 protein compared to Rhamm−/− tumors. Values are the Mean and SEM of n = 3 tumors. *p < 0.05. E STAT1 staining of Wildtype and Rhamm−/− primary tumors and lung metastases. Wildtype tumors contain more STAT1-positive cells in tumors and the tumor microenvironment than Rhamm−/− counterparts. F Rhamm mRNA expression was knocked down by siRNA in the Py8119 MMTV-PyMT mammary tumor cell line as described in Methods, and the effect on RHAMM and STAT1 protein expression was quantified using Western blots. Rhamm knock-down strongly reduces STAT1 protein expression

Rhamm-loss impairs STING-IFN signaling that senses DNA damage

IFN signaling is complexly associated with experimental and clinical tumorigenesis [42]. Robust IFN signaling is cytotoxic and linked to improved outcome in both luminal B and TNBC/basal-like breast cancer molecular subtypes [43, 44]. However, chronic, low activation of IFN signaling provides pro-survival advantages to cancer cells [45]. Notably, a subset of IFN-regulated genes that constitute an IFN-related DNA damage signature (IRDS), which includes STAT1[46], provides protection against DNA damage sensed by cGAS-STING in breast cancer cell lines [47]. Eighteen of these IRDS genes are regulated by RHAMM expression and are co-expressed with RHAMM in breast cancer datasets (Additional file 2: Table S5). Since blocking IFN signaling increases metastasis but not primary tumor growth in MMTV-PyMT mice [48] and forced RHAMM expression activates cGAS-STING in a BRCA1 mouse model of breast cancer susceptibility [14], we assessed if Rhamm−/− cells have an aberrant cGAS-STING-IFN signaling, which provides a survival advantage in lung tissue. The effects of Rhamm-loss on STING-activated interferon signaling and tumor cell survival were further investigated using MMTV-PyMT tumor cells and the human basal-like MDA-MB-231 breast tumor cell line.

The consequence of RHAMM knockdown by siRNA to STING-mediated MMTV-PyMT tumor cell line survival was assessed using the STING agonist DMXAA. DMXAA significantly reduces RHAMM+ve tumor cell survival but not RHAMM−ve tumor cells (Fig. 5A). Deletion of genomic RHAMM in the human basal-like MDA-MB-231 breast tumor cell line ablates RHAMM protein expression (Tarullo SE, 2023. J. Pathol. 2023 Jul;260(3):289–303) and significantly reduces activation of STAT1 by the STING agonist, G10, as detected by decreased STAT1 nuclear localization relative to RHAMM+/+ mock-transfected cells (Fig. 5B). Further, G10 significantly decreases survival of RHAMM+/+ but not RHAMM−/− MDA-MB-231 cells as detected by AlamarBlue (Fig. 5C) and cleaved CASPASE3 (Fig. 5D). These results identify an impaired response to STING-induced cell death as an intrinsic property of Rhamm−ve mammary tumor cells.

Fig. 5
figure 5

Rhamm-loss blunts response to STING-activated interferon signaling. A Py1189 MMTV-PyMT tumor cells were exposed to the mouse STING agonist, MDXAA, and the effect on tumor cell survival was quantified with AlamarBlue as described in Methods. MDXAA reduces the survival of Rhamm+/+ Py1189 tumor cells but has no significant effect when Rhamm is knocked down. Values are the Median and SD, n = 8 replicates. **p < 0.01, ****p < 0.00001, ******p < 0.000001. B STAT1 activity in the STING agonist, G10-stimulated RHAMM+/+ and RHAMM−/− MDA-MB-231 tumor cells were measured by quantifying nuclear STAT1 protein as described in Methods. Results show that RHAMM-loss significantly reduces STAT1 activity. Values are the median and SD. n = 30 cells/genotype. ****p < 0.0001. C, D RHAMM+/+ and RHAMM−/− MDA-MB-231 tumor cells were exposed to the human STING agonist G10, and survival was quantified using AlamarBlue (C) and cleaved Casapse3 (D) as described in Methods. Results show that similar to the MMTV-PyMT tumor cells, RHAMM-loss in MDA-MB-231 cells generates resistance to G10-induced cell death. Values are the Median and Std. Dev. n = 6 replicates. *p < 0.05, **p < 0.01

The Rhamm −/− lung microenvironment contains high levels of ROS and TGFB

We next compared lung and mammary fat pad microenvironments to identify properties that might provide a selective survival advantage for STING-resistant Rhamm−/− tumor cells in lung tissue. ROS/NOS exposure results in DNA damage, stimulation of cGAS-STING, which senses this damage, and a consequent promotion of IFN-mediated cell apoptosis [42, 49]. We therefore confirmed that exposure to H2O2 at levels that damage DNA (e.g., Fig. 3D), activate interferon signaling as detected by increasing STAT1 nuclear localization in RHAMM+/+ tumor cells (Fig. 6A). ROS/NOS levels (Fig. 6B) were next shown to be significantly higher in lung than in mammary tissue of both genotypes (Fig. 6B). Consistent with RNA-seq data, IFN protein levels are collectively similar in the primary and metastatic tumors of both Rhamm−/− and Rhamm+/+ (Wildtype) genotypes (Fig. 6C, primary tumors shown). Since elevated ROS levels are not specific to the Rhamm−/− genotype, we examined other properties of the lung microenvironment that might contribute to the enhanced metastasis observed in Rhamm−/− lungs. We compared the protein levels of TGF-beta1 in Wildtype and Rhamm−/− tumor-bearing lungs since this cytokine is a metastasis promoter in breast cancer [50, 51] that is linked to RHAMM functions, DNA damage, and regulation of interferon signaling [52]. TGFB mRNA expression does not differ in primary Rhamm−/− vs. Wildtype tumors as detected by RNA-seq but protein levels are significantly increased in tumor-bearing Rhamm−/− lungs relative to Wildtype counterparts (Fig. 6D). Collectively, these results reveal ROS levels as a key difference in lung vs. mammary fat pad microenvironments that likely affect lung colonization and further identify an elevated expression of the immune-suppressive TGFB cytokine specifically detected in Rhamm−/− lung tissue that might favor increased growth and survival of tumor cells with RHAMM expression-loss.

Fig. 6
figure 6

The tumor microenvironments of lung and mammary fat pads differ, but are minimally altered by Rhamm-loss. A MDA-MB-231 RHAMM+/+ tumor cells were exposed to H2O2, and the effect on STAT1 activation was quantified as the number of cells/field with nuclear STAT1 immunofluorescence. H2O2 significantly increases STAT1 activation. Values are the Mean and S.E.M. n = 5 fields. B Reactive oxygen species (ROS/NOS) levels in Wildtype and Rhamm−/− primary and metastatic tumors were quantified using ELISA as described in Methods. Results show that levels are similar in tumors of both genotypes. Notably ROS levels are higher in lung vs. primary tumors from both genotypes. Values are the Mean and S.E.M. of n = 4 biological replicates and 3 technical replicates *p < 0.05, **p < 0.01. C IFNG protein levels were measured by antibody array in Wildtype and Rhamm−/− lungs with metastases as described in Methods. IFNG protein is highly variable in both genotypes but not significantly different. Values are the median and Std. Dev. n = 3 replicates. D TGFB protein was quantified by IHC in primary and lung tumors. TGFB levels are significantly higher in Rhamm−/− lungs than in Wildtype counterparts. Values are the Mean and S.E.M. of n = 5 lungs/genotype

Rhamm −/− tumor cells are resistant to ROS and TGFB-induced cell death

We simulated a ROS/TGFB-high lung microenvironment by exposing RHAMM+ve and RHAMM−ve MMTV-PyMT tumor cells to combinations of ROS, TGFB and STING agonists then quantified the consequence to tumor cell survival. The survival of RHAMM+ve and RHAMM−ve tumor cells is not differentially affected in standard culture conditions of low ROS and TGFB levels (Fig. 7A). In contrast, exposure to TGFB combined with either H2O2 or STING agonist DMXAA significantly decreases the survival of RHAMM+ve but not RHAMM−ve tumor cells (Fig. 7B, C). Thus, tumor cells with RHAMM expression-loss are impaired in their ability to sense DNA damage, which is associated with an impaired response to STING-agonists and with a down-regulation of IFN pathway and target genes. This aberrant response provides a survival advantage in a high ROS and TGFB microenvironment.

Fig. 7
figure 7

Rhamm-loss alters the response of tumor cells to the lung microenvironment. AC The PyMT 8119 cell line was exposed to culture medium alone, TGFB alone or TGFB together with DMXAA or H2O2 then cell survival quantified as described in Methods. A Survival of RHAMM+ve and RHAMM−ve tumor cells is similar in culture medium alone. B, C The addition of TGFB with (B) DMXAA or (C) H2O2 reduce RHAMM+ve but not RHAMM−ve tumor cell survival. The line indicates the Mean and S.E.M. of n = 8 samples. *p < 0.05, ****p < 0.0001

These collective results predict that survival of wildtype lung colonies would also benefit by low or absent Rhamm expression, although not to the extent of a Rhamm−/− lung microenvironment. To assess this possibility, tumor-bearing, Wildtype lung tissue sections were stained for RHAMM and STAT1. As shown in Fig. 8, RHAMM protein expression is low in lung colonies compared to primary tumors and STAT1 protein expression is also blunted (Fig. 4D) indicating impaired IFN signaling, consistent with the evidence that RHAMM expression-loss provides protection from STING-IFN induced apoptosis. These collective data identify a tumor-intrinsic, RHAMM-dependent mechanism for promoting tumor cell survival that favors metastatic lung colonization despite a DNA damage burden.

Fig. 8
figure 8

RHAMM expression is heterogeneous in wildtype primary tumors of MMTV-PyMT mice but is reduced in lung metastases. A, B Wildtype primary MMTV-PyMT tumors express RHAMM in a heterogeneous manner. CF show the varying levels of RHAMM expression in lung colonies of wildtype mice. RHAMM expression is highest in small lung colonies (C), and low (D) to negative (E, F) in larger colonies. RHAMM is consistently expressed in the lung microenvironment (CF)

Discussion

Our results identify a novel tumor-intrinsic mechanism resulting from Rhamm-loss that provides a survival advantage in the lung microenvironment and that associates with clonal dominance. We identify specific properties of the lung microenvironment (high ROS, TGFB levels), which combine with a Rhamm-dependent tumor-intrinsic impairment of DNA damage sensing to provide a selective survival/growth advantage in this tissue but not in the mammary gland microenvironment.

Advanced genomic and diagnostic technologies identify metastasis as an extremely complex process that can result from single or multiple clones arising early or late in the genetic evolution of primary tumors [53]. Tumor clonal genetic heterogeneity is considered to be a driving force in both successful metastatic colonization and treatment resistance. Each step in the colonization of extraneous tissues exerts selective pressure, which results in a genetic and epigenetic heterogeneity that is distinct from the primary tumor. However, the factors driving this selection process are still poorly understood. Therapeutic strategies remain largely based upon analyses of primary tumors [53], and this knowledge gap contributes to the present restricted ability to eradicate and/or manage metastases. Our results identify one mechanism, loss of RHAMM signaling, that associates with reduced genetic diversity and, paradoxically, clonal amplification of lung colonies that are resistant to ROS-mediated DNA damage.

Positive selection, detected by SNV similarities, is clearly evident in Rhamm/− primary mammary tumors even though phenotypic properties are not detectably altered. Selection is strongly amplified in lung tissue metastatic tumors, as indicated by an almost 100-fold enrichment in shared SNVs, which predicts that the lung microenvironment creates a bottleneck limiting expansion of metastatic clones. The enhanced tumor-intrinsic survival capability in the lung is linked to loss of RHAMM expression and blunting of STING/IFN signaling in both Rhamm−/− and Wildtype genotypes. Interestingly, RHAMM expression-loss in Wildtype lung metastatic colonies is not associated with the increased SNV homogeneity observed in Rhamm−/− lung tumor colonies, suggesting that aberrant STING-IFN survival responses alone are not responsible for clonal homogeneity. We analyzed genic SNVs for clues as to mechanisms that might complement this Rhamm−/− intrinsic tumor cell phenotype. It is remarkable that while Wildtype lung metastases share few genic SNVs (Additional file 2: Table S3), and these mutant genes do not group onto the same signalling pathways (data not shown), Rhamm−/− lung metastases collectively share an extraordinary 125 genic SNVs. None of these mutations associate with IFN or TGFB-regulated signalling as queried by KEGG, IPA or GSEA data sets (IPA shown, Additional file 2: Table S6), but in silico analyses predict a role for these mutations in cancer and wound repair. Wound and cancer microenvironments share many immune and stromal fingerprints, suggesting that critical changes in the microenvironment of Rhamm−/− primary tumors may contribute to positive selection of homogeneous clones. Another consideration is evidence that RHAMM expression contributes to maintenance of embryonic stem cell pluripotency [20] and is a marker for a subpopulation of renewing tumor stem cells [54,55,56], which are thought to drive tumor cell heterogeneity. Rhamm-loss may increase clonal homogeneity, in part, by restricting pluripotency of these tumor cell subsets. Further multiplexed and spatial platform analyses are required to identify the tumor-intrinsic and host mechanisms that are regulated by RHAMM, and that underlie the clonal selection and expansion of mammary tumor cells.

Evidence is emerging that the oncogenicity of many dominant oncogenic driver genes is context- and tissue-dependent [36, 48, 49]. For example, high expression of CD44, like RHAMM, is linked to breast cancer initiation yet both proteins can perform metastasis suppressor functions in experimental models of breast cancer susceptibility [26]. Both genes are subject to alternative splicing of mRNA, and selective expression of isoforms can have different tumorigenic effects. For example, expression of one of two RHAMM isoforms is oncogenic in an experimental model of pancreatic cancer [9]. Our results probe the duality in oncogenic functions of RHAMM by providing mechanistic insight into the conditions favoring the metastasis suppressor over promoter functions using RHAMM as an example. This type of contextual and mechanistic information is critical for successful targeting of many multifunctional proteins typified by context-dependent oncogenic and suppressor properties [6, 57]. Thus, while a large body of data predict that RHAMM acts as a tumor promoter in breast cancer [13, 15], our results uncover an unexpected role for RHAMM in suppressing metastasis under the very specific conditions of whole animal Rhamm-loss. In addition to identifying and providing context for the dual oncogenic and tumor suppressor functions of RHAMM, our results raise the possibility of using RHAMM expression as a biomarker for sensitivity to interferon therapy, which, as an example, can impact immune therapy effectiveness in tumor management [58, 59].

In summary, our conclusions from the present study are that RHAMM contributes to primary tumor and metastasis clonal heterogeneity by regulating tumor cell survival via STING/IFNG/STAT1 signaling in microenvironments characterized by a high potential for DNA damage. When RHAMM is lost, these properties do not detectably impact primary tumor initiation or growth but favor survival and growth of homogeneous metastatic clones in lung tissue that are under positive selection as a result of host Rhamm-loss.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

CNV:

Copy number variant

ddPCR:

Digital droplet PCR

GO:

Gene ontology

KEGG:

Kyoto encyclopedia of genes and genomes

MDGA:

Mouse diversity genotyping array

PCA:

Principal component analysis

PCR:

Polymerase chain reaction

RHAMM:

Receptor for HA mediated motility

SNP:

Single nucleotide polymorphism

SNV:

Single nucleotide variant

References

  1. Adams S, Gatti-Mays ME, Kalinsky K, Korde LA, Sharon E, Amiri-Kordestani L, et al. Current Landscape of Immunotherapy in Breast Cancer: A Review. JAMA Oncol. 2019;5(8):1205–14.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Garufi G, Palazzo A, Paris I, Orlandi A, Cassano A, Tortora G, et al. Neoadjuvant therapy for triple-negative breast cancer: potential predictive biomarkers of activity and efficacy of platinum chemotherapy, PARP- and immune-checkpoint-inhibitors. Expert Opin Pharmacother. 2020;21(6):687–99.

    Article  CAS  PubMed  Google Scholar 

  3. Liao Z, Chua D, Tan NS. Reactive oxygen species: a volatile driver of field cancerization and metastasis. Mol Cancer. 2019;18(1):65.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Nicolini A, Barak V, Biava P, Ferrari P, Rossi G, Carpi A. The use of immunotherapy to treat metastatic breast cancer. Curr Med Chem. 2019;26(6):941–62.

    Article  CAS  PubMed  Google Scholar 

  5. Chambers AF, Groom AC, MacDonald IC. Dissemination and growth of cancer cells in metastatic sites. Nat Rev Cancer. 2002;2(8):563–72.

    Article  CAS  PubMed  Google Scholar 

  6. Khan I, Steeg PS. Metastasis suppressors: functional pathways. Lab Invest. 2018;98(2):198–210.

    Article  CAS  PubMed  Google Scholar 

  7. Harihar S, Ray S, Narayanan S, Santhoshkumar A, Ly T, Welch DR. Role of the tumor microenvironment in regulating the anti-metastatic effect of KISS1. Clin Exp Metastasis. 2020;37(2):209–23.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Tolg C, Messam BJ, McCarthy JB, Nelson AC, Turley EA. Hyaluronan functions in wound repair that are captured to fuel breast cancer progression. Biomolecules. 2021;11(11):1551–15517.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Choi S, Wang D, Chen X, Tang LH, Verma A, Chen Z, et al. Function and clinical relevance of RHAMM isoforms in pancreatic tumor progression. Mol Cancer. 2019;18(1):92.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Buttermore ST, Hoffman MS, Kumar A, Champeaux A, Nicosia SV, Kruk PA. Increased RHAMM expression relates to ovarian cancer progression. J Ovarian Res. 2017;10(1):66.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Soliman F, Ye L, Jiang W, Hargest R. Targeting hyaluronic acid and peritoneal dissemination in colorectal cancer. Clin Colorectal Cancer. 2022;21(2):e126–34.

    Article  PubMed  Google Scholar 

  12. Liu M, Tolg C, Turley E. Dissecting the dual nature of hyaluronan in the tumor microenvironment. Front Immunol. 2019;10:947.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Messam BJ, Tolg C, McCarthy JB, Nelson AC, Turley EA. RHAMM is a multifunctional protein that regulates cancer progression. Int J Mol Sci. 2021;22(19):10313–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Mateo F, He Z, Mei L, de Garibay GR, Herranz C, Garcia N, et al. Modification of BRCA1-associated breast cancer risk by HMMR overexpression. Nat Commun. 2022;13(1):1895.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Hinneh JA, Gillis JL, Moore NL, Butler LM, Centenera MM. The role of RHAMM in cancer: Exposing novel therapeutic vulnerabilities. Front Oncol. 2022;12: 982231.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Tolg C, Liu M, Cousteils K, Telmer P, Alam K, Ma J, et al. Cell-specific expression of the transcriptional regulator RHAMM provides a timing mechanism that controls appropriate wound re-epithelialization. J Biol Chem. 2020;295(16):5427–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Dunsch AK, Hammond D, Lloyd J, Schermelleh L, Gruneberg U, Barr FA. Dynein light chain 1 and a spindle-associated adaptor promote dynein asymmetry and spindle orientation. J Cell Biol. 2012;198(6):1039–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Maxwell CA, McCarthy J, Turley E. Cell-surface and mitotic-spindle RHAMM: moonlighting or dual oncogenic functions? J Cell Sci. 2008;121(Pt 7):925–32.

    Article  CAS  PubMed  Google Scholar 

  19. Choudhary M, Zhang X, Stojkovic P, Hyslop L, Anyfantis G, Herbert M, et al. Putative role of hyaluronan and its related genes, HAS2 and RHAMM, in human early preimplantation embryogenesis and embryonic stem cell characterization. Stem Cells. 2007;25(12):3045–57.

    Article  CAS  PubMed  Google Scholar 

  20. Jiang J, Mohan P, Maxwell CA. The cytoskeletal protein RHAMM and ERK1/2 activity maintain the pluripotency of murine embryonic stem cells. PLoS ONE. 2013;8(9):e73548.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Bahrami SB, Tolg C, Peart T, Symonette C, Veiseh M, Umoh JU, et al. Receptor for hyaluronan mediated motility (RHAMM/HMMR) is a novel target for promoting subcutaneous adipogenesis. Integr Biol (Camb). 2017;9(3):223–37.

    Article  CAS  PubMed  Google Scholar 

  22. Pibuel MA, Poodts D, Diaz M, Hajos SE, Lompardia SL. The scrambled story between hyaluronan and glioblastoma. J Biol Chem. 2021;296:100549.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Fluck MM, Schaffhausen BS. Lessons in signaling and tumorigenesis from polyomavirus middle T antigen. Microbiol Mol Biol Rev. 2009;73(3):542–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Valdes-Mora F, Salomon R, Gloss BS, Law AMK, Venhuizen J, Castillo L, et al. Single-cell transcriptomics reveals involution mimicry during the specification of the basal breast cancer subtype. Cell Rep. 2021;35(2):108945.

    Article  CAS  PubMed  Google Scholar 

  25. Christenson JL, Butterfield KT, Spoelstra NS, Norris JD, Josan JS, Pollock JA, et al. MMTV-PyMT and derived Met-1 mouse mammary tumor cells as models for studying the role of the androgen receptor in triple-negative breast cancer progression. Horm Cancer. 2017;8(2):69–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Lopez JI, Camenisch TD, Stevens MV, Sands BJ, McDonald J, Schroeder JA. CD44 attenuates metastatic invasion during breast cancer progression. Cancer Res. 2005;65(15):6755–63.

    Article  CAS  PubMed  Google Scholar 

  27. Tolg C, Hamilton SR, Nakrieko KA, Kooshesh F, Walton P, McCarthy JB, et al. Rhamm−/− fibroblasts are defective in CD44-mediated ERK1,2 motogenic signaling, leading to defective skin wound repair. J Cell Biol. 2006;175(6):1017–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Gibby K, You WK, Kadoya K, Helgadottir H, Young LJ, Ellies LG, et al. Early vascular deficits are correlated with delayed mammary tumorigenesis in the MMTV-PyMT transgenic mouse following genetic ablation of the NG2 proteoglycan. Breast Cancer Res. 2012;14(2):R67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Locke ME, Milojevic M, Eitutis ST, Patel N, Wishart AE, Daley M, et al. Genomic copy number variation in Mus musculus. BMC Genomics. 2015;16:497.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Yang H, Ding Y, Hutchins LN, Szatkiewicz J, Bell TA, Paigen BJ, et al. A customized and versatile high-density genotyping array for the mouse. Nat Methods. 2009;6(9):663–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Butler JL, Osborne Locke ME, Hill KA, Daley M. HD-CNV: hotspot detector for copy number variants. Bioinformatics. 2013;29(2):262–3.

    Article  CAS  PubMed  Google Scholar 

  32. Nik-Zainal S, Davies H, Staaf J, Ramakrishna M, Glodzik D, Zou X, et al. Landscape of somatic mutations in 560 breast cancer whole-genome sequences. Nature. 2016;534(7605):47–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Wang J, Li D, Shen W, Sun W, Gao R, Jiang P, et al. RHAMM inhibits cell migration via the AKT/GSK3beta/Snail axis in luminal A subtype breast cancer. Anat Rec (Hoboken). 2020;303(9):2344–56.

    Article  CAS  PubMed  Google Scholar 

  34. Owyong M, Chou J, van den Bijgaart RJ, Kong N, Efe G, Maynard C, et al. MMP9 modulates the metastatic cascade and immune landscape for breast cancer anti-metastatic therapy. Life Sci Alliance. 2019;2(6).

  35. Tolg C, Poon R, Fodde R, Turley EA, Alman BA. Genetic deletion of receptor for hyaluronan-mediated motility (Rhamm) attenuates the formation of aggressive fibromatosis (desmoid tumor). Oncogene. 2003;22(44):6873–82.

    Article  CAS  PubMed  Google Scholar 

  36. Nedvetzki S, Gonen E, Assayag N, Reich R, Williams RO, Thurmond RL, et al. RHAMM, a receptor for hyaluronan-mediated motility, compensates for CD44 in inflamed CD44-knockout mice: a different interpretation of redundancy. Proc Natl Acad Sci USA. 2004;101(52):18081–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Wang C, Thor AD, Moore DH 2nd, Zhao Y, Kerschmann R, Stern R, et al. The overexpression of RHAMM, a hyaluronan-binding protein that regulates ras signaling, correlates with overexpression of mitogen-activated protein kinase and is a significant parameter in breast cancer progression. Clin Cancer Res. 1998;4(3):567–76.

    CAS  PubMed  Google Scholar 

  38. Attalla S, Taifour T, Bui T, Muller W. Insights from transgenic mouse models of PyMT-induced breast cancer: recapitulating human breast cancer progression in vivo. Oncogene. 2021;40(3):475–91.

    Article  CAS  PubMed  Google Scholar 

  39. Elzamly S, Badri N, Padilla O, Dwivedi AK, Alvarado LA, Hamilton M, et al. Epithelial-mesenchymal transition markers in breast cancer and pathological responseafter neoadjuvant chemotherapy. Breast Cancer (Auckl). 2018;12:1178223418788074.

    PubMed  Google Scholar 

  40. Ross C, Szczepanek K, Lee M, Yang H, Qiu T, Sanford JD, et al. The genomic landscape of metastasis in treatment-naive breast cancer models. PLoS Genet. 2020;16(5):e1008743.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Banik S, Rakshit S, Sarkar K. The role of STAT1 in T helper cell differentiation during breast cancer progression. J Breast Cancer. 2021;24(3):253–65.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Gocher AM, Workman CJ, Vignali DAA. Interferon-gamma: teammate or opponent in the tumour microenvironment? Nat Rev Immunol. 2022;22(3):158–72.

    Article  CAS  PubMed  Google Scholar 

  43. Heimes AS, Madjar K, Edlund K, Battista MJ, Almstedt K, Gebhard S, et al. Prognostic significance of interferon regulating factor 4 (IRF4) in node-negative breast cancer. J Cancer Res Clin Oncol. 2017;143(7):1123–31.

    Article  CAS  PubMed  Google Scholar 

  44. Heimes AS, Hartner F, Almstedt K, Krajnak S, Lebrecht A, Battista MJ, et al. Prognostic significance of interferon-gamma and its signaling pathway in early breast cancer depends on the molecular subtypes. Int J Mol Sci. 2020;21(19):7178–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Cheon H, Wang Y, Wightman SM, Jackson MW, Stark GR. How cancer cells make and respond to interferon-I. Trends Cancer. 2023;9(1):83–92.

    Article  CAS  PubMed  Google Scholar 

  46. Weichselbaum RR, Ishwaran H, Yoon T, Nuyten DS, Baker SW, Khodarev N, et al. An interferon-related gene signature for DNA damage resistance is a predictive marker for chemotherapy and radiation for breast cancer. Proc Natl Acad Sci USA. 2008;105(47):18490–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Cheon H, Holvey-Bates EG, McGrail DJ, Stark GR. PD-L1 sustains chronic, cancer cell-intrinsic responses to type I interferon, enhancing resistance to DNA damage. Proc Natl Acad Sci USA. 2021;118(47).

  48. Rautela J, Baschuk N, Slaney CY, Jayatilleke KM, Xiao K, Bidwell BN, et al. Loss of host type-I IFN signaling accelerates metastasis and impairs NK-cell antitumor function in multiple models of breast cancer. Cancer Immunol Res. 2015;3(11):1207–17.

    Article  CAS  PubMed  Google Scholar 

  49. Du H, Xu T, Cui M. cGAS-STING signaling in cancer immunity and immunotherapy. Biomed Pharmacother. 2021;133:110972.

    Article  CAS  PubMed  Google Scholar 

  50. Mirzaei S, Saghari S, Bassiri F, Raesi R, Zarrabi A, Hushmandi K, et al. NF-kappaB as a regulator of cancer metastasis and therapy response: a focus on epithelial-mesenchymal transition. J Cell Physiol. 2022;237(7):2770–95.

    Article  CAS  PubMed  Google Scholar 

  51. Song X, Wei C, Li X. The signaling pathways associated with breast cancer bone metastasis. Front Oncol. 2022;12:855609.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Guerin MV, Regnier F, Feuillet V, Vimeux L, Weiss JM, Bismuth G, et al. TGFbeta blocks IFNalpha/beta release and tumor rejection in spontaneous mammary tumors. Nat Commun. 2019;10(1):4131.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Hunter KW, Amin R, Deasy S, Ha NH, Wakefield L. Genetic insights into the morass of metastatic heterogeneity. Nat Rev Cancer. 2018;18(4):211–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Nakano M, Taguchi R, Kikushige Y, Isobe T, Miyawaki K, Mizuno S, et al. RHAMM marks proliferative subpopulation of human colorectal cancer stem cells. Cancer Sci. 2023. https://doi.org/10.1111/cas.15795.

    Article  PubMed  Google Scholar 

  55. Pibuel MA, Poodts D, Molinari Y, Diaz M, Amoia S, Byrne A, et al. The importance of RHAMM in the normal brain and gliomas: physiological and pathological roles. Br J Cancer. 2023;128(1):12–20.

    Article  CAS  PubMed  Google Scholar 

  56. Shigeishi H, Biddle A, Gammon L, Emich H, Rodini CO, Gemenetzidis E, et al. Maintenance of stem cell self-renewal in head and neck cancers requires actions of GSK3beta influenced by CD44 and RHAMM. Stem Cells. 2013;31(10):2073–83.

    Article  CAS  PubMed  Google Scholar 

  57. Stepanenko AA, Vassetzky YS, Kavsan VM. Antagonistic functional duality of cancer genes. Gene. 2013;529(2):199–207.

    Article  CAS  PubMed  Google Scholar 

  58. Shiravand Y, Khodadadi F, Kashani SMA, Hosseini-Fard SR, Hosseini S, Sadeghirad H, et al. Immune checkpoint inhibitors in cancer therapy. Curr Oncol. 2022;29(5):3044–60.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Williams JB, Li S, Higgs EF, Cabanov A, Wang X, Huang H, et al. Tumor heterogeneity and clonal cooperation influence the immune selection of IFN-gamma-signaling mutant cancer cells. Nat Commun. 2020;11(1):602.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Mouse genotyping arrays were performed at the London Regional Genomics Centre RNA-seq. was performed at the Toronto Genomics Centre

Funding

Work was supported by funding from the Breast Cancer Society of Canada, London Health Sciences Catalyst grants (ET), Cancer Research Society (ET) and NSERC (NSERC-DG R3511A12. KH). ACN is supported by the American Cancer Society (132574-CSDG-18–139-01-CSM). J.B.M. is supported by the Atwater Fund, Elsa Pardee and Chairman's Fund Professor in Cancer Research.

Author information

Authors and Affiliations

Authors

Contributions

CT conceived and designed experiments, performed IHC, IF, Western blots and qRT-PCR. She interpreted results that were involved in writing as well as revising the manuscript. MM conceived and designed experiments, analyzed and interpreted genotyping array data and was involved in writing and revising the manuscript. FWQ conceived and designed experiments, analyzed and interpreted genotyping array data and was involved in writing and revising the manuscript. HAP conceived and designed experiments, analyzed and interpreted RNA-seq data and was involved in writing and revising the manuscript. JM designed and performed animal experiments, collected and interpreted data. MP designed, generated and characterized CRISPR/CAS9 cells. ACN, JBM, KAH, EAT obtained funding, conceived experiments, interpreted results, wrote and edited the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Kathleen A. Hill or Eva A. Turley.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Supplementary materials and methods.

Additional file 2.

Supplementary figures and tables.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tolg, C., Milojevic, M., Qi, F.W. et al. RHAMM regulates MMTV-PyMT-induced lung metastasis by connecting STING-dependent DNA damage sensing to interferon/STAT1 pro-apoptosis signaling. Breast Cancer Res 25, 74 (2023). https://doi.org/10.1186/s13058-023-01652-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13058-023-01652-1

Keywords

  • RHAMM
  • HMMR
  • Lung metastasis
  • STING
  • Interferon
  • apoptosis