Three complete mitochondrial genomes of Blaberidae (Insecta: Blattodea) (Gromphadorhina portentosa, Panchlora nivea, Blaptica dubia) and one complete mt genome of Blattidae (Insecta: Blattodea) (Shelfordella lateralis) were sequenced to further understand the characteristics of cockroach mitogenomes and reconstruct the phylogenetic relationship of Blattodea. The gene order and orientation of these four cockroach genomes were similar to known cockroach mt genomes, and contained 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes and one control region. The mt genomes of Blattodea exhibited a characteristics of a high A+T composition (70.7%-74.3%) and dominant usage of the TAA stop codon. The AT content of the whole mt genome, PCGs and total tRNAs in G. portentosa was the lowest in known cockroaches. The presence of a 71-bp intergenic spacer region between trnQ and trnM was a unique feature in B. dubia, but absent in other cockroaches, which can be explained by the duplication/random loss model. Based on the nucleotide and amino acid datasets of the 13 PCGs genes, neighbor-joining (NJ), maximum parsimony (MP), maximum likelihood (ML) and bayesian inference (BI) analyses were used to rebuild the phylogenetic relationship of cockroaches. All phylogenetic analyses consistently placed Isoptera as the sister cluster to Cryptocercidae of Blattodea. Ectobiidae and Blaberidae (Blaberoidea) formed a sister clade to Blattidae. Corydiidae is a sister clade of all the remaining cockroach species with a high value in NJ and MP analyses of nucleotide and amino acid datasets, and ML and BI analyses of the amino acid dataset.

The heat shock 90/70 organizing protein (Hop), also known as Sti-1 (stress-induced protein-1), is a co-chaperone that usually mediates the interaction of Hsp90 and Hsp70 and has been extensively characterized in mammals and plants. However, its role in insects remains unknown. In the present study, we isolated and characterized a Hop homologue gene from Frankliniella occidentalis (Fohop). The Fohop contains a 1659bp ORF encoding a protein of 552 amino acids with a caculated molecular mass of approximately 62.25kDa, which displays a reasonable degree of identity with the known Hops and shares several canonical motifs, including three tetratricopeptide repeated motif domains (TPR1, TPR2A and TPR2B) and two aspartic acid-proline (DP) repeat motifs (DP1 and DP2). As in other hops, Fohop contains introns, but the number and the position are quite variable. The mRNA expression patterns indicated that Fohop was constitutively expressed throughout the developmental stages, but was obviously upregulated by heat stress both in larvae and adults. Our studies imply that Hop, as in other Hsps, may play an important role in heat shock response of F. occidentalis.

The mayfly species Siphluriscus chinensis (Siphluriscidae) has valuable structures useful for phylogeny reconstruction, given its putative basal position within the Ephemeroptera. Here its nearly complete mitochondrial genome is sequenced. We built phylogenetic trees through multiple analytical strategies with some other insect mitogenomes. Structurally, the obtained mitochondrial genome of S. chinensis is 16,616 bp in length,(1) containing 37 genes and an extra trnK-like (trnK2 (AAA)) gene. The 12 PCGs start with typical ATN codons, except the nad1 gene which starts with an unnormalized TTG. Like other known mayfly mitogenomes, the strand bias has negative AT-skew and negative GC-skew. Phylogenetically, our topologies suggest that Odonata is the basally diverged clade in Pterygota; Ephemeroptera is the sister group of the Neoptera; and S. chinensis is indeed the most basal mayfly branch.

The Plecoptera (stoneflies) is a hemimetabolous order of insects, whose larvae are usually used as indicators for fresh water biomonitoring. Herein, we describe the complete mitochondrial (mt) genome of a stonefly species, namely Acroneuria hainana Wu belonging to the family Perlidae. This mt genome contains 13 PCGs, 22 tRNA-coding genes and 2 rRNA-coding genes that are conserved in most insect mt genomes, and it also has the identical gene order with the insect ancestral gene order. However, there are three special initiation codons of ND1, ND5 and COI in PCGs: TTG, GTG and CGA, coding for L, V and R, respectively. Additionally, the 899-bp control region, with 73.30% A+T content, has two long repeated sequences which are found at the 3'-end closing to the tRNA(Ile) gene. Both of them can be folded into a stem-loop structure, whose adjacent upstream and downstream sequences can be also folded into stem-loop structures. It is presumed that the four special structures in series could be associated with the D-loop replication. It might be able to adjust the replication speed of two replicate directions.

Complete mitochondrial genomes for two corn rootworm species, Diabrotica virgifera virgifera (16,747 bp) and Diabrotica barberi (16,632; Insecta: Coleoptera: Chrysomelidae), were assembled from Illumina HiSeq2000 read data. Annotation indicated that the order and orientation of 13 protein coding genes (PCGs), and 22 tRNA and 2 rRNA sequences were in typical of insect mitochondrial genomes. Non-standard nad4 and cox3 stop codons were composed of single T nucleotides and likely completed by adenylation, and atypical TTT start codons was predicted for both D. v. virgifera and D. barberinad1 genes. The D. v. virgifera and D. barberi haplotypes showed 819 variable nucleotide positions within PCG regions (7.36% divergence), which suggest that speciation may have occurred ~3.68 million years ago assuming a linear rate of short-term substitution. Phylogenetic analyses of Coleopteran MtD genome show clustering based on family level, and may have the capacity to resolve the evolutionary history within this Order of insects.

In the present study, an in silico analysis was performed to identify transposable element (TE) fragments inserted in Cyps with functions associated with resistance to insecticides and developmental regulation as well as in neighboring genes in two sibling species, Drosophila melanogaster and Drosophila simulans. The Cyps associated with insecticide resistance and their neighboring non-Cyp genes have accumulated a greater number of TE fragments than the other Cyps or a random sample of genes, predominantly in the 5'-flanking regions. Most of the insertions were due to DNA transposons, with DNAREP1 fragments being the most common. These fragments carry putative binding sites for transcription factors, which reinforces the hypothesis that DNAREP1 may influence gene regulation and play a role in the adaptation of the Drosophila species.

The blister beetle is an important resource insect due to its defensive substance cantharidin, which was widely used in pharmacology and plant protection. We determined the complete mitochondrial genome of Epicauta chinensis Laporte (Coleoptera: Tenebrionoidae: Meloidae). The circular genome is 15,717 bp long, encoding 13 protein-coding genes (PCGs), two ribosomal RNAs and 22 tRNAs and containing a A+T-rich region with gene arrangement identical to other Coleopteran species. Twelve PCGs start with typical ATN codon, while ATP8 gene initiate with GTT for first report in Insecta. All PCGs terminate with conventional stop codon TAA or TAG. All tRNAs in E. chinensis are predicted to fold into typical cloverleaf secondary structure, except tRNA-Ser(AGN), in which the dihydrouracil arm (DHU arm) could not form stable stem-loop structure. The secondary structure of lrRNA and srRNA comprises 48 helices and 32 helices respectively. The 1101 bp A+T-rich region contains a 15 bp poly-T stretch and microsatellite-like repeats rather than large tandem repetitive sequences. Phylogenetic analysis, based on 13 PCGs of 45 Coleopteran species, show that E. chinensis grouped with Tenebrionidae species. It also support the topology of (((Chrysomelidae+Curculionoidea)+(Cucujoidea+Cleroidea))+Tenebrionoidea) within Cucujiformia.

Next-generation sequencing has dramatically fostered insect mitogenomic research in recent years. However, studies on the insect mitochondrial genome (mitogenome) assembly mainly rely on the sequencing data from total DNA, which is not cost-effective as a huge data from nuclear DNA are wasted. Besides, many mitogenomic studies require genomic information from individual organisms, whereas the DNA yield from small individual insects is too low to meet the sequencing requirements. Here, we describe a strategy for a high enrichment of insect mitochondrial DNA (mtDNA) using rolling circle amplification (RCA) technique. This strategy consists of standard DNA extraction, RCA enrichment, next-generation sequencing and mitogenome assembly. We have evaluated the performance of this strategy on nine insect species representing eight families of insecta, three other invertebrates, and even two vertebrate specimens. Results show that our strategy is especially suitable for insects, which allows almost all tested insect mtDNA contents to reach 80% and above. A further examination of enrichment efficiency of our strategy among different taxa shows that it is also applicable to other invertebrates and even some vertebrates such as Rhacophorus and ptyas species, although its enrichment efficiency in these groups is lower than that of insects. After treatment with our strategy, small flux sequencing data can realize the assembly of mitogenome with deep coverage, providing a solid base for subsequent mitogenome-based studies.