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

The Genetic Variability of Present-Day Bulgarians Captures Ancient and Recent Ancestral Contributions

OBJECTIVES

Thanks to its pivotal crossroad position, Bulgaria played a fundamental key role during all the migration processes that interested the continent through time. While the genetic variability of the country has been deeply investigated using uniparental markers, previous genome-wide autosomal-based surveys mainly consisted of wider-range analyzes on Europe and the whole Balkan Peninsula. Here, we specifically focused on the Bulgarian population to recapitulate the main patterns of genomic variation and the major events that shaped the present-day genetic landscape.

METHODS

A total of 112 samples from seven highly representative areas of present-day Bulgaria were collected and genotyped for approximately 720 K genome-wide SNPs, and integrated with previously generated genomic data from wide modern and ancient reference panels to explore fine-scale relationship patterns and detail ancestral contributions.

RESULTS

In addition to the combination of ancient ancestries related to the early Mesolithic hunter-gatherers, Neolithic farmers, and Bronze Age Steppe pastoralists, both haplotype-based analyzes on modern populations and the comparisons with ancient genomes suggest the contribution of population processes that have occurred after the Roman rule and during the Medieval period in shaping the current Bulgarian genetic pool.

CONCLUSIONS

Our results align with previous evidence highlighting the impact that some historical events may have had not only in contributing to the ethnical and socio-cultural richness of present-day populations, but also in participating in the formation of the current genomic landscape. By providing new data from modern highly-representative samples, this study integrates further research to provide a comprehensive overview of the genetic history of Bulgaria.

Phylogeography of Y-chromosome haplogroup I-P37.2 in Serbian population groups originating from distinct parts of the Balkan Peninsula

Genetic structure of the contemporary Serbian population was shaped by a long history of turbulent historical and demographical events. The most important migrations of Serbs towards present day Serbia, in the recent history, occurred between the 15th to the 18th century from the regions of Old Herzegovina and Kosovo and Metohija. Previous haplogroup analysis revealed wide spectrum of main haplogroups, among which haplogroup I-P37.2 was the most frequent one in Serbian population groups originating from the Balkan Peninsula. Within this study 464 Serbian males samples from three geographical regions in the Balkan Peninsula inhabited by Serbs: present-day Serbia, regions of Old Herzegovina and Kosovo and Metohija, previously classified as haplogroup I-P37.2, were genotyped using the 22 Y-single nucleotide polymorphisms (Y-SNPs) in order to determine deeper phylogenetic and phylogeographic analysis of haplogroup I-P37.2. Based on SNP typing all samples in the Old Herzegovina and present-day Serbia dataset and 122 out of 128 samples from Kosovo and Metohija were assigned to haplogroup I-L621. Further SNP typing revealed very similar haplogroup distribution in all datasets, with the predominant haplogroup being I-PH908, followed by haplogroup I-Z17855. Analysis within haplogroup I-PH908 distinguished haplogroup I-FT14506 as the most frequent in the Kosovo and Metohija dataset, while haplogroup I-FT16449 was the most frequent in the Old Herzegovina dataset. In the present-day Serbia dataset, occurrence of haplogroups I-FT14506 and I-FT16449 was almost equal, comprising 40.2 % and 34.4 %, respectively. Low level of differentiation, within haplogroup I-PH908, was detected between all datasets, with the lowest one detected between present-day Serbia and Old Herzegovina datasets and highest one between Kosovo and Metohija and Old Herzegovina datasets. Furthermore, median-joining network analysis and shared haplotypes statistics revealed closer genetic relationship between Old Herzegovina and present-day Serbia haplotypes. Results obtained within this study support the thesis that migrations from historical region of Old Herzegovina and geographical region of Kosovo and Metohija, had great contribution on the present-day Serbian population genetic structure. Furthermore, here presented results, gave insight into geographic distribution of detected haplogroups I-Z17855, I-Y4460, I-PH908, I-Y5596, I-Y4882, I-FT14506, I-FT16449 and I-A5913 and analyzed SNPs, enabling further improvement of the geographic resolution of paternal ancestry inference.

A Review of Genetic Diversity Based on the Y Chromosome in the Romanian Population

Y chromosome analysis is used in a number of practical applications, including investigations of criminal cases, establishment of paternity, searching for missing persons, studies on human migration, evolutionary research, and historical and genealogical investigations. Questions about the origin of individual ethnic groups are addressed not only through archaeological, linguistic, and ethnographic methods but also through molecular genetics methods. The study of genetic diversity in Romania is particularly interesting from several perspectives because Romania, located in Southeast Europe, is distinguished by the fact that the Carpathians and the Danube served as natural barriers against the migrations of peoples for centuries, thus influencing the genetic mixture of the population. This is relevant for understanding the history and formation of ethnic groups in the region. In addition, many ethnic minorities live in Romania, which adds an additional dimension of genetic and cultural diversity. This article aims to provide an updated picture of the genetic diversity in Romania and to highlight the significant studies carried out among the Romanian population. By analyzing the articles published in the Web of Science, Scopus, or PubMed databases, which explore genetic diversity using the Y chromosome, the aim is to better understand the current genetic panorama in Romania.

Distribution of Y-chromosome haplogroups in Serbian population groups originating from historically and geographically significant distinct parts of the Balkan Peninsula

Our study enrolled 1200 Serbian males originating from three geographical regions in the Balkan Peninsula inhabited by Serbs: present-day Serbia, regions of Old Herzegovina and Kosovo and Metohija. These samples were genotyped using the combination of 23 Y-chromosomal short tandem repeats (Y-STRs) loci and 17 Ychromosomal single nucleotide polymorphisms (Y-SNPs) loci for the haplotype and haplogroup analysis in order to characterize in detail Y chromosome flow in the recent history. Serbia's borders have changed through history, forcing Serbs constantly to migrate to different regions of Balkan Peninsula. The most significant migration waves in the recent history towards present-day Serbia occurred from the regions of Old- Herzegovina and Kosovo and Metohija that lie in the south-west/south. High haplotype diversity and discrimination capacity were observed in all three datasets, with the highest number of unique haplotypes (381) and discrimination capacity (0.97) detected in the samples originating from the present-day Serbia. Haplogroup composition didn't differ significantly among datasets, with three dominant haplogroups (I-M170, E-P170 and R-M198), and haplogroup I-M170 being the most frequent in all three datasets. Haplogroup E-P170 was the second most dominant in the dataset originating from geographical region of Kosovo and Metohija, whereas haplogroup R-M198 was the second most prevalent in the dataset from historical region of Old Herzegovina. Based on the phylogenetic three for haplogroup I constructed within this study, haplogroup I2a1-P37.2 was the most dominant within all three datasets, especially in the dataset from historical region of Old Herzegovina, where 182 out of 400 samples were derived for SNP P37.2. Genetic distances between three groups of samples, evaluated by the Fst and Rst statistical values, and further visualized through multidimensional scaling plot, showed great genetic similarity between datasets from Old Herzegovina and present-day Serbia. Genetic difference in the haplogroup distribution and frequency between datasets from historical region of Old Herzegovina and from geographical region of Kosovo and Metohija was confirmed with highest Fst and Rst vaules. In this study we have distinguished genetic structure, diversity and haplogroup frequencies within 1200 Serbian males from three datasets, relationships among them as well as with other Balkan and European populations, which is useful for studying recent demographic history.

Genetic insights into the paternal admixture history of Chinese Mongolians via high-resolution customized Y-SNP SNaPshot panels

The Mongolian people, one of the Mongolic-speaking populations, are native to the Mongolian Plateau in North China and southern Siberia. Many ancient DNA studies recently reported extensive population transformations during the Paleolithic to historic periods in this region, while little is known about the paternal genetic legacy of modern geographically different Mongolians. Here, we genotyped 215 Y-chromosomal single nucleotide polymorphisms (Y-SNPs) and 37 Y-chromosomal short tandem repeats (Y-STRs) among 679 Mongolian individuals from Hohhot, Hulunbuir, and Ordos in North China using the AGCU Y37 kit and our developed eight Y-SNP SNaPshot panels (including two panels first reported herein). The C-M130 Y-SNP SNaPshot panel defines 28 subhaplogroups, and the N/O/Q complementary Y-SNP SNaPshot panel defines 30 subhaplogroups of N1b-F2930, N1a1a1a1a3-B197, Q-M242, and O2a2b1a1a1a4a-CTS4658, which improved the resolution our developed Y-SNP SNaPshot panel set and could be applied for dissecting the finer-scale paternal lineages of Mongolic speakers. We found a strong association between Mongolian-prevailing haplogroups and some observed microvariants among the newly generated Y-STR haplotype data, suggesting the possibility of haplogroup prediction based on the distribution of Y-STR haplotypes. We identified three main ancestral sources of the observed Mongolian-dominant haplogroups, including the local lineage of C2*-M217 and incoming lineages from other regions of southern East Asia (O2*-M122, O1b*-P31, and N1*-CTS3750) and western Eurasia (R1*-M173). We also observed DE-M145, D1*-M174, C1*-F3393, G*-M201, I-M170, J*-M304, L-M20, O1a*-M119, and Q*-M242 at relatively low frequencies (< 5.00%), suggesting a complex admixture history between Mongolians and other incoming Eurasians from surrounding regions. Genetic clustering analyses indicated that the studied Mongolians showed close genetic affinities with other Altaic-speaking populations and Sinitic-speaking Hui people. The Y-SNP haplotype/haplogroup-based genetic legacy not only revealed that the stratification among geographically/linguistically/ethnically different Chinese populations was highly consistent with the geographical division and language classification, but also demonstrated that patrilineal genetic materials could provide fine-scale genetic structures among geographically different Mongolian people, suggesting that our developed high-resolution Y-SNP SNaPshot panels have the potential for forensic pedigree searches and biogeographical ancestry inference.

Phylogeographic review of Y chromosome haplogroups in Europe

The Y chromosome has been widely explored for the study of human migrations. Due to its paternal inheritance, the Y chromosome polymorphisms are helpful tools for understanding the geographical distribution of populations all over the world and for inferring their origin, which is really useful in forensics. The remarkable historical context of Europe, with numerous migrations and invasions, has turned this continent into a melting pot. For this reason, it is interesting to study the Y chromosome variability and how it has contributed to improving our knowledge of the distribution and development of European male genetic pool as it is today. The analysis of Y lineages in Europe shows the predominance of four haplogroups, R1b-M269, I1-M253, I2-M438 and R1a-M420. However, other haplogroups have been identified which, although less frequent, provide significant evidence about the paternal origin of the populations. In addition, the study of the Y chromosome in Europe is a valuable tool for revealing the genetic trace of the different European colonizations, mainly in several American countries, where the European ancestry is mostly detected by the presence of the R1b-M269 haplogroup. Therefore, the objective of this review is to compile the studies of the Y chromosome haplogroups in current European populations, in order to provide an outline of these haplogroups which facilitate their use in forensic studies.

The association of gr/gr deletion in the Y chromosome and impaired spermatogenesis in Bulgarian males: a pilot study

Background

The microdeletions of the Y chromosome are associated with a decreased number of sperm cells in the ejaculate and male infertility. One such deletion is the gr/gr, which leads to reduction of the alleles in the DAZ gene. In order to evaluate its role, 30 patients with sperm count below 5 × 106/mL and 30 normal fertile men were genotyped for microdeletions in Azoospermia Factor (AZF) region, including gr/gr deletion.

Results

Twenty-one men (70.00%) from the males with impaired fertility had a normal genotype, seven men (23.33%) were carriers of a gr/gr deletion, and two men (6.67%) demonstrated a deletion of the AZFc region. From the control group, there were 22 men (73.33%) with no deletion, seven (23.33%) carried a gr/gr deletion, and one man (3.33%) carried a possible deletion in the b1-b4 region of the Y chromosome. The gr/gr deletion of the Y chromosome showed a statistically nonsignificant difference (p = 0.94, alpha < 0.05) with the control group (chi-square, degrees of freedom 0.006).

Conclusion

The results from the research show no association between this deletion and male infertility in the studied population. Despite the small sample size, the high frequency of the gr/gr deletion, found in the control group, suggests that this variant may not be associated with impaired spermatogenesis.

Y-chromosome and Surname Analyses for Reconstructing Past Population Structures: The Sardinian Population as a Test Case

Many anthropological, linguistic, genetic and genomic analyses have been carried out to evaluate the potential impact that evolutionary forces had in shaping the present-day Sardinian gene pool, the main outlier in the genetic landscape of Europe. However, due to the homogenizing effect of internal movements, which have intensified over the past fifty years, only partial information has been obtained about the main demographic events. To overcome this limitation, we analyzed the male-specific region of the Y chromosome in three population samples obtained by reallocating a large number of Sardinian subjects to the place of origin of their monophyletic surnames, which are paternally transmitted through generations in most of the populations, much like the Y chromosome. Three Y-chromosome founding lineages, G2-L91, I2-M26 and R1b-V88, were identified as strongly contributing to the definition of the outlying position of Sardinians in the European genetic context and marking a significant differentiation within the island. The present distribution of these lineages does not always mirror that detected in ancient DNAs. Our results show that the analysis of the Y-chromosome gene pool coupled with a sampling method based on the origin of the family name, is an efficient approach to unravelling past heterogeneity, often hidden by recent movements, in the gene pool of modern populations. Furthermore, the reconstruction and comparison of past genetic isolates represent a starting point to better assess the genetic information deriving from the increasing number of available ancient DNA samples.

Ancient human mitochondrial genomes from Bronze Age Bulgaria: new insights into the genetic history of Thracians

One of the best documented Indo-European civilizations that inhabited Bulgaria is the Thracians, who lasted for more than five millennia and whose origin and relationships with other past and present-day populations are debated among researchers. Here we report 25 new complete mitochondrial genomes of ancient individuals coming from three necropolises located in different regions of Bulgaria – Shekerdja mogila, Gabrova mogila and Bereketska mogila – dated to II-III millennium BC. The identified mtDNA haplogroup composition reflects the mitochondrial variability of Western Eurasia. In particular, within the ancient Eurasian genetic landscape, Thracians locate in an intermediate position between Early Neolithic farmers and Late Neolithic-Bronze Age steppe pastoralists, supporting the scenario that the Balkan region has been a link between Eastern Europe and the Mediterranean since the prehistoric time. Spatial Principal Component Analysis (sPCA) performed on Thracian and modern mtDNA sequences, confirms the pattern highlighted on ancient populations, overall indicating that the maternal gene pool of Thracians reflects their central geographical position at the gateway of Europe.

Analysis of the human Y-chromosome haplogroup Q characterizes ancient population movements in Eurasia and the Americas

BACKGROUND

Recent genome studies of modern and ancient samples have proposed that Native Americans derive from a subset of the Eurasian gene pool carried to America by an ancestral Beringian population, from which two well-differentiated components originated and subsequently mixed in different proportion during their spread in the Americas. To assess the timing, places of origin and extent of admixture between these components, we performed an analysis of the Y-chromosome haplogroup Q, which is the only Pan-American haplogroup and accounts for virtually all Native American Y chromosomes in Mesoamerica and South America.

RESULTS

Our analyses of 1.5 Mb of 152 Y chromosomes, 34 re-sequenced in this work, support a "coastal and inland routes scenario" for the first entrance of modern humans in North America. We show a major phase of male population growth in the Americas after 15 thousand years ago (kya), followed by a period of constant population size from 8 to 3 kya, after which a secondary sign of growth was registered. The estimated dates of the first expansion in Mesoamerica and the Isthmo-Colombian Area, mainly revealed by haplogroup Q-Z780, suggest an entrance in South America prior to 15 kya. During the global constant population size phase, local South American hints of growth were registered by different Q-M848 sub-clades. These expansion events, which started during the Holocene with the improvement of climatic conditions, can be ascribed to multiple cultural changes rather than a steady population growth and a single cohesive culture diffusion as it occurred in Europe.

CONCLUSIONS

We established and dated a detailed haplogroup Q phylogeny that provides new insights into the geographic distribution of its Eurasian and American branches in modern and ancient samples.

Reconstructing the genetic history of Italians: new insights from a male (Y-chromosome) perspective

BACKGROUND

Due to its central and strategic position in Europe and in the Mediterranean Basin, the Italian Peninsula played a pivotal role in the first peopling of the European continent and has been a crossroad of peoples and cultures since then.

AIM

This study aims to gain more information on the genetic structure of modern Italian populations and to shed light on the migration/expansion events that led to their formation.

SUBJECTS AND METHODS

High resolution Y-chromosome variation analysis in 817 unrelated males from 10 informative areas of Italy was performed. Haplogroup frequencies and microsatellite haplotypes were used, together with available data from the literature, to evaluate Mediterranean and European inputs and date their arrivals.

RESULTS

Fifty-three distinct Y-chromosome lineages were identified. Their distribution is in general agreement with geography, southern populations being more differentiated than northern ones.

CONCLUSIONS

A complex genetic structure reflecting the multifaceted peopling pattern of the Peninsula emerged: southern populations show high similarity with those from the Middle East and Southern Balkans, while those from Northern Italy are close to populations of North-Western Europe and the Northern Balkans. Interestingly, the population of Volterra, an ancient town of Etruscan origin in Tuscany, displays a unique Y-chromosomal genetic structure.

Ancient and recent Middle Eastern maternal genetic contribution to North Africa as viewed by mtDNA diversity in Tunisian Arab populations

OBJECTIVES

Through previous mitochondrial DNA studies, the Middle Eastern maternal genetic contribution to Tunisian populations appears limited. In fact, most of the studied communities were cosmopolitan, or of Berber or Andalusian origin. To provide genetic evidence for the actual contribution of Middle Eastern mtDNA lineages to Tunisia, we focused on two Arab speaking populations from Kairouan and Wesletia known to belong to an Arab genealogical lineage.

MATERIALS AND METHODS

A total of 114 samples were sequenced for the mtDNA HVS-I and HVS-II regions. Using these data, we evaluated the distribution of Middle Eastern haplogroups in the study populations, constructed interpolation maps, and established phylogenetic networks allowing estimation of the coalescence time for three specific Middle Eastern subclades (R0a, J1b, and T1).

RESULTS

Both studied populations displayed North African genetic structure and Middle Eastern lineages with a frequency of 12% and 28.12% in Kairouan and Wesletia, respectively. TMRCA estimates for haplogroups T1a, R0a, and J1b in Tunisian Arabian samples were around 15 000 YBP, 9000 to 5000 YBP, and 960 to 600 YBP, respectively.

CONCLUSIONS

The Middle Eastern maternal genetic contribution to Tunisian populations, as to other North African populations, occurred mostly in deep prehistory. They were brought in different migration waves during the Upper Paleolithic, probably with the expansion of Iberomaurusian culture, and during Epipaleolithic and Early Neolithic periods, which are concomitant with the Capsian civilization. Middle Eastern lineages also came to Tunisia during the recent Islamic expansion of the 7th CE and the subsequent massive Bedouin migration during the 11th CE.

Whole Y-chromosome sequences reveal an extremely recent origin of the most common North African paternal lineage E-M183 (M81)

E-M183 (E-M81) is the most frequent paternal lineage in North Africa and thus it must be considered to explore past historical and demographical processes. Here, by using whole Y chromosome sequences from 32 North African individuals, we have identified five new branches within E-M183. The validation of these variants in more than 200 North African samples, from which we also have information of 13 Y-STRs, has revealed a strong resemblance among E-M183 Y-STR haplotypes that pointed to a rapid expansion of this haplogroup. Moreover, for the first time, by using both SNP and STR data, we have provided updated estimates of the times-to-the-most-recent-common-ancestor (TMRCA) for E-M183, which evidenced an extremely recent origin of this haplogroup (2,000-3,000 ya). Our results also showed a lack of population structure within the E-M183 branch, which could be explained by the recent and rapid expansion of this haplogroup. In spite of a reduction in STR heterozygosity towards the West, which would point to an origin in the Near East, ancient DNA evidence together with our TMRCA estimates point to a local origin of E-M183 in NW Africa.

Clinical utility of the low-density Infinium QC genotyping Array in a genomics-based diagnostics laboratory

Background

With 15,949 markers, the low-density Infinium QC Array-24 BeadChip enables linkage analysis, HLA haplotyping, fingerprinting, ethnicity determination, mitochondrial genome variations, blood groups and pharmacogenomics. It represents an attractive independent QC option for NGS-based diagnostic laboratories, and provides cost-efficient means for determining gender, ethnic ancestry, and sample kinships, that are important for data interpretation of NGS-based genetic tests.

Methods

We evaluated accuracy and reproducibility of Infinium QC genotyping calls by comparing them with genotyping data of the same samples from other genotyping platforms, whole genome/exome sequencing. Accuracy and robustness of determining gender, provenance, and kinships were assessed.

Results

Concordance of genotype calls between Infinium QC and other platforms was above 99%. Here we show that the chip’s ancestry informative markers are sufficient for ethnicity determination at continental and sometimes subcontinental levels, with assignment accuracy varying with the coverage for a particular region and ethnic groups. Mean accuracies of provenance prediction at a regional level were varied from 81% for Asia, to 89% for Americas, 86% for Africa, 97% for Oceania, 98% for Europe, and 100% for India. Mean accuracy of ethnicity assignment predictions was 63%. Pairwise concordances of AFR samples with the samples from any other super populations were the lowest (0.39–0.43), while the concordances within the same population were relatively high (0.55–0.61). For all populations except African, cross-population comparisons were similar in their concordance ranges to the range of within-population concordances (0.54–0.57). Gender determination was correct in all tested cases.

Conclusions

Our results indicate that the Infinium QC Array-24 chip is suitable for cost-efficient, independent QC assaying in the settings of an NGS-based molecular diagnostic laboratory; hence, we recommend its integration into the standard laboratory workflow. Low-density chips can provide sample-specific measures for variant call accuracy, prevent sample mix-ups, validate self-reported ethnicities, and detect consanguineous cases. Integration of low-density chips into QC procedures aids proper interpretation of candidate sequence variants. To enhance utility of this low-density chip, we recommend expansion of ADME and mitochondrial markers. Inexpensive Infinium-like low-density human chips have a potential to become a “Swiss army knife” among genotyping assays suitable for many applications requiring high-throughput assays.

Electronic supplementary material

The online version of this article (10.1186/s12920-017-0297-7) contains supplementary material, which is available to authorized users.

A study of the Bodrogköz population in north-eastern Hungary by Y chromosomal haplotypes and haplogroups

We have determined the distribution of Y chromosomal haplotypes and haplogroups in population samples from one of the most important areas in north-eastern Hungary from many villages in the Bodrogköz. The Bodrogköz region was chosen due to its isolated nature, because this area was a moorland encircled by the Tisza, Bodrog, and Latorca Rivers and inhabitants of this part of Hungary escaped from both Tatar and Ottoman invasions, which decimated the post-Hungarian Conquest populations in many parts of the country. Furthermore, in the first half of the tenth century, this region served as the Palatial Centre and burial grounds of the Hungarian tribes. It has thus been assumed that the present population in this area is likely to be more similar to the population that lived in the Conquest period. We analysed male-specific markers, 23 Y-STRs and more than 30 Y-SNPs, that reflect the past and recent genetic history. We found that the general haplogroup distribution of the samples showed high genetic similarity to non-Bodrogköz Hungarians and neighbouring populations, despite its sheltered location and historical record. We were able to classify the Y-chromosomal haplogroups into four large groups based on STR mutation events: pre-Roman/Roman ancient lineage, Finno-Ugric speakers arriving into the Carpathian Basin, Migration period admixture, and post-Hungarian Conquest admixture. It is clear that a significantly larger database with deep haplogroup resolution, including ancient DNA data, is required to strengthen this research.

Genetic heritage of Croatians in the Southeastern European gene pool-Y chromosome analysis of the Croatian continental and Island population

OBJECTIVES

The research objective of this study is to enlarge and deepen the Y chromosome research on the Croatian population and enable additional insights into the population diversity and historic events that shaped the current genetic landscape of Croatia and Southeastern Europe (SEE).

MATERIALS AND METHODS

A high-resolution phylogenetic and phylogeographic analysis of 66 biallelic (SNPs) and 17 microsatellite (STRs) markers of the Y chromosome was performed using 720 Croatian samples. The obtained results were placed in a wider European context by comparison with ∼4450 samples from a number of other European populations.

RESULTS

A high diversity of haplogroups was observed in the overall Croatian sample, and all typical European Y chromosome haplogroups with corresponding clinal patterns were observed. Three distinct genetic signals were identifiable in the Croatian paternal gene pool - I2a1b-M423, R1a1a1b1a*-M558, and E1b1b1a1b1a-V13 haplogroups.

DISCUSSION

The analyses of the dominant and autochthonous I2a1b-M423 lineage (>30%) suggest that SEE had a significant role in the Upper Paleolithic, the R1a1a1b1a*-M558 lineage (19%) represents a signal from present day Slavic populations of Central Europe in the Croatian population, and the phylogeography of the E1b1b1a1b1a-V13 clade (around 9%) implies cultural diffusion of agriculture into Europe via the Balkan Peninsula. Am. J. Hum. Biol., 2016. © 2016 Wiley Periodicals, Inc. Am. J. Hum. Biol. 28:837-845, 2016. © 2016Wiley Periodicals, Inc.

Ancient Human Migration after Out-of-Africa

The serial founder model of modern human origins predicts that the phylogeny of ancestries exhibits bifurcating, tree-like behavior. Here, we tested this prediction using three methods designed to investigate gene flow in autosome-wide genotype data from 3,528 unrelated individuals from 163 global samples. Specifically, we investigated whether Cushitic ancestry has an East African or Middle Eastern origin. We found evidence for non-tree-like behavior in the form of four migration events. First, we found that Cushitic ancestry is a mixture of ancestries closely related to Arabian ancestry and Nilo-Saharan or Omotic ancestry. We found evidence for additional migration events in the histories of: 1) Indian and Arabian ancestries, 2) Kalash ancestry, and 3) Native American and Northern European ancestries. These findings, based on analysis of ancestry of present-day humans, reveal migration in the distant past and provide new insights into human history.

Y-chromosome phylogeographic analysis of the Greek-Cypriot population reveals elements consistent with Neolithic and Bronze Age settlements

Background

The archeological record indicates that the permanent settlement of Cyprus began with pioneering agriculturalists circa 11,000 years before present, (ca. 11,000 y BP). Subsequent colonization events followed, some recognized regionally. Here, we assess the Y-chromosome structure of Cyprus in context to regional populations and correlate it to phases of prehistoric colonization.

Results

Analysis of haplotypes from 574 samples showed that island-wide substructure was barely significant in a spatial analysis of molecular variance (SAMOVA). However, analyses of molecular variance (AMOVA) of haplogroups using 92 binary markers genotyped in 629 Cypriots revealed that the proportion of variance among the districts was irregularly distributed. Principal component analysis (PCA) revealed potential genetic associations of Greek-Cypriots with neighbor populations. Contrasting haplogroups in the PCA were used as surrogates of parental populations. Admixture analyses suggested that the majority of G2a-P15 and R1b-M269 components were contributed by Anatolia and Levant sources, respectively, while Greece Balkans supplied the majority of E-V13 and J2a-M67. Haplotype-based expansion times were at historical levels suggestive of recent demography.

Conclusions

Analyses of Cypriot haplogroup data are consistent with two stages of prehistoric settlement. E-V13 and E-M34 are widespread, and PCA suggests sourcing them to the Balkans and Levant/Anatolia, respectively. The persistent pre-Greek component is represented by elements of G2-U5(xL30) haplogroups: U5*, PF3147, and L293. J2b-M205 may contribute also to the pre-Greek strata. The majority of R1b-Z2105 lineages occur in both the westernmost and easternmost districts. Distinctively, sub-haplogroup R1b- M589 occurs only in the east. The absence of R1b- M589 lineages in Crete and the Balkans and the presence in Asia Minor are compatible with Late Bronze Age influences from Anatolia rather than from Mycenaean Greeks.

Electronic supplementary material

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

Exploring the Y Chromosomal Ancestry of Modern Panamanians

Geologically, Panama belongs to the Central American land-bridge between North and South America crossed by Homo sapiens >14 ka ago. Archaeologically, it belongs to a wider Isthmo-Colombian Area. Today, seven indigenous ethnic groups account for 12.3% of Panama’s population. Five speak Chibchan languages and are characterized by low genetic diversity and a high level of differentiation. In addition, no evidence of differential structuring between maternally and paternally inherited genes has been reported in isthmian Chibchan cultural groups. Recent data have shown that 83% of the Panamanian general population harbour mitochondrial DNAs (mtDNAs) of Native American ancestry. Considering differential male/female mortality at European contact and multiple degrees of geographical and genetic isolation over the subsequent five centuries, the Y-chromosome Native American component is expected to vary across different geographic regions and communities in Panama. To address this issue, we investigated Y-chromosome variation in 408 modern males from the nine provinces of Panama and one indigenous territory (the comarca of Kuna Yala). In contrast to mtDNA data, the Y-chromosome Native American component (haplogroup Q) exceeds 50% only in three populations facing the Caribbean Sea: the comarca of Kuna Yala and Bocas del Toro province where Chibchan languages are spoken by the majority, and the province of Colón where many Kuna and people of mixed indigenous-African-and-European descent live. Elsewhere the Old World component is dominant and mostly represented by western Eurasian haplogroups, which signal the strong male genetic impact of invaders. Sub-Saharan African input accounts for 5.9% of male haplotypes. This reflects the consequences of the colonial Atlantic slave trade and more recent influxes of West Indians of African heritage. Overall, our findings reveal a local evolution of the male Native American ancestral gene pool, and a strong but geographically differentiated unidirectional sex bias in the formation of local modern Panamanian populations.

The maternal perspective for five Slovenian regions: The importance of regional sampling

BACKGROUND

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

AIM

The aims of this study were to characterize the Slovenian mitochondrial gene pool at the micro-geographic level and to compare it with surrounding populations.

SUBJECTS AND METHODS

A total of 402 individuals from five Slovenian regions were analysed in this study by typing HVR I, HVR II and coding region polymorphisms of mtDNA.

RESULTS

Analysis revealed 47 haplogroups and sub-haplogroups, the most common of which were H*, H1, J1c, T2 and U5a. Intra-population comparisons revealed a sharp gradient of the J1c haplogroup between Slovenian regions, with a peak frequency of 24.5% being observed in the population of the Littoral Region.

CONCLUSION

The sharp gradient of the J1c haplogroup between Slovenian regions is in line with the archaeological horizon known as Impressed Ware culture and could, therefore, represent a genetic trace of the early Neolithic expansion route along the East Adriatic coastal region.

Letter to the editor: Genetics and the archaeology of ancient Israel

This letter is a call for DNA testing on ancient skeletal materials from the southern Levant to begin a database of genetic information of the inhabitants of this crossroads region. In this region, during the Iron I period traditionally dated to circa 1200-1000 BCE, archaeologists and biblical historians view the earliest presence of a group that called itself Israel. They lived in villages in the varied hill countries of the region, contemporary with urban settlements in the coastal plains, inland valleys, and central hill country attributed to varied indigenous groups collectively called Canaanite. The remnants of Egyptian imperial presence in the region lasted until around 1150 BCE, postdating the arrival of an immigrant group from the Aegean called the Philistines circa 1175 BCE. The period that follows in the southern Levant is marked by the development of territorial states throughout the region, circa 1000-800 BCE. These patrimonial kingdoms, including the United Kingdom of Israel and the divided kingdoms of northern Israel and Judah, coalesced varied peoples under central leadership and newly founded administrative and religious bureaucracies. Ancient DNA testing will give us a further refined understanding of the individuals who peopled the region of the southern Levant throughout its varied archaeological and historic periods and provide scientific data that will support, refute, or nuance our sociohistoric reconstruction of ancient group identities. These social identities may or may not map onto genetic data, but without sampling of ancient DNA we may never know. A database of ancient DNA will also allow for comparisons with modern DNA samples collected throughout the greater region and the Mediterranean littoral, giving a more robust understanding of the long historical trajectories of regional human genetics and the genetics of varied ancestral groups of today's Jewish populations and other cultural groups in the modern Middle East and Mediterranean.

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.

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.