Development and Validation of a GBAÔ-based mtDNA Sequence Evaluation System for Forensic Testing
John V. Planz1, Vicki M. Fish1, John M. Rader1,
Matthew W. DuPont1, Joseph E. Warren2, Peter Gill3, and
Robert C. Giles1. 1GeneScreen, Inc., 2600 Stemmons Fwy, Suite 133,
Dallas, Texas 75207
2Tarrant County Medical Examiners Office, Fort Worth, Texas
3Forensic Science Service, Priory House, Gooch Street North, Birmingham, United
Kingdom.
× Ø × Ø × Ø × Ø × Ø × Ø × Ø × Ø × Ø × Ø × Ø × Ø × Ø × Ø × Ø
ABSTRACT
In recent years, emphasis has increased for the forensic typing of poorly preserved and limited samples. Routine mitochondrial DNA (mtDNA) typing of single hairs, skeletal remains and extremely decomposed materials has been greatly facilitated by the advancements of automated DNA sequencing following PCR amplification. Although this methodology provides very precise information regarding the DNA sequence of the two hypervariable regions of the mtDNA control region, the process itself is extremely time and labor consumptive considering the low number of informative loci that are actually encountered. In an attempt to streamline testing of mtDNA sequences, specific testing of the more highly variable base positions within the control region have been developed, which include an SSO-based test strip method and a microsequencing assay. Using the flexible platform of Genetic Bit Analysis (GBA
Ô) we have developed a rapid screening system for several variable mtDNA control region variants which circumvents the hybridization problems of the ASO system and equipment requirements of microsequencing, while providing an extremely robust and inexpensive test with excellent throughput potential. GBAÔ relies on dideoxy chain terminator chemistry followed by an ELISA-like colorimetric detection system. The initial development of the mtDNA array discussed here consists of ten specific base positions, mirroring the markers developed by the British Forensic Science Service for microsequencing, reporting a discriminating power (Pm) of 0.063 for Caucasians and 0.032 for Afro-Caribbeans. Additional markers, currently under development, significantly increase the value of this test on North American populations. Results of GBAÔ testing were verified through DNA sequencing, and environmental insult samples were tested to establish the robustness of the testing system. The GBAÔ system independently assays for all four of the possible nucleotides at each of the variable positions. As a result of this testing methodology, heteroplasmic sites are clearly and unambiguously identified, allowing the full value of these variants to be realized. The GBAÔ -based mtDNA testing platform allows for the rapid and inexpensive screening of mtDNA markers for routine forensic casework and has the added benefit of high throughput for sample analysis in mass disaster circumstances.INTRODUCTION
The analysis of intraspecific mitochondrial DNA (mtDNA) variation has seen increasing attention in recent years. Initially, studies of both coding and non-coding mtDNA regions were used to develop hypotheses concerning human evolution and the dispersal of human populations throughout the world. With these studies came the refinement of analytical methodologies, originating from RFLP restriction site analysis and leading ultimately to full DNA sequencing of regions under examination. The development of the PCR process and its rapid incorporation into DNA sequence analysis facilitated mtDNA analyses on less than pristine materials such as mummified remains and other ancient materials. As these studies became more commonplace and routine, application of the techniques began to interest the forensic community. Forensic uses of mtDNA methods erupted initially in high profile studies, such as the identification of the Romanov family remains in Russia and the search for Butch Cassidy and the Sundance Kid.
As the desire to incorporate mtDNA analyses into forensic case work increased, methods to increase the ease and throughput of mtDNA typing were developed. Most notable of these was the development of an SSO probe hybridization assay which interrogated nine regions of the mtDNA control region with 23 probes (1). The application of this method was extensively used to characterize the extent of population heterogeneity expressed in Asians, Europeans and sub-Saharan Africans (2-4). Gel-based minisequencing assays were also developed to expedite mtDNA sequence analysis for crime laboratory and databasing operations (5,6). The minisequencing format interrogates nucleotide substitutions at 13 polymorphic positions in the HV1 and HV2 regions of the mtDNA control region through the use of migration modifying sequencing primers. Random match probabilities from this system were reported as 0.054 and 0.026 for British Caucasians and Afro-Caribbeans, respectively. Comparable results were reported for these populations with the SSO assay (1) and a reverse dot-blot typing kit adapted from the SSO system (7).
This work describes the design and implementation of a new testing platform, GBA
Ô (8), for determining mtDNA haplotypes for forensic investigations. GBAÔ, named by analogy to the storage of computer data in bits, is a non-radioactive single-base sequencing method performed on a solid-phase substrate. The GBAÔ system involves the following steps; 1) PCR amplification of the target DNA sequence containing the variant base using one exonuclease-resistant primer per set of primer pairs; 2) Generation of a single stranded target by exonuclease digestion; 3) Capture of the target strand by hybridization to a pre-synthesized GBAÔ primer covalently attached to the substrate; 4) Extension of the GBAÔ primer with labeled chain terminators through single base extension; and 5) Detection of the extension product by colorimetric or direct fluorescent assay.METHODS
Organic DNA extractions were performed on 513 blood samples from Caucasians (n=103), African Americans (n=109), Hispanics ("eastern" south Florida, n=116; "western" south Texas, n=108), Asians (n=47), and Native Americans (n=31). Process groups of 78 individual DNA samples were prepared in 96-well plates to facilitate automated sample handling during PCR setup and GBA
Ô hybridization and detection.Ten mtDNA loci were developed for GBA
Ô testing (Figure 1) that were amplified in two amplicons consisting of the HV1 and HV2 regions (Figure 1). Amplification primers corresponding to mtDNA heavy strand sequences were prepared with the first four 5 nucleotides linked with phosphorothioated bonds. Amplifications were conducted with PE Taq DNA polymerase under typical buffer conditions and 30 cycles. Following amplification, the amplified product was digested with T7 exonuclease (5' to 3' exonuclease) added directly to the PCR reaction wells by an ICN multidrop. The enzyme digests the unprotected (non-phosphorothioated) DNA strand, yielding >99% single stranded product after 1 hour incubation at room temperature. Single stranded template DNA is then hybridized to 96-well plate bound GBAÔ capture oligonucleotides specific to the locus being interrogated. GBAÔ plates are divided into separate plates by locus each having space for 78 samples and controls. GBAÔ hybridization salts and an aliquot of the amplified DNA samples from the PCR plates was dispensed into each reaction well by a Biomek fluid handler. The GBAÔ plates also contain wells carrying control templates, reagent controls and negative controls.Hybridization of the single-stranded PCR DNA to the capture oligonucleotide was carried out by adding an equal volume of 3M NaCl/50 mM EDTA (10ul) to the PCR reactions, and incubating the template in the presence of the plate-bound GBA
Ô primer for 15 minutes at room temperature. The plate was then washed three times with TNTW buffer (10mM Tris-HCl pH 7.5, 150 mM NaCl, 0.05% Tween® 20) using a standard microtiter plate washer (Dynatech), to remove unhybridized templates.Locus interrogation was conducted with an extension mix containing E. coli DNA polymerase, buffer and biotinylated and fluorosceinated acyclo chain terminators was added to the washed GBA
Ô plates by the Biomek. Although the majority of SNPs so far developed have displayed biallelic representation, for forensic test development each possible nucleotide is interrogated, thus two copies of each GBAÔ plate are assembled allowing for two color detection of all four bases in two parallel reactions for each SNP being evaluated.After the optimal extension time (1-5 minutes) at room temperature, the plates were washed with TNTW on an ICN stacking microtiter plate washer. The PCR template was then removed by denaturation following addition of 50ul of 0.2 N NaOH, and washing three times with TNTW buffer using ICN dispensing and washing robots. Incorporated fluorosceinated and biotinylated nucleotides were detected by independent enzyme-linked assays. Each well was incubated with 30ul of anti-fluoroscein-alkaline phosphatase for 30 minutes for the detection of fluorosceinated nucleotides and incorporated biotinylated nucleotides were detected by incubation with 20ul of anti-biotin-HRP for 10min. After washing 6 times with TNTW buffer, 100 ul of the substrate p-nitrophenyl phosphate(PNPP) at a final concentration of 1.5 mg/ml in 1 M diethanolamine, pH10 or 100ul of o-phenylenediamine (OPD, 1 mg/ml in 0.1 M Citric acid, pH4.5) containing 0.012% hydrogen peroxide was added to each well and allowed to incubate for 30 minutes.
Raw data (in OD units) were collected on a standard microplate reader at 405nm for fluorosceinated nucleotides and at 620nm for biotinylated nucleotides. The data were analyzed by a computer algorithm which compares the individual data points to plate bound calibration controls. The locus data were then assembled into haplotypes of each individual tested. Haplotype frequencies, population statistics and haplotype distributions were examined using Microsoft Excel. Estimates of genetic diversity were calculated for each population using
h =(1-
Sx2)n/(n-1),where
Sx2 is the sum of squares of haplotype frequencies and n is the sample size (9). Estimates of the probability of a random match were calculated simply as p=Sx2 following formula 2 of Stoneking et al. (1). Individuals exhibiting heteroplasmy at any of the ten loci were omitted from statistical calculations and are summarized independently.RESULTS
Five hundred and thirteen individuals were tested for the ten GBA
Ô loci developed in this study. Sixty-seven hapoltypes were detected among the 504 individuals that did not exhibit heteroplasmy at any locus. The frequency of the most common haplotype across all populations was approximately 20 percent. The ten-locus panel at present is capable of detecting approximately 90 percent of the genetic diversity in the regions being assessed. Statistical data, delineated by population, are found in Table 1. The results obtained with the ten locus panel described here are concordant with the data obtained by the minsequencing array of Tully et al., (6) and the overall distribution of haplotypes (Figure 2) is in agreement with population data presented by Melton et al. (3, 4) although our data represent fewer genetic loci. The inclusion of additional loci to the GBAÔ panel will increase the power of discrimination of the test to the 0.99 level desired.Heteroplasmy was detected in 16 individuals of the 513 tested and was limited to three loci: H00146 (7/513), H00154 (7/513), and H00075 (2/513). The GBA platform allowed for very clear and unambiguous detection of the occurrence of heteroplasmy in the database samples analyzed. In a preliminary investigation into the occurrence of individual heteroplasmy, a set of hairs collected from one individual were individually tested, both root and shaft, for the presence of multiple mtDNA types. The same mtDNA haplotype was detected in all replicates during this test. An expanded test of this nature will be part of the forensic validation studies, in addition to contamination and mixture tests, that are will be conducted when the final panel of loci has been chosen.
Forensic validation of the GBA
Ô-based kit is scheduled for completion by February of 1998. A panel consisting of 12 mtDNA loci will be tested on simulated casework samples and proficiency test materials following the guidelines of TWGDAM and DAB. Following validation studies, the GBAÔ testing platform will be organized into a kit that will be distributed to alpha and beta testing sites in the first quarter of 1998.
REFERENCES
1. Stoneking, M., Hedgecock, D., Higuchi, R.G., Vigilant, L., and Erlich, H. Population variation of human mtDNA control region sequences detected b enzymatic amplification and sequence specific oligonucleotide probes. Am. J. Hum. Genet. 1991; 48:370-382.
2. Melton, T. and Stoneking, M. Extent of heterogeneity in mitochondrial DNA of ethnic Asian populations. J. Forensic Sci. 1996; 41:591-602.
3. Melton, T., Wilson, M., Batzer, M., and Stoneking, M.. Extent of heterogeneity in
mitochondrial DNA of European populations.
J. Forensic Sci. 1997; 42:437-446.
4. Melton, T., Ginther, C., Sensabaugh, G., Soodyall, H., and Stoneking, M. Extent of
heterogeneity in mitochondrial DNA of Sub-Saharan African populations. J. Forensic Sci.
1997; 42:582-592.
5. Sullivan, K., Tully, G., Alliston-Greiner, R, Hopwood, A., Bark, J.E., and Gill, P. A two stage strategy for the automated analysis of mitochondrial DNA. In: Proceedings of the International Society of Forensic Haemogenetics; 1995; pp. 11-13.
6. Tully, G., Sullivan, K.M., Nixon, P., Stones, R.E., and Gill, P. Rapid detection of mitochondrial sequence polymorphisms using multiplex solid-phase fluorescent minisequencing. Genomics 1996; 34:107-113.
7. Reynolds, R., Clark, B., Marschak, T., Valaro, J., and Rubin, S. Mitochondrial DNA typing using sequence-specific oligonucleotide probes. In: Proceedings of the International Society of Forensic Haemogenetics; 1996; pp. 17-23.
8. Nikiforov, T.T., Rendle, R. B., Goelet, P., Rogers, Y., Kotewicz, M.L., Anderson, S. Trainor, G.L., and Knapp, M.R. Genetic Bit Analysis
Ô: a solid phase method for typing single nucleotide polymorphisms. Nuc. Acids. Res. 1994, 22:4167-4175.9. Tajima, F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 1989; 123:585-595.
Figure 1. mDNA loci in the HV1 and HV2 regions assayed by Genetic Bit Analysis

Figure 2. Distribution of mtDNA haplotypes in six populations using a ten-locus GBA panel.

Table 1. Population statistics for ten mtDNA loci interrogated through GBA
Ô methodology.| Population |
n |
# of types |
Frequency of most common type |
Genetic Diversity |
Probability of Identity within populations |
Probability of Identity between populations |
Power of Discrimination |
Am. Indian |
31 |
9 |
48.39 |
.725 |
29.87 |
5.63 |
0.701 |
Asain |
47 |
13 |
27.66 |
.873 |
14.53 |
2.64 |
0.854 |
Caucasian |
102 |
28 |
23.0 |
.907 |
10.15 |
1.41 |
0.899 |
Black |
106 |
28 |
22.11 |
.906 |
10.24 |
2.33 |
0.898 |
W. Hispanic |
102 |
20 |
30.69 |
.840 |
16.04 |
2.43 |
0.840 |
E. Hispanic |
116 |
31 |
19.83 |
.918 |
8.25 |
1.31 |
0.920 |
Total |
504 |
67 |
20.44 |
.922 |
7.79 |
- |
0.922 |
Go to proceedings home page