Digital holographic microscopy (DHM) is a volumetric technique that allows imaging samples 50-100X thicker than brightfield microscopy at comparable resolution, with focusing performed post-processing. Here DHM is used for identifying, counting, and tracking microorganisms at very low densities and compared with optical density measurements, plate count, and direct count.
Accurately detecting and counting sparse bacterial samples has many applications in the food, beverage, and pharmaceutical processing industries, in medical diagnostics, and for life detection by robotic missions to other planets and moons of the solar system. Currently, sparse bacterial samples are counted by culture plating or epifluorescence microscopy. Culture plates require long incubation times (days to weeks), and epifluorescence microscopy requires extensive staining and concentration of the sample. Here, we demonstrate how to use off-axis digital holographic microscopy (DHM) to enumerate bacteria in very dilute cultures (100-104 cells/mL). First, the construction of the custom DHM is discussed, along with detailed instructions on building a low-cost instrument. The principles of holography are discussed, and a statistical model is used to estimate how long videos should be to detect cells, based on the optical performance characteristics of the instrument and the concentration of the bacterial solution (Table 2). Video detection of cells at 105, 104, 103, and 100 cells/mL is demonstrated in real time using un-reconstructed holograms. Reconstruction of amplitude and phase images is demonstrated using an open-source software package.
Determination of accurate bacterial counts in very dilute samples is crucial in many applications: a few examples are water and food quality analysis1,2,3; detection of pathogens in blood, cerebrospinal fluid, or sputum4,5; production of pharmaceutical products, including sterile water6; and environmental community analysis in oligotrophic environments such as the open ocean and sediments7,8,9. There is also increasing interest in detection of possible extant microbial life on the icy moons of Jupiter and Saturn, particularly Europa10,11 and Enceladus12,13,14, which are known to have subsurface liquid oceans. Because no mission since Viking in 1978 has attempted to find extant life on another planet, there has been limited development of technologies and instruments for bacterial identification and counting during space missions15.
Traditional methods of plate count find only culturable cells, which can represent a minority of species in environmental strains, sometimes <1%16. Plates require days or weeks of incubation for maximum success, depending upon the strain. Epifluorescence microscopy has largely replaced plate counts as the gold standard for rapid and accurate microbial enumeration. Nucleic-acid-labeling fluorescent dyes such as 4′,6-diamidino-2-phenylindole dihydrochloride (DAPI), SYBR Green, or acridine orange that bind to nucleic acids are the typical dyes used17,18,19, though many studies use fluorescent indicators of Gram sign20,21,22,23,24. Using these methods without pre-concentration steps leads to limits of detection (LoDs) of ~105 cells per mL. Improvements in LoD are possible using filtration. A liquid sample is vacuum-filtered onto a membrane, usually polycarbonate and ideally black to reduce background. Low-background dyes such as the DNA stains mentioned above may be applied directly to the filter25. For accurate counting by eye, ~105 cells are required per filter, which means that for samples more dilute than ~105 cells per mL, significant sample volumes must be collected and filtered. Laser-scanning devices have been developed in order to systematically explore all regions of the filter and thus reduce the number of cells required for counting, pushing limits of detection down to ~102 cells per mL26. However, these are not available in most laboratories, and require sophisticated hardware as well as software that permit expert confirmation that observed particles are bacteria and not debris.
For reference, adults with sepsis usually begin showing symptoms at <100 cells/mL of blood, and infants at <10 cells/mL. A blood draw from an adult takes 10 mL, and from an infant, 1 mL. PCR-based methods are inhibited by the presence of human and non-pathogenic flora DNA and by PCR-inhibiting components in the blood27,28. Despite a variety of emerging techniques, cultures remain the gold standard for the diagnosis of bloodstream infections, especially in more rural areas or developing nations. For detection of life on other planets, thermodynamic calculations can estimate the energy budget for life and thus the expected possible biomass. 1 – 100 cells/mL are expected to be thermodynamically reasonable on Europa29. It can be readily seen from these numbers that detection of very small numbers of cells in large amounts of aqueous solution is an important unsolved problem.
In this paper, we demonstrate detection of Serratia marcescens and Shewanella oneidensis (wild-type and non-motile mutant) at concentrations of 105, 104, 103, and 100 cells/mL using an off-axis digital holographic microscope (DHM). The key advantage of DHM over traditional light microscopy is the simultaneous imaging of a thick sample volume at high resolution—in this implementation, the sample chamber was 0.8 mm thick. These sample chambers were constructed by the soft-lithography of polydimethylsiloxane (PDMS) from a precision-machined aluminum mold with a tolerance of ± 50 µm. This represents an approximately 100-fold improvement in depth of field over high-power light microscopy. DHM also provides quantitative phase information, allowing for measurements of optical path length (product of refractive index and thickness). DHM and similar techniques have been used for monitoring bacterial and yeast cell cycle and calculation of bacterial dry mass30,31,32; scattering differences may even be used to differentiate bacterial strains33.
The instrument we use is custom-built specifically for use with microorganisms, as previously published34,35, and its design and construction are demonstrated and discussed. Aqueous solutions are continuously supplied to a 0.25 µL volume sample chamber via syringe pump; the flow rate is determined by camera frame rate in order to ensure imaging of the entire sample volume. A statistical calculation predicts the number of sample volumes that must be imaged in order to detect a significant number of cells at a given concentration.
For cell-detection applications, reconstruction of the holograms into amplitude and phase images was not required; analysis was performed on the raw hologram. This saves significant computational resources and disk space: a 500 Mb hologram video will be 1 – 2 Tb when reconstructed. However, we do discuss reconstruction through depth of the sample to confirm that the holograms represent the desired species. An important feature of DHM is its ability to monitor both intensity and phase of the images. Organisms that are nearly transparent in intensity (such as most biological cells) appear clearly in phase. As it is a label-free technique, no dyes are used. This is an advantage for possible space flight applications, since dyes may not survive the conditions of a mission and—more importantly—cannot be assumed to work with extraterrestrial organisms, which may not use DNA or RNA for encoding. It is also an advantage for work in extreme environments such as the Arctic and Antarctic, where dyes may be difficult to bring to the remote location and may degrade upon storage. Reconstruction of images into phase and amplitude is performed using an open-source software package that we have made available on GitHub (SHAMPOO) or using ImageJ.
1. Growth and Enumeration of Bacteria
NOTE: This is applicable to almost any bacterial strain grown in the appropriate medium36. In our example, we use three strains: Serratia marcescens as a common, easy identifiable lab strain; and a smaller, highly motile environmental strain, Shewanella oneidensis MR-1. To compare detection of motile vs. non-motile cells, a non-motile Shewanella mutant, Δ FlgM, is also used for comparison37. All strains are grown in lysogeny broth (LB).
2. Preparation of Highly Dilute Samples for DHM
3. Recording DHM Videos
4. Calculation of Cell Density and Limits of Detection
5. Image Reconstruction to Amplitude and Phase
The results should indicate the ability to detect living and dead bacteria at very low levels by DHM. The number of bacteria counted should be consistent with the results obtained using the Petroff-Hauser counting chamber and plate counts. Standard statistical methods provide information about the accuracy of the different detection methods at various bacterial concentrations.
Figure 1 shows the Petroff-Hausser counting chamber used as well as the microstructure of the counting chamber that is used in the direct enumeration of cells. Figure 2A shows a typical enriched culture plate that has been allowed sufficient time for S. marcescens to grow colonies (16 hours at room temperature). Figure 2B shows growth curves for all the strains used. While the exact spectrophotometer measurements vary, the relationship between OD600 and colony count should be linear within the reliable range of the spectrophotometer; we measured ~OD600 = 0.1 as ~108/mL. The direct Petroff-Hausser counts and colony counts should agree to within 10% until the cells begin to die as seen by turnover in the growth curve.
Figure 3 shows the design of the custom sample chambers. Custom microfluidic chambers were used for this experiment to allow the continuous flow of sample through the field of view of the holographic microscope. Standard microfluidic techniques and materials were used in the construction of the sample chambers. Blunt tipped stainless steel needles with male Luer-Lok fittings were inserted into the microfluidic device to provide inlet and outlet ports to the sample chamber. The microfluidic channel geometry was established by the soft lithography of polydimethylsiloxane (PDMS) using a machined aluminum master mold, which was compressed, via aluminum frames, between two optical quality glass windows. The PDMS establishes flow channel geometries as well as secures hypodermic needles, which provide fluid transport to and from the flow channels. These sample chambers use PDMS solely to establish flow channel geometry and thus contain no PDMS in the optical path of the instrument. Note that due to the optical nature of an off-axis holographic microscope, two micro-channels were molded parallel to each other in the sample chamber. This is intended to match optical path lengths between both arms of the interferometer (referred to as the 'specimen' and 'reference' arms). The reference channel should only contain sterile media, and should not be subject to flow.
Figure 4 shows a schematic and images of the digital holographic microscope in its laboratory implementation. It is operated by the software provided with the camera. Figure 5 shows representative DHM data for bacterial enumeration and their correspondence with predicted numbers of cells per field of view. From the raw holograms, without reconstruction, it is possible to see and count the cells by median subtraction and manual tracking. This greatly reduces the computational burden compared to reconstructing all z-planes, as the raw holograms show the cells throughout the depth of the chamber. The appearance of cultures at different concentrations and their correspondence with plate count and Petroff-Hauser count are indicated. Total cell counts are determined by averaging the number of cells seen in each sample volume over the total number of volumes imaged. The limit of detection is determined by the total number of sample volumes imaged. In this experiment, we set a maximum limit on this value, with a maximum total of 10 mL of sample imaged in 0.25 µL increments (a total of 40,000 images per sample). The recordings are stopped before the maximum total volume is achieved if cells are clearly observed in each frame over a series of 20 frames. In a previous paper, we reported a formula for estimating the number of sample volumes needed to detect bacteria at a 50% confidence level at concentrations ranging from 1 – 105 cells/mL (Table 1). Based upon this estimate, an examination of all 40,000 frames should allow for detection of cells at 1/mL, though this calculation is computationally intensive.
Figure 6 shows the use of reconstruction software (Fiji)39. A single hologram is opened and the software "OD-Numerical Propagation-Numerical Diffraction" is run. The imaging parameters are given in the dialog box; amplitude, intensity, and phase may all be reconstructed. A single z plane may be chosen, or batch reconstructions may be used to reconstruct the entire depth.
Figure 7 shows reconstructions of amplitude and phase images at different depths in a selected sample. The Fourier mask used to generate the reconstructions is shown, along with a median-subtracted intensity image, a phase image, a maximum projection through the phase stack, and a volume projection of the intensity stack.
Figure 1: Direct counting of bacterial cells. Bacterial concentration in cells/mL is calculated by counting the total number of bacteria in the counting chamber and scaling that number accounting for the fact that the Petroff-Hauser chamber only contains 20 nL of sample. (A) The Petroff-Hausser counting chamber. (B) Image of counting chamber microstructure under 10X phase contrast, showing boxes of different sizes useful for enumerating cells of different sizes. (C) Appearance of S. marcescens within the counting chamber smallest boxes; 40X phase contrast. (D) Method for counting that excludes cells falling on the upper and right edges (N), but not the lower and left edges (Y). Please click here to view a larger version of this figure.
Figure 2: Counting cells by colony growth and optical density. (A) S. marcescens incubated at room temperature for 16 hours. At the correct dilution, culture plates should have 20 – 200 colonies for accurate counting and statistics. The plate shown is somewhat too crowded for reliable counting. (B) Growth curves for the 3 strains used. An optical density of 0.1 corresponds to ~108 cells/mL, but exact values vary by instrument. Direct Petroff-Hauser counts agree with plate counts to within 10% within the linear growth range. Error bars represent the standard deviation of 5 independent samples. Please click here to view a larger version of this figure.
Figure 3: Sample chamber design. (A) Microfluidic schematic (specimen and reference microchannels labeled). (B) Fully assembled sample chamber. (C) Exploded CAD view showing the constituent elements of the sample chamber. Please click here to view a larger version of this figure.
Figure 4: Schematic and images of the digital holographic microscope. (A) Schematic showing four main elements: the source, the sample (specimen path is labeled Spec. and reference path is labeled Ref.), the microscope, and the sensor. (B) Solid model of the hardware. The fiber-fed source assembly is at the bottom, and the imaging camera is at the top. The microscope optics – comprised of the two aspheric lenses and the relay lens – are contained within the 300 mm long lens tube. The three-axis stage between the source the microscope optics provides easy manual manipulation of the specimen under study. (C) Photograph of the instrument in the laboratory. The scale bar in the bottom portion of the image represents 1 inch. The schematics are re-drawn from Wallace et al.34. Please click here to view a larger version of this figure.
Figure 5: Growth curves and representative DHM data. (A) Predicted cells per field of view based upon a sample volume of 365 µm x 365 µm x 1 mm. The optimal range for DHM enumeration is shown in the rectangle, and real samples are indicated as (B) and (C). (B) DHM hologram of a Shewanella culture at 8 x 105 cells/mL measured by Petroff-Hauser and 7.7 x 105 cells/mL as measured by plate count; this image shows 110 cells, an excellent fit to the prediction. (C) DHM hologram of a Shewanella culture at 2.0 ±0.1 x 105 cells/mL measured by Petroff-Hauser and 2.0 ±0.3 x 105 cells/mL as measured by plate count; this image shows 27 cells. Scale bar = 35 µm. Please click here to view a larger version of this figure.
Figure 6: Performing reconstructions. (A) Open a hologram in Fiji. (B) Run Numerical Propagation and enter the wavelength and image size. Choose a z distance close to 0 to start. (C) Draw a Fourier mask when prompted to select the real or virtual image. (D) Bath processing allows for reconstructions throughout the volume. Please click here to view a larger version of this figure.
Figure 7: Appearance of reconstructions. Images of Serratia culture at ~106 cells/mL. (A) Raw hologram. (B) Fourier transform. The area inside the circle is the area to be used; the rest is masked. (C) Amplitude reconstruction at a single z-plane. (D) Phase reconstruction at a single z-plane; the scale bar shows the phase shift in degrees. (E) Maximum intensity projection through all z-planes in the phase stack. (F) Volume view of the full amplitude stack (365 x 365 x 800 µm). Scale bar = 35 µm. Please click here to view a larger version of this figure.
Concentration [cells per mL] | Frames Needed |
1.00E+00 | 6,686 |
1.00E+01 | 669 |
1.00E+02 | 67 |
1.00E+03 | 7 |
1.00E+04 | 1 |
1.00E+05 | 1 |
Table 1: Minimum number of required sample volumes to detect bacteria by DHM (detection confidence of 50%).
Property | Value | Unit |
Operating Wavelength | 405 | nm |
Objective Numerical Aperture | 0.3 | – |
System Magnification | 20 | – |
Lateral Resolution | 0.7 | μm |
Sample Imaging Volume | 360 x 360 x 800 | μm x μm x μm |
Instrument Length | 400 | mm |
Table 2: Optical and performance specifications for the digital holographic microscope.
Numerical reconstruction of holograms: For the numerical reconstruction of holograms, the angular spectrum method (ASM) is used. This involves the convolution of the hologram with the Green's Function for the DHM. The complex wavefront of the image at a particular focal plane can be calculated by employing the Fourier Convolution Theorem as follows:
(1)
Where is the Fourier Transform operator,
is the hologram matrix, and
is the Green's Function of the DHM, which is defined as:
(2)
Where is the distance, along the optical axis, of the desired focal plane to be reconstructed,
is the illumination wavelength,
and
are the number of pixels in the x and y directions, respectively, and
and
are the pixel sizes (at the detector plane) in the x and y directions, respectively40.
From the complex wave front, the intensity can be calculated by the magnitude squared of , and the phase can be obtained by the arctangent of the quotient between the imaginary and real parts of
.
Modifications and troubleshooting
Experiments involving the enumeration of sparse bacterial concentrations are highly susceptible to contamination, false positives, and false negatives. For this reason, the growth and enumeration of bacteria, as well as the preparation of highly dilute samples for DHM are critical steps in this experiment. Caution must be taken to avoid contamination in this experiment by carefully controlling the environment where the bacteria are grown and enumerated. Proper sterilization techniques, including autoclaving, help prevent against contamination as well as false positives. To further prevent against false positives, bacteria were selected for this experiment that exhibit relatively unique characteristics (e.g., S. marcescens has a very distinct red color when cultured). To prevent against false negatives, the bacteria chosen for this experiment were selected such that they would be able to grow on enriched culture plates, as well as have characteristic shapes and/or motility patterns to be identified via DHM. It is important to periodically test for contamination and/or unexpected cell growth by plating the samples and quantifying the bacteria by plate count.
In conducting this experiment, a data size of roughly 2.4 GB per 1 mL of sample is expected. This number, of course, is dependent on the camera and file format used. The data size reported is by the use of a 4 Mpx detector and a standard TIFF file format. The post-processing of the data requires roughly 6 min of computation per 1 mL of sample. This processing time was observed while using an Intel i7 8-core processor at 3.5 GHz and 32 GB of RAM. These calculations can be parallelized to the GPU for greatly reduced processing times, but was not done here.
Limitations of the technique
While using DHM to detect very low concentrations of bacteria is possible, this method is non-specific. Morphology alone is not a conclusive means of identifying bacterial strains. Diagnostic applications will present specific challenges due to the presence of non-target cells and aggregates. However, the advantage of the dual-beam instrument presented here is that it is able to image bacteria in relatively turbid samples. The ability to detect bacteria in clinical samples, and the degree of pre-processing of the samples that is required, remain to be determined.
Critical steps in the protocol
Step 2: Obtaining accurate dilutions is essential to quantitative enumeration. Any dilutions made should be backed up by plate count. Step 3: The choice of flow rate and corresponding frame rate must be made carefully to ensure that all cells are seen. Adjustment of the pump parameters should be carefully tuned before performing a limit of detection experiment. Step 4: For very low bacterial concentrations, it is crucial to capture enough data to be confident of detection. It may be necessary to reconstruct some data to be certain that cells, and not debris, are being seen. In extreme cases, more than 10 mL of solution may be required. Non-motile strains are more likely to be ambiguous than motile strains. Step 5: If there is any ambiguity in the data, reconstruction can help demonstrate that holograms represent cells and not debris. High-resolution reconstructions will show classic cell morphologies, and motility can clearly be seen in tracks of motile strains.
The authors have nothing to disclose.
The authors acknowledge the Gordon and Betty Moore Foundation Grants 4037 and 4038 to the California Institute of Technology for funding this work.
Bacto Yeast | BD Biosciences | 212750 | |
Bacto Tryptone | BD Biosciences | 211705 | |
Sodium chloride | Sigma-Aldrich | 7710 | Many options for purchase |
Bacto Agar | BD Biosciences | 214010 | |
10 cm Petri dishes | VWR | 10053-704 | |
15 mL culture tubes | Falcon (Corning Life Sciences) | 352002 | Loose-capped |
Petroff-Hauser chamber | Electron Microscopy Sciences | 3920 | |
10 mL syringes | BD Biosciences | 309604 | Luer-Lok tip not necessary |
Male Luer to 1/16” barbed fitting | McMaster-Carr | 51525K291 | |
Autoclavable 1/16” ID PVC tubing for flow | Nalgene | 8000-0004 | Sold by length, purchase accordingly |
Syringe pump | Harvard Apparatus | PHD 2000 | |
Sample Chamber | Custom | n/a | See Materials Section |
Holographic Microscope | Custom | n/a | See Wallace et al. |
Open-source software | Custom | https://github.com/bmorris3/shampoo |