Gene rearrangement in the mitogenome of whiteflies (Hemiptera: Aleyrodinae) along with their phylogeny and characterization of complete mitogenome of Aleurodicus rugioperculatus

Background Mitochondrial genome rearrangements have been used for defining evolutionary relationships, but there have been incidences of convergences at different taxonomic levels. To gain new insights into whitefly mitogenome evolution, gene arrangement and phylogeny, we sequenced the complete mitogenome of Aleurodicus rugioperculatus (Aleyrodidae: Aleurodicinae) and conducted a comparative analysis with the previously published mitogenomes across the family Aleyrodidae. Methods and Results The complete mitogenome of Aleurodicus rugioperculatus was generated by Next generation sequencing method. It is 15,060 bp circular molecule with 86.5 A + T and 5.5% G + C content. It contains 37 genes (13 PCGs, two rRNAs, 22 tRNAs) and a non-coding control region (CR). Comparative analysis of codon usage indicated that the subfamily Aleyrodinae have weaker bias than Aleurodicinae. Bayesian Inference (BI) and Maximum Likelihood (ML) phylogenetic analyses yielded similar topologies supporting the monophyly of Aleyrodinae and Aleurodicinae. The gene order of 13 whiteflies was compared with ancestor to examine the plesiomorphies, synapomorphies and autapomorphies. We identified five gene blocks (I-V) in the whitefly ancestor that are shared plesiomorphies retained in different whitefly lineages. Gene block I is conserved in all whiteflies except three species. Conversely, we detected 83 derived gene boundaries within the family. Mapping these gene boundaries onto a phylogenetic tree revealed that 16 were symplesiomorphies for two subfamilies; 9 were synapomorphies between the species, and 28 autapomorphies for individual species. Conclusions Comparative analyses of gene order of whiteflies revealed the derived gene boundaries which can be further investigated with more mitogenome data to examine the genome evolution in whiteflies.