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.

Massively parallel sequencing holds great promise for new possibilities in the field of forensic genetics, enabling simultaneous analysis of multiple markers as well as offering enhanced short tandem repeat allele resolution. A challenge in forensic DNA analysis is that the samples often contain low amounts of DNA in a background that may interfere with downstream analysis. PCR inhibition mechanisms of some relevant molecules have been studied applying e.g. real-time PCR and digital PCR. However, a detailed understanding of the effects of inhibitory molecules on forensic MPS, including mechanisms and ways to relieve inhibition, is missing. In this study, the effects of two well-characterized PCR inhibitors, humic acid and hematin, have been studied using the ForenSeq DNA Signature Prep kit. Humic acid and hematin resulted in lowered read numbers as well as specific negative effects on certain markers. Quality control of libraries with Fragment analyzer showed that increasing amounts of inhibitors caused a lowered amplicon quantity and that the larger amplicons were more likely to drop out. Further, the inhibitor tolerance could be improved 5-10 times by addition of bovine serum albumin in the initial PCR. On the contrary to the samples with inhibitors, low-template samples resulted in lowered read numbers for all markers. This difference strengthened the conclusion that the inhibitors have a negative effect on the DNA polymerase activity in the initial PCR. Additionally, a common capillary gel electrophoresis-based STR kit was shown to handle at least 200 times more inhibitors than the ForenSeq DNA Signature Prep kit. This suggests that there is room for improvement of the PCR components to ensure analytical success for challenging samples, which is needed for a broad application of MPS for forensic STR analysis.

A total of 1327 platinum-quality mitochondrial DNA haplotypes from United States (U.S.) populations were generated using a robust, semi-automated next-generation sequencing (NGS) workflow with rigorous quality control (QC). The laboratory workflow involved long-range PCR to minimize the co-amplification of nuclear mitochondrial DNA segments (NUMTs), PCR-free library preparation to reduce amplification bias, and high-coverage Illumina MiSeq sequencing to produce an average per-sample read depth of 1000 × for low-frequency (5%) variant detection. Point heteroplasmies below 10% frequency were confirmed through replicate amplification, and length heteroplasmy was quantitatively assessed using a custom read count analysis tool. Data analysis involved a redundant, dual-analyst review to minimize errors in haplotype reporting with additional QC checks performed by EMPOP. Applying these methods, eight sample sets were processed from five U.S. metapopulations (African American, Caucasian, Hispanic, Asian American, and Native American) corresponding to self-reported identity at the time of sample collection. Population analyses (e.g., haplotype frequencies, random match probabilities, and genetic distance estimates) were performed to evaluate the eight datasets, with over 95% of haplotypes unique per dataset. The platinum-quality mitogenome haplotypes presented in this study will enable forensic statistical calculations and thereby support the usage of mitogenome sequencing in forensic laboratories.