FaceBase 3: analytical tools and FAIR resources for craniofacial and dental research

Tendon-associated gene expression precedes osteogenesis in mid-palatal suture establishment

Orthodontic maxillary expansion relies on intrinsic mid-palatal suture mechanobiology to induce guided osteogenesis, yet establishment of the mid-palatal suture within the continuous secondary palate and causes of maxillary insufficiency remain poorly understood. In contrast, advances in cranial suture research hold promise to improve surgical repair of prematurely fused cranial sutures in craniosynostosis to potentially restore the obliterated signaling environment and ensure continual success of the intervention. We hypothesized that mid-palatal suture establishment is governed by shared principles with calvarial sutures and involves functional linkage between expanding primary ossification centres with the midline mesenchyme. We characterized establishment of the mid-palatal suture from late embryonic to early postnatal timepoints. Suture establishment was visualized using histological techniques and multimodal transcriptomics. We identified that mid-palatal suture formation depends on a spatiotemporally controlled signalling milieu in which tendon-associated genes play a significant role. We mapped relationships between extracellular matrix-encoding gene expression, tenocyte markers, and novel suture patency candidate genes. We identified similar expression patterns in FaceBase-deposited scRNA-seq datasets from cranial sutures. These findings demonstrate shared biological principles for suture establishment, providing further avenues for future development and understanding of maxillofacial interventions.

A common cis-regulatory variant impacts normal-range and disease-associated human facial shape through regulation of PKDCC during chondrogenesis

Genome-wide association studies (GWAS) identified thousands of genetic variants linked to phenotypic traits and disease risk. However, mechanistic understanding of how GWAS variants influence complex morphological traits and can, in certain cases, simultaneously confer normal-range phenotypic variation and disease predisposition, is still largely lacking. Here, we focus on rs6740960, a single nucleotide polymorphism (SNP) at the 2p21 locus, which in GWAS studies has been associated both with normal-range variation in jaw shape and with an increased risk of non-syndromic orofacial clefting. Using in vitro derived embryonic cell types relevant for human facial morphogenesis, we show that this SNP resides in an enhancer that regulates chondrocytic expression of PKDCC - a gene encoding a tyrosine kinase involved in chondrogenesis and skeletal development. In agreement, we demonstrate that the rs6740960 SNP is sufficient to confer chondrocyte-specific differences in PKDCC expression. By deploying dense landmark morphometric analysis of skull elements in mice, we show that changes in Pkdcc dosage are associated with quantitative changes in the maxilla, mandible, and palatine bone shape that are concordant with the facial phenotypes and disease predisposition seen in humans. We further demonstrate that the frequency of the rs6740960 variant strongly deviated among different human populations, and that the activity of its cognate enhancer diverged in hominids. Our study provides a mechanistic explanation of how a common SNP can mediate normal-range and disease-associated morphological variation, with implications for the evolution of human facial features.

Dynamic enhancer landscapes in human craniofacial development

The genetic basis of human facial variation and craniofacial birth defects remains poorly understood. Distant-acting transcriptional enhancers control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development. However, a lack of accurate maps of the genomic locations and cell type-resolved activities of craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combine histone modification, chromatin accessibility, and gene expression profiling of human craniofacial development with single-cell analyses of the developing mouse face to define the regulatory landscape of facial development at tissue- and single cell-resolution. We provide temporal activity profiles for 14,000 human developmental craniofacial enhancers. We find that 56% of human craniofacial enhancers share chromatin accessibility in the mouse and we provide cell population- and embryonic stage-resolved predictions of their in vivo activity. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

Single-cell multiomics decodes regulatory programs for mouse secondary palate development

Perturbations in gene regulation during palatogenesis can lead to cleft palate, which is among the most common congenital birth defects. Here, we perform single-cell multiome sequencing and profile chromatin accessibility and gene expression simultaneously within the same cells (n = 36,154) isolated from mouse secondary palate across embryonic days (E) 12.5, E13.5, E14.0, and E14.5. We construct five trajectories representing continuous differentiation of cranial neural crest-derived multipotent cells into distinct lineages. By linking open chromatin signals to gene expression changes, we characterize the underlying lineage-determining transcription factors. In silico perturbation analysis identifies transcription factors SHOX2 and MEOX2 as important regulators of the development of the anterior and posterior palate, respectively. In conclusion, our study charts epigenetic and transcriptional dynamics in palatogenesis, serving as a valuable resource for further cleft palate research.

Proteomic Analyses Discern the Developmental Inclusion of Albumin in Pig Enamel: A New Model for Human Enamel Hypomineralization

Excess albumin in enamel is a characteristic of the prevalent developmental dental defect known as chalky teeth or molar hypomineralization (MH). This study uses proteomic analyses of pig teeth to discern between developmental origin and post-eruptive contamination and to assess the similarity to hypomineralized human enamel. Here, the objective is to address the urgent need for an animal model to uncover the etiology of MH and to improve treatment. Porcine enamel is chalky and soft at eruption; yet, it hardens quickly to form a hard surface and then resembles human teeth with demarcated enamel opacities. Proteomic analyses of enamel from erupted teeth, serum, and saliva from pigs aged 4 (n = 3) and 8 weeks (n = 2) and human (n = 4) molars with demarcated enamel opacities show alpha-fetoprotein (AFP). AFP expression is limited to pre- and perinatal development and its presence in enamel indicates pre- or perinatal inclusion. In contrast, albumin is expressed after birth, indicating postnatal inclusion into enamel. Peptides were extracted from enamel and analyzed by nano-liquid chromatography-tandem mass spectrometry (nanoLC-MS/MS) after tryptic digestion. The mean total protein number was 337 in the enamel of all teeth with 13 different unique tryptic peptides of porcine AFP in all enamel samples but none in saliva samples. Similarities in the composition, micro-hardness, and microstructure underscore the usefulness of the porcine model to uncover the MH etiology, cellular mechanisms of albumin inclusion, and treatment for demarcated opacities.

Integrative analysis of transcriptome dynamics during human craniofacial development identifies candidate disease genes

Craniofacial disorders arise in early pregnancy and are one of the most common congenital defects. To fully understand how craniofacial disorders arise, it is essential to characterize gene expression during the patterning of the craniofacial region. To address this, we performed bulk and single-cell RNA-seq on human craniofacial tissue from 4-8 weeks post conception. Comparisons to dozens of other human tissues revealed 239 genes most strongly expressed during craniofacial development. Craniofacial-biased developmental enhancers were enriched +/- 400 kb surrounding these craniofacial-biased genes. Gene co-expression analysis revealed that regulatory hubs are enriched for known disease causing genes and are resistant to mutation in the normal healthy population. Combining transcriptomic and epigenomic data we identified 539 genes likely to contribute to craniofacial disorders. While most have not been previously implicated in craniofacial disorders, we demonstrate this set of genes has increased levels of de novo mutations in orofacial clefting patients warranting further study.

Cell Type- and Tissue-specific Enhancers in Craniofacial Development

The genetic basis of craniofacial birth defects and general variation in human facial shape remains poorly understood. Distant-acting transcriptional enhancers are a major category of non-coding genome function and have been shown to control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development1-3. However, a lack of accurate maps of the genomic location and cell type-specific in vivo activities of all craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combined histone modification and chromatin accessibility profiling from different stages of human craniofacial development with single-cell analyses of the developing mouse face to create a comprehensive catalogue of the regulatory landscape of facial development at tissue- and single cell-resolution. In total, we identified approximately 14,000 enhancers across seven developmental stages from weeks 4 through 8 of human embryonic face development. We used transgenic mouse reporter assays to determine the in vivo activity patterns of human face enhancers predicted from these data. Across 16 in vivo validated human enhancers, we observed a rich diversity of craniofacial subregions in which these enhancers are active in vivo. To annotate the cell type specificities of human-mouse conserved enhancers, we performed single-cell RNA-seq and single-nucleus ATAC-seq of mouse craniofacial tissues from embryonic days e11.5 to e15.5. By integrating these data across species, we find that the majority (56%) of human craniofacial enhancers are functionally conserved in mice, providing cell type- and embryonic stage-resolved predictions of their in vivo activity profiles. Using retrospective analysis of known craniofacial enhancers in combination with single cell-resolved transgenic reporter assays, we demonstrate the utility of these data for predicting the in vivo cell type specificity of enhancers. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

Comparing ant morphology measurements from microscope and online AntWeb.org 2D z-stacked images

Unprecedented technological advances in digitization and the steadily expanding open-access digital repositories are yielding new opportunities to quickly and efficiently measure morphological traits without transportation and advanced/expensive microscope machinery. A prime example is the AntWeb.org database, which allows researchers from all over the world to study taxonomic, ecological, or evolutionary questions on the same ant specimens with ease. However, the reproducibility and reliability of morphometric data deduced from AntWeb compared to traditional microscope measurements has not yet been tested. Here, we compared 12 morphological traits of 46 Temnothorax ant specimens measured either directly by stereomicroscope on physical specimens or via the widely used open-access software tpsDig utilizing AntWeb digital images. We employed a complex statistical framework to test several aspects of reproducibility and reliability between the methods. We estimated (i) the agreement between the measurement methods and (ii) the trait value dependence of the agreement, then (iii) compared the coefficients of variation produced by the different methods, and finally, (iv) tested for systematic bias between the methods in a mixed modeling-based statistical framework. The stereomicroscope measurements were extremely precise. Our comparisons showed that agreement between the two methods was exceptionally high, without trait value dependence. Furthermore, the coefficients of variation did not differ between the methods. However, we found systematic bias in eight traits: apart from one trait where software measurements overestimated the microscopic measurements, the former underestimated the latter. Our results shed light on the fact that relying solely on the level of agreement between methods can be highly misleading. In our case, even though the software measurements predicted microscope measurements very well, replacing traditional microscope measurements with software measurements, and especially mixing data collected by the different methods, might result in erroneous conclusions. We provide guidance on the best way to utilize virtual specimens (2D z-stacked images) as a source of morphometric data, emphasizing the method's limitations in certain fields and applications.

Micro-CT analysis reveals the changes in bone mineral density in zebrafish craniofacial skeleton with age

Bone has multiple functions in animals, such as supporting the body for mobility. The zebrafish skeleton is composed of craniofacial and axial skeletons. It shares a physiological curvature and consists of a similar number of vertebrae as humans. Bone degeneration and malformations have been widely studied in zebrafish as human disease models. High-resolution imaging and different bone properties such as density and volume can be obtained using micro-computed tomography (micro-CT). This study aimed to understand the possible changes in the structure and bone mineral density (BMD) of the vertebrae and craniofacial skeleton with age (4, 12 and 24 months post fertilisation [mpf]) in zebrafish. Our data showed that the BMD in the vertebrae and specific craniofacial skeleton (mandibular arch, ceratohyal and ethmoid plate) of 12 and 24 mpf fish were higher than that of the 4 mpf fish. In addition, we found the age-dependent increase in BMD was not ubiquitously observed in facial bones, and such differences were not correlated with bone type. In summary, such additional information on the craniofacial skeleton could help in understanding bone development throughout the lifespan of zebrafish.

Spatiotemporal MicroRNA-Gene Expression Network Related to Orofacial Clefts

Craniofacial structures change dynamically in morphology during development through the coordinated regulation of various cellular molecules. However, it remains unclear how these complex mechanisms are regulated in a spatiotemporal manner. Here we applied natural cubic splines to model gene and microRNA (miRNA) expression from embryonic day (E) 10.5 to E14.5 in the proximal and distal regions of the maxillary processes to identify spatiotemporal patterns of gene and miRNA expression, followed by constructing corresponding regulatory networks. Three major groups of differentially expressed genes (DEGs) were identified, including 3,927 temporal, 314 spatial, and 494 spatiotemporal DEGs. Unsupervised clustering further resolved these spatiotemporal DEGs into 8 clusters with distinct expression patterns. Interestingly, we found 2 clusters of differentially expressed miRNAs: 1 had 80 miRNAs monotonically decreasing and the other had 97 increasing across developmental stages. To evaluate the phenotypic relevance of these DEGs during craniofacial development, we integrated data from the CleftGeneDB database and constructed the regulatory networks of genes related to orofacial clefts. Our analysis revealed 2 hub miRNAs, mmu-miR-325-3p and mmu-miR-384-5p, that repressed cleft-related genes Adamts3, Runx2, Fgfr2, Acvr1, and Edn2, while their expression increased over time. On the contrary, 2 hub miRNAs, mmu-miR-218-5p and mmu-miR-338-5p, repressed cleft-related genes Pbx2, Ermp1, Snai1, Tbx2, and Bmi1, while their expression decreased over time. Our experiments indicated that these miRNA mimics significantly inhibited cell proliferation in mouse embryonic palatal mesenchymal (MEPM) cells and O9-1 cells through the regulation of genes associated with cleft palate and validated the role of our regulatory networks in orofacial clefts. To facilitate interactive exploration of these data, we developed a user-friendly web tool to visualize the gene and miRNA expression patterns across developmental stages, as well as the regulatory networks (https://fyan.shinyapps.io/facebase_shiny/). Taken together, our results provide a valuable resource that serves as a reference map for future research in craniofacial development.

FaceBase: A Community-Driven Hub for Data-Intensive Research

The FaceBase Consortium, funded by the National Institute of Dental and Craniofacial Research of the National Institutes of Health, was established in 2009 with the recognition that dental and craniofacial research are increasingly data-intensive disciplines. Data sharing is critical for the validation and reproducibility of results as well as to enable reuse of data. In service of these goals, data ought to be FAIR: Findable, Accessible, Interoperable, and Reusable. The FaceBase data repository and educational resources exemplify the FAIR principles and support a broad user community including researchers in craniofacial development, molecular genetics, and genomics. FaceBase demonstrates that a model in which researchers "self-curate" their data can be successful and scalable. We present the results of the first 2.5 y of FaceBase's operations as an open community and summarize the data sets published during this period. We then describe a research highlight from work on the identification of regulatory networks and noncoding RNAs involved in cleft lip with/without cleft palate that both used and in turn contributed new findings to publicly available FaceBase resources. Collectively, FaceBase serves as a dynamic and continuously evolving resource to facilitate data-intensive research, enhance data reproducibility, and perform deep phenotyping across multiple species in dental and craniofacial research.

A Roadmap for the Human Oral and Craniofacial Cell Atlas

Oral and craniofacial tissues are uniquely adapted for continuous and intricate functioning, including breathing, feeding, and communication. To achieve these vital processes, this complex is supported by incredible tissue diversity, variously composed of epithelia, vessels, cartilage, bone, teeth, ligaments, and muscles, as well as mesenchymal, adipose, and peripheral nervous tissue. Recent single cell and spatial multiomics assays-specifically, genomics, epigenomics, transcriptomics, proteomics, and metabolomics-have annotated known and new cell types and cell states in human tissues and animal models, but these concepts remain limitedly explored in the human postnatal oral and craniofacial complex. Here, we highlight the collaborative and coordinated efforts of the newly established Oral and Craniofacial Bionetwork as part of the Human Cell Atlas, which aims to leverage single cell and spatial multiomics approaches to first understand the cellular and molecular makeup of human oral and craniofacial tissues in health and to then address common and rare diseases. These powerful assays have already revealed the cell types that support oral tissues, and they will unravel cell types and molecular networks utilized across development, maintenance, and aging as well as those affected in diseases of the craniofacial complex. This level of integration and cell annotation with partner laboratories across the globe will be critical for understanding how multiple variables, such as age, sex, race, and ancestry, influence these oral and craniofacial niches. Here, we 1) highlight these recent collaborative efforts to employ new single cell and spatial approaches to resolve our collective biology at a higher resolution in health and disease, 2) discuss the vision behind the Oral and Craniofacial Bionetwork, 3) outline the stakeholders who contribute to and will benefit from this network, and 4) outline directions for creating the first Human Oral and Craniofacial Cell Atlas.

Sharing Begins at Home: How Continuous and Ubiquitous FAIRness Can Enhance Research Productivity and Data Reuse

The broad sharing of research data is widely viewed as critical for the speed, quality, accessibility, and integrity of science. Despite increasing efforts to encourage data sharing, both the quality of shared data and the frequency of data reuse remain stubbornly low. We argue here that a significant reason for this unfortunate state of affairs is that the organization of research results in the findable, accessible, interoperable, and reusable (FAIR) form required for reuse is too often deferred to the end of a research project when preparing publications-by which time essential details are no longer accessible. Thus, we propose an approach to research informatics in which FAIR principles are applied continuously, from the inception of a research project and ubiquitously, to every data asset produced by experiment or computation. We suggest that this seemingly challenging task can be made feasible by the adoption of simple tools, such as lightweight identifiers (to ensure that every data asset is findable), packaging methods (to facilitate understanding of data contents), data access methods, and metadata organization and structuring tools (to support schema development and evolution). We use an example from experimental neuroscience to illustrate how these methods can work in practice.

Marsupials and Multi-Omics: Establishing New Comparative Models of Neural Crest Patterning and Craniofacial Development

Studies across vertebrates have revealed significant insights into the processes that drive craniofacial morphogenesis, yet we still know little about how distinct facial morphologies are patterned during development. Studies largely point to evolution in GRNs of cranial progenitor cell types such as neural crest cells, as the major driver underlying adaptive cranial shapes. However, this hypothesis requires further validation, particularly within suitable models amenable to manipulation. By utilizing comparative models between related species, we can begin to disentangle complex developmental systems and identify the origin of species-specific patterning. Mammals present excellent evolutionary examples to scrutinize how these differences arise, as sister clades of eutherians and marsupials possess suitable divergence times, conserved cranial anatomies, modular evolutionary patterns, and distinct developmental heterochrony in their NCC behaviours and craniofacial patterning. In this review, I lend perspectives into the current state of mammalian craniofacial biology and discuss the importance of establishing a new marsupial model, the fat-tailed dunnart, for comparative research. Through detailed comparisons with the mouse, we can begin to decipher mammalian conserved, and species-specific processes and their contribution to craniofacial patterning and shape disparity. Recent advances in single-cell multi-omics allow high-resolution investigations into the cellular and molecular basis of key developmental processes. As such, I discuss how comparative evolutionary application of these tools can provide detailed insights into complex cellular behaviours and expression dynamics underlying adaptive craniofacial evolution. Though in its infancy, the field of "comparative evo-devo-omics" presents unparalleled opportunities to precisely uncover how phenotypic differences arise during development.

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

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

Technological scouting of bi-material face masks: simulation of adherence using 3D Facial Norms

During the COVID-19 pandemic started in March 2020, the need for personal protective equipment rapidly grew as it became mandatory. The availability of a set of faces can be of great utility in designing a face mask with proper adherence and comfortability in wearing and breathing. A 3D geometry of a face with user-defined anthropometric measures was generated with Blender, a powerful development tool for creating 3D images. Using 3D Facial Norms, a free online database, it was possible to compute the mean anthropometric measures for the age groups of 17-20, 20-30, and 30-40 years old and then generate the respective faces for both genders. The adherence of an innovative face mask was then simulated with the reverse engineering software considering both the face mask and the faces rigid.

Genetic models and approaches to study orofacial clefts

INTRODUCTION

Orofacial clefts (OFCs) are common craniofacial birth defects with heterogeneous phenotype and etiology. Geneticists have applied nearly every available method and technology to further our understanding of the genetic architectures of OFCs.

OBJECTIVE

This review describes the evidence for a genetic etiology in OFCs, statistical genetic approaches employed to identify genetic causes, and how the results have shaped our current understanding of the genetic architectures of syndromic and nonsyndromic OFCs.

CONCLUSION

There has been rapid progress towards elucidating the genetic architectures of OFCs due to the availability of large collections of DNA samples from cases, controls, and families with OFCs and the consistent adoption of new methodologies and novel statistical approaches as they are developed. Genetic studies have identified rare and common variants influencing risk of OFCs in both Mendelian and complex forms of OFCs, blurring the distinctions traditional categories used in genetic studies and clinical medicine.

Single-cell analysis identifies a key role for Hhip in murine coronal suture development

Craniofacial development depends on formation and maintenance of sutures between bones of the skull. In sutures, growth occurs at osteogenic fronts along the edge of each bone, and suture mesenchyme separates adjacent bones. Here, we perform single-cell RNA-seq analysis of the embryonic, wild type murine coronal suture to define its population structure. Seven populations at E16.5 and nine at E18.5 comprise the suture mesenchyme, osteogenic cells, and associated populations. Expression of Hhip, an inhibitor of hedgehog signaling, marks a mesenchymal population distinct from those of other neurocranial sutures. Tracing of the neonatal Hhip-expressing population shows that descendant cells persist in the coronal suture and contribute to calvarial bone growth. In Hhip^-/- coronal sutures at E18.5, the osteogenic fronts are closely apposed and the suture mesenchyme is depleted with increased hedgehog signaling compared to those of the wild type. Collectively, these data demonstrate that Hhip is required for normal coronal suture development.

The Chromatin Regulator Ankrd11 Controls Palate and Cranial Bone Development

Epigenetic and chromatin regulation of craniofacial development remains poorly understood. Ankyrin Repeat Domain 11 (ANKRD11) is a chromatin regulator that has previously been shown to control neural stem cell fates via modulation of histone acetylation. ANKRD11 gene variants, or microdeletions of the 16q24.3 chromosomal region encompassing the ANKRD11 gene, cause KBG syndrome, a rare autosomal dominant congenital disorder with variable neurodevelopmental and craniofacial involvement. Craniofacial abnormalities include a distinct facial gestalt, delayed bone age, tooth abnormalities, delayed fontanelle closure, and frequently cleft or submucosal palate. Despite this, the dramatic phenotype and precise role of ANKRD11 in embryonic craniofacial development remain unexplored. Quantitative analysis of 3D images of KBG syndromic subjects shows an overall reduction in the size of the middle and lower face. Here, we report that mice with heterozygous deletion of Ankrd11 in neural crest cells (Ankrd11nchet) display a mild midfacial hypoplasia including reduced midfacial width and a persistent open fontanelle, both of which mirror KBG syndrome patient facial phenotypes. Mice with a homozygous Ankrd11 deletion in neural crest cells (Ankrd11ncko) die at birth. They show increased severity of several clinical manifestations described for KBG syndrome, such as cleft palate, retrognathia, midfacial hypoplasia, and reduced calvarial growth. At E14.5, Ankrd11 expression in the craniofacial complex is closely associated with developing bony structures, while expression at birth is markedly decreased. Conditional deletion of Ankrd11 leads to a reduction in ossification of midfacial bones, with several ossification centers failing to expand and/or fuse. Intramembranous bones show features of delayed maturation, with bone remodeling severely curtailed at birth. Palatal shelves remain hypoplastic at all developmental stages, with a local reduction in proliferation at E13.5. Our study identifies Ankrd11 as a critical regulator of intramembranous ossification and palate development and suggests that Ankrd11nchet and Ankrd11ncko mice may serve as pre-clinical models for KBG syndrome in humans.

A landmark-free morphometrics pipeline for high-resolution phenotyping: application to a mouse model of Down syndrome

Characterising phenotypes often requires quantification of anatomical shape. Quantitative shape comparison (morphometrics) traditionally uses manually located landmarks and is limited by landmark number and operator accuracy. Here, we apply a landmark-free method to characterise the craniofacial skeletal phenotype of the Dp1Tyb mouse model of Down syndrome and a population of the Diversity Outbred (DO) mouse model, comparing it with a landmark-based approach. We identified cranial dysmorphologies in Dp1Tyb mice, especially smaller size and brachycephaly (front-back shortening), homologous to the human phenotype. Shape variation in the DO mice was partly attributable to allometry (size-dependent shape variation) and sexual dimorphism. The landmark-free method performed as well as, or better than, the landmark-based method but was less labour-intensive, required less user training and, uniquely, enabled fine mapping of local differences as planar expansion or shrinkage. Its higher resolution pinpointed reductions in interior mid-snout structures and occipital bones in both the models that were not otherwise apparent. We propose that this landmark-free pipeline could make morphometrics widely accessible beyond its traditional niches in zoology and palaeontology, especially in characterising developmental mutant phenotypes.