Genome digging: insight into the mitochondrial genome of Homo

Background

A fraction of the Neanderthal mitochondrial genome sequence has a similarity with a 5,839-bp nuclear DNA sequence of mitochondrial origin (numt) on the human chromosome 1. This fact has never been interpreted. Although this phenomenon may be attributed to contamination and mosaic assembly of Neanderthal mtDNA from short sequencing reads, we explain the mysterious similarity by integration of this numt (mtAncestor-1) into the nuclear genome of the common ancestor of Neanderthals and modern humans not long before their reproductive split.

Principal Findings

Exploiting bioinformatics, we uncovered an additional numt (mtAncestor-2) with a high similarity to the Neanderthal mtDNA and indicated that both numts represent almost identical replicas of the mtDNA sequences ancestral to the mitochondrial genomes of Neanderthals and modern humans. In the proteins, encoded by mtDNA, the majority of amino acids distinguishing chimpanzees from humans and Neanderthals were acquired by the ancestral hominins. The overall rate of nonsynonymous evolution in Neanderthal mitochondrial protein-coding genes is not higher than in other lineages. The model incorporating the ancestral hominin mtDNA sequences estimates the average divergence age of the mtDNAs of Neanderthals and modern humans to be 450,000–485,000 years. The mtAncestor-1 and mtAncestor-2 sequences were incorporated into the nuclear genome approximately 620,000 years and 2,885,000 years ago, respectively.

Conclusions

This study provides the first insight into the evolution of the mitochondrial DNA in hominins ancestral to Neanderthals and humans. We hypothesize that mtAncestor-1 and mtAncestor-2 are likely to be molecular fossils of the mtDNAs of Homo heidelbergensis and a stem Homo lineage. The d_N/d_S dynamics suggests that the effective population size of extinct hominins was low. However, the hominin lineage ancestral to humans, Neanderthals and H. heidelbergensis, had a larger effective population size and possessed genetic diversity comparable with those of chimpanzee and gorilla.     

mtDNA variation in the South African Kung and Khwe-and their genetic relationships to other African populations

The mtDNA variation of 74 Khoisan-speaking individuals (Kung and Khwe) from Schmidtsdrift, in the Northern Cape Province of South Africa, was examined by high-resolution RFLP analysis and control region (CR) sequencing. The resulting data were combined with published RFLP haplotype and CR sequence data from sub-Saharan African populations and then were subjected to phylogenetic analysis to deduce the evolutionary relationships among them. More than 77% of the Kung and Khwe mtDNA samples were found to belong to the major mtDNA lineage, macrohaplogroup L* (defined by a HpaI site at nucleotide position 3592), which is prevalent in sub-Saharan African populations. Additional sets of RFLPs subdivided macrohaplogroup L* into two extended haplogroups-L1 and L2-both of which appeared in the Kung and Khwe. Besides revealing the significant substructure of macrohaplogroup L* in African populations, these data showed that the Biaka Pygmies have one of the most ancient RFLP sublineages observed in African mtDNA and, thus, that they could represent one of the oldest human populations. In addition, the Kung exhibited a set of related haplotypes that were positioned closest to the root of the human mtDNA phylogeny, suggesting that they, too, represent one of the most ancient African populations. Comparison of Kung and Khwe CR sequences with those from other African populations confirmed the genetic association of the Kung with other Khoisan-speaking peoples, whereas the Khwe were more closely linked to non-Khoisan-speaking (Bantu) populations. Finally, the overall sequence divergence of 214 African RFLP haplotypes defined in both this and an earlier study was 0.364%, giving an estimated age, for all African mtDNAs, of 125,500-165,500 years before the present, a date that is concordant with all previous estimates derived from mtDNA and other genetic data, for the time of origin of modern humans in Africa.     

Mapping of mitochondrial DNA of individual sheep and goats: Rapid evolution in the D loop region

Mitochondrial DNA (mtDNA) from sheep and goat was compared by restriction endonuclease analysis and heteroduplex mapping in the electron microscope. The fragment patterns produced by endonuclease Hae III from three individual sheep and two goat mtDNAs all differed from each other. The three sheep mtDNAs had identical Eco RI and Hind III fragments, but the two goat mtDNA patterns differed from each other as well as from sheep mtDNA. We estimate that each sheep mtDNA differs from each other by 0.5–1% of its nucleotide sequences, the two goat mtDNAs by 1–2%, and there is a 6–11% sequence difference between sheep and goat mtDNAs. We have mapped the Eco RI and Hind III sites of goat and sheep mtDNA and determined the positions of the D loop, which marks the replication origin, relative to the restriction map. The D loops are at homologous positions on the mtDNAs from both species, but the goat D loop is only 75% as long as the sheep D loop. Regions with a high degree of sequence divergence occur at both ends of the D loop. We suggest that a duplication of about 150 base pairs has occurred in the region where the sheep and goat D loops differ in length. We discuss mtDNA evolution in terms of divergence of isolated “mitochondrial DNA clones.”     

A high frequency of mtDNA polymorphisms in HeLa cell sublines

The complete mtDNA sequences from the uncloned "founder" HeLa cells and from five sublines have been determined. These sequences all carry a common "core" of 38 single basepair alterations relative to the revised Cambridge Reference Sequence (CRS). The HeLa mitochondrial genome is of African descent and it is a member of the African L3 haplogroup. The sequence of the HeLa mtDNA resolves the uncertainty surrounding the mosaic composition of the original CRS for human mtDNA. Most importantly, we detected a total of eight polymorphisms that have arisen in the mtDNA coding region of different HeLa sublines. These observations suggest that HeLa mtDNA has a high rate of sequence divergence, relative to the phylogenetically-derived divergence rate for mtDNAs in the human population, which results from a relaxation of negative selection against the fixation of deleterious mutations. Furthermore, this high frequency of polymorphisms in HeLa mtDNA may reflect a process similar to the accumulation of somatic mtDNA mutations in human cancers. Preliminary analysis of single-cell derived subclone lines revealed the occurrence of another polymorphism and provided evidence for a large number of mtDNA segregation units.     

A complete Neandertal mitochondrial genome sequence determined by high-throughput sequencing

A complete mitochondrial (mt) genome sequence was reconstructed from a 38,000 year-old Neandertal individual with 8341 mtDNA sequences identified among 4.8 Gb of DNA generated from approximately 0.3 g of bone. Analysis of the assembled sequence unequivocally establishes that the Neandertal mtDNA falls outside the variation of extant human mtDNAs, and allows an estimate of the divergence date between the two mtDNA lineages of 660,000 +/- 140,000 years. Of the 13 proteins encoded in the mtDNA, subunit 2 of cytochrome c oxidase of the mitochondrial electron transport chain has experienced the largest number of amino acid substitutions in human ancestors since the separation from Neandertals. There is evidence that purifying selection in the Neandertal mtDNA was reduced compared with other primate lineages, suggesting that the effective population size of Neandertals was small.     

Mitochondrial nucleoids maintain genetic autonomy but allow for functional complementation

Mitochondrial DNA (mtDNA) is packaged into DNA-protein assemblies called nucleoids, but the mode of mtDNA propagation via the nucleoid remains controversial. Two mechanisms have been proposed: nucleoids may consistently maintain their mtDNA content faithfully, or nucleoids may exchange mtDNAs dynamically. To test these models directly, two cell lines were fused, each homoplasmic for a partially deleted mtDNA in which the deletions were nonoverlapping and each deficient in mitochondrial protein synthesis, thus allowing the first unequivocal visualization of two mtDNAs at the nucleoid level. The two mtDNAs transcomplemented to restore mitochondrial protein synthesis but were consistently maintained in discrete nucleoids that did not intermix stably. These results indicate that mitochondrial nucleoids tightly regulate their genetic content rather than freely exchanging mtDNAs. This genetic autonomy provides a molecular mechanism to explain patterns of mitochondrial genetic inheritance, in addition to facilitating therapeutic methods to eliminate deleterious mtDNA mutations.     

Mitogenomic perspectives into iguanid phylogeny and biogeography: Gondwanan vicariance for the origin of Madagascan oplurines

Complete or nearly complete nucleotide sequences of mitochondrial genomes (mtDNAs) were determined from eight species which, together with previous mtDNA data for two other taxa, cover most subfamilies of Iguanidae sensu lato. These iguanid mtDNAs were found to be rather conservative with respect to gene arrangements and molecular evolutionary rates, which contrasts with mtDNAs of Acrodonta (Agamidae and Chamaeleonidae) in which several gene rearrangements and highly accelerated molecular evolutionary rates have been known. Phylogenetic analyses consistently suggested the earliest shoot-off of a Malagasy subfamily Oplurinae and an affinity of Polychrotinae and Tropidurinae sensu stricto. However, even with the ample molecular characters derived from complete mtDNA sequences, phylogenetic relationships between iguanid subfamilies were poorly resolved in general, presumably due to the rapid ancient cladogenesis. Divergence time estimation without assuming the molecular clock suggested the Late Triassic/Early Jurassic divergence of Iguanidae from acrodonts and the Middle/Late Jurassic divergence of Oplurinae from the other iguanids. Together with geological and paleontological evidence, these results led us to propose Gondwanan vicariance for the origin of Malagasy oplurines without invoking a land bridge connection between South America/Antarctica and drifting Madagascar/India.     

mtDNA and the origin of Caucasians: identification of ancient Caucasian-specific haplogroups, one of which is prone to a recurrent somatic duplication in the D-loop region.

mtDNA sequence variation was examined in 175 Caucasians from the United States and Canada by PCR amplification and high-resolution restriction-endonuclease analysis. The majority of the Caucasian mtDNAs were subsumed within four mtDNA lineages (haplogroups) defined by mutations that are rarely seen in Africans and Mongoloids. The sequence divergence of these haplogroups indicates that they arose early in Caucasian radiation and gave raise to modern European mtDNAs. Although ancient, none of these haplogroups is old enough to be compatible with a Neanderthal origin, suggesting that Homo sapiens sapiens displaced H. s. neanderthaliensis, rather than mixed with it. The mtDNAs of one of these haplogroups have a unique homoplasmic insertion between nucleotide pair (np) 573 and np 574, within the D-loop control region. This insertion makes these mtDNAs prone to a somatic mutation that duplicates a 270-bp portion of the D-loop region between np 309 and np 572. This finding suggests that certain nonpathogenic mtDNA mutations could predispose individuals to mtDNA rearrangements.     

Genetic relationships among the shads (Alosa) revealed by mitochondrial DNA analysis

We surveyed restriction site differences in mitochondrial DNA. (mtDNA) among five species of shad ( Alosa ) from North America and Europe. Allis shad, Alosa alosa and twaite shad, Alosa fallax shared two divergent genotype groups, suggesting that the two forms are either a single species, or are distinct species that have hybridized. Phenetic and cladistic analyses of the relationships among the mitochondrial genotypes defined two groups of shad, corresponding to the subgenera, Alosa and Pomolobus . The mean estimated sequence divergence between the mtDNAs of these two groups of shad was 6.5%. Taken in conjunction with fossil data, this divergence estimate suggests that the rate of mtDNA divergence between the two subgenera has been almost 10-fold lower than the 'conventional' clock calibration for mtDNA.     

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 populations and implications for detection of mitochondrial DNA mutations of pathological significance

Haplotype and phylogenetic analyses of normal mitochondrial DNAs (mtDNAs) have allowed a clarification of several controversial issues concerning the origin of humans, the time and colonization pattern of the various regions of the world, and the genetic relationships of modern human populations. More recently, the same type of analyses has also been applied to mtDNA disease studies. A review of these studies indicates that exhaustive screenings of normal mtDNA variation in all human populations associated with haplotype and phylogenetic analyses are essential if we are to understand the etiology of mitochondrial pathologies.     

Nuclear-mitochondrial epistasis and drosophila aging: introgression of Drosophila simulans mtDNA modifies longevity in D. melanogaster nuclear backgrounds

Under the mitochondrial theory of aging, physiological decline with age results from the accumulated cellular damage produced by reactive oxygen species generated during electron transport in the mitochondrion. A large body of literature has documented age-specific declines in mitochondrial function that are consistent with this theory, but relatively few studies have been able to distinguish cause from consequence in the association between mitochondrial function and aging. Since mitochondrial function is jointly encoded by mitochondrial (mtDNA) and nuclear genes, the mitochondrial genetics of aging should be controlled by variation in (1) mtDNA, (2) nuclear genes, or (3) nuclear-mtDNA interactions. The goal of this study was to assess the relative contributions of these factors in causing variation in Drosophila longevity. We compared strains of flies carrying mtDNAs with varying levels of divergence: two strains from Zimbabwe (<20 bp substitutions between mtDNAs), strains from Crete and the United States (approximately 20-40 bp substitutions between mtDNAs), and introgression strains of Drosophila melanogaster carrying mtDNA from Drosophila simulans in a D. melanogaster Oregon-R chromosomal background (>500 silent and 80 amino acid substitutions between these mtDNAs). Longevity was studied in reciprocal cross genotypes between pairs of these strains to test for cytoplasmic (mtDNA) factors affecting aging. The intrapopulation crosses between Zimbabwe strains show no difference in longevity between mtDNAs; the interpopulation crosses between Crete and the United States show subtle but significant differences in longevity; and the interspecific introgression lines showed very significant differences between mtDNAs. However, the genotypes carrying the D. simulans mtDNA were not consistently short-lived, as might be predicted from the disruption of nuclear-mitochondrial coadaptation. Rather, the interspecific mtDNA strains showed a wide range of variation that flanked the longevities seen between intraspecific mtDNAs, resulting in very significant nuclear x mtDNA epistatic interaction effects. These results suggest that even "defective" mtDNA haplotypes could extend longevity in different nuclear allelic backgrounds, which could account for the variable effects attributable to mtDNA haplogroups in human aging.     

Geographic mitochondrial DNA variation in the rock hyrax, Procavia capensis.

Restriction-fragment-length polymorphisms in mitochondrial DNA (mtDNA) were used to evaluate geographic population genetic structure in the rock hyrax, Procavia capensis, a species which occurs widely, though restricted to rocky habitat, throughout South Africa. Ten restriction endonucleases were employed to assay mtDNAs from 55 specimens representing 10 localities. Haplotypes showed strong geographic patterning, and estimates of nucleotide sequence divergence indicate two major clades thought to be dispersing along separate routes. The divergence time of approximately 2 Myr between clades is relatively high for intraspecific variation. We speculate that the marked genetic break distinguishing the northwestern populations from those constituting the south/central clade may be indicative of two species in what has conventionally been regarded as P. capensis.     

Mitochondrial DNA evidence supports northeast Indian origin of the aboriginal Andamanese in the Late Paleolithic

In view of the geographically closest location to Andaman archipelago,Myanmar was suggested to be the origin place of aboriginal Andamanese.However,for lacking any genetic information from this region,which has prevented to resolve the dispute on whether the aboriginal Andamanese were originated from mainland India or Myanmar.To solve this question and better understand the origin of the aboriginal Andamanese,we screened for haplogroups M31(from which Andaman-specific lineage M31a1 branched off)and M32 among 846mitochondrial DNAs(mtDNAs)sampled across Myanmar.As a result,two Myanmar individuals belonging to haplogroup M31 were identified,and completely sequencing the entire mtDNA genomes of both samples testified that the two M31 individuals observed in Myanmar were probably attributed to the recent gene flow from northeast India populations.Since no root lineages of haplogroup M31 or M32 were observed in Myanmar,it is unlikely that Myanmar may serve as the source place of the aboriginal Andamanese.To get further insight into the origin of this unique population,the detailed phylogenetic and phylogeographic analyses were performed by including additional 7 new entire mtDNA genomes and 113 M31 mtDNAs pinpointed from South Asian populations,and the results suggested that Andaman-specific M31a1 could in fact trace its origin to northeast India.Time estimation results further indicated that the Andaman archipelago was likely settled by modern humans from northeast India via the land-bridge which connected the Andaman archipelago and Myanmar around the Last Glacial Maximum(LGM),a scenario in well agreement with the evidence from linguistic and palaeoclimate studies.     

Phylogeny of east Asian mitochondrial DNA lineages inferred from complete sequences

The now-emerging mitochondrial DNA (mtDNA) population genomics provides information for reconstructing a well-resolved mtDNA phylogeny and for discerning the phylogenetic status of the subcontinentally specific haplogroups. Although several major East Asian mtDNA haplogroups have been identified in studies elsewhere, some of the most basal haplogroups, as well as numerous minor subhaplogroups, were not yet determined or fully characterized. To fill the lacunae, we selected 48 mtDNAs from >2,000 samples across China for complete sequencing that cover virtually all (sub)haplogroups discernible to date in East Asia. This East Asian mtDNA phylogeny can henceforth serve as a solid basis for phylogeographic analyses of mtDNAs, as well as for studies of mitochondrial diseases in East and Southeast Asia.     

No evidence of Neandertal mtDNA contribution to early modern humans

The retrieval of mitochondrial DNA (mtDNA) sequences from four Neandertal fossils from Germany, Russia, and Croatia has demonstrated that these individuals carried closely related mtDNAs that are not found among current humans. However, these results do not definitively resolve the question of a possible Neandertal contribution to the gene pool of modern humans since such a contribution might have been erased by genetic drift or by the continuous influx of modern human DNA into the Neandertal gene pool. A further concern is that if some Neandertals carried mtDNA sequences similar to contemporaneous humans, such sequences may be erroneously regarded as modern contaminations when retrieved from fossils. Here we address these issues by the analysis of 24 Neandertal and 40 early modern human remains. The biomolecular preservation of four Neandertals and of five early modern humans was good enough to suggest the preservation of DNA. All four Neandertals yielded mtDNA sequences similar to those previously determined from Neandertal individuals, whereas none of the five early modern humans contained such mtDNA sequences. In combination with current mtDNA data, this excludes any large genetic contribution by Neandertals to early modern humans, but does not rule out the possibility of a smaller contribution.     

Genome structure and gene content in protist mitochondrial DNAs.

Although the collection of completely sequenced mitochondrial genomes is expanding rapidly, only recently has a phylogenetically broad representation of mtDNA sequences from protists (mostly unicellular eukaryotes) become available. This review surveys the 23 complete protist mtDNA sequences that have been determined to date, commenting on such aspects as mitochondrial genome structure, gene content, ribosomal RNA, introns, transfer RNAs and the genetic code and phylogenetic implications. We also illustrate the utility of a comparative genomics approach to gene identification by providing evidence that orfB in plant and protist mtDNAs is the homolog of atp8 , the gene in animal and fungal mtDNA that encodes subunit 8 of the F0portion of mitochondrial ATP synthase. Although several protist mtDNAs, like those of animals and most fungi, are seen to be highly derived, others appear to be have retained a number of features of the ancestral, proto-mitochondrial genome. Some of these ancestral features are also shared with plant mtDNA, although the latter have evidently expanded considerably in size, if not in gene content, in the course of evolution. Comparative analysis of protist mtDNAs is providing a new perspective on mtDNA evolution: how the original mitochondrial genome was organized, what genes it contained, and in what ways it must have changed in different eukaryotic phyla.     

The Mitochondrial Genomes of Two Nematodes, Caenorhabditis Elegans and Ascaris Suum

The nucleotide sequences of the mitochondrial DNA (mtDNA) molecules of two nematodes, Caenorhabditis elegans [13,794 nucleotide pairs (ntp)], and Ascaris suum (14,284 ntp) are presented and compared. Each molecule contains the genes for two ribosomal RNAs (s-rRNA and l-rRNA), 22 transfer RNAs (tRNAs) and 12 proteins, all of which are transcribed in the same direction. The protein genes are the same as 12 of the 13 protein genes found in other metazoan mtDNAs: Cyt b, cytochrome b; COI-III, cytochrome c oxidase subunits I-III; ATPase6, Fo ATPase subunit 6; ND1-6 and 4L, NADH dehydrogenase subunits 1-6 and 4L: a gene for ATPase subunit 8, common to other metazoan mtDNAs, has not been identified in nematode mtDNAs. The C. elegans and A. suum mtDNA molecules both include an apparently noncoding sequence that contains runs of AT dinucleotides, and direct and inverted repeats (the AT region: 466 and 886 ntp, respectively). A second, apparently noncoding sequence in the C. elegans and A. suum mtDNA molecules (109 and 117 ntp, respectively) includes a single, hairpin-forming structure. There are only 38 and 89 other intergenic nucleotides in the C. elegans and A. suum mtDNAs, and no introns. Gene arrangements are identical in the C. elegans and A. suum mtDNA molecules except that the AT regions have different relative locations. However, the arrangement of genes in the two nematode mtDNAs differs extensively from gene arrangements in all other sequenced metazoan mtDNAs. Unusual features regarding nematode mitochondrial tRNA genes and mitochondrial protein gene initiation codons, previously described by us, are reviewed. In the C. elegans and A. suum mt-genetic codes, AGA and AGG specify serine, TGA specifies tryptophan and ATA specifies methionine. From considerations of amino acid and nucleotide sequence similarities it appears likely that the C. elegans and A. suum ancestral lines diverged close to the time of divergence of the cow and human ancestral lines, about 80 million years ago.     

Investigation of the Origin and Spread of a Mammalian Transposable Element Based on Current Sequence Diversity

Almost half the human genome consists of mobile DNA elements, and their analysis is a vital part of understanding the human genome as a whole. Many of these elements are ancient and have persisted in the genome for tens or hundreds of millions of years, providing a window into the evolution of modern mammals. The Golem family have been used as model transposons to highlight computational analyses which can be used to investigate these elements, particularly the use of molecular dating with large transposon families. Whole-genome searches found Golem sequences in 20 mammalian species. Golem A and B subsequences were only found in primates and squirrel. Interestingly, the full-length Golem, found as a few copies in many mammalian genomes, was found abundantly in horse. A phylogenetic profile suggested that Golem originated after the eutherian-metatherian divergence and that the A and B subfamilies originated at a much later date. Molecular dating based on sequence diversity suggests an early age, of 175 Mya, for the origin of the family and that the A and B lineages originated much earlier than expected from their current taxonomic distribution and have subsequently been lost in some lineages. Using publically available data, it is possible to investigate the evolutionary history of transposon families. Determining in which organisms a transposon can be found is often used to date the origin and expansion of the families. However, in this analysis, molecular dating, commonly used for determining the age of gene sequences, has been used, reducing the likelihood of errors from deleted lineages.     

Mitochondrial DNA diversity and population differentiation in southern East Asia

Mitochondrial DNA (mtDNA) polymorphism has been studied systematically in the Han, Tibeto-Buman, and Hmong-Mien ethnic families of southern East Asia. Only two families in this region, Daic and Austro-Asiatic, were still uninvestigated. Daic is a major ethnic family in South China and Southeast Asia and has a long history. To study mtDNA polymorphism within this family, all the Daic populations of China and some of Vietnam (774 individuals from 30 populations) were typed by HVS-1 region sequencing and by PCR-RFLP assays. The observed high Southern type frequencies (B, F, M7, R) confirmed Daic as a typical Southern group. mtDNAs of other populations (126 individuals from 14 populations) from Austro-Asiatic ethnic families neighboring the Daic were also typed. Networks of mtDNA haplogroups in South China were traced from these new data and those from the literature. Ethnic families share many haplogroups, indicating their common origin. However, the two largest families in South China, Daic, and Hmong-Mien, polarized into several ethnic family specific haplogroups. Haplogroup ages were estimated in the networks of high-frequency haplogroups (B, F, M7, R), and they were found to originate about 50,000 years ago. In contrast, ethnic family specific haplogroups all originated around 20,000 years ago. We therefore conclude that modern humans have lived in South China for a long time, inside-ethnogenesis was a rather late event, and frequent inmixing was taking place throughout. MtDNA data of Daic, Austro-Asiatic and other populations in South China has therefore proven pivotal for studying the human history of East Asia.