The mitochondrial genomes of Muscidifurax similadanacus, M. sinesensilla, Nasonia vitripennis, and Pachycrepoideus vindemmiae were sequenced to better understand the structural evolution of Pteromalidae mitogenomes. These newly sequenced mitogenomes all contained 37 genes. Nucleotide composition was AT-biased and the majority of the protein-coding genes exhibited a negative AT skew. All 13 protein-coding genes (PCGs) initiated with the standard start codon of ATN, excepted for nad1 of N. vitripennis, which started with TTG, and terminated with a typical stop codon TAA/TAG or an incomplete stop codon T. All transfer RNA (tRNA) genes were predicted to fold into the typical clover-leaf secondary structures, except for trnS1, which lacks the DHU arm in all species. In P. vindemmiae, trnR and trnQ lack the DHU arm and TΨC arm, respectively. Although most genes evolved under a strong purifying selection, the Ka/Ks value of the atp8 gene of P. vindemmiae was greater than 1, indicating putative positive selection. A novel transposition of trnR in P. vindemmiae was revealed, which was the first of this kind to be reported in Pteromalidae. Two kinds of datasets (PCG12 and AA) and two inference methods (maximum likelihood and Bayesian inference) were used to reconstruct a phylogenetic hypothesis for the newly sequenced mitogenomes of Pteromalidae and those deposited in GenBank. The topologies obtained recovered the monophyly of the three subfamilies included. Pachyneurinae and Pteromalinae were recovered as sister families, and both appeared sister to Sycophaginae. The pairwise breakpoint distances of mitogenome rearrangements were estimated to infer phylogeny among pteromalid species. The topology obtained was not totally congruent with those reconstructed using the ML and BI methods.
Keywords: Pteromalidae; gene rearrangement; mitogenome; parasitic lifestyles; phylogenetic position.