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

Abstract

Abstract 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. Graphical Abstract