Assessment of cortical and sub-cortical function in neonates by electrophysiological monitoring

The aim of this thesis was the assessment of cortical and sub-cortical function in neonates by electrophysiological monitoring, i.e. to evaluate the function of the neonatal cortex and brainstem through quantitative analysis of signals readily available in the NICU. These signals include amplitude-integrated EEG (aEEG), EEG, ECG and ABP. Important factors influencing neonatal cerebral function are the immaturity of preterm infants, and pathological conditions such as stroke, asphyxia and hemorrhage. Electrophysiological monitoring may be a valuable tool to monitor the development of cortical and sub-cortical function during the intensive care period. However, currently our knowledge of functional maturation and related signal changes is mostly based on visual pattern recognition. Since there are no exact reference values, maturation is difficult to assess and interpretations may vary between investigators. Digital monitoring recordings provide an opportunity towards quantitative signal analysis. This approach may lead to a more objective and uniform description of cortical and sub-cortical function. In the introductory Chapter 1, a brief overview of the mechanisms of perinatal encephalopathy is given, clinical diagnostic protocols are discussed and potential treatment strategies are highlighted. Furthermore, the origin and characteristics of the neonatal monitoring signals analyzed in this thesis are explained. The First part of this thesis focuses on traditional analysis of cortical function, quantifying maturational changes in multichannel aEEG and EEG recordings and quantifying the effects of anti-epileptic drugs on the EEG, to provide an objective and reproducible method for assessing brain maturation and long-term prognosis. In Chapter 2, multi-channel aEEG recordings were quantitatively analyzed and assessed for regional differences in maturation. We investigated 40 preterm infants (PMA 27–37 weeks) with normal follow-up at 24 months of age, at a median postnatal age of 8 days, using 4-h EEG recordings according to the international 10–20 system reduced montage. Nine channels were selected and converted to aEEG registrations. For each aEEG registration, lower margin amplitude (LMA), upper margin amplitude (UMA) and bandwidth (UMA–LMA) were calculated. In all channels, LMA showed a strong positive correlation with PMA. Below 32 weeks of PMA, LMA was =5 µV. Linear regression analysis showed a maximum LMA difference between channels of approximately 2 and 1 µV at 27 and 37 weeks of PMA, respectively. The lowest are LMA values in the occipital channel and the highest values are in centro-occipital channels. In the frontal, centro-temporal and centro-occipital channels, UMA and bandwidth changed with PMA. No differences in LMA, UMA and bandwidth were found between hemispheres. Skewness of LMA values strongly correlated with PMA, with positive skewness related to an immature brain (PMA=32 weeks) and negative skewness to a maturing (PMA>32 weeks) brain. In conclusion, we detected symmetric increase of aEEG characteristics, indicating symmetric brain maturation of the left and right hemispheres. In Chapter 3 we focused on automated power spectral analysis of the EEG in preterm infants, to identify maturation related changes of spectral measures. Weekly 4-hour EEG recordings (10–20 montage) were performed in 18 preterm infants with GA <32 weeks and normal neurological follow-up at 2 years of age, resulting in 79 recordings studied from 27+4 to 36+3 weeks PMA. Automated spectral analysis was performed on the 4-h EEG recordings. The frequency spectrum was divided in d1 (0.5-1 Hz), d2 (1-4 Hz), ¿ (4-8 Hz), a (8-13 Hz) and ß (13-30 Hz) band. Absolute and relative power of each frequency band and spectral edge frequency were calculated. Maturational changes in spectral measures were observed most clearly in the centro-temporal channels. With advancing PMA, absolute powers of d1, d2 and ¿ decreased. With advancing PMA, relative power of d1 decreased and relative powers of a and ß increased, respectively. In conclusion, with maturation, spectral analysis of the EEG showed a decrease in total power and a significant shift from the lower to the higher frequencies. In Chapter 4 we studied limitations of the EEG in a clinical environment. Assessment of the neonatal EEG may be hampered by drug-specific changes in electrocortical activity. To quantify effects of a loading dose of midazolam and lidocaine on the EEG frequency spectrum of full-term neonates with perinatal arterial ischemic stroke, 12 full-term infants with multi-channel aEEG and EEG recordings were retrospectively selected. During recording, midazolam and/or lidocaine were administered as anti-epileptic drug. We performed spectral analysis on 4-hour EEG segments around the loading dose. The frequency spectrum was divided in d (1-4 Hz), ¿ (4-8 Hz), a (8-13 Hz) and ß (13-30 Hz) bands. Midazolam induced immediate suppression of the aEEG background pattern for 30-60 minutes. Spectral EEG analysis showed decreased total and absolute frequency band powers. Relative d power decreased, ¿ power increased while a and ß powers remained constant. Lidocaine induced no aEEG background pattern suppression. Total and absolute EEG band powers were unchanged. Relative d power decreased, ¿ and a power increased and ß power remained constant. Effects were more pronounced in the stroke-affected hemisphere. In conclusion, both drugs induced a shift from low to higher frequency electrocortical activity. Additionally, midazolam reduced total EEG power. These spectral changes differ from those seen in adult studies. The Second part of this thesis focuses on novel methods for analysis of cortical function, i.e. quantifying changes in cortical function related to maturation and pathological condition through extraction and analysis of EEG transients, and through localization of the electrical sources of cortical activity. In Chapter 5, a method for automated detection of bursts, inter-burst-intervals (IBI), and continuous patterns in the EEG was developed in MATLAB® to aid with prognosis and stratification of clinical treatment for preterm infants. Results were evaluated for preterm infants with normal neurological follow-up at 2 years. The detection algorithm for burst, IBI and continuous pattern is based on selection by amplitude, time span, number of channels and numbers of active electrodes. Annotations of two neurophysiologists were used to determine threshold values. The training-set consisted of EEG recordings of 4 preterm infants with PMA of 29-34 weeks. Optimal threshold values were based on overall highest sensitivity. For evaluation, both observers verified detections in an independent dataset of 4 EEG recordings with comparable PMA. Algorithm performance was assessed by calculation of sensitivity and positive predictive value. Results of algorithm evaluation: sensitivity values of 90 ± 6, 80 ± 9 and 97 ± 5 % for burst, IBI and continuous patterns, respectively. Corresponding positive predictive values were 88 ± 8, 96 ± 3 and 85 ± 15 %, respectively. In conclusion, the algorithm showed high sensitivity and positive predictive values for bursts, IBI and continuous patterns in preterm EEG. In Chapter 6 we quantified the EEG burst frequency spectrum of preterm infants by automated analysis and described the topography of maturational change in spectral parameters. With this new method, we aim to increase sensitivity for maturational changes in cortical activity transients compared to non-specific analysis as presented in Chapters 2 and 3. Therefore, 18 preterm infants <32 weeks gestation and normal neurological follow-up at 2 years, underwent weekly 4-hour EEG recordings (10-20 system). The recordings (n=77) represent a large variability in PMA (28 - 36 weeks). We applied the automated burst detection algorithm presented in Chapter 5, and subsequently performed spectral analysis on burst activity. The frequency spectrum was divided in d1 (0.5-1 Hz), d2 (1-4 Hz), ¿ (4-8 Hz), a (8-13 Hz) and ß (13-30 Hz) band. Spectral parameters were evaluated as function of PMA by regression analysis. Results were interpolated and topographically visualized. The majority of spectral parameters show significant change with PMA. Highest correlation is found for d and ¿ band. Absolute band powers decrease with increasing PMA, while relative a and ß powers increase. Maturational change is largest in frontal and temporal region. In conclusion, the topographic distribution of maturational changes in spectral parameters corresponds with studies showing ongoing gyration and postnatal white matter maturation in frontal and temporal lobes. In Chapter 7 it was hypothesized that combining EEG with MRI could give more information about the development of possible brain damage. The study consisted of continuous EEG monitoring of neonates at risk of hypoxia-ischemia (HI), alternated by MRI examinations on predefined time points. Included were 8 term newborns with suspected neonatal encephalopathy caused by perinatal HI. The EEG was measured between day 1 and 8 after birth. MRI data was acquired between day 1 and 9, within 24 hrs from an observed EEG event. Preliminary results indicated an agreement between abnormalities in diffusion weighted MR images and EEG signal changes at corresponding electrode locations. In Chapter 8 we extended the preliminary results presented in Chapter 7 by exploring the feasibility of locating the sources of cortical activity as measured by scalp EEG in neonates. Three term neonates, representative examples of perinatal asphyxia and stroke, were selected for swLORETA analysis. A head model was created based on generic MRI data. Resulting 3-dimensional current source distributions showed good conformity with neonatal MRI images obtained within a 24-hour period from the EEG measurements. This visualization method may provide valuable information for longitudinal monitoring of neonates at risk, for guiding in-detail MRI examinations or for assessing treatment efficacy. The Third part of this thesis focuses on models for analysis of sub-cortical function, and describes modeling of the autonomic cardiovascular regulation of a preterm infant and assessment of the influence of asphyxia on autonomic cardiovascular regulation. In Chapter 9 we presented a model of the autonomic cardiovascular system of a preterm infant of 28 weeks gestation, with a birth weight of 1000 g and a closed ductus arteriosus by the end of the first week. This model describes the complex interactions between heart rate, blood pressure and respiration. The hemodynamic part of the model consists of a pulsatile heart and several vascular compartments, and is regulated by a baroreflex control system. The model is relatively simple to allow for a mathematical analysis of the dynamics, yet sufficiently complex to provide a realistic representation of the underlying physiology. Baroreflex related fluctuations were calculated by transfer function analysis between SBP and R-R fluctuations, using fast Fourier transform, estimated in the low frequency (LF, 0.04-0.15 Hz) band. LF transfer gain (baroreflex sensitivity, BRS, ms/mmHg) and phase (s) were assumed to reflect the baroreflex mediated cardiovascular fluctuations. The model provides (beat-to-beat) values of R-R interval and blood pressure that resemble realistic signals of preterm infants, and is validated with experimental heart rate and blood pressure data obtained in preterm infants. In Chapter 10 we monitored longitudinal changes of baroreflex activity in a fetal sheep model by analysis of the transfer function between low-frequency R-R and systolic blood pressure (SBP) fluctuations after umbilical cord occlusion (UCO) induced asphyxia. Fourteen preterm fetal lambs (UCO, n=5 and controls, n=9) were instrumented at 102 days gestation (term: 146 days) during caesarian section with a femoral arterial catheter, amniotic pressure catheter, 3-lead ECG and vascular occluder around the umbilical cord (UCO group). At 4 days after surgery, intra-uterine asphyxia was achieved by rapidly filling the vascular occluder for 25-min. Baroreflex related fluctuations were calculated by transfer function analysis. Before occlusion, no differences in R-R interval, SBP, and transfer gain and phase were observed between UCO and controls. Median BRS and delay were 7.8 (IQR, 5.6-11.4) and 4 (IQR, 3-6), and 6.2 (IQR, 5.6–8.4) ms/mmHg and 6 (IQR, 3-7) s, respectively. Three days after UCO, BRS gradually decreased in the UCO group, resulting in significantly lower values by the end of the week (UCO median 4.8, IQR 3.1-6.7 versus controls 7.7, IQR 6.7-10.6 ms/mmHg). In UCO group, delay increased to a median value of 7 (IQR, 4-8) s and differed from the controls. In conclusion, intra-uterine asphyxia results in long-lasting effects of baroreflex mediated fluctuations, limiting the possibility to buffer changes in SBP by adapting fetal heart rate. Finally, in Chapter 11 the results of all studies within this thesis are discussed and further research directions are proposed

Validation of a preterm infant cardiovascular system model under baroreflex control with heart rate and blood pressure data

In this paper we present an autonomic cardiovascular model of a preterm infant of 28 weeks of gestation with a birth weight of 1000 g and a closed ductus arteriosus by the end of the first week, that is capable of describing the complex interactions between heart rate, blood pressure and respiration. The hemodynamic model consists of a pulsatile heart and several vascular compartments, and is regulated by a baroreflex control system. The model is relatively simple to allow for a mathematical analysis of the dynamics but sufficiently complex to provide a realistic representation of the underlying physiology. The model provides (beat-to-beat) values of R-R interval and blood pressure that resemble realistic signals of preterm infants. The model is validated with experimental data obtained in preterm infants

Validation of a preterm infant cardiovascular system model under baroreflex control with heart rate and blood pressure data

\u3cp\u3eIn this paper we present an autonomic cardiovascular model of a preterm infant of 28 weeks of gestation with a birth weight of 1000 g and a closed ductus arteriosus by the end of the first week, that is capable of describing the complex interactions between heart rate, blood pressure and respiration. The hemodynamic model consists of a pulsatile heart and several vascular compartments, and is regulated by a baroreflex control system. The model is relatively simple to allow for a mathematical analysis of the dynamics but sufficiently complex to provide a realistic representation of the underlying physiology. The model provides (beat-to-beat) values of R-R interval and blood pressure that resemble realistic signals of preterm infants. The model is validated with experimental data obtained in preterm infants.\u3c/p\u3

EEG source localization in full-term newborns with hypoxic-ischemia

\u3cp\u3eThe aim of this study was to evaluate EEG source localization by standardized weighted low-resolution brain electromagnetic tomography (swLORETA) for monitoring of fullterm newborns with hypoxic-ischemic encephalopathy, using a standard anatomic head model. Three representative examples of neonatal hypoxic-ischemia were included. The method was validated with MRI data. Hypoxic-ischemic areas, visible on MRI, correlated well with swLORETA current density distributions. In addition, neonatal seizure activity may be localized. The calculated current density distributions provide easy-to-interpret localized information about neonatal brain function, which may enable detailed longitudinal monitoring and potential assessment of treatment efficacy.\u3c/p\u3

Generieke verdiepende RI&E MRI, voor werken met MRI: Een generieke risico-inventarisatie en -evaluatie voor werken met MRI in ziekenhuizen en Universitair Medische Centra in het kader van Arbowetgeving elektromagnetische velden

Dit rapport heeft als doel het opstellen van een generieke verdiepende risico-inventarisatie en -evaluatie (RI&E) van toepassingen van magnetic resonance imaging (MRI) ten behoeve van de volksgezondheid in algemene ziekenhuizen en universitair medische centra. Aanleiding voor deze RI&E is een wijziging van het arbeidsomstandighedenbesluit, een uitwerking van Europese richtlijn 2013/35/EG, in juli 2016. Voor de opzet van de RI&E is gekeken naar specifieke risico’s van werken met MRI, vergelijkbaar met de huidige arbocatalogus MRI van de NFU. Elektromagnetische velden (EMV) vormen de directe aanleiding voor deze RI&E. Teneinde de aanwezigheid van EMV bij toepassing van MRI toe te lichten, zijn de basisprincipes van MRI kort uitgelegd. Hieruit is ook de verdeling van EMV over de werkruimtes en de relatie van EMV tot de actieniveaus en grenswaarden bepaald. Voor de verschillende groepen werknemers die met MRI werken is, ten gevolge van grote inter-subject variatie, in globale termen een bepaling van het niveau van blootstelling aan EMV gemaakt, met een inschatting van frequentie en duur ervan. De (indirecte) effecten van EMV op de gezondheid zijn beschreven. Daarnaast zijn ook de risico’s van projectielwerking van ferromagnetische voorwerpen, geluid, werken met cryogene vloeistoffen en ethische risico’s benoemd. Voor al deze risico’s worden beheersmaatregelen beschreven, waarmee het gewenste beschermingsniveau voor medewerkers bereikt kan worden. Deze beheersmaatregelen worden gepresenteerd per risico en zijn geordend volgens de RI&E systematiek: afscherming van de bron, collectieve maatregelen, individuele maatregelen en persoonlijke bescherming. Voor kwetsbare werknemers, zijnde medewerkers met implantaten en zwangere medewerkers, worden specifieke beheersmaatregelen benoemd. Een voorbeeld hiervan is screening van medewerkers op contra-indicaties vóór betreding van de gecontroleerde zone. Een zwangere werknemer mag de gecontroleerde ruimte van een MRI scanner alleen betreden als er op dat moment niet gescand wordt, uitzonderingen op deze regel is mogelijk in overleg met de MR veiligheidsdeskundige. Dit beleid is gebaseerd op voorzorg, niet op bewezen schadelijkheid. Om als instelling rond MRI een veilige werksituatie te handhaven voor zowel medewerkers als patiënten is het van belang om te zorgen dat helder is wie verantwoordelijk is voor beleidsvorming en voor uitvoering ervan op de werkvloer, en dat lokaal expertise aanwezig is. Vanwege noodzakelijke gespecialiseerde kennis is het wenselijk om in lijn met het standpunt van de Europese federatie van organisaties van klinische fysica (EFOMP) hiertoe specifieke MR functionarissen te benoemen: een toezichthoudend medewerker MR veiligheid (MR safety officer) en een MR veiligheidsdeskundige (MR safety expert). Voor deze functionarissen worden taken beschreven en vereisten op gebied van kennis en opleiding benoemd. Door het aanstellen van deze functionarissen wordt de in de Arbowet vereiste toegang tot een deskundige die kan bijstaan bij het uitvoeren van de RI&E en het plan van aanpak gerealiseerd. Tenslotte is beschreven hoe de omstandigheden de overschrijding van grenswaarden uit het Arbobesluit EMV rechtvaardigen voor medewerkers bij de MRI. Hiervoor worden, voor de verschillende grenswaarden, de procedures benoemd waarbij er risico bestaat dat deze grenswaarden overschreden worden. Deze procedures komen af en toe voor, en vormen slechts een fractie van de MRI-gerelateerde werkzaamheden die medewerkers uitvoeren. Deze verdiepende RI&E laat zien dat overschrijdingen van de grenswaarden bij gebruik van MRI beperkt zijn en dat de risico’s voor medewerkers met de beschreven beheersmaatregelen zodanig zijn dat de werknemers bij gebruik van MRI goed beschermd zijn. Deze RI&E kan de basis vormen voor een beleidsdocument voor afdelingen, arbodiensten en medezeggenschapsraden van ziekenhuizen waar met MRI wordt gewerkt. Ook kan het de basis vormen voor een update van de huidige arbocatalogus MRI van de NFU en aanpassing van de arbocatalogus niet-ioniserende straling van de NVZ

EEG based time points for longitudinal MRI studies in neonates with hypoxia-ischemia

\u3cp\u3eThis study consisted of continuous EEG monitoring of neonates at risk of hypoxia ischemia (HI), alternated by MRI examinations on predefined time points. It is hypothesized that combining EEG with MRI could give more information about the development of possible brain damage. Included were 8 term newborns with suspected neonatal encephalopathy caused by perinatal HI. The EEG was intermittently measured between day 1 and 8 after birth. MRI data was acquired between day 1 and 9, within 24 hrs from an observed EEG event. Preliminary results indicate an agreement between abnormalities in diffusion weighted MR images and signal fluctuations at corresponding electrode locations.\u3c/p\u3

Automatic detection of burst synchrony in preterm infants

\u3cp\u3eElectroencephalographic characteristics are useful in assessment of the functional status of specific neuronal connections relative to postmenstrual age. Interhemispheric burst synchrony (IBS) is a measure of the functional connectivity between the hemispheres in the maturing preterm brain. An algorithm was developed to assess IBS and was used in a prospective, longitudinal EEG study on 18 very preterm infants (&lt; 32 weeks gestational age) with normal follow-up at 2 years of age. The preterm infants underwent weekly 4-hour multi-channel EEG recordings, resulting in n = 77 EEGs. After automated detection of bursts, the algorithm defines the start and end of interhemispheric synchronous burst activity, based on selection criteria found in literature. The algorithm was designed to emulate visual inspection, providing objective results in an automated manner. This approach may be applied in clinical use and open novel avenues to automated analysis in EEG monitoring and, moreover, it may facilitate assessment of the functional status of interhemispheric connections. As such, assessment of low interhemispheric synchrony may be associated with brain injury.\u3c/p\u3

Functional connectivity in preterm infants derived from EEG coherence analysis

Objective : To quantify the neuronal connectivity in preterm infants between homologous channels of both hemispheres. Methods : EEG coherence analysis was performed on serial EEG recordings collected from preterm infants with normal neurological follow-up. The coherence spectrum was divided in frequency bands: dnewborn(0–2 Hz), ¿newborn(2–6 Hz), anewborn(6–13 Hz), ßnewborn(13–30 Hz). Coherence values were evaluated as a function of gestational age (GA) and postnatal maturation. Results : All spectra show two clear peaks in the dnewborn and ¿newborn-band, corresponding to the delta and theta EEG waves observed in preterm infants. In the dnewborn-band the peak magnitude coherence decreases with GA and postnatal maturation for all channels. In the ¿newborn-band, the peak magnitude coherence decreases with GA for all channels, but increases with postnatal maturation for the frontal polar channels. In the ßnewborn-band a modest magnitude coherence peak was observed in the occipital channels, which decreases with GA. Conclusions : Interhemispherical connectivity develops analogously with electrocortical maturation: signal intensities at low frequencies decrease with GA and postnatal maturation, but increase at high frequencies with postnatal maturation. In addition, peak magnitude coherence is a clear trend indicator for brain maturation. Significance : Coherence analysis can aid in the clinical assessment of the functional connectivity of the infant brain with maturation

Effects of midazolam and lidocaine on spectral properties of the EEG in full-term neonates with stroke

\u3cp\u3eAssessment of the neonatal EEG may be hampered by drug-specific changes in electrocortical activity. To quantify effects of a loading dose of midazolam and lidocaine on the EEG frequency spectrum of full-term neonates with perinatal arterial ischemic stroke (PAIS), 11 full-term infants underwent multi-channel amplitude-integrated EEG (aEEG) and EEG recordings. During recording, midazolam and/or lidocaine were administered as anti-epileptic drug. Retrospectively, we performed spectral analysis on 4-h EEG segments around the loading dose. The frequency spectrum was divided in δ (1-4 Hz), θ (4-8 Hz), α (8-13 Hz) and β (13-30 Hz) bands. Midazolam induced immediate suppression of the aEEG background pattern for 30-60 min. Spectral EEG analysis showed decreased total and absolute frequency band powers. Relative δ power decreased, θ power increased while α and β powers remained constant. Lidocaine induced no aEEG background pattern suppression. Total and absolute EEG band powers were unchanged. Relative δ power decreased, θ and α power increased and β power remained constant. Effects of lidocaine were more pronounced in the stroke-affected hemisphere. In conclusions, both drugs induced a shift from low to higher frequency electrocortical activity. Additionally, midazolam reduced total EEG power. These spectral changes differ from those seen in adult studies.\u3c/p\u3