Hemispheric asymmetries in cortical grey matter of gyri and sulci in modern human populations from South America.

Structural asymmetries of brain regions associated with lateralised functions have been extensively studied. However, there are fewer morphometric analyses of asymmetries of the gyri and sulci of the entire cortex. The current study assessed cortical asymmetries in a sample of healthy adults (N = 175) from an admixed population from South America. Grey matter volume and surface area of 66 gyri and sulci were quantified on T1 magnetic resonance images. The departure from zero of the differences between left and right hemispheres (L-R), a measure of directional asymmetry (DA), the variance of L-R, and an index of fluctuating asymmetry (FA) were evaluated for each region. Significant departures from perfect symmetry were found for most cortical gyri and sulci. Regions showed leftward asymmetry at the population level in the frontal lobe and superior lateral parts of the parietal lobe. Rightward asymmetry was found in the inferior parietal, occipital, frontopolar, and orbital regions, and the cingulate (anterior, middle, and posterior-ventral). Despite this general pattern, several sulci showed the opposite DA compared to the neighbouring gyri, which remarks the need to consider the neurobiological differences in gyral and sulcal development in the study of structural asymmetries. The results also confirm the absence of DA in most parts of the inferior frontal gyrus and the precentral region. This study contributes with data on populations underrepresented in the databases used in neurosciences. Among its findings, there is agreement with previous results obtained in populations of different ancestry and some discrepancies in the middle frontal and medial parietal regions. A significant DA not reported previously was found for the volume of long and short insular gyri and the central sulcus of the insula, frontomarginal, transverse frontopolar, paracentral, and middle and posterior parts of the cingulate gyrus and sulcus, gyrus rectus, occipital pole, and olfactory sulcus, as well as for the volume and area of the transverse collateral sulcus and suborbital sulcus. Also, several parcels displayed significant variability in the left-right differences, which can be partially attributable to developmental instability, a source of FA. Moreover, a few gyri and sulci displayed ideal FA with non-significant departures from perfect symmetry, such as subcentral and posterior cingulate gyri and sulci, inferior frontal and fusiform gyri, and the calcarine, transverse collateral, precentral, and orbital sulci. Overall, these results show that asymmetries are ubiquitous in the cerebral cortex.

Classifying high-dimensional phenotypes with ensemble learning.

Classification is a fundamental task in biology used to assign members to a class. While linear discriminant functions have long been effective, advances in phenotypic data collection are yielding increasingly high-dimensional datasets with more classes, unequal class covariances, and non-linear distributions. Numerous studies have deployed machine learning techniques to classify such distributions, but they are often restricted to a particular organism, a limited set of algorithms, and/or a specific classification task. In addition, the utility of ensemble learning or the strategic combination of models has not been fully explored.We performed a meta-analysis of 33 algorithms across 20 datasets containing over 20,000 high-dimensional shape phenotypes using an ensemble learning framework. Both binary (e.g., sex, environment) and multi-class (e.g., species, genotype, population) classification tasks were considered. The ensemble workflow contains functions for preprocessing, training individual learners and ensembles, and model evaluation. We evaluated algorithm performance within and among datasets. Furthermore, we quantified the extent to which various dataset and phenotypic properties impact performance.We found that discriminant analysis variants and neural networks were the most accurate base learners on average. However, their performance varied substantially between datasets. Ensemble models achieved the highest performance on average, both within and among datasets, increasing average accuracy by up to 3% over the top base learner. Higher class R2 values, mean class shape distances, and between- vs. within-class variances were positively associated with performance, whereas higher class covariance distances were negatively associated. Class balance and total sample size were not predictive.Learning-based classification is a complex task driven by many hyperparameters. We demonstrate that selecting and optimizing an algorithm based on the results of another study is a flawed strategy. Ensemble models instead offer a flexible approach that is data agnostic and exceptionally accurate. By assessing the impact of various dataset and phenotypic properties on classification performance, we also offer potential explanations for variation in performance. Researchers interested in maximizing performance stand to benefit from the simplicity and effectiveness of our approach made accessible via the R package pheble.

MusMorph, a database of standardized mouse morphology data for morphometric meta-analyses.

Complex morphological traits are the product of many genes with transient or lasting developmental effects that interact in anatomical context. Mouse models are a key resource for disentangling such effects, because they offer myriad tools for manipulating the genome in a controlled environment. Unfortunately, phenotypic data are often obtained using laboratory-specific protocols, resulting in self-contained datasets that are difficult to relate to one another for larger scale analyses. To enable meta-analyses of morphological variation, particularly in the craniofacial complex and brain, we created MusMorph, a database of standardized mouse morphology data spanning numerous genotypes and developmental stages, including E10.5, E11.5, E14.5, E15.5, E18.5, and adulthood. To standardize data collection, we implemented an atlas-based phenotyping pipeline that combines techniques from image registration, deep learning, and morphometrics. Alongside stage-specific atlases, we provide aligned micro-computed tomography images, dense anatomical landmarks, and segmentations (if available) for each specimen (N = 10,056). Our workflow is open-source to encourage transparency and reproducible data collection. The MusMorph data and scripts are available on FaceBase ( www.facebase.org , https://doi.org/10.25550/3-HXMC ) and GitHub ( https://github.com/jaydevine/MusMorph ).

Morphological correspondence between brain and endocranial surfaces in mice exposed to undernutrition during development.

The morphological changes of the brain and the skull are highly integrated as a result of shared developmental pathways and different types of interactions between them. Shared developmental trajectories between these two structures might be influenced by genetic and environmental factors. Although the effect of environmental factors on neural and craniofacial traits has been extensively studied, less is known about the specific impact of stressful conditions on the coordinated variation between these structures. Here, we test the effect of early nutrient restriction on morphological correspondence between the brain and the endocast. For this purpose, mice exposed to protein or calorie-protein restriction during gestation and lactation were compared with a control group in which dams were fed standard food ad libitum. High-resolution images were obtained after weaning to describe brain and endocranial morphology. By magnetic resonance imaging (MRI), brain volumes were obtained and endocasts were segmented from skull reconstructions derived from micro-computed tomography (microCT). Brain and endocranial volumes were compared to assess the correspondence in size. Shape changes were analyzed using a set of landmarks and semilandmarks on 3D surfaces. Results indicated that brain volume is relatively less affected by undernutrition during development than endocast volume. Shape covariation between the brain and the endocast was found to be quite singular for protein-restricted animals. Procrustes distances were larger between the brain and the endocast of the same specimens than between brains or endocasts of different animals, which means that the greatest similarity is by type of structure and suggests that the use of the endocast as a direct proxy of the brain at this intraspecific scale could have some limitations. In the same line, patterns of brain shape asymmetry were not directly estimated from endocranial surfaces. In sum, our findings indicate that morphological variation and association between the brain and the endocast is modulated by environmental factors and support the idea that head morphogenesis results from complex processes that are sensitive to the pervasive influence of nutrient intake.

Influence of accessory sulci of the frontoparietal operculum on gray matter quantification.

Introduction: The perisylvian region is the cortical core of language and speech. Several accessory sulci have been described in this area, whose presence could modify the results of the automatic quantification of gray matter by popularly used software. This study aimed to assess the expression of accessory sulci in the frontoparietal operculum (FPO) and to evaluate their influence on the gray matter volume estimated by an automatic parcellation of cortical gyri and sulci. Methods: Brain MRI scans of 100 healthy adult volunteers were visually analyzed. The existence of the triangular and diagonal sulci, and the number of accessory sulci in the frontoparietal operculum, were assessed on T1 images. Also, the gray matter volume of gyri and sulci was quantified by an automatized parcellation method. Interhemispheric differences in accessory sulci were evaluated with Chi-square and Wilcoxon paired tests. The effects of the hemisphere, sex, age, total intracranial volume, and accessory sulci on morphometric variables were assessed by linear models. Results: These sulci were found in more than half of the subjects, mostly in the left hemisphere, and showed a significant effect on the gray matter content of the FPO. In particular, the volume of the inferior frontal sulcus, pars opercularis of the inferior frontal gyrus, horizontal ramus of the lateral sulcus, angular gyrus, and postcentral gyrus showed a significant influence on the presence of accessory sulci. Discussion: The prevalence of tertiary sulci in the FPO is high, although their meaning is not yet known. Therefore, they should be considered to reduce the risk of misclassifications of normal variation.

Association between shape changes and bone remodeling patterns in the middle face during ontogeny in South American populations.

The morphology of facial bones is modeled by processes of bone formation and resorption induced by interactions between tissues and compensatory responses. However, the role of remodeling patterns on the morphological changes within and among populations has been scarcely explored. Here, we assess the association between facial shape and the underlying bone cell activity throughout the ontogeny in two Amerindian populations that represent the extremes of craniofacial variation in South America. The sample comprises 71 individuals (36 adults and 35 subadults) representing hunter-gatherers from Patagonia and horticulturists from Northwest Argentina. We analyzed the shape and size of the zygomatic and the maxilla, and compared them with the patterns of bone formation and resorption. Bone formation and resorption were described by quantitative histological analysis of bone surfaces. Morphological changes were described by landmarks and semilandmarks digitized on 3D surfaces obtained from CT images. The results from multivariate statistics analysis show that the patterns of bone remodeling are associated with variation in the morphology of the middle face. We found a similar pattern of facial shape variation along the ontogenetic trajectory in the two samples, and a similar trend in the amount of formation and resorption activities across ages. The main differences between samples were found in the distribution of the areas of bone formation and resorption, possibly associated with mechanical bone response to masticatory loading. These findings provide clues about the processes and mechanisms of bone development involved in the facial morphological differentiation in human populations from southern South America.

Endocranial asymmetry in New World monkeys: a comparative phylogenetic analysis of morphometric data.

Brain lateralization is a widespread phenomenon although its expression across primates is still controversial due to the reduced number of species analyzed and the disparity of methods used. To gain insight into the diversification of neuroanatomical asymmetries in non-human primates we analyze the endocasts, as a proxy of external brain morphology, of a large sample of New World monkeys and test the effect of brain size, home range and group sizes in the pattern and magnitude of shape asymmetry. Digital endocasts from 26 species were obtained from MicroCT scans and a set of 3D coordinates was digitized on endocast surfaces. Results indicate that Ateles, Brachyteles, Callicebus and Cacajao tend to have a rightward frontal and a leftward occipital lobe asymmetry, whereas Aotus, Callitrichinae and Cebinae have either the opposite pattern or no directional asymmetry. Such differences in the pattern of asymmetry were associated with group and home range sizes. Conversely, its magnitude was significantly associated with brain size, with larger-brained species showing higher inter-hemispheric differences. These findings support the hypothesis that reduction in inter-hemispheric connectivity in larger brains favors the lateralization and increases the structural asymmetries, whereas the patterns of shape asymmetry might be driven by socio-ecological differences among species.

Lack of head sparing following third-trimester caloric restriction among Tanzanian Maasai.

The reduction of food intake during pregnancy is part of many cultural and religious traditions around the world. The impact of such practices on fetal growth and development are poorly understood. Here, we examined the patterns of diet intake among Maasai pregnant women and assessed their effect on newborn morphometrics. We recruited 141 mother-infant pairs from Ngorongoro Conservation Area (NCA) in Northern Tanzania and quantified dietary intake and changes in maternal diet during pregnancy. We obtained measurements of body weight (BW) and head circumference (HC) at birth. We found that Maasai women significantly reduced their dietary intake during the third trimester, going from an average of 1601 kcal/day during the first two trimesters to 799 kcal/day in the final trimester. The greatest proportion of nutrient reduction was in carbohydrates. Overall, 40% of HC Z-scores of the NCA sample were more than 2 standard deviations below the WHO standard. Nearly a third of neonates classify as low birth weight (< 2500g). HC was smaller relative to BW in this cohort than predicted using the WHO standard. This contrasts markedly to a Tanzanian birth cohort obtained at the same time in an urban context in which only 12% of infants exhibited low weight, only two individuals had HC Z-scores < 2 and HC's relative to birth weight were larger than predicted using the WHO standards. The surprising lack of head sparing in the NCA cohort suggests that the impact of third trimester malnutrition bears further investigation in both animal models and human populations, especially as low HC is negatively associated with long term health outcomes.

Modeling the effect of brain growth on cranial bones using finite-element analysis and geometric morphometrics.

PURPOSE: Brain expansion during ontogeny has been identified as a key factor for explaining the growth pattern of neurocranial bones. However, the dynamics of this relation are only partially understood and a detailed characterization of integrated morphological changes of the brain and the neurocranium along ontogeny is still lacking. The aim of this study was to model the effect of brain growth on cranial bones by means of finite-element analysis (FEA) and geometric morphometric techniques. METHODS: First, we described the postnatal changes in brain size and shape by digitizing coordinates of 3D semilandmarks on cranial endocasts, as a proxy of brain, segmented from CT-scans of an ontogenetic sample. Then, two scenarios of brain growth were simulated: one in which brain volume increases with the same magnitude in all directions, and other that includes the information on the relative expansion of brain regions obtained from morphometric analysis. RESULTS: Results indicate that in the first model, in which a uniform pressure is applied, the largest displacements were localized in the sutures, especially in the anterior and posterior fontanels, as well as the metopic suture. When information of brain relative growth was introduced into the model, displacements were also concentrated in the lambda region although the values along both sides of the neurocranium (parietal and temporal bones) were larger than under the first scenario. CONCLUSION: In sum, we propose a realistic approach to the use of FEA based on morphometric data that offered different results to more simplified models.

Region-specific changes in Mus musculus brain size and cell composition under chronic nutrient restriction.

Nutrition is one of the most influential environmental factors affecting the development of different tissues and organs. It is suggested that under nutrient restriction the growth of the brain is spared as a result of the differential allocation of resources from other organs. However, it is not clear whether this sparing occurs brain-wide. Here, we analyzed morphological changes and cell composition in different regions of the offspring mouse brain after maternal exposure to nutrient restriction during pregnancy and lactation. Using high-resolution magnetic resonance imaging, we found that brain regions were differentially sensitive to maternal protein restriction and exhibited particular patterns of volume reduction. The cerebellum was reduced in absolute and relative volume, while cortex volume was relatively preserved. Alterations in cell composition (examined by the isotropic fractionator method) and organization of white matter (measured by diffusor tensor images) were also region specific. These changes were not related to the metabolic rate of the regions and were only partially explained by their specific growth trajectories. This study is a first step towards understanding the mechanisms of regional brain sparing at microstructural and macrostructural levels resulting from undernutrition.

Integration and the Developmental Genetics of Allometry.

Allometry refers to the ways in which organismal shape is associated with size. It is a special case of integration, or the tendency for traits to covary, in that variation in size is ubiquitous and evolutionarily important. Allometric variation is so commonly observed that it is routinely removed from morphometric analyses or invoked as an explanation for evolutionary change. In this case, familiarity is mistaken for understanding because rarely do we know the mechanisms by which shape correlates with size or understand their significance. As with other forms of integration, allometric variation is generated by variation in developmental processes that affect multiple traits, resulting in patterns of covariation. Given this perspective, we can dissect the genetic and developmental determinants of allometric variation. Our work on the developmental and genetic basis for allometric variation in craniofacial shape in mice and humans has revealed that allometric variation is highly polygenic. Different measures of size are associated with distinct but overlapping patterns of allometric variation. These patterns converge in part on a common genetic basis. Finally, environmental modulation of size often generates variation along allometric trajectories, but the timing of genetic and environmental perturbations can produce deviations from allometric patterns when traits are differentially sensitive over developmental time. These results question the validity of viewing allometry as a singular phenomenon distinct from morphological integration more generally.

Variation in facial bone growth remodeling in prehistoric populations from southern South America.

OBJECTIVES: To assess the intraspecific variation in bone remodeling patterns in modern humans, we studied two populations from southern South America that represent the extremes of morphological variation in this region. We particularly analyzed the ontogenetic changes in the patterns of bone growth remodeling and compared the patterns between samples. MATERIALS AND METHODS: We obtained high-resolution casts of the periosteal surface of the upper and middle face of subadults (n = 36) and adult (n = 36) individuals from a sample of hunter-gatherers from Patagonia and a sample of horticulturists from Northwest Argentina. The areas of bone formation and resorption were registered using an incident-light microscope. We then estimated the average bone remodeling map by sample and age, and performed principal component analysis and multivariate regressions to assess the extension and distribution of these areas across ontogeny and between samples. RESULTS: We found that the remodeling pattern of the glabella, supraorbital arch, frontal process of the maxilla, and a large part of the zygomatic bone is relatively constant in subadults and adults from both sample with a clear predominance of bone formation. In contrast, the middle face is characterized by the spatial alternation between formation and resorption areas, and greater variation with age and between samples. The main differences were found in areas related to chewing and muscle insertions. CONCLUSIONS: Our study provides the first evidence of interpopulation variation in bone growth remodeling and suggests that biomechanical factors can influence the observed patterns. It also underlines the need to account for ecological factors in within and between species comparisons.

Mitochondrial DNA Diversity and Evolutionary History of Native Human Populations of Argentinean Northwest Patagonia.

The genetic composition of Amerindian descendants from Patagonia has long been a focus of interest, although the information available is still scarce for many geographic areas. Here, we report the first analysis of the variation in the mitochondrial DNA (mtDNA) control region for an area of northwestern Patagonia, the North of Neuquén, with the aim of studying the processes and historical events that modeled the evolutionary history of these human groups. We analyzed 113 individuals from two localities of northern Neuquén, along with 6 from southern Neuquén and 223 previously published mtDNA sequences from neighboring areas in Argentina and Chile. We estimated the haplotypic variation and spatial structure of molecular variability. Amerindian subhaplogroups predominate in the two samples from northern Neuquén (n = 70), with D1g and C1b13 the most represented, although in different proportions. These samples exhibit Amerindian mtDNA haplotypes similar to the variants from neighboring areas. Most of haplotype variability was within group; variation among groups was relatively low and scarcely associated with geographical space. The most frequent subhaplogroups in northern Neuquén are characteristic of native populations from Patagonia and Chilean Araucanía, and probably originated in the region during the Late Pleistocene or Early Holocene. However, the spatial variation of mtDNA haplotypes departs from a latitudinal pattern and suggests differential levels of gene flow among areas during the Late Holocene, with moderate levels across the North of Neuquén as well as between this area and neighboring populations from Chile, the South of Neuquén, and Río Negro.

Integrating data on bone modeling and morphological ontogenetic changes of the maxilla in modern humans.

The aim of this work is to assess the association between the patterns of bone modeling and the changes in shape and size of the maxilla along ontogeny in modern humans. The sample analyzed includes an ontogenetic series of 30 individuals from an archeological site from Pampa Grande, northwest of Argentina. The areas of bone resorption and formation were described by histological analysis of bone surfaces and then quantified using spatial statistics. Morphological changes were analyzed by geometric morphometric methods using landmarks and semilandmarks digitized on 3D surfaces obtained from CT-scans. The regression of bone modeling maps on the centroid size shows no significant association between both variables neither in subadult nor adult individuals. On the contrary, the results of the partial least squares analysis shows a strong association between the shape changes in the maxilla with changes in the pattern of bone modeling in both groups of age, subadults and adults. Overall, this study contributes to the understanding of the mechanisms and processes that model maxillary morphology during growth.

Brain Structural Networks in Mouse Exposed to Chronic Maternal Undernutrition.

Brain structural connectivity is known to be altered in cases of intrauterine growth restriction and premature birth, although the specific effect of maternal nutritional restriction, a common burden in human populations, has not been assessed yet. Here we analyze the effects of maternal undernutrition during pregnancy and lactation by establishing three experimental groups of female mice divided according to their diet: control (Co), moderate calorie-protein restriction (MCP) and severe protein restriction (SP). Nutritionally restricted dams gained relatively less weight during pregnancy and the body weight of the offspring was also affected by maternal undernutrition, showing global growth restriction. We performed magnetic resonance imaging (MRI) of the offspring's brains after weaning and analyzed their connectivity patterns using complex graph theory. In general, changes observed in the MCP group were more subtle than in SP. Results indicated that brain structures were not homogeneously affected by early nutritional stress. In particular, the growth of central brain regions, such as the temporo-parietal cortex, and long integrative myelinated tracts were relatively preserved, while the frequency of short tracts was relatively reduced. We also found a differential effect on network parameters: network degree, clustering, characteristic path length and small-worldness remained mainly unchanged, while the rich-club index was lower in nutritionally restricted animals. Rich-club decrease reflects an impairment in the structure by which brain regions with large number of connections tend to be more densely linked among themselves. Overall, the findings presented here support the hypothesis that chronic nutritional stress produces long-term changes in brain structural connectivity.

Body size and allometric variation in facial shape in children.

OBJECTIVES: Morphological integration, or the tendency for covariation, is commonly seen in complex traits such as the human face. The effects of growth on shape, or allometry, represent a ubiquitous but poorly understood axis of integration. We address the question of to what extent age and measures of size converge on a single pattern of allometry for human facial shape. METHODS: Our study is based on two large cross-sectional cohorts of children, one from Tanzania and the other from the United States (N = 7,173). We employ 3D facial imaging and geometric morphometrics to relate facial shape to age and anthropometric measures. RESULTS: The two populations differ significantly in facial shape, but the magnitude of this difference is small relative to the variation within each group. Allometric variation for facial shape is similar in both populations, representing a small but significant proportion of total variation in facial shape. Different measures of size are associated with overlapping but statistically distinct aspects of shape variation. Only half of the size-related variation in facial shape can be explained by the first principal component of four size measures and age while the remainder associates distinctly with individual measures. CONCLUSIONS: Allometric variation in the human face is complex and should not be regarded as a singular effect. This finding has important implications for how size is treated in studies of human facial shape and for the developmental basis for allometric variation more generally.

Facial shape manifestations of growth faltering in Tanzanian children.

Variation in the shape of the human face and in stature is determined by complex interactions between genetic and environmental influences. One such environmental influence is malnourishment, which can result in growth faltering, usually diagnosed by means of comparing an individual's stature with a set of age-appropriate standards. These standards for stature, however, are typically ascertained in groups where people are at low risk for growth faltering. Moreover, genetic differences among populations with respect to stature are well established, further complicating the generalizability of stature-based diagnostic tools. In a large sample of children aged 5-19 years, we obtained high-resolution genomic data, anthropometric measures and 3D facial images from individuals within and around the city of Mwanza, Tanzania. With genome-wide complex trait analysis, we partitioned genetic and environmental variance for growth outcomes and facial shape. We found that children with growth faltering have faces that look like those of older and taller children, in a direction opposite to the expected allometric trajectory, and in ways predicted by the environmental portion of covariance at the community and individual levels. The environmental variance for facial shape varied subtly but significantly among communities, whereas genetic differences were minimal. These results reveal that facial shape preserves information about exposure to undernourishment, with important implications for refining assessments of nutritional status in children and the developmental-genetics of craniofacial variation alike.

A quantitative approach for analysing bone modelling patterns from craniofacial surfaces in hominins.

Bone size and shape arise throughout ontogeny as a result of the coordinated activity of osteoblasts and osteoclasts, responsible for bone deposition and resorption, and growth displacements. The modelling processes leave specific microstructural features on the bone surface, which can be used to infer the mechanisms shaping craniofacial traits in extinct and extant species. However, the analysis of bone surfaces from fossils and archaeological samples faces some difficulties related to the bone loss caused by taphonomic factors, and the lack of formal methods for estimating missing information and comparing the patterns of bone modelling among several specimens and samples. The present study provides a new approach for the quantitative analysis of bone formation and resorption patterns obtained from craniofacial surfaces. First, interpolation techniques were used to estimate missing data on high-resolution replicas of the left maxilla in a sample of sub-adult and adult modern humans and sub-adult fossil hominins. The performance of this approach was assessed by simulating variable amounts of missing data. Then, we applied measures of dispersion and central tendency to represent the variation and average pattern of bone modelling within samples. The spatial interpolation resulted in reliable estimations of the type of cell activity (deposition or resorption) in the missing areas, even when large extensions of the bone surface were lost. The quantification of the histological data allowed us to integrate the information of different specimens and depict the areas with higher and lower variation in the bone modelling pattern of the maxilla among specimens. Overall, the main advantages of the quantitative approach used here for generating bone modelling patterns are the high replicability and the possibility of incorporating variation among specimens into the comparisons among samples.

Morphological asymmetries of mouse brain assessed by geometric morphometric analysis of MRI data.

Mammalian brain has repeated structures at both sides of the median plane, although some asymmetries have been described even under normal conditions. Characterizing normal patterns of asymmetry in mouse brain is important to recognize features that depart from expected ranges in the most widely used mammalian model. Analyses on brain morphology based on magnetic resonance image (MRI) have largely focused on volumes while less is known about shape asymmetry. We introduce a flexible protocol based on geometric morphometrics to assess patterns of asymmetry in shape and size of mouse brain from microMRI scans. After systematic digitization of landmarks and semilandmarks, we combine multivariate methods for statistical analyses with visualization tools to display the results. No preliminary treatment of the images (e.g. space normalization) is needed to collect data on MRI slices and visual representations improve the interpretation of the results. Results indicated that the protocol is highly repeatable. Asymmetry was more evident for shape than for size. Particularly, fluctuating asymmetry accounted for more variation than directional asymmetry in all brain regions. Since this approach can detect subtle shape variation between sides, it is a promising methodology to explore morphological changes in the brain of model organisms and can be applied in future studies addressing the effect of genetic and environmental factors on brain morphology.