Substitution bias, rapid saturation, and the use of mtDNA for nematode systematics

Only relatively recently have researchers turned to molecular methods for nematode phylogeny reconstruction. Thus, we lack the extensive literature on evolutionary patterns and phylogenetic usefulness of different DNA regions for nematodes that exists for other taxa. Here, we examine the usefulness of mtDNA for nematode phylogeny reconstruction and provide data that can be used for a priori character weighting or for parameter specification in models of sequence evolution. We estimated the substitution pattern for the mitochondrial ND4 gene from intraspecific comparisons in four species of parasitic nematodes from the family Trichostrongylidae (38-50 sequences per species). The resulting pattern suggests a strong mutational bias toward A and T, and a lower transition/transversion ratio than is typically observed in other taxa. We also present information on the relative rates of substitution at first, second, and third codon positions and on relative rates of saturation of different types of substitutions in comparisons ranging from intraspecific to interordinal. Silent sites saturate extremely quickly, presumably owing to the substitution bias and, perhaps, to an accelerated mutation rate. Results emphasize the importance of using only the most closely related sequences in order to infer patterns of substitution accurately for nematodes or for other taxa having strongly composition-biased DNA. ND4 also shows high amino acid polymorphism at both the intra- and interspecific levels, and in higher level comparisons, there is evidence of saturation at variable amino acid sites. In general, we recommend using mtDNA coding genes only for phylogenetics of relatively closely related nematode species and, even then, using only nonsynonymous substitutions and the more conserved mitochondrial genes (e.g., cytochrome oxidases). On the other hand, the high substitution rate in genes such as ND4 should make them excellent for population genetics studies, identifying cryptic species, and resolving relationships among closely related congeners when other markers show insufficient variation.      

Nuclear insertions of mitochondrial origin: Database updating and usefulness in cancer studies

Nuclear insertions of mitochondrial origin (NUMTs) can be useful tools in evolution and population studies. However, due to their similarity to mitochondrial DNA (mtDNA), NUMTs may also be a source of contamination in mtDNA studies. The main goal of this work is to present a database of NUMTs, based on the latest version of the human genome—GRCh37 draft. A total of 755 insertions were identified. There are 33 paralogous sequences with over 80% sequence similarity and of a greater length than 500 bp. The non-identical positions between paralogous sequences are listed for the first time. As an application example, the described database is used to evaluate the impact of NUMT contamination in cancer studies. The evaluation reveals that 220 positions from 256 with zero hits in the current mtDNA phylogeny could in fact be traced to one or more nuclear insertions of mtDNA. This is due to they are located in non-identical positions between mtDNA and nuclear DNA (nDNA). After in silico primer validation of each revised cancer study, risk of co-amplification between mtDNA and nDNA was detected in some cases, whereas in others no risk of amplification was identified. This approach to cancer studies clearly proves the potential of our NUMT database as a valuable new tool to validate mtDNA mutations described in different contexts. Moreover, due to the amount of information provided for each nuclear insertion, this database should play an important role in designing evolutionary, phylogenetic and epidemiological studies.     

Molecular phylogenetic systematics of twelve species of Acipenseriformes based on mtDNA ND4L -ND4 gene sequence analysis

Acipenseriformes is an endangered primitive fish group, which occupies a special place in the history of ideas concerning fish evolution, even in vertebrate evolution. However, the classification and evolution of the fishes have been debated. The mitochondrial DNA (mtDNA) ND4L and partial ND4 genes were first sequenced in twelve species of the order Acipenseriformes, including endemic Chinese species. The following points were drawn from DNA sequences analysis: (i) the two species of Huso can be ascribed to Acipenser; (ii) A. dabryanus is the mostly closely related to A. sinensis, and most likely the landlocked form of A. sinensis; (iii) genus Acipenser in trans-Pacific region might have a common origin; (iv) mtDNA ND4L and ND4 genes are the ideal genetic markers for phylogenetic analysis of the order Acipenseriformes.     

云南猕猴mtDNA控制区全序列测定

目的测定云南猕猴线粒体DNA控制区全序列,对其进行鉴定及进化分析。方法利用PCR技术扩增猕猴线粒体DNA控制区全序列,结合GenBank中下载的猕猴参考序列（AY612638）,采用多个生物学软件对序列碱基组成、同源性、转换／颠换比等遗传信息进行分析,并基于邻接法（NJ）和最小进化法（ME）构建系统进化树。结果云南猕猴线粒体DNA控制区全长为（1084-1089）bp,A、T、G和c四种碱基平均含量分别为29．9％、26．9％、12．3％和30．9％,A＋T含量（56．8％）高于G＋C含量（43．2％）。所分析序列间的同源性为91．5％-99．5％,平均核苷酸变异率为4．5％,变异类型包括转换、颠换、插入和缺失4种形式,转换／颠换比值平均为26．1。进化树显示云南猕猴存在两个平行进化的姐妹分支。结论本研究获得了云南猕猴mtDNA控制区全序列,为猕猴进化关系研究及mtDNA控制区功能研究奠定基础。     

Rates of nucleotide substitution and mammalian nuclear gene evolution. Approximate and maximum-likelihood methods lead to different conclusions

Rates and patterns of synonymous and nonsynonymous substitutions have important implications for the origin and maintenance of mammalian isochores and the effectiveness of selection at synonymous sites. Previous studies of mammalian nuclear genes largely employed approximate methods to estimate rates of nonsynonymous and synonymous substitutions. Because these methods did not account for major features of DNA sequence evolution such as transition/transversion rate bias and unequal codon usage, they might not have produced reliable results. To evaluate the impact of the estimation method, we analyzed a sample of 82 nuclear genes from the mammalian orders Artiodactyla, Primates, and Rodentia using both approximate and maximum-likelihood methods. Maximum-likelihood analysis indicated that synonymous substitution rates were positively correlated with GC content at the third codon positions, but independent of nonsynonymous substitution rates. Approximate methods, however, indicated that synonymous substitution rates were independent of GC content at the third codon positions, but were positively correlated with nonsynonymous rates. Failure to properly account for transition/transversion rate bias and unequal codon usage appears to have caused substantial biases in approximate estimates of substitution rates.     

Mitochondrial DNA-like sequences in the nuclear genome of the opossum genus Didelphis (Marsupialia: Didelphidae).

This study reports the occurrence of mtDNA-like sequences in the nuclear genome of the opossum genus Didelphis (Didelphidae, Marsupialia). A specific primer pair designed to amplify a region encompassing a 3' terminal 118 bp region of the cytochrome b gene, the Thr and Pro tRNA genes, and a 489 bp region of the D-loop of the D. virginiana mtDNA, was used in highly stringent PCR reactions. These PCR reactions resulted in several fragments per individual varying in size from 259 bp to 1 kb. The sequencing of some of these fragments showed the occurrence of paralogous mtDNA-like sequences among the PCR amplified fragments. Analyses of qualitative aspects of these sequences, their transition/transversion ratios, and phylogenetic relationships were conclusive in showing the occurrence of mtDNA-like sequences in the nuclear genome of the genus Didelphis. Comparisons and phylogenetic analysis of orthologous mtDNA from the four Didelphis species and paralogous nuclear sequences suggested that mtDNA migration to the nuclear genome occurred more than once in Didelphis evolution.     

Sequence evolution of mitochondrial tRNA genes and deep-branch animal phylogenetics

Mitochondrial DNA sequences are often used to construct molecular phylogenetic trees among closely related animals. In order to examine the usefulness of mtDNA sequences for deep-branch phylogenetics, genes in previously reported mtDNA sequences were analyzed among several animals that diverged 20–600 million years ago. Unambiguous alignment was achieved for stem-forming regions of mitochondrial tRNA genes by virtue of their conservative secondary structures. Sequences derived from stem parts of the mitochondrial tRNA genes appeared to accumulate much variation linearly for a long period of time: nearly 100 Myr for transition differences and more than 350 Myr for transversion differences. This characteristic could be attributed, in part, to the structural variability of mitochondrial tRNAs, which have fewer restrictions on their tertiary structure than do nonmitochondrial tRNAs. The tRNA sequence data served to reconstruct a well-established phylogeny of the animals with 100% bootstrap probabilities by both maximum parsimony and neighbor joining methods. By contrast, mitochondrial protein genes coding for cytochrome b and cytochrome oxidase subunit I did not reconstruct the established phylogeny or did so only weakly, although a variety of fractions of the protein gene sequences were subjected to tree-building. This discouraging phylogenetic performance of mitochondrial protein genes, especially with respect to branches originating over 300 Myr ago, was not simply due to high randomness in the data. It may have been due to the relative susceptibility of the protein genes to natural selection as compared with the stem parts of mitochondrial tRNA genes. On the basis of these results, it is proposed that mitochondrial tRNA genes may be useful in resolving deep branches in animal phylogenies with divergences that occurred some hundreds of Myr ago. For this purpose, we designed a set of primers with which mtDNA fragments encompassing clustered tRNA genes were successfully amplified from various vertebrates by the polymerase chain reaction.Abbreviations AA stem amino acid-acceptor stem - AC stem anticodon stem - COI cytochrome oxidase subunit I - cytb cytochrome b - D stem dihydrouridine stem - MP maximum parsimony - mtDNA mitochondrial DNA - Myr million years - NJ neighbor joining - PCR polymerase chain reaction - Ti transition - T stem tC stem - Tv transversion Correspondence to: Y. Kumazawa     

Lineage-specific selection in human mtDNA: lack of polymorphisms in a segment of MTND5 gene in haplogroup J

Human mitochondrial DNA (mtDNA) is a nonrecombining genome that codes for 13 subunits of the mitochondrial oxidative phosphorylation system, 2 rRNAs, and 22 tRNAs. Mutations have accumulated sequentially in mtDNA lineages that diverged tens of thousands of years ago. The genes in mtDNA are subject to different functional constraints and are therefore expected to evolve at different rates, but the rank order of these rates should be the same in all lineages of a phylogeny. Previous studies have indicated, however, that specific regions of mtDNA may have experienced different histories of selection in different lineages, possibly because of lineage-specific interactions or environmental factors such as climate. We report here on a survey for lineage-specific patterns of nucleotide polymorphism in human mtDNA. We calculated molecular polymorphism indices and neutrality tests for classes of functional sites and genes in 837 human mtDNA sequences, compared the results between continent-specific mtDNA lineages, and used two sliding window methods to identify differences in the patterns of polymorphism between haplogroups. A general correlation between nucleotide position and the level of nucleotide polymorphism was identified in the coding region of the mitochondrial genome. Nucleotide diversity in the protein-coding sequence of mtDNA was generally not much higher than that found for many genes in nuclear DNA. A comparison of nonsynonymous/synonymous rate ratios in the 13 protein-coding genes suggested differences in the relative levels of selection between haplogroups, including the European haplogroup clusters. Interestingly, a segment of the MTND5 gene was found to be almost void of segregating sites and nonsynonymous mutations in haplogroup J, which has been associated with susceptibility to certain complex diseases. Our results suggest that there are haplogroup-specific differences in the intensity of selection against particular regions of the mitochondrial genome, indicating that some mutations may be non-neutral within specific phylogenetic lineages but neutral within others.     

重庆市26个南亚果实蝇种群mtDNA 16S rRNA基因部分序列及其系统进化

对重庆市26个南亚果实蝇Bactrocera(Zeugodacus)tau(Walker)种群线粒体16S rRNA基因进行测序,获得长约350bp片段的序列。对获得的序列分析表明,A,T,C,G平均含量分别为35·0%,41·3%,7·2%,16·5%,其中保守位点数342个,变异位点数5个,简约信息位点2个,自裔位点2个,所有碱基转换总数为136,替换总数为50。利用MEGA2·1软件重建系统发生树,发现其中21个南亚果实蝇种群未出现分化,另外有5个南亚果实蝇种群出现了分化,但遗传分化程度小。     

Current Intraspecific Dynamics of Sequence Evolution Differs from Long-Term Trends and Can Account for the AT-Richness of Honeybee Mitochondrial DNA

The very high AT content of hymenopteran mtDNA has warranted speculation about nucleotide substitution processes in this group. Here we investigate the pattern of honeybee, Apis mellifera, mtDNA nucleotide polymorphisms inferred from phylogeny in terms of differences between the ATPase6, COI, COII, COIII, cytochrome b, and ND2 genes and strand asymmetry in mutation rates. The observed transition/transversion ratios and the distribution of nonsynonymous substitutions between regions differed significantly. The pattern of differences between genes leading to these heterogeneities (the ATPase6 and COIII genes group apart from the rest) differed markedly from that predicted on the basis of long-term evolutionary change and may indicate differences between current and long-term dynamics of sequence evolution. Also, there is strong strand asymmetry in substitutions, which probably results in a mutability of G and C sufficiently high to account for the AT-richness of honeybee mtDNA. Received: 21 October 1998 / Accepted: 27 January 1999     

Molecular phylogenetic systematics of twelve species of Acipenseriformes based on mtDNA ND4L -ND4 gene sequence analysis

Acipenseriformes is an endangered primitive fish group, which occupies a special place in the history of ideas concerning fish evolution, even in vertebrate evolution. However, the classification and evolution of the fishes have been debated. The mitochondrial DMA (mtDNA) ND4L and partial A7D4 genes were first sequenced in twelve species of the order Acipenseriformes, including endemic Chinese species. The following points were drawn from DNA sequences analysis: (i) the two species of Huso can be ascribed to Acipenser; (ii) A. dabryanus is the mostly closely related to A. sinensis, and most likely the landlocked form of A. sinensis; (iii) genus Acipenser in trans-Pacific region might have a common origin; (iv) mtDNA ND4L and ND4 genes are the ideal genetic markers for phylogenetic analysis of the order Acipenseriformes.     

鲤和鲫线粒体(mtDNA)全基因组分析

为研究鲤(Cyprinus carpio)和鲫(Carassius auratus) mtDNA上的基因分布特征,丰富鱼类mtDNA数据库。本研究通过PCR扩增和测序,获得了鲤和鲫mtDNA基因组序列,经比对分析表明,鲤和鲫的mtDNA序列全长均为16 596 bp,共有13个蛋白质编码基因、22个tRNA基因、2个rRNA基因和1个D-Loop区。鲤和鲫的mtDNA序列中(A+T)百分含量分别为57.2%和57.1%,其碱基组成均具有一定的A/T碱基偏向性。N-J系统发育树分析表明,鲤和鲫与琵琶湖鮈亲缘关系最近,与达氏深水尾魟亲缘关系最远。本研究揭示了鲤和鲫mtDNA遗传变异特征及基因分布规律,为鱼类遗传资源保护及开发利用提供了参考依据。     

Acclimatization of Pink Salmon Oncorhynchus gorbuscha Walbaum in the European North: mtDNA Restriction Data

Pink salmon spawners introduced into the White Sea basin (the Umba River) were compared to the spawners from the basin of the Sea of Okhotsk (the Ola River) using restriction analysis of two fragments of mitochondrial DNA (mtDNA). One of the fragments included genes ND5/ND6, the other, the cytochrome b gene and the D-loop. It was found that mtDNA variation and diversity at the earlier examined nuclear allozyme genes significantly decreased in the odd broodline of pink salmon 8 years after the introduction. The haplotype diversity in the even broodline was considerably lower than in the odd broodline exhibiting virtually no change two generations after the introduction. Based on the results obtained, a possible role of these changes in adaptation of White Sea pink salmon from the odd broodline to the new environment is discussed.     

Somatic mutations of mitochondrial genome in hepatocellular carcinoma

Somatic mutations have been identified in mitochondrial DNA (mtDNA) of various human primary cancers. However, their roles in the pathophysiology of cancers are still unclear. In our previous study, high frequency of somatic mutations was found in the D-loop region of mtDNA of hepatocellular carcinomas (HCCs). In the present study, we examined 44 HCCs and corresponding non-cancerous liver tissues, and identified 13 somatic mutations in the coding region of mtDNAs from 11 HCC samples (11/44, 25%). Among the 13 mtDNA mutations, six mutations (T6787C, G7976A, A9263G, G9267A, A9545G and A11708G) were homoplasmic while seven mutations (956delC, T1659C, G3842A, G5650A, 11032delA, 12418insA and a 66 bp deletion) were heteroplasmic. Moreover, the G3842A transition created a premature stop codon and the 66 bp deletion could omit 22 amino acid residues in the NADH dehydrogenase (ND) subunit 1 (ND1) gene. The 11032delA and 12418insA could result in frame-shift mutation in the ND4 and ND5 genes, respectively. The T1659C transition in tRNA^Val gene and G5650A in tRNA^Ala gene were reported to be clinically associated with some mitochondrial disorders. In addition, the T6787C (cytochrome c oxidase subunit I, COI), G7976A (COII), G9267A (COIII) and A11708G (ND4) mutations could result in amino acid substitutions in the highly conserved regions of the affected mitochondrial genes. These mtDNA mutations (10/13, 76.9%) have the potential to cause mitochondrial dysfunction in HCCs. Taken these results together, we suggest that there may be a higher frequency of mtDNA mutations in HCC than in normal liver tissues from the same individuals.     

Novel features of metazoan mtDNA revealed from sequence analysis of three mitochondrial DNA segments of the land snail Albinaria turrita (Gastropoda: Clausiliidae)

The mitochondrial DNA (mtDNA) size of the terrestrial gastropod Albinaria turrita was determined by restriction enzyme mapping and found to be approximately 14.5 kb. Its partial gene content and organization were examined by sequencing three cloned segments representing about one-fourth of the mtDNA molecule. Complete sequences of cytochrome c oxidase subunit II (COII), and ATPase subunit 8 (ATPase8), as well as partial sequences of cytochrome c oxidase subunit I (COI), NADH dehydrogenase subunit 6 (ND6), and the large ribosomal RNA (IrRNA) genes were determined. Nine putative tRNA genes were also identified by their ability to conform to typical mitochondrial tRNA secondary structures. An 82-nt sequence resembles a noncoding region of the bivalve Mytilus edulis, even though it might contain a tenth tRNA gene with an unusual 5-nt overlap with another tRNA gene. The genetic code of Albinaria turrita appears to be the same as that of Drosophila and Mytilus edulis. The structures of COI and COII are conservative, but those of ATPase8 and ND6 are diversified. The sequenced portion of thelrRNA gene (1,079 nt) is characterized by conspicuous deletions in the 5 and 3 ends; this gene represents the smallest coelomate IrRNA gene so far known. Sequence comparisons of the identified genes indicate that there is greater difference between Albinaria and Mytilus than between Albinaria and Drosophila. An evolutionary analysis, based on COII sequences, suggests a possible nonmonophyletic origin of molluskan mtDNA. This is supported also by the absence of the ATPase8 gene in the mtDNA of Mytilus and nematodes, while this gene is present in the mtDNA of Albinaria and Cepaea nemoralis and in all other known coelomate metazoan mtDNAs.     

Nucleotide sequence of nine protein-coding genes and 22 tRNAs in the mitochondrial DNA of the sea starPisaster ochraceus

Summary We have cloned and sequenced over 9 kb of the mitochondrial genome from the sea starPisaster ochraceus. Within a continuous 8.0-kb fragment are located the genes for NADH dehydrogenase subunits 1, 2, 3, and 4L (ND1, ND2, ND3, and ND4L), cytochrome oxidase subunits I, II, and III (COI, COII, and COIII), and adenosine triphosphatase subunits 6 and 8 (ATPase 6 and ATPase 8). This large fragment also contains a cluster of 13 tRNA genes between ND1 and COI as well as the genes for isoleucine tRNA between ND1 and ND2, arginine tRNA between COI and ND4L, lysine tRNA between COII and ATPase 8, and the serine (UCN) tRNA between COIII and ND3. The genes for the other five tRNAs lie outside this fragment. The gene for phenylalanine tRNA is located between cytochrome b and the 12S ribosomal genes. The genes for tRNA^glu and tRNA^thr are 3 to the 12S ribosomal gene. The tRNAs for histidine and serine (AGN) are adjacent to each other and lie between ND4 and ND5. These data confirm the novel gene order in mitochondrial DNA (mtDNA) of sea stars and delineate additional distinctions between the sea star and other mtDNA molecules.     

Complete nucleotide sequence and gene rearrangement of the mitochondrial genome of the bell-ring frog, Buergeria buergeri (family Rhacophoridae)

In this study we determined the complete nucleotide sequence (19,959 bp) of the mitochondrial DNA of the rhacophorid frog Buergeria buergeri. The gene content, nucleotide composition, and codon usage of B. buergeri conformed to those of typical vertebrate patterns. However, due to an accumulation of lengthy repetitive sequences in the D-loop region, this species possesses the largest mitochondrial genome among all the vertebrates examined so far. Comparison of the gene organizations among amphibian species (Rana, Xenopus, salamanders and caecilians) revealed that the positioning of four tRNA genes and the ND5 gene in the mtDNA of B. buergeri diverged from the common vertebrate gene arrangement shared by Xenopus, salamanders and caecilians. The unique positions of the tRNA genes in B. buergeri are shared by ranid frogs, indicating that the rearrangements of the tRNA genes occurred in a common ancestral lineage of ranids and rhacophorids. On the other hand, the novel position of the ND5 gene seems to have arisen in a lineage leading to rhacophorids (and other closely related taxa) after ranid divergence. Phylogenetic analysis based on nucleotide sequence data of all mitochondrial genes also supported the gene rearrangement pathway.     

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.     

<Emphasis Type="Italic">Gallus gallus</Emphasis> aggrecan gene-based phylogenetic analysis of selected avian taxonomic groups

Mitochondrial DNA (mtDNA) sequences remain the most widely used for phylogenetic analysis in birds. A major limitation of mtDNA sequences, however, is that mitochondria genes are inherited as a single linkage group. Here we describe the use of a 540-bp DNA sequence corresponding to the G3 domain of Gallus gallus nuclear aggrecan gene (AGC1) for phylogenetic analysis of the main groups of Galliformes including Phasianidae, Numididae, and Odontophoridae. We also included species from Cracidae and Megapodiidae which are considered by some as Craciformes and others, including here as Galliformes. The uncorrected sequence divergence of the G3 fragments ranges from 1 among the grouses to 36% between some of the distant groups within Galliformes. These sequences contain 39–48% AT nucleotides and the ratios of transition versus transversion are above 1.5 in majority of the comparisons. Using G3 sequences from an Anseriform, Oxyura jamaicensis, as out-groups, phylogenetic trees were obtained using maximum parsimony and distance algorithms and bootstrap analyses. These trees were consistent with those described using Avian sarcoma and leucosis virus gag genes and those from amino acid sequences of hemoglobin and lysozyme c. Our data also support relationships among Galliformes which were defined using mtDNA sequences. In addition to the general support of the five main families of Galliformes, our data are also consistent with previous work that showed Francolinus africanus and Gallus gallus are in the same clade and that Tetraoninae is a well-supported monophyletic subfamily within Phasianidae. The results presented here suggest that the AGC1 sequences meet the criterion of novel nuclear DNA sequences that can be used to help resolve the relationships among Galliformes.     

Neutral and Non-Neutral Evolution of Drosophila Mitochondrial DNA

To test hypotheses of neutral evolution of mitochondrial DNA (mtDNA), nucleotide sequences were determined for 1515 base pairs of the NADH dehydrogenase subunit 5 (ND5) gene in the mitochondrial DNA of 29 lines of Drosophila melanogaster and 9 lines of its sibling species Drosophila simulans. In contrast to the patterns for nuclear genes, where D. melanogaster generally exhibits much less nucleotide polymorphism, the number of segregating sites was slightly higher in a global sample of nine ND5 sequences in D. melanogaster (s = 8) than in the nine lines of D. simulans (s = 6). When compared to variation at nuclear loci, the mtDNA variation in D. melanogaster does not depart from neutral expectations. The ND5 sequences in D. simulans, however, show fewer than half the number of variable sites expected under neutrality when compared to sequences from the period locus. While this reduction in variation is not significant at the 5% level, HKA tests with published restriction data for mtDNA in D. simulans do show a significant reduction of variation suggesting a selective sweep of variation in the mtDNA in this species. Tests of neutral evolution based on the ratios of synonymous and replacement polymorphism and divergence are generally consistent with neutral expectations, although a significant excess of amino acid polymorphism within both species is localized in one region of the protein. The rate of mtDNA evolution has been faster in D. melanogaster than in D. simulans and the population structure of mtDNA is distinct in these species. The data reveal how different rates of mtDNA evolution between species and different histories of neutral and adaptive evolution within species can compromise historical inferences in population and evolutionary biology.