Select      
 Site search   
  Home    Publications

Browse Publications
 
List of Contents

Quantitative Neurogenetics & QTL Mapping

Genetics of Myopia

Control of Neuron Number and Stereology

Growth Cones and Dying Axons

Retina Development and Visual System Mutants

Grant Application

U.S. Patent

Abstracts


Need Help?
Help with Publications
Help with Nervenet
Contact Us

     
Note to the Reader This is a preprint of a paper now in press at the Journal of Neuroscience (Spring 2001). Please cite this work as: Lu L, Airey DC, Williams RW (2001) Complex trait analysis of the hippocampus: Mapping and biometric analysis of two novel gene loci with specific effects on hippocampal structure in mice. J Neurosci 21:3503-3514. The definitive print and html versions are available on the Journal of Neuroscience web site.

Print Friendly
Complex trait analysis of the hippocampus:
Mapping and biometric analysis of two novel gene loci with specific effects on hippocampal structure in mice


Lu Lu, David C. Airey, Robert W. Williams
Center for Neuroscience and Department of Anatomy and Neurobiology, University of Tennessee, 855 Monroe Avenue, Memphis, Tennessee 38163 USA
Email questions and comments to lulu@nb.utmem.edu

 

Contents

Introduction
Methods
     Volumetric Analysis
     Stereology
     Regression Analysis
     Genotyping

Results
     Table 1: Hippocampal Weight Data
     Table 2: Correlation Matrix
     Figure 2: Regression Plots
     Figure 3: Relations among Subregions
     Figure 4: Interval Maps


Primary Data Files
     BXD strain data in Map Manager text format: New data set
     BXD data for individual mice: Sex, age, weight, etc.
     F2 genotype and phenotype data in Map Manager text format
     F2 data for individual mice: Sex, age, weight, etc.

Discussion
     Structure and Function
     Behavioral QTLs and Hipp1
     Specificity of QTL Action
     Candidate Genes
     Sex Difference
     Adult Neurogenesis and Hippocampal QTLs

References

 

 

ABSTRACT

Notable differences in hippocampal structure are associated with intriguing differences in development and behavioral capabilities. In this study we explore genetic and environmental factors that modulate hippocampal size, structure, and cell number using sets of C57BL/6J (B6) and DBA/2J (D2) mice; their F1 and F2 intercrosses (n = 180); and 35 lines of BXD recombinant inbred (RI) strains. Hippocampal weights of the parental strains differ by approximately 20%. Estimates of granule cell number also differs by 20%. Weights of RI strains range from 21 to 31 mg; those of F2 mice from 23 to 36 mg (bilateral weights). Volumes and granule cell numbers are well correlated (r = 0.7–0.8). An important component of the structural variation is associated with age and sex differences. The hippocampus gains weight with age (~0.24 mg/month from 2 to 10 months), and hippocampi of males typically weigh 0.55 mg more than those of females.
     Heritability of variation is ~50% and half of this genetic variation is generated by two quantitative trait loci (QTLs) that we have mapped to chromosome 1 (Hipp1: genome-wide p < 0.005, between 65 and 100 cM) and to chromosome 5 (Hipp2, p < 0.05, between 15 and 40 cM). These are the first gene loci known to produce normal variation in forebrain structure in mammals. Hipp1 and Hipp2 individually modulate hippocampal weight by 1.0 to 2.0 mg, an effect size greater than that generated by age or sex. The Hipp gene loci modulate neuron number in the dentate gyrus, collectively shifting the population up or down by as much as 200,000 cells. Candidate genes for these loci include Rxrg and Fgfr3.

Abstract in Chinese
Our thanks to Dr. Baoyan Hu for this translation.
 

 

Introduction

Structural variation in hippocampus can be substantial. Among normal humans hippocampal volume varies roughly two-fold—from 6 to 11 cm3 (Csernansky et al., 1998; De Bellis et al., 2000). Classic biometric studies by Richard and Cynthia Wimer and their colleagues demonstrated that pervasive differences in hippocampal structure among inbred strains of mice are due to a complex mixture of genetic, environmental, and maternal effects (Wimer and Wimer, 1971; Barber et al., 1974; Wimer et al., 1976; Wimer et al., 1978; Wimer and Wimer, 1982). This well characterized variation in mice is now especially intriguing for two reasons. First, the discovery, characterization, and manipulation of neuronal and glial stem cells in rodent forebrain (Altman and Das, 1965; Bayer et al.,1982; Kaplan and Bell, 1984; Stanfield and Trice, 1988; Cameron et al., 1993; Gould, 1994) has invigorated research directed at isolating factors that modulate rates of stem cell cycling and proliferation (Palmer et al., 1997; Kuhn et al., 1997; Kempermann et al., 1997ab, 1998). Gage and colleagues have shown marked differences in the dynamics of cell proliferation in the dentate gyrus among inbred strains, including the two key strains, C57BL/6 and DBA/2, that we have used in the present study (Kempermann et al., 1998; Kempermann, personal communication).
     The second complementary reason to be interested in strain differences is that recent advances in quantitative and molecular genetics make it practical to map, characterize, and clone gene loci that influence a wide range of heritable neuroanatomical and behavioral traits (Williams, 1998; Rikke and Johnson, 1998; Williams, 2000). This forward genetic approach initially involves associating differences in alleles at defined chromosomal positions with differences in a trait—in our case the size, architecture, and cellular composition of the mouse hippocampus. A strong association between phenotypes and genotypes indicates the presence and position of a quantitative trait locus, or QTL (Darvasi, 1998; Williams 2000). Complex trait analysis has proved to be highly effective for neuroanatomical and behavioral traits. For example, by taking advantage of differences between several strains of mice, we recently succeeded in mapping a set of gene loci that specifically modulate numbers of neurons in the mouse retina (Williams et al., 1998; Strom and Williams, 1998; Strom 1999). The deluge of sequence data, combined with greatly improved QTL mapping strategies, promises to revolutionize the quantitative genetic analysis of CNS development (Darvasi, 1998; for progress on mouse sequencing efforts see, www.ncbi.nlm.nih.gov/genome/seq/MmProgress.shtml). In the case of the hippocampus, it is now possible to use QTL methods to highlight a crucial subset of genes and alleles that modulate hippocampal structure and function at all stages of life.
     In this study we use biometric and genetic techniques to explore the basis of variation in the size and structure of the hippocampus. Our work begins with a systematic multiple linear regression analysis of the effects of age, sex, body and brain weight on the adult hippocampus. This work is of interest in its own right, but in the context of gene mapping studies is a prelude to interval mapping. We have succeeded in identifying and verifying two gene loci that control the size of the hippocampus. These two loci are among the first genes associated with normal variation in forebrain anatomy in any vertebrate. We have examined the cytoarchitecture of the hippocampus and granule cell density in the dentate gyrus to assess the scope of action of these novel gene loci.

 

Material and Methods

Two complementary groups of mice were used in this study. The first consisted of a set of 35 BXD/Ty recombinant inbred (RI) strains. All of these strains were used both for gene mapping and for the biometric analysis of the hippocampus and its parts. The BXD RI strains were generated by crossing C57BL/6J (B6) and DBA/2J (D2) parental strains in the mid-1970s (BXD1 through 32) and 1990s (BXD33 through 42) by Benjamin Taylor and colleagues (1989, 1999). RI strains are completely inbred lines derived from brother-sister matings starting from an F2 intercross. They are advantageous for complex trait analysis for several reasons explained in Williams (1998, 2000). First, numerous individuals with identical genomes can be phenotyped and this greatly improves the precision with which phenotypes such as hippocampal weight can be measured. The key parameter used during QTL mapping becomes a strain mean rather than an individual value. We typically averaged values from eight animals for weight data and five animals for volumetric data (both sides were analyzed in both cases), with the goal of reducing standard error of the mean to well under 1 mg. In fact, the coefficient of error (SE/mean) averages only 1.8% whereas the coefficient of variation (SD/mean) averages 3.7% (Table 1). A related advantage of RI strains is that they can be shared among many investigators over a period of many years. The large strain differences in hippocampal structure and cell density that we report can now be followed up by developmental, physiological, pharmacological, and behavioral studies. From the point of view of gene mapping, one drawback of RI strains is that only a limited number of strains are available for analysis, a factor that decreases the power to detect QTLs with small effects. For this reason we have supported our RI work with a matched analysis of an F2 intercross.
     The second group of mice that we used consisted of reciprocal F1 and F2 intercrosses between B6 and D2, the same strains used to make the BXD RI strains. Thus these two groups share precisely the same sets of parents and parental alleles. The F2 animals were generated at the University of Tennessee by intercrossing both BDF1 and DBF1 mice as described in Zhou and Williams (1999); 105 of the animals were BDF2s and 75 of the animals were DBF2s. An F2 intercross is often used to complement the analysis of recombinant inbred strains, and in this study we have used the F2 to confirm and refine the location of QTLs that modulate hippocampal weight. An F2 intercross has the advantage over the RI set of a much larger number of genomes to correlate with phenotypes. The disadvantage to this cross is the effort required to genotype each individual mouse. In our study 180 animals were genotyped at 128 microsatellite marker loci as describe below.

Animal husbandry and age

Mice were maintained at 20–24 °C on a 14/10h light-dark cycle in a pathogen-free colony at the University of Tennessee. Most animals were fed a 5% fat Agway Prolab 3000 rat and mouse chow and given tap water in glass bottles. The average age of BXD/Ty animals was 80 days; that of F2 mice was 98 days (F2 data in text format for cases). Parental strains and the set of BXD/Ty strains were obtained from the Jackson Laboratory (Bar Harbor, Maine) from 1994 through 1999.

Fixation

Mice were deeply anesthetized with Avertin (1.25% 2,2,2-tribromoethanol and 0.8% tert-pentyl alcohol in water, 0.5–1.0 ml ip). Most mice were perfused through the heart with 0.1 M phosphate buffered saline followed by 5% glutaraldehyde and 1.0% paraformaldehyde in 0.1 M phosphate buffer, and then by 2.5% glutaraldehyde and 2.0% paraformaldehyde in 0.1 M buffer.

Dissection

Fixed brains were bisected along the midline (Fig. 1). Left and right hippocampal regions were dissected under a dissecting microscope by inserting fine blunt forceps into the ventricular cavity just dorsal to the hippocampus and removing overlying cortex and callosum. The surface of the hippocampus and dentate gyrus was used to guide removal of cortex along the septotemporal axis. The exposed hippocampus and dentate gyrus was pulled free of the hemisphere in a ventral-to-dorsal direction. The dissection includes a small part of the subiculum adjacent to CA1 and occasionally a small strand of the fimbria. As shown in the Results, technical error and inconsistency was very low as judged by right left differences witiin cases. The most dorsal-anterior aspect of each hippocampus was then trimmed free of septum and dorsal fornix, rolled quickly in tissue paper, and immediately weighed to the nearest 0.1 mg. The first author made all dissections. The following data were obtained for virtually all cases: sex, age at death, body weight, brain weight, cerebellar weight, and left and right hippocampal weights. Original data files are available at www.nervenet.org/databases/hippocampus.html.

Figure 1: Dissection of the Mouse Hippocampus
 

Figure 1. An example of the dissected hippocampi used to generate the weight data for mapping Hipp1 and Hipp2. This dissection is from the left hemisphere and is oriented properly with respect to the septotemporal (S and T) axis of the inset mouse brain. The internal anatomy of the hippocampus is referenced by five insets showing coronal Nissl stained sections (a-e). Scale bar is 1 mm. [Click on figure for a larger version.]

 

Volumetric measurement of the hippocampus

A second independent set of 31 BXD strains and the two parental strains that are part of the Mouse Brain Library (www.mbl.org) was studied to assess specificity of QTL action, and to determine if QTLs affecting hippocampal size have regional effects limited to the dentate gyrus, the hippocampus proper, the pyramidal cell layer, or the granular cell layer. Images of the serial sections through the entire hippocampus from 154 BXD cases (average of 5 per strain), 8 C57BL/6J cases, and 6 DBA/2J cases were downloaded from www.mbl.org and analyzed in NIH Image. Images were not available for four BXD strains 6, 16, 21, 37. Serial section images have a resolution of 4.5 µrs of the hippocampus, excluding the subiculum, but including the fimbria and dentate gyrus (see Figure 1, a-e) were traced manually. The dentate gyrus was traced separately. The area of the set of sections was multiplied by the section interval to estimate volume. Finally, areas of the pyramidal cell layer (CA1–CA3) and the granule cell layer were measured bilaterally for all cases. These layers have a high cell packing density and can be reliably defined by a thresholding operation in NIH Image. To improve uniformity of the data set the first author conducted all thresholding. The reliability of thresholding was tested by repeated measures analysis. The correlation of duplicated estimates is 0.95. Shrinkage among cases in the Mouse Brain Library is variable, but known. To correct for variance due to shrinkage, we divided hippocampal volumes by the total brain volume and then multiplied by the brain volume expected from the known brain weight, assuming a brain density of one mg per mm3 of fixed tissue. The post-processing volume of the total brain was measured by point counting as described in Williams (2000). All volume data have been corrected for shrinkage, case-by-case. The correlation between weight and volume from two different sets of BXD animals is 0.61 (df = 31, p < 0.01). Volumetric data used in the results are group means based on 4 to 6 cases per stain.

 

Stereological analysis of cell density in the dentate gyrus of BXD strains

Mean cell density in the dentate gyrus was determined bilaterally in an average of five cases of each of 33 BXD strains. For this work we used tissue from the Mouse Brain Library described in the preceding paragraph. All tissue was embedded in celloidin and cut at 30 µm in coronal or horizontal planes by G.D. Rosen and colleagues. Precise volumetric shrinkage for each case had been previously computed as described at www.mbl.org and this has made it possible to compute the expected density for fixed brain tissue prior to processing (Table 1; see web movie at www.nervenet.org/mbl/movies/Z-hilus.html). Granule cells, glial cells, and endothelial cells were counted separately using the 3D counting procedure of Williams and Rakic (1988). Glial cells and endothelial cells make up only a very small fraction of the total cell count (<10%) and were relatively easy to distinguish. The volume of the count box was fixed at 17 x 18 x 30 µm and did not include an upper guard space. The z-axis was monitored using a calibrated rotary shaft encoder as described in Williams and Rakic (1988). This modification of the 3D-counting protocol avoids potential bias introduced by differential z-axis shrinkage (Hatton and von Bartheld, 1999; von Bartheld 1999). The sample volume was counted using a contrast-enhanced video differential interference contrast system and a 1.25 NA 100x planachromat objective. Final magnification on the monitor was 4000x. All cell density estimates are corrected case-by-case for volumetric shrinkage.
     We computed the mean neuron cell packing density in the granule cell layer at a consistent location 2.0 to 2.3 mm posterior to bregma, at a level that corresponds to the central part of the dorsal lateral geniculate nucleus (Fig. 1B). The counting boxes were located in the dorsomedial apex. For each RI strain the mean cell density estimate is based on ten counts from five individuals. These density estimates have a very low coefficients of error (3% or 60,000 cells/mm3) across cases within strains. Given this single point sampling protocol, it was of interest to assess how representative values obtained at minus 2.0 AP are of the entire dentate gyrus. For a subset of five cases (three C57BL/6J and two DBA/2J cases) we explored this question by carrying out high density sampling of sections spaced at 300 µm intervals in coronal and horizontal planes. An average of 70 fields were counted using a systematic random protocol (Howard and Reed, 1998). Cell densities tend to be higher in the rostral and dorsal regions than in the caudal temporal region. The difference between extreme poles is approximately 40%. The gradient of cell density across the granule cell layer is more modest. The marginal part the granule cell layer (both molecular and hilar sides) has a density that is approximately 10–20% lower than that of the most central part the granule cell layer. The sampling sites that we selected for analysis is intermediate in both position and in mean cell density and yield good estimates of total granule cell number. For example, comprehensive estimates of granule cell number from three C57BL/6J mice obtained using a dense sampling grid (60 to 80 samples per side) are 484,000 ± 16,000, 465,000 ± 16,000, and 435,000 ± 16,000. This gives a unilateral average of 443,000 granule cells. In comparison, the estimate given in Table 1 that is based on counts from 5 cases and 10 bilateral samples taken at AP –2 is 443,000 with a standard error of about 15,000 (note that the estimates in Table 1 are given for the sum of both sides and for comparison were divided by two). Our estimate for this strain is close to that recently obtained by Abusaad and colleagues (493,000 based on six 9-week-old cases also of both sexes). Similarly, estimates for the two comprehensive estimates of the other parent, DBA/2J, were 360,000 ± 11,000 and 312,000 ± 8,000. This compares to a mean of 349,000 cells in Table 1. The agreement is satisfactory and estimates listed in Table 1 are of sufficient accuracy to gauge effects of the two Hipp loci on dentate granule cell populations.

 

Regression analysis

Complex trait analysis tests the strength of the relationships between genotypes and phenotypes. The advantages and restrictions of this biometric approach have been reviewed extensively (Tanksley, 1993; Lander and Schork, 1994; Lynch and Walsh, 1998). One issue that is important in our analysis of the hippocampus is the specificity of gene effects. We did not want to inadvertently map genes that control body or brain weight. For this reason we have carried out linear regression analysis to explore and statistically control for covariance between hippocampal weight and other variables using standard techniques explained more fully elsewhere (Sokal and Rohlf, 1995; Williams, 2000). For example, strain BXD5 has an unusually large brain (544 mg), and not surprisingly has an unusually large hippocampus (30 mg). However, after normalizing for brain weight with regression, a process that involves computing residual weight measure, the hippocampus of this strain is not at all exceptional. Computing residuals is preferable to the use of ratios to correct for size differences (Lynch and Walsh, 1998; Bishop and Wahlsten,1999). Controlled parameters in our multiple linear regression analysis include brain weight, sex, age, and body weight. Interactions and second-order terms were not significant. The standard assumptions required from regression analysis were met. Statistics and statistical tests were computed using Data Desk (Data Description Inc, www.datadesk.com).

 

Genotyping and QTL mapping

In essence, the QTL mapping involves categorizing cases (genetic individuals) based on their genotypes (e.g., BB, BD, or DD) at defined chromosomal markers (e.g., microsatellite loci) and comparing these groups with a quantitative variable; in our case hippocampal measurements (Table 1). A gene locus that affects the hippocampus will be recognized when the variation in phenotype matches the variation in genotype (see www.nervenet.org/papers/brainrev99.html for several examples). QTL mapping generally proceeds from analysis at defined loci (single marker analysis), to an analysis at positions inferred between loci (simple interval mapping), and then finally to positions inferred between loci but with statistical control for background loci of interest (composite interval mapping). Map Manager implements both simple and composite interval mapping methods described by Haley and Knott (1992), and evaluates genotype-phenotype associations with likelihood statistics and permutation tests. Genome-wide significance levels for assessing the confidence of the linkage statistics is estimated by comparing the peak likelihood ratio statistic (LRS) of correctly ordered data sets with LRSs computed for 10,000 permutations (Churchill and Doerge, 1994). Permutation tests are a commonly accepted method for determining the probability of the effect occurring by chance. LRS scores can be converted to logarithm of odds ratios (LOD scores) by dividing by 4.6.
     Genomic DNA from all F2 mice was extracted from their spleens using a high-salt procedure (Laird et al., 1991; see PCR details ). A set of 128 microsatellite loci distributed across all autosomes and the X chromosome were typed in the F2 progeny using a modified version of the protocol of Love et al. (1990) and Dietrich and colleagues (1992). Each 10 m l PCR reaction contained 1
X PCR buffer, 1.92 mM MgCl2, 0.25 units of Taq DNA polymerase, 0.2 mM of each deoxynucleotide, 132 nM of the primers, and 50 ng of genomic DNA. The microsatellite primer pairs were purchased from Research Genetics (Huntsville, AL, www.resgen.com). A loading dye (60% sucrose, 1.0 mM cresol red) was added to the reaction before the PCR (Routman et al., 1994). PCRs were carried out in 96-well microtiter plates. We used a high-stringency touchdown protocol in which the annealing temperature was lowered progressively from 60 °C to 50 °C in 2 °C steps over the first 6 cycles (Don et al., 1991). After 30 cycles, PCR products were run on 2.5% Metaphor agarose gels (FMC Inc., Rockland ME), stained with ethidium bromide, and photographed. Genotypes were entered into Microsoft Excel 98 and transferred to Map Manager QT and QTX for mapping and permutation analysis. A nonredundant set of loci typed for all 35 BXD/Ty strains was used for analysis (Williams et al., 1998; Zhou and Williams, 1999). Approximately 200 genotypes were taken from Taylor et al., (1999), and an additional 630 new genotypes were generated for this analysis in our laboratory using the same methods described above for the F2. The new genotypes are available from the Informatics Center for Mouse Neurogenetics (www.nervenet.org/MMfiles/MMlist.html). Map Manager implements both simple and composite interval mapping methods described by Haley and Knott (1992). Genome-wide significance levels for assessing the confidence of the linkage statistics was estimated by comparing the peak likelihood ratio statistic (LRS) of correctly ordered data sets with LRSs computed for 10,000 permutations (Churchill and Doerge, 1994). LRS scores can be converted to LOD scores simply by dividing by 4.6.


 

Results

The results are divided into two major sections. The first is a biometric analysis of normal variation in the structure and size of the hippocampus as a function of sex, age, brain weight, and body weight. We initially focus on hippocampal weight data collected for a large number of cases, including all BXD strains (Table 1), 180 F2 progeny, two parental strains, and the reciprocal F1 hybrids. The analysis of weight data is followed by an analysis of absolute and relative volumes of four parts of the hippocampus and a stereological analysis of granule cells in the dentate gyrus of BXD strains. The second section summarizes QTL mapping results and the effects of two gene loci on hippocampal weight, structure, and granule cell number. The analysis of effect specificity of QTLs relies on the BXD strains because a number of individuals of a defined genotype could be examined to ensure accurate trait means. The multiple regression analysis ensures that the QTLs that we have mapped have comparatively selective (although not exclusive) effects on size, structure, and cell populations of the hippocampus.

 

Normal Variation in Size, Structure, and Cell Number

Hippocampal weights of parental strains

Hippocampal weights of the parental strains B6 and D2 are 28.4 ± 0.7 and 23.0 ± 0.3, respectively (Table 1). This 19% difference is highly significant (t(20) = 17.0, p < 0.001) and corresponds to a 20% difference in number of granule cells in the dentate gyrus of these two strains (Table 1). Total brain weight of B6 is also approximately 20% greater than that of D2 (496 ± 6 mg versus 415 ± 4 mg) whereas body weight is approximately 28% greater at 75 days (24.7 versus 19.3 gm). If hippocampal weights of these strains are adjusted to account for differences in brain and body weight, then the strain difference is reduced to 1 mg: 27.4 ± 0.3 versus 26.4 ± 0.3 (Table 1, column 3). This difference is still statistically significant (t(24) = 17.0, p = 0.04). The hippocampal weight of the reciprocal F1 hybrids—BDF1 and DBF1—are 27.8 ± 0.3 and 28.2 ± 0.4 mg, respectively, an insignificant difference from each other or from the B6 parent.

 

Table 1. Strain Variation in Hippocampal Structure

 

Brain weight and hippocampal weight

Variation in brain weight is the single most important predictor of variation in hippocampal weight (r = 0.73, Table 2, Fig 2A), and 50% of the variance in hippocampal size among individual BXD mice can be accounted for by the regression equation: hippocampal weight = 4.08 + 0.05(brain weight in mg). Brain weight includes the hippocampus and to compare fully independent variables we recomputed the relationship after subtracting the hippocampus. Although the slope is almost precisely the same, the amount of variance explained by the equation is reduced by 5.0% (Fig. 2B, Table 2). Among the F2 sample, variation in brain weight is also the most important predictor of variation in hippocampal weight: 57% of variance in hippocampal weight is accounted for by brain weight, whereas 51% is accounted for by brain minus hippocampus.

Figure 2A: Correlation of brain weight and hippocampal weight
 

Figure 2B: Correlation of brain weight and hippocampal weight
 

Figure 2A and 2B. Brain weight effects on hippocampal weight (bilateral). (A) 53% of the variation in hippocampal weight is explained by brain weight. (B) 47% of the variation in hippocampal weight is explained by brain minus hippocampal weight. There are no significant differences in brain weight among our sample of mice. In this and other parts of Figure 2 males and females are coded as blue crosses and red circles, respectively.

 

Body weight and hippocampal weight

Body weight and hippocampal weight are correlated among BXD mice (r = 0.42, Table 2, Fig 2C). A 1 g increase in body weight is associated with a 0.19 ± 0.03 mg increase in hippocampal weight (Fig. 2C). However, there is no significant correlation when brain weight is used as a cofactor in a multiple regression analysis (Fig. 3B). In F2 mice body weight correlates even more weakly with hippocampus (r = 0.23), and variation in body weight accounts for only 5% of variation in hippocampal weight. A 1 g increase in body weight is associated with a 0.10 ± 0.03 mg increase in hippocampal weight.

Figure 2C: Body Weight and Hippocampal Weight
 

Figure 2C: Body Weight and Hippocampal Weight
 

Figure 2C. Body weight effects on the hippocampal weight. Note the large sex difference in body weight. (C) Body weight is significantly correlated with hippocampal weight. (D) X and Y axes represent body and hippocampus residuals respectively, after correction for covariation with brain weight.


 

Age and hippocampal weight

The weight of the hippocampus increases markedly with age. Among BXD individuals that ranged in age from 30 to 300 days, the slope of this increase is 3.4 mg (~15%) for a 10-fold increase in age. Mice are sexually mature by 50 days of age, and over the next 200 days the summed weight of the hippocampi increases by 2.5 mg, or close to 10%. Variation in age among BXD mice accounts for 9% of the variance in hippocampal weight.
     In the BXD animals the weights of both the forebrain and the hippocampus increase with age. However, the effect of age on the hippocampal weight is not explained completely by a general increase in forebrain weight. The slope of the regression of hippocampus against age remains significant (t(247) = 4.0, p < 0.0001) but decreases from 3.9 to 2.2 ± 0.5 mg /log unit age when forebrain weight is added as a cofactor (Fig. 2F). For the sample of F2 progeny there is a significant upward trend in hippocampal weight between 75 and 150 days that amounts to ~0.02 mg/day (R2 = 3%, p < 0.05). In contrast to the BXD sample, there is no greater proportional increase in hippocampal weight of the F2 sample than in the remainder of the forebrain.

Figure 2e: Body Weight and Hippocampal Weight
 

Figure 2f: Body Weight and Hippocampal Weight
 

Figure 2E and F. Age effects on the hippocampus. X and Y axes represent the logarithm of age and hippocampus residuals respectively. (E) Age is significantly correlated with hippocampal weight. (F) Age remains significantly correlated with hippocampal weight after accounting for variation in forebrain minus hippocampus. m is regression line for male; f is regression line for female.

 

Sex and hippocampal weight

Hippocampi of male mice typically weigh 0.50 to 0.60 mg more than those of female mice in the BXD sample and F1 and F2 intercrosses (p < 0.01). Despite differences in body weight, male and female mice have almost precisely the same brain weight—422.7 ± 3.9 versus 422.8 ± 3.0 mg, respectively, for the BXD sample. The difference in hippocampal weight is therefore a relatively specific CNS sex difference. Although statistically significant, the 2% mean difference in hippocampal weight between sexes is relatively trivial given the other numerous sources of variance. This difference should not be characterized as a sexual dimorphism—the overlap between sexes is vast, and sex accounts for only 2% of the total variance in hippocampal weight.

 

Comparison of right and left sides

Measured differences in weights of right and left hippocampi are due to biological differences and technical error. The mean difference between the two sides averages 0.3 mg. Following a correction for small n, this corresponds to a right-left coefficient of variation of merely 2% (Gurland and Tripathi correction; Sokal and Rohlf, 1995). This value sets an upper limit on the magnitude of variation generated by developmental noise and the magnitude of error introduced by fixation and dissection technique. The mean weights of right and left hippocampi across all BXD and F2 mice differ by only 0.05 mg. The right side is on average 0.4% heavier. This tiny difference does reach statistical significance (paired t(4,30) = 2.4; p = 0.012), but it is possible that a slight technical bias is introduced during dissection.

 

Table 2. Correlation matrix of ten parameters for individual BXD animals

 

Multiple regression analysis

To map QTLs that have relatively specific or intense effects on the hippocampus—as opposed to whole brain weight and other extraneous variables—we corrected all original hippocampal weight data by multiple linear regression using the equation

hippocampal weight (bilateral, fixed in mg) = 4.57 + 1.05 (logarithm of age in days) – 0.004(body weight in gm) + 0.048 (brain weight – hippocampal weight in mg) – 0.56 (if female)

This equation accounts for ~50% of the variance among BXD cases (F(4,246) = 59.5). We refer to these regression-corrected values as adjusted hippocampal weights (Table 1, column 3). The adjusted weights represent the weight of the hippocampus after removing predictable effects of sex, brain size, age, and body size. These corrected values were used for mapping QTLs that have comparatively specific effects on the hippocampus. The adjusted hippocampal weight averages 26.1 ± 0.4 mg and ranges from a low of 23.5 mg in BXD32 to a high of 30.8 mg in BXD40. This range extends below the adjusted value of the parental strain DBA/2J (26.4 mg) and far above the value of parental strain C57BL/6J (27.4 mg). A similar multiple linear regression model was used to create a set of adjusted values for the F2 mice (F2 data in text format for all cases):

hippocampal weight (fixed in mg) = 2.81– 1.35 (logarithm of age in days) + 0.017(body weight in gm) + 0.063 (brain weight – hippocampal weight in mg) – 0.52 (if female)

This equation accounts for ~53% of the variance among the F2 sample. The mean of the adjusted weights for the F2 mice was 26.7 ± 0.1. The cumulative probability density plot for the F2 sample after adjustment is unimodal and slightly steeper (leptokurtic) than a normal distribution.

Heritability of hippocampal weight variation

An analysis of variance of BXD mice indicates that the broad-sense heritability for variation in hippocampal weight computed using the method of Hegmann and Possidente (1981) ranges from 44 to 51%. The corresponding estimate of broad sense heritability derived by comparing variance of F1 and F2 animals is 35%. Using the Castle-Wright equation (Wright, 1978) we estimate that a minimum of 3 to 7 loci contribute to the increased variance in F2 progeny compared to their F1 parents. As we show below, we have been able to map two of these factors or QTLs.

Volumetric analysis of the hippocampus

The volume of the hippocampus and several major components were estimated to screen for interesting strain differences as well as to assess the specificity of gene effects. Regions that were measured include the entire hippocampal formation (a measure that corresponds to hippocampal weight; see Methods section), the hippocampus proper (a measurement that excludes the dentate gyrus), the pyramidal cell layer of CA1 through CA3, the dentate gyrus, and the granule cell layer of the dentate gyrus. We find significant differences among BXD and parental strains for each hippocampal component: the total hippocampal complex (R2 = 54%, F(32,135) = 5, p < 0.0001), the hippocampus proper (R2 = 53.6%, F(32,135) = 4.9, p < 0.0001), the dentate gyrus (R2 = 53%, F(32,135) = 4.8, p < 0.0001), the granule cell layer (R2 = 44%, F(32,135) = 3.3, p < 0.0001), and the pyramidal layer (R2 = 51%, F(32,135) = 4.4, p < 0.0001). We also find that the volume of each part correlates significantly with the volume of the total hippocampus. Correlations range from 0.76 to 0.93 (Fig 3). An interesting aspect of this work is that the relative size of the four components differ appreciably among strains. For example, the volume of the dentate gyrus in strains that have almost precisely the same total hippocampal volume ranges from 3.7 mm3 to 4.5 mm3 (Fig 3A, see BXD2 and BXD9). Similarly, the pyramidal cell volume ranges from 1.1 mm3 to 1.7 mm3 among strains with essentially the same hippocampal volume (Fig 3B).

Figure 3A: Structural Correlations
 

Figure 3A. Relations between volume of hippocampus and volume of the dentate gyrus.

Figure 3B: Structural Correlations
 

Figure 3B. Relations between volume of hippocampus and volume of the pyramidal cell layer, including CA1 through CA3.

Figure 3C: Structural Correlations
 

Figure 3C. Relations between volume of hippocampus and the volume of the granule cell layer of the dentate gryus.

Figure 3D: Structural Correlations
 

Figure 3D. Relation between granular layer volume and granular cell density. All volumetric data are corrected for differential shrinkage.

 

Stereological analysis of the dentate granule cell layer

The average packing density of granule cells in the dentate gyrus of BXD RI strains ranges from 750,000 to 1,050,000 cells/mm3. Standard errors of these estimates average 60,000 cells/mm3. Variation in density is not correlated with the total volume of whole brain, hippocampus, dentate gyrus, or even the volume of the granule cell layer itself (Fig. 3D). However, there is a good correspondence between estimates of granule cell number in the two parental strains and the weight of the hippocampus (Table 1). Similarly, there is a significant correlation between estimates of granule cell number and dentate gyrus volume and total hippocampal volume (r = 0.77 and 0.70, respectively). Bilateral estimates of cell number range from just over 1,000,000 cells in BXD1, BXD8, BXD14, and BXD19 to under 700,000 in BXD2, BXD27, BXD29, and BXD34 (Table 1). The magnitude of this difference is as great as differences in retinal ganglion cell number among the same strains of mice (Williams et al., 1998).

 

QTL Analysis of the Mouse Hippocampus

Mapping hippocampal QTLs using BXD mice

Among the BXD strains we detected an excellent match between variation in hippocampal weight and the distribution of B and D alleles at the marker D1Mit145 on distal Chr 1 (Table 1). The average corrected hippocampal weights for 17 BXD strains with a B/B genotype at this locus was 27.05 ± 0.35 mg, whereas that for 18 strains with a D/D genotype was 25.09 ± 0.23 mg. A single B allele in this interval on Chr 1 therefore has an additive effect of approximately +1.0 mg on hippocampal weight. The correlation between hippocampal weight and alleles at D1Mit145 is 0.63 (Table 1), suggesting that as much as 40% of the genetic variance and 15–20% of the total phenotypic variance is generated by a QTL on Chr 1.
     The association between differences in hippocampal weight and alleles on Chr 1 is strong and has a likelihood ratio statistic (LRS) of 19.9 (genome-wide p < 0.05), equivalent to a LOD score of 4.2 (Fig. 4A) between the loci D1Mit145 (89 cM) and D1Mit113 (92 cM). The 2-LOD confidence interval, the chromosomal region in which the QTL is located with a confidence of approximately 95%, extends from D1Mit446 at 70 cM (LRS or 4.5) to D1Mit291 at 101 cM (LRS of 7.8). We have named this locus Hippocampus 1. The abbreviated form or symbol for this QTL is Hipp1.
     Hipp1 is a robust QTL in the sense that it is detected even using the original hippocampal weight values that do not factor out differences in brain weight. However, the linkage statistics are about 10-fold weaker with a peak LRS of 14.7 at 72 cM on Chr 1 near D1Mit194.

 

Figure 4A: Structural Correlations
Figure 4B: Structural Correlations
Figure 4C: Structural Correlations
Figure 4D: Structural Correlations
Figure 4E: Structural Correlations
Figure 4F: Structural Correlations
 

Figure 4. Interval maps of QTLs on Chr 1 (A, B) and Chr 5 (C, D, E, F) in BXD and F2. In each of the six plots, the left axis and bold black line represent values of the likelihood ratio statistic (LRS) computed at 1 cM intervals for hippocampal weight (residuals). The right axis and the thin line represent values for the additive effect of the substitution of a single D allele with a B allele. The X axis represents the entire genetic lengths of Chr 1 (A, B) and Chr 5 (C, D, E, F). Positions of several markers used for mapping both QTLs (e.g., D1Mit200, D5Mit356) are labeled on the X axes. The approximate 2-LOD confidence band is represented by a gray horizontal bar. Plots A, B, C, D are simple interval maps whereas E, F are composite interval maps that control for the Hipp1 interval (D1Mit145 region) of Chr 1.

We have recently generated high density maps for BXD strains that are available as text files for use with Map Manager. These files include some of the phenotype data listed in Table 1 and can be use to verify and extend the main results of this QTL analysis. Please note that Figure 4 was prepared with a more limited set of about 340 BXD genotypes and the LRS plots that are obtained with the new data will have different contours. However, the results and statistics are not affected.


BXD data in Map Manager text format. 36 progeny. Alleles are coded B, D, H, and U (for untyped).
F2 data in Map Manager text format. 183 progeny. Alleles are coded B, D, H, and U.

 


     We controlled for variation generated by the Hipp1 interval on Chr 1 and searched for secondary QTLs affecting hippocampal weight. This procedure uncovered an interval on Chr 5 that is flanked proximally by D5Mit352 (20 cM) and distally by D5Mit356 (41 cM). The LRS score in this interval peaks at 11.0 (Fig. 4E). The genome-wide p is 0.24 and falls short of the level needed to declare a QTL. We subsequently verified the position of this Chr 5 QTL and that on Chr 1 with the F2 intercross.
     The composite mapping procedure was then reversed: we controlled for the Chr 5 interval and remapped Chr 1. LRS values on Chr 1 were still very high, but now in two locations—the original interval near D1Mit145 (LRS of 21.0) and a slightly more proximal region centered at D1Mit45 (59 cM; LRS of 21.6).

Mapping with the F2 Intercross

The hippocampi of F2 progeny were analyzed to extend and validate the analysis of BXD strains. A linkage between weight of the hippocampus was again discovered on distal Chr 1, with a peak LRS of 34.0 between marker D1Mit57 and D1Mit145 (Fig. 4B). This linkage statistic is highly significant with a genome-wide p of ~ 0.0001. The 2-LOD confidence interval on Chr 1 is approximately 25 cM wide and extends from 70 to 95 cM and bracketed by the markers D1Mit103 and D1Mit356. Individuals with B/B, B/D, and D/D genotypes at D1Mit145, the key marker locus identified in the BXD data set, have mean adjusted hippocampal weights of 27.43 ± 0.20 mg, 26.72 ± 0.11 mg, and 26.16 ± 0.16 mg, respectively. Hipp1 is responsible for ~14% of the total phenotypic variance in hippocampal weight. The additive effect of a single B-to-D allele substitution is about 0.64 mg. This is a slightly more modest effect than that noted in the BXD strains at D1Mit145 (27.04 and 25.08 mg for B/B and D/D genotypes).
     The Hipp2 locus on proximal Chr 5 was also confirmed using the F2 intercross. This Chr 5 interval has LRS scores of 12.6 at D5Mit345 and 11.8 at D5Mit352 (Fig. 4D). When variation associated with Hipp1 is controlled the LRS at D5Mit352 increases to between 15 and 20 depending on the particular set of distal Chr 1 loci used as cofactors. For example, using only D1Mit145 gives the most conservative LRS of 15.9 (Fig. 4F), whereas controlling for all markers between 50 and 100 cM (D1Mit80, 387, 103, 57, 145, and 356) raises the composite score to 19.9. Mean adjusted hippocampal weights of the B/B, B/D, and D/D genotypes at D5Mit346 are approximately 27.3, 26.6, and 26.5 mg (ANOVA F(2,174) = 6.0, p = 0.003). The B allele at Hipp2 is fully dominant. When data from BXD and F2 sets are combined, Hipp2 is most likely to map between 15 and 40 cM on Chr 5 (genome-wide p < 0.05).

Test of epistasis between Hipp1 and Hipp2

As judged by their effects in the BXD strains, Hipp1 and Hipp2 alleles sum almost linearly. Those strains with B alleles at both loci (B/B-B/B) have hippocampi that weigh 3.1 mg (12%) more than those of BXD strains with D alleles. The predicted summed effect is 3.5 mg. Residual hippocampal weights for the four possible two-locus genotypes are 1.7 mg (B/B-B/B at Hipp1 and Hipp2), 0.27 (B/B-D/D), –0.25 (D/D-B/B) and –1.4 (D/D-D/D). The same general pattern characterizes F2 animals. Those with B alleles on both Chr 1 and Chr 5 have hippocampi that weigh more than those of F2 animals with D alleles (Table 3). The double heterozygote is within 0.1 mg of the mean of all cases. The predicted linear sum of B alleles at the two loci is 2.1 mg. There is evidence of non-linear interaction between Hipp1 and Hipp2 in the F2 cross—animals with the D/D-B/D two-locus genotype at D1Mit57 and D5Mit346 have a lower mean weight (a residual of —1.13 mg) than even the D/D-D/D animals (—0.51 mg). This suggests that the heterozygous Hipp2 genotype may interact epistatically with Hipp1 to modulate hippocampal size. Another way of expressing this is that the D allele at Hipp2 is overdominant when combined with a homozygous D genotype at Hipp1.

Table 3: Potential interactions between Hipp1 and Hipp2

 

D1Mit57

D5Mit346

N

Mean Residual*

B

B

15

1.057

B

H

17

0.509

B

D

11

0.668

H

B

13

0.459

H

H

47

-0.117

H

D

18

-0.312

D

B

15

-0.013

D

H

21

-1.133

D

D

12

-0.505

*Residual weight values in milligrams for both sides.
Genotypes for the two loci: B = B/B; H = B/D; D = D/D
N is the number of F2 individuals with a particular two-locus genotype.

 

Specificity of QTL action

Hipp1 is defined as a QTL that modulates total hippocampal weight, but it is possible that this locus has more intense effects on one or more parts of the hippocampus, such as the dentate gyrus. To test this possibility we mapped volumetric data from the independent set of BXD data. The top four panels of figure 5 (A–D) demonstrate that volumetric data also show linkage to Chr 1. There is concordance in the locations of peak LRS scores for hippocampal volume (A), pyramidal cell layer volume (B), granule cell layer volume (C), dentate gyrus volume (D), and mossy fiber area (E, our reanalysis of data from Lassalle, et al., 1999). The important point is that in all of these maps the LRS peaks between 65 and 85 cM, the same interval defined as the location of Hipp1 in figure 4 using our more extensive hippocampal weight data set. This is also true of the untransformed and uncorrected cell counts from the dentate gyrus. Peak LRS values for all of these traits have point-wise probabilities less than 0.05, and less then 0.001 in the case of mossy fiber area. There is a similar concordance between the position of Hipp2 on Chr 5 for weight (Fig. 6) and the peak LRS for volume of total hippocampus, the dentate gyrus, the granular layer, and the pyramidal layer. This suggests that Hipp2 also has broad effects on hippocampus structure.

Figure 5: Hipp1 on Chr 1


 

Figure 5. Interval maps for five different hippocampal traits on maps of Chr 1 derived from the BXD strains. The X axes represents the entire genetic length of Chr 1. Y axes represents values for the likelihood ratio statistic (LRS) computed at 1-cM intervals. A, hippocampal volume; B, pyramidal cell layer volume (CA1-CA3); C. Granule cell layer volume of the dentate gyrus; D. Volume of the dentate gyrus (excludes part of the hilus); E. cross-section area of the mossy fiber projection (from Lassalle et al, 1999).

Figure 6: Hipp2 on Chr 5
 

Figure 6. Interval map of QTL on Chr 5. X axes represents the entire genetic length of Chr 5. Y axes represents values for the likelihood ratio statistic (LRS) computed at 1 cM intervals. Conventions as in figure 7.

 


 
      The relation between the two Hipp loci and granule cell number indicates that these two intervals are responsible for generating differences of just over 100,000 cells per dentate gyrus (200,000 bilaterally). Two markers in the Hipp1 and Hipp2 intervals—D1Mit218 and D5Mit197—collectively account for 28% of the variance in granule cell number (Table 1) among the BXD strains means (F(30,2) = 7.1, p = 0.007 and LS of 7.5 at D1Mit218, p = 0.039 at D5Mit197). It is therefore highly likely that Hipp1 has broad effects on all parts of the hippocampus that we measured, affecting both volume and neuron numbers.
     The bottom panels in figure 5 and 6 (5E and 6E) illustrate one of the unique advantages of gene mapping using recombinant inbred strains. We have remapped data published by Lassalle and colleagues (1999) on the area of the mossy fiber projection in BXD strains. Our remapping demonstrates that Hipp1, and possibly Hipp2, act to control the volume of the mossy fiber projection, a finding not unexpected given the wide range in the granule cell populations among BXD strains.

 


 

Increasing specificity by controlling for covariation with cerebellar weight

In our analysis of the F2 intercross we discovered that hippocampal residuals were negatively correlated (r = –0.15) with cerebellar weight residuals that had been computed using similar procedures (Airey et al., 1999). The slope of this relation is modest, and after all corrections have been made there is merely a 0.05 mg increase in hippocampal weight for a 1 mg increase in cerebellar weight (p = 0.04, R2 = 2%). Nonetheless, we decided to include cerebellar weight as a cofactor in computing adjusted hippocampus weights for the F2 set. These latter values were used to map QTLs with improved specificity for the hippocampus in the F2 set. Mapping these cerebellum-corrected scores revealed linkage to a somewhat more proximal region of Chr 1, circa 50–60 cM. The LRS peaks at 27.6 at D1Mit80 located at ~53 cM. A more distal peak that corresponds to Hipp1 is still located near D1Mit57 and D1Mit145 and has a peak LRS of 28.2. Both the more proximal and more distal intervals are significant and have genome-wide p < 0.0005. QTLs defined using F2 intercrosses are often mapped with poor precision (Darvasi, 1998), and it is possible that only a single QTL affecting hippocampal weight maps to Chr 1. However, it is also possible that two linked QTLs are located on Chr 1 between 50 and 100 cM, and that these QTLs are collectively responsible for the genetic variation in hippocampal weight and the high LRS scores on distal Chr 1. This suggestion is supported by the BXD mapping results that are also consistent with a pair of loci located at about 60 and 90 cM. High resolution mapping with a tenth generation advanced intercross that we have recently generated will resolve this intriguing problem.


 

Discussion

Synopsis

We have mapped gene loci with pronounced effects on hippocampal structure to chromosomes 1 and 5 in mouse. Hipp1 and Hipp2 are the first loci known to modulate normal variation in the size of any part of the mammalian forebrain. While effects of QTLs are relatively modest in comparison to the dramatic effects noted in some knockouts and spontaneous mutations, these QTLs have effects that exceed those generated by sex, age, or body weight. It is possible, even likely, that 10–15% shifts in hippocampal weight and neuron number might have functional consequences. This work represents a key step to characterizing genes that normally modulate proliferation, growth, and maturation of the hippocampus. The volumetric and stereological analyses of the Hipp loci demonstrate that these loci have widespread effects on hippocampal structure, including an effect on at least one well-defined hippocampal cell population. The biometric analysis also suggests that age, and to a lesser extent sex, are important factors influencing hippocampus size. A surprising 10% increase in hippocampal weight characterizes sexually mature mice between 50 and 150 days of age.

 

Relations between structure and function

A series of studies has confirmed the existence of heritable differences in the development of the hippocampus (Symons et al., 1988; Lipp et al., 1989). This work has also shown a rough correspondence between hippocampal structural variation and learning ability among different strains of mice, particularly spatial and contextual learning ability (Crusio et al., 1986; Lipp et al., 1989). We find that the hippocampus of C57BL/6J is approximately 5 mg (20%) heavier than that of DBA/2J. At the level of single gene loci, the hippocampus of F2 progeny that have inherited both Hipp1 alleles from C57BL/6J are 1.3 mg heavier than those of progeny that have inherited alleles from DBA/2J. Upchurch and colleagues (1989) examined the inheritance of spatial learning and found that C57BL/6 mice are significantly better at spatial learning than DBA/2 mice. Several studies have also shown that the mossy fiber projection to CA3 is more extensive in C57BL/6 than in DBA/2J mice (Crusio et al., 1989; Schopke et al., 1991). We now add a 20% difference in granule cell number to this list of differences between the parental stains. The involvement of the projection from dentate gyrus to CA3 in learning and memory has been demonstrated in a variety of tasks (Schwegler et al., 1983, Roullet et al., 1990; Bertholet et al., 1991; Flint et al., 1995). It is plausible that allelic differences at Hipp1 and Hipp2 are functionally related to hippocampal-dependent behavioral differences. Of particular interest, Hipp1 maps to the same part of Chr 1 known to harbor one of more QTLs for contextual learning in mice, a hippocampal dependent behavior (Owen et al., 1997; Wehner et al., 1997; Calderone et al., 1997).

 

Behavioral QTLs associated with Hipp1

In recent quantitative genetic studies of contextual conditioning, a learned behavior that depends on the hippocampus, a set of behavioral QTLs have been identified using precisely the same crosses—BXD and B6D2F2—that we have used (Owen et al., 1997; Wehner et al., 1997). One of the most significant and well-replicated QTLs identified in these studies is located on Chr 1 at 80 ± 10 cM (Flint et al., 1995; Calderon et al., 1997; Owen et al., 1997; Wehner et al., 1997). The chance that the position of our strongest morphometric QTL, Hipp1, would match that of the single best behavioral QTL is small—this common 20 cM interval makes up less than 2% of the mouse genome. The D allele at the learning locus on distal Chr 1 is strongly dominant and is associated with greater responsivity to a fearful context. In contrast, alleles at Hipp1 behave in an almost perfectly additive manner and the D allele in this case is associated with a smaller hippocampus. These differences argue either that Hipp1 and the contextual learning QTL are different loci or effects of alleles on anatomical and behavioral traits are not linearly related.

 

Specificity or selectivity of QTL action

What assurance do we have that the QTLs we have mapped have specific effects on hippocampus and not widespread effects on many parts of the CNS? As a first line of defense we mapped hippocampal residualsÑthe signal that remains after controlling for major factors, particularly brain weight, age, and sex. As a second line of defense we have directly mapped QTLs that modulate total brain weight and compared these QTLs with those that modulate hippocampus. Major brain weight QTLs map to Chrs 6, 7, 11, and 14 (Strom, 1999; Williams 2000). Ongoing analysis of the F2 intercross used in the present study indicates that there is a brain weight QTL on Chr 1 and it is possible that Hipp1a has especially intense effects on the hippocampus and less intense and possibly widespread effects on other forebrain derivatives. QTLs will often have graded effects on different CNS components. This leads to our third line of defense: mapping QTLs that control other CNS regions in the same crosses. Hipp1a and Hipp5a do not correspond to any of the olfactory bulb or cerebellar QTLs we have mapped in BXD and the F2 crosses (Williams et al., 2001, Airey et al., 2001). Hipp1a and a cerebellar QTL, Cbs1a, both map to distal Chr 1 but B6 and D2 alleles at these two QTLs have opposite effect polarities arguing against a common gene.
     Regional volumetric and stereological analyses of the hippocampus indicates that Hipp1 and Hipp2 both influence multiple parts of the hippocampus. The correlation between genotype at either locus and phenotypic measures in BXD mice was significant for all the hippocampal traits, including total hippocampus weight, volume of the hippocampus proper, volumes of the dentate gyrus and the granule and pyramidal cell layers, and granule cell population. Even an axon-projection specific phenotypeÑthe area of the mossy fiber projection to CA3 as measured by Timm's stain, significantly correlates to the Hipp1 locus (Lassalle et al., 1999). Thus the QTLs at these loci seem to affect adult hippocampus size selectively, but have shared effects that modulate major parts of the hippocampus.

 

Candidate genes for Hipp1 and Hipp2

Hipp1 maps in the same region as the retinoid x receptor gamma (Rxrg) gene (88 cM) on distal chromosome 1. During brain development, RXRg and retinoid-binding proteins are expressed in CA1–CA3 and the hilus (Zetterstrom et al., 1999). There is evidence that allelic variants at this gene modulate hippocampal development. RXRg is expressed in the adult murine hippocampus (Zetterstrom et al., 1999) and in Rxrg knockout mice, maze learning performance is compromised (Chiang et al., 1998), a finding that suggests that RXRg is at least partly involved in hippocampal spatial learning and memory. It will be interesting to determine 1. if this gene is polymorphic between C57BL/6J and DBA/2J, and 2. if the Rxrg knockout has an effect on hippocampal weight, volume of hippocampal regions, or on neuron number.
      Gene expression profiling is already beginning to identify polymorphic candidate genes in intervals that harbor QTLs. Sandberg, Yasuda, Barlow and colleagues (2000) compared mRNA expression levels of 13,000 genes in the hippocampus of C57BL/6Tac and 129S6/SvEvTac. Their work highlighted an inward rectifying potassium channel gene (Kcnj9 or Girk3) that maps to distal Chr 1 (94 cM) and that is expressed at much lower levels in C57BL/6 than in 129S6/SvEvTac.  Kcnj9 is probably too distal to be a prime candidate for Hipp1, but this example illustrates the growing power of QTL analysis when combined with gene array comparisions between parental strains.
     Both BXD and F2 data sets indicate that Hipp2 maps between 15 and 40 cM on Chr 5. This region contains two candidates—Fgfr3 at 20 cM and a cluster of GABA receptor genes at 40 cM. These genes are expressed in the hippocampus during development and at maturity (Asai et al., 1993; Peters et al., 1993; Mohler et al., 1990; Gaiarsa et al., 1995). FGFR3 and two of its ligands (FGF1 and FGF2) are expressed in hippocampus during development. FGF-dependent mechanisms influence a wide variety of CNS traits, including the survival, growth, and differentiation of hippocampal neurons (Walicke et al., 1986; Sasaki et al., 1992). GABA receptors are potential candidates because GABA is now known to have a trophic role in morphological development of hippocampal neurons (Barbin et al., 1993). With such a roughly defined QTL interval there are of course numerous other genes, many presumably still entirely unknown that may influence hippocampal development.

 

Hippocampal size and sex differences

Polygamous male voles traverse large home ranges in search of mates (Jacobs et al., 1990). Polygamous kangaroo rats also exhibit a sex difference in home range size (Jacobs et al., 1994). Hippocampi of males are typically 10 to 15% larger than those of females. In these polygamous species, an increase in the size of the hippocampus is associated with superior spatial ability (Sherry et al., 1992; Jacobs et al., 1990). We found a smaller sex difference in laboratory mice that conforms to the pattern seen in wild rodent species—the hippocampus of males is about 0.5 mg heavier (bilateral), a difference of only 2%. However, this difference may be highly dependent on strain background. As much as a 15% difference in granule cell number was discovered between males and females of strain LG/J, whereas no significant sex difference was discovered in C58/J mice (Wimer and Wimer, 1985, 1989).

 

Age-related changes in hippocampal weight and adult hippocampal neurogenesis

Neurogenesis occurs in the dentate gyrus of the hippocampus throughout the life of a rodent (Kaplan et al., 1984; Kuhn et al., 1996). Kempermann et al. (1997b) reported that 3-month-old mice produce at least 1 new neuron per 2000 granule cells per day. Assuming complete additivity rather than turnover, this is equivalent to a 10% gain over 200 days. In mice, total granule cell number increases into midlife and then reaches a stable plateau (Kempermann et al., 1998). We have discovered a 10% gain in hippocampal weight over a 200-day period in adult mice. This correspondence is probably fortuitous, particularly so since the dentate gyrus—the only region with notable adult cell production—makes up a small fraction (~20%) of the hippocampal weight.

Kempermann and colleagues (1997b) have shown that the production of new neurons in adults differs significantly among strains of mice, including those we have use in this study. Hipp1 and Hipp2 are both interesting candidates that may influence neurogenesis. It would be feasible to test whether allelic variants at these QTLs modulate levels or kinetics of neurogenesis in the adult mouse. The analysis would rely on identifying 10 to 20 mice homozygous for either B or D alleles at the Hipp1 and Hipp2 intervals and phenotyping these genetically defined animals as adults. Just as is true of Mendelian mutations, an analysis of the functional role of allelic variants at QTLs can precede cloning. In fact, functional and developmental studies of QTLs (Strom and Williams, 1999) can greatly aid in narrowing the set of candidate genes that subsequently need to be analyzed in detail.

 

Acknowledgments

This research project was support by a grant from the National Institute of Neurological Disorders and Stroke (R01 NS35485). The authors thank Drs. Guomin Zhou, Jing Gu, and Xiyun Peng for their assistance in generating, processing, genotyping F2 and BXD mice. We thank Dr. Anand Kulkarni, Toppy Malasri, and David Seecharan for their assistance in the stereological analysis of the dentate gyrus. We thank Dr. Wim Crucio for helpful comments on this work.  

 

References

Abusaad I, MacKay D, Zhao J, Stanford P, Collier DA, Everall IP (1999) Stereological estimation of the total number of neurons in the murine hippocampus using the optical disector. J Comp Neurol 408:560–566.

Airey DC, Lu L, Williams RW (2001) Genetic control of the mouse cerebellum: Identification of quantitative trait loci modulating size and architecture. J Neurosci XX:XXXÐXXX.

Altman J, Das GD (1965) Autoradiographic and histological evidence of postnatal hippocampal neurogenesis in rats. J Comp Neurol 124:319–335.

Asai T, Wanaka A, Kato H, Masana Y, Seo M, Tohyama M (1993) Differential expression of two members of FGF receptor gene family, FGFR-1 and FGFR-2 mRNA, in the adult rat central nervous system. Brain Res Mol Brain Res 17:174–178.

Barber RP, Vaughn JE, Wimer RE, Wimer CC (1974) Genetically-associated variations in the distribution of dentate granule cell synapses upon the pyramidal cell dendrites in mouse hippocampus. J Comp Neurol 156:417–434.

Barbin G, Pollard H, Gaiarsa JL, Benari Y (1993) Involvement of GABAA Receptors in the outgrowth of cultured hippocampal neurons. Neurosci Lett 152:150–154.

Bayer SA (1982) Changes in the total number of dentate granule cells in juvenile and adult rats: a correlated volumetric and 3H-thymidine autoradiographic study. Exp Brain Res 46:315–323.

Belknap JK, Hitzemann R, Crabbe JC, Phillips TJ, Buck KJ, Williams RW (2001) QTL analysis and genome-wide mutagenesis in mice: complementray genetic approaches to the dissection of complex traits. Behav Genet 31:XXXÐXXX.

Bertholet JY, Crusio WE (1991) Spatial and non-spatial spontaneous alternation and hippocampus mossy fiber distribution in nine inbred mouse strains. Behav Brain Res 43:197–202.

Bishop KM, Wahlsten D (1999) Sex and species differences in mouse and rat forebrain commissures depend on the method of adjusting for brain size. Brain Res 815:358–366.

Caldarone B, Saavedra C, Tartaglia K, Wehner JM, Dudek BC, Flaherty L (1997) Quantitative trait loci analysis affecting contextual conditioning in mice. Nat Gen 17:335–337.

Cameron HA, Woolley CS, McEwen BS, Gould (1993) Differentiation of newly born neurons and glia in the dentate gyrus of the adult rat. Neurosci 56:337–344.

Chiang MY, Misner D, Kempermann G, Schikorski T, Giguere V, Sucov HM, Gage FH, Stevens CF, Evans RM (1998) An essential role for retinoid receptors RARb and RXRg in long-term potentiation and depression. Neuron 21:1353–1361.

Churchill GA, Doerge RW (1994) Empirical threshold values for quantitative trait mapping. Genetics 138:963–971.

Crusio WE, Genthner-Grimm G, Schwegler H (1986) A quantitative-genetic analysis of hippocampal variation in the mouse. J Neurogen 3:203–214.

Crusio WE, Schwegler H, van Abeelen JHF (1989) Behavioral responses of novelty and structural variation of hippocampus in mice: II. Multivariate genetic analysis. Behavior Brain Res 32:81–88.

Csernansky JG, Joshi S, Wang L, Haller JW, Gado M, Miller P, Grenander U, Miller MI (1998) Hippocampal morphometry in schizophrenia by high dimensional brain mapping. Proc Natl Acad Sci USA 95:11406–11411.

Darvasi (1998) Experimental strategies for the genetic dissection of complex traits in animals. Nat Gen 18:19–24.

De Bellis MD, Clark DB, Beers SR, Soloff PH, Boring AM, Hall J, Kersh A, Keshavan MS (2000) Hippocampal volume in adolescent-onset alcohol disorders. Am J Psych 157:737–744.

Dietrich WF, Katz H, Lincoln SE (1992) A genetic map of the mouse suitable for typing in intraspecific crosses. Genetics 131:423–447.

Don RH, Cox PT, Wainwright BJ, Baker K, Mattick JS (1991) ‘Touchdown’ PCR to circumvent spurious priming during gene amplification. Nucleic Acids Res 19:4008.

Flint J, Corley R, DeFries JC, Fulker DW, Gray JA, Miller S, Collins AC (1995) A simple genetic basis for a complex psychological trait in laboratory mice. Science 269:1432–1435.

Gaiarsa JL, McLean H, Congar P, Leinekugel X, Khazipov R, Tseeb V, Ben-Ari Y (1995) Postnatal maturation of gamma-aminobutyric acid A and B-mediated inhibition in the CA3 hippocampal region of the rat. J Neurobiol 26:339–349.

Gould E, Caneron HA, McEwen BS (1994) Blockade of NMDA receptors increases cell death and birth in the developing rat dentate gyrus. J Comp Neurol 22:551–565.

Haley CS and Knott SA (1992) A simple regression method for mapping quantitative trait loci in line crosses using flanking markers. Heredity 69:315–324.

Hatton WJ, von Bartheld CS (1999) Analysis of cell death in the trochlear nucleus of chick embryos: calibration of the optical disector counting technique reveals systematic bias. J Comp Neurol 409:169–186.

Hegmann JP, Possidente B (1981) Estimating genetic correlations from inbred strains. Behav Genet 11:103–114.

Jacobs LF, Gaulin SJC, Sherry DF, Hoffman GE (1990) Evolution of spatial cognition: sex-specific patterns of spatial behavior predict hippocampal size. Proc Natl Acad Sci USA 87:6349–6352.

Jacobs LF, Spencer WD (1994) Natural space-use patterns and hippocampal size in kangaroo rats. Brain Behav Evol 44:125–132.

Kaplan MS, Bell DH (1984) Mitotic neuroblasts in the 9-day-old and 11-month-old rodent hippocampus. J Neurosci 4:1429–1441.

Kempermann G, Kuhn HG, Gage FH (1997a) More hippocampal neurons in adult mice living in an enriched environment. Nature 386:493–495.

Kempermann G, Kuhn HG, Gage FH (1997b) Genetic influence on neurogenesis in the dentate gyrus of adult mice. Proc Natl Acad Sci USA 94:10409–10414.

Kempermann G, Kuhn HG, Gage FH (1998) Experience-induced neurogenesis in the senescent dentate gyrus J Neurosci 18:3206–3212.

Kuhn HG, Dickinson-Anson H, Gage FH (1996) Neurogenesis in the dentate gyrus of the adult rat: age-related decrease of neuronal progenitor proliferation. J Neurosci 16:2027–2033.

Kuhn HG, Winkler J, Kempermann G, Thal LJ, Gage FH (1997) Epidermal growth factor and fibroblast growth factor-2 have different effects on neural progenitors in the adult rat brain. J Neurosci 17:5820–5829.

Kulkarni A, Airey DC, Williams RW (2000) Genetic architecture of the mouse retinogeniculate system: a QTL analysis of numerical matching. Soc Neurosci Abst 26: XXXX.

Lassalle JM, Halley H, Milhaud JM, Roullet P (1999) Genetic architecture of the hippocampus mossy fiber subfields in the BXD RI mouse strain series: a preliminary QTL analysis. Behav Genet 29:273–282.

Laird PW, Zijderveld A, Linders K, Rudnicki M, Jaenisch R, Berns A (1991) Simplified mammalian DNA isolation procedure. Nucleic Acids Res 19:4293.

Lipp HP, Schwegler H, Crusio WE, Wolfer DP, Leisinger-Trigona MC, Heimrich B, Driscoll P (1989) Using genetically-defined rodent strains for the identification of hippocampal traits relevant for two-way avoidance behavior: a non-invasive approach. Experientia 45:845–859.

Lipp HP, Wolfer DP (1998) Genetically modified mice and cognition. Curr Opin Neurobiol 8:272–280.

Love JM, Knight AM, McAleer MA, Todd JA (1990) Towards construction of a high resolution map of the mouse genome using PCR-analyzed microsatellites. Nucleic Acids Res 18:4123–4130.

Manly K (1993) A Macintosh program for storage and analysis of experimental genetic mapping data. Mamm Genome, 4: 303–313. in the adult mouse brain. Proc Natyl Acad Sci USA

Manly KF, Olson JM (1999) Overview of QTL mapping software and introduction to Map Manager QT. Mamm Genome, 10:327–334.

Mohler H, Malherbe P, Draguhn A, Richards JG (1990) GABAA Receptors: structural requirements and sites of gene expression in mammalian brain. Neurochem Res 15:199–207.

Owen EH, Logue SF, Rasmussen DL, Wehner JM (1997) Assessment of learning by the Morris water task and fear conditioning in inbred mouse strains and F1 hybrids: implications of genetic background for single gene mutations and quantitative trait loci analyses. Neurosci 80:1087–1099.

Palmer TD, Takahashi J, Gage FH (1997) The adult rat hippocampus contains primordial neural stem cells. Mol Cell Neurosci 8:389–404.

Peters K, Ornitz D, Werner S, Williams L (1993) Unique expression pattern of the FGF receptor 3 gene during mouse organogenesis. Devel Biol 155:423–430

Rikke BA, Johnson TE (1998) Towards the cloning of genes underlying murine QTLs. Mamm Genome 9:963–968.

Rosen GD, Williams RW (2000) Stereological and quantitative genetic analysis of the mouse caudate nucleus. Soc Neurosci Abst 26:311.

Roullet P, Lassalle JM (1990) Genetic variation, hippocampal mossy fibers distribution, novelty reactions and spatial representation in mice. Behav Brain Res 41:61–69.

Routman E, Cheverud J (1994) A rapid method of scoring simple sequence repeat polymorphisms with agarose gel electrophoresis. Mamm Genome 5:187–188.

Sandberg R, Yasuda R, Pankratz DG, Carter TA, Del Rio JO, Wodicka L, Mayford M, Lockhart DJ, Barlow C (2000) Regional and strain-specific gene expression mapping in the adult mouse brain. Proc Natl Acad Sci USA 97: 11039–11043.

Sasaki K, Oomiya Y, Suzuki K, Hanal K, Yagi H (1992) Acidic fibroblast growth factor prevents death of hippocampal CA1 pyramidal cells following ischemia. Neurochem Int 21:397–402.

Schopke R, Wolfer DP, Lipp HP, Leisinger-Trigona MC (1991) Swimming navigation and structural variations of the infrapyramidal mossy fibers in the hippocampus of the mouse. Hippocampus 1:315–328.

Schwegler H, Lipp HP (1983) Hereditary covariations of neuronal circuitry and behavior: correlations between the proportions of hippocampal synaptic fields in the regio inferior and two-way avoidance in mice and rats. Behav Brain Res 7:1–38.

Sherry DF, Jacobs LF, Gaulin SJ (1992) Spatial memory and adaptive specialization of the hippocampus. Trends Neurosci 15:298–303.

Sokal RR, Rohlf FJ (1995) Biometry: the principles and practice of statistics in biological research. 3rd Ed. New York: WH Freeman.

Stanfield BB, Trice JE (1988) Evidence that granule cells generated in the dentate gyrus of adult rats extend axonal projections. Exp Brain Res 72:399–406.

Strom RC, Williams RW (1998) Cell production and cell death in the generation of variation in neuron number. J Neurosci 18:9948–9953.

Strom RC (1999) Genetic analysis of variation in neuron number. Dissertation, University of Tennessee, Memphis.

Symons JP, Davis RE, Marriott JG (1988) Water-maze learning and effects of cholinergic drugs in mouse strains with high and low hippocampal pyramidal cell counts. Life Sciences 42:375–383.

Tanksley SD (1993) Mapping polygenes. Annu Rev Genet 27:205–233.

Taylor BA (1989) Recombinant inbred strains. In: Genetic variants and strains of the laboratory mouse 2nd Ed (Lyon ML, Searle AG, eds), pp 773–796. Oxford: Oxford UP.

Upchurch M, Wehner JM (1989) Inheritance of spatial learning ability in inbred mice: a classical genetic analysis. Behavior Neurosci 103:1251–1258.

von Bartheld CS (1999) Systematic bias in an "unbiased" neuronal counting technique. Anat Rec (New Anat) 257:119–120.

Walicke P, Cowan WM, Ueno N, Baird A, Guillemin R (1986) Fibroblast growth factor promotes survival of dissociated hippocampal neurons and enhances neurite extension. Proc Natl Acad Sci USA 83:3012–3016.

Wehner JM, Radcliffe RA, Rosmann ST, Christensen SC, Rasmussen DL, Fulker DW, Wiles M (1997) Quantitative trait locus analysis of contextual fear conditioning in mice. Nat Gen 17:331-334.

Williams RW (1998) Neuroscience meets quantitative genetics: Using morphometric data to map genes that modulate CNS architecture. In: Morrison J, Hof P (eds) Short course in quantitative neuroanatomy. Society of Neuroscience, Washington DC, pp 66—78.

Williams RW (2000) Mapping genes that modulate mouse brain development: a quantitative genetic approach. In: Mouse brain development. (Goffinet A, Rakic P, eds), pp 21–49. Berlin: Springer.

Williams RW, Strom RC, Goldowitz D (1998a) Natural variation in neuron number in mice is linked to a major quantitative trait locus on Chr 11. J Neurosci 18:138–146.

Williams RW, Rakic P (1988) Three-dimensional counting: an accurate and direct method to estimate numbers of cells in sectioned material. J Comp Neurol 278:344–352, and 281:355.

Williams RW, Airey DC, Kulkarni A, Zhou G, Lu L (2000) Genetic dissection of the olfactory bulb of mice: QTLs on chromosomes 4, 6, 11, and 17 modulate bulb size. Behavior Genetics, in press.

Wimer CC, Wimer RE (1971) Some behavioral differences associated with relative size of hippocampus in the mouse. J Comp Physiol Psychol 76:57–65.

Wimer CC, Wimer RE, Wimer JS (1983) An association between granule cell density in the dentate gyrus and two-way avoidance conditioning in the house mouse. Behav Neurosci 97:844–855.

Wimer RE, Wimer CC (1982) A biometrical–genetic analysis of granule cell number in the area dentata of house mice. Dev Brain Res 2:129–140.

Wimer RE, Wimer CC, Vaughn JE, Barber RP, Balvanz BA, Chernow CR (1976) The genetic organization of neuron number in Ammon’s horns of house mice. Brain Res 118:219–243.

Wimer RE, Wimer CC, Vaughn JE, Barber RP, Balvanz BA, Chernow CR (1978) The genetic organization of neuron number in the granule cell layer of the area dentata in house mice. Brain Res 157:105–122.

Wimer RE, Wimer C (1985) Three sex dimorphisms in the granule cell layer of the hippocampus in house mice. Brain Res 328:105–109.

Wimer CC, Wimer RE (1989) On the sources of strain and sex differences in granule cell number in the dentate area of house mice. Devel Brain Res 48:167–176.

Wright S (1978) Evolution and the genetics of populations, Vol 4. Variability within and among natural populations. Chicago: University of Chicago.

Zetterstrom RH, Lindqvist E, Urquiza AM, Tomac A, Eriksson U, Perlmann T, Olson L (1999) Role of retinoids in the CNS: differential expression of retinoid binding proteins and receptors and evidence for presence of retinoic acid. Eur J Neurosci 11:407–416.

Zhou G, Williams RW (1999) Eye1 and Eye2: Gene loci that modulate eye size, lens weight, and retinal area in mouse. Invest Ophthalmol Vis Sci 40:817–825


 

Since 12 March 2000
Revised 20 Feb 2001

 


   


Neurogenetics at University of Tennessee Health Science Center

Print Friendly Top of Page

Home Page  |  Genome DBs  |  Phenome DBs  |  Publications  |  People & Associates
Mouse Brain Library  |  Related Sites  |  Complextrait.org

Nervenet.org  |   MBL.ORG

Robert W. Williams | Alex Williams © 2002, Nervenet.org modify this page