Transcriptomics and Cell Transformation Assay: an Integrated Approach to Evaluate the Effects of Low Dose Ionizing Radiation
Stefania Perdichizzi1, Maria Grazia Mascolo1, Francesca Rotondo1, Monica Vaccari1, Sandro Grilli2 and Annamaria Colacci1*
1Agency for Prevention, Environment and Energy, Italy
2Department of Experimental Pathology-Cancer Research Section, School of Medicine University of Bologna, Italy
*Corresponding author: Annamaria Colacci, Adjunct Professor of Environmental Hygiene, Head Center for Environmental Health and Prevention, Agency for Prevention, Environment and Energy (Arpae) Emilia-Romagna, University of Bologna, Italy
Article History
Received: June 28, 2021 Accepted: July 19, 2021 Published: July 23, 2021
Citation: Perdichizzi S, Mascolo MG, Rotondo F, et al. Transcriptomics and Cell Transformation Assay: an Integrated Approach to Evaluate the Effects of Low Dose Ionizing Radiation. Int J Onco Radiother. 2021;2(1):12‒19. DOI: 10.51626/ijor.2021.02.00005
Abstract
Background and aims: Ionizing radiation (IR) are a well-known carcinogenic agent, acting through genotoxic mechanisms. In the last years, great attention has been paid to the effects of IR at low doses and to the non-monotonic dose-response curve for IR exposures. To improve the knowledge of IR-mediated effects and possibly identify biomarkers for IR effects, we combined the Cell Transformation Assay (CTA) with transcriptomics, to correlate cytotoxicity and transformation endpoints with the modulation of gene profiles after IR exposure.
Methods: BALB/c3T3 cells were exposed to ionizing radiation ranging from 0.25Gy and 6Gy. Irradiated cells were seeded for the CTA 20h later. At the same time, RNA was extracted for microarray experiments. The cell clonal survival was significantly increased in 0.25Gy IR exposed cells, while the 3Gy dose strongly inhibited cellular growth. Cell transformation was observed only at the highest dose (3Gy).
Results: Cell’s transformation was observed at 1.5, 2 and 3Gy doses. The 0.25Gy dose, which was able to induce an increment of clonal efficiency, did not induce cell transformation. The gene expression profile, which was obtained by comparing cells treated with the highest tested dose of 3Gy with the cells exposed to the lowest, not transforming, dose of 0.25Gy, identified several genes related to mitotic cell cycle and cholesterol biosynthesis.
Conclusion: Our study showed that the up-regulation of genes belonging to the Spindle Assembly Checkpoint and mitosis progression could support the transforming ability of the 3Gy BALB/c3T3 exposed cells, probably through the involvement of genomic instability. Gene transcripts involved into cholesterol biosynthesis appear to be critical, as well. All these transcripts may be regarded as potential biomarkers of IR effects.
Keywords: Ionizing radiation; BALB/c3T3 transformation assay; Transcriptomics; Integrated testing approach; Gene modulation; Genome instability; Cholesterol biosynthesis; Mitochondrial dysfunction
Abbreviations
ACE: Absolute Clonal Efficiency; APC/C-Cdc20: Anaphase Promoting Complex/Cyclosome; AR: Adaptive Response; Aurka: Aurora Kinase A; Calm-1: Calmodulin 1; Cenpa: Centromere Protein A; Cdkn1a; p21: Cyclin-Dependent Kinase Inhibitor 1A; CIN: Chromosome Instability; CTA: Cell Transformation Assay; CVD: Cardiovascular Disease; EASE: Expression Analysis Systematic Explorer analysis; FDPS: Farnesyl Diphosphate Synthase; Gadd45a: DNA-Damage-Inducible 45 Alpha; Gapdh: Glyceraldehyde-3-Phosphate Dehydrogenase; GO: Gene Ontology; Hmgcr: HMG-CoA Reductase; HRS: Hyper-Radiosensitivity; IR: Ionizing Radiation; IRR: Radioresistance; PBS: Phosphate-Buffered Saline; Pmaip1; Noxa: Phorbol-12-Myristate-13-Acetate-Induced Protein 1; RCE: Relative Clonal Efficiency; SAC: Spindle Assembly Checkpoint; TF: Transformation Frequency
Introduction
Ionizing radiation (IR) represents an important source of both occupational and general population exposure. General population can be exposed through the environment, due to naturally occurring radioactive materials. Additional environmental exposure to IR can caused by the increased utilization of industrial processes requiring radiation sources, or as the consequence of disasters, involving nuclear plants. Medical exposure from diagnostic procedure, using x-rays or nuclear tracers, represents a large source of exposure, especially in developed countries with easy access to high levels of medical care [1,3]. Recent literature reports have risen concern about possible adverse effects related to the increased use of IR for diagnostic and therapeutic procedures [4]. As X- rays and -radiation are known to be human carcinogens, IR potential effects on humans have been of concern for over half a century. Several studies on IR-induced damage pointed out at a notable variety of genetic lesions by IR, including DNA double strand breaks, point mutations, chromosomal aberrations, and deletions [5-7]. Cellular oncogenes were shown to be targets for the effects of IR during in vivo carcinogenesis and in vitro cell transformation [8, 9]. However, the molecular mechanisms sustaining the response to IR, especially at low doses, are not fully understood as yet [10, 11].
It is well known that -radiations are able to alter the expression of stress-responsive genes at doses as low as 10 cGy and to induce large-scale changes in gene expression by activating several unexpected signaling pathways [12,15]. Given the complexity of the radiation response, it is evident that conventional approaches pointing at limited number of genes may not be satisfactory. Microarray technology allows the simultaneous analysis of thousands genes in a single experiment, providing new opportunities to understand the molecular events involved in the response to stressors. Modulated genes have been considered as potential biomarkers of IR exposure [16,18]. Studies in radio-sensitive mouse organs or cell lines were performed to identify differentially expressed genes related to IR exposure [19,23]. Even if these studies paved the way for the discovery of new potential IR biomarkers in human peripheral blood cells, the high variability in time or doses of exposure in both experimental and observational studies as well as the heterogeneity of the study design and final outcomes limit the applicability of these results [16,18].
Therefore, the possibility to identify markers of IR exposure or effect is still a challenge [24,25]. In recent years, the scientific community encouraged the use of experimental models based on the integration of in vitro methodologies with toxicogenomics technologies. This combined approach may provide information at the cellular and molecular level and promote a shift from in vivo costly and time-consuming animal studies to short term in vitro assays [26].
The cell transformation assay (CTA) can represent an alternative to animal models, when used in an integrated approach to support the weight of evidence in carcinogenicity evaluation [27, 28]. The cellular processes involved in in vitro cell transformation closely resembles in vivo carcinogenesis and occur as a result of comprehensive cellular responses to direct and indirect damage to DNA [29]. In particular, the in vitro transformation assay in BALB/c3T3 cells a system that benefits from a convenient protocol and high predictability of mammalian carcinogenicity [30-33]. The test is based on the malignant transformation of BALB/c3T3 embryonic mouse fibroblasts, which results in the induction of morphologically aberrant foci, formed by cells that lost contact inhibition and piled up with random orientation [31]. However, this assay alone does not provide information on the molecular mechanisms sustaining the carcinogenesis process, limiting its use in the regulatory context.
We propose here an approach that conjugates toxicogenomic technologies and in vitro BALB/c3T3 CTA to identify gene expression profiles supporting IR-induced transformation and possibly discriminate the response to low and high IR doses at the molecular level.
Whilst the integrated experimental approaches concerning the implementation of BALB/c3T3 transformation assay with toxicogenomic has been adopted to examine the effects of chemicals acting through genotoxic or non-genotoxic mechanisms in the carcinogenesis process, no similar studies are available for IR [34,35].
Our aim was to identify specific biological signatures of IR exposure in the BALB/c3T3 model by investigating modulated genes in the early steps of the oncotransformation. To address this issue, the effects of two IR doses, 0.25 and 3Gy, were investigated. Results show a different response according to the dose of exposure. Toxic and transforming effects were observed only in 3Gy exposed cells Differentially expressed genes between 3 and 0.25Gy treated cells are mainly related to the cell cycle network. These results show that the use of integrated testing strategies, including transcription analysis, provide further information on the response associated with different doses of ionizing radiation. The gene list we found would be of possible use in the monitoring of IR exposure.
Materials and Methods
Cells
The original stock of BALB/c3T3 cells, clone A31, was obtained from the American Type Culture Collection, Rockville, MD, USA. Working cultures were expanded from the original cryopreserved stock. Cells were grown in Dulbecco’s modified Eagle’s medium (D-MEM, Gibco BRL, Glasgow, Scotland) supplemented with 10% New-born Calf Serum (NCS, Gibco BRL). Only sub-confluent cells were used and the target cells were not maintained beyond the third passage after thawing. Cultures were maintained in a humidified incubator with an atmosphere of 5% CO2 in air at 37°C.
Cell irradiation
Exponentially growing BALB/c3T3 cells, were exposed to 137Cs radiation source (IBL 437C, 0.66 MeV, Dose-rate 221c Gy/min), to obtain an exposure range from 0.25 to 6Gy, then maintained in culture for 20 hours before seeding.
Cytotoxicity test
In order to calculate the number of cells surviving after radiation treatment, -irradiated cells were seeded at a concentration of 250 cells/60mm dish in 5 dishes for each treatment and plates were incubated at 37°C in a 5% CO2 humidified atmosphere. After 48h, cells were washed with phosphate-buffered saline (PBS) and fresh culture medium was added. Cells were maintained in culture for 10 days, with bi-weekly medium changes, then were fixed with methanol, stained with 10% aqueous Giemsa and scored for colony formation. Only colonies containing more than 50 cells were counted. Untreated BALB/c3T3 cells were carried out as negative controls. Results were expressed as the mean number of colonies/plate (± Standard Error, S.E.), absolute clonal efficiency (ACE.), i.e., the fraction of cells which survived after treatment with respect to the number of seeded cells, and relative clonal efficiency (RCE.), which estimates the percentage of reduction of cell clonal efficiency in treated groups compared to non-exposed control cells.
Cell Transformation Assay
Exponentially growing BALB/c3T3 cells were g-irradiated (0.25-3Gy) to estimate the occurrence of transformed foci. Cells were seeded at a concentration of 3 x 104 cells/60 mm dish, 20h after the irradiation, in 10 dishes for each exposure. Plates were incubated at 37°C in a 5% CO2 humidified atmosphere for 48 h. Then, cells were washed with PBS. Unexposed BALB/c3T3 cells were carried out as negative controls. Cells were maintained in culture for 5 weeks, with bi-weekly medium changes, then were fixed with methanol, stained with 10% aqueous Giemsa and scored for the formation of transformed foci. Only foci considered as positive, i.e., not smaller than 1 mm deeply basophilic, showing a dense layer formation and a random orientation of cells at the edges, were counted. Data was reported as transformation frequency (TF), calculated on the cells that survived after radiation exposure [30].
Microarray: RNA Extraction and Hybridization
Total RNA was isolated 20 h after IR exposure by using TRIzol Reagent (Invitrogen, San Diego, CA, USA) followed by purification with Rneasy affinity column (Qiagen, Valencia, CA, USA), according to the manufacturer’s instructions. RNA quality was analysed with Agilent Bioanalyzer 2100 (RNA 6000 Nano Lab Chip, Agilent Technologies, Palo Alto, CA, USA).
cDNA was synthesized from 20µg of total RNA and directly labelled with Cy3-dCTP and Cy5-dCTP, following the manufacturer protocol optimized for the use with Agilent oligo microarray Kit. Labelled cDNAs from the control and treatment reactions were combined and purified by QIA quick spin column (Qiagen), then hybridized into the Mouse (V2) oligo microarray slide according to the Agilent 60-mer oligo microarray processing protocol. Slides were scanned with Agilent dual laser Microarray Scanner, including both Cy-3 and Cy-5 channels. Scanned images were analysed by Feature Extraction software 8.1 (Agilent Technologies) to derive raw intensity data used in the next analysis step.
Microarray Experimental Design and Data Analysis
Two doses were chosen to perform microarray experiments on the basis of results from the preliminary cytotoxicity test: the lowest assayed dose 0.25Gy, which fostered the cell clonal efficiency, and the intermediate dose 3Gy, responsible for significant cytotoxic effects, although still compatible with cell survival. We used four slides for each IR-dose to have four technical replicates, two of which were performed in dye-swap. Slide’s control and saturated features were filtered out. Filtered raw intensity data sets (median green and red signal) were then log2-transformed and normalized intra-array with Joint-LOWESS algorithm. Transformed data were analysed by MAANOVA (Micro Array Analysis Of Variance) data analysis package of R programming environments. A fixed-effect linear ANOVA model was chosen to fit transformed data (variance due to different sources, including array, dye, gene, variety and their combined effects). We then tested a null hypothesis of no differential expression using F statistics computed on the James-Stein shrinkage estimates of the error variance. To avoid any assumption on error distribution, we computed p values for hypothesis tests via permutation methods (1000 permutations were carried out). The false-discovery rate controlling method was used to correct significance estimate for multiple testing hypothesis. MAANOVA analysis returned a list of 2523 differential expressed genes in the 3Gy sample versus 0.25Gy dose treatment. Gene modulations were calculated as the mean log ratio difference between the two sets of data.
Quantitative Real-Time PCR (qPCR)
An amount of 1µg of total RNA was reverse-transcribed at 50°C for 30 in 30µl reaction mixture containing SuperScriptTM III Reverse Transcriptase, RNaseOUT Recombinant Ribonuclease Inhibitor, oligo(dT) (2.5µM), random hexamers (2.5ng/µl), MgCl2 (10mM), and NTPs (SuperScriptTM III First-Strand Synthesis SuperMix for qRT-PCR, Invitrogen). Then, 2.5 ng cDNA were amplified with the appropriate primer pair (200nM) by using SYBR GreenER qPCR SuperMix for iCycler (Invitrogen), in the same cycling conditions (50°C for 2, 95°C for 8, 45 cycles at 95°C, 15 and at 60°C, 60). The primers for aurora kinase A (Aurka), calmodulin 1 (Calm-1), centromere protein A (Cenpa), phorbol-12-myristate-13-acetate-induced protein 1 (Pmaip1, Noxa), were designed with Primer 3 software [36], while those for the growth arrest and DNA-damage-inducible 45 alpha (Gadd45a), cyclin-dependent kinase inhibitor 1A (Cdkn1a, p21) and glyceraldehyde-3-phosphate dehydrogenase (Gapdh) were obtained from the Quantitative PCR Primer Database [37]. Primers sequences were used at a concentration of 200nM. Sequences are listed in Table 1. All amplifications were run at least in triplicate with an iCycler IQ Multicolor Real Time PCR detection system (Biorad, Hemel Hempstead, UK). Reactions were characterized at the time point when the fluorescence signal of the PCR product was first detected (quantification Cycle). Gene target Cq values were normalized to Gapdh and mRNA levels for each gene calculated as 2-Cq, by using the non-exposed sample as calibrator.
Table 1: Forward and reverse sequences of Quantitative Real-Time PCR primers.
Results
Cytotoxicity test
A preliminary cytotoxicity test was performed, in order to evaluate the effects of a wide range of doses, ranging from 0.25Gy to 6Gy, on the ability of BALB/c3T3 to form colonies and to select the working doses to be used in the CTA. At the lowest tested dose (0.25Gy), the statistically significant increase of the cell clonal efficiency was observed. As the IR dose increased, the clonal efficiency linearly decreased, with the almost complete inhibition of the colony formation at 6Gy of exposure. The calculation of ACE. and RCE. further confirmed this biphasic cell response to the tested IR doses [Table 2].
Table 2: Mean number of colonies, Absolute Clonal Efficiency (ACE) and Relative Clonal Efficiency (RCE) of BALB/c 3T3 cells after IR exposure.
aData are reported as a mean of five replicates.
bSignificantly different (p<0.01) from controls at the t-Student test.
cSignificantly different (p<0.01) from controls at the z-test.
dSignificantly different (p<0.01) from controls at the Chi-square test.
Transformation assay
Within the tested range of doses (0.25-3Gy), statistically significant increase of the transformation frequency (TF) (p<0.01, Poisson test) was observed only at 1.5, 2 and 3Gy doses. The 0.25Gy dose, which was able to induce an increment of clonal efficiency, did not induce cell transformation, as compared to the negative controls. Results are shown in Figure 1.
Figure 1: Transformation frequency of Balb/c 3T3 cells at each tested dose. **p-value < 0.01, by the Poisson test after comparison with untreated cells.
Microarray Biological Interpretation
To investigate the biological processes that resulted from the list of differentially expressed genes between 3 and 0.25Gy IR exposures, we used the Expression Analysis Systematic Explorer (EASE) analysis [38]. EASE is a powerful tool for the interpretation of microarray data, that gives gene ontology (GO) annotations of gene lists with statistics adjusted for multiple comparisons using a conservative variant of Fisher exact probability (EASE score), that penalizes the significance of categories supported by few genes and favors more robust biological themes [39]. We selected the main biological processes, according to an EASE score threshold of 0.01, corrected for multiple comparisons with the Bonferroni method, (p-value < 0.01). The results are summarized in Table 3. We obtained several categories within biological processes sustained by the list of up-modulated genes, but not by the down-expressed genes.
The most significant GO term that distinguished the 3Gy exposure from the lowest dose treatment was mitotic cell cycle, GO: A0000278, together with its parent category cell cycle (GO: A0007049, (EASE score <0.01, and p-value < 0.01). This term was enriched by genes involved in the spindle assembly checkpoint (SAC), like Bub1b, Bub3 and Mad2l1, that have the role of preventing premature separation of sister chromatids and mitotic exit throughout a metaphase block. Other important genes belonging to this gene category were Kif22, Kif23 and Kif2c, which encode for KIF proteins, a superfamily of microtubule-based molecular motors involved in intracellular organelle transport and chromosomes movement during cell division. Moreover, Kif2c seems to be required to coordinate the onset of sister centromere separation during anaphase chromosome segregation.
Table 3: Main biological processes returned by EASE analysis of genes ratios obtained by direct comparison between 3Gy and 0.25Gy irradiation.
Mitotic cell cycle GO term was also supported by several genes specifically attributed to S phase of mitotic cell cycle, (GO: A0000084), (Dut, Mcm2, Mcm 3, Mcm 5, Mcm 6, Mcm 7, Nfic, Nfix, and Top2a) and M phase of mitotic cell cycle, (GO: A0000087), (Anapc5, Ccna2, Ccnb2, Cdc20, Fzr1 and Pttg1) showing an EASE score 0.01, but a p-value >0.05. All these genes seem to promote cell cycle progression, and some are overexpressed in tumors.
Within the biological processes significantly modulated by 3Gy versus 0.25Gy dose treatments, EASE also identified gene networks for biosynthesis (GO: A0009058, EASE score <0.01 and p <0.01), with its child GO terms sterol metabolism, (GO: A0016125), steroid biosynthesis, (GO: A000669), and sterol biosynthesis, (GO: A0016126), which clearly converge on cholesterol biosynthesis, (GO: A0006695). Indeed, all the genes belonging to this GO term codify for enzymes that have a key role in catalyzing the most important steps of cholesterol formation and operate both in the pre-squalenic phase (Idi1, Fdps, Hmgcr, Hmgcs1, Hmgcs2, MVD), and in the squalenic phase (Dhcr24, Dhcr7, Ebp). The up-modulation of HMG-CoA reductase (Hmgcr), which represents the rate-limiting enzyme for cholesterol synthesis, was also found.
Full details of the roles of all genes mentioned in this article can be conveniently obtained via the hyperlinks to NCBI Entrez in the Gene Section [40].
Quantitative RT-PCR Validation
Real-Time PCR analysis was performed on selected modulated genes (Aurka, Cenpa, Calm1, Gadd45a, Pmaip1, Cdkn1a) returned from MAANOVA analysis of microarray data. The validation was performed across each sample tested in our study by using biological replicates. Indeed, cells were exposed to each IR dose used in microarray experiments and the RNA was re-extracted and tested in Real-Time PCR.
The relative gene expression was determined by comparing the mRNA level obtained from cells exposed to 3Gy and 0.25Gy with the mRNA level of the untreated cells, by calculating the 2-Cq, where Cq= Cq gene-Cq GAPDH and Cq: Cq exposed – Cq control. We used a Cq t-test to establish the statistical significance of each comparison between exposed and not exposed cells. Statistically significant (p<0.05) Real-Time PCR expression data were substantially consistent with microarray results [Table 4].
Table 4: Real-Time PCR validation of six modulated genes detected by microarray. Underlined modulations indicate the statistically significant values obtained in Real-Time PCR that confirm microarray results.
NE = not exposed.
*p-value <0. 05 at the t- test after comparison between exposed cells/nonexposed cells in Real-Time PCR.
**p-value <0.01 at the t-test after comparison between exposed cells/non-exposed cells in Real-Time PCR
Discussion
In this work, we integrated the in vitro BALB/c3T3 CTA with a global gene expression approach to obtain further information on the molecular events sustaining dose-related cytotoxic and transforming processes induced by IR. We focused on two IR doses, 0.25Gy and 3Gy, which led to completely different outcomes in the BALB/c3T3 cell system. Indeed, a significant increase in RCE and the absence of cell transformation was found in cells exposed at 0.25Gy. Conversely, in response to the 3Gy dose treatment, cell growth was reduced by about 40% while 15-fold increase of TF was observed. These data substantially confirmed the transforming ability of high-dose IR on BALB/c3T3 cells as previously reported by other authors [41,42]. Unlike the well-defined response to the high-dose irradiations, the biological effects at low doses are considerably more complex. Mammalian cells have been reported to show low-dose hyper-radiosensitivity (HRS) at doses below 20-30cGy [43,44], while increased radioresistance (IRR), was observed at doses comprised in the 20-70cGy range in vitro [43]. HRS is related to increased cell death, due to the failure of cell machinery to activate the early G2/M checkpoint and, consequently, DNA repair. Conversely, IRR is associated with increased clonal survival, following the induction of a minimum damage, which is able to activate a protective response [43]. The lowest tested dose in our model (0.25Gy) is included in the range of doses, which can induce either HRS or IRR. Cell survival at this dose was increased by almost 3 times with respect to non-irradiated cell survival. This result is consistent with possible not monotone dose-response relationship (NMDR), often related to hormetic effects, which have been described to occur at this level of IR doses and ascribed to an adaptive response (AR) [42].
The AR is commonly known as a dose-related phenomenon characterized by the ability of pre-exposure at low doses to sustain the cells response to subsequent exposure to radiation high doses. However, an AR can also be observed after a single low dose exposure, not followed by exposure to high doses, triggered by substances released in the culture medium [43]. Some authors have demonstrated that even spontaneous cell transformation can be reduced by low doses of radiation, due to the increase of clonal survival [45,46]. It was previous reported that the rate of neoplastic transformation in non-tumorigenic Hela x human skin fibroblasts, treated with IR doses between 100-300mGy, was at the same level as that of the untreated cells [46]. This observation supports our results.
To have a better comprehension of the different outcomes elicited by the exposure to 0.25Gy and 3Gy IR doses, microarray experiments were performed. The list of statistically different transcripts obtained from MAANOVA analysis of the 3Gy vs 0.25Gy dose treatments underwent EASE biological interpretation.
One of the most significant biological processes belongs to the gene category of mitotic cell cycle, which was supported by several genes related to the spindle assembly checkpoint (SAC) or to mitotic progression. SAC is the major cell cycle control mechanism in mitosis, ensuring the accurate chromosome segregation by delaying anaphase, until all chromosomes have properly attached to spindle microtubules. The mechanism responsible for SAC silencing is not completely understood. It seems to be the consequence of the occupancy of all kinetochores by microtubules and, probably, of the tension from bipolar attachment [47,48]. Among the key components and regulators of SAC, Bub1b, Bub3 and Mad2l1 were overexpressed in the 3Gy vs 0.25Gy data set. The products of these genes, in particular Bub1b, inhibit the protein ubiquitin ligase Anaphase Promoting Complex/Cyclosome (APC/C-Cdc20), by targeting its coactivator protein Cdc20 for degradation, until all chromosomes are correctly attached to all microtubules [48,51], Bub1b is also the core SAC protein responsible for the recruitment of other SAC components for assembling to kinetocores [51]. Defects in the SAC machinery or its silencing, and the consequent chromosome miss-aggregation, cause cell death in most cases, and also malignant transformation [52]. Mad2l and Bub1b were found to be overexpressed in several cancers. The association of their overexpression with tumorigenesis and chromosome instability (CIN) was supported by the increasing non-disjunction events and aneuploidy [53,54]. These genes have been also considered prognostic markers of tumor progression [55]. Beside these SAC genes, the 3Gy vs 0.25Gy data set also showed the up-regulation of Cdc20, whose degradation is the initial response to SAC activation. Without this early step in the checkpoint establishment, consisting, in Cdc20 deregulation, the existing levels of CDC20 could reach the threshold for APC/C activity causing cells to progress through mitosis [56,58]. Moreover, CDC20 overexpression was found to induce mice tumorigenesis and is associated with human cancer [59,61]. Therefore, in our model, CDC20 may abrogate SAC, possibly leading to CIN, which may represent the critical step in the TF increase observed after 3Gy IR exposure, compared to the lowest dose. This hypothesis is further confirmed by the overexpression of Ccna2 and Ccnb2 genes, which was detected in our study, and was reported to be associated with CIN in colorectal cancer cell lines [62,63].
IR is not only associated with cancer. Studies on high-dose exposures in Japanese population, from the A-bomb follow-out, showed that IR are equally responsible for both cancer and cardiovascular disease (CVD). Moreover, IR doses between 0.5 and 2Gy were associated with an elevated risk of stroke and heart disease, comparable with the risk of cancer, and equally contributing to one third of deaths due to irradiation [64]. Some studies showed the role of IR low doses in stroke and CVD, as reported by Zielinski et al. [65] who analyzed a Canadian cohort exposed to a range of doses between 5-400mSv [65]. Cancer radiotherapy has been associated with cardiovascular diseases, especially in patients affected with breast and lung cancers [66]. Indeed, the effects of radiotherapy in cancer cells are mediated by the generation of reactive oxygen species (ROS), which activate p53-mediated apoptosis. Heart damage has been ascribed to the high production of ROS by mitochondria in cardiomyocytes and endothelial cells after radiation treatment, leading to mitochondrial dysfunction [66]. The activation of NOXA-mediated pathway has been described as an essential event in IR-mediated mitochondrial dysfunction [67]. NOXA is one of the genes modulated in our model, supporting the hypothesis that NOXA activation represents an early key event in the response to IR-induced damage and that our model is able to highlight molecular mechanisms related to adverse outcomes. Indeed, the epidemiological findings alone, without any further information of the key events at the cellular and molecular level, are not sufficient to give a clear evidence for a direct association between the exposure to IR low doses and the adverse outcome of CVD or heart disease [68,69]. Cardiovascular disease due to IR exposure was associated with hypercholesterolemia and specifically to an increase of low-density lipoprotein cholesterol [70,72]. Interestingly, our results showed the modulation of different GO categories converging on cholesterol biosynthesis in the 3Gy vs 0.25Gy data set. The overexpression of HMG-CoA Reductase, the rate limiting enzyme in the synthesis of total cholesterol, pointed at the accumulation of cholesterol together with other genes like the farnesyl diphosphate synthase (FDPS). FDPS, catalyzes the sequential condensation of three 5-carbon isoprene units to form farnesyl pyrophosphate, and the squalenic phase transcripts Dhcr7, Dhcr24 and Ebp, whose defects are associated with syndromes characterized by cholesterol deficiency [73,75].
Beside the role of hypercholesterolemia in sustaining CVD, cholesterol may also play a role in cancer. Indeed, since cholesterol is an essential component of lipid rafts, a more elevated amount of this molecule is likely going to alter raft-dependent signaling in tumors cells, favoring cell survival mechanisms [76]. It can be speculated that cholesterol biosynthesis may be the target of the IR and play a role in different outcomes according to the dose of exposure.
In conclusion, we suggest that the cytotoxicity test and the cell transformation assay on BALB/c3T3 cells have the potential to screen low and high IR exposures and represent an interesting model system for the evaluation of specific biomarkers of effect, when integrated with –omics technologies. Indeed, the use of transcriptional profiling highlighted a set of gene transcripts specifically triggered by the 3Gy vs 0.25Gy dose and defined specific biological processes that could discriminate the response to high doses of IR compared to low levels of exposure. Moreover, the results are not limited only to the carcinogenesis process but could interestingly include other pathologies.
Acknowledgement
Authors would like to thank Dr. Daniele Quercioli for his help in the statistical analysis
Conflict of Interest Statement
The authors declare no conflicts of interest.
References
- Roy L, Gruel G, Vaurijoux A (2009) Cell response to ionising radiation analysed by gene expression patterns. Ann Ist super sanita 45(3): 272-277.
- WHO (2012) Preliminary dose estimation from the nuclear accident after the 2011 Great East Japan earthquake and tsunami. Geneve, Switzerland.
- Walsh L, Zhang W, Shore RE, Auvinen A, Laurier D, et al. (2014) A framework for estimating radiation-related cancer risks in Japan from the 2011 Fukushima nuclear accident. Radiat Res 182(5): 556-572.
- http://www.who.int/ionizing_radiation/en/
- Grosovsky AJ, de Boer JG, de Jong PJ, Drobetsky EA, Glickman BW (1988) Base substitutions, frameshifts, and small deletions constitute ionizing radiation-induced point mutations in mammalian cells. Proc Nati Acad Sci U S A 85(1): 185-188.
- Little JB (2000) Radiation carcinogenesis. Carcinogenesis 21(3): 397-404.
- Nikitaki Z, Hellweg CE, Georgakilas AG, Ravanat JL (2015) Stress-induced DNA damage biomarkers: applications and limitations. Front Chem 3: 35.
- Garte SJ, Burns FJ (1991) Oncogenes and radiation carcinogenesis. Environ Health Perspect 93: 45-9.
- Santoro M, Carlomagno F (2013) Oncogenic rearrangements driving ionizing radiation-associated human cancer. J clin invest 123(11): 4566-4568.
- Sishc BJ, Nelson CB, McKenna MJ, Battaglia CL, Herndon A, et al. (2015) Telomeres and Telomerase in the Radiation Response: Implications for Instability, Reprograming, and Carcinogenesis. Front Oncol 5: 257.
- Shin SC, Lee KM, Kang YM, Kim K, Lim SA, et al. (2011) Differential expression of immune-associated cancer regulatory genes in low- versus high-dose-rate irradiated AKR/J mice. Genomics 97(6): 358-363.
- Amundson SA, Lee RA, Koch Paiz CA, Bittner ML, Meltzer P, et al. (2003) Differential responses of stress genes to low dose-rate gamma irradiation. Mol Cancer Res 1(6): 445-452.
- Rieger KE, Chu G (2004) Portrait of transcriptional responses to ultraviolet and ionizing radiation in human cells. Nucleic Acids Res 32(16): 4786-4803.
- Gurung A, Uddin F, Hill RP, Ferguson PC, Alman BA (2009) Beta-catenin is a mediator of the response of fibroblasts to irradiation. Am J pathol 174(1): 248-255.
- Spitz DR, Hauer Jensen M (2014) Ionizing radiation-induced responses: where free radical chemistry meets redox biology and medicine. Antioxid Redox Signal 20(9): 1407-14009.
- Amundson SA, Grace MB, McLeland CB, Epperly MW, Yeager A, et al. (2004) Human in vivo radiation-induced biomarkers: gene expression changes in radiotherapy patients. Cancer Res 64(18): 6368-6371.
- Fachin AL, Mello SS, Sandrin Garcia P, Junta CM, Ghilardi Netto T, et al. (2009) Gene expression profiles in radiation workers occupationally exposed to ionizing radiation. J Radiat Res 50(1): 61-71.
- Morandi E, Severini C, Quercioli D, Perdichizzi S, Mascolo MG, et al. (2009) Gene expression changes in medical workers exposed to radiation. Radiat Res 172(4): 500-508.
- Roudkenar MH, Li L, Baba T, Kuwahara Y, Nakagawa H, et al. (2008) Gene expression profiles in mouse liver cells after exposure to different types of radiation. J Radiat Res 49(1): 29-40.
- Dai JM, Sun DC, Lin RX, Yang J, Lou S, et al. (2006) Microarray analysis of differentially expressed genes in mouse bone marrow tissues after ionizing radiation. Int J Radiat Biol 82(7): 511-521.
- Taki K, Wang B, Nakajima T, Wu J, Ono T, et al. (2009) Microarray analysis of differentially expressed genes in the kidneys and testes of mice after long-term irradiation with low-dose-rate gamma-rays. J Radiat Res 50(3): 241-252.
- Bang HS, Choi MH, Kim CS, Choi SJ (2016) Gene expression profiling in undifferentiated thyroid carcinoma induced by high-dose radiation. J Radiat Res 57(3): 238-249.
- Sokolov M, Nguyen V, Neumann R (2015) Comparative Analysis of Whole-Genome Gene Expression Changes in Cultured Human Embryonic Stem Cells in Response to Low, Clinical Diagnostic Relevant, and High Doses of Ionizing Radiation Exposure. Int J Mol Sci 16(7): 14737-14748.
- Lourenço J, Mendo S, Pereira R (2016) Radioactively contaminated areas: Bioindicator species and biomarkers of effect in an early warning scheme for a preliminary risk assessment. J H Mater 317: 503-542.
- Chauhan V, Kuo B, McNamee JP, Wilkins RC, Yauk CL (2016) Transcriptional benchmark dose modeling: Exploring how advances in chemical risk assessment may be applied to the radiation field. Environ Mol Mutagen 57(8): 589-604.
- Annys E, Billington R, Clayton R, Bremm KD, Graziano M, et al. (2014) Advancing the 3Rs in regulatory toxicology – Carcinogenicity testing: Scope for harmonisation and advancing the 3Rs in regulated sectors of the European Union. Regul Toxicol pharmacol. 69(2): 234-242.
- Jacobs MN, Colacci A, Louekari K, Luijten M, Hakkert BC, et al. (2016) International regulatory needs for development of an IATA for non-genotoxic carcinogenic chemical substances. ALTEX 33(4): 359-392.
- Jacobs MN, Colacci A, Corvi R, Vaccari M, Aguila MC, et al. (2020) Chemical carcinogen safety testing: OECD expert group international consensus on the development of an integrated approach for the testing and assessment of chemical non-genotoxic carcinogens. Arch Toxicol 94(8): 2899-2923.
- Combes R, Balls M, Curren R, Fischbach M, Fusenig N, et al. (1999) Cell transformation assays as predictors of human carcinogenicity. Altern Lab Anim 27(5): 745-767.
- Mascolo MG, Perdichizzi S, Rotondo F, Morandi E, Guerrini A, et al. (2010) BALB/c 3T3 cell transformation assay for the prediction of carcinogenic potential of chemicals and environmental mixtures. Toxicol In Vitro 24(4): 1292-1300.
- Sakai A (2007) BALB/c 3T3 cell transformation assays for the assessment of chemical carcinogenicity. AATEX. 2007,14(SI): 367-373.
- Vasseur P, Lasne C (2007) OECD Detailed review paper on cell transformation assay for detection of chemicals carcinogens. Mutant Res 744(11): 8-11.
- IARC/NCI/EPA Working Group (1985) Cellular and Molecular Mechanisms of Cell Transformation and Standardization of Transformation Assays of Established Cell Lines for the Prediction of Carcinogenic Chemicals: Overview and Recommended Protocols. Cancer Research 45(5): 2395-2399.
- Ao L, Liu JY, Gao LH, Liu SX, Yang MS, et al. (2008) Differential expression of genes associated with cell proliferation and apoptosis induced by okadaic acid during the transformation process of BALB/c 3T3 cells. Toxicol In Vitro 22(1): 116-127.
- Maeshima H, Ohno K, Tanaka Azuma Y, Nakano S, Yamada T (2009) Identification of tumor promotion marker genes for predicting tumor promoting potential of chemicals in BALB/c 3T3 cells. Toxicol In Vitro 23(1):148-157.
- http://bioinfo.ut.ee/primer3/
- http://web.ncifcrf.gov/rtp/gel/primerdb/
- https://david.ncifcrf.gov/
- Hosack DA, Dennis G, Sherman BT, Lane HC, Lempicki RA (2003) Identifying biological themes within lists of genes with EASE. Genome Biol 4(10): R70.
- https://www.ncbi.nlm.nih.gov/gene/
- Little JB (1979) Quantitative studies of radiation transformation with the A31-11 mouse BALB/3T3 cell line. Cancer Res 39(5): 1474-1480.
- Tauchi H, Sawada S (1993) Suppression of gamma- and neutron-induced neoplastic transformation by ascorbic acid in Balb/c 3T3 cells. Int J Radiat Biol 63(3): 369-374.
- Bonner WM (2004) Phenomena leading to cell survival values which deviate from linear-quadratic models. Mutat Res 568(1): 33-39.
- Marples B (2004) Is low-dose hyper-radiosensitivity a measure of G2-phase cell radiosensitivity? Cancer Metastasis Rev 23(3-4): 197-207.
- Azzam EI, de Toledo SM, Raaphorst GP, Mitchel RE (1996) Low-dose ionizing radiation decreases the frequency of neoplastic transformation to a level below the spontaneous rate in C3H 10T1/2 cells. Radiat Res 146(4): 369-373.
- Redpath JL, Elmore E (2006) Radiation-induced neoplastic transformation in vitro, hormesis and risk assessment. Dose Response 5(2): 123-130.
- Fuller BG, Stukenberg PT (2009) Cell division: righting the check. Curr Biol 19(14): R550-553.
- Etemad B, Kops GJ (2016) Attachment issues: kinetochore transformations and spindle checkpoint silencing. Curr Opin Cell Biol 39: 101-108.
- Rao CV, Yamada HY, Yao Y, Dai W (2009) Enhanced genomic instabilities caused by deregulated microtubule dynamics and chromosome segregation: a perspective from genetic studies in mice. Carcinogenesis 30(9): 1469-1474.
- Kulukian A, Han JS, Cleveland DW (2009) Unattached kinetochores catalyze production of an anaphase inhibitor that requires a Mad2 template to prime Cdc20 for BubR1 binding. Dev Cell 16(1): 105-17.
- Yu L, Shang ZF, Abdisalaam S, Lee KJ, Gupta A, et al. (2016) Tumor suppressor protein DAB2IP participates in chromosomal stability maintenance through activating spindle assembly checkpoint and stabilizing kinetochore-microtubule attachments. Nucleic Acids Res 44(18): 8842-8854.
- Burum Auensen E, DeAngelis PM, Schjølberg AR, Røislien J, Mjåland O, et al. (2008) Reduced level of the spindle checkpoint protein BUB1B is associated with aneuploidy in colorectal cancers. Cell Prolif 41(4): 645-659.
- Sotillo R, Hernando E, Díaz-Rodríguez E, Teruya Feldstein J, Cordón Cardo C, et al. (2007) Mad2 overexpression promotes aneuploidy and tumorigenesis in mice. Cancer cell 11(1): 9-23.
- Yamamoto Y, Matsuyama H, Chochi Y, Okuda M, Kawauchi S, et al. (2007) Overexpression of BUBR1 is associated with chromosomal instability in bladder cancer. Cancer Genet Cytogenet 174(1): 42-47.
- Bargiela Iparraguirre J, Prado Marchal L, Pajuelo Lozano N, Jiménez B, Perona R, et al. (2014) Mad2 and BubR1 modulates tumourigenesis and paclitaxel response in MKN45 gastric cancer cells. Cell Cycle. 13(22): 3590-3601.
- Ge S, Skaar JR, Pagano M (2009) APC/C- and Mad2-mediated degradation of Cdc20 during spindle checkpoint activation. Cell Cycle 8(1): 167-171.
- Alfieri C, Chang L, Zhang Z, Yang J, Maslen S, et al. (2016) Molecular basis of APC/C regulation by the spindle assembly checkpoint. Nature 536(7617): 431-436.
- Di Fiore B, Wurzenberger C, Davey NE, Pines J (2016) The Mitotic Checkpoint Complex Requires an Evolutionary Conserved Cassette to Bind and Inhibit Active APC/C. Mol Cell 64(6): 1144-1153.
- Mondal G, Sengupta S, Panda CK, Gollin SM, Saunders WS et al (2007) Overexpression of Cdc20 leads to impairment of the spindle assembly checkpoint and aneuploidization in oral cancer. Carcinogenesis 28(1): 81-92.
- Li M, Fang X, Wei Z, York JP, Zhang P (2009) Loss of spindle assembly checkpoint-mediated inhibition of Cdc20 promotes tumorigenesis in mice. J Cell Biol 185(6): 983-994.
- Marucci G, Morandi L, Magrini E, Farnedi A, Franceschi E, et al. (2008) Gene expression profiling in glioblastoma and immunohistochemical evaluation of IGFBP-2 and CDC20. Virchows Arch 453(6): 599-6090.
- Dunican DS, McWilliam P, Tighe O, Parle McDermott A, Croke DT (2002) Gene expression differences between the microsatellite instability (MIN) and chromosomal instability (CIN) phenotypes in colorectal cancer revealed by high-density cDNA array hybridization. Oncogene 21(20): 3253-3257.
- Sarafan Vasseur N, Lamy A, Bourguignon J, Le Pessot F, Hieter P, et al. (2002) Overexpression of B-type cyclins alters chromosomal segregation. Oncogene 21(13):2051-2057.
- Shimizu Y, Kodama K, Nishi N, Kasagi F, Suyama A, et al. (2010) Radiation exposure and circulatory disease risk: Hiroshima and Nagasaki atomic bomb survivor data, 1950-2003. BMJ 340: b5349.
- Zielinski JM, Ashmore PJ, Band PR, Jiang H, Shilnikova NS, et al. (2009) Low dose ionizing radiation exposure and cardiovascular disease mortality: cohort study based on Canadian national dose registry of radiation workers. Int J Occup Med Environ Health 22(1): 27-33.
- Fu X, Tang J, Wen P, Huang Z, M Najafi (2021) Redox interactions-induced cardiac toxicity in cancer therapy. Arch Biochem Biophys 708: 108952.
- Young Woo Seo, Jin Na Shin, Kang Hee Ko, Jong Hee Cha, Jae Yoon Park, et al. (2003) The Molecular Mechanism of Noxa-induced Mitochondrial Dysfunction in p53-Mediated Cell Death. J Biol Chem 278(48): 48292-48299.
- Baselet B, Rombouts C, Benotmane AM, Baatout S, Aerts A (2016) Cardiovascular diseases related to ionizing radiation: The risk of low-dose exposure (Review). Int J Mol Med 38(6): 1623-1641.
- Boerma M, Sridharan V, Mao XW, Nelson GA, Cheema AK et al (2016) Effects of ionizing radiation on the heart. Mutat Res 770(Pt B): 319-327.
- Baker JE, Fish BL, Su J, Haworth ST, Strande JL, et al. (2009) 10 Gy total body irradiation increases risk of coronary sclerosis, degeneration of heart structure and function in a rat model. Int J Radiat Biol 85(12): 1089-1100.
- Little MP, Gola A, Tzoulaki I (2009) A model of cardiovascular disease giving a plausible mechanism for the effect of fractionated low-dose ionizing radiation exposure. PLoS Comput Biol 5(10): e1000539.
- Kumarathasan P, Vincent R, Blais E, Saravanamuthu A, Gupta P, et al. (2013) Cardiovascular changes in atherosclerotic ApoE-deficient mice exposed to Co60 (γ) radiation. PloS one 8(6): e65486.
- Platt FM, Wassif C, Colaco A, Dardis A, Lloyd Evans E, et al. (2014) Disorders of cholesterol metabolism and their unanticipated convergent mechanisms of disease. Annu Rev Genomics Hum Genet 15: 173-194.
- Zhou Y, Zhang X, Chen L, Wu J, Dang H, et al. (2008) Expression profiling of hepatic genes associated with lipid metabolism in nephrotic rats. Am J Physiol Renal Physiol 295(3): F662-671.
- http://omim.org/
- Freeman MR, Solomon KR (2004) Cholesterol and prostate cancer. J Cell Biochem 91(1): 54-69.