Main

Understanding neurodevelopmental alterations in very preterm children (≤32 wk GA) is of increasing importance due to improved survival rates in recent decades, together with the prevalence of cognitive, behavioral, and motor deficits in this group (1) and the increasing incidence of developmental difficulties in this vulnerable population (2). Even in the absence of major neurological impairment and when intelligence is broadly normal, children born very preterm frequently exhibit selective cognitive deficits at school age, such as in executive function and visual processing (37). Children, adolescents, and adults born very preterm have also been shown to activate different neural networks during cognitive and perceptual processing (810). Considerable advances have been made in relating structural brain alterations associated with very preterm birth to developmental outcome (1113). As previous imaging research on brain development in children born very prematurely has predominantly focused on structure, however, very little is known about cortical oscillatory activity in this population. This represents an important question as cortical oscillations play a central role in the generation of cognition and perception (1416) and age-related changes in cortical oscillations reflect cognitive development (17,18). The generation of cortical oscillatory activity relevant to cognition and perception has been shown to depend critically on thalamic input (19).

Several lines of research have indicated altered development of thalamocortical systems in children born very preterm. White matter damage is the prevalent form of brain injury in this vulnerable population, and advanced imaging techniques have identified abnormal white matter connectivity in preterm children who do not exhibit brain injury on conventional MRI (20,21). White matter alterations in very preterm children include abnormal development of thalamocortical connections, even in the absence of focal lesions as identified by conventional MR imaging (22), and when developmental outcome is normal (23). The 24–32 wk gestational period corresponding to very premature birth is also characterized by the prominence of the transient cortical subplate zone, which reaches four to five times the size and thickness of the cortical plate (24), and plays a critical role in the development of thalamocortical connections (25,26). Moreover, investigation of resting state network dynamics using hemodynamic measures of functional connectivity (fcMRI) has demonstrated that very preterm infants scanned at term equivalent age exhibit reduced functional connectivity between cortex and thalamus, relative to full-term controls (27).

Previous magnetoencephalographic (MEG) research has established that disordered thalamocortical interactions are commonly associated with a slowing of the 10 Hz peak cortical oscillation from the alpha-band (8–14 Hz) toward the slower theta (4–7 Hz) frequency range (28). Animal research has demonstrated that thalamic input and thalamocortical interactions are critical for the generation of alpha- and theta-band oscillations in the cerebral cortex and that disruption of thalamic and thalamocortical systems can induce deceleration of cortical alpha rhythms toward the theta-band (2931). In humans, slowing of peak oscillatory frequency has been identified in several neurological and neuropsychiatric disorders including Alzheimer's disease, neurogenic pain, Parkinson's disease, and schizophrenia (28,3236). This oscillatory slowing is often accompanied by reduced alpha-band activity (33,34,37,38), and studies using implanted electrodes and fcMRI have confirmed the involvement of disordered thalamocortical interactions in these pathological conditions (3741). The relationship between the development of corticothalamic connectivity and cortical alpha rhythms is further underlined by recent evidence of correlations between white matter architecture in corticothalamic systems and the parameters of alpha oscillations recorded using EEG (42).

Evidence indicating altered development of thalamocortical systems in children born very prematurely, in light of previous work implicating disordered thalamocortical activity in the slowing and reduction of cortical alpha rhythms, suggests that very preterm children may exhibit a slowing of peak oscillatory frequency and reduced alpha power. To test this hypothesis, we recorded resting neuromagnetic activity from school age children who were born very prematurely without major neurological or intellectual impairment and age-matched full-term controls and analyzed the power spectrum of spontaneous neuromagnetic oscillations.

METHODS

Subjects.

Eleven school age children born very preterm (≤32 wk; mean gestation 30.29 wk; SD = 2.39 wk; range = 26–32 wk), mean age 7.53 y (range = 7.28–7.92 y; SD = 0.18 y) were recruited as part of a longitudinal study on the neurocognitive development of very preterm children (43,44). Eleven age-matched full-term control children, mean age 7.54 y (range = 7.35–7.87 y; SD = 0.15 y), were recruited from a combination of the longitudinal study and from the community at school age. Very preterm children had a mean birth weight of 1471 g (range = 755-2030 g; SD = 527 g). The very preterm group comprised six girls and five boys, and the full-term group consisted of three girls and eight boys. Individual cognitive assessment was carried out using the Wechsler Intelligence Scale for Children (45), with mean full scale intelligence quotient (FSIQ) of 97.18 (15.03) for the very preterm and 108.73 (16.89) for the full-term groups, which did not differ significantly (p = 0.11). None of the children in the very preterm group had significant brain injury (periventricular leukomalacia or grade III-IV intraventricular hemorrhage) evident on neonatal cranial ultrasound (46). No child had any known neurological illness or major sensory, motor, or intellectual impairment. Table 1 displays clinical and birth characteristics of children in our very preterm group. As several of the full-term controls were recruited from outside the longitudinal cohort, characterization of neonatal information was not possible for this group. Informed consent was obtained from each subject and their parent. This study was approved by the Clinical Research Ethics Board of the University of British Columbia and the Research Ethics Board of the Children's & Women's Health Centre of BC, and conforms with the conventions set out in the Declaration of Helsinki.

Table 1 Summary of clinical and birth characteristics of the very preterm group

MEG recording.

Two minutes of spontaneous eyes open data were recorded from each subject using a 151-channel whole-head MEG system (CTF systems; Port Coquitlam, Canada). Children were supine, viewed a “happy face” which was projected onto a screen 40 cm above their eyes, and were monitored by a research assistant in the recording chamber to ensure that subjects maintained open eyes during recording. Data were digitized continuously at 1200 Hz and stored for offline analysis. Eyes closed resting state data were not recorded in addition to eyes open data due to limitations imposed by neuroimaging in special child populations, given that children were also expected to perform a cognitive task in the MEG, described elsewhere (47) and in the discussion. Fiducial coils were attached at the nasion and at left and right preauricular locations, and each coil was energized at a distinct high narrow-band frequency.

MEG analysis.

To standardize head location relative to sensors, between and within subjects, dipolar source solutions were computed for each of the fiducial coils 30 times per second, thereby creating a continuous record of head position during MEG recording. MEG data were then aligned to a common position by performing an inverse solution, data rotation, and forward solution 30 times per second (48). Energized high-frequency activity emitted by the fiducial coils during recording was then removed using notch filtering and data were down sampled to 300 Hz. The record of ocular and nonocular artifacts was removed from MEG recordings using a principal component analysis-based procedure (49). Data were then transformed from axial gradiometer to planar gradiometer sensor space (50) using routines implemented in the FieldTrip software package (http://www.ru.nl/fcdonders/fieldtrip). This method has previously been used to accurately assess the topography of cortical oscillations in sensor-level analyses using 151-channel CTF MEG systems (51). This conversion was performed so that signals correspond more directly to cortical activity directly underlying each sensor, as planar gradiometers exhibit more spatially restricted lead fields than their axial counterparts.

To determine the spectral density of cortical activity, data from each sensor were filtered at 1-Hz intervals, using methods established in previous studies for the accurate measurement of spectral activity (5153), from 6 to 60 Hz (passband equals filtered frequency ±0.05 Hz). This frequency range was chosen as previous results indicate that the slowing of alpha rhythms in neurological and neuropsychiatric conditions produces a peak in the upper theta or lower alpha frequency range (28,32,3638) and to distinguish this peak from lower frequency theta activity which is distinct from the alpha peak. Power was then calculated for each data point in the 120-s filtered time series for each analyzed frequency and sensor and then averaged across all time points. This provided a single value for each subject representing the average power of neuromagnetic oscillations at each frequency and sensor during the 2-min recording section. To examine the topography of oscillatory changes, data were averaged across sensors in within-regional sensor groupings and midline electrodes were excluded (Fig. 1). These regional sensor groupings were selected because they i) have been successfully used to reveal alterations in spontaneous oscillatory activity in clinical populations using a 151-channel CTF system (54,55), ii) conform to regional groupings designated by the system manufacturer, iii) correspond roughly to underlying cortical regions, and iv) allow exploration of topography while limiting the number of statistical comparisons. This averaging produced power values at each analyzed frequency within each regional sensor grouping for each subject. Alpha power was calculated by averaging power across 8–14 Hz values for each sensor grouping for each subject, as this frequency range encapsulates the alpha-band. One-tailed t tests were used to analyze group differences in peak oscillatory frequency and oscillatory power within each regional sensor grouping, as a priori hypotheses regarding the direction of effects were based on oscillatory changes in other special populations with analogous physiological alterations.

Figure 1
figure 1

The sensor montage and regional groupings. Gray circles represent sensors used in the analysis of resting MEG by cortical sector, black circles denote midline sensors excluded from regional analyses, and black lines denote the boundaries of cortical regions in the sensor level analysis. Letters indicate cortical regions: L, left; R, right; F, frontal; C, central; T, temporal; P, parietal; O, occipital.

RESULTS

MEG results.

Inspection of global spectral activity (averaged across all 151 sensors) indicated a reduction of the alpha peak, together with a slowing of peak frequency, in very preterm children, relative to full-term controls (Fig. 2). To investigate the topography of observed alterations in the spectral density of resting MEG in very preterm children and perform statistical tests for group differences, we investigated group differences in peak oscillatory frequency and alpha power (8–14 Hz) within each regional sensor grouping. This revealed significantly slower peak oscillatory frequency in very preterm children (Fig. 3A) which was concentrated and statistically significant over bilateral frontal cortex (p = 0.039 left; p = 0.0031 right). All analyzed regions trended toward oscillatory slowing except for the left occipital area (Fig. 3A). Alpha-band power in very preterm children was significantly reduced over bilateral frontal (p = 0.047 left; p = 0.019 right) and temporal cortex (p = 0.015 left; p = 0.0028 right) and trended toward reduced activity in all other analyzed regions (Fig. 3B). Additional statistical information regarding peak oscillatory frequency and resting alpha power within each analyzed region is presented in Table 2.

Figure 2
figure 2

Global power spectrum of spontaneous MEG activity for very preterm children and full-term controls averaged over all 151 sensors. The blue line represents the global power spectrum for the full-term controls and the red line represents that of the very preterm children.

Figure 3
figure 3

(A) Shift of peak oscillatory frequency in very preterm children, relative to full-term controls, in each cortical region. The magnitude of frequency shifts is denoted by color (red, −1 to −2 Hz; orange, −0.5 to −1 Hz; yellow, 0 to −0.5 Hz; green, 0 to +0.5 Hz) and statistically significant group differences are labeled (*p < 0.05; **p < 0.005). (B) Reduction of alpha (8–14 Hz) power in children born very preterm, relative to full-term controls, in each cortical region. Colors represent the magnitude of group differences (red, −50 to −75 femptotesla (fT); orange, −25 to −50 fT; and yellow, 0 to −25 fT) and statistically significant differences are indicated (*p < 0.05; **p < 0.02; ¶p < 0.003).

Table 2 Mean and standard deviation of peak oscillatory frequency and alpha (8–14 Hz) power in femptotesla (fT) within each regional sensor grouping for both the preterm and full-term groups

DISCUSSION

We demonstrate slowing of resting peak oscillatory frequency and reduced alpha-band power in very preterm children, providing the first evidence of altered spontaneous neuromagnetic activity in this vulnerable population. Slowing of the alpha peak toward the theta frequency range has been identified in a variety of pathological conditions in adults (28,3236) and is often accompanied by a reduction in alpha-band power (33,34,37,38). Our findings provide a novel demonstration of a mild slowing of peak oscillatory frequency in a clinical population without major neurological or intellectual impairment. The findings of this study add to accumulating evidence that children born very preterm who escape significant neonatal brain injury and have neurocognitive function in the broadly normal range, nonetheless express altered cortical activity relevant to cognition. Such observations are consistent with previous work showing that very preterm infants express altered event-related potentials in response to visual, auditory, and somatosensory stimulation (56,57) and more recent evidence that such responses are altered in school age children born prematurely (58) even in the absence of major impairment (59). The results of this study build upon recent functional imaging results showing altered cortical network dynamics in individuals born very prematurely (810) by demonstrating altered expression of cortical oscillations known to be critical for the generation of cognitive processes (1416) in this population. Using a partially overlapping group of school age children born very prematurely drawn from the same longitudinal cohort (43,44), we previously demonstrated reduced inter-regional alpha-band phase locking during visual short-term memory retention in these children together with increased interhemishperic long-range theta-band synchronization (47), suggesting that reduction and slowing of alpha oscillations may play a critical role in the altered neurocognitive development of very preterm children. This outlook is underscored by findings indicating that alpha rhythms play an important role in cognition and perception (60,61) and age-related changes in alpha responses have been demonstrated during cognitive and perceptual processing in children (18,6264).

During the neonatal period corresponding to premature birth, there is a rapid shift toward relatively more high-frequency EEG activity, which includes decreasing theta oscillations (65), indicating that neural mechanisms underlying cortical rhythms are undergoing a critical developmental phase in the perinatal period corresponding to premature birth (66). Although the density of spectral power recorded from the cortex of very preterm infants using EEG becomes very similar to that of full-term infants at term age (67), a progressive shift toward higher frequency oscillations continues throughout normal childhood development (68). This change proceeds along a posterior-to-anterior axis (69), and progressive increases in the expression of alpha rhythms are observed throughout childhood (68). Accordingly, oscillatory slowing and reduced alpha-band power in very preterm children observed in this study could be conceived as a delay in the development of normal brain maturation. However, the finding that the ratio of low-frequency to high-frequency oscillations is increased in young adults born with extremely low birth weight and very prematurely (70) suggests that slowed oscillatory activity is either permanent or extremely long lasting and is embedded among networks across a wider frequency scale. The relationship between neonatal experience and altered cortical activity in school age children born very preterm is an important question for future research, as it is possible that altered neuromagnetic oscillations in these children are due, at least in part, to regulatory challenges associated with exposure to a stressful extrauterine environment (3,71) during a period wherein neural systems relevant to the generation of cortical oscillatory activity are undergoing major development (65,66,72,73).

CONCLUSION

We demonstrate that the peak oscillatory frequency of resting cortical oscillations is slowed in school age children born very preterm, even in the absence of major neurological or intellectual impairment. This constitutes the first evidence of altered spontaneous neuromagnetic activity in this vulnerable population. Slowing of the alpha rhythms toward the theta-band has been previously identified in several neurological and neuropsychiatric populations and has been attributed to disordered thalamocortical dynamics. The results of this study suggest that altered thalamocortical development in very preterm children may contribute to altered expression of oscillatory brain activity in this vulnerable population.