Application of SNPs with low minor allele frequencies in missing person identification (MPI) through kinship analysis of DNA mixtures.

The need to identify a missing person (MP) through kinship analysis of DNA samples found at a crime scene has become increasingly prevalent. DNA samples from MPs can be severely degraded, contain little DNA and mixed with other contributors, which often makes it difficult to apply conventional methods in practice. This study developed a massively parallel sequencing-based panel that contains 1661 single-nucleotide polymorphisms (SNPs) with low minor allele frequencies (MAFs) (averaged at 0.0613) in the Chinese Han population, and the strategy for relationship inference from DNA mixtures comprising different numbers of contributors (NOCs) and of varying allele dropout probabilities. Based on the simulated dataset and genotyping results of 42 artificial DNA mixtures (NOC = 2-4), it was observed that the present SNP panel was sufficient for balanced mixtures when referenced to the closest relatives (parents/offspring and full siblings). When the mixture profiles suffered from dropout, incorrect assignments were markedly associated with relatedness, NOC and the dropout level. We, therefore, indicate that SNPs with low MAFs could be reliably interpreted for MP identification through the kinship analysis of complex DNA mixtures. Further studies should be extended to more possible scenarios to test the feasibility of this present approach.

Comparing preservation substrates under field conditions for efficient DNA recovery in bone.

Often bones are the only biological material left for the identification of human remains. As situations may occur where samples need to be stored for an extended period without access to cooling, appropriate storage of the bone samples is necessary for maintaining the integrity of DNA for profiling. To simulate DNA preservation under field conditions, pig rib bones were used to evaluate the effects of bone cleaning, buffer composition, storage temperature, and time on DNA recovery from bone samples. Bones were stored in three different buffers: TENT, solid sodium chloride, and ethanol-EDTA, at 20 °C and 35 °C for 10, 20, and 30 days. Bones were subsequently dried and ground to powder. DNA was extracted and quantified. Results show that temperature and storage time have no significant influence on DNA yield. DNA recovery from bones stored in solid sodium chloride or ethanol-EDTA was significantly higher compared to bones stored in TENT, and grinding of bones was facilitated by the extent of dehydration in solid sodium chloride and ethanol-EDTA compared to TENT. Overall, solid sodium chloride was found to be superior over ethanol-EDTA; when it comes to transportation, dry material such as salt eliminates the risk of leaking; it is non-toxic and in contrast to ethanol not classified as dangerous goods. Based on this study's results, we recommend NaCl as a storage substrate for forensic samples in cases where no cooling/freezing conditions are available.

Easykin: a flexible and user-friendly online tool for forensic kinship testing and missing person identification.

In forensic kinship testing and missing person identification, it is a fundamental question to choose the most informative reference relatives, select appropriate genotyping systems, and evaluate the weight of evidence comprehensively. Despite that several useful tools have been developed, they have not addressed these questions satisfactorily. In this paper, we develop a flexible and user-friendly online tool, Easykin, to address the aforementioned issues. It has some promising features: (i) Pedigrees can be constructed easily and presented intuitively with just a few mouse clicks. (ii) System power can be estimated before testing based on certain set of markers and reference relatives. (iii) The pruning function of EasyKin enables users to choose appropriate subsets of available references. (iv) Parameters at a specific LR for a single case may ease evidence interpretation. (v) The user interface (UI) is an HTML-based dashboard, which is friendly to both professional and non-professional users and can be used anytime and anywhere. Here, we presented three common cases as examples to demonstrate how kinship testing and missing person identification can be improved with EasyKin. In conclusion, this tool provides a one-stop solution for forensic use, that is, instructing users to choose appropriate kits and reference relatives before testing, calculating LR in the testing, and providing parameters for data interpretation after testing. EasyKin is freely available at https://forensicsysu.shinyapps.io/EasyKin/ .

Comparison of DNA preservation between adult and non-adult ancient skeletons.

Studies evaluating DNA preservation in non-adults, or comparing preservation in adults and non-adults, are very rare. This study compares the preservation of DNA in the skeletal remains of adults and non-adults. It compares the quality and quantity of DNA recovered from different skeletal elements of adults and non-adults, and from non-adults of different age classes. In addition, the preservation of DNA in males and females is compared. Bone DNA preservation was estimated by measuring nuclear DNA concentration and its degradation, and through STR typing success. The study analyzed 29 adult skeletons and 23 non-adult skeletons from the Ljubljana-Polje archeological site, dating from the seventeenth to nineteenth century, and up to four skeletal elements (petrous bone, femur, calcaneus, and talus) were included. After full demineralization extraction, the PowerQuant System and the PowerPlex ESI 17 Fast System (Promega) were used for qPCR and STR typing, respectively. The results showed that, among the four bone types analyzed, only the petrous bone proved to be a suitable source of DNA for STR typing of non-adult skeletal remains, and DNA yield is even higher than in the adult petrous bone, which can be attributed to the higher DNA degradation observed in the adult petrous bone. In adult skeletons, petrous bones and tali produced high STR amplification success and low DNA yield was observed in adult femurs. The results of this study are applicable for the sampling strategy in routine forensic genetics cases for solving identification cases, including badly preserved non-adult and also adult skeletons.

Comparison of DNA preservation between ribs and vertebrae.

The choice of skeletal element types and their intra-bone parts is important because of differences in DNA preservation, and this must be considered when sampling bones for DNA testing. When incomplete skeletons are found, ribs and vertebrae have been shown to be the most suitable for genetic identification of bones from the torso. This study compares the preservation of DNA between 12th thoracic vertebrae and first ribs to determine which bone type is more suitable for genetic typing. The study analyzed 35 12th thoracic vertebrae and 29 first ribs from one mass grave from the Second World War with commingled skeletal remains excavated. Bone DNA preservation was estimated by measuring nuclear DNA concentration and its degradation and through short tandem repeat (STR) typing success. Previous studies performed on aged skeletal remains have shown that the DNA content of the first ribs and 12th thoracic vertebrae has high intra-bone variability, and this was considered when sampling the bones. After full demineralization extraction, the PowerQuant System (Promega) was used to measure the quantity and quality of DNA, and the GlobalFiler kit (Applied Biosystems) was used for STR typing. The results showed that DNA yield and degradation and STR typing success exhibited no statistically significant difference between first ribs and 12th thoracic vertebrae, and there was no intra-individual difference when comparing only paired bones from the same individuals. Consequently, with intra-bone DNA variability considered, the first ribs or the 12th thoracic vertebrae can be selected when sampling to genetically identify the skeletal remains of highly degraded torsos. HIGHLIGHTS: The first ribs and thoracic vertebrae are the most suitable bones for sampling from the torso. The proximal part of first rib and posterior vertebral column of the 12th thoracic vertebrae yielded the most DNA. The first ribs were compared with the 12th thoracic vertebrae, and the sampling process considered intra-bone DNA variability. The quality and quantity of nuclear DNA and success of STR typing were measured. The first ribs yielded the same DNA yields as well as STR typing success as the 12th thoracic vertebrae. When only the torso is present, it is not of high importance whether the first ribs or the 12th thoracic vertebrae are collected.

Intra-bone nuclear DNA variability and STR typing success in Second World War first ribs.

DNA sampling and typing are used for identifying missing persons or war victims. In recent forensic studies, little focus has been placed on determining intra-bone variability within a single skeletal element. When dealing with aged human bones, complete skeletal remains are rarely present. In cases in which only the torso is available, studies have shown that ribs are one of the most appropriate samples, but intra-bone variability has not yet been studied. A higher degree of remodeling was found to contribute to higher DNA yield in the parts of the skeletal element where the most strain is concentrated. This study explores intra-bone variability in proximal, middle, and distal parts of the first human rib by determining the quantity and quality of DNA using the PowerQuant System (Promega) and autosomal STR typing success using the PowerPlex ESI 17 Fast System (Promega). Thirty first ribs from a single Second World War mass grave were sampled. No variation in DNA degradation was observed across the individual rib. The highest quantity of DNA was measured in the proximal part of the first rib, and in all ribs except three, full or almost full genetic profiles were obtained. Thus, when only the torso is present in archaeological or medico-legal cases, first ribs are recommended to be collected if possible, and the proximal or vertebral ends should be sampled for genetic analysis.

The Missing Person problem through the lens of information theory.

Missing person cases typically require a genetic kinship test to determine the relationship between an unidentified individual and the relatives of the missing person. When not enough genetic evidence has been collected the lack of statistical power of these tests might lead to unreliable results. This is particularly true when just a few distant relatives are available for genotyping. In this contribution, we considered a Bayesian network approach for kinship testing and proposed several information theoretic metrics in order to quantitatively evaluate the information content of pedigrees. We show how these statistics are related to the widely used likelihood ratio values and could be employed to efficiently prioritize family members in order to optimize the statistical power in missing person problems. Our methodology seamlessly integrates with Bayesian modeling approaches, like the GENis platform that we have recently developed for high-throughput missing person identification tasks. Furthermore, our approach can also be easily incorporated into Elston-Stewart forensic frameworks. To facilitate the application of our methodology, we have developed the forensIT package, freely available on CRAN repository, which implements all the methodologies described in our manuscript.

A fast and highly efficient automated DNA extraction method from small quantities of bone powder from aged bone samples.

An efficient extraction method is important for obtaining high-quality DNA from degraded aged bone samples. An automated full-demineralization method using the EDTA and DNA Investigator Kit (Qiagen) combined with Qiagen's biorobots was optimized in our laboratory in the past to extract the DNA from 500 mg of aged bone samples. The purpose of this research was to further improve the method with the aim of reducing the required sample material, shortening the extraction time, and achieving higher throughput. To process extremely small samples, the amount of bone powder was reduced to 75 mg, EDTA was replaced with reagents from the Bone DNA Extraction Kit (Promega), and decalcification was shortened from overnight to 2.5 h. Instead of 50 ml tubes, 2 ml tubes were used, which allows higher throughput. The DNA Investigator Kit (Qiagen) and EZ1 Advanced XL biorobot (Qiagen) was used for DNA purification. A comparison between both extraction methods was made on 29 Second World War bones and 22 archaeological bone samples. The differences between both methods were explored by measuring nuclear DNA yield and STR typing success. After cleaning the samples, 500 mg of bone powder was processed using EDTA, and 75 mg of powder from the same bone was processed using the Bone DNA Extraction Kit (Promega). DNA content and DNA degradation were determined using PowerQuant (Promega), and the PowerPlex ESI 17 Fast System (Promega) was used for STR typing. The results showed that the full-demineralization protocol using 500 mg of bone was efficient for Second World War and archaeological samples, and the partial-demineralization protocol using 75 mg of bone powder was only efficient for the Second World War bones. The improved extraction method-for which significantly lower amounts of bone powder can be used, the extraction process is faster, and higher throughput of bone samples is possible-is applicable for genetic identification of relatively well-preserved aged bone samples in routine forensic analyses.

Building an isoscape based on tooth enamel for human provenance estimation in Brazil.

In this study, we present a correlation between δ18OC values of carbonate in tooth enamel samples from the modern Brazilian population and the available δ18ODW data for the meteoric water from the Global Network of Isotopes in Precipitation (GNIP). Tooth enamel from 119 Brazilian individuals from five different regions of the country were analyzed. The δ18OC isoscape obtained is in good agreement with the isoscape based on regional meteoric and drinking water. The regression matrix obtained for the δ18O values of the carbonate tooth enamel and meteoric water was used to build an isoscape using the regression-kriging approach. Our data show that Brazil can be divided in two main regions with respect to the δ18O values of the carbonate tooth enamel: (1) the most easterly part of the northeast region, which is characterized by a warm and dry climate and (2) the remainder of the country, stretching from the Amazon rain forest to the more southernly regions. The data herein reported can be used for forensic purposes related to human identification.

Identification of missing persons through kinship analysis by microhaplotype sequencing of single-source DNA and two-person DNA mixtures.

In forensic applications, there is an increasing demand for the analysis of DNA profiles arising from missing person identification (MPI) cases. A specific DNA profile may originate from a single source or more than one contributor (i.e., a DNA mixture). When direct references are not available, indirect relative references can be used to identify missing persons by kinship analysis. As a novel kind of multiallelic marker, microhaplotypes have proven promising for relatedness determination and mixture deconvolution. Herein, we developed a large panel of 185 microhaplotype markers and demonstrated its application in different scenarios of relationship inference through a simulation study and real pedigree analysis, combined with probabilistic genotyping models for data interpretation. Based on single-source profiles, it was shown that the present microhaplotype panel was sufficient for pairwise close relative testing (parent/child, full-sibling and 2nd-degree relative). For more distant relatives (3rd-degree relatives), there was a clear improvement when data from one well-chosen extra relative were available. We further sought to evaluate the theoretical systematic effectiveness and actual performance of microhaplotype markers in identifying the contribution of a missing pedigree member to a two-person mixture (as a minor donor). It was observed that 100% correct assignments were made in the balanced mixtures (with no dropout) when referenced to close relatives. When the mixture profiles suffered from dropout, incorrect assignments of minor donors were markedly associated with relatedness and the dropout level. Meanwhile, the studied scenarios generally exhibited zero or very low false-positive rates, indicating a low probability of incorrectly assigning an unrelated contributor as a close relative of the reference. Our results indicate that microhaplotype data can be reliably interpreted for identifying missing persons through kinship analysis based on DNA profiles of single-source samples or two-person mixtures. Furthermore, this study could be extended to more complex scenarios, such as determining the relatedness of contributors in (or among) mixed DNA profiles, if combined with different statistical frameworks.

Making decisions in missing person identification cases with low statistical power.

The present work proposes a general strategy for dealing with missing person identification cases through DNA-database search. Our main example is the identification of abducted children in the last civic-dictatorship of Argentina, known as the "Missing Grandchildren of Argentina". Particularly we focus on those pedigrees where few, or only distant relatives of the missing person are available, resulting in low statistical power. For such complex cases we provide a statistical method for selecting a likelihood ratio (LR) threshold for each pedigree based on error rates. Furthermore, we provide an open-source user friendly software for computing LR thresholds and error rates. The strategy described in the paper could be applied to other large-scale cases of DNA-based identification hampered by low statistical power.

Comparison of nuclear DNA yield and STR typing success in Second World War petrous bones and metacarpals III.

DNA yield and STR typing success differ among skeletal element types and within individual bones. Consequently, it is necessary to identify the skeletal elements and their intra-skeletal parts that will most likely yield utilizable and informative endogenous DNA for human identification of skeletal remains. The petrous portion of the temporal bone has been shown to be the most suitable skeletal part for sampling archaeological skeletons, and it has also been used successfully in some forensic cases. When all representative bone types were analyzed for three complete Second World War skeletons, metatarsals and metacarpals yielded more DNA than petrous bones (which generated full profiles even if the DNA yield was lower) and, among almost 200 Second World War metatarsals and metacarpals analyzed, metacarpals III were found to be the highest-yielding bones. To further improve the sampling strategy in DNA analysis of aged skeletal remains, a comparison between 26 petrous bones and 30 metacarpals III from Second World War skeletons was made considering intra-bone DNA yield variability. In metacarpals III only epiphyses were sampled, and in petrous bones only the dense part within the otic capsule was used. To exclude the influence of taphonomic issues as much as possible, petrous bones and metacarpals III from a single Second World War mass grave were examined. The difference between petrous bones and metacarpals III was explored by measuring nuclear DNA yield and success of STR typing. After cleaning the samples, full demineralization extraction was used to decalcify 0.5 g of powdered bone. PowerQuant (Promega) was used to determine DNA content and DNA degradation rates, and STR typing was performed using the PowerPlex ESI 17 Fast System (Promega). Metacarpals III produced the same DNA yields and STR typing success as petrous bones with no intra-individual difference observed in concentration of DNA, degradation rate, percentage of successfully amplified alleles, and intensity of electrophoretic signals. Sampling and processing of metacarpal III epiphyses is consequently recommended for genetic identification of highly degraded skeletal remains in routine forensic casework and in buried non-commingled aged skeletal remains whenever metacarpals III are preserved. Metacarpals III are easy to sample and less prone to contamination with modern DNA because no saw is needed for sampling in comparison to the petrous portion of the temporal bone. The data obtained in this study further improve the sampling strategy for genetic identification of Second World War skeletal remains in Slovenia.

Identification of a Slovenian prewar elite couple killed in the Second World War.

Genetic identification of a Slovenian prewar elite couple killed in 1944 was performed by typing autosomal and Y-chromosomal STRs, and phenotypic HIrisPlex SNPs for hair and eye color prediction were analyzed for the female skeleton using next-generation sequencing (NGS) technology. The clandestine grave containing the couple's skeletal remains was found in 2015 and only the partial remains were found. Living distant relatives could be found only for the male victim. Because of a lack of comparative reference samples, it was not possible to identify the female victim through autosomal and mitochondrial DNA typing. However, the possibility of comparison of eye and hair color with a painting exhibited in the City Museum of Ljubljana by the prominent Slovenian painter Ivana Kobilca existed. Nuclear DNA obtained from the samples was quantified using the PowerQuant System, and then STR typing was carried out with different autosomal and Y-STR kits. From 0.09-9.36 ng DNA/g of powder was obtained from teeth and bones analyzed. Complete autosomal and Y-STR profiles made it possible to identify the male skeleton via comparison with two nephews. For the female victim, predicted eye and hair color was compared to colors on the painting. Kobilca's painting confirms the genetically predicted eye and hair color. After more than seventy years, the skeletal remains of the couple were handed over to their relatives, who buried the victims with dignity in a family grave.

Intra-bone nuclear DNA variability and STR typing success in Second World War 12th thoracic vertebrae.

Bones are an important source of DNA for identification in forensic medicine, especially when the remains are skeletonized, which is the case when dealing with victims of the Second World War. Often the amount of bone available for sampling is limited, and therefore it is crucial to sample the bone segment with the highest adequate DNA quantity for identification. Studies performed on all representative skeletal element types of the human body showed that the amount of DNA obtained from different skeletal elements of different body regions varies greatly. When bones from torso were analyzed, thoracic vertebrae outperformed other vertebrae (cervical and lumbar) and, alongside the first ribs, were among the most appropriate bone elements for identification purposes. It was also shown that the quantity of DNA varies significantly within a single bone type. This study focused on exploring intra-bone DNA variability between five parts of 12th thoracic vertebrae (laminae + spinous process, pedicles + transverse processes, and corpus right, left, and middle). The research was based on the theory that the distribution of body weight and consequently bone remodeling, as well as the ratio between cancellous and cortical bone, contribute to different quantities of DNA in different parts of vertebra sampled. The vertebrae were cleaned and cut into five parts, and each part was completely ground to obtain homogenous bone powder. Half a gram of powder from each part was decalcified using a full demineralization extraction method. The DNA was purified in a Biorobot EZ1 machine (Qiagen). DNA quantity and quality were determined using the PowerQuant System (Promega) and autosomal STR typing success using the GlobalFiler Amplification Kit (Applied Biosystems). Thirty-five 12th thoracic vertebrae were sampled from a single Second World War mass grave. The best results with the highest DNA quantity were found in laminae and the spinous process, and among them all vertebrae analyzed yielded full STR profiles except three, where only a few dropouts occurred. The second-ranked bone part was the pedicles and transverse processes. The comparison of DNA degradation in the vertebral segments analyzed does not show statistically significant differences. Considering our research, when only the torso is available for identification, the 12th thoracic vertebra should be collected and the vertebral arch should be sampled for genetic analyses.

High DNA yield from metatarsal and metacarpal bones from Slovenian Second World War skeletal remains.

DNA yield varies by anatomical region, and the selection of bone types that yield maximum recovery of DNA is important to maximize the success of human identification of skeletal remains. The goal of our study was to explore inter- and intra-individual variation in DNA content by measuring nuclear DNA quantity and quality and autosomal STR typing success to determine the most promising skeletal elements for bone sampling. To exclude the influence of taphonomic issues as much as possible, three complete male skeletons from a single Second World War mass grave were examined and all representative skeletal element types of the human body were analyzed. Forty-eight different types of bones from the head, torso, arm, leg, hand, and foot were sampled from each skeleton, 144 bones altogether. The samples were cleaned, and half a gram of bone powder was decalcified using a full demineralization extraction method. The DNA was purified in a Biorobot EZ1 (Qiagen). DNA content and rates of DNA degradation were determined with the PowerQuant (Promega), and the Investigator ESSplex SE QS (Qiagen) was used for STR typing. The highest-yielding bones mostly produced the most complete STR profiles. Among the skeletal elements containing on average the most DNA and producing the most complete profiles in all three skeletons examined were metacarpals, metatarsals, and the petrous portion of the temporal bone. Metatarsals and metacarpals can easily be sampled without using a saw, thus reducing potential DNA contamination. Skeletons from the Second World War can be used as a model for poorly preserved skeletal remains, and the results of the investigation can be applied for genetic identification of highly degraded skeletal remains in routine forensic casework. Although the research was limited to only three skeletons found in a unique mass grave, the data obtained could contribute to sampling strategies for identifying old skeletal remains. More Second World War skeletons will be analyzed in the future to investigate inter-bone variation in the preservation of DNA.

Strategy for STR typing of bones from the Second World War combining CE and NGS technology: A pilot study.

The genetic identification of skeletal remains found in Second World War mass graves is complicated because of the poor quality of the samples. The aim of this study was to set up a workflow for STR typing of such samples combining PCR/CE and PCR/NGS technologies. To this end, 57 DNA samples from an equal number of 75-year-old femurs were studied. After a first round of PCR typing using GlobalFiler CE, 42 samples yielded a full profile and were therefore submitted to our standard workflow. The 15 samples that yielded no or a limited number (2-17/21) of autosomal STR markers as well four bone control samples that provided a full profile with the conventional PCR/CE test were typed in duplicate by the GlobalFiler NGS kit. Despite the degradation of the samples, which resulted in lower coverage and a lower % of on-target reads, reliable sequencing data were obtained from 16/19 samples. The use of a threshold of 30× for the locus call led to a consensus profile (cp) of 20-31/31 STR autosomal loci in 10 samples and to a cp of 8-10/31 loci in two samples, whereas the four control samples yielded a cp of 26-31/31 loci. Finally, the data of the NGS typing were combined with those of the CE typing. This last task allowed us to recover (on average) three alleles per sample and to increase the number of the heterozygous patterns in 37 cases. In total, the combined approach proposed here made possible the genetic typing of 65-100% of the autosomal STR markers in 10/15 (66.6 %) skeletal remains that yielded no or very poor results with the conventional PCR/CE approach. However, because several artefacts (such as allelic drop-out and allelic drop-in) were scored, the risk of mistyping cannot be neglected.

Whole-genome sequencing of human remains to enable genealogy DNA database searches - A case report.

Recently a number of high profile crime cases (e.g. the "Golden State Killer") have successfully been solved or given new leads with the use of genome wide DNA data in combination with pairwise matching from individuals present in genealogy DNA databases. Such databases will primarily involve distant relatives which in turn require a large amount of genetic information, in the range of several hundred thousand to millions of SNPs, to be genotyped. While it nowadays is fairly straightforward to obtain such as data from high quality and high quantity DNA, it is still a challenge for degraded DNA of low quantity such in the case of forensic samples. Here we present a successful effort in obtaining genome-wide genotype data from human remains. The goal was to get investigative leads in order to identify the remains of an unknown male ("the Ekeby man") that was found murdered in the south of Sweden in 2003. Whole-genome sequencing was performed on DNA originating from a bone sample. Three replicates of libraries were prepared using ThruPLEX DNA-seq Kit (Takara) which were sequenced on a HiSeq X instrument (Illumina). A mean coverage of 30X was obtained when the sequencing reads were mapped to a human reference genome. Following further bioinformatic processing, allele calling, quality checks and filtering to match the genealogy DNA database SNPs, genotypes for approximately one million SNPs were established. The resulting SNP genotypes were then used to search for relatives in the genealogy DNA database GEDmatch (www.gedmatch.com). A candidate list of relatives was obtained which was further processed using traditional genealogy methods in order to get leads about the identity of the unknown. In summary, this report shows how whole-genome sequencing successfully can be applied on forensic samples to create the SNP genotypes required for searches in genealogy DNA databases for the purpose of generating leads to identify missing or unknown persons, including perpetrators and victims.

Identifying victims of the largest Second World War family massacre in Slovenia.

The killings during the Second World War, with nearly one hundred thousand victims, is one of the greatest losses of life in Slovenia's modern history. This article presents the genetic identification of the victims of the largest family massacre that occurred in Slovenia, in which 10 members of the same family were killed. Seven of them were buried in a hidden mass grave and only two children survived. In 2015 and 2016, two graves were found and three incomplete female skeletons and at least three incomplete male skeletons were exhumed. A total of 12 bones and teeth were analysed and compared to two living relatives. Extracted DNA was quantified using the PowerQuant kit, and various autosomal and Y-STR kits were used for STR typing. Up to 2.7 ng DNA/g of powder was acquired from the samples analysed. We managed to obtain nuclear DNA for successful STR typing from seven bones and one molar. From the female grave, autosomal profiles were obtained only from one skeleton, and from the male grave from five out of six femurs. The relationships between the males were additionally confirmed by analyses of Y-STRs. STR profiles made possible the identification of four family members; one of the aunts from the female grave, and two uncles and the father of the surviving children, who were used as family references, from the male grave. The product rule was used to calculate a combined likelihood ratio for autosomal and Y-STRs, and statistical analyses showed high confidence of correct identification with posterior probability (PP) greater than 99.9 % for three out of four victims identified. For identifying the aunt, the PP obtained after ESI-17 and NGM STR typing was too low. To increase the PP, the next-generation sequencing Precision ID GlobalFiler NGS STR Panel was used and, after the analysis of additional STR loci, the statistical analysis showed a PP greater than 99.9 %, indicating that a sufficient number of genetic markers had been investigated in identifying the skeletal remains of the aunt. An elimination database containing the genetic profiles of all individuals that had been in contact with the bones was created to ensure traceability in case of contamination, and no matches were found. After more than 70 years, the skeletal remains were returned to the surviving children, who buried their relatives in a family grave.

Prioritising family members for genotyping in missing person cases: A general approach combining the statistical power of exclusion and inclusion.

Missing person identification typically involves genetic matching of a person of interest against relatives of the missing person. In cases with few available relatives, exhumations or other substantial efforts may be necessary in order to secure adequate statistical power. We propose a simulation approach for solving prioritisation problems arising in such cases. Conditioning on the already typed individuals we estimate the power of each alternative, both to detect the true person, and to exclude false candidates. Graphical summaries of the simulations are given in complementary power plots, facilitating interpretation and decision making. Through a series of examples originating from the well-known Missing grandchildren of Argentina we demonstrate that our method may untangle complex prioritisation problems and other power-related questions. In particular we offer novel insights in recent cases where only children of the potential match are available for testing. We also show that X-chromosomal markers may give high statistical power in missing person identification, but that this requires careful selection of relatives for genotyping. All simulations, power calculations and plots are done with the R package forrel.

Highly efficient automated extraction of DNA from old and contemporary skeletal remains.

We optimised the automated extraction of DNA from old and contemporary skeletal remains using the AutoMate Express system and the PrepFiler BTA kit. 24 Contemporary and 25 old skeletal remains from WWII were analysed. For each skeleton, extraction using only 0.05 g of powder was performed according to the manufacturer's recommendations (no demineralisation - ND method). Since only 32% of full profiles were obtained from aged and 58% from contemporary casework skeletons, the extraction protocol was modified to acquire higher quality DNA and genomic DNA was obtained after full demineralisation (FD method). The nuclear DNA of the samples was quantified using the Investigator Quantiplex kit and STR typing was performed using the NGM kit to evaluate the performance of tested extraction methods. In the aged DNA samples, 64% of full profiles were obtained using the FD method. For the contemporary skeletal remains the performance of the ND method was closer to the FD method compared to the old skeletons, giving 58% of full profiles with the ND method and 71% of full profiles using the FD method. The extraction of DNA from only 0.05 g of bone or tooth powder using the AutoMate Express has proven highly successful in the recovery of DNA from old and contemporary skeletons, especially with the modified FD method. We believe that the results obtained will contribute to the possibilities of using automated devices for extracting DNA from skeletal remains, which would shorten the procedures for obtaining high-quality DNA from skeletons in forensic laboratories.