The Y-chromosome tree bursts into leaf: 13,000 high-confidence SNPs covering the majority of known clades

Y Chromosome Sequences Reveal a Short Beringian Standstill, Rapid Expansion, and early Population structure of Native American Founders

The Americas were the last inhabitable continents to be occupied by humans, with a growing multidisciplinary consensus for entry 15-25 thousand years ago (kya) from northeast Asia via the former Beringia land bridge [1-4]. Autosomal DNA analyses have dated the separation of Native American ancestors from the Asian gene pool to 23 kya or later [5, 6] and mtDNA analyses to ∼25 kya [7], followed by isolation ("Beringian Standstill" [8, 9]) for 2.4-9 ky and then a rapid expansion throughout the Americas. Here, we present a calibrated sequence-based analysis of 222 Native American and relevant Eurasian Y chromosomes (24 new) from haplogroups Q and C [10], with four major conclusions. First, we identify three to four independent lineages as autochthonous and likely founders: the major Q-M3 and rarer Q-CTS1780 present throughout the Americas, the very rare C3-MPB373 in South America, and possibly the C3-P39/Z30536 in North America. Second, from the divergence times and Eurasian/American distribution of lineages, we estimate a Beringian Standstill duration of 2.7 ky or 4.6 ky, according to alternative models, and entry south of the ice sheet after 19.5 kya. Third, we describe the star-like expansion of Q-M848 (within Q-M3) starting at 15 kya [11] in the Americas, followed by establishment of substantial spatial structure in South America by 12 kya. Fourth, the deep branches of the Q-CTS1780 lineage present at low frequencies throughout the Americas today [12] may reflect a separate out-of-Beringia dispersal after the melting of the glaciers at the end of the Pleistocene.

The role of matrilineality in shaping patterns of Y chromosome and mtDNA sequence variation in southwestern Angola

Southwestern Angola is a region characterized by contact between indigenous foragers and incoming food-producers, involving genetic and cultural exchanges between peoples speaking Kx'a, Khoe-Kwadi, and Bantu languages. Although present-day Bantu speakers share a patrilocal residence pattern and matrilineal principle of clan and group membership, a highly stratified social setting divides dominant pastoralists from marginalized groups that subsist on alternative strategies and have previously been thought to have pre-Bantu origins. Here, we compare new high-resolution sequence data from 2.3 Mb of the male-specific region of the Y chromosome (MSY) from 170 individuals with previously reported mitochondrial DNA (mtDNA) genomes, to investigate the population history of seven representative southwestern Angolan groups (Himba, Kuvale, Kwisi, Kwepe, Twa, Tjimba, !Xun), and to study the causes and consequences of sex-biased processes in their genetic variation. We found no clear link between the formerly Kwadi-speaking Kwepe and pre-Bantu eastern African migrants, and no pre-Bantu MSY lineages among Bantu-speaking groups, except for small amounts of "Khoisan" introgression. We therefore propose that irrespective of their subsistence strategies, all Bantu-speaking groups of the area share a male Bantu origin. Additionally, we show that in Bantu-speaking groups, the levels of among-group and between-group variation are higher for mtDNA than for MSY. These results, together with our previous demonstration that the matriclanic systems of southwestern Angolan Bantu groups are genealogically consistent, suggest that matrilineality strongly enhances both female population sizes and interpopulation mtDNA variation.

Genetic diversities and phylogenetic analyses of three Chinese main ethnic groups in southwest China: A Y-Chromosomal STR study

Short tandem repeats (STRs) located on the Y chromosome with the properties of male-specific inheritance and haploidy are widely used in forensics to analyze paternal genealogies and match male trace donors to evidence. Besides, Y-chromosomal haplotypes play an important role in providing breathtaking insights into population genetic history. However, the genetic diversity and forensic characteristics of Y-STRs in Guizhou main ethnic groups (Hans, Miaos and Bouyeis) remain uncharacterized. Here, we obtained Y-chromosomal 23-marker haplotypes in three Guizhou populations and submitted the first batch of Y-STR haplotype data to the YHRD. The HD in the aforementioned three populations are 0.99990, 0.99983, and 0.99979, respectively, and DC values are 0.9902, 0.9908, and 0.97959, respectively. Subsequently, genetic differentiation between our newly studied populations and reference groups along ethnic/administrative divisions, as well as national/continental boundaries were investigated via AMOVA, MDS, and phylogenetic relationship reconstruction. Significant genetic differentiations from our subjects and other groups are identified in ethnically, linguistically and geographically diverse populations, including most prominently Tibetans and Uyghurs among 30 mainland Chinese populations, Taiwanese groups and others among 58 Asian populations, as well as African groups and others among 89 worldwide populations. Qiannan Bouyei has a close genetic relationship with Guangxi Zhuang, and Zunyi Han and Qiandongnan Miao have close genetic affinity with Hunan Han and Guizhou Shui, respectively. Collectively, this new-generation Y-STR amplification system can be used as a supplementary tool in forensic identification and male parentage testing and even pedigree search.

Recombination hotspots in an extended human pseudoautosomal domain predicted from double-strand break maps and characterized by sperm-based crossover analysis

The human X and Y chromosomes are heteromorphic but share a region of homology at the tips of their short arms, pseudoautosomal region 1 (PAR1), that supports obligate crossover in male meiosis. Although the boundary between pseudoautosomal and sex-specific DNA has traditionally been regarded as conserved among primates, it was recently discovered that the boundary position varies among human males, due to a translocation of ~110 kb from the X to the Y chromosome that creates an extended PAR1 (ePAR). This event has occurred at least twice in human evolution. So far, only limited evidence has been presented to suggest this extension is recombinationally active. Here, we sought direct proof by examining thousands of gametes from each of two ePAR-carrying men, for two subregions chosen on the basis of previously published male X-chromosomal meiotic double-strand break (DSB) maps. Crossover activity comparable to that seen at autosomal hotspots was observed between the X and the ePAR borne on the Y chromosome both at a distal and a proximal site within the 110-kb extension. Other hallmarks of classic recombination hotspots included evidence of transmission distortion and GC-biased gene conversion. We observed good correspondence between the male DSB clusters and historical recombination activity of this region in the X chromosomes of females, as ascertained from linkage disequilibrium analysis; this suggests that this region is similarly primed for crossover in both male and female germlines, although sex-specific differences may also exist. Extensive resequencing and inference of ePAR haplotypes, placed in the framework of the Y phylogeny as ascertained by both Y microsatellites and single nucleotide polymorphisms, allowed us to estimate a minimum rate of crossover over the entire ePAR region of 6-fold greater than genome average, comparable with pedigree estimates of PAR1 activity generally. We conclude ePAR very likely contributes to the critical crossover function of PAR1.

95% of our genome is contained in 22 pairs of chromosomes shared by all humans. However, women and men differ in their sex chromosomes: while women have two X chromosomes, men have an X and a smaller, sex-determining Y chromosome. To ensure correct partition of X and Y into sperm, genetic exchange (crossover) must occur between these very different chromosomes in a short, shared region. The location of the boundary of this region was thought to have been conserved since before the divergence from old world monkeys at least 27 million years ago, but recently it has been shown that some human males carry an extended version on their Y chromosomes, thanks to the transposition of a piece of DNA from the X chromosome. Here, we asked if genetic exchange occurs in this newly extended region. To do this, we used previously published information that signposted the positions within the X chromosome segment which exhibit the hallmarks of crossover initiation. We then sought direct evidence of crossover in the sperm of men carrying the extension. This work showed that the signposts were accurate, pointing to frequent crossover in this novel shared sex-chromosomal domain.

Genetic Diversity on the Sex Chromosomes

Levels and patterns of genetic diversity can provide insights into a population’s history. In species with sex chromosomes, differences between genomic regions with unique inheritance patterns can be used to distinguish between different sets of possible demographic and selective events. This review introduces the differences in population history for sex chromosomes and autosomes, provides the expectations for genetic diversity across the genome under different evolutionary scenarios, and gives an introductory description for how deviations in these expectations are calculated and can be interpreted. Predominantly, diversity on the sex chromosomes has been used to explore and address three research areas: 1) Mating patterns and sex-biased variance in reproductive success, 2) signatures of selection, and 3) evidence for modes of speciation and introgression. After introducing the theory, this review catalogs recent studies of genetic diversity on the sex chromosomes across species within the major research areas that sex chromosomes are typically applied to, arguing that there are broad similarities not only between male-heterogametic (XX/XY) and female-heterogametic (ZZ/ZW) sex determination systems but also any mating system with reduced recombination in a sex-determining region. Further, general patterns of reduced diversity in nonrecombining regions are shared across plants and animals. There are unique patterns across populations with vastly different patterns of mating and speciation, but these do not tend to cluster by taxa or sex determination system.

Understanding 6th-century barbarian social organization and migration through paleogenomics

Despite centuries of research, much about the barbarian migrations that took place between the fourth and sixth centuries in Europe remains hotly debated. To better understand this key era that marks the dawn of modern European societies, we obtained ancient genomic DNA from 63 samples from two cemeteries (from Hungary and Northern Italy) that have been previously associated with the Longobards, a barbarian people that ruled large parts of Italy for over 200 years after invading from Pannonia in 568 CE. Our dense cemetery-based sampling revealed that each cemetery was primarily organized around one large pedigree, suggesting that biological relationships played an important role in these early medieval societies. Moreover, we identified genetic structure in each cemetery involving at least two groups with different ancestry that were very distinct in terms of their funerary customs. Finally, our data are consistent with the proposed long-distance migration from Pannonia to Northern Italy.

The Longobards invaded and conquered much of Italy after the fall of the Western Roman Empire. Here, the authors sequence and analyze ancient genomic DNA from 63 samples from two cemeteries associated with the Longobards and identify kinship networks and two distinct genetic and cultural groups in each.

Mediterranean Y-chromosome 2.0-why the Y in the Mediterranean is still relevant in the postgenomic era

CONTEXT

Due to its unique paternal inheritance, the Y-chromosome has been a highly popular marker among population geneticists for over two decades. Recently, the advent of cost-effective genome-wide methods has unlocked information-rich autosomal genomic data, paving the way to the postgenomic era. This seems to have announced the decreasing popularity of investigating Y-chromosome variation, which provides only the paternal perspective of human ancestries and is strongly influenced by genetic drift and social behaviour.

OBJECTIVE

For this special issue on population genetics of the Mediterranean, the aim was to demonstrate that the Y-chromosome still provides important insights in the postgenomic era and in a time when ancient genomes are becoming exponentially available.

METHODS

A systematic literature search on Y-chromosomal studies in the Mediterranean was performed.

RESULTS

Several applications of Y-chromosomal analysis with future opportunities are formulated and illustrated with studies on Mediterranean populations.

CONCLUSIONS

There will be no reduced interest in Y-chromosomal studies going from reconstruction of male-specific demographic events to ancient DNA applications, surname history and population-wide estimations of extra-pair paternity rates. Moreover, more initiatives are required to collect population genetic data of Y-chromosomal markers for forensic research, and to include Y-chromosomal data in GWAS investigations and studies on male infertility.

Paternal origin of Paleo-Indians in Siberia: insights from Y-chromosome sequences

The expansion of modern humans to the American continent after the Last Glacial Maximum led the way to the present-day distribution of American aborigines. Recent advances in autosomal DNA research and expanded testing of mtDNA lineages has provided a clearer picture of the number and timing of founding lineages. However, both autosomal DNA and mtDNA research have provided unresolved competing theories between the short-term and the long-term models of the Beringian standstill hypothesis. Further, the source of founding paternal lineages of American aborigines and their relationship with ancient Siberia populations remains ambiguous. In this study, we reanalyzed a 7.0 Mbp region of 132 paternal Y-chromosome sequences, including 39 newly reported ones, of male samples from American aborigines and Eurasian populations. Among Eurasian samples, we identified Y-chromosome branches that are most closely related to known American aborigine founding lineages, that is, Q1-L804 links to Q1-M3, Q1-L330 links to Q1-Z780, Q1-M120 links to Q1-B143, and C2-F1756 links to C2-P39. The revised phylogenetic tree and age estimates indicate a narrow timeframe (~15.3-14.3 kya) for the upper time limit of human entry to the American continent. Our analysis suggests that the in situ differentiation of Q-M242 in Central Eurasia and South Siberia region gave rise to numerous sub-lineages older than 15.3 kya, and the founding of Paleo-Indian paternal lineages is part of the great Q1-L53 diffusion throughout the Eurasia after the Last Glacial Maximum. The results of our study will assist in future studies of the history of modern populations in Eurasia and the Americas.

Carriers of mitochondrial DNA macrohaplogroup L3 basal lineages migrated back to Africa from Asia around 70,000 years ago

Background

The main unequivocal conclusion after three decades of phylogeographic mtDNA studies is the African origin of all extant modern humans. In addition, a southern coastal route has been argued for to explain the Eurasian colonization of these African pioneers. Based on the age of macrohaplogroup L3, from which all maternal Eurasian and the majority of African lineages originated, the out-of-Africa event has been dated around 60-70 kya. On the opposite side, we have proposed a northern route through Central Asia across the Levant for that expansion and, consistent with the fossil record, we have dated it around 125 kya. To help bridge differences between the molecular and fossil record ages, in this article we assess the possibility that mtDNA macrohaplogroup L3 matured in Eurasia and returned to Africa as basal L3 lineages around 70 kya.

Results

The coalescence ages of all Eurasian (M,N) and African (L3 ) lineages, both around 71 kya, are not significantly different. The oldest M and N Eurasian clades are found in southeastern Asia instead near of Africa as expected by the southern route hypothesis. The split of the Y-chromosome composite DE haplogroup is very similar to the age of mtDNA L3. An Eurasian origin and back migration to Africa has been proposed for the African Y-chromosome haplogroup E. Inside Africa, frequency distributions of maternal L3 and paternal E lineages are positively correlated. This correlation is not fully explained by geographic or ethnic affinities. This correlation rather seems to be the result of a joint and global replacement of the old autochthonous male and female African lineages by the new Eurasian incomers.

Conclusions

These results are congruent with a model proposing an out-of-Africa migration into Asia, following a northern route, of early anatomically modern humans carrying pre-L3 mtDNA lineages around 125 kya, subsequent diversification of pre-L3 into the basal lineages of L3, a return to Africa of Eurasian fully modern humans around 70 kya carrying the basal L3 lineages and the subsequent diversification of Eurasian-remaining L3 lineages into the M and N lineages in the outside-of-Africa context, and a second Eurasian global expansion by 60 kya, most probably, out of southeast Asia. Climatic conditions and the presence of Neanderthals and other hominins might have played significant roles in these human movements. Moreover, recent studies based on ancient DNA and whole-genome sequencing are also compatible with this hypothesis.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1211-4) contains supplementary material, which is available to authorized users.

A finely resolved phylogeny of Y chromosome Hg J illuminates the processes of Phoenician and Greek colonizations in the Mediterranean

In order to improve the phylogeography of the male-specific genetic traces of Greek and Phoenician colonizations on the Northern coasts of the Mediterranean, we performed a geographically structured sampling of seven subclades of haplogroup J in Turkey, Greece and Italy. We resequenced 4.4 Mb of Y-chromosome in 58 subjects, obtaining 1079 high quality variants. We did not find a preferential coalescence of Turkish samples to ancestral nodes, contradicting the simplistic idea of a dispersal and radiation of Hg J as a whole from the Middle East. Upon calibration with an ancient Hg J chromosome, we confirmed that signs of Holocenic Hg J radiations are subtle and date mainly to the Bronze Age. We pinpointed seven variants which could potentially unveil star clusters of sequences, indicative of local expansions. By directly genotyping these variants in Hg J carriers and complementing with published resequenced chromosomes (893 subjects), we provide strong temporal and distributional evidence for markers of the Greek settlement of Magna Graecia (J2a-L397) and Phoenician migrations (rs760148062). Our work generated a minimal but robust list of evolutionarily stable markers to elucidate the demographic dynamics and spatial domains of male-mediated movements across and around the Mediterranean, in the last 6,000 years.

A phylogenetic framework facilitates Y-STR variant discovery and classification via massively parallel sequencing

•

23 Y-chromosomal STRs (PPY23) reanalysed by massively parallel sequencing.

•

Phylogeny-based approach captures wide range of sequence variants in 100 samples.

•

STR variants described in phase with their flanking sequences.

•

Phylogenetic framework clarifies allele nomenclature and mutation processes.

Short tandem repeats on the male-specific region of the Y chromosome (Y-STRs) are permanently linked as haplotypes, and therefore Y-STR sequence diversity can be considered within the robust framework of a phylogeny of haplogroups defined by single nucleotide polymorphisms (SNPs). Here we use massively parallel sequencing (MPS) to analyse the 23 Y-STRs in Promega’s prototype PowerSeq™ Auto/Mito/Y System kit (containing the markers of the PowerPlex® Y23 [PPY23] System) in a set of 100 diverse Y chromosomes whose phylogenetic relationships are known from previous megabase-scale resequencing. Including allele duplications and alleles resulting from likely somatic mutation, we characterised 2311 alleles, demonstrating 99.83% concordance with capillary electrophoresis (CE) data on the same sample set. The set contains 267 distinct sequence-based alleles (an increase of 58% compared to the 169 detectable by CE), including 60 novel Y-STR variants phased with their flanking sequences which have not been reported previously to our knowledge. Variation includes 46 distinct alleles containing non-reference variants of SNPs/indels in both repeat and flanking regions, and 145 distinct alleles containing repeat pattern variants (RPV). For DYS385a,b, DYS481 and DYS390 we observed repeat count variation in short flanking segments previously considered invariable, and suggest new MPS-based structural designations based on these. We considered the observed variation in the context of the Y phylogeny: several specific haplogroup associations were observed for SNPs and indels, reflecting the low mutation rates of such variant types; however, RPVs showed less phylogenetic coherence and more recurrence, reflecting their relatively high mutation rates. In conclusion, our study reveals considerable additional diversity at the Y-STRs of the PPY23 set via MPS analysis, demonstrates high concordance with CE data, facilitates nomenclature standardisation, and places Y-STR sequence variants in their phylogenetic context.

The peopling of the last Green Sahara revealed by high-coverage resequencing of trans-Saharan patrilineages

BACKGROUND

Little is known about the peopling of the Sahara during the Holocene climatic optimum, when the desert was replaced by a fertile environment.

RESULTS

In order to investigate the role of the last Green Sahara in the peopling of Africa, we deep-sequence the whole non-repetitive portion of the Y chromosome in 104 males selected as representative of haplogroups which are currently found to the north and to the south of the Sahara. We identify 5,966 mutations, from which we extract 142 informative markers then genotyped in about 8,000 subjects from 145 African, Eurasian and African American populations. We find that the coalescence age of the trans-Saharan haplogroups dates back to the last Green Sahara, while most northern African or sub-Saharan clades expanded locally in the subsequent arid phase.

CONCLUSIONS

Our findings suggest that the Green Sahara promoted human movements and demographic expansions, possibly linked to the adoption of pastoralism. Comparing our results with previously reported genome-wide data, we also find evidence for a sex-biased sub-Saharan contribution to northern Africans, suggesting that historical events such as the trans-Saharan slave trade mainly contributed to the mtDNA and autosomal gene pool, whereas the northern African paternal gene pool was mainly shaped by more ancient events.

Extensive geographical and social structure in the paternal lineages of Saudi Arabia revealed by analysis of 27 Y-STRs

Saudi Arabia's indigenous population is organized into patrilineal descent groups, but to date, little has been done to characterize its population structure, in particular with respect to the male-specific region of the Y chromosome. We have used the 27-STR Yfiler^® Plus kit to generate haplotypes in 597 unrelated Saudi males, classified into five geographical regions (North, South, Central, East and West). Overall, Yfiler^® Plus provides a good discrimination capacity of 95.3%, but this is greatly reduced (74.7%) when considering the reduced Yfiler^® set of 17 Y-STRs, justifying the use of the expanded set of markers in this population. Comparison of the five geographical divisions reveals striking differences, with low diversity and similar haplotype spectra in the Central and Northern regions, and high diversity and similar haplotype spectra in the East and West. These patterns likely reflect the geographical isolation of the desert heartland of the peninsula, and the proximity to the sea of the Eastern and Western areas, and consequent historical immigration. We predicted haplogroups from Y-STR haplotypes, testing the performance of prediction by using a large independent set of Saudi Arabian Y-STR + Y-SNP data. Prediction indicated predominance (71%) of haplogroup J1, which was significantly more common in Central, Northern and Southern groups than in East and West, and formed a star-like expansion cluster in a median-joining network with an estimated age of ∼2800 years. Most of our 597 participants were sampled within Saudi Arabia itself, but ∼16% were sampled in the UK. Despite matching these two groups by home sub-region, we observed significant differences in haplotype and predicted haplogroup constitutions overall, and for most sub-regions individually. This suggests social structure influencing the probability of leaving Saudi Arabia, correlated with different Y-chromosome compositions. The UK-recruited sample is an inappropriate proxy for Saudi Arabia generally, and caution is needed when considering expatriate groups as representative of country of origin. Our study shows the importance of geographical and social structuring that may affect the utility of forensic databases and the interpretation of Y-STR profiles.

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.

Defining Y-SNP variation among the Flemish population (Western Europe) by full genome sequencing

Y-chromosomal single nucleotide polymorphisms (Y-SNPs) represent a powerful tool in forensic research and casework, especially for inferring paternal ancestry of unknown perpetrators and unidentified bodies. However, the wealth of recently discovered Y-SNPs, the 'jungle' of different evolutionary lineage trees and nomenclatures, and the lack of population-wide data of many phylogenetically mapped Y-SNPs, limits the use of Y-SNPs in routine forensic approaches. Recently, a concise reference phylogeny of the human Y chromosome, the 'Minimal Reference Y-tree', was introduced aiming to provide a stable phylogeny with optimal global discrimination capacity by including the most resolving Y-SNPs. Here, we obtained a representative sample of 270 whole-genome sequences (WGS) to grasp the Y-SNP variation within the autochthonous Flemish population (Belgium, Western Europe) according to this reference Y-tree. The high quality of the Y-SNP calling was guaranteed for the WGS sample as well as its representativeness for the Flemish population based on the comparison of the main haplogroup frequencies with those from earlier studies on Flanders and the Netherlands. The 270 Flemish Y chromosomes were assigned to 98 different sub-haplogroups of the Minimal Reference Y-tree, showing its high potential of discrimination and confirming the spectrum of evolutionary lineages within Western Europe in general and within Flanders in particular. The full database with all Y-SNP calls of the Flemish sample is public available for future updates including forensic and population genetic studies. New initiatives to categorise Y-SNP variation in other populations according to the reference phylogeny of the Y chromosome are highly encouraged for forensic applications. Recommendations to realise such future population sample sets are discussed based on this study.

Genetic diversity in Svaneti and its implications for the human settlement of the Highland Caucasus

OBJECTIVES

In this study, we characterized genetic diversity in the Svans from northwestern Georgia to better understand the phylogeography of their genetic lineages, determine whether genetic diversity in the highland South Caucasus has been shaped by language or geography, and assess whether Svan genetic diversity was structured by regional residence patterns.

MATERIALS AND METHODS

We analyzed mtDNA and Y-chromosome variation in 184 individuals from 13 village districts and townlets located throughout the region. For all individuals, we analyzed mtDNA diversity through control region sequencing, and, for males, we analyzed Y-chromosome diversity through SNP and STR genotyping. The resulting data were compared with those for populations from the Caucasus and Middle East.

RESULTS

We observed significant mtDNA heterogeneity in Svans, with haplogroups U1-U7, H, K, and W6 being common there. By contrast, ∼78% of Svan males belonged to haplogroup G2a, with the remainder falling into four other haplogroups (J2a1, I2, N, and R1a). While showing a distinct genetic profile, Svans also clustered with Caucasus populations speaking languages from different families, suggesting a deep common ancestry for all of them. The mtDNA data were not structured by geography or linguistic affiliation, whereas the NRY data were influenced only by geography.

DISCUSSION

These patterns of genetic variation confirm a complex set of geographic sources and settlement phases for the Caucasus highlands. Such patterns may also reflect social and cultural practices in the region. The high frequency and antiquity of Y-chromosome haplogroup G2a in this region further points to its emergence there.

Population resequencing of European mitochondrial genomes highlights sex-bias in Bronze Age demographic expansions

Interpretations of genetic data concerning the prehistory of Europe have long been a subject of great debate, but increasing amounts of ancient and modern DNA data are now providing new and more informative evidence. Y-chromosome resequencing studies in Europe have highlighted the prevalence of recent expansions of male lineages, and focused interest on the Bronze Age as a period of cultural and demographic change. These findings contrast with phylogeographic studies based on mitochondrial DNA (mtDNA), which have been interpreted as supporting expansions from glacial refugia. Here we have undertaken a population-based resequencing of complete mitochondrial genomes in Europe and the Middle East, in 340 samples from 17 populations for which Y-chromosome sequence data are also available. Demographic reconstructions show no signal of Bronze Age expansion, but evidence of Paleolithic expansions in all populations except the Saami, and with an absence of detectable geographical pattern. In agreement with previous inference from modern and ancient DNA data, the unbiased comparison between the mtDNA and Y-chromosome population datasets emphasizes the sex-biased nature of recent demographic transitions in Europe.

Next Generation Sequencing Plus (NGS+) with Y-chromosomal Markers for Forensic Pedigree Searches

There is high demand for forensic pedigree searches with Y-chromosome short tandem repeat (Y-STR) profiling in large-scale crime investigations. However, when two Y-STR haplotypes have a few mismatched loci, it is difficult to determine if they are from the same male lineage because of the high mutation rate of Y-STRs. Here we design a new strategy to handle cases in which none of pedigree samples shares identical Y-STR haplotype. We combine next generation sequencing (NGS), capillary electrophoresis and pyrosequencing under the term ‘NGS+’ for typing Y-STRs and Y-chromosomal single nucleotide polymorphisms (Y-SNPs). The high-resolution Y-SNP haplogroup and Y-STR haplotype can be obtained with NGS+. We further developed a new data-driven decision rule, FSindex, for estimating the likelihood for each retrieved pedigree. Our approach enables positive identification of pedigree from mismatched Y-STR haplotypes. It is envisaged that NGS+ will revolutionize forensic pedigree searches, especially when the person of interest was not recorded in forensic DNA database.

The Y chromosome: a blueprint for men's health?

The Y chromosome has long been considered a 'genetic wasteland' on a trajectory to completely disappear from the human genome. The perception of its physiological function was restricted to sex determination and spermatogenesis. These views have been challenged in recent times with the identification of multiple ubiquitously expressed Y-chromosome genes and the discovery of several unexpected associations between the Y chromosome, immune system and complex polygenic traits. The collected evidence suggests that the Y chromosome influences immune and inflammatory responses in men, translating into genetically programmed susceptibility to diseases with a strong immune component. Phylogenetic studies reveal that carriers of a common European lineage of the Y chromosome (haplogroup I) possess increased risk of coronary artery disease. This occurs amidst upregulation of inflammation and suppression of adaptive immunity in this Y lineage, as well as inferior outcomes in human immunodeficiency virus infection. From structural analysis and experimental data, the UTY (Ubiquitously Transcribed Tetratricopeptide Repeat Containing, Y-Linked) gene is emerging as a promising candidate underlying the associations between Y-chromosome variants and the immunity-driven susceptibility to complex disease. This review synthesises the recent structural, experimental and clinical insights into the human Y chromosome in the context of men's susceptibility to disease (with a particular emphasis on cardiovascular disease) and provides an overview of the paradigm shift in the perception of the Y chromosome.

Y Chromosome Uncovers the Recent Oriental Origin of Modern Stallions

The Y chromosome directly reflects male genealogies, but the extremely low Y chromosome sequence diversity in horses has prevented the reconstruction of stallion genealogies [1, 2]. Here, we resolve the first Y chromosome genealogy of modern horses by screening 1.46 Mb of the male-specific region of the Y chromosome (MSY) in 52 horses from 21 breeds. Based on highly accurate pedigree data, we estimated the de novo mutation rate of the horse MSY and showed that various modern horse Y chromosome lineages split much later than the domestication of the species. Apart from few private northern European haplotypes, all modern horse breeds clustered together in a roughly 700-year-old haplogroup that was transmitted to Europe by the import of Oriental stallions. The Oriental horse group consisted of two major subclades: the Original Arabian lineage and the Turkoman horse lineage. We show that the English Thoroughbred MSY was derived from the Turkoman lineage and that English Thoroughbred sires are largely responsible for the predominance of this haplotype in modern horses.

The Connection of the Genetic, Cultural and Geographic Landscapes of Transoxiana

We have analyzed Y-chromosomal variation in populations from Transoxiana, a historical region covering the southwestern part of Central Asia. We studied 780 samples from 10 regional populations of Kazakhs, Uzbeks, Turkmens, Dungans, and Karakalpaks using 35 SNP and 17 STR markers. Analysis of haplogroup frequencies using multidimensional scaling and principal component plots, supported by an analysis of molecular variance, showed that the geographic landscape of Transoxiana, despite its distinctiveness and diversity (deserts, fertile river basins, foothills and plains) had no strong influence on the genetic landscape. The main factor structuring the gene pool was the mode of subsistence: settled agriculture or nomadic pastoralism. Investigation of STR-based clusters of haplotypes and their ages revealed that cultural and demic expansions of Transoxiana were not closely connected with each other. The Arab cultural expansion introduced Islam to the region but did not leave a significant mark on the pool of paternal lineages. The Mongol expansion, in contrast, had enormous demic success, but did not impact cultural elements like language and religion. The genealogy of Muslim missionaries within the settled agricultural communities of Transoxiana was based on spiritual succession passed from teacher to disciple. However, among Transoxianan nomads, spiritual and biological succession became merged.

Human Y-chromosome variation in the genome-sequencing era

The properties of the human Y chromosome - namely, male specificity, haploidy and escape from crossing over - make it an unusual component of the genome, and have led to its genetic variation becoming a key part of studies of human evolution, population history, genealogy, forensics and male medical genetics. Next-generation sequencing (NGS) technologies have driven recent progress in these areas. In particular, NGS has yielded direct estimates of mutation rates, and an unbiased and calibrated molecular phylogeny that has unprecedented detail. Moreover, the availability of direct-to-consumer NGS services is fuelling a rise of 'citizen scientists', whose interest in resequencing their own Y chromosomes is generating a wealth of new data.

Carriers of mitochondrial DNA macrohaplogroup R colonized Eurasia and Australasia from a southeast Asia core area

BACKGROUND

The colonization of Eurasia and Australasia by African modern humans has been explained, nearly unanimously, as the result of a quick southern coastal dispersal route through the Arabian Peninsula, the Indian subcontinent, and the Indochinese Peninsula, to reach Australia around 50 kya. The phylogeny and phylogeography of the major mitochondrial DNA Eurasian haplogroups M and N have played the main role in giving molecular genetics support to that scenario. However, using the same molecular tools, a northern route across central Asia has been invoked as an alternative that is more conciliatory with the fossil record of East Asia. Here, we assess as the Eurasian macrohaplogroup R fits in the northern path.

RESULTS

Haplogroup U, with a founder age around 50 kya, is one of the oldest clades of macrohaplogroup R in western Asia. The main branches of U expanded in successive waves across West, Central and South Asia before the Last Glacial Maximum. All these dispersions had rather overlapping ranges. Some of them, as those of U6 and U3, reached North Africa. At the other end of Asia, in Wallacea, another branch of macrohaplogroup R, haplogroup P, also independently expanded in the area around 52 kya, in this case as isolated bursts geographically well structured, with autochthonous branches in Australia, New Guinea, and the Philippines.

CONCLUSIONS

Coeval independently dispersals around 50 kya of the West Asia haplogroup U and the Wallacea haplogroup P, points to a halfway core area in southeast Asia as the most probable centre of expansion of macrohaplogroup R, what fits in the phylogeographic pattern of its ancestor, macrohaplogroup N, for which a northern route and a southeast Asian origin has been already proposed.

The study of human Y chromosome variation through ancient DNA

High throughput sequencing methods have completely transformed the study of human Y chromosome variation by offering a genome-scale view on genetic variation retrieved from ancient human remains in context of a growing number of high coverage whole Y chromosome sequence data from living populations from across the world. The ancient Y chromosome sequences are providing us the first exciting glimpses into the past variation of male-specific compartment of the genome and the opportunity to evaluate models based on previously made inferences from patterns of genetic variation in living populations. Analyses of the ancient Y chromosome sequences are challenging not only because of issues generally related to ancient DNA work, such as DNA damage-induced mutations and low content of endogenous DNA in most human remains, but also because of specific properties of the Y chromosome, such as its highly repetitive nature and high homology with the X chromosome. Shotgun sequencing of uniquely mapping regions of the Y chromosomes to sufficiently high coverage is still challenging and costly in poorly preserved samples. To increase the coverage of specific target SNPs capture-based methods have been developed and used in recent years to generate Y chromosome sequence data from hundreds of prehistoric skeletal remains. Besides the prospects of testing directly as how much genetic change in a given time period has accompanied changes in material culture the sequencing of ancient Y chromosomes allows us also to better understand the rate at which mutations accumulate and get fixed over time. This review considers genome-scale evidence on ancient Y chromosome diversity that has recently started to accumulate in geographic areas favourable to DNA preservation. More specifically the review focuses on examples of regional continuity and change of the Y chromosome haplogroups in North Eurasia and in the New World.

Toward a consensus on SNP and STR mutation rates on the human Y-chromosome

The mutation rate on the Y-chromosome matters for estimating the time-to-the-most-recent-common-ancestor (TMRCA, i.e. haplogroup age) in population genetics, as well as for forensic, medical, and genealogical studies. Large-scale sequencing efforts have produced several independent estimates of Y-SNP mutation rates. Genealogical, or pedigree, rates tend to be slightly faster than evolutionary rates obtained from ancient DNA or calibrations using dated (pre)historical events. It is, therefore, suggested to report TMRCAs using an envelope defined by the average aDNA-based rate and the average pedigree-based rate. The current estimate of the "envelope rate" is 0.75-0.89 substitutions per billion base pairs per year. The available Y-SNP mutation rates can be applied to high-coverage data from the entire X-degenerate region, but other datasets may demand recalibrated rates. While a consensus on Y-SNP rates is approaching, the debate on Y-STR rates has continued for two decades, because multiple genealogical rates were consistent with each other but three times faster than the single evolutionary estimate. Applying Y-SNP and Y-STR rates to the same haplogroups recently helped to clarify the issue. Genealogical and evolutionary STR rates typically provide lower and upper bounds of the "true" (SNP-based) age. The genealogical rate often-but not always-works well for haplogroups less than 7000 years old. The evolutionary rate, although calibrated using recent events, inflates ages of young haplogroups and deflates the age of the entire Y-chromosomal tree, but often provides reasonable estimates for intermediate ages (old haplogroups). Future rate estimates and accumulating case studies should further clarify the Y-SNP rates.

Y-chromosomal sequences of diverse Indian populations and the ancestry of the Andamanese

We present 42 new Y-chromosomal sequences from diverse Indian tribal and non-tribal populations, including the Jarawa and Onge from the Andaman Islands, which are analysed within a calibrated Y-chromosomal phylogeny incorporating South Asian (in total 305 individuals) and worldwide (in total 1286 individuals) data from the 1000 Genomes Project. In contrast to the more ancient ancestry in the South than in the North that has been claimed, we detected very similar coalescence times within Northern and Southern non-tribal Indian populations. A closest neighbour analysis in the phylogeny showed that Indian populations have an affinity towards Southern European populations and that the time of divergence from these populations substantially predated the Indo-European migration into India, probably reflecting ancient shared ancestry rather than the Indo-European migration, which had little effect on Indian male lineages. Among the tribal populations, the Birhor (Austro-Asiatic-speaking) and Irula (Dravidian-speaking) are the nearest neighbours of South Asian non-tribal populations, with a common origin in the last few millennia. In contrast, the Riang (Tibeto-Burman-speaking) and Andamanese have their nearest neighbour lineages in East Asia. The Jarawa and Onge shared haplogroup D lineages with each other within the last ~7000 years, but had diverged from Japanese haplogroup D Y-chromosomes ~53000 years ago, most likely by a split from a shared ancestral population. This analysis suggests that Indian populations have complex ancestry which cannot be explained by a single expansion model.

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.

A Working Model of the Deep Relationships of Diverse Modern Human Genetic Lineages Outside of Africa

A major topic of interest in human prehistory is how the large-scale genetic structure of modern populations outside of Africa was established. Demographic models have been developed that capture the relationships among small numbers of populations or within particular geographical regions, but constructing a phylogenetic tree with gene flow events for a wide diversity of non-Africans remains a difficult problem. Here, we report a model that provides a good statistical fit to allele-frequency correlation patterns among East Asians, Australasians, Native Americans, and ancient western and northern Eurasians, together with archaic human groups. The model features a primary eastern/western bifurcation dating to at least 45,000 years ago, with Australasians nested inside the eastern clade, and a parsimonious set of admixture events. While our results still represent a simplified picture, they provide a useful summary of deep Eurasian population history that can serve as a null model for future studies and a baseline for further discoveries.

Detecting past male-mediated expansions using the Y chromosome

Males and females display biological differences that lead to a higher variance of offspring number in males, and this is frequently exacerbated in human societies by mating practices, and possibly by past socio-cultural circumstances. This implies that the genetic record might contain the imprint of past male-mediated expansions, which can be investigated by analysing the male-specific region of the Y chromosome (MSY). Here, we review studies that have used MSY data to infer such expansions. Sets of short-tandem repeats define haplotypes of very low average frequencies, but in a few cases, high-frequency haplotypes are observed, forming the core of descent clusters. Estimates of the ages of such clusters, together with geographical information, have been used to propose powerful historical founders, including Genghis Khan, although without direct supporting evidence. Resequencing of multi-megabase segments of MSY has allowed the construction of detailed phylogenies in which branch lengths are proportional to time, leading to the identification of lineage expansions in the last few millennia as well as the more distant past. Comparisons with maternally-inherited mitochondrial DNA sequence data allow the male specificity of some of these expansions to be demonstrated. These include expansions in Europe in the last ~5000 years that may be associated with a cultural shift during the Bronze Age, as well as expansions elsewhere in the world for which explanations from archaeological evidence are not yet clear.

Forensic use of Y-chromosome DNA: a general overview

The male-specific part of the human Y chromosome is widely used in forensic DNA analysis, particularly in cases where standard autosomal DNA profiling is not informative. A Y-chromosomal gene fragment is applied for inferring the biological sex of a crime scene trace donor. Haplotypes composed of Y-chromosomal short tandem repeat polymorphisms (Y-STRs) are used to characterise paternal lineages of unknown male trace donors, especially suitable when males and females have contributed to the same trace, such as in sexual assault cases. Y-STR haplotyping applied in crime scene investigation can (i) exclude male suspects from involvement in crime, (ii) identify the paternal lineage of male perpetrators, (iii) highlight multiple male contributors to a trace, and (iv) provide investigative leads for finding unknown male perpetrators. Y-STR haplotype analysis is employed in paternity disputes of male offspring and other types of paternal kinship testing, including historical cases, as well as in special cases of missing person and disaster victim identification involving men. Y-chromosome polymorphisms are applied for inferring the paternal bio-geographic ancestry of unknown trace donors or missing persons, in cases where autosomal DNA profiling is uninformative. In this overview, all different forensic applications of Y-chromosome DNA are described. To illustrate the necessity of forensic Y-chromosome analysis, the investigation of a prominent murder case is described, which initiated two changes in national forensic DNA legislation both covering Y-chromosome use, and was finally solved via an innovative Y-STR dragnet involving thousands of volunteers after 14 years. Finally, expectations for the future of forensic Y-chromosome DNA analysis are discussed.

Y chromosome palindromes and gene conversion

The presence of large and near-identical inverted repeat sequences (called palindromes) is a common feature of the constitutively haploid sex chromosomes of different species. Despite the fact palindromes originated in a non-recombining context, they have evolved a strong recombinational activity in the form of abundant arm-to-arm gene conversion. Their independent appearance in different species suggests they can have a profound biological significance that has yet to be fully clarified. It has been theorized that natural selection may have favored palindromic organization of male-specific genes and that the establishment of intra-palindrome gene conversion has strong adaptive significance. Arm-to-arm gene conversion allows the efficient removal of deleterious mutations, increases the fixation rate of beneficial mutations and has played an important role in modulating the equilibrium between gene loss and acquisition during Y chromosome evolution. Additionally, a palindromic organization of duplicates could favor the formation of unusual chromatin structures and could optimize the use of gene conversion as a mechanism to maintain the structural integrity of male-specific genes. In this review, we describe the structural features of palindromes on mammalian sex chromosomes and summarize different hypotheses regarding palindrome evolution and the functional benefits of arm-to-arm gene conversion on the unique haploid portion of the nuclear genome.

The Y chromosomes of the great apes

The great apes (orangutans, gorillas, chimpanzees, bonobos and humans) descended from a common ancestor around 13 million years ago, and since then their sex chromosomes have followed very different evolutionary paths. While great-ape X chromosomes are highly conserved, their Y chromosomes, reflecting the general lability and degeneration of this male-specific part of the genome since its early mammalian origin, have evolved rapidly both between and within species. Understanding great-ape Y chromosome structure, gene content and diversity would provide a valuable evolutionary context for the human Y, and would also illuminate sex-biased behaviours, and the effects of the evolutionary pressures exerted by different mating strategies on this male-specific part of the genome. High-quality Y-chromosome sequences are available for human and chimpanzee (and low-quality for gorilla). The chromosomes differ in size, sequence organisation and content, and while retaining a relatively stable set of ancestral single-copy genes, show considerable variation in content and copy number of ampliconic multi-copy genes. Studies of Y-chromosome diversity in other great apes are relatively undeveloped compared to those in humans, but have nevertheless provided insights into speciation, dispersal, and mating patterns. Future studies, including data from larger sample sizes of wild-born and geographically well-defined individuals, and full Y-chromosome sequences from bonobos, gorillas and orangutans, promise to further our understanding of population histories, male-biased behaviours, mutation processes, and the functions of Y-chromosomal genes.

Phylogeography of human Y-chromosome haplogroup Q3-L275 from an academic/citizen science collaboration

Background

The Y-chromosome haplogroup Q has three major branches: Q1, Q2, and Q3. Q1 is found in both Asia and the Americas where it accounts for about 90% of indigenous Native American Y-chromosomes; Q2 is found in North and Central Asia; but little is known about the third branch, Q3, also named Q1b-L275. Here, we combined the efforts of population geneticists and genetic genealogists to use the potential of full Y-chromosome sequencing for reconstructing haplogroup Q3 phylogeography and suggest possible linkages to events in population history.

Results

We analyzed 47 fully sequenced Y-chromosomes and reconstructed the haplogroup Q3 phylogenetic tree in detail. Haplogroup Q3-L275, derived from the oldest known split within Eurasian/American haplogroup Q, most likely occurred in West or Central Asia in the Upper Paleolithic period. During the Mesolithic and Neolithic epochs, Q3 remained a minor component of the West Asian Y-chromosome pool and gave rise to five branches (Q3a to Q3e), which spread across West, Central and parts of South Asia. Around 3–4 millennia ago (Bronze Age), the Q3a branch underwent a rapid expansion, splitting into seven branches, some of which entered Europe. One of these branches, Q3a1, was acquired by a population ancestral to Ashkenazi Jews and grew within this population during the 1st millennium AD, reaching up to 5% in present day Ashkenazi.

Conclusions

This study dataset was generated by a massive Y-chromosome genotyping effort in the genetic genealogy community, and phylogeographic patterns were revealed by a collaboration of population geneticists and genetic genealogists. This positive experience of collaboration between academic and citizen science provides a model for further joint projects. Merging data and skills of academic and citizen science promises to combine, respectively, quality and quantity, generalization and specialization, and achieve a well-balanced and careful interpretation of the paternal-side history of human populations.

Electronic supplementary material

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

Signatures of human European Palaeolithic expansion shown by resequencing of non-recombining X-chromosome segments

Human genetic diversity in Europe has been extensively studied using uniparentally inherited sequences (mitochondrial DNA (mtDNA) and the Y chromosome), which reveal very different patterns indicating sex-specific demographic histories. The X chromosome, haploid in males and inherited twice as often from mothers as from fathers, could provide insights into past female behaviours, but has not been extensively investigated. Here, we use HapMap single-nucleotide polymorphism data to identify genome-wide segments of the X chromosome in which recombination is historically absent and mutations are likely to be the only source of genetic variation, referring to these as phylogeographically informative haplotypes on autosomes and X chromosome (PHAXs). Three such sequences on the X chromosome spanning a total of ~49 kb were resequenced in 240 males from Europe, the Middle East and Africa at an average coverage of 181 ×. These PHAXs were confirmed to be essentially non-recombining across European samples. All three loci show highly homogeneous patterns across Europe and are highly differentiated from the African sample. Star-like structures of European-specific haplotypes in median-joining networks indicate past population expansions. Bayesian skyline plots and time-to-most-recent-common-ancestor estimates suggest expansions pre-dating the Neolithic transition, a finding that is more compatible with data on mtDNA than the Y chromosome, and with the female bias of X-chromosomal inheritance. This study demonstrates the potential of the use of X-chromosomal haplotype blocks, and the utility of the accurate ascertainment of rare variants for inferring human demographic history.

Microreact: visualizing and sharing data for genomic epidemiology and phylogeography

Visualization is frequently used to aid our interpretation of complex datasets. Within microbial genomics, visualizing the relationships between multiple genomes as a tree provides a framework onto which associated data (geographical, temporal, phenotypic and epidemiological) are added to generate hypotheses and to explore the dynamics of the system under investigation. Selected static images are then used within publications to highlight the key findings to a wider audience. However, these images are a very inadequate way of exploring and interpreting the richness of the data. There is, therefore, a need for flexible, interactive software that presents the population genomic outputs and associated data in a user-friendly manner for a wide range of end users, from trained bioinformaticians to front-line epidemiologists and health workers. Here, we present Microreact, a web application for the easy visualization of datasets consisting of any combination of trees, geographical, temporal and associated metadata. Data files can be uploaded to Microreact directly via the web browser or by linking to their location (e.g. from Google Drive/Dropbox or via API), and an integrated visualization via trees, maps, timelines and tables provides interactive querying of the data. The visualization can be shared as a permanent web link among collaborators, or embedded within publications to enable readers to explore and download the data. Microreact can act as an end point for any tool or bioinformatic pipeline that ultimately generates a tree, and provides a simple, yet powerful, visualization method that will aid research and discovery and the open sharing of datasets.

Human Y Chromosome Haplogroup N: A Non-trivial Time-Resolved Phylogeography that Cuts across Language Families

The paternal haplogroup (hg) N is distributed from southeast Asia to eastern Europe. The demographic processes that have shaped the vast extent of this major Y chromosome lineage across numerous linguistically and autosomally divergent populations have previously been unresolved. On the basis of 94 high-coverage re-sequenced Y chromosomes, we establish and date a detailed hg N phylogeny. We evaluate geographic structure by using 16 distinguishing binary markers in 1,631 hg N Y chromosomes from a collection of 6,521 samples from 56 populations. The more southerly distributed sub-clade N4 emerged before N2a1 and N3, found mostly in the north, but the latter two display more elaborate branching patterns, indicative of regional contrasts in recent expansions. In particular, a number of prominent and well-defined clades with common N3a3'6 ancestry occur in regionally dissimilar northern Eurasian populations, indicating almost simultaneous regional diversification and expansion within the last 5,000 years. This patrilineal genetic affinity is decoupled from the associated higher degree of language diversity.

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.

Refining the Y chromosome phylogeny with southern African sequences

The recent availability of large-scale sequence data for the human Y chromosome has revolutionized analyses of and insights gained from this non-recombining, paternally inherited chromosome. However, the studies to date focus on Eurasian variation, and hence the diversity of early-diverging branches found in Africa has not been adequately documented. Here, we analyze over 900 kb of Y chromosome sequence obtained from 547 individuals from southern African Khoisan- and Bantu-speaking populations, identifying 232 new sequences from basal haplogroups A and B. We identify new clades in the phylogeny, an older age for the root, and substantially older ages for some individual haplogroups. Furthermore, while haplogroup B2a is traditionally associated with the spread of Bantu speakers, we find that it probably also existed in Khoisan groups before the arrival of Bantu speakers. Finally, there is pronounced variation in branch length between major haplogroups; in particular, haplogroups associated with Bantu speakers have significantly longer branches. Technical artifacts cannot explain this branch length variation, which instead likely reflects aspects of the demographic history of Bantu speakers, such as recent population expansion and an older average paternal age. The influence of demographic factors on branch length variation has broader implications both for the human Y phylogeny and for similar analyses of other species.

New native South American Y chromosome lineages

Many single-nucleotide polymorphisms (SNPs) in the non-recombining region of the human Y chromosome have been described in the last decade. High-coverage sequencing has helped to characterize new SNPs, which has in turn increased the level of detail in paternal phylogenies. However, these paternal lineages still provide insufficient information on population history and demography, especially for Native Americans. The present study aimed to identify informative paternal sublineages derived from the main founder lineage of the Americas-haplogroup Q-L54-in a sample of 1841 native South Americans. For this purpose, we used a Y-chromosomal genotyping multiplex platform and conventional genotyping methods to validate 34 new SNPs that were identified in the present study by sequencing, together with many Y-SNPs previously described in the literature. We updated the haplogroup Q phylogeny and identified two new Q-M3 and three new Q-L54*(xM3) sublineages defined by five informative SNPs, designated SA04, SA05, SA02, SA03 and SA29. Within the Q-M3, sublineage Q-SA04 was mostly found in individuals from ethnic groups belonging to the Tukanoan linguistic family in the northwest Amazon, whereas sublineage Q-SA05 was found in Peruvian and Bolivian Amazon ethnic groups. Within Q-L54*, the derived sublineages Q-SA03 and Q-SA02 were exclusively found among Coyaima individuals (Cariban linguistic family) from Colombia, while Q-SA29 was found only in Maxacali individuals (Jean linguistic family) from southeast Brazil. Furthermore, we validated the usefulness of several published SNPs among indigenous South Americans. This new Y chromosome haplogroup Q phylogeny offers an informative paternal genealogy to investigate the pre-Columbian history of South America.Journal of Human Genetics advance online publication, 31 March 2016; doi:10.1038/jhg.2016.26.