Comparative Analysis of the Development of Acquired Radioresistance in Canine and Human Mammary Cancer Cell Lines
Mark Gray1*† 
Arran K. Turnbull2,3† 
James Meehan2† 
Carlos Martínez-Pérez2,3 
Charlene Kay2,3 
Lisa Y. Pang1 David J. Argyle1- 1The Royal (Dick) School of Veterinary Studies and Roslin Institute, University of Edinburgh, Edinburgh, United Kingdom
- 2Translational Oncology Research Group, Institute of Genetics and Molecular Medicine, Western General Hospital, University of Edinburgh, Edinburgh, United Kingdom
- 3Breast Cancer Now Edinburgh Research Team, Institute of Genetics and Molecular Medicine, Western General Hospital, University of Edinburgh, Edinburgh, United Kingdom
Research using in vitro canine mammary cancer cell lines and naturally-occurring canine mammary tumors are not only fundamental models used to advance the understanding of cancer in veterinary patients, but are also regarded as excellent translational models of human breast cancer. Human breast cancer is commonly treated with radiotherapy; however, tumor response depends on both innate radiosensitivity and on tumor repopulation by cells that develop radioresistance. Comparative canine and human studies investigating the mechanisms of radioresistance may lead to novel cancer treatments that benefit both species. In this study, we developed a canine mammary cancer (REM-134) radioresistant (RR) cell line and investigated the cellular mechanisms related to the development of acquired radioresistance. We performed a comparative analysis of this resistant model with our previously developed human breast cancer radioresistant cell lines (MCF-7 RR, ZR-751 RR, and MDA-MB-231 RR), characterizing inherent differences through genetic, molecular, and cell biology approaches. RR cells demonstrated enhanced invasion/migration capabilities, with phenotypic evidence suggestive of epithelial-to-mesenchymal transition. Similarities were identified between the REM-134 RR, MCF-7 RR, and ZR-751 RR cell lines in relation to the pattern of expression of both epithelial and mesenchymal genes, in addition to WNT, PI3K, and MAPK pathway activation. Following the development of radioresistance, transcriptomic data indicated that parental MCF-7 and ZR-751 cell lines changed from a luminal A classification to basal/HER2-overexpressing (MCF-7 RR) and normal-like/HER2-overexpressing (ZR-751 RR). These radioresistant subtypes were similar to the REM-134 and REM-134 RR cell lines, which were classified as HER2-overexpressing. To our knowledge, our study is the first to generate a canine mammary cancer RR cell line model and provide a comparative genetic and phenotypic analysis of the mechanisms of acquired radioresistance between canine and human cancer cell lines. We demonstrate that the cellular processes that occur with the development of acquired radioresistance are similar between the human and canine cell lines; our results therefore suggest that the canine model is appropriate to study both human and canine radioresistant mammary cancers, and that treatment strategies used in human medicine may also be applicable to veterinary patients.
Introduction
Naturally-occurring mammary tumors are the most frequently diagnosed cancer in bitches. These neoplasms represent ~50% of all canine tumors (1), of which 50% are malignant (2–4). Due to similarities in clinical features, relative age of onset, risk factors and tumor biology, canine mammary tumors (CMT) represent an excellent comparative and translational model for human breast cancer (HBC) (5–7). Gene expression profiles of primary (8) and metastatic (9) CMT have also shown similarities with HBC profiles, providing evidence that canine models can be utilized to help understand the genetic mechanisms of carcinogenesis in humans and dogs (4, 10, 11).
HBC is routinely managed with surgery followed by combinations of adjuvant chemotherapy, radiotherapy (RT), endocrine therapy or targeted therapy. RT is a commonly used breast cancer treatment; estimates indicate that curative or palliative RT can benefit ~83% of breast cancer patients (12). The use of adjuvant whole-breast RT after breast-conserving surgery has also been shown to deliver regional disease control and overall survival rates comparable with patients receiving a mastectomy, while also reducing side effects and improving cosmetic outcomes (13–15). Although 5-year survival rates for breast cancer patients following RT are ~80%, 30% will subsequently develop local recurrence and/or metastasis. Unfortunately, the vast majority of these patients have a poor prognosis and die within 5-years of disease progression (16). Cancer cells that possess intrinsic radioresistance, or develop acquired resistance during RT, can repopulate the tumor site after treatment. This can lead to treatment failures with the development of tumor recurrence and/or metastatic disease. Multiple factors are associated with the ability of cancer cells to develop acquired radioresistance, including signaling pathway dysregulation (e.g., EGFR/PI3K/AKT/mTOR), activation of DNA damage repair mechanisms, the existence of cancer stem cells, alterations in cancer metabolism and epithelial-to-mesenchymal transition (EMT). Hypoxic tumor microenvironments can also drive cancer cells to adopt an aggressive, treatment resistant phenotype (17). In order to develop treatment strategies to overcome/target the clinical issue of radioresistance, we require a detailed understanding of the mechanisms underlying acquired radioresistance.
In HBC, accurate disease staging is a mandatory requirement before beginning definitive treatment; classification systems that provide both predictive and prognostic information are commonly used to inform patient treatment (18). Histological grading with assessment of human epidermal growth factor 2 (HER2), estrogen receptor (ER), and progesterone receptor (PR) status provides an indication of the drivers of the disease and influences endocrine and/or targeted therapy use (19, 20). HBC can also be classified into intrinsic molecular subtypes that can predict treatment responses, overall survival and disease-free survival (21–24). These subtypes include: normal breast-like, luminal A, luminal B, HER2+, basal and “claudin-low” (21, 22, 24, 25). CMT can be solely epithelial (simple carcinoma or adenoma) or mesenchymal (fibrosarcoma, fibroadenoma, sarcoma, or osteosarcoma) in origin; however, a combination of mesenchymal and epithelial (carcinosarcoma or benign mixed tumors) or myoepithelial and epithelial tissues (complex adenoma or complex carcinoma) can also occur (26). Canine inflammatory mammary carcinoma is a rare CMT subtype that carries a very poor prognosis (27, 28). In contrast to human medicine, the diagnosis of CMT relies mainly on this histological grading system without undergoing specific receptor status evaluation or molecular subtype classification. The lack of integration of these diagnostic tests into the treatment decision-making process is largely because surgery alone is the main treatment option for dogs.
The effects of estrogen on CMT development have been previously established from a study that assessed the incidence of CMT in spayed and intact bitches. This study showed that only 0.5% of dogs spayed prior to their first season developed mammary tumors. However, levels rose to 8% and 26% when bitches were neutered following their first or second seasons, respectively. Additionally, no preventative effect was observed on the risk of CMT development when bitches were spayed following their second season (29). Although estrogen is associated with CMT development, the utilization of endocrine therapies, such as the selective ER modulator tamoxifen, has been reported in the literature only a limited number of times (30, 31). Conflicting results produced from these studies, combined with severe side effects associated with the use of anti-oestrogens in dogs (vulvar oedema, vaginal discharge, pyometra, and retinitis), have restricted their use in the treatment of CMT (31, 32). RT is infrequently employed in CMT treatment but could be used to improve regional disease control when cancer-free margins cannot be obtained at surgery, or as a palliative therapy for inflammatory mammary carcinomas or non-resectable tumors. Additional studies to elucidate the future role of RT in CMT treatment regimens are needed due to the impact it has in the management of HBC patients (33).
Compared to chemoresistance, there is an inadequate understanding of the mechanisms underlying radioresistance; this partly results from a lack of model systems for both human and veterinary applications. To begin to address this issue we have recently developed and characterized 3 novel HBC radioresistant (RR) cell lines (34). These included hormone-dependent ER+/PR+/HER2− (MCF-7 and ZR-751) and hormone-independent ER−/PR−/HER2− (MDA-MB-231) cell lines that represented different molecular HBC subtypes. In this new study we now go on to develop a novel CMT radioresistant cell line (REM-134), investigate the mechanisms of CMT radioresistance and perform a comparative analysis with our HBC RR cell line models. Although multiple CMT cell lines have previously been generated (35–38), the REM-134 cell line was the first to be developed, derived from a canine mammary carcinoma (39), a commonly occurring disease subtype with a poor prognosis. Even though it is commonly used in research, we chose to use the REM-134 cell line as it has yet to be fully characterized, especially in terms of its hormone receptor status. Its use in this study would therefore provide valuable information. Through genotypic, phenotypic, and functional analysis of our developed RR cell lines, we provide evidence that the cellular processes that occur with the development of acquired radioresistance are similar between human and canine cell lines; this suggests that not only is our canine model appropriate to study both canine and human radioresistant mammary cancers, but that treatment strategies used in human medicine may also be applicable to veterinary patients. To our knowledge, our study is the first to generate a CMT radioresistant model and provide a comparative analysis of the mechanisms of acquired radioresistance between human and canine breast cancer cell lines. We are also the first to report the use of canine multicellular tumor spheroids (MTS) originating from RR cells in immunohistochemical analysis and functional assays.
Materials and Methods
Cell Culture
Unless otherwise stated, reagents used for cell culture were acquired from Gibco Thermo Fisher Scientific (Paisley, UK). The HBC cell lines MCF-7, ZR-751 and MDA-MB-231 were cultured using Dulbecco's modified Eagle's medium (DMEM) with 10% fetal calf serum, 50 U ml−1 penicillin and 50 mg ml−1 streptomycin. The CMT cell line REM-134 was cultured in DMEM (high glucose) with the same additions. All cells were incubated at 37°C in a humidified atmosphere with 5% CO2. The HBC cell lines were obtained from the American Type Culture Collection (LGC Standards, Teddington, UK); the REM-134 cell line was a kind gift from Professor R.W. Else (College of Veterinary Medicine, University of Edinburgh, UK). All cell lines were authenticated at Health England (Porton Down, Salisbury, UK) by short tandem repeat (STR) profiling. All experiments were performed using cell lines maintained at low passage numbers with new cells acquired from frozen stocks after 10 passages.
Irradiation of Cells and Development of Radioresistant Cell Lines
A Faxitron cabinet X-ray system 43855D (Faxitron X-ray Corporation, IL, USA) was used to irradiate cells. Radioresistant cells were established from parental cell lines through exposure to weekly doses of radiation. Beginning with an initial dose of 2 Gy, cells were irradiated weekly with incremental doses of 0.5 Gy for 12-weeks. After this development period cells were maintained with additional doses of 5 Gy given every week. Cells were routinely passaged 24 h after each radiation dose.
Sulforhodamine B Proliferation (SRB) Assay
Cells were seeded into 96 well plates (500 cells/well). After 24 h of incubation, cells were irradiated and then fixed up to 120 h post-treatment by adding 50 μl 25% trichloracetic acid (Sigma-Aldrich, UK) at 4°C for 1 h. Wells were washed in dH2O and dried. Fifty microliter SRB dye [0.4% SRB dissolved in 1% glacial acetic acid (VWR International)] was added to the wells and the cells were incubated for 30 min. The cells were then washed 4 times in 1% glacial acetic acid and incubated for 60 min after the addition of 150 μl 10 mM Tris-NaOH buffer (pH 10.5). A Biohit BP800 spectrophotometer (Biohit Ltd, UK) and Wallac 1,420 Manager program (PerkinElmer, UK) were used to analyze optical density at 540 nm.
Colony Formation Assay
Cells were irradiated with varying doses 24 h after seeding into 75 mm plates (1,000 cells/plate). 1,9-dimethyl-methylene blue zinc chloride double salt (Sigma-Aldrich, UK) was used to fix and stain the cells between 10 and 14 days after seeding, when colonies consisting of at least 50 cells/colony had formed in the untreated control group. The plating efficiencies and survival fractions for treatment and control colony formation (CF) groups were calculated (40).
Scratch (Migratory) Assays
Cells were seeded into 6 well plates at densities that led to 100% confluence after 24 h. Scratch assays were carried out as previously described (41) following replacement of the media routinely used for cell culture with 0.1% serum-supplemented DMEM. Phase contrast images were taken (Axiovert DS100, x5 objective) up to 48 h post-scratch. The area lacking migrating cells, expressed as a % of the initial scratched area, was calculated at each time point using FIJI software.
Formation of Multicellular Tumor Spheroids
A single cell suspension of ~15,000,000 cells was placed into a spinner flask (Cellcontrol Spinner Flask, Integra, Switzerland) containing 100 ml of standard DMEM. The flask was placed onto a magnetic stirrer platform (Cellspin, Integra, Switzerland) with MTS forming over a period of 7 days under routine incubation conditions. Hypoxyprobe-1 (Hypoxyprobe, USA) was used to detect hypoxic areas within MTS. MTS were incubated with 100 μM hypoxyprobe-1 for 1 h before fixation.
3D Invasion Assay Using Multicellular Tumor Spheroids
MTS were transferred into the wells of a 24-well plate with 500 μl of collagen mix (0.22 M NaOH (Sigma-Aldrich, UK), fetal calf serum, 10x DMEM (Sigma-Aldrich, UK), cell matrix type 1-A (Alphalabs) and ice cold 0.1% filtered acetic acid at concentrations of 10, 10, 10, 25, and 45% respectively). The plates were incubated at 37°C for 1 h to allow collagen polymerization, after which 500 μl of routine DMEM was added. Phase-contrast images were taken at 24 h time-points up to 96 h (Axiovert DS100, x5 objective). Invasion was measured at each time point with a FIJI macro developed by Matthew Pearson (IGMM Advanced Imaging Resource, University of Edinburgh), with invasion expressed as a % of the initial MTS area.
Protein Isolation and Detection
Whole cell lysates were prepared as described previously (42). Equal protein amounts were separated by electrophoresis using sodium dodecyl sulfate (SDS) polyacrylamide gels. Proteins were transferred to Immobilon-P transfer membranes (Millipore, UK). Membranes were blocked using Odyssey Blocking Buffer (LI-COR Biosciences, UK) diluted 1:1 with PBS for 1 h. Membranes were then incubated overnight at 4°C with the appropriate primary antibodies (Table 1). IRDye 680LT (Li-Cor 926-68021, 1:10,000) and IRDye 800CW (Li-Cor, 926-32210, 1:10,000) secondary antibodies were added, and signals detected using a Li-Cor Odyssey Imager.
Table 1. Primary antibodies used for western blotting (WB), immunocytochemistry (ICC), and immunohistochemistry (IHC).
Immunohistochemistry
MTS were fixed in 4% formaldehyde (Genta Medical, UK) for 24 h. After fixation, the MTS were placed in 2% agarose and processed using the Thermo Scientific Excelsior AS Tissue Processor (Thermo Scientific, UK) and embedded in paraffin. Blocks containing MTS were cut with a Leica RM2235 rotary microtome (Leica Microsystems Ltd, UK); 4 μm sections were placed on SuperFrost Plus glass slides (Thermo Scientific, UK) and dried for 4 h at 53°C.
MTS were deparaffinized and rehydrated before undergoing antigen retrieval (Table 1). Incubation with 3% H2O2 solution (Dako, UK) was performed for 10 min to block endogenous peroxidase activity. Total Protein Block (Dako, UK) was used for 10 min to block non-specific antibody staining with the MTS subsequently incubated with primary antibodies for 1 h (Table 1). Envision labeled polymer (Dako, UK) was added for 30 min, after which DAB and substrate buffer (1:50) (Dako, UK) were added to visualize protein staining. Haematoxylin was used to counterstain the MTS, which were then dehydrated and mounted with coverslips using DXP mountant (Sigma-Aldrich, UK). A NanoZoomer ER slide scanner (Hamamatsu Photonics, UK) was used to scan all slides. The stained MTS were viewed using NanoZoomer Digital Pathology software. Image analysis was performed using QuPath version 0.1.2 (43).
RNA Extraction and Whole-Transcriptome Gene Expression Analysis
Cells were seeded in triplicate into 75 mm plates (3,000,000 cells/plate). After 24 h of incubation, pellets containing up to 10,000,000 cells were collected by trypsinisation, snap-frozen on dry ice and stored at −70°C for RNA extraction. The RNeasy Mini Kit using QIAshredder technology (UK Qiagen, Ltd) was used to extract RNA from the cells. Total RNA was purified from animal cells using spin technology, as per the manufacturer's protocol. The NanoDropTM Spectrophotometer ND1000 (Thermo Fischer Scientific) was used to quantify the RNA and assess for contaminants. RNA quality was evaluated by producing an RNA integrity number (RIN) for each sample (Agilent Bioanalyzer); each of the samples had RIN values above 9.7 (Supplementary Table 1). Lexogen QuantSeq 3' FWD sequencing technology was used to produce full genome expression read-counts on an Illumina flow cell; this was scanned with the Illumina HiScanSQ system (Edinburgh Clinical Research Facility's Genetic Core, University of Edinburgh). Next generation sequencing reads were generated toward the poly(A) tail with read 1 directly reflecting the mRNA sequence. The BlueBee high-performance next generation sequencing analysis software was used to pre-process the FASTQ files; this implements poly(A) tail trimming and alignment to the Genome Reference Consortium Human genome build 38 reference genome (for human cell lines) and the Canis lupus familiaris reference genome (for canine cell lines) using the Spliced Transcripts Alignment to a Reference (STAR) algorithm (44).
Data were filtered, removing all genes with <5 reads per sample in at least 90% of samples. Cell lines were mapped to species-specific Ensembl gene identifiers and cross-species matches were determined using Ensembl BioMART (45). In total 17,243 genes were mapped to human Ensembl gene identifiers and 13,703 were mapped to canine Ensembl gene identifiers. Following cross-species mapping, 9,692 genes were identified as common to both datasets. No significant difference in variance of expression was observed across all genes between both datasets (Supplementary Figure 1). Before the analysis took place, data were log2 transformed and quantile normalized in R (Bioconductor) software and packages (46). Heatmap and cluster analysis were implemented with the TM4 MeV (multiple experiment viewer) software (47). Heatmap clustering was performed using Pearson correlation with average linkage. Correction for batch effects was performed to integrate gene expression data produced in this study with public datasets; this was achieved using R with the ComBat package, as described previously (48, 49). Gene enrichment analysis was performed using the DAVID Functional Annotation Bioinformatics Microarray Analysis tool (50). Hierarchical clustering of both the parental and RR cell lines was accomplished using a published list of genes; the expression profile of the genes within this list denotes the breast cancer intrinsic subtypes (luminal A, luminal B, normal-like, basal, and HER2) (24). The genefu R package (51) was used to assign samples to the differing intrinsic subtypes. This package applies a Single Sample Predictor (SSP) algorithm which is a nearest-centroid classifier. Centroids signifying the molecular subtypes of breast cancer were identified through hierarchical clustering using the same intrinsic gene list that was used for cluster analysis within this study. All datasets generated and/or analyzed within this study are available in the NCBI's Gene Expression Omnibus (52) and are accessible through GEO Series accession number GSE149988.
Statistical Analysis
Two-way ANOVA with Holm-Šídák multiple comparisons test was employed in SRB, CF, and migration/invasion assays. Unpaired (two tailed) t-test was employed in the IHC analysis. P < 0.05 were deemed statistically significant. Data are shown as mean ± SEM. Statistical analysis was performed and graphs were generated with GraphPad Prism 8.
Results
Development and Confirmation of Acquired Radioresistance in Canine and Human Breast Cancer Cell Lines
Parental cell line intrinsic radiosensitivity and confirmation of acquired radioresistance in the developed cell lines were investigated through CF and SRB assays. Using the survival fraction of cells that were given a 2 Gy dose of radiation (SF2, a recognized experimental measure of radiosensitivity), a range of intrinsic radiosensitivities was found to be present in the parental cell lines, with the REM-134 cell line showing significantly greater radioresistance compared with the human cell lines. Within the parental human cell lines SF2 was not associated with molecular subtype. The ability of RR cells to form colonies was significantly higher than that of parental cells when subjected to a single radiation dose of up to 2 Gy, confirming the acquisition of radioresistance (Figure 1A). Significantly less inhibition of proliferation was also observed in RR cells in comparison to their parental cells when exposed to radiation doses of up to 10 Gy (Figure 1B). REM-134 RR and REM-134 RR cells that had not been exposed to radiation for 6 months (REM-134 rr) showed similar levels of radioresistance; this suggested that the changes involved in the acquisition of the RR phenotype are maintained over a long period of time (Figure 1C).
Figure 1. Confirmation of radioresistance using colony formation and SRB assays. (A) Colony formation assays comparing MCF-7, ZR-751, MDA-MB-231, and REM-134 cell lines with their derived RR cell lines at 10–14 days post-radiation treatment. (B) SRB assays comparing MCF-7, ZR-751, MDA-MB-231, and REM-134 cell lines with their derived RR cell lines at 120 h post-radiation treatment. (C) SRB assays comparing REM-134, REM-134 RR, and REM-134 rr cell lines at 120 h post-radiation treatment. The REM-134 rr is a radioresistant cell line that had not been radiated for 6 months (24 passages) before the experiment (2-way ANOVA with Holm-Šídák multiple comparisons test; data expressed as mean ± SEM, n = 3, ****p ≤ 0.0001; ***p ≤ 0.001; **p ≤ 0.01; *p ≤ 0.05).
Gene Expression Analysis Identifies Differences Between the Canine Parental and Radioresistant Cell Lines
Using an unsupervised analysis, a large number of genes were found to be inherently differentially expressed (using DESeq2—R, Bioconductor package) between the REM-134 and REM-134 RR cell lines (Figure 2). Higher expression of genes enriched for EMT, cell adhesion/motility, and response to hypoxia (cluster 1), with lower expression of genes involved in steroid biosynthesis, HIPPO signaling, focal adhesion, negative regulators of cell motility and proliferation (cluster 2) were evident in the REM-134 RR cell line compared to its parental cell line. Cluster analysis using these differentially expressed genes was performed using all human and canine parental and RR cell lines. The relative expression of key differentially expressed genes between the REM-134 and REM-134 RR cells associated with EMT (Cluster 1: BMP2, WNT5A, and SNAI1) and HIPPO signaling (Cluster2: WNT6, BMP4, FZD4, SNAI2) across all cell lines are shown in Supplementary Figure 2.
Figure 2. Global gene expression changes between the canine and human parental and RR cell lines. Heatmap showing log2 mean-centered gene expression profiles between parental and RR cell lines in respect of differentially expressed genes (false discovery rate = 0.01). Differences between the REM-134 and REM-134 RR cell lines are shown in the heatmap on the left, and the expression of the same differentially expressed genes across all human and canine parental and RR cells are shown in the heatmap on the right. Heatmap clustering was carried out using Pearson correlation with average linkage; red = higher expression, black = no change, green = lower expression. The gene list and the order in which they appear in the heatmap are shown in Supplementary Table 2.
Overall, results showed that the MDA-MB-231, MDA-MB-231 RR, MCF-7, and ZR-751 cell lines clustered together with similar expression patterns across all genes. The MCF-7 RR, ZR-751 RR, REM-134, and REM-134 RR cell lines formed a separate cluster. Interestingly, the MCF-7 RR and ZR-751 RR cells were found to have higher expression of cluster 1 genes compared to their parental lines, undergoing similar radiation-induced changes in expression to the REM-134 cells, with little change in expression of cluster 2 genes. These data suggest a common link between the expression of genes in cluster 1 and the development of radioresistance in the ER+ HBC and REM-134 cell lines. Cluster 2 expression appears to be specific to the REM-134 cell line; the expression of these genes does not appear to change in the MDA-MB-231 cells, suggesting that a different mechanism may underlie the development of radioresistance in this triple-negative model (Figure 2 and Supplementary Table 2). Analysis of these global gene changes provided the basis for further investigation of the key enriched pathways associated with the response to radiation.
Canine Radioresistant Cells Have Lower Expression of Cell Cycle-Associated Genes and Decreased Proliferation
2D and 3D cell models were used to assess the proliferative capabilities of parental and RR cells. MTS were the 3D model system used; IHC showed the presence of hypoxic areas both within and surrounding a central necrotic core, with a proliferating layer of cells around the periphery. These are recognized characteristics of MTS; therefore, their presence validated the use of the REM-134 MTS within this study (Figure 3). SRB assays were used to evaluate proliferation rates of cells in 2D cultures (Figure 4A). A range of proliferation rates were found to be present in the parental cell lines, with the REM-134 cell line showing a significantly greater proliferation rate compared with the human cell lines. Comparisons between the parental and RR cell lines identified lower rates of proliferation in MCF-7 RR, MDA-MB-231 RR, and REM-134 RR cells compared to their respective parental cell lines. The opposite was seen for the ZR-751 cell line, where the ZR-751 RR showed higher proliferation rates compared to the parental cell line. Gene expression data showed that both MCF-7 RR and ZR-75-1 RR cell lines exhibited lower expression levels of genes related to cell cycle regulation and G1/S-phase transition, including MCM4, MCM6, and cyclin D1 compared to their parental lines (Figure 4B and Supplementary Table 3). Expression levels of these genes were similar across the parental and RR MDA-MB-231 cell lines, both of which clustered together. In accordance with the SRB results, overall expression levels of these proliferation genes were higher in the REM-134 parental cells compared to all others. Lower expression levels of some genes were observed in the RR compared to the parental REM-134 cells, but overall expression remained higher than in the ER+ RR cell lines. Subsequent investigation of proliferation was carried out through IHC; MTS were stained for the proliferation marker Ki67 (Figure 4C) (both parental and RR MDA-MB-231 cells fail to generate MTS that withstood IHC processing). Quantitative IHC analysis identified that there was a lower percentage of Ki67-positive cells in the MCF-7 RR, ZR-751 RR, and REM-134 RR MTS compared with MTS formed from their parental cells. Overall, these results suggested that the RR cells have lower basal proliferation rates compared to their parental cells.
Figure 3. REM-134 MTS show hypoxic and proliferative gradients. The hypoxyprobe-1 compound was used to detect hypoxic regions within REM-134 MTS, while ki67 was used identify proliferating cells.
Figure 4. Radioresistant cell lines have modified basal proliferation rates relative to their parental cells. (A) SRB assays showing differences in proliferation rates between MCF-7, ZR-751, MDA-MB-231, and REM-134 cell lines and their derived RR cell lines grown in 2D cultures (2-way ANOVA with Holm-Šídák multiple comparisons test; data expressed as mean ± SEM, n = 3, ****p ≤ 0.0001; ***p ≤ 0.001). (B) Heatmap showing log2 mean-centered gene expression profiles of proliferation genes in parental and RR cell lines showing key G1/S phase regulators taken from the KEGG database cell cycle pathway (55); red = higher expression, black = no change, green = lower expression. Heatmap clustering was carried out using Pearson correlation with average linkage. The gene list is shown in Supplementary Table 3. (Ci) IHC of MTS stained for Ki67 using MCF-7, ZR-751, and REM-134 parental and RR cell lines. (Cii) Quantitative analysis of the % of cells with Ki67 staining (unpaired, two tailed t-test; data expressed as mean ± SEM, n = 3, ****p ≤ 0.0001).
Canine Radioresistant Cells Have Increased Invasion and Migration Potential
Morphological changes that occurred with the acquisition of radioresistance were identified through H&E staining of cells grown in 2D cultures. The human parental ER+ cell lines (MCF-7 and ZR-751) exhibited an epithelial-like morphology, comprising of tightly packed cells that form cobblestone-like monolayers, typical of luminal subtypes. However, their RR derivatives showed mesenchymal-type characteristics, with cells gaining a spindle-shaped morphology that contacted neighboring cells through focal points rather than the entire cellular circumference. The parental human ER− (MDA-MB-231) and the canine (REM-134) cell lines exhibited a mesenchymal-like phenotype, typical of basal and HER2-overexpressing subtypes; morphological changes in their RR derivatives were not as obvious as those observed in the RR cell lines derived from ER+ cells (Figure 5A).
Figure 5. REM-134 and REM-134 RR exhibit signs of EMT, with gene expression profiles similar to the human RR cell lines derived from ER+ cells lines. (A) H&E staining of cells grown in 2D cultures detailing the morphological differences between parental and RR cell lines. (B) IHC staining of EMT markers (vimentin, E-cadherin, and N-cadherin) in REM-134 and REM-134 RR MTS. (C) Heatmap showing log2 mean-centered gene expression profiles in respect of a published cancer cell EMT-signature (53); red = higher expression, black = no change, green = lower expression. The gene list and the order in which they appear in the heatmap are shown in Supplementary Table 4.
These observed changes in cell morphology provided evidence that the cells were undergoing EMT. We therefore assessed the protein expression levels of EMT markers through IHC. MTS formed from MCF-7 RR and ZR-751 RR cells demonstrated higher expression levels of N-cadherin, vimentin, and SNAIL, in addition to a partial downregulation of E-cadherin. The MDA-MB-231 cell line exhibited high levels of vimentin and SNAIL along with low E-cadherin and N-cadherin expression; no differences between the parental and RR MDA-MB-231 cell lines were identified (Supplementary Figure 3). Although no differences between the REM-134 and REM-134 RR were identified, the pattern of protein expression was similar to that seen in the MCF-7 RR and ZR-751 RR cell lines, with all cell lines showing expression of E-cadherin and vimentin (Figure 5B and Supplementary Figure 3). Gene expression analysis looking at the expression patterns for genes in a published EMT signature produced similar results (53). This study produced a pan-cancer EMT-associated gene expression signature by merging bioinformatic expression data from both The Cancer Cell Line Encyclopedia and The Cancer Genome Atlas (the lists of genes used in our study are provided in Supplementary Table 4). In our analysis, both MCF-7 and ZR-751 parental cell lines exhibited expression patterns consistent with an epithelial genotype, while the parental and RR MDA-MB-231 cells displayed higher expression of mesenchymal genes. The MCF-7 RR, ZR-751 RR, REM-134, and REM-134 RR cell lines all demonstrated a diverse expression pattern, with comparatively higher expression levels of both mesenchymal and epithelial genes, suggestive of a hybrid or transitional phenotype (Figure 5C). The expression levels of key EMT genes (BMP2, WNT5A, and SNAI1) are shown across all cell lines in Supplementary Figure 2. Expression levels of these key genes were found to be increased in the ER+ HBC RR and REM-134 RR cell lines compared to their parental lines. A corresponding decrease in expression was observed in key members of the HIPPO tumor suppressor pathway (WNT6, BMP4, FZD4, and SNAI2) in ER+ HBC RR and REM-134 RR cell lines compared to their parental lines.
Following the identification of cellular changes suggestive of EMT, we investigated the invasive and migratory characteristics of the cell lines. 2D migration assay results demonstrated that the RR cell lines all had significantly increased migratory ability in comparison to their parental cells. Similarly, using 3D invasion assays, the MCF-7 RR, ZR-751 RR, and REM-134 RR cells had increased invasive potential compared to their parental cells (Figures 6A,B).
Figure 6. Radioresistant cell lines have increased migration and invasion potential. (A) Images of 2D migration and 3D MTS invasion assays comparing the parental and the derived RR cell lines. (B) Graphs exhibiting the migration (Bi) and invasion assay (Bii) results. For the migration assays the relative migratory distance was calculated at each time point up to 48 h and expressed as a % area devoid of cells based on the initial scratched area at day 0. Invasion was assessed up to 96 h post-seeding. Area of MTS at each time point was calculated and expressed as a % of initial MTS area at day 0 (2-way ANOVA with Holm-Šídák multiple comparisons test; data expressed as mean ± SEM, n = 3, ****p ≤ 0.0001; ***p ≤ 0.001; *p ≤ 0.05.
Canine Radioresistant Cell Lines Exhibit Enhanced WNT Signaling
WNT signaling was studied due to its apparent role in radioresistance and EMT (54); this pathway was investigated using WNT signaling pathway gene expression signatures and WNT signaling downstream targets, both acquired from the KEGG database (55) (genes lists are provided in Supplementary Table 5). The ZR-751 and MCF-7 parental cells were found to have lower expression of WNT pathway and target genes compared to the other cell lines. The MCF-7 RR, ZR-751 RR, REM-134, and REM-134 RR cell lines clustered together with a similar pattern of gene expression consistent with WNT target gene activation. Interestingly, the pattern of expression of WNT signaling pathway members was found to be different between the ER+ HCB RR cell lines and both the REM-134 parental and RR lines, potentially suggesting that different WNT signaling pathways in these models may be responsible for the downstream WNT target activation. When the REM-134 and REM-134 RR cell lines were analyzed separately for key differentially expressed (false discovery rate = 0.01) members of the WNT pathway represented in the dataset, clear differences were observed, with the REM-134 RR cell line showing overall higher expression of frizzled family members 1/3/6 and lower expression of frizzled members 2/4, in addition to around 4-fold higher expression of WNT5a (Figure 7A). These results were further investigated using IHC in REM-134 and REM-134 RR MTS. Results showed significantly increased expression of WNT5a and pan-frizzled in the REM-134 RR MTS compared to the parental MTS (Figure 7B). Increased WNT5a expression was also identified in MCF-7 RR and ZR-751 RR MTS in comparison to their respective parental cell lines (Supplementary Figure 4).
Figure 7. WNT signaling is increased in MCF-7 RR, ZR-751 RR, and REM-134 RR cell lines. (A) Heatmap showing log2 mean-centered gene expression profiles between parental and RR cell lines in respect of the WNT signaling pathway (left heatmap) and WNT target genes (right heatmap). WNT5a and represented FRIZZLED genes are shown separately for REM-134 and REM-134 RR cell lines (center heatmap). Genes taken from the KEGG pathway database (55); red = higher expression, black = no change, green = lower expression. The gene lists are shown in Supplementary Table 5. Heatmap clustering was carried out using Pearson correlation with average linkage. (B) IHC of pan-Frizzled and WNT5a expression in REM-134 and REM-134 RR MTS with quantitative analysis of the % of positively stained cells (unpaired, two tailed t-test; data expressed as mean ± SEM, n = 3, ****p ≤ 0.0001).
Canine Radioresistant Cell Lines Maintain Their Original Intrinsic Breast Cancer Subtype
Receptor status in the REM-134 and REM-134 RR cell lines was investigated through IHC using MTS. No change in receptor expression was identified with the acquisition of radioresistance, with both the REM-134 and REM-134 RR cell lines classified as ER−/PR−/HER2+ (Figure 8A and Table 2). Although no change was identified in receptor expression between the parental and RR MDA-MB-231 cell lines (both classified as ER−/PR−/HER2−), a difference was seen in the MCF-7 and ZR-751 cell lines. Both MCF-7 and ZR-751 cell lines were classified as ER+/PR+/HER2−; however, their RR derivatives lost ER and PR expression, becoming ER−/PR−/HER2− (Supplementary Figure 5). Further investigation of cell line classification was performed through integration of the gene expression data from this study with a public gene expression dataset (GSE50811) of 51 breast cancer cell lines. Both the REM-134 and REM-134 RR cell lines clustered tightly and were classified, by correlation to centroids, as belonging to the HER2-overexpressing subtype. They also clustered near the MCF-7 RR and ZR-751 RR cell lines, which were classified as basal/HER2-overexpressing and as normal-like/HER2-overexpressing, respectively. The classification of the MCF-7 RR and ZR-751 RR was different to their parental lines; as anticipated, both of the parental cell lines were classified as luminal A. Predictably, the parental and RR MDA-MB-231 cell lines clustered near each other and were classified as the basal breast cancer subtype (Figure 8B and Table 2).
Figure 8. Breast cancer receptor status, ER signaling and intrinsic breast cancer subtype in the canine and human parental and RR cell lines. (Ai) IHC of ERα, HER2 and PR expression in REM-134 and REM-134 RR MTS. (Aii) Quantitative analysis of the % of ERα, PR and HER2 positively stained cells (unpaired, two tailed t-test; data expressed as mean ± SEM, n = 3). (B) The data generated within this study was integrated with a public gene expression dataset (GSE50811) of 51 breast cancer cell lines. Hierarchical clustering of parental and RR cell lines was based on the Sørlie 2003 intrinsic genes (24); subgroup classification was based on correlation to centroids and was performed using the genefu packaging in R; red = basal, dark blue = luminal A, light blue = luminal B, purple = HER2-overexpressing, green = normal-like. (C) Heatmap showing the log2 mean-centered expression profile of a published ER signaling gene signature (56); red = higher expression, black = no change, green = lower expression. Heatmap clustering was carried out using Pearson correlation with average linkage.
Table 2. Immunohistochemical and molecular subtype classification of the canine and human parental and RR cell lines.
ER signaling was investigated using a published ER signaling gene expression signature (56). As expected, the human ER+ parental cell lines (MCF-7 and ZR-751) were characterized by high expression of all of these genes, while in comparison their RR derivatives, together with the REM-134 and REM-134 RR cell lines, were found overall to have lower expression levels of these genes (Figure 8C).
Canine Parental and Radioresistant Cell Lines Show PI3K and MAPK Activity
After we identified that the canine cell lines expressed HER2, signal transduction pathways downstream of the HER/ERBB tyrosine-kinase receptor family were further evaluated. MAPK pathway activity was assessed using a gene expression signature that had previously been published (57). Using a combination of FOXO-regulated genes (which, as downstream targets of inhibition by PI3K, have the opposite pattern of expression to PI3K activity) (58) and genes obtained from the KEGG pathway database (55), PI3K activity was assessed. Results from supervised gene expression analysis were consistent with PI3K and MAPK signaling activity in the REM-134, REM-134 RR, MCF-7 RR and ZR-751 RR cell lines, whereas the MCF-7 and ZR-751 cell lines exhibited inactive PI3K and MAPK activity. Inactive PI3K and active MAPK signaling were also observed in the parental and RR MDA-MB-231 cell lines (Figures 9A,B). Western blot analysis of untreated lysates from all cell lines confirmed the presence of phosphorylated ERK1/2 (MAPK activation) and phosphorylated AKT (PI3K activation) (Supplementary Figure 6).
Figure 9. MAPK and PI3K pathway activity in the canine and human parental and RR cell lines. (A) Heatmap showing log2 mean-centered gene expression profiles between parental and RR cell lines in respect to a MAPK pathway activity gene signature (57). Heatmap clustering was carried out using Pearson correlation with average linkage, red = higher expression, black = no change, green = lower expression. (B) Heatmap showing log2 mean-centered gene expression profiles between parental and RR cell lines in respect of the PI3K pathway [associated genes taken from the PI3K KEGG pathway (upper heatmap) and FOXO-regulated genes (lower heatmap) (55)]. Heatmap clustering was carried out using Pearson correlation with average linkage, red = higher expression, black = no change, green = lower expression.
Discussion
Radiotherapy is a frequently used curative and palliative treatment for a wide range of human and canine tumors. Unfortunately, intrinsic and acquired radioresistance can significantly limit its efficacy and ultimately leads to local recurrence, disease progression or metastasis. In this study, we developed a canine mammary cancer radioresistant cell line and investigated the cellular mechanisms related to the development of acquired radioresistance. We subsequently performed a comparative analysis of this resistant model with our previously developed HBC radioresistant cell lines, characterizing their inherent differences through genetic, molecular and cell biology approaches.
Intrinsic radiosensitivities of the panel of cell lines was first investigated through CF assays. The REM-134 cell line showed significantly greater radioresistance than that of the human cell lines. Intrinsic radioresistance of the 3 human cell lines was not related to subtype, with the ER+ and ER− cell lines showing a similar response to doses of radiation up to 2 Gy. Clinical research studies have identified that HBC subtype is associated with tumor radiosensitivity. One such study investigated invasive breast cancer local recurrence rates following breast-conserving surgery with subsequent adjuvant RT. Their results identified recurrence rates of 0.8% for luminal A, 1.5% for luminal B, 7.1% for basal and 8.4% for HER2-overexpressing HBC subtypes (59). Other studies indicate that triple negative breast cancers and HER2-overexpressing breast cancers that are treated with post-mastectomy RT have increased risks of locoregional recurrence and metastasis, along with significantly reduced overall survival (59–61). However, our results are in accordance with an in vitro cell line study which showed that there was no association between HBC subtype and their intrinsic radiosensitivity (62).
Following the 12-week radiation exposure protocol, radioresistance development was verified using CF and SRB proliferation assays. All developed RR cell lines showed greater resistance to a single fractionated radiation dose compared to their parental cells. These results validated their utility as an in vitro model system to characterize their resistant phenotype and examine the mechanisms associated with the development of acquired radioresistance. Our protocol was able to generate RR models that exhibited significant differences in CF ability compared to parental cells when exposed to 2 Gy, a standard treatment dose routinely used in HBC patients. This is an important consideration for any resistance development protocol, as differences seen between radiosensitive and radioresistant cells at doses comparable to those used in the clinic will produce more translational data. Our development protocol therefore has significant advantages over others used in the literature, which only managed to generate HBC RR cell lines that had significantly different CF ability compared to their respective parental cells at considerably higher radiation doses (8–10 Gy) (63). Our protocol was also able to produce a radioresistant phenotype that was maintained in REM-134 RR cells that had not received radiation for 6 months. This result is similar to that observed with our previously developed MCF-7 RR cell line (34); this indicates that the acquisition of radioresistance in both our human and canine RR cell line models was not transient. Our results differ from those of other studies which were unable to generate stable HBC radioresistant models (64). These differing results are likely due to the use of different radioresistance development protocols in each study and highlights the need to perform frequent CF and/or SRB proliferation assays to verify maintenance of the radioresistant phenotype.
Following the generation of the REM-134 RR cell line, we evaluated the ability of both the parental and RR cells to generate MTS for use within the study. The HBC MTS used in this study have been previously validated for use as 3D tumor models (34). Using the spinner flask method, we successfully generated MTS from both canine cell lines. This result, to our knowledge, is the first time that REM-134 MTS have been generated. MTS reproduce various aspects of the in vivo tumor microenvironment, including low oxygen levels, creating necrotic, peri-necrotic, and hypoxic regions. Proliferative gradients found within cancers are also present in MTS (65). These characteristics were present within the REM-134 MTS produced in this study. Hypoxyprobe-1, a chemical that has been used to detect hypoxic regions in various cancer models (66–68), was used in our MTS model. Hypoxyprobe-1 staining showed that low oxygen regions were present predominantly in the more central areas of MTS, both within and surrounding the necrotic core. Proliferating cells, identified through Ki67 staining, were largely located around the periphery.
Transcriptomic data was initially used to distinguish genes whose expression was significantly changed between the REM-134 and REM-134 RR cell lines following the generation of the RR models. Enriched pathways identified from this analysis, which had previously been linked with radioresistance in studies using human samples, were taken forward for further investigation; these included proliferation, EMT and WNT signaling pathways.
Differences in proliferation between the cell lines were investigated through SRB assays, transcriptomic and IHC analysis. Although the SRB data showed that the ZR-751 RR cell line had an increased rate of proliferation compared to its parental cell line, all of the other analyses, including the gene expression profiles and IHC for the other cell lines, suggested that both the human and canine RR cell lines had reduced proliferation rates. These results are in line with a previous study which identified reduced proliferation rates in human prostate cancer RR cell lines (69). Clinical effects of RT are largely due to the direct and indirect DNA damage it causes. In rapidly-dividing cells there is little time for DNA damage to be repaired by processes such as non-homologous end joining and homologous recombination. If a damaged cell enters cell division with uncorrected DNA damage, then the cell will likely die due to various radiation-induced cell death mechanisms such as mitotic catastrophe, apoptosis, autophagy, necrosis, or senescence. Lower proliferation rates may therefore provide these RR cells with greater time to repair these sites of DNA damage.
Among CMT patients, ~50% of carcinomas will metastasize to local lymph nodes, leading to distant metastases (lung and bone) and death (70, 71). A similar situation occurs in HBC patients; studies indicate that ~7% of patients present with metastatic tumors, while ~20% of patients that are initially diagnosed with local disease subsequently go on to develop metastatic disease (72). Studies that investigate the molecular mechanisms involved in cancer metastasis are therefore of significant clinical importance. In our study, all our RR cell lines showed increased migration and invasion abilities, suggestive of a more aggressive phenotype; this phenotypic change could indicate that RR cells have greater locally invasive and metastatic potentials, factors which are poor prognostic indicators for clinical patients.
Using an initial unsupervised analysis, downregulation of genes involved in HIPPO signaling was identified in the REM-134 RR cell line. This pathway regulates cellular proliferation and apoptosis and has multiple family members that act as tumor suppressors (73). Activation of the HIPPO signaling pathway has also been shown to antagonize WNT signaling, whereas inhibition of HIPPO signaling (causing YAP hypophosphorylation and nuclear localization) can induce EMT (74, 75). Downregulation of HIPPO signaling may therefore be involved in the increase in WNT signaling and EMT seen in the REM-134 RR cell line.
Malignant cellular transformation that leads to loss of epithelial morphology, reduced cellular contact and increased cell migration/invasion is an important feature of EMT (76–78) and is associated with poor prognosis (79). Significant morphological changes were observed in the MCF-7 RR and ZR-751 RR cell lines compared with their parental cells, with the RR cells gaining a more mesenchymal phenotype. The morphology of the MDA-MB-231 and REM-134 cell lines was quite different to that of the parental ER+ cell lines, with the former exhibiting a typical mesenchymal phenotype. Although no significant differences were seen in their RR derivatives, their morphology was similar to the MCF-7 RR and ZR-751 RR cell lines. IHC and gene expression analysis were used to investigate these results further. Firstly, using IHC to investigate established EMT breast cancer biomarkers (80), we demonstrated that the MCF-7 RR and ZR-751 RR cell lines exhibited downregulated E-cadherin expression and upregulated vimentin, N-cadherin and SNAIL expression compared to their parental cells. In contrast to these results, the parental and RR MDA-MB-231 cell lines exhibited low E-cadherin and high vimentin expression. Interestingly, the REM-134 and REM-134 RR cell lines showed a hybrid/intermediate epithelial-mesenchymal phenotype, with E-cadherin and vimentin both being expressed. These results were supported by identifying gene expression profiles in our cell lines through the use of a previously published cancer cell EMT-signature (53). Similarities were seen in the MCF-7 RR, ZR-751 RR, REM-134, and REM-134 RR cell lines, which all showed a mixed expression pattern with high expression of mesenchymal and epithelial genes. These results again suggested that these RR cell lines possess a hybrid/transitional phenotype that traverses the epithelial and mesenchymal states (77). Previous studies have suggested that a hybrid EMT state is linked with the presence of stem-like properties, increased cellular plasticity and enhanced migratory/metastatic abilities. It is thought that this state may represent a cellular survival response to stressful environments (81, 82). A multitude of signaling cascades can stimulate EMT in non-cancerous and cancerous cells, including various receptor tyrosine kinase pathways, Wnt-β-catenin and Notch signaling (83, 84). RT can also activate EMT through increasing the expression of TGFβ (85, 86). In our cell lines, WNT signaling was activated in the MCF-7 RR, ZR-751 RR, REM-134, and REM-134 RR cell lines, suggestive of a potential mechanism through which the cells underwent EMT. These results also showed that the acquisition of radioresistance in the human ER+ cell lines led to the development of a phenotype similar to that of the canine cell lines.
Previous research studies have investigated the value of CMT as a metastatic model for HBC. Metastatic CMT, in comparison to non-metastasizing tumors, have been shown to exhibit upregulation of genes associated with cell cycle regulation, DNA damage repair, extracellular matrix remodeling, proteasomal degradation and protein folding, while genes involved in cellular differentiation, growth factor signaling and actin organization are downregulated. Of these differentially expressed canine genes, 25% were discovered to be linked to HBC (9). Comparable results to these reported at gene level have also been observed at the intracellular protein level. One study detected 21 proteins (predominantly associated with cell adhesion, extracellular matrix remodeling and hypoxic resistance) that were differentially expressed in canine mammary carcinomas which were classified as either metastasizing or non-metastasizing. The majority (19/21) of these proteins were linked with metastasis or malignancy in a range of human cancers, of which 9 had comparable expression patterns to that seen in HBC patients (87). The partly overlapping transcriptome and proteome of metastatic CMT and HBC indicates that there must be similar pathways/mechanisms involved in breast carcinogenesis and pathogenesis between the two species. These studies also demonstrate that metastatic CMT are an appropriate translational model for metastatic HBC. Similarities between the REM-134 and REM-134 RR transcriptomic data produced in this study and the results from the previous studies discussed here provides further evidence of the value of using these cell line models as a metastatic model of human and canine disease.
HBC is typically graded and characterized through IHC with analysis of expression levels of various receptors such as HER2, ER and PR. Additionally, gene expression profiling has been successfully employed to classify breast tumors into luminal A, luminal B, HER2-overexpressing, basal and normal-like intrinsic subtypes (21, 22). As previously mentioned, these varying subtypes exhibit differing inherent sensitivities to RT, indicate prognosis and can influence which patients receive endocrine and/or targeted therapies (61, 88). To characterize the REM-134 and REM-134 RR cell lines within the context of HBC, we investigated the expression of HER2, ERα, and PR and performed molecular profiling. IHC showed that both parental and RR REM-134 cells expressed HER2, with these cell lines also classifying as HER2-overexpressing through their transcriptional profiles; receptor expression and subtype classification did not change with the acquisition of radioresistance. Similarities again were seen in the RR models produced from ER+ cell lines, with a change from luminal A for the MCF-7 and ZR-751 cells to a non-luminal classification for their RR derivatives. The MCF-7 RR cells correlated with the basal/HER2-overexpressing subtypes, while the ZR-751 RR cells correlated with the normal-like/HER2-overexpressing subtypes. Prior to the introduction of HER2-targeted therapies, HER2-overexpressing breast cancers carried a high locoregional recurrence risk and poor overall prognosis (22–24, 89). Luminal A breast cancers typically have an excellent response to RT and endocrine treatments (89); a move away from this subtype classification, as seen in our cell lines, would be consistent with an aggressive, treatment-resistant phenotype. These results indicate that acquired radioresistance can be associated with cellular plasticity, and that gene expression changes can lead to an alteration in molecular subtype. Our results also demonstrate that the RR cells derived from ER+ cell lines again show significant similarities to the canine cell lines.
HER2 expression occurs in ~15–30% of all HBC patients (90–92). HER2 protein overexpression has been determined to be both predictive for tumor response to HER2-targeted treatments and prognostic for disease outcome in human patients (93). Several studies have suggested that ~35% of malignant CMT have either HER2 gene or protein expression (94–97) and that HER2 expression is associated with histological grade, proliferation index and tumor size (94, 96, 98). Human and canine genome sequencing has identified significant homology between the HER2 antigens in both species and although further research is needed, results such as these suggest that human-based immunotherapies (e.g., pertuzumab or trastuzumab) or human developed tyrosine kinase inhibitors, could be successfully employed in HER2+ CMT (99). The REM-134 cell line and our developed REM-134 RR model would therefore be good in vitro models to investigate HER2 signaling in CMT.
In support of the REM-134 HER2-overexpressing classification, we showed that REM-134 cells exhibited lower expression of ER-driven genes in comparison to the human ER+ MCF-7 and ZR-751 cell lines. Again, the REM-134 and REM-134 RR cell lines clustered closely with the MCF-7 RR and ZR-751 RR cell lines. In HBC, HER2 expression and ER activity have been shown to have an inverse relationship, with HER2 overexpression being associated with reduced sensitivity to endocrine therapies (100–102). Because of the change in subtype classification and the suggestion of a shift away from ER signaling in the MCF-7 RR and ZR-751 RR cell lines, a situation akin to that seen in the canine cell lines, downstream signaling pathways of the HER/ERBB family were investigated. The phosphatidylinositol-3-kinase (PI3K)/protein kinase B (AKT) pathway is commonly hyperactivated in various cancer types and leads to cellular responses related to survival, proliferation and metabolism (103, 104). RT can also activate the PI3K/AKT pathway, which is associated with intrinsic radioresistance, proliferation, and resistance to hypoxic environments (105, 106). Similarly, activation of the mitogen-activated protein kinase (MAPK) pathway is thought to be a cytoprotective response which can allow cancer cells to repopulate the tumor during fractionated RT (107–111). In this study we found that MCF-7 RR, ZR-751 RR, REM-134, and REM-134 RR gene expression signatures corresponded with activation of the PI3K and MAPK pathways; these pathways may therefore play a significant role in the development of radioresistance in both human and canine patients.
The development of novel RR CMT cell lines opens up the possibility of future in vivo studies. Xenograft studies using orthotopic (mammary fat pad) or subcutaneous implantation of paired sensitive and RR REM-134 cells could be used for a multitude of pre-clinical research opportunities. Studies evaluating the effectiveness of radiosensitizing agents, the development of metastatic models or the detection of serum-based biomarkers of radiation response are all achievable with the generation of RR cell lines such as ours.
Conclusion
This study is the first to report the development and characterization of a novel canine mammary cancer RR cell line which was used as a comparative model for HBC. The generation of new radioresistant models is important, as these will aid the understanding of the molecular mechanisms that drive the development of radioresistance. Similarities in terms of EMT, WNT signaling, estrogen regulation, HER signaling, and subtype classification were identified between the RR cell lines derived from the human ER+ cell lines and the canine parental and RR cell lines. These results suggest that the mechanisms involved in the acquisition of radioresistance may be similar in the 2 species. As we continue to appreciate the significant similarities between human and canine mammary tumors, comparative studies will become more important for the investigation of carcinogenesis in both species. We believe that comparative studies of resistant disease will be fundamental for future research, leading to the development of novel treatment strategies that are equally applicable to both human and veterinary patients.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/geo/, GSE149988.
Author Contributions
DA secured funding for this research. MG conceptualized, performed, analyzed, interpreted the laboratory work, wrote the majority of the manuscript, and composed the figures, with significant contributions from AT and JM. AT performed all bioinformatic analysis. CM-P uploaded data to GEO. Critical revisions were made by MG, AT, JM, CM-P, CK, LP, and DA. All authors read and approved the final manuscript.
Funding
This work was supported by the Breast Cancer Institute Fund from the Edinburgh & Lothians Health Foundation.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors gratefully acknowledge the Edinburgh Clinical Research Facility's Genetic Core for their expertise and assistance in this body of work.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2020.00439/full#supplementary-material
Supplementary Figure 1. (A) Venn diagram showing overlap in gene mappings between human and canine Ensembl gene IDs. (B) Violin plot comparing the variance for all matched genes between human and canine datasets. Statistical comparison was performed using two-tailed Mann-Whitney test (non-significant).
Supplementary Figure 2. Log2 mean-centered gene expression profiles of key EMT and Hippo associated genes across all cell lines (selected genes from heatmap in Figure 2). Genes were selected from differential gene expression analysis comparing global gene expression between REM-134 parental and RR cell lines. Heatmap clustering was carried out using Pearson correlation with average linkage based on all differentially expressed genes; red = higher expression, black = no change, green = lower expression.
Supplementary Figure 3. IHC and ICC staining of EMT markers (vimentin, E-cadherin, N-cadherin, and SNAIL) in MCF-7, ZR-751, and MDA-MB-231 parental and RR cell lines.
Supplementary Figure 4. IHC staining for WNT5a in MCF-7 and ZR-751 parental and RR MTS.
Supplementary Figure 5. IHC and ICC staining of hormone receptors (ER, PR, and HER2) in MCF-7, ZR-751, and MDA-MB-231 parental and RR cell lines.
Supplementary Figure 6. (A) Combined western blot showing total and phosphorylated ERK1 and ERK2 and total and phosphorylated pan-AKT in parental and RR cell lines. Splicing was performed in the images shown in part (b) to remove redundant lanes to produce this figure. (B) Original western blot images from radiation time course experiments in all 8 cell lines. Bands shown in part (a) are highlighted here, which represent untreated control samples.
Supplementary Table 1. RIN values for all samples used in gene expression analysis.
Supplementary Table 2. List of genes differentially expressed between the REM-134 and REM-134 RR cell lines.
Supplementary Table 3. List of G1/S phase regulatory genes taken from the KEGG database cell cycle pathway.
Supplementary Table 4. Epithelial and mesenchymal gene lists from a published EMT signature.
Supplementary Table 5. WNT pathway signaling and target genes taken from the KEGG database.
References
1. Moe L. Population-based incidence of mammary tumours in some dog breeds. J Reprod fertil Suppl. (2001) 57:439–43.
2. Benjamin S, Lee A, Saunders W. Classification and behavior of canine mammary epithelial neoplasms based on life-span observations in beagles. Vet Pathol. (1999) 36:423–36. doi: 10.1354/vp.36-5-423
3. Fidler I, Brodey R. The biological behavior of canine mammary neoplasms. J Am Vet Med Assoc. (1967) 151:1311–8.
4. Nguyen F, Peña L, Ibisch C, Loussouarn D, Gama A, Rieder N, et al. Canine invasive mammary carcinomas as models of human breast cancer. Part 1: natural history and prognostic factors. Breast Cancer Res Treat. (2018) 167:635–48. doi: 10.1007/s10549-017-4548-2
5. Pinho SS, Carvalho S, Cabral J, Reis CA, Gärtner F. Canine tumors: a spontaneous animal model of human carcinogenesis. Transl Res. (2012) 159:165–72. doi: 10.1016/j.trsl.2011.11.005
6. Rivera P, Von Euler H. Molecular biological aspects on canine and human mammary tumors. Vet Pathol. (2011) 48:132–46. doi: 10.1177/0300985810387939
7. Gray M, Meehan J, Martínez-Pérez C, Kay C, Turnbull AK, Morrison LR, et al. Naturally-occurring canine mammary tumors as a translational model for human breast cancer. Front Oncol. (2020) 10:617. doi: 10.3389/fonc.2020.00617
8. Uva P, Aurisicchio L, Watters J, Loboda A, Kulkarni A, Castle J, et al. Comparative expression pathway analysis of human and canine mammary tumors. BMC Genomics. (2009) 10:135–55. doi: 10.1186/1471-2164-10-135
9. Klopfleisch R, Lenze D, Hummel M, Gruber AD. Metastatic canine mammary carcinomas can be identified by a gene expression profile that partly overlaps with human breast cancer profiles. BMC Cancer. (2010) 10:618–29. doi: 10.1186/1471-2407-10-618
10. Queiroga FL, Raposo T, Carvalho MI, Prada J, Pires I. Canine mammary tumours as a model to study human breast cancer: most recent findings. In Vivo. (2011) 25:455–65.
11. Abadie J, Nguyen F, Loussouarn D, Peña L, Gama A, Rieder N, et al. Canine invasive mammary carcinomas as models of human breast cancer. Part 2: immunophenotypes and prognostic significance. Breast Cancer Res Treat. (2018) 167:459–68. doi: 10.1007/s10549-017-4542-8
12. Delaney G, Jacob S, Featherstone C, Barton M. The role of radiotherapy in cancer treatment: estimating optimal utilization from a review of evidence-based clinical guidelines. Cancer. (2005) 104:1129–37. doi: 10.1002/cncr.21324
13. Onitilo AA, Engel JM, Stankowski RV, Doi SA. Survival comparisons for breast conserving surgery and mastectomy revisited: community experience and the role of radiation therapy. Clin Med Res. (2015) 13:65–73. doi: 10.3121/cmr.2014.1245
14. Cao J, Olson R, Tyldesley S. Comparison of recurrence and survival rates after breast-conserving therapy and mastectomy in young women with breast cancer. Curr Oncol. (2013) 20:593–601. doi: 10.3747/co.20.1543
15. Poortmans P. Evidence based radiation oncology: breast cancer. Radiother Oncol. (2007) 84:84–101. doi: 10.1016/j.radonc.2007.06.002
16. Allemani C, Sant M, Weir HK, Richardson LC, Baili P, Storm H, et al. Breast cancer survival in the US and Europe: A CONCORD high-resolution study. Int J Cancer. (2013) 132:1170–81. doi: 10.1002/ijc.27725
17. Barker HE, Paget JTE, Khan AA, Harrington KJ. The tumour microenvironment after radiotherapy: mechanisms of resistance and recurrence. Nat Rev Cancer. (2015) 15:409–25. doi: 10.1038/nrc3958
18. Malhotra GK, Zhao X, Band H, Band V. Histological, molecular and functional subtypes of breast cancers. Cancer Biol Ther. (2010) 10:955–60. doi: 10.4161/cbt.10.10.13879
19. Gradishar WJ, Anderson BO, Balassanian R, Blair SL, Burstein HJ, Cyr A, et al. Invasive breast cancer version 1.2016, NCCN clinical practice guidelines in oncology. J Natl Compr Cancer Netw. (2016) 14:324–54. doi: 10.6004/jnccn.2016.0037
20. Harris L, Fritsche H, Mennel R, Norton L, Ravdin P, Taube S, et al. American society of clinical oncology 2007 update of recommendations for the use of tumor markers in breast cancer. J Clin Oncol. (2007) 25:5287–312. doi: 10.1200/JCO.2007.14.2364
21. Perou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, et al. Molecular portraits of human breast tumours. Nature. (2000) 406:747–52. doi: 10.1038/35021093
22. Sørlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, et al. Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci USA. (2001) 98:10869–74. doi: 10.1073/pnas.191367098
23. Sotiriou C, Neo SY, McShane LM, Korn EL, Long PM, Jazaeri A, et al. Breast cancer classification and prognosis based on gene expression profiles from a population-based study. Proc Natl Acad Sci USA. (2003) 100:10393–8. doi: 10.1073/pnas.1732912100
24. Sørlie T, Tibshirani R, Parker J, Hastie T, Marron J, Nobel A, et al. Repeated observation of breast tumor subtypes in independent gene expression data sets. Proc Natl Acad Sci USA. (2003) 100:8418–23. doi: 10.1073/pnas.0932692100
25. Prat A, Parker JS, Karginova O, Fan C, Livasy C, Herschkowitz JI, et al. Phenotypic and molecular characterization of the claudin-low intrinsic subtype of breast cancer. Breast Cancer Res. (2010) 12:68–75. doi: 10.1186/bcr2635
26. Goldschmidt M, Peña L, Rasotto R, Zappulli V. Classification and grading of canine mammary tumors. Vet Pathol. (2011) 48:117–31. doi: 10.1177/0300985810393258
27. Alenza MDP, Tabanera E, Peña L. Inflammatory mammary carcinoma in dogs: 33 cases (1995–1999). J Am Vet Med Assoc. (2001) 219:1110–4. doi: 10.2460/javma.2001.219.1110
28. Peña L, Perez-Alenza MD, Rodriguez-Bertos A, Nieto A. Canine inflammatory mammary carcinoma: histopathology, immunohistochemistry and clinical implications of 21 cases. Breast Cancer Res Treat. (2003) 78:141–8. doi: 10.1023/A:1022991802116
29. Schneider R, Dorn CR, Taylor D. Factors influencing canine mammary cancer development and postsurgical survival. J Natl Cancer Inst. (1969) 43:1249–61.
30. Kitchell B, Fidel J. Tamoxifen as a potential therapy for canine mammary carcinoma. Proc Vet Can Soc. (1992) 91:91–4.
31. Morris JS, Dobson JM, Bostock DE. Use of tamoxifen in the control of canine mammary neoplasia. Vet Rec. (1993) 133:539–42. doi: 10.1136/vr.133.22.539
32. Tavares WLF, Lavalle GE, Figueiredo MS, Souza AG, Bertagnolli AC, Viana FAB, et al. Evaluation of adverse effects in tamoxifen exposed healthy female dogs. Acta Vet. Scand. (2010) 52:67–74. doi: 10.1186/1751-0147-52-67
33. Andrew Novosad C. Principles of treatment for mammary gland tumors. Clin Tech Small Anim Pract. (2003) 18:107–9. doi: 10.1053/svms.2003.36625
34. Gray M, Turnbull AK, Ward C, Meehan J, Martínez-Pérez C, Bonello M, et al. Development and characterisation of acquired radioresistant breast cancer cell lines. Radiat Oncol. (2019) 14:64–83. doi: 10.1186/s13014-019-1268-2
35. HellmÉn E. Characterization of four in vitro established canine mammary carcinoma and one atypical benign mixed tumor cell lines. In Vitro Cell Dev Biol Anim. (1992) 28:309–19. doi: 10.1007/BF02877054
36. Wolfe LG, Smith BB, Toivio-Kinnucan MA, Sartin EA, Kwapien RP, Henderson RA, et al. Biologic properties of cell lines derived from canine mammary carcinomas23. J Natl Cancer Inst. (1986) 77:783–92. doi: 10.1093/jnci/77.3.783
37. van der Burg B, van Selm-Miltenburg AJP, van Maurik P, Rutteman GR, Misdorp W, de Laat SW, et al. Isolation of autonomously growing dog mammary tumor cell lines cultured in medium supplemented with serum treated to inactivate growth factors. J Natl Cancer Inst. (1989) 81:1545–51. doi: 10.1093/jnci/81.20.1545
38. Visan S, Balacescu O, Berindan-Neagoe I, Catoi C. In vitro comparative models for canine and human breast cancers. Clujul Med. (2016) 89:38. doi: 10.15386/cjmed-519
39. Else RW, Norval M, Neill WA. The characteristics of a canine mammary carcinoma cell line, REM 134. Br J Cancer. (1982) 46:675–81. doi: 10.1038/bjc.1982.254
40. Franken NA, Rodermond HM, Stap J, Haveman J, Van Bree C. Clonogenic assay of cells in vitro. Nat Protoc. (2006) 1:2315. doi: 10.1038/nprot.2006.339
41. Liang CC, Park AY, Guan JL. In vitro scratch assay: a convenient and inexpensive method for analysis of cell migration in vitro. Nat Protoc. (2007) 2:329–33. doi: 10.1038/nprot.2007.30
42. Meehan J, Ward C, Turnbull A, Bukowski-Wills J, Finch AJ, Jarman EJ, et al. Inhibition of pH regulation as a therapeutic strategy in hypoxic human breast cancer cells. Oncotarget. (2017) 8:42857–75. doi: 10.18632/oncotarget.17143
43. Bankhead P, Loughrey MB, Fernandez JA, Dombrowski Y, McArt DG, Dunne PD, et al. QuPath: open source software for digital pathology image analysis. Sci Rep. (2017) 7:16878–85. doi: 10.1038/s41598-017-17204-5
44. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. (2013) 29:15–21. doi: 10.1093/bioinformatics/bts635
45. Kinsella RJ, Kähäri A, Haider S, Zamora J, Proctor G, Spudich G, et al. Ensembl BioMarts: a hub for data retrieval across taxonomic space. Database. (2011) 2011:bar030. doi: 10.1093/database/bar030
46. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. (2004) 5:80–96. doi: 10.1186/gb-2004-5-10-r80
47. Howe E, Holton K, Nair S, Schlauch D, Sinha R, Quackenbush J. MeV: multiexperiment viewer. In: Ochs MF, Casagrande JT, Davuluri RV, editors. Biomedical Informatics for Cancer Research. Boston, MA: Springer (2010). p. 267–77.
48. Leek JT, Johnson WE, Parker HS, Jaffe AE, Storey JD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. (2012) 28:882–3. doi: 10.1093/bioinformatics/bts034
49. Turnbull AK, Kitchen RR, Larionov AA, Renshaw L, Dixon JM, Sims AH. Direct integration of intensity-level data from Affymetrix and Illumina microarrays improves statistical power for robust reanalysis. BMC Med Genomics. (2012) 5:5–35. doi: 10.1186/1755-8794-5-35
50. Sherman BT, Tan Q, Collins JR, Alvord WG, Roayaei J, Stephens R, et al. The DAVID Gene Functional Classification Tool: a novel biological module-centric algorithm to functionally analyze large gene lists. Genome Biol. (2007) 8:183. doi: 10.1186/gb-2007-8-9-r183
51. Gendoo DM, Ratanasirigulchai N, Schröder MS, Paré L, Parker JS, Prat A, et al. Genefu: an R/Bioconductor package for computation of gene expression-based signatures in breast cancer. Bioinformatics. (2015) 32:1097–9. doi: 10.1093/bioinformatics/btv693
52. Edgar R, Domrachev M, Lash AE. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. (2002) 30:207–10. doi: 10.1093/nar/30.1.207
53. Rokavec M, Kaller M, Horst D, Hermeking H. Pan-cancer EMT-signature identifies RBM47 down-regulation during colorectal cancer progression. Sci Rep. (2017) 7:4687–702. doi: 10.1038/s41598-017-04234-2
54. Theys J, Jutten B, Habets R, Paesmans K, Groot AJ, Lambin P, et al. E-Cadherin loss associated with EMT promotes radioresistance in human tumor cells. Radiother Oncol. (2011) 99:392–7. doi: 10.1016/j.radonc.2011.05.044
55. Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. (1999) 27:29–34. doi: 10.1093/nar/27.1.29
56. Verhaegh W, van Ooijen H, Inda MA, Hatzis P, Versteeg R, Smid M, et al. Selection of personalized patient therapy through the use of knowledge-based computational models that identify tumor-driving signal transduction pathways. Cancer Res. (2014) 74:2936–45. doi: 10.1158/0008-5472.CAN-13-2515
57. Creighton CJ, Hilger AM, Murthy S, Rae JM, Chinnaiyan AM, El-Ashry D. Activation of mitogen-activated protein kinase in estrogen receptor α-positive breast cancer cells in vitro induces an in vivo molecular phenotype of estrogen receptor α-negative human breast tumors. Cancer Res. (2006) 66:3903–11. doi: 10.1158/0008-5472.CAN-05-4363
58. van Ooijen H, Hornsveld M, Dam-de Veen C, Velter R, Dou M, Verhaegh W, et al. Assessment of functional phosphatidylinositol 3-kinase pathway activity in cancer tissue using forkhead box O target gene expression in a knowledge-based computational model. Am. J. Pathol. (2018) 188:1956–72. doi: 10.1016/j.ajpath.2018.05.020
59. Nguyen PL, Taghian AG, Katz MS, Niemierko A, Abi Raad RF, Boon WL, et al. Breast cancer subtype approximated by estrogen receptor, progesterone receptor, and HER-2 is associated with local and distant recurrence after breast-conserving therapy. J Clin Oncol. (2008) 26:2373–8. doi: 10.1200/JCO.2007.14.4287
60. Stål O, Sullivan S, Wingren S, Skoog L, Rutqvist L, Carstensen J, et al. c-erbB-2 expression and benefit from adjuvant chemotherapy and radiotherapy of breast cancer. Eur J Cancer. (1995) 31:2185–90. doi: 10.1016/0959-8049(95)00344-4
61. Kyndi M, Sørensen FB, Knudsen H, Overgaard M, Nielsen HM, Overgaard J. Estrogen receptor, progesterone receptor, HER-2, and response to postmastectomy radiotherapy in high-risk breast cancer: the Danish breast cancer cooperative group. J Clin Oncol. (2008) 26:1419–26. doi: 10.1200/JCO.2007.14.5565
62. Speers C, Zhao S, Liu M, Bartelink H, Pierce LJ, Feng FY. Development and Validation of a novel radiosensitivity signature in human breast cancer. Clin Cancer Res. (2015) 21:3667–77. doi: 10.1158/1078-0432.CCR-14-2898
63. Smith L, Qutob O, Watson MB, Beavis AW, Potts D, Welham KJ, et al. Proteomic identification of putative biomarkers of radiotherapy resistance: a possible role for the 26S proteasome? Neoplasia. (2009) 11:1194–207. doi: 10.1593/neo.09902
64. Li Z, Xia L, Lee LM, Khaletskiy A, Wang J, Wong JY, et al. Effector genes altered in MCF-7 human breast cancer cells after exposure to fractionated ionizing radiation. Radiat Res. (2001) 155:543–53. doi: 10.1667/0033-7587(2001)155[0543:EGAIMH]2.0.CO;2
65. Hirschhaeuser F, Menne H, Dittfeld C, West J, Mueller-Klieser W, Kunz-Schughart LA. Multicellular tumor spheroids: an underestimated tool is catching up again. J Biotechnol. (2010) 148:3–15. doi: 10.1016/j.jbiotec.2010.01.012
66. Raleigh J, Chou S, Arteel G, Horsman M. Comparisons among pimonidazole binding, oxygen electrode measurements, and radiation response in C3H mouse tumors. Radiat Res. (1999) 151:580–9. doi: 10.2307/3580034
67. Bussink J, Kaanders JH, Strik AM, van der Kogel AJ. Effects of nicotinamide and carbogen on oxygenation in human tumor xenografts measured with luminescense based fiber-optic probes. Radiother Oncol. (2000) 57:21–30. doi: 10.1016/S0167-8140(00)00275-9
68. Ljungkvist AS, Bussink J, Rijken PF, Raleigh JA, Denekamp J, Van Der Kogel AJ. Changes in tumor hypoxia measured with a double hypoxic marker technique. Int J Radiat Oncol Biol Phys. (2000) 48:1529–38. doi: 10.1016/S0360-3016(00)00787-2
69. Chang L, Ni J, Beretov J, Wasinger VC, Hao J, Bucci J, et al. Identification of protein biomarkers and signaling pathways associated with prostate cancer radioresistance using label-free LC-MS/MS proteomic approach. Sci Rep. (2017) 7:41834–48. doi: 10.1038/srep41834
70. Gilbertson S, Kurzman I, Zachrau R, Hurvitz A, Black M. Canine mammary epithelial neoplasms: biologic implications of morphologic characteristics assessed in 232 dogs. Vet Pathol. (1983) 20:127–42. doi: 10.1177/030098588302000201
71. Rosol TJ, Tannehill-Gregg SH, LeRoy BE, Mandl S, Contag CH. Animal models of bone metastasis. Cancer. (2003) 97:748–57. doi: 10.1002/cncr.11150
72. O'Shaughnessy J. Extending survival with chemotherapy in metastatic breast cancer. Oncologist. (2005) 10:20–9. doi: 10.1634/theoncologist.10-90003-20
73. Atkins M, Potier D, Romanelli L, Jacobs J, Mach J, Hamaratoglu F, et al. An ectopic network of transcription factors regulated by hippo signaling drives growth and invasion of a malignant tumor model. Curr Biol. (2016) 26:2101–13. doi: 10.1016/j.cub.2016.06.035
74. Imajo M, Miyatake K, Iimura A, Miyamoto A, Nishida E. A molecular mechanism that links Hippo signalling to the inhibition of Wnt/beta-catenin signaling. EMBO J. (2012) 31:1109–1022. doi: 10.1038/emboj.2011.487
75. Wei C, Wang Y, Li X. The role of Hippo signal pathway in breast cancer metastasis. Onco Targets Ther. (2018) 11:2185–93. doi: 10.2147/OTT.S157058
76. Pinto CA, Widodo E, Waltham M, Thompson EW. Breast cancer stem cells and epithelial mesenchymal plasticity – implications for chemoresistance. Cancer Lett. (2013) 341:56–62. doi: 10.1016/j.canlet.2013.06.003
77. Lee JM, Dedhar S, Kalluri R, Thompson EW. The epithelial–mesenchymal transition: new insights in signaling, development, and disease. J Cell Biol. (2006) 172:973–81. doi: 10.1083/jcb.200601018
78. Thompson EW, Newgreen DF. Carcinoma invasion and metastasis: a role for epithelial-mesenchymal transition? Cancer Res. (2005) 65:5991–5. doi: 10.1158/0008-5472.CAN-05-0616
79. Sommers CL, Heckford SE, Skerker JM, Worland P, Torri JA, Thompson EW, et al. Loss of epithelial markers and acquisition of vimentin expression in adriamycin-and vinblastine-resistant human breast cancer cell lines. Cancer Res. (1992) 52:5190–7.
80. Liu F, Gu LN, Shan BE, Geng CZ, Sang MX. Biomarkers for EMT and MET in breast cancer: an update. Oncol Lett. (2016) 12:4869–76. doi: 10.3892/ol.2016.5369
81. Liao TT, Yang M-H. Hybrid epithelial/mesenchymal state in cancer metastasis: clinical significance and regulatory mechanisms. Cells. (2020) 9:623–30. doi: 10.3390/cells9030623
82. Sha Y, Haensel D, Gutierrez G, Du H, Dai X, Nie Q. Intermediate cell states in epithelial-to-mesenchymal transition. Phys Biol. (2019) 16:021001. doi: 10.1088/1478-3975/aaf928
83. Yook JI, Li XY, Ota I, Fearon ER, Weiss SJ. Wnt-dependent regulation of the E-cadherin repressor snail. J Biol Chem. (2005) 280:11740–8. doi: 10.1074/jbc.M413878200
84. Thiery JP, Acloque H, Huang RY, Nieto MA. Epithelial-mesenchymal transitions in development and disease. Cell. (2009) 139:871–90. doi: 10.1016/j.cell.2009.11.007
85. Zhou YC, Liu JY, Li J, Zhang J, Xu YQ, Zhang HW, et al. Ionizing radiation promotes migration and invasion of cancer cells through transforming growth factor-beta–mediated epithelial–mesenchymal transition. Int J Radiat Oncol Biol Phys. (2011) 81:1530–7. doi: 10.1016/j.ijrobp.2011.06.1956
86. Andarawewa KL, Erickson AC, Chou WS, Costes SV, Gascard P, Mott JD, et al. Ionizing radiation predisposes nonmalignant human mammary epithelial cells to undergo transforming growth factor beta induced epithelial to mesenchymal transition. Cancer Res. (2007) 67:8662–78. doi: 10.1158/0008-5472.CAN-07-1294
87. Klopfleisch R, Klose P, Weise C, Bondzio A, Multhaup G, Einspanier R, et al. Proteome of metastatic canine mammary carcinomas: similarities to and differences from human breast cancer. J Proteome Res. (2010) 9:6380–91. doi: 10.1021/pr100671c
88. Wang Y, Yin Q, Yu Q, Zhang J, Liu Z, Wang S, et al. A retrospective study of breast cancer subtypes: the risk of relapse and the relations with treatments. Breast Cancer Res Treat. (2011) 130:489–98. doi: 10.1007/s10549-011-1709-6
89. Dai X, Li T, Bai Z, Yang Y, Liu X, Zhan J, et al. Breast cancer intrinsic subtype classification, clinical use and future trends. Am J Cancer Res. (2015) 5:2929–43. doi: 10.1371/journal.pone.0124964
90. Worzfeld T, Swiercz JM, Looso M, Straub BK, Sivaraj KK, Offermanns S. ErbB-2 signals through Plexin-B1 to promote breast cancer metastasis. J Clin Invest. (2012) 122:1296–305. doi: 10.1172/JCI60568
91. Gutierrez C, Schiff R. HER2: biology, detection, and clinical implications. Arch Pathol Lab Med. (2011) 135:55–62. doi: 10.1043/2010-0454-RAR.1
92. Ross JS, Fletcher JA. HER-2/neu (c-erb-B2) gene and protein in breast cancer. Am J Clin Pathol. (1999) 112 (1 Suppl. 1):53–67.
93. Ross JS, Fletcher JA, Linette GP, Stec J, Clark E, Ayers M, et al. The Her-2/neu gene and protein in breast cancer 2003: biomarker and target of therapy. Oncologist. (2003) 8:307–25. doi: 10.1634/theoncologist.8-4-307
94. Dutra A, Granja N, Schmitt F, Cassali G. c-erbB-2 expression and nuclear pleomorphism in canine mammary tumors. Braz J Med Biol Res. (2004) 37:1673–81. doi: 10.1590/S0100-879X2004001100013
95. Hsu WL, Huang HM, Liao JW, Wong ML, Chang SC. Increased survival in dogs with malignant mammary tumours overexpressing HER-2 protein and detection of a silent single nucleotide polymorphism in the canine HER-2 gene. Vet J. (2009) 180:116–23. doi: 10.1016/j.tvjl.2007.10.013
96. Muhammadnejad A, Keyhani E, Mortazavi P, Behjati F, Haghdoost IS. Overexpression of her-2/neu in malignant mammary tumors; translation of clinicopathological features from dog to human. Asian Pac J Cancer Prev. (2012) 13:6415–21. doi: 10.7314/APJCP.2012.13.12.6415
97. Ahern T, Bird R, Bird A, Wolfe L. Expression of the oncogene c-erbB-2 in canine mammary cancers and tumor-derived cell lines. Am J Vet Res. (1996) 57:693–6.
98. Silva I, Dias A, Bertagnolli A, Cassali G, Ferreira E. Analysis of EGFR and HER-2 expressions in ductal carcinomas in situ in canine mammary glands. Arquivo Brasileiro Med Vet Zootec. (2014) 66:763–8. doi: 10.1590/1678-41626128
99. Singer J, Weichselbaumer M, Stockner T, Mechtcheriakova D, Sobanov Y, Bajna E, et al. Comparative oncology: ErbB-1 and ErbB-2 homologues in canine cancer are susceptible to cetuximab and trastuzumab targeting. Mol Immunol. (2012) 50:200–9. doi: 10.1016/j.molimm.2012.01.002
100. Nicholson RI, McClelland RAJ, Gee MW, Manning DL, Cannon P, Robertson JFR, et al. Epidermal growth factor receptor expression in breast cancer: association with response to endocrine therapy. Breast Cancer Res Treat. (1994) 29:117–125. doi: 10.1007/BF00666187
101. Nicholson RI, McClelland RA, Finlay P, Eaton CL, Gullick WJ, Dixon AR, et al. Relationship between EGF-R, c-erbB-2 protein expression and Ki67 immunostaining in breast cancer and hormone sensitivity. Eur J Cancer. (1993) 29:1018–23. doi: 10.1016/S0959-8049(05)80215-1
102. Wright C, Nicholson S, Angus B, Sainsbury J, Farndon J, Cairns J, et al. Relationship between c-erbB-2 protein product expression and response to endocrine therapy in advanced breast cancer. Br J Cancer. (1992) 65:118–21. doi: 10.1038/bjc.1992.22
103. Shimura T, Noma N, Oikawa T, Ochiai Y, Kakuda S, Kuwahara Y, et al. Activation of the AKT/cyclin D1/Cdk4 survival signaling pathway in radioresistant cancer stem cells. Oncogenesis. (2012) 1:12–20. doi: 10.1038/oncsis.2012.12
104. Shimura T, Kakuda S, Ochiai Y, Nakagawa H, Kuwahara Y, Takai Y, et al. Acquired radioresistance of human tumor cells by DNA-PK/AKT/GSK3β-mediated cyclin D1 overexpression. Oncogene. (2010) 29:4826–37. doi: 10.1038/onc.2010.238
105. Igor V, Charles LS. The phosphatidylinositol 3-Kinase–AKT pathway in human cancer. Nat Rev Cancer. (2002) 2:489–501. doi: 10.1038/nrc839
106. Liang K, Jin W, Knuefermann C, Schmidt M, Mills GB, Ang KK, et al. Targeting the phosphatidylinositol 3-kinase/Akt pathway for enhancing breast cancer cells to radiotherapy. Mol Cancer Ther. (2003) 2:353–60.
107. Kavanagh BD, Dent P, Schmidt-Ullrich RK, Chen P, Mikkelsen RB. Calcium-dependent stimulation of mitogen-activated protein kinase activity in A431 cells by low doses of ionizing radiation. Radiat Res. (1998) 149:579–87. doi: 10.2307/3579904
108. Suy S, Anderson WB, Dent P, Chang E, Kasid U. Association of Grb2 with Sos and Ras with Raf-1 upon gamma irradiation of breast cancer cells. Oncogene. (1997) 15:53–61. doi: 10.1038/sj.onc.1201165
109. Kavanagh BD, Lin PS, Chen P, Schmidt-Ullrich RK. Radiation-induced enhanced proliferation of human squamous cancer cells in vitro: a release from inhibition by epidermal growth factor. Clin Cancer Res. (1995) 1:1557–62.
110. Schmidt-Ullrich RK, Mikkelsen RB, Dent P, Todd DG, Valerie K, Kavanagh BD. Radiation-induced proliferation of the human A431 squamous carcinoma cells is dependent on EGFR tyrosine phosphorylation. Oncogene. (1997) 15:1191–7. doi: 10.1038/sj.onc.1201275
Keywords: canine breast cancer models, human breast cancer, radioresistance, global gene analysis, characterization of radioresistant cell lines, comparative oncology
Citation: Gray M, Turnbull AK, Meehan J, Martínez-Pérez C, Kay C, Pang LY and Argyle DJ (2020) Comparative Analysis of the Development of Acquired Radioresistance in Canine and Human Mammary Cancer Cell Lines. Front. Vet. Sci. 7:439. doi: 10.3389/fvets.2020.00439
Received: 12 May 2020; Accepted: 16 June 2020;
Published: 23 July 2020.
Edited by:
Laura Pena, Complutense University of Madrid, SpainReviewed by:
Laura Cavicchioli, University of Padova, ItalyAlejandro Suarez-Bonnet, Royal Veterinary College (RVC), United Kingdom
Copyright © 2020 Gray, Turnbull, Meehan, Martínez-Pérez, Kay, Pang and Argyle. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mark Gray, mark.gray@.ed.ac.uk
†These authors have contributed equally to this work