Abstract
Retinal ganglion cell numbers in adult mice vary from 40,000 to 80,000. Much of this variation and the prominent bimodality of strain averages are generated by allelic variants at the neuron number control 1 (Nnc1) locus on chromosome 11. TheNnc1 locus may modulate either ganglion cell production or the severity of ganglion cell death. Here we have determined what the relative contributions of these two processes are to variation in adult cell number by estimating total ganglion cell production in 10 strains of mice (A/J, BALB/cJ, BXD32, C57BL/6J, CAST/Ei, CARL/ChGo, CE/J, C3H/HeSnJ, DBA/2J, and LP/J). These strains have adult populations that range from 45,000 to 76,000 (data available athttp://qtl.ml.org). We estimated cell production by counting ganglion cell axons after ganglion cell neurogenesis but before the onset of significant cell death. Total cell production ranges from 131,000 to 224,000, and most of the variation in adult ganglion cell number is explained by this significant variation in cell production. In contrast, the percentage of cell death is relatively uniform in most strains (∼69% cell loss). The exceptions are BXD32, a strain that has an extremely high adult cell population, and Mus caroli (CARL/ChGo), a wild southeast Asian species that is distantly related to laboratory strains. In BXD32 and M. caroli, ∼62% of the population dies. Our analysis indicates that substitutions of single alleles at the Nnc1 locus are responsible for production differences of ∼8000 ganglion cells.
Numbers of retinal ganglion cells range from 50,000 in nocturnal rodents to several million in diurnal birds and primates (Rager and Rager, 1978; Rakic and Riley, 1983;Finlay and Pallas, 1989; Williams et al., 1996). Variation is also marked within species; numbers range from 0.7 to 1.5 million in humans (Curcio and Allen, 1990) and from 40,000 to 80,000 in mice (Williams et al., 1996). In mice, the distribution of strain averages is bimodal, with distinct modes centered at 55,000 and 64,000. We have demonstrated that this variation is primarily genetic, with a heritability of ∼0.8 (Williams et al., 1996). We subsequently discovered that most of the bimodality in strain averages is controlled by a major quantitative trait locus, neuron number control 1 (Nnc1). This quantitative trait locus is located on chromosome (Chr) 11 and is closely linked to three genes known to influence retinal development—the retinoic acid receptor α, neuregulin, and the thyroid hormone receptor α (Williams et al., 1998a). Nnc1may influence ganglion cell number by modulating either cell production or the severity of natural cell death.
Our principal objectives in this study were to determine the developmental mechanism by which allelic variants at Nnc1control neuron number and to assess the relative contributions of cell production and cell death to the consistent and pronounced differences in population size among strains of mice. We examined ganglion cell production in 10 strains—four selected from the low mode, four from the high mode, and one each from the high and low extreme (see Fig. 1). We specifically included strains C57BL/6J and DBA/2J, the parental strains of the set of recombinant inbred mice used to mapNnc1. Mechanisms generating neuron number differences between these two strains can be more confidently assigned toNnc1. Total ganglion cell production was estimated by counting retinal ganglion cell axons at birth, a time at which ganglion cell generation in mice is complete (Dräger, 1985) but before significant ganglion cell death has begun (Linden and Pinto, 1985).
MATERIALS AND METHODS
Animals. As illustrated in Figure1, strains of mice were chosen primarily to represent the two major modes in ganglion cell number (Williams et al., 1996). Three standard inbred strains were selected from the low mode (C57BL/6J, A/J, and LP/J) and four standard inbred strains were selected from the high mode (BALB/cJ, C3H/HeSnJ, CE/J, and DBA/2J). All of these strains were obtained from the Jackson Laboratory (Bar Harbor, ME). In addition, we selected two strains, CAST/Ei and BXD32, that have exceptionally low and high ganglion cell numbers, respectively. BXD32 was obtained from Dr. Benjamin Taylor at the Jackson Laboratory. CAST/Ei is an inbred strain derived from Mus musculus castaneus that we obtained from Dr. Eva Eicher at the Jackson Laboratory. Finally, we studied an outbred sample of Mus caroli that we refer to as CARL/ChGo, a strain that falls into the low ganglion cell mode. CARL/ChGo is a partially inbred strain ofM. caroli given to us by Dr. Dan Goldowitz at the University of Tennessee. Both CARL/ChGo and CAST/Ei are representatives of wild species endemic to southeast Asia. All mice were mated in our colony to produce neonates. The day of birth was designated postnatal day 0 (P0).
Tissue preparation. We anesthetized neonates by placing them on ice for several minutes. Neonates were then perfused transcardially with 0.1 m PBS (0.9%), followed by fixative (2.5% glutaraldehyde and 2.0% paraformaldehyde in 0.1 mphosphate buffer). Midorbital segments of optic nerves were dissected, osmicated, and embedded in Spurr’s resin. Nerves were thin-sectioned, placed on formvar-coated grids, and stained with lead citrate and uranyl acetate.
Estimating ganglion cell number. We estimated ganglion cell numbers by counting axons in optic nerve cross-sections (Williams et al., 1996). Previous studies have demonstrated that axon counts are reliable estimates of ganglion cell number; bifurcating axons, retinopetal axons, and retinoretinal axons are comparatively rare in mammals even during development [Perry et al. (1983); Chalupa et al. (1984); Lia et al. (1986); Williams et al. (1986); Rice et al. (1995), see their Table 2]. Nerves were photographed in a grid pattern at ∼30,000× using a JEOL EX2000II electron microscope. High and low magnifications were calibrated for each case by photographing a grid replica (2160 lines/mm; EMS, Fort Washington, PA). Unmyelinated axons were easily identified (Fig.2). Axons were counted directly on negatives within a 63 mm × 86 mm counting frame. Total axon estimates were calculated by multiplying the mean axon density by the total area of the optic nerve.
We counted necrotic axons in neonatal optic nerves from two strains belonging to the high mode and two strains belonging to the low mode. We did this by systematically scanning the entire optic nerve cross-section for necrotic axons at 15,000×. The criteria for distinguishing necrotic axons are those described by Williams et al. (1986). We also searched for growth cones in the sample of photographs used for counting axons and by scanning several optic nerves at high magnification (>40,000×).
Random numbers for Monte Carlo simulations were generated in Microsoft Excel 98. These numbers were drawn from a normal distribution with seed parameters (mean ± SD) taken from our experimental ganglion cell distributions (see Fig. 4 legend for more details).
RESULTS
The retinal ganglion cell population at birth ranged from 131,000 to 224,000 (Table 1). The mean for all 46 cases is 182,500 ± 4400 (± SE). This value is almost three times higher than the average for an equally diverse sample of adult mice (Williams et al., 1996). On average we counted five neonates per strain. The coefficient of variation within strains averaged 8.2%, only slightly higher than the 7.2% value obtained for adult mice (Williams et al., 1996). The small increase is probably caused by the technical difficulty of counting axons before they are myelinated. Given the anticipated variation in the stage of maturation of sets of neonatal mice killed at birth, this coefficient of variation is low and suggests that the ganglion cell population within a strain is comparatively stable at this stage of development. The average coefficient of error (the SE divided by the sample mean) averaged 4.5% in neonates and 2.5% in adults. These values provide an assessment of the reliability of adult and neonatal ganglion cell counts.
Cell production
If strain differences in adult ganglion cell numbers result from differences in the number of neurons that are generated, then at birth each strain should have a population that is approximately threefold higher than its adult mean. The slope of the regression should be close to 1/3, and the correlation should be high. This is what we found. The slope of a free regression for the 10 strains is 0.26 ± 0.07 (Fig. 3). Forcing the regression line through the origin produces the expected slope of 1/3 with an excellent fit (Fig. 3, inset). The positivey-intercept (11,600 adult cells) in the free regression may result from sampling error or nonlinearity of cell death or may indicate a basal level of cell production. The correlation coefficient of the free regression in Figure 3 is 0.81, and the corresponding coefficient of determination (r2) is 0.66. Thus, two-thirds of the variance in adult cell number can be readily explained by strain differences in cell genesis.
We were particularly interested in understanding the process that produces the bimodality of adult strain averages, and for this reason we also restricted the analysis to the eight strains belonging to high and low modes (Fig. 3, dark line). The coefficient of determination for this subset of points is 0.77, indicating that the bimodality is generated primarily by differences in ganglion cell production. The remaining “unexplained” variance must result from strain differences in the severity of cell death, developmental noise, and technical error.
Our statistical analysis is complicated by two factors. First, the parameters plotted in Figure 3 are not formally independent because total cell production cannot be less than the adult population. Second, the distribution of adult values is far from normal (Fig. 1). Conventional statistical estimates are therefore difficult to interpret. To address these problems, we performed Monte Carlo simulations to test cell production and cell death models using seed parameters taken from the adult distribution. We also subtracted the adult population from the neonatal population to insure independence between the parameters (Fig.4A). Figure 4,B and C, shows the outcomes of two typical Monte Carlo simulations in which we plot adult cell number against the number of lost cells. The first model (Fig. 4B) assumes that all differences in adult cell number are caused by matched differences in cell production and that cell death is strictly proportional to cell production. The second model (Fig. 4C) assumes that all differences among adult strains are caused by variation in the severity of cell death and that at birth all strains have approximately the same cell population (∼180,400 ± 18,400 cells). In the cell production simulation (Fig. 4B), the regression slope is +1.2, whereas in the cell death simulation (Fig. 4C), the slope is −1.1. Our actual data set (Fig. 4A) with its slope of +1.5 strongly supports a cell production model.
The 10 inbred strains were divided into high (BALB/cJ, C3H/HeSnJ, CE/J, BXD32, and DBA/2J) and low (C57BL/6J, A/J, CAST/Ei, CARL/ChGo, and LP/J) groups. Mean adult ganglion cell numbers for these groups are 66,800 ± 2700 and 50,900 ± 1900, respectively. There are highly significant differences in ganglion cell production between these groups, with means of 202,700 ± 7800 and 158,100 ± 10,600, respectively (t test, p < 0.001). In contrast, there is no significant difference in the percentage of ganglion cell loss between high and low groups, with mean percentages of cell loss relative to neonatal values of 66.9 and 67.5%, respectively (p = 0.42).
Nnc1 was mapped using recombinant inbred strains generated from the parental strains DBA/2J and C57BL/6J. For this reason a comparison between these two strains is especially germane in discovering how Nnc1 modulates ganglion cell number. The severity of cell death was closely matched between DBA/2J and C57BL/6J (69 and 70%, respectively). In contrast, DBA/2J produces ∼16,400 more cells than does C57BL/6J. This result, together with our previous finding of additive gene action (Williams et al., 1998a), indicates that the substitution of a single allele at Nnc1 is associated with a production difference of ∼8000 cells.
Cell death
With the exception of strains BXD32, CARL/ChGo, and BALB/cJ, the average percentage of cell death among strains is relatively uniform (Table 1, 69 ± 1.2%). Although the percentage of cell death is relatively uniform, the absolute magnitude of ganglion cell death is variable among strains and is highly correlated with production values (Table 1). There are some notable exceptions to this generality. The percentage of cell death in BXD32 and CARL/ChGo is significantly lower than that in other strains (t test, p < 0.05, with Bonferroni correction). Estimates of ganglion cell production are similar in BALB/cJ and C57BL/6J, yet these strains have adult populations that differ by ∼8000 cells (Table 1). A slight reduction in the severity of cell death in BALB/cJ (65% loss) seems to account for the relatively high cell number of this strain at maturity. In this instance, the marked strain difference in adult population size results primarily from variation in the severity of cell death. Differences in cell death can also compensate for differences in cell production. For example, CARL/ChGo produces an average of 131,000 ganglion cells, 20,000–40,000 fewer cells than A/J and LP/J produce, respectively, yet all three strains have closely matched adult populations (Table 1).
Necrotic axons and growth cones
The validity of our quantitative analysis depends on the assurance with which we can estimate total ganglion cell production in mice. If much cell loss occurs before birth or much cell addition occurs after birth, then production estimates based on axon counts in the optic nerve at P0 will be too low. To eliminate the possibility that significant cell death occurs prenatally, we counted necrotic axons in neonatal optic nerves from strains belonging to the high and low modes using criteria described by Williams et al. (1986). Necrotic axons are relatively easy to see, and it was practical to count all sites of necrosis in single optic nerve cross-sections. Necrotic axons at P0 make up 0.02 and 0.05% of the fiber population in cases selected from the low strains, A/J and C57BL/6J, respectively, whereas they make up 0.07 and 0.09% of the population in cases selected from the high strains, BXD32 and C3H/HeSnJ, respectively. The fact that a somewhat higher incidence of necrosis was noted in nerves taken from the high strains makes it very unlikely that variation in early axon loss accounts for differences between adult values. Growth cones were exceeding rare in all material, and fewer than five profiles among all cases met even a relatively lax criteria for these structures (Williams et al., 1986; Williams et al., 1991; Colello and Guillery, 1992).
Specificity of strain differences
Do strain differences in retinal ganglion cell number correspond to differences in total brain weight, or are differences among strains specific to the ganglion cell population? The correlation of ganglion cell number and brain weight across individual mice is 0.37, but when strain averages are used, the correlation rises to 0.75. This suggests that approximately one-half of the variance in neonatal ganglion cell number can be explained directly or indirectly by differences in brain weight. As assessed by quantitative DNA analysis, brain weight differences among neonatal mice are primarily attributable to differences in total cell number (Zamenhof and Marthens, 1976). Thus, mechanisms modulating ganglion cell number may have common effects on cell number in the other parts of the CNS. The correlation between strain averages of adult brain weight and ganglion cell number for the same strains is only 0.51. Given the wide confidence intervals of correlations computed with low numbers of cases, the difference between the adult and neonatal correlations (0.51 and 0.75, respectively) may be caused by sampling error. But it is also conceivable that strain variation in cell death decreases an initially high correlation between brain weight and retinal ganglion cell number. In any case, the cellular specificity of the strain differences is likely to be low, and we expect differences in numerous other neuronal cell populations to be closely matched with the differences we find in ganglion cell number.
DISCUSSION
Synopsis
Our analysis demonstrates that most of the variation in adult ganglion cell number among strains of mice can be traced to differences in cell production. Allelic variants at the Nnc1 locus on Chr 11 (Williams et al., 1996, 1998a) generate the pronounced bimodality in ganglion cell population size by modulating ganglion cell production.
Generation and death of retinal ganglion cells
Generation of retinal ganglion cells in mice begins on embryonic day 11 (E11) and lasts until just before birth (Dräger, 1985). There is a short delay between neurogenesis and the time at which ganglion cell axons extend into the optic nerve (Colello and Guillery, 1992). This delay could deflate estimates of total cell production. However, very few ganglion cells are produced after E18 (Dräger, 1985), and as anticipated from the work of Colello and Guillery (1992), we did not observe growth cones in neonatal optic nerves. It is therefore unlikely that our estimates of total production are biased downward by late ganglion cell generation.
In contrast, ganglion cell death begins at, or just before, birth, peaks between postnatal days 4–6, and is essentially complete by P12 (Linden and Pinto, 1985). At the peak of cell loss, between 5,000 and 10,000 ganglion cells are eliminated per day (Williams et al., 1990). However, fewer than 2000 cells are lost on the day of birth in mice, consistent with our observation of very few necrotic axons, <300 per nerve. At this rate it is improbable that more than a total of 10,000 ganglion cells are lost prenatally. In chickens there is as much as a 3 d delay between the onset of ganglion cell degeneration in the retina and the elimination of axons in the optic nerve (Rager and Rager, 1978). If there is a similar delay in mice, then our axon counts should more directly reflect production numbers. Nonetheless, our estimates of total production may be biased downward slightly by the early loss of ganglion cell axons. However, the magnitude of this error is sufficiently small (∼10,000 cells) that we did not think this loss warranted correction.
Mechanism generating differences in ganglion cell production
We recently mapped a gene, Nnc1, that is responsible for more than one-half of the genetic variance in ganglion cell number in mice and that generates the pronounced bimodality that we discovered among strain averages (Williams et al., 1998a). Nnc1 is the first locus known to control normal variation in cell number in the vertebrate CNS. The thyroid hormone receptor α gene (Thra) is a superb candidate gene. Thra maps within 1–2 cM of Nnc1 on chromosome 11 (Montgomery et al., 1997) and is expressed within the developing chick retina (Sjöberg et al., 1992). The ligand of THRA, triiodothyronine, is known to influence retinal ganglion cell fate determination (Hoskins, 1985) and retinal maturation rate (Macaione et al., 1984), and hypothyroidism during retinal development results in decreased cell density in the ganglion cell layer (Hoskins, 1985; Navagantes et al., 1996).
Nnc1, possibly Thra, could influence ganglion cell production by affecting (1) the numbers of retinal progenitor cells, (2) the pathways of cell determination, or (3) the kinetics of progenitor cell proliferation. Genetically determined differences in numbers of multipotent retinal progenitors would have consistent effects on the number of many retinal cell types. However, a comparison of horizontal cell and ganglion cell numbers for six strains demonstrates that ratios in these early-generated cell types are not always matched (Williams et al., 1998b). We have also examined other retinal cell populations and in preliminary work have found a weak negative correlation between ganglion cells and photoreceptors (Williams et al., 1998a). This suggests that there may be a reciprocal relationship between the generation of early- and late-generated retinal cell types. An example of this type of reciprocal relationship is found when the Notch signaling pathway is perturbed. This results in a shift in the ratio of early- to late-generated retinal cell types (Dorsky et al., 1997). The idea of temporally regulated competence is supported by a confluence of work (Watanabe and Raff, 1990; Anchan et al., 1991; Williams and Goldowitz, 1992; Guillemot and Joyner, 1993;Cepko et al., 1996; Alexiades and Cepko, 1997).
The progressive slowing of the cell cycle and its eventual cessation result in part from decreased availability of key exogenous factors (Jacobson, 1991; Alexiades and Cepko, 1996). Studies have identified multiple mitogenic factors for retinal progenitors: FGF, TGFα (Lillien and Cepko, 1992), TGFβ, epidermal growth factor (Anchan et al., 1991), and IGF-1 (Hernández-Sánchez et al., 1995; Frade et al., 1996). Interestingly, the addition of the THRA ligand triiodothyronine to cultured fetal rat hypothalamus cells stimulates the release of IGF-1 into the culture medium (Binoux et al., 1985). Finally, genetic variants of Nnc1 could alter the proliferation kinetics of progenitors that give rise to ganglion cells. Rates of mitosis may be influenced via inhibitory molecules and pathways. One extremely interesting example is dopa, a tyrosine metabolite that normally has inhibitory effects on cell genesis in retina (Ilia and Jeffery, 1996). The absence of dopa in albino rats leads to an anomalous upregulation of ganglion cell production followed by an increase in the severity of cell death (Ilia and Jeffery, 1998). Nnc1 could have effects within any of these mitotic regulatory pathways.
Variation in retinal ganglion cell death
The severity of cell death is close to 68–70% in most strains of mice. However, there are three exceptional strains with less severe loss. Three to nine percent fewer cells are lost in BXD32, CARL/ChGo, and BALB/cJ. BXD32 is particularly interesting because it has the highest adult population (75,800 ± 1900) among the 60 strains we have now examined. Yet at birth BXD32 has an unexceptional number, 199,500, that is lower than that of other strains. Clearly, one or more genes controlling rates of ganglion cell death are responsible for the high adult cell number in this strain. It would be feasible to map a cell death gene by crossing BXD32 to a strain with similar ganglion cell production but higher cell death.
Variation in the severity of cell death may result from differences in titers of neurotrophic factors. The neurotrophins (BDNF and neurotrophin-3/4) have been found to increase survival of retinal ganglion cells in chicken and rat (Rosa et al., 1993; Ma et al., 1998). Neuregulin, found on the cell surface and as a secreted protein, can also increase survival of neonatal rat retinal ganglion cells in culture (Bermingham-McDonogh et al., 1996). Differences in the time of expression or concentration of these neurotrophic factors, their receptors, or components within their signaling pathways could produce variation in the severity of naturally occurring ganglion cell death.
Nnc1 controls cell production
We have shown that as much as 77% of the variation among adult strains results from differences in the production of ganglion cells. The percentage of cell death in high and low groups does not differ significantly (66.9 and 67.5%, respectively). We conclude that variation in adult ganglion cell number among inbred mouse strains results predominantly from differences in cell production. Comparison of our data with the Monte Carlo simulations (Fig. 4) corroborates this conclusion.
A major motivation for undertaking the present study was to determine how and when allelic variants at Nnc1 influence the size of the ganglion cell population. Collectively, our results strongly indicate that Nnc1 modulates ganglion cell number by influencing cell production, and because ganglion cell production occurs before birth, our results indicate a time frame for the action of Nnc1.
Footnotes
This research was supported by grants from the National Institutes of Health to R.W. (National Institute of Neurological Disorders and Stroke Grant R01 NS35485 and NEI Grant EY08868). R.C.S. was supported in part by the United States Public Health Service Training Grant RNS-07323. We thank K. Troughton for technical help and D. Goldowitz for providing us with M. caroli pups.
Correspondence should be addressed to Dr. Robert W. Williams, Department of Anatomy and Neurobiology, 855 Monroe Avenue, Memphis, TN 38163.