Papers by Benedikt Hallgrimsson

Bite force performance from wild derived mice has undetectable heritability despite having heritable morphological components
bioRxiv (Cold Spring Harbor Laboratory), Dec 3, 2021
Fitness-related traits tend to have low heritabilities. Conversely, morphology tends to be highly... more Fitness-related traits tend to have low heritabilities. Conversely, morphology tends to be highly heritable. Yet, many fitness-related performance traits such as running speed or bite force depend critically on morphology. Craniofacial morphology correlates with bite performance in several groups including rodents. However, within species, this relationship is less clear, and the genetics of performance, morphology and function are rarely analyzed in combination. Here, we use a half-sib design in outbred wild-derived Mus musculus to study the morphology-bite force relationship and determine whether there is additive genetic (co-)variance for these traits. Results suggest that bite force has undetectable additive genetic variance and heritability in this sample, while morphological traits related mechanically to bite force exhibit varying levels of heritability. The most heritable traits include the length of the mandible which relates to bite force. Despite its correlation with morphology, realized bite force was not heritable, which suggests it is less responsive to selection in comparison to its morphological determinants. We explain this paradox with a non-additive, many-to-one mapping hypothesis of heritable change in complex traits. We furthermore propose that performance traits could evolve if pleiotropic relationships among the determining traits are modified.

A Deep Invertible 3-D Facial Shape Model for Interpretable Genetic Syndrome Diagnosis
IEEE Journal of Biomedical and Health Informatics, Jul 1, 2022
One of the primary difficulties in treating patients with genetic syndromes is diagnosing their c... more One of the primary difficulties in treating patients with genetic syndromes is diagnosing their condition. Many syndromes are associated with characteristic facial features that can be imaged and utilized by computer-assisted diagnosis systems. In this work, we develop a novel 3D facial surface modeling approach with the objective of maximizing diagnostic model interpretability within a flexible deep learning framework. Therefore, an invertible normalizing flow architecture is introduced to enable both inferential and generative tasks in a unified and efficient manner. The proposed model can be used (1) to infer syndrome diagnosis and other demographic variables given a 3D facial surface scan and (2) to explain model inferences to non-technical users via multiple interpretability mechanisms. The model was trained and evaluated on more than 4700 facial surface scans from subjects with 47 different syndromes. For the challenging task of predicting syndrome diagnosis given a new 3D facial surface scan, age, and sex of a subject, the model achieves a competitive overall top-1 accuracy of 71%, and a mean sensitivity of 43% across all syndrome classes. We believe that invertible models such as the one presented in this work can achieve competitive inferential performance while greatly increasing model interpretability in the domain of medical diagnosis.
A Conditional Normalizing Flow Model of Syndromic 3D Facial Shape

Development, Mar 12, 2021
Characterising phenotypes often requires quantification of anatomical shape. Quantitative shape c... more 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 landmarkfree 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 landmarkbased 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 midsnout 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.

Privacy, bias and the clinical use of facial recognition technology: A survey of genetics professionals
American Journal of Medical Genetics Part C: Seminars in Medical Genetics, Feb 7, 2023
Facial recognition technology (FRT) has been adopted as a precision medicine tool. The medical ge... more Facial recognition technology (FRT) has been adopted as a precision medicine tool. The medical genetics field highlights both the clinical potential and privacy risks of this technology, putting the discipline at the forefront of a new digital privacy debate. Investigating how geneticists perceive the privacy concerns surrounding FRT can help shape the evolution and regulation of the field, and provide lessons for medicine and research more broadly. Five hundred and sixty‐two genetics clinicians and researchers were approached to fill out a survey, 105 responded, and 80% of these completed. The survey consisted of 48 questions covering demographics, relationship to new technologies, views on privacy, views on FRT, and views on regulation. Genetics professionals generally placed a high value on privacy, although specific views differed, were context‐specific, and covaried with demographic factors. Most respondents (88%) agreed that privacy is a basic human right, but only 37% placed greater weight on it than other values such as freedom of speech. Most respondents (80%) supported FRT use in genetics, but not necessarily for broader clinical use. A sizeable percentage (39%) were unaware of FRT's lower accuracy rates in marginalized communities and of the mental health effects of privacy violations (62%), but most (76% and 75%, respectively) expressed concern when informed. Overall, women and those who self‐identified as politically progressive were more concerned about the lower accuracy rates in marginalized groups (88% vs. 64% and 83% vs. 63%, respectively). Younger geneticists were more wary than older geneticists about using FRT in genetics (28% compared to 56% “strongly” supported such use). There was an overall preference for more regulation, but respondents had low confidence in governments' or technology companies' ability to accomplish this. Privacy views are nuanced and context‐dependent. Support for privacy was high but not absolute, and clear deficits existed in awareness of crucial FRT‐related discrimination potential and mental health impacts. Education and professional guidelines may help to evolve views and practices within the field.
AlignModel: A Tool for Interactive Manual Registration of Models in 3D Space

bioRxiv (Cold Spring Harbor Laboratory), Aug 11, 2021
Structural variants (SVs) can affect protein-coding sequences as well as gene regulatory elements... more Structural variants (SVs) can affect protein-coding sequences as well as gene regulatory elements. However, SVs disrupting protein-coding sequences that also function as cis-regulatory elements remain largely uncharacterized. Here, we show that craniosynostosis patients with SVs containing the Histone deacetylase 9 (HDAC9) protein-coding sequence are associated with disruption of TWIST1 regulatory elements that reside within HDAC9 sequence. Based on SVs within the HDAC9-TWIST1 locus, we defined the 3' HDAC9 sequence (~500Kb) as a critical TWIST1 regulatory region, encompassing craniofacial TWIST1 enhancers and CTCF sites. Deletions of either Twist1 enhancers (eTw5-7 Δ/Δ) or Ctcf site (Ctcf Δ/Δ) within the Hdac9 proteincoding sequence in mice led to decreased Twist1 expression and altered anterior\posterior limb expression patterns of Shh pathway genes. This decreased Twist1 expression results in a smaller sized and asymmetric skull and polydactyly that resembles Twist1 +/− mouse phenotype. Chromatin conformation analysis revealed that the Twist1 promoter region interacts with Hdac9 sequences that encompass Twist1 enhancers and a Ctcf site and that interactions depended on the presence of both regulatory regions. Finally, a large inversion of the entire Hdac9 sequence (Hdac9 INV/+) in mice that does not disrupt Hdac9 expression but repositions Twist1 regulatory elements showed decreased Twist1 expression and led to a craniosynostosis-like phenotype and polydactyly. Thus, our study elucidated essential components of TWIST1 transcriptional machinery that reside within the HDAC9 sequence, suggesting that SVs, encompassing proteincoding sequence, such as HDAC9, could lead to a phenotype that is not attributed to its protein function but rather to a disruption of the transcriptional regulation of a nearby gene, such as TWIST1.

Development, 2020
The FaceBase Consortium was established by the National Institute of Dental and Craniofacial Rese... more The FaceBase Consortium was established by the National Institute of Dental and Craniofacial Research in 2009 as a ‘big data’ resource for the craniofacial research community. Over the past decade, researchers have deposited hundreds of annotated and curated datasets on both normal and disordered craniofacial development in FaceBase, all freely available to the research community on the FaceBase Hub website. The Hub has developed numerous visualization and analysis tools designed to promote integration of multidisciplinary data while remaining dedicated to the FAIR principles of data management (findability, accessibility, interoperability and reusability) and providing a faceted search infrastructure for locating desired data efficiently. Summaries of the datasets generated by the FaceBase projects from 2014 to 2019 are provided here. FaceBase 3 now welcomes contributions of data on craniofacial and dental development in humans, model organisms and cell lines. Collectively, the Fac...

Clinical diagnosis of syndromes benefits strongly from objective facial phenotyping. This study i... more Clinical diagnosis of syndromes benefits strongly from objective facial phenotyping. This study investigates facial dysmorphism of genetic syndromes by building and investigating a low-dimensional metric space referred to as the clinical face phenotypic space (CFPS). As a facial matching tool for clinical genetics, such CFPS can enhance clinical diagnosis. It helps to interpret facial dysmorphisms of a subject by placing them within the space of known dysmorphisms. In this paper, a triplet loss-based autoencoder developed by geometric deep learning (GDL) is trained using multi-task learning, which combines supervised and unsupervised learning approaches. Experiments are designed to illustrate the following properties of CFPSs that can aid clinicians in narrowing down their search space: A CFPS can 1) classify and cluster syndromes accurately, 2) generalize to novel syndromes, and 3) preserve the relatedness of genetic diseases, meaning that clusters of phenotypically similar disorde...

bioRxiv (Cold Spring Harbor Laboratory), Mar 18, 2022
Background: Asymmetries in craniofacial anomalies are commonly observed. With respect to the faci... more Background: Asymmetries in craniofacial anomalies are commonly observed. With respect to the facial skeleton, the left side is more commonly and/or severely affected than the right. Such asymmetries complicate treatment options. Mechanisms underlying variation in disease severity between individuals as well as within individuals (asymmetries) are still relatively unknown. Results: Developmental reductions in Fibroblast growth factor 8 (Fgf8) have a dosage dependent effect on jaw size, shape, and symmetry. Further, Fgf8 mutants have directionally asymmetric jaws with the left side being more affected than the right. Defects in lower jaw development begin with an early disruption to Meckel's cartilage, which is discontinuous and appears as two separate condensations in Fgf8 mutants. All skeletal elements associated with the proximal condensation are dysmorphic in the mutants, which is exemplified by a malformed and mis-oriented malleus. At later stages, Fgf8 mutants exhibit syngnathia, which falls into 2 broad categories: bony fusion of the maxillary and mandibular alveolar ridges and zygomatico-mandibular fusion. All of these morphological defects exhibit both inter-and intra-individual variation. Conclusions: We hypothesize that these asymmetries are linked to asymmetries in heart development resulting in higher levels of Fgf8 on the right side of the face during development, which may buffer the right side to mild developmental perturbations. This mutant mouse is a good model for future investigations of mechanisms underlying human syngnathia and facial asymmetry.
Genome Research, Jun 16, 2022
Evolutionary Biology-new York, Dec 28, 2016
the dermatocranium may contribute to a passive defence against predation. We hypothesize that the... more the dermatocranium may contribute to a passive defence against predation. We hypothesize that the complexity in dermatocranial shape demonstrated here for P. hernandesi indicates parcellation of shape variance, which may contribute to explanations of the pronounced dermatocranial disparity exhibited by the species of Phrynosoma.

PLOS Genetics, Aug 19, 2021
Facial morphology is highly variable, both within and among human populations, and a sizable port... more Facial morphology is highly variable, both within and among human populations, and a sizable portion of this variation is attributable to genetics. Previous genome scans have revealed more than 100 genetic loci associated with different aspects of normal-range facial variation. Most of these loci have been detected in Europeans, with few studies focusing on other ancestral groups. Consequently, the degree to which facial traits share a common genetic basis across diverse sets of humans remains largely unknown. We therefore investigated the genetic basis of facial morphology in an East African cohort. We applied an openended data-driven phenotyping approach to a sample of 2,595 3D facial images collected on Tanzanian children. This approach segments the face into hierarchically arranged, multivariate features that capture the shape variation after adjusting for age, sex, height, weight, facial size and population stratification. Genome scans of these multivariate shape phenotypes revealed significant (p < 2.5 × 10 −8) signals at 20 loci, which were enriched for active chromatin elements in human cranial neural crest cells and embryonic craniofacial tissue,

Evolutionary Biology-new York, Nov 9, 2022
Fitness-related traits tend to have low heritabilities. Conversely, morphology tends to be highly... more Fitness-related traits tend to have low heritabilities. Conversely, morphology tends to be highly heritable. Yet, many fitnessrelated performance traits such as running speed or bite force depend critically on morphology. Craniofacial morphology correlates with bite performance in several groups including rodents. However, within species, this relationship is less clear, and the genetics of performance, morphology and function are rarely analyzed in combination. Here, we use a halfsib design in outbred wild-derived Mus musculus to study the morphology-bite force relationship and determine whether there is additive genetic (co-)variance for these traits. Results suggest that bite force has undetectable additive genetic variance and heritability in this sample, while morphological traits related mechanically to bite force exhibit varying levels of heritability. The most heritable traits include the length of the mandible which relates to bite force. Despite its correlation with morphology, realized bite force was not heritable, which suggests it is less responsive to selection in comparison to its morphological determinants. We explain this paradox with a non-additive, many-to-one mapping hypothesis of heritable change in complex traits. We furthermore propose that performance traits could evolve if pleiotropic relationships among the determining traits are modified.
Exploring the Distribution and Orientation of Cell Proliferation as Drivers of Mouse Facial Development
The FASEB Journal, May 1, 2022
Perpetual Right to Exclusive Use of Java 3D Digital Atlas for Human Body for Medical Research and Development

Building Generic Anatomical Models Using Virtual Cutting and Iterative Registration
Using 3D generic models to statistically analyze trends in biological structure changes is an imp... more Using 3D generic models to statistically analyze trends in biological structure changes is an important tool in morphometrics research. Therefore, 3D generic models built for a range of populations are in high demand. However, due to the complexity of biological structures and the limited views of them that medical images can offer, it is still an exceptionally difficult task to quickly and accurately create 3D generic models (a model is a 3D graphical representation of a biological structure) based on medical image stacks (a stack is an ordered collection of 2D images). We show that the creation of a generic model that captures spatial information exploitable in statistical analyses is facilitated by coupling our generalized segmentation method to existing automatic image registration algorithms. The method of creating generic 3D models consists of the following processing steps: (i) scanning subjects to obtain image stacks; (ii) creating individual 3D models from the stacks; (iii) interactively extracting sub-volume by cutting each model to generate the sub-model of interest; (iv) creating image stacks that contain only the information pertaining to the sub-models; (v) iteratively registering the corresponding new 2D image stacks; (vi) averaging the newly created sub-models based on intensity to produce the generic model from all the individual sub-models. After several registration procedures are applied to the image stacks, we can create averaged image stacks with sharp boundaries. The averaged 3D model created from those image stacks is very close to the average representation of the population. The image registration time varies depending on the image size and the desired accuracy of the registration. Both volumetric data and surface model for the generic 3D model are created at the final step. Our method is very flexible and easy to use such that anyone can use image stacks to create models and retrieve a sub-region from it at their ease. Java-based implementation allows our method to be used on various visualization systems including personal computers, workstations, computers equipped with stereo displays, and even virtual reality rooms such as the CAVE Automated Virtual Environment. The technique allows biologists to build generic 3D models of their interest quickly and accurately.
A 4D Perspective on Cortical Bone Remodelling by Micro-CT

Characterising phenotypes often requires quantification of anatomical shapes. Quantitative shape ... more Characterising phenotypes often requires quantification of anatomical shapes. Quantitative shape comparison (morphometrics) traditionally uses anatomical landmarks and is therefore limited by the number of landmarks and operator accuracy when landmarks are located manually. Here we apply a landmark-free method to characterise the craniofacial skeletal phenotype of the Dp1Tyb mouse model of Down syndrome (DS), validating it against a landmark-based approach. We identify cranial dysmorphologies in Dp1Tyb mice, especially smaller size and brachycephaly (front-back shortening) homologous to the human phenotype. The landmark-free phenotyping was less labour-intensive and required less user training than the landmark-based method. It also enabled mapping of local differences as planar expansion or shrinkage. This higher resolution and local mapping pinpointed reductions in interior midsnout structures and occipital bones in this DS model that were not as apparent using a traditional landmark-based method. This approach could make morphometrics widelyaccessible beyond traditional niches in zoology and palaeontology, especially in characterising mutant phenotypes. .
Comparing 2D and 3D representations for face-based genetic syndrome diagnosis
European Journal of Human Genetics
Uploads
Papers by Benedikt Hallgrimsson