CoChat
Sign Up
PreprintCited by 144Epigenetics and DNA MethylationPleistocene-Era Hominins and ArchaeologyForensic and Genetic Research

The DOI matches its Crossref record; 86 of the 86 references it cites were found in a registry or an index.

Genes Affecting Vocal and Facial Anatomy Went Through Extensive Regulatory Divergence in Modern Humans

Open access
Gokhman, David, Nissim-Rafinia, Malka, Agranat-Tamir, Lily, Housman, Genevieve, García-Pérez, Raquel, Lizano, Esther, Cheronet, Olivia, Mallick, Swapan, Nieves-Colón, Maria A., Li, Heng, Alpaslan-Roodenberg, Songül, Novak, Mario, Gu, Hongcang, Osinski, Jason M., Ferrando-Bernal, Manuel, Gelabert, Pere, Lipende, Iddi, Mjungu, Deus, Kondova, Ivanela, Bontrop, Ronald, Kullmer, Ottmar, Weber, Gerhard, Shahar, Tal, Dvir-Ginzberg, Mona, Faerman, Marina, Quillen, Ellen E., Meissner, Alexander, Lahav, Yonatan, Kandel, Leonid, Liebergall, Meir, Prada, María E., Vidal, Julio M., Gronostajski, Richard M., Stone, Anne C., Yakir, Benjamin, Lalueza-Fox, Carles, Pinhasi, Ron, Reich, David, Marques-Bonet, Tomas, Meshorer, Eran, Carmel, Liran·2017·Nature Communications
DOI 10.1101/106955PMC PMC7055320
Open this paper in CoChat — free accountOpen PDF

Abstract

Summary: Regulatory changes are broadly accepted as key drivers of phenotypic divergence. However, identifying regulatory changes that underlie human-specific traits has proven very challenging. Here, we use 63 DNA methylation maps of ancient and present-day humans, as well as of six chimpanzees, to detect differentially methylated regions that emerged in modern humans after the split from Neanderthals and Denisovans. We show that genes affecting the face and vocal tract went through particularly extensive methylation changes. Specifically, we identify widespread hypermethylation in a network of face- and voice-affecting genes ( SOX9 , ACAN , COL2A1 , NFIX and XYLT1 ). We propose that these repression patterns appeared after the split from Neanderthals and Denisovans, and that they might have played a key role in shaping the modern human face and vocal tract.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1101/106955

Cites

86 of 86 references verified0 citation errors · 0 not checkable

Statistics

Figures

Fig. 1 Variability filtering and lineage assignment. a Methylation levels across AMH, Denisovan, Neanderthal, and chimpanzee samples in DMR#278 (chr4:38,014,896–38,016,197). This is an example of a lineage-specific DMR, defined as a locus in which all samples of a group are found outside the range of methylation in the other groups. Chimpanzee samples were used during the following step of lineage assignment. b A putative limb-specific DMR (chr3:14,339,371–14,339,823) which was removed from the analysis, as it does not comply with our definition of lineage-specific DMRs. Femur, toe, and finger samples are hypermethylated compared to other skeletal elements. Toe and finger are found at the bottom range of limb samples, suggesting some variation in this locus within limb samples too. c A putative sex-specific DMR (chr3:72,394,336–72,396,901) which was removed from the analysis. Males are hypermethylated compared to females. d Lineage assignment using chimpanzee samples. Only DMRs that passed the previous variability filtering steps were analyzed. Each bar at the tree leaves represents a locus in a sample. Methylation levels of the locus in each sample are marked with red (methylated) and green (unmethylated). The lineage where the methylation change has likely occurred (by parsimony) is marked by a star. Branch lengths are not scaled.
Fig. 1 Variability filtering and lineage assignment. a Methylation levels across AMH, Denisovan, Neanderthal, and chimpanzee samples in DMR#278 (chr4:38,014,896–38,016,197). This is an example of a lineage-specific DMR, defined as a locus in which all samples of a group are found outside the range of methylation in the other groups. Chimpanzee samples were used during the following step of lineage assignment. b A putative limb-specific DMR (chr3:14,339,371–14,339,823) which was removed from the analysis, as it does not comply with our definition of lineage-specific DMRs. Femur, toe, and finger samples are hypermethylated compared to other skeletal elements. Toe and finger are found at the bottom range of limb samples, suggesting some variation in this locus within limb samples too. c A putative sex-specific DMR (chr3:72,394,336–72,396,901) which was removed from the analysis. Males are hypermethylated compared to females. d Lineage assignment using chimpanzee samples. Only DMRs that passed the previous variability filtering steps were analyzed. Each bar at the tree leaves represents a locus in a sample. Methylation levels of the locus in each sample are marked with red (methylated) and green (unmethylated). The lineage where the methylation change has likely occurred (by parsimony) is marked by a star. Branch lengths are not scaled.CC BY · PubMed Central
Fig. 2 Genes affecting voice and face are the most over-represented within AMH-derived DMRs. a The number of DMRs that emerged along each of the human branches. Divergence times are in thousands of years ago (kya). b Distribution of median distances (turquoise) of DMRs to randomized single nucleotide changes that separate AMHs from archaic humans and chimpanzees. Genomic positions of single nucleotide changes were allocated at random. This was repeated 10,000 times. Red arrow marks the observed distance of DMRs, showing that they tend to be significantly closer to AMH-derived single nucleotide changes than expected by chance. This suggests that some of these sequence changes might be associated with the changes in methylation. c A heat map representing the level of enrichment of each anatomical part within the AMH-derived DMRs. Only body parts that are significantly enriched (FDR < 0.05) are colored. Three skeletal parts are significantly over-represented: the face, pelvis, and larynx (voice box, marked with arrows). d Enrichment levels of anatomical parts within the most significant (top quartile, Q statistic) AMH-derived DMRs, showing a more pronounced enrichment of genes affecting vocal and facial anatomy.
Fig. 2 Genes affecting voice and face are the most over-represented within AMH-derived DMRs. a The number of DMRs that emerged along each of the human branches. Divergence times are in thousands of years ago (kya). b Distribution of median distances (turquoise) of DMRs to randomized single nucleotide changes that separate AMHs from archaic humans and chimpanzees. Genomic positions of single nucleotide changes were allocated at random. This was repeated 10,000 times. Red arrow marks the observed distance of DMRs, showing that they tend to be significantly closer to AMH-derived single nucleotide changes than expected by chance. This suggests that some of these sequence changes might be associated with the changes in methylation. c A heat map representing the level of enrichment of each anatomical part within the AMH-derived DMRs. Only body parts that are significantly enriched (FDR < 0.05) are colored. Three skeletal parts are significantly over-represented: the face, pelvis, and larynx (voice box, marked with arrows). d Enrichment levels of anatomical parts within the most significant (top quartile, Q statistic) AMH-derived DMRs, showing a more pronounced enrichment of genes affecting vocal and facial anatomy.CC BY · PubMed Central
Fig. 3 The extent of differential methylation is highest among genes associated with the voice. a Within each lineage, the fraction of differentially methylated CpGs was computed as the number of derived CpGs per 100 kb centered around the middle of each DMR. DMRs were ranked according to the fraction of derived CpG positions in their vicinity. DMRs in genes associated with the voice are marked with red lines. In AMHs, DMRs in voice-affecting genes tend to be ranked significantly higher. Although known voice-associated genes comprise less than 2% of the genome, three of the top five AMH-derived DMRs, and all top five skeleton-related AMH-derived DMRs are in genes known to associate with the voice. In archaic-derived DMRs and in simulated DMRs, voice-associated genes do not show higher ranking compared to the rest of the DMGs. t-test P-values are shown for each group. b The fraction of differentially methylated CpGs along the five chromosomes containing ACAN, SOX9, COL2A1, XYLT1, and NFIX. In each of these chromosomes, the most extensive changes are found within these genes. All five genes control facial projection and the development of the larynx.
Fig. 3 The extent of differential methylation is highest among genes associated with the voice. a Within each lineage, the fraction of differentially methylated CpGs was computed as the number of derived CpGs per 100 kb centered around the middle of each DMR. DMRs were ranked according to the fraction of derived CpG positions in their vicinity. DMRs in genes associated with the voice are marked with red lines. In AMHs, DMRs in voice-affecting genes tend to be ranked significantly higher. Although known voice-associated genes comprise less than 2% of the genome, three of the top five AMH-derived DMRs, and all top five skeleton-related AMH-derived DMRs are in genes known to associate with the voice. In archaic-derived DMRs and in simulated DMRs, voice-associated genes do not show higher ranking compared to the rest of the DMGs. t-test P-values are shown for each group. b The fraction of differentially methylated CpGs along the five chromosomes containing ACAN, SOX9, COL2A1, XYLT1, and NFIX. In each of these chromosomes, the most extensive changes are found within these genes. All five genes control facial projection and the development of the larynx.CC BY · PubMed Central
Fig. 4 Hypermethylation of SOX9, ACAN, and COL2A1 in AMHs. a Methylation levels in the AMH-derived DMRs in SOX9, ACAN, and COL2A1. AMH samples are marked with red lines, archaic human samples are marked with blue lines and chimpanzee samples are marked with gray lines. The distribution of methylation across 52 AMH samples (450K methylation arrays) is presented in red. b SOX9 and its upstream regulatory elements. AMH-derived DMRs are marked with red rectangles. Previously identified putative enhancers are marked with circles: human craniofacial (orange), human craniofacial, chimpanzee-biased (blue), human skeletal (yellow), human non-skeletal (dark gray), and mouse (light gray). Numbers above skeletal enhancers show the difference in mean bone methylation between AMHs and archaic humans (top) and between AMHs and chimpanzee (bottom). Across almost all SOX9 enhancers, AMHs are hypermethylated compared to archaic humans and the chimpanzee.
Fig. 4 Hypermethylation of SOX9, ACAN, and COL2A1 in AMHs. a Methylation levels in the AMH-derived DMRs in SOX9, ACAN, and COL2A1. AMH samples are marked with red lines, archaic human samples are marked with blue lines and chimpanzee samples are marked with gray lines. The distribution of methylation across 52 AMH samples (450K methylation arrays) is presented in red. b SOX9 and its upstream regulatory elements. AMH-derived DMRs are marked with red rectangles. Previously identified putative enhancers are marked with circles: human craniofacial (orange), human craniofacial, chimpanzee-biased (blue), human skeletal (yellow), human non-skeletal (dark gray), and mouse (light gray). Numbers above skeletal enhancers show the difference in mean bone methylation between AMHs and archaic humans (top) and between AMHs and chimpanzee (bottom). Across almost all SOX9 enhancers, AMHs are hypermethylated compared to archaic humans and the chimpanzee.CC BY · PubMed Central
Fig. 5 DNA methylation and expression pattern of NFIX in modern and archaic humans. a Methylation levels along NFIX, color-coded from green (unmethylated) to red (methylated). Methylation levels around the two AMH-derived DMRs (#24 and #167) are shown in the zoomed-in panels. These two DMRs represent the regions where the most significant methylation changes are observed, but hypermethylation of NFIX in AMHs can be seen throughout the entire gene body. Chimpanzee and present-day samples were smoothed using the same sliding window as in ancient samples to allow easier comparison. The inferred schematic regulatory evolution of NFIX is shown using a phylogenetic tree to the left of the top panel. Star marks the shift in methylation from unmethylated (green) to methylated (red). b Methylation levels in DMRs #167 and #24 vs. expression levels of NFIX across 22 AMH tissues (gray). In both DMRs, higher methylation is significantly associated with lower expression of NFIX. Ust’-Ishim, Bone1 and Bone2 methylation levels (red) are plotted against mean NFIX expression across 13 osteoblast lines. Neanderthal and Denisovan methylation levels (green) are plotted against their predicted expression levels, based on the extrapolated regression line (dashed). Standard errors are marked with dotted lines. The expression levels of NFIX in Neanderthal and Denisovan shown in these graphs are extrapolated (green dots). c Box plots of expression levels of NFIX, COL2A1, ACAN, SOX9, and XYLT1 in 89 AMH samples (red) and four mouse samples (green) from appendicular bones (limbs and pelvis). Central line shows mean, box borders show 25th and 75th percentiles, whiskers extend to the most extreme data points not considered outliers, crosses show outliers. Expression levels were converted to percentiles based on the level of gene expression compared to the rest of the genome in each sample. t-test P-values are shown for each gene.
Fig. 5 DNA methylation and expression pattern of NFIX in modern and archaic humans. a Methylation levels along NFIX, color-coded from green (unmethylated) to red (methylated). Methylation levels around the two AMH-derived DMRs (#24 and #167) are shown in the zoomed-in panels. These two DMRs represent the regions where the most significant methylation changes are observed, but hypermethylation of NFIX in AMHs can be seen throughout the entire gene body. Chimpanzee and present-day samples were smoothed using the same sliding window as in ancient samples to allow easier comparison. The inferred schematic regulatory evolution of NFIX is shown using a phylogenetic tree to the left of the top panel. Star marks the shift in methylation from unmethylated (green) to methylated (red). b Methylation levels in DMRs #167 and #24 vs. expression levels of NFIX across 22 AMH tissues (gray). In both DMRs, higher methylation is significantly associated with lower expression of NFIX. Ust’-Ishim, Bone1 and Bone2 methylation levels (red) are plotted against mean NFIX expression across 13 osteoblast lines. Neanderthal and Denisovan methylation levels (green) are plotted against their predicted expression levels, based on the extrapolated regression line (dashed). Standard errors are marked with dotted lines. The expression levels of NFIX in Neanderthal and Denisovan shown in these graphs are extrapolated (green dots). c Box plots of expression levels of NFIX, COL2A1, ACAN, SOX9, and XYLT1 in 89 AMH samples (red) and four mouse samples (green) from appendicular bones (limbs and pelvis). Central line shows mean, box borders show 25th and 75th percentiles, whiskers extend to the most extreme data points not considered outliers, crosses show outliers. Expression levels were converted to percentiles based on the level of gene expression compared to the rest of the genome in each sample. t-test P-values are shown for each gene.CC BY · PubMed Central
Fig. 6 NFIX down-regulation is associated with modern human-derived traits. a Vocal anatomy of chimpanzee and AMH. The vocal tract is the cavity from the lips to the larynx (marked by dashed lines). In AMHs, the flattening of the face together with the descent of the larynx led to approximately 1:1 proportions of the horizontal and vertical portions of the vocal tract. b Craniofacial features of the Neanderthal, healthy AMH, and AMH with Marshall-Smith or Malan syndromes. Each box shows a phenotype that occurs in the Marshall-Smith/Malan syndromes (i.e., when NFIX is partially or completely inactive). The righthand side of each box shows the observed phenotypes of individuals with the syndromes (left), healthy AMHs (middle) and Neanderthals (right). In most phenotypes, the observed phenotypes match the expected phenotypes based on NFIX expression. Binomial test P-value is shown.
Fig. 6 NFIX down-regulation is associated with modern human-derived traits. a Vocal anatomy of chimpanzee and AMH. The vocal tract is the cavity from the lips to the larynx (marked by dashed lines). In AMHs, the flattening of the face together with the descent of the larynx led to approximately 1:1 proportions of the horizontal and vertical portions of the vocal tract. b Craniofacial features of the Neanderthal, healthy AMH, and AMH with Marshall-Smith or Malan syndromes. Each box shows a phenotype that occurs in the Marshall-Smith/Malan syndromes (i.e., when NFIX is partially or completely inactive). The righthand side of each box shows the observed phenotypes of individuals with the syndromes (left), healthy AMHs (middle) and Neanderthals (right). In most phenotypes, the observed phenotypes match the expected phenotypes based on NFIX expression. Binomial test P-value is shown.CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Show all 86 references

Research works printed but not linked: 1

  1. Aiello, L. & Dean, C . An Introduction to Human Evolutionary Anatomy. (Elsevier, 2002).

Cited by

See all 144

Related

See all

Metadata

Authors
Gokhman, D., Nissim-Rafinia, M., Agranat-Tamir, L., Housman, G., García-Pérez, R., Lizano, E., Cheronet, O., Mallick, S., Nieves-Colón, M. A., Li, H., Alpaslan-Roodenberg, S., Novak, M., Gu, H., Osinski, J. M., Ferrando-Bernal, M., Gelabert, P., Lipende, I., Mjungu, D., Kondova, I., Bontrop, R., Kullmer, O., Weber, G., Shahar, T., Dvir-Ginzberg, M., Faerman, M., Quillen, E. E., Meissner, A., Lahav, Y., Kandel, L., Liebergall, M., Prada, M. E., Vidal, J. M., Gronostajski, R. M., Stone, A. C., Yakir, B., Lalueza-Fox, C., Pinhasi, R., Reich, D., Marques-Bonet, T., Meshorer, E., & Carmel, L.
Affiliations
  • Hebrew University of Jerusalem (IL)
  • Arizona State University (US)
  • Consejo Superior de Investigaciones Científicas (ES)
  • University of Vienna (AT)
  • Broad Institute (US)
  • Harvard University (US)
  • University College Dublin (IE)
  • University at Buffalo, State University of New York (US)
  • The Jane Goodall Institute (TZ)
  • Biomedical Primate Research Centre (NL)
  • Senckenberg Research Institute and Natural History Museum Frankfurt/M (DE)
  • Shaare Zedek Medical Center (IL)
  • Texas Biomedical Research Institute (US)
  • Hadassah Medical Center (IL)
  • Junta de Castilla y León (ES)
Corresponding author
  • Gokhman, David
  • Meshorer, Eran
  • Carmel, Liran
ORCID
Venue
Nature Communications
Publisher
openRxiv
Published
2017-02-08
Volume
11
Issue
1
Pages
1189-1189
ISSN
2041-1723 (ISSN-L)
DOI
10.1101/106955
PMC
PMC7055320
OpenAlex
W3009791849
Cited by
144
References
86
Citations by year
2026: 37 · 2025: 23 · 2024: 13 · 2023: 14 · 2022: 21 · 2021: 18 · 2020: 10 · 2019: 6 · 2018: 1 · 2017: 1
FWCI
5.11
Type
journal article
Language
en
Open access
Licence
Locations
  • Nature Communications (publishedVersion, cc-by)
  • bioRxiv (Cold Spring Harbor Laboratory) (acceptedVersion, cc-by-nc-nd)
  • Europe PMC (PubMed Central) (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • Repositori digital de la UPF (Universitat Pompeu Fabra) (publishedVersion, cc-by)
  • LA Referencia (Red Federada de Repositorios Institucionales de Publicaciones Científicas) (submittedVersion, cc-by)
  • RePEc: Research Papers in Economics (submittedVersion)
  • Dipòsit Digital de Documents de la UAB (Universitat Autònoma de Barcelona) (submittedVersion, cc-by)
  • DIGITAL.CSIC (Spanish National Research Council (CSIC)) (submittedVersion, cc-by)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
  • Utrecht University Repository (Utrecht University) (submittedVersion, other-oa)
  • MPG.PuRe (Max Planck Society) (submittedVersion, other-oa)
Funders
  • National Science Foundation (MH106874, 1052693)
  • Howard Hughes Medical Institute
  • Paul G. Allen Family Foundation
  • Leakey Foundation
  • Nacey Maggioncalda Foundation
  • New York State Stem Cell Science (C030133)
  • Centres de Recerca de Catalunya
  • Clore Israel Foundation
  • Generalitat de Catalunya (BFU2017-86471-P)
  • Israel Science Foundation (ISF 1430/13)
  • Max-Planck-Gesellschaft
  • National Institutes of Health (MH106874, P51OD011133)
  • School of Human Evolution and Social Change, Arizona State University
  • National Institute of Mental Health (U01MH106874)
  • National Institute of Allergy and Infectious Diseases (R01AI058715)
  • Division of Materials Research
  • FP7 Ideas: European Research Council (263441)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Epigenetics and DNA Methylation
  • Social Sciences › Anthropology › Pleistocene-Era Hominins and Archaeology
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Genetics › Forensic and Genetic Research
Keywords
DNA methylation, differentially methylated regions, human evolution, Neanderthals, Denisovans, human-specific traits, vocal tract anatomy, SOX9, aggrecan, COL2A1, NFIX, XYLT1, hypermethylation, epigenetic regulation, chimpanzees
MeSH terms
DNA, Ancient, Adult, Aged, Animals, Cells, Cultured, Child, Pan troglodytes, Face, Female, Humans, Larynx, Male, Middle Aged, Phenotype, Phonation, Tongue, Vocal Cords, Vocalization, Animal, Evolution, Molecular, DNA Methylation, Chondrocytes, Genetic Speciation, Gene Regulatory Networks, Neanderthals, Primary Cell Culture

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review

How to cite

Cite CoChat for the checks on this page, and the paper for its own findings.

This page
CoChat. (2026). Citation check: Genes Affecting Vocal and Facial Anatomy Went Through Extensive Regulatory Divergence in Modern Humans. https://app.cochat.ai/papers/10.1101/106955
The paper
Gokhman, D., Nissim-Rafinia, M., Agranat-Tamir, L., Housman, G., García-Pérez, R., Lizano, E., Cheronet, O., Mallick, S., Nieves-Colón, M. A., Li, H., Alpaslan-Roodenberg, S., Novak, M., Gu, H., Osinski, J. M., Ferrando-Bernal, M., Gelabert, P., Lipende, I., Mjungu, D., Kondova, I., . . . Carmel, L. (2017). Genes Affecting Vocal and Facial Anatomy Went Through Extensive Regulatory Divergence in Modern Humans. Nature Communications. https://doi.org/10.1101/106955