Reduced-median-network analysis of complete mitochondrial DNA coding-region sequences for the major African, Asian, and European haplogroups

The evolution of the human mitochondrial genome is characterized by the emergence of ethnically distinct lineages or haplogroups. Nine European, seven Asian (including Native American), and three African mitochondrial DNA (mtDNA) haplogroups have been identified previously on the basis of the presence or absence of a relatively small number of restriction-enzyme recognition sites or on the basis of nucleotide sequences of the D-loop region. We have used reduced-median-network approaches to analyze 560 complete European, Asian, and African mtDNA coding-region sequences from unrelated individuals to develop a more complete understanding of sequence diversity both within and between haplogroups. A total of 497 haplogroup-associated polymorphisms were identified, 323 (65%) of which were associated with one haplogroup and 174 (35%) of which were associated with two or more haplogroups. Approximately one-half of these polymorphisms are reported for the first time here. Our results confirm and substantially extend the phylogenetic relationships among mitochondrial genomes described elsewhere from the major human ethnic groups. Another important result is that there were numerous instances both of parallel mutations at the same site and of reversion (i.e., homoplasy). It is likely that homoplasy in the coding region will confound evolutionary analysis of small sequence sets. By a linkage-disequilibrium approach, additional evidence for the absence of human mtDNA recombination is presented here.     

Phylogeographic differentiation of mitochondrial DNA in Han Chinese

To characterize the mitochondrial DNA (mtDNA) variation in Han Chinese from several provinces of China, we have sequenced the two hypervariable segments of the control region and the segment spanning nucleotide positions 10171-10659 of the coding region, and we have identified a number of specific coding-region mutations by direct sequencing or restriction-fragment-length-polymorphism tests. This allows us to define new haplogroups (clades of the mtDNA phylogeny) and to dissect the Han mtDNA pool on a phylogenetic basis, which is a prerequisite for any fine-grained phylogeographic analysis, the interpretation of ancient mtDNA, or future complete mtDNA sequencing efforts. Some of the haplogroups under study differ considerably in frequencies across different provinces. The southernmost provinces show more pronounced contrasts in their regional Han mtDNA pools than the central and northern provinces. These and other features of the geographical distribution of the mtDNA haplogroups observed in the Han Chinese make an initial Paleolithic colonization from south to north plausible but would suggest subsequent migration events in China that mainly proceeded from north to south and east to west. Lumping together all regional Han mtDNA pools into one fictive general mtDNA pool or choosing one or two regional Han populations to represent all Han Chinese is inappropriate for prehistoric considerations as well as for forensic purposes or medical disease studies.     

Maternal footprints of Southeast Asians in North India

We have analyzed 7,137 samples from 125 different caste, tribal and religious groups of India and 99 samples from three populations of Nepal for the length variation in the COII/tRNA(Lys) region of mtDNA. Samples showing length variation were subjected to detailed phylogenetic analysis based on HVS-I and informative coding region sequence variation. The overall frequencies of the 9-bp deletion and insertion variants in South Asia were 1.9 and 0.6%, respectively. We have also defined a novel deep-rooting haplogroup M43 and identified the rare haplogroup H14 in Indian populations carrying the 9-bp deletion by complete mtDNA sequencing. Moreover, we redefined haplogroup M6 and dissected it into two well-defined subclades. The presence of haplogroups F1 and B5a in Uttar Pradesh suggests minor maternal contribution from Southeast Asia to Northern India. The occurrence of haplogroup F1 in the Nepalese sample implies that Nepal might have served as a bridge for the flow of eastern lineages to India. The presence of R6 in the Nepalese, on the other hand, suggests that the gene flow between India and Nepal has been reciprocal.     

Phylogenetic network of the mtDNA haplogroup U in Northern Finland based on sequence analysis of the complete coding region by conformation-sensitive gel electrophoresis

Mutations in mtDNA have accumulated sequentially, and maternal lineages have diverged to form population-specific genotypes. Classification of the genotypes has been made based on differences found in restriction fragment analysis of the coding region or in the sequence of the hypervariable segment I. Both methods have shortcomings, as the former may not detect all the important polymorphisms and the latter makes use of a segment containing hypervariable nucleotide positions. Here, we have used conformation-sensitive gel electrophoresis (CSGE) to detect polymorphisms within the coding region of mtDNA from 22 Finns belonging to haplogroup U. Sixty-three overlapping PCR fragments covering the entire coding region were analyzed by CSGE, and the fragments that differed in their migration pattern were sequenced. CSGE proved to be a sensitive and specific method for identifying mtDNA substitutions. The phylogenetic network of the 22 coding-region sequences constituted a perfect tree, free of homoplasy, and provided several previously unidentified common polymorphisms characterizing subgroups of U. After contrasting this data with that of hypervariable segment I, we concluded that position 16192 seems to be prone to recurrent mutations and that position 16270 has experienced a back mutation. Interestingly, all 22 samples were found to belong to subcluster U5, suggesting that this subcluster is more frequent in Finns than in other European populations. Complete sequence data of the mtDNA yield a more reliable phylogenetic network and a more accurate classification of the haplogroups than previous ones. In medical genetics, such networks may help to decide between a rare polymorphism and a pathogenic mutation; in population genetics, the networks may enable more detailed analyses of population history and mtDNA evolution.     

Phylogeny of mitochondrial DNA macrohaplogroup N in India, based on complete sequencing: implications for the peopling of South Asia

To resolve the phylogeny of the autochthonous mitochondrial DNA (mtDNA) haplogroups of India and determine the relationship between the Indian and western Eurasian mtDNA pools more precisely, a diverse subset of 75 macrohaplogroup N lineages was chosen for complete sequencing from a collection of >800 control-region sequences sampled across India. We identified five new autochthonous haplogroups (R7, R8, R30, R31, and N5) and fully characterized the autochthonous haplogroups (R5, R6, N1d, U2a, U2b, and U2c) that were previously described only by first hypervariable segment (HVS-I) sequencing and coding-region restriction-fragment-length polymorphism analysis. Our findings demonstrate that the Indian mtDNA pool, even when restricted to macrohaplogroup N, harbors at least as many deepest-branching lineages as the western Eurasian mtDNA pool. Moreover, the distribution of the earliest branches within haplogroups M, N, and R across Eurasia and Oceania provides additional evidence for a three-founder-mtDNA scenario and a single migration route out of Africa.     

High resolution analysis and phylogenetic network construction using complete mtDNA sequences in sardinian genetic isolates

For mitochondrial phylogenetic analysis, the best result comes from complete sequences. We therefore decided to sequence the entire mitochondrial DNA (mtDNA) (coding and D-loop regions) of 63 individuals selected in 3 small Ogliastra villages, an isolated area of eastern Sardinia: Talana, Urzulei, and Perdasdefogu. We studied at least one individual for each of the most frequent maternal genealogical lineages belonging to haplogroups H, V, J, K, T, U, and X. We found in our 63 samples, 172 and 69 sequence changes in the coding and in the D-loop region, respectively. Thirteen out of 172 sequence changes in the coding region are novel. It is our hypothesis that some of them are characteristic of the Ogliastra region and/or Sardinia. We reconstructed the phylogenetic network of the 63 complete mtDNA sequences for the 3 villages. We also drew a network including a large number of European sequences and calculated various indices of genetic diversity in Ogliastra. It appears that these small populations remained extremely isolated and genetically differentiated compared with other European populations. We also identified in our samples a never previously described subhaplogroup, U5b3, which seems peculiar to the Ogliastra region.     

人类线粒体DNA变异的检测方法和思路

基于线粒体DNA（mtDNA)的研究对于人群源流迁移、线粒体相关疾病病因的探讨和法医鉴定等具有重意义,就检测人线粒体突变的一些常用方法,如RFLP、SSO和控制区测序等作一小结和归纳,并重点介绍目前mtDNA突变的筛选方法和思路,另外,还总结了近年来对人mtDNA方面的研究结果,对世界人群中主要单倍型类群（haplogroup)特征变异位点和相应的酶切检测引物作了归纳。     

Mitochondrial DNA variation in two South Siberian Aboriginal populations: implications for the genetic history of North Asia

The mtDNAs of 76 individuals representing the aboriginal populations of South Siberia, the Tuvinians and Buryats, were subjected to restriction fragment length polymorphism (RFLP) analysis and control region hypervariable segment I (HVS-I) sequencing, and the resulting data were combined with those available for other Siberian and East Asian populations and subjected to statistical and phylogenetic analysis. This analysis showed that the majority of the Tuvinian and Buryat mtDNAs (94.4% and 92.5%, respectively) belong to haplogroups A, B, C, D, E, F, and M*, which are characteristic of Mongoloid populations. Furthermore, the Tuvinians and Buryats harbor four Asian- and Native American-specific haplogroups (A-D) with frequencies (72.2% and 55%, respectively) exceeding those reported previously for Mongolians, Chinese, and Tibetans. They represent, therefore, the populations that are most closely related to New World indigenous groups. Despite their geographical proximity, the Tuvinians and Buryats shared no HVS-I sequences in common, although individually they shared such sequences with a variety of other Siberian and East Asian populations. In addition, phylogenetic and principal component analyses data of mtDNA sequences show that the Tuvinians clustered more closely with Turkic-speaking Yakuts, whereas the Mongolic-speaking Buryats clustered closer to Korean populations. Furthermore, HVS-I sequences, comprising one-fourth of the Buryat lineages and characterized by the only C-to-T transition at nucleotide position 16223, were identified as different RFLP haplotypes (B, C, D, E, M*, and H). This finding appears to indicate the putative ancestral state of the 16223T HVS-I sequences to Mongoloid macrohaplogroup M, at least. Finally, the results of nucleotide diversity analysis in East Asian and Siberian populations suggest that Central and East Asia were the source areas from which the genetically heterogeneous Tuvinians and Buryats first emerged.     

The complete nucleotide sequence of the mitochondrial genome of Tetraodon nigroviridis.

The fresh water pufferfish Tetraodon nigroviridis is a model organism for studying evolution of genome and gene functions, but its mitochondrial genome (mtDNA) sequence is still not available. We determined the complete nucleotide sequence of its mtDNA using shotgun sequencing. The T. nigroviridis mtDNA was 16,462 bp, and contained 13 protein coding genes, 22 tRNAs, 2 rRNAs and a major non-coding region. The gene order was identical to the common type of vertebrate mtDNA, whereas the G + C content in the sense strand was 46.9%, much higher than most other fish species. One hundred and three SNPs were detected in the control region of the mtDNA of 35 individuals, a majority of SNPs were detected in the 5' end of the control region. A phylogenetic study including 21 fish species was performed on concatenated amino acid sequences of 12 protein coding genes, and revealed that the T. nigroviridis was clustered with Fugu rubripes into a group. The complete mtDNA sequence and SNPs in its control region will be useful in studying fish evolution, in differentiating different Tetraodon species and in analyzing genetic diversity within T. nigroviridis.     

Mitochondrial DNA lineages of elite Ethiopian athletes

Previous studies have hypothesised that mitochondrial DNA (mtDNA) polymorphisms may influence aerobic performance. The matrilineal inheritance and accumulation of polymorphisms in mtDNA means that mtDNA haplogroups, characterised by key polymorphisms, are often represented at different frequencies in different populations. The present study aimed to compare the mtDNA haplogroup distribution of elite Ethiopian athletes relative to the general Ethiopian population. The haplogroup distribution of 76 endurance athletes (E), members of the Ethiopian national athletics team, was compared to 108 members of the general Ethiopian population (C). DNA was extracted from buccal swabs and haplogroups assigned by sequencing part of the hypervariable sequence (HVS-I), followed by analysis of key coding-region polymorphisms. A high proportion of African 'L' haplogroups was found in athletes and controls (C=53%; E=55%). Haplogroup distribution of endurance runners did not differ from that of C (P=0.63). Elite Ethiopian athletes are not a mitochondrially distinct group relative to the Ethiopian population. It appears that environment and, perhaps, polymorphisms in the nuclear genome are more important determinants of Ethiopian running success than mtDNA polymorphisms.     

The emerging limbs and twigs of the East Asian mtDNA tree

We determine the phylogenetic backbone of the East Asian mtDNA tree by using published complete mtDNA sequences and assessing both coding and control region variation in 69 Han individuals from southern China. This approach assists in the interpretation of published mtDNA data on East Asians based on either control region sequencing or restriction fragment length polymorphism (RFLP) typing. Our results confirm that the East Asian mtDNA pool is locally region-specific and completely covered by the two superhaplogroups M and N. The phylogenetic partitioning based on complete mtDNA sequences corroborates existing RFLP-based classification of Asian mtDNA types and supports the distinction between northern and southern populations. We describe new haplogroups M7, M8, M9, N9, and R9 and demonstrate by way of example that hierarchically subdividing the major branches of the mtDNA tree aids in recognizing the settlement processes of any particular region in appropriate time scale. This is illustrated by the characteristically southern distribution of haplogroup M7 in East Asia, whereas its daughter-groups, M7a and M7b2, specific for Japanese and Korean populations, testify to a presumably (pre-)Jomon contribution to the modern mtDNA pool of Japan.     

Mitochondrial D-Loop Variation in Persian Multiple Sclerosis Patients: K and A Haplogroups as a Risk Factor!!

Multiple Sclerosis (MS) is a multifocal demyelinating central nervous system disorder in which interplay between genes and the environment are supposed to be involved. Mitochondrial DNA (mtDNA) has the only non-coding regions at the displacement loop (D-loop) region that contains two hypervariable segments (HVS-I and HVS-II) with high polymorphism. mtDNA has already been fully sequenced and many subsequent publications have showed polymorphic sites, haplogroups and haplotypes. Haplogroups could have important implications to understand association between mutability of the mitochondrial genome and disease. To assess relationship between mtDNA haplogroups and MS, we have sequenced the mtDNA HVS-I in 54 MS patients and 100 control subjects. We have found that haplogroups A and K are significantly more abundant in MS patients (P=0.042 for haplogroup A and P=0.0005 for haplogroup K). Thus, these two haplogroups might act synergistically to increase the penetrance of MS disease.     

The Genographic Project public participation mitochondrial DNA database

The Genographic Project is studying the genetic signatures of ancient human migrations and creating an open-source research database. It allows members of the public to participate in a real-time anthropological genetics study by submitting personal samples for analysis and donating the genetic results to the database. We report our experience from the first 18 months of public participation in the Genographic Project, during which we have created the largest standardized human mitochondrial DNA (mtDNA) database ever collected, comprising 78,590 genotypes. Here, we detail our genotyping and quality assurance protocols including direct sequencing of the mtDNA HVS-I, genotyping of 22 coding-region SNPs, and a series of computational quality checks based on phylogenetic principles. This database is very informative with respect to mtDNA phylogeny and mutational dynamics, and its size allows us to develop a nearest neighbor-based methodology for mtDNA haplogroup prediction based on HVS-I motifs that is superior to classic rule-based approaches. We make available to the scientific community and general public two new resources: a periodically updated database comprising all data donated by participants, and the nearest neighbor haplogroup prediction tool.     

Biological relationship between Central and South American Chibchan speaking populations: evidence from mtDNA

We examined mitochondrial DNA (mtDNA) haplogroup and haplotype diversity in 188 individuals from three Chibchan (Kogi, Arsario, and Ijka) populations and one Arawak (Wayuú) group from northeast Colombia to determine the biological relationship between lower Central American and northern South American Chibchan speakers. mtDNA haplogroups were obtained for all individuals and mtDNA HVS-I sequence data were obtained for 110 samples. Resulting sequence data were compared to 16 other Caribbean, South, and Central American populations using diversity measures, neutrality test statistics, sudden and spatial mismatch models, intermatch distributions, phylogenetic networks, and a multidimensional scaling plot. Our results demonstrate the existence of a shared maternal genetic structure between Central American Chibchan, Mayan populations and northern South American Chibchan-speakers. Additionally, these results suggest an expansion of Chibchan-speakers into South America associated with a shift in subsistence strategies because of changing ecological conditions that occurred in the region between 10,000-14,000 years before present.     

Phylogenetic network and physicochemical properties of nonsynonymous mutations in the protein-coding genes of human mitochondrial DNA

Theories on molecular evolution predict that phylogenetically recent nonsynonymous mutations should contain more non-neutral amino acid replacements than ancient mutations. We analyzed 840 complete coding-region human mitochondrial DNA (mtDNA) sequences for nonsynonymous mutations and evaluated the mutations in terms of the physicochemical properties of the amino acids involved. We identified 465 distinct missense and 6 nonsense mutations. 48% of the amino acid replacements changed polarity, 26% size, 8% charge, 32% aliphaticity, 13% aromaticity, and 44% hydropathy. The reduced-median networks of the amino acid changes revealed relatively few differences between the major continent-specific haplogroups, but a high variation and highly starlike phylogenies within the haplogroups. Some 56% of the mutations were private, and 25% were homoplasic. Nonconservative changes were more common than expected among the private mutations but less common among the homoplasic mutations. The asymptotic maximum of the number of nonsynonymous mutations in European mtDNA was estimated to be 1,081. The results suggested that amino acid replacements in the periphery of phylogenetic networks are more deleterious than those in the central parts, indicating that purifying selection prevents the fixation of some alleles.     

Genetic polymorphism of mitochondrial DNA HVS-I and HVS-Ⅱ of Chinese Tu ethnic minority group

We analyzed the two hypervariable segments HVS-Ⅰ and HVS-Ⅱ of 108 Chinese Tu ethnic minority group samples for forensic and population genetics purposes.Comparing with Anderson sequence,79 polymorphic loci in HVS-Ⅰ and 40 in HVS-Ⅱ were found in Chi-nese Tu ethnic minority group mtDNA sequences,and 90 and 64 haplotypes were then defined.Haplotype diversity and the mean pair-wise differences were 0.9903±0.0013 and 5.7785 in HVS-Ⅰ,and 0.9777±0.0013 and 3.5819 in HVS-Ⅱ,respectively.By analyzing the hypervariable domain from nucleotide 1,6180 to 1,6193 in HVS-Ⅰ,we defined some new types of sequence variations.We also compared the relationship between Tu population and other populations using mtDNA HVS-Ⅰ sequences.According to Rst genetic distances,the phylogenetic tree showed that the Tu population,the Xi'an Han population,the Chinese Korean,and the Mongol ethnic group were in a clade.This indicated a close genetic relationship between them.There were far relations between the Tu population and other Chinese southern Han populations,Siberian,European,African,and other foreign populations.The results suggest that Tu population has a multi-origin and has also merged with other local populations.     

African Haplogroup L mtDNA sequences show violations of clock-like evolution

A set of 96 complete mtDNA sequences that belong to the three major African haplogroups (L1, L2, and L3) was analyzed to determine if mtDNA has evolved as a molecular clock. Likelihood ratio tests (LRTs) were carried out with each of the haplogroups and with combined haplogroup sequence sets. Evolution has not been clock-like, neither for the coding region nor for the control region, in combined sets of African haplogroup L mtDNA sequences. In tests of individual haplogroups, L2 mtDNAs showed violations of a molecular clock under all conditions and in both the control and coding regions. In contrast, haplogroup L1 and L3 sequences, both for the coding and control regions, show clock-like evolution. In clock tests of individual L2 subclades, the L2a sequences showed a marked violation of clock-like evolution within the coding region. In addition, the L2a and L2c branch lengths of both the coding and control regions were shorter relative to those of the L2b and L2d sequences, a result that indicates lower levels of sequence divergence. Reduced median network analyses of the L2a sequences indicated the occurrence of marked homoplasy at multiple sites in the control region. After exclusion of the L2a and L2c sequences, African mtDNA coding region evolution has not significantly departed from a molecular clock, despite the results of neutrality tests that indicate the mitochondrial coding region has evolved under nonneutral conditions. In contrast, control region evolution is clock-like only at the haplogroup level, and it thus appears to have evolved essentially independently from the coding region. The results of the clock tests, the network analyses, and the branch length comparisons all caution against the use of simple mtDNA clocks.     

Phenotype of FAECB (Facility for Automated Experiments in Cell Biology) Chinese hamster ovary mutants with minimal UV-sensitivity

Mitochondrial DNA (mtDNA) haplogroup U, defined by the polymorphism 12308A>G, may constitute a risk factor for an occipital stroke in migraine. We therefore identified 14 patients with an occipital stroke and with 12308A>G. We determined complete mtDNA coding region sequence for the patients and for population controls by conformation sensitive gel electrophoresis (CSGE) and direct sequencing. Sequence information was used to construct a phylogenetic network of mtDNA haplogroups U and K, which was found to be composed of subclusters U2, U4, U5 and a new subcluster U7, as well as cluster K. Five patients with a migrainous stroke belonged to subcluster U5 (P=0.006; Fisher's exact test). Many unique mutations were found among the patients with an occipital stroke including two tRNA mutations that have previously been suggested to be pathogenic. Analysis of mtDNA sequences by CSGE and comparison of the sequences through phylogenetic analysis greatly enhances the identification of mtDNA clusters in population and detection of mtDNA mutations in patients.     

Maximum-likelihood estimation of site-specific mutation rates in human mitochondrial DNA from partial phylogenetic classification

The mitochondrial DNA hypervariable segment I (HVS-I) is widely used in studies of human evolutionary genetics, and therefore accurate estimates of mutation rates among nucleotide sites in this region are essential. We have developed a novel maximum-likelihood methodology for estimating site-specific mutation rates from partial phylogenetic information, such as haplogroup association. The resulting estimation problem is a generalized linear model, with a nonstandard link function. We develop inference and bias correction tools for our estimates and a hypothesis-testing approach for site independence. We demonstrate our methodology using 16,609 HVS-I samples from the Genographic Project. Our results suggest that mutation rates among nucleotide sites in HVS-I are highly variable. The 16,400–16,500 region exhibits significantly lower rates compared to other regions, suggesting potential functional constraints. Several loci identified in the literature as possible termination-associated sequences (TAS) do not yield statistically slower rates than the rest of HVS-I, casting doubt on their functional importance. Our tests do not reject the null hypothesis of independent mutation rates among nucleotide sites, supporting the use of site-independence assumption for analyzing HVS-I. Potential extensions of our methodology include its application to estimation of mutation rates in other genetic regions, like Y chromosome short tandem repeats.     

Do Haplogroups H and U Act to Increase the Penetrance of Alzheimer’s Disease?

1. Alzheimer’s disease (AD) is the most common form of dementia in the elderly in which interplay between genes and the environment is supposed to be involved. Mitochondrial DNA (mtDNA) has the only noncoding regions at the displacement loop (D-loop) region that contains two hypervariable segments (HVS-I and HVS-II) with high polymorphism. mtDNA has already been fully sequenced and many subsequent publications have shown polymorphic sites, haplogroups, and haplotypes. Haplogroups could have important implications to understand the association between mutability of the mitochondrial genome and the disease. 2. To assess the relationship between mtDNA haplogroup and AD, we sequenced the mtDNA HVS-I in 30 AD patients and 100 control subjects. We could find that haplogroups H and U are significantly more abundant in AD patients (P = 0.016 for haplogroup H and P = 0.0003 for haplogroup U), Thus, these two haplogroups might act synergistically to increase the penetrance of AD disease.