Skip to main content

Accumulation of DNA methylation alterations in paediatric glioma stem cells following fractionated dose irradiation



Radiation is an important therapeutic tool. However, radiotherapy has the potential to promote co-evolution of genetic and epigenetic changes that can drive tumour heterogeneity, formation of radioresistant cells and tumour relapse. There is a clinical need for a better understanding of DNA methylation alterations that may follow radiotherapy to be able to prevent the development of radiation-resistant cells.


We examined radiation-induced changes in DNA methylation profiles of paediatric glioma stem cells (GSCs) in vitro. Five GSC cultures were irradiated in vitro with repeated doses of 2 or 4 Gy. Radiation was given in 3 or 15 fractions. DNA methylation profiling using Illumina DNA methylation arrays was performed at 14 days post-radiation. The cellular characteristics were studied in parallel.


Few fractions of radiation did not result in significant accumulation of DNA methylation alterations. However, extended dose fractionations changed DNA methylation profiles and induced thousands of differentially methylated positions, specifically in enhancer regions, sites involved in alternative splicing and in repetitive regions. Radiation induced dose-dependent morphological and proliferative alterations of the cells as a consequence of the radiation exposure.


DNA methylation alterations of sites with regulatory functions in proliferation and differentiation were identified, which may reflect cellular response to radiation stress through epigenetic reprogramming and differentiation cues.


Paediatric high-grade glioma, including glioblastoma, is one of the most aggressive types of brain tumours. The molecular subtype, determined by DNA methylation profiling, is strongly associated with clinical features [1]. In contrast to adult glioblastoma where most cases have extensive accumulated genetic aberrations, paediatric glioblastoma is characterised by epigenetic deregulation of functions involved in developmental biology [2].

Radiation is an important therapeutic tool for children diagnosed with brain tumours and used as standard treatment for children from the age of 4–5. However, radioresistance remains a major obstacle for radiotherapy treatment [3]. Genetic and epigenetic variabilities, which are currently not fully understood, influence, both directly and indirectly, the cellular response to radiation. Adverse late effects from radiation therapy is common but differ markedly on an individual basis, with one third of patients suffering from severe cognitive disabilities and one third having a seemingly normal development [4].

The biochemical processes that occur within a cell following irradiation are well known, but alterations of the epigenome by radiation have not been thoroughly investigated. Radiation strikes and disrupts the DNA molecule directly by ionization and indirectly by the production of free radicals. Both these actions take place very quickly, in less than a second [5]. However, the biological effects become apparent in hours or days. Ionizing radiation disrupts chemical bonds and impacts specific parts of the DNA molecule, which challenges the cell to continue replication and cell division. If the damage is sufficient to kill the cell, death can take place days or weeks later. Crucial cells as haematopoietic stem cells are the definitive dose-limiting factor of all radiotherapy treatment plans [6]. Still, the dose must be sufficient to eradicate all tumour cells with the potential to repopulate the tumour.

DNA methylation influences phenotypic variation including radiation response. As the cellular radio-response aims to protect the cell, specific DNA methylation alterations may define mechanisms mediating cell survival after radiation exposure and hence constitute a target for radiosensitization [7]. Profiling of altered DNA methylation patterns following irradiation may identify reversible stress adaptations strategies, as cell lineage-specific programs controlled by epigenetic mechanisms could be one way for the cells to resist stress and survive treatment [8, 9]. Cancer stem cells, with their self-renewal capacity and ability to maintain an undifferentiated state, are reported to be particularly resistant against radiation and are thus crucial for brain tumour repopulation and progression [10, 11]. Understanding the cellular and molecular mechanisms of these cancer stem cells involved in the response and resistance to radiotherapy is important to improve therapy.

The overall aim of this study was to identify radiation-induced cytosine methylation alterations; initially (3 FDIR) and after an accumulated methylation drift over repeated radiation dose fractionations (15 FDIR). Identifying targets and dynamics of DNA methylation alterations following radiation has the potential to reveal biological consequences. For this, we used patient-derived paediatric high-grade glioma stem cells (GSCs) grown under conditions that preserve the signatures of the originating tumour and maintain the tumour-initiating capacity of the cells [12]. Genome-wide DNA methylation maps of pre- and post-irradiated GSC cultures were generated to determine the effects radiation has on epigenomic plasticity and stability.


Patients and cell cultures

Cell cultures established from five paediatric patients (three males and two females aged 4–12) with high-grade glioma tumours were used. Cultures originated from pre-radiotherapy tissues except for one that was derived from a tumour relapse post-radiotherapy. All patients received conventional doses of 1.8 Gy in daily fractions (5 days a week) and a total dose of 54 Gy as part of their clinical management. Patients, tumours and cell cultures have been fully described previously [12]. Samples were collected with written informed consent according to the ethical approval (Dnr 604-12). Cells were cultured on laminin-coated (0.5 μg/cm2, Sigma-Aldrich, Sweden, or Biolamina, Sundbyberg, Sweden) plastics (Costar, Washington) in stem cell media (DMEM-F12 supplemented with B27 (Gibco, Thermo Fisher, Waltham, Massachusetts), N2 (Gibco) and EGF (10 ng/ml, Peprotech, Stockholm, Sweden)). Fresh media was replaced every fourth day, and cells were passaged typically every week, 1:2–1:5, using Accutase (Gibco). All established cell cultures were confirmed negative for mycoplasma contamination.


GSCs were irradiated with repeated doses of 2, 4 or 8 Gy (Fig. 1) at a dose rate of 1.7 Gy/min delivered by a RS 2000 Biological Research Irradiator (Rad Source Technologies, Buford, GA). A 0.3-mm Cu filter for low energy radiation filtration at 160 kV × 25 mA was used. As controls, unirradiated cells from the corresponding passages were handled outside the cell incubator in parallel with irradiated cells. The dose and number of fractions of radiation treatment in this study were determined by proliferation assays and DMP analysis based on preliminary experiments. The aim was to analyse cells that resume cell cycle progression following radiation and to avoid early non-stable DNA methylation effects and early cell death that lack interest for the generation of radiation-resistant phenotypes. Dose-response experiments were performed initially, and 2 Gy was chosen as the dose presented if not otherwise specified. Cells were irradiated for 3 or 15 numbers of fraction dose irradiation (FDIR). The cells that recovered after 15 FDIR were re-irradiated for another 3 fractions (15 + 3) to study this specific population of cells. Cells were assayed at either 72 h, 14 days or 21 days post-radiation.

Fig. 1.
figure 1

Schematic representation of the model system and experimental setup of fractionated dose irradiation schemes. a Paediatric high-grade glioma stem cell cultures established pre- or post-radiotherapy (30 fractions irradiation of 1.8 Gy each) were used in the study. b Cells were exposed to 3–18 fractionated dose irradiations 2–8 Gy. To monitor the stability during culturing, GU-pBT-07 cells were cultured and analysed at passage 12 and 22. These sub-cultured cells from early (p12) and later passages (p22) were irradiated and used as biological replicates

Cell viability and cellular characterisation

Morphological features were monitored during culturing and at passaging, comparing radiation-exposed cells with control cells. Representative images were captured over time with an inverted microscope equipped with phase-contrast optics (Leica) using a black and white CCD video camera (Lumenera, Ottawa, Ontario, Canada). The Via1-Cassette™ was used together with the NucleoCounter® NC-3000™ for determination of cell count, cell size and viability. For adherent cells, the LIVE/DEAD® Cell Imaging Kit (Invitrogen) was used to identify cells with intact or permeable membranes. Hoechst (1:2000; Thermo Fisher) was used as a counterstain for cell nuclei. Images were acquired with Operetta High-Content Imaging System, and data processing was performed using the Harmony software version 4.1 (Perkin Elmer, Waltham, Massachusetts). FDIR cells and the corresponding unirradiated cell population were analysed for phenotypic growth pattern and expression of stem cell markers with immunocytochemistry (ICC) method. Cells were seeded in 96-well plates. Proliferation was assessed with the Click-iT EdU AlexaFluor 488 Imaging Kit (C10337, Thermo Fisher, Waltham, Massachusetts). Briefly, cells were incubated during 4 or 24 h with EdU (5 μmol/L final concentration) prior to fixation with paraformaldehyde 4% (Histolab, Västra Frölunda, Sweden), in phosphate-buffered saline (PBS, 0.1 M, pH 7.4) for 10 min, at 14 days after irradiation. After washing quickly twice and once for 30 min with PBS, the cells were permeabilised 30 min with PBST (PBS with 0.1% Triton X-100). Cells were washed once with PBST and incubated with block solution (3% goat serum (Fisher Scientific, Hampton, New Hampshire) and 1% BSA (bovine serum albumin; Sigma, Saint-Louis, Missouri) in PBS for 1 h at room temperature. Primary antibodies used were as follows: mouse anti-α-tubulin (1/1000; T6199, Sigma) and anti-GFAP (1/1000; G3893, Sigma). All primary antibodies were diluted in block solution. The cells were washed with PBS and incubated 1 h at room temperature with the secondary antibodies. Secondary antibodies were chosen according to the primary antibodies: AlexaFluor594-coupled anti-mouse antibody (A21125, Invitrogen, Carlsbad, CA), AlexaFluor647-coupled anti-rabbit antibody (A121244, Invitrogen) and AlexaFluor594-coupled anti-rabbit antibody (A11012, Invitrogen). All secondary antibodies were diluted at 1/1000 in block solution or fluorescence antibody diluent (FAD; FAD901L, Biocare Medical, Pacheco, CA) for the OLIG2 secondary antibody. Cells were washed with PBS and counterstained with DAPI (1/5000; D9542, Sigma) in PBS for 5 min.

DNA sample preparation and DNA methylation profiling

DNA was extracted using the DNeasy Blood & Tissue Kit (Qiagen, Hilden, Germany) according to the protocol provided by the supplier. The DNA was quantified with a Qubit fluorometer (Life Technologies, Carlsbad, CA), and 500 ng was used for bisulfite modification with the EZ DNA methylation kit (D5001, Zymo Research, Orange, CA). The bisulfite-modified DNA was applied to the Infinium HumanMethylation450 BeadChips (Illumina, San Diego, CA) which determine the methylation levels of more than 450,000 CpG sites. The fluorescence signals were measured from the BeadArrays using an Illumina BeadStation GX scanner and further processed with Bead studio to obtain raw data for further downstream analysis.

Data and statistical analysis

Array data were analysed using the R package ChAMP [13]. Probes with multiple and nonspecific mapping from their intended target were removed before normalisation. Each CpG site was assigned a score called a “β-value”, which corresponds to the ratio between the fluorescence signal of the methylated allele (C) and the sum of the fluorescent signals of the methylated (C) and unmethylated (T) alleles. DNA methylation levels were depicted as beta values in the interval 0–1. Differential methylation with delta beta below 0.2 was omitted before further analysis. Chromosome coordinates refer to human genome assembly hg19. Methylation age was determined with Horvath age [14]. The GISTIC module was used to identify regions of the genome that were significantly amplified or deleted ( Graphs were made in R, Pandas Seaborn and Graph Pad. The Gene Set Enrichment Analysis (GSEA) software was used to test for the over-representation of Gene Ontology (GO) biological processes and molecular functions ( Annotation of Infinium DNA Methylation BeadChip probes was provided by the HumanMethylation450 v1.2 Manifest.

The R package Repetitive Element Methylation Prediction (REMP) version 1.10.0 [15] was used to calculate methylation β-values of Alu elements.


In this study, we profiled DNA methylation alterations that follow after few (n = 3) and many (n = 15 or 15 + 3) fractions of radiation. For this, we used previously generated GSC cultures [12] established from patients before (n = 4) and after treatment (n = 1). These well-characterised cells have been shown to be stable during culturing and maintain their DNA methylation profiles and DNA copy number alterations (CNA). The expression of stem cell proteins as well as in vivo tumour-initiating capacity of these cells [16] demonstrates their resemblance to cancer stem cells or progenitor cells. These cells hence provide a relevant and robust model system for functional studies on paediatric brain tumours. The overall scheme of the study is depicted in Fig. 1.

Proliferation is reduced and cellular morphology is altered following radiation

We first determined whether irradiation changed the cellular morphology of the cells using immunocytochemistry (ICC; Fig. 2). An exponential reduction in cell number as well as alterations in the morphology of the remaining cells (reshaped nuclei, increased polyploidy and senescence morphology) was seen with increasing doses of irradiation. The proportion of cells with active replication, estimated with EdU incorporation, indicated that 2- and 4-Gy-treated cells retained proliferative capacity, while only a small proportion of 8-Gy-treated cells were actively proliferating 72 h post-IR (Fig. 2). In all but one GSC culture, the number of cells with enlarged nuclei and altered morphology was increased (Fig. 2). Surviving cells were able to establish cell colonies when single cells were seeded, demonstrating that the cells maintained the capacity to repopulate (data not shown).

Fig. 2
figure 2

Radiation-induced cellular and nuclear morphological changes. Immunocytochemistry showing cellular morphology with α-tubulin-stained cytoskeleton (orange), DNA-synthesizing cells by EdU incorporation (green) and nuclear staining with DAPI (blue) at normal culture conditions (0 Gy) and exposed to radiation (4 and 8 Gy). Scale bars 100 μm

Short-term doses of radiation do not alter genome-wide DNA methylation levels of GSCs

We next profiled DNA methylation of cells following three fractions of radiation to assess the DNA methylation instability. Supervised hierarchical clustering of the top 1000 most variable sites showed that the individual GCS cultures clustered separately, independent of the irradiation treatment, demonstrating that the genome-wide methylation profiles did not alter substantially (Fig. 3a). DNA methylation can be used to accurately predict biological age. As the predicted age can be altered by external stimuli, we analysed whether the methylation age of the cells was altered with radiation. For this, we used the Horvath’s epigenetic clock, a tissue-independent clock that aims to predict chronological age using methylation data [14]. We detected very small alterations in the predicted DNA methylation age following three FDIR (Fig. 3b).

Fig. 3
figure 3

Genome-wide and site-specific methylation alterations following three FDIR of 2 Gy. a Unsupervised clustering and b DNA methylation age calculated with Horvath’s DNA methylation age calculator show that cells maintain their DNA methylation profiles of specific sites with radiation. Error bars represent standard error of mean, n = 2. c Number of differentially methylated positions (DMPs) detected in each cell line. d Venn diagrams depicting the overlapping DMPs between the cell lines. e A selected CpG site showing a change from unmethylated to methylated state following radiation at both early and late passages and at both 2 and 4 Gy

Few site-specific DNA methylation alterations following three irradiation fractions

Since the global profile was stable, we next analysed whether there were site-specific alterations with radiation. For this, we identified the number of differentially methylated positions (DMPs) that occurred with radiation for the five GCS cultures respectively (Fig. 3c). In all cases, a low number of DMPs were identified (n = 10–23). Few of these were detected in several of the individual GSC cultures (Fig. 3d). One of the identified DMPs, cg26027669, in the growth factor independent 1 transcriptional repressor gene (GFI1), showed large, consistent alterations in methylation frequency with both 2 and 4 Gy radiation in one of the GSC cultures, GU-pBT-07 (Fig. 3e). This site was unmethylated in the unirradiated cultures (methylation beta value < 20%) and highly methylated (> 80%) in the irradiated cultures and could hence constitute a tumour-specific marker of radiation exposure, independent of dose. To validate this alteration, we repeated the experiment using cells at a higher passage which identified the same DNA methylation alteration.

Extended dose fractionations alter DNA methylation profiles and induce thousands of DMPs

As we detected few alterations, and to more mimic the clinical situation, we irradiated the most radioresistant GSC culture GU-pBT-07, with repeated fractions (n = 15). The 15 FDIR treatment was used to evaluate the population of cells that survived and repopulated during treatment. These cells were then irradiated for another three fractions (15 + 3) to study the response of the radioresistant cells. The change in DNA methylation profile following repeated fractions of irradiation was evident when the data was visualised in an unsupervised clustering (Fig. 4a). DNA methylation age calculations based on Horvath’s tissue independent CpG sites showed accelerated age following repeated irradiation fractions (Fig. 4b). We next analysed DMPs following the repeated fractions of 2 Gy irradiation and detected a 1000-fold increase in the number of sites compared to the three fraction set up (Fig. 4c). Comparing the number of DMPs between the 15 and 15 + 3 fractions showed that the majority of sites were consistently detected in both treatment groups (data not shown). The DMPs with the largest difference in methylation value (delta beta) are presented in Additional file 1: Table S1. To determine whether the shifted DNA methylation profiles following 15 and 15 + 3 FDIR was a consequence of the radiation exposure and not general cell culture artefacts, we cultured unirradiated cells for the same time as the duration of the irradiations. This resulted in very few changes, demonstrating that the alterations following 15 and 15 + 3 fractions were indeed due to the irradiation and not a culture artefact (data not shown). In addition, we inferred chromosomal copy number status from the array data which verified the chromosomal stability of the cells during the experiment (Additional file 2: Figure S1).

Fig. 4
figure 4

Alterations in DNA methylation levels following 3, 15 and 15 + 3 fractional doses. a Unsupervised clustering of methylation beta values. b DNA methylation age. c Number of hypomethylated (white) and hypermethylated (grey) DMPs observed in the IR samples compared with unirradiated samples. Error bars represent standard error of mean. To address the comparability of the cells between different passages, experiments were performed at both mid and late passages for both control and 3 FDIR

The differentially methylated positions are enriched in enhancer and R-DMR regions

By analysing the chromosomal distribution of the DMPs following 15 and 15 + 3 FDIR, we detected enrichment at specific chromosomes, for example at chromosome 18 (Fig. 5a). The location of the DMPs was enriched in enhancer regions, regions of repetitive DNA elements and in regions previously demonstrated to be differentially methylated in cancer (C-DMRs) and during reprogramming (R-DMRs) [17] (Fig. 5b). Less sites, in relation to the representation of sites on the array, were detected in the CpG islands and around the transcription start sites. We also noted an overrepresentation of sites annotated as involved in alternative splicing among the identified DMPs (1.5-fold enrichment, BH-corrected p value of 9.7 × 10−4). Pathway analysis using the DMPs identified relevant terms such as neurogenesis, development and cell differentiation (Fig. 5c).

Fig. 5
figure 5

Distribution of DMPs in unirradiated versus 15 + 3 FDIR samples. a Chromosomal location. b Genomic features. c Gene Ontology enrichment analysis (cellular compartments, molecular functions and biological processes). d hsa-miR genes. e Genes with higher number of DMP than expected. f Differentially methylated regions (DMRs) implicated in known radiation response functions

MicroRNAs (miRNAs) typically repress gene expression, and modulation of specific miRNAs may provide therapeutic targets. The methylation arrays include 35,699 CpG sites annotated to the promoter regions of non-coding RNA genes. The majority of the 82 identified DMPs in regions associated with miRNAs increased in methylation following radiation, and only eight showed decreased methylation levels. Several miRNAs that have previously been related to stem cells, invasion, migration and radiation response were identified as altered (Fig. 5d). The top hypomethylated miRNA, after irradiation, was miR-525. This miRNA has previously been demonstrated to be upregulated in response to radiation [18].

We further analysed genes enriched with DMPs and found among these, for example, two histone genes (HIST1H1D and H2AFJ), a gene involved in DNA damage repair (NHEDC1) and a gene in stem cell regulation (OLIG2) (Fig. 5e). Moreover, two HSP70 heat shock protein genes, HSPA1A and HSPA1B, were identified (Fig. 5f). Among the 100 CpG sites with the largest radiation-induced difference in methylation, three sites in the gene SHROOM3, encoding for the actin-binding protein and cell morphology regulator, were identified (Additional file 1: Table S1). The Sphingosine-1-Phosphate Transporter (SP1) gene MFSD2B, a regulator of self-renewal and differentiation in neural and cancer progenitor cells [19], was also differentially methylated at multiple CpG sites (n = 6). Other genes with many altered sites included the glial fibrillary acidic protein GFAP, which encodes one of the major intermediate filament proteins of mature astrocytes, and the glutathione metabolism gene GSTP1.

Frequent hypomethylation of repetitive elements occurs following radiation

The majority of methylated CpG sites in the genome are located within regions of repetitive DNA elements, and we therefore specifically analysed probes in these regions. Methylation β-values of repetitive Alu elements were predicted using the R package REMP. These analyses identified global hypomethylation in Alu elements following irradiation (Additional file 2: Figure S2A). The Pearson correlation in methylation beta values between untreated cells at different passages, as well as between untreated cells and cells irradiated with three fractions, was high (0.99). Following 15 or 15 + 3 fractions, the correlation decreased with an increased number of repetitive base pairs in the probe sequence down to 0.95 (Additional file 2: Figure S2B). This indicates that there is an increased methylation variability with an increasing number of repetitive base pairs. Our results thus confirm previous studies showing that irradiation results in hypomethylation of repetitive elements [20].

Fractionated irradiation alters cellular properties

Cell states are determined by DNA methylation patterns, and we thus explored if the large alterations in the methylome that was induced with fractionated irradiation affected the characteristics of the cells. Irradiation resulted in a higher resistance of the radiation-enriched cells following 15 + 3 fractions compared to 3 fractions (Fig. 6a), and both 3 and 15 + 3 irradiated cells showed increased cell diameter (Fig. 6a, b). Extensive alterations between unirradiated cells and cells irradiated with 3 or 15 + 3 fractions were detected, including reorganisation of the cytoskeleton, altered proliferation rate and growth patterns (Fig. 6c). Colonies that formed after single-cell seeding showed an inter-population variation but a high homogeneity within colonies in terms of cellular morphology, proliferation rate and expression of neural and progenitor cell marker proteins (Fig. 6d). The colony-forming ability following 21 days after irradiation decreased after 3 fractions; however, colonies obtained from radiation-enriched cells (15 + 3 fractions) were dense and larger in size compared to unirradiated cells (Fig. 6d).

Fig. 6
figure 6

Radiation-induced morphological changes with increasing number of radiation fractions of 2 Gy. a Bright-field enlarged cells. b Percent of viable cells and cell diameter. c Immunocytochemistry showing cellular morphology with α-tubulin-stained cytoskeleton (orange), DNA-synthesizing cells by EdU incorporation (green) and nuclear staining with DAPI (blue). d Single-cell seeding after irradiation shows reduced clonogenic capacity after 3 fractions compared to unirradiated and 15 + 3 irradiated cells. e The frequency and staining pattern of GFAP varied between treatment groups

As we detected methylation alterations that could be associated with alternative splicing of the astrocytic gene GFAP, we further evaluated this finding on the protein level using ICC. GFAP was distributed more evenly in the cytoplasm after irradiation (Fig. 6e). This pattern has previously been linked to the different isoforms of GFAP [21, 22], suggesting that the differential methylation that we detected could be involved in the regulation.


Cellular response to radiation and intrinsic radiosensitivity are highly influenced by the function of tumour suppressor genes, often epigenetically silenced in tumours [23]. In this study, we therefore aimed to identify loci of DNA methylation alterations in paediatric GSCs following radiation in order to reveal alterations in sites or regions which could mediate the effects of radiation.

The effects of irradiation on gene expression have been documented in some studies, however resulting in very few consistent and informative clinical useful findings [24, 25]. In addition, discrepancy between protein levels and number of RNA molecules is frequently observed [26]. DNA methylation, on the other hand, is a stable mechanism which specifies cell identity and shows tumour entity-specific profiles. DNA methylation response following radiation treatment has previously been reported but not for stem cells or paediatric tumour cells [27,28,29]. DNA methylation in prostate cancer cells was recently reported to be stable following a single dose of irradiation with few stably altered CpG sites up to 14 days following exposure [30].

Paediatric high-grade gliomas originate from cells along the neural stem cell lineages [31]. These cells have the capacity to generate a variety of different cell types: activated stem cells, progenitor cells and mature cells such as astrocytes, oligodendrocytes and neurons [32]. The DNA methylation pattern is strictly regulated during differentiation and development allowing each cell type to acquire its required phenotype to meet its specific functions and needs. In this study, we used patient-derived GSCs to better reflect the plasticity of cells seen within tumours in vivo. We started with three fractions of irradiation which identified few altered loci (DMPs; > 20% difference in methylation). However, individual DMPs with large alterations could be consistently identified. For example, a CpG site in the gene GFI1 showed a large difference in methylation level. The change from unmethylated to methylated indicates that the majority of cells altered methylation status at this site following radiation. This alteration was consistently found in several experiments and with different doses of irradiation and could therefore be of interest for elucidating cell type-specific early radiation responses. The methylation state of this site has previously been shown to alter in a leukaemia cell line undergoing differentiation [33].

To more mimic the clinical treatment protocols, we next performed many fractions (15 or 15 + 3) of irradiation, which resulted in large DNA methylation alterations, not detected in unirradiated cells cultured in parallel experiments. Genome-wide alterations, an accelerated DNA methylation age and a high number of detected DMPs (~ 10,000), most of which were hypermethylated, were detected. The genomic regions that were most frequently targets of the detected radiation-induced DNA methylation alterations were regions previously described to be altered in reprogramming (R-DMRs), cancer (C-DMRs) and in enhancers [17]. Alterations were also detected in CpG island shores, regions known to be variable across cell types and variably methylated during cell differentiation [34]. Regions with fewer alterations than expected included promoter regions and CpG islands.

Inhibition of specific miRNAs has been shown to contribute to radiation resistance in various cancer cells and to control the balance between proliferation and differentiation following irradiation [35]. The miR-525 contained the most hypomethylated site among the miRNA sites annotated on the array. Inhibition of miR-525-3p has previously been shown to elevate radiosensitivity, while overexpression of the precursor miR-525-3p conferred radioresistance in endothelial-like cells and various tumour cell lines [18]. Two DMPs were located within miR-10b, and previous studies have shown that high gene expression of this miRNA is associated with disease progression and self-renewal [36]. Functional studies have suggested that miR-10b may function as a tumour suppressor and is silenced by DNA hypermethylation in gastric cancer [37].

Apart from alterations in miRNAs, we also detected many DMPs in protein-coding genes. Several of these were located in the genomic regions of the two histone genes H2AFJ and HIST1H1D. Histone variant H2A.J accumulates in senescent cells and promotes inflammatory gene expression [38] which implies that the DNA methylation alterations could affect cell fate. The linker histone HIST1H1D interacts with linker DNA between nucleosomes and functions in the compaction of chromatin into higher order structures [39]. We additionally identified many altered sites with large delta beta in the Heat Shock Protein Family A (Hsp70) Member 1A and 1B. Previous studies have shown that the HSPA1A promoter is methylated in cell lines with low gene expression of HSPA1A with expression being restored following treatment with the DNA demethylating agent 5-aza-2′-deoxycytidine [40]. Since heat shock proteins are an attractive target for anticancer therapy, further studies are needed. Also in the SP1 gene MFSD2B, many DMPs were identified. SP1 was recently shown to promote tumour growth by inhibiting apoptosis and conferring chemo- and radiation resistance to cancer cells [41]. Furthermore, radiation treatment was found to recruit SP1 transcription factor in meningioma cells [42]. We also detected promoter hypermethylation of the glutathione metabolism gene GSTP1 in the irradiated cells which is in agreement with a previous study on irradiated leukocytes [43]. The alterations that we identified in these genes could potentially confer adaptation to radio-induced stress.

Knowledge generated from recent studies indicates that changes in DNA methylation in cancer cells impact the organisation of the genome in the nucleus rather than directly affecting gene expression [44]. Repetitive elements represent a large portion of the human genome. DNA methylation, together with other chromatin modifications, is most often associated with silencing of transposable elements (TEs) within repetitive regions [45]. The Alu element is commonly found in repetitive regions. Recently, many environmental factors were shown to alter Alu methylation, which could lead to chromosomal instability resulting in chromosomal rearrangements [46]. We found that loci in long repetitive regions were significantly enhanced in DMPs following 15 and 15 + 3 FDIR, demonstrating that also irradiation can induce such alterations. This is in agreement with earlier studies showing decreased DNA methylation in repetitive elements with radiation in various models [47,48,49,50].

To link the DNA methylation findings to cellular properties, we studied the cells with ICC. Altered morphology following radiation was apparent, and biological processes like apical cytoskeletal protein processes within membrane regions were enriched among radiation-induced DMPs. Cells following many fractions of irradiation formed dense colonies with altered morphology in contrast to few fractions where the opposite was seen.

In this study, we could not distinguish between if the identified alterations were due to a selection of cells with specific profiles or to altered DNA methylation levels in cells, since the DNA methylation arrays cannot be performed on individual cells. The irradiation induced methylation alterations, with the majority below 30%, indicating selection of cells with specific profiles. However, sites with large alterations, as were identified for example in the GFI1 gene, are more likely to have resulted from an active methylation alteration of the cells, due to the short time of the acquired alterations. The alterations could theoretically also result from a genetic instability following radiation exposure. The analysis of CNAs did however not detect relative increase or decrease in DNA copy levels, demonstrating that the genome was stable during the experiments.


On the epigenetic level, at least three categories of elements were identified as associated with long-term radiation exposure: (1) repetitive sequences, (2) enhancers and (3) sites that distinguish cell types and differentiation states. Although the identification of changes in DNA methylation levels does not prove functional consequences, it demonstrates the need for additional studies that focus on detecting possible functional implications of such accumulation of alterations. The alterations that accompanied the generated radioresistant cells were stable and mark the cell identity and characteristics of these cells. Cellular characteristics showed growth selection from typical mixed morphology to more homogeneous populations.

Identifying DNA methylation marks that drive radiation-induced changes in cell state provides candidates that can potentially be targeted therapeutically. This would offer a way to interfere with the generation of the radiation-resistant phenotype.

Availability of data and materials

Data sets will be accessible in GEO upon acceptance of this manuscript.



Arthrobacter luteus


Gibco B-27 Serum Free Supplement


Benjamini and Hochberg procedure


Charge-coupled device


Regions of DNA where a cytosine nucleotide is followed by a guanine nucleotide




Dulbecco’s modified Eagle medium


Differentially methylated positions


Differentially methylated regions




Epidermal growth factor


Ham’s F-12 nutrient mixture


Fractionated dose irradiation


Glioma stem cells




Phosphate-buffered saline


  1. Gusyatiner O, Hegi ME. Glioma epigenetics: from subclassification to novel treatment options. Semin Cancer Biol. 2018;51:50–8.

    Article  CAS  PubMed  Google Scholar 

  2. Carén H, Pollard SM, Beck S. The good, the bad and the ugly: epigenetic mechanisms in glioblastoma. Mol Asp Med. 2013;34(4):849–62.

    Article  CAS  Google Scholar 

  3. Peitzsch C, Kurth I, Ebert N, Dubrovska A, Baumann M. Cancer stem cells in radiation response: current views and future perspectives in radiation oncology. Int J Radiat Biol. 2019;95(7):900–11.

    Article  CAS  PubMed  Google Scholar 

  4. Krasin MJ, Constine LS, Friedman DL, Marks LB. Radiation-related treatment effects across the age spectrum: differences and similarities or what the old and young can learn from each other. Semin Radiat Oncol. 2010;20(1):21–9.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Cannan WJ, Pederson DS. Mechanisms and consequences of double-strand DNA break formation in chromatin. J Cell Physiol. 2016;231(1):3–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Lombardi G, Rumiato E, Bertorelle R, Saggioro D, Farina P, Della Puppa A, et al. Clinical and genetic factors associated with severe hematological toxicity in glioblastoma patients during radiation plus temozolomide treatment: a prospective study. Am J Clin Oncol. 2015;38(5):514–9.

    Article  CAS  PubMed  Google Scholar 

  7. Zang L, Kondengaden SM, Che F, Wang L, Heng X. Potential epigenetic-based therapeutic targets for glioma. Front Mol Neurosci. 2018;11:408.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Zhu X, Wang Y, Tan L, Fu X. The pivotal role of DNA methylation in the radio-sensitivity of tumor radiotherapy. Cancer Med. 2018;7(8):3812–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Dabin J, Fortuny A, Polo SE. Epigenome maintenance in response to DNA damage. Mol Cell. 2016;62(5):712–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Testa U, Castelli G, Pelosi E. Genetic abnormalities, clonal evolution, and cancer stem cells of brain tumors. Medical Sciences. 2018;6:4.

    Article  CAS  Google Scholar 

  11. Kim M, Costello J. DNA methylation: an epigenetic mark of cellular memory. Exp Mol Med. 2017;49(4):e322.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Wenger A, Larsson S, Danielsson A, Elbaek KJ, Kettunen P, Tisell M, et al. Stem cell cultures derived from pediatric brain tumors accurately model the originating tumors. Oncotarget. 2017;8(12):18626–39.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Morris TJ, Butcher LM, Feber A, Teschendorff AE, Chakravarthy AR, Wojdacz TK, et al. ChAMP: 450k chip analysis methylation pipeline. Bioinformatics. 2014;30(3):428–30.

    Article  CAS  PubMed  Google Scholar 

  14. Horvath S. DNA methylation age of human tissues and cell types. Genome Biol. 2013;14(10):R115.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Zheng Y, Joyce BT, Liu L, Zhang Z, Kibbe WA, Zhang W, et al. Prediction of genome-wide DNA methylation in repetitive elements. Nucleic Acids Res. 2017;45(15):8697–711.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Larsson S, Wenger A, Dosa S, Sabel M, Kling T, Caren H. Cell line-based xenograft mouse model of paediatric glioma stem cells mirrors the clinical course of the patient. Carcinogenesis. 2018.

  17. Feinberg AP. Epigenetic stochasticity, nuclear structure and cancer: the implications for medicine. J Intern Med. 2014;276(1):5–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Kraemer A, Barjaktarovic Z, Sarioglu H, Winkler K, Eckardt-Schupp F, Tapio S, et al. Cell survival following radiation exposure requires miR-525-3p mediated suppression of ARRB1 and TXN1. PLoS One. 2013;8(10):e77484.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Cimino PJ, Kim Y, Wu HJ, Alexander J, Wirsching HG, Szulzewsky F, et al. Increased HOXA5 expression provides a selective advantage for gain of whole chromosome 7 in IDH wild-type glioblastoma. Genes Dev. 2018;32(7-8):512–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Gong M, Shi W, Qi J, Shao G, Shi Z, Wang J, et al. Alu hypomethylation and MGMT hypermethylation in serum as biomarkers of glioma. Oncotarget. 2017;8(44):76797–806.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Thomsen R, Daugaard TF, Holm IE, Nielsen AL. Alternative mRNA splicing from the glial fibrillary acidic protein (GFAP) gene generates isoforms with distinct subcellular mRNA localization patterns in astrocytes. PLoS One. 2013;8(8):e72110-e.

    Article  CAS  Google Scholar 

  22. Stassen O, van Bodegraven EJ, Giuliani F, Moeton M, Kanski R, Sluijs JA, et al. GFAPdelta/GFAPalpha ratio directs astrocytoma gene expression towards a more malignant profile. Oncotarget. 2017;8(50):88104–21.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Llinàs-Arias P, Esteller M. Epigenetic inactivation of tumour suppressor coding and non-coding genes in human cancer: an update. Open Biol. 2017;7(9):170152.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Ghandhi SA, Smilenov LB, Elliston CD, Chowdhury M, Amundson SA. Radiation dose-rate effects on gene expression for human biodosimetry. BMC Med Genet. 2015;8:22.

    Google Scholar 

  25. Yang M, Sun Y, Xiao C, Ji K, Zhang M, He N, et al. Integrated analysis of the altered lncRNAs and mRNAs expression in 293 T cells after ionizing radiation exposure. Int J Mol Sci. 2019;20:12.

    Google Scholar 

  26. Liu Y, Beyer A, Aebersold R. On the dependency of cellular protein levels on mRNA abundance. Cell. 2016;165(3):535–50.

    Article  CAS  PubMed  Google Scholar 

  27. Maierhofer A, Flunkert J, Dittrich M, Muller T, Schindler D, Nanda I, et al. Analysis of global DNA methylation changes in primary human fibroblasts in the early phase following X-ray irradiation. PLoS One. 2017;12(5):e0177442.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Antwih DA, Gabbara KM, Lancaster WD, Ruden DM, Zielske SP. Radiation-induced epigenetic DNA methylation modification of radiation-response pathways. Epigenetics. 2013;8(8):839–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Miousse IR, Kutanzi KR, Koturbash I. Effects of ionizing radiation on DNA methylation: from experimental biology to clinical applications. Int J Radiat Biol. 2017;93(5):457–69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Sutton LP, Jeffreys SA, Phillips JL, Taberlay PC, Holloway AF, Ambrose M, et al. DNA methylation changes following DNA damage in prostate cancer cells. Epigenetics. 2019:1–14.

  31. Alcantara Llaguno S, Sun D, Pedraza AM, Vera E, Wang Z, Burns DK, et al. Cell-of-origin susceptibility to glioblastoma formation declines with neural lineage restriction. Nat Neurosci. 2019;22(4):545–55.

    Article  PubMed  CAS  Google Scholar 

  32. Lu QR, Qian L, Zhou X. Developmental origins and oncogenic pathways in malignant brain tumors. Wiley Interdiscip Rev Dev Biol. 2019;8(4):e342.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  33. Severin J, Waterhouse AM, Kawaji H, Lassmann T, van Nimwegen E, Balwierz PJ, et al. FANTOM4 EdgeExpressDB: an integrated database of promoters, genes, microRNAs, expression dynamics and regulatory interactions. Genome Biol. 2009;10(4):R39.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  34. Irizarry RA, Ladd-Acosta C, Wen B, Wu Z, Montano C, Onyango P, et al. The human colon cancer methylome shows similar hypo- and hypermethylation at conserved tissue-specific CpG island shores. Nat Genet. 2009;41(2):178–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Liu Z, Huang S. Inhibition of miR-191 contributes to radiation-resistance of two lung cancer cell lines by altering autophagy activity. Cancer Cell Int. 2015;15(1):16.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Sheedy P, Medarova Z. The fundamental role of miR-10b in metastatic cancer. Am J Cancer Res. 2018;8(9):1674–88.

    CAS  PubMed  PubMed Central  Google Scholar 

  37. Li Z, Lei H, Luo M, Wang Y, Dong L, Ma Y, et al. DNA methylation downregulated mir-10b acts as a tumor suppressor in gastric cancer. Gastric cancer : official journal of the International Gastric Cancer Association and the Japanese Gastric Cancer Association. 2015;18(1):43–54.

    Article  CAS  Google Scholar 

  38. Contrepois K, Coudereau C, Benayoun BA, Schuler N, Roux P-F, Bischof O, et al. Histone variant H2A.J accumulates in senescent cells and promotes inflammatory gene expression. Nat Commun. 2017;8:14995.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Gilbert N, Thomson I, Boyle S, Allan J, Ramsahoye B, Bickmore WA. DNA methylation affects nuclear organization, histone modifications, and linker histone binding but not chromatin compaction. J Cell Biol. 2007;177(3):401–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Qi W, White MC, Choi W, Guo C, Dinney C, McConkey DJ, et al. Inhibition of inducible heat shock protein-70 (hsp72) enhances bortezomib-induced cell death in human bladder cancer cells. PLoS One. 2013;8(7):e69509.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Kumar R, Liu APY, Orr BA, Northcott PA, Robinson GW. Advances in the classification of pediatric brain tumors through DNA methylation profiling: from research tool to frontline diagnostic. Cancer. 2018.

  42. Velpula KK, Gogineni VR, Nalla AK, Dinh DH, Rao JS. Radiation-induced hypomethylation triggers urokinase plasminogen activator transcription in meningioma cells. Neoplasia (New York, NY). 2013;15(2):192–203.

    Article  CAS  Google Scholar 

  43. Kuzmina NS, Lapteva NS, Rusinova GG, Azizova TV, Vyazovskaya NS, Rubanovich AV. Gene hypermethylation in blood leukocytes in humans long term after radiation exposure - validation set. Environ Pollut (Barking, Essex : 1987). 2018;234:935-942.

    Article  CAS  PubMed  Google Scholar 

  44. Madakashira BP, Sadler KC. DNA methylation, nuclear organization, and cancer. Front Genet. 2017;8:76.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  45. Deniz Ö, Frost JM, Branco MR. Regulation of transposable elements by DNA modifications. Nat Rev Genet. 2019;20(7):417–31.

    Article  CAS  PubMed  Google Scholar 

  46. Lewis CW, Golsteyn RM. Cancer cells that survive checkpoint adaptation contain micronuclei that harbor damaged DNA. Cell Cycle (Georgetown, Tex). 2016;15(22):3131-3145.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Miousse IR, Chang J, Shao L, Pathak R, Nzabarushimana E, Kutanzi KR, et al. Inter-strain differences in LINE-1 DNA methylation in the mouse hematopoietic system in response to exposure to ionizing radiation. Int J Mol Sci. 2017;18:7.

    Article  CAS  Google Scholar 

  48. Elliott EN, Sheaffer KL, Kaestner KH. The ‘de novo’ DNA methyltransferase Dnmt3b compensates the Dnmt1-deficient intestinal epithelium. eLife. 2016;5:e12975.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Koturbash I, Miousse IR, Sridharan V, Nzabarushimana E, Skinner CM, Melnyk SB, et al. Radiation-induced changes in DNA methylation of repetitive elements in the mouse heart. Mutat Res. 2016;787:43–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Prior S, Miousse IR, Nzabarushimana E, Pathak R, Skinner C, Kutanzi KR, et al. Densely ionizing radiation affects DNA methylation of selective LINE-1 elements. Environ Res. 2016;150:470–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


We thank the Bergö group, Institute of Medicine, Sahlgrenska Academy, University of Gothenburg, for the use of radiation equipment.


This work was supported by the Swedish Cancer Society, the Swedish Children’s Cancer Society, the Swedish Research Council and the Swedish state under the agreement between the Swedish government and the county councils, the ALF-agreement. AD was supported by the Swedish Radiation Safety Authority. Open access funding provided by University of Gothenburg.

Author information

Authors and Affiliations



AD and HC designed the study. AD performed the experimental work with help from KB. AD performed data analysis with input from TK. MT provided clinical data and clinical input to the study. AD and HC wrote the manuscript. All authors approved the final version of the manuscript.

Corresponding authors

Correspondence to Anna Danielsson or Helena Carén.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Regional Ethical Review Board of Gothenburg (Dnr 604-12). Tumour samples were obtained after signed informed consent from the parents or legal guardians of children who underwent surgery at the Sahlgrenska University Hospital.

Consent for publication

Consent is granted through the ethical permission.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1

. The DMPs with the largest difference in methylation value.

Additional file 2: Figure S1

. CNA of potential tumour driver and tumour suppressor genes. Figure S2. Alterations in DNA methylation of CpG sites associated with repetitive sequences following 15 and 15+3 FDIR of the cell line GU-pBT-07. (a) DNA methylation beta values of Alu elements decreased following 15 and 15+3 FDIR compared to unirradiated and 3 FDIR; (b) Pearson correlation between unirradiated and irradiated cells decreased with increased number of repetitive bp within the probe sequences.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Danielsson, A., Barreau, K., Kling, T. et al. Accumulation of DNA methylation alterations in paediatric glioma stem cells following fractionated dose irradiation. Clin Epigenet 12, 26 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: