M.V. Lareu

Inferring ancestral origin using a single multiplex assay of ancestry-informative marker SNPs

Tests that infer the ancestral origin of a DNA sample have considerable potential in the development of forensic tools that can help to guide crime investigation. We have developed a single-tube 34-plex SNP assay for the assignment of ancestral origin by choosing ancestry-informative markers (AIMs) exhibiting highly contrasting allele frequency distributions between the three major population-groups. To predict ancestral origin from the profiles obtained, a classification algorithm was developed based on maximum likelihood. Sampling of two populations each from African, European and East Asian groups provided training sets for the algorithm and this was tested using the CEPH Human Genome Diversity Panel. We detected negligible theoretical and practical error for assignments to one of the three groups analyzed with consistently high classification probabilities, even when using reduced subsets of SNPs. This study shows that by choosing SNPs exhibiting marked allele frequency differences between population-groups a practical forensic test for assigning the most likely ancestry can be achieved from a single multiplexed assay.

Hierarchical analysis of 30 Y-chromosome SNPs in European populations

Analysis of Y-chromosome haplogroups defined by binary polymorphisms, has became a standard approach for studying the origin of modern human populations and for measuring the variability between them. Furthermore, the simplicity and population specificity of binary polymorphisms allows inferences to be drawn about the population origin of any male sample of interest for forensic purposes. From the 245 binary polymorphisms that can be analysed by PCR described in the Y Chromosome Consortium tree, we have selected 30 markers. The set of 30 has been grouped into 4 multiplexes in order to determine the most frequent haplogroups in Europe, using only 1 or 2 multiplexes. In this way, we avoid typing unnecessary SNPs to define the final haplogroup saving effort and cost, since we only need to type 9 SNPs in the best case and in the worst case, no more than 17 SNPs to define the haplogroup. The selected method for allele discrimination was a single base extension reaction using the SNaPshot multiplex kit. A total of 292 samples from 8 different districts of Galicia (northwest Spain) were analysed with this strategy. No significant differences were detected among the different districts, except for the population from Mariña Lucense, which showed a distant haplogroup frequency but not higher Φst values.

Resolving relationship tests that show ambiguous STR results using autosomal SNPs as supplementary markers

When using a standard battery of STRs for relationship testing a small proportion of analyses can give ambiguous results - where the claimed relationship cannot be confirmed by a high enough paternity index or excluded with fully incompatible genotypes. The majority of such cases arise from unknowingly testing a brother of the true father and observing only a small number of exclusions that can each be interpreted as one- or two-step mutations. Although adding extra STRs might resolve a proportion of cases, there are few properly validated extra STRs available, while the commonly added hypervariable SE33 locus is four times more mutable than average, increasing the risk of ambiguous results. We have found SNPs in large multiplexes are much more informative for both low initial probabilities or ambiguous exclusions and at the same time provide a more reliable genotyping approach for the highly degraded DNA encountered in many identification cases. Eight relationship cases are outlined where the addition of SNP data resolved analyses that had remained ambiguous even with extended STR typing. In addition we have made simulations to ascertain the frequency of failing to obtain exclusions or conclusive probabilities of paternity with different marker sets when a brother of the true father is tested. Results indicate that SNPs are statistically more efficient than STRs in resolving cases that distinguish first-degree relatives in deficient pedigrees.

Analysis of global variability in 15 established and 5 new European Standard Set (ESS) STRs using the CEPH human genome diversity panel

The CEPH human genome diversity cell line panel (CEPH-HGDP) of 51 globally distributed populations was used to analyze patterns of variability in 20 core human identification STRs. The markers typed comprised the 15 STRs of Identifiler, one of the most widely used forensic STR multiplexes, plus five recently introduced European Standard Set (ESS) STRs: D1S1656, D2S441, D10S1248, D12S391 and D22S1045. From the genotypes obtained for the ESS STRs we identified rare, intermediate or off-ladder alleles that had not been previously reported for these loci. Examples of novel ESS STR alleles found were characterized by sequence analysis. This revealed extensive repeat structure variation in three ESS STRs, with D12S391 showing particularly high variability for tandem runs of AGAT and AGAC repeat units. The global geographic distribution of the CEPH panel samples gave an opportunity to study in detail the extent of substructure shown by the 20 STRs amongst populations and between their parent population groups. An assessment was made of the forensic informativeness of the new ESS STRs compared to the loci they will replace: CSF1PO, D5S818, D7S820, D13S317 and TPOX, with results showing a clear enhancement of discrimination power using multiplexes that genotype the new ESS loci. We also measured the ability of Identifiler and ESS STRs to infer the ancestry of the CEPH-HGDP samples and demonstrate that forensic STRs in large multiplexes have the potential to differentiate the major population groups but only with sufficient reliability when used with other ancestry-informative markers such as single nucleotide polymorphisms. Finally we checked for possible association by linkage between the two ESS multiplex STRs closely positioned on chromosome-12: vWA and D12S391 by examining paired genotypes from the complete CEPH data set.

Investigation of the STR locus HUMTH01 using PCR and two electrophoresis formats: UK and Galician Caucasian population surveys and usefulness in paternity investigations

Two hundred unrelated Caucasians from the UK and 210 from Galicia (NW Spain) have been genotyped for the HUMTH01 locus using polymerase chain reaction (PCR) amplification followed by high sieving agarose electrophoresis (UK samples) and polyacrylamide electrophoresis (Galician samples). Allele and genotype frequencies obtained from both populations were in close agreement to those seen by other workers and chi 2 tests of the two populations showed that both were in Hardy-Weinberg equilibrium. The potential usefulness of HUMTH01 in paternity investigations was analysed by constructing, within each population, false family trios, each with a non-father, in order to estimate the actual exclusion rate. The observed rate of exclusion was in close agreement with the expected rate for both populations. Examination of the mother-child pairs in the false families showed a common allele in every case and no evidence of mutation or non-Mendelian inheritance was observed. The HUMTH01 locus shows an informative polymorphism and the production and analysis of population databases is easily achieved. The use of PCR followed by electrophoresis in either gel format appears to provide a quick and straightforward method for investigating the HUMTH01 locus.

Heteroplasmy in mtDNA and the weight of evidence in forensic mtDNA analysis: a case report

Abstract Mitochondrial DNA (mtDNA) sequencing has been validated as a useful tool for forensic analysis. However, there are several aspects of the analysis which need to be considered in order to evaluate the value of the evidence. One of these aspects is related to heteroplasmy which is the state when two or more mtDNA populations occur in a single individual, cell or mitochondrion. In this report a case is described where the mtDNA profile of the blood sample of a raped woman was compared with the mtDNA profile of a single hair found in the suspect’s car. The results obtained show differences in sequence between different portions of the hair and the victim’s sequence. These differences are related to various heteroplasmy events. The concordance between the hair sample and the potential source (victim) of this sample is questionable and the strength of the evidence depends on how the sequence information is interpreted by the expert. The discussion of the results emphasises the necessity to evaluate heteroplasmic events in routine forensic work.

Extending STR markers in Y chromosome haplotypes

Abstract. Two multiplex reactions were developed to amplify 16 Y-STRs (DYS19, DYS385, DYS389 I and II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, GATA A7.1, GATA A7.2, GATA A10, GATA C4, GATA H4). Here we extend previous population studies done in a sample from northern Portugal for the GATA A7.1, GATA A7.2, GATA C4 and GATA H4 loci. A total of 199 different haplotypes identified by the 16 Y-STR markers were observed in a sample of 208 male individuals, of which 190 were unique and 9 were found twice. The overall haplotype diversity was 0.9996. The haplotype diversity of the Y-STR set composed of the 8 new markers is higher than the Y-STR core set included in the Y-STR haplotype reference database. Sequence structure of new alleles for GATA C4 and GATA H4 is reported. The usefulness of the inclusion of this new set of Y-STRs in forensic casework was also assessed. The increase in haplotype diversity with the addition of any new Y-STR marker to the 8 Y-STR core set is dependent not only on the gene diversity (positively) but also (negatively) on the degree of gametic association between the markers and the haplotypes previously defined. For instance, in our sample the addition of the DYS437, DYS438 and GATA A7.2 to a 13-locus set increased haplotype diversity only by 0.0001.

Robustness of the Y STRs DYS19, DYS389 I and II, DYS390 and DYS393: optimization of a PCR pentaplex

Various technical methods were investigated with the aim of developing a multiplex system to amplify five Y-chromosome STR loci in the same PCR reaction: DYS393, DYS19, DYS390, DYS389 I and DYS389 II. A sequenced allelic ladder was constructed with previously sequenced alleles including the most common ones. A number of reamplification conditions of the allelic ladders were tested. The pentaplex was evaluated for typing using two different platforms (ABI and ALF) with promising results. However, in degraded samples non-specific artifacts were observed in the DYS393 system in the same range of sizes as the real alleles. This system can also be typed in females under relatively low stringency conditions in the PCR amplification, making this system prone to errors in critical samples. This lack of specificity can be reduced by increasing the stringency of the PCR conditions. The DYS19 ladder cannot be reamplified as stutters appear after a few reamplifications. These stutters are probably due to a 2 bp slippage induced by the presence of a TA repeat stretch in the PCR amplified fragments. Non-specific products were also noted in the DYS389 I and DYS389 II amplification, although out of the range of other alleles in this pentaplex. This newly constructed pentaplex has proved to be very useful in population genetic studies because all five Y STR markers can be loaded in the same lane of a gel with other Y STR singleplex or multiplexes. The usefulness of Y-chromosome STRs in criminal casework is especially evident in analyzing azoospermic individuals.

SEQUENCE VARIATION OF A HYPERVARIABLE SHORT TANDEM REPEAT AT THE D1S1656 LOCUS

Abstract A short tandem repeat at the D1S1656 locus was sequenced in 45 selected alleles and 13 different alleles were found which were designated according to the total number of repeats. This STR is a compound hypervariable STR consisting of blocks of (TAGR) repeats with a basic sequence structure (TAGA)4(TGA)0-1(TAGA)6-16- (TAGG)0-1(TG)5. The presence of a TGA, probably due to an A deletion in the fifth TAGA repeat leads to intermediate a.3 alleles. Population data showed that this is a highly polymorphic STR with a heterozygosity of more than 0.89. This fact together with its simple structure and small size (129–168 bp) makes this STR one of the most interesting DNA polymorphisms for forensic and genetic purposes.

Further development of forensic eye color predictive tests.

In forensic analysis predictive tests for external visible characteristics (or EVCs), including inference of iris color, represent a potentially useful tool to guide criminal investigations. Two recent studies, both focused on forensic testing, have analyzed single nucleotide polymorphism (SNP) genotypes underlying common eye color variation (Mengel-From et al., Forensic Sci. Int. Genet. 4:323 and Walsh et al., Forensic Sci. Int. Genet. 5:170). Each study arrived at different recommendations for eye color predictive tests aiming to type the most closely associated SNPs, although both confirmed rs12913832 in HERC2 as the key predictor, widely recognized as the most strongly associated marker with blue and brown iris colors. Differences between these two studies in identification of other eye color predictors may partly arise from varying approaches to assigning phenotypes, notably those not unequivocally blue or dark brown and therefore occupying an intermediate iris color continuum. We have developed two single base extension assays typing 37 SNPs in pigmentation-associated genes to study SNP-genotype based prediction of eye, skin, and hair color variation. These assays were used to test the performance of different sets of eye color predictors in 416 subjects from six populations of north and south Europe. The presence of a complex and continuous range of intermediate phenotypes distinct from blue and brown eye colors was confirmed by establishing eye color populations compared to genetic clusters defined using Structure software. Our study explored the effect of an expanded SNP combination beyond six markers has on the ability to predict eye color in a forensic test without extending the SNP assay excessively - thus maintaining a balance between the tests predictive value and an ability to reliably type challenging DNA with a multiplex of manageable size. Our evaluation used AUC analysis (area under the receiver operating characteristic curves) and naïve Bayesian likelihood-based classification approaches. To provide flexibility in SNP-based eye color predictive tests in forensic applications we modified an online Bayesian classifier, originally developed for genetic ancestry analysis, to provide a straightforward system to assign eye color likelihoods from a SNP profile combining additional informative markers from the predictors analyzed by our study plus those of Walsh and Mengel-From. Two advantages of the online classifier is the ability to submit incomplete SNP profiles, a common occurrence when typing challenging DNA, and the ability to handle physically linked SNPs showing independent effect, by allowing the user to input frequencies from SNP pairs or larger combinations. This system was used to include the submission of frequency data for the SNP pair rs12913832 and rs1129038: indicated by our study to be the two SNPs most closely associated to eye color.