Assessing the recovery of Y chromosome microsatellites with population genomic data using Papio and Theropithecus genomes.

Y chromosome markers can shed light on male-specific population dynamics but for many species no such markers have been discovered and are available yet, despite the potential for recovering Y-linked loci from available genome sequences. Here, we investigated how effective available bioinformatic tools are in recovering informative Y chromosome microsatellites from whole genome sequence data. In order to do so, we initially explored a large dataset of whole genome sequences comprising individuals at various coverages belonging to different species of baboons (genus: Papio) using Y chromosome references belonging to the same genus and more distantly related species (Macaca mulatta). We then further tested this approach by recovering Y-STRs from available Theropithecus gelada genomes using Papio and Macaca Y chromosome as reference sequences. Identified loci were validated in silico by a) comparing within-species relationships of Y chromosome lineages and b) genotyping male individuals in available pedigrees. Each STR was selected not to extend in its variable region beyond 100 base pairs, so that loci can be developed for PCR-based genotyping of non-invasive DNA samples. In addition to assembling a first set of Papio and Theropithecus Y-specific microsatellite markers, we released TYpeSTeR, an easy-to-use script to identify and genotype Y chromosome STRs using population genomic data which can be modulated according to available male reference genomes and genomic data, making it widely applicable across taxa.

Genomic variation in baboons from central Mozambique unveils complex evolutionary relationships with other Papio species.

BACKGROUND: Gorongosa National Park in Mozambique hosts a large population of baboons, numbering over 200 troops. Gorongosa baboons have been tentatively identified as part of Papio ursinus on the basis of previous limited morphological analysis and a handful of mitochondrial DNA sequences. However, a recent morphological and morphometric analysis of Gorongosa baboons pinpointed the occurrence of several traits intermediate between P. ursinus and P. cynocephalus, leaving open the possibility of past and/or ongoing gene flow in the baboon population of Gorongosa National Park. In order to investigate the evolutionary history of baboons in Gorongosa, we generated high and low coverage whole genome sequence data of Gorongosa baboons and compared it to available Papio genomes. RESULTS: We confirmed that P. ursinus is the species closest to Gorongosa baboons. However, the Gorongosa baboon genomes share more derived alleles with P. cynocephalus than P. ursinus does, but no recent gene flow between P. ursinus and P. cynocephalus was detected when available Papio genomes were analyzed. Our results, based on the analysis of autosomal, mitochondrial and Y chromosome data, suggest complex, possibly male-biased, gene flow between Gorongosa baboons and P. cynocephalus, hinting to direct or indirect contributions from baboons belonging to the "northern" Papio clade, and signal the presence of population structure within P. ursinus. CONCLUSIONS: The analysis of genome data generated from baboon samples collected in central Mozambique highlighted a complex set of evolutionary relationships with other baboons. Our results provided new insights in the population dynamics that have shaped baboon diversity.

Ancient DNA at the edge of the world: Continental immigration and the persistence of Neolithic male lineages in Bronze Age Orkney.

Orkney was a major cultural center during the Neolithic, 3800 to 2500 BC. Farming flourished, permanent stone settlements and chambered tombs were constructed, and long-range contacts were sustained. From ∼3200 BC, the number, density, and extravagance of settlements increased, and new ceremonial monuments and ceramic styles, possibly originating in Orkney, spread across Britain and Ireland. By ∼2800 BC, this phenomenon was waning, although Neolithic traditions persisted to at least 2500 BC. Unlike elsewhere in Britain, there is little material evidence to suggest a Beaker presence, suggesting that Orkney may have developed along an insular trajectory during the second millennium BC. We tested this by comparing new genomic evidence from 22 Bronze Age and 3 Iron Age burials in northwest Orkney with Neolithic burials from across the archipelago. We identified signals of inward migration on a scale unsuspected from the archaeological record: As elsewhere in Bronze Age Britain, much of the population displayed significant genome-wide ancestry deriving ultimately from the Pontic-Caspian Steppe. However, uniquely in northern and central Europe, most of the male lineages were inherited from the local Neolithic. This suggests that some male descendants of Neolithic Orkney may have remained distinct well into the Bronze Age, although there are signs that this had dwindled by the Iron Age. Furthermore, although the majority of mitochondrial DNA lineages evidently arrived afresh with the Bronze Age, we also find evidence for continuity in the female line of descent from Mesolithic Britain into the Bronze Age and even to the present day.

Mutations along human chromosomes: How randomly scattered are they?

The diversity of mutations in human chromosomes is nowadays very well documented. The mutations characterize populations in the world as well as genetic causes of diseases. In the approach that we follow, we study the patterns of gaps between mutations by means of the rescaled range analysis and the fractal dimension estimates. The results for chromosomes 1 to 22 and X indicate the existence of the so-called Hurst phenomenon in all of them. The interpretation of this outcome entails the presence of long-range correlations and we propose an explanation based on the genomic feature dubbed linkage disequilibrium, a nonrandom association of alleles at different loci. An unexpected outcome is the noteworthy uniform reduction in the Hurst phenomenon when considering the centimorgan metric instead of base position units. By contrast, such uniform reduction is not observed with the fractal dimension values.

Biomolecular insights into North African-related ancestry, mobility and diet in eleventh-century Al-Andalus.

Historical records document medieval immigration from North Africa to Iberia to create Islamic al-Andalus. Here, we present a low-coverage genome of an eleventh century CE man buried in an Islamic necropolis in Segorbe, near Valencia, Spain. Uniparental lineages indicate North African ancestry, but at the autosomal level he displays a mosaic of North African and European-like ancestries, distinct from any present-day population. Altogether, the genome-wide evidence, stable isotope results and the age of the burial indicate that his ancestry was ultimately a result of admixture between recently arrived Amazigh people (Berbers) and the population inhabiting the Peninsula prior to the Islamic conquest. We detect differences between our sample and a previously published group of contemporary individuals from Valencia, exemplifying how detailed, small-scale aDNA studies can illuminate fine-grained regional and temporal differences. His genome demonstrates how ancient DNA studies can capture portraits of past genetic variation that have been erased by later demographic shifts-in this case, most likely the seventeenth century CE expulsion of formerly Islamic communities as tolerance dissipated following the Reconquista by the Catholic kingdoms of the north.

A Geometrical Framework for f-Statistics.

A detailed derivation of the f-statistics formalism is made from a geometrical framework. It is shown that the f-statistics appear when a genetic distance matrix is constrained to describe a four population phylogenetic tree. The choice of genetic metric is crucial and plays an outstanding role as regards the tree-like-ness criterion. The case of lack of treeness is interpreted in the formalism as the presence of population admixture. In this respect, four formulas are given to estimate the admixture proportions. One of them is the so-called [Formula: see text]-ratio estimate and we show that a second one is related to a known result developed in terms of the fixation index [Formula: see text]. An illustrative numerical simulation of admixture proportion estimates is included. Relationships of the formalism with coalescence times and pairwise sequence differences are also provided.

Evolutionary distances corrected for purifying selection and ancestral polymorphisms.

Evolutionary distance formulas that take into account effects due to ancestral polymorphisms and purifying selection are obtained on the basis of the full solution of Jukes-Cantor and Kimura DNA substitution models. In the case of purifying selection two different methods are developed. It is shown that avoiding the dimensional reduction implicitly carried out in the conventional model solving is instrumental to incorporate the quoted effects into the formalism. The problem of estimating the numerical values of the model parameters, as well as those of the correction terms, is not addressed.