Angularly resolved elastic light scattering is an established technique for probing the average size of organelles in biological tissue and cellular ensembles. Focusing of the incident light to illuminate no more than one cell at a time restricts the minimum forward-scattering angle θmin that can be detected. Series of simulated single-cell angular-scattering patterns have been generated to explore how size estimates vary as a function of θmin. At a setting of θmin = 20 deg, the size estimates hop unstably between multiple minima in the solution space as simulated noise (mimicking experimentally observed levels) is varied. As θmin is reduced from 20 deg to 10 deg, the instability vanishes, and the variance of estimates near the correct answer also decreases. The simulations thus suggest that robust Mie theory fits to single-cell scattering at 785 nm excitation require measurements down to at least 15 deg. Notably, no such instability was observed at θmin = 20 deg for narrow bead distributions. Accurate sizing of traditional calibration beads is, therefore, insufficient proof that an angular-scattering system is capable of robust analysis of single cells. Experimental support for the simulation results is also presented using measurements on cells fixed with formaldehyde. |
1.IntroductionLight scattering measurement techniques are useful for studying dynamic biological processes because they are nondestructive, noninvasive, and do not require exogenous labels. Angular elastic scattering, one such technique, can be used in order to obtain morphological information from biological samples, such as the mean size of cell organelles, which can yield important information about the health and behavior of the cells. Angular domain elastic scattering of cell suspensions and tissues has been used to diagnose precancerous conditions and observe subcellular changes in organelles during apoptosis.1–4 The angular scattering properties of the sample depend upon parameters such as size and refractive index. By comparing angular scattering patterns to a theoretical model, such scatterer characteristics can be determined. Mie theory5 is a commonly used elastic scattering model for spherical scatterers with a diameter on the order of the incident wavelength. The scattered intensity depends on the scatterer size, refractive index of both the scatterer and the surrounding medium, the incident wavelength, the polarization of the incident light, and the scattering angle. For angular scattering, the results of Mie theory are particularly sensitive to small changes in scatterer size.6 In the 1980s, generalized Lorenz-Mie theory (GLMT) was developed to describe elastic scattering using incident waves that are Gaussian rather than planar.7–10 Scattering models other than Mie theory can be used to analyze angular scattering, including the T-matrix method11 and complex focus theory.12 Angularly resolved forward scattering over broad ranges (from a few degrees to as high as 90 deg) with fine resolution (on the scale of 1 deg) has been modeled using Mie theory to estimate organelle size distributions in ensembles of cells. Mourant et al. calculated scatterer size distribution ranging from tens of nm to several microns in epithelial cells. Wilson et al.13 used this approach to detect size distribution shifts associated with photodynamic therapy that targeted specific organelles (mitochondria and lysosomes, separately).14,15 Both of these groups used scanning goniometers to obtain their angularly resolved data. The motivation of the work described below is to extend this wide-angular-range, high-angular-resolution, Mie-based approach to the analysis of single cells. Forward-directed angular scattering is collected using an angular scattering microscope design inspired by Cottrell et al.,16 rather than goniometric detection. Detecting a two-dimensional scattering pattern with fine angular resolution over a large angular range can be accomplished quickly, and it is a more light-efficient method of collecting the scattering from single cells. Because we also aim to measure the scattering from a single cell that may have a neighbor, it is necessary to focus the excitation beam onto the sample plane. Angular scattering measurements of single cells over time would allow for the determination of differences in structure and behavior between individual cells. For the work presented in this paper, the excitation beam is focused to a subcellular spot size to ensure that only a single cell is illuminated. Due to the focusing of the beam to a spot, a nontrivial portion of forward-scattering solid angle is obscured by unscattered light. This contrasts with studies that excite with a more collimated beam and correspondingly illuminate a larger region. As such, a trade-off exits between the spot diameter and the measurable scattering angles. For example, a Gaussian spot size of diameter (defined using the intensity width) using 785-nm light restricts angular measurements to 11 deg and higher. To have measurements down to 1 deg unobscured by unscattered light, the spot diameter has to increase to . In addition to the limit imposed by the focused beam size, other practical effects further confine the range of analyzable scattering angles. If a physical mask is placed in the Fourier plane to prevent unscattered light from reaching the detector, flare around the edges of the mask, which can have sharp edges with some irregularities, can obscure additional angles. Furthermore, measuring the scattering from single cells introduces several challenges. These include lower scattering intensities that must compete with the system background noise and speckle with larger grain size which can influence the fits to the data. Although the effect of the speckle grain can be reduced using a temporally incoherent excitation beam, broadband illumination washes out the Mie scattering features in the angular plane.17 In this paper, we present a study of the effect of angular range on the stability and accuracy of Mie theory fits to scattering from organelle-sized beads and from single cells. Simulated data permit investigation of behavior at the noise-free limit as well as at signal-to-noise levels typical of single-cell measurements on our current apparatus. Although several studies in the literature have extracted organelle size parameters of multiple-cell samples, to our knowledge no group to date has assigned uncertainties to size predictions derived from the angular scattering from a broad distribution of sizes similar to that of cell organelles. Here we also investigate the effect that angular range has on the uncertainty using computational methods. Experimental data from formalin-fixed single mammalian cells are also presented. 2.Methods2.1.SimulationsSimulations were performed to test the effect of the detected angular range on the accuracy and precision of the fits to scattering from cell organelles. Scattering from an ensemble of spheres was simulated using GLMT using a linearly polarized excitation beam with a diameter of and a simple addition of intensities. This method ignores any interference effects between spheres, which is a common approximation in the literature.14,18–20 Three different size distributions were simulated: a narrow log-normal distribution of polystyrene beads with a mean diameter of and a standard deviation of , a single bead with a diameter of , and a broad log-normal distribution of cell organelles with a mean diameter of and a standard deviation of . The first two distributions correspond to polystyrene beads that are used for routine calibration. The final distribution is a mitochondria-dominated distribution derived by Wilson et al.15 from angular scattering analysis of squamous cell carcinoma suspensions. For this study, we assume a simplified model of the cell, in which the scattered light originates from a single population of scatterers (mitochondria, reported as the dominant light scatterers in the 10-deg to 60-deg forward-scattering range in some mammalian cell experiments20) so that the effect of minimum scattering angle alone could be isolated. (Comments about multiple population distributions are offered in Sec. 4.) For each case, noise was added to the simulated scattering patterns that mimicked the observed noise in experimentally acquired scattering data. For the two bead cases, noise proportional to the square root of signal was used and scaled empirically to match the percentage seen in our experimental bead measurements (taken with a 50-ms exposure time and power of about 20 mW at the sample) to resemble shot noise. For single cell measurements, the small excitation beam diameter and multiple scatterers (organelles) in the sample result in speckle grains larger than single pixels of the detector array. Therefore, for the mitochondrial distribution, speckle noise was simulated by first adding white noise proportional to the signal square root, and then low-pass filtering the data. This resulted in scattering patterns with noise resembling speckle grain sizes and amplitudes observed in our single-cell experiments. Size distributions for each case and the resulting simulated scattering pattern with noise are shown in Fig. 1. The simulated data were then modeled using Mie theory and a single-log-normal size distribution with mean size and width distribution . As described in Sec. 2.3, a least square error metric function given by Eq. (3) was minimized to identify the parameter pair for the best theoretical fit. Predictions were performed 100 times with different simulated noises. The uncertainties in and in were defined as the standard deviation of the predictions from the 100 noise simulations. The standard error of each uncertainty was calculated by21 where .To test the influence of angular range, size estimates were obtained as a function of the minimum detection angle , ranging from 5 to 20 deg. The maximum angle was held constant at 60 deg. Additional simulations were performed, in which was fixed at 10 deg and the effect of decreasing the maximum angle was studied. Finally, in order to decouple the effects of angular range and speckle noise on the uncertainty in the size estimates, simulations were performed, in which the speckle grain size was also varied to more closely simulate the scattering from multiple cells at once. This allows us to compare the uncertainty in the estimates not only for various angular ranges, but also for single-cell and multicell scattering measurements. 2.2.Experiments2.2.1.Optical systemAngular scattering data were acquired as described previously.22 The angular scattering microscope diagram is shown in Fig. 2. Briefly, a collimated 785-nm diode laser source (I0785SU0050PA-TK, Innovative Photonic Solutions, Monmouth Junction, New Jersey) is linearly polarized and then directed through a condenser lens that produces a spot size with a focusing half-angle of and a power of at the sample. Elastically scattered light is collected by a 1.3 NA oil-immersion microscope objective (-PFO, Nikon Instruments, Tokyo, Japan). A 4f lens system images the back aperture of the objective onto a glass coverslip with an evaporated gold disc that rejects scattering angles of , thus preventing unscattered light from reaching the detector. A second 4f system then relays this Fourier image, referred to here as a scattergram, to a CCD array (, pixel pitch, Luca-S DL-658M-TIL, Andor Technology, Belfast, Northern Ireland). Images of the scattering pattern were acquired by summing five images with 10-ms exposure times. 2.2.2.Bead sample preparationBead samples were prepared for measurements by diluting one drop from the tip of the vial of polystyrene beads from a known size distribution in 10 mL of water. Approximately 0.5 mL of the bead dilution was placed on a coverslip in a sample chamber. The water was allowed to evaporate overnight, which adhered the beads to the glass coverslip due to electrostatic interactions. For measurements, the sample was rehydrated. 2.2.3.Cell culture and preparationSCC7 mouse squamous cell carcinoma (ATCC, Manassas, Virginia, courtesy of the Foster Lab, University of Rochester Medical Center) were grown in an adherent monolayer in a Petri dish designed to promote cell adherence (Falcon standard tissue culture dishes, 08-772E, Fisher Scientific) and containing RPMI 1640 medium (Gibco, ThermoFisher Scientific) with 10% fetal bovine serum (Gibco) and antibacterial/antifungal agents (Plasmocin and Primocin, InvivoGen, San Diego, California). Cells were lifted from the dish using 0.25% trypsin (Gibco) and seeded into a new dish with round glass coverslips (12-546-1, Fisherbrand, Fisher Scientific, Pittsburgh, Pennsylvania) and media. The cells adhered to the coverslips, which were then removed from the dish and placed into an Attofluor cell chamber (A7816, Invitrogen, Carlsbad, California). For fixation, the coverslip in the cell chamber was washed with phosphate-buffered saline (PBS) and the cells were then fixed by filling the chamber with 4% formaldehyde for 15 min at room temperature. The formaldehyde was aspirated and the cells were washed twice with PBS. The chamber was then filled with phenol red-free RPMI 1640 cell culture medium for cell measurements. The chamber was placed in the sample plane of the system in a stage-top incubator (Stable “Z,” Bioptechs, Butler, Pennsylvania) which provided heat from the objective lens and the perimeter of the chamber and allowed for the delivery of humidified 5% to the cells. Although this study used only fixed cells, the temperature, humidity, and settings were the same as for live cells in order to replicate all sources of uncontrolled variation. Prior comparisons between fixed and unfixed SCC7 indicated that the fixation process had little to no influence upon the measured scattergram. In this study, angular scattergrams were measured from one fixed cell repeatedly for 35 min. 2.3.Data FittingEach two-dimensional angular scattering pattern contained scattering signal as a function of polar angle (where represents forward-directed unscattered light) and azimuthal angle . The data in each pattern were azimuthally averaged to produce a single vector of scattered intensity versus polar angle. The maximum analyzable angle (60 deg) is determined by the numerical aperture of our microscope objective. The minimum analyzable angle () will be investigated here. A theoretical model then determined the size distribution of Mie-scattering spheres that best explained the data. The same approach was used for fitting simulated and experimental data. First, a lookup table of azimuthally averaged scattered intensity versus and sphere diameter was calculated using GLMT. ranged from 5 to 60 deg in steps of 0.275 deg, and from 10 nm to in steps of 10 nm. Then inspired by the work of Wilson et al.,13 the ensemble of scatters was modeled as a single log-normal size distribution defined by a mean , a standard deviation , and an overall scaling parameter (unitless, proportional to the total number of scatterers). Defining as the scattering dependence for , the best scaled fit for a given pair was determined by linear least-squares optimization of where is a linear offset to compensate for experimental baseline drift.The MATLAB downhill simplex function fminsearch was then used to estimate the best pair by minimizing a -like error metric given by Because angular scattering is much stronger at smaller values of , the second sum was included to boost the influence at the larger angles. The relative weight of 10 for this sum was selected empirically. Estimates were forbidden which resulted in the ratio of tail of the cross-section weighted size distribution to the peak being larger than a value of 0.1. 3.Results3.1.Simulation Results3.1.1.Narrow bead distributionsThe uncertainty in the predicted mean diameter for the narrow simulated bead distribution and the simulated single five micron bead was determined as a function of between 5 deg and 20 deg over 100 iterations of speckle noise. At the various angles, the uncertainty ranged from 4 to 17 nm for bead fits, and 6 to 13 nm for scattering, in neither case with any pattern versus angle. These uncertainties are on the order of our Mie theory lookup table diameter step size of 10 nm. Figure 3 shows representative maps of [the inverse of the error metric in Eq. (3)] as a function of the two log-normal size population parameters and . is normalized in each plot to a maximum of 1, which corresponds to the best fit (smallest value of .) Maps are shown for the two simulated bead diameters and for values of both 20 deg and 10 deg. Because of physical assertions about the size distributions (e.g., requiring the distribution weight to be close to zero at ), a zone of forbidden predictions exists near the bottom of the plots (shown in dark blue). The error maps have step sizes of 100 nm in the dimension; on this scale, the best fit is confined to a single pixel. As noted above, the fluctuation in the selected value over the 100 simulations is on the scale of 10 nm. The error maps also indicate that in three of the four cases, the fit correctly estimates the population width to be as close to zero as possible. The only exception occurs at , where the best fit for the beads is . Switching to eliminates this error. Overall, the error metric maps for the two different minimum angles are similar for both the large and small bead diameter cases. This further supports the argument that does not have a significant impact on the uncertainty of the fits to narrow bead distributions. 3.1.2.Mitochondria distributionWhen the population width is broader, angular scattering becomes less sensitive to changes in . Although the effect is clear when the narrow-width pattern has multiple peaks, e.g., see Figure 11.6 in Ref. 6), it is less obvious for monotonically decreasing scattering. To illustrate the effect, Mie theory plots for a bead and for the mitochondria-like size distribution ( and ) are shown in Fig. 4. In each plot, the scattering for a 100-nm increase in is also shown, with the two curves scaled to match at . For the bead, a 100-nm change in diameter causes a visible change in the angular dependence both above and below 20 deg, which explains why the parameter maps in Fig. 3 have sharp maxima regardless of . For the mitochondrial distribution, however, the scattering changes are much smaller. The inset, showing the absolute value of the residual, shows that changes are now more concentrated at angles below 20 deg. This suggests that, unlike with narrow bead distributions, using instead of 20 deg should confer greater sensitivity to . When maps similar to those in Fig. 3 are generated for the simulated mitochondrial population ( and ) show that has a significant effect on the stability of the size predictions. Example error metric maps for simulated single cell scattering are shown in Fig. 5. In both cases the value of has a local maximum near and . For , there are also local maxima for near 3 and with comparable values. When the simulated noise was varied 100 times, the global best fit lay in these incorrect lobes of solution space 40% of the time. The extra lobes are absent for , as shown in this figure, and the primary lobe shrinks. We found that values of of were sufficient to prevent global maxima occurring outside of the primary lobe. To emphasize the greater robustness at lower angles, Fig. 5(c) (bottom) shows two fits corresponding to the local maxima in space marked by asterisks in panel a. Both fits follow the data to within the noise. When those two fits are extended to in panel d, the fit from the incorrect lobe (asterisk 2) deviates drastically from the simulated data at the added angles. Even if solutions from the other lobes are rejected (e.g., by restricting the range of accepted values), Fig. 5’s error maps also reveal that is more sharply peaked around its maximum for . The dashed line in (a) traces the contour of fits with , i.e., a moderately good fit relative to the maximum of 1.0. This same contour line is overlaid on the error map, showing that the new area of fits with is now smaller. For these simulated cell distributions, the uncertainty of the predicted and as a function of was determined by calculating the standard deviation of the predictions as the speckle noise was simulated 100 times, and is plotted in Fig. 6. The uncertainty in the predicted mean drops as is decreased, until it plateaus around 14 deg. The uncertainty in the predicted population width continues to decrease as is decreased lower than 14 deg (data not shown). These results illustrate the importance of collecting angles below 20 deg when performing fits, in which and are both free parameters. At the other end of the angular range, varying from 30 deg to 65 deg produced no significant change in the uncertainty of the predictions for . Note that the error maps shown in Fig. 5 are idiosyncratic to the particular speckle noise generated for one simulation. At , according to the upper uncertainty plot of Fig. 6, another noise iteration would produce a map where the best fit has a value up to different from the previous one. For the error map of Fig. 5, points that are horizontally displaced from the best fit (asterisk 1) have values close to 0.7. This suggests a new way to read the error map: that all points with values of 0.7 are pairs that could be optimal for another instance of speckle noise on the same underlying Mie pattern. We will return to this line of reasoning when discussing experimental measurements of a fixed cell. These simulations were repeated using a smaller speckle grain, as would occur if there were contributions over a larger field of view from multiple cells. The results were qualitatively the same: secondary lobes appeared in the map for , and this phenomenon ceases to exist at . Additionally, as shown in Fig. 6, the uncertainty in significantly improved versus reduction in until around 15 deg. Unsurprisingly, the absolute uncertainty at low angles is lower for the smaller speckle grain size. From these results, we can conclude that measuring lower scattering angles is crucial not only for single-cell measurements, but also for multicell measurements. 3.2.Experimental Results3.2.1.Polystyrene bead measurementsAs expected from our simulations of beads, theory fits to experimental scattering measurements from single beads did not result in striking differences in the predicted diameter as a function of . Representative plots for one and one diameter bead are shown in Fig. 7. Plausible size predictions are obtained when the beads’ azimuthally averaged scattering patterns are fit over the angular range of 20 deg to 60 deg [Fig. 7(a)]. Figure 7(b) shows similar fits and size predictions when fits are performed using . Estimates of for 20 single- beads are summarized by box plots in Fig. 8. Angular scattering from 20 different single- beads was fit using values of 10 deg and 20 deg. The average predicted did not differ significantly when was reduced from 20 deg to 10 deg. However, the of the distribution was predicted more accurately using . Because only a single bead at a time was measured, the predicted for each fit should be close to zero. Because our fitting algorithm cannot choose a width of zero, the minimum possible width prediction was the 10-nm step size of the theory lookup table. A Wilcoxon rank-sum test showed the distribution of values for the analysis to be significantly different (), and according to Fig. 8 clearly larger, which is less accurate. It is worth noting that the single beads were sized accurately only if the width of the size distribution was forced to be zero. When the width was a free parameter, regardless of , the fits resulted in predicted means that were smaller than the manufacturer specifications (typically about 0.8 to ) and large distribution widths (100 to 200 nm). This is most likely because of lower signal-to-noise due to the reduced scattering signal from the smaller diameter bead. 3.2.2.Single-cell measurementsThe experimental scattering pattern, azimuthal average, and resulting best Mie theory fit from a single fixed cell are shown in Fig. 9. The theory curve is fit to the data to within the noise, whose amplitude is dictated by the large speckle grains. Figure 10 shows the corresponding error metric map for a typical cell as a function of size distribution mean and width for and 10 deg, similar to those shown in the simulation results. From the experimental error maps, we see that the best fits cluster at population widths that are noticeably larger than 0. This effect is more pronounced for . No single-scatterer size is able to model the intensity falloff of the single-cell data over the entire angular fitting range. Our single-cell scattering data are only well to within the noise when a broad distribution of sizes is modeled, as one would expect a priori. The experimental single-cell map for (Fig. 10, left) exhibits multiple local maxima in the map. This is evocative of the results of the simulated organelle scattering (cf., Fig. 5). The map has only one local maximum lobe, again in parallel with simulations. Recall from the simulations that the true value for a single-log-normal population could be assigned a value as low as 0.7 due to speckle noise. The region of values greater than 0.7 for the experimental results (outlined in the plot) contains estimates with values in the range 0.8 to ; two such points in that regime are indicated by asterisks. In addition, the estimates of for the fixed SCC7 cell fluctuated less over 35 min even when the solution space was confined near the correct lobe. Using , the estimated drifted slowly between 300 and 600 nm over the 35 min. When analysis included angles down to 10 deg, estimates stayed steady to within 30 nm. 4.Discussion4.1.BeadsMeasurements of single beads or bead ensembles with known size distributions are commonly used to calibrate an angular scattering system and ensure that good quality data and fits are being obtained.16,17,22,23 For single beads or narrow size distributions of scatterers, analyzing the scattering down to only 20 deg is sufficient for estimating mean diameter. In our simulations, changing to did not significantly change the uncertainty in the predicted sizes for either large () or small () diameter beads (cf., Fig. 3). This matches our experimental experience with bead measurements (cf., Fig. 7). Even for single beads, however, there was a benefit to incorporating lower angles. Fits to single- bead scattering using consistently returned estimates near zero as was correct, while analysis using did not (cf., Fig. 8). 4.2.CellsThe results above for beads invite the thought that might produce similarly acceptable results for cell organelle analysis. As noted by Bohren and Huffman,6 however, the heterogeneity in broader distributions blurs the features in the scattering patterns. When we fit the simulated scattering from cell organelle-like size distributions, results for were markedly inferior in two ways to those obtained at . First, there were multiple local maxima in solution space for , and realistic levels of measurement noise were sufficient to vary which location was the global maximum; this effect vanished for . Second, in the neighborhood of the true values, the map was less sharply peaked for than for . For our empirically observed level of speckle noise, there appears to be a crucial turning point in the cell scattering curve at , as is shown in Fig. 6. As a result, the fits to these lower angles are not simply more accurate because more scattering angles are included, but because there is a shape change in the scattering versus angle curve specifically at the low angles that must be included to obtain accurate and stable fits. This behavior was not seen in the size predictions from narrower size distributions, as the effect ultimately arises from the broad distribution of organelle sizes in the cell. Experimental measurement specifications that are sufficient for sizing narrow bead distributions may be insufficient for analyzing broader distributions. This highlights that validating an angular scattering system using narrow bead distributions, a common step in the angular scattering literature, is necessary but not sufficient to ensure that the system will work for cells, i.e., for broad distributions of organelle sizes. The results indicate that speckle grain size influences the prediction accuracy. Single cells have a maximum distance between scatterers dictated by the size of the cell, leading to larger speckle grains than if multiple cells contribute to the measurement. In turn, the simulation results in Fig. 6 show that the larger speckle grains associated with single cells lead to a larger uncertainty in estimates of , even when the lowest angles are incorporated. For single cells, the simulation reported an uncertainty of ; for a grain size approximately five times smaller, it dropped to . The likely reason is that a larger speckle grain pulls a larger amount of the lowest-angle data in the same direction, and perturbations at the lowest analyzed angles affect the fit parameters most radically. Even when is used, the single-cell angular scattering data given here can be modeled by Mie theory to within the noise using a single-log-normal size distribution. Previous work by Mourant et al.13 and Wilson and Foster20 using goniometers has modeled angular scattering from cell suspensions to two or more populations of scatterers, and Wilson et al. specifically illustrated that a single population was inadequate to fit their angular data. There are at least two differences in the cell measurements performed here. First, our microscope-based instrument measured down to 10 deg, as compared to 3 deg in Ref. 18. Had our single cells been measured down to 3 deg, a single-population fit might not have been adequate. Second, the speckle grains in a single cell measurement distort the data more. They are larger both in angular extent (due to the reduced area over which the scattering is emitted, as discussed above) and in amplitude (because there is no averaging from multiple cells or from cell motion in suspension). Efforts are underway in our group to mitigate both the angular limitations and the speckle effects. Compared to the simulation results in Fig. 5(b), the region of best fits () in the SCC7 cell’s error map [Fig. 10(b)] was pushed to lower and values overall. The direction of this shift is unsurprising, given that the simulation contained solely a mitochondrially inspired population distribution with a mean of (from Ref. 15), while the experimental measurement of an intact cell included smaller organelles. The best-fit region contained values as high as , which is a plausible intermediate value between mitochondria and smaller organelles, given that the current fits are limited to a single population as discussed in the previous paragraph. In addition to reducing the speckle amplitude and the minimum measured scattering angle, more work is needed to establish both (a) a range of estimates obtained for single SCC7 cells by angular scattering and (b) reference information about the actual organelle size distribution within such cells. This study showed that the size estimates for a fixed cell were dramatically more stable over 35 min when angles below 20 deg were included in the analysis. For , the stability was at the 30-nm level. This does not mean that the estimates were accurate. We believe that when angular scattering is performed on single adherent cells using , the limit on the accuracy is imposed by the speckle. In this case, because the cell was fixed, the speckle remained largely the same throughout the 35 min, hence the best Mie fit kept nearly the same parameters. For the case, the measurement-to-measurement variations due to shot noise and system drift were enough to make the Mie fits unstable. With the added angular information between 10 deg and 20 deg, the data analysis stabilized, but we emphasize that reducing the speckle in single-cell measurements is the essential next step for replicating the accuracy of multicell organelle sizing analysis in single cells using angular scattering. As previously mentioned, the aim of this work is to extend wide-range/high-resolution angular scattering plus Mie analysis to the single-cell regime to obtain scatterer size distributions. There are other methods of deriving size distributions in single cells from variations in angular scattering. Zhang and Boustany24 used optical scatter image ratio (OSIR) to assign scatterer sizes to subdomains of a cell image and thus obtain scatterer size distribution estimates for the cell as a whole. In that work, two wide-field images of cells were acquired at different numerical apertures. A two-element vector of angular information is thus built up at each pixel, sufficient for Mie theory to estimate a sphere diameter at each location, from which the size distribution emerges. Scatterer motion and objects larger than the pixel size introduce error into the estimation process. Here, by comparison, a large number () of angular values is obtained from the illuminated area. A benefit is that this enables multiparameter fitting to the angular data, extracting not just a sphere size but other moments of the size distribution. On the other hand, the multiparameter fit must assume a functional form (here, a log-normal distribution) whose validity will vary by case, especially for cells with smaller numbers of significant scatterers. Another downside is that all spatial information within the cell is lost by the nature of the Fourier-domain intensity measurement. We note that angular scattering from single transilluminated cells can also be obtained by Fourier-transform of a complex-field image.25–27 To our knowledge, however, inversions of such angular data by Mie theory to obtain size distributions of nonnuclear organelle populations have not been reported. DisclosuresThe authors have no relevant financial interests in this article and no potential conflicts of interest to disclose. AcknowledgmentsThe authors thank Tim Baran and Thomas Foster for helpful conversations and support with cell culturing. This work was partially supported by NSF CBET award No. 1402345 (Biophotonics). ReferencesL. T. Perelman et al.,
“Observation of periodic fine structure in reflectance from biological tissue: a new technique for measuring nuclear size distribution,”
Phys. Rev. Lett., 80
(3), 627
–630
(1998). https://doi.org/10.1103/PhysRevLett.80.627 PRLTAO 0031-9007 Google Scholar
V. Backman et al.,
“Polarized light scattering spectroscopy for quantitative measurement of epithelial cellular structures in situ,”
IEEE J. Sel. Top. Quantum Electron., 5
(4), 1019
–1026
(1999). https://doi.org/10.1109/2944.796325 IJSQEN 1077-260X Google Scholar
Y. L. Kim et al.,
“Simultaneous measurement of angular and spectral properties of light scattering for characterization of tissue microarchitecture and its alteration in early precancer,”
J. Sel. Top. Quantum Electron., 9
(2), 243
–256
(2003). https://doi.org/10.1109/JSTQE.2003.814183 Google Scholar
N. N. Boustany et al.,
“BCL-xL-dependent light scattering by apoptotic cells,”
Biophys. J., 87
(6), 4163
–4171
(2004). https://doi.org/10.1529/biophysj.104.048736 BIOJAU 0006-3495 Google Scholar
G. Mie,
“Beiträge zur optik trüber medien, speziell kolloidaler metallösungen,”
Annalen der Physik, 330
(3), 377
–445
(1908). https://doi.org/10.1002/(ISSN)1521-3889 ANPYA2 0003-3804 Google Scholar
C. F. Bohren and D. R. Huffman, Absorption and Scattering by Small Particles, Wiley Interscience, New York
(1983). Google Scholar
G. Gouesbet, G. Grehan and B. Maheu,
“Scattering of a Gaussian beam by a Mie scatter center using a Bromwich formalism,”
J. Opt., 16
(2), 83
–93
(1985). https://doi.org/10.1088/0150-536X/16/2/004 JOOPDB 0150-536X Google Scholar
G. Gouesbet, B. Maheu and G. Grehan,
“The order of approximation in a theory of the scattering of a Gaussian beam by a Mie scatter center,”
J. Opt., 16
(5), 239
–247
(1985). https://doi.org/10.1088/0150-536X/16/5/007 JOOPDB 0150-536X Google Scholar
G. Grehan, B. Maheu and G. Gouesbet,
“Scattering of laser beams by Mie scatter centers: numerical results using a localized approximation,”
Appl. Opt., 25
(19), 3539
–3548
(1986). https://doi.org/10.1364/AO.25.003539 APOPAI 0003-6935 Google Scholar
G. Gouesbet, G. Grehan and B. Maheu,
“Localized interpretation to compute all the coefficients GMN in the generalized Lorenz-Mie theory,”
J. Opt. Soc. Am. A, 7
(6), 998
–1007
(1990). https://doi.org/10.1364/JOSAA.7.000998 JOAOD6 0740-3232 Google Scholar
M. G. Giacomelli et al.,
“Application of the T-matrix method to determine the structure of spheroidal cell nuclei with angle-resolved light scattering,”
Opt. Lett., 33
(21), 2452
–2454
(2008). https://doi.org/10.1364/OL.33.002452 OPLEDP 0146-9592 Google Scholar
N. J. Moore and M. A. Alonso,
“Closed form formula for Mie scattering of nonparaxial analogues of Gaussian beams,”
Opt. Express, 16
(8), 5926
–5933
(2008). https://doi.org/10.1364/OE.16.005926 OPEXFF 1094-4087 Google Scholar
J. R. Mourant et al.,
“Polarized angular dependent spectroscopy of epithelial cells and epithelial cell nuclei to determine the size scale of scattering structures,”
J. Biomed. Opt., 7
(3), 378
–387
(2002). https://doi.org/10.1117/1.1483317 JBOPFO 1083-3668 Google Scholar
J. D. Wilson et al.,
“Light scattering from intact cells reports oxidative-stress-induced mitochondrial swelling,”
Biophys. J., 88 2929
–2938
(2005). https://doi.org/10.1529/biophysj.104.054528 BIOJAU 0006-3495 Google Scholar
J. D. Wilson, W. J. Cottrell and T. H. Foster,
“Index-of-refraction-dependent subcellular light scattering observed with organelle-specific dyes,”
J. Biomed. Opt., 12
(1), 014010
(2007). https://doi.org/10.1117/1.2437765 JBOPFO 1083-3668 Google Scholar
W. J. Cottrell, J. D. Wilson and T. H. Foster,
“Microscope enabling multimodality imaging, angle-resolved scattering, and scattering spectroscopy,”
Opt. Lett., 32
(16), 2348
–2350
(2007). https://doi.org/10.1364/OL.32.002348 OPLEDP 0146-9592 Google Scholar
D. W. Shipp,
“Holographic angular domain elastic scattering of single biological cells,”
Institute of Optics, University of Rochester,
(2014). Google Scholar
J. R. Mourant et al.,
“Mechanisms of light scattering from biological cells relevant to noninvasive optical-tissue diagnostics,”
Appl. Opt., 37
(16), 3586
–3593
(1998). https://doi.org/10.1364/AO.37.003586 APOPAI 0003-6935 Google Scholar
J. R. Mourant et al.,
“Light scattering from cells: the contribution of the nucleus and the efects of prolierative status,”
J. Biomed. Opt., 5
(2), 131
–137
(2000). https://doi.org/10.1117/1.429979 JBOPFO 1083-3668 Google Scholar
J. D. Wilson and T. H. Foster,
“Mie theory interpretations of light scattering from intact cells,”
Opt. Lett., 30
(18), 2442
–2444
(2005). https://doi.org/10.1364/OL.30.002442 OPLEDP 0146-9592 Google Scholar
B. Harding, C. Tremblay and D. Cousineau,
“Standard errors: a review and evaluation of standard error estimators using Monte Carlo simulations,”
Quant. Methods Psychol., 10
(2), 107
–123
(2014). https://doi.org/10.20982/tqmp.10.2.p107 Google Scholar
Z. J. Smith and A. J. Berger,
“Validation of an integrated Raman- and angular-scattering microscopy system on heterogeneous bead mixtures and single human immune cells,”
Appl. Opt., 48
(10), D109
–D120
(2009). https://doi.org/10.1364/AO.48.00D109 APOPAI 0003-6935 Google Scholar
J. W. Pyhtila, R. N. Graf and A. Wax,
“Determining nuclear morphology using an improved angle-resolved low coherence interferometry system,”
Opt. Express, 11
(25), 3473
–3484
(2003). https://doi.org/10.1364/OE.11.003473 OPEXFF 1094-4087 Google Scholar
J. Zheng and N. N. Boustany,
“Alterations in the characteristic size distributions of subcellular scatterers at the onset of apoptosis: effect of Bcl-xL and Bax/Bak,”
J. Biomed. Opt., 15
(4), 045002
(2010). https://doi.org/10.1117/1.3462933 JBOPFO 1083-3668 Google Scholar
M. Kalashnikov et al.,
“Assessing light scattering of intracellular organelles in single intact living cells,”
Opt. Express, 17
(22), 19674
–19681
(2009). https://doi.org/10.1364/OE.17.019674 OPEXFF 1094-4087 Google Scholar
H. Ding et al.,
“Fourier transform light scattering of inhomogeneous and dynamic structures,”
Phys. Rev. Lett., 101
(23), 238102
(2008). https://doi.org/10.1103/PhysRevLett.101.238102 PRLTAO 0031-9007 Google Scholar
W. Choi et al.,
“Field-based angle-resolved light-scattering study of single live cells,”
Opt. Lett., 3314 1596
–1598
(2008). https://doi.org/10.1364/OL.33.001596 OPLEDP 0146-9592 Google Scholar
|