Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Biomass measurements of single neurites in vitro using optical wavefront microscopy

Open Access Open Access

Abstract

Quantitative phase microscopies (QPMs) enable label-free, non-invasive observation of living cells in culture, for arbitrarily long periods of time. One of the main benefits of QPMs compared with fluorescence microscopy is the possibility to measure the dry mass of individual cells or organelles. While QPM dry mass measurements on neural cells have been reported this last decade, dry mass measurements on their neurites has been very little addressed. Because neurites are tenuous objects, they are difficult to precisely characterize and segment using most QPMs. In this article, we use cross-grating wavefront microscopy (CGM), a high-resolution wavefront imaging technique, to measure the dry mass of individual neurites of primary neurons in vitro. CGM is based on the simple association of a cross-grating positioned in front of a camera, and can detect wavefront distortions smaller than a hydrogen atom (∼0.1 nm). In this article, an algorithm for dry-mass measurement of neurites from CGM images is detailed and provided. With objects as small as neurites, we highlight the importance of dealing with the diffraction rings for proper image segmentation and accurate biomass measurements. The high precision of the measurements we obtain using CGM and this semi-manual algorithm enabled us to detect periodic oscillations of neurites never observed before, demonstrating the sufficient degree of accuracy of CGM to capture the cell dynamics at the single neurite level, with a typical precision of 2%, i.e., 0.08 pg in most cases, down to a few fg for the smallest objects.

© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Quantitative phase microscopies (QPMs) are optical imaging techniques capable of mapping the phase of a light wave after passing the object plane of a microscope [1]. Most of the applications of QPMs are related to biology, for the observation of cells in vitro. In this context, the measured phase image $\varphi $ reads

$$\varphi ({x,y} )= \frac{{2\pi }}{\lambda }({n - {n_0}} )h({x,y} )$$
$$\varphi ({x,y} )= \frac{{2\pi }}{\lambda }\delta \ell ({x,y} )$$
where h is the thickness profile of the imaged object, n its refractive index (supposed here to be uniform for the sake of simplicity) and ${n_0}$ the refractive index of the environment. In Eq. (2), $({n - {n_0}} )h$ has been replaced with $\delta \ell $, called the optical path difference (OPD) image, or equivalently the optical wavefront profile.

QPMs have the great advantage, compared with fluorescence microscopies, to be label-free and to not require the use of UV-blue light. For these reasons, they are non-invasive, do not suffer from photobleaching, and can be conducted for arbitrarily long periods of time, without decay of the signal and without perturbing the imaged cells. Moreover, the quantitative nature of the measurements provides QPMs with the ability to measure the bio-mass of the imaged cells [2]. Indeed, the refractive index is in direct relation with dry mass density and the total dry mass of an imaged object can be calculated by summing all the pixels values $\delta {\ell _{ij}}$ of the OPD image:

$$\delta m = {\rho _0}\mathrm{\int\!\!\!\int }\delta \ell ({x,y} )\textrm{d}x\textrm{d}y = {\rho _0}\,{p^2}\mathop \sum \limits_{\textrm{pixels}} \delta {\ell _{ij}}$$
where $\delta m$ is the dry mass of the entire cell, p the pixel size of the image at the sample plane, and ${\rho _0}$ is the dry mass volumetric density. For the sake of simplicity, we consider here that the background is uniform and set to $\delta \ell = 0$. In practice, when the background is non-uniform, image treatment may be required, as done in Ref. [3]. Equation (3) can be recast using the optical volume difference $\delta V$:
$$\delta m = {\rho _0}\delta V$$

An interesting physical quantity to deal with is the dry mass areal density $\delta \sigma $ (dry mass per unit area), directly proportional to the OPD image [4]:

$$\delta \sigma = {\rho _0}\delta \ell $$

Although the use of ${\rho _0}$ (a volumetric dry mass density) seems more natural in this context, the QPM community rather uses the so-called refractive index increment $\alpha = \rho _0^{ - 1}$, as originally introduced in 1954, leading to expressions such as $\delta m = {\alpha ^{ - 1}}\delta V$. The $\alpha $ value has been considered rather constant for biological media, ranging from $0.18$ to $0.21$ µm3·pg-1. We used in this study $\alpha = 0.18$ µm3·pg-1, i.e., ${\rho _0} = 5.5$ pg·µm-3 or g·ml-1.

Although the idea of measuring biomass using light dates from 1952 [5], one had to wait around 50 years for the development of CCD cameras and effective QPMs to observe the first optical measurements of dry mass at high resolution, independently by Popescu et al. [6] and Rappaz et al. [7] in 2008 and 2009. Since then, dry mass measurements using QPM have been reported by several groups for the study of mammalian cells [8,9,3,10,11,12], cancer cells [4,13,14,15], red blood cells [16], bacteria [17], yeast [18], and neurons [19,20]. More recently, masses of single molecules were reported using iScat by the group of Kukura, a technique called mass photometry [21].

Cross-grating wavefront microscopy is a QPM based on the use of a two-dimensional grating (aka a cross-grating) positioned around 1 mm from the sensor of a camera [22]. This simple association enables the mapping of the wavefront profile $\delta \ell $ of a light wave (Eq. (2)). When implemented on a microscope, this wavefront imaging technique can be used as a QPM technique, according to Eq. (2) that establishes the equivalence between phase and wavefront profiles [23,24]. Interestingly, because the dry mass directly depends on $\delta \ell $, CGM does not require the knowledge of the illumination wavelength to measure dry masses [3,12,20], unlike the other QPM techniques based on interferometry, which rather measure $\varphi $. Moreover, unlike Shack-Hartmann wavefront sensing, which has been very rarely used in bio-imaging because of its low spatial resolution [25], CGM is a wavefront imaging technique that can reach a diffraction-limited spatial resolution. In 2020, the group of Zangle reported the first observation of neural cells using CGM [20]. Using a home-made Matlab segmentation algorithm, they performed ensemble measurements of the dry masses of the somas and neurites separately, averaged over the field of view of a microscope. This work was the occasion to highlight the ability of QPM to follow the developments of cells in vitro for days, without perturbation or loss of the signal. However, this work only reported ensemble measurements.

In this article, we report the observation of neural cells in vitro using CGM, and present an image segmentation method for the precise estimation of the dry mass of individual neurites, from their birth until they connect. We highlight the importance of taking into account diffraction rings to achieve quantitative measurements, and the benefit of using CGM over other existing QPMs for the analysis of such tenuous biological objects. We evidence a precision of 2% on the dry mass measurements, i.e., around 0.08 pg. As a means to demonstrate the interest of conducting experiments at the single cell level, we present the observation of periodic oscillations of neurite outgrowth, in the pg range.

2. Material and methods

2.1 Cross-grating wavefront microscopy

All the experiments were conducted on a home-made microscope (Fig. 1), endowed with a 60${\times} $, 1.25 NA oil immersion objective lens (Olympus UPLFLN60XOI), a 200 mm tube lens (leading to a 66${\times} $ total magnification), a Köhler illumination using a mounted LED at 625 nm (Thorlabs M625L2-C1), and a heating stage (Tokai Hit, WSKMX) maintaining the cells at 37°C and regulating CO2 at 5%. For CGM measurements, we used a Zyla 5.5 sCMOS camera with a quadriwave lateral shearing interferometry (QLSI) cross-grating positioned at 0.86 mm from the camera. The pitch of the grating was $\Gamma = 39$ µm, 6 times the camera pixel size p. This configuration gives a zeta-factor [26] equal to $\zeta = \Gamma /2p = 3$. Another common value of the zeta-factor is 4, but higher values of $\zeta $ reduce the spatial resolution.

 figure: Fig. 1.

Fig. 1. (a) Schematic of the experimental setup. (b) Representation of the cross-grating used in CGM, characterized by a checkerboard pattern of $0 - \mathrm{\pi }$ phase shifts.

Download Full Size | PDF

To postprocess the raw images (called interferograms), we used a home-made algorithm that we detail in Ref. [24], and provide on a github repository [27]. Unlike more common algorithms, our algorithm does not yield a reduction of the number of pixels of the image by a factor of ${\zeta ^2}$. The number of pixels of the OPD image is the same as the number of pixels of the camera sensor ($2560 \times 2160$), because we perform a demodulation in the Fourier space, instead of an image crop by a factor of $\zeta $. This common ${\zeta ^2}$-reduction in the number of pixels enables faster computation, and thus faster live acquisition. Note that it does not yield a reduction of the resolution, only of the image definition, i.e. of the number of pixels (Fig. S4). However, avoiding this image reduction is a benefit for data visualization, especially on objects as tenuous and small as neurites. To provide a visual comparison of these low- and high-definition algorithms, Fig. 2 shows the same image of a neurite processed with both types of algorithms. Figure 2 is also the occasion to show that OPD images can be considered as dry mass density images as well, according to Eq. (5).

 figure: Fig. 2.

Fig. 2. (a) OPD image of a neurite processed with a full-pixel algorithm. (b) OPD image from the same measurement, processed with the more common pixel-reduction algorithm. The color bar shows the equivalence between OPD and dry mass density.

Download Full Size | PDF

2.2 Cell preparation

All procedures were in agreements with the guidelines established by the European Animal Care and Use Committee (86/609/CEE) and was approved by the local ethics committee (agreement G13O555). We used primary neural cell cultures from embryonic rat hippocampus. Pregnant female Wistar rats were obtained from Janvier labs. The plug date of the female rats was marked E0. All embryos staged at E18 from the female rats were used in the experiments. For cell culture and observation, we used 35-mm µ-Dishes, high Glass Bottom, 170 µm thick (Ibidi, 81158). µ-dishes were first incubated for 24 h with a PEI solution (1 ml, 20 µg/ml) and then rinsed 3 times with mQ water, leading to a PEI coating favouring cell adhesion. 150,000 to 200,000 cells were plated on the µ-dishes, with 1 ml of culture medium composed of Neurobasal/B27/L-glutamin with antifungal/antibiotics (Gibco 15240062), and let 1h30min in the incubator at 37°C, 5% CO2, before being placed on the microscope stage. The cells were gently washed before imaging to reduce the number of debris commonly found in primary cultures, and to ease the measurements of dry mass of neurites.

3. Results

3.1 Image processing method for single neurite dry-mass measurement

To measure the dry mass of a single neurite using Eq. (3), the OPD background value has to be precisely set. Indeed, CGM primarily measures wavefront gradients and the wavefront profile is obtained by integration of the gradients, to within a constant. For proper estimation of the dry mass, this constant has to be set such that the OPD of the background area is zero, defined as the part of the image where there is no object. However, this background may not be perfectly uniform, especially for long neurites covering a large part of the field of view, which makes precise dry mass estimation not as straightforward as what Eq. (3) suggests.

To circumvent this problem, the neurite to be measured was segmented, defining an area $\mathcal{S}$, and the background value $\delta {\ell _0}$ was calculated as the averaged OPD value over a few-pixel thick boundary (${\mathcal{S}}_0$) outlining the segmented area:

$$ \delta \ell_0=\frac{1}{{S}_{0}} \iint_{\mathcal{S}_{0}} \delta \ell(x, y) \mathrm{d} x \mathrm{d} y $$
$$\delta V = \mathop {\int\!\!\!\int }\limits_{\mathcal{S}} ({\delta \ell ({x,y} )- \delta {\ell_0}} )\textrm{d}x\textrm{d}y$$

The complexity of the problem comes down to automatically determine $\mathcal{S}$. There exist many image segmentation methods but for objects as tenuous as neurites, and surrounded by complex objects such a soma and debris, it was hard to find an automated algorithm that could approach human performance. Machine learning (ML) could be an effective method for this task, but ML algorithms need to be trained on a set of trustful data obtained by more common human-guided algorithms anyway.

We developed a semi-manual segmentation algorithm based on a magic-wand selection. The four-step procedure is described in Fig. 3, and the related Matlab code is provided in Suppl. Mater. and on github [28]. Here is a description of the four steps of the segmentation algorithm illustrated with a particular neuron image (Fig. 3(a)):

  • 1. Exclusion. The soma and any other unwanted objects at the vicinity of the neurite (like debris) have to be manually outlined to avoid their inclusion during the integration step (Fig. 3(b), red area). This step has to be performed mostly for the first image, because all these objects (soma and debris) usually do not move (or very little) during the acquisition. Outlining the soma is important to avoid the natural capture of its boundary by the magic wand, as represented in Fig. 3(d),(e). For most cells, the demarcation between the soma and the neurite is straightforward to define. However, in some cases, the demarcation with the soma is less abrupt, or neurites can arise from a lamellipodium (bottom cell of Fig. 4), making the onset of the neurite less easy to define. In such a case, the exclusion line can be drawn arbitrarily further along the neurite. Part of the dry mass may be missed, but any missing constant offset won't affect measurements of growth dynamics and velocity, which is what essentially matters.
  • 2. Magic-wand selection. The neurite of interest has to be clicked, ideally on the highest pixel value (Fig. 3(b), cross marker). Several clicks on several locations on the neurites may be necessary to properly capture the whole neurite.
  • 3. Adjustment. The magic-wand sensitivity has to be adjusted so that the whole neurite is properly captured, which defines the initial segmentation area $\mathcal{S}$ (Fig. 3(c)).
  • 4. Inflation. The initial area is dilated by a factor of f, ranging typically from 0.5 to 2 (Fig. 3(f),(g),h). The optical volume as a function of f, $\delta V(f )$, defined by Eq. (9), normally reaches a plateau from $f\mathrm{\ \mathbin{\lower.3ex\hbox{$\buildrel> \over {\smash{\scriptstyle\sim}\vphantom{_x}}$}}\ }1.3$, corresponding to the true value of the dry mass (Fig. 3(i)).
    $$\delta \ell_0(f)=\frac{1}{S_0(f)} \iint_{\mathcal{S}_0(\mathrm{f})} \delta \ell(x, y) \mathrm{d} x \mathrm{d} y$$
    $$ \delta V(f)=\iint_{\mathcal{S}(f)}\left(\delta \ell(x, y)-\delta \ell_0(f)\right) \textrm{d}x\textrm{d}y $$
    This inflation of the area above the apparent outline of the neurites actually captures the diffraction rings that we barely see in the image and that would not be captured by standard segmentation algorithms. However, they notably contribute in the total dry mass. To ascertain that diffraction are actually responsible for this missing dry mass, numerical simulations have been conducted using IF-DDA [29] and shown in SI (Fig. S6).
  • 5. Plateau selection. To obtain more accurate measurement, the retained value of $\delta V$ does not only correspond to one selected value on the plateau, but to the average value over the plateau. For this purpose, two f values are selected, ${f_\textrm{1}}$ and ${f_\textrm{2}}$, such that the dry mass equals.
    $$\delta m = \frac{{{\rho _0}}}{{{f_2} - {f_1}}}\int \delta V(f )\textrm{d}f$$
    Above a certain value of f, around $f = 2$ in Fig. 3(i), the estimated value of $\delta m$ becomes unstable due to debris and even tends to diverge due to image noise. One neurite segmentation (steps 1 to 5) takes around 30 s.

 figure: Fig. 3.

Fig. 3. (a) OPD image of a neuron and its unique neurite. (b) Same OPD image as (a). The red area corresponds to the manual exclusion of the soma. The cross indicates where the magic wand was subsequently applied. (c) Area $\mathcal{S}$ (orange area) captured by the magic wand after adjustment of the magic-wand sensitivity. (d) Representation of a magic-wand application, without soma removal, and (e) corresponding area $\mathcal{S}$ captured by the same magic-wand application. A problem of capture of the soma boundary is evidenced. (f, g, h) OPD images of the same neurite along with segmented areas ${\mathcal{S}}(f )$ for 3 values of the inflation factor f, namely $0.6$, $1$, and $1.35$. (i) Dry mass of the neurite measured as a function of the inflation factor f, featuring a plateau above $f = 1.3$. Two dry mass values are highlighted: $\delta m = 2.44$ pg, corresponding to the properly estimated value, and $\delta m = 1.79$ pg, an underestimated value corresponding to $f = 1$, highlighting the importance of inflating the initial segmented area.

Download Full Size | PDF

 figure: Fig. 4.

Fig. 4. (a) OPD image of neural cells, corresponding to time $t = 14.5$ h. The 5 neurites are indicated and numbered. (b – f) Dry mass and length as a function of time for the 5 neurites indicated in (a).

Download Full Size | PDF

This algorithm has been tested and intensively used by 3 different users. Figure S3 plots dry mass measurements performed by two different users on the exact same series of data, demonstrating a very good reproducibility.

The interest of measuring the OV on a plateau of $\delta V(f )$, rather than at $f = 1$, is not only to obtain a value closer to reality. Note that the true dry mass is anyway biased by the arbitrary choice of the $\mathrm{\alpha }$ coefficient within the range $0.18$ to $0.21$ µm3·pg-1. The other, and more important, interest of avoiding $f = 1$ is rather that the slope of $\delta V(f )$ at $f = 1$ is abrupt, leading to a much larger variation of the measured OV, and thus of the dry mass. To illustrate this statement, Figure S2c plots the measured dry masses over 30 successive frames of a neurite of fixed neurons with the two modalities: at $f = 1$ and on the plateau. While the standard deviation of the measurement is 1.3% on the plateau, it is 11% at $f = 1$.

In Ref. [30], The authors introduce a segmentation algorithm to measure dry masses of red blood cells from QPM images. To quantify the efficiency of the measurements, they introduced a metrics, named the segmentation accuracy, defined as “the absolute value of correlation between segmented RBCs image and reference image, which is manually obtained”. We cannot apply such a metrics because the segmentations we perform are already obtained semi-manually, i.e., they are already visually checked individually. However, our algorithm could be the basis for the generation of a ground truth for training a deep learning algorithm for neurite segmentation. In this case, such a metrics could be advantageously used.

3.2 Illustration of the numerical procedure on experimental CGM images

To illustrate the efficiency and the degree of accuracy of the algorithm, Fig. 4 presents dry mass measurements of neurites stemming from a 15-h long image acquisition, one image per minute, of two nearby neurons. Figure 4(a) displays the OPD image of the field of view at the end of the data acquisition, where all the neurites of interest are indicated. The full movie is provided in Suppl. Mater (Visualization 1 and Fig. S7 for a time lapse of 16 images). Of particular interest are neurites 1 and 5, which eventually connect, just at the end of the acquisition. Figure 4(b)–(f) plots both the dry mass and the length of the neurites.

All the neurites are a few pg in mass. Rare are the dry mass measurements reaching such small values. Common dry mass measurements using QPM rather report values around a few 100s of pg (mass of eukaryotic cells). The precision of neurite dry mass measurements was 0.08 pg, calculated as the standard deviation of series of dry mass measurements for which the dry mass of the neurite was constant (see Fig. S1). The precision of other QPMs is rather in the 1 pg range [31,32]. CGM recently confirmed its ability to yield highly precise OPD measurements, better than other conventional QPM technique such as DHM [33]. We recently dedicated an article on the precision and accuracy of CGM, and demonstrated that the OPD precision can typically reach 0.1 nm, i.e. a fraction of hydrogen atom, corresponding to a precision in $\delta \sigma $ of 0.55 fg·µm-1 (Eq. (5)) [26]. This high precision comes from the fact that CGM does not require the use of a reference arm, like other QPMs based on interferometry measurements for which it makes the setup highly sensitive to environmental perturbation, such as air flow, mechanical drift, thermal variation, etc. The common-path configuration also ensures a perfectly flat background. No image treatment has been performed in our data, except the correction of a tilt on the OPD images. With CGM, all the optical metrology happens within a distance of around 1 mm, between the camera sensor and the cross-grating, a volume that can be made particularly solid, robust and isolated from any perturbation. Thanks to this high precision, we recently could characterize for the first time the optical properties of small objects and 1-atom-thick 2D-materials using QPM [34,35].

In the case of the data presented in Fig. 4, this high precision, along with the magic-wand algorithm, enabled us to observe fine variations of neurite masses, in the sub-pg range. Most of the time, the evolution looks chaotic, with succeeding periods of activity, inactivity, regression, growth. Interestingly, neurite lengths do not necessarily correlate with their dry mass. In many occasions, neurites thicken without elongation, and reciprocally, change length without mass increase.

Of particular interest with this set of data is neurite 1, which featured a periodic oscillation of its dry mass, a behavior that has never been reported in the literature. This behavior is not systematic, only this neurite shows this effect in this set of data. During these oscillations, the neurite length also demonstrates an oscillation, correlated with the dry mass, but much weaker. To determine the oscillation frequency, we conducted a short-time Fourier transform analysis: from the dry mass profile $\delta m(t )$ of neurite 1, the baseline $\delta {m_0}(t )$ is estimated with a low-pass filter and then subtracted (green line in Fig. 5(a)). This operation is aimed to remove the low spatial frequencies of the signal that would have masked the signal of interest in Fig. 5(b). Then, a spectrogram (modulus of the short-time Fourier transform with a Gaussian window of size 4 h) of this baseline-free signal is computed to provide a time-frequency representation (Fig. 5(b)). Finally, the ridge of the time-frequency pattern is calculated (red dotted line if Fig. 5(b)) giving an estimation of the instantaneous frequency of the neurite oscillation as a function of time. We can see that the oscillation frequencies lie around $1$ to $1.5$ h-1 at the beginning of the acquisition, and then splits, or maybe becomes more chaotic from $t \approx 10$ h. Let us notice that this time corresponds to neurite 5 arriving is close proximity to neurite 1 (see Visualization 1 in Suppl. Info.). We show another event of neurite oscillation in Suppl. Mater (Fig. S5), with the same period of around 1 h.

 figure: Fig. 5.

Fig. 5. (a) Dry mass of neurite 1 as a function of time $\delta m(t )$, sampled at $30$ h-1 (green line) and its baseline $\delta {m_0}(t )$ (black line). Subtraction of these two profiles (yellow line). (b) Spectrogram (estimated instantaneous frequency of $\delta m$ as a function of time) calculated with a time window of 4 h. The dotted red line is the calculated ridge of the image, giving the estimated frequency of the oscillation as a function of time.

Download Full Size | PDF

To date, we did not elucidate the origin of this oscillation, although we believe it may be correlated with the dynamics of actin waves. Because the aim of this article is primarily to demonstrate the possibility of CGM to investigate single neurites with high precision, and give the tool to do it, we will leave this question for a future work.

4. Discussion and summary

We show in this article that cross-grating wavefront microscopy (CGM) can be used to monitor the dynamics of single neurites in vitro, via measurements of dry mass with a sub-pg sensitivity. For this purpose, we developed a semi-manual segmentation algorithm, suited for neurite segmentation. Previously reported CGM measurements on neurons only investigated ensemble measurements of soma and neurite dry masses, an approach that naturally erases any particular behaviour at the single cell level. We show that imaging such tenuous objects implies two requirements: First, the neurite segmentation has to take into account the diffraction rings that are barely visible on the OPD images in order to achieve precise measurements. Second, a highly sensitive QPM has to be used, the neurite OPD being only around 10 to 20 nm. CGM performs very well, with its typical noise amplitude of 0.1 nm, and because it yields images with a particular flat background over the field of view, being insensitive to environmental perturbation, unlike common QPM technique based on interferometry. CGM is also simple (no modification of the microscope), cost-effective (only involves the addition of a cross-grating). Moreover, CGM offers easy fluorescence/phase correlation measurements because it can be easily implemented on any standard fluorescence microscope. For all these reasons, CGM has the potential for becoming a major QPM technique in the near future, and in particular in the field of neurobiology.

Funding

Association Nationale de la Recherche et de la Technologie.

Acknowledgments

The authors thank Christophe Leterrier and Fanny Boroni-Rueda for providing the primary neural cells, and for helpful discussion; Frédéric Galland and Maëlle Bénéfice for helpful discussion; and the ANRT for financial support.

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Supplemental document

See Supplement 1 for supporting content.

References

1. Y. Park, C. Depeursinge, and G. Popescu, “Quantitative phase imaging in biomedicine,” Nat. Photonics 12(10), 578–589 (2018). [CrossRef]  

2. T. A. Zangle and M. A. Teitell, “Live-cell mass profiling: an emerging approach in quantitative biophysics,” Nat. Methods 11(12), 1221–1228 (2014). [CrossRef]  

3. S. Aknoun, J. Savatier, P. Bon, F. Galland, L. Abdeladim, B. Wattellier, and S. Monneret, “Living cell dry mass measurement using quantitative phase imaging with quadriwave lateral shearing interferometry: an accuracy and sensitivity discussion,” J. Biomed. Opt. 20(12), 126009 (2015). [CrossRef]  

4. P. Girshovitz and N. T. Shaked, “Generalized cell morphological parameters based on interferometric phase microscopy and their application to cell life cycle characterization,” Biomed. Opt. Express 3(8), 1757–1773 (2012). [CrossRef]  

5. R. Barer, “Interference microscopy and mass determination,” Nature 169(4296), 366–367 (1952). [CrossRef]  

6. G. Popescu, Y. Park, N. Lue, C. Best-Popescu, L. Deflores, R. R. Dasari, M. S. Feld, and K. Badizadegan, “Optical imaging of cell mass and growth dynamics,” American Journal of Physiology-Cell Physiology 295(2), C538–C544 (2008). [CrossRef]  

7. B. Rappaz, E. Cano, T. Colomb, J. Kühn, C. Depeursinge, V. Simanis, P. J. Magistretti, and P. Marquet, “Noninvasive characterization of the fission yeast cell cycle by monitoring dry mass with digital holographic microscopy,” J. Biomed. Opt. 14(3), 034049 (2009). [CrossRef]  

8. K. L. Cooper, S. Oh, Y. Sung, R. R. Dasari, M. W. Kirschner, and C. J. Tabin, “Multiple phases of chondrocyte enlargement underlie differences in skeletal proportions,” Nature 495(7441), 375–378 (2013). [CrossRef]  

9. Y. Sung, A. Tzur, S. Oh, W. Choi, V. Li, R. R. Dasari, Z. Yaqoob, and M. W. Kirschner, “Size Homeostasis in Adherent Cells Studied by Synthetic Phase Microscopy,” Proc. Natl. Acad. Sci. 110(41), 16687–16692 (2013). [CrossRef]  

10. M. E. Kandel, W. Lu, J. Liang, O. Aydin, T. A. Saif, and G. Popescu, “Cell-to-Cell Influence on Growth in Large Populations,” Biomed. Opt. Express 10(9), 4664 (2019). [CrossRef]  

11. P. A. Sandoz, C. Tremblay, F. G. van der Goot, and M. Frechin, “Image-Based Analysis of Living Mammalian Cells Using Label-Free 3D Refractive Index Maps Reveals New Organelle Dynamics and Dry Mass Flux,” PLoS Biol 17(12), e3000553 (2019). [CrossRef]  

12. S. Aknoun, M. Yonnet, Z. Djabari, F. Graslin, M. Taylor, T. Pourcher, B. Wattellier, and P. Pognonec, “Quantitative Phase Microscopy for Non-Invasive Live Cell Population Monitoring,” Sci. Rep. 11(1), 4409 (2021). [CrossRef]  

13. O. Tolde, A. Gandalovičová, A. Křížová, P. Veselý, R. Chmelík, D. Rosel, and J. Brábek, “Quantitative Phase Imaging Unravels New Insight into Dynamics of Mesenchymal and Amoeboid Cancer Cell Invasion,” Sci. Rep. 8(1), 12020 (2018). [CrossRef]  

14. V. Ayyappan, A. Chang, C. Zhang, S. K. Paidi, R. Bordett, T. Liang, I. Barman, and R. Pandey, “Identification and Staging of B-Cell Acute Lymphoblastic Leukemia Using Quantitative Phase Imaging and Machine Learning,” ACS Sens. 5(10), 3281–3289 (2020). [CrossRef]  

15. Y. Su, R. Fu, W. Du, H. Yang, L. Ma, X. Luo, R. Wang, X. Lin, X. Jin, X. Shan, W. Lv, and G. Huang, “Label-Free and Quantitative Dry Mass Monitoring for Single Cells during In Situ Culture,” Cells 10(7), 1635 (2021). [CrossRef]  

16. K. G. Phillips, S. L. Jacques, and O. J. T. McCarty, “Measurement of Single Cell Refractive Index, Dry Mass, Volume, and Density Using a Transillumination Microscope,” Phys. Rev. Lett. 109(11), 118105 (2012). [CrossRef]  

17. M. Mir, Z. Wang, Z. Shen, M. Bednarz, R. Bashir, I. Golding, S. G. Prasanth, and G. Popescu, “Optical Measurement of Cycle-Dependent Cell Growth,” Proc. Natl. Acad. Sci. U.S.A. 108(32), 13124–13129 (2011). [CrossRef]  

18. D. Midtvedt, E. Olsén, F. Höök, and G. D. M. Jeffries, “Label-Free Spatio-Temporal Monitoring of Cytosolic Mass, Osmolarity, and Volume in Living Cells,” Nat. Commun. 10(1), 340 (2019). [CrossRef]  

19. P. Cintora, J. Arikkath, M. Kandel, G. Popescu, and C. Best-Popescu, “Cell Density Modulates Intracellular Mass Transport in Neural Networks,” Cytometry 91(5), 503–509 (2017). [CrossRef]  

20. S. Pradeep, T. Tasnim, H. Zhang, and T. A. Zangle, “Simultaneous Measurement of Neurite and Neural Body Mass Accumulation via Quantitative Phase Imaging,” Analyst 146(4), 1361–1368 (2021). [CrossRef]  

21. G. Young, N. Hundt, D. Cole, A. Fineberg, J. Andrecka, A. Tyler, A. Olerinyova, A. Ansari, E. G. Marklund, M. P. Collier, S. A. Chandler, O. Tkachenko, J. Allen, M. Crispin, N. Billington, Y. Takagi, J. R. Sellers, C. Eichmann, P. Selenko, L. Frey, R. Riek, M. R. Galpin, W. B. Struwe, J. L. P. Benesch, and P. Kukura, “Quantitative Mass Imaging of Single Molecules,” Science 360(6387), 423–427 (2018). [CrossRef]  

22. J. Primot and N. Guerineau, “Extended Hartmann Test Based on the Pseudoguiding Property of a Hartmann Mask Completed by a Phase Chessboard,” Appl. Opt. 39(31), 5715–5720 (2000). [CrossRef]  

23. P. Bon, G. Maucort, B. Wattellier, and S. Monneret, “Quadriwave Lateral Shearing Interferometry for Quantitative Phase Microscopy of Living Cells,” Opt. Express 17(15), 13080–13094 (2009). [CrossRef]  

24. G. Baffou, “Quantitative Phase Microscopy Using Quadriwave Lateral Shearing Interferometry (QLSI): Principle, Terminology, Algorithm and Grating Shadow Description,” J. Phys. D: Appl. Phys. 54(29), 294002 (2021). [CrossRef]  

25. H. Gong, T. E. Agbana, P. Pozzi, O. Soloviev, M. Verhaegen, and G. Vdovin, “Optical Path Difference Microscopy with a Shack–Hartmann Wavefront Sensor,” Opt. Lett. 42(11), 2122 (2017). [CrossRef]  

26. B. Marthy and G. Baffou, “Cross-Grating Phase Microscopy (CGM): In Silico Experiment (Insilex) Algorithm, Noise and Accuracy,” Opt. Commun. 521, 128577 (2022). [CrossRef]  

27. G., Baffou, CGM Process, Github, 2022, Github.Com/Baffou/CGMprocess.

28. G. Baffou, CGM_magicWandSegmentation, Github, 2022, Github.Com/Baffou/CGM_magicWandSegmentation.

29. P. C. Chaumet, D. Sentenac, G. Maire, M. Rasedujjaman, T. Zhang, and A. Sentenac, “IFDDA, an Easy-to-Use Code for Simulating the Field Scattered by 3D Inhomogeneous Objects in a Stratified Medium: Tutorial,” J. Opt. Soc. Am. A 38(12), 1841 (2021). [CrossRef]  

30. F. Yi, I. Moon, B. Javidi, D. Boss, and P. Marquet, “Automated Segmentation of Multiple Red Blood Cells with Digital Holographic Microscopy,” J. Biomed. Opt. 18(2), 026006 (2013). [CrossRef]  

31. K. Park, L. J. Millet, N. Kim, H. Li, X. Jin, G. Popescu, N. R. Aluru, K. J. Hsia, and R. Bashir, “Measurement of Adherent Cell Mass and Growth,” Proc. Natl. Acad. Sci. 107(48), 20691–20696 (2010). [CrossRef]  

32. G. Popescu, K. Park, M. Mir, and R. Bashir, “New Technologies for Measuring Single Cell Mass,” Lab Chip 14(4), 646–652 (2014). [CrossRef]  

33. C. Allier, L. Hervé, O. Mandula, P. Blandin, Y. Usson, J. Savatier, S. Monneret, and S. Morales, “Quantitative Phase Imaging of Adherent Mammalian Cells: A Comparative Study,” Biomed. Opt. Express 10(6), 2768 (2019). [CrossRef]  

34. S. Khadir, P. Bon, D. Vignaud, E. Galopin, N. McEvoy, D. McCloskey, S. Monneret, and G. Baffou, “Optical Imaging and Characterization of Graphene and Other 2D Materials Using Quantitative Phase Microscopy,” ACS Photonics 4(12), 3130–3139 (2017). [CrossRef]  

35. S. Khadir, D. Andrén, P. C. Chaumet, S. Monneret, N. Bonod, M. Käll, A. Sentenac, and G. Baffou, “Full optical characterization of single nanoparticles using quantitative phase imaging,” Optica 7(3), 243 (2020). [CrossRef]  

Supplementary Material (2)

NameDescription
Supplement 1       Supplementary Material file gathering all the supplementary figures and comments.
Visualization 1       Movie related to Fig. 4

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (5)

Fig. 1.
Fig. 1. (a) Schematic of the experimental setup. (b) Representation of the cross-grating used in CGM, characterized by a checkerboard pattern of $0 - \mathrm{\pi }$ phase shifts.
Fig. 2.
Fig. 2. (a) OPD image of a neurite processed with a full-pixel algorithm. (b) OPD image from the same measurement, processed with the more common pixel-reduction algorithm. The color bar shows the equivalence between OPD and dry mass density.
Fig. 3.
Fig. 3. (a) OPD image of a neuron and its unique neurite. (b) Same OPD image as (a). The red area corresponds to the manual exclusion of the soma. The cross indicates where the magic wand was subsequently applied. (c) Area $\mathcal{S}$ (orange area) captured by the magic wand after adjustment of the magic-wand sensitivity. (d) Representation of a magic-wand application, without soma removal, and (e) corresponding area $\mathcal{S}$ captured by the same magic-wand application. A problem of capture of the soma boundary is evidenced. (f, g, h) OPD images of the same neurite along with segmented areas ${\mathcal{S}}(f )$ for 3 values of the inflation factor f, namely $0.6$, $1$, and $1.35$. (i) Dry mass of the neurite measured as a function of the inflation factor f, featuring a plateau above $f = 1.3$. Two dry mass values are highlighted: $\delta m = 2.44$ pg, corresponding to the properly estimated value, and $\delta m = 1.79$ pg, an underestimated value corresponding to $f = 1$, highlighting the importance of inflating the initial segmented area.
Fig. 4.
Fig. 4. (a) OPD image of neural cells, corresponding to time $t = 14.5$ h. The 5 neurites are indicated and numbered. (b – f) Dry mass and length as a function of time for the 5 neurites indicated in (a).
Fig. 5.
Fig. 5. (a) Dry mass of neurite 1 as a function of time $\delta m(t )$, sampled at $30$ h-1 (green line) and its baseline $\delta {m_0}(t )$ (black line). Subtraction of these two profiles (yellow line). (b) Spectrogram (estimated instantaneous frequency of $\delta m$ as a function of time) calculated with a time window of 4 h. The dotted red line is the calculated ridge of the image, giving the estimated frequency of the oscillation as a function of time.

Equations (10)

Equations on this page are rendered with MathJax. Learn more.

φ ( x , y ) = 2 π λ ( n n 0 ) h ( x , y )
φ ( x , y ) = 2 π λ δ ( x , y )
δ m = ρ 0 δ ( x , y ) d x d y = ρ 0 p 2 pixels δ i j
δ m = ρ 0 δ V
δ σ = ρ 0 δ
δ 0 = 1 S 0 S 0 δ ( x , y ) d x d y
δ V = S ( δ ( x , y ) δ 0 ) d x d y
δ 0 ( f ) = 1 S 0 ( f ) S 0 ( f ) δ ( x , y ) d x d y
δ V ( f ) = S ( f ) ( δ ( x , y ) δ 0 ( f ) ) d x d y
δ m = ρ 0 f 2 f 1 δ V ( f ) d f
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.