A major Y-chromosome haplogroup R1b Holocene era founder effect in Central and Western Europe

Massive migration from the steppe was a source for Indo-European languages in Europe

We generated genome-wide data from 69 Europeans who lived between 8,000-3,000 years ago by enriching ancient DNA libraries for a target set of almost 400,000 polymorphisms. Enrichment of these positions decreases the sequencing required for genome-wide ancient DNA analysis by a median of around 250-fold, allowing us to study an order of magnitude more individuals than previous studies and to obtain new insights about the past. We show that the populations of Western and Far Eastern Europe followed opposite trajectories between 8,000-5,000 years ago. At the beginning of the Neolithic period in Europe, ∼8,000-7,000 years ago, closely related groups of early farmers appeared in Germany, Hungary and Spain, different from indigenous hunter-gatherers, whereas Russia was inhabited by a distinctive population of hunter-gatherers with high affinity to a ∼24,000-year-old Siberian. By ∼6,000-5,000 years ago, farmers throughout much of Europe had more hunter-gatherer ancestry than their predecessors, but in Russia, the Yamnaya steppe herders of this time were descended not only from the preceding eastern European hunter-gatherers, but also from a population of Near Eastern ancestry. Western and Eastern Europe came into contact ∼4,500 years ago, as the Late Neolithic Corded Ware people from Germany traced ∼75% of their ancestry to the Yamnaya, documenting a massive migration into the heartland of Europe from its eastern periphery. This steppe ancestry persisted in all sampled central Europeans until at least ∼3,000 years ago, and is ubiquitous in present-day Europeans. These results provide support for a steppe origin of at least some of the Indo-European languages of Europe.

Ladakh, India: the land of high passes and genetic heterogeneity reveals a confluence of migrations

Owing to its geographic location near the longitudinal center of Asia, Ladakh, the land of high passes, has witnessed numerous demographic movements during the past millenniums of occupation. In an effort to view Ladakh's multicultural history from a paternal genetic perspective, we performed a high-resolution Y-chromosomal survey of Ladakh, within the context of Y haplogroup and haplotype distributions of 41 Asian reference populations. The results of this investigation highlight the rich ethnic and genetic diversity of Ladkah which includes genetic contributions from disparate regions of the continent including, West, East, South and Central Asia. The phylogenetic signals from Ladakh are consistent with the Indo-Aryans' occupation during the Neolithic age and its historic connection with Tibet, as well as the East-West gene flow associated with the Silk Road.

Surnames and Y-chromosomal markers reveal low relationships in Southern Spain

A sample of 416 males from western and eastern Andalusia has been jointly analyzed for surnames and Y-chromosome haplogroups and haplotypes. The observed number of different surnames was 222 (353 when the second surname of the Spanish system of naming is considered). The great majority of recorded surnames have a Castilian-Leonese origin, while Catalan or Basque surnames have not been found. A few Arab-related surnames appear but none discernible of Sephardic-Jewish descent. Low correlation among surnames with different population frequencies and Y-chromosome markers, at different levels of genetic resolution, has been observed in Andalusia. This finding could be explained mainly by the very low rate of monophyletic surnames because of the historical process of surname ascription and the resulting high frequencies of the most common Spanish surnames. The introduction of surnames in Spain during the Middle Ages coincided with Reconquest of the territories under Islamic rule, and Muslims and Jews progressively adopted the present male line surname system. Sampled surnames and Y-chromosome lineages fit well a power-law distribution and observed isonymy is very close to that of the general population. Besides, our data and results show that the reliability of the isonymy method should be questioned because of the high rate of polyphyletic surnames, even in small geographic regions and autochthonous populations. Random isonymy would be consistently dependent of the most common surname frequencies in the population.

Sousse: extreme genetic heterogeneity in North Africa

The male genetic landscape of the territory currently known as Tunisia is hampered by the scarcity of data, especially from cosmopolitan areas such as the coastal city of Sousse. In order to alleviate this lacuna, 220 males from Sousse were examined, for the first time, for more than 50 Y-chromosome single-nucleotide polymorphisms (Y-SNPs) markers and compared with 3099 individuals from key geographically targeted locations in North Africa, Europe and the Near East. The paternal lineages observed belong to a common set of Y haplogroups previously described in North Africa. In addition to the prominent autochthonous North African E-M81 haplogroup which is exclusively represented by its subclade E-M183 (44.55% of Y-chromosomes), a number of Near Eastern Neolithic lineages including E-M78, J-M267 and J-M172 account for 39% of the Y-chromosomes detected. Principal component analysis based on haplogroup frequencies, multidimensional scaling based on Rst genetic distances and analyses of molecular variance using both Y-chromosome short tandem repeat haplotypes and Y-SNP haplogroup data revealed that the Tunisian and North African groups, as a whole, are intra- and inter-specific diverse with Sousse being highly heterogeneous.

The role of nutrition and genetics as key determinants of the positive height trend

The aim of this study was to identify the most important variables determining current differences in physical stature in Europe and some of its overseas offshoots such as Australia, New Zealand and USA. We collected data on the height of young men from 45 countries and compared them with long-term averages of food consumption from the FAOSTAT database, various development indicators compiled by the World Bank and the CIA World Factbook, and frequencies of several genetic markers. Our analysis demonstrates that the most important factor explaining current differences in stature among nations of European origin is the level of nutrition, especially the ratio between the intake of high-quality proteins from milk products, pork meat and fish, and low-quality proteins from wheat. Possible genetic factors such as the distribution of Y haplogroup I-M170, combined frequencies of Y haplogroups I-M170 and R1b-U106, or the phenotypic distribution of lactose tolerance emerge as comparably important, but the available data are more limited. Moderately significant positive correlations were also found with GDP per capita, health expenditure and partly with the level of urbanization that influences male stature in Western Europe. In contrast, male height correlated inversely with children's mortality and social inequality (Gini index). These results could inspire social and nutritional guidelines that would lead to the optimization of physical growth in children and maximization of the genetic potential, both at the individual and national level.

Different waves and directions of Neolithic migrations in the Armenian Highland

Background

The peopling of Europe and the nature of the Neolithic agricultural migration as a primary issue in the modern human colonization of the globe is still widely debated. At present, much uncertainty is associated with the reconstruction of the routes of migration for the first farmers from the Near East. In this context, hospitable climatic conditions and the key geographic position of the Armenian Highland suggest that it may have served as a conduit for several waves of expansion of the first agriculturalists from the Near East to Europe and the North Caucasus.

Results

Here, we assess Y-chromosomal distribution in six geographically distinct populations of Armenians that roughly represent the extent of historical Armenia. Using the general haplogroup structure and the specific lineages representing putative genetic markers of the Neolithic Revolution, haplogroups R1b1a2, J2, and G, we identify distinct patterns of genetic affinity between the populations of the Armenian Highland and the neighboring ones north and west from this area.

Conclusions

Based on the results obtained, we suggest a new insight on the different routes and waves of Neolithic expansion of the first farmers through the Armenian Highland. We detected at least two principle migratory directions: (1) westward alongside the coastline of the Mediterranean Sea and (2) northward to the North Caucasus.

Electronic supplementary material

The online version of this article (doi:10.1186/s13323-014-0015-6) contains supplementary material, which is available to authorized users.

Simultaneous analysis of hundreds of Y-chromosomal SNPs for high-resolution paternal lineage classification using targeted semiconductor sequencing

SNPs from the non-recombining part of the human Y chromosome (Y-SNPs) are informative to classify paternal lineages in forensic, genealogical, anthropological, and evolutionary studies. Although thousands of Y-SNPs were identified thus far, previous Y-SNP multiplex tools target only dozens of markers simultaneously, thereby restricting the provided Y-haplogroup resolution and limiting their applications. Here, we overcome this shortcoming by introducing a high-resolution multiplex tool for parallel genotyping-by-sequencing of 530 Y-SNPs using the Ion Torrent PGM platform, which allows classification of 432 worldwide Y haplogroups. Contrary to previous Y-SNP multiplex tools, our approach covers branches of the entire Y tree, thereby maximizing the paternal lineage classification obtainable. We used a default DNA input amount of 10 ng per reaction but preliminary sensitivity testing revealed positive results from as little as 100 pg input DNA. Furthermore, we demonstrate that sample pooling using barcodes is feasible, allowing increased throughput for lower per-sample costs. In addition to the wetlab protocol, we provide a software tool for automated data quality control and haplogroup classification. The unique combination of ultra-high marker density and high sensitivity achievable from low amounts of potentially degraded DNA makes this new multiplex tool suitable for a wide range of Y-chromosome applications.

Testing Central and Inner Asian admixture among contemporary Hungarians

Historically, the Carpathian Basin was the final destination for many nomadic peoples who migrated westward from Inner and Central Asia towards Europe. Proto-Hungarians (Steppe Magyars) were among those who came from the East, the Eurasian Steppe in the early middle ages. In order to detect the paternal genetic contribution from nomadic Steppe tribes, we tested 966 samples from Central Asian (Uzbekistan, Kazakhstan), Inner Asian (Mongolians and Buryats in Mongolia) and Hungarian-speaking European (Hungarian, Sekler and Csango) populations. We constructed median-joining networks of certain haplogroups in Hungarian-speaking European, and Altaic-speaking Central and Inner Asian populations. We estimated that the possible paternal genetic contribution from the above described populations among contemporary Hungarian speaking populations ranged between 5% and 7.4%. It is lowest among Hungarians from Hungary (5.1%), while higher among Hungarian-speaking groups in Romania, notably Sekler (7.4%) and Csango (6.3%). However, these results represent only an upper limit. Actual Central/Inner Asian admixture might be somewhat lower as some of the related lineages may have come from a common third source. The main haplogroups responsible for the Central/Inner Asian admixture among Hungarians are J2*-M172 (xM47, M67, M12), J2-L24, R1a-Z93; Q-M242 and E-M78. Earlier studies showed very limited Uralic genetic influence among Hungarians, and based on the present study, Altaic/Turkic genetic contribution is also not significant, although significantly higher than the Uralic one. The conclusion of this study is that present-day Hungarian speakers are genetically very similar to neighbouring populations, isolated Hungarian speaking groups having relatively higher presence of Central and Inner Asian genetic elements. At the same time, the reliable historical and genetic conclusions require an extension of the study to a significantly larger database with deep haplogroup resolution, including ancient DNA data.

New clustering methods for population comparison on paternal lineages

The goal of this study is to show two new clustering and visualising techniques developed to find the most typical clusters of 18-dimensional Y chromosomal haplogroup frequency distributions of 90 Western Eurasian populations. The first technique called "self-organizing cloud (SOC)" is a vector-based self-learning method derived from the Self Organising Map and non-metric Multidimensional Scaling algorithms. The second technique is a new probabilistic method called the "maximal relation probability" (MRP) algorithm, based on a probability function having its local maximal values just in the condensation centres of the input data. This function is calculated immediately from the distance matrix of the data and can be interpreted as the probability that a given element of the database has a real genetic relation with at least one of the remaining elements. We tested these two new methods by comparing their results to both each other and the k-medoids algorithm. By means of these new algorithms, we determined 10 clusters of populations based on the similarity of haplogroup composition. The results obtained represented a genetically, geographically and historically well-interpretable picture of 10 genetic clusters of populations mirroring the early spread of populations from the Fertile Crescent to the Caucasus, Central Asia, Arabia and Southeast Europe. The results show that a parallel clustering of populations using SOC and MRP methods can be an efficient tool for studying the demographic history of populations sharing common genetic footprints.

Human paleogenetics of Europe--the known knowns and the known unknowns

The number of ancient human DNA studies has drastically increased in recent years. This results in a substantial record of mitochondrial sequences available from many prehistoric sites across Western Eurasia, but also growing Y-chromosome and autosomal sequence data. We review the current state of research with specific emphasis on the Holocene population events that likely have shaped the present-day genetic variation in Europe. We reconcile observations from the genetic data with hypotheses about the peopling and settlement history from anthropology and archaeology for various key regions, and also discuss the data in light of evidence from related disciplines, such as modern human genetics, climatology and linguistics.

The place of Slovakian paternal diversity in the clinal European landscape

BACKGROUND

Several demographic events have been postulated to explain the contemporaneous structure of European genetic diversity. First, an initial settlement of the continent by anatomically modern humans; second, the re-settlement of northern latitudes after the Last Glacial Maximum; third, the demic diffusion of Neolithic farmers from the Near East; and, fourth, several historical events such as the Slavic migration.

AIM

The aim of this study was to provide a more integrated picture of male-specific genetic relationships of Slovakia within the broader pan-European genetic landscape.

SUBJECTS AND METHODS

This study analysed a new Y-chromosome data-set (156 individuals) for both SNP and STR polymorphisms in population samples from five different Slovakian localities.

RESULTS

It was found that the male diversity of Slovakia is embedded in the clinal pattern of the major R1a and R1b clades extending over the continent and a similar pattern of population structure is found with Y-specific SNP or STR variation.

CONCLUSION

The highly significant correlation between the results based on fast evolving STRs on one hand and slow evolving SNPs on the other hand suggests a recent timeframe for the settlement of the area.

At the southeast fringe of the Bantu expansion: genetic diversity and phylogenetic relationships to other sub-Saharan tribes

Here, we present 12 loci paternal haplotypes (Y-STR profiles) against the backdrop of the Y-SNP marker system of Bantu males from the Maputo Province of Southeast Africa, a region believed to represent the southeastern fringe of the Bantu expansion. Our Maputo Bantu group was analyzed within the context of 27 geographically relevant reference populations in order to ascertain its genetic relationship to other Bantu and non Bantu (Pygmy, Khoisan and Nilotic) sub-equatorial tribes from West and East Africa. This study entails statistical pair wise comparisons and multidimensional scaling based on YSTR Rst distances, network analyses of Bantu (B2a-M150) and Pygmy (B2b-M112) lineages as well as an assessment of Y-SNP distribution patterns. Several notable findings include the following: 1) the Maputo Province Bantu exhibits a relatively close paternal affinity with both east and west Bantu tribes due to high proportion of Bantu Y chromosomal markers, 2) only traces of Khoisan (1.3%) and Pygmy (1.3%) markers persist in the Maputo Province Bantu gene pool, 3) the occurrence of R1a1a-M17/M198, a member of the Eurasian R1a-M420 branch in the population of the Maputo Province, may represent back migration events and/or recent admixture events, 4) the shared presence of E1b1b1-M35 in all Tanzanian tribes examined, including Bantu and non-Bantu groups, in conjunction with its nearly complete absence in the West African populations indicate that, in addition to a shared linguistic, cultural and genetic heritage, geography (e.g., east vs. west) may have impacted the paternal landscape of sub-Saharan Africa, 5) the admixture and assimilation processes of Bantu elements were both highly complex and region-specific.

•

Maputo Bantus exhibit close affinities with other West and East African Bantus.

•

Traces of Khoisan and Pygmy markers persist in the Maputo Province Bantus.

•

R1a1a-M17/M198 in the Maputo Province may represent back or recent migration.

•

Linguistic, cultural and genetic heritages are reflected in Maputo's gene pool.

•

Admixture and assimilation processes of Bantu elements were region-specific.

Population distribution and ancestry of the cancer protective MDM2 SNP285 (rs117039649)

The MDM2 promoter SNP285C is located on the SNP309G allele. While SNP309G enhances Sp1 transcription factor binding and MDM2 transcription, SNP285C antagonizes Sp1 binding and reduces the risk of breast-, ovary- and endometrial cancer. Assessing SNP285 and 309 genotypes across 25 different ethnic populations (>10.000 individuals), the incidence of SNP285C was 6-8% across European populations except for Finns (1.2%) and Saami (0.3%). The incidence decreased towards the Middle-East and Eastern Russia, and SNP285C was absent among Han Chinese, Mongolians and African Americans. Interhaplotype variation analyses estimated SNP285C to have originated about 14,700 years ago (95% CI: 8,300 – 33,300). Both this estimate and the geographical distribution suggest SNP285C to have arisen after the separation between Caucasians and modern day East Asians (17,000 - 40,000 years ago). We observed a strong inverse correlation (r = -0.805; p < 0.001) between the percentage of SNP309G alleles harboring SNP285C and the MAF for SNP309G itself across different populations suggesting selection and environmental adaptation with respect to MDM2 expression in recent human evolution. In conclusion, we found SNP285C to be a pan-Caucasian variant. Ethnic variation regarding distribution of SNP285C needs to be taken into account when assessing the impact of MDM2 SNPs on cancer risk.

Y-SNP L1034: limited genetic link between Mansi and Hungarian-speaking populations

Genetic studies noted that the Hungarian Y-chromosomal gene pool significantly differs from other Uralic-speaking populations. Hungarians show very limited or no presence of haplogroup N-Tat, which is frequent among other Uralic-speaking populations. We proposed that some genetic links need to be observed between the linguistically related Hungarian and Mansi populations.This is the first attempt to divide haplogroup N-Tat into subhaplogroups by testing new downstream SNP markers L708 and L1034. Sixty Northern Mansi samples were collected in Western Siberia and genotyped for Y-chromosomal haplotypes and haplogroups. We found 14 Mansi and 92 N-Tat samples from 7 populations. Comparative results showed that all N-Tat samples carried the N-L708 mutation. Some Hungarian, Sekler, and Uzbek samples were L1034 SNP positive, while all Mongolians, Buryats, Khanty, Finnish, and Roma samples yielded a negative result for this marker. Based on the above, L1034 marker seems to be a subgroup of N-Tat, which is typical for Mansi and Hungarian-speaking ethnic groups so far. Based on our time to most recent common ancestor data, the L1034 marker arose 2,500 years before present. The overall frequency of the L1034 is very low among the analyzed populations, thus it does not necessarily mean that proto-Hungarians and Mansi descend from common ancestors. It does provide, however, a limited genetic link supporting language contact. Both Hungarians and Mansi have much more complex genetic population history than the traditional tree-based linguistic model would suggest.

Traces of medieval migrations in a socially stratified population from Northern Italy. Evidence from uniparental markers and deep-rooted pedigrees

Social and cultural factors had a critical role in determining the genetic structure of Europe. Therefore, socially stratified populations may help to focus on specific episodes of European demographic history. In this study, we use uniparental markers to analyse the genetic structure of Partecipanza in San Giovanni in Persiceto (Northern Italy), a peculiar institution whose origins date back to the Middle Ages and whose members form the patrilineal descent of a group of founder families. From a maternal point of view (mtDNA), Partecipanza is genetically homogeneous with the rest of the population. However, we observed a significant differentiation for Y-chromosomes. In addition, by comparing 17 Y-STR profiles with deep-rooted paternal pedigrees, we estimated a Y-STR mutation rate equal to 3.90 * 10(-3) mutations per STR per generation and an average generation duration time of 33.38 years. When we used these values for tentative dating, we estimated 1300-600 years ago for the origins of the Partecipanza. These results, together with a peculiar Y-chromosomal composition and historical evidence, suggest that Germanic populations (Lombards in particular) settled in the area during the Migration Period (400-800 AD, approximately) and may have had an important role in the foundation of this community.

Male lineage strata of Brazilian population disclosed by the simultaneous analysis of STRs and SNPs

Brazil has a large territory divided in five geographical regions harboring highly diverse populations that resulted from different degrees and modes of admixture between Native Americans, Europeans and Africans. In this study, a sample of 605 unrelated males was genotyped for 17 Y-STRs and 46 Y-SNPs aiming a deep characterization of the male gene pool of Rio de Janeiro and its comparison with other Brazilian populations. High values of Y-STR haplotype diversity (0.9999±0.0001) and Y-SNP haplogroup diversity (0.7589±0.0171) were observed. Population comparisons at both haplotype and haplogroup levels showed significant differences between Brazilian South Eastern and Northern populations that can be explained by differences in the proportion of African and Native American Y chromosomes. Statistical significant differences between admixed urban samples from the five regions of Brazil were not previously detected at haplotype level based on smaller size samples from South East, which emphasizes the importance of sample size to detected population stratification for an accurate interpretation of profile matches in kinship and forensic casework. Although not having an intra-population discrimination power as high as the Y-STRs, the Y-SNPs are more powerful to disclose differences in admixed populations. In this study, the combined analysis of these two types of markers proved to be a good strategy to predict population sub-structure, which should be taken into account when delineating forensic database strategies for Y chromosome haplotypes.

Genetic diversity in Puerto Rico and its implications for the peopling of the Island and the West Indies

Puerto Rico and the surrounding islands rest on the eastern fringe of the Caribbean's Greater Antilles, located less than 100 miles northwest of the Lesser Antilles. Puerto Ricans are genetic descendants of pre-Columbian peoples, as well as peoples of European and African descent through 500 years of migration to the island. To infer these patterns of pre-Columbian and historic peopling of the Caribbean, we characterized genetic diversity in 326 individuals from the southeastern region of Puerto Rico and the island municipality of Vieques. We sequenced the mitochondrial DNA (mtDNA) control region of all of the samples and the complete mitogenomes of 12 of them to infer their putative place of origin. In addition, we genotyped 121 male samples for 25 Y-chromosome single nucleotide polymorphism and 17 STR loci. Approximately 60% of the participants had indigenous mtDNA haplotypes (mostly from haplogroups A2 and C1), while 25% had African and 15% European haplotypes. Three A2 sublineages were unique to the Greater Antilles, one of which was similar to Mesoamerican types, while C1b haplogroups showed links to South America, suggesting that people reached the island from the two distinct continental source areas. However, none of the male participants had indigenous Y-chromosomes, with 85% of them instead being European/Mediterranean and 15% sub-Saharan African in origin. West Eurasian Y-chromosome short tandem repeat haplotypes were quite diverse and showed similarities to those observed in southern Europe, North Africa and the Middle East. These results attest to the distinct, yet equally complex, pasts for the male and female ancestors of modern day Puerto Ricans.

Human Y-chromosome SNP characterization by multiplex amplified product-length polymorphism analysis

We designed an allele-specific amplification protocol to optimize Y-chromosome SNP typing, which is an unavoidable step for defining the phylogenetic status of paternal lineages. It allows the simultaneous highly specific definition of up to six mutations in a single reaction by amplification fragment length polymorphism (AFLP) without the need of specialized equipment, at a considerably lower cost than that based on single-base primer extension (SNaPshot™) technology or PCR-RFLP systems, requiring as little as 0.5 ng DNA and compatible with the small fragments characteristic of low-quality DNA. By designation of two primers recognizing the derived and ancestral state for each SNP, which can be differentiated by size by the addition of a noncomplementary nucleotide tail, we could define major Y clades E, F, K, R, Q, and subhaplogroups R1, R1a, R1b, R1b1b, R1b1c, J1, J2, G1, G2, I1, Q1a3, and Q1a3a1 through amplification fragments that ranged between 60 and 158bp.

Phylogenetic applications of whole Y-chromosome sequences and the Near Eastern origin of Ashkenazi Levites

Previous Y-chromosome studies have demonstrated that Ashkenazi Levites, members of a paternally inherited Jewish priestly caste, display a distinctive founder event within R1a, the most prevalent Y-chromosome haplogroup in Eastern Europe. Here we report the analysis of 16 whole R1 sequences and show that a set of 19 unique nucleotide substitutions defines the Ashkenazi R1a lineage. While our survey of one of these, M582, in 2,834 R1a samples reveals its absence in 922 Eastern Europeans, we show it is present in all sampled R1a Ashkenazi Levites, as well as in 33.8% of other R1a Ashkenazi Jewish males and 5.9% of 303 R1a Near Eastern males, where it shows considerably higher diversity. Moreover, the M582 lineage also occurs at low frequencies in non-Ashkenazi Jewish populations. In contrast to the previously suggested Eastern European origin for Ashkenazi Levites, the current data are indicative of a geographic source of the Levite founder lineage in the Near East and its likely presence among pre-Diaspora Hebrews.

Population genetics studies continue to debate whether Ashkenazi Levites originated in Europe or the Near East. Here, Rootsi et al. use whole Y-chromosome DNA sequences to unravel the phylogenetic origin of the Ashkenazi Levite and suggest an origin for the Levite founder lineage in the Near East.

Association between Y haplogroups and autosomal AIMs reveals intra-population substructure in Bolivian populations

For the correct evaluation of the weight of genetic evidence in a forensic context, databases must reflect the structure of the population, with all possible groups being represented. Countries with a recent history of admixture between strongly differentiated populations are usually highly heterogeneous and sub-structured. Bolivia is one of these countries, with a high diversity of ethnic groups and different levels of admixture (among Native Americans, Europeans and Africans) across the territory. For a better characterization of the male lineages in Bolivia, 17 Y-STR and 42 Y-SNP loci were genotyped in samples from La Paz and Chuquisaca. Only European and Native American Y-haplogroups were detected, and no sub-Saharan African chromosomes were found. Significant differences were observed between the two samples, with a higher frequency of European lineages in Chuquisaca than in La Paz. A sample belonging to haplogroup Q1a3a1a1-M19 was detected in La Paz, in a haplotype background different from those previously found in Argentina. This result supports an old M19 North-south dispersion in South America, possibly via two routes. When comparing the ancestry of each individual assessed through his Y chromosome with the one estimated using autosomal AIMs, (a) increased European ancestry in individuals with European Y chromosomes and (b) higher Native American ancestry in the carriers of Native American Y-haplogroups were observed, revealing an association between autosomal and Y-chromosomal markers. The results of this study demonstrate that a sub-structure does exist in Bolivia at both inter- and intrapopulation levels, a fact which must be taken into account in the evaluation of forensic genetic evidence.

The phylogenetic and geographic structure of Y-chromosome haplogroup R1a

R1a-M420 is one of the most widely spread Y-chromosome haplogroups; however, its substructure within Europe and Asia has remained poorly characterized. Using a panel of 16 244 male subjects from 126 populations sampled across Eurasia, we identified 2923 R1a-M420 Y-chromosomes and analyzed them to a highly granular phylogeographic resolution. Whole Y-chromosome sequence analysis of eight R1a and five R1b individuals suggests a divergence time of ∼25,000 (95% CI: 21,300-29,000) years ago and a coalescence time within R1a-M417 of ∼5800 (95% CI: 4800-6800) years. The spatial frequency distributions of R1a sub-haplogroups conclusively indicate two major groups, one found primarily in Europe and the other confined to Central and South Asia. Beyond the major European versus Asian dichotomy, we describe several younger sub-haplogroups. Based on spatial distributions and diversity patterns within the R1a-M420 clade, particularly rare basal branches detected primarily within Iran and eastern Turkey, we conclude that the initial episodes of haplogroup R1a diversification likely occurred in the vicinity of present-day Iran.

An unbiased resource of novel SNP markers provides a new chronology for the human Y chromosome and reveals a deep phylogenetic structure in Africa

Sequence diversity and the ages of the deepest nodes of the MSY phylogeny remain largely unexplored due to the severely biased collection of SNPs available for study. We characterized 68 worldwide Y chromosomes by high-coverage next-generation sequencing, including 18 deep-rooting ones, and identified 2386 SNPs, 80% of which were novel. Many aspects of this pool of variants resembled the pattern observed among genome-wide de novo events, suggesting that in the MSY, a large proportion of newly arisen alleles has survived in the phylogeny. Some degree of purifying selection emerged in the form of an excess of private missense variants. Our tree recapitulated the previously known topology, but the relative lengths of major branches were drastically modified and the associated node ages were remarkably older. We found significantly different branch lengths when comparing the rare deep-rooted A1b African lineage with the rest of the tree. Our dating results and phylogeography led to the following main conclusions: (1) Patrilineal lineages with ages approaching those of early AMH fossils survive today only in central-western Africa; (2) only a few evolutionarily successful MSY lineages survived between 160 and 115 kya; and (3) an early exit out of Africa (before 70 kya), which fits recent western Asian archaeological evidence, should be considered. Our experimental design produced an unbiased resource of new MSY markers informative for the initial formation of the anatomically modern human gene pool, i.e., a period of our evolution that had been previously considered to be poorly accessible with paternally inherited markers.

Y-chromosomal DNA analysis in French male lineages

French population, despite of its crucial geographic location for repopulation movements of Europe across time, it has been insufficiently characterized at the genetic level, especially for Y-chromosomal DNA variation. In order to make a genetic structure characterization, we have analyzed the Y-chromosome diversity of 558 male individuals, scattered along 7 different French regions: Alsace (Strasbourg), Auvergne (Clermont-Ferrand), Bretagne (Rennes), Île-de-France (Paris), Midi-Pyrénées (Toulouse), Nord-Pas-de-Calais (Lille) and Provence-Alpes-Côte d'Azur (Marseille). A total of 17 Y-chromosome STRs and 27 Y-chromosome SNPs were genotyped for each individual. Even though we find that most of the individual populations in France were not differentiated from each other, Bretagne population shows population substructure, an important fact to be considered when establishing general population databases.

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.

The paternal perspective of the Slovenian population and its relationship with other populations

BACKGROUND

The Slovenian territory is geographically positioned between the Alps, the Adriatic Sea, the Pannonian basin and the Dinaric Mountains and, as such, has served as a passageway for different populations over different periods of time. Turbulent historic events and the diverse geography of the region have produced a diverse contemporary population whose genetic analysis could provide insight into past demographic events.

AIM

The aim of this study was to analyse Y-chromosome biallelic and STR markers in a Slovenian population from five different regions.

SUBJECTS AND METHODS

A total of 42 Y-chromosomal biallelic markers and 17 Y-STRs were genotyped in 399 individuals from five different Slovenian regions.

RESULTS

The analysis of Y-chromosome markers revealed 29 different haplogroups in the Slovenian population, with the most common being R1a1a, R1b, I2a1 and I1. Analysis of the genetic affiliations between different populations revealed strong affiliations of the Slovenian gene pool with West Slavic populations.

CONCLUSION

Analysis of Y-chromosomal markers in five Slovenian regions revealed a diverse genetic landscape. Slovenian population display close genetic affiliations with West Slavic populations. The homogenous genetic strata of the West Slavic populations and the Slovenian population suggest the existence of a common ancestral Slavic population in central European region.

Uniparental genetic heritage of belarusians: encounter of rare middle eastern matrilineages with a central European mitochondrial DNA pool

Ethnic Belarusians make up more than 80% of the nine and half million people inhabiting the Republic of Belarus. Belarusians together with Ukrainians and Russians represent the East Slavic linguistic group, largest both in numbers and territory, inhabiting East Europe alongside Baltic-, Finno-Permic- and Turkic-speaking people. Till date, only a limited number of low resolution genetic studies have been performed on this population. Therefore, with the phylogeographic analysis of 565 Y-chromosomes and 267 mitochondrial DNAs from six well covered geographic sub-regions of Belarus we strove to complement the existing genetic profile of eastern Europeans. Our results reveal that around 80% of the paternal Belarusian gene pool is composed of R1a, I2a and N1c Y-chromosome haplogroups – a profile which is very similar to the two other eastern European populations – Ukrainians and Russians. The maternal Belarusian gene pool encompasses a full range of West Eurasian haplogroups and agrees well with the genetic structure of central-east European populations. Our data attest that latitudinal gradients characterize the variation of the uniparentally transmitted gene pools of modern Belarusians. In particular, the Y-chromosome reflects movements of people in central-east Europe, starting probably as early as the beginning of the Holocene. Furthermore, the matrilineal legacy of Belarusians retains two rare mitochondrial DNA haplogroups, N1a3 and N3, whose phylogeographies were explored in detail after de novo sequencing of 20 and 13 complete mitogenomes, respectively, from all over Eurasia. Our phylogeographic analyses reveal that two mitochondrial DNA lineages, N3 and N1a3, both of Middle Eastern origin, might mark distinct events of matrilineal gene flow to Europe: during the mid-Holocene period and around the Pleistocene-Holocene transition, respectively.

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.

Updating the Y-chromosomal phylogenetic tree for forensic applications based on whole genome SNPs

The Y-chromosomal phylogenetic tree has a wide variety of important forensic applications and therefore it needs to be state-of-the-art. Nevertheless, since the last 'official' published tree many publications reported additional Y-chromosomal lineages and other phylogenetic topologies. Therefore, it is difficult for forensic scientists to interpret those reports and use an up-to-date tree and corresponding nomenclature in their daily work. Whole genome sequencing (WGS) data is useful to verify and optimise the current phylogenetic tree for haploid markers. The AMY-tree software is the first open access program which analyses WGS data for Y-chromosomal phylogenetic applications. Here, all published information is collected in a phylogenetic tree and the correctness of this tree is checked based on the first large analysis of 747 WGS samples with AMY-tree. The obtained result is one phylogenetic tree with all peer-reviewed reported Y-SNPs without the observed recurrent and ambiguous mutations. Nevertheless, the results showed that currently only the genomes of a limited set of Y-chromosomal (sub-)haplogroups is available and that many newly reported Y-SNPs based on WGS projects are false positives, even with high sequencing coverage methods. This study demonstrates the usefulness of AMY-tree in the process of checking the quality of the present Y-chromosomal tree and it accentuates the difficulties to enlarge this tree based on only WGS methods.

Y-chromosome diversity in modern Bulgarians: new clues about their ancestry

To better define the structure and origin of the Bulgarian paternal gene pool, we have examined the Y-chromosome variation in 808 Bulgarian males. The analysis was performed by high-resolution genotyping of biallelic markers and by analyzing the STR variation within the most informative haplogroups. We found that the Y-chromosome gene pool in modern Bulgarians is primarily represented by Western Eurasian haplogroups with ∼ 40% belonging to haplogroups E-V13 and I-M423, and 20% to R-M17. Haplogroups common in the Middle East (J and G) and in South Western Asia (R-L23*) occur at frequencies of 19% and 5%, respectively. Haplogroups C, N and Q, distinctive for Altaic and Central Asian Turkic-speaking populations, occur at the negligible frequency of only 1.5%. Principal Component analyses group Bulgarians with European populations, apart from Central Asian Turkic-speaking groups and South Western Asia populations. Within the country, the genetic variation is structured in Western, Central and Eastern Bulgaria indicating that the Balkan Mountains have been permeable to human movements. The lineage analysis provided the following interesting results: (i) R-L23* is present in Eastern Bulgaria since the post glacial period; (ii) haplogroup E-V13 has a Mesolithic age in Bulgaria from where it expanded after the arrival of farming; (iii) haplogroup J-M241 probably reflects the Neolithic westward expansion of farmers from the earliest sites along the Black Sea. On the whole, in light of the most recent historical studies, which indicate a substantial proto-Bulgarian input to the contemporary Bulgarian people, our data suggest that a common paternal ancestry between the proto-Bulgarians and the Altaic and Central Asian Turkic-speaking populations either did not exist or was negligible.

Introducing the Algerian mitochondrial DNA and Y-chromosome profiles into the North African landscape

North Africa is considered a distinct geographic and ethnic entity within Africa. Although modern humans originated in this Continent, studies of mitochondrial DNA (mtDNA) and Y-chromosome genealogical markers provide evidence that the North African gene pool has been shaped by the back-migration of several Eurasian lineages in Paleolithic and Neolithic times. More recent influences from sub-Saharan Africa and Mediterranean Europe are also evident. The presence of East-West and North-South haplogroup frequency gradients strongly reinforces the genetic complexity of this region. However, this genetic scenario is beset with a notable gap, which is the lack of consistent information for Algeria, the largest country in the Maghreb. To fill this gap, we analyzed a sample of 240 unrelated subjects from a northwest Algeria cosmopolitan population using mtDNA sequences and Y-chromosome biallelic polymorphisms, focusing on the fine dissection of haplogroups E and R, which are the most prevalent in North Africa and Europe respectively. The Eurasian component in Algeria reached 80% for mtDNA and 90% for Y-chromosome. However, within them, the North African genetic component for mtDNA (U6 and M1; 20%) is significantly smaller than the paternal (E-M81 and E-V65; 70%). The unexpected presence of the European-derived Y-chromosome lineages R-M412, R-S116, R-U152 and R-M529 in Algeria and the rest of the Maghreb could be the counterparts of the mtDNA H1, H3 and V subgroups, pointing to direct maritime contacts between the European and North African sides of the western Mediterranean. Female influx of sub-Saharan Africans into Algeria (20%) is also significantly greater than the male (10%). In spite of these sexual asymmetries, the Algerian uniparental profiles faithfully correlate between each other and with the geography.

AMY-tree: an algorithm to use whole genome SNP calling for Y chromosomal phylogenetic applications

Background

Due to the rapid progress of next-generation sequencing (NGS) facilities, an explosion of human whole genome data will become available in the coming years. These data can be used to optimize and to increase the resolution of the phylogenetic Y chromosomal tree. Moreover, the exponential growth of known Y chromosomal lineages will require an automatic determination of the phylogenetic position of an individual based on whole genome SNP calling data and an up to date Y chromosomal tree.

Results

We present an automated approach, ‘AMY-tree’, which is able to determine the phylogenetic position of a Y chromosome using a whole genome SNP profile, independently from the NGS platform and SNP calling program, whereby mistakes in the SNP calling or phylogenetic Y chromosomal tree are taken into account. Moreover, AMY-tree indicates ambiguities within the present phylogenetic tree and points out new Y-SNPs which may be phylogenetically relevant. The AMY-tree software package was validated successfully on 118 whole genome SNP profiles of 109 males with different origins. Moreover, support was found for an unknown recurrent mutation, wrong reported mutation conversions and a large amount of new interesting Y-SNPs.

Conclusions

Therefore, AMY-tree is a useful tool to determine the Y lineage of a sample based on SNP calling, to identify Y-SNPs with yet unknown phylogenetic position and to optimize the Y chromosomal phylogenetic tree in the future. AMY-tree will not add lineages to the existing phylogenetic tree of the Y-chromosome but it is the first step to analyse whole genome SNP profiles in a phylogenetic framework.

Paleo-Balkan and Slavic contributions to the genetic pool of Moldavians: insights from the Y chromosome

Moldova has a rich historical and cultural heritage, which may be reflected in the current genetic makeup of its population. To date, no comprehensive studies exist about the population genetic structure of modern Moldavians. To bridge this gap with respect to paternal lineages, we analyzed 37 binary and 17 multiallelic (STRs) polymorphisms on the non-recombining portion of the Y chromosome in 125 Moldavian males. In addition, 53 Ukrainians from eastern Moldova and 54 Romanians from the neighboring eastern Romania were typed using the same set of markers. In Moldavians, 19 Y chromosome haplogroups were identified, the most common being I-M423 (20.8%), R-M17* (17.6%), R-M458 (12.8%), E-v13 (8.8%), R-M269* and R-M412* (both 7.2%). In Romanians, 14 haplogroups were found including I-M423 (40.7%), R-M17* (16.7%), R-M405 (7.4%), E-v13 and R-M412* (both 5.6%). In Ukrainians, 13 haplogroups were identified including R-M17 (34.0%), I-M423 (20.8%), R-M269* (9.4%), N-M178, R-M458 and R-M73 (each 5.7%). Our results show that a significant majority of the Moldavian paternal gene pool belongs to eastern/central European and Balkan/eastern Mediterranean Y lineages. Phylogenetic and AMOVA analyses based on Y-STR loci also revealed that Moldavians are close to both eastern/central European and Balkan-Carpathian populations. The data correlate well with historical accounts and geographical location of the region and thus allow to hypothesize that extant Moldavian paternal genetic lineages arose from extensive recent admixture between genetically autochthonous populations of the Balkan-Carpathian zone and neighboring Slavic groups.

Contemporary paternal genetic landscape of Polish and German populations: from early medieval Slavic expansion to post-World War II resettlements

Homogeneous Proto-Slavic genetic substrate and/or extensive mixing after World War II were suggested to explain homogeneity of contemporary Polish paternal lineages. Alternatively, Polish local populations might have displayed pre-war genetic heterogeneity owing to genetic drift and/or gene flow with neighbouring populations. Although sharp genetic discontinuity along the political border between Poland and Germany indisputably results from war-mediated resettlements and homogenisation, it remained unknown whether Y-chromosomal diversity in ethnically/linguistically defined populations was clinal or discontinuous before the war. In order to answer these questions and elucidate early Slavic migrations, 1156 individuals from several Slavic and German populations were analysed, including Polish pre-war regional populations and an autochthonous Slavic population from Germany. Y chromosomes were assigned to 39 haplogroups and genotyped for 19 STRs. Genetic distances revealed similar degree of differentiation of Slavic-speaking pre-war populations from German populations irrespective of duration and intensity of contacts with German speakers. Admixture estimates showed minor Slavic paternal ancestry (~20%) in modern eastern Germans and hardly detectable German paternal ancestry in Slavs neighbouring German populations for centuries. BATWING analysis of isolated Slavic populations revealed that their divergence was preceded by rapid demographic growth, undermining theory that Slavic expansion was primarily linguistic rather than population spread. Polish pre-war regional populations showed within-group heterogeneity and lower STR variation within R-M17 subclades compared with modern populations, which might have been homogenised by war resettlements. Our results suggest that genetic studies on early human history in the Vistula and Oder basins should rely on reconstructed pre-war rather than modern populations.

MLVA based classification of Mycobacterium tuberculosis complex lineages for a robust phylogeographic snapshot of its worldwide molecular diversity

Multiple-locus variable-number tandem repeat analysis (MLVA) is useful to establish transmission routes and sources of infections for various microorganisms including Mycobacterium tuberculosis complex (MTC). The recently released SITVITWEB database contains 12-loci Mycobacterial Interspersed Repetitive Units – Variable Number of Tandem DNA Repeats (MIRU-VNTR) profiles and spoligotype patterns for thousands of MTC strains; it uses MIRU International Types (MIT) and Spoligotype International Types (SIT) to designate clustered patterns worldwide. Considering existing doubts on the ability of spoligotyping alone to reveal exact phylogenetic relationships between MTC strains, we developed a MLVA based classification for MTC genotypic lineages. We studied 6 different subsets of MTC isolates encompassing 7793 strains worldwide. Minimum spanning trees (MST) were constructed to identify major lineages, and the most common representative located as a central node was taken as the prototype defining different phylogenetic groups. A total of 7 major lineages with their respective prototypes were identified: Indo-Oceanic/MIT57, East Asian and African Indian/MIT17, Euro American/MIT116, West African-I/MIT934, West African-II/MIT664, M. bovis/MIT49, M.canettii/MIT60. Further MST subdivision identified an additional 34 sublineage MIT prototypes. The phylogenetic relationships among the 37 newly defined MIRU-VNTR lineages were inferred using a classification algorithm based on a bayesian approach. This information was used to construct an updated phylogenetic and phylogeographic snapshot of worldwide MTC diversity studied both at the regional, sub-regional, and country level according to the United Nations specifications. We also looked for IS6110 insertional events that are known to modify the results of the spoligotyping in specific circumstances, and showed that a fair portion of convergence leading to the currently observed bias in phylogenetic classification of strains may be traced back to the presence of IS6110. These results shed new light on the evolutionary history of the pathogen in relation to the history of peopling and human migration.

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.

Pasture names with Romance and Slavic roots facilitate dissection of Y chromosome variation in an exclusively German-speaking alpine region

The small alpine district of East Tyrol (Austria) has an exceptional demographic history. It was contemporaneously inhabited by members of the Romance, the Slavic and the Germanic language groups for centuries. Since the Late Middle Ages, however, the population of the principally agrarian-oriented area is solely Germanic speaking. Historic facts about East Tyrol's colonization are rare, but spatial density-distribution analysis based on the etymology of place-names has facilitated accurate spatial mapping of the various language groups' former settlement regions. To test for present-day Y chromosome population substructure, molecular genetic data were compared to the information attained by the linguistic analysis of pasture names. The linguistic data were used for subdividing East Tyrol into two regions of former Romance (A) and Slavic (B) settlement. Samples from 270 East Tyrolean men were genotyped for 17 Y-chromosomal microsatellites (Y-STRs) and 27 single nucleotide polymorphisms (Y-SNPs). Analysis of the probands' surnames revealed no evidence for spatial genetic structuring. Also, spatial autocorrelation analysis did not indicate significant correlation between genetic (Y-STR haplotypes) and geographic distance. Haplogroup R-M17 chromosomes, however, were absent in region A, but constituted one of the most frequent haplogroups in region B. The R-M343 (R1b) clade showed a marked and complementary frequency distribution pattern in these two regions. To further test East Tyrol's modern Y-chromosomal landscape for geographic patterning attributable to the early history of settlement in this alpine area, principal coordinates analysis was performed. The Y-STR haplotypes from region A clearly clustered with those of Romance reference populations and the samples from region B matched best with Germanic speaking reference populations. The combined use of onomastic and molecular genetic data revealed and mapped the marked structuring of the distribution of Y chromosomes in an alpine region that has been culturally homogeneous for centuries.

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.

Ancient migratory events in the Middle East: new clues from the Y-chromosome variation of modern Iranians

Knowledge of high resolution Y-chromosome haplogroup diversification within Iran provides important geographic context regarding the spread and compartmentalization of male lineages in the Middle East and southwestern Asia. At present, the Iranian population is characterized by an extraordinary mix of different ethnic groups speaking a variety of Indo-Iranian, Semitic and Turkic languages. Despite these features, only few studies have investigated the multiethnic components of the Iranian gene pool. In this survey 938 Iranian male DNAs belonging to 15 ethnic groups from 14 Iranian provinces were analyzed for 84 Y-chromosome biallelic markers and 10 STRs. The results show an autochthonous but non-homogeneous ancient background mainly composed by J2a sub-clades with different external contributions. The phylogeography of the main haplogroups allowed identifying post-glacial and Neolithic expansions toward western Eurasia but also recent movements towards the Iranian region from western Eurasia (R1b-L23), Central Asia (Q-M25), Asia Minor (J2a-M92) and southern Mesopotamia (J1-Page08). In spite of the presence of important geographic barriers (Zagros and Alborz mountain ranges, and the Dasht-e Kavir and Dash-e Lut deserts) which may have limited gene flow, AMOVA analysis revealed that language, in addition to geography, has played an important role in shaping the nowadays Iranian gene pool. Overall, this study provides a portrait of the Y-chromosomal variation in Iran, useful for depicting a more comprehensive history of the peoples of this area as well as for reconstructing ancient migration routes. In addition, our results evidence the important role of the Iranian plateau as source and recipient of gene flow between culturally and genetically distinct populations.

Ancestral modal Y-STR haplotype shared among Romani and South Indian populations

One of the primary unanswered questions regarding the dispersal of Romani populations concerns the geographical region and/or the Indian caste/tribe that gave rise to the proto-Romani group. To shed light on this matter, 161 Y-chromosomes from Roma, residing in two different provinces of Serbia, were analyzed. Our results indicate that the paternal gene pool of both groups is shaped by several strata, the most prominent of which, H1-M52, comprises almost half of each collection's patrilineages. The high frequency of M52 chromosomes in the two Roma populations examined may suggest that they descend from a single founder that has its origins in the Indian subcontinent. Moreover, when the Y-STR profiles of haplogroup H derived individuals in our Roma populations were compared to those typed in the South Indian emigrants from Malaysia and groups from Madras, Karnataka (Lingayat and Vokkaliga castes) and tribal Soligas, sharing of the two most common haplotypes was observed. These similarities suggest that South India may have been one of the contributors to the proto-Romanis. European genetic signatures (i.e., haplogroups E1b1b1a1b-V13, G2a-P15, I-M258, J2-M172 and R1-M173), on the other hand, were also detected in both groups, but at varying frequencies. The divergent European genetic signals in each collection are likely the result of differential gene flow and/or admixture with the European host populations but may also be attributed to dissimilar endogamous practices following the initial founder effect. Our data also support the notion that a number of haplogroups including G2a-P15, J2a3b-M67(xM92), I-M258 and E1b1b1-M35 were incorporated into the proto-Romani paternal lineages as migrants moved from northern India through Southwestern Asia, the Middle East and/or Anatolia into the Balkans.