Stoneflies comprise an ancient group of insects, but the phylogenetic position of Plecoptera and phylogenetic relations within Plecoptera have long been controversial, and more molecular data is required to reconstruct precise phylogeny. Herein, we present the complete mitogenome of a stonefly, Suwallia teleckojensis, which is 16146 bp in length and consists of 13 protein-coding genes (PCGs), 2 ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs) and a control region (CR). Most PCGs initiate with the standard start codon ATN. However, ND5 and ND1 started with GTG and TTG. Typical termination codons TAA and TAG were found in eleven PCGs, and the remaining two PCGs (COII and ND5) have incomplete termination codons. All transfer RNA genes (tRNAs) have the classic cloverleaf secondary structures, with the exception of tRNA^Ser(AGN), which lacks the dihydrouridine (DHU) arm. Secondary structures of the two ribosomal RNAs were shown referring to previous models. A large tandem repeat region, two potential stem-loop (SL) structures, Poly N structure (2 poly-A, 1 poly-T and 1 poly-C), and four conserved sequence blocks (CSBs) were detected in the control region. Finally, both maximum likelihood (ML) and Bayesian inference (BI) analyses suggested that the Capniidae was monophyletic, and the other five stonefly families form a monophyletic group. In this study, S. teleckojensis was closely related to Sweltsa longistyla, and Chloroperlidae and Perlidae were herein supported to be a sister group.

The family-level relationships within Plecoptera have been a focused area of research for a long time. Its higher classification remains unstable, and the phylogenetic relationships within Plecoptera should be re-examined. Here, we sequenced and analyzed two complete mitochondrial genomes (mitogenomes) of Paraleuctra cercia and Perlomyia isobeae of the family Leuctridae. We reconstructed the phylogeny of Plecoptera based on 13 protein-coding genes (PCGs) from published stoneflies. Our results showed that the Bayesian inference and maximum-likelihood tree had similar topological structures except for the positions of two families, Peltoperlidae and Scopuridae. The Plecoptera is divided into two clades, the suborder Antarctoperlaria and the suborder Arctoperlaria. The two suborders subsequently formed two groups, Eusthenioidea and Gripopterygoidea, and Euholognatha and Systellognatha, which is consistent with the results of morphological studies. In addition, the Leuctridae is the earliest branch within the superfamily Nemouroidea. But the monophyly of Perloidea and Pteronarcyoidea are still not well supported.

The mitochondrial genome of Etrocorema hochii (Wu, 1938), the first representative of Etrocorema, was sequenced and annotated in this study. The complete circular mitochondrial genome (mitogenome) was 15,854 bp in length using Illumina sequencing and harboured 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNA), two ribosomal RNA genes (rRNA), and one control region in the typical stonefly gene order. The nucleotide composition of this genome is A (36.0%), T (31.4%), C (20.9%), and G (11.7%). All PCGs use the typical start codon ATN, except ND1, ND2, and COI, which start with TTG, GTG, and ACC, respectively. Only the ND5 genes terminated with single T--, and the rest use the stop codon TAA/TAG. We used this new nucleotide information under Bayesian and maximum-likelihood analyses to explore the relationships of Pelidae. Phylogenetic analyses showed that E. hochii is closely related to the clade (Togoperla + Kamimuria), which is widely accepted view.

The genus-level relationships within the subfamily Amphinemurinae have been controversial, although attempts have been made based on morphology and limited molecular data. With the establishment of two new genera, the phylogenetic relationships within Amphinemurinae should be re-examined. In this study, the complete mitochondrial genome (mitogenome) of Malenka flexura of the genus Malenka was firstly sequenced and analyzed. The phylogeny of Amphinemurinae was also reconstructed using 13 proteincoding genes (PCGs) from previously published stoneflies. This mitogenome was 15,744 bp long and encoded the typical 37 genes, as well as a putative control region. The gene arrangement of M. flexura mitogenome is identical with the putative ancestral mitogenome in Drosophila yakuba. Most PCGs used standard ATN as start codons and TAA/TAG as termination codons. All tRNA genes exhibited the typical cloverleaf secondary structure, except for tRNA^Ser(AGN), whose dihydrouridine (DHU) arm was lacking. Some structural elements in the control region were founded, such as tandem repeat regions, stemloop structures, polyN stretch and microsatellite structure, etc. Phylogenetic analyses of sequenced Amphinemurinae mitogenomes unsupported the sister relationship of Amphinemura and Malenka. Finally, the phylogenetic analyses inferred a relationship within Amphinemurinae: Amphinemura + (Malenka + (Protonemura + (Indonemoura + (Sphaeronemoura + Mesonemoura)))).

Currently, the classification system of 2 subfamilies within Nemouridae has been widely accepted. However, monophyly of 2 subfamilies has not been well supported by molecular evidence. To date, only mitogenomes from genus Nemoura of the subfamily Nemourinae were used in previous phylogenetic studies and produced conflicting results with morphological studies. Herein, we analyzed mitogenomes of 3 Nemourinae species to reveal their mitogenomic characteristics and to examine genus-level classification among Nemouridae. In this study, the genome organization of 3 mitogenomes is highly conserved in gene order, nucleotide composition, codon usage, and amino acid composition. In 3 Nemourinae species, there is a high variation in nucleotide diversity among the 13 protein-coding genes (PCGs). The K_a/K_s values for all PCGs were far lower than 1, indicating that these genes were evolving under purifying selection. The phylogenetic analyses highly support Nemurella as the sister group to Ostrocerca. Meanwhile, Nemoura is recovered as the sister group of Malenka; they are grouped with other Amphinemurinae and emerged from a paraphyletic Nemourinae. More molecular data from different taxonomic groups are needed to understand stoneflies phylogeny and evolution.