The first mathematically unbiased estimates of neocortical cell numbers are presented from the developing pig brain, including a full description of tissue processing and optimal sampling for application of the stereological optical fractionator method in this species. The postnatal development of neocortical neurons and glial cells from the experimental Göttingen minipig was compared with the postnatal development of neocortical neurons in the domestic pig. A significant postnatal development was observed in the Göttingen minipig brain for both neuronal (28%; P=0.01) and glial cells (87%; P<0.01). A corresponding postnatal development of neurons was not detected in the domestic pig brain. The reason for this strain difference is not known. The mean total number of neocortical neurons is 324 million in the adult Göttingen minipig compared with 432 million in the domestic pig. The glial-to-neuron cell ratio is around 2.2 in the adult Göttingen minipig. Based on these results, the domestic pig seems to be a more suitable model for evaluating the effects of developmental insults on human brain growth and neuronal development than the Göttingen minipig.

The Göttingen minipig and the domestic pig are increasingly used as non-primate models in basic experimental studies of neurological diseases. The gyrated pig brain is more similar to the primate brain than lissencephalic brains from small laboratory animals. The pig is affordable, it is easily handled and its use may potentially avoid some of the ethical considerations concerning the use of primates as laboratory animals.

The pig has previously been considered a potential animal model for evaluating the effects of developmental insults on human brain growth and development (Dickerson and Dobbing,1966; Book and Bustad,1974; Dobbing and Sands,1979; Pond et al.,2000). In general, the mammalian brain appears to go through a momentary period of rapid growth, exemplified by a characteristic growth rate curve when the percent of adult brain mass is plotted against age(Dobbing and Sands, 1979). One of the most important interspecies differences seems to be the complexity of the final product as well as the timing of the brain growth spurt in relation to birth (Dobbing, 1974). The timing of the brain growth spurt can be used to roughly categorize mammalian species as prenatal, perinatal or postnatal developers(Dobbing and Sands, 1973). In a comparison of seven mammalian species, it was demonstrated that the pig brain, like that of humans, develops perinatally, with a brain growth spurt extending from midgestation to early postnatal life(Dickerson and Dobbing, 1966; Dobbing and Sands, 1979; Pond et al., 2000). This is in contrast to other mammalian species, e.g. the brain of guinea pig, sheep and monkey, which has a prenatal growth spurt, or the brain of rat and rabbit,which develops postnatally (Dobbing and Sands, 1979). The development of the pig brain is also considered more similar to the human brain with respect to myelination, compositions and electrical activity (Dickerson and Dobbing,1966; Fang et al.,2005; Flynn, 1984;Pampliglione, 1971; Thibault and Margulies, 1998).

The traditional view of the primate neocortex is that it is structurally stable and that neurogenesis and synapse formation occur during early development (Bourgeois et al.,1994; Rakic,1985b). Quantitative studies based on DNA quantification in the human brain have indicated that the major phase of neuronal multiplication occurs during the first half of gestation, prior to the numerically larger but slower phase of glial multiplication, which continues into the first postnatal years (Dobbing and Sands,1973; Dobbing,1974). A two-phased growth pattern has similarly been observed in a stereological study on total cell numbers in the developing human fetal forebrain (Samuelsen et al.,2003). A clear cell discrimination has, however, not been performed in early fetal life.

In an attempt to further evaluate the pig as a potential model for human brain development and to provide a quantitative structural basis for comparative and experimental studies, a number of quantitative examinations on the neonate and adult pig brain have been initiated. Quantitative data are obtained using assumption-free stereological methods. The methods are designed to describe quantitative parameters without assumptions about shape, size,orientation and distribution of cells in the reference space and are based on established procedures for systematic, uniformly random sampling that are superior in precision compared with results obtained by independent sampling(Gundersen and Jensen, 1987). In the present paper, neocortical cell numbers were obtained using the optical fractionator method (Gundersen,1986; West et al.,1991). A detailed description of tissue processing and optimal sampling procedures for application in the pig brain is presented.

Experimental animals

The Göttingen minipig

The Göttingen minipig was developed in 1961–1962 at the Institute of Animal Breeding and Genetics of the University of Göttingen(Germany). The present characteristics of the Göttingen minipig, as a small, white miniature pig with good fertility and stable genetics, were obtained as a result of cross-breeding the Minnesota minipig with the Vietnamese potbelly pig and the German Landrace(Bollen and Ellegaard, 1997; Damm Jorgensen, 1998). The Göttingen minipig in Denmark originates from caesarean sections performed on German sows in the early 1990s. The offspring was transferred to a barrier facility and kept in a non-contaminated environment. Only three populations of limited population size exist in the world today.

The Göttingen minipig is an outbred animal, with less than 5%in-breeding (Glodek, 1986). The newborn Göttingen minipig has a body mass of ∼350–450 g. Boars become sexually mature at an age of 3–4 months, weighing 6–8 kg, while sows are sexually mature at an age of 4–5 months, weighing 7–9 kg. The gestational period is 112–114 days and the average litter size is 5–6 animals. The adult mass of the Göttingen minipig, at an age of two years, is 35–40 kg(Bollen and Ellegaard, 1997; Damm Jorgensen, 1998).

The Göttingen minipig of today is not gnotobiotic but is kept in barriers to avoid bacterial contamination and thereby minimize the risk of any influence on research caused by microbiological organisms. Health monitoring is carried out twice a year, based on FELASA (Federation of European Laboratory Animal Science Associations) guidelines.

The domestic pig

The domestic pigs used in the present study are all crossbreeds between Danish Landrace and Yorkshire. This combination is normally used for production of mother animals. Danish Landrace is a long and lean breed known for high fertility and good motherhood. The litter size is around 13.5 piglets per litter. The Yorkshire pig has a high lean meat percentage, high daily gain and good meat quality. The mothering characteristics are good, with a litter size of around 12 piglets per litter. The National Committee for pig production mainly organizes The Danish Pig breeding program.

The newborn Landrace/Yorkshire domestic pig has a body mass of∼1.3–1.9 kg. Boars become sexually mature at an age of 6 months,while sows are sexually mature at an age of 7 months. The gestational period is 114 days and the average litter size is 10–14 animals. The adult mass of the domestic pig, at an age of two years, is 200–300 kg(Bollen et al., 2000). In Denmark, most of the domestic pigs are SPF-tested (Specific Pathogen Free) but are not submitted to the same strict control procedures as the Göttingen minipigs.

Experimental set-up

A total of 10 Göttingen minipigs and 12 domestic pigs were used in the study. All Göttingen minipigs [five neonate polts (1–2 days old)and five adult sows (1.5–3 years old)] were perfusion fixated by a procedure approved by the Danish Animal Research Inspectorate. The pigs were anesthetized with an intramuscular injection (1 ml per 10 kg body mass) of a mixture of 6.5 ml Narcoxyl®Vet (20 mg ml–1; Intervet,Denmark), 1.5 ml Ketaminol®Vet (100 mg ml–1; Intervet,Denmark) and 2.5 ml Methadone DAK (19 mg ml–1; Nycomed,Denmark) added to one bottle of Zoletil®50 Vet (Virbac Laboratories,France) without additional solvent. The pigs were then placed in a supine position on the operation table and supplied with a lethal dose of Pentobarbital (1 ml kg–1 body mass, 200 mg ml–1) before intervention. The deeply `anesthetized' pigs received a midsternal incision followed by a sternal split. The left cardiac ventricle was punctured by an infusion cannula, the right auricle was cut open, and a transcardial perfusion with 0.5–2.5 l saline followed by 2.5–7.5 l of 4% paraformaldehyde in 0.15 mol l–1Sørensen phosphate buffer (pH 7.4) was executed. The procedure took on average 10–15 min. The brain was removed and postfixed for an additional 24 h at 5°C in 1% paraformaldehyde.

The domestic pigs [six neonate sows (1 day old) and six adult sows(3–4 years old)] were euthanized with a lethal injection of Pentobarbital i.v. or with carbon dioxide in a Danish pig slaughterhouse. Following death, the brains were carefully removed from the skull and fixed by immediate immersion in 4% formaldehyde (0.15 mol l–1phosphate buffer, pH 7.4) for at least two weeks.

Tissue processing

The cerebral hemispheres were divided medially through the corpus callosum,and the left or right hemisphere was selected randomly before dehydration and infiltration in paraffin (Fig. 1A,B). The hemispheres were placed in a container with the midsagittal surface facing down, embedded in blocks of paraffin and sectioned into coronal series of 40 μm-thick sections on a Leica microtome(Fig. 1B,C). Satisfactory adhesion of the thick sections to glass slides was obtained using Superfrost Plus glass slides primed with a chromalun–gelatin solution(Feinstein et al., 1996). A predetermined fraction of the 40 μm-thick sections was sampled using moist filter paper and carefully rolled onto 40°C preheated primed slides covered by a thin layer of distilled water. The sections were dried overnight at 37°C and a final selection of slides was subsequently stained with Giemsa.

Stereological techniques

Estimates of neocortical cell numbers were obtained using the optical fractionator method (Fig. 1). The optical fractionator combines two stereological principles: the optical disector for counting particles and the systematic uniform sampling scheme of the fractionator (Gundersen,1986; West et al.,1991). The optical fractionator involves counting particles with optical disectors in a uniform random systematic sample that constitutes a known volume fraction of the region being analyzed. This is achieved by counting cells on a known fraction of sections (ssf; Fig. 1D), under a known fraction of the sectional area of the region (asf; Fig. 1E) in a known fraction of the thickness of a section (hsf; Fig. 1G). The total numbers of neocortical cells are estimated by multiplying the numbers of cells counted(∑Q) with the inverse sampling fractions. Finally, the bilateral cell number can be estimated by multiplying the unilateral number by two.
\[\ N_{\mathrm{total}}=(1{/}ssf)(1{/}asf)(1{/}hsf){\ }{\Sigma}Q^{-}{\ }2.\]
(1)
Fig. 1.

The optical fractionator sampling scheme. The cerebral hemispheres were divided medially (A), and the left or right hemisphere was selected systematically randomly before dehydration and infiltration in paraffin (B). The hemispheres were sectioned exhaustively into 40 μm-thick coronal sections (C), from which a predetermined fraction, the section sampling fraction (ssf), was sampled systematically randomly (D) and subsequently stained with Giemsa. Optical disectors were positioned systematically randomly on each of the sampled sections at regular predetermined x,y-positions (E,F). The area of the counting frame of the disector to the area associated with each x,y-step represents the area sampling fraction (asf). Counts were performed in all optical disectors hitting the structure of interest (E,F). Cellular nuclei were counted by moving the counting frame through a continuous stack of thin optical planes inside the section (G). The height of the disector(hdis) relative to the height of the section (t)represents the height sampling fraction (hsf). The disector is protected by upper and lower guard zones.

Fig. 1.

The optical fractionator sampling scheme. The cerebral hemispheres were divided medially (A), and the left or right hemisphere was selected systematically randomly before dehydration and infiltration in paraffin (B). The hemispheres were sectioned exhaustively into 40 μm-thick coronal sections (C), from which a predetermined fraction, the section sampling fraction (ssf), was sampled systematically randomly (D) and subsequently stained with Giemsa. Optical disectors were positioned systematically randomly on each of the sampled sections at regular predetermined x,y-positions (E,F). The area of the counting frame of the disector to the area associated with each x,y-step represents the area sampling fraction (asf). Counts were performed in all optical disectors hitting the structure of interest (E,F). Cellular nuclei were counted by moving the counting frame through a continuous stack of thin optical planes inside the section (G). The height of the disector(hdis) relative to the height of the section (t)represents the height sampling fraction (hsf). The disector is protected by upper and lower guard zones.

The optical disector may be considered as a three-dimensional probe,generated with the aid of a microscope with a high numerical aperture (NA=1.4)and an oil immersion objective, in which it is possible to observe thin focal planes in relatively thick sections (Fig. 1G). A counting frame with `exclusion' and `inclusion' lines is superimposed on the magnified image of the tissue on a computer screen, and the orientation in the z-axis is measured with a digital microcator with a precision of 0.5 μm. The purpose of `exclusion' and `inclusion'lines of the counting frame is to exclude edge effects arising from sub-sampling (Gundersen,1978). Similarly, upper and lower guard zones protect the counting frame to prevent bias as a consequence of loss of cells close to the section surfaces (Fig. 1G). All cells that come into focus within the frame and are not in focus at the uppermost position are counted as the focal plane is moved through the section. All cell counts were obtained using a BH-2 Olympus microscope and CAST-GRID software(Olympus, Ballerup, Denmark).

Counting criteria

The cells were identified as neurons if they had a combination of dendritic processes, a Giemsa-positive cytoplasm, a clearly defined nucleus with a pale and homogeneous nucleoplasm and a dark and condense centrally located nucleolus. The nucleus was used as the counting item, and around 250 neuronal nuclei were counted per brain(

\({\Sigma}\mathrm{Q}_{\mathrm{neu}}^{-}\)
⁠; Table 1). The glial cells were usually smaller and identified by the absence of a Giemsa-positive cytoplasm,the presence of heterochromatin clumps within the ovoid or irregularly shaped nucleus and the lack of a clearly identifiable nucleolus. Also here, the nucleus was used as the counting item, and an average of 360 glial nuclei were counted per brain(
\({\Sigma}\mathrm{Q}_{\mathrm{glial}}^{-}\)
; Table 1). No differentiation was made between astrocytes, oligodendrocytes or microglia. Endothelial cells were easily recognized by their dark and elongated nucleus and were excluded from all counts.

Table 1.

The fractionator sampling parameters, mean values

BA (μm)1/ssf∑secta(frame) (μm2)Step size (x,y) (μm)1/asf∑dishdis (μm)tQ- (μm)1/hsf
\({\Sigma}\mathrm{Q}_{\mathrm{neu}}^{-}\)
\({\Sigma}\mathrm{Q}_{\mathrm{glial}}^{-}\)
Göttingen minipig             
    Neonate 40 50 11-13 749 1600 3417 172 15 39.1 2.6 253 383 
    Adult 40 100 11-12 2497 2600 2707 168 20 37.7 1.9 289 337 
Domestic pig             
    Neonate 40 100 8-9 999 1800 3242 141 15 38.4 2.6 260 
    Adult 40 200 7-8 2706 2800 2896 137 25 37.9 1.5 235 
BA (μm)1/ssf∑secta(frame) (μm2)Step size (x,y) (μm)1/asf∑dishdis (μm)tQ- (μm)1/hsf
\({\Sigma}\mathrm{Q}_{\mathrm{neu}}^{-}\)
\({\Sigma}\mathrm{Q}_{\mathrm{glial}}^{-}\)
Göttingen minipig             
    Neonate 40 50 11-13 749 1600 3417 172 15 39.1 2.6 253 383 
    Adult 40 100 11-12 2497 2600 2707 168 20 37.7 1.9 289 337 
Domestic pig             
    Neonate 40 100 8-9 999 1800 3242 141 15 38.4 2.6 260 
    Adult 40 200 7-8 2706 2800 2896 137 25 37.9 1.5 235 

Abbreviations: BA, block advance; ssf, section sampling fraction;∑sect, number of sections; a(frame), area of counting frame;∑dis, number of disectors hitting neocortex; hdis,disector height; tQ-, q-weighted section thickness; hsf, height sampling fraction;

\({\Sigma}\mathrm{Q}_{\mathrm{neu}}^{-}\)
⁠, number of neurons counted;
\({\Sigma}\mathrm{Q}_{\mathrm{glial}}^{-}\)
, number of glial cells counted.

The a(frame) for glial cell estimations in adult Göttingen minipigs was 50% of the original frame.

The section sampling fraction

A sampling scheme was designed based on systematic uniform random sampling(SURS) to ensure that all parts of the neocortex had an equal probability of being sampled. Based on a pilot study, every 50th or 100th section was sampled from the neonate and adult Göttingen minipig brain, respectively, whereas every 100th or 200th section was sampled from the domestic pig(Table 1). The first section was selected randomly using a random number within the sampled section period. In the event of poor technical quality, e.g. due to the presence of folds or breaks, the subsequent section was sampled instead. Extra sections were also sampled to circumvent potential damage of sections during staining procedures. These deviations are of no consequence to the results, provided the overall sampling scheme is maintained. The section fraction represents the section sampling fraction (ssf) and provided 11–13 sections for final quantitative analyses in the Göttingen minipig and 7–9 sections in the domestic pig.

The area sampling fraction

In each of the sampled sections, counts of neurons or glial cells were made with optical disectors at regular predetermined x,y positions in the neocortex. The neocortex was defined as the isocortex and mesocortex. Again,the first disector was positioned randomly within the first x,yinterval by the CAST-GRID software. The area of the counting frame of the disector, a(frame), is known relative to the area associated with each x,y-step. The area of the sampling fraction (asf) is accordingly:
\[\ asf=a(frame){\ }{/}{\ }a(x,y\mathrm{-step}).\]
(2)
In a pilot study, the step length between disectors was adjusted to provide approximately 150 disectors in each brain. Subsequently, the frame area and height of the disector were adjusted to obtain 1–2 particles per disector (see below).

The height sampling fraction

The height of the disector (hdis) should be known relative to the thickness (t) of the section. To prevent bias as a consequence of loss of cells, the disector is guarded by upper and lower guard zones. The size of the guard zones depends upon the size of the particles counted. For neocortical neurons it should usually be at least 5 μm in the top and bottom of the section. To compensate for potential deformation of the sections in the z-axis, the height sampling fraction (hsf)depends on the Q weighted mean section thickness(Q):
\[\ hsf=h_{\mathrm{dis}}{\ }{/}{\ }{\bar{t}}_{\mathrm{Q}^{-}},\]
(3)
\[\ {\bar{t}}_{\mathrm{Q}^{-}}={\Sigma}(t_{\mathrm{i}}q_{\mathrm{i}}^{-}){\ }{/}{\ }{\Sigma}(q_{\mathrm{i}}^{-}),\]
(4)
where ti is the local section thickness centrally in the ith counting frame with a disector count of
\(q_{{\bar{\mathrm{i}}}}\)
(Dorph-Petersen et al.,2001).

Volume estimates

The systematic uniform random placement of disectors in the fractionator design was used to estimate volumes in accordance with the unbiased principles of the Cavalieri-estimator:
\[\ V_{\mathrm{neo}}=2{\ }{\Sigma}P{\ }a(p){\ }t{\ }k,\]
(5)
where ∑P is the number of frame upper-corner-points hitting the neocortex, a(p) is the x,y-step area, t is the block or microtome advance and k is the section sampling fraction. Note that estimated volumes are not used for the estimation of total cell number; the fractionator estimates of total cell number are independent of the containing volume and its shrinkage and deformation.

Statistical analyses and estimation of precision

The differences between means were analyzed using an unpaired two-tailed Student's t-test. Group variability is shown in parentheses as the coefficient of variation (CV=s.d./mean)(Table 2). The precision of the estimate of the total cell number in each subject was estimated as the coefficient of error (CE=s.e.m./mean), caused by sampling error related to the counting noise, the systematic uniform random sampling and variances in section thickness (Table 3). The precision of an individual estimate is related to the uniformity of the distribution of particles being counted and the amount of sampling that has been performed (Table 3). The sampling was considered optimal when the observed variance of the individual estimate, CE2, was less than half the observed interindividual variance, CV2.

Table 2.

Litter parameters and major estimated quantities

NBody mass (kg)Brain mass (g)Neocortex volume (cm3)Mean CE volumeNeocortical neurons (×106)Mean CE neuronsNeocortical glial cells (×106)Mean CE glial cells
Göttingen minipig          
    Neonate 0.56 (0.081) 27.8 (0.037) 1.75 (0.057) 0.049 252.5 (0.11) 0.065 381.9 (0.10) 0.053 
    Adult 37.9 (0.11) 79.0 (0.077) 9.01 (0.11) 0.052 323.8 (0.093) 0.064 714.2 (0.12) 0.057 
    P-value  — — — — 0.01 — <0.01 — 
Domestic pig          
    Neonate 214.5 (0.30) 29.4 (0.16) 3.64 (0.098) 0.060 424.8 (0.13) 0.067 — — 
    Adult 212-217 134.0 (0.11) 17.2 (0.083) 0.051 432.1 (0.046) 0.068 — — 
    P-value  — — — — 0.76 — — — 
NBody mass (kg)Brain mass (g)Neocortex volume (cm3)Mean CE volumeNeocortical neurons (×106)Mean CE neuronsNeocortical glial cells (×106)Mean CE glial cells
Göttingen minipig          
    Neonate 0.56 (0.081) 27.8 (0.037) 1.75 (0.057) 0.049 252.5 (0.11) 0.065 381.9 (0.10) 0.053 
    Adult 37.9 (0.11) 79.0 (0.077) 9.01 (0.11) 0.052 323.8 (0.093) 0.064 714.2 (0.12) 0.057 
    P-value  — — — — 0.01 — <0.01 — 
Domestic pig          
    Neonate 214.5 (0.30) 29.4 (0.16) 3.64 (0.098) 0.060 424.8 (0.13) 0.067 — — 
    Adult 212-217 134.0 (0.11) 17.2 (0.083) 0.051 432.1 (0.046) 0.068 — — 
    P-value  — — — — 0.76 — — — 

Mean

\(\mathrm{CE}=\sqrt{\mathrm{mean}(\mathrm{CE}^{2})}\)
⁠. For calculation of CE for cell number and volume, see Table 3. The inter-individual variability in each group is shown in parentheses as coefficient of variability (CV=s.d./mean). Volume estimates are from shrunken paraffin-embedded tissue and should not be compared with fresh brain volumes.

Table 3.

How to estimate the coefficient of error (CE) of the estimate in an individual animal

Section no.Q-A
\(\mathrm{Q}_{\mathrm{i}}^{-}{\times}\mathrm{Q}_{\mathrm{i}}^{-}\)
B
\(\mathrm{Q}_{\mathrm{i}}^{-}{\times}\mathrm{Q}_{\mathrm{i}+1}^{-}\)
C
\(\mathrm{Q}_{\mathrm{i}}^{-}{\times}\mathrm{Q}_{\mathrm{i}+2}^{-}\)
14 196 224 308 
16 256 352 288 
22 484 396 506 
18 324 414 216 
23 529 276 552 
12 144 288 312 
24 576 624 528 
26 676 572 884 
22 484 726 748 
10 33 1089 1122 462 
11 34 1156 476  
12 14 196   
Σ 258 6110 5470 4778 
VARNOISEa 258    
VARSURSb 1.89    
CE(t)c 0.007    
CE(ΣQ-)d 0.062    
CE(N)e 0.063    
Section no.Q-A
\(\mathrm{Q}_{\mathrm{i}}^{-}{\times}\mathrm{Q}_{\mathrm{i}}^{-}\)
B
\(\mathrm{Q}_{\mathrm{i}}^{-}{\times}\mathrm{Q}_{\mathrm{i}+1}^{-}\)
C
\(\mathrm{Q}_{\mathrm{i}}^{-}{\times}\mathrm{Q}_{\mathrm{i}+2}^{-}\)
14 196 224 308 
16 256 352 288 
22 484 396 506 
18 324 414 216 
23 529 276 552 
12 144 288 312 
24 576 624 528 
26 676 572 884 
22 484 726 748 
10 33 1089 1122 462 
11 34 1156 476  
12 14 196   
Σ 258 6110 5470 4778 
VARNOISEa 258    
VARSURSb 1.89    
CE(t)c 0.007    
CE(ΣQ-)d 0.062    
CE(N)e 0.063    

The precision, CE, of a fractionator estimate is a function of three independent factors: the noise variance (VARNOISE), the systematic uniform random sampling variance (VARSURS) and the variance attributable to variations in section thickness [CE(t)](see also West and Gundersen, 1990; West et al., 1991; West et al., 1996).

From the calculations above (a—e), it is obvious that the CE(N) for all practical purposes may be considered a function of the NOISE variance. If one wants to reduce the CEs of the individual estimates further, this would, in this example, be achieved most effectively by sampling more on the sections already in the series.

a

VARNOISE is the uncertainty' in the estimate that comes from disector counts within a section and is equal toΣQ-.

b

VARSURS is the `uncertainty' in the estimate that arises due to sampling between sections, i.e. because repeated estimates based on different sets of sections may vary. The VARSURS is calculated using a prediction model that takes into account the systematic nature of the sampling (Gundersen et al.,1999): VARSURS= [3(A—Noise)-4B+C]/240. For calculation of A, B and C, see columns above. When one uses more than 5-10 sections for a biological estimate, the SURS variance is usually negligible relative to the Noise variance. The denominator is a constant(Gundersen et al., 1999).

c

CE(t)=[s.d.()/(√n)](1/),where is the mean section thickness in each section measured with the digital microcator, n is the number of sections and is the mean section thickness between all sections. The CE(t) usually contributes less than 1% to the overall estimator variance and can be ignored in most studies where section homogeneity is high.

d

The total sampling variance, CE(ΣQ-), is calculated as:

\(\sqrt{\mathrm{VAR}_{\mathrm{NOISE}}+\mathrm{VAR}_{\mathrm{SURS}}}{/}{\Sigma}\mathrm{Q}^{-}\)
⁠.

e

The total CE for the final estimate, CE(N), is eventually calculated as:

\(\sqrt{\mathrm{CE}^{2}({\Sigma}\mathrm{Q}^{-})+\mathrm{CE}^{2}(t)}\)
⁠.

Fig. 2.

The postnatal development of neocortical cell numbers in the Göttingen minipig and the domestic pig brain. Total neuron number (A) and total glial cell number (B) in the Göttingen minipig. Total neuron number in the domestic pig (C).

Fig. 2.

The postnatal development of neocortical cell numbers in the Göttingen minipig and the domestic pig brain. Total neuron number (A) and total glial cell number (B) in the Göttingen minipig. Total neuron number in the domestic pig (C).

The total number of neocortical neurons in the Göttingen minipig brain increases from ∼253 million at birth to ∼324 million in adulthood(Table 2, Fig. 2A). This significant(P=0.01) 28% difference demonstrates a pronounced postnatal development of neurons in the Göttingen minipig brain. A significant postnatal development is also observed for neocortical glial cells(P<0.01), increasing from ∼382 million in the neonate to∼714 million glial cells in the adult Göttingen minipig, an 87%difference (Table 2; Fig. 2B). The glial-to-neuron ratio changes accordingly from 1.5 to 2.2. The total brain mass increases almost threefold from a mean of 27.8 g at birth to 79.0 g as adult. Meanwhile,the proportional increase of neocortex volume is more than 500%(Table 2).

A corresponding postnatal development of neocortical neurons is not observed in the domestic pig brain. The domestic pig has ∼425 million neocortical neurons at birth and ∼432 million in adulthood(Table 2; Fig. 2C), which is ∼35%more than in the adult Göttingen minipig neocortex (P<0.01). The number of neocortical glial cells was not estimated. The total brain mass increases from a mean of 29.4 g at birth to 134.0 g as adult. The postnatal proportional increase in neocortex volume is, by comparison, close to 500%(Table 2).

With the principal aim to evaluate the pig as a potential animal model for human brain development, a comparative study has been performed to evaluate the postnatal development of neocortical cell numbers in the experimental Göttingen minipig and the domestic pig brain. It is demonstrated that the cellular development in the Göttingen minipig brain is incomplete at term and that both neocortical neuron and glial cell numbers increase significantly from time of birth to adulthood. By contrast, the neocortical development of neurons in the domestic pig seems to be fully established at birth. An explanation for this strain difference is not at hand.

A rapid development of the brain prior to the general body growth seems to be characteristic for all mammals. Differences between species are related to the gestational time of major cellular multiplication, the developmental stage at birth and the complexity of the final product(Dobbing, 1974). At the time of birth, the human brain constitutes around 23% of its adult mass, and in the first two postnatal years it increases rapidly to around 75% of its adult mass(Dobbing, 1974). In comparison, the brain of the closely related macaque monkey constitutes almost 65% of its adult mass at birth (Dobbing,1974). The pig brain constitutes around 25% of its adult mass at birth and seems in this aspect to be more similar to that of humans(Dickerson and Dobbing, 1966). The pig is also considered a perinatal brain developer, like human, and several studies have shown a good correspondence to the developing human brain with respect to myelination, compositions and electrical activity(Dickerson and Dobbing, 1966;Pampliglione, 1971; Fang et al.,2005; Flynn, 1984; Thibault and Margulies, 1998). Accordingly, the pig has been considered an appropriate model for human brain development (Dickerson and Dobbing,1966; Book and Bustad,1974; Pond et al.,2000).

The present study demonstrates different developmental patterns in the neocortex of two strains of pigs. A significant postnatal development of neuron and glial cells is observed in the Göttingen minipig, whereas the adult number of neurons is established at birth in the domestic pig. This difference in cell number is not represented by a corresponding difference in relative brain mass of the neonates. The brain constitutes 35% of its adult mass in the Göttingen minipig and 22% in the domestic pig, whereas the relative growth of neocortical volume is close to 500% in both strains. The differences can by explained neither by differences in the stereological sampling or counting procedures (see below) nor the different tissue processing. Both fixation methods provided adequate histological sections with no clear difference in cellular morphology. No clinical data corroborate the observed difference; the gestational period is the same and birth mass and behavior is more or less identical in both strains. However, the results do substantiate that strain differences should be considered in future experimental studies using the pig brain and exemplifies the need for a full designation of the specific strain used.

Based on the results, the domestic pig seems to be a more proper model for evaluating the effects of developmental insults on the human brain than the Göttingen minipig. Even though recent developments have made it possible to unambiguously demonstrate that new neurons are added to the adult primate neocortex as well as to other mammal brains(Gould et al., 1999; Gould et al., 2001; Gould and Gross, 2002), it is still generally accepted that most neurogenesis in monkeys(Rakic and Sidman, 1968; Sidman and Rakic, 1973; Rakic, 1974; Rakic, 1978; Rakic, 1985a; Rakic, 1985b; Rakic, 1988) and humans(Dobbing and Sands, 1973; Dobbing, 1974; Samuelsen et al., 2003) is completed at midgestation or at least before term. The rate of neurons added to the neocortex in adulthood has no relative influence on the total neocortical neuron number when compared with the rate of multiplication during early development (Gould et al.,2001). Furthermore, recent stereological results from our laboratory reveal that the total number of neurons in the cortical plate of human newborns equals the total number in adults(Larsen et al., in press). These results definitively refute previous results demonstrating a major postnatal neurogenesis in humans (Shankle et al., 1998; Shankle et al.,1999). Regrettably, glial cells were not estimated in the domestic pig. Results from one neonate and one adult domestic pig (not presented) do indicate a prenatal development in glial cell number similar to the postnatal increase observed in the Göttingen minipig. Further studies on several postnatal ages are considered valuable in order to describe the growth slope of neuronal increase in the Göttingen minipig.

The Göttingen minipig and the domestic pig have also been considered as useful non-primate models for a number of human neurological diseases(McClellan, 1968; Douglas, 1972). However,although several neuroanatomical studies have been performed, the documentation of pig brain anatomy, connectivity and function is still incomplete. Quantitative information of cell numbers based on systematic sampling procedures has been limited to subcortical areas, e.g. the hippocampus (Holm and West,1994) and the subthalamic nucleus(Larsen et al., 2004). Here,we present the total number of neocortical neurons in two strains of pigs. In the adult Göttingen minipig, the neocortex contains 324 million neurons,whereas the domestic pig brain contains 432 million neocortical neurons. This 33% strain difference in neocortical neuron number was not an unexpected finding considering the general relationship between body size and neuron number (for a review, see Williams and Herrup, 1988). The total number of neocortical neurons has also been estimated in a number of other species. There are ∼3 million neocortical neurons in the mouse brain(Bonthius et al., 2004), 21 million in the rat (Korbo et al.,1990), 12.8 billion in the minke whale (Nina Eriksen and Bente Pakkenberg, unpublished data) and 19–23 billion in the human neocortex(Pakkenberg and Gundersen,1997). When compared to these species, one of the most valuable findings in the pig brain is the rather low coefficient of variation (CV) from which the true biological variance can be estimated to be less than 10% (see the additivity of variances below). The low biological variance seems to be a reproducible finding for the pig brain(Holm and West, 1994; Larsen et al., 2004; Jelsing et al., 2005a; Jelsing et al., 2005b) and supports the continued use of pigs in neurotoxicologic studies, since perturbations may be detected with great sensitivity.

Stereological design

Estimates of neocortical cell numbers were obtained using the optical fractionator method, which is efficient and independent of any tissue shrinkage or expansion that may take place during any stage of tissue preparation (Gundersen, 1986; West et al., 1991). It was considered optimal on paraffin sections in which shrinkage during processing is significant. Provided that mounted sections maintain a sufficient depth to accommodate optical disectors, the optical fractionator can also be applied on other preparations, e.g. cryostat or vibratome sections.

Sampling parameters were optimized to obtain a high efficiency in terms of precision and effort. A total count of ∼100–150 cells in 75–100 disectors distributed systematically randomly on 5–10 sections is usually enough to obtain an estimate with a precision appropriate for most biological structures (Pakkenberg and Gundersen, 1988; Korbo et al., 1990; West,1993). However, because the variance in the experimental groups appeared to be rather low, it was decided to increase the efficiency of the fractionator by intensifying the sampling of cells to an average of 250 particles per brain in 140–170 disectors. An additional number of sections were sampled from the Göttingen minipig. This was done to ensure enough sections for forthcoming studies of specific glial subpopulations, of which some have a more heterogeneous distribution in the brain. The efficiency of the fractionator sampling was evaluated from variance analysis of relative variances and estimator CE. These two measures are related through the basic equation (the additivity of variances):
\[\ \mathrm{CV}_{\mathrm{obs}}^{2}=\mathrm{CV}_{\mathrm{biol}}^{2}+\mathrm{CE}^{2}.\]
(6)
The sampling was considered optimal when the variance of the individual estimate, CE2, was less than half the observed interindividual variance, CV2, because the real inherent biological variability in the cohorts of pigs (CV2biol) then contributes most to the observed relative variance. Subsequently, the estimator CE was optimized in relation to variances from sampling with disectors within a section and sampling between sections (Table 3).

Even though the efficiency and mathematical unbiasedness of the optical fractionator method for estimating cells in the pig brain neocortex are evident, some fundamental requirements have to be fulfilled for a proper application (West, 1993). A first requirement is that the whole structure is accessible; secondly, one must be confident that all cells of interest can be identified unambiguously and that penetration of staining is complete throughout the thickness of the section. In the present study, a complete penetration of staining was optimized beforehand by registering the z-distribution of all cells. However, difficulty in distinguishing glial cells from neurons appeared when stained with the modified Giemsa method. This was especially true in the neonate brains, where the neuronal density was high and morphometric differences between cells were less pronounced. The lack of clear criteria for distinguishing neurons and glial cells in cortical regions has previously been a major problem for stereologists(Braendgaard et al., 1990; Davanlou and Smith, 2004) and may partly explain the somewhat lower hemispheric cell number published in a screening procedure of young Göttingen minipig(Jelsing et al., 2005a). Recently, a number of immunohistochemical markers have been evaluated in the pig brain (Lyck et al., 2006),and the combination of immunohistochemistry and stereology may provide a better approach for quantifying neurons and glial cell populations in future quantitative studies of the pig brain.

The authors would like to thank Skaerbaek Pig Slaughterhouse (Skaerbaek,Denmark) for donating the brains from adult domestic pigs. Thanks to the technical staff at the Research Laboratory for Stereology and Neuroscience,especially Susanne Sørensen and Hans Jørgen Jensen. The study received financial support from the Lundbeck Foundation and the Gangsted Foundation.

Bollen, P. and Ellegaard, L. (
1997
). The Göttingen minipig in pharmacology and toxicology.
Pharmacol. Toxicol.
80
,
3
-4.
Bollen, P. J. A., Hansen, A. K. and Rasmussen, H. J.(
2000
).
The Laboratory Swine.
Boca Raton:CRC Press.
Bonthius, D. J., McKim, R., Koele, L., Harb, H., Karacay, B.,Mahoney, J. and Pantazis, N. J. (
2004
). Use of frozen sections to determine neuronal number in the murine hippocampus and neocortex using the optical disector and optical fractionator.
Brain Res. Brain Res. Protoc.
14
,
45
-57.
Book, S. A. and Bustad, L. K. (
1974
). The fetal and neonatal pig in biomedical research.
J. Anim. Sci.
38
,
997
-1002.
Bourgeois, J. P., Goldman-Rakic, P. S. and Rakic, P.(
1994
). Synaptogenesis in the prefrontal cortex of rhesus monkeys.
Cereb. Cortex
4
,
78
-96.
Braendgaard, H., Evans, S. M., Howard, C. V. and Gundersen, H. J. (
1990
). The total number of neurons in the human neocortex unbiasedly estimated using optical disectors.
J. Microsc.
157
,
285
-304.
Damm Jorgensen, K. (
1998
). Teratogenic activity of tretionin in the Göttingen minipig.
Scand. J. Lab. Anim. Sci.
25
,
235
-243.
Davanlou, M. and Smith, D. F. (
2004
). Unbiased stereological estimation of different cell types in rat cerebral cortex.
Image Anal. Stereol.
23
,
1
-12.
Dickerson, J. W. T. and Dobbing, J. (
1966
). Prenatal and postnatal growth and development of the central nervous system of the pig.
Proc. R. Soc. Lond. B Biol. Sci.
166
,
384
-395.
Dobbing, J. (
1974
). The later development of the brain and its vulnerability. In
Scientific Foundation of Paediatrics
(ed. J. A. Davis and J. Dobbing), pp.
565
-577. Philadelphia: Saunders.
Dobbing, J. and Sands, J. (
1973
). Quantitative growth and development of human brain.
Arch. Dis. Child
.
48
,
757
-767.
Dobbing, J. and Sands, J. (
1979
). Comparative aspects of the brain growth spurt.
Early Hum. Dev.
3
,
79
-83.
Dorph-Petersen, K. A., Nyengaard, J. R. and Gundersen, H. J.(
2001
). Tissue shrinkage and unbiased stereological estimation of particle number and size.
J. Microsc.
204
,
232
-246.
Douglas, W. R. (
1972
). Of pigs and men and research: a review of applications and analogies of the pig, Sus Scrofa, in human medical research.
Space Life Sci.
3
,
226
-234.
Fang, M., Li, J., Gong, X., Antonio, G., Lee, F., Kwong, W. H.,Wai, S. M. and Yew, D. T. (
2005
). Myelination of the pig's brain: a correlated MRI and histological study.
Neurosignals
14
,
102
-108.
Feinstein, R. E., Westergren, E., Bucht, E., Sjoberg, H. E. and Grimelius, L. (
1996
). Estimation of the C-cell numbers in rat thyroid glands using the optical fractionator.
J. Histochem. Cytochem.
44
,
997
-1003.
Flynn, T. J. (
1984
). Developmental changes of myelin-related lipids in brain of miniature swine.
Neurochem. Res.
9
,
935
-945.
Glodek, P. (
1986
). Breeding program and population standards of the Göttingen miniature swine. In
Swine in Biomedical Research
(ed. M. E. Tumbleson),pp.
23
-28. New York: Plenum Press.
Gould, E. and Gross, C. G. (
2002
). Neurogenesis in adult mammals: some progress and problems.
J. Neurosci.
22
,
619
-623.
Gould, E., Reeves, A. J., Graziano, M. S. and Gross, C. G.(
1999
). Neurogenesis in the neocortex of adult primates.
Science
286
,
548
-552.
Gould, E., Vail, N., Wagers, M. and Gross, C. G.(
2001
). Adult-generated hippocampal and neocortical neurons in macaques have a transient existence.
Proc. Natl. Acad. Sci. USA
98
,
10910
-10917.
Gundersen, H. J. (
1978
). Estimators of the number of objects per area unbiased by edge effects.
Microsc. Acta
81
,
107
-117.
Gundersen, H. J. (
1986
). Stereology of arbitrary particles. A review of unbiased number and size estimators and the presentation of some new ones, in memory of William R. Thompson.
J. Microsc.
143
,
3
-45.
Gundersen, H. J. and Jensen, E. B. (
1987
). The efficiency of systematic sampling in stereology and its prediction.
J. Microsc.
147
,
229
-263.
Gundersen, H. J., Jensen, E. B., Kieu, K. and Nielsen, J.(
1999
). The efficiency of systematic sampling in stereology– reconsidered.
J. Microsc.
193
,
199
-211.
Holm, I. E. and West, M. J. (
1994
). Hippocampus of the domestic pig: a stereological study of subdivisional volumes and neuron numbers.
Hippocampus
4
,
115
-125.
Jelsing, J., Olsen, A. K., Cumming, P., Gjedde, A., Hansen, A. K., Arnfred, S., Hemmingsen, R. and Pakkenberg, B.(
2005a
). A volumetric screening procedure for the Göttingen minipig brain.
Exp. Brain Res.
162
,
428
-435.
Jelsing, J., Rostrup, E., Markenroth, K., Paulson, O. B.,Gundersen, H. J., Hemmingsen, R. and Pakkenberg, B.(
2005b
). Assessment of in vivo MR imaging compared to physical sections in vitro – a quantitative study of brain volumes using stereology.
Neuroimage
26
,
57
-65.
Jones, N. A., Field, T., Davalos, M. and Hart, S.(
2004
). Greater right frontal EEG asymmetry and nonemphathic behavior are observed in children prenatally exposed to cocaine.
Int. J. Neurosci.
114
,
459
-480.
Kjellmer, I., Thordstein, M. and Wennergren, M.(
1992
). Cerebral function in the growth-retarded fetus and neonate.
Biol. Neonate
62
,
265
-270.
Korbo, L., Pakkenberg, B., Ladefoged, O., Gundersen, H. J.,Arlien-Soborg, P. and Pakkenberg, H. (
1990
). An efficient method for estimating the total number of neurons in rat brain cortex.
J. Neurosci. Methods
31
,
93
-100.
Larsen, C. C., Larsen, K. B., Bogdanovic, N., Laursen, H.,Græm, N., Samuelsen, G. B. and Pakkenberg, B. (in press). Total number of cells in the human newborn telencephalic wall.
Neuroscience
Larsen, M., Bjarkam, C. R., Ostergaard, K., West, M. J. and Sorensen, J. C. (
2004
). The anatomy of the porcine subthalamic nucleus evaluated with immunohistochemistry and design-based stereology.
Anat. Embryol. (Berl)
.
208
,
239
-247.
Lyck, L., Jelsing, J., Jensen, P. S., Lambertsen, K. L.,Pakkenberg, B. and Finsen, B. (
2006
). Immunohistochemical visualization of neurons and specific glial cells for stereological application in the porcine neocortex.
J. Neurosci. Methods
152
,
229
-242.
McClellan, R. O. (
1968
). Applications of swine in biomedical research.
Lab. Anim. Care
18
,
120
-126.
Pakkenberg, B. and Gundersen, H. J. (
1988
). Total number of neurons and glial cells in human brain nuclei estimated by the disector and the fractionator.
J. Microsc.
150
,
1
-20.
Pakkenberg, B. and Gundersen, H. J. (
1997
). Neocortical neuron number in humans: effect of sex and age.
J. Comp. Neurol.
384
,
312
-320.
Pampiglione, G. (
1971
). Some aspects of development of cerebral function in mammals.
Proc. R. Soc. Med.
64
,
429
-435.
Pond, W. G., Boleman, S. L., Fiorotto, M. L., Ho, H., Knabe, D. A., Mersmann, H. J., Savell, J. W. and Su, D. R.(
2000
). Perinatal ontogeny of brain growth in the domestic pig.
Proc. Soc. Exp. Biol. Med.
223
,
102
-108.
Rakic, P. (
1974
). Neurons in rhesus monkey visual cortex: systematic relation between time of origin and eventual disposition.
Science
183
,
425
-427.
Rakic, P. (
1978
). Neuronal migration and contact guidance in the primate telencephalon.
Postgrad. Med. J
.
54
,
25
-40.
Rakic, P. (
1985a
). DNA synthesis and cell division in the adult primate brain.
Ann. N. Y. Acad. Sci.
457
,
193
-211.
Rakic, P. (
1985b
). Limits of neurogenesis in primates.
Science
227
,
1054
-1056.
Rakic, P. (
1988
). Specification of cerebral cortical areas.
Science
241
,
170
-176.
Rakic, P. and Sidman, R. L. (
1968
). Supravital DNA synthesis in the developing human and mouse brain.
J. Neuropathol. Exp. Neurol.
27
,
246
-276.
Rodning, C., Beckwith, L. and Howard, J.(
1989
). Prenatal exposure to drugs: behavioral distortions reflecting CNS impairment?
Neurotoxology
10
,
629
-634.
Samuelsen, G. B., Larsen, K. B., Bogdanovic, N., Laursen, H.,Graem, N., Larsen, J. F. and Pakkenberg, B. (
2003
). The changing number of cells in the human fetal forebrain and its subdivisions: a stereological analysis.
Cereb. Cortex
13
,
115
-122.
Shankle, W. R., Landing, B. H., Rafii, M. S., Schiano, A., Chen,J. M. and Hara, J. (
1998
). Evidence for a postnatal doubling of neuron number in the developing human cerebral cortex between 15 months and 6 years.
J. Theor. Biol.
191
,
115
-140.
Shankle, W. R., Rafii, M. S., Landing, B. H. and Fallon, J. H. (
1999
). Approximate doubling of numbers of neurons in postnatal human cerebral cortex and in 35 specific cytoarchitectural areas from birth to 72 months.
Pediatr. Dev. Pathol.
2
,
244
-259.
Sidman, R. L. and Rakic, P. (
1973
). Neuronal migration, with special reference to developing human brain: a review.
Brain Res.
62
,
1
-35.
Thibault, K. L. and Margulies, S. S. (
1998
). Age-dependent material properties of the porcine cerebrum: effect on pediatric inertial head injury criteria.
J. Biomech.
31
,
1119
-1126.
West, M. J. (
1993
). New stereological methods for counting neurons.
Neurobiol. Aging
14
,
275
-285.
West, M. J. and Gundersen, H. J. (
1990
). Unbiased stereological estimation of the number of neurons in the human hippocampus.
J. Comp. Neurol.
296
,
1
-22.
West, M. J., Slomianka, L. and Gundersen, H. J.(
1991
). Unbiased stereological estimation of the total number of neurons in the subdivisions of the rat hippocampus using the optical fractionator.
Anat. Rec.
231
,
482
-497.
West, M. J., Ostergaard, K., Andreassen, O. A. and Finsen,B. (
1996
). Estimation of the number of somatostatin neurons in the striatum: an in situ hybridization study using the optical fractionator method.
J. Comp. Neurol.
370
,
11
-22.
Williams, R. W. and Herrup, K. (
1988
). The control of neuron number.
Annu. Rev. Neurosci.
11
,
423
-453.