Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Neonatal Pain-Related Stress Predicts Cortical Thickness at Age 7 Years in Children Born Very Preterm

  • Manon Ranger,

    Affiliations Pediatrics, University of British Columbia, Vancouver, British Columbia, Canada, Developmental Neurosciences and Child Health, Child and Family Research Institute, Vancouver, British Columbia, Canada

  • Cecil M. Y. Chau,

    Affiliations Developmental Neurosciences and Child Health, Child and Family Research Institute, Vancouver, British Columbia, Canada, BC Mental Health and Addictions Research Institute, Vancouver, British Columbia, Canada

  • Amanmeet Garg,

    Affiliation Engineering Science, Simon Fraser University, Burnaby, British Columbia, Canada

  • Todd S. Woodward,

    Affiliations BC Mental Health and Addictions Research Institute, Vancouver, British Columbia, Canada, Psychiatry, University of British Columbia, Vancouver, British Columbia, Canada

  • Mirza Faisal Beg,

    Affiliation Engineering Science, Simon Fraser University, Burnaby, British Columbia, Canada

  • Bruce Bjornson,

    Affiliations Pediatrics, University of British Columbia, Vancouver, British Columbia, Canada, Developmental Neurosciences and Child Health, Child and Family Research Institute, Vancouver, British Columbia, Canada

  • Kenneth Poskitt,

    Affiliations Developmental Neurosciences and Child Health, Child and Family Research Institute, Vancouver, British Columbia, Canada, Radiology, University of British Columbia, Vancouver, British Columbia, Canada

  • Kevin Fitzpatrick,

    Affiliation Developmental Neurosciences and Child Health, Child and Family Research Institute, Vancouver, British Columbia, Canada

  • Anne R. Synnes,

    Affiliations Pediatrics, University of British Columbia, Vancouver, British Columbia, Canada, Developmental Neurosciences and Child Health, Child and Family Research Institute, Vancouver, British Columbia, Canada, BC Children’s and Women’s Hospitals, Vancouver, British Columbia, Canada

  • Steven P. Miller,

    Affiliations Pediatrics, University of British Columbia, Vancouver, British Columbia, Canada, Developmental Neurosciences and Child Health, Child and Family Research Institute, Vancouver, British Columbia, Canada, Pediatrics, Hospital for Sick Children and University of Toronto, Toronto, Ontario, Canada

  • Ruth E. Grunau

    rgrunau@cw.bc.ca

    Affiliations Pediatrics, University of British Columbia, Vancouver, British Columbia, Canada, Developmental Neurosciences and Child Health, Child and Family Research Institute, Vancouver, British Columbia, Canada, BC Children’s and Women’s Hospitals, Vancouver, British Columbia, Canada

Abstract

Background

Altered brain development is evident in children born very preterm (24–32 weeks gestational age), including reduction in gray and white matter volumes, and thinner cortex, from infancy to adolescence compared to term-born peers. However, many questions remain regarding the etiology. Infants born very preterm are exposed to repeated procedural pain-related stress during a period of very rapid brain development. In this vulnerable population, we have previously found that neonatal pain-related stress is associated with atypical brain development from birth to term-equivalent age. Our present aim was to evaluate whether neonatal pain-related stress (adjusted for clinical confounders of prematurity) is associated with altered cortical thickness in very preterm children at school age.

Methods

42 right-handed children born very preterm (24–32 weeks gestational age) followed longitudinally from birth underwent 3-D T1 MRI neuroimaging at mean age 7.9 yrs. Children with severe brain injury and major motor/sensory/cognitive impairment were excluded. Regional cortical thickness was calculated using custom developed software utilizing FreeSurfer segmentation data. The association between neonatal pain-related stress (defined as the number of skin-breaking procedures) accounting for clinical confounders (gestational age, illness severity, infection, mechanical ventilation, surgeries, and morphine exposure), was examined in relation to cortical thickness using constrained principal component analysis followed by generalized linear modeling.

Results

After correcting for multiple comparisons and adjusting for neonatal clinical factors, greater neonatal pain-related stress was associated with significantly thinner cortex in 21/66 cerebral regions (p-values ranged from 0.00001 to 0.014), predominately in the frontal and parietal lobes.

Conclusions

In very preterm children without major sensory, motor or cognitive impairments, neonatal pain-related stress appears to be associated with thinner cortex in multiple regions at school age, independent of other neonatal risk factors.

Introduction

Altered brain development is evident in children born very preterm (24–32 weeks gestational age) in early infancy [1][3] and at school-age [4], [5], including reduced gray and white matter volumes in infancy, childhood and adolescence [6][10], compared to term-born peers. Normal brain growth involves changes in cortical thickness reflecting cellular maturational changes related to myelination and synaptic pruning [11]. Recent findings from our group, in a separate cohort, showed delayed microstructural development of the cortical gray matter in very preterm neonates at term equivalent [12]. Children born preterm show altered cortical thickness in childhood and adolescence [7], [13][18]. Specifically, thinner cortex has been reported in superior and temporal, middle frontal, anterior cingulate cortex, supramaginal, precuneus, and post central regions when compared to term born controls at school age [7], [14]. In adolescents born preterm, thinner cortex has been reported in the enthorhinal, temporal, and parietal regions [13], [16], [18]. Thicker cortex has been reported in preterm children with periventricular leukomalacia, identified on neonatal ultrasound and still evident on magnetic resonance (MR) images at school age [14]. Little is known about the etiology and possible clinical risk factors underlying regional cortical thinning in children born preterm who do not have major brain injury.

Early-life adversity in animal models, including maternal separation [19][22] and pain exposure [23][27], induces long-term changes in brain and behavior (reviews [28], [29]). Infants born very preterm are exposed to weeks to months of hospitalization in the neonatal intensive care unit (NICU) during a period of rapid neuronal proliferation and cell differentiation (e.g. differentiation of subplate neurons, formation of synapses, selective pruning of neuronal processes and synapses). Repeated procedures inducing pain-related stress place very preterm infants at particular risk due to the very rapid brain development and programming of stress systems while they are in the NICU [30], [31]. Given our previous findings that pain-related stress was associated with atypical brain development in the neonatal period [8], our aim was to evaluate whether neonatal pain-related stress, quantified as the number of skin-breaking procedures adjusted for confounding clinical factors, is associated with variations in cortical thickness at age 7 years in a different cohort of children born very preterm.

We hypothesized that greater neonatal pain-related stress would be related to reduced cortical thickness, after adjusting for confounding neonatal clinical factors associated with premature birth and neonatal care such as gestational age, number of days on mechanical ventilation, early severity of illness, surgery, infection, and morphine exposure.

Materials and Methods

1. Study Design and Participants

Participants were part of a larger longitudinal study of long-term effects of neonatal pain-related stress on neurodevelopment of children born very preterm (24–32 weeks gestation) e.g. [31], [32], who were admitted to the level III NICU at British Columbia’s Women’s Hospital between 2000 and 2004. Out of 106 children seen at age 7 years, 43 refused magnetic resonance imaging (MRI), leaving 63 school-age children, 2 were not scanned, 61 underwent MRI. From the 61 preterm children that underwent MRI scans at age 7 years, 13 were excluded due to poor quality MR images due to movement artifact and one for missing neonatal clinical data. In addition, two children with periventricular leukomalacia (PVL) and/or intraventricular hemorrhage (IVH) grade 3 or 4 on neonatal ultrasound and confirmed on MR at age 7 years (author KP) were excluded, and three left handed children were excluded to eliminate the effect of handedness on cortical thickness asymmetries [33], [34]. We included in our study three very preterm infants diagnosed on neonatal ultrasound with IVH grade 1 or 2 (two infants with IVH grade 1) and six children that showed minimal to moderate white matter injury (≤3 lesions) on MR scans at school-age, one of which had IVH grade 1 on neonatal ultrasound. None of the cortical thickness measures for these children were outliers compared to the rest of the sample, therefore they were included in the study. All of the children scanned had an IQ score above 70 on the Wechsler Intelligence Scale for Children –4th Ed [35], and none had a major sensory or motor impairment. The final study sample comprised 42 right-handed children (38% boys).

The study was approved by the Clinical Research Ethics Board of the University of British Columbia and the British Columbia Children’s and Women’s Research Ethics Board. Written informed consent was obtained from parents and assent from children.

2. Procedures

Magnetic resonance imaging.

MRI was performed using a standard 12 channel head coil on a Siemens 1.5 Tesla Avanto (Berlin, Germany) with VB 16 software. The following images were acquired: a 3D T1 weighted SPGR sequence 18/9.2/256/1 mm/0/256×256 (TR/TE/FOV/Thickness/Gap/Matrix), axial FSE T2 4030/90/220/3 mm/0.1 mm/512×354, axial FLAIR 8900/87/5 mm/1 mm/256×154 and a 12 direction DTI sequence 7800/82/256/2 mm/0/128×128 using B values of 700 and 1000. All imaging sessions were performed without sedation. On the study day, each child first had a session in a mock scanner to acclimatize to the noise and feeling of being in a MRI scanner, followed by the actual study scan. Children were instructed to remain still and watched a video during the sessions that lasted approximately 30 minutes.

An experienced pediatric neuroradiologist (KP), blinded to the child’s medical history, assessed the MR scans for ventriculomegaly, cerebellar hemorrhage and severity of white matter injury, as previously described [36]. No child had a severe brain injury at school age (i.e. no cerebellar hemorrhage, ventriculomegaly or severe white matter injury [i.e. >3 lesions or 2 with 5% hemisphere involved]).

Cortical thickness measures.

Regional cortical thickness was calculated using a topology aware technique originally proposed by Gibson et al. [37] which modifies the normally accepted method of calculating the thickness as the length of the streamlines obtained from the well-established and validated Laplacian streamlines method [38], [39]. This implementation uses FreeSurfer surfaces to enforce separation of volumetric domains, which is important where adjacent gyri/sulci may abut. In the original implementation, the abutting gyri/sulci can lead to overestimation of cortical thickness, which is avoided in this custom implementation, thereby giving more accurate thickness in regions where topological issues such as touching gyri occur, and no difference in thickness calculated in other regions. For detecting group-wise maps of statistically significant thickness changes, FreeSurfer and Laplacian-streamlines based volumetric methods lead to similar results [40]. The coupled surface method establishes the inner surface initially and propagates that mesh to another surface thus preserving the corresponding points incorporating the topology dictated by the inner and the outer surfaces of the gray matter in the volumetric Laplacian streamline thickness computation as described in [41] as the arc length of the streamlines of the solution of Laplace’s equation on the cortical mantle.

Cortical thickness was calculated for each individual brain using the segmentation output from the preprocessing step, where the output obtained is the pial surface for each subject where each vertex has a thickness value associated with it. The surfaces used to define the outer and inner boundaries in our methods were constructed using an automated cortical surface reconstruction method developed by Dale and Fischl [42], [43]. The 66 brain regions of interest (ROIs) used in this study were parcellated by FreeSurfer using the Desikan-Killiany Atlas, an automated labeling system developed and validated by Desikan et al. [44].

Clinical data collection.

Medical and nursing chart review of neonatal data from birth to term equivalent was carried out by highly trained neonatal research nurses. Data collected included, but was not limited to, birth weight, gestational age (GA), number of days on mechanical ventilation and/or oscillation, illness severity on day 1 (Score for Neonatal Acute Physiology [SNAP]- II [45]), number of surgeries, presence of culture proven infection, and cumulative dose of morphine. The cumulative dose of morphine was calculated (intravenous dose plus converted oral dose) as the average daily dose adjusted for daily body weight, multiplied by the number of days the drug was given, as we have used previously [8], [31]. We quantified neonatal pain-related stress as the number of skin-breaking procedures (e.g., heel lance, peripheral intravenous or central line insertion, chest-tube insertion, tape removal, and nasogastric tube insertion) during the stay in the NICU, as previously used [8], [31], [46]. Each attempt at a procedure was counted as one skin-break; all nursing staff in our NICU have been trained to precisely record each attempt.

3. Statistical Analysis

Principal component analysis (PCA) is a method of data reduction of large sets of related variables by reducing them to a few vectors of weightings that best explain the variance, while losing as little information as possible. Each extracted vector, referred to as a “component”, accounts for a portion of the total variance in the data: the first component accounts for the largest amount of variance, with each successive component accounting for a smaller amount of the total variance. Constrained principal component analysis (CPCA) combines multivariate multiple regression and PCA, and allows examination of the component structure of the variance in a set of dependent variables that is specifically predicted by a set of predictor variables [47][49]. CPCA is two-step process, referred to as the external and internal analysis. The external analysis consists of a multivariate least squares multiple regression of the dependent measures on the independent measures, producing predicted and residual scores for each dependent measure. In the present study, the matrix of predicted scores reflects the variance in cortical thickness that can be accounted for by the neonatal clinical variables, and the residual matrix reflects the variance that cannot be accounted for by the neonatal clinical variables. Seven neonatal clinical predictors were included: number of skin-breaking procedures, gestational age, number of days on mechanical ventilation, illness severity on day 1 (SNAP-II), number of surgeries, culture proven infection, and cumulative morphine exposure adjusted for daily weight. The internal analysis consists of PCAs on each of the aforementioned matrices. The resulting component solutions (overall, predicted, and residual solutions) are examined to determine which dimensions of the cortical thickness can be explained by the neonatal clinical variables. To determine the particular neonatal clinical predictors that are related to each of the components extracted from the predicted solution, correlations were computed between each component score and the set of independent variables (i.e., neonatal clinical data). The number of components retained for each PCA was determined by inspection of scree plots. All PCA solutions were separately rotated using varimax with Kaiser normalization. Solutions were bootstrapped 1000 times by Monte Carlo methods to produce confidence intervals and p-values [50], [51]. Computations for CPCA and bootstrapping were carried out using MATLAB version 7.6 (The MathWorks, 2008, Natick, MA).

Then, to further examine the individual contribution of each neonatal predictor on the measured variations of cortical thickness, we conducted a generalized linear (GENLIN) modeling analysis, adjusting for neonatal clinical factors (GA, SNAP-II at day 1, culture proven infection, number of days on mechanical ventilation, number of surgeries, cumulative morphine exposure adjusted for daily weight, and number of skin-breaking procedures). GENLIN provides an extension of general linear models and relaxes the requirement of equality or constancy of variances that is required in traditional linear models. GENLIN was carried out using the Statistical Package for Social Sciences (SPSS) version 16.0 (IBM, Somers, NY). False discovery rate (FDR 5%) [52] was employed to correct for multiple comparisons between brain regions for CPCA and GENLIN.

Results

Demographic and neonatal clinical data are presented in Table 1. Children were scanned at a median age of 7.8 years (interquartile range [IQR] 7.7–8). During their NICU stay, these children underwent a median of 74 (IQR 45–136) skin-breaking procedures.

Since child age may contribute to variations in cortical development [11], we performed correlations between age at MR scan and cortical thickness and no statistically significant associations were found, after Bonferroni correction for multiple comparisons. Given that gender differences in cortical thickness have been reported in adolescence [53], we conducted independent t-tests comparing cortical thickness in boys and girls, finding no statistically significant differences in any of the 66 brain regions assessed, after Bonferroni correction for multiple comparisons. It was noted that without correction, cortical thickness differed by gender only in 1/66 regions.

1. Constrained Principal Component Analysis (CPCA)

Given the limited number of participants, the high number of brain regions examined, and correlated neonatal clinical factors, we first examined the relationships between the neonatal clinical variables and cortical thickness in all brain regions. This permitted unbiased handling and interpretation of our data.

Table 2 shows the distribution of variance of the overall, predicted and residual solutions of the CPCA. Three components were extracted from the predicted solutions (eigenvalues 14.77, 2.46, and 1.56 respectively), which corresponds to 13.4% of the overall variance for component 1, 7.5% and 6.9% for components 2 and 3 respectively. The external analysis showed that the 7 clinical predictors accounted for 32.9% of the overall variance in cortical thickness.

thumbnail
Table 2. Constrained principal component analysis results to explain variance in cortical thickness in relation to neonatal clinical factors.

https://doi.org/10.1371/journal.pone.0076702.t002

Component loadings and their bootstrapped confidence intervals and p-values in each brain regions of interest (ROIs) are summarized in Table 3. Both the first and second components were characterized by negative loadings on all brain regions (i.e. decreased cortical thickness), while the third component had positive loadings (i.e. increased cortical thickness). The dominant loadings on the first component were distributed prominently on various frontal and central brain regions. On the second component, left supramarginal, bilateral inferior parietal, and right lateral occipital ROIs had the largest loadings. Bilateral caudal anterior cingulate and left lingual ROIs had the largest loadings for the third component.

thumbnail
Table 3. Brain regions with loadings for the 3 components.

https://doi.org/10.1371/journal.pone.0076702.t003

The variance in cortical thickness explained by the 7 neonatal clinical factors within each component is shown in Table 4. The external analysis in CPCA showed that pain-related stress was the strongest predictor from our set of 7 neonatal factors. The number of skin-breaking procedure loaded highly on the first two components, infection loaded significantly only on component 1, while surgery significantly loaded uniquely and predominately on component 2. Morphine exposure and duration of mechanical ventilation significantly loaded on all three components but dominated the 3rd component. Severity of illness on day 1 (SNAP-II) loaded significantly on all 3 components. Component 1 suggests that the number of skin-breaking procedures during NICU care, illness severity on day 1, culture-proven infection, gestational age, cumulative morphine exposure, and days on mechanical ventilation are related to reduced cortical thickness, primarily in frontal and parietal regions. Component 2 in addition, captured predominately surgery, but not infection, and is related to reduced cortical thickness, primarily in parietal, temporal and occipital regions. Component 3, contrary to the two other components, was positively associated with cortical thickness and primarily reflected cumulative morphine exposure and number of days on mechanical ventilation. Therefore, we interpreted component 1 as procedural and inflammation related pain, component 2 as surgery related pain, and component 3 as primarily morphine exposure. Since both components 1 and 2 had negative loadings on all brain regions, this implied that higher procedural/inflammatory and surgical pain were associated with a reduction in cortical thickness predominately in the frontal and parietal regions for component 1, and parietal, temporal and occipital regions for component 2. In contrast, component 3 showed a positive association with cortical thickness, predominately within the anterior cingulate cortex ROIs, suggesting thicker cortex in association with higher cumulative morphine exposure in these regions.

thumbnail
Table 4. Constrained principal component analysis loadings for the 7 neonatal clinical factors to explain variance in cortical thickness.

https://doi.org/10.1371/journal.pone.0076702.t004

To check whether birth weight (BW) might provide important additional information beyond GA in explaining variations in cortical thickness, we added BW to the set of neonatal factors in a further CPCA analysis. With the inclusion of BW and GA, the overall predictable variance accounted for by the neonatal factors increased by only 1.8%. Given the low sample size and the high correlation between GA and BW (r = 0.76, p<0.0001) we only retained GA in further analyses.

2. Generalized Linear Modeling (GENLIN)

Neonatal clinical factors were inspected for normality, then when necessary were log transformed (neonatal skin-breaking procedures and number of days on mechanical ventilation). Using the findings from the CPCA, a GENLIN model for each of the 66 brain regions was conducted, corrected for multiple comparisons (i.e. 5% FDR). After adjusting for neonatal clinical factors (i.e. gestational age, SNAP-II at day 1, infection, number of days on mechanical ventilation, number of surgery, cumulative morphine exposure), greater number of neonatal skin-breaking procedures was significantly associated with reduced cortical thickness in 21 out of the 66 brain regions assessed (p-values ranged from 0.00001 to 0.014; Table 5), p-value threshold for significance adjusted for multiple comparisons with a FDR correction set at 5%. Most significant (p≤0.001) thinning in relation to neonatal pain-related stress was found in bilateral postcentral, superior frontal, rostral middle frontal, left hemisphere precentral and pars orbitalis, as well as the right hemisphere supramarginal region. Of all the neonatal clinical factors, pain-related stress (adjusted for neonatal clinical factors) was the most robust independent predictor of regional variation in cortical thickness. The association between neonatal pain-related stress (adjusted for confounders) and right hemisphere postcentral and left pars orbitalis cortical thickness are shown as examples (Figure 1).

thumbnail
Figure 1. Cortical thickness at age 7 years in relation to neonatal pain-related stress adjusted for clinical confounders.

Cortical thickness in two of the brain regions where most significant (p≤0.0001) thinning in relation to neonatal pain-related stress was found. A) Scatter plot of right hemisphere postcentral cortical thickness in relation to neonatal pain-related stress (number of skin-breaking procedures log transformed) adjusted for gestational age, severity of illness on day 1, number of culture proven infection, number of days on mechanical ventilation, number of surgeries, and cumulative daily morphine dose in milligrams adjusted for daily body weight. B) Scatter plot of left hemisphere pars orbitalis cortical thickness in relation to neonatal pain-related stress (same as in A).

https://doi.org/10.1371/journal.pone.0076702.g001

thumbnail
Table 5. GENLIN results for the 21/66 cortical regions significantly associated with neonatal pain-related stress (adjusted for clinical confounders).

https://doi.org/10.1371/journal.pone.0076702.t005

After adjusting for neonatal clinical factors including neonatal pain-related stress, and correcting for multiple comparisons (FDR 5%), greater morphine exposure was associated with increased cortical thickness in bilateral caudal anterior cingulate (p = 0.000 right; p = 0.003 left), right posterior cingulate (p = 0.002), as well as in the left lingual (p = 0.002) and middle temporal regions (p = 0.002) (Table 6). In our NICU, very preterm infants only received morphine if they underwent mechanical ventilation (n = 25). To confirm these GENLIN results, we reran the analysis on the 25 infants exposed to mechanical ventilation. Morphine remained significant only in the right caudal anterior and posterior cingulate cortices in the subset of ventilated very preterm children.

thumbnail
Table 6. GENLIN results for the 5/66 cortical regions significantly associated with neonatal cumulative morphine exposure (adjusted for clinical confounders).

https://doi.org/10.1371/journal.pone.0076702.t006

Discussion

Our primary finding was that greater neonatal pain-related stress (adjusted for multiple clinical factors related to preterm birth and NICU care) was associated with lower cortical thickness at school age in children born very preterm. Specifically, in a conservative two-step approach, first using constrained principal component analysis followed by generalized linear modeling to control for clinical confounders related to prematurity, we found that pain-related stress was the strongest predictor of the variation in overall cortical thickness. Precisely, after applying a 5% false discovery rate to control for overall Type-I error, and adjusting for neonatal clinical factors, greater exposure to neonatal skin-breaking procedures was the predominant consistent neonatal factor related to thinner cortex in 21/66 areas, specifically in frontal, parietal, and temporal regions, independent of the neonatal confounding factors.

Consistent with our present study, we recently found in a different cohort of preterm infants, after similarly controlling for multiple confounding neonatal clinical factors, that greater exposure to neonatal pain-related stress was associated with altered brain microstructural development (subcortical gray matter and white matter) from early in life to term equivalent age [8]. It is known that during the last trimester of gestation and early postnatally, major axonal development takes place in the cerebrum, and damage to white matter tracts (i.e. decreased myelination) and subcortical structures leads to detrimental effects on neuronal migration and cortical development [30], [54]. Thus, in two independent cohorts, our work provides converging evidence that higher exposure to neonatal pain-related stress may contribute to altered brain development early in life and later at school-age in children born very preterm.

Prematurity has been associated with short and long-term abnormalities in cortical thickness compared to healthy term born controls [7], [13][18]. Lax et al. [7] showed lower overall mean cortical thickness in 25 very preterm 9 year-old children compared to 32 full-term controls, with significant cortical thinning in bilateral precuneus and right supplemental motor cortex, anterior insula, superior temporal, and postcentral regions. Compared to 22 control children born at full-term, Zubiaurre-Elorza et al. [14] found reduced cortical thickness in left middle and superior frontal regions, supramarginal, and post central gyri in 14 preterm 9 year-old children with no evidence of PVL. In the present study, we found lower cortical thickness was related to greater neonatal pain exposure in these same regions among children born very preterm. Thus it appears we have identified a potential factor that may contribute to the etiology of reduced cortical thickness in very preterm children. Moreover, our current findings of an association between neonatal pain-related stress (adjusted for clinical confounders) and cortical thickness are consistent with our previous reports on short-term [8] and long-term [55] brain development. Importantly we excluded children with severe brain injury (i.e. PVL and/or IVH grade 3–4) in the newborn period, and/or major neurodevelopmental, motor or sensory impairments, therefore our findings are not explained by severe brain injury. Testosterone levels play a role in cortical thickness in adolescence [53]. In the present study, in 7 year olds born very preterm we did not find differences in cortical thickness between boys and girls in preliminary analyses. Gender was not included as a predictor in our statistical modeling due to limited sample size, and remains for future studies in preterm pre-adolescents.

Thinner cortex in specific regions in our sample of very preterm children may be related to delayed neuronal or glial maturation or cell death, or both. Normal cortical maturation shows an increase in cortical thickness during childhood (up to age 11 years) in most of the lateral frontal, temporal, parietal and occipital regions, a reduction during adolescence, followed by stabilization in adulthood [56]. Hence it is possible that normally occurring cortical growth is delayed in children born very preterm. Recent findings in a premature fetal sheep model showed the importance of cell dysmaturation rather than cell loss (especially neuronal) during cortical growth and microstructure development [54]. The potential to normalize the maturation of these cells remains unknown. Although mechanisms underlying cortical thinning in preterm school-age children remain unclear, findings in rat pup pain models [57][59] stress the long-term damaging effects of repetitive neonatal pain exposure on cortical and subcortical neuronal development. However, since we did not compare cortical thickness of very preterm children with a term born control group, it is not possible at this time to conclude that relative reduction in cortical thickness among the very preterm children reflects delayed or accelerated cortical development.

Moreover, volumetric MRI and diffusion tensor imaging (DTI) analyses have demonstrated decreased neuronal volumes (e.g. in cerebellum, thalamus, cerebral cortex) and possible axonal disturbances in very low birth-weight neonates from term-equivalent age to later in adulthood [30]. Sensorimotor, premotor, temporal and parieto-occipital regions are generally reported as the cortical regions showing the most significant volumetric reductions in preterm born school-age children and adolescents [6], [9], [60][62]. Smith and colleagues [2] found that greater exposure to stressful procedures (e.g. heel lance/venipuncture, intubation/extubation, diaper change) in the NICU was associated with reduced brain size in the frontal and parietal regions in preterm neonates assessed at term equivalent age. In addition, in that study, functional connectivity MRI and DTI measures showed that alteration in brain microstructure and functional connectivity within the temporal lobes were related to greater stress exposure. Thus, signifying of the importance of early adverse stressful and painful experience converge.

In response to heel lance, Slater et al. [63] showed that greater neuronal activity was evident in premature neonates compared to age-matched term-born counterparts. This increased neuronal excitation could be detrimental to the preterm neonate’s immature and rapidly developing neural circuitry by altering apoptosis (programmed cell death) and neuronal survival [64]. However, our understanding of how neonatal pain and stress exposure impacts neurodevelopment in these immature infants is still emerging and basic animal research is crucial to determine mechanisms.

To our knowledge, only a few studies in the neonatal rodent model have specifically examined pain and altered brain development. Early inflammatory pain in neonatal rodent pups induced increased cortical and subcortical neuronal activation [59], increased hippocampal gene expression [58], and widespread cell death [57], [59], thus modifying both the structure and function of the developing brain [57]. Persistent inflammatory pain (formalin injections) or repetitive pain (saline injections) may induce major neuronal apoptosis and altered expression of neurodevelopmentally important proteins in the rat pup brain during the first week of life [57]. The most degenerative cells were found in the lamina II of both frontal and parietal cortex. However, there were differences in long-term effects on the adult brain depending on the timing and type of neonatal pain [57]. Given that the majority of the neonatal pain-related stress procedures occur within the first weeks of the NICU stay, especially in those who undergo the most procedures [8], addressing the management of procedural pain-related stress to protect the developing brain in very preterm infants is pressing. Moreover, opioids seem to have differing effects in the presence or absence of pain. In the Duhrsen et al. study [57] and in previous findings in neonatal rat models [23], [65], pre-emptive morphine only protected the neonatal pup against long-term brain and/or behavioral changes when inflammatory pain was present. In the present study, neonatal cumulative morphine exposure was not related to variation in cortical thickness in 20 of the 21 brain regions significantly adversely associated with pain-related stress exposure. In contrast, greater neonatal morphine was associated with thicker cortex in the cingulate cortices, regions known to be involved in mediating opioid analgesia [66][69], as well as in processing emotional aspects of pain [66], [70][72]. However, this relationship remained significant only in 2/66 regions when examined in the subset of children who had been exposed to mechanical ventilation and to morphine as a neonate. Thus this finding may reflect Type-I error.

Routine treatment of pain with morphine is no longer advocated in ventilated preterm infants in the NICU [73], since in the short-term continuous morphine infusions lead to longer duration of mechanical ventilation and longer time to reach full enteral feeding [74], and concerns remain about long-term effects of analgesia and sedation on the developing brain [73][78]. Our present findings of an association between higher exposure to neonatal morphine and greater cortical thickness in the cingulate regions should be considered in the context of the lack of significant effect in the rest of the cortex. It remains for future research to verify whether these findings imply delayed or accelerated cortical development in very preterm infants at school age. Randomized controlled clinical trials of short [74], [75] and long-term effects [76], [79], [80] of pre-emptive morphine infusion for mechanically ventilated neonates have shown no beneficial effects in the short-term outcomes [74], [75]. Longer-term cognitive/motor and behavioral outcomes in relation to morphine exposure have reported mixed results [76], [79], with the most recent findings suggesting a possible neuroprotective effect of pre-emptive morphine infusion on some aspects of parent report of executive functions at 8–9 years [80]. However, in that study, morphine was not found to be beneficial on teacher reports or on individual testing of child executive functions, thus the results of that study [80] remain tentative.

Importantly, in the present study, severity of illness on day 1 (SNAP-II) also contributed to lower cortical thickness at school age. We have recently found, in a separate cohort, that illness severity of the very preterm neonate during the first 24 hours of life together with pain-related stress exposure contributed to slower microstructural development of the corticospinal tract up to term equivalent age [81]. Together, these findings suggest the importance of protecting the developing brain during the first days of life.

During their NICU care, preterm neonates are exposed to multiple factors that may alter the developing brain, and teasing out specific pain-related effects is challenging. Prenatal and post-natal clinical factors and treatments might interact or may lead to similar end points, which makes them difficult to isolate [82]. Nonetheless, one strength of the present study is extensive control for multiple confounding neonatal clinical factors in our statistical modeling. Unique to the present study of this cohort is the use of constrained principal component analysis as a first step which allowed an unbiased examination of the contribution of 7 pre-selected neonatal clinical factors on cortical thickness variation in 66 brain regions of interest. Our group has shown that neonatal pain-related stress remained a robust predictor of white and gray matter deficits in the short-term [8]. On brain imaging at school age in the present study, the number of skin-breaking procedures during NICU care stood out as the strongest predictor of cortical thickness. Interestingly, we found that greater number of skin-breaks was highly related to bilateral thinner cortex of the postcentral region where the somatosensory cortices lie. Given the known plasticity of the developing pain system [64], associated sensitization, and findings of altered fMRI response to pain at school age [83] it appears highly possible that repeated early noxious nociceptive stimulation may cause rewiring of the somatosensory neural circuitry leading to wider spread detrimental effects on cortical development. Further studies are needed to examine cortical thickness alterations and neonatal pain-related stress exposure in relation to pain threshold and specific domains of cognitive and motor functioning at school age.

This study has limitations in that there could be other confounding factors associated with early pain-related stress exposure and/or prematurity, and a relatively limited sample size. Moreover, the lack of a control group of term born children is a limitation of this study. Advances in prenatal ultrasound have permitted reliable dating of GA. Currently GA is viewed as the single best index of physiological maturity in the preterm newborn, and is thus more widely used as a predictor of long-term outcomes than BW [2], [7]. In contrast, BW reflects a combination of immaturity and potential growth retardation. Examining pain-related stress in infants born small for gestational age (i.e. lower BW than expected relative to GA) in relation to long-term brain development remains to be examined in future studies. Important strengths of our study are that we controlled for seven neonatal factors that capture the most salient neonatal clinical features, the strict exclusion criteria, and a false discovery rate of 5%, thus a conservative approach that reduced the chance of finding significant relationships between neonatal pain-related stress and cortical thickness at school age.

Conclusion

In very preterm children without severe brain injury in the newborn period, and free from major sensory, motor or neurodevelopmental impairments, we showed that greater exposure to neonatal pain-related stress was associated with thinner cortex in multiple regions at school age, independent of other neonatal risk factors. Important to further advancing our understanding of the relationship between pain-related stress and alterations in cortical development will be the examination of corticospinal tracts, white matter and sub-cortical gray matter structures in school-age children born very preterm. Finally, it remains for future studies to evaluate to what extent these pain-related stress brain alterations are related to cognitive, motor, and behavioral outcomes in these children.

Acknowledgments

The authors would like to acknowledge the patients and families who generously gave their time for this project. We thank Ivan Cepeda and Gisela Gosse for coordinating the study and Amanda Degenhardt and Katia Jitlina for help in data collection, and Dr. Yoshio Takane for his support on statistical analysis.

Author Contributions

Conceived and designed the experiments: REG SPM. Performed the experiments: CMYC KP BB KF AG MFB. Analyzed the data: MR CMYC TSW REG. Contributed reagents/materials/analysis tools: MFB AG. Wrote the paper: MR REG CMYC. Contributed to the interpretation of the data: ARS SPM REG. Provided critical review of the manuscript: ARS SPM REG. Approved the final manuscript: MR REG CMYC AG TSW MFB BB KP KF ARS SPM.

References

  1. 1. Rogers CE, Anderson PJ, Thompson DK, Kidokoro H, Wallendorf M, et al. (2012) Regional Cerebral Development at Term Relates to School-Age Social–Emotional Development in Very Preterm Children. J Am Acad Child Adolesc Psychiatry 51: 181–191.
  2. 2. Smith GC, Gutovich J, Smyser C, Pineda R, Newnham C, et al. (2011) Neonatal intensive care unit stress is associated with brain development in preterm infants. Ann Neurol 70: 541–549.
  3. 3. Miller SP, Vigneron DB, Henry RG, Bohland MA, Ceppi-Cozzio C, et al. (2002) Serial quantitative diffusion tensor MRI of the premature brain: development in newborns with and without injury. J Magn Reson Imaging 16: 621–632.
  4. 4. Doesburg SM, Ribary U, Herdman AT, Cheung T, Moiseev A, et al. (2010) Altered Long-Range Phase Synchronization and Cortical Activation in Children Born Very Preterm. IFMBE Proc 29: 250–253.
  5. 5. Doesburg SM, Ribary U, Herdman AT, Moiseev A, Cheung T, et al. (2011) Magnetoencephalography reveals slowing of resting peak oscillatory frequency in children born very preterm. Pediatr Res 70: 171–175.
  6. 6. Nosarti C, Al-Asady MHS, Frangou S, Stewart AL, Rifkin L, et al. (2002) Adolescents who were born very preterm have decreased brain volumes. Brain 125: 1616–1623.
  7. 7. Lax I, Duerden E, Lin S, Mallar Chakravarty M, et al.. (2012) Neuroanatomical consequences of very preterm birth in middle childhood. Brain Struct Funct: 1–11.
  8. 8. Brummelte S, Grunau RE, Chau V, Poskitt KJ, Brant R, et al. (2012) Procedural pain and brain development in premature newborns. Ann Neurol 71: 385–396.
  9. 9. Nosarti C, Giouroukou E, Healy E, Rifkin L, Walshe M, et al. (2008) Grey and white matter distribution in very preterm adolescents mediates neurodevelopmental outcome. Brain 131: 205–217.
  10. 10. Nosarti C, Mechelli A, Herrera A, Walshe M, Shergill SS, et al. (2011) Structural covariance in the cortex of very preterm adolescents: A voxel-based morphometry study. Hum Brain Mapp 32: 1615–1625.
  11. 11. Sowell ER, Thompson PM, Leonard CM, Welcome SE, Kan E, et al. (2004) Longitudinal mapping of cortical thickness and brain growth in normal children. J Neurosci 24: 8223–8231.
  12. 12. Vinall J, Grunau RE, Brant R, Chau V, Poskitt KJ, et al. (2013) Slower postnatal growth is associated with delayed cerebral cortical maturation in preterm newborns. Sci Transl Med 5: 168ra8.
  13. 13. Nagy Z, Lagercrantz H, Hutton C (2011) Effects of preterm birth on cortical thickness measured in adolescence. Cereb Cortex 21: 300–306.
  14. 14. Zubiaurre-Elorza L, Soria-Pastor S, Junque C, Sala-Llonch R, Segarra D, et al. (2012) Cortical thickness and behavior abnormalities in children born preterm. PLoS One 7: e42148.
  15. 15. Phillips JP, Montague EQ, Aragon M, Lowe JR, Schrader RM, et al. (2011) Prematurity Affects Cortical Maturation in Early Childhood. Pediatr Neurol 45: 213–219.
  16. 16. Martinussen M, Fischl B, Larsson HB, Skranes J, Kulseng S, et al. (2005) Cerebral cortex thickness in 15-year-old adolescents with low birth weight measured by an automated MRI-based method. Brain 128: 2588–2596.
  17. 17. Hatfield T, Wing DA, Buss C, Head K, Muftuler LT, et al. (2011) Magnetic resonance imaging demonstrates long-term changes in brain structure in children born preterm and exposed to chorioamnionitis. Am J Obstet Gynecol 205: 384.
  18. 18. Skranes J, Lohaugen GCC, Evensen KA, Indredavik MS, Haraldseth O, et al. (2012) Entorhinal cortical thinning affects perceptual and cognitive functions in adolescents born preterm with very low birth weight (VLBW). Early Hum Dev 88: 103–109.
  19. 19. Pryce CR, Bettschen D, Bahr NI, Feldon J (2001) Comparison of the effects of infant handling, isolation, and nonhandling on acoustic startle, prepulse inhibition, locomotion, and HPA activity in the adult rat. Behav Neurosci 115: 71–83.
  20. 20. Maccari S, Morley-Fletcher S (2007) Effects of prenatal restraint stress on the hypothalamus-pituitary-adrenal axis and related behavioural and neurobiological alterations. Psychoneuroendocrinology 32 Suppl 1S10–S15.
  21. 21. Brummelte S, Pawluski JL, Galea LA (2006) High post-partum levels of corticosterone given to dams influence postnatal hippocampal cell proliferation and behavior of offspring: A model of post-partum stress and possible depression. Horm Behav 50: 370–382.
  22. 22. Champagne FA, Meaney MJ (2006) Stress during gestation alters postpartum maternal care and the development of the offspring in a rodent model. Biol Psychiatry 59: 1227–1235.
  23. 23. Bhutta AT, Rovnaghi C, Simpson PM, Gossett JM, Scalzo FM, et al. (2001) Interactions of inflammatory pain and morphine in infant rats: long-term behavioral effects. Physiology & Behavior 73: 51–58.
  24. 24. Walker CD, Xu Z, Rochford J, Johnston CC (2008) Naturally occurring variations in maternal care modulate the effects of repeated neonatal pain on behavioral sensitivity to thermal pain in the adult offspring. Pain 140: 167–176.
  25. 25. de Medeiros CB, Fleming AS, Johnston CC, Walker CD (2009) Artificial rearing of rat pups reveals the beneficial effects of mother care on neonatal inflammation and adult sensitivity to pain. Pediatr Res 66: 272–277.
  26. 26. Low LA, Fitzgerald M (2012) Acute Pain and a Motivational Pathway in Adult Rats: Influence of Early Life Pain Experience. PLoS ONE 7: e34316.
  27. 27. Anand KJ, Coskun V, Thrivikraman KV, Nemeroff CB, Plotsky PM (1999) Long-term behavioral effects of repetitive pain in neonatal rat pups. Physiol Behav 66: 627–637.
  28. 28. Lupien SJ, McEwen BS, Gunnar MR, Heim C (2009) Effects of stress throughout the lifespan on the brain, behaviour and cognition. Nat Rev Neurosci 10: 434–445.
  29. 29. McEwen BS (2007) Physiology and neurobiology of stress and adaptation: central role of the brain. Physiol Rev 87: 873–904.
  30. 30. Volpe JJ (2009) Brain injury in premature infants: a complex amalgam of destructive and developmental disturbances. Lancet Neurol 8: 110–124.
  31. 31. Grunau RE, Whitfield MF, Petrie-Thomas J, Synnes AR, Cepeda IL, et al. (2009) Neonatal pain, parenting stress and interaction, in relation to cognitive and motor development at 8 and 18 months in preterm infants. Pain 143: 138–146.
  32. 32. Grunau RE, Haley DW, Whitfield MF, Weinberg J, Yu W, et al. (2007) Altered basal cortisol levels at 3, 6, 8 and 18 months in infants born at extremely low gestational age. J Pediatr 150: 151–156.
  33. 33. Hamilton LS, Narr KL, Luders E, Szeszko PR, Thompson PM, et al. (2007) Asymmetries of cortical thickness: effects of handedness, sex, and schizophrenia. Neuroreport 18: 1427–1431.
  34. 34. Toga AW, Thompson PM (2003) Mapping brain asymmetry. Nat Rev Neurosci 4: 37–48.
  35. 35. Wechsler D (2003) Wechsler Intelligence Scale for Children (WISC-IV). San Antonio, TX: Psychological Corporation.
  36. 36. Miller SP, Ferriero DM, Leonard C, Piecuch R, Glidden DV, et al. (2005) Early brain injury in premature newborns detected with magnetic resonance imaging is associated with adverse early neurodevelopmental outcome. J Pediatr 147: 609–616.
  37. 37. Gibson E, Wang L, Beg M (2009) Cortical thickness measurement using Eulerian PDEs and surface-based global topological information. Org Human Brain Mapping, 15th Ann Meeting.
  38. 38. Jones SE, Buchbinder BR, Aharon I (2000) Three-dimensional mapping of cortical thickness using Laplace’s Equation. Hum Brain Mapp 11: 12–32.
  39. 39. Lerch JP, Evans AC (2005) Cortical thickness analysis examined through power analysis and a population simulation. NeuroImage 24: 163–173.
  40. 40. Clarkson MJ, Cardoso MJ, Ridgway GR, Modat M, Leung KK, et al. (2011) A comparison of voxel and surface based cortical thickness estimation methods. NeuroImage 57: 856–865.
  41. 41. Yezzi AJ Jr, Prince JL (2003) An Eulerian PDE approach for computing tissue thickness. IEEE Trans Med Imaging 22: 1332–1339.
  42. 42. Dale AM, Fischl B, Sereno MI (1999) Cortical surface-based analysis. I. Segmentation and surface reconstruction. NeuroImage 9: 179–194.
  43. 43. Fischl B, Sereno MI, Dale AM (1999) Cortical surface-based analysis. II: Inflation, flattening, and a surface-based coordinate system. NeuroImage 9: 195–207.
  44. 44. Desikan RS, Segonne F, Fischl B, Quinn BT, Dickerson BC, et al. (2006) An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. NeuroImage 31: 968–980.
  45. 45. Richardson DK, Corcoran JD, Escobar GJ, Lee SK (2001) SNAP-II and SNAPPE-II: Simplified newborn illness severity and mortality risk scores. J Pediatr 138: 92–100.
  46. 46. Grunau RE, Holsti L, Haley DW, Oberlander T, Weinberg J, et al. (2005) Neonatal procedural pain exposure predicts lower cortisol and behavioral reactivity in preterm infants in the NICU. Pain 113: 293–300.
  47. 47. Hunter M, Takane Y (2002) Constrained principal component analysis: Various applications. J Educ Behav Stat 27: 41–81.
  48. 48. Takane Y, Shibayama T (1991) Principal component analysis with external information on both subjects and variables. Psychometrika 56: 97–120.
  49. 49. Takane Y, Hunter M (2001) Constrained principal component analysis: A comprehensive theory. Applicable Algebra in Engineering, Communication, and Computing 12: 391–419.
  50. 50. Efron B, Tibshirani RJ (1994) An Introduction to the Bootstrap. Boca Raton, FL: Chapman & Hall/CRC.
  51. 51. Shao J, Tu D (1995) The Jacknife and Bootstrap. New York, NY: Springer.
  52. 52. Benjamini Y, Hochberg Y (1995) Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J R Stat Soc Series B Stat Methodol 57: 289–300.
  53. 53. Bramen JE, Hranilovich JA, Dahl RE, Chen J, Rosso C, et al. (2012) Sex Matters during Adolescence: Testosterone-Related Cortical Thickness Maturation Differs between Boys and Girls. PLoS ONE 7: e33850.
  54. 54. Dean JM, McClendon E, Hansen K, Azimi-Zonooz A, Chen K, et al. (2013) Prenatal Cerebral Ischemia Disrupts MRI-Defined Cortical Microstructure Through Disturbances in Neuronal Arborization. Sci Transl Med 5: 168ra7.
  55. 55. Doesburg SM, Chau CM, Cheung TPL, Moiseev A, Ribary U, et al.. (2013) Neonatal pain-related stress, functional cortical activity and visual-perceptual abilities in school-age children born at extremely low gestational age. Pain [Epub ahead of print] doi: 10.1016/j.pain.2013.04.009.
  56. 56. Shaw P, Kabani NJ, Lerch JP, Eckstrand K, Lenroot R, et al. (2008) Neurodevelopmental trajectories of the human cerebral cortex. J Neurosci 28: 3586–3594.
  57. 57. Duhrsen L, Simons SHP, Dzietko M, Genz K, Bendix I, et al. (2013) Effects of Repetitive Exposure to Pain and Morphine Treatment on the Neonatal Rat Brain. Neonatology 103: 35–43.
  58. 58. Juul SE, Beyer RP, Bammler TK, Farin FM, Gleason CA (2011) Effects of neonatal stress and morphine on murine hippocampal gene expression. Pediatr Res 69: 285–292.
  59. 59. Anand KJ, Garg S, Rovnaghi CR, Narsinghani U, Bhutta AT, et al. (2007) Ketamine reduces the cell death following inflammatory pain in newborn rat brain. Pediatr Res 62: 283–290.
  60. 60. Kesler SR, Ment LR, Vohr B, Pajot SK, Schneider KC, et al. (2004) Volumetric analysis of regional cerebral development in preterm children. Pediatr Neurol 31: 318–325.
  61. 61. Kesler SR, Reiss AL, Vohr B, Watson C, Schneider KC, et al. (2008) Brain Volume Reductions within Multiple Cognitive Systems in Male Preterm Children at Age Twelve. J Pedr 152: 513–520.
  62. 62. Peterson BS, Vohr B, Staib LH, Cannistraci CJ, Dolberg A, et al. (2000) Regional brain volume abnormalities and long-term cognitive outcome in preterm infants. JAMA 284: 1939–1947.
  63. 63. Slater R, Fabrizi L, Worley A, Meek J, Boyd S, et al. (2010) Premature infants display increased noxious-evoked neuronal activity in the brain compared to healthy age-matched term-born infants. NeuroImage 52: 583–589.
  64. 64. Anand KJ (2000) Effects of perinatal pain and stress. Prog Brain Res 122: 117–129.
  65. 65. Laprairie JL, Johns ME, Murphy AZ (2008) Preemptive morphine analgesia attenuates the long-term consequences of neonatal inflammation in male and female rats. Pediatr Res 64: 625–630.
  66. 66. Apkarian AV, Bushnell MC, Treede RD, Zubieta JK (2005) Human brain mechanisms of pain perception and regulation in health and disease. Eur J Pain 9: 463–484.
  67. 67. Casey KL, Svensson P, Morrow TJ, Raz J, Jone C, et al. (2000) Selective Opiate Modulation of Nociceptive Processing in the Human Brain. J Neurophysiol 84: 525–533.
  68. 68. Treede RD, Apkarian AV (2008) Nociceptive Processing in the Cerebral Cortex. In: Allan I, Akimichi K, Gordon MS, Gerald W, Thomas DA et al.., editors. The Senses: A Comprehensive Reference. New York: Academic Press. pp. 669–697.
  69. 69. Schlaepfer TE, Strain EC, Greenberg BD, Preston KL, Lancaster E, et al. (1998) Site of Opioid Action in the Human Brain: Mu and Kappa Agonists’ Subjective and Cerebral Blood Flow Effects. Am J Psychiatry 155: 470–473.
  70. 70. Hodge CJ Jr, Apkarian AV (1990) The spinothalamic tract. Crit Rev Neurobiol 5: 363–397.
  71. 71. Treede RD, Kenshalo DR, Gracely RH, Jones AKP (1999) The cortical representation of pain. Pain 79: 105–111.
  72. 72. Derbyshire SWG, Osborn J (2007) Modeling Pain Circuits: How Imaging May Modify Perception. Neuroimaging Clin N A 17: 485–493.
  73. 73. Bellu R, de Waal KA, Zanini R (2008) Opioids for neonates receiving mechanical ventilation. Cochrane Database Syst Rev CD004212.
  74. 74. Anand KJS, Hall RW, Desai N, Shephard B, Bergqvist LL, et al. (2004) Effects of morphine analgesia in ventilated preterm neonates: primary outcomes from the NEOPAIN randomised trial. The Lancet 363: 1673–1682.
  75. 75. Simons SH, van DM, van Lingen RA, Roofthooft D, Duivenvoorden HJ, et al. (2003) Routine morphine infusion in preterm newborns who received ventilatory support: a randomized controlled trial. JAMA 290: 2419–2427.
  76. 76. de Graaf J, van Lingen RA, Simons SHP, Anand KJS, Duivenvoorden HJ, et al. (2011) Long-term effects of routine morphine infusion in mechanically ventilated neonates on childrens functioning: Five-year follow-up of a randomized controlled trial. Pain 152: 1391–1397.
  77. 77. Johnston CC, Collinge JM, Henderson SJ, Anand KJ (1997) A cross-sectional survey of pain and pharmacological analgesia in Canadian neonatal intensive care units. Clin J Pain 13: 308–312.
  78. 78. Roze J, Denizot S, Carbajal R, Ancel P, Kaminski M, et al. (2008) Prolonged sedation and/or analgesia and 5-year neurodevelopment outcome in very preterm infants: Results from the epipage cohort. Arch Pediatr Adolesc Med 162: 728–733.
  79. 79. Ferguson SA, Ward WL, Paule MG, Hall RW, Anand KJS (2012) A pilot study of preemptive morphine analgesia in preterm neonates: Effects on head circumference, social behavior, and response latencies in early childhood. Neurotoxicol Teratol 34: 47–55.
  80. 80. de Graaf J, van Lingen RA, Valkenburg AJ, Weisglas-Kuperus N, Groot Jebbink L, et al. (2013) Does neonatal morphine use affect neuropsychological outcomes at 8 to 9years of age? Pain 154: 449–458.
  81. 81. Zwicker JG, Grunau RE, Adams E, Chau V, Brant R, et al. (2013) Score for Neonatal Acute Physiology II and Neonatal Pain Predict Corticospinal Tract Development in Premature Newborns. Pediatr Neurol 48: 123–129.
  82. 82. Grunau RE, Holsti L, Peters JWB (2006) Long-term consequences of pain in human neonates. Semin Fetal Neonatal Med 11: 268–275.
  83. 83. Hohmeister J, Kroll A, Wollgarten-Hadamek I, Zohsel K, Demirakça S, et al. (2010) Cerebral processing of pain in school-aged children with neonatal nociceptive input: An exploratory fMRI study. Pain 150: 257–267.
  84. 84. Tukey JW (1977) Exploratory Data Analysis. Don Mills, Ontario: Addison-Wesley.