Comprehensive mutation analysis of 17 Y-chromosomal short tandem repeat polymorphisms included in the AmpFlSTR Yfiler PCR amplification kit

Bibliometric analysis of kinship analysis from 1960 to 2023: global trends and development

Kinship analysis is a crucial aspect of forensic genetics. This study analyzed 1,222 publications on kinship analysis from 1960 to 2023 using bibliometric analysis techniques, investigating the annual publication and citation patterns, most productive countries, organizations, authors and journals, most cited documents and co-occurrence of keywords. The initial publication in this field occurred in 1960. Since 2007, there has been a significant increase in publications, with over 30 published annually except for 2010. China had the most publications (n = 213, 17.43%), followed by the United States (n = 175, 14.32%) and Germany (n = 89, 7.28%). The United States also had the highest citation count. Sichuan University in China has the largest number of published articles. The University of Leipzig and the University of Cologne in Germany exhibit the highest total citation count and average citation, respectively. Budowle B was the most prolific author and Kayser M was the most cited author. In terms of publications, Forensic Science International- Genetics, Forensic Science International, and International Journal of Legal Medicine were the most prolific journals. Among them, Forensic Science International-Genetics boasted the highest h-index, citation count, and average citation rate. The most frequently cited publication was "Van Oven M, 2009, Hum Mutat", with a total of 1,361 citations. The most frequent co-occurrence keyword included "DNA", "Loci", "Paternity testing", "Population", "Markers", and "Identification", with recent interest focusing on "Kinship analysis", "SNP" and "Inference". The current research is centered around microhaplotypes, forensic genetic genealogy, and massively parallel sequencing. The field advanced with new DNA analysis methods, tools, and genetic markers. Collaborative research among nations, organizations, and authors benefits idea exchange, problem-solving efficiency, and high-quality results.

Microsatellites' mutation modeling through the analysis of the Y-chromosomal transmission: Results of a GHEP-ISFG collaborative study

The Spanish and Portuguese Speaking Working Group of the International Society for Forensic Genetics (GHEP-ISFG) organized a collaborative study on mutations of Y-chromosomal short tandem repeats (Y-STRs). New data from 2225 father-son duos and data from 44 previously published reports, corresponding to 25,729 duos, were collected and analyzed. Marker-specific mutation rates were estimated for 33 Y-STRs. Although highly dependent on the analyzed marker, mutations compatible with the gain or loss of a single repeat were 23.2 times more likely than those involving a greater number of repeats. Longer alleles (relatively to the modal one) showed to be nearly twice more mutable than the shorter ones. Within the subset of longer alleles, the loss of repeats showed to be nearly twice more likely than the gain. Conversely, shorter alleles showed a symmetrical trend, with repeat gains being twofold more frequent than reductions. A positive correlation between the paternal age and the mutation rate was observed, strengthening previous findings. The results of a machine learning approach, via logistic regression analyses, allowed the establishment of algebraic formulas for estimating the probability of mutation depending on paternal age and allele length for DYS389I, DYS393 and DYS627. Algebraic formulas could also be established considering only the allele length as predictor for DYS19, DYS389I, DYS389II-I, DYS390, DYS391, DYS393, DYS437, DYS439, DYS449, DYS456, DYS458, DYS460, DYS481, DYS518, DYS533, DYS576, DYS626 and DYS627 loci. For the remaining Y-STRs, a lack of statistical significance was observed, probably as a consequence of the small effective size of the subsets available, a common difficulty in the modeling of rare events as is the case of mutations. The amount of data used in the different analyses varied widely, depending on how the data were reported in the publications analyzed. This shows a regrettable waste of produced data, due to inadequate communication of the results, supporting an urgent need of publication guidelines for mutation studies.

Improving the strategy to identify historical military remains: a literature review and Y-STR meta-analysis

The identification of historical military remains by Unrecovered War Casualties-Army (UWC-A) currently relies on Y-chromosome Short Tandem Repeat (Y-STR) testing when maternal relatives are not available, or when a mitochondrial DNA match does not provide sufficient certainty of identification. However, common Y-STR profiles (using Yfiler™) between sets of remains or families often prevent identification. To resolve these cases, an investigation of additional Y-DNA markers is needed for their potential inclusion into the DNA identification strategy. The number of genetic transmissions between missing soldiers and their living relatives needs to be considered to avoid false exclusions between paternal relatives. Analysis of 236 World War I/II (WWI/II) era pairs of relatives identified up to seven genetic transmissions between WWII soldiers and their living relatives, and nine for WWI. Previous Y-STR meta-analyses were published approximately 10 years ago when rapidly mutating markers were relatively new. This paper reports a contemporary literature review and meta-analysis of 35 studies (which includes 23 studies not previously used in meta-analysis) and 23 commonly used Y-STR's mutation rates to inform the inclusion of additional loci to UWC-A's DNA identification strategy. Meta-analysis found mutation data for a given Y-STR locus could be pooled between studies and that the mutation rates were significantly different between some loci (at P < 0.05). Based on this meta-analysis, we have identified two additional markers from PowerPlex® Y23 for potential inclusion in UWC-A's identification strategy. Further avenues for potential experimental exploration are discussed.

Key points

From 236 UWC-A pairs of relatives, we observed up to nine genetic transmissions between WWI soldiers and their living relatives, and seven for WWII.MedCalc® software for meta-analysis utilizing the Freeman-Tukey transformation was run, which analysed 35 published studies and 23 commonly used loci. Previous Y-STR mutation rate meta-analyses are now 10 years old; this paper includes 23 studies that were not included in previous meta-analyses.Through meta-analysis, we identify two markers from PowerPlex® Y23 for potential inclusion in UWC-A's historical remains identification strategy (alongside Yfiler™). We discuss potential next steps for experimental exploration of additional Y-DNA markers.

Delineating the dispersal of Y-chromosome sub-haplogroup O2a2b-P164 among Austronesian-speaking populations

This article reports on an exploration of the Y-chromosome sub-haplogroup O2a2b-P164 in Austronesian-speaking populations. Moderate to high abundance of the P 164 mutation is seen in the West Pacific including the Amis of Formosa (36%) and the Filipinos of Mindanao (50%) as well as in the Kiritimati of Micronesia (70%), and Tonga and Samoa of West Polynesia (54% and 33%, respectively), and it drops to low frequencies in populations of East Polynesia. The communities of Polynesia and Micronesia exhibit considerable inter- and intra-population haplotype sharing suggesting extensive population affinity. The observed affinities, as well as the ages and diversity values within the P 164 sub-haplogroup among Austronesian-speaking populations signal an ancestral migration route and relationships that link the Amis of Taiwan with distant communities in West and East Polynesia, Micronesia, and the Maori of New Zealand. High resolution sequencing of the Austronesian Y chromosome indicate that the P 164 lineage originated about 19,000 ya and then split into three branches separating the Ami aborigines, Southeast Asian and Polynesian/Micronesian populations about 4700 ya, roughly coinciding with the initiation of the Austronesian diaspora. The Y-chromosomes of all the Polynesian and Micronesian population examined belong to the new FT 257096 haplogroup.

A posteriori parameters from paternity tests of a Mexican laboratory with the powerplex fusion system

Population studies regarding Human identification (HID) systems report a priori forensic parameters, but rarely they describe a posteriori parameters from concluded paternity tests. We analyzed data from the PowerPlex® Fusion System in 1503 paternity tests from a Mexican laboratory for five years (2016-2020). The motherless duo paternity tests (89.8%) were more frequent than the standard trio tests (10.2%). A notable increase in motherless tests was noted regarding our previous report (89.8% vs 77.3%), probably explained by the COVID-19 pandemic. The estimated exclusion frequency in Mexico ranged from 30.1 (trio) to 32.1% (duo). For paternity exclusions, we report the number of mismatches and the frequency at which each STR was involved. The PowerPlex® Fusion system showed more than five mismatches in 100% of the standard trio tests excluding paternity, and the majority of motherless-duo tests (98.1%). In positive paternity tests, PowerPlex® Fusion offered a higher combined paternity index (PI) (average 1.18 E + 10) regarding HID systems with 15 and 20 STRs, even without the inclusion of the Y-linked locus DYS391 to the kinship interpretation. Individual and global STR mutation rates were estimated from 17 paternal mutations (μ = 0.0017), the majority involving a single-step mutation (94.11%). Five independent null alleles were detected, most of them involving the Penta E locus (80%), which suggests caution to the users working with DNA databases or kinship analysis, to avoid false exclusions with Penta E. In brief, our results provide a better overview of a posteriori informativeness offered by the PowerPlex® Fusion system for paternity testing in Mexico.

The sequence of the repetitive motif influences the frequency of multistep mutations in Short Tandem Repeats

Microsatellites, or Short Tandem Repeats (STRs), are subject to frequent length mutations that involve the loss or gain of an integer number of repeats. This work aimed to investigate the correlation between STRs' specific repetitive motif composition and mutational dynamics, specifically the occurrence of single- or multistep mutations. Allelic transmission data, comprising 323,818 allele transfers and 1,297 mutations, were gathered for 35 Y-chromosomal STRs with simple structure. Six structure groups were established: ATT, CTT, TCTA/GATA, GAAA/CTTT, CTTTT, and AGAGAT, according to the repetitive motif present in the DNA leading strand of the markers. Results show that the occurrence of multistep mutations varies significantly among groups of markers defined by the repetitive motif. The group of markers with the highest frequency of multistep mutations was the one with repetitive motif CTTTT (25% of the detected mutations) and the lowest frequency corresponding to the group with repetitive motifs TCTA/GATA (0.93%). Statistically significant differences (α = 0.05) were found between groups with repetitive motifs with different lengths, as is the case of TCTA/GATA and ATT (p = 0.0168), CTT (p < 0.0001) and CTTTT (p < 0.0001), as well as between GAAA/CTTT and CTTTT (p = 0.0102). The same occurred between the two tetrameric groups GAAA/CTTT and TCTA/GATA (p < 0.0001) - the first showing 5.7 times more multistep mutations than the second. When considering the number of repeats of the mutated paternal alleles, statistically significant differences were found for alleles with 10 or 12 repeats, between GATA and ATT structure groups. These results, which demonstrate the heterogeneity of mutational dynamics across repeat motifs, have implications in the fields of population genetics, epidemiology, or phylogeography, and whenever STR mutation models are used in evolutionary studies in general.

Bidirectional dispersals during the peopling of the North American Arctic

It is thought that Paleo-Inuit were the first people that settled the American Arctic about 5000 BP (before the present) from a migration that crossed Beringia from Northeast Asia. It is theorized that this group initially migrated to the North Slopes of Alaska and subsequently expanded eastward, eventually reaching Greenland. A second circumpolar dispersal of Neo-Inuit from the North Slopes associated with the Thule-Inuk culture has been postulated to have extended eastward around 800 BP, totally replacing the original Paleo-Inuit without admixing. Although generally accepted, this migration scenario is incompatible with previously reported indications of east to west gene flow across the American Arctic. Here we report on the Y-chromosome haplogroup and Y-STR diversity of the four circumpolar populations of the Tuva Republic (N = 24), Northeast Siberia (N = 9), Bethel, Alaska (N = 40), and Barrow, Alaska (N = 31). Four haplogroup lineages (Q-NWT01, Q-M3, Q-M346, and Q-M120) were detected, Q-NWT01 and Q-M3 being the most abundant at 11.11 and 66.67% in Northeast Siberia, 32.50 and 65.00% in Bethel, and 67.74 and 32.26% in Barrow, respectively. The same samples genotyped for Y-chromosome SNPs were typed for 17 Y-STYR loci using the AmpFlSTR Yfiler system. Age estimates and diversity values for the Q-NWT01 and Q-M3 mutations suggest extensive movement of male individuals along the entire longitudinal stretch of the American circumpolar region. Throughout the entire region, Q-M3 exhibits a west to east decreasing gradient in age and diversity while Q-NWT01 indicates the opposite with older TMRCA and higher diversity values running from east to west with the most recent estimates in Canada and Alaska. The high age and diversity values in Greenland are congruent with an origin of the Q-NWT01 mutation in the east of the circumpolar range about 2000-3000 ya. This scenario is incompatible with a complete biological replacement starting about 700 BP of Paleo-Inuit like the Dorset by the Thule-Inuit (Neo-Inuit), as is currently thought, and more parsimonious with gene flow carrying the NWT01 mutation from a pre-Thule population to the ancestors of the present-day Inuit.

Development and validation of a new multiplex for upgrading Y-STRs population databases from 12 to 23 markers and its forensic casework application

Y chromosomal short tandem repeats (Y-STRs) are used in forensic investigations as a useful complementary tool to autosomal markers. The ongoing development of new kits with an increasing number of markers makes it necessary to update populations typed in the Y-STR Haplotype Reference Database to reach at least 23 Y-STRs. A novel Y-STR multiplex panel was developed to offer a cost-efficient alternative to update Y-STR haplotypes from 12 to 23 loci. This panel includes the eleven markers, DYS448, DYS456, DYS458, DYS635, Y-GATA H4, DYS576, DYS481, DYS549, DYS533, DYS570 and DYS643, as well as DYS385a/b for traceability purpose. Developmental validation of this panel was conducted following the recommendations of the Scientific Working Group on DNA Analysis Methods (SWGDAM), showing high sensitivity, tolerance to common inhibitors as well as high species specificity. It was efficient for degraded DNA samples and for detection of male mixtures. When applying it for extending the current data of the Ibiza population, both the discrimination capacity and the haplotype diversity increased from 0.5952 to 0.9048 and from 0.9808 to 0.9977, respectively. Together, the study demonstrates the suitability of this panel in forensic casework.

Genetic diversity and forensic application of Y-filer STRs in four major ethnic groups of Pakistan

17 Y-chromosomal STRs which are part of the Yfiler Amplification Kit were investigated in 493 unrelated Pakistani individuals belonging to the Punjabi, Sindhi, Baloch, and Pathan ethnic groups. We have assessed the forensic parameters and population genetic structure for each group. Among the 493 unrelated individuals from four ethnic groups (128 Baloch, 122 Pathan, 108 Punjabi, and 135 Sindhi), 82 haplotypes were observed with haplotype diversity (HD) of 0.9906 in Baloch, 102 haplotypes with HD value of 0.9957 in Pathans, 80 haplotypes with HD value of 0.9924 in Punjabi, and 105 haplotypes with HD value of 0.9945 in the Sindhi population. The overall gene diversity for Baloch, Pathan, Punjabi, and Sindhi populations was 0.6367, 0.6479, 0.6657, and 0.6112, respectively. The results had shown us that Pakistani populations do not have a unique set of genes but share the genetic affinity with regional (Central Asia and Northern India) populations. The observed low gene diversity (heterozygosity) values may be because of endogamy trends and this observation is equally supported by the results of forensic parameters which are mostly static across 4 combinations (minimal STRs, extended 11 Y-STRs, Powerplex 12 Y System, and Yfiler 17 Y-STRs) of STRs in these four populations.

Mutation Rate Analysis of RM Y-STRs in Deep-Rooted Multi-Generational Punjabi Pedigrees from Pakistan

Y chromosome short tandem repeat polymorphisms (Y-STRs) are important in many areas of human genetics. Y chromosomal STRs, being normally utilized in the field of forensics, exhibit low haplotype diversity in consanguineous populations and fail to discriminate among male relatives from the same pedigree. Rapidly mutating Y-STRs (RM Y-STRs) have received much attention in the past decade. These 13 RM Y-STRs have high mutation rates (>10−2) and have considerably higher haplotype diversity and discrimination capacity than conventionally used Y-STRs, showing remarkable power when it comes to differentiation in paternal lineages in endogamous populations. Previously, we analyzed two to four generations of 99 pedigrees with 1568 pairs of men covering one to six meioses from all over Pakistan and 216 male relatives from 18 deep-rooted endogamous Sindhi pedigrees covering one to seven meioses. Here, we present 861 pairs of men from 62 endogamous pedigrees covering one to six meioses from the Punjabi population of Punjab, Pakistan. Mutations were frequently observed at DYF399 and DYF403, while no mutation was observed at DYS526a/b. The rate of differentiation ranged from 29.70% (first meiosis) to 80.95% (fifth meiosis), while overall (first to sixth meiosis) differentiation was 59.46%. Combining previously published data with newly generated data, the overall differentiation rate was 38.79% based on 5176 pairs of men related by 1−20 meioses, while Yfiler differentiation was 9.24% based on 3864 pairs. Using father−son pair data from the present and previous studies, we also provide updated RM Y-STR mutation rates.

Contrasting maternal and paternal genetic histories among five ethnic groups from Khyber Pakhtunkhwa, Pakistan

Northwest Pakistan has served as a point of entry to South Asia for different populations since ancient times. However, relatively little is known about the population genetic history of the people residing within this region. To better understand human dispersal in the region within the broader history of the subcontinent, we analyzed mtDNA diversity in 659 and Y-chromosome diversity in 678 individuals, respectively, from five ethnic groups (Gujars, Jadoons, Syeds, Tanolis and Yousafzais), from Swabi and Buner Districts, Khyber Pakhtunkhwa Province, Pakistan. The mtDNAs of all individuals were subject to control region sequencing and SNP genotyping, while Y-chromosomes were analyzed using 54 SNPs and 19 STR loci. The majority of the mtDNAs belonged to West Eurasian haplogroups, with the rest belonging to either South or East Asian lineages. Four of the five Pakistani populations (Gujars, Jadoons, Syeds, Yousafzais) possessed strong maternal genetic affinities with other Pakistani and Central Asian populations, whereas one (Tanolis) did not. Four haplogroups (R1a, R1b, O3, L) among the 11 Y-chromosome lineages observed among these five ethnic groups contributed substantially to their paternal genetic makeup. Gujars, Syeds and Yousafzais showed strong paternal genetic affinities with other Pakistani and Central Asian populations, whereas Jadoons and Tanolis had close affinities with Turkmen populations from Central Asia and ethnic groups from northeast India. We evaluate these genetic data in the context of historical and archeological evidence to test different hypotheses concerning their origins and biological relationships.

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.

Comparative Y-chromosome analysis among Cypriots in the context of historical events and migrations

Y-chromosome analysis provides valuable information regarding the migration patterns of male ancestors, ranging from the Paleolithic age to the modern humans. STR and SNP genotyping analysis provides data regarding the genetic and geographical ancestry of the populations studied. This study focused on the analysis of the Y-chromosome in Maronite Cypriots and Armenian Cypriots, who came to the island as a result of different historical events. The aim was to provide information on the paternal genetic ancestry of Maronites and Armenians of Cyprus and investigate any affinity with the Greek Cypriots and Turkish Cypriots of the island. Since there is limited information in the current literature, we proceeded and used 23 Y-chromosome STRs and 28 Y-chromosome SNPs to genotype 57 Maronite Cypriots and 56 Armenian Cypriots, which were then compared to data from 344 Greek Cypriots and 380 Turkish Cypriots. All samples were assigned to eight major Y-haplogroups but the most frequent haplogroup among all Cypriots is haplogroup J in the major subclade J2a-L559. The calculated pairwise genetic distances between the populations show that Armenian Cypriots are genetically closer to Greek and Turkish Cypriots compared to Maronite Cypriots. Median Joining Network analysis in 17 Y-STR haplotypes of all Cypriots assigned to J2a-L559, revealed that Cypriots share a common paternal ancestor, prior to the migration of the Armenians and Maronites to Cyprus, estimated in the Late Bronze Age and Early Iron Age.

The Y-STR landscape of coastal southeastern Han: Forensic characteristics, haplotype analyses, mutation rates, and population genetics

The Y-STR landscape of Coastal Southeastern Han (CSEH) living in Chinese southeast areas (including Guangdong, Fujian, and Zhejiang provinces) is still unclear. We investigated 62 Y-STR markers in a reasonably large number of 1021 unrelated males and 1027 DNA-confirmed father-son pairs to broaden the genetic backgrounds of CSEH. In total, 85 null alleles, 121 off-ladder alleles, and 95 copy number variants were observed, and 1012 distinct haplotypes were determined with the overall HD and DC values of 0.999974 and 0.9912. We observed 369 mutations in 76 099 meiotic transfers, and the average estimated Y-STR mutation rate was 4.85 × 10^-3 (95% CI, 4.4 × 10^-3 -5.4 × 10^-3 ). The Spearman correlation analyses indicated that GD values (R² = 0.6548) and average allele sizes (R² = 0.5989) have positive correlations with Y-STR mutation rates. Our RM Y-STR set including 8 candidate RM Y-STRs, of which DYS534, DYS630, and DYS713 are new candidates in CSEH, distinguished 18.52% of father-son pairs. This study also clarified the population structures of CSEH which isolated in population-mixed South China relatively. The strategy, SM Y-STRs for familial searching and RM Y-STRs for individual identification regionally, could be applicable based on enough knowledge of the Y-STR mutability of different populations.

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.

Mutational data and population profiling of 23 Y-STRs in three Brazilian populations

Y-chromosomal STRs are important markers in forensic genetics, due to some peculiar characteristics. The absence of recombination makes them a useful tool to infer kinship in complex cases involving distant paternal relatives, or to infer paternal bio-geographic ancestry. The presence of a single copy, being transmitted from father to son, allow tracing mutational events in Y-STRs without ambiguity. For the statistical interpretation of forensic evidences based on Y-STR profiles, it is necessary to have estimates on both mutation rates and haplotype frequencies. In this work, 407 father-son duos from São Paulo and Rio de Janeiro states and 204 unrelated individuals from Manaus were analyzed. Haplotype frequencies and mutation rates for the Y-STRs from the PowerPlex Y23 commercial kit were estimated. Thirty-six mutations were observed in 15 of the 22 Y-STRs analyzed, for an average mutation rate of 3.84 × 10^-3 (95 % CI 2.69 × 10^-3 to 5.32 × 10^-3). All mutations in GAAA repeats occurred in alleles with 13 or more uninterrupted units. Mutations in GATA repeats were observed in alleles with 9-17 uninterrupted units. An analysis carried out in different father's age groups showed an increase of 2.48 times the mutation rate in the age group of 40-50 years, when compared to the 20-30 age group, in agreement with the described for autosomal STRs. A high haplotype diversity was found in the three Brazilian populations. Pairwise genetic distance analysis (F_ST) showed no significant differences between the three populations in this study, which were also close to populations with strong European influence. The highest distances among the Brazilian populations were with São Gabriel da Cachoeira, which has a high Native American ancestry.

Recovery & identification of human Y-STR DNA from immatures of chrysomya albiceps (Diptera: Calliphoridae). Simulation of sexual crime investigation involving victim corpse in state of decay

The number of sexual crimes in Brazil, as in several other countries, is very high. In many of these crimes the women raped are murdered and their bodies are found days later, in an advanced state of decomposition, with intense cadaverous fauna. Forensic Entomology studies insects and other arthropods that can be used in the expert analysis of various types of crimes. Diptera, the order of insects that comprises the two-winged or true flies, represents one of the largest known groups of insects and is the principal source of cadaveric entomofauna. Members of its Calliphoridae family are observed in cadavers in all phases of decomposition. The retrieval and identification of human Y-STR DNA from the gastrointestinal tract of Calliphoridae species Chrysomya albiceps maggots and pupae can provide a good tool for the gathering of evidence in sexual crime investigations involving rape and death, in which the abandoned victim's body is found in a putrefied state. In this study, the animal model used was a female pig, Sus scrofa, which was sacrificed in a forested area with three shots from a 0.40 calibre Taurus pistol, and inoculated with semen to its anal and vaginal regions, simulating rape and homicide. During decomposition, 20-80 maggots were collected every 24 h and preserved in 70 % alcohol, totalling 289 maggots and 157 pupae (446 immatures) over a period of 14 days (336 h) of decomposition. Each maggot was then dissected for removal of the digestive tract, which was placed in extraction buffer. The molecular phase proceeded with extraction, quantification, amplification and capillary electrophoresis of samples, testing 16 STR loci of the Y chromosome. It was possible to establish a partial Y-STR DNA profile, with the amplification of up to eight sites, by considering a combination of the samples taken at hours 144 h, 168 h, 192 h, 216 h, 240 h, 288 h, 312 h and 336 h.

Analysis of Mutation Rate of 17 Y-Chromosome Short Tandem Repeats Loci Using Tanzanian Father-Son Paired Samples

Hundred unrelated father-son buccal swab sample pairs collected from consented Tanzanian population were examined to establish mutation rates using 17 Y-STRs loci DYS19, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS385a, DYS385b, DYS437, DYS438, DYS439, DYS448, DYS456, DYS458, DYS635, and Y-GATA-H4 of the AmpFlSTRYfiler kit used in forensics and paternity testing. Prior to 17 Y-STRs analysis, father-son pair biological relationships were confirmed using 15 autosomal STRs markers and found to be paternally related. A total of four single repeat mutational events were observed between father and sons. Two mutations resulted in the gain of a repeat and the other two resulted in a loss of a repeat in the son. All observed mutations occurred at tetranucleotide loci DYS389II, DYS385a, and DYS385b. The locus specific mutation rate varied between 0 and 1.176 x10-3 and the average mutation rate of 17Y-STRs loci in the present study was 2.353x10-3 (6.41x10-4 - 6.013x10-3) at 95% CI. Furthermore the mean fathers' age with at least one mutation at son's birth was 32 years with standard error of 2.387 while the average age of all fathers without mutation in a sampled population at son's birth was 26.781 years with standard error of 0.609. The results shows that fathers' age at son's birth may have an effect on Y-STRs mutation rate analysis, though this age difference was statistically not significant using unpaired samples t-test (p = 0.05). As a consequence of observed mutation rates in this study, the precise and reliable understanding of mutation rate at Y-chromosome STR loci is necessary for a correct evaluation and interpretation of DNA typing results in forensics and paternity testing involving males. The criterion for exclusion in paternity testing should be defined, so that an exclusion from paternity has to be based on exclusion constellations at a minimum of two 17 Y-STRs loci.

Discriminating power of rapidly mutating Y-STRs in deep rooted endogamous pedigrees from Sindhi population of Pakistan

Rapidly mutating Y-STRs (RM Y-STRs) have been paid much attention in recent years. The 13 RM Y-STRs (DYF387S1, DYF399S1, DYF403S1a/b, DYF404S1, DYS449, DYS518, DYS526I/II, DYS547, DYS570, DYS576, DYS612, DYS626, and DYS627) have been proved to have substantially higher haplotype diversity and discrimination capacity than conventionally used Y-STRs indicating the considerable power in paternal lineage differentiation in endogamous populations, separation of which is usually impossible with standard Y-STRs. In current study, we analyzed the RM Y-STRs and PowerPlex® Y23 System in 216 male relatives from 18 deep rooted endogamous Sindhi families from Pakistan. Mutations were frequently observed at DYF399S1, DYS449, DYS518DYS547 and DYF403S1b2 loci, which are known to mutate more rapidly than other RM Y-STRs. Overall differentiation rate with RM Y-STRs was as high as 32.88%, while those with PowerPlex® Y23 System and AmpFℓSTR® Yfiler™ kit were 6.85% and 3.65% respectively. The differentiation rate of RM Y-STRs was 29.22% and 26.03% higher than those of AmpFlSTR® Yfiler™ kit and PowerPlex® Y23 System, respectively.

Haplotypic polymorphisms and mutation rate estimates of 22 Y-chromosome STRs in the Northern Chinese Han father-son pairs

Y chromosome Short tandem repeats (Y-STRs) analysis has been widely used in forensic identification, kinship testing, and population evolution. An accurate understanding of haplotype and mutation rate will benefit these applications. In this work, we analyzed 1123 male samples from Northern Chinese Han population which including 578 DNA-confirmed father-son pairs at 22 Y-STRs loci. A total of 537 haplotypes were observed and the overall haplotype diversity was calculated as 1.0000 ± 0.0001. Except that only two haplotypes were observed twice, all the rest of the 535 were unique. Furthermore, totally 47 mutations were observed during 13,872 paternal meiosis. The mutation rate for each locus estimates ranged from 0.0 to 15.6 × 10-3 with an average mutation rate 3.4 × 10-3 (95% CI 2.5-4.5 × 10-3). Among the 22 loci, DYS449, DYS389 II and DYS458 are the most prone to mutations. This study adds to the growing data on Y-STR haplotype diversity and mutation rates and could be very useful for population and forensic genetics.

Inconsistent genotyping call at DYS389 locus and implications for interpretation

The male-specific Y chromosome short tandem repeat (STR) locus is used widely in forensic case, which are useful molecular tool to providing the biological evidence for male/female mixture and paternal lineage cases. The Y-STR analysis has been greatly facilitated by advent of commercial multiplex kit. However, even with well-designed robust multiplex kit, abnormal genotyping profile may be observed when encountering with mutations, such as deletion/duplication within the target region or mutation at the primer binding site. In this study, a single-allele shift by five nucleotides for the DYS389I marker between the AmpFlSTR^® Yfiler^® and Yfiler^® Plus PCR amplification kits while the same allele count for DYS389II was observed in eight unrelated Chinese male individuals. After further investigations by re-amplified with three additional multiplex kits, sanger, and next-generation sequencing, the discordance was finally proven caused by existing rare mutation in those sample, which contained two adjacent SNPs only one base apart in the sequence. This paper describes the molecular basis of the discordance at DYS389I genotyping between different commercial multiplex kits and could provide available information for enhancing of interpretation of abnormal Y-STR genotyping in forensic practice.

Y-chromosomal analysis of Greek Cypriots reveals a primarily common pre-Ottoman paternal ancestry with Turkish Cypriots

Genetics can provide invaluable information on the ancestry of the current inhabitants of Cyprus. A Y-chromosome analysis was performed to (i) determine paternal ancestry among the Greek Cypriot (GCy) community in the context of the Central and Eastern Mediterranean and the Near East; and (ii) identify genetic similarities and differences between Greek Cypriots (GCy) and Turkish Cypriots (TCy). Our haplotype-based analysis has revealed that GCy and TCy patrilineages derive primarily from a single gene pool and show very close genetic affinity (low genetic differentiation) to Calabrian Italian and Lebanese patrilineages. In terms of more recent (past millennium) ancestry, as indicated by Y-haplotype sharing, GCy and TCy share much more haplotypes between them than with any surrounding population (7-8% of total haplotypes shared), while TCy also share around 3% of haplotypes with mainland Turks, and to a lesser extent with North Africans. In terms of Y-haplogroup frequencies, again GCy and TCy show very similar distributions, with the predominant haplogroups in both being J2a-M410, E-M78, and G2-P287. Overall, GCy also have a similar Y-haplogroup distribution to non-Turkic Anatolian and Southwest Caucasian populations, as well as Cretan Greeks. TCy show a slight shift towards Turkish populations, due to the presence of Eastern Eurasian (some of which of possible Ottoman origin) Y-haplogroups. Overall, the Y-chromosome analysis performed, using both Y-STR haplotype and binary Y-haplogroup data puts Cypriot in the middle of a genetic continuum stretching from the Levant to Southeast Europe and reveals that despite some differences in haplotype sharing and haplogroup structure, Greek Cypriots and Turkish Cypriots share primarily a common pre-Ottoman paternal ancestry.

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.

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

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

Forensic use of Y-chromosome DNA: a general overview

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

Mayans: a Y chromosome perspective

In spite of the wealth of available cultural and archeological information as well as general interest in the Mayans, little is known about their genetics. In this study, for the first time, we attempt to alleviate this lacuna of knowledge by comprehensively investigating the Y chromosome composition of contemporary Mayan populations throughout their domain. To accomplish this, five geographically targeted and ethnically distinct Mayan populations are investigated using Y-SNP and Y-STR markers.

FINDINGS

overall, the Mayan populations as a group are highly homogeneous, basically made up of only two autochthonous haplogroups, Q1a2a1a1*-M3 and Q1a2a1*-L54. Although the Y-STR data illustrates diversity, this diversity, for the most part, is uniformly distributed among geographically distant Mayan populations. Similar haplotypes among populations, abundance of singletons and absence of population partitioning within networks among Mayan populations suggest recent population expansion and substantial gene flow within the Mayan dominion, possibly due to the development of agriculture, the establishment of interacting City-State systems and commerce.

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.

Analysis of Rapidly Mutating Y Chromosome Short Tandem Repeats (RM Y-STRs)

Rapidly mutating (RM) Y-STRs were recently identified and can help to differentiate between paternally related males in paternity and forensic casework. Normal Y-STRs are often used in casework due to their paternal inheritance, which can help to resolve kinship cases, and identify male components in male/female mixtures; however, more discriminating profiles are obtained if rapidly mutating Y-STRs are used. Previously two or three multiplex PCRs have been used to amplify 13 RM Y-STRs; here, an assay amplifying these 13 markers in a single multiplex PCR is described. Commercially available male control DNA samples have been genotyped during the validation of this assay, thus providing a tool for calibrating genotyping results. It is expected that the assay will provide a niche tool for genotyping casework samples.

Genetic structure of the Kuwaiti population revealed by paternal lineages

OBJECTIVE

We analyzed the Y-chromosome haplogroup diversity in the Kuwaiti population to gain a more complete overview of its genetic landscape.

METHOD

A sample of 117 males from the Kuwaiti population was studied through the analysis of 22 Y-SNPs. The results were then interpreted in conjunction with those of other populations from the Middle East, South Asia, North and East Africa, and East Europe.

RESULTS

The analyzed markers allowed the discrimination of 19 different haplogroups with a diversity of 0.7713. J-M304 was the most frequent haplogroup in the Kuwaiti population (55.5%) followed by E-M96 (18%). They revealed a genetic homogeneity between the Kuwaiti population and those of the Middle East (FST  = 6.1%, P-value < 0.0001), although a significant correlation between genetic and geographic distances was found (r = 0.41, P-value = 0.009). Moreover, the nonsignificant pairwise FST genetic distances between the Kuwait population on the one hand and the Arabs of Iran and those of Sudan on the other, corroborate the hypothesis of bidirectional gene flow between Arabia and both Iran and Sudan.

CONCLUSION

Overall, we have revealed that the Kuwaiti population has experienced significant gene flow from neighboring populations like Saudi Arabia, Iran, and East Africa. Therefore, we have confirmed that the population of Kuwait is genetically coextensive with those of the Middle East.

New clues to the evolutionary history of the main European paternal lineage M269: dissection of the Y-SNP S116 in Atlantic Europe and Iberia

The dissection of S116 in more than 1500 individuals from Atlantic Europe and the Iberian Peninsula has provided important clues about the controversial evolutionary history of M269. First, the results do not point to an origin of M269 in the Franco-Cantabrian refuge, owing to the lack of sublineage diversity within M269, which supports the new theories proposing its origin in Eastern Europe. Second, S116 shows frequency peaks and spatial distribution that differ from those previously proposed, indicating an origin farther west, and it also shows a high frequency in the Atlantic coastline. Third, an outstanding frequency of the DF27 sublineage has been found in Iberia, with a restricted distribution pattern inside this peninsula and a frequency maximum in the area of the Franco-Cantabrian refuge. This entire panorama indicates an old arrival of M269 into Western Europe, because it has generated at least two episodes of expansion in the Franco-Cantabrian area. This study demonstrates the importance of continuing the dissection of the M269 lineage in different European populations because the discovery and study of new sublineages can adjust or even completely revise the theories about European peopling, as has been the case for the place of origin of M269.

Ladakh, India: the land of high passes and genetic heterogeneity reveals a confluence of migrations

Owing to its geographic location near the longitudinal center of Asia, Ladakh, the land of high passes, has witnessed numerous demographic movements during the past millenniums of occupation. In an effort to view Ladakh's multicultural history from a paternal genetic perspective, we performed a high-resolution Y-chromosomal survey of Ladakh, within the context of Y haplogroup and haplotype distributions of 41 Asian reference populations. The results of this investigation highlight the rich ethnic and genetic diversity of Ladkah which includes genetic contributions from disparate regions of the continent including, West, East, South and Central Asia. The phylogenetic signals from Ladakh are consistent with the Indo-Aryans' occupation during the Neolithic age and its historic connection with Tibet, as well as the East-West gene flow associated with the Silk Road.

Three hundred years of low non-paternity in a human population

When cuckoldry is frequent we can expect fathers to withhold investment in offspring that may not be theirs. Human paternal investment can be substantial and is in line with observations from tens of thousands of conceptions that suggest that cuckoldry is rare in humans. The generality of this claim seems to be in question as the rate of cuckoldry varies across populations and studies have mostly been on Western populations. Two additional factors complicate our conclusions, (1) current estimates of the rate of cuckoldry in humans may not reflect our past behaviour as adultery can be concealed by the use of contraceptives; and (2) it is difficult to obtain samples that are random with respect to their paternity certainty. Studies that combine genealogies with Y-chromosome haplotyping are able to circumvent some of these problems by probing into humans' historical behaviour. Here we use this approach to investigate 1273 conceptions over a period of 330 years in 23 families of the Afrikaner population in South Africa. We use haplotype frequency and diversity and coalescent simulations to show that the male population did not undergo a severe bottleneck and that paternity exclusion rates are high for this population. The rate of cuckoldry in this Western population was 0.9% (95% confidence interval 0.4-1.5%), and we argue that given the current data on historical populations we have to conclude that, at least for Western human populations, cuckoldry rate is probably in the range of 1%.

From Arabia to Iberia: A Y chromosome perspective

At different times during recent human evolution, northern Africa has served as a conduit for migrations from the Arabian Peninsula. Although previous researchers have investigated the possibility of the Strait of Gibraltar as a pathway of migration from North Africa to Iberia, we now revisit this issue and theorize that although the Strait of Gibraltar, at the west end of this corridor, has acted as a barrier for human dispersal into Southwest Europe, it has not provided an absolute seal to gene flow. To test this hypothesis, here we use the spatial frequency distributions, STR diversity and expansion time estimates of Y chromosome haplogroups J1-P58 and E-M81 to investigate the genetic imprints left by the Arabian and Berber expansions into the Iberian Peninsula, respectively. The data generated indicate that Arabian and Berber genetic markers are detected in Iberia. We present evidence that suggest that Iberia has received gene flow from Northwest Africa during and prior to the Islamic colonization of 711A.D. It is interesting that the highest frequencies of Arabia and Berber markers are not found in southern Spain, where Islam remained the longest and was culturally most influential, but in Northwest Iberia, specifically Galicia. We propose that Moriscos' relocations to the north during the Reconquista, the migration of cryptic Muslims seeking refuge in a more lenient society and/or more geographic extensive pre-Islamic incursions may explain the higher frequencies and older time estimates of mutations in the north of the Peninsula. These scenarios are congruent with the higher diversities of some diagnostic makers observed in Northwest Iberia.

A novel multiplex assay amplifying 13 Y-STRs characterized by rapid and moderate mutation rate

As microsatellites located on Y chromosome mutate with different rates, they may be exploited in evolutionary studies, genealogical testing of a variety of populations and even, as proven recently, aid individual identification. Currently available commercial Y-STR kits encompass mostly low to moderately mutating loci, making them a perfect choice for the first two applications. Some attempts have been made so far to utilize Y-STRs to provide a discriminatory tool for forensic purposes. Although all 13 rapidly mutating Y-STRs were already multiplexed, no single assay based on single-copy markers allowing at least a portion of close male relatives to be differentiated from one another is available. To fill in the blanks, we constructed and validated an assay comprised of single-copy Y-STR markers only with a mutation rate ranging from 8×10(-3) to 1×10(-2). Performance of the resulting combination of nine RM Y-STRs and four moderately mutating ones was tested on 361 father-son pairs and 1326 males from 9 populations revealing an overall mutation rate of 1.607×10(-1) for the assay as a whole. Application of the proposed 13 Y-STR set to differentiation of haplotypes present among homogenous population of Buryats resulted in a threefold increase of discrimination as compared with 10 Y-STRs from the PowerPlex(®) Y.