A western route of prehistoric human migration from Africa into the Iberian Peninsula.

Being at the western fringe of Europe, Iberia had a peculiar prehistory and a complex pattern of Neolithization. A few studies, all based on modern populations, reported the presence of DNA of likely African origin in this region, generally concluding it was the result of recent gene flow, probably during the Islamic period. Here, we provide evidence of much older gene flow from Africa to Iberia by sequencing whole genomes from four human remains from northern Portugal and southern Spain dated around 4000 years BP (from the Middle Neolithic to the Bronze Age). We found one of them to carry an unequivocal sub-Saharan mitogenome of most probably West or West-Central African origin, to our knowledge never reported before in prehistoric remains outside Africa. Our analyses of ancient nuclear genomes show small but significant levels of sub-Saharan African affinity in several ancient Iberian samples, which indicates that what we detected was not an occasional individual phenomenon, but an admixture event recognizable at the population level. We interpret this result as evidence of an early migration process from Africa into the Iberian Peninsula through a western route, possibly across the Strait of Gibraltar.

The genetics of an early Neolithic pastoralist from the Zagros, Iran.

The agricultural transition profoundly changed human societies. We sequenced and analysed the first genome (1.39x) of an early Neolithic woman from Ganj Dareh, in the Zagros Mountains of Iran, a site with early evidence for an economy based on goat herding, ca. 10,000 BP. We show that Western Iran was inhabited by a population genetically most similar to hunter-gatherers from the Caucasus, but distinct from the Neolithic Anatolian people who later brought food production into Europe. The inhabitants of Ganj Dareh made little direct genetic contribution to modern European populations, suggesting those of the Central Zagros were somewhat isolated from other populations of the Fertile Crescent. Runs of homozygosity are of a similar length to those from Neolithic farmers, and shorter than those of Caucasus and Western Hunter-Gatherers, suggesting that the inhabitants of Ganj Dareh did not undergo the large population bottleneck suffered by their northern neighbours. While some degree of cultural diffusion between Anatolia, Western Iran and other neighbouring regions is possible, the genetic dissimilarity between early Anatolian farmers and the inhabitants of Ganj Dareh supports a model in which Neolithic societies in these areas were distinct.

Ancient Ethiopian genome reveals extensive Eurasian admixture throughout the African continent.

Characterizing genetic diversity in Africa is a crucial step for most analyses reconstructing the evolutionary history of anatomically modern humans. However, historic migrations from Eurasia into Africa have affected many contemporary populations, confounding inferences. Here, we present a 12.5× coverage ancient genome of an Ethiopian male ("Mota") who lived approximately 4500 years ago. We use this genome to demonstrate that the Eurasian backflow into Africa came from a population closely related to Early Neolithic farmers, who had colonized Europe 4000 years earlier. The extent of this backflow was much greater than previously reported, reaching all the way to Central, West, and Southern Africa, affecting even populations such as Yoruba and Mbuti, previously thought to be relatively unadmixed, who harbor 6 to 7% Eurasian ancestry.

Paging through history: parchment as a reservoir of ancient DNA for next generation sequencing.

Parchment represents an invaluable cultural reservoir. Retrieving an additional layer of information from these abundant, dated livestock-skins via the use of ancient DNA (aDNA) sequencing has been mooted by a number of researchers. However, prior PCR-based work has indicated that this may be challenged by cross-individual and cross-species contamination, perhaps from the bulk parchment preparation process. Here we apply next generation sequencing to two parchments of seventeenth and eighteenth century northern English provenance. Following alignment to the published sheep, goat, cow and human genomes, it is clear that the only genome displaying substantial unique homology is sheep and this species identification is confirmed by collagen peptide mass spectrometry. Only 4% of sequence reads align preferentially to a different species indicating low contamination across species. Moreover, mitochondrial DNA sequences suggest an upper bound of contamination at 5%. Over 45% of reads aligned to the sheep genome, and even this limited sequencing exercise yield 9 and 7% of each sampled sheep genome post filtering, allowing the mapping of genetic affinity to modern British sheep breeds. We conclude that parchment represents an excellent substrate for genomic analyses of historical livestock.

A paleogenomic perspective on evolution and gene function: new insights from ancient DNA.

The publication of partial and complete paleogenomes within the last few years has reinvigorated research in ancient DNA. No longer limited to short fragments of mitochondrial DNA, inference of evolutionary processes through time can now be investigated from genome-wide data sampled as far back as 700,000 years. Tremendous insights have been made, in particular regarding the hominin lineage. With rare exception, however, a paleogenomic perspective has been mired by the quality and quantity of recoverable DNA. Though conceptually simple, extracting ancient DNA remains challenging, and sequencing ancient genomes to high coverage remains prohibitively expensive for most laboratories. Still, with improvements in DNA isolation and declining sequencing costs, the taxonomic and geographic purview of paleogenomics is expanding at a rapid pace. With improved capacity to screen large numbers of samples for those with high proportions of endogenous ancient DNA, paleogenomics is poised to become a key technology to better understand recent evolutionary events.

A phantom extinction? New insights into extinction dynamics of the Don-hare Lepus tanaiticus.

The Pleistocene to Holocene transition was accompanied by a worldwide extinction event affecting numerous mammalian species. Several species such as the woolly mammoth and the giant deer survived this extinction wave, only to go extinct a few thousand years later during the Holocene. Another example for such a Holocene extinction is the Don-hare, Lepus tanaiticus, which inhabited the Russian plains during the late glacial. After being slowly replaced by the extant mountain hare (Lepus timidus), it eventually went extinct during the middle Holocene. Here, we report the phylogenetic relationship of L. tanaiticus and L. timidus based on a 339-basepair (bp) fragment of the mitochondrial D-loop. Phylogenetic tree- and network reconstructions do not support L. tanaiticus and L. timidus being different species. Rather, we suggest that the two taxa represent different morphotypes of a single species and the extinction of 'L. tanaiticus' represents the disappearance of a local morphotype rather than the extinction of a species.

Patterns of nucleotide misincorporations during enzymatic amplification and direct large-scale sequencing of ancient DNA.

Whereas evolutionary inferences derived from present-day DNA sequences are by necessity indirect, ancient DNA sequences provide a direct view of past genetic variants. However, base lesions that accumulate in DNA over time may cause nucleotide misincorporations when ancient DNA sequences are replicated. By repeated amplifications of mitochondrial DNA sequences from a large number of ancient wolf remains, we show that C/G-to-T/A transitions are the predominant type of such misincorporations. Using a massively parallel sequencing method that allows large numbers of single DNA strands to be sequenced, we show that modifications of C, as well as to a lesser extent of G, residues cause such misincorporations. Experiments where oligonucleotides containing modified bases are used as templates in amplification reactions suggest that both of these types of misincorporations can be caused by deamination of the template bases. New DNA sequencing methods in conjunction with knowledge of misincorporation processes have now, in principle, opened the way for the determination of complete genomes from organisms that became extinct during and after the last glaciation.

Nuclear insertions help and hinder inference of the evolutionary history of gorilla mtDNA.

Numts are fragments of mitochondrial DNA (mtDNA) that have been translocated to the nucleus, where they can persist while their mitochondrial counterparts continue to rapidly evolve. Thus, numts represent 'molecular fossils' useful for comparison with mitochondrial variation, and are particularly suited for studies of the fast-evolving hypervariable segment of the mitochondrial control region (HV1). Here we used information from numts found in western gorillas (Gorilla gorilla) and eastern gorillas (Gorilla beringei) to estimate that these two species diverged about 1.3 million years ago (Ma), an estimate similar to recent calculations for the divergence of chimpanzee and bonobo. We also describe the sequence of a gorilla numt still possessing a segment lost from all contemporary gorilla mtDNAs. In contrast to that sequence, many numts of the HV1 are highly similar to authentic mitochondrial organellar sequences, making it difficult to determine whether purported mitochondrial sequences truly derive from that genome. We used all available organellar HV1 and corresponding numt sequences from gorillas in a phylogenetic analysis aimed at distinguishing these two types of sequences. Numts were found in several clades in the tree. This, in combination with the fact that only a limited amount of the extant variation in gorillas has been sampled, suggests that categorization of new sequences by the indirect means of phylogenetic comparison would be prone to uncertainty. We conclude that for taxa such as gorillas that contain numerous numts, direct approaches to the authentication of HV1 sequences, such as amplification strategies relying upon the circularity of the mtDNA molecule, remain necessary.

Mutations induced by ancient DNA extracts?

We have investigated whether some factor in ancient DNA extracts induces site-specific mutations in modern DNA. We find no evidence for higher mutation rates when extracts from three different Pleistocene mammals are added to modern DNA than when water or extraction blanks are added. We also fail to find evidence that any such factor affects ancient DNA sequences determined from the same extracts. This as well as the patterns of nucleotide substitutions seen in DNA sequences determined from hundreds of other specimens leads us to doubt that a previously unknown mutagenic factor can be a general feature of extracts from old tissues.

DNA sequences from multiple amplifications reveal artifacts induced by cytosine deamination in ancient DNA.

We show that DNA molecules amplified by PCR from DNA extracted from animal bones and teeth that vary in age between 25 000 and over 50 000 years carry C-->T and G-->A substitutions. These substitutions can reach high proportions among the molecules amplified and are due to the occurrence of modified deoxycytidine residues in the template DNA. If the template DNA is treated with uracil N-glycosylase, these substitutions are dramatically reduced. They are thus likely to result from deamination of deoxycytidine residues. In addition, 'jumping PCR', i.e. the occurrence of template switching during PCR, may contribute to these substitutions. When DNA sequences are amplified from ancient DNA extracts where few template molecules initiate the PCR, precautions such as DNA sequence determination of multiple clones derived from more than one independent amplification are necessary in order to reduce the risk of determination of incorrect DNA sequences. When such precautionary measures are taken, errors induced by damage to the DNA template are unlikely to be more frequent than approximately 0.1% even under the unlikely scenario where each amplification starts from a single template molecule.

Paternity and relatedness in wild chimpanzee communities.

The genetic structure of three contiguous wild chimpanzee communities in West Africa was examined to determine the extent to which the community, the mixed-sex social unit of chimpanzees, represents a closed reproductive unit. An analysis of paternity for 41 offspring resulted in 34 cases of paternity assignment to an adult male belonging to the same community. Among the 14 offspring for which all potential within-community fathers have been tested, one likely case of extra-group paternity (EGP) has been identified, suggesting an incidence of EGP of 7%. This more extensive analysis contradicts a previous genetic study of the Tai chimpanzees that inferred 50% extra-group fathers. We suggest, based on direct comparison of results for 33 individuals at 1 microsatellite locus and direct comparison of paternity assignments for 11 offspring, that the error rate in the previous study was too high to produce accurate genotypes and assignments of paternity and hence caused the false inference of a high rate of EGP. Thus, the community is the primary but not exclusive unit for reproduction in wild chimpanzees, and females do not typically reproduce with outside males. Despite the inferred low level of gene flow from extra-community males, relatedness levels among the community males are not significantly higher than among community females, and the distribution of genetic relationships within the group suggests that, rather than a primarily male-bonded social structure, the group is bonded through relationships between males and females. Kinship may explain cooperative behaviors directed against other communities, but is unlikely to explain the high levels of affiliation and cooperation seen for male within-community interactions.

Exclusion of linkage to the CDL1 gene region on chromosome 3q26.3 in some familial cases of Cornelia de Lange syndrome.

Cornelia de Lange Syndrome (CdLS) is a complex developmental disorder consisting of characteristic facial features, limb abnormalities, hirsutism, ophthalmologic involvement, gastroesophageal dysfunction, hearing loss, as well as growth and neurodevelopmental retardation. Most cases of CdLS appear to be sporadic. Familial cases are rare and indicate autosomal dominant inheritance. Several individuals with CdLS have been reported with chromosomal abnormalities, suggesting candidate genomic regions within which the causative gene(s) may lie. A CdLS gene location (CDL1) has been assigned to 3q26.3 based on phenotypic overlap with the duplication 3q syndrome (critical region 3q26.2-q27) and the report of a CdLS individual with a balanced de novo t(3;17)(q26.3;q23.1). It has been postulated that a gene within the dup3q critical region results in the CdLS when deleted or mutated. We have performed a linkage analysis to the minimal critical region for the dup3q syndrome (that encompasses the translocation breakpoint) on chromosome 3q in 10 rare familial cases of CdLS. Nineteen markers spanning a region of approximately 40 Mb (37 cM) were used. Results of a multipoint linkage analysis demonstrated total lod-scores that were negative across the chromosome 3q26-q27 region. In 4/10 families, lod-scores were less than -2 in the 2 cM region encompassing the translocation, while in the remaining 6/10 families, lod-scores could not exclude linkage to this region. These studies indicate that in some multicase families, the disease gene does not map to the CDL1 region at 3q26.3.

Ancient DNA.

DNA that has been recovered from archaeological and palaeontological remains makes it possible to go back in time and study the genetic relationships of extinct organisms to their contemporary relatives. This provides a new perspective on the evolution of organisms and DNA sequences. However, the field is fraught with technical pitfalls and needs stringent criteria to ensure the reliability of results, particularly when human remains are studied.

A molecular analysis of ground sloth diet through the last glaciation.

DNA was extracted from five coprolites, excavated in Gypsum Cave, Nevada and radiocarbon dated to approximately 11 000, 20 000 and 28 500 years BP. All coprolites contained mitochondrial DNA sequences identical to a DNA sequence determined from a bone of the extinct ground sloth Nothrotheriops shastensis. A 157-bp fragment of the chloroplast gene for the large subunit of the ribulosebisphosphate carboxylase (rbcL) was amplified from the boluses and several hundred clones were sequenced. In addition, the same DNA fragment was sequenced from 99 plant species that occur in the vicinity of Gypsum Cave today. When these were compared to the DNA sequences in GenBank, 69 were correctly (two incorrectly) assigned to taxonomic orders. The plant sequences from the five coprolites as well as from one previously studied coprolite were compared to rbcL sequences in GenBank and the contemporary plant species. Thirteen families or orders of plants that formed part of the diet of the Shasta ground sloth could be identified, showing that the ground sloth was feeding on trees as well as herbs and grasses. The plants in the boluses further indicate that the climate 11 000 years BP was dryer than 20 000 and 28 500 years BP. However, the sloths seem to have visited water sources more frequently at 11 000 BP than at earlier times.

Modulation of volume regulated anion current by I(Cln).

I(Cln), a cytosolic protein associated with a nucleotide-sensitive chloride current, may be involved in the regulation of a volume-regulated anion current (VRAC) associated with hypotonic cell swelling. We have determined the nucleic acid sequences of I(Cln) from human tsA201a, colonic (T84) and myeloma (RPMI 8826) cell lines. The amino acid sequences are highly homologous (>/=99%) to each other but less homologous to I(Cln) protein from other species. Using the whole-cell patch clamp technique, we examined the effect of I(Cln) protein expression levels on VRAC properties during a hyposmotic challenge. Overexpression of T84 or RPMI 8226-derived I(Cln) protein in tsA201a cells results in a more than 9-fold increase in the rate of VRAC activation over control values, while having no effect on VRAC inactivation properties. Underexpression of endogenous I(Cln) protein in tsA201a cells using antisense oligonucleotides results in a more than 180-fold decrease in VRAC activation rate as compared to control values. These results indicate that I(Cln) protein expression modulates VRAC activation but not inactivation.

Molecular coproscopy: dung and diet of the extinct ground sloth Nothrotheriops shastensis.

DNA from excrements can be amplified by means of the polymerase chain reaction. However, this has not been possible with ancient feces. Cross-links between reducing sugars and amino groups were shown to exist in a Pleistocene coprolite from Gypsum Cave, Nevada. A chemical agent, N-phenacylthiazolium bromide, that cleaves such cross-links made it possible to amplify DNA sequences. Analyses of these DNA sequences showed that the coprolite is derived from an extinct sloth, presumably the Shasta ground sloth Nothrotheriops shastensis. Plant DNA sequences from seven groups of plants were identified in the coprolite. The plant assemblage that formed part of the sloth's diet exists today at elevations about 800 meters higher than the cave.