Skip to main content
  • Research article
  • Open access
  • Published:

Long-term effects of cranial irradiation and intrathecal chemotherapy in treatment of childhood leukemia: a MEG study of power spectrum and correlated cognitive dysfunction



Prophylaxis to prevent relapses in the central nervous system after childhood acute lymphoblastic leukemia (ALL) used to consist of both intrathecal chemotherapy (CT) and cranial irradiation (CRT). CRT was mostly abolished in the eighties because of its neurotoxicity, and replaced with more intensive intrathecal CT. In this study, a group of survivors treated with CRT before 1983 and another group treated without CRT thereafter are investigated 20–25 years later, giving a much stronger perspective on long-term quality of life than previous studies. The outcomes will help to better understand these groups’ current needs and will aid in anticipating late effects of prophylactic CRT that is currently applied for other diseases. This study evaluates oscillatory neuronal activity in these long-term survivors. Power spectrum deviations are hypothesized to correlate with cognitive dysfunction.


Resting state eyes-closed magnetoencephalography (MEG) recordings were obtained from 14 ALL survivors treated with CT + CRT, 18 treated with CT alone and 35 controls. Relative spectral power was calculated in the δ, θ, α1, α2, β and γ frequency bands. The Amsterdam Neuropsychological Tasks (ANT) program was used to assess cognition in the executive functions domain. MEG data and ANT scores were correlated.


In the CT + CRT group, relative θ power was slightly increased (p = 0.069) and α2 power was significantly decreased (p = 0.006). The CT + CRT group performed worse on various cognitive tests. A deficiency in visuomotor accuracy, especially of the right hand, could be clearly associated with the deviating regional θ and α2 powers (0.471 < r < 0.697). A significant association between decreased regional α2 power and less attentional fluctuations was found for CT + CRT patients as well as controls (0.078 < r < 0.666). Patients treated with CT alone displayed a power spectrum similar to controls, except for a significantly increased level of left frontal α2 power (p = 0.030).


The tendency towards global slowing of brain oscillatory activity, together with the fact that dementia has been reported as a late effect of CRT and the neuropsychological deficiencies currently present, suggest that the irradiated brain might be aging faster and could be at risk for early‐onset dementia. The CT group showed no signs of early aging.

Peer Review reports


Acute lymphoblastic leukemia (ALL) is the most common malignancy diagnosed in children, representing almost one third of all childhood cancers. There is a peak incidence at 2–5 years of age. Over the last decades, cure rates for ALL patients have gone from less than 5% to over 80% due to improved treatment protocols [1]. An essential part of ALL treatment is central nervous system (CNS) prophylaxis, mostly performed by cranial irradiation (CRT) and/or intrathecal chemotherapy (CT). In the Netherlands, CRT was abolished in 1983 [2].

As survival rates for ALL continue to rise, the focus of clinicians is increasingly shifting to the late effects of therapy and related quality of life [3]. CRT patients have a 13.6 times higher chance of developing secondary neoplasms in the irradiated field within 20 years following treatment [4]. Late neurotoxic effects of CRT are also well described. Many studies have reported a causative link between prophylactic CRT and long-term neurocognitive deficits among ALL survivors, for example memory problems, cognitive slowing and attention deficits [58]. Cognitive effects beyond 10 years after treatment are largely unknown [9]. This study provides data from patients, on average 25 years post childhood treatment, giving a much stronger perspective on long-term quality of life than previous studies. These outcomes not only help to better understand these groups’ current needs, but can also aid in anticipating late effects of currently applied prophylactic CRT that is still part of treatment for ALL in some countries, but also for brain tumours and small cell lung cancer.

Neurocognitive late effects of chemotherapy alone are relatively subtle compared to the effects of CRT. However, a significant percentage of children that received chemotherapy alone also suffer from neurocognitive deficits, particularly in the domain of executive functioning (EF) [1013]. The underlying mechanisms of neurotoxicity, caused by either CRT or CT, are not yet fully understood. Structural imaging has revealed smaller white matter volumes in this population, which correlated with cognitive dysfunction [9]. Functional neuroimaging methods like functional magnetic resonance imaging (fMRI), electroencephalography (EEG) or magnetoencephalography (MEG) have rarely been applied in late effect studies of cancer treatment, while these methods are very suitable to study subtle brain dysfunction. Among the available functional imaging techniques, MEG has the unique property of combining high temporal resolution with good spatial resolution. The MEG signal is quantified in spectral power, which indicates the strength of magnetic induction fields generated by oscillatory neuronal activity and is decomposed into frequency bands [14]. Relative power is the percentage of power in any band compared with the sum of the power of all frequency bands. A power spectrum that displays a shift towards relatively more power in the lower frequency bands is generally considered to be a pathological sign of slowing brain activity and is associated with deteriorated cognition [15, 16]. Also, increased δ or θ power is known as a nonspecific sign of brain pathology [17]. In this paper we will describe relative spectral power changes assessed by MEG during resting state with eyes closed. To our knowledge, this is the first study to use MEG in this population. We hypothesize that the pattern of oscillatory activity will deviate from that of controls and that the deviations will correlate with cognitive dysfunction.

In 1989, Tucker et al. conducted a study of long-term neurocognitive effects (11.5 ± 7 years after treatment) of CT + CRT in ALL and Non Hodgkin Lymphoma (NHL) survivors treated at age 15–69 [18]. They compared the effects between a group treated at 25 years of age or younger and a group treated at an older age. In the younger group, abnormalities were more outspoken, indicating higher susceptibility to neurotoxicity at a younger age. Using EEG, they found excess θ activity and slightly slowed α. The abnormalities were interpreted as minor, but relatively more severe in the younger group, reflecting subtle cerebral cortical dysfunction. Therefore, an association between worse outcome and younger age at treatment is hypothesized.



Seventy‐five eligible ALL survivors were identified from patient records of the VU University Medical Center (VUMC) and the Academic Medical Center (AMC) Amsterdam. They were sent a letter explaining the study and were subsequently contacted by phone. Unfortunately, 45% were unwilling to participate or met one or more of the exclusion criteria (i.e. use of centrally acting drugs, active psychiatric disease or symptoms, pre-existing CNS disorders, metal parts in the body or pregnancy). Survivors who were willing to participate were requested to recruit a control (sibling, partner or friend, n = 44). The ALL survivors were treated according to two different protocols developed by the Dutch Childhood Leukemia Study Group (DCLSG). The first group, CT + CRT (n = 18) was treated according to DCLSG protocol ALL-5 (1979–1984), which consisted of standard dose CT (intravenously administered vincristine, prednisone, L-asparaginase, daunorubicine, intrathecal injections of methotrexate (MTX) and prednisolone) and CRT (2500 Gy). The second group, CT (n = 23), was treated according to DCLSG protocol ALL-6 (1984–1988), which consisted of standard dose CT (intravenous vincristine, dexamethasone, L-asparaginase, intrathecal MTX and prednisolone), without CRT. Duration of treatment was approximately 2 years.

Two CT + CRT patients were excluded from data analysis because meningiomas were discovered during assessment. Additionally, data from two CT + CRT patients, five CT patients and nine controls could not be used due to artefacts discussed in the magnetoencephalography section below. None of the CT patients or controls presented with CNS disorders. Eventually, data from 14 CT + CRT patients, 18 CT patients and 35 controls were used in the analysis for this paper. Characteristics of these participants are reported in Table 1.

Table 1 Characteristics of the included subjects

The ethical principles of the Helsinki Declaration were followed and approval was obtained from the Medical Ethics Committee of the VU University Medical Center (n° 2006/200). All participants signed an informed consent form.


MEG was recorded for ten minutes while subjects were inside a magnetically shielded room (Vacuumschmelze GmbH, Germany) using a 151 channel whole-head MEG system (CTF systems Inc., Canada) [19]. A third order software gradient was used with a recording passband of 0.25-125 Hz and a sample frequency of 625 Hz in 49 subjects (assessment around 2009), and a sample frequency of 312.5 Hz in 18 subjects (assessment around 2007) [20]. Both sample frequencies occurred in each of the subject groups. Head position was monitored. MEG was recorded during a no-task, eyes-closed resting state condition. For this study, 149 of the 151 channels could be used. Additionally, 7 channels appeared to produce no signal for some of the subjects and were appropriately excluded from the calculation of the mean global and regional powers per subject. MEG recordings were converted to ASCII files and from these files four artefact free epochs of 13 s per subject were carefully selected by visual inspection by the two first authors (based on consensus). The MEG data were filtered in the following frequency bands: δ (0.5-4 Hz), θ (4–8 Hz), α1 (8–10 Hz), α2 (10–12 Hz), β (13–30 Hz) and γ (30–50 Hz). Hypothetically, the lower and upper α bands are associated with different cognitive processes and many studies have found differential effects in these subbands [21, 22]. There is no such rationale for subdividing the other frequency bands. Data in the δ band were largely corrupted by movement artefacts like breathing. Therefore, all subjects with artefacts in the δ band were excluded and group differences were investigated. Within this artefact‐free dataset (CT + CRT n = 10, CT n = 11, controls n = 12), the groups did not differ on δ power. Therefore, the δ band was disregarded and more subjects with artefact‐free data in the other frequency bands were included. Eventually, only 16 subjects (CT + CRT n = 2, CT n = 5, CON n = 9) were excluded due to movement artefacts in the other frequency bands, mostly caused by eye blinks and heart beats.

Relative instead of absolute powers were analyzed, because absolute powers are partially determined by the distance between the MEG sensor and the measured neural substrate. This would result in misleading outcomes, because part of our population has a smaller head size due to CRT [23]. The data were processed using CTF DataEditor and BrainWave 0.8.80 [19, 24]. Relative band powers were computed for four 13 s epochs for each subject in the five frequency bands, and averaged for each subject. MEG channels were clustered based upon the approximate underlying cortical areas (see Figure 1). Regional means were calculated averaging relative powers of the available channels per cluster, over four epochs.

Figure 1
figure 1

Schematic projection of MEG channels. Head drawing to illustrate the clustering of channels per region. Midline channels were disregarded. L = left, R = right, F = frontal, C = central, P = parietal, O = occipital, T = temporal.

Neurocognitive assessment

The Amsterdam Neuropsychological Tasks (ANT) program was used to assess EF [25]. The computerized ANT provides for highly standardized assessments and automated recordings of speed and accuracy of information processing, attention processes and working memory. It has proven to be helpful in defining neurocognitive deficit profiles in various clinical domains associated with a generally diffuse impact on the brain, such as phenylketonuria, multiple sclerosis, neurofibromatosis, and middle late effects of childhood ALL [12, 13, 2631]. Based on these studies, outcome parameters of baseline response speed, complex visual information processing, sustained attention, work pace and attentional fluctuations, cognitive flexibility (set shifting and inhibition), visuomotor skills, and visuospatial sequential working memory were selected for assessment of this study’s population. Short task descriptions are given in Table 2, and in more detail in Additional file 1. The reliability and validity of these tasks have extensively been described elsewhere [32]. Task parameters that discriminated between the groups were selected for correlation with MEG data.

Table 2 Amsterdam Neuropsychological Tasks program: description of subtests

Statistical analyses

Individuals’ global and regional means of relative power were entered in SPSS software (version 16.0; SPSS Inc, Chicago, IL) for statistical analyses. All powers were transformed according to the following logarithmic function [33]:

power ln = e log power 1 power

Global relative power differences between subject groups were analyzed using a Multivariate General Linear Model (GLM) with group as between-subjects factor and the frequency bands as dependent variables. Age was used as a covariate when significantly correlated with the dependent variable within the healthy control group. Wilks’ Lambda corrected F and p-values are reported. For each frequency band, regional power differences between groups were analyzed using a MANOVA for ten regions (bilateral central, frontal, parietal, occipital and temporal power). Differences between groups were tested with simple contrasts comparing each patient group with controls. Partial eta squared (ηp 2) was computed to estimate effect sizes (weak effect: ηp² ≈ 0.03; moderate: ηp² ≈ 0.06; large: ηp² ≥ 0.14) [34]. Cohen’s d was used for simple contrasts (small effect: 0.2 < d < 0.3; medium: 0.3 < d < 0.8; large: d > 0.8).

Correlations between regional powers and discriminative cognitive variables

The ANT variables that discriminated between patients and controls were correlated with MEG regional power outcomes in the discriminating frequency bands. Pearson’s correlations (r) and their p-values were calculated between regional powers and these neuropsychological variables, controlling for age. When r in the patient group or in the control group was significant, the r-values were compared. First, they were transformed into z'-values according to the Fisher transformation [35]. The difference |z*| between r1 within the patient group and r2 within controls was calculated and significance of that difference was determined by comparing |z*| to normal z-values [35].

Linear regression models of powers predicting cognitive performance

Separately for each frequency band, regional means of relative powers that differed significantly between groups were entered into linear regression analysis predicting cognitive parameters that discriminated between groups, together with age. The backward method was used, using a probability (F) of 0.05 for entry and 0.10 for removal.


Global relative power

After elog transformation, separate inspection per group revealed normal distributions within all frequency bands, except for α1 in controls (SW statistic: 0.926 (p = 0.022)). Based on these results, the elog transformation was considered sufficiently effective to allow for parametric statistical testing.

The groups differed significantly in age (p = 0.001) due to the consecutive time periods the treatment protocols were applied. The correlation analyses within the control group between the powers and age resulted in a significant outcome for θ_ln (r = −0.560, p = 0.0005), so age was used as a covariate in the multivariate power analyses. The GLM results for the mean global relative powers are shown in Table 3. In the CT + CRT group, θ power was slightly increased compared to controls and α2 power was significantly decreased. Within the CT group, global power levels did not differ from controls. To facilitate interpretation, Figure 2 displays the untransformed mean global relative power values per group for all frequency bands. For additional graphs of transformed and untransformed global relative powers, see Additional file 2.

Table 3 Global relative power differences
Figure 2
figure 2

Mean global relative power values. X-axis: Frequency band; y-axis: Mean global relative power (%). Error bars represent one standard deviation.

Regional relative power

For each frequency band, a separate MANOVA was used to test for regional power differences between groups in the ten regions shown in Figure 1. The multivariate tests did not reach significance, with just a trend for the group difference on regional α2 powers (p = 0.067). Group contrasts versus controls within the θ band revealed increased powers for the CT + CRT group in the RC (p = 0.038) and RP (p = 0.029) regions and trends for increases in the LO (p = 0.066) and RO (p = 0.051) regions. This group also displayed two trends for increased levels of α1 power in the LO (p = 0.081) and RO (p = 0.059) regions. Within the α2 band, the CT + CRT group displayed decreased power levels in all regions, except for LF. The CT group only displayed a significant increase of α2 power (p = 0.030) in the left frontal region compared to controls. For a table of all statistics on the multivariate tests and contrasts, see Additional file 3. Significant contrast results for the CT + CRT group versus controls are displayed in Figure 3. The directions of the deviations are colour coded per region (red for increases, blue for decreases). No significant regional Group x Hemisphere interactions were found.

Figure 3
figure 3

Visualization of regional power differences between the CT + CRT group and controls. Differences concern the θ and α2 frequency bands. The directions of the deviations are colour coded per region (red for increases, blue for decreases).

The time that has passed since treatment could hypothetically allow for either recovery or worsening of late effects. Therefore, time since diagnosis was correlated with powers in all frequency bands, separately within each patient group. There were no significant correlations with either global or regional powers. Hypothetically, a younger age could mean more vulnerability to neurotoxicity of treatment. However, no effect of age at diagnosis was found on either global or regional powers.

Group differences on cognition

The cognitive variables are explained in Additional file 1. Age did not correlate significantly with any of the cognitive parameters within the control group and was therefore not used as a covariate. The CT + CRT group performed worse than controls on measures of cognitive flexibility (percentage of errors (P_flex, p = 0.002), attentional fluctuations during sustained attention (SD_sa, p = 0.062), visuomotor accuracy (D_pu, p = 0.049; S_pu, p = 0.075; Dr_pu, p = 0.015; Sr_pu, p = 0.029) and sequential visuospatial working memory (Nit_vs, p = 0.065; Nitco_vs, p = 0.023; Nit-Nitco_vs, p = 0.036). A table of statistics, including effect sizes, can be found in Additional file 4. When the CT group was compared with controls, trends were found for worse sequential visuospatial working memory (Nitco_vs, p = 0.092) and accuracy of inhibition (P_inhib, p = 0.059).

Correlations between regional powers and discriminative cognitive variables

Correlations were calculated within the composite group of CT + CRT patients and controls to investigate general patterns. Significantly differentiating regional powers were correlated with differentiating cognitive variables. Increased levels of RC and RP θ power correlated significantly with worse performance on visuomotor control (D_pu, Dr_pu and Sr_pu). Decreased α2 powers in the eight significant regions of the CT + CRT group (LO, LP, LT, RC, RF, RO, RP and RT) were significantly correlated with better performance on attentional fluctuations (SD_sa) and with worse performance on visuomotor control (D_pu and Dr_pu). Correlations within the CT + CRT group and within controls were separately calculated and then compared (see Additional file 5). There was a trend for the difference between CT + CRT patients and controls in correlation coefficients of attentional fluctuations with right frontal α2 powers. This correlation was absent in the CT + CRT group (r = 0.078), while it was significantly positive in controls (r = 0.549). Correlations between visuomotor accuracy of the right hand (Dr_pu) and regional α2 powers differed significantly in seven regions (not in RF) between CT + CRT patients (where there were significantly negative correlations with r ranging from −0.508 to −0.697)), and controls (where the correlations were not significant and ranging from −0.129 to 0.142). Roughly the same pattern was seen in the θ band, but with positive correlations within the CT + CRT group between Dr_pu and RC θ (r = 0.630) and RP θ (r = 0.599).

To verify whether the increased LF α2 power level in the CT group was associated with the mild deficiencies in accuracy of inhibition and sequential visuospatial working memory, these outcomes were correlated, but no significant results were found. The significant correlations between α2 and attentional fluctuations observed in controls were absent in the CT group.

To explore the possibility that increased levels of θ power reflect the results of compensatory activity, the data were searched for associations between increased regional θ powers and better cognitive performance on any of the neuropsychological tasks. Two trends were found for increased left occipital θ power in association with less attentional fluctuations (r = −0.509, p = 0.075) and better right hand visuomotor accuracy on the tracking task (r = −0.509, p = 0.075). A description of the tracking task can be found in Additional file 1.

Linear regression models

The differentiating regional powers were entered in linear regression analysis for each group separately, and for the composite group of CT + CRT patients and controls. Of the three pursuit variables, visuomotor accuracy of the preferred hand (Dr_pu) is best explained by regional θ power. Within the CT + CRT group, 39.6% variance on Dr_pu could be explained by increased RC θ power. No significant model could be reached for controls.

Attentional fluctuations (SD_sa) could be significantly explained by α2 powers in one or two regions, regardless of group composition. More variance could be explained in the CT + CRT group (R2 = 0.653), than in controls (R2 = 0.345). D_pu and Dr_pu could be significantly explained by α2 powers in multiple regions. Similarly, more variance could be explained in CT + CRT patients (R2 = 0.912 (D_pu) and 0.727 (Dr_pu)) than in controls (R2 = 0.173 (D_pu) and 0.327 (Dr_pu)). An overview of linear regression outcomes is given in Additional file 6.


The global power spectrum of CRT‐treated patients reflects a slight increase of θ activity and a significant decrease of α2 activity. Patients treated with CT alone display a power spectrum similar to controls, except for a significantly increased level of left frontal α2 power.

Late effects of cranial irradiation

The increase in θ and decrease in α2 indicate a tendency towards global slowing of brain oscillatory activity, a finding that is consistent with other recent MEG studies in low grade glioma patients and Alzheimer patients [15, 16]. Increased θ power is also known as a nonspecific sign of brain pathology [17]. Our results confirm Tucker’s findings from 1989 of increased θ and decreased α power in irradiated patients and are also in line with results from a review by Klimesch, published in 1999 [18, 36]. He concluded that α power is lowered and θ power enhanced in subjects with a variety of different neurological disorders. Also in healthy people in the late part of their lifespan this pattern is familiar and is supposed to reflect some degree of normal aging in the form of vascular or fibrillary degeneration [17]. Taken together with the fact that dementia is frequently reported as a late effect of whole brain irradiation and the neuropsychological deficiencies currently present, these findings suggest that the irradiated brain might be aging faster and could be at risk for early‐onset dementia [18, 3740]. Long-term follow‐up of this patient group is therefore increasingly important. Longitudinal research is also recommended.

High percentages of variance of visuomotor accuracy, especially of the right hand, are explained by increased θ and decreased α2 powers in the CT + CRT group. This is a strong indication that these deviant power levels are underlying this particular deficiency. Unexpectedly, decreased regional α2 powers are significantly associated with less attentional fluctuations in both irradiated patients and controls. In general, α is suppressed by stimuli (e.g. light) and mental activities [17]. Rihs et al. (2009) describe how α power increases and decreases co-vary with visual cortex excitability related to anticipatory visual attention processes [4143]. This illustrates that patterns of α power fluctuations are difficult to decipher and it is almost impossible to interpolate between resting state and attentional task data.

The CT + CRT group demonstrates decreased α2 levels and worse attentional fluctuations. Within the right frontal region of irradiated patients, where the α2 power level is decreased, the association with attentional fluctuations is absent, while there is a strong positive association in controls. This suggests that especially the significantly decreased α2 power in the right frontal region underlies the deficiency in sustained attention in the CT + CRT group. Lower α2 is not helping sustained attention here. The α2 associations in the occipital regions are slightly, although not significantly, stronger than in controls, which could hypothetically indicate that the occipital lobes are trying to compensate for decreased α2 power in the right frontal region.

Additionally, recent research on θ activity, e.g. by Voytek et al. (2010) [44], suggests that increased θ power may reflect compensatory activity for a damaged brain region that is being challenged, rather than reflecting pathology. In our study, subtle implications were found for compensatory θ activity in the left occipital lobe. Combined, the increase in left occipital θ power and the decrease in occipital α2 powers both being associated with less attentional fluctuations, suggest that attentional fluctuations in irradiated patients may be partly counterbalanced by compensatory activity in the occipital lobes. One might arguably postulate that the observed trends for elevated occipital α1 power represent nothing more than a side effect of working with relative powers instead of absolute. These slightly increased levels were not associated with any of the cognitive deficiencies in the irradiated patients, although reports have been made about associations between an increase in lower α power and disturbed working memory [15].

Late effects of chemotherapy

The patient group treated with CT alone showed no signs of early aging. The subtlety of the cognitive deficiencies and limited spectral power changes are consistent with overall findings that the patients treated with CT only do not perform significantly worse than controls. The significant reduction of late effects makes CT a superior alternative to CRT, especially since there was no trade-off in terms of mortality or recurrence [1, 2]. Event-free survival at 10 years after treatment with ALL-6 was even 30% better than after treatment with ALL-5 and CNS relapse rates improved from 12.9% to 1.1% [45]. Meanwhile, suspected neurotoxic effects of CT might be masked by compensatory mechanisms. It remains unclear whether the elevated level of left frontal α2 power should be interpreted as compensatory or pathological. If lowered α2 is regarded as compensatory, increased α2 in the CT patients should be regarded as pathological. However, increased α2 power and the mild cognitive deficiencies observed in the CT patients were not associated, which contradicts the pathology hypothesis. In CT patients, α2 power is not correlated with attentional fluctuations, so the absence of attentional fluctuations combined with increased α2 is not necessarily inconsistent with the hypothesis of lower α2 levels being compensatory. Also, the possibility of increased α2 being compensatory in other cognitive domains cannot be ruled out.

The time that passed since end of treatment (20–25 years) and the age of the patients at time of treatment (an age of relatively high plasticity of the brain) could have allowed for alternative brain strategies to develop. For the CT patients, this could have resulted in mostly normal cognitive performance. For the irradiated group, however, neurotoxicity is much more severe. Compensation might exist, yet appears to be inadequate, which emphasizes this group’s need for long-term follow‐up and interventions.

Methodological considerations

Although there were good reasons for analyzing relative powers instead of absolute, this also creates a limitation. A decrease in one frequency band might cause an increase in another band, although this would not be the observed pattern in absolute measures. On the other hand, consistency with the literature and significant associations with cognitive deficiencies support our findings in the CT + CRT group. Another possible limitation of the study is the relatively small sample size which raises the question whether this set of subjects is representative for the whole population. Confirmation by future studies is needed.

No corrections were applied for multiple comparisons, although commonly the Bonferroni method would be applied in order to lower the chances of false discovery of significant results (the type I error rate). In a publication by Perneger et al. (1998), a view widely held by epidemiologists is expressed that this method creates more problems than it solves and the best way of dealing with multiple comparisons is by simply describing what tests of significance have been performed, and why [46]. This way, findings are not interpreted differently according to how many other tests were performed. Also, type I errors cannot be prevented without increasing the number of type II errors, which are equally worth avoiding. In the present study, pre-established hypotheses were tested in our correlation analyses, instead of just correlating everything, and more than 50% of these hypotheses were confirmed. In addition, the correlation results were supported by the regression results. This, combined with reporting effect sizes in addition to significance levels, should provide the reader with sufficient information to base an interpretation of the results on.


Power spectrum analysis of MEG registrations is able to demonstrate abnormal patterns of resting-state oscillatory brain activity in irradiated long-term survivors of ALL. The suspicion of a faster aging brain after CRT warrants careful long-term follow‐up and screening for early‐onset dementia. In future ALL therapy, CRT should be completely abolished.


This study was supported by a grant to LS, AV and CB from the Dutch Cancer Society (UL 2006–3630).



Acute lymphoblastic leukemia




Cranial radiation therapy




Amsterdam Neuropsychological Tasks program


Central nervous system


Executive functioning


Dutch Childhood Leukemia Study Group




  1. Pui C-H, Evans WE: Treatment of acute lymphoblastic leukemia. N Engl J Med. 2006, 354: 166-178. 10.1056/NEJMra052603.

    Article  CAS  PubMed  Google Scholar 

  2. Veerman AJ, Hählen K, Kamps WA, Van Leeuwen EF, De Vaan GA, Solbu G, Suciu S, Van Wering ER, Van der Does-Van der Berg A: High cure rate with a moderately intensive treatment regimen in non- high-risk childhood acute lymphoblastic leukemia. Results of protocol ALL VI from the Dutch Childhood Leukemia Study Group. J Clin Oncol. 1996, 14 (Suppl 3): 911-918.

    CAS  PubMed  Google Scholar 

  3. Cole PD, Kamen BA: Delayed neurotoxicity associated with therapy for children with acute lymphoblastic leukemia. Ment Retard Dev Disabil Res Rev. 2006, 12: 174-183. 10.1002/mrdd.20113.

    Article  PubMed  Google Scholar 

  4. Hijiya N, Hudson MM, Lensing S, Zacher M, Onciu M, Behm FG, Razzouk BI, Ribeiro RC, Rubnitz JE, Sandlund JT, Rivera GK, Evans WE, Relling MV, Pui CH: Cumulative incidence of secondary neoplasms as a first event after childhood acute lymphoblastic leukemia. JAMA. 2007, 297 (Suppl 11): 1207-1215.

    Article  CAS  PubMed  Google Scholar 

  5. Brouwers P, Poplack D: Memory and learning sequelae in long-term survivors of acute lymphoblastic leukemia: association with attention deficits. Am J Pediatr Hematol Oncol. 1990, 12: 174-181. 10.1097/00043426-199022000-00009.

    Article  CAS  PubMed  Google Scholar 

  6. Butler RW, Copeland DR: Neuropsychological effects of central nervous system prophylactic treatment in childhood leukemia: methodological considerations. J Pediatr Psychol. 1993, 18: 319-338. 10.1093/jpepsy/18.3.319.

    Article  CAS  PubMed  Google Scholar 

  7. Butler RW, Copeland DR: Attentional processes and their remediation in children treated for cancer: A literature review and the development of a therapeutic approach. J Int Neuropsychol Soc. 2002, 8: 115-124. 10.1017/S1355617702811110.

    Article  PubMed  Google Scholar 

  8. Precourt S, Robaey P, Lamothe I, Lassonde M, Sauerwein HC, Moghrabi A: Verbal cognitive functioning and learning in girls treated for acute lymphoblastic leukemia by chemotherapy with or without cranial irradiation. Dev Neuropsychol. 2002, 21: 173-195. 10.1207/S15326942DN2102_4.

    Article  PubMed  Google Scholar 

  9. Reddick WE, Shan ZY, Glass JO, Helton S, Xiong X, Wu S, Bonner MJ, Howard SC, Christensen R, Khan RB, Pui CH, Mulhern RK: Smaller white-matter volumes are associated with larger deficits in attention and learning among long-term survivors of acute lymphoblastic leukemia. Cancer. 2006, 106: 941-949. 10.1002/cncr.21679.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Kingma A, Van Dommelen RI, Mooyaart EL, Wilmink JT, Deelman BG, Kamps WA: No major cognitive impairment in young children with acute lymphoblastic leukemia using chemotherapy only: A prospective longitudinal study. Journal of Pediatric Hematology Oncology. 2002, 24: 106-114. 10.1097/00043426-200202000-00010.

    Article  Google Scholar 

  11. Copeland DR, Moore BD, Francis DJ, Jaffe N, Culbert SJ: Neuropsychologic effects of chemotherapy on children with cancer: a longitudinal study. J Clin Oncol. 1996, 14 (Suppl 10): 2826-2835.

    CAS  PubMed  Google Scholar 

  12. Mennes M, Stiers P, Vandenbussche E, Vercruysse G, Uyttebroeck A, De Meyer G, Van Gool S: Attention and information processing in survivors of childhood acute lymphoblastic leukemia treated with chemotherapy only. Pediatric Blood Cancer. 2004, 43: 1-8. 10.1002/pbc.20105.

    Article  Google Scholar 

  13. Buizer AI, de Sonneville LMJ, van den Heuvel-Eibrink MM, Veerman AJ: Chemotherapy and attentional dysfunction in survivors of childhood acute lymphoblastic leukemia: Effect of treatment intensity. Pediatric Blood and Cancer. 2005, 45: 281-290. 10.1002/pbc.20397.

    Article  PubMed  Google Scholar 

  14. Stam CJ: Use of magnetoencephalography (MEG) to study functional brain networks in neurodegenerative disorders. J Neurol Sci. 2010, 289: 128-134. 10.1016/j.jns.2009.08.028.

    Article  CAS  PubMed  Google Scholar 

  15. Bosma I, Stam CJ, Douw F, Bartolomei F, Heimans JJ, van Dijk BW, Postma TJ, Klein M, Reijneveld JC: The influence of low-grade glioma on resting state oscillatory brain activity: a magnetoencephalography study. J Neurooncol. 2008, 88: 77-85. 10.1007/s11060-008-9535-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. De Haan W, Stam CJ, Jones BF, Zuiderwijk IM, van Dijk BW, Scheltens P: Resting-state oscillatory brain dynamics in Alzheimer disease. J Clin Neurophysiol. 2008, 25: 187-193. 10.1097/WNP.0b013e31817da184.

    Article  PubMed  Google Scholar 

  17. Niedermeyer E, da Silva FH L: Electroencephalography: Basic Principles, Clinical Applications and Related Fields. 1993, Baltimore: Lippincott Williams & Wilkins

    Google Scholar 

  18. Tucker J, Prior PF, Green CR, Ede GMV, Stevenson JF, Gawler J, Jamal GA, Charlesworth M, Thakkar CM, Patel P, Lister TA: Minimal neuropsychological sequelae following prophylactic treatment of the central nervous system in adult leukeamia and lymphoma. Br J Cancer. 1989, 60: 775-780. 10.1038/bjc.1989.358.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Analysis Software CTF MEG. Version 5.0.1. Coquitlam, Canada: VSM MedTech Systems Inc, []

  20. Vrba J, Robinson SE: Signal processing in magnetoencephalography. Methods. 2001, 25 (Suppl 2): 249-271.

    Article  CAS  PubMed  Google Scholar 

  21. Klimesch W: Memory processes, brain oscillations and EEG sychronization. Int J Psychophysiol. 1996, 24: 61-100. 10.1016/S0167-8760(96)00057-8.

    Article  CAS  PubMed  Google Scholar 

  22. Klimesch W, Freunberger R, Sauseng P, Gruber W: A short review of slow phase synchronization and memory: evidence for control processes in different memory systems?. Brain Res. 2008, 1235: 31-44.

    Article  CAS  PubMed  Google Scholar 

  23. Appleton RE, Farrell K, Zaide J, Rogers P: Decline in head growth and cognitive impairment in survivors of acute lymphoblastic leukaemia. Arch Dis Child. 1990, 65: 530-534. 10.1136/adc.65.5.530.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Stam CJ: BrainWave: A Java based application for functional connectivity and network analysis. []

  25. De Sonneville LMJ: Amsterdam Neuropsychological Tasks: a computer-aided assessment program. Computers in Psychology: Cognitive ergonomics, clinical assessment and computer-assisted learning, Volume 6. Edited by: den Brinker BPLM, Beek PJ, Brand AN, Maarse SJ, Mulder LJM. 1999, Lisse: Swets & Zeitlinger, 187-203.

    Google Scholar 

  26. Burgard P, Rey F, Rupp A, Abadie V, Rey J: Neuropsychologic functions of early treated patients with phenylketonuria, on and off diet: results of a cross-national and cross-sectional study. Pediatr Res. 1997, 41: 368-374. 10.1203/00006450-199703000-00011.

    Article  CAS  PubMed  Google Scholar 

  27. Huijbregts S, de Sonneville L, Licht R, Sergeant J, van Spronsen FJ: Inhibition of Prepotent Responding and Attentional Flexibility in Treated Phenylketonuria. Dev Neuropsychol. 2002, 22: 481-499. 10.1207/S15326942DN2202_4.

    Article  PubMed  Google Scholar 

  28. Huijbregts SCJ, de Sonneville LMJ, Licht R, van Spronsen FJ, Verkerk PH, Sergeant JA: Sustained attention and inhibition of cognitive interference in treated phenylketonuria: associations with concurrent and lifetime phenylalanine concentrations. Neuropsychologia. 2002, 40: 7-15. 10.1016/S0028-3932(01)00078-1.

    Article  CAS  PubMed  Google Scholar 

  29. De Sonneville LMJ, Boringa JB, Reuling IEW, Lazeron RHC, Adèr HJ, Polman CH: Information Processing Characteristics in Subtypes of Multiple Sclerosis. Neuropsychologia. 2002, 14: 1751-1765.

    Article  Google Scholar 

  30. Rowbotham I, Pit-ten Cate IM, Sonuga-Barke EJS, Huijbregts SCJ: Cognitive control in adolescents with neurofibromatosis type 1. Neuropsychology. 2009, 23: 50-60.

    Article  CAS  PubMed  Google Scholar 

  31. Huijbregts S, Swaab H, de Sonneville L: Cognitive and Motor Control in Neurofibromatosis Type I: Influence of Maturation and Hyperactivity-Inattention. Dev Neuropsychol. 2010, 35: 737-751. 10.1080/87565641.2010.508670.

    Article  PubMed  Google Scholar 

  32. Gunther T, Herpertz-Dahlmann B, Konrad K: Reliability of attention and verbal memory tests with normal children and adolescents - clinical implications. Zeitschrift für Kinder- und Jugendpsychiatrie und Psychotherapie. 2005, 33: 169-179. 10.1024/1422-4917.33.3.169.

    Article  PubMed  Google Scholar 

  33. Gasser T, Bacher P, Mocks J: Transformations towards the normal distribution of broad band spectral parameters of the EEG. Electroencephalogr Clin Neurophysiol. 1982, 53: 119-124. 10.1016/0013-4694(82)90112-2.

    Article  CAS  PubMed  Google Scholar 

  34. Stevens J: Applied multivariate statistics for the social sciences. 1986, London: Lawrence Erlbaum Associates

    Google Scholar 

  35. Neter J, Kutner MH, Wasserman W, Nachtsheim CJ: Applied Linear Statistical Models. 1996, New York: WCB McGraw-Hill, 4

    Google Scholar 

  36. Klimesch W: EEG alpha and theta oscillations reflect cognitive and memory performance: a review and analysis. Brain Res Rev. 1999, 29: 169-195. 10.1016/S0165-0173(98)00056-3.

    Article  CAS  PubMed  Google Scholar 

  37. Bleyer WA: Neurologic sequelae of methotrexate and ionizing radiation: a new classification. Cancer treatment reports. 1981, 65: 89-98.

    PubMed  Google Scholar 

  38. D’Ambrosio DJ, Cohen RB, Glass J, Konski A, Buyyounouski MK, Feigenberg SJ: Unexpected dementia following prophylactic cranial irradiation for small cell lung cancer: case report. J Neurooncol. 2007, 85: 77-79. 10.1007/s11060-007-9384-5.

    Article  PubMed  Google Scholar 

  39. Einhorn LH: The case against prophylactic cranial irradiation in limited small cell lung cancer. Semin Radiat Oncol. 1995, 5 (1): 57-60. 10.1016/S1053-4296(05)80011-6.

    Article  PubMed  Google Scholar 

  40. Pedersen AG, Kristjansen PEG, Hansen HH: Prophylactic cranial irradiation and small cell lung cancer. Cancer Treatment Reviews. 1988, 15: 85-103. 10.1016/0305-7372(88)90019-9.

    Article  CAS  PubMed  Google Scholar 

  41. Rihs TA, Michel CM, Thut G: A bias for posterior alpha-band power suppression versus enhancement during shifting versus maintenance of spatial attention. NeuroImage. 2009, 44: 190-199. 10.1016/j.neuroimage.2008.08.022.

    Article  PubMed  Google Scholar 

  42. Romei V, Brodbeck V, Michel C, Amedi A, Pascual-Leone A, Thut G: Spontaneous fluctuations in posterior alpha-band EEG activity reflect variability in excitability of human visual areas. Cereb Cortex. 2008, 18: 2010-2018.

    Article  PubMed  Google Scholar 

  43. Romei V, Rihs T, Brodbeck V, Thut G: Resting electroencephalogram alpha power over posterior sites indexes baseline visual cortex excitability. NeuroReport. 2008, 19: 203-208. 10.1097/WNR.0b013e3282f454c4.

    Article  PubMed  Google Scholar 

  44. Voytek B, Davis M, Yago E, Barcelo F, Vogel EK, Knight RT: Dynamic neuroplasticity after human prefrontal cortex damage. Neuron. 2010, 68 (3): 401-408. 10.1016/j.neuron.2010.09.018.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Schrappe M, Stanulla M: Current Treatment Approaches in Childhood Acute Lymphoblastic Leukemia. 2010, International Society of Paediatric Oncology Education Bookp, 25-38. []

    Google Scholar 

  46. Pergener T: What’s wrong with Bonferroni adjustments. BMJ. 1998, 316: 1236-1238. 10.1136/bmj.316.7139.1236.

    Article  Google Scholar 

Pre-publication history

Download references


We thank H.J.H. van der Pal, R. Heinen, M. van Overveld (Academic Medical Center, Amsterdam), C.J. Stam, G.W.G.A. de Vos, J.P.A. Verbunt, H.J.A. Smaling, P. J. Ris, N.C.G. Sijsma and A.M.E. Linger (VU University Medical Center, Amsterdam) for their contributions.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Ilse Schuitema.

Additional information

Competing interests

The company Sonares BV commercially distributes the ANT program that has been used in this study. Dr. LMJ de Sonneville is director of this firm. The other authors have nothing to declare.

Authors’ contributions

All authors participated in the discussion and design of the study. LS, AV and CB applied for the grant. Data collection and analysis was done by MD, IS en BD. MD and IS wrote the manuscript. All authors read and approved the final manuscript.

Marita Daams, Ilse Schuitema contributed equally to this work.

Electronic supplementary material


Additional file 1: ANTprogram. Description of subtests of the Amsterdam Neuropsychological Tasks program (ANT). (PDF 154 KB)


Additional file 2: Global relative powers. Mean global relative powers per group. Four visualisations of the same data. A: Bar chart of the raw power values, equal to Figure 2 in the main text. B: Bar chart of the log transformed power values. C: Line chart of the raw power values. D: Line chart of the log transformed power values. Error bars represent one standard deviation. (PDF 455 KB)


Additional file 3: Regional relative power differences. Statistics of the multivariate tests and group contrasts of the regional relative powers in each frequency band. (MANOVA 1: theta, MANOVA 2: alpha 1, MANOVA 3: alpha 2, MANOVA 4: beta, MANOVA 5: gamma). Group differences were tested with age as a covariate. (PDF 297 KB)


Additional file 4: Neuropsychological outcome. Overview of group differences on all neuropsychological variables assessed with the Amsterdam Neuropsychological Tasks (ANT) program. Separate ANOVA’s were applied per variable, with simple contrasts between each patient group and controls. (PDF 210 KB)


Additional file 5: Correlations. Overview of differences between correlation coefficients of regional power and ANT variables, calculated separately within the CT + CRT group and in controls. (PDF 557 KB)


Additional file 6: Linear regression models. Overview of linear regression models within the separate (composite) groups. R2 indicates the explained variance of ANT variables by regional powers; p is the significance of each linear regression model. (PDF 190 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Daams, M., Schuitema, I., van Dijk, B.W. et al. Long-term effects of cranial irradiation and intrathecal chemotherapy in treatment of childhood leukemia: a MEG study of power spectrum and correlated cognitive dysfunction. BMC Neurol 12, 84 (2012).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: