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.     

Assessing non-LUS stutter in DNA sequence data

Forensic DNA analysis is among the most well-recognized and well-developed forensic disciplines. The field’s use of DNA markers known as short tandem repeats (STRs) offer a robust means of discriminating individuals while also introducing challenges to the analysis. One of these challenges, stutter, is the result of a non-biological artifact introduced during PCR. The formation and amplification of these stutter products can occur at rates as high as 15–20% of the parent allele. The challenge inherent in this process is differentiating stutter artifacts from true alleles, particularly in the presence of a minor contributor. Traditionally, DNA profiles are obtained using capillary electrophoresis (CE), where amplified DNA fragments are separated by size, not sequence, and the identification of stutter is performed on a locus-specific level. The use of CE-based fragment data rather than sequence-based data, has limited the community’s understanding of the precise behavior of stutter. Massively parallel sequencing (MPS) data provides an opportunity to better characterize stutter, permitting a more accurate means of detecting both size- or longest uninterrupted stretch (LUS)-based stutter but also allele and motif-specific stutter characteristics. This study sheds light on the value of characterizing motif- and allele-specific stutter, including non-LUS stutter, when using MPS methods. Analysis and characterization of stutter sequences was performed using data generated from 539 samples amplified with the ForenSeq and PowerSeq 46GY library preparation kit and sequenced on the Illumina MiSeq FGx. Assessment of non-LUS stutter begins with calculating stutter rates for all potential stutter products at a given locus (and allele), additionally, the occurrence of these discrete stutter products were quantified. Results show that although the LUS sequence stutters at a higher rate than non-LUS motifs, the non-LUS stutter products do occur at detectable levels and potentially influence sequence-based mixture analysis. The data indicate that the stutter from one motif or allele can be distinguished from another motif or allele based on their unique stutter rates; however, the number of stutter products from each motif or allele may similarly make up the overall pool of stutter products. Motif- and allele-specific stutter models provide the most comprehensive analysis of sequence stutter rates and provide the ability to differentiate stutter sequences more accurately from true allele stutter. This information provides a foundation for including the characterization of non-LUS stutter products when analyzing DNA profiles, specifically mixtures with potential low-level contributors.     

DNA Commission of the International Society for Forensic Genetics: Revised and extended guidelines for mitochondrial DNA typing

The DNA Commission of the International Society of Forensic Genetics (ISFG) regularly publishes guidelines and recommendations concerning the application of DNA polymorphisms to the question of human identification. Previous recommendations published in 2000 addressed the analysis and interpretation of mitochondrial DNA (mtDNA) in forensic casework. While the foundations set forth in the earlier recommendations still apply, new approaches to the quality control, alignment and nomenclature of mitochondrial sequences, as well as the establishment of mtDNA reference population databases, have been developed. Here, we describe these developments and discuss their application to both mtDNA casework and mtDNA reference population databasing applications. While the generation of mtDNA for forensic casework has always been guided by specific standards, it is now well-established that data of the same quality are required for the mtDNA reference population data used to assess the statistical weight of the evidence. As a result, we introduce guidelines regarding sequence generation, as well as quality control measures based on the known worldwide mtDNA phylogeny, that can be applied to ensure the highest quality population data possible. For both casework and reference population databasing applications, the alignment and nomenclature of haplotypes is revised here and the phylogenetic alignment proffered as acceptable standard. In addition, the interpretation of heteroplasmy in the forensic context is updated, and the utility of alignment-free database searches for unbiased probability estimates is highlighted. Finally, we discuss statistical issues and define minimal standards for mtDNA database searches.     

Developmental validation of STRmix™, expert software for the interpretation of forensic DNA profiles

In 2015 the Scientific Working Group on DNA Analysis Methods published the SWGDAM Guidelines for the Validation of Probabilistic Genotyping Systems [1]. STRmix™ is probabilistic genotyping software that employs a continuous model of DNA profile interpretation. This paper describes the developmental validation activities of STRmix™ following the SWGDAM guidelines. It addresses the underlying scientific principles, and the performance of the models with respect to sensitivity, specificity and precision and results of interpretation of casework type samples. This work demonstrates that STRmix™ is suitable for its intended use for the interpretation of single source and mixed DNA profiles.     

A comparison of stochastic variation in mixed and unmixed casework and synthetic samples

Understanding the behaviour of mixed DNA profiles is of paramount importance in forensic DNA analysis. Key parameters are those of heterozygote balance and mixture proportion and its variability. These parameters have been previously explored as a function of the average peak height of the active alleles in single source and mixed samples derived from pristine DNA. Here we report a comparison of this data with data obtained from casework samples. This allows an assessment of the difference in the distribution of heterozygote balance between mixed and single source stains and between casework mixtures and synthetic mixtures constructed from pristine DNA.     

Presentation Of DNA Data in Cases of Disputed Paternity

Flip-open style cell phones were investigated for the potential to produce quality genetic profiles that could be used in forensic casework. Swabs were taken of the outside/back and the inside ear speaker of ten flip-phones on two occasions – prior to and seven days after cleaning with 95% ethanol. Buccal swabs were collected as exemplars. The samples were amplified using the AmpFlSTR ProfilerPlus PCR Kit for 35 cycles and STR profiles were generated using an ABI Prism 310 Genetic Analyzer and GeneMapper ID analysis software v3.2. The phone profiles were compared to the references and to each other, to assess the quality of the profiles. The completeness of the profiles varied greatly, even within an experimental condition. There was no significant difference in the percentage correct alleles or in the number of drop-in alleles in the DNA profiles for the outside/inside locations or for the pre/post-cleaning times. The findings of this study demonstrate the need for collecting multiple samples from discrete locations on a cell phone if such evidence is encountered in forensic cases.     

RNA/DNA co-analysis from blood stains--results of a second collaborative EDNAP exercise

A second collaborative exercise on RNA/DNA co-analysis for body fluid identification and STR profiling was organized by the European DNA Profiling Group (EDNAP). Six human blood stains, two blood dilution series (5-0.001 μl blood) and, optionally, bona fide or mock casework samples of human or non-human origin were analyzed by the participating laboratories using a RNA/DNA co-extraction or solely RNA extraction method. Two novel mRNA multiplexes were used for the identification of blood: a highly sensitive duplex (HBA, HBB) and a moderately sensitive pentaplex (ALAS2, CD3G, ANK1, SPTB and PBGD). The laboratories used different chemistries and instrumentation. All of the 18 participating laboratories were able to successfully isolate and detect mRNA in dried blood stains. Thirteen laboratories simultaneously extracted RNA and DNA from individual stains and were able to utilize mRNA profiling to confirm the presence of blood and to obtain autosomal STR profiles from the blood stain donors. The positive identification of blood and good quality DNA profiles were also obtained from old and compromised casework samples. The method proved to be reproducible and sensitive using different analysis strategies. The results of this collaborative exercise involving a RNA/DNA co-extraction strategy support the potential use of an mRNA based system for the identification of blood in forensic casework that is compatible with current DNA analysis methodology.     

Sequence polymorphisms of forensic Y-STRs revealed by a 68-plex in-house massively parallel sequencing panel

Sequence polymorphisms of Y chromosome short tandem repeat (Y-STR) markers can be unveiled using next generation sequencing (NGS). Compared to capillary electrophoresis, NGS has the advantage of distinguishing between some alleles of the same length. Here, a 68-plex in-house panel covering 67 Y-STR loci and the sex determinant Amelogenin locus, was developed. The accuracy of this panel was 100% concordant with three standard reference samples. The sensitive was as low as 250 pg. A total of 466 length-based alleles, 806 sequence-based alleles, and 149 haplotypes were observed across 149 Chinese Han individuals. The total haplotype diversity and discrimination capacity was 1.0000 in detected samples. The DYS710 locus possessed the highest diversity by sequence among these Y-STRs, with 109 sequence-based alleles observed. Micro-variant alleles with the same length were observed in 39 Y-STR loci, with their sequence variations mainly attributable to repeat pattern variations. While the number of sequence-based alleles identified for DYS447, DYS449, DYS710, DYS720 and DYF387S1a/b was approximately three times that of their length-based alleles, flanking sequence variations were observed in 18 alleles. In addition, 201 sequence-based alleles in 42 loci were newly discovered. This significantly expanded the knowledge of human Y-STR sequence polymorphisms. Collectively, the 68-plex panel provided reliable Y-STR results as well as higher resolution for paternal lineage analysis.     

Suspension fluorescence in situ hybridization (S-FISH) combined with automatic detection and laser microdissection for STR profiling of male cells in male/female mixtures

Laser microdissection is a valuable tool for isolating specific cells from mixtures, such as male cells in a mixture with female cells, e.g., in cases of sexual assault. These cells can be stained with Y-chromosome-specific probes. We developed an automatic screening method to detect male cells after fluorescence in situ hybridization in suspension (S-FISH). To simulate forensic casework, the method was tested on female saliva after cataglottis (a kiss involving tongue-to-tongue contact) and on licking traces (swabs of dried male saliva on female skin) even after drying. After isolation of the detected cells, short tandem repeat profiling was performed. Full DNA profiles could consistently be obtained from as little as ten buccal cells. Isolation of five cells resulted in a mean of 98% (SD of 3.4%) of the alleles detected, showing that the developed S-FISH staining had no significant negative influence on DNA recovery and can be used in forensic casework.     

Development and validation of a fast and automated DNA identification line

The importance of DNA evidence for gaining investigative leads demands a fast workflow for forensic DNA profiling performed in large volumes. Therefore, we developed software solutions for automated DNA profile analysis, contamination check, major donor inference, DNA database (DDB) comparison and reporting of the conclusions. This represents the Fast DNA IDentification Line (FIDL) and this study describes its development, validation and implementation in criminal casework at the authors’ institute. This first implementation regards single donor profiles and major contributors to mixtures. The validation included testing of the software components on their own and examination of the performance of different DDB search strategies. Furthermore, end-to-end testing was performed under three conditions: (1) testing of scenarios that can occur in DNA casework practice, (2) tests using three months of previous casework data, and (3) testing in a casework production environment in parallel to standard casework practices. The same DNA database candidates were retrieved by this automated line as by the manual workflow. The data flow was correct, results were reproducible and robust, results requiring manual analysis were correctly flagged, and reported results were as expected. Overall, we found FIDL valid for use in casework practice in our institute. The results from FIDL are automatically reported within three working days from receiving the trace sample. This includes the time needed for registration of the case, DNA extraction, quantification, polymerase chain reaction and capillary electrophoresis. FIDL itself takes less than two hours from intake of the raw CE data to reporting. Reported conclusions are one of five options: (1) candidate retrieved from DDB, (2) no candidate retrieved from DDB, (3) high evidential value with regards to reference within the case, (4) results require examination of expert, or (5) insufficient amount of DNA obtained to generate a DNA profile. In our current process, the automated report is sent within three working days and a complete report, with confirmation of the FIDL results, and signed by a reporting officer is sent at a later time. The signed report may include additional analyses regarding e.g. minor contributors. The automated report with first case results is quickly available to the police enabling them to act upon the DNA results prior to receiving the full DNA report. This line enables a uniform and efficient manner of handling large numbers of traces and cases and provides high value investigative leads in the early stages of the investigation.     

On the exact distribution of the numbers of alleles in DNA mixtures

When more than one individual contributes biological material to a forensic stain, the resulting DNA type is termed a DNA mixture. DNA mixtures occur frequently in forensic genetic casework, and in recent years, much research has been devoted to this subject. This paper presents a derivation of the exact distribution of the number of alleles for any number of profiles and investigated loci. The per locus number of observed alleles is of interest as it indicates the plausible range on the number of contributors. Hence, by specifying a prior distribution on the number of contributors, the locus distribution may be used to assess the number of contributors. Furthermore, the total number of alleles across all loci is used by some forensic geneticists to estimate the probability that an allele has failed to be detected (allelic drop-out).     

Allele frequencies of 31 autosomal short tandem repeat (auSTR) loci obtained using the Precision ID GlobalFiler™ NGS STR Panel v2 in 322 individuals from the Japanese population

In human identification methods that target short tandem repeats (STRs), massively parallel sequencing (MPS) technology has made it possible to genotype at the level of the specific sequence itself. This allows for the detection of repeat unit variants and single nucleotide polymorphisms (SNPs) adjacent to the STRs. Using the GlobalFiler™ NGS STR Panel v2, Ion S5, and Converge software, this study constructed a Japanese database of 31 autosomal STRs (auSTRs) and two sex markers from 322 individuals. After excluding some sequence errors and stutters, a total of 31 novel alleles were identified. Additionally, using the allele frequencies of 31 auSTR loci, the match probabilities for the length-based and sequence-based data were calculated to be 1.433 × 10⁻³⁴ and 9.163 × 10⁻³⁸, respectively. These values are at least nine orders of magnitude higher than that obtained from 21 auSTR loci in the Japanese population using the conventional capillary electrophoresis method. The database generated in this study is expected to be implemented in forensic practice and used to solve difficult casework.     

New application for haplotype-specific extraction: Separation of mitochondrial DNA mixtures

A particularly challenging aspect of forensic casework is the accurate interpretation of samples that contain DNA from more than one contributor, since DNA analysis of specific identification markers often results in ambiguous mixed profiles. Physically separating the various haploid components of a diploid or polyploid DNA sample or of a specimen containing DNA from multiple individuals can resolve this problem. We adapted and combined previously developed methods of haplotype-specific extraction (HSE) to separate and analyze mitochondrial DNA (mtDNA) in samples containing DNA from two individuals. After optimizing conditions and probes, we performed HSE to separate mitochondrial DNA mixtures, which manifest as mixed nucleotide calls at single base positions. We were able to subsequently identify the DNA of the two individuals by sequencing. Based on our findings from this proof-of-concept experiment, this novel assay will be useful for distinguishing among the mtDNA of individuals in mixed DNA samples.     

Increasing amplification success of forensic DNA samples using multiple displacement amplification

Multiple displacement amplification (MDA) is capable of amplifying nanogram template amounts of DNA with high accuracy to generate micrograms of representative product. Although MDA is able to amplify small template amounts (<100 pg), increased levels of preferential amplification and allelic dropout are observed. The use of molecular crowders (polyethylene glycol 400) can decrease the amplification bias and increase amplification success. We describe the use of standard and crowded MDA on low-concentration casework samples originating from blood, semen, saliva, hair, and trace DNA. While standard MDA produced no significant increases in STR genotyping success, the use of crowded MDA yielded an average increase of 15% in the number of alleles, with a significant decrease in the amplification bias, resulting in clearer, more-complete profiles from forensically relevant samples. Presented at the 18th International Symposium on the Forensic Sciences organised by the Australian and New Zealand Forensic Science Society (ANZFS), Fremantle, Western Australia, Australia 2–7 April 2006.     

FDSTools: A software package for analysis of massively parallel sequencing data with the ability to recognise and correct STR stutter and other PCR or sequencing noise

Massively parallel sequencing (MPS) is on the advent of a broad scale application in forensic research and casework. The improved capabilities to analyse evidentiary traces representing unbalanced mixtures is often mentioned as one of the major advantages of this technique. However, most of the available software packages that analyse forensic short tandem repeat (STR) sequencing data are not well suited for high throughput analysis of such mixed traces. The largest challenge is the presence of stutter artefacts in STR amplifications, which are not readily discerned from minor contributions. FDSTools is an open-source software solution developed for this purpose. The level of stutter formation is influenced by various aspects of the sequence, such as the length of the longest uninterrupted stretch occurring in an STR. When MPS is used, STRs are evaluated as sequence variants that each have particular stutter characteristics which can be precisely determined. FDSTools uses a database of reference samples to determine stutter and other systemic PCR or sequencing artefacts for each individual allele. In addition, stutter models are created for each repeating element in order to predict stutter artefacts for alleles that are not included in the reference set. This information is subsequently used to recognise and compensate for the noise in a sequence profile. The result is a better representation of the true composition of a sample. Using Promega Powerseq™ Auto System data from 450 reference samples and 31 two-person mixtures, we show that the FDSTools correction module decreases stutter ratios above 20% to below 3%. Consequently, much lower levels of contributions in the mixed traces are detected. FDSTools contains modules to visualise the data in an interactive format allowing users to filter data with their own preferred thresholds.     

Presence and potential of cell free DNA in different types of forensic samples

Extracellular or cell free DNA has been found to exist in many biological media such as blood and saliva. To check whether cell free DNA is present in the supernatant which is normally discarded during several DNA extraction processes, such as Chelex^® extraction, DNA profiles of cell pellet and concentrated supernatant from 30 artificial case like samples and from 100 real forensic samples were compared. Presence of cell free DNA was shown in all investigated sample types. Moreover, in some samples additional alleles, not detected during analysis of the cell pellet, were detected, offering valuable information which would normally have been discarded together with the supernatant. The results presented here indicate that cell free DNA deserves further consideration since it has the potential to increase the DNA yield in forensic casework samples in general and in contact traces in particular.     

Massively parallel sequencing for the forensic scientist – sequencing archived amplified products of AmpFlSTR Identifiler and PowerPlex Y multiplex kits to capture additional information

In the last few years the cost and ease of massively parallel sequencing (MPS) has reduced dramatically to the point that it can now be considered as a tool for use in forensic case work. An important consideration for the implementation of any new forensic technology is the ability to remain compatible with previous technology. During this study we sequenced the amplicons of two commercial forensic short tandem repeat (STR) multiplexes AmpFlSTR Identifiler and PowerPlex Y using the Illumina MiSeq and Ion PGM Sequencer (Life Technologies) and characterised the sequence data from a forensic perspective. Using the MPS data from both platforms we determined the STR genotypes of forensic samples and found previously undocumented sequence variation in seven STR alleles. By characterising features of the DNA sequence profiles, such as stutter and locus imbalance we identified areas for future development that will be needed prior to casework implementation. The rapid development of this technology has meant many in the forensic community have been ‘left behind’. We also provide an explanation, for forensic scientists, of what is happening at the different stages of the MPS workflow, from library preparation through to bioinformatics, and how this may affect the results.     

Massively parallel sequence data of 31 autosomal STR loci from 496 Spanish individuals revealed concordance with CE-STR technology and enhanced discrimination power

This study reports Short Tandem Repeat (STR) sequence-based allele data from 496 Spanish individuals across 31 autosomal STR (auSTR) loci included in the Precision ID GlobalFiler™ NGS STR Panel v2: D12S391, D13S317, D8S1179, D21S11, D3S1358, D5S818, D1S1656, D2S1338, vWA, D2S441, D5S2800, D7S820, D16S539, D6S474, D12ATA63, D4S2408, D6S1043, D19S433, D14S1434, CSF1PO, D10S1248, D18S51, D1S1677, D22S1045, D2S1776, D3S4529, FGA, Penta D, Penta E, TH01 and TPOX. The sequence of each allele was aligned to the reference sequence GRCh37 (hg19) and formatted according to the guidance of the International Society for Forensic Genetics. A subset of 221 samples was evaluated for testing concordance with allele calls derived from CE-based analysis using PowerPlex Fusion 6C, and there was 99.95% allele concordance. Twenty-five out of 31 auSTR loci showed an increased number of alleles due to repeat region sequence variation and/or single nucleotide polymorphisms (SNP) residing in the flanking regions. A total of 18 loci showed increased observed heterozygosity due to sequence variation; the loci exhibiting the greatest increase were: D13S317 (12% points), D5S818 (10% points), D8S1179 (7% points), D3S1358 (7% points), and D21S11 (6% points). The combined match probability decreased from 2.022E-24 (length-based data) to 1.042E-27 (sequence-based data) for the 20 CODIS core STR loci. The combined match probability (sequence-based data) for the 31 STR loci studied was 4.777E-40. The combined typical paternity index increased from 1.118E + 12 to 8.179E + 13 using length and sequence-based data, respectively. This Spanish population study performed in the framework of the EU-funded DNASEQEX project is expected to provide STR sequence-based allele frequencies for forensic casework and support implementation of massively parallel sequencing (MPS) technology in forensic laboratories.     

Short tandem repeat typing on the 454 platform: Strategies and considerations for targeted sequencing of common forensic markers

To investigate the feasibility of next generation sequencing technology (NGS) for the multiplex detection and sequence production of short tandem repeats (STRs) from degraded and low DNA quantity samples, standard polymerase chain reaction amplification methods were used to enrich for commonly employed STR markers. Samples were amplified with two multiplexing strategies: a multiplex containing thirteen miniSTR markers and a series of multiplexes containing four miniSTR markers each. Each sample multiplex was barcoded with a sample-specific multiplex identifier for subsequent parallel tagged sequencing on the GS Junior System (454 Life Sciences, a Roche company, Branford, CT). Sequencing results from over fifty DNA extracts representing both pristine samples and low-quality evidentiary specimens reflected known genotypes and were consistent across multiple extracts and/or amplifications of the same sample. Furthermore, the NGS data revealed sequence information not available with standard capillary electrophoresis-based detection alone. For the population samples tested, a total of 152 single nucleotide polymorphisms or insertions/deletions were identified in over 935 recovered alleles, averaging one polymorphism for every six recovered alleles. For three of the loci, the sequence information doubled the number of alleles detected via traditional STR typing by fragment analysis. In addition, twenty-eight of these variants were only seen once within our dataset, highlighting the potential for discrimination among individuals. These additional data are likely to be particularly valuable in missing persons and disaster victim identification cases for which only partial profiles may be recovered and/or only distant kin are available as references. And, considering the opportunity to target only small amplicons with NGS, this type of STR typing will allow for greater information recovery from challenging casework samples. While our results highlight the potential of new technologies for recovering discriminatory genetic information from evidentiary specimens, our data also reveal the complexities of NGS-based STR typing, both in terms of the laboratory assays themselves as well as the downstream data processing and analysis.     

D1S80 subtyping by PCR-RFLP: new nomenclature and further characterization

In a previous study, we screened out sequence variations within alleles at the D1S80 locus of a Japanese population using PCR-RFLP with EcoRII as a restriction enzyme. In the present study, through analyzing the alleles in a Chinese population, we were able to demonstrate four new electrophoretic band patterns that were complementary to the Japanese data. After summarizing the band patterns and sequencing results of these two populations, we established a new nomenclature for the PCR-RFLP band patterns, closely relating them to their corresponding sequences so that the new types could be designated easily and accurately. After PCR-RFLP subtyping, nineteen alleles in the Chinese population were revealed to have a total of thirty-seven subtypes. The discrimination power of this locus in the Chinese population was elevated from 0.974 to 0.988, and the Hardy-Weinberg equilibrium in this population showed no deviation when checked. Two samples typed as homozygotes 24/24 and 30/30 were identified to be actually heterozygous according to their band patterns. The result was supported by the sequencing analysis of the two samples in which overlapping of eight and eleven repeat units, respectively, were revealed. The heterozygosity was thus elevated from 0.85 to 0.87. The present study proved that PCR-RFLP was an effective method for subtyping D1S80 alleles in the Chinese.