Characterizing stutter variants in forensic STRs with massively parallel sequencing

Development and validation of YARN: A novel SE-400 MPS kit for East Asian paternal lineage analysis

Y-chromosomal short tandem repeat polymorphisms (Y-STRs) and Y-chromosomal single nucleotide polymorphisms (Y-SNPs) are valuable genetic markers used in paternal lineage identification and population genetics. Currently, there is a lack of an effective panel that integrates Y-STRs and Y-SNPs for studying paternal lineages, particularly in East Asian populations. Hence, we developed a novel Y-chromosomal targeted panel called YARN (Y-chromosome Ancestry and Region Network) based on multiplex PCR and a single-end 400 massive parallel sequencing (MPS) strategy, consisting of 44 patrilineage Y-STRs and 260 evolutionary Y-SNPs. A total of 386 reactions were validated for the effectiveness and applicability of YARN according to SWGDAM validation guidelines, including sensitivity (with a minimum input gDNA of 0.125 ng), mixture identification (ranging from 1:1-1:10), PCR inhibitor testing (using substances such as 50 μM hematin, 100 μM hemoglobin, 100 μM humic acid, and 2.5 mM indigo dye), species specificity (successfully distinguishing humans from other animals), repeatability study (achieved 100% accuracy), and concordance study (with 99.91% accuracy for 1121 Y-STR alleles). Furthermore, we conducted a pilot study using YARN in a cohort of 484 Han Chinese males from Huaiji County, Zhaoqing City, Guangdong, China (GDZQHJ cohort). In this cohort, we identified 52 different Y-haplogroups and 73 different surnames. We found weak to moderate correlations between the Y-haplogroups, Chinese surnames, and geographical locations of the GDZQHJ cohort (with λ values ranging from 0.050 to 0.340). However, when we combined two different categories into a new independent variable, we observed stronger correlations (with λ values ranging from 0.617 to 0.754). Overall, the YARN panel, which combines Y-STR and Y-SNP genetic markers, meets forensic DNA quality assurance guidelines and holds potential for East Asian geographical origin inference and paternal lineage analysis.

Ultrasensitive sequencing of STR markers utilizing unique molecular identifiers and the SiMSen-Seq method

Massively parallel sequencing (MPS) is increasingly applied in forensic short tandem repeat (STR) analysis. The presence of stutter artefacts and other PCR or sequencing errors in the MPS-STR data partly limits the detection of low DNA amounts, e.g., in complex mixtures. Unique molecular identifiers (UMIs) have been applied in several scientific fields to reduce noise in sequencing. UMIs consist of a stretch of random nucleotides, a unique barcode for each starting DNA molecule, that is incorporated in the DNA template using either ligation or PCR. The barcode is used to generate consensus reads, thus removing errors. The SiMSen-Seq (Simple, multiplexed, PCR-based barcoding of DNA for sensitive mutation detection using sequencing) method relies on PCR-based introduction of UMIs and includes a sophisticated hairpin design to reduce unspecific primer binding as well as PCR protocol adjustments to further optimize the reaction. In this study, SiMSen-Seq is applied to develop a proof-of-concept seven STR multiplex for MPS library preparation and an associated bioinformatics pipeline. Additionally, machine learning (ML) models were evaluated to further improve UMI allele calling. Overall, the seven STR multiplex resulted in complete detection and concordant alleles for 47 single-source samples at 1 ng input DNA as well as for low-template samples at 62.5 pg input DNA. For twelve challenging mixtures with minor contributions of 10 pg to 150 pg and ratios of 1-15% relative to the major donor, 99.2% of the expected alleles were detected by applying the UMIs in combination with an ML filter. The main impact of UMIs was a substantially lowered number of artefacts as well as reduced stutter ratios, which were generally below 5% of the parental allele. In conclusion, UMI-based STR sequencing opens new means for improved analysis of challenging crime scene samples including complex mixtures.

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.

Whole Mitochondrial Genome Detection and Analysis of Two- to Four-Generation Maternal Pedigrees Using a New Massively Parallel Sequencing Panel

Mitochondrial DNA (mtDNA) is an effective genetic marker in forensic practice, especially for aged bones and hair shafts. Detection of the whole mitochondrial genome (mtGenome) using traditional Sanger-type sequencing is laborious and time-consuming. Additionally, its ability to distinguish point heteroplasmy (PHP) and length heteroplasmy (LHP) is limited. The application of massively parallel sequencing in mtDNA detection helps researchers to study the mtGenome in-depth. The ForenSeq mtDNA Whole Genome Kit, which contains a total of 245 short amplicons, is one of the multiplex library preparation kits for the mtGenome. We used this system to detect the mtGenome in the blood samples and hair shafts of thirty-three individuals from eight two-generation pedigrees, one three-generation pedigree, and one four-generation pedigree. High-quality sequencing results were obtained. Ten unique mtGenome haplotypes were observed in the mothers from the ten pedigrees. A total of 26 PHPs were observed using the interpretation threshold of 6%. Eleven types of LHPs in six regions were evaluated in detail. When considering homoplasmic variants only, consistent mtGenome haplotypes were observed between the twice-sequenced libraries and between the blood and hair shafts from the same individual and among maternal relatives in the pedigrees. Four inherited PHPs were observed, and the remainder were de novo/disappearing PHPs in the pedigrees. Our results demonstrate the effective capability of the ForenSeq mtDNA Whole Genome Kit to generate the complete mtGenome in blood and hair shafts, as well as the complexity of mtDNA haplotype comparisons between different types of maternal relatives when heteroplasmy is considered.

Characterizing the amplification of STR markers in multiplex polymerase chain displacement reaction using massively parallel sequencing

Polymerase chain displacement reaction (PCDR) showed advantages in forensic low-template DNA analysis with improved amplification efficiency, higher allele detection capacity, and lower stutter artifact than PCR. However, characteristics of STR markers after PCDR amplification remain unclarified for the limited resolving power of capillary electrophoresis (CE). This issue can be addressed by massively parallel sequencing (MPS) technology with higher throughput and discriminability. Here, we developed a multiplex PCDR system including 24 STRs and amelogenin. In addition, a PCR reference was established for comparison. After amplification, products were subjected to PCR-free library construction and sequenced on the Illumina NovaSeq system. We implemented a sequence-matching pipeline to separate different amplicon types of PCDR products from the combination of primers. In the sensitivity test, the PCDR multiplex obtained full STR profiles with as low as 125 pg 2800M control DNA. Based on that, single-source DNA samples were tested. First, highly concordant genotypes were observed among the PCDR multiplex, the PCR reference, and CE-based STR kits. Next, read counts of different PCDR amplicon types were investigated, showing a relative abundance of 78:12:12:1 for the shortest amplicon S, the two medium amplicons M1 and M2, and the longest amplicon L. We also analyzed the stutter artifacts for distinct amplicon types, and the results revealed the reduction of N - 1 and N - 2 contraction stutters, and the increase of N + 1 and N + 2 elongation stutters in PCDR samples. Moreover, we confirmed the feasibility of PCDR for amplifying degraded DNA samples and unbalanced DNA mixtures. Compared to the previous proof of principle study, our work took a further step to characterize the complete profile of STR markers in the PCDR context. Our results suggested that the PCDR-MPS workflow is an effective approach for forensic STR analysis. Corresponding findings in this study may help the development of PCDR-based assays and probabilistic methods in future studies.

Phosphoryl guanidine oligonucleotides as primers for RNA-dependent DNA synthesis using murine leukemia virus reverse transcriptase

Modern approaches to the detection and analysis of low-copy-number RNAs are often based on the use of RNA-dependent DNA polymerases, for example, in reverse-transcription PCR. The accuracy and eff iciency of cDNA synthesis in the reverse-transcription reaction catalyzed by reverse transcriptase (RNA-dependent DNA polymerase) signif icantly affect the correctness of the results of PCR diagnostic assays and/or RNA sequencing. In this regard, many studies are focused on the optimization of the reverse-transcription reaction, including the search for more perfect primers necessary to obtain a full-length DNA copy of RNA under study. The best-known completely uncharged analogs of oligonucleotides – morpholine oligonucleotides and peptide nucleic acids – cannot be substrates for enzymes that process nucleic acids. The aim of this work was to conduct a pilot study of uncharged phosphoryl guanidine oligodeoxyribonucleotides (PGOs) as primers for mouse leukemia virus reverse transcriptase (MMLV H-). Specif ic features of elongation of partially and completely uncharged PGO primers were investigated. It was demonstrated that PGOs can be elongated eff iciently, e. g., in the presence of a fragment of human ribosomal RNA having complex spatial structure. It was shown that the proportion (%) of abortive elongation products of a PGO primer depends on buffer ionic strength, nucleotide sequence of the primer, and the presence and location of phosphoryl guanidine groups in the primer. The results indicate the suitability of PGOs, including completely electroneutral ones, as primers for reverse-transcription PCR, thereby opening up new prospects for the creation of experimental models for the analysis of highly structured RNA.

Kinship assignment with the ForenSeq™ DNA Signature Prep Kit: Sources of error in simulated and real cases

Kinship recognition between anonymous DNA samples is becoming a relevant issue in forensics, more so with the increasing number of DNA profiles in databanks. Also, NGS-based genotyping is being increasingly used in routine personal identification, to simultaneously type large numbers of markers of different kind. In the present work, we explored computationally and experimentally the performance of the ForenSeq™ DNA Signature Prep Kit in identifying the true relationship between two anonymous samples, distinguishing it from other possible relationships. We analyzed with Familias R series of 10,000 pairs with 9 different simulated relationships, corresponding to different degrees of autosomal sharing. For each pair we obtained likelihood ratios for five kinship hypotheses vs. unrelatedness, and used their ranking to identify the preferred relationship. We also typed 21 subjects from two pedigrees, representing from parent-child to 4th cousins relationships. As expected, the power for identifying the true relationship decays in the order of autosomal sharing. Parent-child and full siblings can be robustly identified against other relationships. For half-siblings the chance of reaching a significant conclusion is already small. For more distant relationships the proportion of cases correctly and significantly identified is 10% or less. Bidirectional errors in kinship attribution include the suggestion of relatedness when this does not exist (10-50%), and the suggestion of independence in pairs of individuals more than 4 generations apart (25-60%). The real cases revealed a relevant effect of genotype miscalling at some loci, which could only be partly avoided by modulating the analysis parameters. In conclusion, with the exception of first degree relatives, the kit can be useful to inform additional investigations, but does not usually provide probatory results.

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.

A proof-of-principle study on implementing polymerase chain displacement reaction (PCDR) to improve forensic low-template DNA analysis

Polymerase chain reaction (PCR) plays an important role in forensic DNA analysis. However, the amplification of low-template DNA (LTDNA) samples usually encounters unsatisfactory results for the limited efficiency of PCR, which would interfere with the subsequent profile interpretation. Polymerase chain displacement reaction (PCDR) is a highly-efficient technique characterized by combining PCR and strand displacement reaction into a single PCDR cycle. This study explored the feasibility of PCDR for improving forensic LTDNA analysis. STR markers commonly used in forensic genetics were subjected to PCDR amplification and capillary electrophoresis detection. The results of singleplex reactions indicated that PCDR surpassed original PCR in efficiency for STR amplification. The average peak height of alleles in PCDR profiles was linearly correlated to the number of outer primers adopted for initiating the strand displacement process. Further, we assessed the multiplexing potential of PCDR by incorporating 17 STRs included in the expanded CODIS core loci and Amelogenin gene into a multiplex PCDR system. For pristine DNA templates ranged from 200 pg to 12.5 pg, the multiplex PCDR system consistently exhibited higher allele peak height as well as less allele dropout compared to the multiplex PCR references. Meanwhile, a significant reduction of stutter ratio was extensively observed in PCDR profiles. We also tested mock casework samples to verify the practical ability of multiplex PCDR for LTDNA detection. With DNA input varying from 48.1 pg to 6.6 pg, the multiplex PCDR system consistently obtained more allelic information than multiplex PCR methods. Our data collectively suggested that it is feasible to apply PCDR in forensic LTDNA analysis.

Validation of a 6-Dye Short Tandem Repeat System: A Dry Kit With Lyophilized Amplification Reagent

The SureID®S6 system used a lyophilized pellet as the amplification reagent to enable multiplexing of sex-determining marker Amelogenin, 21 autosomal short tandem repeats (STRs), and one Y-STR. To assess the performance, reliability, and limitation of the dry amplification system, the validation studies including PCR condition, reproducibility, sizing and precision, analytical threshold calculation, sensitivity and stochastic threshold calculation, species specificity, stability, mixture, case sample, and population and concordance were conducted according to the Scientific Working Group on DNA Analysis Methods (SWGDAM) Validation Guidelines. Experimental data suggested that the optimal range of total input DNA was from 125 to 500 pg; the appropriate analytical threshold was 80 relative fluorescence units (RFUs) while the stochastic threshold was 260 RFUs; for the stability studies, SureID®S6 system could resist against less than 500 μmol/L of hematin, 100 ng/μl of humic acid, 4 mM of indigotin, 800 mM of tannic acid, and 800 mM of calcium ion. Population and concordance studies using 500 unrelated individuals showed that the combined probability of discrimination (CPD) and cumulative probability of exclusion (CPE) values were 0.999999999999 and 0.999999998416, respectively. The genotypes for the same sample were concordant with the previously validated HUAXIA™ Platinum kit. The validation results demonstrated that the SureID®S6 system could be used for forensic applifications.

Pairwise kinship testing with microhaplotypes: Can advancements be made in kinship inference with these markers?

Kinship testing based on genetic relatedness is one of the major tasks in forensic genetics. Although short tandem repeats (STRs) are the "gold standard" biomarkers for relationship testing, microhaplotypes (MHs) have also emerged as viable options for kinship elucidation. In this work, the kinship testing efficiency of 54 highly polymorphic MHs was studied in two extended families consisting of parent-offspring, full siblings, grandparent-grandchildren, uncle/aunt-nephew/nieces, and first cousins. In addition, ten-thousand pairs of different degrees of relationships were simulated using various datasets including 54 MHs, 27 STRs plus 94 single nucleotide polymorphisms (SNPs) that were included in the ForenSeq DNA Signature Prep Kit (ForenSeq), 54 MHs plus loci in ForenSeq, and different subsets of 417-published MHs. The panels' system effectiveness in the kinship analysis were accessed by likelihood ratio distributions. The results showed that 54 MHs could be used in first-degree relationship testing with high reliability. The effectiveness of 54 MHs was slightly lower than ForenSeq but only by a narrow margin. Both 54 MHs and ForenSeq were not sufficient for distant relationship testing, and approximately 200 microhaplotypes with an average expected heterozygosity (He) = 0.79 were enough to determine second-degree relationships, but a panel of 417 MHs with an average He = 0.72 was not sufficient to first cousins testing according to the simulation analysis. In conclusion, 54 MHs could be used to serve as supplement markers for kinship testing; and well-established STR markers plus well-performing microhaplotype markers may become collective tools in forensic applications, though an enlarged pool of forensic markers is needed for distant relationship testing.

The forensic landscape and the population genetic analyses of Hainan Li based on massively parallel sequencing DNA profiling

Due to the formation of the Qiongzhou Strait by climate change and marine transition, Hainan island was isolated from the mainland southern China during the Last Glacial Maximum. Hainan island, located at the southernmost part of China and separated from the Leizhou Peninsula by the Qiongzhou Strait, laid on one of the modern human northward migration routes from Southeast Asia to East Asia. The Hlai language-speaking Li minority, the second largest population after Han Chinese in Hainan island, is the direct descendants of the initial migrants in Hainan island and has unique ethnic properties and derived characteristics; however, the forensic-associated studies on Hainan Li population are still insufficient. Hence, 136 Hainan Li individuals were genotyped in this study using the MPS-based ForenSeq™ DNA Signature Prep Kit (DNA Primer Set A, DPMA) to characterize the forensic genetic polymorphism landscape, and DNA profiles were obtained from 152 different molecular genetic markers (27 autosomal STRs, 24 Y-STRs, 7 X-STRs, and 94 iiSNPs). A total of 419 distinct length variants and 586 repeat sequence sub-variants, with 31 novel alleles (at 17 loci), were identified across the 58 STR loci from the DNA profiles of Hainan Li population. We evaluated the forensic characteristics and efficiencies of DPMA, demonstrating that the STRs and iiSNPs in DPMA were highly polymorphic in Hainan Li population and could be employed in forensic applications. In addition, we set up three datasets, which included the genetic data of (i) iiSNPs (27 populations, 2640 individuals), (ii) Y-STRs (42 populations, 8281 individuals), and (iii) Y haplogroups (123 populations, 4837 individuals) along with the population ancestries and language families, to perform population genetic analyses separately from different perspectives. In conclusion, the phylogenetic analyses indicated that Hainan Li, with a southern East Asia origin and Tai-Kadai language-speaking language, is an isolated population relatively. But the genetic pool of Hainan Li influenced by the limited gene flows from other Tai-Kadai populations and Hainan populations. Furthermore, the establishment of isolated population models will be beneficial to clarify the exquisite population structures and develop specific genetic markers for subpopulations in forensic genetic fields.

Noninvasive Prenatal Paternity Testing with a Combination of Well-Established SNP and STR Markers Using Massively Parallel Sequencing

Cell-free fetal DNA (cffDNA) from maternal plasma has made it possible to develop noninvasive prenatal paternity testing (NIPPT). However, most studies have focused on customized single nucleotide polymorphism (SNP) typing systems and few have used conventional short tandem repeat (STR) markers. Based on massively parallel sequencing (MPS), this study used a widely-accepted forensic multiplex assay system to evaluate the effect of noninvasive prenatal paternity testing with a combination of well-established SNP and STR markers. Using a ForenSeq DNA Signature Prep Kit, NIPPT was performed in 17 real parentage cases with monovular unborn fetuses at 7 to 24 gestational weeks. Different analytical strategies for the identification of paternally inherited allele (PIA) were developed to deal with SNPs and STRs. Combined paternity index (CPI) for 17 real trios as well as 272 unrelated trios was calculated. With the combination of SNPs and A-STRs, 82.35% (14/17), 88.24% (15/17), 94.12% (16/17), and 94.12% (16/17) of real trios could be accurately determined when the likelihood ratio (LR) threshold for paternity inclusion was set to 10,000, 1000, 100, and 10, respectively. This reveals that simultaneous surveys of SNP and STR markers included in the ForenSeq DNA Signature Prep Kit offer a promising method for NIPPT using MPS technology.

Modeling allelic analyte signals for aSTRs in NGS DNA profiles

We describe an adaption of Bright et al.'s work modeling peak height variability in CE-DNA profiles to the modeling of allelic aSTR (autosomal short tandem repeats) read counts from NGS-DNA profiles, specifically for profiles generated from the ForenSeq™ DNA Signature Prep Kit, DNA Primer Mix B. Bright et al.'s model consists of three key components within the estimation of total allelic product-template, locus-specific amplification efficiencies, and degradation. In this work, we investigated the two mass parameters-template and locus-specific amplification efficiencies-and used MLE (maximum likelihood estimation) and MCMC (Markov chain Monte Carlo) methods to obtain point estimates to calculate the total allelic product. The expected read counts for alleles were then calculated after proportioning some of the expected stutter product from the total allelic product. Due to preferential amplicon selection introduced by the sample purification beads, degradation is difficult to model from the aSTR outputs alone. Improved modeling of the locus-specific amplification efficiencies may mask the effects of degradation. Whilst this model could be improved by introducing locus specific variances in addition to locus specific priors, our results demonstrate the suitability of adapting Bright et al.'s allele peak height model for NGS-DNA profiles. This model could be incorporated into continuous probabilistic interpretation approaches for mixed DNA profiles.

Developmental validation of the MGIEasy Signature Identification Library Prep Kit, an all-in-one multiplex system for forensic applications

Analyzing genetic markers in nuclear and mitochondrial genomes is helpful in various forensic applications, such as individual identifications and kinship analyses. However, most commercial kits detect these markers separately, which is time-consuming, laborious, and more error-prone (mislabelling, contamination, ...). The MGIEasy Signature Identification Library Prep Kit (hereinafter "MGIEasy identification system"; MGI Tech, Shenzhen, China) has been designed to provide a simple, fast, and robust way to detect appropriate markers in one multiplex PCR reaction: 52 autosomal STRs, 27 X-chromosomal STRs, 48 Y-chromosomal STRs, 145 identity-informative SNPs, 53 ancestry-informative SNPs, 29 phenotype-informative SNPs, and the hypervariable regions of mitochondrial DNA (mtDNA). Here, we validated the performance of MGIEasy identification system following the guidelines of the Scientific Working Group on DNA Analysis Methods (SWGDAM), assessing species specificity, sensitivity, mixture identification, stability under non-optimal conditions (degraded samples, inhibitor contamination, and various substrates), repeatability, and concordance. Libraries prepared using MGIEasy identification system were sequenced on a MGISEQ-2000 instrument (MGI Tech). MGIEasy-derived STR, SNP, and mtDNA genotypes were highly concordant with CE-based STR genotypes (99.79%), MiSeq FGx-based SNP genotypes (99.78%), and Sanger-based mtDNA genotypes (100%), respectively. This system was strongly human-specific, resistant to four common PCR inhibitors, and reliably amplified both low quantities of DNA (as low as 0.125 ng) and degraded DNA (~ 150 nt). Most of the unique alleles from the minor contributor were detected in 1:10 male-female and male-male mixtures; some minor Y-STR alleles were even detected in 1:1000 male-female mixtures. MGIEasy also successfully directly amplified markers from blood stains on FTA cards, filter papers, and swabs. Thus, our results demonstrated that MGIEasy identification system was suitable for use in forensic analyses due to its robust and reliable performance on samples of varying quality and quantity.

Massively parallel sequencing of STRs using a 29-plex panel reveals stutter sequence characteristics

Massively parallel sequencing of forensic STRs simultaneously provides length-based genotypes and core repeat sequences as well as flanking sequence variations. Here, we report primer sequences and concentrations of a next-generation sequencing (NGS)-based in-house panel covering 28 autosomal STR loci (CSF1PO, D1GATA113, D1S1627, D1S1656, D1S1677, D2S441, D2S1776, D3S3053, D5S818, D6S474, D6S1017, D6S1043, D8S1179, D9S2157, D10S1435, D11S4463, D13S317, D14S1434, D16S539, D18S51, D18S853, D20S482, D20S1082, D22S1045, FGA, TH01, TPOX, and vWA) and the sex determinant locus Amelogenin. Preliminary evaluation experiments showed that the panel yielded intralocus- and interlocus-balanced sequencing data with a sensitivity as low as 62.5 pg input DNA. A total of 203 individuals from Yunnan Bai population were sequenced with this panel. Comparative forensic genetic analyses showed that sequence-based matching probability of this 29-plex panel reached 2.37 × 10^-29 , which was 23 times lower than the length-based data. Compound stutter sequences of eight STRs were compared with parental alleles. For seven loci, repeat motif insertions or deletions occurred in the longest uninterrupted repeat sequences (LUS). However, LUS and non-LUS stutters co-existed in the locus D6S474 with different sequencing depth ratios. These results supplemented our current knowledge of forensic STR stutters, and provided a sound basis for DNA mixture deconvolution.