The peopling of Europe and the cautionary tale of Y chromosome lineage R-M269

Genetic genealogy reveals true Y haplogroup of House of Bourbon contradicting recent identification of the presumed remains of two French Kings

Genetic analysis strongly increases the opportunity to identify skeletal remains or other biological samples from historical figures. However, validation of this identification is essential and should be done by DNA typing of living relatives. Based on the similarity of a limited set of Y-STRs, a blood sample and a head were recently identified as those belonging respectively to King Louis XVI and his paternal ancestor King Henry IV. Here, we collected DNA samples from three living males of the House of Bourbon to validate the since then controversial identification of these remains. The three living relatives revealed the Bourbon's Y-chromosomal variant on a high phylogenetic resolution for several members of the lineage between Henry IV and Louis XVI. This 'true' Bourbon's variant is different from the published Y-STR profiles of the blood as well as of the head. The earlier identifications of these samples can therefore not be validated. Moreover, matrilineal genealogical data revealed that the published mtDNA sequence of the head was also different from the one of a series of relatives. This therefore leads to the conclusion that the analyzed samples were not from the French kings. Our study once again demonstrated that in order to realize an accurate genetic identification of historical remains DNA typing of living persons, who are paternally or maternally related with the presumed donor of the samples, is required.

Multiplex genotyping assays for fine-resolution subtyping of the major human Y-chromosome haplogroups E, G, I, J, and R in anthropological, genealogical, and forensic investigations

Inherited DNA polymorphisms located within the nonrecombing portion of the human Y chromosome provide a powerful means of tracking the patrilineal ancestry of male individuals. Recently, we introduced an efficient genotyping method for the detection of the basal Y-chromosome haplogroups A to T, as well as an additional method for the dissection of haplogroup O into its sublineages. To further extend the use of the Y chromosome as an evolutionary marker, we here introduce a set of genotyping assays for fine-resolution subtyping of haplogroups E, G, I, J, and R, which make up the bulk of Western Eurasian and African Y chromosomes. The marker selection includes a total of 107 carefully selected bi-allelic polymorphisms that were divided into eight hierarchically organized multiplex assays (two for haplogroup E, one for I, one for J, one for G, and three for R) based on the single-base primer extension (SNaPshot) technology. Not only does our method allow for enhanced Y-chromosome lineage discrimination, the more restricted geographic distribution of the subhaplogroups covered also enables more fine-scaled estimations of patrilineal bio-geographic origin. Supplementing our previous method for basal Y-haplogroup detection, the currently introduced assays are thus expected to be of major relevance for future DNA studies targeting male-specific ancestry for forensic, anthropological, and genealogical purposes.

A comparison of Y-chromosomal lineage dating using either resequencing or Y-SNP plus Y-STR genotyping

We have compared phylogenies and time estimates for Y-chromosomal lineages based on resequencing ∼9 Mb of DNA and applying the program GENETREE to similar analyses based on the more standard approach of genotyping 26 Y-SNPs plus 21 Y-STRs and applying the programs NETWORK and BATWING. We find that deep phylogenetic structure is not adequately reconstructed after Y-SNP plus Y-STR genotyping, and that times estimated using observed Y-STR mutation rates are several-fold too recent. In contrast, an evolutionary mutation rate gives times that are more similar to the resequencing data. In principle, systematic comparisons of this kind can in future studies be used to identify the combinations of Y-SNP and Y-STR markers, and time estimation methodologies, that correspond best to resequencing data.

Uniparental markers in Italy reveal a sex-biased genetic structure and different historical strata

Located in the center of the Mediterranean landscape and with an extensive coastal line, the territory of what is today Italy has played an important role in the history of human settlements and movements of Southern Europe and the Mediterranean Basin. Populated since Paleolithic times, the complexity of human movements during the Neolithic, the Metal Ages and the most recent history of the two last millennia (involving the overlapping of different cultural and demic strata) has shaped the pattern of the modern Italian genetic structure. With the aim of disentangling this pattern and understanding which processes more importantly shaped the distribution of diversity, we have analyzed the uniparentally-inherited markers in ∼900 individuals from an extensive sampling across the Italian peninsula, Sardinia and Sicily. Spatial PCAs and DAPCs revealed a sex-biased pattern indicating different demographic histories for males and females. Besides the genetic outlier position of Sardinians, a North West–South East Y-chromosome structure is found in continental Italy. Such structure is in agreement with recent archeological syntheses indicating two independent and parallel processes of Neolithisation. In addition, date estimates pinpoint the importance of the cultural and demographic events during the late Neolithic and Metal Ages. On the other hand, mitochondrial diversity is distributed more homogeneously in agreement with older population events that might be related to the presence of an Italian Refugium during the last glacial period in Europe.

Increasing phylogenetic resolution still informative for Y chromosomal studies on West-European populations

Many Y-chromosomal lineages which are defined in the latest phylogenetic tree of the human Y chromosome by the Y Chromosome Consortium (YCC) in 2008 are distributed in (Western) Europe due to the fact that a large number of phylogeographic studies focus on this area. Therefore, the question arises whether newly discovered polymorphisms on the Y chromosome will still be interesting to study Western Europeans on a population genetic level. To address this question, the West-European region of Flanders (Belgium) was selected as study area since more than 1000 Y chromosomes from this area have previously been genotyped at the highest resolution of the 2008 YCC-tree and coupled to in-depth genealogical data. Based on these data the temporal changes of the population genetic pattern over the last centuries within Flanders were studied and the effects of several past gene flow events were identified. In the present study a set of recently reported novel Y-SNPs were genotyped to further characterize all those Flemish Y chromosomes that belong to haplogroups G, R-M269 and T. Based on this extended Y-SNP set the discrimination power increased drastically as previous large (sub-)haplogroups are now subdivided in several non-marginal groups. Next, the previously observed population structure within Flanders appeared to be the result of different gradients of independent sub-haplogroups. Moreover, for the first time within Flanders a significant East-West gradient was observed in the frequency of two R-M269 lineages, and this gradient is still present when considering the current residence of the DNA donors. Our results thus suggest that an update of the Y-chromosomal tree based on new polymorphisms is still useful to increase the discrimination power based on Y-SNPs and to study population genetic patterns in more detail, even in an already well-studied region such as Western Europe.

Uniparental markers of contemporary Italian population reveals details on its pre-Roman heritage

Background

According to archaeological records and historical documentation, Italy has been a melting point for populations of different geographical and ethnic matrices. Although Italy has been a favorite subject for numerous population genetic studies, genetic patterns have never been analyzed comprehensively, including uniparental and autosomal markers throughout the country.

Methods/Principal Findings

A total of 583 individuals were sampled from across the Italian Peninsula, from ten distant (if homogeneous by language) ethnic communities — and from two linguistic isolates (Ladins, Grecani Salentini). All samples were first typed for the mitochondrial DNA (mtDNA) control region and selected coding region SNPs (mtSNPs). This data was pooled for analysis with 3,778 mtDNA control-region profiles collected from the literature. Secondly, a set of Y-chromosome SNPs and STRs were also analyzed in 479 individuals together with a panel of autosomal ancestry informative markers (AIMs) from 441 samples. The resulting genetic record reveals clines of genetic frequencies laid according to the latitude slant along continental Italy – probably generated by demographical events dating back to the Neolithic. The Ladins showed distinctive, if more recent structure. The Neolithic contribution was estimated for the Y-chromosome as 14.5% and for mtDNA as 10.5%. Y-chromosome data showed larger differentiation between North, Center and South than mtDNA. AIMs detected a minor sub-Saharan component; this is however higher than for other European non-Mediterranean populations. The same signal of sub-Saharan heritage was also evident in uniparental markers.

Conclusions/Significance

Italy shows patterns of molecular variation mirroring other European countries, although some heterogeneity exists based on different analysis and molecular markers. From North to South, Italy shows clinal patterns that were most likely modulated during Neolithic times.

Towards improvements in the estimation of the coalescent: implications for the most effective use of Y chromosome short tandem repeat mutation rates

Over the past two decades, many short tandem repeat (STR) microsatellite loci on the human Y chromosome have been identified together with mutation rate estimates for the individual loci. These have been used to estimate the coalescent age, or the time to the most recent common ancestor (TMRCA) expressed in generations, in conjunction with the average square difference measure (ASD), an unbiased point estimator of TMRCA based upon the average within-locus allele variance between haplotypes. The ASD estimator, in turn, depends on accurate mutation rate estimates to be able to produce good approximations of the coalescent age of a sample. Here, a comparison is made between three published sets of per locus mutation rate estimates as they are applied to the calculation of the coalescent age for real and simulated population samples. A novel evaluation method is developed for estimating the degree of conformity of any Y chromosome STR locus of interest to the strict stepwise mutation model and specific recommendations are made regarding the suitability of thirty-two commonly used Y-STR loci for the purpose of estimating the coalescent. The use of the geometric mean for averaging ASD and across loci is shown to improve the consistency of the resulting estimates, with decreased sensitivity to outliers and to the number of STR loci compared or the particular set of mutation rates selected.

A calibrated human Y-chromosomal phylogeny based on resequencing

We have identified variants present in high-coverage complete sequences of 36 diverse human Y chromosomes from Africa, Europe, South Asia, East Asia, and the Americas, representing eight major haplogroups. After restricting our analysis to 8.97 Mb of the unique male-specific Y sequence, we identified 6662 high-confidence variants, including single-nucleotide polymorphisms (SNPs), multi-nucleotide polymorphisms (MNPs), and indels. We constructed phylogenetic trees using these variants, or subsets of them, and recapitulated the known structure of the tree. Assuming a male mutation rate of 1 × 10⁻⁹ per base pair per year, the time depth of the tree (haplogroups A3-R) was ∼101,000–115,000 yr, and the lineages found outside Africa dated to 57,000–74,000 yr, both as expected. In addition, we dated a striking Paleolithic male lineage expansion to 41,000–52,000 yr ago and the node representing the major European Y lineage, R1b, to 4000–13,000 yr ago, supporting a Neolithic origin for these modern European Y chromosomes. In all, we provide a nearly 10-fold increase in the number of Y markers with phylogenetic information, and novel historical insights derived from placing them on a calibrated phylogenetic tree.

The genetic landscape of Equatorial Guinea and the origin and migration routes of the Y chromosome haplogroup R-V88

Human Y chromosomes belonging to the haplogroup R1b1-P25, although very common in Europe, are usually rare in Africa. However, recently published studies have reported high frequencies of this haplogroup in the central-western region of the African continent and proposed that this represents a 'back-to-Africa' migration during prehistoric times. To obtain a deeper insight into the history of these lineages, we characterised the paternal genetic background of a population in Equatorial Guinea, a Central-West African country located near the region in which the highest frequencies of the R1b1 haplogroup in Africa have been found to date. In our sample, the large majority (78.6%) of the sequences belong to subclades in haplogroup E, which are the most frequent in Bantu groups. However, the frequency of the R1b1 haplogroup in our sample (17.0%) was higher than that previously observed for the majority of the African continent. Of these R1b1 samples, nine are defined by the V88 marker, which was recently discovered in Africa. As high microsatellite variance was found inside this haplogroup in Central-West Africa and a decrease in this variance was observed towards Northeast Africa, our findings do not support the previously hypothesised movement of Chadic-speaking people from the North across the Sahara as the explanation for these R1b1 lineages in Central-West Africa. The present findings are also compatible with an origin of the V88-derived allele in the Central-West Africa, and its presence in North Africa may be better explained as the result of a migration from the south during the mid-Holocene.

In the name of the migrant father--analysis of surname origins identifies genetic admixture events undetectable from genealogical records

Patrilineal heritable surnames are widely used to select autochthonous participants for studies on small-scale population genetic patterns owing to the unique link between the surname and a genetic marker, the Y-chromosome (Y-chr). Today, the question arises as to whether the surname origin will be informative on top of in-depth genealogical pedigrees. Admixture events that happened in the period after giving heritable surnames but before the start of genealogical records may be informative about the additional value of the surname origin. In this context, an interesting historical event is the demic migration from French-speaking regions in Northern France to the depopulated and Dutch-speaking region Flanders at the end of the sixteenth century. Y-chr subhaplogroups of individuals with a French/Roman surname that could be associated with this migration event were compared with those of a group with autochthonous Flemish surnames. Although these groups could not be differentiated based on in-depth genealogical data, they were significantly genetically different from each other. Moreover, the observed genetic divergence was related to the differences in the distributions of main Y-subhaplogroups between contemporary populations from Northern France and Flanders. Therefore, these results indicate that the surname origin can be an important feature on top of in-depth genealogical results to select autochthonous participants for a regional population genetic study based on Y-chromosomes.

Discovery of Western European R1b1a2 Y chromosome variants in 1000 genomes project data: an online community approach

The authors have used an online community approach, and tools that were readily available via the Internet, to discover genealogically and therefore phylogenetically relevant Y-chromosome polymorphisms within core haplogroup R1b1a2-L11/S127 (rs9786076). Presented here is the analysis of 135 unrelated L11 derived samples from the 1000 Genomes Project. We were able to discover new variants and build a much more complex phylogenetic relationship for L11 sub-clades. Many of the variants were further validated using PCR amplification and Sanger sequencing. The identification of these new variants will help further the understanding of population history including patrilineal migrations in Western and Central Europe where R1b1a2 is the most frequent haplogroup. The fine-grained phylogenetic tree we present here will also help to refine historical genetic dating studies. Our findings demonstrate the power of citizen science for analysis of whole genome sequence data.

Mitochondrial DNA signals of late glacial recolonization of Europe from near eastern refugia

Human populations, along with those of many other species, are thought to have contracted into a number of refuge areas at the height of the last Ice Age. European populations are believed to be, to a large extent, the descendants of the inhabitants of these refugia, and some extant mtDNA lineages can be traced to refugia in Franco-Cantabria (haplogroups H1, H3, V, and U5b1), the Italian Peninsula (U5b3), and the East European Plain (U4 and U5a). Parts of the Near East, such as the Levant, were also continuously inhabited throughout the Last Glacial Maximum, but unlike western and eastern Europe, no archaeological or genetic evidence for Late Glacial expansions into Europe from the Near East has hitherto been discovered. Here we report, on the basis of an enlarged whole-genome mitochondrial database, that a substantial, perhaps predominant, signal from mitochondrial haplogroups J and T, previously thought to have spread primarily from the Near East into Europe with the Neolithic population, may in fact reflect dispersals during the Late Glacial period, ∼19-12 thousand years (ka) ago.

Temporal differentiation across a West-European Y-chromosomal cline: genealogy as a tool in human population genetics

The pattern of population genetic variation and allele frequencies within a species are unstable and are changing over time according to different evolutionary factors. For humans, it is possible to combine detailed patrilineal genealogical records with deep Y-chromosome (Y-chr) genotyping to disentangle signals of historical population genetic structures because of the exponential increase in genetic genealogical data. To test this approach, we studied the temporal pattern of the 'autochthonous' micro-geographical genetic structure in the region of Brabant in Belgium and the Netherlands (Northwest Europe). Genealogical data of 881 individuals from Northwest Europe were collected, from which 634 family trees showed a residence within Brabant for at least one generation. The Y-chr genetic variation of the 634 participants was investigated using 110 Y-SNPs and 38 Y-STRs and linked to particular locations within Brabant on specific time periods based on genealogical records. Significant temporal variation in the Y-chr distribution was detected through a north-south gradient in the frequencies distribution of sub-haplogroup R1b1b2a1 (R-U106), next to an opposite trend for R1b1b2a2g (R-U152). The gradient on R-U106 faded in time and even became totally invisible during the Industrial Revolution in the first half of the nineteenth century. Therefore, genealogical data for at least 200 years are required to study small-scale 'autochthonous' population structure in Western Europe.