Multi-laboratory validation of DNAxs including the statistical library DNAStatistX

This study describes a multi-laboratory validation of DNAxs, a DNA eXpert System for the data management and probabilistic interpretation of DNA profiles [1], and its statistical library DNAStatistX to which, besides the organising laboratory, four laboratories participated. The software was modified to read multiple data formats and the study was performed prior to the release of the software to the forensic community. The first exercise explored all main functionalities of DNAxs with feedback on user-friendliness, installation and general performance. Next, every laboratory performed likelihood ratio (LR) calculations using their own dataset and a dataset provided by the organising laboratory. The organising laboratory performed LR calculations using all datasets. The datasets were generated with different STR typing kits or analysis systems and consisted of samples varying in DNA amounts, mixture ratios, number of contributors and drop-out level. Hypothesis sets had the correct, under- and over-assigned number of contributors and true and false donors as person of interest. When comparing the results between laboratories, the LRs were foremost within one unit on log10 scale. The few LR results that deviated more had differences for the parameters estimated by the optimizer within DNAStatistX. Some of these were indicated by failed iteration results, others by a failed model validation, since unrealistic hypotheses were included. When these results that do not meet the quality criteria were excluded, as is in accordance with interpretation guidelines, none of the analyses in the different laboratories yielded a different statement in the casework report. Nonetheless, changes in software parameters were sought that minimized differences in outcomes, which made the DNAStatistX module more robust. Overall, the software was found intuitive, user-friendly and valid for use in multiple laboratories.

A global analysis of Y-chromosomal haplotype diversity for 23 STR loci

In a worldwide collaborative effort, 19,630 Y-chromosomes were sampled from 129 different populations in 51 countries. These chromosomes were typed for 23 short-tandem repeat (STR) loci (DYS19, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS385ab, DYS437, DYS438, DYS439, DYS448, DYS456, DYS458, DYS635, GATAH4, DYS481, DYS533, DYS549, DYS570, DYS576, and DYS643) and using the PowerPlex Y23 System (PPY23, Promega Corporation, Madison, WI). Locus-specific allelic spectra of these markers were determined and a consistently high level of allelic diversity was observed. A considerable number of null, duplicate and off-ladder alleles were revealed. Standard single-locus and haplotype-based parameters were calculated and compared between subsets of Y-STR markers established for forensic casework. The PPY23 marker set provides substantially stronger discriminatory power than other available kits but at the same time reveals the same general patterns of population structure as other marker sets. A strong correlation was observed between the number of Y-STRs included in a marker set and some of the forensic parameters under study. Interestingly a weak but consistent trend toward smaller genetic distances resulting from larger numbers of markers became apparent.

Standing at the gateway to Europe--the genetic structure of Western balkan populations based on autosomal and haploid markers

Contemporary inhabitants of the Balkan Peninsula belong to several ethnic groups of diverse cultural background. In this study, three ethnic groups from Bosnia and Herzegovina - Bosniacs, Bosnian Croats and Bosnian Serbs - as well as the populations of Serbians, Croatians, Macedonians from the former Yugoslav Republic of Macedonia, Montenegrins and Kosovars have been characterized for the genetic variation of 660 000 genome-wide autosomal single nucleotide polymorphisms and for haploid markers. New autosomal data of the 70 individuals together with previously published data of 20 individuals from the populations of the Western Balkan region in a context of 695 samples of global range have been analysed. Comparison of the variation data of autosomal and haploid lineages of the studied Western Balkan populations reveals a concordance of the data in both sets and the genetic uniformity of the studied populations, especially of Western South-Slavic speakers. The genetic variation of Western Balkan populations reveals the continuity between the Middle East and Europe via the Balkan region and supports the scenario that one of the major routes of ancient gene flows and admixture went through the Balkan Peninsula.

RNA/DNA co-analysis from human menstrual blood and vaginal secretion stains: results of a fourth and fifth collaborative EDNAP exercise

The European DNA Profiling Group (EDNAP) organized a fourth and fifth collaborative exercise on RNA/DNA co-analysis for body fluid identification and STR profiling. The task was to identify dried menstrual blood and vaginal secretion stains using specific RNA biomarkers, and additionally test 3 housekeeping genes for their suitability as reference genes. Six menstrual blood and six vaginal secretion stains, two dilution series (1/4-1/64 pieces of a menstrual blood/vaginal swab) and, optionally, bona fide or mock casework samples of human or non-human origin were analyzed by 24 participating laboratories, using RNA extraction or RNA/DNA co-extraction methods. Two novel menstrual blood mRNA multiplexes were used: MMP triplex (MMP7, MMP10, MMP11) and MB triplex (MSX1, LEFTY2, SFRP4) in conjunction with a housekeeping gene triplex (B2M, UBC, UCE). Two novel mRNA multiplexes and a HBD1 singleplex were used for the identification of vaginal secretion: Vag triplex (MYOZ1, CYP2B7P1 and MUC4) and a Lactobacillus-specific Lacto triplex (Ljen, Lcris, Lgas). The laboratories used different chemistries and instrumentation and all were able to successfully isolate and detect mRNA in dried stains. The simultaneous extraction of RNA and DNA allowed for positive identification of the tissue/fluid source of origin by mRNA profiling as well as a simultaneous identification of the body fluid donor by STR profiling, also from old and compromised casework samples. The results of this and the previous collaborative RNA exercises support RNA profiling as a reliable body fluid identification method that can easily be combined with current STR typing technology.

A single-nucleotide polymorphism (SNP) multiplex system: the association of five SNPs with human eye and hair color in the Slovenian population and comparison using a Bayesian network and logistic regression model

AIM

To analyze two phenotype characteristics--eye and hair color--using single-nucleotide polymorphisms (SNPs) and evaluate their prediction accuracy in Slovenian population.

METHODS

Twelve SNPs (OCA2 - rs1667394, rs7170989, rs1800407, rs7495174; HERC2 - rs1129038, rs12913832; MC1R - rs1805005, rs1805008; TYR - rs1393350; SLC45A2 - rs16891982, rs26722; SLC24A5 - rs1426654) were used for the development of a single multiplex assay. The single multiplex assay was based on SNaPshot chemistry and capillary electrophoresis. In order to evaluate the accuracy of the prediction of eye and hair color, we used the logistic regression model and the Bayesian network model, and compared the parameters of both.

RESULTS

The new single multiplex assay displayed high levels of genotyping sensitivity with complete profiles generated from as little as 62 pg of DNA. Based on a prior evaluation of all SNPs in a single multiplex, we focused on the five most statistically significant in our population in order to investigate the predictive value. The two prediction models performed reliably without prior ancestry information, and revealed very good accuracy for both eye and hair color. Both models determined the highest predictive value for rs12913832 (P<0.0001), while the other four SNPs (rs1393350, rs1800407, rs1805008, and rs7495174) showed additional association for color prediction.

CONCLUSION

We developed a sensitive and reliable single multiplex genotyping assay. More samples from different populations should be analyzed before this assay could be used as one of the supplemental tools in tracing unknown individuals in more complicated crime investigations.

Forensic DNA databases in Western Balkan region: retrospectives, perspectives, and initiatives

The European Network of Forensic Science Institutes (ENFSI) recommended the establishment of forensic DNA databases and specific implementation and management legislations for all EU/ENFSI members. Therefore, forensic institutions from Bosnia and Herzegovina, Serbia, Montenegro, and Macedonia launched a wide set of activities to support these recommendations. To assess the current state, a regional expert team completed detailed screening and investigation of the existing forensic DNA data repositories and associated legislation in these countries. The scope also included relevant concurrent projects and a wide spectrum of different activities in relation to forensics DNA use. The state of forensic DNA analysis was also determined in the neighboring Slovenia and Croatia, which already have functional national DNA databases. There is a need for a ‘regional supplement’ to the current documentation and standards pertaining to forensic application of DNA databases, which should include regional-specific preliminary aims and recommendations.

Identification of skeletal remains of Communist Armed Forces victims during and after World War II: combined Y-chromosome (STR) and MiniSTR approach

AIM

To report on the use of STR, Y-STRs, and miniSTRs typing methods in the identification of victims of revolutionary violence and crimes against humanity committed by the Communist Armed Forces during and after World War II in which bodies were exhumed from mass and individual graves in Slovenia.

METHODS

Bone fragments and teeth were removed from human remains found in several small and closely located hidden mass graves in the Skofja Loka area (Lovrenska Grapa and Zolsce) and 2 individual graves in the Ljubljana area (Podlipoglav), Slovenia. DNA was isolated using the Qiagen DNA extraction procedure optimized for bone and teeth. Some DNA extracts required additional purification, such as N-buthanol treatment. The QuantifilerTM Human DNA Quantification Kit was used for DNA quantification. Initially, PowerPlex 16 kit was used to simultaneously analyze 15 short tandem repeat (STR) loci. The PowerPlex S5 miniSTR kit and AmpF/STR MiniFiler PCR Amplification Kit was used for additional analysis if preliminary analysis yielded weak partial or no profiles at all. In 2 cases, when the PowerPlex 16 profiles indicated possible relatedness of the remains with reference samples, but there were insufficient probabilities to call the match to possible male paternal relatives, we resorted to an additional analysis of Y-STR markers. PowerPlex Y System was used to simultaneously amplify 12 Y-STR loci. Fragment analysis was performed on an ABI PRISM 310 genetic analyzer. Matching probabilities were estimated using the DNA-View software.

RESULTS

Following the Y-STR analysis, 1 of the "weak matches" previously obtained based on autosomal loci, was confirmed while the other 1 was not. Combined standard STR and miniSTR approach applied to bone samples from 2 individual graves resulted in positive identifications. Finally, using the same approach on 11 bone samples from hidden mass grave Zolosce, we were able to obtain 6 useful DNA profiles.

CONCLUSION

The results of this study, in combination with previously obtained results, demonstrate that Y-chromosome testing and mini-STR methodology can contribute to the identification of human remains of victims of revolutionary violence from World War II.

Validation of SRY marker for forensic casework analysis

Determining the gender of the source of forensic DNA evidence is based on the amelogenin test. However, at times the assay may not be indicative of gender assignment, because of deletions at the amelogenin site. Previously, we described successful coamplification of a marker residing within the SRY gene with the short tandem repeat markers from two commercially available human identification kits. The study herein addresses the validation of primers for the target SRY gene regarding specificity, sensitivity, and robustness. Among 115 unrelated male Slovenians no null allele was observed. Repeatable and reliable results were obtained from as little as 25 pg of template DNA, indicating a high sensitivity of detection for the assay. No polymerase chain reaction product was observed even at a concentration of 10 ng/microL of template female DNA. Additionally, the male specific marker could be detected in mixed male and female samples down to a ratio of 1:16.

DNA identification of skeletal remains from the World War II mass graves uncovered in Slovenia

AIM

To present the joint effort of three institutions in the identification of human remains from the World War II found in two mass graves in the area of Skofja Loka, Slovenia.

METHODS

The remains of 27 individuals were found in two small and closely located mass graves. The DNA was isolated from bone and teeth samples using either standard phenol/chloroform alcohol extraction or optimized Qiagen DNA extraction procedure. Some recovered samples required the employment of additional DNA purification methods, such as N-buthanol treatment. Quantifiler Human DNA Quantification Kit was used for DNA quantification. PowerPlex 16 kit was used to simultaneously amplify 15 short tandem repeat (STR) loci. Matching probabilities were estimated using the DNA View program.

RESULTS

Out of all processed samples, 15 remains were fully profiled at all 15 STR loci. The other 12 profiles were partial. The least successful profile included 13 loci. Also, 69 referent samples (buccal swabs) from potential living relatives were collected and profiled. Comparison of victims' profile against referent samples database resulted in 4 strong matches. In addition, 5 other profiles were matched to certain referent samples with lower probability.

CONCLUSION

Our results show that more than 6 decades after the end of the World War II, DNA analysis may significantly contribute to the identification of the remains from that period. Additional analysis of Y-STRs and mitochondrial DNA (mtDNA) markers will be performed in the second phase of the identification project.

Population Data at Two Short Tandem Repeat Loci D2S1338 and D19S433 in the Sample of Multinational Bosnia and Herzegovina Residents

POPULATION

We have analyzed the distribution of allele frequencies at two short tandem repeats loci (D2S1338 and D19S433) in a multinational sample of Bosnia and Herzegovina (B&H) residents. A total of 110 unrelated male and female individuals (Caucasians) from different regions of B&H were sampled for the analysis. We ensured that the sample reflected approximate proportional participation of the three main ethnic groups in the population of B&H (Bosniacs-Muslim [45%], Serbs [34%], Croats [21%]).

Allele frequencies for 15 short tandem repeat loci in a representative sample of Bosnians and Herzegovinians

Allele frequencies for the 15 STR loci (D3S1358, TH01, D21S11, D18S51, Penta E, D5S818, D13S317, D7S820, D16S539, CSF1PO, Penta D, vWA, D8S1179, TPOX, FGA) included in the PowerPlex 16 kit were obtained from a multiethnic sample of 100 unrelated individuals born in Bosnia and Herzegovina.

The peopling of modern Bosnia-Herzegovina: Y-chromosome haplogroups in the three main ethnic groups

The variation at 28 Y-chromosome biallelic markers was analysed in 256 males (90 Croats, 81 Serbs and 85 Bosniacs) from Bosnia-Herzegovina. An important shared feature between the three ethnic groups is the high frequency of the "Palaeolithic" European-specific haplogroup (Hg) I, a likely signature of a Balkan population re-expansion after the Last Glacial Maximum. This haplogroup is almost completely represented by the sub-haplogroup I-P37 whose frequency is, however, higher in the Croats (approximately 71%) than in Bosniacs (approximately 44%) and Serbs (approximately 31%). Other rather frequent haplogroups are E (approximately 15%) and J (approximately 7%), which are considered to have arrived from the Middle East in Neolithic and post-Neolithic times, and R-M17 (approximately 14%), which probably marked several arrivals, at different times, from eastern Eurasia. Hg E, almost exclusively represented by its subclade E-M78, is more common in the Serbs (approximately 20%) than in Bosniacs (approximately 13%) and Croats (approximately 9%), and Hg J, observed in only one Croat, encompasses approximately 9% of the Serbs and approximately 12% of the Bosniacs, where it shows its highest diversification. By contrast, Hg R-M17 displays similar frequencies in all three groups. On the whole, the three main groups of Bosnia-Herzegovina, in spite of some quantitative differences, share a large fraction of the same ancient gene pool distinctive for the Balkan area.

Comparative Study of Genetic Variation at 15 STR Loci in Three Isolated Populations of the Bosnian Mountain Area

Fifteen autosomal STR loci (D3S1358, TH01, D21S11, D18S51, Penta E, D5S818, D13S317, D7S820, D16S539, CSF1PO, Penta D, VWA, D8S1179, TPOX, and FGA) were studied in three geographically close but isolated populations from the Bosnian mountain area. The three villages are Bobovica, Dejcici, and Lukomir. DNA was obtained from 83 individuals, and the allele frequencies and genetic diversity among the three sample groups were compared. In addition, seven of the STR loci (CSF1PO, D13S317, D3S1358, D5S818, D7S820, FGA, TH01) were used in a comparative population analysis of the Bjelasnica-Treskavica region and the Adriatic islands of Brac, Hvar, and Korcula. Although the sample sizes are relatively small, the observed variation within any of the small isolated populations is high and comparable to less isolated groups. In addition, even though the populations are geographically isolated, the STR data are similar among the populations. The most significant frequency differences were observed at the TH01 locus. Although the specific allele distributions in any untyped population cannot be determined a priori, we find support for a high degree of diversity for the STR loci in most populations. In addition, the multiple locus profile is highly informative not only for various population studies but also for forensic studies, even when specific population data are not available.

Mitochondrial DNA variability in Bosnians and Slovenians

Mitochondrial DNA variability in two Slavonic-speaking populations of the northwestern Balkan peninsula, Bosnians (N = 144) and Slovenians (N = 104), was studied by hypervariable segments I and II (HVS I and II) sequencing and restriction fragment-length polymorphism (RFLP) analysis of the mtDNA coding region. The majority of the mtDNA detected in Southern Slavonic populations falls into the common West Eurasian mitochondrial haplogroups (e.g., H, pre-V, J, T, U, K, I, W, and X). About 2% of the Bosnian mtDNAs encompass East Eurasian and African lineages (e.g., M and L1b, respectively). The distribution of mtDNA subclusters in Bosnians, Slovenians and the neighbouring European populations reveals that the common genetic substratum characteristic for Central and Eastern European populations (such as Germans, Poles, Russians and Finns) penetrates also South European territories as far as the Western Balkans. However, the observed differentiation between Bosnian and Slovenian mtDNAs suggests that at least two different migration waves of the Slavs may have reached the Balkans in the early Middle Ages.

Analysis of DNA evidence recovered from epithelial cells in penile swabs

In the rape case presented here, no semen, hair, or other biological evidence were left by the perpetuator at the crime scene or on the victim. The alleged assailant was arrested soon after the crime. A classical stain recovery technique using cotton swab moistened with sterile water was taken for recovering potential female epithelial cells and leukocytes deposited on the alleged assailant's penis during sexual assault. The organic method used for DNA extraction was quantified according to the slot-blot procedure and amplified at 9 and 15 polymorphic loci. Penile swab revealed a DNA profile of mixed origin. In addition to the suspect's DNA profile, DNA contribution from the victim was identified as a minor component in the mixture. Frequency of the profile resulted in a value of 5 x 10(-14) for the multiplex systems AmpFlSTR Plus and 2.5 x 10(-18) for the multiplex system PowerPlex 16, taking into account only non-overlapping alleles between the suspect and the victim from the minor component in the DNA mixture. Moreover, three additional alleles were observed at D21S11 locus by use of PowerPlex and STR SGM plus primers, which could not belong to the suspect. The victim's DNA profile showed the same three-banded genotype at this locus. The same pattern was detected when the victim's saliva or blood were used as reference samples. Our laboratory finding was consistent with the police report that the victim was a person with Down syndrome, a human genetic disease mainly resulting from trisomy (triplication) of the 21 chromosome.

STR data for the AmpFlSTR SGM plus from Slovenia

Allele frequencies for the 10 STRs loci included in the AmpFISTR SGM Plus kit (PE Applied Biosystem) were obtained from a sample of 321 unrelated individuals born in Slovenia.

The analysis of three short tandem repeat (STR) loci in the Slovene population by multiplex PCR

Allele frequencies for three tetrameric short tandem repeat (STR) loci D3S1358, HUMVWA, and HUMFGA were determined in a Slovene Caucasian population sample. DNA samples from a total of 221 Slovenes were amplified by multiplex PCR using the commercial kit AmpFISTR Blue (Perkin-Elmer). Separation and detection of the amplified STR fragments were carried out using a 377 automated genetic analyzer (Applied Biosystem Division/Perkin Elmer). Seven alleles at the D3S1358 locus, 8 alleles at the HUMVWA31A locus, and 13 alleles at the HUMFGA locus were observed. A deviation from Hardy-Weinberg equilibrium was observed, only at the HUMVWA31A locus (p = 0.045, exact test). The departure at this locus was not significant after Bonferroni correction. There were no detectable departures between pairwise comparisons of the loci. The combined power of discrimination for all three loci is 0.9998, and the power of exclusion is 0.9526. The observed allele frequencies for the loci D3S1358, HUMVWA31A, and HUMFGA are similar to those in European and U.S. Caucasian populations.

The Slovenian population data on the PCR based loci HLA-DQA1, LDLR, GYPA, HBGG, D7S8, GC, and D1S80

Allele frequencies for the loci HLA-DQA1, LDLR, GYPA, HBGG, D7S8, GC, and D1S80 were determined for a sample population of unrelated individuals from Slovenia. All loci meet Hardy-Weinberg expectations, except the loci GYPA (p = 0.041) and D1S80 (p = 0.009). There is little evidence for association of alleles among the seven loci. Only one out of 21 pairwise comparisons demonstrated departures from independence (HLA-DQA1/HBGG, p = 0.008). The allelic frequency data generally are similar to that of U.S. Caucasians.

Improved purification of steroid 1:2-dehydrogenase from Nocardia opaca and partial characterization of its cloned gene sequence

We have purified a steroid-inducible 1:2-dehydrogenase from Nocardia opaca. The final enzyme preparation was purified 120-fold with a recovery of 38%. The N-terminal amino acid sequence was determined to be: Met-Gln-Asp-Trp-Thr-Ser-Glu-(Cys)-Asp-Val-Leu-Val-Val-Gly-. From the genomic library of Nocardia opaca in the plasmid pUC19, a clone designated as pSTD23 containing a 0.9 kb KpnI-PstI fragment was found to hybridize with an oligonucleotide probe corresponding to the first six amino acids from the N-terminal of the purified protein. The nucleotide sequence of the upstream region and a part of the structural domain were determined. The sequence of the first 56 amino acids of the steroid 1:2-dehydrogenase from Nocardia opaca as deduced from its gene sequence showed a 58% homology with the corresponding gene from Pseudomonas testosteroni, and the conservative sequences in the FAD-binding domain were also determined.