Mitochondrial diversity and phylogeographic structure of Chinese domestic goats

China has numerous native domestic goat breeds, but so far there has been no extensive study on genetic diversity, population demographic history, and origin of Chinese goats. Here, we examined the genetic diversity and phylogeographic structure of Chinese domestic goats by determining a 481-bp fragment of the first hypervariable region of mitochondrial DNA (mtDNA) control region from 368 individuals representing 18 indigenous breeds. Phylogenetic analyses revealed that there were four mtDNA lineages (A-D) identified in Chinese goats, in which lineage A was predominant, lineage B was moderate, and lineages C and D were at low frequency. These results further support the multiple maternal origins of domestic goats. The pattern of genetic variation in goat mtDNA sequences indicated that the two larger lineages A and B had undergone population expansion events. In a combined analysis of previously reported sequences and our sequences belonging to lineage B, we detected two subclades, in which one was unique to eastern Asia and another was shared between eastern and southern Asia. A larger genetic variation in eastern Asia than southern Asia and the pattern of phylogeographic variation in lineage B suggest that at least one subclade of lineage B originated from eastern Asia. There was no significant geographical structuring in Chinese goat populations, which suggested that there existed strong gene flow among goat populations caused by extensive transportation of goats in history.     

Variation in mitochondrial COII gene sequences among two species of Japanese knotweed-boring moths, Ostrinia latipennis and O. ovalipennis (Lepidoptera: Crambidae)

The Ostrinia latipennis group contains two species, O. latipennis (Warren) and O. ovalipennis Ohno. These two species commonly utilize perennial knotweeds (Fallopia spp.) as their host plants, which are serious invasive weeds in Europe and North America. Ostrinia latipennis is widely distributed across north-east Asia including Japan whereas O. ovalipennis is restricted to north Japan (Hokkaido Is.) and highland areas of central Japan (Nagano Prefecture in Honshu Is.). To estimate the phylogenetic relatedness and geographical differentiation of the two species, mitochondrial COII gene sequences were determined for specimens covering their distribution ranges in Japan. The uncorrected sequence divergence between O. latipennis and O. ovalipennis was 0.6-0.7%, supporting a close relationship. According to the standard molecular clock proposed for arthropod mtDNA, the two species are speculated to have diverged about 0.3 Myr ago. A single COII gene haplotype was found in O. latipennis irrespective of collection locality. In contrast, two haplotypes were found in O. ovalipennis, and their frequencies were significantly different between the Hokkaido and Honshu populations. The patterns of geographical variation in the COII gene within the two species were in agreement with previously reported patterns of geographical differentiation in morphology of the two species in Japan. The present results support the hypothesis that gene flow among local populations of O. ovalipennis has been limited by geographical isolation.     

The origins of the Polynesians: an interpretation from mitochondrial lineage analysis.

Using mitochondrial lineage analysis of 1,178 individuals from Polynesia, the western Pacific, and Taiwan, we show that the major prehistoric settlement of Polynesia was from the west and involved two or possibly three genetically distinct populations. The predominant lineage group, accounting for 94% of Polynesian mtDNA, shares a 9-bp COII/tRNA(Lys) intergenic deletion and characteristic control region transition variants, compared to the Cambridge reference sequence. In Polynesia, the diversity of this group is extremely restricted, while related lineages in Indonesia, the Philippines, and Taiwan are increasingly diverse. This suggests a relatively recent major eastward expansion into Polynesia, perhaps originating from Taiwan, in agreement with archeological and linguistic evidence, but which experienced one or more severe population bottlenecks. The second mitochondrial lineage group, accounting for 3.5% of Polynesian mtDNA haplotypes, does not have the 9-bp deletion and its characterized by an A-C transversional variant at nt position 16265. Specific oligonucleotides for this variant were used to select individuals from the population sample who, with other sequences, show that the Polynesian lineages were part of a diverse group in Vanuatu and Papua New Guinea. The very low overall diversity of both lineage groups in Polynesia suggests there was severe population restriction during the colonization of remote Oceania. A third group, represented by only four individuals (0.6%) in Polynesia but also present in the Philippines, shares variants at nt positions 16172 and 16304. Two Polynesians had unrelated haplotypes matching published sequences from native South Americans, which may be the first genetic evidence of prehistoric human contact between Polynesia and South America.     

Founding Amerindian mitochondrial DNA lineages in ancient Maya from Xcaret, Quintana Roo

Ancient DNA from the bone remains of 25 out of 28 pre-Columbian individuals from the Late Classic-Postclassic Maya site of Xcaret, Quintana Roo, was recovered, and mitochondrial DNA (mtDNA) was amplified by using the polymerase chain reaction. The presence of the four founding Amerindian mtDNA lineages was investigated by restriction analysis and by direct sequencing in selected individuals. The mtDNA lineages A, B, and C were found in this population. Eighty-four percent of the individuals were lineage A, whereas lineages B and C were present at low frequencies, 4% and 8%, respectively. Lineage D was absent from our sample. One individual did not possess any of the four lineages. Six skeletons out of 7 dated from the Late Classic period were haplotype A, whereas 11 skeletons out of 16 dated from the Postclassic period were also haplotype A. The distribution of mtDNA lineages in the Xcaret population contrasts sharply with that found in ancient Maya from Copán, which lack lineages A and B. On the other hand, our results resemble more closely the frequencies of mtDNA lineages found in contemporary Maya from the Yucatán Peninsula and in other Native American contemporary populations of Mesoamerican origin.     

Spatial structure of lemming populations (Dicrostonyx groenlandicus) fluctuating in density

The pattern and scale of the genetic structure of populations provides valuable information for the understanding of the spatial ecology of populations, including the spatial aspects of density fluctuations. In the present paper, the genetic structure of periodically fluctuating lemmings (Dicrostonyx groenlandicus) in the Canadian Arctic was analysed using mitochondrial DNA (mtDNA) control region sequences and four nuclear microsatellite loci. Low genetic variability was found in mtDNA, while microsatellite loci were highly variable in all localities, including localities on isolated small islands. For both genetic markers the genetic differentiation was clear among geographical regions but weaker among localities within regions. Such a pattern implies gene flow within regions. Based on theoretical calculations and population census data from a snap-trapping survey, we argue that the observed genetic variability on small islands and the low level of differentiation among these islands cannot be explained without invoking long distance dispersal of lemmings over the sea ice. Such dispersal is unlikely to occur only during population density peaks.     

Mitochondrial sequence variation in the Guahibo Amerindian population from Venezuela

New data were obtained on mitochondrial DNA (mtDNA) from Guahibo from Venezuela, a group so far not studied using molecular data. A population sample (n = 59) was analyzed for mtDNA variation in two control-region hypervariable segments (HV1 and HV2) by sequencing. The presence or absence of a 9-bp polymorphism in the COII/tRNA(Lys) region was studied by direct amplification and electrophoretic identification. Thirty-eight variable sites were detected in regions HV1 and HV2, defining 26 mtDNA lineages; 23.7% of these were present in a single individual. The 9-bp deletion was found in 3.39% of individuals. Nucleotide and haplotype diversities were relatively high compared with other New World populations. The identified sequence haplotypes were classified into four major haplogroups (A-D) according to previous studies, with high frequencies for A (47.46%) and C (49.15%), low frequency for B (3.39%), and an absence of D.     

Limits of Bayesian skyline plot analysis of mtDNA sequences to infer historical demographies in Pacific herring (and other species)

A previous analysis of Pacific herring mitochondrial (mt) DNA with Bayesian skyline plots (BSPs) was interpreted to reflect population growth in the late Pleistocene that was preceded by population stability over several hundred thousand years. Here we use an independent set of mtDNA control region (CR) sequences and simulations to test these hypotheses. The CR haplotype genealogy shows three deeply divided lineages, A, B and C, with divergences ranging from d=1.6% to 1.9% and with similar genetic diversities (h=0.95, 0.96, 0.94; Θ(π)=0.011, 0.012, 0.014, respectively). Lineage A occurs almost exclusively in the NW Pacific and Bering Sea, but lineages B and C are co-distributed in the Northeastern Pacific. This distribution points to a historical allopatric separation between A and B-C across the North Pacific during Pleistocene glaciations. The origins of B and C are uncertain. One hypothesis invokes long-term isolation of lineage C in the Sea of Cortez, but the present-day lack of geographical segregation from lineage B argues for lineage sorting to explain the deep divergence between B and C. BSPs depict rapid population growth in each lineage, but the timing of this growth is uncertain, because of questions about an appropriate molecular clock calibration. We simulated historical demographies under a Pleistocene climate model using observed genetic parameters. BSPs for these sequences showed rapid population growth after the Last Glacial Maximum (LGM) 18-20kyears ago and a flat population history during previous climate fluctuations. Population declines during the LGM appear to have erased signals of previous population fluctuations.     

Mitochondrial DNA phylogeny of the Old-World monkey tribe Papionini.

The evolution of the Old World monkey tribe Papionini, composed of macaques, baboons, mandrills, drills, and mangabeys, was examined using mitochondrial DNA (mtDNA) sequence data on the cytochrome oxidase subunit II gene. When analyzed cladistically, these data support a baboon clade of savannah (Papio) plus gelada (Theropithecus) baboons, as well as a clade containing drill (Mandrillus) plus mangabey (Cerocebus) genera. This result stands in opposition to most morphological phylogenies, which break up the baboon clade by placing Papio and Mandrillus as sister taxa and Theropithecus as a more distantly related lineage. Analyses of COII gene sequences also suggest that the papionin ancestral stock divided into two lineages, one leading to macaques and the other to the purely African genera. From a molecular evolutionary perspective, the papionin COII gene sequences reveal a pattern of amino acid replacements concentrated in the regions spanning the mitochondrial membrane.     

Genetic differentiation and cryptic speciation in natural populations of Drosophila lacertosa

Drosophila lacertosa, an Oriental member of the robusta species group in the virilis-repleta radiation, has a wide distribution from northern India throughout China to the Far East. Phylogenetic analyses of mitochondrial ND2 gene sequences revealed two genetically significantly diverged lineages with 2.1% DNA sequence divergence. These two lineages are largely allopatric: one is mainly found in southwest China (lineage A), whereas the other ranges central and east of China to northeast Asia (lineage B). The geographic distributions of these two clades narrowly overlap near the provincial boundaries of Yunnan and Guizhou or Guangxi, probably as a consequence of population expansion and secondary contact. These two lineages were estimated to have diverged about 0.56 Mya, and the estimated time of expansion was approximately 70,000 years ago for lineage A population, and 110,000 years ago for lineage B population. The present molecular analysis, together with nearly complete reproductive isolation between those two lineages, strongly suggests that the two lineages might be incipient species.     

Mitochondrial DNA diversity in an apomictic Daphnia complex from the Canadian High Arctic

Cyclic parthenogenesis is the ancestral mode of reproduction in the cladoceran crustacean, Daphnia pulex, but some populations have made the transition to obligate parthenogenesis and this is the only mode of reproduction known to occur in arctic populations. Melanism and polyploidy are also common in arctic populations of this species. Prior allozyme studies of arctic D. pulex revealed substantial levels of clonal diversity on a regional scale. Clonal groupings based on cluster analysis of allozyme genotypes do not conform to groupings based on the presence/absence of melanin or on ploidy level. In order to further elucidate genetic relationships among arctic D. pulex clones, mitochondrial DNA (mtDNA) variation was examined in 31 populations from two Canadian high-arctic sites. The data were also compared to a previous study of mtDNA variation in populations from a Canadian low-arctic site. Cladistic analysis of restriction site variation of the entire mitochondrial genome and nucleotide sequence variation of the mitochondrial control region was used to construct genetic relationships among mitochondrial genotypes. Three distinct mitochondrial lineages were detected. One lineage was associated with diploid, nonmelanic clones and is the same as the lineage that is found in temperate populations of D. pulex. The other two lineages (A & B) were associated with polyploid, melanic clones. Sequence divergence between the A and B lineages was 2.4%. Sequence divergence between D. pulex and either of these two lineages exceeded 3%. It is suggested that the melanic, polyploid clones are hybrids between males of D. pulex (and/or a closely related congener, D. pulicaria) and females of either of two ancestral melanic species that have mitochondrial lineages A and B. Geographic patterns of mitochondrial diversity in ‘melanic’ lineage B support the hypothesis of an high-arctic refuge for the ancestral species during the last glacial period.     

American Indian prehistory as written in the mitochondrial DNA: a review.

Native Americans have been divided into three linguistic groups: the reasonably well-defined Eskaleut and Nadene of northern North America and the highly heterogeneous Amerind of North, Central, and South America. The heterogeneity of the Amerinds has been proposed to be the result of either multiple independent migrations or a single ancient migration with extensive in situ radiation. To investigate the origin and interrelationship of the American Indians, we examined the mitochondrial DNA (mtDNA) variation in 87 Amerinds (Pima, Maya, and Ticuna of North, Central, and South America, respectively), 80 Nadene (Dogrib and Tlingit of northwest North America and Navajo of the southwest North America), and 153 Asians from 7 diverse populations. American Indian mtDNAs were found to be directly descended from five founding Asian mtDNAs and to cluster into four lineages, each characterized by a different rare Asian mtDNA marker. Lineage A is defined by a HaeIII site gain at np 663, lineage B by a 9-bp deletion between the COII and tRNA(Lys) genes, lineage C by a HincII site loss at np 13259, and lineage D by an AluI site loss at np 5176. The North, Central, and South America Amerinds were found to harbor all four lineages, demonstrating that the Amerinds originated from a common ancestral genetic stock. The genetic variation of three of the four Amerind lineages (A, C, and D) was similar with a mean value of 0.084%, whereas the sequence variation in the fourth lineage (B) was much lower, raising the possibility of an independent arrival. By contrast, the Nadene mtDNAs were predominantly from lineage A, with 27% of them having a Nadene-specific RsaI site loss at np 16329. The accumulated Nadene variation was only 0.021%. These results demonstrate that the Amerind mtDNAs arose from one or maybe two Asian migrations that were distinct from the migration of the Nadene and that the Amerind populations are about four times older than the Nadene.     

Independent maternal origin of Chinese swamp buffalo (Bubalus bubalis)

To obtain more knowledge on the origin and genetic diversity of the swamp buffalo (Bubalus bubalis) in China, the complete mitochondrial D-loop sequences of 119 samples representing seven native types were compared. Two mitochondrial DNA (mtDNA) lineages (lineages A and B) were determined for the Chinese swamp buffalo. Examination of the diversity patterns suggest that lineage A has undergone a population expansion event. Divergence of lineages A and B was estimated at 18,000 years ago. Combined analyses of mtDNA sequences from Chinese, Indian, Brazilian/Italian and Southeast Asian/Australian buffalo samples showed independent domestication events in the swamp buffalo from China and the river buffalo from the India subcontinent. The spread of swamp and river buffalo from China and India respectively to mainland Southeast Asia suggests that Southeast Asia is a hybrid zone for buffalo. Our data support the hypothesis of the evolution of domesticated swamp and river buffalo from ancestral swamp-like animals. These ancestral animals were extensively distributed across mainland Asia and most likely are represented today by the wild Asian buffalo (Bubalus arnee).     

Incongruence in the pattern and timing of intra-specific diversification in bronze frogs and bullfrogs (Ranidae)

We compare patterns of lineage divergence in mitochondrial DNA (mtDNA) sequences of two protein-encoding mitochondrial genes (cyt b and ND2) in two ecologically similar, co-distributed, and closely related ranid frogs (Rana clamitans and Rana catesbeiana), that are geographically widespread, and frequently syntopic. We identified three lineages in R. clamitans, separated by 0.5% to 2.1% net corrected sequence divergence, comparable to two R. catesbeiana lineages separated by 0.6%. The geographic pattern of lineage distribution differed notably between the two species. In R. clamitans, we found a Coastal Plain-Appalachian (CPA) lineage restricted to south and east of the Appalachian Mountains and a widespread lineage that encompassing nearly all the sampled range. A third distinct and divergent lineage was detected in one location in the southwest portion of the range (Louisiana). This pattern contrasts with the east-west pattern in R. catesbeiana, and reflects possible differences in refugial dynamics and patterns of range expansion. Although both species have undergone range expansion and population growth, coalescent reconstruction of N(e) reflects larger lineages but more recent divergence in R. clamitans relative to R. catesbeiana, reflecting significant differences in population history or divergent patterns of molecular evolution at mtDNA.     

Behavioural Reproductive Isolation in a Rotifer Hybrid Zone

A hybrid zone between two Brachionus plicatilis rotifer mitochondrial DNA (mtDNA) lineages was recently described in the Iberian Peninsula between a pond (Santed 2) and a lake (Gallocanta). The patterns of mitochondrial and nuclear genetic variation observed suggested that gene flow is mainly male-mediated from the lake to the pond. Here we test two hypotheses: (a) that male-mediated gene flow occurs through assortative mating between individuals from these ponds, (b) that behavioural isolation occurs between the two mtDNA lineages. We isolated, reared and genotyped rotifer clones from resting eggs collected in the sediments of these and two other distant ponds. We devised a quick, inexpensive RFLP method to discriminate between B. plicatilis and its sibling species B. ‘Manjavaeas’ and between both mtDNA B. plicatilis lineages. Behavioural no-choice tests using new-born, virgin males and females were performed between five clones. B. ‘Manjavacas’ and B. plicatilis were reproductively isolated. B. plicatilis clones did not show evidence of reproductive isolation, regardless of their mtDNA lineage, except Santed 2 males, which discriminated strongly against Gallocanta females. These results could help to explain the discrepancies between mitochondrial and nuclear genetic variation reported in the two populations.     

Comparative phylogeography of the two pink salmon broodlines: an analysis based on a mitochondrial DNA genealogy

Over most of their natural northern Pacific Ocean range, pink salmon (Oncorhynchus gorbuscha) spawn in a habitat that was repeatedly and profoundly affected by Pleistocene glacial advances. A strictly two-year life cycle of pink salmon has resulted in two reproductively isolated broodlines, which spawn in alternating years and evolved as temporal replicates of the same species. To study the influence of historical events on phylogeographical and population genetic structure of the two broodlines, we first reconstructed a fine-scale mtDNA haplotype genealogy from a sample of 80 individuals and then determined the geographical distribution of the major genealogical assemblages for 718 individuals sampled from nine Alaskan and eastern Asian even- and nine odd-year pink salmon populations. Analysis of restriction site states in seven polymerase chain reaction (PCR)-amplified mtDNA regions (comprising 97% of the mitochondrial genome) using 13 endonucleases resolved 38 haplotypes, which clustered into five genealogical lineages that differed from 0.065 to 0.225% in net sequence divergence. The lineage sorting between broodlines was incomplete, which suggests a recent common ancestry. Within each lineage, haplotypes exhibited star-like genealogies indicating recent population growth. The depth of the haplotype genealogy is shallow ( approximately 0.5% of nucleotide sequence divergence) and probably reflects repeated decreases in population size due to Pleistocene glacial advances. Nested clade analysis (NCA) of geographical distances showed that the geographical distribution observed for mitochondrial DNA (mtDNA) haplotypes resulted from alternating influences of historical range expansions and episodes of restricted dispersal. Analyses of molecular variance showed weak geographical structuring of mtDNA variation, except for the strong subdivision between Asian and Alaskan populations within the even-year broodline. The genetic similarities observed among and within geographical regions probably originated from postglacial recolonizations from common sources rather than extensive gene flow. The phylogeographical and population genetic structures differ substantally between broodlines. This can be explained by stochastic lineage sorting in glacial refugia and perhaps different recolonization routes in even- and odd-year broodlines.     

Extensive Variation and Sub-Structuring in Lineage A mtDNA in Indian Sheep: Genetic Evidence for Domestication of Sheep in India

Previous studies on mitochondrial DNA analysis of sheep from different regions of the world have revealed the presence of two major- A and B, and three minor- C, D and E maternal lineages. Lineage A is more frequent in Asia and lineage B is more abundant in regions other than Asia. We have analyzed mitochondrial DNA sequences of 330 sheep from 12 different breeds of India. Neighbor-joining analysis revealed lineage A, B and C in Indian sheep. Surprisingly, multidimensional scaling plot based on F_ST values of control region of mtDNA sequences showed significant breed differentiation in contrast to poor geographical structuring reported earlier in this species. The breed differentiation in Indian sheep was essentially due to variable contribution of two major lineages to different breeds, and sub- structuring of lineage A, possibly the latter resulting from genetic drift. Nucleotide diversity of this lineage was higher in Indian sheep (0.014 ± 0.007) as compared to that of sheep from other regions of the world (0.009 ± 0.005 to 0.01 ± 0.005). Reduced median network analysis of control region and cytochrome b gene sequences of Indian sheep when analyzed along with available published sequences of sheep from other regions of the world showed that several haplotypes of lineage A were exclusive to Indian sheep. Given the high nucleotide diversity in Indian sheep and the poor sharing of lineage A haplotypes between Indian and non-Indian sheep, we propose that lineage A sheep has also been domesticated in the east of Near East, possibly in Indian sub-continent. Finally, our data provide support that lineage B and additional lineage A haplotypes of sheep might have been introduced to Indian sub-continent from Near East, probably by ancient sea trade route.     

Molecular systematics and phylogeography of Cebus capucinus (Cebidae, Primates) in Colombia and Costa Rica by means of the mitochondrial COII gene

We propose the first molecular systematic hypothesis for the origin and evolution of Cebus capucinus based on an analysis of 710 base pairs (bp) of the cytochrome c oxidase subunit II (COII) mitochondrial gene in 121 C. capucinus specimens sampled in the wild. The animals came from the borders of Guatemala and Belize, Costa Rica, and eight different departments of Colombia (Antioquia, Chocó, Sucre, Bolivar, Córdoba, Magdalena, Cauca, and Valle del Cauca). Three different and significant haplotype lineages were found in Colombia living sympatrically in the same departments. They all presented high levels of gene diversity but the third Colombian gene pool was determined likely to be the most ancestral lineage. The second Colombian mitochondrial (mt) haplogroup is likely the source of origin of the unique Central America mt haplogroup that was detected. Our molecular population genetics data do not agree with the existence of two well-defined subspecies in Central America (limitaneus and imitator). This Central America mt haplogroup showed significantly less genetic diversity than the Colombian mt haplogroups. All the C. capucinus analyzed showed evidence of historical population expansions. The temporal splits among these four C. capucinus lineages were related to the completion of the Panamanian land bridge as well as to climatic changes during the Quaternary Period.     

Ancient sex-specific extension of the cytochrome c oxidase II gene in bivalves and the fidelity of doubly-uniparental inheritance

Bivalves of the families Mytilidae and Unionidae show a unique mode of mitochondrial DNA inheritance called doubly uniparental inheritance. In addition to receiving the maternally transmitted mtDNA lineage, males receive a separate mtDNA genome from their fathers. This system is sometimes compromised, in that female genomes are occasionally recruited into the male cycle of inheritance. These masculinization events are common in the Mytilidae but have not been reported in the Unionidae. In order to estimate the age of the male and the female lineages in the Unionidae and to look for evidence of masculinization, we sequenced the junction between the cytochrome c oxidase II gene and the cytochrome c oxidase I gene. The unionid male and female lineages diverged approximately 450 MYA. There is no evidence for masculinization during this period, suggesting that there are taxon-specific differences in the rate of masculinization. Coincidentally, a 200-codon extension of the COII gene is present in the male genome of the Unionidae and may be responsible for the absence of masculinization.     

Nuclear loci and coalescent methods support ancient hybridization as cause of mitochondrial paraphyly between gadwall and falcated duck (Anas spp.)

Many species have mitochondrial DNA lineages that are phylogenetically intermixed with other species, but studies have rarely tested the cause of such paraphyly. In this study, we tested two hypotheses that could explain mitochondrial paraphyly of Holarctic gadwalls (Anas strepera) with respect to Asian falcated ducks (A. falcata). First, hybridization could have resulted in falcated duck mitochondrial DNA (mtDNA) introgressing into the gadwall gene pool. Second, gadwalls and falcated ducks could have diverged so recently that mtDNA lineages have not sorted to reciprocal monophyly. We used coalescent analyses of three independent loci to distinguish between these two hypotheses. Two lines of evidence support introgression. First, analyses of the three loci combined show that some introgression is necessary to explain current genetic diversity in gadwalls. Second, we generated alternative predictions regarding time since divergence estimated from mtDNA: falcated ducks and gadwalls would have diverged between 65,000 and 700,000 years before present (ybp) under the introgression hypothesis and between 11,000 and 76,000 ybp under the incomplete lineage sorting hypothesis. The two independent nuclear introns indicated that these species diverged between 210,000 and 5,200,000 ybp, which did not overlap the predicted time for incomplete lineage sorting. These analyses also suggested that ancient introgression ( approximately 14,000 ybp) has resulted in the widespread distribution and high frequency of falcated-like mtDNA (5.5% of haplotypes) in North America. This is the first study to use a rigorous quantitative framework to reject incomplete lineage sorting as the cause of mitochondrial paraphyly.