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

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.

Comparative analysis of mitochondrial genomes among the family Peltoperlidae (Plecoptera: Systellognatha) and phylogenetic implications

Nowadays, the position of Peltoperlidae in Systellognatha has been resolved based on morphological analyses. However, there are different opinions based on molecular data. To date, only three peltoperlid mitogenomes are available, and more sampling is needed to obtain precise phylogenetic relationships. In this study, we obtained the complete mitogenomes of Cryptoperla kawasawai (15,832 bp) and Peltoperlopsis sagittata (15,756 bp). Our results show that gene content, gene order, DmTTF binding site, nucleotide composition, codon usage, ribonucleic acid (RNA) structure, and structural elements in the control region are highly conserved in peltoperlids. Heatmap analysis of codon usage shows that the AT-rich codons UUA, AUU, UUU, and AUA were commonly used codons in the Peltoperlidae. Evolutionary rate analyses of protein-coding genes reveal that different genes have been subject to different rates of molecular evolution correlated with the GC content. All tRNA genes in peltoperlid mitogenomes have a canonical cloverleaf secondary structure except for trnS1, whose dihydrouridine arm simply forms a loop. The control region of the family has several distinct structural characteristics and has the potential to serve as effective phylogenetic markers. Phylogenetic analyses support the monophyly of Perloidea, but the monophyly of Pteronarcyoidea is still not supported. The Peltoperlidae is placed as the earliest branch within the Systellognatha, and the estimated phylogenetic relationship is: Peltoperlidae + {(Styloperlidae + Pteronarcyidae) + [Perlidae + (Chloroperlidae + Perlodidae)]}. Our results provide new insight into the phylogeny of this group.

First Comprehensive Analysis of Both Mitochondrial Characteristics and Mitogenome-Based Phylogenetics in the Subfamily Eumeninae (Hymenoptera: Vespidae)

The subfamily Eumeninae plays a significant role in the biological control of agricultural pests. However, the characteristics of eumenine mitogenomes that are important molecular markers for phylogenetics are not clearly revealed. Here, 52 eumenine mitogenomes are newly sequenced and annotated, and the phylogenetic relationships of the subfamily are comprehensively analyzed based on 87 vespid mitogenomes. Through the comparative analysis of the 54 eumenine mitogenomes, the gene compositions of about one half of the 54 species match with ancestral insect mitogenome, and remaining others contain two trnM which are highly similar, with 51.86% (Eumenes tripunctatus) to 90.65% (Pseumenes nigripectus) sequence identities, which is unique among the reported mitogenomes of the family Vespidae. Moreover, the translocation trnL1 upstream of nad1 is a common rearrangement event in all eumenine mitogenomes. The results of phylogenetic analyses support the paraphyly of the subfamily Eumeninae and the tribe Odynerini, respectively, and the monophyly of the tribe Eumenini, and verify that the tribe Zethini is a valid subfamily Zethinae. In this study, the relationships between some genera such as Allorhynchium and Pararrhynchium or the taxonomic status of the subgenera such as Eremodynerus and Dirhynchium are found to be confusing and there should be further inquiry with more samples.

Mitochondrial composition of and diffusion limiting factors of three social wasp genera Polistes, Ropalidia, and parapolybia (Hymenoptera: Vespidae)

Background

Social wasps Polistes, Ropalidia, and Parapolybia, belonging to the subfamily Polistinae, have obviously different distribution patterns, yet the factors leading to this difference remain unknown.

Results

The 17 newly sequenced mitogenomes of Polistes, Ropalidia, and Parapolybia contain 37 genes, and there are obvious differences among the compositions of the three genera. The monophyly of the genus Polistes and a monophyletic Ropalidiini: (Ropalidia + Parapolybia) are concordant with previous morphological analysis of the subfamily Polistinae. Our inferred divergence time demonstrates Polistes (at around 69 Ma) was diverged earlier than Ropalidia and Parapolybia (at around 61 Ma). The rearrangement of both trnY and trnL1 are shared by all the Polistinae. In addition, the unique rearrangement of TDRL derived at 69 Ma is detected in Polistes, and Ropalidia contains a Reversal which may derive at 61 Ma. Hereafter, the possibility is elaborated that Polistes originated in Aisa and then dispersed from Africa to South America, and Polistes and Ropalidia spread from Southeast Asia to Australia. At last, continental drift and Quaternary Ice Ages are inferred to be two main limiting factors in the current distributions of the three genera.

Conclusions

Obvious differences occur in the mitochondrial composition of Polistes, Ropalidia, and Parapolybia. According to the reconstructed time-calibrated framework, it is inquired that the continental drifts and the climate are mainly diffusion limiting factors of the three genera.

Supplementary information

The online version contains supplementary material available at 10.1186/s12862-022-02017-6.

A phylogenomic perspective on the evolutionary history of the stonefly genus Suwallia (Plecoptera: Chloroperlidae) revealed by ultraconserved genomic elements

Evolutionary biologists have long sought to disentangle phylogenetic relationships among taxa spanning the tree of life, an increasingly important task as anthropogenic influences accelerate population declines and species extinctions, particularly in insects. Phylogenetic analyses are commonly used to identify unique evolutionary lineages, to clarify taxonomic designations of the focal taxa, and to inform conservation decisions. Advances in DNA sequencing techniques have increasingly facilitated the ability of researchers to apply genomic methods to phylogenetic analyses, even for non-model organisms. Stoneflies are non-model insects that are important bioindicators of the quality of freshwater habitats and landscape disturbance as they spend the immature stages of their life cycles in fresh water, and the adult stages in terrestrial environments. Phylogenetic relationships within the stonefly genus Suwallia (Insecta: Plecoptera: Chloroperlidae) are poorly understood, and have never been assessed using molecular data. We used DNA sequence data from genome-wide ultraconserved element loci to generate the first molecular phylogeny for the group and assess its monophyly. We found that Palearctic and Nearctic Suwallia do not form reciprocally monophyletic clades, and that a biogeographic history including dispersal, vicariance, and founder event speciation via jump dispersal best explains the geographic distribution of this group. Our results also strongly suggest that Neaviperla forcipata (Neave, 1929) is nested within Suwallia, and the concept of the genus Suwallia should be revised to include it. Thus, we formally propose a new taxonomic combination wherein Neaviperla forcipata (Neave, 1929) is reclassified as Suwallia forcipata (Neave, 1929). Moreover, some Suwallia species (e.g., S. amoenacolens, S. kerzhneri, S. marginata, S. pallidula, and S. starki) exhibit pronounced cryptic diversity that is worthy of further investigation. These findings provide a first glimpse into the evolutionary history of Suwallia, improve our understanding of stonefly diversity in the tribe Suwallini, and highlight areas where additional research is needed.

Two Complete Mitochondrial Genomes From Leuctridae (Plecoptera: Nemouroidea): Implications for the Phylogenetic Relationships Among Stoneflies

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.

Characterization of the Complete Mitochondrial Genome of Harpalus sinicus and Its Implications for Phylogenetic Analyses

In this study, we report the complete mitochondrial genome of Harpalus sinicus (occasionally named as the Chinese ground beetle) which is the first mitochondrial genome for Harpalus. The mitogenome is 16,521 bp in length, comprising 37 genes, and a control region. The A + T content of the mitogenome is as high as 80.6%. A mitochondrial origins of light-strand replication (OL)-like region is found firstly in the insect mitogenome, which can form a stem-loop hairpin structure. Thirteen protein-coding genes (PCGs) share high homology, and all of them are under purifying selection. All tRNA genes (tRNAs) can be folded into the classic cloverleaf secondary structures except tRNA-Ser (GCU), which lacks a dihydrouridine (DHU) stem. The secondary structure of two ribosomal RNA genes (rRNAs) is predicted based on previous insect models. Twelve types of tandem repeats and two stem-loop structures are detected in the control region, and two stem-loop structures may be involved in the initiation of replication and transcription. Additionally, phylogenetic analyses based on mitogenomes suggest that Harpalus is an independent lineage in Carabidae, and is closely related to four genera (Abax, Amara, Stomis, and Pterostichus). In general, this study provides meaningful genetic information for Harpalus sinicus and new insights into the phylogenetic relationships within the Carabidae.

Tick mitochondrial genomes: structural characteristics and phylogenetic implications

Ticks are obligate blood-sucking arachnid ectoparasites from the order Acarina, and many are notorious as vectors of a wide variety of zoonotic pathogens. However, the systematics of ticks in several genera is still controversial. The mitochondrial genome (mt-genome) has been widely used in arthropod phylogeny, molecular evolution and population genetics. With the development of sequencing technologies, an increasing number of tick mt-genomes have been sequenced and annotated. To date, 63 complete tick mt-genomes are available in the NCBI database, and these genomes have become an increasingly important genetic resource and source of molecular markers in phylogenetic studies of ticks in recent years. The present review summarizes all available complete mt-genomes of ticks in the NCBI database and analyses their characteristics, including structure, base composition and gene arrangement. Furthermore, a phylogenetic tree was constructed using mitochondrial protein-coding genes (PCGs) and ribosomal RNA (rRNA) genes from ticks. The results will provide important clues for deciphering new tick mt-genomes and establish a foundation for subsequent taxonomic research.

The complete mitochondrial genome of a stonefly species, Suwallia bimaculata (Plecoptera: Chloroperlidae)

The stonefly Suwallia bimaculata belongs to the Chloroperlidae. The mitogenome of S. bimaculata was sequenced and annotated, the new representative of the complete mitogenome of the genus Suwallia. The entire genome of S. bimaculata is 16,125 bp totally with an A + T content of 68.5%, consisting of 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNAs, and a 1,210 bp control region. All genes have the similar locations and strands with that of other published species of Plecoptera. Most PCGs use typical start/stop codon, whereas ND5 use GTG as start codon. And, COII and ND5 stopped with incomplete terminaton signal T. Phylogenetic analysis was performed using the Bayesian (BI) method and Maximum-likelihood (ML) methods based on 13 protein-coding genes and two ribosomal RNAs showed that S. bimaculata was the sister group to Suwallia teleckojensis and the clade Suwallia was closely to four perlodid species. More chloroperlid data was needed for further study of phylogeny in Chloroperlidae.

Tracking outbreak populations of the pepper weevil Anthonomus eugenii (Coleoptera; Curculionidae) using complete mitochondrial genomes

The pepper weevil, Anthonomus eugenii, is a major pest on Capsicum species. Apart from natural spread, there is a risk of spread via international pepper trade. In the Netherlands, a pepper weevil outbreak occurred in 2012 and affected six greenhouses producing different sweet pepper varieties. The following year, a pepper weevil outbreak occurred in Italy. To trace the origin of the Dutch outbreak and to establish if the Dutch and Italian outbreaks were linked, we determined the mitogenomes of A. eugenii specimens collected at outbreak locations, and compared these with specimens from the native area, and other areas where the pest was introduced either by natural dispersal or via trade. The circular 17,257 bp A. eugenii mitogenome comprises thirteen mitochondrial genes typically found in insect species. Intra-species variation of these mitochondrial genes revealed four main mitochondrial lineages encompassing 41 haplotypes. The highest diversity was observed for specimens from its presumed native area (i.e. Mexico). The Dutch outbreak specimens represented three highly similar haplotypes, suggesting a single introduction of the pest. The major Dutch haplotype was also found in two specimens from New Jersey. As the Netherlands does not have pepper trade with New Jersey, it is likely that the specimens sampled in New Jersey and those sampled in the Netherlands originate from a shared source that was not included in this study. In addition, our analysis shows that the Italian and Dutch outbreaks were not linked. The mitochondrial genome is a useful tool to trace outbreak populations and the methodology presented in this paper could prove valuable for other invasive pest species, such as the African fruit moth Thaumatotibia leucotreta and emerald ash borer Agrilus planipennis.

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 first mitochondrial genome from Scopuridae (Insecta: Plecoptera) reveals structural features and phylogenetic implications

This study determined the first complete mitochondrial genome (mitogenome) of a stonefly, Scopura longa (Plecoptera: Scopuridae), and reconstructed a phylogeny based on two datasets of mitogenomes in eighteen available stoneflies to examine the relationships among Plecoptera. The complete mitogenome of S. longa is a circular molecule of 15,798bp in size. It contains 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a control region (CR). Most PCGs used standard ATN start codons and TAN stop codons. All tRNAs could be folded as typical cloverleaf secondary structures except tRNA^Ser(AGN), with the dihydrouridine (DHU) arm missing. Tandem repeat regions, two potential stem-loop (SL) structures, Poly T structures and G+C-rich regions are detected in the control region. Finally, the phylogenetic relationships among the families within the Arctoperlaria were reconstructed. The topological structures of the two trees were almost identical. The present phylogenetic analysis shows that S. longa belongs to the infraorder Euholognatha. The monophyly of each family is generally well supported based on nucleotide sequences. The Pteronarcyidae is sister to ((Peltoperlidae+Styloperlidae)+(Perlidae+(Perlodidae+Chloroperlidae))).