key: cord-0016241-zbb8wemp authors: Rantakari, Krista; Rinta-Koski, Olli-Pekka; Metsäranta, Marjo; Hollmén, Jaakko; Särkkä, Simo; Rahkonen, Petri; Lano, Aulikki; Lauronen, Leena; Nevalainen, Päivi; Leskinen, Markus J.; Andersson, Sture title: Early oxygen levels contribute to brain injury in extremely preterm infants date: 2021-03-22 journal: Pediatr Res DOI: 10.1038/s41390-021-01460-3 sha: a83663a1d618ae1d90b541b6f80fee44cba05caf doc_id: 16241 cord_uid: zbb8wemp BACKGROUND: Extremely low gestational age newborns (ELGANs) are at risk of neurodevelopmental impairments that may originate in early NICU care. We hypothesized that early oxygen saturations (SpO(2)), arterial pO(2) levels, and supplemental oxygen (FiO(2)) would associate with later neuroanatomic changes. METHODS: SpO(2), arterial blood gases, and FiO(2) from 73 ELGANs (GA 26.4 ± 1.2; BW 867 ± 179 g) during the first 3 postnatal days were correlated with later white matter injury (WM, MRI, n = 69), secondary cortical somatosensory processing in magnetoencephalography (MEG-SII, n = 39), Hempel neurological examination (n = 66), and developmental quotients of Griffiths Mental Developmental Scales (GMDS, n = 58). RESULTS: The ELGANs with later WM abnormalities exhibited lower SpO(2) and pO(2) levels, and higher FiO(2) need during the first 3 days than those with normal WM. They also had higher pCO(2) values. The infants with abnormal MEG-SII showed opposite findings, i.e., displayed higher SpO(2) and pO(2) levels and lower FiO(2) need, than those with better outcomes. Severe WM changes and abnormal MEG-SII were correlated with adverse neurodevelopment. CONCLUSIONS: Low oxygen levels and high FiO(2) need during the NICU care associate with WM abnormalities, whereas higher oxygen levels correlate with abnormal MEG-SII. The results may indicate certain brain structures being more vulnerable to hypoxia and others to hyperoxia, thus emphasizing the role of strict saturation targets. IMPACT: This study indicates that both abnormally low and high oxygen levels during early NICU care are harmful for later neurodevelopmental outcomes in preterm neonates. Specific brain structures seem to be vulnerable to low and others to high oxygen levels. The findings may have clinical implications as oxygen is one of the most common therapies given in NICUs. The results emphasize the role of strict saturation targets during the early postnatal period in preterm infants. Extremely low gestational age newborns (ELGANs), born before 28 weeks of gestational age (GA), are at risk of brain injury and later neurodevelopmental complications. [1] [2] [3] [4] [5] [6] [7] [8] Many of these complications are thought to originate during their early care, when they are exposed to various factors, such as oxygen, which is one of the most common therapies given in neonatal intensive care units (NICUs). 1, 2, 7, [9] [10] [11] [12] [13] [14] Inappropriately low oxygen levels are, e.g., associated with increased mortality and impaired neurodevelopment, 9, 15, 16 and higher levels with retinopathy and lung injury. 16 However, despite intensive research, the long-term effects of early oxygen levels in preterm infants are incompletely characterized. 7, [10] [11] [12] Medical data from electronic health record systems and monitoring devices are increasingly available and methods for analyzing large data sets have improved, making analysis of NICU information system data feasible. As a part of our multimethodological study of ELGANs, we correlated oxygen-related parameters from the first 3 postnatal days of life with neuroimaging (magnetic resonance imaging (MRI)) and magnetoencephalography (MEG) at term equivalent age (TEA) and neurodevelopmental outcomes at 2 years of corrected age. The main outcomes were white matter (WM) injury in brain MRI and changes in secondary cortical somatosensory processing measured by MEG-SII, as these abnormalities may originate during the neonatal phase. 8, [17] [18] [19] [20] [21] [22] [23] Neurodevelopmental outcome was assessed by Hempel neurological examination and the Griffiths Mental Developmental Scales (GMDS). We hypothesized that, in ELGANs aberrant oxygen saturations (SpO 2 ), lower arterial partial pressures of oxygen (pO 2 ), and higher oxygen demand (as measured by fraction of inspired oxygen (FiO 2 )) during the first 3 postnatal days would correlate with pathologic changes in neuroanatomical structures and unfavorable neurodevelopmental outcomes. The original patient group consisted of 82 ELGANs, who were born before 28 gestational weeks and treated in level IV NICU at Children's Hospital, Helsinki, Finland. They were recruited for a multimethodological study between 5/2006 and 9/2008 and the parents provided signed informed consents. All clinical decisions were made by existing treatment protocols. The SpO 2 target range was 90-95%, the blood transfusion range was hematocrit (HCT) <40% when on ventilator care and/or needing supplemental oxygen, and HCT <30% for more stable infants. No delayed cord clamping was performed at birth. Six infants died during the NICU period and three infants did not participate in the follow-up examinations. Consequently, the final study population included 73 ELGANs (Fig. 1) . The clinical characteristics of the infants are presented in Table 1A NICU biosignal data and oxygen-related parameters Data gathered during the NICU phase, including patient monitor data, data from ventilators and other medical devices, observation variables, laboratory results, and diagnoses were stored in Centricity Critical Care information system (GE Healthcare, Chicago, IL). Monitor and ventilator data including SpO 2 , FiO 2 , and heart rate (HR) were recorded automatically and stored as 2min averages of median values for 10-s intervals. The observation variables include both background information (e.g., GA, mode of delivery, sex) and clinical measurements (e.g., SpO 2 , FiO 2 ). Manual readings were entered into the database by staff and represent the value at the time of reading. The partial pressures of oxygen (pO 2 ) and carbon dioxide (pCO 2 ) in arterial blood were measured using routine blood gas analyzers. Brain MRI and determination of WM abnormalities Brain MRI (1.5 T) including T2-weighted axial and T1-weighted 3D sagittal images was performed at TEA. Two experienced neuroradiologists classified the MRI images according to Woodward et al. with modifications. 17 WM was classified based on five variables: WM signal abnormality, periventricular WM volume loss, cystic abnormalities, ventricular dilatation, and thinning of the corpus callosum. The thinning of corpus callosum had two grades (grade 1 was normal and grade 2 had thinning of the corpus callosum). All the other variables had three grades, grade 1 being normal and grade 3 having most severe abnormalities. Based on the total score, the WM MRI findings were classified as abnormal if the total score was ≥7. In a subgroup analysis, comparisons were performed between patients with normal WM and those with severe WM changes having total scores of ≥10. MEG and sensory evoked magnetic fields MEG analysis was performed at TEA as previously described. 18 MEG recording was performed using a whole-head adult-sized helmet-shaped sensor array consisting of 306 independent channels: 204 gradiometers and 102 magnetometers (Vectorview, Elekta Neuromag Oy, Helsinki, Finland). Electroencephalography and electro-oculography were recorded for sleep stage monitoring. The sensory evoked magnetic fields were elicited by tactile stimulus to the tip of the index finger by a thin elastic membrane expanded by an air pressure pulse delivered through a plastic tube (Somatosensory Stimulus Generator, 4-D NeuroImaging Inc., San Diego, CA) and the interstimulus interval was 2 s. The MEG data were analyzed as previously described. 18, 19 The SII response, peaking at about 200 ms after tactile stimulation, was defined as abnormal when it was absent after both contralateral and ipsilateral stimulation in at least one hemisphere (right or left). 18 Neurodevelopmental assessments An experienced child neurologist performed the clinical neurodevelopmental assessments at 2 years of corrected age according to a structured Hempel neurological examination 24 and GMDS. 25, 26 The Hempel neurological examination consisted of five functional domains: posture and muscle tone, gross motor function, fine motor function, reflexes, and visual system. For the present study, infants with any dysfunctional domain(s) were classified as abnormal and the patients with normal findings as normal. In a subgroup analysis, comparisons were performed between patients with normal Hempel assessments and those with major neurologic impairment (cerebral palsy (CP)). The general developmental quotient (DQ) of GMDS was based on five subscales: locomotor, personal-social, hearing-language, eye-hand coordination, and performance. DQ and subscale quotients (SQs) were calculated on the basis of the raw scores and the corrected age of the child. DQ and SQ scores <−1 SD were defined as abnormal and the cut-off points for impairment were for general DQ 88.7, locomotor 84.3, personal-social 84.8, hearing-language 84.6, eye-hand coordination 84.3, and for performance 84.4. In addition to separate SQ analysis, we also compared patients having subscale impairments in either eye-hand coordination or performance scores or both (abnormal group) with patients having normal scores in both these two subscales (i.e., eye-hand coordination and performance; normal group). Data are presented as n (%) or mean (SD) or median (95% CI). Statistically significant p-values are in bold. Early oxygen levels contribute to brain injury in extremely preterm infants K Rantakari et al. The methods of big data analysis were used, not due to the number of the patients (final n = 73) but rather due to large and complex NICU database, e.g., altogether over million SpO 2 measurements. Time series data (measurements with associated time stamps) for oxygen-related parameters were extracted from Centricity for Critical Care information system and analyzed at the Aalto University Department of Computer Science, Espoo, Finland using PostgreSQL database engine. SpO 2 and FiO 2 data were extracted from data logged automatically via equipment interfaces. Data for pO 2 and pCO 2 analysis came from laboratory information system via Centricity for Critical Care information system interface. Supplementary oxygen time series was created by combining automatically logged FiO 2 from ventilator interface when available (mechanical ventilation, synchronized nasal ventilation) and manually entered FiO 2 for patients on nasal continuous positive airway pressure or high flow nasal cannulas, where no equipment interface was available. The time series data were preprocessed by removing out-of-range values caused by, e.g., missing or misplaced sensors and monitoring equipment drifting out of calibration. The data were correlated with WM injury, MEG-SII abnormalities (both at TEA), and Hempel neurological examination and GMDS quotients (at 2 years of corrected age). Student's T test was used for statistical comparisons. A p value <0.05 was considered significant. The study design is shown in Fig. 1 and the clinical characteristics of the ELGANs (final n = 73) in Table 1A , B. WM (at TEA) was classified as abnormal in 38%, MEG-SII (at TEA) in 36%, Hempel neurological examination (at 2 years of corrected age) in 61%, and GMDS quotients in 41% of the patients. These relatively large percentages of infants in the abnormal groups are likely to be explained by the study design with classifying infants even with mild alterations in the main abnormal groups. In the subgroup analysis, there were 6 infants with severe WM changes (8.7%) and 6 with major neurologic impairment in Hempel assessments (CP, 9.1%). In comparison with infants with normal WM, those with WM abnormalities were more often males, had tendency toward lower Apgar scores, had more often suffered from respiratory distress syndrome (RDS), and had more intraventricular hemorrhage (IVH) findings in ultrasound examinations (Table 1A) . The ELGANs with severe WM changes additionally suffered more from bronchopulmonary dysplasia (BPD) (data not shown (DNS)). The patients with abnormal Hempel neurological examination were more often males (Table 1B) , and the ones with major neurologic impairment (CP) also had more IVH (DNS) than the infants with normal Hempel assessment. The infants with low GMDS scores, in turn, had more retinopathy of prematurity than those with normal scores (Table 1B) . The infants with abnormal WM, MEG-SII, Hempel, or low GMDS did not significantly differ from those without these abnormalities in terms of GA, birth weight, being small for gestational age (SGA), being twins, or other parameters shown in Table 1 . In addition, none of these groups statistically significantly differed from the other groups in regard to age, parity, and the smoking status of the mother; being exposed to preterm premature rupture of membranes, chorioamnionitis, pre-eclampsia, gestational diabetes, cesarean section; being extubated by day 3; or treatment for patent ductus arteriosus (DNS). The infants with abnormal WM and/or MEG-SII had tendencies toward having received more postnatal dexamethasone than the ones with normal imaging, but findings were not statistically significant (abnormal WM 19% vs normal WM 9%, p = 0.07 and MEG-SII 14 vs 2%, p = 0.09). Severe WM changes and abnormal MEG-SII were associated with adverse neurodevelopment The severe WM changes at TEA were associated with abnormalities in GMDS and Hempel assessments at two years of corrected age (all p < 0.05, DNS). The abnormal WM group that included also minor WM changes was not significantly associated with the 2year neurodevelopmental assessments (GMDS or Hempel, DNS). The abnormal MEG-SII findings (at TEA) were significantly associated with worse outcomes in 2-year GMDS, which is consistent with our previous findings. 18 Here, with a larger study group than in our previous study, 18 Early oxygen levels contribute to brain injury in extremely preterm infants K Rantakari et al. subscales were statistically significant: the infants with abnormal MEG-SII had significantly lower GMDS SQ scores in locomotor, personal and social, hearing and language, eye-hand coordination, and performance (all p < 0.05, DNS). The patients with abnormal MEG-SII findings tended also to have poorer Hempel outcomes, although the finding was not statistically significant (p = 0.07). Oxygen saturations and partial pressures of oxygen during the first 3 days of life in ELGANs with later WM, MEG-SII, and neurodevelopmental abnormalities During the first 3 days of life, the infants who at TEA displayed WM abnormalities showed lower average SpO 2 (92.8% ± 0.3 vs 93.8% ± 0.3; p < 0.05. DNS), had more SpO 2 <90% (p < 0.05) and <85% (p < 0.05), and less SpO 2 >95% (p < 0.05) than the patients with normal WM (Fig. 2a) The infants with abnormal MEG-SII, in turn, had higher average SpO 2 (94.5% ± 0.4 vs 92.8% ± 0.4 p < 0.05, DNS) and had more SpO 2 >95% (p < 0.05) and >98% (p < 0.05), as well as fewer SpO 2 <90% (p < 0.05), than the patients with normal MEG-SII findings (Fig. 2b) . They also stayed less within current oxygen saturation target ranges (90-95%, p < 0.05). Furthermore, their cumulative times (percentage of time during the first 3 days) when having higher SpO 2 (>95% and >98%) were greater. The proportional times (percentage of time) in patients with abnormal and normal MEG-SII when having SpO 2 <85% were 1.7 vs 3.1% (NS), SpO 2 <90% 9.2 vs 17.6 (p < 0.05), SpO 2 90-95% 40.1 vs 51.6% (p < 0.05), SpO 2 >95% 47.6 vs 26.5% (p < 0.05), and SpO 2 >98% 17.1 vs 6.8% (p < 0.05), respectively. The ELGANs with low GMDS SQ scores in eye-hand coordination or performance or both had significantly fewer low SpO 2 (i.e., had higher SpO 2 ) than their controls with better scores (DNS). The SpO 2 levels in other abnormal GMDS subscales did not significantly differ from the ones with normal findings (DNS). The infants with major neurologic impairment in Hempel assessment (CP) had more SpO 2 measurements <85% and higher cumulative times when having SpO 2 <85% than the ones with Early oxygen levels contribute to brain injury in extremely preterm infants K Rantakari et al. normal Hempel examinations (all p < 0.05, DNS). The oxygen saturations in those abnormal Hempel groups that included also the minor impairments did not significantly differ from the group with normal assessments (DNS). Consistently with the SpO 2 findings, the infants with WM abnormalities had lower arterial pO 2 levels and patients with abnormal MEG-SII had higher arterial pO 2 levels than the infants with normal WM and MEG-SII, respectively, during the first 3 days of life (Fig. 3) . Also, the ELGANs with low SQ scores in locomotor, eye-hand coordination, or performance had significantly more high and fewer low pO 2 levels (p < 0.05, DNS) than the ELGANs with normal scores (DNS). The arterial pO 2 levels of the other abnormal GMDS SQs or Hempel examinations did not significantly differ from the ones with normal examinations (DNS). Patients with abnormal WM findings and major neurologic impairment in Hempel assessment (CP) had higher and patients with abnormal MEG-SII or low GMDS had lower need for supplemental oxygen (FiO 2 ) The patients with WM abnormalities had spent more time (percentage of time) with supplemental oxygen (FiO 2 >21%) and had required higher average FiO 2 than the ones with normal WM (Fig. 4) . The findings with severe WM abnormalities were similar (p < 0.01, DNS). Also the proportional times with FiO 2 >30% (p < 0.01) and >70% (p < 0.05) were higher in infants with abnormal WM (Fig. 4c) and with severe WM outcomes (p < 0.01, DNS) than in those with normal WM. The ELGANs with major neurologic impairment in Hempel assessment (CP) also had a higher need for FiO 2 , i.e., they had higher proportional times with FiO 2 >70% (p < 0.01) and >90% (p < 0.05) than the ones with normal Hempel assessments (DNS). In turn, the infants with abnormal MEG-SII had spent less time with supplemental oxygen and received less FiO 2 than the patients with normal MEG-SII (Fig. 4a, b) . Also the group with lower GMDS scores spent significantly less time with FiO 2 >50% than the ones with better scores (p < 0.01; DNS). To investigate the role of excessive oxygen administration, we analyzed the use of supplemental oxygen when the saturations were >95% or >98%. During the first 3 days of life, the cumulative times when patients received supplemental oxygen while having SpO 2 >95% or >98% were not statistically significantly different between any of the study groups (DNS). Bradycardia (HR <100 bpm) and low oxygen saturations (SpO 2 ) As bradycardia, in addition to oxygen-related factors, or in combination with hypoxia (apneas), is potentially associated with later neurodevelopmental impairments, we next studied the role of HR, when being <100 bpm. During the first 3 days of life, the ELGANs with WM abnormalities tended to have more bradycardia than those with normal WM. In patients with abnormal and normal WM, the proportions of HR values <100 bpm relative to all HR values (percentage of total measurements) were 0.10 and 0.05% (p = 0.06), and the cumulative times (percentage of time during the first 3 days) when having HR <100 min were 0.18 and 0.15% (NS), respectively. The infants with abnormal WM tended also to suffer more from the combination of bradycardia and hypoxia as measured by the cumulative time (percentage of time during the first 3 days) when having HR <100 bpm and SpO 2 <85% (p = 0.05) or HR <100 bpm and SpO 2 <90% (n = 0.08) than the patients with normal WM (Fig. 5a) . The differences between normal and abnormal MEG-SII, Hempel, and GMDS (DQ) groups with HR data were not significant (DNS). Early oxygen levels contribute to brain injury in extremely preterm infants K Rantakari et al. In arterial blood gas analysis, the infants with WM abnormalities had higher average levels of pCO 2 (Fig. 5b) than those with normal WM. Consistently, when analyzing different pCO 2 levels, the patients with abnormal WM had more high level and fewer low level pCO 2 values than those with normal WM (Fig. 5b) . Regarding pCO 2 levels, the patients with abnormal MEG-SII, low GMDS scores, or abnormal Hempel neurological examination did not statistically differ from those with normal findings (DNS). Due to advancements in neonatal intensive care and increased survival of ELGANs, the prediction and prevention of later complications have become crucially important. [4] [5] [6] [7] [8] [21] [22] [23] In the present study, we demonstrate that the ELGANs with later WM abnormalities exhibit lower actually achieved SpO 2 during their first 3 days of life than the infants with normal WM. In turn, the results are different with infants having abnormal MEG-SII responses, who have higher early SpO 2 than ELGANs with normal MEG-SII. We also show that both these abnormalities, i.e., severe WM changes and abnormal MEG-SII at TEA, correlate with adverse neurodevelopment at 2 years of corrected age. The role of WM injury in unfavorable development is in accordance with previous publications. 8, 17, [20] [21] [22] [23] The SII response, in turn, is considered as an indicator of higher intracortical processing, and the absence of this response may reflect overall reduced corticocortical connectivity. 19 The findings with WM and oxygen-linked factors are also in agreement with previous studies showing relations between low oxygen and WM injury, 1, 20, 21, 27 although hyperoxia has also been suggested to be harmful to WM. 28 In terms of oxygen parameters and MEG-SII responses in infants, there are no previous data and the present findings are novel. Consistent with the SpO 2 findings, the achieved oxygen levels in arterial samples (pO 2 ) in ELGANs with WM abnormalities were lower and with MEG-SII and GMDS abnormalities higher than in their normal controls. Thus, despite technical, limit setting, and averaging challenges with pulse oximetry, 14, 29 and discontinuous data sampling with arterial pO 2 levels, obtaining similar results with both methods strengthens the results. The other limitations Early oxygen levels contribute to brain injury in extremely preterm infants K Rantakari et al. of the study include the number of patients, i.e., although the amount of data from the database was big, the number of patients (final study population of 73 patients) may not have been powerful enough for discovering all clinically significant factors. Furthermore, retrospective investigations always have limitations. In terms of updates in later devices, we consider the data obtained from the monitoring devices of our study compatible, as, e.g., the limited setting and averaging challenges still exist with the updated pulse oximetry. The possible role of compliance with preestablished SpO 2 ranges by the stuff was not determined, but one could assume that the potential effects of compliance would be equally targeted to all different groups and thus compliance would not explain the observed differences between the patient groups. The mechanisms underlying the opposite findings in early oxygen levels in patients with later WM injury and with abnormal MEG-SII responses are unclear. Only two of the patients had both abnormal WM and abnormal MEG-SII responses, whereas the others had either abnormal WM or MEG-SII. This suggests that the noxae causing these two conditions are likely to be different. Indeed, certain cell types and tissue structures have been shown to be more sensitive to hypoxia and others to hyperoxia. Hypoxia can, e.g., lead to proinflammatory and oxidative stress, 1 increased production of hypoxia-inducible factors, 30 and WM injury. 21 Hyperoxia, in turn, has an impact on oxidative stress, changes in blood flow, disruption on neural plasticity and myelination, and may contribute to the encephalopathy of prematurity. 1, 28, 31 Thus, the present results may indicate certain brain structures being more vulnerable to hypoxia and others to higher oxygen levels. The administration of supplemental oxygen (FiO 2 ) seemed, at first glance, to have been clinically appropriate, as the patients with lower SpO 2 (abnormal WM) were supplied with higher FiO 2 and the patients with higher SpO 2 (abnormal MEG-SII) with less supplemental oxygen (Fig. 4) . Moreover, the cumulative times of inappropriate oxygen administration did not differ between the groups. However, the administration of appropriate FiO 2 in the NICU may not have been dynamic enough as the actually achieved oxygen saturations were relatively often out-of-thetarget range (Fig. 2) . Having said that, certain high SpO 2 levels were not iatrogenically caused as some infants had high SpO 2 at room air, i.e., without supplementary oxygen. Nevertheless, those occasions may still have been harmful, as even at room air the oxygen exposure is higher than that encountered in corresponding fetal period in utero. 2, 10, 32 Other than oxygen-associated factors are undoubtedly likely to be involved in causing neurological impairments in ELGANs. In the present study, the infants with WM abnormalities had lower Apgar scores, more RDS, BPD, and IVH, and thus seem to have been in more serious general condition during the NICU period, than those with normal WM. Moreover, they were more often male, in whom the antioxidant defense maturation is shown to be delayed specifically when born preterm. 32 All these factors may have affected the later outcomes of these ELGANs. The SGA or twin infants, in turn, were not significantly overrepresented in the abnormal groups, but the study is limited by the number of patients and may not have been powerful enough to discover all clinically significant factors. There are several additional factors that may play roles in inducing neurodevelopmental impairments in ELGANs. Indeed, in our study, the partial pressures of carbon dioxide in patients with abnormal WM were higher than in those with normal WM. This is in accordance with previous studies showing the involvement of hypercarbia and pCO 2 fluctuations in adverse neurological outcomes. [33] [34] [35] [36] [37] Furthermore, vascularization and circulation are candidates for playing important roles. For example, vascularization of specific regions and the different regulation of vascular tone by oxygen and CO 2 have been demonstrated with different regions of retina. 38 Moreover, although in this study the findings with bradycardia (HR <100 bpm ± low SpO 2 ) did not reach statistical significance, there were trends, and apneas (low SpO 2 and low HR) as well as intermittent hypoxemia or hyperoxemia, systemic circulatory factors (e.g., HR, blood pressure), and hemoglobin levels may play important roles and need to be evaluated in the future. Additionally, we wish to correlate the early NICU data with the neurodevelopmental outcome at later ages, such as pre-school and school ages. Naturally, prospective investigations are needed to confirm the effects of early oxygen parameters and other related factors in later neurodevelopment in ELGANs. In conclusion, lower SpO 2 and pO 2 levels and higher FiO 2 need during the first 3 days of life in ELGANs associate with later WM abnormalities and adverse neurodevelopment. The infants with abnormal MEG-SII, which is associated with later low GMDS, show the opposite findings, i.e., they have higher early SpO 2 and pO 2 levels, than the infants with more favorable MEG-SII and GMDS outcomes. The results indicate that of brain structures WM are more vulnerable to hypoxia, whereas others, such as cortical areas, are more sensitive to hyperoxia. Although the findings have to be confirmed in further studies, they emphasize the importance of maintaining strict saturation targets during the early postnatal period. Cell therapy for the preterm infant: promise and practicalities Oxygen metabolism and oxygenation of the newborn Early oxygen levels contribute to brain injury in extremely preterm infants Early neurodevelopmental outcomes of extremely preterm infants Neurodevelopmental outcomes in early childhood Neurologic consequences of preterm birth EPICure Study Group. Neurologic and developmental disability at six years of age after extremely preterm birth The evolution of modern respiratory care for preterm infants Neuroimaging and neurodevelopmental outcome in extremely preterm infants Patterns of oxygenation, mortality, and growth status in the surfactant positive pressure and oxygen trial cohort Lower versus higher oxygen concentrations titrated to target oxygen saturations during resuscitation of preterm infants at birth Effects of targeting lower versus higher arterial oxygen saturations on death or disability in preterm infants Oxygen targeting in infants born extremely preterm who are small for gestational age: a need for heightened vigilance Impact of the neonatal resuscitation program-recommended low oxygen strategy on outcomes of infants born preterm Committee on Fetus and Newborn. Oxygen targeting in extremely low birth weight infants Cerebral oxygenation is associated with neurodevelopmental outcome of preterm children at age 2 to 3 years Optimal oxygenation of extremely low birth weight infants: a meta-analysis and systematic review of the oxygen saturation target studies Neonatal MRI to predict neurodevelopmental outcomes in preterm infants Cortical somatosensory processing measured by magnetoencephalography predicts neurodevelopment in extremely low-gestationalage infants Evaluation of somatosensory cortical processing in extremely preterm infants at term with MEG and EEG Impaired oligodendrocyte maturation in preterm infants: potential therapeutic targets White matter injury in the preterm infant: pathology and mechanisms MRI findings at term-corrected age and neurodevelopmental outcomes in a large cohort of very preterm infants Neurodevelopmental outcomes in preterm infants with white matter injury using a new MRI classification The Neurological Examination for Toddler Age Griffiths Entwicklungsskalen zur Beurteilung der Entwicklung in den ersten beiden Lebensjahren Atypical sensory processing is common in extremely low gestational age children Mild intermittent hypoxemia in neonatal mice causes permanent neurofunctional deficit and white matter hypomyelination Bendix, I. & Felderhoff-Mueser, U. Hyperoxia and the immature brain Oxygen targeting in preterm infants: a physiological interpretation Gene expression profiles in preterm infants on continuous long-term oxygen therapy suggest reduced oxidative stress-dependent signaling during hypoxia Oxidative stress as a primary risk factor for brain damage in preterm newborns Antenatal steroids and antioxidant enzyme activity in preterm infants: influence of gender and timing Incidence of hypocapnia, hypercapnia, and acidosis and the associated risk of adverse events in preterm neonates Carbon dioxide fluctuations are associated with changes in cerebral oxygenation and electrical activity in infants born preterm Newborn blood gas derangements of children born extremely preterm and neurocognitive dysfunctions at age 10 years PaCO2 in surfactant, positive pressure, and oxygenation randomised trial (SUPPORT) Carbon dioxide as a drug in neonatology Metabolic adaptation and neuroprotection differ in the retina and choroid in a piglet model of acute postnatal hypoxia We are deeply indebted to Marita Suni, RN, for caring for the very prematurely born infants during the MEG measurements. We also express our gratitude to all our subjects and their families for preparing the way for this study. Finally, we thank the personnel of the maternity ward of the Department of Obstetrics and Gynecology and the Neonatal Intensive Care Unit at the Helsinki University Hospital and the neonatal ward at Jorvi Hospital for their cooperation. This study was supported by Päivikki and Sakari Sohlberg Foundation (to K.R.), Foundation for Pediatric Research in Open access funding provided by University of Helsinki including Helsinki University Central Hospital. Competing interests: The authors declare no competing interests. Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.