ORIGINAL RESEARCH article
Measuring Absolute Cell Volume Using Quantitative-Phase Digital Holographic Microscopy and a Low-Cost, Open-Source, and 3D-Printed Flow Chamber
- 1Département de Physique, de Génie Physique et d'Optique, Université Laval, Quebec, QC, Canada
- 2Centre de Recherche CERVO, Université Laval, Quebec, QC, Canada
- 3Centre d'Optique, Photonique et Laser, Université Laval, Quebec, QC, Canada
- 4Département de Psychiatrie et Neurosciences, Université Laval, Quebec, QC, Canada
Although cell volume is a critical cell parameter and acute changes of its regulation mechanisms can reveal pathophysiological states, its accurate measurement remains difficult. On the other hand, custom-made devices are now possible through the application of 3D printing that meet the needs of certain specific applications. Therefore, in this article we describe the development of a low-cost, open-source, and 3D-printed flow chamber optimized to perform non-invasive accurate measurements of important cell biophysical parameters. This measurement development includes the absolute cell volume, mean cell thickness, and whole-cell refractive index using quantitative-phase digital holographic microscopy. Specifically, this advancement makes it possible to develop a flow chamber that preserves cell viability while exhibiting a fast washout process for the two-liquid decoupling procedure—an experimental procedure allowing to separately measure the intracellular refractive index and cell thickness from the quantitative-phase signal—characterized by a homogeneous return of the extracellular refractive index over the entire imaging slit toward the known value of the last perfusion solution. Such a fast washout process has been extensively characterized by monitoring the time required to obtain a stabilization of the quantitative-phase signal following the perfusion of the two liquids. A more rapid quantitative-phase signal stabilization time results in fewer cell changes, including important cell position and shape modifications. These changes likely take place over the washout process time and ultimately result in important artifacts to all cell biophysical measurements. The proposed flow chamber has been used to perform measurements of the whole-cell refractive index, mean cell thickness, and absolute cell volume of three cell types in culture during resting state, and hypo-osmotic and hyper-osmotic challenges. For a typical cell body, these biophysical parameters are measured with a precision of 0.0006, 100 nm, and 50 μm3, respectively. Finally, the design files of the flow chamber will be shared with the scientific community through an open-source model.
Animal cell membranes are highly permeable to water. Movements of water across the membranes are in large part dictated by osmotic pressure gradients. Consequently, even at constant extracellular osmolarity, volume constancy of any mammalian cell is permanently challenged by its own activity including, for example, movements of different ionic species across the cell membrane or transport of osmotically active substances involved in cell metabolism all of which induce osmotic water movements [1, 2]. On the other hand, a variety of vital cell functions, including enzyme activity, gene expression mechanisms, second messenger cascades, and hormone releases are all extremely dependent on a homeostatic intracellular ionic environment [3, 4]. Therefore, cell volume is one of the critical cell parameters and is closely regulated. Accordingly, cells have developed numerous efficient volume regulatory mechanisms involving multiple ion transport processes including channels, cotransporters, and exchangers which operate to maintain a constant ionic balance and hydrostatic pressure gradient. To do so, osmotically active solutes (mainly inorganic ions such as Na+, K+, Cl−, or organic osmolytes) are gained or lost to regulate the cell volume .
Currently, reliable and non-invasive measurements of the absolute cell volume Vc are difficult to achieve . Although current and common flow cytometry techniques offer tremendous high-throughput capability and can be seen as the gold standard, they are dedicated to the measurement of Vc using large samples of isolated cells in suspension . Nevertheless, studying volume homeostasis and its regulatory mechanisms requires the ability to be able to measure Vc within the environmental conditions in which cells grow, migrate, and communicate with one another. In this case, there is no methodological consensus to be able to achieve such Vc measurements. Among the few suitable techniques, dye exclusion microscopy, either based on absorption  or fluorescence , use the equivalent displacement of stained liquid created by a cell to measure its thickness, but this requires extracellular labels. This thickness measurement is used in conjunction with the projected cell surface (in this instance, measured in the corresponding microscopy image) to calculate Vc.
Alternatively, quantitative-phase imaging (QPI) consists of an ensemble of microscopy techniques offering the possibility of performing label-free, non-invasive, high-resolution, and real-time imaging of transparent microscopic specimens such as living cells in culture [10, 11]. QPI measures the minute phase shift, namely the quantitative-phase signal (QPS), that transparent specimens induce on the diffracted wave-front, as long as the observed sample is a weakly-scattering object, i.e., differing from the surrounding medium only by a slight difference of refractive index [12–17]. The QPS, proportional to the optical path length (OPL), thus depends on both the intracellular refractive index (IRI) and the cell thickness. According to the seminal work of Barer , this IRI dependency has allowed the development of QPI as an efficient technique to quantitatively monitor cellular dry mass dynamics, particularly in connection with the cell growth cycle [19–22]. Different QPI approaches based on interferometric principles have also succeeded to detect cell membrane fluctuations at the nanometric scale, especially with respect to red blood cells [23–25].
Despite these fruitful applications in simple cases, the information on the intrinsic mix of cell content and morphology, the QPS, and its variations, remain generally difficult to interpret in cell biology . Measuring the IRI and the cell thickness separately represents thus the key scientific step leading to the calculation of various quantitative cellular biophysical parameters, including Vc [27, 28], the osmotic membrane water permeability , the refractive index of transmembrane water and solute flux , allowing the study of specific physiological and pathophysiological cellular processes [30–33]. Specifically, Vc is calculated from both the cell thickness and the projected cell surface, as in the case of dye exclusion microscopy.
Optical diffraction tomography makes it possible to obtain, also from the QPS, both the refractive index distribution and the morphology of living cells in 3D [34, 35]. However, this unique 3D capability comes at the price of complex imaging setups and relies heavily on the approximations used in the algorithms of the reconstruction process which requires considerable expertise to be used.
However, less complicated experimental procedures have been developed to extract both the IRI and the corresponding cell thickness from the QPS [26, 36–42], which greatly promotes their application. Conceptually, there are two main classes of such procedures. The first involves the extraction of the IRI by cell thickness evaluation either by shape approximation using geometrical models, direct measurement with complementary imaging modalities or physical confinement in the axial direction using microstructures with a priori known geometry. While the shape approximation and the physical confinement are more suited for suspended rather than adherent cells, the direct measurement with a complementary imaging modality is usually conducted separately, making the observation of dynamic biological processes cumbersome. The second group of procedures solves the IRI and the cell thickness coupling problem by achieving two different QPS measurements resulting in a system of independent linear equations with two unknowns having a unique solution. Among these techniques is the two-liquid decoupling approach which is particularly attractive for adherent cells as it does not require any cell morphology assumptions, nor axial constraining, generally resulting in better physiological environmental conditions than in those in the first class of procedures. Used in conjunction with quantitative-phase digital holographic microscopy (QP-DHM), it provides accurate measurements of the average refractive index and the mean thickness of a cellular body . However, this technique has several drawbacks related to its washout process arising from the requirement to sequentially perfuse two different equi-osmolar physiological solutions with different, but a priori, known refractive indices.
Recent developments in affordable and precise 3D printers, as well as the availability of computer aided-design (CAD) software through additive manufacturing, offer a great opportunity to create custom-made and innovative devices that meet the needs of a specific application. For researchers, 3D printing is a versatile and inexpensive technique that allows rapid prototyping, and appreciable improvements of previous versions and variations of a design . Within the context of live-cell imaging, 3D printing is therefore ideal to create microscope accessories tailored to fit the precise needs of QP-DHM, and with QPI in a more general manner.
To our knowledge, the authors have created the first, low-cost, open-source, and 3D-printed flow chamber with tailor-made flow field characteristics that produces an efficient washout process. Within the framework of the two-liquid decoupling approach, such a washout process is characterized by an early stage during which the QPS is temporally unstable, followed by a period of QPS stabilization corresponding to a state where the refractive index of the second solution is distributed in a spatially uniform way, throughout the flow chamber. From a fluid mechanics point of view, the use of a steady-state and a fully-developed laminar flow (which reduces lateral mixing) may be beneficial in this situation.
The design and scientific particularities of this 3D-printed flow chamber will be presented in this paper. First, the mechanical design of the do-it-yourself flow chamber is presented. Second, the optimization of the geometry of the flow chamber to obtain such a flow of fluid (resulting from numerical simulations), will be explained. In addition, the presence of a steady-state and fully-developed laminar flow has been confirmed experimentally for the geometry of the flow chamber. In such a situation, the washout time required to alternate between the two liquids has been extensively experimentally characterized. The rapidity at which the QPS stabilizes (an indication of the homogeneity of the refractive index of the perfusion medium) has been assessed. As a validation to demonstrate how this proposed 3D-printed flow chamber works, it has been used to perform decoupling experiments during the resting state of three common cell types in culture, namely primary neurons, HEK293T, and NIH3T3 cells. Specifically, both and as well as their corresponding Vc, have been systematically calculated for each cell type. Furthermore, the cell viability rate inside the flow chamber has been compared to a control condition to determine its compatibility with live cell imaging. Finally, the changes of these biophysical parameters have been quantitatively evaluated during hypo-osmotic and hyper-osmotic challenges.
Theory and Principles
The cell-phase signal of each (x, y)-pixel of a quantitative-phase image of QP-DHM can be expressed as:
where λ0 is the free-space wavelength of the light source, z is the axial coordinate, hc(x, y) is the cell thickness corresponding to the OPL of each pixel, nc(x, y, z) is the 3D distribution of the IRI, is the mean value of the IRI along hc(x, y) (otherwise known as the integral IRI), nm is the refractive index of the perfusion medium, and H is the thickness of the flow chamber. Similarly, the phase signal corresponding to the background of a quantitative-phase image can be formulated as:
Monitoring and subsequently reusing this background value are important because this compensates for mechanical and thermal instabilities of the setup during the experiment. In addition, by considering limited background regions around the cell under study, effective variations in the chamber height are not an issue . By subtracting Equation (3) from Equation (2), the QPS of each (x, y)-pixel  can be expressed as:
Although temporal dependency is not explicitly indicated in Equation (4), it remains valid for each quantitative-phase image of a time-lapse. The dual dependence expressed by Equation (4) (known as the refractive index-thickness coupling problem ) results in a counterintuitive and challenging interpretation of the QPS.
In order to obtain these two contributions independently, a two-liquid decoupling procedure  was performed. In this method, two iso-osmotic and iso-tonic perfusion media having two different but a priori known refractive indices (namely nm1 and nm2) were sequentially perfused. The extracellular refractive index of the perfusion medium varied from nm1 to nm2 with a transient period during which the QPS is not defined, followed by a stabilization period. Meanwhile, a time-lapse of quantitative-phase images was recorded. This two-liquid decoupling approach created a system of two equations with two unknowns that can be solved as:
where the subscripts 1 and 2 refer to two different QPS stabilized states corresponding to the two different perfusion media . From there, the calculation of and can be performed by averaging, for each state i = 1, 2, Δϕci(x, y) over all the segmented pixels contained in a single-cell body, namely . Therefore, we obtain:
From the automated segmentation, the projected surface of the cell, Sc, is obtained which makes it possible to calculate the absolute cell volume , Vc, as follows:
where Ap is the area of a magnified pixel, and Nc is the number of pixels corresponding to the cell body. However, the following working assumptions must be fulfilled to obtain accurate results. It is assumed that and do not exhibit significant variations during the delay caused by the washout process. This notably requires the use of two equi-osmolar perfusion media to avoid any osmotic challenges. Furthermore, in order to minimize the effect of cell changes that take place over time and affecting either or (or both), the washout process allowing a return of the extracellular refractive indices over the entire imaging slit, toward the known values of either nm1 or nm2, must be as short as possible.
Materials and Methods
Fabrication of the Flow Chamber
The mechanical design of the flow chamber was drawn with SolidWorks (2016, Dassault Systèmes). The flow chamber was 3D-printed with high-strength PETG plastic (Guideline, Taulman3D) using the Ultimaker 2 Extended+ printer. The base and the lid were printed with a 0.4 mm nozzle (Micro Swiss) with light infill. Both the base and the lid each have four circular holes meant for inserting the rare-earth magnets (5862K961, McMaster-Carr). These magnets are glued with fast-adhesive glue (5551T72, McMaster-Carr). The parafilm gasket contained in the base was cut with scissors from a sheet of parafilm (M, BEMIS). The sample holder was printed with a 0.25 mm nozzle (Micro Swiss) with 100 and 30% infill overlap. The 30% infill overlap of the central element generally prevents the internal leaks; however, it must be treated with acetone to ensure complete waterproofness considering the pressure level while the fluid is flowing . Only the top (the roughest surface) of the sample holder shall be exposed to the acetone (270725, Sigma-Aldrich) vapors for the advised 15 min. Following this treatment, a clean 18 mm coverslip is pressed in the groove made for this purpose, resulting in its sealing and fixing. Pressure was applied and held on the coverslip for about 1 min, with the first 10 s being decisive. On each side of the central element, stainless steel microtubes (304F10125X006SL, MicroGroup) are partially inserted and glued with fast adhesive. Cutting the tubes is a delicate process since an inappropriate technique results in it collapsing. A sharp blade (LB, OLFA) is used to roll the microtube while applying some pressure. After several back and forth cycles, the microtubes can be cut without altering their cavity. A detailed step-by-step protocol for the fabrication of the flow chamber can be found in Appendix A.
Mounting the Flow Chamber
The chamber is intended to reduce the handling time and ensure the well-being of the studied cells in keeping with the reality of cellular imaging. First, the parafilm gasket is applied to the base. Then, a small amount of perfusion liquid is injected with a connected syringe in the upside-down sample holder which has a fixed-top window. The 18 mm coverslip containing the cell culture (the bottom coverglass), is next placed upright directly on the bottom mount and wetted with a drop of temperate medium. Then, the filled central element is rotated back and positioned on the base. The chamber is then closed by pressure from only the lid and the magnets. Finally, any air bubbles are quickly purged by gently injecting medium via the syringe fixed to the sample holder.
Computational Fluid Dynamics
Simulations were made with the computational fluid dynamics (CFD) module COMSOL Multiphysics (5.3a, COMSOL) and steady-state solutions were sought. A top-view cross-section was used to define the geometry for the simulation and assuming a 2D representation to be a satisfactory approximation. To ensure closer resemblance with the experimental setup, microtubes were included in the geometry, using their inner diameter as the entrance size. To input a laminar flow at the inlet of the flow chamber, the entrance length was calculated for each flow rate with the expression for the Reynolds number in a circular pipe, resulting in values from 1 to 2 mm. Laminar flow rates were varied from 0.2 to 5.1 mL/min at the inlet. A no-slip boundary condition was used for the edges as well as a no-pressure condition for the outlet of the flow chamber. For all the simulations, the fluid flows from right to left.
Micro-particle Image Velocimetry
The experimental flow inside the imaging slit was monitored with a micro-particle image velocimetry (μPIV) setup composed of a machine vision camera (Grasshopper3, FLIR) and a variable magnification zoom lens system (MVL-7000, NAVITAR). The camera was equipped with a monochrome 2.3-megapixel Sony IMX174 CMOS sensor with a pixel size of 5.86 μm and a maximum acquisition rate of 163 fps. The zoom was a macro lens featuring a variable focal length from 18 to 108 mm. The maximum magnification is ~1.1X at its minimum working distance of 130 mm. A volume of 200 μL was taken from a 1% w/v stock solution of polystyrene microspheres with a 104 μm average diameter (PPB-1000-5, Spherotech). Then, the microspheres were suspended in 10 mL of a 1X PBS solution containing BSA at 1 mg/mL. The solution of suspended microspheres was infused with a perfusion pump (PPS-1100, Scientifica) at equivalent continuous flow rates ranging between 0.2 and 5.1 mL/min. In total, approximately three videos per flow rate, each lasting 3 s, were recorded for all four sample holders. The TrackMate plug-in of ImageJ  was used to follow the trajectory of the beads in the μPIV data to obtain the experimental streamlines inside the imaging slit.
Primary rat neurons were obtained from dissected brain hippocampus from neonatal rats (P0-P1), as described previously [44, 48], then plated on 18-mm borosilicate coverglass coated poly-D-lysine (GG-18-1.5, Neuvitro) and cultured for 2 weeks in Neurobasal media supplemented with B27 (LS21103049, Gibco™) to obtain mature neurons forming effective networks. Two cell lines, NIH3T3 embryonic mouse fibroblast cells (ATCC® CRL-1658™) and human fetal HEK293T cells (ATCC® CRL-11268™) were also used in this study. Briefly, NIH3T3 and HEK293T cells were seeded 12–24 h before imaging on 18-mm coverglass coated with matrigel (C354277, Corning™) at a density of 25,000 cells or 50,000 cells per well, respectively. Both cell lines were kept in high-glucose DMEM (10-017-CV, Corning™) supplemented with 10% bovine calf serum (97068-085, VWR) or 5% of premium grade fetal bovine serum (10158-358, VWR), respectively.
Neurons Viability Assay
A total of 10 cell cultures of rat primary neurons were prepared in two experimental groups. The five cultures of the control group were kept in a petri dish (08-772B, Fisher Scientific) and immersed in mannitol-based physiological solution without perfusion for 30 min at room temperature. The five cell cultures of the test group were placed in the flow chamber and were perfused at 2 mL/min with the same physiological solution at room temperature for 30 min. All coverslips were then stained for 3 min with a 0.04% Trypan Blue solution (MT25900CI, Fisher Scientific). Fourteen images at a 20X-magnification randomly positioned along the central line of the imaging slit, and two at a 4X-magnification centered in the imaging slit have been recorded per coverslip with a transmitted light microscope (EVOS FL, Invitrogen). Negatively- and positively-stained cells on each image were identified and manually counted with ImageJ (v1.52p, NIH) to assess the number of living and dead cells. Overall cell viability rates were calculated on the randomly positioned 20X-magnification images as they are representative of the whole cell culture, whereas the 4X-magnification images were meant to quantify cell viability in a position-dependent manner since the edges of the imaging slit are visible. In this case, images were separated into four segments parallel to the length of the imaging slit, and the living and dead cells were then separately counted in each segment.
Quantitative-Phase Digital Holographic Microscopy
A digital holographic microscope (T-1003, Lyncée Tec) based on a Mach-Zehnder interferometer was used. Within the microscope, the beam produced by a 666-nm laser diode is separated into two, forming an object wave and a reference wave. The object wave, scattered by the sample, interferes with the reference wave, producing an off-axis hologram recorded by a digital camera. The microscope objective magnification is 20X with a free working distance of 590 μm, and a NA of 0.70 (HC PL APO, Leica). The camera is equipped with a monochrome 2.3-megapixel Sony IMX174 CMOS sensor with a pixel size of 5.86 μm and capable of recording at 164 fps (acA1920-155um, Basler). The digital holograms were integrated over a period of ~100 μs and acquired at a frequency of 4 Hz. The integration time was adjusted carefully at the beginning of each imaging session to avoid any saturation in the holograms. The off-axis geometry (meaning that holograms are recorded with the object and reference waves having a small angle between their directions of propagation) allows the reconstruction of quantitative-phase images from a single recorded hologram . The reconstruction of the quantitative-phase image from the hologram is numerically computed using proprietary software (Koala, Lyncée Tec) and the resulting field-of-view is 215 μm by 215 μm. The reconstruction algorithm consists of a simulation of the illumination of the recorded hologram by a digital reference wave . This is followed by a numerical correction of the wave-front aberration induced by the optical components, the objective, and the off-axis geometry of the microscope using a reference-conjugated hologram . As for background flattening, it has been performed by means of polynomials .
Two-Liquid Decoupling Approach
The two-liquid decoupling approach requires sequentially perfusing two distinct physiological solutions with the same osmolarity and different, but a priori, known refractive indices. Concretely, a small amount of NaCl (25 mM) of an iso-osmotic and isotonic medium with an osmolarity appropriate for the cell under study (i.e., 220 mOsm/kg for neurons or 340 mOsm/kg for HEK293T and NIH3T3 cells) was replaced by either 50 mM of mannitol or 50 mM of Nycodenz to obtain the two different perfusion media . In practice, the refractive index and osmolarity of every physiological solution has been measured with an Abbe refractometer (Abbemat MW, Anton Paar) at 666 nm and a freezing-point osmometer (Micro-Sample model 3320, Advanced Instruments), respectively, prior to the QP-DHM session. The refractive index of the mannitol-based perfusion medium is ~1.3335, whereas for the Nycodenz-based physiological solution, it is 1.3380. Finally, consecutive perfusion of these physiological solutions, using a perfusion gravity system linked to an 8-port manual distribution valve , was performed either successively for three times, or before and after the cells were submitted to an osmotic challenge. For all the physiological solutions, the flow rate of the perfusion was set to 2 mL/min.
Three different cell types, primary hippocampal neurons, NIH3T3 cells, and HEK293T cells mounted in at least four different sample holders, were submitted to hypo-osmotic or hyper-osmotic challenges. In these experiments, the primary neurons were perfused with HEPES-buffered physiological solution with an osmolarity similar to Neurobasal + B27® medium containing (in mM) 75 NaCl, 5.4 KCl, 1.8 CaCl2, 0.8 MgCl2, 0.9 NaH2PO4, 10 HEPES, 5 dextrose (D-glucose), 50 mannitol, or Nycodenz (pH 7.4; 225-230 mOsm). The other cell cultures (HEK293T and NIH3T3) were maintained in DMEM which has a higher osmolarity (330 mOsm), and thus were perfused with an appropriate iso-osmotic solution (in mM) 135 NaCl, 3 KCl, 5 dextrose, 10 HEPES, 3 CaCl2, and 1.8 MgCl2, 50 mannitol or Nycodenz (pH 7.4; 330–335 mOsm). The hypo-osmotic media were prepared by removing 25 mM of NaCl from the iso-osmotic solutions, and the hyper-osmotic media by adding 25 mM of NaCl to the iso-osmotic solutions, resulting in refractive index alterations of −0.0003 and 0.0003, respectively, from their iso-osmotic counterpart. All cells were perfused with the adequate mannitol iso-osmotic medium and the QP-DHM recording started when the QPS was stable, typically after 2 min. Usually, the baseline was recorded for ~1 min, and then replaced by the Nycodenz iso-osmotic solution required for decoupling the resting state of the cells for ~30 s. The cells were then submitted to the osmotic challenge by either perfusing for ~5 min with either the hypo-osmotic or the hyper-osmotic solution containing Nycodenz. Finally, after the osmotic challenge was completed, a second decoupling sequence was performed by changing the medium for the mannitol-based solution but having the same osmolarity as that used for the challenge, for ~1 min.
Post-processing of Quantitative-Phase Images
Rectifications are commonly applied during the numerical reconstruction of the holograms to correct phase aberrations and image distortions. However, some irregularities may nevertheless persist in the quantitative-phase images. The reference-conjugated hologram and polynomials (both determined at the beginning of the experiment) gradually degrade and lose some of their effectiveness over time. Therefore, a spline smoothing method was used to iteratively determine an illumination function that when subtracted from the quantitative-phase image, effectively offsets the background fluctuations. This correction was individually applied to each quantitative-phase image in a time series as shown in reference . Thereafter, a global threshold-based segmentation strategy, implementing a minimum cross-entropy algorithm, was used to classify each pixel of the quantitative-phase image as either belonging to the background or to the foreground. In addition, an intensity-based method was used to distinguish between individual objects in the foreground that are touching one another and therefore not properly delineated as two entities by the segmentation alone. Finally, the segmented cells were numerically labeled throughout the time-lapse of quantitative-phase images to perform a cell-by-cell analysis.
Results and Discussion
As previously mentioned, an effective and accurate two-liquid decoupling approach is required to obtain a rapid and homogeneous return of the extracellular refractive index toward the known values of either nm1 or nm2 over the entire imaging slit of the flow chamber. Correspondingly, it is well-known in fluid mechanics that a flow system in which the pipe diameter changes abruptly from one size to another gives rise to what is referred to as vena contracta. The loss coefficient for a sudden size variation is a function of the geometry, with a sharp-edged change having more losses. Accordingly, the flow is said to separate from the sharp corner. The separated flow location is then the host of a swirl that traps the liquid and creates a counter-current flow, or an eddy. The size of such an eddy is flow rate dependent; in a faster flow, it spans a larger zone. In the context of the two-liquid decoupling approach (where two different physiological solutions are sequentially perfused), these dead volumes and unstirred layers are detrimental to obtaining a fast-homogeneous return of the extracellular refractive index toward either nm1 or nm2. Most available perfusion chambers have abrupt size expansions, e.g., a square-edged entrance and exit. A concrete example is the junctions between the microtubes for the perfusion and the imaging slit. These locations give rise to important losses causing separated flows containing eddies even at very low flow rates as confirmed by the CFD simulation shown in Figure 1. This simulation corresponds to a flow rate of 2 mL/min at the entrance of a rectangular imaging slit with a square-edged entrance and exit of same size as the one developed and described hereafter. Thus, a rectangular imaging slit with a square-edged entrance and exit does not allow for a fast and homogeneous washout. However, from a fluid mechanics point of view, obtaining a steady-state and fully-developed laminar flow is expected to be beneficial. Therefore, the geometry of the flow chamber has been optimized to provide such a flow field.
Figure 1. Separated flow and eddies. Flow simulation of a 2 mL/min flow rate through a rectangular imaging slit without diffusers. Streamlines are shown in black and are overlaid on the velocity field with values indicated by the color bar. The formation of a vena contracta and eddies are made evident with the separation of the flow. The flow direction is specified with the black arrow.
The exploded view of the proposed flow chamber is shown in Figure 2A. The flow chamber was 3D-printed with a plastic filament intended for medical and clinical applications since the raw material itself is biocompatible. The raw material is listed with the Food and Drug Administration and has complied with specific biocompatibility tests, namely ISO 10993, including cell cytotoxicity . The sample holder contains the channel for the flow and is encapsulated by two 18-mm coverslips. The top coverglass acts as an observation window while the bottom coverslip mounts the cell culture. The parafilm gasket is placed underneath the bottom coverglass to ensure waterproofness. Stainless steel microtubes with 0.4-mm internal diameters and 0.8-mm outer diameters, are partially inserted up to the entrance of the channel in the central element to convey the perfusion medium. The base plate holds the sample holder. The entire chamber is closed using the lid and four rare-earth magnets. For sizing purposes, a mounted flow chamber (top) and a sample holder (bottom) are shown in Figure 2B.
Figure 2. Mechanical design. (A) An exploded CAD-model view representing the flow chamber and detailing its constituents. (B) Photograph showing the assembled flow chamber (top) and a sample holder (bottom) as well as a ruler, for sizing purposes. (C,D) CAD models of the sample holder showing, respectively, a vertical and horizontal cross-section of the channel highlighting the diffusers where the outlet of the stainless steel microtubes will be attached to the imaging slit.
The central element, as shown in Figures 2C,D and being, respectively, a vertical and horizontal cross-section of the channel, is the key component of the flow chamber and is described in more detail in this paragraph. The channel changes smoothly from a circular cross-section of 0.4 mm in diameter (at the outlet of the stainless steel microtubes), to a rectangular cross-section 1.2-mm high and 2.0-mm wide at the level of the imaging slit using pyramidal diffusers. The opening and closing angles of these pyramidal diffusers, highlighted in Figures 2C,D, have been optimized and this aspect will be discussed in Section Steady-State and Fully-Developed Laminar Flow. The 9.0-mm long imaging slit has an area of 18.0 mm2 and an internal volume of 28.8 μl. The internal volume has been voluntarily minimized to enable fast washouts. Concretely, the area of the imaging slit is ~400 times that of a typical field of view for a microscope objective having a 20x-magnification. The thickness of the sample holder cannot be further reduced because the microtubes are inserted into it. Based on a daily usage in the laboratory, the lifetime of the central element is ~3 months after which leaks may occur or a degradation of the fluidic performances may be experienced. While the lid and base outlast the lifetime of the sample holder, the gasket is designed for single use. The fabrication cost of the entire flow chamber is 12 USD, of which 0.86 USD is for the plastic itself. Printing all flow chamber parts takes ~4 h. Gluing the tubes to the sample holder and the magnets to the mounts takes 15 min. The acetone treatment takes ~30 min to complete. It is advised to wait at least 2 h after the treatment to test its water tightness. The entire process takes ~7 h, 6 of which are spent waiting. The design files of the do-it-yourself flow chamber will be shared with the scientific community through an open-source model for rapid dissemination and maximum impact. The technical drawings will be available on the web page of the authors' research group upon publication.
Steady-State and Fully-Developed Laminar Flow
A simple way to reduce the losses mentioned above is to round off the entrance and exit regions of the imaging slit which connects it to the channels. A geometry such as this, shaped to decelerate the fluid, is called a diffuser. There is an optimum angle for which the loss coefficient is minimum, which is usually around 10°. The relatively small value of this angle for minimal losses results in a long diffuser and is an indication of the fact that it is difficult to efficiently decelerate a fluid . Thus, 6.5-mm long pyramidal diffusers with a total opening and closing angle of 11.45° have been added at the ends of the imaging slit as shown in Figures 2C,D. Such an angle, close to the optimal one, was chosen to have both minimal losses and an imaging slit of a size compatible with the resolution of the 3D printer. To evaluate the proposed design, CFD simulations have been computed for 16 flow rate values ranging from 0.2 to 5.1 mL/min. As an indication, the results shown in Figures 3A–C, are zoomed in on the imaging slit and the fluid flows from right to left. For synthesis purposes, only 3 of the 16 flow rates are illustrated. The field of streamlines (the black solid lines) for 1.0, 2.0, and 2.8 mL/min, superimposed on their respective velocity distributions and color coded, are shown in Figures 3A–C (top). The simulations show that fluid flows in constant-velocity parallel layers along a path parallel to the wall of the imaging slit. Adjacent layers exhibit slightly different speeds according to a velocity gradient perpendicular to the streamlines. This behavior is typical of a steady-state and fully-developed laminar flow. However, at a relatively high flow rate of 2.8 mL/min, there is a small region at the entrance where separated flow appears. For the specific application presented here, it has been noted that the Reynolds number is not a sufficiently sensitive metric to predict the presence of vortices and backward flows. As an example, at a flow rate of 5.1 mL/min, the flow chamber is completely filled with eddies, has revealed by the streamlines. But, the maximum corresponding Reynolds number is 12, well below the threshold value of 2100 from which the flow supposedly stops to be laminar, and this by a factor of ~200 (data not shown). Therefore, it is inappropriate to optimize the geometry and the sizing of the flow chamber based on this measure.
Figure 3. Steady-state and fully-developed laminar flow. (A–C) Comparison of CFD flow simulation (top) and experimental μPIV (bottom) for the imaging slit of the flow chamber with pyramidal diffusers at three different flow rates: 1.0, 2.0, and 2.8 mL/min, respectively. (D) Laminar length for an imaging slit without diffuser (green line, triangles) were obtained solely with simulations. Laminar length within the imaging slit with diffusers as a function of flow rate obtained with CFD simulations (blue line, circles) and experimental μPIV (red line, squares). Arrow indicates the flow rates shown in (A–C).
To experimentally confirm the validity of the design, μPIV experiments have been conducted using the same 16 values of flow rates as for the CFD simulations as shown in Figures 3A–C (bottom). Again, only the same three flow rates are represented. Each yellow solid line corresponds to the trajectory of a single particle inside the imaging slit flow field, therefore providing an experimental way to obtain the streamlines inside the actual flow chamber. The experimentally retrieved streamlines correspond reasonably well with those resulting from the numerical simulations. For the flow rates of 1.0 and 2.0 mL/min, the streamlines are relatively straight and parallel while for a high flow rate of 2.8 mL/min the two areas of separated flow at the entrance are observed. This stresses that the actual flow of fluid inside the imaging slit is steady-state, fully-developed, and laminar.
Since the regions of separated flow contain an eddy that traps the liquid and creates a counter-current flow, they are detrimental to obtaining a washout process characterized by a rapid stabilization of the refractive index values around nm1 or nm2. To quantify the extent of this phenomenon, the reattachment length Lr (defined as the distance from the entrance of the imaging slit at which the flow resumes in the positive flow direction all over the cross-section) has been measured. In order to obtain a metric that represents the laminar length Ll expressed as a percentage of the total length Lt of the imaging slit, and therefore free from separated flow regions, Ll is defined as follows: Ll = 100(Lt − Lr)/Lt. For the design having pyramidal diffusers, Ll has been measured on both the simulated (blue line, circles) and experimental data (red line, squares), and for all the studied flow rates, as shown in Figure 3D. An arrow demonstrates the three flow rates shown in Figures 3A–C. The analysis shows that the imaging slit is entirely usable, Ll = 100%, for flows up to 2 mL/min after which Ll decreases monotonically. For comparison purposes, an Ll based on CFD simulations for a rectangular imaging slit with a square-edged entrance and exit, having the same size than the proposed design comprising pyramidal diffusers, is plotted in Figure 3D. In this case, Ll is already 90% at 0.2 mL/min, and it decreases, monotonically as a function of the flow rate. It is important to note that the rectangular design contains separated flows throughout its length, Ll = 0%, at the typical flow rate for cellular imaging of 2 mL/min. At this flow rate, the proposed design comprising pyramidal diffusers is still entirely laminar.
Another important aspect of the flow chamber design is its compatibility with live cell imaging. A cell viability assay has been performed on neurons; this cell type is known to be very sensitive to environmental conditions. As shown in Figure 4A, the percentage of living neurons in the control group is 84.74 ± 2.26% while the viability in the test group is 81.71 ± 4.91%. Both viability rates are close, but the slight trend is likely to result from the additional manipulation of assembling the chamber i.e., risking a brief exposure to the air when the cells are momentarily removed from their physiological solution and again when closing the chamber. Nevertheless, this extra step is required with all perfusion chambers whether homemade or commercial in nature. However, this slight decrease in viability rate does not preclude imaging healthy cell cultures. Since neurons are more sensitive to environmental conditions than most other cell types, higher viability rates are expected for HEK293T and NIH3T3 cells. Furthermore, cells lying in the different segments of the imaging slit, as seen on Figure 4B, exhibit similar viability rates, namely 80.74 ± 7.28%, 83.68 ± 5.37%, 85.59 ± 4.17%, and 81.97 ± 7.74%, from segment one through four, respectively. Therefore, there is no significant variation in cell viability between the distinct regions of the imaging slit. Cells close to the edges of the imaging slit are not preferentially damaged nor do they obtain any benefit with respect to their viability.
Figure 4. Neurons viability assay. (A) Average percentage of living neurons observed in the control experiment and in the imaging slit of the flow chamber. (B) Average percentage of living neurons in each segment of the imaging slit. Inset: top-view visualization of the separation of the four segments; lines one and four are positioned near the edges of the imaging slit.
Quantitative-Phase Signal Stabilization Time
The acquisition protocol of the two-liquid decoupling approach, repeated twice, is illustrated in Figure 5. Rectangles on the top illustrates the perfusion sequence of the various equi-osmolar solutions, namely mannitol- and Nycodenz-based. The black solid line represents the mean cellular QPS of the HEK293T cell delineated by the green contour line, , obtained through automated segmentation, shown as insets. The is sensitive to the alternate flow at 2 mL/min of the first physiological solution having a refractive index of nm1 followed by the second perfusion medium having a refractive index of nm2. For each of the two-liquid perfusion sequences, exhibits a transient state during which its value is not defined, followed by an evolution period toward a plateau. The time to obtain a stabilized (namely the QPS stabilization time) reflects the time to obtain a complete washout of the flow chamber. Since the refractive index of Nycodenz is higher than that of mannitol, nm2 > nm1, the Nycodenz phase plateaus (71.0 and 70.5°) are lower than the mannitol plateaus (86.6 and 85.6°). Accordingly, this highlights the fact that the perfusion medium bathing the cells is an active parameter of the quantitative-phase image formation process, as described by Equation (4). Coherently, the quantitative-phase images (shown as insets in Figure 5) corresponding to the Nycodenz phase plateaus, are darker than the mannitol plateaus especially in the cytoplasmic region, as confirmed by the color bars expressed in degrees. These quantitative-phase images correspond to the time-points defined by the arrowheads.
Figure 5. QPS stabilization time. Graph showing the mean cellular QPS of one HEK293T cell during the two-solution decoupling approach. The QPS is averaged over the cell body delineated by the green contour line, shown in the insets. The quantitative-phase images correspond to the time-points as defined by the respective arrowheads. Rectangles across the top illustrate the perfusion sequence of the mannitol- (dark gray) and Nycodenz-based (light gray) solutions. The mean cellular QPS as well as the color bars of the insets are expressed in degrees and the scale bars are 5 μm.
Experimentally, the typical fluctuations (on the scale of a few seconds), for which the automated segmentation algorithm is stable, are ~0.2°. Such phase fluctuations, although dependent on the coherent noise, result mainly from the natural cell micro-movements and deformations namely the cell-mediated noise . Considering these cell body phase fluctuations, the typical time to obtain an appropriate stabilization, to perform a two-liquid decoupling procedure with the 3D-printed flow chamber, is ~10.09 ± 3.46 and 12.79 ± 3.57 s from mannitol to Nycodenz, and vice versa, respectively. This determination results from an extensive analysis of 37 turnovers, randomly positioned over the entire imaging slit of the four different sample holders. Compared to the originally reported stabilization time of 30–60 s for the mannitol-Nycodenz transition, this represents a reduction by at least a factor of three to six [16, 25]. The fabrication parameter having the most restrictive tolerance to obtain the reported performances on the washout process is the concentric alignment of the axes of the channels to insert the microtubes. This tolerance is met by the micrometric resolution of the 3D printer described in this paper. Although QP-DHM has coherent noise originating from the degree coherence of the light source required to perform off-axis interferometry, from experimental observations, there are no significant speckles originating specifically from the potential remaining roughness of the internal walls of the imaging slit.
Keeping the QPS stabilization time as short as possible is important considering that even more significant cell changes, including important cell position and shape modifications involving thickness and/or refractive index, are very likely to take place over time causing important artifacts on , and Vc measurements. In general, the Nycodenz-mannitol QPS stabilization times are longer by a factor of ~1.3. In addition, the variability in the QPS stabilization time for a given pair of liquid is governed (for now) by the reproducibility of the fabrication process and aging of the sample holder. The rapidity of the QPS stabilization time is an indication of the homogeneity of the refractive index of the perfusion medium inside the imaging slit.
Thus, the presence of a steady-state and fully-developed laminar flow improves the washout process of the solutions resulting in a faster QPS stabilization time corresponding to the restauration of the known values nm1or nm2 of the two perfusion solutions. This makes the a priori knowledge of the refractive indices of the perfusion media a more reliable approximation than the one in situ. Moreover, a flow chamber with such a homogenous flow of fluid is well-suited to perform pharmacological experiments.
Measurements of the Whole-Cell Refractive Index and Absolute Cell Volume
As previously mentioned, the decoupling procedure measures the integral IRI over the OPL of each pixel and the corresponding cell thickness hc(x, y) separately from Δϕc(x, y). Thanks to the mean value theorem for integrals, , such a cell thickness measurement remains valid even if exhibits some variations along the OPL of each pixel, as long as the observed sample is a weakly-scattering object such as a living cell, i.e., one that does not significantly deviate optical rays from their incident direction. However, as mentioned above, during the washout processes, both cell-mediated noise as well as coherent noise (contributing to the major part of the phase fluctuations), preclude the possibility to obtain a precise value for hc(x, y) and measurements at the level of individual pixels. Nonetheless, the spatial phase averaging over a cell body  significantly decreases both effects leading to the precise measurement of and , and consequently of Vc according to Equations (7)–(9). These biophysical parameters, along with their respective errors obtained by the technique of partial derivatives, were calculated for populations of three cell types in culture including primary neurons, HEK293T cells, and NIH3T3 cells.
The fluctuations of the mean cellular QPS represent the main factor limiting the precision of and measurements. Typical fluctuations of 0.2°, as seen in Figure 5, used in conjunction with the technique of partial derivatives on Equations (7)–(8) lead to values of = 0.0003 and = 70 nm. Such a high precision on and measurements comes at the price of losing the spatial resolution of these two biophysical parameters. Interestingly, and , shown in Figure 6A, do not exhibit a significant decrease as a function of Sc. This is consistent with the fact that the main effect of the spatial averaging on the cell-mediated noise is already reached for surfaces ranging from 100 to 400 μm2 . Typically, across all the cells of the populations of the three cell types studied, the precisions at the single-cell level of and are 0.0006 and 100 nm, respectively. The measurements of Vc as a function of Sc are shown in Figure 7A. It is observed that the precision of Vc is drastically affected by Sc. This is detailed in Figure 6B presenting the absolute and relative errors of Vc as a function of Sc obtained with the technique of partial derivatives on Equation (9). While δVc increases for large cellular surfaces, the relative error decreases. According to Equation (9) and given the fact that does not significantly depend on Sc as shown in Figure 6A, it turns out that the influence of Sc on Vc should come from our capacity of measuring Sc with precision. This is consistently illustrated in Figure 6C showing the behavior of both δSc and as a function of Sc. Considering the SNR of , it is understood that the precision of Sc is mainly limited by the lateral resolution of the quantitative-phase images dictated by the NA of the microscope objective, which is 0.7 in this case. Therefore, the decrease of as a function of Sc reflects the ratio between the cell perimeter and its surface. The tendency of a behavior inversely proportional to the radius of the cell surface comes from the presence of a significant number of nearly round cells which is specially the case for the smallest cells of our populations. Figure 7B shows as a function of Vc. Interestingly, this highlights the fact that cells with smaller Vc tend to have a higher . This observation is consistent with a ratio between the solid and the cytosolic portion, which is higher for smaller cells .
Figure 6. Error of biophysical measurements. Absolute and relative errors of numerous biophysical measurements as a function of Sc for three cell types in culture, including primary neurons, and HEK293T and NIH3T3 cells. Data points plotted with darker squares are in relation to the left axis and data points plotted with lighter shade triangles are in relation to the right axis. (A) Plots of the absolute error (left) and absolute error (right), (B) absolute (left) and relative (right) Vc error, and (C) absolute (left) and relative (right) Sc error.
Figure 7. Biophysical parameters. Various biophysical measurements for three cell types in culture, including primary neurons, and HEK293T and NIH3T3 cells. (A) Graph of Vc as a function of Sc, (B) Scatter of as a function of Vc.
Table 1 shows the means and the standard deviations of the distributions of a set of biophysical parameters measured over populations of the three different cell types studied. Mean values and standard deviations, for cell populations, of whole-cell refractive index and average cell thickness of 1.3672 ± 0.0070 and 5.87 ± 2.38 μm for neurons (n = 64), 1.3637 ± 0.0086 and 5.27 ± 1.46 μm for HEK293T cells (n = 202), 1.3660 ± 0.0038 and 3.01 ± 0.96 μm for NIH3T3 cells (n = 14) were obtained, respectively, for resting state. Mean values and standard deviations, for cell populations, of absolute cell volumes V of 1,210 ± 890, 1,460 ± 840, and 2,200 ± 480 μm3 were obtained for the populations of neurons, and HEK293T and NIH3T3 cells, respectively. The standard deviations reported in Table 1 relate to the spread of the distribution of the different biophysical parameters, reflecting the heterogeneity across the cell populations of a given cell type and not the precision of a single-cell measurement. Interestingly, in the conditions of cell culture used here, NIH3T3 cells are thinner by a factor of two, but have greater V because of their larger Sc. In addition, it seems that the distribution of is rather narrow while for V it is much broader. It is interesting to note that the absolute values of the cell biophysical parameters presented here are consistent with existing literature, even when considering various measurement techniques [26, 29, 36, 55].
Table 1. Measured biophysical parameters on populations of primary neurons, HEK293T, and NIH3T3 cells during their resting state.
The membranes of animal cells are highly permeable to water. Water movement across membranes is in large part dictated by osmotic pressure gradients. During a hypo-osmotic challenge, a net water influx occurs, causing a cellular swelling. On the contrary, in a hyper-osmotic challenge, there is a water efflux, leading to cell shrinkage. Changing the osmolarity of the extracellular medium to perform osmotic challenges implies the use of at least two solutions, a number that is doubled in the case of assessing both as well as Vc through the two-liquid decoupling approach. Therefore, these experiments offer an excellent paradigm to confirm the efficiency and performance of the proposed flow chamber.
The responses of two cell types in culture (namely the neurons and HEK293T cells), to osmotic challenges, have been studied. First, a typical response of the mean cellular QPS (black solid line) for a HEK293T cell during a hypo-osmotic challenge is shown in Figure 8. The mean cellular QPS results from an averaging performed over the cell body delineated by the green contour lines, , as shown in insets. These quantitative-phase images correspond to the time-points defined by the arrowheads. Rectangles on the top illustrate the perfusion sequence of the various perfusion solutions, namely the iso-osmotic mannitol (blue), iso-osmotic Nycodenz (red), hypo-osmotic Nycodenz (light red), and hypo-osmotic mannitol (light blue). As described in Section Two-Liquid Decoupling Approach, the strength and duration of the hypo-osmotic stress are 50 mOsm and 5 min, respectively. The vertical gray dashed line illustrates the onset of the hypo-osmotic challenge.
Figure 8. Hypo-osmotic challenge. A graph of the mean cellular QPS of one HEK293T cell during a hypo-osmotic challenge is shown. The QPS is averaged over the cell body delineated by the green contour line, shown in the insets. The quantitative-phase images correspond to the time-points defined by the respective arrowheads. Rectangles across the top illustrate the perfusion sequence of the various solutions, namely the iso-osmotic mannitol (blue), iso-osmotic Nycodenz (red), hypo-osmotic Nycodenz (light red), and hypo-osmotic mannitol (light blue). The vertical gray dashed line indicates the onset of the hypo-osmotic challenge. The mean cellular QPS as well as the color bars of the insets are expressed in degrees and the scale bars are 5 μm.
As discussed in Section Quantitative-Phase Signal Stabilization Time, the first two plateaus of are used to obtain the and values during the cell resting state; in this specific case their values are 1.3737 ± 0.0007 and 4.50 ± 0.12 μm, respectively. This given , used in conjunction with the corresponding Sc of 247 ± 5 μm2, translate to a Vc of 1,110 ± 30 μm3. As soon as the hypo-osmotic Nycodenz solution is perfused, is rapidly decreased to finally reach a stabilized state after around 60 s of challenge, corresponding to an equilibrium between intra- and extra-cellular osmotic pressures. Next, the perfusion of the hypo-osmotic mannitol medium leads to the fourth plateau. These last two plateaus of are used to calculate = 1.3670 ± 0.0006 and = 5.26 ± 0.11 μm values corresponding to the new steady-state induced by the hypo-osmotic challenge, namely a cellular swelling. For this particular cell, the hypo-osmotic challenge effectively mediated a decrease of , given by = −0.007 ± 0.001, and an increase of , corresponding to = 0.8 ± 0.2 μm. The hypo-osmotic challenge caused a cell volume increase of ΔVc = Vchypo−Vc = 190 ± 70 μm3 and a swelling factor β = Vchypo/Vc = 1.17 ± 0.06.
For a hyper-osmotic challenge, the experiment is virtually the same as in Figure 8 except that the two last perfusion media were different, namely hyper-osmotic Nycodenz and hyper-osmotic mannitol. Similarly, biophysical data are obtained from the values of the different plateaus of . In summary, Table 2 shows the hypo-osmotic and hyper-osmotic induced changes of a set of biophysical parameters expressed as the mean value and standard deviation calculated over population of two cell types in culture, including neurons (n = 16 and n = 9, respectively) and HEK293T cells (n = 16 and n = 25, respectively). Here also, the standard deviations reported in Table 2 relate to the range of the values of the different biophysical parameters, reflecting the heterogeneity across the cell population of a given cell type and not the precision of a single-cell measurement.
Table 2. Measured variations of biophysical parameters induced by osmotic challenges on populations of primary neurons and HEK293T cells.
The data presented in Table 2 suggests that the biophysical parameter is sufficient to quantify water movements over a cell population. The hypo-osmotic challenge induces a decrease of at the level of a population of cells. This reduction is consistent with a hypo-osmotic water influx, resulting in a dilution of the intracellular protein concentration—the cellular component that largely determines the value of . The contradictory fall of results from a lateral dilution of the dry mass concentration. The explanation is also valid for the case of a hyper-osmotic challenge. In addition, both the and the ΔV as well as the β values obtained here are consistent with the existing literature [26, 27, 29]. The measurements of these biophysical data for both osmotic conditions reinforce our understanding of the cellular mechanisms of these phenomena.
In conclusion, a low-cost, open-source, and do-it-yourself 3D-printed flow chamber optimized to efficiently perform dissociated measurements of and , resulting in accurate calculations of the critical cell parameter Vc with QP-DHM and with QPI in a more general manner, has been developed. The flow chamber has been 3D-printed with a plastic filament for which the raw material itself is biocompatible. The entire flow chamber costs <12 USD to fabricate, of which a mere 0.86 USD is for the plastic itself. The compatibility of the flow chamber with live cell imaging has been demonstrated by quantifying its cell viability rate when compared to a control condition. The assay revealed similar cell viability rates in both conditions, namely 81.71 ± 4.91% for the flow chamber and 84.74 ± 2.26% for the control condition. The geometry of the flow chamber has been customized to have a steady-state and fully-developed laminar flow which optimizes the washout process. The presence of such a flow of fluid inside the imaging slit has been verified both numerically and experimentally. Specifically, following the two-liquid decoupling sequence, a state where the refractive index of the solution is homogeneously distributed over the entire imaging slit can be rapidly obtained. This can be assessed by experimentally measuring the time to obtain a stabilized QPS. The mean QPS stabilization time of 37 turnovers resulting from four different sample holders has been quantified and found to be 10.09 ± 3.46 and 12.79 ± 3.57 s, depending on the given sequence of liquids used. Such a washout optimization has permitted to measure with high reliability the important cell biophysical parameters including , , and Vc and results in a faster stabilization time, thereby resulting in reduced cell changes and the resulting artifacts. For a typical cell body having a Vc of 1,250 μm3 and a Sc of 300 μm2, , , and Vc values are usually measured with a precision of 0.0006, 100 nm, and 50 μm3, respectively. Furthermore, the variations of , , and Vc have been quantitatively calculated during hypo-osmotic and hyper-osmotic challenges. Hypo-osmotic challenges of 50 mOsm on cell populations have shown globally a dilution of , on an average of −0.0115, as well as an increase of , typically of 1.50 μm, and a gain of V, characteristically of 290 μm3. As expected, changes in the opposite directions, reflecting a reversible process, have been observed in the case of hyper-osmotic challenges. Finally, the design files of the do-it-yourself flow chamber will be shared with the scientific community through an open-source model for rapid dissemination and maximum impact.
Data Availability Statement
The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.
EB and PM designed the study, wrote the paper with help from SL, ÉR-P and PL, and co-supervised the research. SL, ÉR-P, and PL conducted experiments. SL and PL processed the data. ÉR-P analyzed the data.
This work was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC) (RGPIN-2018-06198), the Canada Excellence Research Chairs (CERC) program and the Canada Foundation for Innovation (CFI) (36689, 34265).
Conflict of Interest
PM declares the presence of a potential conflict of interest as a co-founder of Lyncée Tec, a company that commercializes digital holographic microscopes. However, the study has been performed independently of Lyncée Tec in academic laboratories.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors would like to acknowledge the contribution all former interns and staff members, namely Gabriel Anctil, Vincent Roy, Alyson Bernatchez, Valerie Watters, Anne-Sophie Poulin-Girard, Chloé Martel, Jean-Michel Mugnes, and Jean-Xavier Giroux for their contribution in the development of the flow chamber and post-processing of quantitative-phase images. The authors would also like to acknowledge the services of CMC Microsystems which facilitated this research.
4. Lang F, Ritter M, Gamper N, Huber S, Fillon S, Tanneur V, et al. Cell volume in the regulation of cell proliferation and apoptotic cell death. Cell Physiol Biochem. (2000) 10:417–28. doi: 10.1159/000016367
8. Gregg JL, McGuire KM, Focht DC, Model MA. Measurement of the thickness and volume of adherent cells using transmission-through-dye microscopy. Pflugers Arch Eur J Physiol. (2010) 460:1097–104. doi: 10.1007/s00424-010-0869-2
12. Marquet P, Rappaz B, Magistretti PJ, Cuche E, Emery Y, Colomb T, et al. Digital holographic microscopy: a noninvasive contrast imaging technique allowing quantitative visualization of living cells with subwavelength axial accuracy. Opt Lett. (2005) 30:468–70. doi: 10.1364/OL.30.000468
15. Bon P, Maucort G, Wattellier B, Monneret S. Quadriwave lateral shearing interferometry for quantitative phase microscopy of living cells. Opt Express. (2009) 17:13080–94. doi: 10.1364/OE.17.013080
16. Moon I, Jaferzadeh K, Ahmadzadeh E, Javidi B. Automated quantitative analysis of multiple cardiomyocytes at the single-cell level with three-dimensional holographic imaging informatics. J Biophotonics. (2018) 11:e201800116. doi: 10.1002/jbio.201800116
17. Chmelik R, Krizova A, Kollarova V, Vesely P, Collakova J, Slaba M, et al. Coherence-controlled holographic microscopy enabled recognition of necrosis as the mechanism of cancer cells death after exposure to cytopathic turbid emulsion. J Biomed Opt. (2015) 20:111213. doi: 10.1117/1.JBO.20.11.111213
21. Rappaz B, Cano E, Colomb T, Kühn J, Depeursinge C, Simanis V, et al. Noninvasive characterization of the fission yeast cell cycle by monitoring dry mass with digital holographic microscopy. J Biomed Opt. (2009) 14:034049. doi: 10.1117/1.3147385
23. Rappaz B, Barbul A, Hoffmann A, Boss D, Korenstein R, Depeursinge C, et al. Spatial analysis of erythrocyte membrane fluctuations by digital holographic microscopy. Blood Cell Mol Dis. (2009) 42:228–32. doi: 10.1016/j.bcmd.2009.01.018
24. Boss D, Hoffmann A, Rappaz B, Depeursinge C, Magistretti PJ, Van de Ville D, et al. Spatially-resolved eigenmode decomposition of red blood cells membrane fluctuations questions the role of ATP in flickering. PLoS ONE. (2012) 7:e40667. doi: 10.1371/journal.pone.0040667
25. Javidi B, Markman A, Rawat S, O'Connor T, Anand A, Andemariam B. Sickle cell disease diagnosis based on spatio-temporal cell dynamics analysis using 3D printed shearing digital holographic microscopy. Opt Express. (2018) 26:13614–27. doi: 10.1364/OE.26.013614
26. Rappaz B, Marquet P, Cuche E, Emery Y, Depeursinge C, Magistretti PJ. Measurement of the integral refractive index and dynamic cell morphometry of living cells with digital holographic microscopy. Opt Express. (2005) 13:9361–73. doi: 10.1364/OPEX.13.009361
27. Curl CL, Bellair CJ, Harris PJ, Allman BE, Roberts A, Nugent KA, et al. Single cell volume measurement by quantitative phase microscopy (QPM): a case study of erythrocyte morphology. Cell Physiol Biochem. (2006) 17:193–200. doi: 10.1159/000094124
28. Rappaz B, Barbul A, Emery Y, Korenstein R, Depeursinge C, Magistretti PJ, et al. Comparative study of human erythrocytes by digital holographic microscopy, confocal microscopy, and impedance volume analyzer. Cytom A. (2008) 73:895–903. doi: 10.1002/cyto.a.20605
29. Boss D, Kühn J, Jourdain P, Depeursinge C, Magistretti PJ, Marquet P. Measurement of absolute cell volume, osmotic membrane water permeability, and refractive index of transmembrane water and solute flux by digital holographic microscopy. J Biomed Opt. (2013) 18:036007. doi: 10.1117/1.JBO.18.3.036007
30. Jourdain P, Allaman I, Rothenfusser K, Fiumelli H, Marquet P, Magistretti PJ. L-Lactate protects neurons against excitotoxicity: implication of an ATP-mediated signaling cascade. Sci Rep. (2016) 6:21250. doi: 10.1038/srep21250
31. Jourdain P, Becq F, Lengacher S, Boinot C, Magistretti PJ, Marquet P. The human CFTR protein expressed in CHO cells activates aquaporin-3 in a cAMP-dependent pathway: study by digital holographic microscopy. J Cell Sci. (2014) 127:546–56. doi: 10.1242/jcs.133629
32. Jourdain P, Boss D, Rappaz B, Moratal C, Hernandez MC, Depeursinge C, et al. Simultaneous optical recording in multiple cells by digital holographic microscopy of chloride current associated to activation of the ligand-gated chloride channel GABAA receptor. PLoS ONE. (2012) 7:e51041. doi: 10.1371/journal.pone.0051041
33. Jourdain P, Pavillon N, Moratal C, Boss D, Rappaz B, Depeursinge C, et al. Determination of transmembrane water fluxes in neurons elicited by glutamate ionotropic receptors and by the cotransporters KCC2 and NKCC1: a digital holographic microscopy study. J Neurosci. (2011) 31:11846–54. doi: 10.1523/JNEUROSCI.0286-11.2011
35. Yang SA, Yoon J, Kim K, Park YK. Measurements of morphological and biophysical alterations in individual neuron cells associated with early neurotoxic effects in Parkinson's disease. Cytom A. (2017) 91:510–8. doi: 10.1002/cyto.a.23110
36. Curl CL, Bellair CJ, Harris T, Allman BE, Harris PJ, Stewart AG, et al. Refractive index measurement in viable cells using quantitative phase-amplitude microscopy and confocal microscopy. Cytom A. (2005) 65:88–92. doi: 10.1002/cyto.a.20134
37. Kemper B, Carl D, Schnekenburger J, Bredebusch I, Schafer M, Domschke W, et al. Investigation of living pancreas tumor cells by digital holographic microscopy. J Biomed Opt. (2006) 11:034005. doi: 10.1117/1.2204609
39. Kemper B, Kosmeier S, Langehanenberg P, von Bally G, Bredebusch I, Domschke W, et al. Integral refractive index determination of living suspension cells by multifocus digital holographic phase contrast microscopy. J Biomed Opt. (2007) 12:054009. doi: 10.1117/1.2798639
40. Rappaz B, Charrière F, Depeursinge C, Magistretti PJ, Marquet P. Simultaneous cell morphometry and refractive index measurement with dual-wavelength digital holographic microscopy and dye-enhanced dispersion of perfusion medium. Opt Lett. (2008) 33:744–6. doi: 10.1364/OL.33.000744
41. Przibilla S, Dartmann S, Vollmer A, Ketelhut S, Greve B, von Bally G, et al. Sensing dynamic cytoplasm refractive index changes of adherent cells with quantitative phase microscopy using incorporated microspheres as optical probes. J Biomed Opt. (2012) 17:0970011. doi: 10.1117/1.JBO.17.9.097001
42. Dardikman G, Nygate YN, Barnea I, Turko NA, Singh G, Javidi B, et al. Integral refractive index imaging of flowing cell nuclei using quantitative phase microscopy combined with fluorescence microscopy. Biomed Opt Express. (2018) 9:88–92. doi: 10.1364/BOE.9.001177
43. Przybytek A, Gubanska I, Kucinska-Lipka J, Janik H. Polyurethanes as a potential medical-grade filament for use in fused deposition modeling 3D printers – a brief review. Fibres Text East Eur. (2018) 6:120–5. doi: 10.5604/01.3001.0012.5168
44. Lévesque SA, Mugnes JM, Bélanger E, Marquet P. Sample and substrate preparation for exploring living neurons in culture with quantitative-phase imaging. Methods. (2018) 136:90–107. doi: 10.1016/j.ymeth.2018.02.001
45. Dardikman G, Shaked NT. Review on methods of solving the refractive index–thickness coupling problem in digital holographic microscopy of biological cells. Opt Commun. (2018) 422:8–16. doi: 10.1016/j.optcom.2017.11.084
46. Bélanger E, Lévesque SA, Anctil G, Poulin-Girard A-S, Marquet P. Low-cost production and sealing procedure of mechanical parts of a versatile 3D-printed perfusion chamber for digital holographic microscopy of primary neurons in culture. In: Proceedings SPIE 10074, Quantitative Phase Imaging III. San Francisco, CA (2017). p. 100741S. doi: 10.1117/12.2253036
47. Tinevez JY, Perry N, Schindelin J, Hoopes GM, Reynolds GD, Laplantine E, et al. TrackMate: an open and extensible platform for single-particle tracking. Methods. (2017) 115:80–90. doi: 10.1016/j.ymeth.2016.09.016
48. Nault F, De Koninck P. Dissociated hippocampal cultures. In: Doering LC, editor. Protocols 799 for Neural Cell Culture. Totowa, NJ: Humana Press (2009). p. 137–59. doi: 10.1007/978-1-60761-292-6_8
49. Cuche E, Marquet P, Depeursinge C. Simultaneous amplitude-contrast and quantitative phase-contrast microscopy by numerical reconstruction of Fresnel off-axis holograms. Appl Opt. (1999) 38:6994–7001. doi: 10.1364/AO.38.006994
51. Colomb T, Kühn J, Charrière F, Depeursinge C, Marquet P, Aspert N. Total aberrations compensation in digital holographic microscopy with a reference conjugated hologram. Opt Express. (2006) 14:4300–6. doi: 10.1364/OE.14.004300
52. Colomb T, Cuche E, Charrière F, Kühn J, Aspert N, Montfort F, et al. Automatic procedure for aberration compensation in digital holographic microscopy and applications to specimen shape compensation. Appl Opt. (2006) 45:851–63. doi: 10.1364/AO.45.000851
54. Olsen H, Kinne-Saffran E, Kinne RKH, Wehner F, Tinel H. Cell volume regulation: osmolytes, osmolyte transport, and signal transduction. Rev Physiol Biochem Pharmacol. (2003) 148:1–80. doi: 10.1007/s10254-003-0009-x
55. Liu PY, Chin LK, Ser W, Chen HF, Hsieh CM, Lee CH, et al. Cell refractive index for cell biology and disease diagnosis: past, present and future. Lab Chip. (2016) 16:634–44. doi: 10.1039/C5LC01445J
Flow Chamber Fabrication Protocol
Fabrication of the Lid and Base
• Slice .stl files produced with SolidWorks of the lid and base in Cura, flat face down on the print bed.
◦ Use 0.4 mm nozzle and 30% infill
• Print with Taulman3D Guideline and remove from print bed.
• Glue magnets in the 4 slots of both parts.
◦ Make sure all magnets placed in the lid have inverted polarities to those placed in the base.
• Wait for the glue to dry.
Fabrication of the Sample Holder
• Slice sample holder .stl file produced with SolidWorks in Cura, flat face down on the print bed.
◦ Use 0.25 mm nozzle and full infill
◦ Use custom settings to increase infill overlap percentage to 30%
■ This setting can be added via: Preferences - > Settings - > Infill
◦ Reduce print speed to 30 mm/s
◦ Add brim for better adhesion to the print bed
• Change nozzle on printer to 0.25 mm and precisely calibrate the glass print bed.
• Print with Taulman3D Guideline using their recommendation for nozzle and bed temperatures.
• Wait for the bed to cool almost completely to remove.
• Cut micro-tubes to have two segments of ~2.5 cm.
◦ With small gage wire, verify the openings of the microtube segments.
• Glue micro-tubes in slots on each side of the sample holder.
◦ To prevent closing the tubes with the glue, pass small gage wire through them during the procedure.
• Wait for glue to dry.
• Treat with acetone as described in the reference of the main text.
• Press top coverslip on the sample holder directly after acetone treatment.
◦ The top of the holder can be distinguished by its small round indentation.
◦ This coverslip might break or fall during extensive use and should be replaced. This replacement can be greased upon the surface around the slit.
• Mount the flow chamber and test water tightness. The steps to mount the chamber are described in the next section (use clean coverslip instead of cell culture).
Mounting the Chamber With Cell Cultures
• Position the holder with the flattest face upwards.
• Slowly fill the channel with cell or perfusion medium. This will prevent the cells from dying when the coverslip is added.
• Place coverslip with cells on the sample holder.
• Clean off excess cell medium.
• Press the parafilm gasket onto the coverslip and holder, and make sure it sticks to the surface.
• Place sample holder in the base.
Place the lid slowly to avoid breaking the coverslips with the force of the magnets.
Keywords: quantitative-phase imaging, quantitative-phase digital holographic microscopy, flow chamber, whole-cell refractive index, mean cell thickness, absolute cell volume, biophysics
Citation: Bélanger E, Lévesque SA, Rioux-Pellerin É, Lavergne P and Marquet P (2019) Measuring Absolute Cell Volume Using Quantitative-Phase Digital Holographic Microscopy and a Low-Cost, Open-Source, and 3D-Printed Flow Chamber. Front. Phys. 7:172. doi: 10.3389/fphy.2019.00172
Received: 28 February 2019; Accepted: 17 October 2019;
Published: 06 November 2019.
Edited by:Pietro Ferraro, Italian National Research Council (CNR), Italy
Reviewed by:Emiliano Descrovi, Politecnico di Torino, Italy
Juergen Schnekenburger, University of Münster, Germany
Teresa Cacace, University of Campania Luigi Vanvitelli, Italy
Copyright © 2019 Bélanger, Lévesque, Rioux-Pellerin, Lavergne and Marquet. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Erik Bélanger, firstname.lastname@example.org