Multiple independent origins of duplicated mitochondrial control regions indicate an apomorphy in the Thysanoptera (Insecta)

The mitochondrial genome (mitogenome) of thrips is characterized by the presence of control region (CR) duplication. However, the evolution pattern of duplicated CRs in thrips is still unclear. In this study, the multiple independent origins of duplicated CR indicated that the CR duplication was not an ancestral state for Thysanoptera. The macroevolutionary pattern suggested that the earliest CR duplication event occurred in the middle Cretaceous (94.85 Ma) coincided with rearrangement events forming the ancestors of Aeolothripidae, but much later than that forming the ancestors of the suborder Terebrantia. The mitogenome with duplicated CRs showed a higher rate of gene rearrangement. The sequence similarity of the CR copies and divergence time were negatively correlated, indicating age-related deterioration of mitochondrial function. No significant differences were found in the mitochondrial DNA, the P₁₂₃ and P_4FD between the single and multiple-CR charactered mitogenomes, which suggested that the duplicated CRs may not affect the replication process in thrip mitogenome. The mitogenomes with duplicated CRs (mean: 0.0088 subs/s/my) show a significantly increased evolutionary rate than that with a single one (mean: 0.0058 subs/s/my). However, it seems that this higher evolutionary rate did not have adaptive mechanisms in Terebrantia. We speculated that the duplicated CRs may cause a more intense production of energy by mitochondria, and an accelerated mutation and substitution rate is expected in such mitogenomes. Our study provided new insights into the presence of CR duplications and their evolution in the mitogenomes of thrips.