Phylogenetic position of the Pentastomida and (pan)crustacean relationships

Gaps in parasitological research in the molecular era

We live in the age of molecular biology and '-omics', and molecular methods have opened up unimagined possibilities for biological research, including parasitology. However, too one-sided a focus on new approaches can lead to major gaps as less 'cool' topics are neglected. Selected areas of research are briefly discussed to highlight the gaps caused by the current excessive focus on molecular and '-omics' methods. It is crucial to combine both 'classical' and modern methods without neglecting the complexity of the interactions of parasites with their hosts and the environment (One Health concept), which is even more urgent in today's rapidly changing world. Parasitologists should be more involved in field studies and multidisciplinary assessment of parasites.

Complete Mitogenome sequencing of the fish louse Argulus japonicus (Crustacea: Branchiura): Comparative analyses and phylogenetic implications

The fish louse Argulus japonicus, a branchiuran crustacean of the Argulidae family, is attracting increasing attention because of its parasitic tendencies and significant health threats to global fish farming. The mitogenomes can yield a foundation for studying epidemiology, genetic diversity, and molecular ecology and therefore may be used to assist in the surveillance and control of A. japonicus. In this study, we sequenced and assembled the complete mitogenome of A. japonicus to shed light on its genetic and evolutionary blueprint. Our investigation indicated that the 15,045-bp circular genome of A. japonicus encodes 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), and 2 ribosomal RNAs (rRNAs) with significant AT and GC skews. Comparative genomics provided an evolutionary scenario for the genetic diversity of 13 PCGs: all were under purifying selection, with cox1 and nad6 having the lowest and highest evolutionary rates, respectively. Genome-wide phylogenetic trees established a close relationship between species of the families Argulidae (Arguloida) and Armilliferidae (Porocephalida) within Crustacea, and further, A. japonicus and Argulus americanus were determined to be more closely related to each other than to others within the family Argulidae. Single PCG-based phylogenies supported nad1 and nad6 as the best genetic markers for evolutionary and phylogenetic studies for branchiuran crustaceans due to their similar phylogenetic topologies with those of genome-based phylogenetic analyses. To sum up, these comprehensive mitogenomic data of A. japonicus and related species refine valuable marker resources and should contribute to molecular diagnostic methods, epidemiological investigations, and ecological studies of the fish ectoparasites in Crustacea.

The mitochondrial genome of the pentastome parasite Raillietiella orientalis Hett, 1915 (Raillietiellida; Raillietiellidae) with notes on its phylogenetic position

In this study we sequenced and annotated the complete mitochondrial genome of the invasive reptile parasite Raillietiella orientalis using Illumina DNA sequencing. The length of the mitogenome was 15,320 bp and had a GC content of 33.1%. The mitogenome contained 13 protein-coding genes, two ribosomal RNA genes, and 22 tRNA genes, the order of which was diagnostic of Pancrustacean mitogenomes. A phylogenetic tree constructed from the 13 protein-coding genes of R. orientalis and 26 other Pancrustacean mitogenomes supported the placement of R. orientalis as part of the monophyletic subclass Pentastomida within the Maxillopoda and sister to the subclass Branchiura.

Conflicts in Mitochondrial Phylogenomics of Branchiopoda, with the First Complete Mitogenome of Laevicaudata (Crustacea: Branchiopoda)

Conflicting phylogenetic signals are pervasive across genomes. The potential impact of such systematic biases may be reduced by phylogenetic approaches accommodating for heterogeneity or by the exclusive use of homoplastic sites in the datasets. Here, we present the complete mitogenome of Lynceus grossipedia as the first representative of the suborder Laevicaudata. We employed a phylogenomic approach on the mitogenomic datasets representing all major branchiopod groups to identify the presence of conflicts and concordance across the phylogeny. We found pervasive phylogenetic conflicts at the base of Diplostraca. The homogeneity of the substitution pattern tests and posterior predictive tests revealed a high degree of compositional heterogeneity among branchiopod mitogenomes at both the nucleotide and amino acid levels, which biased the phylogenetic inference. Our results suggest that Laevicaudata as the basal clade of Phyllopoda was most likely an artifact caused by compositional heterogeneity and conflicting phylogenetic signal. We demonstrated that the exclusive use of homoplastic site methods combining the application of site-heterogeneous models produced correct phylogenetic estimates of the higher-level relationships among branchiopods.

The Mitochondrial Genome of the Globally Invasive Barnacle Megabalanus coccopoma Darwin 1854 (Crustacea: Balanomorpha): Rearrangement and Phylogenetic Consideration within Balanomorpha

Megabalanus coccopoma (Darwin, 1854) is a globally invasive species in Balanomorpha (Crustacea). This species is a model organism for studying marine pollution and ecology. However, its mitogenome remains unknown. The mitogenome sequencing of M. coccopoma is completed in the present study. It has a 15,098 bp in length, including 13 protein-coding genes (PCGs), 2 ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs), along with a putative regulatory area. A substantial A+T bias was observed in the genome composition (68.2%), along with a negative AT (0.82) and GC (−0.136) skew. Compared to the gene sequence of the ground model of pan-crustacea, 13 gene clusters (or genes), such as 10 tRNAs and 3 PCGs, were observed in a different order. This was in line with the previously observed large-scale gene rearrangements of Balanomorpha. Among the 37 genes, the gene cluster (M-nad2-W-cox1-L2-cox2-D-atp8-atp6-cox3-G- nad3-R-N-A-E-S1) Balanomorpha was conserved. Furthermore, phylogeny analysis indicated that the existing Balanomorpha species family was divided into nine rearrangement patterns, supporting the polyphyly of Balanoidea.

The complete mitochondrial genome of the Arctic fairy shrimp Branchinectapaludosa (Müller, 1788) (Anostraca, Branchinectidae) from Sirius Passet, North Greenland

Here we report the complete mitochondrial genome of the Arctic fairy shrimp, Branchinectapaludosa (Müller, 1788) (Anostraca, Branchinectidae), which was collected in the High Arctic of North Greenland. A complete 16,059 bp mitochondrion of B.paludosa was sequenced and assembled with the Illumina next generation sequencing platform. The B.paludosa mitogenome contains 13 PCGs, 22 tRNAs and 2 rRNA genes that are commonly observed in most metazoans and shows the conserved gene arrangement pattern of Anostraca. Our results of the phylogenomic analysis are consistent with the previous phylogenetic relationship, based on nuclear 18S ribosomal DNA. The B.paludosa mitogenome will be useful for understanding the geographical distribution and phylogenetic relationship of anostracans.

The complete mitochondrial genome of the Antarctic fairy shrimp Branchinectagaini Daday, 1910 (Branchiopoda, Anostraca, Branchinectidae)

The complete mitochondrial genome of Antarctic fairy shrimp Branchinectagaini Daday, 1910 was sequenced, assembled and annotated using next-generation sequencing technology. The mitogenome of B.gaini is circular at 15,536 bp in length, consisting of 13 protein-coding genes, 23 tRNAs, two rRNAs and two major non-coding regions. In particular, there are two tRNA^Gly genes and one non-coding region between these two tRNA^Gly genes. A phylogenetic tree was constructed using concatenated amino acid sequences of 13 protein-coding genes. It reveals that B.gaini is clustered with the Anostraca group within the Branchiopoda clade. This study helps us understand the evolution of Anostraca.

A Rearrangement of the Mitochondrial Genes of Centipedes (Arthropoda, Myriapoda) with a Phylogenetic Analysis

Due to the limitations of taxon sampling and differences in results from the available data, the phylogenetic relationships of the Myriapoda remain contentious. Therefore, we try to reconstruct and analyze the phylogenetic relationships within the Myriapoda by examining mitochondrial genomes (the mitogenome). In this study, typical circular mitogenomes of Mecistocephalus marmoratus and Scolopendra subspinipes were sequenced by Sanger sequencing; they were 15,279 bp and 14,637 bp in length, respectively, and a control region and 37 typical mitochondrial genes were annotated in the sequences. The results showed that all 13 PCGs started with ATN codons and ended with TAR codons or a single T; what is interesting is that the gene orders of M. marmoratus have been extensively rearranged compared with most Myriapoda. Thus, we propose a simple duplication/loss model to explain the extensively rearranged genes of M. marmoratus, hoping to provide insights into mitogenome rearrangement events in Myriapoda. In addition, our mitogenomic phylogenetic analyses showed that the main myriapod groups are monophyletic and supported the combination of the Pauropoda and Diplopoda to form the Dignatha. Within the Chilopoda, we suggest that Scutigeromorpha is a sister group to the Lithobiomorpha, Geophilomorpha, and Scolopendromorpha. We also identified a close relationship between the Lithobiomorpha and Geophilomorpha. The results also indicate that the mitogenome can be used as an effective mechanism to understand the phylogenetic relationships within Myriapoda.

Comparative analysis of mitochondrial genomes among the family Peltoperlidae (Plecoptera: Systellognatha) and phylogenetic implications

Nowadays, the position of Peltoperlidae in Systellognatha has been resolved based on morphological analyses. However, there are different opinions based on molecular data. To date, only three peltoperlid mitogenomes are available, and more sampling is needed to obtain precise phylogenetic relationships. In this study, we obtained the complete mitogenomes of Cryptoperla kawasawai (15,832 bp) and Peltoperlopsis sagittata (15,756 bp). Our results show that gene content, gene order, DmTTF binding site, nucleotide composition, codon usage, ribonucleic acid (RNA) structure, and structural elements in the control region are highly conserved in peltoperlids. Heatmap analysis of codon usage shows that the AT-rich codons UUA, AUU, UUU, and AUA were commonly used codons in the Peltoperlidae. Evolutionary rate analyses of protein-coding genes reveal that different genes have been subject to different rates of molecular evolution correlated with the GC content. All tRNA genes in peltoperlid mitogenomes have a canonical cloverleaf secondary structure except for trnS1, whose dihydrouridine arm simply forms a loop. The control region of the family has several distinct structural characteristics and has the potential to serve as effective phylogenetic markers. Phylogenetic analyses support the monophyly of Perloidea, but the monophyly of Pteronarcyoidea is still not supported. The Peltoperlidae is placed as the earliest branch within the Systellognatha, and the estimated phylogenetic relationship is: Peltoperlidae + {(Styloperlidae + Pteronarcyidae) + [Perlidae + (Chloroperlidae + Perlodidae)]}. Our results provide new insight into the phylogeny of this group.

Insights into the Evolution of Aphid Mitogenome Features from New Data and Comparative Analysis

The complete mitochondrial genomes and their rearrangement patterns can provide useful information for inferring evolutionary history of organisms. Aphids are one of the insect groups with some unique mitogenome features. In this study, to examine whether some features in aphid mitogenomes are independent species-specific evolutionary events or clade-specific events at certain taxonomic levels, we sequenced three new aphid mitogenomes (Hormaphidinae: Ceratovacuna keduensis, Pseudoregma panicola; Lachninae: Nippolachnus piri) and compared them with all known aphid mitogenomes. The three mitogenomes are 16,059–17,033 bp in length, with a set of 37 typical mitochondrial genes, a non-coding control region and a tandem repeat region. The gene orders of them are all highly rearranged. Within the subfamily Hormaphidinae, the presence of repeat region and mitogenome rearrangement in Cerataphidini species but not in the other two tribes indicate that these may be Cerataphidini-specific features. The same gene rearrangement pattern in the two Lachninae species, N. piri (Tuberolachnini) and Stomaphis sinisalicis (Stomaphidini), supports that this feature should be at least derived from the common ancestor of two tribes. Overall, our data and analyses provide new insights into the evolutionary patterns of gene rearrangement and repeat region in aphid mitogenomes, and further corroborate the potential role of gene rearrangement in elucidating the evolutionary history of different insect lineages.

Linguatula serrata in an imported dog in Germany: Single-case or emerging disease?

Linguatula serrata is a worm-like parasite with zoonotic potential that inhabits the nasal cavities of canids. Although most cases of linguatulosis are associated with unspecific and rather mild respiratory symptoms, cases of unusual infestations and severe courses in both animals and humans have been reported. In central and northern Europe, the pathogen used to appear only sporadically, however, within the last few years the number of detections has increased noticeably. In July 2020 an approximately nine-month-old dog, imported from Romania, was presented in a veterinary practice in Gotha, central Germany, due to persistent worsening cough. Despite antibiotic treatment the tussis became more severe until the dog expectorated multiple worm-like structures. Three of these specimens were sent to the Institute of Parasitology (Faculty of Veterinary Medicine, University of Leipzig, Leipzig) for morphological and genetic species identification. The latter was based on a 1000-bp fragment of the mitochondrial cytochrome c oxidase subunit I gene (cox1) and the complete nuclear 18S rRNA gene. The dog presented in this study suffered from a severe respiratory impairment caused by worm-like parasites inhabiting its upper respiratory tract. The detected parasites were morphologically identified as female specimens of the so-called tongue-worm L. serrata, which was confirmed by pairwise alignment and phylogenetic analysis of the produced sequences. We report an unusually severe case of L. serrata infection in an imported dog and discuss the spread of this potentially dangerous parasite in central and northern Europe.

Raillietiella hemidactyli (Pentastomida: Raillietiellidae) in Tarentola mauritanica geckoes: A new zoonotic parasite for Europe

Pentastomida is a subclass of parasitic arthropods, related to crustaceans, which develop in the respiratory tract of vertebrates (i.e., fishes, amphibians, reptiles, birds and mammals). Within this group of parasites, Raillietiella spp. adults develop in the lungs of lizards, snakes and toads, whereas larval stages in insects (e.g., cockroaches), which are intermediate hosts. Lizards were captured under the frame of a study on reptile zoonotic parasites. Feces of the collected animals were examined and pentastomids were diagnosed in Tarentola mauritanica geckoes (1.2%; 3/259) from Linosa island. Adult forms of Railietiella hemidactyli pentastomids were morphologically characterized and molecularly identified through 18S rDNA amplification and sequencing. Positive animals had adult forms of R. hemidactyli pentastomids in the lungs as well as embryonated eggs in feces. Raillietiella was herein identified for the first time in synanthropic geckoes in a confined population of one of the southernmost islands of Italy, representing the first report of this zoonotic pentastomid in synanthropic and invasive reptiles in Europe. Further studies should focus on the prevalence of pentastomids on synanthropic reptiles in other Italian regions to assess the zoonotic risk of infection and to warn veterinarians and physicians about the risk they may represent for several species of hosts, including dogs, cats and humans.

Chromosome-level genome assembly, annotation, and phylogenomics of the gooseneck barnacle Pollicipes pollicipes

BACKGROUND

The barnacles are a group of >2,000 species that have fascinated biologists, including Darwin, for centuries. Their lifestyles are extremely diverse, from free-swimming larvae to sessile adults, and even root-like endoparasites. Barnacles also cause hundreds of millions of dollars of losses annually due to biofouling. However, genomic resources for crustaceans, and barnacles in particular, are lacking.

RESULTS

Using 62× Pacific Biosciences coverage, 189× Illumina whole-genome sequencing coverage, 203× HiC coverage, and 69× CHi-C coverage, we produced a chromosome-level genome assembly of the gooseneck barnacle Pollicipes pollicipes. The P. pollicipes genome is 770 Mb long and its assembly is one of the most contiguous and complete crustacean genomes available, with a scaffold N50 of 47 Mb and 90.5% of the BUSCO Arthropoda gene set. Using the genome annotation produced here along with transcriptomes of 13 other barnacle species, we completed phylogenomic analyses on a nearly 2 million amino acid alignment. Contrary to previous studies, our phylogenies suggest that the Pollicipedomorpha is monophyletic and sister to the Balanomorpha, which alters our understanding of barnacle larval evolution and suggests homoplasy in a number of naupliar characters. We also compared transcriptomes of P. pollicipes nauplius larvae and adults and found that nearly one-half of the genes in the genome are differentially expressed, highlighting the vastly different transcriptomes of larvae and adult gooseneck barnacles. Annotation of the genes with KEGG and GO terms reveals that these stages exhibit many differences including cuticle binding, chitin binding, microtubule motor activity, and membrane adhesion.

CONCLUSION

This study provides high-quality genomic resources for a key group of crustaceans. This is especially valuable given the roles P. pollicipes plays in European fisheries, as a sentinel species for coastal ecosystems, and as a model for studying barnacle adhesion as well as its key position in the barnacle tree of life. A combination of genomic, phylogenetic, and transcriptomic analyses here provides valuable insights into the evolution and development of barnacles.

Comparative Mitogenomes of Two Coreamachilis Species (Microcoryphia: Machilidae) along with Phylogenetic Analyses of Microcoryphia

The order Microcoryphia, commonly known as bristletails, is considered as the most primitive one among living insects. Within this order, two species, Coreamachilis coreanus and C. songi (Machilidae: Machilinae), display the following contrasting reproductive strategies: parthenogenesis occurs in C. coreanus, whereas sexual reproduction is found in C. songi. In the present study, the complete mitogenomes of C. coreanus and C. songi were sequenced to compare their mitogenome structure, analyze relationships within the Microcoryphia, and assess adaptive evolution. The length of the mitogenomes of C. coreanus and C. songi were 15,578 bp and 15,570 bp, respectively, and the gene orders were those of typical insects. A long hairpin structure was found between the ND1 and 16S rRNA genes of both species that seem to be characteristic of Machilinae and Petrobiinae species. Phylogenetic assessment of Coreamachilis was conducted using BI and ML analyses with concatenated nucleotide sequences of the 13 protein-coding genes. The results showed that the monophyly of Machilidae, Machilinae, and Petrobiinae was not supported. The genus Coreamachilis (C. coreanus and C. songi) was a sister clade to Allopsontus helanensis, and then the clade of ((C. coreanus + C. songi) + A. helanensis) was a sister clade to A. baii, which suggests that the monophyly of Allopsontus was not supported. Positive selection analysis of the 13 protein-coding genes failed to reveal any positive selection in C. coreanus or C. songi. The long hairpin structures found in Machilinae and Petrobiinae were highly consistent with the phylogenetic results and could potentially be used as an additional molecular characteristic to further discuss relationships within the Microcoryphia.

Mitochondrial Genomes from Two Specialized Subfamilies of Reduviidae (Insecta: Hemiptera) Reveal Novel Gene Rearrangements of True Bugs

Reduviidae, a hyper-diverse family, comprise 25 subfamilies with nearly 7000 species and include many natural enemies of crop pests and vectors of human disease. To date, 75 mitochondrial genomes (mitogenomes) of assassin bugs from only 11 subfamilies have been reported. The limited sampling of mitogenome at higher categories hinders a deep understanding of mitogenome evolution and reduviid phylogeny. In this study, the first mitogenomes of Holoptilinae (Ptilocnemus lemur) and Emesinae (Ischnobaenella hainana) were sequenced. Two novel gene orders were detected in the newly sequenced mitogenomes. Combined 421 heteropteran mitogenomes, we identified 21 different gene orders and six gene rearrangement units located in three gene blocks. Comparative analyses of the diversity of gene order for each unit reveal that the tRNA gene cluster trnI-trnQ-trnM is the hotspot of heteropteran gene rearrangement. Furthermore, combined analyses of the gene rearrangement richness of each unit and the whole mitogenome among heteropteran lineages confirm Reduviidae as a 'hot-spot group' of gene rearrangement in Heteroptera. The phylogenetic analyses corroborate the current view of phylogenetic relationships between basal groups of Reduviidae with high support values. Our study provides deeper insights into the evolution of mitochondrial gene arrangement in Heteroptera and the early divergence of reduviids.

Morphology, ultrastructure, genomics, and phylogeny of Euplotes vanleeuwenhoeki sp. nov. and its ultra-reduced endosymbiont "Candidatus Pinguicoccus supinus" sp. nov

Taxonomy is the science of defining and naming groups of biological organisms based on shared characteristics and, more recently, on evolutionary relationships. With the birth of novel genomics/bioinformatics techniques and the increasing interest in microbiome studies, a further advance of taxonomic discipline appears not only possible but highly desirable. The present work proposes a new approach to modern taxonomy, consisting in the inclusion of novel descriptors in the organism characterization: (1) the presence of associated microorganisms (e.g.: symbionts, microbiome), (2) the mitochondrial genome of the host, (3) the symbiont genome. This approach aims to provide a deeper comprehension of the evolutionary/ecological dimensions of organisms since their very first description. Particularly interesting, are those complexes formed by the host plus associated microorganisms, that in the present study we refer to as "holobionts". We illustrate this approach through the description of the ciliate Euplotes vanleeuwenhoeki sp. nov. and its bacterial endosymbiont "Candidatus Pinguicoccus supinus" gen. nov., sp. nov. The endosymbiont possesses an extremely reduced genome (~ 163 kbp); intriguingly, this suggests a high integration between host and symbiont.

The complete mitochondrial genome of the pentastomid Linguatula arctica (Pentastomida) from reindeer (Rangifer tarandus) in Northern Norway

Here, we present the first complete mitochondrial genome of the pentastomid Linguatula arctica collected from the nasal passages of a reindeer (Rangifer tarandus) in Norway. The full length mitochondrial genome of L. arctica, which measures 14,789 bp in length, contains 13 protein-coding genes, 2 ribosomal RNA genes and 22 transfer RNA genes. A clear A + T bias is observed in the mitogenome of L. arctica with an overall base composition of 32.6% A, 27.5% T, 32.8% C, and 7,1% G., and a GC content of 39.9%. The gene arrangement is identical to that of previously described pentastomid mitogenomes.

Detection of a Larva of Armillifer armillatus in a Potto (Perodicticus potto) from the Republic of the Congo

We determined the complete sequence of the mitochondrial DNA (mtDNA) of a parasite discovered between the subcutaneous tissue and the peritoneum of an African nocturnal non-human primate (NHP). The parasite and host sequences were obtained by a combination of Sanger sequencing and nanopore MinION techniques. Analyses of mtDNA gene arrangements and sequences unambiguously showed that the parasite investigated was the pentastomid Armillifer armillatus, also commonly named the tongue worm. The full-length mitochondrial genome of A. armillatus, measuring 16,706 bp in length, contains 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes, an arrangement identical to that of previously described pentastomid mitochondrial genomes. We describe here the second full mitochondrial genome of A. armillatus to date. To identify the NHP host, maximum likelihood phylogenetic analyses of a 441-bp fragment on the 12S rDNA gene and of a 1,140-bp fragment of the mitochondrial cytochrome b strongly support clustering with the African lorisid Perodicticus potto, a species that has rarely been reported as an intermediate host of this parasite.

Morphological and molecular description of Armillifer moniliformis larvae isolated from Sri Lankan brown palm civet (Paradoxurus montanus)

We report Armillifer moniliformis species infecting the endemic Sri Lankan brown palm civet (Paradoxurus montanus) from the Knuckles Range Forest Conservation Area, Sri Lanka. Larval stages of A. moniliformis were found during the postmortem of three civet cats found dead. Morphological studies were done by a light microscope and a scanning electron microscope (SEM). Histopathological examination was conducted using tissue samples obtained from the liver. For the molecular analysis, DNA was extracted from the isolated third-stage larvae. The NADH dehydrogenase subunit 5 (ND5) and the second internal transcribed spacer region (ITS-2), a portion of the large subunit nuclear ribosomal DNA (28S), a portion of 18S ribosomal rRNA gene (18S), and cytochrome c oxidase subunit 1 gene (COX1) were amplified using polymerase chain reaction (PCR). Excysted third-stage larvae were observed in the lungs, omentum, the pleural cavity, the abdominal cavity, and the surface of the spleen and the pericardium. Around 88 third-stage larvae were isolated from three civet cats. First-stage larvae in the liver were surrounded by outer fibrous layer over the inner germinal layer and filled with clear fluid. Slight hemorrhage, leukocyte infiltration, and mild hepatocellular degeneration in the liver were observed. The SEM examination indicated the unique oral apparatus comprises the oval-shaped mouth opening in between two pairs of curved, retractile hamuli. The sequences obtained for ND5, ITS-2, 28S, 18S, and COX1 were 301, 382, 325, 414, and 644 bp in length respectively. Morphology, sequence similarity search, sequence alignment, and phylogenetic analysis identified this parasite as A. moniliformis.

Unexpected Infection with Armillifer Parasites

Visceral pentastomiasis is usually found incidentally during surgery. We describe a case of visceral pentastomiasis discovered during inguinoscrotal hernia surgery for a man from Benin, Africa. Because surgical removal of nymphs is needed for symptomatic patients only, this patient's asymptomatic pentastomiasis was not treated and he recovered from surgery uneventfully.

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

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

Co-infections in Visceral Pentastomiasis, Democratic Republic of the Congo

Results of PCR and histology indicate this infection is endemic to this country.

Snakeborne Armillifer pentastomiasis is an emerging human parasitic infection in rural tropical areas where snake meat is eaten. After a series of severe ocular A. grandis larval infections and anecdotal abdominal infection in Sankuru District, Democratic Republic of the Congo, during 2014–2015, we systematically investigated possible pentastomid etiology in patients who underwent surgery in the region. Histologic and molecular analyses by established pentastomid 18S rDNA- and newly developed Armillifer-specific cytochrome oxidase PCRs revealed larval pentastomid lesions in 3.7% of patients. Some persons had A. armillatus and A. grandis co-infections. Another pentastomid larva, Raillietiella sp., was molecularly detected in 1 patient who had concomitant A. grandis and A. armillatus infection. The PCRs used were suitable for detecting pentastomid species even in highly necrotic tissues. Phylogenetic analyses of Armillifer cytochrome oxidase genes detected multiple local strains.

More limbs on the tree: mitogenome characterisation and systematic position of ‘living fossil’ species Neoglyphea inopinata and Laurentaeglyphea neocaledonica (Decapoda : Glypheidea : Glypheidae)

Glypheids first appeared in the Lower Triassic period and were believed to be extinct until specimens of Neoglyphea inopinata Forest & Saint Laurent and Laurentaeglyphea neocaledonica Richer de Forges were described in 1975 and 2006, respectively. The finding of extant species has meant that molecular data can now be used to complement morphological and fossil-based studies to investigate the relationships of Glypheidea within the Decapoda. However, despite several molecular studies, the placement of this infraorder within the decapod phylogenetic tree is not resolved. One limitation is that molecular resources available for glypheids have been limited to a few nuclear and mitochondrial gene fragments. Many of the more recent large-scale studies of decapod phylogeny have used information from complete mitogenomes, but have excluded the infraorder Glypheidea due to the unavailability of complete mitogenome sequences. Using next-generation sequencing, we successfully sequenced and assembled complete mitogenome sequences from museum specimens of N. inopinata and L. neocaledonica, the only two extant species of glypheids. With these sequences, we constructed the first decapod phylogenetic tree based on whole mitogenome sequences that includes Glypheidea as one of 10 decapod infraorders positioned within the suborder Pleocyemata. From this, the Glypheidea appears to be a relatively derived lineage related to the Polychelida and Astacidea. Also in our study, we conducted a survey on currently available decapod mitogenome resources available on National Center for Biotechnology Information (NCBI) and identified infraorders that would benefit from more strategic and expanded taxonomic sampling.

Armillifer-Infected Snakes Sold at Congolese Bushmeat Markets Represent an Emerging Zoonotic Threat

African pythons (Pythonidae) and large vipers (Bitis spp.) act as definitive hosts for Armillifer armillatus and Armillifer grandis parasites (Crustacea: Pentastomida) in the Congo Basin. Since the proportion of snakes in bushmeat gradually increases, human pentastomiasis is an emerging zoonotic disease. To substantiate the significance of this threat, we surveyed snakes offered for human consumption at bushmeat markets in the Kole district, Democratic Republic of the Congo, for the presence of adult pentastomids. In Bitis vipers (n = 40), Armillifer spp. infestations exhibited an 87.5% prevalence and 6.0 median intensity. Parasite abundance covaried positively with viper length, but not with body mass. In pythons (n = 13), Armillifer spp. exhibited a 92.3% prevalence and 3.5 median intensity. The positive correlations between parasite abundance and python length or mass were statistically nonsignificant. Ninety-one percent of A. grandis were discovered in vipers and 97% of infected vipers hosted A. grandis, whereas 81% of A. armillatus specimens were found in pythons and 63% of infected pythons hosted A. armillatus. Thus, challenging the widespread notion of strict host specificity, we found 'reversed' infections and even a case of coinfection. In this study, we also gathered information about the snake consumption habits of different tribal cultures in the area. Infective parasite ova likely transmit to humans directly by consumption of uncooked meat, or indirectly through contaminated hands, kitchen tools or washing water.

The complete mitochondrial genome of the pentastomid Armillifer grandis (Pentastomida) from the Democratic Republic of Congo

We present the first complete mitochondrial genome of the pentastomid Armillifer grandis (Arthropoda: Pentastomida) collected from the lungs of a rhinoceros viper (Bitis nasicornis) in the Democratic Republic of Congo. The full length mitochondrial genome of Armillifer grandis, which measures 16,073 bp in length, contains 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. A clear A + T bias is observed in the mitogenome of Armillifer grandis with an overall base composition of 34.6% A, 29.4% T, 29% C, and 6.9% G, and a GC content of 35.9%. The gene arrangement is identical to that of previously described pentastomid mitogenomes.

Comparison of Three Different DNA Extraction Methods for Linguatula serrata as a Food Born Pathogen

BACKGROUND

One of the most important items in molecular characterization of food-borne pathogens is high quality genomic DNA. In this study, we investigated three protocols and compared their simplicity, duration and costs for extracting genomic DNA from Linguatula serrata.

METHODS

The larvae were collected from the sheep's visceral organs from the Yazd Slaughterhouse during May 2013. DNA extraction was done in three different methods, including commercial DNA extraction kit, Phenol Chloroform Isoamylalcohol (PCI), and salting out. Extracted DNA in each method was assessed for quantity and quality using spectrophotometery and agarose gel electrophoresis, respectively.

RESULTS

The less duration was regarding to commercial DNA extraction kit and then salting out protocol. The cost benefit one was salting out and then PCI method. The best quantity was regarding to PCI with 72.20±29.20 ng/μl, and purity of OD260/OD280 in 1.76±0.947. Agarose gel electrophoresis for assessing the quality found all the same.

CONCLUSION

Salting out is introduced as the best method for DNA extraction from L. seratta as a food-borne pathogen with the least costand appropriate purity. Although, the best purity was regarding to PCI but PCI is not safe as salting out. In addition, the duration of salting out was less than PCI. The least duration was seen in commercial DNA extraction kit, but it is expensive and therefore is not recommended for developing countries where consumption of offal is common.

Analysis of the complete mitochondrial genome sequence of Palinura homarus

The complete mitochondrial genome sequence of Palinura homarus was obtained using PCR amplification and walking sequencing (GenBank accession no. JN_542716). The complete mitochondrial genome of P. homarus was 15,665 bp long and showed significant AT bias (67% AT content, 33% GC content). The A + T-rich region included copy-related control information and a poly (dT) structure that related to replication and transcription. In this study, the gene arrangement was consistent with other Palinura mitochondrial genomes and the sequence was strikingly similar to Panulirus ornatus, which would be useful in species identification and natural resources conservation.

Complete Mitochondrial Genome of a Tongue Worm Armillifer agkistrodontis

Armillifer agkistrodontis (Ichthyostraca: Pantastomida) is a parasitic pathogen, only reported in China, which can cause a zoonotic disease, pentastomiasis. A complete mitochondrial (mt) genome was 16,521 bp comprising 13 protein-coding genes (PCGs), 22 tRNA genes, 2 rRNA genes, and 1 non-coding region (NCR). A phylogenetic tree drawn with the concatenated amino acid sequences of the 6 conserved PCGs (atp6, cox1-3, and nad2) showed that A. agkistrodontis and Armillifer armillatus constituted a clade Pentastomida which was a sister group of the Branchiura. The complete mt genome sequence of A. agkistrodontis provides important genetic markers for both phylogenetic and epidemiological studies of pentastomids.

Armillifer armillatus infestation in Human; public health scenario of a snake parasite: a report of three cases

We report cases of Armillifer Armillatus infestation in three Nigerian adults within two and half years in our health facility. The first patient was a 70 year old farmer and a regular consumer of snake meat for over 50 years. He presented in February, 2014 for follow-up visit as he was a known systemic hypertensive patient. He was incidentally discovered to have multiple comma-shaped calcific lesions in the lungs and liver on a chest radiograph. These were better demonstrated on abdominal ultrasound and computed tomographic scans. He was asymptomatic. The second patient was a 42 year old male civil servant who presented in December 2015 with dry cough and right loin pain for five and three days respectively. His past medical history revealed that he had been treated previously for pneumonia. He has never eaten snake meat but consumed Alligator (Amphibious reptile) for many years but stopped about 12 years ago. Similar calcific lesions were also noted in his liver and lung parenchyma on chest radiograph and abdominal ultrasound scan. The third patient was an 80 year old man who presented in April, 2014 with dizziness and diminished urine output of one day duration. He was a farmer who has been consuming snake meat for many years, and has been on management for systemic arterial hypertension and prostatic hypertrophy. Chest radiograph and abdomino-pelvic ultrasound incidentally revealed multiple comma-shaped calcific lesions in the lungs and liver. The liver function test parameters were all within normal limits but the electrolytes were deranged and he was anaemic with a Packed Cell Volume of 27%. A diagnosis of Armillifer Armillatus infestation was made in these patients, and they were conservatively managed with Mebendazole. The third case was catherized and the deranged electrolytes were corrected. The first patient was lost to follow-up, whiles the second and third had no remarkable symptoms on subsequent follow-up visits.

Molecular characterization of Argulus bengalensis and Argulus siamensis (Crustacea: Argulidae) infecting the cultured carps in West Bengal, India using 18S rRNA gene sequences

The present study characterized Argulus spp. infecting the cultured carps using 18S rRNA gene sequences, estimated the genetic similarity among Argulus spp. and established their phylogenetic relationship. Of the 320 fish samples screened, 34 fish (10.6%) had Argulus infection. The parasitic frequency index (PFI) was observed to be high (20%) in Hypophthalmichthys molitrix and Labeo bata. The frequency of infection was high in September (PFI: 17%) and October (PFI: 12.9%). The 18S rRNA sequences of five A. bengalensis (KF583878, KF192316, KM016968, KM016969, and KM016970) and one A. siamensis (KF583879) of this study showed genetic heterogeneity and exhibited 77-99% homology among the 18S rRNA gene sequences of Argulus spp. of NCBI GenBank database. Among the Indian Argulus spp. the sequence homology was 87-100%. Evolutionary pair-wise distances between Indian Argulus spp. and other Argulus spp. ranged from 0 to 20.20%. In the phylogenetic tree, all the crustaceans were clustered together as a separate clade with two distinct lineages. The lineage-1 comprised exclusive of Branchiura (Argulus spp.). All Argulus bengalensis clustered together and A. siamensis (KF583879) was closely related to Argulus sp. JN558648. The results of the present study provided baseline data for future work on population structure analysis of Indian Argulus species.

The complete mitochondrial genome of the semiterrestrial crab, Chiromantes neglectum (Eubrachyura: Grapsoidea: Sesarmidae)

The complete mitogenome of the semiterrestrial crab Chiromantes neglectum was sequenced. It contained the entire set of 37 genes. The gene order was basically identical to pancrustacean ground pattern, except for thetrnH and trnQ genes. Phylogenetic inferences based on protein-coding genes (PCGs) provide strong evidence that places C. neglectum within an intermingled 'Grapsoidea & Ocypodoidea' clade.

Finding of pentastomes of genus Reighardia (Pentastomida) in the Belcher's gull (Larus belcheri)

This report describes the finding of Reighardia sp. (Pentastomida) infecting the air sac of two Belcher's gulls (Larus belcheri) found dead on the beaches of Pucusana, a district in southern Lima, Peru. Three pentastomes were collected from two Belcher's gulls. Then, they were morphologically and molecular analyzed. Molecular characterization of the parasite was achieved by amplifying a fragment of the small subunit ribosomal RNA gene (SSU rRNA). Based on both morphological and molecular data the pentastomes were identified as pentastomes of the genus Reighardia. This is the first report showing that the Belcher's gull is a new natural definitive host for this pentastome.

Fossils of parasites: what can the fossil record tell us about the evolution of parasitism?

Parasites are common in many ecosystems, yet because of their nature, they do not fossilise readily and are very rare in the geological record. This makes it challenging to study the evolutionary transition that led to the evolution of parasitism in different taxa. Most studies on the evolution of parasites are based on phylogenies of extant species that were constructed based on morphological and molecular data, but they give us an incomplete picture and offer little information on many important details of parasite-host interactions. The lack of fossil parasites also means we know very little about the roles that parasites played in ecosystems of the past even though it is known that parasites have significant influences on many ecosystems. The goal of this review is to bring attention to known fossils of parasites and parasitism, and provide a conceptual framework for how research on fossil parasites can develop in the future. Despite their rarity, there are some fossil parasites which have been described from different geological eras. These fossils include the free-living stage of parasites, parasites which became fossilised with their hosts, parasite eggs and propagules in coprolites, and traces of pathology inflicted by parasites on the host's body. Judging from the fossil record, while there were some parasite-host relationships which no longer exist in the present day, many parasite taxa which are known from the fossil record seem to have remained relatively unchanged in their general morphology and their patterns of host association over tens or even hundreds of millions of years. It also appears that major evolutionary and ecological transitions throughout the history of life on Earth coincided with the appearance of certain parasite taxa, as the appearance of new host groups also provided new niches for potential parasites. As such, fossil parasites can provide additional data regarding the ecology of their extinct hosts, since many parasites have specific life cycles and transmission modes which reflect certain aspects of the host's ecology. The study of fossil parasites can be conducted using existing techniques in palaeontology and palaeoecology, and microscopic examination of potential material such as coprolites may uncover more fossil evidence of parasitism. However, I also urge caution when interpreting fossils as examples of parasites or parasitism-induced traces. I point out a number of cases where parasitism has been spuriously attributed to some fossil specimens which, upon re-examination, display traits which are just as (if not more) likely to be found in free-living taxa. The study of parasite fossils can provide a more complete picture of the ecosystems and evolution of life throughout Earth's history.

Molecular Diagnosis of Abdominal Armillifer grandis Pentastomiasis in the Democratic Republic of Congo

Pentastomiasis is an emerging snake-borne parasitic zoonosis in the tropics. We describe a molecular and morphological study to diagnose a cluster of asymptomatic abdominal human infections caused by Armillifer grandis. The findings may indicate a silent epidemic in a rural area where severe symptomatic ocular cases with the same parasite species have recently surfaced. Molecular diagnostics are of increasing importance when patient material from remote areas cannot be thoroughly examined locally for logistic reasons.

Mitochondrial genome of the intertidal acorn barnacle Tetraclita serrata Darwin, 1854 (Crustacea: Sessilia): Gene order comparison and phylogenetic consideration within Sessilia

The complete mitochondrial genome of the intertidal barnacle Tetraclita serrata Darwin, 1854 (Crustacea: Maxillopoda: Sessilia) is presented. The genome is a circular molecule of 15,200 bp, which encodes 13 PCGs, 2 ribosomal RNA genes, and 22 transfer RNA genes. All non-coding regions are 591 bp in length, with the longest one speculated as the control region (389 bp), which is located between srRNA and trnK. The overall A+T content of the mitochondrial genome of T. serrata is 65.4%, which is lowest among all the eight mitochondrial genomes reported from sessile barnacles. There are variations of initiation and stop codons in the reported sessile barnacle mitochondrial genomes. Large-scale gene rearrangements are found in these genomes as compared to the pancrustacean ground pattern. ML and Bayesian analyses of all 15 complete mitochondrial genomes available from Maxillopoda lead to identical phylogenies. The phylogenetic tree based on mitochondrial PCGs shows that Argulus americanus (Branchiura) cluster with Armillifer armillatus (Pentastomida), distinct from all ten species from Cirripedia. Within the order Sessilia, Amphibalanus amphitrite (Balanidae) clusters with Striatobalanus amaryllis (Archaeobalanidae), and Nobia grandis (Pyrgomatidae). However, the two Megabalanus (Balanidae) are separated from the above grouping, resulting in non-monophyly of the family Balanidae. Moreover, the two Megabalanus have large-scale rearrangements as compared to the gene order shared by former three species. Therefore, both phylogenetic analysis using PCG sequences and gene order comparison suggest that Balanidae is not a monophyletic group. Given the limited taxa and moderate support values of the internal branches, the non-monophyly of the family Balanidae requires further verification.

Complete mitochondrial genome of Epigonichthys cultellus (Cephalochordata: Branchiostomatidae)

The systematic position of the amphioxus species with dextral gonads distributed in the southeast Pacific and Indian Oceans has remained to be clarified due to the adoption of different names by different researchers. Mitochondrial (mt) DNA is generally considered to be a powerful molecular marker in taxonomic studies. For a reliable systematics of the amphioxi collected from the South China Sea, we sequenced the complete mtDNA from a single specimen and compared it with those of the other eight amphioxus species. The present mtDNA genome contains 13 protein-coding genes, 22 tRNA genes, and two rRNA genes, with the same gene order as those in Branchiostoma and Epigonichthys, which, however, is different from that in Asymmetron. Based on our morphological data (including measurements of some characters) and the features of the mt genome, together with the distribution records of the dextral-gonad amphioxi, we conclude that the present mtDNA sequence most likely represents that of E. cultellus.

Control region sequences indicate that multiple externae represent multiple infections by Sacculina carcini (Cirripedia: Rhizocephala)

The rhizocephalan barnacle, Sacculina carcini, is a common parasite of the European shore crab, Carcinus maenas, in which it causes significant detrimental physical and behavioral modifications. In the vast majority of cases, the external portion of the parasite is present in the form of a single sac-like externa; in rare cases, double or even triple externae may occur on the same individual host. Here, we use a highly variable DNA marker, the mitochondrial control region (CR), to investigate whether multiple externae in S. carcini represent infection by multiple parasites or asexual cloning developed by a single parasite individual. Sequences for multiple externae from C. maenas hosts from the Danish inlet, Limfjorden, and from the mud flates at Roscoff, France, were compared. In almost all cases, double or triple externae from an individual host yielded different haplotypes. In the few cases where identical haplotypes were identified from externae on a multiple-infected host, this always represented the most commonly found haplotype in the population. This indicates that in Sacculina carcini, the presence of multiple externae on a single host reflects infection by different individual parasites. A haplotype network of CR sequences also suggests a degree of geographical partitioning, with no shared haplotypes between the Limfjorden and Roscoff. Our data represent the first complete CR sequences for a rhizocephalan, and a unique gene order was also revealed. Although the utility of CR sequences for population-level work must be investigated further, the CR has proved a simple to use and highly variable marker for studies of S. carcini and can easily be applied to a variety of studies in this important parasite.

First insights on the mitochondrial genetic variability of Lightiella magdalenina (Crustacea), the sole Mediterranean cephalocarid species

BACKGROUND

Here we report the first insight into the mitochondrial (Cytochrome c Oxidase subunit I - COI and Cytochrome b - Cyt b) genetic variation of the only Mediterranean cephalocarid Lightiella magdalenina.

FINDINGS

COI sequences provide a scenario of low intraspecific variability, while significant genetic divergence occurs between L. magdalenina and L. incisa. Interestingly, Cyt b sequences reveal a higher degree of intraspecific variability, with no shared haplotypes between the sites considered.

CONCLUSIONS

In the future, COI and Cyt b molecular markers could be used as valuable tools to shed new light into the extant species within the genus Lightiella thus providing molecular support to the taxonomical identifications carried out on a morphological basis.

The mitochondrial genome of the plant bug Apolygus lucorum (Hemiptera: Miridae): Presently known as the smallest in Heteroptera

The complete mitochondrial (mt) genome of the plant bug, Apolygus lucorum, an important cotton pest, has been sequenced and annotated in this study. The entire circular genome is 14 768 bp in size and represents the smallest in presently known heteropteran mt genomes. The mt genome is encoding for two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, 13 protein coding genes and a control region, and the order, content, codon usage and base organization show similarity to a great extent to the hypothetical ancestral model. All protein coding genes use standard initiation codons ATN. Conventional stop codons TAA and TAG have been assigned to the most protein coding genes; however, COIII, ND4 and ND5 genes show incomplete terminator signal (T). All tRNA genes possess the typical clover leaf structure, but the dihydrouridine arm of tRNA(Ser(AGN)) only forms a simple loop. Secondary structure models of rRNA genes are generally in accordance with the former models, although some differences exist in certain parts. Three intergenic spacers have never been found in sequenced mt genomes of Heteroptera. The phylogenetic study based on protein coding genes is largely congruent with previous phylogenetic work. Both Bayesian inference and maximum likelihood analyses highly support the sister relationship of A. lucorum and Lygus lineolaris, and Miridae presents a sister position to Anthocoridae.

Prospects of using DNA barcoding for species identification and evaluation of the accuracy of sequence databases for ticks (Acari: Ixodida)

Ticks are important vectors of disease and parasites of livestock. Species identification of ticks has been traditionally based on morphological characters, which is usually limited by the condition of samples and little variation among specimens, so a rapid and reliable identification method is needed. DNA barcoding uses a standard fragment of the mitochondrial gene cytochrome oxidase c subunit I (COI) to identify species and has been successfully used in many taxa. In this study, we applied DNA barcoding to tick species. K2P distances showed that most interspecific divergences exceed 8%, while intraspecific distances were usually lower than 2%. However, intraspecific distances of 12 species were unexpectedly high. ABGD grouping results demonstrated that sequences of these species should be divided into 2 or more groups. And some exceptional clustering occurred among sequences of Hyalomma marginatum, Hy. truncatum, and Hy. dromedarii, Amblyomma testudinarium and A. pattoni, Rhipicephalus sanguineus and R. pumilio, Haemaphysalis parva and Ha. concinna, Ixodes asanumai and I. nipponensis. Additionally, 226 unnamed sequences were assigned to known species or constituted different groups, and K2P distances of all these groups were less than 2%. In conclusion, our study demonstrated that DNA barcoding is a useful tool for the identification of tick species, and further work is needed to reveal ambiguous species delimitation in some problematic genera.

Imported Armillifer pentastomiasis: report of a symptomatic infection in The Netherlands and mini-review

We report a case of symptomatic visceral Armillifer pentastomiasis in a 23-year-old female Liberian immigrant to The Netherlands. The patient was referred to the gynecologist because of lower abdominal pain. During laparotomy, multiple adhesions were seen in the lower pelvis and a hydrosalpinx with an encapsulated Armillifer nymph, most likely Armillifer armillatus, was found. Key features of the parasite's cuticle which facilitate the diagnosis of pentastomiasis, are presented. Symptomatic pentastomiasis is uncommon, and most cases are diagnosed incidentally during surgery for other reasons, or at autopsy. With regard to increasing international migration, other imported pentastomiasis cases to Europe and North America are reviewed, and more cases are likely to be seen in the future.

The complete mitochondrial genomes of three parasitic nematodes of birds: a unique gene order and insights into nematode phylogeny

Background

Analyses of mitochondrial (mt) genome sequences in recent years challenge the current working hypothesis of Nematoda phylogeny proposed from morphology, ecology and nuclear small subunit rRNA gene sequences, and raise the need to sequence additional mt genomes for a broad range of nematode lineages.

Results

We sequenced the complete mt genomes of three Ascaridia species (family Ascaridiidae) that infest chickens, pigeons and parrots, respectively. These three Ascaridia species have an identical arrangement of mt genes to each other but differ substantially from other nematodes. Phylogenetic analyses of the mt genome sequences of the Ascaridia species, together with 62 other nematode species, support the monophylies of seven high-level taxa of the phylum Nematoda: 1) the subclass Dorylaimia; 2) the orders Rhabditida, Trichinellida and Mermithida; 3) the suborder Rhabditina; and 4) the infraorders Spiruromorpha and Oxyuridomorpha. Analyses of mt genome sequences, however, reject the monophylies of the suborders Spirurina and Tylenchina, and the infraorders Rhabditomorpha, Panagrolaimomorpha and Tylenchomorpha. Monophyly of the infraorder Ascaridomorpha varies depending on the methods of phylogenetic analysis. The Ascaridomorpha was more closely related to the infraorders Rhabditomorpha and Diplogasteromorpha (suborder Rhabditina) than they were to the other two infraorders of the Spirurina: Oxyuridorpha and Spiruromorpha. The closer relationship among Ascaridomorpha, Rhabditomorpha and Diplogasteromorpha was also supported by a shared common pattern of mitochondrial gene arrangement.

Conclusions

Analyses of mitochondrial genome sequences and gene arrangement has provided novel insights into the phylogenetic relationships among several major lineages of nematodes. Many lineages of nematodes, however, are underrepresented or not represented in these analyses. Expanding taxon sampling is necessary for future phylogenetic studies of nematodes with mt genome sequences.

Evolution and phylogeny of the mud shrimps (Crustacea: Decapoda) revealed from complete mitochondrial genomes

BACKGROUND

The evolutionary history and relationships of the mud shrimps (Crustacea: Decapoda: Gebiidea and Axiidea) are contentious, with previous attempts revealing mixed results. The mud shrimps were once classified in the infraorder Thalassinidea. Recent molecular phylogenetic analyses, however, suggest separation of the group into two individual infraorders, Gebiidea and Axiidea. Mitochondrial (mt) genome sequence and structure can be especially powerful in resolving higher systematic relationships that may offer new insights into the phylogeny of the mud shrimps and the other decapod infraorders, and test the hypothesis of dividing the mud shrimps into two infraorders.

RESULTS

We present the complete mitochondrial genome sequences of five mud shrimps, Austinogebia edulis, Upogebia major, Thalassina kelanang (Gebiidea), Nihonotrypaea thermophilus and Neaxius glyptocercus (Axiidea). All five genomes encode a standard set of 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes and a putative control region. Except for T. kelanang, mud shrimp mitochondrial genomes exhibited rearrangements and novel patterns compared to the pancrustacean ground pattern. Each of the two Gebiidea species (A. edulis and U. major) and two Axiidea species (N. glyptocercus and N. thermophiles) share unique gene order specific to their infraorders and analyses further suggest these two derived gene orders have evolved independently. Phylogenetic analyses based on the concatenated nucleotide and amino acid sequences of 13 protein-coding genes indicate the possible polyphyly of mud shrimps, supporting the division of the group into two infraorders. However, the infraordinal relationships among the Gebiidea and Axiidea, and other reptants are poorly resolved. The inclusion of mt genome from more taxa, in particular the reptant infraorders Polychelida and Glypheidea is required in further analysis.

CONCLUSIONS

Phylogenetic analyses on the mt genome sequences and the distinct gene orders provide further evidences for the divergence between the two mud shrimp infraorders, Gebiidea and Axiidea, corroborating previous molecular phylogeny and justifying their infraordinal status. Mitochondrial genome sequences appear to be promising markers for resolving phylogenetic issues concerning decapod crustaceans that warrant further investigations and our present study has also provided further information concerning the mt genome evolution of the Decapoda.