Slow mitochondrial DNA sequence evolution in the Anthozoa (Cnidaria)

Mitochondrial genes have been used extensively in population genetic and phylogeographical analyses, in part due to a high rate of nucleotide substitution in animal mitochondrial DNA (mtDNA). Nucleotide sequences of anthozoan mitochondrial genes, however, are virtually invariant among conspecifics, even at third codon positions of protein-coding sequences. Hence, mtDNA markers are of limited use for population-level studies in these organisms. Mitochondrial gene sequence divergence among anthozoan species is also low relative to that exhibited in other animals, although higher level relationships can be resolved with these markers. Substitution rates in anthozoan nuclear genes are much higher than in mitochondrial genes, whereas nuclear genes in other metazoans usually evolve more slowly than, or similar to, mitochondrial genes. Although several mechanisms accounting for a slow rate of sequence evolution have been proposed, there is not yet a definitive explanation for this observation. Slow evolution and unique characteristics may be common in primitive metazoans, suggesting that patterns of mtDNA evolution in these organisms differ from that in other animal systems.     

Animal mitochondrial DNA as a genetic marker in population and evolutionary biology

Animal mitochondrial DNA (mtDNA) is playing an increasingly important role as a genetic marker in population and evolutionary biology. The popularity of this molecule derives, in part, from the relative ease with which clearly homologous sequences can be isolated and compared. Simple sequence organization, maternal inheritance and absence of recombination make mtDNA an ideal marker for tracing maternal genealogies. Rapid rate of sequence divergence (at least in vertebrates) allows discrimination of recently diverged lineages. Studies of mtDNAs from a diversity of animal groups have revealed significant variation among taxa in mtDNA sequence dynamics, gene order and genome size. They have also provided important insights into population structure, geographic variation, zoogeography and phylogeny.     

Mitochondrial pseudogenes: evolution's misplaced witnesses

Nuclear copies of mitochondrial DNA (mtDNA) have contaminated PCR-based mitochondrial studies of over 64 different animal species. Since the last review of these nuclear mitochondrial pseudogenes (Numts) in animals, Numts have been found in 53 of the species studied. The recent evidence suggests that Numts are not equally abundant in all species, for example they are more common in plants than in animals, and also more numerous in humans than in Drosophila. Methods for avoiding Numts have now been tested, and several recent studies demonstrate the potential utility of Numt DNA sequences in evolutionary studies. As relics of ancient mtDNA, these pseudogenes can be used to infer ancestral states or root mitochondrial phylogenies. Where they are numerous and selectively unconstrained, Numts are ideal for the study of spontaneous mutation in nuclear genomes.     

A novel mitochondrial DNA-like sequence in the human nuclear genome.

We describe here a nuclear mitochondrial DNA-like sequence (numtDNA) that is nearly identical in sequence to a continuous 5842 bp segment of human mitochondrial DNA (mtDNA) that spans nucleotide positions 3914 to 9755. On the basis of evolutionary divergence among modern primates, this numtDNA molecule appears to represent mtDNA from a hominid ancestor that has been translocated to the nuclear genome during the recent evolution of humans. This numtDNA sequence harbors synonymous and nonsynonymous nucleotide substitutions relative to the authentic human mtDNA sequence, including an array of substitutions that was previously found in the cytochrome c oxidase subunit 1 and 2 genes. These substitutions were previously reported to occur in human mtDNA, but subsequently contended to be present in a nuclear pseudogene sequence. We now demonstrate their exclusive association with this 5842-bp numtDNA, which we have characterized in its entirety. This numtDNA does not appear to be expressed as a mtDNA-encoded mRNA. It is present in nuclear DNA from human blood donors, in human SH-SY5Y and A431 cell lines, and in rho(0) SH-SY5Y and rho(0) A431 cell lines that were depleted of mtDNA. The existence of human numtDNA sequences with great similarities to human mtDNA renders the amplification of pure mtDNA from cellular DNA very difficult, thereby creating the potential for confounding studies of mitochondrial diseases and population genetics.     

Mitochondrial genome sequence evolution in Chlamydomonas

The mitochondrial genomes of the Chlorophyta exhibit significant diversity with respect to gene content and genome compactness; however, quantitative data on the rates of nucleotide substitution in mitochondrial DNA, which might help explain the origin of this diversity, are lacking. To gain insight into the evolutionary forces responsible for mitochondrial genome diversification, we sequenced to near completion the mitochondrial genome of the chlorophyte Chlamydomonas incerta, estimated the evolutionary divergence between Chlamydomonas reinhardtii and C. incerta mitochondrial protein-coding genes and rRNA-coding regions, and compared the relative evolutionary rates in mitochondrial and nuclear genes. Synonymous and nonsynonymous substitution rates do not differ significantly between the mitochondrial and nuclear protein-coding genes. The mitochondrial rRNA-coding regions, however, are evolving much faster than their nuclear counterparts, and this difference might be explained by relaxed functional constraints on the mitochondrial translational apparatus due to the small number of proteins synthesized in Chlamydomonas mitochondria. Substitution rates at synonymous sites in a nonstandard mitochondrial gene (rtl) and at intronic and synonymous sites in nuclear genes expressed at low levels suggest that the mutation rate is similar in these two genetic compartments. Potential evolutionary forces shaping mitochondrial genome evolution in Chlamydomonas are discussed.     

Mitochondrial DNA sequences of primates: Tempo and mode of evolution

Summary We cloned and sequenced a segment of mitochondrial DNA from human, chimpanzee, gorilla, orangutan, and gibbon. This segment is 896 bp in length, contains the genes for three transfer RNAs and parts of two proteins, and is homologous in all 5 primates. The 5 sequences differ from one another by base substitutions at 283 positions and by a deletion of one base pair. The sequence differences range from 9 to 19% among species, in agreement with estimates from cleavage map comparisons, thus confirming that the rate of mtDNA evolution in primates is 5 to 10 times higher than in nuclear DNA. The most striking new finding to emerge from these comparisons is that transitions greatly outnumber transversions. Ninety-two percent of the differences among the most closely related species (human, chimpanzee, and gorilla) are transitions. For pairs of species with longer divergence times, the observed percentage of transitions falls until, in the case of comparisons between primates and non-primates, it reaches a value of 45. The time dependence is probably due to obliteration of the record of transitions by multiple substitutions at the same nucleotide site. This finding illustrates the importance of choosing closely related species for analysis of the evolutionary process. The remarkable bias toward transitions in mtDNA evolution necessitates the revision of equations that correct for multiple substitutions at the same site. With revised equations, we calculated the incidence of silent and replacement substitutions in the two protein-coding genes. The silent substitution rate is 4 to 6 times higher than the replacement rate, indicating strong functional constraints at replacement sites. Moreover, the silent rate for these two genes is about 10% per million years, a value 10 times higher than the silent rate for the nuclear genes studied so far. In addition, the mean substitution rate in the three mitochondrial tRNA genes is at least 100 times higher than in nuclear tRNA genes. Finally, genealogical analysis of the sequence differences supports the view that the human lineage branched off only slightly before the gorilla and chimpanzee lineages diverged and strengthens the hypothesis that humans are more related to gorillas and chimpanzees than is the orangutan.Abbreviations mtDNA mitochondrial DNA - bp base pair - URF unidentified reading frame     

On shortcomings of using mtDNA sequence divergence for the systematics of hares (genus Lepus): An example from cape hares

Despite the well-known fact that evolutionary patterns of single genes or sequences are not necessarily paralleled by organismic evolution, using mitochondrial sequence divergence for inferring phylogenetic relationships among hare species is not uncommon. Earlier studies reported interspecific introgression in the genus Lepus that may slur systematic conclusions drawn exclusively from mtDNA data. We examined pairwise divergence in partial mtHV-1 sequences and nuclear gene pools based on microsatellite and allozyme loci separately in a South and North African population of cape hares, Lepus capensis, and did not find significant correspondence on the individual level within either population. Also, the former population had a significantly higher level of mtHV-1 sequence divergence compared to the latter, but levels of individual nuclear gene pool divergence did not differ significantly between the two populations. Hence, the absence of correspondence between mtHV-1 sequence divergence and nuclear gene pool divergence among individuals within a population might be independent of the population-specific level of differentiation among individuals. Our results complement earlier findings in hares, and strongly recommend to include nuclear gene pool evidence for systematic inferences within this genus, even under the absence of mitochondrial introgression.     

Complete Mitochondrial Genomes and Eutherian Evolution

Recent large-scale nuclear DNA phylogenies have supported unconventional interordinal relationships among modern eutherians as well as divergence dates (100 mya) that substantially predate the first appearance of fossils from modern eutherians near the Cretaceous/Cenozoic (K/T) boundary (65-70 mya). For comparison to the nuclear data, I analyzed 12 complete mitochondrial DNA (mtDNA) protein-coding genes (10,677 bp) from 53 eutherian taxa, using maximum-likelihood methods to estimate model parameters (GTR + I + ) and to optimize topology and branch-length estimates. Although closely resembling the nuclear DNA trees, the mtDNA maximum-likelihood tree is just one of seven statistically indistinguishable ( lnL 1.747) trees, each suggesting different evolutionary relationships. This 53-taxon data set and another including 56 taxa provide no statistically significant support for a monophyletic afrotherian clade. In fact, these mitochondrial DNA sequences fail to support the monophyly of three putative eutherian divisions suggested by the nuclear data (Afrotheria, Laurasiatheria or Euarchontoglires). By comparison to well-supported branches describing relationships among families, those describing interordinal relationships are extremely short and only tenuously supported. Neither these sequences, nor sequences simulated under a known tree, fully resolve any interordinal relationship. Even simulated sequences that are twice as long (22kb) as mtDNA protein-coding genes are too short and too saturated to resolve the deepest and shortest interordinal relationships. Further, the mammalian mtDNA sequences appear to depart significantly from molecular-clock and quartet dating assumptions. Unlike recent nuclear DNA studies, I find that mtDNA genes, by themselves, are inadequate to describe relationships or divergence times at the base of the eutherian tree.     

A Method for Calibrating Molecular Clocks and Its Application to Animal Mitochondrial DNA

A generalized least-squares procedure is introduced for the calibration of molecular clocks and applied to the complete mitochondrial DNA sequences of 13 animal species. The proposed technique accounts for both nonindependence and heteroscedasticity of molecular-distance data, problems that have not been taken into to account in such analyses in the past. When sequence-identity data are transformed to account for multiple substitutions/site, the molecular divergence scales linearly with time, but with substantially more variation in the substitution rate than expected under a Poisson model. Significant levels of divergence are predicted at zero divergence time for most loci, suggesting high levels of site-specific heterozygosity among mtDNA molecules establishing in sister taxa. For nearly all loci, the baseline heterozygosity is lower and the substitution rate is higher in mammals relative to other animals. There is considerable variation in the evolutionary rate among loci but no compelling evidence that the average rate of mtDNA evolution is elevated with respect to that of nuclear DNA. Using the observed patterns of interspecific divergence, empirical estimates are derived for the mean coalescence times of organelles colonizing sister taxa.     

The microevolution of the Galápagos marine iguana Amblyrhynchus cristatus assessed by nuclear and mitochondrial genetic analyses

Marine iguanas may have inhabited the Galápagos archipelago and its former, now sunken islands for more than 10 million years (Myr). It is therefore surprising that morphological and immunological data indicate little evolutionary divergence within the genus. We utilized mitochondrial DNA (mtDNA) sequence analyses and nuclear DNA fingerprinting to re-evaluate the level and pattern of genetic differentiation among 22 marine iguana populations from throughout the archipelago. Both genetic marker systems detect a low level of within-genus divergence, but they show contrasting levels of geographical subdivision among the populations. The mitochondrial gene pools of populations from different regions of the archipelago are isolated, and the mtDNA pattern appears to follow the sequence in which the islands were colonized by marine iguanas. Conversely, the nuclear DNA study indicates substantial interpopulational gene exchange, and the geographical distribution of the nuclear markers seems to be determined by isolation by distance among the populations. The natural history of marine iguanas suggests that the contrasting nuclear and mitochondrial DNA patterns result from an asymmetric migration behaviour of the two sexes, with higher (active and passive) interisland dispersal for males than females. Separate genetic analyses for the sexes appear to support this hypophesis. Based on these findings, a scenario is proposed that explains the marine iguanas' low genetic divergence, notwithstanding their long evolutionary history in the Galápagos archipelago.     

Conflicting patterns of mitochondrial and nuclear DNA diversity in Phylloscopus warblers

Molecular variation is often used to infer the demographic history of species, but sometimes the complexity of species history can make such inference difficult. The willow warbler, Phylloscopus trochilus, shows substantially less geographical variation than the chiffchaff, Phylloscopus collybita, both in morphology and in mitochondrial DNA (mtDNA) divergence. We therefore predicted that the willow warbler should harbour less nuclear DNA diversity than the chiffchaff. We analysed sequence data obtained from multiple samples of willow warblers and chiffchaffs for the mtDNA cytochrome b gene and four nuclear genes. We confirmed that the mtDNA diversity among willow warblers is low (pi = 0.0021). Sequence data from three nuclear genes (CHD-Z, AFLP-WW1 and MC1R) not linked to the mitochondria demonstrated unexpectedly high nucleotide diversity (pi values of 0.0172, 0.0141 and 0.0038) in the willow warbler, on average higher than the nucleotide diversity for the chiffchaff (pi values of 0.0025, 0.0017 and 0.0139). In willow warblers, Tajima's D analyses showed that the mtDNA diversity, but not the nuclear DNA diversity, has been reduced relative to the neutral expectation of molecular evolution, suggesting the action of a selective sweep affecting the maternally inherited genes. The large nuclear diversity seen within willow warblers is not compatible with processes of neutral evolution occurring in a population with a constant population size, unless the long-term effective population size has been very large (N(e) > 10(6)). We suggest that the contrasting patterns of genetic diversity in the willow warbler may reflect a more complex evolutionary history, possibly including historical demographic fluctuations or historical male-biased introgression of nuclear genes from a differentiated population of Phylloscopus warblers.     

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.     

Assessing Nuclear and Mitochondrial DNA Sequence Variation Within Steinernema (Rhabditida: Steinernematidae)

DNA sequence analysis was used to characterize the nuclear ribosomal DNA ITS1 region and a portion of the COII and 16S rDNA genes of the mitochondrial genome from Steinernema entomopathogenic nematodes. Nuclear ITS1 nucleotide divergence among seven Steinernema spp. ranged from 6 to 22%, and mtDNA divergence among five species ranged from 12 to 20%. No intraspecific variation was observed among three S. feltiae strains. Phylogenetic analysis of both nuclear and mitochondrial DNA sequences confirms the existing morphological relationships of several Steinernema species. Both the rDNA ITS1 and mtDNA sequences were useful for resolving relationships among Steinernema taxa.     

Frequent assimilation of mitochondrial DNA by grasshopper nuclear genomes

Multiple copies of mitochondrial-like DNA were found in the brown mountain grasshopper, Podisma pedestris (Orthoptera: Acrididae), paralogous to COI and ND5 regions. The same was discovered using the ND5 regions of nine other grasshopper species from four separate subfamilies (Podisminae, Calliptaminae, Cyrtacanthacridinae, and Gomphocerinae). The extra ND5-like sequences were shown to be nuclear in the desert locust, Schistocerca gregaria (Cyrtacanthacridinae), and probably so in P. pedestris and an Italopodisma sp. (Podisminae). Eighty-seven different ND5-like nuclear mitochondrial pseudogenes (Numts) were sequenced from 12 grasshopper individuals. Different nuclear mitochondrial pseudogenes, if descended from the same mitochondrial immigrant, will have diverged from each other under no selective constraints because of their loss of functionality. Evidence of selective constraints in the differences between any two Numt sequences (e.g., if most differences are at third positions of codons) implies that they have separate mitochondrial origins. Through pairwise comparisons of pseudogene sequences, it was established that there have been at least 12 separate mtDNA integrations into P. pedestris nuclear genomes. This is the highest reported rate of horizontal transfer between organellar and nuclear genomes within a single animal species. The occurrence of numerous mitochondrial pseudogenes in nuclear genomes derived from separate integration events appears to be a common phenomenon among grasshoppers. More than one type of mechanism appears to have been involved in generating the observed grasshopper Numts.     

The problems of molecular phylogenetics with the example of squamate reptiles: Mitochondrial DNA markers

The review considers the current problems of molecular phylogenetics based on mitochondrial and chromosomal DNA sequences. The emphasis is placed on mtDNA markers, which are widely employed in reconstructing molecular evolution, but often without a critical analysis of the physiological and biochemical features of mitochondria that affect the adequacy and reliability of the results. In addition to the factors that make mtDNA-based phylogenies difficult to interpret (unrecognized hybridization and introgression events, ancestral polymorphism, and nuclear paralogs of mtDNA sequences), attention is paid to the nonneutrality and unequal mutation rates of mtDNA genes and their fragments, violations of uniparental inheritance of mitochondria, recombination events, natural heteroplasmy, and mtDNA haplotypic diversity. These factors may influence the congruence of phylogenetic inferences and trees constructed for the same organisms with different mtDNA markers or with mitochondrial and nuclear markers. The review supports the viewpoint that mitochondrial genes and their fragments fail to provide reliable evolutionary markers when considered without a thorough study of the environmental conditions and life of the taxa. The influence of external conditions on the metabolism and physiology of mitochondria cannot be taken into account in full nor modeled well enough for phylogenetic applications. It is assumed that mtDNA is valuable as a phylogenetic marker primarily because its complete sequence may be analyzed to identify the apomorphic and synmorphic properties of a taxon and to search for informative nuclear paralogs of mtDNA for phylogeographical studies and estimations of relative evolution times.     

Free from mitochondrial DNA: Nuclear genes and the inference of species trees among closely related darter lineages (Teleostei: Percidae: Etheostomatinae)

Investigations into the phylogenetics of closely related animal species are dominated by the use of mitochondrial DNA (mtDNA) sequence data. However, the near-ubiquitous use of mtDNA to infer phylogeny among closely related animal lineages is tempered by an increasing number of studies that document high rates of transfer of mtDNA genomes among closely related species through hybridization, leading to substantial discordance between phylogenies inferred from mtDNA and nuclear gene sequences. In addition, the recent development of methods that simultaneously infer a species phylogeny and estimate divergence times, while accounting for incongruence among individual gene trees, has ushered in a new era in the investigation of phylogeny among closely related species. In this study we assess if DNA sequence data sampled from a modest number of nuclear genes can resolve relationships of a species-rich clade of North American freshwater teleost fishes, the darters. We articulate and expand on a recently introduced method to infer a time-calibrated multi-species coalescent phylogeny using the computer program ^*BEAST. Our analyses result in well-resolved and strongly supported time-calibrated darter species tree. Contrary to the expectation that mtDNA will provide greater phylogenetic resolution than nuclear gene data; the darter species tree inferred exclusively from nuclear genes exhibits a higher frequency of strongly supported nodes than the mtDNA time-calibrated gene tree.     

Mitochondrial DNA in the sea urchin Arbacia lixula: evolutionary inferences from nucleotide sequence analysis.

From the stirodont Arbacia lixula we determined the sequence of 5,127 nucleotides of mitochondrial DNA (mtDNA) encompassing 18 tRNAs, two complete coding genes, parts of three other coding genes, and part of the 12S ribosomal RNA (rRNA). The sequence confirms that the organization of mtDNA is conserved within echinoids. Furthermore, it underlines the following peculiar features of sea urchin mtDNA: the clustering of tRNAs, the short noncoding regulatory sequence, and the separation by the ND1 and ND2 genes of the two rRNA genes. Comparison with the orthologous sequences from the camarodont species Paracentrotus lividus and Strongylocentrotus purpuratus revealed that (1) echinoids have an extra piece on the amino terminus of the ND5 gene that is probably the remnant of an old leucine tRNA gene; (2) third-position codon nucleotide usage has diverged between A. lixula and the camarodont species to a significant extent, implying different directional mutational pressures; and (3) the stirodont-camarodont divergence occurred twice as long ago as did the P. lividus-S. purpuratus divergence.     

Structure and chromosomal distribution of human mitochondrial pseudogenes

Nuclear mitochondrial pseudogenes (Numts) have been found in the genome of many eukaryote species, including humans. Using a BLAST approach, we found 1105 DNA sequences homologous to mitochondrial DNA (mtDNA) in the August 2001 Goldenpath human genome database. We assembled these sequences manually into 286 pseudogenes on the basis of single insertion events and constructed a chromosomal map of these Numts. Some pseudogenes appeared highly modified, containing inversions, deletions, duplications, and displaced sequences. In the case of four randomly selected Numts, we used PCR tests on cells lacking mtDNA to ensure that our technique was free from genome-sequencing artifacts. Furthermore, phylogenetic investigation suggested that one Numt, apparently inserted into the nuclear genome 25-30 million years ago, had been duplicated at least 10 times in various chromosomes during the course of evolution. Thus, these pseudogenes should be very useful in the study of ancient mtDNA and nuclear genome evolution.     

Seventeen new complete mtDNA sequences reveal extensive mitochondrial genome evolution within the Demospongiae

Two major transitions in animal evolution--the origins of multicellularity and bilaterality--correlate with major changes in mitochondrial DNA (mtDNA) organization. Demosponges, the largest class in the phylum Porifera, underwent only the first of these transitions and their mitochondrial genomes display a peculiar combination of ancestral and animal-specific features. To get an insight into the evolution of mitochondrial genomes within the Demospongiae, we determined 17 new mtDNA sequences from this group and analyzing them with five previously published sequences. Our analysis revealed that all demosponge mtDNAs are 16- to 25-kbp circular molecules, containing 13-15 protein genes, 2 rRNA genes, and 2-27 tRNA genes. All but four pairs of sampled genomes had unique gene orders, with the number of shared gene boundaries ranging from 1 to 41. Although most demosponge species displayed low rates of mitochondrial sequence evolution, a significant acceleration in evolutionary rates occurred in the G1 group (orders Dendroceratida, Dictyoceratida, and Verticillitida). Large variation in mtDNA organization was also observed within the G0 group (order Homosclerophorida) including gene rearrangements, loss of tRNA genes, and the presence of two introns in Plakortis angulospiculatus. While introns are rare in modern-day demosponge mtDNA, we inferred that at least one intron was present in cox1 of the common ancestor of all demosponges. Our study uncovered an extensive mitochondrial genomic diversity within the Demospongiae. Although all sampled mitochondrial genomes retained some ancestral features, including a minimally modified genetic code, conserved structures of tRNA genes, and presence of multiple non-coding regions, they vary considerably in their size, gene content, gene order, and the rates of sequence evolution. Some of the changes in demosponge mtDNA, such as the loss of tRNA genes and the appearance of hairpin-containing repetitive elements, occurred in parallel in several lineages and suggest general trends in demosponge mtDNA evolution.     

Speciation, introgressive hybridization and nonlinear rate of molecular evolution in flycatchers

Evolutionary history of Muscicapidae flycatchers is inferred from nuclear and mitochondrial DNA (mtDNA) sequence comparisons and population genetic analysis of nuclear and mtDNA markers. Phylogenetic reconstruction based on sequences from the two genomes yielded similar trees with respect to the order at which the species split off. However, the genetic distances fitted a nonlinear, polynomial model reflecting diminishing divergence rate of the mtDNA sequences compared to the nuclear DNA sequences. This could be explained by Haldane's rule because genetic isolation might evolve more rapidly on the mitochondrial rather than the nuclear genome in birds. This is because hybrid sterility of the heterogametic sex (females) would predate that of the homogametic sex (males), leading to sex biased introgression of nuclear genes. Analyses of present hybrid zones of pied (Ficedula hypoleuca) and collared flycatchers (F. albicollis) may indicate a slight sexual bias in rate of introgression, but the introgression rates were too low to allow proper statistical analyses. It is suggested, however, that the observed deviation from linearity can be explained by a more rapid mutational saturation of the mtDNA sequences than of the nuclear DNA sequences, as supported by analyses of third codon position transversions at two protein coding mtDNA genes. A phylogeographic scenario for the black and white flycatcher species is suggested based on interpretation of the genetic data obtained. Four species appear to have diverged from a common ancestor relatively simultaneously during the Pleistocene. After the last glaciation period, pied and collared flycatchers expanded their breeding ranges and eventually came into secondary contact in Central and Eastern Europe and on the Baltic Isles.