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.     

inka1b expression in the head mesoderm is dispensable for facial cartilage development

Inka box actin regulator 1 (Inka1) is a novel protein identified in Xenopus and is found in vertebrates. While Inka1 is required for facial skeletal development in Xenopus and zebrafish, it is dispensable in mice despite its conserved expression in the cranial neural crest, indicating that Inka1 function in facial skeletal development is not conserved among vertebrates. Zebrafish bears two paralogs of inka1 (inka1a and inka1b) in the genome, with the biological roles of inka1b barely known. Here, we analyzed the expression and function of inka1b during facial skeletal development in zebrafish. inka1b was expressed sequentially in the head mesoderm adjacent to the pharyngeal pouches essential for facial skeletal development at the stage of arch segmentation. However, a loss-of-function mutation in inka1b displayed normal head development, including the pouches and facial cartilages. The normal head of inka1b mutant fish was unlikely a result of the genetic redundancy of inka1b with inka1a, given the distinct expression of inka1a and inka1b in the cranial neural crest and head mesoderm, respectively, during craniofacial development. Our findings suggest that the inka1b expression in the head mesoderm might not be essential for head development in zebrafish.     

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.     

A genome-wide association study of facial morphology identifies novel genetic loci in Han Chinese

The human face is a heritable surface with many complex sensory organs. In recent years, many genetic loci associated with facial features have been reported in different populations, yet there is a lack of studies on the Han Chinese population. Here, we report a genome-wide association study of 3 D normal human faces of 2,659 Han Chinese with autosegment phenotypes of facial morphology. We identify singlenucleotide polymorphisms(SNPs) encompassing four genomic regions showing significant associations with different facial regions, including SNPs in DENND1 B associated with the chin, SNPs among PISRT1 associated with eyes, SNPs between DCHS2 and SFRP2 associated with the nose, and SNPs in VPS13 B associated with the nose. We replicate 24 SNPs from previously reported genetic loci in different populations, whose candidate genes are DCHS2, SUPT3 H, HOXD1, SOX9, PAX3, and EDAR. These results provide a more comprehensive understanding of the genetic basis of variation in human facial morphology.     

3D facial phenotyping by biometric sibling matching used in contemporary genomic methodologies

The analysis of contemporary genomic data typically operates on one-dimensional phenotypic measurements (e.g. standing height). Here we report on a data-driven, family-informed strategy to facial phenotyping that searches for biologically relevant traits and reduces multivariate 3D facial shape variability into amendable univariate measurements, while preserving its structurally complex nature. We performed a biometric identification of siblings in a sample of 424 children, defining 1,048 sib-shared facial traits. Subsequent quantification and analyses in an independent European cohort (n = 8,246) demonstrated significant heritability for a subset of traits (0.17–0.53) and highlighted 218 genome-wide significant loci (38 also study-wide) associated with facial variation shared by siblings. These loci showed preferential enrichment for active chromatin marks in cranial neural crest cells and embryonic craniofacial tissues and several regions harbor putative craniofacial genes, thereby enhancing our knowledge on the genetic architecture of normal-range facial variation.     

Turning the tables of sex distinction in craniofacial identification: Why females possess thicker facial soft tissues than males,not vice versa

Males are universally reported to possess larger facial soft‐tissue thickness (FSTT) than females, however, this observation oversimplifies the raw data yielding an underpowered assessment of FSTT sex‐patterning where: differences are small (η² < 5%) and inconsistent (females are routinely larger than males at the cheeks). Here we investigate body‐size normalized data to assess whether more general and improved understanding of FSTT sex‐variation in humans is possible. FSTTs were measured in 52 healthy living Australians aged 18 to 30 years using B‐mode ultrasound. Participants' stature and body mass were also measured. Sex differences were calculated before and after normalization by the aforementioned body‐composition variables. Methods were repeated in three other independent samples to evaluate reproducibility: 100 American Whites and 60 American Blacks measured by B‐mode ultrasound; and 50 Turkish residents measured by regular supine CT. Compared to raw mean differences (F < M, by ?6%), females displayed much thicker FSTTs than males when normalized for body mass (F > M, by +16%). Consequently, while the sexes share similar raw values, females possess much larger FSTTs for their relatively lighter bodies. The relative FSTT difference was 2.7× larger than the raw mean difference. Sex differences in FSTT are of larger magnitude and reversed direction in mass normalized data. Contrary to popular thought, females possess much larger FSTTs than males owing to their generically lighter bodies (?18 kg). These data patterns help explain why the pooling of sex‐categorized FSTT does not jeopardize the sex‐difference—it is encoded more strongly in terms relative to body mass.     

Cranial variation in prehistoric human skeletal remains from the Marianas

Nonmetric cranial variation and facial flatness of the Pacific and circum-Pacific populations are investigated. The peoples of the Marianas, eastern Polynesia and Hawaii form a cluster and show affinities in terms of nonmetric cranial variation with the Southeast and East Asians rather than with the Jomon-Ainu, a view which is widely supported by others. Facial flatness analysis also indicates that Polynesians have different patterns of facial prominence as compared with the Jomon-Ainu. These results increase the difficulty of accepting the Jomon-Pacific cluster proposed by Brace and his coworkers. Although genetic and nonmetric cranial variation reveal relatively close relationships, the Mariana skeletons are markedly different in facial flatness and limb bone morphology from those of Polynesians. Am J Phys Anthropol 104:399–410, 1997. © 1997 Wiley-Liss, Inc.     

Comparative population genomics of latitudinal variation in Drosophila simulans and Drosophila melanogaster

Examples of clinal variation in phenotypes and genotypes across latitudinal transects have served as important models for understanding how spatially varying selection and demographic forces shape variation within species. Here, we examine the selective and demographic contributions to latitudinal variation through the largest comparative genomic study to date of Drosophila simulans and Drosophila melanogaster, with genomic sequence data from 382 individual fruit flies, collected across a spatial transect of 19 degrees latitude and at multiple time points over 2 years. Consistent with phenotypic studies, we find less clinal variation in D. simulans than D. melanogaster, particularly for the autosomes. Moreover, we find that clinally varying loci in D. simulans are less stable over multiple years than comparable clines in D. melanogaster. D. simulans shows a significantly weaker pattern of isolation by distance than D. melanogaster and we find evidence for a stronger contribution of migration to D. simulans population genetic structure. While population bottlenecks and migration can plausibly explain the differences in stability of clinal variation between the two species, we also observe a significant enrichment of shared clinal genes, suggesting that the selective forces associated with climate are acting on the same genes and phenotypes in D. simulans and D. melanogaster.     

Functional characterization of powdery mildew resistance gene MlIW172, a new Pm60 allele and its allelic variation in wild emmer wheat

Wild emmer wheat (Triticum dicoccoides, WEW) is an immediate progenitor of both the cultivated tetraploid and hexaploid wheats and it harbors rich genetic diversity against powdery mildew caused by Blumeria graminis f. sp. tritici (Bgt). A powdery mildew resistance gene MlIW172 originated from WEW accession IW172 (G-797-M) is fine mapped in a 0.048 centimorgan (cM) genetic interval on 7AL, corresponding to a genomic region spanning 233 kb, 1 Mb and 800 kb in Chinese Spring, WEW Zavitan, and T. urartu G1812, respectively. MlIW172 encodes a typical NLR protein NLR^IW172 and physically locates in an NBS-LRR gene cluster. NLR^IW172 is subsequently identified as a new allele of Pm60, and its function is validated by EMS mutagenesis and transgenic complementation. Haplotype analysis of the Pm60 alleles reveals diversifications in sequence variation in the locus and presence and absence variations (PAV) in WEW populations. Four common single nucleotide variations (SNV) are detected between the Pm60 alleles from WEW and T. urartu, indicative of speciation divergence between the two different wheat progenitors. The newly identified Pm60 alleles and haplotypes in WEW are anticipated to be valuable for breeding powdery mildew resistance wheat cultivars via marker-assisted selection.     

Inter‐ and intraspecific variation in mycotoxin tolerance: A study of four Drosophila species

Many mycophagous Drosophila species have adapted to tolerate high concentrations of mycotoxins, an ability not reported in any other eukaryotes. Although an association between mycophagy and mycotoxin tolerance has been established in many Drosophila species, the genetic mechanisms of the tolerance are unknown. This study presents the inter‐ and intraspecific variation in the mycotoxin tolerance trait. We studied the mycotoxin tolerance in four Drosophila species from four separate clades within the immigrans‐tripunctata radiation from two distinct locations. The effect of mycotoxin treatment on 20 isofemale lines per species was studied using seven gross phenotypes: survival to pupation, survival to eclosion, development time to pupation and eclosion, thorax length, fecundity, and longevity. We observed interspecific variation among four species, with D. falleni being the most tolerant, followed by D. recens, D. neotestacea, and D. tripunctata, in that order. The results also revealed geographical variation and intraspecific genetic variation in mycotoxin tolerance. This report provides the foundation for further delineating the genetic mechanisms of the mycotoxin tolerance trait.     

Drift,selection and adaptive variation in small populations of a threatened rattlesnake

An important goal of conservation genetics is to determine if the viability of small populations is reduced by a loss of adaptive variation due to genetic drift. Here, we assessed the impact of drift and selection on direct measures of adaptive variation (toxin loci encoding venom proteins) in the eastern massasauga rattlesnake (Sistrurus catenatus), a threatened reptile that exists in small isolated populations. We estimated levels of individual polymorphism in 46 toxin loci and 1,467 control loci across 12 populations of this species, and compared the results with patterns of selection on the same loci following speciation of S. catenatus and its closest relative, the western massasauga (S. tergeminus). Multiple lines of evidence suggest that both drift and selection have had observable impacts on standing adaptive variation. In support of drift effects, we found little evidence for selection on toxin variation within populations and a significant positive relationship between current levels of adaptive variation and long‐ and short‐term estimates of effective population size. However, we also observed levels of directional selection on toxin loci among populations that are broadly similar to patterns predicted from interspecific selection analyses that pre‐date the effects of recent drift, and that functional variation in these loci persists despite small short‐term effective sizes. This suggests that much of the adaptive variation present in populations may represent an example of “drift debt,” a nonequilibrium state where present‐day levels of variation overestimate the amount of functional genetic diversity present in future populations.     

High genetic and epigenetic variation of transposable elements: Potential drivers to rapid adaptive evolution for the noxious invasive weed Mikania micrantha

Why invasive species can rapidly adapt to novel environments is a puzzling question known as the genetic paradox of invasive species. This paradox is explainable in terms of transposable elements (TEs) activity, which are theorized to be powerful mutational forces to create genetic variation. Mikania micrantha, a noxious invasive weed, in this sense provides an excellent opportunity to test the explanation. The genetic and epigenetic variation of 21 invasive populations of M. micrantha in southern China have been examined by using transposon display (TD) and transposon methylation display (TMD) techniques to survey 12 TE superfamilies. Our results showed that M. micrantha populations maintained an almost equally high level of TE‐based genetic and epigenetic variation and they have been differentiated into subpopulations genetically and epigenetically. A similar positive spatial genetic and epigenetic structure pattern was observed within 300 m. Six and seven TE superfamilies presented significant genetic and epigenetic isolation by distance (IBD) pattern. In total, 59 genetic and 86 epigenetic adaptive TE loci were identified. Of them, 51 genetic and 44 epigenetic loci were found to correlate with 25 environmental variables (including precipitation, temperature, vegetation coverage, and soil metals). Twenty‐five transposon‐inserted genes were sequenced and homology‐based annotated, which are found to be involved in a variety of molecular and cellular functions. Our research consolidates the importance of TE‐associated genetic and epigenetic variation in the rapid adaptation and invasion of M. micrantha.     

Phenotypic variation of seed traits of Haloxylon ammodendron and its affecting factors

The typical desert-inhabiting shrub Haloxylon ammodendron (C. A. Mey.) Bunge is an important species that is endangered in Northwest China. In this study, variance, correlation and cluster analysis were employed to examine phenotypic variations of H. ammodendron seed traits from 14 populations in China to investigate their relationships with environmental and genetic factors. The average coefficient of variation was higher among regions (26.9%) than within regions (19.8%). Correlation analysis revealed that the diameter/thickness and fruit wing length were significantly correlated with environmental factors. A cluster dendrograms designed based on seed traits was consistent with the finding from environmental and genetic data obtained using Ward's method and the unweighted pair-group method with arithmetic means. A Mantel test revealed correlation between seed traits and environmental data (r = 0.53, P < 0.01) was greater than that between seed traits and genetic data (r = 0.25, P < 0.05). Thus, seed traits associated with morphology, mass, and germination varied significantly among the populations, are strongly affected by climatic and geographical factors, and were moderately correlated with genetic diversity. These results provide insights for designing appropriate conservation strategies for H. ammodendron.     

Cranial variation in geographically widespread dwarf gerbil Gerbillus nanus (Gerbillinae,Rodentia) populations: Isolation by distance versus adaptation to local environments

The dwarf gerbil (Gerbillus nanus) is broadly distributed in Asia, with a range that encompasses altitudinally diverse terrain, including two major mountain ranges. Previous studies have shown this species to be generally varied across its geographic range, both genetically and morphologically. Physical barriers (e.g. mountains) and geographic distance (i.e. isolation by distance [IBD]) are expected to reduce dispersal rates, and consequently could lead to cranial morphological differentiation among populations. Adaptation to local environments is also expected to lead to cranial morphological differentiation among populations. Here, I test these hypotheses by examining variation in cranial shape and size across the geographic distribution of G. nanus using geometric morphometric analysis. Based on a sample of 473 specimens from throughout its distribution, G. nanus populations do not seem to show biologically meaningful variation in cranial shape. Cranial size, on the other hand, did show geographic variation—yet, this variation does not seem to show strong patterns of IBD nor adaptation to local environments, which could indicate that the geographic variation in the cranial size of G. nanus populations may be accounted for by factors unexamined in this study.     

Experimental evidence for ecological selection on genome variation in the wild

Understanding natural selection's effect on genetic variation is a major goal in biology, but the genome‐scale consequences of contemporary selection are not well known. In a release and recapture field experiment we transplanted stick insects to native and novel host plants and directly measured allele frequency changes within a generation at 186 576 genetic loci. We observed substantial, genome‐wide allele frequency changes during the experiment, most of which could be attributed to random mortality (genetic drift). However, we also documented that selection affected multiple genetic loci distributed across the genome, particularly in transplants to the novel host. Host‐associated selection affecting the genome acted on both a known colour‐pattern trait as well as other (unmeasured) phenotypes. We also found evidence that selection associated with elevation affected genome variation, although our experiment was not designed to test this. Our results illustrate how genomic data can identify previously underappreciated ecological sources and phenotypic targets of selection.     

The contribution of additive genetic variation to personality variation: heritability of personality

Individual animals frequently exhibit repeatable differences from other members of their population, differences now commonly referred to as ‘animal personality’. Personality differences can arise, for example, from differences in permanent environmental effects―including parental and epigenetic contributors―and the effect of additive genetic variation. Although several studies have evaluated the heritability of behaviour, less is known about general patterns of heritability and additive genetic variation in animal personality. As overall variation in behaviour includes both the among-individual differences that reflect different personalities and temporary environmental effects, it is possible for personality to be largely genetically influenced even when heritability of behaviour per se is quite low. The relative contribution of additive genetic variation to personality variation can be estimated whenever both repeatability and heritability are estimated for the same data. Using published estimates to address this issue, we found that approximately 52% of animal personality variation was attributable to additive genetic variation. Thus, while the heritability of behaviour is often moderate or low, the heritability of personality is much higher. Our results therefore (i) demonstrate that genetic differences are likely to be a major contributor to variation in animal personality and (ii) support the phenotypic gambit: that evolutionary inferences drawn from repeatability estimates may often be justified.     

Spatial variation of biochemical and ecological phenotypes in Drosophila: Electrophoretic and quantitative variation

Electrophoretic variation at three enzyme loci-alcohol dehydrogenase (Adh), glycerophosphate dehydrogenase (Gpdh), triosephosphate isomerase (Tpi)- is compared in Australian Drosophila melanogaster populations at three levels of spatial heterogeneity; among breeding sites, within populations, and between populations at the geographic level. Heterogeneity at the breeding site level greatly exceeds that among adults within populations, indicating greater intermixing at the mobile adult stage than at the developmentally immature and less migratory larval stage. Heterogeneity at the microspatial level is large relative to the geographic level at two of these loci. Spatial patterns of variation in ecological phenotypes are also considered. It is argued that electrophoretic variants may contribute little to an understanding of this quantitative variation, and that a more useful approach in ecological genetics is to consider ecological phenotypes as primary data.     

Genetic diversity and structure of the noxious alien grass Praxelis clematidea in southern China

Praxelis clematidea (Asteraceae), a plant species native to South America, is a noxious weed in southern China. We examined the genetic variation and population structure of 12 populations (76 individuals) of P. clematidea from Fujian, Guangdong, and Hainan Provinces in China using inter-simple sequence repeat (ISSR) analysis. From an initial set of 69 ISSR primers, 10 were selected which yielded 80 reproducible bands. Polymorphic bands (P) were 100%, Shannon's information index (I) was 0.4226, and Nei's gene diversity (H) was 0.2791. We infer that the high levels of genetic diversity exhibited by P. clematidea may have contributed to its invasiveness. Gene flow among populations was 2.4930, which has led to homogenization. The coefficient of population differentiation (Gst = 0.1671) indicated low levels of genetic variation among populations and high levels of genetic polymorphism within populations. There was a negative correlation between population elevation and genetic diversity, while there was a significant positive correlation between genetic distance and geographic distance based on a Mantel test (r = 0.5820, P < 0.01). Some populations from different provinces clustered together in principal coordinate and UPGMA analyses indicating that human-mediated events may have contributed to the dispersal of the species.