- Research article
- Open Access
The genetic control of neocortex volume and covariation with neocortical gene expression in mice
BMC Neurosciencevolume 10, Article number: 44 (2009)
The size of the cerebral cortex varies widely within human populations, and a large portion of this variance is modulated by genetic factors. The discovery and characterization of these genes and their variants can contribute to an understanding of individual differences in brain development, behavior, and disease susceptibility. Here we use unbiased stereological techniques to map quantitative trait loci (QTLs) that modulate the volume of neocortex.
We estimated volumes bilaterally in an expanded set of BXD recombinant inbred strains (n = 56 strains and 223 animals) taken from the Mouse Brain Library http://www.mbl.org. We generated matched microarray data for the cerebral cortex in the same large panel of strains and in parental neonates to efficiently nominate and evaluate candidate genes. Volume of the neocortex varies widely, and is a heritable trait. Genome-wide mapping of this trait revealed two QTLs – one on chromosome (Chr) 6 at 88 ± 5 Mb and another at Chr 11 (41 ± 8 Mb). We generated both neonatal and adult neocortical gene expression databases using microarray technology. Using these databases in combination with other bioinformatic tools we have identified positional candidates on these QTL intervals.
This study is the first to use the expanded set of BXD strains to map neocortical volume, and we found that normal variation of this trait is, at least in part, genetically modulated. These results provide a baseline from which to assess the genetic contribution to regional variation in neocortical volume, as well as other neuroanatomic phenotypes that may contribute to variation in regional volume, such as proliferation, death, and number and packing density of neurons
The cerebral cortex is among the most complicated structures in the brain of mammals, and comprises a large but highly variable fraction of the total brain volume. The volume of the human cortex varies by as much as 60% among normal adults. Much of this variation is genetic [1, 2] and has been linked with measures of intelligence [3, 4], as well as with differences in susceptibility to disorders including developmental dyslexia , anxiety-related personality traits , and schizophrenia [7, 8]. Dissecting those genetic variants that modulate normal variation of the cerebral cortex could have an impact on our understanding of cortical development, normal function, and the etiology of several pervasive diseases.
In the present study we use BXD recombinant inbred (RI) strains of mice to investigate the genetic basis of normal variation in the size of the neocortex. This genetic reference population, which is both genetically diverse and phenotypically well-characterized, provides an experimental system to test relations and interactions between normal variation among a potentially unlimited variety of traits . Thus, the scientific community has been able use these mice to systematically accumulate a vast amount of morphometric, behavioral, and physiologic data for each strain. This in turn enables scientists with different interests and expertise to test and verify entire systems of traits, their covariance, and their genetic causes [10–13].
Recently, Beatty and Laughlin  identified a QTL for neocortical volume on chromosome (Chr) 11 using images of BXD cases taken from the Mouse Brain Library (MBL; http://www.mbl.org). Dong et al.  also used the same BXD reference population – and many of the same cases and images from the MBL – and found suggestive QTLs for gray matter volume on Chrs 2, 8, 16, and 19, but not on Chr 11. These studies were constrained by the relatively small number of strains (n = 34) and cases. The recent addition of many new BXD strains and cases to the MBL  has significantly increased the utility of this RI set and improved both the power and precision of QTL mapping. In this experiment, we used unbiased stereology to estimate the neocortical volume in both right and left hemispheres in 54 BXD RI lines, and both parental strains. We mapped QTLs modulating neocortical volume to intervals on Chrs 6 and 11. We exploited a variety of bioinformatic resources, as well as our own neonatal and adult databases of gene expression in the neocortex, to identify potential candidate genes within these QTLs.
Neocortical volume is highly variable
All estimates of volume are fully corrected for case-by-case differences in shrinkage and should be considered close to the original size of these regions in well-fixed tissue. Bilateral neocortical volume of individual mice has a remarkably wide range – from 61.6 to 162.2 mm3 (mean ± SEM = 102.5 ± 0.9 mm3), and is normally distributed (Fig 1A). Strain averages extend from a low of 73.4 ± 4.5 mm3 in BXD30 to a high of 126.9 ± 3.8 mm3 in BXD5 (Fig 1B). The right and left hemisphere neocortical volumes of individual mice ranged from 30.4–82.3 mm3 and from 31.2–79.9 mm3, with means of 51.5 ± 0.45 mm3 and 51.1 ± 0.45 mm3, respectively.
There is no significant left/right asymmetry in the neocortex
We computed an asymmetry coefficient δ asymm by subtracting the volume of the left from right hemisphere and dividing by one half of the total volume. To determine whether there were overall asymmetric biases irrespective of direction, we calculated an absolute asymmetry coefficient (δ abs ) by taking the absolute value of δ asymm . We computed the distribution of both δ asymm and δ abs for the neocortex (mean ± SEM = -0.007 ± 0.003 and 0.039 ± 0.002 mm3, respectively). ANOVA with strain as the independent measure and δ asymm and δ abs as the dependent measures revealed no significant effects of strain (F53,169 < 1, NS in both cases). These results suggest that there is no significant left/right neocortical asymmetry in this population.
Neocortical volume is heritable, but asymmetry is not
In order to determine the heritability of neocortical volume, we computed an ANOVA with strain as the independent measure and bilateral volume as the dependent measure. We found a significant effect of strain (F55,167 = 3.2, P < .001). The strain main effect accounts for 51% of the variance, which provides a reasonable upper bound on the fraction of variance that might be explained by additive and epistatic interactions. (Dominance effects cannot be measured using RI strains because there are no heterozygous genotypes or heterozygote phenotypes to analyze.)
There were no significant main effects for δ asymm (F55,167 < 1, NS) or δ abs (F55,167 = 1.1; explained variance = 22–27%). These results strongly support the notion that there are significant differences between strains on the volumetric measures, but that strain differences in asymmetry are not apparent.
We also computed heritability (h2) using inbred strain data and the adjustment method of Hegmann and Possidente . Neocortical volume is a moderately heritable trait (h2 = 0.29). The broad sense estimate of heritability, which takes into account the large sample sizes of genetically identical members of RI strains, yielded an h2 of 0.67. The heritability factor for δ asymm and δ abs was 0.10 and 0.02, respectively. This low heritability strongly suggests that neither of these traits could be profitably mapped.
Mapping neocortical volume
We mapped bilateral neocortical volume and detected two loci (Fig 2A) with closely matched likelihood ratio statistics (LRS) linkage scores but with opposite allelic effects (that is, B6 vs D2 alleles contributing to greater volume). The first locus on Chr 6 peaks between 88 and 92 Mb (LRS = 13.3, LOD = 2.9). This region corresponds to human Chr 3 at ≈ 134 Mb (3p25.1). The second locus on Chr 11 peaks between 35 and 40 Mb (LRS = 12.7, LOD = 2.6, corresponding to human Chr 5 at ≈ 158 Mb, 3q35.1).
In order to map variation related to volume of neocortex rather than possible confounding covariates, we performed multiple regression to remove all covariance associated with differences in age, sex, plane of section, strain epoch (the original Taylor BXD set versus the newer UTHSC BXD set; see , for an explanation of the history of BXD strains), and non-neocortical brain weight (brain weight – neocortical weight). There were significant effects of non-neocortical brain weight (F1,209 = 10.8, P < .01) and epoch (F1,209 = 8.2, P < .01; Taylor strains = 103.5 ± 1.1 mm3; UTHSC strains = 100.3 ± 1.4). Effects of sex (F1,209 < 1, NS; females = 101.1 ± 1.2 mm3; males = 102.6 ± 1.3 mm3), plane of section (F1,209 = 2.1, NS; horizontal = 101.2 ± 1.3 mm3; coronal = 103.5 ± 1.2 mm3), or age (F1,209 = < 1, NS) were modest and did not reach statistical significance. We regressed neocortical volume for each subject against age, sex, plane of section, epoch, and non-neocortical brain weight and calculated residuals. The use of all variables, even those that do not reach an alpha level of .05, is warranted and removes that variable as a potential confound in mapping. These residuals were used to compute adjusted strain means.
When we re-mapped the variation in adjusted neocortical volume, we detected the same two loci as the original values (Fig 2B), but with a modest increase in the strength of linkage (Chr 6 at 88 ± 5 Mb LRS = 14.5, LOD = 3.1; Chr 11 at 37 ± 5 Mb LRS = 14.4, LOD = 3.1). The D2 allele has a positive effect of about 4 mm3 per allele at the Chr 6 locus, whereas the B6 allele has a similar positive effect at the Chr 11 locus (Additional File 1). We did not detect any epistatic interactions that could account for adjusted neocortical volume. Composite interval mapping of neocortical volume did not reveal any significant or consistent secondary loci for this trait.
Candidate gene analysis
The Chr 6 and Chr 11 loci can be subdivided into large blocks that have common haplotypes in B6, D2, and BXD strains. These blocks are essentially identical by descent (except for any recent mutations) and have low densities of single nucleotide polymorphisms (SNPs). Such regions are less likely to contain polymorphisms that modulate the volume of forebrain derivatives. In contrast, several large blocks have dissimilar haplotypes and much higher densities of SNPs (Figs 2C, D). Genes within these regions have a higher prior probability of containing functional polymorphisms.
Chr 11 candidate analysis
We combined data on SNP density with LRS values to rank positional candidate genes. Of ≈ 70 genes in the Chr 11 interval (34–49 Mb), 17 were located in regions that were highly polymorphic (1 to 6 SNPs per kb). Members of this subset are good positional candidates, irrespective of any role that they may have in brain development or adult forebrain structure (Additional File 2). An additional 9 candidates were located in regions with low SNP densities but had either missense mutations or SNPs in their promoters or UTRs. The restrictive criteria for inclusion were increased by requiring that strong candidates also have moderate to high expression in the neocortex at some stage of development. To apply this filter, we extracted data on expression of all of these genes in the neocortex from the adult and P1 databases. This subset of 26 genes was ranked using data on missense SNPs in these genes and by evidence of local regulatory variation that controls their own expression – so-called cis-QTLs. Those genes with SNPs in Illumina microarray probes that produced high cis-QTLs were excluded from the analysis. Fifteen genes met these criteria: Slit3, Odz2, Gabrg2, Gabra1, Gabrb2, Atp10b, Pttg1, Slu7, Ttc1, Ublcp1, Ebf1, Adam19, Cyfip2, Thg1l, and Itk (Table 1).
Chr 6 candidate analysis
The Chr 6 interval identified here was previously associated with striatal volume, and the identification of candidate genes (using the criteria outlined above), was identical to those recently described . A subset of 10 genes met all the criteria: Htra2, Tia1, Mxd1, Anxa4, Aak1, Nfu1, Nup210, Hdac11, Fbln2, and Slc25a26 (Table 2).
Gene Expression correlations with neocortical volume
We asked which subset of the positional candidates at both QTLs also had differences in expression across BXD strains that co-varied well with our trait. We produced a correlation matrix based on the subset of 47 strains for which we have data of all types. This matrix was then thresholded at a Pearson r = |0.4| and plotted as a graph. There were 7 genes – Gabra1, Gabrb2, Pttg1, Slu7, Ublcp1, Ebf1, and Thg1l – that co-varied with neocortical volume on Chr 11 (Fig 3). We further found confirmatory evidence from the Allen Brain Atlas and GENSAT that these genes are expressed either in the adult brain, the developing brain, or both. In marked contrast, none of the Chr 6 gene expression profiles co-vary with neocortical volume.
Trait correlation analysis
One advantage of employing the BXD RI set is the ability to test for covariation with other phenotypes that have also been studied in this population. We correlated neocortical volume (adjusted and non-adjusted) with the BXD Published Phenotypes Database of GN, which contains a large number of behavioral, anatomic, and physiologic phenotypes gathered from BXD RI sets. Alpha levels were adjusted after permuting neocortical volume and adjusted neocortical volume 20 times each, determining the top 50 correlations with the BXD Published Phenotypes Database for each permutation, and determining the computed alpha level of the top 5% of all correlations. From this analysis, it was determined that computed correlations by GN with P < .001 were significant at an adjusted alpha level of .05.
Both neocortical and adjusted neocortical volume correlated highly with a large number of morphological traits, including volume of dorsal thalamus, dorsal striatum, lateral geniculate nucleus, hippocampus, basolateral amygdala, fimbria, trigeminal tract, nucleus accumbens, and corpus callosum. There were no significant correlations with non-morphological traits.
We systematically explored sources of genetic variation in the size of the neocortex. Using unbiased stereological estimations of this trait in the BXD RI genetic reference population, we found that neocortical volume is highly heritable and variable. There was a 2.6-fold difference between the neocortical volume of the smallest and largest individuals in this study, and a 1.7-fold difference when strain means were compared. The heritability of this trait is consistent with that reported for similar neuroanatomic traits including olfactory bulb, hippocampus, striatum, neocortex, and thalamus [12, 14, 19–24]. Even with a large sample size and a powerful paired t-test design, we did not find any significant asymmetry between the left and right hemispheres of mice. Neocortical volume was associated with QTLs on Chrs 6 and 11. Large microarray data sets were used to efficiently nominate small numbers of viable candidate genes, one or more of which may modulate the volume of the neocortex, and these results supported candidates on the Chr 11 interval as modulating neocortical volume.
There is a lack of anatomic asymmetry in the neocortex
Because previous research has found evidence of volumetric asymmetry in the rodent brain [25–31], we attempted to assess this variable in the BXD RI set. We found that there was little evidence of asymmetry in the neocortex. This was true both for measures that assessed directional (right vs. left) asymmetry and those that assessed magnitude (without considering directionality). These results are consistent with previous reports that failed to reveal neocortical asymmetries in the mouse brain . The results reported here suggest that there are no population level anatomic asymmetries in the neocortex, but do not discount the possibility that anatomic asymmetries at the individual level may correlate with functional asymmetries.
QTLs for neocortical volume
Using images from the MBL, Beatty and Laughlin  reported significant QTLs for neocortical volume on Chr 11 and for non-neocortical brain volume on Chr 19. A more recent study, which again used images from the MBL , did not find any statistically significant QTLs for gray matter volume, although there were suggestive QTLs on Chrs 2, 8, 16, and 19.
We find some support in this study for the Chr 11 QTL reported by Beatty and Laughlin . Although the Chr 11 interval identified by Beatty and Laughlin is more proximal than that described in the current study, the intervals do overlap. In some ways, this was not surprising given that both Beatty and Laughlin and we analyzed brains from the MBL. Yet there were a number of methodological differences that could have worked against this. Our estimations of cerebral volume were performed using unbiased stereology on histologic tissue that was corrected for section thickness, whereas Beatty and Laughlin used planimetric techniques on images from the MBL to estimate volume. Despite these differences in methodology, there was a highly significant correlation (r = 0.78, P < .001) between the neocortical volume estimates of Beatty and Laughlin and our own for the 36 strains common to both studies, and the correlation between the regressed cerebral volume estimates – the variable used to map these traits – is also quite high (r = 0.68,P < .001). This was encouraging given that we regressed out the effects of brain weight (minus the weight of the cortex itself), age, plane of section, epoch, and sex, whereas Beatty and Laughlin regressed out brain weight and logarithm of age only.
There were differences, however, between our results and those published previously. None of the early data sets detect a QTL on Chr 6, and we find no support for any of the QTLs identified by Dong et al. . The most obvious difference among these studies is the inclusion of the new BXD strains . QTL detection and map precision improves with number of strains, whereas increasing the number of subjects within a strain has only marginal effects [33, 34]. We have increased the number of strains by over 50% (34 to 54) in this experiment, and it is likely that our ability to detect QTLs that modulate a smaller percent of the variance of neocortical volume is improved.
Positional candidate genes
Before considering positional candidate genes, it may be useful to review factors that are likely to impact normal variation in brain volume. Cerebral cortical volume can be modified by a variety of developmental factors, as well as those that appear during adulthood. Thus, genes modulating cell proliferation, differentiation, and/or death may contribute to the variation in cerebral cortical volume reported here. In adulthood, cerebral cortical volume can be modified by neurodegenerative changes, some of which may be under genetic control. In the current experiment, age was not a significant predictor of cerebral cortical volume, and none of the mice would be qualified as senescent. When considering potential candidate genes, we therefore concentrated our efforts on those genes known to modulate factors in brain development.
We have recently reported a QTL modulating striatal volume on the same interval on Chr 6 as that in the current report . A detailed analysis of the genes lying within that chromosome (see Table 2) identified a number of potential candidate genes. Analysis of gene expression strain distribution profiles revealed that a number of these genes correlated well with striatal volume. In contrast, none correlated with cerebral cortical volume in the current study. In a complementary manner, none of the Chr 11 candidates reported here correlated with striatal volume.
There are a number of genes in the Chr 11 interval whose functions suggest that they could play a role in modulating neocortical volume. Slit3 is a well-characterized gene that is part of the Slit family of ligands for the Robo receptors , and which are important for neuronal migration, axon outgrowth, and dendritic development [e.g., ]. Moreover, it has recently been demonstrated that this gene is upregulated in the caudal ganglionic eminence, which gives rise to neocortical GABAergic neurons . This is particularly intriguing given the inclusion of three GABA receptor subtype genes (Gabrg2, Gabra1, and Gabrb2) in the QTL interval. Odz2 is expressed in the mouse brain quite early during development , although its specific role is not yet known. Cyfip2 has been shown to be an interactor with FMRP (fragile × mental retardation protein), which has been implicated in Fragile × syndrome, and is highly expressed in neurons throughout the brain .
Of the candidate genes on Chr 11, we found seven whose expression levels were highly correlated with neocortical volume. All are strong cis-QTLs; that is, a sequence variant in or very close to the gene itself controls its expression (Fig 3). Many of these genes have been documented to have roles during brain development. GABA receptors (see above) are obvious candidates given their roles in brain development in general and neocortical functioning in specific [e.g., ]. Pttg1 is upregulated during development of the mouse telencephalon . Likewise, Ublcp1 is upregulated in fast growing or tumor tissues perhaps through its activity as a regulator of the phosphorylation state of RNA Polymerase II . The transcription factor Ebf1 plays a pivotal role in ensuring that a set of thalamic axons reaches the neocortex as well as in the development of striatal neurons [43, 44]. Though the last two transcripts, Slu7 and Thg1l, have not yet been implicated in neural growth, their molecular functions suggest a developmental role. The homolog of Slu7 is well characterized in S. cerevisae as essential during the second catalytic step in the pre-mRNA splicing process . Additionally, Thg1l is important for tRNA priming for protein translation .
Mechanisms that control numbers of neurons and glial cells assigned during development to the neocortex are crucial in understand human behavior, human evolution, and the basis of individual differences. This study can be considered an essential step in a more refined regional genetic anatomy of the several dozen distinct cytoarchitectonic regions that together make up the entire neocortex. The cellular demographics of the neocortex are a strong determinant of functional repertoire. This is most obvious when comparisons are made at a coarse taxonomic level comparing, for example, the massive auditory cortex of echo-locating bats and the complex, highly evolved whisker fields of mice and rats, and the expanded prefrontal cortical regions of humans. But even within a species, pronounced differences in cytoarchitecture are associated with function [47–49]. Now that we have a baseline, we can begin the more painstaking genetic studies of discrete regions of the mouse neocortex and determine if there are QTLs with smaller but regionally intense effects on, for example, barrel cortex, primary visual cortex, or temporal regions.
A number of research strategies that can further refine the QTL interval are available. Higher density maps of the B6 and D2 parental strains will allow a more refined screen for sequence variants. Knockouts of the candidate genes are another potentially productive avenue of investigation. Although none of the Chr 11 candidates have been knocked out, Pttg1, Ebf1, Ublcp1, Slu7, Thg1l, and Gabrb2 have been targeted by the Knockout Mouse Project http://www.komp.org. Additionally, there are a number of forward genetic methods that can be employed to confirm our results. The interval identified here using the BXD RI line could be confirmed in other recombinant inbred lines, including the AXB, BXA, and LXS. The Collaborative Cross – an 8-strain cross to generate hundreds of RI lines, each with more recombinations than existing RI strains – will enhance precision of QTL mapping [50–52]. Recombinant Inbred Intercross (RIX) schemes – whereby selective generation of F1 progeny from recombinant inbred parents allows for the examination of known heterozygous and homozygous intervals – extends the number of genomes that can be phenotyped by n(n-1)/2, where n is the number of original recombinant inbred strains (in our case, 54). The RIX lines have known genomes and quantifiable phenotypes and, therefore, will contribute to mapping traits to their genetic determinants .
This study is the first to use the expanded set of BXD strains to map neocortical volume. With the enhanced precision garnered from the addition of over 20 strains, we found that neocortical volume is a heritable trait with large inter-strain variability. By quantifying cerebral cortical volume for both the right and left hemispheres, we were able to conclude that asymmetry is most likely not genetically modulated, at least in this group of inbred mice. We identified two QTLs, one on Chr 6 and another on Chr 11, that are likely to contain gene variants that modulate neocortical volume. From correlational analysis to gene expression data, we identified a number of candidates on the Chr 11 interval that may modulate this trait. Future research exploiting this expanded BXD genetic reference population will allow us to further dissect other important neuroanatomic phenotypes that may contribute to variation in regional volume, such as proliferation, death, and number and packing density of neurons.
All experiments were approved by the Institutional Animal Care and Use Committees at Beth Israel Deaconess Medical Center and University of Tennessee Health Science Center.
All histologic data for this study were obtained from The Mouse Brain Library (MBL) – a physical and Internet resource that contains high-resolution images of histologically processed slides from over 2900 adult mouse brains http://www.mbl.org with roughly balanced numbers of male and female specimens . The ages ranged from 21–694 days of age (mean ± SEM = 103 ± 5), with most of the cases ranging from 50 – 120 days. Mice were obtained from either the Jackson Laboratory (Bar Harbor, ME) or the University of Tennessee Health Science Center (UTHSC) as detailed previously . All procedures were approved by institutional animal care and use committees and conform to NIH guidelines for humane treatment of animals. Mice were deeply anesthetized with Avertin (0.8 ml i.p.) and transcardially perfused with 0.9%saline, followed by fixative (1.25% glutaraldehyde/1.0% paraformaldehyde in phosphate buffer), and their brains removed and weighed. After variable post-fixation times, the brains were embedded in 12% celloidin and sliced in either a coronal or horizontal plane at a width of approximately 30 μm. Actual section thickness was determined by direct examination of 10 sections for each brain using an ×100 oil immersion objective and a z-axis micrometer.
Estimation of neocortical volume
The volume of the neocortex was estimated in 223 mice (114 female and 109 male) from 54 BXD RI strains and the two parental strains (C57BL/6J and DBA/2J, abbreviated B6 and D2, respectively) by one of us (SMG) using a computer controlled microscope (Nikon E800, Nikon, Inc., Melville, NY) and Stereo Investigator (MBF Biosciences, Williston, VT). For most strains, 4 mice were used with the exception of BXD2 (N = 3), BXD15 (N = 5), and BXD83 (N = 3). Neocortex was parcellated by the criteria of the Allen Brain Atlas (mouse.brain-map.org), and includes agranular insular (AId, AIv, AIp), auditory (AUDv, AUDd, AUDv), cingulate (ACAd, ACAv), ectorhinal (ECT), gustatory (GU), infralimbic (ILA), motor (MOs, MOp), orbital (ORBl, ORBvl, ORBm), parietal (PPLp), prelimbic (PL), perirhinal (PERI), retrosplenial (RSPv, RSPd, RSPagl), somatosensory (SSp, SSs), temporal (TEa), visceral (VISC), and visual (VISam, VISp, VISal) cortices. Volume was estimated by point counting using Cavalieri's method. Grid spacing for horizontal (N = 110) and coronal (N = 113) sections was 800 μm and 400 μm, respectively. For each brain, two interleaved 1-in-10 series of sections were examined, representing every fifth section. For horizontal sections, volume was estimated by measuring every tenth section on each slide, whereas every thirtieth section was measured for coronal sections (six to ten sections per slide). Volume was independently estimated for both the right and left hemispheres, and the sum of these measures estimated total neocortical volume. In cases where there were missing or damaged sections, a piece-wise parabolic estimation was used . Final volume estimates were individually corrected for histological shrinkage by determining the previously computed ratio between the brain volume at fixation (brain weight) and that after processing. Neocortical volume was blindly re-measured on 10 slides to assess intra-observer reliability. The experimenter was blind with respect to strain and sex.
Intra-observer reliability was high for estimation of neocortical volume. The percentage difference between the original and repeated volume estimations ranged from 0–5% and the average difference was 0.42%. A correlation coefficient between the two measurements was highly significant (r = 0.99) indicating that technical error at this level of the analysis contributes little to case variation or strain variation. A paired t-test confirmed that the difference between the first and second estimations was not significant (t < 1, NS).
Data were analyzed using standard ANOVA and multiple regression techniques (JMP, SAS Institute, Cary, NC). QTL analysis was performed using the WebQTL module of GeneNetwork (GN, http://www.genenetwork.org). This online resource includes all known morphometric data for the BXD strains, several neocortical transcriptome data sets, high density marker maps based on approximately 3795 fully informative markers distributed on all chromosomes except Chr Y , and a database containing ≈ 8.3 million SNPs taken from dbSNP . WebQTL incorporates three common mapping methods: (1) simple interval mapping, (2) composite interval mapping, and (3) a scan for two-locus epistatic interactions .
We exploited an expression data set available at http://www.genenetwork.org that was generated with support from the CHDI Foundation. This data set estimated steady-state mRNA levels in the neocortex of 73 strains of adult (P60 ± 5) mice, of which 44 are common to the morphometric data described here. Strains for which we have morphometric data but no array data include BXD22, BXD24, BXD25, BXD30, BXD35, BXD48, BXD63, BXD83, BXD85 and BXD92. We specifically used the GeneNetwork database named "HQF BXD Neocortex ILM6.1 (Feb08) RankInv," which can be accessed in the main search page by selecting Species = Mouse, Group = BXD, Type = Neocortex mRNA. Data were generated using the Illumina Sentrix Mouse-6.1 microarray. This array estimates expression for a great majority of mouse genes with confirmed protein products and consists of sets of ≈ 46,000 unique 50-nucleotide-long probe sequences. Like other array data in GeneNetwork , the original Illumina bead array data (rank invariant transform) were logged and re-centered to a mean of 8 units and a standard deviation of 2 units – essentially a Z transform of the data. For complete metadata on the adult neocortex data set, including quality control procedures, error-checking, and normalization see http://www.genenetwork.org/dbdoc/HQFNeoc_0208_RankInv.html.
Comparing adult morphometric data with adult cortical expression data is potentially insufficient since it is possible that the transient expression of a small number of genes during development may leave a lasting imprint on differences in volume maturity [but see ]. For this reason, we generated a companion expression data set of neonatal (P1) neocortical expression data for the two parental strains, B6 and D2 (NIH Neuroscience Microarray Consortium http://arrayconsortium.tgen.org, rosen-illu-mouse- 588443). Data were generated using the Illumina MouseWG-6 v2 Expression Beadchip. We have exploited this companion developmental data to 1) test whether genes with expression differences that co-vary with neocortical volume at maturity are also expressed during a key stage of development in the neocortex, and 2) to test whether any genes in QTL intervals have high expression only during development.
All genome coordinates in this paper are given using the mouse genome assembly of February 2006 (UCSC Genome Browser release mm8, NCBI Build 36). These position coordinates differ slightly (usually less than 1–3 Mb) from mm9 and NCBI Build 37.1.
To evaluate candidate genes and to study other molecular co-variates of neocortical volume we correlated our estimated volumes with the P1 and P60 neocortical transcriptome data sets. Co-variation networks were constructed using on-line tools in GN. We additionally correlated our neocortical volume data set with a database of over 1000 previously published and unpublished BXD traits http://www.genenetwork.org/dbdoc/BXDPublish.html.
On-line data access
Phenotypes for the BXD strains generated as part of this study are all available at Genenetwork.org using the accession numbers:
Neocortical volume (mm3): GN BXD Phenotypes Trait ID: 10995
Neocortical volume adjusted (mm3): GN BXD Phenotypes Trait ID: 10997
Cerebral cortex volume (mm3, Beatty and Laughlin, 2006): GN BXD Phenotypes Trait ID: 10992
Neocortex expression data exploited in this study are directly accessible at http://www.genenetwork.org (adult) or arrayconsortium.tgen.org (P1 data set).
Mouse Brain Library
Quantitative Trait Locus
single nucleotide polymorphism.
Pennington BF, Filipek PA, Lefly D, Chhabildas N, Kennedy DN, Simon JH, Filley CM, Galaburda A, DeFries JC: A twin MRI study of size variations in human brain. J Cogn Neurosci. 2000, 12 (1): 223-232.
Baaré WF, Hulshoff Pol HE, Boomsma DI, Posthuma D, de Geus EJ, Schnack HG, van Haren NE, van Oel CJ, Kahn RS: Quantitative genetic modeling of variation in human brain morphology. Cereb Cortex. 2001, 11 (9): 816-824.
Thoma RJ, Yeo RA, Gangestad SW, Halgren E, Sanchez NM, Lewine JD: Cortical volume and developmental instability are independent predictors of general intellectual ability. Intelligence. 2005, 33 (1): 27-38.
Haier RJ, Jung RE, Yeo RA, Head K, Alkire MT: Structural brain variation and general intelligence. Neuroimage. 2004, 23 (1): 425-433.
Sandu AL, Specht K, Beneventi H, Lundervold A, Hugdahl K: Sex-differences in grey-white matter structure in normal-reading and dyslexic adolescents. Neurosci Lett. 2008, 438 (1): 80-84.
Yamasue H, Abe O, Suga M, Yamada H, Inoue H, Tochigi M, Rogers M, Aoki S, Kato N, Kasai K: Gender-common and -specific neuroanatomical basis of human anxiety-related personality traits. Cereb Cortex. 2008, 18 (1): 46-52.
Szeszko PR, Hodgkinson CA, Robinson DG, Derosse P, Bilder RM, Lencz T, Burdick KE, Napolitano B, Betensky JD, Kane JM, et al.: DISC1 is associated with prefrontal cortical gray matter and positive symptoms in schizophrenia. Biol Psychol. 2008, 79 (1): 103-110.
Reiss AL, Faruque F, Naidu S, Abrams M, Beaty T, Bryan RN, Moser H: Neuroanatomy of Rett syndrome: a volumetric imaging study. Ann Neurol. 1993, 34 (2): 227-234.
Rosen GD, La Porte NT, Diechtiareff B, Pung CJ, Nissanov J, Gustafson C, Bertrand L, Gefen S, Fan Y, Tretiak O, et al.: Informatics center for mouse genomics: The dissection of complex traits of the nervous system. Neuroinformatics. 2003, 1 (4): 327-342.
Chesler EJ, Lu L, Shou S, Qu Y, Gu J, Wang J, Hsu HC, Mountz JD, Baldwin NE, Langston MA, et al.: Complex trait analysis of gene expression uncovers polygenic and pleiotropic networks that modulate nervous system function. Nat Genet. 2005, 37 (3): 233-242.
Chesler EJ, Lu L, Wang J, Williams RW, Manly KF: WebQTL: rapid exploratory analysis of gene expression and genetic networks for brain and behavior. Nat Neurosci. 2004, 7 (5): 485-486.
Mozhui K, Hamre KM, Holmes A, Lu L, Williams RW: Genetic and structural analysis of the basolateral amygdala complex in BXD recombinant inbred mice. Behav Genet. 2007, 37 (1): 223-243.
Yang RJ, Mozhui K, Karlsson RM, Cameron HA, Williams RW, Holmes A: Variation in mouse basolateral amygdala volume is associated with differences in stress reactivity and fear learning. Neuropsychopharmacology. 2008, 33 (11): 2595-2604.
Beatty J, Laughlin R: Genomic regulation of natural variation in cortical and noncortical brain volume. BMC Neuroscience. 2006, 7 (1): 16.
Dong H, Martin MV, Colvin J, Ali Z, Wang L, Lu L, Williams RW, Rosen GD, Csernansky JG, Cheverud JM: Quantitative trait loci linked to thalamus and cortex gray matter volumes in BXD recombinant inbred mice. Heredity. 2007, 99 (1): 62-69.
Peirce JL, Lu L, Gu J, Silver LM, Williams RW: A new set of BXD recombinant inbred lines from advanced intercross populations in mice. BMC Genet. 2004, 5 (1): 7.
Hegmann JP, Possidente B: Estimating genetic correlations from inbred strains. Behav Genet. 1981, 11 (2): 103-114.
Rosen GD, Pung CJ, Owens CB, Caplow J, Kim H, Mozhui K, Lu L, Williams RW: Genetic modulation of striatal volume by loci on Chrs 6 and 17 in BXD recombinant inbred mice. Genes Brain Behav. 2009, 8 (3): 296-308.
Williams RW, Lu L, Kulkarnik A, Zhou G, Airey DC: Genetic dissection of the olfactory bulbs of mice: QTLs on chromosomes 4, 6, 11, and 17 modulate bulb size. Behav Genet. 2001, 31: 61-77.
Yang J, Zhu J, Williams RW: Mapping the genetic architecture of complex traits in experimental populations. Bioinformatics. 2007, 23 (12): 1527-1536.
Lu L, Airey DC, Williams RW: 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. 2001, 21 (10): 3503-3514.
Martin MV, Dong H, Vallera D, Lee D, Lu L, Williams RW, Rosen GD, Cheverud JM, Csernansky JG: Independent quantitative trait loci influence ventral and dorsal hippocampal volume in recombinant inbred strains of mice. Genes Brain Behav. 2006, 5 (8): 614-623.
Peirce JL, Chesler EJ, Williams RW, Lu L: Genetic architecture of the mouse hippocampus: identification of gene loci with selective regional effects. Genes Brain Behav. 2003, 2 (4): 238-252.
Rosen GD, Williams RW: Complex trait analysis of the mouse striatum: Independent QTLs modulate volume and neuron number. BMC Neuroscience. 2001, 2 (1): 5.
Rosen GD, Sherman GF, Galaburda AM: Interhemispheric connections differ between symmetrical and asymmetrical brain regions. Neuroscience. 1989, 33: 525-533.
Rosen GD, Sherman GF, Galaburda AM: Ontogenesis of neocortical asymmetry: A [3H]thymidine study. Neuroscience. 1991, 41 (2–3): 779-790.
Zilles K, Dabringhaus A, Geyer S, Amunts K, Qü M, Schleicher A, Gilissen E, Schlaug G, Steinmetz H: Structural asymmetries in the human forebrain and the forebrain of non-human primates and rats. Neurosci Biobehav Rev. 1996, 20 (4): 593-605.
Rosen GD, Sherman GF, Galaburda AM: Neuronal subtypes and anatomic asymmetry: Changes in neuronal number and cell-packing density. Neuroscience. 1993, 56 (4): 833-839.
Verstynen T, Tierney R, Urbanski T, Tang A: Neonatal novelty exposure modulates hippocampal volumetric asymmetry in the rat. NeuroReport. 2001, 12 (14): 3019-3022.
Caparelli-Daquer EM, Schmidt SL: Morphological brain asymmetries in male mice with callosal defects due to prenatal gamma irradiation. Int J Dev Neurosci. 1999, 17 (1): 67-77.
Rosen GD, Sherman GF, Mehler C, Emsbo K, Galaburda AM: The effect of developmental neuropathology on neocortical asymmetry in New Zealand Black mice. Int J Neurosci. 1989, 45: 247-254.
Rosen GD, Sherman GF, Emsbo K, Mehler C, Galaburda AM: The midsagittal area of the corpus callosum and total neocortical volume differ in three inbred strains of mice. Exp Neurol. 1990, 107: 271-276.
Belknap JK: Effect of within-strain sample size on QTL detection and mapping using recombinant inbred mouse strains. Behav Genet. 1998, 28 (1): 29-38.
Crusio WE: A note on the effect of within-strain sample sizes on QTL mapping in recombinant inbred strain studies. Genes Brain Behav. 2004, 3 (4): 249-251.
Chedotal A: Slits and their receptors. Adv Exp Med Biol. 2007, 621: 65-80.
Whitford KL, Marillat V, Stein E, Goodman CS, Tessier-Lavigne M, Chedotal A, Ghosh A: Regulation of cortical dendrite development by Slit-Robo interactions. Neuron. 2002, 33 (1): 47-61.
Willi-Monnerat S, Migliavacca E, Surdez D, Delorenzi M, Luthi-Carter R, Terskikh AV: Comprehensive spatiotemporal transcriptomic analyses of the ganglionic eminences demonstrate the uniqueness of its caudal subdivision. Mol Cell Neurosci. 2008, 37 (4): 845-856.
Zhou XH, Brandau O, Feng K, Oohashi T, Ninomiya Y, Rauch U, Fassler R: The murine Ten-m/Odz genes show distinct but overlapping expression patterns during development and in adult brain. Gene Expr Patterns. 2003, 3 (4): 397-405.
Schenck A, Bardoni B, Langmann C, Harden N, Mandel JL, Giangrande A: CYFIP/Sra-1 controls neuronal connectivity in Drosophila and links the Rac1 GTPase pathway to the fragile × protein. Neuron. 2003, 38 (6): 887-898.
Levitt P, Harvey JA, Friedman E, Simansky K, Murphy EH: New evidence for neurotransmitter influences on brain development. Trends Neurosci. 1997, 20 (6): 269-274.
Tarabykin V, Britanova O, Fradkov A, Voss A, Katz LS, Lukyanov S, Gruss P: Expression of PTTG and prc1 genes during telencephalic neurogenesis. Mech Dev. 2000, 92 (2): 301-304.
Zheng H, Ji C, Gu S, Shi B, Wang J, Xie Y, Mao Y: Cloning and characterization of a novel RNA polymerase II C-terminal domain phosphatase. Biochem Biophys Res Commun. 2005, 331 (4): 1401-1407.
Lobo MK, Karsten SL, Gray M, Geschwind DH, Yang XW: FACS-array profiling of striatal projection neuron subtypes in juvenile and adult mouse brains. Nat Neurosci. 2006, 9 (3): 443-452.
Garel S, Yun K, Grosschedl R, Rubenstein JL: The early topography of thalamocortical projections is shifted in Ebf1 and Dlx1/2 mutant mice. Development. 2002, 129 (24): 5621-5634.
Frank D, Guthrie C: An essential splicing factor, SLU7, mediates 3' splice site choice in yeast. Genes Dev. 1992, 6 (11): 2112-2124.
Jahn D, Pande S: Histidine tRNA guanylyltransferase from Saccharomyces cerevisiae. II. Catalytic mechanism. J Biol Chem. 1991, 266 (34): 22832-22836.
Airey DC, Wu F, Guan M, Collins CE: Geometric morphometrics defines shape differences in the cortical area map of C57BL/6J and DBA/2J inbred mice. BMC Neurosci. 2006, 7: 63.
Heimel JA, Hermans JM, Sommeijer JP, Levelt CN: Genetic control of experience-dependent plasticity in the visual cortex. Genes Brain Behav. 2008, 7 (8): 915-923.
Jan TA, Lu L, Li CX, Williams RW, Waters RS: Genetic analysis of posterior medial barrel subfield (PMBSF) size in somatosensory cortex (SI) in recombinant inbred strains of mice. BMC Neurosci. 2008, 9: 3.
Churchill GA, Airey DC, Allayee H, Angel JM, Attie AD, Beatty J, Beavis WD, Belknap JK, Bennett B, Berrettini W, et al.: The Collaborative Cross, a community resource for the genetic analysis of complex traits. Nat Genet. 2004, 36 (11): 1133-1137.
Valdar W, Flint J, Mott R: Simulating the collaborative cross: power of quantitative trait loci detection and mapping resolution in large sets of recombinant inbred strains of mice. Genetics. 2006, 172 (3): 1783-1797.
Chesler E, Miller D, Branstetter L, Galloway L, Jackson B, Philip V, Voy B, Culiat C, Threadgill D, Williams R, et al.: The Collaborative Cross at Oak Ridge National Laboratory: developing a powerful resource for systems genetics. Mammalian Genome. 2008, 19 (6): 382-389.
Threadgill DW, Hunter KW, Williams RW: Genetic dissection of complex and quantitative traits: from fantasy to reality via a community effort. Mamm Genome. 2002, 13 (4): 175-178.
Rosen GD, Harry JD: Brain volume estimation from serial section measurements: A comparison of methodologies. J Neurosci Methods. 1990, 35: 115-124.
Shifman S, Bell JT, Copley RR, Taylor MS, Williams RW, Mott R, Flint J: A high-resolution single nucleotide polymorphism genetic map of the mouse genome. PLoS biology. 2006, 4 (12): e395.
Frazer KA, Eskin E, Kang HM, Bogue MA, Hinds DA, Beilharz EJ, Gupta RV, Montgomery J, Morenzoni MM, Nilsen GB, et al.: A sequence-based variation map of 8.27 million SNPs in inbred mouse strains. Nature. 2007, 448 (7157): 1050-1053.
Wang J, Williams RW, Manly KF: WebQTL: Web-based complex trait analysis. Neuroinformatics. 2003, 1 (4): 299-308.
Zapala MA, Hovatta I, Ellison JA, Wodicka L, Del Rio JA, Tennant R, Tynan W, Broide RS, Helton R, Stoveken BS, et al.: Adult mouse brain gene expression patterns bear an embryologic imprint. Proc Natl Acad Sci USA. 2005, 102 (29): 10357-10362.
This work was supported, in part, by NS052397 from NINDS to GDR, by DA021131 from NIDA, NIAAA, and NIMH to GDR and RWW, U01AA014425 to LL, and by U01AA13499, U01CA105417, and U24 RR021760 for support of GeneNetwork. SMG was supported by the Research Science Institute of the Center for Excellence in Education and the Harvard College Research Program. The CHDI Foundation supported the generation of expression data for neocortex. We thank the NIH Microarray Consortium and the Yale/NIH Microarray Center for performing the P1 Illumina microarray analyses. We thank Dr. Xusheng Wang for help annotating the content of the Illumina array. We thank Stephanie Chin, Bradford Craigen, Annelise Joseph, Monica Larimer, Christopher Pung, and Nataliya Ryzhenko for their help in generating the neocortical transcriptome databases.
SMG was involved in the design or the experiment, carried out the data acquisition, participated in the data analysis and in the drafting of the manuscript. LL provided the animals used to generate portions of the adult transcriptome database and participated in editing the manuscript. RWW curated the GeneNetwork databases, participated in the data analysis and in the editing of the manuscript. GDR conceived the experiment, supervised the data analysis and drafting of the manuscript. All authors read and approved the final manuscript.