Reconstructing Denisovan Anatomy Using DNA Methylation Maps

Advances in skeletal genomics research across tissues and cells

Phenotypic variation within the skeleton has biological, behavioral, and biomedical functional implications for individuals and species. Thus, it is critical to understand how genomic, environmental, and mediating regulatory factors combine and interact to drive skeletal trait development and evolution. Recent research efforts to clarify these mechanisms have been made possible by expanded collections of genomic and phenotypic data from in vivo skeletal tissues, as well as the development of relevant in vitro skeletal cell culture systems. This review outlines this current work and recommends that continued exploration of this complexity should include an increased focus on how interactions between genomic and physiologically relevant contexts contribute to skeletal trait variation at population and evolutionary scales.

Craniometric variation and the ancestry of modern humans

OBJECTIVES

Ancient and contemporary DNA provide information about geographic variation in the ancestry of present-day humans. All living populations have ancestry from early Homo sapiens originating in sub-Saharan Africa. Populations of Eurasian descent also have a small amount of Neandertal ancestry. This study examines whether craniometric distances between recent modern human samples reflect this geographic variation in ancestry. Among recent modern humans, Eurasians are expected to be more similar to Neandertals, whereas both sub-Saharan Africans and Eurasians are expected to be equidistant from early H. sapiens.

MATERIALS AND METHODS

Data on 33 craniometric traits from 2524 recent modern humans were compared with data from the literature for Neandertals and early H. sapiens. Mahalanobis distances were computed for each modern specimen to both the Neandertal and early H. sapiens means. These distances were examined for differences between recent humans from sub-Saharan Africa (N = 373) and those of Eurasian descent (N = 2151).

RESULTS

Eurasians as a group are significantly closer than sub-Saharan Africans to Neandertals. There is no significant difference between the distances of sub-Saharan Africans and Eurasians to early H. sapiens.

DISCUSSION

The differences between sub-Saharan Africans and Eurasians for both Neandertals and early H. sapiens are as expected. Although there has been geographic differentiation among recent modern humans, including differences in Neandertal admixture, these differences have not affected overall similarity of recent modern sub-Saharan Africans and Eurasians to the earliest samples of H. sapiens.

A multi-ancestry GWAS meta-analysis of facial features and its application in predicting archaic human features

Facial morphology, a complex trait influenced by genetics, holds great significance in evolutionary research. However, due to limited fossil evidence, the facial characteristics of Neanderthals and Denisovans have remained largely unknown. In this study, we conducted a large-scale multi-ethnic meta-analysis of the genome-wide association study (GWAS), including 9674 East Asians and 10,115 Europeans, quantitatively assessing 78 facial traits using 3D facial images. We identified 71 genomic loci associated with facial features, including 21 novel loci. We developed a facial polygenic score (FPS) that enables the prediction of facial features based on genetic information. Interestingly, the distribution of FPSs among populations from diverse continental groups exhibited relevant correlations with observed facial features. Furthermore, we applied the FPS to predict the facial traits of seven Neanderthals and one Denisovan using ancient DNA and aligned predictions with the fossil records. Our results suggested that Neanderthals and Denisovans likely shared similar facial features, such as a wider but shorter nose and a wider endocanthion distance. The decreased mouth width was characterized specifically in Denisovans. The integration of genomic data and facial trait analysis provides valuable insights into the evolutionary history and adaptive changes in human facial morphology.

Denisovan admixture facilitated environmental adaptation in Papua New Guinean populations

Neandertals and Denisovans, having inhabited distinct regions in Eurasia and possibly Oceania for over 200,000 y, experienced ample time to adapt to diverse environmental challenges these regions presented. Among present-day human populations, Papua New Guineans (PNG) stand out as one of the few carrying substantial amounts of both Neandertal and Denisovan DNA, a result of past admixture events with these archaic human groups. This study investigates the distribution of introgressed Denisovan and Neandertal DNA within two distinct PNG populations, residing in the highlands of Mt Wilhelm and the lowlands of Daru Island. These locations exhibit unique environmental features, some of which may parallel the challenges that archaic humans once confronted and adapted to. Our results show that PNG highlanders carry higher levels of Denisovan DNA compared to PNG lowlanders. Among the Denisovan-like haplotypes with higher frequencies in highlander populations, those exhibiting the greatest frequency difference compared to lowlander populations also demonstrate more pronounced differences in population frequencies than frequency-matched nonarchaic variants. Two of the five most highly differentiated of those haplotypes reside in genomic areas linked to brain development genes. Conversely, Denisovan-like haplotypes more frequent in lowlanders overlap with genes associated with immune response processes. Our findings suggest that Denisovan DNA has provided genetic variation associated with brain biology and immune response to PNG genomes, some of which might have facilitated adaptive processes to environmental challenges.

Predicting the direction of phenotypic difference

Predicting phenotypes from genomic data is a key goal in genetics, but for most complex phenotypes, predictions are hampered by incomplete genotype-to-phenotype mapping. Here, we describe a more attainable approach than quantitative predictions, which is aimed at qualitatively predicting phenotypic differences. Despite incomplete genotype-to-phenotype mapping, we show that it is relatively easy to determine which of two individuals has a greater phenotypic value. This question is central in many scenarios, e.g., comparing disease risk between individuals, the yield of crop strains, or the anatomy of extinct vs extant species. To evaluate prediction accuracy, i.e., the probability that the individual with the greater predicted phenotype indeed has a greater phenotypic value, we developed an estimator of the ratio between known and unknown effects on the phenotype. We evaluated prediction accuracy using human data from tens of thousands of individuals from either the same family or the same population, as well as data from different species. We found that, in many cases, even when only a small fraction of the loci affecting a phenotype is known, the individual with the greater phenotypic value can be identified with over 90% accuracy. Our approach also circumvents some of the limitations in transferring genetic association results across populations. Overall, we introduce an approach that enables accurate predictions of key information on phenotypes - the direction of phenotypic difference - and suggest that more phenotypic information can be extracted from genomic data than previously appreciated.

A regulatory variant impacting TBX1 expression contributes to basicranial morphology in Homo sapiens

Changes in gene regulatory elements play critical roles in human phenotypic divergence. However, identifying the base-pair changes responsible for the distinctive morphology of Homo sapiens remains challenging. Here, we report a noncoding single-nucleotide polymorphism (SNP), rs41298798, as a potential causal variant contributing to the morphology of the skull base and vertebral structures found in Homo sapiens. Screening for differentially regulated genes between Homo sapiens and extinct relatives revealed 13 candidate genes associated with basicranial development, with TBX1, implicated in DiGeorge syndrome, playing a pivotal role. Epigenetic markers and in silico analyses prioritized rs41298798 within a TBX1 intron for functional validation. CRISPR editing revealed that the 41-base-pair region surrounding rs41298798 modulates gene expression at 22q11.21. The derived allele of rs41298798 acts as an allele-specific enhancer mediated by E2F1, resulting in increased TBX1 expression levels compared to the ancestral allele. Tbx1-knockout mice exhibited skull base and vertebral abnormalities similar to those seen in DiGeorge syndrome. Phenotypic differences associated with TBX1 deficiency are observed between Homo sapiens and Neanderthals (Homo neanderthalensis). In conclusion, the regulatory divergence of TBX1 contributes to the formation of skull base and vertebral structures found in Homo sapiens.

Reconstructing DNA methylation maps of ancient populations

Studying premortem DNA methylation from ancient DNA (aDNA) provides a proxy for ancient gene activity patterns, and hence valuable information on evolutionary changes in gene regulation. Due to statistical limitations, current methods to reconstruct aDNA methylation maps are constrained to high-coverage shotgun samples, which comprise a small minority of available ancient samples. Most samples are sequenced using in-situ hybridization capture sequencing which targets a predefined set of genomic positions. Here, we develop methods to reconstruct aDNA methylation maps of samples that were not sequenced using high-coverage shotgun sequencing, by way of pooling together individuals to obtain a DNA methylation map that is characteristic of a population. We show that the resulting DNA methylation maps capture meaningful biological information and allow for the detection of differential methylation across populations. We offer guidelines on how to carry out comparative studies involving ancient populations, and how to control the rate of falsely discovered differentially methylated regions. The ability to reconstruct DNA methylation maps of past populations allows for the development of a whole new frontier in paleoepigenetic research, tracing DNA methylation changes throughout human history, using data from thousands of ancient samples.

More than a decade of genetic research on the Denisovans

Denisovans, a group of now extinct humans who lived in Eastern Eurasia in the Middle and Late Pleistocene, were first identified from DNA sequences just over a decade ago. Only ten fragmentary remains from two sites have been attributed to Denisovans based entirely on molecular information. Nevertheless, there has been great interest in using genetic data to understand Denisovans and their place in human history. From the reconstruction of a single high-quality genome, it has been possible to infer their population history, including events of admixture with other human groups. Additionally, the identification of Denisovan DNA in the genomes of present-day individuals has provided insights into the timing and routes of dispersal of ancient modern humans into Asia and Oceania, as well as the contributions of archaic DNA to the physiology of present-day people. In this Review, we synthesize more than a decade of research on Denisovans, reconcile controversies and summarize insights into their population history and phenotype. We also highlight how our growing knowledge about Denisovans has provided insights into our own evolutionary history.

A Neanderthal/Denisovan GLI3 variant contributes to anatomical variations in mice

Changes in genomic structures underlie phenotypic diversification in organisms. Amino acid-changing mutations affect pleiotropic functions of proteins, although little is known about how mutated proteins are adapted in existing developmental programs. Here we investigate the biological effects of a variant of the GLI3 transcription factor (GLI3R1537C) carried in Neanderthals and Denisovans, which are extinct hominins close to modern humans. R1537C does not compromise protein stability or GLI3 activator-dependent transcriptional activities. In contrast, R1537C affects the regulation of downstream target genes associated with developmental processes. Furthermore, genome-edited mice carrying the Neanderthal/Denisovan GLI3 mutation exhibited various alterations in skeletal morphology. Our data suggest that an extinct hominin-type GLI3 contributes to species-specific anatomical variations, which were tolerated by relaxed constraint in developmental programs during human evolution.

Resurrecting the alternative splicing landscape of archaic hominins using machine learning

Alternative splicing contributes to adaptation and divergence in many species. However, it has not been possible to directly compare splicing between modern and archaic hominins. Here, we unmask the recent evolution of this previously unobservable regulatory mechanism by applying SpliceAI, a machine-learning algorithm that identifies splice-altering variants (SAVs), to high-coverage genomes from three Neanderthals and a Denisovan. We discover 5,950 putative archaic SAVs, of which 2,186 are archaic-specific and 3,607 also occur in modern humans via introgression (244) or shared ancestry (3,520). Archaic-specific SAVs are enriched in genes that contribute to traits potentially relevant to hominin phenotypic divergence, such as the epidermis, respiration and spinal rigidity. Compared to shared SAVs, archaic-specific SAVs occur in sites under weaker selection and are more common in genes with tissue-specific expression. Further underscoring the importance of negative selection on SAVs, Neanderthal lineages with low effective population sizes are enriched for SAVs compared to Denisovan and shared SAVs. Finally, we find that nearly all introgressed SAVs in humans were shared across the three Neanderthals, suggesting that older SAVs were more tolerated in human genomes. Our results reveal the splicing landscape of archaic hominins and identify potential contributions of splicing to phenotypic differences among hominins.

Three-dimensional geometric morphometric study of the Xuchang 2 cranium

Results of traditional metric and nonmetric assessments suggest that the Xuchang hominin shares features with Neanderthals. To comprehensively compare the nuchal morphology of XC 2 to those of the genus Homo, we conduct a three-dimensional geometric morphometric study with 35 cranial landmarks and surface semilandmarks of XC 2, Homo erectus, Middle Pleistocene humans, Neanderthals, and early and recent modern humans. Results reveal that the centroid size of XC 2 is larger than that of early and recent modern humans and can only be compared to that of Middle Pleistocene humans and H. erectus. Early and recent modern humans share a nuchal morphology distinct from archaic hominins (Ngandong H. erectus, Middle Pleistocene humans, and Neanderthals), except for SM 3, Sangiran 17, and Asian and African H. erectus. Although Ngandong specimens differ from the other H. erectus, it is unclear whether this represents a temporal or spatial trend in the process of evolution of this species. The nuchal morphological resemblance between Middle Pleistocene humans and Neanderthals may be attributed to similar cranial architecture and cerebellar shape. The great nuchal morphological variation shared by recent modern humans may indicate a particular developmental pattern. In conclusion, the nuchal morphology of different human groups is highly variable and may be caused by different factors including brain globularization and developmental plasticity. XC 2 shares similar nuchal morphology with Middle Pleistocene humans and Neanderthals, but these results are insufficient to fully resolve the taxonomic status of XC 2.

DNA methylation-based profiling of horse archaeological remains for age-at-death and castration

Age profiling of archaeological bone assemblages can inform on past animal management practices, but is limited by the fragmentary nature of the fossil record and the lack of universal skeletal markers for age. DNA methylation clocks offer new, albeit challenging, alternatives for estimating the age-at-death of ancient individuals. Here, we take advantage of the availability of a DNA methylation clock based on 31,836 CpG sites and dental age markers in horses to assess age predictions in 84 ancient remains. We evaluate our approach using whole-genome sequencing data and develop a capture assay providing reliable estimates for only a fraction of the cost. We also leverage DNA methylation patterns to assess castration practice in the past. Our work opens for a deeper characterization of past husbandry and ritual practices and holds the potential to reveal age mortality profiles in ancient societies, once extended to human remains.

Beringia and the peopling of the Western Hemisphere

Did Beringian environments represent an ecological barrier to humans until less than 15 000 years ago or was access to the Americas controlled by the spatial-temporal distribution of North American ice sheets? Beringian environments varied with respect to climate and biota, especially in the two major areas of exposed continental shelf. The East Siberian Arctic Shelf ('Great Arctic Plain' (GAP)) supported a dry steppe-tundra biome inhabited by a diverse large-mammal community, while the southern Bering-Chukchi Platform ('Bering Land Bridge' (BLB)) supported mesic tundra and probably a lower large-mammal biomass. A human population with west Eurasian roots occupied the GAP before the Last Glacial Maximum (LGM) and may have accessed mid-latitude North America via an interior ice-free corridor. Re-opening of the corridor less than 14 000 years ago indicates that the primary ancestors of living First Peoples, who already had spread widely in the Americas at this time, probably dispersed from the NW Pacific coast. A genetic 'arctic signal' in non-arctic First Peoples suggests that their parent population inhabited the GAP during the LGM, before their split from the former. We infer a shift from GAP terrestrial to a subarctic maritime economy on the southern BLB coast before dispersal in the Americas from the NW Pacific coast.

Evolution of Homo in the Middle and Late Pleistocene

The Middle and Late Pleistocene is arguably the most interesting period in human evolution. This broad period witnessed the evolution of our own lineage, as well as that of our sister taxon, the Neanderthals, and related Denisovans. It is exceptionally rich in both fossil and archaeological remains, and uniquely benefits from insights gained through molecular approaches, such as paleogenetics and paleoproteomics, that are currently not widely applicable in earlier contexts. This wealth of information paints a highly complex picture, often described as 'the Muddle in the Middle,' defying the common adage that 'more evidence is needed' to resolve it. Here we review competing phylogenetic scenarios and the historical and theoretical developments that shaped our approaches to the fossil record, as well as some of the many remaining open questions associated with this period. We propose that advancing our understanding of this critical time requires more than the addition of data and will necessitate a major shift in our conceptual and theoretical framework.

The contribution of Neanderthal introgression to modern human traits

Neanderthals, our closest extinct relatives, lived in western Eurasia from 400,000 years ago until they went extinct around 40,000 years ago. DNA retrieved from ancient specimens revealed that Neanderthals mated with modern human contemporaries. As a consequence, introgressed Neanderthal DNA survives scattered across the human genome such that 1-4% of the genome of present-day people outside Africa are inherited from Neanderthal ancestors. Patterns of Neanderthal introgressed genomic sequences suggest that Neanderthal alleles had distinct fates in the modern human genetic background. Some Neanderthal alleles facilitated human adaptation to new environments such as novel climate conditions, UV exposure levels and pathogens, while others had deleterious consequences. Here, we review the body of work on Neanderthal introgression over the past decade. We describe how evolutionary forces shaped the genomic landscape of Neanderthal introgression and highlight the impact of introgressed alleles on human biology and phenotypic variation.

Decoding the Human Face: Progress and Challenges in Understanding the Genetics of Craniofacial Morphology

Variations in the form of the human face, which plays a role in our individual identities and societal interactions, have fascinated scientists and artists alike. Here, we review our current understanding of the genetics underlying variation in craniofacial morphology and disease-associated dysmorphology, synthesizing decades of progress on Mendelian syndromes in addition to more recent results from genome-wide association studies of human facial shape and disease risk. We also discuss the various approaches used to phenotype and quantify facial shape, which are of particular importance due to the complex, multipartite nature of the craniofacial form. We close by discussing how experimental studies have contributed and will further contribute to our understanding of human genetic variation and then proposing future directions and applications for the field.

A Middle Pleistocene Denisovan molar from the Annamite Chain of northern Laos

The Pleistocene presence of the genus Homo in continental Southeast Asia is primarily evidenced by a sparse stone tool record and rare human remains. Here we report a Middle Pleistocene hominin specimen from Laos, with the discovery of a molar from the Tam Ngu Hao 2 (Cobra Cave) limestone cave in the Annamite Mountains. The age of the fossil-bearing breccia ranges between 164–131 kyr, based on the Bayesian modelling of luminescence dating of the sedimentary matrix from which it was recovered, U-series dating of an overlying flowstone, and U-series–ESR dating of associated faunal teeth. Analyses of the internal structure of the molar in tandem with palaeoproteomic analyses of the enamel indicate that the tooth derives from a young, likely female, Homo individual. The close morphological affinities with the Xiahe specimen from China indicate that they belong to the same taxon and that Tam Ngu Hao 2 most likely represents a Denisovan.

Evidence for the presence of Homo during the Middle Pleistocene is limited in continental Southeast Asia. Here, the authors report a hominin molar from Tam Ngu Hao 2 (Cobra Cave), dated to 164–131 kyr. They use morphological and paleoproteomic analysis to show that it likely belonged to a female Denisovan.

Predicting Archaic Hominin Phenotypes from Genomic Data

Ancient DNA provides a powerful window into the biology of extant and extinct species, including humans' closest relatives: Denisovans and Neanderthals. Here, we review what is known about archaic hominin phenotypes from genomic data and how those inferences have been made. We contend that understanding the influence of variants on lower-level molecular phenotypes-such as gene expression and protein function-is a promising approach to using ancient DNA to learn about archaic hominin traits. Molecular phenotypes have simpler genetic architectures than organism-level complex phenotypes, and this approach enables moving beyond association studies by proposing hypotheses about the effects of archaic variants that are testable in model systems. The major challenge to understanding archaic hominin phenotypes is broadening our ability to accurately map genotypes to phenotypes, but ongoing advances ensure that there will be much more to learn about archaic hominin phenotypes from their genomes. Expected final online publication date for the Annual Review of Genomics and Human Genetics, Volume 23 is October 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Evolution of cranial capacity revisited: A view from the late Middle Pleistocene cranium from Xujiayao, China

The Late Middle Pleistocene hominin fossils from the Xujiayao site in northern China have been closely studied in light of their morphological variability. However, all previous studies have focused on separated cranial fragments. Here, we report the first reconstruction of a fairly complete posterior cranium, Xujiayao 6 (XJY 6), confidently dated to ∼200-160 ka, which facilitated an assessment of its overall cranial size. XJY 6 was reconstructed from three of the original fragments-the PA1486 (No.7/XJY 6a) occipital bone, PA1490 (No.10/XJY 6b) right parietal bone, and PA1498 (No.17/XJY 15) left temporal bone-which originated from the same young adult individual. The XJY 6 endocranial capacity, estimated by measuring endocranial volume, was estimated using multiple regression formulae derived from ectocranial and endocranial measurements on select samples of Pleistocene hominins and recent modern humans. The results indicate that the larger pooled sample of both Pleistocene and recent modern humans was more robust for the endocranial capacity estimate. Based on the pooled sample using the ectocranial and endocranial measurements, we conservatively estimate the XJY 6 endocranial volume to be ∼1700 cm³ with a 95% confidence interval of 1555-1781 cm³. This is close to Xuchang 1, which dates to 125-105 ka and whose endocranial volume is ∼1800 cm³. Thus, XJY 6 provides the earliest evidence of a brain size that falls in the upper range of Neanderthals and modern Homo sapiens. XJY 6, together with Xuchang 1, Homo floresiensis, Homo luzonensis, and Homo naledi, challenge the general pattern that brain size gradually increases over geological time. This study also finds that hominin brain size expansion occurred at different rates across time and space.

Human brain evolution: Emerging roles for regulatory DNA and RNA

Humans diverge from other primates in numerous ways, including their neuroanatomy and cognitive capacities. Human-specific features are particularly prominent in the cerebral cortex, which has undergone an expansion in size and acquired unique cellular composition and circuitry. Human-specific gene expression is postulated to explain neocortical anatomical differences across evolution. In particular, noncoding regulatory loci are strongly linked to human traits, including progenitor proliferation and cortical size. In this review, we highlight emerging noncoding elements implicated in human cortical evolution, including roles for regulatory DNA and RNA. Further, we discuss the association of human-specific genetic changes with neurodevelopmental diseases.

Progress and Prospects in Epigenetic Studies of Ancient DNA

Development of technologies for high-throughput whole-genome sequencing and improvement of sample preparation techniques made it possible to study ancient DNA (aDNA) from archaeological samples over a million year old. The studies of aDNA have shed light on the history of human migration, replacement of populations, interbreeding of Cro-Magnons with Neanderthals and Denisovans, evolution of human pathogens, etc. Equally important is the possibility to investigate epigenetic modifications of ancient genomes, which has allowed to obtain previously inaccessible information on gene expression, nucleosome positioning, and DNA methylation. Analysis of methylation status of certain genomic sites can predict an individual's age at death and reconstruct some phenotypic features, as it was done for the Denisovan genome, and even to elucidate unfavorable environmental factors that had affected this archaic individual. In this review, we discuss current progress in epigenetic studies of aDNA, including methodological approaches and promising research directions in this field.

The cis-regulatory effects of modern human-specific variants

The Neanderthal and Denisovan genomes enabled the discovery of sequences that differ between modern and archaic humans, the majority of which are noncoding. However, our understanding of the regulatory consequences of these differences remains limited, in part due to the decay of regulatory marks in ancient samples. Here, we used a massively parallel reporter assay in embryonic stem cells, neural progenitor cells, and bone osteoblasts to investigate the regulatory effects of the 14,042 single-nucleotide modern human-specific variants. Overall, 1791 (13%) of sequences containing these variants showed active regulatory activity, and 407 (23%) of these drove differential expression between human groups. Differentially active sequences were associated with divergent transcription factor binding motifs, and with genes enriched for vocal tract and brain anatomy and function. This work provides insight into the regulatory function of variants that emerged along the modern human lineage and the recent evolution of human gene expression.

Human-chimpanzee fused cells reveal cis-regulatory divergence underlying skeletal evolution

Gene regulatory divergence is thought to play a central role in determining human-specific traits. However, our ability to link divergent regulation to divergent phenotypes is limited. Here, we utilized human-chimpanzee hybrid induced pluripotent stem cells to study gene expression separating these species. The tetraploid hybrid cells allowed us to separate cis- from trans-regulatory effects, and to control for non-genetic confounding factors. We differentiated these cells into cranial neural crest cells (CNCCs), the primary cell type giving rise to the face. We discovered evidence of lineage-specific selection on the hedgehog signaling pathway, including a human-specific 6-fold down-regulation of EVC2 (LIMBIN), a key hedgehog gene. Inducing a similar down-regulation of EVC2 substantially reduced hedgehog signaling output. Mice and humans lacking functional EVC2 show striking phenotypic parallels to human-chimpanzee craniofacial differences, suggesting that the regulatory divergence of hedgehog signaling may have contributed to the unique craniofacial morphology of humans.

Visual Depictions of Our Evolutionary Past: A Broad Case Study Concerning the Need for Quantitative Methods of Soft Tissue Reconstruction and Art-Science Collaborations

Flip through scientific textbooks illustrating ideas about human evolution or visit any number of museums of natural history and you will notice an abundance of reconstructions attempting to depict the appearance of ancient hominins. Spend some time comparing reconstructions of the same specimen and notice an obvious fact: hominin reconstructions vary in appearance considerably. In this review, we summarize existing methods of reconstruction to analyze this variability. It is argued that variability between hominin reconstructions is likely the result of unreliable reconstruction methods and misinterpretation of available evidence. We also discuss the risk of disseminating erroneous ideas about human evolution through the use of unscientific reconstructions in museums and publications. The role an artist plays is also analyzed and criticized given how the aforementioned reconstructions have become readily accepted to line the halls of even the most trusted institutions. In conclusion, improved reconstruction methods hold promise for the prediction of hominin soft tissues, as well as for disseminating current scientific understandings of human evolution in the future.

Harnessing epigenetics to study human evolution

Recent advances in ancient DNA extraction and high-throughput sequencing technologies enabled the high-quality sequencing of archaic genomes, including the Neanderthal and the Denisovan. While comparisons with modern humans revealed both archaic-specific and human-specific sequence changes, in the absence of gene expression information, understanding the functional implications of such genetic variations remains a major challenge. To study gene regulation in archaic humans, epigenetic research comes to our aid. DNA methylation, which is highly correlated with transcription, can be directly measured in modern samples, as well as reconstructed in ancient samples. This puts DNA methylation as a natural basis for comparative epigenetics between modern humans, archaic humans and nonhuman primates.

Archaic hominin genomics provides a window into gene expression evolution

Differences in gene expression are thought to account for most phenotypic differences within and between species. Consequently, gene expression is a powerful lens through which to study divergence between modern humans and our closest evolutionary relatives, the Neanderthals and Denisovans. Such insights complement biological knowledge gleaned from the fossil record, while also revealing general features of the mode and tempo of regulatory evolution. Because of the degradation of ancient RNA, gene expression profiles of archaic hominins must be studied by indirect means. As such, conclusions drawn from these studies are often laden with assumptions about the genetic architecture of gene expression, the complexity of which is increasingly apparent. Despite these challenges, rapid technical and conceptual advances in the fields of ancient genomics, functional genomics, statistical genomics, and genome engineering are revolutionizing understanding of hominin gene expression evolution.

Differential DNA methylation of vocal and facial anatomy genes in modern humans

Changes in potential regulatory elements are thought to be key drivers of phenotypic divergence. However, identifying changes to regulatory elements that underlie human-specific traits has proven very challenging. Here, we use 63 reconstructed and experimentally measured DNA methylation maps of ancient and present-day humans, as well as of six chimpanzees, to detect differentially methylated regions that likely emerged in modern humans after the split from Neanderthals and Denisovans. We show that genes associated with face and vocal tract anatomy went through particularly extensive methylation changes. Specifically, we identify widespread hypermethylation in a network of face- and voice-associated 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.

More Arrows in the Ancient DNA Quiver: Use of Paleoepigenomes and Paleomicrobiomes to Investigate Animal Adaptation to Environment

Whether and how epigenetic mechanisms and the microbiome play a role in mammalian adaptation raised considerable attention and controversy, mainly because they have the potential to add new insights into the Modern Synthesis. Recent attempts to reconcile neo-Darwinism and neo-Lamarckism in a unified theory of molecular evolution give epigenetic mechanisms and microbiome a prominent role. However, supporting empirical data are still largely missing. Because experimental studies using extant animals can hardly be done over evolutionary timescales, we propose that advances in ancient DNA techniques provide a valid alternative. In this piece, we evaluate 1) the possible roles of epigenomes and microbiomes in animal adaptation, 2) advances in the retrieval of paleoepigenome and paleomicrobiome data using ancient DNA techniques, and 3) the plasticity of either and interactions between the epigenome and the microbiome, while emphasizing that it is essential to take both into account, as well as the underlying genetic factors that may confound the findings. We propose that advanced ancient DNA techniques should be applied to a wide range of past animals, so novel dynamics in animal evolution and adaption can be revealed.

Predicted Archaic 3D Genome Organization Reveals Genes Related to Head and Spinal Cord Separating Modern from Archaic Humans

High coverage sequences of archaic humans enabled the reconstruction of their DNA methylation patterns. This allowed comparing gene regulation between human groups, and linking such regulatory changes to phenotypic differences. In a previous work, a detailed comparison of DNA methylation in modern humans, archaic humans, and chimpanzees revealed 873 modern human-derived differentially methylated regions (DMRs). To understand the regulatory implications of these DMRs, we defined differentially methylated genes (DMGs) as genes that harbor DMRs in their promoter or gene body. While most of the modern human-derived DMRs could be linked to DMGs, many others remained unassigned. Here, we used information on 3D genome organization to link ~70 out of the remaining 288 unassigned DMRs to genes. Combined with the previously identified DMGs, we reinforce the enrichment of these genes with vocal and facial anatomy, and additionally find significant enrichment with the spinal column, chin, hair, and scalp. These results reveal the importance of 3D genomic organization in understanding gene regulation by DNA methylation.