Efficacy of extended kinship analyses utilizing commercial STR kit in establishing personal identification

Unprecedented fidelity and specificity have afforded DNA testing its long reigning status as the gold standard for establishing personal identification. While the method itself is flawless, forensic experts have undoubtedly stumbled across challenging cases in which no reference samples for an unknown person (UP) are available for comparison. In such cases, experts often must resort to an assortment of kinship analyses-primarily those involving alleged parents or children of a UP-to establish personal identification. The present study derives likelihood ratio (LR) distributions from an extensive series of kinship simulations and places actual data, obtained from 120 cases in which personal identification of a UP was established via kinship analyses, to a comprehensive comparison in order to evaluate the efficacy of kinship assessments in establishing personal identification. A commercially available AmpFlSTR Identifiler kit was used to obtain DNA profiles. UP DNAs were extracted and isolated from fingernail (n=87), cardiac blood (24), carpal bone (7) and tooth (2). Buccal cells were procured from alleged kin (AK) for subsequent kinship analyses. In 72 cases 1-3 alleged children were available for comparison; in 46 cases, one or both alleged parents were available; and in the final 2 cases (involving a pair of bodies discovered together in a dwelling), their alleged children were typed for comparison. For each case a LR was calculated based on the DNA typing results. Interestingly, we found that the median LR observed in the actual cases virtually mirrored those of the simulations. With exception to 2 cases in which a silent allele was observed at D19S433, biological relatives showed a LR greater than 100 and in these cases, kinship between the UP and AK were further supported by additional forms of evidence. We show here that in the vast majority of identification cases where direct reference samples are unavailable for a UP, kinship analyses referring to alleged parents/children and using 15 standard loci is more than capable of establishing the identification of a UP. However, discretion should be advised for silent alleles which-albeit rare-are known to occur at loci such as D19S433, along with other mutations which could render a deceivingly reduced LR.     

Next-generation sequencing analysis of off-ladder alleles due to migration shift caused by sequence variation at D12S391 locus

In short tandem repeat (STR) analysis, length polymorphisms are detected by capillary electrophoresis (CE). At most STR loci, mobility shift due to sequence variation in the repeat region was thought not to affect the typing results. In our recent population studies of 1501 Japanese individuals, off-ladder calls were observed at the D12S391 locus using PowerPlex Fusion in nine samples for allele 22, one sample for allele 25, and one sample for allele 26. However, these samples were typed as ordinary alleles within the bins using GlobalFiler. In this study, next-generation sequencing analysis using MiSeq was performed for the D12S391 locus from the 11 off-ladder samples and 33 other samples, as well as the allelic ladders of PowerPlex Fusion and GlobalFiler. All off-ladder allele 22 in the nine samples had [AGAT]₁₁[AGAC]₁₁ as a repeat structure, while the corresponding allele was [AGAT]₁₅[AGAC]₆[AGAT] for the PowerPlex Fusion ladder, and [AGAT]₁₃[AGAC]₉ for the GlobalFiler ladder. Overall, as the number of [AGAT] in the repeat structure decreased at the D12S391 locus, the peak migrated more slowly using PowerPlex Fusion, the reverse strand of which was labeled, and it migrated more rapidly using GlobalFiler, the forward strand of which was labeled. The allelic ladders of both STR kits were reamplified with our small amplicon D12S391 primers and their mobility was also examined. In conclusion, off-ladder observations of allele 22 at the D12S391 locus using PowerPlex Fusion were mainly attributed to a relatively large difference of the repeat structure between its allelic ladder and off-ladder allele 22.     

Massively parallel sequencing data of 31 autosomal STR loci obtained using the Precision ID GlobalFiler NGS STR Panel v2 for 82 Japanese population samples

Allele frequencies for 31 autosomal short tandem repeat (STR) loci (CSF1PO, D10S1248, D12ATA63, D12S391, D13S317, D14S1434, D16S539, D18S51, D19S433, D1S1656, D1S1677, D21S11, D22S1045, D2S1338, D2S1776, D2S441, D3S1358, D3S4529, D4S2408, D5S2800, D5S818, D6S1043, D6S474, D7S820, D8S1179, FGA, Penta D, Penta E, TH01, TPOX, and vWA) were obtained using Precision ID GlobalFiler NGS STR Panel v2 from 82 unrelated individuals sampled from the Japanese population. Autosomal STR alleles designated by NGS and conventional capillary electrophoresis were found to be concordant except at D2S441 allele 9.     

Developmental validation of the HomyGene19+14Y System

The HomyGene19+14Y System (HG19+14Y) is a PCR-based amplification kit that enables typing of 18 autosomal short tandem repeat (STR) loci (i.e., CSF1PO, D2S1338, D3S1358, D5S818, D6S1043, D7S820, D8S1179, D12S391, D13S317, D16S539, D18S51, D19S433, D21S11, FGA, Penta E, TPOX, TH01, vWA), 14 widely used Y chromosome STR (Y-STR) loci (Y_GATA_H4, DYS385a/b, DYS389I/II, DYS390, DYS391, DYS392, DYS393, DYS438, DYS439, DYS456, DYS458, DYS635), and amelogenin. This multiplex system was designed for the simultaneous analysis of amelogenin-Y allele mutation, single-source searches, kinship (including familial searching), mixture profiles, international data sharing, and other forensic applications. In this study, the multiplex system was validated for sensitivity, specificity, DNA mixtures, stability, precision, stutter, reproducibility, parallel tests, PCR-based conditions, and population analysis according to the Scientific Working Group on DNA Analysis Methods (SWGDAM) developmental validation guidelines. A total of 212 alleles were detected for the 18 autosomal STR loci among 528 Guangdong Han individuals, and 431 haplotypes were found for 14 Y-STRs among 452 unrelated males. The combined match probability (CMP) of the HG19+14Y System was calculated as 2.39 × 10⁻²⁹. All the validation results showed that the HG19+14Y System would be a robust, reliable, highly polymorphic, and informative forensic kit.

     

Null alleles and sequence variations at primer binding sites of STR loci within multiplex typing systems

Rare variants are widely observed in human genome and sequence variations at primer binding sites might impair the process of PCR amplification resulting in dropouts of alleles, named as null alleles. In this study, 5 cases from routine paternity testing using PowerPlex^®21 System for STR genotyping were considered to harbor null alleles at TH01, FGA, D5S818, D8S1179, and D16S539, respectively. The dropout of alleles was confirmed by using alternative commercial kits AGCU Expressmarker 22 PCR amplification kit and AmpFℓSTR^®. Identifiler^® Plus Kit, and sequencing results revealed a single base variation at the primer binding site of each STR locus. Results from the collection of previous reports show that null alleles at D5S818 were frequently observed in population detected by two PowerPlex^® typing systems and null alleles at D19S433 were mostly observed in Japanese population detected by two AmpFℓSTR™ typing systems. Furthermore, the most popular mutation type appeared the transition from C to T with G to A, which might have a potential relationship with DNA methylation. Altogether, these results can provide helpful information in forensic practice to the elimination of genotyping discrepancy and the development of primer sets.     

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.     

Long D13S317 variant allele: A cautionary case report

Exclusion at locus D13S317 between a father and child was observed in a kinship case, which at first glance appeared as a silent allele. However, a closer inspection using three commercial kits showed that the observed pattern is due to a long variant allele overlapping with different loci in different kits. Sequencing of the variant revealed a duplication within D13S317, that had also created an additional binding site for short amplicon reverse primer. If ignored, genotyping results including this variant may be wrongly interpreted.     

Sibling assessment based on likelihood ratio and total number of shared alleles using 21 short tandem repeat loci included in the GlobalFiler™ kit

Sibling assessment using the 15 autosomal short tandem repeat (STR) loci included in the Identifiler® kit can be difficult when comparing an unidentified party to an alleged sibling. Therefore, we investigated the likelihood ratio (LR) and the total number of shared alleles (TNSA) for sibship determination using the 21 autosomal STR loci included in the GlobalFiler™ kit.We computationally generated the genotypes of 10,000 sibling pairs and 10,000 unrelated pairs based on previously reported allele frequencies of the 15 Identifiler loci and the remaining 6 GlobalFiler loci. The LR and the TNSA were then calculated in each pair using the 15 and 21 loci. Next, these calculations were applied to 22 actual sibling pairs.LR values ⩾10,000 were observed in 48% of the sibling pairs using the 15 loci and in 80% of the sibling pairs using the 21 loci. The TNSA distribution between siblings and unrelated pairs was more divergent in GlobalFiler than in Identifiler. TNSA values ⩾20 were found only in true siblings in Identifiler, while TNSA values ⩾24 in GlobalFiler. In Identifiler, all pairs with TNSA ⩾24 had LR values ⩾10,000 and the same was true in GlobalFiler for TNSA ⩾29. Therefore, increasing the number of loci is very efficient for sibship determination.The LR is most reliable for determining sibship. However, TNSA values may be useful for the preliminary method of LR values because LR value demonstrated a significantly positive correlation with TNSA value in both Identifiler and GlobalFiler.     

Characterisation of variant alleles in the STR systems D2S1338, D3S1358 and D19S433

We have observed three hitherto undescribed off-ladder alleles at three widely used STR loci. These were isolated, sequenced and designated as follows: allele 10 (D2S1338, one case), allele 21 (D3S1358, two cases) and allele 6.2 (D19S433, six cases). These sequences are described in comparison to non-variant alleles, and their implications for the semi-automated STR analysis will be discussed.     

Allele frequencies for 22 autosomal short tandem repeat loci obtained by PowerPlex Fusion in a sample of 1501 individuals from the Japanese population

Allele frequencies for 22 autosomal short tandem repeat loci (D3S1358, D1S1656, D2S441, D10S1248, D13S317, Penta E, D16S539, D18S51, D2S1338, CSF1PO, Penta D, TH01, vWA, D21S11, D7S820, D5S818, TPOX, D8S1179, D12S391, D19S433, FGA, and D22S1045) were obtained from 1501 unrelated individuals sampled from the Japanese population.     

Genetic polymorphism studies on 22 autosomal STR loci of the PowerPlex Fusion System in Bangladeshi population

Genetic polymorphism of 22 autosomal STR loci included in PowerPlex® Fusion System (D3S1358, D1S1656, D2S441, D10S1248, D13S317, Penta E, D16S539, D18S51, D2S1338, CSF1PO, Penta D, TH01, vWA, D21S11, D7S820, D5S818, TPOX, D8S1179, D12S391, D19S433, FGA and D22S1045) was studied in 188 unrelated Bangladeshi Bengali individuals. Allele frequencies and forensic efficiency parameters such as, the power of discrimination (PD), observed and expected heterozygosity (Ho & He), polymorphism information content (PIC), probability of match (PM), power of exclusion (PE) and typical paternity index was calculated for the loci. The combined PM and PE for all 22 STR loci were calculated to be 5.29 × 10⁻²⁷ and 0.99999999945 respectively. The dataset indicated the usefulness of these loci in personal identification, parentage testing and complex kinship analysis in Bangladeshi population. A neighbor-joining tree was constructed based on pair-wise Nei’s genetic distance by comparing allele frequency data for the 22 loci with six other populations. The analysis showed that Bangladeshi population lies closer to a clade consisting Japan, the Philippines and East Timot populations.     

Population genetics of the hypervariable locus D12S391 in Koreans

The hypervariable short tandem repeat (STR) locus D12S391 was investigated in a Korean population and 34 fragments were sequenced to confirm the structure of alleles. From these sequenced fragments an allelic ladder containing 13 sequenced alleles was constructed. From 595 unrelated Koreans, 14 alleles were detected and one variant allele 19.3 was observed. The observed heterozygosity was 0.795 and no deviation from Hardy-Weinberg equilibrium was observed in the Korean population (p = 0.606). The allele frequency distribution in the Korean population was not similar to other racial or ethnic groups except for Egyptians, Yemenis, Japanese and Caucasoids from the Rhine area. No mutations were observed in the 702 meioses from 144 Korean families. This study demonstrates that the STR locus D12S391 is a useful tool for forensic identification and parentage testing. Received: 15 September 1999 / Accepted: 18 December 1999     

D20S161 data for three ethnic populations and forensic validation

In order to evaluate the forensic applicability of the STR locus D20S161 and construct a preliminary database, the genotype distributions and allele frequencies in five populations from three main ethnic groups were investigated, including Germans, Slovakians, African Americans, Japanese and Chinese. A total of 512 samples from unrelated individuals and 85 confirmed father/mother/¶child triplets were analyzed by PCR and allele determination was carried out by comparison with a sequenced human allelic ladder. The results showed that D20S161 typing was both precise and reliable. A total of 7 alleles was found in these populations and no evidence of deviation from Hardy-Weinberg equilibrium was observed. Pairwise comparisons between populations showed that there were significant differences in the distributions of the allele frequencies among the three main ethnic groups. No mutation events were observed from the confirmed father/mother/¶child triplets. With a maximum likelihood method, the mutation rate was indirectly estimated as 2.5 × 10^–5. These results suggest that D20S161 is a useful marker for forensic casework and paternity analysis.     

STR allele sequence variation: Current knowledge and future issues

This article reviews what is currently known about short tandem repeat (STR) allelic sequence variation in and around the twenty-four loci most commonly used throughout the world to perform forensic DNA investigations. These STR loci include D1S1656, TPOX, D2S441, D2S1338, D3S1358, FGA, CSF1PO, D5S818, SE33, D6S1043, D7S820, D8S1179, D10S1248, TH01, vWA, D12S391, D13S317, Penta E, D16S539, D18S51, D19S433, D21S11, Penta D, and D22S1045. All known reported variant alleles are compiled along with genomic information available from GenBank, dbSNP, and the 1000 Genomes Project. Supplementary files are included which provide annotated reference sequences for each STR locus, characterize genomic variation around the STR repeat region, and compare alleles present in currently available STR kit allelic ladders. Looking to the future, STR allele nomenclature options are discussed as they relate to next generation sequencing efforts underway.     

Effect of linkage between vWA and D12S391 in kinship analysis

Ideally for use in forensic analyses, genetic markers on the same chromosome should be more than 50 Mb in physical distance to ensure full recombination and thus independent inheritance. The forensic community has given attention to two STR markers, D12S391 and vWA, that are 6.3 megabases (Mb) apart on chromosome 12. Recent studies have shown no significant linkage disequilibrium between vWA and D12S391 in U.S. and worldwide populations, although genetic linkage has been identified. It is important to evaluate the impact of linkage effects on kinship analysis. In this study, we aimed to determine a more precise measurement of the recombination frequency between vWA and D12S391 based on a larger number of informative meiosis than has been studied previously. We estimated the recombination frequency (θ) to 0.089 (95% CI 0.044–0.158). Using pedigrees simulated under specific kinship scenarios where recombination was expected to affect the likelihood ratio (LR), we evaluated the impact on LR values of including or ignoring linkage between vWA and D12S391. For all pedigree scenarios considered, on average, LR values ignoring linkage were slightly underestimated than when linkage was considered. However, in the incest scenario considered, LR values could be overestimated up to 25–30 times when linkage was ignored. We demonstrate that the effect of ignoring linkage in the likelihood ratio calculation can be considerable. These results suggest that linkage should be considered during kinship analysis when vWA and D12S391 are tested for pedigrees where a recombination could impact the LR value.     

Identification of a rare off-ladder allele of the D13S325 locus during paternity testing

This study demonstrates an unusual rare allele of D13S325 that was falsely categorized as an allele of D12S391 under the STRtyper™-10F/G system. The parentage cases with these rare alleles were analyzed using the Sinofiler™ system and singleplex amplification system, and the alleles of D13S325 extracted from the electrophoresis gel were sequenced. 5 Cases with the rare alleles misread as allele 20 of D12S391 were identified in total 2618 cases (including 3200 unrelated parents). This rare allele was designated as allele 5.1 of D13S325 based on its DNA sequence. Its frequency in the Chinese population was 1.6 × 10⁻³. Because the rare allele 5.1 of D13S325 locus tends to be incorrectly labeled in the STRtyper™-10F/G system, particular attention should be paid when the system is used in paternity testing, personal identification, and DNA database comparisons.     

Analysis of eight STR loci in two Hungarian populations

A multiplex reaction for the eight STR loci D3S1358, FGA, D8S1179, D21S11, D18S51, D5S818, D13S317, D7S820 was used to generate allele frequency databases for two Hungarian population samples, Caucasians from the Budapest area and Romanies from Baranya county. During the analysis two intermediate-sized alleles and a sequence variant allele were observed at the D7S820 locus. All three types of allelic variants were found to have modifications in the same block of a (T)₉ stretch located within the 3′ flanking region of each allele, which may indicate a possible higher mutation rate of this (T)₉ block. For the loci D3S1358 and D7S820 the Romany population database showed departures from Hardy-Weinberg equilibrium. The forensic efficiency values for the Romany population were slightly different from those found in the Hungarian Caucasian population. Comparing the allele frequency values by G-statistic, calculating the F_ST indices and with the pair-wise comparisons of inter-population variance, the two Hungarian populations could be distinguished using data from the eight STR loci. Received: 18 May 1999 / Received in revised form: 26 August 1999     

Polymorphisms at fifteen tetrameric short tandem repeat loci in three ethnic populations of Bengal, India

The study presents allele frequency data at 15 tetrameric short tandem repeat (STR) loci (D3S1358, THO1, D21S11, D18S51, D5S818, D13S317, D7S820, D16S539, CSF1PO, vWA, D8S1179, TPOX, D2S1338, D19S433 and FGA) in three ethnic populations--Mahishya, Bauri and Namasudra of Bengal to evaluate their utility in Forensic testing and understanding population structure and dynamics. A total of 169 individuals were studied from the selected populations. On an average the combined power of discrimination and power of exclusion in these groups was found 0.97 and 0.99, respectively. The allele distribution pattern shows possible genetic admixture between these ethnic groups which could be attributed to their close geographical proximity and occupying almost similar position in the social hierarchy. This study suggests that the 13 Combined DNA Index System (CODIS) markers and two added markers named D2S1338, D19S433 are highly informative and therefore suitable in matching biological specimen in human identification and population genetic study.     

A silent allele in the locus D19S433 contained within the AmpFℓSTR® Identifiler™ PCR Amplification Kit

We present two cases where a single locus mismatch was found in the locus D19S433 using the AmpF?STR® Identifiler? PCR Amplification Kit (Applied Biosystems) (Identifiler Kit) during paternity and maternity tests. This mismatch differed from the mismatch pattern where there is usually a one repeat difference. We designed forward and reverse primers so that they were positioned further away from the primer set contained in the Identifiler Kit. The results showed the existence of a silent allele 13 in both families, due to a point mutation that changed guanine to adenine at 32 nucleotides downstream from the 3′ end of the AAGG repeat sequences in all four members. A single locus mismatch due to a silent allele may occur in any locus using any kit. Accordingly, we should pay attention to this silent allele when carrying out human identification and parentage analysis.     

Genetic variation of microsatellite markers D1S117, D6S89, D11S35, APOC2, and D21S168 in the Spanish population

Summary We have used PCR amplification to analyse the allele frequency, distribution and heterozygosity of 5 microsatellite markers (D1S117, D6S89, D11S35, APOC2, and D21S168), in a sample of 100 unrelated Spanish individuals. The loci tested exhibit wide allelic variability having 7-17 alleles, PIC (polymorphic information content) between 0.79 and 0.86, and heterozygosity between 0.81 and 0.86. D1S117 and D21S168 have unimodal distribution, APOC2 has 4 common alleles which account for 71% of the total variation, D11S35 has a bimodal distribution and D6S89 is trimodal. The allelic distribution observed for each locus is in agreement with slippage and mispairing as the main mechanisms involved in the evolution of microsatellite alleles. Multiplex amplification of loci D6S89 and APOC2 was possible due to their non-overlapping allele sizes. The rapidity with which microsatellites can be analysed, and the accurate determination of alleles, make these markers very powerful tools for genetic typing. The information obtained for loci D1S117, D6S89, D11S35, APOC2, and D21S168, provides a basis for their use for DNA typing and paternity analysis in the Spanish population.