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.

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.

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.

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.

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.

Massively parallel sequencing of forensic STRs: Considerations of the DNA commission of the International Society for Forensic Genetics (ISFG) on minimal nomenclature requirements.

The DNA Commission of the International Society for Forensic Genetics (ISFG) is reviewing factors that need to be considered ahead of the adoption by the forensic community of short tandem repeat (STR) genotyping by massively parallel sequencing (MPS) technologies. MPS produces sequence data that provide a precise description of the repeat allele structure of a STR marker and variants that may reside in the flanking areas of the repeat region. When a STR contains a complex arrangement of repeat motifs, the level of genetic polymorphism revealed by the sequence data can increase substantially. As repeat structures can be complex and include substitutions, insertions, deletions, variable tandem repeat arrangements of multiple nucleotide motifs, and flanking region SNPs, established capillary electrophoresis (CE) allele descriptions must be supplemented by a new system of STR allele nomenclature, which retains backward compatibility with the CE data that currently populate national DNA databases and that will continue to be produced for the coming years. Thus, there is a pressing need to produce a standardized framework for describing complex sequences that enable comparison with currently used repeat allele nomenclature derived from conventional CE systems. It is important to discern three levels of information in hierarchical order (i) the sequence, (ii) the alignment, and (iii) the nomenclature of STR sequence data. We propose a sequence (text) string format the minimal requirement of data storage that laboratories should follow when adopting MPS of STRs. We further discuss the variant annotation and sequence comparison framework necessary to maintain compatibility among established and future data. This system must be easy to use and interpret by the DNA specialist, based on a universally accessible genome assembly, and in place before the uptake of MPS by the general forensic community starts to generate sequence data on a large scale. While the established nomenclature for CE-based STR analysis will remain unchanged in the future, the nomenclature of sequence-based STR genotypes will need to follow updated rules and be generated by expert systems that translate MPS sequences to match CE conventions in order to guarantee compatibility between the different generations of STR data.

The new Y Chromosome Haplotype Reference Database.

After opening the first version of an internet-accessible worldwide reference database of Y chromosome profiles 14 years ago and six years after the last major relaunch the new YHRD 4.0 repository and website has been rolled-out. By November 2014 about 136k 9-locus haplotypes, among these 84k 17-locus haplotypes, 25k 23-locus haplotypes and 15k Y SNP profiles from 917 sampling locations in 128 countries have been submitted by more than 250 institutes and laboratories. In geographic terms, about 39% of the YHRD samples are from Europe, 32% from Asia, 16% from South America, 6% from North America, 4% from Africa and 2% from Oceania/Australia. Worldwide collaboration is the driving force for the rapid growth of the database and this, in turn, allows the evaluation and implementation of enhanced interpretation tools (variable frequency estimators, LR-based mixture and kinship analysis, Y-SNP-based ancestry assessment).

No shortcut solution to the problem of Y-STR match probability calculation.

Match probability calculation is deemed much more intricate for lineage genetic markers, including Y-chromosomal short tandem repeats (Y-STRs), than for autosomal markers. This is because, owing to the lack of recombination, strong interdependence between markers is likely, which implies that haplotype frequency estimates cannot simply be obtained through the multiplication of allele frequency estimates. As yet, however, the practical relevance of this problem has not been studied in much detail using real data. In fact, such scrutiny appears well warranted because the high mutation rates of Y-STRs and the possibility of backward mutation should have worked against the statistical association of Y-STRs. We examined haplotype data of 21 markers included in the PowerPlex(®)Y23 set (PPY23, Promega Corporation, Madison, WI) originating from six different populations (four European and two Asian). Assessing the conditional entropies of the markers, given different subsets of markers from the same panel, we demonstrate that the PowerPlex(®)Y23 set cannot be decomposed into smaller marker subsets that would be (conditionally) independent. Nevertheless, in all six populations, >94% of the joint entropy of the 21 markers is explained by the seven most rapidly mutating markers. Although this result might render a reduction in marker number a sensible option for practical casework, the partial haplotypes would still be almost as diverse as the full haplotypes. Therefore, match probability calculation remains difficult and calls for the improvement of currently available methods of haplotype frequency estimation.

A global analysis of Y-chromosomal haplotype diversity for 23 STR loci.

In a worldwide collaborative effort, 19,630 Y-chromosomes were sampled from 129 different populations in 51 countries. These chromosomes were typed for 23 short-tandem repeat (STR) loci (DYS19, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS385ab, DYS437, DYS438, DYS439, DYS448, DYS456, DYS458, DYS635, GATAH4, DYS481, DYS533, DYS549, DYS570, DYS576, and DYS643) and using the PowerPlex Y23 System (PPY23, Promega Corporation, Madison, WI). Locus-specific allelic spectra of these markers were determined and a consistently high level of allelic diversity was observed. A considerable number of null, duplicate and off-ladder alleles were revealed. Standard single-locus and haplotype-based parameters were calculated and compared between subsets of Y-STR markers established for forensic casework. The PPY23 marker set provides substantially stronger discriminatory power than other available kits but at the same time reveals the same general patterns of population structure as other marker sets. A strong correlation was observed between the number of Y-STRs included in a marker set and some of the forensic parameters under study. Interestingly a weak but consistent trend toward smaller genetic distances resulting from larger numbers of markers became apparent.

Y-chromosomal analysis identifies the skeletal remains of Swiss national hero Jörg Jenatsch (1596-1639).

Jörg Jenatsch was a Swiss defender of independence and a fighter for liberty in the 17th century. With the help of three living male members of the Jenatsch family, we successfully identified a skeleton exhumed from Chur cathedral as the remains of Jörg Jenatsch. Our conclusion was based upon complete Y-STR and Y-SNP profiles that could be generated by replicate analyses of a bone sample available to us. The skeleton and the three living family members carried the same Y-SNP haplogroup, but were discordant at three of 23 Y-STR loci. This notwithstanding, conservative biostatistical evaluation of the data suggests that the Chur skeleton is at least 20 times more likely than not to be Jörg Jenatsch.

Multiple recurrent mutations at four human Y-chromosomal single nucleotide polymorphism sites in a 37 bp sequence tract on the ARSDP1 pseudogene.

The male-specific region of the human Y chromosome (MSY) is passed down clonally from father to son and mutation is the single driving force for Y-chromosomal diversification. The geographical distribution of MSY variation is non-random. Therefore, Y-chromosomal single nucleotide polymorphisms (Y-SNPs) are of forensic interest, as they can be utilized, e.g. for deducing the bio-geographical origin of biological evidence. This extra information can complement short tandem repeat data in criminal investigations. For forensic applications, however, any targeted marker has to be unequivocally interpretable. Here, we report findings for 17 samples from a population study comprising specimens from ∼3700 men living in Tyrol (Austria), indicating apparent homoplasic mutations at four Y-SNP loci on haplogroup R-M412/L51/S167, R-U152/S28, and L-M20 Y chromosomes. The affected Y-SNPs P41, P37, L202, and L203 mapped to a 37bp region on Yq11.21. Observing in multiple phylogenetic contexts up to four homoplasic mutations within such a short sequence tract is unlikely to result from a series of independent parallel mutations. Hence, we rather propose X-to-Y gene conversion as a more likely scenario. Practical implications arising from markers exhibiting paralogues on the Y chromosome or sites with a high propensity to recurrent mutation for database searches are addressed.

Continent-wide decoupling of Y-chromosomal genetic variation from language and geography in native South Americans.

Numerous studies of human populations in Europe and Asia have revealed a concordance between their extant genetic structure and the prevailing regional pattern of geography and language. For native South Americans, however, such evidence has been lacking so far. Therefore, we examined the relationship between Y-chromosomal genotype on the one hand, and male geographic origin and linguistic affiliation on the other, in the largest study of South American natives to date in terms of sampled individuals and populations. A total of 1,011 individuals, representing 50 tribal populations from 81 settlements, were genotyped for up to 17 short tandem repeat (STR) markers and 16 single nucleotide polymorphisms (Y-SNPs), the latter resolving phylogenetic lineages Q and C. Virtually no structure became apparent for the extant Y-chromosomal genetic variation of South American males that could sensibly be related to their inter-tribal geographic and linguistic relationships. This continent-wide decoupling is consistent with a rapid peopling of the continent followed by long periods of isolation in small groups. Furthermore, for the first time, we identified a distinct geographical cluster of Y-SNP lineages C-M217 (C3*) in South America. Such haplotypes are virtually absent from North and Central America, but occur at high frequency in Asia. Together with the locally confined Y-STR autocorrelation observed in our study as a whole, the available data therefore suggest a late introduction of C3* into South America no more than 6,000 years ago, perhaps via coastal or trans-Pacific routes. Extensive simulations revealed that the observed lack of haplogroup C3* among extant North and Central American natives is only compatible with low levels of migration between the ancestor populations of C3* carriers and non-carriers. In summary, our data highlight the fact that a pronounced correlation between genetic and geographic/cultural structure can only be expected under very specific conditions, most of which are likely not to have been met by the ancestors of native South Americans.

Estimating trace-suspect match probabilities for singleton Y-STR haplotypes using coalescent theory.

Estimation of match probabilities for singleton haplotypes of lineage markers, i.e. for haplotypes observed only once in a reference database augmented by a suspect profile, is an important problem in forensic genetics. We compared the performance of four estimators of singleton match probabilities for Y-STRs, namely the count estimate, both with and without Brenner's so-called 'kappa correction', the surveying estimate, and a previously proposed, but rarely used, coalescent-based approach implemented in the BATWING software. Extensive simulation with BATWING of the underlying population history, haplotype evolution and subsequent database sampling revealed that the coalescent-based approach is characterized by lower bias and lower mean squared error than the uncorrected count estimator and the surveying estimator. Moreover, in contrast to the two count estimators, both the surveying and the coalescent-based approach exhibited a good correlation between the estimated and true match probabilities. However, although its overall performance is thus better than that of any other recognized method, the coalescent-based estimator is still computation-intense on the verge of general impracticability. Its application in forensic practice therefore will have to be limited to small reference databases, or to isolated cases of particular interest, until more powerful algorithms for coalescent simulation have become available.

Y-chromosomal STR analysis in the Pashtun population of Southern Afghanistan.

Afghanistan is a landlocked country in the heart of Asia and since the dawn of humankind Afghanistan has faced centuries of turmoil, strife, conflict, warfare, distress, social unrest, difficult climate, harsh terrain and due to its unique geostrategic position in Eurasia which has historically attracted commerce and conflict. It is an important stop along the Silk Road, connecting the far eastern civilizations to the western world. A 5000-year history of constant invasion. Afghanistan has been repeatedly invaded and conquered by rulers and super powers, neighboring interference in this conflict-tattered land for centuries yet rarely leading to the conquest of this rugged and challenging terrain nation. Afghans are not only shepherds, farmers and nomads but also intense fighters and fierce warriors. Currently very limited genetic studies have been performed in Afghan populations. 17 Y chromosomal short tandem repeats (Y-STRs) were analyzed in 125 unrelated Pashtun (in hindi: Pathan) males residing in the Kandahar region of Southern Afghanistan. A total of 92 unique haplotypes were observed. The predominant haplotype reached a frequency of 9.6%. The haplotype diversity was 0.987 and the discrimination capacity 73.6%. Analysis of molecular variance (AMOVA) reveals a considerable regional stratification within the country as well as between different Pashtun (Pathan) groups from Afghanistan, Pakistan and India.

Y-STR Frequency Surveying Method: A critical reappraisal.

Reasonable formalized methods to estimate the frequencies of DNA profiles generated from lineage markers have been proposed in the past years and were discussed in the forensic community. Recently, collections of population data on the frequencies of variations in Y chromosomal STR profiles have reached a new quality with the establishment of the comprehensive neatly quality-controlled reference database YHRD. Grounded on such unrivalled empirical material from hundreds of populations studies the core assumption of the Haplotype Frequency Surveying Method originally described 10 years ago can be tested and improved. Here we provide new approaches to calculate the parameters used in the frequency surveying method: a maximum likelihood estimation of the regression parameters (r(1), r(2), s(1) and s(2)) and a revised Frequency Surveying framework with variable binning and a database preprocessing to take the population sub-structure into account. We found good estimates for 11 metapopulations using both approaches and demonstrate that the statistical basis of the method is well supported and independent of the population under study. The results of the estimation process are reliable and robust if the underlying datasets are large and representative and show small average and pairwise genetic distances.

Geostatistical inference of main Y-STR-haplotype groups in Europe.

We examined the multifarious genetic heterogeneity of Europe and neighboring regions from a geographical perspective. We created composite maps outlining the estimated geographical distribution of major groups of genetically similar individuals on the basis of forensic Y-chromosomal markers. We analyzed Y-chromosomal haplotypes composed of 7 highly polymorphic STR loci, genotyped for 33,010 samples, collected at 249 sites in Europe, Western Asia and North Africa, deposited in the YHRD database (www.yhrd.org). The data set comprised 4176 different haplotypes, which we grouped into 20 clusters. For each cluster, the frequency per site was calculated. All geostatistical analysis was performed with the geographic information system GRASS-GIS. We interpolated frequency values across the study area separately for each cluster. Juxtaposing all 20 interpolated surfaces, we point-wisely screened for the highest cluster frequencies and stored it in parallel with the respective cluster label. We combined these two types of data in a composite map. We repeated this procedure for the second highest frequencies in Europe. Major groups were assigned to Northern, Western and Eastern Europe. North Africa built a separate region, Southeastern Europe, Turkey and Near East were divided into several regions. The spatial distribution of the groups accounting for the second highest frequencies in Europe overlapped with the territories of the largest countries. The genetic structure presented in the composite maps fits major historical geopolitical regions and is in agreement with previous studies of genetic frequencies, validating our approach. Our genetic geostatistical approach provides, on the basis of two composite maps, detailed evidence of the geographical distribution and relative frequencies of the most predominant groups of the extant male European population, examined on the basis of forensic Y-STR haplotypes. The existence of considerable genetic differences among geographic subgroups in Europe has important consequences for the statistical inference in forensic Y-STR haplotype analyses.