Specimens
We processed postmortem brain samples from 50 individuals across nine primate genera. Subjects included platyrrhines (cotton-top tamarins, common marmosets, owl monkeys, and tufted capuchins), cercopithecids (pig-tailed and rhesus macaques), African great apes (common chimpanzees and bonobos), and humans (see Table 1 for details). Sexes were balanced as much as possible. All individuals were adult, free of gross neuropathology. Nonhuman primates had been housed in accordance with their home institution’s animal care and use standards. We obtained the human brain specimens from the Northwestern University Alzheimer’s Disease Center Brain Bank. All individuals were non-geriatric and free of neuropathological disease. We received whole brains or sections for nonhuman primates from the National Chimpanzee Brain Resource, The Great Ape Neuroscience Project, the Oregon National Primate Research Center, the Washington National Primate Research Center, the New England Regional Primate Research Center (Harvard University), Alpha Genesis, and various zoos. Brain samples were fixed by immersion in 10% buffered formalin for 7-10 days, then transferred to a solution of 0.1 M phosphate-buffered saline (PBS, pH 7.4) containing 0.1% sodium azide and stored at 4º C until further processing.
Sample Processing
All samples were cryoprotected in a graded series of sucrose solutions (10%, 20%, and 30%) prior to sectioning. The brain samples were rapidly frozen with dry ice and then sectioned at 40 µm using a freezing sliding microtome (SM2000R, Leica, Chicago, IL). Each section was placed into an individual centrifuge tube containing freezer storage solution (30% each distilled water, ethylene glycol, and glycerol and 10% 0.244 M PBS), numbered sequentially, and then stored at -20° C until histological or immunohistochemical processing. Every tenth section was Nissl-stained to visualize the regions of interest for immunohistochemical staining and stereological quantification of neuron cell density.
Figures 1 and 2 show the two areas of interest for the present study: the NAcc, which comprises most of the ventral striatum and in this study we conservatively defined as the portion of the ventromedial striatum that is ventral and medial to the inferior border of the internal capsule in sections rostral to the appearance of the anterior commissure ; the VP, which was traced in sections ventral and adjacent to the anterior commissure with readily discernable borders (Paxinos et al. 2000; 2012, Ding et al. 2016; Haber and McFarland, 1999). We did not differentiate between the shell and core region of the NAcc in this study as TH does not provide clear boundaries between these regions (e.g., McCollum et al., 2016).
Immunohistochemistry
Sections that spanned the NAcc and VP for each subject were stained for tyrosine hydroxylase (TH), the rate-limiting enzyme for catecholamine synthesis, using the avidin-biotin-peroxidase method as described previously (Raghanti et al. 2016; Raghanti et al. 2009; Raghanti et al. 2008). Briefly, sections were pretreated for antigen retrieval by incubating in 0.05% citraconic acid (pH 7.4) at 86°C in a water bath for 30 minutes. Sections were then rinsed, and endogenous peroxidase was quenched using a solution of 75% methanol, 2.5% hydrogen peroxide (30%), and 22.5% distilled water for 20 minutes at room temperature. Sections were pre-blocked in a solution of 4% normal goat serum, 0.6% Triton X-100 detergent, 90.4% PBS, and 5% bovine serum albumin. Following this, sections were incubated in a rabbit anti-TH polyclonal antibody (Millipore, Bedford, MA, AB152, RRID 390204) diluted to 1:1,000 for 24 hours at room temperature followed by 24 hours at 4° C. Sections were then incubated in a biotinylated secondary antibody (1:200) followed by the avidin-peroxidase complex (PK-6100, Vector Laboratories, Burlingame, CA) for 1 hour at room temperature. A 3,3’ -diaminobenzidine-peroxidase (DAB) substrate with nickel solution enhancement was used as the chromogen (SK-4100, Vector Laboratories).
Data Acquisition
Quantitative analyses were performed using computer-assisted stereology on an Olympus BX-51 photomicroscope system equipped with a digital camera and StereoInvestigator software version 10 (MBF Bioscience, Williston, VT, USA). Subsampling techniques were performed for each species to determine appropriate sampling parameters. TH-immunoreactive (ir) axons were quantified using the SpaceBalls probe at 100x (N.A. 1.4) magnification under Koehler illumination using a hemisphere set at 7µm with a 2% guard zone. Section thickness was measured at every fifth sampling site with an average mounted section thickness for immunostained sections being 23.9 ± 8.2 µm. Sampling grids were optimized for each specimen, with an average grid of 496 x 401 µm (range: 106 x 66 – 1632 x 1092). Axons were marked where they intersected the outline of the hemisphere. Axon length density (ALv) was calculated as the total fiber length divided by the planimetric measurement of the reference volume sampled. Based on subsampling results, three sections per region per individual were quantified to reliably provide coefficients of error of less than 0.10. The section sampling interval ranged from 1 in 5 for smaller brains to 1 in 20 for larger brains. The mean number of sampling sites in the NAcc per individual was 29 ± 5 (s.d.) with an average of 131 ± 38 axon intersections counted per region. The mean number of sampling sites in the VP per individual was 28 ± 6 with an average of 72 ± 40 axon intersections counted per region. The Gundersen coefficient of error (CE, m = 1) was 0.06 ± 0.03.
Neuron density (Nv) and glial density (Gv) were assessed in adjacent Nissl-stained sections using the optical fractionator probe. While neuron densities were quantified for the purpose of standardizing the comparative data, glia densities were also collected at the same time. The ratio of glia to neurons (G/N) has been used to examine cell type composition across species and brain region (Sherwood et al. 2006). The NAcc and VP were outlined at 4x magnification and neurons and glia were manually counted under a 40x objective (N.A. 0.75) with a counting frame of 50 x 50 µm. As there is a significantly large range in the size of the specimens in this comparative sample, sampling grids also varied in size with an average of 550 x 459 µm (range: 78 x 111 – 1280 x 1702). The optical disector height was 7 µm with a 2% guard zone. Section thickness was measured at every 5th sampling site and the average mounted section thickness for Nissl-stained sections was 13.4 ± 4.9 µm. The Gundersen CE (m = 1) was 0.07 ± 0.003 for Nv and 0.06 ± 0.01 for Gv. Neurons were identified by the presence of a large, lightly stained nucleus and a distinct nucleolus, accompanied by lightly stained dendritic processes. Glia cells do not possess a visible nucleolus or dendritic processes. Due to the difficulty in differentiating astrocytes, oligodendrocytes, and microglia in Nissl-stained sections, all glia were included. Nv and Gv were calculated as the sum of neurons or glia counted, respectively, within the sum of optical disectors divided by the product and volume of the disector (e.g., Sherwood et al. 2006).
Because the specimens included in this study are derived from species with vastly different brain sizes, we used the ratio of ALv/Nv for comparative analysis among species. Using neuron density as a denominator provides a variable that accounts for the fact that axons are innervating neurons regardless of brain size. Additionally, glia densities (Gv) and glia-to-neuron ratios were evaluated for each region across species.
Statistical Analysis
We used IBM SPSS software,TIBCO Statistica Academicand R to analyze the data. Among-species differences were evaluated using analysis of variance (ANOVA) in the NAcc and VP separately. The dependent variables were TH ALv/Nv, Nv, Gv, and glia-to-neuron ratio (G/N). A Brown-Forsythe correction was applied when Levene’s test for homogeneity of variance was significant. Significant results were evaluated using Fisher’s Least Significant Difference (LSD) post hoc tests. Prior to among-species analyses, separate independent sample T-tests were used to test for differences between sexes in each species. The results from these analyses were non-significant (p > 0.05 for all), therefore sexes were pooled.
Additional analyses to determine allometric scaling relationships with brain weight were performed using R studio software (R core Team, Vienna, Austria, version 4.0.2). Species mean data were evaluated using phylogenetic generalized least squares (PGLS) to calculate scaling slopes while accounting for the covariance structure of evolutionary relatedness. In these analyses, Aotus spp. individuals were pooled with Aotus trivirgatus. Phylogenetic analysis of covariance (pANCOVA) was used to test whether human values were significantly different from what would be expected for their brain size (Smaers and Rohlf 2016). PGLS and pANCOVA were performed for NAcc and VP data separately. The dependent variables were TH ALv/Nv, Nv, Gv, and glia-to-neuron (G/N) ratio while the dependent variable was brain weight. The level of significance (α) was set at 0.05 for all statistical tests.