Modern humans did not admix with Neanderthals during their range expansion into Europe

The process by which the Neanderthals were replaced by modern humans between 42,000 and 30,000 before present is still intriguing. Although no Neanderthal mitochondrial DNA (mtDNA) lineage is found to date among several thousands of Europeans and in seven early modern Europeans, interbreeding rates as high as 25% could not be excluded between the two subspecies. In this study, we introduce a realistic model of the range expansion of early modern humans into Europe, and of their competition and potential admixture with local Neanderthals. Under this scenario, which explicitly models the dynamics of Neanderthals' replacement, we estimate that maximum interbreeding rates between the two populations should have been smaller than 0.1%. We indeed show that the absence of Neanderthal mtDNA sequences in Europe is compatible with at most 120 admixture events between the two populations despite a likely cohabitation time of more than 12,000 y. This extremely low number strongly suggests an almost complete sterility between Neanderthal females and modern human males, implying that the two populations were probably distinct biological species.     

A Haplotype at STAT2 Introgressed from Neanderthals and Serves as a Candidate of Positive Selection in Papua New Guinea

Signals of archaic admixture have been identified through comparisons of the draft Neanderthal and Denisova genomes with those of living humans. Studies of individual loci contributing to these genome-wide average signals are required for characterization of the introgression process and investigation of whether archaic variants conferred an adaptive advantage to the ancestors of contemporary human populations. However, no definitive case of adaptive introgression has yet been described. Here we provide a DNA sequence analysis of the innate immune gene STAT2 and show that a haplotype carried by many Eurasians (but not sub-Saharan Africans) has a sequence that closely matches that of the Neanderthal STAT2. This haplotype, referred to as N, was discovered through a resequencing survey of the entire coding region of STAT2 in a global sample of 90 individuals. Analyses of publicly available complete genome sequence data show that haplotype N shares a recent common ancestor with the Neanderthal sequence (∼80 thousand years ago) and is found throughout Eurasia at an average frequency of ∼5%. Interestingly, N is found in Melanesian populations at ∼10-fold higher frequency (∼54%) than in Eurasian populations. A neutrality test that controls for demography rejects the hypothesis that a variant of N rose to high frequency in Melanesia by genetic drift alone. Although we are not able to pinpoint the precise target of positive selection, we identify nonsynonymous mutations in ERBB3, ESYT1, and STAT2—all of which are part of the same 250 kb introgressive haplotype—as good candidates.     

The Date of Interbreeding between Neandertals and Modern Humans

Comparisons of DNA sequences between Neandertals and present-day humans have shown that Neandertals share more genetic variants with non-Africans than with Africans. This could be due to interbreeding between Neandertals and modern humans when the two groups met subsequent to the emergence of modern humans outside Africa. However, it could also be due to population structure that antedates the origin of Neandertal ancestors in Africa. We measure the extent of linkage disequilibrium (LD) in the genomes of present-day Europeans and find that the last gene flow from Neandertals (or their relatives) into Europeans likely occurred 37,000–86,000 years before the present (BP), and most likely 47,000–65,000 years ago. This supports the recent interbreeding hypothesis and suggests that interbreeding may have occurred when modern humans carrying Upper Paleolithic technologies encountered Neandertals as they expanded out of Africa.     

Ancient structure in Africa unlikely to explain neanderthal and non-african genetic similarity

Neanderthals have been shown to share more genetic variants with present-day non-Africans than Africans. Recent admixture between Neanderthals and modern humans outside of Africa was proposed as the most parsimonious explanation for this observation. However, the hypothesis of ancient population structure within Africa could not be ruled out as an alternative explanation. We use simulations to test whether the site frequency spectrum, conditioned on a derived Neanderthal and an ancestral Yoruba (African) nucleotide (the doubly conditioned site frequency spectrum [dcfs]), can distinguish between models that assume recent admixture or ancient population structure. We compare the simulations to the dcfs calculated from data taken from populations of European, Chinese, and Japanese descent in the Complete Genomics Diversity Panel. Simulations under a variety of plausible demographic parameters were used to examine the shape of the dcfs for both models. The observed shape of the dcfs cannot be explained by any set of parameter values used in the simulations of the ancient structure model. The dcfs simulations for the recent admixture model provide a good fit to the observed dcfs for non-Africans, thereby supporting the hypothesis that recent admixture with Neanderthals accounts for the greater similarity of Neanderthals to non-Africans than Africans.     

Dynamics between Cancer Cell Subpopulations Reveals a Model Coordinating with Both Hierarchical and Stochastic Concepts

Tumors are often heterogeneous in which tumor cells of different phenotypes have distinct properties. For scientific and clinical interests, it is of fundamental importance to understand their properties and the dynamic variations among different phenotypes, specifically under radio- and/or chemo-therapy. Currently there are two controversial models describing tumor heterogeneity, the cancer stem cell (CSC) model and the stochastic model. To clarify the controversy, we measured probabilities of different division types and transitions of cells via in situ immunofluorescence. Based on the experiment data, we constructed a model that combines the CSC with the stochastic concepts, showing the existence of both distinctive CSC subpopulations and the stochastic transitions from NSCCs to CSCs. The results showed that the dynamic variations between CSCs and non-stem cancer cells (NSCCs) can be simulated with the model. Further studies also showed that the model can be used to describe the dynamics of the two subpopulations after radiation treatment. More importantly, analysis demonstrated that the experimental detectable equilibrium CSC proportion can be achieved only when the stochastic transitions from NSCCs to CSCs occur, indicating that tumor heterogeneity may exist in a model coordinating with both the CSC and the stochastic concepts. The mathematic model based on experimental parameters may contribute to a better understanding of the tumor heterogeneity, and provide references on the dynamics of CSC subpopulation during radiotherapy.     

Coalescence Times and FST Under a Skewed Offspring Distribution Among Individuals in a Population

Estimates of gene flow between subpopulations based on F_ST (or N_ST) are shown to be confounded by the reproduction parameters of a model of skewed offspring distribution. Genetic evidence of population subdivision can be observed even when gene flow is very high, if the offspring distribution is skewed. A skewed offspring distribution arises when individuals can have very many offspring with some probability. This leads to high probability of identity by descent within subpopulations and results in genetic heterogeneity between subpopulations even when Nm is very large. Thus, we consider a limiting model in which the rates of coalescence and migration can be much higher than for a Wright–Fisher population. We derive the densities of pairwise coalescence times and expressions for F_ST and other statistics under both the finite island model and a many-demes limit model. The results can explain the observed genetic heterogeneity among subpopulations of certain marine organisms despite substantial gene flow.     

Archaic admixture in the human genome

One of the enduring questions in the evolution of our species surrounds the fate of 'archaic' forms of Homo. Did Neanderthals go extinct without interbreeding with modern humans 25-40 thousand years ago or are their genes present among modern-day Europeans? Recent work suggests that Neanderthals and an as yet unidentified archaic African population contributed to at least 5% of the modern European and West African gene pools, respectively. Extensive sequencing of Neanderthal and other archaic human nuclear DNA has the potential to answer this question definitively within the next few years.     

Evaluating sample allocation and effort in detecting population differentiation for discrete and continuously distributed individuals

One of the most pressing issues in spatial genetics concerns sampling. Traditionally, substructure and gene flow are estimated for individuals sampled within discrete populations. Because many species may be continuously distributed across a landscape without discrete boundaries, understanding sampling issues becomes paramount. Given large-scale, geographically broad conservation efforts, researchers are looking for guidance as to the trade-offs between sampling more individuals within a population versus few individuals scattered across more populations. Here, we conducted simulations that address these issues. We first established two archetypical patterns of dispersion: (1) individuals within discrete populations, and (2) continuously distributed individuals with limited dispersal. We used genotypes generated from a spatially-explicit, individual-based program and simulated genetic structure in individuals from nine different population sizes across a landscape that either had barriers to movement (defining discrete populations) or isolation-by-distance patterns (defining continuously distributed individuals). Then, given each pattern of dispersion, we allocated samples across four different sampling strategies for each of the nine population sizes in various configurations for sampling more individuals within a population versus fewer individuals scattered across more populations. We assessed the population genetic substructure with both the population-based metric, F _ST, and an individual-based metric, D _PS regardless of the true pattern of dispersion to allow us to better understand the effect of incorrectly matching the metric and the distribution (e.g., F _ST with continuously distributed individuals, and vice versa). We show that sampling many subpopulations (or sampling areas), thus sampling fewer individuals per subpopulation, overestimates measures of population subdivision with the population-based metric for both patterns of dispersion. In contrast, using the individual-based metric gives the opposite results: sampling too few subpopulations, and many individuals per subpopulation, produces an underestimate of the strength of isolation-by-distance. By comparing all results, we were able to suggest a strong predictive model of a chosen genetic structure metric for elucidating the sampling design trade-offs given each pattern of dispersion and configuration on the landscape.     

Evaluating Neanderthal Genetics and Phylogeny

The retrieval of Neanderthal (Homo neanderthalsensis) mitochondrial DNA is thought to be among the most significant ancient DNA contributions to date, allowing conflicting hypotheses on modern human (Homo sapiens) evolution to be tested directly. Recently, however, both the authenticity of the Neanderthal sequences and their phylogenetic position outside contemporary human diversity have been questioned. Using Bayesian inference and the largest dataset to date, we find strong support for a monophyletic Neanderthal clade outside the diversity of contemporary humans, in agreement with the expectations of the Out-of-Africa replacement model of modern human origin. From average pairwise sequence differences, we obtain support for claims that the first published Neanderthal sequence may include errors due to postmortem damage in the template molecules for PCR. In contrast, we find that recent results implying that the Neanderthal sequences are products of PCR artifacts are not well supported, suffering from inadequate experimental design and a presumably high percentage (>68%) of chimeric sequences due to “jumping PCR” events. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Martin Kreitman]     

A deterministic model of admixture and genetic introgression: the case of Neanderthal and Cro-Magnon

There is an ongoing debate in the field of human evolution about the possible contribution of Neanderthals to the modern human gene pool. To study how the Neanderthal private alleles may have spread over the genes of Homo sapiens, we propose a deterministic model based on recursive equations and ordinary differential equations. If the Neanderthal population was large compared to the Homo sapiens population at the beginning of the contact period, we show that genetic introgression should have been fast and complete meaning that most of the Neanderthal private alleles should be found in the modern human gene pool in case of ancient admixture. In order to test/reject ancient admixture from genome-wide data, we incorporate the model of genetic introgression into a statistical hypothesis-testing framework. We show that the power to reject ancient admixture increases as the ratio, at the time of putative admixture, of the population size of Homo sapiens over that of Neanderthal decreases. We find that the power to reject ancient admixture might be particularly low if the population size of Homo sapiens was comparable to the Neanderthal population size.     

Genetic Origins of Lactase Persistence and the Spread of Pastoralism in Africa

In humans, the ability to digest lactose, the sugar in milk, declines after weaning because of decreasing levels of the enzyme lactase-phlorizin hydrolase, encoded by LCT. However, some individuals maintain high enzyme amounts and are able to digest lactose into adulthood (i.e., they have the lactase-persistence [LP] trait). It is thought that selection has played a major role in maintaining this genetically determined phenotypic trait in different human populations that practice pastoralism. To identify variants associated with the LP trait and to study its evolutionary history in Africa, we sequenced MCM6 introns 9 and 13 and ∼2 kb of the LCT promoter region in 819 individuals from 63 African populations and in 154 non-Africans from nine populations. We also genotyped four microsatellites in an ∼198 kb region in a subset of 252 individuals to reconstruct the origin and spread of LP-associated variants in Africa. Additionally, we examined the association between LP and genetic variability at candidate regulatory regions in 513 individuals from eastern Africa. Our analyses confirmed the association between the LP trait and three common variants in intron 13 (C-14010, G-13907, and G-13915). Furthermore, we identified two additional LP-associated SNPs in intron 13 and the promoter region (G-12962 and T-956, respectively). Using neutrality tests based on the allele frequency spectrum and long-range linkage disequilibrium, we detected strong signatures of recent positive selection in eastern African populations and the Fulani from central Africa. In addition, haplotype analysis supported an eastern African origin of the C-14010 LP-associated mutation in southern Africa.     

DNA fingerprinting analysis of subpopulations of the red squirrel (Sciurus vulgaris, L. 1758)

The relationships within two squirrel groups of 7 and 10 individuals living in small isolated forest areas (»20 ha) were studied by DNA fingerprinting. The study area was located in northwestern Germany, close to the town of Osnabrück. The degrees of relationship determined by the multi locus probes 33.15 and MZ1.3 were >20–40% higher in the two subpopulations compared with a group of unrelated squirrels. The band-sharing rates in the latter group varied between 6 and 57% (33.15) and 17 and 63% (MZ1.3), whereas the squirrels living in the two forest islands had band-sharing rates between 13 and 80%, depending on the probe used. Although the forest islands were isolated by fields and roads, there seemed to be a small genetic exchange between the squirrels of the forest islands and those of the next adjacent larger woodland area.     

The effect of trauma on Neanderthal culture: A mathematical analysis

Traumatic lesions are often observed in ancient skeletal remains. Since ancient medical technology was immature, severely traumatized individuals may have frequently lost the physical ability for cultural skills that demand complex body movements. I develop a mathematical model to analyze the effect of trauma on cultural transmission and apply it to Neanderthal culture using Neanderthal fossil data. I estimate from the data that the proportion of adult individuals who suffered traumatic injuries before death was approximately 0.79–0.94, in which 0.37–0.52 were injured severely and 0.13–0.19 were injured before adulthood. Assuming that every severely traumatized individual and a quarter to a half of the other traumatized individuals lost the capacity for a cultural skill that demands complex control of the traumatized body part, I estimate that if an upper limb is associated with a cultural skill, each individual had to communicate closely with at least 1.5–2.6 individuals during adulthood to maintain the skill in Neanderthal society, and if a whole body is associated, at least 3.1–11.5 individuals were necessary. If cultural transmissions between experts and novices were inaccurate, or if low frequency skills easily disappeared from the population due to random drift, more communicable individuals were necessary. Since the community size of Neanderthals was very small, their high risk of injury may have inhibited the spread of technically difficult cultural skills in their society. It may be important to take this inhibition into consideration when we study Neanderthal culture and the replacement of Neanderthals by modern humans.     

Genome digging: insight into the mitochondrial genome of Homo

Background

A fraction of the Neanderthal mitochondrial genome sequence has a similarity with a 5,839-bp nuclear DNA sequence of mitochondrial origin (numt) on the human chromosome 1. This fact has never been interpreted. Although this phenomenon may be attributed to contamination and mosaic assembly of Neanderthal mtDNA from short sequencing reads, we explain the mysterious similarity by integration of this numt (mtAncestor-1) into the nuclear genome of the common ancestor of Neanderthals and modern humans not long before their reproductive split.

Principal Findings

Exploiting bioinformatics, we uncovered an additional numt (mtAncestor-2) with a high similarity to the Neanderthal mtDNA and indicated that both numts represent almost identical replicas of the mtDNA sequences ancestral to the mitochondrial genomes of Neanderthals and modern humans. In the proteins, encoded by mtDNA, the majority of amino acids distinguishing chimpanzees from humans and Neanderthals were acquired by the ancestral hominins. The overall rate of nonsynonymous evolution in Neanderthal mitochondrial protein-coding genes is not higher than in other lineages. The model incorporating the ancestral hominin mtDNA sequences estimates the average divergence age of the mtDNAs of Neanderthals and modern humans to be 450,000–485,000 years. The mtAncestor-1 and mtAncestor-2 sequences were incorporated into the nuclear genome approximately 620,000 years and 2,885,000 years ago, respectively.

Conclusions

This study provides the first insight into the evolution of the mitochondrial DNA in hominins ancestral to Neanderthals and humans. We hypothesize that mtAncestor-1 and mtAncestor-2 are likely to be molecular fossils of the mtDNAs of Homo heidelbergensis and a stem Homo lineage. The d_N/d_S dynamics suggests that the effective population size of extinct hominins was low. However, the hominin lineage ancestral to humans, Neanderthals and H. heidelbergensis, had a larger effective population size and possessed genetic diversity comparable with those of chimpanzee and gorilla.     

Neanderthal and Denisova genetic affinities with contemporary humans: Introgression versus common ancestral polymorphisms

Analyses of the genetic relationships among modern humans, Neanderthals and Denisovans have suggested that 1–4% of the non-Sub-Saharan African gene pool may be Neanderthal derived, while 6–8% of the Melanesian gene pool may be the product of admixture between the Denisovans and the direct ancestors of Melanesians. In the present study, we analyzed single nucleotide polymorphism (SNP) diversity among a worldwide collection of contemporary human populations with respect to the genetic constitution of these two archaic hominins and Pan troglodytes (chimpanzee). We partitioned SNPs into subsets, including those that are derived in both archaic lineages, those that are ancestral in both archaic lineages and those that are only derived in one archaic lineage. By doing this, we have conducted separate examinations of subsets of mutations with higher probabilities of divergent phylogenetic origins. While previous investigations have excluded SNPs from common ancestors in principal component analyses, we included common ancestral SNPs in our analyses to visualize the relative placement of the Neanderthal and Denisova among human populations. To assess the genetic similarities among the various hominin lineages, we performed genetic structure analyses to provide a comparison of genetic patterns found within contemporary human genomes that may have archaic or common ancestral roots. Our results indicate that 3.6% of the Neanderthal genome is shared with roughly 65.4% of the average European gene pool, which clinally diminishes with distance from Europe. Our results suggest that Neanderthal genetic associations with contemporary non-Sub-Saharan African populations, as well as the genetic affinities observed between Denisovans and Melanesians most likely result from the retention of ancient mutations in these populations.     

Development of genetic differentiation and postzygotic isolation in experimental metapopulations of spider mites

We studied the development of genetic differentiation and postzygotic isolation in experimental metapopulations of the two-spotted spider mite, Tetranychus urticae Koch. A genetically diverse starter population was made by allowing six inbred sublines to interbreed. Then three migration patterns were tested: no migration, or one or three immigrants per subpopulation per generation. Variations in four traits were investigated: allozymes, acaricide resistance, diapause, and hatchability. In the allozymes, acaricide resistance, and diapause, genetic variation among subpopulations became high in metapopulations with no migration, but not in the others, which showed that one immigrant is enough to prevent genetic differentiation. Hatchability, which was decreased by interbreeding among the six sublines, gradually recovered in succeeding generations. In metapopulations with no migration, hatchability was reduced again after in-migration at the 15th generation. Different karyotypes or coadapted gene complexes can survive in different subpopulations by genetic drift, and both Wolbachia-infected and -noninfected subpopulations may be selected, which would lead to postzygotic isolation between isolated subpopulations. Our results indicate that sampling effects such as genetic drift or stochastic loss of Wolbachia produce postzygotic isolation in laboratory populations of spider mite.     

Strongyloides stercoralis genotypes in humans in Cambodia

Little is known about the genetic variability of the soil-transmitted nematode, Strongyloides stercoralis, in humans. We sequenced portions of the small subunit rDNA (SSU), including the hyper variable regions (HVR) I and IV from S. stercoralis larvae derived from individuals living in a rural setting in Cambodia. We identified three polymorphic positions, including a previously reported one within the HVR I. HVR IV was invariable. Six different SSU alleles existed in our sample. Although different genotypes of S. stercoralis were found in the same individuals, no heterozygous larvae were found. This indicates that there is no or very little interbreeding between the different genotypes. Further studies are needed to examine if this is because sexual reproduction, which is facultative, is rare in our study area's S. stercoralis population or because what is considered to be S. stercoralis today is actually a complex of closely related species or subspecies.     

Comparison of dental measurement systems for taxonomic assignment of Neanderthal and modern human lower second deciduous molars

Traditional morphometric approaches for taxonomic assignment of Neanderthal and modern human dental remains are mainly characterized by caliper measurements of tooth crowns. Several studies have recently described differences in dental tissue proportions and enamel thickness between Neanderthal and modern human teeth. At least for the lower second deciduous molar (dm₂), a three-dimensional lateral relative enamel thickness index has been proposed for separating the two taxa. This index has the advantage over other measurements of being applicable to worn teeth because it ignores the occlusal aspect of the crown. Nevertheless, a comparative evaluation of traditional crown dimensions and lateral dental tissue proportion measurements for taxonomic assignment of Neanderthal and modern human dm₂s has not yet been performed.In this study, we compare various parameters gathered from the lateral aspects of the crown. These parameters include crown diameters, height of the lateral wall of the crown (lateral crown height = LCH), lateral enamel thickness, and dentine volume of the lateral wall, including the volume of the coronal pulp chamber (lateral dentine plus pulp volume = LDPV), in a 3D digital sample of Neanderthal and modern human dm₂s to evaluate their utility in separating the two taxa.The LDPV and the LCH allow us to discriminate between Neanderthals and modern humans with 88.5% and 92.3% accuracy, respectively. Though our results confirm that Neanderthal dm₂s have lower relative enamel thickness (RET) index compared with modern humans (p = 0.005), only 70% of the specimens were correctly classified on the basis of the RET index. We also emphasize that results of the lateral enamel thickness method depend on the magnitude of the interproximal wear. Accordingly, we suggest using the LCH or the LDPV to discriminate between Neanderthal and modern human dm₂s. These parameters are more independent of interproximal wear and loss of lateral enamel.