The complete mitochondrial genome of Amphinemura bulla Shimizu, 1997 (Plecoptera: Nemouridae) from Japan

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.

The First Complete Mitochondrial Genome of Genus Isocapnia (Plecoptera: Capniidae) and Phylogenetic Assignment of Superfamily Nemouroidea

Capniidae are a family of stoneflies, also known as snow flies, who emerge in winter. The phylogeny of Capniidae is widely accepted to be based on morphological analysis. Until now, only five Capniidae mitochondrial genomes have been sequenced so far. In addition, sampling is required to determine an accurate phylogenetic association because the generic classification of this family is still controversial and needs to be investigated further. In this study, the first mitogenome of genus Isocapnia was sequenced with a length of 16,200 bp and contained 37 genes, including a control region, two rRNAs, 22 tRNAs, and 13 PCGs, respectively. Twelve PCGs originated with the common start codon ATN (ATG, ATA, or ATT), while nad5 used GTG. Eleven PCGs had TAN (TAA or TAG) as their last codon; however, cox1 and nad5 had T as their final codon due to a shortened termination codon. All tRNA genes demonstrated the cloverleaf structure, which is distinctive for metazoans excluding the tRNASer1 (AGN) that missed the dihydrouridine arm. A Phylogenetic analysis of the superfamily Nemouroidea was constructed using thirteen PCGs from 32 formerly sequenced Plecoptera species. The Bayesian inference and maximum likelihood phylogeny tree structures derived similar results across the thirteen PCGs. Our findings strongly supported Leuctridae + ((Capniidae + Taeniopterygidae) + (Nemouridae + Notonemouridae)). Ultimately, the best well-supported generic phylogenetic relationship within Capniidae is as follows; (Isocapnia + (Capnia + Zwicknia) + (Apteroperla + Mesocapnia)). These findings will enable us to better understand the evolutionary relationships within the superfamily Nemouroidea and the generic classification and mitogenome structure of the family Capniidae.

The mitochondrial genome analysis of a stonefly, Nemoura longicercia Okamoto, 1922 (Plecoptera: Nemouridae)

To better understand the diversity and phylogeny of Nemouridae, we sequenced and annotated the complete mitochondrial genome (mitogenome) of Nemoura longicercia Okamoto, 1922. The entire mitochondrial genome of N. longicercia was 15,728 bp. It contained 37 representative genes (22 transfer RNA genes (tRNAs), 13 protein-coding genes (PCGs), and two ribosomal RNA genes (rRNAs)) and a control region. In this genome, 23 genes were found on the heavy strand (H-strand) and 14 genes were on the light strand (L-strand). All PCGs began with the initiation codon ATN. Eleven PCGs stopped with the termination codon TAA or TAG, while COII and ND5 genes stopped with incomplete codon T. Bayesian’s inference (BI) and maximum-likelihood (ML) methods generated the identical tree topology across the PCGR dataset (13 PCGs plus two rRNAs). The monophyly of each subfamily was well-supported. This analysis supported the clade N. longicercia plus N. papilla as sister taxon to the clade N. meniscata plus N. nankinensis.

Mitochondrial Genome of Strophopteryx fasciata (Plecoptera: Taeniopterygidae), with a Phylogenetic Analysis of Nemouroidea

Taeniopterygidae is a medium-sized family of stoneflies. The phylogeny of Taeniopterygidae was widely accepted based on the morphological analyses. However, there are different opinions based on molecular data. To date, only two taeniopterygid mitochondrial genomes (mitogenomes) were available, and more sampling is needed to obtain precise phylogenetic relationships. In this research, the Strophopteryx fasciata mitogenome was sequenced and analyzed. The complete mitogenome of S. fasciata was 15,527 bp in length and contained 37 genes and a non-coding control region. Among taeniopterygid mitogenomes, the length variation was minimal in protein-coding genes (PCGs), transfer RNA genes (tRNAs) and ribosomal RNA genes (rRNAs), but very different in the control region. Similar to mitogenomes of other taeniopterygid species, the S. fasciata mitogenome was consistently AT biased and displayed positive AT- and negative GC-skews of the whole mitogenome. Most PCGs used ATN as the start codon and TAA/TAG as the stop codon. The stop codons were far less variable than the start codons in taeniopterygid mitogenomes. All Ka/Ks ratios were less than 1, indicating the presence of purifying selection in these genes. The secondary structures of transfer and ribosomal RNA genes of S. fasciata mitogenome are highly conserved with other taeniopterygid species. In the control region of the S. fasciata mitogenome, some essential elements (tandem repeats, stem–loop structures, and poly−N stretch, etc.) were observed. Two phylogenetic trees were inferred from Bayesian inference (BI) and Maximum Likelihood (ML) methods generated the identical topology across the PCGR dataset. The relationships of five families in Nemouroidea were recovered as Leuctridae + ((Capniidae + Taeniopterygidae) + (Nemouridae + Notonemouridae)). These results will help us understand the mitogenome structure of taeniopterygid species and the evolutionary relationship within Plecoptera.

Complete Mitochondrial Genome of Malenka flexura (Plecoptera: Nemouridae) and Phylogenetic Analysis

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