Geometric morphometrics defines shape differences in the cortical area map of C57BL/6J and DBA/2J inbred mice
© Airey et al; licensee BioMed Central Ltd. 2006
Received: 12 June 2006
Accepted: 15 September 2006
Published: 15 September 2006
We previously described planar areal differences in adult mouse visual, somatosensory, and neocortex that collectively discriminated C57BL/6J and DBA/2J inbred strain identity. Here we use a novel application of established methods of two-dimensional geometric morphometrics to examine shape differences in the cortical area maps of these inbred strains.
We used Procrustes superimposition to align a reliable set of landmarks in the plane of the cortical sheet from tangential sections stained for the cytochrome oxidase enzyme. Procrustes superimposition translates landmark configurations to a common origin, scales them to a common size, and rotates them to minimize an estimate of error. Remaining variation represents shape differences. We compared the variation in shape between C57BL/6J and DBA/2J relative to that within each strain using a permutation test of Goodall's F statistic. Significant differences in shape in the posterior medial barrel subfield (PMBSF), as well as differences in shape across primary sensory areas, characterize the cortical area maps of these common inbred, isogenic strains.
C57BL/6J and DBA/2J have markedly different cortical area maps, in both size and shape. These differences suggest polymorphism in genetic factors underlying cortical specification, even between common isogenic strains. Comparing cortical phenotypes between normally varying inbred mice or between genetically modified mice can identify genetic contributions to cortical specification. Geometric morphometric analysis of shape represents an additional quantitative tool for the study of cortical development, regardless of whether it is studied from phenotype to gene or gene to phenotype.
Similar differences may be present between standard inbred strains of mice. Such differences, although likely more subtle than those induced by Hamasaki and colleagues, suggest polymorphism in genes underlying morphogenesis of the cerebral cortex. Because of the extensive changes in the cortical area map of transgenic Emx2 mice, the quantitative descriptions used by Hamasaki et al. were adequate in portraying the visible differences in area and shape . When differences in the cortical area map are more subtle, then more sensitive, multivariate statistical methods are required. We previously demonstrated significant differences in the size of both V1 (12%) and the PMBSF (10%) between the inbred strains C57BL/6J and DBA/2J . A multivariate logit model predicted strain identity with 90% correct classification rate . We were unable, however, to demonstrate any reliable differences in landmark positions in the cortical area map using the methods described by Hamasaki et al., . In this paper, we use established statistical methods to analyze shape, called geometric morphometrics (see [7–9] for introductions to shape analysis, and Methods for more details).
The analysis of shape has been successfully used to investigate neuroanatomical differences, but mostly within the context of human brain imaging data. For example, using anatomical landmarks from two dimensional magnetic resonance image (MRI) scans of the human midsagittal plane, Gharaibeh et al., , examined shape changes in the configuration of landmarks following the first episode of schizophrenia. Bookstein et al., , used other geometric morphometric methods to describe shape differences in the corpus callosum of fetal alcohol affected human brains, again from MRI scans. Computational anatomy methods have also been developed to describe shape changes in three dimensional data (e.g., [12, 13]). The lissencephalic mouse cortex, in stained tangential sections, lends itself well to the application of established geometric morphometric methods for landmarks in two dimensions. There is a long history of using tangential cortical tissue sections in neuroanatomy, but the application of geometric morphometric methods to this histological preparation is novel.
Using two sets of landmark configurations, one local (PMBSF) and one more global (V1, S1, A1), we sought to identify differences in shape of primary sensory cortex between C57BL/6J and DBA/2J inbred strains. We collected landmark data from cytochrome oxidase stained sections of flattened cortex and compared landmark configurations using the Procrustes transformation, a standard method in the analysis of shape (see Methods). Landmarks included the positions of barrels in the PMBSF and minimal landmarks for the triangular shapes of primary visual cortex (V1) and primary somatosensory cortex (S1), and the position of the centroid of primary auditory cortex (A1) (see Figure 2). A single landmark for A1 is used, because A1 is ovoid in shape and lacks reliable outline landmarks, as in V1 or S1. Similarly, centroids were used as landmarks for each ovoid shaped barrel. The discovery of significant differences between these strains encourages further use of geometric morphometrics to identify genetic or transcript correlates of shape, either by classical trait QTL mapping (quantitative trait locus mapping) or expression QTL mapping [18, 19]. Geometric morphometrics may also find use in describing the consequences of experimental genetic manipulations affecting the mouse cortical area map. Regardless of whether the cortical area map is studied from phenotype to gene or gene to phenotype, the analysis of cortical area map shape will help increase our understanding of cerebral cortical morphogenesis.
Local differences: PMBSF
Deformation grids and vector plots
Thin-plate spline deformation grids for the Procrustes transformed barrel landmarks are shown for C57BL/6J, the consensus (average) of the two strains, and for DBA/2J (Figure 4A). Immediately noticeable from the general outline of the grids and spacing between barrel landmarks, is an anteroposterior expansion in C57BL/6J (taller grid) and the corresponding anteroposterior compression in DBA/2J (shorter grid). Differences in anteroposterior spacing between barrel landmarks is particularly apparent in barrel rows A and B. Less obvious, but still noticeable from the general grid outline, are mediolateral effects. The overall shape of the C57BL/6J grid in 4A is wider posteriorly, with the opposite effect apparent in DBA/2J – in other words, there is posterior lateral expansion and anterior lateral compression in the PMBSF of C57BL/6J. In addition to general anteroposterior and mediolateral differences in PMBSF shape, there are also focal anteroposterior and mediolateral effects. Along the A-P axis, there is straightening of the anterior ends of rows C, D, and E in DBA/2J (see lower left of the DBA/2J panel in 4A versus lower left of C57BL/6J). Along the M-L axis, rows A, B, and C are relatively compressed in C57BL/6J, while rows C, D, and E are relatively compressed in DBA/2J. Each of these general and focal effects described in the deformation grids can actually be discerned in the scatterplot, though it requires looking back and forth between the plots. Conversely, the comments made above for the scatterplot also are evident in the deformation grids. An alternative representation of the noted shape changes above is shown in 4B. These plots represent the change in barrel position by vectors.
Uniform components of shape
The plots in Figure 4A and 4B represent shape changes from the total variation shown in Figure 3. Shape variation can be decomposed into uniform and non-uniform components , and both the deformation grids and the vector plots in 4A or 4B illustrate the combined effects of uniform and non-uniform components of shape. Uniform components of shape do not bend the parallel grid lines in deformation plots such as those shown in Figure 4A. Two kinds of uniform shape differences are compression and shearing. Compression and shearing, for example, would make a rectangle more narrow or would make a rectangle into a parallelogram, respectively (see Methods for additional comments).
Figure 4C represents the uniform component of shape change for C57BL/6J and DBA/2J mean configurations. Regression of the partial warps for the X and Y uniform shape variables against the independent variables (here, strain) gives separate tests. For PMBSF, these test probabilities were 0.0005 and 0.0008, respectively. In addition to the anteroposterior compression noted above (the height of the DBA/2J parallelogram is less than in C57BL/6J), shearing is also evident in these plots (the longer diagonal of the DBA/2J plot lies on the A-P axis, whereas that of C57BL/6J lies on the M-L axis).
Global differences: V1, S1, A1
Separation by strain is evident in the lateral and rostral V1 landmarks, and the most lateral S1 landmark. Other landmarks show no clear separation between strains.
Deformation grids and vector plots
Uniform components of shape
Despite the presence of significant uniform shape components in the local PMBSF, and the presence of significant shape differences across V1, S1, and A1, there were no significant uniform shape components in this second landmark configuration.
The 3 landmarks for the triangular shape of V1 were analyzed separately. Approximately 29.7% of variation in shape of V1 is explained by strain. Using a permutation test to confirm strain differences in shape variation, 1 in 10,000 permuted data sets had a Goodall's F statistic as great or greater than the F = 12.3117 (df = 2, 58) for the observed data, giving a significance probability of 0.0001. V1 appears elongated both in the anteroposterior and mediolateral dimensions in C57BL/6J relative to DBA/2J.
Cortical arealization initiates from early secreted ligands and graded distribution of transcription factors, which trigger the activation of distinctive area-specific morphogenetic programs, and result in an ordered topology of species-typical primary sensory and motor cortical areas (see [1–4] for review). Manipulation of transcription factor gradients have altered caudal and rostral cortical allotment to specific cortical fields. For example, Hamasaki et al.,  caused a caudal over expression of Emx2. Homzygote ne-Emx2 mice, compared to wild-type mice, showed a gradient of change in cortical area size from the caudal to the rostral poles. Transgenic ne-Emx2 mice showed a 52% increase in V1, a 25% decrease in PMBSF, and a 36% decrease in motor cortex. There was also a rostrolateral shift in the positions of S1 and A1. Hamasaki et al.,  have been the only group to quantitate in some detail the anatomical effects of their genetic manipulations, although their methods are unrelated to the geometric morphometric descriptions of shape in two dimensions used in the present paper.
The application of shape statistics to sections of the flattened cortical sheet is an obvious but powerful innovation that has the potential to improve and unify attempts to detect and describe effects of genetic manipulations on cortical arealization. In addition to improving phenotypic description of genetic manipulations putatively affecting cortex arealization, shape analysis of the cortical area map enables forward genetic efforts to isolate quantitative trait loci (QTLs). This approach would correlate measured quantitative phenotypic variation in cortex to genetic variation between mice or lines of mice (e.g., [21, 18]). Shape statistics has been applied successfully to isolate a large number of QTLs for mouse mandible shape . Cortical morphogenesis has been predominantly studied using experimental (reverse) genetic strategies. Relatively few genes have been identified that play a role in specifying the cortical area map, and our understanding of the genetic network active during cortical specification is poor. Genes that have been identified, however, are well supported [2, 4, 1, 3]. Forward genetic approaches to cortical morphogenesis would provide a complementary strategy. Although the conversion of QTLs to identified genes is difficult, advances in forward genetic methods have demonstrated success [23–25].
Another approach is to use gene microarray technologies to identify spatial or temporal differences in gene expression in the developing cortex. Use of microarrays in genetically manipulated mice , or mice of differing genetic backgrounds [26, 27], is a promising discovery approach to understanding the network of transcription factors that must coordinately define the cortical area map. Sansom et al., , used microarrays to characterize the early cortical transcriptome of mice lacking a receptor important in Fgf signaling, and identified and validated the novel target gene Mest. Given differences in V1 area and shape between C57BL/6J and DBA/2J strains, microarrays could be used to test for differences between strains in the transcriptome from caudal, embryonic cortex. Although the work of Sansom et al., , is a model study, it lacked direct correlation of the transcriptome differences to cortical map phenotypes.
An approach that integrates genetic, transcript, and phenotypic variation would use the genetic reference population of BXD recombinant inbred lines . This panel of genetically varied lines derived from C57BL/6J and DBA/2J strains can be phenotyped for cortical area and shape differences (as in  and this paper) and also characterized for anterior-posterior early cortical gene expression differences between lines (as in  or ). Since both parental inbred strains are fully sequenced, the full complement of genetic variation is archival in the BXD and need not be re-genotyped (see ). Phenotype to genotype (classical trait QTLs), phenotype to transcript, and transcript to genotype (expression QTLs) relationships could then all be evaluated, genome-wide . Because we have now shown that the parental strains C57BL/6J and DBA/2J markedly differ in both quantitative measures of cortex area size  and shape, this assures variation in the derivative BXD lines, and provides an empirical basis for using the BXD panel to study cortical development.
C57BL/6J and DBA/2J have markedly different cortical area maps, in both size and shape. These differences suggest polymorphism in genetic factors underlying cortical specification, even between common isogenic strains. Comparing cortical phenotypes between normally varying inbred mice or between genetically modified mice can identify genetic contributions to cortical specification. Geometric morphometrics represents an additional quantitative tool for the study of cortical development, regardless of whether it is studied from phenotype to gene or gene to phenotype.
Mice were purchased from Jackson Laboratories at 4–6 weeks of age, housed on a 12:12 light:dark cycle in same sex groups in standard laboratory animal cages (5 animals per cage). All experimental procedures were performed in accordance with the Guidelines for the Care and Use of Laboratory Animals published by the National Institutes of Health (publication 86–23) and the Vanderbilt University Animal Care and Use Committee. Mice were provided a standardized diet and clean water ad libitum. Mice at this age are not compromised by known visual and or auditory sensorineural deficits common to older animals of these strains. None of the mice used had visible body or facial wounds.
At 6–8 weeks of age (14 day age span), young adult mice were brought to complete anesthesia with a sodium pentobarbital overdose (100 mg/kg) injected intraperitoneally (IP), and then transcardially perfused with 0.1 M phosphate buffered 0.9% saline wash followed by 3% buffered paraformaldehyde fixative. Intact brains were removed from the skull, and the cortex was dissected free of the underlying white matter. Dissected cortices were flattened between glass slides and transferred to 30% sucrose for 12–18 hours. Cortices were sectioned parallel to the cortical surface at a thickness of 70–80 μm on a freezing, sliding microtome, stained for cytochrome oxidase according to the method of Wong-Riley , mounted on glass slides, air dried, and coverslipped. Thicker sections encourage all cortical areas to be within a single section, and are less prone to deformation or tearing.
A scale bar and the outlines of cortical regions of interest (ROIs) were drawn under a light microscope with a camera lucida attachment. Regions of interest included neocortex (C), visual cortex (V1), auditory cortex (A1), somatosensory cortex (S1), and select barrels of the posterior medial barrel subfield (PMBSF). Barrels drawn in rows A, B, C, D, and E were standardized to 5, 4, 6, 7, and 8 barrels, respectively (alpha, beta, gamma, and delta barrels were also included). Digital scans of the drawings were imported into a computer and landmark coordinates acquired with NIH ImageJ software after scale bar calibration . All data were acquired blind to animal and strain identity, and the order of acquisition was independent of animal or strain.
Two landmark configurations were generated for this study, the first to assess local shape differences in the PMBSF, the second to assess global shape differences across the cortical sheet (see Figure 2). Landmarks are points that can be located precisely (repeatedly) on the anatomical structure under study and display a one-to-one correspondence among the specimens included in the study. The first landmark configuration consisted of 34 barrel centroids. Centroids are the average coordinates of a set of coordinates. In ImageJ a centroid is the average X and Y from all pixel coordinates within a region of interest. Due to missing data in some animals for barrels alpha, beta, gamma, delta, A4–A5, B4, and D7, the final data set consisted of 26 landmarks for barrels A1–A3, B1–B3, C1–C6, D1–D6, and E1–E8, in 13 C57BL/6J mice (4 males, 9 females) and 12 DBA/2J mice (8 males, 4 females). The methods of shape analysis used in this paper required that all landmarks be present in all cortices measured. The second landmark configuration consisted of 7 landmarks, 3 that reduced V1 to a triangle, 3 that reduced S1 to a triangle, and 1 that was the centroid of A1. For this configuration, 16 C57BL/6J mice (6 males, 10 females) were compare to 15 DBA/2J mice (9 males, 6 females). In this paper sex effects in the cortical area maps were not a focus, and sex was determined not to be a significant factor by Goodall's F tests performed separately on each strain for each landmark configuration (each P > 0.05). Li et al., , also recently looked at the barrels of C57BL/6J, DBA/2J, and a small sample of derivative BXD lines, and found no sex differences in PMBSF area, although shape was not examined. Li et al., also found no differences in PMBSF area over a 21 day age spanning the same ages sampled in this study. In our data, over a 14 day age span, we found no significant difference that we could attribute to age, in either strain for either landmark configuration (Goodall's F tests, each P > 0.05). No substantive differences in equipment or procedures confounded the significant strain differences we report in this paper.
Reliability of landmark configurations was assessed by acquiring (drawing) each landmark configuration in triplicate for each animal. The first set of drawings was discarded to avoid training effects. The second and third sets were tested for between set differences, using shape analysis. None were found for either landmark configuration (Goodall's F tests, P > 0.05). The final data set used for shape analysis was formed by averaging over the second and third drawing sets, as well as hemisphere within animal, if two hemispheres were drawn. If two hemispheres were drawn, they were appropriately flipped to the same orientation. In this paper we do not consider average differences in left and right hemispheric cortical maps. Studies in rats  and mice  found no left and right hemispheric differences in PMBSF area, although shape was not examined. Before shape analysis, final data sets were preanalyzed by the program tpsSmall  to confirm variation was appropriate (small enough) for shape analysis.
The collection of complex multivariate statistical methods used to analyze shape are organized around Kendall's definition of shape . Shape is defined as all the geometric information that remains when location, scale and rotational effects are filtered out from an object or landmark configuration. Shape analysis in this study made use of generalized Procrustes analysis (GPA) to superimpose landmark configurations and remove nonshape variation in the landmark coordinates [34, 20, 35]. The algorithm of this method involves three transformations [7, 9]. First, the centroid of each configuration is translated to the origin by subtracting centroid coordinates from each landmark. Second, the configuration is then scaled by dividing each landmark coordinate by the centroid size of that configuration. The centroid size is the square root of the summed squared distances of each landmark from the centroid of the landmark configuration. Third, with respect to a given configuration, another configuration is rotated to minimize the summed squared distances between labeled landmarks. In tpsRegr, this is done iteratively: one reference configuration is chosen and every other configuration is fitted to it, but from the second round on, the average of the coordinates in the previous fit is used as a reference. Final superimposition is not influenced by which landmark configuration is chosen to begin the procedure.
A statistical test adapted to the coordinates produced by Procrustes fit is Goodall's F test , which compares the difference in mean shape between two samples relative to the shape variation found within the samples. To determine statistical significance we employed a permutation test based on Goodall's F statistic. In this test, the data are permuted (randomized) and Goodall's F statistic is calculated. This is done 10,000 times. The proportion of Goodall's F statistics from permuted data sets as great or greater than the Goodall's statistic on the original data set is given as the significance probability. Use of permutation relaxes some of the restrictive assumptions of Goodall's F test. Goodall's F test only considers the total amount of shape variation, and does not consider the directionality of the variation. With small samples (relative to the number of landmark coordinates), as in this study, this is a useful property. Indeed, the sample size in this study is too small to use alternative MANOVA multivariate statistics to test for shape differences in the PMBSF, given the number of barrel landmarks and the number of mice measured.
In addition to scatter plots of Procrustes transformed landmark data, broken down by symbol to represent strain (e.g., Figure 3 above), thin-plate spline deformation grid plots are used (e.g., Figure 4 above). A thin-plate spline deformation grid represents a smooth interpolation, mapping the coordinates of one landmark configuration into another. These plots are analogous to an elastic sheet of graph paper that is stretched to map one morph to another, and are in the spirit of the historically famous efforts of D'Arcy W. Thompson . Deformation grids are useful for visualizing changes between landmarks over an entire configuration. Aside from deformation grids, vector plots are also provided. These represent the change in landmark positions by a vector with size proportional to the observed differences.
Although the statistical methods used in this paper are complex, reliable software for geometric morphometrics is readily and freely available and makes these methods accessible to neuroscientists. The primary software used in this study was tpsRegr by statistician F. James Rohlf . The strain differences reported above were confirmed using two additional software titles that implement Procrustes fit and Goodall's F test: TwoGroup by H. D. Sheets , and shapes, an R package by Ian Dryden .
We thank Dr. Jon Kaas, Dr. Elaine Sanders-Bush, and Dr. Randy Blakely for encouragement, lab space, and materials. We thank Jon Kaas, Ford Ebner, and Erin Henry for manuscript comments. We thank F. James Rohlf and Chris P. Klingenberg for guiding a revision. We thank Dennis O'Leary for permission to adapt and use figures from . We thank Elsevier for permission to use figures from . This work was supported by NIH NS16446 to Jon Kaas and NIH MH34007 to Elaine Sanders-Bush, and a VUMC NIH training grant (T32 MH065215).
- O'Leary D, Nakagawa Y: Patterning centers, regulatory genes and extrinsic mechanisms controlling arealization of the neocortex. Curr Opin Neurobiol. 2002, 12: 14-25. 10.1016/S0959-4388(02)00285-4.View ArticlePubMedGoogle Scholar
- Grove E, Fukuchi-Shimogori T: Generating the cerebral cortical area map. Annu Rev Neurosci. 2003, 26: 355-380. 10.1146/annurev.neuro.26.041002.131137.View ArticlePubMedGoogle Scholar
- Rash B, Grove E: Area and layer patterning in the developing cerebral cortex. Curr Opin Neurobiol. 2006, 16: 25-34. 10.1016/j.conb.2006.01.004.View ArticlePubMedGoogle Scholar
- Mallamaci A, Stoykova A: Gene networks controlling early cerebral cortex arealization. Eur J Neurosci. 2006, 23 (4): 847-856. 10.1111/j.1460-9568.2006.04634.x.View ArticlePubMedGoogle Scholar
- Hamasaki T, Leingartner A, Ringstedt T, O'Leary D: EMX2 regulates sizes and positioning of the primary sensory and motor areas in neocortex by direct specification of cortical progenitors. Neuron. 2004, 43 (3): 359-372. 10.1016/j.neuron.2004.07.016.View ArticlePubMedGoogle Scholar
- Airey D, Robbins A, Enzinger K, Wu F, Collins C: Variation in the cortical area map of C57BL/6J and DBA/2J inbred mice predicts strain identity. BMC Neurosci. 2005, 6: 18-10.1186/1471-2202-6-18.PubMed CentralView ArticlePubMedGoogle Scholar
- Dryden IL, Mardia KV: Statistical shape analysis Volume Wiley series in probability and statistics Probability and statistics. 1998, Chichester; New York: John Wiley & SonsGoogle Scholar
- Klingenberg C: Morphometrics and the role of the phenotype in studies of the evolution of developmental mechanisms. Gene. 2002, 287 (1–2): 3-10. 10.1016/S0378-1119(01)00867-8.View ArticlePubMedGoogle Scholar
- Zelditch ML, Swiderski DL, Sheets DH, Fink WL: Geometric Morphometrics for Biologists. 2004, Academic PressGoogle Scholar
- Gharaibeh W, Rohlf F, Slice D, DeLisi L: A geometric morphometric assessment of change in midline brain structural shape following a first episode of schizophrenia. Biol Psychiatry. 2000, 48 (5): 398-405. 10.1016/S0006-3223(00)00916-1.View ArticlePubMedGoogle Scholar
- Bookstein F, Sampson P, Streissguth A, Connor P: Geometric morphometrics of corpus callosum and subcortical structures in the fetal-alcohol-affected brain. Teratology. 2001, 64: 4-32. 10.1002/tera.1044.View ArticlePubMedGoogle Scholar
- Styner M, Gerig G, Lieberman J, Jones D, Weinberger D: Statistical shape analysis of neuroanatomical structures based on medial models. Med Image Anal. 2003, 7 (3): 207-220. 10.1016/S1361-8415(02)00110-X.View ArticlePubMedGoogle Scholar
- Davatzikos C: Voxel-based morphometric analysis using shape transformations. Int Rev Neurobiol. 2005, 66: 125-146.View ArticlePubMedGoogle Scholar
- Li C, Wei X, Lu L, Peirce J, Williams R, Waters R: Genetic analysis of barrel field size in the first somatosensory area (SI) in inbred and recombinant inbred strains of mice. Somatosens Mot Res. 2005, 22 (3): 141-150. 10.1080/08990220500262182.View ArticlePubMedGoogle Scholar
- Riddle D, Purves D: Individual variation and lateral asymmetry of the rat primary somatosensory cortex. J Neurosci. 1995, 15 (6): 4184-4195.PubMedGoogle Scholar
- Xu Y, Sari Y, Zhou F: Selective serotonin reuptake inhibitor disrupts organization of thalamocortical somatosensory barrels during development. Brain Res Dev Brain Res. 2004, 150 (2): 151-161. 10.1016/j.devbrainres.2003.02.001.View ArticlePubMedGoogle Scholar
- Van der Loos H, Dorfl J, Welker E: Variation in pattern of mystacial vibrissae in mice. A quantitative study of ICR stock and several inbred strains. J Hered. 1984, 75 (5): 326-336.PubMedGoogle Scholar
- Airey D, Lu L, Shou S, Williams R: Genetic sources of individual differences in the cerebellum. Cerebellum. 2002, 1 (4): 233-240. 10.1080/147342202320883542.View ArticlePubMedGoogle Scholar
- Chesler E, Lu L, Shou S, Qu Y, Gu J, Wang J, Hsu H, Mountz J, Baldwin N, Langston M, Threadgill D, Manly K, Williams R: Complex trait analysis of gene expression uncovers polygenic and pleiotropic networks that modulate nervous system function. Nat Genet. 2005, 37 (3): 233-242. 10.1038/ng1518.View ArticlePubMedGoogle Scholar
- Rohlf F, Bookstein F: Computing the uniform component of shape variation. Syst Biol. 2003, 52: 66-69.View ArticlePubMedGoogle Scholar
- Airey D, Lu L, Williams R: Genetic control of the mouse cerebellum: identification of quantitative trait loci modulating size and architecture. J Neurosci. 2001, 21 (14): 5099-5109.PubMedGoogle Scholar
- Klingenberg C, Leamy L, Cheverud J: Integration and modularity of quantitative trait locus effects on geometric shape in the mouse mandible. Genetics. 2004, 166 (4): 1909-1921. 10.1534/genetics.166.4.1909.PubMed CentralView ArticlePubMedGoogle Scholar
- Flint J, Valdar W, Shifman S, Mott R: Strategies for mapping and cloning quantitative trait genes in rodents. Nat Rev Genet. 2005, 6 (4): 271-286. 10.1038/nrg1576.View ArticlePubMedGoogle Scholar
- Leygraf A, Hohoff C, Freitag C, Willis-Owen S, Krakowitzky P, Fritze J, Franke P, Bandelow B, Fimmers R, Flint J, Deckert J: Rgs 2 gene polymorphisms as modulators of anxiety in humans?. J Neural Transm. 2006Google Scholar
- Yalcin B, Willis-Owen S, Fullerton J, Meesaq A, Deacon R, Rawlins J, Copley R, Morris A, Flint J, Mott R: Genetic dissection of a behavioral quantitative trait locus shows that Rgs2 modulates anxiety in mice. Nat Genet. 2004, 36 (11): 1197-1202. 10.1038/ng1450.View ArticlePubMedGoogle Scholar
- Sansom S, Hebert J, Thammongkol U, Smith J, Nisbet G, Surani M, McConnell S, Livesey F: Genomic characterisation of a Fgf-regulated gradient-based neocortical protomap. Development. 2005, 132 (17): 3947-3961. 10.1242/dev.01968.View ArticlePubMedGoogle Scholar
- Funatsu N, Inoue T, Nakamura S: Gene expression analysis of the late embryonic mouse cerebral cortex using DNA microarray: identification of several region- and layer-specific genes. Cereb Cortex. 2004, 14 (9): 1031-1044. 10.1093/cercor/bhh063.View ArticlePubMedGoogle Scholar
- Peirce J, Lu L, Gu J, Silver L, Williams R: A new set of BXD recombinant inbred lines from advanced intercross populations in mice. BMC Genet. 2004, 5: 7-10.1186/1471-2156-5-7.PubMed CentralView ArticlePubMedGoogle Scholar
- WebQTL. [http://www.genenetwork.org]
- Wong-Riley M: Changes in the visual system of monocularly sutured or enucleated cats demonstrable with cytochrome oxidase histochemistry. Brain Res. 1979, 171: 11-28. 10.1016/0006-8993(79)90728-5.View ArticlePubMedGoogle Scholar
- Image Processing and Analysis in Java. [http://rsb.info.nih.gov/ij/]
- SUNY Stoneybrook Morphometrics. [http://life.bio.sunysb.edu/morph/index.html]
- Kendall D: The diffusion of shape. Advances in Applied Probability. 1977, 9: 428-430. 10.2307/1426091.View ArticleGoogle Scholar
- Rohlf F, Slice D: Extensions of the Procrustes method for the optimal superimposition of landmarks. Systematic Zoology. 1990, 39: 40-59. 10.2307/2992207.View ArticleGoogle Scholar
- Rohlf F: Bias and error in estimates of mean shape in geometric morphometrics. J Hum Evol. 2003, 44 (6): 665-683. 10.1016/S0047-2484(03)00047-2.View ArticlePubMedGoogle Scholar
- Goodall C: Procrustes methods in the statistical analysis of shape. Journal of the Royal Statistical Society. 1991, 53 (2): 285-339.Google Scholar
- Thompson DW: On growth and form. 1917, Cambridge [Eng.]: University pressView ArticleGoogle Scholar