Non-uniform phenotyping of D12S391 resolved by second generation sequencing

Developmental validation of the ForenSeq MainstAY kit, MiSeq FGx sequencing system and ForenSeq Universal Analysis Software

For human identification purposes, forensic genetics has primarily relied upon a core set of autosomal (and to a lesser extent Y chromosome) short tandem repeat (STR) markers that are enriched by amplification using the polymerase chain reaction (PCR) that are subsequently separated and detected using capillary electrophoresis (CE). While STR typing conducted in this manner is well-developed and robust, advances in molecular biology that have occurred over the last 15 years, in particular massively parallel sequencing (MPS) [1-7], offer certain advantages as compared to CE-based typing. First and foremost is the high throughput capacity of MPS. Current bench top high throughput sequencers enable larger batteries of markers to be multiplexed and multiple samples to be sequenced simultaneously (e.g., millions to billions of nucleotides can be sequenced in one run). Second, compared to the length-based CE approach, sequencing STRs increases discrimination power, enhances sensitivity of detection, reduces noise due to instrumentation, and improves mixture interpretation [4,8-23]. Third, since detection of STRs is based on sequence and not fluorescence, amplicons can be designed that are shorter in length and of similar lengths among loci, where possible, which can improve amplification efficiency and analysis of degraded samples. Lastly, MPS offers a single format approach that can be applied to analysis of a wide variety of genetic markers of forensic interest (e.g., STRs, mitochondrial DNA, single nucleotide polymorphisms, insertion/deletions). These features make MPS a desirable technology for casework [14,15,24,25-48]. The developmental validation of the ForenSeq MainstAY library preparation kit with the MiSeq FGx Sequencing System and ForenSeq Universal Software is reported here to assist with validation of this MPS system for casework [49]. The results show that the system is sensitive, accurate and precise, specific, and performs well with mixtures and mock case-type samples.

A new digital approach to SNP encoding for DNA identification

Identification of individuals has become an urgent problem for mankind. In the last three decades, STR-based DNA identification has actively evolved along with traditional biometric methods. Nonetheless, single-nucleotide polymorphisms (SNPs) are now of great interest and a number of relevant SNP panels have been proposed for DNA identification. Here, a simple approach to SNP data digitization that can provide assigning a unique genetic identification number (GIN) to each person is proposed. The key points of this approach are as follows: 1) SNP data are digitized as whole 4-bit boxes in the most convenient binary format, where character "1" (YES) is assigned to revealed nucleotides, and character "0" (NO) to missing nucleotides after SNP-typing; 2) all SNPs should be considered tetra-allelic. Calculations showed that a 72-plex SNP panel is enough to provide the population with unique GINs, which can be represented in digital (binary or hexadecimal) or graphic (linear or two-dimensional) formats. Simple software for SNP data processing and GINs creation in any format was written. It is likely that the national and global GIN databases will facilitate the solution of problems related to identification of individuals or human biological materials. The proposed approach may be extended to other living organisms as well.

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.

Massively Parallel Sequencing of Forensic STRs Using the Ion Chef™ and the Ion S5™ XL Systems

Next-generation sequencing (NGS) has been used to genotype forensic short tandem repeat (STR) markers for individual identification and kinship analysis. STR data from several NGS platforms have been published, but forensic application trials using the Ion S5™ XL system have not been reported. In this work, we report sensitivity, reproducibility, mixture, simulated degradation, and casework sample data on the Ion Chef™ and S5™ XL systems using an early access 25-plex panel. Sensitivity experiments showed that over 97% of the alleles were detectable with down to 62 pg input of genomic DNA. In mixture studies, alleles from minor contributors were correctly assigned at 1:9 and 9:1 ratios. NGS successfully gave 12 full genotype results from 13 challenging casework samples, compared with five full results using the CE platform. In conclusion, the Ion Chef™ and the Ion S5™ XL systems provided an alternative and promising approach for forensic STR genotyping.

Next-generation sequencing analysis of off-ladder alleles due to migration shift caused by sequence variation at D12S391 locus

In short tandem repeat (STR) analysis, length polymorphisms are detected by capillary electrophoresis (CE). At most STR loci, mobility shift due to sequence variation in the repeat region was thought not to affect the typing results. In our recent population studies of 1501 Japanese individuals, off-ladder calls were observed at the D12S391 locus using PowerPlex Fusion in nine samples for allele 22, one sample for allele 25, and one sample for allele 26. However, these samples were typed as ordinary alleles within the bins using GlobalFiler. In this study, next-generation sequencing analysis using MiSeq was performed for the D12S391 locus from the 11 off-ladder samples and 33 other samples, as well as the allelic ladders of PowerPlex Fusion and GlobalFiler. All off-ladder allele 22 in the nine samples had [AGAT]11[AGAC]11 as a repeat structure, while the corresponding allele was [AGAT]15[AGAC]6[AGAT] for the PowerPlex Fusion ladder, and [AGAT]13[AGAC]9 for the GlobalFiler ladder. Overall, as the number of [AGAT] in the repeat structure decreased at the D12S391 locus, the peak migrated more slowly using PowerPlex Fusion, the reverse strand of which was labeled, and it migrated more rapidly using GlobalFiler, the forward strand of which was labeled. The allelic ladders of both STR kits were reamplified with our small amplicon D12S391 primers and their mobility was also examined. In conclusion, off-ladder observations of allele 22 at the D12S391 locus using PowerPlex Fusion were mainly attributed to a relatively large difference of the repeat structure between its allelic ladder and off-ladder allele 22.

Exploring the efficacy of paternity and kinship testing based on single nucleotide polymorphisms

Short tandem repeats (STRs) are conventional genetic markers typically used for paternity and kinship testing. As supplementary markers of STRs, single nucleotide polymorphisms (SNPs) have less discrimination power but broader applicability to degraded samples. The rapid improvement of next-generation sequencing (NGS) and multiplex amplification technologies also make it possible now to simultaneously identify dozens or even hundreds of SNP loci in a single pool. However, few studies have been endeavored to kinship testing based on SNP loci. In this study, we genotyped 90 autosomal human identity SNP loci with NGS, and investigated their testing efficacies based on the likelihood ratio model in eight pedigree scenarios involving paternity, half/full-sibling, uncle/nephew, and first-cousin relationships. We found that these SNPs might be sufficient to discriminate paternity and full-sibling, but impractical for more distant relatives such as uncle and cousin. Furthermore, we conducted an in silico study to obtain the theoretical tendency of how testing efficacy varied with increasing number of SNP loci. For each testing battery in a given pedigree scenario, we obtained distributions of logarithmic likelihood ratio for both simulated relatives and unrelated controls. The proportion of the overlapping area between the two distributions was defined as a false testing level (FTL) to evaluate the testing efficacy. We estimated that 85, 127, 491, and 1,858 putative SNP loci were required to discriminate paternity, full-sibling, half-sibling/uncle-nephew, and first-cousin (FTL, 0.1%), respectively. To test a half-sibling or nephew, an additional uncle relative could be included to decrease the required number of putative SNP loci to ∼320 (FTL, 0.1%). As a systematic computation of paternity and kinship testing based only on SNPs, our results could be informative for further studies and applications on paternity and kinship testing using SNP loci.

Introduction of the Python script STRinNGS for analysis of STR regions in FASTQ or BAM files and expansion of the Danish STR sequence database to 11 STRs

This work introduces the in-house developed Python application STRinNGS for analysis of STR sequence elements in BAM or FASTQ files. STRinNGS identifies sequence reads with STR loci by their flanking sequences, it analyses the STR sequence and the flanking regions, and generates a report with the assigned SNP-STR alleles. The main output file from STRinNGS contains all sequences with read counts above 1% of the total number of reads per locus. STR sequences are automatically named according to the nomenclature used previously and according to the repeat unit definitions in STRBase (http://www.cstl.nist.gov/strbase/). The sequences are named with (1) the locus name, (2) the length of the repeat region divided by the length of the repeat unit, (3) the sequence(s) of the repeat unit(s) followed by the number of repeats and (4) variations in the flanking regions. Lower case letters in the main output file are used to flag sequences with previously unknown variations in the STRs. SNPs in the flanking regions are named by their "rs" numbers and the nucleotides in the SNP position. Data from 207 Danes sequenced with the Ion Torrent™ HID STR 10-plex that amplified nine STRs (CSF1PO, D3S1358, D5S818, D7S820, D8S1179, D16S539, TH01, TPOX, vWA), and Amelogenin was analysed with STRinNGS. Sequencing uncovered five common SNPs near four STRs and revealed 20 new alleles in the 207 Danes. Three short homopolymers in the D8S1179 flanking regions caused frequent sequencing errors. In 29 of 3726 allele calls (0.8%), sequences with homopolymer errors were falsely assigned as true alleles. An in-house developed script in R compensated for these errors by compiling sequence reads that had identical STR sequences and identical nucleotides in the five common SNPs. In the output file from the R script, all SNP-STR haplotype calls were correct. The 207 samples and six additional samples were sequenced for D3S1358, D12S391, and D21S11 using the 454 GS Junior platform in this and a previous work. Overall, next generation sequencing (NGS) of the 11 STRs lowered the mean match probability 386 times and increased the typical paternity indexes (i.e. the geometric mean) for trios and duos 47 and 23 times, respectively, compared to the traditional PCR-CE typing of the same population.

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.

A novel set of DIP-STR markers for improved analysis of challenging DNA mixtures

Currently available molecular biology tools allow forensic scientists to characterize DNA evidence found at crime scenes for a large variety of samples, including those of limited quantity and quality, and achieve high levels of individualization. Yet, standard forensic markers provide limited or no results when applied to mixed DNA samples where the contributors are present in very different proportions (unbalanced DNA mixtures). This becomes an issue mostly for the analysis of trace samples collected on the victim or from touched objects. To this end, we recently proposed an innovative type of genetic marker, named DIP-STR that relies on pairing deletion/insertion polymorphisms (DIP) with standard short tandem repeats (STR). This novel compound marker allows detection of the minor DNA contributor in a DNA mixture of any gender and cellular origin with unprecedented resolution (beyond a DNA ratio of 1:1000). To provide a novel analytical tool useful in practice to common forensic laboratories, this article describes the first set of 10 DIP-STR markers selected according to forensic technical standards. The novel DIP-STR regions are short (between 146 and 271 bp), include only highly polymorphic tri-, tetra- and pentanucleotide tandem repeats and are located on different chromosomes or chromosomal arms to provide statistically independent results. This novel set of DIP-STR can target the amplification of 0.03-0.1 ng of DNA when mixed with a 1000-fold excess of major DNA. DIP-STR relative allele frequencies are estimated based on a survey of 103 Swiss individuals. Finally, this study provides an estimate of the occurrence of informative alleles and a calculation of the corresponding random match probability of the detected minor DIP-STR genotype assessed across 10,506 pairwise conceptual mixtures.

Massively parallel sequencing of complete mitochondrial genomes from hair shaft samples

Though shed hairs are one of the most commonly encountered evidence types, they are among the most limited in terms of DNA quantity and quality. As a result, DNA testing has historically focused on the recovery of just about 600 base pairs of the mitochondrial DNA control region. Here, we describe our success in recovering complete mitochondrial genome (mtGenome) data (∼16,569bp) from single shed hairs. By employing massively parallel sequencing (MPS), we demonstrate that particular hair samples yield DNA sufficient in quantity and quality to produce 2-3kb mtGenome amplicons and that entire mtGenome data can be recovered from hair extracts even without PCR enrichment. Most importantly, we describe a small amplicon multiplex assay comprised of sixty-two primer sets that can be routinely applied to the compromised hair samples typically encountered in forensic casework. In all samples tested here, the MPS data recovered using any one of the three methods were consistent with the control Sanger sequence data developed from high quality known specimens. Given the recently demonstrated value of complete mtGenome data in terms of discrimination power among randomly sampled individuals, the possibility of recovering mtGenome data from the most compromised and limited evidentiary material is likely to vastly increase the utility of mtDNA testing for hair evidence.

Second generation sequencing of three STRs D3S1358, D12S391 and D21S11 in Danes and a new nomenclature for sequenced STR alleles

Second generation sequencing (SGS) may revolutionize the field of forensic STR typing. Two of the essential requirements for implementation of an SGS based approach for forensic investigations are (1) establishment of adequate frequency databases and (2) adoption of a new STR nomenclature. We report the STR sequences and allele frequencies of three STR loci: D3S1358, D12S391 and D21S11 in 197 unrelated Danes. We used a new STR nomenclature that depicts the locus name used in forensic genetics, the length of the repeat region divided by the repeat length (typically 4 nucleotides) and detailed sequence information of possible sub-repeats and SNPs within the amplified fragment.