- Research article
- Open Access
Estrogen receptors genotypes and canine mammary neoplasia
BMC Veterinary Researchvolume 15, Article number: 325 (2019)
Estrogens are essential for the development and proper function of several hormone-dependent organs. There are, however, several lines of evidence associating estrogens with mammary carcinogenesis. A marked individual genetic variability concerning estrogens biosynthesis, metabolism and mechanism of action was recognized and associated with human breast cancer susceptibility, clinical features and progression. Although some genetic variations in canine ESR1 gene were reported, their influence in clinicopathological features and progression of canine mammary tumors has not been fully evaluated. This study aims to assess the influence of SNPs in ESR1 gene (rs397512133, rs397510462, rs851327560, rs397510612, rs852887655, rs852684753 and rs852398698) in canine mammary tumors characteristics and progression. A group of 155 non-neutered bitches with mammary tumors was included in the study. Follow-up information was assessed 24 months after surgery.
Genetic profiles associated with a later onset of mammary tumors and less aggressive clinicopathological features, namely smaller tumor size (≤ 3 cm) with extensive tubular differentiation and low canine-adapted prognostic index (vet-NPI), were identified in this study.
Our data suggest that the ESR1 genetic profile may help on the decision regarding the selection of individual tailored preventive measures against canine mammary tumors development, such as early neutering.
Estrogens are crucial for normal development and function of several organs and systems, namely the mammary gland. In women, estrogens play a pivotal role in the development of the mammary branching ductal-alveolar system in puberty, throughout the menstrual cycle and also during pregnancy [1,2,3,4,5].
Due to their pro-proliferative and anti-apoptotic affects, estrogens have also been implicated in human breast cancer development and progression [6,7,8]. Furthermore, some of the intermediate compounds derived from estrogens metabolism have a well-recognized genotoxic action [7, 9, 10]. Indeed, several conditions related to increased or prolonged exposure of the mammary tissue to estrogens, such as early menarche, late menopause, oral contraception or hormone replacement therapy, constitute well-known risk factors for human breast cancer .
There is a large body of evidence linking estrogens to mammary carcinogenesis, also in canine species. Most mammary tumors are reported in females and the few cases documented in males are related to estrogen-secretor testicular neoplasms [12, 13]. Besides, canine contraceptive hormonal therapy has long been associated to an increased risk of mammary tumors . On the other hand, the protective effect of ovariectomy against the development of mammary neoplasia has been advocated for decades [15,16,17]. Moreover, levels of serum steroid hormones were reported to be higher in dogs with mammary carcinomas than in normal ones [18, 19].
Estrogens bind to estrogen receptors (ER) found in the normal canine mammary gland. Several studies reported changes in the expression pattern of ER in canine mammary gland in the course of neoplastic transformation and progression. An underexpression of ER has been documented in canine malignant mammary tumors, compared to benign neoplasms and to the normal mammary tissue [20,21,22,23]. This feature reinforces the importance of estrogens in canine mammary carcinogenesis. Furthermore, the decreased ER expression has been related to larger tumor size and lymph node metastasis, suggesting that ER status can be regarded as a marker with predictive and prognostic value in canine mammary tumors [21, 22].
Over the last decades, a considerable individual genetic variability concerning estrogens biosynthesis, metabolism and mechanism of action was recognized in humans. This individual genetic background is considered a significant contributor to breast cancer susceptibility, allowing the identification of subpopulations of women with higher breast cancer risk [24,25,26,27]. It has also been related to specific breast cancer clinical features, as well as to the clinical course of the disease [28,29,30]. In canine species data regarding the genetic profile related to estrogens and mammary tumor risk is not completely understood. Some ESR1 genetic differences were described between different dog breeds known to be at high and at low risk of mammary tumor development; furthermore, an association between ESR1 variation and the susceptibility to mammary tumors was described in a cohort of English Springer Spaniels [31, 32]. However, in a recent investigation, our group could not confirm a relationship between ESR1 genetic profile and the risk of development of mammary tumors . On the other hand, genetic variations in canine COMT gene (which encodes catechol-O-methyltransferase, an enzyme involved in estrogens metabolism through inactivation of carcinogenic catechol estrogens) has not been proved to influence susceptibility to canine mammary tumors. Nevertheless, COMT genetic variation has been linked to the age of onset of canine mammary carcinomas, to the development of high-grade mammary carcinomas, vascular invasion and recurrences [34,35,36]. To date, and to the best of the author’s knowledge, the influence of ESR1 genetic profile in clinicopathological features and progression of canine mammary tumors has not been fully assessed.
The aim of this study is to investigate the association between seven single nucleotide polymorphisms (SNPs) in canine ESR1 gene (rs397512133, rs397510462, rs851327560, rs397510612, rs852887655, rs852684753 and rs852398698) and clinicopathological features of canine mammary tumors and the clinical outcome of the disease.
One hundred and fifty-five non-neutered bitches were included in this study. The mean age of the whole population was 10.1 years-old (7–18 years-old): 9.8 years for dogs with benign tumors (n = 56; 36.1%) and 10.3 years for dogs with at least one malignant tumor (n = 99; 63.9%). Multiple tumours were diagnosed in 69.0% of the cases. Based on the Nottingham histological grading method 25/77 (32.5%) carcinomas were graded I, 36/77 (46.8%) were graded II and 16/77 (20.8%) were graded III. Vascular invasion and lymph node metastasis were observed in 14/96 (14.6%) and in 20/85 (23.5%) of cases, respectively. Two-year follow-up data was available for 88 dogs with malignant tumors. Of those, 44/88 (50%) were alive at the end of the follow-up period, while 21/88 (23.9%) died due to progression of the disease. Animals lost to follow-up (n = 16/88; 18.2%) and animals that died from causes not related to the mammary neoplasia (n = 7/88; 8%) were censored.
A significant association was found between the animal’s age at the time of the tumor diagnosis and SNPs rs397512133, rs397510462, rs851327560 and rs397510612. In fact, carriers of the variant allele for these SNPs developed mammary tumors later than wild type dogs (p = 0.014; p = 0.005; p = 0.007 and p = 0.008, respectively) (Fig. 1).
Table 1 presents results regarding the association between the studied SNPs and clinicopathological parameters, namely tumor number and size, histological classification (benign/malignant), mode of growth, NHG histological grade, NHG grading parameters (tubule formation, nuclear pleomorphism and mitotic index), NHG total score, vet-NPI, vascular invasion and lymph node metastases.
Carriers of variant allele for SNPs rs397512133, rs397510462, rs851327560 and rs397510612 developed smaller size carcinomas (≤ 3 cm) than wild type animals (p = 0.014; p = 0.020; p = 0.016; p = 0.020, respectively). The NHG grading parameters (tubule formation, nuclear pleomorphism and mitotic counts) were evaluated and a statistically significant relationship was observed between SNPs rs852887655, rs852684753 and rs852398698 and tubule formation (p = 0.039; p = 0.014; p = 0.026, respectively). The majority of carcinomas scored 3 for tubule formation were found in wild type animals. In fact, only 7.1% (rs852887655), 22.2% (rs852684753) and 25.0% (rs852398698) of the carcinomas scored 3 for this parameter corresponded to variant allele carriers. Furthermore, a statistically significant association was found between rs851327560 and the vet-NPI (p = 0.023). Most cases of vet-NPI > 4 (80.6%) corresponded to wild type animals, while only 19.4% of those cases were observed in variant allele carriers.
No significant associations could be established between any of the SNPs considered and the number of tumors, histological classification (benign/malignant tumors) and pattern of tumor growth, NGH histological grade, nuclear pleomorphism, mitotic counts, vascular invasion and lymph node metastasis. Furthermore, none of the SNPs considered were related to OS.
In this study genetic variations of the canine ESR1 were associated with the development of less aggressive canine mammary tumors.
Our results demonstrated that carriers of the variant allele for SNPs rs397512133, rs397510462, rs851327560 and rs397510612 tended to develop mammary neoplasia later in life than wild type dogs. This result finds parallel in data from a previous investigation of our group demonstrating that SNP in canine COMT gene were associated with the age of onset of mammary tumors . According to that study, variant allele carriers for SNP rs853046495 (also known as SNP COMT G482A) presented a threefold likelihood of developing mammary tumors after 9 years of age, when compared to wild type animals. Furthermore, a significantly longer waiting time of onset of malignant disease was observed in variant allele carriers than in wild type animals.
The presence of the variant allele for SNPs rs397512133, rs397510462, rs851327560 and rs397510612 was also associated with the development of small size malignant tumors (≤ 3 cm). Tumor size has long been considered an important prognostic factor in canine mammary neoplasia, with tumors larger than 3 cm being associated with short survival times [37,38,39,40].
Our data also demonstrated that variant allele genotypes for SNPs rs852887655, rs852684753 and rs852398698 were associated with carcinomas with a high percentage of tubular arrangements, a feature related to well differentiated and low-malignancy neoplasias, and usually associated with good prognosis .
Despite the fact that none of the SNPs studied was related to NHG histological grade, carriers of the variant allele for rs851327560 were significantly associated with low vet-NPI (≤4). This recently described index combines well-recognized prognostic factors for canine mammary tumors, namely tumor size, NHG histological grade and vascular/lymph node invasion. Previous studies from our group demonstrated that this canine-adapted index is associated with disease free-interval and OS in bitches [42, 43].
Taken together, our results allow the identification of a subgroup of dogs that tend to develop mammary tumors at older ages and with less aggressive clinicopathological features (small size tumors, with extensive tubular differentiation and low vet-NPI). Based on these findings, we hypothesize that carriers of the variant allele for ESR1 SNPs (rs397512133, rs397510462, rs851327560, rs397510612, rs852887655, rs852684753 and rs852398698) may possess receptors less sensitive to estrogen binding, resulting in a mammary tissue less responsive to the hormone, thus being more protected from its carcinogenic action. It is conceivable that in those dogs, a longer period of exposure to estrogens would be required to achieve carcinogenic levels, which could explain the development of tumors in older age. Despite the fact that SNPs assessed in this study are synonymous or intronic, and do not involve amino acid changes, they can induce alternative splicing, promote changes in RNA stability or structure, interfere with the speed and accuracy of transcription and/or translation and disturb protein folding [44,45,46,47]. Even if there is no change in protein production, as it happens in non-synonymous, the allegedly subtle characteristics produced by synonymous and intronic SNPs may explain the diversity of individual reaction to hormonal stimulus exhibited by different animals.
There are several epidemiological, biochemical and toxicological lines of evidence associating estrogens with mammary carcinogenesis [9, 10, 48, 49]. However, estrogens are also essential for the development and proper function of several hormone-dependent organs [50,51,52]. In fact, there are various detrimental effects related to neutering, such as urinary incontinence, musculoskeletal disorders and development of different types of neoplasia [53,54,55,56,57,58,59,60,61,62,63,64]. Besides, several sources of bias were identified in previous studies that link neutering to a decreased risk of mammary neoplasia. A more recent revision of the data previously available was performed, demonstrating that the scientific support for this evidence is weak . In this vein, a growing debate on whether and when spaying should be recommended has emerged. Several researchers have questioned the value of early spaying, taking into account the secondary effects of this procedure in the animal’s health and quality of life in the medium and long term [54, 60, 62]. In fact, the relationship between gonadectomy, cancer development, overall lifespan and cancer-related death is complex. Besides, biological and molecular mechanisms involved in this process are poorly understood. The conflicting data obtained from different studies may be related to the individual genetic profile of the animal, namely to SNPs.
It is important to identify subsets of animals with increased risk or more susceptible to the development of aggressive mammary tumors that, in spite of the medium-long term secondary effects resulting from neutering, could benefit from that intervention. Similarly, recognition of the animals in which the deleterious consequences from neutering overlap their protective effect against the development of mammary neoplasia is of uttermost importance, avoiding unnecessary surgical interventions. In this sense, further research must be conducted considering other genetic variants and involving a larger number of animals, in order to assess the importance of the individual genetic profile in gonadectomy effects.
Our data suggest that canine ESR1 genetic profile may constitute a rational basis for the evaluation of the cost-benefit ratio related to early spaying. However, additional broad investigation is needed to expand these findings and to clarify the relevance of the genetic profile, assisting in the selection of the animals that could benefit from neutering as an individualized preventive strategy against the development of canine mammary neoplasia.
The study was conducted involving 155 non neutered bitches with histologically confirmed mammary tumors collected from the Veterinary Pathology Laboratory (ICBAS-University of Porto). All animals were treated with radical unilateral and/or partial mastectomy. Owners provided consent for surgery with curative intents as well as for the use of the material for research purposes. This protocol was approved by the Ethics Committee of the Institute of Biomedical Sciences Abel Salazar, University of Porto (P151/2016).
After surgery, mammary specimens were immediately fixed in a 10% buffered formalin solution and routinely processed for histopathological examination. For each case, clinicopathological features including age at the time of the diagnosis, tumor number and size (corresponding to the largest diameter measured by a pathologist (ACS) during trimming) were recorded. The histological diagnosis was established by consensus of 3 pathologists (ACS, MS and PDP) in a multi-head microscope, using the criteria of the World Health Organization for the classification of mammary tumors of dogs and cats . Each malignant tumor was assessed for the mode of growth and the presence of vascular invasion and regional lymph node metastasis, as previously described . In the subgroup of animals with multiple malignant tumors, a reference lesion was assigned for the statistical study, according to criteria previously reported . Histological grading was performed based on the Nottingham histological grading method -NHG . A veterinary adaptation of the human Nottingham Prognostic Index (vet-NPI) was also computed, as previously reported .
Follow-up data was obtained by consulting the medical records and by contact with the referring veterinarian. Disease-specific overall survival (OS) was calculated from the time of diagnosis to the date of the animal’s death/euthanasia due to the neoplastic disease. Animals that died or were euthanized for unrelated causes and those that were lost to follow-up were censored, respectively, at the time of death and at the data of their last clinical examination. Euthanasia was performed only in terminal stage of the disease. Necropsy examination was performed upon the owner consent.
Genomic DNA was extracted from peripheral blood samples (obtained by standard venipuncture) using High Pure PCR Template preparation kit (Roche). The DNA quality was evaluated by measuring the optical density and the quantity was assessed employing the NanoDrop 1000 Spectrophotometer. SNP genotyping was performed using MassARRAY iPLEX Gold Technology at the Unidade de Genómica/Serviço de Genotipagem do Instituto Gulbenkian de Ciência. This technology for SNP genotyping consists of an initial PCR reaction, followed by multiplexed primer extension (using mass-modified dideoxynucleotide terminators of an oligonucleotide primer) which anneals immediately upstream of the polymorphic site of interest. Additionally, the MALDI-TOF (matrix assisted laser desorption/ionization - time of flight) mass spectrometry allows the recognition of the SNP allele by the different mass of the extended primer [67, 68].
Seven canine ESR1 (chromosome 1) SNPs were assessed as referred in Table 2 and Fig. 2. Statistical analysis of data was performed using the computer software SPSS for Windows (version 25). Chi-square analysis (or Fisher’s exact test, when appropriated) was used to evaluate the significance of the relationship between ESR1 SNPs and the categorical variables. Cumulative risk curves were computed using Kaplan-Meier product-limit estimates method, with Log-rank (Mantel-Cox) tests being used to estimate the differences in risk of canine mammary tumors development according to each genetic profile. A 5% level was considered to define statistical significance.
Availability of data and materials
The datasets used and analyzed for this study are available from the corresponding author on reasonable request.
Estrogen receptor 1
Matrix assisted laser desorption/ionization - time of flight
Nottingham histologic grade
Single nucleotide polymorphism
Veterinary adapted Nottingham Prognostic Index
Hilton N, Clarke L, Graham D. Estrogen and progesterone signalling in the normal breast and its implications for cancer development. Mol Cell Endocrinol. 2018;466:2–14.
Macias H, Hinck L. Mammary gland development. Wiley Interdiscip Rev Dev Biol. 2012;1(4):533–57.
Arendt LM, Kuperwasser C. Form and function: how estrogen and progesterone regulate the mammary epithelial hierarchy. J Mammary Gland Biol Neoplasia. 2015;20(1–2):9–25.
Stingl J. Estrogen and progesterone in normal mammary gland development and in cancer. Horm Cancer. 2011;2(2):85–90.
Sternlicht MD, Kouros-Mehr H, Lu P, Werb Z. Hormonal and local control of mammary branching morphogenesis. Differentiation. 2006;74(7):365–81.
Parl FF, Dawling S, Roodi N, Crooke PS. Estrogen metabolism and breast cancer: a risk model. Ann N Y Acad Sci. 2009;1155:68–75.
Santen R, Cavalieri E, Rogan E, Russo J, Guttenplan J, Ingle J, Yue W. Estrogen mediation of breast tumor formation involves estrogen receptor-dependent, as well as independent, genotoxic effects. Ann N Y Acad Sci. 2009;1155:132–40.
Chen F-P, Chien M-H, Chen H-Y, Huang T-S, Ng Y-T. Effects of estradiol and progestogens on human breast cells: regulation of sex steroid receptors. Taiwanese J Obstet Gynecol. 2013;52(3):365–73.
Cavalieri E, Frenkel K, Liehr JG, Rogan E, Roy D. Chapter 4: estrogens as endogenous genotoxic agents—DNA adducts and mutations. JNCI Monographs. 2000;2000(27):75–94.
Liehr JG. Is estradiol a genotoxic mutagenic carcinogen? Endocr Rev. 2000;21(1):40–54.
Dall GV, Britt KL. Estrogen effects on the mammary gland in early and late life and breast Cancer risk. Front Oncol. 2017;7:110.
Saba CF, Rogers KS, Newman SJ, Mauldin GE, Vail DM: Mammary gland tumors in male dogs. J Vet Intern Med. 2007;21(5):1056-9.
Kwon S-C, Yoo D-Y, Ko M, Lee K-Y, Kwak H-H, Park I-C, Hwang I-K, Choi J-H, Chung J-Y. Mammary gland tumors in a male cocker spaniel. Acta Vet Scand. 2017;59:20.
Rutteman GR. Contraceptive steroids and the mammary gland: is there a hazard?--insights from animal studies. Breast Cancer Res Treat. 1992;23(1–2):29–41.
Schneider R, Dorn CR, Taylor DO. Factors influencing canine mammary cancer development and postsurgical survival. J Natl Cancer Inst. 1969;43(6):1249–61.
Sorenmo KU, Shofer FS, Goldschmidt MH. Effect of spaying and timing of spaying on survival of dogs with mammary carcinoma. J Vet Intern Med. 2000;14(3):266–70.
Kristiansen VM, Pena L, Diez Cordova L, Illera JC, Skjerve E, Breen AM, Cofone MA, Langeland M, Teige J, Goldschmidt M, et al. Effect of ovariohysterectomy at the time of tumor removal in dogs with mammary carcinomas: a randomized controlled trial. J Vet Intern Med. 2016;30(1):230–41.
Queiroga FL, Perez-Alenza MD, Silvan G, Pena L, Lopes C, Illera JC. Role of steroid hormones and prolactin in canine mammary cancer. J Steroid Biochem Mol Biol. 2005;94(1–3):181–7.
Illera JC, Perez-Alenza MD, Nieto A, Jimenez MA, Silvan G, Dunner S, Pena L. Steroids and receptors in canine mammary cancer. Steroids. 2006;71(7):541–8.
Rutteman GR, Misdorp W, Blankenstein MA, van den Brom WE. Oestrogen (ER) and progestin receptors (PR) in mammary tissue of the female dog: different receptor profile in non-malignant and malignant states. Br J Cancer. 1988;58(5):594–9.
Nieto A, Pena L, Perez-Alenza MD, Sanchez MA, Flores JM, Castano M. Immunohistologic detection of estrogen receptor alpha in canine mammary tumors: clinical and pathologic associations and prognostic significance. Vet Pathol. 2000;37(3):239–47.
de Las Mulas JM, Millan Y, Dios R. A prospective analysis of immunohistochemically determined estrogen receptor alpha and progesterone receptor expression and host and tumor factors as predictors of disease-free period in mammary tumors of the dog. Vet Pathol. 2005;42(2):200–12.
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 (Athens, Greece). 2011;25(3):455–65.
Martin A-M, Weber BL. Genetic and hormonal risk factors in breast Cancer. JNCI: J Natl Cancer Inst. 2000;92(14):1126–35.
Bugano DD, Conforti-Froes N, Yamaguchi NH, Baracat EC. Genetic polymorphisms, the metabolism of estrogens and breast cancer: a review. Eur J Gynaecol Oncol. 2008;29(4):313–20.
Evans DG, Howell SJ, Howell A. Personalized prevention in high risk individuals: Managing hormones and beyond. Breast (Edinburgh, Scotland). 2018;39:139–47.
Sun MY, Du HY, Zhu AN, Liang HY, de Garibay GR, Li FX, Li M, Yang XX. Genetic polymorphisms in estrogen-related genes and the risk of breast cancer among Han Chinese women. Int J Mol Sci. 2015;16(2):4121–35.
Larsen MJ, Thomassen M, Gerdes A-M, Kruse TA. Hereditary breast Cancer: clinical, pathological and molecular characteristics. Breast Cancer : Basic and Clinical Research. 2014;8:145–55.
Veronesi A, Cd G, Magri MD, Lombardi D, Zanetti M, Scuderi C, Dolcetti R, Viel A, Crivellari D, Bidoli E, et al. Familial breast cancer: characteristics and outcome of BRCA 1–2 positive and negative cases. BMC Cancer. 2005;5:70.
van’t Veer LJ, Dai H, van de Vijver MJ, He YD, AAM H, 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:530.
Borge KS, Borresen-Dale AL, Lingaas F. Identification of genetic variation in 11 candidate genes of canine mammary tumour. Vet Comp Oncol. 2011;9(4):241–50.
Borge KS, Melin M, Rivera P, Thoresen SI, Webster MT, von Euler H, Lindblad-Toh K, Lingaas F. The ESR1 gene is associated with risk for canine mammary tumours. BMC Vet Res. 2013;9:69.
Canadas A, Santos M, Nogueira A, Assis J, Gomes M, Lemos C, Medeiros R, Dias-Pereira P. Canine mammary tumor risk is associated with polymorphisms in RAD51 and STK11 genes. J Vet Diagn Invest. 2018;30(5):733–8.
Canadas A, Santos M, Pinto R, Medeiros R, Dias-Pereira P. Catechol-o-methyltransferase genotypes are associated with progression and biological behaviour of canine mammary tumours. Vet Comp Oncol. 2018;16(4):664–9.
Dias Pereira P, Lopes CC, Matos AJ, Pinto D, Gartner F, Lopes C, Medeiros R. Influence of catechol-O-methyltransferase (COMT) genotypes on the prognosis of canine mammary tumors. Vet Pathol. 2009;46(6):1270–4.
Dias Pereira P, Lopes CC, Matos AJ, Pinto D, Gartner F, Lopes C, Medeiros R. Estrogens metabolism associated with polymorphisms: influence of COMT G482a genotype on age at onset of canine mammary tumors. Vet Pathol. 2008;45(2):124–30.
Sorenmo KU, Kristiansen VM, Cofone MA, Shofer FS, Breen AM, Langeland M, Mongil CM, Grondahl AM, Teige J, Goldschmidt MH. Canine mammary gland tumours; a histological continuum from benign to malignant; clinical and histopathological evidence. Vet Comp Oncol. 2009;7(3):162–72.
Yamagami T, Kobayashi T, Takahashi K, Sugiyama M. Prognosis for canine malignant mammary tumors based on TNM and histologic classification. J Vet Med Sci. 1996;58(11):1079–83.
Perez Alenza MD, Pena L, Nieto AI, Castano M. Clinical and pathological prognostic factors in canine mammary tumors. Ann Ist Super Sanita. 1997;33(4):581–5.
Ferreira E, Bertagnolli AC, Cavalcanti MF, Schmitt FC, Cassali GD. The relationship between tumour size and expression of prognostic markers in benign and malignant canine mammary tumours. Vet Comp Oncol. 2009;7(4):230–5.
Fisher ER, Sass R, Fisher B. Pathologic findings from the national surgical adjuvant breast project. Correlations with concordant and discordant estrogen and progesterone receptors. Cancer. 1987;59(9):1554–9.
Santos M, Correia-Gomes C, Marcos R, Santos A, De Matos A, Lopes C, Dias-Pereira P. Value of the Nottingham histological grading parameters and Nottingham prognostic index in canine mammary carcinoma. Anticancer Res. 2015;35(7):4219–27.
Canadas A, Franca M, Pereira C, Vilaca R, Vilhena H, Tinoco F, Silva MJ, Ribeiro J, Medeiros R, Oliveira P, Dias-Pereira P, Santos M. Canine mammary tumors: comparison of classification and grading methods in a survival study. Vet Pathol. 2019;56:208-19.
Cooper DN. Functional intronic polymorphisms: buried treasure awaiting discovery within our genes. Hum Genomics. 2010;4(5):284–8.
Ramirez-Bello J, Jimenez-Morales M. Functional implications of single nucleotide polymorphisms (SNPs) in protein-coding and non-coding RNA genes in multifactorial diseases. Gac Med Mex. 2017;153(2):238–50.
Hunt R, Sauna ZE, Ambudkar SV, Gottesman MM, Kimchi-Sarfaty C. Silent (synonymous) SNPs: should we care about them? Methods Mol Biol (Clifton, NJ). 2009;578:23–39.
Fernandez-Calero T, Cabrera-Cabrera F, Ehrlich R, Marin M. Silent polymorphisms: can the tRNA population explain changes in protein properties? Life (Basel). 2016;6(1).
Yager JD, Liehr JG. Molecular mechanisms of estrogen carcinogenesis. Annu Rev Pharmacol Toxicol. 1996;36:203–32.
Gail MH. Twenty-five years of breast cancer risk models and their applications. J Natl Cancer Inst. 2015;107(5). https://doi.org/10.1093/jnci/djv042.
Hughes ZA, Liu F, Marquis K, Muniz L, Pangalos MN, Ring RH, Whiteside GT, Brandon NJ. Estrogen receptor neurobiology and its potential for translation into broad spectrum therapeutics for CNS disorders. Curr Mol Pharmacol. 2009;2(3):215–36.
Xiao J, Wang NL, Sun B, Cai GP. Estrogen receptor mediates the effects of pseudoprotodiocsin on adipogenesis in 3T3-L1 cells. Am J Physiol Cell Physiol. 2010;299(1):C128–38.
Swedenborg E, Power KA, Cai W, Pongratz I, Ruegg J. Regulation of estrogen receptor beta activity and implications in health and disease. Cell Mol Life Sci : CMLS. 2009;66(24):3873–94.
Arlt S, Wehrend A, Reichler IM. Neutering of female dogs - old and new insights into pros and cons. Tierarztl Prax Ausg K Kleintiere Heimtiere. 2017;45(4):253–63.
Hart BL, Hart LA, Thigpen AP, Willits NH. Long-term health effects of neutering dogs: comparison of Labrador retrievers with Golden retrievers. PLoS One. 2014;9(7):e102241.
Prymak C, McKee LJ, Goldschmidt MH, Glickman LT. Epidemiologic, clinical, pathologic, and prognostic characteristics of splenic hemangiosarcoma and splenic hematoma in dogs: 217 cases (1985). J Am Vet Med Assoc. 1988;193(6):706–12.
Ru G, Terracini B, Glickman LT. Host related risk factors for canine osteosarcoma. Vet J (London, England : 1997). 1998;156(1):31–9.
Veronesi MC, Rota A, Battocchio M, Faustini M, Mollo A. Spaying-related urinary incontinence and oestrogen therapy in the bitch. Acta Vet Hung. 2009;57(1):171–82.
Villamil JA, Henry CJ, Hahn AW, Bryan JN, Tyler JW, Caldwell CW. Hormonal and sex impact on the epidemiology of canine lymphoma. J Cancer Epidemiol. 2009;2009:591753.
Ware WA, Hopper DL. Cardiac tumors in dogs: 1982-1995. J Vet Intern Med. 1999;13(2):95–103.
Beauvais W, Cardwell JM, Brodbelt DC. The effect of neutering on the risk of urinary incontinence in bitches - a systematic review. J Small Anim Pract. 2012;53(4):198–204.
Beauvais W, Cardwell JM, Brodbelt DC. The effect of neutering on the risk of mammary tumours in dogs-a systematic review. J Small Anim Pract. 2012;53(6):314–22.
Torres de la Riva G, Hart BL, Farver TB, Oberbauer AM, Messam LL, Willits N, Hart LA. Neutering dogs: effects on joint disorders and cancers in golden retrievers. PloS one. 2013;8(2):e55937.
Hoffman JM, Creevy KE, Promislow DE. Reproductive capability is associated with lifespan and cause of death in companion dogs. PLoS One. 2013;8(4):e61082.
Kent MS, Burton JH, Dank G, Bannasch DL, Rebhun RB. Association of cancer-related mortality, age and gonadectomy in golden retriever dogs at a veterinary academic center (1989-2016). PLoS One. 2018;13(2):e0192578.
Misdorp W, Else RW, Hellmen E, Lipscomb TP. World Health Organization international histological classification of tumors of domestic animals vol. VII. 2nd ed. Washington, D.C: Armed Forces Institute of Pathology; 1999.
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–10.
Oeth P, del Mistro G, Marnellos G, Shi T, van den Boom D. Qualitative and quantitative genotyping using single base primer extension coupled with matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MassARRAY). Methods Mol Biol (Clifton, NJ). 2009;578:307–43.
Gabriel S, Ziaugra L, Tabbaa D. SNP genotyping using the Sequenom MassARRAY iPLEX platform. Curr Protoc Human Genet. 2009;2(2):12.
The authors would like to express their gratitude to Ricardo Pinto for his contribution in the gene map construction. We also thank to Veterinary Pathology Laboratory and Immunogenetic Laboratory (ICBAS, University of Porto) for technical support; to the clinicians Miguel França, Cristina Pereira, Jorge Ribeiro, Flora Tinoco, Raquel Vilaça, Maria João Silva and Hugo Vilhena; to Helena Frias of Vila Nova de Gaia Kennel and Aanifeira animal association; to the animals and owners involved in the study.
No funding was used in this study.
Ethics approval and consent to participate
This protocol was approved by the Ethics Committee of the Institute of Biomedical Sciences Abel Salazar, University of Porto (P151/2016).
Informed consent was obtained from all owners of animals that participate in this study. Euthanasia, required by the owners or proposed by the veterinary physician, was performed only in terminal stage of the disease using pentobarbital.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.