Characterization of the complete mitochondrial genomes of six horseflies (Diptera: Tabanidae)

The family Tabanidae (Insecta: Diptera) is one of the economically most important group of haematophagous insects, causing millions of livestock deaths per year. However, current knowledge on the mitochondrial genomes (mitogenomes) from this family is limited. Additional tabanid mitogenomes characterization is of utmost importance for their identification, epidemiologic and phylogenetic studies. We sequenced the mt genomes of six horseflies with an Illumina platform and their phylogenetic relationship was conducted with other infraorder Tabanomorpha members with available mt genome datasets. All six newly sequenced mitogenomes were typical 37-gene circular structures retaining the gene order of Tabanomorpha. The trnQ, trnM and trnA were highly conserved among the six mitogenomes (identity = 100%). The TΨC arm and variable loop regions were relatively more variable compared to the amino acid receptor arm, anticodon arm and DHU arm of the tRNAs. Among 13 protein-coding genes (PCGs) of tabanids mitogenomes, the highest nucleotide diversity was detected in atp8, cox1, cox3, nad6 and cytb (0.1 for each). In addition, atp8 genes exhibited the highest evolutionary rate (ω = 0.24) among 13 PCGs. The interspecies K2P genetic distances among some Tabanus spp. across the mitogenome was greater (0.08) than intergeneric genetic distance between T. amaenus and Atylotus miser (0.07). Phylogenetic analyses revealed non-monophyletic relationships among horseflies of the genus Tabanus. The present study showed mt gene order is highly conserved within Tabanus species. Our mito-phylogenomic analysis supports the paraphyly of the genus Tabanus. The new data provide novel genetic markers for studies of population genetics and systematics of horseflies.