Genetics of murine craniofacial morphology: diallel analysis of the eight founders of the Collaborative Cross

Palatal segment contributions to midfacial anterior-posterior growth

Anterior-posterior (A-P) elongation of the palate is a critical aspect of integrated midfacial morphogenesis. Reciprocal epithelial-mesenchymal interactions drive secondary palate elongation that is coupled to the periodic formation of signaling centers within the rugae growth zone (RGZ). However, the relationship between RGZ driven morphogenetic processes, the differentiative dynamics of underlying palatal bone mesenchymal precursors, and the segmental organization of the upper jaw has remained enigmatic. A detailed ontogenetic study of these relationships is important, because palatal segment growth is a critical aspect of normal midfacial growth, can be modified to produce dysmorphology, and is a likely basis for evolutionary differences in upper jaw morphology. Variation in palatal-segment specific growth may also underlie known differences in palatal segment proportions between inbred mouse strains. We completed a combined whole mount gene expression and morphometric analysis of normal murine palatal growth dynamics and their association with palatal segment elongation and resulting upper jaw morphology. Our results demonstrated that the first formed palatal ruga (ruga 1), found just posterior to the RGZ, maintained an association with important nasal, neurovascular and palatal structures throughout early midfacial development; suggesting that these features are positioned at a proximal source of embryonic midfacial directional growth. Our detailed characterization of midfacial morphogenesis revealed a one-to-one relationship between palatal segments and upper jaw bones during the earliest stages of palatal elongation. Growth of the maxillary anlage within the anterior secondary palate is uniquely coupled to RGZ-driven morphogenesis that more than doubles the length of this palatal segment prior to palatal shelf fusion. Our results also demonstrate that the future maxillary-palatine suture, approximated by the position ruga 1 and consistently associated with the palatine anlage, forms predominantly via the posterior differentiation of the maxilla within the expanding anterior secondary palate. Our complementary ontogenetic comparison of three inbred mouse strains identified small but significant strain-specific differences in early embryonic palatal segment contributions to the upper jaw. Although early palatal segment specific growth is not primarily responsible for adult differences in upper jaw morphology between these strains, our ontogenetic series of measurements provide a useful foundation for understanding the impact of background genetic effects on facial shape and elongation. In combination, our results provide a novel and particularly detailed picture of the earliest spatiotemporal dynamics of intramembranous midfacial skeletal specification and differentiation within the context of the surrounding palatal segment A-P elongation and associated rugae formation.

Morphometric analysis of the size-adjusted linear dimensions of the skull landmarks revealed craniofacial dysmorphology in Mid1-cKO mice

BACKGROUND

The early craniofacial development is a highly coordinated process involving neural crest cell migration, proliferation, epithelial apoptosis, and epithelial-mesenchymal transition (EMT). Both genetic defects and environmental factors can affect these processes and result in orofacial clefts. Mutations in MID1 gene cause X-linked Opitz Syndrome (OS), which is a congenital malformation characterized by craniofacial defects including cleft lip/palate (CLP). Previous studies demonstrated impaired neurological structure and function in Mid1 knockout mice, while no CLP was observed. However, given the highly variable severities of the facial manifestations observed in OS patients within the same family carrying identical genetic defects, subtle craniofacial malformations in Mid1 knockout mice could be overlooked in these studies. Therefore, we propose that a detailed morphometric analysis should be necessary to reveal mild craniofacial dysmorphologies that reflect the similar developmental defects seen in OS patients.

RESULTS

In this research, morphometric study of the P0 male Mid1-cKO mice were performed using Procrustes superimposition as well as EMDA analysis of the size-adjusted three-dimensional coordinates of 105 skull landmarks, which were collected on the bone surface reconstructed using microcomputed tomographic images. Our results revealed the craniofacial deformation such as the increased dimension of the frontal and nasal bone in Mid1-cKO mice, in line with the most prominent facial features such as hypertelorism, prominent forehead, broad and/or high nasal bridge seen in OS patients.

CONCLUSION

 While been extensively used in evolutionary biology and anthropology in the last decades, geometric morphometric analysis was much less used in developmental biology. Given the high interspecies variances in facial anatomy, the work presented in this research suggested the advantages of morphometric analysis in characterizing animal models of craniofacial developmental defects to reveal phenotypic variations and the underlining pathogenesis.

Merging morphological and genetic evidence to assess hybridization in Western Eurasian late Pleistocene hominins

Previous scientific consensus saw human evolution as defined by adaptive differences (behavioural and/or biological) and the emergence of Homo sapiens as the ultimate replacement of non-modern groups by a modern, adaptively more competitive group. However, recent research has shown that the process underlying our origins was considerably more complex. While archaeological and fossil evidence suggests that behavioural complexity may not be confined to the modern human lineage, recent palaeogenomic work shows that gene flow between distinct lineages (for example, Neanderthals, Denisovans, early H. sapiens) occurred repeatedly in the late Pleistocene, probably contributing elements to our genetic make-up that might have been crucial to our success as a diverse, adaptable species. Following these advances, the prevailing human origins model has shifted from one of near-complete replacement to a more nuanced view of partial replacement with considerable reticulation. Here we provide a brief introduction to the current genetic evidence for hybridization among hominins, its prevalence in, and effects on, comparative mammal groups, and especially how it manifests in the skull. We then explore the degree to which cranial variation seen in the fossil record of late Pleistocene hominins from Western Eurasia corresponds with our current genetic and comparative data. We are especially interested in understanding the degree to which skeletal data can reflect admixture. Our findings indicate some correspondence between these different lines of evidence, flag individual fossils as possibly admixed, and suggest that different cranial regions may preserve hybridization signals differentially. We urge further studies of the phenotype to expand our ability to detect the ways in which migration, interaction and genetic exchange have shaped the human past, beyond what is currently visible with the lens of ancient DNA.

The genetic basis of neurocranial size and shape across varied lab mouse populations

Brain and skull tissues interact through molecular signalling and mechanical forces during head development, leading to a strong correlation between the neurocranium and the external brain surface. Therefore, when brain tissue is unavailable, neurocranial endocasts are often used to approximate brain size and shape. Evolutionary changes in brain morphology may have resulted in secondary changes to neurocranial morphology, but the developmental and genetic processes underlying this relationship are not well understood. Using automated phenotyping methods, we quantified the genetic basis of endocast variation across large genetically varied populations of laboratory mice in two ways: (1) to determine the contributions of various genetic factors to neurocranial form and (2) to help clarify whether a neurocranial variation is based on genetic variation that primarily impacts bone development or on genetic variation that primarily impacts brain development, leading to secondary changes in bone morphology. Our results indicate that endocast size is highly heritable and is primarily determined by additive genetic factors. In addition, a non-additive inbreeding effect led to founder strains with lower neurocranial size, but relatively large brains compared to skull size; suggesting stronger canalization of brain size and/or a general allometric effect. Within an outbred sample of mice, we identified a locus on mouse chromosome 1 that is significantly associated with variation in several positively correlated endocast size measures. Because the protein-coding genes at this locus have been previously associated with brain development and not with bone development, we propose that genetic variation at this locus leads primarily to variation in brain volume that secondarily leads to changes in neurocranial globularity. We identify a strain-specific missense mutation within Akt3 that is a strong causal candidate for this genetic effect. Whilst it is not appropriate to generalize our hypothesis for this single locus to all other loci that also contribute to the complex trait of neurocranial skull morphology, our results further reveal the genetic basis of neurocranial variation and highlight the importance of the mechanical influence of brain growth in determining skull morphology.

Modelling the shape of the pig scapula

BACKGROUND

The shape of pig scapula is complex and is important for sow robustness and health. To better understand the relationship between 3D shape of the scapula and functional traits, it is necessary to build a model that explains most of the morphological variation between animals. This requires point correspondence, i.e. a map that explains which points represent the same piece of tissue among individuals. The objective of this study was to further develop an automated computational pipeline for the segmentation of computed tomography (CT) scans to incorporate 3D modelling of the scapula, and to develop a genetic prediction model for 3D morphology.

RESULTS

The surface voxels of the scapula were identified on 2143 CT-scanned pigs, and point correspondence was established by predicting the coordinates of 1234 semi-landmarks on each animal, using the coherent point drift algorithm. A subsequent principal component analysis showed that the first 10 principal components covered more than 80% of the total variation in 3D shape of the scapula. Using principal component scores as phenotypes in a genetic model, estimates of heritability ranged from 0.4 to 0.8 (with standard errors from 0.07 to 0.08). To validate the entire computational pipeline, a statistical model was trained to predict scapula shape based on marker genotype data. The mean prediction reliability averaged over the whole scapula was equal to 0.18 (standard deviation = 0.05) with a higher reliability in convex than in concave regions.

CONCLUSIONS

Estimates of heritability of the principal components were high and indicated that the computational pipeline that processes CT data to principal component phenotypes was associated with little error. Furthermore, we showed that it is possible to predict the 3D shape of scapula based on marker genotype data. Taken together, these results show that the proposed computational pipeline closes the gap between a point cloud representing the shape of an animal and its underlying genetic components.

A comprehensive and comparative phenotypic analysis of the collaborative founder strains identifies new and known phenotypes

The collaborative cross (CC) is a large panel of mouse-inbred lines derived from eight founder strains (NOD/ShiLtJ, NZO/HILtJ, A/J, C57BL/6J, 129S1/SvImJ, CAST/EiJ, PWK/PhJ, and WSB/EiJ). Here, we performed a comprehensive and comparative phenotyping screening to identify phenotypic differences and similarities between the eight founder strains. In total, more than 300 parameters including allergy, behavior, cardiovascular, clinical blood chemistry, dysmorphology, bone and cartilage, energy metabolism, eye and vision, immunology, lung function, neurology, nociception, and pathology were analyzed; in most traits from sixteen females and sixteen males. We identified over 270 parameters that were significantly different between strains. This study highlights the value of the founder and CC strains for phenotype-genotype associations of many genetic traits that are highly relevant to human diseases. All data described here are publicly available from the mouse phenome database for analyses and downloads.

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.

Hybridization in human evolution: Insights from other organisms

During the late Pleistocene, isolated lineages of hominins exchanged genes thus influencing genomic variation in humans in both the past and present. However, the dynamics of this genetic exchange and associated phenotypic consequences through time remain poorly understood. Gene exchange across divergent lineages can result in myriad outcomes arising from these dynamics and the environmental conditions under which it occurs. Here we draw from our collective research across various organisms, illustrating some of the ways in which gene exchange can structure genomic/phenotypic diversity within/among species. We present a range of examples relevant to questions about the evolution of hominins. These examples are not meant to be exhaustive, but rather illustrative of the diverse evolutionary causes/consequences of hybridization, highlighting potential drivers of human evolution in the context of hybridization including: influences on adaptive evolution, climate change, developmental systems, sex-differences in behavior, Haldane's rule and the large X-effect, and transgressive phenotypic variation.

A Diallel of the Mouse Collaborative Cross Founders Reveals Strong Strain-Specific Maternal Effects on Litter Size

Reproductive success in the eight founder strains of the Collaborative Cross (CC) was measured using a diallel-mating scheme. Over a 48-month period we generated 4,448 litters, and provided 24,782 weaned pups for use in 16 different published experiments. We identified factors that affect the average litter size in a cross by estimating the overall contribution of parent-of-origin, heterosis, inbred, and epistatic effects using a Bayesian zero-truncated overdispersed Poisson mixed model. The phenotypic variance of litter size has a substantial contribution (82%) from unexplained and environmental sources, but no detectable effect of seasonality. Most of the explained variance was due to additive effects (9.2%) and parental sex (maternal vs. paternal strain; 5.8%), with epistasis accounting for 3.4%. Within the parental effects, the effect of the dam's strain explained more than the sire's strain (13.2% vs. 1.8%), and the dam's strain effects account for 74.2% of total variation explained. Dams from strains C57BL/6J and NOD/ShiLtJ increased the expected litter size by a mean of 1.66 and 1.79 pups, whereas dams from strains WSB/EiJ, PWK/PhJ, and CAST/EiJ reduced expected litter size by a mean of 1.51, 0.81, and 0.90 pups. Finally, there was no strong evidence for strain-specific effects on sex ratio distortion. Overall, these results demonstrate that strains vary substantially in their reproductive ability depending on their genetic background, and that litter size is largely determined by dam's strain rather than sire's strain effects, as expected. This analysis adds to our understanding of factors that influence litter size in mammals, and also helps to explain breeding successes and failures in the extinct lines and surviving CC strains.

The effect of automated landmark identification on morphometric analyses

Morphometric analysis of anatomical landmarks allows researchers to identify specific morphological differences between natural populations or experimental groups, but manually identifying landmarks is time-consuming. We compare manually and automatically generated adult mouse skull landmarks and subsequent morphometric analyses to elucidate how switching from manual to automated landmarking will impact morphometric analysis results for large mouse (Mus musculus) samples (n = 1205) that represent a wide range of 'normal' phenotypic variation (62 genotypes). Other studies have suggested that the use of automated landmarking methods is feasible, but this study is the first to compare the utility of current automated approaches to manual landmarking for a large dataset that allows the quantification of intra- and inter-strain variation. With this unique sample, we investigated how switching to a non-linear image registration-based automated landmarking method impacts estimated differences in genotype mean shape and shape variance-covariance structure. In addition, we tested whether an initial registration of specimen images to genotype-specific averages improves automatic landmark identification accuracy. Our results indicated that automated landmark placement was significantly different than manual landmark placement but that estimated skull shape covariation was correlated across methods. The addition of a preliminary genotype-specific registration step as part of a two-level procedure did not substantially improve on the accuracy of one-level automatic landmark placement. The landmarks with the lowest automatic landmark accuracy are found in locations with poor image registration alignment. The most serious outliers within morphometric analysis of automated landmarks displayed instances of stochastic image registration error that are likely representative of errors common when applying image registration methods to micro-computed tomography datasets that were initially collected with manual landmarking in mind. Additional efforts during specimen preparation and image acquisition can help reduce the number of registration errors and improve registration results. A reduction in skull shape variance estimates were noted for automated landmarking methods compared with manual landmarking. This partially reflects an underestimation of more extreme genotype shapes and loss of biological signal, but largely represents the fact that automated methods do not suffer from intra-observer landmarking error. For appropriate samples and research questions, our image registration-based automated landmarking method can eliminate the time required for manual landmarking and have a similar power to identify shape differences between inbred mouse genotypes.

Creating diversity in mammalian facial morphology: a review of potential developmental mechanisms

Mammals (class Mammalia) have evolved diverse craniofacial morphology to adapt to a wide range of ecological niches. However, the genetic and developmental mechanisms underlying the diversification of mammalian craniofacial morphology remain largely unknown. In this paper, we focus on the facial length and orofacial clefts of mammals and deduce potential mechanisms that produced diversity in mammalian facial morphology. Small-scale changes in facial morphology from the common ancestor, such as slight changes in facial length and the evolution of the midline cleft in some lineages of bats, could be attributed to heterochrony in facial bone ossification. In contrast, large-scale changes of facial morphology from the common ancestor, such as a truncated, widened face as well as the evolution of the bilateral cleft possessed by some bat species, could be brought about by changes in growth and patterning of the facial primordium (the facial processes) at the early stages of embryogenesis.

Craniomandibular form and body size variation of first generation mouse hybrids: A model for hominin hybridization

Hybridization occurs in a number of mammalian lineages, including among primate taxa. Analyses of ancient genomes have shown that hybridization between our lineage and other archaic hominins in Eurasia occurred numerous times in the past. However, we still have limited empirical data on what a hybrid skeleton looks like, or how to spot patterns of hybridization among fossils for which there are no genetic data. Here we use experimental mouse models to supplement previous studies of primates. We characterize size and shape variation in the cranium and mandible of three wild-derived inbred mouse strains and their first generation (F1) hybrids. The three parent taxa in our analysis represent lineages that diverged over approximately the same period as the human/Neanderthal/Denisovan lineages and their hybrids are variably successful in the wild. Comparisons of body size, as quantified by long bone measurements, are also presented to determine whether the identified phenotypic effects of hybridization are localized to the cranium or represent overall body size changes. The results indicate that hybrid cranial and mandibular sizes, as well as limb length, exceed that of the parent taxa in all cases. All three F1 hybrid crosses display similar patterns of size and form variation. These results are generally consistent with earlier studies on primates and other mammals, suggesting that the effects of hybridization may be similar across very different scenarios of hybridization, including different levels of hybrid fitness. This paper serves to supplement previous studies aimed at identifying F1 hybrids in the fossil record and to introduce further research that will explore hybrid morphologies using mice as a proxy for better understanding hybridization in the hominin fossil record.

Developmental constraint through negative pleiotropy in the zygomatic arch

BACKGROUND

Previous analysis suggested that the relative contribution of individual bones to regional skull lengths differ between inbred mouse strains. If the negative correlation of adjacent bone lengths is associated with genetic variation in a heterogeneous population, it would be an example of negative pleiotropy, which occurs when a genetic factor leads to opposite effects in two phenotypes. Confirming negative pleiotropy and determining its basis may reveal important information about the maintenance of overall skull integration and developmental constraint on skull morphology.

RESULTS

We identified negative correlations between the lengths of the frontal and parietal bones in the midline cranial vault as well as the zygomatic bone and zygomatic process of the maxilla, which contribute to the zygomatic arch. Through gene association mapping of a large heterogeneous population of Diversity Outbred (DO) mice, we identified a quantitative trait locus on chromosome 17 driving the antagonistic contribution of these two zygomatic arch bones to total zygomatic arch length. Candidate genes in this region were identified and real-time PCR of the maxillary processes of DO founder strain embryos indicated differences in the RNA expression levels for two of the candidate genes, Camkmt and Six2.

CONCLUSIONS

A genomic region underlying negative pleiotropy of two zygomatic arch bones was identified, which provides a mechanism for antagonism in component bone lengths while constraining overall zygomatic arch length. This type of mechanism may have led to variation in the contribution of individual bones to the zygomatic arch noted across mammals. Given that similar genetic and developmental mechanisms may underlie negative correlations in other parts of the skull, these results provide an important step toward understanding the developmental basis of evolutionary variation and constraint in skull morphology.

A population level atlas of Mus musculus craniofacial skeleton and automated image-based shape analysis

Laboratory mice are staples for evo/devo and genetics studies. Inbred strains provide a uniform genetic background to manipulate and understand gene-environment interactions, while their crosses have been instrumental in studies of genetic architecture, integration and modularity, and mapping of complex biological traits. Recently, there have been multiple large-scale studies of laboratory mice to further our understanding of the developmental basis, evolution, and genetic control of shape variation in the craniofacial skeleton (i.e. skull and mandible). These experiments typically use micro-computed tomography (micro-CT) to capture the craniofacial phenotype in 3D and rely on manually annotated anatomical landmarks to conduct statistical shape analysis. Although the common choice for imaging modality and phenotyping provides the potential for collaborative research for even larger studies with more statistical power, the investigator (or lab-specific) nature of the data collection hampers these efforts. Investigators are rightly concerned that subtle differences in how anatomical landmarks were recorded will create systematic bias between studies that will eventually influence scientific findings. Even if researchers are willing to repeat landmark annotation on a combined dataset, different lab practices and software choices may create obstacles for standardization beyond the underlying imaging data. Here, we propose a freely available analysis system that could assist in the standardization of micro-CT studies in the mouse. Our proposal uses best practices developed in biomedical imaging and takes advantage of existing open-source software and imaging formats. Our first contribution is the creation of a synthetic template for the adult mouse craniofacial skeleton from 25 inbred strains and five F1 crosses that are widely used in biological research. The template contains a fully segmented cranium, left and right hemi-mandibles, endocranial space, and the first few cervical vertebrae. We have been using this template in our lab to segment and isolate cranial structures in an automated fashion from a mixed population of mice, including craniofacial mutants, aged 4-12.5 weeks. As a secondary contribution, we demonstrate an application of nearly automated shape analysis, using symmetric diffeomorphic image registration. This approach, which we call diGPA, closely approximates the popular generalized Procrustes analysis (GPA) but negates the collection of anatomical landmarks. We achieve our goals by using the open-source advanced normalization tools (ANT) image quantification library, as well as its associated R library (ANTsR) for statistical image analysis. Finally, we make a plea to investigators to commit to using open imaging standards and software in their labs to the extent possible to increase the potential for data exchange and improve the reproducibility of findings. Future work will incorporate more anatomical detail (such as individual cranial bones, turbinals, dentition, middle ear ossicles) and more diversity into the template.

Integration of Brain and Skull in Prenatal Mouse Models of Apert and Crouzon Syndromes

The brain and skull represent a complex arrangement of integrated anatomical structures composed of various cell and tissue types that maintain structural and functional association throughout development. Morphological integration, a concept developed in vertebrate morphology and evolutionary biology, describes the coordinated variation of functionally and developmentally related traits of organisms. Syndromic craniosynostosis is characterized by distinctive changes in skull morphology and perceptible, though less well studied, changes in brain structure and morphology. Using mouse models for craniosynostosis conditions, our group has precisely defined how unique craniosynostosis causing mutations in fibroblast growth factor receptors affect brain and skull morphology and dysgenesis involving coordinated tissue-specific effects of these mutations. Here we examine integration of brain and skull in two mouse models for craniosynostosis: one carrying the FGFR2c C342Y mutation associated with Pfeiffer and Crouzon syndromes and a mouse model carrying the FGFR2 S252W mutation, one of two mutations responsible for two-thirds of Apert syndrome cases. Using linear distances estimated from three-dimensional coordinates of landmarks acquired from dual modality imaging of skull (high resolution micro-computed tomography and magnetic resonance microscopy) of mice at embryonic day 17.5, we confirm variation in brain and skull morphology in Fgfr2c^C342Y/+ mice, Fgfr2^+/S252W mice, and their unaffected littermates. Mutation-specific variation in neural and cranial tissue notwithstanding, patterns of integration of brain and skull differed only subtly between mice carrying either the FGFR2c C342Y or the FGFR2 S252W mutation and their unaffected littermates. However, statistically significant and substantial differences in morphological integration of brain and skull were revealed between the two mutant mouse models, each maintained on a different strain. Relative to the effects of disease-associated mutations, our results reveal a stronger influence of the background genome on patterns of brain-skull integration and suggest robust genetic, developmental, and evolutionary relationships between neural and skeletal tissues of the head.

The Interaction of Genetic Background and Mutational Effects in Regulation of Mouse Craniofacial Shape

Inbred genetic background significantly influences the expression of phenotypes associated with known genetic perturbations and can underlie variation in disease severity between individuals with the same mutation. However, the effect of epistatic interactions on the development of complex traits, such as craniofacial morphology, is poorly understood. Here, we investigated the effect of three inbred backgrounds (129X1/SvJ, C57BL/6J, and FVB/NJ) on the expression of craniofacial dysmorphology in mice (Mus musculus) with loss of function in three members of the Sprouty family of growth factor negative regulators (Spry1, Spry2, or Spry4) in order to explore the impact of epistatic interactions on skull morphology. We found that the interaction of inbred background and the Sprouty genotype explains as much craniofacial shape variation as the Sprouty genotype alone. The most severely affected genotypes display a relatively short and wide skull, a rounded cranial vault, and a more highly angled inferior profile. Our results suggest that the FVB background is more resilient to Sprouty loss of function than either C57 or 129, and that Spry4 loss is generally less severe than loss of Spry1 or Spry2 While the specific modifier genes responsible for these significant background effects remain unknown, our results highlight the value of intercrossing mice of multiple inbred backgrounds to identify the genes and developmental interactions that modulate the severity of craniofacial dysmorphology. Our quantitative results represent an important first step toward elucidating genetic interactions underlying variation in robustness to known genetic perturbations in mice.

Angiotensin receptor blockade mediated amelioration of mucopolysaccharidosis type I cardiac and craniofacial pathology

Mucopolysaccharidosis type I (MPS IH) is a lysosomal storage disease (LSD) caused by inactivating mutations to the alpha-L-iduronidase (IDUA) gene. Treatment focuses on IDUA enzyme replacement and currently employed methods can be non-uniform in their efficacy particularly for the cardiac and craniofacial pathology. Therefore, we undertook efforts to better define the pathological cascade accounting for treatment refractory manifestations and demonstrate a role for the renin angiotensin system (RAS) using the IDUA^-/- mouse model. Perturbation of the RAS in the aorta was more profound in male animals suggesting a causative role in the observed gender dimorphism and angiotensin receptor blockade (ARB) resulted in improved cardiac function. Further, we show the ability of losartan to prevent shortening of the snout, a common craniofacial anomaly in IDUA^-/- mice. These data show a key role for the RAS in MPS associated pathology and support the inclusion of losartan as an augmentation to current therapies.

Development Shapes a Consistent Inbreeding Effect in Mouse Crania of Different Line Crosses

Development translates genetic variation into a multivariate pattern of phenotypic variation, distributing it among traits in a nonuniform manner. As developmental processes are largely shared within species, this suggests that heritable phenotypic variation will be patterned similarly, in spite of the different segregating alleles. To investigate developmental effect on the variational pattern in the shape of the mouse skull across genetically differentiated lines, we employed the full set of reciprocal crosses (a.k.a. diallel) between eight inbred mouse strains of the Collaborative Cross Project. We used geometric morphometrics and multivariate analysis to capture cranial size and shape changes in 8 parentals and their 54 F1 crosses. The high heterozygosity generated in the F1 crosses allowed us to compare the multivariate deviations of the F1 phenotypes from the expected midparental phenotypes in different haplotype combinations. In contrast to body weight, we found a high degree of nonadditive deviation in craniofacial shape. Whereas the phenotypic and genetic divergence of parental strains manifested in high dimensionality of additive effects, the nonadditive deviations exhibited lesser dimensionality and in particular a strikingly coherent direction in shape space. We interpret this finding as evidence for a strong structuring effect of a relatively small set of developmental processes on the mapping of genetic to phenotypic variation.

Resources for Systems Genetics

A key characteristic of systems genetics is its reliance on populations that vary to a greater or lesser degree in genetic complexity-from highly admixed populations such as the Collaborative Cross and Diversity Outcross to relatively simple crosses such as sets of consomic strains and reduced complexity crosses. This protocol is intended to help investigators make more informed decisions about choices of resources given different types of questions. We consider factors such as costs, availability, and ease of breeding for common scenarios. In general, we recommend using complementary resources and minimizing depth of resampling of any given genome or strain.

Human Facial Shape and Size Heritability and Genetic Correlations

The human face is an array of variable physical features that together make each of us unique and distinguishable. Striking familial facial similarities underscore a genetic component, but little is known of the genes that underlie facial shape differences. Numerous studies have estimated facial shape heritability using various methods. Here, we used advanced three-dimensional imaging technology and quantitative human genetics analysis to estimate narrow-sense heritability, heritability explained by common genetic variation, and pairwise genetic correlations of 38 measures of facial shape and size in normal African Bantu children from Tanzania. Specifically, we fit a linear mixed model of genetic relatedness between close and distant relatives to jointly estimate variance components that correspond to heritability explained by genome-wide common genetic variation and variance explained by uncaptured genetic variation, the sum representing total narrow-sense heritability. Our significant estimates for narrow-sense heritability of specific facial traits range from 28 to 67%, with horizontal measures being slightly more heritable than vertical or depth measures. Furthermore, for over half of facial traits, >90% of narrow-sense heritability can be explained by common genetic variation. We also find high absolute genetic correlation between most traits, indicating large overlap in underlying genetic loci. Not surprisingly, traits measured in the same physical orientation (i.e., both horizontal or both vertical) have high positive genetic correlations, whereas traits in opposite orientations have high negative correlations. The complex genetic architecture of facial shape informs our understanding of the intricate relationships among different facial features as well as overall facial development.

Dissecting the Effect of Genetic Variation on the Hepatic Expression of Drug Disposition Genes across the Collaborative Cross Mouse Strains

A central challenge in pharmaceutical research is to investigate genetic variation in response to drugs. The Collaborative Cross (CC) mouse reference population is a promising model for pharmacogenomic studies because of its large amount of genetic variation, genetic reproducibility, and dense recombination sites. While the CC lines are phenotypically diverse, their genetic diversity in drug disposition processes, such as detoxification reactions, is still largely uncharacterized. Here we systematically measured RNA-sequencing expression profiles from livers of 29 CC lines under baseline conditions. We then leveraged a reference collection of metabolic biotransformation pathways to map potential relations between drugs and their underlying expression quantitative trait loci (eQTLs). By applying this approach on proximal eQTLs, including eQTLs acting on the overall expression of genes and on the expression of particular transcript isoforms, we were able to construct the organization of hepatic eQTL-drug connectivity across the CC population. The analysis revealed a substantial impact of genetic variation acting on drug biotransformation, allowed mapping of potential joint genetic effects in the context of individual drugs, and demonstrated crosstalk between drug metabolism and lipid metabolism. Our findings provide a resource for investigating drug disposition in the CC strains, and offer a new paradigm for integrating biotransformation reactions to corresponding variations in DNA sequences.

Genetic structure of phenotypic robustness in the collaborative cross mouse diallel panel

Developmental stability and canalization describe the ability of developmental systems to minimize phenotypic variation in the face of stochastic micro-environmental effects, genetic variation and environmental influences. Canalization is the ability to minimize the effects of genetic or environmental effects, whereas developmental stability is the ability to minimize the effects of micro-environmental effects within individuals. Despite much attention, the mechanisms that underlie these two components of phenotypic robustness remain unknown. We investigated the genetic structure of phenotypic robustness in the collaborative cross (CC) mouse reference population. We analysed the magnitude of fluctuating asymmetry (FA) and among-individual variation of cranial shape in reciprocal crosses among the eight parental strains, using geometric morphometrics and a diallel analysis based on a Bayesian approach. Significant differences among genotypes were found for both measures, although they were poorly correlated at the level of individuals. An overall positive effect of inbreeding was found for both components of variation. The strain CAST/EiJ exerted a positive additive effect on FA and, to a lesser extent, among-individual variance. Sex- and other strain-specific effects were not significant. Neither FA nor among-individual variation was associated with phenotypic extremeness. Our results support the existence of genetic variation for both developmental stability and canalization. This finding is important because robustness is a key feature of developmental systems. Our finding that robustness is not related to phenotypic extremeness is consistent with theoretical work that suggests that its relationship to stabilizing selection is not straightforward.

Craniofacial shape transition across the house mouse hybrid zone: implications for the genetic architecture and evolution of between-species differences

Craniofacial shape differences between taxa have often been linked to environmental adaptation, e.g., new food sources, or have been studied in the context of domestication. Evidence for the genetic basis of such phenotypic differences to date suggests that between-species as well as between-population variation has an oligogenic basis, i.e., few loci of large effect explain most of the variation. In mice, it has been shown that within-population craniofacial variation has a highly polygenic basis, but there are no data regarding the genetic basis of between-species differences in natural populations. Here, we address this question using a phenotype-focused approach. Using 3D geometric morphometrics, we phenotyped a panel of mice derived from a natural hybrid zone between Mus musculus domesticus and Mus mus musculus and quantify the transition of craniofacial shape along the hybridization gradient. We find a continuous shape transition along the hybridization gradient and unaltered developmental stability associated with hybridization. This suggests that the morphospace between the two subspecies is continuous despite reproductive isolation and strong barriers to gene flow. We show that quantitative changes in overall genome composition generate quantitative changes in craniofacial shape; this supports a highly polygenic basis for between-species craniofacial differences in the house mouse. We discuss our findings in the context of oligogenic versus polygenic models of the genetic architecture of morphological traits.

Does 3D Phenotyping Yield Substantial Insights in the Genetics of the Mouse Mandible Shape?

We describe the application of high-resolution 3D microcomputed tomography, together with 3D landmarks and geometric morphometrics, to validate and further improve previous quantitative genetic studies that reported QTL responsible for variation in the mandible shape of laboratory mice using a new backcross between C57BL/6J and A/J inbred strains. Despite the increasing availability of 3D imaging techniques, artificial flattening of the mandible by 2D imaging techniques seems at first an acceptable compromise for large-scale phenotyping protocols, thanks to an abundance of low-cost digital imaging systems such as microscopes or digital cameras. We evaluated the gain of information from considering explicitly this additional third dimension, and also from capturing variation on the bone surface where no precise anatomical landmark can be marked. Multivariate QTL mapping conducted with different landmark configurations (2D vs. 3D; manual vs. semilandmarks) broadly agreed with the findings of previous studies. Significantly more QTL (23) were identified and more precisely mapped when the mandible shape was captured with a large set of semilandmarks coupled with manual landmarks. It appears that finer phenotypic characterization of the mandibular shape with 3D landmarks, along with higher density genotyping, yields better insights into the genetic architecture of mandibular development. Most of the main variation is, nonetheless, preferentially embedded in the natural 2D plane of the hemi-mandible, reinforcing the results of earlier influential investigations.

The old and new face of craniofacial research: How animal models inform human craniofacial genetic and clinical data

The craniofacial skeletal structures that comprise the human head develop from multiple tissues that converge to form the bones and cartilage of the face. Because of their complex development and morphogenesis, many human birth defects arise due to disruptions in these cellular populations. Thus, determining how these structures normally develop is vital if we are to gain a deeper understanding of craniofacial birth defects and devise treatment and prevention options. In this review, we will focus on how animal model systems have been used historically and in an ongoing context to enhance our understanding of human craniofacial development. We do this by first highlighting "animal to man" approaches; that is, how animal models are being utilized to understand fundamental mechanisms of craniofacial development. We discuss emerging technologies, including high throughput sequencing and genome editing, and new animal repository resources, and how their application can revolutionize the future of animal models in craniofacial research. Secondly, we highlight "man to animal" approaches, including the current use of animal models to test the function of candidate human disease variants. Specifically, we outline a common workflow deployed after discovery of a potentially disease causing variant based on a select set of recent examples in which human mutations are investigated in vivo using animal models. Collectively, these topics will provide a pipeline for the use of animal models in understanding human craniofacial development and disease for clinical geneticist and basic researchers alike.