Increasing the discrimination power of ancestry- and identity-informative SNP loci within the ForenSeq™ DNA Signature Prep Kit

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.

Developmental validation of the ForenSeq® Kintelligence kit, MiSeq FGx® sequencing system and ForenSeq Universal Analysis Software

Forensic Investigative Genetic Genealogy, a recent sub discipline of forensic genomics, leverages the high throughput and sensitivity of detection of next generation sequencing and established genetic and genealogical approaches to support the identification of human remains from missing persons investigations and investigative lead generation in violent crimes. To facilitate forensic DNA evidence analysis, the ForenSeq® Kintelligence multiplex, consisting of 10,230 SNPs, was developed. Design of the ForenSeq Kintelligence Kit, the MiSeq FGx® Sequencing System and the ForenSeq Universal Analysis Software is described. Developmental validation in accordance with SWGDAM guidelines and forensic quality assurance standards, using single source samples, is reported for the end-to-end workflow from library preparation to data interpretation. Performance metrics support the conclusion that more genetic information can be obtained from challenging samples compared to other commercially available forensic targeted DNA assays developed for capillary electrophoresis (CE) or other current next generation sequencing (NGS) kits due to the higher number of markers, the overall shorter amplicon sizes (97.8% <150 bp), and kit design. Data indicate that the multiplex is robust and fit for purpose for a wide range of quantity and quality samples. The ForenSeq Kintelligence Kit and the Universal Analysis Software allow transfer of the genetic component of forensic investigative genetic genealogy to the operational forensic laboratory.

Performance and characterization of 94 identity-informative SNPs in Northern Han Chinese using ForenSeq ™ DNA signature prep kit

Target and flanking region (FR) variation at 94 identity-informative SNPs (iSNPs) are investigated in 635 Northern Han Chinese using the ForenSeq DNA Signature Prep Kit on the MiSeq FGx Forensic Genomics System. The dataset presents the following performance characteristics (average values): ≥60% bases with a quality score of 20 or higher (%≥ Q20); >700 × of depth of coverage (DoC) from both Sample Details Reports and Flanking Region Reports; >80% of effective reads; ≥60% of allele coverage ratio (ACR); and ≥70% of inter-locus balance, while some stable low-performance characteristics are also observed: low DoC at rs1736442, rs1031825, rs7041158, rs338882, rs2920816, rs1493232, rs719366, and rs2342747; high noise at rs891700; and imbalanced ACR at rs6955448 and rs338882. The average amplicon length is 69 bp, suitable for detecting degraded samples. Bioinformatic concordance achieves 99.99% between the ForenSeq Universal Analysis Software (UAS) and the Integrative Genomic Viewer (IGV) inspection. Discordance results from flanking region deletions of rs10776839, rs8078417, rs2831700, and rs1454361. Due to FR variants within amplicons detected by massively parallel sequencing (MPS), the increases in the number of unique alleles, effective alleles (Ae), and observed heterozygosity (Hobs) are 46.81%, 4.51%, and 3.29%, respectively. Twelve FR variants are first reported to dbSNP, such as rs1252699848, rs1665500714, rs1771121532, rs2097285015, rs1851671415, rs2045669877, rs2046758811, rs2044248635, rs1251308240, rs1968822112, rs1981638299, and rs1341756746. All 94 iSNPs from target and amplicon data are in Hardy-Weinberg equilibrium (HWE) and independent within autosomes. As expected, forensic parameters from the amplicon data increase significantly on the combined power of discrimination (CPD = 1 - 3.9876 × 10-38) and the combined power of exclusion (CPE = 1 - 6.6690 × 10-8). Additionally, the power of the system effectiveness (CPD = 1 - 6.7054 × 10-72 and CPE = 1 - 4.4719 × 10-20) with sequence-based 27 autosomal STRs and 94 iSNP amplicons in combination is substantially improved compared to one type of marker alone. In conclusion, we have established a traditional length-based and current sequence-based reference database with 58 STRs and 94 iSNPs in the Northern Han Chinese population. We hope these data can serve as a solid reference and foundation for forensic practice.

Population genetic analyses of Eastern Chinese Han nationality using ForenSeq™ DNA Signature Prep Kit

Currently, the most commonly used method for human identification and kinship analysis in forensic genetics is the detection of length polymorphism in short tandem repeats (STRs) using polymerase chain reaction (PCR) and capillary electrophoresis (CE). However, numerous studies have shown that considerable sequence variations exist in the repeat and flanking regions of the STR loci, which cannot be identified by CE detection. Comparatively, massively parallel sequencing (MPS) technology can capture these sequence differences, thereby enhancing the identification capability of certain STRs. In this study, we used the ForenSeq™ DNA Signature Prep Kit to sequence 58 STRs and 94 individual identification SNPs (iiSNPs) in a sample of 220 unrelated individuals from the Eastern Chinese Han population. Our aim is to obtain MPS-based STR and SNP data, providing further evidence for the study of population genetics and forensic applications. The results showed that the MPS method, utilizing sequence information, identified a total of 486 alleles on autosomal STRs (A-STRs), 97 alleles on X-chromosome STRs (X-STRs), and 218 alleles on Y-chromosome STRs (Y-STRs). Compared with length polymorphism, we observed an increase of 260 alleles (157, 31, and 72 alleles on A-STRs, X-STRs, and Y-STRs, respectively) across 36 STRs. The most substantial increments were observed in DYF387S1 and DYS389II, with increases of 287.5% and 250%, respectively. The most increment in the number of alleles was found at DYF387S1 and DYS389II (287.5% and 250%, respectively). The length-based (LB) and sequence-based (SB) combined random match probability (RMP) of 27 A-STRs were 6.05E-31 and 1.53E-34, respectively. Furthermore, other forensic parameters such as total discrimination power (TDP), cumulative probability of exclusion of trios (CPEtrio), and duos (CPEduo) were significantly improved when using the SB data, and informative data were obtained for the 94 iiSNPs. Collectively, these findings highlight the advantages of MPS technology in forensic genetics, and the Eastern Chinese Han genetic data generated in this study could be used as a valuable reference for future research in this field.

Development of a coding SNP panel for tracking the origin of whole-exome sequencing samples

Whole-exome sequencing (WES) is widely used to diagnose complex genetic diseases and rare conditions. The implementation of a robust and effective quality control system for sample identification and tracking throughout the WES process is essential. We established a multiplex panel that included 22 coding single-nucleotide polymorphism (cSNP) loci. The personal identification and paternity identification abilities of the panel were evaluated, and a preliminary validation of the practical feasibility of the panel was conducted in a clinical WES case. These results indicate that the cSNP panel could be a useful tool for sample tracking in WES.

Recombulator-X: A fast and user-friendly tool for estimating X chromosome recombination rates in forensic genetics

Genetic markers (especially short tandem repeats or STRs) located on the X chromosome are a valuable resource to solve complex kinship cases in forensic genetics in addition or alternatively to autosomal STRs. Groups of tightly linked markers are combined into haplotypes, thus increasing the discriminating power of tests. However, this approach requires precise knowledge of the recombination rates between adjacent markers. The International Society of Forensic Genetics recommends that recombination rate estimation on the X chromosome is performed from pedigree genetic data while taking into account the confounding effect of mutations. However, implementations that satisfy these requirements have several drawbacks: they were never publicly released, they are very slow and/or need cluster-level hardware and strong computational expertise to use. In order to address these key concerns we developed Recombulator-X, a new open-source Python tool. The most challenging issue, namely the running time, was addressed with dynamic programming techniques to greatly reduce the computational complexity of the algorithm. Compared to the previous methods, Recombulator-X reduces the estimation times from weeks or months to less than one hour for typical datasets. Moreover, the estimation process, including preprocessing, has been streamlined and packaged into a simple command-line tool that can be run on a normal PC. Where previous approaches were limited to small panels of STR markers (up to 15), our tool can handle greater numbers (up to 100) of mixed STR and non-STR markers. In conclusion, Recombulator-X makes the estimation process much simpler, faster and accessible to researchers without a computational background, hopefully spurring increased adoption of best practices.

US Population Data for 94 Identity-Informative SNP Loci

The US National Institute of Standards and Technology (NIST) analyzed a set of 1036 samples representing four major US population groups (African American, Asian American, Caucasian, and Hispanic) with 94 single nucleotide polymorphisms (SNPs) used for individual identification (iiSNPs). The compact size of iiSNP amplicons compared to short tandem repeat (STR) markers increases the likelihood of successful amplification with degraded DNA samples. Allele frequencies and relevant forensic statistics were calculated for each population group as well as the aggregate population sample. Examination of sequence data in the regions flanking the targeted SNPs identified additional variants, which can be combined with the target SNPs to form microhaplotypes (multiple phased SNPs within a short-read sequence). Comparison of iiSNP performance with and without flanking SNP variation identified four amplicons containing microhaplotypes with observed heterozygosity increases of greater than 15% over the targeted SNP alone. For this set of 1036 samples, comparison of average match probabilities from iiSNPs with the 20 CODIS core STR markers yielded an estimate of 1.7 × 10^-38 for iiSNPs (assuming independence between all 94 SNPs), which was four orders of magnitude lower (more discriminating) than STRs where internal sequence variation was considered, and 10 orders of magnitude lower than STRs using established capillary electrophoresis length-based genotypes.

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.

Forensic identity SNPs: Characterisation of flanking region variation using massively parallel sequencing

Single nucleotide polymorphisms (SNPs) can be analysed for identity or kinship applications in forensic genetics to either provide an adjunct to traditional STR typing or as a stand-alone approach. The advent of massively parallel sequencing technology (MPS) has provided a useful opportunity to more easily deploy SNP typing in a forensic context, given the ability to simultaneously amplify a large number of markers. Furthermore, MPS also provides valuable sequence data for the targeted regions, which enables the detection of any additional variation seen in the flanking regions of amplicons. In this study we genotyped 977 samples across five UK-relevant population groups (White British, East Asian, South Asian, North-East African and West African) for 94 identity-informative SNP markers using the ForenSeq DNA Signature Prep Kit. Examination of flanking region variation allowed for the identification of 158 additional alleles across all populations studied. Here we present allele frequencies for all 94 identity-informative SNPs, both including and excluding the flanking region sequence of these markers. We also present information on the configuration of these SNPs in the ForenSeq DNA Signature Prep Kit, including performance metrics for the markers and investigation of bioinformatic and chemistry-based discordances. Overall, the inclusion of flanking region variation in the analysing workflow for these markers reduced the average combined match probability 2175 times across all populations, with a maximum reduction of 675,000-fold in the West African population. The gain due to flanking region-based discrimination increased the heterozygosity of some loci above that of some of the least useful forensic STR loci; thus demonstrating the benefit of enhanced analysis of currently targeted SNP markers for forensic applications.

Application of a custom haplotype caller to analyze sequence-based data of 56 microhaplotypes

Microhaplotypes (microhaps) are recently introduced markers that aim to complement the limitations of conventional forensic markers such as short tandem repeats (STRs) and single nucleotide polymorphisms (SNPs). With the potential of microhaps in forensics becoming clearer through massively parallel sequencing (MPS), MPS-based studies on microhaps are being actively reported. However, simpler workflow schemes for the generation and analysis of MPS data are still required to facilitate the practical application of MPS in forensics. In this study, we developed an in-house MPS panel that simultaneously amplifies 56 microhaps and a custom haplotype caller, Visual Microhap. The developed tool works on a web browser and provides four analysis options to extract SNP-based haplotypes from sequence-based data obtained by STRait Razor 3.0. To demonstrate the utility of the MPS panel and data analysis workflow scheme, we also analyzed 56 microhaps of 286 samples from four populations (African-American, Caucasian, Hispanic, and Korean). The average effective number of alleles (A_e) for the four groups was 3.45, ranging from 1.74 to 6.98. Forensic statistical parameters showed that this microhap panel is more powerful than conventional autosomal STRs for human identification. Meanwhile, the 56-plex panel mostly comprised microhaps with high Ae; however, the four populations were grossly distinguishable from each other by cluster analysis. Consequently, the developed in-house MPS panel for 56 microhaps and the adopted workflow using open-source tools can increase the utility of microhap MPS in forensic research and practice.

Peptide ancestry informative markers in uterine neoplasms from women of European, African, and Asian ancestry

Characterization of ancestry-linked peptide variants in disease-relevant patient tissues represents a foundational step to connect patient ancestry with disease pathogenesis. Nonsynonymous single-nucleotide polymorphisms encoding missense substitutions within tryptic peptides exhibiting high allele frequencies in European, African, and East Asian populations, termed peptide ancestry informative markers (pAIMs), were prioritized from 1000 genomes. In silico analysis identified that as few as 20 pAIMs can determine ancestry proportions similarly to >260K SNPs (R² = 0.99). Multiplexed proteomic analysis of >100 human endometrial cancer cell lines and uterine leiomyoma tissues combined resulted in the quantitation of 62 pAIMs that correlate with patient race and genotype-confirmed ancestry. Candidates include a D451E substitution in GC vitamin D-binding protein previously associated with altered vitamin D levels in African and European populations. pAIMs will support generalized proteoancestry assessment as well as efforts investigating the impact of ancestry on the human proteome and how this relates to the pathogenesis of uterine neoplasms.

Evaluation of a custom QIAseq targeted DNA panel with 164 ancestry informative markers sequenced with the Illumina MiSeq

Introduction of new methods requires meticulous evaluation before they can be applied to forensic genetic case work. Here, a custom QIAseq Targeted DNA panel with 164 ancestry informative markers was assessed using the MiSeq sequencing platform. Concordance, sensitivity, and the capability for analysis of mixtures were tested. The assay gave reproducible and nearly concordant results with an input of 10 and 2 ng DNA. Lower DNA input led to an increase in both locus and allele drop-outs, and a higher variation in heterozygote balance. Locus or allele drop-outs in the samples with less than 2 ng DNA input were not necessarily associated with the overall performance of a locus. Thus, the QIAseq assay will be difficult to implement in a forensic genetic setting where the sample material is often scarce and of poor quality. With equal or near equal mixture ratios, the mixture DNA profiles were easily identified by an increased number of imbalanced heterozygotes. For more skewed mixture ratios, the mixture DNA profiles were identified by an increased noise level. Lastly, individuals from Great Britain and the Middle East were investigated. The Middle Eastern individuals showed a greater affinity with South European populations compared to North European populations.

Evaluation of ArmedXpert software tools, MixtureAce and Mixture Interpretation, to analyze MPS-STR data

Massively parallel sequencing (MPS) technologies have revolutionized studies of genomic variations and transformed DNA analysis in multiple fields. Assays based on MPS must be capable of discriminating variations introduced by the method, i.e. artifacts from true polymorphisms. In PCR-MPS methods targeting microsatellite markers, artifacts can arise from PCR mis-incorporation, PCR strand slippage (stutter), and sequencing error. Reliable detection of artifacts in mixed DNA samples is a significant challenge that must be addressed in forensic DNA analysis. The ArmedXpert (NicheVision) software tools, MixtureAce™ and Mixture Interpretation, can analyze MPS data by categorizing sequence reads in alleles, stutter, and non-stutter artifacts and analyzing autosomal STR loci of mixed samples. In this study, we evaluated the ArmedXpert tools for the analysis of STR profiles of single-sourced and mixed samples generated by the ForenSeq™ DNA Signature Prep kit (Verogen). Data from eight experimental runs (240 samples) were analyzed: one benchmark run, two runs testing sensitivity with down to 50 pg DNA input, one run testing artificially degraded samples and DNA derived from bones, blood cards and teeth, as well as four runs with mixed DNA samples of varying ratios, sex, and different number of contributors (two to six). The MixtureAce stutter thresholds were initially set following the recommendations from Verogen, plus a non-stutter artifact threshold was set at 5% of allele read counts. A benchmark run, of 30 samples, plus two controls, containing 2310 total alleles, revealed over 5000 artifacts, above an analytical threshold of 10. A total of 4869 artifacts were correctly classified, while 435 were mis-classified as alleles due to exceedance of initial threshold settings. False positives must be resolved by an analyst, which can be time consuming. Stutter thresholds were adjusted based on the benchmark data and the samples were re-tested, resulting in only 57 false positive allele calls. The revised settings were then used in the analysis of the remaining seven experimental runs. Results show that MixtureAce can accurately classify artifacts and alleles when laboratory-specific threshold settings are used. The Mixture Interpretation tool was applied on two- and three-person mixtures. This tool utilized the analyzed data from MixtureAce to calculate, based on the number of alleles at a locus and their read counts, possible deconvolution outcomes with their respective ratios.

High-throughput DNA sequencing of environmentally insulted latent fingerprints after visualization with nanoscale columnar-thin-film technique

The goals of this study were to (a) ascertain human identity capabilities of DNA obtained from latent fingerprints that have been first environmentally insulted and then developed by the deposition of a columnar thin film (CTF), and (b) to determine if the CTF process and material are detrimental to single nucleotide polymorphism (SNP) analysis. Fingerprints were deposited on five different types of substrates and aged for one day, 7 days or 30 days while being environmentally insulted under one of the four conditions: 16.6 °C and 60% relative humidity (RH) (Condition A), 24.5 °C and 60% RH (Condition B), 35 °C and 67% RH (Condition C) and a cold condition (Condition D). Then CTF technique was then on 59% of these fingerprints. DNA samples from 805 fingerprints were extracted, quantified, subjected to manual library preparation using the Precision ID Identity Panel, and underwent high-throughput sequencing. The Ion S5™ platform was employed to sequence 124 SNP amplicons. SNPs were successfully sequenced from 802/805 samples. Total read depth was consistent across environmental conditions, and majority of samples had 100% profile completeness and 100% concordance. Anecdotally, libraries that were amplified with a higher cycle number had more 'Major Allele Frequency' flags compared to samples amplified with 23 cycle numbers, possibly due to stochastic effects. Neither the substrates nor the CTF process and materials inhibit downstream DNA analysis. DNA of low quality and quantity from the chosen samples can be sequenced using the Precision ID Identity Panel on the Ion S5™ platform which performed well, however, a different approach may be needed if spurious alleles are suspected.

Reverse complement-PCR, an innovative and effective method for multiplexing forensically relevant single nucleotide polymorphism marker systems

DNA analyses from challenging samples such as touch evidence, hairs and skeletal remains push the limits of the current forensic DNA typing technologies. Reverse complement PCR (RC-PCR) is a novel, single-step PCR target enrichment method adapted to amplify degraded DNA. The sample preparation process involves a limited number of steps, decreasing the labor required for library preparation and reducing the possibility of contamination due to less sample manipulation. These features of the RC-PCR make the technology a unique application to successfully target single nucleotide polymorphisms (SNPs) in fragmented and low copy number DNA and yield results from samples in which no or limited data are obtained with standard DNA typing methods. The developed RC-PCR short amplicon 85 SNP-plex panel is a substantial improvement over the previously reported 27-plex RC-PCR multiplex that will provide higher discrimination power for challenging DNA sample analyses.

A targeted ancestry informative InDels panel on capillary electrophoresis for ancestry inference in Asian populations

CE is the primary methodology used in forensic DNA typing. Alleles of commonly used types of genetic markers could be separated and detected via CE based on dye color and migration time. Insertion/deletion (InDel) is an ideal genetic marker for forensic DNA analysis due to their abundance in the human genome, low mutation rate, availability of their allele types via CE, and elimination of stutter peaks. Moreover, InDels could be used as ancestry informative markers since allele frequencies of InDels is different among geographically separated populations. Several ancestry informative insertion/deletion panels have been established based on CE platform to achieve the intercontinental populations distinction. However, improvements to differentiate intracontinental populations is few. In this study, 21 InDels with fixation index (F_ST ) > 0.15 were selected and assembled into one ancestry informative insertion/deletion panel. Using well-designed primers, those 21 InDels could be amplified successfully and genotyped on the CE platform accurately and completely. The panel showed a large F_ST distance distinction among the ten Asian populations. Using clustering analysis, ten Asian populations were classified into three subgroups: East Asian, Southeast Asian, and South Asian subgroups. To evaluate the panel's capability in ancestry inference, a validation experiment was undertaken with 319 individuals from four geographically separated populations in China. Four Chinese populations were classified into different ancestry subgroups and 81.8% test individuals' ancestry could be inferred correctly. Our result showed that development of high ancestry informative InDels panel based on CE platform is a potential for individual ancestry inference among intracontinental populations.

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.

Autosomal STR and SNP characterization of populations from the Northeastern Peruvian Andes with the ForenSeq™ DNA Signature Prep Kit

Autosomal DNA data from Peru for human identity testing purposes are scarce in the scientific literature, which hinders obtaining an appropriate portrait of the genetic variation of the resident populations. In this study we genetically characterize five populations from the Northeastern Peruvian Andes (Chachapoyas, Awajún, Wampís, Huancas and Cajamarca). Autosomal short tandem repeat (aSTR) and identity informative single nucleotide polymorphism (iiSNP) data from a total of 233 unrelated individuals are provided, and forensic genetic parameters are calculated for each population and for the combined set Northeastern Peruvian Andes. After correction for multiple testing in the whole dataset of the Northeastern Peruvian Andes, the only departure from Hardy-Weinberg equilibrium was observed in locus rs2111980. Twenty one out of 27 aSTR loci exhibited an increased number of alleles due to sequence variation in the repeat motif and flanking regions. For iiSNPs 33% of the loci displayed flanking region variation. The combined random match probability (RMP), assuming independence of all loci (aSTRs and iiSNPs), in the Chachapoyas, the population with the largest samples size (N = 172), was 8.14 × 10^-62 for length-based data while for sequence-based was 4.15 × 10^-67. In the merged dataset (Northeastern Peruvian Andes; N = 233), the combined RMP when including all markers were 2.96 × 10^-61 (length-based) and 3.21 × 10^-66 (sequence-based). These new data help to fill up some of the gaps in the genetic canvas of South America and provide essential length- and sequence-based background information for other forensic genetic studies in Peru.

Evaluation of OpenArray™ as a Genotyping Method for Forensic DNA Phenotyping and Human Identification

A custom plate of OpenArray™ technology was evaluated to test 60 single-nucleotide polymorphisms (SNPs) validated for the prediction of eye color, hair color, and skin pigmentation, and for personal identification. The SNPs were selected from already validated subsets (Hirisplex-s, Precision ID Identity SNP Panel, and ForenSeq DNA Signature Prep Kit). The concordance rate and call rate for every SNP were calculated by analyzing 314 sequenced DNA samples. The sensitivity of the assay was assessed by preparing a dilution series of 10.0, 5.0, 1.0, and 0.5 ng. The OpenArray™ platform obtained an average call rate of 96.9% and a concordance rate near 99.8%. Sensitivity testing performed on serial dilutions demonstrated that a sample with 0.5 ng of total input DNA can be correctly typed. The profiles of the 19 SNPs selected for human identification reached a random match probability (RMP) of, on average, 10⁻⁸. An analysis of 21 examples of biological evidence from 8 individuals, that generated single short tandem repeat profiles during the routine workflow, demonstrated the applicability of this technology in real cases. Seventeen samples were correctly typed, revealing a call rate higher than 90%. Accordingly, the phenotype prediction revealed the same accuracy described in the corresponding validation data. Despite the reduced discrimination power of this system compared to STR based kits, the OpenArray™ System can be used to exclude suspects and prioritize samples for downstream analyses, providing well-established information about the prediction of eye color, hair color, and skin pigmentation. More studies will be needed for further validation of this technology and to consider the opportunity to implement this custom array with more SNPs to obtain a lower RMP and to include markers for studies of ancestry and lineage.

STRait Razor Online: An enhanced user interface to facilitate interpretation of MPS data

Since 2013, STRait Razor has enabled analysis of massively parallel sequencing (MPS) data from various marker systems such as short tandem repeats, single nucleotide polymorphisms, insertion/deletions, and mitochondrial DNA. In this paper, STRait Razor Online (SRO), available at https://www.unthsc.edu/straitrazor, is introduced as an interactive, Shiny-based user interface for primary analysis of MPS data and secondary analysis of STRait Razor haplotype pileups. This software can be accessed from any common browser via desktop, tablet, or smartphone device. SRO is available also as a standalone application and open-source R script available at https://github.com/ExpectationsManaged/STRaitRazorOnline. The local application is capable of batch processing of both fastq files and primary analysis output. Processed batches generate individual report folders and summary reports at the locus- and haplotype-level in a matter of minutes. For example, the processing of data from ∼700 samples generated with the ForenSeq Signature Preparation Kit from allsequences.txt to a final table can be performed in ∼40 min whereas the Excel-based workbooks can take 35-60 h to compile a subset of the tables generated by SRO. To facilitate analysis of single-source, reference samples, a preliminary triaging system was implemented that calls potential alleles and flags loci suspected of severe heterozygote imbalance. When compared to published, manually curated data sets, 98.72 % of software-assigned allele calls without manual interpretation were consistent with curated data sets, 0.99 % loci were presented to the user for interpretation due to heterozygote imbalance, and the remaining 0.29 % of loci were inconsistent due to the analytical thresholds used across the studies.

Massively parallel sequencing and STR analysis from partial bloody fingerprints enhanced with columnar thin films

Fingerprint enhancement often includes either physical or chemical approaches, such as fingerprint powder or cyanoacrylate fuming, to improve the quality of a fingerprint for visualization and analysis. However, these methods become more complex when fingerprints are partial bloody, and these procedures may interfere with downstream DNA analysis. Columnar thin film (CTF) deposition is a type of nanotechnology that utilizes an evaporant material to enhance a fingerprint under low-pressure conditions. Short tandem repeat (STR) analysis is the traditional method employed in crime laboratories. When DNA is of poor quality and quantity, like that often obtained from fingerprints, little to no genetic information may be obtained. Single nucleotide polymorphisms (SNPs) may be used to glean additional information when STR analysis fails. In this pilot study, 100 partial bloody fingerprints were collected from two donors and deposited on five different crime scene substrates, in which half were enhanced with CTFs and were graded for quality by an IAI-certified latent fingerprint examiner. CTF-developed fingerprints, on average, had higher grades compared to non-developed partial bloody fingerprints. STR analysis using Fusion 6C was performed to assess inhibition from the evaporant materials, in which no inhibition was observed. Sequencing of SNPs using the Precision ID Identity Panel was also employed, in which genetic information that could not be obtained from STRs was acquired with SNPs. Various sample types (i.e. pristine, low quality, and contaminated) utilized in this project demonstrated the acceptable performance of the Precision ID Identity Panel.

Identification of sequence polymorphisms at 58 STRs and 94 iiSNPs in a Tibetan population using massively parallel sequencing

Massively parallel sequencing (MPS) has rapidly become a promising method for forensic DNA typing, due to its ability to detect a large number of markers and samples simultaneously in a single reaction, and sequence information can be obtained directly. In the present study, two kinds of forensic genetic markers, short tandem repeat (STR) and identity-informative single nucleotide polymorphism (iiSNP) were analyzed simultaneously using ForenSeq DNA Signature Prep Kit, a commercially available kit on MPS platform. A total of 152 DNA markers, including 27 autosomal STR (A-STR) loci, 24 Y chromosomal STR (Y-STR) loci, 7 X chromosomal STR (X-STR) loci and 94 iiSNP loci were genotyped for 107 Tibetan individuals (53 males and 54 females). Compared with length-based STR typing methods, 112 more A-STR alleles, 41 more Y-STR alleles, and 24 more X-STR alleles were observed at 17 A-STRs, 9 Y-STRs, and 5 X-STRs using sequence-based approaches. Thirty-nine novel sequence variations were observed at 20 STR loci. When the flanking regions were also analyzed in addition to target SNPs at the 94 iiSNPs, 38 more alleles were identified. Our study provided an adequate genotype and frequencies data of the two types of genetic markers for forensic practice. Moreover, we also proved that this panel is highly polymorphic and informative in Tibetan population, and should be efficient in forensic kinship testing and personal identification cases.

Validation of novel forensic DNA markers using multiplex microhaplotype sequencing

Microhaplotypes (MH) are comprised of multiple single nucleotide polymorphisms (SNPs) that are located within 300 bases of genomic sequence. Improved tools are needed to facilitate broader application of microhaplotypes in a diverse range of populations and forensic settings. We designed an assay for multiplex sequencing of 90 microhaplotypes (mMHseq) that include 46 MH loci with high Effective Number of Alleles (A_e) from previous studies [1], and 44 high A_e MH loci containing between four to fourteen SNPs that were identified from the 1000 Genomes (1KG) Project. The unique design of mMHseq integrates a novel method for multiplex amplification from small DNA amounts, and multiplex sequencing of 48 samples in a single MiSeq run to detect all relevant MH variation. Assay performance was evaluated in a cohort of 156 individuals from seven different world populations from Africa, Asia, and Europe. Three of those populations from East Africa (Chagga, Sandawe, and Zaramo) and one from Eastern Europe (Adygei) had sufficient individuals sequenced by the assay to be included in statistical analyses with the 26 1KG populations. For those 30 populations the mean global average A_e was 5.08 (range: 2.7-11.54) and mean informativeness for biogeographic variation (I_n) was 0.30 (range: 0.08-0.70). Eighty-five novel SNPs were detected in 58 of the 90 microhaplotypes. Open-source, web-based software was developed to visualize haplotype phase data for each microhaplotype and individual. Our approach for multiplex microhaplotype sequencing can be customized and expanded as novel loci are being discovered.

Sequenced-based French population data from 169 unrelated individuals with Verogen's ForenSeq DNA signature prep kit

Massively Parallel Sequencing (MPS) applied to forensic genetics allows the simultaneous analysis of hundreds of genetic markers and the access to full amplicon sequences which help to increase available allele diversity. Meanwhile, sequence variation within the repeat regions represents the majority of the allele diversity, flanking regions adjacent to the repeat core provide an additional degree of variation. The forensic genetics community needs access to population data, from relevant parts of the world that contain this new sequence diversity in order to perform statistical calculations. In this study, we report sequence-based Short Tandem Repeat (STR) and identity Single Nucleotide Polymorphism (iSNPs) allele data for 169 French individuals across 58 STRs and 92 SNPs included in the Verogen ForenSeq DNA Signature Prep kit. 42 STRs out of 58 showed an increased number of alleles due to sequence variation in the repeat motif and/or the flanking regions. D9S1122 showed the largest overall gain with an increase of observed heterozygosities of almost 25 %. The combined match probability combining 27 autosomal STRs and 91 identity SNPs was 1.11E-69. Sequence-based allele frequencies included in this publication will help forensic laboratories to increase the power of discrimination for identification, kinship analysis and mixture interpretation.

An innovative panel containing a set of insertion/deletion loci for individual identification and its forensic efficiency evaluations in Chinese Hui ethnic minority

BACKGROUND

Individual identification is one of the most important tasks in the field of forensic genetics. Insertion/Deletion (InDel) polymorphism marker has been a promising marker for individual identification. However, a part of InDel loci in commonly used commercial kit show low polymorphisms in Chinese populations.

METHODS

We evaluated a panel of 35 InDel loci constructed previously for individual identifications in Hui group. Subsequently, population data of three Chinese populations from 1,000 Genomes Project database were used to evaluate individual identification performance of these 35 InDels. Forensic parameters, such as heterozygosity, power of exclusion, match probability and power of discrimination, were calculated to evaluate the forensic efficiency of these loci in Hui group. The heatmap of insertion allelic frequencies, Nei's genetic distances, pairwise fixation index values, principal component analyses and admixture analyses were used to analyze the genetic differentiations and structure between Hui group and other populations.

RESULTS

In studied Hui group, besides rs3054057, polymorphism information content values of the remaining loci were greater than 0.3. Values of expected heterozygosity of these loci were close to 0.5. The combined power of discrimination and power of exclusion values were 0.99999999999999659609 and 0.998682, respectively. Analyses of population genetics revealed that Chinese Hui group had closer genetic relationships with East Asian populations than other intercontinental populations.

CONCLUSION

The forensic statistical analyses revealed these loci showed relatively high genetic polymorphisms in Chinese Hui group, and could be served as a useful tool for individual identifications in Hui group. Population genetic evaluations indicated that Chinese Hui group had close genetic relationships with East Asian populations.

Cadaveric blood cards: Assessing DNA quality and quantity and the utility of STRs for the individual estimation of trihybrid ancestry and admixture proportions

As part a body donation program, blood samples were collected and stored on untreated (non-FTA) blood cards. The blood cards were evaluated in terms of DNA preservation and STR typing success with resulting profiles assessed with special consideration given to profile matching for positive identification and biogeographic ancestry estimation. While STR profiles were successfully generated for all samples, results indicate that the time interval between date of death and sample collection have an impact on DNA quantity and quality. There is a statistically significant decrease in relative fluorescent unit (RFU) values with increasing time interval between date of death and sample collection, indicating degradation in the blood card samples related to the post-mortem interval prior to sample collection. The STR profiles were used to estimate ancestry and admixture using the program STRUCTURE, demonstrating utility of these markers beyond individual identification purposes, with caveats for application based on population history.

Massively parallel sequencing of autosomal STRs and identity-informative SNPs highlights consanguinity in Saudi Arabia

While many studies have been undertaken of Middle Eastern populations using autosomal STR profiling by capillary electrophoresis, little has so far been published from this region on the forensic use of massively parallel sequencing (MPS). Here, we carried out MPS of 27 autosomal STRs and 91 identity-informative SNPs (iiSNPs) with the Verogen ForenSeq™ DNA Signature Prep Kit on a representative sample of 89 Saudi Arabian males, and analysed the resulting sequence data using Verogen's ForenSeq Universal Analysis Software (UAS) v1.3 and STRait Razor v3.0. This revealed sequence variation in the composition of complex STR arrays, and SNPs in their flanking regions, which raised the number of STR alleles from 238 distinct length variants to 357 sequence sub-variants. Similarly, between one and three additional polymorphic sites were observed within the amplicons of 37 of the 91 iiSNPs, forming up to six microhaplotypes per locus. These further enhance discrimination compared to the biallelic target SNP data presented by the primary UAS interface. In total, we observed twenty-two STR alleles previously unrecognised in the STRait Razor v3.0 default allele list, along with nine SNPs flanking target iiSNPs that were not highlighted by the UAS. Sequencing reduced the STR-based random match probability (RMP) from 2.62E-30 to 3.49E-34, and analysis of the iiSNP microhaplotypes reduced RMP from 9.97E-37 to 6.83E-40. The lack of significant linkage disequilibrium between STRs and target iiSNPs allowed the two marker types to be combined using the product rule, yielding a RMP of 2.39E-73. Evidence of consanguinity was apparent from both marker types. While TPOX was the only locus displaying a significant deviation from Hardy-Weinberg equilibrium, 23 out of 27 STRs and 63 out of 91 iiSNPs showed fewer than expected heterozygotes, demonstrating an overall homozygote excess probably reflecting the high frequency of first-cousin marriages in Saudi Arabia. We placed our data in a global context by considering the same markers in the Human Genome Diversity Panel (HGDP), revealing that the Saudi sample was typical of Middle Eastern populations, with a higher level of inbreeding than is seen in most European, African and Central/South Asian populations, correlating with known patterns of endogamy. Given reduced levels of diversity within endogamous groups, the ability to combine the discrimination power of both STRs and SNPs offers significant benefits in the analysis of forensic evidence in Saudi Arabia and the Middle East region more generally.

Utility of the Ion S5™ and MiSeq FGx™ sequencing platforms to characterize challenging human remains

Often in missing persons' and mass disaster cases, the samples remaining for analysis are hard tissues such as bones, teeth, nails, and hair. These remains may have been exposed to harsh environmental conditions, which pose challenges for downstream genotyping. Short tandem repeat analysis (STR) via capillary electrophoresis (CE) is still the gold standard for DNA typing; however, a newer technology known as massively parallel sequencing (MPS) could improve upon our current techniques by typing different and more markers in a single analysis, and consequently improving the power of discrimination. In this study, bone and tooth samples exposed to a variety of DNA insults (cremation, embalming, decomposition, thermal degradation, and fire) were assessed and sequenced using the Precision ID chemistry and a custom AmpliSeq™ STR and iiSNP panel on the Ion S5™ System, and the ForenSeq DNA Signature Prep Kit on the MiSeq FGx™ system, as well as the GlobalFiler™ PCR Amplification Kit on the 3500™ Genetic Analyzer. The results demonstrated that using traditional CE-based genotyping performed as expected, producing a partial or full DNA profile for all samples, and that both sequencing chemistries and platforms were able to recover sufficient STR and SNP information from a majority of the same challenging samples. Run metrics including profile completeness and mean read depth produced good results with each system, considering the degree of damage of some samples. Most sample insults (except decomposed) produced similar numbers of alleles for both MPS systems. Comparable markers produced full concordance between the two platforms.

Identity informative SNP associations in the UK Biobank

Single nucleotide polymorphisms (SNPs) are amenable to genotyping DNA from degraded, inhibited, and/or ancient substrates due to their relatively small amplicon size. Though they have clear advantages over traditional short tandem repeat (STR) typing for specific casework scenarios, the advances in massively parallel sequencing (MPS) have drastically increased the utility of this marker type. The biallelic nature of SNPs makes them individually less informative than STRs due to limited heterozygosity; however, in sufficiently large multiplexes, identity informative SNPs (iiSNPs) may produce combined random match probabilities comparable to STR typing. Multiple MPS library preparation kits now include iiSNPs and similar to STRs, these loci have been rigorously characterized during multiplex development. The relative accessibility of genome-wide association study (GWAS) summary statistics enables re-investigation of forensically relevant targets in high-quality datasets. Here, 4085 GWASs from the UK Biobank European datasets (UKB; 787 ≤ N ≤ 361,194) were mined for iiSNPs typed by the ForenSeq DNA Signature Prep Kit (Verogen). Seven iiSNPs had genome-wide association (p ≤ 5 × 10^-8) with 17 phenotypes in UKB Europeans. Most notably, these relationships involve two outwardly visible characteristics: standing height (rs907100; β = 0.011, p = 1.35 × 10^-10) and hair/balding patterns (rs2399332; β = -0.009, p = 3.83 × 10^-8). The remaining associations involve red blood cell characteristics and measures of lung function. Though these traits are highly polygenic and the individual SNP effects described here have been refuted empirically, we describe the importance and ease of exploring high-quality, freely accessible data to continuously and robustly characterize new and existing forensically relevant loci.