Sequence polymorphism of the mitochondrial DNA control region in the Slovenian population

The forensic application of mitochondrial DNA (mtDNA) typing requires large and regionally well-defined databases. To expand the database for forensic identification purposes in Slovenia, the mtDNA control region sequences of the hypervariable regions HVI and HVII were determined in a population of 129 maternally unrelated Slovenians, using a fluorescent-based capillary electrophoresis sequencing method. A total of 111 different haplotypes resulting from 124 polymorphic positions (80 polymorphic positions in HVI and 44 in HVII) were found. Of these, 101 mtDNA types were unique, 6 haplotypes were shared by 2 individuals, 1 haplotype by 3 individuals, 2 haplotypes by 4 individuals, and the most common haplotype was found in 5 individuals. The most frequent haplotypes in the Slovenian population ,263(G), 315.1(C) and 263(G), 309.1(C), 315.1(C) are also the most common in other European populations. The data support the concept that these haplotypes may represent a common European mtDNA sequence types. The sequence poymorphisms were compared to the databases of west Austria and central Italy and the HVI and HVII sequence matching probabilities within and between populations were calculated. It is 1.1–4.5 times more likely to find a sequence match in a random pair of Slovenians than in a random Slovenian-Italian pair and in a random Slovenian-Austrian pair. The length heteroplasmy in the homopolymeric C-stretch regions located at nucleotide positions 16184–16193 in HVI and at positions 303–315 in HVII was observed in 17% and 8% of individuals, respectively. A statistical estimate of the results for this population showed the random match probability and the genetic diversity of 1.16% and 0.996, respectively.Electronic Supplementary Material Supplementary material is available in the online version of this article at http://dx.doi.org/10.1007/s00414-003-0394-3     

Mitochondrial DNA sequence analysis of a native Bolivian population

Mitochondrial DNA analysis is very useful for the interpretation of the history of human migration and to estimate the frequency of a haplotype in the forensic context. From a human settlement perspective, La Paz area is greatly interesting since the first planned city of the region is located there. Samples from 110 individuals from La Paz were studied analysing the polymorphisms in the D-loop, hypervariable region I (HVI) and hypervariable region II (HVII) in order to verify the genetic diversity. The aim of this study was to start the creation of a population database in order to obtain the genetic interpopulation variability and classify haplotypes into characteristic haplogroups of South America. A total of 97 different haplotypes were identified, 90 being unique, expressed by 122 polymorphic nucleotide positions. Nucleotide and sequence diversity were estimated to be 0.015 ± 0.0075 and 0.996, respectively. Haplogroup distribution in the samples was 57.27% B4, 19.09% C1, 10.00% A2, 3.64% D1, 2.73% D4h3, 1.82% H, and 0.91% for each of the haplogroups A4, B4c1a, CZ, D4J, M7a and M8/N9b. The rate of length heteroplasmy was 36.36% in HVI and 52.73% in HVII. Phylogenetic analysis reveals proximity to the Korean, Chilean aboriginal, Japanese and Australian populations. The estimated genetic variability of the studied population was high, suggesting an early settlement.     

Sequence polymorphisms of the mitochondrial DNA control region and phylogenetic analysis of mtDNA lineages in the Japanese population

Sequence polymorphisms of the hypervariable mitochondrial DNA (mtDNA) regions HVI and HVII, and coding region polymorphisms were investigated in 211 unrelated individuals from the Japanese population. Sequence comparison of the HVI and HVII regions led to the identification of 169 mitochondrial haplotypes defined by 147 variable positions. Among them 145 types were observed in only 1 individual; the other 24 types were shared by 2 or more individuals. The gene diversity was estimated at 0.9961, and the probability of two randomly selected individuals from the population having identical mtDNA types was 0.86%. We also established phylogenetic haplogroups in the Japanese population based on the coding and control region polymorphisms and compared the haplotypes with those in other Japanese, Korean and Chinese populations. As a result, three new subhaplogroups, G4a, G4b, and N9b, and several haplotypes specific for the Japanese and Korean populations were identified. The present database can be used not only for personal identification but also as an aid for geographic or phenotype (race) estimation in forensic casework in Japan.Electronic Supplementary Material Supplementary material is available for this article if you access the article at     

HVI and HVII mitochondrial DNA data in Apaches and Navajos

Most mtDNA studies on Native Americans have concentrated on hypervariable region I (HVI) sequence data. Mitochondrial DNA haplotype data from hypervariable regions I and II (HVI and HVII) have been compiled from Apaches (N=180) and Navajos (N=146). The inclusion of HVII data increases the amount of information that can be obtained from low diversity population groups. Less mtDNA variation was observed in the Apaches and Navajos than in major population groups. The majority of the mtDNA sequences were observed more than once; only 17.8% (32/180) of the Apache sequences and 25.8% of the Navajo sequences were observed once. Most of the haplotypes in Apaches and Navajos fall into the A and B haplogroups. Although a limited number of haplogroups were observed, both sample populations exhibit sufficient variation for forensic mtDNA typing. Genetic diversity was 0.930 in the Apache sample and 0.963 in the Navajo sample. The random match probability was 7.48% in the Apache sample and 4.40% in the Navajo sample. The average number of nucleotide differences between individuals in a database is 9.0 in the Navajo sample and 7.7 in the Apache sample. The data demonstrate that mtDNA sequencing can be informative in forensic cases where Native American population data are used.     

Mitochondrial DNA population data of HV1 and HV2 sequences from Japanese individuals

Mitochondrial DNA sequences of the hypervariable regions HV1 and HV2 were determined for 1204 unrelated Japanese individuals. A total of 741 different mtDNA haplotypes were found, 157 of which were seen in multiple individuals. Twenty-seven of these individuals showed point heteroplasmy. The most frequent haplotype (16223T-16362C-73G-263G-315.1C) was found in 31 individuals and the second most frequent haplotypes (16129A-16223T- 16362C-73G-152C-263G-309.1C-315.1C) was found in 24 individuals. The haplotypes diversity and random match probability were calculated to be 0.9969 and 0.0040, respectively.     

Forensic mitochondrial coding region analysis for increased discrimination using pyrosequencing technology

Analysis of mitochondrial DNA (mtDNA) is very useful when nuclear DNA analysis fails due to degradation or insufficient amounts of DNA in forensic analysis. However, mtDNA analysis has a lower discrimination power compared to what can be obtained by nuclear DNA (nDNA) analysis, potentially resulting in multiple individuals showing identical mtDNA types in the HVI/HVII region. In this study, the increase in discrimination by analysis of mitochondrial coding regions has been evaluated for identical or similar HVI/HVII sequences. A pyrosequencing-based system for coding region analysis, comprising 17 pyrosequencing reactions performed on 15 PCR fragments, was utilised. This assay was evaluated in 135 samples, resulting in an average read length of 81 nucleotides in the pyrosequencing analysis. In the sample set, a total of 52 coding region SNPs were identified, of which 18 were singletons. In a group of 60 samples with 0 or 1 control region difference from the revised Cambridge reference sequence (rCRS), only 12 samples could not be resolved by at least two differences using the pyrosequencing assay. Thus, the use of this pyrosequencing-based coding region assay has the potential to substantially increase the discriminatory power of mtDNA analysis.     

Sequence polymorphism of mitochondrial DNA in Japanese individuals from Gifu Prefecture

Sequence polymorphisms of the hypervariable region HV1 in mitochondrial DNA (mtDNA) were analyzed in a sample of 137 unrelated Japanese individuals living in Gifu Prefecture (central region of Japan) using polymerase chain reaction amplification and direct sequencing. Eighty-two different haplotypes resulting from 81 variable sites were found in the mtDNA HV1 region between positions 16061 and 16450. The most frequent haplotype (16223T, 16362C) was shared by ten individuals. The genetic diversity and the genetic identity were 0.985 and 0.022, respectively. The C-stretch region located around position 16189 was observed in 23.4% of this population sample. Sequence heteroplasmy at the position 16103 (A/G) was found in one individual.     

Sequence polymorphism of the mitochondrial DNA hypervariable regions I and II in 205 Singapore Malays

Mitochondrial DNA sequences of the hypervariable regions HV1 and HV2 were analyzed in 205 unrelated ethnic Malays residing in Singapore as an initial effort to generate a database for forensic identification purposes. Sequence polymorphism was detected using PCR and direct sequencing analysis. A total of 152 haplotypes was found containing 152 polymorphisms. Out of the 152 haplotypes, 115 were observed only once and 37 types were seen in multiple individuals. The most common haplotype (16223T, 16295T, 16362C, 73G, 146C, 199C, 263G, and 315.1C) was shared by 7 (3.41%) individuals, two haplotypes were shared by 4 individuals, seven haplotypes were shared by 3 individuals, and 27 haplotypes by 2 individuals. Haplotype diversity and random match probability were estimated to be 0.9961% and 0.87%, respectively.     

Investigation of length heteroplasmy in mitochondrial DNA control region by massively parallel sequencing

Accurate sequencing of the control region of the mitochondrial genome is notoriously difficult due to the presence of polycytosine bases, termed C-tracts. The precise number of bases that constitute a C-tract and the bases beyond the poly cytosines may not be accurately defined when analyzing Sanger sequencing data separated by capillary electrophoresis. Massively parallel sequencing has the potential to resolve such poor definition and provides the opportunity to discover variants due to length heteroplasmy. In this study, the control region of mitochondrial genomes from 20 samples was sequenced using both standard Sanger methods with separation by capillary electrophoresis and also using massively parallel DNA sequencing technology. After comparison of the two sets of generated sequence, with the exception of the C-tracts where length heteroplasmy was observed, all sequences were concordant. Sequences of three segments 16184–16193, 303–315 and 568–573 with C-tracts in HVI, II and III can be clearly defined from the massively parallel sequencing data using the program SEQ Mapper. Multiple sequence variants were observed in the length of C-tracts longer than 7 bases. Our report illustrates the accurate designation of all the length variants leading to heteroplasmy in the control region of the mitochondrial genome that can be determined by SEQ Mapper based on data generated by massively parallel DNA sequencing.     

Full mtGenome reference data: Development and characterization of 588 forensic-quality haplotypes representing three U.S. populations

Though investigations into the use of massively parallel sequencing technologies for the generation of complete mitochondrial genome (mtGenome) profiles from difficult forensic specimens are well underway in multiple laboratories, the high quality population reference data necessary to support full mtGenome typing in the forensic context are lacking. To address this deficiency, we have developed 588 complete mtGenome haplotypes, spanning three U.S. population groups (African American, Caucasian and Hispanic) from anonymized, randomly-sampled specimens. Data production utilized an 8-amplicon, 135 sequencing reaction Sanger-based protocol, performed in semi-automated fashion on robotic instrumentation. Data review followed an intensive multi-step strategy that included a minimum of three independent reviews of the raw data at two laboratories; repeat screenings of all insertions, deletions, heteroplasmies, transversions and any additional private mutations; and a check for phylogenetic feasibility. For all three populations, nearly complete resolution of the haplotypes was achieved with full mtGenome sequences: 90.3–98.8% of haplotypes were unique per population, an improvement of 7.7–29.2% over control region sequencing alone, and zero haplotypes overlapped between populations. Inferred maternal biogeographic ancestry frequencies for each population and heteroplasmy rates in the control region were generally consistent with published datasets. In the coding region, nearly 90% of individuals exhibited length heteroplasmy in the 12418-12425 adenine homopolymer; and despite a relatively high rate of point heteroplasmy (23.8% of individuals across the entire molecule), coding region point heteroplasmies shared by more than one individual were notably absent, and transversion-type heteroplasmies were extremely rare. The ratio of nonsynonymous to synonymous changes among point heteroplasmies in the protein-coding genes (1:1.3) and average pathogenicity scores in comparison to data reported for complete substitutions in previous studies seem to provide some additional support for the role of purifying selection in the evolution of the human mtGenome. Overall, these thoroughly vetted full mtGenome population reference data can serve as a standard against which the quality and features of future mtGenome datasets (especially those developed via massively parallel sequencing) may be evaluated, and will provide a solid foundation for the generation of complete mtGenome haplotype frequency estimates for forensic applications.     

Forensic evaluation of mtDNA in a population from south west Switzerland

The polymorphism of the two hypervariable segments (HVI and HVII) of the control region of mtDNA was analyzed in a population of 154 unrelated individuals from south west Switzerland using a fluorescent based capillary electrophoresis sequencing method. In our population data of 154 random individuals, 137 mtDNA types were observed. Of these, 124 sequences were observed only in one individual whereas 10 sequences were observed in 2 individuals, 2 sequences in 3 individuals and 1 sequence in 4 individuals. The probability of two unrelated individuals having the same sequence was 0.84%. The results were compared with four other Caucasian populations. Furthermore, the usefulness of the mtDNA sequencing was tested, for exclusion and inclusion, in 18 forensic cases including 69 evidence samples and 44 reference samples. Despite the fact that 55% of the evidence samples yielded a negative result for the nuclear DNA with the human dot quantitation system, the success rate of the mtDNA sequencing was 71.0%. This validation study proves the great usefulness and sensitivity of the mtDNA sequencing technique using nested PCR and fluorescent capillary electrophoresis.     

Different methods to determine length heteroplasmy within the mitochondrial control region

The first and second hypervariable regions of the human mitochondrial DNA control region contain two homopolymeric stretches of cytosine (nt 16184–16193 and nt 303–315, respectively). According to the Cambridge reference sequence these homopolymeric stretches are interrupted by thymine (T), at positions 16189 and 310, respectively. Monotonous runs of the same base have been suggested to be hot spots for mutations, probably caused by replication slippage, resulting in length heteroplasmy. This paper describes a rapid method based on restriction cleavage of labelled PCR products encompassing the homopolymeric tract in HVII to quantify the relative proportions of different length variants present in an individual. To compare the accuracy of this method, cloned PCR products from several heteroplasmic individuals have been additionally sequenced.     

Identification of single nucleotide polymorphisms within the mtDNA genome of the domestic dog to discriminate individuals with common HVI haplotypes

We sequenced the entire ～16 kb canine mitochondrial genome (mtGenome) of 100 unrelated domestic dogs (Canis lupus familiaris) and compared these to 246 published sequences to assess hypervariable region I (HVI) haplotype frequencies. We then used all available sequences to identify informative single nucleotide polymorphisms (SNPs) outside of the control region for use in further resolving mtDNA haplotypes corresponding to common HVI haplotypes. Haplotype frequencies in our data set were highly correlated with previous ones (e.g., F(ST)=0.02, r=0.90), suggesting the total data set reasonably reflected the broader dog population. A total of 128 HVI haplotypes was represented. The 10 most common HVI haplotypes (n=184 dogs) represented 53.3% of the sample. We identified a total 71 SNPs in the mtGenomes (external to the control region) that resolved the 10 most common HVI haplotypes into 63 mtGenome subhaplotypes. The random match probability of the dataset based solely on the HVI sequence was 4%, whereas the random match probability of the mtGenome subhaplotypes was <1%. Thus, the panel of 71 SNPs identified in this study represents a useful forensic tool to further resolve the identity of individual dogs from mitochondrial DNA (mtDNA).     

Mitochondrial heteroplasmy among maternally related individuals

The second hypervariable segment of the human mtDNA control region contains a homopolymeric tract of cytidines between nucleotides (nt) 303 and 315, interrupted by a thymidine at position 310, according to the Cambridge reference sequence. By direct sequencing, some individuals show blurred sequence chromatograms in this region which are not caused by a sequencing artefact but by high levels of length heteroplasmy. With respect to this length heteroplasmy ten maternally related individuals and two unrelated probands were examined. The relative proportions of length variants in the homopolymeric tract in selected individuals were determined by cloning and sequencing of multiple independent clones. All ten family members examined were heteroplasmic while the proportions of each genotype varied widely in different individuals. The size of a possible mitochondrial bottleneck during embryonic development of the offspring is discussed with respect to the changes in mitochondrial haplotypes within mother-offspring pairs. Our data are consistent with both slow and rapid segregation of mtDNAs between the generations, which would implicate a tight as well as a wide bottleneck. Therefore, a common bottleneck size in all individuals from this lineage seems to be very unlikely. Received: 25 February 1999 / Received in revised form: 7 June 1999     

Mitochondrial DNA sequences for 118 individuals from northeastern Spain

A population database was generated from 118 unrelated Caucasoid individuals living in Spain. Sequence polymorphisms of the mitochondrial DNA (mtDNA) control region, hypervariable regions I and II (HVRI and HVRII) were determined using the polymerase chain reaction (PCR) and direct sequencing. A total of 102 different sequences were found as defined by 105 variable positions. The most common sequence occurred six times, and this sequence is also the most frequent in other European populations such as Austria, Germany and Britain. The mean pair-wise difference for the two HVR regions taken together was 7.74. The study revealed that transitions made up the majority of the variations (88%), whereas we observed a significantly lower frequency of transversions (8%). Also one individual in this study was observed with two positions of heteroplasmy at nucleotides 150 (C/T) and 153 (G/A). A statistical estimate of the results for this population showed a genetic diversity of 0.99. The probability of two random individuals showing identical mtDNA haplotypes is 1.3%. In order to use the mtDNA analysis in forensic casework, we consider that it is of crucial importance to know the frequency of the different sequences of mtDNA, and this data base study could be a useful tool to statistically evaluate the results. Received: 12 July 1999 / Accepted: 11 April 2000     

A new database of mitochondrial DNA hypervariable regions I and II sequences from 162 Japanese individuals

A database of mitochondrial DNA (mtDNA) hypervariable region 1 (HV1) and region 2 (HV2) sequences of the mtDNA control region was established from 162 unrelated Japanese individuals. The random match probability and the genetic diversity for this database were 0.96% and 0.997, respectively. Length heteroplasmy in the C-stretch regions located around position 16189 in HV1 and 310 in HV2 was observed in 37% and 38% of the samples, respectively. A strategy using internal sequencing primers was devised to obtain confirmed sequences in these length heteroplasmic individuals. This database, combined with other mtDNA sequence databases from the Japanese population, will permit the significance of mtDNA match results to be properly reported in mtDNA typing casework in Japan.Received: 10 August 2000 / Accepted: 7 January 2001An erratum to this article can be found at     

Mitochondrial DNA control region diversity in hairs and body fluids of monozygotic triplets

Length heteroplasmy of the homopolymeric cytosine stretch in the hypervariable region II of the mitochondrial D-loop was investigated in blood, buccal cells and hair shafts of monozygotic triplets. The proportions of length heteroplasmy were determined by cloning and sequencing of multiple independent clones. Blood and buccal cells showed an accumulation of molecules with one and two insertions of cytosine residues in relation to the Cambridge Reference Sequence (CRS). The results did not show statistically significant differences between blood and buccal cells of one and the same individual and also not between the three monozygotic brothers. In the hair samples a loss of cytosine residues was established in all three monozygotic individuals compared to blood and buccal cells, suggesting that this must be a regular process. Furthermore, the hair shaft samples showed significant differences between the frequencies of 7, 8 or 9 Cs in the poly C region comparing the three individuals (p<0.008) and in addition there were highly significant differences (p<0.0001) when comparing the results for six different hairs of each individual separately. From these results it can be assumed that besides a common genetic bottleneck during embryonic development, a post-embryonic bottleneck seems to exist in each hair follicle.     

Sequence polymorphisms of mtDNA HV1, HV2, and HV3 regions in the Malay population of Peninsular Malaysia

The uniparentally inherited mitochondrial DNA (mtDNA) is in the limelight for the past two decades, in studies relating to demographic history of mankind and in forensic kinship testing. In this study, human mtDNA hypervariable segments 1, 2, and 3 (HV1, HV2, and HV3) were analyzed in 248 unrelated Malay individuals in Peninsular Malaysia. Combined analyses of HV1, HV2, and HV3 revealed a total of 180 mtDNA haplotypes with 149 unique haplotypes and 31 haplotypes occurring in more than one individual. The genetic diversity was estimated to be 99.47%, and the probability of any two individuals sharing the same mtDNA haplotype was 0.93%. The most frequent mtDNA haplotype (73, 146, 150, 195, 263, 315.1C, 16140, 16182C, 16183C, 16189, 16217, 16274, and 16335) was shared by 11 (4.44%) individuals. The nucleotide diversity and mean of pair-wise differences were found to be 0.036063 ± 0.020101 and 12.544022 ± 6.230486, respectively.     

Molecular genetic identification of skeletal remains from the Second World War Konfin I mass grave in Slovenia

This paper describes molecular genetic identification of one third of the skeletal remains of 88 victims of postwar (June 1945) killings found in the Konfin I mass grave in Slovenia. Living relatives were traced for 36 victims. We analyzed 84 right femurs and compared their genetic profiles to the genetic material of living relatives. We cleaned the bones, removed surface contamination, and ground the bones into powder. Prior to DNA isolation using Biorobot EZ1 (Qiagen), the powder was decalcified. The nuclear DNA of the samples was quantified using the real-time polymerase chain reaction method. We extracted 0.8 to 100 ng DNA/g of bone powder from 82 bones. Autosomal genetic profiles and Y-chromosome haplotypes were obtained from 98% of the bones, and mitochondrial DNA (mtDNA) haplotypes from 95% of the bones for the HVI region and from 98% of the bones for the HVII region. Genetic profiles of the nuclear and mtDNA were determined for reference persons. For traceability in the event of contamination, we created an elimination database including genetic profiles of the nuclear and mtDNA of all persons that had been in contact with the skeletal remains. When comparing genetic profiles, we matched 28 of the 84 bones analyzed with living relatives (brothers, sisters, sons, daughters, nephews, or cousins). The statistical analyses showed a high confidence of correct identification for all 28 victims in the Konfin I mass grave (posterior probability ranged from 99.9% to more than 99.999999%).     

Sequence analysis of mitochondrial DNA hypervariable region I in Thai individuals

Nucleotide sequence analysis of hypervariable region I (HVRI) in human mitochondrial DNA (mtDNA) was investigated in 100 unrelated Thai individuals. A total of 85 variable sites and 423 base substitutions, which consisted of 390 nucleotide transitions and 33 nucleotide transversions were found. The following nucleotide substitutions were found: 48% at 16,223, 31% at 16,304, 30% at 16,332, and 26% at 16,129, respectively. Transition from T to C (43.7%) was the most frequent substitution. The nucleotide insertions were found at two sites with T at position 16,188 and C at position 16,194. Eighty-two haplotypes were investigated of which 72 haplotypes were unique. The most frequent haplotypes (16,108T-16,129A-16,162G-16,172C-16,304C and 16,260T-16,298C-16,355T-16,362C) were observed. From position 16,180 to 16,193, thirteen patterns of polycytosine or C-stretch were observed, whereas 68 Thai individuals were found to be similar to the references. The genetic diversity, random match probability, and discrimination power were estimated to be 0.9943, 0.0156, and 0.9844, respectively.