Mitochondrial genome of the stonefly Kamimuria wangi (Plecoptera: Perlidae) and phylogenetic position of plecoptera based on mitogenomes

Mitochondrial genomes of Nemourinae species (Plecoptera: Nemouridae) and the phylogenetic implications

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 Ka/Ks 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.

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.

Molecular phylogeny inferred from the mitochondrial genomes of Plecoptera with Oyamia nigribasis (Plecoptera: Perlidae)

In this study, the mitochondrial genome of the stonefly, Oyamia nigribasis Banks, 1920 (Plecoptera: Perlidae), was sequenced and compared with the mtDNA genomes of 38 other stoneflies and two Ephemerae. The O. nigribasis mitogenome is a circular 15,923 bp molecule that encodes a large, noncoding control region (CR) and 37 typical mtDNA genes; these include 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and two ribosomal RNA genes (rRNAs), respectively. Most of the PCGs initiated with ATN and terminated with TAN. The dihydrouridine (DHU) arm of tRNASer (AGN) was missing, whereas the other 21 tRNAs all exhibited the typical cloverleaf secondary structure. Stem-loop (SL) structures and tandem repeats were identified in the CR. Phylogenetic analyses using Bayesian inference and maximum likelihood were undertaken to determine relationships between stoneflies. Results indicated that the Antarctoperlaria, which contains Gripopterygidae, was absolutely separated from Arctoperlaria; this finding agrees with morphology. Finally, the overall relationships could be summarized as follows ((((Notonemouridae + Nemouridae) + Leuctridae) + (Scopuridae + (Capniidae + Taeniopterygidae))) + (((Perlodidae + Chloroperlidae) + Perlidae) + (Pteronarcyidae + (Peltoperlidae + Styloperlidae))) + ((Diamphipnoidae + Eustheniidae) + Gripopterygidae)).

The mitochondrial genome of the stonefly Togoperla limbata Pictet (Plecoptera: Perlidae)

We sequenced the complete mitochondrial genomes of Togoperla limbata which is the unique species distributed in Japan to provide supplementary data for future study in phylogenetic studies of Perlidae. This complete mitogenome of T. limbata is 15,915 bp long, including 37 typical genes and a control region as other stoneflies. Through the data analysis, the A + T content of the whole mitogenome, PCGs, tRNAs, rRNAs, and control region was 63.7%, 61.3%, 68.7%, 66.8%, and 74.8%. Most PCGs used typical start or stop codon except COI, ND5, and ND1 used the exceptional start codon and the COII and ND5 used the single T as the stop codon. Phylogenetic tree was constructed by using Bayesian inference (BI) and maximum-likelihood (ML) analysis based on the sequences of the 13 PCGs and two rRNAs and the result showed that subfamily Perlinae was a monophyletic group and the clade Togoperla and Paragnetina had a closer relationship.

The mitochondrial genome analysis of Paragnetina indentata (Plecoptera: Perlidae)

In this study, the entire mitochondrial genome of Paragnetina indentata was sequenced. The circular mitochondrial genome was the first mitochondrial genome representing the genus Paragnetina, which consists of 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, 22 transfer RNA genes, and 1 control region with the length of 15,885 bp. In the complete sequence, the A + T content for 64.1% (A for 33.6%, T for 30.5%, C for 23.3%, and G for 12.7%). Among the 13 PCGs, 11 start codons is ATN and the start codon of ND1, ND2, and ND5 genes is TTG, GTG, and GTG, respectively. In addition, 11 of the PCGs used conservative termination codon TAA or TAG, except for COII and ND5 which terminated by the single T. By using the Bayesian inference (BI) and maximum-likelihood (ML) methods, the phylogenetic relationship showed that P. indentata was closely related to Togoperla sp. and the species of Perlinae were clustered in a clade.

The mitochondrial genome from the stonefly: Claassenia sp. (Plecoptera: Perlidae)

Claasseniini is the apparent sister group to all the other groups of Perlidae, and the small tribe Claasseniini has only one genus includes 12 species in the world. Only five mitochondrial genomes in the Perlinae have previously been reported. In this study, we sequenced the mitochondrial genome of Claassenia sp. Wu in genus Claassenia for the first time. The total length of the mitochondrial genome is 15,774 bp, and the DNA base composition is 34.0% A, 27.4% T, 25.6% C and 13.0% G. The whole genome contains 37 genes and one control region, which is similar to other published stoneflies. All PCG started with ATN except CCG for COI gene, TTG for ND1 gene, GTG for ND4 and ND5 genes. The termination codon of 11 PCGs used the normal complete codon TAA or TAG, and the COII and ND5 genes were terminated by a single T. All tRNAs except tRNASer(AGN) fold into a classic clover structure. The phylogenetic tree by Maximum likelihood (ML) and Bayesian (BI) methods supported that Claassenia sp. was the basal taxa of the Perlidae clade. This study provides new data for further mitochondrial genome study on the Perlidae.

Characterization of the complete mitochondrial genome of a stonefly species, Kamimuria klapaleki (Plecoptera: Pelidae)

We have sequenced and analyzed complete mitochondrial genome (mitogenomes) of the Kamimuria klapaleki, the third mitochondrial genome in the genus Kamimuria, which belongs to the family Pelidae in this paper. The mitogenome of K. klapaleki is circular with the length of 16,077 bp, which possessed 37 genes and a control region like other stonefly. The A + T content of the whole mitogenome was 67.1%. All PCGs were 11,217 bp in length, accounting for 65.3% of the content of A + T. The A + T content of lrRNA and srRNA were 72.3% and 66.4%. The highest A + T content was in the control region (76.2%). There are 15 gene overlaps and 9 gene intergenic spacers in this mitochondrial genome. In addition, we performed phylogenetic analysis by using the Bayesian (BI) and Maximum Likelihood (ML) methods based on the concatenated data set of PCGs from 11 species in Pelidae and two species in Styloperlidae (outgroups). The clade K. chungnanshana + K. wangi was a sister group to K. klapaleki, which is consistent with the traditional morphological classification.

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.

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).

Comparative analysis of the Liriomyza chinensis mitochondrial genome with other Agromyzids reveals conserved genome features

Liriomyza chinensis is a serious pest of onions in many countries, especially in East Asia. We sequenced the complete mitochondrial genome of this species and compared it with five other Agromyzidae species. The L. chinensis mitogenome is a double-stranded 16,175 bp circular molecule with an A + T content of 78.3%. It contains 37 genes and a control region as do the sequenced Liriomyza species. The mitogenomes of L. chinensis and other Agromyzidae species showed a clear bias in nucleotide composition with a positive AT-skew. Most PCGs used standard ATN as start codons, and TAN as termination codons. The tRNAs exhibited the typical clover-leaf structure, except for tRNASer^(AGN) and the two rRNA genes are conserved with those of other Agromyzids. The L. chinensis mitogenome control region included several conserved regions, including a poly-T, two (TA)n and one poly-A stretch, which are considered important replication and transcription. The 13 PCGs were used to study the phylogeny of L. chinensis and five related Agromyzids. Analysis by maximum likelihood, Bayesian inference and genetic distance suggest congruent phylogenetic relationships in Liriomyza spp. and provide a useful supplement to taxonomic classification by morphology.

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.

The Complete Mitogenome of the Wood-Feeding Cockroach Cryptocercus meridianus (Blattodea: Cryptocercidae) and Its Phylogenetic Relationship among Cockroach Families

In this study, the complete mitochondrial genome of Cryptocercus meridianus was sequenced. The circular mitochondrial genome is 15,322 bp in size and contains 13 protein-coding genes, two ribosomal RNA genes (12S rRNA and 16S rRNA), 22 transfer RNA genes, and one D-loop region. We compare the mitogenome of C. meridianus with that of C. relictus and C. kyebangensis. The base composition of the whole genome was 45.20%, 9.74%, 16.06%, and 29.00% for A, G, C, and T, respectively; it shows a high AT content (74.2%), similar to the mitogenomes of C. relictus and C. kyebangensis. The protein-coding genes are initiated with typical mitochondrial start codons except for cox1 with TTG. The gene order of the C. meridianus mitogenome differs from the typical insect pattern for the translocation of tRNA-Ser^AGN, while the mitogenomes of the other two Cryptocercus species, C. relictus and C. kyebangensis, are consistent with the typical insect pattern. There are two very long non-coding intergenic regions lying on both sides of the rearranged gene tRNA-Ser^AGN. The phylogenetic relationships were constructed based on the nucleotide sequence of 13 protein-coding genes and two ribosomal RNA genes. The mitogenome of C. meridianus is the first representative of the order Blattodea that demonstrates rearrangement, and it will contribute to the further study of the phylogeny and evolution of the genus Cryptocercus and related taxa.

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.

Progress, pitfalls and parallel universes: a history of insect phylogenetics

The phylogeny of insects has been both extensively studied and vigorously debated for over a century. A relatively accurate deep phylogeny had been produced by 1904. It was not substantially improved in topology until recently when phylogenomics settled many long-standing controversies. Intervening advances came instead through methodological improvement. Early molecular phylogenetic studies (1985-2005), dominated by a few genes, provided datasets that were too small to resolve controversial phylogenetic problems. Adding to the lack of consensus, this period was characterized by a polarization of philosophies, with individuals belonging to either parsimony or maximum-likelihood camps; each largely ignoring the insights of the other. The result was an unfortunate detour in which the few perceived phylogenetic revolutions published by both sides of the philosophical divide were probably erroneous. The size of datasets has been growing exponentially since the mid-1980s accompanied by a wave of confidence that all relationships will soon be known. However, large datasets create new challenges, and a large number of genes does not guarantee reliable results. If history is a guide, then the quality of conclusions will be determined by an improved understanding of both molecular and morphological evolution, and not simply the number of genes analysed.

Sharp switches between regular and swinger mitochondrial replication: 16S rDNA systematically exchanging nucleotides A<->T+C<->G in the mitogenome of Kamimuria wangi

Swinger DNAs are sequences whose homology with known sequences is detected only by assuming systematic exchanges between nucleotides. Nine symmetric (X<->Y, i.e. A<->C) and fourteen asymmetric (X->Y->Z, i.e. A->C->G) exchanges exist. All swinger DNA previously detected in GenBank follow the A<->T+C<->G exchange, while mitochondrial swinger RNAs distribute among different swinger types. Here different alignment criteria detect 87 additional swinger mitochondrial DNAs (86 from insects), including the first swinger gene embedded within a complete genome, corresponding to the mitochondrial 16S rDNA of the stonefly Kamimuria wangi. Other Kamimuria mt genome regions are "regular", stressing unanswered questions on (a) swinger polymerization regulation; (b) swinger 16S rDNA functions; and (c) specificity to rDNA, in particular 16S rDNA. Sharp switches between regular and swinger replication, together with previous observations on swinger transcription, suggest that swinger replication might be due to a switch in polymerization mode of regular polymerases and the possibility of swinger-encoded information, predicted in primordial genes such as rDNA.

Comparison of the complete mitochondrial genome of the stonefly Sweltsa longistyla (Plecoptera: Chloroperlidae) with mitogenomes of three other stoneflies

The complete mitochondrial genome of the stonefly, Sweltsa longistyla Wu (Plecoptera: Chloroperlidae), was sequenced in this study. The mitogenome of S. longistyla is 16,151bp and contains 37 genes including 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a large non-coding region. S. longistyla, Pteronarcys princeps Banks, Kamimuria wangi Du and Cryptoperla stilifera Sivec belong to the Plecoptera, and the gene order and orientation of their mitogenomes were similar. The overall AT content for the four stoneflies was below 72%, and the AT content of tRNA genes was above 69%. The four genomes were compact and contained only 65-127bp of non-coding intergenic DNAs. Overlapping nucleotides existed in all four genomes and ranged from 24 (P. princeps) to 178bp (K. wangi). There was a 7-bp motif ('ATGATAA') of overlapping DNA and an 8-bp motif (AAGCCTTA) conserved in three stonefly species (P. princeps, K. wangi and C. stilifera). The control regions of four stoneflies contained a stem-loop structure. Four conserved sequence blocks (CSBs) were present in the A+T-rich regions of all four stoneflies.