From next generation sequencing to now generation sequencing in forensics

Solving the twin paradox-forensic strategies to identify the identical twins

Identical twins are also called monozygotic twins which originate from the same zygote that possesses the same genetic make-up. To discriminate between identical monozygotic twins, short tandem repeats has not been found effective, therefore, various techniques, including next-generation sequencing (NGS), are applied. Monozygotic twins can be identified through germ line genomes, through speech using deep learning networks, and through epigenetic analysis. Fingerprint analysis has also been used to distinguish between identical twins, as human beings have unique fingerprints. Two distinct levels of fingerprint are used to distinguish between monozygotic twins based upon the differences in the minutiae points. Examination of the methylation pattern of the genome has an enormous potential to differentiate between identical twins, as the methylation of DNA occurs uniquely to each individual. This article offers an insight into the latest methods and techniques used for the differentiation between the identical twins.

Applications and Performance of Precision ID GlobalFiler NGS STR, Identity, and Ancestry Panels in Forensic Genetics

Short Tandem Repeat (STR) testing via capillary electrophoresis is undoubtedly the most popular forensic genetic testing method. However, its low multiplexing capabilities and limited performance with challenging samples are among the factors pushing scientists towards new technologies. Next-generation sequencing (NGS) methods overcome some of these limitations while also enabling the testing of Single-Nucleotide Polymorphisms (SNPs). Nonetheless, these methods are still under optimization, and their adoption into practice is limited. Among the available kits, Thermo Fisher Scientific (Waltham, MA, USA) produces three Precision ID Panels: GlobalFiler NGS STR, Identity, and Ancestry. A clear review of these kits, providing information useful for the promotion of their use, is, however, lacking. To close the gap, a literature review was performed to investigate the popularity, applications, and performance of these kits. Following the PRISMA guidelines, 89 publications produced since 2015 were identified. China was the most active country in the field, and the Identity Panel was the most researched. All kits appeared robust and useful for low-quality and low-quantity samples, while performance with mixtures varied. The need for more population data was highlighted, as well as further research surrounding variables affecting the quality of the sequencing results.

An exploratory view into allelic drop-out of sequenced autosomal STRs

As massively parallel sequencing is implemented in forensic genetics, an understanding of sequence data must accompany these advancements, that is, accurate modeling of data for proper statistical analysis. Allelic drop-out, a common stochastic effect seen in genetic data, is often modeled in statistical analysis of STR results. This proof-of-concept study sequenced several serial dilutions of a standard sample ranging from 4 ng to 7.82 pg to evaluate allelic drop-out trends on a select panel of autosomal STRs using the ForenSeq™ DNA Signature Prep Kit, Primer Set A on the Illumina MiSeq FGx. Parameters assessed included locus, profile, and run specific information. A majority of the allelic drop-out occurred in DNA concentrations less than 31.25 pg. Statistical results indicated a need for locus-specific modeling based on STR descriptors, like simple versus compound repeat patterns. No correlation was seen between average read count of scored alleles and allelic drop-out at a locus. A statistical correlation was observed between the amount of allelic drop-out and the starting amount of DNA in a sample, average read count of a sample, and total read count generated on a flow cell. This study supports using common allelic drop-out factors used in fragment length analysis on sequenced STRs while including additional locus, sample, and run specific information. Results demonstrate multiple factors that can be considered when developing probability of allelic drop-out models for sequenced autosomal STRs including locus-specific analysis, total read count of a profile, and total read count sequenced on a flow cell.

A novel 193-plex MPS panel integrating STRs and SNPs highlights the application value of forensic genetics in individual identification and paternity testing

Massively parallel sequencing (MPS) has emerged as a promising technology for targeting multiple genetic loci simultaneously in forensic genetics. Here, a novel 193-plex panel was designed to target 28 A-STRs, 41 Y-STRs, 21 X-STRs, 3 sex-identified loci, and 100 A-SNPs by employing a single-end 400 bp sequencing strategy on the MGISEQ-2000™ platform. In the present study, a series of validations and sequencing of 1642 population samples were performed to evaluate the overall performance of the MPS-based panel and its practicality in forensic application according to the SWGDAM guidelines. In general, the 193-plex markers in our panel showed good performance in terms of species specificity, stability, and repeatability. Compared to commercial kits, this panel achieved 100% concordance for standard gDNA and 99.87% concordance for 14,560 population genotypes. Moreover, this panel detected 100% of the loci from 0.5 ng of DNA template and all unique alleles at a 1:4 DNA mixture ratio (0.2 ng minor contributor), and the applicability of the proposed approach for tracing and degrading DNA was further supported by case samples. In addition, several forensic parameters of STRs and SNPs were calculated in a population study. High CPE and CPD values greater than 0.9999999 were clearly demonstrated and these results could be useful references for the application of this panel in individual identification and paternity testing. Overall, this 193-plex MPS panel has been shown to be a reliable, repeatable, robust, inexpensive, and powerful tool sufficient for forensic practice.

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.

Perceptions of forensic scientists on statistical models, sequence data, and ethical implications for DNA evidence evaluations: A qualitative assessment

With the introduction of next generation sequencing (NGS) technology in the forensic field, it will be of interest to assess if forensic scientists feel equipped to interpret and present DNA evidence for sequence data. Here, we describe perceptions of sixteen U.S.-based forensic scientists on statistical models, sequence data, and ethical implications for DNA evidence evaluations. To get an in-depth understanding of the current situation, we used a qualitative research approach with a cross-sectional study design. Semi-structured interviews (N = 16) were conducted with U.S. forensic scientists working with DNA evidence. Open-ended interview questions were used to explore participants' views and needs surrounding the use of statistical models and sequence data for forensic purposes. We conducted a conventional content analysis using ATLAS. ti software and employed a second coder to ensure reliability of our results. Eleven themes emerged: 1) a statistical model that maximizes the value of the evidence is preferred; 2) a high-level understanding of the statistical model used is generally sufficient; 3) transparency is key in minimizing the risk of creating black boxes; 4) training and education should be an ongoing effort; 5) the effectiveness of presenting results in court can be improved; 6) NGS has the potential to become revolutionary; 7) some hesitations surrounding the use of sequence data remain; 8) there is a need for a concrete plan to alleviate barriers to the implementation of sequencing techniques; 9) ethics plays a major part in the role of a forensic scientist; 10) ethical barriers for sequence data depend on the application; 11) DNA evidence has its limitations. The results of this study give insight into the perceptions of forensic scientists regarding the use of statistical models and sequence data, providing valuable information in the move towards implementing sequencing methods for DNA evidence evaluations.

Accuracy of Eye and Hair Color Prediction in Mexican Mestizos from Monterrey City Based on ForenSeqTM DNA Signature Prep

Forensic genomic systems allow simultaneously analyzing identity informative (iiSNPs), ancestry informative (aiSNPs), and phenotype informative (piSNPs) genetic markers. Among these kits, the ForenSeq DNA Signature prep (Verogen) analyzes identity STRs and SNPs as well as 24 piSNPs from the HIrisPlex system to predict the hair and eye color. We report herein these 24 piSNPs in 88 samples from Monterrey City (Northeast, Mexico) based on the ForenSeq DNA Signature prep. Phenotypes were predicted by genotype results with both Universal Analysis Software (UAS) and the web tool of the Erasmus Medical Center (EMC). We observed predominantly brown eyes (96.5%) and black hair (75%) phenotypes, whereas blue eyes, and blond and red hair were not observed. Both UAS and EMC showed high performance in eye color prediction (p ≥ 96.6%), but a lower accuracy was observed for hair color prediction. Overall, UAS hair color predictions showed better performance and robustness than those obtained with the EMC web tool (when hair shade is excluded). Although we employed a threshold (p > 70%), we suggest using the EMC enhanced approach to avoid the exclusion of a high number of samples. Finally, although our results are helpful to employ these genomic tools to predict eye color, caution is suggested for hair color prediction in Latin American (admixed) populations such as those studied herein, principally when no black color is predicted.

Targeted DNA methylation analysis and prediction of smoking habits in blood based on massively parallel sequencing

Tobacco smoking is a frequent habit sustained by > 1.3 billion people in 2020 and the leading preventable factor for health risk and premature mortality worldwide. In the forensic context, predicting smoking habits from biological samples may allow broadening DNA phenotyping. In this study, we aimed to implement previously published smoking habit classification models based on blood DNA methylation at 13 CpGs. First, we developed a matching lab tool based on bisulfite conversion and multiplex PCR followed by amplification-free library preparation and targeted paired-end massively parallel sequencing (MPS). Analysis of six technical duplicates revealed high reproducibility of methylation measurements (Pearson correlation of 0.983). Artificially methylated standards uncovered marker-specific amplification bias, which we corrected via bi-exponential models. We then applied our MPS tool to 232 blood samples from Europeans of a wide age range, of which 90 were current, 71 former and 71 never smokers. On average, we obtained 189,000 reads/sample and 15,000 reads/CpG, without marker drop-out. Methylation distributions per smoking category roughly corresponded to previous microarray analysis, showcasing large inter-individual variation but with technology-driven bias. Methylation at 11 out of 13 smoking-CpGs correlated with daily cigarettes in current smokers, while solely one was weakly correlated with time since cessation in former smokers. Interestingly, eight smoking-CpGs correlated with age, and one displayed weak but significant sex-associated methylation differences. Using bias-uncorrected MPS data, smoking habits were relatively accurately predicted using both two- (current/non-current) and three- (never/former/current) category model, but bias correction resulted in worse prediction performance for both models. Finally, to account for technology-driven variation, we built new, joint models with inter-technology corrections, which resulted in improved prediction results for both models, with or without PCR bias correction (e.g. MPS cross-validation F₁-score > 0.8; 2-categories). Overall, our novel assay takes us one step closer towards the forensic application of viable smoking habit prediction from blood traces. However, future research is needed towards forensically validating the assay, especially in terms of sensitivity. We also need to further shed light on the employed biomarkers, particularly on the mechanistics, tissue specificity and putative confounders of smoking epigenetic signatures.

Sequence-based allelic variations and frequencies for 22 autosomal STR loci in the Lebanese population

This is the first study that characterizes the sequence-based allelic variations of 22 autosomal Short Tandem Repeat (aSTR) loci in a population dataset collected from Lebanon. Genomic DNA extracts from 195 unrelated Lebanese individuals were amplified with PowerSeq 46GY System Prototype. Targeted amplicons were subjected to DNA library preparation and sequenced on the Verogen MiSeq FGx Sequencing System. Raw FASTQ data files were processed by STRait Razor v3. Sequence strings were annotated according to the considerations of the DNA Commission of the International Society for Forensic Genetics (ISFG) and tabulated herein with their respective allelic frequencies and GeneBank accession and version numbers. The sequenced Lebanese dataset resulted in 429 distinct allelic sequences as compared to the 236 alleles identified by length only. The increase in the number of alleles was observed at 18 out of 22 aSTR loci and was attributed to the sequence variations residing in both the STR repeat motifs and flanking regions. The study uncovered 25 novel aSTR allelic sequences across 12 loci for which GenBank records did not previously exist in the STRSeq BioProject, PRJNA380127. For a concordance check, the length-based allelic calls derived from the full sequences were compared to those genotyped using capillary electrophoresis (CE) methods. Population genetic parameters relevant to the evaluation of forensic DNA evidence were assessed for the sequence-based data and compared to the parameters generated from the length-based information. Using the sequence-based data, Analysis of MOlecular VAriance (AMOVA), genetic distances, and population genetic structure were evaluated for 1231 individuals sampled from the Lebanese and four U.S. populations (African American, Asian, Caucasian, and Hispanic). The results were tabulated and visualized in a population tree, multidimensional scaling scatter plots, and bar plots. This newly established sequence-based database for the Lebanese population can be beneficial for extending NGS applicability to casework or paternity testing and assessing the strength of evidence for NGS-STR profiles. The described novel sequence variants at certain loci can further help in the effort to characterize the sequence diversity of STR markers from different populations around the world.

Applying Unique Molecular Indices with an Extensive All-in-One Forensic SNP Panel for Improved Genotype Accuracy and Sensitivity

One of the major challenges in forensic genetics is being able to detect very small amounts of DNA. Massively parallel sequencing (MPS) enables sensitive detection; however, genotype errors may exist and could interfere with the interpretation. Common errors in MPS-based analysis are often induced during PCR or sequencing. Unique molecular indices (UMIs) are short random nucleotide sequences ligated to each template molecule prior to amplification. Applying UMIs can improve the limit of detection by enabling accurate counting of initial template molecules and removal of erroneous data. In this study, we applied the FORCE panel, which includes ~5500 SNPs, with a QIAseq Targeted DNA Custom Panel (Qiagen), including UMIs. Our main objective was to investigate whether UMIs can enhance the sensitivity and accuracy of forensic genotyping and to evaluate the overall assay performance. We analyzed the data both with and without the UMI information, and the results showed that both genotype accuracy and sensitivity were improved when applying UMIs. The results showed very high genotype accuracies (>99%) for both reference DNA and challenging samples, down to 125 pg. To conclude, we show successful assay performance for several forensic applications and improvements in forensic genotyping when applying UMIs.

A global snapshot of current opinions of next-generation sequencing technologies usage in forensics

The future of forensic DNA testing is being shaped by the research and usage of next-generation systems, which have increased the multiplexing capabilities of the field and the type and amount of genetic data that can be utilized for investigations. The NGS adoption for casework has been slow, albeit the plethora of data that has been published. This study evaluated the current opinions on sequencing in forensics. A 20-question online-survey focusing on NGS knowledge, training, and usage was distributed to 6001 forensic DNA researchers and practitioners worldwide. A total of 367 responses were obtained from all continents (North/South America (69.8%), Europe (21.2%), Asia (5.5%), Oceania (2.5%), and Africa (1%)). The respondents consisted of 50% practitioners, 31% researchers, and 19% both. Of these, 38% already own a next-gen sequencing instrument, and 13% are planning to purchase one. Overall, there exists an extensive knowledge on next-gen sequencing within the forensic community, including among laboratories that have not yet implemented this high-throughput technology in their workflows. Current usage focuses primarily on SNP analysis for investigative leads and mitochondrial DNA analysis while future applications included both STR and SNP testing applied to general casework. The major overall concerns respondents have for implementing a sequencing instrument include limited funding, staffing, lack of time, and the cost-effectiveness of providing this service. Specific technical concerns that the respondents had are the lack of training, statistical applications, bioinformatics support, and of rigorous guidelines and recommendations. Most of the respondents do believe there will be a technology shift from using CE only to the use of NGS on casework in 5-10 years. In addition, around 66% of respondents believe that it is moderately to very likely that the court will accept sequencing analysis. Sixteen percent fell in the middle, and the remaining 15% believe it is more unlikely, with 3% of respondents believing it is very unlikely. In conclusion, this work outlines current analytical challenges experienced by the global forensic DNA community and addresses different strategies for the implementation of next-gen sequencing technologies in casework.

Identification and characterization of novel DIP-STRs from whole-genome sequencing data

In an attempt to enhance forensic DNA mixture deconvolution several alternative DNA typing approaches have been developed. Among these, DIP-STR compound markers can resolve extremely unbalanced two-source DNA mixtures of same-or-opposite sex donors, up to a 1:1000 minor:major DNA ratio. A forensic set of 10 markers was validated for casework and a larger set of 23 DIP-STRs has proven suitable to biogeographic ancestry inference and for prenatal paternity testing. Yet, to promote the widespread use of this original approach, more markers and multiplex panels need to be developed. To this end, here we describe an extended set of forensic DIP-STRs identified using currently available whole-genome sequencing datasets. Complete lists of Indels and STRs were obtained from reported frequencies of genetic variants of 76,156 genomes. About 3000 identified DIP-STRs candidates were shorter than 200 bp and 500 showed high haplotype variability estimated using the genotypes of individuals homozygous for the DIP or the STR. Here, we present 23 additional DIP-STRs validated for sensitivity, specificity and Swiss population variability. Finally, a set of 30 markers comprising seven previously validated ones is proposed for the prospective development of a forensic DIP-STR multiplex panel.

Developmental validation of STRmix™ NGS, a probabilistic genotyping tool for the interpretation of autosomal STRs from forensic profiles generated using NGS

We describe the developmental validation of the probabilistic genotyping software - STRmix™ NGS - developed for the interpretation of forensic DNA profiles containing autosomal STRs generated using next generation sequencing (NGS) also known as massively parallel sequencing (MPS) technologies. Developmental validation was carried out in accordance with the Scientific Working Group on DNA Analysis Methods (SWGDAM) Guidelines for the Validation of Probabilistic Genotyping Systems and the International Society for Forensic Genetics (ISFG) recommendations and included sensitivity and specificity testing, accuracy, precision, and the interpretation of case-types samples. The results of developmental validation demonstrate the appropriateness of the software for the interpretation of profiles developed using NGS technology.

Recent advances in forensic biology and forensic DNA typing: INTERPOL review 2019-2022

This review paper covers the forensic-relevant literature in biological sciences from 2019 to 2022 as a part of the 20th INTERPOL International Forensic Science Managers Symposium. Topics reviewed include rapid DNA testing, using law enforcement DNA databases plus investigative genetic genealogy DNA databases along with privacy/ethical issues, forensic biology and body fluid identification, DNA extraction and typing methods, mixture interpretation involving probabilistic genotyping software (PGS), DNA transfer and activity-level evaluations, next-generation sequencing (NGS), DNA phenotyping, lineage markers (Y-chromosome, mitochondrial DNA, X-chromosome), new markers and approaches (microhaplotypes, proteomics, and microbial DNA), kinship analysis and human identification with disaster victim identification (DVI), and non-human DNA testing including wildlife forensics. Available books and review articles are summarized as well as 70 guidance documents to assist in quality control that were published in the past three years by various groups within the United States and around the world.

BigFiRSt: A Software Program Using Big Data Technique for Mining Simple Sequence Repeats From Large-Scale Sequencing Data

Background

Simple Sequence Repeats (SSRs) are short tandem repeats of nucleotide sequences. It has been shown that SSRs are associated with human diseases and are of medical relevance. Accordingly, a variety of computational methods have been proposed to mine SSRs from genomes. Conventional methods rely on a high-quality complete genome to identify SSRs. However, the sequenced genome often misses several highly repetitive regions. Moreover, many non-model species have no entire genomes. With the recent advances of next-generation sequencing (NGS) techniques, large-scale sequence reads for any species can be rapidly generated using NGS. In this context, a number of methods have been proposed to identify thousands of SSR loci within large amounts of reads for non-model species. While the most commonly used NGS platforms (e.g., Illumina platform) on the market generally provide short paired-end reads, merging overlapping paired-end reads has become a common way prior to the identification of SSR loci. This has posed a big data analysis challenge for traditional stand-alone tools to merge short read pairs and identify SSRs from large-scale data.

Results

In this study, we present a new Hadoop-based software program, termed BigFiRSt, to address this problem using cutting-edge big data technology. BigFiRSt consists of two major modules, BigFLASH and BigPERF, implemented based on two state-of-the-art stand-alone tools, FLASH and PERF, respectively. BigFLASH and BigPERF address the problem of merging short read pairs and mining SSRs in the big data manner, respectively. Comprehensive benchmarking experiments show that BigFiRSt can dramatically reduce the execution times of fast read pairs merging and SSRs mining from very large-scale DNA sequence data.

Conclusions

The excellent performance of BigFiRSt mainly resorts to the Big Data Hadoop technology to merge read pairs and mine SSRs in parallel and distributed computing on clusters. We anticipate BigFiRSt will be a valuable tool in the coming biological Big Data era.

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.

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

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

Developments in forensic DNA analysis

The analysis of DNA from biological evidence recovered in the course of criminal investigations can provide very powerful evidence when a recovered profile matches one found on a DNA database or generated from a suspect. However, when no profile match is found, when the amount of DNA in a sample is too low, or the DNA too degraded to be analysed, traditional STR profiling may be of limited value. The rapidly expanding field of forensic genetics has introduced various novel methodologies that enable the analysis of challenging forensic samples, and that can generate intelligence about the donor of a biological sample. This article reviews some of the most important recent advances in the field, including the application of massively parallel sequencing to the analysis of STRs and other marker types, advancements in DNA mixture interpretation, particularly the use of probabilistic genotyping methods, the profiling of different RNA types for the identification of body fluids, the interrogation of SNP markers for predicting forensically relevant phenotypes, epigenetics and the analysis of DNA methylation to determine tissue type and estimate age, and the emerging field of forensic genetic genealogy. A key challenge will be for researchers to consider carefully how these innovations can be implemented into forensic practice to ensure their potential benefits are maximised.

CSYseq: The first Y-chromosome sequencing tool typing a large number of Y-SNPs and Y-STRs to unravel worldwide human population genetics

Male-specific Y-chromosome (chrY) polymorphisms are interesting components of the DNA for population genetics. While single nucleotide polymorphisms (Y-SNPs) indicate distant evolutionary ancestry, short tandem repeats (Y-STRs) are able to identify close familial kinships. Detailed chrY analysis provides thus both biogeographical background information as paternal lineage identification. The rapid advancement of high-throughput massive parallel sequencing (MPS) technology in the past decade has revolutionized genetic research. Using MPS, single-base information of both Y-SNPs as Y-STRs can be analyzed in a single assay typing multiple samples at once. In this study, we present the first extensive chrY-specific targeted resequencing panel, the 'CSYseq', which simultaneously identifies slow mutating Y-SNPs as evolution markers and rapid mutating Y-STRs as patrilineage markers. The panel was validated by paired-end sequencing of 130 males, distributed over 65 deep-rooted pedigrees covering 1,279 generations. The CSYseq successfully targets 15,611 Y-SNPs including 9,014 phylogenetic informative Y-SNPs to identify 1,443 human evolutionary Y-subhaplogroup lineages worldwide. In addition, the CSYseq properly targets 202 Y-STRs, including 81 slow, 68 moderate, 27 fast and 26 rapid mutating Y-STRs to individualize close paternal relatives. The targeted chrY markers cover a high average number of reads (Y-SNP = 717, Y-STR = 150), easy interpretation, powerful discrimination capacity and chrY specificity. The CSYseq is interesting for research on different time scales: to identify evolutionary ancestry, to find distant family and to discriminate closely related males. Therefore, this panel serves as a unique tool valuable for a wide range of genetic-genealogical applications in interdisciplinary research within evolutionary, population, molecular, medical and forensic genetics.

Massively parallel sequencing and capillary electrophoresis of a novel panel of falcon STRs: Concordance with minisatellite DNA profiles from historical wildlife crime

Birds of prey have suffered persecution for centuries through trapping, shooting, poisoning and theft from the wild to meet the demand from egg collectors and falconers; they were also amongst the earliest beneficiaries of DNA testing in wildlife forensics. Here we report the identification and characterisation of 14 novel tetramer, pentamer and hexamer short tandem repeat (STR) markers which can be typed either by capillary electrophoresis or massively parallel sequencing (MPS) and apply them to historical casework samples involving 49 peregrine falcons, 30 of which were claimed to be the captively bred offspring of nine pairs. The birds were initially tested in 1994 with a multilocus DNA fingerprinting probe, a sex test and eight single-locus minisatellite probes (SLPs) demonstrating that 23 birds were unrelated to the claimed parents. The multilocus and SLP approaches were highly discriminating but extremely time consuming and required microgram quantities of high molecular weight DNA and the use of radioisotopes. The STR markers displayed between 2 and 21 alleles per locus (mean = 7.6), lengths between 140 and 360 bp, and heterozygosities from 0.4 to 0.93. They produced wholly concordant conclusions with similar discrimination power but in a fraction of the time using a hundred-fold less DNA and with standard forensic equipment. Furthermore, eleven of these STRs were amplified in a single reaction and typed using MPS on the Illumina MiSeq platform revealing eight additional alleles (three with variant repeat structures and five solely due to flanking SNPs) across four loci. This approach gave a random match probability of < 1E-9, and a parental pair false inclusion probability of < 1E-5, with a further ten-fold reduction in the amount of DNA required (~3 ng) and the potential to analyse mixed samples. These STRs will be of value in monitoring wild populations of these key indicator species as well as for testing captive breeding claims and establishing a database of captive raptors. They have the potential to resolve complex cases involving trace, mixed and degraded samples from raptor persecution casework representing a significant advance over the previously applied methods.

A worldwide map of swine short tandem repeats and their associations with evolutionary and environmental adaptations

Background

Short tandem repeats (STRs) are genetic markers with a greater mutation rate than single nucleotide polymorphisms (SNPs) and are widely used in genetic studies and forensics. However, most studies in pigs have focused only on SNPs or on a limited number of STRs.

Results

This study screened 394 deep-sequenced genomes from 22 domesticated pig breeds/populations worldwide, wild boars from both Europe and Asia, and numerous outgroup Suidaes, and identified a set of 878,967 polymorphic STRs (pSTRs), which represents the largest repository of pSTRs in pigs to date. We found multiple lines of evidence that pSTRs in coding regions were affected by purifying selection. The enrichment of trinucleotide pSTRs in coding sequences (CDS), 5′UTR and H3K4me3 regions suggests that trinucleotide STRs serve as important components in the exons and promoters of the corresponding genes. We demonstrated that, compared to SNPs, pSTRs provide comparable or even greater accuracy in determining the breed identity of individuals. We identified pSTRs that showed significant population differentiation between domestic pigs and wild boars in Asia and Europe. We also observed that some pSTRs were significantly associated with environmental variables, such as average annual temperature or altitude of the originating sites of Chinese indigenous breeds, among which we identified loss-of-function and/or expanded STRs overlapping with genes such as AHR, LAS1L and PDK1. Finally, our results revealed that several pSTRs show stronger signals in domestic pig—wild boar differentiation or association with the analysed environmental variables than the flanking SNPs within a 100-kb window.

Conclusions

This study provides a genome-wide high-density map of pSTRs in diverse pig populations based on genome sequencing data, enabling a more comprehensive characterization of their roles in evolutionary and environmental adaptation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00631-4.

Progress in the implementation of massively parallel sequencing for forensic genetics: results of a European-wide survey among professional users

A European-wide online survey was conducted to generate an overview on the state-of-the-art using massively parallel sequencing (MPS) platforms for forensic DNA analysis and DNA phenotyping among forensic practitioners in Europe. The survey was part of the dissemination activities of the “VISible Attributes through GEnomics – VISAGE” Horizon 2020 funded European research project [30], in preparation of a series of educational training activities. A total of 105 replies from 32 European countries representing participants from police, governmental, academic, and private laboratories providing professional services in the field of forensic genetics were included in the final analysis. Of these, 73% already own an MPS platform or are planning to acquire one within the next 1–2 years. One-third of the participants have already carried out MPS-based STR sequencing, identity, or ancestry SNP typing. A total of 23–40% of participants are planning to explore all FDP applications showing the overall very high interest in using MPS for the whole range of forensic MPS markers and applications. About 50% of the participants have previously gathered experience using forensic DNA phenotyping (FDP) markers based on conventional (i.e., not MPS-based) DNA typing methods. A total of 55% of the participants have attended training on the general use of MPS technology, but 36% have received no training whatsoever. Accordingly, 90% have expressed high or medium interest to attend training on the analysis and interpretation of DNA phenotyping data for predicting appearance, ancestry, and age. The results of our survey will provide valuable information for organizing relevant training workshops on all aspects of MPS-based DNA phenotyping for the forensic genetics scientific community.

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.

Parallel sequencing of 87 STR and 294 SNP markers using the prototype of the SifaMPS panel on the MiSeq FGx™ system

Massively parallel sequencing (MPS), or next generation sequencing (NGS), is a promising methodology for the detection of short tandem repeats (STRs) and single nucleotide polymorphisms (SNPs) in forensic genetics. Here, the prototype SifaMPS Panel is designed to simultaneously target 87 STRs and 294 SNPs with forensic interest in a single multiplex in conjunction with the TruSeq™ Custom Amplicon workflow and MiSeq FGx™ System. Two in-house python scripts are adopted for the fastq-to-genotype interpretation of MPS data concerning STR and SNP, respectively. In the present study, by sequencing 50 Chinese Hans and many other DNA samples involved in validation studies, system parameters including the depth of coverage (DoC), heterozygote balance (Hb) and sequence coverage ratios (SCRs), as well as different forensic parameters of STRs and SNPs in a population study, were calculated to evaluate the overall performance of this new panel and its practicality in forensic application. In general, except for two STRs (DYS505 and DYS449) and one SNP (rs4288409) that performed poorly, the other 85 STRs and 293 SNPs in our panel had good performance that could strengthen efficiency for human identification and paternity testing. In addition, discordant STR genotype results between those generated from capillary electrophoresis (CE) and from the MPS platform were clearly illustrated, and these results could be a useful reference for applying these particular non-CODIS STRs in forensic practice.

Application of a probabilistic genotyping software to MPS mixture STR data is supported by similar trends in LRs compared with CE data

The interpretation of short tandem repeat (STR) profiles can be challenging when, for example, alleles are masked due to allele sharing among contributors and/or when they are subject to drop-out, for instance from sample degradation. Mixture interpretation can be improved by increasing the number of STRs and/or loci with a higher discriminatory power. Both capillary electrophoresis (CE, 6-dye) and massively parallel sequencing (MPS) provide a platform for analysing relatively large numbers of autosomal STRs. In addition, MPS enables distinguishing between sequence variants, resulting in enlarged discriminatory power. Also, MPS allows for small amplicon sizes for all loci as spacing is not an issue, which is beneficial with degraded DNA. Altogether, MPS has the potential to increase the weights of evidence for true contributors to (complex) DNA profiles. In this study, likelihood ratio (LR) calculations were performed using STR profiles obtained with two different MPS systems and analysed using different settings: 1) MPS PowerSeq™ Auto System profiles analysed using FDSTools equipped with optimized settings such as noise correction, 2) ForenSeq™ DNA Signature Prep Kit profiles analysed using the default settings in the Universal Analysis Software (UAS), and 3) ForenSeq™ DNA Signature Prep Kit profiles analysed using FDSTools empirically adapted to cope with one-directional reads and provisional, basic settings. The LR calculations used genotyping data for two- to four-person mixtures varying for mixture proportion, level of drop-out and allele sharing and were generated with the continuous model EuroForMix. The LR results for the over 2000 sets of propositions were affected by the variation for the number of markers and analysis settings used in the three approaches. Nevertheless, trends for true and non-contributors, effects of replicates, assigned number of contributors, and model validation results were comparable for the three MPS approaches and alike the trends known for CE data. Based on this analogy, we regard the probabilistic interpretation of MPS STR data fit for forensic DNA casework. In addition, guidelines were derived on when to apply LR calculations to MPS autosomal STR data and report the corresponding results.

Challenges with co-amplification of microbial DNA in interpretation of STR profiles obtained from human skeletal remains

STR-based DNA analysis is still the main tool for human DNA identification in most forensic DNA laboratories. DNA typing of aged human skeletal elements faces well-known interpretation challenges characteristic of degraded and low copy number DNA samples. Analyzing tens of thousands of human bone and teeth samples, we found that the occasional presence of artefactual peaks of presumed microbial origin adds another layer of complexity to the interpretation of STR profiles. In this paper, we present our approach and suggest guidelines for identifying and distinguishing non-human peaks, developed over the last 18 years. Additionally, we report a compendium of artefact peaks of presumed microbial origin recorded in human STR profiles obtained from bone and teeth samples, originating from Iraq, Chile, Maldives, Brazil and Western Balkans. Our experience has shown that these artefacts are not uncommon in bone STR testing, suggesting the possibility of occurrence in other forensic contexts, particularly trace DNA samples. Raising awareness among the forensic DNA community and accounting for this phenomenon is important for accurate STR interpretation.

Monozygotic twins: Identical or distinguishable for science and law?

Monozygotic twins, also known as monovular twins, share an identical genetic heritage because they are two individuals who derive from the same zygote. For this reason, they have been considered indistinguishable. They represent a limit for the application of markers and analytical methods that are routinely used in forensic science because analyses of DNA fragments (short tandem repeats analysed by capillary electrophoresis) are unable to distinguish monozygotic twins. The recent introduction of ultra-deep next generation sequencing in forensic genetics, also known as massively parallel sequencing, has made it possible to identify a number of genetic variations through genome sequencing (such as copy number variations, single nucleotide polymorphisms and DNA methylation) that make it possible to distinguish monozygotic twins. Here, we present a case of ascertaining biological paternity, in which the alleged father had a monozygotic twin brother. This case led to the examination of international law in similar cases in which the only available biological evidence derives from classical forensic genetic analysis, performed with short tandem repeat (autosomal and/or gonosomal) capillary electrophoresis and the probative value, if recognised, of the next generation sequencing technology in the courtroom.

Blood group typing from whole-genome sequencing data

Many questions can be explored thanks to whole-genome data. The aim of this study was to overcome their main limits, software availability and database accuracy, and estimate the feasibility of red blood cell (RBC) antigen typing from whole-genome sequencing (WGS) data. We analyzed whole-genome data from 79 individuals for HLA-DRB1 and 9 RBC antigens. Whole-genome sequencing data was analyzed with software allowing phasing of variable positions to define alleles or haplotypes and validated for HLA typing from next-generation sequencing data. A dedicated database was set up with 1648 variable positions analyzed in KEL (KEL), ACKR1 (FY), SLC14A1 (JK), ACHE (YT), ART4 (DO), AQP1 (CO), CD44 (IN), SLC4A1 (DI) and ICAM4 (LW). Whole-genome sequencing typing was compared to that previously obtained by amplicon-based monoallelic sequencing and by SNaPshot analysis. Whole-genome sequencing data were also explored for other alleles. Our results showed 93% of concordance for blood group polymorphisms and 91% for HLA-DRB1. Incorrect typing and unresolved results confirm that WGS should be considered reliable with read depths strictly above 15x. Our results supported that RBC antigen typing from WGS is feasible but requires improvements in read depth for SNV polymorphisms typing accuracy. We also showed the potential for WGS in screening donors with rare blood antigens, such as weak JK alleles. The development of WGS analysis in immunogenetics laboratories would offer personalized care in the management of RBC disorders.

Massively parallel sequencing of sex-chromosomal STRs in Saudi Arabia reveals patrilineage-associated sequence variants

Massively parallel sequencing (MPS) of forensic STRs has the potential to reveal additional allele diversity compared to conventional capillary electrophoresis (CE) typing strategies, but population studies are currently relatively few in number. The Verogen ForenSeq™ DNA Signature Prep Kit includes both Y-STRs and X-STRs among its targeted loci, and here we report the sequences of these loci, analysed using Verogen's ForenSeq™ Universal Analysis Software (UAS) v1.3 and STRait Razor v3.0, in a representative sample of 89 Saudi Arabian males. We identified 56 length variants (equivalent to CE alleles) and 75 repeat sequence sub-variants across the six X-STRs analysed; equivalent figures for the set of 24 Y-STRs were 147 and 192 respectively. We also observed two flanking sequence variants for the X-, and six for the Y-STRs. Recovery of sequence data and concordance with CE data (where available) across the tested loci was good, though rare flanking variation affected interpretation and allele calling at DYF387S1 and DXS7132. Examination of flanking sequences of the Y-STRs revealed five SNPs (L255, M4790, BY7692, Z16708 and S17543) previously shown to define specific haplogroups by Y-chromosome sequencing. These define Y-haplogroups in 62 % of our sample, a proportion that increases to 91 % when haplogroup-associated repeat-sequence motifs are also considered. A population-level comparison of the Saudi Arabian X-STRs with a global sample showed our dataset to be part of a large cluster of populations of West Eurasian and Middle Eastern origin.

Past, present and future contributions of evolutionary biology to wildlife forensics, management and conservation

Successfully implementing fundamental concepts into concrete applications is challenging in any given field. It requires communication, collaboration and shared will between researchers and practitioners. We argue that evolutionary biology, through research work linked to conservation, management and forensics, had a significant impact on wildlife agencies and department practices, where new frameworks and applications have been implemented over the last decades. The Quebec government's Wildlife Department (MFFP: Ministère des Forêts, de la Faune et des Parcs) has been proactive in reducing the “research–implementation” gap, thanks to prolific collaborations with many academic researchers. Among these associations, our department's outstanding partnership with Dr. Louis Bernatchez yielded significant contributions to harvest management, stocking programmes, definition of conservation units, recovery of threatened species, management of invasive species and forensic applications. We discuss key evolutionary biology concepts and resulting concrete examples of their successful implementation that derives directly or indirectly from this successful partnership. While old and new threats to wildlife are bringing new challenges, we expect recent developments in eDNA and genomics to provide innovative solutions as long as the research–implementation bridge remains open.

Ancestry Prediction Comparisons of Different AISNPs for Five Continental Populations and Population Structure Dissection of the Xinjiang Hui Group via a Self-Developed Panel

Ancestry informative markers are genetic markers that show distinct genetic divergences among different populations. These markers can be utilized to discern population substructures and estimate the ancestral origins of unknown individuals. Previously, we developed a multiplex system of 30 ancestry informative single nucleotide polymorphism (AISNP) loci to facilitate ancestral inferences in different continental populations. In the current study, we first compared the ancestry resolutions of the 30 AISNPs and the other previously reported AISNP panels for African, European, East Asian, South Asian and American populations. Next, the genetic components of the Xinjiang Hui group were further explored in comparison to these continental populations based on the 30 AISNPs. Genetic divergence analyses of the 30 AISNPs in these five continental populations revealed that most of the AISNPs showed high genetic differentiations between these populations. Ancestry analysis comparisons of the 30 AISNPs and other published AISNPs revealed that these 30 AISNPs had comparable efficiency to other AISNP panels. Genetic relationship analyses among the studied Hui group and other continental populations demonstrated that the Hui group had close genetic affinities with East Asian populations and might share the genetic ancestries with East Asian populations. Overall, the 30 AISNPs can be used to predict the bio-geographical origins of different continental populations. Moreover, the obtained genetic data of 30 AISNPs in the Hui group can further enrich the extant reference data, which can be used as reference data for ancestry analyses of the Hui group.

PCR inhibition in qPCR, dPCR and MPS-mechanisms and solutions

DNA analysis has seen an incredible development in terms of instrumentation, assays and applications over the last years. Massively parallel sequencing (MPS) and digital PCR are now broadly applied in research and diagnostics, and quantitative PCR is used for more and more practises. All these techniques are based on in vitro DNA polymerization and fluorescence measurements. A major limitation for successful analysis is the various sample-related substances that interfere with the analysis, i.e. PCR inhibitors. PCR inhibition affects library preparation in MPS analysis and skews quantification in qPCR, and some inhibitors have been found to quench the fluorescence of the applied fluorophores. Here, we provide a deeper understanding of mechanisms of specific PCR inhibitors and how these impact specific analytical techniques. This background knowledge is necessary in order to take full advantage of modern DNA analysis techniques, specifically for analysis of samples with low amounts of template and high amounts of background material. The classical solution to handle PCR inhibition is to purify or dilute DNA extracts, which leads to DNA loss. Applying inhibitor-tolerant DNA polymerases, either single enzymes or blends, provides a more straightforward and powerful solution. This review includes mechanisms of specific PCR inhibitors as well as solutions to the inhibition problem in relation to cutting-edge DNA analysis.

Interpol review of forensic biology and forensic DNA typing 2016-2019

This review paper covers the forensic-relevant literature in biological sciences from 2016 to 2019 as a part of the 19th Interpol International Forensic Science Managers Symposium. The review papers are also available at the Interpol website at: https://www.interpol.int/content/download/14458/file/Interpol%20Review%20Papers%202019.pdf.

Utility of genetic variation in coat color genes to distinguish wild, domestic and hybrid South American camelids for forensic and judicial applications

A molecular genetic protocol for distinguishing pure and hybrid South American camelids was developed to provide strong, quantifiable, and unbiased species identification. We detail the application of the approach in the context of a criminal case in the Andes Mountains of central Chile where the defendants were alleged to have illegally hunted three wild guanacos (Lama guanicoe), as opposed to hybrid domestic llama (Lama glama)/wild guanaco crosses, which are unregulated. We describe a workflow that differentiates among wild, domestic and hybrid South American camelids (Lama versus Vicugna) based on mitochondrial cytochrome b genetic variation (to distinguish between Lama and Vicugna), and MC1R and exon 4 variation of the ASIP gene (to differentiate wild from domestic species). Additionally, we infer the population origin and sex of each of the three individuals from a panel of 15 autosomal microsatellite loci and the presence or absence of the SRY gene. Our analyses strongly supported the inference that the confiscated carcasses corresponded with 2 male and 1 female guanacos that were hunted illegally. Statistical power analyses suggested that there was an extremely low probability of misidentifying domestic camelids as wild camelids (an estimated 0 % Type I error rate), or using more conservative approached a 1.17 % chance of misidentification of wild species as domestic camelids (Type II error). Our case report and methodological and analytical protocols demonstrate the power of genetic variation in coat color genes to identify hybrids between wild and domestic camelid species and highlight the utility of the approach to help combat illegal wildlife hunting and trafficking.

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.

Detection and analysis of the cause of false-tetra-allelic patterns of locus D10S1435 at the sequence level

A number of artifacts produced in forensic DNA typing make the interpretation more complicated and even lead to typing errors. Here, we reported the cause of false-tetra-allelic patterns of STR locus D10S1435 at the sequence level. To confirm the true genotyping, the sample with four allelic peaks was re-amplified and sequenced. The amplicon sequences of D10S1435, D20S482, D6S1017, and D10S1248 loci were analyzed by software BioXM and RNAstructure. We successfully reproduced the four-peak phenomenon by adding various concentration of magnesium chloride into the loading mixtures to simulate the suboptimal electrophoresis conditions. The false four allelic peaks may be caused by the specific nucleotide sequence of locus D10S1435 which tends to form secondary structures under the suboptimal electrophoresis conditions. The relatively high GC content and extremely uneven distribution give the amplicon a potency to resist complete denaturation at the phase of sample preparation and a tendency to form intra- and intermolecular secondary structures during post-injection.