A comprehensive evolutionary classification of proteins encoded in complete eukaryotic genomes

Background

Sequencing the genomes of multiple, taxonomically diverse eukaryotes enables in-depth comparative-genomic analysis which is expected to help in reconstructing ancestral eukaryotic genomes and major events in eukaryotic evolution and in making functional predictions for currently uncharacterized conserved genes.

Results

We examined functional and evolutionary patterns in the recently constructed set of 5,873 clusters of predicted orthologs (eukaryotic orthologous groups or KOGs) from seven eukaryotic genomes: Caenorhabditis elegans, Drosophila melanogaster, Homo sapiens, Arabidopsis thaliana, Saccharomyces cerevisiae, Schizosaccharomyces pombe and Encephalitozoon cuniculi. Conservation of KOGs through the phyletic range of eukaryotes strongly correlates with their functions and with the effect of gene knockout on the organism's viability. The approximately 40% of KOGs that are represented in six or seven species are enriched in proteins responsible for housekeeping functions, particularly translation and RNA processing. These conserved KOGs are often essential for survival and might approximate the minimal set of essential eukaryotic genes. The 131 single-member, pan-eukaryotic KOGs we identified were examined in detail. For around 20 that remained uncharacterized, functions were predicted by in-depth sequence analysis and examination of genomic context. Nearly all these proteins are subunits of known or predicted multiprotein complexes, in agreement with the balance hypothesis of evolution of gene copy number. Other KOGs show a variety of phyletic patterns, which points to major contributions of lineage-specific gene loss and the 'invention' of genes new to eukaryotic evolution. Examination of the sets of KOGs lost in individual lineages reveals co-elimination of functionally connected genes. Parsimonious scenarios of eukaryotic genome evolution and gene sets for ancestral eukaryotic forms were reconstructed. The gene set of the last common ancestor of the crown group consists of 3,413 KOGs and largely includes proteins involved in genome replication and expression, and central metabolism. Only 44% of the KOGs, mostly from the reconstructed gene set of the last common ancestor of the crown group, have detectable homologs in prokaryotes; the remainder apparently evolved via duplication with divergence and invention of new genes.

Conclusions

The KOG analysis reveals a conserved core of largely essential eukaryotic genes as well as major diversification and innovation associated with evolution of eukaryotic genomes. The results provide quantitative support for major trends of eukaryotic evolution noticed previously at the qualitative level and a basis for detailed reconstruction of evolution of eukaryotic genomes and biology of ancestral forms.     

Repeated regions in mitochondrial genomes: Distribution, origin and evolutionary significance

All complete or nearly complete mitochondrial genomes of Metazoa (2819) have been subject to bioinformatic analysis to investigate the distribution and features of repeated and palindromic sequences. Repeats are ubiquitous, with 29.9% of genomes containing at least one and 1.95% of total genome length being repeated. Repeat boundaries were tested for the presence of secondary structure motifs, consensus sequences or small repeats, features generally reported as associated with duplications. No significant relationship was detected, suggesting the non ubiquitousness of such features. A mechanism related to gene conversion is proposed to explain the origin of small interspersed repeats.     

Structural variation of mitochondrial genomes sheds light on evolutionary history of soybeans

The architecture and genetic diversity of mitogenome (mtDNA) are largely unknown in cultivated soybean (Glycine max), which is domesticated from the wild progenitor, Glycine soja, 5000 years ago. Here, we de novo assembled the mitogenome of the cultivar ‘Williams 82’ (Wm82_mtDNA) with Illumina PE300 deep sequencing data, and verified it with polymerase chain reaction (PCR) and Southern blot analyses. Wm82_mtDNA maps as two autonomous circular chromosomes (370 871-bp Chr-m1 and 62 661-bp Chr-m2). Its structure is extensively divergent from that of the mono-chromosomal mitogenome reported in the landrace ‘Aiganhuang’ (AGH_mtDNA). Synteny analysis showed that the structural variations (SVs) between two genomes are mainly attributed to ectopic and illegitimate recombination. Moreover, Wm82_mtDNA and AGH_mtDNA each possess six and four specific regions, which are absent in their counterparts and likely result from differential sequence-loss events. Mitogenome SV was further studied in 39 wild and 182 cultivated soybean accessions distributed world-widely with PCR/Southern analyses or a comparable in silico analysis. The results classified both wild and cultivated soybeans into five cytoplasmic groups, named as GSa–GSe and G1–G5; ‘Williams 82’ and ‘Aiganhuang’ belong to G1 and G5, respectively. Notably, except for members in GSe and G5, all accessions carry a bi-chromosomal mitogenome with a common Chr-m2. Phylogenetic analyses based on mtDNA structures and chloroplast gene sequences both inferred that G1–G3, representing >90% of cultigens, likely inherited cytoplasm from the ancestor of domestic soybean, while G4 and G5 likely inherited cytoplasm from wild soybeans carrying GSa- and GSe-like cytoplasm through interspecific hybridization, offering new insights into soybean cultivation history.     

A comprehensive analysis of Trypanosoma brucei mitochondrial proteome

The composition of the large, single, mitochondrion (mt) of Trypanosoma brucei was characterized by MS (2‐D LC‐MS/MS and gel‐LC‐MS/MS) analyses. A total of 2897 proteins representing a substantial proportion of procyclic form cellular proteome were identified, which confirmed the validity of the vast majority of gene predictions. The data also showed that the genes annotated as hypothetical (species specific) were overpredicted and that virtually all genes annotated as hypothetical, unlikely are not expressed. By comparing the MS data with genome sequence, 40 genes were identified that were not previously predicted. The data are placed in a publicly available web‐based database (www.TrypsProteome.org). The total mitochondrial proteome is estimated at 1008 proteins, with 401, 196, and 283 assigned to the mt with high, moderate, and lower confidence, respectively. The remaining mitochondrial proteins were estimated by statistical methods although individual assignments could not be made. The identified proteins have predicted roles in macromolecular, metabolic, energy generating, and transport processes providing a comprehensive profile of the protein content and function of the T. brucei mt.     

Comparative analysis of mitochondrial genomes of Nirvanini and Evacanthini (Hemiptera: Cicadellidae) reveals an explicit evolutionary relationship

Mitogenomes of five leafhopper species, Chudania hellerina and Concaveplana rufolineata in Nirvanini, Carinata rufipenna, Evacanthus danmainus and E. heimianus representing Evacanthini, were sequenced. The lengths of these five mitogenomes range from 15,044 (C. hellerina) to 15,680 bp (E. heimianus). All five mitogenomes exhibit similar base composition, gene size and codon usage of protein-coding genes. All 22 tRNA genes have typical cloverleaf secondary structures, except for trnS1 (AGN) which appears to lack the dihydrouridine arm. The two included Nirvanini species employ the anticodon TCT instead of the commonly used GCT in trnS1 (AGN). Genes nad2, atp8 and nad6 were highly variable while cox1 and cob showed the lowest nucleotide diversity. Phylogenetic analyses of two concatenated nucleotide datasets, incorporating the newly sequenced taxa and other available membracoid mitogenomes, recovered each included leafhopper subfamily as monophyletic with evacanthine tribes Nirvanini and Evacanthini forming monophyletic sister clades. A relationship among Evacanthinae, Cicadellinae and Typhlocybinae received moderate branch support.     

Towards a comprehensive picture of alloacceptor tRNA remolding in metazoan mitochondrial genomes

Remolding of tRNAs is a well-documented process in mitochondrial genomes that changes the identity of a tRNA. It involves a duplication of a tRNA gene, a mutation that changes the anticodon and the loss of the ancestral tRNA gene. The net effect is a functional tRNA that is more closely related to tRNAs of a different alloacceptor family than to tRNAs with the same anticodon in related species. Beyond being of interest for understanding mitochondrial tRNA function and evolution, tRNA remolding events can lead to artifacts in the annotation of mitogenomes and thus in studies of mitogenomic evolution. Therefore, it is important to identify and catalog these events. Here we describe novel methods to detect tRNA remolding in large-scale data sets and apply them to survey tRNA remolding throughout animal evolution. We identify several novel remolding events in addition to the ones previously mentioned in the literature. A detailed analysis of these remoldings showed that many of them are derived from ancestral events.     

Metazoan OXPHOS gene families: evolutionary forces at the level of mitochondrial and nuclear genomes

Mitochondrial and nuclear DNAs contribute to encode the whole mitochondrial protein complement. The two genomes possess highly divergent features and properties, but the forces influencing their evolution, even if different, require strong coordination. The gene content of mitochondrial genome in all Metazoa is in a frozen state with only few exceptions and thus mitochondrial genome plasticity especially concerns some molecular features, i.e. base composition, codon usage, evolutionary rates. In contrast the high plasticity of nuclear genomes is particularly evident at the macroscopic level, since its redundancy represents the main feature able to introduce genetic material for evolutionary innovations. In this context, genes involved in oxidative phosphorylation (OXPHOS) represent a classical example of the different evolutionary behaviour of mitochondrial and nuclear genomes. The simple DNA sequence of Cytochrome c oxidase I (encoded by the mitochondrial genome) seems to be able to distinguish intra- and inter-species relations between organisms (DNA Barcode). Some OXPHOS subunits (cytochrome c, subunit c of ATP synthase and MLRQ) are encoded by several nuclear duplicated genes which still represent the trace of an ancient segmental/genome duplication event at the origin of vertebrates.     

Analysis of the evolutionary forces shaping mitochondrial genomes of a Neotropical malaria vector complex

Many vectors of human malaria belong to complexes of morphologically indistinguishable cryptic species. Here we report the analysis of the newly sequenced complete mitochondrial DNA molecules from six recognized or putative species of one such group, the Neotropical Anopheles albitarsis complex. The molecular evolution of these genomes had been driven by purifying selection, particularly strongly acting on the RNA genes. Directional mutation pressure associated with the strand-asynchronous asymmetric mtDNA replication mechanism may have shaped a pronounced DNA strand asymmetry in the nucleotide composition in these and other Anopheles species. The distribution of sequence polymorphism, coupled with the conflicting phylogenetic trees inferred from the mitochondrial DNA and from the published white gene fragment sequences, indicates that the evolution of the complex may have involved ancient mtDNA introgression. Six protein coding genes (nad5, nad4, cox3, atp6, cox1 and nad2) have high levels of sequence divergence and are likely informative for population genetics studies. Finally, the extent of the mitochondrial DNA variation within the complex supports the notion that the complex consists of a larger number of species than until recently believed.     

The complete mitochondrial genomes of Schistosoma haematobium and Schistosoma spindale and the evolutionary history of mitochondrial genome changes among parasitic flatworms

Complete mitochondrial genome sequences for the schistosomes Schistosoma haematobium and Schistosoma. spindale have been characterized. S. haematobium is the causative agent of urinary schistosomiasis in humans and S. spindale uses ruminants as its definitive host; both are transmitted by freshwater snail intermediate hosts. Results confirm a major gene order rearrangement among schistosomes in all traditional Schistosoma species groups other than Schistosoma japonicum; i.e., species groups S. mansoni, S. haematobium, and S. indicum. These data lend support to the 'out of Asia' (East and Southeast Asia) hypothesis for Schistosoma. The gene order change involves translocation of atp6-nad2-trnA and a rearrangement of nad3-nad1 relative to other parasitic flatworm mt genomes so far sequenced. Gene order and tRNA secondary structure changes (loss and acquisition of the DHU and/or TPsiC arms of trnC, trnF, and trnR) between mitochondrial genomes of these and other (digenean and cestode) flatworms were inferred by character mapping onto a phylogeny estimated from nuclear small subunit rRNA gene sequences of these same species, in order to find additional rare genomic changes suitable as synapomorphies. Denser and wider taxon sampling of mt genomes across the Platyhelminthes will validate these putative characters.     

Pattern and timing of diversification of Cetartiodactyla (Mammalia, Laurasiatheria), as revealed by a comprehensive analysis of mitochondrial genomes

The order Cetartiodactyla includes cetaceans (whales, dolphins and porpoises) that are found in all oceans and seas, as well as in some rivers, and artiodactyls (ruminants, pigs, peccaries, hippos, camels and llamas) that are present on all continents, except Antarctica and until recent invasions, Australia. There are currently 332 recognized cetartiodactyl species, which are classified into 132 genera and 22 families. Most phylogenetic studies have focused on deep relationships, and no comprehensive time-calibrated tree for the group has been published yet. In this study, 128 new complete mitochondrial genomes of Cetartiodactyla were sequenced and aligned with those extracted from nucleotide databases. Our alignment includes 14,902 unambiguously aligned nucleotide characters for 210 taxa, representing 183 species, 107 genera, and all cetartiodactyl families. Our mtDNA data produced a statistically robust tree, which is largely consistent with previous classifications. However, a few taxa were found to be para- or polyphyletic, including the family Balaenopteridae, as well as several genera and species. Accordingly, we propose several taxonomic changes in order to render the classification compatible with our molecular phylogeny. In some cases, the results can be interpreted as possible taxonomic misidentification or evidence for mtDNA introgression. The existence of three new cryptic species of Ruminantia should therefore be confirmed by further analyses using nuclear data. We estimate divergence times using Bayesian relaxed molecular clock models. The deepest nodes appeared very sensitive to prior assumptions leading to unreliable estimates, primarily because of the misleading effects of rate heterogeneity, saturation and divergent outgroups. In addition, we detected that Whippomorpha contains slow-evolving taxa, such as large whales and hippos, as well as fast-evolving taxa, such as river dolphins. Our results nevertheless indicate that the evolutionary history of cetartiodactyls was punctuated by four main phases of rapid radiation during the Cenozoic era: the sudden occurrence of the three extant lineages within Cetartiodactyla (Cetruminantia, Suina and Tylopoda); the basal diversification of Cetacea during the Early Oligocene; and two radiations that involve Cetacea and Pecora, one at the Oligocene/Miocene boundary and the other in the Middle Miocene. In addition, we show that the high species diversity now observed in the families Bovidae and Cervidae accumulated mainly during the Late Miocene and Plio-Pleistocene.     

Complete mitochondrial genome of compactin-producing fungus Penicillium solitum and comparative analysis of Trichocomaceae mitochondrial genomes

We determined the complete mitochondrial genome sequence of the compactin-producing fungus Penicillium solitum strain 20-01. The 28?601-base pair circular-mapping DNA molecule encodes a characteristic set of mitochondrial proteins and RNA genes and is intron-free. All 46 protein- and RNA-encoding genes are located on one strand and apparently transcribed in one direction. Comparative analysis of this mtDNA and previously sequenced but unannotated mitochondrial genomes of several medically and industrially important species of the Aspergillus/Penicillium group revealed their extensive similarity in terms of size, gene content and sequence, which is also reflected in the almost perfect conservation of mitochondrial gene order in Penicillium and Aspergillus. Phylogenetic analysis based on concatenated mitochondrial protein sequences confirmed the monophyletic origin of Eurotiomycetes.     

Angiosperm mitochondrial genomes and mutations

Flowering plants harbor the largest mitochondrial genomes reported so far. At present, the nucleotide sequences of 15 mitochondrial genomes from seven angiosperm species are available, making detailed comparative analysis feasible. The gene content is variable among the species, but the most striking feature is the fluidity of intergenic regions, where species-specific sequences predominate. Additionally, angiosperm mitochondrial genomes, even within a species, show a remarkable amount of rearrangement. We also review mitochondrial mutants in angiosperms from a genomic viewpoint, and discuss how they have arisen. The involvement of nuclear genes in mitochondrial genome stability and organization is currently being revealed through the analysis of mutants.     

A comparison of three fission yeast mitochondrial genomes

The fission yeasts are members of the fungal order Schizosaccharomycetales, a candidate deep-diverging group within Ascomycota. Although a great deal of molecular information is available from Schizosaccharomyces pombe, a model eukaryote, very little is available from other members of this group. In order to better characterize mitochondrial genome evolution in this fungal lineage, the mitochondrial DNA (mtDNA) of two additional fission yeasts, Schizosaccharomyces octosporus and Schizosaccharomyces japonicus var. japonicus, was sequenced. Whereas the mtDNA of S.pombe is only 19 431 bp, the mtDNA of S.octosporus is 44 227 bp, and that of S.japonicus var. japonicus is over 80 kb. The size variation of these mtDNAs is due largely to non-coding regions. The gene content in the latter two mtDNAs is almost identical to that of the completely sequenced S.pombe mtDNA, which encodes 25 tRNA species, the large and small mitochondrial ribosomal RNAs (rnl and rns), the RNA component of mitochondrial RNaseP (rnpB), mitochondrial small subunit ribosomal protein 3 (rps3), cytochrome oxidase subunits 1, 2 and 3 (cox1, cox2 and cox3) and ATP-synthase subunits 6, 8 and 9 (atp6, atp8 and atp9). However, trnI2(cau) (C modified to lysidine) is absent in the S.octosporus mtDNA, as are corresponding ATA codons in its protein-coding genes, and rps3 and rnpB are not found in the mtDNA of S.japonicus var. japonicus. The mtDNA of S.octosporus contains five double hairpin elements, the first report of these elements in an ascomycete. This study provides further evidence in favor of the mobility of these elements, and supports their role in mitochondrial genome rearrangement. The results of our phylogenetic analysis support the monophyly of the Schizosaccharomycetales, but question their grouping within the Archiascomycota.     

Phylogenetic analysis of diprotodontian marsupials based on complete mitochondrial genomes

Australidelphia is the cohort, originally named by Szalay, of all Australian marsupials and the South American Dromiciops. A lot of mitochondria and nuclear genome studies support the hypothesis of a monophyly of Australidelphia, but some familial relationships in Australidelphia are still unclear. In particular, the familial relationships among the order Diprotodontia (koala, wombat, kangaroos and possums) are ambiguous. These Diprotodontian families are largely grouped into two suborders, Vombatiformes, which contains Phascolarctidae (koala) and Vombatidae (wombat), and Phalangerida, which contains Macropodidae, Potoroidae, Phalangeridae, Petauridae, Pseudocheiridae, Acrobatidae, Tarsipedidae and Burramyidae. Morphological evidence and some molecular analyses strongly support monophyly of the two families in Vombatiformes. The monophyly of Phalangerida as well as the phylogenetic relationships of families in Phalangerida remains uncertain, however, despite searches for morphological synapomorphy and mitochondrial DNA sequence analyses. Moreover, phylogenetic relationships among possum families (Phalangeridae, Petauridae, Pseudocheiridae, Acrobatidae, Tarsipedidae and Burramyidae) as well as a sister group of Macropodoidea (Macropodidae and Potoroidae) remain unclear. To evaluate familial relationships among Dromiciops and Australian marsupials as well as the familial relationships in Diprotodontia, we determined the complete mitochondrial sequence of six Diprotodontian species. We used Maximum Likelihood analyses with concatenated amino acid and codon sequences of 12 mitochondrial protein genomes. Our analysis of mitochondria amino acid sequence supports monophyly of Australian marsupials+Dromiciops and monophyly of Phalangerida. The close relatedness between Macropodidae and Phalangeridae is also weakly supported by our analysis.     

In silico analysis of SSRs in mitochondrial genomes of plants

Simple sequence repeats (SSRs) or microsatellites constitute a countable portion of genomes. However, the significance of SSRs in organelle genomes has not been completely understood. The availability of organelle genome sequences allows us to understand the organization of SSRs in their genic and intergenic regions. In the current study we surveyed the patterns of SSRs in mitochondrial genomes of different taxa of plants. A total of 16 mitochondrial genomes, from algae to angiosperms, have been considered to analyze the pattern of simple sequence repeats present in them. Based on study, the mononucleotide repeats of A/T were found to be more prevalent in mitochondrial genomes over other repeat types. The dinucleotides repeats, TA/AT, were the second most numerous, whereas tri-, tetra-, and pentanucleotide repeats were in less number and present in intronic or intergenic portions only. Mononucleotide repeats prevailed in protein-coding exonic portions of all organisms. These results indicates that microsatellite pattern in mitochondrial genomes is different from nuclear genomes and also focuses on organization and diversity at SSR locuses in mitochondrial genomes. This is the novel report of microsatellite polymorphism in plant mitochondrion on whole genome level.     

Four new mitochondrial genomes and the increased stability of evolutionary trees of mammals from improved taxon sampling

We have sequenced four new mitochondrial genomes to improve the stability of the tree for placental mammals; they are two insectivores (a gymnure, Echinosorex gymnurus and Formosan shrew Soriculus fumidus); a Formosan lesser horseshoe bat (Rhinolophus monoceros); and the New Zealand fur seal (Arctocephalus forsteri). A revision to the hedgehog sequence (Erinaceus europaeus) is also reported. All five are from the Laurasiatheria grouping of eutherian mammals. On this new data set there is a strong tendency for the hedgehog and its relative, the gymnure, to join with the other Laurasiatherian insectivores (mole and shrews). To quantify the stability of trees from this data we define, based on nuclear sequences, a major four-way split in Laurasiatherians. This ([Xenarthra, Afrotheria], [Laurasiatheria, Supraprimates]) split is also found from mitochondrial genomes using either protein-coding or RNA (rRNA and tRNA) data sets. The high similarity of the mitochondrial and nuclear-derived trees allows a quantitative estimate of the stability of trees from independent data sets, as detected from a triplet Markov analysis. There are significant changes in the mutational processes within placental mammals that are ignored by current tree programs. On the basis of our quantitative results, we expect the evolutionary tree for mammals to be resolved quickly, and this will allow other problems to be solved.     

The complete mitochondrial genomes of sixteen ardeid birds revealing the evolutionary process of the gene rearrangements

Background

The animal mitochondrial genome is generally considered to be under selection for both compactness and gene order conservation. As more mitochondrial genomes are sequenced, mitochondrial duplications and gene rearrangements have been frequently identified among diverse animal groups. Although several mechanisms of gene rearrangement have been proposed thus far, more observational evidence from major taxa is needed to validate specific mechanisms. In the current study, the complete mitochondrial DNA of sixteen bird species from the family Ardeidae was sequenced and the evolution of mitochondrial gene rearrangements was investigated. The mitochondrial genomes were then used to review the phylogenies of these ardeid birds.

Results

The complete mitochondrial genome sequences of the sixteen ardeid birds exhibited four distinct mitochondrial gene orders in which two of them, named as “duplicate tRNA^Glu–CR” and “duplicate tRNA^Thr–tRNA^Pro and CR”, were newly discovered. These gene rearrangements arose from an evolutionary process consistent with the tandem duplication - random loss model (TDRL). Additionally, duplications in these gene orders were near identical in nucleotide sequences within each individual, suggesting that they evolved in concert. Phylogenetic analyses of the sixteen ardeid species supported the idea that Ardea ibis, Ardea modesta and Ardea intermedia should be classified as genus Ardea, and Ixobrychus flavicollis as genus Ixobrychus, and indicated that within the subfamily Ardeinae, Nycticorax nycticorax is closely related to genus Egretta and that Ardeola bacchus and Butorides striatus are closely related to the genus Ardea.

Conclusions

The duplicate tRNAThr–CR gene order is found in most ardeid lineages, suggesting this gene order is the ancestral pattern within these birds and persisted in most lineages via concerted evolution. In two independent lineages, when the concerted evolution stopped in some subsections due to the accumulation of numerous substitutions and deletions, the duplicate tRNAThr–CR gene order was transformed into three other gene orders. The phylogenetic trees produced from concatenated rRNA and protein coding genes have high support values in most nodes, indicating that the mitochondrial genome sequences are promising markers for resolving the phylogenetic issues of ardeid birds when more taxa are added.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-573) contains supplementary material, which is available to authorized users.     

Twenty-fold difference in evolutionary rates between the mitochondrial and plastid genomes of species with secondary red plastids

Within plastid-bearing species, the relative rates of evolution between mitochondrial and plastid genomes are poorly studied, but for the few lineages in which they have been explored, including land plants and green algae, the mitochondrial DNA mutation rate is nearly always estimated to be lower than or equal to that of the plastid DNA. Here, we show that in protists from three distinct lineages with secondary, red algal-derived plastids, the opposite is true: their mitochondrial genomes are evolving 5-30 times faster than their plastid genomes, even when the plastid is nonphotosynthetic. These findings have implications for understanding the origins and evolution of organelle genome architecture and the genes they encode.     

Four new avian mitochondrial genomes help get to basic evolutionary questions in the late cretaceous

Good phylogenetic trees are required to test hypotheses about evolutionary processes. We report four new avian mitochondrial genomes, which together with an improved method of phylogenetic analysis for vertebrate mt genomes give results for three questions in avian evolution. The new mt genomes are: magpie goose (Anseranas semipalmata), an owl (morepork, Ninox novaeseelandiae); a basal passerine (rifleman, or New Zealand wren, Acanthisitta chloris); and a parrot (kakapo or owl-parrot, Strigops habroptilus). The magpie goose provides an important new calibration point for avian evolution because the well-studied Presbyornis fossils are on the lineage to ducks and geese, after the separation of the magpie goose. We find, as with other animal mitochondrial genomes, that RY-coding is helpful in adjusting for biases between pyrimidines and between purines. When RY-coding is used at third positions of the codon, the root occurs between paleognath and neognath birds (as expected from morphological and nuclear data). In addition, passerines form a relatively old group in Neoaves, and many modern avian lineages diverged during the Cretaceous. Although many aspects of the avian tree are stable, additional taxon sampling is required.