|
1.IntroductionThe near-infrared (NIR) optical properties of tissue, namely absorption and reduced scattering coefficients, can provide information on a variety of tissue physiology processes. Wavelength-dependent absorption is widely used to quantify the concentration of biologically important chromophores, such as hemoglobin, water, and lipid. To date, scattering properties are not fully understood. Basically, the wavelength dependence of scattering is associated with tissue structure. In the past decades, increasing use of light either as a diagnostic or therapeutic tool has demanded a better knowledge of the optical properties of human tissues.1 2 3 4 The properties of the head, especially the neonatal brain, are not particularly well known. In recent years, great effort has been devoted to evaluating the optical properties of brain tissue. Using integrating-sphere measurements, van der Zee, Essenpreis, and Delpy5 determined the optical properties of postmortem neonate and adult human brain tissue over the wavelength range of 500 to 1000 nm. Yaroslavsky et al.6 reported the optical properties of brain tissue in vitro from 360 to 1100 nm. However, due to significant alternation that may occur after the death of the organism, to obtain reliable information, it is important to estimate the optical properties in vivo. The time-resolved method has proven to be a valuable tool to assess tissue optical properties in vivo noninvasively. Recently, in vivo optical properties of piglet brain7 and adult forehead8 9 were measured by time-resolved systems. Hebden et al.10 obtained the first 3-D images of a premature infant brain using a 32-channel time-resolved imaging system. However, these measurements were obtained with a complex time-resolved system not suitable for clinical use. A frequency-domain apparatus, on the contrary, is relatively inexpensive and can easily be made portable, which are two important considerations for a clinically useful instrument. Fantini et al.11 and Hueber et al.12 reported the optical properties of newborn piglet brain using a frequency-domain multidistance spectroscopy. Bevilacqua et al.13 measured the local optical properties of normal and malignant human brain tissue in vivo during brain surgery with short source-detector separations. More recently, Choi et al.14 published the measurements of absorption and reduced scattering coefficients through the forehead on 30 adult volunteers using a multidistance frequency-domain method. Up to now, however, published data on in vivo optical properties of neonatal brain in large populations, to the best of our knowledge, is relatively scarce. Recently, we have developed a portable, multiwavelength, frequency-domain near-infrared spectroscopy (NIRS) instrument for quantitative, noninvasive measurement of tissue optical properties.15 16 A calibration procedure was employed to quantify the absorption and reduced scattering coefficients with single source-detector separation. This instrument has been validated by extensive tissue equivalent phantom and animal experimental studies.17 18 We investigated the optical properties of 23 infant heads using our frequency-domain NIRS instrument. The physiological parameters were also derived from the absorption coefficients and compared with those available in the literature. 2.Theory and Method2.1.AlgorithmThe propagation of photons in tissue is well described by the radiation transport theory, and its simplification to the highly scattering regime results in the diffusion equation.19 20 21 In frequency-domain systems, the light source is intensity modulated by a sinusoidal wave. The light energy propagates as a damped spherical wave outward from the source. This macroscopic scalar wave is often called the diffuse photon density wave (DPDW).22 23 24 DPDW dispersion is highly dependent on the optical properties of the tissue. Experimentally, we measured the reflectance, which is the number of photons crossing the tissue boundary per unit time per unit area at distance ρ. Pogue and Patterson25 derived the expression for the ac component of reflectance R(ρ,ω,μa,μs ′,cn) by application of the extrapolated boundary condition to the diffusion equation. where and s is a constant that depends on the refractive index mismatch (s=0.4258 for a tissue-air boundary).Note that the measured reflectance R(ρ,ω,μa,μs ′,cn) is complex. The phase of the ac component of the reflectance is given by and the amplitude is The measured amplitude and phase on the infant heads were first normalized by the amplitude and phase measured on a reference phantom, then fitted in Eq. (1) (Levenberg-Marquardt nonlinear fitting) to estimate the absorption coefficient μa and reduced scattering coefficient μs ′.2.2.Derivation of Physiological ParametersWe assume that the chromophores contributing to μa in the brain tissue are principally oxy- and deoxyhemoglobin: where μa λ is the absorption coefficient at wavelength λ, ε chrom λ is the extinction coefficient of a given chromophore at wavelength λ, and [Hb] and [HbO 2] are the concentration of hemoglobin and oxygenated hemoglobin, respectively.In this study, the absorption coefficients obtained at two wavelengths (i.e., 788 and 832 nm) were employed to calculate the [Hb] and [HbO 2]. ε Hb λ and ε HbO 2 λ were obtained from the literature.26 27 Two physiologically important parameters, namely tissue oxygen saturation (StO 2) and total hemoglobin concentration (THC), could be derived from [Hb] and [HbO 2]. 3.Instrumentation and Calibration3.1.Frequency-Domain NIRS InstrumentThis study has been performed with a homodyne frequency-domain NIRS instrument (see Fig. 1), which was described in detail elsewhere.15 16 Briefly, this frequency-domain NIRS instrument uses three laser diodes at 788-, 814-, and 832-nm wavelength. The three wavelengths are time multiplexed at a rate of 1 Hz. The light is modulated at 140 MHz and directed to the forehead of neonates by a 1-mm-diam fiber guide. The diffused light is collected to the photomultiplier tube (PMT, Hamamatsu, H7683-02) by another 3-mm-diam fiber bundle, placed several centimeters apart from the source fiber. The amplitude and phase shift of the detected photon density wave can be further derived from the output of the in-phase and quadrature demodulator (Mini-Circuit, MIQY-140D). 3.2.Instrument CalibrationIt is noted that the measured phase and amplitude of photon density waves include contribution from both the sample and the instrument.16 28 The measured phase is the sum of the sample and instrument phase, while the measured amplitude is the product of the sample and instrument amplitude. Two kinds of calibration procedure can be used to remove the instrument response from the raw data: the multidistance method and use of a standard reference phantom calibration. The first approach works well in a macroscopically homogeneous turbid medium, while the latter is more effective for heterogeneous in vivo measurements.28 In this study, a standard reference phantom with known optical properties was employed to remove the instrument response from the raw data. The phantom was made of a pliable RTV silicon resin. A calculated amount of TiO 2 powder and carbon black powder were added to the resin to obtain the expected reduced scattering coefficient and absorption coefficient, respectively. The optical properties of the calibration phantom had been previously determined by measurements using the multidistance method.23 29 To match the shape of the infant heads, the calibration phantom was cylindrically molded with a radius of curvature of 5.5 cm. 4.Experimental Protocol23 infants were enrolled in this study. They had a median gestional age of 39 (range 35 to 42.5) weeks, a body weight of 3200 (2400 to 4600) g, and a postnatal age of 5 (1 to 17) days. Two of them were preterm infants with a gestional age of 35 and 36 weeks, respectively. The subjects were admitted due some diseases but were all in recovery at the time of study and had normal cranial ultrasound scans: mild jaundice (n=12), pneumonia (n=3), hypoglycemia (n=2), and healthy infants (n=6). Informed parental consent was obtained before each investigation. The subjects were in a supine, comfortable position during the entire experiment. Measurements were performed on each subject by gently placing the probe on the forehead. Just enough pressure was applied to ensure optical contact between skin and probe. At least two measurements were performed successively at a given location to exclude the problem of optical coupling. Sometimes we changed the position a little bit for better contact. The source-detector separation of the probe was fixed at 4 cm, which means the frequency-domain NIRS instrument views a banana-shaped tissue volume between the emitter and detector located approximately 1 to 3 cm beneath the skin.3 9 With the probe on the forehead, this instrument appears to monitor both gray and white matter in the frontal neocortex. The calibration procedures were performed right before or after the measurements. The top and bottom sides of the cylindrical silicon resin phantom were covered with aluminum foil to simulate the semi-infinite boundary condition. 5.Results and Discussion5.1.Optical PropertiesFigure 2 shows the histograms of the fitted absorption coefficient μa and reduced scattering coefficient μs ′ in 23 infant foreheads at wavelengths of 788 and 832 nm, respectively. The vertical axis in each histogram is the number of subjects. Mean and standard deviation of the fitted optical properties are listed in Table 1. Table 1
Note that the absorption coefficients of the infant foreheads are lower than the values reported in adults.8 9 14 This difference may partly be due to the brain itself, and partly due to the influence of the superficial layer. NIRS monitors the tissue beneath the optical fiber, which includes scalp (skin, subcutaneous fat, muscle), skull, cerebrospinal fluid (CSF), and the brain. The measurement of an adult head may be contaminated by these layered structures. In infants, however, the thin extracerebral tissues (<5 mm) do not significantly contaminate the measurement of optical properties from the surface of the head.30 The broadly distributed μs ′ in Fig. 2 may be due to the large intersubject variation. In other words, the variation of μs ′ should be taken into account when comparing data among subjects. However, the current commercially available continuous-wave NIRS cannot separate the effect of scattering from the attenuation of light and thus simply assumes a constant μs ′. Frequency-domain NIRS allows separation of μa and μs ′, which make it possible to provide absolute concentration of Hb and HbO 2. 5.2.Physiological ParametersHistograms of the physiological properties extracted from μa values are summarized in Fig. 3. Note that we have ignored the background absorption and used absorption coefficients at two wavelengths to derive the hemoglobin concentration. Despite the simplistic assumptions involved, the physiological properties derived from the absorption coefficients (Table 1) agree reasonably well with the results available in the literature. Cooper et al.31 measured 19 newborn infants with second differential spectroscopy. Their preliminary results gave a mean deoxyhemoglobin concentration of 14.6±4.0 μM, which is consistent with our result (16.0±5.1 μM). NIRS measurements of absolute cerebral hemoglobin concentrations in neonates by Wolf et al.32 found an average tissue oxygen saturation (StO 2) of about 66±4%, while we obtain a somewhat lower value of 58.7±11.2%. Wyatt et al.33 reported the mean cerebral blood volume of the infant brain (2.2±0.4 ml per 100 g of tissue). This corresponds to an estimated total hemoglobin concentration (THC) of 36±11 μM of tissue, which agrees well with our average of 39.7±9.8 μM. 5.3.Reproducibility of the Optical MeasurementsThe reproducibility of our optical measurements were assessed by successive reposition of the optical probe on the same spot. 15 infants were measured just after they were fed. They were all quiet during the measurements. In fact, most of them were asleep. Since optical coupling is one of the main sources of reproducibility errors, all measurements were performed by the same experimenter using a consistent pressure. Five measurements were performed consecutively per infant. Reproducibility is based on the maximum difference in the parameters as a result of repeated repositioning of the optical probe on the same spot.11 The reproducibility errors estimate the uncertainty in the absolute readings. The average results of 15 infants are listed in Table 2. It suggests that the results were repeatable in quiet situations. However, an optical probe that could tolerate extensive movements is still desirable, especially in the study of the infants. Table 2
6.DiscussionThe main clinical interest in NIRS was initially for monitoring brain oxygenation in newborn babies. NIRS has been used for perinatal application since 1985, and hundreds of studies have been performed.34 35 Used in a research context, NIRS has improved understanding of cerebral circulation. However, NIRS is used infrequently to guide clinical care because of quantitation problems related to the biophysics of measurement. The improved NIRS technologies, such as time- and frequency-resolved NIRS enables experimenters to separate tissue scattering effects from tissue absorption effects, and thus extract more quantitative information about tissue chromophores. In this work, a novel frequency-domain method was employed to quantify the optical properties of infant brains. However, several points should be considered during the interpretation of the experimental data. First, the measurements were obtained with a single source-detector separation (4 cm). The NIRS instrument views a banana-shaped tissue volume including the scalp, skull, cerebrospinal fluid (CSF), and the brain. In newborn infants, the thin extracerebral tissues make it possible to detect more information from both gray and white matter in the frontal neocortex.30 In adults, however, the interference from the thicker extracranial tissues cannot be ignored.36 To obtain depth information on different layers, the multiple-distance method could be introduced in adult brain applications.14 37 Secondly, a simplified model of light propagation (homogeneity and semi-infinite boundary condition) was employed to recover the optical properties of tissue from the experimental data. The neonatal head, however, was both highly inhomogeneous and curved. The difference between the actual anatomy and model may induce systematic errors. Analytical solutions of the diffusion equation for more accurate sphere and cylinder geometries are also available in the literature.25 38 Unfortunately, they are time consuming and not suitable for real-time application. Finally, we assumed that oxy- and deoxyhemoglobin are the primary chromophores in the brain when deriving the physiological properties from μa at two wavelengths (788 and 832 nm). In brain tissue, hemoglobin, water, and cytochrome-c oxidase are known chromophores in the near-infrared region.27 The accuracy of quantitative hemoglobin concentration may be hindered by ignoring the nonhemoglobin absorption. There are also some practical limitations to obtaining the quantitative measurements of infant brains, which have not been discussed here. The main concern of this work is to investigate the potential of this novel frequency-domain NIRS in pediatric applications. The baseline of the optical and physiological properties of infant brains was obtained. Work is in progress to improve the instrument and to design more complex measurement protocols for assessment of infant brain injury and neurodevelopment. 7.ConclusionsWe present determination of in vivo optical properties of 23 neonatal foreheads using a portable frequency-domain NIRS instrument. A calibration phantom of known optical properties is employed to quantify the absorption coefficient and reduced scattering coefficient with a single source-detector separation. The ability to separate tissue scattering effects from tissue absorption effects makes it possible to obtain quantitative measurements. Physiological parameters derived from the measured absorption coefficients are consistent with those in the literature. The factors affecting the interpretation of the experimental data are also discussed, which suggests considerable work remains to assess the utility of this approach in clinical practice. REFERENCES
A. Villringer
and
B. Chance
,
“Non-invasive optical spectroscopy and imaging of human brain function,”
Trends Neurosci. , 20
(10), 435
–442
(1997). Google Scholar
T. H. Pham
,
O. Coquoz
,
J. B. Fishkin
,
E. Anderson
, and
B. J. Tromberg
,
“Broad bandwidth frequency domain instrument for quantitative tissue optical spectroscopy,”
Rev. Sci. Instrum. , 71
(6), 2500
–2513
(2000). Google Scholar
E. Gratton
,
S. Fantini
,
M. A. Franceschini
,
G. Gratton
, and
M. Fabiani
,
“Measurement of scattering and absorption changes in muscle and brain,”
Philos. Trans. R. Soc. London, Ser. B , 352 727
–735
(1997). Google Scholar
M. Wolf
,
U. Wolf
,
J. H. Choi
,
V. Toronov
,
L. A. Paunescu
,
A. Michalos
, and
E. Gratton
,
“Fast cerebral functional signal in the 100-ms range detected in the visual cortex by frequency-domain near-infrared spectrophotometry,”
Psychophysiology , 40
(4), 521
–528
(2003). Google Scholar
P. van der Zee
,
M. Essenpreis
, and
D. T. Delpy
,
“Optical properties of brain tissue,”
Proc. SPIE , 1888 454
–465
(1993). Google Scholar
A. N. Yaroslavsky
,
P. C. Schulze
,
I. V. Yaroslavsky
,
R. Schober
,
F. Ulrich
, and
H. J. Schwarzmaier
,
“Optical properties of selected native and coagulated human brain tissues in vitro in the visible and near infrared spectral range,”
Phys. Med. Biol. , 47 2059
–2073
(2002). Google Scholar
A. Sassaroli
,
F. Martelli
,
Y. Tanikawa
,
K. Tanaka
,
R. Araki
,
Y. Onodera
, and
Y. Yamada
,
“Time-resolved measurements of in vivo optical properties of piglet brain,”
Opt. Rev. , 7
(5), 420
–425
(2000). Google Scholar
S. J. Matcher
,
M. Cope
, and
D. T. Delpy
,
“In vivo measurements of the wavelength dependence of tissue-scattering coefficients between 760 and 900 nm measured with time-resolved spectroscopy,”
Appl. Opt. , 36
(1), 386
–396
(1997). Google Scholar
A. Torricelli
,
A. Pifferi
,
P. Taroni
,
E. Giambattistelli
, and
R. Cubeddu
,
“In vivo optical characterization of human tissues from 610 to 1010 nm by time-resolved reflectance spectroscopy,”
Phys. Med. Biol. , 46 2227
–2237
(2001). Google Scholar
J. C. Hebden
,
A. Gibson
,
R. M. Yusof
,
N. Everdell
,
E. M. C. Hillman
,
D. T. Delpy
,
S. R. Arridge
,
T. Austin
,
J. H. Meek
, and
J. S. Wyatt
,
“Three-dimensional optical tomography of the premature infant brain,”
Phys. Med. Biol. , 47 4155
–4166
(2002). Google Scholar
S. Fantini
,
D. Hueber
,
M. A. Franceschini
,
E. Gratton
,
W. Rosenfeld
,
P. G. Stubblefield
,
D. Maulk
, and
M. R. Stankovic
,
“Non-invasive optical monitoring of the newborn piglet brain using continous-wave and frequency-domain spectroscopy,”
Phys. Med. Biol. , 44 1543
–1563
(1999). Google Scholar
D. M. Hueber
,
M. A. Franceschini
,
H. Y. Ma
,
Q. Zhang
,
J. R. Ballesteros
,
S. Fantini
,
D. Wallace
,
V. Ntziachristos
, and
B. Chance
,
“Non-invasive and quantitative near-infrared haemoglobin spectrometry in the piglet brain during hypoxic stress, using a frequency-domain multidistance instrument,”
Phys. Med. Biol. , 46 41
–62
(2001). Google Scholar
F. Bevilacqua
,
D. Piguet
,
P. Marquet
,
J. D. Gross
,
B. J. Tromberg
, and
C. Depeursinge
,
“In vivo local determination of tissue optical properties: applications to human brain,”
Appl. Opt. , 38
(22), 4939
–4950
(1999). Google Scholar
J. Choi
,
M. Wolf
,
V. Toronov
,
U. Wolf
,
C. Polzonetti
,
D. Hueber
,
L. P. Safonova
,
R. Gupta
,
A. Michalos
,
W. Mantulin
, and
E. Gratton
,
“Noninvasive determination of the optical properties of adult brain: near-infrared spectroscopy approach,”
J. Biomed. Opt. , 9
(1), 221
–229
(2004). Google Scholar
Y. S. Yang
,
H. L. Liu
,
X. D. Li
, and
B. Chance
,
“Low-cost frequency-domain photon migration instrument for tissue spectroscopy, oximetry, and imaging,”
Opt. Eng. , 36
(5), 1562
–1569
(1997). Google Scholar
B. Chance
,
M. Cope
,
E. Gratton
,
N. Ramanujam
, and
B. Tromberg
,
“Phase measurement of light absorption and scatter in human tissue,”
Rev. Sci. Instrum. , 69
(10), 3457
–3481
(1998). Google Scholar
H. Y. Ma
,
Q. Xu
,
J. R. Ballesteros
,
V. Ntziachristos
,
Q. Zhang
, and
B. Chance
,
“Quantitative study of hypoxia stress in piglet brain by IQ phase modulation oximeter,”
Proc. SPIE , 3597 642
–649
(1999). Google Scholar
Y. Chen
,
R. T. Dharmesh
,
X. Intes
, and
B. Chance
,
“Correlation between near-infrared spectroscopy and magnetic resonance imaging of rat brain oxygenation modulation,”
Phys. Med. Biol. , 48
(4), 417
–427
(2003). Google Scholar
M. S. Patterson
,
B. Chance
, and
B. C. Wilson
,
“Time resolved reflectance and transmittance for the non-invasive measurement of tissue optical properties,”
Appl. Opt. , 26
(12), 2331
–2336
(1989). Google Scholar
M. S. Patterson
,
J. D. Moulton
,
B. C. Wilson
,
K. W. Berndt
, and
J. R. Lakowicz
,
“Frequency-domain reflectance for the determination of the scattering and absorption properties of tissue,”
Appl. Opt. , 30
(31), 4474
–4476
(1991). Google Scholar
R. C. Haskell
,
L. O. Svaasand
,
T. T. Tsay
,
T. C. Feng
,
M. S. McAdams
, and
B. J. Tromberg
,
“Boundary conditions for the diffusion equation in radiative transfer,”
J. Opt. Soc. Am. A , 11
(10), 2727
–2741
(1994). Google Scholar
M. A. O’Leary
,
D. A. Boas
,
B. Chance
, and
A. G. Yodh
,
“Refraction of diffuse photon density waves,”
Phys. Rev. Lett. , 69 2658
(1992). Google Scholar
J. B. Fishkin
and
E. Gratton
,
“Propagation of photon-density waves in strongly scattering media containing an absorbing semi-infinite plane bounded by a straight edge,”
J. Opt. Soc. Am. A , 10
(1), 127
–140
(1993). Google Scholar
B. J. Tromberg
,
L. O. Svaasand
,
T. T. Tsay
, and
R. C. Haskell
,
“Properties of photon density waves in multiple-scattering media,”
Appl. Opt. , 32 607
–616
(1993). Google Scholar
B. W. Pogue
and
M. S. Patterson
,
“Frequency-domain optical-absorption spectroscopy of finite tissue volumes using diffusion-theory,”
Phys. Med. Biol. , 39 1157
–1180
(1994). Google Scholar
S. J. Matcher
,
C. E. Elwell
,
C. E. Cooper
,
M. Cope
, and
D. T. Delpy
,
“Performance comparison of several published tissue near-infrared sepctroscopy algorithms,”
Anal. Biochem. , 227 54
–68
(1995). Google Scholar
B. J. Tromberg
,
O. Coquoz
,
J. B. Fishkin
,
T. Pham
,
E. Anderson
,
J. Butler
,
M. Cahn
,
J. D. Gross
,
V. Venugopalan
, and
D. Pham
,
“Non-invasive measurements of breast tissue optical properties using frequency-domain photon migration,”
Philos. Trans. R. Soc. London, Ser. B , 352 661
–668
(1997). Google Scholar
S. Fantini
,
M. A. Franceschini
, and
E. Gratton
,
“Semi-infinite-geometry boundary problem for light migration in highly scattering media: a frequency-domain study in the diffusion approximation,”
J. Opt. Soc. Am. B , 11
(10), 2128
–2138
(1994). Google Scholar
C. D. Kurth
and
W. S. Thayer
,
“A multiwavelength frequency-domain near-infrared cerebral oximeter,”
Phys. Med. Biol. , 44 727
–740
(1999). Google Scholar
C. E. Cooper
,
C. E. ElWell
,
J. H. Meek
,
S. J. Matcher
,
J. S. Wyatt
,
M. Cope
, and
D. T. Delpy
,
“The noninvasive measurement of absolute cerebral deoxyhemoglobin concentration and mean optical path length in the neonatal brain by second derivative near infrared spectroscopy,”
Pediatr. Res. , 39
(1), 32
–38
(1996). Google Scholar
M. Wolf
,
P. Evans
,
H. U. Bucher
,
V. Dietz
,
M. Keel
,
R. Strebel
, and
K. von Siebenthal
,
“The measurement of absolute cerebral hemoglobin concentration in adults and neonates,”
Adv. Exp. Med. Biol. , 428 219
–227
(1997). Google Scholar
J. S. Wyatt
,
M. Cope
,
D. T. Delpy
,
C. E. Richardson
,
A. D. Edwards
,
S. Wray
, and
E. O. R. Reynolds
,
“Quantitation of cerebral blood volume in human infants by near-infrared spectroscopy,”
J. Appl. Phys. , 68
(3), 1086
–1091
(1990). Google Scholar
J. E. Brazy
,
V. Darrell
,
M. D. Lewis
,
M. H. Mitnick
, and
F. F. Jobsis
,
“Nonivasive monitoring of cerebral oxygenation in preterm infants: preliminary observations,”
Pediatrics , 75 217
–225
(1985). Google Scholar
S. E. Nicklin
,
I. A. A. Hassan
,
Y. A. Wickramasinghe
, and
S. A. Spencer
,
“The light still shines, but not that brightly? The current status of perinatal near infrared spectroscopy,”
Arch. Dis. Child , 88
(4), 263
–268
(2003). Google Scholar
E. Okada
,
M. Firbank
,
M. Schweiger
,
S. R. Arridge
,
M. Cope
, and
D. T. Delpy
,
“Theoretical and experimental investigation of near-infrared light propagation in a model of the adult head,”
Appl. Opt. , 36
(1), 21
–31
(1997). Google Scholar
M. A. Franceschini
,
S. Fantini
,
L. A. Paunescu
,
J. S. Maier
, and
E. Gratton
,
“Influence of a superficial layer in the quantitative spectroscopic study of strongly scattering media,”
Appl. Opt. , 37
(31), 7447
–7458
(1998). Google Scholar
S. R. Arridge
,
M. Cope
, and
D. T. Delpy
,
“The theoretical basis for the determination of optical pathlengths in tissue: temporal and frequency analysis,”
Phys. Med. Biol. , 37 1531
–1560
(1992). Google Scholar
|