Functional Differences between Mitochondrial Haplogroup T and Haplogroup H in HEK293 Cybrid Cells

Background

Epidemiological case-control studies have revealed associations between mitochondrial haplogroups and the onset and/or progression of various multifactorial diseases. For instance, mitochondrial haplogroup T was previously shown to be associated with vascular diseases, including coronary artery disease and diabetic retinopathy. In contrast, haplogroup H, the most frequent haplogroup in Europe, is often found to be more prevalent in healthy control subjects than in patient study groups. However, justifications for the assumption that haplogroups are functionally distinct are rare. Therefore, we attempted to compare differences in mitochondrial function between haplogroup H and T cybrids.

Methodology/Principal Findings

Mitochondrial haplogroup H and T cybrids were generated by fusion of HEK293 cells devoid of mitochondrial DNA with isolated thrombocytes of individuals with the respective haplogroups. These cybrid cells were analyzed for oxidative phosphorylation (OXPHOS) enzyme activities, mitochondrial DNA (mtDNA) copy number, growth rate and susceptibility to reactive oxygen species (ROS). We observed that haplogroup T cybrids have higher survival rate when challenged with hydrogen peroxide, indicating a higher capability to cope with oxidative stress.

Conclusions/Significance

The results of this study show that functional differences exist between HEK293 cybrid cells which differ in mitochondrial genomic background.     

Mitochondrial DNA associations with East Asian metabolic syndrome

Mitochondrial dysfunction has repeatedly been reported associated with type 2 diabetes mellitus (T2DM) and metabolic syndrome (MS), as have mitochondrial DNA (mtDNA) tRNA and duplication mutations and mtDNA haplogroup lineages. We identified 19 Taiwanese T2DM and MS pedigrees from Taiwan, with putative matrilineal transmission, one of which harbored the pathogenic mtDNA tRNA^Leu(UUR) nucleotide (nt) 3243A>G mutation on the N9a3 haplogroup background. We then recruited three independent Taiwanese cohorts, two from Taipei (N?=?498, mean age 52 and N?=?1002, mean age 44) and one from a non-urban environment (N?=?501, mean age 57). All three cohorts were assessed for an array of metabolic parameters, their mtDNA haplogroups determined, and the haplogroups correlated with T2DM/MS phenotypes. Logistic regression analysis revealed that mtDNA haplogroups D5, F4, and N9a conferred T2DM protection, while haplogroups F4 and N9a were risk factors for hypertension (HTN), and F4 was a risk factor for obesity (OB). Additionally, the 5263C>T (ND2 A165V) variant commonly associated with F4 was associated with hypertension (HTN). Cybrids were prepared with macro-haplogroup N (defined by variants m.ND3 10398A (114T) and m.ATP6 8701A (59T)) haplogroups B4 and F1 mtDNAs and from macro-haplogroup M (variants m.ND3 10398G (114A) and m.ATP6 8701G (59A)) haplogroup M9 mtDNAs. Additionally, haplogroup B4 and F1 cybrids were prepared with and without the mtDNA variant in ND1 3394T>C (Y30H) reported to be associated with T2DM. Assay of mitochondria complex I in these cybrids revealed that macro-haplogroup N cybrids had lower activity than M cybrids, that haplogroup F cybrids had lower activity than B4 cybrids, and that the ND1 3394T>C (Y30H) variant reduced complex I on both the B4 and F1 background but with very different cumulative effects. These data support the hypothesis that functional mtDNA variants may contribute to the risk of developing T2DM and MS.     

Unexpected patterns of mitochondrial DNA variation among Native Americans from the southeastern United States

Mitochondrial DNA (mtDNA) haplogroups were determined by restriction fragment length polymorphism-typing for 66 individuals from four southeastern North American populations, and the HVS I portion of the mtDNA control region was sequenced in 48 of these individuals. Although populations from the same geographic region usually exhibit similar haplogroup frequency distributions (Lorenz and Smith [1996] Am. J. Phys. Anthropol. 101:307-323; Malhi et al. [2001] Hum. Biol. 73:17-55), those from the Southeast instead exhibit haplogroup frequency distributions that differ significantly from one another. Such divergent haplogroup frequency distributions are unexpected for the Muskogean-speaking southeastern populations, which share many sociocultural traits, speak closely related languages, and have experienced extensive admixture both with each other and with other eastern North American populations. Independent origins, genetic isolation from other Native American populations due to matrilocality, differential admixture, or a genetic bottleneck could be responsible for this heterogeneous distribution of haplogroup frequencies. Within a given haplogroup, however, the HVS I sequences from the four Muskogean-speaking populations appear relatively similar to one another, providing evidence for close relationships among them and for reduced diversity within haplogroups in the Southeast. Given additional archaeological, linguistic, and ethnographic evidence, these results suggest that a genetic bottleneck associated with the historical population decline is the most plausible explanation for such patterns of mtDNA variation.     

Mitochondrial Haplogroups Modify the Risk of Developing Hypertrophic Cardiomyopathy in a Danish Population

Hypertrophic cardiomyopathy (HCM) is a genetic disorder caused by mutations in genes coding for proteins involved in sarcomere function. The disease is associated with mitochondrial dysfunction. Evolutionarily developed variation in mitochondrial DNA (mtDNA), defining mtDNA haplogroups and haplogroup clusters, is associated with functional differences in mitochondrial function and susceptibility to various diseases, including ischemic cardiomyopathy. We hypothesized that mtDNA haplogroups, in particular H, J and K, might modify disease susceptibility to HCM. Mitochondrial DNA, isolated from blood, was sequenced and haplogroups identified in 91 probands with HCM. The association with HCM was ascertained using two Danish control populations. Haplogroup H was more prevalent in HCM patients, 60% versus 46% (p = 0.006) and 41% (p = 0.003), in the two control populations. Haplogroup J was less prevalent, 3% vs. 12.4% (p = 0.017) and 9.1%, (p = 0.06). Likewise, the UK haplogroup cluster was less prevalent in HCM, 11% vs. 22.1% (p = 0.02) and 22.8% (p = 0.04). These results indicate that haplogroup H constitutes a susceptibility factor and that haplogroup J and haplogroup cluster UK are protective factors in the development of HCM. Thus, constitutive differences in mitochondrial function may influence the occurrence and clinical presentation of HCM. This could explain some of the phenotypic variability in HCM. The fact that haplogroup H and J are also modifying factors in ischemic cardiomyopathy suggests that mtDNA haplotypes may be of significance in determining whether a physiological hypertrophy develops into myopathy. mtDNA haplotypes may have the potential of becoming significant biomarkers in cardiomyopathy.     

Decreased Reactive Oxygen Species Production in Cells with Mitochondrial Haplogroups Associated with Longevity

Mitochondrial DNA (mtDNA) is highly polymorphic, and its variations in humans may contribute to individual differences in function. Zhang and colleagues found a strikingly higher frequency of a C150T transition in the D-loop of mtDNA from centenarians and twins of an Italian population, and also demonstrated that this base substitution causes a remodeling of the mtDNA 151 replication origin in human leukocytes and fibroblasts [1]. The C150T transition is a polymorphism associated with several haplogroups. To determine whether haplogroups that carry the C150T transition display any phenotype that may be advantageous for longevity, we analyzed cybrids carrying or not the C150T transition. These cybrids were obtained by fusing cytoplasts derived from human fibroblasts with human mtDNA-less cells (ρ⁰ cells). We chose for cybrid construction and analysis haplogroup-matched pairs of fibroblast strains containing or not the C150T transition. In particular, we used, as one pair of mtDNA donors, a fibroblast strain of the U3a haplogroup, carrying the C150T transition and a strain of the U-K2 haplogroup, without the C150T transition, and as another pair, fibroblasts of the J2b haplogroup, carrying the C150T transition and of the J1c haplogroup, without the C150T transition. We have found no association of respiratory capacity, mtDNA level, mitochondrial gene expression level, or growth rate with the presence of the C150T transition. However, we have found that the cybrids with haplogroups that include the C150T transition have in common a lower reactive oxygen species (ROS) production rate than the haplogroup-matched cybrids without that transition. Thus, the lower ROS production rate may be a factor in the increased longevity associated with the U and the J2 haplogroups. Of further interest, we found that cybrids with the U3a haplogroup exhibited a higher respiration rate than the other cybrids examined.     

Mitochondrial DNA inherited variants are associated with successful aging and longevity in humans.

Mitochondrial DNA (mtDNA) is characterized by high variability, maternal inheritance, and absence of recombination. Studies of human populations have revealed ancestral associated polymorphisms whose combination defines groups of mtDNA types (haplogroups) that are currently used to reconstruct human evolution lineages. We used such inherited mtDNA markers to compare mtDNA population pools between a sample of individuals selected for successful aging and longevity (212 subjects older than 100 years and in good clinical condition) and a sample of 275 younger individuals (median age 38 years) carefully matched as to sex and geographic origin (northern and southern Italy). All nine haplogroups that are typical of Europeans were found in both samples, but male centenarians emerged in northern Italy as a particular sample: 1) mtDNA haplogroup frequency distribution was different between centenarians and younger individuals (P=0.017 by permutation tests); and 2) the frequency of the J haplogroup was notably higher in centenarians than in younger individuals (P=0.0052 by Fisher exact test). Since haplogroups are defined on the basis of inherited variants, these data show that mtDNA inherited variability could play a role in successful aging and longevity.     

Experimental assessment of bioenergetic differences caused by the common European mitochondrial DNA haplogroups H and T

Studies of both survival after sepsis and sperm motility in human populations have shown significant associations with common European mitochondrial DNA haplogroups, and have led to proposals that mitochondria bearing haplogroup H have different bioenergetic capacities than those bearing haplogroup T. However, the validity of such associations assumes that there are no non-random influences of nuclear genes or other factors. Here, we removed the effect of any differences in nuclear genes by constructing transmitochondrial cybrids harbouring mitochondria with either haplogroup H or haplogroup T in cultured A549 human lung carcinoma cells with identical nuclear backgrounds. We compared the bioenergetic capacities and coupling efficiencies of mitochondria isolated from these cells, and of mitochondria retained within the cells, as a critical experimental test of the hypothesis that these haplogroups affect mitochondrial bioenergetics. We found that there were no functionally-important bioenergetic differences between mitochondria bearing these haplogroups, using either isolated mitochondria or mitochondria within cells.     

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.     

Oxidative phosphorylation differences between mitochondrial DNA haplogroups modify the risk of Leber's hereditary optic neuropathy

Leber's hereditary optic neuropathy is a maternally inherited optic atrophy caused by mitochondrial DNA point mutations. Previous epidemiological studies have shown that individuals from mitochondrial genetic backgrounds (haplogroups) J/Uk and H have a higher and a lower risk, respectively, of suffering this disorder. To analyze the bases of these associations at cellular and molecular levels, functional studies with cybrids provide high quality evidence. Cybrids from haplogroup J contain less mitochondrial deoxyribonucleic acid (DNA) and ribonucleic acid (RNA) and synthesize a smaller amount of mitochondrial DNA-encoded polypeptides than those from haplogroup H. Haplogroup J cybrids also display lower oxygen consumption, mitochondrial inner membrane potential and total adenosine-5'-triphosphate (ATP) levels. Moreover, mitochondrial DNA levels correlate with many parameters of the oxidative phosphorylation system. These results suggest that the mitochondrial DNA amount determines oxidative phosphorylation capacity and, along with other recently published observations, support the possibility that mitochondrial DNA levels may be responsible for the bias of the disorder toward males, for the incomplete penetrance of mutations causing Leber's hereditary optic neuropathy and for the association of the disease with particular mitochondrial DNA haplogroups.     

mtDNA evidence: genetic background associated with related populations at high risk for esophageal cancer between Chaoshan and Taihang Mountain areas in China

There are three major geographic regions in China known for their high incidences of esophageal cancer (EC): the Taihang Mountain range of north-central China, the Minnan area of Fujian province, and the Chaoshan plain of Guangdong province. Historically, waves of great population migrations from north-central China through coastal Fujian to the Chaoshan plain were recorded. To study the genetic relationship among the related EC high-risk populations, we analyzed mitochondrial DNA (mtDNA) haplogroups based on 30 EC patients from Chaoshan and used control samples from the high-risk populations, including 48, 73, and 89 subjects from the Taihang, Fujian, and Chaoshan areas, respectively. The principal component of all haplogroups, correlation analysis of haplogroup frequency distributions between populations, and haplogroup D network analysis showed that compared with other Chinese populations, populations in the three studied areas are genetically related. The highest haplogroup frequency shared by all studied populations was haplogroup D, with much higher frequency in the Chaoshan area EC patients. The majority of haplogroup D individuals among the Chaoshan area EC patients belonged to subhaplogroups D4a and D5a, with the total frequency of these two haplogroups significantly higher than that in the high-risk population in the same area (chi(2)=9.017, p<0.01). In conclusion, EC high-risk populations in these three areas share a similar matrilineal genetic background, and D4a and D5a might be candidate genetic markers for screening populations susceptible to EC in the Chaoshan area. Ours is the first report to show the association between mtDNA haplogroups (D4a and D5a) and esophageal cancer.     

Major mitochondrial DNA haplotype heterogeneity in highland and lowland Amerindian populations from Bolivia

This study provides the frequencies of four mitochondrial DNA (mtDNA) haplogroups of 233 native South Amerindians in eight populations living in the Beni Department of Bolivia, including six populations not previously studied. Linguistically, these populations belong to the three principal South Amerindian language stocks, Andean, Equatorial-Tucanoan, and Ge-Pano-Carib. Frequency analyses under geographic, historic, linguistic, and genetic configurations using the theta statistic of Weir (Weir 1990) and analysis of molecular variance (AMOVA) show similar results. Results are also similar when phenetic cluster is used. Aymara belongs almost exclusively to haplogroup B, Quechua- and Moseten-speaking tribes belong to haplogroups A and B, but the first tribe presents high frequencies of haplogroup B. Yuracare, Trinitario, and Ignaciano exhibit high frequencies of A, B, and C haplogroups, and the Movima present a large proportion of haplogroup C. There is some correspondence between mtDNA haplogroup frequencies and language affiliation and historical connections, but less so with geographic aspects. The present study provides a context for understanding the relationship between different Amerindian populations living in a multiethnic area of Bolivia.     

Analysis of mtDNA variation in African populations reveals the most ancient of all human continent-specific haplogroups.

mtDNA sequence variation was examined in 140 Africans, including Pygmies from Zaire and Central African Republic (C.A.R.) and Mandenkalu, Wolof, and Pular from Senegal. More than 76% of the African mtDNAs (100% of the Pygmies and 67.3% of the Senegalese) formed one major mtDNA cluster (haplogroup L) defined by an African-specific HpaI site gain at nucleotide pair (np) 3592. Additional mutations subdivided haplogroup L into two subhaplogroups, each encompassing both Pygmy and Senegalese mtDNAs. A novel 12-bp homoplasmic insertion in the intergenic region between tRNA(Tyr) and cytochrome oxidase I (COI) genes was also observed in 17.6% of the Pygmies from C.A.R. This insertion is one of the largest observed in human mtDNAs. Another 25% of the Pygmy mtDNAs harbored a 9-bp deletion between the cytochrome oxidase II (COII) and tRNA(Lys) genes, a length polymorphism previously reported in non-African populations. In addition to haplogroup L, other haplogroups were observed in the Senegalese. These haplogroups were more similar to those observed in Europeans and Asians than to haplogroup L mtDNAs, suggesting that the African mtDNAs without the HpaI np 3592 site could be the ancestral types from which European and Asian mtDNAs were derived. Comparison of the intrapopulation sequence divergence in African and non-African populations confirms that African populations exhibit the largest extent of mtDNA variation, a result that further supports the hypothesis that Africans represent the most ancient human group and that all modern humans have a common and recent African origin. The age of the total African variation was estimated to be 101,000-133,000 years before present (YBP), while the age of haplogroup L was estimated at 98,000-130,000 YBP. These values substantially exceed the ages of all Asian- and European-specific mtDNA haplogroups.     

Mitochondrial DNA variation in human metabolic rate and energy expenditure

The role of climate in driving selection of mtDNA as Homo sapiens migrated out of Africa into Eurasia remains controversial. We evaluated the role of mtDNA variation in resting metabolic rate (RMR) and total energy expenditure (TEE) among 294 older, community-dwelling African and European American adults from the Health, Aging and Body Composition Study. Common African haplogroups L0, L2 and L3 had significantly lower RMRs than European haplogroups H, JT and UK with haplogroup L1 RMR being intermediate to these groups. This study links mitochondrial haplogroups with ancestry-associated differences in metabolic rate and energy expenditure.     

The Background of Mitochondrial DNA Haplogroup J Increases the Sensitivity of Leber's Hereditary Optic Neuropathy Cells to 2,5-Hexanedione Toxicity

Leber''s hereditary optic neuropathy (LHON) is a maternally inherited blinding disease due to mitochondrial DNA (mtDNA) point mutations in complex I subunit genes, whose incomplete penetrance has been attributed to both genetic and environmental factors. Indeed, the mtDNA background defined as haplogroup J is known to increase the penetrance of the 11778/ND4 and 14484/ND6 mutations. Recently it was also documented that the professional exposure to n-hexane might act as an exogenous trigger for LHON. Therefore, we here investigate the effect of the n-hexane neurotoxic metabolite 2,5-hexanedione (2,5-HD) on cell viability and mitochondrial function of different cell models (cybrids and fibroblasts) carrying the LHON mutations on different mtDNA haplogroups. The viability of control and LHON cybrids and fibroblasts, whose mtDNAs were completely sequenced, was assessed using the MTT assay. Mitochondrial ATP synthesis rate driven by complex I substrates was determined with the luciferine/luciferase method. Incubation with 2,5-HD caused the maximal loss of viability in control and LHON cells. The toxic effect of this compound was similar in control cells irrespective of the mtDNA background. On the contrary, sensitivity to 2,5-HD induced cell death was greatly increased in LHON cells carrying the 11778/ND4 or the 14484/ND6 mutation on haplogroup J, whereas the 11778/ND4 mutation in association with haplogroups U and H significantly improved cell survival. The 11778/ND4 mutation on haplogroup U was also more resistant to inhibition of complex I dependent ATP synthesis by 2,5-HD. In conclusion, this study shows that mtDNA haplogroups modulate the response of LHON cells to 2,5-HD. In particular, haplogroup J makes cells more sensitive to its toxic effect. This is the first evidence that an mtDNA background plays a role by interacting with an environmental factor and that 2,5-HD may be a risk element for visual loss in LHON. This proof of principle has broad implications for other neurodegenerative disorders such as Parkinson''s disease.     

Mitochondrial DNA (mtDNA) haplogroups in 1526 unrelated individuals from 11 Departments of Colombia

The frequencies of four mitochondrial Native American DNA haplogroups were determined in 1526 unrelated individuals from 11 Departments of Colombia and compared to the frequencies previously obtained for Amerindian and Afro-Colombian populations. Amerindian mtDNA haplogroups ranged from 74% to 97%. The lowest frequencies were found in Departments on the Caribbean coast and in the Pacific region, where the frequency of Afro-Colombians is higher, while the highest mtDNA Amerindian haplogroup frequencies were found in Departments that historically have a strong Amerindian heritage. Interestingly, all four mtDNA haplogroups were found in all Departments, in contrast to the complete absence of haplogroup D and high frequencies of haplogroup A in Amerindian populations in the Caribbean region of Colombia. Our results indicate that all four Native American mtDNA haplogroups were widely distributed in Colombia at the time of the Spanish conquest.     

Respiratory function in cybrid cell lines carrying European mtDNA haplogroups: implications for Leber's hereditary optic neuropathy

The possibility that some combinations of mtDNA polymorphisms, previously associated with Leber's hereditary optic neuropathy (LHON), may affect mitochondrial respiratory function was tested in osteosarcoma-derived transmitochondrial cytoplasmic hybrids (cybrids). In this cellular system, in the presence of the same nuclear background, different exogenous mtDNAs are used to repopulate a parental cell line previously devoid of its original mtDNA. No detectable differences in multiple parameters exploring respiratory function were observed when mtDNAs belonging to European haplogroups X, H, T and J were used. Different possible explanations for the previously established association between haplogroup J and LHON 11778/ND4 and 14484/ND6 pathogenic mutations are discussed, including the unconventional proposal that mtDNA haplogroup J may exert a protective rather than detrimental effect.     

Uniparental genetic markers in South Amerindians

A comprehensive review of uniparental systems in South Amerindians was undertaken. Variability in the Y-chromosome haplogroups were assessed in 68 populations and 1,814 individuals whereas that of Y-STR markers was assessed in 29 populations and 590 subjects. Variability in the mitochondrial DNA (mtDNA) haplogroup was examined in 108 populations and 6,697 persons, and sequencing studies used either the complete mtDNA genome or the highly variable segments 1 and 2. The diversity of the markers made it difficult to establish a general picture of Y-chromosome variability in the populations studied. However, haplogroup Q1a3a* was almost always the most prevalent whereas Q1a3* occurred equally in all regions, which suggested its prevalence among the early colonizers. The STR allele frequencies were used to derive a possible ancient Native American Q-clade chromosome haplotype and five of six STR loci showed significant geographic variation. Geographic and linguistic factors moderately influenced the mtDNA distributions (6% and 7%, respectively) and mtDNA haplogroups A and D correlated positively and negatively, respectively, with latitude. The data analyzed here provide rich material for understanding the biological history of South Amerindians and can serve as a basis for comparative studies involving other types of data, such as cultural data.     

Asian affinities and continental radiation of the four founding Native American mtDNAs.

The mtDNA variation of 321 individuals from 17 Native American populations was examined by high-resolution restriction endonuclease analysis. All mtDNAs were amplified from a variety of sources by using PCR. The mtDNA of a subset of 38 of these individuals was also analyzed by D-loop sequencing. The resulting data were combined with previous mtDNA data from five other Native American tribes, as well as with data from a variety of Asian populations, and were used to deduce the phylogenetic relationships between mtDNAs and to estimate sequence divergences. This analysis revealed the presence of four haplotype groups (haplogroups A, B, C, and D) in the Amerind, but only one haplogroup (A) in the Na-Dene, and confirmed the independent origins of the Amerinds and the Na-Dene. Further, each haplogroup appeared to have been founded by a single mtDNA haplotype, a result which is consistent with a hypothesized founder effect. Most of the variation within haplogroups was tribal specific, that is, it occurred as tribal private polymorphisms. These observations suggest that the process of tribalization began early in the history of the Amerinds, with relatively little intertribal genetic exchange occurring subsequently. The sequencing of 341 nucleotides in the mtDNA D-loop revealed that the D-loop sequence variation correlated strongly with the four haplogroups defined by restriction analysis, and it indicated that the D-loop variation, like the haplotype variation, arose predominantly after the migration of the ancestral Amerinds across the Bering land bridge.     

Microevolution in prehistoric Andean populations: chronologic mtDNA variation in the desert valleys of northern Chile

Archeological evidence suggests that the iconographic and technological developments that took place in the highlands around Lake Titicaca in the Central Andean region had an influence on the cultural elaborations of the human groups in the valleys and the Pacific coast of northern Chile. In a previous communication, we were able to show, by means of a distance analysis, that a craniofacial differentiation accompanied the process of cultural evolution in the valleys (Rothhammer and Santoro [2001] Lat. Am. Antiq. 12:59-66). Recently, numerous South Amerindian mtDNA studies were published, and more accurate molecular techniques to study ancient mtDNA are available. In view of these recent developments, we decided 1) to study chronological changes of ancient mtDNA haplogroup frequencies in the nearby Lluta, Azapa, and Camarones Valleys, 2) to identify microevolutionary forces responsible for such changes, and 3) to compare ancient mtDNA haplogroup frequencies with previous data in order to validate craniometrical results and to reconstruct the biological history of the prehistoric valley groups in the context of their interaction with culturally more developed highland populations. From a total of 97 samples from 83 individuals, 68 samples (61 individuals) yielded amplifications for the fragments that harbor classical mtDNA markers. The haplogroup distribution among the total sample was as follows: 26.2%, haplogroup A; 34.4%, haplogroup B; 14.8%, haplogroup C; 3.3%, haplogroup D; and 21.3%, other haplogroups. Haplogroup B tended to increase, and haplogroup A to decrease during a 3,900-year time interval. The sequence data are congruent with the haplogroup analysis. In fact, the sequencing of hypervariable region I of 30 prehistoric individuals revealed 43 polymorphic sites. Sequence alignment and subsequent phylogenetic tree construction showed two major clusters associated with the most common restriction haplogroups. Individuals belonging to haplogroups C and D tended to cluster together with nonclassical lineages.