Complete mitochondrial genome of Penicillidia dufourii (Diptera: Hippoboscoidea: Nycteribiidae) and phylogenetic relationship.

Penicillidia dufourii (Westwood 1834) is a specialized parasite categorized under family Nycteribiidae that prefers to parasitize the body surface of various bats under the genus Myotis. Many species of the family Nycteribiidae are carriers of various pathogens; however, research on P. dufourii remains scarce, and studies on its molecular identification and population genetic structure are still lacking. In this study, the complete mitochondrial genome of P. dufourii was elucidated for the first time using Illumina sequencing. The mitochondrial genome is 15,354 bp in size and encodes approximately 37 genes, including 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 1 control region. Analysis of 13 protein-coding genes revealed that UUA, UCA, CGA, and GGA were the most common codons, while nad4L had the fastest evolutionary rate and cox1 the slowest. Phylogenetic analysis based on the mitochondrial genome indicated that P. dufourii is clustered with other species of the family Nycteribiidae and is most closely related to Nycteribia parvula and Phthiridium szechuanum.

The first mitogenome of the subfamily Stenoponiinae (Siphonaptera: Ctenophthalmidae) and implications for its phylogenetic position.

Fleas are the most important insect vectors that parasitize warm-blooded animals and are known vectors of zoonotic pathogens. A recent study showed that Stenoponia polyspina parasitizing Eospalax baileyi in Zoige County have carried Bartonella spp. and Spotted fever group Rickettsia (SFGR). Accurate identification and differentiation of fleas are essential for prevention and control of zoonotic pathogens. To understand phylogenetic relationship of the subfamily Stenoponiinae, we described morphological characteristics of S. polyspina and sequenced its mitogenome with 14,933 bp in size and high A + T content (~ 79%). The S. polyspina mitogenome retained the ancestral pattern of mitochondrial gene arrangement of arthropods without rearrangement. The start codons of 13 protein-coding genes (PCGs) are traditional ATN and the stop codons are TAA or TAG. Anticodon loops of all tRNA genes were 7 bp except for trnL2 and trnD had anticodon loops with 9 bp and the abnormal anticodon loops may be associated with frameshifting mutation. Genetic distance and Ka/Ks ratios indicated that all 13 PCGs of S. polyspina were subjected to purifying selection, with cox1 at the slowest rate and atp8 at the fastest rate. The mitogenomes of 24 species representing 7 families in the order Siphonaptera were selected to reconstruct phylogenetic tree based on concatenated nucleotide sequences of two datasets (PCGRNA matrix and PCG12RNA matrix) using Bayesian inference (BI) and Maximum likelihood (ML) methods. Phylogenetic tree supported that the superfamilies Ceratophylloidea, Vermipsylloidea, Pulicoidea were monophyletic and the superfamily Hystrichopsylloidea was paraphyletic. The family Ctenophthalmidae was monophyletic in PCGRNA-ML (codon partition) and paraphyletic in the remain trees. S. polyspina belongs to the subfamily Stenoponiinae was closely more related to the subfamily Rhadinopsyllinae. This paper explored phylogenetic position of diverse clades within the order Siphonaptera based on morphological and mitogenome data of S. polyspina. Our research enriched NCBI database of the order Siphonaptera.