We describe the use of a standard optical microscope to perform quantitative measurements of cellular mass, volume, and density through a combination of bright field and differential interference contrast imagery.
We describe the use of a standard optical microscope to perform quantitative measurements of mass, volume, and density on cellular specimens through a combination of bright field and differential interference contrast imagery. Two primary approaches are presented: noninterferometric quantitative phase microscopy (NIQPM), to perform measurements of total cell mass and subcellular density distribution, and Hilbert transform differential interference contrast microscopy (HTDIC) to determine volume. NIQPM is based on a simplified model of wave propagation, termed the paraxial approximation, with three underlying assumptions: low numerical aperture (NA) illumination, weak scattering, and weak absorption of light by the specimen. Fortunately, unstained cellular specimens satisfy these assumptions and low NA illumination is easily achieved on commercial microscopes. HTDIC is used to obtain volumetric information from through-focus DIC imagery under high NA illumination conditions. High NA illumination enables enhanced sectioning of the specimen along the optical axis. Hilbert transform processing on the DIC image stacks greatly enhances edge detection algorithms for localization of the specimen borders in three dimensions by separating the gray values of the specimen intensity from those of the background. The primary advantages of NIQPM and HTDIC lay in their technological accessibility using “off-the-shelf” microscopes. There are two basic limitations of these methods: slow z-stack acquisition time on commercial scopes currently abrogates the investigation of phenomena faster than 1 frame/minute, and secondly, diffraction effects restrict the utility of NIQPM and HTDIC to objects from 0.2 up to 10 (NIQPM) and 20 (HTDIC) μm in diameter, respectively. Hence, the specimen and its associated time dynamics of interest must meet certain size and temporal constraints to enable the use of these methods. Excitingly, most fixed cellular specimens are readily investigated with these methods.
The use of optical microscopy is now ubiquitous in the investigation of cellular organisms. Due to their low endogenous absorbance and weak scattering properties over the visible optical spectrum, cells do not strongly affect the amplitude of optical waves traversing them and thus appear semitransparent when imaged with standard bright field microscopes. Cellular specimens do, however, slow down optical waves traveling through them in a manner that is linearly related to the amount of local mass density in a particular region of space through which the light travels. Utilization of this heterogeneous time lag or “phase” profile of optical waves transmitted through microscopic specimens was first described in 1935 by Frits Zernike1 and experimentally realized by Zernikein 19422. Zernike was awarded the Nobel prize in 1953 for this achievement. Zeiss commercialized this modality in 19453. In 1955, Smith and Nomarski would present their initial work on the use4 and theory5 of differential interference contrast (DIC) microscopy, a modality that uses the spatial gradient of phase as a contrast mechanism. DIC was commercialized by Zeiss in 1965in close collaboration with Nomarski6. In 1981, two laboratories demonstrated the first recorded live cell DIC imagery with the incorporation of video cameras into the optics train of the DIC microscope7,8. The era of live cell imaging was born.
Since this time, the execution of both phase contrast and DIC on commercial microscopes has largely been unchanged. These methods are principally utilized by biologists to produce images of cells for qualitative purposes: monitoring of morphology, tracking of subcellular structures, and investigations of membrane dynamics9. These techniques are qualitative in their “off-the-shelf” configuration as both phase and DIC images are arbitrary functions of the light source intensity, the illumination optics settings, and CCD camera gain, gamma, and exposure settings.
A small legion of physicists and optical engineers has endeavored to make commercial imaging modalities quantitative. Among the first efforts were two letters to Nature in 1952 and 1953 in which the physician-turned-biophysicist Robert Barer demonstrated the use of phase microscopy to determine cellular dry mass of cells by estimating phase shifts through these cell types using a commercially available phase microscope10,11. The field has developed a multitude of techniques over the ensuing years based around three basic label-free contrast mechanisms: phase microscopy10,11, DIC microscopy12-17, and bright field18-22 to determine optical path-length, phase, mass density, refractive index, and cellular volume.
In parallel, a large collection of custom optical instruments have also been developed since the 1950s, and have made far-reaching optical measurements ranging from applications in parasite growth23, to documenting the cell cycle24 to investigating membrane dynamics of red blood cells25. In particular, the past ten years has seen a wealth of label-free quantitative microscopy in the form of diffraction phase microscopy26, tomographic phase microscopy27, digital holographic microscopy28, phase sensitive optical coherence microscopy29, spatial light interference microscopy30, Hilbert phase microscopy31, and quantitative phase microscopy32. Despite their collective successes, these instruments have not been disseminated to the larger field of biological researchers owing, mostly, to their complex instrumentation and computational requirements.
Herein we describe the use of a standard optical microscope to perform quantitative measurements of mass, volume, and density on cellular specimens through a combination of bright field and DIC imagery. Two primary approaches are presented: noninterferometric quantitative phase microscopy (NIQPM), to perform measurements of total cell mass and subcellular density distribution, and Hilbert transform differential interference contrast microscopy (HTDIC), to determine volume. The primary advantages of NIQPM and HTDIC lay in their technological accessibility. The imaging conditions required for their successful execution are within the scope of normal operation of most commercially available microscopes. Additionally, the post-processing algorithms are stable, quick, and robust – having been implemented in MATLAB using fast Fourier transform (FFT) based algorithms whenever possible.
NIQPM is a method to reconstruct phase and the axially integrated mass density of cellular specimens from bright field imagery. Summation of this axially integrated mass density over the area of the specimen gives the total dry mass content of the specimen. The NIQPM protocol is based on the experimental foundations laid by Paganin and Nugent18,19 – in which it was demonstrated that the phase profile of a cell could be reconstructed from through-focus bright field imagery of the sample – and the theoretical work of Frank, Altmeyer, and Wernicke20 – on solving the paraxial wave models in an efficient FFT-based manner. The connection of phase to the dry mass density is based upon the work by Barer10,11 and Popescu33.
Volumetric information can be obtained from through-focus DIC imagery under high NA illumination conditions that enable optical sectioning of the specimen along the optical axis. Hilbert transform processing on the DIC image stacks greatly enhances edge detection algorithms for localization of the specimen borders in three dimensions by separating the gray values of the specimen intensity from those of the background. This work originates with Arinson et al.34 although we have introduced both Fourier filtering methods to enhance contrast and a Sobel-based edge detection method for automated volumetric analysis of the sample. We have also validated HTDIC previously on polystyrene spheres ranging in size from the diffraction limit up to 20 μm in diameter36.
While both NIQPM and HTDIC are technologically accessible owing to their development on commercial microscopes, the methods are fundamentally limited by the hardware configuration of the microscopes themselves. The primary limitations of these techniques are two-fold: slow z-stack acquisition time on commercial scopes, due to translation of the entire sample stage as opposed to just the objective lens, currently limits the investigation of phenomena faster than roughly 1 frame/minute, and secondly, diffraction effects restrict the utility of NIQPM and HTDIC to objects ranging in size from 0.2 up to 10 and 20 μm in diameter, respectively. Hence, the specimen and its associated time dynamics of interest must meet certain size and temporal constraints to enable the use of these methods on typical “off-the-shelf” instruments. Excitingly, most fixed cellular specimens are readily investigated with these methods.
An overview of the NIQPM and HTDIC protocols are given in Figure 1. In Figure 2 we illustrate optimal and suboptimal through-focus imaging under both low and high NA illumination conditions for both bright field and DIC imagery. Figures 3 and 4 demonstrate the parameter dependence of the NIQPM algorithm highlighting successful and unsuccessful implementations. Figure 5 illustrates the steps involved in the HTDIC image-processing algorithm and demonstrates the optimal implementation of the algorithm to determine cellular volume.
1. Microscope Specifications
To carry out imaging in the correct fashion the microscope should have the following specifications:
2. Differential Interference Contrast (DIC) Z-stack Acquisition
3. Bright Field (BF) Z-stack Acquisition
4. Exporting Z-stack Images
5. Volume Measurements
6. Mass Measurements
Correct sample illumination during through-focus image acquisition is critical to the successful implementation of the NIQPM and HTDIC algorithms. In Figure 2 we illustrate low and high NA illumination under both DIC and bright field contrast for a polystyrene sphere and the human colorectal adenocarcinoma cell line SW620. Figures 2A, 2C, 2I, and 2K demonstrate optimal imaging for NIQPM. Figures 2F, 2H, 2N, and 2P demonstrate optimal imaging for HTDIC.
Figures 3 and 4 demonstrate the parameter dependence of the NIQPM algorithm highlighting both successful and unsuccessful implementations. In Figure 3 we explore the phase profile of a 4.8 μm diameter polystyrene sphere – the expected profile is known theoretically and can thus be compared directly to the NIQPM reconstruction. Figures 2I-L present the best reconstruction obtainable for the sphere. Diffraction effects at both the boundaries of the sphere and inside prevent stable reconstruction at all points on the sphere. The central region of the sphere can be captured with NIQPM with a % error of 1-5%, Figure 2L.
Cellular specimens, whose phase properties are not known a priori, can be reconstructed using NIQPM in conjunction with a procedure to compare a “pseudo DIC” image – computed from the reconstructed phase – to an actual DIC image, determined from a direct optical measurement, of the cell. NIQPM has one free parameter – the plane in the bright field image stack at which the calculation is centered. This central focal plane should be adjusted until the pseudo DIC and true DIC images look as similar as possible. Figures 4E-G demonstrate an out-of-focus pseudo DIC image, an optimal pseudo DIC image, and the corresponding DIC image of the cell taken with an illumination NA of 0.9. Interestingly, the best phase map and corresponding pseudo DIC image do not necessarily correspond to an in focus bright field image, viz. Figures 1A and 1B.
Lastly, Figure 5 illustrates the steps involved in the HTDIC image-processing algorithm. From the DIC through-focus imagery under NA = 0.9 illumination, Figures 5A and 5D, the Hilbert transform is preformed to remove the bas relief of the DIC images, Figures 5B and 5E. This comes with some blurring along the optical axis that can be removed from high pass Fourier filtering, Figures 5C and 5F. These final images are easily segmented to determine the area in each cross sectional plane of the specimen to infer the total cellular volume.
Figure 1. NIQPM and HTDIC Workflow. (1.) Microscopic specimens such as cells should be mounted on microscope slides with cover glass affixed over the sample using Fluoromount G. (2.) Through-focus images acquired under both DIC and bright field contrast with a standard “off-the-shelf” microscope form the input to the image processing algorithms. (3.) Post-processing of images in MATLAB to determine cell volume from DIC and cellular mass distribution from bright field imagery. (4.) Quantitative endpoint metrics: heat maps and bar graphs. Please click here to view a larger version of this figure.
Figure 2. Through-focus DIC and bright field imagery of polystyrene spheres and SW620 cell lines. The optimal imaging conditions for NIQPM are bright field contrast with 0.1 NA illumination. For HTDIC DIC, 0.9 NA illumination is optimal. (A, B) En face bright field imagery of 4.8 μm diameter polystyrene sphere subject to 0.1 NA and 0.9 NA illumination, respectively. (C, D) Cross sectional bright field imagery of 4.8 μm diameter polystyrene sphere subject to 0.1 NA and 0.9 NA illumination, respectively. (E, F) En face DIC imagery of 4.8 μm diameter polystyrene sphere subject to 0.1 NA and 0.9 NA illumination, respectively. (G, H) Cross sectional DIC imagery of 4.8 μm diameter polystyrene sphere subject to 0.1 NA and 0.9 NA illumination, respectively. (I, J) En face bright field imagery of SW620 colorectal cancer cell line subject to 0.1 and 0.9 NA illumination, respectively. (K, L) Cross sectional bright field imagery of SW620 cell subject to 0.1 and 0.9 NA illumination, respectively. (M, N) En face DIC imagery of SW620 cell subject to 0.1 NA and 0.9 NA illumination, respectively. (O, P) Cross sectional DIC imagery of SW620 cell subject to 0.1 NA and 0.9 NA illumination. Please click here to view a larger version of this figure.
Figure 3. Optimal and suboptimal phase reconstructions:4.8 μm diameter polystyrene sphere. Each row investigates the focal plane dependence of the en face bright field intensity (first column) and the reconstructed phase profile as the focal plane begins at the midpoint through the sphere, z = 0, and moves past this plane in 1 μm steps. (A, E, I, M) En face bright field intensity, (B, F, J, N) corresponding phase reconstruction, (C, G, K, O) comparison of the reconstructed phase profile along the diagonal of the phase map, blue circles, theoretical phase profile, in red, and manufacturer’s error tolerance, black lines, and (D, H, L, P) percent (%) error of the reconstructed phase with respect to the theoretical phase profile. Please click here to view a larger version of this figure.
Figure 4. Optimal and suboptimal phase reconstructions: SW620 cell. The two columns compare suboptimal and optimal bright field image planes (A, B) from the z-stack to center the phase calculation. The corresponding phase reconstructions (C, D) are used to generate a pseudo DIC image (E, F). These are compared to the true DIC image, NA = 0.9 illumination, in (G). The pseudo DIC image best matching the DIC image determines the correct input images from the bright field image stack to use in phase reconstruction. Finally, phase is mapped to the projected mass density (H, I). Please click here to view a larger version of this figure.
Figure 5. Optimal DIC contrast enhancement using HTDIC: polystyrene spheres and SW620 cell. (A) Cross-sectional DIC imagery of a 4.8 polystyrene sphere, (B) corresponding HTDIC image, (C) Fourier filtered HTDIC image. The axial dimension of the sphere images have been scaled to account for refractive index mismatch of the specimen (1.597) and the mounting media (1.4). (D-F) demonstrate the same image types for an SW620 cell, however, with no refractive index mismatch correction owing to the weak index contrast of the specimen. No thresholding has been performed on these images. Please click here to view a larger version of this figure.
In general, NIQPM is a diffraction-limited technique validated on optical path-lengths ranging from 0.25-44.7. Validation was performed on n = 1.596 polystyrene spheres ranging in diameter from 0.11-9.8 μm suspended in Fluoromount G (data not shown). Cells possess optical path-lengths that range from nearly 0-7.
When measuring the density distribution of a specimen one may find that the pseudo DIC image looks fine while the density map possesses unwanted background contributions. This is due to noise in the bright field imagery of the sample used as the input to the NIQPM algorithm. Two possible modifications to the NIQPM method can be used to eliminate the unwanted background contributions: firstly, one may use longer exposure times to enhance signal-to-noise in the bright field imagery. This is especially critical when imaging optically thin specimens. Alternatively, one may reduce the axial spacing of the z-stack acquisition from 0.1-0.05 μm and average 2-3 planes together before inputting them into the NIQPM algorithm.
HTDIC based volume determination has been validated on polystyrene spheres ranging in diameter from 1-20 μm and has been cross-validated with confocal fluorescence microscopy35. The technique overestimates the volume of objects below the diffraction limit of the system due to diffraction effects associated with the point spread function of the optics.
Troubleshooting the HTDIC volume determination maybe required for small or thin specimens. The primary source of error comes from image segmentation using the built-in MATLAB Sobel-based edge detection used to determine the borders of the cell in the cross sectional images. The value of “threshold” in section 4 of the HTDIC program is critical in obtaining the best segmentation. Results tend to vary in a nonlinear way – with small changes in “threshold” drastically changing the enclosed area resulting from the edge detection. The current program presents the user with options to perform image segmentation from either the DIC z-stack, the HTDIC z-stack, or the F-HTDIC z-stack.
DIC alone is often better for large specimens while HTDIC and F-HTDIC work for smaller specimens (<5 μm in diameter).
A mask can enhance the ability of the method to segment the image. Section 2 allows the user to determine a mask to apply to the image cube. We recommend the use of a rectangular mask that truncates the horizontal (x-extent) of the cube while preserving the vertical (y-extent) of the cube.
In summary, NIQPM and HTDIC are technologically accessible quantitative imaging modalities that can be performed on standard “off-the-shelf” optical microscopes, in contrast to most existing methods. Central to these techniques is through-focus imagery of the sample under the appropriate imaging conditions: low NA illumination for NIQPM, and high NA illumination for HTDIC. The methods presented here can be generalized for use on systems other than those demonstrated here. The primary criteria required for these procedures is a microscope for which the user has control over the z-movement of the stage and the ability to acquire images as the focal position of the objective lens is varied. The methods presented are suitable for imaging fixed cells or slowly moving biological specimens.
The authors have nothing to disclose.
This work was supported by grants from the National Institutes of Health (U54CA143906 to K.G.P., O.J.T.M and R01HL101972 to O.J.T.M.) and an Oregon Medical Research Foundation Early Clinical Investigator Award (K.G.P.). O.J.T.M. is an American Heart Association Established Investigator (13EIA12630000). We thank Dr. Eric Anderson of the Knight Cancer Institute for preparing cell samples used in this work.
Zeiss Axio Imager 2 microscope | Carl Zeiss MicroImaging GmbH, Germany | Axio Imager D2 | Microscope |
Green filter (λ = 540 ± 25 nm) | Chroma Technology Corp., Bellows Falls, Vermont | D540/25x | Green filter |
SlideBook 5.5 software | Intelligent Imaging Innovations, Denver, Colorado | Image acquistion software | |
Polystyrene microspheres | Bangs Laboratory, Inc., Fishers, IN | PS06N | Polystyrene spheres |
Fluoromount-G | SouthernBiotech, Birmingham, Alabama | 0100-01 | Mounting media |