Autosomal SNP typing of forensic samples with the GenPlex™ HID System: results of a collaborative study

TIA: algorithms for development of identity-linked SNP islands for analysis by massively parallel DNA sequencing

Background

Single nucleotide polymorphisms (SNPs) located within the human genome have been shown to have utility as markers of identity in the differentiation of DNA from individual contributors. Massively parallel DNA sequencing (MPS) technologies and human genome SNP databases allow for the design of suites of identity-linked target regions, amenable to sequencing in a multiplexed and massively parallel manner. Therefore, tools are needed for leveraging the genotypic information found within SNP databases for the discovery of genomic targets that can be evaluated on MPS platforms.

Results

The SNP island target identification algorithm (TIA) was developed as a user-tunable system to leverage SNP information within databases. Using data within the 1000 Genomes Project SNP database, human genome regions were identified that contain globally ubiquitous identity-linked SNPs and that were responsive to targeted resequencing on MPS platforms. Algorithmic filters were used to exclude target regions that did not conform to user-tunable SNP island target characteristics. To validate the accuracy of TIA for discovering these identity-linked SNP islands within the human genome, SNP island target regions were amplified from 70 contributor genomic DNA samples using the polymerase chain reaction. Multiplexed amplicons were sequenced using the Illumina MiSeq platform, and the resulting sequences were analyzed for SNP variations. 166 putative identity-linked SNPs were targeted in the identified genomic regions. Of the 309 SNPs that provided discerning power across individual SNP profiles, 74 previously undefined SNPs were identified during evaluation of targets from individual genomes. Overall, DNA samples of 70 individuals were uniquely identified using a subset of the suite of identity-linked SNP islands.

Conclusions

TIA offers a tunable genome search tool for the discovery of targeted genomic regions that are scalable in the population frequency and numbers of SNPs contained within the SNP island regions. It also allows the definition of sequence length and sequence variability of the target region as well as the less variable flanking regions for tailoring to MPS platforms. As shown in this study, TIA can be used to discover identity-linked SNP islands within the human genome, useful for differentiating individuals by targeted resequencing on MPS technologies.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2133-2) contains supplementary material, which is available to authorized users.

Development of a SNP-based panel for human identification for Indian populations

The widely employed short tandem repeat (STR)-based panels for forensic human identification (HID) have limitations while dealing with challenging forensic samples involving DNA degradation, resulting in dropping-out of higher molecular weight alleles/loci. To address this issue, bialleic markers like single nucleotide polymorphisms (SNPs) and insertion-deletions (indels), which can be scored even when the template DNA is heavily degraded (<100bp), have been suggested as alternative markers for HID testing. Recent studies have highlighted their utility in forensic HID and several panels based on biallelic markers have been described for worldwide populations. However, there has been very little information about the behavior of such DNA markers in Indian populations, which is known to possess great genetic diversity. This study describes a two-step approach for designing a SNP-based panel consisting of 70 SNPs for HID testing in Indian populations. In the first step, candidate SNPs were shortlisted from public databases by screening them for several criteria including allelic distribution, genomic location, potential phenotypic expression or functionality and species specificity. The second step involved genotyping the shortlisted SNPs in various Indian populations followed by shortlisting of the best performers for identity-testing. Starting with 592,652 SNPs listed in Human660W-Quad Beadchip (Illumina Inc.), we shortlisted 275 candidate SNPs for identity-testing and genotyped them in 462 unrelated individuals from different population groups in India. Post genotyping and statistical analyses based on biogeographic regions, 206 SNPs demonstrated desired allelic distribution (Heterozygosity≥0.4 and F_ST≤0.02), from which 2-4 widely separated (>20 Mb apart) SNPs from each chromosome were finally selected to construct a panel of 70 SNPs. This panel on average possessed match probability 10e-29 and probability of paternity of 0.99999997, which was orders of magnitude higher than most of the currently employed STR-based chemistries and SNP-based panels that were proposed previously for HID testing. For comparison purpose, genotyping previously reported SNPs for HID in our samples led us to conclude that the panel developed in this study is much more efficient and robust and better suited for the Indian populations.

Massively parallel sequencing of forensically relevant single nucleotide polymorphisms using TruSeq™ forensic amplicon

The TruSeq™ Forensic Amplicon library preparation protocol, originally designed to attach sequencing adapters to chromatin-bound DNA for chromatin immunoprecipitation sequencing (TruSeq™ ChIP-Seq), was used here to attach adapters directly to amplicons containing markers of forensic interest. In this study, the TruSeq™ Forensic Amplicon library preparation protocol was used to detect 160 single nucleotide polymorphisms (SNPs), including human identification SNPs (iSNPs), ancestry, and phenotypic SNPs (apSNPs) in 12 reference samples. Results were compared with those generated by a second laboratory using the same technique, as well as to those generated by whole genome sequencing (WGS). The genotype calls made using the TruSeq™ Forensic Amplicon library preparation protocol were highly concordant. The protocol described herein represents an effective and relatively sensitive means of preparing amplified nuclear DNA for massively parallel sequencing (MPS).

A 27-locus STR assay to meet all United States and European law enforcement agency standards

Different national and international agencies have selected specific STR sets for forensic database use. To enhance database comparison across national and international borders, a 27-locus multiplex system was developed comprising all 15 STR loci of the European standard set, the current 13 STR loci of the CODIS core, the proposed 22 STR loci of the expanded CODIS core, 4 additional commonly used STR loci, and the amelogenin locus. Development required iterative primer design to resolve primer-related artifacts, amplicon sizing, and locus-to-locus balance issues. The 19.5-min assay incorporated newly developed six-dye chemistry analyzed using a novel microfluidic electrophoresis instrument capable of simultaneous detection and discrimination of 8 or more fluorescent dyes. The 27-locus multiplex offers the potential for a new international STR standard permitting laboratories in any jurisdiction to use a single reaction to determine profiles for loci they typically generate plus an expanded common STR profiling set of global interest.

Evaluation of the iPLEX® Sample ID Plus Panel designed for the Sequenom MassARRAY® system. A SNP typing assay developed for human identification and sample tracking based on the SNPforID panel

Sequenom launched the first commercial SNP typing kit for human identification, named the iPLEX(®) Sample ID Plus Panel. The kit amplifies 47 of the 52 SNPs in the SNPforID panel, amelogenin and two Y-chromosome SNPs in one multiplex PCR. The SNPs were analyzed by single base extension (SBE) and Matrix Assisted Laser Desorption/Ionization-Time of Flight Mass Spectrometry (MALDI-TOF MS). In this study, we evaluated the accuracy and sensitivity of the iPLEX(®) Sample ID Plus Panel by comparing the typing results of the iPLEX(®) Sample ID Plus Panel with those obtained with our ISO 17025 accredited SNPforID assay. The average call rate for duplicate typing of any one SNPs in the panel was 90.0% when the mass spectra were analyzed automatically with the MassARRAY(®) TYPER 4.0 genotyping software in real time. Two reproducible inconsistencies were observed (error rate: 0.05%) at two different SNP loci. In addition, four inconsistencies were observed once. The optimal amount of template DNA in the PCR was ≥10ng. There was a relatively high risk of allele and locus drop-outs when ≤1ng template DNA was used. We developed an R script with a stringent set of "forensic analysis parameters" based on the peak height and the signal to noise data exported from the TYPER 4.0 software. With the forensic analysis parameters, all inconsistencies were eliminated in reactions with ≥10ng DNA. However, the average call rate decreased to 69.9%. The iPLEX(®) Sample ID Plus Panel was tested on 10 degraded samples from forensic case-work. Two samples could not be typed, presumably because the samples contained PCR and SBE inhibitors. The average call rate was generally lower for degraded DNA samples and the number of inconsistencies higher than for pristine DNA. However, none of the inconsistencies were reproduced and the highest match probability for the degraded samples typed with the panel was 1.7E-9 using the stringent forensic analysis parameters. Although the relatively low sensitivity of the iPLEX(®) Sample ID Plus Panel makes it inappropriate for typing of trace samples from crime scenes, the panel may be interesting for relationship testing and for identification of e.g. samples in biobanks because of the low reagent costs, the limited hands-on time of the iPLEX(®) assay and the automatic analysis of the mass spectra.

Revision of the SNPforID 34-plex forensic ancestry test: Assay enhancements, standard reference sample genotypes and extended population studies

A revision of an established 34 SNP forensic ancestry test has been made by swapping the under-performing rs727811 component SNP with the highly informative rs3827760 that shows a near-fixed East Asian specific allele. We collated SNP variability data for the revised SNP set in 66 reference populations from 1000 Genomes and HGDP-CEPH panels and used this as reference data to analyse four U.S. populations showing a range of admixture patterns. The U.S. Hispanics sample in particular displayed heterogeneous values of co-ancestry between European, Native American and African contributors, likely to reflect in part, the way this disparate group is defined using cultural as well as population genetic parameters. The genotyping of over 700 U.S. population samples also provided the opportunity to thoroughly gauge peak mobility variation and peak height ratios observed from routine use of the single base extension chemistry of the 34-plex test. Finally, the genotyping of the widely used DNA profiling Standard Reference Material samples plus other control DNAs completes the audit of the 34-plex assay to allow forensic practitioners to apply this test more readily in their own laboratories.

Kinship Analysis with Diallelic SNPs - Experiences with the SNPforID Multiplex in an ISO17025 Accreditated Laboratory

BACKGROUND: The mutation rate of single nucleotide polymorphisms (SNPs) is estimated to be 100,000 times lower than that of short tandem repeats (STRs), which makes SNPs very suitable for relationship testing. The SNPforID multiplex assay was the first SNP typing assay that was a real alternative to the commonly used STR kits in kinship and crime case work and the first SNP assay to be validated in a forensic laboratory accredited according to the ISO17025 standard. METHODS: A total of 54 crime case samples were typed with the SNPforID multiplex assay. 30 samples from relationship cases were sequenced in selected SNP loci. RESULTS: It was demonstrated that mixtures were easily detected with the SNPforID assay by analyzing the signal strengths of the detected alleles. Unusual imbalances in signal strengths that were observed in a few individuals could be explained by unexpected SNPs in one of the primer binding sites. A complicated relationship case with four closely related individuals is presented. CONCLUSION: Mixtures can be detected with bi-allelic SNPs. The SNPforID assay is a very useful supplement to the STR kits in relationship testing.

A 48-plex autosomal SNP GenPlex™ assay for human individualization and relationship testing

SNPs are being increasingly used by forensic laboratories. Different platforms have been developed for SNP typing. We describe the GenPlex™ HID system protocol, a new SNP-typing platform developed by Applied Biosystems where 48 of the 52 SNPforID SNPs and amelogenin are included. The GenPlex™ HID system protocol has been successfully tested by a number of forensic laboratories using both ordinary and forensic samples.