The complete mitochondrial genomes of the sea lily Gymnocrinus richeri and the feather star Phanogenia gracilis: Signature nucleotide bias and unique nad4L gene rearrangement within crinoids

Mitochondrial Genomic Landscape: A Portrait of the Mitochondrial Genome 40 Years after the First Complete Sequence

Notwithstanding the initial claims of general conservation, mitochondrial genomes are a largely heterogeneous set of organellar chromosomes which displays a bewildering diversity in terms of structure, architecture, gene content, and functionality. The mitochondrial genome is typically described as a single chromosome, yet many examples of multipartite genomes have been found (for example, among sponges and diplonemeans); the mitochondrial genome is typically depicted as circular, yet many linear genomes are known (for example, among jellyfish, alveolates, and apicomplexans); the chromosome is normally said to be “small”, yet there is a huge variation between the smallest and the largest known genomes (found, for example, in ctenophores and vascular plants, respectively); even the gene content is highly unconserved, ranging from the 13 oxidative phosphorylation-related enzymatic subunits encoded by animal mitochondria to the wider set of mitochondrial genes found in jakobids. In the present paper, we compile and describe a large database of 27,873 mitochondrial genomes currently available in GenBank, encompassing the whole eukaryotic domain. We discuss the major features of mitochondrial molecular diversity, with special reference to nucleotide composition and compositional biases; moreover, the database is made publicly available for future analyses on the MoZoo Lab GitHub page.

Gene expression analysis of three homeobox genes throughout early and late development of a feather star Anneissia japonica

Crinoids are considered as the most basal extant echinoderms. They retain aboral nervous system with a nerve center, which has been degraded in the eleutherozoan echinoderms. To investigate the evolution of patterning of the nervous systems in crinoids, we examined temporal and spatial expression patterns of three neural patterning-related homeobox genes, six3, pax6, and otx, throughout the development of a feather star Anneissia japonica. These genes were involved in the patterning of endomesodermal tissues instead of the ectodermal neural tissues in the early planktonic stages. In the stages after larval attachment, the expression of these genes was mainly observed in the podia and the oral nervous systems instead of the aboral nerve center. Our results indicate the involvement of these three genes in the formation of oral nervous system in the common ancestor of the echinoderms and suggest that the aboral nerve center is not evolutionally related to the brain of other bilaterians.

Complete mitochondrial genome of the crinoid echinoderm, Florometra species (Echinodermata, Crinoidea)

In this study, we report the sequence of the mitochondrial genome (mitogenome) of the crinoid echinoderm, Florometra species (Echinodermata, Crinoidea). The complete mitogenome of Florometra sp. was 15,792 base pairs long and was composed of 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and three regions of unassigned sequence (UAS) including one putative control region. Unique nucleotide composition, a clear positive bias for T with an apparent loss of C in PCGs as observed in the Crinoidea mitogenomes, was also seen in the Florometra sp. mitogenome (45% T, 12% C). Phylogenetic analysis with the concatenated nucleotide sequences of entire PCGs of echinoderms confirmed that Florometra sp. is highly related to F. serratissima within the family Crinoidea.

Characterization and phylogenetic analysis of the complete mitochondrial genome of Stephnometra indica (Pelmatozoa: Crinoidea)

The crinoid Stephnometra indica inhabits coral reefs and surrounding waters in the western Pacific Ocean. In this study, the complete mitochondrial genome of S. indica (15,878 bp) was confirmed to contain 13 protein-coding genes, 22 tRNA genes, and 2 rRNA genes. The gene order of S. indica was identical to that of previously published Crinoidea. The nucleotide composition of the S. indica mitogenome was biased toward A + T nucleotide 73.68%. There were four unassigned regions (UASs) found in S. indica and all four UAS regions were AT-biased. The phylogenetic relationship illustrated that S. indica is more closely related to Antedon mediterranea within the family Crinoidea.

The first complete mitochondrial genome of the Mariana Trench Freyastera benthophila (Asteroidea: Brisingida: Brisingidae) allows insights into the deep-sea adaptive evolution of Brisingida

Starfish (phylum Echinodermata) are ecologically important and diverse members of marine ecosystems in all of the world's oceans, from the shallow water to the hadal zone. The deep sea is recognized as an extremely harsh environment on earth. In this study, we present the mitochondrial genome sequence of Mariana Trench starfish Freyastera benthophila, and this study is the first to explore in detail the mitochondrial genome of a deep-sea member of the order Brisingida. Similar to other starfish, it contained 13 protein-coding genes, two ribosomal RNA genes, and 22 transfer RNA genes (duplication of two tRNAs: trnL and trnS). Twenty-two of these genes are encoded on the positive strand, while the other 15 are encoded on the negative strand. The gene arrangement was identical to those of sequenced starfish. Phylogenetic analysis showed the deep-sea Brisingida as a sister taxon to the traditional members of the Asteriidae. Positive selection analysis indicated that five residues (8 N and 16 I in atp8, 47 D and 196 V in nad2, 599 N in nad5) were positively selected sites with high posterior probabilities. Compared these features with shallow sea starfish, we predict that variation specifically in atp8, nad2, and nad5 may play an important role in F. benthophila's adaptation to deep-sea environment.

A complete logical approach to resolve the evolution and dynamics of mitochondrial genome in bilaterians

Investigating how recombination might modify gene order during the evolution has become a routine part of mitochondrial genome analysis. A new method of genomic maps analysis based on formal logic is described. The purpose of this method is to 1) use mitochondrial gene order of current taxa as datasets 2) calculate rearrangements between all mitochondrial gene orders and 3) reconstruct phylogenetic relationships according to these calculated rearrangements within a tree under the assumption of maximum parsimony. Unlike existing methods mainly based on the probabilistic approach, the main strength of this new approach is that it calculates all the exact tree solutions with completeness and provides logical consequences as highly robust results. Moreover, this method infers all possible hypothetical ancestors and reconstructs character states for all internal nodes of the trees. We started by testing our method using the deuterostomes as a study case. Then, with sponges as an outgroup, we investigated the evolutionary history of mitochondrial genomes of 47 bilaterian phyla and emphasised the peculiar case of chaetognaths. This pilot work showed that the use of formal logic in a hypothetico-deductive background such as phylogeny (where experimental testing of hypotheses is impossible) is very promising to explore mitochondrial gene order in deuterostomes and should be applied to many other bilaterian clades.

Multiple reversals of strand asymmetry in molluscs mitochondrial genomes, and consequences for phylogenetic inferences

Strand asymmetry in nucleotide composition is a remarkable feature of animal mitochondrial genomes. The strand-specific bias in the nucleotide composition of the mtDNA has been known to be highly problematic for phylogenetic analyses. Here, the strand asymmetry was compared across 140 mollusc species and analyzed for a mtDNA fragment including twelve protein-coding genes. The analyses show that almost all species in Gastropoda (except Heterobranchia) and all species in Bivalvia present reversals of strand bias. The skew values on individual genes for all codon positions (P₁₂₃), third codon positions (P₃), and fourfold redundant third codon positions (P_4FD) indicated that CG skews are the best indicators of strand asymmetry. The differences in the patterns of strand asymmetry significantly influenced the amino acid composition of the encoded proteins. These biases are most striking for the amino acids Valine, Cysteine, Asparagine and Threonines, which appear to have evolved asymmetrical exchanges in response to shifts in nucleotide composition. Molluscs with strong variability of genome architectures (ARs) are usually characterized by a reversal of the usual strand bias. Phylogenetic analyses show that reversals of asymmetric mutational constraints have consequences on the phylogenetic inferences, as taxa characterized by reverse strand bias (Heterobranchia and Bivalvia) tend to group together due to long-branch attraction (LBA) artifacts. Neutral Transitions Excluded (NTE) model did not overcome the problem of heterogeneous biases present in molluscs mt genomes, suggested it may not be appropriate for molluscs mt genome data. Further refinement phylogenetic models may help us better understand internal relationships among these diverse organisms.

Fixed, free, and fixed: the fickle phylogeny of extant Crinoidea (Echinodermata) and their Permian-Triassic origin

Although the status of Crinoidea (sea lilies and featherstars) as sister group to all other living echinoderms is well-established, relationships among crinoids, particularly extant forms, are debated. All living species are currently placed in Articulata, which is generally accepted as the only crinoid group to survive the Permian-Triassic extinction event. Recent classifications have recognized five major extant taxa: Isocrinida, Hyocrinida, Bourgueticrinina, Comatulidina and Cyrtocrinida, plus several smaller groups with uncertain taxonomic status, e.g., Guillecrinus, Proisocrinus and Caledonicrinus. Here we infer the phylogeny of extant Crinoidea using three mitochondrial genes and two nuclear genes from 59 crinoid terminals that span the majority of extant crinoid diversity. Although there is poor support for some of the more basal nodes, and some tree topologies varied with the data used and mode of analysis, we obtain several robust results. Cyrtocrinida, Hyocrinida, Isocrinida are all recovered as clades, but two stalked crinoid groups, Bourgueticrinina and Guillecrinina, nest among the featherstars, lending support to an argument that they are paedomorphic forms. Hence, they are reduced to families within Comatulida. Proisocrinus is clearly shown to be part of Isocrinida, and Caledonicrinus may not be a bourgueticrinid. Among comatulids, tree topologies show little congruence with current taxonomy, indicating that much systematic revision is required. Relaxed molecular clock analyses with eight fossil calibration points recover Articulata with a median date to the most recent common ancestor at 231-252mya in the Middle to Upper Triassic. These analyses tend to support the hypothesis that the group is a radiation from a small clade that passed through the Permian-Triassic extinction event rather than several lineages that survived. Our tree topologies show various scenarios for the evolution of stalks and cirri in Articulata, so it is clear that further data and taxon sampling are needed to recover a more robust phylogeny of the group.

Current Status of Echinoderm Genome Analysis - What do we Know?

Echinoderms have long served as model organisms for a variety of biological research, especially in the field of developmental biology. Although the genome of the purple sea urchin Strongylocentrotus purpuratus has been sequenced, it is the only echinoderm whose whole genome sequence has been reported. Nevertheless, data is rapidly accumulating on the chromosomes and genomic sequences of all five classes of echinoderms, including the mitochondrial genomes and Hox genes. This blossoming new data will be essential for estimating the phylogenetic relationships among echinoderms, and also to examine the underlying mechanisms by which the diverse morphologies of echinoderms have arisen.

Potential pitfalls of modelling ribosomal RNA data in phylogenetic tree reconstruction: evidence from case studies in the Metazoa

Background

Failure to account for covariation patterns in helical regions of ribosomal RNA (rRNA) genes has the potential to misdirect the estimation of the phylogenetic signal of the data. Furthermore, the extremes of length variation among taxa, combined with regional substitution rate variation can mislead the alignment of rRNA sequences and thus distort subsequent tree reconstructions. However, recent developments in phylogenetic methodology now allow a comprehensive integration of secondary structures in alignment and tree reconstruction analyses based on rRNA sequences, which has been shown to correct some of these problems. Here, we explore the potentials of RNA substitution models and the interactions of specific model setups with the inherent pattern of covariation in rRNA stems and substitution rate variation among loop regions.

Results

We found an explicit impact of RNA substitution models on tree reconstruction analyses. The application of specific RNA models in tree reconstructions is hampered by interaction between the appropriate modelling of covarying sites in stem regions, and excessive homoplasy in some loop regions. RNA models often failed to recover reasonable trees when single-stranded regions are excessively homoplastic, because these regions contribute a greater proportion of the data when covarying sites are essentially downweighted. In this context, the RNA6A model outperformed all other models, including the more parametrized RNA7 and RNA16 models.

Conclusions

Our results depict a trade-off between increased accuracy in estimation of interdependencies in helical regions with the risk of magnifying positions lacking phylogenetic signal. We can therefore conclude that caution is warranted when applying rRNA covariation models, and suggest that loop regions be independently screened for phylogenetic signal, and eliminated when they are indistinguishable from random noise. In addition to covariation and homoplasy, other factors, like non-stationarity of substitution rates and base compositional heterogeneity, can disrupt the signal of ribosomal RNA data. All these factors dictate sophisticated estimation of evolutionary pattern in rRNA data, just as other molecular data require similarly complicated (but different) corrections.

The enigmatic mitochondrial genome of Rhabdopleura compacta (Pterobranchia) reveals insights into selection of an efficient tRNA system and supports monophyly of Ambulacraria

Background

The Hemichordata comprises solitary-living Enteropneusta and colonial-living Pterobranchia, sharing morphological features with both Chordata and Echinodermata. Despite their key role for understanding deuterostome evolution, hemichordate phylogeny is controversial and only few molecular data are available for phylogenetic analysis. Furthermore, mitochondrial sequences are completely lacking for pterobranchs. Therefore, we determined and analyzed the complete mitochondrial genome of the pterobranch Rhabdopleura compacta to elucidate deuterostome evolution. Thereby, we also gained important insights in mitochondrial tRNA evolution.

Results

The mitochondrial DNA of Rhabdopleura compacta corresponds in size and gene content to typical mitochondrial genomes of metazoans, but shows the strongest known strand-specific mutational bias in the nucleotide composition among deuterostomes with a very GT-rich main-coding strand. The order of the protein-coding genes in R. compacta is similar to that of the deuterostome ground pattern. However, the protein-coding genes have been highly affected by a strand-specific mutational pressure showing unusual codon frequency and amino acid composition. This composition caused extremely long branches in phylogenetic analyses. The unusual codon frequency points to a selection pressure on the tRNA translation system to codon-anticodon sequences of highest versatility instead of showing adaptations in anticodon sequences to the most frequent codons. Furthermore, an assignment of the codon AGG to Lysine has been detected in the mitochondrial genome of R. compacta, which is otherwise observed only in the mitogenomes of some arthropods. The genomes of these arthropods do not have such a strong strand-specific bias as found in R. compacta but possess an identical mutation in the anticodon sequence of the tRNA_Lys.

Conclusion

A strong reversed asymmetrical mutational constraint in the mitochondrial genome of Rhabdopleura compacta may have arisen by an inversion of the replication direction and adaptation to this bias in the protein sequences leading to an enigmatic mitochondrial genome. Although, phylogenetic analyses of protein coding sequences are hampered, features of the tRNA system of R. compacta support the monophyly of Ambulacraria. The identical reassignment of AGG to Lysine in two distinct groups may have occurred by convergent evolution in the anticodon sequence of the tRNA_Lys.

Mutation and selection on the wobble nucleotide in tRNA anticodons in marine bivalve mitochondrial genomes

Background

Animal mitochondrial genomes typically encode one tRNA for each synonymous codon family, so that each tRNA anticodon essentially has to wobble to recognize two or four synonymous codons. Several factors have been hypothesized to determine the nucleotide at the wobble site of a tRNA anticodon in mitochondrial genomes, such as the codon-anticodon adaptation hypothesis, the wobble versatility hypothesis, the translation initiation and elongation conflict hypothesis, and the wobble cost hypothesis.

Principal Findings

In this study, we analyzed codon usage and tRNA anticodon wobble sites of 29 marine bivalve mitochondrial genomes to evaluate features of the wobble nucleotides in tRNA anticodons. The strand-specific mutation bias favors G and T on the H strand in all the 29 marine bivalve mitochondrial genomes. A bias favoring G and T is also visible in the third codon positions of protein-coding genes and the wobble sites of anticodons, rejecting that codon usage bias drives the wobble sites of tRNA anticodons or tRNA anticodon bias drives the evolution of codon usage. Almost all codon families (98.9%) from marine bivalve mitogenomes support the wobble versatility hypothesis. There are a few interesting exceptions involving tRNA^Trp with an anticodon CCA fixed in Pectinoida species, tRNA^Ser with a GCU anticodon fixed in Mytiloida mitogenomes, and the uniform anticodon CAU of tRNA^Met translating the AUR codon family.

Conclusions/Significance

These results demonstrate that most of the nucleotides at the wobble sites of tRNA anticodons in marine bivalve mitogenomes are determined by wobble versatility. Other factors such as the translation initiation and elongation conflict, and the cost of wobble translation may contribute to the determination of the wobble nucleotide in tRNA anticodons. The finding presented here provides valuable insights into the previous hypotheses of the wobble nucleotide in tRNA anticodons by adding some new evidence.

New views on strand asymmetry in insect mitochondrial genomes

Strand asymmetry in nucleotide composition is a remarkable feature of animal mitochondrial genomes. Understanding the mutation processes that shape strand asymmetry is essential for comprehensive knowledge of genome evolution, demographical population history and accurate phylogenetic inference. Previous studies found that the relative contributions of different substitution types to strand asymmetry are associated with replication alone or both replication and transcription. However, the relative contributions of replication and transcription to strand asymmetry remain unclear. Here we conducted a broad survey of strand asymmetry across 120 insect mitochondrial genomes, with special reference to the correlation between the signs of skew values and replication orientation/gene direction. The results show that the sign of GC skew on entire mitochondrial genomes is reversed in all species of three distantly related families of insects, Philopteridae (Phthiraptera), Aleyrodidae (Hemiptera) and Braconidae (Hymenoptera); the replication-related elements in the A+T-rich regions of these species are inverted, confirming that reversal of strand asymmetry (GC skew) was caused by inversion of replication origin; and finally, the sign of GC skew value is associated with replication orientation but not with gene direction, while that of AT skew value varies with gene direction, replication and codon positions used in analyses. These findings show that deaminations during replication and other mutations contribute more than selection on amino acid sequences to strand compositions of G and C, and that the replication process has a stronger affect on A and T content than does transcription. Our results may contribute to genome-wide studies of replication and transcription mechanisms.

Comparative mitogenomics of Braconidae (Insecta: Hymenoptera) and the phylogenetic utility of mitochondrial genomes with special reference to Holometabolous insects

BACKGROUND

Animal mitochondrial genomes are potential models for molecular evolution and markers for phylogenetic and population studies. Previous research has shown interesting features in hymenopteran mitochondrial genomes. Here, we conducted a comparative study of mitochondrial genomes of the family Braconidae, one of the largest families of Hymenoptera, and assessed the utility of mitochondrial genomic data for phylogenetic inference at three different hierarchical levels, i.e., Braconidae, Hymenoptera, and Holometabola.

RESULTS

Seven mitochondrial genomes from seven subfamilies of Braconidae were sequenced. Three of the four sequenced A+T-rich regions are shown to be inverted. Furthermore, all species showed reversal of strand asymmetry, suggesting that inversion of the A+T-rich region might be a synapomorphy of the Braconidae. Gene rearrangement events occurred in all braconid species, but gene rearrangement rates were not taxonomically correlated. Most rearranged genes were tRNAs, except those of Cotesia vestalis, in which 13 protein-coding genes and 14 tRNA genes changed positions or/and directions through three kinds of gene rearrangement events. Remote inversion is posited to be the result of two independent recombination events. Evolutionary rates were lower in species of the cyclostome group than those of noncyclostomes. Phylogenetic analyses based on complete mitochondrial genomes and secondary structure of rrnS supported a sister-group relationship between Aphidiinae and cyclostomes. Many well accepted relationships within Hymenoptera, such as paraphyly of Symphyta and Evaniomorpha, a sister-group relationship between Orussoidea and Apocrita, and monophyly of Proctotrupomorpha, Ichneumonoidea and Aculeata were robustly confirmed. New hypotheses, such as a sister-group relationship between Evanioidea and Aculeata, were generated. Among holometabolous insects, Hymenoptera was shown to be the sister to all other orders. Mecoptera was recovered as the sister-group of Diptera. Neuropterida (Neuroptera + Megaloptera), and a sister-group relationship with (Diptera + Mecoptera) were supported across all analyses.

CONCLUSIONS

Our comparative studies indicate that mitochondrial genomes are a useful phylogenetic tool at the ordinal level within Holometabola, at the superfamily within Hymenoptera and at the subfamily level within Braconidae. Variation at all of these hierarchical levels suggests that the utility of mitochondrial genomes is likely to be a valuable tool for systematics in other groups of arthropods.

Nervous system development of two crinoid species, the sea lily Metacrinus rotundus and the feather star Oxycomanthus japonicus

Nervous system development in echinoderms has been well documented, especially for sea urchins and starfish. However, that of crinoids, the most basal group of extant echinoderms, has been poorly studied due to difficulties in obtaining their larvae. In this paper, we report nervous system development from two species of crinoids, from hatching to late doliolaria larvae in the sea lily Metacrinus rotundus and from hatching to cystidean stages after settlement in the feather star Oxycomanthus japonicus. The two species showed a similar larval nervous system pattern with an extensive anterior larval ganglion. The ganglion was similar to that in sea urchins which is generally regarded as derived. In contrast with other echinoderm and hemichordate larvae, synaptotagmin antibody 1E11 failed to reveal ciliary band nerve tracts. Basiepithelial nerve cells formed a net-like structure in the M. rotundus doliolaria larvae. In O. japonicus, the larval ganglion was still present 1 day after settlement when the adult nervous system began to appear inside the crown. Stalk nerves originated from the crown and extended down the stalk, but had no connections with the remaining larval ganglion at the base of the stalk. The larval nervous system was not incorporated into the adult nervous system, and the larval ganglion later disappeared. The aboral nerve center, the dominant nervous system in adult crinoids, was formed at the early cystidean stage, considerably earlier than previously suggested. Through comparisons with nervous system development in other ambulacraria, we suggest the possible nervous system development pattern of the echinoderm ancestor and provide new implications on the evolutionary history of echinoderm life cycles.

Accurate and efficient reconstruction of deep phylogenies from structured RNAs

Ribosomal RNA (rRNA) genes are probably the most frequently used data source in phylogenetic reconstruction. Individual columns of rRNA alignments are not independent as a consequence of their highly conserved secondary structures. Unless explicitly taken into account, these correlation can distort the phylogenetic signal and/or lead to gross overestimates of tree stability. Maximum likelihood and Bayesian approaches are of course amenable to using RNA-specific substitution models that treat conserved base pairs appropriately, but require accurate secondary structure models as input. So far, however, no accurate and easy-to-use tool has been available for computing structure-aware alignments and consensus structures that can deal with the large rRNAs. The RNAsalsa approach is designed to fill this gap. Capitalizing on the improved accuracy of pairwise consensus structures and informed by a priori knowledge of group-specific structural constraints, the tool provides both alignments and consensus structures that are of sufficient accuracy for routine phylogenetic analysis based on RNA-specific substitution models. The power of the approach is demonstrated using two rRNA data sets: a mitochondrial rRNA set of 26 Mammalia, and a collection of 28S nuclear rRNAs representative of the five major echinoderm groups.

The mitochondrial genome structure of Xenoturbella bocki (phylum Xenoturbellida) is ancestral within the deuterostomes

Background

Mitochondrial genome comparisons contribute in multiple ways when inferring animal relationships. As well as primary sequence data, rare genomic changes such as gene order, shared gene boundaries and genetic code changes, which are unlikely to have arisen through convergent evolution, are useful tools in resolving deep phylogenies. Xenoturbella bocki is a morphologically simple benthic marine worm recently found to belong among the deuterostomes. Here we present analyses comparing the Xenoturbella bocki mitochondrial gene order, genetic code and control region to those of other metazoan groups.

Results

The complete mitochondrial genome sequence of Xenoturbella bocki was determined. The gene order is most similar to that of the chordates and the hemichordates, indicating that this conserved mitochondrial gene order might be ancestral to the deuterostome clade. Using data from all phyla of deuterostomes, we infer the ancestral mitochondrial gene order for this clade. Using inversion and breakpoint analyses of metazoan mitochondrial genomes, we test conflicting hypotheses for the phylogenetic placement of Xenoturbella and find a closer affinity to the hemichordates than to other metazoan groups. Comparative analyses of the control region reveal similarities in the transcription initiation and termination sites and origin of replication of Xenoturbella with those of the vertebrates. Phylogenetic analyses of the mitochondrial sequence indicate a weakly supported placement as a basal deuterostome, a result that may be the effect of compositional bias.

Conclusion

The mitochondrial genome of Xenoturbella bocki has a very conserved gene arrangement in the deuterostome group, strikingly similar to that of the hemichordates and the chordates, and thus to the ancestral deuterostome gene order. Similarity to the hemichordates in particular is suggested by inversion and breakpoint analysis. Finally, while phylogenetic analyses of the mitochondrial sequences support a basal deuterostome placement, support for this decreases with the use of more sophisticated models of sequence evolution.

Arachnid relationships based on mitochondrial genomes: asymmetric nucleotide and amino acid bias affects phylogenetic analyses

Phylogenetic analyses based on mitochondrial DNA have yielded widely differing relationships among members of the arthropod lineage Arachnida, depending on the nucleotide coding schemes and models of evolution used. We enhanced taxonomic coverage within the Arachnida greatly by sequencing seven new arachnid mitochondrial genomes from five orders. We then used all 13 mitochondrial protein-coding genes from these genomes to evaluate patterns of nucleotide and amino acid biases. Our data show that two of the six orders of arachnids (spiders and scorpions) have experienced shifts in both nucleotide and amino acid usage in all their protein-coding genes, and that these biases mislead phylogeny reconstruction. These biases are most striking for the hydrophobic amino acids isoleucine and valine, which appear to have evolved asymmetrical exchanges in response to shifts in nucleotide composition. To improve phylogenetic accuracy based on amino acid differences, we tested two recoding methods: (1) removing all isoleucine and valine sites and (2) recoding amino acids based on their physiochemical properties. We find that these methods yield phylogenetic trees that are consistent in their support of ancient intraordinal divergences within the major arachnid lineages. Further refinement of amino acid recoding methods may help us better delineate interordinal relationships among these diverse organisms.

Two circular chromosomes of unequal copy number make up the mitochondrial genome of the rotifer Brachionus plicatilis

The monogonont rotifer Brachionus plicatilis is an emerging model system for a diverse array of questions in limnological ecosystem dynamics, the evolution of sexual recombination, cryptic speciation, and the phylogeny of basal metazoans. We sequenced the complete mitochondrial genome of B. plicatilis sensu strictu NH1L and found that it is composed of 2 circular chromosomes, designated mtDNA-I (11,153 bp) and mtDNA-II (12,672 bp). Hybridization to DNA isolated from mitochondria demonstrated that mtDNA-I is present at 4 times the copy number of mtDNA-II. The only nucleotide similarity between the 2 chromosomes is a 4.9-kbp region of 99.5% identity including a transfer RNA (tRNA) gene and an extensive noncoding region that contains putative D-loop and control sequence. The mtDNA-I chromosome encodes 4 proteins (ATP6, COB, NAD1, and NAD2), 13 tRNAs, and the large and small subunit ribosomal RNAs; mtDNA-II encodes 8 proteins (COX1-3, NAD3-6, and NAD4L) and 9 tRNAs. Gene order is not conserved between B. plicatilis and its closest relative with a sequenced mitochondrial genome, the acanthocephalan Leptorhynchoides thecatus, or other sequenced mitochondrial genomes. Polymerase chain reaction assays and Southern hybridization to DNA from 18 strains of Brachionus suggest that the 2-chromosome structure has been stable for millions of years. The novel organization of the B. plicatilis mitochondrial genome into 2 nearly equal chromosomes of 4-fold different copy number may provide insight into the evolution of metazoan mitochondria and the phylogenetics of rotifers and other basal animal phyla.

Evolution of mitochondrial gene orders in echinoderms

A comprehensive analysis of the mitochondrial gene orders of all previously published and two novel Antedon mediterranea (Crinoidea) and Ophiura albida (Ophiuroidea) complete echinoderm mitochondrial genomes shows that all major types of rearrangement operations are necessary to explain the evolution of mitochondrial genomes. In addition to protein coding genes we include all tRNA genes as well as the control region in our analysis. Surprisingly, 7 of the 16 genomes published in the GenBank database contain misannotations, mostly unannotated tRNAs and/or mistakes in the orientation of tRNAs, which we have corrected here. Although the gene orders of mt genomes appear very different, only 8 events are necessary to explain the evolutionary history of echinoderms with the exception of the ophiuroids. Only two of these rearrangements are inversions, while we identify three tandem-duplication-random-loss events and three transpositions.

Ribosomal RNA genes and deuterostome phylogeny revisited: More cyclostomes, elasmobranchs, reptiles, and a brittle star

This is an expanded study of the relationships among the deuterostome animals based on combined, nearly complete 28S and 18S rRNA genes (>3925 nt.). It adds sequences from 20 more taxa to the approximately 45 sequences used in past studies. Seven of the new taxa were sequenced here (brittle star Ophiomyxa, lizard Anolis, turtle Chrysemys, sixgill shark Hexanchus, electric ray Narcine, Southern Hemisphere lamprey Geotria, and Atlantic hagfish Myxine for 28S), and the other 13 were from GenBank and the literature (from a chicken, dog, rat, human, three lungfishes, and several ray-finned fishes, or Actinopterygii). As before, our alignments were based on secondary structure but did not account for base pairing in the stems of rRNA. The new findings, derived from likelihood-based tree-reconstruction methods and by testing hypotheses with parametric bootstrapping, include: (1) brittle star joins with sea star in the echinoderm clade, Asterozoa; (2) with two hagfishes and two lampreys now available, the cyclostome (jawless) fishes remain monophyletic; (3) Hexanchiform sharks are monophyletic, as Hexanchus groups with the frilled shark, Chlamydoselachus; (4) turtle is the sister taxon of all other amniotes; (5) bird is closer to the lizard than to the mammals; (6) the bichir Polypterus is in a monophyletic Actinopterygii; (7) Zebrafish Danio is the sister taxon of the other two teleosts we examined (trout and perch); (8) the South American and African lungfishes group together to the exclusion of the Australian lungfish. Other findings either upheld those of the previous rRNA-based studies (e.g., echinoderms and hemichordates group as Ambulacraria; orbitostylic sharks; batoids are not derived from any living lineage of sharks) or were obvious (monophyly of mammals, gnathostomes, vertebrates, echinoderms, etc.). Despite all these findings, the rRNA data still fail to resolve the relations among the major groups of deuterostomes (tunicates, Ambulacraria, cephalochordates and vertebrates) and of gnathostomes (chondrichthyans, lungfishes, coelacanth, actinopterygians, amphibians, and amniotes), partly because tunicates and lungfishes are rogue taxa that disrupt the tree. Nonetheless, parametric bootstrapping showed our RNA-gene data are only consistent with these dominant hypotheses: (1) deuterostomes consist of Ambulacraria plus Chordata, with Chordata consisting of tunicates and 'vertebrates plus cephalochordates'; and (2) lungfishes are the closest living relatives of tetrapods.