Citation:

Wolf phase tomography (WPT) of transparent structures using partially coherent illumination


  • Light: Science & Applications  9, Article number: (2020)
More Information
  • Corresponding author:
    Gabriel Popescu (gpopescu@illinois.edu)
  • Received: 10 April 2020
    Revised: 30 July 2020
    Accepted: 03 August 2020
    Published online: 19 August 2020

doi: https://doi.org/10.1038/s41377-020-00379-4

  • In 1969, Emil Wolf proposed diffraction tomography using coherent holographic imaging to extract 3D information from transparent, inhomogeneous objects. In the same era, the Wolf equations were first used to describe the propagation correlations associated with partially coherent fields. Combining these two concepts, we present Wolf phase tomography (WPT), which is a method for performing diffraction tomography using partially coherent fields. WPT reconstruction works directly in the space–time domain, without the need for Fourier transformation, and decouples the refractive index (RI) distribution from the thickness of the sample. We demonstrate the WPT principle using the data acquired by a quantitative-phase-imaging method that upgrades an existing phase-contrast microscope by introducing controlled phase shifts between the incident and scattered fields. The illumination field in WPT is partially spatially coherent (emerging from a ring-shaped pupil function) and of low temporal coherence (white light), and as such, it is well suited for the Wolf equations. From three intensity measurements corresponding to different phase-contrast frames, the 3D RI distribution is obtained immediately by computing the Laplacian and second time derivative of the measured complex correlation function. We validate WPT with measurements of standard samples (microbeads), spermatozoa, and live neural cultures. The high throughput and simplicity of this method enables the study of 3D, dynamic events in living cells across the entire multiwell plate, with an RI sensitivity on the order of 10−5.
  • 加载中
  • [1] Liu, P. Y. et al. Cell refractive index for cell biology and disease diagnosis: past, present and future. Lab. Chip. 16, 634–644 (2016).
    [2] Chen, X. & Korotkova, O. Optical beam propagation in soft anisotropic biological tissues. OSA Contin. 1, 1055–1067 (2018). doi: 10.1364/OSAC.1.001055
    [3] Chen, X., Li, J. & Korotkova, O. Light scintillation in soft biological tissues. Waves Random Complex Media 30, 481–489 (2020). doi: 10.1080/17455030.2018.1530814
    [4] Zhu, R. Y. et al. Correlation-induced spectral changes in tissues. Opt. Lett. 36, 4209–4211 (2011). doi: 10.1364/OL.36.004209
    [5] Takabayashi, M. et al. Disorder strength measured by quantitative phase imaging as intrinsic cancer marker in fixed tissue biopsies. PLoS ONE 13, e0194320 (2018). doi: 10.1371/journal.pone.0194320
    [6] Dunn, G. A. & Zicha, D. Dynamics of fibroblast spreading. J. Cell Sci. 108, 1239–1249 (1995). http://www.ncbi.nlm.nih.gov/pubmed/7622607
    [7] Gerber, S. A. et al. Absolute quantification of proteins and phosphoproteins from cell lysates by tandem MS. Proc. Natl Acad. Sci. USA 100, 6940–6945 (2003). doi: 10.1073/pnas.0832254100
    [8] Barer, R. Interference microscopy and mass determination. Nature 169, 366–367 (1952). doi: 10.1038/169541a0
    [9] Davies, H. G. & Wilkins, M. H. F. Interference microscopy and mass determination. Nature 169, 541 (1952). doi: 10.1038/169541a0
    [10] Wang, Z. et al. Tissue refractive index as marker of disease. J. Biomed. Opt. 16, 116017 (2011). doi: 10.1117/1.3656732
    [11] Giannios, P. et al. Complex refractive index of normal and malignant human colorectal tissue in the visible and near-infrared. J. Biophotonics 10, 303–310 (2017). doi: 10.1002/jbio.201600001
    [12] Subramanian, H. et al. Partial-wave microscopic spectroscopy detects subwavelength refractive index fluctuations: an application to cancer diagnosis. Opt. Lett. 34, 518–520 (2009). doi: 10.1364/OL.34.000518
    [13] Uttam, S. et al. Early prediction of cancer progression by depth-resolved nanoscale mapping of nuclear architecture from unstained tissue specimens. Cancer Res. 75, 4718–4727 (2015). doi: 10.1158/0008-5472.CAN-15-1274
    [14] Park, Y. et al. Refractive index maps and membrane dynamics of human red blood cells parasitized by Plasmodium falciparum. Proc. Natl Acad. Sci. USA 105, 13730–13735 (2018). http://www.pnas.org/content/105/37/13730.abstract
    [15] Bista, R. K. et al. Quantification of nanoscale nuclear refractive index changes during the cell cycle. J. Biomed. Opt. 16, 070503 (2011). doi: 10.1117/1.3597723
    [16] Yanina, I. Y., Lazareva, E. N. & Tuchin, V. V. Refractive index of adipose tissue and lipid droplet measured in wide spectral and temperature ranges. Appl. Opt. 57, 4839–4848 (2018). doi: 10.1364/AO.57.004839
    [17] Kim, D. et al. Holotomography: refractive index as an intrinsic imaging contrast for 3-D label-free live cell imaging. Preprint at https://doi.org/10.1101/106328 (2018).
    [18] Popescu, G. Quantitative Phase Imaging of Cells and Tissues (McGraw-Hill Education, New York, 2011).
    [19] Schürmann, M. et al. Three-dimensional correlative single-cell imaging utilizing fluorescence and refractive index tomography. J. Biophotonics 11, e201700145 (2018). doi: 10.1002/jbio.201700145
    [20] Zhang, Q. N. et al. Quantitative refractive index distribution of single cell by combining phase-shifting interferometry and AFM imaging. Sci. Rep. 7, 2532 (2017). doi: 10.1038/s41598-017-02797-8
    [21] Park, Y. K., Depeursinge, C. & Popescu, G. Quantitative phase imaging in biomedicine. Nat. Photonics 12, 578–589 (2018). doi: 10.1038/s41566-018-0253-x
    [22] Depeursinge, A. et al. Three-dimensional solid texture analysis in biomedical imaging: review and opportunities. Med. Image Anal. 18, 176–196 (2014). doi: 10.1016/j.media.2013.10.005
    [23] Merola, F. et al. Tomographic flow cytometry by digital holography. Light.: Sci. Appl. 6, e16241 (2017). doi: 10.1038/lsa.2016.241
    [24] Robles, F. E. et al. Molecular imaging true-colour spectroscopic optical coherence tomography. Nat. Photonics 5, 744–747 (2011). doi: 10.1038/nphoton.2011.257
    [25] Choi, W. et al. Tomographic phase microscopy. Nat. Methods 4, 717–719 (2007). doi: 10.1038/nmeth1078
    [26] Cotte, Y. et al. Marker-free phase nanoscopy. Nat. Photonics 7, 113–117 (2013). doi: 10.1038/nphoton.2012.329
    [27] Lee, K. R. & Park, Y. K. Quantitative phase imaging unit. Opt. Lett. 39, 3630–3633 (2014). doi: 10.1364/OL.39.003630
    [28] Wang, Z. et al. Spatial light interference microscopy (SLIM). Opt. Express 19, 1016–1026 (2011). doi: 10.1364/OE.19.001016
    [29] Nguyen, T. H. et al. Gradient light interference microscopy for 3D imaging of unlabeled specimens. Nat. Commun. 8, 210 (2017). doi: 10.1038/s41467-017-00190-7
    [30] Bhaduri, B. et al. Diffraction phase microscopy with white light. Opt. Lett. 37, 1094–1096 (2012). doi: 10.1364/OL.37.001094
    [31] Kuś, A. et al. Tomographic phase microscopy of living three-dimensional cell cultures. J. Biomed. Opt. 19, 046009 (2014). http://europepmc.org/abstract/med/24723114
    [32] Rappaz, B. et al. Measurement of the integral refractive index and dynamic cell morphometry of living cells with digital holographic microscopy. Opt. Express 13, 9361–9373 (2005). doi: 10.1364/OPEX.13.009361
    [33] Lue, N. et al. Live cell refractometry using microfluidic devices. Opt. Lett. 31, 2759–2761 (2006). doi: 10.1364/OL.31.002759
    [34] Wolf, E. Three-dimensional structure determination of semi-transparent objects from holographic data. Opt. Commun. 1, 153–156 (1969). doi: 10.1016/0030-4018(69)90052-2
    [35] Chen, M. et al. Multi-layer Born multiple-scattering model for 3D phase microscopy. Optica 7, 394–403 (2020). doi: 10.1364/OPTICA.383030
    [36] B., V. et al. Integrated dual-tomography for refractive index analysis of free-floating single living cell with isotropic superresolution. Sci. Rep. 8, 5943 (2018).
    [37] Lue, N. et al. Synthetic aperture tomographic phase microscopy for 3D imaging of live cells in translational motion. Opt. Express 16, 16240–16246 (2008). doi: 10.1364/OE.16.016240
    [38] Kim, T. et al. White-light diffraction tomography of unlabelled live cells. Nat. Photonics 8, 256–263 (2014). doi: 10.1038/nphoton.2013.350
    [39] Phillips, K. G., Jacques, S. L. & McCarty, O. J. T. Measurement of single cell refractive index, dry mass, volume, and density using a transillumination microscope. Phys. Rev. Lett. 109, 118105 (2012). doi: 10.1103/PhysRevLett.109.118105
    [40] Paganin, D. & Nugent, K. A. Noninterferometric phase imaging with partially coherent light. Phys. Rev. Lett. 80, 2586–2589 (1998). doi: 10.1103/PhysRevLett.80.2586
    [41] Leblanc-Hotte, A. et al. On-chip refractive index cytometry for whole-cell deformability discrimination. Lab Chip 19, 464–474 (2019). doi: 10.1039/C8LC00938D
    [42] Liang, X. J. et al. Determining refractive index of single living cell using an integrated microchip. Sens. Actuators A: Phys. 133, 349–354 (2007). doi: 10.1016/j.sna.2006.06.045
    [43] Wang, Z. et al. Spatial light interference tomography (SLIT). Opt. Express 19, 19907–19918 (2011). doi: 10.1364/OE.19.019907
    [44] Kim, K. et al. Real-time visualization of 3-D dynamic microscopic objects using optical diffraction tomography. Opt. Express 21, 32269–32278 (2013). doi: 10.1364/OE.21.032269
    [45] Macias-Garza, F., Diller, K. R. & Bovik, A. C. Missing cone of frequencies and low-pass distortion in three-dimensional microscopic images. Optical Eng. 27, 276461 (1988). http://adsabs.harvard.edu/abs/1988OptEn..27..461M
    [46] Habaza, M. et al. Tomographic phase microscopy with 180° rotation of live cells in suspension by holographic optical tweezers. Opt. Lett. 40, 1881–1884 (2015). doi: 10.1364/OL.40.001881
    [47] Habaza, M. et al. Rapid 3D refractive-index imaging of live cells in suspension without labeling using Dielectrophoretic cell rotation. Adv. Sci. 4, 1600205 (2017). doi: 10.1002/advs.201600205
    [48] Charrière, F. et al. Living specimen tomography by digital holographic microscopy: morphometry of testate amoeba. Opt. Express 14, 7005–7013 (2006). doi: 10.1364/OE.14.007005
    [49] Memmolo, P. et al. 3D morphometry of red blood cells by digital holography. Cytom. Part A 85, 1030–1036 (2014). doi: 10.1002/cyto.a.22570
    [50] Chen, X. & Korotkova, O. Phase structuring of 2D complex coherence states. Opt. Lett. 44, 2470–2473 (2019). doi: 10.1364/OL.44.002470
    [51] Mandel, L. & Wolf, E. Coherence properties of optical fields. Rev. Mod. Phys. 37, 231–287 (1965). doi: 10.1103/RevModPhys.37.231
    [52] Kandel, M. E. et al. Real-time halo correction in phase contrast imaging. Biomed. Opt. Express 9, 623–635 (2018). doi: 10.1364/BOE.9.000623
    [53] Shan, M. G., Kandel, M. E. & Popescu, G. Refractive index variance of cells and tissues measured by quantitative phase imaging. Opt. Express 25, 1573–1581 (2017). doi: 10.1364/OE.25.001573
    [54] Ryu, D. H. et al. Deep learning-based optical field screening for robust optical diffraction tomography. Sci. Rep. 9, 15239 (2019). doi: 10.1038/s41598-019-51363-x
    [55] Wang, H. D. et al. Deep learning enables cross-modality super-resolution in fluorescence microscopy. Nat. Methods 16, 103–110 (2019). doi: 10.1038/s41592-018-0239-0
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Figures(6)

Article Metrics

Article views(309) PDF downloads(95) Citation(0) Citation counts are provided from Web of Science. The counts may vary by service, and are reliant on the availability of their data.

Wolf phase tomography (WPT) of transparent structures using partially coherent illumination

  • Quantitative Light Imaging Laboratory, Beckman Institute for Advanced Science and Technology, Department of Electrical and Computer Engineering, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
  • Corresponding author:

    Gabriel Popescu, gpopescu@illinois.edu

doi: https://doi.org/10.1038/s41377-020-00379-4

Abstract: In 1969, Emil Wolf proposed diffraction tomography using coherent holographic imaging to extract 3D information from transparent, inhomogeneous objects. In the same era, the Wolf equations were first used to describe the propagation correlations associated with partially coherent fields. Combining these two concepts, we present Wolf phase tomography (WPT), which is a method for performing diffraction tomography using partially coherent fields. WPT reconstruction works directly in the space–time domain, without the need for Fourier transformation, and decouples the refractive index (RI) distribution from the thickness of the sample. We demonstrate the WPT principle using the data acquired by a quantitative-phase-imaging method that upgrades an existing phase-contrast microscope by introducing controlled phase shifts between the incident and scattered fields. The illumination field in WPT is partially spatially coherent (emerging from a ring-shaped pupil function) and of low temporal coherence (white light), and as such, it is well suited for the Wolf equations. From three intensity measurements corresponding to different phase-contrast frames, the 3D RI distribution is obtained immediately by computing the Laplacian and second time derivative of the measured complex correlation function. We validate WPT with measurements of standard samples (microbeads), spermatozoa, and live neural cultures. The high throughput and simplicity of this method enables the study of 3D, dynamic events in living cells across the entire multiwell plate, with an RI sensitivity on the order of 10−5.

Introduction
  • The refractive index (RI) is a fundamental physical property that determines how light interacts with a medium in terms of scattering, governed by its real part, and absorption, through its imaginary part1-5. In biological applications, the RI distribution correlates strongly with cellular properties, such as dry mass and chemical concentrations6-9. Tissue RI can also act as an intrinsic marker for cancer diagnosis10, 11. Nanoscale morphological changes in cells and tissues can be revealed by RI maps12, 13. For example, it has been shown that cancer tissue exhibits higher RI variances than normal tissue10, 11. The RI can also be used to study biological dynamics, including cellular transport and mitosis14, 15, and can be used for phenotypic screening and cellular monitoring16, 17. To obtain the RI distribution of cells and tissues from the measured field properties in different imaging modalities, one must go beyond the typical quantities measured in phase imaging and solve an inverse scattering problem. A condition for this problem to yield unique solutions is to measure the full information about the scattered field, meaning both the amplitude and phase. Interferometric microscopy provides a method for phase retrieval in weakly scattered samples such as cells and tissues18-20.

    Quantitative-phase imaging (QPI) has emerged as a growing field focused on extracting the phase distributions of an imaging field and exploiting that information for biomedical applications18, 21-27. White-light-based methods, such as spatial light interference microscopy (SLIM)28, gradient light interferometric microscopy29, and (white-light) diffraction-phase microscopy30, can render phase images of live cells without the speckle noise typically associated with coherent illumination21. As a result, the spatial sensitivity to pathlength changes is very high. The optical pathlength measurement depends on both the RI and the thickness of the sample31. To estimate the 2D (axially averaged) RI from the optical pathlength, the thickness distribution of the structures needs to be known or decoupled from the optical pathlength32, 33. However, the accuracy is low due to the geometrical optics approximation, and the results provide only a 2D map of the longitudinally averaged RI.

    For inferring the 3D RI distribution from QPI data, several approaches have been proposed based on solving the deterministic wave equations34, 35. One of them is the filtered back-projection algorithm, which uses the Fourier diffraction theorem and the first-order Born or Rytov approximation25. It connects the object function with the Fourier transform of the projection. The reconstruction of the RI distribution is obtained by combining the frequency bands with respect to different angles36. One can achieve this by rotating the illumination angles or measuring a set of image fields at successive points across a cell with focused beam illumination, also known as synthetic aperture tomography-phase microscopy37. However, the Fourier diffraction theorem34 assumes plane-wave illumination, and is only an approximation for partially coherent fields. To obtain a more accurate solution to the inverse problem, the coherence properties must be taken into account. White-light diffraction tomography (WDT) uses the temporal correlation and instrument response to perform deconvolution on the complex field data to extract the 3D scattering potential38. However, WDT requires a priori knowledge of the instrument impulse response (or transfer function), which is often limited. At the same time, deconvolution operations are usually time-consuming and sensitive to noise.

    The transport of the intensity equation can connect the RI to the intensity of bright-field images39, 40. However, this method is only applicable under a paraxial approximation (see, e.g., Section 12.1 in ref. 18). Integrated microchips have been used to measure RI information, combining an external cavity laser, microlenses, and microfluidic channels into a monolithic device41, 42. Such a device can determine the average RI of a single live cell in real time, but cannot render the RI distribution. 3D RI distributions are usually constructed through axial scanning (z scanning)29, 38, 43 or projections from different angles (computed tomography)25, 26, 44. To increase the axial resolution in 3D reconstruction, efforts have been devoted toward alleviating the incomplete frequency coverage of imaging systems, or "the missing cone problem"45. Illumination angle scanning and rotation of the sample can help fill in a missing cone region. Cells can be rotated by optical tweezers or dielectrophoretic forces in microfluidics46-49. However, these methods involve more complicated procedures.

    In this paper, we propose a fast 3D RI construction method based on the Wolf equations for propagating correlations of partially coherent light50, 51. This approach, referred to as Wolf phase tomography (WPT), involves minimal computational steps, and renders high-resolution RI tomograms without time-consuming deconvolution operations. WPT decouples the RI distribution from the thickness of the sample in the space–time domain directly, without the need for Fourier transformation. We demonstrate that from three independent intensity measurements corresponding to each phase shift, the RI distribution can be reconstructed directly from the Laplacian and second time derivative of the complex correlation functions. We demonstrate WPT with standard polystyrene beads, fixed spermatozoa, and dynamic live-cell imaging over many hours. Interestingly, we find that WPT is able to extract intrinsic RI changes in live cells with a sensitivity on the order of 10−5, which can indicate cell viability in screening applications.

Results
  • WPT relies on a commercial phase-contrast microscope upgraded with a spatial light modulator (SLM) conjugated to the pupil plane. In our implementation, this hardware is provided by a SLIM module (SLIM Pro, Phi Optics), as shown in Fig. 1a. Figure 1b–d shows the temporal spectrum and autocorrelation properties of the illumination (white-light) field. In addition to the π/2 phase shift between the incident and scattered fields introduced by the objective-phase ring, the SLIM module provides further phase shifts with π/2 increments. At the camera plane, we record three intensity images, corresponding to each phase shift, as illustrated in Fig. 1e, namely

    Fig. 1  Optical setup and working principle of Wolf phase tomography (WPT).

    a WPT optical setup. The spatial light interference microscopy (SLIM) add-on module is mounted to the output camera port of a phase-contrast microscope. This module shifts the phase between the incident and scattered field with π/2 increments using the SLM. The interference patterns are recorded by the camera. b Spectrum, c autocorrelation, and d second-order derivative of the autocorrelation of the halogen source measured by the spectrometer. e Three phase-shifted frames of 2-μm polystyrene beads are acquired using the SLIM module (×63/1.4 NA)
    $$ I_d({\bf{r}}) = I_i({\bf{r}}) + I_S({\bf{r}}) + 2\Re [{\Gamma} _{is}\left( {\left\langle \omega \right\rangle \tau _d + {\Delta} \phi ({\bf{r}})} \right)] $$ (1)

    where $\left\langle \omega \right\rangle \tau _d = - d\pi /2$, d = 1, 2, 3, $\left\langle \omega \right\rangle$ is the central frequency of the incident field, $\Re$ stands for the real part, ${\Delta} \phi$ is the phase difference between the incident field Ui and scattered field US, and ${\Gamma} _{pq}({\bf{r}}_1, {\bf{r}}_2, \tau) = \langle U_p^ \ast ({\bf{r}}_1, t)U_q({\bf{r}}_2, t + \tau)\rangle _t$, $p, q = \{ i, s\}.$ From these three frames, we solve for $\Re \left[ {{\Gamma} _{is}(\left\langle \omega \right\rangle \tau _d + {\Delta} \phi ({\bf{r}}))} \right]$. Based on partially coherent light propagation, governed by the Wolf equations51, the RI of the object can be obtained by (see the full derivation in Supplementary Note 1)

    $$ n({\bf{r}}) = \sqrt {\frac{{m({\bf{r}}) - n_0^2\left[ {1 - g({\bf{r}})} \right]}}{{1 + g({\bf{r}})}}} $$ (2)

    In Eq. (2), the functions m and g are defined as

    $$ m({\bf{r}}) = \left. {\frac{{c^2\left( {\nabla ^2\Re [{\Gamma} _{is}({\bf{r}},{\bf{r}},\tau )] + \zeta ({\bf{r}})} \right)}}{{\frac{{\partial ^2\Re [{\Gamma} _{is}({\bf{r}},{\bf{r}},\tau )]}}{{\partial \tau ^2}}}}} \right|_{\tau = - \pi /\left\langle \omega \right\rangle } $$ (3$a$)
    $$ \zeta ({\bf{r}}) = - 2\Re \mathop {\int}\limits_0^\infty {\langle \nabla U_i^ \ast ({\bf{r}},\omega ) \cdot \nabla U_s({\bf{r}},\omega )\rangle e^{i\omega \pi /\left\langle \omega \right\rangle }d\omega } $$ (3$b$)
    $$ g({\bf{r}}) = \left. {\frac{{\frac{{\partial ^2\Re [{\Gamma} _{ii}({\bf{r}},{\bf{r}},\tau )]}}{{\partial \tau ^2}}}}{{\frac{{\partial ^2\Re [{\Gamma} _{is}({\bf{r}},{\bf{r}},\tau )]}}{{\partial \tau ^2}}}}} \right|_{\tau = - \pi /\left\langle \omega \right\rangle } $$ (3$c$)

    where r = (x, y, z) is the spatial coordinate, n0 is the RI of the background media, and c is the speed of light in vacuum. The detailed steps for calculating the terms in Eqs. (3a)–(3c) are given in Supplementary Note 1. The term in Eq. (3b) does not substantially contribute to the final RI and can be omitted for faster construction (see the discussion in Supplementary Note 2). Figure 1b describes the normalized spectrum of a halogen source measured by the spectrometer (ocean optics). The real part of the normalized autocorrelation $\Re [{\Gamma} _{ii}({\bf{r}}, {\bf{r}}, \omega _0\tau)]$ is obtained by taking the Fourier transform of the spectrum (see Fig. 1c). To retrieve the temporal correlation function quantitatively, we normalized the Γii(r, r, 0) value from the spectrometer data to the background intensity from the camera, and corrected it with the spectrally dependent quantum efficiency of the camera. Thus, we ensured that the autocorrelations Γii(r, r, 0) measured by the two different devices have the same value. The second-order time derivative of Γii(r, r, τ) is depicted in Fig. 1d. The Laplacian in Eq. (3a) is calculated using three images with a first-order finite-difference approximation. The z component of the Laplacian was computed using three axially distributed frames separated by a distance that matches the xy pixel sampling, and is much smaller than the diffraction spot. For example, for a ×40/0.75 NA objective, this distance is 0.14 μm, while the diffraction-limited resolution is 0.4 μm. The second-order derivatives in Eqs. (3a) and (3c) are calculated in MATLAB using three phase-shifted frames. Smaller phase shifts would give more accurate derivatives. However, the contrast between different frames would greatly decrease; thus, the signal-to-noise ratio would decrease as well, resulting in a lower accuracy for the derivative. Therefore, to increase the signal-to-noise ratio and accuracy, we keep the phase increment at π/2. This algorithm requires 40 ms to reconstruct the RI map at one z position with a 3-megapixel field of view (MATLAB, i7-8650U CPU).

  • To validate the capability of WPT in extracting the RI distribution, we imaged 2-μm polystyrene microsphere (Polysciences Inc.) z stacks with an RI value of 1.59 at the central wavelength. The beads are suspended in immersion oil (Zeiss) with an RI value of 1.518. Figure 2a shows the three frames corresponding to the different phase shifts of the polystyrene beads. For these experiments, we use a ×63/1.4 NA objective. The real parts of the correlation function Γis at different time-lapse values are illustrated in Fig. 2b. The RI distribution of the microspheres for each z slice is reconstructed via Eq. (2). The 3D rendering of the RI distribution described in Fig. 2c was obtained in AMIRA (Thermo Fisher Scientific). The reconstructed RI value of the microspheres agrees well with the expected value of 1.59 at the central wavelength. The halo artifacts associated with phase-contrast and SLIM images were removed from the final RI maps using our previously reported algorithm52.

    Fig. 2  Wolf phase tomography (WPT) on standard samples.

    a Three phase-shifted frames of 2-μm polystyrene beads suspended in oil are imaged by spatial light interference microscopy (SLIM) using a ×63/1.4 NA objective. b The real part of the correlation function at three different time delays is obtained by solving Eq. (1). c 3D refractive index (RI) tomogram of the 2-μm polystyrene bead
  • A 3D rendering of a bovine sperm cell is displayed in Fig. 3a (Supplementary Video 1). In the sperm head, the acrosome and the nucleus can be identified with RI values between 1.35 and 1.37. The centriole and mitochondria-rich midpiece of the sperm cell yield high RI values (Fig. 3b). The tail of the sperm has an RI value of 1.35, and the axial filament inside the tail, with a slightly higher RI value of 1.36, can be recognized. The end piece of the sperm has the lowest RI value, ~1.34.

    Fig. 3  Wolf phase tomography (WPT) of sperm cells.

    a 3D refractive index (RI) tomogram of a spermatozoon (×40/0.75 NA objective). b xy-plane projection view. The nucleus, acrosome, centriole, and axial filament of the sperm cell are indicated by white arrows. c xz-plane projection view. d Histogram of the RI of the sperm cell
  • Applying the WPT principle, the three frames of hippocampal neurons and their correlation functions are depicted in Fig. 4a and b. The reconstructed RI distribution and 3D rendering of the neurons (Supplementary Video 2) are displayed in Fig. 4c and d. The more detailed structures of individual hippocampal neurons (Supplementary Videos 3 and 4) are illustrated in Fig. 4e and f. The rendering in this case used two colormaps, as shown in Fig. 4e and f. The neuron dendrites have an RI value of ~1.34, while the cell body ranges from 1.35 to 1.38, with a nucleolus of 1.39–1.4. The axon can be recognized in Fig. 4f, as the morphology shows a longer and thinner filamentous structure.

    Fig. 4  Wolf phase tomography (WPT) of neurons.

    a Three phase-shifted frames of hippocampal neurons (×40/0.75 NA objective). b The real part of the correlation function at three different time lags is solved with Eq. (1). c Refractive index (RI) map of hippocampal neurons. df 3D rendering of RI tomograms of hippocampal neurons. e, f Two colormaps are used as indicated to enhance the dendrites and axons. The axon is indicated with a red arrow
  • Due to its high throughput, low phototoxicity, absence of photobleaching, and easy sample preparation, WPT is capable of studying real-time volumetric biological events in living cells. We imaged the growth and proliferation of hippocampal neurons over the course of several days in six-well plates typical in phenotypic screening applications. The RI distribution of the whole well of neurons is displayed in Fig. 5a (Supplementary Video 5). One tile zoom-in of the whole well and its distribution of RI is shown in Fig. 5b (Supplementary Video 6). Figure 5c describes the averages of the RI values within this tile versus time. The average RI values increase with time due to neuron growth. Figure 5c illustrates the average RI of the whole tile, including the neurons and background. As the neurons grow, more pixels in the region of interest appear with higher RIs; thus, the average RI becomes larger. Another point worth mentioning is that the range of the y axis in Fig. 5c is from 1.34045 to 1.34070. Thus, due to the averaging over the large field of view, the change in the RI value detected by our system is at the fifth decimal place, indicating the high sensitivity of WPT. Figure 5d shows that the variance of the RI for this tile increases with time as well53. Note that the range of RI variance values is on the order of 10−6, which is detectable due to the sensitivity conferred by the common-path stability and lack of speckles in SLIM.

    Fig. 5  Dynamic Wolf phase tomography (WPT) of live cells across multiwell plates.

    a The refractive index (RI) map across a whole well of living hippocampal neurons (×10/0.3 NA objective) is composed of 20 × 21 mosaic tiles, each 214 × 204 µm2 in area. b Enlarged RI map of the purple box in (a) with the average (c) and variance (d) of the RI versus time. e Enlarged RI map of the area in the red box in (b) with the average (f) and variance (g) of the RI. The green arrow indicates the increase in the RI when the two neurons separated and their dendrites appeared, and the red arrow shows the decrease in RI when the two neurons died. h Enlarged RI map of the area in the yellow box in (b) with the average (i) and variance (j) of the RI. The green arrow indicates the change in the RI when the dendrites appeared, and the red arrow indicates the decrease in RI when the neuron died

    Figure 5e is the enlarged image of the area in the red box in Fig. 5b containing two neurons. The neurons spread out into two regions at ~t = 16 h, continued growing until ~t = 53 h, and then died. We can see that both the average and variance of the RI show three different stages (Fig. 5f, g). One significant change in the average and variance of the RI appeared when the two neurons separated (red arrows). Another change is visible when the two neurons died (green arrows). The death event was accompanied by a decrease in the mean RI, likely due to the membrane permeability, which allowed for water influx.

    Figure 5h is a magnified image of the area in the yellow box in Fig. 5b containing one neuron. The neuron dendrites started to appear at ~t = 13 h, resulting in a jump in the average RI (Fig. 5i). The neuron continued growing until ~t = 62 h and then died, leading to a decrease in the average RI. Some oscillations in the variance (Fig. 5j) of the RI appeared before the neuron died and exhibited a clear change after the neuron died. Figure 6 demonstrates the capability of WPT for 3D real-time live-cell imaging. The changes in the morphology of the neuron can be recognized at different time frames.

    Fig. 6  Time-lapse Wolf phase tomography (WPT) of live neurons.

    a Refractive index (RI) map of a live hippocampal neuron imaged with a ×40/0.75 NA objective. be 3D RI tomograms of the hippocampal neuron 12-min apart
Discussion
  • In summary, we proposed a new high-throughput RI tomography method, WPT, based on the correlation propagation of partially coherent light. We demonstrated the capability of WPT with tomographic reconstructions of standard polystyrene beads, spermatozoa, and hippocampal neurons. Our method builds on the coherence theory pioneered by Emil Wolf by combining the Wolf equations and diffraction tomography to perform the reconstruction directly in the space–time domain, without the need for Fourier transformation. WPT decouples the RI distribution from the thickness of the object by calculating the Laplacian and the second-order time derivative of the complex correlation functions. High-resolution tomograms of RI distributions are acquired using a z stack of three phase-shifted intensity frames. As a result, the tomographic reconstruction is very fast, requiring only 40 ms per z slice. In our implementation, the total acquisition time is 180 ms per z slice. With more advanced SLMs, z-scanning stages, and cameras, ~40 ms of the total acquisition time can be achieved per z slice. WPT has a high RI sensitivity, on the order of 10−5, which is useful as an intrinsic marker for live-cell monitoring. We illustrated this ability by imaging dynamic live cells over many hours. As a label-free method, WPT is nondestructive and is not limited by the photobleaching and phototoxicity commonly associated with fluorescence microscopy.

    In a larger context, WPT highlights the advantage of partially coherent illumination, phase shifting, and phase-contrast geometry. The white-light illumination and common-path interferometry allow for speckle-free and nanometer-path-length stability. As phase shifting is performed in the pupil plane, in the time domain, our reconstruction preserves the diffraction-limited resolution of the microscope without introducing coherent artifacts (such as residual fringes)54.

    Although the mathematical model derived in WPT does not use approximation, the ability to access the information of the correlation functions is limited experimentally to weakly scattering samples. It is also worth mentioning that since we use white light with a very short coherence length (~1–2 µm), long optical paths are cut off via coherence gating. This implies that within the optical section, scattering is characterized well by the Born approximation, even though the entire specimen might produce multiple scattering. However, WPT is still limited by the "missing cone problem" due to its finite numerical aperture. Therefore, WPT can also adopt other methods, such as cell rotation and illumination angle scanning, to achieve better axial resolution and depth sectioning. As an alternative, new advances in deep learning appear promising in addressing this issue by frequency extrapolation55. We envision that WPT will find important applications in material and life sciences, such as studying adherent cell growth with the segmentation of a nucleus compared to the whole cell, cell classification based on RI, histopathology for cancer diagnosis based on RI10, and 3D tracking of collagen fibers.

Materials and methods
  • The SLIM add-on module is mounted to the output camera port of a commercial phase-contrast microscope. The module contains an SLM (Meadowlark, XY Series) and a camera (Hamamatsu, ORCA-Flash 4.0 V2). Measurements were conducted using an inverted microscope (Zeiss, Axio Observer Z1) with a halogen light source (Zeiss, HAL 100). Cells were imaged with an incubation system under ×63/1.4 NA, ×40/0.75 NA, and 10×0.3 NA objectives with matching phase-contrast illumination. The sampling was uniform in the x, y, and z directions. For example, for a ×40/0.75 NA objective, the distance is 0.14 μm, which is smaller than the diffraction-limited resolution of 0.4 μm.

    The 40-ms reconstruction per frame is faster than the SLIM image acquisition rate of 180 ms, which requires 30 ms for SLM stabilization (Meadowlark XY Series), 10 ms for exposure (Hamamatsu, V2 Orca Flash), and 60 ms for z scanning. We expect the technique to be able to be implemented on faster hardware without modification.

  • The hippocampal neurons were prepared as follows: primary neurons were harvested from dissected hippocampi of Sprague–Dawley rat embryos. The hippocampi were dissociated with the enzyme to obtain hippocampal neurons. The hippocampal neurons were then plated onto a six-well plate that was precoated with poly-D-lysine (0.1 mg/ml, Sigma-Aldrich). The hippocampal neurons were initially incubated with plating medium containing 86.55% MEM Eagle's with Earle's BSS (Lonza), 10% fetal bovine serum (refiltered, heat-inactivated, Thermo Fisher), 0.45% of 20% (wt./vol.) glucose, 1 × 100 mM sodium pyruvate (100×, Sigma-Aldrich), 1 × 200 mM glutamine (100×, Sigma-Aldrich), and 1× penicillin/streptomycin (100×, Sigma-Aldrich) to facilitate the attachment of neurons (300 cells/mm2). After 3 h of incubation in an incubator (37 ℃ and 5% CO2), the plating medium was aspirated and replaced with maintenance medium containing NeurobasalTM growth medium supplemented with B-27 (Thermo Fisher), 1% 200 mM glutamine (Thermo Fisher), and 1% penicillin/streptomycin (Thermo Fisher) at 37 ℃ in the presence of 5% CO2. The hippocampal neurons in Fig. 5 were grown for 2 days in vitro, and dynamic images were taken for 4 days. The hippocampal neurons in Fig. 6 were grown for 14 days in vitro, and snapshots were taken every 12 min. The hippocampal neurons in Fig. 4 were fixed. The sperm cells in Fig. 3 were fixed in 10% paraformaldehyde.

Acknowledgements
  • This work is supported by the National Science Foundation (CBET0939511 STC, NRT-UtB 1735252), the National Institute of General Medical Sciences (GM129709), and the National Cancer Institute (CA238191).

Author contributions
  • G.P. and X.C. proposed the idea. X.C. solved the inverse problem and derived the model. M.E.K. wrote the control software. X.C., M.E.K., and C.H. performed the imaging. X.C. wrote the reconstruction algorithm and analyzed the data. Y.L. prepared the hippocampal neurons. X.C., M.E.K., and G.P. wrote the paper. G.P. supervised the work.

Data availability
  • The data that support the findings in this study are available from the corresponding author upon reasonable request.

Code availability
  • The code and computer algorithms that support the findings of this study are available from the corresponding author upon reasonable request.

Conflict of interest
  • G.P. has a financial interest in Phi Optics, Inc., a company developing quantitative-phase-imaging technology for materials and life science applications. The remaining authors declare that they have no conflict of interest.

Supplementary information
Reference (55)

Catalog

    /

    DownLoad:  Full-Size Img PowerPoint
    Return
    Return