Time resolved reflectance and transmittance for the non-invasive measurement of tissue optical properties

Theoretical optimal modulation frequencies for scattering parameter estimation and ballistic photon filtering in diffusing media

The efficiency of using intensity modulated light for the estimation of scattering properties of a turbid medium and for ballistic photon discrimination is theoretically quantified in this article. Using the diffusion model for modulated photon transport and considering a noisy quadrature demodulation scheme, the minimum-variance bounds on estimation of parameters of interest are analytically derived and analyzed. The existence of a variance-minimizing optimal modulation frequency is shown and its evolution with the properties of the intervening medium is derived and studied. Furthermore, a metric is defined to quantify the efficiency of ballistic photon filtering which may be sought when imaging through turbid media. The analytical derivation of this metric shows that the minimum modulation frequency required to attain significant ballistic discrimination depends only on the reduced scattering coefficient of the medium in a linear fashion for a highly scattering medium.

Generalized Beer-Lambert model for near-infrared light propagation in thick biological tissues

The attenuation of near-infrared (NIR) light intensity as it propagates in a turbid medium like biological tissue is described by modified the Beer–Lambert law (MBLL). The MBLL is generally used to quantify the changes in tissue chromophore concentrations for NIR spectroscopic data analysis. Even though MBLL is effective in terms of providing qualitative comparison, it suffers from its applicability across tissue types and tissue dimensions. In this work, we introduce Lambert-W function-based modeling for light propagation in biological tissues, which is a generalized version of the Beer–Lambert model. The proposed modeling provides parametrization of tissue properties, which includes two attenuation coefficients μ0 and η. We validated our model against the Monte Carlo simulation, which is the gold standard for modeling NIR light propagation in biological tissue. We included numerous human and animal tissues to validate the proposed empirical model, including an inhomogeneous adult human head model. The proposed model, which has a closed form (analytical), is first of its kind in providing accurate modeling of NIR light propagation in biological tissues.

NIR Spectroscopic Detection of Breast Cancer

Near infrared spectroscopy (NIRS) utilizes intrinsic optical absorption signals of blood, water, and lipid concentration available in the NIR window (600–1000 nm) as well as a developing array of extrinsic organic compounds to detect and localize cancer. This paper reviews optical cancer detection made possible through high tumor-tissue signal-to-noise ratio (SNR) and providing biochemical and physiological data in addition to those obtained via other methods.

NIRS detects cancers in vivo through a combination of blood volume and oxygenation from measurements of oxy- and deoxy-hemoglobin giving signals of tumor angiogenesis and hypermetabolism. The Chance lab tends towards CW breast cancer systems using manually scannable detectors with calibrated low pressure tissue contact. These systems calculate angiogenesis and hypermetabolism by using a pair of wavelengths and referencing the mirror image position of the contralateral breast to achieve high ROC/AUC. Time domain and frequency domain spectroscopy were also used to study similar intrinsic breast tumor characteristics such as high blood volume. Other NIRS metrics are water-fat ratio and the optical scattering coefficient. An extrinsic FDA approved dye, ICG, has been used to measure blood pooling with extravasation, similar to Gadolinium in MRI.

A key future development in NIRS will be new Molecular Beacons targeting cancers and fluorescing in the NIR window to enhance in vivo tumor-tissue ratios and to afford biochemical specificity with the potential for effective photodynamic anti-cancer therapies.

Spatial and Spectral Information in Optical Mammography

This article reviews our research activities in the area of optical mammography and relates them to the historical developments and the current state and trends in the field. The guiding threads for this article are the roles played in optical mammography by spatial and spectral information. The first feature, spatial information, is limited by the diffusive nature of light propagation but can take advantage of the exceptionally high optical contrast featured by blood vessels and blood-rich areas in the breast. We describe a method to correct for edge effects, a spatial second-derivative algorithm, and a two-dimensional phased-array approach that enhance the image contrast, the spatial resolution, and the depth discrimination in optical mammograms. The second feature, spectral information, is the most powerful and unique capability of optical mammography, and allows for functional measurements associated with hemoglobin concentration and oxygenation, water concentration, lipids content, and the wavelength dependence of tissue scattering. We present oxygenation-index images obtained from multi-wavelength optical data that point to the diagnostic potential of oxygenation information in optical mammography. The optimization of the spatial and spectral information in optical mammography has the potential to create a role for this imaging modality in the detection and monitoring of breast cancer.

T2-prepared velocity selective labelling: A novel idea for full-brain mapping of oxygen saturation

BACKGROUND AND AIM

Disturbances in cerebral oxygenation saturation (SO₂) have been linked to adverse outcome in adults, children, and neonates. In intensive care, the cerebral SO₂ is increasingly being monitored by Near-InfraRed Spectroscopy (NIRS). Unfortunately NIRS has a limited penetration depth. The "modified T₂-prepared Blood Imaging of Oxygen Saturation" (T₂-BIOS) MR sequence provides a step towards full brain SO₂ measurement.

MATERIALS AND METHODS

Tissue SO₂, and venous SO₂ (S_vO₂) were obtained simultaneously by T₂-BIOS during a respiratory challenge in ten healthy volunteers. These two measures were compared to SO₂ that was obtained by a single probe MR-compatible NIRS setup, and to cerebral blood flow and venous SO₂ that were obtained by arterial spin labelling and T₂-TRIR, respectively.

RESULTS

SO_2-T2-BIOS and SO_2-NIRS had a mean bias of -4.0% (95% CI -21.3% to 13.3%). S_vO_2-T2-BIOS correlated with SO_2-NIRS (R²=0.41, p=0.002) and S_vO_2-T2-TRIR (R²=0.87, p=0.002). In addition, SO_2-NIRS correlated with S_vO_2-T2-TRIR (R²=0.85, p=0.003) Frontal cerebral blood flow correlated with SO_2-T2-BIOS (R²=0.21, p=0.04), but was not significant in relation to SO_2-NIRS.

DISCUSSION/CONCLUSION

Full brain SO₂ assessment by any technique may help validating NIRS and may prove useful in guiding the clinical management of patient populations with cerebral injury following hypoxic-ischaemic events. The agreement between NIRS and T₂-BIOS provides confidence in measuring cerebral SO₂ by either technique. As it stands now, the T₂-BIOS represents a novel idea and future work will focus on improvements to make it a reliable tool for SO₂ assessment.

Mathematical model to interpret localized reflectance spectra measured in the presence of a strong fluorescence marker

Quantification of multiple fluorescence markers during neurosurgery has the potential to provide complementary contrast mechanisms between normal and malignant tissues, and one potential combination involves fluorescein sodium (FS) and aminolevulinic acid-induced protoporphyrin IX (PpIX). We focus on the interpretation of reflectance spectra containing contributions from elastically scattered (reflected) photons as well as fluorescence emissions from a strong fluorophore (i.e., FS). A model-based approach to extract μa and μ′s in the presence of FS emission is validated in optical phantoms constructed with Intralipid (1% to 2% lipid) and whole blood (1% to 3% volume fraction), over a wide range of FS concentrations (0 to 1000  μg/ml 1000  μg/ml ). The results show that modeling reflectance as a combination of elastically scattered light and attenuation-corrected FS-based emission yielded more accurate tissue parameter estimates when compared with a nonmodified reflectance model, with reduced maximum errors for blood volume (22% versus 90%), microvascular saturation (21% versus 100%), and μs′ (13% versus 207%). Additionally, quantitative PpIX fluorescence sampled in the same phantom as FS showed significant differences depending on the reflectance model used to estimate optical properties (i.e., maximum error 29% versus 86%). These data represent a first step toward using quantitative optical spectroscopy to guide surgeries through simultaneous assessment of FS and PpIX.

Reproducibility of parameters of postocclusive reactive hyperemia measured by diffuse optical tomography

The application of near-infrared spectroscopy (NIRS) to assess microvascular function has shown promising results. An important limitation when using a single source-detector pair, however, is the lack of depth sensitivity. Diffuse optical tomography (DOT) overcomes this limitation using an array of sources and detectors that allow the reconstruction of volumetric hemodynamic changes. This study compares the key parameters of postocclusive reactive hyperemia measured in the forearm using standard NIRS and DOT. We show that while the mean parameter values are similar for the two techniques, DOT achieves much better reproducibility, as measured by the intraclass correlation coefficient (ICC). We show that DOT achieves high reproducibility for muscle oxygen consumption (ICC: 0.99), time to maximal HbO2 (ICC: 0.94), maximal HbO2 (ICC: 0.99), and time to maximal HbT (ICC: 0.99). Absolute reproducibility as measured by the standard error of measurement is consistently smaller and close to zero (ideal value) across all parameters measured by DOT compared to NIRS. We conclude that DOT provides a more robust characterization of the reactive hyperemic response and show how the availability of volumetric hemodynamic changes allows the identification of areas of temporal consistency, which could help characterize more precisely the microvasculature.

Optical properties of drying wood studied by time-resolved near-infrared spectroscopy

We measured the optical properties of drying wood with the moisture contents ranging from 10% to 200%. By using time-resolved near-infrared spectroscopy, the reduced scattering coefficient μ_s' and absorption coefficient μ_a were determined independent of each other, providing information on the chemical and structural changes, respectively, of wood on the nanometer scale. Scattering from dry pores dominated, which allowed us to determine the drying process of large pores during the period of constant drying rate, and the drying process of smaller pores during the period of decreasing drying rate. The surface layer and interior of the wood exhibit different moisture states, which affect the scattering properties of the wood.

The role of perfusion in the oxygen extraction capability of skin and skeletal muscle

Oxygen extraction (OE) by all cells is dependent on an adequate supply of oxygen in proximal blood vessels and the cell's need and ability to uptake that oxygen. Here the role of blood flow in regulating OE in skin and skeletal muscle was investigated in lean and obese men. OE was derived by two optical reflectance spectroscopy techniques: 1) from the rate of fall in mean blood saturation during a 4 min below knee arterial occlusion, and thus no blood flow, in calf skin and skeletal muscle and 2) in perfused, unperturbed skin, using the spontaneous falls in mean blood saturation induced by vasomotion in calf and forearm skin of 24 subjects, 12 lean and 12 obese. OE in perfused skin was significantly higher in lean compared with obese subjects in forearm (Mann-Whitney, P < 0.004) and calf (P < 0.001) and did not correlate with OE in unperfused skin (ρ = -0.01, P = 0.48). With arterial occlusion and thus no blood flow, skin OE in lean and obese subjects no longer differed (P = 0.23, not significant). In contrast in skeletal muscle with arterial occlusion and no blood flow, the difference in OE between lean and obese subjects occurred, with obese subjects exhibiting significantly higher OE (P < 0.012). The classic model of metabolic blood flow regulation to support oxygen extraction is evident in perfused skin; OE is perturbed without blood flow and reduced in obesity. In resting skeletal muscle other mechanism(s), independent of blood flow, are implicated in oxygen extraction.

Time-resolved subtraction method for measuring optical properties of turbid media

Near-infrared spectroscopy is a noninvasive optical method used primarily to monitor tissue oxygenation due to the absorption properties of hemoglobin. Accurate estimation of hemoglobin concentrations and other light absorbers requires techniques that can separate the effect of absorption from the much greater effect of light scattering. One of the most advanced methods is time-resolved near-infrared spectroscopy (TR-NIRS), which measures the absorption and scattering coefficients of a turbid medium by modeling the recorded distribution time of flight of photons. A challenge with TR-NIRS is that it requires accurate characterization of the dispersion caused by the system. In this study, we present a method for circumventing this problem by applying statistical moment analysis to two time-of-flight distributions measured at separated source-detector distances. Simulations based on analytical models and Monte Carlo code, and tissue-mimicking phantoms, were used to demonstrate its accuracy for source-detector distances typically used in neuroimaging applications. The simplicity of the approach is well suited to real-time applications requiring accurate quantification of the optical properties of a turbid medium.

Finite-size effect in light transmission through highly forward scattering media at grazing angles

We present a theoretical study of light transmission through a disordered medium with large (compared to the light wavelength) inhomogeneities. Both numerical integration and analytic treatments of the radiative transfer equation are performed. An effect of the single-scattering phase function on the total transmittance is found in a subdiffusion thickness range. The effect reveals itself at grazing angles of incidence and originates from small-angle multiple scattering of light. A simple analytic formula for the total transmittance is derived. Our results are in good agreement with data of independent numerical calculations.

Determination of true optical absorption and scattering coefficient of wooden cell wall substance by time-of-flight near infrared spectroscopy

The true absorption coefficient (μa) and reduced scattering coefficient (μ´s) of the cell wall substance in Douglas fir were determined using time-of-flight near infrared spectroscopy. Samples were saturated with hexane, toluene or quinolone to minimize the multiple reflections of light on the boundary between pore-cell wall substance in wood. μ´s exhibited its minimum value when the wood was saturated with toluene because the refractive index of toluene is close to that of the wood cell wall substance. The optical parameters of the wood cell wall substance calculated were μa = 0.030 mm(-1) and μ´s= 18.4 mm(-1). Monte Carlo simulations using these values were in good agreement with the measured time-resolved transmittance profiles.

Effect of blood flow restriction on tissue oxygenation during knee extension

PURPOSE

Time-resolved near-infrared spectroscopy was used to quantify tissue oxy- and deoxyhemoglobin concentrations ([HbO2] and [HbR]) and O2 saturation (stO2) in the oblique fibers of the vastus medialis muscle and brain prefrontal cortex during knee extension with and without blood flow restriction (BFR).

METHODS

Six young healthy males performed three sets of knee extensions on a dynamometer (50% one-repetition maximum) separated by 90-s rest periods in three conditions: 1) until fatigue without BFR (fatigue), 2) until fatigue with BFR (100 mm Hg cuff constriction around thigh (BFR)), 3) same number of repetitions from condition 2 without BFR (matched). Each condition was performed on a separate visit.

RESULTS

BFR was associated with higher [HbR] at the oblique fibers of the vastus medialis muscle (rest 1: 57.8 (BFR) vs 35.0 μM (matched); P < 0.0001) and a significantly lower stO2 during recovery periods between sets (7.5%-11.2 % lower than non-BFR conditions for rest 1 and 2, P < 0.0001). Using a piecewise linear spline method, a spike in [HbR] was observed before the onset of HbR clearance during recovery, causing HbR clearance to begin at a higher concentration (81 (BFR) vs 62 μM (matched), P = 0.029). [HbO2] kinetics during recovery were also affected by BFR, with longer duration (BFR, 51 s; matched, 31 s; P = 0.047) but lower rate of increase (BFR, 58 μM·min; matched, 89 μM·min; P = 0.004) during recovery. In the prefrontal cortex, BFR was associated with increased [HbR], diminished increase in [HbO2], and higher subjective exertion.

CONCLUSIONS

These findings yield insight into possible physiological mechanisms of BFR and suggest a role of time-resolved near-infrared spectroscopy in monitoring and optimization of BFR exercise on an individual basis.

Effect of a thin superficial layer on the estimate of hemodynamic changes in a two-layer medium by time domain NIRS

In order to study hemodynamic changes involved in muscular metabolism by means of time domain fNIRS, we need to discriminate in the measured signal contributions coming from different depths. Muscles are, in fact, typically located under other tissues, e.g. skin and fat. In this paper, we study the possibility to exploit a previously proposed method for analyzing time-resolved fNIRS measurements in a two-layer structure with a thin superficial layer. This method is based on the calculation of the time-dependent mean partial pathlengths. We validated it by simulating venous and arterial arm cuff occlusions and then applied it on in vivo measurements.

Toward noninvasive assessment of flap viability with time-resolved diffuse optical tomography: a preclinical test on rats

The noninvasive assessment of flap viability in autologous reconstruction surgery is still an unmet clinical need. To cope with this problem, we developed a proof-of-principle fully automatized setup for fast time-gated diffuse optical tomography exploiting Mellin-Laplace transform to obtain three-dimensional tomographic reconstructions of oxy- and deoxy-hemoglobin concentrations. We applied this method to perform preclinical tests on rats inducing total venous occlusion in the cutaneous abdominal flaps. Notwithstanding the use of just four source-detector couples, we could detect a spatially localized increase of deoxyhemoglobin following the occlusion (up to 550 μM in 54 min). Such capability to image spatio-temporal evolution of blood perfusion is a key issue for the noninvasive monitoring of flap viability.

Diffuse Optical Characterization of the Healthy Human Thyroid Tissue and Two Pathological Case Studies

The in vivo optical and hemodynamic properties of the healthy (n = 22) and pathological (n = 2) human thyroid tissue were measured non-invasively using a custom time-resolved spectroscopy (TRS) and diffuse correlation spectroscopy (DCS) system. Medical ultrasound was used to guide the placement of the hand-held hybrid optical probe. TRS measured the absorption and reduced scattering coefficients (μ_a, μ_s′) at three wavelengths (690, 785 and 830 nm) to derive total hemoglobin concentration (THC) and oxygen saturation (StO₂). DCS measured the microvascular blood flow index (BFI). Their dependencies on physiological and clinical parameters and positions along the thyroid were investigated and compared to the surrounding sternocleidomastoid muscle. The THC in the thyroid ranged from 131.9 μM to 144.8 μM, showing a 25–44% increase compared to the surrounding sternocleidomastoid muscle tissue. The blood flow was significantly higher in the thyroid (BFI_thyroid = 16.0 × 10^-9 cm²/s) compared to the muscle (BFI_muscle = 7.8 × 10^-9 cm²/s), while StO₂ showed a small (StO_{2, muscle} = 63.8% to StO_{2, thyroid} = 68.4%), yet significant difference. Two case studies with thyroid nodules underwent the same measurement protocol prior to thyroidectomy. Their THC and BFI reached values around 226.5 μM and 62.8 × 10^-9 cm²/s respectively showing a clear contrast to the nodule-free thyroid tissue as well as the general population. The initial characterization of the healthy and pathologic human thyroid tissue lays the ground work for the future investigation on the use of diffuse optics in thyroid cancer screening.

Hemodynamic Measurements of the Human Adult Head in Transmittance Mode by Near-Infrared Time-Resolved Spectroscopy

Using a near-infrared time-resolved spectroscopy (TRS) system, we measured the human head in transmittance mode to obtain the optical properties and the hemodynamic changes of deep brain tissues in seven healthy adult volunteers during hyperventilation. For six out of seven volunteers, we obtained the optical signals with sufficient intensity within 10 sec. of sampling. We confirmed that it is possible to non-invasively measure the hemodynamic changes of the human head during hyperventilation, even in the transmittance measurements by the developed TRS system. These results showed that the level of deoxygenated hemoglobin was significantly increased, and the level of oxygenated and total hemoglobin and tissue oxygen saturation were also significantly decreased during hyperventilation. We expect that this TRS technique will be applied to clinical applications for measuring deep brain tissues and deep biological organs.

Interferometric Near-Infrared Spectroscopy (iNIRS) for determination of optical and dynamical properties of turbid media

We introduce and implement interferometric near-infrared spectroscopy (iNIRS), which simultaneously extracts optical and dynamical properties of turbid media through analysis of a spectral interference fringe pattern. The spectral interference fringe pattern is measured using a Mach-Zehnder interferometer with a frequency-swept narrow linewidth laser. Fourier analysis of the detected signal is used to determine time-of-flight (TOF)-resolved intensity, which is then analyzed over time to yield TOF-resolved intensity autocorrelations. This approach enables quantification of optical properties, which is not possible in conventional, continuous-wave near-infrared spectroscopy (NIRS). Furthermore, iNIRS quantifies scatterer motion based on TOF-resolved autocorrelations, which is a feature inaccessible by well-established diffuse correlation spectroscopy (DCS) techniques. We prove this by determining TOF-resolved intensity and temporal autocorrelations for light transmitted through diffusive fluid phantoms with optical thicknesses of up to 55 reduced mean free paths (approximately 120 scattering events). The TOF-resolved intensity is used to determine optical properties with time-resolved diffusion theory, while the TOF-resolved intensity autocorrelations are used to determine dynamics with diffusing wave spectroscopy. iNIRS advances the capabilities of diffuse optical methods and is suitable for in vivo tissue characterization. Moreover, iNIRS combines NIRS and DCS capabilities into a single modality.

Fabrication and application of heterogeneous printed mouse phantoms for whole animal optical imaging

This work demonstrates the usefulness of 3D printing for optical imaging applications. Progress in developing optical imaging for biomedical applications requires customizable and often complex objects for testing and evaluation. There is therefore high demand for what have become known as tissue-simulating "phantoms." We present a new optical phantom fabricated using inexpensive 3D printing methods with multiple materials, allowing for the placement of complex inhomogeneities in complex or anatomically realistic geometries, as opposed to previous phantoms, which were limited to simple shapes formed by molds or machining. We use diffuse optical imaging to reconstruct optical parameters in 3D space within a printed mouse to show the applicability of the phantoms for developing whole animal optical imaging methods. This phantom fabrication approach is versatile, can be applied to optical imaging methods besides diffusive imaging, and can be used in the calibration of live animal imaging data.

Evaluation of venous return in lower limb by passive ankle exercise performed by PHARAD

This paper presents evaluation of venous return, i.e., blood flow volume of vein (BF), in the lower limb after passive exercise performed by our developed "parallel link type human ankle rehabilitation assistive device (PHARAD)". The PHARAD can perform complex passive exercises (plantar flexion/dorsiflexion, inversion/eversion, adduction/abduction, and combination of these motions) by reproducing input motions of a foot plate that is attached to a sole of foot. The passive exercise can be performed for not only rehabilitation but also prevention of deep vein thrombosis (DVT). In this study, we measured the concentration of Total hemoglobin (Total-Hb) using multi-channel near infra-red spectroscopy (NIRS)-based tissue oximeters and calculated a gradient of Total-Hb during a venous occlusion. We defined the gradient as BF and evaluated BF after 3 min passive exercise performed by the PHARAD comparing to BF of resting. Seven healthy young adult people were recruited for the experiment and we assessed passive exercise, active exercise, and walking. Experimental results show that BF after the passive exercises significantly increases compare to BF of resting and this indicates that passive exercises performed by the PHARAD increases BF and has a potential to prevent DVT.

Monte Carlo Modeling of Photon Propagation Reveals Highly Scattering Coral Tissue

Corals are very efficient at using solar radiation, with photosynthetic quantum efficiencies approaching theoretical limits. Here, we investigated potential mechanisms underlying such outstanding photosynthetic performance through extracting inherent optical properties of the living coral tissue and skeleton in a massive faviid coral. Using Monte Carlo simulations developed for medical tissue optics it is shown that for the investigated faviid coral, the coral tissue was a strongly light scattering matrix with a reduced scattering coefficient of μ_s' = 10 cm^-1 (at 636 nm). In contrast, the scattering coefficient of the coral skeleton was μ_s' = 3.4 cm^-1, which facilitated the efficient propagation of light to otherwise shaded coral tissue layers, thus supporting photosynthesis in lower tissues. Our study provides a quantification of coral tissue optical properties in a massive faviid coral and suggests a novel light harvesting strategy, where tissue and skeletal optics act in concert to optimize the illumination of the photosynthesizing algal symbionts embedded within the living coral tissue.

Finite element modeling of light propagation in turbid media under illumination of a continuous-wave beam

Spatially resolved spectroscopy provides a means for measuring the optical properties of biological tissues, based on analytical solutions to diffusion approximation for semi-infinite media under the normal illumination of an infinitely small light beam. The method is, however, prone to error in measurement because the actual boundary condition and light beam often deviate from that used in deriving the analytical solutions. It is therefore important to quantify the effect of different boundary conditions and light beams on spatially resolved diffuse reflectance in order to improve the measurement accuracy of the technique. This research was aimed at using finite element method (FEM) to model light propagation in turbid media, subjected to normal illumination by a continuous-wave beam of infinitely small or finite size. Three types of boundary conditions [i.e., partial current (PCBC), extrapolated (EBC), and zero (ZBC)] were evaluated and compared against Monte Carlo (MC) simulations, since MC could provide accurate fluence rate and diffuse reflectance. The effect of beam size was also investigated. Overall results showed that FEM provided results as accurate as those of the analytical method when an appropriate boundary condition was applied. ZBC did not give satisfactory results in most cases. FEM-PCBC yielded a better fluence rate at the boundary than did FEM-EBC, while they were almost identical in predicting diffuse reflectance. Results further showed that FEM coupled with EBC effectively simulated spatially resolved diffuse reflectance under the illumination of a finite size beam. A large beam introduced more error, especially within the region of illumination. Research also confirmed an earlier finding that a light beam of less than 1 mm diameter should be used for estimation of optical parameters. FEM is effective for modeling light propagation in biological tissues and can be used for improving the optical property measurement by the spatially resolved technique.

Effect of the chest wall on the measurement of hemoglobin concentrations by near-infrared time-resolved spectroscopy in normal breast and cancer

BACKGROUND

Optical imaging and spectroscopy using near-infrared light have great potential in the assessment of tumor vasculature. We previously measured hemoglobin concentrations in breast cancer using a near-infrared time-resolved spectroscopy system. The purpose of the present study was to evaluate the effect of the chest wall on the measurement of hemoglobin concentrations in normal breast tissue and cancer.

METHODS

We measured total hemoglobin (tHb) concentration in both cancer and contralateral normal breast using a near-infrared time-resolved spectroscopy system in 24 female patients with breast cancer. Patients were divided into two groups based on menopausal state. The skin-to-chest wall distance was determined using ultrasound images obtained with an ultrasound probe attached to the spectroscopy probe.

RESULTS

The apparent tHb concentration of normal breast increased when the skin-to-chest wall distance was less than 20 mm. The tHb concentration in pre-menopausal patients was higher than that in post-menopausal patients. Although the concentration of tHb in cancer tissue was statistically higher than that in normal breast, the contralateral normal breast showed higher tHb concentration than cancer in 9 of 46 datasets. When the curves of tHb concentrations as a function of the skin-to-chest wall distance in normal breast were applied for pre- and post-menopausal patients separately, all the cancer lesions plotted above the curves.

CONCLUSIONS

The skin-to-chest wall distance affected the measurement of tHb concentration of breast tissue by near-infrared time-resolved spectroscopy. The tHb concentration of breast cancer tissue was more precisely evaluated by considering the skin-to-chest wall distance.

Cerebral metabolic rate of oxygen (CMRO2) assessed by combined Doppler and spectroscopic OCT

A method of measuring cortical oxygen metabolism in the mouse brain that uses independent quantitative measurements of three key parameters: cerebral blood flow (CBF), arteriovenous oxygen extraction (OE), and hemoglobin concentration ([HbT]) is presented. Measurements were performed using a single visible light spectral/Fourier domain OCT microscope, with Doppler and spectroscopic capabilities, through a thinned-skull cranial window in the mouse brain. Baseline metabolic measurements in mice are shown to be consistent with literature values. Oxygen consumption, as measured by this method, did not change substantially during minor changes either in the fraction of inspired oxygen (FiO2) or in the fraction of inspired carbon dioxide (FiCO2), in spite of larger variations in oxygen saturations. This set of experiments supports, but does not prove, the validity of the proposed method of measuring brain oxygen metabolism.

Fast single photon avalanche photodiode-based time-resolved diffuse optical tomography scanner

Resolution in diffuse optical tomography (DOT) is a persistent problem and is primarily limited by high degree of light scatter in biological tissue. We showed previously that the reduction in photon scatter between a source and detector pair at early time points following a laser pulse in time-resolved DOT is highly dependent on the temporal response of the instrument. To this end, we developed a new single-photon avalanche photodiode (SPAD) based time-resolved DOT scanner. This instrument uses an array of fast SPADs, a femto-second Titanium Sapphire laser and single photon counting electronics. In combination, the overall instrument temporal impulse response function width was 59 ps. In this paper, we report the design of this instrument and validate its operation in symmetrical and irregularly shaped optical phantoms of approximately small animal size. We were able to accurately reconstruct the size and position of up to 4 absorbing inclusions, with increasing image quality at earlier time windows. We attribute these results primarily to the rapid response time of our instrument. These data illustrate the potential utility of fast SPAD detectors in time-resolved DOT.

An extended analytical approach for diffuse optical imaging

In this work, we introduce an analytical method to solve the diffusion equation in a cylindrical geometry. This method is based on an integral approach to derive the Green's function for specific boundary conditions. Using our approach, we obtain comprehensive analytical solutions with the Robin boundary condition for diffuse optical imaging in both two and three dimensions. The solutions are expressed in terms of the optical properties of tissue and the amplitude and position of the light source. Our method not only works well inside the tissue but provides very accurate results near the tissue boundaries as well. The results obtained by our method are first compared with those obtained by a conventional analytical method then validated using numerical simulations. Our new analytical method allows not only implementation of any boundary condition for a specific problem but also fast simulation of light propagation making it very suitable for iterative image reconstruction algorithms.

Pressure modulation algorithm to separate cerebral hemodynamic signals from extracerebral artifacts

We introduce and validate a pressure measurement paradigm that reduces extracerebral contamination from superficial tissues in optical monitoring of cerebral blood flow with diffuse correlation spectroscopy (DCS). The scheme determines subject-specific contributions of extracerebral and cerebral tissues to the DCS signal by utilizing probe pressure modulation to induce variations in extracerebral blood flow. For analysis, the head is modeled as a two-layer medium and is probed with long and short source-detector separations. Then a combination of pressure modulation and a modified Beer-Lambert law for flow enables experimenters to linearly relate differential DCS signals to cerebral and extracerebral blood flow variation without a priori anatomical information. We demonstrate the algorithm's ability to isolate cerebral blood flow during a finger-tapping task and during graded scalp ischemia in healthy adults. Finally, we adapt the pressure modulation algorithm to ameliorate extracerebral contamination in monitoring of cerebral blood oxygenation and blood volume by near-infrared spectroscopy.

Evaluation of light scattering and absorption properties of in vivo rat liver using a single-reflectance fiber probe during preischemia, ischemia-reperfusion, and postmortem

Diffuse reflectance spectroscopy (DRS) has been extensively used for characterization of biological tissues as a noninvasive optical technique to evaluate the optical properties of tissue. We investigated a method for evaluating the reduced scattering coefficient μ(s)', the absorption coefficient μ(a), the tissue oxygen saturation StO₂, and the reduction of heme aa3 in cytochrome c oxidase CcO of in vivo liver tissue using a single-reflectance fiber probe with two source-collector geometries. We performed in vivo recordings of diffuse reflectance spectra for exposed rat liver during the ischemia-reperfusion induced by the hepatic portal (hepatic artery, portal vein, and bile duct) occlusion. The time courses of μ a at 500, 530, 570, and 584 nm indicated the hemodynamic change in liver tissue as well as StO₂. Significant increase in μ(a)(605)/μ(a)(620) during ischemia and after euthanasia induced by nitrogen breathing was observed, which indicates the reduction of heme aa3, representing a sign of mitochondrial energy failure. The time courses of μ(s)' at 500, 530, 570, and 584 nm were well correlated with those of μ(a), which also reflect the scattering by red blood cells. On the other hand, at 700 and 800 nm, a temporary increase in μ(s)' and an irreversible decrease in μ(s)' were observed during ischemia-reperfusion and after euthanasia induced by nitrogen breathing, respectively. The change in μ(s)' in the near-infrared wavelength region during ischemia is indicative of the morphological changes in the cellular and subcellular structures induced by the ischemia, whereas that after euthanasia implies the hepatocyte vacuolation. The results of the present study indicate the potential application of the current DRS system for evaluating the pathophysiological conditions of in vivo liver tissue.

Multigrid-based reconstruction algorithm for quantitative photoacoustic tomography

This paper proposes a multigrid inversion framework for quantitative photoacoustic tomography reconstruction. The forward model of optical fluence distribution and the inverse problem are solved at multiple resolutions. A fixed-point iteration scheme is formulated for each resolution and used as a cost function. The simulated and experimental results for quantitative photoacoustic tomography reconstruction show that the proposed multigrid inversion can dramatically reduce the required number of iterations for the optimization process without loss of reliability in the results.

Explicit solutions of the radiative transport equation in the P3 approximation

PURPOSE

Explicit solutions of the monoenergetic radiative transport equation in the P3 approximation have been derived which can be evaluated with nearly the same computational effort as needed for solving the standard diffusion equation (DE). In detail, the authors considered the important case of a semi-infinite medium which is illuminated by a collimated beam of light.

METHODS

A combination of the classic spherical harmonics method and the recently developed method of rotated reference frames is used for solving the P3 equations in closed form.

RESULTS

The derived solutions are illustrated and compared to exact solutions of the radiative transport equation obtained via the Monte Carlo (MC) method as well as with other approximated analytical solutions. It is shown that for the considered cases which are relevant for biomedical optics applications, the P3 approximation is close to the exact solution of the radiative transport equation.

CONCLUSIONS

The authors derived exact analytical solutions of the P3 equations under consideration of boundary conditions for defining a semi-infinite medium. The good agreement to Monte Carlo simulations in the investigated domains, for example, in the steady-state and time domains, as well as the short evaluation time needed suggests that the derived equations can replace the often applied solutions of the diffusion equation for the homogeneous semi-infinite medium.

Breast tissue composition and its dependence on demographic risk factors for breast cancer: non-invasive assessment by time domain diffuse optical spectroscopy

BACKGROUND

Breast tissue composition is recognized as a strong and independent risk factor for breast cancer. It is a heritable feature, but is also significantly affected by several other elements (e.g., age, menopause). Nowadays it is quantified by mammographic density, thus requiring the use of ionizing radiation. Optical techniques are absolutely non-invasive and have already proved effective in the investigation of biological tissues, as they are sensitive to tissue composition and structure.

METHODS

Time domain diffuse optical spectroscopy was performed at 7 wavelengths (635-1060 nm) on 200 subjects to derive their breast tissue composition (in terms of water, lipid and collagen content), blood parameters (total hemoglobin content and oxygen saturation level), and information on the microscopic structure (scattering amplitude and power). The dependence of all optically-derived parameters on age, menopausal status, body mass index, and use of oral contraceptives, and the correlation with mammographic density were investigated.

RESULTS

Younger age, premenopausal status, lower body mass index values, and use of oral contraceptives all correspond to significantly higher water, collagen and total hemoglobin content, and lower lipid content (always p < 0.05 and often p < 10-4), while oxygen saturation level and scattering parameters show significant dependence only on some conditions. Even when age-adjusted groups of subjects are compared, several optically derived parameters (and in particular always collagen and total hemoglobin content) remain significantly different.

CONCLUSIONS

Time domain diffuse optical spectroscopy can probe non-invasively breast tissue composition and physiologic blood parameters, and provide information on tissue structure. The measurement is suitable for in vivo studies and monitoring of changes in breast tissue (e.g., with age, lifestyle, chemotherapy, etc.) and to gain insight into related processes, like the origin of cancer risk associated with breast density.

Fast method of retrieving the asymmetry factor and scattering albedo from the maximum time-resolved reflectance of participating media

This research presents a parametric study of the time-resolved hemispherical reflectance of a semi-infinite plane-parallel slab of homogeneous, nonemitting, absorbing, and anisotropic scattering medium exposed to a collimated Gaussian pulse. The one-dimensional transient radiative transfer equation was solved by using the finite volume method. The internal reflection at the medium-air interface caused by the mismatch of the refractive indices was considered. In particular, this work focused on the maximum diffuse hemispherical reflectance. Three different optical regions were identified according to the dimensionless pulsewidth βct_p. The correlation between the normalized maximum hemispherical reflectance and βct_p was conformed to the Boltzmann function. The coefficients in the correlating functions of the match and mismatch refractive index cases were fitted as polynomial fitting functions of the single scattering albedo ω and Henyey-Greenstein asymmetric factor g. Thus, ω and g can be simultaneously reconstructed by the semi-empirical correlations without solving the forward model. In conclusion, the proposed method can potentially retrieve the asymmetry factor and single scattering albedo of participating media accurately and efficiently.

Fast silicon photomultiplier improves signal harvesting and reduces complexity in time-domain diffuse optics

We present a proof of concept prototype of a time-domain diffuse optics probe exploiting a fast Silicon PhotoMultiplier (SiPM), featuring a timing resolution better than 80 ps, a fast tail with just 90 ps decay time-constant and a wide active area of 1 mm². The detector is hosted into the probe and used in direct contact with the sample under investigation, thus providing high harvesting efficiency by exploiting the whole SiPM numerical aperture and also reducing complexity by avoiding the use of cumbersome fiber bundles. Our tests also demonstrate high accuracy and linearity in retrieving the optical properties and suitable contrast and depth sensitivity for detecting localized inhomogeneities. In addition to a strong improvement in both instrumentation cost and size with respect to legacy solutions, the setup performances are comparable to those of state-of-the-art time-domain instrumentation, thus opening a new way to compact, low-cost and high-performance time-resolved devices for diffuse optical imaging and spectroscopy.

Polymer-encapsulated organic nanoparticles for fluorescence and photoacoustic imaging

Polymer encapsulated organic nanoparticles have recently attracted increasing attention in the biomedical field because of their unique optical properties, easy fabrication and outstanding performance as imaging and therapeutic agents. Of particular importance is the polymer encapsulated nanoparticles containing conjugated polymers (CP) or fluorogens with aggregation induced emission (AIE) characteristics as the core, which have shown significant advantages in terms of tunable brightness, superb photo- and physical stability, good biocompatibility, potential biodegradability and facile surface functionalization. In this review, we summarize the latest advances in the development of polymer encapsulated CP and AIE fluorogen nanoparticles, including preparation methods, material design and matrix selection, nanoparticle fabrication and surface functionalization for fluorescence and photoacoustic imaging. We also discuss their specific applications in cell labeling, targeted in vitro and in vivo imaging, blood vessel imaging, cell tracing, inflammation monitoring and molecular imaging. We specially focus on strategies to fine-tune the nanoparticle property (e.g. size and fluorescence quantum yield) through precise engineering of the organic cores and careful selection of polymer matrices. The review also highlights the merits and limitations of these nanoparticles as well as strategies used to overcome the limitations. The challenges and perspectives for the future development of polymer encapsulated organic nanoparticles are also discussed.

Human brown adipose tissue assessed by simple, noninvasive near-infrared time-resolved spectroscopy

OBJECTIVE

Human brown adipose tissue (BAT) activity has been typically evaluated by (18) F-fluorodeoxyglucose (FDG)-positron emission tomography (PET) combined with computed tomography (CT). However, FDG-PET/CT has serious limitations (e.g., radiation and cold exposure). This study evaluated BAT density using near-infrared time-resolved spectroscopy (NIRTRS ), a simple and noninvasive method of measuring the indices of tissue hemoglobin concentration [total-Hb] and mitochondrial density (µs ').

METHODS

The NIRTRS parameters at 760, 800, and 830 nm in the supraclavicular region potentially containing BAT were evaluated. First, the NIRTRS parameters were compared at 27 °C and during a 2-h cold exposure (19 °C) in 18 men. Then, NIRTRS parameters at 27 °C were compared with mean standardized uptake values (SUVmean ) assessed by FDG-PET/CT after the 2-h cold exposure (19 °C) in 29 men.

RESULTS

There was no significant difference between the NIRTRS parameters at 27 °C and 19°C. The [total-Hb] and µs ' were significantly correlated to SUVmean (r = 0.73 and r = 0.64, respectively). A receiver operating characteristic analysis revealed that the sensitivity (75.0-82.4%), specificity (91.7-100%), and accuracy (82.8-86.2%) of the NIRTRS parameters were all good to determine the NIRTRS reliability.

CONCLUSIONS

Our novel NIRTRS method is noninvasive and simple and can reliably assess human BAT density in the supraclavicular region.

Towards next-generation time-domain diffuse optics for extreme depth penetration and sensitivity

Light is a powerful tool to non-invasively probe highly scattering media for clinical applications ranging from oncology to neurology, but also for molecular imaging, and quality assessment of food, wood and pharmaceuticals. Here we show that, for a paradigmatic case of diffuse optical imaging, ideal yet realistic time-domain systems yield more than 2-fold higher depth penetration and many decades higher contrast as compared to ideal continuous-wave systems, by adopting a dense source-detector distribution with picosecond time-gating. Towards this aim, we demonstrate the first building block made of a source-detector pair directly embedded into the probe based on a pulsed Vertical-Cavity Surface-Emitting Laser (VCSEL) to allow parallelization for dense coverage, a Silicon Photomultiplier (SiPM) to maximize light harvesting, and a Single-Photon Avalanche Diode (SPAD) to demonstrate the time-gating capability on the basic SiPM element. This paves the way to a dramatic advancement in terms of increased performances, new high impact applications, and availability of devices with orders of magnitude reduction in size and cost for widespread use, including quantitative wearable imaging.

Study of the effect introduced by an integrating sphere on the temporal profile characterization of short laser pulses propagating through a turbid medium

When a nanosecond laser pulse is transmitted through a highly scattering material, its irradiance decreases as it propagates; this is because of the spatial and temporal pulse profile stretching owing to multiple scattering events. Although the effect of temporal distortion is much less significant than that of the spatial distortion for applications where the laser beam is focused on a subsurface target (writing of waveguides, for example), it becomes significant for applications where the laser pulse must attain certain temporal width after the beam propagated is collimated through a turbid medium (photoacoustic tomography, for example). The objective of this work is to determine the transfer function associated to an integrating sphere measurement of the temporal intensity profile involving turbid media samples. The transfer function is found to be related to the geometrical characteristics of the integrating sphere and the optical properties of the turbid media. This procedure opens a new possibility for optical property characterization and enables the use of an integrating sphere for time-dependent intensity measurements.

Novel analytical solution for the radiance in an anisotropically scattering medium

We report on a novel analytical solution of the three-dimensional radiative transport equation for the case of an infinitely extended anisotropically scattering medium that is illuminated by an isotropic point light source. The resulting expression for the radiance can be evaluated efficiently and accurately and exhibits significant improvements with respect to the convergence and the numerical stability compared to the solutions found in the literature so far. The equations obtained were successfully verified by comparisons with the Monte Carlo method.

Computationally effective solution of the inverse problem in time-of-flight spectroscopy

Photon time-of-flight (PTOF) spectroscopy enables the estimation of absorption and reduced scattering coefficients of turbid media by measuring the propagation time of short light pulses through turbid medium. The present investigation provides a comparison of the assessed absorption and reduced scattering coefficients from PTOF measurements of intralipid 20% and India ink-based optical phantoms covering a wide range of optical properties relevant for biological tissues and dairy products. Three different models are used to obtain the optical properties by fitting to measured temporal profiles: the Liemert-Kienle model (LKM), the diffusion model (DM) and a white Monte-Carlo (WMC) simulation-based algorithm. For the infinite space geometry, a very good agreement is found between the LKM and WMC, while the results obtained by the DM differ, indicating that the LKM can provide accurate estimation of the optical parameters beyond the limits of the diffusion approximation in a computational effective and accurate manner. This result increases the potential range of applications for PTOF spectroscopy within industrial and biomedical applications.

In situ fiber-optical monitoring of cytosolic calcium in tissue explant cultures

We present a fluorescence-lifetime based method for monitoring cell and tissue activity in situ, during cell culturing and in the presence of a strong autofluorescence background. The miniature fiber-optic probes are easily incorporated in the tight space of a cell culture chamber or in an endoscope. As a first application we monitored the cytosolic calcium levels in porcine tracheal explant cultures using the Calcium Green-5N (CG5N) indicator. Despite the simplicity of the optical setup we are able to detect changes of calcium concentration as small as 2.5 nM, with a monitoring time resolution of less than 1 s.

Impact of ischemic preconditioning on functional sympatholysis during handgrip exercise in humans

Repeated bouts of ischemia followed by reperfusion, known as ischemic preconditioning (IPC), is found to improve exercise performance. As redistribution of blood from the inactive areas to active skeletal muscles during exercise (i.e., functional sympatholysis) is important for exercise performance, we examined the hypothesis that IPC improves functional sympatholysis in healthy, young humans. In a randomized study, 15 healthy young men performed a 10-min resting period, dynamic handgrip exercise at 10% maximal voluntary contraction (MVC), and 25% MVC. This protocol was preceded by IPC (IPC; 4 × 5-min 220-mmHg unilateral occlusion) or a sham intervention (CON; 4 × 5-min 20-mmHg unilateral occlusion). Near-infrared spectroscopy was used to assess changes in oxygenated hemoglobin and myoglobin in skeletal muscle (HbO₂ + MbO₂) in response to sympathetic activation (via cold pressor test (CPT)) at baseline and during handgrip exercise (at 10% and 25%). In resting conditions, HbO₂ + MbO₂ significantly decreased during CPT (−11.0 ± 1.0%), which was significantly larger during the IPC-trial (−13.8 ± 1.2%, P = 0.006). During handgrip exercise at 10% MVC, changes in HbO₂ + MbO₂ in response to the CPT were blunted after IPC (−8.8 ± 1.5%) and CON (−8.3 ± 0.4%, P = 0.593). During handgrip exercise at 25% MVC, HbO₂ + MbO₂ in response to the CPT increased (2.0 ± 0.4%), whereas this response was significantly larger when preceded by IPC (4.2 ± 0.6%, P = 0.027). Collectively, these results indicate that IPC-induced different vascular changes at rest and during moderate exercise in response to sympathetic activation. This suggests that, in healthy volunteers, exposure to IPC may alter tissue oxygenation during sympathetic stimulation at rest and during exercise.