Novel genetic loci affecting facial shape variation in humans
Abstract
The human face represents a combined set of highly heritable phenotypes, but knowledge on its genetic architecture remains limited, despite the relevance for various fields. A series of genome-wide association studies on 78 facial shape phenotypes quantified from 3-dimensional facial images of 10,115 Europeans identified 24 genetic loci reaching study-wide suggestive association (p<5x10-8), among which 17 were previously unreported. A follow-up multi-ethnic study in additional 7,917 individuals confirmed 10 loci including 6 unreported ones (padjusted<2.1x10-3). A global map of derived polygenic face scores assembled facial features in major continental groups consistent with anthropological knowledge. Analyses of epigenomic datasets from cranial neural crest cells revealed abundant cis-regulatory activities at the face-associated genetic loci. Luciferase reporter assays in neural crest progenitor cells highlighted enhancer activities of several face-associated DNA variants. These results substantially advance our understanding of the genetic basis underlying human facial variation and provide candidates for future in-vivo functional studies.
Data availability
GWAS meta-analysis summary statistics data of the significantly associated SNPs are provided with the paper in the supplementary file 1. In addition, GWAS meta-analysis summary statistics of all SNPs and all facial phenotypes, including for each SNP the effect allele, non-effect allele and for each phenotype the effect size alligned to the effect allele with standard error and p-value, are made publically available via figshare under https://doi.org/10.6084/m9.figshare.10298396 (updated file). Moreover, after the paper is accepted for publication, we will upload to the EBI GWAS Catalogue the complete summary statistics of all SNPs (same information as on figshare now) into the GWAS Catalogue. We included this information and the website links in the Material and Method section.
Article and author information
Author details
Funding
European Union Horizon 2020 Research and Innovation Programme (740580 (VISAGE))
- Manfred Kayser
Levelhulm Trust (F/07 134/DF)
- Andrés Ruiz-Linares
National Natural Science Foundation of China (91631307)
- Sijia Wang
National Natural Science Foundation of China (91731303)
- Shu-Hua Xu
National Natural Science Foundation of China (30890034)
- Li Jin
Australian NHMRC
- Nicholas G Martin
Australian NHMRC Fellowship (APP1103623)
- Sarah E Medland
National Natural Science Foundation of China (31771388)
- Shu-Hua Xu
National Natural Science Foundation of China (315014)
- Shu-Hua Xu
National Natural Science Foundation of China (31711530331)
- Shu-Hua Xu
National Natural Science Foundation of China (31271338)
- Li Jin
National Science Foundation of China (91651507)
- Fan Liu
National Institute of Dental and Craniofacial Research (R01-DE027023)
- Seth M Weinberg
National Institute of Dental and Craniofacial Research (R01-DE016148)
- Seth M Weinberg
National Institute of Dental and Craniofacial Research (X01-HG007821)
- Seth M Weinberg
Netherlands Organization of Scientific Research (911-03-012)
- M Arfan Ikram
National Key R&D Program of China (2017YFC083501)
- Fan Liu
Strategic Priority Reserach Program Chinese Academy of Sciences (XDC010400100)
- Fan Liu
China Scholarship Council (PhD Fellowship)
- Ziyi Xiong
Netherlands Organization of Scientific Research (1750102005011)
- M Arfan Ikram
Wellcome Trust
- Timothy D Spector
Medical Research Council (102215/2/13/2)
- Evie Stergiakouli
National Institute of Dental and Craniofacial Research (U01-DE20078)
- Seth M Weinberg
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Reviewing Editor
- Andrew P Morris, University of Liverpool, United Kingdom
Ethics
Human subjects: All cohort participants gave informed consent and consent to publish. The different cohort studies involved have been approved by their local ethics committees and in part higher institutions such as ministries, as described in the Material and Method section. Protocol numbers can be found for each cohort in the Materials and Methods section.
Version history
- Received: July 3, 2019
- Accepted: November 22, 2019
- Accepted Manuscript published: November 25, 2019 (version 1)
- Accepted Manuscript updated: November 26, 2019 (version 2)
- Version of Record published: December 11, 2019 (version 3)
Copyright
© 2019, Xiong et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 5,600
- views
-
- 764
- downloads
-
- 59
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Genetics and Genomics
Runs-of-homozygosity (ROH) segments, contiguous homozygous regions in a genome were traditionally linked to families and inbred populations. However, a growing literature suggests that ROHs are ubiquitous in outbred populations. Still, most existing genetic studies of ROH in populations are limited to aggregated ROH content across the genome, which does not offer the resolution for mapping causal loci. This limitation is mainly due to a lack of methods for the efficient identification of shared ROH diplotypes. Here, we present a new method, ROH-DICE (runs-of-homozygous diplotype cluster enumerator), to find large ROH diplotype clusters, sufficiently long ROHs shared by a sufficient number of individuals, in large cohorts. ROH-DICE identified over 1 million ROH diplotypes that span over 100 single nucleotide polymorphisms (SNPs) and are shared by more than 100 UK Biobank participants. Moreover, we found significant associations of clustered ROH diplotypes across the genome with various self-reported diseases, with the strongest associations found between the extended human leukocyte antigen (HLA) region and autoimmune disorders. We found an association between a diplotype covering the homeostatic iron regulator (HFE) gene and hemochromatosis, even though the well-known causal SNP was not directly genotyped or imputed. Using a genome-wide scan, we identified a putative association between carriers of an ROH diplotype in chromosome 4 and an increase in mortality among COVID-19 patients (p-value = 1.82 × 10−11). In summary, our ROH-DICE method, by calling out large ROH diplotypes in a large outbred population, enables further population genetics into the demographic history of large populations. More importantly, our method enables a new genome-wide mapping approach for finding disease-causing loci with multi-marker recessive effects at a population scale.
-
- Chromosomes and Gene Expression
- Genetics and Genomics
Members of the diverse heterochromatin protein 1 (HP1) family play crucial roles in heterochromatin formation and maintenance. Despite the similar affinities of their chromodomains for di- and tri-methylated histone H3 lysine 9 (H3K9me2/3), different HP1 proteins exhibit distinct chromatin-binding patterns, likely due to interactions with various specificity factors. Previously, we showed that the chromatin-binding pattern of the HP1 protein Rhino, a crucial factor of the Drosophila PIWI-interacting RNA (piRNA) pathway, is largely defined by a DNA sequence-specific C2H2 zinc finger protein named Kipferl (Baumgartner et al., 2022). Here, we elucidate the molecular basis of the interaction between Rhino and its guidance factor Kipferl. Through phylogenetic analyses, structure prediction, and in vivo genetics, we identify a single amino acid change within Rhino’s chromodomain, G31D, that does not affect H3K9me2/3 binding but disrupts the interaction between Rhino and Kipferl. Flies carrying the rhinoG31D mutation phenocopy kipferl mutant flies, with Rhino redistributing from piRNA clusters to satellite repeats, causing pronounced changes in the ovarian piRNA profile of rhinoG31D flies. Thus, Rhino’s chromodomain functions as a dual-specificity module, facilitating interactions with both a histone mark and a DNA-binding protein.