Characterization of the complete mitochondrial genome of the swine kidney worm Stephanurus dentatus (Nematoda: Syngamidae) and phylogenetic implications

The complete mitochondrial genome of the rodent flea Nosopsyllus laeviceps: genome description, comparative analysis, and phylogenetic implications

BACKGROUND

Fleas are one of the most common and pervasive ectoparasites worldwide, comprising at least 2500 valid species. They are vectors of several disease-causing agents, such as Yersinia pestis. Despite their significance, however, the molecular genetics, biology, and phylogenetics of fleas remain poorly understood.

METHODS

We sequenced, assembled, and annotated the complete mitochondrial (mt) genome of the rodent flea Nosopsyllus laeviceps using next-generation sequencing technology. Then we combined the new mitogenome generated here with mt genomic data available for 23 other flea species to perform comparative mitogenomics, nucleotide diversity, and evolutionary rate analysis. Subsequently, the phylogenetic relationship within the order Siphonaptera was explored using the Bayesian inference (BI) and maximum likelihood (ML) methods based on concentrated data for 13 mt protein-coding genes.

RESULTS

The complete mt genome of the rodent flea N. laeviceps was 16,533 base pairs (bp) in a circular DNA molecule, containing 37 typical genes (13 protein-coding genes, 22 transfer RNA [tRNA] genes, and two ribosomal RNA [rRNA] genes) with one large non-coding region (NCR). Comparative analysis among the order Siphonaptera showed a stable gene order with no gene arrangement, and high AT content (76.71-83.21%) with an apparent negative AT and GC skew except in three fleas Aviostivalius klossi bispiniformis, Leptopsylla segnis, and Neopsylla specialis. Moreover, we found robust evidence that the cytochrome c oxidase subunit 1 (cox1) gene was the most conserved protein-coding gene (Pi = 0.15, non-synonymous/synonymous [Ka/Ks] ratio = 0.13) of fleas. Phylogenomic analysis conducted using two methods revealed different topologies, but both results strongly indicated that (i) the families Ceratophyllidae and Leptopsyllidae were paraphyletic and were the closest to each other, and (ii) the family Ctenophthalmidae was paraphyletic.

CONCLUSIONS

In this study, we obtained a high-quality mt genome of the rodent flea N. laeviceps and performed comparative mitogenomics and phylogeny of the order Siphonaptera using the mt database. The results will enrich the mt genome data for fleas, lay a foundation for the phylogenetic analysis of fleas, and promote the evolutionary analysis of Siphonaptera.

Molecular phylogeny of the family Rhabdiasidae (Nematoda: Rhabditida), with morphology, genetic characterization and mitochondrial genomes of Rhabdias kafunata and R. bufonis

BACKGROUND

The family Rhabdiasidae (Nematoda: Rhabditida) is a globally distributed group of nematode parasites, with over 110 species parasitic mainly in amphibians and reptiles. However, the systematic position of the family Rhabdiasidae in the order Rhabditida remains unsolved, and the evolutionary relationships among its genera are still unclear. Moreover, the present knowledge of the mitochondrial genomes of rhabdiasids remains limited.

METHODS

Two rhabdiasid species: Rhabdias kafunata Sata, Takeuchi & Nakano, 2020 and R. bufonis (Schrank, 1788) collected from the Asiatic toad Bufo gargarizans Cantor (Amphibia: Anura) in China, were identified based on morphology (light and scanning electron microscopy) and molecular characterization (sequencing of the nuclear 28S and ITS regions and mitochondrial cox1 and 12S genes). The complete mitochondrial genomes of R. kafunata and R. bufonis were also sequenced and annotated for the first time. Moreover, phylogenetic analyses based on the amino acid sequences of 12 protein-coding genes (PCGs) of the mitochondrial genomes were performed to clarify the systematic position of the family Rhabdiasidae in the order Rhabditida using maximum likelihood (ML) and Bayesian inference (BI). The phylogenetic analyses based on the 28S + ITS sequences, were also inferred to assess the evolutionary relationships among the genera within Rhabdiasidae.

RESULTS

The detailed morphology of the cephalic structures, vulva and eggs in R. kafunata and R. bufonis was revealed using scanning electron microscopy (SEM) for the first time. The characterization of 28S and ITS regions of R. kafunata was reported for the first time. The mitogenomes of R. kafunata and R. bufonis are 15,437 bp and 15,128 bp long, respectively, and both contain 36 genes, including 12 PCGs (missing atp8). Comparative mitogenomics revealed that the gene arrangement of R. kafunata and R. bufonis is different from all of the currently available mitogenomes of nematodes. Phylogenetic analyses based on the ITS + 28S data showed Neoentomelas and Kurilonema as sister lineages, and supported the monophyly of Entomelas, Pneumonema, Serpentirhabdias and Rhabdias. Mitochondrial phylogenomic results supported Rhabdiasidae as a member of the superfamily Rhabditoidea in the suborder Rhabditina, and its occurrance as sister to the family Rhabditidae.

CONCLUSIONS

The complete mitochondrial genome of R. kafunata and R. bufonis were reported for the first time, and two new gene arrangements of mitogenomes in Nematoda were revealed. Mitogenomic phylogenetic results indicated that the family Rhabdiasidae is a member of Rhabditoidea in Rhabditina, and is closely related to Rhabditidae. Molecular phylogenies based on the ITS + 28S sequence data supported the validity of Kurilonema, and showed that Kurilonema is sister to Neoentomelas. The present phylogenetic results also indicated that the ancestors of rhabdiasids seem to have initially infected reptiles, then spreading to amphibians.

The complete chloroplast genome of Holarrhena pubescens and its phylogenetic analysis

Holarrhena pubescens Wall. ex G. Don, 1837 is an important medicinal plant belonging to the Holarrhena genus in the Apocynaceae family. In this study, the complete chloroplast (cp) genome sequence of H. pubescens was sequenced using the Illumina NovaSeq platform. The cp genome of H. pubescens was 160,108 bp in length with 37.21% overall GC content. The cp genome of H. pubescens containing a large single-copy region (LSC, 88,685 bp), a small single-copy region (SSC, 18,671 bp), and a pair of inverted repeat regions (SSC, 26,376 bp). The cp genome encoded 129 genes, including 84 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Phylogenetic analysis based on complete protein coding genes sequences revealed that H. pubescens was closest to Beaumontia murtonii.

A Novel Mitochondrial Genome Fragmentation Pattern in the Buffalo Louse Haematopinus tuberculatus (Psocodea: Haematopinidae)

Sucking lice are obligate ectoparasites of mammalian hosts, causing serious public health problems and economic losses worldwide. It is well known that sucking lice have fragmented mitochondrial (mt) genomes, but many remain undetermined. To better understand patterns of mt genome fragmentation in the sucking lice, we sequenced the mt genome of the buffalo louse Haematopinus tuberculatus using next-generation sequencing (NGS). The mt genome of H. tuberculatus has ten circular minichromosomes containing a total of 37 genes. Each minichromosome is 2.9-5.0 kb long and carries one to eight genes plus one large non-coding region. The number of mt minichromosomes of H. tuberculatus (ten) is different from those of congeneric species (horse louse H. asini, domestic pig louse H. suis and wild pig louse H. apri) and other sucking lice. Two events (gene translocation and merger of mt minichromosome) are observed in Haematopinus. Compared to other studies, our phylogeny generated from mt genome datasets showed a different topology, suggesting that inclusion of data other than mt genomes would be required to resolve phylogeny of sucking lice. To our knowledge, this is the first report of a ten mt minichromosomes genome in sucking lice, which opens a new outlook into unexplored mt genome fragmentation patterns in sucking lice.