Massively parallel sequencing of the entire control region and targeted coding region SNPs of degraded mtDNA using a simplified library preparation method

Comparative analysis of SNaPshot and massively parallel sequencing for body fluid-specific DNA methylation markers

The identification of tissue-specific differentially methylated regions has significantly contributed to the field of forensic genetics, particularly in body fluid identification crucial for linking evidence to crimes. Among the various approaches to analyzing DNA methylation, the SNaPshot assay has been popularly studied in numerous researches. However, there is a growing interest in exploring alternative methods such as the use of massively parallel sequencing (MPS), which can process a large number of samples simultaneously. This study compares SNaPshot and MPS multiplex assays using nine cytosine-phosphate-guanine markers for body fluid identification. As a result of analyzing 112 samples, including blood, saliva, vaginal fluid, menstrual blood, and semen, both methods demonstrated high sensitivity and specificity, indicating their reliability in forensic investigations. A total of 92.0% samples were correctly identified by both methods. Although both methods accurately identified all blood, saliva, and semen samples, some vaginal fluid samples showed unexpected methylation signals at nontarget loci in addition to the target loci. In the case of menstrual blood samples, due to their complexity, independent typing criteria were applied, and successful menstrual blood typing was possible, whereas a few samples showed profiles similar to vaginal fluid. The MPS method worked better in vaginal fluid samples, and the SNaPshot method performed better in menstrual blood samples. This study offers valuable insights into body fluid identification based on the characteristics of the SNaPshot and MPS methods, which may help in more efficient forensic applications.

A Customized Human Mitochondrial DNA Database (hMITO DB v1.0) for Rapid Sequence Analysis, Haplotyping and Geo-Mapping

The field of mitochondrial genomics has advanced rapidly and has revolutionized disciplines such as molecular anthropology, population genetics, and medical genetics/oncogenetics. However, mtDNA next-generation sequencing (NGS) analysis for matrilineal haplotyping and phylogeographic inference remains hindered by the lack of a consolidated mitogenome database and an efficient bioinformatics pipeline. To address this, we developed a customized human mitogenome database (hMITO DB) embedded in a CLC Genomics workflow for read mapping, variant analysis, haplotyping, and geo-mapping. The database was constructed from 4286 mitogenomes. The macro-haplogroup (A to Z) distribution and representative phylogenetic tree were found to be consistent with published literature. The hMITO DB automated workflow was tested using mtDNA-NGS sequences derived from Pap smears and cervical cancer cell lines. The auto-generated read mapping, variants track, and table of haplotypes and geo-origins were completed in 15 min for 47 samples. The mtDNA workflow proved to be a rapid, efficient, and accurate means of sequence analysis for translational mitogenomics.

Internal validation and improvement of mitochondrial genome sequencing using the Precision ID mtDNA Whole Genome Panel

With the recent advances in next-generation sequencing (NGS), mitochondrial whole-genome sequencing has begun to be applied to the field of the forensic biology as an alternative to the traditional Sanger-type sequencing (STS). However, experimental workflows, commercial solutions, and output data analysis must be strictly validated before being implemented into the forensic laboratory. In this study, we performed an internal validation for an NGS-based typing of the entire mitochondrial genome using the Precision ID mtDNA Whole Genome Panel (Thermo Fisher Scientific) on the Ion S5 sequencer (Thermo Fisher Scientific). Concordance, repeatability, reproducibility, sensitivity, and heteroplasmy detection analyses were assessed using the 2800 M and 9947A standard control DNA as well as typical casework specimens, and results were compared with conventional Sanger sequencing and another NGS sequencer in a different laboratory. We discuss the strengths and limitations of this approach, highlighting some issues regarding noise thresholds and heteroplasmy detection, and suggesting solutions to mitigate these effects and improve overall data interpretation. Results confirmed that the Precision ID Whole mtDNA Genome Panel is highly reproducible and sensitive, yielding useful full mitochondrial DNA sequences also from challenging DNA specimens, thus providing further support for its use in forensic practice.

Massively parallel sequencing of 25 autosomal STRs including SE33 in four population groups for forensic applications

The introduction of massively parallel sequencing (MPS) in forensic investigation enables sequence-based large-scale multiplexing beyond size-based analysis using capillary electrophoresis (CE). For the practical application of MPS to forensic casework, many population studies have provided sequence data for autosomal short tandem repeats (STRs). However, SE33, a highly polymorphic STR marker, has little sequence-based data because of difficulties in analysis. In this study, 25 autosomal STRs were analyzed, including SE33, using an in-house MPS panel for 350 samples from four populations (African-American, Caucasian, Hispanic, and Korean). The barcoded MPS library was generated using a two-step PCR method and sequenced using a MiSeq System. As a result, 99.88% genotype concordance was obtained between length- and sequence-based analyses. In SE33, the most discordances (eight samples, 0.08%) were observed because of the 4 bp deletion between the CE and MPS primer binding sites. Compared with the length-based CE method, the number of alleles increased from 332 to 725 (2.18-fold) for 25 autosomal STRs in the sequence-based MPS method. Notably, additional 129 unique alleles, a 4.15-fold increase, were detected in SE33 by identifying sequence variations. This population data set provides sequence variations and sequence-based allele frequencies for 25 autosomal STRs.

Development and Verification of an Economical Method of Custom Target Library Construction

Although technological advances have greatly reduced the cost of DNA sequencing, sample preparation time and reagent costs remain the limiting factors for many studies. Based on low-cost targeted amplification, we developed an economical method for custom target library construction based on DNA nanoball (DNB) technology and two-step polymerase chain reaction (PCR). Here, we refer to this method as the two-step PCR, which was compared to traditional multiplex PCR methods in three aspects, data quality, efficiency, and specificity to humans. The results confirmed that two-step PCR reduces to finishing 128 sequencing libraries in only 2 h 24 min 59 s of the total PCR time and at a data utilization rate of 0.44 at a cost of approximately $1.70 per sample for targeted sequencing via the two-step PCR. The replacement of traditional multiplex PCR methods with this strategy makes the sample preparation process before sequencing relatively more cost-effective and further reduces the cost of next-generation sequencing (NGS). This method may also be free from the interference of other species and the limitations of sample type and DNA content. These findings reveal possibilities for broad applications of this approach in forensic research.

A review of bioinformatic methods for forensic DNA analyses

Short tandem repeats, single nucleotide polymorphisms, and whole mitochondrial analyses are three classes of markers which will play an important role in the future of forensic DNA typing. The arrival of massively parallel sequencing platforms in forensic science reveals new information such as insights into the complexity and variability of the markers that were previously unseen, along with amounts of data too immense for analyses by manual means. Along with the sequencing chemistries employed, bioinformatic methods are required to process and interpret this new and extensive data. As more is learnt about the use of these new technologies for forensic applications, development and standardization of efficient, favourable tools for each stage of data processing is being carried out, and faster, more accurate methods that improve on the original approaches have been developed. As forensic laboratories search for the optimal pipeline of tools, sequencer manufacturers have incorporated pipelines into sequencer software to make analyses convenient. This review explores the current state of bioinformatic methods and tools used for the analyses of forensic markers sequenced on the massively parallel sequencing (MPS) platforms currently most widely used.

Investigation into the sequence structure of 23 Y chromosomal STR loci using massively parallel sequencing

Next-generation sequencing (NGS) can produce massively parallel sequencing (MPS) data for many targeted regions with a high depth of coverage, suggesting its successful application to the amplicons of forensic genetic markers. In the present study, we evaluated the practical utility of MPS in Y-chromosome short tandem repeat (Y-STR) analysis using a multiplex polymerase chain reaction (PCR) system. The multiplex PCR system simultaneously amplified 24 Y-chromosomal markers, including the PowerPlex^® Y23 loci (DYS19, DYS385ab, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, DYS448, DYS456, DYS458, DYS481, DYS533, DYS549, DYS570, DYS576, DYS635, DYS643, and YGATAH4) and the M175 marker with the small-sized amplicons ranging from 85 to 253bp. The barcoded libraries for the amplicons of the 24 Y-chromosomal markers were produced using a simplified PCR-based library preparation method and successfully sequenced using MPS on a MiSeq^® System with samples from 250 unrelated Korean males. The genotyping concordance between MPS and the capillary electrophoresis (CE) method, as well as the sequence structure of the 23 Y-STRs, were investigated. Three samples exhibited discordance between the MPS and CE results at DYS385, DYS439, and DYS576. There were 12 Y-STR loci that showed sequence variations in the alleles by a fragment size determination, and the most varied alleles occurred in DYS389II with a different sequence structure in the repeat region. The largest increase in gene diversity between the CE and MPS results was in DYS437 at +34.41%. Single nucleotide polymorphisms (SNPs), insertions, and deletions (indels) were observed in the flanking regions of DYS481, DYS576, and DYS385, respectively. Stutter and noise ratios of the 23 Y-STRs using the developed MPS system were also investigated. Based on these results, the MPS analysis system used in this study could facilitate the investigation into the sequences of the 23 Y-STRs in forensic genetics laboratories.