Abstract
Recently, the distribution of radioactivity among a population of cells labeled with ^{210}Po was shown to be well described by a lognormal (LN) distribution function (J Nucl Med. 2006;47:1049–1058) with the aid of autoradiography. To ascertain the influence of Poisson statistics on the interpretation of the autoradiographic data, the present work reports on a detailed statistical analysis of these earlier data. Methods: The measured distributions of αparticle tracks per cell were subjected to statistical tests with Poisson, LN, and Poissonlognormal (PLN) models. Results: The LN distribution function best describes the distribution of radioactivity among cell populations exposed to 0.52 and 3.8 kBq/mL of ^{210}Pocitrate. When cells were exposed to 67 kBq/mL, the PLN distribution function gave a better fit; however, the underlying activity distribution remained lognormal. Conclusion: The present analysis generally provides further support for the use of LN distributions to describe the cellular uptake of radioactivity. Care should be exercised when analyzing autoradiographic data on activity distributions to ensure that Poisson processes do not distort the underlying LN distribution.
Autoradiography has been used for decades to quantify the distribution of radioactivity in tissues after administration of radiochemicals (1–6). Recently, αparticle track autoradiography was used to determine the distribution of ^{210}Pocitrate among a clonal cell population of Chinese hamster V79 lung fibroblasts (7). Analysis of the data revealed a lognormal (LN) distribution of cellular activity, a distribution that can have a substantial effect on the biologic response of the cell population (7,8). To arrive at this distribution, the number of αparticle tracks per cell was scored in more than 1,000 cells. However, it was necessary to limit the number of counted tracks per cell to fewer than 10 to ensure accurate scoring. The statistical uncertainties associated with measuring these low numbers of tracks per cell can affect the measuredtrack distribution and thereby potentially distort the underlying activity distribution among the cell population (9,10).
As pointed out by Kvinnsland et al. (9), it is the statistical nature of radioactive decay that can influence the measured distribution of αparticle tracks, particularly at the low numbers of tracks per cell studied by Neti et al. (7). Therefore, it may be necessary to “peel off” the distribution associated with the statistics of radioactive decay to arrive at the underlying activity distribution. It has been hypothesized that in radioactive decay, all atoms are identical and independent, and the chance for an atom to disintegrate during a given time interval is the same for all time intervals of equal size. Accordingly, after the discovery of radioactivity, initial experiments suggested that the decay process, and particle counting associated with its observation, is well approximated by simple Poisson processes. Subsequently, during a study of αparticle emission rates, a quantum selfinterference phenomenon was proposed to account for 1/f fluctuations (11). Experiments showed that αemission is in fact not a simple Poisson process, and particle counting is not adequately described by Poisson statistics (12). However, another report states that 1/f fluctuations are not present in the decay of ^{210}Po at decay frequencies in excess of 10^{−6} Hz (13). The mean track numbers for the range of ^{210}Pocitrate concentrations discussed in the present report fall in the window of approximately 0.16–1.6 × 10^{−6} Hz. Consequently, it is unlikely that 1/f fluctuations play a significant role in interpreting our experimental αparticle track data (7). Nevertheless, Poisson statistics may play a role; it is necessary, therefore, to examine whether the LN distribution obtained from our ^{210}Po αparticle track data is influenced by the Poisson distribution associated with its decay. Accordingly, in this report we subject our αparticle track data to statistical tests with Poisson, LN, and Poissonlognormal (PLN) distribution functions (14).
MATERIALS AND METHODS
Summary of Published Experimental Data and Analyses
The methods used to obtain the αparticle track data statistically analyzed in the current study are explained in detail in an earlier report (7). Briefly, cultured cells were labeled with ^{210}Pocitrate, washed, and coated with autoradiographic emulsion. Decays were allowed to accumulate for various times to obtain scoreable track data that covered the entire range of cellular activities encountered in a given labeled cell population. Accordingly, the emulsion was developed for 2 or 3 different exposure times. For each emulsion exposure time, 500–1,200 cells were examined, and the number of tracks in each cell was recorded. Cells with 0–9 tracks per cell were scored as such. Cells with more than 9 tracks could not be accurately counted and were simply scored as greater than 9 tracks. The entire process was repeated for ^{210}Pocitrate–labeling concentrations of 0.52, 3.8, and 67 kBq/mL, which resulted in mean cellular activities of 0.054, 0.12, and 1.8 mBq per cell. The numbers of tracks per cell from the short exposure times were decaycorrected to the longest exposure time and then normalized. All datasets for a given labeling concentration were then combined into a single convolved track distribution with the longest exposure time data providing the distribution for the lower multiplicity of convolved tracks and the short exposure time data providing the distribution of the higher multiplicity. A nonlinear leastsquares (LS) fit of the resulting convolved track distribution was performed with an LN function (7).
Experimental Mean Number of αParticle Tracks per Cell
The arithmetic mean number of tracks per cell <n> can be calculated for the track data corresponding to the shortest emulsion exposure times for each labeling concentration. These data have no cells with more than 9 tracks, and therefore the mean is obtained by simply tallying the total number of tracks and dividing by the number of cells. A different approach was required to calculate <n> for the datasets corresponding to longer emulsion exposure times that contained cells with more than 9 tracks. Regardless of the mathematic nature of the track distribution, the total number of tracks n_{total} recorded in the cell population at time point t_{1} can be projected to a later time point t_{2} according to Equation 1,Eq. 1where T_{p} is the physical halflife of the radionuclide (T_{p} = 138 d for ^{210}Po). In the same manner, <n(t_{1})> (shortest exposure time with all cells having <9 tracks) can be used to calculate <n(t_{2})> (longer exposure time with some cells having >9 tracks).
Statistical Distributions and Analyses
Poisson Distribution.
If each cell in the population had the same activity, one would anticipate a Poisson distribution of tracks per cell. A 1/f distribution may be expected according to one study (12); however, another report indicates that 1/f fluctuations are not present in the decay of ^{210}Po for frequencies in excess of 10^{−6} Hz (13), a frequency range that encompasses our data. Thus, in this report, the Poisson distribution is assumed to adequately describe the statistical distribution associated with the decay of ^{210}Po. After finding <n> for each set of track distribution data, we calculated the Poisson probability of n discrete tracks per cell, given byEq. 2For noninteger values of n, which may result because of decay correction (Eq. 1), truncated values are used to calculate the probability.
LN Distribution.
The LN probability density function is given byEq. 3where μ_{n} is the scale parameter, σ is the shape parameter, and g is a constant. The scale parameter is related to the mean according to μ_{n} = ln <n> − σ^{2}/2. The properties of this function and its use in analysis of the distribution of tracks per cell are given in a previous study (7). Detailed information on LN distributions can be found in other studies (15,16) as well.
PLN Distribution.
The PLN compound probability of obtaining a realization n and all its possible Poisson realizations k is given in an earlier study (17),Eq. 4
χ̂^{2} Analyses.
To quantitatively assess which of the above statistical distributions best describes the observed track distribution data, a reduced χ^{2} analysis () was performed. This quantity is given by Equation 5,Eq. 5where n_{i} is the number of tracks per cell for the ith experimental data, n_{expected,i} is the expected value based on the assumed distribution, ν is the degree of freedom (ν = i_{max} − p − 1), and p is the number of parameters in the assumed distribution.
LS Fit.
LS fits of the data to an LN distribution function (Eq. 3) were performed with software (Sigmaplot; Systat Software). The LS fit of the LN distribution function to the data corresponding to emulsion exposure times where more than 9 tracks were observed was performed by 2 methods. In the first method, <n> was calculated by decaycorrecting <n(t_{1})> to obtain <n(t_{2})> and then constraining <n(t_{2})> to obtain σ. In the second method, <n> and σ were allowed to vary. No constraints were placed on the LS fit of the LN distribution function to the convolved data.
RESULTS
Statistical Analysis for Cells Exposed to 0.52 kBq/mL
To study the effect of Poisson statistics on our analyses, the 3 discrete raw datasets from Figure 3A in an earlier study (7) (0.52 kBq/mL) are first revisited. The 3 sets of track distributions were acquired from autoradiographs that were developed at 7, 26, and 52 d, as shown in Figures 1A–1C of the current study. Each set of track distribution data includes the number of cells scored with 0–9 tracks per cell as well as the number of cells with an unscoreable number of tracks (>9 tracks). The t = 7 d data have a maximum of 6 tracks in any given cell, and no cells have more than 9 tracks. The mean value for the t = 7 d data is 1.0 tracks per cell. The values of <n> for t = 26 d and t = 52 d cannot be obtained directly from the respective datasets because some of the cells have greater than 9 tracks. Thus, the projected <n> values obtained with Eq. 1 are 3.6 for the 26d data and 6.7 for the 52d data. Using these <n> values, the Poisson probabilities for each of 0–9 discrete tracks per cell are calculated at t = 7, 26, and 52 d. These probabilities can be used to assess the extent to which our experimental data represent a Poisson distribution. As shown in Figure 1A, the t = 7 d experimental data follow the general trend of a Poisson distribution, which suggests that Poisson statistics may have some impact on these data. The t = 26 d (Fig. 1B) and t = 52 d (Fig. 1C) datasets do not follow a Poisson distribution, although this does not rule out a Poisson component in these distributions.
To assess the possible Poisson component, we return to the t = 52 d dataset. This is the primary dataset, covers the range of 0–9 tracks in our convolved data, and primarily dictates the shape of the probability distribution. For the t = 52 d data, the values for the Poisson, LN, and PLN distributions were 149, 3.2, and 6.5. The latter 2 values were obtained via a minimization procedure with the σ. The values of σ were 0.81 and 1.1 for the LN and PLN distributions, respectively. The lowest value was obtained for the LN distribution as shown in Table 1, suggesting that Poisson corrections are not required for the t = 52 d data. Furthermore, the σ obtained with the LN distribution corresponding to the discrete t = 52 d data was the same as that obtained by LS fit of the convolved data, which included the t = 7, 26, and 52 d data.
Now that the LN distribution has been firmly established with the primary t = 52 d data, we return to the ancillary t = 7 and 26 d data where the Poisson distribution may play a more significant role. These data represent cells from the tail of the LN distribution and therefore contain little information about the overall shape of the entire activity distribution in the cell population. Accordingly, although an LN distribution with σ = 0.8 fits these data extremely well, it is not appropriate to apply this distribution to this subset of data. For the same reason, the PLN distribution cannot be used to assess the significance of the Poisson distribution in these data. This makes it difficult to tease the possible Poisson influence out of the underlying distributions in the t = 7 and 26 d datasets. With this limitation in mind, it is worth revisiting the convolved data where these distributions are applicable. The 52d data were used for the low multiplicity of convolved tracks per cell (1–9 tracks), the 26d data were used for the middle multiplicity (5–15 tracks), and the 7d data were used for high multiplicity (10–32 tracks). Within the overlapping regions (Fig. 1D), considerable agreement exists between the data points from independent datasets. Furthermore, an analysis of the convolved data gives 7.7 × 10^{7}, 3.6, and 12 for the Poisson, LN, and PLN. The difference between the convolved data and the 3 predicted probability curves can be seen in the inset of Figure 1D at very small probabilities. Again, the data are best described by the LN distribution (Table 1). This suggests that the convolved t = 7 d and 26 d data follow the expected trend of the LN distribution, and any Poisson component in these data has a negligible impact on the overall distribution.
Statistical Analysis for Cells Exposed to 3.8 kBq/mL
An approach similar to the one described earlier in this study was also applied to the 2 discrete datasets that were obtained at a labeling concentration of 3.8 kBq/mL. The 2 sets of track distributions were acquired from autoradiographs that were developed at t = 0.67 and 4 d, as shown in Figures 2A and 2B. In this case, the t = 4 d dataset is the primary dataset, covers the range of 0–9 tracks of convolved data, and primarily dictates the shape of the probability distribution. The <n> for the t = 0.67 d data, where all cells have fewer than 9 tracks per cell, is 1.2 tracks per cell. The mean number of tracks per cell for the 4d data can be obtained with Equation 1, which results in <n> = 7.3. The Poisson probabilities for each of the 0–9 discrete tracks per cell are calculated for <n> = 1.2 (t = 0.67 d) and <n> = 7.3 (t = 4 d) (Figs. 2A and 2B). The values for the Poisson distribution are 96 for the 0.67d and 25 for the 4d datasets. The values for all 3 distribution functions are given in Table 1. On minimization of , the values of σ were 0.64 and 0.60 for the LN and PLN, respectively. The lowest value was obtained for the LN distribution (Table 1), thereby again suggesting that Poisson corrections are not required for the t = 4 d data. Furthermore, the σ obtained for the LN distribution corresponding to the t = 4 d data was the same as that obtained by an LS fit to the convolved data. Agreement exists between the data points within the overlapping regions of the convolved datasets (Fig. 2C). The arithmetic mean (<n>) value of 6.1 for the convolved dataset was obtained from an LS fit of the convolved data to the LN distribution. Furthermore, analysis of the convolved data gives 1.2 × 10^{7}, 2.8, and 2.1 for the Poisson, LN, and PLN. The difference between the convolved data and the 3 predicted probability curves is shown in the inset of Figure 2C. Again, the data are best described by the LS fit of the convolved data to an LN distribution function (Table 1, ).
Statistical Analysis for Cells Exposed to 67 kBq/mL
Finally, the results are somewhat different at the highest labeling concentration of 67 kBq/mL. In this case, the 2 sets of track distributions were acquired from autoradiographs that were developed at t = 0.25 and 1 d (Figs. 3A and 3B). The <n> for the t = 0.25 d data is 1.4 tracks per cell and, using Equation 1, <n> = 5.7 for the 1d data. The Poisson probabilities for each of 0–9 discrete tracks per cell are calculated with Equation 2 for each <n>. As shown in Figure 3A, the t = 0.25 d dataset does not follow a Poisson distribution (). However, the t = 1 d data are better described by the Poisson distribution () than by the LN distribution (). The PLN distribution function provides the best fit to the data, suggesting a significant Poisson contribution (). The values of σ were 0.54 and 0.28 for the LN and PLN, respectively. The σ obtained for the LN distribution corresponding to t = 1 d data was the same as that obtained by an LS fit for the convolved data (Table 1). Again, within the overlapping regions of the convolved data considerable agreement exists between the data points from independent datasets (Fig. 3C). Furthermore, analysis of the convolved data gives 9.6 × 10^{6}, 8.4, and 1.7 for the Poisson, LN, and PLN. A comparison between the convolved data and the 3 predicted probability curves is provided in the inset of Figure 3C. The data are much better described by the LS fit of the convolved data to an LN distribution function (Table 1, ) than to a Poisson distribution (), but the data are best explained by the PLN model ().
DISCUSSION
As described above, a variety of approaches were used to analyze our ^{210}Po αparticle track distribution data. These involve statistical analyses of both the raw track data and the convolved track data using Poisson, LN, and PLN functions. The results, summarized in Table 1, indicate that the Poisson distribution alone does not adequately describe the observed track distributions as it would if the cellular activity was distributed equally among all the cells in the population. In fact, with the exception of the data corresponding to the highest concentration of ^{210}Po, a pure LN distribution generally provides the best description of the data. In the case of the highest concentration (67 kBq/mL), the data are best described by the PLN distribution, which suggests that the Poisson statistics associated with radioactive decay play a significant role in the analysis of this set of autoradiographic data. Nevertheless, the underlying distribution of activity in this cell population is well described by an LN distribution.
It is interesting to more closely examine the conditions under which Poisson statistics play a significant role in interpreting our αparticle track data. Figures 1C, 2B, and 3B show that the autoradiographic track data appear to become more Poissonian as the concentration of ^{210}Pocitrate in the culture medium is increased. This is supported by the values for the Poisson, LN, and PLN distributions. For the 1d dataset (67 kBq/mL), the corresponding values are 1.1, 5.6, and 0.66 (Table 1). The reason for the increased Poisson influence at higher concentrations of ^{210}Pocitrate may be related to a decrease in the time interval for the observation (i.e., shortemulsion exposure time). However, the values for Poisson fits to the shortinterval 0.25d (67 kBq/mL) dataset and the 0.67d (3.8 kBq/mL) dataset are 91 and 96 (data are not shown in Table 1). These high values indicate a poor fit, which suggests that the time interval may not be a significant factor.
It is also possible that the Poisson distribution plays a bigger role with increasing concentration of ^{210}Pocitrate because of a reduction in nonuniformity in the distribution of radioactivity among the cell population. To further pursue this possibility, the value of the LN σ obtained by the various methods described above are graphically compared in Figure 4. The dashed line connects what we believe are the best values for σ, where best is defined as the lowest (Table 1). On the basis of the linearity of the 3 data points on a logarithmic scale, the σ is exponentially dependent on the extracellular concentration. A decreasing value of σ implies more uniform distribution of radioactivity among the cell population. This implies that the nonuniformity of the cellular uptake of ^{210}Pocitrate decreases with increasing extracellular concentration (7,8). Figure 6 in an earlier study (7) shows that this will result in a more exponential survival curve rather than the saturating doseresponse curves that are expected for LN distributions with large σ values.
In our earlier article, we used the convolution approach to analyze the experimental αparticle track data (7). Kvinnsland et al. expressed concern about the effect of the convolution approach on the impact of Poisson statistics on the resulting distribution (9). The overlapping regions of convolved data from independent datasets (Figs. 1D, 2C, and 3C) suggest that the data can be convolved. The σ obtained by LS fit of the primary dataset (1–9 tracks, t_{2} > t_{1}) and the convolved data to an LN function are also similar for the cases of 0.52 and 3.8 kBq/mL, which supports the convolution approach as well. However, as shown in Figure 4, the convolution approach is not satisfactory for the case of 67 kBq/mL datasets. In this instance, as described above, the underlying LN distribution must be extracted from the measured distribution using the PLN function. Thus, caution must be exercised when implementing the convolution approach. The σ values obtained directly from the 52d data and the 4d data are very similar to those values obtained from the convolved data. Therefore, provided that one has a single autoradiograph that covers most of the track distribution (e.g., Figs. 1C or 2B), one may be able to obtain σ without resorting to convolution.
To test the potential use of LN distributions to describe the cellular uptake of other radiochemicals within a population of cells, the data of Lehman and Scavia are revisited (18). Their cells were labeled with ^{33}P, an emitter of βparticles with a mean energy of 76.9 keV, under unstirred and stirred conditions, and the frequency distribution of autoradiographic grains per cell was measured. Lehman and Scavia noted that the distribution of background grains per cell was adequately described by a Poisson distribution, but the experimental treatments (stirred and unstirred mixtures) were not explained by a simple Poisson distribution (18). Consequently, they used a Neyman's Type A distribution, which consists of a composite of 2 Poisson processes (one process that describes graincounting statistics and another process that pertains to the randomness of radioactive decay (19)). They fit their data to this distribution and obtained values of 0.68 and 1.7 for the stirred and unstirred conditions, respectively. To assess the capacity of the LN distribution to describe their data, the arithmetic mean (<n>) was calculated for each mixture, and an analysis was conducted for Poisson, LN, and PLN. Furthermore, an unconstrained LS fit of the data to an LN function was performed. The results of these fits show that the LN, PLN, and LS fit to LN generally yield similar values of and LN σ (Table 2). Thus, the underlying distribution of radioactivity among the cell population is well described by an LN distribution.
The analyses above support the conclusion that the distribution of radioactivity in a cell population can often be well represented by an LN distribution. LN distributions have been shown to describe a variety of natural processes. Bulmer (20) fitted the species abundance data, which showed that truncated, grouped LN distributions provided a satisfactory fit over the logarithmic and PLN models, omitting zero class. Other models exist that can be useful for explaining data with a priori knowledge of the underlying distribution, such as multivariate PLN distributions (21), bivariate Poisson distributions (19), and so on. However, LN distributions are ubiquitous and have been observed across many fields (22,23).
Finally, this exercise was undertaken to demonstrate the statistical significance of the LN distribution that was obtained for ^{210}Pocitrate using an αparticle track autoradiographic approach (7). The reliability of this approach stems from its use of emitted radiations to ascertain the distribution. Other experimental techniques, such as flow cytometry (24,25), can also be used to derive the distribution. Because of the ubiquitous presence of LN distributions, many investigators studying radiobiologic responses to radiopharmaceuticals and other radiochemicals may find this distribution useful to fold into their doseresponse models (8,24,26). Its implementation is facilitated by several factors. First, it is an analytic function that is described by only 2 parameters (σ, μ). Second, the LN probability density function is provided in standard subroutine libraries (e.g., National Algorithm Group) for computational purposes.
CONCLUSION
Statistical tests show that the LN distribution function is favored as a general form to describe the distribution of cellular uptake of ^{210}Po among a population of cells exposed to the same concentration of ^{210}Pocitrate. Furthermore, this distribution is likely to be applicable to a variety of radiopharmaceuticals.
Acknowledgments
This work was supported in part by U.S. Public Health Service grant R01CA83838.
Footnotes

COPYRIGHT © 2008 by the Society of Nuclear Medicine, Inc.
References
 Received for publication November 5, 2007.
 Accepted for publication December 6, 2007.