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

Genomic insights into Neolithic founding paternal lineages around the Qinghai-Xizang Plateau using integrated YanHuang resource

Indigenous populations of the Qinghai-Xizang Plateau exhibit unique high-altitude adaptations, especially within Tibeto-Burman (TB) groups. However, the paternal genetic heritage of eastern Plateau regions remains less explored. We present one integrative Y chromosome dataset of 9,901 modern and ancient individuals, including whole Y chromosome sequences from 1,297 individuals and extensive Y-SNP/STR genotype data. We reveal the Paleolithic common origin and following divergence of Qinghai-Xizang Plateau ancestors from East Asian lowlands, marked by subsequent isolation and Holocene expansion involving local hunter-gatherers and millet-farming communities. We identified two key TB-related founding lineages, D-Z31591 and O-CTS4658, which underwent significant expansions around 5,000 years ago on the Qinghai-Xizang Plateau and its eastern Tibetan-Yi Corridor. The genetic legacy of these TB lineages highlights crucial migration pathways linking the Plateau and lowland southwestern China. Our findings align paternal genetic structures with East Asian geography and linguistic groups, underscoring the utility of Y chromosome analyses in unraveling complex paternal histories.

Development and validation of YARN: A novel SE-400 MPS kit for East Asian paternal lineage analysis

Y-chromosomal short tandem repeat polymorphisms (Y-STRs) and Y-chromosomal single nucleotide polymorphisms (Y-SNPs) are valuable genetic markers used in paternal lineage identification and population genetics. Currently, there is a lack of an effective panel that integrates Y-STRs and Y-SNPs for studying paternal lineages, particularly in East Asian populations. Hence, we developed a novel Y-chromosomal targeted panel called YARN (Y-chromosome Ancestry and Region Network) based on multiplex PCR and a single-end 400 massive parallel sequencing (MPS) strategy, consisting of 44 patrilineage Y-STRs and 260 evolutionary Y-SNPs. A total of 386 reactions were validated for the effectiveness and applicability of YARN according to SWGDAM validation guidelines, including sensitivity (with a minimum input gDNA of 0.125 ng), mixture identification (ranging from 1:1-1:10), PCR inhibitor testing (using substances such as 50 μM hematin, 100 μM hemoglobin, 100 μM humic acid, and 2.5 mM indigo dye), species specificity (successfully distinguishing humans from other animals), repeatability study (achieved 100% accuracy), and concordance study (with 99.91% accuracy for 1121 Y-STR alleles). Furthermore, we conducted a pilot study using YARN in a cohort of 484 Han Chinese males from Huaiji County, Zhaoqing City, Guangdong, China (GDZQHJ cohort). In this cohort, we identified 52 different Y-haplogroups and 73 different surnames. We found weak to moderate correlations between the Y-haplogroups, Chinese surnames, and geographical locations of the GDZQHJ cohort (with λ values ranging from 0.050 to 0.340). However, when we combined two different categories into a new independent variable, we observed stronger correlations (with λ values ranging from 0.617 to 0.754). Overall, the YARN panel, which combines Y-STR and Y-SNP genetic markers, meets forensic DNA quality assurance guidelines and holds potential for East Asian geographical origin inference and paternal lineage analysis.

Male Pedigree Toolbox: A Versatile Software for Y-STR Data Analyses

Y-chromosomal short tandem repeats (Y-STRs) are widely used in forensic, genealogical, and population genetics. With the recent increase in the number of rapidly mutating (RM) Y-STRs, an unprecedented level of male differentiation can be achieved, widening and improving the applications of Y-STRs in various fields, including forensics. The growing complexity of Y-STR data increases the need for automated data analyses, but dedicated software tools are scarce. To address this, we present the Male Pedigree Toolbox (MPT), a software tool for the automated analysis of Y-STR data in the context of patrilineal genealogical relationships. The MPT can estimate mutation rates and male relative differentiation rates from input Y-STR pedigree data. It can aid in determining ancestral haplotypes within a pedigree and visualize the genetic variation within pedigrees in all branches of family trees. Additionally, it can provide probabilistic classifications using machine learning, helping to establish or prove the structure of the pedigree and the level of relatedness between males, even for closely related individuals with highly similar haplotypes. The tool is flexible and easy to use and can be adjusted to any set of Y-STR markers by modifying the intuitive input file formats. We introduce the MPT software tool v1.0 and make it publicly available with the goal of encouraging and supporting forensic, genealogical, and other geneticists in utilizing the full potential of Y-STRs for both research purposes and practical applications, including criminal casework.

Contrasting genetic trajectories of endangered and expanding red fox populations in the western U.S

As anthropogenic disturbances continue to drive habitat loss and range contractions, the maintenance of evolutionary processes will increasingly require targeting measures to the population level, even for common and widespread species. Doing so requires detailed knowledge of population genetic structure, both to identify populations of conservation need and value, as well as to evaluate suitability of potential donor populations. We conducted a range-wide analysis of the genetic structure of red foxes in the contiguous western U.S., including a federally endangered distinct population segment of the Sierra Nevada subspecies, with the objectives of contextualizing field observations of relative scarcity in the Pacific mountains and increasing abundance in the cold desert basins of the Intermountain West. Using 31 autosomal microsatellites, along with mitochondrial and Y-chromosome markers, we found that populations of the Pacific mountains were isolated from one another and genetically depauperate (e.g., estimated Ne range = 3–9). In contrast, red foxes in the Intermountain regions showed relatively high connectivity and genetic diversity. Although most Intermountain red foxes carried indigenous western matrilines (78%) and patrilines (85%), the presence of nonindigenous haplotypes at lower elevations indicated admixture with fur-farm foxes and possibly expanding midcontinent populations as well. Our findings suggest that some Pacific mountain populations could likely benefit from increased connectivity (i.e., genetic rescue) but that nonnative admixture makes expanding populations in the Intermountain basins a non-ideal source. However, our results also suggest contact between Pacific mountain and Intermountain basin populations is likely to increase regardless, warranting consideration of risks and benefits of proactive measures to mitigate against unwanted effects of Intermountain gene flow.

Improving the regional Y-STR haplotype resolution utilizing haplogroup-determining Y-SNPs and the application of machine learning in Y-SNP haplogroup prediction in a forensic Y-STR database: A pilot study on male Chinese Yunnan Zhaoyang Han population

Improving the resolution of the current widely used Y-chromosomal short tandem repeat (Y-STR) dataset is of great importance for forensic investigators, and the current approach is limited, except for the addition of more Y-STR loci. In this research, a regional Y-DNA database was investigated to improve the Y-STR haplotype resolution utilizing a Y-SNP Pedigree Tagging System that includes 24 Y-chromosomal single nucleotide polymorphism (Y-SNP) loci. This pilot study was conducted in the Chinese Yunnan Zhaoyang Han population, and 3473 unrelated male individuals were enrolled. Based on data on the male haplogroups under different panels, the matched or near-matching (NM) Y-STR haplotype pairs from different haplogroups indicated the critical roles of haplogroups in improving the regional Y-STR haplotype resolution. A classic median-joining network analysis was performed using Y-STR or Y-STR/Y-SNP data to reconstruct population substructures, which revealed the ability of Y-SNPs to correct misclassifications from Y-STRs. Additionally, population substructures were reconstructed using multiple unsupervised or supervised dimensionality reduction methods, which indicated the potential of Y-STR haplotypes in predicting Y-SNP haplogroups. Haplogroup prediction models were built based on nine publicly accessible machine-learning (ML) approaches. The results showed that the best prediction accuracy score could reach 99.71% for major haplogroups and 98.54% for detailed haplogroups. Potential influences on prediction accuracy were assessed by adjusting the Y-STR locus numbers, selecting Y-STR loci with various mutabilities, and performing data processing. ML-based predictors generally presented a better prediction accuracy than two available predictors (Nevgen and EA-YPredictor). Three tree models were developed based on the Yfiler Plus panel with unprocessed input data, which showed their strong generalization ability in classifying various Chinese Han subgroups (validation dataset). In conclusion, this study revealed the significance and application prospects of Y-SNP haplogroups in improving regional Y-STR databases. Y-SNP haplogroups can be used to discriminate NM Y-STR haplotype pairs, and it is important for forensic Y-STR databases to develop haplogroup prediction tools to improve the accuracy of biogeographic ancestry inferences.

CSYseq: The first Y-chromosome sequencing tool typing a large number of Y-SNPs and Y-STRs to unravel worldwide human population genetics

Male-specific Y-chromosome (chrY) polymorphisms are interesting components of the DNA for population genetics. While single nucleotide polymorphisms (Y-SNPs) indicate distant evolutionary ancestry, short tandem repeats (Y-STRs) are able to identify close familial kinships. Detailed chrY analysis provides thus both biogeographical background information as paternal lineage identification. The rapid advancement of high-throughput massive parallel sequencing (MPS) technology in the past decade has revolutionized genetic research. Using MPS, single-base information of both Y-SNPs as Y-STRs can be analyzed in a single assay typing multiple samples at once. In this study, we present the first extensive chrY-specific targeted resequencing panel, the 'CSYseq', which simultaneously identifies slow mutating Y-SNPs as evolution markers and rapid mutating Y-STRs as patrilineage markers. The panel was validated by paired-end sequencing of 130 males, distributed over 65 deep-rooted pedigrees covering 1,279 generations. The CSYseq successfully targets 15,611 Y-SNPs including 9,014 phylogenetic informative Y-SNPs to identify 1,443 human evolutionary Y-subhaplogroup lineages worldwide. In addition, the CSYseq properly targets 202 Y-STRs, including 81 slow, 68 moderate, 27 fast and 26 rapid mutating Y-STRs to individualize close paternal relatives. The targeted chrY markers cover a high average number of reads (Y-SNP = 717, Y-STR = 150), easy interpretation, powerful discrimination capacity and chrY specificity. The CSYseq is interesting for research on different time scales: to identify evolutionary ancestry, to find distant family and to discriminate closely related males. Therefore, this panel serves as a unique tool valuable for a wide range of genetic-genealogical applications in interdisciplinary research within evolutionary, population, molecular, medical and forensic genetics.

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.

Improved Models of Coalescence Ages of Y-DNA Haplogroups

Databases of commercial DNA-testing companies now contain more customers with sequenced DNA than any completed academic study, leading to growing interest from academic and forensic entities. An important result for both these entities and the test takers themselves is how closely two individuals are related in time, as calculated through one or more molecular clocks. For Y-DNA, existing interpretations of these clocks are insufficiently accurate to usefully measure relatedness in historic times. In this article, I update the methods used to calculate coalescence ages (times to most-recent common ancestor, or TMRCAs) using a new, probabilistic statistical model that includes Y-SNP, Y-STR and ancilliary historical data, and provide examples of its use.

Genetic reconstruction and phylogenetic analysis by 193 Y-SNPs and 27 Y-STRs in a Chinese Yi ethnic group

Yi is the seventh-largest ethnic group in China and features mountainous regional characteristics. The Liangshan Yi Autonomous Prefecture is the largest Yi agglomeration with isolated geographical conditions, profoundly impeding genetic communication. Here, we investigated 427 unrelated males of Liangshan from 193 Y-chromosome single nucleotide polymorphisms (Y-SNPs) and 27 Y-chromosome short tandem repeats (Y-STRs) to reveal the genetic structure and paternal phylogeny of the group. The haplogroup diversity reached 0.9169 with 46 different subhaplogroups by 193 Y-SNPs, and the haplotype diversity reached 0.9999 by 27 Y-STR loci. Multidimensional scaling (MDS), N-J tree, and Network were constructed to decipher and visualize the genetic relations between Yi and worldwide groups. Our results revealed: (1) the Network by Y-STRs and Y-SNPs showed the haplogroup D1a1a-M15 in the Liangshan Yi population was a ramification of Tibetan groups' expansion from west to east on the plateau; (2) the haplogroup distribution and the mismatch mutation analysis indicated the haplogroup O2a2b1a1a1a4a2-Z25929 of Liangshan Yi derived from manifold Southeast Asian immigrants; (3) a high-resolution Y-SNPs panel is vital to depict accurate paternal derivations and build an integrated and refining genetic structure of ethnic groups.

Y chromosome analysis for common surnames in the Japanese male population

For many years of Japan's long history, Japanese surnames have been handed down patrilineally. This study investigated relations between major surnames and Y chromosomal polymorphism among the Japanese male population. To analyze genetic phylogeny in namesakes, the Y-single nucleotide polymorphism (SNP) plus Y-short tandem repeat (STR) approach was employed. A haplogroup based on SNPs and haplotypes at 17 STR loci were typed in 567 unrelated volunteers recruited in Kanagawa, Japan. Samples covered 27 common surnames such as Satoh and Suzuki, each name having 10-55 bearers. Significant difference was found for SNP haplogroup compositions and a multidimensional scaling plot using STR haplotypes in several surname groups. By contrast, these common surnames displayed wide diversity with phylogenetic networks, suggesting that no genetic drift event has occurred in their history. In all, 22 descent clusters were found, as judgcriteria ed by ad hoc of groups within five mutational steps in the 15 STR loci with the same haplogroup. The times of the most recent common ancestor ranged from 279 to over 2577 years. According to the approximate millennium span of Japanese surname history, descent criteria are expected to be reasonable for grouping within four step-neighbors. High heterogeneity of common surnames resembles that observed for England and Spain, but not for Ireland. Our results highlight that common Japanese surnames consist of descent clusters and many singletons, reflecting a mixture of long-term bearers and short-term bearers among the population. The genetic study of this population revealed characteristic features of Japanese surnames.

Exploitation of a novel slowly mutating Y-STRs set and evaluation of slowly mutating Y-STRs plus Y-SNPs typing strategy in forensic genetics and evolutionary research

The Y-chromosome short tandem repeats (Y-STRs) loci with different mutation rates existing in the Y chromosome non-recombination region (NRY) allow to be applied in human forensics, genealogical researches, historical investigations and evolutionary studies. Currently, there is a high demand for pedigree search to narrow the scope of crime investigations. However, the commonly used Y-STRs kits generally contain Y-STRs with high mutation rates that could cause individuals from the same pedigree to display different haplotypes. Herein, we put forward a new strategy of Slowly Mutating (SM) Y-STRs plus Y-SNPs typing, which could not only improve the resolution and accuracy of pedigree search, but also be applicable to evolutionary research. First, we developed a nine SM Y-STRs assay by evaluating their mutation rates in 210 pedigrees. Then the gene diversity and efficiency of the SM Y-STRs and 172 Y-SNPs sets were investigated by 2304 unrelated males from 24 populations. Furthermore, network and time estimation analyses were performed to evaluate the new strategy's capability to reconstruct phylogenetic tree and reliability to infer the time to the most recent common ancestor (TMRCA). The nine SM Y-STRs assay even had a higher resolution and a comparable capacity of revealing population genetic differentiation compared to 172 Y-SNPs system. This new strategy could optimize the phylogenetic tree generated by commonly used Y-STR panels and obtain a quite consistent time estimations with the published dating.

Autosomal genetics and Y-chromosome haplogroup L1b-M317 reveal Mount Lebanon Maronites as a persistently non-emigrating population

Currently, there are 18 different religious communities living in Lebanon. While evolving primarily within Lebanon, these communities show a level of local isolation as demonstrated previously from their Y-haplogroup distributions. In order to trace the origins and migratory patterns that may have led to the genetic isolation and autosomal clustering in some of these communities we analyzed Y-chromosome STR and SNP sample data from 6327 individuals, in addition to whole genome autosomal sample data from 609 individuals, from Mount Lebanon and other surrounding communities. We observed Y chromosome L1b Levantine STR branching that occurred around 5000 years ago. Autosomal DNA analyses suggest that the North Lebanese Mountain Maronite community possesses an ancestral Fertile Crescent genetic component distinct from other populations in the region. We suggest that the Levantine L1b group split from the Caucasus ancestral group around 7300 years ago and migrated to the Levant. This event was distinct from the earlier expansions from the Caucasus region that contributed to the wider Levantine populations. Differential cultural adaption by populations from the North Lebanese Mountains are clearly aligned with the L1b haplotype STR haplogroup clusters, indicating pre-existing and persistent cultural barriers marked by the transmission of L1b lineages. Our findings highlight the value of uniparental haplogroups and STR haplotype data for elucidating biosocial events among these populations.

Profiling of 17 Y-STR loci in Mazandaran and Gilan provinces of Iran

Background/aim

The Y-chromosome mainly consists of heterochromatin regions that have a father-to-son inheritance. Short tandem repeat polymorphic (STRP) markers distributed all over the chromosome provide the opportunity for investigations in forensic medicine and ancestral lineage studies. Due to the existence of wide varieties of geographical and ethnic groups in Iran, studying Y-STRP markers is necessary for further applications. Here we investigated the provinces of Mazandaran and Gilan for the first time.

Materials and methods

Samples included 119 and 90 unrelated males from Mazandaran and Gilan, respectively. Using a PCR amplification kit, 17 Y-STRP markers were amplified and genotyping was conducted by capillary electrophoresis. Allele frequency, haplotype diversity (HD), and haplotype discrimination capacity (DC) were calculated. The populations were compared together and to neighboring countries including Afghanistan and Azerbaijan by F_ST index.

Results

A total of 204 unique haplotypes were observed. No uniqueness was observed between the two provinces. HD was 0.9993 and 0.9998 in Mazandaran and Gilan, respectively. DC was 0.9666 and 0.9888 for Mazandaran and Gilan, respectively. DYS385b and DYS391 had the most and least polymorphic content in both provinces, respectively. There was not a significant difference between these two provinces (F_ST = 0.0006 and P = 0.00) and neighboring countries.

Conclusion

The results highlight the effectiveness of these Y-STRP markers for male discrimination in the north of Iran. Using additional markers along with extended sample size would provide a better opportunity for removing matched haplotypes and introducing the best polymorphic markers in this specific population.

Comparison of Markov Chain Monte Carlo Software for the Evolutionary Analysis of Y-Chromosomal Microsatellite Data

The evolutionary analysis of genetic data is an important subject of modern bioscience, with practical applications in diverse fields. Parameters of interest in this context include effective population sizes, mutation rates, population growth rates and the times to most recent common ancestors. Studying Y-chromosomal microsatellite data, in particular, has proven useful to unravel the recent patrilineal history of Homo sapiens populations. We compared the individual analysis options and technical details of four software tools that are widely used for this purpose, namely BATWING, BEAST, IMa2 and LAMARC, all of which use Bayesian coalescent-based Markov chain Monte Carlo (MCMC) methods for parameter estimation. More specifically, we simulated datasets for either eight or 20 hypothetical Y-chromosomal microsatellites, assuming a mutation rate of 0.0030 per generation and a constant or exponentially increasing population size, and used these data to evaluate the parameter estimation capacity of each tool. The datasets comprised between 100 and 1000 samples. In addition to runtime, the practical utility of the tools of interest can also be expected to depend critically upon the convergence behavior of the actual MCMC implementation. In fact, we found that runtime increased, and convergence rate decreased, with increasing sample size as expected. BATWING performed best with respect to runtime and convergence behavior, but only supports simple evolutionary models. As regards the spectrum of evolutionary models covered, and also in terms of cross-platform usability, BEAST provided the greatest flexibility. Finally, IMa2 and LAMARC turned out best to incorporate elaborate migration models in the analysis process.

Forensic characteristics and genetic analysis of both 27 Y-STRs and 143 Y-SNPs in Eastern Han Chinese population

Y-chromosome short tandem repeat (Y-STR) and Y-chromosome single nucleotide polymorphism (Y-SNP) frequency distributions provide resources for assessment of male population stratification among world-wide populations. Currently, the Y-STR Haplotype Reference Database (YHRD) contains numerous Y-chromosome haplotype profiles from various populations and countries around the world. However, for many of the recently discovered and already phylogenetically mapped Y-SNPs, the population data are scarce. Herein, the typing of 27 Y-STRs (Yfiler Plus) and 143 Y-SNPs (self-designed Y-SNP panel) was performed on 1269 unrelated males from 11 Han Chinese populations. Haplogroup O-M175 was the most predominant haplogroup in our Han Chinese data, ranging from 67.34% (Henan Han) to 93.16% (Guangdong Han). The highest haplogroup diversity (0.967056) was observed in Heilongjiang Han, with a discrimination capacity (DC) value of 0.3723. The number of alleles at single-copy loci varied from 2 for DYS391 (Guangdong Han) to 16 for DYS518 (Henan Han). For the majority of the populations (8/11), both the haplotype diversity and DC values are 1.0000. Furthermore, genetic differentiations were observed between Northern and Southern Han Chinese. These genetic differences were mainly reflected in haplogroup distribution and frequency, and they were confirmed by statistical analysis.

Estimating Y-Str Mutation Rates and Tmrca Through Deep-Rooting Italian Pedigrees

In the population genomics era, the study of Y-chromosome variability is still of the greatest interest for several fields ranging from molecular anthropology to forensics and genetic genealogy. In particular, mutation rates of Y-chromosomal Short Tandem Repeats markers (Y-STRs) are key parameters for different interdisciplinary applications. Among them, testing the patrilineal relatedness between individuals and calculating their Time of Most Recent Common Ancestors (TMRCAs) are of the utmost importance. To provide new valuable estimates and to address these issues, we typed 47 Y-STRs (comprising Yfiler, PowerPlex23 and YfilerPlus loci, the recently defined Rapidly Mutating [RM] panel and 11 additional markers often used in genetic genealogical applications) in 135 individuals belonging to 66 deep-rooting paternal genealogies from Northern Italy. Our results confirmed that the genealogy approach is an effective way to obtain reliable Y-STR mutation rate estimates even with a limited number of samples. Moreover, they showed that the impact of multi-step mutations and backmutations is negligible within the temporal scale usually adopted by forensic and genetic genealogy analyses. We then detected a significant association between the number of mutations within genealogies and observed TMRCAs. Therefore, we compared observed and expected TMRCAs by implementing a Bayesian procedure originally designed by Walsh (2001) and showed that the method yields a good performance (up to 96.72%), especially when using the Infinite Alleles Model (IAM).

Relating Clans Ao and Aisin Gioro from northeast China by whole Y-chromosome sequencing

The Y-chromosome haplogroup C2b1a3a2-F8951 is the paternal lineage of the Aisin Gioro clan, the most important brother branch of the famous Mongolic-speaking population characteristic haplogroup C2*-Star Cluster (C2b1a3a1-F3796). However, investigations on its internal phylogeny are still limited. In this study, we used whole Y-chromosome sequencing to update its phylogenetic tree. In the revised tree, C2b1a3a2-F8951 and C2*-Star Cluster differentiated 3852 years ago (95% CI = 3295-4497). Approximately 3558 years ago (95% CI = 3013-4144), C2b1a3a2-F8951 was divided into two main subclades, C2b1a3a2a-F14753 and C2b1a3a2b-F5483. Currently, samples of C2b1a3a2-F8951 were mainly from the House of Aisin Gioro clan, the Ao family from Daur and some individuals mainly from northeast China. Although other haplogroups are also found in the Ao family, including C2b1a2-M48, C2b1a3a1-F3796, C2a1b-F845, and N1c-M178, the haplogroup C2b1a3a2-F8951 is still the most distinct genetic component. For haplogroup C2b1a3a2-F8951, the time of the most recent common ancestor of the House of Aisin Gioro clan and the Ao family were both very late, just a few hundred years ago. Some family-specific Y-SNPs of the House of Aisin Gioro and the Ao family were also discovered. This revision evidently improved the resolving power of Y-chromosome phylogeny in northeast Asia, deepening our understanding of the origin of these two families, even the Mongolic-speaking population.

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.

Dispersals of the Siberian Y-chromosome haplogroup Q in Eurasia

The human Y-chromosome has proven to be a powerful tool for tracing the paternal history of human populations and genealogical ancestors. The human Y-chromosome haplogroup Q is the most frequent haplogroup in the Americas. Previous studies have traced the origin of haplogroup Q to the region around Central Asia and Southern Siberia. Although the diversity of haplogroup Q in the Americas has been studied in detail, investigations on the diffusion of haplogroup Q in Eurasia and Africa are still limited. In this study, we collected 39 samples from China and Russia, investigated 432 samples from previous studies of haplogroup Q, and analyzed the single nucleotide polymorphism (SNP) subclades Q1a1a1-M120, Q1a2a1-L54, Q1a1b-M25, Q1a2-M346, Q1a2a1a2-L804, Q1a2b2-F1161, Q1b1a-M378, and Q1b1a1-L245. Through NETWORK and BATWING analyses, we found that the subclades of haplogroup Q continued to disperse from Central Asia and Southern Siberia during the past 10,000 years. Apart from its migration through the Beringia to the Americas, haplogroup Q also moved from Asia to the south and to the west during the Neolithic period, and subsequently to the whole of Eurasia and part of Africa.

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.

Paternal phylogeographic structure of the brown bear (Ursus arctos) in northeastern Asia and the effect of male-mediated gene flow to insular populations

BACKGROUND

Sex-biased dispersal is widespread among mammals, including the brown bear (Ursus arctos). Previous phylogeographic studies of the brown bear based on maternally inherited mitochondrial DNA have shown intraspecific genetic structuring around the northern hemisphere. The brown bears on Hokkaido Island, northern Japan, comprise three distinct maternal lineages that presumably immigrated to the island from the continent in three different periods. Here, we investigate the paternal genetic structure across northeastern Asia and assess the connectivity among and within intraspecific populations in terms of male-mediated gene flow.

RESULTS

We analyzed paternally inherited Y-chromosomal DNA sequence data and Y-linked microsatellite data of 124 brown bears from Hokkaido, the southern Kuril Islands (Kunashiri and Etorofu), Sakhalin, and continental Eurasia (Kamchatka Peninsula, Ural Mountains, European Russia, and Tibet). The Hokkaido brown bear population is paternally differentiated from, and lacked recent genetic connectivity with, the continental Eurasian and North American populations. We detected weak spatial genetic structuring of the paternal lineages on Hokkaido, which may have arisen through male-mediated gene flow among natal populations. In addition, our results suggest that the different dispersal patterns between male and female brown bears, combined with the founder effect and subsequent genetic drift, contributed to the makeup of the Etorofu Island population, in which the maternal and paternal lineages show different origins.

CONCLUSIONS

Brown bears on Hokkaido and the adjacent southern Kuril Islands experienced different maternal and paternal evolutionary histories. Our results indicate that sex-biased dispersal has played a significant role in the evolutionary history of the brown bear in continental populations and in peripheral insular populations, such as on Hokkaido, the southern Kuril Islands, and Sakhalin.

Patrilineal background of the She minority population from Chaoshan Fenghuang Mountain, an isolated mountain region, in China

The She ethnic minority population is distributed in southern China. The origin of the She population has been controversial. The purpose of this work was to investigate the genomic diversity of She. The Chaoshan She population living in the Chaoshan Fenghuang mountain is a relatively isolated population. We detected 14 Y chromosome biallelic markers (Y-SNPs) and 6 Y chromosome short tandem repeat (Y-STR) loci in Chaoshan She people. Y-SNP analysis showed the Chaoshan She was closely related to the Chaoshan Hakka, Chaoshanese, Tujia and Gaoshan national minority. Compared with the Fujian She, the Chaoshan She maintained a more southern native genetic structure. Y-STR analysis revealed the Chaoshan She population was more closely related to the Hakka population than the other Hans. We concluded the Chaoshan She population had a closer genetic relationship with the southern national minority and Hakka Han and it may be representative of She ancestors' patrilineal genetic structure.

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.

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.

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.

Genetic trail for the early migrations of Aisin Gioro, the imperial house of the Qing dynasty

The House of Aisin Gioro, the imperial clan of Qing dynasty (1644-1911), affected the history of China and the formation of Manchu ethnicity greatly. However, owing to the lack of historical records and archeological evidences, the origin of the House of Aisin Gioro remains ambiguous. To clarify the origin of Aisin Gioro clan, we conducted whole Y-chromosome sequencing on three samples and Y-single-nucleotide polymorphism (Y-SNP) genotyping on other four samples beside those reported in previous work. We confirmed that the paternal lineage of the Aisin Gioro clan belongs to haplogroup C3b1a3a2-F8951, a brother branch of C3*-Star Cluster (currently named as C3b1a3a1-F3796, once linked to Genghis Khan), which is quite different from the predominant lineage C3c-M48 in other Tungusic-speaking populations. We also determined a series of unique Y-SNP markers for the Aisin Gioro clan. Diversity analyses of haplogroup C3b1a3a2-F8951 revealed the early migration of the ancestors of the Aisin Gioro clan from the middle reaches of Amur River to their later settlement in southeastern Manchuria. Hence, our results suggest that the Aisin Gioro clan may be descendants of ancient populations in Transbaikal region and closely related to origin of current Daur populations. Our research indicated that detailed research of stemma and deep sequencing of Y chromosomes are helpful to explore the prehistoric activities of populations lacking historical records and archeological evidences.

Mutation Rates and Discriminating Power for 13 Rapidly-Mutating Y-STRs between Related and Unrelated Individuals

Rapidly Mutating Y-STRs (RM Y-STRs) were recently introduced in forensics in order to increase the differentiation of Y-chromosomal profiles even in case of close relatives. We estimate RM Y-STRs mutation rates and their power to discriminate between related individuals by using samples extracted from a wide set of paternal pedigrees and by comparing RM Y-STRs results with those obtained from the Y-filer set. In addition, we tested the ability of RM Y-STRs to discriminate between unrelated individuals carrying the same Y-filer haplotype, using the haplogroup R-M269 (reportedly characterised by a strong resemblance in Y-STR profiles) as a case study. Our results, despite confirming the high mutability of RM Y-STRs, show significantly lower mutation rates than reference germline ones. Consequently, their power to discriminate between related individuals, despite being higher than the one of Y-filer, does not seem to improve significantly the performance of the latter. On the contrary, when considering R-M269 unrelated individuals, RM Y-STRs reveal significant discriminatory power and retain some phylogenetic signal, allowing the correct classification of individuals for some R-M269-derived sub-lineages. These results have important implications not only for forensics, but also for molecular anthropology, suggesting that RM Y-STRs are useful tools for exploring subtle genetic variability within Y-chromosomal haplogroups.

Population-Scale Sequencing Data Enable Precise Estimates of Y-STR Mutation Rates

Short tandem repeats (STRs) are mutation-prone loci that span nearly 1% of the human genome. Previous studies have estimated the mutation rates of highly polymorphic STRs by using capillary electrophoresis and pedigree-based designs. Although this work has provided insights into the mutational dynamics of highly mutable STRs, the mutation rates of most others remain unknown. Here, we harnessed whole-genome sequencing data to estimate the mutation rates of Y chromosome STRs (Y-STRs) with 2-6 bp repeat units that are accessible to Illumina sequencing. We genotyped 4,500 Y-STRs by using data from the 1000 Genomes Project and the Simons Genome Diversity Project. Next, we developed MUTEA, an algorithm that infers STR mutation rates from population-scale data by using a high-resolution SNP-based phylogeny. After extensive intrinsic and extrinsic validations, we harnessed MUTEA to derive mutation-rate estimates for 702 polymorphic STRs by tracing each locus over 222,000 meioses, resulting in the largest collection of Y-STR mutation rates to date. Using our estimates, we identified determinants of STR mutation rates and built a model to predict rates for STRs across the genome. These predictions indicate that the load of de novo STR mutations is at least 75 mutations per generation, rivaling the load of all other known variant types. Finally, we identified Y-STRs with potential applications in forensics and genetic genealogy, assessed the ability to differentiate between the Y chromosomes of father-son pairs, and imputed Y-STR genotypes.

Deep Roots for Aboriginal Australian Y Chromosomes

Australia was one of the earliest regions outside Africa to be colonized by fully modern humans, with archaeological evidence for human presence by 47,000 years ago (47 kya) widely accepted [1, 2]. However, the extent of subsequent human entry before the European colonial age is less clear. The dingo reached Australia about 4 kya, indirectly implying human contact, which some have linked to changes in language and stone tool technology to suggest substantial cultural changes at the same time [3]. Genetic data of two kinds have been proposed to support gene flow from the Indian subcontinent to Australia at this time, as well: first, signs of South Asian admixture in Aboriginal Australian genomes have been reported on the basis of genome-wide SNP data [4]; and second, a Y chromosome lineage designated haplogroup C(∗), present in both India and Australia, was estimated to have a most recent common ancestor around 5 kya and to have entered Australia from India [5]. Here, we sequence 13 Aboriginal Australian Y chromosomes to re-investigate their divergence times from Y chromosomes in other continents, including a comparison of Aboriginal Australian and South Asian haplogroup C chromosomes. We find divergence times dating back to ∼50 kya, thus excluding the Y chromosome as providing evidence for recent gene flow from India into Australia.

Antiquity and diversity of aboriginal Australian Y-chromosomes

OBJECTIVE

Understanding the origins of Aboriginal Australians is crucial in reconstructing the evolution and spread of Homo sapiens as evidence suggests they represent the descendants of the earliest group to leave Africa. This study analyzed a large sample of Y-chromosomes to answer questions relating to the migration routes of their ancestors, the age of Y-haplogroups, date of colonization, as well as the extent of male-specific variation.

METHODS

Knowledge of Y-chromosome variation among Aboriginal Australians is extremely limited. This study examined Y-SNP and Y-STR variation among 657 self-declared Aboriginal males from locations across the continent. 17 Y-STR loci and 47 Y-SNPs spanning the Y-chromosome phylogeny were typed in total.

RESULTS

The proportion of non-indigenous Y-chromosomes of assumed Eurasian origin was high, at 56%. Y lineages of indigenous Sahul origin belonged to haplogroups C-M130*(xM8,M38,M217,M347) (1%), C-M347 (19%), K-M526*(xM147,P308,P79,P261,P256,M231,M175,M45,P202) (12%), S-P308 (12%), and M-M186 (0.9%). Haplogroups C-M347, K-M526*, and S-P308 are Aboriginal Australian-specific. Dating of C-M347, K-M526*, and S-P308 indicates that all are at least 40,000 years old, confirming their long-term presence in Australia. Haplogroup C-M347 comprised at least three sub-haplogroups: C-DYS390.1del, C-M210, and the unresolved paragroup C-M347*(xDYS390.1del,M210).

CONCLUSIONS

There was some geographic structure to the Y-haplogroup variation, but most haplogroups were present throughout Australia. The age of the Australian-specific Y-haplogroups suggests New Guineans and Aboriginal Australians have been isolated for over 30,000 years, supporting findings based on mitochondrial DNA data. Our data support the hypothesis of more than one route (via New Guinea) for males entering Sahul some 50,000 years ago and give no support for colonization events during the Holocene, from either India or elsewhere.

Great ape Y Chromosome and mitochondrial DNA phylogenies reflect subspecies structure and patterns of mating and dispersal

The distribution of genetic diversity in great ape species is likely to have been affected by patterns of dispersal and mating. This has previously been investigated by sequencing autosomal and mitochondrial DNA (mtDNA), but large-scale sequence analysis of the male-specific region of the Y Chromosome (MSY) has not yet been undertaken. Here, we use the human MSY reference sequence as a basis for sequence capture and read mapping in 19 great ape males, combining the data with sequences extracted from the published whole genomes of 24 additional males to yield a total sample of 19 chimpanzees, four bonobos, 14 gorillas, and six orangutans, in which interpretable MSY sequence ranges from 2.61 to 3.80 Mb. This analysis reveals thousands of novel MSY variants and defines unbiased phylogenies. We compare these with mtDNA-based trees in the same individuals, estimating time-to-most-recent common ancestor (TMRCA) for key nodes in both cases. The two loci show high topological concordance and are consistent with accepted (sub)species definitions, but time depths differ enormously between loci and (sub)species, likely reflecting different dispersal and mating patterns. Gorillas and chimpanzees/bonobos present generally low and high MSY diversity, respectively, reflecting polygyny versus multimale–multifemale mating. However, particularly marked differences exist among chimpanzee subspecies: The western chimpanzee MSY phylogeny has a TMRCA of only 13.2 (10.8–15.8) thousand years, but that for central chimpanzees exceeds 1 million years. Cross-species comparison within a single MSY phylogeny emphasizes the low human diversity, and reveals species-specific branch length variation that may reflect differences in long-term generation times.

Genetic Diversity in the Lesser Antilles and Its Implications for the Settlement of the Caribbean Basin

Historical discourses about the Caribbean often chronicle West African and European influence to the general neglect of indigenous people's contributions to the contemporary region. Consequently, demographic histories of Caribbean people prior to and after European contact are not well understood. Although archeological evidence suggests that the Lesser Antilles were populated in a series of northward and eastern migratory waves, many questions remain regarding the relationship of the Caribbean migrants to other indigenous people of South and Central America and changes to the demography of indigenous communities post-European contact. To explore these issues, we analyzed mitochondrial DNA and Y-chromosome diversity in 12 unrelated individuals from the First Peoples Community in Arima, Trinidad, and 43 unrelated Garifuna individuals residing in St. Vincent. In this community-sanctioned research, we detected maternal indigenous ancestry in 42% of the participants, with the remainder having haplotypes indicative of African and South Asian maternal ancestry. Analysis of Y-chromosome variation revealed paternal indigenous American ancestry indicated by the presence of haplogroup Q-M3 in 28% of the male participants from both communities, with the remainder possessing either African or European haplogroups. This finding is the first report of indigenous American paternal ancestry among indigenous populations in this region of the Caribbean. Overall, this study illustrates the role of the region's first peoples in shaping the genetic diversity seen in contemporary Caribbean populations.

Phylogeographic Refinement and Large Scale Genotyping of Human Y Chromosome Haplogroup E Provide New Insights into the Dispersal of Early Pastoralists in the African Continent

Haplogroup E, defined by mutation M40, is the most common human Y chromosome clade within Africa. To increase the level of resolution of haplogroup E, we disclosed the phylogenetic relationships among 729 mutations found in 33 haplogroup DE Y-chromosomes sequenced at high coverage in previous studies. Additionally, we dissected the E-M35 subclade by genotyping 62 informative markers in 5,222 samples from 118 worldwide populations. The phylogeny of haplogroup E showed novel features compared with the previous topology, including a new basal dichotomy. Within haplogroup E-M35, we resolved all the previously known polytomies and assigned all the E-M35* chromosomes to five new different clades, all belonging to a newly identified subhaplogroup (E-V1515), which accounts for almost half of the E-M35 chromosomes from the Horn of Africa. Moreover, using a Bayesian phylogeographic analysis and a single nucleotide polymorphism-based approach we localized and dated the origin of this new lineage in the northern part of the Horn, about 12 ka. Time frames, phylogenetic structuring, and sociogeographic distribution of E-V1515 and its subclades are consistent with a multistep demic spread of pastoralism within north-eastern Africa and its subsequent diffusion to subequatorial areas. In addition, our results increase the discriminative power of the E-M35 haplogroup for use in forensic genetics through the identification of new ancestry-informative markers.

Evaluating the Y chromosomal STR dating in deep-rooting pedigrees

Background

Y chromosomal short tandem repeat (STR) has been used in time estimations for single nucleotide polymorphism (SNP) lineages or eminent persons. But to choose which mutation rate and estimation method in the Y chromosome dating is controversial, since different rates and methods can result in several-fold deviation.

Findings

We used two deep-rooting pedigrees with full records and reliable dates to directly evaluate the Y chromosomal STR mutation rates and dating methods. We found that the Y chromosomal genealogical mutation rates (OMRB and lmMR) in BATWING method can give the best-fit estimation for historical lineage dating.

Conclusions

This study validated a very efficient and reliable way for genealogy and historical anthropology researches.

Electronic supplementary material

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

A novel multiplex assay for simultaneously analysing 13 rapidly mutating Y-STRs

A multiplex polymerase chain reaction (PCR) assay (RM-Yplex) was developed which is capable of simultaneously amplifying 13 recently introduced rapidly mutating Y-STR markers (RM Y-STRs). This multiplex assay is expected to aid human identity testing in forensic and other applications to improve differentiating unrelated males and allow separating related males. The 13 RM Y-STR markers included in the multiplex are: DYF387S1, DYF399S1, DYF403S1ab, DYF404S1, DYS449, DYS518, DYS526ab, DYS547, DYS570, DYS576, DYS612, DYS626 and DYS627. This study reflects the proof of concept to analyse all currently known RM Y-STRs simultaneously and describes the optimization of the multiplex assay. The RM-Yplex assay generated complete RM Y-STR profiles down to 62.5pg of male template DNA, and from male-female DNA mixtures at all ratios tested. We herewith introduce and make available for widespread use in forensic and anthropological studies, an effective and sensitive single multiplex assay for simultaneous genotyping of 13 RM Y-STRs.

Inference of sex-specific expansion patterns in human populations from Y-chromosome polymorphism

Studying the current distribution of genetic diversity in humans has important implications for our understanding of the history of our species. We analyzed a set of linked STR and SNP loci from the paternally inherited Y chromosome to infer the past demography of 55 African and Eurasian populations, using both the parametric and nonparametric coalescent-based methods implemented in the BEAST application. We inferred expansion events in most sedentary farmer populations, while we found constant effective population sizes for both nomadic hunter-gatherers and seminomadic herders. Our results differed, on several aspects, from previous results on mtDNA and autosomal markers. First, we found more recent expansion patterns in Eurasia than in Africa. This discrepancy, substantially stronger than the ones found with the other kind of markers, may result from a lower effective population size for men, which might have made male-transmitted markers more sensitive to the out-of-Africa bottleneck. Second, we found expansion signals only for sedentary farmers but not for nomadic herders in Central Asia, while these signals were found for both kind of populations in this area when using mtDNA or autosomal markers. Expansion signals in this area may result from spatial expansion processes and may have been erased for the Y chromosome among the herders because of restricted male gene flow.

24 Y-chromosomal STR haplotypic polymorphisms for Chinese Uygur ethnic group and its phylogenic analysis with other Chinese groups

The Uygur ethnic minority is the largest ethnic group in the Xinjiang Uygur Autonomous Region of China, and is a precious resource for the study of ethnogeny and forensic biology. Previous studies have focused on the genetic background of the Uygur group, however, the patrilineal descent of the group is still unclear. In this study, we investigated the genetic diversity of 24 Y-STR loci in the Uygur group and analyzed the population differentiations as well as the genetic relationships between the Uygur group and other previously reported populations using 17 Y-filer loci. According to haplotypic analysis of the 24 Y-STR loci in 109 Uygur individuals, 104 different haplotypes were obtained, 99 of which were unique. The haplotypic diversity and discrimination capacity of these 24 Y-STR loci in Uygur group were 0.9992 and 0.9541, respectively. An additional 7 loci (DYS388, DYS444, DYS447, DYS449, DYS522, and DYS527a,b) showed high genetic diversity and improved the overall discrimination capacity of the 24 Y-STR system. Pairwise Fst and neighbor-joining analysis showed that the Uygur group was genetically close to the Han populations from different regions.

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

Many studies of human populations have used the male-specific region of the Y chromosome (MSY) as a marker, but MSY sequence variants have traditionally been subject to ascertainment bias. Also, dating of haplogroups has relied on Y-specific short tandem repeats (STRs), involving problems of mutation rate choice, and possible long-term mutation saturation. Next-generation sequencing can ascertain single nucleotide polymorphisms (SNPs) in an unbiased way, leading to phylogenies in which branch-lengths are proportional to time, and allowing the times-to-most-recent-common-ancestor (TMRCAs) of nodes to be estimated directly. Here we describe the sequencing of 3.7 Mb of MSY in each of 448 human males at a mean coverage of 51×, yielding 13,261 high-confidence SNPs, 65.9% of which are previously unreported. The resulting phylogeny covers the majority of the known clades, provides date estimates of nodes, and constitutes a robust evolutionary framework for analyzing the history of other classes of mutation. Different clades within the tree show subtle but significant differences in branch lengths to the root. We also apply a set of 23 Y-STRs to the same samples, allowing SNP- and STR-based diversity and TMRCA estimates to be systematically compared. Ongoing purifying selection is suggested by our analysis of the phylogenetic distribution of nonsynonymous variants in 15 MSY single-copy genes.

Mitochondrial and Y chromosome haplotype motifs as diagnostic markers of Jewish ancestry: a reconsideration

Several authors have proposed haplotype motifs based on site variants at the mitochondrial genome (mtDNA) and the non-recombining portion of the Y chromosome (NRY) to trace the genealogies of Jewish people. Here, we analyzed their main approaches and test the feasibility of adopting motifs as ancestry markers through construction of a large database of mtDNA and NRY haplotypes from public genetic genealogical repositories. We verified the reliability of Jewish ancestry prediction based on the Cohen and Levite Modal Haplotypes in their "classical" 6 STR marker format or in the "extended" 12 STR format, as well as four founder mtDNA lineages (HVS-I segments) accounting for about 40% of the current population of Ashkenazi Jews. For this purpose we compared haplotype composition in individuals of self-reported Jewish ancestry with the rest of European, African or Middle Eastern samples, to test for non-random association of ethno-geographic groups and haplotypes. Overall, NRY and mtDNA based motifs, previously reported to differentiate between groups, were found to be more represented in Jewish compared to non-Jewish groups. However, this seems to stem from common ancestors of Jewish lineages being rather recent respect to ancestors of non-Jewish lineages with the same "haplotype signatures." Moreover, the polyphyly of haplotypes which contain the proposed motifs and the misuse of constant mutation rates heavily affected previous attempts to correctly dating the origin of common ancestries. Accordingly, our results stress the limitations of using the above haplotype motifs as reliable Jewish ancestry predictors and show its inadequacy for forensic or genealogical purposes.

The place of Slovakian paternal diversity in the clinal European landscape

BACKGROUND

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

AIM

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

SUBJECTS AND METHODS

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

RESULTS

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

CONCLUSION

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

Human paternal and maternal demographic histories: insights from high-resolution Y chromosome and mtDNA sequences

Background

Comparisons of maternally-inherited mitochondrial DNA (mtDNA) and paternally-inherited non-recombining Y chromosome (NRY) variation have provided important insights into the impact of sex-biased processes (such as migration, residence pattern, and so on) on human genetic variation. However, such comparisons have been limited by the different molecular methods typically used to assay mtDNA and NRY variation (for example, sequencing hypervariable segments of the control region for mtDNA vs. genotyping SNPs and/or STR loci for the NRY). Here, we report a simple capture array method to enrich Illumina sequencing libraries for approximately 500 kb of NRY sequence, which we use to generate NRY sequences from 623 males from 51 populations in the CEPH Human Genome Diversity Panel (HGDP). We also obtained complete mtDNA genome sequences from the same individuals, allowing us to compare maternal and paternal histories free of any ascertainment bias.

Results

We identified 2,228 SNPs in the NRY sequences and 2,163 SNPs in the mtDNA sequences. Our results confirm the controversial assertion that genetic differences between human populations on a global scale are bigger for the NRY than for mtDNA, although the differences are not as large as previously suggested. More importantly, we find substantial regional variation in patterns of mtDNA versus NRY variation. Model-based simulations indicate very small ancestral effective population sizes (<100) for the out-of-Africa migration as well as for many human populations. We also find that the ratio of female effective population size to male effective population size (N_f/N_m) has been greater than one throughout the history of modern humans, and has recently increased due to faster growth in N_f than N_m.

Conclusions

The NRY and mtDNA sequences provide new insights into the paternal and maternal histories of human populations, and the methods we introduce here should be widely applicable for further such studies.

Inferring population structure and demographic history using Y-STR data from worldwide populations

The Y chromosome is one of the best genetic materials to explore the evolutionary history of human populations. Global analyses of Y chromosomal short tandem repeats (STRs) data can reveal very interesting world population structures and histories. However, previous Y-STR works tended to focus on small geographical ranges or only included limited sample sizes. In this study, we have investigated population structure and demographic history using 17 Y chromosomal STRs data of 979 males from 44 worldwide populations. The largest genetic distances have been observed between pairs of African and non-African populations. American populations with the lowest genetic diversities also showed large genetic distances and coancestry coefficients with other populations, whereas Eurasian populations displayed close genetic affinities. African populations tend to have the oldest time to the most recent common ancestors (TMRCAs), the largest effective population sizes and the earliest expansion times, whereas the American, Siberian, Melanesian, and isolated Atayal populations have the most recent TMRCAs and expansion times, and the smallest effective population sizes. This clear geographic pattern is well consistent with serial founder model for the origin of populations outside Africa. The Y-STR dataset presented here provides the most detailed view of worldwide population structure and human male demographic history, and additionally will be of great benefit to future forensic applications and population genetic studies.

Genetic structure of Qiangic populations residing in the western Sichuan corridor

The Qiangic languages in western Sichuan (WSC) are believed to be the oldest branch of the Sino-Tibetan linguistic family, and therefore, all Sino-Tibetan populations might have originated in WSC. However, very few genetic investigations have been done on Qiangic populations and no genetic evidences for the origin of Sino-Tibetan populations have been provided. By using the informative Y chromosome and mitochondrial DNA (mtDNA) markers, we analyzed the genetic structure of Qiangic populations. Our results revealed a predominantly Northern Asian-specific component in Qiangic populations, especially in maternal lineages. The Qiangic populations are an admixture of the northward migrations of East Asian initial settlers with Y chromosome haplogroup D (D1-M15 and the later originated D3a-P47) in the late Paleolithic age, and the southward Di-Qiang people with dominant haplogroup O3a2c1*-M134 and O3a2c1a-M117 in the Neolithic Age.