Non-invasive Optical Measurement of Cerebral Metabolism and Hemodynamics in Infants

Summary

We combined frequency-domain near-infrared spectroscopy measures of cerebral hemoglobin oxygenation with diffuse correlation spectroscopy measures of cerebral blood flow index to estimate an index of oxygen metabolism. We tested the utility of this measure as a bedside screening tool to evaluate the health and development of the newborn brain.

Abstract

Perinatal brain injury remains a significant cause of infant mortality and morbidity, but there is not yet an effective bedside tool that can accurately screen for brain injury, monitor injury evolution, or assess response to therapy. The energy used by neurons is derived largely from tissue oxidative metabolism, and neural hyperactivity and cell death are reflected by corresponding changes in cerebral oxygen metabolism (CMRO₂). Thus, measures of CMRO₂ are reflective of neuronal viability and provide critical diagnostic information, making CMRO₂ an ideal target for bedside measurement of brain health.

Brain-imaging techniques such as positron emission tomography (PET) and single-photon emission computed tomography (SPECT) yield measures of cerebral glucose and oxygen metabolism, but these techniques require the administration of radionucleotides, so they are used in only the most acute cases.

Continuous-wave near-infrared spectroscopy (CWNIRS) provides non-invasive and non-ionizing radiation measures of hemoglobin oxygen saturation (SO₂) as a surrogate for cerebral oxygen consumption. However, SO₂ is less than ideal as a surrogate for cerebral oxygen metabolism as it is influenced by both oxygen delivery and consumption. Furthermore, measurements of SO₂ are not sensitive enough to detect brain injury hours after the insult ^1,2, because oxygen consumption and delivery reach equilibrium after acute transients ³. We investigated the possibility of using more sophisticated NIRS optical methods to quantify cerebral oxygen metabolism at the bedside in healthy and brain-injured newborns. More specifically, we combined the frequency-domain NIRS (FDNIRS) measure of SO₂ with the diffuse correlation spectroscopy (DCS) measure of blood flow index (CBF_i) to yield an index of CMRO₂ (CMRO_2i) ^4,5.

With the combined FDNIRS/DCS system we are able to quantify cerebral metabolism and hemodynamics. This represents an improvement over CWNIRS for detecting brain health, brain development, and response to therapy in neonates. Moreover, this method adheres to all neonatal intensive care unit (NICU) policies on infection control and institutional policies on laser safety. Future work will seek to integrate the two instruments to reduce acquisition time at the bedside and to implement real-time feedback on data quality to reduce the rate of data rejection.

Introduction

The FDNIRS device is a customized frequency-domain system from ISS Inc. with two identical sets of 8 laser diodes emitting at eight wavelengths ranging from 660 to 830 nm, and two photomultiplier tube (PMT) detectors. Sources and detectors are modulated at 110 MHz and 110 MHz plus 5 kHz, respectively, to achieve heterodyne detection ⁶. Each laser diode is turned on for 10 msec in sequence, for a 160 msec total acquisition time per cycle. Sources and detectors are coupled to fiber optics and arranged in a row in an optical probe. The arrangement of fibers on the probe is such that it produces four different source-detector separations. By measuring transmitted light (amplitude attenuation and phase shift) at multiple distances, we can quantify the absorption (μa) and scattering (μs') coefficients of the tissue under observation. From the absorption coefficients at multiple wavelengths, we then estimate the absolute values of oxygenated (HbO) and deoxygenated (HbR) hemoglobin concentrations ⁷, cerebral blood volume (CBV) and hemoglobin oxygen saturation (SO₂).

The DCS device is a home-built system similar to the one developed by Drs. Arjun Yodh and Turgut Durduran at the University of Pennsylvania ^8,9. The DCS system that consists of a solid-state, long coherence laser at 785 nm, four photon-counting avalanche photodiode (APD) detectors (EG&G Perkin Elmer SPCM-AQRH) featuring low dark counts (<50 counts/sec) and a high quantum yield (>40% at 785 nm), and a four channel, 256-bin multi-tau correlator, with 200 nsec resolution. With the DCS we measure microvascular blood flow in cerebral cortex by quantifying the temporal intensity fluctuations of multiply scattered light that arises from Doppler shifts produced by moving red blood cells. The technique, similar to laser Doppler blood flowmetry (i.e. they are Fourier Transform analogs), measures an autocorrelation function of the intensity fluctuations of each detector channel computed by a digital correlator over a delay time range of 200 nsec - 0.5 sec. The correlator computes the temporal intensity auto-correlation of the light re-emerging from tissue. We then fit the diffusion correlation equation to the measured autocorrelation function, acquired sequentially, about once per second, to obtain the blood flow index (CBF_i) ^10,11. DCS measures of blood flow changes have been extensively validated ^12,13. By combining the FDNIRS measures of SO₂ with the DCS measures of CBF_i, we achieve an estimate of cerebral oxygen metabolism (CMRO_2i).

Protocol

1. Preparation for Bedside Measures

1. The FDNIRS and the DCS systems are compact and easy to move on a small cart to the infant's bedside in the hospital (Figure 1).
2. After moving the cart with the devices to the bedside, turn on the systems and connect the optical probe to the FDNIRS and DCS devices. Ensure that two experimenters are present for every measurement: one to manage the instruments and computer, and one to hold the probe.
3. Choose the appropriate probe according to the infant's postmenstrual age (PMA). The optical probe with FDNIRS source-detector separations of 1, 1.5, 2 and 2.5 cm is used for infants <37 wks PMA and the probe with FDNIRS separations 1.5, 2, 2.5 and 3 cm is used for older infants (Figure 2-A). The choice of shorter source-detector separations is dictated by preterm infants' small size and larger head curvature. When using a larger probe with a preterm infant, the relatively smaller size of the baby's head and its significant curvature together impede effective contact between the infant's head and all sources and detectors. For this reason, the probe with FDNIRS source-detector separations of 1, 1.5, 2 and 2.5 cm is fitting for use with preterm infants. Our research has verified that the chosen source-detector separations are sufficient to measure optical properties of the cerebral cortex of both preterm and term ¹⁴. DCS source and detector fibers are arranged in a row parallel to the FDNIRS fibers with source-detector distances of 1.5 (one detector) and 2 cm (three detectors) in both premature and term infants probes.
4. Sanitize the optical probes with a Sani-cloth disinfecting wipe and insert the probe and fibers into a single-use polypropylene plastic sleeve.

2. FDNIRS Gain Settings and Calibration

1. Open the FDNIRS Graphical User Interface (GUI) and select the program settings file corresponding to the probe and calibration block being used.
2. To adjust detector gains, gently place the probe on an area of the subject's head without hair (preferably the left side of the forehead) and maintain it in the same position without applying any pressure. Turn on sources and detectors and adjust PMT voltage until the amplitude of any of the source lasers reaches 20,000 counts. 32,000 counts is the maximum digitization of the analog to digital acquisition card, and gains need to be set below that threshold to avoid saturation during data acquisition. The gains should be set in the frontal area because this region generally has the lowest absorption of laser light and is therefore most prone to saturation.
3. Turn off the sources and detectors and return the probe to the calibration block. The lasers need to be turned off when moving the probe for eye safety; the detectors need to be turned off because PMTs are very sensitive and exposure to any bright light increases background noise and may permanently damage them.
4. With the probe back on the calibration block, use the neutral density (ND) filter if any source-detector pair saturates. Different ND filters may be selected due to optimizing gains in infants with different skin tones Hold the probe still for 16 sec while running the calibration procedure. Since we do not physically move one source to different distances from a single detector to achieve a multi-distance scheme, but instead use four combinations of two independent sources and two independent detectors, we need to calibrate for the different power of the two sources and the different gains of the two detectors. By measuring a calibration block of known optical properties, we estimate the amplitude and phase correction factors needed to recover the absorption and scattering coefficients of the calibration block.
5. After calibration, acquire 16 more sec of data on the block and visually assess the adequacy of the calibration with an in-house MATLAB GUI. The measured μa and μs' should match the actual coefficients of the calibration block at all wavelengths. Recalibrate if the fit is poor.
6. If detector gains need to be changed, or source and detector fibers need to be disconnected during measurements, repeat the calibration procedure of the FDNIRS device.
7. At the end of the measurement session, acquire another 16 sec of data on the calibration block to verify whether the calibration was maintained during measurements on the subject. If the calibration has not been maintained, take a second calibration at the end of the measurement and apply to the acquired data.

3. DCS Settings

1. Open the in-house DCS data acquisition GUI and load the settings file corresponding to the optical probe being used.
2. Before starting measurements, verify that the laser power of the DCS source is appropriate for skin exposure by measuring the laser power of the DCS source with a power meter and checking the spot size with an IR viewing card (the laser emits at 785 nm, which is not visible). The DCS laser power is ~60 mW and coupled to a relatively small diameter fiber (400-600 μm). To meet ANSI standards for skin exposure, the light at the probe must be attenuated and diffused across a large area. This is achieved by covering the tip of the fiber with a 3 mm diameter white Teflon sheet (Figure 2-A). The Teflon is highly scattering and widely diffuses the laser beam. At the bedside, ensure that the laser power at the probe is less than 25 mW and the spot size is larger than 3 mm in diameter. As for the FDNIRS, always turn off sources and detectors when moving the optical probe.
3. DCS detection is photon-counting and there is no APD gain adjustment as is required for the FDNIRS device. A flag in the acquisition software indicates if too much light is detected, in which case light coupling to either source or detector fibers needs to be reduced by turning the fiber connectors. Adequate light detection is on the range of 200,000-4,000,000 detected photons (corresponding to -26~0 dB on the computer display). Avoid excessive room light to reduce background noise.
4. The DCS does not require calibration to measure CBF_i. Blood flow is proportional to the time it takes to lose correlation. A solid block does not suffice to check signal quality because there are no moving scattering particles to cause decay. An experimenter's arm instead shows decay -- the faster the blood flow, the steeper the decay.

4. Data Acquisition

1. While FDNIRS and DCS measurements can be done quickly in sequence, first measure all locations with one device and then repeat the same progression with the other device, using independent acquisition software corresponding to each.
2. Measure seven locations covering frontal, temporal and parietal areas, according to a 10-20 system (Fp1, FpZ, FP2, C3, C4, P3, P4), in sequence (Figure 2-B). Part the hair along the source-detector line and place the probe on that area of the head.
3. Turn on FDNIRS lasers and detectors and check the signal quality: amplitude counts should be between 2,000 and 20,000 and phase shifts SNR <2 degrees. If outside of these ranges, reposition the probe, ensuring hair is parted and all sources and detectors are in contact with the skin.
4. Acquire data for 16 sec. Repeat measurements up to three times in each location (Figure 2-C), parting the hair and repositioning the probe in a slightly different spot for each acquisition. This is done to minimize effect of local inhomogeneities such as hair and superficial large vessels and to provide values representative of a region, rather than a single spot.
5. Turn on the DCS laser and detectors and acquire data for 10 sec. Reposition the probe and repeat the acquisitions (as with the FDNIRS measures).
6. Turn off all lasers when moving the probe between locations. Data collection in all seven locations is not always possible. Discontinue measurements if the subject manifests any sign of distress or motion. Retry the acquisition if possible. EEG electrodes or respiratory equipment may also preclude measurements in some locations.

5. Measure of Systemic Parameters

1. For the calculation of CMRO_2i, two systemic parameters, arterial oxygenation (SaO₂) and hemoglobin in the blood (HGB), must be acquired. HGB is also needed to calculate CBV. While conventional pulse oximetry provides measures of SaO₂, HGB is conventionally measured with a blood test. A new pulse oximeter, developed by Masimo Corporation, is able to measure HGB non-invasively by using multiple wavelengths. The device is FDA-approved for infants >3 kg, and allows for a quick bedside measure of both SaO₂ and HGB.
2. Record SaO₂ and HGB using a Masimo pulse oximeter (Pronto spot check pulse co-oximeter). For these measurements, attach an adhesive single-use sensor to the big toe of the baby's foot. HGB will be displayed on the monitor within a few seconds.
3. When it is not possible to use the Masimo pulse co-oximeter, measure SaO₂ with other FDA-approved pulse oximeters. HGB can be either recovered from patients' clinical charts or estimated using normative values.

6. Data Analysis

1. Open an in-house post-processing data analysis GUI using MATLAB. This software not only estimates all hemodynamic parameters, but also uses the redundancy of data to automatically assess measurement quality and constrain results.
2. Automatic objective criteria for quality control consist of discarding data for FDNIRS if: R2 < 0.98 for the model fit of the experimental data, p-value > 0.02 for the Pearson product moment correlation coefficient between the eight measured absorption coefficients and the hemoglobin fit (Figure 3-A), p-value > 0.02 for the linear fit of the reduced scattering coefficients versus wavelength (Figure 3-B) ¹⁵. If more than 33 percent of the data merits discarding, the whole set is discarded. For the DCS, data is discarded if: the tail of the fitting curve differs from 1 by more than 0.02, the cumulative variation among the 3 first points of the curve is more than 0.1, or the average value of the 3 first points is more than 1.6 (Figure 4). If more than 50 percent of curves are discarded, or the fit values have a coefficient of variation > 15 percent, the whole set is discarded¹⁵.
3. Calculate absolute HbO and HbR concentrations by fitting the absorption coefficients at all wavelengths, using literature values for Hb extinction coefficients ¹⁶ and a 75 percent concentration of water in tissue ¹⁷. Derive total hemoglobin concentration (HbT = HbO + HbR) and SO₂ (HbO / HbT) from HbO and HbR concentrations.
4. Estimate cerebral blood volume using the equation described in Ijichi et al ¹⁸. CBV= (HbT × MW_Hb) / (HGB × D_bt), where MW_Hb = 64,500 [g/Mol] is the molecular weight of Hb, and D_bt = 1.05 g/ml is the brain tissue density.
5. Calculate CBF_i by fitting the measured temporal autocorrelation functions to the diffusion correlation equation. The detailed theoretical framework to calculate CBF_i is in Boas et al. and Boas and Yodh ^10,11. In the equations, use individual absorption coefficients measured from FDNIRS and an average of the scattering coefficients across the whole population.
6. Calculate the index of cerebral oxygen consumption by using the FDNIRS measure of SO₂ and the DCS measure of CBFi with the following equation: CMRO_2i = (HGB × CBFi × (SaO₂ - SO₂)) / (4 × MW_Hb × β)¹⁵, where the factor 4 reflects the four O₂ molecules bound to each hemoglobin and β is the percent contribution of the venous compartment to the hemoglobin oxygenation measurement ¹⁹.

Results

In the past five years we have demonstrated the feasibility and clinical utility of the proposed method. In particular, we have shown CMRO₂ to be more representative of brain health and development than SO₂.

In a cross-sectional study on more than 50 healthy infants, we found that while CBV is more than double during the first year of life, SO₂ remains constant ⁴ (Figure 5). In a study on 70 healthy newborns we also found that SO₂ is constant across brain regions while CMRO_2i, CBV and CBF are higher in temporal and parietal regions than in the frontal region (Figure 6)²⁰, which is consistent with PET glucose uptake findings ²¹. In both of our studies, the constant SO₂, within a 60-70 percent range indicates that oxygen delivery closely matches local consumption, while CBV, CBF and CMRO₂ are more tightly coupled with neural development.

To verify that CMRO_2i is a better screening tool than SO₂ in detecting neonatal brain injury, we measured brain injured infants during the acute phase ⁵, and (in a few infants) during the chronic phase several months after injury. Results in Figure 7 show how SO₂ is not significantly altered by brain injury in both early (1-15 days after insult) and chronic (months after injury) stages, while CMRO_2i is significantly different than normal during both the acute and chronic stages. Specifically, CMRO_2i is elevated during the acute phase because of seizure activity after brain injury, and lower than normal during the chronic phase due to neuronal loss.

Infants with hypoxic ischemic injuries are currently treated with therapeutic hypothermia (TH) to lower brain metabolism and reduce damage after the hypoxic insult. Therapeutic hypothermia is maintained for three days and we have been able to monitor 11 infants during treatment (Figure 8). We found that CMRO_2i significantly decreases to levels below normal during TH, and this decrease seems to be related to response to therapy and developmental outcome. These preliminary results suggest that the FDNIRS-DCS method may be able to guide and optimize hypothermia therapy.

Figure 1. Picture of the cart with the FDNIRS and DCS devices. The two instruments are compact enough to fit on a small cart that can be moved to the infant's bedside in the NICU.

Figure 2. (A) Optical probe configuration. (B) The measurement location scheme. (C) A photo of a typical FDNIRS-DCS measurement on an infant.

Figure 3. Representative examples of good and bad fit of measured (A) absorption coefficients and the hemoglobin fit (B) scattering coefficients and the linear fit. P-value > 0.02 refers to a bad fit. Click here to view larger figure.

Figure 4. A representative example of good and bad fit of an autocorrelation function of the intensity fluctuations computed by a correlator over a delay time range of 200 nsec - 0.5 sec. In the bad fit figure the tail of the fitting curve differs from 1 by more than 0.02 and the variation of the 3 first points is more than 0.1. Click here to view larger figure.

Figure 5. Changes in CBV and SO₂ across frontal, temporal and parietal cortical regions in infants from birth to one year of age.

Figure 6. CBF, SO₂, CBV and CMRO_2i of the frontal, temporal and parietal regions in 70 healthy newborns.

Figure 7. Examples of abnormal oxygen consumption and normal SO₂ after brain injury in infants. Brain injury is marked by changes in CMRO₂ with respect to normal while SO₂ is not significantly different from normal. Please note that in these two figures, CMRO₂ was calculated using the Grubb relationship, because the DCS measure was not available at the time of those measurements.

Figure 8. rCMRO₂ of 11 infants during therapeutic hypothermia vs. age-matched healthy controls. Oxygen metabolism is strongly reduced in all infants with hypothermia therapy.

Discussion

We demonstrated a quantitative measurement of cerebral hemodynamic and metabolism with FDNIRS and DCS in the neonatal population. The probe configuration is optimized for measuring neonatal cerebral cortex ¹⁴. Blood flow changes measured by DCS have been extensively validated against other techniques in animal and human studies ^22,23. By using a direct DCS measure of blood flow, we are able to reduce the variance in the calculation of CMRO_2i ²⁴. The variance from repeated measures was also smaller than the changes between brain regions and with age ²⁰.

From our previous results, CBFi and CMRO_2i showed significant changes with PMA in healthy preterm neonates. The measure of CMRO_2i is better able to detect brain damage than the measure of SO₂. This suggests that combined measures of vascular and metabolic parameters serve as more robust biomarkers of neonatal brain health and development than oxygen saturation alone. Technical improvements will focus on the integration of two instruments to reduce acquisition time 35-40% per session and the implementation of real-time feedback on data quality to reduce the frequency of discarded measures. In the near future, this system can be delivered to clinical end-users as a novel bedside monitor of altered cerebral oxygen metabolism. By measuring trajectories of CMRO₂ over time may also increase clinical significance and predict outcomes. This tool could ultimately make a significant contribution towards improved management of neonatal brain injury.

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
Imagent	ISS		FDNIRS
DCS laser fibers	Thorlabs	FT400	DCS component
DCS detector fiber	Thorlabs	780HP	DCS component
DCS laser	CrystaLaser	DL785-070-S	DCS component
DCS detector	Pacer International	SPCM-AQRH-14-FC	DCS component
DCS Correlator	Correlator.com	Flex05-8ch	DCS component
Pronto co-oximeter	Masimo		HGB and SaO2 monitor
NOVA	OPHIR	7Z01500	Laser power meter
Sensor card	Newport	F-IRC1-S	IR viewer
Neutral Density filter	Kodak	NT54-453