Application of a time-resolved optical brain imager for monitoring cerebral oxygenation during carotid surgery

Investigating cortical hypoxia in multiple sclerosis via time-domain near-infrared spectroscopy

OBJECTIVES

Hypoperfusion and tissue hypoxia have been implicated as contributory mechanisms in the neuropathology of multiple sclerosis (MS). Our objective has been to study cortical oxygenation in vivo in patients with MS and age-matched controls.

METHODS

A custom, multiwavelength time-domain near-infrared spectroscopy system was developed for assessing tissue hypoxia from the prefrontal cortex. A cross-sectional case-control study was undertaken assessing patients with secondary progressive MS (SPMS) and age-matched controls. Co-registered magnetic resonance imaging was used to verify the location from which near-infrared spectroscopy data were obtained through Monte Carlo simulations of photon propagation. Additional clinical assessments of MS disease severity were carried out by trained neurologists. Linear mixed effect models were used to compare cortical oxygenation between cases and controls, and against measures of MS severity.

RESULTS

Thirty-three patients with secondary progressive MS (median expanded disability status scale 6 [IQR: 5-6.5]; median age 53.0 [IQR: 49-58]) and 20 age-matched controls were recruited. Modeling of photon propagation confirmed spectroscopy data were obtained from the prefrontal cortex. Patients with SPMS had significantly lower cortical hemoglobin oxygenation compared with controls (-6.0% [95% CI: -10.0 to -1.9], P = 0.004). There were no significant associations between cortical oxygenation and MS severity.

INTERPRETATION

Using an advanced, multiwavelength time-domain near-infrared spectroscopy system, we demonstrate that patients with SPMS have lower cortical oxygenation compared with controls.

Two-layered blood-lipid phantom and method to determine absorption and oxygenation employing changes in moments of DTOFs

Near-infrared spectroscopy (NIRS) is an established technique for measuring tissue oxygen saturation (StO₂), which is of high clinical value. For tissues that have layered structures, it is challenging but clinically relevant to obtain StO₂ of the different layers, e.g. brain and scalp. For this aim, we present a new method of data analysis for time-domain NIRS (TD-NIRS) and a new two-layered blood-lipid phantom. The new analysis method enables accurate determination of even large changes of the absorption coefficient (Δµ_a) in multiple layers. By adding Δµ_a to the baseline µ_a, this method provides absolute µ_a and hence StO₂ in multiple layers. The method utilizes (i) changes in statistical moments of the distributions of times of flight of photons (DTOFs), (ii) an analytical solution of the diffusion equation for an N-layered medium, (iii) and the Levenberg-Marquardt algorithm (LMA) to determine Δµ_a in multiple layers from the changes in moments. The method is suitable for NIRS tissue oximetry (relying on µ_a) as well as functional NIRS (fNIRS) applications (relying on Δµ_a). Experiments were conducted on a new phantom, which enabled us to simulate dynamic StO₂ changes in two layers for the first time. Two separate compartments, which mimic superficial and deep layers, hold blood-lipid mixtures that can be deoxygenated (using yeast) and oxygenated (by bubbling oxygen) independently. Simultaneous NIRS measurements can be performed on the two-layered medium (variable superficial layer thickness, L), the deep (homogeneous), and/or the superficial (homogeneous). In two experiments involving ink, we increased the nominal µ_a in one of two compartments from 0.05 to 0.25 cm^-1, L set to 14.5 mm. In three experiments involving blood (L set to 12, 15, or 17 mm), we used a protocol consisting of six deoxygenation cycles. A state-of-the-art multi-wavelength TD-NIRS system measured simultaneously on the two-layered medium, as well as on the deep compartment for a reference. The new method accurately determined µ_a (and hence StO₂) in both compartments. The method is a significant progress in overcoming the contamination from the superficial layer, which is beneficial for NIRS and fNIRS applications, and may improve the determination of StO₂ in the brain from measurements on the head. The advanced phantom may assist in the ongoing effort towards more realistic standardized performance tests in NIRS tissue oximetry. Data and MATLAB codes used in this study were made publicly available.

Multi-laboratory performance assessment of diffuse optics instruments: the BitMap exercise

SIGNIFICANCE

Multi-laboratory initiatives are essential in performance assessment and standardization-crucial for bringing biophotonics to mature clinical use-to establish protocols and develop reference tissue phantoms that all will allow universal instrument comparison.

AIM

The largest multi-laboratory comparison of performance assessment in near-infrared diffuse optics is presented, involving 28 instruments and 12 institutions on a total of eight experiments based on three consolidated protocols (BIP, MEDPHOT, and NEUROPT) as implemented on three kits of tissue phantoms. A total of 20 synthetic indicators were extracted from the dataset, some of them defined here anew.

APPROACH

The exercise stems from the Innovative Training Network BitMap funded by the European Commission and expanded to include other European laboratories. A large variety of diffuse optics instruments were considered, based on different approaches (time domain/frequency domain/continuous wave), at various stages of maturity and designed for different applications (e.g., oximetry, spectroscopy, and imaging).

RESULTS

This study highlights a substantial difference in hardware performances (e.g., nine decades in responsivity, four decades in dark count rate, and one decade in temporal resolution). Agreement in the estimates of homogeneous optical properties was within 12% of the median value for half of the systems, with a temporal stability of <5  %   over 1 h, and day-to-day reproducibility of <3  %  . Other tests encompassed linearity, crosstalk, uncertainty, and detection of optical inhomogeneities.

CONCLUSIONS

This extensive multi-laboratory exercise provides a detailed assessment of near-infrared Diffuse optical instruments and can be used for reference grading. The dataset-available soon in an open data repository-can be evaluated in multiple ways, for instance, to compare different analysis tools or study the impact of hardware implementations.

Lawrence K. Incorporating early and late-arriving photons to improve the reconstruction of cerebral hemodynamic responses acquired by time-resolved near-infrared spectroscopy

Significance: Despite its advantages in terms of safety, low cost, and portability, functional near-infrared spectroscopy applications can be challenging due to substantial signal contamination from hemodynamics in the extracerebral layer (ECL). Time-resolved near-infrared spectroscopy (tr NIRS) can improve sensitivity to brain activity but contamination from the ECL remains an issue. This study demonstrates how brain signal isolation can be further improved by applying regression analysis to tr data acquired at a single source–detector distance.

Aim: To investigate if regression analysis can be applied to single-channel trNIRS data to further isolate the brain and reduce signal contamination from the ECL.

Approach: Appropriate regressors for trNIRS were selected based on simulations, and performance was evaluated by applying the regression technique to oxygenation responses recording during hypercapnia and functional activation.

Results: Compared to current methods of enhancing depth sensitivity for trNIRS (i.e., higher statistical moments and late gates), incorporating regression analysis using a signal sensitive to the ECL significantly improved the extraction of cerebral oxygenation signals. In addition, this study demonstrated that regression could be applied to trNIRS data from a single detector using the early arriving photons to capture hemodynamic changes in the ECL.

Conclusion: Applying regression analysis to trNIRS metrics with different depth sensitivities improves the characterization of cerebral oxygenation signals.

Depth-selective data analysis for time-domain fNIRS: moments vs. time windows

Time-domain measurements facilitate the elimination of the influence of extracerebral, systemic effects, a key problem in functional near-infrared spectroscopy (fNIRS) of the adult human brain. The analysis of measured time-of-flight distributions of photons often relies on moments or time windows. However, a systematic and quantitative characterization of the performance of these measurands is still lacking. Based on perturbation simulations for small localized absorption changes, we compared spatial sensitivity profiles and depth selectivity for moments (integral, mean time of flight and variance), photon counts in time windows and their ratios for different time windows. The influence of the instrument response function (IRF) was investigated for all measurands and for various source-detector separations. Variance exhibits the highest depth selectivity among the moments. Ratios of photon counts in different late time windows can achieve even higher selectivity. An advantage of moments is their robustness against the shape of the IRF and instrumental drifts.

Characteristic of Mayer Waves in Electrophysiological, Hemodynamic and Vascular Signals

We evaluated the properties of oscillations in the Mayer waves (MW) frequency range (∼0.1Hz) detected in blood pressure, heart rate variability, cerebral blood oxygenation changes and evolution of electroencephalographic (EEG) rhythms to elucidate the mechanisms of MW generation. We examined the persistence of MW in different signals and stability of their oscillations on the level of individual MW waveforms, which was achieved by applying matching pursuit (MP). MP yields adaptive time-frequency approximation of signal's structures in terms of frequency, amplitude, time occurrence, and time-span. The number of waveforms contributing to 95% of the energy of the signals was vastly different for the time series, but the average number of waveforms conforming to the MW criteria was almost the same (3.5 ± 0.4 per 120s epoch). In all the investigated signals, MW had the same distributions of frequency and the number of cycles. We show that the MW energy ratios in different signals varied strongly, p < 0.001. The highest percentage of MW energy was observed in blood pressure signals, heart rate variability, and reduced hemoglobin, in contrast to brain signals and oxygenated hemoglobin. The percentage of MW energy was related to the strength of causal influence exerted by them on the other signals. Our results indicate existence of a common mechanism of MW generation and support the hypothesis of MW generation in the baroreflex loop; however, they do not exclude the action of a central pacemaker.

Depth-resolved assessment of changes in concentration of chromophores using time-resolved near-infrared spectroscopy: estimation of cytochrome-c-oxidase uncertainty by Monte Carlo simulations

Time-resolved near-infrared spectroscopy (TR-NIRS) measurements can be used to recover changes in concentrations of tissue constituents ( Δ C ) by applying the moments method and the Beer-Lambert law. In this work we carried out the error propagation analysis allowing to calculate the standard deviations of uncertainty in estimation of the Δ C . Here, we show the process of choosing wavelengths for the evaluation of hemodynamic (oxy-, deoxyhemoglobin) and metabolic (cytochrome-c-oxidase (CCO)) responses within the brain tissue as measured with an in-house developed TR-NIRS multi-wavelength system, which measures at 16 consecutive wavelengths separated by 12.5 nm and placed between 650 and 950 nm. Data generated with Monte Carlo simulations on three-layered model (scalp, skull, brain) for wavelengths range from 650 to 950 nm were used to carry out the error propagation analysis for varying choices of wavelengths. For a detector with a spectrally uniform responsivity, the minimal standard deviation of the estimated changes in CCO within the brain layer, σ Δ C CCO brain  = 0.40 µM, was observed for the 16 consecutive wavelengths from 725 to 912.5 nm. For realistic a detector model, i.e. the spectral responsivity characteristic is considered, the minimum, σ Δ C CCO brain  = 0.47 µM, was observed at the 16 consecutive wavelengths from 688 to 875 nm. We introduce the method of applying the error propagation analysis to data as measured with spectral TR-NIRS systems to calculate uncertainty of recovery of tissue constituents concentrations.

Influence of intra-abdominal pressure on the amplitude of fluctuations of cerebral hemoglobin concentration in the respiratory band

An intra-abdominal pressure (IAP) is correlated with cerebral perfusion, in a mechanism of reducing venous outflow. The elevated intra-abdominal pressure leads to an increase in the intracranial pressure and a decrease in the cerebral perfusion pressure. We studied the relationship between the IAP and the cerebral oxygenation with the use of the near infrared spectroscopy technique during a gynecological surgery. The changes in hemoglobin concentrations were analyzed in the time-frequency domain in the frequency band related to respiration. The measurements were carried out in 15 subjects who underwent laparoscopic surgery. During the laparoscopy, the intra-abdominal cavity was insufflated with CO₂, which caused a controlled increase in the IAP. It was observed that the amplitudes of respiration-related waves present in hemoglobin concentration signals show an increase of 1.5 to 8.5 times during elevation of the IAP by 15 mmHg.

Self-calibrating time-resolved near infrared spectroscopy

Time-resolved near infrared spectroscopy is considered to be a gold standard technique when measuring absolute values of tissue optical properties, as it provides separable and independent information about both tissue absorption and scattering. However, time-resolved instruments require an accurate characterization by measuring the instrument response function in order to decouple the contribution of the instrument itself from the measurement. In this work, a new approach to the methodology of analysing time-resolved data is presented where the influence of instrument response function is eliminated from the data and a self-calibrating analysis is proposed. The proposed methodology requires an instrument to provide at least two wavelengths and allows spectral parameters recovery (optical properties or constituents concentrations and reduced scatter amplitude and power). Phantom and in-vivo data from two different time-resolved systems are used to validate the accuracy of the proposed self-calibrating approach, demonstrating that parameters recovery compared to the conventional curve fitting approach is within 10% and benefits from introducing a spectral constraint to the reconstruction problem. It is shown that a multi-wavelength time-resolved data can be used for parameters recovery directly without prior calibration (instrument response function measurement).

Clinical Brain Monitoring with Time Domain NIRS: A Review and Future Perspectives

Near-infrared spectroscopy (NIRS) is an optical technique that can measure brain tissue oxygenation and haemodynamics in real-time and at the patient bedside allowing medical doctors to access important physiological information. However, despite this, the use of NIRS in a clinical environment is hindered due to limitations, such as poor reproducibility, lack of depth sensitivity and poor brain-specificity. Time domain NIRS (or TD-NIRS) can resolve these issues and offer detailed information of the optical properties of the tissue, allowing better physiological information to be retrieved. This is achieved at the cost of increased instrument complexity, operation complexity and price. In this review, we focus on brain monitoring clinical applications of TD-NIRS. A total of 52 publications were identified, spanning the fields of neonatal imaging, stroke assessment, traumatic brain injury (TBI) assessment, brain death assessment, psychiatry, peroperative care, neuronal disorders assessment and communication with patient with locked-in syndrome. In all the publications, the advantages of the TD-NIRS measurement to (1) extract absolute values of haemoglobin concentration and tissue oxygen saturation, (2) assess the reduced scattering coefficient, and (3) separate between extra-cerebral and cerebral tissues, are highlighted; and emphasize the utility of TD-NIRS in a clinical context. In the last sections of this review, we explore the recent developments of TD-NIRS, in terms of instrumentation and methodologies that might impact and broaden its use in the hospital.

Frequency analysis of oscillations in cerebral hemodynamics measured by time domain near infrared spectroscopy

In this paper, we propose the application of time-domain near-infrared spectroscopy to the assessment of oscillations in cerebral hemodynamics. These oscillations were observed in the statistical moments of the distributions of time of flight of photons (DTOFs) measured on the head. We analyzed the zeroth and second centralized moments of DTOFs (total number of photons and variance) to obtain their spectra to provide parameters for the frequency components of microcirculation, which differ between the extracerebral and intracerebral layers of the head. Analysis of these moments revealed statistically significant differences between a control group of healthy subjects and a group of patients with severe neurovascular disorders, which is a promising result for the assessment of cerebral microcirculation and cerebral autoregulation mechanisms.

Causal Coupling Between Electrophysiological Signals, Cerebral Hemodynamics and Systemic Blood Supply Oscillations in Mayer Wave Frequency Range

The aim of the study was to assess causal coupling between neuronal activity, microvascular hemodynamics and blood supply oscillations in the Mayer wave frequency range. An electroencephalogram, cerebral blood oxygenation changes, an electrocardiogram and blood pressure were recorded during rest and during a movement task. Causal coupling between them was evaluated using directed transfer function, a measure based on the Granger causality principle. The multivariate autoregressive model was fitted to all the signals simultaneously, which made it possible to construct a complete scheme of interactions between the considered signals. The obtained pattern of interactions in the resting state estimated in the 0.05-0.15 Hz band revealed a predominant influence of blood pressure oscillations on all the other variables. Reciprocal connections between blood pressure and heart rate variability time series indicated the presence of feedback loops between these signals. During movement, the pattern of connections did not change dramatically. The number of connections decreased, but the couplings between blood pressure and heart rate variability signal were not significantly changed, and the strong influence of the decreased blood hemoglobin concentration on the oxygenated hemoglobin concentration persisted. For the first time our results provided a comprehensive scheme of interactions between electrical and hemodynamic brain signals, heart rate and blood pressure oscillations. Persistent reciprocal connections between blood pressure and heart rate variability time series suggest possible feedforward and feedback coupling of cardiovascular variables which may lead to the observed oscillations in Mayer wave range.

Multiwavelength time-resolved near-infrared spectroscopy of the adult head: assessment of intracerebral and extracerebral absorption changes

An optical technique based on diffuse reflectance measurement combined with indocyanine green (ICG) bolus tracking is extensively tested as a method for the clinical assessment of brain perfusion at the bedside. We report on multiwavelength time-resolved diffuse reflectance spectroscopy measurements carried out on the head of a healthy adult during the intravenous administration of a bolus of ICG. Intracerebral and extracerebral changes in absorption were estimated from an analysis of changes in statistical moments (total number of photons, mean time of flight and variance) of the distributions of times of flight (DTOF) of photons recorded simultaneously at 16 wavelengths from the range of 650-850 nm using sensitivity factors estimated by diffusion approximation based on a layered model of the studied medium. We validated the proposed method in a series of phantom experiments and in-vivo measurements. The results obtained show that changes in the concentration of the ICG can be assessed as a function of time of the experiment and depth in the tissue. Thus, the separation of changes in ICG concentration appearing in intra- and extracerebral tissues can be estimated from optical data acquired at a single source-detector pair of fibers/fiber bundles positioned on the surface of the head.

Confirmation of brain death using optical methods based on tracking of an optical contrast agent: assessment of diagnostic feasibility

We aimed to determine whether optical methods based on bolus tracking of an optical contrast agent are useful for the confirmation of cerebral circulation cessation in patients being evaluated for brain death. Different stages of cerebral perfusion disturbance were compared in three groups of subjects: controls, patients with posttraumatic cerebral edema, and patients with brain death. We used a time-resolved near-infrared spectroscopy setup and indocyanine green (ICG) as an intravascular flow tracer. Orthogonal partial least squares-discriminant analysis (OPLS-DA) was carried out to build statistical models allowing for group separation. Thirty of 37 subjects (81.1%) were classified correctly (8 of 9 control subjects, 88.9%; 13 of 15 patients with edema, 86.7%; and 9 of 13 patients with brain death, 69.2%; p < 0.0001). Depending on the combination of variables used in the OPLS-DA model, sensitivity, specificity, and accuracy were 66.7–92.9%, 81.8–92.9%, and 77.3–89.3%, respectively. The method was feasible and promising in the demanding intensive care unit environment. However, its accuracy did not reach the level required for brain death confirmation. The potential usefulness of the method may be improved by increasing the depth of light penetration, confirming its accuracy against other methods evaluating cerebral flow cessation, and developing absolute parameters for cerebral perfusion.

Time-resolved near infrared light propagation using frequency domain superposition

Time-resolved temporal point spread function (TPSF) measurement of near infrared spectroscopic (NIRS) data allows the estimation of absorption and reduced scattering properties of biological tissues. Such analysis requires an iterative calculation of the theoretical TPSF curve using mathematical and computational models of the domain being imaged which are computationally complex and expensive. In this work, an efficient methodology for representing the TPSF data using a superposition of cosines calculated in frequency domain is presented. The proposed method is outlined and tested on finite element realistic models of the human neck and head. Using an adult head model containing ~140k nodes, the TPSF calculation at each node for one source is accelerated from 3.11 s to 1.29 s within an error limit of ± 5% related to the time domain calculation method.

Human skull translucency: post mortem studies

Measurements of optical translucency of human skulls were carried out. An incandescent light source and a CCD camera were used to measure the distribution of light transmitted through the skull in 10 subjects post-mortem. We noticed that intra-individual differences in optical translucency may be up to 100 times but inter-individual translucency differences across the skull reach 10⁵ times. Based on the measurement results, a "theoretical" experiment was simulated. Monte-Carlo calculations were used in order to evaluate the influence of the differences in optical translucency of the skull on results of NIRS measurements. In these calculations a functional stimulation was done, in which the oxyhemoglobin and deoxyhemoglobin concentrations in the brain cortex change by 5μM and -5μM respectively. The maximal discrepancies between assumed hemoglobin concentration changes and hemoglobin concentration changes estimated with Monte-Carlo simulation may reach 50% depending of the translucency of the skull.

Brain Oxygenation Monitoring

A mismatch between cerebral oxygen supply and demand can lead to cerebral hypoxia/ischemia and deleterious outcomes. Cerebral oxygenation monitoring is an important aspect of multimodality neuromonitoring. It is increasingly deployed whenever intracranial pressure monitoring is indicated. Although there is a large body of evidence demonstrating an association between cerebral hypoxia/ischemia and poor outcomes, it remains to be determined whether restoring cerebral oxygenation leads to improved outcomes. Randomized prospective studies are required to address uncertainties about cerebral oxygenation monitoring and management. This article describes the different methods of monitoring cerebral oxygenation, their indications, evidence base, limitations, and future perspectives.

Separation of superficial and cerebral hemodynamics using a single distance time-domain NIRS measurement

In functional near-infrared spectroscopy (fNIRS) superficial hemodynamics can mask optical signals related to brain activity. We present a method to separate superficial and cerebral absorption changes based on the analysis of changes in moments of time-of-flight distributions and a two-layered model. The related sensitivity factors were calculated from individual optical properties. The method was validated on a two-layer liquid phantom. Absorption changes in the lower layer were retrieved with an accuracy better than 20%. The method was successfully applied to in vivo data and compared to the reconstruction of homogeneous absorption changes.

Systematic review of near-infrared spectroscopy determined cerebral oxygenation during non-cardiac surgery

Near-infrared spectroscopy (NIRS) is used to monitor regional cerebral oxygenation (rScO₂) during cardiac surgery but is less established during non-cardiac surgery. This systematic review aimed (i) to determine the non-cardiac surgical procedures that provoke a reduction in rScO₂ and (ii) to evaluate whether an intraoperative reduction in rScO₂ influences postoperative outcome. The PubMed and Embase database were searched from inception until April 30, 2013 and inclusion criteria were intraoperative NIRS determined rScO₂ in adult patients undergoing non-cardiac surgery. The type of surgery and number of patients included were recorded. There was included 113 articles and evidence suggests that rScO₂ is reduced during thoracic surgery involving single lung ventilation, major abdominal surgery, hip surgery, and laparoscopic surgery with the patient placed in anti-Tredelenburg's position. Shoulder arthroscopy in the beach chair and carotid endarterectomy with clamped internal carotid artery (ICA) also cause pronounced cerebral desaturation. A >20% reduction in rScO₂ coincides with indices of regional and global cerebral ischemia during carotid endarterectomy. Following thoracic surgery, major orthopedic, and abdominal surgery the occurrence of postoperative cognitive dysfunction (POCD) might be related to intraoperative cerebral desaturation. In conclusion, certain non-cardiac surgical procedures is associated with an increased risk for the occurrence of rScO₂. Evidence for an association between cerebral desaturation and postoperative outcome parameters other than cognitive dysfunction needs to be established.

Multi-channel medical device for time domain functional near infrared spectroscopy based on wavelength space multiplexing

We have designed a compact dual wavelength (687 nm, 826 nm) multi-channel (16 sources, 8 detectors) medical device for muscle and brain imaging based on time domain functional near infrared spectroscopy. The system employs the wavelength space multiplexing approach to reduce wavelength cross-talk and increase signal-to-noise ratio. System performances have been tested on homogeneous and heterogeneous tissue phantoms following specifically designed protocols for photon migration instruments. Preliminary in vivo measurements have been performed to validate the instrument capability to monitor hemodynamic parameters changes in the arm muscle during arterial occlusion and in the adult head during a motor task experiment.

The effect of temporal impulse response on experimental reduction of photon scatter in time-resolved diffuse optical tomography

New fast detector technology has driven significant renewed interest in time-resolved measurement of early photons in improving imaging resolution in diffuse optical tomography and fluorescence mediated tomography in recent years. In practice, selection of early photons results in significantly narrower instrument photon density sensitivity functions (PDSFs) than the continuous wave case, resulting in a better conditioned reconstruction problem. In this work, we studied the quantitative impact of the instrument temporal impulse response function (TIRF) on experimental PDSFs in tissue mimicking optical phantoms. We used a multimode fiber dispersion method to vary the system TIRF over a range of representative literature values. Substantial disagreement in PDSF width--by up to 40%--was observed between experimental measurements and Monte Carlo (MC) models of photon propagation over the range of TIRFs studied. On average, PDSFs were broadened by about 0.3 mm at the center plane of the 2 cm wide imaging chamber per 100 ps of the instrument TIRF at early times. Further, this broadening was comparable on both the source and detector sides. Results were confirmed by convolution of instrument TIRFs with MC simulations. These data also underscore the importance of correcting imaging PDSFs for the instrument TIRF when performing tomographic image reconstruction to ensure accurate data-model agreement.