Abstract
Background/Aim: Neutrons are used as a type of high linear energy transfer (LET) radiation and they have stronger carcinogenic effects compared to low LET radiation. We sought to clarify the features of mammary carcinomas for which the incidence increases when these were exposed to neutron radiation. Materials and Methods: We compared mammary carcinomas from female Sprague-Dawley rats irradiated at 7 weeks of age with 0.485 Gy neutron beams or 0.5-Gy γ rays, with carcinomas of non-irradiated rats. Tumors were classified into luminal and non-luminal subtypes based on immunohistochemistry, while their copy number aberrations were determined using microarrays. Results: Neutrons and γ rays significantly increased the incidence of luminal carcinomas. The carcinomas in the three groups contained multiple aberrations affecting 46 genes for which mutations have been reported in human breast cancer. Conclusion: Neutrons and γ rays increase the incidence of luminal mammary carcinoma in rats, probably via genetic aberrations similar to those found in human breast cancer patients.
Exposure to radiation increases the risk of breast cancer, as revealed by epidemiological studies of atomic bomb survivors and patients who have undergone radiotherapy (1). Neutrons can be used as a type of high linear energy transfer (LET) radiation and produce more complex DNA damage with stronger carcinogenic effects compared to low LET radiation (2). In clinical settings, neutrons are generated as a secondary type of radiation during proton therapy (3), However, the carcinogenic effects of neutrons in humans have not yet been reported.
Cancer is caused by the accumulation of genomic and epigenomic alterations, such as gene mutations, chromosomal and DNA copy number aberrations in addition to DNA methylation (4, 5). Rats develop mammary carcinomas that are similar to human breast cancer with regard to their hormone dependence and pathology (6) and thus have been used as a model of radiation-induced breast cancer (7). Moreover, various animal experiments have revealed that neutrons have a strong effect on mammary carcinogenesis (8). Other studies have revealed various genomic and epigenomic alterations in radiation-induced cancers (9, 10). Although copy number aberrations are common in many types of tumors, there are only a few findings on genetic alterations of cancers induced by neutrons (11).
Breast cancer is a heterogeneous disease associated with genomic alterations and is classified into subtypes such as luminal, triple-negative, and human epidermal growth factor receptor 2 (HER2)-enriched (12). A few studies have investigated the subtypes of breast cancer that occur after X-ray exposure (13-15), however, the results are inconsistent. Furthermore, although other studies have assessed the incidence of mammary carcinoma after exposure to neutrons (2), none have described the subtypes of cancer advancing following radiation, while it is unclear which subtype(s) contributes to the increased risk of breast cancer caused by it.
We have previously conducted a series of experiments to quantify the incidence of neutron-induced mammary carcinomas in rats (16, 17). In our present study, we further analysed the mammary carcinomas from those experiments to understand the subtype(s) of these carcinomas for which the incidence was increased following neutron exposure. We also clarified the DNA copy number aberrations of these mammary carcinomas to elucidate the mechanisms of neutron-induced carcinogenesis.
Materials and Methods
Tissue samples. Tumor samples were obtained from our previous study (8), which was approved by the Institutional Animal Care and Use Committee of National Institute of Radiological Sciences. Briefly, in that previous study, 7-week-old female Jcl:SD rats (Clea Japan Inc., Tokyo, Japan) were whole-body irradiated with 9Be(d,nγ)10B fast neutrons (mean energy, 2 MeV) or 137Cs γ rays, and mammary carcinoma development was assessed by weekly palpation until the age of 90 weeks. All of the carcinomas available from the neutron-irradiated (0.485 Gy, n=16), γ-ray–irradiated (0.5 Gy, n=14), and the matched untreated rats (n=28), excluding the ones from rats found dead, were analysed.
Immunostaining and classification of tumor subtypes. Formalin-fixed, paraffin-embedded tissue sections were immunostained with antibodies against the estrogen receptor (ER) α (clone 6F11, Leica Biosystems, Nussloch, Germany; dilution, 1:400), the progesterone receptor (PgR; clone SP42, Acris Antibodies, Herford, Germany; 1:400), HER2 (a mixture of clones e2-4001 and 3B5, Thermo Fisher Scientific, Waltham, MA, USA; 1:100), and the proliferation marker Ki-67 (clone SP6, Spring Bioscience, Pleasanton, CA, USA; 1:200) as previously described (17). Each slide was scanned using a NanoZoomer XR C12000-03 slide scanner (Hamamatsu Photonics, Hamamatsu, Japan), and all images were saved in NDPI format. Images of 6 regions of interest (ROI) per slide were captured independently by two researchers (i.e., 12 ROIs for each carcinoma) at 40× magnification using the NDP.view2 software (Hamamatsu Photonics) and were stored in JPEG format. The percentage of antigen-positive epithelial tumor cells was determined using the Tissue Studio version 3.6.1 software (Definiens, Munich, Germany) for automated scoring and averaging over the 12 ROIs. For ERα and PgR, the cutoff for positivity of a tumor was set at 1% (12). For Ki-67, the cutoff was set at the median percentage (14.8%) (12). For HER2, tumors were classified as previously described (12). Luminal and non-luminal tumor subtypes were defined according to the standard classification criteria of human breast cancer (12), i.e., luminal subtype is defined as positive for either ERα or PgR and negative for HER2, and non-luminal subtype defined as any other.
Laser microdissection and DNA extraction. Neoplastic epithelial cells were collected using a laser microdissection system (MMI CellCut, Molecular Machines & Industries, Munich, Germany). Briefly, cryosections (20 μm thick) were prepared from frozen tissues embedded in an optimal cutting temperature compound (Sakura Finetek Japan, Tokyo, Japan) on membrane slides (MMI Membrane Slides, RNA free, Molecular Machines & Industries). Following fixation in 2-propanol for 45 seconds, sections were stained with hematoxylin and eosin. Dehydration was then achieved through a graded ethanol series and xylene for 15 seconds each. Neoplastic epithelial cells were subsequently dissected from the frozen sections and were collected. Since our preliminary examination indicated that 50 ng of double-stranded DNA was retrievable from 1 mm2 of sectioned epithelial cells,, a total of 5 mm2 or more of epithelial cells were dissected for the analyses described below. DNA was extracted from dissected epithelial cell sections as well as from normal liver sections stained in the same manner as the carcinoma sections, which served as the source of reference DNA from the individual rats. DNA extraction was performed using a QIAmp DNA Micro kit (Qiagen, Hilden, Germany) and isolated DNA was collected in 20 μL water (UltraPure DNase/RNase-Free Distilled Water, Thermo Fisher Scientific).
Microarray-based comparative genomic hybridization (aCGH). Double-stranded DNA (50 ng) was fragmented and amplified using the GenomePlex Whole Genome Amplification kit (Sigma-Aldrich, St. Louis, MO, USA). The amplified DNAs from normal liver and mammary carcinoma samples were labelled with cyanine 3– and cyanine 5–dUTP, respectively, and was purified using a purification column (Agilent Technologies, Santa Clara, CA, USA). Labelled DNA was hybridized with probes of microarrays (SurePrint G3 Rat CGH, 4×180K; Agilent Technologies) at 67°C with rotation at 20 rpm for 24 hours, and was then washed with Wash Buffers 1 and 2 (Agilent Technologies). The resolution of microarrays was ~17.5 kb (as an average) with 155,049 probes. Microarray scanning was performed using the Agilent G2505C microarray scanner. Fluorescence intensity values were obtained from the scanned images with the Agilent Feature Extraction software (ver. 10.7.3.1) and were analysed using the Agilent Genomic Workbench software 7.0.4.0.
Statistical analysis. Statistical analysis was performed using the statistical software R with the aid of a graphical user interface EZR (Jichi Medical University, Saitama, Japan) (18). Comparison among three groups was performed using the Kruskal-Wallis test. Comparison between two groups was performed using the Mann-Whitney's U-test. A p-value of <0.05 was considered statistically significant.
Results
Neutrons and γ rays increase the incidence of luminal carcinoma. We selected two groups of rats from our previous studies (8, 16): i) irradiated with neutrons (0.485 Gy, n=24) and ii) with γ rays (0.5 Gy, n=24) as well as a group of matched, non-irradiated rats (n=76). Table I summarizes the incidence, tumor number, hazard ratio (i.e., ratio of the probabilities of having a new carcinoma per unit time), and the age at which mammary carcinomas were first detected in these groups. Exposure to either neutrons or γ-rays significantly increased the incidence and the hazard ratio, whereas neutron exposure also significantly reduced the age at which the carcinoma was first detected, as compared to the non-irradiated group (Table I).
To reveal the subtypes of mammary carcinomas in these groups of rats, we performed immunohistochemical staining for ERα, PgR, HER2, and Ki-67 in all available tumors. As expected, signals for ERα, PgR, and Ki-67 were observed in the nucleus, whereas the HER2 signal was in the plasma membrane (Figure 1A-D). The percentage of cells that were positive for each antigen did not differ significantly between groups (Figure 1E). Next, we classified all available mammary carcinomas with respect to being luminal or non-luminal, based on the immunohistochemistry results. The fractions of luminal carcinomas in the neutron- and γ-ray–irradiated groups tended to be higher compared to the non-irradiated group (Table II). The hazard ratio of luminal mammary carcinomas, which were palpable, was also significantly higher in the neutron- and γ-ray–irradiated groups compared to the non-irradiated group (Table III and Figure 2A). Neutron exposure also increased the hazard ratio of palpable non-luminal carcinomas at a marginally significant level and significantly reduced the age at first palpation of luminal and non-luminal carcinomas compared to the non-irradiated group (Table III and Figure 2B), although the number of non-luminal carcinomas was relatively small in the neutron-irradiated group. Taken together, the increased incidence of mammary carcinomas in the neutron-irradiated rats mostly involved luminal carcinomas.
Rat mammary carcinomas include copy number changes of genes associated with human breast cancer. We next performed an aCGH analysis of 16 mammary carcinomas (i) 7 from the non-irradiated, ii) 4 from the neutron-irradiated, and iii) 5 from the γ-ray–irradiated groups) to detect DNA copy number aberrations. We identified 99 aberrations: i) 42 in the non-irradiated: ii) 27 in the neutron-irradiated, and iii) 30 in the γ-ray–irradiated groups, none of which displayed any significant differences between them (Figure 3). Given that ionizing radiation is known to cause DNA strand breaks (19), we then focused on partial-chromosome copy number aberrations, excluding whole-chromosome aberrations, because the latter is more likely related to abnormal chromosome separation during mitosis (20, 21). The results indicated no significant difference in the number of partial-chromosome aberrations among the non-irradiated, neutron-irradiated, and γ-ray–irradiated groups (Table IV).
We subsequently examined whether these aberrations included genes important for the genesis of breast cancer. For this analysis, we focused on 92 rat orthologues of driver genes of human breast cancer, which have been previously reported in a study using next-generation sequencing (22). Our analysis indicated that 46 of the 92 genes were affected due to 11 partial- and 22 whole-chromosome aberrations in 7 carcinomas (data not shown). We further classified these 46 genes into: i) 13 proto-oncogenes, ii) 17 tumor suppressor genes, and iii) 16 ambivalent (i.e., genes thath can be both) genes, based on information from the COSMIC database (http://cancer.sanger.ac.uk/cosmic). We excluded contradictory copy number changes, such as gain of tumor suppressors and loss of proto-oncogenes and we ultimately obtained 26 aberrations affecting 44 genes, out of which 10 were partial-chromosome aberrations (affected genes are listed in Table V). Of note, many of these genes (Erbb2, Erbb3, Pten, Nras, Notch1, Notch2, Smad4, and Apc) encode signal transducers, and their orthologues are mutated in human breast cancer. There were no differences in the distribution of these genes among our three groups (Table V).
These results revealed that the rat mammary carcinomas we analysed had multiple DNA copy number aberrations—either partial- or whole-chromosome – and these aberrations were not associated with tumor etiology (i.e., spontaneous, neutron-induced, or γ-ray–induced) even though they affected regulatory genes that undergo mutations in human breast cancer.
Discussion
In the present study, we analysed the subtype and copy number aberrations of rat mammary carcinomas to reveal a mechanism of carcinogenesis induced by neutrons and γ rays. We, herein, revealed that these exposures increased the incidence of luminal mammary carcinoma. Although a few studies have investigated the subtype of radiation-induced breast cancer in humans, there has been no consistent evidence that radiation preferentially induces specific subtypes (13-15). Notably, the subtype of breast cancer is influenced by various factors such as the number of births and obesity, and differences in these reproductive and lifestyle factors might underlie the inconsistency among the reported results (23). On the other hand, rats mainly develop hormone receptor–positive mammary carcinomas in both spontaneous and chemical-induced models (6, 17). Our present data concerning subtype frequency are consistent with published data (17).
Our aCGH analysis revealed copy number aberrations in mammary carcinomas of rats. The number of aberrations we detected here is greater than what we reported previously (24) possibly due to the increased number of cancer cells we isolated using laser microdissection. We found that these aberrations were present in a variety of genes that are frequently affected in human breast cancer (22), thus supporting the idea that these aberrations play important functional roles in the development of carcinomas in rats. On the other hand, we did not find any significant difference among the carcinomas in all our groups with regards to the frequency, the number, or the identity of genes for which copy number aberrations were apparent. This observation suggests that the molecular mechanism of mammary carcinogenesis induced by neutron or γ-ray irradiation is similar to that of spontaneous tumor development, at least with respect to copy number aberrations. This does not indicate, however, whether the type of radiation has any influence on carcinogenesis, as the age at which tumors were detected was earlier and the hazard ratio was greater for rats in the irradiated groups compared to those in ethe non-irradiated group. Thus, irradiation may have increased the probability of the occurrence of these copy number aberrations, or it has shortened the onset period. In the present study we did not analyze other types of mutations except for copy number aberrations, thus, further analysis is warranted, using more comprehensive techniques, such as next-generation sequencing.
Our results indicate that the increased incidence of rat mammary carcinoma caused by exposure to neutrons or γ rays could primarily account for an increased incidence of luminal carcinoma. Finally, these results also suggest that copy number aberrations affecting genes associated with human breast cancer are involved in the mechanism of tumorigenesis, providing information on previously unknown aspects of breast cancer development following exposure to neutrons.
Acknowledgements
This work was supported, in part, by a Grant-in-Aid for Young Scientists (B) JP26740022 and for Scientific Research (A) JP15H01834 from the Japan Society for the Promotion of Science KAKENHI. We thank all the technical staff at the Department of Radiation Effects Research, National Institutes of Radiological Sciences, and National Institute for Quantum and Radiological Science and Technology.
Footnotes
Authors' Contributions
HM, KD, TI, MN, YN, MT and TM designed the study and performed experiments. HM, KD, TI and MT analyzed the data and drafted the manuscript. MF, YS and KS supervised the work and helped draft and revise the manuscript. All authors read and approved the final manuscript.
This article is freely accessible online.
Conflicts of Interest
The Authors declare that there are no confricts of interest.
- Received December 26, 2018.
- Revision received January 25, 2019.
- Accepted February 1, 2019.
- Copyright© 2019, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved