Mitochondrial DNA typing screens with control region and coding region SNPs

Mitochondrial DNA (mtDNA) analysis has found an important niche in forensic DNA typing. It is used with highly degraded samples or low-copy number materials such as might be found from shed hair or bones exposed to severe environmental conditions. The primary advantage of mtDNA is that it is present in high copy number within cells and therefore more likely to be recovered from highly degraded specimens. A major disadvantage to traditional forensic mtDNA analysis is that it is time-consuming and labor-intensive to generate and review the 610 nucleotides of sequence information commonly targeted in hypervariable regions I and II (HVI and HVII) of the control region. In addition, common haplotypes exist in HVI/HVII mtDNA sequences that can reduce the ability to differentiate two unrelated samples. In this report we describe the utility of two newly available screening assays for rapid exclusion of non-matching samples. The LINEAR ARRAY mtDNA HVI/HVII Region-Sequencing Typing Kit (Roche Applied Science, Indianapolis, IN) was used to type 666 individuals from U.S. Caucasian, African American, and Hispanic groups. Processing of the LINEAR ARRAY probe panels "mito strips" was automated on a ProfiBlot workstation. Observable variation in 666 individuals is reported and frequencies of the mitotypes within and between populations are presented. Samples exhibiting the most common Caucasian mitotype were subdivided with a multiplexed amplification and detection assay using eleven single nucleotide polymorphisms in the mitochondrial genome. These types of screening assays should enable more rapid evaluation of forensic casework samples such that only samples not excluded would be subjected to further characterization through full HVI/HVII mtDNA sequence analysis.     

Typing the 1.1 kb control region of human mitochondrial DNA in Japanese individuals

This study presents a reliable method that uses high-fidelity long-range PCR and optimized primers to assess polymorphism and to genotype human mitochondrial DNA (mtDNA). This method was used to analyze polymorphic sites in the human mtDNA control region, including hypervariable regions I, II, and III (HVI, HVII, and HVIII), from 124 unrelated Japanese individuals. In HVI, HVII, and HVIII, 80, 37, and 14 polymorphic sites were identified, respectively, excluding those in the homopolymeric cytosine stretch (C-stretch) regions. The region between HVI and HVII also contained 15 polymorphic sites. On the other hand, C-stretch length heteroplasmy in HVI or HVII was observed in 66 of 124 Japanese individuals (53%), which is much higher than in Caucasian populations. The variants in the C-stretch regions were characterized by counting the number of heteroplasmic peaks split from the single peak in homoplasmic sequences (i.e., 16244G and 16255G in HVI and 285G in HVII). Including the C-stretch length heteroplasmy, the 124 Japanese mtDNA samples were classified into 116 distinct haplotypes. The random match probability and the genetic diversity were estimated to be 0.95% and 0.998581, respectively, indicating that the method presented here has higher discrimination than the conventional method for mtDNA typing using HVI and HVII. [Correction added after publication 30 January 2007: in the preceding sentence random match probability and genetic diversity estimates were corrected from 0.95 and 0.998581%, respectively, to 0.95% and 0.998581, respectively.] The haplogroups and their frequencies observed in this study (i.e., D4; 13.7%, M7a1; 11.3%, D4a; 9.7% and M7b2; 8.9%) were similar to those observed in other studies of Japanese mtDNA polymorphism. The method described here is suitable for forensic applications, as shown by successful analysis of tissues from highly putrefied remains of an infant, which allowed maternal relationship to be determined via mtDNA haplotyping.     

Is it possible to differentiate mtDNA by means of HVIII in samples that cannot be distinguished by sequencing the HVI and HVII regions?

The mitochondrial control region includes three so-called hypervariable (HV) regions, in which the polymorphic positions show a particularly high frequency. According to a population study of 200 unrelated individuals from Germany, HVI (positions 16,024-16,365, according to Anderson) showed 88 variable positions in a total length of 342 bp (26%) and HVII (positions 73-340) displayed 65 mutable sites in 268 bp (24%). HVIII (positions 438-574) exhibited a slightly lower variability, with 25 polymorphic sites within 137 bp (18%), but contrasted clearly with the background, which showed variability rates of only 7% (positions 16,366-16,569, 1-72) and 3% (positions 341-437), respectively. At present, the displacement (D)-loop database in Magdeburg comprises 904 sequences of the mitochondrial HVI region and HVII region from Germans, Austrians and Swiss. By means of this material, the extent to which the mtDNA sequences that do not differ in the HVI and HVII regions can be differentiated by additionally sequencing HVIII was investigated.     

Detection of sequence variation in the HVII region of the human mitochondrial genome in 689 individuals using immobilized sequence-specific oligonucleotide probes

We have developed a rapid, immobilized probe-based assay for the detection of sequence variation in the hypervariable segment II (HVII) of the mitochondrial DNA (mtDNA) control region. Using a panel of 17 sequence-specific oligonucleotide (SSO) probes immobilized on nylon membrane strips, we typed 689 individuals from four population groups. The genetic diversity value for each population was calculated from the frequency data, and the frequencies of distinct "mitotypes" in each group were determined. We performed DNA sequence analysis of 129 samples to characterize the sequences associated with "blanks" (absence of probe signals) and weak probe signals. Out of 689 samples, we observed five heteroplasmic samples (excluding the variable C-stretch beginning at position 303) using the immobilized SSO probe panel. The SSO probe strips were used for the analysis of shed hairs and bloodstains from several criminal cases in Sweden, one of which is described here. We conclude that this mtDNA typing system is useful for human identification and significantly decreases casework turnaround time.     

Different informativeness of the three hypervariable mitochondrial DNA regions in the population of Bologna (Italy)

Mitochondrial DNA (mtDNA) sequence variations at hypervariable regions HVI, HVII and HVIII were analysed in 100 unrelated Italians from Bologna. The values of the statistical parameters are in agreement with the range of European populations. We suggest that the less informative HVIII region may be useful to distinguish HVI-HVII identical sequences in forensic analysis especially when nuclear DNA cannot be investigated.     

Forensic casework analysis using the HVI/HVII mtDNA linear array assay

The mitochondrial hypervariable regions I and II have proven to be a useful target for analysis of forensic materials, in which the amount of DNA is limited or highly degraded. Conventional mitochondrial DNA (mtDNA) sequencing can be time-consuming and expensive, limitations that can be minimized using a faster and less expensive typing assay. We have evaluated the exclusion capacity of the linear array mtDNA HVI/HVII region-sequence typing assay (Roche Applied Science) in 16 forensic cases comprising 90 samples. Using the HVI/HVII mtDNA linear array, 56% of the samples were excluded and thus less than half of the samples require further sequencing due to a match or inconclusive results. Of all the samples that were excluded by sequence analysis, 79% could be excluded using the HVI/HVII linear array alone. Using the HVI/HVII mtDNA linear array assay, we demonstrate the potential to decrease sequencing efforts substantially and thereby reduce the cost and the turn-around time in casework analysis.     

Characterization of human control region sequences of the African American SWGDAM forensic mtDNA data set

The scientific working group on DNA analysis Methods (SWGDAM) mitochondrial DNA (mtDNA) population data set is used to infer the relative rarity of control region mtDNA profiles obtained from evidence samples and of profiles used for identification of missing persons. In this study, the African American haplogroup patterns in the SWGDAM data were analyzed in a phylogenetic context to determine relevant single nucleotide polymorphisms (SNPs) and to describe haplogroup distributions for Africans observed in these data sets. Over 200 SNPs (n=217) were observed in the African American data set (n=1148). These SNPs ranged from having 1-39 changes in the phylogenetic tree, with sites 152 and 16519 being the most variable. On average there were 5.8 changes for a character on the tree. The most variable sites (with 19 or more changes each) observed included 16093, 16129, 16189, 16311, 16362, 16519, 146, 150, 152, 189, and 195. These rapidly changing sites are consistent with other published analyses. Only 34 SNPs are needed to identify all clusters containing 10 or more individuals in the African American data set. The results show that the African American SWGDAM mtDNA data set contains variation consistent with that described in continental African populations. Thirteen of the 18 haplogroups previously observed in African populations were observed and include: L1a, L1b, L1c, L2a, L2b, L2c, L3b, L3d, L3e1, L3e2, L3e3, L3e4 and L3f. Haplogroup L2a is the most commonly observed cluster (18.8%) in the African American data set. The next most common haplogroups in the African American data set include the clusters L1c (11.0%), L1b (9.1%), L3e2 (9.0%) and L3b (8.1%). Approximately 8% of the haplogroups observed within African Americans were common in European Caucasians or East Asians; these were H (n=32), J (n=4), K (n=5), T (n=2), U5 (n=6), U6 (n=9 also known from North Africa), A (n=12), B (n=7), C (n=4), and M (n=16), respectively. The European Caucasian and East Asian haplogroups are expected due to admixture between individuals with recent ancestry in Western Eurasia and sub-Saharan Africa. The genetic characterization of these relevant data sets is fully consistent with other published mtDNA genetic variation. The sequence diversity observed in this data set makes it a valuable tool for forensic applications.     

Sequence polymorphism of mitochondrial D-loop DNA in the Taiwanese Han population

In order to demonstrate the sequence diversity of mitochondrial D-loop DNA in the Taiwanese Han population, we established a database of 155 unrelated individuals. For each individual, the complete 980bp DNA region from the 5' end of HVI to 3' end of HVII segment was sequenced. In these 155 sequence data, 149 different haplotypes were observed, amongst these haplotypes, 144 were unique, 4 were found in 2 individuals and 1 was found in 3 individuals. When compare to the Anderson sequence, 144 transitions, 24 transversions, 5 insertions and 5 deletions were found. Eight positions exhibited more than one polymorphic sequence, six exhibited two variants while two exhibited three variants. Over the 1024bp that was analysed, pairwise differences between the sequences were 11.35+/-3.53bp. The sequence and nucleotide diversity were 0.9994 and 0.0116, respectively. The probability of two individuals randomly matching over the entire control region was 0.007. The diversity in the mitochondrial D-loop indicates the value of this locus for casework within Taiwan.     

Analysis of eight mtDNA coding region polymorphisms for characterization of the female lineages ancestry in Alagoas,Brazil

In this work we evaluated eight mitochondrial DNA coding region polymorphisms (8281–8289d, 1736, 13263, 4883, 3594, 10873, 10400, and 12705) in order to identify lineages of African, Amerindian and European origin in Alagoas, Northeastern Brazil. Seven of these polymorphisms were detected by multiplex minisequencing with SnaPshot. For validation purposes a sample previously characterized by sequencing of the HVI and HVII regions was analyzed. Good agreement was observed between the phylogenetic results obtained in this work and the HVI and HVII regions sequencing. The results obtained suggest that the set of markers evaluated in this study can be used in ancestry studies of Brazilian populations.     

A rape case solved by mitochondrial DNA mixture analysis

In cases of stains that contain mixed DNA from different contributors, analyzing mitochondrial DNA (mtDNA) requires the use of cloning techniques. We developed an efficient cloning technique that was applied in a rape case. After a differential lysis-based DNA extraction from vaginal swabs, hypervariable region I and II (HVI, HVII) amplicons obtained from the male fraction were cloned. Although we mainly found the victim's haplotype, we were able to detect the suspect's haplotype in two clones for HVI and in one clone for HVII. As the midpiece of the flagellum, which contains mitochondria, can be lost during the differential lysis, we also investigated the female fraction by cloning to evaluate the proportion of victim/suspect mtDNA. Unfortunately, only clones presenting the victim's haplotype were found. This case highlights the need for an optimal differential lysis protocol to enrich the male fraction not only with nuclear but also mitochondrial DNA.     

SSCP法mtDNA分型的法医学应用

目的探讨SSCP法mtDNA分型在法医鉴定实践中的意义。方法PCR扩增mtDNAHV-I和HV-II的序列多态性片段,直接用SSCP技术分型。对70个湖北汉族真三联家系和140例随机个体进行检测,比较家系中mtDNA的单倍型SSCP图谱。统计分析SSCP法用于两个高变区在母系确定、个人识别中的意义和鉴别能力。结果在70个家系中,母亲和孩子的HV-I、HV-IISSCP带型完全相同;家系中父亲与母子的HV-I图谱不同占98.57%;HV-IISSCP图谱不同占97.13%。140例随机个体的HVI、HVII区分别检出21、16种单倍型,GD值分别为0.9556、0.9356。结论mtDNA-SSCP分型在嫌疑人筛查和母系亲缘关系推断中有实际应用价值。     

Variability of the Entire Mitochondrial DNA Control Region in a Human Isolate from the Pas Valley (Northern Spain)

Abstract: In this study, we analyzed the entire mtDNA control region in 61 unrelated individuals from the Pas Valley (Cantabria), a human isolate from northern Spain, to evaluate the suitability of this analysis to increase the power of discrimination of this locus for forensic purposes in human isolates. Low values obtained for the diversity parameters confirmed the relative isolation of this human group. The main findings of this study indicated that even the analysis of the entire mtDNA control region may have important limitations for use in forensic casework when dealing with human isolates: none of the 44 individuals who exhibited identical HVI‐HVII haplotypes could be further differentiated by analysis of segment HVIII. Nevertheless, analysis of the entire mtDNA control region proved to be useful to determine the ancestry of the samples examined, by contributing to the confirmation, and, on occasion, even to the refinement of the haplogroup assignment.     

Twelve short tandem repeat loci Y chromosome haplotypes: genetic analysis on populations residing in North America

A total of 2443 male individuals, previously typed for the 13 CODIS STR loci, distributed across the five North American population groups African American, Asian, Caucasian, Hispanic, and Native American were typed for the Y-STR loci DYS19, DYS385a/b, DYS389I/II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438 and DYS439 using the PowerPlex Y System. All population samples were highly polymorphic for the 12 Y-STR loci with the marker DYS385a/b being the most polymorphic across all sample populations. The Native American population groups demonstrated the lowest genetic diversity, most notably at the DYS393 and DYS437 loci. Almost all of the 12-locus haplotypes observed in the sample populations were represented only once in the database. Haplotype diversities were greater than 99.6% for the African Americans, Caucasians, Hispanics, and Asians. The Native Americans had the lowest haplotype diversities (Apaches, 97.0%; Navajo, 98.1%). Population substructure effects were greater for Y-haplotypes, compared with that for the autosomal loci. For the apportionment of variance for the 12 Y-STRs, the within sample population variation was the largest component (>98% for each major population group and approximately 97% in Native Americans), and the variance component contributed by the major population groups was less than the individual component, but much greater than among sample populations within a major group (11.79% versus 1.02% for African Americans/Caucasians/Hispanics and 15.35% versus 1.25% for all five major populations). When each major population is analyzed individually, the R(ST) values were low but showed significant among group heterogeneity. In 692 confirmed father-son pairs, 14 mutation events were observed with the average rate of 1.57x10(-3)/locus/generation (a 95% confidence bound of 0.83x10(-3) to 2.69x10(-3)). Since the Y-STR loci reside on the non-recombining region of the Y chromosome, the counting method is one approach suggested for conveying an estimate of the rarity of the Y-haplotype. Because the Y-STR loci are not all in disequilibrium to the same extent, the counting method is a very conservative approach. The data also support that autosomal STR frequencies can be multiplied by the upper bound frequency estimate of a Y-haplotype in the individual population group or those pooled into major population groups (i.e., Caucasian, African American, Hispanic, and Asian). These analyses support use of the haplotype population data for estimating Y-STR profile frequencies for populations residing in North America.     

The EDNAP mitochondrial DNA population database (EMPOP) collaborative exercises: organisation, results and perspectives

This paper presents an overview of the organisation and the results of the collaborative exercises (CE) of the European DNA Profiling (EDNAP) Group's mitochondrial DNA population database project (EMPOP). The aim of the collaborative exercises was to determine whether uniformity of mtDNA sequencing results could be achieved among different laboratories. These were asked to sequence either the complete mtDNA control region or the two hypervariable regions HVI (16024-16365) and HVII (73-340) from DNA extracts, buccal swabs or bloodstains, proceeding in accordance with the protocol and strategies used in each individual laboratory. The results of the collaborative exercises were employed to identify possible sources of errors that could arise during the analysis and interpretation of mtDNA profiles. These findings were taken as a basis to tentatively make suitable arrangements for the construction of a high quality mtDNA database. One hundred fifty mtDNA profiles were submitted to the evaluating laboratory, and disaccording profiles were classified into four groups corresponding to the source of error: clerical errors, sample mix-ups, contaminations and discrepancies with respect to the mtDNA nomenclature. Overall, 14 disaccording haplotypes (16 individual errors) were observed. The errors included 10 clerical errors, 3 interpretation problems, 2 cases of sample mix-up and 1 case of point heteroplasmic mixture, where the 2 sequencing reactions brought inconsistent base calls. This corresponds to an error rate of 10.7% in a virtual mtDNA database consisting of the collaborative exercise results. However, this estimate is still conservative compared to conclusions drawn by authors of meanwhile numerous publications critically reviewing published mtDNA population databases. Our results and earlier published concerns strongly emphasize the need for appropriate safety regulations when mtDNA profiles are compiled for database purposes in order to accomplish the high standard required for mtDNA databases that are used in the forensic context.     

Optimization of a duplex amplification and sequencing strategy for the HVI/HVII regions of human mitochondrial DNA for forensic casework

A duplex primer set for the amplification of mitochondrial DNA HVI and HVII control regions was evaluated for the optimization of a DNA sequencing protocol suitable for forensic casework. HVI and HVII products, with the absence of non-specific products, could be detected by agarose gel electrophoresis when as little as 0.5 and 0.1pg of DNA were amplified for 34 and 38 cycles, respectively. Because HVI and HVII amplicons are co-synthesized in the duplex PCR, fewer steps are required (lessening the risk of cross contamination events) and more frugal use of precious extracted DNA samples is possible, both desirable features for forensic casework. The ABI Prism BigDyetrade mark version 1.1 chemistry provided high quality sequencing data, with little or no background noise and uniform peak heights, outcomes that favored reliable detection of heteroplasmy, particularly at early sequence reads (<40 bases). Optimal compromise between sensitivity and sequence accuracy in the absence of noise was achieved starting at 150 mitochondrial genome copies. The protocol is effective (no sequence errors) with highly degraded DNA (average detectable template size of 200bp). Dual artificial template mixtures with the minor component at 15% suggests that heteroplasmy should be detected at this level with confidence.     

Optimization of Human mtDNA Control Region Sequencing for Forensic Applications

Sequencing mitochondrial DNA hypervariable regions I and II (HVI and HVII) is useful in forensic missing person and unidentified remains cases. Improvements in ease and sensitivity of testing will yield results from more samples in a timely fashion. Routinely, amplification of HVI and HVII is followed by Sanger sequencing using the BigDye^® Terminator v3.1 Cycle Sequencing kit (Applied Biosystems) using 4 μL of ready reaction mix (RRM). Each sequencing reaction is then purified through column filtration before capillary electrophoresis. Using lower amounts of RRM (2 μL or 1 μL) and purification using BigDye^® XTerminator^? (Applied Biosystems) instead of columns showed no loss of sequence length and increased the quality and the sensitivity of testing, allowing HVI and HVII typing from mitochondrial genome equivalent to 125 fg of nuclear DNA, or 100 pg of HVI/HVII amplicons. Using this methodology, testing can be completed in 1 day, and the cost of testing is reduced.     

High‐Resolution Melting (HRM) of Hypervariable Mitochondrial DNA Regions for Forensic Science

Forensic strategies commonly are proceeding by analysis of short tandem repeats (STRs); however, new additional strategies have been proposed for forensic science. Thus, this article standardized the high‐resolution melting (HRM) of DNA for forensic analyzes. For HRM, mitochondrial DNA (mtDNA) from eight individuals were extracted from mucosa swabs by DNAzol reagent, samples were amplified by PCR and submitted to HRM analysis to identify differences in hypervariable (HV) regions I and II. To confirm HRM, all PCR products were DNA sequencing. The data suggest that is possible discriminate DNA from different samples by HRM curves. Also, uncommon dual‐dissociation was identified in a single PCR product, increasing HRM analyzes by evaluation of melting peaks. Thus, HRM is accurate and useful to screening small differences in HVI and HVII regions from mtDNA and increase the efficiency of laboratory routines based on forensic genetics.     

Necrophagous Caterpillars Provide Human mtDNA Evidence*

Abstract: Decomposition of large mammalian carcasses is greatly accelerated through the action of insects. Specialized feeders capable of digesting keratin and collagen found in skin, hair, and tendons and ligaments are attracted to corpses in late stages of dry decomposition and include Tinea pellionella, the casemaking clothes moth, and Tineola bisselliella, the webbing clothes moth (Lepidoptera; Tineidae). Until now, details of the caterpillar behavior as necrophagous insects were vague. Here, we detail the behavior of each species and document the incorporation of human hair into the portable larval shelters constructed by the caterpillars of T. pellionella. Hair of the decedent used as building material for caterpillar shelters provided enough starting template to amplify and sequence the HVI and HVII sections of the control region (mtDNA) of the decedent.     

Application of Fordisc 3.0 to Explore Differences Among Crania of North American and South African Blacks and Whites,

Using discriminant function analysis, classification accuracies for ancestry and sex in white and black South Africans were compared using North American (FDB), African groups in Howells (HDB), and South African (SADB) databases in FORDISC 3.0. (FD3). Twenty‐four standard linear measures were collected from a total of 86 black and 101 white crania obtained from the Pretoria Bone Collection. White and black South Africans classified 73% correctly in FDB, 55% correctly in HDB, and 71% correctly in SADB. The percentage of atypical cases was higher with FDB than SADB. In all three databases, misclassification occurred more with sex than ancestry revealing differences in sexual dimorphism between population groups. Broad ancestral differences may explain low misclassification rates for ancestry. FD3, with a modern South African reference sample, can assist South African anthropologists to standardize methodology and to justify procedures for estimating ancestry.     

Geographic Variation in Zygomaxillary Suture Morphology and its Use in Ancestry Estimation

Angled/curved zygomaxillary suture coding is widely employed in cranial assessments of ancestry. However, the efficacy of this method has not been extensively evaluated across diverse populations. In this study, zygomaxillary suture morphology was assessed on a total of 411 human crania from six populations (European, Native American, African, Asian, Arctic Circle, and Aboriginal Australian) using a novel 3D coordinate landmark method. Our results indicate a predominance of angled sutures among native peoples of the Arctic and North America (85–86%), a prevalence of curved sutures among Africans and Aboriginal Australians (77–81%), and essentially equal proportions of both configurations in Asians and Europeans (50–56%). Statistically, angled/curved coding generally discriminates poorly between groups, except when populations with antithetically high frequencies of the two configurations (e.g., African vs. Native American) are compared. Moreover, comparisons across previous studies reveal conflicting frequencies for many populations, further suggesting limited utility of this trait in ancestry estimation.