A comparison of adjustment methods to test the robustness of an STR DNA database comprised of 24 European populations

DNA databases of an important tropical timber tree species Shorea leprosula (Dipterocarpaceae) for forensic timber identification

International timber trade communities are increasingly demanding that timber in the wood supply chain be sourced from sustainably harvested forests and certified plantations. This is to combat illegal logging activities to prevent further depletion of our precious forests worldwide. Hence, timber tracking tools are important to support law enforcement officials in ensuring only sustainably harvested timbers are traded in the market. In this study, we developed chloroplast DNA (cpDNA) and simple sequence repeat (SSR) databases as tracking tools for an important tropical timber tree species, Shorealeprosula from Peninsular Malaysia. A total of 1410 individual trees were sampled from 44 natural populations throughout Peninsular Malaysia. Four cpDNA regions were used to generate a cpDNA haplotype database, resulting in a haplotype map comprising 22 unique haplotypes derived from 28 informative intraspecific variable sites. This cpDNA database can be used to trace the origin of an unknown log at the regional level. Ten SSR loci were used to develop the SSR allele frequency database. Bayesian cluster analysis divided the 44 populations into two genetic clusters corresponding to Region A and Region B. Based on conservativeness evaluation of the SSR databases for individual identification, the coancestry coefficients (θ) were adjusted to 0.1900 and 0.1500 for Region A and B, respectively. These databases are useful tools to complement existing timber tracking systems in ensuring only legally sourced timbers are allowed to enter the wood supply chain.

DNA databases of a CITES listed species Aquilaria malaccensis (Thymelaeaceae) as the tracking tools for forensic identification and chain of custody certification

Aquilaria malaccensis (Thymelaeaceae) is the main source of high-grade agarwood in Southeast Asia. Aggressive collections and trade activities over the past decades have put great pressure on the natural stands and raised concerns over the long-term survival potential of A. malaccensis. Tracking and authentication of agarwood require method with a high degree of accuracy. Therefore, this study aimed to develop DNA databases of A. malaccensis as the tracking tools at species, population and individual levels for forensic identification and chain of custody certification. Using two cpDNA (rbcL and matK) and an rDNA (ITS2) markers, species identification database of Aquilaria was developed to distinguish A. malaccensis from A. hirta, A. microcarpa, A. beccariana, A. crassna, A. sinensis and A. rostrata. In addition, based on 35 populations of A. malaccensis throughout Peninsular Malaysia, cpDNA haplotype and STR allele frequency databases were developed for population and individual identification. A haplotype distribution map based on 29 haplotypes derived from seven cpDNA showed that the A. malaccensis in Peninsular Malaysia can be associated to Kedah-Perak and Kelantan-Johor regions. Similarly, genetic relatedness and Bayesian clustering analyses based on 10 STR markers also divided the 35 populations into two main genetic clusters, corresponding to Kedah-Perak and Kelantan-Johor regions. The STR allele frequency databases were established and characterized according to these two regions. To determine the performance of the STR allele frequency databases for population identification, independent self-assignment tests showed that the percentage of individuals correctly assigned into the origin population was 93.88% in Kedah-Perak and 90.29% in Kelantan-Johor. For the STR allele frequency databases to be used for individual identification, conservativeness tests showed that the θ should be adjusted to 0.250 and 0.200 in the Kedah-Perak and Kelantan-Johor databases, respectively. To ensure consistency in allele calling for the dinucleotide repeat loci across different electrophoretic platforms or laboratories, allelic ladders have been developed for the 10 STR loci. Two case studies are presented of how these databases were used to track A. malaccensis to the origin population and stump. These databases are ready to be used to provide admissible forensic evidence for legal proceedings against the illegal harvesters of agarwood and for agarwood certification to meet the consumer country regulations.

Sequenced-based paternity analysis to improve breeding and identify self-incompatibility loci in intermediate wheatgrass (Thinopyrum intermedium)

KEY MESSAGE

Paternity assignment and genome-wide association analyses for fertility were applied to a Thinopyrum intermedium breeding program. A lack of progeny between combinations of parents was associated with loci near self-incompatibility genes. In outcrossing species such as intermediate wheatgrass (IWG, Thinopyrum intermedium), polycrossing is often used to generate novel recombinants through each cycle of selection, but it cannot track pollen-parent pedigrees and it is unknown how self-incompatibility (SI) genes may limit the number of unique crosses obtained. This study investigated the potential of using next-generation sequencing to assign paternity and identify putative SI loci in IWG. Using a reference population of 380 individuals made from controlled crosses of 64 parents, paternity was assigned with 92% agreement using Cervus software. Using this approach, 80% of 4158 progeny (n = 3342) from a polycross of 89 parents were assigned paternity. Of the 89 pollen parents, 82 (92%) were represented with 1633 unique full-sib families representing 42% of all potential crosses. The number of progeny per successful pollen parent ranged from 1 to 123, with number of inflorescences per pollen parent significantly correlated to the number of progeny (r = 0.54, p < 0.001). Shannon's diversity index, assessing the total number and representation of families, was 7.33 compared to a theoretical maximum of 8.98. To test our hypothesis on the impact of SI genes, a genome-wide association study of the number of progeny observed from the 89 parents identified genetic effects related to non-random mating, including marker loci located near putative SI genes. Paternity testing of polycross progeny can impact future breeding gains by being incorporated in breeding programs to optimize polycross methodology, maintain genetic diversity, and reveal genetic architecture of mating patterns.

The STRidER Report on Two Years of Quality Control of Autosomal STR Population Datasets

STRidER, the STRs for Identity ENFSI Reference Database, is a curated, freely publicly available online allele frequency database, quality control (QC) and software platform for autosomal Short Tandem Repeats (STRs) developed under the endorsement of the International Society for Forensic Genetics. Continuous updates comprise additional STR loci and populations in the frequency database and many further STR-related aspects. One significant innovation is the autosomal STR data QC provided prior to publication of datasets. Such scrutiny was lacking previously, leaving QC to authors, reviewers and editors, which led to an unacceptably high error rate in scientific papers. The results from scrutinizing 184 STR datasets containing >177,000 individual genotypes submitted in the first two years of STRidER QC since 2017 revealed that about two-thirds of the STR datasets were either being withdrawn by the authors after initial feedback or rejected based on a conservative error rate. Almost no error-free submissions were received, which clearly shows that centralized QC and data curation are essential to maintain the high-quality standard required in forensic genetics. While many errors had minor impact on the resulting allele frequencies, multiple error categories were commonly found within single datasets. Several datasets contained serious flaws. We discuss the factors that caused the errors to draw the attention to redundant pitfalls and thus contribute to better quality of autosomal STR datasets and allele frequency reports.

A geographical traceability system for Merbau (Intsia palembanica Miq.), an important timber species from peninsular Malaysia

To inform product users about the origin of timber, the implementation of a traceability system is necessary for the forestry industry. In this study, we developed a comprehensive genetic database for the important tropical timber species Merbau, Intsia palembanica, to trace its geographic origin within peninsular Malaysia. A total of 1373 individual trees representing 39 geographically distinct populations of I. palembanica were sampled throughout peninsular Malaysia. We analyzed the samples using a combination of four chloroplast DNA (cpDNA) markers and 14 short tandem repeat (STR) markers to establish both cpDNA haplotype and STR allele frequency databases. A haplotype map was generated through cpDNA sequencing for population identification, resulting in six unique haplotypes based on 10 informative intraspecifically variable sites. Subsequently, an STR allele frequency database was developed from 14 STRs allowing individual identification. Bayesian cluster analysis divided the individuals into two genetic clusters corresponding to the northern and southern regions of peninsular Malaysia. Tests of conservativeness showed that the databases were conservative after the adjustment of the θ values to 0.2000 and 0.2900 for the northern (f = 0.0163) and southern (f = 0.0285) regions, respectively. Using self-assignment tests, we observed that individuals were correctly assigned to populations at rates of 40.54-94.12% and to the identified regions at rates of 79.80-80.62%. Both the cpDNA and STR markers appear to be useful for tracking Merbau timber originating from peninsular Malaysia. The use of these forensic tools in addition to the existing paper-based timber tracking system will help to verify the legality of the origin of I. palembanica and to combat illegal logging issues associated with the species.

DNA commission of the International society for forensic genetics: Assessing the value of forensic biological evidence - Guidelines highlighting the importance of propositions: Part I: evaluation of DNA profiling comparisons given (sub-) source propositions

The interpretation of evidence continues to be one of the biggest challenges facing the forensic community. This is the first of two papers intended to provide advice on difficult aspects of evaluation and in particular on the formulation of propositions. The scientist has a dual role: investigator (crime-focused), where often there is no suspect available and a database search may be required; evaluator (suspect-focused), where the strength of evidence is assessed in the context of the case. In investigative mode, generally the aim is to produce leads regarding the source of the DNA. Sub-source level propositions will be adequate to help identify potential suspects who can be further investigated by the authorities. Once in evaluative mode, given the defence version of events of the person of interest, it may become necessary to consider alternatives that go beyond the source of the DNA (i.e., to consider activity level propositions). In the evaluation phase, it is crucial that formulation of propositions allows the assessment of all the results that will help with the issue at hand. Propositions should therefore be precise (indication of the number of contributors, information on the relevant population etc.), be about causes, not effects (e.g. a 'matching' DNA profile) and to avoid bias, must not be findings-led. This means that ideally, propositions should be decided based on the case information and before the results of the comparisons are known. This paper primarily reflects upon what has been coined as "sub-source level propositions". These are restricted to the evaluation of the DNA profiles themselves, and help answer the issue regarding the source of the DNA. It is to be emphasised that likelihood ratios given sub-source level propositions cannot be carried over to a different level - for example, activity level propositions, where the DNA evidence is put into the context of the alleged activities. This would be highly misleading and could give rise to miscarriages of justice; this will be discussed in a second paper. The value of forensic results depends not only on propositions, but also on the type of results (e.g. allelic designations, peak heights, replicates) and upon the model used: it is therefore important to discuss these aspects. Finally, since communication is key to help understanding by courts, we will explore how to convey the value of the results and explain the importance of avoiding the practice of transposing the conditional.

Evaluating DNA evidence in a genetically complex population

In forensic genetics, the likelihood ratio (LR), measuring the value of DNA profile evidence, is computed from a database of allele frequencies. Here, we address the choice of database and adjustments for population structure and sample size in the context of Brazil. The Brazilian population underwent a complex process of colonization, migration and mating, which created an admixed genetic composition that makes it difficult to obtain an appropriate database for a given case. National databases are now available, as well as databases for many Brazilian states. However, those databases are not statistically random samples, and state boundaries may not accurately reflect the sub-structuring of genetic diversity. We compared the LR calculated using the relevant state-specific database with the statistics calculated when a national database and when international databases were used. We evaluated two methods of adjustment for population structure, due to Wright [13] and Balding and Nichols [14]. We also considered two adjustments for database sample size: the Balding size bias correction [15] and a minimum allele frequency [16]. Our results show that the use of a national database with the Balding and Nichols adjustment and θ = 0.002 generated lower LR values than did the state-specific database in more than 50% of the profiles simulated using the state-based allele frequencies, while θ = 0.01 produced lower LRs for more than 90% of the profiles. We conclude that the utilization of a national database for Brazilian cases can be justified in association with the appropriate adjustment for population structure.

Geographic origin and individual assignment of Shorea platyclados (Dipterocarpaceae) for forensic identification

The development of timber tracking methods based on genetic markers can provide scientific evidence to verify the origin of timber products and fulfill the growing requirement for sustainable forestry practices. In this study, the origin of an important Dark Red Meranti wood, Shorea platyclados, was studied by using the combination of seven chloroplast DNA and 15 short tandem repeats (STRs) markers. A total of 27 natural populations of S. platyclados were sampled throughout Malaysia to establish population level and individual level identification databases. A haplotype map was generated from chloroplast DNA sequencing for population identification, resulting in 29 multilocus haplotypes, based on 39 informative intraspecific variable sites. Subsequently, a DNA profiling database was developed from 15 STRs allowing for individual identification in Malaysia. Cluster analysis divided the 27 populations into two genetic clusters, corresponding to the region of Eastern and Western Malaysia. The conservativeness tests showed that the Malaysia database is conservative after removal of bias from population subdivision and sampling effects. Independent self-assignment tests correctly assigned individuals to the database in an overall 60.60−94.95% of cases for identified populations, and in 98.99−99.23% of cases for identified regions. Both the chloroplast DNA database and the STRs appear to be useful for tracking timber originating in Malaysia. Hence, this DNA-based method could serve as an effective addition tool to the existing forensic timber identification system for ensuring the sustainably management of this species into the future.

The factor of 10 in forensic DNA match probabilities

An update was performed of the classic experiments that led to the view that profile probability assignments are usually within a factor of 10 of each other. The data used in this study consist of 15 Identifiler loci collected from a wide range of forensic populations. Following Budowle et al. [1], the terms cognate and non-cognate are used. The cognate database is the database from which the profiles are simulated. The profile probability assignment was usually larger in the cognate database. In 44%-65% of the cases, the profile probability for 15 loci in the non-cognate database was within a factor of 10 of the profile probability in the cognate database. This proportion was between 60% and 80% when the FBI and NIST data were used as the non-cognate databases. A second experiment compared the match probability assignment using a generalised database and recommendation 4.2 from NRC II (the 4.2 assignment) with a proxy for the matching proportion developed using subpopulation allele frequencies and the product rule. The findings support that the 4.2 assignment has a large conservative bias. These results are in agreement with previous research results.

Forensic timber identification: a case study of a CITES listed species, Gonystylus bancanus (Thymelaeaceae)

Illegal logging and smuggling of Gonystylus bancanus (Thymelaeaceae) poses a serious threat to this fragile valuable peat swamp timber species. Using G. bancanus as a case study, DNA markers were used to develop identification databases at the species, population and individual level. The species level database for Gonystylus comprised of an rDNA (ITS2) and two cpDNA (trnH-psbA and trnL) markers based on a 20 Gonystylus species database. When concatenated, taxonomic species recognition was achieved with a resolution of 90% (18 out of the 20 species). In addition, based on 17 natural populations of G. bancanus throughout West (Peninsular Malaysia) and East (Sabah and Sarawak) Malaysia, population and individual identification databases were developed using cpDNA and STR markers respectively. A haplotype distribution map for Malaysia was generated using six cpDNA markers, resulting in 12 unique multilocus haplotypes, from 24 informative intraspecific variable sites. These unique haplotypes suggest a clear genetic structuring of West and East regions. A simulation procedure based on the composition of the samples was used to test whether a suspected sample conformed to a given regional origin. Overall, the observed type I and II errors of the databases showed good concordance with the predicted 5% threshold which indicates that the databases were useful in revealing provenance and establishing conformity of samples from West and East Malaysia. Sixteen STRs were used to develop the DNA profiling databases for individual identification. Bayesian clustering analyses divided the 17 populations into two main genetic clusters, corresponding to the regions of West and East Malaysia. Population substructuring (K=2) was observed within each region. After removal of bias resulting from sampling effects and population subdivision, conservativeness tests showed that the West and East Malaysia databases were conservative. This suggests that both databases can be used independently for random match probability estimation within respective regions. The reliability of the databases was further determined by independent self-assignment tests based on the likelihood of each individual's multilocus genotype occurring in each identified population, genetic cluster and region with an average percentage of correctly assigned individuals of 54.80%, 99.60% and 100% respectively. Thus, after appropriate validation, the genetic identification databases developed for G. bancanus in this study could support forensic applications and help safeguard this valuable species into the future.

Worldwide F(ST) estimates relative to five continental-scale populations

We estimate the population genetics parameter (also referred to as the fixation index) from short tandem repeat (STR) allele frequencies, comparing many worldwide human subpopulations at approximately the national level with continental-scale populations. is commonly used to measure population differentiation, and is important in forensic DNA analysis to account for remote shared ancestry between a suspect and an alternative source of the DNA. We estimate comparing subpopulations with a hypothetical ancestral population, which is the approach most widely used in population genetics, and also compare a subpopulation with a sampled reference population, which is more appropriate for forensic applications. Both estimation methods are likelihood-based, in which is related to the variance of the multinomial-Dirichlet distribution for allele counts. Overall, we find low values, with posterior 97.5 percentiles when comparing a subpopulation with the most appropriate population, and even for inter-population comparisons we find . These are much smaller than single nucleotide polymorphism-based inter-continental estimates, and are also about half the magnitude of STR-based estimates from population genetics surveys that focus on distinct ethnic groups rather than a general population. Our findings support the use of up to 3% in forensic calculations, which corresponds to some current practice.

Genotyping and interpretation of STR-DNA: Low-template, mixtures and database matches-Twenty years of research and development

The introduction of Short Tandem Repeat (STR) DNA was a revolution within a revolution that transformed forensic DNA profiling into a tool that could be used, for the first time, to create National DNA databases. This transformation would not have been possible without the concurrent development of fluorescent automated sequencers, combined with the ability to multiplex several loci together. Use of the polymerase chain reaction (PCR) increased the sensitivity of the method to enable the analysis of a handful of cells. The first multiplexes were simple: 'the quad', introduced by the defunct UK Forensic Science Service (FSS) in 1994, rapidly followed by a more discriminating 'six-plex' (Second Generation Multiplex) in 1995 that was used to create the world's first national DNA database. The success of the database rapidly outgrew the functionality of the original system - by the year 2000 a new multiplex of ten-loci was introduced to reduce the chance of adventitious matches. The technology was adopted world-wide, albeit with different loci. The political requirement to introduce pan-European databases encouraged standardisation - the development of European Standard Set (ESS) of markers comprising twelve-loci is the latest iteration. Although development has been impressive, the methods used to interpret evidence have lagged behind. For example, the theory to interpret complex DNA profiles (low-level mixtures), had been developed fifteen years ago, but only in the past year or so, are the concepts starting to be widely adopted. A plethora of different models (some commercial and others non-commercial) have appeared. This has led to a confusing 'debate' about the 'best' to use. The different models available are described along with their advantages and disadvantages. A section discusses the development of national DNA databases, along with details of an associated controversy to estimate the strength of evidence of matches. Current methodology is limited to searches of complete profiles - another example where the interpretation of matches has not kept pace with development of theory. STRs have also transformed the area of Disaster Victim Identification (DVI) which frequently requires kinship analysis. However, genotyping efficiency is complicated by complex, degraded DNA profiles. Finally, there is now a detailed understanding of the causes of stochastic effects that cause DNA profiles to exhibit the phenomena of drop-out and drop-in, along with artefacts such as stutters. The phenomena discussed include: heterozygote balance; stutter; degradation; the effect of decreasing quantities of DNA; the dilution effect.

Confidence interval of the likelihood ratio associated with mixed stain DNA evidence

Likelihood ratios are necessary to properly interpret mixed stain DNA evidence. They can flexibly consider alternate hypotheses and can account for population substructure. The likelihood ratio should be seen as an estimate and not a fixed value, because the calculations are functions of allelic frequency estimates that were estimated from a small portion of the population. Current methods do not account for uncertainty in the likelihood ratio estimates and are therefore an incomplete picture of the strength of the evidence. We propose the use of a confidence interval to report the consequent variation of likelihood ratios. The confidence interval is calculated using the standard forensic likelihood ratio formulae and a variance estimate derived using the Taylor expansion. The formula is explained, and a computer program has been made available. Numeric work shows that the evidential strength of DNA profiles decreases as the variation among populations increases.

A simple method for establishing concordance between short-tandem-repeat allele frequency databases

BACKGROUND

Current paternity and forensic accreditation standards do not require concordance to be established between short-tandem-repeat allele frequency databases representing the same population. The current statistical methods for evaluating databases do not establish concordance. Although acceptable under current forensic statistical methods, databases representing the same population may have sufficient variation to influence the outcome of a nondirect relationship testing result (i.e., siblingship). Hence there is a need for a quantitative method to determine concordance between databases.

STUDY DESIGN AND METHODS

Local allele frequency databases were generated for major US ethnic groups. Statistical analysis was performed as recommended by international forensic standards. A new method was developed and used for evaluating concordance between the locally developed and published databases. Smaller deviation values signify greater concordance between compared databases. The results were quantitatively confirmed against data obtained from a multidimensional scaling analysis system (SPSS, SPSS, Inc.).

RESULTS

The locally developed database deviated from the mean of the published data by approximately 0.1073 for Caucasians, 0.1341 for US Hispanics, and 0.1287 for African Americans. Upon pairwise comparison, the published databases deviated from one another by 0.1210 in US Caucasians, 0.1457 in US Hispanics, and 0.1228 in African Americans. The local database was observed to be more concordant to published databases than some of the published data to themselves.

CONCLUSION

The developed database was shown to be concordant with four previously published databases using the new method. The results were confirmed by comparison to qualitative data obtained from the multidimensional scaling analysis system.

Determination of the variables affecting mixed MiniFiler™ DNA profiles

The interpretation of mixed DNA profiles presents additional challenges for the forensic scientist. There has been a broad based call for transparency in the process of interpretation of all evidence including mixed DNA profiles. This interpretation is greatly facilitated by a sound understanding of the variability in peak heights for the two peaks of a heterozygote, in the sizes of stutter peaks and in the variability in peak heights across loci. This study examines single source and mixed DNA profiles to assess this variability. The relative variability in peak height between the two peaks of a heterozygote and in the peak heights across loci becomes greater as the peaks themselves become smaller. This is consistent with findings from other multiplexes. This variability appears larger in the MiniFiler™ system at 30 cycles than, for example, in the Identifiler™ system at 28 cycles and this difference is largely explained by the two extra cycles of amplification. Stutter peaks appear no larger in the MiniFiler™ system at 30 cycles than in the Identifiler™ system at 28 cycles.

A universal strategy to interpret DNA profiles that does not require a definition of low-copy-number

In this paper we critically examine the causes of the underlying confusion that relates to the issue of low-template (LT) DNA profile interpretation. Firstly, there is much difficulty in attempting to distinguish between LT-DNA vs. conventional DNA because there is no discrete 'cut-off' point that can be reasonably defined or evaluated. LT-DNA is loosely characterised by drop-out (where alleles may be missing) and drop-in (where additional alleles may be present). We have previously described probabilistic methods that can be used to incorporate these phenomena using likelihood ratio (LR) principles. This is preferred to the random man not excluded (RMNE) method, because we cannot identify a coherent way forward within the restrictions provided by this framework. Most LT-DNA profiles are interpreted using a 'consensus' profile method, we called this the 'biological model', where only those alleles that are duplicated in consecutive tests are reported. We recognise that there is an increased need for probabilistic models to take precedence over the biological model. These models are required for all kinds of DNA profiles, not just those that are believed to be low-template. We also recognise that there is a need for education and training if the methods we recommend are to be widely introduced.

Object-oriented Bayesian networks for paternity cases with allelic dependencies

This study extends the current use of Bayesian networks by incorporating the effects of allelic dependencies in paternity calculations. The use of object-oriented networks greatly simplify the process of building and interpreting forensic identification models, allowing researchers to solve new, more complex problems. We explore two paternity examples: the most common scenario where DNA evidence is available from the alleged father, the mother and the child; a more complex casewhere DNA is not available from the alleged father, but is available from the alleged father's brother. Object-oriented networks are built, using HUGIN, for each example which incorporate the effects of allelic dependence caused by evolutionary relatedness.

Populationsgenetik autosomaler Polymorphismen

Für die Sammlung valider populationsgenetischer Daten zu den in der forensischen Genetik gebräuchlichen autosomalen DNA-Polymorphismen müssen eine Reihe von Kriterien erfüllt werden. Diese umfassen zunächst die untersuchten polymorphen Marker und ihre Typisierungsverfahren, da nur auf der Basis einer gesicherten Unterscheidung der Allele sowie einer verbindlichen Nomenklatur die Vergleichbarkeit und Reproduzierbarkeit der gesammelten Daten gegeben ist. Die genetische Validierung umfasst den Nachweis des Mendel-Erbgangs sowie die Überprüfung der Mutationsrate. Die Selektion der Probanden für die Bestimmung der Allelhäufigkeiten, die für eine Bevölkerungsgruppe repräsentativ sein sollen, muss zufällig erfolgen, und die Daten müssen auf das Vorliegen des Hardy-Weinberg-Gleichgewichts geprüft werden. Eine allgemein anerkannte Definition von Populationen im forensisch-genetischen Kontext ist Gegenstand der aktuellen Diskussion. Daher beruht die Auswahl der Probanden aus pragmatischen Gründen primär auf ihrer geografischen Herkunft. Zusätzlich wird empfohlen, bei der biostatistischen Bewertung populationsgenetische Maßzahlen wie Fst und Θ als Korrekturen für mögliche Inhomogenitäten in der untersuchten Populationsstichprobe zu verwenden.

Scientific standards for studies in forensic genetics

Forensic molecular genetics has evolved from a rapidly developing field with changing technologies into a highly recognized and generally accepted forensic science, leading to the establishment of national DNA databases with DNA profiles from suspects and convicted offenders. DNA evidence has taken a central role by carrying a significant weight for convictions, as well as by excluding innocent suspects early on in a criminal investigation. Due to this impact on the criminal justice system, guidelines for research in forensic genetics have been introduced already since many years. The most important issues regarding the selection and definition of typing systems both for paternity testing and for forensic identification, the criteria for technical and biostatistical validation, as well as the use of mitochondrial DNA analysis are summarized and discussed.

Allele frequencies for 70 autosomal SNP loci with U.S. Caucasian, African-American, and Hispanic samples

189 samples from 3 different U.S. sample groups Caucasian (74), African American (71) and Hispanic (44) were typed for 70 autosomal genetic markers. These 70 markers are bi-allelic (C/T) short nucleotide polymorphisms (SNPs). For each sample, the 70 SNP markers were typed in 11 unique 6-plexes and a single 4-plex PCR. A total of 10 of the 210 tests (70 loci x 3 populations) for Hardy-Weinberg equilibrium indicated a statistically significant result. In order to evaluate the minimum number of SNP loci needed to distinguish all 189 samples from one another, we ranked the loci according to their levels of observed heterozygosity and p-values obtained upon testing for Hardy-Weinberg equilibrium. The top 12 loci according to these ranking criteria were tabulated along with the number of unique genotypes observed when combining subsequent SNP markers. The 12 selected SNPs possessed an observed heterozygosity of >0.45 in all three populations examined and thus would be expected to exhibit more differences between samples. All of the 189 samples in this study were individualized with a subset of 12 SNP loci. However, it is likely that the addition of more than 12 SNP loci will be required to resolve larger sets of unrelated individuals from one another. By way of comparison, in these same 189 individuals all but one pair is resolved from one another with three of the traditional short tandem repeat (STR) loci possessing the highest heterozygosity values (D2S1338, D18S51, and FGA) run with the Identifiler kit. The final pair of unrelated samples could be resolved with the combination of 4 STR loci: D2S1338, D18S51, FGA, and VWA.

A graphical simulation model of the entire DNA process associated with the analysis of short tandem repeat loci

The use of expert systems to interpret short tandem repeat DNA profiles in forensic, medical and ancient DNA applications is becoming increasingly prevalent as high-throughput analytical systems generate large amounts of data that are time-consuming to process. With special reference to low copy number (LCN) applications, we use a graphical model to simulate stochastic variation associated with the entire DNA process starting with extraction of sample, followed by the processing associated with the preparation of a PCR reaction mixture and PCR itself. Each part of the process is modelled with input efficiency parameters. Then, the key output parameters that define the characteristics of a DNA profile are derived, namely heterozygote balance (Hb) and the probability of allelic drop-out p(D). The model can be used to estimate the unknown efficiency parameters, such as π_extraction. ‘What-if’ scenarios can be used to improve and optimize the entire process, e.g. by increasing the aliquot forwarded to PCR, the improvement expected to a given DNA profile can be reliably predicted. We demonstrate that Hb and drop-out are mainly a function of stochastic effect of pre-PCR molecular selection. Whole genome amplification is unlikely to give any benefit over conventional PCR for LCN.

Validation of a 21-locus autosomal SNP multiplex for forensic identification purposes

A single nucleotide polymorphism (SNP) multiplex has been developed to analyse highly degraded and low copy number (LCN) DNA template, i.e. <100 pg, for scenarios including mass disaster identification. The multiplex consists of 20 autosomal non-coding loci plus Amelogenin for sex determination, amplified in a single tube PCR reaction and visualised on the Applied Biosystems 3100 capillary electrophoresis (CE) system. Allele-specific primers tailed with shared universal tag sequences were designed to speed multiplex design and balance the amplification efficiencies of all loci through the use of a single reverse and two differentially labelled allele denoting forward universal primers. As the multiplex is intended for use with samples too degraded for conventional profiling, a computer program was specifically developed to aid interpretation. Critical factors taken into account by the software include empirically determined extremes of heterozygous imbalance (Hb) and the drop-out threshold (Ht) defined as the maximum peak height of a surviving heterozygous allele, where its partner may have dropped out. The discrimination power of the system is estimated at 1 in 4.5 million, using a White Caucasian population database. Comparisons using artificially degraded samples profiled with both the SNP multiplex and AMPFISTR SGM plus (Applied Biosystems) demonstrated a greater likelihood of obtaining a profile using SNPs for certain sample types. Saliva stains degraded for 147 days generated an 81% complete SNP profile whilst short tandem repeats (STRs) were only 18% complete; similarly blood degraded for 243 days produced full SNP profiles but only 9% with STRs. Reproducibility studies showed concordance between SNP profiles for different sample types, such as blood, saliva, semen and hairs, for the same individual, both within and between different DNA extracts.