skip to main content


Title: 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.  more » « less
Award ID(s):
1711847
NSF-PAR ID:
10300498
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; « less
Editor(s):
Tishkoff, Sarah A.
Date Published:
Journal Name:
PLOS Genetics
Volume:
17
Issue:
8
ISSN:
1553-7404
Page Range / eLocation ID:
e1009695
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Deciphering the genetic basis of vertebrate craniofacial variation is a longstanding biological problem with broad implications in evolution, development, and human pathology. One of the most stunning examples of craniofacial diversification is the adaptive radiation of birds, in which the beak serves essential roles in virtually every aspect of their life histories. The domestic pigeon (Columba livia) provides an exceptional opportunity to study the genetic underpinnings of craniofacial variation because of its unique balance of experimental accessibility and extraordinary phenotypic diversity within a single species. We used traditional and geometric morphometrics to quantify craniofacial variation in an F2laboratory cross derived from the straight‐beaked Pomeranian Pouter and curved‐beak Scandaroon pigeon breeds. Using a combination of genome‐wide quantitative trait locus scans and multi‐locus modeling, we identified a set of genetic loci associated with complex shape variation in the craniofacial skeleton, including beak shape, braincase shape, and mandible shape. Some of these loci control coordinated changes between different structures, while others explain variation in the size and shape of specific skull and jaw regions. We find that in domestic pigeons, a complex blend of both independent and coupled genetic effects underlie three‐dimensional craniofacial morphology.

     
    more » « less
  2. Abstract

    Understanding the molecular basis of repeated evolution improves our ability to predict evolution across the tree of life. Only since the last decade has high‐throughput sequencing enabled comparative genome scans to thoroughly examine the repeatability of genetic changes driving repeated phenotypic evolution. The Asian corn borer (ACB),Ostrinia furnacalis(Guenée), and the European corn borer (ECB),Ostrinia nubilalis(Hübner), are two closely related moths displaying repeatable phenological adaptation to a wide range of climates on two separate continents, largely manifesting as changes in the timing of diapause induction and termination across latitude. Candidate genes underlying diapause variation in North American ECB have been previously identified. Here, we sampled seven ACB populations across 23 degrees of latitude in China to elucidate the genetic basis of diapause variation and evolutionary mechanisms driving parallel clinal responses in the two species. Using pooled whole‐genome sequencing (Pool‐seq) data, population genomic analyses revealed hundreds of single nucleotide polymorphisms (SNP) whose allele frequencies covaried with mean diapause phenotypes along the cline. Genes involved in circadian rhythm were over‐represented among candidate genes with strong signatures of spatially varying selection. Only one of two circadian clock genes associated with diapause evolution in ECB showed evidence of reuse in ACB (period [per]), butperalleles were not shared between species nor with their outgroup, implicating independent mutational paths. Nonetheless, evidence of adaptive introgression was discovered at putative diapause loci located elsewhere in the genome, suggesting that de novo mutations and introgression might both underlie the repeated phenological evolution.

     
    more » « less
  3. Wittkopp, Patricia (Ed.)
    Abstract Investigating closely related species that rapidly evolved divergent feeding morphology is a powerful approach to identify genetic variation underlying variation in complex traits. This can also lead to the discovery of novel candidate genes influencing natural and clinical variation in human craniofacial phenotypes. We combined whole-genome resequencing of 258 individuals with 50 transcriptomes to identify candidate cis-acting genetic variation underlying rapidly evolving craniofacial phenotypes within an adaptive radiation of Cyprinodon pupfishes. This radiation consists of a dietary generalist species and two derived trophic niche specialists—a molluscivore and a scale-eating species. Despite extensive morphological divergence, these species only diverged 10 kya and produce fertile hybrids in the laboratory. Out of 9.3 million genome-wide SNPs and 80,012 structural variants, we found very few alleles fixed between species—only 157 SNPs and 87 deletions. Comparing gene expression across 38 purebred F1 offspring sampled at three early developmental stages, we identified 17 fixed variants within 10 kb of 12 genes that were highly differentially expressed between species. By measuring allele-specific expression in F1 hybrids from multiple crosses, we found that the majority of expression divergence between species was explained by trans-regulatory mechanisms. We also found strong evidence for two cis-regulatory alleles affecting expression divergence of two genes with putative effects on skeletal development (dync2li1 and pycr3). These results suggest that SNPs and structural variants contribute to the evolution of novel traits and highlight the utility of the San Salvador Island pupfish system as an evolutionary model for craniofacial development. 
    more » « less
  4. Abstract

    Adaptive radiations are often characterized by the rapid evolution of traits associated with divergent feeding modes. For example, the evolutionary history of African cichlids is marked by repeated and coordinated shifts in skull, trophic, fin and body shape. Here, we seek to explore the molecular basis for fin shape variation in Lake Malawi cichlids. We first described variation within an F2mapping population derived by crossing two cichlid species with divergent morphologies including fin shape. We then used this population to genetically map loci that influence variation in this trait. We found that the genotype–phenotype map for fin shape is largely distinct from other morphological characters including body and craniofacial shape. These data suggest that key aspects of fin, body and jaw shape are genetically modular and that the coordinated evolution of these traits in cichlids is more likely due to common selective pressures than to pleiotropy or linkage. We next combined genetic mapping data with population‐level genome scans to identifywnt7aaandcol1a1as candidate genes underlying variation in the number of pectoral fin ray elements. Gene expression patterns across species with different fin morphologies and small molecule manipulation of the Wnt pathway during fin development further support the hypothesis that variation at these loci underlies divergence in fin shape between cichlid species. In all, our data provide additional insights into the genetic and molecular mechanisms associated with morphological divergence in this important adaptive radiation.

     
    more » « less
  5. Abstract

    The genetic architecture of phenotypic traits can affect the mode and tempo of trait evolution. Human‐altered environments can impose strong natural selection, where successful evolutionary adaptation requires swift and large phenotypic shifts. In these scenarios, theory predicts that adaptation is due to a few adaptive variants of large effect, but empirical studies that have revealed the genetic architecture of rapidly evolved phenotypes are rare, especially for populations inhabiting polluted environments.Funduluskillifish have repeatedly evolved adaptive resistance to extreme pollution in urban estuaries. Prior studies, including genome scans for signatures of natural selection, have revealed some of the genes and pathways important for evolved pollution resistance, and provide context for the genotype–phenotype association studies reported here. We created multiple quantitative trait locus (QTL) mapping families using progenitors from four different resistant populations, and using RAD‐seq genetically mapped variation in sensitivity (developmental perturbations) following embryonic exposure to a model toxicant PCB‐126. We found that one to two large‐effect QTL loci accounted for resistance to PCB‐mediated developmental toxicity. QTLs harbored candidate genes that govern the regulation of aryl hydrocarbon receptor (AHR) signaling. One QTL locus was shared across all populations and another was shared across three populations. One QTL locus showed strong signatures of recent natural selection in the corresponding wild population but another QTL locus did not. Some candidate genes for PCB resistance inferred from genome scans in wild populations were identified as QTL, but some key candidate genes were not. We conclude that rapidly evolved resistance to the developmental defects normally caused by PCB‐126 is governed by few genes of large effect. However, other aspects of resistance beyond developmental phenotypes may be governed by additional loci, such that comprehensive resistance to PCB‐126, and to the mixtures of chemicals that distinguish urban estuaries more broadly, may be more genetically complex.

     
    more » « less