Quantitative Optical Microscopy: Measurement of Cellular Biophysical Features with a Standard Optical Microscope

Summary

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.

Abstract

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.

Introduction

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 Zernike¹ and experimentally realized by Zernikein 1942². Zernike was awarded the Nobel prize in 1953 for this achievement. Zeiss commercialized this modality in 1945³. In 1955, Smith and Nomarski would present their initial work on the use⁴ and theory⁵ 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 Nomarski⁶. 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 microscope^7,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 dynamics⁹. 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 microscope^10,11. The field has developed a multitude of techniques over the ensuing years based around three basic label-free contrast mechanisms: phase microscopy^10,11, DIC microscopy^12-17, and bright field^18-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 growth²³, to documenting the cell cycle²⁴ to investigating membrane dynamics of red blood cells²⁵. In particular, the past ten years has seen a wealth of label-free quantitative microscopy in the form of diffraction phase microscopy²⁶, tomographic phase microscopy²⁷, digital holographic microscopy²⁸, phase sensitive optical coherence microscopy²⁹, spatial light interference microscopy³⁰, Hilbert phase microscopy³¹, and quantitative phase microscopy³². 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 Nugent^18,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 Wernicke²⁰ – 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 Barer^10,11 and Popescu³³.

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.³⁴ 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 diameter³⁶.

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.

Protocol

1. Microscope Specifications

To carry out imaging in the correct fashion the microscope should have the following specifications:

1. Possess both differential interference contrast (DIC) and brightfield (BF) contrast.
2. Have computer-controlled z-axis movement.
3. Have an adjustable aperture stop to vary the condenser lens numerical aperture. An aperture with a graded rule or electronic readout is required to know the value of the numerical aperture. The numerical aperture should range from 0.1 for NIQPM up to 0.9 (or higher) for HTDIC.
4. The microscope should have a narrow band color filter to be utilized for bright field imaging. This filter is required to fix the refractive increment, a wavelength-dependent facture used to convert phase to mass density for NIQP.

2. Differential Interference Contrast (DIC) Z-stack Acquisition

1. Open the SlideBook software and create a new slide for image collection.
2. Next, open the Focus Window. Under the Filter Set section, select DIC. On the Scope tab, under the Condenser section, adjust the Aperture slide bar to the furthest right position (open all the way). This provides high numerical aperture illumination and enhances optical sectioning of the specimen.
3. Focus on the sample using a DIC objective. Adjust the halogen lamp intensity until the sample is easily visible. To ensure the camera is not being saturated, select the Camera tab in the Focus Window and check that the histogram of pixel intensities is within the dynamic range of the camera.
4. Open the Image Capture window. In the Capture Type section, check the 3D box. In the 3D Capture section, select the Range around center and Return to current location boxes and specify the Z-direction range, number of planes, and step size. The Z-direction range should incorporate the entire sample.
5. In the Filter Set section of the Image Capture window, check the DIC box and specify exposure time. In the Image Information section, name the image (optional), and select Start to initiate the imaging.

3. Bright Field (BF) Z-stack Acquisition

1. Once the microscope has finished collecting the DIC Z-stack of the sample, open the Focus Window and select Open under the Filter Set section. Adjust the Aperture slide bar to the furthest left position (closed all the way) to provide low NA illumination.
2. Change the microscope light-path filter to green. Adjust the halogen lamp intensity until the sample is visible. Ensure there is no saturation of the camera by selecting the Camera tab in the Focus Window and checking that the histogram of pixel intensities is within the specified range.
3. Open the Image Capture window. The 3D Capture settings from the DIC Z-stack acquisition will be displayed.
4. In the Filter Set section of the Image Capture window, check the OPEN box and specify exposure time. In the Image Information section, name the image (optional), and select Start to initiate Z-stack image acquisition.

4. Exporting Z-stack Images

1. Open a Z-stack capture. Change the view to 100%. Adjust the histogram of the pixel values to be displayed between 0 and the maximum pixel value of the camera (4,095 for 12 bit cameras, 65,535 for 16 bit.).
2. Select View > Export > TIFF Series. This will export the Z-stack as a series of TIFF images (one for each Z plane). Save the TIFF series in a separate folder with a name that has an underscore at the end (i.e. "DIC Stack 1_"). Repeat this with each DIC or BF Z-stack.

5. Volume Measurements

1. Open the HTDIC MATLAB program entitled "JoVE_HTDIC_v1.m".
2. Under Section 0 of the HTDIC program, update the dependencies directory variable. Copy and paste the directory containing the hilbert_transform_dic.m and sobel_edge_detect.m files from explorer (PC) in between the single quotes following "dependencies_directory =". ExecuteSection 0 of the JoVE_HTDIC_v1.m program.
3. In Section 1, update the "images_directory". Again, copy and paste the directory containing the through-focus images in .tif form in between the single quotes. ONLY RUN SECTION 1 ONCE.
4. Alignment and rotation of images for Hilbert transform Run Section 2 of the code. A dialog box titled "Define HTDIC Parameters" will appear. Five numbers are required from the user: the focal plane number where the DIC image of the sample is in focus, lateral resolution (μm/pixel in the image), the axial resolution (this is 0.1 for image acquisition in 0.1 μm axial steps), the rotation angle of the DIC image needed to perform the Hilbert transform, typical values are 45 and -135. Lastly, enter the region of interest size, this parameter defines the length of a side of a box that will subsequently appear - typical value is 400. Click Ok.
5. An image of the DIC focal plane specified by the focal plane number will appear with a blue box. Position the box over the feature of interest, e.g. a cell. The box need not be square. Once the box has been positioned over the desired region, double click inside the box.
6. Another figure will now appear: this figure contains a cropped and rotated image of the region of interest selected in the previous step. The contrast of the image should be such that dark features appear on the left while bright features appear on the right. Drag the blue box over the region of interest, reshape as necessary.
7. Section 3 generates a mask for the z-stack cube. There are two types of masks available: a rectangular mask to create a rectangle around the cell and a free-hand tool to outline the cell of interest by hand.
8. To generate a rectangular mask, uncomment line 167 and comment line 170 (comment a line out by placing a "%" at the beginning of that line). Execute Section 3 of the program. Click in the image and drag the mouse to begin defining the rectangular mask. Double click on the box to accept it.
9. To generate a hand-drawn mask, comment line 167 and uncomment line 170 and execute Section 3. Click and draw the desired mask with the mouse. Double click on the mask to accept it.
10. Run Section 4to construct the Hilbert transformed DIC image stack and the corresponding DIC image stack of the area of interest. The mask constructed in Section 3 will be applied to the Hilbert transform DIC stack and the regular DIC stack when "maskON" is set to 1 at line 179. Setting "maskON" to 0 will construct image stacks without applying the mask.
11. Run Section 5 to optimize the image segmentation of the xz cross sectional images of the region of interest. Figure 500, produced by the program, appears displaying three different types of contrast. The success of the algorithm to find the borders of the cell is reliant on a combination of the mask used and the value of "threshold" at line 229 of the program. Begin with a value of 0.5. Adjust the value of threshold and rerun this section of the program until the proper outlining is achieved in one of the columns.
12. If the outlining was best in column 1, use Section 6 to determine the volume using DIC image segmentation. If column 2 gave optimal results, run Section 7 to determine cell volume from Hilbert-transformed DIC imagery. If column 3 gave the optimal results, run Section 8 to determine volume using Fourier-filtered Hilbert-transformed DIC imagery.
13. The measured volume of the sample, reported in cubic μm (fl) is presented in the title of Figure 600, 700, or 800 produced by the program depending on which volume measurement is chosen.

6. Mass Measurements

1. Open the NIQPM MATLAB program entitled "JoVE_NIQPM_v1.m".
2. Under Section 0 of the NIQPM program, update the location of three directories needed to run the program. These are the "dependencies_directory," the "brightfield_directory," and the "dic_directory."
3. Next, run Section 1. This section of code generates a bright field image cube of the full field of the entire set of images. ONLY RUN SECTION 1 ONCE.
4. Next, run Section 2. A dialog box titled "Define NIQPM Parameters" will appear. Four numbers are required from the user: the focal plane number where the bright field image of the sample is in focus, lateral resolution (μm/pixel in the image), the axial resolution (this is 0.1 for image acquisition in 0.1 μm axial steps), and the region of interest size, this parameter defines the length of a side of a box that will subsequently appear - a typical value is 200. Click Ok.
5. An image of the bright field focal plane specified by the focal plane number will appear with a blue box.
   1. If the image is not in focus, double click in the blue box and rerun Section 2, be sure to adjust the focal plane number in the dialog box.
   2. If the image is in focus, drag the blue box around the image and resize it by selecting the nodes of the box and dragging them as necessary. Position the box around the feature of interest, e.g. a cell. The box need not be square. Double click inside the box to accept it.
6. Next, run Section 3 to construct a stack of bright field images cropped to the region of interest.
7. To generate the phase map, pseudo DIC image, and comparisons to bright field imagery and the true DIC image, run Section 4.
8. With the pseudo DIC and true DIC images as similar as possible the mass density map, total mass, and histogram of the cell density can be determined by running Section 5a or 5b. Section 5a carries out automated boarder detection of the field - optimize "threshold" at line 300 - typical values vary between 0.1-1. Rerun this section as needed to optimize the value of threshold. Section 5b carries out mass determination, etc. after the user outlines the cell of interest.

Results

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 demons...

Discussion

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. Thi...

Materials

Name	Company	Catalog Number	Comments
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