Massively parallel sequencing of human skeletal remains in Vietnam using the precision ID mtDNA control region panel on the Ion S5™ system

Simultaneous sequencing of 102 Y-STRs on Ion Torrent ™ GeneStudio ™ S5 System

The Precision ID NGS System from Thermo Fisher Scientific is a mainstream next-generation sequencing (NGS) platform used in forensic laboratories to detect almost all commonly used forensic markers, except for Y-chromosomal short tandem repeats (Y-STRs). This study aimed to: 1) develop a Y-STR panel compatible with the automatic workflow of the NGS system using Ion AmpliSeq Technology, 2) evaluate the panel performance following the SWGDAM guidelines, and 3) explore the possibility of using a combination workflow to detect autosomal STRs and Y-STRs (AY-STR NGS workflow). The GrandFiler Y-STR Panel was successfully designed using the 'separating' and 'merging' strategies, including 102 Y-STRs and Amelogenin with an average amplicon length of 133 bp. It is a mega Y-STR multiplex system in which up to 16 samples can be sequenced simultaneously on an Ion 530 ™ Chip. Developmental validation studies of the performance of the NGS platform, species specificity, reproducibility, concordance, sensitivity, degraded samples, case-type samples, and mixtures were conducted to unequivocally determine whether the GrandFiler Y-STR Panel is suitable for real scenarios. The newly developed Y-STR panel showed compelling run metrics and NGS performance, including 92.47% bases with ≥ Q20, 91.80% effective reads, 2106 × depth of coverage (DoC), and 97.09% inter-locus balance. Additionally, it showed high specificity for human males and 99.40% methodological and bioinformatical concordance, generated complete profiles at ≥ 0.1 ng input DNA, and recovered more genetic information from severely degraded and diverse case samples. Although the outcome when used on mixtures was not as expected, more genetic information was obtained compared to that from capillary electrophoresis (CE) methods. The AY-STR NGS workflow was established by combining the GrandFiler Y-STR Panel with the Precision ID GlobalFiler ™ NGS STR Panel v2 at a 2:1 concentration ratio. The combination workflow on NGS performance, reproducibility, concordance, and sensitivity was as stable as the single Y-STR NGS workflow, providing more options for forensic scientists when dealing with different case scenarios. Overall, the GrandFiler Y-STR Panel was confirmed as the first to effectively detect a large number of Y-STR markers on the Precision ID NGS System, which is compatible with 51 Y-STRs in commercial CE kits and 51 Y-STRs in commercial NGS kits and the STRBase. The panel is as robust, reliable, and sensitive as current CE/NGS kits, and is suitable for solving real cases, especially for severely degraded samples (degradation index > 10).

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.

Comparison of Two Variant Analysis Programs for Next-Generation Sequencing Data of Whole Mitochondrial Genome

BACKGROUND

With advance of next-generation sequencing (NGS) techniques, the need for mitochondrial DNA analysis is increasing not only in the forensic area, but also in medical fields.

METHODS

Two commercial programs, Converge Software (CS) and Torrent Variant Caller for variant calling of NGS data, were compared with a considerable amount of sequence data of 50 samples with a homogeneous ethnicity.

RESULTS

About 2,300 variants were identified and the two programs showed about 90% of consistency. CS, a dedicated analysis program for mitochondrial DNA, showed some advantages for forensic use. By additional visual inspection, several causes of discrepancy in variant calling results were identified. Application of different notation rules for mitochondrial sequence and the minor allele frequency close to detection threshold were the two most significant reasons.

CONCLUSION

With prospective improvement of each program, researchers and practitioners should be aware of characteristics of the analysis program they use and prepare their own strategies to determine variants.

Mitochondrial DNA control region typing from highly degraded skeletal remains by single-multiplex next-generation sequencing

Poor nuclear DNA preservation from highly degraded skeletal remains is the most limiting factor for the genetic identification of individuals. Mitochondrial DNA (mtDNA) typing, and especially of the control region (CR), using next-generation sequencing (NGS), enables retrieval of valuable genetic information in forensic contexts where highly degraded human skeletal remains are the only source of genetic material. Currently, NGS commercial kits can type all mtDNA-CR in fewer steps than the conventional Sanger technique. The PowerSeq CRM Nested System kit (Promega Corporation) employs a nested multiplex-polymerase chain reaction (PCR) strategy to amplify and index all mtDNA-CR in a single reaction. Our study analyzes the success of mtDNA-CR typing of highly degraded human skeletons using the PowerSeq CRM Nested System kit. We used samples from 41 individuals from different time periods to test three protocols (M1, M2, and M3) based on modifications of PCR conditions. To analyze the detected variants, two bioinformatic procedures were compared: an in-house pipeline and the GeneMarker HTS software. The results showed that many samples were not analyzed when the standard protocol (M1) was used. In contrast, the M3 protocol, which includes 35 PCR cycles and longer denaturation and extension steps, successfully recovered the mtDNA-CR from highly degraded skeletal samples. Mixed base profiles and the percentage of damaged reads were both indicators of possible contamination and can provide better results if used together. Furthermore, our freely available in-house pipeline can provide variants concordant with the forensic software.

The first maternal genetic study of hunter-gatherers from Vietnam

The current limitation of ancient DNA data from Vietnam led to the controversy surrounding the prehistory of people in this region. The combination of high heat and humidity damaged ancient bones that challenged the study of human evolution, especially when using DNA as study materials. So far, only 4 k years of history have been recorded despite the 65 k years of history of anatomically modern human occupations in Vietnam. Here we report, to our knowledge, the oldest mitogenomes of two hunter-gatherers from Vietnam. We extracted DNA from the femurs of two individuals aged 6.2 k cal BP from the Con Co Ngua (CCN) site in Thanh Hoa, Vietnam. This archeological site is the largest cemetery of the hunter-gatherer population in Southeast Asia (SEA) that was discovered, but their genetics have not been explored until the present. We indicated that the CCN haplotype belongs to a rare haplogroup that was not detected in any present-day Vietnamese individuals. Further matrilineal analysis on CCN mitogenomes showed a close relationship with ancient farmers and present-day populations in SEA. The mitogenomes of hunter-gatherers from Vietnam debate the "two layers" model of peopling history in SEA and provide an alternative solution for studying challenging ancient human samples from Vietnam.

DNA and protein analyses of hair in forensic genetics

Hair is one of the most common pieces of biological evidence found at a crime scene and plays an essential role in forensic investigation. Hairs, especially non-follicular hairs, are usually found at various crime scenes, either by natural shedding or by forcible shedding. However, the genetic material in hairs is usually highly degraded, which makes forensic analysis difficult. As a result, the value of hair has not been fully exploited in forensic investigations and trials. In recent years, with advances in molecular biology, forensic analysis of hair has achieved remarkable strides and provided crucial clues in numerous cases. This article reviews recent developments in DNA and protein analysis of hair and attempts to provide a comprehensive solution to improve forensic hair analysis.

Genetic Background of Kirgiz Ethnic Group From Northwest China Revealed by Mitochondrial DNA Control Region Sequences on Massively Parallel Sequencing

The mitochondrial DNA (mtDNA) has been used to trace population evolution and apply to forensic identification due to the characteristics including lack of recombination, higher copy number and matrilineal inheritance comparing with nuclear genome DNA. In this study, mtDNA control region sequences of 91 Kirgiz individuals from the Northwest region of China were sequenced to identify genetic polymorphisms and gain insight into the genetic background of the Kirgiz ethnic group. MtDNA control region sequences of Kirgiz individuals presented relatively high genetic polymorphisms. The 1,122 bp sequences of mtDNA control region could differ among unrelated Kirgiz individuals, which suggested the mtDNA control region sequences have a good maternal pedigree tracing capability among different Kirgiz individuals. The neutrality test, mismatch distribution, Bayesian phylogenetic inference, Bayesian skyline analysis, and the median network analyses showed that the Kirgiz group might occurred population expansion, and the expansion could be observed at about ∼53.41 kilo years ago (kya) when ancestries of modern humans began to thrive in Eurasia. The pairwise population comparisons, principal component analyses, and median network analyses were performed based on haplogroup frequencies or mtDNA control region sequences of 5,886 individuals from the Kirgiz group and the 48 reference populations all over the world. And the most homologous haplotypes were found between Kirgiz individuals and the East Asian individuals, which indicated that the Kirgiz group might have gene exchanges with the East Asian populations.

Development and validation of a novel 133-plex forensic STR panel (52 STRs and 81 Y-STRs) using single-end 400 bp massive parallel sequencing

Short tandem repeats (STRs) are the preferred genetic markers in forensic DNA analysis, routinely measured by capillary electrophoresis (CE) method based on the fragment length features. While, the massive parallel sequencing (MPS) technology could simultaneously target a large number of intriguing forensic STRs, bypassing the intrinsic limitations of amplicon size separation and accessible fluorophores in CE, which is efficient and promising for enabling the identification of forensic biological evidence. Here, we developed a novel MPS-based Forensic Analysis System Multiplecues SetB Kit of 133-plex forensic STR markers (52 STRs and 81 Y-STRs) and one Y-InDel (M175) based on multiplex PCR and single-end 400 bp sequencing strategy. This panel was subjected to developmental validation studies according to the SWGDAM Validation Guidelines. Approximately 2185 MPS-based reactions using 6 human DNA standards and 8 male donors were conducted for substrate studies (filter paper, gauze, cotton swab, four different types of FTA cards, peripheral venous blood, saliva, and exfoliated cells), sensitivity studies (from 2 ng down to 0.0625 ng), mixture studies (two-person DNA mixtures), PCR inhibitor studies (seven commonly encountered PCR inhibitors), species specificity studies (11 non-human species), and repeatability studies. Results of concordance studies (413 Han males and 6 human DNA standards) generated by STRait Razor and in-house Python scripts indicated 99.98% concordance rate in STR calling relative to CE for STRs between 41,900 genotypes at 100 STR markers. Moreover, the limitations of present studies, the nomenclature rules and forensic MPS applications were also described. In conclusion, the validation studies based on ~ 2200 MPS-based and ~ 2500 CE-based DNA profiles demonstrated that the novel MPS-based panel meets forensic DNA quality assurance guidelines with robust, reliable, and reproducible performance on samples of various quantities and qualities, and the STR nomenclature rules should be further regulated to integrate the inconformity between MPS-based and CE-based methods.