A preliminary assessment of the ForenSeq™ FGx System: next generation sequencing of an STR and SNP multiplex

Population genetic analyses of Eastern Chinese Han nationality using ForenSeq™ DNA Signature Prep Kit

Currently, the most commonly used method for human identification and kinship analysis in forensic genetics is the detection of length polymorphism in short tandem repeats (STRs) using polymerase chain reaction (PCR) and capillary electrophoresis (CE). However, numerous studies have shown that considerable sequence variations exist in the repeat and flanking regions of the STR loci, which cannot be identified by CE detection. Comparatively, massively parallel sequencing (MPS) technology can capture these sequence differences, thereby enhancing the identification capability of certain STRs. In this study, we used the ForenSeq™ DNA Signature Prep Kit to sequence 58 STRs and 94 individual identification SNPs (iiSNPs) in a sample of 220 unrelated individuals from the Eastern Chinese Han population. Our aim is to obtain MPS-based STR and SNP data, providing further evidence for the study of population genetics and forensic applications. The results showed that the MPS method, utilizing sequence information, identified a total of 486 alleles on autosomal STRs (A-STRs), 97 alleles on X-chromosome STRs (X-STRs), and 218 alleles on Y-chromosome STRs (Y-STRs). Compared with length polymorphism, we observed an increase of 260 alleles (157, 31, and 72 alleles on A-STRs, X-STRs, and Y-STRs, respectively) across 36 STRs. The most substantial increments were observed in DYF387S1 and DYS389II, with increases of 287.5% and 250%, respectively. The most increment in the number of alleles was found at DYF387S1 and DYS389II (287.5% and 250%, respectively). The length-based (LB) and sequence-based (SB) combined random match probability (RMP) of 27 A-STRs were 6.05E-31 and 1.53E-34, respectively. Furthermore, other forensic parameters such as total discrimination power (TDP), cumulative probability of exclusion of trios (CPEtrio), and duos (CPEduo) were significantly improved when using the SB data, and informative data were obtained for the 94 iiSNPs. Collectively, these findings highlight the advantages of MPS technology in forensic genetics, and the Eastern Chinese Han genetic data generated in this study could be used as a valuable reference for future research in this field.

Forensic analysis and sequence variation of 133 STRs in the Hakka population

Introduction: Short Tandem Repeats (STRs) are highly valuable genetic markers in forensic science. However, the conventional PCR-CE technique has limitations, and the emergence of massively parallel sequencing (MPS) technology presents new opportunities for STR analysis. Yet, there is limited research on Chinese population diversity using MPS. Methods: In this study, we obtained genotype data for 52 A-STRs and 81 Y-STRs from the Hakka population in Meizhou, Guangdong, China, using the Forensic Analysis System Multiplecues SetB Kit on the MGISEQ-2000 platform. Results: Our findings demonstrate that these 133 STRs are highly efficient for forensic applications within the Meizhou Hakka population. Statistical analysis revealed Hobs values ranging from 0.61306 to 0.91083 and Hexp values ranging from 0.59156 to 0.91497 for A-STRs based on length polymorphism. For sequence polymorphism, Hobs values ranged from 0.61306 to 0.94586, and Hexp values fluctuated between 0.59156 and 0.94487. The CPE values were 1-5.0779620E-21 and 1-3.257436E-24 for length and sequence polymorphism, respectively, while the CPD values were 1-1.727007E-59 and 1-5.517015E-66, respectively. Among the 80 Y-STR loci, the HD values for length and sequence polymorphism were 0.99764282 and 0.99894195, respectively. The HMP values stood at 0.00418102 and 0.00288427, respectively, and the DC values were 0.75502742 and 0.83363803, respectively. For the 52 A-STR loci, we identified 554 and 989 distinct alleles based on length and sequence polymorphisms, respectively. For the 81 Y-STR loci, 464 and 652 unique alleles were detected at the length and sequence level, respectively. Population genetic analysis revealed that the Meizhou Hakka population has a close kinship relationship with the Asian populations THI and KOR based on length polymorphism data of A-STRs. Conversely, based on length polymorphism data of Y-STRs, the Meizhou Hakka population has the closest kinship relationship with the Henan Han population. Discussion: Overall, the variation information of repeat region sequences significantly enhances the forensic identification efficacy of STR genetic markers, providing an essential database for forensic individual and paternity testing in this region. Additionally, the data generated by our study will serve as a vital resource for research into the genetic structure and historical origins of the Meizhou Hakka population.

High-resolution genotyping of 58 STRs in 635 Northern Han Chinese with MiSeq FGx ® Forensic Genomics System

Sequence polymorphisms were characterized at 27 autosomal STRs (A-STRs), 7 X chromosomal STRs (X-STRs), and 24 Y chromosomal STRs (Y-STRs) in 635 Northern Han Chinese with the ForenSeq DNA Signature Prep Kit on the MiSeq FGx Forensic Genomics System. Since repeat region (RR) and flanking region (FR) variation can be detected by massively parallel sequencing (MPS), the increase in the number of unique alleles and the average of gene diversity was 78.18% and 3.51% between sequence and length, respectively. A total of 74 novel RR variants were identified at 33 STRs compared with STRSeq and previous studies, and 13 FR variants (rs1770275883, rs2053373277, rs2082557941, rs1925525766, rs1926380862, rs1569322793, rs2051848492, rs2051848696, rs2016239814, rs2053269960, rs2044518192, rs2044536444, and rs2089968964) were first submitted to dbSNP. Also, 99.94% of alleles were concordant between the ForenSeq DNA Signature Prep Kit and commercial CE kits. Discordance resulted from the low performance at D22S1045 and occasionally at DYS392, flanking region deletions at D7S820 and DXS10074, and the strict alignment algorithm at DXS7132. Null alleles at DYS505 and DYS448 and multialleles at DYS387S1a/b, DYS385a/b, DYS448, DYS505, DXS7132, and HPRTB were validated with other MPS and CE kits. Thus, a high-resolution sequence-based (SB) and length-based (LB) allele frequencies dataset from Northern Han Chinese has been established already. As expected, forensic parameters increased significantly on combined power of discrimination (PD) and combined power of exclusion (PE) at A-STRs, mildly on combined PD and combined mean exclusion chance (MEC) at X-STRs, and barely on discrimination capacity (DC) at Y-STRs. Additionally, MiSeq FGx quality metrics and MPS performance were evaluated in this study, which presented the high-quality of the dataset at 20 consecutive runs, such as ≥ 60% bases with a quality score of 20 or higher (%≥ Q20), > 60% of effective reads, > 2000 × of depth of coverage (DoC), ≥ 60% of allele coverage ratio (ACR) or heterozygote balance, ≥ 70% of inter-locus balance, and ≤ 0.4 of the absolute value of observed minus expected heterozygosity (|H_exp - H_obs|). In conclusion, MiSeq FGx can help us generate a high-resolution and high-quality dataset for human identification and population genetic studies.

From crime scene to courtroom: A review of the current bioanalytical evidence workflows used in rape and sexual assault investigations in the United Kingdom

Sexual assault casework requires the collaboration of multiple agency staff to formalise an investigative pipeline running from crime scene to court. While the same could be said of many other forensic investigations, few require the additional support of health care staff and the combined forensic involvement of body-fluid examiners, DNA experts and analytical chemists. The sheer amount of collaborative effort between agencies is laid out through a detailed examination of the investigative workflow from crime scene to courtroom with each step in the pipelines detailed and discussed. Beginning with a review of sexual assault legislation in the United Kingdom this article details how sexual assault investigations are initiated by police and supported by sexual assault referral centre (SARC) staff who are often the first responders providing primary healthcare and patient support to victims while simultaneously collecting and assessing forensic evidence. Detailing the myriad of evidential material that can be documented and collected at the SARC, the review identifies and categorises key forensic tests to first detect and identify body-fluids recovered from evidence through to the secondary analysis of DNA to help identify the suspect. This review also focusses on the collection and analysis of biological material used to support the allegation that the sexual activity was non-consensual and provides a breakdown of common marks and trauma as well as a review of common analytical methods used to infer Drug Facilitated Sexual Assault (DFSA). The culmination of the investigative pipeline is discussed by reviewing the Rape and Serious Sexual Assault (RASSO) workflow used by the Crown Prosecution Service before providing our thoughts on the future of forensic analysis and possible changes to the described workflows.

Validation and beyond: Next generation sequencing of forensic casework samples including challenging tissue samples from altered human corpses using the MiSeq FGx system

The proceeding developments in next generation sequencing (NGS) technologies enable increasing discrimination power for short tandem repeat (STR) analyses and provide new possibilities for human identification. Therefore, the growing relevance and demand in forensic casework display the need for reliable validation studies and experiences with challenging DNA samples. The presented validation of the MiSeq FGx system and the ForenSeq™ DNA Signature Prep Kit (1) investigated sensitivity, repeatability, reproducibility, concordance, pooling variations, DNA extraction method variances, DNA mixtures, degraded, and casework samples and (2) optimized the sequencing workflow for challenging samples from human corpses by testing additional PCR purification, pooling adjustments, and adapter volume reductions. Overall results indicate the system's reliability in concordance to traditional capillary electrophoresis (CE)‐based genotyping and reproducibility of sequencing data. Genotyping success rates of 100% were obtained down to 62.5 pg DNA input concentrations. Autosomal STR (aSTR) profiles of artificially degraded samples revealed significantly lower numbers of locus and allelic dropouts than CE. However, it was observed that the system still exposed drawbacks when sequencing highly degraded and inhibited samples from human remains. Due to the lack of studies evaluating the sequencing success of samples from decomposed or skeletonised corpses, the presented optimisation studies provide valuable recommendations such as an additional PCR purification, an increase in library pooling volumes, and a reduction of adapter volumes for samples with concentrations ≥31.2 pg. Thus, this research highlights the importance of all‐encompassing validation studies for implementing novel technologies in forensic casework and presents recommendations for challenging samples.

Predicting eye and hair colour in a Norwegian population using Verogen's ForenSeq™ DNA signature prep kit

Prediction of eye and hair colour from DNA can be an important investigative tool in forensic cases if conventional DNA profiling fails to match DNA from any known suspects or cannot obtain a hit in a DNA database. The HIrisPlex model for simultaneous eye and hair colour predictions was developed for forensic usage. To genotype a DNA sample, massively parallel sequencing (MPS) has brought new possibilities to the analysis of forensic DNA samples. As part of an in-house validation, this study presents the genotyping and predictive performance of the HIrisPlex SNPs in a Norwegian study population, using Verogen's ForenSeq™ DNA Signature Prep Kit on the MiSeq FGx system and the HIrisPlex webtool. DNA-profiles were successfully typed with DNA input down to 125 pg. In samples with DNA input < 125 pg, false homozygotes were observed with as many as 92 reads. Prediction accuracies in terms of AUC were high for red (0.97) and black (0.93) hair colours, as well as blue (0.85) and brown (0.94) eye colours. The AUCs for blond (0.72) and brown (0.70) hair colour were considerably lower. None of the individuals was predicted to have intermediate eye colour. Therefore, the error rates of the overall eye colour predictions were 37% with no predictive probability threshold (pmax) and 26% with a probability threshold of 0.7. We also observed that more than half of the incorrect predictions were for individuals carrying the rs12913832 GG genotype. For hair colour, 65% of the individuals were correctly predicted when using the highest probability category approach. The main error was observed for individuals with brown hair colour that were predicted to have blond hair. Utilising the prediction guide approach increased the correct predictions to 75%. Assessment of phenotype-genotype associations of eye colours using a quantitative eye colour score (PIE-score), revealed that rs12913832 AA individuals of Norwegian descent had statistically significantly higher PIE-score (less brown eye colour) than individuals of non-northern European descent. To our knowledge, this has not been reported in other studies. Our study suggests that careful assessment of the target population prior to the implementation of forensic DNA phenotyping to case work is beneficial.

Comprehensive Insights Into Forensic Features and Genetic Background of Chinese Northwest Hui Group Using Six Distinct Categories of 231 Molecular Markers

The Hui minority is predominantly composed of Chinese-speaking Islamic adherents distributed throughout China, of which the individuals are mainly concentrated in Northwest China. In the present study, we employed the length and sequence polymorphisms-based typing system of 231 molecular markers, i.e., amelogenin, 22 phenotypic-informative single nucleotide polymorphisms (PISNPs), 94 identity-informative single nucleotide polymorphisms (IISNPs), 24 Y-chromosomal short tandem repeats (Y-STRs), 56 ancestry-informative single nucleotide polymorphisms (AISNPs), 7 X-chromosomal short tandem repeats (X-STRs), and 27 autosomal short tandem repeats (A-STRs), into 90 unrelated male individuals from the Chinese Northwest Hui group to comprehensively explore its forensic characteristics and genetic background. Total of 451 length-based and 652 sequence-based distinct alleles were identified from 58 short tandem repeats (STRs) in 90 unrelated Northwest Hui individuals, denoting that the sequence-based genetic markers could pronouncedly provide more genetic information than length-based markers. The forensic characteristics and efficiencies of STRs and IISNPs were estimated, both of which externalized high polymorphisms in the Northwest Hui group and could be further utilized in forensic investigations. No significant departure from the Hardy-Weinberg equilibrium (HWE) expectation was observed after the Bonferroni correction. Additionally, four group sets of reference population data were exploited to dissect the genetic background of the Northwest Hui group separately from different perspectives, which contained 26 populations for 93 IISNPs, 58 populations for 17 Y-STRs, 26 populations for 55 AISNPs (raw data), and 109 populations for 55 AISNPs (allele frequencies). As a result, the analyses based on the Y-STRs indicated that the Northwest Hui group primarily exhibited intimate genetic relationships with reference Hui groups from Chinese different regions except for the Sichuan Hui group and secondarily displayed close genetic relationships with populations from Central and West Asia, as well as several Chinese groups. However, the AISNP analyses demonstrated that the Northwest Hui group shared more intimate relationships with current East Asian populations apart from reference Hui group, harboring the large proportion of ancestral component contributed by East Asia.

Exploring STR sequencing for forensic DNA intelligence databasing using the Austrian National DNA Database as an example

Here, we present the results from a population study that evaluated the performance of massively parallel sequencing (MPS) of short tandem repeats (STRs) with a particular focus on DNA intelligence databasing purposes. To meet this objective, 247 randomly selected reference samples, earlier being processed with conventional capillary electrophoretic (CE) STR sizing from the Austrian National DNA Database, were reanalyzed with the PowerSeq 46Y kit (Promega). This sample set provides MPS-based population data valid for the Austrian population to increase the body of sequence-based STR variation. The study addressed forensically relevant parameters, such as concordance and backward compatibility to extant amplicon-based genotypes, sequence-based stutter ratios, and relative marker performance. Of the 22 autosomal STR loci included in the PowerSeq 46GY panel, 99.98% of the allele calls were concordant between MPS and CE. Moreover, 25 new sequence variants from 15 markers were found in the Austrian dataset that are yet undescribed in the STRSeq online catalogue and were submitted for inclusion. Despite the high degree of concordance between MPS and CE derived genotypes, our results demonstrate the need for a harmonized allele nomenclature system that is equally applicable to both technologies, but at the same time can take advantage of the increased information content of MPS. This appears to be particularly important with regard to database applications in order to prevent false exclusions due to varying allele naming based on different analysis platforms and ensures backward compatibility.

Evaluation of ArmedXpert software tools, MixtureAce and Mixture Interpretation, to analyze MPS-STR data

Massively parallel sequencing (MPS) technologies have revolutionized studies of genomic variations and transformed DNA analysis in multiple fields. Assays based on MPS must be capable of discriminating variations introduced by the method, i.e. artifacts from true polymorphisms. In PCR-MPS methods targeting microsatellite markers, artifacts can arise from PCR mis-incorporation, PCR strand slippage (stutter), and sequencing error. Reliable detection of artifacts in mixed DNA samples is a significant challenge that must be addressed in forensic DNA analysis. The ArmedXpert (NicheVision) software tools, MixtureAce™ and Mixture Interpretation, can analyze MPS data by categorizing sequence reads in alleles, stutter, and non-stutter artifacts and analyzing autosomal STR loci of mixed samples. In this study, we evaluated the ArmedXpert tools for the analysis of STR profiles of single-sourced and mixed samples generated by the ForenSeq™ DNA Signature Prep kit (Verogen). Data from eight experimental runs (240 samples) were analyzed: one benchmark run, two runs testing sensitivity with down to 50 pg DNA input, one run testing artificially degraded samples and DNA derived from bones, blood cards and teeth, as well as four runs with mixed DNA samples of varying ratios, sex, and different number of contributors (two to six). The MixtureAce stutter thresholds were initially set following the recommendations from Verogen, plus a non-stutter artifact threshold was set at 5% of allele read counts. A benchmark run, of 30 samples, plus two controls, containing 2310 total alleles, revealed over 5000 artifacts, above an analytical threshold of 10. A total of 4869 artifacts were correctly classified, while 435 were mis-classified as alleles due to exceedance of initial threshold settings. False positives must be resolved by an analyst, which can be time consuming. Stutter thresholds were adjusted based on the benchmark data and the samples were re-tested, resulting in only 57 false positive allele calls. The revised settings were then used in the analysis of the remaining seven experimental runs. Results show that MixtureAce can accurately classify artifacts and alleles when laboratory-specific threshold settings are used. The Mixture Interpretation tool was applied on two- and three-person mixtures. This tool utilized the analyzed data from MixtureAce to calculate, based on the number of alleles at a locus and their read counts, possible deconvolution outcomes with their respective ratios.

Strategy for STR typing of bones from the Second World War combining CE and NGS technology: A pilot study

The genetic identification of skeletal remains found in Second World War mass graves is complicated because of the poor quality of the samples. The aim of this study was to set up a workflow for STR typing of such samples combining PCR/CE and PCR/NGS technologies. To this end, 57 DNA samples from an equal number of 75-year-old femurs were studied. After a first round of PCR typing using GlobalFiler CE, 42 samples yielded a full profile and were therefore submitted to our standard workflow. The 15 samples that yielded no or a limited number (2-17/21) of autosomal STR markers as well four bone control samples that provided a full profile with the conventional PCR/CE test were typed in duplicate by the GlobalFiler NGS kit. Despite the degradation of the samples, which resulted in lower coverage and a lower % of on-target reads, reliable sequencing data were obtained from 16/19 samples. The use of a threshold of 30× for the locus call led to a consensus profile (cp) of 20-31/31 STR autosomal loci in 10 samples and to a cp of 8-10/31 loci in two samples, whereas the four control samples yielded a cp of 26-31/31 loci. Finally, the data of the NGS typing were combined with those of the CE typing. This last task allowed us to recover (on average) three alleles per sample and to increase the number of the heterozygous patterns in 37 cases. In total, the combined approach proposed here made possible the genetic typing of 65-100% of the autosomal STR markers in 10/15 (66.6 %) skeletal remains that yielded no or very poor results with the conventional PCR/CE approach. However, because several artefacts (such as allelic drop-out and allelic drop-in) were scored, the risk of mistyping cannot be neglected.

A first for forensic genetics in Africa: successful identification of skeletal remains from the marine environment using massively parallel sequencing

In unrelated circumstances, two young adult males allegedly went missing off the coast of Cape Town, South Africa, within two months of each other. Weeks after the second disappearance, a decomposed human lower limb was recovered from a beach in Cape Town, followed by a washed-up decomposed hand three days later. An item of female clothing was found with the remains, and preliminary analysis of the skeleton indicated a female, leading to confusion regarding the possible identity of the decedent. Consequently, DNA analyses were requested to determine the biological sex of the remains, and whether the two sets of remains originated from the same individual. Various samples were collected, including bone, nails and swabs of soft tissue. DNA quantity and quality varied between sample types, with better results obtained from metacarpal bone and swab lysates. DNA profiling revealed a male sex, which suggests cognitive bias may have played a role in initial sex estimations. In addition, massively parallel sequencing confidently matched the two sets of remains (random match probability: 1 in 2.70 x 10³¹). These results were a first for Africa where massively parallel sequencing was successfully used and assisted in the identification of human remains, thus, affording closure to the next-of-kin. Moreover, this constitutes the first global report where soft tissue lysates from a marine decomposition case yielded full DNA profiles with a massively parallel sequencing approach.

Identification of sequence polymorphisms at 58 STRs and 94 iiSNPs in a Tibetan population using massively parallel sequencing

Massively parallel sequencing (MPS) has rapidly become a promising method for forensic DNA typing, due to its ability to detect a large number of markers and samples simultaneously in a single reaction, and sequence information can be obtained directly. In the present study, two kinds of forensic genetic markers, short tandem repeat (STR) and identity-informative single nucleotide polymorphism (iiSNP) were analyzed simultaneously using ForenSeq DNA Signature Prep Kit, a commercially available kit on MPS platform. A total of 152 DNA markers, including 27 autosomal STR (A-STR) loci, 24 Y chromosomal STR (Y-STR) loci, 7 X chromosomal STR (X-STR) loci and 94 iiSNP loci were genotyped for 107 Tibetan individuals (53 males and 54 females). Compared with length-based STR typing methods, 112 more A-STR alleles, 41 more Y-STR alleles, and 24 more X-STR alleles were observed at 17 A-STRs, 9 Y-STRs, and 5 X-STRs using sequence-based approaches. Thirty-nine novel sequence variations were observed at 20 STR loci. When the flanking regions were also analyzed in addition to target SNPs at the 94 iiSNPs, 38 more alleles were identified. Our study provided an adequate genotype and frequencies data of the two types of genetic markers for forensic practice. Moreover, we also proved that this panel is highly polymorphic and informative in Tibetan population, and should be efficient in forensic kinship testing and personal identification cases.

Classification of STR allelic variation using massively parallel sequencing and assessment of flanking region power

The application of massively parallel sequencing (MPS) to forensic genetics has led to improvements in multiple aspects of DNA analysis, however, additional complexities are concurrently associated with these advances. In relation to short tandem repeat (STR) typing, the move to sequence rather than length-based methodologies has highlighted the extent to which previous allelic variation was masked - both within and outside of the repeat regions (the flanking regions). In order to fully implement MPS for autosomal STR analysis, sequence-based allelic frequencies must be available, and concordance with previous typing techniques needs to be assessed. In this work, a series of samples (n = 1007) from five different population groups were genotyped using the MiSeq FGx™ Forensic Genomics System. Results were compared to those obtained using capillary electrophoresis (CE), and sequence variation has been characterised both within and outside STR repeat regions, with allelic frequencies provided for all variants observed within this database. Analysing and characterising flanking region sequence is currently less straightforward than studying repeat region variation alone, and the added value of doing so remains largely unexplored - this paper provides data to show that the gain in polymorphism achieved when analysing flanking regions is less than might be expected. In the White British population for example, including the sequence variation within repeat regions of 26 autosomal STRs made the average combined random match probability (RMP) over 700 times lower than with length-based alleles alone. Including the sequence variation within the flanking regions only resulted in a combined RMP that was a further 4 times lower.

Massive parallel sequencing in forensics: advantages, issues, technicalities, and prospects

In the last decade, next-generation sequencing (NGS) technology, alternatively massive parallel sequencing (MPS), was applied to all fields of biological research. Its introduction to the field of forensics was slower, mainly due to lack of accredited sequencers, kits, and relatively higher sequencing error rates as compared with standardized Sanger sequencing. Currently, a majority of the problematic issues have been solved, which is proven by the body of reports in the literature. Here, we discuss the utility of NGS sequencing in forensics, emphasizing the advantages, issues, the technical aspects of the experiments, commercial solutions, and the potentially interesting applications of MPS.

Inter-laboratory study on standardized MPS libraries: evaluation of performance, concordance, and sensitivity using mixtures and degraded DNA

We present results from an inter-laboratory massively parallel sequencing (MPS) study in the framework of the SeqForSTRs project to evaluate forensically relevant parameters, such as performance, concordance, and sensitivity, using a standardized sequencing library including reference material, mixtures, and ancient DNA samples. The standardized library was prepared using the ForenSeq DNA Signature Prep Kit (primer mix A). The library was shared between eight European laboratories located in Austria, France, Germany, The Netherlands, and Sweden to perform MPS on their particular MiSeq FGx sequencers. Despite variation in performance between sequencing runs, all laboratories obtained quality metrics that fell within the manufacturer’s recommended ranges. Furthermore, differences in locus coverage did not inevitably adversely affect heterozygous balance. Inter-laboratory concordance showed 100% concordant genotypes for the included autosomal and Y-STRs, and still, X-STR concordance exceeded 83%. The exclusive reasons for X-STR discordances were drop-outs at DXS10103. Sensitivity experiments demonstrated that correct allele calling varied between sequencing instruments in particular for lower DNA amounts (≤ 125 pg). The analysis of compromised DNA samples showed the drop-out of one sample (FA10013B01A) while for the remaining three degraded DNA samples MPS was able to successfully type ≥ 87% of all aSTRs, ≥ 78% of all Y-STRs, ≥ 68% of all X-STRs, and ≥ 92% of all iSNPs demonstrating that MPS is a promising tool for human identity testing, which in return, has to undergo rigorous in-house validation before it can be implemented into forensic routine casework.

Massively parallel sequencing of autosomal STRs and identity-informative SNPs highlights consanguinity in Saudi Arabia

While many studies have been undertaken of Middle Eastern populations using autosomal STR profiling by capillary electrophoresis, little has so far been published from this region on the forensic use of massively parallel sequencing (MPS). Here, we carried out MPS of 27 autosomal STRs and 91 identity-informative SNPs (iiSNPs) with the Verogen ForenSeq™ DNA Signature Prep Kit on a representative sample of 89 Saudi Arabian males, and analysed the resulting sequence data using Verogen's ForenSeq Universal Analysis Software (UAS) v1.3 and STRait Razor v3.0. This revealed sequence variation in the composition of complex STR arrays, and SNPs in their flanking regions, which raised the number of STR alleles from 238 distinct length variants to 357 sequence sub-variants. Similarly, between one and three additional polymorphic sites were observed within the amplicons of 37 of the 91 iiSNPs, forming up to six microhaplotypes per locus. These further enhance discrimination compared to the biallelic target SNP data presented by the primary UAS interface. In total, we observed twenty-two STR alleles previously unrecognised in the STRait Razor v3.0 default allele list, along with nine SNPs flanking target iiSNPs that were not highlighted by the UAS. Sequencing reduced the STR-based random match probability (RMP) from 2.62E-30 to 3.49E-34, and analysis of the iiSNP microhaplotypes reduced RMP from 9.97E-37 to 6.83E-40. The lack of significant linkage disequilibrium between STRs and target iiSNPs allowed the two marker types to be combined using the product rule, yielding a RMP of 2.39E-73. Evidence of consanguinity was apparent from both marker types. While TPOX was the only locus displaying a significant deviation from Hardy-Weinberg equilibrium, 23 out of 27 STRs and 63 out of 91 iiSNPs showed fewer than expected heterozygotes, demonstrating an overall homozygote excess probably reflecting the high frequency of first-cousin marriages in Saudi Arabia. We placed our data in a global context by considering the same markers in the Human Genome Diversity Panel (HGDP), revealing that the Saudi sample was typical of Middle Eastern populations, with a higher level of inbreeding than is seen in most European, African and Central/South Asian populations, correlating with known patterns of endogamy. Given reduced levels of diversity within endogamous groups, the ability to combine the discrimination power of both STRs and SNPs offers significant benefits in the analysis of forensic evidence in Saudi Arabia and the Middle East region more generally.

Is MPS always the answer? Use of two PCR-based methods for Y-chromosomal haplotyping in highly and moderately degraded bone material

Forensic and population genetics often rely on Y-chromosomal studies. Whether it is a human identification case, trace evidence examination or phylogenetic analysis, a Y-STR haplotype is an important tool in the hands of law enforcement agencies. A common obstacle in achieving satisfactory results in all of the above mentioned circumstances, is low DNA quantity and quality within samples obtained. In this study we have examined Y-STR haplotypes in 75 bone material samples, coming from different time periods. For this purpose we have chosen YFiler Plus PCR Amplification Kit (ThermoFisher Scientific) and ForenSeq Signature DNA Prep Kit (Verogen Inc.), which use two different allele calling technologies - capillary electrophoresis and Massively Parallel Sequencing respectively. Full profiles were obtained from DNA extracts with as little as 0.1896 ng (Degradation Index 1.3) (ForenSeq) and 0.0591 ng (Degradation Index 26.8) (YFiler Plus) DNA input. The results that we present in this paper show differences in amplification rates between common markers in both kits. The differences strictly reflect mean amplicon length of markers. This, however, does not seem to influence Y-haplogroup estimation results noticeably. In one sample a discordance occurred between haplotypes obtained with both methods, where a 24 allele was called in DYS390 marker by capillary electrophoresis, while for the same sample in this locus a 23 allele was shown with MPS. A reason for this is yet to be investigated. The sequence analysis revealed a significant variation between isometric alleles, especially within repetitive regions of studied Y-STR markers.

Impact of DNA degradation on massively parallel sequencing-based autosomal STR, iiSNP, and mitochondrial DNA typing systems

Biological samples, including skeletal remains exposed to environmental insults for extended periods of time, exhibit increasing levels of DNA damage and fragmentation. Human forensic identification methods typically use a combination of mitochondrial (mt) DNA sequencing and short tandem repeat (STR) analysis, which target segments of DNA ranging from 80 to 500 base pairs (bps). Larger templates are often unavailable as skeletal samples age and the associated DNA degrades. Single-nucleotide polymorphism (SNP) loci target shorter templates and may serve as a solution to the problem. Recently developed assays for STR and SNP analysis using a massively parallel sequencing approach, such as the ForenSeq kit (Verogen, San Diego, CA), offer a means for generating results from degraded samples as they target templates down to 60 to 170 bps. We performed a modeling study that demonstrates that SNPs can increase the significance of an identification when analyzing DNA down to an average size of 100 bps for input amounts between 0.375 and 1 ng of nuclear DNA. Observations from this study were then compared with human skeletal material results (n = 14, ninth to eighteenth centuries), which further demonstrated the utility of the ForenSeq kit for degraded samples. The robustness of the Promega PowerSeq™ Mito System was also tested with human skeletal remains (n = 70, ninth to eighteenth centuries), resulting in successful coverage of 99.29% of the mtDNA control region at 50× coverage or more. This was accompanied by modifications to a mainstream DNA extraction technique for skeletal remains that improved recovery of shorter templates.

The impact of common PCR inhibitors on forensic MPS analysis

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.

Massively parallel sequencing techniques for forensics: A review

DNA sequencing, starting with Sanger's chain termination method in 1977 and evolving into the next generation sequencing (NGS) techniques of today that employ massively parallel sequencing (MPS), has become essential in application areas such as biotechnology, virology, and medical diagnostics. Reflected by the growing number of articles published over the last 2-3 years, these techniques have also gained attention in the forensic field. This review contains a brief description of first, second, and third generation sequencing techniques, and focuses on the recent developments in human DNA analysis applicable in the forensic field. Relevance to the forensic analysis is that besides generation of standard STR-profiles, DNA repeats can also be sequenced to look for polymorphisms. Furthermore, additional SNPs can be sequenced to acquire information on ancestry, paternity or phenotype. The current MPS systems are also very helpful in cases where only a limited amount of DNA or highly degraded DNA has been secured from a crime scene. If enough autosomal DNA is not present, mitochondrial DNA can be sequenced for maternal lineage analysis. These developments clearly demonstrate that the use of NGS will grow into an indispensable tool for forensic science.

An automated independent workflow for the analysis of massively parallel sequence data from forensic SNP assays

Illumina and Thermo Fisher Scientific have developed assays that permit the sequencing of forensically relevant single nucleotide polymorphisms (SNPs), along with software to determine the associated genotypes. Currently there is no method to either independently confirm the genotypes determined using the manufacturer's software, or to compare genotypes and quality metrics among samples processed using both platforms. This paper outlines an automated workflow developed in CLC Genomics Workbench that permits accurate, fast and independent analysis of SNP sequence data from either platform. To facilitate the straightforward comparison of genotypes generated from both the manufacturer's software and the independent CLC analysis, a Python script was written. Data for a total of 323 forensically relevant ancestry, identity and phenotypic SNPs can be analyzed, and the resulting genotypes, coverage, quality flags and major allele frequencies are easily compared across samples and platforms.

Sequencing of 231 forensic genetic markers using the MiSeq FGx™ forensic genomics system - an evaluation of the assay and software

The MiSeq FGx™ Forensic Genomics System types 231 genetic markers in one multiplex polymerase chain reaction (PCR) assay. The markers include core forensic short tandem repeats (STRs) as well as identity, ancestry and phenotype informative short nucleotide polymorphisms (SNPs). In this work, the MiSeq FGx™ Forensic Genomics System was evaluated by analysing reproducibility, sensitivity, mixture identification and forensic phenotyping capabilities of the assay. Furthermore, the genotype calling of the ForenSeq™ Universal Analysis Software was verified by analysing fastq.gz files from the MiSeq FGx™ platform using the softwares STRinNGS and GATK. Overall, the performance of the MiSeq FGx™ Forensic Genomics System was high. However, locus and allele drop-outs were relatively frequent at six loci (two STRs and four human identification SNPs) due to low read depth or skewed heterozygote balances, and the stutter ratios were larger than those observed with conventional STR genotyping methods. The risk of locus and allele drop-outs increased dramatically when the amount of DNA in the first PCR was lower than 250 pg. Two-person 50:1 mixtures were identified as mixtures, whereas 100:1 and 1 000:1 mixtures were not. Y-chromosomal short tandem repeats (Y-STRs) alleles were detected in the 100:1 and 1 000:1 female/male mixtures. The ForenSeq™ Universal Analysis Software provided the data analyst with useful alerts that simplified the analysis of the large number of markers. Many of the alerts were due to user-defined, locus-specific criteria. The results shown here indicated that the default settings should be altered for some loci. Also, recommended changes to the assay and software are discussed.

Qualitative and quantitative assessment of Illumina's forensic STR and SNP kits on MiSeq FGx™

Massively parallel sequencing (MPS) is a powerful tool transforming DNA analysis in multiple fields ranging from medicine, to environmental science, to evolutionary biology. In forensic applications, MPS offers the ability to significantly increase the discriminatory power of human identification as well as aid in mixture deconvolution. However, before the benefits of any new technology can be employed, a thorough evaluation of its quality, consistency, sensitivity, and specificity must be rigorously evaluated in order to gain a detailed understanding of the technique including sources of error, error rates, and other restrictions/limitations. This extensive study assessed the performance of Illumina’s MiSeq FGx MPS system and ForenSeq™ kit in nine experimental runs including 314 reaction samples. In-depth data analysis evaluated the consequences of different assay conditions on test results. Variables included: sample numbers per run, targets per run, DNA input per sample, and replications. Results are presented as heat maps revealing patterns for each locus. Data analysis focused on read numbers (allele coverage), drop-outs, drop-ins, and sequence analysis. The study revealed that loci with high read numbers performed better and resulted in fewer drop-outs and well balanced heterozygous alleles. Several loci were prone to drop-outs which led to falsely typed homozygotes and therefore to genotype errors. Sequence analysis of allele drop-in typically revealed a single nucleotide change (deletion, insertion, or substitution). Analyses of sequences, no template controls, and spurious alleles suggest no contamination during library preparation, pooling, and sequencing, but indicate that sequencing or PCR errors may have occurred due to DNA polymerase infidelities. Finally, we found utilizing Illumina’s FGx System at recommended conditions does not guarantee 100% outcomes for all samples tested, including the positive control, and required manual editing due to low read numbers and/or allele drop-in. These findings are important for progressing towards implementation of MPS in forensic DNA testing.

Investigation of the STR loci noise distributions of PowerSeq™ Auto System

Aim

To characterize the noise and stutter distribution of 23 short tandem repeats (STRs) included in the PowerSeq^TM Auto System.

Methods

Raw FASTQ files were analyzed using STRait Razor v2s to display alleles and coverage. The sequence noise was divided into several categories: noise at allele position, noise at -1 repeat position, and artifact. The average relative percentages of locus coverage for each noise, stutter, and allele were calculated from the samples used for this locus noise analysis.

Results

Stutter products could be routinely observed at the -2 repeat position, -1 repeat position, and +1 repeat position of alleles. Sequence noise at the allele position ranged from 10.22% to 28.81% of the total locus coverage. At the allele position, individual noise reads were relatively low.

Conclusion

The data indicate that noise generally will be low. In addition, the PowerSeq^TM Auto System could capture nine flanking region single nucleotide polymorphisms (SNPs) that would not be observed by other current kits for massively parallel sequencing (MPS) of STRs.

Population and performance analyses of four major populations with Illumina's FGx Forensic Genomics System

The MiSeq FGx Forensic Genomics System (Illumina) enables amplification and massively parallel sequencing of 59 STRs, 94 identity informative SNPs, 54 ancestry informative SNPs, and 24 phenotypic informative SNPs. Allele frequency and population statistics data were generated for the 172 SNP loci included in this panel on four major population groups (Chinese, African Americans, US Caucasians, and Southwest Hispanics). Single-locus and combined random match probability values were generated for the identity informative SNPs. The average combined STR and identity informative SNP random match probabilities (assuming independence) across all four populations were 1.75E-67 and 2.30E-71 with length-based and sequence-based STR alleles, respectively. Ancestry and phenotype predictions were obtained using the ForenSeq™ Universal Analysis System (UAS; Illumina) based on the ancestry informative and phenotype informative SNP profiles generated for each sample. Additionally, performance metrics, including profile completeness, read depth, relative locus performance, and allele coverage ratios, were evaluated and detailed for the 725 samples included in this study. While some genetic markers included in this panel performed notably better than others, performance across populations was generally consistent. The performance and population data included in this study support that accurate and reliable profiles were generated and provide valuable background information for laboratories considering internal validation studies and implementation.

Comparison between magnetic bead and qPCR library normalisation methods for forensic MPS genotyping

Massively parallel sequencing (MPS) is fast approaching operational use in forensic science, with the capability to analyse hundreds of DNA identity and DNA intelligence markers in multiple samples simultaneously. The ForenSeq™ DNA Signature Kit on MiSeq FGx™ (Illumina) workflow can provide profiles for autosomal short tandem repeats (STRs), X chromosome and Y chromosome STRs, identity single nucleotide polymorphisms (SNPs), biogeographical ancestry SNPs and phenotype (eye and hair colour) SNPs from a sample. The library preparation procedure involves a series of steps including target amplification, library purification and library normalisation. This study highlights the comparison between the manufacturer recommended magnetic bead normalisation and quantitative polymerase chain reaction (qPCR) methods. Furthermore, two qPCR chemistries, KAPA® (KAPA Biosystems) and NEBNext® (New England Bio Inc.), have also been compared. The qPCR outperformed the bead normalisation method, while the NEBNext® kit obtained higher genotype concordance than KAPA®. The study also established an MPS workflow that can be utilised in any operational forensic laboratory.