Measuring Near-Infrared Spectroscopy Derived Cerebral Autoregulation in Neonates: From Research Tool Toward Bedside Multimodal Monitoring.

Introduction: Cerebral autoregulation (CAR), the ability of the human body to maintain cerebral blood flow (CBF) in a wide range of perfusion pressures, can be calculated by describing the relation between arterial blood pressure (ABP) and cerebral oxygen saturation measured by near-infrared spectroscopy (NIRS). In literature, disturbed CAR is described in different patient groups, using multiple measurement techniques and mathematical models. Furthermore, it is unclear to what extent cerebral pathology and outcome can be explained by impaired CAR. Aim and methods: In order to summarize CAR studies using NIRS in neonates, a systematic review was performed in the PUBMED and EMBASE database. To provide a general overview of the clinical framework used to study CAR, the different preprocessing methods and mathematical models are described and explained. Furthermore, patient characteristics, definition of impaired CAR and the outcome according to this definition is described organized for the different patient groups.
Results: Forty-six articles were included in this review. Four patient groups were established: preterm infants during the transitional period, neonates receiving specific medication/treatment, neonates with congenital heart disease and neonates with hypoxic-ischemic encephalopathy (HIE) treated with therapeutic hypothermia. Correlation, coherence and transfer function (TF) gain are the mathematical models most frequently used to describe CAR. The definition of impaired CAR is depending on the mathematical model used. The incidence of intraventricular hemorrhage in preterm infants is the outcome variable most frequently correlated with impaired CAR. Hypotension, disease severity, dopamine treatment, injury on magnetic resonance imaging (MRI) and long term outcome are associated with impaired CAR. Prospective interventional studies are lacking in all research areas. Discussion and conclusion: NIRS derived CAR measurement is an important research tool to improve knowledge about central hemodynamic fluctuations during the transitional period, cerebral pharmacodynamics of frequently used medication (sedatives-inotropes) and cerebral effects of specific therapies in neonatology. Uniformity regarding measurement techniques and mathematical models is needed. Multimodal monitoring databases of neonatal intensive care patients of multiple centers, together with identical outcome parameters are needed to compare different techniques and make progress in this field. Real-time bedside monitoring of CAR, together with conventional monitoring, seems a promising technique to improve individual patient care.

Introduction

The transitional period in neonates is an extremely vulnerable phase prone to hemodynamic instability (i.e., hypotension, cyanosis, shock, ischemia, reperfusion injury) at high risk for cerebral ischemic and/or hemorrhagic lesions. At the same moment interventions such as intubation, surfactant administration, cooling in hypoxic-ischemic encephalopathy (HIE), treatment with inotropes, surgery and sedation can increase or decrease this risk. Maintaining adequate brain tissue oxygenation, a stable balance between cerebral oxygen delivery and extraction, is one of the major goals in neonatology because brain damage by ischemic or hemorrhagic lesions will often lead to impaired neurodevelopmental outcome (1). Ninety percent of the intraventricular hemorrhages (IVH) occur in the first 72 h after birth, suggesting that brain circulation is especially vulnerable in this period (2). To monitor and improve brain tissue oxygenation, arterial blood pressure (ABP) is often used as a surrogate measurement for cerebral blood flow (CBF). However, this classical early-goal directed therapy -increasing ABP if below a cut-off value- may not be adequate nor ideal in the healthy preterm or sick neonate in the transitional phase (3). Several mechanisms play a role in maintaining cerebral oxygenation, ABP being only one factor in this very complicated physiological system (4) (Figure 1, with permission). In this introduction, the determinants of cerebral oxygen delivery are discussed. Then, the concept of static and dynamic cerebral autoregulation (CAR) is explained with an emphasis on its assessment using near-infrared spectroscopy (NIRS).

Figure 1

Model for brain circulation. Cerebral autoregulation is only a single element in the interaction between blood processes and vascular smooth muscle processes. Multiple factors play a role in brain hemodynamics regulation. This figure shows the complex interaction between the different processes in brain circulation. Furthermore, techniques to study these processes in a non-invasive way are reported (with permission (4)). EEG, (amplitude-integrated) electroencephalogram; CSF, cerebrospinal fluid; fMRI, functional magnetic resonance imaging; NIRS, near-infrared spectroscopy; PET, positron emission tomography.

Cerebral oxygen delivery is determined by CBF and blood oxygen content. CBF itself is the result of the gradient between cerebral perfusion pressure (CPP) and cerebrovascular resistance (CVR). CPP is determined by ABP and intracranial pressure (ICP).

CVR reflects the varying tone of smooth muscle cells in the wall of arteries. One of the important factors influencing this smooth muscle tone is ABP. The myogenic reflex will cause the vessel to constrict or dilate if CPP increases or decreases, respectively. This reflex leads to the classical concept of CAR, the mechanism in which CBF is maintained stable regardless of changes in ICP and first described in humans by Lassen (5). If ICP is stable, CPP can be replaced by ABP. In this way, changes in CBF can be measured for a range of ABP values to determine CAR. However, multiple factors, apart from changes in ABP, can influence smooth muscle tone. The most relevant is the chemical influence of pCO2 on the muscle tone, known as CO2 vaso-reactivity, but also effects of NO, calcium and physical stimuli have been described (6, 7). Furthermore, different tissue processes (i.e., functional activation, autonomic neural activity, among others) cause alteration in regional CBF. This is known as blood flow metabolism coupling and is more extensively described in the paper by Huneau et al. (8). In this review article we will focus on the classical definition of flow-pressure CAR.

A good overview of the first flow-pressure CAR studies is provided by Weindling et al and Greisen (9, 10). In this review it is indicated that static CAR can be studied as a steady-state response: looking at the relationship between CBF and CPP (ABP) without considering the time course of changes in flow following changes in pressure. Impaired CAR (CBF being strongly correlated with ABP) in preterm infants is described using 133Xe clearance or plethysmography after which it was hypothesized that loss of CAR plays a decisive part in the pathogenesis of brain lesions in the neonate (11–14). To study dynamic flow-pressure CAR, the continuous measurement of changes in CPP, and thus ABP, is mandatory. Spontaneous rhythmic oscillations in ABP and CBF are used to describe the CAR mechanisms within different frequency bands. The low and very low frequency bands are dominated by myogenic, neurogenic and metabolic regulatory factors. Sometimes, the lack of a large variation in CPP and CBF makes it necessary to challenge the regulatory mechanisms in order to be able to measure autoregulation reliably. Therefore, some studies have induced changes in ABP (15) while other use the changes in ABP due to pathology, in order to be able to assess CAR. This might seem as a limiting factor for the continuous monitoring of CAR. However, the more stable the clinical situation, the less problems regarding CAR are expected.

The dynamical aspect of CAR can be measured using Doppler flowmetry or NIRS. Doppler flowmetry uses ultrasound to measure blood flow velocity in a specific vessel, measuring the immediate effect of changes in blood pressure on the CBF velocity, while NIRS reflects the effect of changes in ABP over a longer period of time. For a further review of the Doppler technique to measure CAR, we refer to Panerai et al. (16). With NIRS, cerebral oxygen saturation and cerebral fractional tissue oxygen extraction (cFTOE) can be measured (17, 18). NIRS is based on the specific absorption of near-infrared light by oxyhemoglobin (HbO2) and deoxyhemoglobin (HHb) in cerebral tissue, reflecting the mixed venous-capillary-arterial oxygenation state of the parenchyma directly beneath the light emitting sensor. It has been shown that under constant arterial oxygen saturation (SaO2) and a constant brain metabolism, NIRS derived cerebral oxygen saturation can be used as a surrogate measurement for CBF. This was first described by Tyszczuk et al, describing normal CBF even during periods of low ABP (19).

The use of NIRS to study CAR was first documented by Tsuji et al, using the correlation between the hemoglobin difference (HbD) and mean ABP, in periods of stable oxygen saturation, in preterm infants. They were the first to correlate impaired CAR with brain lesions, mainly IVH, using the classical NIRS technology with the Beer-Lambert law (20). The introduction of spatially resolved spectroscopy, using different receptors in the light emitting sensor, led to measurements less prone to movement artifacts and thus more easy to use in clinical situations. Wong et al. were the first to describe impaired CAR, assessed by means of the brain tissue oxygenation index (TOI) in the sickest infants with the highest clinical risk index for babies (CRIB) scores (21).

The dynamic approach might become the reference method for clinical assessment of CAR. However, the major disadvantage is the lack of standardization in the measurement. The dynamics in time and magnitude between changes in ABP and NIRS derived cerebral oxygenation can be assessed by different mathematical models. These models range from the simple time-domain correlations between cerebral oxygenation and ABP, which assumes autoregulation is a simple linear process, to more complex techniques, based on continuous wavelet transforms (CWT), to describe the stochastic, non-stationary dynamic nature of CAR.

Since measurement of cerebral oxygenation is easy and non-invasive in neonates, several groups have been studying the added value of measuring CAR in different patient groups, using different NIRS instruments, different mathematical models and different outcome parameters in the last 2 decades. The aim of this review is to provide a general overview of the clinical framework to study CAR by describing and explaining the different preprocessing methods and mathematical models. Furthermore, we will critically summarize the available literature.

Methods

Between July 2015 and 20 November 2017, we conducted an extensive search of the literature to identify clinical studies measuring NIRS-derived CAR in neonates. This search was conducted using the PRISMA checklist (http://prisma-statement.org) in the PUBMED and EMBASE database. The search strategy included the terms (cerebrum OR cerebral OR brain) AND (homeostasis OR autoregulation) AND (infant, newborn OR infant OR newborn infant OR neonates) AND (spectroscopy, near-infrared OR spectroscopy AND near-infrared OR NIRS OR near infrared spectroscopy OR near AND infrared AND spectroscopy) without language restriction. Inclusion criteria were (1) clinical study in (2) neonates measuring quantitatively (3) flow-pressure CAR using (4) NIRS methodology. We excluded studies or parts of certain methods/results measuring other parameters (i.e., heart rate, cardiac output) to evaluate CAR. Also, abstracts only, case reports, and articles limited to methodological questions without short or long term outcome were excluded. LT screened the article titles and abstracts to determine whether they met the inclusion criteria. Then, LT reviewed full text articles to assess for eligibility. Any articles presenting doubts or inconsistencies were fully reviewed by GN and AC until a decision was reached on their inclusion or exclusion. By using cross-references, additional eligible articles were added by hand searching.

We collected study type and data about patient numbers, gestational age (GA), postmenstrual age (PMA) or day of life (DOL) at start of the study. Furthermore, data about the mathematical model, initial sample frequency of data extraction, epoch length and duration of measurement per patient was extracted. The proposed definition of impaired CAR and the outcome determined by the authors was defined. Different patient groups according to pathology were identified to summarize the results.

Results

The result of the systematic search is presented in an adapted PRISMA flow diagram (Figure 2). Included studies are organized in 4 different patient groups: (1) preterm infants during the transitional period, (2) neonates receiving medication/treatment, (3) neonates with congenital heart diseases (CHD), and (4) neonates with HIE treated with therapeutic hypothermia. An overview of the different studies per patient group is provided in Table 1.

Figure 2

Adapted PRISMA flow diagram.

Table 1

Overview of included studies.

Author	-Patient number-GA/PMA in weeks in mean (SD) or median (range)/*DOL-Study type	NIRS instrument	-Mathematical model-Sample frequency-Epoch length-Mean or median duration	Definition impaired CAR	Outcome
PRETERMS DURING THE TRANSITIONAL PERIOD
Schat et al. (26)	-n = 28-case: 27.9 (26.3–34.7); control: 27.4 (25.6–34.7)-Case control	INVOS 5100C	-COR between MAP and FTOE-0.003 Hz-5 min-48 h	Significant negative COR coefficient	No COR between presence of impaired CAR and NEC
Riera et al. (34)	-n = 54-27 (1.9)-Cohort	NIRO-200nx	-(p)PD COH between MAP TOI in 0.003–0.04 Hz band-2 Hz-30 min-9.5 h	-Threshold PDCMAP>>TOI for low SVC flow was 0.554 = PDCMAP>>TOI classifier-Low SVC flow as surrogate for cerebral hypoperfusion	-PDCMAP>>TOI predicted low SVC flow-PDCMAP>>TOI classifier associated with % of time of MABP <GA-5 mmHg-pPDC MAP>>TOI predicted cardiovascular support, severe IVH and death
Vesoulis et al. (69)	-n = 62-25.4 (1.3)-Cohort	Foresight	-Log transformation of TF gain between MAP and SctO2 in 0.08–0.12 Hz band-0.5 Hz-20 min-68 h	Stronger dampening (more negative TF gain coefficient) is better autoregulation	-Greater dampening independently associated with advancing GA, BW and chorioamnionitis-Less dampening independently associated with African-American race, IVH during first week
Stammwitz et al. (70)	-n = 31-272/7 (261/7-322/7)-Cohort	Critikon Cerebral Oxygenation Monitor 2001	-COH between MAP and tHb/OI in 0–0.01 Hz band-2 Hz-≥12 min-8 h	No cut off. Hypothesis that high COH indicate a coordination of physiological sub-systems and thus are a sign of health	Low COH in the first 24 h were associated with IVH ≥ 3, death and MDI
Binder-Heschl et al. (39)	-n = 46-case: 33.4 (1.9); control: 33.3 (1.3)-Case control	INVOS5100 C	-COR between (invasive n = 3) MABP and crSO2-0.13 Hz-1 h-24 h	No definition	Very weak COR between MABP and crSO2 suggesting intact CAR during borderline hypotension
Eriksen et al. (67)	-n = 60-26.6 (1.3)-Cohort	NIRO 300	-COx (moving linear COR), regression coefficient vs. COH, TF gain between MAP and OI in 0.003–0.04 Hz band-2 Hz-10 min-2.3 h	COx ≥ 0.4 and COH ≥ 0.5	-COR between TF gain and regression coefficient was weak (r = 0.245) but significant after exclusion of outliers-COx is more robust, TF gain also increases if MAP and OI are in counterphase
Riera et al. (72)	-n = 54-27 (1.9)-Cohort	NIRO 200 NX	-(p)BiAR-COH, (p)COH between MABP and TOI in 0.003–0.04 Hz band-2 Hz-30 min-9.5 h	-Threshold BiAR-COH and COH for low SVC flow was 0.58 and 0.52 respectively-Low SVC flow as surrogate for cerebral hypoperfusion	-BiAR-COH better in predicting low SVC flow with compared to COH,-pBiAR-COH but not pCOH was associated with IVH grade 3–4 and predicted mortality
Verhagen et al. (30)	-n = 25-29.1 (25.4–31.7)-Cohort	INVOS 4100–5100	-COR between MABP and rcSO2/FTOE-0.003 Hz-5 min-24 h	Statistically significant positive and negative COR between rcSO2-MABP and FTOE/MABP, respectively	Identification of absent CAR in 40% of patients, no correlations between absent CAR and clinical variables except higher hemoglobin levels
Caicedo et al. (71)	-n = 9- <32-Case control	INVOS 4100–5100	-BPRSA between MABP and rScO2-1 Hz-?-?	No definition	Presence of non-linear relations between the variables. In addition, the BPRSA curves from the control subjects converge faster to zero than the curves for the subjects with IVH gr 3–4.
Alderliesten et al. (42)	-n = 90-246/7-31-Case control	INVOS 4100–5100	-COR between MABP and rScO2-1 Hz-1 min-Cases: 48 h; controls: 47 h	COR > 0.5	-More time with impaired autoregulation before and after detection of PIVH compared to controls
Hahn et al. (66)	-n = 60-27 (1.3)-Cohort	NIRO 300	-COH, TF gain between MAP and OI in 0.003–0.04 and 0.04–0.1 Hz band-2 Hz-10 min-2.3 h	COH ≥ 0.45–0.47	-Negative association between TF gain and MAP-No association between impaired CAR and antenatal or postnatal signs of inflammation, IVH or mortality
Caicedo et al. (63)	-n = 42-28.1 (2.27)-Cohort	Criticon Cerebral RedOx monitor	-TF gain, phase between MABP and HbD in 0.003–0.02, 0.02–0.05 and 0.05–0.1 Hz band-0.333 Hz-20 min-72 h	No definition	-Significant higher TF gain in normal compared to abnormal (IVH, PVL, death, abnormal MDI and/or PDI) population in 0.05–0.1 Hz band
Wong et al. (73)	-n = 32-26.3 (1.5)-Cohort	NIRO 200	-COH and TF gain between MABP and TOI in 0.003–0.02 Hz band-6 Hz-20 min-5 × 20 min on day 1-2-3	COH ≥0.5	-High COH and TF gain at low BPV in unstable children with brain injury.-Significant association between maximum COH and BPV in stable children.
Zhang et al. (60)	-n = 17-26 (2)-Cohort	NIRO 300	-COH, TF gain and phase between MAP and HbO2/HHb/HbD/TOI in 0.02–0.04, 0.04–0.15, 0.15–0.25 Hz band-1 kHz-10 min-10 min	COH ≥ 0.5 or max COH if <0.5	-Multiple testing-Strongest relation was found between COH MAP-HHb in 0.04–0.15 Hz band and CRIB-II
Caicedo et al. (57)	-n = 33 and 20-28.9 (1.8) and 28.4 (3.5)-Cohort	NIRO 300 and INVOS 4100	-COR, (PA)COH, between MABP and HbD/TOI/rSO2 in 0.003–0.1 Hz band-0.333 Hz-20 min-50–70 and 6–9 h	COR/(PA)COH>0.5CPRT COR/(PA)COH: % time with impaired CAR	No significant correlation with CRIB/MDI/PDI/Griffith score and mean or CPRT COR/(PA)COH
Gilmore et al. (44)	-n = 23-26.7 (1.4)-Cohort	Foresight	-COx (moving linear COR) between MABP and SctO2-0.5 Hz-5 min-3.2 days	COx > 0.5	Impaired autoregulation was associated with low MABP but not with IVH.
Hahn et al. (56)	-n = 22-27.5 (24.1–29.4)-Cohort	NIRO 300	-COH between MAP and OI in 0.003–0.04 and 0.04–0.1 Hz band-2 Hz-10 min-2.1 h	Threshold COH with simulationCOHST = COH minus ThresholdCOHCOHST ≥ 0 implies impaired CAR	-Precision of COH to measure CAR is improved when the magnitude of variability in ABP is taken into account
De Smet et al. (32)	-n = 10, 10 and 10-281/7 (21/7), 292/7(12/7) and 285/7 (32/7)-Cohort	Critikon Cerebral Oxygentation monitor 2001, INVOS 4100 and NIRO 300	-(PA)COH between MABP and HbD/rSO2/TOI in 0.0033–0.04 Hz band-1.677, 1 and 10 Hz-10.15 and 12.5 min-72 h	(PA)COH>0.5PPI: % epochs with impaired CARCPRT: % time with impaired CAR	-High PCOH values are better indicators of poor clinical outcome (MDI <84, PDI <84, Apgar <7) than COH-CPRT and PPI are better indicators of poor clinical outcome than mean score values
O'Leary et al. (64)	-n = 88-26 (23–30)-Cohort	NIRO 500	-COH, TF gain between MAP and HbD in 0.05–0.25, 0.25–0.5 and 0.5–1.0 Hz bands-2 Hz-10 min-75.2 h	COH > 0.69	High TF gain was significantly associated with IVH in 0.05–0.25 Hz band
De Smet et al. (40)	-n = 20-28.7 (24–39)-Cohort	NIRO 300	-COR, (PA)COH between MAP and HbD/TOI in 0–0.01 Hz band-0.2 Hz-30 min-?	COR, (PA)COH > 0.5CPRT: % time with impaired CAR	-TOI may be used for the calculation of cerebral autoregulation.-CPRT generates a measurement of the autoregulation impairment proportional to COR and PA(COH)
Wong et al. (21)	-n = 24-26 (2.3)-Cohort	NIRO 300	-COH, TF gain between MAP and TOI in 0.003–0.02, 0.02–0.05, 0.05–0.1 Hz band-1 Hz-20 min-52 min	COH ≥ 0.5	-High COH and high TF gain were found in sickest infants in 0.003–0.02 Hz band-CRIB score best predictor of COH-COH ≥ 0.5 predictive for mortality
Soul et al. (59)	-n = 90-26.5 (23–30)-Cohort	NIRO 500	-COH between MAP and HbD in 0–0.04 Hz band-2 Hz-10 min-17.2 h	-COH ≥0.77-PPI: % epochs with impaired CAR	-Pressure passive cerebral circulation associated with GA and BW, hypotension, maternal hemodynamic factors.-Pressure passivity in 87/90 patients, with mean PPI of 20.3% (range 0–48.6)
Lemmers et al. (29)	-n = 38-case 28.6 (1.32); control 29.3 (1.74)-Case control	INVOS 4100	-COR between MAPB and ScO2/FTOE-10 Hz-15 min-420 min	-COR MABP/ScO2 > 0.5-COR MABP/FTOE <-0.5	More 15 min periods with impaired autoregulation in RDS in comparison with no RDS
Morren et al. (65)	-n = ?-?-Cohort	NIRO 300	-COR, COH, CPC between MAP and HbD in 0–0.01 Hz band-0.2 Hz-30 min-?	No definition	CPC and COR are better measures to detect impaired autoregulation than COH analysis.
Tsuji et al. (20)	-n = 32-27.1 (2.5)-Cohort	NIRO 500	-COH between MAP and HbD in 0–0.01, 0.01–0.05 and 0.05–0.1 Hz band-2 Hz-30 min-207 min	COH > 0.5	-Impaired autoregulation in 0–0.01 Hz band was observed in 53% of patients and in 80% of patients with IVH grade 3/4 or PVL
NEONATES RECEIVING MEDICATION/TREATMENT
Li et al. (75)	-n = 44-case:29.5 (1.3); control: 29.3 (1.6)-Case control	MC-2030CCerebral oximeter	-COR between MAP and ScO2-0.1 Hz-5 min-20 min	Deviation from baseline of COR coefficient suggests less effective CAR	Longer lasting impaired CAR with surfactant administration with INSURE compared to LISA method
Alderliesten et al. (74)	-n = 132-case: 292/7 (256/7-314/7); control: 293/7 (255/7-314/7)-Case control	INVOS 4100–5100	-COR between MABP and rScO2-1 Hz-15 min-72 h	% time with COR > 0.5	Impaired CAR associated with treatment with higher doses of dopamine compared to no blood pressure support
Eriksen et al. (45)	-n = 60-case: 26.2(1.5); control: 26.7(1.2)-Case control	NIRO 300	-COx (moving COR) between MAP and OI-2 Hz-10 min-2.3 h	COx > 0	Impaired CAR associated with dopamine treatment compared to no dopamine treatment
Baerts et al. (58)	-n = 18-case: 27.2 (252/7-294/7) control: 27.4 (250/7-295/7)-Case control	INVOS 4100–5100	-COR in the very slow frequency range (1/60 HZ) between MABP and rScO2-1 Hz-15 min-1 h	COR > 0.5 during 10% or more of time	No difference in CAR between offsprings of mothers treated with indomethacine and controls
Caicedo et al. (62)	-n = 56-29 (24.7–31.9)-Case control	INVOS 4100	-COR, COH and TF gain between MABP and rScO2 in 0.003–0.02, 0.02–0.05 and 0.05–0.1 Hz band-1 Hz-15 min-72 h	High TF gain	Higher TF gain in offsprings of mothers treated with labetalol during 1 day of life in 0.003–0.02 and 0.02–0.05 Hz band compared to controls
Kooi et al. (31)	-n = 14-26.7 w-Cohort	INVOS 5100 C	-COR between changes in MABP and cFTOE-?-10–30 min-2 h	Increase in MABP of 2 mmHg combined with decrease of cFTOE of 5%	Unable to define subgroup of infants lacking CAR after volume treatment
Papademetriou et al. (24)	-n = 6-*day 3–16-Interventional	Hitachi ETG-100	-CWT, WCC between MAP and HbO2 in 0.06–0.13, 0.13–0.25 and 0.25–1 Hz band-5 Hz-10 min-70 min	WCC > 0.5	-Loss of CAR at low ECMO flow-Right hemisphere more susceptible to low flow
Chock et al. (41)	-n = 40-26 (1)-Case control	INVOS 5100	-COR between MAP and rSo2-0.2 Hz-20 min-26 h	-COR > 0.5-PPI:% epochs with impaired CAR-max COR	-PPI was significantly higher 2 h after ductal ligation compared with control and indomethacin PDA treatment-Dopamine use was associated with max COR, independent of PDA treatment strategy
Wagner et al. (15)	-n = 24 (11 neonates)-neonates >36 weeks; other 2 month-15 year-Interventional	NIRO 500	-ΔHb_DiffΔMABP, ΔHb_TotalΔMABP after phenylephrine bolus-2 Hz-5–60 s-240 s	ARIHbDiff/ARIHbTotal >0	-Significant correlations between ARI using cerebral HB signals and direct CBF measures-Studying dynamic CAR using bolus phenylephrine is faster and has better signal-to-noise ratio compared to use of spontaneous blood pressure fluctuations
Munro et al. (53)	- n = 17-26 (0.4)-Case control	NIRO 500	-Linear regression of CBF vs. MAP-?-30 min-4 short measurements	Identification of breakpoint of MAP at which the residual sums of squares reaches a minimum	-CBF is autoregulated above 29 mmHg in extreme prematures-No evidence of autoregulation in dopamine treated infants
NEONATES WITH CONGENITAL HEART DISEASE
Smith et al. (51)	-n = 64-*day 9 (4–30)-Cohort	Reflectance spectroscopy monitor	-HVx (moving COR) between ABP and blood volume index-240 Hz-300 s-414 min	COR > 0	-Hypothermia was associated with hypotension, dysautoregulation and increased cerebral oximetry but collinearity between 3 variables during CPB
Votava-Smith et al. (76)	-n = 24-38.6 (0.9)-Cohort	FORE-SIGHT	-COH between MAP and SctO2 in 0.003–0.04 Hz band-0.5 Hz-20 min-23.4 h	COH > 0.58% of epochs/patient with abnormal CAR	-All subjects had epochs with impaired CAR, with mean 15.3% (3.5–56%) during first days of life-Lack of sedative medication, low hemoglobin, low MAP and greater FTOE were associated with impaired CAR
Brady et al. (43)	-n = 54-*25 months (0–222)-Cohort	INVOS	-COx (moving COR) between MABP and rSO2-60 Hz-0.25 Hz-300 s-?	COx > 0.4, sorting by MABP, to determine LLA	-Broad range of individual LLA during CPB-CPB was associated with both hypotension and impairment of CAR
Bassan et al. (55)	-n = 43-*day 9 (2–210)-Cohort	NIRO 500	-COH between ΔMAP and ΔHbD in 0–0.1 Hz band-2 Hz-5 min-105 min	-COH > 0.5-PPI: % epochs with impaired CAR	In early postoperative phase after CPB, higher end- tidal CO2 and higher MAP variability increased odds of impaired CAR
NEONATES WITH HYPOXIC-ISCHEMIC ENCEPHALOPATHY TREATED WITH THERAPEUTIC HYPOTHERMIA
Chavez-Valdez et al. (47)	-n = 75-386/7 (16/7)-Cohort	INVOS 5100	-HVx (moving COR) between MAP and rTHb-100 Hz-300 s-57.8 h	-HVx>0-% time, max deviation and AUC below MAPopt	Impaired CAR in HIE and therapeutic hypothermia correlated with cardiopulmonary injury and sex
Lee et al. (77)	-n = 64-*day 1-Cohort	INVOS 5100	-HVx (moving COR) between MAP and rTHb-100 Hz-300 s-58.1 h	-HVx > 0-% time, max deviation and AUC below MAPopt	Impaired CAR during and after therapeutic hypothermia correlated with neurologic injury on MRI
Tian et al. (25)	-n = 9-39(2)-Cohort	INVOS 4100–5100	-CWT between MAP and SctO2-0.033 Hz-not applicable- <72 h	Significant in-phase and anti-phase coherence between blood pressure and SctO2	Impaired CAR correlated with MRI severity score and clinical outcome
Tekes et al. (48)	-n = 27->35 w, *day 1-Cohort	INVOS5100	-HVx (moving COR) between MAP and rTHb-100 Hz-300 s-43 h	-HVx>0-% time, max deviation and AUC below MAPopt	Impaired CAR during hypothermia and rewarming correlated with ADC scalars in specific anatomic regions on MRI
Massaro et al. (68)	-n = 36-38.6 (1.7) and 39.2 (1.3)-Cohort	NIRO 200	-COH, TF gain between MAP and HbD in 0.05–0.25 Hz band-1 kHz-10 min-73 h	COH>0.384PPI: % epochs COH>0.384TF gain within pressure-passive epochs	Impaired CAR during hypothermia and rewarming correlated with MRI severity score or death
Burton et al. (49)	-n = 19-38.9 (1.5)-Cohort	INVOS 5100	-HVx (moving COR) between MAP and rTHb-100 Hz-300 s-43 h	-HVx>0-% time, max deviation and AUC below MAPopt	Impaired CAR during rewarming correlated with 2-year neurodevelopmental outcome
Howlett et al. (50)	-n = 24-39.2(1.5)-Cohort	INVOS5100	-HVx (moving COR) between MAP and rTHb-100 Hz-300 s-43.4 h	-HVx>0-% time and max deviation below MAPopt	Impaired CAR during rewarming correlated with MRI injury severity in specific anatomic regions

First, we will present the clinical framework to study CAR by describing the different preprocessing methods and mathematical models, based on the literature search and identified clinical studies. Secondly, the definition of impaired CAR is discussed. Finally, the different study outcomes organized by the 4 identified patient groups are summarized.

Clinical framework

The clinical framework with general setup for the assessment of CAR is presented in Figure 3. The patient is connected to different monitors in order to acquire the relevant data for CAR assessment (Figure 3A). The main goal of this monitoring system is to indicate the status of the CAR mechanisms (Figure 3B). Recent advances enable the use of multimodal monitoring technologies, where different signals can be acquired at the same time (Figure 3C). In this context, different surrogates for CPP and CBF can be obtained. Additionally, other systemic parameters, which may influence CAR assessment can be measured. When the signals are obtained, different preprocessing algorithms are used in order to retrieve a reliable assessment of CAR (Figure 3D). Once the data are ready to be processed, different mathematical models exist in order to assess CAR (Figure 3E). The output from these models is then confronted to the state of the neonate, and its prognostic value is assessed. In the following sections we will describe the preprocessing methods and mathematical models more in detail.

Figure 3

Clinical framework to study cerebral flow-pressure autoregulation status using multimodal monitoring. Hereby we propose a setup for determination of flow-pressure CAR in a NICU patient with typical age-appropriate monitoring (A,B). In a multimodal setup (C), invasive ABP and non-invasive vital parameters (SaO2, HR, CO2, Temperature), combined with non-invasive measurement of cerebral oxygenation, are collected continuously in a time-stamped method. If ICP is stable, ABP is a surrogate measurement for CPP. If SaO2 is stable, NIRS derived cerebral oxygenation is a surrogate measurement for CBF. The next step is preprocessing of the data, where the data is down sampled and filtered. Artifact removal and correction for SaO2 is applied (D). Afterwards, several mathematical models can be applied (E). The derived scores provide information about the status of the CAR mechanisms in the patient. Currently, analysis is done offline but real-time bedside information about the CAR status of the patient might be of interest to adapt treatment (written informed parental consent was obtained for publication of this image). BiAR-COH, bivariate autoregressive spectral coherence; BPRSA, bivariate phase rectified signal averaging; CAR, cerebral autoregulation; CBF, cerebral blood flow; CO2, carbon dioxide; COH, coherence; COR, correlation; COx, cerebral oximetry index; CPP, cerebral perfusion pressure; CWT, continuous wavelet transform; FTOE, fractional tissue oxygen extraction; HbD, hemoglobin difference; HbO2, ogygenated hemoglobin; HbT, total hemoglobin; HR, heart rate; HVx, hemoglobin volume index; ICP, intracranial pressure; M/S/DABP, mean/systolic/diastolic arterial blood pressure; NICU, neonatal intensive care unit; NIRS, near-infrared spectroscopy; rScO2, regional cerebral tissue oxygen saturation; SaO2, arterial oxygen saturation; TF, transfer function; TOI, tissue oxygenation index.

Preprocessing

Data preparation

Different methods exist to transfer the requested parameters from the bedside monitor to an off-line system for processing. Even with different manufacturers, similar processing of the often large data files is necessary. To compare studies, detailed description of the different steps used to transfer the data is mandatory. Data is acquired with a given sampling frequency fs. The sampling frequency is extremely variable in the different studies and ranges from 100 to 0.03 Hz (Table 1).

Different filtering techniques are applied to the data. Generally, low-pass filters are used to remove high frequency oscillations. It is particularly common that the systemic data and the NIRS data are acquired at different sampling frequencies. In those cases, data is normally filtered with an anti-aliasing filter and down sampled in order to have a common sampling frequency of all the measurements. Different studies report different final sampling frequency values (22–24).

Artifact removal

When overlooking large data files, sudden, non-physiological changes in baseline or excessive variance can be due to artifacts. The main source of artifacts in the ABP can be due to disconnection of the sensor or blood withdrawal from the arterial catheter. In the case of NIRS, the main source of artifacts is due to the movement or replacement of the sensor, which causes a sudden offset for the measurements. Visual or automated artifact removal is possible and both methods are described in the different studies. However, due to the multifactorial nature of artifacts, automated methods are mostly accompanied by visual inspection. Once artifacts have been detected, they can be corrected by linear interpolation (25) or simply eliminated (20, 26) for further analysis. An example of automated artifact removal is proposed by Scholkmann et al. This method uses a moving window and detects the artifact by detecting sudden changes in the standard deviation of the signals, and corrects it by using spline interpolation of the affected segment (27). Other methods using moving windows are also described (28).

Correction for SaO2

Arterial oxygen saturation (SaO2) has a major influence on NIRS derived cerebral oxygenation, leading to a hypoxic (low oxygen content) cerebral desaturation but not necessarily an ischemic (low blood flow) cerebral desaturation. Different methods have been used in order to correct for the influence of third variables, such as the variability in SaO2. Most prevalent is the exclusion of data with variability in SaO2 larger than 5%. Therefore, conclusions are based on patients during stable SaO2 and extrapolation toward impairment in CAR during desaturation episodes is not possible. Several groups have used cFTOE instead of brain oxygenation to correct for changes in SaO2 (26, 29–31). However, whether this is a valid measurement technique is yet to be defined. De Smet et al. proposed the use of partial coherence (PACOH) in order to correct for variations in SaO2 on the NIRS signals (32). Another technique, partial directed coherence (PDC) has been developed by Baccala et al. to take out the influence of a third signal (33). Although used by Riera et al., this correction was not applied since they excluded the SaO2 signals in their analysis (34). Caicedo et al. proposed the use of oblique sub-space projections (ObSP) (35, 36). ObSP makes use of sub-space system identification that uses input-output observations of the system in order to produce a mathematical model that can explain the measured output. Furthermore, ObSP is able to decouple the linked dynamics between the different underlying subsystems in order to decompose the observed output in terms of the partial contributions of each input variable. This set of signals can be used to define scores for the assessment of the coupling between systemic and brain hemodynamic variables. In essence a correction for SaO2 is provided by eliminating its contribution in the observed NIRS signal, which makes the residual component suited for the assessment of CAR even when changes in SaO2 are present.

Mathematical models

Several mathematical models have been used to assess CAR (37, 38). All of these methodologies try to quantify the relationship between ABP and CBF. These scores are then used to assess the status of the CAR mechanisms in the infants. These models can be divided in 2 main groups, linear and non-linear. The linear models can be further segmented in time-domain, frequency-domain, and non-stationary methods and are discussed below. However, visual inspection of the CAR curve indicates that the CAR mechanism should be nonlinear. Taking this into account several groups have studied CAR using nonlinear models with transcranial Doppler as a surrogate for CBF. However, clinical studies using nonlinear models for CAR assessment in neonates, when NIRS is used as a surrogate for CBF, are not identified during this systematic research.

Time domain

Among the temporal analysis of CAR, correlation (COR) and linear regression are most commonly used. Correlations between the NIRS-derived cerebral oxygenation and ABP measurement are often used to determine flow-pressure CAR. A graphical representation for the COR method can be seen in Figure 4. Different study groups use different correlation methods. Correlations can be carried out in predefined time epochs with or without a predefined threshold (22, 29–31, 39–42). In patients with CHD, a moving correlation coefficient method has been proposed by Brady et al. and named cerebral oximetry index (COx) (43). In preterm infants, COx was used to study the effect of dopamine on CAR (44, 45). A comparable method is used in patients with HIE using HVx (hemoglobin volume index), defining an optimal blood pressure MAPopt from which deviations are correlated with outcome (46–50). Smith et al. used HVx in patients with CHD to study CAR during cardiopulmonary bypass (51). A similar moving correlation coefficient but with the use of Doppler was previously described in adults by Czosnyka et al. (52).

Figure 4

Schematic representation of the correlation. The COR coefficient is a measurement of the linear relationship between two variables. A large and positive COR coefficient indicates impaired CAR, while a small or negative COR coefficient indicates intact CAR. CAR, cerebral autoregulation; CBF, cerebral blood flow; COR, correlation; MABP, mean arterial blood pressure; ρ, correlation coefficient.

Additionally, some groups have explored the use of linear regression for the analysis of CAR. Tyszczuk et al. used COR coefficient and multiple linear regression to study CBF, MABP and other variables. However, no quantitative measurement of CAR was attempted (19). In contrast, Munro et al used linear regression in order to quantify the relationship between MABP and CBF. Based on this analysis they constructed curves for CAR and concluded that a breakpoint exists at 30 mmHg (53).

Frequency domain

Frequency domain analysis explores the relation between 2 signals in specific frequency bands. A major advantage is that it considers the fact that CAR may be composed of responses with different temporal properties, thereby taking into account the effect of time lag between changes in ABP and cerebral NIRS. Determination of the frequencies of interest is based on observed spontaneous oscillations and the cerebrovascular transit time and has an impact on the sample frequency. Different parameter choices like the length of the epochs, the frequency band and others will have an important influence on the calculation as described by Caicedo et al. (23).

Transfer Function (TF) analysis involves the computation of the signal power spectral density (PSD). The PSD of a signal is a representation of how much each particular frequency component contributes to form a given signal. In other terms, each measured signal can be considered as a sum of sinusoids with different frequencies, amplitudes and phases, the PSD represents how much power of each frequency is contained in the measured signal. Different methods exist to compute PSD. The most common is based on the Welch averaged periodogram (54). By measuring the PSD of the ABP and CBF, and their cross-power spectral density (CPSD), the coherence, gain and phase can be computed.

Coherence:

The coherence (COH) function describes the linearity of the relation between two signals in the frequency domain. In the framework of system theory, COH can be used as a measure to indicate if a linear relationship exists between the input and output of the system, which within the context of CAR assessment uses measurements of ABP as input, and NIRS-derived measurements of CBF as output, considering the CAR mechanism as the system to be identified. COH can be seen as the frequency based analog of COR in time. High COH values can be anticipated when CAR is impaired and a threshold for high COH can be defined or calculated according to the length of the windows. COH is used in different studies, with or without correction for SaO2 in specific time frames of 10–15 min or longer (40, 55–60). A graphical representation of the COH method is presented in Figure 5.

Figure 5

Schematic representation of the coherence. The COH is a measure of the linear dependencies between 2 signals. In the figure there are three panels, corresponding to three different conditions for the relation between the input: MABP (gray line), and the output: CBF (black line). Each panel is divided in an upper figure, representing the time course of the signals, a middle figure, representing their PSD, and a lower figure, representing the COH values in the region of interest. Since the input is sinusoidal, we consider the COH value as the value provided in the plot at that specified frequency. In the left panel, the output is contaminated with some noise, however it can be seen that the COH value is large, since the output contains a sinusoid of the same frequency as the input. In the middle panel, the sinusoid has been reduced in amplitude and more noise has been added to the output signal, however, as observed in the figure, the COH value is still large, since the output contains a sinusoid at that specified frequency. In the right panel, no sinusoid has been added to the output and only noise is considered for the computation of the COH. In this case it can be seen that the COH at the specified frequency is low. Within the framework of CAR a low COH value represents intact CAR, while a large COH is associated to impaired CAR. CAR, cerebral autoregulation; CBF, cerebral blood flow; COH, coherence; MABP, mean arterial blood pressure; PSD, power spectral density.

Gain:

TF analysis can be used to study the linear relation between physiological signals. In the framework of system theory and CAR, a TF uses measurements of ABP as input variables and NIRS-derived CBF as output. In this way, when using a TF, a model for the spectral behavior of the CAR mechanism can be found. This process assumes that CAR can be described by a linear, stationary system. The TF of a system can provide two features: the TF gain and the TF phase. TF gain represents the relationship in magnitude between the input and the output, while TF phase represents their temporal relation (time-shifts). A graphical representation of TF gain can be seen in Figure 6. Within the context of CAR, TF gain reflects the change in HbD/rScO2 caused by 1 mmHg change in ABP. The TF gain provides a measure of the magnitude of pressure passivity. A low TF gain would indicate that although CAR was not perfect, at least the magnitude of changes in CBF was small or moderate. This reflects the hypothesis that the plateau of the CAR curve is not flat but increases moderately (61). Similarly, a high TF gain would indicate that even moderate changes in ABP were associated with large changes in CBF. TF gain values have been analyzed over the complete set of measurements (62, 63) or only within coherent epochs (21, 60, 64, 66–68). TF analysis with logarithmic transformation of the gain coefficient to provide the amplitude of the dampening response is described by Vesoulis and measured in decibels (dB) where 0 dB represents no transformation and-10 dB represents a 10-fold reduction in power (69).

Figure 6

Schematic representation of transfer function gain and phase. In the framework of CAR as a system theory, the changes in ABP are considered as input and the changes in NIRS derived CBF as output. The CAR mechanisms are thus considered as a system that will be modeled by its TF. The TF analysis provides two outputs, the gain and the phase. The gain represents the magnitude of the relationship between the variables in the frequency domain, while the phase represents their relative shift in time. Assuming we have as input a pure sinusoidal change in MABP, as indicated in the figure, with this approach we will have a pure sinusoidal change in CBF at the same frequency. In this case the gain at that particular frequency can be interpreted as the ratio between the amplitude of the change in CBF and the change in MABP. A large gain indicates that the change in MABP produces a large change in CBF, while a small gain will indicate that the change in MABP will produce only a small change in CBF. On the other hand, the phase represents the time shift between the two signals. This time shift can be positive or negative. In the context of CAR analysis, a large gain and reduced phase indicate impaired CAR mechanism, while a low gain and large phase indicate intact CAR. CAR, cerebral autoregulation; CBF, cerebral blood flow; (M)ABP, mean arterial blood pressure.

Phase:

In simple terms, the TF phase represents the time lag between the two signals. When using TF phase as a measure for CAR capacity, a high gain with a reduced phase has been associated to impaired autoregulation (60).

Non-stationary methods

Classical methods for the assessment of CAR, such as COR, COH, and TF gain, are often based on the assumption that changes in ABP and cerebral hemodynamics are stationary, i.e., assuming the statistical properties of these signals do not change with time. However, under pathophysiological conditions, but also during stable situations, ABP and CBF can behave non-stationary. Also, time information is averaged out by using TF gain. Thus, other methods are needed to describe the non-stationary aspects of CAR.

Continuous wavelet transform:

The CWT can be used to construct a time-frequency representation of a signal with a very good tradeoff between time and frequency resolution. It makes no assumption about the stationarity of input signals. Therefore, spectral analysis using wavelets provides a mathematical framework for the analysis of nonstationary effects in cerebral hemodynamics, thus overcoming the restrictions intrinsic to earlier methods (25). Tian et al. used rScO2 and described the setting of the neurovascular unit of the group of Chalak, together with EEG measurements (25). Tachtsidis et al. used HbO2 measurements, obtained with a multichannel NIRS device, to analyze the relationship between systemic variables and brain hemodynamics in ECMO patients (24). Using CWT several features or scores can be extracted, for instance compared to ordinary cross-correlations, Wavelet Cross-Correlation (WCC) is also a measure of similarity between two time series, but localized in frequencies, and taking into account the non-stationary nature of physiological measurements.

Bi-variate phase rectified signal averaging:

Bi-variate phase rectified signal averaging (BPRSA) is a method that is used in order to describe the response of a signal to changes in other signals. BPRSA identifies some points of interest in one signal, called the anchor points. These points are traditionally defined as increments or decrements of the signal. Then a segment around each anchor point is extracted from the other signal. All the segments are aligned with the anchor point positions in the middle, and averaged in order to produce the BPRSA response. Since the signal has been averaged over a large amount of segments, the influence of the non-stationarities as well as the noise on the estimation of the BPRSA curve is reduced. The BPRSA curve then represents how one signal reacts to increments or decrements of another variable, which can be of value for the study of the regulatory mechanism. In the context of CAR, Caicedo et al. used increments in MABP and HR as anchor points to study their effect on rScO2 and by studying the BPRSA curve they hypothesized that it was possible to assess CAR. They found differences in the obtained BPRSA curves between control subjects and subjects with an IVH (71).

Bivariate autoregressive coherence:

COH analysis, as presented before, is a very useful tool to detect linear dependencies among signals. However, in the case of third signals influencing both measurements, or in case that the interest lies in identifying the direction of influence, COH analysis is not the adequate tool. For this reason Riera et al. have used the bivariate autoregressive coherence (BiAR-COH). In short, BiAR-COH is able to produce a more reliable estimation of the linear dependencies among the signals, and indicate the directionality of the coupling. In order to do this, BiAR-COH first models the ABP as well as the TOI as follows: . When computing the COH assuming that ABP leads the changes in TOI, COHABP → TOI, it makes more sense to use the component of TOI that is affected by ABP, TOIABP−, instead of using the original ABP and TOI measurements. In this way it can be evaluated if changes in TOI can be predicted linearly by changes in ABP. BiAR-COH is computed as a ratio between the cross-power spectrum of , and it is normalized using the sum of this cross-power spectrum and the power spectrum of . Riera et al. found that this method has a good prognostic value to identify neonates at risk of brain hypoperfusion and adverse outcomes (72). Furthermore, they proposed the PDC as a valuable complementary analysis to BiAR-COH (34).

Corrections for variability in ABP:

One of the main drawbacks of different models for the assessment of CAR is that they require large changes in ABP and CBF in order to produce reliable scores. In long recordings it is possible to encounter these large changes. However, changes in ABP are not regular and thus non-stationary. Hahn et al. has indicated that weighting measurements with large in favor of those with small variations in ABP, increases the precision in the assessment of CAR (66, 67). Also other groups have taken the ABP variability into account for the assessment of CAR (73).

Definition of impaired autoregulation

CAR is defined as a stable CBF during changes in CPP.

The CAR curve describes a plateau with a stable CBF between a lowest (minimum autoregulation pressure) and a highest (maximum autoregulation pressure) ABP. When the ABP decreases further, an almost linear relation between ABP and CBF is described possibly leading to ischemic lesions, while a further increase in ABP at the end of the plateau will give a linear increase in CBF leading to an overflow with possible IVH as a consequence. The plateau is not flat but slightly increasing explaining the low but significant gain during periods of autoregulation (61). The CAR curve can change in form by changes in CO2, PO2, acidosis, NO and other parameters causing vasodilation or vasoconstriction.

A loss of CAR means that the patient is -at that specific moment- outside the plateau of the curve. In contrast with earlier visions, recent observations suggest that this can change in individual patients according to changing parameters like ABP, CO2, sepsis, … Therefore, increasingly, studies describe the percentage of time of impaired CAR. However, it is important to notice that a loss of CAR does not necessarily means that the patient will develop complications, but is regarded as a higher risk for complications. Consequently, the clinician will aim for a stable CAR during treatment, within possible limits.

The definition of impaired CAR requires the definition of parameters that characterize where the individuals are located within the CAR curve or a specific cut-off point. Depending on the mathematical model used to study CAR, definitions can change. CAR can be seen as an on-off phenomenon. Different authors assess this phenomenon by identifying when the changes in ABP are linearly related to the changes in CBF. This is done by indicating when there is maximum COH, or COH above a certain threshold (20, 21, 59). However, others denote impaired CAR as a curve with a large slope in the autoregulation plateau and changes in the infliction points (53). Other groups use the change in CBF divided by the change in ABP or the autoregulation index (ARI) of static autoregulation (15). A wide range of definitions exist and thus, different ways of describing the autoregulation curve and its borders are used.

As mentioned before, some methods make use of thresholds in order to identify the instants of time when CAR is impaired. Different thresholds have been described using mathematical models, adult studies, animal studies, Monte Carlo simulations or visual inspection of data. Apart from cut-off values for COH and COR, also expressed as critical value score (CVS) (20, 40, 55, 57, 59), the threshold can also be expressed as amount of pressure passivity per total study time or per number of epochs. CPRT (critical percentage of recording time) (40, 57) and pressure passivity index (PPI) of COR,COH or TF gain is used by different groups (41, 57, 59, 68). The identifications of the lower infliction point is estimated by defining a lower limit of pressure autoregulation (LLA) by Brady et al. (43). The ABP associated with the most negative COR coefficient (less impaired CAR) is defined as MAPopt in HIE studies (46–50).

Impaired cerebral autoregulation and study outcome

Most studies in CAR try to relate the scores used for the assessment of CAR with clinical outcome. We will concisely present the different study outcomes for the different patient groups earlier described.

Preterms during the transitional period

This is the main research area for CAR in neonates and 25 articles were identified describing both short and long term outcome parameters. Starting in 2000, Tsuji described impaired CAR in 53% of ventilated preterms with a strong relationship between this impairment and severe IVH or periventricular leukomalacia (PVL) (20). The association between IVH and impaired CAR was confirmed by different groups (34, 42, 64, 69, 71–73), but not by others (30, 44, 63, 66).

A good correlation was described between impaired CAR and sickest infants, hypotension, and a higher CRIB-score. Different studies described the relation between low ABP and impaired CAR in preterms. Munro et al. described a cut-off value at 30 mmHg (53). Soul et al found a good correlation with hypotension, defined as less than the 10th percentile for PMA and postnatal age (59). They did not find a good correlation with blood pressure variability. Also, Hahn et al. described an association between impaired CAR and hypotension (defined as mean ABP in mmHg minus GA in weeks), but not with blood pressure variability (66). Gilmore et al. and Wong et al. described a correlation with lower ABP although no cut-off value was provided (21, 44).

Regarding long term outcome, impaired CAR was associated with a higher mortality and worse neurodevelopmental outcome by different groups (21, 34, 44, 59, 60, 63, 72). All studies were retrospective and no controlled prospective studies were performed.

Neonates receiving medication/treatment

Studies have explored whether treatment with volume/inotropes influences CAR. Treatment with dopamine is associated with impaired autoregulation in 4 studies (41, 45, 53, 74).

Neonatal effects on CAR of maternal drug treatment during the last days of pregnancy have been investigated (58, 62). Also, effects of certain postnatal drug treatments on CAR are described: surfactant administration in RDS (75), patent ductus arteriosus (PDA) treatment with indomethacin and/or surgery (41).

Extracorporeal membrane oxygenation (ECMO) treatment also disrupts CAR in neonates. The mechanisms are multifactorial, including heparinisation, as well as hemodynamic instability. In this context, it has been found that during low ECMO flow values, CAR is impaired in neonates undergoing ECMO weaning (24). In this study was also shown that the right hemisphere is more susceptible to disrupted CAR.

Patients with CHD

A single study has identified impaired CAR in 24 neonates with CHD in the preoperative phase. An association with lack of sedative medication, low hemoglobin, low mean ABP and greater FTOE was found (76). The group of Brady et al. has described an association with both hypotension and impairment of CAR during CPB. A cut-off value of 42±7 mmHg in pediatric patients undergoing cardiac surgery was described (43). However, hypothermia was identified as a confounding factor since ABP, temperature and CAR were collinear in a cohort studied later (51). In the early postoperative phase after CPB, higher end-tidal CO2, higher mean ABP variability but not lower ABP increased the odds of impaired CAR (55).

Patients with HIE treated with therapeutic hypothermia

A selected patient group to study CAR are infants with HIE after a perinatal ischemic event. Therapeutic hypothermia decreases the metabolic rate of the brain to preserve the remaining brain tissue as much as possible (78). By defining the optimal MAP (MAPopt) as the MAP where CAR is most robust, impaired CAR is identified in % of time, maximum deviation from and area under the curve below MAPopt. Retrospective studies where ABP measurements were set against this MAPopt during hypothermia, rewarming and after the therapy showed an association between the impaired CAR and cardiopulmonary injury, subjective and objective injury on magnetic resonance imaging (MRI) and long term outcome (46–50). Other studies using different definitions have correlated impaired CAR with lesions on MRI and outcome (25, 68).

Discussion

In this manuscript we provide an extensive overview of the measurement of CAR in neonates using NIRS.

Clinical framework of CAR measurement

Different studies have shown that it is possible to measure CAR in neonates. No clear evidence for the best clinical framework has been provided. This may be caused due to the lack of a gold standard for its assessment. The lack of standardization on preprocessing methods and mathematical models used, length of window of analysis and thresholds of impaired CAR hinders the reproducibility of studies that have related CAR assessment with clinical outcome.

One of the first questions arising concerns the sample frequency. According to the Nyquist theorem, the sampling frequency must be at least twice the highest intrinsic frequency in the signal. To measure CAR, spontaneous cardiovascular oscillations are used. Some studies use the intrinsic variations of ABP within the low to very low frequency (frequency ≥0.0067 Hz, which corresponds to oscillations of 150 s) which represent the composite influence of autonomic, myogenic and cellular control mechanisms (79). Others use oscillations in the higher frequency range. The oscillations with a frequency component around 0.1 Hz are related to the Mayer waves and result from an oscillation of the sympathetic vasomotor tone. Several theories have been proposed to explain the constancy of the Mayer waves frequency (80). To use these waves to measure CAR, sampling frequencies should be at least 0.2 Hz (fs ≥ 0.2 Hz, 1 value every 5 s). By using a lower signal sample frequency, information about the autoregulatory capacity of the brain around these higher frequencies will be lost. Several groups avoid measuring CAR using this 0.1 Hz frequency since the sensitivity of NIRS to brain activation is diminished by these physiological fluctuations that arises from scalp, skull and brain. These systemic changes may result in false positives by mimicking the brain hemodynamic response or result in false negatives by attenuating it (81). However, the question remains at which frequencies CAR can be studied best and which influences are measured.

Soul et al. and Bassan et al. were the first to measure bedside during days but often this involved a fellow being constantly present to make corrections in case of movement artifacts (55, 59). Therefore, the signal preprocessing, specifically artifact detection and removal, still needed further development in order to provide a reliable solution for the bedside monitoring of CAR. It has been shown that when assessing CAR, small artifacts can be easily truncated or interpolated (82). However, if large artifacts are present and large interpolations are made, the question arises if the results represent the real mechanisms under study.

In this context some efforts have been made in order to evaluate whether is possible to use some technologies that are more robust for CAR assessment. Caicedo et al, have shown that TOI and rScO2, can be used instead of HbD and HbT for the study of CAR. TOI and rScO2, derived from spatially resolved spectroscopy, are less prone to artifact movements than HbD and HbT (22).

Comparisons between mathematical models

Few mathematical models are validated in animal studies. The COH and TF gain method to measure CAR is validated in piglets (83). A good correlation with TF gain was found when the COH was larger than 0.47. However, TF gain or COH alone were less optimal parameters, suggesting both parameters need to be used together to obtain a good assessment of CAR. The COx has been validated in piglets with hypotension (84) and correlates with transcranial Doppler-derived measurements of flow-pressure CAR in adult patients (43). Eriksen et al. compared COx (time domain) vs. COH and TF gain (frequency domain). They reported that time-domain analysis appeared more robust compared with COH function analysis. TF gain also increases when ABP and OI are in counter phase (67). This can be solved by including TF phase values, which provide an indication of the temporal shift between the signals in the frequencies of interest. One of the main drawbacks of the TF gain model is that it requires the presence of large variations, in a large range of frequencies, in MABP and CBF, in order to obtain reliable estimations of the spectral content of the signals. On the other hand, in the COx model, as with other COR models, the phase shift between the signals is not taken into account. These analyses assume that changes in ABP are immediately reflected in changes in CBF. Also, both COH and TF gain and COx models assume linearity, while this is not always true in biological signals. More studies comparing COH and TF gain vs. COx will elucidate the value of each model.

Another important issue concerns the precision of different details concerning the mathematical models. A comparison between selected details was published by Caicedo et al. The epoch length (segment on which the scores are calculated), overlapping percentage (different parameters in the COR and the frequency based methods), and sub-window length (segment used in the Welch method) all influenced importantly the final scores (23). This has been noticed by the Cerebral Autoregulation Network (CARNET, http://www.car-net.org/), which reported that TF gain and phase, which was used for CAR assessment, using ABP as input and transcranial Doppler measurements as output, change significantly between different centers. This was mainly due to the differences in preprocessing, selection of parameters for the estimation of PSD, among others (85). This phenomenon is also noted in the included studies causing difficulties in comparing study methods, definition of impaired CAR and outcome.

To overcome the non-stationary nature of biological signals, Chalak et al. focused on what they call the neurovascular unit approach, which involves the continuous monitoring of different modalities to identify potential mechanisms of dysfunctional CBF regulation by using the CWT. This model, compared to more traditional models, has the advantage to allow identification of the spectral content of the signal, together with the changes in time. Since CAR is a dynamic process, CWT seems to be more suited for its analysis. This strategy is proposed as a neuromonitoring tool in HIE to determine dynamic CAR (25, 86).

Study outcome in CAR research

From a physiological point of view, impaired CAR is not by definition associated with pathology. Theoretically, infants with less adapted mechanisms to overcome the change in ABP are at risk for complications and can possible be identified by using CAR monitoring. If impaired CAR is identified, not all patients will develop complications.

Both short and long term outcome parameters are used in the studies described in this review. In general, the incidence of IVH is the outcome variable most frequently correlated with impaired CAR in preterm infants. Multiple animal and clinical studies have demonstrated that disturbed CBF is associated with IVH and PVL (1). This is reflected in studies where impaired CAR is associated with IVH using different mathematical methods (20, 42, 64, 69, 71–73). Alderliesten et al described a detailed time frame with impaired CAR and hyperperfusion before the occurrence of IVH. Other possible pathophysiological mechanisms are low SVC flow (34, 72) or low cardiac output (9, 19, 87). However, this association of IVH and impaired CAR is not confirmed in other studies (41, 44, 59). It is unknown whether the used method lacks precision or that in the subgroup studied, the impaired CAR did not necessarily provoke brain injury or other complications.

Correlation between impaired CAR and long term outcome is difficult due to other confounders (i.e., PMA, bronchopulmonary dysplasia, infection, medication and parental factors). This might be the reason why correlation with long term outcome is weak or absent.

CAR assessment during pharmacodynamic research is of special interest, since many drugs used in neonatal care are prescribed off-label. Identification of the drugs hemodynamic impact, both central and peripheral, is of utmost importance. Impairment of CAR after dopamine therapy is described in 4 previously mentioned studies. Also, in an extensive meta-analysis, increases in CBF following dopamine infusion were greater in hypotensive than normotensive preterm infants suggesting that dopamine does not exert a selective vasodilatory effect in the cerebral circulation and that the ABP of infants treated for hypotension was below the lower elbow of the CAR curve (88). However, impaired CAR might have developed by other confounders in this small patient group. A placebo-controlled randomized trial with dopamine in extremely premature infants is currently recruiting, to elucidate the effect of dopamine on long term outcome (89). A large patient group with measurements of invasive ABP and cerebral oxygenation will be available to study CAR during hypotension and treatment with dopamine or placebo.

Limitations and future research

Due to the very heterogeneous study methods and outcome parameters, a meta-analysis was not performed. Due to small patient groups, corrections for other known variables defining long term outcome are not performed in most of the studies. At this stage, prospective interventional studies are lacking in all research areas.

After the identification of the patient at risk, continuous bedside CAR measurement can improve knowledge about cerebral pharmacodynamics of frequently used medication (sedatives-inotropes) and cerebral effects of specific therapies. However, most CAR studies are done offline, and at present, offline analysis is the standard. Online analysis to guide care, with bedside treatment adaptation if CAR is failing, may be possible in the future. In this way, physiology will be displayed on the multimodal monitor, hand in hand with the clinical examination.

A limitation for the continuous monitoring of CAR is the presence of an arterial line to monitor ABP invasively. A promising technique is to use heart rate instead of ABP to describe CAR. Describing the functional activation of the brain, heart rate is used as a surrogate marker for neurogenic activity (35, 60, 71, 90, 91). Da Costa et al. used moving correlation between heart rate and TOI to determine optimal ABP (92). Cardiac output is another possible variable that can be used as a surrogate for APB (9, 93).

In the future, multimodal monitoring databases of NICU-patients of multiple centers, together with identical outcome parameters are needed to compare different mathematical models and make progress in this field. To achieve this, easy, cheap, user-friendly NIRS-equipment with validated neonatal sensors is necessary.

Conclusion

NIRS derived CAR measurement is an important research tool to improve knowledge about central hemodynamic fluctuations during the transitional period, cerebral pharmacodynamics of frequently used medication (sedatives-inotropes) and cerebral effects of specific therapies in neonatology. Uniformity regarding measurement techniques and mathematical models is needed. Multimodal monitoring databases of neonatal intensive care patients of multiple centers, together with identical outcome parameters are needed to compare different techniques and make progress in this field. Real-time bedside monitoring of CAR, together with conventional monitoring, seems a promising technique to improve individual patient care.

Author contributions

LT screened the article titles and abstracts to determine whether they met the inclusion criteria. Then, LT reviewed full text articles to assess for eligibility. Any articles presenting doubts or inconsistencies were fully reviewed by GN and AC until a decision was reached on their inclusion. The first draft was made by LT, GN, and AC and all other authors contributed to the final article.

Conflict of interest statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.