Familias 3 - Extensions and new functionality

Evaluation of the usefulness of insertion-null markers in critical skeletal remains

Forensic DNA analysis in compromised skeletal remains may pose challenges due to DNA degradation, often resulting in partial or negative autosomal STRs profiles. To address this issue, alternative approaches such as mitochondrial DNA or SNPs typing may be employed; however, they are labour-intensive and costly. Insertion-null alleles (INNULs), short interspersed nuclear elements, have been suggested as a valuable tool for human identification in challenging samples due to their small amplicon size. A commercial kit including 20 INNULs markers along with amelogenin (InnoTyper® 21) has been developed. This study assesses its utility using degraded skeletal remains, comparing the results obtained (the number of detected alleles, RFU values, PHR, and the number of reportable markers) to those obtained using GlobalFiler™. Subsequently, the random match probability of the two profiles for each sample was determined using Familias version 3 to evaluate the power of discrimination of the results obtained from each kit. In every sample, InnoTyper® 21 yielded more alleles, higher RFU values, and a greater number of reportable loci. However, in most cases, both profiles were similarly informative. In conclusion, InnoTyper® 21 serves as a valuable complement to the analysis of challenging samples in cases where a poor or negative profile was obtained.

An overview of autosomal STRs and identity SNPs in a Norwegian population using massively parallel sequencing

In recent years, probabilistic genotyping software has been adapted for the analysis of massively parallel sequencing (MPS) forensic data. Likelihood ratios (LR) are based on allele frequencies selected from populations of interest. This study provides an outline of sequence-based (SB) allele frequencies for autosomal short tandem repeats (aSTRs) and identity single nucleotide polymorphisms (iSNPs) in 371 individuals from Southern Norway. 27 aSTRs and 94 iSNPs were previously analysed with the ForenSeq™ DNA Signature Prep Kit (Verogen). The number of alleles with frequencies less than 0.05 for sequenced-based alleles was 4.6 times higher than for length-based alleles. Consistent with previous studies, it was observed that sequence-based data (both with and without flanks) exhibited higher allele diversity compared to length-based (LB) data; random match probabilities were lower for SB alleles confirming their advantage to discriminate between individuals. Two alleles in markers D22S1045 and Penta D were observed with SNPs in the 3´ flanking region, which have not been reported before. Also, a novel SNP with a minor allele frequency (MAF) of 0.001, was found in marker TH01. The impact of the sample size on minor allele frequency (MAF) values was studied in 88 iSNPs from Southern Norway (n = 371). The findings were then compared to a larger Norwegian population dataset (n = 15,769). The results showed that the smaller Southern Norway dataset provided similar results, and it was a representative sample. Population structure was analyzed for regions within Southern Norway; FST estimates for aSTR and iSNPs did not indicate any genetic structure. Finally, we investigated the genetic differences between Southern Norway and two other populations: Northern Norway and Denmark. Allele frequencies between these populations were compared, and we found no significant frequency differences (p-values > 0.0001). We also calculated the pairwise FST values per marker and comparisons between Southern and Northern Norway showed small differences. In contrast, the comparisons between Southern Norway and Denmark showed higher FST values for some markers, possibly driven by distinct alleles that were present in only one of the populations. In summary, we propose that allele frequencies from each population considered in this study could be used interchangeably to calculate genotype probabilities.

Nails as optimal source of DNA for molecular identification of 5 decomposed bodies recovered from seawater: from Y-ancestry to personal identification

Molecular identification of extremely compromised human remains in forensic field is usually performed from DNA typing of bones, which are a difficult sample to work with. Moreover, autosomal STR profiles do not always result in the identification of the donor due to lack of comparisons or non-hit throughout database searching. An attempt to overcome these issues is represented by fingernails as an alternative DNA source and Y-STRs typing to infer both geographical and familial ancestry of the unknown donor. In this study, we analyzed both 24 autosomal and 27 Y-chromosome STRs from unidentified human remains (UHRs) of five males recovered from the water near the southwestern coast of Sardinia by the Italian Harbor Master's Office. Nail clippings provided an optimal source of autologous DNA for molecular identification in a very short time, producing complete autosomal and Y-STR profiles even under conditions of high body degradation. Unfortunately, no match neither compatibility occurred using autosomal STRs (aSTRs), initially. Upon analyzing the Y-haplotypes, we found out they had already been observed in northern Africa, providing us important investigative leads. This prompted the International Criminal Police Organization (INTERPOL) to provide us with references of alleged relatives that were then confirmed to be related. The use of fingernails represents an excellent DNA source especially for genetic identification of decomposed bodies recovered in seawater environment. Notably, DNA extracted from nails gave high-quality Y-STR haplotypes by which predicting paternal ancestry of the unidentified donors may result fundamental in the forensic investigative context.

Determining the effects of genetic linkage when using a combination of STR and SNP loci for kinship testing

The pedigree likelihood ratio (LR) can be used for determining kinship in the forensic kinship testing. LR can be obtained by analyzing the DNA data of Short tandem repeat (STR) and single nucleotide polymorphism (SNP) loci. With the advancement of biotechnology, increasing number of genetic markers have been identified, thereby expanding the pedigree range of kinship testing. Moreover, some of the loci are physically closer to each other and genetic linkage between loci is inevitable. LRs can be calculated by accounting for linkage or ignoring linkage (LRlinkage and LRignore, respectively). GeneVisa is a software for kinship testing (www.genevisa.net) and adopts the Lander-Green algorithm to deal with genetic linkage. Herein, we used the simulation program of the software GeneVisa to investigate the effects of genetic linkage on 1st-degree, 2nd-degree, and 3rd-degree kinship testing. We used this software to simulate LRlinkage and LRignore values based on 43 STRs and 134 SNPs in commercial kits by using the allele frequency rate and genetic distance data of the European population. The effects of linkage on LR distribution and LRs of routine cases were investigated by comparing the LRlinkage values with the LRignore values. Our results revealed that the linkage effect on LR distributions is small, but the effect on LRs of routine cases may be large. Moreover, the results indicated that the discriminatory power of genetic markers for kinship testing can be improved by accounting for linkage.

Joint application of A-InDels and miniSTRs for forensic personal, full and half sibling identifications, and genetic differentiation analyses in two populations from China

BACKGROUND

Previously, a novel multiplex system of 64 loci was constructed based on capillary electrophoresis platform, including 59 autosomal insertion/deletions (A-InDels), two Y-chromosome InDels, two mini short tandem repeats (miniSTRs), and an Amelogenin gene. The aim of this study is to evaluate the efficiencies of this multiplex system for individual identification, paternity testing and biogeographic ancestry inference in Chinese Hezhou Han (CHH) and Hubei Tujia (CTH) groups, providing valuable insights for forensic anthropology and population genetics research.

RESULTS

The cumulative values of power of discrimination (CDP) and probability of exclusion (CPE) for the 59 A-InDels and two miniSTRs were 0.99999999999999999999999999754, 0.99999905; and 0.99999999999999999999999999998, 0.99999898 in CTH and CHH groups, respectively. When the likelihood ratio thresholds were set to 1 or 10, more than 95% of the full sibling pairs could be identified from unrelated individual pairs, and the false positive rates were less than 1.2% in both CTH and CHH groups. Biogeographic ancestry inference models based on 35 populations were constructed with three algorithms: random forest, adaptive boosting and extreme gradient boosting, and then 10-fold cross-validation analyses were applied to test these three models with the average accuracies of 86.59%, 84.22% and 87.80%, respectively. In addition, we also investigated the genetic relationships between the two studied groups with 33 reference populations using population statistical methods of FST, DA, phylogenetic tree, PCA, STRUCTURE and TreeMix analyses. The present results showed that compared to other continental populations, the CTH and CHH groups had closer genetic affinities to East Asian populations.

CONCLUSIONS

This novel multiplex system has high CDP and CPE in CTH and CHH groups, which can be used as a powerful tool for individual identification and paternity testing. According to various genetic analysis methods, the genetic structures of CTH and CHH groups are relatively similar to the reference East Asian populations.

Witness for resolution: post-conflict quadratic affiliation in semi-free ranging pigs

In social mammals, post-conflict resolution can involve the reunion of former opponents (reconciliation), spontaneous/solicited post-conflict affiliation of a third party with either opponent (triadic contacts), and affiliation between other individuals (hereafter bystanders; quadratic contacts). Quadratic contacts-possibly informing complex cognitive abilities-have been neglected in post-conflict studies. We investigated quadratic affiliation in semi-free ranging pigs Sus scrofa, at the ethical farm Parva-Domus (Cavagnolo, Italy). Kinship was known. We collected behavioral data on adult pigs (n = 104) via video recordings (43 h) followed by video analyses. Affiliative and anxiety behaviors between bystanders were collected under post-conflict (PC; following a conflict between non-bystanders) and matched-control (MC; no conflict) conditions. Quadratic affiliation was present in pigs, as bystanders affiliated more in PC than MC, and such affiliation was followed by a decrease in the anxiety behaviors of both the interacting bystanders. Thus, quadratic contacts may be partly aimed at reducing one's own anxiety (intrinsic regulation). Quadratic affiliation was highest between closely related bystanders, which suggests that such affiliation may be most effective when close kin is involved. Quadratic affiliation was lowest after reconciliation and spontaneous triadic contacts. This suggests that direct peacemaking between opponents and spontaneous triadic contacts with close kin may most likely replace quadratic affiliation. Hence, pigs can be influenced by the negative events that affect other pigs-but not themselves-and their response may be modulated by social factors. Such non-random quadratic affiliation may point toward the presence of elements of social appraisal abilities in pigs.

Using two self-developed InDel panels to explore forensic traits and ancestral components in the Hui group

To address the challenges faced by forensic examiners in degraded DNA analysis, we have developed two different panels for various forensic applications, encompassing an AIM-InDel panel for ancestry inference and a Multi-InDel panel for individual identification, respectively. Herein, the efficiencies of these two panels were tested in the Chinese Hui group. By calculating forensic parameters and simulating family relationships, we verified that the Multi-InDel panel could be an effective tool for individual identification, paternity testing, and could assist in kinship identification in the Hui group. For full siblings, the true positive rate of kinship discrimination was 96.553%, when the threshold of log10LR was 1. The cumulative probability of matching as well as cumulative probability of exclusion were 3.8117 × 10-26 and 0.999999722, respectively. Meanwhile, we found that the AIM-InDel panel was effective for bio-geographic ancestry inference, with the vast majority of loci contributing significantly to distinguish East Asian, African, and European populations. By studying the population structure of the Hui ethnic minority, the genetic distance to the Beijing Han population was the closest among the 26 reference populations, which was similar to previous reports. In summary, we have developed two panels which can be effectively applied to the Hui group for individual identification, parentage testing and bio-geographic ancestry inference.

Targeted sequencing of high-density SNPs provides an enhanced tool for forensic applications and genetic landscape exploration in Chinese Korean ethnic group

BACKGROUND

In this study, we present a NGS-based panel designed for sequencing 1993 SNP loci for forensic DNA investigation. This panel addresses unique challenges encountered in forensic practice and allows for a comprehensive population genetic study of the Chinese Korean ethnic group. To achieve this, we combine our results with datasets from the 1000 Genomes Project and the Human Genome Diversity Panel.

RESULTS

We demonstrate that this panel is a reliable tool for individual identification and parentage testing, even when dealing with degraded DNA samples featuring exceedingly low SNP detection rates. The performance of this panel for complex kinship determinations, such as half-sibling and grandparent-grandchild scenarios, is also validated by various kinship simulations. Population genetic studies indicate that this panel can uncover population substructures on both global and regional scales. Notably, the Han population can be distinguished from the ethnic minorities in the northern and southern regions of East Asia, suggesting its potential for regional ancestry inference. Furthermore, we highlight that the Chinese Korean ethnic group, along with various Han populations from different regional areas and certain northern ethnic minorities (Daur, Tujia, Japanese, Mongolian, Xibo), exhibit a higher degree of genetic affinities when examined from a genomic perspective.

CONCLUSION

This study provides convincing evidence that the NGS-based panel can serve as a reliable tool for various forensic applications. Moreover, it has helped to enhance our knowledge about the genetic landscape of the Chinese Korean ethnic group.

Review: Computational analysis of human skeletal remains in ancient DNA and forensic genetics

Degraded DNA is used to answer questions in the fields of ancient DNA (aDNA) and forensic genetics. While aDNA studies typically center around human evolution and past history, and forensic genetics is often more concerned with identifying a specific individual, scientists in both fields face similar challenges. The overlap in source material has prompted periodic discussions and studies on the advantages of collaboration between fields toward mutually beneficial methodological advancements. However, most have been centered around wet laboratory methods (sampling, DNA extraction, library preparation, etc.). In this review, we focus on the computational side of the analytical workflow. We discuss limitations and considerations to consider when working with degraded DNA. We hope this review provides a framework to researchers new to computational workflows for how to think about analyzing highly degraded DNA and prompts an increase of collaboration between the forensic genetics and aDNA fields.

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/ .

Adapting an established Ampliseq microhaplotype panel to nanopore sequencing through direct PCR

We have adapted an established Ampliseq microhaplotype panel for nanopore sequencing with the Oxford Nanopore Technologies (ONT) system, as a cost-effective and highly scalable solution for forensic genetics applications. For this purpose, we designed a protocol combining direct PCR amplification from unextracted DNA with ONT library construction and sequencing using the MinION device and workflow. The analysis of reference samples at input amounts of 5-10 ng of DNA demonstrates stable coverage patterns, allele balance, and strand bias, reaching profile completeness and concordance rates of ∼95%. Similar levels were achieved when using direct-PCR from blood, buccal and semen swabs. Dilution series results indicate sensitivity is maintained down to 250 pg of input DNA, and informative profiles are produced down to 62.5 pg. Finally, we demonstrated the forensic utility of the nanopore workflow by analyzing two third degree pedigrees that showed low likelihood ratio values after the analysis of an extended panel of 38 STRs, achieving likelihood ratios 2-3 orders of magnitude higher when testing with the MinION-based haplotype data.

Privacy and ethical challenges of the Amelogenin sex test in forensic paternity/kinship analysis: Insights from a 13-year case history

The Amelogenin sex test included in forensic DNA typing kits has the potential to identify congenital conditions such as differences/disorders of sex development (DSD). It can also reveal mismatches between genotypic sex and gender marker in identity documents of transgender persons who obtained legal gender recognition. In a 13-year case history of paternity/kinship tests, involving n = 962 females and n = 1001 males, two mismatches between Amelogenin sex test (male) and gender marker (female), and three cases of chromosomal DSD (Klinefelter syndrome) were observed. The concrete risk of observing Amelogenin anomalies, their potential causes, and the context in which they occur (forensic, i.e. non-medical) mean that laboratory operators are called to strike a complex balance between privacy interests and individual health rights when providing preliminary information and reporting Amelogenin incidental findings. This case history argues for the need of a more responsible approach towards the Amelogenin sex test in the forensic community.

Applications of 1993 single nucleotide polymorphism loci in forensic pairwise kinship identifications and inferences

Kinship testing plays critical roles in criminal investigations, missing person searches, civil disputes, as well as identifying disaster victims. The existing commonly used short tandem repeat (STR) loci have limited effectiveness in the identification of second-degree and more distant kinships. In this study, a total of 1993 SNP loci of 119 Chinese Han individuals from eight families were sequenced on the MGISEQ-2000RS platform. The system powers of this panel for kinship identifications were evaluated based on both the likelihood ratio (LR) and identical by state (IBS) methods. The results indicated that this panel could be used as an effective tool to kinship analyses including paternity testing, full sibling testing, second-degree kinships, and first cousin kinship analyses. Both the LR and IBS methods could be applied in distinguishing first-degree and second-degree pairs from unrelated individuals. Based on the 1993 SNP loci, LR>1000 and LR<0.001 are recommended as the thresholds of identifying first-cousin kinships from unrelated individuals, and the system power of such thresholds was 0.9470. Besides, kinship coefficients for different kinship pairs were estimated and then were used to predict the kinships for pairwise individuals. This panel performs an effective kinship inference power for the predictions of first-degree, second-degree kinships and unrelated individual pairs, while presenting low sensitivity in the prediction of first-cousin kinships.

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

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

Incorporating non-genetic evidence in large scale missing person searches: A general approach beyond filtering

The search for missing persons implies several steps, from the preliminary investigation that involves collecting background data related to the case to the genetic kinship testing. Despite its crucial importance in identifications, only some approaches mathematically formalize the possibility of using preliminary investigation data. In some cases, a filtering strategy is applied, which implies selecting a subset of possible victims where some non-genetic variables perfectly match those of the missing. Through a Bayesian approach, we propose a mathematical model for computing the prior odds based on non-genetic variables usually collected during the preliminary investigation, such as biological sex, hair colour, and age. We use computational simulations to show how to incorporate these prior odds in DNA-database searches. Importantly, our results suggest that applying the proposed model leads to better search performance in underpowered cases from the genetic point of view, where few or distant relatives of the missing person are available for genotyping. Furthermore, the results are also helpful when using non-genetic data for prior odds in well-powered cases, where genetic data are enough to reach a reliable conclusion. It performs better than other approaches, such as using non-genetic data for filtering. The software mispitools, freely available on CRAN, implements all described methods (https://CRAN.R-project.org/package=mispitools).

Paternity pseudo-exclusion caused by tetragametic chimerism in a gestational surrogacy case

Chimerism is a condition when an organism is composed of two or more populations of genetically distinct cells. Chimerism often produces curious results of medical and genetic investigations and could be a major cause of false negative conclusions in parentage testing. Here we describe a paternity pseudo-exclusion due to tetragametic chimerism in a gestational surrogacy case originated in a fertility clinic. Initial analysis using a buccal swab from the child and a peripheral blood sample from the father showed paternity exclusion at 6 STR loci. To find out the reason for observed paternal discrepancy father's semen sample used for IVF and samples from other tissues were genotyped. Buccal swabs, semen, hair follicles, nail clippings and cerumen showed identical mixed autosomal STR profiles originated from two genetically different cell lines and for all the 24 informative loci contain paternal obligate alleles. Results of Y-STR profiling of all paternal sample types showed a DNA profile originated from a single man. The mixed profiles obtained for different tissue types suggest that two genetically different cell lines contributed to formation of both the endoderm and the ectoderm of the father. The mesoderm seems to be monoclonal having originated from a genetically homogenous cell line as evidenced by the STR profile of peripheral blood. Such allelic pattern for various tissues suggests that the clonal origin happened at the very early stage of embryonic development. Approaches to minimise the rate of false exclusions in DNA parentage tests due to chimerism are discussed.

Evaluation of microhaplotype panels for complex kinship analysis using massively parallel sequencing

In recent years, microhaplotypes (MHs) have become a research hotspot within the field of forensic genetics. Traditional MHs contain only SNPs that are closely linked within short fragments. Herein, we broaden the concept of general MHs to include short InDels. Complex kinship identification plays an important role in disaster victim identification and criminal investigations. For distant relatives (e.g., 3rd-degree), many genetic markers are required to enhance power of kinship testing. We performed genome-wide screening for new MH markers composed of two or more variants (InDel or SNP) within 220 bp based on the Chinese Southern Han from the 1000 Genomes Project. An NGS-based 67plex MH panel (Panel B) was successfully developed, and 124 unrelated individual samples were sequenced to obtain population genetic data, including alleles and allele frequencies. Of the 67 genetic markers, 65 MHs were, as far as we know, newly discovered, and 32 MHs had effective number of allele (A_e) values greater than 5.0. The average A_e and heterozygosity of the panel were 5.34 and 0.7352, respectively. Next, 53 MHs from a previous study were collected as Panel A (average A_e of 7.43), and Panel C with 87 MHs (average A_e of 7.02) was formed by combining Panels A and B. We investigated the utility of these three panels in kinship analysis (parent-child, full siblings, 2nd-degree, 3rd-degree, 4th-degree, and 5th-degree relatives), with Panel C exhibiting better performance than the two other panels. Panel C was able to separate parent-child, full-sibling, and 2nd-degree relative duos from unrelated controls in real pedigree data, with a small false testing level (FTL) of 0.11% in simulated 2nd-degree duos. For more distant relationships, the FTL was much higher: 8.99% for 3rd-degree, 35.46% for 4th-degree, and 61.55% for 5th-degree. When a carefully chosen extra relative was known, this may enhance the testing power for distant kinship analysis. Two twins from the Q family (2-5 and 2-7) and W family (3-18 and 3-19) shared the same genotypes in all tested MHs, which led to the incorrect conclusion that an uncle-nephew duo was classified as a parent-child duo. In addition, Panel C showed great capacity for excluding close relatives (2nd-degree and 3rd-degree relatives) during paternity tests. Among 18,246 real and 10,000 simulated unrelated pairs, none were misinterpreted as a relative within 2nd-degree at a log₁₀(LR) cutoff of 4. The panels presented herein could provide supplementary power for the analysis of complex kinship.

Case Report: A Novel Homozygous Variant of the SERPINF1 Gene in Rare Osteogenesis Imperfecta Type VI

Osteogenesis imperfecta (OI) is a group of connective tissue disorders with different types of inheritance. OI is characterized by bone fragility and deformities, frequent fractures, low bone-mineral density, and impaired bone micro-architectonics. We described here a case of a one-year-old Tuvan patient with multiple fractures. The disease manifestation occurred first at 12 weeks of age as a shoulder joint bruise, and during the year, the patient sustained 27 fractures. Genetic testing revealed a novel homozygous mutation, c.259_260insCGGCC (p.T87fs), in the SERPINF1 gene. This insertion leads to an open-reading frameshift, and the mutation is not represented in the databases. Mutations in SERPINF1 lead to type VI OI, the clinical picture of which is similar to the disease phenotype manifestation of the patient. Thus, the patient’s diagnosis was established by finding a novel pathogenic variant in the SERPINF1 gene.

Applying Unique Molecular Indices with an Extensive All-in-One Forensic SNP Panel for Improved Genotype Accuracy and Sensitivity

One of the major challenges in forensic genetics is being able to detect very small amounts of DNA. Massively parallel sequencing (MPS) enables sensitive detection; however, genotype errors may exist and could interfere with the interpretation. Common errors in MPS-based analysis are often induced during PCR or sequencing. Unique molecular indices (UMIs) are short random nucleotide sequences ligated to each template molecule prior to amplification. Applying UMIs can improve the limit of detection by enabling accurate counting of initial template molecules and removal of erroneous data. In this study, we applied the FORCE panel, which includes ~5500 SNPs, with a QIAseq Targeted DNA Custom Panel (Qiagen), including UMIs. Our main objective was to investigate whether UMIs can enhance the sensitivity and accuracy of forensic genotyping and to evaluate the overall assay performance. We analyzed the data both with and without the UMI information, and the results showed that both genotype accuracy and sensitivity were improved when applying UMIs. The results showed very high genotype accuracies (>99%) for both reference DNA and challenging samples, down to 125 pg. To conclude, we show successful assay performance for several forensic applications and improvements in forensic genotyping when applying UMIs.

Comprehensive evaluations of individual discrimination, kinship analysis, genetic relationship exploration and biogeographic origin prediction in Chinese Dongxiang group by a 60-plex DIP panel

BACKGROUND

Dongxiang group, as an important minority, resides in Gansu province which is located at the northwest China, forensic detection system with more loci needed to be studied to improve the application efficiency of forensic case investigation in this group.

METHODS

A 60-plex system including 57 autosomal deletion/insertion polymorphisms (A-DIPs), 2 Y chromosome DIPs (Y-DIPs) and the sex determination locus (Amelogenin) was explored to evaluate the forensic application efficiencies of individual discrimination, kinship analysis and biogeographic origin prediction in Gansu Dongxiang group based on the 60-plex genotype results of 233 unrelated Dongxiang individuals. The 60-plex genotype results of 4582 unrelated individuals from 33 reference populations in five different continents were also collected to analyze the genetic background of Dongxiang group and its genetic relationships with other continental populations.

RESULTS

The system showed high individual discrimination power, as the cumulative power of discrimination (CPD), cumulative power of exclusion (CPE) for trio and cumulative match probability (CMP) values were 0.99999999999999999999997297, 0.999980 and 2.7029E^- 24, respectively. The system could distinguish 98.12%, 93.78%, 82.18%, 62.35% and 39.32% of full sibling pairs from unrelated individual pairs, when the likelihood ratio (LR) limits were set as 1, 10, 100, 1000 and 10,000 based on the simulated family samples, respectively. Additionally, Dongxiang group had the close genetic distances with populations in East Asia, especially showed the intimate genetic relationships with Chinese Han populations, which were concluded from the genetic affinities and genetic background analyses of Dongxiang group and 33 reference populations. In terms of the effectiveness of biogeographic origin inference, different artificial intelligent algorithms possessed different efficacies. Among them, the random forest (RF) and extreme gradient boosting (XGBoost) algorithm models could accurately predict the biogeographic origins of 99.7% and 90.59% of three and five continental individuals, respectively.

CONCLUSION

This 60-plex system had good performance for individual discrimination, kinship analysis and biogeographic origin prediction in Dongxiang group, which could be used as a powerful tool for case investigation.

Improved computations for relationship inference using low-coverage sequencing data

Pedigree inference, for example determining whether two persons are second cousins or unrelated, can be done by comparing their genotypes at a selection of genetic markers. When the data for one or more of the persons is from low-coverage next generation sequencing (lcNGS), currently available computational methods either ignore genetic linkage or do not take advantage of the probabilistic nature of lcNGS data, relying instead on first estimating the genotype. We provide a method and software (see familias.name/lcNGS) bridging the above gap. Simulations indicate how our results are considerably more accurate compared to some previously available alternatives. Our method, utilizing a version of the Lander-Green algorithm, uses a group of symmetries to speed up calculations. This group may be of further interest in other calculations involving linked loci.

Forensic efficiencies of individual identification, kinship testing and ancestral inference in three Yunnan groups based on a self-developed multiple DIP panel

Deletion/insertion polymorphism (DIP), as a short insertion/deletion sequence polymorphic genetic marker, has attracted the attention of forensic genetic scientist due to its lack of stutter, short amplicon and abundant ancestral information. In this study, based on a self-developed 43 autosomal deletion/insertion polymorphism (A-DIP) loci panel which could meet the forensic application purposes of individual identification, kinship testing and ancestral inference to some extent, we evaluated the forensic efficiencies of the above three forensic objectives in Chinese Yi, Hani and Miao groups of Yunnan province. The cumulative match probability (CPM) and combined probability of exclusion (CPE) of these three groups were 1.11433E-18, 8.24299E-19, 4.21721E-18; 0.999610217, 0.999629285 and 0.999582084, respectively. Average 96.65% full sibling pairs could be identified from unrelated individual pairs (as likelihood ratios > 1) using this DIP panel, whereas the average false positive rate was 3.69% in three target Yunnan groups. With the biogeographical ancestor prediction models constructed by extreme gradient boosting (XGBoost) and support vector machine (SVM) algorithms, 0.8239 (95% CI 0.7984, 0.8474) of the unrelated individuals could be correctly divided according to the continental origins based on the 43 A-DIPs which were large frequency distribution differentiations among different continental populations. The present results of principal component analysis (PCA), multidimensional scaling (MDS), neighbor joining (NJ) and maximum likelihood (ML) phylogenetic trees and STRUCTURE analyses indicated that these three Yunnan groups had relatively close genetic distances with East Asian populations.

Domestic pigs (Sus scrofa) engage in non-random post-conflict affiliation with third parties: cognitive and functional implications

In social mammals, conflict resolution involves the reunion of former opponents (aggressor and victim) after an aggressive event (reconciliation) or post-conflict triadic contacts with a third party, started by either opponent (solicited-TSC) or spontaneously offered by the third party (unsolicited-TUC). These post-conflict strategies can serve different functions, including consolation (specifically when TUCs reduce the victim's anxiety). We investigated the possible presence and modulating factors of such strategies on semi-free ranging pigs (Sus scrofa; N = 104), housed at the ethical farm Parva Domus (Cavagnolo, Italy). Kinship was known. Reconciliation was present and mainly occurred between weakly related pigs to possibly improve tolerant cohabitation. Triadic contacts (all present except aggressor TSCs) mostly occurred between close kin. TSCs enacted by victims reduced neither their post-conflict anxiety behaviors nor further attacks by the previous aggressor, possibly because TSCs remained largely unreciprocated. TUCs towards aggressors did not reduce aggressor post-conflict anxiety but limited aggression redirection towards third parties. TUCs towards the victim reduced the victim but not the third-party's anxiety. However, TUCs may also provide inclusive fitness benefits to third parties by benefiting close kin. In sum, pigs engaged in non-random solicited/unsolicited triadic contacts, which suggests that pigs might possess socio-emotional regulation abilities to change their own or others' experience and elements of social appraisal, necessary to detect the emotional arousal of relevant others and (in case of TUCs) take the agency to restore homeostasis.

A sequence-based 163plex microhaplotype assay for forensic DNA analysis

Novel genetic marker microhaplotype has led to an upsurge in forensic genetic research. This study established a 163 microhaplotype (MH) multiplex assay based on next-generation sequencing (NGS) and evaluated the assay’s performance and applicability. Our results showed that the 163 MH assay was accurate, repeatable and reliable, and could distinguish between African, European-American, Southern Asia and Eastern Asia populations. Among the 163 MH makers, 48 MHs with Ae > 3.0 in China Eastern Han were selected and confirmed to be highly polymorphic, with a combined power of discrimination of 1–8.26 × 10–44 and the combined power of exclusion in duos and trios of 1–1.26 × 10–8 and 1–8.27 × 10–16, respectively. Moreover, the mixture study demonstrated the realizability of the MHs in deconvoluting mixtures with different proportions of two to five-person. In conclusion, our findings support the use of this MH assay for ancestry inference, human identification, paternity testing and mixture deconvolution in forensic research.

The application of short and highly polymorphic microhaplotype loci in paternity testing and sibling testing of temperature-dependent degraded samples

Paternity testing and sibling testing become more complex and difficult when samples degrade. But the commonly used genetic markers (STR and SNP) cannot completely solve this problem due to some disadvantages. The novel genetic marker microhaplotype proposed by Kidd’s research group combines the advantages of STR and SNP and is expected to become a promising genetic marker for kinship testing in degraded samples. Therefore, in this study, we intended to select an appropriate number of highly polymorphic SNP-based microhaplotype loci, detect them by the next-generation sequencing technology, analyze their ability to detect degraded samples, calculate their forensic parameters based on the collected 96 unrelated individuals, and evaluate their effectiveness in paternity testing and sibling testing by simulating kinship relationship pairs, which were also compared to 15 STR loci. Finally, a short and highly polymorphic microhaplotype panel was developed, containing 36 highly polymorphic SNP-based microhaplotype loci with lengths smaller than 100 bp and Ae greater than 3.00, of which 29 microhaplotype loci could not reject the Hardy-Weinberg equilibrium and linkage equilibrium after the Bonferroni correction. The CPD and CPE of these 29 microhaplotype loci were 1-2.96E-26 and 1-5.45E-09, respectively. No allele dropout was observed in degraded samples incubated with 100°C hot water for 40min and 60min. According to the simulated kinship analysis, the effectiveness at the threshold of 4/−4 reached 98.39% for relationship parent-child vs. unrelated individuals, and the effectiveness at the threshold of 2/−2 for relationship full-sibling vs. unrelated individuals was 93.01%, which was greater than that of 15 STR loci (86.75% for relationship parent-child vs. unrelated individuals and 81.73% for relationship full-sibling vs. unrelated individuals). After combining our 29 microhaplotype loci with other 50 short and highly polymorphic microhaplotype loci, the effectiveness values at the threshold of 2/−2 were 82.42% and 90.89% for relationship half-sibling vs. unrelated individuals and full-sibling vs. half-sibling. The short and highly polymorphic microhaplotype panel we developed may be very useful for paternity testing and full sibling testing in degraded samples, and in combination with short and highly polymorphic microhaplotype loci reported by other researchers, may be helpful to analyze more distant kinship relationships.

Forensic DNA Typing From Femurs and Bones of the Foot: A Study of 3 Cases

ABSTRACT

Evidence has been accumulating in the sense that femur may not always be the best option for DNA typing of skeletal remains. Recent studies have shown that bones of the hands and feet appear to be a superior source of preserved DNA. The current study reanalyzed DNA quantitation, degradation, and short tandem repeat typing in femurs, lateral cuneiforms, and distal foot phalanges. Data from 3 human identification cases involving corpses in an advanced decomposition state were collected. We found that in the studied cases, the femur provided equal or inferior results, recovering 84.9% of true alleles. Lateral cuneiforms (99.2%) and distal foot phalanges (96.8%) yielded higher percentages. In addition, more drop-ins and drop-outs were detected in femurs than cuneiforms and phalanges. This study adds to current findings that advocate for further investigation into bone selection for use in forensic practice. The impacts of our findings are limited by the small number of individuals studied and may not apply to old and degraded bones.

Estimations of Mutation Rates Depend on Population Allele Frequency Distribution: The Case of Autosomal Microsatellites

Microsatellites (or short-tandem repeats (STRs)) are widely used in anthropology and evolutionary studies. Their extensive polymorphism and rapid evolution make them the ideal genetic marker for dating events, such as the age of a gene or a population. This usage requires the estimation of mutation rates, which are usually estimated by counting the observed Mendelian incompatibilities in one-generation familial configurations (typically parent(s)–child duos or trios). Underestimations are inevitable when using this approach, due to the occurrence of mutational events that do not lead to incompatibilities with the parental genotypes (‘hidden’ or ‘covert’ mutations). It is known that the likelihood that one mutation event leads to a Mendelian incompatibility depends on the mode of genetic transmission considered, the type of familial configuration (duos or trios) considered, and the genotype(s) of the progenitor(s). In this work, we show how the magnitude of the underestimation of autosomal microsatellite mutation rates varies with the populations’ allele frequency distribution spectrum. The Mendelian incompatibilities approach (MIA) was applied to simulated parent(s)/offspring duos and trios in different populational scenarios. The results showed that the magnitude and type of biases depend on the population allele frequency distribution, whatever the type of familial data considered, and are greater when duos, instead of trios, are used to obtain the estimates. The implications for molecular anthropology are discussed and a simple framework is presented to correct the naïf estimates, along with an informatics tool for the correction of incompatibility rates obtained through the MIA.

Extended population genetic analysis of 12 X-STRs - Exemplified using a Norwegian population sample

The use of X-chromosomal markers to resolve questions of relatedness has experienced a significant increase during the last years in forensic genetics. Perhaps primarily due to the emergence of commercial kits, but equally important due to an increased awareness of the utility of those markers. The X-chromosomal inheritance pattern entails that some cases, for instance paternal half-sisters, can potentially be resolved using a few X-chromosomal markers alone. For the statistical assessment in kinship cases it is of importance to have relevant population frequency data. In the present study 631 unrelated males from a Norwegian population sample are analyzed. The resulting haplotypes are compared to previously studied population samples and a deeper analysis of the linkage disequilibrium (LD) structure is conducted. We demonstrate that the power to detect LD will be low when few males, say below 300, are analyzed. We use entropy to describe the degree of LD between multiallelic loci and describe how this measure varies between different studied populations. Large population frequency databases have been recommended when using X-chromosomal markers, and we show that by combining reference databases from genetically similar populations, more precise haplotype frequency estimates can be obtained for rare haplotypes which improves the statistical assessment of the weight of evidence. In addition, we promote the use of simulations to assess the utility of STR markers in contrast to standard forensic parameters. Specifically we perform extensive simulations on cases where X-chromosomal markers are important and illustrate how the results can be used to infer the information gained from these markers.

Forensic Feature Exploration and Comprehensive Genetic Insights Into Yugu Ethnic Minority and Northern Han Population via a Novel NGS-Based Marker Set

The MPS technology has expanded the potential applications of DNA markers and increased the discrimination power of the targeted loci by taking variations in their flanking regions into consideration. Here, a collection of nuclear and extranuclear DNA markers (totally six kinds of nuclear genetic markers and mtDNA hypervariable region variations) were comprehensively and systematically assessed for polymorphism detections, further employed to dissect the population backgrounds in the Yugu ethnic group from Gansu province (Yugu) and Han population from the Inner Mongolia Autonomous Region (NMH) of China. The elevated efficiencies of the marker set in separating full sibling and challenging half sibling determination cases in parentage tests (iiSNPs), as well as predicting ancestry origins of unknown individuals from at least four continental populations (aiSNPs) and providing informative characteristic-related clues for Chinese populations (piSNPs) are highlighted in the present study. To sum up, different sets of DNA markers revealed sufficient effciencies to serve as promising tools in forensic applications. Genetic insights from the perspectives of autosomal DNA, Y chromosomal DNA, and mtDNA variations yielded that the Yugu ethnic group was genetically close related to the Han populations of the northern region. But we admit that more reference populations (like Mongolian, Tibetan, Hui, and Tu) should be incorporated to gain a refined genetic background landscape of the Yugu group in future studies.

A multi-dimensional evaluation of the 'NIST 1032' sample set across four forensic Y-STR multiplexes

This manuscript reports Y-chromosomal short tandem repeat (Y-STR) haplotypes for 1032 male U.S. population samples across 30 Y-STR loci characterized by three capillary electrophoresis (CE) length-based kits (PowerPlex Y23 System, Yfiler Plus PCR Amplification Kit, and Investigator Argus Y-28 QS Kit) and one sequence-based kit (ForenSeq DNA Signature Prep Kit): DYF387S1, DYS19, DYS385 a/b, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, DYS448, DYS449, DYS456, DYS458, DYS460, DYS481, DYS505, DYS518, DYS522, DYS533, DYS549, DYS570, DYS576, DYS612, DYS627, DYS635, DYS643, and Y-GATA-H4. The length-based Y-STR haplotypes include six loci that are not reported in the sequence-based kit (DYS393, DYS449, DYS456, DYS458, DYS518, and DYS627), whereas three loci included in the sequence-based kit are not present in length-based kits (DYS505, DYS522, and DYS612). For the latter, a custom multiplex was used to generate CE length-based data, allowing 1032 samples to be evaluated for concordance across the 30 Y-STR loci included in these four commercial Y-STR typing kits. Discordances between typing methods were analyzed further to assess underlying causes such as primer binding site mutations and flanking region insertions/deletions. Allele-level frequency and statistical information is provided for sequenced loci, excluding the multi-copy loci DYF387S1 and DYS385 a/b, for which locus-specific haplotype-level frequencies are provided instead. The resulting data reveals the degree of information gained through sequencing: 88% of sequenced Y-STR loci contain additional sequence-based alleles compared to length-based data, with the DYS389II locus containing the most additional alleles (51) observed by sequencing. Despite these allelic increases, only minimal improvement was observed in haplotype resolution by sequence, with all four commercial kits providing a similar ability to differentiate length-based haplotypes in this sample set. Finally, a subset of 369 male samples were compared to their corresponding additionally sequenced father samples, revealing the sequence basis for the 50 length-based changes observed, and no additional sequence-based mutations. GenBank accession numbers are reported for each unique sequence, and associated records are available in the STRSeq Y-Chromosomal STR Loci National Center for Biotechnology Information (NCBI) BioProject, accession PRJNA380347. Haplotype data is updated in the Y-STR Haplotype Reference Database (YHRD) for the 'NIST 1032' data set to now achieve the level of maximal haplotype of YHRD. All supplementary files including revisions to previously published Y-STR data are available in the NIST Public Data Repository: U.S. population data for human identification markers, DOI 10.18434/t4/1500024.

Kinship assignment with the ForenSeq™ DNA Signature Prep Kit: Sources of error in simulated and real cases

Kinship recognition between anonymous DNA samples is becoming a relevant issue in forensics, more so with the increasing number of DNA profiles in databanks. Also, NGS-based genotyping is being increasingly used in routine personal identification, to simultaneously type large numbers of markers of different kind. In the present work, we explored computationally and experimentally the performance of the ForenSeq™ DNA Signature Prep Kit in identifying the true relationship between two anonymous samples, distinguishing it from other possible relationships. We analyzed with Familias R series of 10,000 pairs with 9 different simulated relationships, corresponding to different degrees of autosomal sharing. For each pair we obtained likelihood ratios for five kinship hypotheses vs. unrelatedness, and used their ranking to identify the preferred relationship. We also typed 21 subjects from two pedigrees, representing from parent-child to 4th cousins relationships. As expected, the power for identifying the true relationship decays in the order of autosomal sharing. Parent-child and full siblings can be robustly identified against other relationships. For half-siblings the chance of reaching a significant conclusion is already small. For more distant relationships the proportion of cases correctly and significantly identified is 10% or less. Bidirectional errors in kinship attribution include the suggestion of relatedness when this does not exist (10-50%), and the suggestion of independence in pairs of individuals more than 4 generations apart (25-60%). The real cases revealed a relevant effect of genotype miscalling at some loci, which could only be partly avoided by modulating the analysis parameters. In conclusion, with the exception of first degree relatives, the kit can be useful to inform additional investigations, but does not usually provide probatory results.

Frequencies of D19S433 silent alleles in a Japanese population of 1501 individuals and their effect on likelihood ratios calculated in kinship tests

Although silent alleles in D19S433 typing using the GlobalFiler PCR Amplification Kit have been reported, the exact frequency of the D19S433 silent alleles in population data of 1501 Japanese individuals, which are widely used for the assessment of Japanese STR typing results, is unclear. In this study, we examined the exact D19S433 silent allele frequency in this population data. We newly observed the G32A variant causing silent alleles at D19S433 in five samples. Combining them with data including 30 samples with the variant reported previously, we determined that the total frequency of the silent alleles (i.e. the frequency of the G32A variant) in the 1501 Japanese samples was 0.0117 (35/3002). Using the D19S433 allele frequency data, we evaluated the effect of presence/absence information for the D19S433 silent allele on kinship tests. Likelihood ratios (LRs) were calculated for both simulated parent-child and full sibling cases, revealing that the LR may change by approximately 10^-2 to 10³ fold when the presence/absence of the D19S433 silent allele is revealed in a kinship test. Therefore, if a sufficiently large or small LR is obtained, there is little need to determine the presence/absence of the D19S433 silent allele in Japanese kinship tests using GlobalFiler. This study will be beneficial for the assessment of Japanese human identification and kinship test results using GlobalFiler.

Evaluation of diallelic STR markers with inter-population allelic database for their usefulness in paternity trios in the Central Indian population

BACKGROUND

Most of the forensic DNA laboratories have migrated to new generation STR kits of 6 dye chemistry with more than 20 autosomal STRs. The population-specific databases of such STR markers are lacking in many regions.

AIM

To evaluate the effect of the inter-population database in 100 paternity trios with no inconsistencies at 23 STRs.

SUBJECTS AND METHODS

100 paternity trios were evaluated considering inter-population allelic frequency values for calculation of Combined Paternity Index (CPI) and Probability of Paternity (POP).

RESULTS

No significant variation (p < 0.05) among the allele frequencies at the interpopulation level was observed. The number of obligate alleles and the likelihood of transferring obligate alleles from the putative father showed a positive correlation (p < 0.005) with Power of Discrimination (PD), Polymorphic Information Content (PIC), Power of Exclusion (PE), Paternity Index (PI), Observed and Expected Heterozygosity (Ho and He), and a statistically significant negative correlation (p < 0.005) with Matching Probability (Pm). The average Combined Paternity Index (CPI) and Probability of Paternity (POP) did not show any statistically significant difference (p < 0.05) at the interpopulation level.

CONCLUSION

The allelic database showed no effect on the CPI and POP in 100 paternity trios. This suggests no urgent need for using population-specific databases for statistical evaluation of paternity trios without inconsistencies.

19th century family saga re-told by DNA recovered from postcard stamps

Old postcards with stamps might help unravelling historical family stories and relationships. By employing ancient DNA recovered from world war I postage stamps, we disprove a family saga of an illegitimate child born in 1887. We developed a protocol to collect DNA from saliva, trapped and protected on the backside of postage stamps glued on postcards. With replicate STR analyses we were able to assemble almost full autosomal and Y-STR profiles of three male, deceased family members. The illegitimate child turned out to be a legitimate child of a later married couple.

Pairwise kinship analysis of 17 pedigrees using massively parallel sequencing

With the tremendous development of massively parallel sequencing (MPS) in the last decade, it has been widely applied in basic science, clinical diagnostics, microbial genomics, as well as forensic genetics. MPS has lots of advantages that may facilitate the kinship analysis. In this study, 243 Chinese Han individuals from 17 families were involved and sequenced using the ForenSeq™ DNA Signature Prep Kit (Verogen, Inc., San Diego, USA), which provided the sequence information of 27 autosomal STRs (A-STRs), 7 X chromosomal STRs (X-STRs), 24 Y chromosomal STRs (Y-STRs) and 94 identity-informative SNPs (iSNPs). A total of 275 pairs of parent-child, 123 pairs of full siblings, 1 pair of twins, 1 pair of half siblings, 158 pairs of grandparent-grandchild, 222 pairs of uncle/aunt-nephew/niece and 121 pairs of first cousins, as well as 701 pairs of unrelated individuals were identified. Using both likelihood ratio (LR) and identical by state (IBS) methods, the kinship analysis was conducted among these relative and non-relative pairs based on the A-STRs and SNPs. As a result, the ForenSeq Signature Kit could solve the analysis of parent-child (t1 = -4, t2 = 4), full siblings (t1 = -2, t2 = 2) and most second-degree kinships (t1 = -1, t2 = 1) using the LR method. When the IBS method was applied, 123 full sibling pairs had a higher average IBS value than other kinship groups in this study. And the IBS method could play a role in the testing of parent-child and full siblings.

Streamlining the decision-making process for international DNA kinship matching using Worldwide allele frequencies and tailored cutoff log10LR thresholds

The identification of human remains belonging to missing persons is one of the main challenges for forensic genetics. Although other means of identification can be applied to missing person investigations, DNA is often extremely valuable to further support or refute potential associations. When reference DNA samples cannot be collected from personal items belonging to a missing person, a direct DNA identification cannot be carried out. However, identifications can be made indirectly using DNA from the missing person's relatives. The ranking of likelihood ratio (LR) values, which measure the fit of a missing person for any given pedigree, is often the first step in selecting candidates in a DNA database. Although implementing DNA kinship matching in a national environment is feasible, many challenges need to be resolved before applying this method to an international configuration. In this study, we present an innovative and intuitive method to perform international DNA kinship matching and facilitate the comparison of DNA profiles when the ancestry is unknown or unsure and/or when different marker sets are used. This straightforward method, which is based on calculations performed with the DNA matching software BONAPARTE, Worldwide allele frequencies and tailored cutoff log₁₀LR thresholds, allows for the classification of potential candidates according to the strength of the DNA evidence and the predicted proportion of adventitious matches. This is a powerful method for streamlining the decision-making process in missing person investigations and DVI processes, especially when there are low numbers of overlapping typed STRs. Intuitive interpretation tables and a decision tree will help strengthen international data comparison for the identification of reported missing individuals discovered outside their national borders.

Precision DNA Mixture Interpretation with Single-Cell Profiling

Wet-lab based studies have exploited emerging single-cell technologies to address the challenges of interpreting forensic mixture evidence. However, little effort has been dedicated to developing a systematic approach to interpreting the single-cell profiles derived from the mixtures. This study is the first attempt to develop a comprehensive interpretation workflow in which single-cell profiles from mixtures are interpreted individually and holistically. In this approach, the genotypes from each cell are assessed, the number of contributors (NOC) of the single-cell profiles is estimated, followed by developing a consensus profile of each contributor, and finally the consensus profile(s) can be used for a DNA database search or comparing with known profiles to determine their potential sources. The potential of this single-cell interpretation workflow was assessed by simulation with various mixture scenarios and empirical allele drop-out and drop-in rates, the accuracies of estimating the NOC, the accuracies of recovering the true alleles by consensus, and the capabilities of deconvolving mixtures with related contributors. The results support that the single-cell based mixture interpretation can provide a precision that cannot beachieved with current standard CE-STR analyses. A new paradigm for mixture interpretation is available to enhance the interpretation of forensic genetic casework.

A Swiss collaborative exercise for Disaster Victim Identification (DVI): Covering situations with different levels of complexity

The identification of victims of a disaster (DVI) requires the collaboration of different specialists. Within a DVI context, DNA analyses often play an important role. Consequently, forensic genetic laboratories should be prepared to cope with DVI situations, as this can involve large-scale DNA profile comparisons. Six forensic genetic laboratories from Switzerland participated in an exercise where supposedly a plane had crashed. The goal of the exercise was to monitor participants use of dedicated software with ground truth cases and to make them aware of the existence of particular situations that may occur in real cases. For assigning the value of the comparison of the DNA profiles, all participating laboratories used the DVI module of Familias v3.2.¹ In addition, one of the 6 laboratories used the Pedigree Searcher from CODIS v7.0. The data (AmpFlSTR® NGM SElect™ profiles) were generated to challenge the participating laboratories: cases with first, second degree biological parents, mutation events, as well as non-paternity cases were included. This study shows that the majority of the participants used the software in an appropriate way. However, a few misleading conclusions were detected for the most challenging situations. These errors belonged to one of the following categories: false pedigree, false association using the higher LR, misleading contextual information (false paternity) and not clustering family members. Specific recommendations are provided in order to reduce misuse of the software and the risk of misinterpretations by using all the relevant information.

Interpretation of DNA data within the context of UK forensic science - investigation

This article is the second part of a review of the interpretation of DNA data in forensic science. The first part describes the evaluation of autosomal profile for criminal trials where an evidential weight is assigned to the profile of a person of interest (POI) and a crime-scene profile. This part describes the state of the art and future advances in the interpretation of forensic DNA data for providing intelligence information during an investigation. Forensic DNA is crucial in the investigative phase of an undetected crime where a POI needs to be identified. A sample taken from a crime scene is profiled using a range of forensic DNA tests. This review covers investigation using autosomal profiles including searching national and international crime and reference DNA databases. Other investigative methodologies described are kinship analysis; familial searching; Y chromosome (Y-STR) and mitochondrial (mtDNA) profiles; appearance prediction and geographic ancestry; forensic genetic genealogy; and body identification. For completeness, the evaluation of Y-STRs, mtDNA and kinship analysis are briefly described. Taken together, parts I and II, cover the range of interpretation of DNA data in a forensic context.

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.

Strategies for pairwise searches in forensic kinship analysis

Testing kinship between pairs of individuals is central to a wide range of applications. We focus on cases where many tests are done jointly. Typical examples include cases where DNA profiles are available from a burial site, a plane crash or a database of convicted offenders. The task is to determine the relationships between DNA profiles or individuals. Our approach generalises previous methods and implementations in several respects. We model general, possibly inbred, pairwise relationships which is important for non-human applications and in archaeological studies of ancient inbred populations. Furthermore, we do not restrict attention to autosomal markers. Some cases, such as distinguishing between maternal and paternal half siblings, can be solved using X-chromosomal markers. When many tests are done, the risk of errors increases. We address this problem by building on the theory of multiple testing and show how optimal thresholds for tests can be determined. We point out that the likelihood ratios in a blind search may be dependent so multiple testing methods and interpretation need to account for this. In addition, we show how a Bayesian approach can be helpful. Our examples, using simulated and real data, demonstrate the practical importance of the methods and implementation is based on freely available software.

Joint DNA-based disaster victim identification

We address computational and statistical aspects of DNA-based identification of victims in the aftermath of disasters. Current methods and software for such identification typically consider each victim individually, leading to suboptimal power of identification and potential inconsistencies in the statistical summary of the evidence. We resolve these problems by performing joint identification of all victims, using the complete genetic data set. Individual identification probabilities, conditional on all available information, are derived from the joint solution in the form of posterior pairing probabilities. A closed formula is obtained for the a priori number of possible joint solutions to a given DVI problem. This number increases quickly with the number of victims and missing persons, posing computational challenges for brute force approaches. We address this complexity with a preparatory sequential step aiming to reduce the search space. The examples show that realistic cases are handled efficiently. User-friendly implementations of all methods are provided in the R package dvir, freely available on all platforms.