On the phylogenetic position of Myzostomida: can 77 genes get it wrong?

First report of the mitogenome of the invasive reef-building polychaete Ficopomatus enigmaticus (Annelida: Serpulidae) and a cryptic lineage from the Japanese Archipelago

BACKGROUND

The mitochondrial genomes (mitogenomes) of the family Serpulidae are characterized by a high nucleotide sequence divergence and a significant number of gene order rearrangements compared with other families within the phylum Annelida. However, only two of 50 genera of serpulids have mitogenomes already sequenced. In this study, we report the first sequencing and assembly of the complete mitogenome of Ficopomatus, thus providing further knowledge on mitochondrial gene sequences of Serpulidae.

METHODS AND RESULTS

A mitogenome of the invasive reef-building polychaete Ficopomatus enigmaticus was amplified by long PCR and sequenced using the Illumina MiSeq System. It comprised 15,853 bp and consisted of 12 protein-coding genes (atp8 was not found), 23 tRNA, and two rRNA genes. The AT and GC skew values were infrequent when compared to annelid mitogenomes but similar to other serpulids sequenced to date (i.e., Spirobranchus and Hydroides). The mitochondrial gene order of F. enigmaticus was highly rearranged compared to other serpulids. To amplify 16S rRNA gene sequences, we developed a 16S rRNA primer set by modifying the universal primer set 16SarL/16SbrH. We detected the 16S rRNA sequence of F. enigmaticus deposited in GenBank erroneously characterized as of serpulid origin. We reported for the first time the presence of two lineages of F. enigmaticus in Japan, which have already been identified in California, Australia, and the Mediterranean.

CONCLUSIONS

The first mitochondrial genome of F. enigmaticus showed a unique gene order rearrangement, corroborating the remarkable diversity in the previously reported mitogenomes of other serpulid species. The presence of the two lineages of F. enigmaticus identified for the first time in Japan represents another case of cryptic invasion. The first 16S rRNA gene sequences of F. enigmaticus obtained in the present study can be used as reference sequences in future DNA metabarcoding studies.

Comparative mitogenomes provide new insights into phylogeny and taxonomy of the subfamily Xenocyprinae (Cypriniformes: Cyprinidae)

Xenocyprinae is a cyprinid subfamily that not only has a discrete geographic distribution but also has a long history dating to the Early Miocene. However, it is controversial whether systematic classification and some species validity of Xenocyprinae exist, as well as its phylogenetic relationships and evolutionary history. In the present study, we first reviewed the description and taxonomic history of Xenocyprinae, and then the complete mitochondrial genome of Distoechodon compressus, an endemic and locally distributed species belonging to Xenocyprinae, was sequenced and annotated. Finally, all the mitogenomes of Xenocyprinae were compared to reconstruct the phylogenetic relationship and estimate the divergence time. The results showed that the mitogenomes are similar in organization and structure with 16618–16630 bp length from 12 mitogenomes of eight species. Phylogenetic analysis confirmed the monology of Xenocyprinae and illustrated three clades within the Xenocyprinae to consist of ambiguous generic classification. Plagiognathops is a valid genus located at the base of the phylogenetic tree. The genus Xenocypris was originally monophyletic, but X. fangi was excluded. Divergence time estimation revealed that the earliest divergence within Xenocyprinae occurred approximately 12.1 Mya when Plagiognathops separated from the primitive Xenocypris. The main two clades (Xenocypris and (Distoechodon + Pseudobrama + X. fangi)) diverged 10.0 Mya. The major divergence of Xenocyprinae species possibly occurred in the Middle to Late Miocene and Late Pliocene, suggesting that speciation and diversifications could be attributed to the Asian monsoon climate. This study clarifies some controversial issues of systematics and provides essential information on the taxonomy and phylogeny of the subfamily Xenocyprinae.

Mitogenome of a stink worm (Annelida: Travisiidae) includes degenerate group II intron that is also found in five congeneric species

Mitogenomes are useful for inferring phylogenetic relationships between organisms. Although the mitogenomes of Annelida, one of the most morphologically and ecologically diverse metazoan groups have been well sequenced, those of several families remain unexamined. This study determined the first mitogenome from the family Travisiidae (Travisia sanrikuensis), analyzed its mitogenomic features, and reconstructed a phylogeny of Sedentaria. The monophyly of the Terebellida + Arenicolida + Travisiidae clade is supported by molecular phylogenetic analysis. The placement of Travisiidae is unclear because of the lack of mitogenomes from closely related lineages. An unexpected intron appeared within the cox1 gene of T. sanrikuensis and in the same positions of five undescribed Travisia spp. Although the introns are shorter (790–1386 bp) than other group II introns, they can be considered degenerate group II introns due to type II intron maturase open reading frames, found in two of the examined species, and motifs characteristic of group II introns. This is likely the first known case in metazoans where mitochondrial group II introns obtained by a common ancestor are conserved in several descendants. Insufficient evolutionary time for intron loss in Travisiidae, or undetermined mechanisms may have helped maintain the degenerate introns.

Annelid Diversity: Historical Overview and Future Perspectives

Annelida is a ubiquitous, common and diverse group of organisms, found in terrestrial, fresh waters and marine environments. Despite the large efforts put into resolving the evolutionary relationships of these and other Lophotrochozoa, and the delineation of the basal nodes within the group, these are still unanswered. Annelida holds an enormous diversity of forms and biological strategies alongside a large number of species, following Arthropoda, Mollusca, Vertebrata and perhaps Platyhelminthes, among the species most rich in phyla within Metazoa. The number of currently accepted annelid species changes rapidly when taxonomic groups are revised due to synonymies and descriptions of a new species. The group is also experiencing a recent increase in species numbers as a consequence of the use of molecular taxonomy methods, which allows the delineation of the entities within species complexes. This review aims at succinctly reviewing the state-of-the-art of annelid diversity and summarizing the main systematic revisions carried out in the group. Moreover, it should be considered as the introduction to the papers that form this Special Issue on Systematics and Biodiversity of Annelids.

The complete mitochondrial genome of the sand bubbler crab Scopimera globosa and its phylogenetic position

The mitochondrial genome has become commonly used for the molecular phylogenetic analysis of animals. Most phylogenetic studies on brachyurans using mitogenome sequences have indicated the paraphyly of superfamilies Grapsoidea and Ocypodoidea but taxon sampling remains limited. The phylogenetic position of Scopimera has been tested in several previous studies using nuclear and/or mitochondrial DNA sequences, but the phylogenetic relationship within the family remains to be resolved. We newly sequenced the complete mitochondrial genome of the sand bubbler crab Scopimera globosa (Ocypodoidea: Dotillidae). Scopimera globosa was clustered with Ilyoplax despite the morphological similarity between Scopimera and Dotilla. The mitochondrial gene order of S. globosa was unique, whereas that of other genera in the family was the same. These results suggest that phylogenetic analysis based on mitogenome sequences and gene order comparison would provide a more robust phylogeny of Dotillidae.

Phylogenomics of Lophotrochozoa with Consideration of Systematic Error

Phylogenomic studies have improved understanding of deep metazoan phylogeny and show promise for resolving incongruences among analyses based on limited numbers of loci. One region of the animal tree that has been especially difficult to resolve, even with phylogenomic approaches, is relationships within Lophotrochozoa (the animal clade that includes molluscs, annelids, and flatworms among others). Lack of resolution in phylogenomic analyses could be due to insufficient phylogenetic signal, limitations in taxon and/or gene sampling, or systematic error. Here, we investigated why lophotrochozoan phylogeny has been such a difficult question to answer by identifying and reducing sources of systematic error. We supplemented existing data with 32 new transcriptomes spanning the diversity of Lophotrochozoa and constructed a new set of Lophotrochozoa-specific core orthologs. Of these, 638 orthologous groups (OGs) passed strict screening for paralogy using a tree-based approach. In order to reduce possible sources of systematic error, we calculated branch-length heterogeneity, evolutionary rate, percent missing data, compositional bias, and saturation for each OG and analyzed increasingly stricter subsets of only the most stringent (best) OGs for these five variables. Principal component analysis of the values for each factor examined for each OG revealed that compositional heterogeneity and average patristic distance contributed most to the variance observed along the first principal component while branch-length heterogeneity and, to a lesser extent, saturation contributed most to the variance observed along the second. Missing data did not strongly contribute to either. Additional sensitivity analyses examined effects of removing taxa with heterogeneous branch lengths, large amounts of missing data, and compositional heterogeneity. Although our analyses do not unambiguously resolve lophotrochozoan phylogeny, we advance the field by reducing the list of viable hypotheses. Moreover, our systematic approach for dissection of phylogenomic data can be applied to explore sources of incongruence and poor support in any phylogenomic data set. [Annelida; Brachiopoda; Bryozoa; Entoprocta; Mollusca; Nemertea; Phoronida; Platyzoa; Polyzoa; Spiralia; Trochozoa.].

Curious bivalves: Systematic utility and unusual properties of anomalodesmatan mitochondrial genomes

Mitogenomic trees for Bivalvia have proved problematic in the past, but several highly divergent lineages were missing from these analyses and increased representation of these groups may yet improve resolution. Here, we add seven new sequences from the Anomalodesmata and one unidentified semelid species (Bryopa lata, Euciroa cf. queenslandica, Laternula elliptica, Laternula truncata, Lyonsia norwegica, Myadora brevis, Tropidomya abbreviata, "Abra" sp.). We show that relationships in a mitogenomic tree for the Class are improved by the addition of seven anomalodesmatans from this highly divergent clade, but are still not completely consistent with relationships recovered in studies of nuclear genes. We suggest that some anomalous relationships (for instance the non-monophyly of Bivalvia) may be partially explained by compositional heterogeneity in the mitogenome and suggest that the addition of more taxa may help resolve both this effect and possible instances of long branch attraction. We also identify several curious features about anomalodesmatan mitogenomes. For example, many protein-coding gene boundaries are poorly defined in marine bivalves, but particularly so in anomalodesmatans, primarily due to non-conserved boundary sequences. The use of transcriptomic and genomic data together enabled better definition of gene boundaries, the identification of possible pseudogenes and suggests that most genes are translated monocistronically, which contrasts with many other studies. We also identified a possible case of gene duplication of ND5 in Myadora brevis (Myochamidae). Mitogenome size in the Anomalodesmata ranges from very small compact molecules, with the smallest for Laternula elliptica (Laternulidae) only 14,622bp, to Bryopa lata (Clavagellidae) which is at least 31,969bp long and may be >40,000bp. Finally, sampled species show a high degree of sequence divergence and variable gene order, although intraspecific variation in Laternula elliptica is very low.

Multiple measures could alleviate long-branch attraction in phylogenomic reconstruction of Cupressoideae (Cupressaceae)

Long-branch attraction (LBA) is a major obstacle in phylogenetic reconstruction. The phylogenetic relationships among Juniperus (J), Cupressus (C) and the Hesperocyparis-Callitropsis-Xanthocyparis (HCX) subclades of Cupressoideae are controversial. Our initial analyses of plastid protein-coding gene matrix revealed both J and C with much longer stem branches than those of HCX, so their sister relationships may be attributed to LBA. We used multiple measures including data filtering and modifying, evolutionary model selection and coalescent phylogenetic reconstruction to alleviate the LBA artifact. Data filtering by strictly removing unreliable aligned regions and removing substitution saturation genes and rapidly evolving sites could significantly reduce branch lengths of subclades J and C and recovered a relationship of J (C, HCX). In addition, using coalescent phylogenetic reconstruction could elucidate the LBA artifact and recovered J (C, HCX). However, some valid methods for other taxa were inefficient in alleviating the LBA artifact in J-C-HCX. Different strategies should be carefully considered and justified to reduce LBA in phylogenetic reconstruction of different groups. Three subclades of J-C-HCX were estimated to have experienced ancient rapid divergence within a short period, which could be another major obstacle in resolving relationships. Furthermore, our plastid phylogenomic analyses fully resolved the intergeneric relationships of Cupressoideae.

Comparative Mitogenomics of Leeches (Annelida: Clitellata): Genome Conservation and Placobdella-Specific trnD Gene Duplication

Mitochondrial DNA sequences, often in combination with nuclear markers and morphological data, are frequently used to unravel the phylogenetic relationships, population dynamics and biogeographic histories of a plethora of organisms. The information provided by examining complete mitochondrial genomes also enables investigation of other evolutionary events such as gene rearrangements, gene duplication and gene loss. Despite efforts to generate information to represent most of the currently recognized groups, some taxa are underrepresented in mitochondrial genomic databases. One such group is leeches (Annelida: Hirudinea: Clitellata). Herein, we expand our knowledge concerning leech mitochondrial makeup including gene arrangement, gene duplication and the evolution of mitochondrial genomes by adding newly sequenced mitochondrial genomes for three bloodfeeding species: Haementeria officinalis, Placobdella lamothei and Placobdella parasitica. With the inclusion of three new mitochondrial genomes of leeches, a better understanding of evolution for this organelle within the group is emerging. We found that gene order and genomic arrangement in the three new mitochondrial genomes is identical to previously sequenced members of Clitellata. Interestingly, within Placobdella, we recovered a genus-specific duplication of the trnD gene located between cox2 and atp8. We performed phylogenetic analyses using 12 protein-coding genes and expanded our taxon sampling by including GenBank sequences for 39 taxa; the analyses confirm the monophyletic status of Clitellata, yet disagree in several respects with other phylogenetic hypotheses based on morphology and analyses of non-mitochondrial data.

The Utility of Genome Skimming for Phylogenomic Analyses as Demonstrated for Glycerid Relationships (Annelida, Glyceridae)

Glyceridae (Annelida) are a group of venomous annelids distributed worldwide from intertidal to abyssal depths. To trace the evolutionary history and complexity of glycerid venom cocktails, a solid backbone phylogeny of this group is essential. We therefore aimed to reconstruct the phylogenetic relationships of these annelids using Illumina sequencing technology. We constructed whole-genome shotgun libraries for 19 glycerid specimens and 1 outgroup species (Glycinde armigera). The chosen target genes comprise 13 mitochondrial proteins, 2 ribosomal mitochondrial genes, and 4 nuclear loci (18SrRNA, 28SrRNA, ITS1, and ITS2). Based on partitioned maximum likelihood as well as Bayesian analyses of the resulting supermatrix, we were finally able to resolve a robust glycerid phylogeny and identified three clades comprising the majority of taxa. Furthermore, we detected group II introns inside the cox1 gene of two analyzed glycerid specimens, with two different insertions in one of these species. Moreover, we generated reduced data sets comprising 10 million, 4 million, and 1 million reads from the original data sets to test the influence of the sequencing depth on assembling complete mitochondrial genomes from low coverage genome data. We estimated the coverage of mitochondrial genome sequences in each data set size by mapping the filtered Illumina reads against the respective mitochondrial contigs. By comparing the contig coverage calculated in all data set sizes, we got a hint for the scalability of our genome skimming approach. This allows estimating more precisely the number of reads that are at least necessary to reconstruct complete mitochondrial genomes in Glyceridae and probably non-model organisms in general.

GPU MrBayes V3.1: MrBayes on Graphics Processing Units for Protein Sequence Data

We present a modified GPU (graphics processing unit) version of MrBayes, called ta(MC)(3) (GPU MrBayes V3.1), for Bayesian phylogenetic inference on protein data sets. Our main contributions are 1) utilizing 64-bit variables, thereby enabling ta(MC)(3) to process larger data sets than MrBayes; and 2) to use Kahan summation to improve accuracy, convergence rates, and consequently runtime. Versus the current fastest software, we achieve a speedup of up to around 2.5 (and up to around 90 vs. serial MrBayes), and more on multi-GPU hardware. GPU MrBayes V3.1 is available from http://sourceforge.net/projects/mrbayes-gpu/.

Error, signal, and the placement of Ctenophora sister to all other animals

Elucidating relationships among early animal lineages has been difficult, and recent phylogenomic analyses place Ctenophora sister to all other extant animals, contrary to the traditional view of Porifera as the earliest-branching animal lineage. To date, phylogenetic support for either ctenophores or sponges as sister to other animals has been limited and inconsistent among studies. Lack of agreement among phylogenomic analyses using different data and methods obscures how complex traits, such as epithelia, neurons, and muscles evolved. A consensus view of animal evolution will not be accepted until datasets and methods converge on a single hypothesis of early metazoan relationships and putative sources of systematic error (e.g., long-branch attraction, compositional bias, poor model choice) are assessed. Here, we investigate possible causes of systematic error by expanding taxon sampling with eight novel transcriptomes, strictly enforcing orthology inference criteria, and progressively examining potential causes of systematic error while using both maximum-likelihood with robust data partitioning and Bayesian inference with a site-heterogeneous model. We identified ribosomal protein genes as possessing a conflicting signal compared with other genes, which caused some past studies to infer ctenophores and cnidarians as sister. Importantly, biases resulting from elevated compositional heterogeneity or elevated substitution rates are ruled out. Placement of ctenophores as sister to all other animals, and sponge monophyly, are strongly supported under multiple analyses, herein.

Phylogenomic analyses uncover origin and spread of the Wolbachia pandemic

Of all obligate intracellular bacteria, Wolbachia is probably the most common. In general, Wolbachia are either widespread, opportunistic reproductive parasites of arthropods or essential mutualists in a single group of filarial nematodes, including many species of medical significance. To date, a robust phylogenetic backbone of Wolbachia is lacking and consequently, many Wolbachia-related phenomena cannot be discussed in a broader evolutionary context. Here we present the first comprehensive phylogenomic analysis of Wolbachia supergroup relationships based on new whole-genome-shotgun data. Our results suggest that Wolbachia has switched between its two major host groups at least twice. The ability of some arthropod-infecting Wolbachia to universally infect and to adapt to a broad range of hosts quickly is restricted to a single monophyletic lineage (containing supergroups A and B). Thus, the currently observable pandemic has likely a single evolutionary origin and is unique within the radiation of Wolbachia strains.

The molecular symplesiomorphies shared by the stem groups of metazoan evolution: can sites as few as 1% have a significant impact on recognizing the phylogenetic position of myzostomida?

Although it is clear that taxon sampling, alignments, gene sampling, tree reconstruction methods and the total length of the sequences used are critical to the reconstruction of evolutionary history, weakly supported or misleading nodes exist in phylogenetic studies with no obvious flaw in those aspects. The phylogenetic studies focusing on the basal part of bilaterian evolution are such a case. During the past decade, Myzostomida has appeared in the basal part of Bilateria in several phylogenetic studies of Metazoa. However, most researchers have entertained only two competing hypotheses about the position of Myzostomida-an affinity with Annelida and an affinity with Platyhelminthes. In this study, dozens of symplesiomorphies were discovered by means of ancestral state reconstruction in the complete 18S and 28S rDNAs shared by the stem groups of Metazoa. By contrastive analysis on the datasets with or without such symplesiomorphic sites, we discovered that Myzostomida and other basal groups are basal lineages of Bilateria due to the corresponding symplesiomorphies shared with earlier lineages. As such, symplesiomorphies account for approximately 1-2% of the whole dataset have an essential impact on phylogenetic inference, and this study reminds molecular systematists of the importance of carrying out ancestral state reconstruction at each site in sequence-based phylogenetic studies. In addition, reasons should be explored for the low support of the hypothesis that Myzostomida belongs to Annelida in the results of phylogenomic studies. Future phylogenetic studies concerning Myzostomida should include all of the basal lineages of Bilateria to avoid directly neglecting the stand-alone basal position of Myzostomida as a potential hypothesis.

TreSpEx-Detection of Misleading Signal in Phylogenetic Reconstructions Based on Tree Information

Phylogenies of species or genes are commonplace nowadays in many areas of comparative biological studies. However, for phylogenetic reconstructions one must refer to artificial signals such as paralogy, long-branch attraction, saturation, or conflict between different datasets. These signals might eventually mislead the reconstruction even in phylogenomic studies employing hundreds of genes. Unfortunately, there has been no program allowing the detection of such effects in combination with an implementation into automatic process pipelines. TreSpEx (Tree Space Explorer) now combines different approaches (including statistical tests), which utilize tree-based information like nodal support or patristic distances (PDs) to identify misleading signals. The program enables the parallel analysis of hundreds of trees and/or predefined gene partitions, and being command-line driven, it can be integrated into automatic process pipelines. TreSpEx is implemented in Perl and supported on Linux, Mac OS X, and MS Windows. Source code, binaries, and additional material are freely available at http://www.annelida.de/research/bioinformatics/software.html.

Evaluating topological conflict in centipede phylogeny using transcriptomic data sets

Relationships between the five extant orders of centipedes have been considered solved based on morphology. Phylogenies based on samples of up to a few dozen genes have largely been congruent with the morphological tree apart from an alternative placement of one order, the relictual Craterostigmomorpha, consisting of two species in Tasmania and New Zealand. To address this incongruence, novel transcriptomic data were generated to sample all five orders of centipedes and also used as a test case for studying gene-tree incongruence. Maximum likelihood and Bayesian mixture model analyses of a data set composed of 1,934 orthologs with 45% missing data, as well as the 389 orthologs in the least saturated, stationary quartile, retrieve strong support for a sister-group relationship between Craterostigmomorpha and all other pleurostigmophoran centipedes, of which the latter group is newly named Amalpighiata. The Amalpighiata hypothesis, which shows little gene-tree incongruence and is robust to the influence of among-taxon compositional heterogeneity, implies convergent evolution in several morphological and behavioral characters traditionally used in centipede phylogenetics, such as maternal brood care, but accords with patterns of first appearances in the fossil record.

Transcriptome data reveal Syndermatan relationships and suggest the evolution of endoparasitism in Acanthocephala via an epizoic stage

The taxon Syndermata comprises the biologically interesting wheel animals ("Rotifera": Bdelloidea + Monogononta + Seisonidea) and thorny-headed worms (Acanthocephala), and is central for testing superordinate phylogenetic hypotheses (Platyzoa, Gnathifera) in the metazoan tree of life. Recent analyses of syndermatan phylogeny suggested paraphyly of Eurotatoria (free-living bdelloids and monogononts) with respect to endoparasitic acanthocephalans. Data of epizoic seisonids, however, were absent, which may have affected the branching order within the syndermatan clade. Moreover, the position of Seisonidea within Syndermata should help in understanding the evolution of acanthocephalan endoparasitism. Here, we report the first phylogenomic analysis that includes all four higher-ranked groups of Syndermata. The analyzed data sets comprise new transcriptome data for Seison spec. (Seisonidea), Brachionus manjavacas (Monogononta), Adineta vaga (Bdelloidea), and Paratenuisentis ambiguus (Acanthocephala). Maximum likelihood and Bayesian trees for a total of 19 metazoan species were reconstructed from up to 410 functionally diverse proteins. The results unanimously place Monogononta basally within Syndermata, and Bdelloidea appear as the sister group to a clade comprising epizoic Seisonidea and endoparasitic Acanthocephala. Our results support monophyly of Syndermata, Hemirotifera (Bdelloidea + Seisonidea + Acanthocephala), and Pararotatoria (Seisonidea + Acanthocephala), rejecting monophyly of traditional Rotifera and Eurotatoria. This serves as an indication that early acanthocephalans lived epizoically or as ectoparasites on arthropods, before their complex lifecycle with arthropod intermediate and vertebrate definite hosts evolved.

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.

Exploiting gene families for phylogenomic analysis of myzostomid transcriptome data

BACKGROUND

In trying to understand the evolutionary relationships of organisms, the current flood of sequence data offers great opportunities, but also reveals new challenges with regard to data quality, the selection of data for subsequent analysis, and the automation of steps that were once done manually for single-gene analyses. Even though genome or transcriptome data is available for representatives of most bilaterian phyla, some enigmatic taxa still have an uncertain position in the animal tree of life. This is especially true for myzostomids, a group of symbiotic (or parasitic) protostomes that are either placed with annelids or flatworms.

METHODOLOGY

Based on similarity criteria, Illumina-based transcriptome sequences of one myzostomid were compared to protein sequences of one additional myzostomid and 29 reference metazoa and clustered into gene families. These families were then used to investigate the phylogenetic position of Myzostomida using different approaches: Alignments of 989 sequence families were concatenated, and the resulting superalignment was analyzed under a Maximum Likelihood criterion. We also used all 1,878 gene trees with at least one myzostomid sequence for a supertree approach: the individual gene trees were computed and then reconciled into a species tree using gene tree parsimony.

CONCLUSIONS

Superalignments require strictly orthologous genes, and both the gene selection and the widely varying amount of data available for different taxa in our dataset may cause anomalous placements and low bootstrap support. In contrast, gene tree parsimony is designed to accommodate multilocus gene families and therefore allows a much more comprehensive data set to be analyzed. Results of this supertree approach showed a well-resolved phylogeny, in which myzostomids were part of the annelid radiation, and major bilaterian taxa were found to be monophyletic.

A molecular phylogeny of bryozoans

We present the most comprehensive molecular phylogeny of bryozoans to date. Our concatenated alignment of two nuclear ribosomal and five mitochondrial genes includes 95 taxa and 13,292 nucleotide sites, of which 8297 were included. The number of new sequences generated during this project are for each gene:ssrDNA (32), lsrDNA (22), rrnL (38), rrnS (35), cox1 (37), cox3 (34), and cytb (44). Our multi-gene analysis provides a largely stable topology across the phylum. The major groups were unambiguously resolved as (Phylactolaemata (Cyclostomata (Ctenostomata, Cheilostomata))), with Ctenostomata paraphyletic. Within Phylactolaemata, (Stephanellidae, Lophopodidae) form the earliest divergent clade. Fredericellidae is not resolved as a monophyletic family and forms a clade together with Plumatellidae, Cristatellidae and Pectinatellidae, with the latter two as sister taxa. Hyalinella and Gelatinella nest within the genus Plumatella. Cyclostome taxa fall into three major clades: i. (Favosipora (Plagioecia, Rectangulata)); ii. (Entalophoroecia ((Diplosolen, Cardioecia) (Frondipora, Cancellata))); and iii. (Articulata ((Annectocyma, Heteroporidae) (Tubulipora (Tennysonia, Idmidronea)))), with suborders Tubuliporina and Cerioporina, and family Plagioeciidae each being polyphyletic. Ctenostomata is composed of three paraphyletic clades to the inclusion of Cheilostomata: ((Alcyonidium, Flustrellidra) (Paludicella (Anguinella, Triticella)) (Hislopia (Bowerbankia, Amathia)) Cheilostomata); Flustrellidra nests within the genus Alcyonidium, and Amathia nests within the genus Bowerbankia. Suborders Carnosa and Stolonifera are not monophyletic. Within the cheilostomes, Malacostega is paraphyletic to the inclusion of all other cheilostomes. Conopeum is the most early divergent cheilostome, forming the sister group to ((Malacostega, Scrupariina, Inovicellina) ((Hippothoomorpha, Flustrina) (Lepraliomorpha, Umbonulomorpha))); Flustrina is paraphyletic to the inclusion of the hippothoomorphs; neither Lepraliomorpha nor Umbonulomorpha is monophyletic. Ascophorans are polyphyletic, with hippothoomorphs grouping separately from lepraliomorphs and umbonulomorphs; no cribrimorphs were included in the analysis. Results are discussed in the light of molecular and morphological evidence. Ancestral state reconstruction of larval strategy in Gymnolaemata revealed planktotrophy and lecithotrophy as equally parsimonious solutions for the ancestral condition. More comprehensive taxon sampling is expected to clarify this result. We discuss the extent of non-bryozoan contaminant sequences deposited in GenBank and their impact on the reconstruction of metazoan phylogenies and those of bryozoan interrelationships.

The complete mitochondrial genome of Flustra foliacea (Ectoprocta, Cheilostomata) - compositional bias affects phylogenetic analyses of lophotrochozoan relationships

Background

The phylogenetic relationships of the lophophorate lineages, ectoprocts, brachiopods and phoronids, within Lophotrochozoa are still controversial. We sequenced an additional mitochondrial genome of the most species-rich lophophorate lineage, the ectoprocts. Although it is known that there are large differences in the nucleotide composition of mitochondrial sequences of different lineages as well as in the amino acid composition of the encoded proteins, this bias is often not considered in phylogenetic analyses. We applied several approaches for reducing compositional bias and saturation in the phylogenetic analyses of the mitochondrial sequences.

Results

The complete mitochondrial genome (16,089 bp) of Flustra foliacea (Ectoprocta, Gymnolaemata, Cheilostomata) was sequenced. All protein-encoding, rRNA and tRNA genes are transcribed from the same strand. Flustra shares long intergenic sequences with the cheilostomate ectoproct Bugula, which might be a synapomorphy of these taxa. Further synapomorphies might be the loss of the DHU arm of the tRNA L(UUR), the loss of the DHU arm of the tRNA S(UCN) and the unique anticodon sequence GAG of the tRNA L(CUN). The gene order of the mitochondrial genome of Flustra differs strongly from that of the other known ectoprocts. Phylogenetic analyses of mitochondrial nucleotide and amino acid data sets show that the lophophorate lineages are more closely related to trochozoan phyla than to deuterostomes or ecdysozoans confirming the Lophotrochozoa hypothesis. Furthermore, they support the monophyly of Cheilostomata and Ectoprocta. However, the relationships of the lophophorate lineages within Lophotrochozoa differ strongly depending on the data set and the used method. Different approaches for reducing heterogeneity in nucleotide and amino acid data sets and saturation did not result in a more robust resolution of lophotrochozoan relationships.

Conclusion

The contradictory and usually weakly supported phylogenetic reconstructions of the relationships among lophotrochozoan phyla based on mitochondrial sequences indicate that these alone do not contain enough information for a robust resolution of the relations of the lophotrochozoan phyla. The mitochondrial gene order is also not useful for inferring their phylogenetic relationships, because it is highly variable in ectoprocts, brachiopods and some other lophotrochozoan phyla. However, our study revealed several rare genomic changes like the evolution of long intergenic sequences and changes in the structure of tRNAs, which may be helpful for reconstructing ectoproct phylogeny.

Loss of different inverted repeat copies from the chloroplast genomes of Pinaceae and cupressophytes and influence of heterotachy on the evaluation of gymnosperm phylogeny

The relationships among the extant five gymnosperm groups--gnetophytes, Pinaceae, non-Pinaceae conifers (cupressophytes), Ginkgo, and cycads--remain equivocal. To clarify this issue, we sequenced the chloroplast genomes (cpDNAs) from two cupressophytes, Cephalotaxus wilsoniana and Taiwania cryptomerioides, and 53 common chloroplast protein-coding genes from another three cupressophytes, Agathis dammara, Nageia nagi, and Sciadopitys verticillata, and a non-Cycadaceae cycad, Bowenia serrulata. Comparative analyses of 11 conifer cpDNAs revealed that Pinaceae and cupressophytes each lost a different copy of inverted repeats (IRs), which contrasts with the view that the same IR has been lost in all conifers. Based on our structural finding, the character of an IR loss no longer conflicts with the "gnepines" hypothesis (gnetophytes sister to Pinaceae). Chloroplast phylogenomic analyses of amino acid sequences recovered incongruent topologies using different tree-building methods; however, we demonstrated that high heterotachous genes (genes that have highly different rates in different lineages) contributed to the long-branch attraction (LBA) artifact, resulting in incongruence of phylogenomic estimates. Additionally, amino acid compositions appear more heterogeneous in high than low heterotachous genes among the five gymnosperm groups. Removal of high heterotachous genes alleviated the LBA artifact and yielded congruent and robust tree topologies in which gnetophytes and Pinaceae formed a sister clade to cupressophytes (the gnepines hypothesis) and Ginkgo clustered with cycads. Adding more cupressophyte taxa could not improve the accuracy of chloroplast phylogenomics for the five gymnosperm groups. In contrast, removal of high heterotachous genes from data sets is simple and can increase confidence in evaluating the phylogeny of gymnosperms.

Statistics and truth in phylogenomics

Phylogenomics refers to the inference of historical relationships among species using genome-scale sequence data and to the use of phylogenetic analysis to infer protein function in multigene families. With rapidly decreasing sequencing costs, phylogenomics is becoming synonymous with evolutionary analysis of genome-scale and taxonomically densely sampled data sets. In phylogenetic inference applications, this translates into very large data sets that yield evolutionary and functional inferences with extremely small variances and high statistical confidence (P value). However, reports of highly significant P values are increasing even for contrasting phylogenetic hypotheses depending on the evolutionary model and inference method used, making it difficult to establish true relationships. We argue that the assessment of the robustness of results to biological factors, that may systematically mislead (bias) the outcomes of statistical estimation, will be a key to avoiding incorrect phylogenomic inferences. In fact, there is a need for increased emphasis on the magnitude of differences (effect sizes) in addition to the P values of the statistical test of the null hypothesis. On the other hand, the amount of sequence data available will likely always remain inadequate for some phylogenomic applications, for example, those involving episodic positive selection at individual codon positions and in specific lineages. Again, a focus on effect size and biological relevance, rather than the P value, may be warranted. Here, we present a theoretical overview and discuss practical aspects of the interplay between effect sizes, bias, and P values as it relates to the statistical inference of evolutionary truth in phylogenomics.

A phylogenomic approach to resolve the arthropod tree of life

Arthropods were the first animals to conquer land and air. They encompass more than three quarters of all described living species. This extraordinary evolutionary success is based on an astoundingly wide array of highly adaptive body organizations. A lack of robustly resolved phylogenetic relationships, however, currently impedes the reliable reconstruction of the underlying evolutionary processes. Here, we show that phylogenomic data can substantially advance our understanding of arthropod evolution and resolve several conflicts among existing hypotheses. We assembled a data set of 233 taxa and 775 genes from which an optimally informative data set of 117 taxa and 129 genes was finally selected using new heuristics and compared with the unreduced data set. We included novel expressed sequence tag (EST) data for 11 species and all published phylogenomic data augmented by recently published EST data on taxonomically important arthropod taxa. This thorough sampling reduces the chance of obtaining spurious results due to stochastic effects of undersampling taxa and genes. Orthology prediction of genes, alignment masking tools, and selection of most informative genes due to a balanced taxa-gene ratio using new heuristics were established. Our optimized data set robustly resolves major arthropod relationships. We received strong support for a sister group relationship of onychophorans and euarthropods and strong support for a close association of tardigrades and cycloneuralia. Within pancrustaceans, our analyses yielded paraphyletic crustaceans and monophyletic hexapods and robustly resolved monophyletic endopterygote insects. However, our analyses also showed for few deep splits that were recently thought to be resolved, for example, the position of myriapods, a remarkable sensitivity to methods of analyses.

Phylogeny of Annelida (Lophotrochozoa): total-evidence analysis of morphology and six genes

BACKGROUND

Annelida is one of the major protostome phyla, whose deep phylogeny is very poorly understood. Recent molecular phylogenies show that Annelida may include groups once considered separate phyla (Pogonophora, Echiurida, and Sipunculida) and that Clitellata are derived polychaetes. SThe "total-evidence" analyses combining morphological and molecular characters have been published for a few annelid taxa. No attempt has yet been made to analyse simultaneously morphological and molecular information concerning the Annelida as a whole.

RESULTS

Phylogenetic relationships within Annelida were analysed on the basis of 93 morphological characters and sequences of six genes (18S, 28S, and 16S rRNA, EF1alpha, H3, COI), altogether, 87 terminals of all annelid "families" and 3,903 informative characters, by Bayesian and maximum-parsimony methods. The analysis of the combined dataset yields the following scheme of relationships: Phyllodocida and Eunicida are monophyletic groups, together probably forming monophyletic Aciculata (incl. Orbiniidae and Parergodrilidae that form a sister group of the Eunicida). The traditional "Scolecida" and "Canalipalpata" are both polyphyletic, forming instead two clades: one including Cirratuliformia and the "sabelloid-spionoid clade" (incl. Sternaspis, Sabellidae-Serpulidae, Sabellariidae, Spionida s.str.), the other ("terebelloid-capitelloid clade") including Terebelliformia, Arenicolidae-Maldanidae, and Capitellidae-Echiurida. The Clitellata and "clitellate-like polychaetes" (Aeolosomatidae, Potamodrilidae, Hrabeiella) form a monophyletic group. The position of the remaining annelid groups is uncertain--the most problematic taxa are the Opheliidae-Scalibregmatidae clade, the Amphinomida-Aberranta clade, Apistobranchus, Chaetopteridae, Myzostomida, the Sipunculida-Dinophilidae clade, and the "core Archiannelida" (= Protodrilidae, Nerillidae, Polygordiidae, Saccocirridae).

CONCLUSION

The combined ("total-evidence") phylogenetic analysis provides a modified view of annelid evolution, with several higher-level taxa, i.e. Phyllodocida, Eunicida, orbinioid-parergodrilid clade (OPC), Cirratuliformia, sabelloid-spionoid clade (SSC), terebelloid-capitelloid clade (TCC), and "Clitellatomorpha". Two unorthodox clades, the "core Archiannelida" and Sipunculida-Dinophilidae, are proposed. Although the deep-level evolutionary relationships of Annelida remain poorly understood, we propose the monophyly of the Aciculata, sister-group relationships between the Eunicida and OPC, between the Cirratuliformia and SSC, and possibly also between the "Clitellatomorpha" and Oweniidae-Pogonophora clades.

Support for the monophyletic origin of Gnathifera from phylogenomics

The monophyletic origin of Spiralia within the metazoan tree of life is supported by many large-scale phylogenomic data. While there is now substantial molecular evidence for Lophotrochozoa being a monophyletic taxon within Spiralia, the phylogenetic affiliations of many other spiralian phyla remain unclear. Here we focus on the question of a monophyletic taxon Gnathifera, which was originally characterized by jaw morphology as comprising the taxa Rotifera, Acanthocephala and Gnathostomulida. Based on a large-scale molecular sequence dataset of 11,146 amino acid residues, we reconstructed phylogenetic trees of spiralian phyla using maximum-likelihood and Bayesian approaches. We obtain the first phylogenomic evidence for the clade Gnathifera, linking Syndermata (Rotifera+Acanthocephala) with Gnathostomulida. Furthermore, our data support recent findings concerning the paraphyly of Eurotatoria.