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

Mitochondrial DNA control region diversity in a population from Parana state-increasing the Brazilian forensic database

The entire mitochondrial DNA (mtDNA) control region (nucleotide position 16024-576) sequences were obtained through Sanger sequencing method for 122 individuals from Parana state, South of Brazil. We observed a total of 108 different haplotypes of which 97 were unique and 11 were shared by more than one individual. The haplogroups were classified according to the updated mtDNA phylogeny, by EMMA (estimating mitochondrial haplogroups using a maximum likelihood approach). Our results revealed the predominance of Amerindian haplogroups with a frequency of 49.2% of the population sample, followed by European lineages with 38.5% and 12.3% of African lineages. Parana population sample set presented a high haplotype diversity (0.9976) and the random match probability was 0.0106. The phylogenetical findings and the diversity indices confirm the high genetic heterogeneity of this population and suggest a high informativeness of mtDNA analyses in forensic cases. The population data will contribute to increase the Brazilian mtDNA database for forensic purposes and it is available through EMPOP (European DNA Profiling Group mitochondrial DNA population database) under the accession number EMP00714.

Detection of somatic mutations in the mitochondrial DNA control region D-loop in brain tumors: The first report in Malaysian patients

Although the role of nuclear-encoded gene alterations has been well documented in brain tumor development, the involvement of the mitochondrial genome in brain tumorigenesis has not yet been fully elucidated and remains controversial. The present study aimed to identify mutations in the mitochondrial DNA (mtDNA) control region D-loop in patients with brain tumors in Malaysia. A mutation analysis was performed in which DNA was extracted from paired tumor tissue and blood samples obtained from 49 patients with brain tumors. The D-loop region DNA was amplified using the PCR technique, and genetic data from DNA sequencing analyses were compared with the published revised Cambridge sequence to identify somatic mutations. Among the 49 brain tumor tissue samples evaluated, 25 cases (51%) had somatic mutations of the mtDNA D-loop, with a total of 48 mutations. Novel mutations that had not previously been identified in the D-loop region (176 A-deletion, 476 C>A, 566 C>A and 16405 A-deletion) were also classified. No significant associations between the D-loop mutation status and the clinicopathological parameters were observed. To the best of our knowledge, the current study presents the first evidence of alterations in the mtDNA D-loop regions in the brain tumors of Malaysian patients. These results may provide an overview and data regarding the incidence of mitochondrial genome alterations in Malaysian patients with brain tumors. In addition to nuclear genome aberrations, these specific mitochondrial genome alterations may also be considered as potential cancer biomarkers for the diagnosis and staging of brain cancers.

Human genetic data reveal contrasting demographic patterns between sedentary and nomadic populations that predate the emergence of farming

Demographic changes are known to leave footprints on genetic polymorphism. Together with the increased availability of large polymorphism data sets, coalescent-based methods allow inferring the past demography of populations from their present-day patterns of genetic diversity. Here, we analyzed both nuclear (20 noncoding regions) and mitochondrial (HVS-I) resequencing data to infer the demographic history of 66 African and Eurasian human populations presenting contrasting lifestyles (nomadic hunter-gatherers, nomadic herders, and sedentary farmers). This allowed us to investigate the relationship between lifestyle and demography and to address the long-standing debate about the chronology of demographic expansions and the Neolithic transition. In Africa, we inferred expansion events for farmers, but constant population sizes or contraction events for hunter-gatherers. In Eurasia, we inferred higher expansion rates for farmers than herders with HVS-I data, except in Central Asia and Korea. Although isolation and admixture processes could have impacted our demographic inferences, these processes alone seem unlikely to explain the contrasted demographic histories inferred in populations with different lifestyles. The small expansion rates or constant population sizes inferred for herders and hunter-gatherers may thus result from constraints linked to nomadism. However, autosomal data revealed contraction events for two sedentary populations in Eurasia, which may be caused by founder effects. Finally, the inferred expansions likely predated the emergence of agriculture and herding. This suggests that human populations could have started to expand in Paleolithic times, and that strong Paleolithic expansions in some populations may have ultimately favored their shift toward agriculture during the Neolithic.

Mitochondrial haplogroup B4 may be a protective factor to oral lichen planus susceptibility in Chinese

OBJECTIVE

Oral lichen planus (OLP) is a common inflammatory disorder with a higher prevalence among women than men. The pathogenesis of OLP is still unclear, and its heredity is not well characterized. Maternal inheritance of mitochondrial DNA (mtDNA) indicates its importance in human ethnic group classification as well disease susceptibility. So, this study was conducted to find whether mtDNA haplogroup associates with OLP susceptibility in Chinese, and determine the influence of C-stretch structure of mtDNA on OLP susceptibility.

SUBJECTS AND METHODS

We sequenced three hypervariable regions of mtDNA from 242 patients with OLP and 237 healthy controls. The association between mtDNA haplogroups and OLP is tested by chi-square test.

RESULTS

Seventeen kinds of haplogroups were identified, and the frequency of haplogroup B4 was significantly higher in control group than that of the OLP group (P = 0.013, OR=0.429). In female samples, B4 declared even more significance (P = 0.003, OR=0.296). Meanwhile, the haplotypes of C-stretch in mtDNA did not have any significant difference between case-control groups.

CONCLUSIONS

Mitochondrial DNA haplogroup B4 might have a protective effect to OLP, and its protective effect resides predominantly in women. However, the association between C-stretch haplotype and OLP susceptibility still needs more patients for evaluation.

Analysis of human mitochondrial DNA polymorphisms in the Japanese population

The highly polymorphic nature and high amplification efficiency of mitochondrial DNA (mtDNA) is valuable for the analysis of biological evidence in forensic casework, such as the identification of individuals and assignment of race/ethnicity. To be useful, a mtDNA polymorphism database for the Japanese population requires an understanding of the range of haplotype variation and phylogenies of mtDNA sequences. To extend current knowledge on the haplotypes in the Japanese population, this study defines new lineages and provides more detail about some of those previously described. We compared the hypervariable regions (HVRs) of 270 healthy, unrelated Japanese individuals and demonstrated 192 haplotypes. Combining HVR1 and HVR2, the genetic diversity was 0.9935, thus providing a high level of identification capability. Haplogroup status was defined for 160 individuals using HVR1, HVR2, and particular coding region polymorphisms; these individuals belonged to 94 haplotypes, four of which were new lineages. The complete mtDNA sequence was also determined from seven individuals.

Hypervariable region polymorphism of mtDNA of recurrent oral ulceration in Chinese

Background

MtDNA haplogroups could have important implication for understanding of the relationship between the mutations of the mitochondrial genome and diseases. Distribution of a variety of diseases among these haplogroups showed that some of the mitochondrial haplogroups are predisposed to disease. To examine the susceptibility of mtDNA haplogroups to ROU, we sequenced the mtDNA HV1, HV2 and HV3 in Chinese ROU.

Methodology/Principal Findings

MtDNA haplogroups were analyzed in the 249 cases of ROU patients and the 237 cases of healthy controls respectively by means of primer extension analysis and DNA sequencing. Haplogroups G1 and H were found significantly more abundant in ROU patients than in healthy persons, while haplogroups D5 and R showed a trend toward a higher frequency in control as compared to those in patients. The distribution of C-stretch sequences polymorphism in mtDNA HV1, HV2 and HV3 regions was found in diversity.

Conclusions/Significance

For the first time, the relationship of mtDNA haplogroups and ROU in Chinese was investigated. Our results indicated that mtDNA haplogroups G1 and H might constitute a risk factor for ROU, which possibly increasing the susceptibility of ROU. Meanwhile, haplogroups D5 and R were indicated as protective factors for ROU. The polymorphisms of C-stretch sequences might being unstable and influence the mtDNA replication fidelity.

A new strategy for the discrimination of mitochondrial DNA haplogroups in Han population

Mitochondrial DNA (mtDNA) haplogroup discrimination is interesting not only for phylogenetic and clinical but also for forensic studies. We discriminated the mtDNA haplogroups of 570 healthy unrelated Han people from Zhejiang Province, Southeast China, by comprehensive analysis mutations of the hypervariable segments-I sequence and diagnostic polymorphisms in mtDNA coding region using real-time polymerase chain reaction (RT-PCR), which was compared with the widely used PCR and restriction fragment length polymorphism (PCR-RFLP) method. The results showed that in superhaplogroup M, haplogroup D was the most common haplotype within this assay to 24.6%, and in the other superhaplotype N, haplogroup B and F were the most common groups. Samples re-identified by PCR-RFLP showed the consistent results that were got with RT-PCR. In conclusion, the RT-PCR strategy appears to be an accurate, reproducible, and sensitive technique for the discrimination of mtDNA haplogroups, especially for mass screenings quickly and economically.

Alterations of length heteroplasmy in mitochondrial DNA under various amplification conditions

There are several areas within mitochondrial DNA that show length heteroplasmy. If the heteroplasmy pattern is unique and consistent for each person, it may be used to support an interpretation of exclusion in identity testing. We investigated whether the length heteroplasmy pattern would be consistent under different amplification conditions. We also determined whether various amplification parameters would affect the homopolymeric cytosine stretches (C-stretch) in HV1. Monoclonal samples tended to be heteroplasmic after amplification. After several repetitions, C-stretch patterns of all samples were inconsistent even under the same amplification conditions. Increased PCR cycles and high template concentrations resulted in a more frequent heteroplasmic tendency. These amplification parameters seem to have little effect if samples are not long enough in C-stretch or total length of the segment from nt 16180 to nt 16193. It is suggested that the pattern of length heteroplasmy cannot be used as an additional polymorphic marker.

Forensic and phylogeographic characterization of mtDNA lineages from northern Thailand (Chiang Mai)

The immigration of diverse ethnic groups over the past centuries from surrounding countries into Thailand left footprints in the genetic composition of Thai mitochondrial DNA (mtDNA) lineages. The entire mtDNA control region (1,122 bp) was typed in 190 unrelated male volunteers from the northern Thailand province of Chiang Mai following highest quality standards. For a more precise haplogroup classification, selected single nucleotide polymorphisms from the mtDNA coding region were genotyped. We found several new, so far undescribed mtDNA lineages. Quasi-median networks were constructed for visualisation of character conflicts. The data were put into population-genetic relationships with other Southeast Asian populations. Although the frequencies of the Thai haplogroups were characteristic for Southeast Asia in terms of haplotype composition and genetic structure, the Thai population was significantly different from other Southeast Asian populations. This necessitates establishing regional databases, especially for forensic applications. The population data have been submitted to the EMPOP database (www.empop.org) and will be available on publication.

Sequence-length variation of mtDNA HVS-I C-stretch in Chinese ethnic groups

The purpose of this study was to investigate mitochondrial DNA (mtDNA) hypervariable segment-I (HVS-I) C-stretch variations and explore the significance of these variations in forensic and population genetics studies. The C-stretch sequence variation was studied in 919 unrelated individuals from 8 Chinese ethnic groups using both direct and clone sequencing approaches. Thirty eight C-stretch haplotypes were identified, and some novel and population specific haplotypes were also detected. The C-stretch genetic diversity (GD) values were relatively high, and probability (P) values were low. Additionally, C-stretch length heteroplasmy was observed in approximately 9% of individuals studied. There was a significant correlation (r=-0.961, P<0.01) between the expansion of the cytosine sequence length in the C-stretch of HVS-I and a reduction in the number of upstream adenines. These results indicate that the C-stretch could be a useful genetic maker in forensic identification of Chinese populations. The results from the Fst and dA genetic distance matrix, neighbor-joining tree, and principal component map also suggest that C-stretch could be used as a reliable genetic marker in population genetics.

Study of the cytochrome b gene sequence in populations of Taiwan

The cytochrome b gene (MTCYB) has been widely used in taxonomic research. In this study, the sequence polymorphism of the MTCYB gene was determined in 417 subjects of eight populations living in Taiwan (Taiwanese Han, indigenous Taiwanese, Tao, mainland Chinese, Filipino, Thai, Vietnamese, and Caucasian). Sequence variation from the revised Cambridge Reference Sequence and genetic distance between these populations were analyzed. There were 108 variable positions with a total of 99 haplotypes. Population-specific positions of MTCYB gene were noted in Tao and Caucasian populations. There were statistically significant differences of genetic distance between Taiwanese Han and Caucasian, between Taiwanese Han and Tao, and between Taiwanese Han and Filipino. A phylogenetic tree presents the genetic distances between these populations. In conclusion, there are sufficient sequence polymorphisms of the MTCYB gene in individuals of different populations, which may be used in the analyses of human ethnic groups in forensic casework.

MtDNA control region sequence polymorphisms and phylogenetic analysis of Malay population living in or around Kuala Lumpur in Malaysia

Control region polymorphisms in the mitochondrial DNA of 124 unrelated individuals from the Malay population living in or around Kuala Lumpur in Malaysia were investigated and phylogenetic haplogroup lineages were determined. The intergenic COII/tRNALys 9-bp deletion, 3010 and 5178 mutations, and several coding region polymorphisms were examined to discriminate some phylogenetic haplogroups. Sequence comparison of the control regions led to the identification of 117 mitochondrial haplotypes, in which 103 types were observed in only one individual and the other nine types were shared by more than two individuals. Gene diversity was estimated to be 0.997. Phylogenetic haplogroup determination revealed that the gene pool of the modern Malay population in Malaysia consisted mainly of southeast Asian, east Asian, unidentified and unique, and aboriginal southeast-specific haplogroups. These results suggest a multi-original nature for the modern Malay population. The present database may help not only in personal identification but also in determining geographic origin in forensic casework in Malaysian, Southeast Asian and East Asian populations.

Genetic studies of eight X-STRs in a Japanese population

We studied eight X-STRs (DXS7132, DXS7423, DXS8378, DXS10074, DXS10101, DXS10134, DXS10135, HPRTB) polymorphism in 494 unrelated Japanese individuals (313 males, 181 females) using Mentype Argus X-8 PCR Amplification Kit. PD of the eight X-STRs ranged from 0.558 (male) to 0.987 (female). Allele frequencies, number of alleles, and PIC were 0.001-0.587, 6-20, and 0.470-0.913, respectively.

Mitochondrial DNA control region polymorphism in the population of Alagoas state, north-eastern Brazil

The sequences of the two hypervariable (HV) segments of the mitochondrial DNA (mtDNA) control region were determined in 167 randomly selected, unrelated individuals living in the state of Alagoas, north-eastern Brazil. One hundred and forty-five different haplotypes, associated with 139 variable positions, were determined. More than 95% of the mtDNA sequences could be allocated to specific mtDNA haplogroups according to the mutational motifs. Length heteroplasmy in the C-stretch HV1 and HV2 regions was observed in 22 and 11%, respectively, of the population sample. The genetic diversity was estimated to be 0.9975 and the probability of two random individuals presenting identical mtDNA haplotypes was 0.0084. The most frequent haplotype was shared by six individuals. All sequences showed high-quality values and phantom mutations were not detected. The diversity revealed in the mitochondrial control region indicates the importance of this locus for forensic casework and population studies within Alagoas, Brazil.

mtDNA haplogroup distribution in Chinese patients with Leber's hereditary optic neuropathy and G11778A mutation

Mitochondrial DNA background has been shown to be involved in the penetrance of Leber's hereditary optic neuropathy (LHON) in western Eurasian populations. To analyze mtDNA haplogroup distribution pattern in Han Chinese patients with LHON and G11778A mutation, we analyzed the mtDNA haplogroups of 41 probands with LHON known to harbor G11778A mutation by sequencing the mtDNA control region hypervariable segments and some coding region polymorphisms. Each mtDNA was classified according to the available East Asian haplogroup system. The haplogroup distribution pattern of LHON sample was then compared to the reported Han Chinese samples. Haplogroups M7, D, B, and A were detected in 11 (26.8%), 10 (24.4%), 8 (19.5%), and 5 (12.2%) LHON families, respectively, and accounted for 82.9% of the total samples examined. For the remaining seven mtDNAs, six belonged to M8a, M10a, C, N9a, F1a, and R11, respectively, and one could only be assigned into macro-haplogroup M. The LHON sample was distinguished from other Han Chinese samples in the principal component map based on haplogroup distribution frequency. Our results show that matrilineal genetic components of Chinese LHON patients with G11778A are diverse and differ from related Han Chinese regional samples. mtDNA background might affect the expression of LHON and the G11778A mutation in Chinese population.

Sequence analysis of mitochondrial DNA hypervariable regions: an approach to personal identification

BACKGROUND

Mitochondrial DNA (mtDNA) sequence analysis of the hypervariable control region has been shown to be an effective tool for personal identification. The high copy and maternal mode of inheritance make mtDNA analysis particularly useful when old samples or degradation of biological samples prohibits the detection of nuclear DNA analysis.

METHODS

This study reports mtDNA polymorphisms in both HV1 and HV2 of the non-coding D-loop region from 30 Iranian persons (10 unrelated families in three sequential maternal generations).

RESULTS

Eighty one polymorphic nucleotide positions were found, 32 in HV1 and 49 in HV2. The sequence of HV1 and HV2 and the occurred polymorphism were completely similar in each family, except for heteroplasmy in five positions in HV2. The average numbers of nucleotide differences between families were 5.2 nucleotides in HV1 and 2.8 nucleotides in HV2.

CONCLUSIONS

We expect eight nucleotide differences in sequence of both HV1 and HV2 in two unknown, unrelated Iranian samples.

Migration rates and genetic structure of two Hungarian ethnic groups in Transylvania, Romania

Transylvania's ethnic mosaic is composed of Romanians, German Saxons and Hungarians. The ethnic groups of the Hungarian minority that settled in Romania show differences in dialects, customs and religious affiliations. In this study entire mtDNA control region sequences from 360 individuals of Hungarian ethnicity from two populations (the Csángó and the Székely), settled in the historical region of Transylvania in Romania, were generated and analyzed following high quality sequencing standards. Phylogenetic analyses were used for haplogroup determination, quasi-median network analyses were applied for the visualization of character conflicts, and median joining reconstructions were used for depicting haplotype structures. Affiliation of haplotypes to major west Eurasian haplogroups was confirmed using coding region SNPs. Gene flow between the two populations was low and biased towards a higher migration rate from the Csángó to the Székely than vice versa. Phylogeographic analyses revealed effects of genetic isolation within the Csángó population, which is, in its genetic structure, clearly different from the Székely population. The pronounced genetic divergence between the two populations is in sharp contrast to the expectation of high genetic similarity due to the close geographic proximity of their native homelands. The population data will be incorporated in the EMPOP database (http://www.empop.org).

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.

Development of species identification tests targeting the 16S ribosomal RNA coding region in mitochondrial DNA

This report describes the development of a species testing system based on the diversity of nucleotide sequences in mitochondrial DNA (mtDNA) among species. Five species, human, cow, pig, dog, and cat, were considered. The partial nucleotide sequences in 16S ribosomal RNA coding region were chosen as the target for discriminating the species. The sequence diversities of this approximately 400 bp long region ranged from 15.7 to 24.1% among the five species. Sequencing analysis of this target on 50 individuals of each species (53 for dogs) revealed that the nucleotide sequences were well preserved within species. Species-specific PCR for each species was also designed, and satisfactory results with regard to both sensitivity and specificity were obtained. A validation study with DNA extracted from bovine bone exposed to the environment revealed that the PCRs designed in this study worked correctly. From the results obtained, it was shown that this testing system could be a good tool for species identification. One successful case report is also demonstrated.

Application of a quasi-median network analysis for the visualization of character conflicts to a population sample of mitochondrial DNA control region sequences from southern Germany (Ulm)

Entire mtDNA control region sequences from 100 individuals in a west Eurasian population sample from southern Germany (around the city of Ulm) were generated and analyzed. The control region was amplified in one piece and sequenced with ten different sequencing primers. Sequence evaluation was performed independently. Phylogenetic analyses were used for quality assurance purposes and for the determination of the haplogroup affiliation of the samples. The sequences were scrutinized performing a quasi-median network analysis. To visualize character conflicts, frequent mutations were filtered, and the reduced data were represented by the torso of their quasi-median network. Character incompatibilities were found to be based on real biological patterns of homoplasy. The population data will be incorporated in the EMPOP database ( http://www.empop.org ).

Updating the East Asian mtDNA phylogeny: a prerequisite for the identification of pathogenic mutations

Knowledge about the world phylogeny of human mitochondrial DNA (mtDNA) is essential not only for evaluating the pathogenic role of specific mtDNA mutations but also for performing reliable association studies between mtDNA haplogroups and complex disorders. In the past few years, the main features of the East Asian portion of the mtDNA phylogeny have been determined on the basis of complete sequencing efforts, but representatives of several basal lineages were still lacking. Moreover, some recently published complete mtDNA sequences did apparently not fit into the known phylogenetic tree and conflicted with the established nomenclature. To refine the East Asian mtDNA tree and resolve data conflicts, we first completely sequenced 20 carefully selected mtDNAs--likely representatives of novel sub-haplogroups--and then, in order to distinguish diagnostic mutations of novel haplogroups from private variants, we applied a 'motif-search' procedure to a large sample collection. The novel information was incorporated into an updated East Asian mtDNA tree encompassing more than 1000 (near-) complete mtDNA sequences. A reassessment of the mtDNA data from a series of disease studies testified to the usefulness of such a refined mtDNA tree in evaluating the pathogenicity of mtDNA mutations. In particular, the claimed pathogenic role of mutations G3316A, T3394C, A4833G and G15497A appears to be most questionable as those initial claims were derived from anecdotal findings rather than e.g. appropriate association studies. Following a guideline based on the phylogenetic knowledge as proposed here could help avoiding similar problems in the future.

Forensic genetic analysis of mitochondrial DNA hypervariable region I/II sequences: An expanded Korean population database

We have analyzed variation of the mitochondrial DNA (mtDNA) hypervariable segments I and II (HVS-I and HVS-II) in 185 randomly chosen individuals from Korea to provide an expanded and reliable Korean database. Combined sequence comparison of HVS-I and HVS-II led to the identification of 167 different haplotypes characterized by 154 variable sites. One hundred and fifty-one of the haplotypes were individual-specific, 14 were found in two individuals and 2 were found in three individuals. A pairwise comparison of the 185 HVS-I/II sequences found an average of 10.11 +/- 4.63 differences between individuals. The random match probability and gene diversity for the combined hypervariable regions were estimated at 0.66% and 0.9988, respectively. Analyzing the expanded database including three previously reported data sets and the present data using haplogroup-based comparisons and comparison with closely related sequences allowed errors to be detected and eliminated, thus considerably improving data quality. Sample division comparisons based on PhiST genetic distance measures revealed no significant population differentiation in the distribution of mtDNA sequence variations between the present data set and a database in The Scientific Working Group on DNA Analysis Methods (SWGDAM), but did indicate differences from other sets of data. Based on the results of mtDNA profiles, almost all of the mtDNA types studied here could be classified into subsets of haplogroups common in east Asia, and show that the Koreans possess lineages from both the southern and the northern haplogroup complexes of east Asian populations. The new data, combined with other mtDNA sequences, demonstrate how useful comparison with closely related mtDNA sequences can be for improving database quality, as well as providing haplotype information for forensic and population genetic analyses in the Korean population.

Sequence polymorphisms of the mtDNA control region in a human isolate: the Georgians from Swanetia

In this work, we analyzed the sequence diversity of the mtDNA control region (HVI and HVII) in a sample of 48 individuals from Swanetia (Georgia), using direct fluorescent-based sequencing methods. We identified 43 different mtDNA haplotypes resulting from 78 polymorphic sites (46 in HVI and 32 in HVII). Most of the variable positions identified in both HVI and HVII were transitions (82.6 and 71.9%, respectively). The frequency of length heteroplasmy in the homopolymeric C-stretch regions was the same for both segments (10.4%). The sequence diversity increased markedly when both hypervariable regions were analyzed jointly (HVI: 0.985, HVII: 0.975, HVI+HVII: 0.994). Accordingly, the probability of two randomly selected sequences matching (random match probability, RMP) decreased from 3.4% (HVI) to 2.6% (HVI+HVII), despite which the RMP values in Georgians remained higher than estimated in most Europeans. This suggests that the variability of maternal lineages tends to be lower in traditional human isolates and, therefore, the potential of discrimination of mtDNA in forensic analysis is more limited in this type of population. The incorporation of HVII data also contributed to the refinement of results regarding the genetic relationships among the samples included in the analyses, which stress the importance of considering HVII in both population and forensic genetics.

A reappraisal of complete mtDNA variation in East Asian families with hearing impairment

In a number of recent studies, we summarized the obvious errors and shortcomings that can be spotted in many (if not most) mitochondrial DNA (mtDNA) data sets published in medical genetics. We have reanalyzed here the complete mtDNA genome data published in various recent reports of East Asian families with hearing impairment, using a phylogenetic approach, in order to demonstrate the persistence of lab-specific mistakes in mtDNA genome sequencing in cases where those caveats were (deliberately) neglected. A phylogenetic reappraisal of complete mtDNAs with mutation A1555G (or G11778A) indeed supports the suggested lack of association between haplogroup background and phenotypic presentation of these mutations in East Asians. In contrast, the claimed pathogenicity of mutation T1095C in Chinese families with hearing impairment seems unsupported, basically because this mutation is rather basal in the mtDNA phylogeny, being specific to haplogroup M11 in East Asia. The roles of other haplogroup specific or associated variants, such as A827G, T961C, T1005C, in East Asian subjects with aminoglycoside-induced and non-syndromic hearing loss are also unclear in view of the known mtDNA phylogeny.

Multistep microsatellite mutation in the maternally transmitted locus D13S317: a case of maternal allele mismatch in the child

Examination of a case of a paternity dispute with 17 autosomal short tandem repeats (STR) loci revealed a mismatch of the maternally transmitted allele at the locus D13S317 in the questioned child. The composition of the alleles of this locus in the mother, questioned child and suspected father was 8/8, 11/11 and 8/11, respectively. The sequence analysis of the regions flanking the locus D13S317 and peak height measurements of the paternal, maternal and child alleles at this locus excluded the possibility of null allele as a cause of the allelic mismatch inherited by the child. The results suggested expansion of the microsatellite repeat motif, TATC by three repeat units as a probable cause for the allelic mismatch in the child. This is a rare case of maternally transmitted multistep microsatellite mutation reported for the first time for this locus in the forensic DNA analysis. The mutation rate at D13S317 locus in maternal and paternal meiosis was 0.04 and 0.14%, respectively, and overall mutation rate was 0.15%. The probability of maternity and paternity were 0.999999 and 0.999999, respectively, for all the 17 autosomal STR loci analyzed. Furthermore, the sequence of two hypervariable regions of mitochondrial DNA, HV1 and HV2 and the maternal alleles of six X chromosome STR loci in the questioned child matched completely with the mother. These results conclusively proved that the mother and suspected father are the biological parents of the questioned child.

Ancient mitochondrial DNA from Malaysian hair samples: Some indications of Southeast Asian population movements

The late Pleistocene and early Holocene population history of Southeast Asia is not well-known. Our study provides new data on mitochondrial DNA (mtDNA) lineages of the aboriginal inhabitants of the Malay Peninsula, and through an extensive comparison to the known mtDNA diversity in Southeast and East Asia, provides some new insights into the origins and historical geography of certain mtDNA lineages in the region. We extracted DNA from hair samples (dating back 100 years) preserved in the Duckworth Collection and belonging to two Peninsular Malaysian individuals identified as "Negrito." Ancient DNA was analyzed by sequencing hypervariable region I (HVS-I) of the mtDNA control region and the mtDNA region V length polymorphism. The results show that the maternal lineages of these individuals belong to a recently defined haplogroup B sub-branch called B4c2. A comparison of mitochondrial haplotypes and haplogroups with those of 10,349 East Asian individuals indicates their very restricted geographical distribution (southwestern China, Southeast Asia Peninsula, and Indonesia). Recalculation of the B4c2 age across all of East Asia ( approximately 13,000 years) and in different subregions/populations suggests its rapid diffusion in Southeast Asia between the end of the Last Glacial Maximum and the Neolithic expansion of the Holocene.

Hypervariable region structure and polymorphism of mtDNA from dental pulp and a family analysis

Nucleotide sequences of the hypervariable region in the D-loop of mitochondrial DNA (mtDNA) were analyzed using DNA extracted from 140 old dental pulp samples. These sequences were compared with the sequence reported by Anderson et al. Nucleotide substitution in the HV1 region was identified at 77 positions. A C-to-T transition at position 16223 (C16223T) was most frequently detected (77.9%). Fourteen types of C-stretch sequence patterns were detected and the same sequence as Anderson had the highest frequency (57.9%). In the HV2 region, base transitions were identified at 56 positions. A263G was identified in all samples. Seven types of C-stretch were detected, but none had the same sequence as Anderson. In the HV3 region, base transitions were identified at 21 positions. T489C was most frequently identified (64.3%). Five types of C-stretch were detected, and the same sequence as Anderson accounted for 92.9%. The 140 samples were classified into 128 kinds by the sequence patterns of the HV region. Next, using the blood and oral mucosa epithelium from 23 subjects comprising four generations in a family line, the hereditary relationship of mtDNA was examined. All mtDNA types of the first-generation mother were infallibly inherited by the fourth generation.

A critical reassessment of the role of mitochondria in tumorigenesis

BACKGROUND

Mitochondrial DNA (mtDNA) is being analyzed by an increasing number of laboratories in order to investigate its potential role as an active marker of tumorigenesis in various types of cancer. Here we question the conclusions drawn in most of these investigations, especially those published in high-rank cancer research journals, under the evidence that a significant number of these medical mtDNA studies are based on obviously flawed sequencing results.

METHODS AND FINDINGS

In our analyses, we take a phylogenetic approach and employ thorough database searches, which together have proven successful for detecting erroneous sequences in the fields of human population genetics and forensics. Apart from conceptual problems concerning the interpretation of mtDNA variation in tumorigenesis, in most cases, blocks of seemingly somatic mutations clearly point to contamination or sample mix-up and, therefore, have nothing to do with tumorigenesis.

CONCLUSION

The role of mitochondria in tumorigenesis remains unclarified. Our findings of laboratory errors in many contributions would represent only the tip of the iceberg since most published studies do not provide the raw sequence data for inspection, thus hindering a posteriori evaluation of the results. There is no precedent for such a concatenation of errors and misconceptions affecting a whole subfield of medical research.

Mosaics of ancient mitochondrial DNA: positive indicators of nonauthenticity

Research into ancient mitochondrial DNA is plagued by contamination, post mortem damage, and other artefacts. The stringent set of controls suggested by Cooper and Poinar a few years ago are, however, rarely followed in practice, and even when applied carefully, these criteria need not be sufficient to guarantee authenticity. The fairly relaxed prerequisites now common for ancient population studies have opened the door for all kinds of contamination and sequencing errors to enter ancient mtDNA data. To reject or question authenticity of particular sequencing results a posteriori, one can follow similar strategies of focused database comparisons that have proven to be effective and successful in the case of flawed modern mtDNA data.

Recent progress in mitochondrial DNA analysis

In this review, we describe the current state of knowledge of mitochondrial genetics of East Asian populations and its application to forensic science. Recent advances in mitochondrial DNA (mtDNA) phylogeny have identified haplogroup-specific single nucleotide polymorphisms (SNPs) and the control region motifs of haplogroups. By analyzing haplogroup-specific SNPs, we can rapidly and accurately connect the mtDNA under study to the relevant haplogroup. Haplogroups are fairly continent- and/or region-specific; therefore, we can infer the ethnic background of that mtDNA. In addition, errors in hypervariable region sequences can be detected by means of haplogroup motif analysis.

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 nt303-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.

Artificial recombination in forensic mtDNA population databases

Artificial recombination of two or more mitochondrial DNA fragments from different samples would constitute a serious cause of error in forensic DNA typing, and yet one can demonstrate that such events have happened in the preparation of several published mtDNA databases. Focussed database searches, phylogenetic analysis, and network representations can highlight mosaic patterns and thus pinpoint sample mix-up. Therefore, we suggest that this approach should be applied to data prior to publication in order to uncover such errors in time.

Mitochondrial DNA control region sequences from Nairobi (Kenya): inferring phylogenetic parameters for the establishment of a forensic database

Large forensic mtDNA databases which adhere to strict guidelines for generation and maintenance, are not available for many populations outside of the United States and western Europe. We have established a high quality mtDNA control region sequence database for urban Nairobi as both a reference database for forensic investigations, and as a tool to examine the genetic variation of Kenyan sequences in the context of known African variation. The Nairobi sequences exhibited high variation and a low random match probability, indicating utility for forensic testing. Haplogroup identification and frequencies were compared with those reported from other published studies on African, or African-origin populations from Mozambique, Sierra Leone, and the United States, and suggest significant differences in the mtDNA compositions of the various populations. The quality of the sequence data in our study was investigated and supported using phylogenetic measures. Our data demonstrate the diversity and distinctiveness of African populations, and underline the importance of establishing additional forensic mtDNA databases of indigenous African populations.

DNA typing of bone specimens—the potential use of the profiler test as a tool for bone identification

Twenty-six bone DNA identification cases are described. The postmortem periods of the studied remains ranged from three days to over 30 years, and the locations where the remains were found varied resulting in a variety of postmortem conditions. Nuclear DNA typing using an AmpFLSTR Profiler kit and mitochondrial DNA (mtDNA) typing of hypervariable regions 1 and 2 (HV1 and HV2) in a control region were performed both with decalcified and non-treated bone powder samples. Decalcification was shown to improve the success of DNA typing. The nucleotide sequences of the HV1 and HV2 regions were successfully determined in all cases examined. Nuclear DNA typing was very successful, more than half of the loci were typed during multiple amplifications (10 loci in one reaction) in 23 cases. Polymerase chain reaction (PCR) inhibition was observed in five cases including three samples that were found buried in soil. This inhibitory effect was identified as the result of unbalanced multiple PCR during the profiler test. These results revealed that DNA typing targeting nuclear DNA is a potentially powerful tool for bone identification.

Identification of Unknown Body Using DNA Analysis and Dental Characteristics in Chest X-ray Photograph

An unknown skeletonized body was identified by DNA analysis and dental information. The body had already been cremated when a candidate for the unknown body was proposed. Therefore, for DNA analysis we used teeth that had been kept for a long time after use for serological examination. We also used a chest X-ray photograph of the candidate and photographs of dentition, as well as dental X-ray photographs taken when the unknown body was found. Because DNA obtained from teeth was highly degraded, we amplified three PCR fragments to determine the 766 bp mitochondrial DNA (mtDNA) sequence including HV1 and HV2. Polymorphism of the ABO locus was also analyzed using small PCR fragments. Although the isolated DNA was contaminated, probably with DNA from a different individual, DNA polymorphisms of mtDNA and the ABO locus could be analyzed. We discuss the reliability of our conclusions from the point of view of the necessity of constructing an accurate mtDNA database. Although a dentist who had treated the teeth of the unknown body could not be found, a chest X-ray photograph for medical diagnosis was very useful in comparing dental characteristics, as it included an image of the frontal part of the lower jaw and upper teeth.

Mitochondrial genome variation in eastern Asia and the peopling of Japan

To construct an East Asia mitochondrial DNA (mtDNA) phylogeny, we sequenced the complete mitochondrial genomes of 672 Japanese individuals (http://www.giib.or.jp/mtsnp/index_e.html). This allowed us to perform a phylogenetic analysis with a pool of 942 Asiatic sequences. New clades and subclades emerged from the Japanese data. On the basis of this unequivocal phylogeny, we classified 4713 Asian partial mitochondrial sequences, with <10% ambiguity. Applying population and phylogeographic methods, we used these sequences to shed light on the controversial issue of the peopling of Japan. Population-based comparisons confirmed that present-day Japanese have their closest genetic affinity to northern Asian populations, especially to Koreans, which finding is congruent with the proposed Continental gene flow to Japan after the Yayoi period. This phylogeographic approach unraveled a high degree of differentiation in Paleolithic Japanese. Ancient southern and northern migrations were detected based on the existence of basic M and N lineages in Ryukyuans and Ainu. Direct connections with Tibet, parallel to those found for the Y-chromosome, were also apparent. Furthermore, the highest diversity found in Japan for some derived clades suggests that Japan could be included in an area of migratory expansion to Continental Asia. All the theories that have been proposed up to now to explain the peopling of Japan seem insufficient to accommodate fully this complex picture.

Identification of Native American founder mtDNAs through the analysis of complete mtDNA sequences: some caveats

In this study, a detailed analysis of both previously published and new data was performed to determine whether complete, or almost complete, mtDNA sequences can resolve the long-debated issue of which Asian mtDNAs were founder sequences for the Native American mtDNA pool. Unfortunately, we now know that coding region data and their analysis are not without problems. To obtain and report reasonably correct sequences does not seem to be a trivial task, and to discriminate between Asian and Native American mtDNA ancestries may be more complex than previously believed. It is essential to take into account the effects of mutational hot spots in both the control and coding regions, so that the number of apparent Native American mtDNA founder sequences is not erroneously inflated. As we report here, a careful analysis of all available data indicates that there is very little evidence that more than five founder mtDNA sequences entered Beringia before the Last Glacial Maximum and left their traces in the current Native American mtDNA pool.

A call for mtDNA data quality control in forensic science

There is increasing evidence that many of the mitochondrial DNA (mtDNA) databases published in the fields of forensic science and molecular anthropology are flawed. An a posteriori phylogenetic analysis of the sequences could help to eliminate most of the errors and thus greatly improve data quality. However, previously published caveats and recommendations along these lines were not yet picked up by all researchers. Here we call for stringent quality control of mtDNA data by haplogroup-directed database comparisons. We take some problematic databases of East Asian mtDNAs, published in the Journal of Forensic Sciences and Forensic Science International, as examples to demonstrate the process of pinpointing obvious errors. Our results show that data sets are not only notoriously plagued by base shifts and artificial recombination but also by lab-specific phantom mutations, especially in the second hypervariable region (HVR-II).

Results of a collaborative study of the EDNAP group regarding mitochondrial DNA heteroplasmy and segregation in hair shafts

A collaborative exercise was carried out by the European DNA Profiling Group (EDNAP) in order to evaluate the distribution of mitochondrial DNA (mtDNA) heteroplasmy amongst the hairs of an individual who displays point heteroplasmy in blood and buccal cells. A second aim of the exercise was to study reproducibility of mtDNA sequencing of hairs between laboratories using differing chemistries, further to the first mtDNA reproducibility study carried out by the EDNAP group. Laboratories were asked to type 2 sections from each of 10 hairs, such that each hair was typed by at least two laboratories. Ten laboratories participated in the study, and a total of 55 hairs were typed. The results showed that the C/T point heteroplasmy observed in blood and buccal cells at position 16234 segregated differentially between hairs, such that some hairs showed only C, others only T and the remainder, C/T heteroplasmy at varying ratios. Additionally, differential segregation of heteroplasmic variants was confirmed in independent extracts at positions 16093 and the poly(C) tract at 302-309, whilst a complete A-G transition was confirmed at position 16129 in one hair. Heteroplasmy was observed at position 16195 on both strands of a single extract from one hair segment, but was not observed in the extracts from any other segment of the same hair. Similarly, heteroplasmy at position 16304 was observed on both strands of a single extract from one hair. Additional variants at positions 73, 249 and the HVII poly(C) region were reported by one laboratory; as these were not confirmed in independent extracts, the possibility of contamination cannot be excluded. Additionally, the electrophoresis and detection equipment used by this laboratory was different to those of the other laboratories, and the discrepancies at position 249 and the HVII poly(C) region appear to be due to reading errors that may be associated with this technology. The results, and their implications for forensic mtDNA typing, are discussed in the light of the biology of hair formation.

Mitochondrial DNA D-loop hypervariable regions: Czech population data

In order to identify polymorphic sites and to find out their frequencies and the frequency of haplotypes, the complete D-loop of mitochondrial DNA (mtDNA) from 93 unrelated Czech Caucasians was sequenced. Sequence comparison showed that 85 haplotypes were found and of these 78 were unique, 6 were observed twice and 1 was observed three times. Genetic diversity (GD) was estimated at 0.999 and the probability of two randomly selected sequences matching (random match probability, RMP) at 1.2%. Additionally these calculations were carried out for hypervariable regions 1, 2 (HV1, HV2), for the area between HV1 and HV2 and for the area of the hypervariable region HV3. The average number of nucleotide differences (ANND) was established to be 10.2 for the complete D-loop. The majority of sequence variations were substitutions, particularly transitions. Deletions were found only in the region where HV3 is situated and insertions in the same place and in poly-C tracts between positions 303 and 315 in HV2. A high degree of length heteroplasmy was found especially in the regions of poly-C tracts between positions 16184 and 16193 in HV1 and between positions 303 and 315 in HV2. Position heteroplasmies were found in two cases.

Sequence analysis of two de novo mutation alleles at the DXS10011 locus

We have detected two unusual alleles at the DXS10011 locus in two paternity trio cases. In one case, one allele of the daughter was found not to have been derived from the mother but the other allele was shared with the father. In the other case, the mother and the son shared no bands. Paternity in both cases was established using conventional polymorphic markers in addition to DNA markers (probabilities: >0.999999). Sequencing showed that the two de novo alleles of the children acquired a single unit (GAAA).

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.

Genetic analysis and ethnic affinities from two Scytho-Siberian skeletons

We extracted DNA from two skeletons belonging to the Sytho-Siberian population, which were excavated from the Sebÿstei site (dating back 2,500 years) in the Altai Republic (Central Asia). Ancient DNA was analyzed by autosomal short tandem repeats (STRs) and by the sequencing of the hypervariable region 1 (HV1) of the mitochondrial DNA (mtDNA) control region. The results showed that these two skeletons were not close relatives. Moreover, their haplogroups were characteristic of Asian populations. Comparison with the haplogroup of 3,523 Asian and American individuals linked one skeleton with a putative ancestral paleo-Asiatic population and the other with Chinese populations. It appears that the genetic study of ancient populations of Central Asia brings important elements to the understanding of human population movements in Asia.