A simple and novel algorithm for time-resolved multiwavelength oximetry

Approach to optimize 3-dimensional brain functional activation image with high resolution: a study on functional near-infrared spectroscopy

In this study, 3-dimensional (3-D) enhanced brain-function-map generation and estimation methodology is presented. Optical signals were modelled in the form of numerical optimization problem to infer the best existing waveform of canonical hemodynamic response function. Inter-channel activity patterns were also estimated. The estimation of activation of inter-channel gap depends on the minimization of generalized cross-validation. 3-D brain activation maps were produced through inverse discrete cosine transform. The proposed algorithm acquired significant results for 3-D functional maps with high resolution, in comparison with that of 2-D functional t-maps. A comprehensive analysis by exhibiting images corresponding to several layers has also been appended.

Near-Infrared Diffuse Optical Imaging for Early Prediction of Breast Cancer Response to Neoadjuvant Chemotherapy: A Comparative Study Using 18F-FDG PET/CT

Diffuse optical spectroscopic imaging (DOSI) is used as an indicator of tumor blood volume quantified by tissue hemoglobin concentrations. We aimed to determine whether early changes in tumor total hemoglobin (tHb) concentration can predict a pathologic complete response (pCR) to neoadjuvant chemotherapy in patients with operable breast cancer, and we compared the predictive value of pCR between DOSI and (18)F-FDG PET combined with CT.

METHODS

Of the 100 patients enrolled, 84 patients were prospectively evaluated for primary objective analysis. Sixty-four of the patients underwent both sequential DOSI scans at baseline after their first and second chemotherapy courses and (18)F-FDG PET/CT at baseline and after their second chemotherapy course. The mean tHb (tHbmean) concentration and SUVmax of the lesion were measured using DOSI and (18)F-FDG PET/CT, respectively, and the percentage change in tHbmean (∆tHbmean) and change in SUVmax (∆SUVmax) were calculated. We compared the diagnostic performances of DOSI and (18)F-FDG PET/CT for predicting pCR via the analysis of the receiver-operating-characteristic curves.

RESULTS

pCR was achieved in 16 patients, and neoadjuvant chemotherapy caused a significant reduction of ∆tHbmean in pCR compared with non-pCR after the 2 chemotherapy courses. When the tentative ∆tHbmean cutoff values after the first and second courses were used, the ability to predict pCR was as follows: 81.2% sensitivity/47.0% specificity and 93.7% sensitivity/47.7% specificity, respectively. Comparison of the diagnostic performances of DOSI and (18)F-FDG PET/CT revealed areas under the curve of 0.69 and 0.75 of ∆tHbmean after the first and second courses, respectively, which were lower than those of ∆SUVmax (0.90).

CONCLUSION

DOSI predicted pCR in patients with breast cancer with moderate accuracy. The diagnostic performance of DOSI was inferior to that of the early metabolic response as monitored by (18)F-FDG PET/CT.

Time-resolved diffuse optical spectroscopy: a differential absorption approach

A method is presented for the estimate of spectral changes in the absorption properties of turbid media from time-resolved diffuse optical spectroscopy. The method relies on the hypothesis of constant scattering over the wavelength range of interest, but no limitations come from the sample size and shape as the method is derived directly from the Beer-Lambert law. The effects of a moderate spectral dependence of the scattering properties and of the non-ideal instrument response function were investigated theoretically, and the results were confirmed experimentally, showing that the method can be profitably applied in cases of practical interest.

Estimation of the absorption coefficients of two-layered media by a simple method using spatially and time-resolved reflectances

Our newly developed method using spatially and time-resolved reflectances can easily estimate the absorption coefficients of each layer in a two-layered medium if the thickness of the upper layer and the reduced scattering coefficients of the two layers are known a priori. We experimentally validated this method using phantoms and examined its possibility of estimating the absorption coefficients of the tissues in human heads. In the case of a homogeneous plastic phantom (polyacetal block), the absorption coefficient estimated by our method agreed well with that obtained by a conventional method. Also, in the case of two-layered phantoms, our method successfully estimated the absorption coefficients of the two layers. Furthermore, the absorption coefficients of the extracerebral and cerebral tissue inside human foreheads were estimated under the assumption that the human heads were two-layered media. It was found that the absorption coefficients of the cerebral tissues were larger than those of the extracerebral tissues.

Progress of near-infrared spectroscopy and topography for brain and muscle clinical applications

This review celebrates the 30th anniversary of the first in vivo near-infrared (NIR) spectroscopy (NIRS) publication, which was authored by Professor Frans Jobsis. At first, NIRS was utilized to experimentally and clinically investigate cerebral oxygenation. Later it was applied to study muscle oxidative metabolism. Since 1993, the discovery that the functional activation of the human cerebral cortex can be explored by NIRS has added a new dimension to the research. To obtain simultaneous multiple and localized information, a further major step forward was achieved by introducing NIR imaging (NIRI) and tomography. This review reports on the progress of the NIRS and NIRI instrumentation for brain and muscle clinical applications 30 years after the discovery of in vivo NIRS. The review summarizes the measurable parameters in relation to the different techniques, the main characteristics of the prototypes under development, and the present commercially available NIRS and NIRI instrumentation. Moreover, it discusses strengths and limitations and gives an outlook into the "bright" future.

Clinical evaluation of time-resolved spectroscopy by measuring cerebral hemodynamics during cardiopulmonary bypass surgery

We developed a three-wavelength time-resolved spectroscopy (TRS) system, which allows quantitative measurement of hemodynamics within relatively large living tissue. We clinically evaluated this TRS system by monitoring cerebral circulation during cardiopulmonary bypass surgery. Oxyhemoglobin, deoxyhemoglobin, total hemoglobin and oxygen saturation (SO(2)) were determined by TRS on the left forehead attached with an optode spacing of 4 cm. We also simultaneously monitored jugular venous oxygen saturation (SjvO(2)) and arterial blood hematocrit (Hct) using conventional methods. The validity and usefulness of the TRS system were assessed by comparing parameters obtained with the TRS and conventional methods. Although the changes in SO(2) were lower than those in SjvO(2), SO(2) obtained by TRS paralleled the fluctuations in SjvO(2), and a good correlation between these values was observed. The only exceptions occurred during the perfusion period. Moreover, there was a good correlation between tHb and Hct values (r(2)=0.63). We concluded that time-resolved spectroscopy reflected the conditions of cerebral hemodynamics of patients during surgical operations.

Cerebral hemodynamics evaluation by near-infrared time-resolved spectroscopy: Correlation with simultaneous positron emission tomography measurements

We compared pharmacologically-perturbed hemodynamic parameters (cerebral blood volume; CBV, and flow; CBF) by acetazolamide administration in six healthy human subjects studied with positron emission tomography (PET) and near-infrared (NIR) time-resolved spectroscopy (TRS) simultaneously to investigate whether NIR-TRS could measure in vivo hemodynamics in the brain tissue quantitatively. Simultaneously with the PET measurements, TRS measurements were performed at the forehead with four different optode spacing from 2 cm to 5 cm. Total hemoglobin and oxygen saturation (SO2) measured by TRS significantly increased after administration of acetazolamide at any optode spacing in all subjects. In PET study, CBV and CBF were estimated in the following three volumes of interest (VOIs) determined on magnetic resonance images, VOI1: scalp and skull, VOI2: gray matter region, VOI3: gray and white matter regions. Acetazolamide treatment elevated CBF and CBV significantly in VOI2 and VOI3 but VOI1. TRS-derived CBV was more strongly correlated with PET-derived counterpart in VOI2 and VOI3 when the optode spacing was above 4 cm, although optical signal from cerebral tissue could be caught with any optode spacing. As to increase of the CBV, 4 cm of optode spacing correlated best with VOI2. To support the result of TRS-PET experiment, we also estimated the contribution ratios of intracerebral tissue to observed absorption change based on diffusion theory. The contribution ratios at 4 cm were estimated as follows: 761 nm: 50%, 791 nm: 72%, 836 nm: 70%. These results demonstrated that NIR-TRS with 4 cm of optode spacing could measure cerebral hemodynamic responses optimally and quantitatively.

A tentative proposal on physiological polymorphism and its experimental approaches

The purpose of this paper is to propose a tentative concept of physiological polymorphism and experimental approaches to it. We think that the concept of physiological polymorphism is to categorize the diversified phenomena into various types with statistical methods and to explain the differences among the categorized types from the viewpoint of their physiological mechanisms. Furthermore, it is necessary to take into consideration the fact that physiological polymorphism is observed as a phenotype, and the phenotype results from a genotype modified by culture and environment. As an experimental approach, we studied the effects of gustatory stimulation by chocolate on the activities of the prefrontal area and found that the activities were increased in some cases and decreased in other cases. Therefore, to begin with, when we divided them into an "increasing group" and a "decreasing group," we found that the increasing group had many subjects of Type B and High anxiety, and the decreasing group had many subjects of Type A and Normal anxiety. By the chi-square test for independence, it was found that the ratio of "increase" and "decrease" was related to the trait anxiety and type A personality, respectively. Next, we divided the activities of the prefrontal area into Type A and Type B, as well as a high anxiety group and a low anxiety group. As a result, the Type B and high anxiety groups showed significantly increased activities, while the Type A and normal anxiety groups showed no changes in their activities. Consequently, this data enabled us to explain the difference in activities of the prefrontal area from the viewpoint of personality characteristics. To conclude, we were able to categorize diversified changes in the prefrontal area into certain types when a gustatory stimulus was applied, and to explain them by using personality characteristics (State Trait Anxiety Inventory, Type A behavioral pattern) that are commonly known for their reflection of genotypes.

In vivo absorption spectroscopy of tumor sensitizers with femtosecond white light

A system based on a femtosecond white-light continuum and a streak camera was used for recordings of the in vivo absorption spectra of the tumor-seeking agent disulphonated aluminum phthalocyanine. Measurements for different drug doses were performed on tumor tissue (muscle-implanted adenocarcinoma) and normal muscle tissue in rats. It was found that the shape of the spectrum is tissue dependent. The peak of the absorption spectrum is blueshifted in tumor tissue as compared with the muscle. Thus the contrast in the drug-related absorption can be altered by up to a factor of 2 from the primary drug molecular-concentration contrast between normal muscle and tumor by the proper selection of the illumination wavelength.

Artefacts in the analysis of temporal response functions measured by photon counting

The least-squares (LS) method in fluorescence decay analyses and in time-domain analyses of the diffuse scattering light for data measured by the time-correlated single photon counting (TCSPC) technique is experimentally evaluated, and the artefact in LS analysis for data with different counting statistics is discussed. In single exponential decay analysis, the error of the decay parameter by the LS method is smaller than 10% of the expected true value when the average number of counts per bin (N/k) is more than 1, and the fitting region covers a period on the order of the decay time. In multi-exponential analysis, the decay parameters are sensitively dependent on the counting statistics. In contrast, the fitting by the maximum likelihood estimation (MLE), assuming Poissonian statistics, greatly reduces such dependence of parameters on the counting statistics. In another application, time-domain diffuse scattering measurements, the LS method is only accurate at N/k > 50 (10% error in the absorption coefficient). In particular, the absorption coefficient is largely dependent on the count. In both examples, the problem of stability in the fitting process by MLE still remains: the convergence of the fitting is critically dependent on the selection of initial guesses of the parameters in contrast to the convergence in the LS method. Thus, a hybrid method using the LS method for the determination of the initial guesses is a practical solution to this problem.

Simple peak shift analysis of time-of-flight data with a slow instrumental response function

Analysis of time-of-flight (TOF) data is sometimes limited by the instrumental response function, and optical parameters are extracted from the observed response curve by several mathematical methods, such as deconvolution. In contrast to this, we demonstrate that a method using shifts of the peak time of the response curve with different source-detector separations can yield the average path length of the light traveling in a tissue-like sample without deconvolution. In addition, combining the intensity information allows us to separate the scattering and absorption coefficients. This simple method is more robust in signal-to-noise ratio than the moment analysis, which also does not require the deconvolution procedure, because the peak position is not significantly dependent on the baseline fluctuation and the contamination of the scattering. The analysis is demonstrated by TOF measurements of an Intralipid solution at 800 nm, and is applied to the measurements at 1.29 microm, where the temporal response of photomultiplier tubes is not sufficiently good.

Photon path distribution and optical responses of turbid media: theoretical analysis based on the microscopic Beer-Lambert law

A concise theoretical treatment has been developed to describe the optical responses of a highly scattering inhomogeneous medium using functions of the photon path distribution (PPD). The treatment is based on the microscopic Beer-Lambert law and has been found to yield a complete set of optical responses by time- and frequency-domain measurements. The PPD is defined for possible photons having a total zigzag pathlength of l between the points of light input and detection. Such a distribution is independent of the absorption properties of the medium and can be uniquely determined for the medium under quantification. Therefore, the PPD can be calculated with an imaginary reference medium having the same optical properties as the medium under quantification except for the absence of absorption. One of the advantages of this method is that the optical responses, the total attenuation, the mean pathlength, etc are expressed by functions of the PPD and the absorption distribution.

Visuospatial imagery is a fruitful strategy for the digit span backward task: a study with near-infrared optical tomography

Our newly developed 64-channel time-resolved optical tomographic imaging system using near-infrared light enables us to obtain a quantitative image of hemoglobin concentration changes associated with neuronal activation in the human brain ¿H. Eda, I. Oda, Y. Ito, Y. Wada, Y. Oikawa, Y. Tsunazawa, M. Takada, Y. Tsuchiya, Y. Yamashita, M. Oda, A. Sassaroll, Y. Yamada, M. Tamura, Multi-channel time-resolved optical tomographic imaging system, Rev. Sci. Instrum., 70 (1999) 3595-3602. Here, we used this optical imaging system to demonstrate that the backward digit span (DB) task activated the dorsolateral prefrontal cortex (DLPFC) of each hemisphere more than the forward digit span (DF) task in healthy adult volunteers, and higher performance of the DB task was closely related to the activation of the right DLPFC. These results suggest that visuospatial imagery is a useful strategy for the DB task. Optical tomography described here is a new modality of neuropsychological studies.

Compact tissue oximeter based on dual-wavelength multichannel time-resolved reflectance

We developed a compact dual-wavelength multichannel tissue oximeter based on the time-correlated single-photon counting (TCSPC) technique. The light sources are two pulsed diode lasers (output wavelengths of 672 and 818 nm, an average power of 1 mW, a pulse duration of 100 ps, and a pulse-repetition rate as high as 80 MHz). The time-resolved reflectance photons are detected by a multianode photomultiplier, and the output signals are redirected by a router to different memory blocks of the TCSPC personal computer board. The system's accuracy in determining the absorption microa and the reduced-scattering micros' coefficients and in reconstructing absorber concentrations in diffusive media was tested on phantoms. Preliminary in vivo tissue-oxygenation measurements were performed on healthy volunteers under different physiological conditions with a minimum acquisition time of 100 ms and an injected power of less than 100 microW.

Optical quantitation of absorbers in variously shaped turbid media based on the microscopic Beer-Lambert law. A new approach to optical computerized tomography

To determine the concentrations of an absorber in variously shaped turbid media such as human tissue, we propose analytical expressions for diffuse re-emission in time and frequency domains, based on the microscopic Beer-Lambert law that holds true when we trace a zigzag photon path in the medium. Our expressions are implicit for the scattering properties, the volume shape, and the source-detector separation. We show that three observables are sufficient to determine the changes in the concentration and the absolute concentrations of an absorber in scattering media as long as the scattering property remains constant. The three observables are: the re-emission, the mean pathlength or group delay, and the extinction coefficient of the absorber. We also show that our equations can be extended to describe photon migration in nonuniform media. The validity of the predictions is confirmed by measuring a tissue-like phantom.

Redox behavior of cytochrome oxidase in the rat brain measured by near-infrared spectroscopy

Using near-infrared spectroscopy, we developed a new approach for measuring the redox state of cytochrome oxidase in the brain under normal blood-circulation conditions. Our algorithm does not require the absorption coefficient of cytochrome oxidase, which differs from study to study. We employed this method for evaluation of effects of changes in oxygen delivery on cerebral oxygenation in rats. When fractional inspired oxygen was decreased in a stepwise manner from 100 to <10%, at which point the concentration of oxygenated hemoglobin ([HbO2]) decreased by approximately 60%, cytochrome oxidase started to be reduced. Increases in arterial PO2 under hyperoxic conditions caused an increase in [HbO2], whereas further oxidation of cytochrome oxidase was not observed. The dissociation of the responses of hemogloblin and cytochrome oxidase was also clearly observed after the injection of epinephrine under severely hypoxic conditions; that is, cytochrome oxidase was reoxidized with increasing blood pressure, whereas hemoglobin oxygenation was not changed. These data indicated that oxygen-dependent redox changes in cytochrome oxidase occur only when oxygen delivery is extremely impaired. This is consistent with the in vitro data of our previous study.

Relationship between time-resolved and non-time-resolved Beer-Lambert law in turbid media

The time-resolved Beer-Lambert law proposed for oxygen monitoring using pulsed light was extended to the non-time-resolved case in a scattered medium such as living tissues with continuous illumination. The time-resolved Beer-Lambert law was valid for the phantom model and living tissues in the visible and near-infrared regions. The absolute concentration and oxygen saturation of haemoglobin in rat brain and thigh muscle could be determined. The temporal profile of rat brain was reproduced by Monte Carlo simulation. When the temporal profiles of rat brain under different oxygenation states were integrated with time, the absorbance difference was linearly related to changes in the absorption coefficient. When the simulated profiles were integrated, there was a linear relationship within the absorption coefficient which was predicted for fractional inspiratory oxygen concentration from 10 to 100% and, in the case beyond the range of the absorption coefficient, the deviation from linearity was slight. We concluded that an optical pathlength which is independent of changes in the absorption coefficient is a good approximation for near-infrared oxygen monitoring.