Review of SNP assays for disaster victim identification: Cost, time, and performance information for decision-makers.

In mass disaster events, forensic DNA laboratories may be called upon to quickly pivot their operations toward identifying bodies and reuniting remains with family members. Ideally, laboratories have considered this possibility in advance and have a plan in place. Compared with traditional short tandem repeat (STR) typing, single nucleotide polymorphisms (SNPs) may be better suited to these disaster victim identification (DVI) scenarios due to their small genomic target size, resulting in an improved success rate in degraded DNA samples. As the landscape of technology has shifted toward DNA sequencing, many forensic laboratories now have benchtop instruments available for massively parallel sequencing (MPS), facilitating this operational pivot from routine forensic STR casework to DVI SNP typing. Herein, we present the commercially available SNP sequencing assays amenable to DVI, we use data simulations to explore the potential for kinship prediction from SNP panels of varying sizes, and we give an example DVI scenario as context for presenting the matrix of considerations: kinship predictive potential, cost, and throughput of current SNP assay options. This information is intended to assist laboratories in choosing a SNP system for disaster preparedness.

Harmonizing the forensic nomenclature for STR loci D6S474 and DYS612.

The autosomal STR D6S474 and the Y-chromosomal STR DYS612 have been reported in multiple ways in the forensic literature, with differences in both the bracketed repeat structures and counting of numerical length-based capillary electrophoresis (CE) alleles. These issues often come to light when STR loci are introduced in commercial assays and results compared with historical publications of allele frequency data, or multiple assays are characterized with reference materials. We review the forensic literature and other relevant information, and provide suggestions for the future treatment of each STR.

Report from the STRAND Working Group on the 2019 STR sequence nomenclature meeting.

This report summarizes topics discussed at the STR sequence nomenclature meeting hosted by the STRAND Working Group in April 2019. Invited attendees for this meeting included researchers known-to-us to be developing STR sequence-based nomenclature schemata, scientific representatives from vendors developing STR sequence bioinformatic methods, DNA intelligence database curators, and academic experts in STR genomics. The goal of this meeting was to provide a forum for individuals developing nomenclature schemata to present and discuss their ideas, encouraging mutual awareness, identification of differences in approaches, opposing aspects, and opportunities for parallelization while some approaches are still under development.

Estimating number of contributors in massively parallel sequencing data of STR loci.

In recent years a number of computer-based algorithms have been developed for the deconvolution of complex DNA mixtures in forensic science. These procedures utilize likelihood ratios that quantify the evidence for a hypothesis for the presence of a person of interest in a DNA profile compared to an alternative hypothesis. Proper operation of these software systems requires an assumption regarding the total number of contributors present in the mixture. Unfortunately, estimates based on counting the number of alleles at a locus can be inaccurate due to the sharing and masking of alleles at individual loci. The effects of allele masking become increasingly severe as the number of contributors increases, rendering estimates about high-order mixtures uncertain. The accuracy of these estimates can be improved by increasing the number of STR markers in panels, and by using highly polymorphic markers. Increasing the number of STR markers from 13 to 20 (expanded CODIS panel) improves the accuracy of allele count-based estimation methods for low-order mixtures, but accuracy for high-order mixtures (> 3 contributors) remains poor due to allele masking. An alternative technique, massively parallel sequencing, holds great potential to improve the accuracy of the estimate of number of contributors due to its ability to detect sequence polymorphisms within alleles. This process results in an expansion of the number of alleles when compared to that obtained using capillary electrophoresis. Here, we show that the detection of these additional sequence-defined alleles in 22-marker panels improves number of contributor estimates in conceptual mixtures of 4 and 5 contributors.

Sequence-based U.S. population data for 27 autosomal STR loci.

This manuscript reports Short Tandem Repeat (STR) sequence-based allele frequencies for 1036 samples across 27 autosomal STR loci: D1S1656, TPOX, D2S441, D2S1338, D3S1358, D4S2408, FGA, D5S818, CSF1PO, D6S1043, D7S820, D8S1179, D9S1122, D10S1248, TH01, vWA, D12S391, D13S317, Penta E, D16S539, D17S1301, D18S51, D19S433, D20S482, D21S11, Penta D, and D22S1045. Sequence data were analyzed by two bioinformatic pipelines and all samples have been evaluated for concordance with alleles derived from CE-based analysis at all loci. Each reported sequence includes high-quality flanking sequence and is properly formatted according to the most recent guidance of the International Society for Forensic Genetics. In addition, GenBank accession numbers are reported for each sequence, and associated records are available in the STRSeq BioProject (https://www.ncbi.nlm.nih.gov/bioproject/380127). The D3S1358 locus demonstrates the greatest average increase in heterozygosity across populations (approximately 10 percentage points). Loci demonstrating average increase in heterozygosity from 10 to 5 percentage points include (in descending order) D9S1122, D13S317, D8S1179, D21S11, D5S818, D12S391, and D2S441. The remaining 19 loci each demonstrate less than 5 percentage point increase in average heterozygosity. Discussion includes the utility of this data in understanding traditional CE results, such as informing stutter models and understanding migration challenges, and considerations for population sampling strategies in light of the marked increase in rare alleles for several of the sequence-based STR loci. This NIST 1036 data set is expected to support the implementation of STR sequencing forensic casework by providing high-confidence sequence-based allele frequencies for the same sample set which are already the basis for population statistics in many U.S. forensic laboratories.

STRSeq: A catalog of sequence diversity at human identification Short Tandem Repeat loci.

The STR Sequencing Project (STRSeq) was initiated to facilitate the description of sequence-based alleles at the Short Tandem Repeat (STR) loci targeted in human identification assays. This international collaborative effort, which has been endorsed by the ISFG DNA Commission, provides a framework for communication among laboratories. The initial data used to populate the project are the aggregate alleles observed in targeted sequencing studies across four laboratories: National Institute of Standards and Technology (N=1786), Kings College London (N=1043), University of North Texas Health Sciences Center (N=839), and University of Santiago de Compostela (N=944), for a total of 4612 individuals. STRSeq data are maintained as GenBank records at the U.S. National Center for Biotechnology Information (NCBI), which participates in a daily data exchange with the DNA DataBank of Japan (DDBJ) and the European Nucleotide Archive (ENA). Each GenBank record contains the observed sequence of a STR region, annotation ("bracketing") of the repeat region and flanking region polymorphisms, information regarding the sequencing assay and data quality, and backward compatible length-based allele designation. STRSeq GenBank records are organized within a BioProject at NCBI (https://www.ncbi.nlm.nih.gov/bioproject/380127), which is sub-divided into: commonly used autosomal STRs, alternate autosomal STRs, Y-chromosomal STRs, and X-chromosomal STRs. Each of these categories is further divided into locus-specific BioProjects. The BioProject hierarchy facilitates access to the GenBank records by browsing, BLAST searching, or ftp download. Future plans include user interface tools at strseq.nist.gov, a pathway for submission of additional allele records by laboratories performing population sample sequencing and interaction with the STRidER web portal for quality control (http://strider.online).

Sequence variation of 22 autosomal STR loci detected by next generation sequencing.

Sequencing short tandem repeat (STR) loci allows for determination of repeat motif variations within the STR (or entire PCR amplicon) which cannot be ascertained by size-based PCR fragment analysis. Sanger sequencing has been used in research laboratories to further characterize STR loci, but is impractical for routine forensic use due to the laborious nature of the procedure in general and additional steps required to separate heterozygous alleles. Recent advances in library preparation methods enable high-throughput next generation sequencing (NGS) and technological improvements in sequencing chemistries now offer sufficient read lengths to encompass STR alleles. Herein, we present sequencing results from 183 DNA samples, including African American, Caucasian, and Hispanic individuals, at 22 autosomal forensic STR loci using an assay designed for NGS. The resulting dataset has been used to perform population genetic analyses of allelic diversity by length compared to sequence, and exemplifies which loci are likely to achieve the greatest gains in discrimination via sequencing. Within this data set, six loci demonstrate greater than double the number of alleles obtained by sequence compared to the number of alleles obtained by length: D12S391, D2S1338, D21S11, D8S1179, vWA, and D3S1358. As expected, repeat region sequences which had not previously been reported in forensic literature were identified.

STR allele sequence variation: Current knowledge and future issues.

This article reviews what is currently known about short tandem repeat (STR) allelic sequence variation in and around the twenty-four loci most commonly used throughout the world to perform forensic DNA investigations. These STR loci include D1S1656, TPOX, D2S441, D2S1338, D3S1358, FGA, CSF1PO, D5S818, SE33, D6S1043, D7S820, D8S1179, D10S1248, TH01, vWA, D12S391, D13S317, Penta E, D16S539, D18S51, D19S433, D21S11, Penta D, and D22S1045. All known reported variant alleles are compiled along with genomic information available from GenBank, dbSNP, and the 1000 Genomes Project. Supplementary files are included which provide annotated reference sequences for each STR locus, characterize genomic variation around the STR repeat region, and compare alleles present in currently available STR kit allelic ladders. Looking to the future, STR allele nomenclature options are discussed as they relate to next generation sequencing efforts underway.

Performance of a next generation sequencing SNP assay on degraded DNA.

Forensic DNA casework samples are often of insufficient quantity or quality to generate full profiles by conventional DNA typing methods. Polymerase chain reaction (PCR) amplification of short tandem repeat (STR) loci is inherently limited in samples containing degraded DNA, as the cumulative size of repeat regions, primer binding regions, and flanking sequence is necessarily larger than the PCR template. Additionally, traditional capillary electrophoresis (CE) assay design further inherently limits shortening amplicons because the markers must be separated by size. Non-traditional markers, such as single nucleotide polymorphisms (SNPs) and insertion deletion polymorphisms (InDels), may yield more information from challenging samples due to their smaller amplicon size. In this study, the performance of a next generation sequencing (NGS) SNP assay and CE-based STR, mini-STR, and InDel assays was evaluated with a series of fragmented, size-selected samples. Information obtained from the NGS SNP assay exhibited higher overall inverse random match probability (1/RMP) values compared to the CE-based typing assays, with particular benefit for fragment sizes ≤ 150 base pairs (bp). The InDel, mini-STR, and NGS SNP assays all had similar percentages of loci with reportable alleles at this level of degradation; however, the relatively fewer number of loci in the InDel and mini-STR assays results in the NGS SNP assay having at least nine orders of magnitude higher 1/RMP values. In addition, the NGS SNP assay and three CE-based assays (two STR and one InDel assay) were tested using a dilution series consisting of 0.5 ng, 0.1 ng, and 0.05 ng non-degraded DNA. All tested assays showed similar percentages of loci with reportable alleles at these levels of input DNA; however, due to the larger number of loci, the NGS SNP assay and the larger of the two tested CE-based STR assays both resulted in considerably higher 1/RMP values than the other assays. These results indicate the potential advantage of NGS SNP assays for forensic analysis of degraded DNA samples.

A 50-SNP assay for biogeographic ancestry and phenotype prediction in the U.S. population.

When an STR DNA profile obtained from crime scene evidence does not match identified suspects or profiles from available databases, further DNA analyses targeted at inferring the possible ancestral origin and phenotypic characteristics of the perpetrator could yield valuable information. Single Nucleotide Polymorphisms (SNPs), the most common form of genetic polymorphisms, have alleles associated with specific populations and/or correlated to physical characteristics. We have used single base primer extension (SBE) technology to develop a 50 SNP assay (composed of three multiplexes) designed to predict ancestry among the primary U.S. populations (African American, East Asian, European American, and Hispanic American/Native American), as well as pigmentation phenotype (eye, hair, and skin color) among European American. We have optimized this assay to a sensitivity level comparable to current forensic DNA analyses, and shown robust performance on forensic-type samples. In addition, we developed a prediction model for ancestry in the U.S. population, based on the random match probability and likelihood ratio formulas already used in forensic laboratories. Lastly, we evaluated the biogeographic ancestry prediction model using a test set, and we evaluated an existing model for eye color with our U.S. sample set. Using these models with recommended thresholds, the 50 SNP assay provided accurate ancestry information in 98.6% of the test set samples, and provided accurate eye color information in 61% of the European samples tested (25% were inconclusive and 14% were incorrect). This method, which uses equipment already available in forensic DNA laboratories, is recommended for use in U.S. forensic casework to provide additional information about the donor of a DNA sample when the STR profile has not been linked to an individual.