Genetic insights into the social organization of Neanderthals.

Genomic analyses of Neanderthals have previously provided insights into their population history and relationship to modern humans1-8, but the social organization of Neanderthal communities remains poorly understood. Here we present genetic data for 13 Neanderthals from two Middle Palaeolithic sites in the Altai Mountains of southern Siberia: 11 from Chagyrskaya Cave9,10 and 2 from Okladnikov Cave11-making this one of the largest genetic studies of a Neanderthal population to date. We used hybridization capture to obtain genome-wide nuclear data, as well as mitochondrial and Y-chromosome sequences. Some Chagyrskaya individuals were closely related, including a father-daughter pair and a pair of second-degree relatives, indicating that at least some of the individuals lived at the same time. Up to one-third of these individuals' genomes had long segments of homozygosity, suggesting that the Chagyrskaya Neanderthals were part of a small community. In addition, the Y-chromosome diversity is an order of magnitude lower than the mitochondrial diversity, a pattern that we found is best explained by female migration between communities. Thus, the genetic data presented here provide a detailed documentation of the social organization of an isolated Neanderthal community at the easternmost extent of their known range.

Initial Upper Palaeolithic humans in Europe had recent Neanderthal ancestry.

Modern humans appeared in Europe by at least 45,000 years ago1-5, but the extent of their interactions with Neanderthals, who disappeared by about 40,000 years ago6, and their relationship to the broader expansion of modern humans outside Africa are poorly understood. Here we present genome-wide data from three individuals dated to between 45,930 and 42,580 years ago from Bacho Kiro Cave, Bulgaria1,2. They are the earliest Late Pleistocene modern humans known to have been recovered in Europe so far, and were found in association with an Initial Upper Palaeolithic artefact assemblage. Unlike two previously studied individuals of similar ages from Romania7 and Siberia8 who did not contribute detectably to later populations, these individuals are more closely related to present-day and ancient populations in East Asia and the Americas than to later west Eurasian populations. This indicates that they belonged to a modern human migration into Europe that was not previously known from the genetic record, and provides evidence that there was at least some continuity between the earliest modern humans in Europe and later people in Eurasia. Moreover, we find that all three individuals had Neanderthal ancestors a few generations back in their family history, confirming that the first European modern humans mixed with Neanderthals and suggesting that such mixing could have been common.

Microstratigraphic preservation of ancient faunal and hominin DNA in Pleistocene cave sediments.

Ancient DNA recovered from Pleistocene sediments represents a rich resource for the study of past hominin and environmental diversity. However, little is known about how DNA is preserved in sediments and the extent to which it may be translocated between archaeological strata. Here, we investigate DNA preservation in 47 blocks of resin-impregnated archaeological sediment collected over the last four decades for micromorphological analyses at 13 prehistoric sites in Europe, Asia, Africa, and North America and show that such blocks can preserve DNA of hominins and other mammals. Extensive microsampling of sediment blocks from Denisova Cave in the Altai Mountains reveals that the taxonomic composition of mammalian DNA differs drastically at the millimeter-scale and that DNA is concentrated in small particles, especially in fragments of bone and feces (coprolites), suggesting that these are substantial sources of DNA in sediments. Three microsamples taken in close proximity in one of the blocks yielded Neanderthal DNA from at least two male individuals closely related to Denisova 5, a Neanderthal toe bone previously recovered from the same layer. Our work indicates that DNA can remain stably localized in sediments over time and provides a means of linking genetic information to the archaeological and ecological records on a microstratigraphic scale.

A humanized version of Foxp2 affects cortico-basal ganglia circuits in mice.

It has been proposed that two amino acid substitutions in the transcription factor FOXP2 have been positively selected during human evolution due to effects on aspects of speech and language. Here, we introduce these substitutions into the endogenous Foxp2 gene of mice. Although these mice are generally healthy, they have qualitatively different ultrasonic vocalizations, decreased exploratory behavior and decreased dopamine concentrations in the brain suggesting that the humanized Foxp2 allele affects basal ganglia. In the striatum, a part of the basal ganglia affected in humans with a speech deficit due to a nonfunctional FOXP2 allele, we find that medium spiny neurons have increased dendrite lengths and increased synaptic plasticity. Since mice carrying one nonfunctional Foxp2 allele show opposite effects, this suggests that alterations in cortico-basal ganglia circuits might have been important for the evolution of speech and language in humans.

Reconstructing the genetic history of late Neanderthals.

Although it has previously been shown that Neanderthals contributed DNA to modern humans, not much is known about the genetic diversity of Neanderthals or the relationship between late Neanderthal populations at the time at which their last interactions with early modern humans occurred and before they eventually disappeared. Our ability to retrieve DNA from a larger number of Neanderthal individuals has been limited by poor preservation of endogenous DNA and contamination of Neanderthal skeletal remains by large amounts of microbial and present-day human DNA. Here we use hypochlorite treatment of as little as 9 mg of bone or tooth powder to generate between 1- and 2.7-fold genomic coverage of five Neanderthals who lived around 39,000 to 47,000 years ago (that is, late Neanderthals), thereby doubling the number of Neanderthals for which genome sequences are available. Genetic similarity among late Neanderthals is well predicted by their geographical location, and comparison to the genome of an older Neanderthal from the Caucasus indicates that a population turnover is likely to have occurred, either in the Caucasus or throughout Europe, towards the end of Neanderthal history. We find that the bulk of Neanderthal gene flow into early modern humans originated from one or more source populations that diverged from the Neanderthals that were studied here at least 70,000 years ago, but after they split from a previously sequenced Neanderthal from Siberia around 150,000 years ago. Although four of the Neanderthals studied here post-date the putative arrival of early modern humans into Europe, we do not detect any recent gene flow from early modern humans in their ancestry.

Nuclear DNA sequences from the Middle Pleistocene Sima de los Huesos hominins.

A unique assemblage of 28 hominin individuals, found in Sima de los Huesos in the Sierra de Atapuerca in Spain, has recently been dated to approximately 430,000 years ago. An interesting question is how these Middle Pleistocene hominins were related to those who lived in the Late Pleistocene epoch, in particular to Neanderthals in western Eurasia and to Denisovans, a sister group of Neanderthals so far known only from southern Siberia. While the Sima de los Huesos hominins share some derived morphological features with Neanderthals, the mitochondrial genome retrieved from one individual from Sima de los Huesos is more closely related to the mitochondrial DNA of Denisovans than to that of Neanderthals. However, since the mitochondrial DNA does not reveal the full picture of relationships among populations, we have investigated DNA preservation in several individuals found at Sima de los Huesos. Here we recover nuclear DNA sequences from two specimens, which show that the Sima de los Huesos hominins were related to Neanderthals rather than to Denisovans, indicating that the population divergence between Neanderthals and Denisovans predates 430,000 years ago. A mitochondrial DNA recovered from one of the specimens shares the previously described relationship to Denisovan mitochondrial DNAs, suggesting, among other possibilities, that the mitochondrial DNA gene pool of Neanderthals turned over later in their history.

The genetic history of Ice Age Europe.

Modern humans arrived in Europe ~45,000 years ago, but little is known about their genetic composition before the start of farming ~8,500 years ago. Here we analyse genome-wide data from 51 Eurasians from ~45,000-7,000 years ago. Over this time, the proportion of Neanderthal DNA decreased from 3-6% to around 2%, consistent with natural selection against Neanderthal variants in modern humans. Whereas there is no evidence of the earliest modern humans in Europe contributing to the genetic composition of present-day Europeans, all individuals between ~37,000 and ~14,000 years ago descended from a single founder population which forms part of the ancestry of present-day Europeans. An ~35,000-year-old individual from northwest Europe represents an early branch of this founder population which was then displaced across a broad region, before reappearing in southwest Europe at the height of the last Ice Age ~19,000 years ago. During the major warming period after ~14,000 years ago, a genetic component related to present-day Near Easterners became widespread in Europe. These results document how population turnover and migration have been recurring themes of European prehistory.

An early modern human from Romania with a recent Neanderthal ancestor.

Neanderthals are thought to have disappeared in Europe approximately 39,000-41,000 years ago but they have contributed 1-3% of the DNA of present-day people in Eurasia. Here we analyse DNA from a 37,000-42,000-year-old modern human from Pestera cu Oase, Romania. Although the specimen contains small amounts of human DNA, we use an enrichment strategy to isolate sites that are informative about its relationship to Neanderthals and present-day humans. We find that on the order of 6-9% of the genome of the Oase individual is derived from Neanderthals, more than any other modern human sequenced to date. Three chromosomal segments of Neanderthal ancestry are over 50 centimorgans in size, indicating that this individual had a Neanderthal ancestor as recently as four to six generations back. However, the Oase individual does not share more alleles with later Europeans than with East Asians, suggesting that the Oase population did not contribute substantially to later humans in Europe.

A mitochondrial genome sequence of a hominin from Sima de los Huesos.

Excavations of a complex of caves in the Sierra de Atapuerca in northern Spain have unearthed hominin fossils that range in age from the early Pleistocene to the Holocene. One of these sites, the 'Sima de los Huesos' ('pit of bones'), has yielded the world's largest assemblage of Middle Pleistocene hominin fossils, consisting of at least 28 individuals dated to over 300,000 years ago. The skeletal remains share a number of morphological features with fossils classified as Homo heidelbergensis and also display distinct Neanderthal-derived traits. Here we determine an almost complete mitochondrial genome sequence of a hominin from Sima de los Huesos and show that it is closely related to the lineage leading to mitochondrial genomes of Denisovans, an eastern Eurasian sister group to Neanderthals. Our results pave the way for DNA research on hominins from the Middle Pleistocene.

Historical biogeography of the leopard (Panthera pardus) and its extinct Eurasian populations.

BACKGROUND: Resolving the historical biogeography of the leopard (Panthera pardus) is a complex issue, because patterns inferred from fossils and from molecular data lack congruence. Fossil evidence supports an African origin, and suggests that leopards were already present in Eurasia during the Early Pleistocene. Analysis of DNA sequences however, suggests a more recent, Middle Pleistocene shared ancestry of Asian and African leopards. These contrasting patterns led researchers to propose a two-stage hypothesis of leopard dispersal out of Africa: an initial Early Pleistocene colonisation of Asia and a subsequent replacement by a second colonisation wave during the Middle Pleistocene. The status of Late Pleistocene European leopards within this scenario is unclear: were these populations remnants of the first dispersal, or do the last surviving European leopards share more recent ancestry with their African counterparts? RESULTS: In this study, we generate and analyse mitogenome sequences from historical samples that span the entire modern leopard distribution, as well as from Late Pleistocene remains. We find a deep bifurcation between African and Eurasian mitochondrial lineages (~ 710 Ka), with the European ancient samples as sister to all Asian lineages (~ 483 Ka). The modern and historical mainland Asian lineages share a relatively recent common ancestor (~ 122 Ka), and we find one Javan sample nested within these. CONCLUSIONS: The phylogenetic placement of the ancient European leopard as sister group to Asian leopards suggests that these populations originate from the same out-of-Africa dispersal which founded the Asian lineages. The coalescence time found for the mitochondrial lineages aligns well with the earliest undisputed fossils in Eurasia, and thus encourages a re-evaluation of the identification of the much older putative leopard fossils from the region. The relatively recent ancestry of all mainland Asian leopard lineages suggests that these populations underwent a severe population bottleneck during the Pleistocene. Finally, although only based on a single sample, the unexpected phylogenetic placement of the Javan leopard could be interpreted as evidence for exchange of mitochondrial lineages between Java and mainland Asia, calling for further investigation into the evolutionary history of this subspecies.

Patterns of coding variation in the complete exomes of three Neandertals.

We present the DNA sequence of 17,367 protein-coding genes in two Neandertals from Spain and Croatia and analyze them together with the genome sequence recently determined from a Neandertal from southern Siberia. Comparisons with present-day humans from Africa, Europe, and Asia reveal that genetic diversity among Neandertals was remarkably low, and that they carried a higher proportion of amino acid-changing (nonsynonymous) alleles inferred to alter protein structure or function than present-day humans. Thus, Neandertals across Eurasia had a smaller long-term effective population than present-day humans. We also identify amino acid substitutions in Neandertals and present-day humans that may underlie phenotypic differences between the two groups. We find that genes involved in skeletal morphology have changed more in the lineage leading to Neandertals than in the ancestral lineage common to archaic and modern humans, whereas genes involved in behavior and pigmentation have changed more on the modern human lineage.

Long-Term Balancing Selection in LAD1 Maintains a Missense Trans-Species Polymorphism in Humans, Chimpanzees, and Bonobos.

Balancing selection maintains advantageous genetic and phenotypic diversity in populations. When selection acts for long evolutionary periods selected polymorphisms may survive species splits and segregate in present-day populations of different species. Here, we investigate the role of long-term balancing selection in the evolution of protein-coding sequences in the Homo-Pan clade. We sequenced the exome of 20 humans, 20 chimpanzees, and 20 bonobos and detected eight coding trans-species polymorphisms (trSNPs) that are shared among the three species and have segregated for approximately 14 My of independent evolution. Although the majority of these trSNPs were found in three genes of the major histocompatibility locus cluster, we also uncovered one coding trSNP (rs12088790) in the gene LAD1. All these trSNPs show clustering of sequences by allele rather than by species and also exhibit other signatures of long-term balancing selection, such as segregating at intermediate frequency and lying in a locus with high genetic diversity. Here, we focus on the trSNP in LAD1, a gene that encodes for Ladinin-1, a collagenous anchoring filament protein of basement membrane that is responsible for maintaining cohesion at the dermal-epidermal junction; the gene is also an autoantigen responsible for linear IgA disease. This trSNP results in a missense change (Leucine257Proline) and, besides altering the protein sequence, is associated with changes in gene expression of LAD1.

Complete mitochondrial genome sequence of a Middle Pleistocene cave bear reconstructed from ultrashort DNA fragments.

Although an inverse relationship is expected in ancient DNA samples between the number of surviving DNA fragments and their length, ancient DNA sequencing libraries are strikingly deficient in molecules shorter than 40 bp. We find that a loss of short molecules can occur during DNA extraction and present an improved silica-based extraction protocol that enables their efficient retrieval. In combination with single-stranded DNA library preparation, this method enabled us to reconstruct the mitochondrial genome sequence from a Middle Pleistocene cave bear (Ursus deningeri) bone excavated at Sima de los Huesos in the Sierra de Atapuerca, Spain. Phylogenetic reconstructions indicate that the U. deningeri sequence forms an early diverging sister lineage to all Western European Late Pleistocene cave bears. Our results prove that authentic ancient DNA can be preserved for hundreds of thousand years outside of permafrost. Moreover, the techniques presented enable the retrieval of phylogenetically informative sequences from samples in which virtually all DNA is diminished to fragments shorter than 50 bp.

Illuminating the base of the annelid tree using transcriptomics.

Annelida is one of three animal groups possessing segmentation and is central in considerations about the evolution of different character traits. It has even been proposed that the bilaterian ancestor resembled an annelid. However, a robust phylogeny of Annelida, especially with respect to the basal relationships, has been lacking. Our study based on transcriptomic data comprising 68,750-170,497 amino acid sites from 305 to 622 proteins resolves annelid relationships, including Chaetopteridae, Amphinomidae, Sipuncula, Oweniidae, and Magelonidae in the basal part of the tree. Myzostomida, which have been indicated to belong to the basal radiation as well, are now found deeply nested within Annelida as sister group to Errantia in most analyses. On the basis of our reconstruction of a robust annelid phylogeny, we show that the basal branching taxa include a huge variety of life styles such as tube dwelling and deposit feeding, endobenthic and burrowing, tubicolous and filter feeding, and errant and carnivorous forms. Ancestral character state reconstruction suggests that the ancestral annelid possessed a pair of either sensory or grooved palps, bicellular eyes, biramous parapodia bearing simple chaeta, and lacked nuchal organs. Because the oldest fossil of Annelida is reported for Sipuncula (520 Ma), we infer that the early diversification of annelids took place at least in the Lower Cambrian.

Fine-scale recombination patterns differ between chimpanzees and humans.

Recombination rates seem to vary extensively along the human genome. Pedigree analysis suggests that rates vary by an order of magnitude when measured at the megabase scale, and at a finer scale, sperm typing studies point to the existence of recombination hotspots. These are short regions (1-2 kb) in which recombination rates are 10-1,000 times higher than the background rate. Less is known about how recombination rates change over time. Here we determined to what degree recombination rates are conserved among closely related species by estimating recombination rates from 14 Mb of linkage disequilibrium data in central chimpanzee and human populations. The results suggest that recombination hotspots are not conserved between the two species and that recombination rates in larger (50 kb) genomic regions are only weakly conserved. Therefore, the recombination landscape has changed markedly between the two species.

Opsins in onychophora (velvet worms) suggest a single origin and subsequent diversification of visual pigments in arthropods.

Multiple visual pigments, prerequisites for color vision, are found in arthropods, but the evolutionary origin of their diversity remains obscure. In this study, we explore the opsin genes in five distantly related species of Onychophora, using deep transcriptome sequencing and screening approaches. Surprisingly, our data reveal the presence of only one opsin gene (onychopsin) in each onychophoran species, and our behavioral experiments indicate a maximum sensitivity of onychopsin to blue-green light. In our phylogenetic analyses, the onychopsins represent the sister group to the monophyletic clade of visual r-opsins of arthropods. These results concur with phylogenomic support for the sister-group status of the Onychophora and Arthropoda and provide evidence for monochromatic vision in velvet worms and in the last common ancestor of Onychophora and Arthropoda. We conclude that the diversification of visual pigments and color vision evolved in arthropods, along with the evolution of compound eyes-one of the most sophisticated visual systems known.

Annotation of primate miRNAs by high throughput sequencing of small RNA libraries.

BACKGROUND: In addition to genome sequencing, accurate functional annotation of genomes is required in order to carry out comparative and evolutionary analyses between species. Among primates, the human genome is the most extensively annotated. Human miRNA gene annotation is based on multiple lines of evidence including evidence for expression as well as prediction of the characteristic hairpin structure. In contrast, most miRNA genes in non-human primates are annotated based on homology without any expression evidence. We have sequenced small-RNA libraries from chimpanzee, gorilla, orangutan and rhesus macaque from multiple individuals and tissues. Using patterns of miRNA expression in conjunction with a model of miRNA biogenesis we used these high-throughput sequencing data to identify novel miRNAs in non-human primates. RESULTS: We predicted 47 new miRNAs in chimpanzee, 240 in gorilla, 55 in orangutan and 47 in rhesus macaque. The algorithm we used was able to predict 64% of the previously known miRNAs in chimpanzee, 94% in gorilla, 61% in orangutan and 71% in rhesus macaque. We therefore added evidence for expression in between one and five tissues to miRNAs that were previously annotated based only on homology to human miRNAs. We increased from 60 to 175 the number miRNAs that are located in orthologous regions in humans and the four non-human primate species studied here. CONCLUSIONS: In this study we provide expression evidence for homology-based annotated miRNAs and predict de novo miRNAs in four non-human primate species. We increased the number of annotated miRNA genes and provided evidence for their expression in four non-human primates. Similar approaches using different individuals and tissues would improve annotation in non-human primates and allow for further comparative studies in the future.

Deep sequencing of small RNAs confirms an annelid affinity of Myzostomida.

Myzostomida comprise a group of marine worms associated mainly with echinoderms since the Carboniferous. Due to their unusual morphology the phylogenetic position in relation to other Lophotrochozoa is discussed since their description. According to different morphological and molecular markers the Myzostomida are either close to Platyzoa or Annelida. Here we investigated small non-coding RNAs of Myzostoma cirriferum to infer the phylogenetic position of myzostomids. Based on transcriptomic data collected by Illumina Deep Sequencing we analyzed the microRNA (miRNA) families occurring in M. cirriferum. Phylogenetic analysis revealed the presence of 13 miRNA-families exclusively shared by Annelida (including Sipuncula) and Myzostomida, as such highly significantly supporting an annelid origin of myzostomids. Furthermore, using a mapping-approach and secondary structure models we predicted several miRNA-candidates unique for myzostomids.

Transcriptional neoteny in the human brain.

In development, timing is of the utmost importance, and the timing of developmental processes often changes as organisms evolve. In human evolution, developmental retardation, or neoteny, has been proposed as a possible mechanism that contributed to the rise of many human-specific features, including an increase in brain size and the emergence of human-specific cognitive traits. We analyzed mRNA expression in the prefrontal cortex of humans, chimpanzees, and rhesus macaques to determine whether human-specific neotenic changes are present at the gene expression level. We show that the brain transcriptome is dramatically remodeled during postnatal development and that developmental changes in the human brain are indeed delayed relative to other primates. This delay is not uniform across the human transcriptome but affects a specific subset of genes that play a potential role in neural development.

Quantifying and reducing cross-contamination in single- and multiplex hybridization capture of ancient DNA.

The use of hybridization capture has enabled a massive upscaling in sample sizes for ancient DNA studies, allowing the analysis of hundreds of skeletal remains or sediments in single studies. Nevertheless, demands in throughput continue to grow, and hybridization capture has become a limiting step in sample preparation due to the large consumption of reagents, consumables and time. Here, we explored the possibility of improving the economics of sample preparation via multiplex capture, that is, the hybridization capture of pools of double-indexed ancient DNA libraries. We demonstrate that this strategy is feasible, at least for small genomic targets such as mitochondrial DNA, if the annealing temperature is increased and PCR cycles are limited in post-capture amplification to avoid index swapping by jumping PCR, which manifests as cross-contamination in resulting sequence data. We also show that the reamplification of double-indexed libraries to PCR plateau before or after hybridization capture can sporadically lead to small, but detectable cross-contamination even if libraries are amplified in separate reactions. We provide protocols for both manual capture and automated capture in 384-well format that are compatible with single- and multiplex capture and effectively suppress cross-contamination and artefact formation. Last, we provide a simple computational method for quantifying cross-contamination due to index swapping in double-indexed libraries, which we recommend using for routine quality checks in studies that are sensitive to cross-contamination.