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.

Around 95% of the male-specific Y-chromosome (chrY) is non-recombining and therefore inherited in a conserved manner from father to son. It can therefore serve as a powerful marker for interdisciplinary genetic-genealogical research as it provides a strong link between genetic information and a family tree or pedigree. While Y-chromosomal short tandem repeats (Y-STRs) discriminate close paternal kinships, single nucleotide polymorphisms (Y-SNPs) enables the identification of far evolutionary ancestry. Unfortunately, an extensive chrY-specific sequencing panel combining a large number of familial Y-STRs and evolutionary Y-SNPs was not yet available. Therefore, chrY is rarely included in research projects and not often linked to a genealogical, history-demographical or life science database. In this way, the importance of chrY still remains not yet fully understood. Massive parallel sequencing (MPS) allows the simultaneous analysis at sequence level of Y-SNPs and Y-STRs with variable mutation rates in a large number of males. However, up until today, no commercial kit is exploiting the full potential that MPS offers on chrY. Therefore, we developed the ‘CSYseq’, which is the first extensive chrY-specific sequencing panel. The CSYseq simultaneously identifies 9,014 slow mutating Y-SNPs to identify evolutionary ancestry, and 202 rapid mutating Y-STRs to investigate paternal relationships. We validated and optimized the panel through the analysis of 130 males distributed over 65 families. This novel MPS panel is useful for biogeographical identity and ancestry analysis, together with Y-chromosome profiling for the identification of patrilineages and discrimination of closely related males. As the CSYseq includes a very diverse set of markers that can be easily interpreted, it is interesting for different interdisciplinary applications within evolutionary, population, molecular, medical and forensic genetics.

Cell survival and DNA damage repair are promoted in the human blood thanatotranscriptome shortly after death

RNA analysis of post-mortem tissues, or thanatotranscriptomics, has become a topic of interest in forensic science due to the essential information it can provide in forensic investigations. Several studies have previously investigated the effect of death on gene transcription, but it has never been conducted with samples of the same individual. For the first time, a longitudinal mRNA expression analysis study was performed with post-mortem human blood samples from individuals with a known time of death. The results reveal that, after death, two clearly differentiated groups of up- and down-regulated genes can be detected. Pathway analysis suggests active processes that promote cell survival and DNA damage repair, rather than passive degradation, are the source of early post-mortem changes of gene expression in blood. In addition, a generalized linear model with an elastic net restriction predicted post-mortem interval with a root mean square error of 4.75 h. In conclusion, we demonstrate that post-mortem gene expression data can be used as biomarkers to estimate the post-mortem interval though further validation using independent sample sets is required before use in forensic casework.

At the Interface of Life and Death: Post-mortem and Other Applications of Vaginal, Skin, and Salivary Microbiome Analysis in Forensics

Microbial forensics represents a promising tool to strengthen traditional forensic investigative methods and fill related knowledge gaps. Large-scale microbiome studies indicate that microbial fingerprinting can assist forensics in areas such as trace evidence, source tracking, geolocation, and circumstances of death. Nevertheless, the majority of forensic microbiome studies focus on soil and internal organ samples, whereas the microbiome of skin, mouth, and especially vaginal samples that are routinely collected in sexual assault and femicide cases remain underexplored. This review discusses the current and emerging insights into vaginal, skin, and salivary microbiome-modulating factors during life (e.g., lifestyle and health status) and after death (e.g., environmental influences and post-mortem interval) based on next-generation sequencing. We specifically highlight the key aspects of female reproductive tract, skin, and mouth microbiome samples relevant in forensics. To fill the current knowledge gaps, future research should focus on the degree to which the post-mortem succession rate and profiles of vaginal, skin, and saliva microbiota are sensitive to abiotic and biotic factors, presence or absence of oxygen and other gases, and the nutrient richness of the environment. Application of this microbiome-related knowledge could provide valuable complementary data to strengthen forensic cases, for example, to shed light on the circumstances surrounding death with (post-mortem) microbial fingerprinting. Overall, this review synthesizes the present knowledge and aims to provide a framework to adequately comprehend the hurdles and potential application of vaginal, skin, and salivary post-mortem microbiomes in forensic investigations.

YMrCA: Improving Y-chromosomal ancestor time estimation for DNA kinship research

The Y-chromosome is a valuable kinship indicator in family history and forensic research. To reconstruct genealogies, the time to the most recent common ancestor (tMRCA) between paternal relatives can be estimated through Y-STR analysis. Existing models are the stepwise mutation model (SMM, only one-step Y-STR changes) and the infinite allele model (IAM, new allele per Y-STR change). In this study, these mutation models and all existing tMRCA calculators were validated through a genetic-genealogy database containing 1,120 biologically related genealogical pairs confirmed by 46 Y-STRs with known tMRCA (18,109 generations). Consistent under- and overestimation and broad confidence intervals were observed, leading to dubious tMRCA estimates. This is because they do not include individual mutation rates or multi-step changes and ignore hidden multiple, back, or parallel modifications. To improve tMRCA estimation, we developed a user-friendly calculator, the "YMrCA", including all previously mentioned mutation characteristics. After extensive validation, we observed that the YMrCA calculator demonstrated a promising performance. The YMrCA yields a significantly higher tMRCA success rate (96%; +20%) and a lower tMRCA error (7; -3) compared to the mutation models and all online tMRCA calculators. Therefore, YMrCA offers the next step towards more objective tMRCA estimation for DNA kinship research.

Ysurnames? The patrilineal Y-chromosome and surname correlation for DNA kinship research

The Y-chromosome is a widely studied and useful small part of the genome providing different applications for interdisciplinary research. In many (Western) societies, the Y-chromosome and surnames are paternally co-inherited, suggesting a corresponding Y-haplotype for every namesake. While it has already been observed that this correlation may be disrupted by a false-paternity event, adoption, anonymous sperm donor or the co-founding of surnames, extensive information on the strength of the surname match frequency (SMF) with the Y-chromosome remains rather unknown. For the first time in Belgium and the Netherlands, we were able to study this correlation using 2,401 males genotyped for 46 Y-STRs and 183 Y-SNPs. The SMF was observed to be dependent on the number of Y-STRs analyzed, their mutation rates and the number of Y-STR differences allowed for a kinship. For a perfect match, the Yfiler® Plus and our in-house YForGen kit gave a similar high SMF of 98%, but for non-perfect matches, the latter could overall be identified as the best kit. The SMF generally increased due to less mismatches when encountering [1] deep Y-subhaplogroups, [2] less frequently occurring surnames, and [3] small geographical distances between relatives. This novel information enabled the design of a surname prediction model based on genetic and geographical distances of a kinship. The prediction model has an area under the curve (AUC) of 0.9 and is therefore useable for DNA kinship priority listing in estimation applications like forensic familial searching.

An interdisciplinary study around the reliquary of the late cardinal Jacques de Vitry

The reliquary of Jacques de Vitry, a prominent clergyman and theologian in the early 13th century, has experienced several transfers over the last centuries, which seriously question the attribution of the remains to the late Cardinal. Uncertainty about the year of his birth poses an additional question regarding his age at death in 1240. The reliquary, located in the Saint Marie d'Oigines church, Belgium, was reopened in 2015 for an interdisciplinary study around his relics as well as the Treasure of Oignies, a remarkable cultural heritage notably built from Jacques de Vitry's donation. Anthropological, isotopic and genetic analyses were performed independently on the remains found in the reliquary. Results of the analyses provided evidence that the likelihood that these remains are those of Jacques de Vitry is very high: the remains belong to the same human male individual and the historical tradition about his age is confirmed. In addition, a separate relic (left tibia) was analysed and found to match with the remains of the reliquary (right tibia). The unique Jacques de Vitry's mitre, made of parchment, was sampled non-destructively and the extracted parchment collagen was analysed by a proteomic method in order to determine the animal species. The results showed that, surprisingly, not all parts of the mitre were made from the same species. All together, these findings are expected to fertilize knowledge carried by historical tradition around the relics of Jacques de Vitry and his related cultural heritage.

The black legend on the Spanish presence in the low countries: Verifying shared beliefs on genetic ancestry

OBJECTIVES

War atrocities committed by the Spanish army in the Low Countries during the 16th century are so ingrained in the collective memory of Belgian and Dutch societies that they generally assume a signature of this history to be present in their genetic ancestry. Historians claim this assumption is a consequence of the so-called "Black Legend" and negative propaganda portraying and remembering Spanish soldiers as extreme sexual aggressors. The impact of the presence of Spaniards during the Dutch Revolt on the genetic variation in the Low Countries has been verified in this study.

MATERIALS AND METHODS

A recent population genetic analysis of Iberian-associated Y-chromosomal variation among Europe is enlarged with representative samples of Dutch (N = 250) and Flemish (N = 1,087) males. Frequencies of these variants are also compared between donors whose oldest reported paternal ancestors lived in-nowadays Flemish-cities affected by so-called Spanish Furies (N = 116) versus other patrilineages in current Flemish territory (N = 971).

RESULTS

The frequencies of Y-chromosomal markers Z195 and SRY2627 decline steeply going north from Spain and the data for the Flemish and Dutch populations fits within this pattern. No trend of higher frequencies of these variants has been found within the well-ascertained samples associated with Spanish Fury cities.

DISCUSSION

Although sexual aggression did occur in the 16th century, these activities did not leave a traceable "Spanish" genetic signature in the autochthonous genome of the Low Countries. Our results support the view that the 'Black Legend' and historical propaganda on sexual aggression have nurtured today's incorrect assumptions regarding genetic ancestry.

Genetic-genealogy approach reveals low rate of extrapair paternity in historical Dutch populations

OBJECTIVE

Evolutionary theory has shown that seeking out extrapair paternity (EPP) can be a viable reproductive strategy for both sexes in pair-bonded species, also in humans. As yet, estimates of the contemporary or historical EPP rate in human population are still rare. In the present study, we estimated the historical EPP rate in the Dutch population over the last 400 years and compared the rate with those obtained for other human populations to determine the evolutionary, cultural, and socio-demographic factors that influence human cuckoldry behavior.

METHODS

We estimated the historical EPP rate for the Dutch population via the "genealogical pair method", in which the EPP rate is derived from Y-chromosome mismatches between pairs of individuals that, based on genealogical evidence, share a common paternal ancestor.

RESULTS

Based on the analysis of 68 representative genealogical pairs, separated by a total of 1013 fertilization events, we estimated that the historical EPP rate for the Dutch population over the last 400 years was 0.96% per generation (95% confidence interval 0.46%-1.76%).

CONCLUSION

The Dutch EPP rate fits perfectly within the range reported for other contemporary and historical populations in Western Europe and was highly congruent with that estimated for neighboring Flanders, despite the socio-economic and religious differences between both populations. The estimated low EPP rate challenges the "dual mating strategy hypothesis" that states that women could obtain fitness benefits by securing investment from one man while cuckolding him to obtain good genes from an affair partner.

Bight of Benin: a Maternal Perspective of Four Beninese Populations and their Genetic Implications on the American Populations of African Ancestry

The understanding of the first movements of the ancestral populations within the African continent is still unclear, particularly in West Africa, due to several factors that have shaped the African genetic pool across time. To improve the genetic representativeness of the Beninese population and to better understand the patterns of human settlement inside West Africa and the dynamics of peopling of the Democratic Republic of Benin, we analyzed the maternal genetic variation of 193 Beninese individuals belonging to Bariba, Berba, Dendi, and Fon populations. Results support the oral traditions indicating that the western neighbouring populations have been the ancestors of the first Beninese populations, and the extant genetic structure of the Beninese populations is most likely the result of admixture between populations from neighbouring countries and native people. The present findings highlight how the Beninese populations contributed to the gene pool of the extant populations of some American populations of African ancestry. This strengthens the hypothesis that the Bight of Benin was not only an assembly point for the slave trade during the Trans-Atlantic Slave Trade but also an important slave trapping area.

Comparing maternal genetic variation across two millennia reveals the demographic history of an ancient human population in southwest Turkey

More than two decades of archaeological research at the site of Sagalassos, in southwest Turkey, resulted in the study of the former urban settlement in all its features. Originally settled in late Classical/early Hellenistic times, possibly from the later fifth century BCE onwards, the city of Sagalassos and its surrounding territory saw empires come and go. The Plague of Justinian in the sixth century CE, which is considered to have caused the death of up to a third of the population in Anatolia, and an earthquake in the seventh century CE, which is attested to have devastated many monuments in the city, may have severely affected the contemporary Sagalassos community. Human occupation continued, however, and Byzantine Sagalassos was eventually abandoned around 1200 CE. In order to investigate whether these historical events resulted in demographic changes across time, we compared the mitochondrial DNA variation of two population samples from Sagalassos (Roman and Middle Byzantine) and a modern sample from the nearby town of Ağlasun. Our analyses revealed no genetic discontinuity across two millennia in the region and Bayesian coalescence-based simulations indicated that a major population decline in the area coincided with the final abandonment of Sagalassos, rather than with the Plague of Justinian or the mentioned earthquake.

Improved age determination of blood and teeth samples using a selected set of DNA methylation markers

Age estimation from DNA methylation markers has seen an exponential growth of interest, not in the least from forensic scientists. The current published assays, however, can still be improved by lowering the number of markers in the assay and by providing more accurate models to predict chronological age. From the published literature we selected 4 age-associated genes (ASPA, PDE4C, ELOVL2, and EDARADD) and determined CpG methylation levels from 206 blood samples of both deceased and living individuals (age range: 0-91 years). This data was subsequently used to compare prediction accuracy with both linear and non-linear regression models. A quadratic regression model in which the methylation levels of ELOVL2 were squared showed the highest accuracy with a Mean Absolute Deviation (MAD) between chronological age and predicted age of 3.75 years and an adjusted R(2) of 0.95. No difference in accuracy was observed for samples obtained either from living and deceased individuals or between the 2 genders. In addition, 29 teeth from different individuals (age range: 19-70 years) were analyzed using the same set of markers resulting in a MAD of 4.86 years and an adjusted R(2) of 0.74. Cross validation of the results obtained from blood samples demonstrated the robustness and reproducibility of the assay. In conclusion, the set of 4 CpG DNA methylation markers is capable of producing highly accurate age predictions for blood samples from deceased and living individuals.

High Y-chromosomal diversity and low relatedness between paternal lineages on a communal scale in the Western European Low Countries during the surname establishment

There is limited knowledge on the biological relatedness between citizens and on the demographical dynamics within villages, towns and cities in pre-17th century Western Europe. By combining Y-chromosomal genotypes, in-depth genealogies and surname data in a strict genetic genealogical approach, it is possible to provide insights into the genetic diversity and the relatedness between indigenous paternal lineages within a particular community at the time of the surname adoption. To obtain these insights, six Flemish communities were selected in this study based on the differences in geography and historical development. After rigorous selection of appropriate DNA donors, low relatedness between Y chromosomes of different surnames was found within each community, although there is co-occurrence of these surnames in each community since the start of the surname adoption between the 14th and 15th century. Next, the high communal diversity in Y-chromosomal lineages was comparable with the regional diversity across Flanders at that time. Moreover, clinal distributions of particular Y-chromosomal lineages between the communities were observed according to the clinal distributions earlier observed across the Flemish regions and Western Europe. No significant indication for genetic differences between communities with distinct historical development was found in the analysis. These genetic results provide relevant information for studies in historical sciences, archaeology, forensic genetics and genealogy.

Towards a consensus Y-chromosomal phylogeny and Y-SNP set in forensics in the next-generation sequencing era

Currently, several different Y-chromosomal phylogenies and haplogroup nomenclatures are presented in scientific literature and at conferences demonstrating the present diversity in Y-chromosomal phylogenetic trees and Y-SNP sets used within forensic and anthropological research. This situation can be ascribed to the exponential growth of the number of Y-SNPs discovered due to mostly next-generation sequencing (NGS) studies. As Y-SNPs and their respective phylogenetic positions are important in forensics, such as for male lineage characterization and paternal bio-geographic ancestry inference, there is a need for forensic geneticists to know how to deal with these newly identified Y-SNPs and phylogenies, especially since these phylogenies are often created with other aims than to carry out forensic genetic research. Therefore, we give here an overview of four categories of currently used Y-chromosomal phylogenies and the associated Y-SNP sets in scientific research in the current NGS era. We compare these categories based on the construction method, their advantages and disadvantages, the disciplines wherein the phylogenetic tree can be used, and their specific relevance for forensic geneticists. Based on this overview, it is clear that an up-to-date reduced tree with a consensus Y-SNP set and a stable nomenclature will be the most appropriate reference resource for forensic research. Initiatives to reach such an international consensus are therefore highly recommended.

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.

Toward male individualization with rapidly mutating y-chromosomal short tandem repeats

Relevant for various areas of human genetics, Y-chromosomal short tandem repeats (Y-STRs) are commonly used for testing close paternal relationships among individuals and populations, and for male lineage identification. However, even the widely used 17-loci Yfiler set cannot resolve individuals and populations completely. Here, 52 centers generated quality-controlled data of 13 rapidly mutating (RM) Y-STRs in 14,644 related and unrelated males from 111 worldwide populations. Strikingly, >99% of the 12,272 unrelated males were completely individualized. Haplotype diversity was extremely high (global: 0.9999985, regional: 0.99836–0.9999988). Haplotype sharing between populations was almost absent except for six (0.05%) of the 12,156 haplotypes. Haplotype sharing within populations was generally rare (0.8% nonunique haplotypes), significantly lower in urban (0.9%) than rural (2.1%) and highest in endogamous groups (14.3%). Analysis of molecular variance revealed 99.98% of variation within populations, 0.018% among populations within groups, and 0.002% among groups. Of the 2,372 newly and 156 previously typed male relative pairs, 29% were differentiated including 27% of the 2,378 father–son pairs. Relative to Yfiler, haplotype diversity was increased in 86% of the populations tested and overall male relative differentiation was raised by 23.5%. Our study demonstrates the value of RM Y-STRs in identifying and separating unrelated and related males and provides a reference database.

Low historical rates of cuckoldry in a Western European human population traced by Y-chromosome and genealogical data

Recent evidence suggests that seeking out extra-pair paternity (EPP) can be a viable alternative reproductive strategy for both males and females in many pair-bonded species, including humans. Accurate data on EPP rates in humans, however, are scant and mostly restricted to extant populations. Here, we provide the first large-scale, unbiased genetic study of historical EPP rates in a Western European human population based on combining Y-chromosomal data to infer genetic patrilineages with genealogical and surname data, which reflect known historical presumed paternity. Using two independent methods, we estimate that over the last few centuries, EPP rates in Flanders (Belgium) were only around 1–2% per generation. This figure is substantially lower than the 8–30% per generation reported in some behavioural studies on historical EPP rates, but comparable with the rates reported by other genetic studies of contemporary Western European populations. These results suggest that human EPP rates have not changed substantially during the last 400 years in Flanders and imply that legal genealogies rarely differ from the biological ones. This result has significant implications for a diverse set of fields, including human population genetics, historical demography, forensic science and human sociobiology.

Controversial identification in a historical case is illustrative of the complexity of DNA typing in forensic research. Response to Charlier et al

The previously published genetic identification of presumptive samples attributed to two French kings, Henri IV and Louis XVI, by Charlier et al., was recently refuted by a genetic genealogic approach. This (provisional) refutation illustrates the difficulties in confirming the identification of historical DNA samples using limited genetic data. Therefore, we want to stress the necessity of including the genetic genealogic approach--which relies on DNA typing of living relatives of the presumptive donor as a confirmed reference--to validate genetic results in historical cases. Moreover, the popularity and broad media coverage of such studies are useful in bringing awareness to the general public, non-DNA forensic experts and lawyers about the complexity of DNA typing in forensic cases.

Seeing the wood for the trees: a minimal reference phylogeny for the human Y chromosome

During the last few decades, a wealth of studies dedicated to the human Y chromosome and its DNA variation, in particular Y-chromosome single-nucleotide polymorphisms (Y-SNPs), has led to the construction of a well-established Y-chromosome phylogeny. Since the recent advent of new sequencing technologies, the discovery of additional Y-SNPs is exploding and their continuous incorporation in the phylogenetic tree is leading to an ever higher resolution. However, the large and increasing amount of information included in the "complete" Y-chromosome phylogeny, which now already includes many thousands of identified Y-SNPs, can be overwhelming and complicates its understanding as well as the task of selecting suitable markers for genotyping purposes in evolutionary, demographic, anthropological, genealogical, medical, and forensic studies. As a solution, we introduce a concise reference phylogeny whereby we do not aim to provide an exhaustive tree that includes all known Y-SNPs but, rather, a quite stable reference tree aiming for optimal global discrimination capacity based on a strongly reduced set that includes only the most resolving Y-SNPs. Furthermore, with this reference tree, we wish to propose a common standard for Y-marker as well as Y-haplogroup nomenclature. The current version of our tree is based on a core set of 417 branch-defining Y-SNPs and is available online at http://www.phylotree.org/Y.

Genetic genealogy reveals true Y haplogroup of House of Bourbon contradicting recent identification of the presumed remains of two French Kings

Genetic analysis strongly increases the opportunity to identify skeletal remains or other biological samples from historical figures. However, validation of this identification is essential and should be done by DNA typing of living relatives. Based on the similarity of a limited set of Y-STRs, a blood sample and a head were recently identified as those belonging respectively to King Louis XVI and his paternal ancestor King Henry IV. Here, we collected DNA samples from three living males of the House of Bourbon to validate the since then controversial identification of these remains. The three living relatives revealed the Bourbon's Y-chromosomal variant on a high phylogenetic resolution for several members of the lineage between Henry IV and Louis XVI. This 'true' Bourbon's variant is different from the published Y-STR profiles of the blood as well as of the head. The earlier identifications of these samples can therefore not be validated. Moreover, matrilineal genealogical data revealed that the published mtDNA sequence of the head was also different from the one of a series of relatives. This therefore leads to the conclusion that the analyzed samples were not from the French kings. Our study once again demonstrated that in order to realize an accurate genetic identification of historical remains DNA typing of living persons, who are paternally or maternally related with the presumed donor of the samples, is required.

Increasing phylogenetic resolution still informative for Y chromosomal studies on West-European populations

Many Y-chromosomal lineages which are defined in the latest phylogenetic tree of the human Y chromosome by the Y Chromosome Consortium (YCC) in 2008 are distributed in (Western) Europe due to the fact that a large number of phylogeographic studies focus on this area. Therefore, the question arises whether newly discovered polymorphisms on the Y chromosome will still be interesting to study Western Europeans on a population genetic level. To address this question, the West-European region of Flanders (Belgium) was selected as study area since more than 1000 Y chromosomes from this area have previously been genotyped at the highest resolution of the 2008 YCC-tree and coupled to in-depth genealogical data. Based on these data the temporal changes of the population genetic pattern over the last centuries within Flanders were studied and the effects of several past gene flow events were identified. In the present study a set of recently reported novel Y-SNPs were genotyped to further characterize all those Flemish Y chromosomes that belong to haplogroups G, R-M269 and T. Based on this extended Y-SNP set the discrimination power increased drastically as previous large (sub-)haplogroups are now subdivided in several non-marginal groups. Next, the previously observed population structure within Flanders appeared to be the result of different gradients of independent sub-haplogroups. Moreover, for the first time within Flanders a significant East-West gradient was observed in the frequency of two R-M269 lineages, and this gradient is still present when considering the current residence of the DNA donors. Our results thus suggest that an update of the Y-chromosomal tree based on new polymorphisms is still useful to increase the discrimination power based on Y-SNPs and to study population genetic patterns in more detail, even in an already well-studied region such as Western Europe.

Updating the Y-chromosomal phylogenetic tree for forensic applications based on whole genome SNPs

The Y-chromosomal phylogenetic tree has a wide variety of important forensic applications and therefore it needs to be state-of-the-art. Nevertheless, since the last 'official' published tree many publications reported additional Y-chromosomal lineages and other phylogenetic topologies. Therefore, it is difficult for forensic scientists to interpret those reports and use an up-to-date tree and corresponding nomenclature in their daily work. Whole genome sequencing (WGS) data is useful to verify and optimise the current phylogenetic tree for haploid markers. The AMY-tree software is the first open access program which analyses WGS data for Y-chromosomal phylogenetic applications. Here, all published information is collected in a phylogenetic tree and the correctness of this tree is checked based on the first large analysis of 747 WGS samples with AMY-tree. The obtained result is one phylogenetic tree with all peer-reviewed reported Y-SNPs without the observed recurrent and ambiguous mutations. Nevertheless, the results showed that currently only the genomes of a limited set of Y-chromosomal (sub-)haplogroups is available and that many newly reported Y-SNPs based on WGS projects are false positives, even with high sequencing coverage methods. This study demonstrates the usefulness of AMY-tree in the process of checking the quality of the present Y-chromosomal tree and it accentuates the difficulties to enlarge this tree based on only WGS methods.

Genetic genealogy comes of age: perspectives on the use of deep-rooted pedigrees in human population genetics

In this article, we promote the implementation of extensive genealogical data in population genetic studies. Genealogical records can provide valuable information on the origin of DNA donors in a population genetic study, going beyond the commonly collected data such as residence, birthplace, language, and self-reported ethnicity. Recent studies demonstrated that extended genealogical data added to surname analysis can be crucial to detect signals of (past) population stratification and to interpret the population structure in a more objective manner. Moreover, when in-depth pedigree data are combined with haploid markers, it is even possible to disentangle signals of temporal differentiation within a population genetic structure during the last centuries. Obtaining genealogical data for all DNA donors in a population genetic study is a labor-intensive task but the vastly growing (genetic) genealogical databases, due to the broad interest of the public, are making this job more time-efficient if there is a guarantee for sufficient data quality. At the end, we discuss the advantages and pitfalls of using genealogy within sampling campaigns and we provide guidelines for future population genetic studies.

Pig domestication and human-mediated dispersal in western Eurasia revealed through ancient DNA and geometric morphometrics

Zooarcheological evidence suggests that pigs were domesticated in Southwest Asia ∼8,500 BC. They then spread across the Middle and Near East and westward into Europe alongside early agriculturalists. European pigs were either domesticated independently or more likely appeared so as a result of admixture between introduced pigs and European wild boar. As a result, European wild boar mtDNA lineages replaced Near Eastern/Anatolian mtDNA signatures in Europe and subsequently replaced indigenous domestic pig lineages in Anatolia. The specific details of these processes, however, remain unknown. To address questions related to early pig domestication, dispersal, and turnover in the Near East, we analyzed ancient mitochondrial DNA and dental geometric morphometric variation in 393 ancient pig specimens representing 48 archeological sites (from the Pre-Pottery Neolithic to the Medieval period) from Armenia, Cyprus, Georgia, Iran, Syria, and Turkey. Our results reveal the first genetic signatures of early domestic pigs in the Near Eastern Neolithic core zone. We also demonstrate that these early pigs differed genetically from those in western Anatolia that were introduced to Europe during the Neolithic expansion. In addition, we present a significantly more refined chronology for the introduction of European domestic pigs into Asia Minor that took place during the Bronze Age, at least 900 years earlier than previously detected. By the 5th century AD, European signatures completely replaced the endemic lineages possibly coinciding with the widespread demographic and societal changes that occurred during the Anatolian Bronze and Iron Ages.

In the name of the migrant father--analysis of surname origins identifies genetic admixture events undetectable from genealogical records

Patrilineal heritable surnames are widely used to select autochthonous participants for studies on small-scale population genetic patterns owing to the unique link between the surname and a genetic marker, the Y-chromosome (Y-chr). Today, the question arises as to whether the surname origin will be informative on top of in-depth genealogical pedigrees. Admixture events that happened in the period after giving heritable surnames but before the start of genealogical records may be informative about the additional value of the surname origin. In this context, an interesting historical event is the demic migration from French-speaking regions in Northern France to the depopulated and Dutch-speaking region Flanders at the end of the sixteenth century. Y-chr subhaplogroups of individuals with a French/Roman surname that could be associated with this migration event were compared with those of a group with autochthonous Flemish surnames. Although these groups could not be differentiated based on in-depth genealogical data, they were significantly genetically different from each other. Moreover, the observed genetic divergence was related to the differences in the distributions of main Y-subhaplogroups between contemporary populations from Northern France and Flanders. Therefore, these results indicate that the surname origin can be an important feature on top of in-depth genealogical results to select autochthonous participants for a regional population genetic study based on Y-chromosomes.

Automated DNA extraction of single dog hairs without roots for mitochondrial DNA analysis

Dogs are intensely integrated in human social life and their shed hairs can play a major role in forensic investigations. The overall aim of this study was to validate a semi-automated extraction method for mitochondrial DNA analysis of telogenic dog hairs. Extracted DNA was amplified with a 95% success rate from 43 samples using two new experimental designs in which the mitochondrial control region was amplified as a single large (± 1260 bp) amplicon or as two individual amplicons (HV1 and HV2; ± 650 and 350 bp) with tailed-primers. The results prove that the extraction of dog hair mitochondrial DNA can easily be automated to provide sufficient DNA yield for the amplification of a forensically useful long mitochondrial DNA fragment or alternatively two short fragments with minimal loss of sequence in case of degraded samples.

Deep into the roots of the Libyan Tuareg: a genetic survey of their paternal heritage

Recent genetic studies of the Tuareg have begun to uncover the origin of this semi-nomadic northwest African people and their relationship with African populations. For centuries they were caravan traders plying the trade routes between the Mediterranean coast and south-Saharan Africa. Their origin most likely coincides with the fall of the Garamantes who inhabited the Fezzan (Libya) between the 1st millennium BC and the 5th century AD. In this study we report novel data on the Y-chromosome variation in the Libyan Tuareg from Al Awaynat and Tahala, two villages in Fezzan, whose maternal genetic pool was previously characterized. High-resolution investigation of 37 Y-chromosome STR loci and analysis of 35 bi-allelic markers in 47 individuals revealed a predominant northwest African component (E-M81, haplogroup E1b1b1b) which likely originated in the second half of the Holocene in the same ancestral population that contributed to the maternal pool of the Libyan Tuareg. A significant paternal contribution from south-Saharan Africa (E-U175, haplogroup E1b1a8) was also detected, which may likely be due to recent secondary introduction, possibly through slavery practices or fusion between different tribal groups. The difference in haplogroup composition between the villages of Al Awaynat and Tahala suggests that founder effects and drift played a significant role in shaping the genetic pool of the Libyan Tuareg.

Mitochondrial analysis of a Byzantine population reveals the differential impact of multiple historical events in South Anatolia

The archaeological site of Sagalassos is located in Southwest Turkey, in the western part of the Taurus mountain range. Human occupation of its territory is attested from the late 12th millennium BP up to the 13th century AD. By analysing the mtDNA variation in 85 skeletons from Sagalassos dated to the 11th-13th century AD, this study attempts to reconstruct the genetic signature potentially left in this region of Anatolia by the many civilizations, which succeeded one another over the centuries until the mid-Byzantine period (13th century BC). Authentic ancient DNA data were determined from the control region and some SNPs in the coding region of the mtDNA in 53 individuals. Comparative analyses with up to 157 modern populations allowed us to reconstruct the origin of the mid-Byzantine people still dwelling in dispersed hamlets in Sagalassos, and to detect the maternal contribution of their potential ancestors. By integrating the genetic data with historical and archaeological information, we were able to attest in Sagalassos a significant maternal genetic signature of Balkan/Greek populations, as well as ancient Persians and populations from the Italian peninsula. Some contribution from the Levant has been also detected, whereas no contribution from Central Asian population could be ascertained.

Micro-geographic distribution of Y-chromosomal variation in the central-western European region Brabant

One of the future issues in the forensic application of the haploid Y-chromosome (Y-chr) is surveying the distribution of the Y-chr variation on a micro-geographical scale. Studies on such a scale require observing Y-chr variation on a high resolution, high sampling efforts and reliable genealogical data of all DNA-donors. In the current study we optimised this framework by surveying the micro-geographical distribution of the Y-chr variation in the central-western European region named Brabant. The Duchy of Brabant was a historical region in the Low Countries containing three contemporary Belgian provinces and one Dutch province (Noord-Brabant). 477 males from five a priori defined regions within Brabant were selected based on their genealogical ancestry (known pedigree at least before 1800). The Y-haplotypes were determined based on 37 Y-STR loci and the finest possible level of substructuring was defined according to the latest published Y-chr phylogenetic tree. In total, eight Y-haplogroups and 32 different subhaplogroups were observed, whereby 70% of all participants belonged to only four subhaplogroups: R1b1b2a1 (R-U106), R1b1b2a2* (R-P312*), R1b1b2a2g (R-U152) and I1* (I-M253*). Significant micro-geographical differentiation within Brabant was detected between the Dutch (Noord-Brabant) vs. the Flemish regions based on the differences in (sub)haplogroup frequencies but not based on Y-STR variation within the main subhaplogroups. A clear gradient was found with higher frequencies of R1b1b2 (R-M269) chromosomes in the northern vs. southern regions, mainly related to a trend in the frequency of R1b1b2a1 (R-U106).

Mutability of Y-chromosomal microsatellites: rates, characteristics, molecular bases, and forensic implications

Nonrecombining Y-chromosomal microsatellites (Y-STRs) are widely used to infer population histories, discover genealogical relationships, and identify males for criminal justice purposes. Although a key requirement for their application is reliable mutability knowledge, empirical data are only available for a small number of Y-STRs thus far. To rectify this, we analyzed a large number of 186 Y-STR markers in nearly 2000 DNA-confirmed father-son pairs, covering an overall number of 352,999 meiotic transfers. Following confirmation by DNA sequence analysis, the retrieved mutation data were modeled via a Bayesian approach, resulting in mutation rates from 3.78 × 10(-4) (95% credible interval [CI], 1.38 × 10(-5) - 2.02 × 10(-3)) to 7.44 × 10(-2) (95% CI, 6.51 × 10(-2) - 9.09 × 10(-2)) per marker per generation. With the 924 mutations at 120 Y-STR markers, a nonsignificant excess of repeat losses versus gains (1.16:1), as well as a strong and significant excess of single-repeat versus multirepeat changes (25.23:1), was observed. Although the total repeat number influenced Y-STR locus mutability most strongly, repeat complexity, the length in base pairs of the repeated motif, and the father's age also contributed to Y-STR mutability. To exemplify how to practically utilize this knowledge, we analyzed the 13 most mutable Y-STRs in an independent sample set and empirically proved their suitability for distinguishing close and distantly related males. This finding is expected to revolutionize Y-chromosomal applications in forensic biology, from previous male lineage differentiation toward future male individual identification.

Identification and sequence analysis of discordant phenotypes between AmpFlSTR SGM Plus and PowerPlex 16

During duplicate analysis of buccal swabs from 1,377 individuals with 2 commercial short tandem repeat (STR) kits, we observed 8 discordant phenotypes with SGM Plus (SGM, second generation multiplex) for the STRs THO1 (2), vWA (4) and D18S51 (2), and 1 discrepancy with PowerPlex 16 for D18S51. One individual even showed two discrepancies (vWA and THO1) for SGM Plus. In each case, the difference observed was due to the non-amplification or allele dropout of the second allele in a heterozygous genotype. Sequence analysis revealed each time the presence of a mutation that probably coincided with the primer-binding site. Primer-binding site mutations for vWA and D18S51 have been reported previously, while the mutation for THO1 (C-to-T substitution at position 1286 of GenBank sequence D00269) is reported here for the first time. While the frequency of these silent alleles remains low (0.58% in our study), it is suggested that appropriate measures should be taken for database comparisons and that allelic dropout should be further investigated by sequence analysis and be reported to the forensic community.

Y-chromosomal STR haplotypes in a Belgian population sample and identification of a micro-variant with a flanking site mutation at DYS19

Allele frequencies and haplotypes for 12 Y-chromosomal STR loci included in the Powerplex System (Promega, Madison, USA) were determined in a sample of 113 unrelated males of Belgian origin. Ninety-nine different haplotypes were observed with an overall haplotype diversity of 0.997.

Belgian population data for 15 STR loci (AmpFlSTR SGM Plus and AmpFlSTR profiler PCR amplification kit)

The allele and genotype distributions for 15 STR loci included in the AmpFlSTR SGM Plus and AmpFlSTR Profiler kits (Applied Biosystems, Foster City, USA) were determined in a sample of 222 unrelated individuals of Belgian origin.

STR analysis of artificially degraded DNA—results of a collaborative European exercise

Degradation of human DNA extracted from forensic stains is, in most cases, the result of a natural process due to the exposure of the stain samples to the environment. Experiences with degraded DNA from casework samples show that every sample may exhibit different properties in this respect, and that it is difficult to systematically assess the performance of routinely used typing systems for the analysis of degraded DNA samples. Using a batch of artificially degraded DNA with an average fragment size of approx. 200 bp a collaborative exercise was carried out among 38 forensic laboratories from 17 European countries. The results were assessed according to correct allele detection, peak height and balance as well as the occurrence of artefacts. A number of common problems were identified based on these results such as strong peak imbalance in heterozygous genotypes for the larger short tandem repeat (STR) fragments after increased PCR cycle numbers, artefact signals and allelic drop-out. Based on the observations, strategies are discussed to overcome these problems. The strategies include careful balancing of the amount of template DNA and the PCR cycle numbers, the reaction volume and the amount of Taq polymerase. Furthermore, a careful evaluation of the results of the fragment analysis and of automated allele calling is necessary to identify the correct alleles and avoid artefacts.