Quantitative trait loci affecting the 3D skull shape and size in mouse and prioritization of candidate genes in-silico

Craniofacial anomalies in schizophrenia-relevant GFAP.HMOX10-12m mice

Subtle craniofacial dysmorphology has been reported in schizophrenia patients. This dysmorphology includes midline facial elongation, frontonasal anomalies and a sexually dimorphic deviation from normal directional asymmetry of the face, with male patients showing reduced and female patients showing enhanced facial asymmetry relative to healthy control subjects. GFAP.HMOX10-12m transgenic mice (Mus musculus) that overexpress heme oxygenase-1 in astrocytes recapitulate many schizophrenia-relevant neurochemical, neuropathological and behavioral features. As morphogenesis of the brain, skull and face are highly interrelated, we hypothesized that GFAP.HMOX10-12m mice may exhibit craniofacial anomalies similar to those reported in persons with schizophrenia. We examined craniofacial anatomy in male GFAP.HMOX10-12m mice and wild-type control mice at the early adulthood age of 6-8 months. We used computer vision techniques for the extraction and analysis of mouse head shape parameters from systematically acquired 2D digital images, and confirmed our results with landmark-based geometric morphometrics. We performed skull bone morphometry using digital calipers to take linear distance measurements between known landmarks. Relative to controls, adult male GFAP.HMOX10-12m mice manifested craniofacial dysmorphology including elongation of the nasal bones, alteration of head shape anisotropy and reduction of directional asymmetry in facial shape features. These findings demonstrate that GFAP.HMOX10-12m mice exhibit craniofacial anomalies resembling those described in schizophrenia patients, implicating heme oxygenase-1 in their development. As a preclinical mouse model, GFAP.HMOX10-12m mice provide a novel opportunity for the study of the etiopathogenesis of craniofacial and other anomalies in schizophrenia and related disorders.

Natural selection on feralization genes contributed to the invasive spread of wild pigs throughout the United States

Despite a long presence in the contiguous United States (US), the distribution of invasive wild pigs (Sus scrofa × domesticus) has expanded rapidly since the 1980s, suggesting a more recent evolutionary shift towards greater invasiveness. Contemporary populations of wild pigs represent exoferal hybrid descendants of domestic pigs and European wild boar, with such hybridization expected to enrich genetic diversity and increase the adaptive potential of populations. Our objective was to characterize how genetic enrichment through hybridization increases the invasiveness of populations by identifying signals of selection and the ancestral origins of selected loci. Our study focused on invasive wild pigs within Great Smoky Mountains National Park, which represents a hybrid population descendent from the admixture of established populations of feral pigs and an introduction of European wild boar to North America. Accordingly, we genotyped 881 wild pigs with multiple high-density single-nucleotide polymorphism (SNP) arrays. We found 233 markers under putative selection spread over 79 regions across 16 out of 18 autosomes, which contained genes involved in traits affecting feralization. Among these, genes were found to be related to skull formation and neurogenesis, with two genes, TYRP1 and TYR, also encoding for crucial melanogenesis enzymes. The most common haplotypes associated with regions under selection for the Great Smoky Mountains population were also common among other populations throughout the region, indicating a key role of putatively selective variants in the fitness of invasive populations. Interestingly, many of these haplotypes were absent among European wild boar reference genotypes, indicating feralization through genetic adaptation.

Comparative Cranial Geometric Morphometrics among Wistar Albino, Sprague Dawley, and WAG/Rij Rat Strains

This research utilizes geometric morphometrics to investigate shape variation in the skull, mandible, and teeth among three rat strains: Wistar Albino (WA), Sprague Dawley (SD), and WAG/Rij (WR). Through the analysis of 48 rats using 2D geometric morphometric techniques, significant differences in their skull morphology were identified. This study indicates a shift from a rectangular to an oval cranial shape across strains, with notable size and morphological variances. Particularly, the WR strain's skull shape significantly differs from the SD and WA strains, suggesting distinct ecological or genetic pathways. Compared to the skull, mandible shape differences are less pronounced, but still significant. The WR strain exhibits a distinct mandible shape, potentially reflecting ecological adaptations like dietary habits. The teeth shape of WR rats is the most distinct. SD rats consistently exhibited larger sizes in both skull and mandible measurements, while WR rats were notably smaller. Interestingly, sexual dimorphism was not statistically significant in skull and teeth sizes, aligning with findings from previous studies. However, the mandible showed clear size differences between sexes, underscoring its potential for adaptive or behavioral studies. In summary, this study provides a comprehensive analysis of morphological variations in rat strains, highlighting the intricate interplay of size, shape, and ecological factors. These findings lay a foundation for deeper explorations into the adaptive, ecological, or genetic narratives influencing rat morphology.

Joint multi-ancestry and admixed GWAS reveals the complex genetics behind human cranial vault shape

The cranial vault in humans is highly variable, clinically relevant, and heritable, yet its genetic architecture remains poorly understood. Here, we conduct a joint multi-ancestry and admixed multivariate genome-wide association study on 3D cranial vault shape extracted from magnetic resonance images of 6772 children from the ABCD study cohort yielding 30 genome-wide significant loci. Follow-up analyses indicate that these loci overlap with genomic risk loci for sagittal craniosynostosis, show elevated activity cranial neural crest cells, are enriched for processes related to skeletal development, and are shared with the face and brain. We present supporting evidence of regional localization for several of the identified genes based on expression patterns in the cranial vault bones of E15.5 mice. Overall, our study provides a comprehensive overview of the genetics underlying normal-range cranial vault shape and its relevance for understanding modern human craniofacial diversity and the etiology of congenital malformations.

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.

A complex genetic architecture underlies mandibular evolution in big mice from Gough Island

Some of the most compelling examples of morphological evolution come from island populations. Alterations in the size and shape of the mandible have been repeatedly observed in murid rodents following island colonization. Despite this pattern and the significance of the mandible for dietary adaptation, the genetic basis of island-mainland divergence in mandibular form remains uninvestigated. To fill this gap, we examined mandibular morphology in 609 F2s from a cross between Gough Island mice, the largest wild house mice on record, and mice from a mainland reference strain (WSB). Univariate genetic mapping identifies 3 quantitative trait loci (QTL) for relative length of the temporalis lever arm and 2 distinct QTL for relative condyle length, 2 traits expected to affect mandibular function that differ between Gough Island mice and WSB mice. Multivariate genetic mapping of coordinates from geometric morphometric analyses identifies 27 QTL contributing to overall mandibular shape. Quantitative trait loci show a complex mixture of modest, additive effects dispersed throughout the mandible, with landmarks including the coronoid process and the base of the ascending ramus frequently modulated by QTL. Additive effects of most shape quantitative trait loci do not align with island-mainland divergence, suggesting that directional selection played a limited role in the evolution of mandibular shape. In contrast, Gough Island mouse alleles at QTL for centroid size and QTL for jaw length increase these measures, suggesting selection led to larger mandibles, perhaps as a correlated response to the evolution of larger bodies.

Testing the accuracy of 3D automatic landmarking via genome-wide association studies

Various advances in 3D automatic phenotyping and landmark-based geometric morphometric methods have been made. While it is generally accepted that automatic landmarking compromises the capture of the biological variation, no studies have directly tested the actual impact of such landmarking approaches in analyses requiring a large number of specimens and for which the precision of phenotyping is crucial to extract an actual biological signal adequately. Here, we use a recently developed 3D atlas-based automatic landmarking method to test its accuracy in detecting QTLs associated with craniofacial development of the house mouse skull and lower jaws for a large number of specimens (circa 700) that were previously phenotyped via a semiautomatic landmarking method complemented with manual adjustment. We compare both landmarking methods with univariate and multivariate mapping of the skull and the lower jaws. We find that most significant SNPs and QTLs are not recovered based on the data derived from the automatic landmarking method. Our results thus confirm the notion that information is lost in the automated landmarking procedure although somewhat dependent on the analyzed structure. The automatic method seems to capture certain types of structures slightly better, such as lower jaws whose shape is almost entirely summarized by its outline and could be assimilated as a 2D flat object. By contrast, the more apparent 3D features exhibited by a structure such as the skull are not adequately captured by the automatic method. We conclude that using 3D atlas-based automatic landmarking methods requires careful consideration of the experimental question.

Relating multivariate shapes to genescapes using phenotype-biological process associations for craniofacial shape

Realistic mappings of genes to morphology are inherently multivariate on both sides of the equation. The importance of coordinated gene effects on morphological phenotypes is clear from the intertwining of gene actions in signaling pathways, gene regulatory networks, and developmental processes underlying the development of shape and size. Yet, current approaches tend to focus on identifying and localizing the effects of individual genes and rarely leverage the information content of high-dimensional phenotypes. Here, we explicitly model the joint effects of biologically coherent collections of genes on a multivariate trait – craniofacial shape – in a sample of n = 1145 mice from the Diversity Outbred (DO) experimental line. We use biological process Gene Ontology (GO) annotations to select skeletal and facial development gene sets and solve for the axis of shape variation that maximally covaries with gene set marker variation. We use our process-centered, multivariate genotype-phenotype (process MGP) approach to determine the overall contributions to craniofacial variation of genes involved in relevant processes and how variation in different processes corresponds to multivariate axes of shape variation. Further, we compare the directions of effect in phenotype space of mutations to the primary axis of shape variation associated with broader pathways within which they are thought to function. Finally, we leverage the relationship between mutational and pathway-level effects to predict phenotypic effects beyond craniofacial shape in specific mutants. We also introduce an online application that provides users the means to customize their own process-centered craniofacial shape analyses in the DO. The process-centered approach is generally applicable to any continuously varying phenotype and thus has wide-reaching implications for complex trait genetics.

MesoNet allows automated scaling and segmentation of mouse mesoscale cortical maps using machine learning

Understanding the basis of brain function requires knowledge of cortical operations over wide spatial scales and the quantitative analysis of brain activity in well-defined brain regions. Matching an anatomical atlas to brain functional data requires substantial labor and expertise. Here, we developed an automated machine learning-based registration and segmentation approach for quantitative analysis of mouse mesoscale cortical images. A deep learning model identifies nine cortical landmarks using only a single raw fluorescent image. Another fully convolutional network was adapted to delimit brain boundaries. This anatomical alignment approach was extended by adding three functional alignment approaches that use sensory maps or spatial-temporal activity motifs. We present this methodology as MesoNet, a robust and user-friendly analysis pipeline using pre-trained models to segment brain regions as defined in the Allen Mouse Brain Atlas. This Python-based toolbox can also be combined with existing methods to facilitate high-throughput data analysis.

Variation in Pleiotropic Hub Gene Expression Is Associated with Interspecific Differences in Head Shape and Eye Size in Drosophila

Revealing the mechanisms underlying the breathtaking morphological diversity observed in nature is a major challenge in Biology. It has been established that recurrent mutations in hotspot genes cause the repeated evolution of morphological traits, such as body pigmentation or the gain and loss of structures. To date, however, it remains elusive whether hotspot genes contribute to natural variation in the size and shape of organs. As natural variation in head morphology is pervasive in Drosophila, we studied the molecular and developmental basis of differences in compound eye size and head shape in two closely related Drosophila species. We show differences in the progression of retinal differentiation between species and we applied comparative transcriptomics and chromatin accessibility data to identify the GATA transcription factor Pannier (Pnr) as central factor associated with these differences. Although the genetic manipulation of Pnr affected multiple aspects of dorsal head development, the effect of natural variation is restricted to a subset of the phenotypic space. We present data suggesting that this developmental constraint is caused by the coevolution of expression of pnr and its cofactor u-shaped (ush). We propose that natural variation in expression or function of highly connected developmental regulators with pleiotropic functions is a major driver for morphological evolution and we discuss implications on gene regulatory network evolution. In comparison to previous findings, our data strongly suggest that evolutionary hotspots are not the only contributors to the repeated evolution of eye size and head shape in Drosophila.

Capturing the rapidly evolving study of adaptation

Research on the genomics of adaptation is rapidly changing. In the last few decades, progress in this area has been driven by methodological advances, not only in the way increasingly large amounts of molecular data are generated (e.g. with high-throughput sequencing), but also in the way these data are analysed. This includes a growing appreciation and quantitative treatment of covariation among units within the same data type (e.g. genes) or across data types (e.g. genes and phenotypes). The development and adoption of more and more integrative tools have resulted in richer and more interesting empirical work. This special issue - comprising methodological, empirical, and review papers - aims to capture a 'snapshot' of this rapidly evolving field. We discuss in particular three important themes in the study of adaptation: the genetic architecture of adaptive variation, protein-coding and regulatory changes, and parallel evolution. We highlight how more traditional key themes in the study of genetic architecture (e.g. the number of loci underlying adaptive traits and the distribution of their effects) are now being complemented by other factors (e.g. how patterns of linkage and number of loci interact to affect the ability to adapt). Similarly, apart from addressing the relative importance of protein-coding and regulatory changes, we now have the tools to look in-depth at specific types of regulatory variation to gain a clearer picture of regulatory networks. Finally, parallel evolution has always been central to the study of adaptation, but now we are often able to address the question of whether - and to what extent - parallelism at the organismal or phenotypic level is matched by parallelism at the genetic level. Perhaps most importantly, we can now determine what mechanisms are driving parallelism (or lack thereof) across levels of biological organization. All these recent methodological developments open up new directions for future studies of adaptive changes across traits, levels of biological organization, demographic contexts and time scales.

The comparative genomic landscape of adaptive radiation in crater lake cichlid fishes

Factors ranging from ecological opportunity to genome composition might explain why only some lineages form adaptive radiations. While being rare, particular systems can provide natural experiments within an identical ecological setting where species numbers and phenotypic divergence in two closely related lineages are notably different. We investigated one such natural experiment using two de novo assembled and 40 resequenced genomes and asked why two closely related Neotropical cichlid fish lineages, the Amphilophus citrinellus species complex (Midas cichlids; radiating) and Archocentrus centrarchus (Flyer cichlid; nonradiating), have resulted in such disparate evolutionary outcomes. Although both lineages inhabit many of the same Nicaraguan lakes, whole-genome inferred demography suggests that priority effects are not likely to be the cause of the dissimilarities. Also, genome-wide levels of selection, transposable element dynamics, gene family expansion, major chromosomal rearrangements and the number of genes under positive selection were not markedly different between the two lineages. To more finely investigate particular subsets of the genome that have undergone adaptive divergence in Midas cichlids, we also examined if there was evidence for 'molecular pre-adaptation' in regions identified by QTL mapping of repeatedly diverging adaptive traits. Although most of our analyses failed to pinpoint substantial genomic differences, we did identify functional categories containing many genes under positive selection that provide candidates for future studies on the propensity of Midas cichlids to radiate. Our results point to a disproportionate role of local, rather than genome-wide factors underlying the propensity for these cichlid fishes to adaptively radiate.

Chronic Cranial Window for Imaging Cortical Activity in Head-Fixed Mice

Chronic cranial window surgery is a critical procedure for in vivo imaging in neuroscience. Here, we describe our surgical protocol with several subtle improvements that increase the success rate significantly. The window allows high-quality imaging in head-fixed behaving mice within the first week after the surgical procedure and remains clear for months. We used this procedure to prepare mice for intrinsic signal imaging and two-photon imaging of layer 6 neurons in visual cortex.

For complete details on the use and execution of this protocol, please refer to Augustinaite and Kuhn (2020).

•

Long-term, high-quality chronic cranial window for imaging applications

•

Instructions for craniotomy and mounting of a cranial window with headplate

•

Instructions for tracer, dye, and/or virus injection during the surgery

•

Specific instructions for surgeries over primary visual cortex (V1) in mice

Sex as a Biological Variable in Preclinical Modeling of Blast-Related Traumatic Brain Injury

Approaches to furthering our understanding of the bioeffects, behavioral changes, and treatment options following exposure to blast are a worldwide priority. Of particular need is a more concerted effort to employ animal models to determine possible sex differences, which have been reported in the clinical literature. In this review, clinical and preclinical reports concerning blast injury effects are summarized in relation to sex as a biological variable (SABV). The review outlines approaches that explore the pertinent role of sex chromosomes and gonadal steroids for delineating sex as a biological independent variable. Next, underlying biological factors that need exploration for blast effects in light of SABV are outlined, including pituitary, autonomic, vascular, and inflammation factors that all have evidence as having important SABV relevance. A major second consideration for the study of SABV and preclinical blast effects is the notable lack of consistent model design—a wide range of devices have been employed with questionable relevance to real-life scenarios—as well as poor standardization for reporting of blast parameters. Hence, the review also provides current views regarding optimal design of shock tubes for approaching the problem of primary blast effects and sex differences and outlines a plan for the regularization of reporting. Standardization and clear description of blast parameters will provide greater comparability across models, as well as unify consensus for important sex difference bioeffects.

Facial shape and allometry quantitative trait locus intervals in the Diversity Outbred mouse are enriched for known skeletal and facial development genes

The biology of how faces are built and come to differ from one another is complex. Discovering normal variants that contribute to differences in facial morphology is one key to untangling this complexity, with important implications for medicine and evolutionary biology. This study maps quantitative trait loci (QTL) for skeletal facial shape using Diversity Outbred (DO) mice. The DO is a randomly outcrossed population with high heterozygosity that captures the allelic diversity of eight inbred mouse lines from three subspecies. The study uses a sample of 1147 DO animals (the largest sample yet employed for a shape QTL study in mouse), each characterized by 22 three-dimensional landmarks, 56,885 autosomal and X-chromosome markers, and sex and age classifiers. We identified 37 facial shape QTL across 20 shape principal components (PCs) using a mixed effects regression that accounts for kinship among observations. The QTL include some previously identified intervals as well as new regions that expand the list of potential targets for future experimental study. Three QTL characterized shape associations with size (allometry). Median support interval size was 3.5 Mb. Narrowing additional analysis to QTL for the five largest magnitude shape PCs, we found significant overrepresentation of genes with known roles in growth, skeletal and facial development, and sensory organ development. For most intervals, one or more of these genes lies within 0.25 Mb of the QTL's peak. QTL effect sizes were small, with none explaining more than 0.5% of facial shape variation. Thus, our results are consistent with a model of facial diversity that is influenced by key genes in skeletal and facial development and, simultaneously, is highly polygenic.

Sex differences in cued fear responses and parvalbumin cell density in the hippocampus following repetitive concussive brain injuries in C57BL/6J mice

There is strong evidence to suggest a link between repeated head trauma and cognitive and emotional disorders, and Repetitive concussive brain injuries (rCBI) may also be a risk factor for depression and anxiety disorders. Animal models of brain injury afford the opportunity for controlled study of the effects of injury on functional outcomes. In this study, male and cycling female C57BL/6J mice sustained rCBI (3x) at 24-hr intervals and were tested in a context and cued fear conditioning paradigm, open field (OF), elevated zero maze and tail suspension test. All mice with rCBI showed less freezing behavior than sham control mice during the fear conditioning context test. Injured male, but not female mice also froze less in response to the auditory cue (tone). Injured mice were hyperactive in an OF environment and spent more time in the open quadrants of the elevated zero maze, suggesting decreased anxiety, but there were no differences between injured mice and sham-controls in depressive-like activity on the tail suspension test. Pathologically, injured mice showed increased astrogliosis in the injured cortex and white matter tracts (optic tracts and corpus callosum). There were no changes in the number of parvalbumin-positive interneurons in the cortex or amygdala, but injured male mice had fewer parvalbumin-positive neurons in the hippocampus. Parvalbumin-reactive interneurons of the hippocampus have been previously demonstrated to be involved in hippocampal-cortical interactions required for memory consolidation, and it is possible memory changes in the fear-conditioning paradigm following rCBI are the result of more subtle imbalances in excitation and inhibition both within the amygdala and hippocampus, and between more widespread brain regions that are injured following a diffuse brain injury.

The genetic architecture of adaptation: convergence and pleiotropy in Heliconius wing pattern evolution

Unravelling the genetic basis of adaptive traits is a major challenge in evolutionary biology. Doing so informs our understanding of evolution towards an adaptive optimum, the distribution of locus effect sizes, and the influence of genetic architecture on the evolvability of a trait. In the Müllerian co-mimics Heliconius melpomene and Heliconius erato some Mendelian loci affecting mimicry shifts are well known. However, several phenotypes in H. melpomene remain to be mapped, and the quantitative genetics of colour pattern variation has rarely been analysed. Here we use quantitative trait loci (QTL) analyses of crosses between H. melpomene races from Peru and Suriname to map, for the first time, the control of the broken band phenotype to WntA and identify a ~100 kb region controlling this variation. Additionally, we map variation in basal forewing red-orange pigmentation to a locus centred around the gene ventral veins lacking (vvl). The locus also appears to affect medial band shape variation as it was previously known to do in H. erato. This adds to the list of homologous regions controlling convergent phenotypes between these two species. Finally we show that Heliconius wing-patterning genes are strikingly pleiotropic among wing pattern traits. Our results demonstrate how genetic architecture can shape, aid and constrain adaptive evolution.

PIN1 is a new therapeutic target of craniosynostosis

Gain-of-function mutations in fibroblast growth factor receptors (FGFRs) cause congenital skeletal anomalies, including craniosynostosis (CS), which is characterized by the premature closure of craniofacial sutures. Apert syndrome (AS) is one of the severest forms of CS, and the only treatment is surgical expansion of prematurely fused sutures in infants. Previously, we demonstrated that the prolyl isomerase peptidyl-prolyl cis-trans isomerase interacting 1 (PIN1) plays a critical role in mediating FGFR signaling and that Pin1+/- mice exhibit delayed closure of cranial sutures. In this study, using both genetic and pharmacological approaches, we tested whether PIN1 modulation could be used as a therapeutic regimen against AS. In the genetic approach, we crossbred Fgfr2S252W/+, a mouse model of AS, and Pin1+/- mice. Downregulation of Pin1 gene dosage attenuated premature cranial suture closure and other phenotypes of AS in Fgfr2S252W/+ mutant mice. In the pharmacological approach, we intraperitoneally administered juglone, a PIN1 enzyme inhibitor, to pregnant Fgfr2S252W/+ mutant mice and found that this treatment successfully interrupted fetal development of AS phenotypes. Primary cultured osteoblasts from Fgfr2S252W/+ mutant mice expressed high levels of FGFR2 downstream target genes, but this phenotype was attenuated by PIN1 inhibition. Post-translational stabilization and activation of Runt-related transcription factor 2 (RUNX2) in Fgfr2S252W/+ osteoblasts were also attenuated by PIN1 inhibition. Based on these observations, we conclude that PIN1 enzyme activity is important for FGFR2-induced RUNX2 activation and craniofacial suture morphogenesis. Moreover, these findings highlight that juglone or other PIN1 inhibitors represent viable alternatives to surgical intervention for treatment of CS and other hyperostotic diseases.

Epistasis regulates the developmental stability of the mouse craniofacial shape

Fluctuating asymmetry is a classic concept linked to organismal development. It has traditionally been used as a measure of developmental instability, which is the inability of an organism to buffer environmental fluctuations during development. Developmental stability has a genetic component that influences the final phenotype of the organism and can lead to congenital disorders. According to alternative hypotheses, this genetic component might be either the result of additive genetic effects or a by-product of developmental gene networks. Here we present a genome-wide association study of the genetic architecture of fluctuating asymmetry of the skull shape in mice. Geometric morphometric methods were applied to quantify fluctuating asymmetry: we estimated fluctuating asymmetry as Mahalanobis distances to the mean asymmetry, correcting first for genetic directional asymmetry. We applied the marginal epistasis test to study epistasis among genomic regions. Results showed no evidence of additive effects but several interacting regions significantly associated with fluctuating asymmetry. Among the candidate genes overlapping these interacting regions we found an over-representation of genes involved in craniofacial development. A gene network is likely to be associated with skull developmental stability, and genes originally described as buffering genes (e.g., Hspa2) might occupy central positions within these networks, where regulatory elements may also play an important role. Our results constitute an important step in the exploration of the molecular roots of developmental stability and the first empirical evidence about its genetic architecture.

Mouse Skull Mean Shape and Shape Robustness Rely on Different Genetic Architectures and Different Loci

The genetic architecture of skull shape has been extensively studied in mice and the results suggest a highly polygenic and additive basis. In contrast few studies have explored the genetic basis of the skull variability. Canalization and developmental stability are the two components of phenotypic robustness. They have been proposed to be emergent properties of the genetic networks underlying the development of the trait itself, but this hypothesis has been rarely tested empirically. Here we use outbred mice to investigate the genetic architecture of canalization of the skull shape by implementing a genome-wide marginal epistatic test on 3D geometric morphometric data. The same data set had been used previously to explore the genetic architecture of the skull mean shape and its developmental stability. Here, we address two questions: (1) Are changes in mean shape and changes in shape variance associated with the same genomic regions? and (2) Do canalization and developmental stability rely on the same loci and genetic architecture and do they involve the same patterns of shape variation? We found that unlike skull mean shape, among-individual shape variance and fluctuating asymmetry (FA) show a total lack of additive effects. They are both associated with complex networks of epistatic interactions involving many genes (protein-coding and regulatory elements). Remarkably, none of the genomic loci affecting mean shape contribute these networks despite their enrichment for genes involved in craniofacial variation and diseases. We also found that the patterns of shape FA and individual variation are largely similar and rely on similar multilocus epistatic genetic networks, suggesting that the processes channeling variation within and among individuals are largely common. However, the loci involved in these two networks are completely different. This in turn underlines the difference in the origin of the variation at these two levels, and points at buffering processes that may be specific to each level.

Quantitative morphometric analysis of adult teleost fish by X-ray computed tomography

Vertebrate models provide indispensable paradigms to study development and disease. Their analysis requires a quantitative morphometric study of the body, organs and tissues. This is often impeded by pigmentation and sample size. X-ray micro-computed tomography (micro-CT) allows high-resolution volumetric tissue analysis, largely independent of sample size and transparency to visual light. Importantly, micro-CT data are inherently quantitative. We report a complete pipeline of high-throughput 3D data acquisition and image analysis, including tissue preparation and contrast enhancement for micro-CT imaging down to cellular resolution, automated data processing and organ or tissue segmentation that is applicable to comparative 3D morphometrics of small vertebrates. Applied to medaka fish, we first create an annotated anatomical atlas of the entire body, including inner organs as a quantitative morphological description of an adult individual. This atlas serves as a reference model for comparative studies. Using isogenic medaka strains we show that comparative 3D morphometrics of individuals permits identification of quantitative strain-specific traits. Thus, our pipeline enables high resolution morphological analysis as a basis for genotype-phenotype association studies of complex genetic traits in vertebrates.

Mapping genetic variants for cranial vault shape in humans

The shape of the cranial vault, a region comprising interlocking flat bones surrounding the cerebral cortex, varies considerably in humans. Strongly influenced by brain size and shape, cranial vault morphology has both clinical and evolutionary relevance. However, little is known about the genetic basis of normal vault shape in humans. We performed a genome-wide association study (GWAS) on three vault measures (maximum cranial width [MCW], maximum cranial length [MCL], and cephalic index [CI]) in a sample of 4419 healthy individuals of European ancestry. All measures were adjusted by sex, age, and body size, then tested for association with genetic variants spanning the genome. GWAS results for the two cohorts were combined via meta-analysis. Significant associations were observed at two loci: 15p11.2 (lead SNP rs2924767, p = 2.107 × 10⁻⁸) for MCW and 17q11.2 (lead SNP rs72841279, p = 5.29 × 10⁻⁹) for MCL. Additionally, 32 suggestive loci (p < 5x10^-6) were observed. Several candidate genes were located in these loci, such as NLK, MEF2A, SOX9 and SOX11. Genome-wide linkage analysis of cranial vault shape in mice (N = 433) was performed to follow-up the associated candidate loci identified in the human GWAS. Two loci, 17q11.2 (c11.loc44 in mice) and 17q25.1 (c11.loc74 in mice), associated with cranial vault size in humans, were also linked with cranial vault size in mice (LOD scores: 3.37 and 3.79 respectively). These results provide further insight into genetic pathways and mechanisms underlying normal variation in human craniofacial morphology.

Genome-wide mapping of global-to-local genetic effects on human facial shape

Genome-wide association scans of complex multipartite traits like the human face typically use preselected phenotypic measures. Here we report a data-driven approach to phenotyping facial shape at multiple levels of organization, allowing for an open-ended description of facial variation while preserving statistical power. In a sample of 2,329 persons of European ancestry, we identified 38 loci, 15 of which replicated in an independent European sample (n = 1,719). Four loci were completely new. For the others, additional support (n = 9) or pleiotropic effects (n = 2) were found in the literature, but the results reported here were further refined. All 15 replicated loci highlighted distinctive patterns of global-to-local genetic effects on facial shape and showed enrichment for active chromatin elements in human cranial neural crest cells, suggesting an early developmental origin of the facial variation captured. These results have implications for studies of facial genetics and other complex morphological traits.

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.

Genetic mapping of molar size relations identifies inhibitory locus for third molars in mice

Molar size in Mammals shows considerable disparity and exhibits variation similar to that predicted by the Inhibitory Cascade model. The importance of such developmental systems in favoring evolutionary trajectories is also underlined by the fact that this model can predict macroevolutionary patterns. Using backcross mice, we mapped QTL for molar sizes controlling for their sequential development. Genetic controls for upper and lower molars appear somewhat similar, and regions containing genes implied in dental defects drive this variation. We mapped three relationship QTLs (rQTL) modifying the control of the mesial molars on the focal third molar. These regions overlap Shh, Sostdc1, and Fst genes, which have pervasive roles in development and should be buffered against new variation. It has theoretically been shown that rQTL produces new variation channeled in the direction of adaptive changes. Our results provide evidence that evolutionary/disease patterns of tooth size variation could result from such a non-random generating process.

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.

A Genome-Wide Association Study Identifies Multiple Regions Associated with Head Size in Catfish

Skull morphology is fundamental to evolution and the biological adaptation of species to their environments. With aquaculture fish species, head size is also important for economic reasons because it has a direct impact on fillet yield. However, little is known about the underlying genetic basis of head size. Catfish is the primary aquaculture species in the United States. In this study, we performed a genome-wide association study using the catfish 250K SNP array with backcross hybrid catfish to map the QTL for head size (head length, head width, and head depth). One significantly associated region on linkage group (LG) 7 was identified for head length. In addition, LGs 7, 9, and 16 contain suggestively associated regions for head length. For head width, significantly associated regions were found on LG9, and additional suggestively associated regions were identified on LGs 5 and 7. No region was found associated with head depth. Head size genetic loci were mapped in catfish to genomic regions with candidate genes involved in bone development. Comparative analysis indicated that homologs of several candidate genes are also involved in skull morphology in various other species ranging from amphibian to mammalian species, suggesting possible evolutionary conservation of those genes in the control of skull morphologies.

Genome-Wide Association Study Reveals Multiple Loci Influencing Normal Human Facial Morphology

Numerous lines of evidence point to a genetic basis for facial morphology in humans, yet little is known about how specific genetic variants relate to the phenotypic expression of many common facial features. We conducted genome-wide association meta-analyses of 20 quantitative facial measurements derived from the 3D surface images of 3118 healthy individuals of European ancestry belonging to two US cohorts. Analyses were performed on just under one million genotyped SNPs (Illumina OmniExpress+Exome v1.2 array) imputed to the 1000 Genomes reference panel (Phase 3). We observed genome-wide significant associations (p < 5 x 10⁻⁸) for cranial base width at 14q21.1 and 20q12, intercanthal width at 1p13.3 and Xq13.2, nasal width at 20p11.22, nasal ala length at 14q11.2, and upper facial depth at 11q22.1. Several genes in the associated regions are known to play roles in craniofacial development or in syndromes affecting the face: MAFB, PAX9, MIPOL1, ALX3, HDAC8, and PAX1. We also tested genotype-phenotype associations reported in two previous genome-wide studies and found evidence of replication for nasal ala length and SNPs in CACNA2D3 and PRDM16. These results provide further evidence that common variants in regions harboring genes of known craniofacial function contribute to normal variation in human facial features. Improved understanding of the genes associated with facial morphology in healthy individuals can provide insights into the pathways and mechanisms controlling normal and abnormal facial morphogenesis.

There is a great deal of evidence that genes influence facial appearance. This is perhaps most apparent when we look at our own families, since we are more likely to share facial features in common with our close relatives than with unrelated individuals. Nevertheless, little is known about how variation in specific regions of the genome relates to the kinds of distinguishing facial characteristics that give us our unique identities, e.g., the size and shape of our nose or how far apart our eyes are spaced. In this paper, we investigate this question by examining the association between genetic variants across the whole genome and a set of measurements designed to capture key aspects of facial form. We found evidence of genetic associations involving measures of eye, nose, and facial breadth. In several cases, implicated regions contained genes known to play roles in embryonic face formation or in syndromes in which the face is affected. Our ability to connect specific genetic variants to ubiquitous facial traits can inform our understanding of normal and abnormal craniofacial development, provide potential predictive models of evolutionary changes in human facial features, and improve our ability to create forensic facial reconstructions from DNA.

Natural history collections as windows on evolutionary processes

Natural history collections provide an immense record of biodiversity on Earth. These repositories have traditionally been used to address fundamental questions in biogeography, systematics and conservation. However, they also hold the potential for studying evolution directly. While some of the best direct observations of evolution have come from long-term field studies or from experimental studies in the laboratory, natural history collections are providing new insights into evolutionary change in natural populations. By comparing phenotypic and genotypic changes in populations through time, natural history collections provide a window into evolutionary processes. Recent studies utilizing this approach have revealed some dramatic instances of phenotypic change over short timescales in response to presumably strong selective pressures. In some instances, evolutionary change can be paired with environmental change, providing a context for potential selective forces. Moreover, in a few cases, the genetic basis of phenotypic change is well understood, allowing for insight into adaptive change at multiple levels. These kinds of studies open the door to a wide range of previously intractable questions by enabling the study of evolution through time, analogous to experimental studies in the laboratory, but amenable to a diversity of species over longer timescales in natural populations.

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 Genetic Basis of Baculum Size and Shape Variation in Mice

The rapid divergence of male genitalia is a preeminent evolutionary pattern. This rapid divergence is especially striking in the baculum, a bone that occurs in the penis of many mammalian species. Closely related species often display diverse baculum morphology where no other morphological differences can be discerned. While this fundamental pattern of evolution has been appreciated at the level of gross morphology, nearly nothing is known about the genetic basis of size and shape divergence. Quantifying the genetic basis of baculum size and shape variation has been difficult because these structures generally lack obvious landmarks, so comparing them in three dimensions is not straightforward. Here, we develop a novel morphometric approach to quantify size and shape variation from three-dimensional micro-CT scans taken from 369 bacula, representing 75 distinct strains of the BXD family of mice. We identify two quantitative trait loci (QTL) that explain ∼50% of the variance in baculum size, and a third QTL that explains more than 20% of the variance in shape. Together, our study demonstrates that baculum morphology may diverge relatively easily, with mutations at a few loci of large effect that independently modulate size and shape. Based on a combination of bioinformatic investigations and new data on RNA expression, we prioritized these QTL to 16 candidate genes, which have hypothesized roles in bone morphogenesis and may enable future genetic manipulation of baculum morphology.

Postnatal Development of the Craniofacial Skeleton in Male C57BL/6J Mice

C57BL/6J is one of the most commonly used inbred mouse strains in biomedical research, including studies of craniofacial development and teratogenic studies of craniofacial malformation. The current study quantitatively assessed the development of the skull in male C57BL/6J mice by using high-resolution 3D imaging of 55 landmarks from 48 male mice over 10 developmental time points from postnatal day 0 to 90. The growth of the skull plateaued at approximately postnatal day 60, and the shape of the skull did not change markedly thereafter. The amount of asymmetry in the craniofacial skeleton seemed to peak at birth, but considerable variation persisted in all age groups. For C57BL/6J male mice, postnatal day 60 is the earliest time point at which the skull achieves its adult shape and proportions. In addition, C57BL/6J male mice appear to have an inherent susceptibility to craniofacial malformation.

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.

Performance of single and multi-atlas based automated landmarking methods compared to expert annotations in volumetric microCT datasets of mouse mandibles

BACKGROUND

Here we present an application of advanced registration and atlas building framework DRAMMS to the automated annotation of mouse mandibles through a series of tests using single and multi-atlas segmentation paradigms and compare the outcomes to the current gold standard, manual annotation.

RESULTS

Our results showed multi-atlas annotation procedure yields landmark precisions within the human observer error range. The mean shape estimates from gold standard and multi-atlas annotation procedure were statistically indistinguishable for both Euclidean Distance Matrix Analysis (mean form matrix) and Generalized Procrustes Analysis (Goodall F-test). Further research needs to be done to validate the consistency of variance-covariance matrix estimates from both methods with larger sample sizes.

CONCLUSION

Multi-atlas annotation procedure shows promise as a framework to facilitate truly high-throughput phenomic analyses by channeling investigators efforts to annotate only a small portion of their datasets.

Mapping of Craniofacial Traits in Outbred Mice Identifies Major Developmental Genes Involved in Shape Determination

The vertebrate cranium is a prime example of the high evolvability of complex traits. While evidence of genes and developmental pathways underlying craniofacial shape determination is accumulating, we are still far from understanding how such variation at the genetic level is translated into craniofacial shape variation. Here we used 3D geometric morphometrics to map genes involved in shape determination in a population of outbred mice (Carworth Farms White, or CFW). We defined shape traits via principal component analysis of 3D skull and mandible measurements. We mapped genetic loci associated with shape traits at ~80,000 candidate single nucleotide polymorphisms in ~700 male mice. We found that craniofacial shape and size are highly heritable, polygenic traits. Despite the polygenic nature of the traits, we identified 17 loci that explain variation in skull shape, and 8 loci associated with variation in mandible shape. Together, the associated variants account for 11.4% of skull and 4.4% of mandible shape variation, however, the total additive genetic variance associated with phenotypic variation was estimated in ~45%. Candidate genes within the associated loci have known roles in craniofacial development; this includes 6 transcription factors and several regulators of bone developmental pathways. One gene, Mn1, has an unusually large effect on shape variation in our study. A knockout of this gene was previously shown to affect negatively the development of membranous bones of the cranial skeleton, and evolutionary analysis shows that the gene has arisen at the base of the bony vertebrates (Eutelostomi), where the ossified head first appeared. Therefore, Mn1 emerges as a key gene for both skull formation and within-population shape variation. Our study shows that it is possible to identify important developmental genes through genome-wide mapping of high-dimensional shape features in an outbred population.

Formation of the face, mandible, and skull is determined in part by genetic factors, but the relationship between genetic variation and craniofacial development is not well understood. We demonstrate how recent advances in mouse genomics and statistical methods can be used to identify genes involved in craniofacial development. We use outbred mice together with a dense panel of genetic markers to identify genetic loci affecting craniofacial shape. Some of the loci we identify are also known from past studies to contribute to craniofacial development and bone formation. For example, the top candidate gene identified in this study, Mn1, is a gene that appeared at a time when animals started to form bony skulls, suggesting that it may be a key gene in this evolutionary innovation. This further suggests that Mn1 and other genes involved in head formation are also responsible for more fine-grained regulation of its shape. Our results confirm that the outbred mouse population used in this study is suitable to identify single genetic factors even under conditions where many genes cooperate to generate a complex phenotype.