Metastatic canine mammary carcinomas can be identified by a gene expression profile that partly overlaps with human breast cancer profiles
BMC Cancer volume 10, Article number: 618 (2010)
Similar to human breast cancer mammary tumors of the female dog are commonly associated with a fatal outcome due to the development of distant metastases. However, the molecular defects leading to metastasis are largely unknown and the value of canine mammary carcinoma as a model for human breast cancer is unclear. In this study, we analyzed the gene expression signatures associated with mammary tumor metastasis and asked for parallels with the human equivalent.
Messenger RNA expression profiles of twenty-seven lymph node metastasis positive or negative canine mammary carcinomas were established by microarray analysis. Differentially expressed genes were functionally characterized and associated with molecular pathways. The findings were also correlated with published data on human breast cancer.
Metastatic canine mammary carcinomas had 1,011 significantly differentially expressed genes when compared to non-metastatic carcinomas. Metastatic carcinomas had a significant up-regulation of genes associated with cell cycle regulation, matrix modulation, protein folding and proteasomal degradation whereas cell differentiation genes, growth factor pathway genes and regulators of actin organization were significantly down-regulated. Interestingly, 265 of the 1,011 differentially expressed canine genes are also related to human breast cancer and, vice versa, parts of a human prognostic gene signature were identified in the expression profiles of the metastatic canine tumors.
Metastatic canine mammary carcinomas can be discriminated from non-metastatic carcinomas by their gene expression profiles. More than one third of the differentially expressed genes are also described of relevance for human breast cancer. Many of the differentially expressed genes are linked to functions and pathways which appear to be relevant for the induction and maintenance of metastatic progression and may represent new therapeutic targets. Furthermore, dogs are in some aspects suitable as a translational model for human breast tumors in order to identify prognostic molecular signatures and potential therapeutic targets.
Canine mammary tumor (CMT) is the most common cancer among female dogs and often becomes fatal due to the development of distant metastases [1–3]. Metastasis to the regional lymph node is an early step in metastasis and one of the most important prognostic factors in the diagnosis of CMT, a criterion that is also valid for human breast cancer [4, 5]. Lymph node metastases of CMT are usually followed by the development of distant metastases, mainly in the lung, ultimately leading to the death of the dog . However, knowledge of the molecular mechanisms contributing to lymph node and distant metastasis is still fragmentary. Despite numerous studies on this issue, significant metastasis-associated and predictable expression patterns of single genes have not been identified in CMT as yet [7–9]. Global gene expression profiles that compare metastasizing versus non-metastasizing CMT are unavailable whereas several studies on human breast cancer found significant metastasis associated expression profiles. The latter studies identified several non-overlapping expression signatures which are related to the development of lymph node and distant metastases and worse prognoses [10–13]. The available studies on global gene expression in CMT compared normal mammary gland, benign and malignant tumors with unknown lymph node status and clinical follow-up [14, 15]. The authors reported that the gene expression profiles of CMT include a gene expression signature associated with neoplastic transformation. Furthermore, comparison of canine and human expression profiles disclosed an overlap of deregulated genes in human and canine mammary tumors . Generally, clinical and molecular features of human and canine bear a likeness in several aspects. Both malignancies are the most common cancer of the female, lymph node metastases indicate a poor prognosis, the hormonal status influences the development of CMT and estrogen receptor (ER), progesterone receptor (PR) and ERBB2 expression patterns do influence the overall survival rate [1, 16–18].
The aim of this study was the identification of gene expression signatures in primary CMT that are associated with early lymph node metastasis. Global mRNA expression profiles obtained from metastatic versus non-metastatic CMT cases were compared and differentially expressed genes were analyzed for their function and their role in pathway activation. Moreover, mining in published literature revealed interesting overlapping features when compared to data derived from gene expression profiles of metastatic human breast carcinomas.
Thirteen simple mammary carcinomas with invasive growth and lymph node metastases at the time of tumor resection and 14 simple carcinomas without lymph node metastases were included in the study (Table 1). Complex carcinomas were excluded from the study to avoid differences in gene expression levels due to differences in mesenchymal/epithelial ratio in the different tissues. None of the patients had a history of progestin treatment or radiographically detectable pulmonary metastases at the time of tumor resection. Distant metastases as the cause of death were determined postoperatively by radiographic detection of metastases (nos. 1-11) or necropsy (nos. 12 and 13). Selection criteria for carcinomas without lymph node metastases included an invasive growth, a negative lymph node status, a histological grade III and a minimal tumor diameter above the average of the lymph node positive tumors (> 2.42 cm). All animals with non-metastatic carcinomas had a overall survival rate of over 24 months except animal no. 22 which developed radiographic detectable lung metastases 8 months after surgery.
The experimental research reported in the manuscript has been performed with the approval of animal welfare authorities and the ethical committee of the Freie Universität Berlin. Surgical excision of tumour biopsies was part of the tumour treatment according to the state of the art treatment and solely to improve the animals welfare. Furthermore, the animals were under full anaesthesia and not exposed to any additional manipulation due to the inclusion in this study. All animal owners received and approved an informed client consent form.
Tissue specimens were fixed in neutral-buffered 4% formalin or snap frozen in liquid nitrogen within 15 minutes after resection and stored at -80°C until further use. Formalin fixed tumor tissue samples were routinely embedded in paraffin and sections of 2-μm were stained with hematoxylin and eosin. Tumor and lymph node histologies were evaluated independently by two board-certified pathologists, following the criteria of the WHO classification of canine mammary tumors and the Nottingham grading system [19, 20]. All 27 tumors were simple carcinomas and characterized by an invasive, mostly solid growth pattern, marked cellular pleomorphism, anisokaryosis and 3 or more mitotic figures per high power field.
Oestrogen receptor alpha (ER), ERBB2 expression was immunohistochemically determined using the ABC-method. In brief, monoclonal mouse anti-human ER specific antibody (1:1000, clone CC4-5, Novocastra, Wetzlar, Germany) and rabbit polyclonal anti-human ERBB2 specific antibody (1:150, cat. no. A0485, Dako, Hamburg, Germany) were diluted in Tris-buffered saline (TBS, 50 mM, pH 7.6) and incubated at 4°C overnight after a blocking step with 50% goat serum in TBS for 30 min at room temperature. Polyclonal goat anti-rabbit IgG (1:200; Vector, England, BA1000) and goat anti-mouse IgG (1:200; BA9200, Vector, Burlingame, USA) were used as secondary antibody. Diaminobenzidine tetrahydrochloride (D8001, Sigma Aldrich, Munich, Germany) was used as chromogen and slides were counterstained with hematoxylin (Merck). Normal canine and human mammary gland were used as positive tissue controls. ER and ERBB2 immunolabeling was evaluated in 10 random 200× magnification fields. Tumors were defined as ER or ERBB2 positive if more than 10% of the cells stained positive with the respective antibody.
Macrodissection of tumor samples and RNA isolation
Macrodissection was performed on all tumor specimens to ensure high tumor cellularity. Only sections with more than 70% carcinoma cells were included in the study as shown by digital image analysis (Scanscope T3, Aperio, Vista, USA; Zeiss Axiovision, Jena, Germany). For mRNA isolation, tissue sections were transferred into 300 μl of lysis buffer (NucleoSpin RNA; Macherey&Nagel, Dueren, Germany) and homogenized (Precellyse 24, Bertin Technology, France). RNA was extracted and purified using a commercial system (NucleoSpin RNA; Macherey&Nagel, Düren, Germany). The RNA quality was controlled using the BioAnalyzer (Agilent Technologies, USA) and only high quality RNA (RIN >8) was used for microarray analyses. For reverse transcription a commercial kit (Iscript cDNA synthesis Kit, Biorad, Germany) was used. 100 ng of total RNA of each sample were used for reverse transcription of each tissue sample and the reactions were performed exactly according to the manufacturers instructions. cDNA was stored in low adhesion tubes at -20°C.
Affymetrix GeneChip hybridization (Canine Genome 2.0 Array) was performed with 2 μg total RNA according to the manufacturer's recommendations. The chips were stained and washed with the GeneChip Fluidics Station 450 and visualized on an Affymetrix GeneChip Scanner 3000. Microarray data were deposited at the Gene Expression Omnibus data repository under the number GSE20718: (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=bbwvhacsemokkfq&acc=GSE20718)
Quantitative RT-PCR (qPCR)
Quantitative PCR was used to validate array gene expression data. To this end, established qPCR assays for eight genes that were identified in the metastatic gene signature were used to compare mRNA expression levels in metastatic and non-metastatic tumors. Primer sequences, product length and sequence accession numbers for the PCR assays for AURKA, ALOX12, BMP6, ERBB4, HEPACAM2, IGFR2, RAD51, TGFBR3 and the housekeeper genes A5B, HPRT and RP32 mRNA expression levels were used before [17, 21–24] and are describe in Additional file 1. Each qPCR reaction had a 20 μl reaction volume containing 5 μl cDNA corresponding to 100 ng input RNA, 500 nM of each forward and reverse primer, 12.5 μl of the Maxima Sybr Green/Rox qPCR Master Mix (Fermentas, Germany) containing 2.5 mM Mg2+ and 5 mM dNTP. Tubes and lids were purchased from Fermentas (Germany). Cycling parameters for all assays were, 95°C for 10 min, 95°C for 30 sec, 58°C for 60 sec, 72°C for 30 sec for 40 cycles. Melting curve ranged from 58°C to 95°C and was read every 0.2°C and hold for 2 sec. The cDNA of all samples were amplified on the same plate for every primer pair to ensure equal amplification conditions. Specificity of amplification products was confirmed by melting curve and sequence analyses. For each sample, results were documented as cycle threshold (CT) set to 100 relative fluorescence unit values of background subtracted qPCR fluorescence kinetics by using the MX Pro Stratagene analysis software, applying the adaptive baseline, amplification based threshold and moving average algorithm enhancement. Relative expression of the target gene (TG) was calculated using a comparative ΔCT-method with multiple housekeepers as previously described[25–27]. The housekeeper genes used were selected from a panel of reference genes (RG) according to the GeNorm algorithm. Data are presented as fold change in gene expression level of the gene of interest (GOI) in lymph node-positive carcinomas (LN+) normalized to the housekeepers and to the similarly normalized GOI expression levels in lymph node-negative carcinomas (LN-). Significance (p < 0.05) of differential gene expression was tested using Shapiro-Wilk-test and t-test.
Affymetrix CEL files were imported into Partek Genomic Suite Software (Version 6.4, Partek Inc., St. Louis, USA) and processed by the implemented gcRMA workflow (median polish probe set summarization, RMA background correction, quantile normalization). Differences in gene expression between samples with and without metastases were analyzed by ANOVA and false discovery rate was controlled by using the q-value method. Differentially expressed genes were selected by applying a filter of q < 0.001 and a fold-change of >1.7 in both directions. Unnamed genes were excluded from the list. Hierarchical clustering of the samples and genes was conducted using Pearson correlation and complete linkage.
To supplement the gene annotations of the differentially expressed genes with functional information, BLAST search and Affymetrix-provided human to canine microarray comparisons were used to map canine genes to their human equivalents as has been shown before . Using the human equivalents as templates, the DAVID database was queried for gene ontology information . To study enriched functional gene families and functional annotation, all down-regulated and all up-regulated genes were submitted separately to DAVID. In the case of redundant probe with a fold change in the same direction only the probe set with the highest fold change was included in further analysis. Selection criteria for DAVID included a medium stringency, ≥ 4 probes within a cluster and an enrichment factor >1.3. In case of multiple appearances of similar gene families or functional annotation terms, the cluster with the higher enrichment factor was selected. Reported associations between the differentially expressed genes and human breast cancer were determined using the MEDGENE database , entering the search term: 'breast neoplasms'.
The GenMAPP software was used for pathway analysis . Analyses of up- or down-regulated gene sets were performed independently. Quality criteria for pathway selection included more than 3 genes, z-score of |>2| and p < 0.05.
Immunohistological phenotype of the tumors
All lymph node positive and lymph node negative tumors were ER negative. In addition, 11/13 lymph node positive tumors and 13/14 lymph node negative tumors were ERBB2 positive.
Identification of differentially expressed genes
Comparison of the global expression of canine mammary carcinomas with or without lymph node metastasis revealed 1011 differentially expressed genes, 744 genes of which were up-regulated metastatic carcinomas whereas only 267 genes were down-regulated (1.7× fold change, false discovery rate q < 0.001, Additional file 2). With the exception of two cases (dogs no. 14, 22), these genes were able to clearly separate the two tumor groups by means of two-dimensional hierarchical clustering (Figure 1). Interestingly, one of the tumors (no. 22) developed radiographic detectable lung metastases 8 months postoperatively. However, the other case (14) did not develop metastases during a 24 month clinical follow-up. qPCR confirmed significant expression differences for all seven of the eight genes analyzed (Table 2).
Up-regulation of cell cycle and DNA-damage repair genes in metastatic canine carcinomas
Functional analysis of the metastasis-associated genes identified an up-regulation of 93 cell cycle associated genes with functions in cell cycle progression (Table 3). Interestingly, the 13 cell cycle checkpoint genes and 27 DNA-damage repair genes in the gene list were exclusively up-regulated and only one of the 6 cyclin genes was down-regulated (Table 3). 49 genes associated with protein metabolism were also up-regulated, including protein folding, ubiquitin dependent proteins and prefoldins, suggesting an increased biosynthesis of the cells. Furthermore, anti-apoptotic and pro-apoptotic genes, members of the RAS signaling cascade were up-regulated (Table 4). Interestingly, transcription factors and methyltransferases were significantly enriched in both up-regulated and down-regulated gene lists indicating that the metastatic phenotype is associated with a different equilibrium of gene activation and deactivation (Table 4).
Down-regulation of growth factor receptor pathways and cell differentiation genes in metastatic canine carcinomas
39 genes with receptor activity were up- and down regulated in metastatic tumors (Table 4, 5). However, transmembrane tyrosine kinase receptor genes were significantly and exclusively down-regulated in metastatic carcinomas (Table 5). Specifically, genes of growth factor signaling pathways epidermal growth factor (EGF), transforming growth factor-beta (TGFB), fibroblast growth factor (FGF), hepatocyte growth factor receptor (MET) and insulin-like growth factor (IGF) I and II were preferentially down-regulated in the metastatic tumors. Besides their growth stimulatory effects, EGF, TGF and FGF signaling also induces cell differentiation in mammary epithelial cells and a down regulation of these pathways and associated proteins has been described before [24, 33–37].
In line with this changed receptor signaling in metastatic carcinomas is therefore the de-regulation of several genes associated with cell growth regulation and cell differentiation in the metastatic tumors (Table 5). Interestingly, statistical analysis by DAVID gene ontology identified a significant enrichment of cell differentiation genes only in the list of down-regulated genes (Tables 4).
Cell-cell adhesion and matrix modulation genes are differentially expressed in metastatic tumors
It is commonly accepted that metastasis requires loss of cell-cell adhesion and expression of extracellular matrix modulators by the tumor cells. We, therefore, specifically screened the gene list for modulators of the extracellular matrix, cell adhesion and angiogenesis genes. Statistical analysis found no significant enrichment of matrix modulator genes in the gene list (Table 4). However, three matrix modulator genes: matrix-metalloproteinase 11, cathepsin L2 and the pro-cancerogenic serin peptidase inhibitor E1, were up-regulated in metastatic tumors (Table 6) and this is in agreement with earlier observations in human breast cancer [38–40]. Similarly, cell adhesion genes were also not significantly enriched in the gene list. However, the 20 adhesion and focal adhesion complex genes present in the differential gene list were predominantly down-regulated, like HEPACAM2, EMILIN1, PCDHGC3 and ASAM but some adhesion genes were also up-regulated, like ICAM2 and the ITGA4 (Table 6). A similar expression pattern for several of these genes has been described for human breast cancer and other cancer types before (Table 6) [23, 41–44].
Finally, analysis of the gene list disclosed no significant enrichment of angiogenesis genes. The 7 angiogenesis genes present in the list were also predominantly down-regulated (e.g. NRP1, JMJD6, CYR61, FGFR1) and again the impact of these genes on the carcinogenesis of breast cancer cells has been described[36, 45, 46].
Comparison with human breast cancer expression profiles
The identified metastasis-associated expression profile of canine mammary tumors has several significant overlaps with expression profiles of metastatic human breast cancer. Literature mining by Medgene disclosed that 25% of the canine gene list has been cited in association with human breast cancer before. However, the mining approach by Medgene is based on information presented in the abstract of the articles and not complete gene expression profiles . It is therefore likely that a comparative meta-analysis of gene expression profiles will identify additional overlapping gene expression patterns. The overlapping genes are mostly associated with cell division, cell growth, kinase activity, translation, DNA integrity checkpoint and transcription regulation.
GenMapp pathway analysis identified a marked overlap of deregulated pathways in the canine metastatic profile and a recent meta-analysis on metastasis-associated gene expression profiles (Table 7) . All 5 significantly up-regulated pathways were Identical to human breast cancer the cell cycle, G1 to S control reactome and DNA replication reactome pathways were up-regulated, whereas striated muscle contraction and IL-6 pathways were significantly down-regulated. (Table 7).
Furthermore, comparison with the van't Veer 70 gene prognostic signature which reliably identifies human breast cancers with metastatic potential also found a significant overlap . However, only 34 of the 70 genes are present on the canine gene array but 23 (68%) of these genes were included in the metastatic gene expression profile of CMT indicating that comparable prognostic gene expression patterns may exist in canine mammary tumors (Table 8).
In the present study we identified a gene expression profile in CMT that is associated with early metastatic spread to the lymph nodes and is able to discriminate carcinomas with similar histological features but divergent metastatic potential. The differential expression profile contains several enriched functional gene classes and has significant overlaps with expression profiles of metastatic human breast cancer.
The similarities found between canine and human mammary tumors in terms of increased proliferation, altered cell differentiation status and decreased cell adhesion were also confirmed by Medgene literature mining. Approximately 25% of the deregulated genes in metastatic canine carcinomas have been cited in association with human breast cancer. In this subset of cited genes, a significant enrichment of genes associated with cell cycle regulation, protein kinases, DNA integrity checkpoint and protein metabolism was observed. In addition, more than 60 percent of the van't Veer 70-gene prognostic signature was identified in our differential gene list. It is therefore likely that gene expression profiles may also predict metastasis in CMT. For instance, one of the dogs without lymph node metastases at the time of investigation developed later distant metastases. Strikingly, this case clustered within the group of metastatic tumors indicating that metastatic potential can be identified in CMT before clinically and morphologically overt metastases. It will be therefore of interest to compare the gene expression profiles of lymph node negative CMT with a comprehensive clinical follow-up to identify potential predictive patterns similar to human breast cancer.
Generally, metastatic CMT were characterized by increased expression of cell division genes. These findings are identical to a recent meta-analysis on expression profiles of metastasizing human breast cancer which found cell cycle and DNA replication pathways most significantly up-regulated [10, 13, 47]. Unexpectedly, this expression profile of increased proliferation was accompanied by significantly decreased growth factor receptor expression. Increased proliferation and malignant transformation of most cancer types is at least in part commonly attributed to aberrant growth factor signaling. In addition, growth factor signaling is also a major differentiation impulse. The decreased growth factor pathway expression in metastatic carcinomas hence may be a possible cause for the altered expression of cell differentiation and unlocked cell growth regulation genes in these CMT. The dysregulation of cell cycle control is therefore a dominant, cross species feature metastatic mammary tumors. Common "driver" mutations of this cell proliferation in the tumor genome of both species are unknown at the moment. The partly overlapping transcriptome of metastatic canine and human carcinomas, however, indicate similar mechanisms of mammary carcinogenesis and a suitable molecular model character of canine mammary tumors.
Genes associated with focal adhesion and regulation of actin cytoskeleton organization were significantly down-regulated in metastasizing CMT. This is in agreement with the observation that loss of cell to cell contact is a major step in the metastatic process and has been observed in canine mammary tumors before [48, 49]. The effect of the down-regulated pathway of actin regulation is somewhat unclear, but disorganization of the actin cytoskeleton is a common feature of malignant breast cancer cells . Of interest qPCR analyses confirmed the differences in mRNA expression levels in metastatic and non-metastatic carcinomas although the fold change in expression difference consistently lower in the qPCR assays when compared with microarray assays for all eight genes analyzed. The eight genes, AURKA, ALOX12, BMP6, ERBB4, HEPACAM2, IGFR2, RAD51, TGFBR3, were selected due to their description in the literature on canine mammary tumors before[17, 21–24, 51]. Their gene expression has been described metastatic carcinomas, metastases, adenomas and normal mammary gland and these studies showed a similar tendency of changes in gene expression between benign and malignant canine mammary tissue types[17, 21–24, 51]. Of note, a comparison of the proteome of the same subset of tumours did not identify similar proteins.
In conclusion, our findings show that metastatic spread of CMT to the lymph nodes is associated with a gene expression profile of increased cell cycle progression, altered cell differentiation and decreased growth factor signaling. Furthermore, several matrix modulators and cell adhesion genes have an altered expression in metastatic tumor cells while angiogenesis gene expression is not differentially regulated in metastatic and non-metastatic CMT. Several major aspects of metastasis-associated gene expression are therefore similar between human and canine mammary tumors. Furthermore, the differential expression profiles in metastatic and non-metastatic carcinomas allow a prediction of clinical behavior and offer several new hypotheses for molecular pathways and networks involved in the metastatic process of canine mammary tumors.
adipocyte-specific adhesion molecule
basic local alignment search tool
canine mammary tumour
cysteine-rich, angiogenic inducer 61
database for annotation, visualization and integrated discovery
epidermal growth factor
elastin microfibril interfacer 1
epidermal growth factor receptor 2
fibroblast growth factor
fibroblast growth factor receptor 1
hepatocyte cellular adhesion molecule
intercellular adhesion molecule
insulin-like growth factor
jumonji domain containing 6
hepatocyte growth factor
protocadherin gamma subfamily C3
RNA integrity number
tTransforming growth factor-beta
World Health Organization.
Benjamin SA, Lee AC, Saunders WJ: Classification and behavior of canine mammary epithelial neoplasms based on life-span observations in beagles. Vet Pathol. 1999, 36 (5): 423-436. 10.1354/vp.36-5-423.
Kitchell BE, Loar AS: Diseases of the mammary glands. Handbook of Small Animal Practice. Edited by: Morgan RV. 1997, WB Saunders, Philadelphia, 615-625. 3
Misdorp W: Tumors of the mammary gland. Tumours in Domestic Animals. Edited by: Meuten DJ. 2002, Iowa State Press, Iowa, 575-606. full_text. 4
Fitzgibbons PL, Page DL, Weaver D, Thor AD, Allred DC, Clark GM, Ruby SG, O'Malley F, Simpson JF, Connolly JL, et al: Prognostic factors in breast cancer. College of American Pathologists Consensus Statement 1999. Archives of pathology & laboratory medicine. 2000, 124 (7): 966-978.
Misdorp W, Hart AA: Prognostic factors in canine mammary cancer. J Natl Cancer Inst. 1976, 56 (4): 779-786.
Misdorp W, Hart AA: Canine mammary cancer. II. Therapy and causes of death. The Journal of small animal practice. 1979, 20 (7): 395-404. 10.1111/j.1748-5827.1979.tb06744.x.
Setoguchi A, Sakai T, Okuda M, Minehata K, Yazawa M, Ishizaka T, Watari T, Nishimura R, Sasaki N, Hasegawa A, et al: Aberrations of the p53 tumor suppressor gene in various tumors in dogs. Am J Vet Res. 2001, 62 (3): 433-439. 10.2460/ajvr.2001.62.433.
Nieto A, Perez-Alenza MD, Del Castillo N, Tabanera E, Castano M, Pena L: BRCA1 expression in canine mammary dysplasias and tumours: relationship with prognostic variables. J Comp Pathol. 2003, 128 (4): 260-268. 10.1053/jcpa.2002.0631.
Muthuswamy SK, Li D, Lelievre S, Bissell MJ, Brugge JS: ErbB2, but not ErbB1, reinitiates proliferation and induces luminal repopulation in epithelial acini. Nat Cell Biol. 2001, 3 (9): 785-792. 10.1038/ncb0901-785.
Thomassen M, Tan Q, Kruse TA: Gene expression meta-analysis identifies metastatic pathways and transcription factors in breast cancer. BMC Cancer. 2008, 8: 394-10.1186/1471-2407-8-394.
Lu X, Lu X, Wang ZC, Iglehart JD, Zhang X, Richardson AL: Predicting features of breast cancer with gene expression patterns. Breast Cancer Res Treat. 2008, 108 (2): 191-201. 10.1007/s10549-007-9596-6.
Weigelt B, Hu Z, He X, Livasy C, Carey LA, Ewend MG, Glas AM, Perou CM, Van't Veer LJ: Molecular portraits and 70-gene prognosis signature are preserved throughout the metastatic process of breast cancer. Cancer Res. 2005, 65 (20): 9155-9158. 10.1158/0008-5472.CAN-05-2553.
van 't Veer LJ, Dai H, van de Vijver MJ, He YD, Hart AA, Mao M, Peterse HL, van der Kooy K, Marton MJ, Witteveen AT, et al: Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002, 415 (6871): 530-536. 10.1038/415530a.
Rao NA, van Wolferen ME, Gracanin A, Bhatti SF, Krol M, Holstege FC, Mol JA: Gene expression profiles of progestin-induced canine mammary hyperplasia and spontaneous mammary tumors. J Physiol Pharmacol. 2009, 60 (Suppl 1): 73-84.
Uva P, Aurisicchio L, Watters J, Loboda A, Kulkarni A, Castle J, Palombo F, Viti V, Mesiti G, Zappulli V, et al: Comparative expression pathway analysis of human and canine mammary tumors. BMC genomics. 2009, 10: 135-10.1186/1471-2164-10-135.
Khanna C, Lindblad-Toh K, Vail D, London C, Bergman P, Barber L, Breen M, Kitchell B, McNeil E, Modiano JF, et al: The dog as a cancer model. Nat Biotechnol. 2006, 24 (9): 1065-1066. 10.1038/nbt0906-1065b.
Klopfleisch R, Klose P, Gruber AD: The combined expression pattern of BMP2, LTBP4, and DERL1 discriminates malignant from benign canine mammary tumors. Vet Pathol. 2010, 47 (3): 446-454. 10.1177/0300985810363904.
Sassi F, Benazzi C, Castellani G, Sarli G: Molecular-based tumour subtypes of canine mammary carcinomas assessed by immunohistochemistry. BMC Vet Res. 2010, 6: 5-10.1186/1746-6148-6-5.
Misdorp W, Else RW, Hellmén E, Lipscomb TP: Histological classification of mammary tumors of the dog and the cat. 1999, Washington: Armed Forces Institute of Pathology, 7:
Elston CW, Ellis IO: Pathological prognostic factors in breast cancer. I. The value of histological grade in breast cancer: experience from a large study with long-term follow-up. Histopathology. 1991, 19 (5): 403-410. 10.1111/j.1365-2559.1991.tb00229.x.
Klopfleisch R, Henning H, Costa da A, Klose P, Gruber AD: Insulin receptor is expressed in normal canine mammary gland and benign adenomas but decreased in metastatic canine mammary carcinomas similar to human breast cancer. Vet Comp Oncol. 2010
Klopfleisch R, Gruber AD: Increased Expression of BRCA2 and RAD51 in Lymph Node Metastases of Canine Mammary Adenocarcinomas. Vet Pathol. 2009, 46 (3): 416-422. 10.1354/vp.08-VP-0212-K-FL.
Klopfleisch R, Klose P, da Costa A, Brunnberg L, Gruber AD: HEPACAM1 and 2 are differentially regulated in canine mammary adenomas and carcinomas and its lymph node metastases. BMC Vet Res. 2010, 6: 15-10.1186/1746-6148-6-15.
Klopfleisch R, Schutze M, Gruber AD: Downregulation of transforming growth factor beta (TGFbeta) and latent TGFbeta binding protein (LTBP)-4 expression in late stage canine mammary tumours. The Veterinary Journal. 2009
Klopfleisch R, Gruber AD: Derlin-1 and stanniocalcin-1 are differentially regulated in metastasizing canine mammary adenocarcinomas. J Comp Pathol. 2009, 141 (2-3): 113-120. 10.1016/j.jcpa.2008.09.010.
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3 (7): RESEARCH0034-10.1186/gb-2002-3-7-research0034.
Etschmann B, Wilcken B, Stoevesand K, von der Schulenburg A, Sterner-Kock A: Selection of reference genes for quantitative real-time PCR analysis in canine mammary tumors using the GeNorm algorithm. Vet Pathol. 2006, 43 (6): 934-942. 10.1354/vp.43-6-934.
Cardin S, Libby E, Pelletier P, Le Bouter S, Shiroshita-Takeshita A, Le Meur N, Leger J, Demolombe S, Ponton A, Glass L, et al: Contrasting gene expression profiles in two canine models of atrial fibrillation. Circ Res. 2007, 100 (3): 425-433. 10.1161/01.RES.0000258428.09589.1a.
Huang da W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4 (1): 44-57. 10.1038/nprot.2008.211.
Hu Y, Hines LM, Weng H, Zuo D, Rivera M, Richardson A, LaBaer J: Analysis of genomic and proteomic data using advanced literature mining. J Proteome Res. 2003, 2 (4): 405-412. 10.1021/pr0340227.
Salomonis N, Hanspers K, Zambon AC, Vranizan K, Lawlor SC, Dahlquist KD, Doniger SW, Stuart J, Conklin BR, Pico AR: GenMAPP 2: new features and resources for pathway analysis. BMC bioinformatics. 2007, 8: 217-10.1186/1471-2105-8-217.
Klopfleisch R, Gruber AD: Differential expression of cell cycle regulators p21, p27 and p53 in metastasizing canine mammary adenocarcinomas versus normal mammary glands. Res Vet Sci. 2009, 87 (1): 91-96. 10.1016/j.rvsc.2008.12.010.
Klopfleisch R, Schutze M, Gruber AD: Loss of p27 expression in canine mammary tumors and their metastases. Res Vet Sci. 2010, 88 (2): 300-303. 10.1016/j.rvsc.2009.08.007.
Benson JR: Role of transforming growth factor beta in breast carcinogenesis. The lancet oncology. 2004, 5 (4): 229-239. 10.1016/S1470-2045(04)01426-3.
Schwertfeger KL: Fibroblast growth factors in development and cancer: insights from the mammary and prostate glands. Current drug targets. 2009, 10 (7): 632-644.
Muraoka-Cook RS, Feng SM, Strunk KE, Earp HS: ErbB4/HER4: role in mammary gland development, differentiation and growth inhibition. Journal of mammary gland biology and neoplasia. 2008, 13 (2): 235-246. 10.1007/s10911-008-9080-x.
Lah TT, Kalman E, Najjar D, Gorodetsky E, Brennan P, Somers R, Daskal I: Cells producing cathepsins D, B, and L in human breast carcinoma and their association with prognosis. Human pathology. 2000, 31 (2): 149-160. 10.1016/S0046-8177(00)80214-2.
Van den Eynden GG, Van Laere SJ, Van der Auwera I, Trinh XB, Van Marck EA, Dirix LY, Vermeulen PB: Two distinct expression patterns of urokinase, urokinase receptor and plasminogen activator inhibitor-1 in colon cancer liver metastases. Int J Cancer. 2009, 125 (6): 1494-1496. 10.1002/ijc.24524.
Nakopoulou L, Panayotopoulou EG, Giannopoulou I, Alexandrou P, Katsarou S, Athanassiadou P, Keramopoulos A: Stromelysin-3 protein expression in invasive breast cancer: relation to proliferation, cell survival and patients' outcome. Mod Pathol. 2002, 15 (11): 1154-1161. 10.1097/01.MP.0000037317.84782.CD.
Bon G, Di Carlo SE, Folgiero V, Avetrani P, Lazzari C, D'Orazi G, Brizzi MF, Sacchi A, Soddu S, Blandino G, et al: Negative regulation of beta4 integrin transcription by homeodomain-interacting protein kinase 2 and p53 impairs tumor progression. Cancer Res. 2009, 69 (14): 5978-5986. 10.1158/0008-5472.CAN-09-0244.
Moh MC, Zhang T, Lee LH, Shen S: Expression of hepaCAM is downregulated in cancers and induces senescence-like growth arrest via a p53/p21-dependent pathway in human breast cancer cells. Carcinogenesis. 2008, 29 (12): 2298-2305. 10.1093/carcin/bgn226.
Yoon KJ, Phelps DA, Bush RA, Remack JS, Billups CA, Khoury JD: ICAM-2 expression mediates a membrane-actin link, confers a nonmetastatic phenotype and reflects favorable tumor stage or histology in neuroblastoma. PloS one. 2008, 3 (11): e3629-10.1371/journal.pone.0003629.
Christgen M, Christgen H, Heil C, Krech T, Langer F, Kreipe H, Lehmann U: Expression of KAI1/CD82 in distant metastases from estrogen receptor-negative breast cancer. Cancer science. 2009, 100 (9): 1767-1771. 10.1111/j.1349-7006.2009.01231.x.
Hirschfeld M, zur Hausen A, Bettendorf H, Jager M, Stickeler E: Alternative splicing of Cyr61 is regulated by hypoxia and significantly changed in breast cancer. Cancer Res. 2009, 69 (5): 2082-2090. 10.1158/0008-5472.CAN-08-1997.
Gagnon ML, Bielenberg DR, Gechtman Z, Miao HQ, Takashima S, Soker S, Klagsbrun M: Identification of a natural soluble neuropilin-1 that binds vascular endothelial growth factor: In vivo expression and antitumor activity. Proc Natl Acad Sci USA. 2000, 97 (6): 2573-2578. 10.1073/pnas.040337597.
Wang Y, Klijn JG, Zhang Y, Sieuwerts AM, Look MP, Yang F, Talantov D, Timmermans M, Meijer-van Gelder ME, Yu J, et al: Gene-expression profiles to predict distant metastasis of lymph-node-negative primary breast cancer. Lancet. 2005, 365 (9460): 671-679.
Gama A, Paredes J, Gartner F, Alves A, Schmitt F: Expression of E-cadherin, P-cadherin and beta-catenin in canine malignant mammary tumours in relation to clinicopathological parameters, proliferation and survival. Vet J. 2008, 177 (1): 45-53. 10.1016/j.tvjl.2007.05.024.
De Matos AJ, Lopes CC, Faustino AM, Carvalheira JG, Rutteman GR, Gartner Mde F: E-cadherin, beta-catenin, invasion and lymph node metastases in canine malignant mammary tumours. Apmis. 2007, 115 (4): 327-334. 10.1111/j.1600-0463.2007.apm_544.x.
Otto AM, Muller CS, Huff T, Hannappel E: Chemotherapeutic drugs change actin skeleton organization and the expression of beta-thymosins in human breast cancer cells. Journal of cancer research and clinical oncology. 2002, 128 (5): 247-256. 10.1007/s00432-002-0332-7.
Klopfleisch R, Schutze M, Gruber AD: RAD51 protein expression is increased in canine mammary carcinomas. Vet Pathol. 2010, 47 (1): 98-101. 10.1177/0300985809353310.
Klopfleisch R, Klose P, Weise C, Bondzio A, Multhaup G, Einspanier R, Gruber AD: The Proteome of Metastatic Canine Mammary Carcinomas: Similarities and Differences to Human Breast Cancer. J Proteome Res. 2010
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/10/618/prepub
We thank L. Brunnberg, H. Linzmann, M. Köhle and T. Lottermoser for providing tumors samples and S. Mende and M. Schärig for excellent technical assistance. Research by RK and ADG is funded be the German Research Council (DFG KL 2240/1).
The authors declare that they have no competing interests.
All authors read and approved the manuscript. RK conducted the experimental work including microarray hybridization. RK and ADG designed the study, performed the histological analysis and the follow-up of the patients. DL, MH and RK accomplished the data analysis and prepared the manuscript.
Electronic supplementary material
Additional file 2: excel data sheet, gene list of differentially expressed genes in metastatic CMT, selection criteria included fold-change >2 or <-2 and q < 0.001(XLS 200 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Klopfleisch, R., Lenze, D., Hummel, M. et al. Metastatic canine mammary carcinomas can be identified by a gene expression profile that partly overlaps with human breast cancer profiles. BMC Cancer 10, 618 (2010). https://doi.org/10.1186/1471-2407-10-618