First Sequenced Mitochondrial Genome from the Phylum Acanthocephala(Leptorhynchoides thecatus) and Its Phylogenetic Position Within Metazoa

Characterization of the complete mitochondrial genomes of the zoonotic parasites Bolbosoma nipponicum and Corynosoma villosum (Acanthocephala: Polymorphida) and the molecular phylogeny of the order Polymorphida

Acanthocephalans of the order Polymorphida mainly parasitic in birds and mammals, are of veterinary, medical and economic importance. However, the evolutionary relationships of its 3 families (Centrorhynchidae, Polymorphidae and Plagiorhynchidae) remain under debate. Additionally, some species of Polymorphida (i.e. Bolbosoma spp. and Corynosoma spp.) are recognized as zoonotic parasites, associated with human acanthocephaliasis, but the mitochondrial genomes for representatives of Bolbosoma and Corynosoma have not been reported so far. In the present study, the complete mitochondrial genomes B. nipponicum and C. villosum (Acanthocephala: Polymorphidae) are reported for the first time, which are 14 296 and 14 241 bp in length, respectively, and both contain 36 genes [including 12 PCGs, 22 tRNA genes and 2 rRNA genes] and 2 non-coding regions (NCR1 and NCR2). The gene arrangement of some tRNAs in the mitogenomes of B. nipponicum and C. villosum differs from that found in all other acanthocephalans, except Polymorphus minutus. Phylogenetic results based on concatenated amino acid (AA) sequences of the 12 protein-coding genes (PCGs) strongly supported that the family Polymorphidae is a sister to the Centrorhynchidae rather than the Plagiorhynchidae, and also confirmed the sister relationship of the genera Bolbosoma and Corynosoma in the Polymorphidae based on the mitogenomic data for the first time. Our present findings further clarified the phylogenetic relationships of the 3 families Plagiorhynchidae, Centrorhynchidae and Polymorphidae, enriched the mitogenome data of the phylum Acanthocephala (especially the order Polymorphida), and provided the resource of genetic data for diagnosing these 2 pathogenic parasites of human acanthocephaliasis.

Characterization of the complete mitochondrial genome of Acanthogyrus (Acanthosentis) bilaspurensis Chowhan, Gupta & Khera, 1987 (Eoacanthocephala: Quadrigyridae), the smallest mitochondrial genome in Acanthocephala, and its phylogenetic implications

The phylum Acanthocephala is an important group of parasites with more than 1,300 species parasitizing intestine of all major vertebrate groups. However, our present knowledge of the mitochondrial genomes of Acanthocephala remains very limited. In the present study, we sequenced and annotated the complete mitochondrial genome of Acanthogyrus (Acanthosentis) bilaspurensis (Gyracanthocephala: Quadrigyridae) for the first time based on the specimens recovered from the intestine of common carp Cyprinus carpio Linnaeus (Cyprinidae) in Pakistan. The mitochondrial genome of A. bilaspurensis is 13,360 bp in size and contains 36 genes, representing the smallest mitogenome of acanthocephalans reported so far. The mitogenome of A. bilaspurensis also has the lowest level of overall A+T contents (59.3%) in the mitogenomes of Eoacanthocephala, and the non-coding region 3 (NCR3) lies between trnS2 and trnI, which is different from all of the other acanthocephalan species. Phylogenetic analyses based on concatenating the amino acid sequences of 12 protein-coding genes using maximum likelihood (ML) and Bayesian inference (BI) methods revealed that the family Pseudoacanthocephalidae is a sister to the Arhythmacanthidae rather than the Cavisomatidae, and the families Rhadinorhynchidae and Cavisomatidae showed sister relationships.

Two-Fold ND5 Genes, Three-Fold Control Regions, lncRNA, and the "Missing" ATP8 Found in the Mitogenomes of Polypedates megacephalus (Rhacophridae: Polypedates)

In prior research on the mitochondrial genome (mitogenome) of Polypedates megacephalus, the one copy of ND5 gene was translocated to the control region (CR) and the ATP8 gene was not found. Gene loss is uncommon among vertebrates. However, in this study, we resequenced the mitogenomes of P. megacephalus from different regions using a "primer bridging" approach with Sanger sequencing technologies, which revealed the "missing" ATP8 gene in P. megacephalus as well as three other previously published Polypedates. The mitogenome of this species was found to contain two copies of the ND5 genes and three copies of the control regions. Furthermore, multiple tandem repeats were identified in the control regions. Notably, we observed that there was no correlation between genetic divergence and geographic distance. However, using the mitogenome, gene expression analysis was performed via RT-qPCR of liver samples and it was thus determined that COIII, ND2, ND4, and ND6 were reduced to 0.64 ± 0.24, 0.55 ± 0.34, 0.44 ± 0.21 and 0.65 ± 0.17, respectively, under low-temperature stress (8 °C) as compared with controls (p < 0.05). Remarkably, the transcript of long non-coding RNA (lncRNA) between positions 8029 and 8612 decreased significantly with exposure to low-temperature stress (8 °C). Antisense ND6 gene expression showed a downward trend, but this was not significant. These results reveal that modulations of protein-coding mitochondrial genes and lncRNAs of P. megacephalus play a crucial role in the molecular response to cold stress.

Phylomitogenomic Analyses Provided Further Evidence for the Resurrection of the Family Pseudoacanthocephalidae (Acanthocephala: Echinorhynchida)

The phylum Acanthocephala is an important monophyletic group of parasites, with adults parasitic in the digestive tracts of all major vertebrate groups. Acanthocephalans are of veterinary, medical, and economic importance due to their ability to cause disease in domestic animals, wildlife, and humans. However, the current genetic data for acanthocephalans are sparse, both in terms of the proportion of taxa surveyed and the number of genes sequenced. Consequently, the basic molecular phylogenetic framework for the phylum is still incomplete. In the present study, we reported the first complete mitochondrial genome from a representative of the family Pseudoacanthocephalidae Petrochenko, 1956. The mitogenome of Pseudoacanthocephalus bufonis (Shipley, 1903) is 14,056 bp in length, contains 36 genes (12 protein-coding genes (PCGs) (lacking atp8), 22 tRNA genes, and 2 rRNA genes (rrnL and rrnS)) and two non-coding regions (NCR1 and NCR2), and displayed the highest GC-skew in the order Echinorhynchida. Phylogenetic results of maximum likelihood (ML) and Bayesian inference (BI) using the amino acid sequences of 12 protein-coding genes in different models provided further evidence for the resurrection of the family Pseudoacanthocephalidae and also supported that the order Echinorhynchida is paraphyletic. A monophyletic clade comprising P. bufonis and Cavisoma magnum suggests a close affinity between Pseudoacanthocephalidae and Cavisomatidae. Our phylogenetic analyses also showed that Polymorphidae has a closer relationship with Centrorhynchidae than Plagiorhynchidae in the monophyletic order Polymorphida.

The mitochondrial genome of Heterosentis pseudobagri (Wang & Zhang, 1987) Pichelin & Cribb, 1999 reveals novel aspects of tRNA genes evolution in Acanthocephala

BACKGROUND

Acanthocephala is a clade of obligate endoparasites whose mitochondrial genomes (mitogenomes) and evolution remain relatively poorly understood. Previous studies reported that atp8 is lacking from acanthocephalan mitogenomes, and that tRNA genes often have nonstandard structures. Heterosentis pseudobagri (Arhythmacanthidae) is an acanthocephalan fish endoparasite for which no molecular data are currently available, and biological information is unavailable in the English language. Furthermore, there are currently no mitogenomes available for Arhythmacanthidae.

METHODS

We sequenced its mitogenome and transcriptome, and conducted comparative mitogenomic analyses with almost all available acanthocephalan mitogenomes.

RESULTS

The mitogenome had all genes encoded on the same strand and unique gene order in the dataset. Among the 12 protein-coding genes, several genes were highly divergent and annotated with difficulty. Moreover, several tRNA genes could not be identified automatically, so we had to identify them manually via a detailed comparison with orthologues. As common in acanthocephalans, some tRNAs lacked either the TWC arm or the DHU arm, but in several cases, we annotated tRNA genes only on the basis of the conserved narrow central segment comprising the anticodon, while the flanking 5' and 3' ends did not exhibit any resemblance to orthologues and they could not be folded into a tRNA secondary structure. We corroborated that these are not sequencing artefacts by assembling the mitogenome from transcriptomic data. Although this phenomenon was not observed in previous studies, our comparative analyses revealed the existence of highly divergent tRNAs in multiple acanthocephalan lineages.

CONCLUSIONS

These findings indicate either that multiple tRNA genes are non-functional or that (some) tRNA genes in (some) acanthocephalans might undergo extensive posttranscriptional tRNA processing which restores them to more conventional structures. It is necessary to sequence mitogenomes from yet unrepresented lineages and further explore the unusual patterns of tRNA evolution in Acanthocephala.

Mitochondrial phylogenomics of Acanthocephala: nucleotide alignments produce long-branch attraction artefacts

BACKGROUND

Classification of the Acanthocephala, a clade of obligate endoparasites, remains unresolved because of insufficiently strong resolution of morphological characters and scarcity of molecular data with a sufficient resolution. Mitochondrial genomes may be a suitable candidate, but they are available for a small number of species and their suitability for the task has not been tested thoroughly.

METHODS

Herein, we sequenced the first mitogenome for the large family Rhadinorhynchidae: Micracanthorhynchina dakusuiensis. These are also the first molecular data generated for this entire genus. We conducted a series of phylogenetic analyses using concatenated nucleotides (NUC) and amino acids (AAs) of all 12 protein-coding genes, three different algorithms, and the entire available acanthocephalan mitogenomic dataset.

RESULTS

We found evidence for strong compositional heterogeneity in the dataset, and Micracanthorhynchina dakusuiensis exhibited a disproportionately long branch in all analyses. This caused a long-branch attraction artefact (LBA) of M. dakusuiensis resolved at the base of the Echinorhynchida clade when the NUC dataset was used in combination with standard phylogenetic algorithms, maximum likelihood (ML) and Bayesian inference (BI). Both the use of the AA dataset (BI-AAs and ML-AAs) and the CAT-GTR model designed for suppression of LBA (CAT-GTR-AAs and CAT-GTR-NUC) at least partially attenuated this LBA artefact. The results support Illiosentidae as the basal radiation of Echinorhynchida and Rhadinorhynchidae forming a clade with Echinorhynchidae and Pomporhynchidae. The questions of the monophyly of Rhadinorhynchidae and its sister lineage remain unresolved. The order Echinorhynchida was paraphyletic in all of our analyses.

CONCLUSIONS

Future studies should take care to attenuate compositional heterogeneity-driven LBA artefacts when applying mitogenomic data to resolve the phylogeny of Acanthocephala.

Complete Mitogenomes of Polypedates Tree Frogs Unveil Gene Rearrangement and Concerted Evolution within Rhacophoridae

New developments in sequencing technology and nucleotide analysis have allowed us to make great advances in reconstructing anuran phylogeny. As a clade of representative amphibians that have radiated from aquatic to arboreal habitats, our understanding of the systematic status and molecular biology of rhacophorid tree frogs is still limited. We determined two new mitogenomes for the genus Polypedates (Rhacophoridae): P. impresus and P. mutus. We conducted comparative and phylogenetic analyses using our data and seven other rhacophorid mitogenomes. The mitogenomes of the genera Polypedates, Buergeria, and Zhangixalus were almost identical, except that the ATP8 gene in Polypedates had become a non-coding region; Buergeria maintained the legacy "LTPF" tRNA gene cluster compared to the novel "TLPF" order in the other two genera; and B. buergeri and Z. dennysi had no control region (CR) duplication. The resulting phylogenetic relationship supporting the above gene rearrangement pathway suggested parallel evolution of ATP8 gene loss of function (LoF) in Polypedates and CR duplication with concerted evolution of paralogous CRs in rhacophorids. Finally, conflicting topologies in the phylograms of 185 species reflected the advantages of phylogenetic analyses using multiple loci.

Molecular Characterization of a New Moniliformis sp. From a Plateau Zokor (Eospalax fontanierii baileyi) in China

In the present study, a new species of the genus Moniliformis species is described taxonomically in the mitochondrial genomic context. The parasite was found in a plateau zokor captured in a high-altitude area of Xiahe County of Gansu Province, China. The mitochondrial (mt) genome length of this new species was 14,066 bp comprising 36 genes and 2 additional non-coding regions (SNR and LNR), without atp8. The molecular phylogeny inferred by the cytochrome c oxidase subunit I gene (cox1) and the18S ribosomal RNA gene (18S rDNA) sequences showed that the parasite as a sister species to other Moniliformis spp. and was named Moniliformis sp. XH-2020. The phylogeny of the concatenated amino acid sequences of the 12 protein-coding genes (PCGs) showed Moniliformis sp. XH-2020 in the same cluster as Macracanthorhynchus hirudinaceus and Oncicola luehei confirming the cox1 and 18S rDNA phylogenetic inference. In addition, the entire mt genome sequenced in this study represents the first in the order Moniliformida, providing molecular material for further study of the phylogeny of the class Archiacanthocephala. Moreover, the species of this class, use arthropods as intermediate hosts and mammals as definitive hosts and are agents of acanthocephaliasis, a zoonosis in humans. Therefore, this study not only expands the host range among potential wild animal hosts for Archiacanthocephalans which is of great ecological and evolutionary significance but also has important significance for the research of zoonotic parasitic diseases.

Morphological and complete mitogenomic characterisation of the acanthocephalan Polymorphus minutus infecting the duck Anas platyrhynchos

Morphological characteristics of the acanthocephalan Polymorphus minutus (Goeze, 1782), which was collected from the duck Anas platyrhynchos Linnaeus in the Czech Republic, are described. The mitochondrial (mt) genome of P. minutus was sequenced, with a total length of 14,149 bp, comprising 36 genes including 12 protein coding genes (PCGs), 22 transfer RNA (tRNA) genes and two ribosomal RNA genes (rrnL and rrnS). This genome is similar to the mt genomes of other syndermatan species. All these genes were encoded on the same DNA strand and in the same orientation. The overall nucleotide composition of the P. minutus mt genome was 38.2% T, 27.3% G, 26.2% A, and 8.3% C. The amino acid sequences of 12 PCGs for mt genomes of 28 platyzoans, including P. minutus, were used for phylogenetic analysis, and the resulting topology recovers P. minutus as sister to Southwellina hispida (Van Cleave, 1925), and the two taxa form a sister clade to Centrorhynchus aluconis (Müller, 1780) and Plagiorhynchus transversus (Rudolphi, 1819), which are all species in the Palaeacanthocephala, thus supporting the monophyly of this class.

The genome, transcriptome, and proteome of the fish parasite Pomphorhynchus laevis (Acanthocephala)

Thorny-headed worms (Acanthocephala) are endoparasites exploiting Mandibulata (Arthropoda) and Gnathostomata (Vertebrata). Despite their world-wide occurrence and economic relevance as a pest, genome and transcriptome assemblies have not been published before. However, such data might hold clues for a sustainable control of acanthocephalans in animal production. For this reason, we present the first draft of an acanthocephalan nuclear genome, besides the mitochondrial one, using the fish parasite Pomphorhynchus laevis (Palaeacanthocephala) as a model. Additionally, we have assembled and annotated the transcriptome of this species and the proteins encoded. A hybrid assembly of long and short reads resulted in a near-complete P. laevis draft genome of ca. 260 Mb, comprising a large repetitive portion of ca. 63%. Numbers of transcripts and translated proteins (35,683) were within the range of other members of the Rotifera-Acanthocephala clade. Our data additionally demonstrate a significant reorganization of the acanthocephalan gene repertoire. Thus, more than 20% of the usually conserved metazoan genes were lacking in P. laevis. Ontology analysis of the retained genes revealed many connections to the incorporation of carotinoids. These are probably taken up via the surface together with lipids, thus accounting for the orange coloration of P. laevis. Furthermore, we found transcripts and protein sequences to be more derived in P. laevis than in rotifers from Monogononta and Bdelloidea. This was especially the case in genes involved in energy metabolism, which might reflect the acanthocephalan ability to use the scarce oxygen in the host intestine for respiration and simultaneously carry out fermentation. Increased plasticity of the gene repertoire through the integration of foreign DNA into the nuclear genome seems to be another underpinning factor of the evolutionary success of acanthocephalans. In any case, energy-related genes and their proteins may be considered as candidate targets for the acanthocephalan control.

Mitochondrial DNA dataset suggest that the genus Sphaerirostris Golvan, 1956 is a synonym of the genus Centrorhynchus Lühe, 1911

Our present genetic data of Acanthocephala, especially the mitochondrial (mt) genomes, remains very limited. In the present study, the nearly complete mt genome sequences of Sphaerirostris lanceoides (Petrochenko, 1949) was sequenced and determined for the first time based on specimens collected from the Indian pond heron Ardeola grayii (Sykes) (Ciconiiformes: Ardeidae) in Pakistan. The mt genome of S. lanceoides is 13 478 bp in size and contains 36 genes, including 12 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs) and two ribosomal RNA genes (rRNAs). Moreover, in order to clarify the phylogenetic relationship of the genera Centrorhynchus and Sphaerirostris, and to test the systematic position of S. lanceoides in the Centrorhynchidae, the phylogenetic analyses were performed using Bayesian inference and maximum likelihood methods, based on concatenated nucleotide sequences of 12 PCGs, rRNAs and tRNAs. The phylogenetic results further confirmed the monophyly of the order Polymorphida and the paraphyly of the order Echinorhynchida in the class Palaeacanthocephala. Our results also challenged the validity of the genus Sphaerirostris (Polymorphida: Centrorhynchidae) and showed a sister relationship between S. lanceoides and S. picae (Rudolphi, 1819).

Characterization of the complete mitogenome of Centrorhynchus clitorideus (Meyer, 1931) (Palaeacanthocephala: Centrorhynchidae), the largest mitochondrial genome in Acanthocephala, and its phylogenetic implications

Species of Centrorhynchus (Polymorphida: Centrorhynchidae) commonly parasitize various falconiform and strigiform birds worldwide. In the present study, the complete mitochondrial (mt) genome sequences of Centrorhynchus clitorideus was sequenced and annotated for the first time based on specimens collected from the little owl Athene noctua (Scopoli) (Strigiformes: Strigidae) in Pakistan. The complete mt genome sequences of C. clitorideus is 15,884 bp in length, and contained 36 genes [two rRNA genes (rrnL and rrnS), 22 tRNA genes and 12 protein-coding genes (PCGs) (lacking atp8)] and two non-coding regions (NCR1 and NCR2), which represents the largest mt genome of acanthocephalan reported so far. In order to assess the systematic position of C. clitorideus and the interrelationship of the family Centrorhynchidae and the other families in order Polymorphida, the phylogenetic tree was constructed using Bayesian inference (BI) based on amino acid sequences of 12 PCGs. Phylogenetic results supported C. clitorideus formed a sister relationship to C. milvus in Centrorhynchidae, which has a sister relationship to the representatives of Polymorphidae + Plagiorhynchidae. Our results revealed the monophyly of Polymorphida and paraphyly of Echinorhynchida in the class Palaeacanthocephala. The validity of the genus Sphaerirostris (Polymorphida: Centrorhynchidae) was also challenged by our phylogenetic results, which seems to be a synonym of Centrorhynchus. Moreover, the present phylogenetic analysis indicated that the family Quadrigyridae and subfamily Pallisentinae (A. cheni and P. celatus) are polyphyletic.

Characterization of the complete mitochondrial genome of Cavisoma magnum () (Acanthocephala: Palaeacanthocephala), first representative of the family Cavisomidae, and its phylogenetic implications

The phylum Acanthocephala is a small group of endoparasites occurring in the alimentary canal of all major lineages of vertebrates worldwide. In the present study, the complete mitochondrial (mt) genome of Cavisoma magnum (Southwell, 1927) (Palaeacanthocephala: Echinorhynchida) was determined and annotated, the representative of the family Cavisomidae with the characterization of the complete mt genome firstly decoded. The mt genome of this acanthocephalan is 13,594 bp in length, containing 36 genes plus two non-coding regions. The positions of trnV and SNCR (short non-coding region) in the mt genome of C. magnum are different comparing to those of the other acanthocephalan species available in GenBank. Phylogenetic analysis based on amino acid sequences of 12 protein-coding genes using Bayesian inference (BI) supported the class Palaeacanthocephala and its included order Polymorphida to be monophyletic, but rejected monophyly of the order Echinorhynchida. Our phylogenetic results also challenged the validity of the genus Sphaerirostris (Polymorphida: Centrorhynchidae). The novel mt genomic data of C. magnum are very useful for understanding the evolutionary history of this group of parasites and establishing a natural classification of Acanthocephala.

The revised complete mitogenome sequence of the tree frog Polypedatesmegacephalus (Anura, Rhacophoridae) by next-generation sequencing and phylogenetic analysis

The mitochondrial genome (mitogenome) sequence of the tree frog Polypedates megacephalus (16,473 bp) was previously reported as having the unusual characteristic of lacking the ND5 gene. In this study, a new mitogenome of P. megacephalus (19,952 bp) was resequenced using the next-generation sequencing (NGS) and standard Sanger sequencing technologies. It was discovered that the ND5 gene was not lost but translocated to the control region (CR) from its canonical location between the ND4 and ND6 genes. In addition, a duplicated control region was found in the new mitogenome of this species. Conservative region identification of the ND5 gene and phylogenetic analysis confirmed that the ND5 gene was located between two control regions. The phylogenetic relationship among 20 related species of anura revealed a rearrangement of the ND5 gene during the evolutionary process. These results also highlighted the advantages of next-generation sequencing. It will not only decrease the time and cost of sequencing, but also will eliminate the errors in published mitogenome databases.

Characterization of the complete mitochondrial genome of Brentisentisyangtzensis Yu & Wu, 1989 (Acanthocephala, Illiosentidae)

The mitogenome of Brentisentisyangtzensis is 13,864 bp in length and has the circular structure typical of metazoans. It contains 36 genes: 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and 12 protein-encoding genes (PCGs). All genes are transcribed from the same strand. Thirteen overlapping regions were found in the mitochondrial genome. The overall A+T content of B.yangtzensis is 68.3% versus 31.7% of G+C content (A = 27.8%, T = 40.5%, C = 9.0%, G = 22.7%). B.yangtzenensis (Illiosentidae) and Leptorhynchoidesthecatus (Rhadinorhynchidae) form a sister clade, showing the relatively close relationship between the Illiosentidae and the Rhadinorhynchidae. The mitochondrial gene arrangements of acanthocephalan species are relatively conserved, with only a few translocations of tRNAs (trnS1, trnS2, trnV, and trnK) detected. An identical gene order was found both in a sister clade (Centrorhynchusaluconis and Plagiorhynchustransversus) and across different classes (B.yangtzensis (Palaeacanthocephala), Acanthosentischeni (Eoacanthocephala) and Macracanthorhynchushirudinaceus (Archiacanthocephala), Oncicolaluehei and L.thecatus (Palaeacanthocephala)). More studies and more sequences of acanthocephalan species are needed to gain a clear understanding of the phylogenetic relationships.

Characterization of the complete mitochondrial genome of Centrorhynchus milvus (Acanthocephala: Polymorphida) and its phylogenetic implications

The phylum Acanthocephala is a small group of obligate parasites of animals. However, the current classifications of Acanthocephala are still under debate. Moreover, our present knowledge of the complete mitochondrial genome of this parasite group remains limited. To fill this knowledge gap, the complete mitochondrial (mt) genome of Centrorhynchus milvusWard, 1956 (Palaeacanthocephala: Polymorphida) was firstly sequenced and determined based on specimens collected from the red kite (Milvus milvus) in Pakistan. The mitochondrial genome of C. milvus is 14,314 bp in length and contains 36 genes, including 12 protein-coding (PCGs) genes, 22 transfer RNA (tRNA) genes and ribosomal RNA (rRNA) genes (rrnL and rrnS). To elucidate the phylogenetic relationships of the four classes of Acanthocephala and the systematic position of C. milvus, phylogenetic analysis based on concatenated amino acid sequences of 12 PCGs was performed using Bayesian inference (BI). The results supported the monophyly of Archiacanthocephala and Palaeacanthocephala with strong support (BPP = 1) and also indicated that Archiacanthocephala is the sister clade to the remaining classes of Acanthocephala (Palaeacanthocephala, Eoacanthocephala and Polyacanthocephala). However, Polyacanthocephala with only one representative species (P. caballeroi) is nested within Eoacanthocephala. Our phylogenetic analysis also confirmed C. milvus as the member of the family Centrorhynchidae with a sister relationship to C. aluconis. Our present mt genomic data are very useful for studying the molecular epidemiology, population genetics and systematics of acanthocephalans.

Characterization of the complete mitochondrial genome of Sphaerirostris picae (Rudolphi, 1819) (Acanthocephala: Centrorhynchidae), representative of the genus Sphaerirostris

The Centrorhynchidae (Acanthocephala: Palaeacanthocephala) is a cosmopolitan family commonly found in various avian and mammalian hosts. Within Centrorhynchidae, species of the genus Sphaerirostris Golvan, 1956 are usually parasitic in the digestive tract of various passerine birds. In the present study, adult specimens of Sphaerirostris picae (Rudolphi, 1819), the type species of this genus, were recovered from the small intestine of Acridotheres tristis (Sturnidae) and Dendrocitta vagabunda (Corvidae) in Pakistan. Molecular data from the nuclear or mitochondrial genome is either very limited or completely absent from this phylogenetically understudied group of acanthocephalans. To fill this knowledge gap, we sequenced and determined the internal transcribed spacers of ribosomal DNA (ITS rDNA) and the complete mitochondrial (mt) genome of S. picae. The ITS rDNA of S. picae was 95.2% similar to that of Sphaerirostris lanceoides which is the only member of the Centrorhynchidae whose ITS rDNA is available in GenBank. The phylogenetic tree based on the amino acid sequences of 12 mt protein-coding genes (PCGs) placed S. picae close to Centrorhynchus aluconis in a monophyletic clade of Polymorphida which also contain members of the families Polymorphidae and Plagiorhynchidae on separate branches. The mt gene arrangement, nucleotide composition and codon usage of 12 PCGs were discussed and compared with those of other acanthocephalan mt genomes. Within the Centrorhynchidae, S. picae and C. aluconis showed 67.7-86.8% similarity in the nucleotide sequences of 12 PCGs and 2 rRNAs, where nad4L is the most conserved gene while atp6 is the least conserved. The similarity in amino acid sequences ranged from 68.1 to 91.8%, where cox1 was recorded as the most conserved gene, while atp6 is highly variable among 12 PCGs. This novel mt genome of S. picae provides genetic resources for further studies of phylogenetics and molecular epidemiology of acanthocephalans.

Monogonont Rotifer, Brachionus calyciflorus, Possesses Exceptionally Large, Fragmented Mitogenome

In contrast to the highly conserved mitogenomic structure and organisation in most animals (including rotifers), the two previously sequenced monogonont rotifer mitogenomes were fragmented into two chromosomes similar in size, each of which possessed one major non-coding region (mNCR) of about 4–5 Kbp. To further explore this phenomenon, we have sequenced and analysed the mitogenome of one of the most studied monogonont rotifers, Brachionus calyciflorus. It is also composed of two circular chromosomes, but the chromosome-I is extremely large (27 535 bp; 3 mNCRs), whereas the chromosome-II is relatively small (9 833 bp; 1 mNCR). With the total size of 37 368 bp, it is one of the largest metazoan mitogenomes ever reported. In comparison to other monogononts, gene distribution between the two chromosomes and gene order are different and the number of mNCRs is doubled. Atp8 was not found (common in rotifers), and Cytb was present in two copies (the first report in rotifers). A high number (99) of SNPs indicates fast evolution of the Cytb-1 copy. The four mNCRs (5.3–5.5 Kb) were relatively similar. Publication of this sequence shall contribute to the understanding of the evolutionary history of the unique mitogenomic organisation in this group of rotifers.

The mitochondrial genome of the egg-laying flatworm Aglaiogyrodactylus forficulatus (Platyhelminthes: Monogenoidea)

BACKGROUND

The rather species-poor oviparous gyrodactylids are restricted to South America. It was suggested that they have a basal position within the otherwise viviparous Gyrodactylidae. Accordingly, it was proposed that the species-rich viviparous gyrodactylids diversified and dispersed from there.

METHODS

The mitochondrial genome of Aglaiogyrodactylus forficulatus was bioinformatically assembled from next-generation illumina MiSeq sequencing reads, annotated, and compared to previously published mitochondrial genomes of other monogenoidean flatworm species.

RESULTS

The mitochondrial genome of A. forficulatus consists of 14,371 bp with an average A + T content of 75.12 %. All expected 12 protein coding, 22 tRNA, and 2 rRNA genes were identified. Furthermore, there were two repetitive non-coding regions essentially consisting of 88 bp and 233 bp repeats, respectively. Maximum Likelihood analyses placed the mitochondrial genome of A. forficulatus in a well-supported clade together with the viviparous Gyrodactylidae species. The gene order differs in comparison to that of other monogenoidean species, with rearrangements mainly affecting tRNA genes. In comparison to Paragyrodactylus variegatus, four gene order rearrangements, i.e. three transpositions and one complex tandem-duplication-random-loss event, were detected.

CONCLUSION

Mitochondrial genome sequence analyses support a basal position of the oviparous A. forficulatus within Gyrodactylidae, and a sister group relationship of the oviparous and viviparous forms.

Mitochondrial genome evolution and tRNA truncation in Acariformes mites: new evidence from eriophyoid mites

The subclass Acari (mites and ticks) comprises two super-orders: Acariformes and Parasitiformes. Most species of the Parasitiformes known retained the ancestral pattern of mitochondrial (mt) gene arrangement of arthropods, and their mt tRNAs have the typical cloverleaf structure. All of the species of the Acariformes known, however, have rearranged mt genomes and truncated mt tRNAs. We sequenced the mt genomes of two species of Eriophyoidea: Phyllocoptes taishanensis and Epitrimerus sabinae. The mt genomes of P. taishanensis and E. sabinae are 13,475 bp and 13,531 bp, respectively, are circular and contain the 37 genes typical of animals; most mt tRNAs are highly truncated in both mites. On the other hand, these two eriophyoid mites have the least rearranged mt genomes seen in the Acariformes. Comparison between eriophyoid mites and other Aacariformes mites showed that: 1) the most recent common ancestor of Acariformes mites retained the ancestral pattern of mt gene arrangement of arthropods with slight modifications; 2) truncation of tRNAs for cysteine, phenylalanine and histidine occurred once in the most recent common ancestor of Acariformes mites whereas truncation of other tRNAs occurred multiple times; and 3) the placement of eriophyoid mites in the order Trombidiformes needs to be reviewed.

Revision of Leptorhynchoides thecatus (Acanthocephala: Illiosentidae), with morphometric analysis and description of six new species

Six new species of Leptorhynchoides from the southeastern United States are described. These new species were once part of the Leptorhynchoides thecatus complex of species that was previously recognized on the basis of DNA sequence data. Multivariate morphometric analysis including discriminant function analysis and decision tree analysis indicated that each of the species is morphologically distinct. Both analyses classified more than 90% of specimens correctly and most misclassifications occurred between members of 2 pairs of species that are morphologically similar. The most discriminating continuous characters were: trunk length, number of longitudinal rows of hooks, length of the longest hook, and testes width. Hook asymmetry and missing hooks on the proboscis were also important taxonomic characters. The discriminant function and the decision tree generated from the data were used to classify new specimens, yielding a 96% and 84% correct classification rate, respectively. The new taxonomic designations account for much of the previously recognized variability in host use, habitat use, and development as determined by survey data. With the addition of these 6 new taxa, 10 species currently are recognized within the genus.

Classification of the acanthocephala

In 1985, Amin presented a new system for the classification of the Acanthocephala in Crompton and Nickol's (1985) book 'Biology of the Acanthocephala' and recognized the concepts of Meyer (1931, 1932, 1933) and Van Cleave (1936, 1941, 1947, 1948, 1949, 1951, 1952). This system became the standard for the taxonomy of this group and remains so to date. Many changes have taken place and many new genera and species, as well as higher taxa, have been described since. An updated version of the 1985 scheme incorporating new concepts in molecular taxonomy, gene sequencing and phylogenetic studies is presented. The hierarchy has undergone a total face lift with Amin's (1987) addition of a new class, Polyacanthocephala (and a new order and family) to remove inconsistencies in the class Palaeacanthocephala. Amin and Ha (2008) added a third order (and a new family) to the Palaeacanthocephala, Heteramorphida, which combines features from the palaeacanthocephalan families Polymorphidae and Heteracanthocephalidae. Other families and subfamilies have been added but some have been eliminated, e.g. the three subfamilies of Arythmacanthidae: Arhythmacanthinae Yamaguti, 1935; Neoacanthocephaloidinae Golvan, 1960; and Paracanthocephaloidinae Golvan, 1969. Amin (1985) listed 22 families, 122 genera and 903 species (4, 4 and 14 families; 13, 28 and 81 genera; 167, 167 and 569 species in Archiacanthocephala, Eoacanthocephala and Palaeacanthocephala, respectively). The number of taxa listed in the present treatment is 26 families (18% increase), 157 genera (29%), and 1298 species (44%) (4, 4 and 16; 18, 29 and 106; 189, 255 and 845, in the same order), which also includes 1 family, 1 genus and 4 species in the class Polyacanthocephala Amin, 1987, and 3 genera and 5 species in the fossil family Zhijinitidae.

A comprehensive analysis of bilaterian mitochondrial genomes and phylogeny

About 2800 mitochondrial genomes of Metazoa are present in NCBI RefSeq today, two thirds belonging to vertebrates. Metazoan phylogeny was recently challenged by large scale EST approaches (phylogenomics), stabilizing classical nodes while simultaneously supporting new sister group hypotheses. The use of mitochondrial data in deep phylogeny analyses was often criticized because of high substitution rates on nucleotides, large differences in amino acid substitution rate between taxa, and biases in nucleotide frequencies. Nevertheless, mitochondrial genome data might still be promising as it allows for a larger taxon sampling, while presenting a smaller amount of sequence information. We present the most comprehensive analysis of bilaterian relationships based on mitochondrial genome data. The analyzed data set comprises more than 650 mitochondrial genomes that have been chosen to represent a profound sample of the phylogenetic as well as sequence diversity. The results are based on high quality amino acid alignments obtained from a complete reannotation of the mitogenomic sequences from NCBI RefSeq database. However, the results failed to give support for many otherwise undisputed high-ranking taxa, like Mollusca, Hexapoda, Arthropoda, and suffer from extreme long branches of Nematoda, Platyhelminthes, and some other taxa. In order to identify the sources of misleading phylogenetic signals, we discuss several problems associated with mitochondrial genome data sets, e.g. the nucleotide and amino acid landscapes and a strong correlation of gene rearrangements with long branches.

Phylogenetic analyses of endoparasitic Acanthocephala based on mitochondrial genomes suggest secondary loss of sensory organs

The metazoan taxon Syndermata (Monogononta, Bdelloidea, Seisonidea, Acanthocephala) comprises species with vastly different lifestyles. The focus of this study is on the phylogeny within the syndermatan subtaxon Acanthocephala (thorny-headed worms, obligate endoparasites). In order to investigate the controversially discussed phylogenetic relationships of acanthocephalan subtaxa we have sequenced the mitochondrial (mt) genomes of Echinorhynchus truttae (Palaeacanthocephala), Paratenuisentis ambiguus (Eoacanthocephala), Macracanthorhynchus hirudinaceus (Archiacanthocephala), and Philodina citrina (Bdelloidea). In doing so, we present the largest molecular phylogenetic dataset so far for this question comprising all major subgroups of Acanthocephala. Alongside with publicly available mt genome data of four additional syndermatans as well as 18 other lophotrochozoan (spiralian) taxa and one outgroup representative, the derived protein-coding sequences were used for Maximum Likelihood as well as Bayesian phylogenetic analyses. We achieved entirely congruent results, whereupon monophyletic Archiacanthocephala represent the sister taxon of a clade comprising Eoacanthocephala and monophyletic Palaeacanthocephala (Echinorhynchida). This topology suggests the secondary loss of lateral sensory organs (sensory pores) within Palaeacanthocephala and is further in line with the emergence of apical sensory organs in the stem lineage of Archiacanthocephala.

A mitogenomic re-evaluation of the bdelloid phylogeny and relationships among the Syndermata

Molecular and morphological data regarding the relationships among the three classes of Rotifera (Bdelloidea, Seisonidea, and Monogononta) and the phylum Acanthocephala are inconclusive. In particular, Bdelloidea lacks molecular-based phylogenetic appraisal. I obtained coding sequences from the mitochondrial genomes of twelve bdelloids and two monogononts to explore the molecular phylogeny of Bdelloidea and provide insight into the relationships among lineages of Syndermata (Rotifera + Acanthocephala). With additional sequences taken from previously published mitochondrial genomes, the total dataset included nine species of bdelloids, three species of monogononts, and two species of acanthocephalans. A supermatrix of these 10-12 mitochondrial proteins consistently recovered a bdelloid phylogeny that questions the validity of a generally accepted classification scheme despite different methods of inference and various parameter adjustments. Specifically, results showed that neither the family Philodinidae nor the order Philodinida are monophyletic as currently defined. The application of a similar analytical strategy to assess syndermate relationships recovered either a tree with Bdelloidea and Monogononta as sister taxa (Eurotatoria) or Bdelloidea and Acanthocephala as sister taxa (Lemniscea). Both outgroup choice and method of inference affected the topological outcome emphasizing the need for sequences from more closely related outgroups and more sophisticated methods of analysis that can account for the complexity of the data.

The complete mitochondrial genome sequence of Oncicola luehei (Acanthocephala: Archiacanthocephala) and its phylogenetic position within Syndermata

In the present study, we determined the complete mitochondrial genome sequence of Oncicola luehei (14,281bp), the first archiacanthocephalan representative and the second complete sequence from the phylum Acanthocephala. The complete genome contains 36 genes including 12 protein coding genes, 22 transfer RNA (tRNA) genes and 2 ribosomal RNA genes (rrnL and rrnS) as reported for other syndermatan species. All genes are encoded on the same strand. The overall nucleotide composition of O. luehei mtDNA is 37.7% T, 29.6% G, 22.5% A, and 10.2% C. The overall A+T content (60.2%) is much lower, compared to other syndermatan species reported so far, due to the high frequency (18.3%) of valine encoded by GTN in its protein-coding genes. Results from phylogenetic analyses of amino acid sequences for 10 protein-coding genes from 41 representatives of major metazoan groups including O. luehei supported monophyly of the phylum Acanthocephala and of the clade Syndermata (Acanthocephala+Rotifera), and the paraphyly of the clade Eurotatoria (classes Bdelloidea+Monogononta from phylum Rotifera). Considering the position of the acanthocephalan species within Syndermata, it is inferred that obligatory parasitism characteristic of acanthocephalans was acquired after the common ancestor of acanthocephalans diverged from its sister group, Bdelloidea. Additional comparison of complete mtDNA sequences from unsampled acanthocephalan lineages, especially classes Polyacanthocephala and Eoacanthocephala, is required to test if mtDNA provides reliable information for the evolutionary relationships and pattern of life history diversification found in the syndermatan groups.

Repeated parallel evolution of minimal rRNAs revealed from detailed comparative analysis

The concept of a minimal ribosomal RNA-containing ribosome, a structure with a minimal set of elements capable of providing protein biosynthesis, is essential for understanding this fundamental cellular process. Nematodes and trypanosomes have minimal mitochondrial rRNAs and detailed reconstructions of their secondary structures indicate that certain conserved helices have been lost in these taxa. In contrast, several recent studies on acariform mites have argued that minimal rRNAs may evolve via shortening of secondary structure elements but not the loss of these elements as shown for trypanosomes and nematodes. Based on extensive structural analysis of chelicerate arthropods, we demonstrate that extremely short rRNAs of acariform mites share certain structural modifications with nematodes and trypanosomes: loss of helices of the GTPase region and divergence in the evolutionarily conserved connecting loop between helices H1648 and H1764 of the large subunit rRNA. These highly concerted parallel modifications indicate that minimal rRNAs were generated under the strong selection that favored or tolerated reductions of helices in particular locations while maintaining the functionality of the rRNA molecules throughout evolution. We also discuss potential evolution of minimal rRNAs and atypical transfer RNAs.

The mitochondrial genomes of two nemerteans, Cephalothrix sp. (Nemertea: Palaeonemertea) and Paranemertes cf. peregrina (Nemertea: Hoplonemertea)

The mitochondrial genome sequences were determined for two species of nemerteans, Cephalothrix sp. (15,800 bp sequenced, near-complete) and Paranemertes cf. peregrina (14,558 bp, complete). As seen in most metazoans, the genomes encode 13 protein, 2 ribosomal RNA and 22 transfer RNA genes. The nucleotide composition is strongly biased toward A and T, as is typical for metazoan mtDNAs. There is also a significant strand skew in the distribution of these nucleotides, with the coding strand being richer in T than A and in G than C. All genes are transcribed in the same direction except for trnP and trnT, which is consistent with that reported for Cephalothrix hongkongiensis and Lineus viridis. Gene arrangement of Cephalothrix sp. is identical to that of C. hongkongiensis, while in P. cf. peregrina it is similar to L. viridis, but differs significantly from the three Cephalothrix species in the position of four protein-coding genes and seven tRNAs. Some protein-coding genes have 3' end stem-loop structures, which may allow mRNA processing without flanking tRNAs. The major non-coding regions observed in the two genomes with potential to form stem-loop structures may be involved in the initiation of replication or transcription. The average Ka/Ks values, varying from 0.12 to 0.89, are markedly different among the 13 mitochondrial protein-coding genes, suggesting that there may exist different selective pressure among mitochondrial genes of nemerteans.

Characterization of a mitochondrial ORF from the gender-associated mtDNAs of Mytilus spp. (Bivalvia: Mytilidae): identification of the "missing" ATPase 8 gene

Bivalve species are characterized by extraordinary variability in terms of mitochondrial (mt) genome size, gene arrangement and tRNA gene number. Many species are thought to lack the mitochondrial protein-coding gene atp8. Of these species, the Mytilidae appears to be the only known taxon with doubly uniparental inheritance of mtDNA that does not possess the atp8 gene. This raises the question as to whether mytilids have completely lost the ATP8 protein, whether the gene has been transferred to the nucleus or whether they possess a highly modified version of the gene/protein that has led to its lack of annotation. In the present study, we re-investigated all complete (or nearly complete) F and M mytilid mt genomes previously sequenced for the presence of conserved open reading frames (ORFs) that might code for ATP8 and/or have other functional importance in these bivalves. We also revised the annotations of all available complete mitochondrial genomes of bivalves and nematodes that are thought to lack atp8 in an attempt to detect it. Our results indicate that a novel mytilid ORF of significant length (i.e., the ORF is >85 amino acids in length), with complete start and stop codons, is a candidate for the atp8 gene: (1) it possesses a pattern of evolution expected for a protein-coding gene evolving under purifying selection (i.e., the 3rd>1st>2nd codon pattern of evolution), (2) it is actively transcribed in Mytilus species, (3) it has one predicted transmembrane helix (as do other metazoan ATP8 proteins), (4) it has conserved functional motifs and (5), comparisons of its amino acid sequence with ATP8 sequences of other molluscan or bivalve species reveal similar hydropathy profiles. Furthermore, our revised annotations also confirmed the mt presence of atp8 in almost all bivalve species and in one nematode species. Our results thus support recognizing the presence of ATPase 8 in most bivalves mt genomes (if not all) rather than the continued characterization of these genomes as lacking this gene.

Eurotatorian paraphyly: Revisiting phylogenetic relationships based on the complete mitochondrial genome sequence of Rotaria rotatoria (Bdelloidea: Rotifera: Syndermata)

BACKGROUND

The Syndermata (Rotifera+Acanthocephala) is one of the best model systems for studying the evolutionary origins and persistence of different life styles because it contains a series of lineage-specific life histories: Monogononta (cyclic parthenogenetic and free-living), Bdelloidea (entirely parthenogenetic and mostly benthic dweller), Seisonidea (exclusively bisexual and epizoic or ectoparasitic), and Acanthocephala (sexual and obligatory endoparasitic). Providing phylogenetic resolution to the question of Eurotatoria (Monogononta and Bdelloidea) monophyly versus paraphyly is a key factor for better understanding the evolution of different life styles, yet this matter is not clearly resolved. In this study, we revisited this issue based on comparative analysis of complete mitochondrial genome information for major groups of the Syndermata.

RESULTS

We determined the first complete mitochondrial genome sequences (15,319 bp) of a bdelloid rotifer, Rotaria rotatoria. In order to examine the validity of Eurotatoria (Monogononta and Bdelloidea) monophyly/paraphyly, we performed phylogenetic analysis of amino acid sequences for eleven protein-coding genes sampled from a wide variety of bilaterian representatives. The resulting mitochondrial genome trees, inferred using different algorithms, consistently failed to recover Monogononta and Bdelloidea as monophyletic, but instead identified them as a paraphyletic assemblage. Bdelloidea (as represented by R. rotatoria) shares most common ancestry with Acanthocephala (as represented by L. thecatus) rather than with monogonont B. plicatilis, the other representative of Eurotatoria.

CONCLUSION

Comparisons of inferred amino acid sequence and gene arrangement patterns with those of other metazoan mtDNAs (including those of acanthocephalan L. thecatus and monogonont B. plicatilis) support the hypothesis that Bdelloidea shares most common ancestry with Acanthocephala rather than with Monogononta. From this finding, we suggest that the obligatory asexuality of bdelloideans may have secondarily derived from some other preexisting condition in earlier lineage of rotifers. Providing a more complete assessment of phylogenetic relationships and inferring patterns of evolution of different types of life styles among Syndermata awaits comparisons requiring mitochondrial genome sequencing of Seisonidea.

Evolution of the mitochondrial genome of Metazoa as exemplified by comparison of congeneric species

The mitochondrial genome (mtDNA) of Metazoa is a good model system for evolutionary genomic studies and the availability of more than 1000 sequences provides an almost unique opportunity to decode the mechanisms of genome evolution over a large phylogenetic range. In this paper, we review several structural features of the metazoan mtDNA, such as gene content, genome size, genome architecture and the new parameter of gene strand asymmetry in a phylogenetic framework. The data reviewed here show that: (1) the plasticity of Metazoa mtDNA is higher than previously thought and mainly due to variation in number and location of tRNA genes; (2) an exceptional trend towards stabilization of genomic features occurred in deuterostomes and was exacerbated in vertebrates, where gene content, genome architecture and gene strand asymmetry are almost invariant. Only tunicates exhibit a very high degree of genome variability comparable to that found outside deuterostomes. In order to analyse the genomic evolutionary process at short evolutionary distances, we have also compared mtDNAs of species belonging to the same genus: the variability observed in congeneric species significantly recapitulates the evolutionary dynamics observed at higher taxonomic ranks, especially for taxa showing high levels of genome plasticity and/or fast nucleotide substitution rates. Thus, the analysis of congeneric species promises to be a valuable approach for the assessment of the mtDNA evolutionary trend in poorly or not yet sampled metazoan groups.

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.

Complete nucleotide sequences of mitochondrial genomes of two solitary entoprocts, Loxocorone allax and Loxosomella aloxiata: implications for lophotrochozoan phylogeny

The complete nucleotide sequences of the mitochondrial (mt) genomes of the entoprocts Loxocorone allax and Loxosomella aloxiata were determined. Both species carry the typical gene set of metazoan mt genomes and have similar organizations of their mt genes. However, they show differences in the positions of two tRNA(Leu) genes. Additionally, the tRNA(Val) gene, and half of the long non-coding region, is duplicated and inverted in the Loxos. aloxiata mt genome. The initiation codon of the Loxos. aloxiata cytochrome oxidase subunit I gene is expected to be ACG rather than AUG. The mt gene organizations in these two entoproct species most closely resemble those of mollusks such as Katharina tunicata and Octopus vulgaris, which have the most evolutionarily conserved mt gene organization reported to date in mollusks. Analyses of the mt gene organization in the lophotrochozoan phyla (Annelida, Brachiopoda, Echiura, Entoprocta, Mollusca, Nemertea, and Phoronida) suggested a close phylogenetic relationship between Brachiopoda, Annelida, and Echiura. However, Phoronida was excluded from this grouping. Molecular phylogenetic analyses based on the sequences of mt protein-coding genes suggested a possible close relationship between Entoprocta and Phoronida, and a close relationship among Brachiopoda, Annelida, and Echiura.

Cryptic speciation and patterns of phenotypic variation of a highly variable acanthocephalan parasite

An investigation of a parasite species that is broadly host- and habitat-specific and exhibits alternative transmission strategies was undertaken to examine intraspecific variability and if it can be attributed to cryptic speciation or environmentally induced plasticity. Specimens of an acanthocephalan parasite, Leptorhynchoides thecatus, collected throughout North America were analysed phylogenetically using sequences of the cytochrome oxidase I gene and the internal transcribed spacer region. Variation in host use, habitat use, and transmission were examined in a phylogenetic context to determine if they were more likely phylogenetically based or due to environmental influences. Results indicated that most of the variation detected can be explained by the presence of cryptic species. The majority of these species have narrow host and microhabitat specificities although one species, which also may comprise a complex of species, exhibits broad host and habitat specificity. Alternate transmission pathways only occurred in two of the cryptic species and correlate with host use patterns. Taxa that mature in piscivorous piscine hosts use a paratenic fish host to bridge the trophic gap between their amphipod intermediate host and piscivorous definitive host. One potential example of environmentally induced variation was identified in three populations of these parasites, which differ on their abilities to infect different host species.

The mitochondrial genome of the screamer louse Bothriometopus (phthiraptera: ischnocera): effects of extensive gene rearrangements on the evolution of the genome

Mitochondrial (mt) genome rearrangement has generally been studied with respect to the phenomenon itself, focusing on their phylogenetic distribution and causal mechanisms. Rearrangements have additional significance through effects on substitution, transcription, and mRNA processing. Lice are an ideal group in which to study the interactions between rearrangements and these factors due to the heightened rearrangement rate within this group. The entire mt genome of the screamer louse Bothriometopus was sequenced and compared to previously sequenced louse genomes. The mt genome is 15,564 bp, circular, and all genes are encoded on the same strand. The gene arrangement differs radically from both other louse species and the ancestral insect. Nucleotide composition is A+T biased, but there is no skew which may be due to reversal of replication direction or a transcriptional effect. Bothriometopus has both tRNA duplication and concerted evolution which has not been observed previously. Eleven of the 13 protein-coding genes have 3' end stem-loop structures which may allow mRNA processing without flanking tRNAs and so facilitate gene rearrangements. There are five candidate control regions capable of forming stem-loop structures. Two are structurally more similar to the control regions of other insect species than those of other lice. Analyses of Bothriometopus demonstrate that louse mt genomes, in addition to being extensively rearranged, differ significantly from most insect species in nucleotide composition biases, tRNA evolution, protein-coding gene structures and putative signaling sites such as the control region. These may be either a cause or a consequence of gene rearrangements.

Visualizing differences in phylogenetic information content of alignments and distinction of three classes of long-branch effects

BACKGROUND

Published molecular phylogenies are usually based on data whose quality has not been explored prior to tree inference. This leads to errors because trees obtained with conventional methods suppress conflicting evidence, and because support values may be high even if there is no distinct phylogenetic signal. Tools that allow an a priori examination of data quality are rarely applied.

RESULTS

Using data from published molecular analyses on the phylogeny of crustaceans it is shown that tree topologies and popular support values do not show existing differences in data quality. To visualize variations in signal distinctness, we use network analyses based on split decomposition and split support spectra. Both methods show the same differences in data quality and the same clade-supporting patterns. Both methods are useful to discover long-branch effects. We discern three classes of long branch effects. Class I effects consist of attraction of terminal taxa caused by symplesiomorphies, which results in a false monophyly of paraphyletic groups. Addition of carefully selected taxa can fix this effect. Class II effects are caused by drastic signal erosion. Long branches affected by this phenomenon usually slip down the tree to form false clades that in reality are polyphyletic. To recover the correct phylogeny, more conservative genes must be used. Class III effects consist of attraction due to accumulated chance similarities or convergent character states. This sort of noise can be reduced by selecting less variable portions of the data set, avoiding biases, and adding slower genes.

CONCLUSION

To increase confidence in molecular phylogenies an exploratory analysis of the signal to noise ratio can be conducted with split decomposition methods. If long-branch effects are detected, it is necessary to discern between three classes of effects to find the best approach for an improvement of the raw data.

The mitochondrial genome of Priapulus caudatus Lamarck (Priapulida: Priapulidae)

We sequenced and annotated the complete mitochondrial (mt) genome of the priapulid Priapulus caudatus in order to provide a source of phylogenetic characters including an assessment of gene order arrangement. The genome was 14,919 bp in its entirety with few, short non-coding regions. A number of protein-coding and tRNA genes overlapped, making the genome relatively compact. The gene order was: cox1, cox2, trnK, trnD, atp8, atp6, cox3, trnG, nad3, trnA, trnR, trnN, rrnS, trnV, rrnL, trnL(yaa), trnL(nag), nad1, -trnS(nga), -cob, -nad6, trnP, -trnT, nad4L, nad4, trnH, nad5, trnF, -trnE, -trnS(nct), trnI, -trnQ, trnM, nad2, trnW, -trnC, -trnY; where '-' indicates genes transcribed on the opposite strand. The gene order, although unique amongst Metazoa, shared the greatest number of gene boundaries and the longest contiguous fragments with the chelicerate Limulus polyphemus. The mt genomes of these taxa differed only by a single inversion of 18 contiguous genes bounded by rrnS and trnS(nct). Other arthropods and nematodes shared fewer gene boundaries but considerably more than the most similar non-ecdysozoan.

A protein extension to shorten RNA: elongated elongation factor-Tu recognizes the D-arm of T-armless tRNAs in nematode mitochondria

Nematode mitochondria possess extremely truncated tRNAs. Of 22 tRNAs, 20 lack the entire T-arm. The T-arm is necessary for the binding of canonical tRNAs and EF (elongation factor)-Tu (thermo-unstable). The nematode mitochondrial translation system employs two different EF-Tu factors named EF-Tu1 and EF-Tu2. Our previous study showed that nematode Caenorhabditis elegans EF-Tu1 binds specifically to T-armless tRNA. C. elegans EF-Tu1 has a 57-amino acid C-terminal extension that is absent from canonical EF-Tu, and the T-arm-binding residues of canonical EF-Tu are not conserved. In this study, the recognition mechanism of T-armless tRNA by EF-Tu1 was investigated. Both modification interference assays and primer extension analysis of cross-linked ternary complexes revealed that EF-Tu1 interacts not only with the tRNA acceptor stem but also with the D-arm. This is the first example of an EF-Tu recognizing the D-arm of a tRNA. The binding activity of EF-Tu1 was impaired by deletion of only 14 residues from the C-terminus, indicating that the C-terminus of EF-Tu1 is required for its binding to T-armless tRNA. These results suggest that C. elegans EF-Tu1 recognizes the D-arm instead of the T-arm by a mechanism involving its C-terminal region. This study sheds light on the co-evolution of RNA and RNA-binding proteins in nematode mitochondria.

An evolutionary 'intermediate state' of mitochondrial translation systems found in Trichinella species of parasitic nematodes: co-evolution of tRNA and EF-Tu

EF-Tu delivers aminoacyl-tRNAs to ribosomes in the translation system. However, unusual truncations found in some animal mitochondrial tRNAs seem to prevent recognition by a canonical EF-Tu. We showed previously that the chromadorean nematode has two distinct EF-Tus, one of which (EF-Tu1) binds only to T-armless aminoacyl-tRNAs and the other (EF-Tu2) binds to D-armless Ser-tRNAs. Neither of the EF-Tus can bind to canonical cloverleaf tRNAs. In this study, by analyzing the translation system of enoplean nematode Trichinella species, we address how EF-Tus and tRNAs have evolved from the canonical structures toward those of the chromadorean translation system. Trichinella mitochondria possess three types of tRNAs: cloverleaf tRNAs, which do not exist in chromadorean nematode mitochondria; T-armless tRNAs; and D-armless tRNAs. We found two mitochondrial EF-Tu species, EF-Tu1 and EF-Tu2, in Trichinella britovi. T.britovi EF-Tu2 could bind to only D-armless Ser-tRNA, as Caenorhabditis elegans EF-Tu2 does. In contrast to the case of C.elegans EF-Tu1, however, T.britovi EF-Tu1 bound to all three types of tRNA present in Trichinella mitochondria. These results suggest that Trichinella mitochondrial translation system, and particularly the tRNA-binding specificity of EF-Tu1, could be an intermediate state between the canonical system and the chromadorean nematode mitochondrial system.

GenDecoder: genetic code prediction for metazoan mitochondria

Although the majority of the organisms use the same genetic code to translate DNA, several variants have been described in a wide range of organisms, both in nuclear and organellar systems, many of them corresponding to metazoan mitochondria. These variants are usually found by comparative sequence analyses, either conducted manually or with the computer. Basically, when a particular codon in a query-species is linked to positions for which a specific amino acid is consistently found in other species, then that particular codon is expected to translate as that specific amino acid. Importantly, and despite the simplicity of this approach, there are no available tools to help predicting the genetic code of an organism. We present here GenDecoder, a web server for the characterization and prediction of mitochondrial genetic codes in animals. The analysis of automatic predictions for 681 metazoans aimed us to study some properties of the comparative method, in particular, the relationship among sequence conservation, taxonomic sampling and reliability of assignments. Overall, the method is highly precise (99%), although highly divergent organisms such as platyhelminths are more problematic. The GenDecoder web server is freely available from http://darwin.uvigo.es/software/gendecoder.html.