Gene Order Data Research Articles

Complete mitochondrial genome of Metathelazia capsulata (Pneumospiruridae) and comparison with other Spiruromorpha species.

Metathelazia capsulata (family Pneumospiruridae) is a lungworm parasitizing the bronchi and bronchioles, described in four species of wild carnivores. Very little molecular data are available on this nematode and none on other species of the Pneumospiruridae family. In this work, we describe for the first time the complete mitogenome (mitochondrial genome) of M. capsulata, being the first described of the family Pneumospiruridae. The mitogenome of M. capsulata has 13,659bp in length, an A + T content of 79.2%. The mitogenome included 12 protein-coding genes (PCGs) (lacking the atp8 gene), 22 tRNA genes, 2 rRNA genes (all the genes are coded by the heavy strand), and an AT-rich region. The PCGs varied in size (232bp-1645bp). Only the tRNA-Trp has the standard cloverleaf secondary structure, while the other 21 do not. The AT-rich region, with a 90.5% A + T content and a length of 389bp, is located between the cox3 and tRNA-Ala genes. Comparison with the mitogenomes of 29 species of Spiruromorpha infraorder, belonging to different families, demonstrates that M. capsulata mitogenome shared the common characteristics of most of them. The phylogeny constructions yielded phylogenies that were in agreement with the obtained previously by using sequences and gene order data of mitogenomes.

Different gene rearrangements of the genus Dardanus (Anomura: Diogenidae) and insights into the phylogeny of Paguroidea

Complete mitochondrial genomes (mitogenomes) can provide useful information for phylogenetic relationships, gene rearrangement, and molecular evolution. In this study, the complete mitogenomes of two hermit crabs, Dardanus arrosor and Dardanus aspersus, were sequenced for the first time and compared with other published mitogenomes of Paguroidea. Each of the two mitogenomes contains an entire set of 37 genes and a putative control region, but they display different gene arrangements. The different arrangements of the two mitogenomes might be the result of transposition, reversal, and tandem duplication/random loss events from the ancestral pancrustacean pattern. Genome sequence similarity analysis reveals the gene rearrangement in 15 Paguroidea mitogenomes. After synteny analysis between the 15 Paguroidea mitogenomes, an obvious rearranged region is found in D. aspersus mitogenome. Across the 13 protein-coding genes (PCGs) tested, COI has the least and ND6 has the largest genetic distances among the 15 hermit crabs, indicating varied evolution rates of PCGs. In addition, the dN/dS ratio analysis shows that all PCGs are evolving under purifying selection. The phylogenetic analyses based on both gene order and sequence data present the monophyly of three families (Paguridae, Coenobitidae, and Pylochelidae) and the paraphyly of the family Diogenidae. Meanwhile, the phylogenetic tree based on the nucleotide sequences of 13 PCGs shows that two Dardanus species formed a sister group with five Coenobitidae species. These findings help to better understand the gene rearrangement and phylogeny of Paguroidea, as well as provide new insights into the usefulness of mitochondrial gene order as a phylogenetic marker.

Structure, gene order, and nucleotide composition of mitochondrial genomes in parasitic lice from Amblycera.

Parasitic lice have unique mitochondrial (mt) genomes characterized by rearranged gene orders, variable genome structures, and less AT content compared to most other insects. However, relatively little is known about the mt genomes of Amblycera, the suborder sister to all other parasitic lice. Comparing among nine different genera (including representative of all seven families), we show that Amblycera have variable and highly rearranged mt genomes. Some genera have fragmented genomes that vary considerably in length, whereas others have a single mt chromosome. Notably, these genomes are more AT-biased than most other lice. We also recover genus-level phylogenetic relationships among Amblycera that are consistent with those reported from large nuclear datasets, indicating that mt sequences are reliable for reconstructing evolutionary relationships in Amblycera. However, gene order data cannot reliably recover these same relationships. Overall, our results suggest that the mt genomes of lice, already know to be distinctive, are even more variable than previously thought.

Revision shock in Pacific oysters taxonomy: the genusMagallana(formerlyCrassostreain part) is well-founded and necessary

AbstractThe description of the genus Magallana provoked taxonomic ‘revision shock’ among Bayne and colleagues, who criticised: (1) the supporting evidence for this taxonomic change, (2) the procedure used and (3) its impact on the nomenclatural stability of oysters. Here, we demonstrate that the description of the genus Magallana fulfils the norms of taxonomic and indexing revisions, and is well-founded on a scientific basis. The clade named Magallana is supported by a robust phylogeny based on comprehensive taxon sampling, independent datasets and varied analytical methods. A new maximum likelihood phylogenetic analysis of oysters, based on mitochondrial gene order data (representing molecular phenotypes above the sequence level), has provided additional support for this classification. These results are fully consistent with all previously published phylogenetic studies, thus providing an unambiguous indication of the stability of the clade that meets all the currently accepted criteria for naming clades as taxa. We show that, while the criterion of ‘morphological diagnosability’ is not applicable at any level of oyster classification, several striking molecular phenotypes are diagnostic of Magallana, both at the nuclear and the mitochondrial genome level. The classification with three genera Magallana, Crassostrea and Talonostrea reflects the evolutionary diversity of Crassostreinae and is consistent with taxonomic ranking criteria adopted for other oyster subfamilies.

Phylogenetic implications of mitogenome rearrangements in East Asian potamiscine freshwater crabs (Brachyura: Potamidae)

We present here a combined mitogenome gene order dataset totalling 62% of the known genera of East Asian potamiscine freshwater crabs that includes first-time mitogenome data for 14 species and gene order data from 37 other species of potamiscines. A surprisingly high number of mitogenome gene order rearrangements were found in the taxa studied (comprising nine different rearrangements and seven major patterns, one of which has two sub-arrangements). Our phylogenetic reconstructions indicate that the mitogenome gene order rearrangements are associated with the evolutionary history of potamiscine lineages. We also used a new Event-based Maximum Parsimony method to reconstruct ancestral gene orders, which takes into consideration gene duplication, pseudogeneticization, and tandem duplication random loss. Furthermore, shared mitogenome gene order patterns were used to inform the taxonomic placement of Sinopotamon parvum, and the cryptic diversity in Potamiscus. The remarkably frequent mitogenome rearrangements in potamiscine freshwater crabs have great potential to contribute to our understanding of the evolutionary history of these highly diverse decapods in East Asia.

The complete mitochondrial genomes of two vent squat lobsters, Munidopsis lauensis and M.verrilli: Novel gene arrangements and phylogenetic implications.

Hydrothermal vents are considered as one of the most extremely harsh environments on the Earth. In this study, the complete mitogenomes of hydrothermal vent squat lobsters, Munidopsis lauensis and M. verrilli, were determined through Illumina sequencing and compared with other available mitogenomes of anomurans. The mitogenomes of M. lauensis (17,483 bp) and M. verrilli (17,636 bp) are the largest among all Anomura mitogenomes, while the A+T contents of M. lauensis (62.40%) and M. verrilli (63.99%) are the lowest. The mitogenomes of M. lauensis and M. verrilli display novel gene arrangements, which might be the result of three tandem duplication–random loss (tdrl) events from the ancestral pancrustacean pattern. The mitochondrial gene orders of M. lauensis and M. verrilli shared the most similarities with S. crosnieri. The phylogenetic analyses based on both gene order data and nucleotide sequences (PCGs and rRNAs) revealed that the two species were closely related to Shinkaia crosnieri. Positive selection analysis revealed that eighteen residues in seven genes (atp8, Cytb, nad3, nad4, nad4l, nad5, and nad6) of the hydrothermal vent anomurans were positively selected sites.

Order, please! Uncertainty in the ordinal-level classification of Chlorophyceae.

BackgroundChlorophyceae is one of three most species-rich green algal classes and also the only class in core Chlorophyta whose monophyly remains uncontested as gene and taxon sampling improves. However, some key relationships within Chlorophyceae are less clear-cut and warrant further investigation. The present study combined genome-scale chloroplast data and rich sampling in an attempt to resolve the ordinal classification in Chlorophyceae. The traditional division into Sphaeropleales and Volvocales (SV), and a clade containing Oedogoniales, Chaetopeltidales, and Chaetophorales (OCC) was of particular interest with the addition of deeply branching members of these groups, as well as the placement of several incertae sedis taxa.MethodsWe sequenced 18 chloroplast genomes across Chlorophyceae to compile a data set of 58 protein-coding genes of a total of 68 chlorophycean taxa. We analyzed the concatenated nucleotide and amino acid datasets in the Bayesian and Maximum Likelihood frameworks, supplemented by analyses to examine potential discordant signal among genes. We also examined gene presence and absence data across Chlorophyceae.ResultsConcatenated analyses yielded at least two well-supported phylogenies: nucleotide data supported the traditional classification with the inclusion of the enigmatic Treubarinia into Sphaeropleales sensu lato. However, amino acid data yielded equally strong support for Sphaeropleaceae as sister to Volvocales, with the rest of the taxa traditionally classified in Sphaeropleales in a separate clade, and Treubarinia as sister to all of the above. Single-gene and other supplementary analyses indicated that the data have low phylogenetic signal at these critical nodes. Major clades were supported by genomic structural features such as gene losses and trans-spliced intron insertions in the plastome.DiscussionWhile the sequence and gene order data support the deep split between the SV and OCC lineages, multiple phylogenetic hypotheses are possible for Sphaeropleales s.l. Given this uncertainty as well as the higher-taxonomic disorder seen in other algal groups, dwelling on well-defined, strongly supported Linnaean orders is not currently practical in Chlorophyceae and a less formal clade system may be more useful in the foreseeable future. For example, we identify two strongly and unequivocally supported clades: Treubarinia and Scenedesminia, as well as other smaller groups that could serve a practical purpose as named clades. This system does not preclude future establishment of new orders, or emendment of the current ordinal classification if new data support such conclusions.

A Median Solver and Phylogenetic Inference Based on Double-Cut-and-Join Sorting.

Genome rearrangement is known as one of the main evolutionary mechanisms on the genomic level. Phylogenetic analysis based on rearrangement played a crucial role in biological research in the past decades, especially with the increasing availability of fully sequenced genomes. In general, phylogenetic analysis aims to solve two problems: small parsimony problem (SPP) and big parsimony problem (BPP). Maximum parsimony is a popular approach for SPP and BPP, which relies on iteratively solving an NP-hard problem, the median problem. As a result, current median solvers and phylogenetic inference methods based on the median problem all face serious problems on scalability and cannot be applied to data sets with large and distant genomes. In this article, we propose a new median solver for gene order data that combines double-cut-and-join sorting with the simulated annealing algorithm. Based on this median solver, we built a new phylogenetic inference method to solve both SPP and BPP problems. Our experimental results show that the new median solver achieves an excellent performance on simulated data sets, and the phylogenetic inference tool built based on the new median solver has a better performance than other existing methods.

Phylogeny analysis from gene-order data with massive duplications

BackgroundGene order changes, under rearrangements, insertions, deletions and duplications, have been used as a new type of data source for phylogenetic reconstruction. Because these changes are rare compared to sequence mutations, they allow the inference of phylogeny further back in evolutionary time. There exist many computational methods for the reconstruction of gene-order phylogenies, including widely used maximum parsimonious methods and maximum likelihood methods. However, both methods face challenges in handling large genomes with many duplicated genes, especially in the presence of whole genome duplication.MethodsIn this paper, we present three simple yet powerful methods based on maximum-likelihood (ML) approaches that encode multiplicities of both gene adjacency and gene content information for phylogenetic reconstruction.ResultsExtensive experiments on simulated data sets show that our new method achieves the most accurate phylogenies compared to existing approaches. We also evaluate our method on real whole-genome data from eleven mammals. The package is publicly accessible at http://www.geneorder.org.ConclusionsOur new encoding schemes successfully incorporate the multiplicity information of gene adjacencies and gene content into an ML framework, and show promising results in reconstruct phylogenies for whole-genome data in the presence of massive duplications.

Comparative analyses of chloroplast genome data representing nine green algae in Sphaeropleales (Chlorophyceae, Chlorophyta).

The chloroplast genomes of green algae are highly variable in their architecture. In this article we summarize gene content across newly obtained and published chloroplast genomes in Chlorophyceae, including new data from nine of species in Sphaeropleales (Chlorophyceae, Chlorophyta). We present genome architecture information, including genome synteny analysis across two groups of species. Also, we provide a phylogenetic tree obtained from analysis of gene order data for species in Chlorophyceae with fully sequenced chloroplast genomes. Further analyses and interpretation of the data can be found in “Chloroplast phylogenomic data from the green algal order Sphaeropleales (Chlorophyceae, Chlorophyta) reveal complex patterns of sequence evolution” (Fučíková et al., In review) [1].

Distinctive Architecture of the Chloroplast Genome in the Chlorodendrophycean Green Algae Scherffelia dubia and Tetraselmis sp. CCMP 881.

The Chlorodendrophyceae is a small class of green algae belonging to the core Chlorophyta, an assemblage that also comprises the Pedinophyceae, Trebouxiophyceae, Ulvophyceae and Chlorophyceae. Here we describe for the first time the chloroplast genomes of chlorodendrophycean algae (Scherffelia dubia, 137,161 bp; Tetraselmis sp. CCMP 881, 100,264 bp). Characterized by a very small single-copy (SSC) region devoid of any gene and an unusually large inverted repeat (IR), the quadripartite structures of the Scherffelia and Tetraselmis genomes are unique among all core chlorophytes examined thus far. The lack of genes in the SSC region is offset by the rich and atypical gene complement of the IR, which includes genes from the SSC and large single-copy regions of prasinophyte and streptophyte chloroplast genomes having retained an ancestral quadripartite structure. Remarkably, seven of the atypical IR-encoded genes have also been observed in the IRs of pedinophycean and trebouxiophycean chloroplast genomes, suggesting that they were already present in the IR of the common ancestor of all core chlorophytes. Considering that the relationships among the main lineages of the core Chlorophyta are still unresolved, we evaluated the impact of including the Chlorodendrophyceae in chloroplast phylogenomic analyses. The trees we inferred using data sets of 79 and 108 genes from 71 chlorophytes indicate that the Chlorodendrophyceae is a deep-diverging lineage of the core Chlorophyta, although the placement of this class relative to the Pedinophyceae remains ambiguous. Interestingly, some of our phylogenomic trees together with our comparative analysis of gene order data support the monophyly of the Trebouxiophyceae, thus offering further evidence that the previously observed affiliation between the Chlorellales and Pedinophyceae is the result of systematic errors in phylogenetic reconstruction.

MLGO: phylogeny reconstruction and ancestral inference from gene-order data.

BackgroundThe rapid accumulation of whole-genome data has renewed interest in the study of using gene-order data for phylogenetic analyses and ancestral reconstruction. Current software and web servers typically do not support duplication and loss events along with rearrangements.ResultsMLGO (Maximum Likelihood for Gene-Order Analysis) is a web tool for the reconstruction of phylogeny and/or ancestral genomes from gene-order data. MLGO is based on likelihood computation and shows advantages over existing methods in terms of accuracy, scalability and flexibility.ConclusionsTo the best of our knowledge, it is the first web tool for analysis of large-scale genomic changes including not only rearrangements but also gene insertions, deletions and duplications. The web tool is available from http://www.geneorder.org/server.php.

Molecular identification key of the family Streptococcaceae

The gene order conservation (GOC) between the species of family Streptococcaceae was analysed. The rate of GOC in the strains belonging to the same species is 70% or more. When we compared different species belonging to the same genus, the rate of GOC was 30-47% while it was below 20% when the species were from different genera. A molecular identification key was established for identifying those genera and species within the family Streptococcaceae which have an already known full genome sequence (24 Streptococcus and 2 Lactococcus species). Identical genome parts of the species belonging to the same genus were used for determination of genera. These are the sections surrounding the replication origin dnaA, the sequence from gene phaB to the gene accA, and the sequence of alr acpS secA. Sections around the genes pepX, leuS and rplM were used for identifying the species. The gene order analysis and data in molecular identification key showed that S. uberis and S. parauberis also belong to the same species, and our suggestion for their new names is S. uberis subsp. uberis and S uberis subsp. parauberis. Based on this data, a new definition of bacterial species is proposed: two isolates belong to the same species if the order of the genes in their genomes is almost identical.

Mitochondrial genomes to the rescue – Diurodrilidae in the myzostomid trap

Diurodrilidae is a taxon of Lophotrochozoa comprising about six, exclusively interstitial species, which are up to 500μm long and dorsoventrally flattened. Traditionally, Diurodrilidae had been regarded as an annelid family. However, recently Diurodrilidae had been excluded from Annelida and been placed in closer relationship to platyzoan taxa based on both morphological and nuclear rRNA data. Since both, Diurodrilidae and platyzoan taxa, exhibit long branches in the molecular analyses, the close relationship might be due to a long branch attraction artifact. The annelid taxon Myzostomida had been trapped in a similar long branch attraction artifact with platyzoan taxa using nuclear rRNA data, but determination of the nearly complete mitochondrial genome of myzostomids revealed their annelid affinity. Therefore, we determined the nearly complete mitochondrial genome of Diurodrilus subterraneus as well as new nuclear rRNA data for D. subterraneus and some platyzoan taxa. All our analyses of nuclear rRNA and mitochondrial sequence and gene order data presented herein clearly place Diurodrilidae within Annelida and with strong nodal support values in some analyses. Therefore, the previously suggested exclusion of Diurodrilidae from Annelida and its close relationship with platyzoan taxa can be attributed to a long branch artifact. Morphological data do not unambiguously support a platyzoan affinity of Diurodrilidae, but instead would also be in line with a progenetic origin of Diurodrilidae within Annelida.

Inferring Phylogenetic Networks from Gene Order Data

Existing algorithms allow us to infer phylogenetic networks from sequences (DNA, protein or binary), sets of trees, and distance matrices, but there are no methods to build them using the gene order data as an input. Here we describe several methods to build split networks from the gene order data, perform simulation studies, and use our methods for analyzing and interpreting different real gene order datasets. All proposed methods are based on intermediate data, which can be generated from genome structures under study and used as an input for network construction algorithms. Three intermediates are used: set of jackknife trees, distance matrix, and binary encoding. According to simulations and case studies, the best intermediates are jackknife trees and distance matrix (when used with Neighbor-Net algorithm). Binary encoding can also be useful, but only when the methods mentioned above cannot be used.

A Cluster-on-a-Chip Architecture for High-Throughput Phylogeny Search

In this paper, we describe an FPGA-based coprocessor architecture that performs a high-throughput branch-and-bound search of the space of phylogenetic trees corresponding to the number of input taxa. Our coprocessor architecture is designed to accelerate maximum-parsimony phylogeny reconstruction for gene-order and sequence data and is amenable to both exhaustive and heuristic tree searches. Our architecture exposes coarse-grain parallelism by dividing the search space among parallel processing elements (PEs) and each PE exposes fine-grain memory parallelism for their lower-bound computation, the kernel computation performed by each PE. Inter-PE communication is performed entirely on-chip. When using this coprocessor for maximum-parsimony reconstruction for gene-order data, our coprocessor achieves a 40X improvement over software in search throughput, corresponding to a 14X end-to-end application improvement when including all communication and systems overheads.

The complete mitochondrial genomes of six heterodont bivalves (Tellinoidea and Solenoidea): variable gene arrangements and phylogenetic implications.

BackgroundTaxonomy and phylogeny of subclass Heterodonta including Tellinoidea are long-debated issues and a complete agreement has not been reached yet. Mitochondrial (mt) genomes have been proved to be a powerful tool in resolving phylogenetic relationship. However, to date, only ten complete mitochondrial genomes of Heterodonta, which is by far the most diverse major group of Bivalvia, have been determined. In this paper, we newly sequenced the complete mt genomes of six species belonging to Heterodonta in order to resolve some problematical relationships among this subclass.Principal FindingsThe complete mt genomes of six species vary in size from 16,352 bp to 18,182. Hairpin-like secondary structures are found in the largest non-coding regions of six freshly sequenced mt genomes, five of which contain tandem repeats. It is noteworthy that two species belonging to the same genus show different gene arrangements with three translocations. The phylogenetic analysis of Heterodonta indicates that Sinonovacula constricta, distant from the Solecurtidae belonging to Tellinoidea, is as a sister group with Solen grandis of family Solenidae. Besides, all five species of Tellinoidea cluster together, while Sanguinolaria diphos has closer relationship with Solecurtus divaricatus, Moerella iridescens and Semele scaba rather than with Sanguinolaria olivacea.Conclusions/SignificanceBy comparative study of gene order rearrangements and phylogenetic relationships of the five species belonging to Tellinoidea, our results support that comparisons of mt gene order rearrangements, to some extent, are a useful tool for phylogenetic studies. Based on phylogenetic analyses of multiple protein-coding genes, we prefer classifying the genus Sinonovacula within the superfamily Solenoidea and not the superfamily Tellinoidea. Besides, both gene order and sequence data agree that Sanguinolaria (Psammobiidae) is not monophyletic. Nevertheless, more studies based on more mt genomes via combination of gene order and phylogenetic analysis are needed to further understand the phylogenetic relationships in subclass Heterodonta.

MitoZoa 2.0: a database resource and search tools for comparative and evolutionary analyses of mitochondrial genomes in Metazoa

The MITOchondrial genome database of metaZOAns (MitoZoa) is a public resource for comparative analyses of metazoan mitochondrial genomes (mtDNA) at both the sequence and genomic organizational levels. The main characteristics of the MitoZoa database are the careful revision of mtDNA entry annotations and the possibility of retrieving gene order and non-coding region (NCR) data in appropriate formats. The MitoZoa retrieval system enables basic and complex queries at various taxonomic levels using different search menus. MitoZoa 2.0 has been enhanced in several aspects, including: a re-annotation pipeline to check the correctness of protein-coding gene predictions; a standardized annotation of introns and of precursor ORFs whose functionality is post-transcriptionally recovered by RNA editing or programmed translational frameshifting; updates of taxon-related fields and a BLAST sequence similarity search tool. Database novelties and the definition of standard mtDNA annotation rules, together with the user-friendly retrieval system and the BLAST service, make MitoZoa a valuable resource for comparative and evolutionary analyses as well as a reference database to assist in the annotation of novel mtDNA sequences. MitoZoa is freely accessible at http://www.caspur.it/mitozoa.

Rec-DCM-Eigen: reconstructing a less parsimonious but more accurate tree in shorter time.

Maximum parsimony (MP) methods aim to reconstruct the phylogeny of extant species by finding the most parsimonious evolutionary scenario using the species' genome data. MP methods are considered to be accurate, but they are also computationally expensive especially for a large number of species. Several disk-covering methods (DCMs), which decompose the input species to multiple overlapping subgroups (or disks), have been proposed to solve the problem in a divide-and-conquer way.We design a new DCM based on the spectral method and also develop the COGNAC (Comparing Orders of Genes using Novel Algorithms and high-performance Computers) software package. COGNAC uses the new DCM to reduce the phylogenetic tree search space and selects an output tree from the reduced search space based on the MP principle. We test the new DCM using gene order data and inversion distance. The new DCM not only reduces the number of candidate tree topologies but also excludes erroneous tree topologies which can be selected by original MP methods. Initial labeling of internal genomes affects the accuracy of MP methods using gene order data, and the new DCM enables more accurate initial labeling as well. COGNAC demonstrates superior accuracy as a consequence. We compare COGNAC with FastME and the combination of the state of the art DCM (Rec-I-DCM3) and GRAPPA . COGNAC clearly outperforms FastME in accuracy. COGNAC –using the new DCM–also reconstructs a much more accurate tree in significantly shorter time than GRAPPA with Rec-I-DCM3.

The complete mitochondrial genome sequence of Whitmania pigra (Annelida, Hirudinea): The first representative from the class Hirudinea

The mitochondrial genome is a significant tool for investigating the evolutionary history of metazoan animals. The currently available mitochondrial genome data in GenBank is limited to understand the detail evolutionary relationship among the metazoan animals, especially in the phylum Annelida. Here we present the mitochondrial gene organization, gene order and codon usage of the leech Whitmania pigra (Annelida), which is the first representative from the class Hirudinea. It is a circular molecule of 14,426 bp, and encodes 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. All 37 genes of W. pigra mitochondrial genome are transcribed from the same strand, which is identical to studied annelids, two echiurans, two sipunculans and many other lophotrochozoans. Five conserved gene clusters can be found in mitochondrial genomes of nine studied annelids, including (1) cox1- N- cox2; (2) cox3- Q- nad6- cob- W- atp6; (3) H- nad5- F- E- P- T- nad4L- nad4; (4) srRNA- V- lrRNA; and (5) nad3- S 1 - nad2. Compared with that of other studied annelids, translocations of transfer RNAs were found in the gene arrangement of W. pigra mitochondrial genome. Phylogenetic analysis strongly support that the species from Hirudinina and Oligochaeta form a monophyletic group Clitellata (BPM = 100, BPP = 100), which is consistent with previous research based on morphological and other molecular data. Both gene order data and amino acid sequences reveal that echiurans are derived annelids and sipunculans should be clustered with annelids and echiurans.