Estimating mean particle volume and number from random sections by sampling profile boundaries

We describe some new shape-independent stereological estimates of particle mean volume and surface area. Finding volumes or surface areas of cell nuclei, from electron micrographs of random thin sections, is a central problem of biological stereology. The well-known point-sampled intercept (PSI) method samples profile interiors to find the volume-weighted mean volume. This can be used in place of the true mean volume, but to do so introduces bias when volumes vary a great deal, as they do in fixed specimens. Jensen and Gundersen quite recently extended the PSI estimator to provide particle surface area, with no bias in the case of uniform surface areas. Here we extend the PSI volume estimator in a different way, sampling profile boundaries rather than their interiors. We obtain a boundary-sampled intercept (BSI) volume estimator, simpler than the PSI surface area estimator, but also unbiased for uniform surface areas. Both of these estimators are attractive, for example, in measuring and counting cell nuclei, where membrane surface area varies less than volume. Furthermore, they have no shape bias whatsoever. This paper also examines the general relationship between boundary- and area-sampled estimates, and we clarify the formal connection between our volume estimator and the PSI surface area estimator. We also calculate and compare their theoretical efficiencies.     

Stereological ratio estimation based on counts from integral test systems

The stereological practice of using integral test systems in the estimation of the fundamental stereological ratios is studied in the light of recent theoretical developments in sampling. The estimation of a ratio is based on counts only, obtained from two, in general different, test sets constituting the integral test system. The ordinary ratio-of-sums estimator based on counts from uniformly positioned integral test systems is compared with two estimators based on non-uniform, weighted sampling. It is shown that the estimators based on weighted sampling are not, in general, unbiased. Furthermore, it is pointed out that the mean-of-ratios estimator based on replicated weighted sampling needs not have smaller MSE than the ordinary ratio-of-sums estimator based on replicated uniform sampling. The fact that the estimation is based on counts as opposed to complete 2-d observations does not necessarily mean a reduction in information. For certain types of stereological ratios, the ordinary ratio-of-sums estimator based on complete observation is shown, counter to intuition, to be less accurate than that based on simple and fast counting.     

Tissue shrinkage and unbiased stereological estimation of particle number and size*

This paper is a review of the stereological problems related to the unbiased estimation of particle number and size when tissue deformation is present. The deformation may occur during the histological processing of the tissue. It is especially noted that the widely used optical disector may be biased by dimensional changes in the z‐axis, i.e. the direction perpendicular to the section plane. This is often the case when frozen sections or vibratome sections are used for the stereological measurements. The present paper introduces new estimators to be used in optical fractionator and optical disector designs; the first is, as usual, the simplest and most robust. Finally, it is stated that when tissue deformation only occurs in the z‐direction, unbiased estimation of particle size with several estimators is possible.     

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.     

Fundamental stereological formulae based on isotropically orientated probes through fixed points with applications to particle analysis

A new set of fundamental stereological formulae based on isotropically orientated probes through fixed points is derived. Volume, surface area and integrals of curvature can be estimated by this method. Some of the estimators require, however, local 3-D information. The estimators can be applied locally to each of a sample of particles whereby moments of particle volume, particle surface area and particle integral of curvature can be estimated. The results are derived from a generalized version of the so-called Blaschke-Petkantschin formula.     

Estimating the surface area of nonconvex particles from central planar sections

In this paper, we present a new surface area estimator in local stereology. This new estimator is called the ‘Morse‐type surface area estimator’ and is obtained using a two‐stage sampling procedure. First a plane section through a fixed reference point of a three‐dimensional structure is taken. In this section plane, a modification of the area tangent count method is used. The Morse‐type estimator generalizes Cruz‐Orive's pivotal estimator for convex objects to nonconvex objects. The advantages of the Morse‐type estimator over existing local surface area estimators are illustrated in a simulation study. The Morse‐type estimator is well suited for computer‐assisted confocal microscopy and we demonstrate its practicability in a biological application: the surface area estimation of the nuclei of giant‐cell glioblastoma from microscopy images. We also present an interactive software that allows the user to efficiently obtain the estimator.     

Vertical LM sectioning and parallel CT scanning designs for stereology: application to human lung

Practical, unbiased stereological methods are described to estimate lung volume and external surface area, and total volume and surface area of relatively large and anisotropic structures (bronchi and arteries) inside the lung. The volume of each of five lung strata was estimated first by fluid displacement and then by computed tomography (CT) using Cavalieri's method; the reliability of CT was assessed through a calibration procedure, and image thresholding criteria for an accurate volume estimation using CT were established. The parallel, perfectly registered CT section images were also used to estimate the external surface area of each stratum by the spatial grid method. Unbiased estimation of internal surface areas in lung is a long-standing problem: since the structures are large and essentially void, large sections are needed; to facilitate identification, thin sections have to be used for light microscopy, and since such structures are anisotropic, the sections should be vertical. A practical stereological design is demonstrated here on an infant lung, which fulfils all these requirements. This study illustrates the potential of using unbiased stereology to characterize infant pulmonary hypoplasia.     

The spatial rotator

This paper presents a new local volume estimator, the spatial rotator, which is based on measurements on a virtual 3D probe, using computer assisted microscopy. The basic design of the probe builds upon the rotator principle which requires only a few manual intersection markings, thus making the spatial rotator fast to use. Since a 3D probe is involved, it is expected that the spatial rotator will be more efficient than the the nucleator and the planar rotator, which are based on measurements in a single plane. An extensive simulation study shows that the spatial rotator may be more efficient than the traditional local volume estimators. Furthermore, the spatial rotator can be seen as a further development of the Cavalieri estimator, which does not require randomization of sectioning or viewing direction. The tissue may thus be sectioned in any arbitrary direction, making it easy to identify the specific tissue region under study. In order to use the spatial rotator in practice, however, it is necessary to be able to identify intersection points between cell boundaries and test rays in a series of parallel focal planes, also at the peripheral parts of the cell boundaries. In cases where over‐ and underprojection phenomena are not negligible, they should therefore be corrected for if the spatial rotator is to be applied. If such a correction is not possible, it is needed to avoid these phenomena by using microscopy with increased resolution in the focal plane.     

On semiautomatic estimation of surface area

In this paper, we propose a semiautomatic procedure for estimation of particle surface area. It uses automatic segmentation of the boundaries of the particle sections and applies different estimators depending on whether the segmentation was judged by a supervising expert to be satisfactory. If the segmentation is correct the estimate is computed automatically, otherwise the expert performs the necessary measurements manually. In case of convex particles we suggest to base the semiautomatic estimation on the so‐called flower estimator, a new local stereological estimator of particle surface area. For convex particles, the estimator is equal to four times the area of the support set (flower set) of the particle transect. We study the statistical properties of the flower estimator and compare its performance to that of two discretizations of the flower estimator, namely the pivotal estimator and the surfactor. For ellipsoidal particles, it is shown that the flower estimator is equal to the pivotal estimator based on support function measurements along four perpendicular rays. This result makes the pivotal estimator a powerful approximation to the flower estimator. In a simulation study of prolate and oblate ellipsoidal particles, the surfactor also performs well for particles which are not extremely elongated. In particular, the surfactor is not very much affected by the singularity in the surfactor formula or by possible inaccuracies in the necessary angle measurements. We also assess the performance of the semiautomatic procedure in a study of somatostatin positive inhibitory interneurons from mice hippocampi. Only 35% of the cells needed to be analysed manually and an important decrease in workload was obtained by using the semiautomatic approach.     

Estimation of subcellular organelle volume from ultrathin sections through centrioles with a discretized version of the vertical rotator

A method for the fast and efficient estimation of the volume (but not surface area) of subcellular organelles is presented. It consists of a rotator/coaxial-section approach based on the Pappus theorem and represents a discretized version of the vertical rotator where, instead of measuring intercept lengths, the points in distance classes are counted. Centrioles serve as a unique reference 'double-point' with constant size allowing unbiased cell selection from the whole population with equal probability and without the disector application. The sandwich-like method of sample preparation allows comparison of control and experimental cases with the same errors caused by overlapping and overprojection. Test experiments demonstrated that the vertical discretized rotator was an efficient and precise tool for the estimation of volume and that a few independent sections of unknown thickness were sufficient for the quantification of one experimental point.     

Confocal microscopy and stereology: estimating volume, number, surface area and length by virtual test probes applied to three-dimensional images

The opportunities of confocal microscopy applied to morphometry of microscopical structures are presented and demonstrated on stereological methods based on evaluation of optical sections within a thick slice and using computer-generated virtual test probes. Such methods, allowing arbitrary orientation of the thick slice, can be used for estimating volume, number, surface area, and length. The methods using spatial grid of points, disector, fakir, and slicer probes are described and illustrated by different examples using our freeware 3DTOOLS software and their variance and applicability are discussed. It is shown that shifted triple or quadruple spatial grids of lines are very efficient for the surface area and volume estimation by the fakir method.     

Stereological estimation of the volume-weighted mean volume of arbitrary particles observed on random sections*

A stereological estimator of the weighted mean volume of particles of arbitrary shape is described. This unbiased estimator is based on simple point-sampling of linear intercept lengths. The complete absence of shape assumptions effectively breaks the long-standing ‘convexity-barrier’: the only requirement here is that individual particles can be unambiguously identified by their profiles on random sections. Practical details of the simple estimation procedure and an example with very irregular particles are reported. Finally, an estimator of the variance of the weighted distribution of particle volume is discussed. This estimator is also valid for particles of arbitrary shape. For any mixture of ellipsoids (spheres, oblates, prolates and triaxial ellipsoids) the estimator is reduced to a simple function of measurements of diameters in the section plane.     

Estimation of variance components of local stereological volume estimators: a pilot study

Local stereological techniques can be used for particle volume estimation based on information collected on a section plane through a reference point of the particle. We present methods for variability estimation of the local stereological volume estimators. This variability arises during the stereological estimation procedure and in the particle population. Both of these components can be estimated separately from planar sections. Our aim is to give a preliminary analysis of the possibility to include the particle structure interaction into the estimation procedure. For this reason, not only the section profiles, but also their locations, have to be recorded. The methods are applied for the sectional data obtained from neurons in the hippocampal brain region subiculum of four 3-month-old male Wistar rats. The proposed procedure enables one to obtain information about particle volume distribution.     

Unbiased volume estimation with coaxial sections: an application to the human bladder

A problem of current medical interest is how to estimate the volume of an object from coaxial sections. Such images can be produced in ultrasound scanning, when the planar scan sections all emanate from a common axis. Two unbiased volume estimators for this sampling regime, previously published without derivation, are derived and presented here in some detail. The estimators are based on an ancient theorem of Pappus of Alexandria (c. ad 320). One of the estimators was used to estimate the volume of urine in the bladder of six human volunteers. The case example was attractive because the corresponding volumes could be voided and directly measured after the scanning session, thus confirming the unbiasedness of the estimator used. Sample size guidelines are also given.     

Direct and efficient stereological estimation of total cell quantities using electron microscopy

Stereological estimation of total subcellular quantities in bioscience is presented in this report. Special emphasis is placed on the use of electron microscopy, which under certain circumstances may be combined with light microscopy. Three strategies based on the Cavalieri principle, the disector and local stereological probes through arbitrarily fixed points for estimation of total quantities inside cells are presented. The quantities comprise (total) number, length, surface area, volume or 3D spatial distribution for organelles as well as total amount of gold particles, various compounds or certain cytochemical markers.     

Stereological estimation of average cell volume in monolayer culture by combined light and electron microscopy

In order to convert stereological ratio estimates into ‘absolute’ values (‘per cell’ values) the average cell volume must be estimated. The present paper describes a stereological method based on well-known point counting procedures for the estimation of average cell volumes in monolayer cultures fixed in situ. This method involves estimation of the average attachment area per cell by light microscopy combined with estimation of the attachment membrane surface density by electron microscopy. There is no need for any assumption as to cellular or nuclear shape. The method has been tested on an established cell line, NHIK 3025, and shows a good accuracy. It has also been used to analyse the volume changes that take place in human monocytes during monolayer culture, demonstrating a 28-fold increase of the average cell volume over 10 days.     

Stereology of isolated objects with the invariator

The invariator is a new stereological design to generate motion invariant test lines in three dimensions on an isotropic plane through a fixed point. The theory has been published recently. The purpose of this paper is to illustrate the application of the invariator on a group of rat brains to estimate brain volume and external surface area. Each brain was first split into its two hemispheres and then embedded into a ball filled with agar following a configuration named the antithetic isector, with the idea of reducing the error variance. After rolling the ball at random it was scanned by magnetic resonance imaging into a stack of parallel systematic sections: this is the isotropic Cavalieri design which, combined with the antithetic isector idea, proves to be very accurate. The invariator used only an equatorial section of the ball, and in the present case the coefficient of error of the volume and surface area estimators of an individual brain was about 30%. As it is design unbiased, the invariator may prove its strength mainly to estimate population means.