Twenty-Seven Y-Chromosome Short Tandem Repeats Analysis of Italian Mummies of the 16th and 18th Centuries: An Interdisciplinary Research

Roccapelago (MO) is a small village located in the Northern Central Apennines, with a population of 31 inhabitants (2014). In 2010, more than 400 individuals dated between the end of the 16th and the 18th century, many of which partially mummified, were discovered in the crypt of the church. This small village, because of its geographical location and surrounding environment, seems to possess the characteristics of a genetic isolate, useful for population genetics and genealogical analyses. Thus, a diachronic study of DNA aimed at investigating the structure and dynamics of the population of Roccapelago over the about 4 centuries, was conducted by analyzing ancient and modern inhabitants of the village. The 14 modern samples were selected by considering both the founder surnames of the village, identified thanks to the study of parish registers, and the grandparent’s criterion. From 25 ancient mummies, morphologically assigned to male individuals, the petrous bone, that harbors high DNA amounts, was selected for the DNA extraction. The quantification and qualitative assessment of total human male DNA were evaluated by a real-time PCR assay using the Quantifiler Trio DNA Quantification Kit and multiplex PCR of 27 Y-chromosome short tandem repeat (Y-STR) markers included in the Yfiler Plus PCR Amplification Kit, with seven rapidly mutating Y-STR loci for improving discrimination of male lineages, was performed to genotype the samples. Y-STRs were analyzed according to the criteria of ancient DNA (aDNA) analysis to ensure that authentic DNA typing results were obtained from these ancient samples. The molecular analysis showed the usefulness of the Y chromosome to identify historically relevant remains and discover patterns of relatedness in communities moving from anthropology to genetic genealogy and forensics.

Introduction

Remains of ancient people represent an invaluable opportunity to understand patterns of past population dynamics, appearance and distribution of phenotypes, familial relationships, sex determination, and spread of diseases. They also constitute the opportunity to apply and test protocols, combining and optimizing the archaeogenetic and the forensic ones, in order to achieve and improve results in both fields.

In the context of studies about past populations, the Italian Apennines kept a precious witness of its past. In fact, during the excavation campaign in the Church of the Conversion of Saint Paul, in Roccapelago (Modena, Northern Italy; Figure 1A), a hidden crypt was discovered. The archaeological excavation was conducted between 2009 and 2011, and led to the recovery of the remains of more than 400 individuals from the inside of the crypt, most of which incomplete or commingled.

FIGURE 1

Figure 1. (A) The geographical location of Roccapelago; (B) an image of the current musealization inside the crypt; (C) graphical section of the stratigraphic units (SU) of the Roccapelago crypt. The SU interested by the sampling are labeled with different colors. Modified from the original illustration made by Alessandra Alvisi.

This crypt was used as a cemetery by the inhabitants of the village of Roccapelago between the 16th and the 18th centuries. The significance of this discovery lies in the large number of individuals retrieved that constitute a rare example of the natural conservation of an entire mountain community spanning 2 centuries. Thanks to such peculiarities, these important remains have been the subject of numerous anthropological and paleopathological studies (Traversari et al., 2016, 2019; Figus et al., 2017; Lugli et al., 2017, 2018; Paudice et al., 2021).

In this study, the mummified remains of Roccapelago were collected in order to examine population dynamics and relatedness between centuries in this small mountain village. For its peculiar isolated environment, Roccapelago population seems to possess the features of a genetic isolate, useful for population genetics, biomedical research, and genealogical analyses. The selection of both ancient and present-day samples was guided by anthropological, archaeological, and historical analyses. Contextually to the anthropological study, essential to acquire biodemographic and anthropometric data, it was also possible to analyze the parish registers of the church. The analysis of these documents was of fundamental importance to identify the surname lines present since the 16th century and currently still represented on the territory. Forensic and archaeogenetic protocols were thus applied to the selected individuals for the genetic investigation taking advantage of the increased analytical systems.

Research on ancient and forensic DNA are related since they deal with similar challenges, although, they are guided from different purposes and usually they rely on different methodological approaches (Hofreiter et al., 2020). Despite the progress made in recent years in the development of technologies and protocols, the selection of a suitable source of endogenous DNA is still crucial for the success of genetic analyses on degraded samples (Gonzalez et al., 2020). In particular, in ancient samples it has been demonstrated that the inner ear cochlear region of the otic capsule of the petrous bone, can lead to endogenous DNA yields up to 50–100-fold higher than those obtained from other bones, such as teeth or long bones (Gamba et al., 2014; Pinhasi et al., 2015).

Y-chromosome short tandem repeats (Y-STRs) are widely used for forensic and anthropological applications, like paternal kinship analysis, familial searching, and genetic identification of human remains (Coble et al., 2009; Haas et al., 2013; Kyser, 2017), but also for discerning the origins of migration routes and demographic processes that occurred in historical and prehistoric times (Jobling and Tyler-Smith, 2003; Larmuseau and Ottoni, 2018). The Yfiler Plus multiplex kit consists of 27 STRs including highly polymorphic loci and seven well characterized rapidly mutating markers (RM-YSTRs) with mutation rates upward of 1% (Gopinath et al., 2016), that can provide an increased power of discrimination and improve the level of paternal lineage differentiation (Olofsson et al., 2015; Rapone et al., 2016). Moreover, due to the new primer design and master mix optimization, the Yfiler Plus is designed for efficient amplification of extracted DNA casework samples, like degraded and low-template DNA.

Combining a multidisciplinary approach and following the criteria of ancient DNA (aDNA) analysis, results of Y-STRs markers are presented here with the aim to undertake a diachronic analysis of population genetics of Roccapelago community from 16th century until today. The aims of this study were to predict Y-chromosome paternal haplogroups and to assess genetic continuity between ancient mummies and the current population of Roccapelago. Moreover, this study aimed at testing and combine archaeogenetic and forensic protocols, criteria and approaches, in order to maximize the yield of endogenous DNA, obtain more likely complete Y-STR profiles and minimize exogenous contamination.

Materials and Methods

The Mummies and the Anthropological Analyses

The large amount of human remains retrieved in the crypt, mostly consisted of commingled and completely skeletonized human remains, comprised some mummified corpses, especially among the most recent stratigraphic units (SU; Figure 1B). Several remains were partially mummified, as a result of a natural process due to the particular microclimatic conditions of the environment. The bodies, piled in a pyramid shape, were divided into three SU: SU 28 has returned very few individuals and was dated between the end of the 16th century and the beginning of the 17th century; SU 26 was placed at the second half of the 17th and the beginning of the 18th century; SU 23, the most recent, dated at the 18th century (Figure 1C). The chronological division was made through the study of the clothing and grave goods. Sex estimation was performed by observing morphological characteristics of the skull (Acsádi and Nemeskéri, 1970). Only individuals classified as males were collected for the purposes of this study. For Y-STR analysis, we chose 25 individuals, entirely or partially skeletonized, which, among the other skeletal elements, presented the skull intact or partially preserved, but containing the portion of pars petrosa. They included two samples from SU 28, 11 from SU 26, and 12 from SU 23 (Table 1). For assessment of age at death, we combined different methods, obtaining the smallest range possible: the dental attrition was evaluated by approach of Lovejoy (1985); methods of Acsádi and Nemeskéri (1970) and Meindl and Lovejoy (1985) were used for judging cranial sutures degree of closure.

TABLE 1

Table 1. List of ancient and modern samples analyzed in this study.

The Selection of Modern Samples

The 14 modern samples, representing about half of the present-day population of Roccapelago, were selected by considering both the founder surnames of the village, identified thanks to the study of parish registers (that became mandatory in Italy for baptisms and marriages in 1563 after the Council of Trent and in 1614 for deaths when their rules of compilation were as well regulated by the Church), and the grandparent’s criterion.

Every single death certificate was recorded on a table in chronological order. Subsequently, the individual surnames were isolated and analyzed in terms of continuity and frequency; finally, it was possible to extract only the surnames that have demonstrated a continuous presence in the area. Collaterally, the analysis highlighted some cognominal lines that have met with real extinctions on the territory, a phenomenon perhaps due to migrations of entire families.

Modern individuals were sampled by means of cotton buccal swabs, after the collection of informed consent. The study was approved by the Bioethical Committee of the University of Bologna (IRB approval date: April 8th, 2013) and was conducted in agreement with the principles of the WMA Declaration of Helsinki.

DNA Extraction

DNA extraction of ancient specimens was carried out in clean-rooms reserved for the analysis of degraded DNA at the Laboratory of ancient DNA (aDNALab) of the Department of Cultural Heritage, University of Bologna, Ravenna Campus. This pre-PCR area of the laboratory is physically isolated and is equipped with positive air pressure with HEPA filters and laminar flow cabinets, reserved for the different phases of the work. All steps were conducted under strict guidelines for the monitoring of the exogenous contamination and reproducibility of data (Cooper and Poinar, 2000; Gilbert et al., 2005; Llamas et al., 2017). Blanks as negative controls were processed in parallel with the samples, in order to monitor for the presence of contaminations.

In order to maximize the results of the study, the sampling was directed at the inner ear cochlear region of the petrous bone that provides good percentages of endogenous DNA (Gamba et al., 2014; Pinhasi et al., 2015). The sampling was conducted on 24 complete skulls and one isolated petrous bone (sample 26-5). For complete skulls, sampling was performed using the Cranial Base Drilling Method (CBDM; Sirak et al., 2017) in order to minimize structural damages to the cranial features. The isolated petrous bones were directly sampled on the dense parts around the cochlea (Orfanou et al., 2021). Between 40 and 91mg of bone powder were collected from each individual (Table 1) and processed with a silica-based protocol (Dabney et al., 2013; Damgaard et al., 2015), slightly modified as in Cilli et al. (2020). Samples were predigested in 1ml of Extraction buffer (0.5M EDTA, 0.25mg/ml Proteinase K, 0.05% Tween20) for 15min at 37°C to increase the endogenous DNA content and then digested overnight in a new aliquot of 1ml of the same buffer. Supernatant was centrifuged and DNA bound into silica columns with PB buffer (Qiagen). Silica columns were finally washed twice with PE buffer (Qiagen) and eluted in 50μl of elution buffer.

The DNA extraction of modern samples was performed by using the QIAmp DNA Mini kit (Qiagen), following manufacturer instructions.

DNA Quantification and Y-STR Genotyping

After the extraction, all samples were quantified with the Quantifiler® Trio DNA Quantification Kit (Applied Biosystems) following the manufacturer’s instructions (Vernarecci et al., 2015). Quantitative PCR (qPCR) was performed in double using QuantStudio 5 Real-Time PCR System (Applied Biosystems) and data were analyzed using the HID Real-Time PCR Analysis Software v1.3 using the default settings. The quantification data obtained from the small autosomal probe of the Quantifiler Trio Kit was used for determining the concentration of aDNA.

Y-chromosome short tandem repeats amplification was performed with the Yfiler® Plus PCR Amplification Kit (Applied Biosystems), including the 27 STRs loci DYS19, DYS385a/b, DYS389I/II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, DYS448, DYS456, DYS458, DYS635, YGATAH4, DYS460, DYS481, DYS533, DYF387S1a/b, DYS449, DYS518, DYS570, DYS576, and DYS627, on a Veriti® 96 Well Thermal Cycler System (Applied Biosystems) and following the manufacturer’s instructions. PCR products were electrophoresed on a SeqStudio Genetic Analyzer (Applied Biosystems) and the fragment analysis was performed by GeneMapper® ID-X Software v1.6 (Applied Biosystems). Only peaks above 50 RFU were considered and replicates were performed for all ancient samples.

Statistical Analyses

Allele frequencies, Y-chromosome diversity indexes and between-population pairwise distances were calculated with the Arlequin software version 3.5 (Excoffier and Lischer, 2010). In particular, the gene diversity (GD) for each locus was calculated using the formula GD=n(1−∑x_i²)/(n−1), with x_i being the frequency of the ith allele and n the sample size of each considered population. The intra-population haplotype diversity (HD) was estimated as HD=n(1−∑p_i²)/(n−1), where n is the number of analyzed samples and p_i the frequency of the ith haplotype. The haplotype match probability (HMP) was calculated using the formula HMP=∑p_i², whereas the discrimination capacity (DC) was obtained by dividing the total number of observed haplotypes for the total number of individuals in each dataset.

The obtained Y-chromosome STR haplotypes were used to infer haplogroups of individual samples by using the Whit Athey’s Haplogroup Predictor software (Athey, 2006) and the Nevgen prediction tool,¹ while the genetic relationships between Y-STR haplotypes were explored by means of a Median Joining (MJ) Network analysis using the Network 10.2 program.²

For all statistical analyses, the locus DYS389II was converted to the DYS389b nomenclature by subtracting the repeat number of the DYS389I from the detected DYS389II allele, and each of both the DYS385 and DYF387S1 multi-locus markers were treated as two separated allelic series.

In order to investigate the possible genetic relationships between ancient and modern individuals and assess genetic continuity in the population of Roccapelago over about 4 centuries, we estimated the Time to the Most Recent Common Ancestor (TMRCA) for each ancient-modern pair by adopting the Walsh (2001) method as implemented in the R script from Boattini et al. (2019). For each tested pair, the input parameters required by the script are (i) the mutation model (either the stepwise mutation model – SMM or the infinite allele model – IAM, originally implemented by Walsh, 2001); (ii) the average mutation rate for the set of considered Y-STRs markers; (iii) the lambda prior hyper-parameter (Walsh, 2001); and (iv) a vector of per couple differences for the analyzed loci (i.e., a file indicating matching “0” and non-matching “1” Y-STR alleles for the IAM model, or the number of repeat differences for each locus in the case of SMM model). Data output for each tested pair returns the posterior distribution of the estimated TMRCA with mean, median, and mode values as well as the computed confidence interval(s). Comparisons between documented genealogical pairs and the TMRCA estimates obtained with the Walsh method in Boattini et al. (2019), showed that the IAM model yielded the best performance compared to SMM. In addition, among the considered summary statistics, the modal value of the posterior distributions was found to be the one best approximating the observed values. Accordingly, the IAM model was therefore chosen also in the present analysis, while the Y-STR mutation rates adopted in the procedure were taken from Ballantyne et al. (2010). Since the Walsh method implemented in Boattini et al. (2019) was originally designed for estimating TMRCA of genealogical pairs of individuals whose life span overlaps, in the present study, where pairs are instead composed by chronologically distant pairs of ancient and modern individuals, we adapted TMRCA to calculate the total number of generations separating them, which corresponds to the double of the originally provided estimate. Finally, the obtained numbers of generations were converted in years before present by using an average generation time of 33years, in agreement with genealogical-based estimates obtained from previous studies on the same geographic area (Boattini et al., 2014, 2019; Sarno et al., 2021).

Results

Authenticity of the Results and Success Rate of Ancient DNA Analysis

The strict criteria followed in this study allowed us to exclude any modern DNA contamination and confirm the reliability of the genotyping results. No contamination was observed in any of the blank extractions or negative PCR controls included in each reaction. Moreover, the data were consistent between replicates of amplifications. In addition, the haplotypes of the samples were different from those of the anthropologists and geneticists involved in the sampling and DNA analyses.

DNA extraction from petrous bone yielded endogenous DNA by qPCR analysis ranging from 0.004 to 0.226ng/μl for the small autosomal target and from 0.001 to 0.175ng/μl for the Y-target. The range of the degradation index (DI) spanned between 2 and 107, except one sample for which the DI could not be determined as the large autosomal target was not detected (Table 2). Furthermore, for two samples (26-12 and 26-53) the Y-target was not determined, and therefore they were not included in the amplification session.

TABLE 2

Table 2. Quantification by real-time PCR and Y-chromosome short tandem repeat (Y-STRs) amplification results from the 25 mummies analyzed.

Yfiler Plus STRs markers were successfully amplified for at least 10 loci in 19 out of 25 ancient DNA samples and, as expected, in all the 14 analyzed modern individuals from Roccapelago (Table 2). In particular, only five of the analyzed ancient specimens (namely samples 23-7, 23-31, 23-46, 26-05, and 26-32) provided complete results for all the 27 loci of the Yfiler Plus; for the remaining genotyped mummies partial Y-chromosome haplotypes were instead obtained consisting of 25 (samples 23-06, 23-41, 26-1, and 26-10), 24 (samples 23-12, 23-44, 26-42, and 26-43), 21 (samples 23-42, 26-50), 16 (samples 23-01, 28-05), 14 (sample 26-49), and 10 (sample 28-02) STR markers, respectively. Not surprisingly, the two human remains from the oldest stratigraphic unit (i.e., SU28) showed the lowest number of genotyped Y-STRs loci. The remaining amplified ancient samples provided less than nine Y-STRs loci and therefore were not considered for the analyses.

All samples provided results for loci DYS385, DYS393, DYS456, DYS458, DYS460, DYS570, and DYS576, while the worst performance was obtained for the largest loci DYS391 and the DYS533. We also observed a duplication (alleles 14 and 15) at the locus DYS19 in sample 26-10 confirmed by PCR replicates. In Supplementary Table S1, all observed Y-STRs haplotypes in ancient and modern individual sample are provided. When comparing haplotypes of present-day individuals sharing the same surnames, one step mutations in two groups of modern samples for three of seven RM Y-STRs included in the kit were found, in particular for DYS627 in two individuals and for DYS570 and DYS518 in four modern samples.

Haplotype Diversity Analysis and Forensic Indexes

Being aware that the presence of some missing loci in ancient individuals may affect the power of discrimination in distinguishing haplotypes, in order to maximize as much as possible the level of male lineage differentiation we considered for statistical analyses only individuals with no more than three missing genotyped loci, thus resulting in 13 ancient and all 14 modern samples left for subsequent computations. On the whole, using the Y-STR Yfiler Plus Kit a total of 20 unique haplotypes were observed among the 27 considered individuals of Roccapelago (Table 3), with two haplotypes instead occurring twice and another haplotype shared among three samples. In particular, two modern individuals (4A and 10A, who also share the same surname), were found to share a complete Yfiler Plus profile. The ancient mummy 23-44 instead matches for its 24 genotyped loci the ones of the modern individual 7A. Finally, also the ancient sample 23-7 and the modern individual 17A shared an identical complete Yfiler Plus haplotype, which also matches the one of the mummy 23-06 for its 25 genotyped loci.

TABLE 3

Table 3. Diversity indexes for the whole analyzed samples from Roccapelago, and for the ancient (RP_A) and modern (RP_M) population sets, separately.

Haplotype diversity and DC for the Yfiler Plus Kit were found to be 0.9858 and 0.8519, respectively in the ancient-plus-modern dataset. Comparable, if not higher values were observed also when the 13 ancient (RP_A: HD=0.9872 and DC=0.9231) and the 14 modern (RP_M: HD=0.9890 and DC=0.9286) individuals were considered as two separated population sets. In all of the cases, the use of the Yfiler Plus Kit was proven to significantly increase the power of discrimination compared to the original Yfiler panel (Table 3).

In particular, the 10 Y-STRs loci added to the 17-plex Yfiler kit, and especially the seven newly-included RM loci, were indeed the ones showing on average the higher values of gene diversity (GD>0.70) in all the considered population sets (Table 4)

TABLE 4

Table 4. Genetic diversity for each Y-STR marker in the considered datasets (RP_T: whole samples from Roccapelago; RP_A: ancient mummies; and RP_M: modern individuals).

Population Relationship Analyses

Y-chromosome short tandem repeat haplotypes were used to predict Y-chromosome paternal haplogroups and to assess genetic continuity between ancient mummies and the current population of Roccapelago.

Y-chromosome haplogroups were successfully inferred for all samples, with high prediction probabilities in almost all of the cases. Lower probabilities were observed for those samples showing more than 10 missing loci (Supplementary Table S1). Ancient DNA samples were found to belong to the Y-chromosome macro lineages I2a, J1, J2b, and R1b. The same haplogroups also characterize the present-day individuals, suggesting a homogenous Y-chromosome composition with respect to the current population of Roccapelago, even when using a Y haplogroup predictor with a higher level of resolution. Accordingly, the degree of differentiation between RP_A and RP_M population sets, evaluated based on both Y-chromosome haplogroup frequencies and Y-STR RSTs genetic distances, resulted low and not significant (FST: 0.00381, p: 0.41877 and RST: 0.01730, p: 0.28433, respectively).

The phylogenetic relationships between modern and ancient haplotypes (including only individuals showing no more than three missing loci as above) were visually represented by means of a MJ Network (Figure 2). As expected, haplotypes clustered based on corresponding haplogroup lineage, also showing little differences in STR profiles between the ancient and modern individuals of Roccapelago within the same haplogroups. In particular, the haplotype shared between the modern 4A and 10A samples, also shows a single different locus with respect to the 25 markers successfully typed in the ancient mummy 23-41. Analogously, all the six J1 ancient samples (23-31, 23-12, 23-46, 26-05, 26-32, 26-43) appear one-step or two-steps neighbors of the modern J1 sample 16. The ancient J2b mummies 23-7 and 23-06 (successfully typed for all the 27 and for 25 loci, respectively) share an identical Y-chromosome haplotype with the modern individual 17A, as outlined also above by the haplotype sharing analysis. Similarly, the 24 genotyped loci of the ancient mummy 23-44 match the haplotype of the modern individual 7A, further showing one or at most two different loci also with respect to the other modern R1b individuals who share the same surname of 7A (namely 15A, 18, and 21A).

FIGURE 2

Figure 2. Median Joining (MJ) network analysis of the Y-chromosome haplotypes based on the 27 Yfiler Plus STRs in the 13 ancient and 14 modern samples from Roccapelago. Present-day individuals are highlighted in white dots, ancient mummies from the stratigraphic units SU23 in light-gray and those from the SU26 in dark-gray.

In order to formally assess the presence of genetic connections among the ancient mummies and the present-day inhabitants of Roccapelago, we estimated the expected TMRCA between ancient-modern pairs (Supplementary Table S2) by exploiting the Bayesian procedure originally designed by Walsh (2001) and implemented by Boattini et al. (2019), as detailed in Materials and Methods. In particular, we looked for pairs for which the CI of the inferred values fall within the chronological range of dating of the ancient mummies, by specifically retaining those (highlighted in bold in the upper panel of Figure 3) whose more stringent 50% CI overlaps with the time window of 300–500years ago. Accordingly, we observed that 11 of 13 mummies indeed showed signs of continuity with modern inhabitants of Roccapelago (Figure 3), thus suggesting that the genetic composition of the population did not significantly change with time. Interestingly, significant TMRCA results emerged for the same couples also highlighted in the Network analysis (i.e., 23-41 with 4A/10A; 23-12, 23-31, 23-46, 26-05, 26-32, and 26-43 with 16A; 23-06 and 23-7 with 17A; and 23-44 with 7A, 21A, 18, and 15A).

FIGURE 3

Figure 3. Barplots representing confidence intervals (black: 50% CI; gray: 95% CI) of estimates of times separating ancient and modern samples. Vertical dashed lines (in red) delineate the chronological window of Roccapelago mummies (~300–500years ago). The upper panel shows pairs with best-matching results (i.e., overlap between mummies chronological window and the 50% CIs of time estimates).

Discussion

In forensic and archaeological applications, DNA degradation represents a challenge during each phase of the analyses, and methodological improvements are necessary to obtain successful results from heavily degraded and low template DNA samples like mummies. Certainly, in the last years, the scientific progress and technological development contributed to overcome the difficulties to obtain more likely complete STR profiles for genetic identification, and at the same time the standard molecular biology approaches have been widely improved to address issues of DNA degradation, contamination, and inhibition (Malmström et al., 2005; Votrubova-Dubska et al., 2016).

Among the methodological innovations that have made possible to overcome these difficulties there is the identification of the best skeletal elements, of which in primis petrous bones (Edson et al., 2009; Gamba et al., 2014; Pinhasi et al., 2015). Even though, this skeletal element can lead to endogenous DNA yields that are up to 50–100-fold higher than those obtained from other bones, recent studies demonstrated that other samples (in situ molars, thoracic vertebrae, talus, and distal phalanx) can yield adequate DNA for most applications in human palaeogenetics (Higgins et al., 2013; Hansen et al., 2017; Parker et al., 2020), especially ear ossicles (Sirak et al., 2020). Important methodological improvements also regarded the development of increasingly efficient extraction protocols and genotyping kits. The petrous bone (the densest part of the temporal bone) is a source of DNA in both larger quantities and of better quality than other bones, and it represents the substrate of choice for DNA extraction and analyses of ancient and/or degraded DNA (Gonzalez et al., 2020). Although, it was originally published in the forensic field (Edson et al., 2009), the petrous bone was not truly implemented in laboratory work by the forensic community, probably due to ethical aspects related to the invasive sampling (Hofreiter et al., 2020). Recently, several studies have however demonstrated the usefulness of this skeletal element to genotype forensic samples (Pilli et al., 2018; Gonzalez et al., 2020), and the improving Y-STR typing for human identification by the YFiler Plus kit in skeletal remains (Ambers et al., 2018). In addition, the increased efficiency of extraction protocols and amplification systems should also be noted (Loreille et al., 2007; Amory et al., 2012; Dabney et al., 2013; Hofreiter et al., 2020).

This study confirmed the value of the otic capsule of the petrous bone as a skeletal element in maximizing DNA yield in order to increase the probability of obtaining complete Y-STR profiles and limit allelic drop-out, but also in preventing contamination with exogenous human DNA (Gonzalez et al., 2020). The genetic material showed different rates of preservation, ranging from a moderate to a significant degradation (reaching 107), as indicated from the values of the DI reported in Table 2. In only one sample the large autosomal target was not detected. Although, it is well-known that physical, biological, and chemical factors are deeply involved in the long-term preservation of DNA (Collins et al., 2002), in our results the typing success reached 76% and PCR negative samples were related to the small detected quantity of Y DNA. For the two samples completely negative for the Y-target quantification, the results could also be attributed to a not certain morphological sex attribution since a minimal quantity of autosomal DNA was detected. Moreover, in ancient samples showing similar DNA quantity, a strong correlation between DI values and typing success was not observed.

Specifically, 19 out of 25 ancient DNA samples showed more than 10 loci typed. A total of five specimens provided a complete profile for all the 27 loci of the Yfiler Plus kit, while for the other samples it was possible to obtain a number of amplified loci from 25 to 10, where the worst results were obtained from the oldest stratigraphic unit (i.e., SU28). Regarding the STR markers, our typing data suggest that for the analysis of skeletal remains a better performance can be obtained for the DYS385a/b, DYS393, DYS456, DYS458, DYS460, DYS570, and DYS576 loci but not for DYS391 and DYS533 loci, possibly due to the new primer design resulting in a larger amplicon size in the YFiler Plus kit compared to Yfiler kit. Indeed, although, a multiplex system with a greater number of markers increases discriminatory power, often aDNA is limited in quantity and compromised quality resulting in reduction or loss of signal especially for larger loci.

Among the modern individuals sharing the same surnames a one-step mutation was observed in at least one RM Y-STR (DYS570, DYS627, and DYS518). This finding supports again the potentially increased mutability of these markers, which would ideally become the first choice for Y-STR testing in forensic casework (excluding family testing) due to their superior value in paternal lineage differentiation as well as male relative separation (Ballantyne et al., 2012; Boattini et al., 2016). The data also confirmed that the YFiler Plus kit offers higher power of discrimination and enhanced chemistry for improved performance with challenging samples, as well as a higher resolution than the YFiler panel, especially thanks to the seven newly-included RM loci that on average showed the higher values of gene diversity (Table 4). This agrees with previously published data (Olofsson et al., 2015; Rapone et al., 2016; Fan et al., 2019). The Yfiler Plus ability to discriminate among closely related males could have useful applications especially in situations of multiple victims from the same family buried together, like a mass grave, or mummies discovered in a crypt (Ambers et al., 2018), by however considering the possibility of false exclusions due to mutation when compared to living relatives in kinship testing.

As concern population relationships analyses, overall our results support a paternal genetic continuity of the analyzed population from the 16th century until today. In general, the predicted Y-chromosome haplogroups reflect lineages commonly observed in Italy (Grugni et al., 2018) and revealed a homogeneous genetic composition through centuries between the ancient and modern samples of Roccapelago, as supported even by both Y-chromosome haplogroup frequencies and RSTs genetic distances. The connection between ancient and modern individuals of Roccapelago is well documented also in the network analysis (Figure 2), where both mummies’ and present-day haplotypes clustered together based on haplogroups assignment. Additionally, searching against the Y-haplotype reference database (YHRD) resulted in non-matches between each ancient sample and any of the other Yfiler Plus reference samples included in the database (encompassing also other Italian groups), thus providing further evidence of a closer proximity between the mummies’ haplotypes and the modern individuals of Roccapelago. With the aim to formally test the genetic connection among the ancient mummies and the present-day inhabitants of Roccapelago, we finally estimated the expected TMRCA between them. This analysis, supported by the network results, confirmed signals of genetic continuity between ancient and modern inhabitants of Roccapelago for 11 out of the 13 mummies here considered.

These findings agree with the historical data and the partial past geographical isolation of this small village. In fact, from the analysis of parish registers it has been highlighted, during centuries, a trend toward exogamic marriages, with external male contributions, mainly from the areas interested by transhumance. However, for the founding surnames, like the ones sampled for this study, a tendency to isonymia and consanguinity has been noted from the parish registers. This probably was due to a social-cultural practice to pass down lands through generations between the same ancient founder families (Traversari, 2020).

In conclusion, a multidisciplinary approach was successfully applied to the mummies of a small mountain community dated to 16th–18th centuries, whose Y-STR haplotypes were analyzed and compared to data obtained from the modern people of the village. We combined forensic genetics and archaeogenetics approaches and protocols with historical, archaeological, and anthropological data, in order to understand population dynamics and relatedness on this peculiar community from the 16th century. This study highlights the need of the contribution from different disciplines to a better and more efficient human identification.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Ethics Statement

The studies involving human participants were reviewed and approved by Bioethical Committee of the University of Bologna (approval date: April 8th, 2013). The patients/participants provided their written informed consent to participate in this study.

Author Contributions

DL and SP conceived the study. EC, MT, and CB designed the experiments. CB, FF, and EC performed the genetic analyses. SS, AB, CB, and EC analyzed the data. MT performed the archaeological, anthropological, and historical analyses. EC, SS, and CB wrote the original draft. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the MIUR-PRIN 20177PJ9XF grant to DL.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We are grateful to all the donors who kindly provided their biological samples making this study possible. We would like to thank Sara Amurri e Martina Montagnoli who assisted during the experimental procedures. Our gratitude also goes to the “Soprintendenza Archeologia, Belle Arti e Paesaggio per la città metropolitana di Bologna e le province di Modena, Reggio Emilia e Ferrara” for allowing us to perform the analysis on these samples.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.720640/full#supplementary-material

Footnotes

1. ^https://www.nevgen.org

2. ^http://www.fluxus-engineering.com

References

Acsádi, G., and Nemeskéri, J. (1970). Determination of Sex and Age at Death From Skeletal Finds in History of Human Life. Span and Mortality. Budapest: Akadémiai Kiadό.

Google Scholar

Ambers, A., Votrubova, J., Vanek, D., Sajantila, A., and Budowle, B. (2018). Improved Y-STR typing for disaster victim identification, missing persons investigations, and historical human skeletal remains. Int. J. Legal Med. 132, 1545–1553. doi: 10.1007/s00414-018-1794-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Amory, S., Huel, R., Bilić, A., Loreille, O., and Parsons, T. J. (2012). Automatable full demineralization DNA extraction procedure from degraded skeletal remains. Forensic Sci. Int. Genet. 6, 398–406. doi: 10.1016/j.fsigen.2011.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Athey, T. W. (2006). Haplogroup prediction from Y-STR values using a Bayesian-allele frequency approach. J. Genet. Gen. 2, 34–39.

Google Scholar

Ballantyne, K. N., Goedbloed, M., Fang, R., Schaap, O., Lao, O., Wollstein, A., et al. (2010). Mutability of Y-chromosomal microsatellites: rates, characteristics, molecular bases, and forensic implications. Am. J. Hum. Genet. 10, 341–353. doi: 10.1016/j.ajhg.2010.08.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Ballantyne, K. N., Keerl, V., Wollstein, A., Choi, Y., Zuniga, S. B., Ralf, A., et al. (2012). A new future of forensic Y-chromosome analysis: rapidly mutating Y-STRs for differentiating male relatives and paternal lineages. Forensic Sci. Int. Genet. 6, 208–218. doi: 10.1016/j.fsigen.2011.04.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Boattini, A., Sarno, S., Bini, C., Pesci, V., Barbieri, C., De Fanti, S., et al. (2016). Mutation rates and discriminating power for 13 rapidly-mutating Y-STRs between related and unrelated individuals. PLoS One 11:e0165678. doi: 10.1371/journal.pone.0165678

PubMed Abstract | CrossRef Full Text | Google Scholar

Boattini, A., Sarno, S., Mazzarisi, A. M., Viroli, C., De Fanti, S., Bini, C., et al. (2019). Estimating Y-Str mutation rates and Tmrca through deep-rooting Italian pedigrees. Sci. Rep. 9:9032. doi: 10.1038/s41598-019-45398-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Boattini, A., Sarno, S., Pedrini, P., Medoro, C., Carta, M., Tucci, S., et al. (2014). Traces of medieval migrations in a socially stratified population from northern Italy. Evidence from uniparental markers and deep-rooted pedigrees. Heredity 114, 155–162. doi: 10.1038/hdy.2014.77

PubMed Abstract | CrossRef Full Text | Google Scholar

Cilli, E., Gabanini, G., Ciucani, M. M., De Fanti, S., Serventi, P., Bazaj, A., et al. (2020). A multifaceted approach towards investigating childbirth deaths in double burials: anthropology, paleopathology and ancient DNA. J. Archaeol. Sci. 122:105219. doi: 10.1016/j.jas.2020.105219

CrossRef Full Text | Google Scholar

Coble, M. D., Loreille, O. M., Wadhams, M. J., Edson, S. M., Maynard, K., Meyer, C. E., et al. (2009). Mystery solved: the identification of the two missing Romanov children using DNA analysis. PLoS One 4:e4838. doi: 10.1371/journal.pone.0004838

PubMed Abstract | CrossRef Full Text | Google Scholar

Collins, M. J., Nielsen-Marsh, C. M., Hiller, J., Smith, C. I., Roberts, J. P., Prigodich, R. V., et al. (2002). The survival of organic matter in bone: a review. Archaeometry 44, 383–394. doi: 10.1111/1475-4754.t01-1-00071

CrossRef Full Text | Google Scholar

Cooper, A., and Poinar, H. N. (2000). Ancient DNA: do it right or not at all. Science 289:1139. doi: 10.1126/science.289.5482.1139b

PubMed Abstract | CrossRef Full Text | Google Scholar

Dabney, J., Knapp, M., Glocke, I., Gansauge, M. T., Weihmann, A., Nickel, B., et al. (2013). Complete mitochondrial genome sequence of a middle Pleistocene cave bear reconstructed from ultrashort DNA fragments. Proc. Natl. Acad. Sci. U. S. A. 110, 15758–15763. doi: 10.1073/pnas.1314445110

PubMed Abstract | CrossRef Full Text | Google Scholar

Damgaard, P. B., Margaryan, A., Schroeder, H., Orlando, L., Willerslev, E., and Allentoft, M. E. (2015). Improving access to endogenous DNA in ancient bones and teeth. Sci. Rep. 5:11184. doi: 10.1038/srep11184

PubMed Abstract | CrossRef Full Text | Google Scholar

Edson, S. M., Christensen, A. F., Barritt, S. M., Meehan, A., Leney, M. D., and Finelli, L. N. (2009). Sampling of the cranium for mitochondrial DNA analysis of human skeletal remains. Forensic Sci. Int. Genet. Suppl. Ser. 2, 269–270. doi: 10.1016/j.fsigss.2009.09.029

CrossRef Full Text | Google Scholar

Excoffier, L., and Lischer, H. (2010). Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and windows. Mol. Ecol. Resour. 10, 564–567. doi: 10.1111/j.1755-0998.2010.02847.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, G. Y., An, Y. R., Peng, C. X., Deng, J. L., Pan, L. P., and Ye, Y. (2019). Forensic and phylogenetic analyses among three Yi populations in Southwest China with 27 Y chromosomal STR loci. Int. J. Legal Med. 133, 795–797. doi: 10.1007/s00414-018-1984-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Figus, C., Traversari, M., Scalise, L. M., Oxilia, G., Vazzana, A., Buti, L., et al. (2017). The study of commingled non-adult human remains: insights from the 16th–18th centuries community of Roccapelago (Italy). J. Archaeol. Sci. Rep. 14, 382–391. doi: 10.1016/j.jasrep.2017.06.023

CrossRef Full Text | Google Scholar

Gamba, C., Jones, E. R., Teasdale, M. D., McLaughlin, R. L., Gonzalez-Fortes, G., Mattiangeli, V., et al. (2014). Genome flux and stasis in a five millennium transect of European prehistory. Nat. Commun. 5:5257. doi: 10.1038/ncomms6257

PubMed Abstract | CrossRef Full Text | Google Scholar

Gilbert, M. T. P., Bandelt, H.-J., Hofreiter, M., and Barnes, I. (2005). Assessing ancient DNA studies. Trends Ecol. Evol. 20, 541–544. doi: 10.1016/j.tree.2005.07.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Gonzalez, A., Cannet, C., Zvénigorosky, V., Geraut, A., Koch, G., Delabarde, T., et al. (2020). The petrous bone: ideal substrate in legal medicine? Forensic Sci. Int. Genet. 47:102305. doi: 10.1016/j.fsigen.2020.102305

PubMed Abstract | CrossRef Full Text | Google Scholar

Gopinath, S., Zhong, C., Nguyen, V., Ge, J., Lagacé, R. E., Short, M. L., et al. (2016). Developmental validation of the Yfiler(®) plus PCR amplification kit: an enhanced Y-STR multiplex for casework and database applications. Forensic Sci. Int. Genet. 24, 164–175. doi: 10.1016/j.fsigen.2016.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Grugni, V., Raveane, A., Mattioli, F., Battaglia, V., Sala, C., Toniolo, D., et al. (2018). Reconstructing the genetic history of Italians: new insights from a male (Y-chromosome) perspective. Ann. Hum. Biol. 45, 44–56. doi: 10.1080/03014460.2017.1409801

PubMed Abstract | CrossRef Full Text | Google Scholar

Haas, C., Shved, N., Rühli, F. J., Papageorgopoulou, C., Purps, J., Geppert, M., et al. (2013). Y-chromosomal analysis identifies the skeletal remains of Swiss national hero Jörg Jenatsch (1596–1639). Forensic Sci. Int. Genet. 7, 610–617. doi: 10.1016/j.fsigen.2013.08.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Hansen, H. B., Damgaard, P. B., Margaryan, A., Stenderup, J., Lynnerup, N., Willerslev, E., et al. (2017). Comparing ancient DNA preservation in petrous bone and tooth cementum. PLoS One 12:e0170940. doi: 10.1371/journal.pone.0170940

PubMed Abstract | CrossRef Full Text | Google Scholar

Higgins, D., Kaidonis, J., Townsend, G., Hughes, T., and Austin, J. J. (2013). Targeted sampling of cementum for recovery of nuclear DNA from human teeth and the impact of common decontamination measures. Investig. Genet. 4:18. doi: 10.1186/2041-2223-4-18

PubMed Abstract | CrossRef Full Text | Google Scholar

Hofreiter, M., Sneberger, J., Pospisek, M., and Vanek, D. (2020). Progress in forensic bone DNA analysis: lessons learned from ancient DNA. Forensic Sci. Int. Genet. 54:102538. doi: 10.1016/j.fsigen.2021.102538

PubMed Abstract | CrossRef Full Text | Google Scholar

Jobling, M. A., and Tyler-Smith, C. (2003). The human Y chromosome: an evolutionary marker comes of age. Nat. Rev. Genet. 4, 598–612. doi: 10.1038/nrg1124

PubMed Abstract | CrossRef Full Text | Google Scholar

Kyser, M. (2017). Forensic use of Y-chromosome DNA: a general overview. Hum. Genet. 136, 621–635. doi: 10.1007/s00439-017-1776-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Larmuseau, M. H. D., and Ottoni, C. (2018). Mediterranean Y-chromosome 2.0-why the Y in the Mediterranean is still relevant in the postgenomic era. Ann. Hum. Biol. 45, 20–33. doi: 10.1080/03014460.2017.1402956

PubMed Abstract | CrossRef Full Text | Google Scholar

Llamas, B., Valverde, G., Fehren-Schmitz, L., Weyrich, L. S., Cooper, A., and Haak, W. (2017). From the field to the laboratory: controlling DNA contamination in human ancient DNA research in the high-throughput sequencing era. Sci. Technol. Archaeol. Res. 3, 1–14. doi: 10.1080/20548923.2016.1258824

CrossRef Full Text | Google Scholar

Loreille, O. M., Diegoli, T. M., Irwin, J. A., Coble, M. D., and Parsons, T. J. (2007). High efficiency DNA extraction from bone by total demineralization. Forensic Sci. Int. Genet. 1, 191–195. doi: 10.1016/j.fsigen.2007.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Lovejoy, C. O. (1985). Dental wear in the Libben population: its functional pattern and role in the determination of adult skeletal age at death. Am. J. Phys. Anthropol. 68, 47–56. doi: 10.1002/ajpa.1330680105

PubMed Abstract | CrossRef Full Text | Google Scholar

Lugli, F., Brunelli, D., Cipriani, A., Bosi, G., Traversari, M., and Gruppioni, G. (2017). C4-plant foraging in northern Italy: stable isotopes, Sr/ca and Ba/ca data of human osteological samples from Roccapelago (16th–18th centuries AD). Archaeometry 59, 1119–1134. doi: 10.1111/arcm.12295

CrossRef Full Text | Google Scholar

Lugli, F., Cipriani, A., Tavaglione, V., Traversari, M., and Benazzi, S. (2018). Transhumance pastoralism of Roccapelago (Modena, Italy) early-modern individuals: inferences from Sr isotopes of hair strands. Am. J. Phys. Anthropol. 167, 470–483. doi: 10.1002/ajpa.23643

PubMed Abstract | CrossRef Full Text | Google Scholar

Malmström, H., Storå, J., Dalén, L., Holmlund, G., and Götherström, A. (2005). Extensive human DNA contamination in extracts from ancient dog bones and teeth. Mol. Biol. Evol. 22, 2040–2047. doi: 10.1093/molbev/msi195

PubMed Abstract | CrossRef Full Text | Google Scholar

Meindl, R. S., and Lovejoy, C. O. (1985). Ectocranial suture closure: a revised method for the determination of skeletal age at death based on the lateral-anterior sutures. Am. J. Phys. Anthropol. 68, 57–66. doi: 10.1002/ajpa.1330680106

PubMed Abstract | CrossRef Full Text | Google Scholar

Olofsson, J. K., Mogensen, H. S., Buchard, A., Børsting, C., and Morling, N. (2015). Forensic and population genetic analyses of Danes, Greenlanders and Somalis typed with the Yfiler® plus PCR amplification kit. Forensic Sci. Int. Genet. 16, 232–236. doi: 10.1016/j.fsigen.2015.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Orfanou, E., Himmel, M., Aron, F., and Haak, W. (2021). Minimally-invasive sampling of pars petrosa (os temporale) for ancient DNA extraction V.2. protocols.io doi: 10.17504/protocols.io.bqd8ms9w

CrossRef Full Text | Google Scholar

Parker, C., Rohrlach, A. B., Friederich, S., Nagel, S., Meyer, M., Krause, J., et al. (2020). A systematic investigation of human DNA preservation in medieval skeletons. Sci. Rep. 10:18225. doi: 10.1038/s41598-020-75163-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Paudice, M., Traversari, M., Biatta, C. M., Peñuela, L., Buffelli, F., Spina, B., et al. (2021). Histopathological evidence of coal-workers’ anthracofibrosis in natural mummies of 16th–18th centuries. Med. Secoli 33, 1–10.

Google Scholar

Pilli, E., Vai, S., Caruso, M. G., D’Errico, G., Berti, A., and Caramelli, D. (2018). Neither femur nor tooth: petrous bone for identifying archaeological bone samples via forensic approach. Forensic Sci. Int. 283, 144–149. doi: 10.1016/j.forsciint.2017.12.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Pinhasi, R., Fernandes, D., Sirak, K., Novak, M., Connell, S., Alpaslan-Roodenberg, S., et al. (2015). Optimal ancient DNA yields from the inner ear part of the human petrous bone. PLoS One 10:e0129102. doi: 10.1371/journal.pone.0129102

PubMed Abstract | CrossRef Full Text | Google Scholar

Rapone, C., D’Atanasio, E., Agostino, A., Mariano, M., Papaluca, M. T., Cruciani, F., et al. (2016). Forensic genetic value of a 27 Y-STR loci multiplex (Yfiler® plus kit) in an Italian population sample. Forensic Sci. Int. Genet. 21, e1–e5. doi: 10.1016/j.fsigen.2015.11.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Sarno, S., Boscolo Agostini, R., De Fanti, S., Ferri, G., Ghirotto, S., Modenini, G., et al. (2021). Y-chromosome variability and genetic history of commons from northern Italy. Am. J. Phys. Anthropol. 175, 665–679. doi: 10.1002/ajpa.24302

PubMed Abstract | CrossRef Full Text | Google Scholar

Sirak, K. A., Fernandes, D. M., Cheronet, O., Harney, E., Mah, M., Mallick, S., et al. (2020). Human auditory ossicles as an alternative optimal source of ancient DNA. Genome Res. 30, 427–436. doi: 10.1101/gr.260141.119

CrossRef Full Text | Google Scholar

Sirak, K. A., Fernandes, D. M., Cheronet, O., Novak, M., Gamarra, B., Balassa, T., et al. (2017). A minimally-invasive method for sampling human petrous bones from the cranial base for ancient DNA analysis. Biotechniques 62, 283–289. doi: 10.2144/000114558

PubMed Abstract | CrossRef Full Text | Google Scholar

Traversari, M. (2020). “Nozze a Roccapelago, tra endogamia, esogamia e isonimia,” in Le mummie di Roccapelago. Un progetto pilota di ricerca interdisciplinare tra archeologia, antropologia, storia e scienze applicate. eds. E. Cilli and M. Traversari (Bologna, IT: IBC), 140–148.

Google Scholar

Traversari, M., Feletti, F., Vazzana, A., Gruppioni, G., and Frelat, M. (2016). Three cases of developmental dysplasia of the hip on partially mummified human remains (Roccapelago, Modena, 18th century): a study of palaeopathological indicators through direct analysis and 3D virtual models. Bull. Mem. Soc. Anthropol. Paris 28, 202–212. doi: 10.1007/s13219-015-0140-7

CrossRef Full Text | Google Scholar

Traversari, M., Figus, C., Petrella, E., Piciucchi, S., Vazzana, A., Cilli, E., et al. (2019). Paleopathological analysis of a probable case of Jarcho-Levin syndrome from the 18th century northern Italy. Med. Hist. 3, 39–45.

Google Scholar

Vernarecci, S., Ottaviani, E., Agostino, A., Mei, E., Calandro, L., and Montagna, P. (2015). Quantifiler® Trio Kit and forensic samples management: a matter of degradation. Forensic Sci. Int. Genet. 16, 77–85. doi: 10.1016/j.fsigen.2014.12.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Votrubova-Dubska, J., Vanek, D., Zikmund, J., Mestek, O., Urbanova, V., Brzobohata, H., et al. (2016). Technical note: efficient removal of a PCR inhibitory agent (vivianite) found on excavated bones. Forensic Sci. Int. 261, 8–13. doi: 10.1016/j.forsciint.2015.12.043

PubMed Abstract | CrossRef Full Text | Google Scholar

Walsh, B. (2001). Estimating the time to the most recent common ancestor for the Y chromosome or mitochondrial DNA for a pair of individuals. Genetics 158, 897–912. doi: 10.1093/genetics/158.2.897

PubMed Abstract | CrossRef Full Text | Google Scholar