Open Access
1 January 2008 Visualization of fast-moving cells in vivo using digital holographic video microscopy
Hongyue Sun, Bing Song, Hongpai Dong, Brian Reid, Michael A. Player, John Watson, Min Zhao
Author Affiliations +
Abstract
Digital in-line holography offers some significant advantages over conventional optical holography and microscopy to image biological specimens. By combining holography with digital video microscopy, an in-line holographic video microscope is developed and is capable of recording spatial 3D holographic images of biological specimens, while preserving the time dimension. The system enables high-speed video recording of fast cell movement, such as the rapid movement of blood cells in the blood stream in vivo. This capability is demonstrated with observations of fast 3-D movement of live cells in suspension cultures in response to a gentle shake to the Petri dish. The experimental and numerical procedures are incorporated with a fast reconstruction algorithm for reconstruction of holographic video frames at various planes (z axis) from the hologram and along the time axis. The current system enables both lateral and longitudinal resolutions down to a few micrometers. Postreconstruction processing of background subtraction is utilized to eliminate noise caused by scattered light, thereby enabling visualization of, for example, blood streams of live Xenopos tadpoles. The combination of digital holography and microscopy offers unique advantages for imaging of fast moving cells and other biological particles in three dimensions in vivo with high spatial and temporal resolution.

1.

Introduction

Recently, significant advances have been made in the technologies for obtaining three-dimensional (3-D) information from imaging biological specimens in two dimensions. Most cell behavioral studies have been carried out with in vitro two-dimensional (2-D) culture models with only a few providing 3-D information. It becomes more and more obvious that in some instances cellular behavior in three dimensions may be significantly different from that in 2-D cultures.1 Various methods have been used to observe and record cell behavior in three dimensions. For instance, migration of green-fluorescence-protein-labeled cells has been studied2 in vivo and tumor cell motility in 3-D matrices that involve “rho signaling” from that in 2-D cultures.3, 4 Using a Hoffman microscope, the in situ behavior of cornea cells is demonstrated.5, 6 However, only relatively slow cell migration can be imaged by these methods and they are generally unsuited to study fast moving cells such as those in blood streams. Although the latter have been studied with other imaging systems, for instance, using a laser scanning ophthalmoscope,7 3-D spatial information was not obtained. When cells move out of the focal plane with high speed, it is impossible to trace them with currently available imaging techniques.

Holography provides biomedical scientists with a unique tool for imaging and analyzing the behavior of living cells in three dimensions in situ. A hologram records the full 3-D information of a scene while retaining the full parallax and perspective of the original scene. It does this on a single planar (2-D) recording medium such as photographic film or an electronic photosensor. The technique is especially attractive for 3-D particle tracking because a 3-D spatial volume can be reconstructed from a single 2-D hologram. If the recording is made using digital techniques, i.e., by recording on an electronic sensor, such as a CCD (charged-coupled device) or CMOS (complementary metal oxide semiconductor) device, and reconstructed numerically within a computer, then holographic videos can be made, thereby retaining the fourth dimension, namely, time. Digital holographic microscopy has a great potential for fast 3-D image acquisition. This has been demonstrated experimentally in a number of papers. For example, in-line digital holography has been used to visualize particle fields 8, 9, 10, 11 in aerosols and biological samples or as a noninvasive contrast imaging tool for observing marine biological samples in 3-D space.12, 13, 14

Recording holograms on electronic sensors, however, limits the resolution attainable to around 6 to 10μm , depending on the pixel size of the sensor, whereas on photographic film about 2μm is possible under favorable conditions. The limitation of spatial resolution can be overcome, to a certain extent, by using a divergent laser beam to illuminate the sample volume and record the hologram on a lensless camera.15 Since the depth of field of conventional microscopy is very limited due to the high numerical aperture necessary for high resolution, digital holographic microscopy (DHM) offers the possibility of focusing on different object layers by numerical reconstruction methods while still preserving high resolution. This technique also proved to be successful for 4-D (space and time) tracking16, 17 of microstructures and organisms using a lensless microscopic system. Algae motion is tracked in a 3-D space along the time coordinate by applying prereconstruction processing, i.e., the sum of subtracted holograms of each of two consecutive video frames. This was successful when dealing with suspended particles with low population density that were moving slowly. Several recent studies18, 19, 20 utilized a Mach-Zehnder interferometer configuration (with a separated reference beam) with a coherent or partially coherent light source to record intensity and phase-contrast digital microscopic holograms of a 3-D object on a glass slide. The authors successfully demonstrated the ability of DHM for resolving a microstructure of the phase object at high spatial resolution, often down to hundreds of nanometers. Sheng 21 successfully increased the depth of field and imaging resolution by combining holography with microscopy. This study mainly focused on examining the spatial distribution and velocity of a dense cloud of particles with an extended depth. Tracking cell movement with high speed in three dimensions in vivo has, however, not been attempted.

To observe fast-moving cells in 3-D space in animal tissues in vivo, we introduce a digital holographic video microscope (DHVM) with the use of coherent light. This achieves holographic imaging of biological specimens in vitro and in vivo over the 3-D volumetric space and in time. In this paper, we describe our DHVM system and demonstrate cell motions in different media. We also demonstrate that fast cell motion in blood streams in vivo can be recorded, visualized, and quantified with our system. The system enables both lateral and depth resolutions at the micrometer level. The experimental and numerical procedures have been incorporated into a program package with a fast reconstruction algorithm, which is capable of reconstructing of holographic video frames at various planes ( z axis) from the hologram and in the time dimension after image acquisition. This capability is interpreted in 3-D mappings, showing the motion of live cells (neutrophils) in suspension cultures and fast movement of blood cells in 3-D blood vessels of a live Xenopos tadpole. We also apply postreconstruction processing of background subtraction to eliminate background noise caused by the scattered light from dense cell population density in live tadpole tissue.

2.

Principle and Implementation of DHVM

A digital/electronic hologram (ehologram) is formed by the interference between waves scattered from an object on coherent illumination and the straight-through wave at the hologram/sensor plane. We took a simplified technical approach to buildup our system. The system [Fig. 1a ] utilizes a single-mode diode laser (continuous wave, 667-nm wavelength, 10-mW power) that delivers the collimated laser beam to illuminate the specimen through an optical fiber attached with a collimator. The microscope generates a magnified real image of the specimen in the space in front of a CCD (Sony ICX-082AL-6) sensor implemented in a Hitachi camera (KP-MA type) together with an associated image grabber (ViewCast Corporation Osprey 210). The sensor dimensions are 8.72mm (horizontal) by 6.52 (vertical) mm with the pixel pitch of 11.6μm (horizontal) by 11.2μm (vertical). The magnified real image and the directly transmitted light through the specimen form an interference pattern in space, i.e., a hologram. So a hologram is observed by varying the observation plane from the in-focus image plane. The interference pattern will then be captured by the CCD imaging sensor at the sensor or hologram plane (x,y) . Such a recorded 2-D ehologram, represented by E(x,y,0) , where (x,y) are the position coordinates of an arbitrary point P1 at the ehologram plane [Fig. 1b], stores the interfering fringe patterns and is postprocessed digitally in a computer. The CCD sensor records the 2-D spatial sampling on the recorded ehologram, which can be implemented for numerical reconstruction.

Fig. 1

Schematics of the DHVM system: (a) the DHVM system and (b) reconstruction scheme of an ehologram.

014007_1_026801jbo1.jpg

Digital hologram reconstruction has been widely implemented recently in a diverse range of scientific and engineering problems. Reconstructing an ehologram involves tracing the beam propagation through the system and calculating the diffraction pattern by means of a numerical representation of the physical process of diffraction. Many algorithms have been employed based on, for example, the Fresnel transform,22, 23, 24 the Fraunhofer transform,25 the wavelet transform,26 and “Fresnelets” with wavelet bases.27 The reconstruction algorithm used in our system is derived from the angular spectrum theory.28 A parallel reconstruction wave propagates through the 2-D ehologram aperture at z=0 [Fig. 1b], the wave is diffracted by the fringe pattern on the ehologram to form the reconstructed field at the observation plane (x,y) . The reconstructed field at an arbitrary point P2 in (x,y) plane is E(x,y,z) . The mathematical representation for 3-D diffraction of the angular spectrum is

E(x,y,z)=F1(Fh(x,y,0)exp{izk[1(qk)2]12}),
where q=qx2+qy2 is the 2-D variable in the Fourier domain and the conjugate of the variable (x,y) in the spatial domain, E(x,y,z) is the reconstructed complex amplitude of the diffraction wavefront in the (x,y) plane at a z distance from the hologram plane. Here Fh(x,y,0) represents the Fourier transform (FT) of the recorded wave interfering distribution at the hologram plane (x,y) at a point P1, k=2πλ , and λ the wavelength of the illuminating light. This is then rewritten in its integral form after taking into account the Fresnel region approximation28 as

Eq. 1

E(x,y,z)=iλzexp(i2πλz)×E(x,y,0)exp{iπλz[(xx)2+(yy)2]}dxdy.
Equation 1 is rearranged as a 2-D convolution:

Eq. 2

E(x,y,z)=E(x,y,0)iπτ2exp[i2(πτλ)2]K(xτ,yτ),
where xτ=xτ ; yτ=yτ ; τ2=λzπ , as the scaled distance parameter; and K(xτ,yτ)=exp[i(xτ2+yτ2)] is the chirp function and is considered as the mother wavelet. Function K(xτ,yτ) in Eq. 2 is not limited in both spatial and frequency domains, it has the scale and shift features of a wavelet family for 2-D function spaces. Because reconstruction involves numerical convolution of the kernel with the hologram intensity, it is necessary to avoid aliasing in the sampled version of the kernel. To satisfy the wavelet admissibility condition a Gaussian envelope function G(x,y)=exp[(x2+y2)δ]2 was introduced by Buraga-Lefebvre 29 to the mother wavelet. Thus, the Fresnel transform in Eq. 1 can be expressed by a wavelet transform:

Eq. 3

E(x,y,z)=E(x,y,0)k(τ)K(τ,p,q),
where k(τ)=(1πτ2)exp{i[2(πτλ)2(π2)]} is the amplitude of the daughter function, K(τ,p,q)=K[(xx)τ,(yy)τ]G(x,y) is the wavelet set; and p=(xx)τ and q=(yy)τ . In our DHVM system, the reconstruction algorithm is based on the wavelet transform in Eq. 3. Its kernel is based on that used by Onural and Kocatepe.30 However, we adjusted the Gaussian envelope function to provide the necessary filtering, by taking δ=πτ2A , with a parameter A selected to give a suitable reduction in the envelope amplitude at a value of (x2+y2) such that the wavelet frequency equals the Nyquist frequency. More detailed discussions on the performance of various algorithms can be found from Dong 31 It is readily verified that for pixel spacing d the envelope amplitude, at the Nyquist frequency, δ=exp(Ad2) , is independent of the reconstruction distance z . However, due to the usage of objective lenses through the microscopy system and the in-air wavelength in the reconstruction algorithm, the actual physical distance of the specimen image will not be identical to the reconstruction distance even if a parallel recording beam is used. Therefore, to make a valid measurement of the specimen, it is essential to calibrate the system transversely (along x and y ) and longitudinally (along z ). The transverse scale is calibrated by using a recorded hologram reconstruction of a standard scale slide for microscopy. Similarly, the longitudinal calibration is carried out by recording a hologram of an object with known thickness (dg) and reconstructing this hologram to find the z interval (in micrometers) of dg ; in our case, a 0.60-mm glass slide was used. This method will inevitably introduce deviations in z distance calculation, mainly because the refractive index (RI) of the cell solution (nc) differs from the RI of the glass (ng) , such that in a given time interval, the wave travels different distances in the cell solution (dc) and the glass. Ideally, the reconstruction program should take into account of the wavelength change when simulating wave propagation in different media. For simplicity, the algorithm presented here uses only one wavelength (667nm) to calibrate the longitudinal distance. The inaccuracy (Δd) is estimated as less than 10% of the distance traveled in dc , by considering that in a given time interval, the wave optical paths through all media are the same: dgng=dsolnsol (e.g., ng1.50 , nc1.35 , Δddc=0.1 ). The scattering effect was not considered in reconstruction for error estimation because the cell size concerned was larger than the laser wavelength (667nm) . However, when imaging submicrometer objects and the diffraction limit is reached, care should be taken in configuring the recording system.

3.

Results and Examples

3.1.

Living Cells in Suspension Culture

Adherent cells on a planar surface have been successfully imaged with holographic techniques in several laboratories.15 In addition to the adherent cells on flat surfaces, with our system cells in suspension can also be recorded. Here we present holograms and reconstructed images of two types of human cells, live umbilical vein endothelial cells and neutrophils, in suspension in culture medium.

Human umbilical vein endothelial cells were collected shortly after trypsin treatment to dislodge them from adhesion. This treatment released the cells from the planar surface into suspension status in a 3-D volumetric space. A typical hologram interference pattern [Fig. 2a ] and its reconstructed images at different planes demonstrates the viability of such a holographic system in imaging cells in suspension culture, i.e., in a 3-D volume. Figures 2b to 2h show a sequence of images reconstructed numerically from the hologram [Fig. 2a, also see Video 1 ] at reconstruction distances z of 100, 115, 120, 135, 150, 179, and 189μm , respectively [Figs. 2b, 2c, 2d, 2e, 2f, 2g, 2h], where the two arrowed cells are in focus at z=120μm . A 2-D in-focus micrograph of the same group of cells is presented [Fig. 2i] for comparison with the holographic images. Details of cell morphology and contrast in the in-focus holographic image appear comparable to or better than the in-focus micrographic image. Holographic images can be reconstructed at any z plane and offer detailed 3-D volumetric information. We notice that two bright dots, indicated by a white arrow in Figs. 2g and 2h, are present in the reconstruction domain at z=179 and z=189μm . These are formed by the sphere-shaped cells [arrowed in Fig. 2d, respectively], which act as small positive lenses, focusing the collimated illumination beam on their focal points, i.e., the bright dots. One can accurately measure the distance between the in focus reconstruction image [e.g., Fig. 2d] and the bright spots along z . The measured distance can be used to obtain cell property information, such as the refractive index of cell (if the cell diameter is known) or the dimension of the cell along z (if its refractive index is known). The distance between the spherical centers of these cells and their focal points [i.e., the effective focal length (EFL) of the cell lens], are obtained: EFL=29 and 69μm for the bottom-left and top-right arrowed cells, respectively, in Fig. 2d.

Fig. 2

Trypsinized human umbilical vein endothelial cells in suspension culture visualized by the DHVM: (a) hologram at z=0 plane; (b) to (f) reconstructions of (a) at z=100 , 115, 120 (two arrowed cells are on focus), 135, and 150μm , respectively; (g) and (h) reconstructions of the focal points formed by the arrowed spherical cells in (d); (i) micrograph for the same cells, arrowhead pointing at the on-focus cells; scale bar is 10μm .

014007_1_026801jbo2.jpg

Video 1

Reconstructed hologram sequence of human umbilical vein endothelial cells in z axis. Steps of z 30, step interval; 5.5μm ; whole frame size, 105×88μm (QuickTime, 3.0MB ). 1

014007_1_026801jbov1.jpg
10.1117/1.2841050.1

3.2.

Rapid Cell Motion in 3-D Suspension Culture

Cells can move rapidly when they are grown in a volume of suspension culture. It is impossible to track many such fast-moving cells simultaneously with conventional microscopy because one can only choose one focal plane at one time point. In DHVM recording, however, it is the interference patterns of all cells in the light path that a hologram stores, so that choosing a focal plane is no longer an issue. The original scenes of moving cells in the recording volume are obtained and extracted to show, for instance, the image focal plane and the 3-D spatial properties in the reconstruction stage. By reconstructing an electronic holographic video (eholovideo) in time sequence, fast-moving cells in a volume can be recorded on an eholovideo and then reconstructed for visualization and quantitative determination of the 3-D distribution over the recording period. To observe cell motion in free space, we chose small nonadherent cell lines (HL60 cells) in suspension culture as our volumetric recording object. HL60 cells are from the promyelocytic leukemia cell line, which normally grow nonadherently in suspension in culture medium. A gentle shake to the Petri dish caused fast movement of the cells in a volume. The movement was recorded in a 3-s eholovideo (original video not shown) using the system shown in Fig. 1a. The eholovideo was recorded at a rate of 10fps (frames per second). A reconstruction video at a certain z distance 698μm against the time axis is shown in Video 2 . In this reconstruction, one individual cell (arrowed) was shown in focus in some of the video frames at the chosen z distance at various time points, but was not visible on the adjacent reconstructed frames, indicating that cells moved in and out of the chosen z plane as well as in the x and y directions (Video 2) along the time axis. By reconstructing each hologram frame in the video in time sequence, the exact locations of moving cells are extracted; therefore, the fast movement of cells in three dimensions can be visualized and quantified successfully. Some eholovideo frames and their reconstructions at the chosen z distance of 698μm are selected to show one individual cell reconstruction at the indicated time points, marked in Figs. 3a, 3b, 3c , corresponding to the movement of the selected HL60 cell in 3s , shown in the 3-D mapping in Figs. 3d, 3e, 3f, respectively. A video of the 3-D trajectory for this selected cell was generated (Video 3 ) in MATLAB to further illustrate this motion schematically. According to the reconstruction coordinates (x,y,z) at each time point, the motion speed is estimated to be 260±26μms1 on average. Imaging such rapid cell motion is almost impossible with other imaging systems currently available, because the speed at which cells move in and out of the observation plane in a volume is too high to be focused by the imaging system.

Fig. 3

Fast-moving neutrophils in suspension culture captured and visualized with the DHVM: (a) to (c) the same in-focus cell, pointed to by an arrowhead, at a given z distance but at various time points of 400, 1300, and 2300ms , respectively; (d) to (f) the corresponding relative locations of the pointed-to cell in (a) to (c) in a 4-D (including time) mapping. Axis dimensions are as indicated in (d) and the scale bar is 20μm .

014007_1_026801jbo3.jpg

Video 2

Eholovideo (digital holographic video) of fast-moving neutrophil cells. This video is a reconstructed eholovideo from a 3-s eholovideo (not shown) recorded at frame rate of 10fps . The focusing plane z was randomly selected to illustrate the movement of neutrophils in time sequence in 3-D space. One of the in-focus cells, pointed to by an arrowhead, is visible only at some certain time points, e.g., at 400, 1300, and 2300ms , indicating that cells move in and out of the chosen focusing plane over the time. Neutrophil cells were grown in suspension culture at the time of recording. Scale bar is 50μm (QuickTime 2.5MB ). 2

014007_1_026801jbov2.jpg
10.1117/1.2841050.2

Video 3

Four-dimensional schematic of the 3-D mapping of a fast-moving neutrophil along the fourth dimension of time (QuickTime, 0.7MB ). 3

014007_1_026801jbov3.jpg
10.1117/1.2841050.3

3.3.

Spatial Distribution of Cells In Vivo

Subpopulations of cells in a living organism have a defined spatial distribution. A single hologram contains the 3-D information that reveals such a spatial distribution of cells at one time point. A large amount of cells in the tail region of a live tadpole are shown in Fig. 4 . From a single recorded hologram [Fig. 4a], we were able to reconstruct the images and quantify the distribution of cells at the moment that the hologram was taken. Fig. 4b and 4c show the reconstructions of Xenopus tadpole tail at z=178 and 58μm for the cell groups marked in white rectangles, respectively, which were enlarged for the convenience of visualization. The in-focus cell groups are randomly chosen from various focal planes through z and plotted in a 3-D mapping to illustrate their spatial relationship [Fig. 4d] and the distribution in a 3-D tissue space at the time point. The arrowhead in Fig. 4d points to the cells located in the enlarged area marked as “group b” in Fig. 4b at the z=58μm plane; the arrow points at the cells in one of the enlarged areas marked as “group c” in Fig. 4c at the z=178μm plane. The whole reconstruction through the z axis is given in Video 4 , in which one can choose to view the in-focus cells at any arbitrary plane through depths.

Fig. 4

Three-dimensional distribution of cell population in vivo captured and visualized in a live Xenopus tadpole by a single digital hologram: (a) hologram; (b) in-focus cell group (group b) at the z=58μm plane; (c) in-focus cell groups (group c) at z=178μm ; where the in-focus cell images (in white rectangles) were enlarged to be visualized; (d) 3-D mapping of the in-focus cell groups, with the arrowhead pointing at “group b” cells in (b) and the arrow sign pointing at “group c” cells in (c). Axis dimensions are as indicated and scale bar is 32μm .

014007_1_026801jbo4.jpg

Video 4

Eholovideo reconstruction of a tissue block in vivo. An ehologram (digital hologram) was recorded [Fig. 4a in the main text] using Xenopus tadpole tail tissues and reconstructed at various focal planes (31 steps) with a step interval of 5.5μm (QuickTime 1.6MB ). 4

014007_1_026801jbov4.jpg
10.1117/1.2841050.4

3.4.

Digital Holographic Video of Fast Cell Movement in Blood Stream In Vivo

Cell movement in blood streams has been reported with video microcopy for decades. However, when cells move out of the focal plane of the objective lens at a speed such as in the blood stream, conventional video microscopy is not able to capture the movement. With our digital holographic video system, we recorded in vivo blood cells in the tadpole blood stream. Fast-moving cells could be tracked and quantified in three dimensions. A 3-s eholovideo (not shown) of the tadpole tail with blood vessels was recorded, in which two blood stream lines are seen when playing back the video. Figure 5a shows the 14th frame ( t=1400ms , z=0 ) where, the two blood streams were marked as “Main stream” and “Stream 2” with red arrows pointing at the blood flow directions. The eholovideo was then reconstructed at z=57μm against the time axis (Video 5 ). The spatial distribution of blood cells in the “Main stream” is tracked by reconstructing holograms at each time point (t) and scanning through the z axis in a range of 30 to 120μm . Fig. 5c and 5e illustrate two examples of such tracked blood cell, marked with a red dashed circle, at z=103μm and t=200ms and z=75μm and t=1100ms , respectively. The arrows indicate the direction of cell movement. For the purpose of visualization, the “Main stream” blood vessel is highlighted by blurring the surrounding area of the stream. The movement of a cell in three dimensions in the stream was schematically illustrated in four dimensions by plotting the tracked blood cell along time axis with an interval of 100ms (Video 6 ). Figure 5d and 5f show the relative locations [corresponding to Figs. 5b and 5d, respectively] of the tracked cell in the “Main stream” in a 4-D representation, the arrows also indicate the direction of the cell motion.

Fig. 5

Visualization of fast-moving cells in blood vessel in vivo using the DHVM: (a) an eholoimage of tadpole tissue, the first frame (t=0) from a 3-s eholovideo recorded with an ×32 lens; (b) the sum of subtracted holograms (at z=75μm ) of a tadpole blood vessel, tracked two blood streams, i.e., “Main stream” and “Stream 2”; (c) and (e) reconstructed cell (marked in a red circle) in blood stream “Main stream” (highlighted) at time points of 200 and 1100ms , and at a z distances of (c) 103 and (e) 75μm , respectively; (d) and (f) schematic mapping of the tracked cell in a 4-D representation, i.e., 3-D space and time course; the time axis dimension is 3000ms and the space dimensions are as indicated, the scale bar is 20μm ; (g) to (l) reconstruction images of an individual cell extracted from Video 8. The individual cell, pointed to and followed by a red arrow, flows in the blood stream “Stream 2,” and is reconstructed at a z=49μm distance and various time points, as indicated in each image. The frame numbers shown are identical to the frame numbers in Video 5 and the whole frame size is 43 (width)×41 (height) μm .

014007_1_026801jbo5.jpg

Video 5

Visualization of tadpole blood vessels in vivo. A live tadpole was used to record a 3-s eholovideo (not shown) with 10fps frame rate. The eholovideo is reconstructed along the time axis at the reconstruction distance 75μm . We see that only several blood cells are in focus at this reconstruction distance. The full frame size is 273 (width) ×204μm (height) (QuickTime 2.7MB ). 5

014007_1_026801jbov5.jpg
10.1117/1.2841050.5

Video 6

Four-dimensional blood cell movement in the tadpole blood vessel. Upper frame, eholovideo of tadpole blood vessel which was recorded at the frame rate of 10fps over a 3-s period. Individual blood cells are detected to show their migration in the blood vessel along the time axis. The main blood stream, marked as “Main stream” in Fig. 5a in the main text, was highlighted as a clear image against a blurred background area. Lower frame, schematic 3-D mapping of the blood cell movement in the “Main stream” (QuickTime 2.1MB ). 6

014007_1_026801jbov6.jpg
10.1117/1.2841050.6

Video 8

Visualization of a single tadpole blood cell flowing in the blood stream in vivo. An individual blood cell in blood stream [marked as “Stream 2” in Fig. 5a in the main text] was recorded in a 3-s eholovideo (digital holographic video, not shown). The recorded eholovideo was then reconstructed along each time point (interval 100ms ) and along z to locate the best focusing plane. This video tracks the in-focus reconstructed cell over a 2400-ms time. The full frame size is 43 (width)×41 (height) μm (mov, 1.3MB ). 8

014007_1_026801jbov8.jpg
10.1117/1.2841050.8

In the reconstruction processing, however, image quality was degraded due to the severe speckle observed. To eliminate speckle caused by the scattered light from the highly populated cells in the living tissue in Video 5, a motion detection technique, i.e., background subtraction, was utilized as a postreconstruction processing. The postreconstruction processing also includes a deshaking tool (Virtualdub 1.6.3). A deshaking processing is essential when a living tissue is chosen as the recording object because the bulk vibration of tissue was inevitable and it is impossible to obtain a “pure” background image for the processing. The resulting background-subtracted video frames (Video 7 ) were formed by subtracting frame I(k) from I(k+3) , where, k=1,2, , FN3 , and FN is the total number of frames in Video 5. For the purpose of visualization, a sum of all subtracted reconstruction images is calculated to show the blood stream tracks [Fig. 5b]. A moving blood cell in “Stream 2” is further traced from this processing and shown in Video 8 , which lasts over 2400ms along time axis with interval of 100ms . Video 8 shows cell movement at the actual recording rate of 10fps and at a z distance of 49μm . Some of the reconstructed frames are selected [Figs. 5g, 5h, 5i, 5j, 5k, 5l] from Video 8 to show the cell coordinators in four dimensions (x,y,z,t) . Therefore, the moving speed of the cell can be estimated by substituting the sensor pixel properties in x and y . This individual blood cell in “Stream 2” flows with an estimated average speed 34μms1 , which is much slower than the cells in the “Main stream” (estimated average speed 195μms1 ) and the cell size is measured approximately 2.8μm in diameter. An overall image resolution of around 1μm in this system is achieved.

Video 7

Postreconstruction processing of the tadpole blood vessels eholovideo. A live tadpole was used to record a 3-s eholovideo (not shown) with a 10fps frame rate. The eholovideo was reconstructed (Video 5) along the time axis at the reconstruction distance 75μm . Video 5 was then processed by a motion detection technique by a deshaking processing and subtracting the k ’th frame from the (k+3) ’th frame, k=1,2,,FN3 , and FN is the total number of frames in Video 5. The full frame size is 273 (width)×204μm (height) (QuickTime 2.1MB ). 7

014007_1_026801jbov7.jpg
10.1117/1.2841050.7

4.

Conclusions

In this paper, in vitro and in vivo fast-moving cells in a 3-D space were filmed nonintrusively using our DHVM system. The cell locations and 3-D spatial distributions of cells at one time point were recreated in a computer by scanning through the recorded volume at various depths for the chosen hologram frame in the recorded eholovideo; while cell motions were recreated by reconstructing a sequence of eholovideo frames at a specific distance of cell focal planes in a time coordinate. In comparison with a conventional microscopic imaging system, an ehologram records a larger area and in a 3-D volume.

Recording fast-moving cells in suspension and in a blood stream in vivo has not been possible with currently available imaging systems because cells move in and out of the microscopic focal plane rapidly in a 3-D volume. With our digital video microscopic system, such rapid movements were recorded and recreated in a volumetric space. The DHVM system we report here is capable of recording video holographic images of cells on a glass slide, in a tissue culture dish, and in blood vessels. Both intensity and phase-contrast digital microscopic holograms of a 3-D object volume can be recorded. The recording of a nonstained live cell distribution in vivo and the movement of cells in a blood stream in vivo offers the exciting possibility of tracking fast cell movement in three dimensions. In comparison with other holographic systems, simplicity in video recording geometry is one of the main advantages of a DHVM system. Computer simulation of diffraction wave and video processing techniques are conveniently applied in the reconstruction of specimen eholovideos. The reference beam in reconstruction can be omitted if a parallel beam is used to illuminate the recording specimen. With selection of objective lenses, high spatial and temporal resolution can be achieved. The disadvantages of a DHVM include a very limited recording volume (several cubic centimeters) and restricted longitudinal image resolution when a highly magnifying lens is selected.

A limitation of our DHVM is the requirement of partial transparency of the specimen and the limited resolution for submicrometer recordings, where the scattering effect must be taken into account when carrying out size and distance measurements. When scattering becomes significant, the recorded hologram will be significantly degraded, to the point that no reference beam can reach the recording plane, in which case, no hologram will be formed and therefore no object image can be reconstructed by any of the mentioned reconstruction algorithms. This problem can be solved if an off-axis recording scheme, i.e., an angularly separated reference beam, is used. However, when cell population density is rather high, too much scattered light will produce severe speckle in a hologram due to the scattered coherence laser light from dense cells. Therefore, further developments on reconstruction algorithms and postreconstruction algorithms, recording geometry, and light sources are necessary.

Acknowledgments

We thank the funding bodies MRC (British Medical Research Council to HS, JW and MZ) and the Welcome Trust (to MZ). The work could not be started without the outcomes from a NERC (British National Environmental Research Council to JW and MAP) supported project (Biomediation of Sediment Erosion using laser holography). Thanks are also given to Drs. Stefan Hoppler and Fei Liu for their help with the tadpoles.

References

1. 

E. Cukierman, R. Pankov, D. R. Stevens, and K. M. Yamada, “Taking cell-matrix adhesions to the third dimension,” Science, 294 1708 –1702 (2001). https://doi.org/10.1126/science.1064829 0036-8075 Google Scholar

2. 

B. Knight, C. Laukaitis, N. Akhtar, N. A. Hotchin, M. Edlund, and A. R. Horwitz, “Visualizing muscle cell migration in situ,” Curr. Biol., 10 576 –585 (2000). 0960-9822 Google Scholar

3. 

E. Sahai and C. J. Marshall, “Differing modes of tumour cell invasion have distinct requirements for Rho/ROCK signalling and extracellular proteolysis,” Nat. Cell Biol., 5 711 –719 (2003). 1465-7392 Google Scholar

4. 

D. J. Webb and A. F. Horwitz, “New dimensions in cell migration,” Nat. Cell Biol., 5 690 –692 (2003). 1465-7392 Google Scholar

5. 

M. Zhao, B. Song, J. Pu, J. V. Forrester, and C. D. McCaig, “Direct visualization of a stratified epithelium reveals that wounds heal by unified sliding of cell sheets,” FASEB J., 17 397 –406 (2003). 0892-6638 Google Scholar

6. 

M. Zhao, B. Song, J. Pu, T. Wada, B. Reid, G. Tai, F. Wang, A. Guo, P. Walczysko, Y. Gu, T. Sasaki, A. Suzuki, J. V. forrester, H. R. Bourne, P. N. devreotes, C. D. McCaig, and J. M. Penninger, “Electrical signals control wound healing through phosphatidylinositol-3-OH kinase-γ and PTEN,” Nature (London), 442 457 –460 (2006). https://doi.org/10.1038/nature04925 0028-0836 Google Scholar

7. 

H. Xu, A. Manivannan, D. Garry, J. Liversidge, P. F. Sharp, J. V. Forrester, and I. J. Crane, “Evaluation of leukocyte dynamic in mouse retinal circulation with scanning laser ophthalmoloscopy (Video report),” Br. J. Ophthamol., 85 765 –890 (2001). 0007-1161 Google Scholar

8. 

E. Malkiel, J. Sheng, J. Katz, and J. R. Strickler, “The three dimensional flow field generated by a feeding calanoid copepod measured using digital holography,” J. Exp. Biol., 206 3657 –3666 (2003). https://doi.org/10.1242/jeb.00586 0022-0949 Google Scholar

9. 

G. Pan and H. Meng, “Digital holography of particle fields: reconstruction by use of complex amplitude,” Appl. Opt., 42 827 –833 (2003). https://doi.org/10.1364/AO.42.000827 0003-6935 Google Scholar

10. 

U. Schnars and W. O. Jüptner, “Digital recording and numerical reconstruction of holograms,” Meas. Sci. Technol., 13 R85 –101 (2002). https://doi.org/10.1088/0957-0233/13/9/201 0957-0233 Google Scholar

11. 

J. H. Milgram and W. Li, “Computational reconstruction of images from holograms,” Appl. Opt., 41 853 –864 (2002). https://doi.org/10.1364/AO.41.000853 0003-6935 Google Scholar

12. 

H. Sun, D. C. Hendry, M. A. Player, and J. Watson, “In situ underwater electronic holographic camera for studies of plankton,” IEEE J. Ocean. Eng., 32 (2), 1 –10 (2007). 0364-9059 Google Scholar

13. 

H. Sun, H. Dong, M. A. Player, J. Watson, D. M. Paterson, and R. G. Perkins, “In-line digital video holography for the study of erosion processes in sediments,” Meas. Tech., 13 L7 –12 (2002). 0543-1972 Google Scholar

14. 

H. Sun, M. A. Player, J. Watson, D. C. Hendry, R. G. Perkins, G. Gust, and D. M. Paterson, “The use of digital/electronic holography for biological applications,” J. Opt. A, Pure Appl. Opt., 7 S399 –S407 (2005). https://doi.org/10.1088/1464-4258/7/6/022 1464-4258 Google Scholar

15. 

W. B. Xu, M. H. Jericho, I. A. Meinertzhagenm, and H. J. Kreuzer, “Digital in-line holography for biological applications,” Proc. Natl. Acad. Sci. U.S.A., 98 11301 –11305 (2001). https://doi.org/10.1073/pnas.191361398 0027-8424 Google Scholar

16. 

L. Xu, X. Peng, J. Miao, and A. K. Asundi, “Studies of digital microscopic holography with applications to microstructure testing,” Appl. Opt., 40 5046 –5051 (2001). https://doi.org/10.1364/AO.40.005046 0003-6935 Google Scholar

17. 

W. Xu, M. H. Jericho, H. J. Kreuzer, and I. A. Meinertzhagen, “Tracking particles in four dimensions with in-line holographic microscopy,” Opt. Lett., 28 164 –166 (2003). https://doi.org/10.1038/nature01433 0146-9592 Google Scholar

18. 

P. Marquet, B. Rappaz, and P. Magistretti, “Digital holographic microscopy: a non-invasive contrast imaging technique allowing quantitative visualization of living cells with subwavelength axial accuracy,” Opt. Lett., 30 468 –470 (2005). https://doi.org/10.1364/OL.30.000468 0146-9592 Google Scholar

19. 

F. Dubois, N. Callens, C. Yourassowsdy, M. Hoyos, P. Kurowski, and O. Monnom, “Digital holographic microscopy with reduced spatial coherence for three-dimensional particle flow analysis,” Appl. Opt., 45 846 –870 (2006). 0003-6935 Google Scholar

20. 

C. J. Mann, L. Yu, and M. K. Kim, “Movies of cellular and sub-cellular motion by digital holographic microscopy,” Biomed. Eng. Online, 5 (21), (2006). 1475-925X Google Scholar

21. 

J. Sheng, E. Malkiel, and J. Katz, “Digital holographic microscope for measuring three-dimensional particle distributions and motions,” Appl. Opt., 45 3893 –3901 (2006). https://doi.org/10.1364/AO.45.003893 0003-6935 Google Scholar

22. 

U. Schnars and W. O. Jüptner, “Direct recording of holograms by a CCD target and numerical reconstruction,” Appl. Opt., 33 2187 –2197 (1994). 0003-6935 Google Scholar

23. 

R. B. Owen and A. A. Zozulya, “In-line digital holographic sensor for monitoring and characterizing marine particulates,” Opt. Eng., 39 2187 –2197 (2000). https://doi.org/10.1117/1.1305542 0091-3286 Google Scholar

24. 

U. Schnars and W. O. Jüptner, “Digital recording and numerical reconstruction of holograms,” Meas. Sci. Technol., 13 R85 –101 (2002). https://doi.org/10.1088/0957-0233/13/9/201 0957-0233 Google Scholar

25. 

X. M. Huang, J. M. Zuo, and J. H. Spence, “Wavefront reconstruction for in-line holograms formed by pure amplitude objects,” Appl. Surf. Sci., 148 229 –234 (1999). https://doi.org/10.1016/S0169-4332(99)00079-3 0169-4332 Google Scholar

26. 

D. Lebrun, S. Belaid, and C. Ozkul, “Hologram reconstruction by use of optical wavelet transform,” Appl. Opt., 38 3730 –374 (1999). 0003-6935 Google Scholar

27. 

M. Liebling, T. Blu, and M. Unser, “New multiresolution wavelet bases for digital holography,” IEEE Trans. Image Process., 12 29 –43 (2003). https://doi.org/10.1109/TIP.2002.806243 1057-7149 Google Scholar

28. 

D. C. Champney, Fourier Transforms and Their Physical Applications, 138 –180 Academic Press, London (1973). Google Scholar

29. 

C. Buraga-Lefebvre, S. Coëtmellec, D. Lebrun, and C. Özkul, “Application of wavelet transform to hologram analysis: three-dimensional location of particles,” Opt. Lasers Eng., 33 409 –421 (2000). https://doi.org/10.1016/S0143-8166(00)00050-6 0143-8166 Google Scholar

30. 

L. Onural and M. Kocatepe, “Family of scaling chirp functions, diffraction, and holography,” IEEE Trans. Signal Process., 43 1568 –1578 (1995). https://doi.org/10.1109/78.398718 1053-587X Google Scholar

31. 

H. Dong, C. Khong, M. A. Player, M. Solan, and J. Watson, “Algorithms and applications for electronically-recorded holography,” 354 –365 (2004). Google Scholar
©(2008) Society of Photo-Optical Instrumentation Engineers (SPIE)
Hongyue Sun, Bing Song, Hongpai Dong, Brian Reid, Michael A. Player, John Watson, and Min Zhao "Visualization of fast-moving cells in vivo using digital holographic video microscopy," Journal of Biomedical Optics 13(1), 014007 (1 January 2008). https://doi.org/10.1117/1.2841050
Published: 1 January 2008
Lens.org Logo
CITATIONS
Cited by 30 scholarly publications and 9 patents.
Advertisement
Advertisement
RIGHTS & PERMISSIONS
Get copyright permission  Get copyright permission on Copyright Marketplace
KEYWORDS
Video

Holograms

Digital holography

Holography

Blood

In vivo imaging

Visualization

Back to Top