Phylogenetic relationships of Limacodidae and insights into the higher phylogeny of Lepidoptera

To determine the systematic status of family Limacodidae within Lepidoptera, the complete mitochondrial genome (mitogenome) of Thosea sinensis (Lepidoptera: Zygaenoidea: Limacodidae) was sequenced. The genome is 15,544 base pairs (bp), including 13 protein-coding genes (PCGs), two ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs), and an AT-rich region. These characteristics are similar to of other lepidopterans. The gene order of T. sinensis is identical to that of Ditrysia lepidopterans. The nucleotide composition of the T. sinensis mitochondrial genome is highly biased toward A + T nucleotides (81.1%) and exhibits negative AT and GC skew. All the other 13 PCGs except cox1 are initiated by ATN codons. All tRNA genes are folded into the typical cloverleaf secondary structure, except for trnS1, which lacked the dihydrouridine (DHU) stem. There are 20 intergenic spacer regions ranging from 1 to 56 bp in length, and two gene overlap regions throughout the entire genome. The AT-rich region includes the ATAGA motif, followed by a 19-bp poly T stretch, a microsatellite-like (AT)10, and a poly-A element. Analysis of phylogenetic relationships indicated that T. sinensis belongs to the Limacodidae, and the monophyly of each lepidopteran family was well supported.