Whole-field Fluorescence Lifetime Imaging With Picosecond ...

1 downloads 0 Views 106KB Size Report
Instruments Ltd., Unit 9, Hall Farm Workshops, South Moreton, Didcot, Oxon. OX11 9AG, U.K. ..... electricity supply with little or no water-cooling. We have.
370

IEEE JOURNAL OF SELECTED TOPICS IN QUANTUM ELECTRONICS, VOL. 4, NO. 2, MARCH/APRIL 1998

Whole-Field Fluorescence Lifetime Imaging with Picosecond Resolution Using Ultrafast 10-kHz Solid-State Amplifier Technology Keith Dowling, Mark J. Dayel, Sam C. W. Hyde, J. C. Dainty, Member, IEEE, Paul M. W. French, Periklis Vourdas, M. John Lever, Anthony K. L. Dymoke-Bradshaw, Jonathon D. Hares, and Paul A. Kellett

Abstract— We report the development of a high temporal resolution whole-field fluorescence lifetime imaging system based on an ultrafast solid-state laser system and a time-gated image intensifier operating at up to 10 kHz. The temporal instrument response is 110 ps and we have imaged (environmentally perturbed) differences in fluorescence lifetime as small as 20 ps. Fluorophores exhibiting single- or double-exponential fluorescence decay profiles are routinely imaged and a near real-time update time of 3 s for the fluorescence lifetime map has been demonstrated using a modest personal computer. We also present provisional fluorescence lifetime images of tissue constituents. This fluorescence lifetime imaging technology is applicable to almost any optical instrument configuration and, when coupled with existing all-solid-state diode-pumped ultrafast laser technology, may yield a potentially inexpensive instrument for in vitro and in vivo biomedical imaging. Index Terms— Biomedical imaging, biomedical microscopy, biomedical optical imaging, fluorescence, laser applications, laser biomedical applications. Fig. 1. Experimental setup for fluorescence lifetime imaging.

I. INTRODUCTION

I

N RECENT YEARS, the advances in all-solid-state laser technology have prompted the investigation of optical imaging in the visible/near-infrared (NIR) spectral region as an alternative to traditional medical diagnostic techniques. The potential benefits of spectroscopic resolution (discrimination), low hazard and low cost make biomedical optical techniques increasingly attractive. Biomedical optical imaging is limited in its efficacy, however by the structural and chemical heterogeneity of biological tissue and the strong scattering which limits the ability to quantify variations in the optical properties of the tissue under investigation. A method commonly employed to increase the optical contrast between tissue types is to use fluorescent marker dyes (fluorophores) which can be designed such that they are selectively absorbed Manuscript received September 29, 1997. This work was supported by the U.K. Engineering and Physical Sciences Research Council (EPSRC) and the U.K. Defence Evaluation and Research Agency (DERA). K. Dowling, M. J. Dayel, S. C. W. Hyde, J. C. Dainty, and P. M. W. French are with the Femtosecond Optics Group, Physics Department, Imperial College, London SW7 2BZ, U.K. P. Vourdas and M. J. Lever are with the Center for Biological and Medical Systems, Imperial College, London SW7 2BY, U.K. A. K. L. Dymoke-Bradshaw, J. D. Hares, and P. A. Kellett are with Kentech Instruments Ltd., Unit 9, Hall Farm Workshops, South Moreton, Didcot, Oxon. OX11 9AG, U.K. Publisher Item Identifier S 1077-260X(98)03776-9.

in the specific area of tissue under investigation. In medical diagnostic techniques, the presence of a particular tissue type (e.g., cancerous tissue) can be established by detecting the emission-wavelength signature of the fluorophore (which will only be seen if the tissue of interest is present). This spectroscopic technique can be combined with optical imaging techniques to produce a “map” of the localization of the fluorophore and hence a map of the tissue under investigation. Problems arise firstly, because biological tissue exhibits spectrally and spatially heterogeneous absorption properties and secondly, it fluoresces itself upon optical excitation (autofluorescence). It is often necessary, therefore, to use a wavelength-ratiometric technique to detect and image a specific fluorophore. Unfortunately, there are not many wavelength-ratiometric fluorophore probes available and many require UV excitation which limits the optical penetration depth and exacerbates the problems associated with tissue autofluorescence. An alternative approach is fluorescence lifetime imaging (FLIM), in which the temporal decay of the fluorescence signal, rather than just its intensity, is measured. Fluorescence lifetime is a signature of a fluorophore that is relatively unaffected by the heterogeneous absorption properties of tissue,

1077–260X/98$10.00  1998 IEEE

DOWLING et al.: WHOLE-FIELD FLUORESCENCE LIFETIME IMAGING

371

(a)

(a)

(b)

(b)

Fig. 2. (a) FLIM map of interleaved Coumarin 314/DASPI samples, as discussed in text. (b) Cross section through center of FLIM map.

making it possible for measurements to be made through greater tissue depths using visible or near-infrared (NIR) excitation. Fluorescence lifetime measurements may also be spectrally resolved to provide more information about the fluorophores and to aid discrimination between them. In biochemical applications of fluorescence imaging, the effect of environment on the process of fluorescence may be used to map chemical or physical changes within a sample. The quantum efficiency of fluorescence is a function of the radiative and nonradiative decay rates. The radiative decay rate is considered constant for a given fluorophore, while the nonradiative decay rate can vary with environment. Unfortunately the quantum efficiency is not easy to determine as it is difficult to measure the exact quantity of fluorophore in a particular region, and to quantify how much pump-light is absorbed. Fluorescence lifetime, however, is also a function of fluorophore environment [1] and, since it may be determined using only relative intensity measurements, its determination does not require knowledge of the fluorophore concentration or excitation flux in the sample. Imaging fluorescence lifetime may therefore provide spatially resolved chemically specific (functional) data about a tissue sample under investigation. Fluorescence lifetime probes already exist for the measurement of, e.g., Ca concentration, O concentration and pH. Nonbiomedical applications of FLIM have also been demonstrated, including determination of impurities in metal samples for nuclear process control [2], and in combustion related studies [3].

Fig. 3. Lifetime dependence of DASPI on viscosity. From l ! r: 20=80– 40=60–60=40–80=20  glycerol/ethanol. (a) FLIM map. (b) Cross section through its center.

In the following sections, we describe the operation of our FLIM system and its performance is evaluated using convenient laser dye fluorophores. We then discuss its application to laser induced autofluorescence in biological tissue constituents and then briefly outline the scope for future developments of this technology.

II. MEASUREMENT OF FLUORESCENCE LIFETIME Fluorescence lifetime may be measured in the frequency domain or the time domain. The most common technique is in the frequency domain where the sample is illuminated with a sinusoidally modulated continuous-wave laser and the fluorescence lifetime determined from the phase change between the excitation and measured fluorescence modulation [4]. This work concerns the time domain in which an ultrashort light pulse is used to excite the fluorophore and the intensity of the fluorescence is then measured as a function of time. Previous demonstrations of this technique have used detectors which have exhibited temporal resolutions of a few nanoseconds e.g., [5]–[7]. We report a fluorescence lifetime imaging system based on a novel time-gated optical intensifier [8], [9] which allows simultaneous measurement of the fluorescence lifetime at all pixels in the field of view with a demonstrated temporal resolution significantly better than alternative FLIM techniques. The detector technology operates up to 10-kHz

372

IEEE JOURNAL OF SELECTED TOPICS IN QUANTUM ELECTRONICS, VOL. 4, NO. 2, MARCH/APRIL 1998

(a)

(b)

(c)

(d)

Fig. 4. (a) Conventional fluorescence microscope image of slide showing collagen, elastin and collagen samples. (b) Fluorescence intensity image recorded just after excitation. FLIM maps correspond to (c) 1 , and (d) 2 of a double exponential fit the fluorescence decays of each pixel. Note nonlinear time scale for 2 .

repetition rate and is well matched to our solid-state ultrafast oscillator-amplifier technology based on Cr:LiSAF [10]. The pulse energies required ( 1 J) have already been demonstrated from an all-solid-state diode-pumped version of this laser system [11] and so this technology promises to provide a more versatile, practical and considerably cheaper alternative to large-frame ultrafast lasers and streak cameras. Using the apparatus shown in Fig. 1, fluorescent samples were illuminated by 10-ps pulses of up to 1 J energy at 415 nm, at a repetition rate of 5 kHz. These were derived from a commercial ultrafast Ti:sapphire laser (Spectra-Physics Tsunami) and amplified in a home-built Cr:LiSAF regenerative amplifier whose output was tuneable from 800 to 880 nm (400–440 nm in the second-harmonic signal). An image of the fluorescent sample was relayed onto the cathode of the time-gated image intensifier (Kentech Instruments Ltd. Gated Optical Imager), for which the gate width was measured to be 110 ps (this measurement includes triggering jitter). The intensifier was triggered by the switching of the regenerative amplifier Pockels cell via an electronic delay (Stanford Research Systems DG535) to set its position relative to the arrival time of the excitation light. By recording images of the fluorescence at different delays after excitation, a temporal profile of the fluorescence intensity was obtained simultaneously for each point in the field of view. Exponentials were fitted to the decay profiles on a pixel by pixel basis, and the spatial distribution of decay time constants (lifetimes) displayed—a FLIM map. In principle, since we record the complete temporal history of the fluorescence, we can fit (using a least-squares algorithm) any order of exponential decay to the fluorescence profiles. We have so far demonstrated FLIM imaging of fluorophores exhibiting single and double exponential decays. The data was acquired, via an intensified charged-coupled device (CCD) camera,

to a personal computer (Cyrix P200 processor). The total time required to acquire data and calculate and display a FLIM map depended on the number of samples and typically took several minutes. For single exponential, and with lower spatial and temporal resolution, it was possible to record FLIM maps with an update time of only three seconds. Further development of the processing algorithms and computer hardware should yield FLIM update rates of less than one second. Fig. 2 shows a typical FLIM map, assuming single exponential decay profiles, of a fluorescent phantom consisting of three pipettes containing a Coumarin 314 solution in ethanol, separated by pipettes containing a DASPI solution in an ethanol/ethylene glycol mixture. The fluorescence lifetime M Coumarin 314 in ethanol was of a solution of 10 independently measured to be 3.45 ns using time-correlated photon counting. If we spatially average over a 1-mm area ( 2000 pixels) of the Coumarin 314 portion of our FLIM maps to derive a characteristic lifetime and average these values obtained for two measurements, each of two different M Coumarin 314 in ethanol, we obtain a samples of 10 fluorescence lifetime of 3.47 ns with a standard deviation of 0.02 ns. The fluorescence lifetime of DASPI is highly sensitive to the solvent properties, making calibration difficult, but an independent measurement of a pure ethanol solution of (10 M) DASPI using a streak camera has yielded a lifetime of 70 ps. Fig. 2 illustrates the excellent dynamic range of our FLIM system. The upper limit of the lifetimes that can be measured is constrained only by the repetition rate of the system, while the lower limit is determined by the instrument response function. Note that the system can detect and image differences in fluorescence lifetime shorter than 10 ps, provided that the fluorescence decays are measured over sufficient dynamic

DOWLING et al.: WHOLE-FIELD FLUORESCENCE LIFETIME IMAGING

373

range [12], although absolute lifetime measurements are not presently determined to better than 100 ps. This latter precision could be improved by using more sophisticated deconvolution techniques. III. TEMPORAL DISCRIMINATION To demonstrate the temporal resolution of the FLIM imaging system presented here, DASPI, a laser dye whose lifetime is a sensitive function of its solvent viscosity [13], was used to provide a fluorescence phantom with an adjustable lifetime. It was dissolved in ethanol mixed with different amounts of glycerol to provide solutions of variable viscosity. FLIM lifetime maps were thus made of a range of samples prepared with different viscosities. As can be seen in the FLIM maps shown in Fig. 3, the measured lifetimes varied from 155 ps, for the least viscous phantom (lhs), to 320 ps for the most viscous solution (rhs). Further improvements in resolution could be made by reducing the system response time. This has recently been reduced to 90 ps by triggering the detection system from the Tsunami oscillator output and thereby reducing the jitter. IV. IMAGING

OF

TISSUE CONSTITUENTS

Any practical biomedical application of FLIM, either in vivo or in vitro, must obviously take account of the autofluorescence of tissue itself. Indeed, it is possible to exploit endogenous fluorescence as a means to detect and study metabolic function and disease, e.g., [14]. To this end we have begun to characterize the fluorescence of biological tissue, in vitro, using the apparatus reported here. Fig. 4(a) shows a conventional fluorescence microscope image of a sample of collagen located between two samples of elastin (excited at 450–490 nm), which are from calcaneous tendon and aorta of rat. Fig. 4(b) shows the time-gated fluorescence intensity image (horizontal illumination) recorded directly after excitation and Fig. 4(c) and (d) shows the FLIM maps obtained by fitting a double exponential decay to the data, giving two lifetimes, and , for each region of interest. Fig. 5 shows the fluorescence decay curves of (a) elastin and (b) collagen, derived from the FLIM data of Fig. 4. These decay curves clearly fit a double exponential profile as shown by the plots of the residuals to single- and double-exponential fits, Fig. 5. Although this preliminary data has not yet been fully analyzed and so should be treated with caution, it appears that the FLIM profiles do provide a means to discriminate between different tissue constituents. Future work will try to establish how different techniques to fix and preserve tissue samples influence their fluorescence lifetime signatures. Before these measurements were made, these samples had been frozen for five weeks since their preparation, during which time some dehydration may have occurred. We are currently investigating the impact of different modes of sample preparation and sample history on the FLIM signature of biological tissue. Ultimately, this work will be directed toward FLIM investigations of diseased and healthy tissue in vivo in the hope of achieving a useful contrast.

(a)

(b) Fig. 5. Fluorescence decay profiles (spatially averaged over samples) of (a) elastin and (b) collagen extracted from freshly excised rat tissue. The upper graphs are fluorescence decay profiles with the calculated exponential fits. The dashed line is a fit to a single-exponential decay, with  1.08 ns for elastin and 0.74 ns for collagen. The solid line is a fit to a double exponential decay, with 1 285 ps, 2 2.16 ns for elastin and for collagen 1 257 ps, 2 1.71 ns. The lower traces show the residuals to each of the fits, again the dashed line represents the residuals from a single exponential fit (the initial value has been cropped to show detail) while the solid line is for a double exponential fit.

=

=

=

=

=

374

IEEE JOURNAL OF SELECTED TOPICS IN QUANTUM ELECTRONICS, VOL. 4, NO. 2, MARCH/APRIL 1998

V. CONCLUSION

AND

FUTURE WORK

We have demonstrated a fluorescence lifetime imaging system, with an excellent temporal dynamic range, based on a time-gated image intensifier and a solid-state regenerative oscillator/amplifier laser system. We have measured an instrument temporal response of 110 ps and have demonstrated the ability to detect (and image) fluorescence lifetime differences of 20 ps. The laser source for this FLIM system may be replaced by an all-solid-state diode-pumped oscillator amplifier system (already demonstrated [11]), potentially leading to a commercially viable instrument. It is hoped that this may be useful for in vitro and in vivo imaging of fluorophore location and environment for biomedical and other applications. Recently, we have also started to characterize the temporal signatures of the endogenous fluorophores of biological tissue and have imaged double exponential fluorescence decays for collagen and elastin. In the near future the all-solid-state diode-pumped oscillator amplifier system will be developed such that the apparatus may be operated from a conventional electricity supply with little or no water-cooling. We have already established that this system requires a total of only three pump diodes providing 500-mW power at 670 nm from 100- m-wide stripe facets. Anticipated reductions in cost and size compared to existing laser and time-resolved detection systems will permit this technology to be located outside laser laboratories and in hospitals/medical research centers. It should be understood that this whole-field time-gated detection system is applicable to almost any optical instrumentation. We intend to investigate operation in conjunction with both microscope and endoscope configurations. ACKNOWLEDGMENT The authors acknowledge invaluable discussions with R. Taylor. The assistance of J. Gray and G. Rumbles of the Chemistry Department, Imperial College, in carrying out time correlated single-photon counting measurements is gratefully acknowledged. K. Dowling acknowledges an EPSRC CASE studentship supported by Imperial Cancer Research at The Royal Marsden NHS Trust. M. J. Dayel acknowledges an EPSRC CASE studentship supported by Scientific Generics plc.

[7] R. Cubeddu, P. Taroni, and G. Valentini, “Time gated imaging system for tumour diagnosis,” Opt. Eng., vol. 32, pp. 320–325, 1993. [8] A. D. Scully, A. J. MacRobert, S. Botchway, P. O’Neill, A. W. Parker, R. B. Ostler, and D. Phillips, “Development of a laser-based fluorescence microscope with subnanosecond time resolution,” J. Fluorescence, vol. 6, pp. 119–125, 1996. [9] K. Dowling, S. C. W. Hyde, J. C. Dainty, P. M. W. French, and J. D. Hares, “2-D fluorescence lifetime imaging using a time-gated image intensifier,” Opt. Commun., vol. 135, pp. 27–31, 1997. [10] S. C. W. Hyde, N. P. Barry, R. Mellish, P. M. W. French, J. R. Taylor, C. J. van der Poel, and A. Valster, “Argon-ion-pumped and diode pumped all solid state femtosecond Cr:LiSrAlF6 regenerative amplifier,” Opt. Lett. vol. 20, pp. 160–162, 1995. [11] R. Mellish, N. P. Barry, S. C. W. Hyde, R. Jones, P. M. W. French, J. R. Taylor, C. J. van der Poel, and A. Valster, “Diode-pumped Cr:LiSAF all-solid-state femtosecond oscillator and regenerative amplifier,” Opt. Lett., vol. 20, pp. 2312–2314, 1995. [12] T. Oida, Y. Sako, and A. Kusumi, “Fluorescence lifetime imaging microscopy (flimscopy),” Biophys. J., vol. 64, pp. 676–685, 1993. [13] W. Sibbett and J. R. Taylor, “Synchroscan streak camera study of potential saturable absorbers in the blue spectral region,” J. Luminescence, vol. 28, pp. 367–375, 1983. [14] B. B. Das, F. Liu, and R. R. Alfano, “Time-resolved fluorescence and photon migration studies in biomedical and model random media,” Rep. Prog. Phys., vol. 60, pp. 227–292, 1997.

Keith Dowling was born in Dublin, Ireland, in 1973. He received the B.Sc. degree in physics from Imperial College, London University, in 1994 and is working toward the Ph.D. degree in the application of ultrafast lasers to fluorescence lifetime imaging in the Physics Department at Imperial College.

Mark J. Dayel was born in 1972. He received the B.Sc. degree in physics in 1994 and the M.Sc. degree in engineering and physical science in medicine, from which he graduated with the Ash Prize, in 1996, both from Imperial College, London University. He started the Ph.D. degree in biomedical applications of fluorescent lifetime imaging in the Physics Department at Imperial College in 1996. He is currently working in the University of California at San Francisco.

Sam C. W. Hyde, for a biography, see this issue, p. 192.

REFERENCES [1] J. R. Lakowicz, Principles of Fluorescence Spectroscopy. New York: Plenum, 1983. [2] C. Moulin, P. Decambox, and L. Trecani, “Direct and fast uranium determination in zirconium by time-resolved laser-induced fluorescence spectrometry,” Analytica Chimica Acta, vol. 321, no. 1, pp. 121–126, 1996. [3] T. Q. Ni and L. A. Melton, “2-dimensional gas-phase temperaturemeasurements using fluorescence lifetime imaging,” Appl. Spectrosc., vol. 50, no. 9, pp. 1112–1116, 1996. [4] H. Szmacinski, J. R. Lakowicz, and M. L. Johnson, “Fluorescence lifetime imaging microscopy: Homodyne technique using high-speed gated image intensifier,” Methods in Enzymol., vol. 240, pp. 723–748, 1994. [5] X. F. Wang, T. Uchida, D. M. Coleman, and S. Minami, “A 2dimensional fluorescence lifetime imaging-system using a gated image intensifier,” Appl. Spectrosc., vol. 45, pp. 360–366, 1991. [6] T. Oida, Y. Sako, and A. Kusumi, “Fluorescence lifetime imaging microscopy (flimscopy)—Methodology development and application to studies of endosome fusion in single cells,” Biophys. J., vol. 64, pp. 676–685, 1993.

J. C. Dainty (M’83), for a biography, see this issue, p. 369.

Paul M. W. French, for a biography, see this issue, p. 192.

Periklis Vourdas was born in Greece in 1971. He studied electronics and computing at the Aristotle University of Thessaloniki, Greece, and received the B.Sc. degree in 1996. He joined Imperial College in 1996 and has recently completed the M.Sc. degree in engineering and physical science in medicine. During his time at Aristotle University, he spent an extended period in Grenoble, France, studying cluster formation by pulsed laser activation.

DOWLING et al.: WHOLE-FIELD FLUORESCENCE LIFETIME IMAGING

M. John Lever was born in Kenton, U.K., in 1943. After graduating in chemistry with chemical pharmacology at the University of Oxford, Oxford. U.K., he remained there to undertake a D.Phil. degree in hyperbaric medicine. Following Post-Doctoral fellowships at the Universities of Pennsylvania, Philadelphia, and Essex, U.K., he joined Imperial College where he is currently Reader in Physiological Mechanics in the Department of Biological & Medical Systems. For many years, his major research interests have been in mechanical behavior of body tissues, particularly in relation to cardiovascular disease. Imaging the transport of atherogenic materials through blood vessel walls led to collaboration with the femtosecond optics group in the Physics Department with the aim of developing novel methods for differentiating tissues and their properties. He is currently Honorary Secretary of the British Microcirculation Society and on the editorial board of the journal, Heart & Vessels. Dr. Lever is a member of the Physiological Society, The European Society for Microcirculation, and the European Vascular Biology Association.

Anthony K. L. Dymoke-Bradshaw was born in London, U.K., in 1952. He received the B.Sc. and Ph.D. degrees, both from Imperial College, University of London, in 1973 and 1980 respectively. In the 1980’s, he worked at Imperial College and Rutherford Appleton Laboratory on thermal conductivity in plasmas and developed techniques for using streak cameras for Thompson scattering. This led to work on plasma beat wave experiments, which eventually demonstrated the production of relativistic beat waves in plasmas. In 1983, he was a Founding Director of Kentech Instruments Ltd., which makes high-voltage fast-pulse generators and several imagers that use this technology. His current interests include the application of fast imaging systems to biomedical imaging. In 1993, Dr. Dymoke-Bradshaw was awarded the Paterson medal and prize jointly with J. Hares.

375

Jonathon D. Hares was born in Moreton-in-Marsh, U.K., in 1955. He received the Ph.D. degree in physics from Imperial College, London University in 1979 and has many years of experience in high speed instrumentation including that used with plasma physics, broadband radar, high-power lasers, lightning simulation, and night vision. He is the technical director of Kentech Instruments Ltd., U.K. Dr. Hares was awarded the Institute of Physics Paterson medal and prize at the Institute of Directors (jointly with A. K. L. Dymoke-Bradshaw) in 1993 for highly innovative design and successful marketing of a high-speed X-ray detector.

Paul A. Kellett was born in Stockport, U.K., in 1958. He received the B.Sc. degree in physics from Oxford University, Oxford, U.K. In his previous employment at Protech he was the development manager and was responsible for the design of computer based systems used for process control, including such sensitive applications as nuclear power stations. He is the Managing Director of Kentech Instruments Ltd., U.K. When he joined Kentech, he was able to combine his electronics experience with his physics training and has designed many of the complex, high-speed instrumentation systems, which have gained Kentech its worldwide reputation.