An interactive atlas of three-dimensional syndromic facial morphology.

Craniofacial phenotyping is critical for both syndrome delineation and diagnosis because craniofacial abnormalities occur in 30% of characterized genetic syndromes. Clinical reports, textbooks, and available software tools typically provide two-dimensional, static images and illustrations of the characteristic phenotypes of genetic syndromes. In this work, we provide an interactive web application that provides three-dimensional, dynamic visualizations for the characteristic craniofacial effects of 95 syndromes. Users can visualize syndrome facial appearance estimates quantified from data and easily compare craniofacial phenotypes of different syndromes. Our application also provides a map of morphological similarity between a target syndrome and other syndromes. Finally, users can upload 3D facial scans of individuals and compare them to our syndrome atlas estimates. In summary, we provide an interactive reference for the craniofacial phenotypes of syndromes that allows for precise, individual-specific comparisons of dysmorphology.

HDAC9 structural variants disrupting TWIST1 transcriptional regulation lead to craniofacial and limb malformations.

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 the HDAC9 sequence. Based on SVs within the HDAC9-TWIST1 locus, we defined the 3'-HDAC9 sequence as a critical TWIST1 regulatory region, encompassing craniofacial TWIST1 enhancers and CTCF sites. Deletions of either Twist1 enhancers (eTw5-7Δ/Δ) or CTCF site (CTCF-5Δ/Δ) within the Hdac9 protein-coding sequence 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 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 elucidates essential components of TWIST1 transcriptional machinery that reside within the HDAC9 sequence. It suggests that SVs encompassing protein-coding sequences 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.

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.

Genome scans of facial features in East Africans and cross-population comparisons reveal novel associations.

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 open-ended 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, consistent with an early developmental origin of the facial variation. Two of these associations were in highly conserved regions showing craniofacial-specific enhancer activity during embryological development (5q31.1 and 12q21.31). Six of the 20 loci surpassed a stricter threshold accounting for multiple phenotypes with study-wide significance (p < 6.25 × 10-10). Cross-population comparisons indicated 10 association signals were shared with Europeans (seven sharing the same associated SNP), and facilitated fine-mapping of causal variants at previously reported loci. Taken together, these results may point to both shared and population-specific components to the genetic architecture of facial variation.

Privacy, bias and the clinical use of facial recognition technology: A survey of genetics professionals.

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.

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

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 FaceBase Consortium, along with other NIH-supported data resources, provide a continuously growing, dynamic and current resource for the scientific community while improving data reproducibility and fulfilling data sharing requirements.

Refining nosology by modelling variation among facial phenotypes: the RASopathies.

BACKGROUND: In clinical genetics, establishing an accurate nosology requires analysis of variations in both aetiology and the resulting phenotypes. At the phenotypic level, recognising typical facial gestalts has long supported clinical and molecular diagnosis; however, the objective analysis of facial phenotypic variation remains underdeveloped. In this work, we propose exploratory strategies for assessing facial phenotypic variation within and among clinical and molecular disease entities and deploy these techniques on cross-sectional samples of four RASopathies: Costello syndrome (CS), Noonan syndrome (NS), cardiofaciocutaneous syndrome (CFC) and neurofibromatosis type 1 (NF1). METHODS: From three-dimensional dense surface scans, we model the typical phenotypes of the four RASopathies as average 'facial signatures' and assess individual variation in terms of direction (what parts of the face are affected and in what ways) and severity of the facial effects. We also derive a metric of phenotypic agreement between the syndromes and a metric of differences in severity along similar phenotypes. RESULTS: CFC shows a relatively consistent facial phenotype in terms of both direction and severity that is similar to CS and NS, consistent with the known difficulty in discriminating CFC from NS based on the face. CS shows a consistent directional phenotype that varies in severity. Although NF1 is highly variable, on average, it shows a similar phenotype to CS. CONCLUSIONS: We established an approach that can be used in the future to quantify variations in facial phenotypes between and within clinical and molecular diagnoses to objectively define and support clinical nosologies.

Fgf8 dosage regulates jaw shape and symmetry through pharyngeal-cardiac tissue relationships.

BACKGROUND: Asymmetries in craniofacial anomalies are commonly observed. In 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 disruption to Meckel's cartilage, which is discontinuous. All skeletal elements associated with the proximal condensation are dysmorphic, exemplified by a malformed and misoriented malleus. At later stages, Fgf8 mutants exhibit syngnathia, which falls into two 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-specimen variation. CONCLUSIONS: We hypothesize that these asymmetries are linked to heart development resulting in higher levels of Fgf8 on the right side of the face, which may buffer the right side to developmental perturbations. This mouse model may facilitate future investigations of mechanisms underlying human syngnathia and facial asymmetry.

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

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

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.

The developmental-genetics of canalization.

Canalization, or robustness to genetic or environmental perturbations, is fundamental to complex organisms. While there is strong evidence for canalization as an evolved property that varies among genotypes, the developmental and genetic mechanisms that produce this phenomenon are very poorly understood. For evolutionary biology, understanding how canalization arises is important because, by modulating the phenotypic variation that arises in response to genetic differences, canalization is a determinant of evolvability. For genetics of disease in humans and for economically important traits in agriculture, this subject is important because canalization is a potentially significant cause of missing heritability that confounds genomic prediction of phenotypes. We review the major lines of thought on the developmental-genetic basis for canalization. These fall into two groups. One proposes specific evolved molecular mechanisms while the other deals with robustness or canalization as a more general feature of development. These explanations for canalization are not mutually exclusive and they overlap in several ways. General explanations for canalization are more likely to involve emergent features of development than specific molecular mechanisms. Disentangling these explanations is also complicated by differences in perspectives between genetics and developmental biology. Understanding canalization at a mechanistic level will require conceptual and methodological approaches that integrate quantitative genetics and developmental biology.

Body mass estimation from footprint size in hominins.

Although many studies relating stature to foot length have been carried out, the relationship between foot size and body mass remains poorly understood. Here we investigate this relationship in 193 adult and 50 juvenile habitually unshod/minimally shod individuals from five different populations-Machiguenga, Daasanach, Pumé, Hadzabe, and Samoans-varying greatly in body size and shape. Body mass is highly correlated with foot size, and can be predicted from foot area (maximum length × breadth) in the combined sample with an average error of about 10%. However, comparisons among populations indicate that body shape, as represented by the body mass index (BMI), has a significant effect on foot size proportions, with higher BMI samples exhibiting relatively smaller feet. Thus, we also derive equations for estimating body mass from both foot size and BMI, with BMI in footprint samples taken as an average value for a taxon or population, estimated independently from skeletal remains. Techniques are also developed for estimating body mass in juveniles, who have relatively larger feet than adults, and for converting between foot and footprint size. Sample applications are given for five Pliocene through Holocene hominin footprint samples from Laetoli (Australopithecus afarensis), Ileret (probable Homo erectus), Happisburgh (possible Homo antecessor), Le Rozel (archaic Homo sapiens), and Barcin Höyük (H. sapiens). Body mass estimates for Homo footprint samples appear reasonable when compared to skeletal estimates for related samples. However, estimates for the Laetoli footprint sample using the new formulae appear to be too high when compared to skeletal estimates for A. afarensis. Based on the proportions of A.L. 288-1, this is apparently a result of relatively large feet in this taxon. A different method using a ratio between body mass and foot area in A.L. 288-1 provides estimates more concordant with skeletal estimates and should be used for A. afarensis.

Automated syndrome diagnosis by three-dimensional facial imaging.

PURPOSE: Deep phenotyping is an emerging trend in precision medicine for genetic disease. The shape of the face is affected in 30-40% of known genetic syndromes. Here, we determine whether syndromes can be diagnosed from 3D images of human faces. METHODS: We analyzed variation in three-dimensional (3D) facial images of 7057 subjects: 3327 with 396 different syndromes, 727 of their relatives, and 3003 unrelated, unaffected subjects. We developed and tested machine learning and parametric approaches to automated syndrome diagnosis using 3D facial images. RESULTS: Unrelated, unaffected subjects were correctly classified with 96% accuracy. Considering both syndromic and unrelated, unaffected subjects together, balanced accuracy was 73% and mean sensitivity 49%. Excluding unrelated, unaffected subjects substantially improved both balanced accuracy (78.1%) and sensitivity (56.9%) of syndrome diagnosis. The best predictors of classification accuracy were phenotypic severity and facial distinctiveness of syndromes. Surprisingly, unaffected relatives of syndromic subjects were frequently classified as syndromic, often to the syndrome of their affected relative. CONCLUSION: Deep phenotyping by quantitative 3D facial imaging has considerable potential to facilitate syndrome diagnosis. Furthermore, 3D facial imaging of "unaffected" relatives may identify unrecognized cases or may reveal novel examples of semidominant inheritance.

Fully Automatic Landmarking of Syndromic 3D Facial Surface Scans Using 2D Images.

3D facial landmarks are known to be diagnostically relevant biometrics for many genetic syndromes. The objective of this study was to extend a state-of-the-art image-based 2D facial landmarking algorithm for the challenging task of 3D landmark identification on subjects with genetic syndromes, who often have moderate to severe facial dysmorphia. The automatic 3D facial landmarking algorithm presented here uses 2D image-based facial detection and landmarking models to identify 12 landmarks on 3D facial surface scans. The landmarking algorithm was evaluated using a test set of 444 facial scans with ground truth landmarks identified by two different human observers. Three hundred and sixty nine of the subjects in the test set had a genetic syndrome that is associated with facial dysmorphology. For comparison purposes, the manual landmarks were also used to initialize a non-linear surface-based registration of a non-syndromic atlas to each subject scan. Compared to the average intra- and inter-observer landmark distances of 1.1 mm and 1.5 mm respectively, the average distance between the manual landmark positions and those produced by the automatic image-based landmarking algorithm was 2.5 mm. The average error of the registration-based approach was 3.1 mm. Comparing the distributions of Procrustes distances from the mean for each landmarking approach showed that the surface registration algorithm produces a systemic bias towards the atlas shape. In summary, the image-based automatic landmarking approach performed well on this challenging test set, outperforming a semi-automatic surface registration approach, and producing landmark errors that are comparable to state-of-the-art 3D geometry-based facial landmarking algorithms evaluated on non-syndromic subjects.

Nonlinear gene expression-phenotype relationships contribute to variation and clefting in the A/WySn mouse.

BACKGROUND: Cleft lip and palate is one of the most common human birth defects, but the underlying etiology is poorly understood. The A/WySn mouse is a spontaneously occurring model of multigenic clefting in which 20% to 30% of individuals develop an orofacial cleft. Recent work has shown altered methylation at a specific retrotransposon insertion downstream of the Wnt9b locus in clefting animals, which results in decreased Wnt9b expression. RESULTS: Using a newly developed protocol that allows us to measure morphology, gene expression, and DNA methylation in the same embryo, we relate gene expression in an individual embryo directly to its three-dimensional morphology for the first time. We find that methylation at the retrotransposon relates to Wnt9b expression and morphology. IAP methylation relates to shape of the nasal process in a manner consistent with clefting. Embryos with low IAP methylation exhibit increased among-individual variance in facial shape. CONCLUSIONS: Methylation and gene expression relate nonlinearly to nasal process morphology. Individuals at one end of a continuum of phenotypic states display a clinical phenotype and increased phenotypic variation. Variable penetrance and expressivity in this model is likely determined both by among-individual variation in methylation and changes in phenotypic robustness along the underlying liability distribution for orofacial clefting.

The effect of automated landmark identification on morphometric analyses.

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

Conditional Creation and Rescue of Nipbl-Deficiency in Mice Reveals Multiple Determinants of Risk for Congenital Heart Defects.

Elucidating the causes of congenital heart defects is made difficult by the complex morphogenesis of the mammalian heart, which takes place early in development, involves contributions from multiple germ layers, and is controlled by many genes. Here, we use a conditional/invertible genetic strategy to identify the cell lineage(s) responsible for the development of heart defects in a Nipbl-deficient mouse model of Cornelia de Lange Syndrome, in which global yet subtle transcriptional dysregulation leads to development of atrial septal defects (ASDs) at high frequency. Using an approach that allows for recombinase-mediated creation or rescue of Nipbl deficiency in different lineages, we uncover complex interactions between the cardiac mesoderm, endoderm, and the rest of the embryo, whereby the risk conferred by genetic abnormality in any one lineage is modified, in a surprisingly non-additive way, by the status of others. We argue that these results are best understood in the context of a model in which the risk of heart defects is associated with the adequacy of early progenitor cell populations relative to the sizes of the structures they must eventually form.

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

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

Hybridization in human evolution: Insights from other organisms.

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

Genomewide Association Study of African Children Identifies Association of SCHIP1 and PDE8A with Facial Size and Shape.

The human face is a complex assemblage of highly variable yet clearly heritable anatomic structures that together make each of us unique, distinguishable, and recognizable. Relatively little is known about the genetic underpinnings of normal human facial variation. To address this, we carried out a large genomewide association study and two independent replication studies of Bantu African children and adolescents from Mwanza, Tanzania, a region that is both genetically and environmentally relatively homogeneous. We tested for genetic association of facial shape and size phenotypes derived from 3D imaging and automated landmarking of standard facial morphometric points. SNPs within genes SCHIP1 and PDE8A were associated with measures of facial size in both the GWAS and replication cohorts and passed a stringent genomewide significance threshold adjusted for multiple testing of 34 correlated traits. For both SCHIP1 and PDE8A, we demonstrated clear expression in the developing mouse face by both whole-mount in situ hybridization and RNA-seq, supporting their involvement in facial morphogenesis. Ten additional loci demonstrated suggestive association with various measures of facial shape. Our findings, which differ from those in previous studies of European-derived whites, augment understanding of the genetic basis of normal facial development, and provide insights relevant to both human disease and forensics.