Assessment of particle retention and clearance in the intrapulmonary conducting airways of hamster lungs with the fractionator1

A modified version of the fractionator was used to estimate the total number of polystyrene microspheres retained in the airways of hamster lungs at two different time points after inhalation. A systematic three-stage subsampling procedure with known sampling fractions was adopted. First, each lung was cut into slices, from which primary disectors were sampled systematically with a known sampling fraction. From each primary disector, smaller sub-disectors were subsampled, and the corresponding sampling fraction was estimated by point counting. Finally, a few particles were counted at the microscopic level in the sub-disectors, and the final estimate of total particle number (which is unbiased irrespective of any tissue deformations) was easily computed as a product of the counted number times the reciprocal of the successive sampling fractions. The error variance of each estimate was assessed from the data using a new estimator. An average of 6% of the deposited particles were retained on the epithelial surface of the intrapulmonary conducting airways shortly after the inhalation, from which at least one-third was already phagocytosed by macrophages. After 24 h, an average of 87% of the particles retained shortly after the inhalation had been cleared. The proportion of particles ingested by macrophages had increased to at least 87%, in three out of four animals studied.     

On the empirical variance of a fractionator estimate

The fractionator consists of several sampling stages with systematic sampling at each stage; data are collected only at the last stage. Therefore, predicting the error variance of a fractionator estimator is a non-trivial problem, because the observations are correlated in a complicated, unknown way. Gundersen proposed to split the material sampled at the first stage into two subsets, and to compute the variance of the pooled estimate empirically using the corresponding pair of observations made in these two subsets. The idea is very effective, but the estimator thus proposed needed some corrections. The purpose of this paper is to present an improved estimator of the coefficient of error of a fractionator estimator using Gundersen's design.     

Precision of the fractionator from Cavalieri designs

A popular procedure to predict the variance of the fractionator consists in splitting the initial collection of fragments into two subsets, in order to use the corresponding particle counts (or any other pertinent measure), in the calculation. The current formula does not account for local or ‘nugget’ errors inherent in the estimation of fragment contents, however. Moreover, it does not account for the fact that the contribution of the variability between fragments or slices should rapidly decrease as the sampling fraction increases. For these reasons, an update to the formula is overdue. It should be stressed, however, that the formula applies to Cavalieri slices designs – its application for arbitrary partition designs is therefore not warranted.     

Estimating the number of complex particles using the ConnEulor principle

An unbiased counting rule for the number of topologically simple objects of any shape, size and distribution in 3D space is a pertinent problem in stereology. Combining the disector principle with the object's 3D Euler number makes possible number estimation, which until now has been obtainable only by exhaustive serial sections. The disector is a set of two sections where the object's profiles in one section are compared with its profiles on the neighbouring section, and the number of new 2D topological events is recorded. In a disector of known volume the sum of topological events is a direct estimate of the disector contribution to the total Euler number, which forms the basis for an ultimate number estimator in 3D, the ConnEulor. The method is illustrated by an electron microscopic study of the number of mitochondria in the exocrine cells of the pancreas.     

Evaluation of variance models for fractionator sampling of trees

We compared the performance of several models for predicting, from small samples, the precision of estimates of the total number of blossoms on fruit trees obtained using a three‐stage fractionator, in which the sampling units were defined by the tree structure: (1) primary branches and stem (2) secondary branches and shoots and (3) flowering buds. The models considered were the semiempirical models of Cruz‐Orive (1990, 2004 ) (CO), a random sample model (SR), a Poisson model (P), successive differences (D) and repeated systematic sampling (R). Procedures that relied upon a single sample and a model of the variance (SR, P, D) were not able to predict the estimator variance because the considered structures strongly violated model assumptions. The semiempirical CO model performed acceptably in some cases where model assumptions were violated to a moderate degree. The repeated systematic sampling procedure, which does not rely upon a model of the variance, usually provided very good predictions when the resampling terms were distributed appropriately across more than one sampling stage.     

Using biased image analysis for improving unbiased stereological number estimation – a pilot simulation study of the smooth fractionator

The smooth fractionator was introduced in 2002. The combination of a smoothing protocol with a computer‐aided stereology tool provides better precision and a lighter workload. This study uses simulation to compare fractionator sampling based on the smooth design, the commonly used systematic uniformly random sampling design and the ordinary simple random sampling design. The smooth protocol is performed using biased information from crude (but fully automatic) image analysis of the fields of view. The different design paradigms are compared using simulation in three different cell distributions with reference to sample size, noise and counting frame position. Regardless of clustering, sample size or noise, the fractionator based on a smooth design is more efficient than the fractionator based on a systematic uniform random design, which is more efficient than a fractionator based on simple random design. The fractionator based on a smooth design is up to four times more efficient than a simple random design.     

Unbiased and efficient estimation of the total number of terminal bronchiolar duct endings in lung: a modified physical disector

A novel modification of the physical disector is described which was used to estimate the total number of terminal bronchiolar duct endings (TBDEs) in human infant lung. TBDEs are closed three-dimensional space curves of complex shape that are inherently difficult to count from histological sections. However, careful consideration of the microanatomy of the terminal duct endings provides us with the opportunity to define a very simple and unbiased counting rule. To apply the rule in practice we also need to determine a suitable disector height. Owing to the complex shape of the TBDE we had no prior knowledge of what disector height would be suitable for counting the TBDE structures. Exhaustive serial sectioning of complete TBDE structures was carried out and showed that any disector height under 90 μm would give unbiased counts. A further empirical study was then undertaken to determine the most efficient disector height. This was found to be 50 μm.
The total number of TBDEs in the upper lobe of the right lung of six human infants aged between 13 and 25 weeks was also estimated. The estimates of numerical density obtained with our modification of the physical disector were multiplied by estimates of lung lobe volume obtained using Cavalieri's Principle. The total number of TBDEs in the lobes ranged from 15 323 to 57 768, with a mean of 40 306. The average coefficient of error of the number estimates was 19%, which was deemed precise enough given the biological coefficient of variation between TBDE number of 36%.     

Estimating the length of a bounded curve in three dimensions using total vertical projections

A new method is proposed to estimate the total length of a bounded, isolated linear feature in three dimensions from ‘total vertical projections’, obtained by rotating the curve about a fixed axis (arbitrarily called ‘vertical’) and projecting it onto a fixed vertical plane. No sections are required. Properly stained and embedded neuron dendrites, mycelial trees, fluorescent cytoskeletal filaments within a cell, etc., are candidate specimens for the method, especially in combination with the new devices for non-invasive three-dimensional microscopy. It is necessary that the specimen curve is rigid (i.e. of constant shape), that its length density is not too high (so that overlapping effects are not important) and that the embedding medium is fairly transparent. Given these requirements, the method can be very accurate and convenient to use, as is exemplified here.     

Determination of the coordination number in disordered packings of equal spheres

A computation of the coordination number in a packing of equal spheres can be performed from the information obtained along a random line through the packing. It involves the distribution of the number per unit length of separators between intercepts of spheres with the line. The number of contacts per unit volume is derived from the slope of the distribution at the origin. The method is validated by computer experiments on four kinds of simulated ordered packings which are intercepted by a set of random lines (instead of a single one) in order to compensate biased results due to the small size of these packings: the comparison is then extended to a series of disordered dense packings of equal spheres. We present some possible applications of this study to stereological and physical problems.     

The efficiency of systematic sampling in stereology and its prediction*

The superior efficiency of systematic sampling at all levels in stereological studies is emphasized and various commonly used ways of implementing it are briefly described. Summarizing recent theoretical and experimental studies a set of very simple estimators of efficiency are presented and illustrated with a variety of biological examples. In particular, a nomogram for predicting the necessary number of points when performing point counting is provided. The very efficient and simple unbiased estimator of the volume of an arbitrary object based on Cavalieri's principle is dealt with in some detail. The efficiency of the systematic fractionating of an object is also illustrated.     

The total number of neurons in the human neocortex unbiasedly estimated using optical disectors

An efficient method is presented for obtaining, in under 4h, an unbiased estimate of the total number of neurons in the human neocortex, with a coefficient of error on the estimate of ～ 5%. The novel sampling scheme used in this study is unbiased and was designed so that only a small amount of neocortical grey matter had to be removed. Hence, the majority of the cerebral grey matter and all the internal grey matter was left intact for further resampling and analysis. Each cerebral hemisphere was subdivided into the four major neocortical regions, sliced coronally at 7-mm intervals and the volume of the neocortex determined using Cavalieri's principle. Uniform sampling of neocortex was performed in the hemisphere followed by regional subsampling with a varying sampling fraction being taken from each region. Neuronal density estimates were made in thick plastic sections using optical disectors. Shrinkage estimates were made in parallel with the number estimates and found to be negligible. The total number of neocortical neurons in the right hemisphere of five normal 80-year-old men was found to be 13·7 × 10⁹ with an inter-individual coefficient of variation of 12%.     

Effect of lubricating oil additives on particle size distribution and total number concentration in diesel engine

This paper investigates the characteristics of particle size distribution in exhaust gas of engine fuelled with pure diesel and with diesel mixed with base oil or with oil additives. The experiments are conducted on a turbocharged diesel engine with fast particulate spectrometer DMS 500 connected to the exhaust pipe. Base oil and two kinds of commonly used lubricating oil additives, antioxidant additives and antifoaming additives, are chosen to be added into the fuel, with the concentrations being 0.5%, 1.0% and 1.5% of fuel weight individually. The particle size distribution is measured under medium load (100 Nm) and full load at different speeds. The results indicate that the existence of base oil or oil additives shows great influence on particle size distribution. Copyright © 2012 John Wiley & Sons, Ltd.     

Estimation of the number of stellate cells in a liver with the smooth fractionator

To better evaluate the activation and proliferative response of hepatic stellate cells (HSC) in hepatic fibrosis, it is essential to have sound quantitative data in non‐pathological conditions. Our aim was to obtain the first precise and unbiased estimate of the total number of HSC in the adult rat, by combining the optical fractionator, in a smooth sampling design, with immunocytochemistry against glial fibrillary acidic protein. Moreover, we wanted to verify whether there was sufficiently relevant specimen inhomogeneity that could jeopardize the high expected estimate precision when using the smooth fractionator design for HSC. Finally, we wanted to address the question of what sampling scheme would be advisable a priori for future studies. Microscopical observations and quantitative data provided no evidence for inhomogeneity of tissue distribution of HSC. Under this scenario, we implemented a baseline sampling strategy estimating the number (N?) of HSC as 207E⁰⁶ (CV = 0.17). The coefficient of error [CE(N?)] was 0.04, as calculated by two formerly proposed approaches. The biological difference among animals contributed ? 95% to the observed variability, whereas methodological variance comprised the remaining 5%. We then carried out a half reduction of sampling effort, at the level of both sections and fields. In either occasion, the CE(N?) values were low (? 0.05) and the biological variance continued to be far more important than methodological variance. We concluded that our baseline sampling (counting 650–1000 cells/rat) would be appropriate to assess the lobular distribution and the N? of HSC. However, if the latter is the only parameter to be estimated, around half of our baseline sampling (counting 250–600 cells/rat) would still generate precise estimates [CE(N?) < 0.1], being in this case more efficient to reduce the number of sections than to reduce the sampled fields.     

Precision of Cavalieri sections and slices with local errors

Cavalieri sections — and more recently Cavalieri slices, especially in combination with non-invasive scanning — are widely used to estimate volumes. Physical Cavalieri slices are also increasingly used to estimate neuron numbers via the optical fractionator. In either case, the prediction of the error variance is important to assess optimal sample sizes. The error variance consists of two components, one due to the variation among the true contents of sections or slices, and the other due to local or ‘nugget’ errors. The latter may arise for instance when estimating section areas by point counting, or when counting discrete particles in slices or disectors. In this paper, a fairly comprehensive set of prediction formulae is presented to separate both variance components.     

The new stereological tools in metallography: estimation of pore size and number in aluminium

A number of simple, unbiased and efficient methods are now available in stereology for estimating the number and size of arbitrary particles or voids in a material, with the only assumption that particles must be identifiable on serial sections or confocal planes through the material. In recent years, these methods have been developed and applied mainly in a biomedical context: this paper reviews and illustrates them with the aid of a metallographic example, namely the pore population of a sand-cast aluminium alloy. Our goal is to convey the fact that stereology is sampling in three dimensions, and therefore its principles remain valid and applicable in no matter what context. The disector, the selector, and an indirect method to estimate the distance between two parallel planes of polish are thereby illustrated. It is also shown how to split the error variance of the estimator of the pore volume fraction (‘porosity’) into the three components due to blocks, sections within blocks and systematic point counting on sections.     

The efficiency of systematic sampling in stereology — reconsidered

In the present paper, we summarize and further develop recent research in the estimation of the variance of stereological estimators based on systematic sampling. In particular, it is emphasized that the relevant estimation procedure depends on the sampling density. The validity of the variance estimation is examined in a collection of data sets, obtained by systematic sampling. Practical recommendations are also provided in a separate section.     

Global spatial sampling with isotropic virtual planes: estimators of length density and total length in thick, arbitrarily orientated sections

Existing design-based direct length estimators require random rotation around at least one axis of the tissue specimen prior to sectioning to ensure isotropy of test probes. In some tissue it is, however, difficult or even impossible to define the region of interest, unless the tissue is sectioned in a specific, nonrandom orientation. Spatial uniform sampling with isotropic virtual planes circumvents the use of physically isotropic or vertical sections. The structure that is contained in a thick physical section is investigated with software-randomized isotropic virtual planes in volume probes in systematically sampled microscope fields using computer-assisted stereological analysis. A fixed volume of 3D space in each uniformly sampled field is probed with systematic random, isotropic virtual planes by a line that moves across the computer screen showing live video images of the microscope field when the test volume is scanned with a focal plane. The intersections between the linear structure and the virtual probes are counted with columns of two dimensional disectors.
Global spatial sampling with sets of isotropic uniform random virtual planes provides a basis for length density estimates from a set of parallel physical sections of any orientation preferred by the investigator, i.e. the simplest sampling scheme in stereology. Additional virtues include optimal conditions for reducing the estimator variance, the possibility to estimate total length directly using a fractionator design and the potential to estimate efficiently the distribution of directions from a set of parallel physical sections with arbitrary orientation.
Other implementations of the basic idea, systematic uniform sampling using probes that have total 3D × 4π freedom inside the section, and therefore independent of the position and the orientation of the physical section, are briefly discussed.     

A heuristic method to determine the number of pallets in a flexible manufacturing system with several pallet types

This article presents a new heuristic method, called PANORAMA, that enables the determination of an appropriate number of pallets of each type that should circulate in a flexible manufacturing system. Production and return-on-investment constraints are considered. Statistical results derived using data from three existing flexible manufacturing systems, as well as randomly generated data, demonstrate the accuracy of the heuristic proposed.