Mitochondrial genomes of Galathealinum, Helobdella, and Platynereis: sequence and gene arrangement comparisons indicate that Pogonophora is not a phylum and Annelida and Arthropoda are not sister taxa

Cost-efficient PCR based DNA barcoding of marine invertebrate specimens with NovaSeq amplicon sequencing

BACKGROUND

The marine environment harbors high biodiversity; however, it is poorly understood. Nucleotide sequence data of all marine organisms should be accumulated before natural and/or anthropogenic environmental changes jeopardize the marine environment. In this study, we report a cost-effective and easy DNA barcoding method. This method can be readily adopted without using library preparation kits. It includes multiplex PCR of short targets, indexing PCR, and outsourcing to a sequencing service using the NovaSeq system.

METHODS AND RESULTS

We targeted four mitochondrial genes [cytochrome c oxidase subunit I (COI), COIII, 16S rRNA (16S), and 12S rRNA (12S)] and three nuclear genes [18S rRNA (18S), 28S rRNA (28S), internal transcribed spacer 2 (ITS2)] in 95 marine invertebrate specimens, which were primarily annelids. The primers, including adapters and indices for NovaSeq sequencing, were newly designed. Two PCR runs were conducted. The 1st PCR amplified specific loci with universal primers and the 2nd added sequencing adapters and indices to the 1st PCR products. The gene sequences obtained from the FASTQ files were subjected to BLAST search and phylogenetic analyses. One run using 95 specimens yielded sequences averaging 2816 bp per specimen for a total length of six loci. Nuclear genes were more successfully assembled compared with mitochondrial genes. A weak but significantly negative correlation was observed between the average length of each locus and success rate of the assembly. Some of the sequences were almost identical to the sequences obtained from specimens collected far from Japan, indicating the presence of potentially invasive species identified for the first time.

CONCLUSIONS

We obtained gene sequences efficiently using next-generation sequencing rather than Sanger sequencing. Although this method requires further optimization to increase the success rate for some loci, it is used as a first step to select specimens for further analyses by determining the specific loci of the targets.

The Mitochondrial Genomes of Siboglinum plumosum and Oligobrachia dogieli (Annelida: Siboglinidae) and Their Phylogenetic Analysis

Frenulates are a group of sedentary Annelida within the family Siboglinidae that inhabit the ocean floor and present a unique challenge for comprehensive molecular and phylogenetic investigations. In this study, we focused on the frenulates, specifically assembling the mitochondrial genomes of Siboglinum plumosum and Oligobrachia dogieli. The phylogenetic reconstruction placed S. plumosum as a sister taxon to S. ekmani, and O. dogieli as a sister taxon to S. fiordicum, supporting the non-monophyletic nature of the genus Siboglinum. Overall, this study supports the phylogeny of the family Siboglinidae while highlighting the need for additional molecular data within frenulates.

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.

Mitochondrial Genome Evolution in Annelida-A Systematic Study on Conservative and Variable Gene Orders and the Factors Influencing its Evolution

The mitochondrial genomes of Bilateria are relatively conserved in their protein-coding, rRNA, and tRNA gene complement, but the order of these genes can range from very conserved to very variable depending on the taxon. The supposedly conserved gene order of Annelida has been used to support the placement of some taxa within Annelida. Recently, authors have cast doubts on the conserved nature of the annelid gene order. Various factors may influence gene order variability including, among others, increased substitution rates, base composition differences, structure of noncoding regions, parasitism, living in extreme habitats, short generation times, and biomineralization. However, these analyses were neither done systematically nor based on well-established reference trees. Several focused on only a few of these factors and biological factors were usually explored ad-hoc without rigorous testing or correlation analyses. Herein, we investigated the variability and evolution of the annelid gene order and the factors that potentially influenced its evolution, using a comprehensive and systematic approach. The analyses were based on 170 genomes, including 33 previously unrepresented species. Our analyses included 706 different molecular properties, 20 life-history and ecological traits, and a reference tree corresponding to recent improvements concerning the annelid tree. The results showed that the gene order with and without tRNAs is generally conserved. However, individual taxa exhibit higher degrees of variability. None of the analyzed life-history and ecological traits explained the observed variability across mitochondrial gene orders. In contrast, the combination and interaction of the best-predicting factors for substitution rate and base composition explained up to 30% of the observed variability. Accordingly, correlation analyses of different molecular properties of the mitochondrial genomes showed an intricate network of direct and indirect correlations between the different molecular factors. Hence, gene order evolution seems to be driven by molecular evolutionary aspects rather than by life history or ecology. On the other hand, variability of the gene order does not predict if a taxon is difficult to place in molecular phylogenetic reconstructions using sequence data or not. We also discuss the molecular properties of annelid mitochondrial genomes considering canonical views on gene evolution and potential reasons why the canonical views do not always fit to the observed patterns without making some adjustments. [Annelida; compositional biases; ecology; gene order; life history; macroevolution; mitochondrial genomes; substitution rates.].

Mitochondrial genome of Leocrates chinensis Kinberg, 1866 (Annelida: Hesionidae)

We report the complete mitochondrial genome of Leocrates chinensis Kinberg, 1866 - the type species of the genus. It is 15061 bp long, and contains 13 protein-coding genes (PCGs), 22 tRNA genes (tRNAs), and 2 rRNA genes (rRNAs), and 1 putative control region. Phylogenetic analysis indicated that L. chinensis was placed as sister to Sirsoe methanicola (BS = 100) of the same family Hesionidae.

Two New Uromunna Species (Isopoda: Asellota: Munnidae) from the Korean Peninsula and Their Phylogenetic Position within Munnoid Groups

Two new Uromunna species, Uromunna mundongensis sp. n. and Uromunna jejuensis sp. n., are described from the Korean Peninsula, representing the first record of the genus in Korean waters. Genetic analyses using two mitochondrial (COI, Cytb) and one ribosomal (18S rRNA) gene allowed us to confirm high pair-wise interspecific differences with the establishment of reliable barcoding gaps of COI (19.9%) and 18S (0.4%) between the new species. Finally, the phylogenetic relationship between the Uromunna species as well as the position of the genus within the munnoid groups are reconstructed using 18S.

A new eyeless species of Nereis (Annelida, Nereididae) from deep-sea sediments of the northern South China Sea

A variety of nereidid species have been reported from the South China Sea, although little is known about the deep-sea species in this area. Recently, two specimens belonging to a novel nereidid polychaete were collected from a sedimentary habitat during an environmental survey to a deep-sea basin where cold seeps occur. This new species, Nereistricirrata sp. nov., is described herein, based on morphological and molecular analyses. The most noteworthy feature is the absence of eyes on the prostomium; it can be distinguished from other eyeless Nereis species by the arrangement of conical paragnaths on the pharynx, the nature of homogomph falcigers and the shape of notopodial lobes in posterior chaetigers. The reconstructed phylogenetic tree, using concatenated sequences of mtCOI, 16S, and 18S rRNA, showed that all Nereis species included in this study form a monophyletic clade with full support. The mtCOI-based interspecific comparisons revealed a high genetic divergence (23.1%-37.3% K2P) from four-eyed Nereis species with the available sequences. This is the first record of an eyeless Nereis species in the South China Sea.

Sequencing and phylogenomics of the complete mitochondrial genome of Allodiplogaster sp. (Rhabditida: Diplogasteridae): A new gene order and its phylogenetic implications

Gene order has been utilized as a phylogenetic signal for many taxa. However, its phylogenetic performance has not been evaluated in Nematoda. As there is only one nematode mitogenome available to date, in the Diplogasteridae family, we sequenced the mitogenome of Allodiplogaster sp. and constructed a phylogeny for Nematoda using this updated mitogenome dataset. We then compared this phylogeny to one constructed using gene order. The complete mitochondrial genome of Allodiplogaster sp. was 13,953 bp in size and included 22 tRNAs, two rRNAs, and 12 protein-coding genes. To assess how Allodiplogaster sp. is related to other nematode species, we used Bayesian inference and maximum likelihood algorithms to construct phylogenetic trees of the Nematoda. We found that: 1) The target species Allodiplogaster sp. is closely related to Allodiplogaster sudhausi. The topology of the mitogenome based phylogeny was nearly identical to previous phylogenies created using 18S rRNA data, except for the placement of the Strongyloididae family. 2) The maximum likelihood tree constructed using gene order was roughly consistent with the mitogenome-based tree at the family level, but not at the species level. 3) Protein-coding genes were ordered differently in Allodiplogaster sp. versus Allodiplogaster sudhausi; this represents the first report of such a reordering in the class Chromadorea in our study. Our study confirms that gene order represents useful phylogenetic information for the Nematoda: the maximum likelihood tree based on gene order provided additional support for the nematode phylogeny constructed using molecular data.

First mitochondrial genomes of Chrysopetalidae (Annelida) from shallow-water and deep-sea chemosynthetic environments

Among Annelida, Chrysopetalidae is an ecologically and morphologically diverse group, which includes shallow-water, deep-sea, free-living, and symbiotic species. Here, the four first mitochondrial genomes of this group are presented and described. One of the free-living shallow-water species Chrysopetalum debile (Chrysopetalinae), one of the yet undescribed free-living deep-sea species Boudemos sp., and those of the two deep-sea bivalve endosymbionts Craseoschema thyasiricola and Iheyomytilidicola lauensis (Calamyzinae). An updated phylogeny of Chrysopetalidae is performed, which supports previous phylogenetic hypotheses within Chrysopetalinae and indicates a complex ecological evolution within Calamyzinae. Additionally, analyses of natural selection pressure in the four mitochondrial genomes and additional genes from the two shallow-water species Bhawania goodei and Arichlidon gathofi were performed. Relaxed selection pressure in the mitochondrion of deep-sea and symbiotic species was found, with many sites under selection identified in the COX3 gene of deep-sea species.

Molecular phylogeny of Nipponacmea (Patellogastropoda, Lottiidae) from Japan: a re-evaluation of species taxonomy and morphological diagnosis

The patellogastropod limpet genus Nipponacmea is widely distributed in Japan and adjacent East Asia. Species identification within Nipponacmea is challenging due to the high variation in shell morphology. In this study, we examined the taxonomy of this genus represented by nine nominal species from 43 localities (including type localities). Results of the molecular phylogenetic analysis revealed that: (1) N.gloriosa, the sole species in this genus inhabiting the subtidal zone, represents the most basal independent branch; (2) the remaining species are divided into two large clades with lower- and higher-apex shell profiles; and (3) the high-apex morphology was derived from the low-apex type. The terminal clades defined using the molecular data were consistent with nine morphospecies and had 100% bootstrap values, strongly supporting the conventional taxonomy of Nipponacmea. Although morphological similarities do not always reflect phylogeny, the set of morphological characters used in the current taxonomy were proven to be adequate for diagnosis. In conclusion, this study provided solid evidence to uphold the monophyly of known species of Nipponacmea in Japan and demonstrated the usefulness of morphological characters for species diagnosis.

Sequencing of the Complete Mitochondrial Genome of Pingus sinensis (Spirurina: Quimperiidae): Gene Arrangements and Phylogenetic Implications

Despite several decades of intensive research on spirurine nematodes, molecular data on some of the main lineages are still absent, which makes taxonomic classification insufficiently resolved. In the present study, we sequenced the first complete mitogenome for the family Quimperiidae, belonging to P. sinensis (Spirurina: Quimperiidae), a parasite living in the intestines of snakehead (Ophiocephalus argus). The circular mitogenome is 13,874 bp long, and it contains the standard nematode gene set: 22 transfer RNAs, 2 ribosomal RNAs and 12 protein-coding genes. There are also two long non-coding regions (NCR), in addition to only 8 other intergenic regions, ranging in size from 1 to 58 bp. To investigate its phylogenetic position and study the relationships among other available Spirurina, we performed the phylogenetic analysis using Bayesian inference and maximum likelihood approaches by concatenating the nucleotide sequences of all 36 genes on a dataset containing all available mitogenomes of the suborder Spirurina from NCBI and compared with gene order phylogenies using the MLGO program. Both supported the closer relationship of Ascaridoidea to Seuratoidea than to Spiruroidea. Pingus formed a sister-group with the Cucullanus genus. The results provide a new insights into the relationships within Spirurina.

Aln2tbl: building a mitochondrial features table from a assembly alignment in fasta format

The sequencing, annotation and analysis of complete mitochondrial genomes is an important research tool in phylogeny and evolution. Starting with the primary sequence, genes/features are generally annotated automatically to obtain preliminary annotations in the form of a feature table. Further manual curation in a graphic alignment editor is nevertheless necessary to revise annotations. As such, the automatically generated feature table is invalidated and has to be modified manually before submission to data banks. We developed aln2tbl.py, a python script that recreates a feature table from a manually refined alignment of genes mapped on the mitochondrial genome in fasta format. The feature table is populated with notes and annotations specific to mitochondrial genomes. The table can be used to create a sqn file to be submitted directly to data banks. In summary, our scripts fills one gap in the available toolbox and, combined with other software, allows the automation of the entire process, from primary sequence to annotated genome submission, even if a manual curation step is conducted in a visual sequence editor.

Molecular phylogenetic assessment of Spirobranchus kraussii-complex (Annelida: Serpulidae) from the Japanese Archipelago

Spirobranchus kraussii (Annelida: Serpulidae) was recognized as being widely distributed both in the Pacific and Atlantic Oceans. However, the sampling records far from its type locality (South Africa) have been questioned. Actually, recent molecular phylogenetic studies showed that S. kraussii contains genetically distinct species. In this study, we performed molecular phylogenetic analyses of S. cf. kraussii collected from Japan using the nucleotide sequences of a mitochondrial gene and two nuclear genes. Three lineages were recovered within Spirobranchus kraussii-complex in Japan, and one (Spirobranchus sp. 6) showed moderate genetic difference (approximately 4%) in the mitochondrial cytb gene sequence from Spirobranchus sp. 1, an undescribed sequenced species from Honshu Island, Japan. However, the nucleotide sequences of the 18S rRNA gene and ITS2 region were nearly indistinguishable. The other lineage was clearly distinct from the other previously sequenced species and is thus considered to be another distinct species of this species complex (Spirobranchus sp. 5). Although detailed morphological assessment of these lineages is necessary to define their taxonomic status, the present study provided further implications for the species diversity within the S. kraussii-complex.

Mitochondrial genomes of two Polydora (Spionidae) species provide further evidence that mitochondrial architecture in the Sedentaria (Annelida) is not conserved

Contrary to the early evidence, which indicated that the mitochondrial architecture in one of the two major annelida clades, Sedentaria, is relatively conserved, a handful of relatively recent studies found evidence that some species exhibit elevated rates of mitochondrial architecture evolution. We sequenced complete mitogenomes belonging to two congeneric shell-boring Spionidae species that cause considerable economic losses in the commercial marine mollusk aquaculture: Polydora brevipalpa and Polydora websteri. The two mitogenomes exhibited very similar architecture. In comparison to other sedentarians, they exhibited some standard features, including all genes encoded on the same strand, uncommon but not unique duplicated trnM gene, as well as a number of unique features. Their comparatively large size (17,673 bp) can be attributed to four non-coding regions larger than 500 bp. We identified an unusually large (putative) overlap of 14 bases between nad2 and cox1 genes in both species. Importantly, the two species exhibited completely rearranged gene orders in comparison to all other available mitogenomes. Along with Serpulidae and Sabellidae, Polydora is the third identified sedentarian lineage that exhibits disproportionally elevated rates of mitogenomic architecture rearrangements. Selection analyses indicate that these three lineages also exhibited relaxed purifying selection pressures.

Extensive gene rearrangements in the mitogenomes of congeneric annelid species and insights on the evolutionary history of the genus Ophryotrocha

BACKGROUND

Annelids are one the most speciose and ecologically diverse groups of metazoans. Although a significant effort has been recently invested in sequencing genomes of a wide array of metazoans, many orders and families within the phylum Annelida are still represented by a single specimen of a single species. The genus of interstitial annelids Ophryotrocha (Dorvilleidae, Errantia, Annelida) is among these neglected groups, despite its extensive use as model organism in numerous studies on the evolution of life history, physiological and ecological traits. To compensate for the paucity of genomic information in this genus, we here obtained novel complete mitochondrial genomes of six Ophryotrocha species using next generation sequencing. In addition, we investigated the evolution of the reproductive mode in the Ophryotrocha genus using a phylogeny based on two mitochondrial markers (COXI and 16S rDNA) and one nuclear fragment (Histone H3).

RESULTS

Surprisingly, gene order was not conserved among the six Ophryotrocha species investigated, and varied greatly as compared to those found in other annelid species within the class Errantia. The mitogenome phylogeny for the six Ophryotrocha species displayed a separation of gonochoric and hermaphroditic species. However, this separation was not observed in the phylogeny based on the COX1, 16S rDNA, and H3 genes. Parsimony and Bayesian ancestral trait reconstruction indicated that gonochorism was the most parsimonious ancestral reproductive mode in Ophryotrocha spp.

CONCLUSIONS

Our results highlight the remarkably high level of gene order variation among congeneric species, even in annelids. This encourages the need for additional mitogenome sequencing of annelid taxa in order to properly understand its mtDNA evolution, high biodiversity and phylogenetic relationships.

Mitogenomic architecture of the multivalent endemic black clam (Villorita cyprinoides) and its phylogenetic implications

The Indian black clam Villorita cyprinoides (Family: Cyrenidae), an extractive commercially exploited species with aquaculture importance contributing more than 70% of clam fishery in India, is endemic to the Indian peninsula. Currently, there is very sparse information, especially on the molecular data of Villorita. The present study aims to provide a comprehensive knowledge of mitogenome architecture and assess the phylogenetic status of Cyrenidae. This has resulted in reporting the first complete mitogenome of V. cyprinoides using next-generation sequencing technology. The A+T circular mitogenome was 15,880 bp long, exhibiting 13 protein-coding genes (PCGs) including ATP8 (absent in several bivalves), 22 transfer RNA, and two ribosomal RNA genes residing in the heavy strand in a clockwise orientation and a gene order akin to Corbicula fluminea. The molecular phylogeny inferred from a concatenated multi-gene sequence [14 mitochondrial (12 PCGs, rrnS and rrnL) and two nuclear genes (Histone H3, 18S rRNA)] from 47 representative species of superorder Imparidentia, clustered V. cyprinoides and Cyrenid clams to a single clade supporting the monophyly of Cyrenidae. The subsequent mitochondrial gene order analysis substantiates the close relationship of V. cyprinoides and C. fluminea, analogous to phylogenetic output. The multilocus tree topology calibrated with verified fossil data deciphered the origin and diversification of Cyrenid clams during late Triassic-early Jurassic. The data derived from this study shall contribute remarkably for further insights on cryptic species identification, molecular characterization of bivalve mitogenomes and mitochondrial evolutionary history of genus Villorita. Moreover, complete mitogenome can aid in potential marker development for assessing the genetic health of black clam populations.

Myxozoan hidden diversity: the case of Myxobolus pseudodispar Gorbunova, 1936

Myxobolus pseudodispar Gorbunova, 1936 (Myxozoa) was originally described as a parasite of common roach, Rutilus rutilus (Linnaeus), with developing stages in muscles and spores disseminated in macrophage centres of different organs and tissues. Later, this parasite was described from several other cyprinids, but with relatively large intraspecific differences based on SSU rDNA gene sequences. Within our long-term study on myxozoan biodiversity, we performed a broad microscopic and molecular screening of various freshwater fish species (over 450 specimens, 36 species) from different localities. We investigated the cryptic species status of M. pseudodispar. Our analysis revealed four new unique SSU rDNA sequences of M. pseudodispar as well as an infection in new fish host species. Myxobolus pseudodispar sequence analysis showed clear phylogenetic grouping according to fish host criterion forming 13 well-recognised clades. Using 1% SSU rDNA-based genetic distance criterion, at least ten new species of Myxobolus Bütschli, 1882 may be recognised in the group of M. pseudodispar sequences. Our analysis showed the paraphyletic character of M. pseudodispar sequences and the statistical tests rejected hypothetical tree topology with the monophyletic status of the M. pseudodispar group. Myxobolus pseudodispar represents a species complex and it is a typical example of myxozoan hidden diversity phenomenon confirming myxozoans as an evolutionary very successful group of parasites with a great ability to adapt to a new hosts with subsequent speciation events.

Characterization of the Complete Mitochondrial Genome of Epicauta impressicornis (Coleoptera: Meloidae) and Its Phylogenetic Implications for the Infraorder Cucujiformia

The complete mitochondrial genome (mitogenome) of Epicauta impressicornis Pic (Coleoptera: Meloidae) was determined. The circular genome is 15,713-bp long, and encodes 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and a control region (CR). The 13 PCGs start with the typical ATN codon and terminate with the typical stop codon TAA (ND2, ND4L, ND6, ATP6, ATP8, and CYTB), TAG (ND1 and ND3), and T- (COX1, COX2, COX3, ND4, and ND5). The two rRNA genes (rrn12S and rrn16S) are encoded on the minority strand. All tRNAs genes except trnS1 (AGN) are predicted to fold into the typical cloverleaf structure. The longest overlap (10 bp) is observed between ATP8 and ATP6. CR mainly harbors a conserved poly-T stretch (15 bp), a short repeat unit (17 bp), some universal microsatellite-like repeats, and a canonical poly-A tail. Phylogenetic analysis using Bayesian inferences and maximum likelihood based on nucleotide and corresponding amino acid sequences of the 13 PCGs showed that E. impressicornis is closely related to E. chinensis, this relationship is and supported within Cucujiformia belonging to Meloidae (Tenebrionoidea). Our results further confirmed the monophyly of Tenebrionoidea, Lymexyloidea, Curculionoidea, Chrysomeloidea, Cucujoidea, Coccinelloidea, and Cleroidea within Cucujiformia, and revealed the sister relationships of (Cleroidea + Coccinelloidea), (Lymexyloidea + Tenebrionoidea), and ((Chrysomeloidea + Cucujoidea) + Curculionoidea). We believe that the complete mitogenome of E. impressicornis will contribute to further studies on molecular bases for the classification and phylogeny of Meloidae or even Cucujiformia.

The complete mitochondrial genome of a parasite at the animal-fungal boundary

BACKGROUND

Sphaerothecum destruens is an obligate intracellular fish parasite which has been identified as a serious threat to freshwater fishes. Taxonomically, S. destruens belongs to the order Dermocystida within the class Ichthyosporea (formerly referred to as Mesomycetozoea), which sits at the animal-fungal boundary. Mitochondrial DNA (mtDNA) sequences can be valuable genetic markers for species detection and are increasingly used in environmental DNA (eDNA) based species detection. Furthermore, mtDNA sequences can be used in epidemiological studies by informing detection, strain identification and geographical spread.

METHODS

We amplified the entire mitochondrial (mt) genome of S. destruens in two overlapping long fragments using primers designed based on the cox1, cob and nad5 partial sequences. The mt-genome architecture of S. destruens was then compared to close relatives to gain insights into its evolution.

RESULTS

The complete mt-genome of Sphaerothecum destruens is 23,939 bp in length and consists of 47 genes including 21 protein-coding genes, 2 rRNA, 22 tRNA and two unidentified open reading frames. The mitochondrial genome of S. destruens is intronless and compact with a few intergenic regions and includes genes that are often missing from animal and fungal mt-genomes, such as, the four ribosomal proteins (small subunit rps13 and 14; large subunit rpl2 and 16), tatC (twin-arginine translocase component C), and ccmC and ccmF (cytochrome c maturation protein ccmC and heme lyase).

CONCLUSIONS

We present the first mt-genome of S. destruens which also represents the first mt-genome for the order Dermocystida. The availability of the mt-genome can assist the detection of S. destruens and closely related parasites in eukaryotic diversity surveys using eDNA and assist epidemiological studies by improving molecular detection and tracking the parasite's spread. Furthermore, as the only representative of the order Dermocystida, its mt-genome can be used in the study of mitochondrial evolution of the unicellular relatives of animals.

Phylogeny, biogeography and systematics of Pacific vent, methane seep, and whale-fall Parougia (Dorvilleidae : Annelida), with eight new species

Dorvilleidae is a diverse group of annelids found in many marine environments and also commonly associated with chemosynthetic habitats. One dorvilleid genus, Parougia, currently has 11 described species, of which two are found at vents or seeps: Parougia wolfi and Parougia oregonensis. Eight new Parougia species are recognised and described in this study from collections in the Pacific Ocean, all from whale-falls, hydrothermal vents, or methane seeps at ~600-m depth or greater. The specimens were studied using morphology and phylogenetic analyses of DNA sequences from mitochondrial (cytochrome c oxidase subunit I, 16S rRNA, and cytochrome b) and nuclear (18S rRNA and histone 3) genes. Six sympatric Parougia spp. were found at Hydrate Ridge, Oregon, while three of the Parougia species occurred at different types of chemosynthetic habitats. Two new species were found over wide geographical and bathymetric ranges. Another dorvilleid genus, Ophryotrocha, has previously been highlighted as diversifying in the deep-sea environment. Our results document the hitherto unknown diversity of another dorvilleid genus, Parougia, at various chemosynthetic environments. http://zoobank.org/urn:lsid:zoobank.org:pub:EC7EBBEA-2FB5-43D6-BE53-1A468B541A5C

Tick mitochondrial genomes: structural characteristics and phylogenetic implications

Ticks are obligate blood-sucking arachnid ectoparasites from the order Acarina, and many are notorious as vectors of a wide variety of zoonotic pathogens. However, the systematics of ticks in several genera is still controversial. The mitochondrial genome (mt-genome) has been widely used in arthropod phylogeny, molecular evolution and population genetics. With the development of sequencing technologies, an increasing number of tick mt-genomes have been sequenced and annotated. To date, 63 complete tick mt-genomes are available in the NCBI database, and these genomes have become an increasingly important genetic resource and source of molecular markers in phylogenetic studies of ticks in recent years. The present review summarizes all available complete mt-genomes of ticks in the NCBI database and analyses their characteristics, including structure, base composition and gene arrangement. Furthermore, a phylogenetic tree was constructed using mitochondrial protein-coding genes (PCGs) and ribosomal RNA (rRNA) genes from ticks. The results will provide important clues for deciphering new tick mt-genomes and establish a foundation for subsequent taxonomic research.

The complete mitochondrial genomes of five longicorn beetles (Coleoptera: Cerambycidae) and phylogenetic relationships within Cerambycidae

Cerambycidae is one of the most diversified groups within Coleoptera and includes nearly 35,000 known species. The relationships at the subfamily level within Cerambycidae have not been convincingly demonstrated and the gene rearrangement of mitochondrial genomes in Cerambycidae remains unclear due to the low numbers of sequenced mitogenomes. In the present study, we determined five complete mitogenomes of Cerambycidae and investigated the phylogenetic relationship among the subfamilies of Cerambycidae based on mitogenomes. The mitogenomic arrangement of all five species was identical to the ancestral Cerambycidae type without gene rearrangement. Remarkably, however, two large intergenic spacers were detected in the mitogenome of Pterolophia sp. ZJY-2019. The origins of these intergenic spacers could be explained by the slipped-strand mispairing and duplication/random loss models. A conserved motif was found between trnS2 and nad1 gene, which was proposed to be a binding site of a transcription termination peptide. Also, tandem repeat units were identified in the A + T-rich region of all five mitogenomes. The monophyly of Lamiinae and Prioninae was strongly supported by both MrBayes and RAxML analyses based on nucleotide datasets, whereas the Cerambycinae and Lepturinae were recovered as non-monophyletic.

The Complete Mitochondrial Genome and Expression Profile of Mitochondrial Protein-Coding Genes in the Bisexual and Parthenogenetic Haemaphysalis longicornis

The tick Haemaphysalis longicornis is widely distributed in eastern Asia, New Zealand and Australia, and is well-known as a vector of multiple zoonotic pathogens. This species exhibits two reproductive strategies, bisexual and obligate parthenogenetic reproduction. Hence, in the current study, the complete mitochondrial genomes of the bisexual and parthenogenetic populations were assembled and analyzed, and the expression of the mitochondrial protein-coding genes was evaluated and compared between the two reproductive populations. The results indicated that the length of the mitochondrial genomes of the two reproductive populations is 14,694 and 14,693 bp in the bisexual and parthenogenetic populations, respectively. The AT content in the mitochondrial genome of the bisexual and obligate parthenogenetic population reached 77.22 and 77.34%, respectively. The phylogenetic tree was constructed combining 13 protein-coding genes, which showed that the genetic distance between the bisexual and parthenogenetic populations was less than that between the subspecies. The expression of the mitochondrial protein-coding genes was quantitatively analyzed at different feeding status for the bisexual and parthenogenetic populations, and the results showed significant differences in the expression patterns of these genes, suggesting that they might trigger specific energy utilization mechanisms due to their different reproductive strategies and environmental pressures.

Reproductive differences among species, and between individuals and cohorts, in the leech genus Helobdella (Lophotrochozoa; Annelida; Clitellata; Hirudinida; Glossiphoniidae), with implications for reproductive resource allocation in hermaphrodites

Leeches and oligochaetes comprise a monophyletic group of annelids, the Clitellata, whose reproduction is characterized by simultaneous hermaphroditism. While most clitellate species reproduce by cross-fertilization, self-fertilization has been described within the speciose genus Helobdella. Here we document the reproductive life histories and reproductive capacities for three other Helobdella species. Under laboratory conditions, both H. robusta and H. octatestisaca exhibit uniparental reproduction, apparently reflecting self-fertility, and suggesting that this trait is ancestral for the genus. However, the third species, H. austinensis, seems incapable of reproduction by self-fertilization, so we inferred its reproductive life history by analyzing reproduction in breeding cohorts. Comparing the reproductive parameters for H. robusta reproducing in isolation and in cohorts revealed that reproduction in cohorts is dramatically delayed with respect to that of isolated individuals, and that cohorts of leeches coordinate their cocoon deposition in a manner that is not predicted from the reproductive parameters of individuals reproducing in isolation. Finally, our comparisons of reproductive capacity for individuals versus cohorts for H. robusta, and between different sizes of cohorts for H. austinensis, reveal differences in resource allocation between male and female reproductive roles that are consistent with evolutionary theory.

A novel gene arrangement among the Stylommatophora by the complete mitochondrial genome of the terrestrial slug Meghimatium bilineatum (Gastropoda, Arionoidea)

Stylommatophora is a main clade of Gastropoda that encompasses approximately 112 gastropod families and may exceed a total of 30,000 species. Twenty-four complete stylommatophoran mitogenomes have been sequenced to date, yet our understanding of mitochondrial evolution in stylommatophorans is still in its infancy. To further expand the set of available mitogenomes, we sequenced the mitogenome of Meghimatium bilineatum (Arionoidea: Philomycidae), a widespread land slug in East Asia. This is the first report on a mitogenome of the superfamily Arionoidea, and indeed on a terrestrial slug. The mitogenome of Meghimatium bilineatum comprises 13,972 bp and exhibits a novel, highly distinctive gene arrangement among the Stylommatophora. Phylogenetic reconstructions based on the sequences of all protein-coding genes consistently recovered Meghimatium bilineatum as sister-group of the Succineidae. A phylogenetic reconstruction based on gene order, however, suggested a highly divergent tree topology, which is less credible when taking into account prior knowledge of stylommatophoran relationships. Our CREx (Common interval Rearrangement Explorer) analysis suggested that three successive events of tandem duplication random loss (TDRL) best explain the evolutionary process of gene order rearrangement in Meghimatium bilineatum from an ancestral stylommatophoran mitogenome. The present example offers new insights into the mechanisms of mitogenome rearrangements in gastropods at large and into the usefulness of mitogenomic gene order as a phylogenetic marker.

Mitochondrial genome of Paruterina candelabraria (Cestoda: Paruterinidae), with implications for the relationships between the genera Cladotaenia and Paruterina

The taxonomic concept for the family Paruterinidae is controversial, especially concerning the position of the genus Cladotaenia, since the latter genus has been placed sometimes in other families, i.e. in the Taeniidae or in the distinct family Cladotaeniidae; finding a solution based on morphological data is difficult and molecular data on paruterinids and related groups are scarce. In this study, the complete mitochondrial (mt) genome sequence of the type-species of the type-genus of the Paruterinidae, Paruterina candelabraria, was determined and annotated. Gene arrangements are identical with those of Cladotaenia vulturi but differing from those of species of the family Taeniidae by the order change between tRNA-SerUCN and tRNA-LeuCUN. Phylogenetic tree was constructed by Bayesian Inference (BI) analysis using the concatenated amino acid sequences of 12 protein-coding genes. The analysis clearly shows that the Paruterinidae and Taeniidae are sister-groups, and Cladotaenia is a sister taxon of Paruterina. This supports the position of the genus Cladotaenia in the family Paruterinidae and reveals the necessity for sequencing additional taxa of the Paruterinidae for better understanding of phylogenetic relationships within the group.

New patellogastropod mitogenomes help counteracting long-branch attraction in the deep phylogeny of gastropod mollusks

Long-branch attraction (LBA) is a well-known artifact in phylogenetic reconstruction. Sparse taxon sampling and extreme heterogeneity of evolutionary rates among lineages generate propitious situations for LBA, even defying probabilistic methods of phylogenetic inference. A clear example illustrating LBA challenges is the difficulty of reconstructing the deep gastropod phylogeny, particularly using mitochondrial (mt) genomes. Previous studies consistently obtained unorthodox phylogenetic relationships due to the LBA between the mitogenomes of patellogastropods (true limpets, represented only by Lottia digitalis), heterobranchs, and outgroup taxa. Here, we use the reconstruction of the gastropod mitogenomic phylogeny as a case exercise to test the effect of key methodological approaches proposed to counteract LBA, including the selection of slow-evolving representatives, the use of different outgroups, the application of site-heterogeneous evolutionary models, and the removal of fast-evolving sites. In this regard, we sequenced three new patellogastropod mt genomes, which displayed shorter branches than the one of Lottia as well as gene organizations more similar to that of the hypothetical gastropod ancestor. Phylogenetic analyses incorporating the mt genomes of Patella ferruginea, Patella vulgata, and Cellana radiata allowed eliminating the artificial clustering of Patellogastropoda and Heterobranchia that had prevailed in previous studies. Furthermore, the use of site-heterogeneous models with certain combinations of lineages within the outgroup allowed eliminating also the LBA between Heterobranchia and the outgroup, and recovering Apogastropoda (i.e., Caenogastropoda + Heterobranchia). Hence, for the first time, we were able to obtain a mitogenomic phylogeny of gastropods that is congruent with both morphological and nuclear datasets.

Complete mitochondrial genome of Benthodytes marianensis (Holothuroidea: Elasipodida: Psychropotidae): Insight into deep sea adaptation in the sea cucumber

Complete mitochondrial genomes play important roles in studying genome evolution, phylogenetic relationships, and species identification. Sea cucumbers (Holothuroidea) are ecologically important and diverse members, living from the shallow waters to the hadal trench. In this study, we present the mitochondrial genome sequence of the sea cucumber Benthodytes marianensis collected from the Mariana Trench. To our knowledge, this is the first reported mitochondrial genome from the genus Benthodytes. This complete mitochondrial genome is 17567 bp in length and consists of 13 protein-coding genes, two ribosomal RNA genes and 22 transfer RNA genes (duplication of two tRNAs: trnL and trnS). Most of these genes are coded on the positive strand except for one protein-coding gene (nad6) and five tRNA genes which are coded on the negative strand. Two putative control regions (CRs) have been found in the B. marianensis mitogenome. We compared the order of genes from the 10 available holothurian mitogenomes and found a novel gene arrangement in B. marianensis. Phylogenetic analysis revealed that B. marianensis clustered with Peniagone sp. YYH-2013, forming the deep-sea Elasipodida clade. Positive selection analysis showed that eleven residues (24 S, 45 S, 185 S, 201 G, 211 F and 313 N in nad2; 108 S, 114 S, 322 C, 400 T and 427 S in nad4) were positively selected sites with high posterior probabilities. We predict that nad2 and nad4 may be the important candidate genes for the further investigation of the adaptation of B. marianensis to the deep-sea environment.

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

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

Conservation of mitochondrial genome arrangements in brittle stars (Echinodermata, Ophiuroidea)

Brittle stars are conspicuous members of benthic ecosystems, fill many ecological niches and are the most speciose of all classes of echinoderms. With high levels of biodiversity, elucidating the evolutionary history of this group is important. Understanding of higher-level relationships within Ophiuroidea has been aided by multilocus nuclear data and DNA barcoding. However, the degree of consistency between mitochondrial and nuclear data within ophiuroids remains unclear and deserves further assessment. In this study, 17 mitochondrial genomes spanning the taxonomic breadth of Ophiuroidea were utilized to explore evolutionary relationships through maximum likelihood analyses, Bayesian inference and comparative assessment of gene order. Our phylogenetic analyses, based on both nucleotide and amino acid residues, support recent findings based on multilocus nuclear data and morphology, in that the brittle star clades Ophintegrida and Euryophiurida were recovered as monophyletic with the latter comprising Euyalida, Ophiuridae and Ophiopyrgidae. Only three different arrangements of the 13 protein coding and 2 ribosomal RNA genes were observed. As expected, tRNA genes were more likely to have undergone rearrangement but the order of all 37 genes was found to be conserved in all sampled Euryalida and Ophiuridae. Both Euryalida and the clade comprised of Ophiuridae and Ophiopyrgidae, each had their own conserved rearrangement of protein coding genes and ribosomal genes, after divergence from their last common ancestor. Euryalida has a rearrangement of the two ribosomal RNA genes, rrnS and rrnL, in contrast to Ophiuridae and Ophiopyrgidae, which had an inversion of the genes nad1, nad2, and cob relative to Ophintegrida. Further, our data support the gene order found in all sampled Euryalida as the most likely ancestral order for all Ophiuroidea.

The mitochondrial genomes of sarcoptiform mites: are any transfer RNA genes really lost?

BACKGROUND

Mitochondrial (mt) genomes of animals typically contain 37 genes for 13 proteins, two ribosomal RNA (rRNA) genes and 22 transfer RNA (tRNA) genes. In sarcoptiform mites, the entire set of mt tRNA genes is present in Aleuroglyphus ovatus, Caloglyphus berlesei, Dermatophagoides farinae, D. pteronyssinus, Histiostoma blomquisti and Psoroptes cuniculi. Loss of 16 mt tRNA genes, however, was reported in Steganacarus magnus; loss of 2-3 tRNA genes was reported in Tyrophagus longior, T. putrescentiae and Sarcoptes scabiei. Nevertheless, convincing evidence for mt gene loss is lacking in these mites.

RESULTS

We sequenced the mitochondrial genomes of two sarcoptiform mites, Histiostoma feroniarum (13,896 bp) and Rhizoglyphus robini (14,244 bp). Using tRNAScan and ARWEN programs, we identified 16 and 17 tRNA genes in the mt genomes of H. feroniarum and R. robini, respectively. The other six mt tRNA genes in H. feroniarum and five mt tRNA genes in R. robini can only be identified manually by sequence comparison when alternative anticodons are considered. We applied this manual approach to other mites that were reported previously to have lost mt tRNA genes. We were able to identify all of the 16 mt tRNA genes that were reported as lost in St. magnus, two of the three mt tRNA genes that were reported as lost in T. longior and T. putrescentiae, and the two mt tRNA genes that were reported as lost in Sa. scabiei. All of the tRNA genes inferred from these manually identified genes have truncation in the arms and mismatches in the stems.

CONCLUSIONS

Our results reveal very unconventional tRNA structures in sarcoptiform mites and do not support the loss of mt tRNA genes in these mites. The functional implication of the drastic structural changes in these tRNA genes remains to be investigated.

Barcoding and multi-locus phylogeography of the globally distributed calcareous tubeworm genus Hydroides Gunnerus, 1768 (Annelida, Polychaeta, Serpulidae)

Hydroides is a large and diverse group of calcareous tubeworms (Serpulidae, Annelida) recognised by a distinctive but variable two-tiered operculum. Despite considerable research using several species of Hydroides as models in ecological and biofouling studies, phylogenetic and biogeographic relationships within the genus are still poorly understood. Using combined mitochondrial (COI, cytochrome b) and nuclear (18S, 28S and ITS) gene markers for 284 individuals of 45 morphospecies of Hydroides, we investigated the global phylogenetic and biogeographic relationships within the genus. Phylogenetic topologies were well supported and indicated high genetic diversity within Hydroides, revealing potential cryptic species. Present results also include the first COI barcoding data enabling rapid and effective species identification of Hydroides on a global scale. Phylogenetic relationships within Hydroides were more concordant with geographical distributions than morphological similarity of their opercula. Molecular divergence estimates suggested the origin and subsequent diversification in the western Tethys Sea followed by a shift of the historical centre of diversity from the Indo-Mediterranean region to the central Indo-Pacific during the last 50 million years. Further studies on population genetics of species consisting of multiple lineages would provide a better understanding on the status of potential cryptic species. Furthermore, paleogeographic studies based on fossil Hydroides tubes would provide evidence to test this biogeographic hypothesis.

Complete mitochondrial genome of bamboo grasshopper, Ceracris fasciata, and the phylogenetic analyses and divergence time estimation of Caelifera (Orthoptera)

The bamboo grasshopper Ceracris fasciata is regarded as a major pest species because of the damage it causes to bamboo, and its classification within the families and subfamilies of the suborder Caelifera remains unclear. Thus, we attempted to resolve these questions using molecular biology methods and analyses. Our results are as follows: (1) the complete mitochondrial genome of C. fasciata is 15,569 bp in length. The mitochondrial genome contains a standard set of 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes and an A + T-rich region in the same order as those of the other analysed Caeliferan species. The putative start codon for the COX1 gene in C. fasciata is ACC, although it is not defined in other genes. The presence of tandem repeats of different sizes in the A + T-rich region may lead to size differences in other mitochondrial genomes. The mitochondrial genome of C. fasciata harbours the typical 37 genes and an A + T-rich region, and it shows similar characteristics to those of other grasshopper species. Characterization of the mitochondrial genome has enriched our knowledge of the mitochondrial genomes of Orthoptera around the world. Therefore, the phylogenetic relationships in Orthoptera can be re-examined. (2) In phylogenetic analyses, the monophyly of Orthoptera and its two suborders (Caelifera and Ensifera) has been consistently recovered based on most of the datasets selected, regardless of the optimal criteria. Our results do not support the monophyly of the subfamily Oedipodinae of Caelifera. We found that Phlaeoba albonema of the Acridinae is sorted into a clade with Ceracris in all our phylogenetic trees, and field experiments show that Phlaeoba always lives with Ceracris in the same ecotopes. Therefore, we suggest that Phlaeoba should be classified as a member of the Oedipodinae. We found that C. fasciata always clustered with Ceracris kiangsu, and both were sisters to Ceracris versicolor. Therefore, the genetic relationship between C. fasciata and C. kiangsu is closer than that between C. fasciata and C. versicolor. (3) The oldest estimated time of divergence of Ensifera in this context was determined to be 146.16 million years ago (Mya), or around the late Jurassic or early Cretaceous. We estimated that katydids (Grylloidea) likely diverged from other groups in the early Cretaceous. According to our divergence time analyses, we concluded that the ancestral Acrididae probably originated in the early Paleogene, and it is likely that the major diversification events happened at the middle Paleogene, well into the next geologic time. We estimated that crickets (Tettigoniidae) likely diverged from other groups in the early Cretaceous. Acrididae and Romaleinae group, Pyrgacrididae and Ommexechidae group, the youngest two clades we observed, were estimated to have diverged 58.79 Mya, between the middle and early Paleogene. C. versicolor is a sister to the group containing C. kiangsu and C. fasciata. First, C. versicolor diverged from the sister group (C. kiangsu + C. fasciata) around 44.81 Mya, and then the C. kiangsu and C. fasciata group separated at 43.04 Mya.

The mitochondrial genome of the oribatid mite Paraleius leontonychus: new insights into tRNA evolution and phylogenetic relationships in acariform mites

Bilaterian mitochondrial (mt) genomes are circular molecules that typically contain 37 genes. To date, only a single complete mitogenome sequence is available for the species-rich sarcoptiform mite order Oribatida. We sequenced the mitogenome of Paraleius leontonychus, another species of this suborder. It is 14,186 bp long and contains 35 genes, including only 20 tRNAs, lacking tRNA^Gly and tRNA^Tyr. Re-annotation of the mitogenome of Steganacarus magnus increased the number of mt tRNAs for this species to 12. As typical for acariform mites, many tRNAs are highly truncated in both oribatid species. The total number of tRNAs and the number of tRNAs with a complete cloverleaf-like structure in P. leontonychus, however, clearly exceeds the numbers previously reported for Sarcoptiformes. This indicates, contrary to what has been previously assumed, that reduction of tRNAs is not a general characteristic for sarcoptiform mites. Compared to other Sarcoptiformes, the two oribatid species have the least rearranged mt genome with respect to the pattern observed in Limulus polyphemus, a basal arachnid species. Phylogenetic analysis of the newly sequenced mt genome and previously published data on other acariform mites confirms paraphyly of the Oribatida and an origin of the Astigmata within the Oribatida.

Phylogeny, evolution and mitochondrial gene order rearrangement in scale worms (Aphroditiformia, Annelida)

Next-generation sequencing (NGS) has become a powerful tool in phylogenetic and evolutionary studies. Here we applied NGS to recover two ribosomal RNA genes (18S and 28S) from 16 species and 15 mitochondrial genomes from 16 species of scale worms representing six families in the suborder Aphroditiformia (Phyllodocida, Annelida), a complex group of polychaetes characterized by the presence of dorsal elytra or scales. The phylogenetic relationship of the several groups of scale worms remains unresolved due to insufficient taxon sampling and low resolution of individual gene markers. Phylogenetic tree topology based on mitochondrial genomes is comparable with that based on concatenated sequences from two mitochondrial genes (cox1 and 16S) and two ribosomal genes (18S and 28S) genes, but has higher statistical support for several clades. Our analyses show that Aphroditiformia is monophyletic, indicating the presence of elytra is an apomorphic trait. Eulepethidae and Aphroditidae together form the sister group to all other families in this suborder, whereas Acoetidae is sister to Iphionidae. Polynoidae is monophyletic, but within this family the deep-sea subfamilies Branchinotogluminae and Macellicephalinae are paraphyletic. Mitochondrial genomes in most scale-worm families have a conserved gene order, but within Polynoidae there are two novel arrangement patterns in the deep-sea clade. Mitochondrial protein-coding genes in polynoids as a whole have evolved under strong purifying selection, but substitution rates in deep-sea species are much higher than those in shallow-water species, indicating that purifying selection is relaxed in deep-sea polynoids. There are positive selected amino acids for some mitochondrial genes of the deep-sea clade, indicating they may involve in the adaption of deep-sea polynoids. Overall, our study (1) provided more evidence for reconstruction of the phylogeny of Aphroditiformia, (2) provided evidence to refute the assumption that mitochondrial gene order in Errantia is conserved, and (3) indicated that the deep-sea extreme environment may have affected the mitochondrial genome evolution rate and gene order arrangement in Polynoidae.

On the genus Spirobranchus (Annelida, Serpulidae) from the northern Red Sea, and a description of a new species

Species of the genus Spirobranchus, commonly known as Christmas tree worms, are abundant throughout tropical Indo-Pacific and Atlantic Oceans. Information on the species inhabiting the Red Sea in general and the Gulf of Eilat (Gulf of Aqaba) in particular, has so far been very limited. Here we present a multigene phylogenetic analysis, examining both mitochondrial (Cyt-b) and nuclear (ITS2 and 18S) markers, to support the presence of four distinct Spirobranchus species in the Gulf of Eilat: S. corniculatus (including three taxa previously regarded as full species: S. gaymardi, S. cruciger, and S. corniculatus), S. cf. tetraceros, S. gardineri and a new species Spirobranchus aloni, likely endemic to the Red Sea (including two morphotypes with slightly different opercular morphology). The results presented here emphasise that the combination of molecular and in-depth morphological evaluation holds great prospects for a better understanding of species divergence and relationships.

Mitochondrial genomes of blister beetles (Coleoptera, Meloidae) and two large intergenic spacers in Hycleus genera

BACKGROUND

Insect mitochondrial genomes (mitogenomes) exhibit high diversity in some lineages. The gene rearrangement and large intergenic spacer (IGS) have been reported in several Coleopteran species, although very little is known about mitogenomes of Meloidae.

RESULTS

We determined complete or nearly complete mitogenomes of seven meloid species. The circular genomes encode 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs) and two ribosomal RNAs (rRNAs), and contain a control region, with gene arrangement identical to the ancestral type for insects. The evolutionary rates of all PCGs indicate that their evolution is based on purifying selection. The comparison of tRNA secondary structures indicates diverse substitution patterns in Meloidae. Remarkably, all mitogenomes of the three studied Hycleus species contain two large intergenic spacers (IGSs). IGS1 is located between trnW and trnC, including a 9 bp consensus motif. IGS2 is located between trnS2 (UCN) and nad1, containing discontinuous repeats of a pentanucleotide motif and two 18-bp repeat units in both ends. To date, IGS2 is found only in genera Hycleus across all published Coleopteran mitogenomes. The duplication/random loss model and slipped-strand mispairing are proposed as evolutionary mechanisms for the two IGSs (IGS1, IGS2). The phylogenetic analyses using MrBayes, RAxML, and PhyloBayes methods based on nucleotide and amino acid datasets of 13 PCGs from all published mitogenomes of Tenebrionoids, consistently recover the monophylies of Meloidae and Tenebrionidae. Within Meloidae, the genus Lytta clusters with Epicauta rather than with Mylabris. Although data collected thus far could not resolve the phylogenetic relationships within Meloidae, this study will assist in future mapping of the Meloidae phylogeny.

CONCLUSIONS

This study presents mitogenomes of seven meloid beetles. New mitogenomes retain the genomic architecture of the Coleopteran ancestor, but contain two IGSs in the three studied Hycleus species. Comparative analyses of two IGSs suggest that their evolutionary mechanisms are duplication/random loss model and slipped-strand mispairing.

Monoplacophoran mitochondrial genomes: convergent gene arrangements and little phylogenetic signal

BACKGROUND

Although recent studies have greatly advanced understanding of deep molluscan phylogeny, placement of some taxa remains uncertain as different datasets support competing class-relationships. Traditionally, morphologists have placed Monoplacophora, a group of morphologically simple, limpet-like molluscs as sister group to all other conchiferans (shelled molluscs other than Polyplacophora), a grouping that is supported by the latest large-scale phylogenomic study that includes Laevipilina. However, molecular datasets dominated by nuclear ribosomal genes support Monoplacophora + Polyplacophora (Serialia). Here, we evaluate the potential of mitochondrial genome data for resolving placement of Monoplacophora.

RESULTS

Two complete (Laevipilina antarctica and Vema ewingi) and one partial (Laevipilina hyalina) mitochondrial genomes were sequenced, assembled, and compared. All three genomes show a highly similar architecture including an unusually high number of non-coding regions. Comparison of monoplacophoran gene order shows a gene arrangement pattern not previously reported; there is an inversion of one large gene cluster. Our reanalyses of recently published polyplacophoran mitogenomes show, however, that this feature is also present in some chiton species. Maximum Likelihood and Bayesian Inference analyses of 13 mitochondrial protein-coding genes failed to robustly place Monoplacophora and hypothesis testing could not reject any of the evaluated placements of Monoplacophora.

CONCLUSIONS

Under both serialian or aculiferan-conchiferan scenarios, the observed gene cluster inversion appears to be a convergent evolution of gene arrangements in molluscs. Our phylogenetic results are inconclusive and sensitive to taxon sampling. Aculifera (Polyplacophora + Aplacophora) and Conchifera were never recovered. However, some analyses recovered Serialia (Monoplacophora + Polyplacophora), Diasoma (Bivalvia + Scaphopoda) or Pleistomollusca (Bivalvia + Gastropoda). Although we could not shed light on deep evolutionary traits of Mollusca we found unique patterns of gene arrangements that are common to monoplacophoran and chitonine polyplacophoran species but not to acanthochitonine Polyplacophora. Complete mitochondrial genome of Laevipilina antarctica.

The complete mitochondrial genome of the tapeworm Cladotaenia vulturi (Cestoda: Paruterinidae): gene arrangement and phylogenetic relationships with other cestodes

BACKGROUND

Tapeworms Cladotaenia spp. are among the most important wildlife pathogens in birds of prey. The genus Cladotaenia is placed in the family Paruterinidae based on morphological characteristics and hosts. However, limited molecular information is available for studying the phylogenetic position of this genus in relation to other cestodes.

METHODS

In this study, the complete mitochondrial (mt) genome of Cladotaenia vulturi was amplified using "Long-PCR" and then sequenced by primer walking. Sequence annotation and gene identification were performed by comparison with published flatworm mt genomes. The phylogenetic relationships of C. vulturi with other cestode species were established using the concatenated amino acid sequences of 12 protein-coding genes with Bayesian Inference and Maximum Likelihood methods.

RESULTS

The complete mitochondrial genome of the Cladotaenia vulturi is 13,411 kb in size and contains 36 genes. The gene arrangement of C. vulturi is identical to those in Anoplocephala spp. (Anoplocephalidae), Hymenolepis spp. (Hymenolepididae) and Dipylidium caninum (Dipylidiidae), but different from that in taeniids owing to the order shift between the tRNA (L1) and tRNA (S2) genes. Phylogenetic analyses based on the amino acid sequences of the concatenated 12 protein-coding genes showed that the species in the Taeniidae form a group and C. vulturi is a sister taxon to the species of the family Taeniidae.

CONCLUSIONS

To our knowledge, the present study provides the first molecular data to support the early proposal from morphological evidence that the Taeniidae is a sister group to the family Paruterinidae. This novel mt genome sequence will be useful for further investigations into the population genetics, phylogenetics and systematics of the family Paruterinidae and inferring phylogenetic relationships among several lineages within the order Cyclophyllidea.

In silico analyses of mitochondrial ORFans in freshwater mussels (Bivalvia: Unionoida) provide a framework for future studies of their origin and function

BACKGROUND

Many species of bivalves exhibit a unique system of mtDNA transmission named Doubly Uniparental Inheritance (DUI). Under this system, species have two distinct, sex-linked mitochondrial genomes: the M-type mtDNA, which is transmitted by males to male offspring and found in spermatozoa, and the F-type mtDNA, which is transmitted by females to all offspring, and found in all tissues of females and in somatic tissues of males. Bivalves with DUI also have sex-specific mitochondrial ORFan genes, (M-orf in the M mtDNA, F-orf in the F mtDNA), which are open reading frames having no detectable homology and no known function. DUI ORFan proteins have previously been characterized in silico in a taxonomically broad array of bivalves including four mytiloid, one veneroid and one unionoid species. However, the large evolutionary distance among these taxa prevented a meaningful comparison of ORFan properties among these divergent lineages. The present in silico study focuses on a suite of more closely-related Unionoid freshwater mussel species to provide more reliably interpretable information on patterns of conservation and properties of DUI ORFans. Unionoid species typically have separate sexes, but hermaphroditism also occurs, and hermaphroditic species lack the M-type mtDNA and possess a highly mutated version of the F-orf in their maternally transmitted mtDNA (named H-orf in these taxa). In this study, H-orfs and their respective proteins are analysed for the first time.

RESULTS

Despite a rapid rate of evolution, strong structural and functional similarities were found for M-ORF proteins compared among species, and among the F-ORF and H-ORF proteins across the studied species. In silico analyses suggest that M-ORFs have a role in transport and cellular processes such as signalling, cell cycle and division, and cytoskeleton organisation, and that F-ORFs may be involved in cellular traffic and transport, and in immune response. H-ORFs appear to be structural glycoproteins, which may be involved in signalling, transport and transcription. Our results also support either a viral or a mitochondrial origin for the ORFans.

CONCLUSIONS

Our findings reveal striking structural and functional similarities among proteins encoded by mitochondrial ORFans in freshwater mussels, and strongly support a role for these genes in the DUI mechanism. Our analyses also support the possibility of DUI systems with elements of different sources/origins and different mechanisms of action in the distantly-related DUI taxa. Parallel situations to the novel mitochondrially-encoded functions of freshwater mussel ORFans present in some other eukaryotes are also discussed.