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⁻⁸) 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⁻¹⁰). 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.     

Morphological integration of soft-tissue facial morphology in Down Syndrome and siblings

Down syndrome (DS), resulting from trisomy of chromosome 21, is the most common live-born human aneuploidy. The phenotypic expression of trisomy 21 produces variable, though characteristic, facial morphology. Although certain facial features have been documented quantitatively and qualitatively as characteristic of DS (e.g., epicanthic folds, macroglossia, and hypertelorism), all of these traits occur in other craniofacial conditions with an underlying genetic cause. We hypothesize that the typical DS face is integrated differently than the face of non-DS siblings, and that the pattern of morphological integration unique to individuals with DS will yield information about underlying developmental associations between facial regions. We statistically compared morphological integration patterns of immature DS faces (N = 53) with those of non-DS siblings (N = 54), aged 6-12 years using 31 distances estimated from 3D coordinate data representing 17 anthropometric landmarks recorded on 3D digital photographic images. Facial features are affected differentially in DS, as evidenced by statistically significant differences in integration both within and between facial regions. Our results suggest a differential affect of trisomy on facial prominences during craniofacial development.     

A genomewide linkage scan for quantitative trait loci influencing the craniofacial complex in baboons (Papio hamadryas spp.)

Numerous studies have detected significant contributions of genes to variation in development, size, and shape of craniofacial traits in a number of vertebrate taxa. This study examines 43 quantitative traits derived from lateral cephalographs of 830 baboons (Papio hamadryas) from the pedigreed population housed at the Southwest National Primate Research Center. Quantitative genetic analyses were conducted using the SOLAR analytic platform, a maximum-likelihood variance components method that incorporates all familial information for parameter estimation. Heritability estimates were significant and of moderate to high magnitude for all craniofacial traits. Additionally, 14 significant quantitative trait loci (QTL) were identified for 12 traits from the three developmental components (basicranium, splanchnocranium, and neurocranium) of the craniofacial complex. These QTL were found on baboon chromosomes (and human orthologs) PHA1 (HSA1), PHA 2 (HSA3), PHA4 (HSA6), PHA11 (HSA12), PHA13 (HSA2), PHA16 (HSA17), and PHA17 (HSA13) (PHA, P. hamadryas; HSA, Homo sapiens). This study of the genetic architecture of the craniofacial complex in baboons provides the groundwork needed to establish the baboon as an animal model for the study of genetic and nongenetic influences on craniofacial variation.     

Effect of craniofacial genotype on the relationship between morphology and feeding performance in cichlid fishes

The relationship between morphology and performance is complex, but important for understanding the adaptive nature of morphological variation. Recent studies have sought to better understand this system by illuminating the interconnectedness of different functional systems; however, the role of genetics is often overlooked. In this study, we attempt to gain insights into this relationship by examining the effect of genotypic variation at putative craniofacial loci on the relationship between morphology and feeding performance in cichlids. We studied two morphologically disparate species, as well as a morphologically intermediate hybrid population. We assessed feeding performance, jaw protrusion, and general facial morphology for each fish. We also genotyped hybrid animals at six previously identified craniofacial loci. Cichlid species were found to differ in facial geometry, kinematic morphology, and performance. Significant correlations were also noted between these variables; however, the explanatory power of facial geometry in predicting performance was relatively poor. Notably, when hybrids were grouped by genotype, the relationship between shape and performance improved. This relationship was especially robust in animals with the specialist allele at sox9b, a well‐characterized regulator of craniofacial development. These data suggest a novel role for genotype in influencing complex relationships between form and function.     

Craniofacial variation,body size and ecological factors in aboriginal populations from central Patagonia (2000–200 years B.P.)

Previous studies have shown that ecological factors had a significant role in shaping the patterns of craniofacial variation among South American populations. Here, we evaluate whether temperature and diet contributed to facial diversification in small geographic areas. Facial size and shape of 9 osteological samples from central Patagonia (Argentina) were described using 2D landmarks and semilandmarks. Data on mean annual temperature, diet composition (δ¹³C and δ¹⁵N values) and femoral head maximum breadth, used as a proxy of body mass, were obtained for each sample. We then tested the association of body mass and the ecological variables with facial morphology using spatial regression techniques and a model selection approach. Akaike Information Criterion produced disparate results for both components of facial morphology. The best model for facial size included temperature and body mass proxy, and accounted for more than 80% of variation in size. Lower temperatures were related to larger facial sizes. Body mass was negatively associated with facial size and showed no relationship with the temperature. This suggests a relatively independent variation of cranial traits and body mass at the spatial scale studied here. Facial shape was not associated with the temperature or diet composition, contrasting with the patterns observed at larger spatial scales. Our results point out that the effect of climatic variables on cranial traits might be a source of morphological differentiation not only at large scales but also in small geographic areas, and that size and shape display a differential preservation of environmental signals.     

Genetic differences among subspecies of the saddle-back tamarin (Saguinus fuscicollis):evidence from hybrids

The subspecies of saddle-back tamarins (Saguinus fuscicollis) are known to be chromatically and morphologically diverse but little is known of the genetic basis for the observed morphological variation. The morphology of first generation subspecific hybrids can be compared to that of the parental subspecies to provide information on the extent and nature of genetic differences in morphology between subspecies. We compare two groups of saddle-back tamarin hybrids (S. f. illigeri × S. f. lagonotus and S. f. illigeri × S. f. leucogenys) to pure-bred members of their parental subspecies. These crosses were examined for heterosis, caused by allele frequency differences between the subspecies in combination with directional dominance. Thirty-nine craniofacial measurements were derived from three-dimensional coordinates of landmarks on 355 adult tamarin skulls. These measurements were corrected for sex differences and differences due to environment (wild-derived vs. laboratory-born) prior to analysis of hybridity. Sex differences were minimal for these traits. Environment had a more significant effect on craniofacial morphology. Laboratory environments produce larger faces but smaller orbits, anterior cranial vaults, and cranial bases. Significant heterosis was found for many individual traits and for the first principal component representing size and size-related shape measurements in the S. f. illigeri × S. f. lagonotus cross. The smaller samples involved in the S. f. illigeri× S. f. leucogenys cross led to a much lower number of statistically significant results, although most traits did display heterosis. Heterosis for craniofacial size was nearly statistically significant. These results suggest that there are large differences in allele frequencies among these subspecies of saddle-back tamarin for genes affecting craniofacial morphology. Based on these data we suggest that these subspecies are likely to be independent, largely isolated, evolutionary units. © 1993 Wiley-Liss, Inc.     

Genetic analysis of morphological variation in Brassica oleracea using molecular markers

A cross between the open-pollinated Brassica oleracea cabbage cultivar Wisconsin Golden Acre and the hybrid broccoli cultivar Packman was used with molecular markers to investigate the genetic control of morphological variation. Twenty-two traits derived from leaf, stem, and flowering measurements were analyzed in 90 F₂ individuals that were also classified for genotype by restriction fragment length polymorphism (RFLP) markers. Seventy-two RFLP loci, which covered the mapped genome at an average of 10 map-unit intervals on all nine linkage groups, were tested individually for associations to phenotypic measurements by single factor ANOVA, and markers with significant associations (P<0.05) were used to develop multilocus models. These data were utilized to describe the location, parental contribution of alleles, magnitude of effect, and the gene action of trait loci. Single marker loci that were significantly associated (P<0.05) with trait measurements accounted for 6.7–42.7% of the phenotypic variation. Multilocus models described as much as 60.1% of the phenotypic variation for a given trait. In some cases, different related traits had common marker-locus associations with similar gene action and genotypic class ranking. The numbers, action, and linkages, of genes controlling traits estimated with marker loci in this population corresponded to estimates based on classical genetic methods from other studies using similar, or similarly-wide, crosses. There was no evidence that genome duplication accounted for a significant portion of multiple genes controlling trait loci over the entire genome, but possible duplications of trait loci were identified for two regions with linked, duplicated marker loci.     

A Genome-Wide Association Study Identifies Five Loci Influencing Facial Morphology in Europeans

Inter-individual variation in facial shape is one of the most noticeable phenotypes in humans, and it is clearly under genetic regulation; however, almost nothing is known about the genetic basis of normal human facial morphology. We therefore conducted a genome-wide association study for facial shape phenotypes in multiple discovery and replication cohorts, considering almost ten thousand individuals of European descent from several countries. Phenotyping of facial shape features was based on landmark data obtained from three-dimensional head magnetic resonance images (MRIs) and two-dimensional portrait images. We identified five independent genetic loci associated with different facial phenotypes, suggesting the involvement of five candidate genes—PRDM16, PAX3, TP63, C5orf50, and COL17A1—in the determination of the human face. Three of them have been implicated previously in vertebrate craniofacial development and disease, and the remaining two genes potentially represent novel players in the molecular networks governing facial development. Our finding at PAX3 influencing the position of the nasion replicates a recent GWAS of facial features. In addition to the reported GWA findings, we established links between common DNA variants previously associated with NSCL/P at 2p21, 8q24, 13q31, and 17q22 and normal facial-shape variations based on a candidate gene approach. Overall our study implies that DNA variants in genes essential for craniofacial development contribute with relatively small effect size to the spectrum of normal variation in human facial morphology. This observation has important consequences for future studies aiming to identify more genes involved in the human facial morphology, as well as for potential applications of DNA prediction of facial shape such as in future forensic applications.     

Use of a natural hybrid zone for genomewide association mapping of craniofacial traits in the house mouse

The identification of the genes involved in morphological variation in nature is still a major challenge. Here, we explore a new approach: we combine 178 samples from a natural hybrid zone between two subspecies of the house mouse (Mus musculus domesticus and Mus musculus musculus), and high coverage of the genome (~ 145K SNPs) to identify loci underlying craniofacial shape variation. Due to the long history of recombination in the hybrid zone, high mapping resolution is anticipated. The combination of genomes from subspecies allows the mapping of both, variation within subspecies and inter‐subspecific differences, thereby increasing the overall amount of causal genetic variation that can be detected. Skull and mandible shape were measured using 3D landmarks and geometric morphometrics. Using principal component axes as phenotypes, and a linear mixed model accounting for genetic relatedness in the mapping populations, we identified nine genomic regions associated with skull shape and 10 with mandible shape. High mapping resolution (median size of significant regions = 148 kb) enabled identification of single or few candidate genes in most cases. Some of the genes act as regulators or modifiers of signalling pathways relevant for morphological development and bone formation, including several with known craniofacial phenotypes in mice and humans. The significant associations combined explain 13% and 7% of the skull and mandible shape variation, respectively. In addition, a positive correlation was found between chromosomal length and proportion of variation explained. Our results suggest a complex genetic architecture for shape traits and support a polygenic model.     

Trisomy 21 and facial developmental instability

The most common live‐born human aneuploidy is trisomy 21, which causes Down syndrome (DS). Dosage imbalance of genes on chromosome 21 (Hsa21) affects complex gene‐regulatory interactions and alters development to produce a wide range of phenotypes, including characteristic facial dysmorphology. Little is known about how trisomy 21 alters craniofacial morphogenesis to create this characteristic appearance. Proponents of the "amplified developmental instability" hypothesis argue that trisomy 21 causes a generalized genetic imbalance that disrupts evolutionarily conserved developmental pathways by decreasing developmental homeostasis and precision throughout development. Based on this model, we test the hypothesis that DS faces exhibit increased developmental instability relative to euploid individuals. Developmental instability was assessed by a statistical analysis of fluctuating asymmetry. We compared the magnitude and patterns of fluctuating asymmetry among siblings using three‐dimensional coordinate locations of 20 anatomic landmarks collected from facial surface reconstructions in four age‐matched samples ranging from 4 to 12 years: 1) DS individuals (n = 55); 2) biological siblings of DS individuals (n = 55); 3) and 4) two samples of typically developing individuals (n = 55 for each sample), who are euploid siblings and age‐matched to the DS individuals and their euploid siblings (samples 1 and 2). Identification in the DS sample of facial prominences exhibiting increased fluctuating asymmetry during facial morphogenesis provides evidence for increased developmental instability in DS faces. We found the highest developmental instability in facial structures derived from the mandibular prominence and lowest in facial regions derived from the frontal prominence. Am J Phys Anthropol 151:49–57, 2013. © 2013 Wiley Periodicals, Inc.     

Morphological differentiation of aboriginal human populations from Tierra del Fuego (Patagonia): implications for South American peopling

This study aims to integrate the craniofacial morphological variation of southern South American populations with the results of mtDNA haplogroup variation, to discuss the South America peopling. Because the causes of morphological differentiation of Fueguian populations are still a controversial subject, the comparison with neutral variation could contribute to elucidate them. Samples of human remains from South America regions were used to analyze the evolutionary relationships. Several craniofacial traits observed in frontal and lateral view were analyzed by means of geometric morphometrics techniques, and the evolutionary relationships based on morphological and molecular data were established in base to ordination analyses. The results from the facial skeleton agree with those obtained from mtDNA haplogroup frequencies, with La Pampa/Chaco samples detached from the Patagonian samples. Hence, the same mechanism that accounts for the pattern of frequency of haplogroups could explain the variation found in facial skeleton among the samples. It is suggested that such geographic pattern of craniofacial and molecular diversity may reflect the effect of genetic drift that occurred in the small founding populations isolated by distance or geographic barriers. Conversely, the results obtained using the traits from the lateral view slightly differ from the molecular results, showing differences between southernmost Patagonian and the other samples. Therefore, mechanisms other than genetic drift (e.g., natural selection) could have acted to shape the pattern observed in some craniofacial structures present in the lateral view, characterized by the fact that the southernmost Patagonian samples display the most robust and dolichocephalic crania.     

Genetic analysis of craniofacial traits in the medaka

Family and twin studies suggest that a substantial genetic component underlies individual differences in craniofacial morphology. In the current study, we quantified 444 craniofacial traits in 100 individuals from two inbred medaka (Oryzias latipes) strains, HNI and Hd-rR. Relative distances between defined landmarks were measured in digital images of the medaka head region. A total of 379 traits differed significantly between the two strains, indicating that many craniofacial traits are controlled by genetic factors. Of these, 89 traits were analyzed via interval mapping of 184 F(2) progeny from an intercross between HNI and Hd-rR. We identified quantitative trait loci for 66 craniofacial traits. The highest logarithm of the odds score was 6.2 for linkage group (LG) 9 and 11. Trait L33, which corresponds to the ratio of head length to head height at eye level, mapped to LG9; trait V15, which corresponds to the ratio of snout length to head width measured behind the eyes, mapped to LG11. Our initial results confirm the potential of the medaka as a model system for the genetic analysis of complex traits such as craniofacial morphology.     

Genetic evidence for facial variation being a composite phenotype of cranial variation and facial soft tissue thickness

Facial and cranial variation represent a multidimensional set of highly correlated and heritable phenotypes. Little is known about the genetic basis explaining this correlation. We develop a software package ALoSFL for simultaneous localization of facial and cranial landmarks from head computed tomography (CT) images, apply it in the analysis of head CT images of 777 Han Chinese women, and obtain a set of phenotypes representing variation in face, skull and facial soft tissue thickness (FSTT). Association analysis of 301 single nucleotide polymorphisms (SNPs) from 191 distinct genomic loci previously associated with facial variation reveals an unexpected larger number of loci showing significant associations (P < 1e–3) with cranial phenotypes than expected under the null (O/E = 3.39), suggesting facial and cranial phenotypes share a substantial proportion of genetic components. Adding FSTT to a SNP-only model shows a large impact in explaining facial variance. A gene ontology analysis reveals that bone morphogenesis and osteoblast differentiation likely underlie our cranial-significant findings. Overall, this study simultaneously investigates the genetic effects on both facial and cranial variation of the same sample, supporting that facial variation is a composite phenotype of cranial variation and FSTT.     

Scaling relationships between craniofacial sexual dimorphism and body mass dimorphism in primates: implications for the fossil record

Craniofacial remains (the most abundant identifiable remains in the fossil record) potentially offer important information about body size dimorphism in extinct species. This study evaluates the scaling relationships between body mass dimorphism and different measures of craniofacial dimorphism, evaluating taxonomic differences in the magnitude and scaling of craniofacial dimorphism across higher taxonomic groups. Data on 40 dimensions from 129 primate species and subspecies demonstrate that few dimensions change proportionally with body mass dimorphism. Primates show general patterns of greater facial vs. neurocranial and orbital dimorphism, and greater dimorphism in lengths as opposed to breadths. Within any species, though, different craniofacial dimensions can yield very different reconstructions of size dimorphism. There are significant taxonomic differences in the relationships between size and craniofacial dimorphism among primate groups that can have a significant impact on reconstructions of body mass dimorphism. Hominoids tend to show lower degrees of facial dimorphism proportional to size dimorphism than other primates. This in turn implies that strong craniofacial dimorphism in Australopithecus africanus could imply very strong body size dimorphism, conflicting with the relatively modest size dimorphism inferred from postcrania. Different methods of estimating the magnitude of size dimorphism from craniofacial measurements yield similar results, and yield comparatively low percent prediction errors for a number of dimensions. However, confidence intervals for most estimates are so large as to render most estimates highly tentative.     

Variation in soft-tissue thicknesses on the human face and their relation to craniometric dimensions

The average thickness of soft tissues on parts of the face is known, but its variation has not been related to cranial morphology. To investigate this relationship, measurements of facial soft-tissue depths and craniometric dimensions were taken on adult, white Australian cadavers (17 male and 23 female). Significant correlations between many soft-tissue depths and craniometric dimensions were found, suggesting a relationship between the amount of soft tissue present on the face and the size of the underlying bony skeleton. Soft-tissue depths were highly positively correlated with each other; craniometric dimensions were correlated but to a lesser extent. Males had thicker soft tissues and larger craniometric dimensions than females; considerable overlap of ranges was also noted. Multiple regression analysis was used to produce equations predicting the soft-tissue depth at specified areas of the face from craniometric dimensions. A subsample of nine cadavers was examined for the effects of tissue embalming. Embalming caused significant initial increases in facial soft-tissue depths. Cadavers embalmed for less than 6 months had soft-tissue depths significantly greater than for fully embalmed cadavers. The evidence that facial soft-tissue thicknesses vary with craniofacial dimensions has implications for forensic identification, facial aesthetic surgery, and approximation of the facial features of extinct individuals.     

Quantitative trait locus analysis of a recombinant inbred line population derived from a Lycopersicon esculentum x Lycopersicon cheesmanii cross

Quantitative trait loci influencing fruit traits were identified by restriction fragment length polymorphism (RFLP) analysis in a population of recombinant inbred lines (RIL) derived from a cross of the cultivated tomato, Lycopersicon esculentum with a related wild species Lycopersicon cheesmanii. One hundred thirty-two polymorphic RFLP loci spaced throughout the tomato genome were scored for 97 F₈ RIL families. Fruit weight and soluble solids were measured in replicated trials during 1991 and 1992. Seed weight was measured in 1992. Significant (P<0.01 level) quantitative trait locus (QTL) associations of marker loci were identified for each trait. A total of 73 significant marker locus-trait associations were detected for the three traits measured. Fifty-three of these associations were for fruit weight and soluble solids, many of which involved marker loci signficantly associated with both traits. QTL with large effects on all three traits were detected on chromosome 6. Greater homozygosity at many loci in the RIL population as compared to F₂ populations and greater genomic coverage resulted in increased precision in the estimation of QTL effects, and large proportions of the total phenotypic variance were explained by marker class variation at significant marker loci for many traits. The RIL population was effective in detecting and discriminating among QTL for these traits previously identified in other investigations despite skewed segregation ratios at many marker loci. Large additive effects were measured at significant marker loci. Lower fruit weight, higher soluble solids, and lower seed weight were generally associated with RFLP alleles from theL. cheesmanii parent.     

Deciphering the genomic architecture of the stickleback brain with a novel multilocus gene‐mapping approach

Quantitative traits important to organismal function and fitness, such as brain size, are presumably controlled by many small‐effect loci. Deciphering the genetic architecture of such traits with traditional quantitative trait locus (QTL) mapping methods is challenging. Here, we investigated the genetic architecture of brain size (and the size of five different brain parts) in nine‐spined sticklebacks (Pungitius pungitius) with the aid of novel multilocus QTL‐mapping approaches based on a de‐biased LASSO method. Apart from having more statistical power to detect QTL and reduced rate of false positives than conventional QTL‐mapping approaches, the developed methods can handle large marker panels and provide estimates of genomic heritability. Single‐locus analyses of an F₂ interpopulation cross with 239 individuals and 15 198, fully informative single nucleotide polymorphisms (SNPs) uncovered 79 QTL associated with variation in stickleback brain size traits. Many of these loci were in strong linkage disequilibrium (LD) with each other, and consequently, a multilocus mapping of individual SNPs, accounting for LD structure in the data, recovered only four significant QTL. However, a multilocus mapping of SNPs grouped by linkage group (LG) identified 14 LGs (1–6 depending on the trait) that influence variation in brain traits. For instance, 17.6% of the variation in relative brain size was explainable by cumulative effects of SNPs distributed over six LGs, whereas 42% of the variation was accounted for by all 21 LGs. Hence, the results suggest that variation in stickleback brain traits is influenced by many small‐effect loci. Apart from suggesting moderately heritable (h² ≈ 0.15–0.42) multifactorial genetic architecture of brain traits, the results highlight the challenges in identifying the loci contributing to variation in quantitative traits. Nevertheless, the results demonstrate that the novel QTL‐mapping approach developed here has distinctive advantages over the traditional QTL‐mapping methods in analyses of dense marker panels.     

Familial resemblance for anthropometric traits in dizygotic twins and siblings

Familial resemblance for fifty anthropometric traits was studied on a sample of 45 MZ, 101 DZ twin pairs and their 125 singleton siblings. Intraclass correlation coefficients were significant for all the traits. However, resemblance within DZ twin pairs was significantly greater than within sibs for 22 variables, showing that the former had a more correlated environment than the latter. The study also showed that head and facial traits were relatively more stable to the environmental factors than the body traits and hence more suitable for cross-cultural comparisons. The study listed measures of girth and skinfold, thickness as the most labile traits.     

PLEIOTROPIC EFFECTS OF INDIVIDUAL GENE LOCI ON MANDIBULAR MORPHOLOGY

The genotypic basis of morphological variation is largely unknown. In this study we examine patterns of pleiotropic effects on mandibular morphology at individual gene loci to determine whether the pleiotropic effects of individual genes are restricted to functionally and developmentally related traits. Mandibular measurements were obtained from 480 mice from the F₂ generation of an intercross between the LG/J and SM/J mouse strains. DNA was also extracted from these animals, and 76 microsatellite loci covering the autosomes were scored. Interval mapping was used to detect chromosomal locations with significant effects on various mandibular measurements. Sets of traits mapping to a common chromosomal region were considered as being affected by a single quantitative trait locus (QTL) for mandibular morphology. Thirty-seven such chromosomal regions were identified spread throughout the autosomes. Gene effects were small to moderate with the allele derived from the LG/J strain typically leading to larger size. When dominance was present, the LG/J allele was typically dominant to the SM/J allele. Most loci affected restricted functional and developmental regions of the mandible. Of the 26 chromosomal regions affecting more than two traits, 50% affect the muscular processes of the ascending ramus, 27% affect the alveolar processes carrying the teeth, and 23% affect the whole mandible. Four additional locations affecting two traits had effects significantly associated with alveolar regions. Pleiotropic effects are typically restricted to morphologically integrated complexes.