A complete pipeline enables haplotyping and phasing macrohaplotype in long sequencing reads for polyploidy samples and a multi-source DNA mixture.

Macrohaplotype combines multiple types of phased DNA variants, increasing forensic discrimination power. High-quality long-sequencing reads, for example, PacBio HiFi reads, provide data to detect macrohaplotypes in multiploidy and DNA mixtures. However, the bioinformatics tools for detecting macrohaplotypes are lacking. In this study, we developed a bioinformatics software, MacroHapCaller, in which targeted loci (i.e., short TRs [STRs], single nucleotide polymorphisms, and insertion and deletions) are genotyped and combined with novel algorithms to call macrohaplotypes from long reads. MacroHapCaller uses physical phasing (i.e., read-backed phasing) to identify macrohaplotypes, and thus it can detect multi-allelic macrohaplotypes for a given sample. MacroHapCaller was validated with data generated from our designed targeted PacBio HiFi sequencing pipeline, which sequenced ∼8-kb amplicon regions harboring 20 core forensic STR loci in human benchmark samples HG002 and HG003. MacroHapCaller also was validated in whole-genome long-read sequencing data. Robust and accurate genotyping and phased macrohaplotypes were obtained with MacroHapCaller compared with the known ground truth. MacroHapCaller achieved a higher or consistent genotyping accuracy and faster speed than existing tools HipSTR and DeepVar. MacroHapCaller enables efficient macrohaplotype analysis from high-throughput sequencing data and supports applications using discriminating macrohaplotypes.

Determining lineages between individuals with high-density mitochondrial and Y-chromosomal single-nucleotide polymorphisms.

Genetic genealogy has been more frequently used in forensic investigations in identifying criminals. However, the current genetic genealogy applications usually do not consider lineage markers (including both Y and mitochondrial deoxyribonucleic acid (DNA)), which is probably because not all distant relatives share the same lineage markers. In addition, there is no study to show how to use lineage markers and what methods or thresholds should be applied in genetic genealogy. In this study, we developed a method to quickly determine if two single-nucleotide polymorphism (SNP) profiles are from the same paternal or material lineages by using a mismatch frequency of the SNPs in Y-chromosomal or mitochondrial DNA. For both Y and mitochondrial SNPs, profile pairs from the same or different lineages can be decided with high accuracies (i.e., 0.380% or 0.157% error rates with Y and mitochondrial DNA, respectively). With kinship coefficient filtering based on autosomal SNPs, the accuracies of determining maternal and paternal lineage can be further improved (i.e., 0.120% or 0.057% error rates with Y and mitochondrial DNA, respectively, using a threshold of kinship coefficient >0). This study shows that lineage markers can be very powerful tools with high accuracies to determine lineages, which could help solve cases and reduce costs for genetic genealogy investigations.

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.

USAT: a bioinformatic toolkit to facilitate interpretation and comparative visualization of tandem repeat sequences.

BACKGROUND: Tandem repeats (TR), highly variable genomic variants, are widely used in individual identification, disease diagnostics, and evolutionary studies. The recent advances in sequencing technologies and bioinformatic tools facilitate calling TR haplotypes genome widely. Both length-based and sequence-based TR alleles are used in different applications. However, sequence-based TR alleles could provide the highest precision in characterizing TR haplotypes. The need to identify the differences at the single nucleotide level between or among TR haplotypes with an easy-use bioinformatic tool is essential. RESULTS: In this study, we developed a Universal STR Allele Toolkit (USAT) for TR haplotype analysis, which takes TR haplotype output from existing tools to perform allele size conversion, sequence comparison of haplotypes, figure plotting, comparison for allele distribution, and interactive visualization. An exemplary application of USAT for analysis of the CODIS core STR loci for DNA forensics with benchmarking human individuals demonstrated the capabilities of USAT. USAT has user-friendly graphic interfaces and runs fast in major computing operating systems with parallel computing enabled. CONCLUSION: USAT is a user-friendly bioinformatics software for interpretation, visualization, and comparisons of TRs.

A prospective cost-benefit analysis for nylon 4N6FLOQSwabs®: example of the process and potential benefits.

Laboratories and their criminal justice systems are confronted with challenges for implementing new technologies, practices, and policies even when there appears to be demonstrative benefits to operational performance. Impacting decisions are the often higher costs associated with, for example, new technologies, limited current budgets, and making hard decisions on what to sacrifice to take on the seemingly better approach. A prospective cost-benefit analysis (CBA) could help an agency better formulate its strategies and plans and more importantly delineate how a relatively small increase to take on, for example, a new technology can have large impact on the system (e.g., the agency, other agencies, victims and families, and taxpayers). To demonstrate the process and potential value a CBA was performed on the use of an alternate and more expensive swab with reported better DNA yield and being certified human DNA free (i.e., nylon 4N6FLOQSwabs®), versus the traditional less costly swab (i.e., cotton swab). Assumptions are described, potential underestimates and overestimates noted, different values applied (for low and modest to high), and potential benefits (monetary and qualitative) presented. The overall outcome is that the cost of using the more expensive technology pales compared with the potential tangible and intangible benefits. This approach could be a guide for laboratories (and associated criminal justice systems) worldwide to support increased funding, although the costs and benefits may vary locally and for different technologies, practices, and policies. With well-developed CBAs, goals of providing the best services to support the criminal justice system and society can be attained.

Population data for 21 STRs for the Salvadoran population using GlobalFiler Express and GlobalFiler STR Amplification Kits.

With the advent of expanded STR (short tandem repeats) typing kits, it was necessary to determine allele frequencies and other appropriate population data parameters for El Salvador. Samples were collected from the central, east, and west regions of the country and typed for 21 forensically relevant STR loci. The data indicate that all loci are highly polymorphic, the three regions are genetically similar, and the population data are similar to those of US Hispanics. The results of this study support that the allele frequency data described herein can be used for statistical calculations for human identity testing in El Salvador.

Enhanced mixture interpretation with macrohaplotypes based on long-read DNA sequencing.

Deconvoluting mixture samples is one of the most challenging problems confronting DNA forensic laboratories. Efforts have been made to provide solutions regarding mixture interpretation. The probabilistic interpretation of Short Tandem Repeat (STR) profiles has increased the number of complex mixtures that can be analyzed. A portion of complex mixture profiles, particularly for mixtures with a high number of contributors, are still being deemed uninterpretable. Novel forensic markers, such as Single Nucleotide Variants (SNV) and microhaplotypes, also have been proposed to allow for better mixture interpretation. However, these markers have both a lower discrimination power compared with STRs and are not compatible with CODIS or other national DNA databanks worldwide. The short-read sequencing (SRS) technologies can facilitate mixture interpretation by identifying intra-allelic variations within STRs. Unfortunately, the short size of the amplicons containing STR markers and sequence reads limit the alleles that can be attained per STR. The latest long-read sequencing (LRS) technologies can overcome this limitation in some samples in which larger DNA fragments (including both STRs and SNVs) with definitive phasing are available. Based on the LRS technologies, this study developed a novel CODIS compatible forensic marker, called a macrohaplotype, which combines a CODIS STR and flanking variants to offer extremely high number of haplotypes and hence very high discrimination power per marker. The macrohaplotype will substantially improve mixture interpretation capabilities. Based on publicly accessible data, a panel of 20 macrohaplotypes with sizes of ~ 8 k bp and the maximum high discrimination powers were designed. The statistical evaluation demonstrates that these macrohaplotypes substantially outperform CODIS STRs for mixture interpretation, particularly for mixtures with a high number of contributors, as well as other forensic applications. Based on these results, efforts should be undertaken to build a complete workflow, both wet-lab and bioinformatics, to precisely call the variants and generate the macrohaplotypes based on the LRS technologies.

Ancestry inference and admixture component estimations of Chinese Kazak group based on 165 AIM-SNPs via NGS platform.

Predicting the biogeographical ancestries of populations and unknown individuals based on ancestry-informative markers (AIMs) has been widely applied in providing DNA clues to criminal investigations, correcting the factor of population stratification in genome-wide association studies (GWAS), and working as the basis of predicting the externally visible characteristics (EVCs) of individuals. The present study chose Chinese Xinjiang Kazak (XJK) group as research object using a 165 AIM-SNPs panel via next generation sequencing (NGS) technology to reveal its ancestral information and genetic background by referencing the populations' data from 1000 Genomes Phase 3. After the Bonferroni correction, there were no significant deviations at the 165 AIM-SNP loci except two loci with homozygote in the studied XJK group. Ancestry information inference and populations genetic analyses were conducted basing on multiplex statistical methods such as forensic statistical parameter analyses, estimation of the success ratios with cross-validation, population tree, principal component analysis (PCA), and genetic structure analysis. The present results revealed that XJK group had the admixed ancestral components of East Asian and European populations with the ratio of about 62:37.

SNP typing using the HID-Ion AmpliSeq™ Identity Panel in a southern Chinese population.

In the present study, 90 autosomal single nucleotide polymorphisms (SNPs) and 34 Y chromosomal SNPs were sequenced simultaneously using HID-Ion AmpliSeq™ Identity Panel on the Ion PGM™ platform for 125 samples in a southern Chinese population. Raw data were analyzed and forensic parameters were calculated. Haplogrouping concordance was also assessed using alternative methods based on Y-SNP haplotypes and Y-STR haplotypes. The results showed that allelic imbalance occurred more frequently with low coverage while several SNPs with high coverage were also observed with poor allelic balance, including rs214955, rs430046, rs7520386, rs876724, rs9171188, rs16981290, and rs2032631. Totally, 21,261 miscalled reads (0.28%) were observed. The rate of allele-specific miscalled reads (ASMRs) was higher than that of allele nonspecific miscalled reads (ANMRs) and associated with genetic diversity of the SNP. The ASMRs of major allele were lower than that of minor allele while there was no difference for ANMRs. The combined discrimination power (CDP) was 1-4.81 × 10-34 and the combined power of exclusion (CPE) was 0.99989 and 0.99999992 for duo and trio paternity testing, respectively. No significant genetic difference was detected between southern and northern Chinese populations. For haplogroup study, O2 was the predominant haplogroup and 97.01% of samples were assigned consistent haplogoups with Y-SNP and Y-STR haplotypes. In conclusion, the AmpliSeq™ Identity Panel was powerful for individual identification and trio paternity testing. ASMRs were associated with the genetic diversity and allele frequency while neither was related for ANMRs. High concordance of haplogrouping assignment can be obtained with Y-STR and Y-SNP haplotypes.

Paternity calculations in a di-spermy case.

In a criminal paternity case, which involved analysis of the product of conception, a rare circumstance was observed. The product of conception was triploidy, apparently due to an egg fertilized by two sperm. Since there is little guidance on how to calculate the probability of the DNA evidence given some basic hypotheses, the formulae were derived and are presented herein. These approaches could provide guidance for similar situations if they arise.

Comparison of preprocessing methods and storage times for touch DNA samples.

AIM: To select appropriate preprocessing methods for different substrates by comparing the effects of four different preprocessing methods on touch DNA samples and to determine the effect of various storage times on the results of touch DNA sample analysis. METHOD: Hand touch DNA samples were used to investigate the detection and inspection results of DNA on different substrates. Four preprocessing methods, including the direct cutting method, stubbing procedure, double swab technique, and vacuum cleaner method, were used in this study. DNA was extracted from mock samples with four different preprocessing methods. The best preprocess protocol determined from the study was further used to compare performance after various storage times. DNA extracted from all samples was quantified and amplified using standard procedures. RESULTS: The amounts of DNA and the number of alleles detected on the porous substrates were greater than those on the non-porous substrates. The performances of the four preprocessing methods varied with different substrates. The direct cutting method displayed advantages for porous substrates, and the vacuum cleaner method was advantageous for non-porous substrates. No significant degradation trend was observed as the storage times increased. CONCLUSION: Different substrates require the use of different preprocessing method in order to obtain the highest DNA amount and allele number from touch DNA samples. This study provides a theoretical basis for explorations of touch DNA samples and may be used as a reference when dealing with touch DNA samples in case work.

A convenient guideline to determine if two Y-STR profiles are from the same lineage.

Y chromosome STR loci are used in forensics primarily for identification purposes by determining the male lineages. The Henan province in China has established a large Y-STR (>200 000 profiles) database for criminal investigations. A large proportion of the Y-STR profiles in the database were generated using either the Applied Biosystems Yfiler(ۛ) or Yfiler(ۛ) Plus PCR Amplification kits. The additional loci in the Yfiler Plus kit as compared to the Yfiler kit results in a concomitant cumulative mutation rate increase across the loci. Therefore, in those cases when two profiles have one to a few mismatched loci, it is difficult to determine if they are from the same lineage. In this study, 7405 unrelated male profiles were manually selected from the database. Analysis showed higher power of discrimination than the corresponding Yfiler haplotypes. Further, the distributions of the number of mismatched loci and the mismatched steps were generated for father-son, grandfather-grandson, uncle-nephew, and cousins (i.e. one, two, three, and four meioses, respectively) by exhaustive pairwise comparison of the unrelated profiles using a dynamic programming approach. The same distributions were generated for unrelated pairs with mutation rates of the loci. With the distributions, the false negative and false positive rates were determined. Two Yfiler profiles with ≤2 mismatched loci or ≤2 steps are more likely from the same lineage than unrelated lineages, and two Yfiler Plus profiles with ≤4 mismatched loci or ≤5 mismatched steps are more likely from the same lineage.

Blind study evaluation illustrates utility of the Ion PGM™ system for use in human identity DNA typing.

AIM: To perform a blind study to assess the capability of the Ion Personal Genome Machine® (PGM™) system to sequence forensically relevant genetic marker panels and to characterize unknown individuals for ancestry and possible relatedness. METHODS: Twelve genomic samples were provided by a third party for blinded genetic analysis. For these 12 samples, the mitochondrial genome and three PGM™ panels containing human identity single nucleotide polymorphisms (SNPs), ancestry informative SNPs, and short tandem repeats (STRs) were sequenced on the PGM™ system and analyzed. RESULTS: All four genetic systems were run and analyzed on the PGM™ system in a reasonably quick time frame. Completeness of genetic profiles, depth of coverage, strand balance, and allele balance were informative metrics that illustrated the quality and reliability of the data produced. SNP genotypes allowed for identification of sex, paternal lineage, and population ancestry. STR genotypes were shown to be in complete concordance with genotypes generated by standard capillary electrophoresis-based technologies. Variants in the mitochondrial genome data provided information on population background and maternal relationships. CONCLUSION: All results from analysis of the 12 genomic samples were consistent with sample information provided by the sample providers at the end of the blinded study. The relatively easy identification of intra-STR allele SNPs offered the potential for increased discrimination power. The promising nature of these results warrants full validation studies of this massively parallel sequencing technology and its further development for forensic data analysis.

Modeling one complete versus triplicate analyses in low template DNA typing.

There are generally two strategies for low template DNA typing: the complete strategy, which uses all available DNA in a single PCR and subsequent typing, and the consensus strategy, in which the biological sample is divided into two or more aliquots and the genotype profile is determined by consensus from these "replicates." In this study, the consensus and complete strategies are compared by a statistical approach in terms of the accuracy of obtaining the correct genotype at a single locus for single source samples. Logistic models were employed to describe the allele drop-out and drop-in events. The parameters of the models were estimated with empirical or hypothetical data. The probabilities of obtaining the true genotype and the chances to observe drop-out and drop-in alleles were estimated and compared for both strategies. Consistent with a previous experimental study, this study found that, with relatively high input DNA (e.g., ≥ 100 pg), the complete strategy performs better than the consensus strategy to obtain the true genotype and the complete strategy will display less dropped out alleles. The consensus strategy had less drop-in alleles for ≤ 100 pg DNA samples. Moreover, the limitations of the logistic models were discussed. Ideal models with better fit of empirical data approximating casework conditions were proposed for future studies.

Population genetic analyses of the STR loci of the AmpFlSTR NGM SElect™ kit for Han population in Fujian Province, China.

Allele frequencies and forensically relevant population statistics of the STR loci in the AmpFlSTR® NGM SElect™ PCR Amplification Kit were estimated for the Han population from Fujian province in China (n = 454). All loci were highly polymorphic and the cumulative match probability was 5.4 × 10(-21). No significant departure from Hardy-Weinberg equilibrium and linkage equilibrium was detected after correction for sampling. The population substructure of Fujian Han population is minor.

Mutation rates at 16 Y-chromosome STRs in the South China Han population.

Nine hundred eighteen DNA-confirmed father-son pairs were typed for 16 Y-chromosomal short tandem repeat (Y-STR) markers by AmpFLSTR® Yfiler™ PCR Amplification kit. In a total of 15,606 allele transmissions, 36 mutations were detected. The average mutation rate across all 16 Y-STR markers was 0.0023 (95 % confidence interval, 0.0016-0.0032). One two-step mutation was found at DYS389II, and all other mutations were single steps. The losses and gains were balanced at all other loci, excluding DYS385 and DYS458, where losses were more frequent than gains. Mutation rates among different Y-STR loci were significantly different (Χ (2) = 69.05, P = 0.000). Mutation rates were correlated with the lengths of the alleles. Alleles with higher number of repeats were more likely to mutate. Mutation rates were also correlated with the gene diversity of the locus (r (2) = 0.565, P = 0.023). Loci with higher gene diversity had higher mutation rates. In addition, the mutation rate of the older father was found to be notably higher than that of the younger father.

Single nucleotide polymorphism typing with massively parallel sequencing for human identification.

The Ion AmpliSeq™ HID single nucleotide polymorphism (SNP) panel, a primer pool of 103 autosomal SNPs and 33 Y-SNPs, was evaluated using the Ion 314™ Chip on the Ion PGM™ Sequencer with four DNA samples. The study focused on the sequencing of DNA at three different initial target quantities, related interpretation issues, and concordance of results with another sequencing platform, i.e., Genome Analyzer IIx. With 10 ng of template DNA, all genotypes at the 136 SNPs were detected. With 1 ng of DNA, all SNPs were detected and one SNP locus in one sample showed extreme heterozygote imbalance on allele coverage. With 100 pg of DNA, an average of 1.6 SNP loci were not detected, and an average of 4.3 SNPs showed heterozygote imbalance. The average sequence coverage was 945-600× at autosomal SNPs and 465-209× at Y-SNPs for 10 ng-100 pg of DNA. The average heterozygote allele coverage ratio was 89.6-61.8 % for 10 ng-100 pg of DNA. At 10 ng of DNA, all genotypes of the 95 SNPs shared between the two different sequencing platforms were concordant except for one SNP, rs1029047. The error was due to the misalignment of a flanking homopolymer. Overall, the data support that genotyping a large battery of SNPs is feasible with massively parallel sequencing. With barcode systems, better allele balance, and specifically designed alignment software, a more comprehensive rapid genotyping and more cost-effective results may be obtained from multiple samples in one analysis than are possible with current typing and capillary electrophoresis systems.

INNULs: A novel design amplification strategy for retrotransposable elements for studying population variation.

OBJECTIVES: Retrotransposable elements (REs), consisting of long interspersed nuclear elements (LINEs) and short interspersed nuclear elements (SINEs), are a group of markers that can be useful for human identity testing. Until now, however, due to the inherent size difference (up to 6 kb in some instances) associated with insertion and null alleles (or INNULs), the use of REs for facilitated population studies has not been sought or practical. The size of the insertion elements (from a few hundred to several thousand bp) has proven to limit their utility as a marker because of the inefficient amplicon yield with PCR. A novel primer design now facilitates INNUL marker testing. A preliminary panel of single-locus markers was developed to evaluate the potential of typing these insertion elements. Nine INNULs (5 Alu and 4 LINEs) were typed in three major North American populations and analyzed for population genetic features. In addition, the variation of each marker among the sample populations provides insight of its potential use as individual identification or ancestral marker. METHODS: INNUL markers were developed into fluorescently labeled single-loci PCR. Nine markers were developed with amplicons that were less than 180 bp in length, and, depending on the locus amplicons of the INNULs, alleles varied in size from 50 to 1 bp. This allele size is noteworthy because the insertion alleles of the 9 loci range in size from 297 to 6,195 bp. The allele distribution of the INNULs was assessed and analyzed in three major North American populations. RESULTS: Upon observation of the distribution of the alleles in three major North American populations, the markers generally met Hardy-Weinberg expectations, and there was little evidence of detectable levels of linkage disequilibrium. Due to varying distributions of the alleles in the major population groups tested, some of the markers might be better suited for use as an individual identification marker, while others are better suited for bio-ancestral studies. CONCLUSIONS: Using the primer design strategy described in our work, SINEs and (for the first time, to our knowledge) LINEs can be utilized as markers for studying population genetic variation that is more amenable to the limitations of the PCR technique. This study lays the foundation for future work of developing a multiplex panel of INNUL markers that can be used as a single-tube assay for human identity testing utilizing small amplicons (<180 bp), which could be useful for ancient or degraded forensic DNA samples.

[Development of Chinese forensic Y-STR DNA database].

Y chromosome is a male-specific paternal inherited chromosome. The STR markers on Y chromosome have been widely used in forensic practices. This article summarizes the characteristics of Y-STR and some factors are considered of selecting appropriate Y-STR markers for Chinese population. The prospects of existing and potential forensic applications of Y-STR profiles are discussed including familial excluding, familial searching, crowd source deducing, mixture sample testing, and kinship identifying. The research, development, verification of Y-STR kit, Y-STR mutation rate, and search software are explored and some suggestions are given.