|
1.IntroductionIn recent years, the interest emerged in the development of new optical methods for individual cell analysis. A significant progress took place in interference methods for evaluation of red blood cells.1–5 These techniques allow for highly accurate measurement of critical parameters of the cell, that is, the refractive index, the geometric volume, and the spectrum of fluctuations. In particular, the refractive index (or refractivity, i.e., the difference between the refractive indices of the object and the environment) is one of the most informative parameters of an individual cell.1–5 Its diagnostic importance has been demonstrated by characteristic changes of erythrocytes infected with Plasmodium falciparum.2,3 The approaches based on holography,4,5 phase contrast, and common-paths method6 have been applied to measure the refractive index by interference microscopy. In so doing, a variety of algorithms, models,1–10 and several methods of image recording can be combined in a single instrument. The method of spatial interference microscopy (SLIM-fluorescence multimodal imaging) developed by Popescu and colleagues11 allows for simultaneous recording of structural elements and dynamic processes in neurons. Common for all interference methods is the representation of the images in the form of a 2-D distribution of optical path difference (OPD) in the image plane. The most important and the unique advantage of the interference microscopy is the measurement of OPD in absolute units of wavelength1–11 with very high precision, which is hardly achievable or impossible with other methods. In addition, the common-path method allowed to significantly decrease the sensitivity of images to external acoustic noise.8 Modern models of interference microscopes allow for measuring OPD based on various modulation methods and algorithms.3–5,8,9 In particular, the measurements were made with a raster OPD-compensation method providing a significant improvement of resolution.9,12,13 Most studies using interference microscopy have been performed with red blood cells,1–5 whereas other biological objects have been seldom addressed.7,10–12 One reason for this preference might be optical homogeneity of erythrocytes. Because these cells lack optically contrast elements, e.g., the nuclei and organelles, erythrocytes may be well characterized by mean refractive index.1–5 However, the mean volume-averaged refractivity is insufficient for characterization of optically heterogeneous mononuclear cells, such as lymphocytes. A significantly more complex algorithm and much higher spatial resolution are required for the analysis of intracellular structures. Furthermore, a problem arises as to how to interpret the phase images of the optically complex object. These issues have been discussed with respect to the refractivity of single mitochondria, chloroplasts, cyanobacteria, and human cells by coherent phase microscopy.7,9,10,12 In this study, we propose an approach to characterize an individual cell, namely the lymphocyte, by its phase portrait, which is a set of physical parameters of subcellular structures. We believe that such a portrait can be supplemented by data obtained with other methods and serve for evaluation of the physiologic state of the cell. 2.Experimental DetailsMeasurements were performed on a microscope Airyscan.7,9,12,13 The cells were illuminated with a collimated beam of light. In the phase images no speckles were observed, in spite of relatively big coherence length of the He-Ne laser (, 1 mW) used for illumination. The OPD values were determined sequentially in each pixel of the phase image by the compensation method. The phase of the reference wave varied linearly and periodically with a movable mirror in the reference arm of the interferometer. The image raster-scanning was performed by a coordinate sensitive photodetector, that is, the dissector image tube LI-620. This method allowed us to get super-diffraction spatial resolution,9,13 and to record cellular processes in real-time.12 The sensitivity to OPD changes limited by noise was ; therefore, the characteristics of the interference microscope did not seriously limit the results of measurements. The suspension of -lymphocytes isolated from peripheral blood of healthy donors was provided by I. Vasilenko (Moscow Institute of Clinical Research). Freshly isolated cells were resuspended in medium 199 (Sigma-Aldruch) supplemented with 5% fetal calf serum (HyClone) and incubated at 37°C for 1 h without (quiescence) or with phytohemagglutinin (PHA; Sigma-Aldrich activation). To prepare the sample for the analysis, a drop of the cell suspension was applied on a polished silicon wafer, covered with a cover slip, and placed on a table of the coherent phase microscope Airyscan. The position of the sample in the chamber is shown schematically in Fig. 1(a). A 15-μm-thick spacer between the silicon substrate and the cover slip prevented cell deformation. Measurements of each sample were conducted within 20 min. A removable screen was used to minimize the effect of laser radiation on cells. Fifty individual cells were analyzed per each sample. 3.Results and Discussion3.1.Concept of Zones and Characteristic Parameters in the Phase Image of the CellNormally, the phase image is described as the function:7–9 in the image plane , where and are the refractive indices of the object and the environment, respectively. The function can be interpreted as a “projection” of the refractive index onto the image plane. Under certain conditions, the inverse problem can be solved when the function is restored from its projection.1–6,14A fundamentally different approach involves the determination of the characteristic parameters of the function . We analyzed this approach using a model of an optically heterogeneous object represented by a system of concentric spherical layers with the diameters , volumes , and the refractive indices . The section of such an object by the vertical plane is shown in Fig. 1(a). For example, OPD in the point is the sum of contribution of two layers: , where in the geometric path the first index “” corresponds to the position of the point, and the second index “” corresponds to the number of the layer. In general, the value of OPD at the boundary of ’th zones can be represented as the sum: where is its refractivity. For example, at the border of the zone (, ), .The OPD “jumps” in the image plane in the proximity of the points that are the projections of the adjacent layer due to the difference (). The projections of the boundaries of these layers in the topogram [Fig. 1(b)] are represented by concentric circles. The area delimited by the circles with the diameters and , is called henceforth the ’th zone; the change of OPD within a zone is assumed to be smooth, that is, without abrupt changes. The number of the zone corresponds to the index of the external ring diameter. The null zone is bound by the diameters ; its area . The algorithm for determining the boundaries will be discussed in Sec. 3. The boundaries of the zones are represented by the circles, even though they generally form a closed contour [Fig. 1(b), dotted]. Figure 1(c) presents two phase thickness profiles and in the orthogonal cross-sections of the topogram. The points with characteristic values of the phase thickness () and diameter () corresponding to the boundaries of the zones are depicted in one of the profiles. Our concept can be summarized as follows:
In line with this concept, we used the following functions and parameters of the phase image of an individual lymphocyte:
The phase volume of the whole cell is the sum of all segments of the phase volumes , where is geometric volume of the ’th organelle. A part of the phase volume of the cell limited by the horizontal plane and the area in the section of a 3-D image of lymphocyte are shown in Fig. 1(e). A more detailed explanation of these functions and parameters is given in Secs. 3.3 and 3.4. Of note, no cell model, or a priori information about cell structure are needed to determine these parameters, except those mentioned in the items , , and . The algorithm for determining the zone boundaries based on the phase thicknesses values () is important, because their accuracy significantly affects the parameters of the phase portrait. The criterion for the zone boundaries in the phase image of the cell is local change in the slope of the phase thickness profile , and/or the patterns of and (Sec. 3). This algorithm is based on the assertion that the boundary between the zones coincides with the projection of the abrupt change of the refractive index onto the image plane in the optically heterogeneous object. An alternative definition of the zone boundaries can be obtained from the refractivity profile function . However, the necessity to use an optical model seems to be disadvantageous. Several indices were used to determine the boundaries between the zones and their parameters. These indices were provisionally identified with organelles using the notation for the phase thickness profile [Fig. 1(c)]: 0 for the outer boundary of the cell; 1. for the border between the peripheral and dense parts of the cytoplasm; 2. for the boundary of the dense cytoplasm with the chondriome, the giant branched mitochondria; 3. for the border between the chondriome and the nucleus; and 4. for the border between the nucleus and the nucleolus. It is noteworthy that the boundaries of organelles in phase images may or may not coincide with the respective boundaries obtained by other optical methods. The values of the area at the border of the ’th zone were used to calculate the area-equivalent diameter () of the circles instead of the contour . The diameters () of the circles are shown in the topogram [Fig. 1(b)] and in the phase thickness profile [Fig. 1(c)]. An optical model of the cell is required to determine the zones and to calculate the refractivity () at their boundaries. The phase thickness at the ’th boundary of the profile in this model can be represented as a summed contribution of the refractivity and the geometric thickness of the layers according to Eq. (2). The solution of the system of Eq. (2) gives the values of refractivity of the layers. Note that there is no need to define the boundaries of the zones to calculate the refractivity profile in the multilayered spherical model.14 We present the measurements of lymphocytes isolated from peripheral blood of healthy donors. Measurements were performed with a coherent phase microscope Airyscan developed by us. Its design and parameters have been described in detail elsewhere.7,9,10,12. 3.2.Phase Portrait of a Quiescent lymphocyteThe points in phase thickness profiles and obtained in the orthogonal diametrical sections of the topogram include the diameters () and the thicknesses () that we identified with the boundaries of zones. As shown in Fig. 1(c), the shapes of the phase thickness profiles and slightly differ, probably due to certain deviation of the internal structure from perfect central symmetry. The profile can be interpreted as follows. The zero level of the phase thickness () [Fig. 1(c)] corresponding to the equivalent circle of diameter , defines the outer boundary of the cell. As follows from the comparison of Fig. 1(b)–1(c), only a peripheral part of the cytoplasm corresponding to the zone 0 significantly contributes to the phase thickness within the interval . The phase thickness tends to increase toward the cell center as determined by continuously increased contribution of subcellular structures: the boundaries of the zone 1 correspond to thicknesses and diameters . It is reasonable to suppose that the increase of the phase thickness from to is due to the contribution of the layer 1 whose refractivity () is quite high. The latter parameter was stable; we identified the layer 1 as the “dense” part of the cytoplasm. The increase of the phase thickness in the zone 2 () is associated with an additional contribution of the chondriome, in the 3d zone () with contribution of the nucleus, and in the 4th zone (), with the nucleolus. The position of characteristic points in the profiles and depends on the symmetry of the structure, the location of the section and is not necessarily single-valued. More reproducible values of can be obtained from the statistical portrait of the topogram [Fig. 1(b)], e.g., from the graph [Fig. 2(a)]. A slow increase in the phase thickness at the boundary of the cell in a range [Fig. 1(c)] corresponds to the area difference [Fig. 2(a)]. It reflects the contribution of the refractivity () of the zone 0 attributed to the peripheral area of the cytoplasm. The point corresponds to a change of the function slope. Its position can be determined more accurately from the intersection of tangents of the adjacent parts of [Fig. 2(a)]. The ordinate at the intersection of the tangents represents the number of pixels in the cell image. This gives the value of the entire cell area and allows to calculate the diameter of an equivalent circle area. The position of the next characteristic point (), which can be regarded as the inner boundary of the zone 0 separating the periphery and a “dense” part of the cytoplasm, can be identified by the point of changed slope () in the phase thickness profile [Fig. 1(c)]. The position of this point can be obtained with even higher accuracy from the graph of function, assuming that the boundary criterion is based on the location of the characteristic point . The ordinate , the diameter of the circle , and the geometric volume were accepted as the characteristics of the boundary of the zone 1. Apparently, the accuracy of the zone boundaries determined by this method should depend on the difference in the slope tangents on both sides of the characteristic points. Given that the refractivity changes monotonously and the optical homogeneity of the cell increases, the boundary between the zones may become less contrast and finally disappear. The values of the phase thicknesses, the area, the diameter, and the volume of geometric layers at the boundaries of corresponding zones were: , , ; , , , ; , , ; , , ; , ; , , , , . The areas of the zones (rings) were determined as the difference in the area values at boards, : (zone 0); (zone 1); (zone 2); (zone 3); (zone 4). The dependence of the phase volume on the phase thickness for a quiescent lymphocyte is shown in Fig. 2(b). The values of the phase volume at the boundaries of the zones were ; ; ; ; and , respectively. The characteristic points in the and graphs coincided, the algorithms for its determination were also identical. The dense part of the lymphocyte (zones 3–5, ) was of the total area (). Accordingly, for the geometric volume of the whole cell , the index ; for the phase volume , the index . These values correspond to the nucleus-to-cytoplasm ratio frequently used as a quantitative parameter in cell morphology. Comparison of the parts of the phase thickness profile in Fig. 1(c) with the values of the characteristic points in Fig. 2(a) suggests the following interpretation of the zones in the topogram [Fig. 1(b)]. Only a peripheral part of the cytoplasm projects and contributes to the phase thickness of the zone 0 with the area limited by the external (, ) and internal (, ) circles. The part of the cell corresponding to a larger slope in the profile and a higher optical density, is projected onto the area of the 1st zone . As noted above, the zone 1 is primarily attributed to the dense part of the cytoplasm whereas the peripheral part contributes to its phase thickness insignificantly. The parts (organelles) with the total phase thickness within the range provide additive contributions to the zone 2. Finally, the areas of the zones 3 and 4 correspond to the interval . The values of the characteristic points of function [Fig. 2(b)] also represent important parameters of the phase portrait and provide valuable information about the morphology of lymphocytes. A low steepness in the range reflects a weak refractivity of the cytoplasm. A higher steepness within may account for a higher refractivity and large geometric volume of the dense structures such as the chondriome and the nucleus. Figure 2(b) presents the values of the phase volumes at the boundaries of the zones. Accordingly, the difference in these values allows for estimation of the phase volumes of the structures projected onto these zones, : ; ; ; ; . An average value of the refractivity can be obtained from the geometrical volume of the peripheral cytoplasm. These values of the phase volumes could provide indirect estimates of the refractivity of organelles. However, their values depend on the accuracy with which the geometric volume is determined in the optical model of the cell. In Sec. 3.3 we show that more reasonable values of refractivity can be obtained by solving the system of Eq. (2). 3.3.Refractivity of OrganellesThe profiles of refractivity in the cross-sections , calculated with the algorithm presented in Ref. 14 are shown in Fig. 2(c). The difference in the shape and position of the maxima of orthogonal sections made in orthogonal directions is probably a consequence of deviation from central symmetry of the optical density distribution of the lymphocyte. Spatial heterogeneity of the structure of the chondriome can be observed by electron15–17 and fluorescence microscopy. The radii in Fig. 2(c) approximately correspond to the zone borders in Fig. 1(b). The interval of radii and a low refractivity should be identified with the zone 0 and the peripheral cytoplasm. The next interval (, ) can be identified with the dense part of the cytoplasm, whereas the interval (, ) can be identified with the chondriome. In our experiments rotenone, a drug that blocks mitochondrial respiration, caused a marked decrease of refractivity.7,10 Similar changes in the zone 3 can also occur in mitogen-activated cells (see Sec. 3.4). Finally, the zone corresponding to the next interval (, ) can be identified as the nucleus, and the interval (, ) with the nucleolus. These values of refractivity, together with previously determined geometrical volumes of the layers, allow for the estimation of the phase volume of the layers. Thus, it becomes possible to determine a variety of parameters in the phase image of an individual cell. Some of these parameters [] do not depend on the assumptions about the cell shape and other details of the optical model. However, the values of the zone refractivity and the refractivity profiles do depend on the adequacy of the optical model for the concrete cell. Therefore, the refractivity of organelles is probably less reliable because this parameter depends on the adequacy of a multilayered model and deviations from the spherical symmetry. Table 1 shows the parameters of the phase portrait of an intact lymphocyte. The values , , and correspond to the zones, not to the zone boundaries. The variability of the parameters calculated after measurements of up to 50 individual cells was relatively small. Commonly the following limits were determined: maximum phase thickness , the area , the phase volume , and mean values of whole cell refractivity . The variability of other parameters tended to stay within the same limits. Table 1Phase portrait of a quiescent T lymphocyte.
Note: Here and in Table 2 shown are mean values of measurements of 50 cells per each sample. Standard deviations were within 10%. 3.4.Phase Portrait of the Mitogen-Stimulated lymphocyteThe plant mitogen PHA is known to be a polyclonal activator of lymphocytes.15–17 In the initial period of activation (within of exposure) PHA stimulates a transition from the “resting” phase into . One may expect that the processes in PHA-activated lymphocyte can be manifested in the phase images (Fig. 3). The phase thickness significantly decreased (from 220 nm to ) after 1 h with PHA. The shape of the phase thickness profile of the lymphocyte changed [Fig. 3(a)]. The profile became more flat, its outer diameter increased up to mainly due to the peripheral areas of the cytoplasm. These changes involved an increase of both the entire cell area to and the proportion of the peripheral cytoplasm (), which is clear from the graph of the function [Fig. 3(b)]. The total area of the dense part of the cell () increased by 13%. The total area of the cell grew bigger due to the increased cytoplasmic area compared with in the quiescent cell [Fig. 2(b)]. The graph [Fig. 3(c)] reveals a 4-fold (up to ) increase in the phase volume of the peripheral part of the cytoplasm and a decrease down to (by ) of the dense part of the cell. Comparison of refractivity profiles in Figs. 3(c) and 2(c) revealed that PHA-activated lymphocyte was more homogeneous, the optical asymmetry virtually disappeared, the nuclear refractivity decreased, and the nucleolus became indistinguishable. Less contrasting characteristic maxima were shifted to the cell periphery. This is consistent with the collapse of the chondriome into smaller mitochondria and their subsequent distribution to cell periphery.15–17 Decondensation of chromatin15–17 associated with activated gene transcription, can be a major cause for the decreased nuclear phase thickness and refractivity. The values of the nucleus-to-cytoplasm index decreased to versus 75% in intact cells; (versus 60%), (versus 76%). Table 2 shows the characteristic phase image parameters of PHA-activated lymphocytes. Most parameters changed compared with intact cells (Table 1). These changes were most clearly manifested in the areas and refractivity of organelles. Together with our earlier findings,7,10,12 the decrease of the refractivity and the phase thickness in micro-objects exposed to a variety of different inhibitors should be considered fairly universal phenomenon. We attributed some of these effects to the changes in the structure of bound water.7,10,12 Table 2Phase portrait of a PHA-activated T lymphocyte.
4.Concluding RemarksThe unique features of phase images of individual cells allowed to localize characteristic points in and functions and quantitatively determine key physical parameters, that is, phase thickness, phase volume, and area without reference to any optical model. Measurements of these parameters in absolute units, the estimation of OPD with a high precision, the possibility to record dynamic processes, and an enhanced spatial resolution provide valuable information about cellular morphology and physiology. Furthermore, reproducible and reliable results can be obtained with a small number of cells that might be critical for clinical situations.18 Simultaneous measurement of the phase thickness and the absorption coefficient allowed to obtain the exact values of these parameters in a single erythrocyte.18 This and other examples3 to mention a few, show fundamental difference between the interference methods and other techniques based on measurements of physical parameters in relative units. Obtaining statistically meaningful results with these methods frequently makes the experiments laborious and time consuming. Our major result is the representation of an individual living cell as a physical object by a phase portrait with potentially big ( to 30) number of important parameters including the refractivity of organelles. A simplified spherical model was sufficient as a first approximation for identification of the elements of phase images with the subcellular structures of the lymphocyte and for quantitative evaluation of their refractivity. In turn, the parameters of the phase images allow making the optical model of the cell more accurate. To our knowledge, the interference methods have not yet been used to analyze morpho-physiological complexity of nucleated cells. This study is the first to utilize the numerical values of the zone parameters in the phase images for characterization of the functional state of individual lymphocytes. The concept of zones is based on a nonobvious assumption that structural elements of cells are characterized by differential refractivity. Further analysis should be focused on the parameters of boundaries between the zones, their widths and slopes. We believe that the parameters of the phase portraits can be used in the design of new biomedical equipment with a high informative potential and suitability for rapid diagnosis. AcknowledgmentsThe authors thank I. Vasilenko and her colleagues for providing donor lymphocytes, and L. Kovalchuk for critical discussions. This work was supported by Grant No. 2.1.2/12944 of Ministry of Education and Science of the Russian Federation. ReferencesY. Parket al.,
“Metabolic remodeling of the human red blood cell membrane,”
Proc. Natl. Acad. Sci. U. S. A., 107
(4), 1289
–1294
(2010). http://dx.doi.org/10.1073/pnas.0910785107 PNASA6 0027-8424 Google Scholar
Y. Parket al.,
“Static and dynamic light scattering of healthy malaria-parasite invaded red blood cells,”
J. Biomed. Opt., 15
(2), 020506
(2010). http://dx.doi.org/10.1117/1.3369966 JBOPFO 1083-3668 Google Scholar
Y. Parket al.,
“Refractive index maps and membrane dynamics of human red blood cells parasitized by Plasmodium falciparum,”
Proc. Natl. Acad. Sci. U. S. A., 105
(37), 13730
–13753
(2008). http://dx.doi.org/10.1073/pnas.0806100105 PNASA6 0027-8424 Google Scholar
B. Rappazet al.,
“Erythrocyte volume and refractive index measurement with digital holographic microscope,”
Proc. SPIE, 6445 644509
(2007). http://dx.doi.org/10.1117/12.700463 PSISDG 0277-786X Google Scholar
B. Rappazet al.,
“Measurements of the integral refractive index and dynamic cell morphometry of living cells with digital holographic microscopy,”
Opt. Express, 13
(23), 9361
–9373
(2005). http://dx.doi.org/10.1364/OPEX.13.009361 OPEXFF 1094-4087 Google Scholar
G. Popescu,
“Quantitative phase imaging of nanoscale cell structure and dynamics,”
Methods in Cell Biology, 90 87
–115 Elsevier Inc., Amsterdam
(2008). Google Scholar
V. Tychinsky,
“The metabolic component of cellular refractivity and its importance for optical cytometry,”
J. Biophotonics, 2
(8–9), 494
–504
(2009). http://dx.doi.org/10.1002/jbio.v2:8/9 JBOIBX 1864-063X Google Scholar
Z. Wanget al.,
“Tissue refractive index as marker of disease,”
J. Biomed. Opt., 16
(11), 116017
(2011). http://dx.doi.org/10.1117/1.3656732 JBOPFO 1083-3668 Google Scholar
V. P. TychinskyA. N. Tikhonov,
“Interference microscopy in cell biophysics. Principles and methodological aspects of coherent phase microscopy,”
Cell Biochem. Biophys., 58
(3), 107
–116
(2010). http://dx.doi.org/10.1007/s12013-010-9114-z CBBIFV 1085-9195 Google Scholar
V. P. TychinskyA. N. Tikhonov,
“Visualization of individual cells and energy-transducing organelles,”
Cell Biochem. Biophys., 58
(3), 117
–128
(2010). http://dx.doi.org/10.1007/s12013-010-9115-y CBBIFV 1085-9195 Google Scholar
Z. Wanget al.,
“Spatial light interference microscopy (SLIM),”
Opt. Express, 19
(2), 1016
–1026
(2011). http://dx.doi.org/10.1364/OE.19.001016 OPEXFF 1094-4087 Google Scholar
V. P. Tychinsky,
“Dynamic phase microscopy: is a “dialog” with the cell possible,”
Physics-Uspekhi, 50
(5), 513
–528
(2007). http://dx.doi.org/10.1070/PU2007v050n05ABEH006222 PHUSEY 1063-7869 Google Scholar
V. P. Tychinsky,
“Super-resolution and singularities in phase images,”
Phys. Uspekhi, 51
(11), 1161
–1169
(2008). http://dx.doi.org/10.1070/PU2008v051n11ABEH006682 PHUSEY 1063-7869 Google Scholar
Y. HiramotoS. ShumeiS. Yoko,
“Refractive index of the protoplasm in sea urchin eggs,”
Develop. Growth Differ., 21
(2), 141
–153
(1979). http://dx.doi.org/10.1111/dgd.1979.21.issue-2 DGDFA5 0012-1592 Google Scholar
V. ManteifelL. BakeevaT. Karu,
“Ultrastructural changes in chondriome of human lymphocytes after irradiation with He-Ne laser: appearance of giant mitochondria,”
J. Photochem. Photobiol., 38
(1), 25
–30
(1997). http://dx.doi.org/10.1016/S1011-1344(96)07426-X JPPBEG 1011-1344 Google Scholar
V. M. ManteifelT. I. Karu,
“Loosening of condensed chromatin in human blood lymphocytes exposed to irradiation with a low-energy He-Ne laser,”
Izv. Akad. Nauk Ser. Biol., 36
(6), 654
–661
(2009). http://dx.doi.org/10.1134/S1062359009060028 IRABEC 1026-3470 Google Scholar
V. ManteifelT. Karu,
“Ultrastructural changes in chondriome of human lymphocytes after irradiation with He-Ne laser: appearance of giant mitochondria,”
J. Photochem. Photobiol., 38 25
–30
(1997). http://dx.doi.org/10.1016/S1011-1344(96)07426-X JPPBEG 1011-1344 Google Scholar
M. MirK. TangellaG. Popescu,
“Blood testing at the single cell level using quantitative phase and amplitude microscopy,”
Biomed. Opt. Express, 2
(12), 3259
–3266
(2011). http://dx.doi.org/10.1364/BOE.2.003259 BOEICL 2156-7085 Google Scholar
|