Next generation sequencing of Y-STRs in father-son pairs and comparison with traditional capillary electrophoresis.

To evaluate the promising advantages of massively parallel sequencing (MPS) in our casework, we analysed a total of 33 Y-chromosomal short tandem repeats (Y-STRs) with traditional capillary electrophoresis (CE) and 25 Y-STRs using the newer MPS technology. We studied the outcome of both technologies in 64 father-son pairs using stock and custom-designed kits. Current MPS technology confirmed the 13 mutational events observed with CE and improved our understanding of the complex nature of STR mutations. By detecting isometric sequence variants between unrelated males, we show that sequencing Y-STRs using MPS can boost discrimination power.

A Y-chromosomal survey of Ecuador's multi-ethnic population reveals new insights into the tri-partite population structure and supports an early Holocene age of the rare Native American founder lineage C3-MPB373.

Ecuador is a multiethnic and pluricultural country with a complex history defined by migration and admixture processes. The present study aims to increase our knowledge on the Ecuadorian Native Amerindian groups and the unique South American Y-chromosome haplogroup C3-MPB373 through the analysis of up to 23 Y-chromosome STRs (Y-STRs) and several Y-SNPs in a sample of 527 Ecuadorians from 7 distinct populations and geographic areas, including Kichwa and non-Kichwa Native Amerindians, Mestizos and Afro-Ecuadorians. Our results reveal the presence of C3-MPB373 both in the Amazonian lowland Kichwa with frequencies up to 28 % and, for the first time, in notable proportions in Kichwa populations from the Ecuadorian highlands. The substantially higher frequencies of C3-MPB373 in the Amazonian lowlands found in Kichwa and Waorani individuals suggest a founder effect in that area. Notably, estimates for the time to the most recent common ancestor (TMRCA) in the range of 7.2-9.0 kya point to an ancient origin of the haplogroup and suggest an early Holocene expansion of C3-MPB373 into South America. Finally, the pairwise genetic distances (RST) separate the Kichwa Salasaka from all the other Native Amerindian and Ecuadorian groups, indicating a so far hidden diversity among the Kichwa-speaking populations and suggesting a more southern origin of this population. In sum, our study provides a more in-depth knowledge of the male genetic structure of the multiethnic Ecuadorian population, as well as a valuable reference dataset for forensic use.

Genetic assessment reveals no population substructure and divergent regional and sex-specific histories in the Chachapoyas from northeast Peru.

Many native populations in South America have been severely impacted by two relatively recent historical events, the Inca and the Spanish conquest. However decisive these disruptive events may have been, the populations and their gene pools have been shaped markedly also by the history prior to the conquests. This study focuses mainly on the Chachapoya peoples that inhabit the montane forests on the eastern slopes of the northern Peruvian Andes, but also includes three distinct neighboring populations (the Jívaro, the Huancas and the Cajamarca). By assessing mitochondrial, Y-chromosomal and autosomal diversity in the region, we explore questions that have emerged from archaeological and historical studies of the regional culture (s). These studies have shown, among others, that Chachapoyas was a crossroads for Coast-Andes-Amazon interactions since very early times. In this study, we examine the following questions: 1) was there pre-Hispanic genetic population substructure in the Chachapoyas sample? 2) did the Spanish conquest cause a more severe population decline on Chachapoyan males than on females? 3) can we detect different patterns of European gene flow in the Chachapoyas region? and, 4) did the demographic history in the Chachapoyas resemble the one from the Andean area? Despite cultural differences within the Chachapoyas region as shown by archaeological and ethnohistorical research, genetic markers show no significant evidence for past or current population substructure, although an Amazonian gene flow dynamic in the northern part of this territory is suggested. The data also indicates a bottleneck c. 25 generations ago that was more severe among males than females, as well as divergent population histories for populations in the Andean and Amazonian regions. In line with previous studies, we observe high genetic diversity in the Chachapoyas, despite the documented dramatic population declines. The diverse topography and great biodiversity of the northeastern Peruvian montane forests are potential contributing agents in shaping and maintaining the high genetic diversity in the Chachapoyas region.

DNA commission of the International Society of Forensic Genetics (ISFG): Recommendations on the interpretation of Y-STR results in forensic analysis.

Forensic genetic laboratories perform a large amount of STR analyses of the Y chromosome, in particular to analyze the male part of complex DNA mixtures. However, the statistical interpretation of evidence retrieved from Y-STR haplotypes is challenging. Due to the uni-parental inheritance mode, Y-STR loci are connected to each other and thus haplotypes show patterns of relationship on the familial and population level. This precludes the treatment of Y-STR loci as independently inherited variables and the application of the product rule. Instead, the dependency structure of Y-STRs needs to be included in the haplotype frequency estimation process affecting also the current paradigm of a random match probability that is in the autosomal case approximated by the population frequency assuming unrelatedness of sampled individuals. Information on the degree of paternal relatedness in the suspect population as well as on the familial network is however needed to interpret Y-chromosomal results in the best possible way. The previous recommendations of the DNA commission of the ISFG on the use of Y-STRs in forensic analysis published more than a decade ago [1] cover the interpretation issue only marginally. The current recommendations address a number of topics (frequency estimators, databases, metapopulations, LR formulation, triage, rapidly mutating Y-STRs) with relevance for the Y-STR statistics and recommend a decision-based procedure, which takes into account legal requirements as well as availability of population data and statistical methods.

Identification and characterization of novel rapidly mutating Y-chromosomal short tandem repeat markers.

Short tandem repeat polymorphisms on the male-specific part of the human Y-chromosome (Y-STRs) are valuable tools in many areas of human genetics. Although their paternal inheritance and moderate mutation rate (~10-3 mutations per marker per meiosis) allow detecting paternal relationships, they typically fail to separate male relatives. Previously, we identified 13 Y-STR markers with untypically high mutation rates (>10-2 ), termed rapidly mutating (RM) Y-STRs, and showed that they improved male relative differentiation over standard Y-STRs. By applying a newly developed in silico search approach to the Y-chromosome reference sequence, we identified 27 novel RM Y-STR candidates. Genotyping them in 1,616 DNA-confirmed father-son pairs for mutation rate estimation empirically highlighted 12 novel RM Y-STRs. Their capacity to differentiate males related by 1, 2, and 3 meioses was 27%, 47%, and 61%, respectively, while for all 25 currently known RM Y-STRs, it was 44%, 69%, and 83%. Of the 647 Y-STR mutations observed in total, almost all were single repeat changes, repeat gains, and losses were well balanced; allele length and fathers' age were positively correlated with mutation rate. We expect these new RM Y-STRs, together with the previously known ones, to significantly improving male relative differentiation in future human genetic applications.

Inter-laboratory study on standardized MPS libraries: evaluation of performance, concordance, and sensitivity using mixtures and degraded DNA.

We present results from an inter-laboratory massively parallel sequencing (MPS) study in the framework of the SeqForSTRs project to evaluate forensically relevant parameters, such as performance, concordance, and sensitivity, using a standardized sequencing library including reference material, mixtures, and ancient DNA samples. The standardized library was prepared using the ForenSeq DNA Signature Prep Kit (primer mix A). The library was shared between eight European laboratories located in Austria, France, Germany, The Netherlands, and Sweden to perform MPS on their particular MiSeq FGx sequencers. Despite variation in performance between sequencing runs, all laboratories obtained quality metrics that fell within the manufacturer's recommended ranges. Furthermore, differences in locus coverage did not inevitably adversely affect heterozygous balance. Inter-laboratory concordance showed 100% concordant genotypes for the included autosomal and Y-STRs, and still, X-STR concordance exceeded 83%. The exclusive reasons for X-STR discordances were drop-outs at DXS10103. Sensitivity experiments demonstrated that correct allele calling varied between sequencing instruments in particular for lower DNA amounts (≤ 125 pg). The analysis of compromised DNA samples showed the drop-out of one sample (FA10013B01A) while for the remaining three degraded DNA samples MPS was able to successfully type ≥ 87% of all aSTRs, ≥ 78% of all Y-STRs, ≥ 68% of all X-STRs, and ≥ 92% of all iSNPs demonstrating that MPS is a promising tool for human identity testing, which in return, has to undergo rigorous in-house validation before it can be implemented into forensic routine casework.

Y-chromosome diversity of the three major ethno-linguistic groups in the Republic of North Macedonia.

A total of 314 individuals representing the three major ethno-linguistic groups (ethnic Macedonians, Albanians and Turks) in the Republic of North Macedonia were analyzed for Y-SNPs and Y-STRs using minisequencing and fragment analysis. The haplogroup composition differed remarkably between the three groups with dominance of haplogroup I2 in ethnic Macedonians (28.1%), E1b in Albanians (35.3%) and J2a (34.9%) in Turks, respectively. The haplotype analysis using the YFilerPlus kit disclosed a significant reduction in diversity values (DC, GD) for the Turkish subgroup compared to the Macedonian and Albanian speaking populations. The Y-STR based population analysis revealed a similarity of ethnic Macedonians with neighboring Serbians and Bulgarians. The same holds true for the Albanian speakers from Macedonia and Albania, whereas the Turkish minority in North Macedonia stands apart from the population in Turkey.

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.

Genetic structure and sex-biased gene flow in the history of southern African populations.

OBJECTIVES: We investigated the genetic history of southern African populations with a special focus on their paternal history. We reexamined previous claims that the Y-chromosome haplogroup E1b1b (E-M293) was brought to southern Africa by pastoralists from eastern Africa, and investigated patterns of sex-biased gene flow in southern Africa. MATERIALS AND METHODS: We analyzed previously published complete mtDNA genome sequences and ∼900 kb of NRY sequences from 23 populations from Namibia, Botswana, and Zambia, as well as haplogroup frequencies from a large sample of southern African populations and 23 newly genotyped Y-linked STR loci for samples assigned to haplogroup E1b1b. RESULTS: Our results support an eastern African origin for Y-chromosome haplogroup E1b1b (E-M293); however, its current distribution in southern Africa is not strongly associated with pastoralism, suggesting more complex demographic events and/or changes in subsistence practices in this region. The Bantu expansion in southern Africa had a notable genetic impact and was probably a rapid, male-dominated expansion. Our finding of a significant increase in the intensity of the sex-biased gene flow from north to south may reflect changes in the social dynamics between Khoisan and Bantu groups over time. CONCLUSIONS: Our study shows that the population history of southern Africa has been complex, with different immigrating groups mixing to different degrees with the autochthonous populations. The Bantu expansion led to heavily sex-biased admixture as a result of interactions between Khoisan females and Bantu males, with a geographic gradient which may reflect changes in the social dynamics between Khoisan and Bantu groups over time.

Inter-laboratory validation study of the ForenSeq™ DNA Signature Prep Kit.

The implementation of massively parallel sequencing (MPS) in forensic science revealed the advantages of the new method compared to the forensic benchmark in DNA-STR analysis, the capillary-electrophoresis (CE): Sequence information and the possibility to multiplex hundreds of markers in one multiplex PCR increase the discrimination power of a forensic (STR-) profile. The EU funded project DNASeqEx (DNA-STR Massive Sequencing & International Information Exchange) aims to evaluate MPS-based materials in their respective developmental stages using the two established platforms MiSeq FGx (Illumina) and Ion S5™ (Thermo Fisher Scientific). As part of this project, we present here an inter-laboratory validation of the Forenseq™ DNA Signature Prep Kit, focussing on STRs included in primer mix A. Our study comprises tests of concordance, reproducibility, sensitivity (1 ng, 500 pg, 250 pg, 125 pg, 63 pg, 31 pg) and mixtures (male-male and male-female at ratios of 1:1, 1:5, 1:10, 1:15, 1:20, 1:100, 1:500, 1:1000). Sequencing results found to be virtually concordant to CE results, to reference profiles and reproducible between duplicates and between both laboratories. We observed first locus drop-outs (LDO) at a DNA input of 63 pg (20 sample pool) and 125 pg (38 sample pool). Alleles were found to be well balanced at a DNA input of 250 pg or more. We found the kit to perform well on moderate mixtures (1:1-1:20).

Current state-of-art of STR sequencing in forensic genetics.

The current state of validation and implementation strategies of massively parallel sequencing (MPS) technology for the analysis of STR markers for forensic genetics use is described, covering the topics of the current catalog of commercial MPS-STR panels, leading MPS-platforms, and MPS-STR data analysis tools. In addition, the developmental and internal validation studies carried out to date to evaluate reliability, sensitivity, mixture analysis, concordance, and the ability to analyze challenged samples are summarized. The results of various MPS-STR population studies that showed a large number of new STR sequence variants that increase the power of discrimination in several forensically relevant loci are also presented. Finally, various initiatives developed by several international projects and standardization (or guidelines) groups to facilitate application of MPS technology for STR marker analyses are discussed in regard to promoting a standard STR sequence nomenclature, performing population studies to detect sequence variants, and developing a universal system to translate sequence variants into a simple STR nomenclature (numbers and letters) compatible with national STR databases.

Enclaves of genetic diversity resisted Inca impacts on population history.

The Inca Empire is claimed to have driven massive population movements in western South America, and to have spread Quechua, the most widely-spoken language family of the indigenous Americas. A test-case is the Chachapoyas region of northern Peru, reported as a focal point of Inca population displacements. Chachapoyas also spans the environmental, cultural and demographic divides between Amazonia and the Andes, and stands along the lowest-altitude corridor from the rainforest to the Pacific coast. Following a sampling strategy informed by linguistic data, we collected 119 samples, analysed for full mtDNA genomes and Y-chromosome STRs. We report a high indigenous component, which stands apart from the network of intense genetic exchange in the core central zone of Andean civilization, and is also distinct from neighbouring populations. This unique genetic profile challenges the routine assumption of large-scale population relocations by the Incas. Furthermore, speakers of Chachapoyas Quechua are found to share no particular genetic similarity or gene-flow with Quechua speakers elsewhere, suggesting that here the language spread primarily by cultural diffusion, not migration. Our results demonstrate how population genetics, when fully guided by the archaeological, historical and linguistic records, can inform multiple disciplines within anthropology.

Revisiting the male genetic landscape of China: a multi-center study of almost 38,000 Y-STR haplotypes.

China has repeatedly been the subject of genetic studies to elucidate its prehistoric and historic demography. While some studies reported a genetic distinction between Northern and Southern Han Chinese, others showed a more clinal picture of small differences within China. Here, we investigated the distribution of Y chromosome variation along administrative as well as ethnic divisions in the mainland territory of the People's Republic of China, including 28 administrative regions and 19 recognized Chinese nationalities, to assess the impact of recent demographic processes. To this end, we analyzed 37,994 Y chromosomal 17-marker haplotype profiles from the YHRD database with respect to forensic diversity measures and genetic distance between groups defined by administrative boundaries and ethnic origin. We observed high diversity throughout all Chinese provinces and ethnicities. Some ethnicities, including most prominently Kazakhs and Tibetans, showed significant genetic differentiation from the Han and other groups. However, differences between provinces were, except for those located on the Tibetan plateau, less pronounced. This discrepancy is explicable by the sizeable presence of Han speakers, who showed high genetic homogeneity all across China, in nearly all studied provinces. Furthermore, we observed a continuous genetic North-South gradient in the Han, confirming previous reports of a clinal distribution of Y chromosome variation and being in notable concordance with the previously observed spatial distribution of autosomal variation. Our findings shed light on the demographic changes in China accrued by a fast-growing and increasingly mobile population.

Male Lineages in Brazil: Intercontinental Admixture and Stratification of the European Background.

The non-recombining nature of the Y chromosome and the well-established phylogeny of Y-specific Single Nucleotide Polymorphisms (Y-SNPs) make them useful for defining haplogroups with high geographical specificity; therefore, they are more apt than the Y-STRs to detect population stratification in admixed populations from diverse continental origins. Different Y-SNP typing strategies have been described to address issues of population history and movements within geographic territories of interest. In this study, we investigated a set of 41 Y-SNPs in 1217 unrelated males from the five Brazilian geopolitical regions, aiming to disclose the genetic structure of male lineages in the country. A population comparison based on pairwise FST genetic distances did not reveal statistically significant differences in haplogroup frequency distributions among populations from the different regions. The genetic differences observed among regions were, however, consistent with the colonization history of the country. The sample from the Northern region presented the highest Native American ancestry (8.4%), whereas the more pronounced African contribution could be observed in the Northeastern population (15.1%). The Central-Western and Southern samples showed the higher European contributions (95.7% and 93.6%, respectively). The Southeastern region presented significant European (86.1%) and African (12.0%) contributions. The subtyping of the most frequent European lineage in Brazil (R1b1a-M269) allowed differences in the genetic European background of the five Brazilian regions to be investigated for the first time.

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.

The complete mitogenome of a 500-year-old Inca child mummy.

In 1985, a frozen mummy was found in Cerro Aconcagua (Argentina). Archaeological studies identified the mummy as a seven-year-old Inca sacrifice victim who lived >500 years ago, at the time of the expansion of the Inca Empire towards the southern cone. The sequence of its entire mitogenome was obtained. After querying a large worldwide database of mitogenomes (>28,000) we found that the Inca haplotype belonged to a branch of haplogroup C1b (C1bi) that has not yet been identified in modern Native Americans. The expansion of C1b into the Americas, as estimated using 203 C1b mitogenomes, dates to the initial Paleoindian settlements (~18.3 thousand years ago [kya]); however, its internal variation differs between Mesoamerica and South America. By querying large databases of control region haplotypes (>150,000), we found only a few C1bi members in Peru and Bolivia (e.g. Aymaras), including one haplotype retrieved from ancient DNA of an individual belonging to the Wari Empire (Peruvian Andes). Overall, the results suggest that the profile of the mummy represents a very rare sub-clade that arose 14.3 (5-23.6) kya and could have been more frequent in the past. A Peruvian Inca origin for present-day C1bi haplotypes would satisfy both the genetic and paleo-anthropological findings.

Validation of a combined autosomal/Y-chromosomal STR approach for analyzing typical biological stains in sexual-assault cases.

DNA testing is an established part of the investigation and prosecution of sexual assault. The primary purpose of DNA evidence is to identify a suspect and/or to demonstrate sexual contact. However, due to highly uneven proportions of female and male DNA in typical stains, routine autosomal analysis often fails to detect the DNA of the assailant. To evaluate the forensic efficiency of the combined application of autosomal and Y-chromosomal short tandem repeat (STR) markers, we present a large retrospective casework study of probative evidence collected in sexual-assault cases. We investigated up to 39 STR markers by testing combinations of the 16-locus NGMSElect kit with both the 23-locus PowerPlex Y23 and the 17-locus Yfiler kit. Using this dual approach we analyzed DNA extracts from 2077 biological stains collected in 287 cases over 30 months. To assess the outcome of the combined approach in comparison to stand-alone autosomal analysis we evaluated informative DNA profiles. Our investigation revealed that Y-STR analysis added up to 21% additional, highly informative (complete, single-source) profiles to the set of reportable autosomal STR profiles for typical stains collected in sexual-assault cases. Detection of multiple male contributors was approximately three times more likely with Y-chromosomal profiling than with autosomal STR profiling. In summary, 1/10 cases would have remained inconclusive (and could have been dismissed) if Y-STR analysis had been omitted from DNA profiling in sexual-assault cases.