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.     

African-American mitochondrial DNAs often match mtDNAs found in multiple African ethnic groups

Background  

Mitochondrial DNA (mtDNA) haplotypes have become popular tools for tracing maternal ancestry, and several companies offer this service to the general public. Numerous studies have demonstrated that human mtDNA haplotypes can be used with confidence to identify the continent where the haplotype originated. Ideally, mtDNA haplotypes could also be used to identify a particular country or ethnic group from which the maternal ancestor emanated. However, the geographic distribution of mtDNA haplotypes is greatly influenced by the movement of both individuals and population groups. Consequently, common mtDNA haplotypes are shared among multiple ethnic groups. We have studied the distribution of mtDNA haplotypes among West African ethnic groups to determine how often mtDNA haplotypes can be used to reconnect Americans of African descent to a country or ethnic group of a maternal African ancestor. The nucleotide sequence of the mtDNA hypervariable segment I (HVS-I) usually provides sufficient information to assign a particular mtDNA to the proper haplogroup, and it contains most of the variation that is available to distinguish a particular mtDNA haplotype from closely related haplotypes. In this study, samples of general African-American and specific Gullah/Geechee HVS-I haplotypes were compared with two databases of HVS-I haplotypes from sub-Saharan Africa, and the incidence of perfect matches recorded for each sample.     

Mitochondrial DNA diversity in Southeast Asian populations

In a previous study of Southeast Asian genetic variation, we characterized mitochondrial DNAs (mtDNAs) from six populations through high-resolution restriction fragment length polymorphism (RFLP) analysis. Our analysis revealed that these Southeast Asian populations were genetically similar to each other, suggesting they had a common origin. However, other patterns of population associations also emerged. Haplotypes from a major founding haplogroup in Papua New Guinea were present in Malaysia; the Vietnamese and Malaysian aborigines (Orang Asli) had high frequencies of haplogroup F, which was also seen in most other Southeast Asian populations; and haplogroup B, defined by the Region V 9-base-pair deletion, was present throughout the region. In addition, the Malaysian and Sabah (Borneo) aborigine populations exhibited a number of unique mtDNA clusters that were not observed in other populations. Unfortunately, it has been difficult to compare these patterns of genetic diversity with those shown in subsequent studies of mtDNA variation in Southeast Asian populations because the latter have typically sequenced the first hypervariable segment (HVS-I) of the control region (CR) sequencing rather than used RFLP haplotyping to characterize the mtDNAs present in them. For this reason, we sequenced the HVS-I of Southeast Asian mtDNAs that had previously been subjected to RFLP analysis, and compared the resulting data with published information from other Southeast Asian and Oceanic groups. Our findings reveal broad patterns of mtDNA haplogroup distribution in Southeast Asia that may reflect different population expansion events in this region over the past 50,000-5,000 years.     

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.     

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.     

MtDNA variation in the Altai-Kizhi population of southern Siberia: a synthesis of genetic variation

The native peoples of Gorno Altai in southern Siberia represent a genetically diverse population and have been of great interest to anthropological genetics. In particular, the southern Altaian population is argued to be the best candidate for the New World ancestral population. In this study we sampled Altai-Kizhi from the southern Altaian village of Mendur-Sokkon, analyzed mtDNA RFLP markers and HVS-I sequences, and compared the results to other published mtDNA data from Derenko et al. (2003) and Shields et al. (1993) encompassing the same region. Because each independent study uses different sampling techniques in characterizing gene pools, in this paper we explore the accuracy and reliability of evolutionary studies on human populations. All the major Native American haplogroups (A, B, C, and D) were identified in the Mendur-Sokkon sample, including a single individual belonging to haplogroup X. The most common mtDNA lineages are C (35.7%) and D (13.3%), which is consistent with the haplogroup profiles of neighboring Siberian groups. The Mendur-Sokkon sample exhibits depressed HVS-I diversity values and neutrality test scores, which starkly differs from the Derenko et al. (2003) data set and more closely resembles the results for neighboring south Siberian groups. Furthermore, the multidimensional scaling plot of DA genetic distances does not cluster the Altai samples, showing different genetic affinities with various Asian groups. The findings underscore the importance of sampling strategy in the reconstruction of evolutionary history at the population level.     

Disuniting uniformity: a pied cladistic canvas of mtDNA haplogroup H in Eurasia

It has been often stated that the overall pattern of human maternal lineages in Europe is largely uniform. Yet this uniformity may also result from an insufficient depth and width of the phylogenetic analysis, in particular of the predominant western Eurasian haplogroup (Hg) H that comprises nearly a half of the European mitochondrial DNA (mtDNA) pool. Making use of the coding sequence information from 267 mtDNA Hg H sequences, we have analyzed 830 mtDNA genomes, from 11 European, Near and Middle Eastern, Central Asian, and Altaian populations. In addition to the seven previously specified subhaplogroups, we define fifteen novel subclades of Hg H present in the extant human populations of western Eurasia. The refinement of the phylogenetic resolution has allowed us to resolve a large number of homoplasies in phylogenetic trees of Hg H based on the first hypervariable segment (HVS-I) of mtDNA. As many as 50 out of 125 polymorphic positions in HVS-I were found to be mutated in more than one subcluster of Hg H. The phylogeographic analysis revealed that sub-Hgs H1*, H1b, H1f, H2a, H3, H6a, H6b, and H8 demonstrate distinct phylogeographic patterns. The monophyletic subhaplogroups of Hg H provide means for further progress in the understanding of the (pre)historic movements of women in Eurasia and for the understanding of the present-day genetic diversity of western Eurasians in general.     

Evidence for mtDNA admixture between the Finns and the Saami

OBJECTIVES: The Finns, and to a more extreme extent the Saami, are genetic outliers in Europe. Despite the close geographical contact between these populations, no major contribution of Saami mtDNA haplotypes to the Finnish population has been detected. METHODS: To examine the extent of maternal gene flow from the Saami into Finnish populations, we determined the mtDNA variation in 403 persons living in four provinces in central and northern Finland. For all of these samples, we assessed the frequencies of mtDNA haplogroups and examined sequence variation in the hypervariable segment I (HVS-I). The resulting data were compared with published information for Saami populations. RESULTS: The frequencies of the mtDNA haplogroups differed between the populations of the four provinces, suggesting a distinction between northern and central Finland. Analysis of molecular variance suggested that the Saami deviated less from the population of northern Finland than from that of central Finland. Five HVS-I haplotypes, including that harboring the Saami motif and the Asian-specific haplogroup Z, were shared between the Finns and the Saami and allowed comparisons between the populations. Their frequency was highest in the Saami and decreased towards central Finland. CONCLUSIONS: The high frequency of certain mtDNA haplotypes considered to be Saami specific in the Finnish population suggests a genetic admixture, which appears to be more pronounced in northern Finland. Furthermore, the presence of haplogroup Z in the Finns and the Saami indicates that traces of Asian mtDNA genotypes have survived in the contemporary populations.     

The Peopling of Korea Revealed by Analyses of Mitochondrial DNA and Y-Chromosomal Markers

Background

The Koreans are generally considered a northeast Asian group because of their geographical location. However, recent findings from Y chromosome studies showed that the Korean population contains lineages from both southern and northern parts of East Asia. To understand the genetic history and relationships of Korea more fully, additional data and analyses are necessary.

Methodology and Results

We analyzed mitochondrial DNA (mtDNA) sequence variation in the hypervariable segments I and II (HVS-I and HVS-II) and haplogroup-specific mutations in coding regions in 445 individuals from seven east Asian populations (Korean, Korean-Chinese, Mongolian, Manchurian, Han (Beijing), Vietnamese and Thais). In addition, published mtDNA haplogroup data (N = 3307), mtDNA HVS-I sequences (N = 2313), Y chromosome haplogroup data (N = 1697) and Y chromosome STR data (N = 2713) were analyzed to elucidate the genetic structure of East Asian populations. All the mtDNA profiles studied here were classified into subsets of haplogroups common in East Asia, with just two exceptions. In general, the Korean mtDNA profiles revealed similarities to other northeastern Asian populations through analysis of individual haplogroup distributions, genetic distances between populations or an analysis of molecular variance, although a minor southern contribution was also suggested. Reanalysis of Y-chromosomal data confirmed both the overall similarity to other northeastern populations, and also a larger paternal contribution from southeastern populations.

Conclusion

The present work provides evidence that peopling of Korea can be seen as a complex process, interpreted as an early northern Asian settlement with at least one subsequent male-biased southern-to-northern migration, possibly associated with the spread of rice agriculture.     

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.     

Characterization of mitochondrial haplogroups in a large population-based sample from the United States

Mitochondrial DNA (mtDNA) haplogroups are valuable for investigations in forensic science, molecular anthropology, and human genetics. In this study, we developed a custom panel of 61 mtDNA markers for high-throughput classification of European, African, and Native American/Asian mitochondrial haplogroup lineages. Using these mtDNA markers, we constructed a mitochondrial haplogroup classification tree and classified 18,832 participants from the National Health and Nutrition Examination Surveys (NHANES). To our knowledge, this is the largest study to date characterizing mitochondrial haplogroups in a population-based sample from the United States, and the first study characterizing mitochondrial haplogroup distributions in self-identified Mexican Americans separately from Hispanic Americans of other descent. We observed clear differences in the distribution of maternal genetic ancestry consistent with proposed admixture models for these subpopulations, underscoring the genetic heterogeneity of the United States Hispanic population. The mitochondrial haplogroup distributions in the other self-identified racial/ethnic groups within NHANES were largely comparable to previous studies. Mitochondrial haplogroup classification was highly concordant with self-identified race/ethnicity (SIRE) in non-Hispanic whites (94.8 %), but was considerably lower in admixed populations including non-Hispanic blacks (88.3 %), Mexican Americans (81.8 %), and other Hispanics (61.6 %), suggesting SIRE does not accurately reflect maternal genetic ancestry, particularly in populations with greater proportions of admixture. Thus, it is important to consider inconsistencies between SIRE and genetic ancestry when performing genetic association studies. The mitochondrial haplogroup data that we have generated, coupled with the epidemiologic variables in NHANES, is a valuable resource for future studies investigating the contribution of mtDNA variation to human health and disease.     

Estimates of Continental Ancestry Vary Widely among Individuals with the Same mtDNA Haplogroup

The association between a geographical region and an mtDNA haplogroup(s) has provided the basis for using mtDNA haplogroups to infer an individual’s place of origin and genetic ancestry. Although it is well known that ancestry inferences using mtDNA haplogroups and those using genome-wide markers are frequently discrepant, little empirical information exists on the magnitude and scope of such discrepancies between multiple mtDNA haplogroups and worldwide populations. We compared genetic-ancestry inferences made by mtDNA-haplogroup membership to those made by autosomal SNPs in ∼940 samples of the Human Genome Diversity Panel and recently admixed populations from the 1000 Genomes Project. Continental-ancestry proportions often varied widely among individuals sharing the same mtDNA haplogroup. For only half of mtDNA haplogroups did the highest average continental-ancestry proportion match the highest continental-ancestry proportion of a majority of individuals with that haplogroup. Prediction of an individual’s mtDNA haplogroup from his or her continental-ancestry proportions was often incorrect. Collectively, these results indicate that for most individuals in the worldwide populations sampled, mtDNA-haplogroup membership provides limited information about either continental ancestry or continental region of origin.     

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.     

Niger-Congo speaking populations and the formation of the Brazilian gene pool: mtDNA and Y-chromosome data

We analyzed sequence variation in the mitochondrial DNA (mtDNA) hypervariable segment I (HVS-I) from 201 Black individuals from two Brazilian cities (Rio de Janeiro and Porto Alegre), and compared these data with published information from 21 African populations. A subset of 187 males of the sample was also characterized for 30 Y-chromosome biallelic polymorphisms, and the data were compared with those from 48 African populations. The mtDNA data indicated that respectively 69% and 82% of the matrilineages found in Rio de Janeiro and Porto Alegre originated from West-Central/Southeast Africa. These estimates are in close agreement with historical records which indicated that most of the Brazilian slaves who arrived in Rio de Janeiro were from West-Central Africa. In contrast to mtDNA, Y-chromosome haplogroup analysis did not allow discrimination between places of origin in West or West-Central Africa. Thus, when comparing these two major African regions, there seems to be higher genetic structure with mtDNA than with Y-chromosome data.     

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 for European mtDNA

The sequence in the first hypervariable segment (HVS-I) of the control region has been used as a source of evolutionary information in most phylogenetic analyses of mtDNA. Population genetic inference would benefit from a better understanding of the variation in the mtDNA coding region, but, thus far, complete mtDNA sequences have been rare. We determined the nucleotide sequence in the coding region of mtDNA from 121 Finns, by conformation-sensitive gel electrophoresis and subsequent sequencing and by direct sequencing of the D loop. Furthermore, 71 sequences from our previous reports were included, so that the samples represented all the mtDNA haplogroups present in the Finnish population. We found a total of 297 variable sites in the coding region, which allowed the compilation of unambiguous phylogenetic networks. The D loop harbored 104 variable sites, and, in most cases, these could be localized within the coding-region networks, without discrepancies. Interestingly, many homoplasies were detected in the coding region. Nucleotide variation in the rRNA and tRNA genes was 6%, and that in the third nucleotide positions of structural genes amounted to 22% of that in the HVS-I. The complete networks enabled the relationships between the mtDNA haplogroups to be analyzed. Phylogenetic networks based on the entire coding-region sequence in mtDNA provide a rich source for further population genetic studies, and complete sequences make it easier to differentiate between disease-causing mutations and rare polymorphisms.     

Genetic analysis on Tuoba Xianbei remains excavated from Qilang Mountain Cemetery in Qahar Right Wing Middle Banner of Inner Mongolia

Sixteen sequences of the hypervariable segment I (HVS-I, 16039-16398) in mtDNA control region from ancient Tuoba Xianbei remains excavated from Qilang Mountain Cemetery were analyzed. In which, 13 haplotypes were found by 25 polymorphic sites. The haplotype diversity and nucleotide diversity were 0.98 and 0.0189, respectively, and the mean of nucleotide number differences was 6.25. Haplogroup analysis indicates these remains mainly belong to haplogroup C (31.25%) and D (43.75%). According to the published data were considered, we can suggest that the Tuoba Xianbei presented a close genetic affinity to Oroqen, Outer Mongolian and Ewenki populations, especially Oroqen.     

Genetic analyses of Xianbei populations about 1,500–1,800 years old

To understand the profile of genetic structure of Xianbei and trace its impacts on the formation and development of the minorities from northern China, we analyzed the sequences of the hypervariable segment I (HVS-I, 16.035–16.398) in mtDNA control region of 17 Tuoba Xianbei remains from Shangdu Dongdajing cemetery (Inner Mongolia). Its haplotype diversity and nucleotide diversity were 0.971 ± 0.032 and 0.0184 ± 0.010, respectively, and the haplogroup status presented 29.5% C, 23.5% D4, 17.6% D5, 17.6%A, 5.9% B and 5.9% G. When the data from Qilang Mountain Tuoba remains and other relevant populations were considered, we found that Dongdajing Tuoba Xianbei presented the closest genetic affinity to Qilang Mountain Tuoba Xianbei. Tuoba Xianbei and Murong Xianbei showed a significant differentiation in the maternal lineages. Tuoba Xianbei may contribute to the gene pool of some northern minorities, and it may mix with Xiongnu in northern China.     

Mitochondrial DNA haplogroup analysis reveals no association between the common genetic lineages and prostate cancer in the Korean population

Mitochondrial DNA (mtDNA) variation has recently been suggested to have an association with various cancers, including prostate cancer risk, in human populations. Since mtDNA is haploid and lacks recombination, specific mutations in the mtDNA genome associated with human diseases arise and remain in particular genetic backgrounds referred to as haplogroups. To assess the possible contribution of mtDNA haplogroup-specific mutations to the occurrence of prostate cancer, we have therefore performed a population-based study of a prostate cancer cases and corresponding controls from the Korean population. No statistically significant difference in the distribution of mtDNA haplogroup frequencies was observed between the case and control groups of Koreans. Thus, our data imply that specific mtDNA mutations/lineages did not appear to have a significant effect on a predisposition to prostate cancer in the Korean population, although larger sample sizes are necessary to validate our results.