Gene arrangement, phylogeny and divergence time estimation of mitogenomes in Thrips

The metazoan mitogenomes usually display conserved gene arrangement while thrips are known for their extensive gene rearrangement, and duplication of the control region.

We sequenced complete mitogenomes of eight species of thrips to determine the gene arrangement, phylogeny and divergence time estimation. All contain 37 genes and one control region, (CR) except four species with two CRs. Duplicated tRNAs were detected in Mycterothrips nilgiriensis and Thrips florum. nad4-nad4L were not found adjacent to each other in Phibalothrips peringueyi and Plicothrips apicalis. Both Bayesian and likelihood phylogenetic analyses of thrips mitogenomes supported the monophyly of two suborders (Terebrantia and Tubulifera) and the two largest families (Phlaeothripidae and Thripidae). Out of seven earlier proposed ancestral gene blocks, six are conserved in Panchaetothripinae, three in Thripinae and two in Phlaeothripidae. Additionally, eight Thrips Gene Blocks were identified, of which, three conserved in Tubulifera, four in Terebrantia, and one only in Aeolothripidae. Forty-two gene boundaries (15 from previous study + 27 new) were identified. The molecular divergence time is estimated for the order Thysanoptera and suggested that these insects may have been diversified from hemipterans in the late Permian period. The most recent ancestors belong to family Thripidae and Phlaeothripidae, which were diversified in upper Cretaceous period and showed higher rates of rearrangement from the ancestral gene order.

The current study is the first largest effort to provide the new insights into the mitogenomic features, gene arrangement, phylogeny and divergence time estimation of thrips belonging to the order Thysanoptera.