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

Label-free imaging of melanoma with nonlinear photothermal microscopy

Open Access Open Access

Abstract

Nonlinear photothermal microscopy, in which the intensity of the pump heating beam is modulated at f and the photothermal signal is extracted from the probe beam with a lock-in amplifier referred to 2f, is applied to the imaging of mouse melanoma without any staining. The pump and probe pulses, with central wavelengths of 488 and 632 nm, and a pulse duration of 100ps, are filtered from a compact commercial supercontinuum fiber laser source. An auto-balanced detector is applied to accumulate the signal and remove the laser noise of the probe. The spatial resolution of the nonlinear photothermal imaging is enhanced by 18% in both theoretical calculations and experiments, compared with a linear photothermal mechanism, and the resolution enhancement is theoretically 42% compared with conventional optical microscopy. This imaging technique shows possibilities for the clinical evaluation of melanoma with a high contrast and spatial resolution.

© 2015 Optical Society of America

Melanomas are among the most commonly occurring cancers and are a clinical challenge to diagnose [1]. Presently, the best method for the clinical evaluation of melanoma is visual inspection by dermoscopy; however, this technique is still far from a reliable diagnosis [2]. The optical diagnosis of melanoma is accessible because the two dominant types of melanin (eumelanin and pheomelanin) have a specific and intrinsic molecular contrast when imaged [3,4]. Several forms of nonfluorescent optical contrast mechanisms can be used to image melanoma without any staining, such as pump–probe [1,3], photoacoustic [58], and photothermal [7,8] microscopy. Photothermal microscopy has been used widely in material and biological research [711] because of the high sensitivity and lack of staining requirements. Compared with fluorescence microscopy, photothermal imaging does not suffer from photobleaching and blinking, which makes it attractive for robust single molecule detection [1214]. Besides the commonly used linear photothermal (LPT) mechanism, nonlinear photothermal (NLPT) microscopy has been used to image silica windows and gold pads with enhanced imaging contrast [15], in which the pump heating beam is chopped at f and the temperature variation of the sample is detected at 2f with a linear temperature sensor. Nonlinear phenomena have also been observed in photothermal-optical-deflection imaging experiments on samples of both high-purity aluminum and aluminum alloys [16]. The NLPT response of tungsten [17], thin solid films and coatings [18], carbon/epoxy composite materials [19], and cracks in solid [20] have also been investigated. The resolution of nonlinear thermal wave microscopes has been discussed, and it has been suggested that nonlinear microscopy is superior to linear microscopy [21]. However, to the best of our knowledge, all of the research on the effect of NLPT has been performed on bulk materials, thin films, and layered-structure samples; there are no experimental results about the NLPT imaging of biological tissues.

There are two possible explanations for the mechanism of the NLPT process. The first attributes the effects to thermal-wave second harmonic generation (SHG) related to the temperature-dependent heat capacity C(T) and thermal conductivity k(T) of the materials [22,23]. In the second explanation, a diffraction theory of continuous-wave photothermal deflection spectroscopy with fundamental and harmonic responses is presented, in which the displacement of the probe beam centroid leads to a set of analytical solutions to the fundamental and SHG, and the harmonics are caused by the diffraction of the probe beams in the mirage range [24]. In this Letter, we give a brief explanation of NLPT using the first mechanism. The conventional form of the heat conduction equation in a medium with temperature dependent heat capacity C(T) and thermal conductivity k(T) is as follows [25]:

CTt·(kT)=Q(r,t),
where Q is the external heat source input per unit volume. C(T) and k(T) can be expanded with the first Taylor expansion as follows:
C(T)C(T0)+[CT|T0]T1C0(1+δ1T1),
k(T)k(T0)+[kT|T0]T1k0(1+δ2T1),
where T0 is the initial temperature of the material, and T1 is the temperature change initiated by an external heating pump beam with an intensity modulated to a frequency of f. Then, Eq. (1) transforms to
(21D0t)T1=1k0Q(r,t)12(δ22δ1D0t)T12,
where D0=k0/C0 is the thermal diffusivity. Both the dependences C=C(T) and k=k(T) contribute to the nonlinear terms in the heat conduction equation. The quadratic nonlinearity in Eq. (4) induces the generation of the second harmonic of the thermal wave. With a one-dimensional simplification of the sample, the fundamental and second harmonic of the thermal wave is as follows [22]:
Tf(z,t)=Af(z)cos(ft+ϕf),
T2f(z,t)=A2f(z)cos(2ft+ϕ2f),
where Af(z) and ϕf are the amplitude and phase of the thermal wave with a fundamental frequency, and A2f(z) and ϕ2f are those for the second harmonic of the thermal wave. The amplitudes of the fundamental and second harmonic of the thermal wave have relationships with the pump intensity Af(z)Ipump and A2f(z)Ipump2, respectively. Assuming the response to be instantaneous, the refractive index caused by the temperature profile is Δn=T1n/T1 and acts as a modulated thermal lens. If Tf+T2f is used to replace T1, then we can obtain
Δn=TfnT1+T2fnT1=ΔnL+ΔnNL,
where ΔnL and ΔnNL are refractive indices induced by the fundamental and second harmonic thermal wave, respectively. The probe scattered by the thermal lens is proportional to Δn and the probe intensity Iprobe [26]. The NLPT signal can then be written as
SIpump2·Iprobe.

The relationship is also proved in [16] experimentally. Compared with the LPT signal, which can be written as SIpump·Iprobe [26], the cubic dependence of the NLPT signal on the laser intensity makes the spatial resolution of NLPT microscopy higher than that of the LPT one. If we assume that the pump and probe beams have the same wavelength and intensity distribution, this spatial resolution enhancement can be estimated to be 18% using scalar diffraction theory. Compared with conventional optical microscopy, the resolution enhancement of NLPT microscopy is 42%.

Figure 1 shows the schematic of the experimental setup. Both the pump and probe pulses, with a 10 nm spectral width, are spectrally filtered from a commercial supercontinuum (SC) fiber laser (WL-SC400-4, Fianium, Southampton, UK) using band-pass filters. The pulse duration of a 10 nm spectral width component of the SC is measured as 100ps and the repetition frequency is tunable in the range of 0.1–40 MHz. The central wavelengths of the pump and probe beams are 488 and 632 nm, respectively. After collinear combination by a short-pass dichroic mirror (DM) with a cutoff wavelength of 560 nm, the two beams are focused into the specimen by a microscope objective with an NA of 0.9 and amplification of 60×. The forward propagating beam is collected and collimated by a condenser lens (Olympus) with an NA tunable from 0 to 1.4. Then, the probe beam is spectrally filtered out by a band-pass filter with a central wavelength of 632 nm. To eliminate the effects of the laser noise of the probe beam, an auto-balanced detector (BD) (Nirvana 2007, Newport) is used to accumulate the signals. The intensity of the pump heating beam is modulated at f with an electro-optic modulator (EOM) (LM0202P, Qioptiq Photonics, Germany), creating an LPT signal at f and an NLPT signal at 2f, as depicted in Fig. 1(b). Two lock-in amplifiers (LIA) (SR844, Stanford Research System, US; Model 7265, Signal Recovery, US) are applied to extract the LPT and NLPT signals from the probe simultaneously, with reference frequencies of f and 2f, respectively. The sample is mounted on a set of piezo stages (P-622.2CL and P-622.ZCL, Physik Instrumente (PI), Germany), with a resolution of better than 1 nm in all three directions.

 figure: Fig. 1.

Fig. 1. (a) Schematic diagram of the experimental setup. EOM, electro-optic modulator; DM, short-wavelength transmitting dichroic mirror with a switching wavelength of 560 nm; BS, beamsplitter; MMF, multi-mode fiber; Filter, bandpass filter with a central wavelength of 632 nm and a bandwidth of 10 nm; BD, auto-balanced detector with a bandwidth of 125 kHz; LIA1,2, lock-in amplifiers. (b) Temporal intensity modulation behavior of the pump and probe pulse trains before and after interacting with the sample. The intensity of the pump pulse is modulated at f with the EOM and the output probe pulse is modulated at both f and f with a modulation depth of ΔIf and ΔI2f, respectively, because of the generation of the second-harmonic of the thermal wave. The LPT and NLPT signals are demodulated from the output probe beam with two LIAs referring to f and 2f, respectively.

Download Full Size | PDF

The sample of mouse skin melanoma is prepared as follows. B16 melanoma cells RCB-1283 (Riken Bio-Resource Center) are cultured in RPMI-1640 supplemented with 10% FBS, 100 units/ml of penicillin, and 100μg/ml of streptomycin (Meiji Seika Kaisha Ltd.) at 37°C in a humidified atmosphere of 95% air and 5% CO2. A 50 μl suspension containing 0.5 million melanoma cells for each head is subcutaneously inoculated into female nude mice (six weeks old) as five aliquots on the dorsal sides from the base of the tail to the neck. Tissues containing the inoculated cells were fixed with 4% paraformaldehyde for three days, embedded in melted paraffin, sliced at 15–20 μm, extended on glass slides, and enclosed with a coverslip.

First, we optimize the whole experimental setup and check the linearity of the detection system. The repetition frequency of the pump and probe pulses is set to its maximum value (40 MHz) for a high signal-to-noise ratio (SNR). The LPT and NLPT signals decrease with the increase of the modulation frequency of EOM, as is shown in Fig. 2(a). The optical noises of the pump and probe pulses increase sharply with the frequency decreases to <20kHz. The modulation frequency of the EOM is set to 30 kHz for the consideration of both the signal and the noise level. The auto-balanced detector is critical for imaging, as it can help to reduce the laser noise of probe pulse by 20–30 dB and obtain images with a higher SNR. The NA of the condenser lens is set within the range of 0.4–0.6 for the maximum LPT and NLPT signals. The time delay between the pump and probe pulses is also adjusted carefully for the highest LPT and NLPT signals (within 1 ns). After optimizing the experimental setup, we check the linearity of the detection system and ensure that the NLPT signal is not due to the saturation and nonlinearity of the detector. As shown in Fig. 2(b), both the LPT and NLPT signals increase linearly with the probe power and no obvious saturation is observed. The LPT signal also shows a linear response to the pump power, shown as a black line in Fig. 2(c). The best fit to the data of NLPT signal is as an equation of the form y=ax2+bx, shown as a blue line in Fig. 2(c). The term bx is because of imperfection of the square wave modulator (EOM in the experiment), which may have generated a small spurious NLPT signal. The corrected curve (without the term bx) for the response of the NLPT signal to the pump power is shown as a red line in Fig. 2(c). To show the fitting more clearly, Fig. 2(d) is shown with vertical (x) and longitudinal (y) axes in the logarithmic scale to the base 10 of pump power and NLPT signal, respectively. The slope of the fitting curve increases from 1 to 1.75 with the pump power increased from 0.4 to 1.0 mW, which means that the NLPT signal is mainly the result of the term bx with low pump power, and the effect of the term ax2 becomes significant at high pump rate.

 figure: Fig. 2.

Fig. 2. (a) Dependence of the LPT and NLPT signals on the modulation frequency of the EOM. Responses of the LPT and NLPT signals to (b) the probe power and (c) the pump power. The pump power is 1 mW in (b) and the probe power is set to 0.5 mW in (c). In (d), the x and y axes are the logarithm to the base 10 of pump power and NLPT signal, respectively. The slopes of the red, green, and blue parts of the fitting curve are 1.00, 1.52, and 1.75, respectively.

Download Full Size | PDF

The LPT and NLPT images of mouse melanoma are obtained simultaneously in the experiment. The pump and probe powers are set to 0.62 and 0.45 mW, respectively, which are in the linear range of the detection system shown in Fig. 2. The pump power is also far below the damage threshold of the sample. To obtain images with a high SNR, the time constant of the LIA and the dwell time of the imaging are set to 10 and 20 ms, respectively. The LPT and NLPT imaging of the mouse melanoma is then performed in a range of 10μm×10μm. The melanin in the melanoma is clearly imaged with SNR>10, as shown in Figs. 3(a) and 3(b). The NLPT signal is about 100 times smaller than the LPT signal, although the ratio is variable with different pump rate and different sample position. As shown in Fig. 3(c), the cross section at the same position [the red line in Figs. 3(a) and 3(b)] indicates that the resolution of LPT and NLPT can achieve 212 and 174 nm, respectively, with a resolution enhancement of 18% for the NLPT microscopy compared with the linear microscopy. The contrast is also higher for the NLPT imaging, as the obscure structures in Fig. 3(a) no longer exist in Fig. 3(b). The contrast enhancement can also be seen in Fig. 3(c), in which the maximum to minimum value of the NLPT curve is larger than that of the LPT curve. The maximum SNR of Fig. 3(b) is >150, indicating that a high-quality image is still possible, even if the dwell time is decreased 10 times [27].

 figure: Fig. 3.

Fig. 3. (a) LPT and (b) NLPT imaging of mouse melanoma with the range of 10μm×10μm and whole pixels of 300×300. The powers of the pump and probe pulses are 0.62 and 0.45 mW, respectively. The line profiles in (c) indicated by the red line in (a) and (b) show that the spatial resolution for the LPT and NLPT imaging are 212 and 174 nm, respectively. The NLPT signal is about 100 times smaller than the LPT signal.

Download Full Size | PDF

The high harmonic demodulation technology is also applied in saturated excitation (SAX) microscopy to obtain sub-diffraction limited resolution [28,29], although the saturated fluorescence in SAX microscopy may induce photobleaching of the fluorescent dye, which will not happen in the label-free NLPT microscopy. The sensitivity of photothermal microscopy is also higher compared with fluorescence microscopy, as discussed in the introduction.

In summary, we demonstrated the label-free imaging of mouse melanoma with both LPT and NLPT microscopy simultaneously. The spatial resolution and contrast of NLPT microscopy are higher than those of LPT microscopy. The pump and probe pulses are filtered from a supercontinuum fiber laser, which makes the experimental setup more compact and cheaper than the Ti:sapphire laser and OPA system. The technique opens further possibilities for the clinical evaluation of stainless melanoma and other bio-tissues with a high resolution and contrast.

References

1. T. E. Matthews, I. R. Piletic, M. A. Selim, M. J. Simpson, and W. S. Warren, Sci. Transl. Med. 71, 71ra15 (2011).

2. H. Kittler, H. Pehamberger, K. Wolff, and M. Binder, Lancet Oncol. 62, 751 (2010).

3. R. Marchesini, A. Bono, and M. Carrara, J. Biomed. Opt. 14, 014027 (2009). [CrossRef]  

4. G. Zonios, A. Dimou, M. Carrara, and R. Marchesini, Photochem. Photobiol. 86, 236 (2010). [CrossRef]  

5. E. I. Galanzha, E. V. Shashkov, P. M. Spring, J. Y. Suen, and V. P. Zharov, Cancer Res. 69, 7926 (2009). [CrossRef]  

6. M. Mehrmohammadi, S. J. Yoon, D. Yeager, and S. Y. Emelianov, Curr. Mol. Imaging 2, 89 (2013). [CrossRef]  

7. V. P. Zharov, Nat. Photonics 5, 110 (2011). [CrossRef]  

8. D. A. Nedosekin, M. A. Juratli, M. Sarimollaoglu, C. L. Moore, N. J. Rusch, M. S. Smeltzer, V. P. Zharov, and E. I. Galanzha, J. Biophotonics 6, 523 (2013). [CrossRef]  

9. D. Lasne, G. A. Blab, F. D. Giorgi, F. Ichas, B. Lounis, and L. Cognet, Opt. Express 15, 14184 (2007). [CrossRef]  

10. L. J. Inglehart, A. Broniaotowski, D. Fournier, A. C. Boccara, and F. Lepoutre, Appl. Phys. Lett. 56, 1749 (1990). [CrossRef]  

11. A. H. Sahraoui, C. Kolinsky, S. Delenclos, A. Daoudi, and J. M. Buisine, J. Appl. Phys. 82, 6209 (1997). [CrossRef]  

12. M. Tokeshi, M. Uchida, A. Hibara, T. Sawada, and T. Kitamori, Anal. Chem. 73, 2112 (2001). [CrossRef]  

13. S. Berciaud, L. Cognet, G. A. Blab, and B. Lounis, Phys. Rev. Lett. 93, 257402 (2004). [CrossRef]  

14. D. Lasne, G. A. Blab, S. Berciaud, M. Heine, L. Groc, D. Choquet, L. Cognet, and B. Lounis, Biophys. J. 91, 4598 (2006). [CrossRef]  

15. Y. N. Rajakarunanayake and H. K. Wickramasinghe, Appl. Phys. Lett. 48, 218 (1986). [CrossRef]  

16. G. C. Wetsel Jr. and J. B. Spicer, Can. J. Phys. 64, 1269 (1986). [CrossRef]  

17. A. Salnick, J. Opsal, A. Rosencwaig, and A. Mandelis, Solid State Commun. 114, 133 (2000). [CrossRef]  

18. V. Gusev, A. Mandelis, and R. Bleiss, Mat. Sci. Eng. B 26, 111 (1994). [CrossRef]  

19. G. Kalogiannakis, D. V. Hemelrijck, and C. Glorieux, NDT & E Int. 39, 236 (2006). [CrossRef]  

20. V. Gusev, A. Mandelis, and R. Bleiss, Appl. Phys. A 57, 229 (1993). [CrossRef]  

21. O. Dka, A. Mikls, and A. Lrincz, Appl. Phys. A 48, 415 (1989). [CrossRef]  

22. V. Gusev, A. Mandelis, and R. Bleiss, Int. J. Thermophys. 14, 321 (1993). [CrossRef]  

23. A. Mandelis, A. Salnick, J. Opsal, and A. Rosencwaig, J. Appl. Phys. 85, 1811 (1999). [CrossRef]  

24. J. Zhao, J. Shen, and C. Hu, Opt. Lett. 27, 1755 (2002). [CrossRef]  

25. H. S. Carslaw and J. C. Jaeger, Conduction of Heat in Solids2nd ed. (Oxford University, 1959).

26. A. Gaiduk, P. V. Ruijgrok, M. Yorulmaz, and M. Orrit, Chem. Sci. 1, 343 (2010). [CrossRef]  

27. J. F. de Boer, B. Cense, B. H. Park, M. C. Pierce, G. J. Tearney, and B. E. Bouma, Opt. Lett. 28, 2067 (2003). [CrossRef]  

28. K. Fujita, M. Kobayashi, S. Kawano, M. Yamanaka, and S. Kawata, Phys. Rev. Lett. 99, 228105 (2007). [CrossRef]  

29. S. Kawano, N. I. Smith, M. Yamanaka, S. Kawata, and K. Fujita, Appl. Phys. Express 4, 042401 (2011). [CrossRef]  

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 (3)

Fig. 1.
Fig. 1. (a) Schematic diagram of the experimental setup. EOM, electro-optic modulator; DM, short-wavelength transmitting dichroic mirror with a switching wavelength of 560 nm; BS, beamsplitter; MMF, multi-mode fiber; Filter, bandpass filter with a central wavelength of 632 nm and a bandwidth of 10 nm; BD, auto-balanced detector with a bandwidth of 125 kHz; LIA1,2, lock-in amplifiers. (b) Temporal intensity modulation behavior of the pump and probe pulse trains before and after interacting with the sample. The intensity of the pump pulse is modulated at f with the EOM and the output probe pulse is modulated at both f and f with a modulation depth of ΔIf and ΔI2f, respectively, because of the generation of the second-harmonic of the thermal wave. The LPT and NLPT signals are demodulated from the output probe beam with two LIAs referring to f and 2f, respectively.
Fig. 2.
Fig. 2. (a) Dependence of the LPT and NLPT signals on the modulation frequency of the EOM. Responses of the LPT and NLPT signals to (b) the probe power and (c) the pump power. The pump power is 1 mW in (b) and the probe power is set to 0.5 mW in (c). In (d), the x and y axes are the logarithm to the base 10 of pump power and NLPT signal, respectively. The slopes of the red, green, and blue parts of the fitting curve are 1.00, 1.52, and 1.75, respectively.
Fig. 3.
Fig. 3. (a) LPT and (b) NLPT imaging of mouse melanoma with the range of 10μm×10μm and whole pixels of 300×300. The powers of the pump and probe pulses are 0.62 and 0.45 mW, respectively. The line profiles in (c) indicated by the red line in (a) and (b) show that the spatial resolution for the LPT and NLPT imaging are 212 and 174 nm, respectively. The NLPT signal is about 100 times smaller than the LPT signal.

Equations (8)

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

CTt·(kT)=Q(r,t),
C(T)C(T0)+[CT|T0]T1C0(1+δ1T1),
k(T)k(T0)+[kT|T0]T1k0(1+δ2T1),
(21D0t)T1=1k0Q(r,t)12(δ22δ1D0t)T12,
Tf(z,t)=Af(z)cos(ft+ϕf),
T2f(z,t)=A2f(z)cos(2ft+ϕ2f),
Δn=TfnT1+T2fnT1=ΔnL+ΔnNL,
SIpump2·Iprobe.
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.