Mitochondrial genome sequencing with short overlapping amplicons on MiSeq FGx system

Forensic efficiency evaluation of a mtDNA whole genome sequencing system constructed with long fragment amplification strategy on DNA nanoball sequencing platform

Mitochondrial DNA (mtDNA) is an important genetic marker for degraded biological sample identification, maternal pedigree tracing, and population genetic structure study owing to its characteristics of high copy number, anti-degradable ring structure, and maternal inheritance. Whole mtDNA genome sequencing is an optimal method for the analysis of mtDNA polymorphism and heterogeneity because it allows for the comprehensive use of maternal genetic information. However, because of lacking quantitative evaluations for sequencing data, the scientific interpretation standards for mtDNA sequencing results of the previously used sequencing systems are often different, and false positive or false negative results are prone to occur when faced with the interference of nuclear genomic DNA, or the heterogeneities of mtDNA sequence and structure. In this study, we evaluated a novel mtDNA whole genome sequencing system using long fragment amplification strategy on the DNA nanoball (DNB) sequencing platform. This system demonstrated high sequencing quality and specific mtDNA sequencing efficiencies on positive control DNA and FTA bloodstain samples, as the average Q20 and Q30 values of the corresponding samples were 97.17 % and 91.93 %; 97.37 % and 92.48 %, respectively. The mean mapping percentages for the reference sequences of whole genome DNA (wgDNA), mtDNA, and nuclear genomic DNA (ngDNA) in the corresponding samples were 99.98 %, 99.97 %, 0.03 %, and 99.91 %, 99.40 %, 0.60 %; respectively. The average error calling rates for the bases A, C, G, and T of the whole mtDNA genome were 0.2519 %, 0.2550 %, 0.2906 %; and 0.2392 %, respectively. The efficacy of heteroplasmy identification was assessed using a set of theoretical sites with predetermined rates. These sites were created by combining the samples with known mtDNA haplotypes in certain proportions. The absolute errors between observed and theoretical heteroplasmy values were 89.59 %, 74.68 %, 50.20 %, 12.65 %, 8.31 %, and 4.85 %, while the theoretical heteroplasmy values were 5 %, 10 %, 20 %, 80 %, 90 %, and 95 %, respectively. The absolute error exhibited relative stability when the mtDNA sequencing depth exceeded 500×. Furthermore, the system sequencing efficiency was also confirmed among different kinds of samples, and these samples included natural samples (e.g., peripheral blood samples preserved on FTA cards for 2 and 11 years, and on filter paper for 6 and 9 years), degraded samples, sensitivity samples, samples derived from various bodily fluids, and maternal pedigree samples. In summary, the whole mtDNA genome sequencing system used for forensic identification demonstrated high performance in analyzing mtDNA sequence information, and showed significant prospects for forensic application and maternal genetic research.

Assessment of ForenSeq mtDNA Whole Genome Kit for forensic application

Mitochondrial DNA (mtDNA) is an indispensable genetic marker in forensic genetics. The emergence and development of massively parallel sequencing (MPS) makes it possible to obtain complete mitochondrial genome sequences more quickly and accurately. The study evaluated the advantages and limitations of the ForenSeq mtDNA Whole Genome Kit in the practical application of forensic genetics by detecting human genomic DNA standards and thirty-three case samples. We used control DNA with different amount to determine sensitivity of the assay. Even when the input DNA is as low as 2.5 pg, most of the mitochondrial genome sequences could still be covered. For the detection of buccal swabs and aged case samples (bloodstains, bones, teeth), most samples could achieve complete coverage of mitochondrial genome. However, when ancient samples and hair samples without hair follicles were sequenced by the kit, it failed to obtain sequence information. In general, the ForenSeq mtDNA Whole Genome Kit has certain applicability to forensic low template and degradation samples, and these results provide the data basis for subsequent forensic applications of the assay. The overall detection process and subsequent analysis are easy to standardize, and it has certain application potential in forensic cases.

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.

Title: Skeletal evidence of the ethnic cleansing actions in the Free City of Danzig (1939-1942) based on the KL Stutthof victims analysis

In the early days of World War II, many of the prominent and influential people of Polish nationality from the Free City of Danzig were arrested by the Germans and sent to the nearby concentration camp KL Stutthof. Nearly a hundred of them died within the next seven months upon their arrival, and were buried in a clandestine mass grave in a nearby forest. However, the exact nature of their death is unknown, as it is unclear what the attitude of the aggressors was toward the victims. We do not know whether there was only one executioner or there were several assassins, nor if the killing methodology was consistent with the other state-regulated executions. The studied material represents the commingled remains of a minimum thirty-four people, possibly all male, aged from under eighteen to over sixty at the time of death. Perimortem traumatic lesions are shown mainly on the skull bones. We asked whether the perimortem trauma lesions visible on the victims' skeletons could be informative on the cause and manner of their death. Our results show the prevalence of the perimortem trauma inflicted by a blunt object are on the parietal bones above the Hat Brim Line (HBL), which is commonly associated with a violent attack. The gunshot trauma was usually localized on the occipital bone or posterior parietal, which could indicate a shot to the back of the head, and this was commonly encountered during executions. No signs of defensive injuries can be explained either by restraining of the hands or by a surprise attack. The abundance and variability of the trauma type can be evident of multiple assailants. Moreover, the multiple impact points detected on several crania prove unnecessary overkill and brutality, which reflects the personal attitudes of the executioners towards the victims.