Purifying selection drove the adaptation of mitochondrial genes along with correlation of gene rearrangements and evolutionary rates in two subfamilies of Whitefly (Insecta: Hemiptera).

Insect mitochondrial genomes (mitogenomes) are usually represented by a conserved gene order. Whiteflies exhibit gene rearrangement in their mitogenomes; however, understanding how nucleotide substitution rates shape gene rearrangement in whiteflies is unclear due to the limited number of mitogenomes. Additionally, the mechanisms by which selection pressure drives adaptations in mitochondrial genes in the two subfamilies of whiteflies are not yet known. Here, we analyzed 18 whitefly mitogenomes, including one newly generated mitogenome, to compare nucleotide substitution rates, selection pressure, and gene arrangements. The newly generated mitogenome is reported along with reannotation of Pealius mori and comparisons to other whitefly mitogenomes. Comparative studies on nucleotide composition of 18 whiteflies revealed the positive GC skewness, confirming the reversal of strand asymmetry. We found 11 rearranged gene orders within two subfamilies of whiteflies with 8-18 breakpoints of gene rearrangements. Members of the subfamily Aleyrodinae exhibit more complex pathways in the evolution of gene order as compared to the subfamily Aleurodicinae. Our findings also revealed that the increase or reduction of nucleotide substitution rates does not have an impact on any of the gene rearrangement scenarios depicting neutral correlation. Selection pressure analysis revealed that the mitogenomes from members of both the subfamilies Aleurodicinae and Aleyrodinae are characterized by intense purifying selection pressure.