18 April 2012 In vivo measurement of the shape of the tissue-refractive-index correlation function and its applicationto detection of colorectal field carcinogenesis
Author Affiliations +
J. of Biomedical Optics, 17(4), 047005 (2012). doi:10.1117/1.JBO.17.4.047005
Polarization-gated spectroscopy is an established method to depth-selectively interrogate the structural properties of biological tissue. We employ this method in vivo in the azoxymethane (AOM)-treated rat model to monitor the morphological changes that occur in the field of a tumor during early carcinogenesis. The results demonstrate a statistically significant change in the shape of the refractive-index correlation function for AOM-treated rats versus saline-treated controls. Since refractive index is linearly proportional to mass density, these refractive-index changes can be directly linked to alterations in the spatial distribution patterns of macromolecular density. Furthermore, we found that alterations in the shape of the refractive-index correlation function shape were an indicator of both present and future risk of tumor development. These results suggest that noninvasive measurement of the shape of the refractive-index correlation function could be a promising marker of early cancer development.
Gomes, Ruderman, DelaCruz, Wali, Roy, and Backman: In vivo measurement of the shape of the tissue-refractive-index correlation function and its applicationto detection of colorectal field carcinogenesis



Colorectal cancer (CRC) is the third leading cause of cancer deaths in Americans despite being highly treatable in its early stages.1 This reflects both the lack of conspicuous symptoms in CRC’s early stages and the lack of compliance with CRC screening guidelines. Colonoscopy is the current gold standard for CRC screening of the at-risk population because of its outstanding accuracy in detecting advanced neoplastic lesions and its simultaneous ability to reduce the CRC incidence by endoscopic removal of these lesions.2,3 However, due to concerns of cost, invasiveness, discomfort with bowel cleansing, and potential complications, the majority of eligible screening participants forgo the colonoscopy exam. Even if compliance could be improved, it would impose heavy financial and resource burdens on the healthcare system.4 Consequently, there have been recent attempts to find more cost-effective and/or less invasive methods of CRC screening. Among the best known are fecal occult blood testing, fecal DNA analysis, sigmoidoscopy, and computed-tomography (CT) colography. Issues of diagnostic sensitivity5,6 and ionizing radiation in the case of CT colography have prevented these methods from gaining widespread acceptance.

The weaknesses of colonoscopy and the current alternatives underscore the need for additional methods of CRC screening. One promising avenue is based on the field effect concept in carcinogenesis: The genetic and environmental milieu that results in a focal lesion should be detectable, at least in some form, throughout the colonic mucosa. For CRC screening, minimally invasive measurements from the rectum could be diagnostic of early carcinogenesis throughout the colon. Our laboratory has previously employed a suite of biophotonic technologies to demonstrate that changes in rectal optical signatures are indicative of premalignant lesions throughout the colon.78.9

One of the technologies we have used is known as polarization-gated spectroscopy. Interrogation of tissue structure with polarized light has found numerous applications in biomedical optics, particularly for its ability to suppress multiply scattered light.10,11 In this paper, we both demonstrate how the shape of the refractive-index correlation function can be extracted from the polarization-gated signals and present in vivo evidence from a well-controlled and established azoxymethane (AOM)-treated rat model of colon carcinogenesis that the refractive-index correlation shape can predict the risk of current and future neoplastic development throughout the colon. Changes in the shape of the refractive-index correlation function characterized by the fractal dimension have been implicated as an important parameter discriminating normal versus dysplastic tissue.1213.14.15 Furthermore, since tissue refractive index is linearly dependent on macromolecular mass density,16 the shape of the refractive-index correlation can be directly employed to assess tissue mass-density organization and structure. Consequently, there has been growing interest in developing models and techniques to measure the spatial correlations of the refractive-index properties of tissue. Our study employs a well-characterized polarization-gated spectroscopy probe17 and the Whittle-Matérn refractive-index correlation model18 to depth-selectively quantify the refractive-index correlation shape of the AOM rat model. Using this methodology, we demonstrate an association between early colon field carcinogenesis and alterations in the refractive-index correlation function shape in vivo.


Materials and Methods


Whittle-Matérn Model of Tissue Refractive-Index Distribution

The Whittle-Matérn correlation family is a three-parameter model that has been used to describe spatially continuous refractive-index fluctuations and is represented by the equation:


where the parameters lc and m denote the index correlation distance, and the shape of refractive-index correlation function, respectively, while Kv(·) represents the modified Bessel function of the second kind and N is a normalization constant. The normalization factor N ensures that Bn(0) is equal to the variance of the refractive-index fluctuations dn2 and has the following form:


where rmin is introduced to bound Bn when r approaches zero such that Bn(r<rmin)=Bn(rmin)=dn2 (Ref. 18). When m is between 1.5 and 2, the shape of the refractive-index correlation function follows a stretched exponential function. A value of m=2 results in an exponential shape while as m, the shape converges to a Gaussian distribution. It has been shown that when klc1 (where k is the wavenumber) and 1<m<2, the reduced scattering coefficient μs is proportional to λ2m4 (Ref. 18). Values of m>2, however, cannot be recovered spectrally in the limit klc1 since the wavelength dependence becomes independent of m in this region. Tissue is believed to have an m value between 1 and 2 and thus the value of m can be retrieved from the wavelength dependence of the signal intensity if the signal intensity has a known functional relationship with μs' according to the following equation:


where I is the signal intensity, f is the functional relationship between μs and I, and c is a constant from the proportionality between μs and λ2m4. We derived the functional relationship f for the polarization-gated signals using polarization-sensitive Monte Carlo (MC) that employed the Whittle-Matérn model.


Polarization-Sensitive MC Simulations Using the Whittle-Matérn Model

Publicly available polarized-light MC algorithms developed by Ramella-Roman et al.19 were used to simulate the propagation and collection of polarized light. The implementation of this algorithm to track photons as functions of polarization, radius, θc, and maximum depth travelled (z) has been described in detail by Turzhitsky et al.17 Reflection, refraction, and alteration of the Stokes vector at the sample/environment interface was computed using Snell’s Law and the Stokes formalism of the Fresnel equations.20 To employ the Whittle-Matérn model described in the previous section, we used a phase function derived from Eq. (1). First, the spectral density Φ at spatial frequency κ can be calculated from the Fourier transform of Bn(r) yielding


Then, following the derivation of Moscoso et al.,21 the phase function for incident light with initial Stokes vector [IoQoUoVo] can be written as


where θ is the angle of scattering and ϕ is the angle of rotation into the scattering plane.

We designed the MC simulations to mimic the geometry of the polarization-gated probe we have used in previous studies. This probe has overlapping illumination-collection areas of radius 400 µm and a collection angle of 14deg.17 The Stokes vectors of each photon packet were tracked and incoherently summed into corresponding bins of radius (r) giving the response to an input pencil beam. Mueller matrix multiplication was used to obtain the output intensity after a polarizer was oriented at 0 or 90 deg (I or I, respectively) with respect to the incident polarization. Convolution was employed to extend the beam from an infinitely narrow source to a circular illumination spot with finite area and radius R. The main objective of this MC study was to determine how the I or I signals behaved as a function of μs and m. Thus, m is varied between 1.1 and 2. μs is varied between 1.75 and 118cm1 for different anisotropy factors (g=0.75, 0.8, 0.85, 0.9, 0.93). By recording the signal intensity of the polarization-gated signals as a function of μs, we developed empirical mathematical relations between the signal intensity, wavelength, and m according to Eq. (2).

To confirm that we could depth-selectively recover m, we also created a two-layer MC model where the first layer had an m=1.1 and the second layer had an m=1.7. We then calculated the value m from the I and I signal intensities. As the I signal probes a deeper penetration depth than the I signal, we expected the I signal to be more sensitive to the m of the top layer and the I signal to be more sensitive to the m of the bottom layer. We also varied the thickness of the top layer to ascertain the precise thickness at which the individual m’s could be distinguished by independent analysis of the I and I signal intensities.


Determination of Length-Scale Sensitivity of the m Measurement

We followed a procedure outlined by Yi et al. to determine the range of length scales that the measurement of m is sensitive to.22 The wavelength dependence of the polarization-gated signals is primarily determined by μs(λ) which in turn can be calculated from the index correlation function. Thus to determine the length-scale sensitivity of the m measurement we perturbed the index correlation function at different length scales and calculated the corresponding change in μs(λ). To determine the inner length scale rmin, the index correlation function was truncated at rmin such that Bn(r<rmin)=Bn(rmin). This perturbed index correlation function was then used to calculate Bn. To do this, the spectral density Φ was taken as the Fourier transform of Bn and the differential scattering cross section σ was calculated as:



Next, the scattering coefficient μs was computed by integrating the differential scattering section over all angles:


The anisotropy factor g is given as the mean cosine of the scattering angle over all directions:


Finally, the reduced scattering coefficient μs can be calculated as μs(1g). These calculations were performed for m=1.4, lc=3μm, dn2=1, and repeated for different rmin and k values. This created a spectrum from which we extracted a measured m value. The error between the measured m value and the input m value was determined as a function of rmin. The inner length scale for which the measurement of m is sensitive to was defined as the minimum rmin value for which the error exceeded 5%. An analogous procedure was performed to determine the maximum length scale the m measurement is sensitive to. The perturbed index correlation function in this case was:


where rmax is the outer length scale. For r>rmax, Bn(r) decays at a power of three to reflect the lack of mass-density fluctuation. The maximum length-scale sensitivity was defined as the maximum rmax for which the error exceeded 5%.


Animal Model

This study utilized the AOM-induced rat carcinogenesis model that we have previously employed to detect morphological alterations associated with field carcinogenesis.15 Animal experimental protocols were approved by the Institutional Animal Care and Use Committee (IACUC) of NorthShore University HealthSystem and Northwestern University. Twenty-four male Fisher 344-strain rats were housed in NorthShore’s animal research facility and provided with a diet of AIN-76A (Harlan Teklad, Madison, WI) rodent chow. Animals were housed in polycarbonate, solid-bottomed cages with Alpha-dry bedding which were changed weekly. Each solid-bottom cage was supplemented with an automatic water system and ventilation ports. Each cage housed two rats in a 12-h light/dark cycle animal room at a temperature between 18 and 26ºC with minimal humidity. Animals were monitored for one week prior to injections as an acclimation period for the rodents. Rats were randomized to either a saline or AOM treatment group. At eight weeks of age, 18 rats received two weekly intraperitoneal injections of 15-mg-per-kg-of-bodyweight 13.4-moles/liter AOM (Sigma, St. Louis, MO) while the remaining six rats instead received intraperitoneal saline injections to serve as a control group.


Animal Colonoscopy and Polarization-Gated Spectroscopy Measurements

At 18 and 40 weeks post initial AOM/saline injection, rats underwent colonoscopy and polarization-gated spectroscopy measurement. Animals were fasted overnight prior to colonoscopic procedures. Rats were anaesthetized using the Veta Mac Lab Animal Anesthesia System (Veta Mac, Rossville, IN). Each animal was anaesthetized individually in an anesthesia induction chamber supplemented with (vaporized) 1–3% isoflorane (Piramal Healthcare Limited, Andhra Pradesh, India) and oxygen mixture.

After complete sedation was achieved, the rat was removed from the chamber and placed on its ventral. A nose cone was placed on the animal to continuously deliver isofluorane and maintain sedation. The colon was then lavaged to remove fecal debris by rectal injection of saline (with a ball-point needle and 50-mL syringe). Karl Storz veterinary (straight) endoscope (Storz, Goleta, CA) was inserted into the rectum and the rat colon was visualized on the video screen. The presence of polyps, their size, and their relative location were noted. Representative endoscopic images of a rat with a normal colon and a rat harboring a polyp are shown in Fig. 1(a) and 1(b), respectively. The scope was then removed and a 2.45-mm diameter fiber-optic polarization-gated spectroscopy probe was inserted into the rectum approximately 5 cm into the animal colon. Five measurements were then taken with the probe from the distal colonic mucosa. Probe measurements were taken in areas devoid of polyps as determined by the colonoscopy visualization. The probe and nose cone were then removed and the rat was monitored until effects of anesthetic were no longer noticeable.

Fig. 1

Endoscopic images of a normal saline control (a) and an AOM-treated rat with a tumor (b).



Polarization-Gated System and Data Analysis

The details of the polarization-gated probe and signal-collection system have been described previously.17 In brief, the probe consists of 3 fibers arranged in an equilateral triangle. One of these fibers delivers a linearly polarized and collimated beam of light onto the tissue surface. The two remaining fibers receive light backscattered from the tissue with one fiber collecting light polarized parallel (I) to the incident beam and the fiber collecting light polarized perpendicular (I) to the incident beam. Each collection fiber is connected to a fiber-optic spectrometer so that the signal intensity can be analyzed as a function of wavelength. The spectrum I(λ) from each tissue site and both polarization channels was analyzed according to the following relationship:


Where εOHb(λ) and εDHb(λ) are the molar extinction spectra due to the presence of oxyhemoglobin (OHb) and deoxyhemoglobin (DHb) within the tissue, α is the product of OHb/DHb concentration and pathlength, and Cpacking(λ) is a correction factor for the inhomogeneous distribution (packing) of hemoglobin into blood vessels:2324.25


where R is the average blood-vessel radius and μa,bl is the absorption coefficient of whole blood. Methods to separate the absorbance spectrum from the scattering spectrum f(m,cλ2m4) have been further described in detail in previous publications.7,17,26 The precise functional form of f(m,cλ2m4) was determined using the MC methods described in Secs. 2.1 and 2.2 and is described in the results section. Once the absorbance spectrum is removed from the signal, the value of m can be fit to the resultant scattering spectrum using the least-squares minimization technique.




MC-Derived Relationship between Polarization-Gated Signal Intensities and μs

In Fig. 2(a), we plot the I signal intensity versus μs for several different anisotropy factors and m=1.5. For g>0.85, the individual curves for different g are within 5%, indicating that the I signal can be written as a function of μs rather than a combination of g and μs. In Fig. 2(b), we plot the same I signal intensity versus μs curves but for different m and a single anisotropy factor (g=0.9). These curves clearly show that the value of m influences the I-signal-intensity-versus-μs-relationship. However, we found that each curve could be fit with a power law such that I(μs)C(m) where the exponent C(m) is a function of m. The power law fit each of the curves well with an R2>0.99 in all cases. The exponent C(m) was found to be a linear function of m with an R2>0.99 as shown in Fig. 2(c). Using the proportionality between μs and λ2m4 derived from the Whittle-Matern model and the linear fit to C(m), we can express the wavelength dependence of the copolarized signal I(λ) as:


We repeated this procedure for the crosspolarized signal I(λ) in Fig. 2(d)2(f). The intensity-versus- μs curves overlapped for all g’s as shown in Fig. 2(d) with percentage deviation due to g of <1%. In Fig. 2(e), we demonstrate that different m values give rise to different crosspolarized-intensity-versus- curves but that each curve can be fit with a power law (R2>0.99). Finally, in Fig. 2(f), we show that the exponent of the power law follows a linear trend with m in the same manner that it did for the copolarized signal. Following this analysis, the wavelength dependence of the crosspolarized signal I(λ) can be written as:


By fitting Eqs. (12) and (13) to the copolarized and crosspolarized signals, the value of m can be extracted from polarization-gated measurements.

Fig. 2

Extraction of the index correlation-function shape (m) from the copolarized intensity (top row) and crosspolarized intensity (bottom row). (a) For g>0.85, the intensity is a function of and independent of g. A power law can be fit to this dependence for different m as shown in (b). The exponent of the power law has a linear relationship with m as shown in (c). (d–f) are analogous to (a–c) but for the crosspolarized intensity rather than copolarized intensity.



Measurement of m from a Dual-Layer MC Model

The analysis so far assumes a homogeneous medium defined by a single m value. In reality, biological tissue is a layered medium most likely defined by more than a single m. Since the copolarized and crosspolarized signals sample different penetration depths (350 µm and 575 µm on average, respectively), it is likely that the m values measured by these individual signals will be different, with the copolarized signal measuring the m value of superficial layers of and the crosspolarized signal measuring the m value of deeper layers. To test this hypothesis, we created a two-layer MC model in which the top layer had an m of 1.1 and the bottom layer an m of 1.7. The optical thickness τ of the first layer, defined as the product of the physical thickness and the scattering coefficient μs, was varied between 0 and 100 while the overall τ of the medium was fixed at 1000. Figure 3 shows the m values measured from the co- and crosspolarized intensities as a function of the optical thickness of the top layer. When τ of the top layer is less than 5, both co- and cross-signals predominantly interrogate the bottom layer (m=1.7). For τ greater than 15, both signals probe the top layer (m=1.1). For 5<τ<15, the copolarized signal probes mainly the top layer and thus m is close to 1.1. In contrast, the crosspolarized signal measures the bottom layer and has an m value closer to 1.7.

Fig. 3

Depth-selective extraction of index correlation shape (m) from a dual-layer Monte Carlo model. The top layer has m=1.1 and bottom layer an m of 1.7. The m values from the copolarized and crosspolarized signals are recorded as the optical thickness (τ) of the top layer is varied. For 5<τ<15, the copolarized m parameter is sensitive to the top layer and the crosspolarized m parameter is sensitive to the bottom layer.



Length-Scale Sensitivity of m Measurements

We determined the minimum and maximum length scale over which our spectral determination of m is sensitive to by perturbing the index correlation function with different rmin and rmax values as outlined in Sec. 2.3. Figure 4(a) and 4(b) display the percent error between the input m value and the measured m value from the perturbed index correlation function for different rmin and rmax, respectively. We can define a percent error threshold whereby if the percent error is greater than this threshold for a particular rmin and rmax, then the length scales to which the measurement is sensitive are between rmin and rmax. Defining a threshold of 5% error suggests that the measurement of the m value is sensitive to length scales between 45 nm and 1.5 µm according to Fig. 4.

Fig. 4

Length-scale sensitivity of index correlation function shape (m) measurement. (a) Percent between measured and input m values after correlation function was truncated at a particular rmin value. (b) Percent error between measured and input m values after index correlation function was perturbed at various rmax values according to Eq. (9).



Tumor Development in Saline/AOM-treated Rats

Rats were assessed at 18 and 40 weeks for tumor development by colonoscopy. The number of rats who had developed tumors by each time point is summarized in Table 1. As expected, none of the saline-treated rats developed tumors over the study time period. 77% of the AOM-treated rats possessed tumors by the end of 40 weeks, with more than half developing tumors by the 18-week time point. This is consistent with earlier studies on tumor development in the AOM-treated rat.27 These results confirm that the saline-treated rats could be reliably regarded as controls and that the AOM treatment was a potent initiator of carcinogenesis.

Table 1

Number of rats that developed tumors at the 18- and 40- week timepoints.

18 weeks40 weeksDid not develop18 weeks40 weeksDid not develop


Measurement of the Refractive-Index Correlation-Function Shape from Salineand AOM Treated Rats

After colonoscopy at the 18-week time point, five polarization-gated probe measurements were taken throughout the visually normal distal mucosa of each saline and AOM-treated rat. For each tissue site, an m value was extracted from the copolarized and crosspolarized signals according to Eqs. (7) and (8). The median m value from these five sites was taken and used to characterize the mucosa of each individual rat. A Welch’s t-test was used to determine statistical significance between groups. The first comparison we made was between saline-treated rats and all rats treated with AOM. Figure 5(a) shows that in rats treated with AOM, there is a statistically significant increase in the m parameter from the visually normal mucosa for both the co- and crosspolarized signals (P value=0.007 and 0.003, respectively). The average±standard deviation m values for saline controls were 1.77±0.04 and 1.64±0.06 from the co- and crosspolarized signals, respectively, while for AOM they were 1.9±0.17 and 1.79±0.16 from the co- and crosspolarized signals, respectively. Figure 5(b) illustrated what the shape of the refractive-index correlation would look like for saline and AOM-treated rats given the m values from the crosspolarized signal. Both curves can be represented by stretched exponentials with the AOM curve decaying more slowly, suggesting that the AOM-rat tissue structure has correlations at longer length scales than the saline-rat tissue structure does. We next investigated whether the m parameter could be a predictor of both concurrent and future risk of neoplasia. In Fig. 5(c), we compare the m parameter values measured at 18 weeks among saline controls, AOM rats that developed tumors by 18 weeks (tumor presence), and AOM rats that had not developed tumors by 18 weeks (no tumor presence). This figure demonstrates that the m parameter mirrors neoplasia risk. For both the copolarized and crosspolarized m parameter, there is an increase in m for rats that had not developed tumors by the 18-week time point but were still at risk of developing tumors because of the AOM treatment (P=0.1 and 0.007, respectively). There is a further increase in m in rats that harbored neoplasia at the time of measurement as shown in Fig. 5(b). These results demonstrate that the m parameter is most sensitive to the presence of neoplasia, but that it is also predictive of the risk of future neoplastic development.

Fig. 5

The shape of the index correlation function is sensitive to concurrent and future risk of neoplasia in the AOM-treated rat. (a) Comparison of the m values from the copolarized and crosspolarized signals for saline versus all AOM-treated rats shows a statistical significant increase in m (P value=0.007 and 0.003). (b). Comparison of the shape of index correlation function shape between saline and AOM-treated rats. (c). Comparison of m values between saline controls and AOM-treated rats who had and had not developed tumors at the 18-week measurement point. The m value is highest for rats harboring a tumor but is also significant for AOM rats that had not yet developed a tumor. Solid circle points are group means and error bars are 95% CI.




The results in Sec. 3.4 demonstrate that the m parameter measured from the wavelength dependence of the backscattered tissue intensity can demarcate tumor risk in the AOM-rat model. Previous studies have employed light reflectance spectroscopy to document alterations in the spectal slope I/λ from tumors compared to normal tissue.2829.30 In addition, prior studies from our laboratory on ex vivo samples used elastic light scattering to demonstrate that the spectral slope was a marker of tumor development in the histologically normal mucosa of the AOM rat and human patients.31 It was assumed that the spectral slope was related to the size distribution of scattering structures. In this paper, we have used the Whittle-Matérn correlation function family and Monte-Carlo simulations to formalize the relationship between the wavelength dependence of the signal intensity and the shape of the index correlation function. This allows one to make quantitative comparisons that directly relate to a physical parameter (mass-density distribution) of the tissue structure. It should be noted that our previous ex vivo light-scattering studies showing a decline in spectral slope with precancerous alterations are consistent with our in vivo evidence of an increase in the m parameter according to Eqs. (1)–(3).

In Sec. 3.2 we demonstrated that the co- and crosspolarized signals could depth-selectively assess the m parameter in a two-layer model given that the optical thickness of the top layer was between approximately 5 and 15τ. The thickness of the top layer in the human colon, the mucosa, is on the order of 0.05 cm while the scattering coefficient is about 200cm1 giving an optical thickness of 10.32,33 In rats, the thickness of the mucosa is 0.02cm,34 giving an assumed mucosal optical thickness of 5. Thus, it is possible to depth-selectively quantify the m parameter for both animal and human colon tissue under the above assumptions. The importance of depth selectivity for accurate assessment of early cancer risk has been demonstrated by several techniques. Using polarization-gated spectroscopy, we have previously shown that early increases in bloody supply associated with carcinogenesis are more prominent at superficial tissue layers.7,26 Low coherence interferometry has been used to demonstrate that the most significant changes in nuclear diameter in the AOM-treated rat occur at a depth of 35 µm.35 In the present study, both the copolarization and crosspolarization m parameter were statistically significant indicators of concurrent neoplasia, suggesting diffuse alterations throughout mucosa and submucosa in the shape of the index correlation function.

The values of m identified in this study ranged from 1.6 to 2, with AOM-treated rats having higher m values than saline controls. These values of m correspond to the index correlation function having the shape of a stretched exponential. Stretched exponential shapes have been observed in alterations associated with ovarian cancer and were linked to scattering objects having similar size.36 The control samples from that study had m<1.5 placing them in the mass fractal regime. Previous studies on rat esophagi also measured control samples to be in the mass fractal regime.13 The precise reason for why our control measurements did not fall into the mass fractal regime is not known, but may have to do with the type of tissue investigated, the depth of the tissue sampled, and the fact that the previous studies mentioned utilized ex vivo biopsies instead of in vivo measurements.

While the results from this study are very promising, there are certain limitations that need to be addressed. The derivation of Eqs. (12) and (13) depends on certain assumptions, the most important being that actual biological tissue conforms to the Whittle-Matérn correlation family. While this correlation family is likely an improvement over functions that treat tissue as discrete scattering particles such as Mie theory, there is a lack of experimental evidence to confirm this. In addition, both the single and dual-layer Monte-Carlo models are simplified models of tissue. In reality, tissue possesses more than two layers and the distribution of optical properties within each layer is not homogeneous. We intend to explore in future studies how modeling a more realistic tissue structure impacts our measurement and interpretation of the index correlation function shape.



In this paper we have presented a method for measuring the shape of the tissue refractive-index correlation function from polarization-gated spectroscopy measurements. We have employed this method to measure the function shape in the AOM-treated rat model of colon carcinogenesis and we have found that changes in the shape were correlated with both present and future risk of neoplasia development. Overall, our study provides in vivo evidence from a well-controlled animal model that depth-selective measurements of the index correlation shape can serve as a biomarker for CRC risk assessment, as well as evidence for the field-effect concept in carcinogenesis.



R. Siegelet al., “Cancer statistics, 2011: the impact of eliminating socioeconomic and racial disparities on premature cancer deaths,” CA Cancer J. Clin.. 61(4), 212–236 (2011).CAMCAM0007-9235http://dx.doi.org/10.3322/caac.v61:4Google Scholar


“Screening for colorectal cancer: U.S. Preventive Services Task Force recommendation statement,” Ann. Intern. Med. 149(9), 627–637 (2008).AIMEAS0003-4819Google Scholar


J. C. van Rijnet al., “Polyp miss rate determined by tandem colonoscopy: a systematic review,” Am. J. Gastroenterol. 101(2), 343–350 (2006).AJGAAR0002-9270http://dx.doi.org/10.1111/ajg.2006.101.issue-2Google Scholar


H. K. RoyV. BackmanM. J. Goldberg, “Colon cancer screening: the good, the bad, and the ugly,” Arch. Intern. Med. 166(20), 2177–2179 (2006).AIMDAP0003-9926http://dx.doi.org/10.1001/archinte.166.20.2177Google Scholar


T. F. Imperialeet al., “Fecal DNA versus fecal occult blood for colorectal-cancer screening in an average-risk population,” N. Engl. J. Med. 351(26), 2704–2714 (2004).NEJMAG0028-4793http://dx.doi.org/10.1056/NEJMoa033403Google Scholar


J. D. Lewiset al., “Detection of proximal adenomatous polyps with screening sigmoidoscopy: a systematic review and meta-analysis of screening colonoscopy,” Arch. Intern. Med. 163(4), 413–420 (2003).AIMDAP0003-9926http://dx.doi.org/10.1001/archinte.163.4.413Google Scholar


A. J. Gomeset al., “Rectal mucosal microvascular blood supply increase is associated with colonic neoplasia,” Clin. Cancer Res. 15(9), 3110–3117 (2009).CCREF41078-0432http://dx.doi.org/10.1158/1078-0432.CCR-08-2880Google Scholar


H. K. Royet al., “Association between rectal optical signatures and colonic neoplasia: potential applications for screening,” Cancer Res. 69(10), 4476–4483 (2009).CNREA80008-5472http://dx.doi.org/10.1158/0008-5472.CAN-08-4780Google Scholar


H. Subramanianet al., “Nanoscale cellular changes in field carcinogenesis detected by partial wave spectroscopy,” Cancer Res. 69(13), 5357–5363 (2009).CNREA80008-5472http://dx.doi.org/10.1158/0008-5472.CAN-08-3895Google Scholar


N. GhoshI. A. Vitkin, “Tissue polarimetry: concepts, challenges, applications, and outlook,” J. Biomed. Opt. 16(11), 110801 (2011).JBOPFO1083-3668http://dx.doi.org/10.1117/1.3652896Google Scholar


V. V. TuchinL. V. WangD. A. Zimnyakov, Optical Polarization in Biomedical Applications, Springer, Berlin (2006).Google Scholar


J. W. BaishR. K. Jain, “Fractals and cancer,” Cancer Res. 60(14), 3683–3688 (2000).CNREA80008-5472Google Scholar


M. Hunteret al., “Tissue self-affinity and polarized light scattering in the born approximation: a new model for precancer detection,” Phys. Rev. Lett. 97(13), 138102 (2006).PRLTAO0031-9007http://dx.doi.org/10.1103/PhysRevLett.97.138102Google Scholar


H. K. Royet al., “Down-regulation of SNAIL suppresses MIN mouse tumorigenesis: modulation of apoptosis, proliferation, and fractal dimension,” Mol. Cancer Ther. 3(9), 1159–1165 (2004).MCTOCF1535-7163Google Scholar


H. K. Royet al., “Four-dimensional elastic light-scattering fingerprints as preneoplastic markers in the rat model of colon carcinogenesis,” Gastroenterology 126(4), 1071–1081 and 1948 (2004).GASTAB0016-5085http://dx.doi.org/10.1053/j.gastro.2004.01.009Google Scholar


R. BarerS. Tkaczyk, “Refractive index of concentrated protein solutions,” Nature 173(4409), 821–822 (1954).NATUAS0028-0836http://dx.doi.org/10.1038/173821b0Google Scholar


V. M. Turzhitskyet al., “Measuring mucosal blood supply in vivo with a polarization-gating probe,” Appl. Opt. 47(32), 6046–6057 (2008).APOPAI0003-6935http://dx.doi.org/10.1364/AO.47.006046Google Scholar


J. D. RogersI. R. CapogluV. Backman, “Nonscalar elastic light scattering from continuous random media in the Born approximation,” Opt. Lett. 34(12), 1891–1893 (2009).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.34.001891Google Scholar


J. C. Ramella-RomanS. A. PrahlS. L. Jacques, “Three Monte Carlo programs of polarized light transport into scattering media: part II,” Opt. Express 13(25), 10392–10405 (2005).OPEXFF1094-4087http://dx.doi.org/10.1364/OPEX.13.010392Google Scholar


F. JaillonH. Saint-Jalmes, “Description and time reduction of a Monte Carlo code to simulate propagation of polarized light through scattering media,” Appl. Opt. 42(16), 3290–3296 (2003).APOPAI0003-6935http://dx.doi.org/10.1364/AO.42.003290Google Scholar


M. MoscosoJ. B. KellerG. Papanicolaou, “Depolarization and blurring of optical images by biological tissue,” J. Opt. Soc. Am. A. Opt. Image Sci. Vis. 18(4), 948–960 (2001).Google Scholar


J. Yiet al., Inverse Scattering Optical Coherence Tomography (ISOCT): Quantifying Sub-Diffractional Tissue Mass Density Correlation Function, PNAS, Northwestern Unviersity, Evanston (2011).Google Scholar


J. C. FinlayT. H. Foster, “Effect of pigment packaging on diffuse reflectance spectroscopy of samples containing red blood cells,” Opt. Lett. 29(9), 965–967 (2004).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.29.000965Google Scholar


L. Svaasandet al., “Therapeutic response during pulsed laser treatment of port-wine stains: dependence on vessel diameter and depth in dermis,” Las. Med. Sci. 10(4), 235–243 (1995).LMSCEZ1435-604Xhttp://dx.doi.org/10.1007/BF02133615Google Scholar


R. L. van VeenW. VerkruysseH. J. Sterenborg, “Diffuse-reflectance spectroscopy from 500 to 1060 nm by correction for inhomogeneously distributed absorbers,” Opt. Lett. 27(4), 246–248 (2002).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.27.000246Google Scholar


H. K. Royet al., “Spectroscopic microvascular blood detection from the endoscopically normal colonic mucosa: biomarker for neoplasia risk,” Gastroenterology 135(4), 1069–1078 (2008).GASTAB0016-5085http://dx.doi.org/10.1053/j.gastro.2008.06.046Google Scholar


M. Kobaek-Larsenet al., “Review of colorectal cancer and its metastases in rodent models: comparative aspects with those in humans,” Comp. Med. 50(1), 16–26 (2000).CMPMDIGoogle Scholar


J. R. Mourantet al., “Spectroscopic diagnosis of bladder cancer with elastic light scattering,” Laser Surg. Med. 17(4), 350–357 (1995).LSMEDI0196-8092http://dx.doi.org/10.1002/(ISSN)1096-9101Google Scholar


U. Utzingeret al., “Reflectance spectroscopy for in vivo characterization of ovarian tissue,” Laser Surg. Med. 28(1), 56–66 (2001).LSMEDI0196-8092http://dx.doi.org/10.1002/(ISSN)1096-9101Google Scholar


I. J. Bigioet al., “Diagnosis of breast cancer using elastic-scattering spectroscopy: preliminary clinical results,” J. Biomed. Opt. 5(2), 221–228 (2000).JBOPFO1083-3668http://dx.doi.org/10.1117/1.429990Google Scholar


H. K. Royet al., “Spectral slope from the endoscopically-normal mucosa predicts concurrent colonic neoplasia: a pilot ex vivo clinical study,” Dis. Colon Rectum. 51(9), 1381–1386 (2008).DICRAG0012-3706http://dx.doi.org/10.1007/s10350-008-9384-3Google Scholar


Z. Huanget al., “Laser-induced autofluorescence microscopy of normal and tumor human colonic tissue,” Int. J. Oncol. 24(1), 59–63 (2004).IJONES1019-6439Google Scholar


H. J. Weiet al., “Differences in optical properties between healthy and pathological human colon tissues using a Ti:sapphire laser: an in vitro study using the Monte Carlo inversion technique,” J. Biomed. Opt. 10(4), 44022 (2005).JBOPFO1083-3668http://dx.doi.org/10.1117/1.1990125Google Scholar


R. Sharmaet al., “Rat intestinal mucosal responses to a microbial flora and different diets,” Gut 36(2), 209–214 (1995).GUTTAK0017-5749http://dx.doi.org/10.1136/gut.36.2.209Google Scholar


F. E. Robleset al., “Detection of early colorectal cancer development in the azoxymethane rat carcinogenesis model with Fourier domain low coherence interferometry,” Biomed. Opt. Express 1(2), 736–745 (2010).BOEICL2156-7085http://dx.doi.org/10.1364/BOE.1.000736Google Scholar


O. Nadiarnykh et al., “Alterations of the extracellular matrix in ovarian cancer studied by Second Harmonic Generation imaging microscopy,” BMC Cancer 10(94) (2010)http://dx.doi.org/10.1186/1471-2407-10-94.BCMACL1471-2407Google Scholar

Andrew J. Gomes, Sarah K. Ruderman, Vadim Backman, Mart Dela Cruz, Ramesh K. Wali, Hemant K. Roy, "In vivo measurement of the shape of the tissue-refractive-index correlation function and its applicationto detection of colorectal field carcinogenesis," Journal of Biomedical Optics 17(4), 047005 (18 April 2012). http://dx.doi.org/10.1117/1.JBO.17.4.047005

Correlation function



Tissue optics

In vivo imaging

Bragg cells


Back to Top