Image reconstruction of oxidized cerebral cytochrome C oxidase changes from broadband near-infrared spectroscopy data

Sabrina Brigadoi; Phong Phan; David Highton; Samuel Powell; Robert J. Cooper; Jeremy C. Hebden; Martin Smith; Ilias Tachtsidis; Clare E. Elwell; Adam P. Gibson

doi:10.1117/1.NPh.4.2.021105

24 May 2017 Image reconstruction of oxidized cerebral cytochrome C oxidase changes from broadband near-infrared spectroscopy data

Sabrina Brigadoi, Phong Phan, David Highton, Samuel Powell, Robert J. Cooper, Jeremy C. Hebden, Martin Smith, Ilias Tachtsidis, Clare E. Elwell, Adam P. Gibson

Author Affiliations +

Neurophotonics, 4(2), 021105 (2017). https://doi.org/10.1117/1.NPh.4.2.021105

Abstract

In diffuse optical tomography (DOT), overlapping and multidistance measurements are required to reconstruct depth-resolved images of oxy- ( HbO 2 ) and deoxy- (HHb) hemoglobin concentration changes occurring in the brain. These can be considered an indirect measure of brain activity, under the assumption of intact neurovascular coupling. Broadband systems also allow changes in the redox state of cytochrome c oxidase (oxCCO) to be measured, which can be an important biomarker when neurovascular coupling is impaired. We used DOT to reconstruct images of Δ [ HbO 2 ] , Δ [ HHb ] , and Δ [ oxCCO ] from data acquired with a broadband system. Four healthy volunteers were measured while performing a visual stimulation task (4-Hz inverting checkerboard). The broadband system was configured to allow multidistance and overlapping measurements of the participants’ visual cortex with 32 channels. A multispectral approach was employed to reconstruct changes in concentration of the three chromophores during the visual stimulation. A clear and focused activation was reconstructed in the left occipital cortex of all participants. The difference between the residuals of the three-chromophore model and of the two-chromophore model (recovering only Δ [ HbO 2 ] and Δ [ HHb ] ) exhibits a spectrum similar to that of oxCCO. These results form a basis for further studies aimed to further optimize image reconstruction of Δ [ oxCCO ] .

1. Introduction

Near-infrared spectroscopy (NIRS) uses light in the red and near-infrared range to monitor concentration changes of oxy- ( ${HbO}_{2}$ ) and deoxy- (HHb) hemoglobin in the brain.¹ Under the assumption of an intact neurovascular coupling, these can be considered as an indirect measure of brain activity. When multiple sources and detectors are used, arranged in overlapping and multidistance channels, depth-resolved images of these cerebral hemoglobin variations can be reconstructed.²^,³ This technique is usually referred to as diffuse optical tomography (DOT). Eggebrecht et al.² have recently shown that by using a high-density DOT system, DOT can be employed to map distributed brain networks and function. Their maps disclosed a strong correspondence with the same maps obtained with functional magnetic resonance imaging (fMRI). Their results are extremely promising since DOT techniques are characterized by noninvasiveness and portability. They can be easily employed to monitor vulnerable subjects, such as neonates,⁴5.^–⁶ patients with implanted devices,² or critically ill patients who need bedside monitoring⁷ and who could not be previously monitored with fMRI. Both functional NIRS (fNIRS) and fMRI rely on the assumption of intact neurovascular coupling, although this assumption is frequently violated in cases of brain pathology.⁸9.10.^–¹¹ Pathology (e.g., acute brain injury, cerebral ischemia, and neurodegenerative disease) can alter the regional hemodynamic response to localized changes in metabolism and in extreme cases can lead to absence or inversion of the response. Although it is well recognized that an inverted fNIRS hemodynamic response [or an fMRI negative blood-oxygen-level dependent (BOLD) response] can be seen in both healthy volunteers and after brain injury, the exact physiological mechanisms underlying these responses remain unclear. Additional information about colocalized cerebral metabolic state, such as that provided by oxCCO, can help elucidate these pathways. Therefore, there is an unmet clinical need for a direct marker of cerebral metabolism.

Cytochrome c oxidase (CCO), the terminal enzyme in the mitochondrial respiratory chain, is responsible for more than 95% of cellular oxygen metabolism.¹ The redox state of CCO can be measured with broadband NIRS¹² and as such represents an opportunity to directly measure cellular oxidative metabolism. The change in the oxidation state of CCO ( $Δ [oxCCO]$ ) reflects the balance between electron acceptance originating from nicotinamide adenine dinucleotide and donation to oxygen to form water, driving the mitochondrial proton electrochemical gradient and $> 90 %$ of adenosine triphosphate production.¹³ The relationship between cerebral oxygen and metabolic substrate delivery and CCO has been extensively investigated in animals,¹⁴ neonatal asphyxia,¹⁵ adult brain injury,¹⁶ and healthy adults,¹⁷ including during functional activation.¹⁸ Importantly, $Δ [oxCCO]$ is concordant with measurements of metabolism, including invasively measured lactate/pyruvate ratio¹⁶^,¹⁹ and magnetic resonance spectroscopy-derived lactate.¹⁵ However, because the concentration of CCO is an order of magnitude lower than that of hemoglobin, an optimized optical technique to measure $Δ [oxCCO]$ using broadband NIRS²⁰^,²¹ is typically required.

Present commercial fNIRS systems do not have the capability to measure $Δ [oxCCO]$ , and in-house systems developed explicitly to measure $Δ [oxCCO]$ provide only a limited number of channels because of hardware cost and complexity, and therefore cannot perform DOT.¹³ Reconstructing images of $Δ [oxCCO]$ is a natural progression of the technique following the development of DOT for displaying $Δ [{HbO}_{2}]$ and $Δ [HHb]$ . Recovering the spatial distribution of this additional chromophore may well facilitate the investigation of cerebral energy status alongside functional activation, particularly in the context of brain pathology.²² However, the small $Δ [oxCCO]$ signal that occurs in the presence of larger changes in $Δ [{HbO}_{2}]$ and $Δ [HHb]$ means that imaging it is a significant challenge.

The recovery of images of concentration changes via DOT constitutes a nonlinear, ill-posed inverse problem.²³24.^–²⁵ This process involves three steps. First, a geometric representation of the target object is required. A subject-specific, multilayer, anatomically accurate head model is preferred, but when the magnetic resonance image of the subject is unavailable, a multilayer atlas head model can be employed.²⁶ The atlas head model can be registered to the subject’s cranial landmarks, and wavelength- and tissue-specific optical properties can be assigned to each tissue type [e.g., scalp, skull, cerebrospinal fluid (CSF), gray matter (GM), and white matter (WM)], thus providing an anatomically meaningful solution space. Second, a model of how light propagates through these head tissues, given an array of sources and detectors, is required, so as to relate changes in optical properties to changes in the acquired signals. This forward model can be constructed either with finite-element method (FEM) approaches (e.g., Toast++²⁷ and Nirfast²⁸) or with Monte Carlo approaches (e.g., MMC²⁹ and MCX³⁰). Third, the solution of the forward model must be inverted so as to solve the inverse problem and reconstruct an image of concentration changes. This step is the most computationally expensive, since it is an underdetermined and ill-posed problem. Corlu et al.³¹ demonstrated a multispectral approach, which directly reconstructs changes in concentration (rather than absorption changes at each wavelength). This approach reduces the number of unknowns and better constrains the inverse problem.

Most DOT studies to date have recovered images of $Δ [{HbO}_{2}]$ , $Δ [HHb]$ , and of changes in total hemoglobin concentration $Δ [HbT]$ (where $Δ [HbT] = Δ [{HbO}_{2}] + Δ [HHb])$ , since these are the most studied chromophores in clinical and neuroscience settings.²^,⁴^,⁵ Other chromophores, such as water and lipids, have been successfully reconstructed with DOT,³¹^,³² mostly for breast imaging applications.

The aim of this work is to provide, for the first time, a proof of concept demonstration of the feasibility of reconstructing images of oxCCO concentration changes, along with $Δ [{HbO}_{2}]$ and $Δ [HHb]$ , in the adult brain. Data were acquired from four healthy volunteers with a broadband system during a visual stimulation paradigm, and images were reconstructed with a multispectral approach. Visual stimulation was employed since it is known to elicit a strong localized activation with a highly repeatable pattern.³³^,³⁴ Limitations and assumptions of the present work will be discussed.

2. Material and Methods

2.1.

Instrumentation

A detailed description of the broadband system used in this study can be found elsewhere.¹⁸ The system has two 50-W halogen light sources, each delivering light to the participant via optical fibers. Only one source can be illuminated at a time. To collect the diffusely reflected light and measure the intensity at each wavelength, the system has eight detector fibers divided into two bundles and two CCD cameras (Pixis 512, Princeton Instruments, Trenton, New Jersey), each connected to a spectrograph. The sampling frequency is 1 Hz.

The optical interface between the system and the participant was based on a three-dimensional (3-D)-printed optode holder. The optode holder incorporated 4 source slots (in the middle row) and 14 detector slots arranged in 2 rows (see Fig. 1).

Fig. 1

Array layout. Sources were positioned in the central row (red dots), while detectors were located in the bottom and top rows (blue dots). Eight channels per source could be recorded, for a total of 32 channels. Note that not all slots could be filled with a source or a detector during a single acquisition, since only one source and eight detectors were available at a time. Each acquisition was repeated four times, with the available source and detector fibers positioned in different slots for each acquisition so as to cover the entire imaging array over four acquisitions in each subject.

2.2.

Participants and Paradigm

Four healthy volunteers (two males and two females; age range 27 to 48 years old) participated in the visual stimulation study after providing written informed consent. Participants were seated in a comfortable chair in a dimly lit room and were asked to focus their attention at the center of the computer monitor, located at a distance of 60 cm in front of the subject. Visual functional activation was achieved with a 4-Hz inverting checkerboard covering the full visual field. The paradigm consisted of 20 s of stimulation followed by 20 s of blank screen repeated over 10 epochs. To monitor the visual cortex, the optode array was fixed horizontally to cover the participant’s left occipital cortex, with the fourth source location positioned over Oz (according to the 10-20 EEG positions) (Fig. 2). The eight nearest detector slots in the array (four above and four below each activated source) were populated with the eight available detector fibers and were used to collect intensity data continuously during the 10 epochs of visual stimulation. The remaining detector slots remained empty during that acquisition. This was repeated for each of the four source slots by translating the fibers over the fixed optode array on the head and repeating the stimulation. A 3-D Patriot^™ Digitizer (Polhemus, Colchester, Vermont) was employed to acquire, for each participant, the coordinates of their cranial landmarks (nasion, inion, preauricular points, and Cz) in addition to the optode locations. The study was approved by University College London's Research Ethics Committee (1133/001).

Fig. 2

Optode holder placed on the left occipital cortex with one source (central row) and eight detectors.

2.3.

Data Analysis

Attenuation changes ( $Δ A_{λ}$ ) were calculated from the intensity changes across each of the wavelengths from 740 to 900 nm collected during the challenge using the following equation:

(1)

Δ A_{λ_{i}} = \log_{10} (\frac{I_{0}}{I_{λ_{i}}}),

where

I_{0}

is an arbitrary reference intensity (set to

10^{5}

) and

I_{λ i}

is the intensity collected by the detector at wavelength

λ_{i}

.

The 400-s trace of attenuation changes in each channel and for each wavelength was bandpass filtered (fifth-order Butterworth filter, cutoff frequencies: 0.01 to 0.25 Hz) and then block-averaged according to the timestamp of the stimulation to produce a 40-s trace of mean attenuation changes.

Channel-wise $Δ [{HbO}_{2}]$ , $Δ [HHb]$ , and $Δ [oxCCO]$ concentrations were computed by applying the Beer–Lambert law as implemented in the UCLn algorithm,²⁰ as is commonly performed in spectroscopy. The wavelength dependence of the differential pathlength factor, as computed by Kolyva et al.,¹⁸ was taken into account when solving the Beer–Lambert law. These channel-wise data were recovered to show the results of the standard analysis employed in spectroscopy and compare them with the recovered images.

To reduce the computational burden of the reconstruction while covering the available near-infrared spectrum, attenuation changes measured at 17 discrete wavelengths (from 740 to 900 nm at intervals of 10 nm) were derived from the measured broadband data to perform the reconstruction.

2.4.

Image Reconstruction Procedure

2.4.1.

Head model

The nonlinear MNI-ICBM152 atlas³⁵ was the basis of our head model, which was built as described in Brigadoi and Cooper³⁶ and Dempsey et al.³⁷ A multilayer tissue mask was created by segmenting the MR image into five tissue types (scalp, skull, CSF, GM, and WM). A high-resolution tetrahedral mesh ( $\sim 2 \times 10^{5}$ nodes and $\sim 10^{6}$ elements) was created with the iso2mesh software.³⁸ The same software was used to create the GM surface mesh ( $\sim 3 \times 10^{4}$ nodes and $\sim 6 \times 10^{4}$ faces), which was used to display the reconstructed images.

The head model was registered to each participant’s head with an affine registration using the measured cranial landmarks as reference. The affine registration was also used to register the GM surface mesh.

2.4.2.

Inverse problem and image reconstruction

Images representing $Δ [{HbO}_{2}]$ , $Δ [HHb]$ , and $Δ [oxCCO]$ were reconstructed using a multispectral approach, which directly generates concentration changes from attenuation data.³¹ The multispectral approach is a standard approach,³¹^,³⁷ which has been shown to improve the accuracy of the reconstructed images as the number of wavelengths increases, as in our case. After linearizing the forward problem,²³ we may relate changes in measured optical attenuation to changes in concentration through the multispectral Jacobian matrix

(2)

[\begin{matrix} Δ A_{λ_{1}}^{1} \\ Δ A_{λ_{1}}^{2} \\ \dots \\ Δ A_{λ_{1}}^{n} \\ Δ A_{λ_{2}}^{1} \\ \dots \\ Δ A_{λ_{17}}^{n} \end{matrix}] = [\begin{matrix} J_{λ_{1}} ϵ_{λ_{1}}^{{HbO}_{2}} & J_{λ_{1}} ϵ_{λ_{1}}^{HHb} & J_{λ_{1}} ϵ_{λ_{1}}^{oxCCO} \\ J_{λ_{2}} ϵ_{λ_{2}}^{{HbO}_{2}} & J_{λ_{2}} ϵ_{λ_{2}}^{HHb} & J_{λ_{2}} ϵ_{λ_{2}}^{oxCCO} \\ \dots \\ J_{λ_{17}} ϵ_{λ_{17}}^{{HbO}_{2}} & J_{λ_{17}} ϵ_{λ_{17}}^{HHb} & J_{λ_{17}} ϵ_{λ_{17}}^{oxCCO} \end{matrix}] [\begin{matrix} Δ C_{{HbO}_{2}} \\ Δ C_{HHb} \\ Δ C_{oxCCO} \end{matrix}],

where

ϵ_{λ}

is the specific absorption coefficient of

{HbO}_{2}

, HHb, and oxCCO at wavelength

λ_{i}

,

J_{λ}

is the wavelength-specific Jacobian matrix (number of channels × number of nodes),

Δ C

is the concentration change of

{HbO}_{2}

, HHb, and oxCCO at each node, and

Δ A

is the attenuation change between active and rest state at each wavelength

λ_{i}

and at each channel (from 1 to n).

Each wavelength-specific Jacobian was computed with the Toast++ software,²⁷ which solves the diffusion approximation via the FEM. The Jacobian was computed on the volumetric head mesh. Optical properties (absorption coefficient and scattering coefficient) were assigned to each tissue type and for each wavelength by fitting all published values for these tissue types.³⁹40.^–⁴¹

The forward model solution was projected by Toast++ onto a $50 \times 60 \times 50$ voxel regular grid for reconstruction, and an intermediate finer grid of $100 \times 120 \times 100$ voxels was used to optimize the mapping between the mesh space and the voxel space. Diffuse boundary sources and detectors were simulated as a Gaussian profile with a 2-mm standard deviation, and Neumann boundary conditions were applied.

To solve the changes in concentration in a least squares-sense, we minimized the form

(3)

E (Δ C) = {‖ Δ A - J Δ C ‖}_{Γ_{e}^{- 1}}^{2} + λ R (Δ C),

where

R (Δ C)

is a suitable regularization functional whose influence is weighted by the hyperparameter

λ

, and

Γ_{e}^{- 1}

is the inverse noise covariance matrix (set here as identity matrix).

The regularization functional $R (Δ C)$ imposes prior knowledge of the solution to overcome the ill-posedness of the inverse problem. In this work, we chose to employ first-order Tikhonov regularization such that

(4)

R (Δ C) = {‖ \nabla (Δ C) ‖}_{Γ_{x}^{- 1}}^{2},

with the inverse parameter covariance matrix

(5)

Γ_{x}^{- 1} = [\begin{matrix} {\begin{matrix} (C_{{HbO}_{2}}^{b}) \end{matrix}}^{- 2} & 0 & 0 \\ 0 & {\begin{matrix} (C_{HHb}^{b}) \end{matrix}}^{- 2} & 0 \\ 0 & 0 & {\begin{matrix} (C_{oxCCO}^{b}) \end{matrix}}^{- 2} \end{matrix}],

in which

C^{b}

is the square matrix (number of voxels × number of voxels) with the background concentration values of the three chromophores: 56 for

{HbO}_{2}

, 24 for HHb, and 12.8 for oxCCO⁴² on their diagonal.

C^{b}

for oxCCO was set to 12.8 (instead of the suggested 6.4), given the higher concentration of oxCCO in the visual cortex.⁴³ Use of the parameter covariance matrix spheres the solution space, thus ensuring that regularization is applied equally to each of the parameters, irrespective of their background concentrations. This is analogous to the normalization usually employed when absorption and scattering coefficients have to be simultaneously recovered.²³^,⁴⁴ We employed first-order Tikhonov regularization to avoid suppressing the solution, thus ensuring an unbiased estimation of the three chromophore concentrations.

The LSQR method⁴⁵ was employed to solve the matrix equations resulting from the minimization. For each reconstruction, the LSQR algorithm was configured to terminate after a maximum of 50 iterations, and with a tolerance of $10^{- 5}$ . The regularization hyperparameter was set to $λ = 10^{- 4}$ by inspection, after examining all images recovered varying $λ$ from $10^{- 2}$ to $10^{- 5}$ , and further confirmed by an L-curve analysis, as the highest point of curvature of the L-curve in each participant. The LSQR method is an iterative method, which has been proven to be analytically equivalent to the more well-known conjugate gradient, but with better numerical properties.⁴⁵ Furthermore, it is a suitable method to solve the inverse problem when the Jacobian matrix is large and sparse, as in our case. By employing the LSQR method and a preconditioning technique,²³ we could also reduce the amount of memory required by the system to solve the inverse problem, when compared to the standard Moore–Penrose inverse technique.²⁴

The reconstructed images were defined on the same regular grid as the Jacobian. A remapping procedure was therefore performed to map back the reconstructed images from the voxel space to the volumetric head mesh. The volumetric head mesh-based image was then projected to the GM surface mesh by assigning a value to each node on the GM boundary surface that was equal to the mean value of all the volumetric mesh node values within a 3-mm radius.

For comparison purposes, the same procedure described above was employed to reconstruct images of $Δ [{HbO}_{2}]$ and $Δ [HHb]$ only. In this case, a two-chromophore model was employed in the reconstruction, instead of the three-chromophore model (with ${HbO}_{2}$ , HHb, and oxCCO) used previously.

2.5.

Residual Computation

Residuals were computed to investigate the spectral changes due to the oxCCO chromophore. Channel-wise attenuation changes were deduced from the concentration changes in the voxel space (i.e., before their mapping to the volumetric mesh and then to the GM surface mesh, which was performed for visualization purposes only). Attenuation changes were determined for both the concentration changes computed using the three-chromophore model and the changes computed using the two-chromophore model by multiplying the concentration changes voxel-wise data with the multispectral Jacobian.

The difference between the three- and the two-chromophore fits was then investigated. Ideally, if the model explains the data, the residuals should be multivariate Gaussian noise with covariance $Γ_{e}$ , i.e., the model accounts for all nonrandom variability in the data. Given that image reconstruction is an ill-posed inverse problem that requires regularization, it is likely that even if the model explains the data, the residuals would not equate to the covariance of the assumed noise model. However, a model that better explains the data will have residuals that tend toward Gaussian noise. If the difference between the residuals of the two models has a nonrandom distribution across wavelengths, with the residuals of the three-chromophore model smaller than those of the two-chromophore models, this would suggest that the two-chromophore model is not accounting for a chromophore that is evident in the data. This analysis of residuals has been previously used in spectroscopy studies recovering $Δ [{HbO}_{2}]$ , $Δ [HHb]$ , and $Δ [oxCCO]$ .¹⁵^,¹⁸

Residuals were computed for each channel at the time of the peak of the expected functional activation (i.e., 20 s).

All analyses and reconstruction procedures were performed in MATLAB^® (Mathworks, Natick, Massachusetts). Figure 3 shows the data analysis and image reconstruction procedures.

Fig. 3

Flow diagram summarizing the main steps of the data analysis and image reconstruction procedures.

3. Results

In Fig. 4, examples of reconstructed $Δ [{HbO}_{2}]$ , $Δ [HHb]$ , and $Δ [oxCCO]$ images on the GM surface mesh at five different time points are displayed for one representative participant (participant 1). The selected time points are 2, 10, 20, 30, and 40 s after stimulus presentation. Figure 5 shows the images obtained at 20-s poststimulus onset for the other three participants.

Fig. 4

Examples of reconstructed $Δ [{HbO}_{2}]$ (first row), $Δ [HHb]$ (second row), and $Δ [oxCCO]$ (third row) images on the GM surface mesh at five different time points for participant 1. The selected time points are 2, 10, 20, 30, and 40 s after stimulus presentation. Images were mapped to the MNI152 GM mesh for visualization purposes.

Fig. 5

Examples of reconstructed $Δ [{HbO}_{2}]$ (first row), $Δ [HHb]$ (second row), and $Δ [oxCCO]$ (third row) images on the GM surface mesh for participant 2, 3, and 4. The reconstruction for 20-s poststimulus onset is displayed.

Figure 6 shows the grand-average channel-wise recovered $Δ [{HbO}_{2}]$ , $Δ [HHb]$ , and $Δ [oxCCO]$ changes obtained using the spectroscopy approach over the whole 740- to 900-nm spectrum. Only some of the channels exhibit a strong activation (defined as an increase in ${HbO}_{2}$ with a concomitant decrease in HHb), which is in line with the focal and localized change in chromophores reconstructed in Figs. 4 and 5.

Fig. 6

Grand-average hemodynamic responses in each channel to the functional activation task. Error bars represent the standard deviation across participants. Only the upper error bar for ${HbO}_{2}$ and the lower error bar for HHb are reported for visualization purposes. $Δ [{HbO}_{2}]$ is reported in red, $Δ [HHb]$ in blue, and $Δ [oxCCO]$ in green. These responses were obtained using the data from the whole broadband spectrum (740 to 900 nm) and the UCLn algorithm.

Figure 7 shows the grand-average difference between the residuals of the two- and three-chromophore models in all channels. Residuals were computed as the differences between the measured data and the back-projected attenuation changes after image reconstruction. The spectra in most of the active channels do not have a random distribution. Instead, they show a shape approximating the spectrum of the oxCCO chromophore, with a broad peak at around 810 nm. In most of the channels that do not show activation, the difference between the residuals of the two models is almost nil. There is some variability in the residuals among participants (as there is in the hemodynamic responses, see Fig. 6) likely due to: (1) each participant can have a different timing of brain activation (all residuals are computed at time point 20 s) and (2) the retinotopic mapping of the visual cortex is different in each person, increasing the chance that each channel is probing slightly different brain visual areas in each participant.⁴⁶

Fig. 7

Grand-average residual difference between the two- and three-chromophore model in each channel and for each wavelength at time point 20. Error bars represent the standard deviation across participants.

4. Discussion

In this study, we explored the feasibility of reconstructing images of $Δ [oxCCO]$ , in addition to $Δ [{HbO}_{2}]$ and $Δ [HHb]$ , using data acquired with a broadband NIRS system. This NIRS system allowed us to measure broadband data from overlapping and multidistance channels. With this configuration, we were able to perform DOT on broadband data and reconstruct images of $Δ [oxCCO]$ for the first time. Data were acquired in four healthy volunteers during a visual stimulation task, which is known to elicit a strong activation, with a highly reproducible response.³³^,³⁴

We evaluated the reconstructed images both qualitatively and quantitatively. All four participants showed a clear and focused activation for the three chromophores in the left occipital cortex during the period of maximum response to the visual stimulation. These results are in line with the channel-wise spectroscopy data, which show a clear response to the task in a subset of channels. Furthermore, by employing a two-chromophore model, we reconstructed from the same data images of $Δ [{HbO}_{2}]$ and $Δ [HHb]$ only. We then computed the residuals for both the two- and the three-chromophore models as the difference between the original attenuation changes and the back-projected attenuation changes. The difference between these residuals provides insight on whether a third chromophore is apparent in the data. The residual difference spectra exhibited a shape similar to that of the oxCCO spectrum, a common finding in NIRS studies, thus suggesting that the functionally induced changes in oxCCO concentration are influencing our data and can be imaged via DOT.

It should be noted that while $Δ [oxCCO]$ , $Δ [{HbO}_{2}]$ , and $Δ [HHb]$ are normally recovered from NIRS data assuming a linear relationship between the changes in concentration and attenuation, in DOT, a nonlinear, ill-posed, underdetermined inverse problem must be solved. Therefore, the residuals in DOT should be interpreted with caution. Solving the inverse problem involves a balance between our trust in the data and our confidence in the prior knowledge. Residuals can be influenced by the fact that the problem is ill-posed and underdetermined, and by uncertainties in the prior knowledge. Indeed, the back-projected attenuation changes are unlikely to be identical to the measurements. Regularization is required to improve the conditioning of the inverse problem and to cope with noise amplification that can occur in the inversion process due to small singular values in the Jacobian matrix. First-order Tikhonov regularization was utilized here to avoid suppressing the mean value of the images and to prevent the build-up of excessive high-frequency components, which cannot be resolved using diffuse optical methods. By sphering the solution space prior to reconstruction, we ensured that the effect of the regularization was applied equally to each concentration value.

The optimal wavelength range for NIRS studies of functional changes in tissue is still under debate, 740 to 900 nm has been investigated,³⁴ with preference for 780 to 900 nm to recover $Δ [oxCCO]$ , $Δ [{HbO}_{2}]$ , and $Δ [HHb]$ ¹⁵^,¹⁸ because of evidence of distortion over the peak of the HHb extinction spectrum between 740 and 780 nm. Arifler et al.²¹ reported the best combinations of wavelengths to recover $Δ [oxCCO]$ , $Δ [{HbO}_{2}]$ , and $Δ [HHb]$ in spectroscopy when the full broadband spectrum cannot be acquired. With a combination of three wavelengths, the mean recovery error versus the gold standard (i.e., using the full spectrum) was 10%, which was reduced to 2% when using a combination of eight wavelengths. In all combinations, the lowest selected wavelength was around 784 to 785 nm. However, previous DOT studies aimed at recovering only $Δ [{HbO}_{2}]$ , and $Δ [HHb]$ found that the best combination of wavelengths contained at least one wavelength below 780 nm, where measurements are more sensitive to HHb.³⁷^,⁴² Further studies should be performed to evaluate which wavelength combination is optimum to recover all three chromophores with DOT and to further investigate this issue. In this study, we selected 17 wavelengths from the broadband spectrum, separated at intervals of 10 nm. This choice was a compromise between the computational burden and the accuracy of the reconstructed image³²^,⁴⁷ when using the multispectral approach proposed by Corlu et al.³¹ The computation time for the inversion process for one frame was $\sim 101 s$ with an Intel^® Core^™ processor i7-4790 at 3.60 GHz.

Previous studies have shown the feasibility of recovering DOT images of $Δ [{HbO}_{2}]$ and $Δ [HHb]$ in real time.²⁵^,⁴⁸ Real-time imaging could provide valuable clinical information, which can be more naturally interpreted by clinicians. Considering the clinical importance of oxCCO, delivering real-time imaging of $Δ [oxCCO]$ , $Δ [{HbO}_{2}]$ , and $Δ [HHb]$ would be a great advance in the optical imaging field. Although from a theoretical perspective this should be readily achievable, there will likely be challenges to overcome when recovering images of all three chromophores in real time, mainly due to the higher number of wavelengths required and the associated increase in computational burden. Further studies should investigate the optimal number of wavelengths that can reliably yield $Δ [oxCCO]$ images, and whether the computational burden of the reconstruction could be reduced to allow images of $Δ [oxCCO]$ , $Δ [{HbO}_{2}]$ , and $Δ [HHb]$ to be delivered in real time, perhaps using a dynamic Kalman filter estimator.⁴⁹^,⁵⁰

This initial study will hopefully promote further research to optimize image reconstruction methods to deliver oxCCO changes and inform instrument development. Further work is required to explore the validity of the assumptions made in this study (e.g., the prior knowledge on the model covariance matrix, the selected wavelengths, the settings employed to solve the inverse problem, etc.), and the image reconstruction would benefit from further validation with simulations, and with phantom and in vivo studies. Ideally, the procedure should be verified using pathological conditions that induce changes in oxCCO but not in ${HbO}_{2}$ or HHb concentration, which would reveal evidence of cross talk during the reconstruction of the three chromophores. Perhaps this could be achieved using an animal model of cardiopulmonary bypass with fluorocarbon perfusion.⁵¹^,⁵²

In this study, we have demonstrated for the first time the reconstruction of images of simultaneous changes in concentration of oxCCO, ${HbO}_{2}$ , and HHb in the healthy adult brain during a visual stimulation task. A clear and focal activation was reconstructed in the left occipital cortex for all three chromophores in all participants. Reconstructing images of $Δ [oxCCO]$ may provide a unique opportunity to investigate cerebral energy status in addition to functional activation, across the cortex, which may have particular implications for studies of brain pathology.

Disclosures

The authors declare no competing financial interests.

Acknowledgments

This work was supported by the Engineering and Physical Sciences Research Council (EPSRC) Grant No. EP/K020315/1. R.J.C. was supported by EPSRC fellowship EP/N025946/1. I.T. was supported by the Wellcome Trust (104580/Z/14/Z).

References

1. F. F. Jöbsis, “Noninvasive, infrared monitoring of cerebral and myocardial oxygen sufficiency and circulatory parameters,” Science 198, 1264–1267 (1977).SCIEAS0036-8075 http://dx.doi.org/10.1126/science.929199 Google Scholar

2. A. T. Eggebrecht et al., “Mapping distributed brain function and networks with diffuse optical tomography,” Nat. Photonics 8, 448–454 (2014).NPAHBY1749-4885 http://dx.doi.org/10.1038/nphoton.2014.107 Google Scholar

3. S. R. Arridge, “Methods in diffuse optical imaging,” Philos. Trans. R. Soc. A 369, 4558–4576 (2011). http://dx.doi.org/10.1098/rsta.2011.0311 Google Scholar

4. H. Singh et al., “Mapping cortical haemodynamics during neonatal seizures using diffuse optical tomography: a case study,” NeuroImage Clin. 5, 256–265 (2014). http://dx.doi.org/10.1016/j.nicl.2014.06.012 Google Scholar

5. M. Chalia et al., “Hemodynamic response to burst-suppressed and discontinuous electroencephalography activity in infants with hypoxic ischemic encephalopathy,” Neurophotonics 3, 031408 (2016). http://dx.doi.org/10.1117/1.NPh.3.3.031408 Google Scholar

6. S. L. Ferradal et al., “Functional imaging of the developing brain at the bedside using diffuse optical tomography,” Cereb. Cortex 26(4), 1558–1568 (2016). http://dx.doi.org/10.1093/cercor/bhu320 Google Scholar

7. M. Smith, “Shedding light on the adult brain: a review of the clinical applications of near-infrared spectroscopy,” Philos. Trans. R. Soc. London A 369, 4452–4469 (2011). http://dx.doi.org/10.1098/rsta.2011.0242 Google Scholar

8. W. B. Baker et al., “Neurovascular coupling varies with level of global cerebral ischemia in a rat model,” J. Cereb. Blood Flow Metab. 33, 97–105 (2013). http://dx.doi.org/10.1038/jcbfm.2012.137 Google Scholar

9. L. M. Carusone et al., “Hemodynamic response changes in cerebrovascular disease: implications for functional MR imaging,” AJNR Am. J. Neuroradiol. 23, 1222–1228 (2002). Google Scholar

10. N. K. Logothetis, “What we can do and what we cannot do with fMRI,” Nature 453, 869–878 (2008). http://dx.doi.org/10.1038/nature06976 Google Scholar

11. J. U. Blicher et al., “Visualization of altered neurovascular coupling in chronic stroke patients using multimodal functional MRI,” J. Cereb. Blood Flow Metab. 32, 2044–2054 (2012). http://dx.doi.org/10.1038/jcbfm.2012.105 Google Scholar

12. C. Cooper et al., “Measurement of cytochrome oxidase redox state by near infrared spectroscopy,” Adv. Exp. Med. Biol. 413, 63–73 (1997).AEMBAP0065-2598 http://dx.doi.org/10.1007/978-1-4899-0056-2 Google Scholar

13. G. Bale, C. E. Elwell and I. Tachtsidis, “From Jöbsis to the present day: a review of clinical near-infrared spectroscopy measurements of cerebral cytochrome-C-oxidase,” J. Biomed. Opt. 21, 091307 (2016).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.JBO.21.9.091307 Google Scholar

14. A. Bainbridge et al., “Brain mitochondrial oxidative metabolism during and after cerebral hypoxia-ischemia studied by simultaneous phosphorus magnetic-resonance and broadband near-infrared spectroscopy,” NeuroImage 102, 173–183 (2014).NEIMEF1053-8119 http://dx.doi.org/10.1016/j.neuroimage.2013.08.016 Google Scholar

15. G. Bale et al., “A new broadband near-infrared spectroscopy system for in-vivo measurements of cerebral cytochrome-c-oxidase changes in neonatal brain injury,” Biomed. Opt. Express 5, 3450–3466 (2014).BOEICL2156-7085 http://dx.doi.org/10.1364/BOE.5.003450 Google Scholar

16. A. Ghosh et al., “Hyperoxia results in increased aerobic metabolism following acute brain injury,” J. Cereb. Blood Flow Metab. (2016). http://dx.doi.org/10.1177/0271678X16679171 Google Scholar

17. C. Kolyva et al., “Cytochrome C oxidase response to changes in cerebral oxygen delivery in the adult brain shows higher brain-specificity than haemoglobin,” NeuroImage 85, 234–244 (2014).NEIMEF1053-8119 http://dx.doi.org/10.1016/j.neuroimage.2013.05.070 Google Scholar

18. C. Kolyva et al., “Systematic investigation of changes in oxidized cerebral cytochrome c oxidase concentration during frontal lobe activation in healthy adults,” Biomed. Opt. Express 3, 2550–2566 (2012).BOEICL2156-7085 http://dx.doi.org/10.1364/BOE.3.002550 Google Scholar

19. M. M. Tisdall et al., “Increase in cerebral aerobic metabolism by normobaric hyperoxia after traumatic brain injury,” J. Neurosurg. 109, 424–432 (2008).JONSAC0022-3085 http://dx.doi.org/10.3171/JNS/2008/109/9/0424 Google Scholar

20. S. J. Matcher et al., “Performance comparison of several published tissue near-infrared spectroscopy algorithms,” Anal. Biochem. 227, 54–68 (1995).ANBCA20003-2697 http://dx.doi.org/10.1006/abio.1995.1252 Google Scholar

21. D. Arifler et al., “Optimal wavelength combinations for near-infrared spectroscopic monitoring of changes in brain tissue hemoglobin and cytochrome C oxidase concentrations,” Biomed. Opt. Express 6, 933–947 (2015).BOEICL2156-7085 http://dx.doi.org/10.1364/BOE.6.000933 Google Scholar

22. P. Phan et al., “Spatial distribution of changes in oxidised cytochrome C oxidase during visual stimulation using broadband near infrared spectroscopy imaging,” Adv. Exp. Med. Biol. 923, 195–201 (2016). http://dx.doi.org/10.1007/978-3-319-38810-6_26 Google Scholar

23. S. Arridge and J. Schotland, “Optical tomography: forward and inverse problems,” Inverse Prob. 25(12) (2009). http://dx.doi.org/10.1088/0266-5611/25/12/123010 Google Scholar

24. A. P. Gibson, J. C. Hebden and S. R. Arridge, “Recent advances in diffuse optical imaging,” Phys. Med. Biol. 50, R1–R43 (2005).PHMBA70031-9155 http://dx.doi.org/10.1088/0031-9155/50/4/R01 Google Scholar

25. S. Brigadoi et al., “Evaluating real-time image reconstruction in diffuse optical tomography using physiologically realistic test data,” Biomed. Opt. Express 6, 4719–4737(2015).BOEICL2156-7085 http://dx.doi.org/10.1364/BOE.6.004719 Google Scholar

26. R. J. Cooper et al., “Validating atlas-guided DOT: a comparison of diffuse optical tomography informed by atlas and subject-specific anatomies,” NeuroImage 62, 1999–2006 (2012).NEIMEF1053-8119 http://dx.doi.org/10.1016/j.neuroimage.2012.05.031 Google Scholar

27. M. Schweiger and S. Arridge, “The Toast++ software suite for forward and inverse modeling in optical tomography,” J. Biomed. Opt. 19, 040801 (2014).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.JBO.19.4.040801 Google Scholar

28. H. Dehghani et al., “Near infrared optical tomography using NIRFAST: algorithm for numerical model and image reconstruction,” Int. J. Numer. Methods Biomed. Eng. 25, 711–732 (2009).CANMER0748-8025 http://dx.doi.org/10.1002/cnm.1162 Google Scholar

29. Q. Fang, “Mesh-based Monte Carlo method using fast ray-tracing in Plücker coordinates,” Biomed. Opt. Express 1, 165–175 (2010).BOEICL2156-7085 http://dx.doi.org/10.1364/BOE.1.000165 Google Scholar

30. Q. Fang and D. A. Boas, “Monte Carlo simulation of photon migration in 3D turbid media accelerated by graphics processing units,” Opt. Express 17, 20178 (2009).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.17.020178 Google Scholar

31. A. Corlu et al., “Diffuse optical tomography with spectral constraints and wavelength optimization,” Appl. Opt. 44, 2082–2093 (2005).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.44.002082 Google Scholar

32. F. Larusson, S. Fantini and E. L. Miller, “Parametric level set reconstruction methods for hyperspectral diffuse optical tomography,” Biomed. Opt. Express 3, 1006–1024 (2012).BOEICL2156-7085 http://dx.doi.org/10.1364/BOE.3.001006 Google Scholar

33. B. W. Zeff et al., “Retinotopic mapping of adult human visual cortex with high-density diffuse optical tomography,” Proc. Natl. Acad. Sci. U.S.A. 104, 12169–12174 (2007). http://dx.doi.org/10.1073/pnas.0611266104 Google Scholar

34. H. R. Heekeren et al., “Noninvasive assessment of changes in cytochrome-C oxidase oxidation in human subjects during visual stimulation,” J. Cereb. Blood Flow Metab. 19, 592–603 (1999). http://dx.doi.org/10.1097/00004647-199906000-00002 Google Scholar

35. V. Fonov et al., “Unbiased average age-appropriate atlases for pediatric studies,” NeuroImage 54, 313–327 (2011).NEIMEF1053-8119 http://dx.doi.org/10.1016/j.neuroimage.2010.07.033 Google Scholar

36. S. Brigadoi and R. J. Cooper, “How short is short? Optimum source-detector distance for short-separation channels in functional near-infrared spectroscopy,” Neurophotonics 2, 025005 (2015). http://dx.doi.org/10.1117/1.NPh.2.2.025005 Google Scholar

37. L. A. Dempsey et al., “Data-driven approach to optimum wavelength selection for diffuse optical imaging,” J. Biomed. Opt. 20, 016003 (2015).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.JBO.20.1.016003 Google Scholar

38. Q. Fang and D. A. Boas, “Tetrahedral mesh generation from volumetric binary and grayscale images,” in IEEE Int. Symp. on Biomedical Imaging: From Nano to Macro, pp. 1142–1145, IEEE (2009). http://dx.doi.org/10.1109/ISBI.2009.5193259 Google Scholar

39. F. Bevilacqua et al., “In vivo local determination of tissue optical properties: applications to human brain,” Appl. Opt. 38, 4939–4950 (1999).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.38.004939 Google Scholar

40. A. Custo et al., “Effective scattering coefficient of the cerebral spinal fluid in adult head models for diffuse optical imaging,” Appl. Opt. 45, 4747–4755 (2006).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.45.004747 Google Scholar

41. G. Strangman, M. A. Franceschini and D. A. Boas, “Factors affecting the accuracy of near-infrared spectroscopy concentration calculations for focal changes in oxygenation parameters,” NeuroImage 18, 865–879 (2003).NEIMEF1053-8119 http://dx.doi.org/10.1016/S1053-8119(03)00021-1 Google Scholar

42. T. Correia, A. Gibson and J. Hebden, “Identification of the optimal wavelengths for optical topography: a photon measurement density function analysis,” J. Biomed. Opt. 15, 056002 (2010).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.3484747 Google Scholar

43. J. C. Horton and D. H. Hubel, “Regular patchy distribution of cytochrome oxidase staining in primary visual cortex of macaque monkey,” Nature 292, 762–764 (1981). http://dx.doi.org/10.1038/292762a0 Google Scholar

44. S. Powell, S. R. Arridge and T. S. Leung, “Gradient-based quantitative image reconstruction in ultrasound-modulated optical tomography: first harmonic measurement type in a linearised diffusion formulation,” IEEE Trans. Med. Imaging 35, 456–467 (2016).ITMID40278-0062 http://dx.doi.org/10.1109/TMI.2015.2478742 Google Scholar

45. C. C. Paige and M. A. Saunders, “LSQR—an algorithm for sparse linear-equations and sparse least-squares,” ACM Trans. Math. Software 8, 43–71 (1982).ACMSCU0098-3500 http://dx.doi.org/10.1145/355984.355989 Google Scholar

46. J. D. Swisher et al., “Visual topography of human intraparietal sulcus,” J. Neurosci. 27, 5326–5337 (2007).JNRSDS0270-6474 http://dx.doi.org/10.1523/JNEUROSCI.0991-07.2007 Google Scholar

47. S. D. Konecky et al., “Hyperspectral optical tomography of intrinsic signals in the rat cortex,” Neurophotonics 2, 045003 (2015). http://dx.doi.org/10.1117/1.NPh.2.4.045003 Google Scholar

48. T. Zhang et al., “Towards real-time detection of seizures in awake rats with GPU-accelerated diffuse optical tomography,” J. Neurosci. Methods 240, 28–36 (2015).JNMEDT0165-0270 http://dx.doi.org/10.1016/j.jneumeth.2014.10.018 Google Scholar

49. S. Prince et al., “Time-series estimation of biological factors in optical diffusion tomography,” Phys. Med. Biol. 48, 1491–1504 (2003).PHMBA70031-9155 http://dx.doi.org/10.1088/0031-9155/48/11/301 Google Scholar

50. S. G. Diamond et al., “Dynamic physiological modeling for functional diffuse optical tomography,” NeuroImage 30, 88–101 (2006).NEIMEF1053-8119 http://dx.doi.org/10.1016/j.neuroimage.2005.09.016 Google Scholar

51. C. E. Cooper and R. Springett, “Measurement of cytochrome oxidase and mitochondrial energetics by near-infrared spectroscopy,” Philos. Trans. R. Soc. B. Biol. Sci. 352, 669–676 (1997).PTRBAE0962-8436 http://dx.doi.org/10.1098/rstb.1997.0048 Google Scholar

52. M. Cope, “The application of near infrared spectroscopy to non-invasive monitoring of cerebral oxygenation in the newborn infant,” PhD Thesis, University of London (1991). Google Scholar

Biography

Sabrina Brigadoi studied bioengineering before receiving a PhD in cognitive science from the University of Padova, Italy. Formerly a research associate in the Biomedical Optics Research Laboratory at the University College London (UCL), she is now a postdoctoral fellow at the University of Padova and honorary research associate at the UCL. Her research interests are focused on advancing the applicability of diffuse optical techniques, with particular interest in signal processing techniques, image reconstruction, and head model development.

Phong Phan is a PhD student in the Biomedical Optics Research Laboratory at the UCL. With a medical background and an interest in engineering, he is focusing on research, using near-infrared spectroscopy and imaging, to investigate the spatial variation in cellular oxygen metabolism in the healthy and injured adult human brains.

Robert J. Cooper studied physics at New College, Oxford, before obtaining his PhD in medical physics from UCL. Formerly a postdoc at the Martinos Center for Biomedical Imaging at Massachusetts General Hospital, he is now an EPSRC early career fellow at the UCL. His research is focused on advancing and translating diffuse optical methods. He is a codirector of neoLAB: an interdisciplinary neonatal neuroimaging research group based in UCL and Cambridge.

Martin Smith is a consultant in the Neurocritical Care at the University College London Hospitals and honorary professor at the UCL. He has a clinical and research interest in monitoring and managing the acutely injured brain. He is a past president of the Society for Neuroscience in Anesthesiology and Critical Care, neuroscience section editor of Intensive Care Medicine, senior editor of Anesthesia and Analgesia, and reviews section editor of the Journal of Neurosurgical Anesthesiology.

Clare E. Elwell gained her BSc and MPhil degrees in medical physics from the University of Exeter, and her PhD from UCL. She is a professor of medical physics and director of the Near-Infrared Spectroscopy Research Group at UCL. She has studied autism and acute brain injury, and leads the Brain Imaging for Global Health project investigating brain development in Africa. She is a president of the Society for Functional Near-Infrared Spectroscopy.

Adam P. Gibson is a professor in medical physics at the UCL. His research includes diffuse optical imaging of breast and brain, radiotherapy, and imaging historical manuscripts.

Biographies for the other authors are not available.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation Download Citation

Sabrina Brigadoi, Sabrina Brigadoi, Phong Phan, Phong Phan, David Highton, David Highton, Samuel Powell, Samuel Powell, Robert J. Cooper, Robert J. Cooper, Jeremy C. Hebden, Jeremy C. Hebden, Martin Smith, Martin Smith, Ilias Tachtsidis, Ilias Tachtsidis, Clare E. Elwell, Clare E. Elwell, Adam P. Gibson, Adam P. Gibson, } "Image reconstruction of oxidized cerebral cytochrome C oxidase changes from broadband near-infrared spectroscopy data," Neurophotonics 4(2), 021105 (24 May 2017). https://doi.org/10.1117/1.NPh.4.2.021105 . Submission: Received: 3 March 2017; Accepted: 1 May 2017

Received: 3 March 2017; Accepted: 1 May 2017; Published: 24 May 2017

JOURNAL ARTICLE
10 PAGES

DOWNLOAD PDF + SAVE TO MY LIBRARY