The first divergence time estimation of the subfamily Stenogastrinae (Hymenoptera: Vespidae) based on mitochondrial phylogenomics

Abstract

In this study, the mitochondrial genomes of three Stenogastrinae species, Eustenogaster scitula, Liostenogaster nitidipennis and Parishnogaster mellyi were sequenced and annotated, and a total of 16 vespid mtgenomes are comparatively analyzed. Our results indicate that codon usage bias is mainly dominated by mutational pressure, and affected only slightly by natural selection. Selective pressure analysis of protein-coding genes (PCGs) shows that the highest evolutionary rate is present in NADH complex I, and the lowest in cox1. Compared with the reported mtgenomes of other Vespidae, in Stenogastrinae, trnH is shifted to a new position. Phylogenetic analyses are performed using Bayesian method and Maximum Parsimony. Phylogenetic analysis further confirms that the Stenogastrinae is the sister group of all remaining Vespidae. Divergence time of Stenogastrinae from other Vespidae is estimated at ~ 166 Mya. Our results also support that eusociality evolved twice in the family Vespidae.

Keywords: Divergence time; Hymenoptera; Mitochondrial genome; Phylogeny; Stenogastrinae.