OSA's Digital Library

Virtual Journal for Biomedical Optics

Virtual Journal for Biomedical Optics

| EXPLORING THE INTERFACE OF LIGHT AND BIOMEDICINE

  • Editor: Gregory W. Faris
  • Vol. 5, Iss. 9 — Jul. 6, 2010

Optics InfoBase > Virtual Journal for Biomedical Optics > Volume 5 > Issue 9 > Page 10538

« Show journal navigation

Measurement of the ratio of forward-propagating to back-propagating second harmonic signal using a single objective

Xiaoxing Han and Edward Brown  »View Author Affiliations


Optics Express, Vol. 18, Issue 10, pp. 10538-10550 (2010)
http://dx.doi.org/10.1364/OE.18.010538


View Full Text Article

Acrobat PDF (1152 KB)


Advanced Search

Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

In this paper, we present a method to determine, for the first time, the ratio of forward-propagating second harmonic (SHG) signal to back-propagating SHG signal (F/B) in vivo on the surface of intact tissue samples without any biopsy or tissue sectioning, using only epidetection (i.e., via a single objective lens). This method has the additional benefit of using the confocal detection apparatus already contained within common commercially available two-photon laser-scanning microscopes, and hence can allow the measurement of the SHG F/B ratio in vivo with minimal purchase of new equipment.

© 2010 OSA

1. Introduction

Recently, second harmonic generation (SHG) has proven to be a useful window into the amount and organization of fibrillar collagen in biological tissues due to its relative specificity and the fact that it is an intrinsic signal [1

R. M. Williams, W. R. Zipfel, and W. W. Webb, “Interpreting second-harmonic generation images of collagen I fibrils,” Biophys. J. 88(2), 1377–1386 (2005). [CrossRef]

–6

S. V. Plotnikov, A. C. Millard, P. J. Campagnola, and W. A. Mohler, “Characterization of the myosin-based source for second-harmonic generation from muscle sarcomeres,” Biophys. J. 90(2), 693–703 (2006). [CrossRef]

]. SHG is a coherent phenomenon, which implies that SHG is sensitive not only to the amplitude of the illumination field but also to its phase. In addition to a spatial resolution that is equal to other imaging techniques (such as two photon excited fluorescence), SHG microscopy can provide information about the sample’s molecular structure. For example, the ratio of the forward-propagating to backward propagating SHG signal (the “F/B ratio”) can help us to understand the axial extent of ordering in collagen fibers [1

R. M. Williams, W. R. Zipfel, and W. W. Webb, “Interpreting second-harmonic generation images of collagen I fibrils,” Biophys. J. 88(2), 1377–1386 (2005). [CrossRef]

–5

A. C. Kwan, D. A. Dombeck, and W. W. Webb, “Polarized microtubule arrays in apical dendrites and axons,” Proc. Natl. Acad. Sci. U.S.A. 105(32), 11370–11375 (2008). [CrossRef] [PubMed]

].

Previously, in vitro measurements of SHG F/B ratios have been used to study collagen fiber ordering in various tissue samples such as rat tail, ovarian cancer biopsies, mouse models of breast cancer, and dermis from mouse models of Osteogenesis Imperfecta (OIM) [1

R. M. Williams, W. R. Zipfel, and W. W. Webb, “Interpreting second-harmonic generation images of collagen I fibrils,” Biophys. J. 88(2), 1377–1386 (2005). [CrossRef]

–5

A. C. Kwan, D. A. Dombeck, and W. W. Webb, “Polarized microtubule arrays in apical dendrites and axons,” Proc. Natl. Acad. Sci. U.S.A. 105(32), 11370–11375 (2008). [CrossRef] [PubMed]

,7

O. Nadiarnykh, R. B. LaComb, M. A. Brewer, and P. J. Campagnola, “Alterations of the extracellular matrix in ovarian cancer studied by second harmonic generation imaging microscopy,” BMC Cancer 10(1), 94 (2010). [CrossRef] [PubMed]

]. The F/B ratio revealed the length scale of ordering in the fibers and in the case of OIM and ovarian cancer, was able to discriminate pathological tissue from healthy tissue [2

R. Lacomb, O. Nadiarnykh, and P. J. Campagnola, “Quantitative second harmonic generation imaging of the diseased state osteogenesis imperfecta: experiment and simulation,” Biophys. J. 94(11), 4504–4514 (2008). [CrossRef] [PubMed]

,7

O. Nadiarnykh, R. B. LaComb, M. A. Brewer, and P. J. Campagnola, “Alterations of the extracellular matrix in ovarian cancer studied by second harmonic generation imaging microscopy,” BMC Cancer 10(1), 94 (2010). [CrossRef] [PubMed]

]. In these measurements a second objective lens was needed to collect forward propagating SHG signal. Hence, the tissue sample had to be dissected from the animal and sectioned to 100um slices to allow signal to reach the second detection lens. For clinical application, such as in endoscopy, it is impossible to put an objective and a PMT detector underneath the tissue sample to collect the forward propagating SHG. The excision and sectioning required to use a second detector for forward propagating SHG also prevents dynamic measurements of collagen ordering over time. Therefore, it becomes highly desirable to develop a new optical system to measure collagen SHG F/B ratio in vivo, using only epi-detection and on intact thick tissue samples.

2. Experimental methods

2.1 Optical setup

The optical setup we are using to determine the SHG F/B ratio in vivo using only epidetection is shown in Fig. 1 below. SHG signal was generated by a Spectra Physics MaiTai Ti:Sapphire laser providing 100fs pulses at 80 MHz and 810 nm. The excitation beam was directed into an Olympus Fluoview F300 scanhead connected to an Olympus BX61WI upright microscope. The focusing objective is an Olympus UMPLFL20XW water immersion lenses (20 × , 0.5 N.A.), with supporting measurements performed using an Olympus LUMPlanFL/IR (40x, 0.8 NA) water immersion lens. This objective was used to focus the excitation beam on the sample and at the same time collect the direct backward-propagating SHG as well as the forward-propagating SHG signal that was subsequently backscattered within the tissue. The SHG signal was collected by the objective, converged by the tube lens, and collimated by the pupil lens. The collimated SHG beam was then de-scanned and focused again on the pinhole plane by the collector lens. All of these lenses are intrinsic to the BX61WI microscope or the Fluoview F300 scanhead. The focal lengths of the pupil lens and the collector lens are 54mm and 185mm respectively, and size of the pinholes on the F300 variable pinhole turret are 60um, 100um, 150um, 200um and 7000um. The SHG beam was separated from the excitation beam by a dichroic mirror (Chroma 670 DCSX) inserted into the F300 scanhead, as well as a band pass filter centered at 405 nm (Chroma HQ405/30m-2P) placed after the pinhole, and detected by the FV300’s intrinsic photomultiplier tubes. The only significant equipment modifications are insertion of the aforementioned dichroic mirror and filter in the appropriate location in the scanhead, and the punching out of one of the pinholes in the pinhole turret to produce one extremely large pinhole setting (see below).

Fig. 1 Experimental setup for in vivo measurement of tumor collagen SHG F/B ratio.

2.2 SHG Back Scattering and Monte Carlo Simulation

Using the apparatus described above, we generate an SHG image using confocal detection of the resultant SHG. While much of the forward-propagating SHG signal will travel into the tissue sample and be lost to our objective lens, a fraction of this signal undergoes multiple scattering events and passes back through the object plane, traveling towards the objective lens. As shown in Fig. 2 , this signal will pass through the object plane at multiple locations. Conversely, at shallow imaging depths the backwards-propagating SHG signal will emanate from the image plane only from the two-photon focal volume with minimal subsequent scatter. When the object plane is imaged onto the confocal pinhole, the spatial distribution of SHG signal on the confocal plane will consist of a sharp central peak due to the backward propagating SHG plus a diffuse signal due to the forwards propagating and subsequently backscattered SHG. As discussed below, our new method consists of repeatedly imaging the sample through a series of different sized pinholes, whereby the shape of this total SHG distribution can be measured and, with suitable calibration, the underlying F/B ratio can be determined.

Fig. 2 Back scattering of forward-propagating SHG and propagation of direct backward propagating SHG

In order to implement this method we must first understand the spatial distribution, in the object plane, of forward propagated light that subsequently backscatters and reaches that plane. 405 nm SHG signal is easily scattered while propagating in the tissue sample. In Legare et al 2007 [8

F. Légaré, C. Pfeffer, and B. R. Olsen, “The role of backscattering in SHG tissue imaging,” Biophys. J. 93(4), 1312–1320 (2007). [CrossRef] [PubMed]

], it was determined that ~21% of the forward-propagating SHG signal was subsequently backscattered and reached the tissue surface, and that the epi-detected image from a 5mm thick Achilles tendon tissue block was actually a combination of the direct backward propagating SHG signal and the forward propagating SHG signal that was subsequently backscattered. In fact for collagen fibers in some tissue samples, such as in mouse breast tumor models, the SHG F/B ratio is very high (~30) [4

X. Han, R. M. Burke, M. L. Zettel, P. Tang, and E. B. Brown, “Second harmonic properties of tumor collagen: determining the structural relationship between reactive stroma and healthy stroma,” Opt. Express 16(3), 1846–1859 (2008). [CrossRef] [PubMed]

], and most SHG signal propagates in the forward direction. In these cases epi-detected SHG microscopy is greatly enhanced due to the subsequent backscattering of the forward propagating SHG.

The propagation of light in turbid media can be well modeled by Monte Carlo simulation of the paths that photons make as they travel through tissue, which are chosen by statistically sampling the probability distributions for step size and angular deflection per scattering event. Using a Monte Carlo code based closely upon that of Wang et al [9

L. Wang and S. L. Jacques, “Hybrid model of Monte Carlo simulation and diffusion theory for light reflectance by turbid media,” J. Opt. Soc. Am. A 10(8), 1746–1752 (1993). [CrossRef]

] (see Appendix), we sent 405nm SHG photons into an infinitely deep scattering tissue in the forward direction. The SHG photons were emitted from a point source and we ignored the direct backward propagating SHG photons. We experimentally measured the angular distribution of the forward propagating SHG emission by imaging the back focal plane of an objective lens which was collecting the forward propagating SHG from a 100 um section of tumor tissue, and found that the forward SHG from collagen fibers is emitted in a rather tight beam confined to +/−15 degrees around the laser axis (data not shown). Hence we used the same initial angular distribution (+/−15 degrees around the axis normal to the tissue-air surface) in the Monte Carlo simulation. The scattering particles are modeled as 10um diameter cells, the refractive indexes inside and outside of the cells are 1.38 and 1.42 respectively, and the scattering and absorption coefficients we used in the simulation are µa = 0.7 cm−1 and µs = 150 cm−1 [10

W. F. Cheong, S. A. Prahl, and A. J. Welch, “A review of the optical properties of biological tissues,” IEEE J. Quantum Electron. 26(12), 2166–2185 (1990). [CrossRef]

]. The probability distribution of scattering angles after every scattering event were calculated according to Mie theory with the cell size and the refractive index inside and outside of the cells, with a 99% chance the scattering angle will fall between +/−6 degrees. We then counted the number of photons which escape from the air-tissue surface in the backward direction and plotted out the steady state distribution of the backscattered photons over the radial distance from the emission point.

The simulation results are shown in Fig. 3 . Figure 3(a) shows the steady state distribution of the forward propagating and subsequently backscattered SHG photons which reach the object plane, over a large range of radial positions i.e. from the center to 3mm away from the emission point. Figure 3(b) shows the same distribution over a short range i.e. from the center to only 50um away from the emission point, which is the length scale of our largest pinhole setting (where the pinhole-bearing foil has been punched out). We can see from Fig. 3(a) that over a large length scale the forward propagating and subsequently back-scattered SHG photon intensity decays exponentially with the distance from emission point, which agrees with previous literature [11

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

]. But from Fig. 3(b) we can see over a very short range of 50 um, the range covered by our pinhole settings and where the diffusion approximation does not hold [11

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

], the intensity of forward propagating and subsequently back scattered SHG photons is constant and does not depend on distance from the focal volume.

Fig. 3 Monte Carlo simulation of forward propagating and subsequently backscattered SHG photons which reach the object plane.(a) Steady state radial distribution of the backscattered SHG photons over a large range of radial position (3mm) (b) The same distribution over a short radial range (50um)

2.3 Fitting model

To measure collagen fiber SHG F/B ratio in vivo with only the epi-detection objective, we will generate a series of collagen fiber SHG images, in the backward channel, through a series of confocal pinholes of different sizes which correspond to a diameter in the object plane ranging from 0.874 um to 102 um (see 3.1 below). As discussed above, the SHG signal we collect after the pinhole is a mixture of both the direct backward propagating SHG and the forward propagating and subsequently backscattered SHG. To analyze our experimental data we model the image in the pinhole plane of the direct backward propagating SHG signal plane as a Gaussian spot [12

W. R. Zipfel, R. M. Williams, and W. W. Webb, “Nonlinear magic: multiphoton microscopy in the biosciences,” Nat. Biotechnol. 21(11), 1369–1377 (2003). [CrossRef] [PubMed]

], while the forward propagating SHG that subsequently backscatters is modeled as a uniform distribution over these length scales (based upon the results of our Monte Carlo simulation, as described above). So the total SHG signal intensity distribution on the object plane can be expressed as:
I SHG(r)=Bexp[−2 ( rω)2]+FC,
(1)
where ω is the e−2 Gaussian spot size of the direct backward propagating SHG, F and B are absolute intensities of forward and backward propagating SHG signals, the parameter C relates the initial forward propagating signal intensity to the average intensity of the uniform distribution of SHG light that reaches the object plane, and is a function of scattering and absorption properties of the underlying tissue. Alternatively, this expression can be written in another way,
I SHG(r)=B[exp[−2 ( rω)2]+ FBC],
(2)
where F/B represents the collagen fiber SHG F/B ratio.

When we generate collagen fiber SHG images through a series of confocal pinholes of different sizes, each pixel on that image represents an integration of the total SHG signal over the pinhole area
I pixel∝ ∫0 2π dθ ∫0R { exp[−2 ( rω)2]+ ( FB)C}rdr,
(3)
where R is the size of the pinhole with respect to Gaussian spot size of the direct backward propagating SHG, i.e. R = rpinhole/ω. If the pixel intensities at various pinhole sizes are normalized to the maximum pixel intensity at the largest pinhole size, the relative pixel intensity is a function of relative pinhole size R:
I rel(R)= ∫0 2π dθ ∫0R { exp[−2 ( rω)2]+ ( FB)C}rdr ∫0 2π dθ ∫0 R max { exp[−2 ( rω)2]+ ( FB)C}rdr,
(4)
where Rmax is the largest pinhole size in our system. With this expression we can plot the relative pixel intensity versus pinhole size and fit the data to produce (F/B)C. In order to determine F/B we must then eliminate C, the fraction of signal which originally propagates in the forward direction but is eventually backscattered by the tissue and reaches the pinhole plane. This can be done by introducing into the sample a reference of known F/B ratio, whose C value is the same. In these experiments we used blue fluorescent polystyrene beads (peak λem = 424nm) sprinkled onto the sample surface. The quantity (F/B)C is then determined for both the calibration beads in the image, as well as the fibers of interest. Next, the correction factor that is required to convert the measured values of bead (F/B)C to the known value of F/B is determined, and the previously unknown value of fiber F/B is determined from the measured (F/B)C using the same correction factor:

measure d collage n SH G F/ B Ratio measure d bead s TPE F F/ B Ratio= rea l collage n SH G F/ B Ratio rea l bead s TPE F F/ B Ratio.
(5)

To illustrate the fitting function given by Eq. (4) in the case of the typical backscattering fraction C = 0.001 and Rmax = 28 we can plot out the relative pixel intensity v.s. pinhole size curves at various collagen fiber SHG F/B ratios in Fig. 4 :

Fig. 4 Relative SHG intensity vs. relative pinhole size curves at different collagen SHG F/B ratios, based upon Eq. (4). The relative SHG intensities equal 1 at a pinhole size of r = Rmax = 28.

From this plot one can see the sharp curvature due to the Gaussian distribution of direct back-propagating SHG (most evident in the steep early rise of the curves with small values of F/B) and the slow and steady rise due to forward-propagating light that is subsequently backscattered and produces a diffuse signal which does not vary with r (most evident in the slow rise of the curves with large values of F/B). This reveals the utility in punching out one of the metal foils containing a pinhole from the pinhole turret of the scanner, as this produces one pinhole of extremely large size (R = 28 in the case of the Fluoview 300) whose radius is far from the initial sharp curvature due to the Gaussian distribution of direct back-propagating SHG, allowing the collection of a great deal of the diffuse forward-propagated signal.

3. Results

3.1 Determining backward propagating SHG spot size on the pinhole plane

In order to produce the required plot of SHG signal versus pinhole size and fit it with our model, we must first determine the sizes of the pinholes, relative to the Gaussian spot size of the direct backward propagating SHG on the pinhole plane. Since we are determining the spot size in the pinhole plane of the direct backward propagating SHG, we must avoid significant backscattering of the forward propagating SHG. Therefore rat tail collagen samples were fresh frozen and sectioned into 10um thin slices, spread out on a cover slip and dried overnight in a refrigerator for good adhesion between sample sections and cover slip. The cover slip was then flipped over with the sample side facing down, submerged in PBS and the excitation beam transmitted through the cover slip. The sample was submerged in saline to minimize the refractive index change and thus reduce the subsequent backscattering of the forward propagating SHG. The saline container is a cup with 4 cm diameter and 4 cm depth. It was painted black to absorb forward propagating SHG that goes through the sample section (see Fig. 5 ).

Fig. 5 Configuration of objective, sample and scatterer holder

To determine the size of the image of the SHG spot in the pinhole plane, we measured the SHG spot directly with a CCD camera (SPOT RT3, SciTech) at the intermediate image plane depicted on Fig. 1. The SHG spot in the pinhole plane is just a magnified version of the spot we captured with the CCD camera with a magnification factor of the ratio of the focal length of the pupil lens and collection lens.

As shown in Fig. 6(a) above, we drew a straight line through the center of the SHG spot captured by CCD camera, and measured SHG intensity along this line. The SHG intensity vs. horizontal location data set was then plotted in Fig. 6(b) as a series of separate points, and fitted to a Gaussian model producing an average spot size in the intermediate image plane of 9.34 pixels or 69 um. Taking into account the magnification factor of the pupil and collector lens, the average pinhole plane spot size of the backward propagating SHG is therefore 236.8 um. We already knew the absolute diameters of the five pinholes on the pinhole turret are 60um 100um 150um 200um to 7000um respectively (7000um corresponds to the location in the pinhole turret where we punched out the pinhole-bearing metal foil). That corresponds to 0.24 ω, 0.4 ω, 0.6 ω, 0.8 ω and 28 ω, where ω stands for the e−2 SHG spot size in the pinhole plane.

Fig. 6 Determination of SHG spot size. (a) CCD captured image of SHG focal spot after the tube lens (b) SHG intensity along the straight line that goes through the center of SHG spot. The red straight line in (a) is the line along which SHG intensities are measured. The red curve in (b) is the Gaussian fit to the intensities we measured in (a)

3.2 Measurement of SHG intensity varied with pinhole size

To evaluate our new method we chose the rat tail tendon, a sample whose SHG properties have been well studied [1

R. M. Williams, W. R. Zipfel, and W. W. Webb, “Interpreting second-harmonic generation images of collagen I fibrils,” Biophys. J. 88(2), 1377–1386 (2005). [CrossRef]

]. A whole rat tail was removed from a previously sacrificed animal (removal from the animal is not necessary for the technique but is convenient for handling the sample). To generate a clear SHG image of the tendon in particular (to match the sample previously studied in the literature) we peeled a thin layer of outer skin off the rat tail at the location of interest and exposed the tendon beneath it. We then put the rat tail on a glass slide, with the exposed collagen fiber facing up and we put another coverslip on top of the collagen fibers to ensure a stable sample. The rat tail and coverslip were then fixed on the glass slide with plastic tape and the collagen fibers were imaged though the coverslip. Note that, unlike previous methods which measure F/B with a second detector for the forward-propagating SHG, the tendons remained within the ~1cm diameter tail.

We prepared whole rat tail samples from 5 separate animals (see Fig. 7 ). On each rat tail we chose 5 image fields. And for each image field we took 5 back detected SHG images with the pinhole size varied from 60um 100um 150um 200um to 7000um, plus one image with no sample in order to quantify the background noise. One of these image sets is shown below.

Fig. 7 In vivo SHG images of rat tail collagen fibers in an intact rat tail. Images 1-5 are SHG images of the same ROI when the size of the pinhole varied from 60um, 100um, 150um, 200um, to 7000um, respectively. The bright spots are blue fluorescent polystyrene beads for calibration. Images are 600 um across.

3.3 Curve fitting and prediction of collagen SHG F/B ratio

In each set of images, we picked 5 small regions of interest (ROIs) around collagen fibers and 3 regions of interest around the calibration beads. ROIs were drawn around fibers that extended over at least ~100 um in the image plane, to ensure that fibers were close to perpendicular to the optical axis. A length 100 um in an optical section of ~12 um thickness (e−2 z diameter of the PSF) corresponds to a maximum angle in the object plane of ~7 degrees. Signal intensities in these collagen and beads ROIs were measured with ImageJ, background was subtracted, and signal intensities were normalized so that the maximum intensity measured with the largest pinhole was set to 1. The average relative SHG intensity of all 25 collagen ROIs (5 ROIs in each of 5 image stacks) and average relative TPEF intensity of all 15 beads ROIs in one rat tail vs. relative pinhole size plot is shown below. Note the separation between the collagen fiber curve and the bead curve due to the different F/B ratios of rat tail collagen and fluorescent calibration beads.

As described above, for one rat tail we measured relative intensity in 25 collagen ROIs and 15 bead ROIs. For each ROI 5 relative intensities were measured with size of the pinhole varying from the smallest to the largest. Average relative SHG intensity of all 25 collagen ROIs and average relative TPEF intensity of all 15 bead ROIs, as shown in Fig. 8 , were considered a set of data for one rat tail. We then fit this set of real measured data to the model given in Eq. (4) and calculated the rat tail collagen fiber SHG F/B ratio in one animal by eliminating C using Eq. (5). This data collection/averaging/curve fitting/calculation process was repeated 5 times to calculate rat tail collagen fiber SHG F/B ratio in 5 animals. The results are listed in Table 1 below.

Fig. 8 Epidetected total collagen SHG intensity (blue triangles), or total bead TPEF intensity (blue squares) vs. pinhole size, fit to the model given by Eq. (4) (red lines). The horizontal axis is pinhole size in units of backward propagating SHG spot size on the pinhole plane i.e. fraction of ω. The vertical axis is normalized SHG intensity. Note that the fitting curves identically equal 1 when the pinhole size reaches R = 28, whose data point is not shown, and that the intensity information for R = 28 provides the normalization value and hence is included in the overall fit.
Table 1  Results of in vivo measured rat tail collagen SHG F/B ratio
Animal IAnimal IIAnimal IIIAnimal IVAnimal VAverage
F/B ratio2.10 ± 0.191.81 ± 0.302.50 ± 0.341.34 ± 0.132.26 ± 0.232.00 ± 0.45

For each animal, means ± 95% confidence intervals are presented, with the data from 25 ROIs being pooled to produce a single curve fit for each animal. The average value for the 5 animals is then presented ± the standard deviation.

To verify the validity of our method using thick samples and only epi-detection, we also directly measured rat tail collagen SHG F/B ratio with both forward and backward detectors and with the sample sectioned in ~10um thin slices, as previously described [4

X. Han, R. M. Burke, M. L. Zettel, P. Tang, and E. B. Brown, “Second harmonic properties of tumor collagen: determining the structural relationship between reactive stroma and healthy stroma,” Opt. Express 16(3), 1846–1859 (2008). [CrossRef] [PubMed]

]. In summary, the SHG image captured in the forward detector was divided by the SHG image captured in the backward detector. We then averaged SHG F/B ratio over all pixels within collagen fibers and ignored all pixels outside the collagen fibers. The scattering of the tissue, though very low in this case, was corrected for by measuring the F/B ratio of the TPEF signals from 10um diameter calibration fluorescent beads.

The measurement results from 5 animals, using forward and backwards detectors, are listed in Table 2 .

Table 2  Results of in vitro measured rat tail collagen SHG F/B ratio
Animal VIAnimal VIIAnimal VIIIAnimal IXAnimal XAverage
F/B ratio1.43 ± 0.801.22 ± 0.272.29 ± 0.491.75 ± 0.502.11 ± 0.861.76 ± 0.45

For each animal, means ± standard deviations are presented from 5 F/B measurements. The average value for the 5 animals is then presented ± the standard deviation.

The overall average F/B ratio from 5 animals is 1.76+/−0.45. This result was not statistically significantly different from our measurement using only the epi-detection objective lens and the model given by Eq. (4) (Student’s t-test p = 0.42) demonstrating that the results from the two methods are in good agreement. To evaluate the new method at a higher NA we also measured the F/B ratio using a 0.8 NA lens, producing a value for F/B (1.37+/−0.28) which is again not statistically significantly different from the direct measurement using forward and backwards detectors (p = 0.17 N = 5).

4. Discussion

In this paper, we present a method to measure the SHG F/B ratio with only an epi-detection objective lens, allowing measurement on thick tissue samples. In our demonstration we used the whole rat tail (~1 cm thick), with the rat removed for convenience, but the technique is equally applicable to a live specimen. By obviating the need for thin sectioning of the sample, this technique provides the opportunity to do time-dependent studies, as well as the possibility of use in an endoscopic setting. This is significant because SHG F/B ratios have been shown to be of interest in discriminating skin with Osteogenesis Imperfecta [2

R. Lacomb, O. Nadiarnykh, and P. J. Campagnola, “Quantitative second harmonic generation imaging of the diseased state osteogenesis imperfecta: experiment and simulation,” Biophys. J. 94(11), 4504–4514 (2008). [CrossRef] [PubMed]

] from normal dermis [2

R. Lacomb, O. Nadiarnykh, and P. J. Campagnola, “Quantitative second harmonic generation imaging of the diseased state osteogenesis imperfecta: experiment and simulation,” Biophys. J. 94(11), 4504–4514 (2008). [CrossRef] [PubMed]

] and SHG F/B ratio measurements have been used to help determine the organization of fibrillar collagen in samples such as breast tumor models [4

X. Han, R. M. Burke, M. L. Zettel, P. Tang, and E. B. Brown, “Second harmonic properties of tumor collagen: determining the structural relationship between reactive stroma and healthy stroma,” Opt. Express 16(3), 1846–1859 (2008). [CrossRef] [PubMed]

], in rat tail tendon [1

R. M. Williams, W. R. Zipfel, and W. W. Webb, “Interpreting second-harmonic generation images of collagen I fibrils,” Biophys. J. 88(2), 1377–1386 (2005). [CrossRef]

], cellulose [3

O. Nadiarnykh, R. B. Lacomb, P. J. Campagnola, and W. A. Mohler, “Coherent and incoherent SHG in fibrillar cellulose matrices,” Opt. Express 15(6), 3348–3360 (2007). [CrossRef] [PubMed]

], mouse Achilles tendon [8

F. Légaré, C. Pfeffer, and B. R. Olsen, “The role of backscattering in SHG tissue imaging,” Biophys. J. 93(4), 1312–1320 (2007). [CrossRef] [PubMed]

], muscle fascia [8

F. Légaré, C. Pfeffer, and B. R. Olsen, “The role of backscattering in SHG tissue imaging,” Biophys. J. 93(4), 1312–1320 (2007). [CrossRef] [PubMed]

], and ovarian cancer [7

O. Nadiarnykh, R. B. LaComb, M. A. Brewer, and P. J. Campagnola, “Alterations of the extracellular matrix in ovarian cancer studied by second harmonic generation imaging microscopy,” BMC Cancer 10(1), 94 (2010). [CrossRef] [PubMed]

] (although see 4.2 Limitations, below). Furthermore, the confocal pinholes used in this technique are available in most commercial two photon microscopes, allowing the measurement of the SHG F/B ratio in intact specimens without addition of new equipment (assuming one is already using a two-photon microscope to study SHG) or extensive equipment modification (except for a single dichroic and filter inserted in the dichroic holder of the scanbox, and the punching out of one of the pinhole settings to produce one very large pinhole).

4.1 Sensitivity

In order to determine the radial intensity distribution of the forward propagating and subsequently backscattered SHG using a Monte Carlo simulation, the tissue scattering parameters we used in the simulation were µa = 0.7 cm−1 and µs = 150 cm−1, which are typical values for 405 nm light [10

W. F. Cheong, S. A. Prahl, and A. J. Welch, “A review of the optical properties of biological tissues,” IEEE J. Quantum Electron. 26(12), 2166–2185 (1990). [CrossRef]

]. However, the scattering properties of tissue vary over a large range among different organs throughout the body, and it is reasonable to wonder if this fitting model is good for all types of tissues from different organs. To determine this we varied the scattering parameters over a large range and repeated the Monte Carlo simulation with the results shown in Fig. 9 and Fig. 10 .

Fig. 9 Monte Carlo simulation of forward propagating SHG that is subsequently backscattered and reaches the object plane, with varying µs. Steady state radial distribution of the backscattered SHG photons when (a) µs = 100 cm−1 (b) µs = 150 cm−1 (c) µs = 200 cm−1 (d) µs = 250 cm−1 (e) µs = 300 cm−1 (f) µs = 350 cm−1 (g) µs = 400 cm−1 (h) µs = 450 cm−1 (i) µs = 500 cm−1. In each case µa = 0.7 cm−1
Fig. 10 Monte Carlo simulation of forward propagating SHG that is subsequently backscattered and reaches the object plane, with varying µa. (a) Steady state radial distribution of the backscattered SHG photons when (a) µa = 0.05 cm−1 (b) µa = 0.1 cm−1 (c) µa = 1 cm−1 (d) µa = 2 cm−1 (e) µa = 3 cm−1 (f) µa = 4 cm−1 (g) µa = 5 cm−1 (h) µa = 6 cm−1 (i) µa = 7 cm−1. In all three cases µs = 500 cm−1

We can see from the simulation results that when µa decreases or µs increases more forward propagating SHG photons reach the object plane after subsequent backscattering; and when µa increases or µs decreases less reach the object plane. However, over a range of scattering parameters in our Monte Carlo simulations that are comparable to values found at these wavelengths in the literature [10

W. F. Cheong, S. A. Prahl, and A. J. Welch, “A review of the optical properties of biological tissues,” IEEE J. Quantum Electron. 26(12), 2166–2185 (1990). [CrossRef]

], varying µa and µs does not alter the distribution of light exiting the tissue over the 0-50 um length scale, and it remains a constant, independent of radial position (a slope of zero is within the 95% confidence interval). We therefore conclude that within this small region of the object plane close to the optical axis, which is the region covered by our pinhole distribution, the radial intensity distribution of the backscattered SHG is not affected by a wide variety of tissue scattering parameters, and this fitting model is therefore applicable for tissue from diverse organs.

4.2 Limitations

This technique has two significant limitations. As shown in Fig. 4, when the F/B ratio increases, the separation between two SHG intensity vs. pinhole size curves decreases. Using a reasonable value for C [8

F. Légaré, C. Pfeffer, and B. R. Olsen, “The role of backscattering in SHG tissue imaging,” Biophys. J. 93(4), 1312–1320 (2007). [CrossRef] [PubMed]

], Fig. 4 shows that when the collagen fiber SHG F/B ratio is significantly more than ~5, the two curves are so close that they are not likely to be distinguishable assuming typical variation in experimental data (such as shown in Fig. 8). This implies that the new method might not be applicable for samples with collagen fiber SHG F/B ratios significantly higher than ~5. To illustrate this, we applied this method to measure the SHG F/B ratio of collagen fibers in intact mouse skeletal muscle fascia (with overlying skin removed to allow access). Using the new technique, we produced a value of the collagen fiber F/B ratio in vivo of 5.05+/−0.81 (N = 5) while after extracting the fascia, mounting it between two coverslips and directly measuring SHG F/B ratio with two objective lenses we produced a value of 5.56+/1.03 (N = 5), which was not statistically significantly different (p = 0.43). We also applied this method to measure the SHG F/B ratio of collagen fibers in 4T1 breast tumor models. According to our previous work the true collagen fiber SHG F/B ratio in 4T1 mouse breast tumor is ~30. We measured collagen fiber SHG F/B ratio in vivo in whole tumor samples from 5 different animals (25 ROIs per animal), and the result is 6.64+/−2.98 (data not shown). As predicted, as the F/B becomes significantly larger than ~5, the result becomes inaccurate and noisy, with a greater relative standard deviation, because of the normal random noise in the experimental data combined with the insensitivity of this new technique to large F/B values. Fortunately, our data and published studies reveal that the F/B ratio in many biological samples is at or below five, including in the rat tail tendon (F/B~1), skeletal muscle fascia (F/B~5), cellulose (F/B~4), ovarian cancer (F/B ~3.4), and dermis with (F/B~3.4) and without (F/B~2.64) Osteogenesis Imperfecta [1

R. M. Williams, W. R. Zipfel, and W. W. Webb, “Interpreting second-harmonic generation images of collagen I fibrils,” Biophys. J. 88(2), 1377–1386 (2005). [CrossRef]

–3

O. Nadiarnykh, R. B. Lacomb, P. J. Campagnola, and W. A. Mohler, “Coherent and incoherent SHG in fibrillar cellulose matrices,” Opt. Express 15(6), 3348–3360 (2007). [CrossRef] [PubMed]

,7

O. Nadiarnykh, R. B. LaComb, M. A. Brewer, and P. J. Campagnola, “Alterations of the extracellular matrix in ovarian cancer studied by second harmonic generation imaging microscopy,” BMC Cancer 10(1), 94 (2010). [CrossRef] [PubMed]

,8

F. Légaré, C. Pfeffer, and B. R. Olsen, “The role of backscattering in SHG tissue imaging,” Biophys. J. 93(4), 1312–1320 (2007). [CrossRef] [PubMed]

]. Care must be taken when this technique produces F/B ratios significantly greater than five.

The second limitation arises from the fact that this technique assumes there is no subsequent scattering of back-propagating SHG signal. Hence it is limited to the surface, or extremely shallow imaging depths, of the sample and is suitable for quantifying F/B ratios on intact thick tissue samples, but not in intact thick tissue samples.

5. Conclusion

In this paper, we present a method to measure SHG F/B ratio suitable for intact tissue samples without sectioning, using just the epidetection objective lens. The method requires minor modification of the dichroics, filters, and confocal pinholes used in most commercial two photon microscopes and hence minimizes the purchase of new equipment. This allows F/B ratio measurements to be done in a dynamic fashion, and offers the possibility of endoscopic measurements. This technique is sensitive to F/B values up to ~5 and within this range are a variety of interesting and clinically relevant materials, tissues, and disease states, including tendon, fascia, cellulose, ovarian cancer, and dermis with and without Osteogenesis Imperfecta. The fact that OIM is within the range of this technique’s sensitivity offers the attractive possibility of a non-invasive optical diagnosis of that disease, using a single objective lens, without biopsy samples being removed from the patient.

Acknowledgements

We thank Dr. Jerome Mertz for helpful conversations. This work is supported by a Department of Defense BCRP Pre-doctoral Traineeship Award (W81XWH-08-1-0323) to Xiaoxing Han, a Department of Defense BCRP Era of Hope Scholar Award (W81XWH05-1-0396), and a Pew Scholar in the Biomedical Sciences Award to Edward Brown III.

References and links

1.

R. M. Williams, W. R. Zipfel, and W. W. Webb, “Interpreting second-harmonic generation images of collagen I fibrils,” Biophys. J. 88(2), 1377–1386 (2005). [CrossRef]

2.

R. Lacomb, O. Nadiarnykh, and P. J. Campagnola, “Quantitative second harmonic generation imaging of the diseased state osteogenesis imperfecta: experiment and simulation,” Biophys. J. 94(11), 4504–4514 (2008). [CrossRef] [PubMed]

3.

O. Nadiarnykh, R. B. Lacomb, P. J. Campagnola, and W. A. Mohler, “Coherent and incoherent SHG in fibrillar cellulose matrices,” Opt. Express 15(6), 3348–3360 (2007). [CrossRef] [PubMed]

4.

X. Han, R. M. Burke, M. L. Zettel, P. Tang, and E. B. Brown, “Second harmonic properties of tumor collagen: determining the structural relationship between reactive stroma and healthy stroma,” Opt. Express 16(3), 1846–1859 (2008). [CrossRef] [PubMed]

5.

A. C. Kwan, D. A. Dombeck, and W. W. Webb, “Polarized microtubule arrays in apical dendrites and axons,” Proc. Natl. Acad. Sci. U.S.A. 105(32), 11370–11375 (2008). [CrossRef] [PubMed]

6.

S. V. Plotnikov, A. C. Millard, P. J. Campagnola, and W. A. Mohler, “Characterization of the myosin-based source for second-harmonic generation from muscle sarcomeres,” Biophys. J. 90(2), 693–703 (2006). [CrossRef]

7.

O. Nadiarnykh, R. B. LaComb, M. A. Brewer, and P. J. Campagnola, “Alterations of the extracellular matrix in ovarian cancer studied by second harmonic generation imaging microscopy,” BMC Cancer 10(1), 94 (2010). [CrossRef] [PubMed]

8.

F. Légaré, C. Pfeffer, and B. R. Olsen, “The role of backscattering in SHG tissue imaging,” Biophys. J. 93(4), 1312–1320 (2007). [CrossRef] [PubMed]

9.

L. Wang and S. L. Jacques, “Hybrid model of Monte Carlo simulation and diffusion theory for light reflectance by turbid media,” J. Opt. Soc. Am. A 10(8), 1746–1752 (1993). [CrossRef]

10.

W. F. Cheong, S. A. Prahl, and A. J. Welch, “A review of the optical properties of biological tissues,” IEEE J. Quantum Electron. 26(12), 2166–2185 (1990). [CrossRef]

11.

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

12.

W. R. Zipfel, R. M. Williams, and W. W. Webb, “Nonlinear magic: multiphoton microscopy in the biosciences,” Nat. Biotechnol. 21(11), 1369–1377 (2003). [CrossRef] [PubMed]

OCIS Codes
(170.0170) Medical optics and biotechnology : Medical optics and biotechnology
(180.0180) Microscopy : Microscopy

ToC Category:
Medical Optics and Biotechnology

History
Original Manuscript: February 23, 2010
Revised Manuscript: April 16, 2010
Manuscript Accepted: April 19, 2010
Published: May 5, 2010

Virtual Issues
Vol. 5, Iss. 9 Virtual Journal for Biomedical Optics

Citation
Xiaoxing Han and Edward Brown, "Measurement of the ratio of forward-propagating to back-propagating second harmonic signal using a single objective," Opt. Express 18, 10538-10550 (2010)
http://www.opticsinfobase.org/vjbo/abstract.cfm?URI=oe-18-10-10538


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. R. M. Williams, W. R. Zipfel, and W. W. Webb, “Interpreting second-harmonic generation images of collagen I fibrils,” Biophys. J. 88(2), 1377–1386 (2005). [CrossRef]
  2. R. Lacomb, O. Nadiarnykh, and P. J. Campagnola, “Quantitative second harmonic generation imaging of the diseased state osteogenesis imperfecta: experiment and simulation,” Biophys. J. 94(11), 4504–4514 (2008). [CrossRef] [PubMed]
  3. O. Nadiarnykh, R. B. Lacomb, P. J. Campagnola, and W. A. Mohler, “Coherent and incoherent SHG in fibrillar cellulose matrices,” Opt. Express 15(6), 3348–3360 (2007). [CrossRef] [PubMed]
  4. X. Han, R. M. Burke, M. L. Zettel, P. Tang, and E. B. Brown, “Second harmonic properties of tumor collagen: determining the structural relationship between reactive stroma and healthy stroma,” Opt. Express 16(3), 1846–1859 (2008). [CrossRef] [PubMed]
  5. A. C. Kwan, D. A. Dombeck, and W. W. Webb, “Polarized microtubule arrays in apical dendrites and axons,” Proc. Natl. Acad. Sci. U.S.A. 105(32), 11370–11375 (2008). [CrossRef] [PubMed]
  6. S. V. Plotnikov, A. C. Millard, P. J. Campagnola, and W. A. Mohler, “Characterization of the myosin-based source for second-harmonic generation from muscle sarcomeres,” Biophys. J. 90(2), 693–703 (2006). [CrossRef]
  7. O. Nadiarnykh, R. B. LaComb, M. A. Brewer, and P. J. Campagnola, “Alterations of the extracellular matrix in ovarian cancer studied by second harmonic generation imaging microscopy,” BMC Cancer 10(1), 94 (2010). [CrossRef] [PubMed]
  8. F. Légaré, C. Pfeffer, and B. R. Olsen, “The role of backscattering in SHG tissue imaging,” Biophys. J. 93(4), 1312–1320 (2007). [CrossRef] [PubMed]
  9. L. Wang and S. L. Jacques, “Hybrid model of Monte Carlo simulation and diffusion theory for light reflectance by turbid media,” J. Opt. Soc. Am. A 10(8), 1746–1752 (1993). [CrossRef]
  10. W. F. Cheong, S. A. Prahl, and A. J. Welch, “A review of the optical properties of biological tissues,” IEEE J. Quantum Electron. 26(12), 2166–2185 (1990). [CrossRef]
  11. T. J. Farrell, M. S. Patterson, and B. Wilson, “A diffusion theory model of spatially resolved, steady-state diffuse reflectance for the noninvasive determination of tissue optical properties in vivo,” Med. Phys. 19(4), 879–888 (1992). [CrossRef] [PubMed]
  12. W. R. Zipfel, R. M. Williams, and W. W. Webb, “Nonlinear magic: multiphoton microscopy in the biosciences,” Nat. Biotechnol. 21(11), 1369–1377 (2003). [CrossRef] [PubMed]

Cited By

 Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited