Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Quantitative interpretations of Visible-NIR reflectance spectra of blood

Open Access Open Access

Abstract

This paper illustrates the implementation of a new theoretical model for rapid quantitative analysis of the Vis-NIR diffuse reflectance spectra of blood cultures. This new model is based on the photon diffusion theory and Mie scattering theory that have been formulated to account for multiple scattering populations and absorptive components. This study stresses the significance of the thorough solution of the scattering and absorption problem in order to accurately resolve for optically relevant parameters of blood culture components. With advantages of being calibration-free and computationally fast, the new model has two basic requirements. First, wavelength-dependent refractive indices of the basic chemical constituents of blood culture components are needed. Second, multi-wavelength measurements or at least the measurements of characteristic wavelengths equal to the degrees of freedom, i.e. number of optically relevant parameters, of blood culture system are required. The blood culture analysis model was tested with a large number of diffuse reflectance spectra of blood culture samples characterized by an extensive range of the relevant parameters.

©2008 Optical Society of America

1. Introduction

A new method for detection of microorganisms in blood culture was recently proposed [1]. The method is founded on the following principles: (1) metabolic activity of microorganisms induces changes in the physical and chemical properties of blood and (2) due to distinct optical properties of blood these changes can be captured with diffuse reflectance spectroscopy [1]. In this publication we present a new model implemented with the proposed method designed for quantitative extraction of blood properties from the diffuse reflectance measurements of blood culture systems. This new model is complete in terms of robust solutions of diffusion and scattering problems and their numerical implementation.

The use of spectroscopy to measure the state of blood and blood-infused tissues has long been appealing for biomedical applications both in vitro and in vivo due to, firstly, the distinctive optical signature of hemoglobin and, secondly, the simple, non-invasive, and nondestructive nature of the measurements [26]. Yet, since blood is a complex fluid, a multitude of blood parameters of both physical (size, shape, and concentration of erythrocytes and other cells (i.e., leukocytes, platelets, etc.)), and chemical (i.e., concentrations of hemoglobin, bilirubin, myoglobin, cytochromes, nucleic acids, etc.) nature and system parameters (media, temperature, pH, etc.) affect the measurements and, thus, quantitative analysis of spectroscopic measurements has been challenging [79].

Numerous approaches (statistical, empirical, numerical, and theoretical) have been undertaken for quantitative analysis of measured spectra. Examples of statistical approaches would be multivariate regression models, i.e. partial least squares or principal component analysis, that employ calibration curves relating the measured spectral amplitudes at characteristic wavelengths or some portions of the spectra to the selected blood parameters [2,3,5,7,9]. Empirical models use arbitrarily derived scaling factors, for example, ratios of spectral amplitudes measured at different wavelengths and/or optical configurations [2,10], to fit desired theoretical approximations [6]. Although these techniques are relatively easy and computationally fast, they are limited by the approximation used, typically require extensive calibration and standardization, and often show poor inter-donor consistency.

Numerical approaches to the solution of Maxwell equations for spectroscopic measurements of blood and tissues such as Monte-Carlo simulations, finite boundary element method, discrete dipole approximation, T-matrix, etc. [6,1114], have been extensively used and have the advantage of providing rigorous solutions that involve no major approximations [6]. However, cumbersome and time-consuming computations are required [5,6,12].

Of the theoretical methods, the photon diffusion models have been known to provide an accurate description of the optical measurements of blood [2,15]. Unfortunately, since in the past the rigorous solution to the photon diffusion problem was considered to be mathematically complex [6], the simplifying empirical approximations were generally used [4,6,1620].

To successfully apply the theory, accurate macroscopic scattering and absorption optical properties are needed. Since these optical properties are the functions of multiple blood parameters, their direct extraction from the measurements, for example, with inverse Monte- Carlo simulations [12,21], have produced a wide range of values [16,2123]. For the same reason simplified empirical approximations for scattering [4,6,1820] and/or neglecting the effect of the particle size on absorption [1922] may lead to erroneous results [2].

The macroscopic optical properties of blood can be computed with a choice of scattering theory, i.e., Lorentz-Mie, Rayleigh-Gans, Fraunhofer, or anomalous diffraction [2326], or a numerical technique [11,14,2630]. The major advantage of the numerical techniques is the capability of rigorous computation of optical properties of realistic non-spherical particle shapes such as the biconcave disk of erythrocytes [28]; yet due to computational complexities they are not yet applicable for the routine analysis of blood samples. Of the scattering theories, Mie theory provides an exact solution for the spherical particle of the size range of blood components [15,24,31]. Since the general solution to Mie theory is only for spheres, its applicability to non-spherical particles, such as erythrocytes, has been criticized [14,28,3234]. However, it can be successfully used for non-spherical randomly oriented cells, as those in agitated blood cultures, when size rather than shape determines light scattering [11,27].

The approach to the deconvolution, i.e. quantitative analysis, of diffuse reflectance measurements in terms of optically relevant blood parameters presented in this paper is based on the thorough wavelength dependent solutions to the photon diffusion problem and Mie scattering theory. We show that the accurate implementation of these theories succeeds in the extraction of realistic estimates of the parameters from the diffuse reflectance spectra of blood cultures. Computations are nearly instantaneous enabling the possibility of real-time continuous monitoring applications. Additionally, our approach requires only the knowledge of the refractive indices of the basic constituents of the samples and, therefore, no calibration or standardization is needed as compared to the empirical and statistical methods.

2. Method

2.1 Measurement

Measurements of diffuse reflectance were conducted with a USB-4000 spectrometer characterized by linear CCD detector, signal to noise ratio of 300:1, and 50 µm aperture, equipped with either standard or angled reflectance QR400-Vis-NIR probe each consisting of one 0.2 mm radius read fiber surrounded by six 0.2 mm radius illumination fibers. The probes’ source-detector distance was 0.4 mm and numerical aperture was 0.22. LS-1 or HL-2000 Tungsten Halogen lamps with dynamic range of 360-2000 nm were used as isotropic light sources. The spectrometer, probes, and light sources were purchased from Ocean Optics Inc. (Dunedin, FL). The optical readings were taken by placing the surface of a reflectance probe at 45° angle to the surface of a blood culture vial. Spectra were acquired with the integration time of 150 msec. The spectral resolution of the raw data was 0.33 nm and was interpolated to 1 nm for data analysis. For the measurements reported in this paper we used standard plastic BacT/ALERT blood culture vials [35] purchased from bioMerieux Inc. (Hazelwood, MO) containing 40 ml of growth media. Approximately 8–10 ml of blood was added for the experiments. Blood samples from healthy individuals were provided by Florida Blood Services (St. Petersburg, FL). Some measured spectra presented in this paper were recorded from blood cultures contaminated with microorganisms [36]. In all cases blood cultures were incubated at 37°C. The vials were continuously agitated to ensure that the solutions were well mixed. The diffuse reflectance R(λ) was calculated by scaling the measured reflected intensity of a sample Isample to that of MgO2 (assuming it to be 100%) and subtracting the dark intensity Idark measured when there was no light according to Eq. (1).

R(λ)=IsampleIdarkIMgO2Idark

2.2 Photon diffusion theory

The described measurement configuration satisfied the assumptions under which the photon diffusion theory provides a robust mathematical description of measured diffuse reflectance. These assumptions are that, first, light is scattered multiple times, second, the scattering is larger than absorption, and, third, the depth of the measured sample is larger than the mean photon free path [2,15,16,37,38].

The detailed derivation of the solution to photon diffusion problem has been previously published [15,37]. The resulting solution expresses the reflected intensity R(λ) as a non-linear function of the measurement geometry (detector aperture rx, incident beam radius a, sourcedetector separation distance b, and depth of the sample d) and optical properties (macroscopic reduced scattering µst(λ) and absorption µa(λ) cross-sections) of a measured sample:

R(λ)=rx2r22r12(Rdiff(r2)Rdiff(r1))

where diffuse reflectance intensity of arbitrary radius r, Rdiff, is defined as

Rdiff(r)=2μst(λ)a2n=1knzncosγnNnζn2(kn2+μt2(λ)){(a22)raI1(ζna)K1(ζnr);r>a(r22)raK1(ζna)·I1(ζnr);r<a}

The variables of the Eqs. (2) and (3) are defined in Table 1.

2.3 Formulation of the optical properties

The macroscopic absorption µa(λ) and scattering µs(λ) cross-sections describe the combined optical properties of multiple scatters and chromophores present in a given sample. The first term of Eq. (4) defining total macroscopic absorption cross-section of the sample describes its dependence on single-particle absorption cross-sections Cabs and particle number density N(p) of N populations of particles and the second term incorporates the absorption contributions from M chromophoric compounds with [Ci] being concentration of ith compound dissolved in the media of the sample. In Eq. (5) the total macroscopic scattering cross-section is a function of single-particle scattering cross-sections Csca, particle number density N(p) and the mean volume V(p) of the Jth particle population [25,39]. We use the multi-scattering approximation suggested by [39] in a form of non-linear dependence of the macroscopic scattering cross-section on the particle concentration. The scattering component enters the reflectance calculations (Eq.(3)) through the reduced macroscopic scattering cross-section µst(λ), which results from the correction of the scattering cross-section µs(λ) for anisotropy of the scattering <µ> as shown in Eq. (6).

μa(λ)=J=1NN(p)JCabsJ+i=1M4πki(λ)λ[Ci]
μs(λ)=J=1NN(p)J(1N(p)JV(p)J)CscaJ
μst(λ)=J=1NN(p)J(1N(P)JV(P)J)CscaJ(1<μJ>)
Tables Icon

Table 1. Definition of parameters used in Eqs. (2 and 3).

The single-particle cross-sections, Cabs and Csca, calculated for each particle population were products of the corresponding absorption Qabs and scattering Qsca efficiencies and average scattering area S of a given particle population (Eq. (7)). The efficiencies were computed as shown in Eqs. (8-10) in terms of the particle size parameter α and the particle relative refractive index m(λ) implicit in the coefficients an and bn, the Mie theory formulas for which can be found in [40]. The asymmetry factor <µ> as formulated in Eq. (11) was also a function of particle size and relative refractive index. The scattering area and size parameter of each particle population were calculated under the assumption of homogeneous spherical scatters characterized by mean volume V(p) (Eq. (12)).

Cabs=Qabs(m(λ),V(p))·SCsca=Qsca(m(λ),V(p))·S
Qsca=2α2n=1(2n+1){an2+bn2}
Qext=2α2n=1(2n+1){Re(an+bn)}
Qabs=QextQsca
<μ>=4α2·Qscan=2[n21n{Re(an1a¯n+bn1b¯n)}+2n1n(n1){Re(an1b¯n1)}]
S=(6π2V(p))234α=(6π2V(p))13λ

The Mie scattering algorithm used was constructed following the general computer code of [41] and includes the refractive index increment computations [42]. The computer code has been extensively tested and compared against well-tested Mie codes available [43,44] over the relevant wavelength, size, and optical property sets of interest. The algorithm has been evaluated previously through quantitative analysis of transmission measurements of biological particulates such as blood components and a variety of microorganisms [45, 46].

The effect of multi-component chemical composition for each particle population was introduced through the additive properties of the complex refractive index m(λ) [46-48]. For a particle consisting of M components the real n(λ) and imaginary k(λ) parts of the particle complex refractive index were calculated as weighted sums of the contributions ωi of each component [4648]:

m(λ)=n(λ)+ik(λ)n0(λ)n(λ)=i=1Mωini(λ)k(λ)=i=1Mωiki(λ)

For example, for erythrocytes three chromophoric components (oxy-, deoxy-, and methemoglobin) were considered (Fig. 1), for each the weighted contribution ωi was a product of its relative fraction fi of the total hemoglobin and hemoglobin concentration within erythrocytes [Hb](e). Distinct characteristic spectral signatures of these forms of hemoglobin can be appreciated with Fig. 1. The refractive indices of water n0(λ) and hemoglobin forms were obtained from the published data [49,50].

 figure: Fig. 1.

Fig. 1. The imaginary parts, k(λ), of refractive indices of oxyhemoglobin (HbO2), deoxyhemoglobin (Hb), and methemoglobin (HbFe(III)).

Download Full Size | PDF

2.4 Analysis of data

The combined solution to photon diffusion approximation and Mie theory constituted the fundamental model that was applied to the diffuse reflectance measurements. With this model the diffuse reflectance spectra were deconvoluted in terms of optically relevant blood culture parameters. In the context of this paper deconvolution (or mathematical interpretation) implies an algorithm-based reconstruction of the measured signal with parameterized mathematical functions by which the values of the parameters are estimated. In practical terms the theoretically predicted spectra were computed as functions of blood culture parameters and then fitted to the corresponding measured spectra in the 500–850 nm wavelength range in the least squares sense. A Nelder-Meade downhill simplex optimization algorithm was used for this purpose [51,52].

The optically relevant blood culture parameters included number density N(p) of particle populations present (erythrocytes and other cells), mean particle volume V(p), the products of these parameters being the volume fraction Vf(p) of a given particle population. Also, total amount of hemoglobin [Hb]T, hemoglobin concentration within erythrocytes [Hb](e), concentration of free hemoglobin in media [CHb], fractions of three hemoglobin forms (oxyf(HbO2), deoxy- f(Hb), and met- f(HbFe(III)), and volume fraction of lysed erythrocytes Vf(e)lysed. The mass balance statements relating some of these parameters are shown in Eqs. (14–16). Firstly, since only three forms of hemoglobin were considered, their fractions always added up to one (Eq. (14)). Second, the total amount of hemoglobin [Hb]T in a given sample was the sum of total hemoglobin within erythrocytes, i.e. a product of the erythrocyte hemoglobin concentration [Hb](e) and volume fraction of erythrocyte population V(e), and total hemoglobin dispersed in the media, i.e. a product of free hemoglobin concentration [CHb] and sample volume V (Eq. (15)). Thirdly, concentration of hemoglobin in the media [CHb] was a function of the fraction of lysed erythrocytes Vf(e)lysed as shown in Eq. (16). The mass balance constraints were implemented using variable transformation techniques [53].

f(HbO2)+f(Hb)+f(HbFe(III))=1
[Hb]T=[Hb](e)·Vf(e)+[CHb]·V
[CHb]=1V[Hb](e)·Vf(e)·Vf(e)lysed

Given these relationships, the independently estimated parameters were number density N(p), mean particle volume V(p), hemoglobin concentration within erythrocytes [Hb](e), fractions of hemoglobin forms fi, and media hemoglobin concentration [CHb]. The erythrocytes’ and media hemoglobin concentrations and the fractions of hemoglobin forms entered the model through the formulation of the complex refractive index of the erythrocyte population in Eq. (13) and media absorption, i.e. second term of Eq. (4). The mean particle volume V (p) was included in the formulation of both single-particle and macroscopic cross-sections (Eqs. (5–12)). The number densities N (p) of particle populations were modeled with Eqs. (4–6). It is evident from Eqs. (2–13) that the model is highly non-linear in terms of the parameters being estimated and as such the uniqueness of the parameters cannot be ensured. Nevertheless, extensive numerical sensitivity analyses demonstrated that the model is indeed sensitive to the parameters listed. Furthermore, by constraining the parameter search to the range of physiologically feasible values, convergence to a unique solution from a wide variety of initial conditions is always achieved. The average errors of the estimation of the blood culture parameters were ±0.1*108 cell ml-1 for number density N (p), ±0.5 µm3 for mean particle volume V(p), ±0.0003 for particle volume fraction Vf(p), ±0.005 for hemoglobin fraction within erythrocytes [Hb](e), ±0.0001 g ml-1 for concentration of free hemoglobin in media [CHb], ±0.0015 for the fractions of hemoglobin forms, and ±0.01 for volume fraction of lysed erythrocytes Vf(e)lysed.

3. Results

Typical sets of wavelength-dependent diffuse reflectance spectra of blood cultures are shown in Fig. 2. These spectral sets were recorded over time from two independent blood cultures of different experimental set-ups. In brief, the spectra plotted in Panel A of Fig. 2 were recorded from a blood culture contaminated with E. coli. The spectra were collected for 12 hours with 2 min interval between the measurements. These spectra show a transition from oxyhemoglobin with characteristic reflectance peak at 660 nm to deoxyhemoglobin with characteristic reflectance peaks at 725 nm and 810 nm over the course of the experiment. The spectra in Panel B of Fig. 2 were recorded from a non-contaminated blood culture for over 120 hours with 17 min interval between the measurements. In this example, a transition from oxyhemoglobin to methemoglobin with characteristic peaks at 630 nm and 690 nm can be observed. These changes in the spectral features were results of metabolic activity of contaminating bacteria (Panel A) and aging process of erythrocytes (Panel B) as discussed in detail in [1].

 figure: Fig. 2.

Fig. 2. Panel A: diffuse reflectance spectra collected during a blood culture experiments with E. coli. Panel B: diffuse reflectance spectra collected during a blood culture experiment with no contaminants.

Download Full Size | PDF

In blood culture experiments erythrocytes are the primary scattering elements, as they are the dominant blood components, and they typically have volume fractions of <0.1 in blood cultures. Hemoglobin is the dominant chromophoric component of erythrocytes and, therefore, determined the major features of the measured spectral curves. Therefore, the evident variability in the recorded spectral shapes of Fig. 2 largely resulted from the changes in the chemical composition of hemoglobin over the measurement time periods. The signature absorption bands of hemoglobin forms, for example, 542 nm and 577 nm of oxyhemoglobin, 558 nm and 760 nm of deoxyhemoglobin, and 627 nm of methemoglobin (Fig. 1), appear as depressions in the diffuse reflectance spectra (Fig. 2). Notice that the characteristic absorption signatures of the hemoglobin forms are enhanced in diffuse reflectance spectra through non-trivial dependence of the reflectance on both scattering and absorption (Eqs. (2–3)). By means of such dependence the differences in the characteristic spectroscopic features of hemoglobin forms enable accurate quantitative assessment of hemoglobin composition. Furthermore, the physiological parameters of the measured samples such as number density and mean volume of erythrocytes (and other cells) can be derived from the measured spectra. In the following discussion we illustrate these claims through a series of characteristic examples. The contributions of the different physiological parameters of the components of blood cultures to the diffuse reflectance signals will be addressed individually.

The first set of the characteristic examples of the deconvolution of the measured reflectance spectra of blood cultures is shown in Fig. 3. It addresses the model’ ability to estimate the fractions of the multiple hemoglobin forms present in a given blood culture sample. The examples shown in Fig. 3 are individual spectra from three blood culture experiments. The spectra in Panels A and B were recorded after inoculation of different blood samples to aerobic and anaerobic blood culture vials, respectively. The spectrum in Panel C was recorded after long-term incubation of an aerobic blood culture. Each of the examples had distinct hemoglobin composition; yet all three main forms of hemoglobin contributed to the measured reflectance signals of all spectra. The estimates of the hemoglobin fractions obtained via deconvolution of the measured spectra are listed in Table 2. The estimates of other blood parameters such as the volume fraction and mean volume of erythrocytes are also given (Table 2). Each subplot of Fig. 3 displays a comparison between a measured and the corresponding theoretically predicted reflectance spectra along with the residuals of the fit.

The following examples address the model’s capability to quantify the optically relevant parameters of blood cultures such as the volume fraction and mean volume of the scatters. These properties determine the magnitudes of macroscopic cross-sections (Eqs. (4–6)) and thus notably contribute to the reflectance signal. Figure 4 illustrates typically observed differences in the shape and magnitude of the reflectance spectra of blood cultures resulting from differences in the mean volume of the scattering elements. In this example the scatter number density and chemical composition of the compared blood cultures are similar (Table 3). However, since the volume fraction of the scatters is directly proportional to their mean volume, this parameter also varied between the measured samples. Figure 4 shows that the changes in both the mean volume and the volume fraction of scattering elements affect the relative amplitude and character of the reflectance spectra.

 figure: Fig. 3.

Fig. 3. The measured (Rmeas) and predicted (Rcalc) diffuse reflectance spectra and residuals (σres) of three different blood cultures.

Download Full Size | PDF

Tables Icon

Table 2. The relevant parameters of blood culture components obtained from the deconvolution of the diffuse reflectance spectra presented in Fig. 2.

 figure: Fig. 4.

Fig. 4. The measured (Rmeas) and predicted (Rcalc) diffuse reflectance spectra and residuals (σres) for blood culture samples characterized by the mean erythrocyte volume of 82 µm3 (A) and the mean erythrocyte volume of 99 µm3 (B).

Download Full Size | PDF

Tables Icon

Table 3. The relevant parameters of blood culture components obtained from the deconvolution of the diffuse reflectance spectra presented in Fig. 4.

Figure 5 demonstrates the spectral differences between two samples characterized by similar chemical composition and the equal volume fractions of scatters but by distinct mean volumes of the scatters (Table 4). In particular, spectrum B was recorded of a blood culture with larger number density of erythrocytes than that of spectrum A. Note, that the effects of the volume fraction and the mean volume of scatters on the reflectance signal are complex through non-linear dependency of macroscopic scattering cross-section on the first parameter (Eq. (5)) and complex formulations of single-particle cross-sections as functions of particle size (Eqs. (7–12))

 figure: Fig. 5.

Fig. 5. The measured (Rmeas) and predicted (Rcalc) diffuse reflectance spectra and residuals (σres) for blood culture samples characterized by the mean erythrocyte volume of 92 µm3 (A) and the mean erythrocyte volume of 77 µm3 (B).

Download Full Size | PDF

Tables Icon

Table 4. The relevant parameters of blood culture components obtained from the deconvolution of the diffuse reflectance spectra presented in Fig. 5.

The previous examples demonstrated good agreement between the measured and reconstructed reflectance spectra and the corresponding deconvolution estimates of blood culture parameters for non-contaminated samples, i.e. containing only typical blood components in culture media. Since erythrocytes are the predominant particle population (>99.99% of all particles by volume) of normal healthy blood, it was reasonable to approximate that those samples had only one particle population. With erythrocytes having the number density of 108 cells ml-1 and therefore occupying up to 10% of liquid phase of a typical blood culture, for any other population to have a measurable contribution to the reflectance signal, comparable particle number density and/or volume fraction are needed.

Although microbial contaminants growing in blood cultures typically have small size (for example, characteristic average diameter of bacteria is about 0.5-2 µm [45,47] and yeast cells are typically 3–10 µm in diameter [47]), they usually reach the number densities of 107 -109 cells ml-1 in positive blood cultures [35,36] and thus can alter the reflectance signal. This outcome is illustrated with Fig. 6, where three spectra recorded at different time steps of the same blood culture experiment with two contaminants (bacterium Staphylococcus epidermidis and yeast Candida glabrata) are presented. Of these spectra, A was recorder prior build-up of microbial biomass, B was recorded when S. epidermidis reached number density of >108 cells ml-1, and C was recorded when the yeast cells reached similar optically relevant concentration. The experiment started with low inocula of both contaminants and because of the significant dissimilarity in the growth rates of these two microorganisms leading to the difference in time needed to reach the optically relevant cell number densities, it was possible to capture the individual contributions of these microorganisms to the reflectance signal. In Fig. 6, the spectral change from spectrum A to spectrum B resulted from S. epidermidis reaching a calculated number density of 0.8·109 cells ml-1 (Table 5). Since the cells of S. epidermidis typically form clusters [36,47], the bacterial mean cell volume was approximated as 3 µm3 (Table 5). Further, the growth of C. glabrata produced additional spectral changes (i.e., the shift from spectrum B to spectrum C in Fig. 6). With approximated mean cell volume of yeast of 33.5 µm3 cell number density for C. glabrata was estimated to be 7·107 cells ml-1 for spectrum C (Table 5). These results suggest that physiological properties of blood culture components other than erythrocytes can also be captured and quantitatively resolved from diffuse reflectance measurements.

 figure: Fig. 6.

Fig. 6. The measured (Rmeas) and predicted (Rcalc) diffuse reflectance spectra and residuals (σres) for blood culture samples characterized by one particle population (erythrocytes) (A), two particle populations (erythrocytes and bacterial cells) (B), and three particle populations (erythrocytes, bacterial cells, and yeast cells) (C).

Download Full Size | PDF

Tables Icon

Table 5. The relevant parameters of blood culture components obtained from the deconvolution of the diffuse reflectance spectra presented in Fig. 6.

Another important process that may to occur in blood cultures is lysis of erythrocytes [36]. Erythrocyte destruction unavoidably induces the changes in the number density and mean volume of the scatters. Since erythrocytes are the primary particle population in blood cultures, the effect of this process on the reflectance signal is dramatic as shown in Fig. 7. The spectra in Fig. 7 were collected from an anaerobic blood culture contaminated with Clostridium perfringens, an anaerobic hemolytic bacterium, at two time points of the experiment. Spectrum A was recorded prior to hemolysis and spectrum B was recorded later into the experiment when hemolysis occurred as a result of metabolic activity of C. perfringens. The observed spectral differences are the collective outcome of the physical changes in the sample, i.e. in the size, volume fraction, and number density of particulates as erythrocytes rupture into pieces, and chemical changes resulted from the increase in the concentration of free hemoglobin in the media (Table 6). Both physical and chemical perturbations of lysis led to the observed decrease in the magnitude of the reflected signal. Scattering decreased with reduced volume fraction of scattering elements and absorption increased with hemoglobin becoming free in solution. Moreover, the reduction of scattering caused the absorptive spectral features in diffuse reflectance to be less pronounced and produced an apparent flattening of the spectrum (Fig. 7). By including mathematical description of these processes in the formulation of macroscopic cross-sections (Eqs. (4–6)) the measured reflectance spectra of samples containing lysed erythrocytes can be suitably predicted and reasonable quantitative estimates of the properties of sample components can be estimated (Table 6).

 figure: Fig. 7.

Fig. 7. Temporal progression of oxyhemoglobin (A), fraction of deoxyhemoglobin (B), fraction of methemoglobin (C), and the mean cell volume (MCV) of erythrocytes (D) along with the residual sum of squares (E) obtained from deconvolution of diffuse reflectance spectra recorded during an on-line blood culture experiment.

Download Full Size | PDF

Tables Icon

Table 6. The relevant parameters of blood culture components obtained from the deconvolution of the diffuse reflectance spectra presented in Fig. 7.

Further, temporal progression of blood culture parameters estimated from continuously measured reflectance spectra can be evaluated. For example, Fig. 8 shows the temporal changes in the fractions of hemoglobin forms and mean volume of erythrocytes estimated from reflectance measurements of non-contaminated blood culture. The measurements were taken every 15 min for over 120 hrs and each measured spectrum was deconvoluted with the model. Tracking the temporal variability in the estimated parameters allows for quantitative assessment of the processes taking place within the blood culture as described in [1,36,54]. The residual sum of squares plotted in panel E of Fig. 8 indicates good agreement between the measured and theoretically predicted spectra and therefore validates the estimated parameter values.

 figure: Fig. 8.

Fig. 8. The measured (Rmeas) and predicted (Rcalc) diffuse reflectance spectra and residuals (σres) for blood culture samples with intact erythrocytes (A) and partially lysed erythrocytes (B).

Download Full Size | PDF

4. Discussion

The solutions to photon diffusion approximation and Mie scattering theory are well-known and have been used for over several decades, but surprisingly their applications to the spectroscopic measurements of blood have been restricted to one or several selected wavelengths [2,15]. Overall, the possibilities offered by the wavelength-dependent measurements have been explored only in few studies [21,22]. Furthermore, the number of reported simulations of the wavelength-dependent spectroscopic behavior of blood in terms of reflectance is very limited [30,55,56]. The approach presented in this paper demonstrates the advantages offered by the non-linear dependence on wavelength of both real and imaginary parts of hemoglobin refractive index that can be exploited in the context of diffuse reflectance measurements of blood cultures. Complete employment of the model in terms of these nonlinear dependencies is vital for mathematical description of scattering and absorption properties of blood.

In this paper we briefly reported on the capabilities of the new approach designed to quantitatively extract the physiological and chemical parameters of the components of blood cultures from diffuse reflectance measurements. This interpretation model has been applied effectively for the analysis of a multitude of reflectance spectra from blood cultures. The examples of the recorded spectra with corresponding theoretical deconvolutions presented in this paper were intended to illustrate how the variability in the optically relevant properties of blood culture components alter the diffuse reflectance spectra and how these properties can be estimated. We used the residuals of the fit between the measured and predicted spectra as a measure of the quality of theoretical predictions. The excellent agreements between the measurements and theory confirmed the validity of our approach. As defined by Eqs. (2–13) the reflectance spectra are functions of the number density, size, and composition as well as the absorption and scattering properties of the particle populations present in a given sample. In terms of blood cultures, such important variables as number density and mean volume of erythrocytes, concentration of hemoglobin within the erythrocytes and dissolved in media, and fractions of the main forms of hemoglobin can be assessed. From these variables, the erythrocyte volume fraction (hematocrit) and the amount of erythrocyte lysis can be derived. Further, we have shown that the composite effects of multiple particle populations can be assessed from diffuse reflectance spectra. The additional populations could be, for instance, microbial contaminants; when their volume fractions exceed 0.01% they can produce measurable changes to the diffuse reflectance signals. Therefore, the mean cell size, number density, and chemical composition (if chromophoric) of the supplementary particle populations can be estimated. The interpretation model was shown to produce comprehensive estimates of these variables.

While a number of comprehensive approaches for theoretical and/or numerical description of reflectance or backscattering measurements has been developed [1521], they have been restricted to single wavelength measurements and even then their computational complexity often require simplifying approximations. From this perspective, the presented approach for spectral deconvolution has an important advantage of being computationally fast and efficient. Our interpretation model is envisioned to be used for on-line and real-time applications, for example, in clinical blood culture analyses as illustrated in Fig. 8.

5. Conclusion

Diffuse reflectance is simple, non-invasive, and non-destructive type of spectroscopic measurements that by its virtues has been extensively exploited for various applications including those in biomedical field. In this paper we report on a new algorithm for accurate quantitative description of the diffuse reflectance spectra from blood cultures in terms of the relevant parameters of the blood culture components. For this method, the rigorous solutions to photon diffusion approximation and Mie scattering theory have been revisited. Adequate mathematical descriptions of multiple scattering, multiple scattering populations, and multiple chromophores present within the scatters as well as in solution were added. In fact, the success of new interpretation model lies in the utilization of non-linear wavelength dependence of the basic optical properties, i.e., refractive indices, of the chromophoric components present in the measured samples. The displayed solution to the diffuse reflectance problem was found to give appropriate estimates of the blood culture parameters as was illustrated with a scope of characteristic examples. Furthermore, the features of this quantitative theoretical deconvolution model eliminate the need for external calibrations (other than scaling for white intensity) and therefore lead to superior interpretation of the measured spectra.

Acknowledgments

The authors would like to acknowledge the support of Claro Scientific, LLC., Florida Blood Services, and of the College of Marine Science, University of South Florida.

References and links

1. J. M. Smith, Y. M. Serebrennikova, D. E. Huffman, G. F. Leparc, and L. H. García-Rubio, “A new method for the detection of microorganisms in blood cultures. Part I: Theoretical analysis and simulation of blood culture processes,” Can. J. Chem. Eng. 86, 947–959 (2008). [CrossRef]  

2. J. M. Steinke and A. P. Shepherd, “Reflectance measurements of hematocrit and oxyhemoglobin saturation,” Am. J. Physiol. 253, H147–H153 (1987). [PubMed]  

3. S. Zhang, B. R. Soller, S. Kaur, K. Perras, and T. J. vander Salm, “Investigation of noninvasive in vivo blood hematocrit measurement using NIR reflectance spectroscopy and partial least-square regression,” Appl. Spectrosc. 54, 294–299 (2000). [CrossRef]  

4. M. Johns, C. A. Giller, and H. Liu, “Determination of hemoglobin oxygen saturation from turbid media using reflectance spectroscopy with small source-detector separations,” Appl. Spectrosc. 55, 1686–1693 (2001). [CrossRef]  

5. M. Meinke, I. Gersonde, M. Friebel, J. Helfmann, and G. Muller, “Chemometric determination of blood parameters using Visible-Near-Infrared spectra,” Appl. Spectrosc. 59, 826–835 (2005). [CrossRef]   [PubMed]  

6. G. Zonios and A. Dimou, “Modeling diffuse reflectance from semi-infinite turbid media: application to the study of skin optical properties,” Opt. Express 14, 8661–8674 (2006). [CrossRef]   [PubMed]  

7. K.R. Denninghoff, R. A. Chipman, and L. W. Hillman, “Oxyhemoglobin saturation measurements by green spectral shift,” Opt. Lett. 31, 924–926 (2006). [CrossRef]   [PubMed]  

8. J. M. Steinke and A. P. Shepherd, “Effects of temperature on optical absorbance spectra of oxy-, carboxy-, and deoxyhemoglobin,” Clin. Chem. 38, 1360–1364 (1992). [PubMed]  

9. K. A. Schenkman and W. A. Ciesielski, “Improved myoglobin saturation measurement made by partial leastsquare analysis of optical reflectance spectra,” Appl. Spectrosc. 56, 1215–1221 (2002). [CrossRef]  

10. A. P. Shepherd, V. T. Randal, J. M. Steinke, and J. L. Schmalzel, “An oximeter for measuring hemoglobin concentration and oxygen content,” Am. J. Physiol. 257, H1705–H1711 (1989). [PubMed]  

11. A. M. K. Enejder, J. Swartling, P. Aruna, and S. Andersson-Engels, “Influence of cell shape and aggregate formation on the optical properties of flowing whole blood,” Appl. Opt. 42, 1384–1394 (2003). [CrossRef]   [PubMed]  

12. M. S. Patterson, B. C. Wilson, and D. R. Wyman, “The propagation of optical radiation in tissue I. Models of radiation transport and their application,” Las. Med. Sci. 6, 155–168 (1991). [CrossRef]  

13. I. Turcu, “Effective phase function for light scattered by blood,” Appl. Opt. 45, 639–647 (2006). [CrossRef]   [PubMed]  

14. S. V. Tsinopoulos, E. J. Sellountos, and D. Polyzos, “Light scattering by aggregated red blood cells,” Appl. Opt. 41, 1408–1417 (2002). [CrossRef]   [PubMed]  

15. L. Reynolds, C. Johnson, and A. Ishimaru, “Diffuse reflectance from a finite blood medium: applications to the modeling of fiber optic catheters,” Appl. Opt. 15, 2059–2067 (1976). [CrossRef]   [PubMed]  

16. R. A. J. Groenhuis, H. A. Ferwerda, and J. J. Ten Bosch, “Scattering and absorption of turbid materials determined from reflection measurements. 1: Theory,” Appl. Opt. 22, 2456–2462 (1983). [CrossRef]   [PubMed]  

17. T. J. Farrell, M. S. Patterson, and B. Wilson, “A diffusion theory model of spatially resolved, steady state diffuse reflectance for noninvasive determination of tissue optical properties in vivo,” Med. Phys. 19, 879–888 (1992). [CrossRef]   [PubMed]  

18. R. Graaff, J. G. Aarnoudse, J. P. Zijp, M. A. Sloot, F. F. M. de Mul, J. Greve, and M. H. Koelink, “Reduced light-scattering properties for mixtures of spherical particles: a simple approximation derived from Mie calculations,” Appl. Opt. 31, 1370–1376 (1992). [CrossRef]   [PubMed]  

19. G. Zonios, L. T. Perelman, V. Backman, R. Manoharan, M. Fitzmaurice, J. Van Dam, and M. S. Feld, “Diffuse reflectance spectroscopy of human adenomatous colon polyps in vivo,” Appl. Opt. 38, 6628–6637 (1999). [CrossRef]  

20. G. Zonios, J. Bykowski, and N. Kollias, “Skin melanin, hemoglobin, and light scattering properties can be quantitatively assessed in vivo using diffuse reflectance spectroscopy,” J. Invest. Dermatol. 117, 1452–1457 (2001). [CrossRef]  

21. A. Roggan, M. Friebel, K. Dorschel, A. Hahn, and G. Muller, “Optical properties of circulating human blood in the wavelength range 400–2500 nm,” J. Biomed. Opt. 4, 36–46 (1999). [CrossRef]  

22. M. Meinke, G. Muller, J. Helfmann, and M. Friebel, “Empirical model functions to calculate hematocritdependent optical properties of human blood,” Appl. Opt. 46, 1742–1753 (2007). [CrossRef]   [PubMed]  

23. D. J. Faber, C. G. Aalders, E. G. Mik, B. A. Hooper, M. J. C. van Gemert, and T. G. van Leeuwen, “Oxygen saturation dependent absorption and scattering of blood,” Phys. Rev. Lett. 93, 028102 (2004). [CrossRef]   [PubMed]  

24. D. H. Tycko, M. H. Metz, E. A. Epstein, and A. Grinbaum, “Flow-cytometric light scattering measurement of red blood cell volume and hemoglobin concentration,” Appl. Opt. 24, 1355–1364 (1985). [CrossRef]   [PubMed]  

25. G. J. Streekstra, A. G. Hoekstra, E. J. Nijhof, and R. M. Heethaar, “Light scattering by red blood cells in ektacytometry: Fraunhofer versus anomalous diffraction,” Appl. Opt. 32, 2266–2272 (1993). [CrossRef]   [PubMed]  

26. M. Hammer, D. Schweitzer, B. Michel, E. Thamm, and A. Kolb, “Single scattering by red blood cells,” Appl. Opt. 37, 7410–7418 (1998). [CrossRef]  

27. A. M. K. Nilsson, P. Alsholm, A. Karlsson, and S. Andersson-Engels, “T-matrix computations of light scattering by red blood cells,” Appl. Opt. 37, 2735–2748 (1998). [CrossRef]  

28. S. V. Tsinopoulos and D. Polyzos, “Scattering of He-Ne laser light by an average-sized red blood cell,” Appl. Opt. 38, 5499–5510 (1999). [CrossRef]  

29. J. He, A. Karlsson, J. Swartling, and S. Andersson-Engels, “Light scattering by multiple red blood cells,” J. Opt. Soc. Am.A 21, 1953–1961 (2004). [CrossRef]  

30. A. Karlsson, J. He, J. Swartling, and S. Andersson -Engels, “Numerical simulations of light scattering by red blood cells,” IEEE Trans. Biomed. Eng. 52, 13–18 (2005). [CrossRef]   [PubMed]  

31. J. M. Steinke and A. P. Shepherd, “Diffusion model of the optical absorbance of whole blood,” J. Opt. Soc. Am. A 5, 813–822 (1988). [CrossRef]   [PubMed]  

32. J. M. Steinke and A. P. Shepherd, “Comparison of Mie theory and the light scattering of red blood cells,” Appl. Opt. 27, 4027–4033 (1988). [CrossRef]   [PubMed]  

33. R. Drezek, A. Dunn, and R. Richards-Kortum, “Light scattering from cells: finite-difference time-domain simulations and goniometric measurements,” Appl. Opt. 38, 3651–3661 (1999). [CrossRef]  

34. R. Drezek, A. Dunn, and R. Richards-Kortum, “A pulsed finite-difference time-domain (FDTD) method for calculating light scattering from biological cells over broad wavelength ranges,” Opt. Express 6, 147–157 (2000). [CrossRef]   [PubMed]  

35. T. C. Thorpe, M. L. Wilson, J. E. Turner, J. L. DiGuiseppi, M. Willert, S. Mirrett, and L. B. Reller, “BacT/Alert: an automated colorimetric microbial detection system,” J. Clin. Microbiol. 28, 1608–1612 (1990). [PubMed]  

36. D. E. Huffman, Y. M. Serebrennikova, J. M. Smith, G. F. Leparc, and L. H. García-Rubio, “A new method for the detection of microorganisms in blood cultures: Application of the quantitative interpretation to aerobic blood cultures,” J. Biomed. Opt. (to be published). [PubMed]  

37. A. Ishimaru, Wave propagation and scattering in random media (Academic Press, 1978).

38. R. Aronson and N. Corngold, “Photon diffusion coefficient in an absorbing medium,” J. Opt. Soc. Am. A 16, 1066–1071 (1999). [CrossRef]  

39. V. Twersky, “Absorption and multiple scattering by biological suspensions,” J. Opt. Soc. Am. 60, 1084–1093 (1970). [CrossRef]   [PubMed]  

40. M. Kerker, The Scattering of Light and other electromagnetic radiation (Academic Press, 1969).

41. C. Bohren and D. R. Huffman, Absorption and Scattering by Small Particles (John Wiley & Sons, 1983).

42. B. H. Zimm and W. B. Dandliker, “Theory of light scattering and refractive index of solutions of large colloidal particles,” J. Phys. Chem. 58, 644–648 (1954). [CrossRef]  

43. W. J. Lentz, “Generating Bessel functions in Mie scattering calculations using continued fractions,” Appl. Opt. 15, 668–671 (1976). [CrossRef]   [PubMed]  

44. W. J. Wiscombe, Mie Scattering Calculations: Advances in Technique and Fast Vector-Speed Computer Codes (National Technical Information Center, 1983), PB 301388.

45. Y. Mattley, G. Leparc, R. Potter, and L. García-Rubio, “Light scattering and absorption model for the quantitative interpretation of human blood platelet spectral data,” Photochem. Photobiol. 71, 610–619 (2000). [CrossRef]   [PubMed]  

46. C. E. Alupoaei and L. H. García-Rubio, “Growth behavior of microorganisms using UV-Vis spectroscopy: Escherichia coli,” Biotech. Bioeng. 86, 163–167 (2004). [CrossRef]  

47. C. E. Alupoaei, J. A. Olivares, and L. H. García-Rubio, “Quantitative spectroscopy analysis of prokaryotic cells: vegetative cells and spores,” Biosens. Bioelectron. 19, 893–908 (2004). [CrossRef]   [PubMed]  

48. C. E. Alupoaei and L. H. García-Rubio, “An interpretation model for the UV-VIS spectra of microorganisms,” Chem. Eng. Comm. 192, 198–218 (2005). [CrossRef]  

49. I. Thormählen, J. Straub, and U. Grigull, “Refractive index of water and its dependence on wavelength temperature, and density,” J. Phys. Chem. Ref. Data 14, 933–946 (1985). [CrossRef]  

50. W. G. Zijistra, A. Buursma, and W. P. Meeuwsen-van der Roest, “Absorption Spectra of Human Fetal and Adult Oxyhemoglobin, De-Oxyhemoglobin, Carboxyhemoglobin, and Methemoglobin,” Clin. Chem. 37, 1633–1638 (1991).

51. J.A. Nelder and R. Mead, “A simplex method for function minimization,” Comp. J. 7, 308–313 (1965).

52. J. L. Kuester and J. H. Mize, Optimization Techniques with Fortran (McGraw-Hill, New York, 1973).

53. M. J. Box, “A comparison of several current optimization methods, and the use of transformations in constrained problems,” Comp. J. 9, 67–77 (1966).

54. Y. M. Serebrennikova, J. M. Smith, D. E. Huffman, L. R. García-Rubio, and Claro Scientific, LLC., 10100 Dr. Martin Luther King Jr. St., N., St. Petersburg, FL, 33716, are preparing a manuscript to be called “A new method for the detection of yeast in blood cultures”.

55. S. J. Matcher, M. Cope, and D. T. Deply, “In vivo measurements of the wavelength dependence of tissuescattering coefficients between 760 and 900 nm measured with time-resolved spectroscopy,” Appl. Opt. 36, 386–396 (1997). [CrossRef]   [PubMed]  

56. V. V. Tuchin, D. M. Zhestkov, A. N. Bashkatov, and E. A. Genina, “Theoretical study of immersion optical clearing of blood in vessels at local hemolysis,” Opt. Express 12, 2966–2971 (2004). [CrossRef]   [PubMed]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (8)

Fig. 1.
Fig. 1. The imaginary parts, k(λ), of refractive indices of oxyhemoglobin (HbO2), deoxyhemoglobin (Hb), and methemoglobin (HbFe(III)).
Fig. 2.
Fig. 2. Panel A: diffuse reflectance spectra collected during a blood culture experiments with E. coli. Panel B: diffuse reflectance spectra collected during a blood culture experiment with no contaminants.
Fig. 3.
Fig. 3. The measured (Rmeas) and predicted (Rcalc) diffuse reflectance spectra and residuals (σres) of three different blood cultures.
Fig. 4.
Fig. 4. The measured (Rmeas) and predicted (Rcalc) diffuse reflectance spectra and residuals (σres) for blood culture samples characterized by the mean erythrocyte volume of 82 µm3 (A) and the mean erythrocyte volume of 99 µm3 (B).
Fig. 5.
Fig. 5. The measured (Rmeas) and predicted (Rcalc) diffuse reflectance spectra and residuals (σres) for blood culture samples characterized by the mean erythrocyte volume of 92 µm3 (A) and the mean erythrocyte volume of 77 µm3 (B).
Fig. 6.
Fig. 6. The measured (Rmeas) and predicted (Rcalc) diffuse reflectance spectra and residuals (σres) for blood culture samples characterized by one particle population (erythrocytes) (A), two particle populations (erythrocytes and bacterial cells) (B), and three particle populations (erythrocytes, bacterial cells, and yeast cells) (C).
Fig. 7.
Fig. 7. Temporal progression of oxyhemoglobin (A), fraction of deoxyhemoglobin (B), fraction of methemoglobin (C), and the mean cell volume (MCV) of erythrocytes (D) along with the residual sum of squares (E) obtained from deconvolution of diffuse reflectance spectra recorded during an on-line blood culture experiment.
Fig. 8.
Fig. 8. The measured (Rmeas) and predicted (Rcalc) diffuse reflectance spectra and residuals (σres) for blood culture samples with intact erythrocytes (A) and partially lysed erythrocytes (B).

Tables (6)

Tables Icon

Table 1. Definition of parameters used in Eqs. (2 and 3).

Tables Icon

Table 2. The relevant parameters of blood culture components obtained from the deconvolution of the diffuse reflectance spectra presented in Fig. 2.

Tables Icon

Table 3. The relevant parameters of blood culture components obtained from the deconvolution of the diffuse reflectance spectra presented in Fig. 4.

Tables Icon

Table 4. The relevant parameters of blood culture components obtained from the deconvolution of the diffuse reflectance spectra presented in Fig. 5.

Tables Icon

Table 5. The relevant parameters of blood culture components obtained from the deconvolution of the diffuse reflectance spectra presented in Fig. 6.

Tables Icon

Table 6. The relevant parameters of blood culture components obtained from the deconvolution of the diffuse reflectance spectra presented in Fig. 7.

Equations (16)

Equations on this page are rendered with MathJax. Learn more.

R ( λ ) = I sample I dark I MgO 2 I dark
R ( λ ) = r x 2 r 2 2 r 1 2 ( R diff ( r 2 ) R diff ( r 1 ) )
R diff ( r ) = 2 μ st ( λ ) a 2 n = 1 k n z n cos γ n N n ζ n 2 ( k n 2 + μ t 2 ( λ ) ) { ( a 2 2 ) ra I 1 ( ζ n a ) K 1 ( ζ n r ) ; r > a ( r 2 2 ) ra K 1 ( ζ n a ) · I 1 ( ζ n r ) ; r < a }
μ a ( λ ) = J = 1 N N ( p ) J C absJ + i = 1 M 4 π k i ( λ ) λ [ C i ]
μ s ( λ ) = J = 1 N N ( p ) J ( 1 N ( p ) J V ( p ) J ) C scaJ
μ st ( λ ) = J = 1 N N ( p ) J ( 1 N ( P ) J V ( P ) J ) C scaJ ( 1 < μ J > )
C abs = Q abs ( m ( λ ) , V ( p ) ) · S C sca = Q sca ( m ( λ ) , V ( p ) ) · S
Q sca = 2 α 2 n = 1 ( 2 n + 1 ) { a n 2 + b n 2 }
Q ext = 2 α 2 n = 1 ( 2 n + 1 ) { Re ( a n + b n ) }
Q abs = Q ext Q sca
< μ > = 4 α 2 · Q sca n = 2 [ n 2 1 n { Re ( a n 1 a ¯ n + b n 1 b ¯ n ) } + 2 n 1 n ( n 1 ) { Re ( a n 1 b ¯ n 1 ) } ]
S = ( 6 π 2 V ( p ) ) 2 3 4 α = ( 6 π 2 V ( p ) ) 1 3 λ
m ( λ ) = n ( λ ) + ik ( λ ) n 0 ( λ ) n ( λ ) = i = 1 M ω i n i ( λ ) k ( λ ) = i = 1 M ω i k i ( λ )
f ( HbO 2 ) + f ( Hb ) + f ( HbFe ( III ) ) = 1
[ Hb ] T = [ Hb ] ( e ) · V f ( e ) + [ C Hb ] · V
[ C Hb ] = 1 V [ Hb ] ( e ) · V f ( e ) · V f ( e ) lysed
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.