Forensic characterization of Brazilian regional populations through massive parallel sequencing of 124 SNPs included in HID ion Ampliseq Identity Panel

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.

Sequenced-based Rwanda population provides insights into demographic history

Rwanda is known as the heart of Africa, reflecting the history of the world. Colonization and genocide have led to Rwanda's existing genetic structure. Herein, we used massively parallel sequencing to analyze 296 loci in 185 Rwandans and constructed a database for Rwandan forensic data for the first time. We found the following results: First, forensic parameters demonstrated that all loci were highly informative and could be used for forensic identification and paternity tests in Rwandans. Second, we found that the differences in genetic background between Rwandans and other African populations were similar but slight, as indicated by the massively parallel sequencing panel. Rwandans belonged to the African population and were inseparable from populations from neighboring countries. Also, Rwandans were closer to the European and American populations because of colonization, war, and other reasons. There was no scientific basis for racial classification established by colonization. Further research still needs to be carried out on more loci and larger Rwandan samples.

Development and validation of a custom panel including 114 InDels using massively parallel sequencing for forensic application

Insertion/deletion polymorphisms (InDels) have particular characteristics, such as a relatively low mutation rate, small amplicon size, and no stutter artifacts when genotyped via the capillary electrophoresis platform. It would be an important complementary tool for individual identification and certain kinship analyses. At present, massively parallel sequencing (MPS) has shown excellent application value in forensic studies. Therefore, in this study, we developed a custom MPS InDel panel that contains 114 InDels [77 autosomal InDels (A-InDels), 32 X-chromosomal InDels (X-InDels), and 5 Y-chromosomal InDels) based on previous studies. To assess this panel's performance, several validation experiments were performed, including sensitivity, inhibitor, degraded DNA testing, species specificity, concordance, repeatability, case-type samples, and population studies. The results showed that the lowest DNA input was 0.25 ng. All genotypes were obtained in the presence of 80 ng/µL humic acid, 2000 µmol/L calcium, 3000 µmol/L EDTA and indigo. In degraded DNA testing, 90% of loci could be detected for 16-day-old formalin-fixed hearts. In addition, this panel has good species specificity. The values of combined power of discrimination and the combined power of exclusion for 77 A-InDels were 1-3.9951 × 10-32 and 1-4.2956 × 10-7 , respectively. The combined mean exclusion chance for 32 X-InDels was 0.99999 in trios and 0.99904 in duos. The validation results indicate that this newly developed MPS multiplex system is a robust tool for forensic applications.

Forensic genetic analysis of single-nucleotide polymorphisms and microhaplotypes in Koreans through next-generation sequencing using precision ID identity panel

BACKGROUND

Forensic DNA analysis has seen remarkable advancements with the advent of Next Generation Sequencing (NGS). In particular, NGS analysis of single nucleotide polymorphisms (SNPs) offers significant advantages in the analysis of challenging samples compared to conventional STR analysis.

OBJECTIVE

This study aimed to investigate the SNPs of the Precision ID Identity Panel, a commercially available NGS panel for personal identification, by generating genetic profiles of 298 Koreans and comparing them with other global populations.

METHODS

A total of 124 SNPs, including 90 autosomal and 34 Y-SNPs, were analyzed using the Precision ID Identity Panel, and forensic parameters, microhaplotypes, and population differences were investigated.

RESULTS

The NGS data were successfully obtained from 298 Koreans. The analysis of forensic parameters exhibited a low combined match probability of 1.532 × 10- 34, which is comparable to that obtained from commonly used STR analysis. Additionally, the microhaplotype analysis revealed that the use of 16 microhaplotypes provided higher discriminatory power compared to single target SNPs. Furthermore, the adoption of microhaplotype data resulted in an increase of over 20% in expected heterozygosity at five loci. Inter-population analysis showed a close genetic relationship between Koreans and individuals from China and Myanmar in East and Southeast Asia, which are geographically adjacent to Korea.

CONCLUSIONS

The results of this study show that the Precision ID Identity panel can be a useful alternative where traditional STR typing is not feasible. Also, the data from our study will be useful as a reference for Koreans in forensic investigations and the prosecution of criminal justice.

MPKin-YSTR: Interpretation of Y chromosome STR haplotypes for missing persons cases

Y chromosome Short Tandem Repeat (STR) haplotypes have been used in assisting forensic investigations primarily for identification and male lineage determination. The current SWGDAM interpretation guidelines for Y-STR typing provide helpful guidance on those purposes but do not address the issue of kinship analysis with Y-STR haplotypes. Because of the high mutation rate of Y-STRs, there are complex missing person cases in which inconsistent Y-STR haplotypes between true paternal lineage relatives will arise and cases with two or more male references in the same lineage and yet differ in their haplotypes. Therefore, more useful methods are needed for interpreting the Y-STR haplotype data. Computational methods and interpretation guidelines have been developed specifically addressing this issue, either using a mismatch-based counting method or a pedigree likelihood ratio method. In this study, a software program, MPKin-YSTR, was developed by implementing those more sophisticated methods. This software should be able to improve the interpretation of complex cases with Y-STR haplotype evidence. Thus, more biological evidence will be interpreted, which in turn will result in more investigation leads to help solve crimes.

Genetic investigation of 124 SNPs in a Myanmar population using the Precision ID Identity Panel and the Illumina MiSeq

Single nucleotide polymorphisms (SNPs) have become popular in forensic genetics as an alternative to short tandem repeats (STRs). The Precision ID Identity Panel (Thermo Fisher Scientific), consisting of 90 autosomal SNPs and 34 Y-chromosomal SNPs, enabled human identification studies on global populations through next-generation sequencing (NGS). However, most previous studies on the panel have used the Ion Torrent platform, and there are few reports on the Southeast Asian population. Here, a total of 96 unrelated males from Myanmar (Yangon) were analyzed with the Precision ID Identity Panel on a MiSeq (Illumina) using an in-house TruSeq compatible universal adapter and a custom variant caller, Visual SNP. The sequencing performance evaluated by locus balance and heterozygote balance was comparable to that of the Ion Torrent platform. For 90 autosomal SNPs, the combined match probability (CMP) was 6.994 × 10^-34, lower than that of 22 PowerPlex Fusion autosomal STRs (3.130 × 10^-26). For 34 Y-SNPs, 14 Y-haplogroups (mostly O2 and O1b) were observed. We found 51 cryptic variations (42 haplotypes) around target SNPs, of which haplotypes corresponding to 33 autosomal SNPs decreased CMP. Interpopulation analysis revealed that the Myanmar population is genetically closer to the East and Southeast Asian populations. In conclusion, the Precision ID Identity Panel can be successfully analyzed on the Illumina MiSeq and provides high discrimination power for human identification in the Myanmar population. This study broadened the accessibility of the NGS-based SNP panel by expanding the available NGS platforms and adopting a robust NGS data analysis tool.

Ancestry resolution of South Brazilians by forensic 165 ancestry-informative SNPs panel

Forensic DNA phenotyping (FDP) includes biogeographic ancestry (BGA) inference and externally visible characteristics (EVCs) prediction directly from an evidential DNA sample as alternatives to provide valuable intelligence when conventional DNA profiling fails to achieve identification. In this context, the application of Massively Parallel Sequencing (MPS) methodologies, which enables simultaneous typing of multiple samples and hundreds of forensic markers, has been gradually implemented in forensic genetic casework. The Precision ID Ancestry Panel (Thermo Fisher Scientific, Waltham, USA) is a forensic multiplex assay consisting of 165 autosomal SNPs designed to provide biogeographic ancestry information. In this work, a sample of 250 individuals from Rio Grande do Sul (RS) State, southern Brazil, apportioned into four main population groups (African-, European-, Amerindian-, and Admixed-derived Gauchos), was evaluated with this panel, to assess the feasibility of this approach in a highly heterogeneous population. Forensic descriptive parameters estimated for each population group revealed that this panel has enough polymorphic and informative SNPs to be used as a supplementary instrument in forensic individual identification and kinship testing regardless of ethnicity. No statistically significant deviation from Hardy-Weinberg equilibrium was observed after Bonferroni correction. However, seven loci pairs displayed linkage disequilibrium in pairwise LD testing (p < 3.70 × 10^-6). Interpopulation comparisons by F_ST analysis, MDS plot, and STRUCTURE analysis among the four RS population groups apart and along with 89 reference worldwide populations demonstrated that Admixed- and African-derived Gauchos present the highest levels of admixture and population stratification, whereas European- and Amerindian-derived exhibit a more homogeneous genetic conformation.

HIV-1 Drug Resistance Assay Using Ion Torrent Next Generation Sequencing and On-Instrument End-to-End Analysis Software

HIV-1 antiretroviral therapy management requires sequencing the protease, reverse transcriptase, and integrase portions of the HIV-1 pol gene. Most resistance testing is performed with Sanger sequencing, which has limited ability to detect minor variants. Next generation sequencing (NGS) platforms enable variant detection at frequencies as low as 1% allowing for earlier detection of resistance and modification of therapy. Implementation of NGS assays in the clinical laboratory is hindered by complicated assay design, cumbersome wet bench procedures, and the complexity of data analysis and bioinformatics. We developed a complete NGS protocol and companion analysis and reporting pipeline using AmpliSeq multiplex PCR, Ion Torrent S5 XL sequencing, and Stanford's HIVdb resistance algorithm. Implemented as a Torrent Suite software plugin, the pipeline runs automatically after sequencing. An optimum variant frequency threshold of 10% was determined by comparing Sanger sequences of archived samples from ViroSeq testing, resulting in a sensitivity of 98.2% and specificity of 99.0%. The majority (91%) of drug resistance mutations were detected by both Sanger and NGS, with 1.7% only by Sanger and 7.3% only by NGS. Variant calls were highly reproducible and there was no cross-reactivity to VZV, HBV, CMV, EBV, and HCV. The limit of detection was 500 copies/mL. The NGS assay performance was comparable to ViroSeq Sanger sequencing and has several advantages, including a publicly available end-to-end analysis and reporting plugin. The assay provides a straightforward path for implementation of NGS for HIV drug resistance testing in the laboratory setting without additional investment in bioinformatics infrastructure and resources.

A New Computational Deconvolution Algorithm for the Analysis of Forensic DNA Mixtures with SNP Markers

Single nucleotide polymorphisms (SNPs) support robust analysis on degraded DNA samples. However, the development of a systematic method to interpret the profiles derived from the mixtures is less studied, and it remains a challenge due to the bi-allelic nature of SNP markers. To improve the discriminating power of SNPs, this study explored bioinformatic strategies to analyze mixtures. Then, computer-generated mixtures were produced using real-world massively parallel sequencing (MPS) data from the single samples processed with the Precision ID Identity Panel. Moreover, the values of the frequency of major allele reads (F_MAR) were calculated and applied as key parameters to deconvolve the two-person mixtures and estimate mixture ratios. Four custom R language scripts (three for autosomes and one for Y chromosome) were designed with the K-means clustering method as a core algorithm. Finally, the method was validated with real-world mixtures. The results indicated that the deconvolution accuracy for evenly balanced mixtures was 100% or close to 100%, which was the same as the deconvolution accuracy of inferring the genotypes of the major contributor of unevenly balanced mixtures. Meanwhile, the accuracy of inferring the genotypes of the minor contributor decreased as its proportion in the mixture decreased. Moreover, the estimated mixture ratio was almost equal to the actual ratio between 1:1 and 1:6. The method proposed in this study provides a new paradigm for mixture interpretation, especially for inferring contributor profiles of evenly balanced mixtures and the major contributor profile of unevenly balanced mixtures.

Haplotype distribution in a forensic full mtDNA genome database of admixed Southern Brazilians and its association with self-declared ancestry and pigmentation traits

BACKGROUND

The advent of massively parallel sequencing (MPS) applications focused on the generation of forensic-quality full mitochondrial genome sequences led to a popularization of the technique on a global scale. However, the lack of forensic-graded population databases has refrained a wider adoption of full genome sequences as the industry standard, despite its better discrimination capacity of individual maternal lineages.

PURPOSE

This work describes a forensic-oriented full mtDNA genome database comprised of 480 samples from a Southern Brazilian population.

METHODS

A collection of mitochondrial sequences were obtained from low-pass, full genome DNA sequencing results. The complete sample set was evaluated regarding haplotype composition and distribution. Summary statistics and forensic parameters were calculated and are presented for the database, with detailed information concerning the impact of removing genetic information in the form of specific variants or increasingly larger genomic regions. Interpopulational analysis comparing haplotypical diversity in Brazilian and 26 worldwide populations was also performed. The association between mitochondrial genetic variability and phenotypic diversity was also evaluated in populations, with self-declared ancestry and three distinct phenotypic pigmentation traits (eyes, skin and hair colors) as parameters.

RESULTS

The presented database can be used to evaluate mitochondrial-related genetic evidence, providing LR values of up to 20,465 for unobserved haplotypes. Haplotype distribution in Southern Brazil seems to be different than the remaining of the country, with a larger contribution of maternal lines with European origin. Despite association can be found between lighter and darker phenotypes or self-declared ancestry and haplotype distribution, prediction models cannot be reliably proposed due to the admixed nature of the Brazilian population.

CONCLUSIONS

The proposed database provides a basis for statistical calculation and frequency estimation of full mitochondrial genomes, and can be part of an integrated, representative, national database comprising most of the genetic diversity of maternal lineages in the country.

Sequence variations, flanking region mutations, and allele frequency at 31 autosomal STRs in the central Indian population by next generation sequencing (NGS)

Capillary electrophoresis-based analysis does not reflect the exact allele number variation at the STR loci due to the non-availability of the data on sequence variation in the repeat region and the SNPs in flanking regions. Herein, this study reports the length-based and sequence-based allelic data of 138 central Indian individuals at 31 autosomal STR loci by NGS. The sequence data at each allele was compared to the reference hg19 sequence. The length-based allelic results were found in concordance with the CE-based results. 20 out of 31 autosomal STR loci showed an increase in the number of alleles by the presence of sequence variation and/or SNPs in the flanking regions. The highest gain in the heterozygosity and allele numbers was observed in D5S2800, D1S1656, D16S539, D5S818, and vWA. rs25768 (A/G) at D5S818 was found to be the most frequent SNP in the studied population. Allele no. 15 of D3S1358, allele no. 19 of D2S1338, and allele no. 22 of D12S391 showed 5 isoalleles each with the same size and with different intervening sequences. Length-based determination of the alleles showed Penta E to be the most useful marker in the central Indian population among 31 STRs studied; however, sequence-based analysis advocated D2S1338 to be the most useful marker in terms of various forensic parameters. Population genetics analysis showed a shared genetic ancestry of the studied population with other Indian populations. This first-ever study to the best of our knowledge on sequence-based STR analysis in the central Indian population is expected to prove the use of NGS in forensic case-work and in forensic DNA laboratories.

Forensic characterization of 124 SNPs in the central Indian population using precision ID Identity Panel through next-generation sequencing

With the advent of next-generation sequencing technology, SNP markers are being explored as a useful alternative to conventional capillary electrophoresis-based STR typing. Low mutation rate and short-sized amplicons are added advantages of SNP markers over the STRs. However, to achieve a sufficient level of discrimination among individuals, a higher number of SNPs need to be characterized simultaneously. Hence, the NGS technique is highly useful to analyze a sufficiently higher number of SNPs simultaneously. Though the technique is in its nascent stage, an attempt has been made to assess its usability in the central Indian population by analyzing 124 SNPs (90 autosomal and 34 Y-chromosome) in 95 individuals. Various quality parameters such as locus balance, locus strand balance, heterozygosity balance, and noise level showed a good quality sequence obtained from the Ion GeneStudio S5 instrument. Obtained frequency of SNP alleles ranged from 0.001 to 0.377 in autosomal SNPs. rs9951171 was found to be the most informative SNP in the studied population with the highest PD and lowest MP value. The cumulative MP of 90 SNPs was found to be 4.76698 × 10^-37. Analysis of 34 Y-chromosome SNPs reveals 11 unique haplogroups in 54 male samples with R1a1 as the most frequent haplogroup found in 22.22% of samples. Interpopulation comparison by FST analysis, PCA plot, and STRUCTURE analysis showed genetic stratification of the studied population suggesting the utility of SNP markers present in the Precision ID Identity Panel for forensic demands of the Indian population.

An Introductory Overview of Open-Source and Commercial Software Options for the Analysis of Forensic Sequencing Data

The top challenges of adopting new methods to forensic DNA analysis in routine laboratories are often the capital investment and the expertise required to implement and validate such methods locally. In the case of next-generation sequencing, in the last decade, several specifically forensic commercial options became available, offering reliable and validated solutions. Despite this, the readily available expertise to analyze, interpret and understand such data is still perceived to be lagging behind. This review gives an introductory overview for the forensic scientists who are at the beginning of their journey with implementing next-generation sequencing locally and because most in the field do not have a bioinformatics background may find it difficult to navigate the new terms and analysis options available. The currently available open-source and commercial software for forensic sequencing data analysis are summarized here to provide an accessible starting point for those fairly new to the forensic application of massively parallel sequencing.

Uniparental disomy of chromosome 21: A statistical approach and application in paternity tests

Considering the overall frequency of paternity investigation cases including mutational events, there is a real possibility that at least a fraction of all inconsistencies reported in paternity cases are caused not by polymerase slippage mutations, but to chromosomic abnormalities, as Uniparental Disomy (UPD). We report here the investigation of a trio paternity case (mother, child and alleged father), with observed inconsistencies that can alternatively be explained by occurrence of maternal uniparental isodisomy of chromosome 21 (miUPD21). A total of 350 short tandem repeat (STR) and single nucleotide polymorphism (SNP) markers were tested, statistically suggesting true biological linkage within the trio. Additionally, we propose miUPD21 explains, with significantly greater probability, the occurrence of detected inconsistencies, when compared to alternative hypothesis of multiple and simultaneous slippage mutations. Similar cases could have their statistical conclusions improved or even altered by including unusual chromosomal segregation patterns in the hypothesis formulation, as well as in mathematical calculations. Such reports of allelic inconsistencies being explained by chromosomal alterations are common in clinical genetics, and such situations might have impact on forensic investigation.

A novel forensic panel of 186-plex SNPs and 123-plex STR loci based on massively parallel sequencing

The MiSeq® FGX Forensic system and the HID-Ion AmpliSeq Panel were previously developed for massively parallel sequencing (MPS) for forensic casework. Among the three major sequencing platforms, BGISEQ-500^TM, which is based on multiple PCRs, is still lacking in forensics. Here, a novel forensic panel was constructed to detect 186 single-nucleotide polymorphisms (SNPs) and 123 short tandem repeats (STRs) with MPS technology on the BGISEQ-500™ platform. First, the library preparation, sequencing process, and data analysis were performed, focusing on the average depth of coverage and heterozygote balance. We calculated the allelic frequencies and forensic parameters of STR and SNP loci in 73 unrelated Chinese Han individuals. In addition, performance was evaluated with accuracy, uniformity, sensitivity, PCR inhibitor, repeatability and reproducibility, mixtures, degraded samples, case-type samples, and pedigree analyses. The results showed that 100% accurate and concordant genotypes can be obtained, and the loci with an abundance in the interquartile range accounted for 92.90% of the total, suggesting reliable uniformity in this panel. We obtained a locus detection rate that was higher than 98.78% from 78 pg of input DNA, and the optimal amount was 1.25-10 ng. The maximum concentrations of hematin and humic acid were 200 and 100 μM, respectively (the ratios of detected loci were 96.52% and 92.41%), in this panel. As a mixture, compared with those of SNPs, minor-contributor alleles of STRs could be detected at higher levels. For the degraded sample, the ratio of detected loci was 98.41%, and most profiles from case-type samples were not significantly different in abundance in our studies. As a whole, this panel showed high-performance, reliable, robust, repeatable, and reproducible results, which are sufficient for paternity testing, individual identification, and use for potentially degraded samples in forensic science.

Full mtDNA genome sequencing of Brazilian admixed populations: A forensic-focused evaluation of a MPS application as an alternative to Sanger sequencing methods

The use of Massive Parallel Sequencing (MPS) techniques have been proposed by the forensic community as an alternative to Sanger sequencing methods in routine forensic casework analysis regarding mitochondrial DNA (mtDNA). Interesting features of MPS include high throughput, ability to simultaneously genotype a significant number of samples by barcoding techniques, processing automation, reduced time and costs, among others. Advantages include the capability of generating full mtDNA genome sequences versus usual techniques, usually limited to hypervariable or control regions exclusively. In this work, 96 reference single-source samples from three different Brazilian cities were subjected to full mtDNA genome sequencing by MPS techniques using an early-access version of Precision ID mtDNA Whole Genome Panel on an Ion Torrent PGM platform (Thermo Fisher Scientific, Waltham, MA, USA). Complete, high-quality sequences were obtained and sequencing performance was evaluated via four different metrics. As a subset of evaluated samples have been previously submitted for Sanger sequencing of the control region, a comparative analysis of both methods' results was conducted in order to compare technique adequacy within a forensic context. Even though this study is one of the first to report full mtDNA genome sequences for Brazilian admixed populations, the observed haplotypes exhibit a predominance of Native American and African maternal lineages in the studied sample set, reproducing results described in the literature for control regions only. Interpopulation analysis among Brazilian and 26 worldwide populations was also carried out. The results indicate that MPS-generated full mtDNA genome sequences may have great utility in forensic real casework applications, with a pronounced gain of genetic information and discrimination power provided by coding region evaluation and the enhanced capacity of heteroplasmies determination. Database construction and other relevant factors concerning implementation of such techniques in Brazilian forensic laboratories are also discussed.