The complete mitochondrial genome and its remarkable secondary structure for a stonefly Acroneuria hainana Wu (Insecta: Plecoptera, Perlidae)

The mitochondrial genome of Neoperla bimaculata (Li et al. 2021) (plecoptera: perlidae) from Tibet of southwest China and its phylogenetic analysis

The complete mitochondrial genome (mitogenome) of Neoperla bimaculata was sequenced and annotated in this study. We found that the mitogenome of N. bimaculata is 15,774 bp in length with an A + T content of 64.3%. It exhibits the classic structure of a mitogenome. Most protein-coding genes (PCGs) of the mitogenome initiate with the standard start codon ATN. Ten PCGs use the standard stop codon TAA/TAG, while the COI, COII, and ND5 genes terminate with a single T nucleotide. Phylogenetic analyses suggested that N. bimaculata, along with two unpublished Neoperla species, formed a cluster within the phylogenetic tree. Our results indicated that the genus Neoperla and Neoperlops were sister groups. Meanwhile, the monophyly of Perlinae and Acroneuriinae was supported in the mitochondrial phylogeny.

Mitochondrial Genomes of the Genus Claassenia (Plecoptera: Perlidae) and Phylogenetic Assignment to Subfamily Perlinae

Mitochondrial genomes of three stoneflies, e.g., Claassenia magna Wu, 1948, Claassenia sp. 2 and Claassenia xucheni Chen, 2019 were sequenced in this study with 15,774, 15,777 and 15,746 bp in length, respectively. Each mitogenome contained 37 genes including 22 tRNAs, two ribosomal RNAs, 13 protein-coding genes (PCGs), and a noncoding control region (CR). In general, standard ATN start and TAN termination codons were evident in the PCGs. Although the dihydrouridine arm was absent in trnSer, the remaining 21 tRNAs displayed the characteristic cloverleaf secondary structure. Stem-loop structures were identified in the CRs of all three mitogenomes, but tandem repeats were only apparent in Claassenia xucheni. The mitogenomes of three Claassenia species were analyzed and compared with mitogenomes in 21 other stoneflies from the Perlidae and three Euholognatha species (Rhopalopsole bulbifera, Capnia zijinshana and Amphinemura longispina) as outgroups. Phylogenetic analyses using maximum likelihood and Bayesian inference. Phylogenetic analysis supported that Claassenia was recovered as the sister group of other Perlinae and Claassenia+Perlinae emerged from the paraphyletic Acroneuriinae. The final results supported that Claassenia was classified into subfamily Perlinae and proposed Claassenia represent a transitional group of the subfamilies Acroneuriinae and Perlinae. This study provided new molecular evidence for exploring the debatable taxonomic position of the genus Claassenia in Perlidae.

Three Complete Mitochondrial Genomes of Orestes guangxiensis, Peruphasma schultei, and Phryganistria guangxiensis (Insecta: Phasmatodea) and Their Phylogeny

Insects of the order Phasmatodea are mainly distributed in the tropics and subtropics and are best known for their remarkable camouflage as plants. In this study, we sequenced three complete mitochondrial genomes from three different families: Orestes guangxiensis, Peruphasma schultei, and Phryganistria guangxiensis. The lengths of the three mitochondrial genomes were 15,896 bp, 16,869 bp, and 17,005 bp, respectively, and the gene composition and structure of the three stick insects were identical to those of the most recent common ancestor of insects. The phylogenetic relationships among stick insects have been chaotic for a long time. In order to discuss the intra- and inter-ordinal relationship of Phasmatodea, we used the 13 protein-coding genes (PCGs) of 85 species for maximum likelihood (ML) and Bayesian inference (BI) analyses. Results showed that the internal topological structure of Phasmatodea had a few differences in both ML and BI trees and long-branch attraction (LBA) appeared between Embioptera and Zoraptera, which led to a non-monophyletic Phasmatodea. Consequently, after removal of the Embioptera and Zoraptera species, we re-performed ML and BI analyses with the remaining 81 species, which showed identical topology except for the position of Tectarchus ovobessus (Phasmatodea). We recovered the monophyly of Phasmatodea and the sister-group relationship between Phasmatodea and Mantophasmatodea. Our analyses also recovered the monophyly of Heteropterygidae and the paraphyly of Diapheromeridae, Phasmatidae, Lonchodidae, Lonchodinae, and Clitumninae. In this study, Peruphasma schultei (Pseudophasmatidae), Phraortes sp. YW-2014 (Lonchodidae), and species of Diapheromeridae clustered into the clade of Phasmatidae. Within Heteropterygidae, O. guangxiensis was the sister clade to O. mouhotii belonging to Dataminae, and the relationship of (Heteropteryginae + (Dataminae + Obriminae)) was recovered.

Increasing 28 mitogenomes of Ephemeroptera, Odonata and Plecoptera support the Chiastomyaria hypothesis with three different outgroup combinations

BACKGROUND

The phylogenetic relationships of Odonata (dragonflies and damselflies) and Ephemeroptera (mayflies) remain unresolved. Different researchers have supported one of three hypotheses (Palaeoptera, Chiastomyaria or Metapterygota) based on data from different morphological characters and molecular markers, sometimes even re-assessing the same transcriptomes or mitochondrial genomes. The appropriate choice of outgroups and more taxon sampling is thought to eliminate artificial phylogenetic relationships and obtain an accurate phylogeny. Hence, in the current study, we sequenced 28 mt genomes from Ephemeroptera, Odonata and Plecoptera to further investigate phylogenetic relationships, the probability of each of the three hypotheses, and to examine mt gene arrangements in these species. We selected three different combinations of outgroups to analyze how outgroup choice affected the phylogenetic relationships of Odonata and Ephemeroptera.

METHODS

Mitochondrial genomes from 28 species of mayflies, dragonflies, damselflies and stoneflies were sequenced. We used Bayesian inference (BI) and Maximum likelihood (ML) analyses for each dataset to reconstruct an accurate phylogeny of these winged insect orders. The effect of outgroup choice was assessed by separate analyses using three outgroups combinations: (a) four bristletails and three silverfish as outgroups, (b) five bristletails and three silverfish as outgroups, or (c) five diplurans as outgroups.

RESULTS

Among these sequenced mitogenomes we found the gene arrangement IMQM in Heptageniidae (Ephemeroptera), and an inverted and translocated tRNA-Ile between the 12S RNA gene and the control region in Ephemerellidae (Ephemeroptera). The IMQM gene arrangement in Heptageniidae (Ephemeroptera) can be explained via the tandem-duplication and random loss model, and the transposition and inversion of tRNA-Ile genes in Ephemerellidae can be explained through the recombination and tandem duplication-random loss (TDRL) model. Our phylogenetic analysis strongly supported the Chiastomyaria hypothesis in three different outgroup combinations in BI analyses. The results also show that suitable outgroups are very important to determining phylogenetic relationships in the rapid evolution of insects especially among Ephemeroptera and Odonata. The mt genome is a suitable marker to investigate the phylogeny of inter-order and inter-family relationships of insects but outgroup choice is very important for deriving these relationships among winged insects. Hence, we must carefully choose the correct outgroup in order to discuss the relationships of Ephemeroptera and Odonata.

The complete mitochondrial genome of Flavoperla biocellata Chu, 1929 (Plecoptera: Perlidae) and the phylogenetic analyses of Plecoptera

Of the roughly 400 species of Perlidae in the world, most species are widely distributed in the northern hemisphere, but a few can be found in South Africa and South America. There are only five species in the genus Flavoperla of the family Perlidae in China. To gain a better understanding of the architecture and evolution of mitochondrial genome in Flavoperla, the entire mitochondrial genome (mitogenome) of a Chinese Flavoperla biocellata Chu, 1929 from family Perlidae (Insecta: Plecoptera) was sequenced. The 15,805-bp long mitochondrial genome of F. biocellata contained 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a putative control region (CR). The gene arrangement of F. biocellata was identical with that of other stoneflies and with the fly Drosophila yakuba. Most PCGs of F. biocellata used the standard ATN start codons and complete TAN termination codons. Twenty-one of the 22 tRNA genes exhibited cloverleaf secondary structures, but the dihydrouridine (DHU) arm of trnSer (AGN) was completely reduced. Phylogenetic analyses with both Bayesian inference (BI) and maximum likelihood methods (ML) generated similar topology, both supporting the monophyly of all stonefly families and the infraorder Systellognatha. The phylogenetic analysis based on mitochondrial genomic data from 30 stonefly species recovered a well-supported tree resolving higher-level relationships within Plecoptera. The northern hemisphere suborder Arctoperlaria divided into two groups, Euholognatha and Systellognatha. The southern hemisphere suborder Antarctoperlaria formed two clades: Eustheniidae+Diamphipnoidae and Austroperlidae+ Gripopterygidae; consistent with relationships proposed based on morphology. The final relationships within Plecoptera were recovered as (((Perlidae+(Perlodidae+Chloroperlidae))+(Pteronarcyidae+(Peltoperlidae+Styloperlidae))) +(Taeniopterygidae+(Capniidae+(Nemouridae+Notonemouridae))))+ (Gripopterygoidae+Eusthenioidae).

The mitochondrial genome analysis of Flavoperla hatakeyamae (Plecoptera: Perlidae)

In this study, we analyzed the entire mitochondrial genome of Flavoperla hatakeyamae, which has 15,730 bp in length. This complete mitochondrial genome possessed 13 protein coding genes (PCGs), 2 ribosomal RNA genes, 22 transfer RNA genes, and 1 control region. The nucleotide composition of the complete sequence of mitochondrial genome was accounting for 35.0%, T for 31.3%, C for 21.2%, G for 12.5%, and A + T for 66.3%. Among the 13 protein-coding genes, the starting codon of 12 protein-coding genes was ATN except the starting codon of COII was ACC. In addition, 11 of the protein-coding genes used conservative termination codon TAA/TAG, except for COII and ND5 which terminated by a single T. Based on Bayesian (BI) and maximum-likelihood (ML) methods, we found that the genus Flavoperla and Niponiella were gathered into one branch and the species of Perlinae were clustered in one clade.

Characterization of the complete mitochondrial genome of a stonefly, Perlesta teaysia (Plecoptera: Perlidae)

In this study, we sequenced the mitochondrial genome of Perlesta teaysia. The complete mitochondrial genome was 16,023 bp in length, including 37 typical genes and a control region. The overall nucleotide composition was biased toward the A/T nucleotides. All the protein-coding genes (PCGs) started with typical ATN codon as the start codon except for ND1 and ND5, which started with TTG and GTG, respectively. In addition, 11 of 13 PCGs possessed the typical stop codon TAA/TAG, whereas ND4 and ND5 terminated with an incomplete stop codon. Phylogenetic analysis suggested that P. teaysia was a sister group to Calineuria stigmatica.

The complete mitochondrial genome analysis of the stonefly, Flavoperla sp. (Plecoptera: Perlidae)

Flavoperla is a small genus in Acroneuriinae, which was established by Chu but was treated as a synonym for a long time until Uchida and Harper recognized Flavoperla as a valid genus. This genus includes 15 species in the world now, including six species from China. Here, we sequenced and annotated the mitochondrial genome of Flavoperla sp. under the next-generation sequence technology, the first representative of Genus Flavoperla. The complete mitochondrial genome of Flavoperla sp. is 15,796 bp in length with an A + T content of 68.3% showed a positive AT-skew (0.037) and a negative GC-skew (-0.262). The gene order and organization of the mitochondrial genome is consistent with other stoneflies. The control region had the highest A + T content. Eleven PCGs started with the typical codon (ATN), the remaining PCGs started with TTG (ND1and ND5). All PCGs terminated with TAA/TAG, except COII and ND5 used single T. The phylogenetic tree by maximum likelihood (ML) and Bayesian (BI) methods supported that Flavoperla sp. was the sister group to Niponiella limbatella.

The characterization of the mitochondrial genome of Calineuria stigmatica (Plecoptera: Perlidae)

The nearly complete sequence of the mitochondrial DNA of Calineuria stigmatica has been completed and annotated in this study. The circular genome is 15,070 bp in length with an A + T content of 61.8% and contains 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, and two ribosomal RNA (rRNA) genes. The control region can only be assembled partially. All PCGs use normal start codon ATN, while COI, ND1, and ND5 use CCG, TTG, and GTG as start codon, respectively. Meanwhile, 10 PCGs use the typical termination codons TAN, except COII, ND4, ND5, which stopped with the incomplete terminaton signal T--. Based on 13 PCGs and two rRNAs using the Bayesian (BI) method supported that C. stigmatica was closely grouped with four other Acroneuriinae species. Our results provide basic data for further study of phylogeny in Plecoptera.

The complete mitochondrial genome of a stonefly species, Etrocorema hochii (Plecoptera: Pelidae)

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.

Comparative mitogenomic analysis of species in the subfamily Amphinemurinae (Plecoptera: Nemouridae) reveal conserved mitochondrial genome organization

The subfamily Amphinemurinae has five genera in China, with each genus of similar morphology. To gain a better understanding of architecture and evolution of mitogenome in Amphinemurinae, mitogenomes of eight species representing four genera (Amphinemura, Indonemoura, Protonemura and Sphaeronemoura) in the subfamily Amphinemurinae were sequenced, and a comparative mitogenomic analysis of five genera (including a published stonefly genus, Mesonemoura) was carried out. By comparative analysis, we found highly conserved genome organization of ten Amphinemurinae species including genome contents, gene order, nucleotide composition, codon usage, amino acid composition, as well as genome asymmetry. GC content was the most significant factor in determining codon bias among organisms. The Ka/Ks values for all PCGs were far lower than 1, indicating that these genes were evolving under purifying selection. We also found some important conserved stem and loop in the cloverleaf structure of tRNAs, and found conserved helices and loops in each domain of the secondary structure of rRNAs. The presence of structural elements in the control region is also discussed. The phylogenetic analyses indicated that within Amphinemurinae, Sphaeronemoura was assigned the sister group of Mesonemoura. Our analyses inferred a relationship within Euholognatha: ((Nemouridae + Notonemouridae) + (Taeniopterygidae + Capniidae) + Scopuridae) + Leuctridae.

Mitochondrial genomes of the stoneflies Mesonemourametafiligera and Mesonemouratritaenia (Plecoptera, Nemouridae), with a phylogenetic analysis of Nemouroidea

In this study, two new mitochondrial genomes (mitogenomes) of Mesonemourametafiligera and Mesonemouratritaenia from the family Nemouridae (Insecta: Plecoptera) were sequenced. The Mesonemourametafiligera mitogenome was a 15,739 bp circular DNA molecule, which was smaller than that of M.tritaenia (15,778 bp) due to differences in the size of the A+T-rich region. Results show that gene content, gene arrangement, base composition, and codon usage were highly conserved in two species. Ka/Ks ratios analyses of protein-coding genes revealed that the highest and lowest rates were found in ND6 and COI and that all these genes were evolving under purifying selection. All tRNA genes in nemourid mitogenomes had a typical cloverleaf secondary structure, except for tRNA^Ser(AGN) which appeared to lack the dihydrouridine arm. The multiple alignments of nemourid lrRNA and srRNA genes showed that sequences of three species were highly conserved. All the A+T-rich region included tandem repeats regions and stem-loop structures. The phylogenetic analyses using Bayesian inference (BI) and maximum likelihood methods (ML) generated identical results. Amphinemurinae and Nemourinae were sister-groups and the family Nemouridae was placed as sister to Capniidae and Taeniopterygidae.

The phylogeny and evolutionary timescale of stoneflies (Insecta: Plecoptera) inferred from mitochondrial genomes

Phylogenetic analysis based on mitochondrial genomic data from 25 stonefly species recovered a well-supported tree resolving higher-level relationships within Plecoptera (stoneflies). The monophyly of both currently recognized suborders was strongly supported, concordant with previous molecular analyses of Plecoptera. The southern hemisphere suborder Antarctoperlaria formed two clades: Eustheniidae + Diamphipnoidae and Austroperlidae + Gripopterygidae; consistent with relationships proposed based on morphology. The largely northern hemisphere suborder Arctoperlaria also divided into two groups, Euholognatha and Systellognatha, each composed of the five families traditionally assigned to each infraorder (the placement Scopuridae by mt genome data remains untested at this time). Within Euholognatha, strong support for the clade Nemouridae + Notonemouridae confirmed the northern origin of the currently southern hemisphere restricted Notonemouridae. Other family level relationships within the Arctoperlaria differ from those recovered by previous morphology and molecular based analyses. A fossil-calibrated divergence estimation suggests the formation of two suborders dates back to the Jurassic (181 Ma), with subsequent diversification of most stonefly families during the Cretaceous. This result confirms the hypothesis that initial divergence between the suborders was driven by the breakup of the supercontinent Pangaea into Laurasia and Gondwanaland (commencing 200 Ma and complete by 150 Ma).

The complete mitochondrial genome of the Caribbean spiny lobster Panulirus argus

Panulirus argus is a keystone species and target of the most lucrative fishery in the Caribbean region. This study reports, for the first time, the complete mitochondrial genome of Panulirus argus (average coverage depth nucleotide-1 = 70×) assembled from short Illumina 150 bp PE reads. The AT-rich mitochondrial genome of Panulirus argus was 15 739 bp in length and comprised 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. A single 801 bp long intergenic space was assumed to be the D-loop. Most of the PCGs were encoded on the H-strand. The gene order observed in the mitochondrial genome of Panulirus argus corresponds to the presumed Pancrustacean ground pattern. KA/KS ratios calculated for all mitochondrial PCGs showed values < 1, indicating that all these PCGs are evolving under purifying selection. A maximum likelihood phylogenetic analysis (concatenated PCGs [n = 13], 154 arthropods) supported the monophyly of the Achelata and other infraorders within the Decapoda. Mitochondrial PCGs have enough phylogenetic informativeness to explore high-level genealogical relationships in the Pancrustacea. The complete mitochondrial genome of the Caribbean spiny lobster Panulirus argus will contribute to the better understanding of meta-population connectivity in this keystone overexploited species.

The mitochondrial genome of Caenis sp. (Ephemeroptera: Caenidae) and the phylogeny of Ephemeroptera in Pterygota

The phylogenetic relationship between Ephemeroptera (mayflies) and Odonata (dragonflies and damselflies) remains hotly debated in the insect evolution community. We sequenced the complete mitochondrial genome of Caenis sp. (Ephemeroptera: Caenidae) to discuss the phylogenetic relationship of Palaeoptera. The mitochondrial genome of Caenis sp. is a circular molecule of 15,254 bp in length containing 37 genes (13 protein-coding genes, 22 tRNAs, and 2 rRNAs), which showed the typical insect mitochondrial gene arrangement. In BI and ML phylogenetic trees using 71 species of 12 orders, our results support the Ephemeroptera as the basal group of winged insects.

Complete Mitochondrial Genome of Suwallia teleckojensis (Plecoptera: Chloroperlidae) and Implications for the Higher Phylogeny of Stoneflies

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 tRNASer(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.

First Mitochondrial Genome from Nemouridae (Plecoptera) Reveals Novel Features of the Elongated Control Region and Phylogenetic Implications

The complete mitochondrial genome (mitogenome) of Nemoura nankinensis (Plecoptera: Nemouridae) was sequenced as the first reported mitogenome from the family Nemouridae. The N. nankinensis mitogenome was the longest (16,602 bp) among reported plecopteran mitogenomes, and it contains 37 genes including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes and two ribosomal RNA (rRNA) genes. Most PCGs used standard ATN as start codons, and TAN as termination codons. All tRNA genes of N. nankinensis could fold into the cloverleaf secondary structures except for trnSer (AGN), whose dihydrouridine (DHU) arm was reduced to a small loop. There was also a large non-coding region (control region, CR) in the N. nankinensis mitogenome. The 1751 bp CR was the longest and had the highest A+T content (81.8%) among stoneflies. A large tandem repeat region, five potential stem-loop (SL) structures, four tRNA-like structures and four conserved sequence blocks (CSBs) were detected in the elongated CR. The presence of these tRNA-like structures in the CR has never been reported in other plecopteran mitogenomes. These novel features of the elongated CR in N. nankinensis may have functions associated with the process of replication and transcription. Finally, phylogenetic reconstruction suggested that Nemouridae was the sister-group of Capniidae.