The complete mitochondrial genome of Tenguna medogensis (Hemiptera: Dictyopharidae)

In this study, we sequenced and annotated the complete mitochondrial genome of Tenguna medogensis Song & Liang, 2007 (Hemiptera: Dictyopharidae). The mitogenome was 15,400 bp in length with a high A + T content (76.95%), containing 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and a control region. All protein-coding genes started with ATC/ATG/ATT, except for nad1 and nad5, which used GTG. The predicted secondary structures of all tRNA genes were typical cloverleaf except for trnS (gcu) and trnV (uac), lacking the dihydrouridine stem. In the maximum-likelihood phylogenetic result of planthoppers, Tenguna was closely related to Metaurus in Dictyopharidae.

Mitogenomic Characterization and Comparative Analysis of Three Egg Parasitoid Wasps Parasitizing Nilaparvata lugens (Stål)

This study reports the first complete mitochondrial genomes of three egg parasitoid wasps parasitizing Nilaparvata lugens-Pseudoligosita nephotetticum, Anagrus frequens, and Anagrus nilaparvatae. Genome sizes ranged from 15,429 to 15,889 bp, with all three mitogenomes displaying strong A + T bias, standard gene content, and characteristic strand asymmetries. While A. frequens and A. nilaparvatae exhibited conserved gene orders, extensive gene rearrangements, including multiple inversions in both protein-coding genes (PCGs) and tRNAs, were observed in P. nephotetticum. Codon usage analyses revealed a preference for codons ending in A or U. The non-synonymous (Ka) to synonymous (Ks) substitution ratio analysis identified signs of positive selection in multiple PCGs, particularly in atp8, nad6, and nad3, suggesting possible adaptive evolution related to host-searching behavior. Secondary structure analyses showed the loss of trnL1 in all Anagrus species, while trnS1 and trnR lacked the DHU arm, indicating possible derived traits in Mymaridae. Phylogenetic analysis was the first time to describe the relationship of the genus Anagrus within Mymaridae from the perspective of 13 protein genes. Furthermore, the grouping of (Pseudoligosita + Megaphragma) + Trichogramma was supporting the distinct evolutionary lineage of Pseudoligosita. This work provides new molecular resources and phylogenetic insight for Chalcidoidea, with implications for parasitoid evolution and biological control strategies.

Phylogenetic Analyses of Bostrichiformia and Characterization of the Mitogenome of Gibbium aequinoctiale (Bostrichiformia Ptinidae)

BACKGROUND

Ptinidae, within the infraorder Bostrichiformia, are a cosmopolitan, ecologically diverse but poorly known group. The phylogeny within Bostrichiformia and the monophyly of Ptinidae and its phylogenetic placement in Bostrichiformia remain contentious.

METHODS

In this research, we determined the entire mitochondrial genome (mitogenome) of Gibbium aequinoctiale, the first representative mitogenome of the subfamily Ptininae, and reconstructed the phylogenetic relationships for Bostrichiformia based on four mitochondrial datasets using maximum likelihood (ML) and Bayesian inference (BI) methods.

RESULTS

The mitogenome of G. aequinoctiale is a circular molecule spanning 17,020 bp and harbors 37 mitochondrial genes and a presumed control region (CR). The mitogenome exhibited a marked preference for the utilization of A and T bases, which was also observed in three kinds of genes and CR. AAT was inferred as the putative candidate initiation codon for cytochrome oxidase subunits 1 (COI). The control region contains three tandem repeats (TDRs) and one poly-thymine stretch (Poly-T) in both coding strands. The phylogenetic results appeared to support the monophyly of four families, Nosodendridae, Derodontidae, Dermestidae, and Bostrichidae, and the basal position of the latter two families within Bostrichiformia. However, the family Ptinidae was not verified as monophyly because of one species diverging from the main lineage. Three families, Dermestidae, Bostrichidae, and Ptinidae, clustered as the major clade in Bostrichiformia, among which Bostrichidae and Ptinidae grouped together as sister groups.

CONCLUSIONS

The present study provides valuable mitochondrial information for Ptinidae and provides novel perspectives on the inner phylogeny within the infraorder Bostrichiformia.

Comparative Analysis of Metopograpsus quadridentatus (Crustacea: Decapoda: Grapsidae) Mitochondrial Genome Reveals Gene Rearrangement and Phylogeny

The taxonomy of the genus Metopograpsus is still a matter of debate because its species have limited morphological differences. Mitochondrial genomes, which are highly informative and easily accessible genetic markers, have been widely used to study molecular evolution and clarify relationships among species. This study aims to investigate the mitochondrial genome of Metopograpsus quadridentatus, a species with unique ecological and evolutionary significance. By analyzing the mitochondrial genome, we seek to address taxonomic uncertainties and provide insights into the evolutionary history of this species. In this study, we sequenced and analyzed the mitochondrial genome of M. quadridentatus to investigate its gene rearrangement patterns and its place within Brachyura. We compared the gene order of M. quadridentatus with that of 40 other Brachyuran species and created phylogenetic trees based on the nucleotide and amino acid sequences of 13 protein-coding genes (PCGs). We found that the mitochondrial gene arrangement of M. quadridentatus is mostly unchanged, similar to the original crustacean pattern, except for the movement of the trnH gene. Notably, the gene orders of several families, such as Eriphiidae, Grapsidae, Camptandriidae, Dotillidae, Plagusiidae, Ocypodidae, and Gecarcinidae, are the same. Phylogenetic analyses consistently placed M. quadridentatus within the genus Metopograpsus and the family Grapsidae, confirming its current taxonomic classification. These results offer important insights into evolutionary relationships and gene order conservation within Brachyura. Our study highlights the significance of mitochondrial genomes in resolving taxonomic uncertainties within the genus Metopograpsus.

Fly in the Ointment: Host-Specificity Challenges for Botanophila turcica, a Candidate Agent for the Biological Control of Saffron Thistle in Australia

In classical biological control of weeds, the risk posed by a candidate agent to close relatives of the target weed in the intended area of release is a key criterion (i.e., candidate agents that demonstrate a high degree of host specificity). In this study, we investigated if the rosette crown-feeding fly Botanophila turcica Hennig (Diptera: Anthomyiidae) could meet this criterion and thus be considered a good candidate to control saffron thistle Carthamus lanatus L. (Asteraceae: Cardueae) in Australia. Previous studies indicated that B. turcica is specific to Ca. lanatus and did not infest the closely related crop, safflower (Carthamus tinctorius L.). However, more recent field observations made in Greece reported that B. turcica infested safflower in cultivated fields. To determine if B. turcica is safe for release as a biocontrol agent, we re-examined the host range of B. turcica by performing new host-specificity testing combined with field surveys carried out in the south of France during two consecutive years. We also investigated the species identity of the flies by comparing DNA sequences (COI barcode region) of specimens collected in France from Ca. lanatus and Centaurea solstitialis L. with those from Greece collected from Ce. solstitialis and Centaurea diffusa Lam. Our COI analyses confirmed that French and Greek samples identified as B. turcica belonged to the same species, while a second group of Greek samples matched B. brunneilinea, indicating two distinct species. Our results also demonstrated that B. turcica has a wider host range than previously suggested. Laboratory testing indicated that Ca. lanatus, Ca. tinctorius, and Ce. solstitialis are suitable for the development of B. turcica. Field surveys also revealed that Ce. diffusa is part of the host range of the fly. Based on the results reported here, B. turcica may have the potential to control both the target weed, Ca. lanatus, and Ce. Solstitialis, but it may also be a threat to safflower, Ca. tinctorius. Further investigations to assess under what conditions B. turcica attacks Ca. tinctorius may help clarify the level of risk to Australian growers.

Mitochondrial Genome of Grapsus albolineatus and Insights into the Phylogeny of Brachyura

Brachyura is among the most diverse groups of crustaceans, with over 7000 described species. Crab mitogenomes are important for understanding molecular evolution and phylogenetic relationships. Grapsus albolineatus exhibits specific rearrangements compared with the Pancrustacean ground pattern and other Brachyura species. The gene arrangement of G. albolineatus is similar to that of ancestral crustaceans, barring that of the translocated trnH gene. In phylogenetic analyses, the Bayesian inference estimation was observed to be superior to the maximum likelihood estimation when the nodal support values were compared. Considering the results of the gene rearrangement pattern and phylogenetic analysis, we speculate that G. albolineatus belongs to Grapsidae. Our comparative study indicated that mitogenomes are a useful phylogenetic tool at the subfamily level within Brachyura. The findings indicate that mitogenomes could be a useful tool for systematics in other Brachyuran species.

The Palaearctic types of Chrysididae (Insecta, Hymenoptera) deposited in the Linsenmaier collection. Part 5. Elampini: genus Holopyga Dahlbom, 1845

An illustrated catalogue of the Palaearctic types of species and subspecies of Holopyga Dahlbom, 1845 described by Walter Linsenmaier is presented. Linsenmaier described 48 species and subspecies of Holopyga and most of the primary types are deposited in the collection of the Natur-Museum in Luzern, Switzerland, with the exception of five holotypes which are deposited in other public and private collections. This catalogue includes the list of the described species, photographs of 41 holotypes and a synthesis of the recent changes in the taxonomy of Holopyga to update Linsenmaier's classification. New synonymies are proposed for the subgenus Chamaeholopyga Linsenmaier, 1987 syn. nov. of Pseudolopyga Krombein, 1969 and for two species: Holopyga unitasculpta Linsenmaier, 1968 syn. nov. of H. solskyi (Radoszkowski, 1877) and H. jurinei Chevrier, 1862 syn. nov. of H. lucida (Lepeletier, 1806). Additionally, a new species, Holopyga dichroica Rosa sp. nov., formally referred to as H. jurinei sensu Linsenmaier 1959a is described.

Molecular Phylogenetic Relationships Based on Mitogenomes of Spider: Insights Into Evolution and Adaptation to Extreme Environments

In this study, we performed a comparative analysis based on a total of 255 spider mitogenomes and four outgroups, of which the mitogenomes of 39 species were assembled de novo, to explore the phylogenetic relationships and the adaptive evolution of mitogenomes. Results showed that Argyroneta aquatica had the longest mitochondrial length and the most pronounced codon preference to be UUA, followed by CCU. Codon usage frequencies were similar between families and codon usage in the mitogenome of spiders was mainly influenced by natural selection pressures rather than G/C mutation bias. Our phylogenetic topology clearly explained the evolutionary relationships among the spiders, and divergence time estimates indicated that the spiders originated in the early Devonian, and that the two clades of Mesothelae and Opisthothelae separated in the late Carboniferous. Ancestral range and trait reconstruction results supported the ancestral origin of spiders to the Devonian Nearctic realm, with the trapdoor being the original trait. Selection analysis detected positive selection signals in the ATP8 gene in Desis jiaxiangi. The ND5 gene is a convergent evolutionary gene between D. jiaxiangi and A. aquatica. Positive selection signals in the ATP8 gene and convergent selection sites in the ND5 gene may facilitate metabolic adaptation to the aquatic environment in two aquatic spiders. In conclusion, our analysis contributes to a better understanding of the taxonomic status, species diversity, mitochondrial characteristics, and environmental adaptations of these spiders.

Phylogeny and divergence time estimation of the subfamily Amphipsyllinae based on the Frontopsylla diqingensis mitogenome

Fleas are primarily parasites of small mammals and serve as essential vectors of the transmission of plague. The subfamily Amphipsyllinae (Siphonaptera: Leptopsyllidae) consists of 182 species across 13 genera, widely distributed worldwide. Only two species of Amphipsyllinae have been sequenced for complete mitogenomes to date. It hinders the taxonomy and evolutionary history studies of fleas. In this study, we first sequenced the Frontopsylla diqingensis mitogenome and performed comparative mitogenomic analyses with the two other species (Frontopsylla spadix and Paradoxopsyllus custodis) in Amphipsyllinae available in the NCBI database. The evolutionary process of Amphipsyllinae was comprehensively analyzed in terms of nucleotide composition, codon usage, nucleotide diversity, tRNA secondary structure, nucleotide skew, phylogeny tree, and divergence time. Nucleotide diversity and tRNAs of three species of fleas of Amphipsyllinae have differences among different species. The effective number of codon (ENC)-plot, neutrality curve, PR2, and correspondence analysis (COA) showed that the codon preference of Amphipsyllinae was influenced mainly by natural selection. For phylogenetic trees and divergence time of the order Siphonaptera, our results showed two concatenated data matrices, namely, PCG: (((Ceratophyllidae + Leptopsyllidae) + ((Vermipsyllidae + Hystrichopsyllidae) + Ctenophthalmidae)) + (Pulicidae + Pygiopsyllidae)); PCGRNA: ((((Ceratophyllidae + Leptopsyllidae) + ((Vermipsyllidae + Hystrichopsyllidae) + Ctenophthalmidae)) + Pulicidae) + Pygiopsyllidae). We concluded that P. custodis and Macrostylophora euteles from GenBank are the same species by phylogenetic trees and sequence alignment, and supported the monophyly of Amphipsyllinae. Amphipsyllinae diverged in the Cenozoic, approximately 73.37-40.32 million years ago (Mya). The majority of the species within the intraordinal divergence into extant lineages occurred after the K-Pg boundary. The common ancestor of the extant order Siphonaptera diverged during the Cretaceous. Our findings supported those of Zhu et al. (1). This study provides new insights into the evolutionary history and taxonomy of the order Siphonaptera.

Toward the Development of the Trojan Female Technique in Pest Insects: Male-Specific Influence of Mitochondrial Haplotype on Reproductive Output in the Seed Beetle Acanthoscelides obtectus

Biocontrol techniques that impair reproductive capacity of insect pests provide opportunities to control the dynamics of their populations while minimizing collateral damage to non-target species and the environment. The Trojan Female Technique, or TFT, is a method of the trans-generational fertility-based population control through the release of females that carry mitochondrial DNA mutations that negatively affect male, but not female, reproductive output. TFT is based on the evolutionary hypothesis that, due to maternal inheritance of mitochondria, mutations which are beneficial or neutral in females but harmful in males can accumulate in the mitochondrial genome without selection acting against them. Although TFT has been theoretically substantiated, empirical work to date has focused only on Drosophila melanogaster populations, while the existence of male-biased mutations and the TFT approach in economically important pest species remain unexplored. Here, we examined the sex-specific effects of three distinct and naturally occurring mitochondrial haplotypes (MG1a, MG1d, and MG3b) on several reproductive and life history traits in the seed beetle Acanthoscelides obtectus. Our results revealed that males harboring the MG3b mitotype exhibited lower early fecundity and fertility, while there were no effects on females or longevity in either sex. Our experiments provide support for the existence of the mitochondrial variant that specifically impairs male reproductive output in pest insects. These results can be harnessed to further develop TFT as a novel form of biocontrol with broad applicability to economic pests and disease vector insects.

New Mitochondrial Genomes of Ithonidae (Neuroptera) and Higher Phylogenetic Implications

Ithonidae (moth lacewings) are an enigmatic, small family of the insect order Neuroptera (lacewings). Its phylogenetic position within Neuroptera and internal subfamily relationships remain unresolved. In this study, the complete mitochondrial genome (mitogenome) of Ithone fulva Tillyard, 1916 representing the first mitogenome of Ithoninae, as well as the complete mitogenome of Rapisma gaoligongensis Liu, Li and Yang, 2018, were newly reported. Molecular phylogenetic trees recovered Ithonidae as the sister group to Chrysopidae + Hemerobiidae. Ithoninae was demonstrated to be the sister group to Polystoechotinae + Rapismatinae across all topologies. Divergence time estimation revealed that Ithonidae originated during the Mid-Triassic. Ithoninae diverged from Polystoechotinae and Rapismatinae in the Late Triassic, while Polystoechotinae diverged from Rapismatinae in the Middle Jurassic. Moreover, more comprehensive samplings, as well as whole genome data, are needed to reconstruct a comprehensive framework of the phylogeny of Ithonidae, as well as Neuropteida.

The first complete mitochondrial genome of Loimia arborea (Polychaeta: Terebellidae) and phylogenetic analysis

In order to understand the molecular insights within the Terebellida, the complete mitochondrial genome of Loimia arborea was sequenced. The mitochondrial genome is 16,023 bp, with 13 protein-coding genes (PCGs), 23 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and a non-coding region (D-loop). Notably, two adjacent copies of methionine tRNA genes (trnMs) were detected. The phylogeny of Terebellida was constructed based on 13 PCGs from 13 species, two main clades were strongly supported, i.e., Cirratuliformia (clade A) and Terebelliformia (clade B). And the monophyly of the Terebellidae was restored.

The partial mitochondrial genome of Semidalis anchoroides Liu & Yang, 1993 (Neuroptera: Coniopterygidae)

Semidalis anchoroides Liu & Yang, 1993 is a small and common insect in southern China. It's known for its small size, being covered by whitewax powder, and having uncomplicated venation. In this study, the mitochondrial genome for S. anchoroides was sequenced and analyzed. The sequenced partial mitogenome is 15,700 bp in length, encoding 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and one partial control region. These findings provide fundamental molecular data, thereby facilitating a more profound understanding of the phylogenetic relationships within the Semidalis species.

The complete mitochondrial genome of Dimorphopterus japonicus (Hidaka, 1959) (Hemiptera, Lygaeoidea) and phylogenetic relationships within the Lygaeoidea superfamily

Dimorphopterus japonicus is a kind of pest which seriously harms sorghum and millet. In this study, we sequenced the mitochondrial whole genome of D. japonicus (Hidaka, 1959), and characterized and analyzed the mitogenome. The D. japonicus genome length is 15,473 bp, and it exhibits a typical high A + T content (77.6%). The mitogenomic structure of D. japonicus is highly conservative and there are no gene rearrangements. The evolutionary rates of the PCGs are in the order of atp8 > nad3 > nad6 > nad2 > nad5 > atp6 > nad4 > nad4L > cox2 > nad1 > cytb > cox3 > cox1. By using the Bayesian inference and maximum likelihood methods, we inferred the phylogenetic relationships within Lygaeoidea and estimated their divergence times based on concatenated mitogenome genes. The stable clades of (Oxycarenidae + Piesmatidae), (Malcidae + Colobathristidae), (Meschiidae + Berytidae), and (Blissidae+(Cymidae + Ninidae)) were consistently recovered in all analyses. Estimated divergence times revealed that the divergence time of the Lygaeoidea was 99.4 Ma (95% HPD: 83.4-124.8 Ma), and most lygaeoid families diverged during the early Cretaceous to late Paleogene. Berytidae is the latest differentiation in the lygaeoid families, and the genus Dimorphopterus was differentiated in 45.9 Ma (95% HPD: 39.8-47.7 Ma). This study is of great significance for reconstructing the phylogeny of Lygaeoidea and providing insights into its evolutionary history.

Phylogenetic and divergence analysis of Pentatomidae, with a comparison of the mitochondrial genomes of two related species (Hemiptera, Pentatomidae)

Pentatomidae, the most diverse family of Pentatomoidea, is found worldwide. Currently, the phylogenetic relationships among Pentatomidae tribes remain unstable, and subfamily divergence has not been estimated. Here, we sequenced and analyzed the complete mitochondrial genomes of two species of Lelia, and studied the phylogenetic relationships among Pentatominae tribes. We also selected three available fossil as the calibration points in the family, and preliminarily discussed the divergence time of Pentatomidae. Trees of Pentatomidae were reconstructed using the Bayesian inference method. Divergence times of Pentatominae were estimated based on the nucleotide sequences of protein-coding genes with a relaxed clock log-normal model in BEASTv.1.8.2. The results showed that the gene arrangements, nucleotide composition, and codon preferences were highly conserved in Lelia. Further, a phylogenetic analysis recovered Eysarcorini, Strachiini, Phyllocephalini, and Menidini as monophyletic with strong support, however, the monophyly of Antestiini, Nezarini, Carpocorini, Pentatomini and Cappaeini were rejected. Moreover, Pentatominae diverged from Pentatomidae soon after the origin of the Cretaceous Period, at approximately 110.38 Ma. This study enriches the mitochondrial genome database of Pentatomidae and provides a reference for further phylogenetic studies, and provides a more accurate estimate of divergence time.

Mitogenome selection shaped the terrestrial adaptation of Grapsidae (Decapoda: Brachyura)

The colonization of aquatic to terrestrial habitats by brachyuran crabs requires genetic innovations as well as morphological adaptations to adapt to terrestrial environments. The genetic basis of such adaptive evolution, however, is largely unknown. This study focuses on terrestrialization in Geograpsus (Grapsidae) the only highly terrestrial genus in this family, which represents a notable example of terrestrial adaptive radiation. Here, we sequenced the mitogenomes of two Geograpsus species and used the mitogenomes of 215 representative crabs to construct phylogenetic and time frameworks that we used to infer terrestrial origins and evolution. Using mitochondrial genomic data, we demonstrated that marine crab ancestors began to settle on land during the early Eocene. Ocean acidification, the Paleocene-Eocene Thermal Maximum (PETM), and mangrove expansion at that time may have driven the diversification and ecological expansion of these terrestrial crabs. Evolutionary analyses reveal strong positive selection signals on monophyletic lineages of Grapsidae, especially the terrestrial species of Geograpsus. Positively selected sites in functionally important regions of ND5 and ND4 may imply enhanced energy metabolism in Grapsidae compared to other crabs, and may have played an important role in their terrestrial adaptation. Overall, our work provides valuable resources and opportunities to reveal the adaptation of crabs to complex terrestrial environments.

The biogeographic patterns of the olive fly and its primary symbiont Candidatus Erwinia dacicola across the distribution area of the olive tree

The olive fly, Bactrocera oleae (Rossi, 1790), is the major insect pest of olives attacking both cultivated and wild olive. Bactrocera oleae carries a primary and vertically transmitted symbiont, the bacterium Candidatus Erwinia dacicola. As any primary symbiont, it plays an important role in the reproduction and lifespan of the fly. The genetic 16S rRNA diversity of the primary symbiont and the mitochondrial haplotype variation of the insect host were simultaneously examined in 54 olive fly populations. The aim was to unravel the biogeographic patterns of this economically relevant host-bacteria interaction across a wide distribution area. Three symbiont haplotypes were identified. The primary symbiont showed a lower haplotype diversity than that of its host, a characteristic indicative of a long-term interaction. A significant genetic and geographic association between host and primary symbiont was observed, with an East-West genetic differentiation pattern in the Mediterranean basin, coinciding with the historical genetic distribution of the olive tree. The study shows promise, informing and aiding the development of future tools for the control of the olive fly.

A new species of the genus Xistra (Orthoptera: Tetrigoidea: Metrodorinae) with comments on the characters of mitochondrial genome

A new species from Jiangxi Province of China, Xistra zhengi Deng, sp. nov. is described and illustrated. Meanwhile, the complete mitochondrical genome of the new species was determined and annotated. It has the typical invertebrate mitochondrial gene arrangement, the size of the sequenced mitogenomes is 18148 bp. The constructed phylogenetic tree showed that the new species was placed in Metrodorinae, and formed a separate clade from other species of Metrodorinae.

Molecular analysis of scats revealed diet and prey choice of grey wolves and Eurasian lynx in the contact zone between the Dinaric Mountains and the Alps

A comprehensive understanding of the dietary habits of carnivores is essential to get ecological insights into their role in the ecosystem, potential competition with other carnivorous species, and their effect on prey populations. Genetic analysis of non-invasive samples, such as scats, can supplement behavioural or microscopic diet investigations. The objective of this study was to employ DNA metabarcoding to accurately determine the prey species in grey wolf (Canis lupus) and Eurasian lynx (Lynx lynx) scat samples collected in the Julian Alps and the Dinaric Mountains, Slovenia. The primary prey of wolves were red deer (Cervus elaphus) (detected in 96% scat samples), European roe deer (Capreolus capreolus) (68%), and wild boar (Sus scrofa) (45%). A smaller portion of their diet consisted of mesocarnivores, small mammals, and domestic animals. In contrast, the lynx diet mostly consisted of European roe deer (82%) and red deer (64%). However, small mammals and domestic animals were also present in lynx diet, albeit to a lesser extent. Our findings indicate that the dietary habits of wolves and lynx are influenced by geographical location. Snapshot dietary analyses using metabarcoding are valuable for comprehending the behaviour and ecology of predators, and for devising conservation measures aimed at sustainable management of both their natural habitats and prey populations. However, to gain a more detailed understanding of wolf and lynx dietary habits and ecological impact, it would be essential to conduct long-term genetic monitoring of their diet.

First mitogenomic characterization of Macromotettixoides (Orthoptera, Tetrigidae), with the descriptions of two new species

Classification of species is commonly based on morphological, molecular, and distribution depending on the taxa. Macromotettixoides Zheng, Wei & Jiang, 2005 and Formosatettix Tinkham, 1937 are both wingless types of Tetrigidae with extremely similar morphological characteristics, and in the current taxonomic system they are placed in two different subfamilies, Metrodorinae and Tetriginae, respectively. It is difficult to clearly identify the species of these two genera by morphological characteristics, and molecular data is often needed to assist identification. Here, the complete mitogenomes of two new species were sequenced and assembled, with that of Macromotettixoidesorthomargina. Molecular data of species of Formosatettix were used to test the monophyly of Macromotettixoides and to re-assess the generic characters, and also to test whether Macromotettixoides belongs to the Asian Metrodorinae or Tetriginae. Furthermore, mitochondrial characteristics were analyzed and the phylogeny of the Tetrigidae reconstructed based on mitochondrial protein-coding genes (PCGs). The results indicated that the two new species were clustered with Macromotettixoides rather than Formosatettix, and the anterior margin of the fastigium and pronotum of the two new species usually had the humeral angle different from that of Formosatettix. Therefore, after integrating morphological and molecular data, the two new species were placed in the genus Macromotettixoides, M.maoershanensissp. nov. and M.brachycornasp. nov. Finally, a phylogenetic reconstruction supported Macromotettixoides being assigned to Tetriginae rather than Metrodorinae, in contrast to the previous classification of this genus.

Complete mitochondrial genome of Stethoconus japonicus (Hemiptera: Miridae): Insights into the evolutionary traits within the family Miridae

The mitochondrial (mt) genome sequence of insects possesses numerous evolutionary traits. To better understand the evolution of mt genomes within the family Miridae, the complete mt genome of the predatory Japanese plant bug Stethoconus japonicus Schumacher was sequenced before undertaking a comparative analysis of all reported plant bug mt genomes. The mt genome of S. japonicus is a closed-circular and double-stranded DNA molecule of 16,274 bp (GenBank: MK341530), which consists of 13 protein-coding genes (PCGs), 2 rRNAs, 22 tRNAs and a putative control region (CR). Consistent with other plant bugs, the mt genome of S. japonicus is strongly AT-biased (73.49 %) with A-skewed (0.202) and C-skewed (-0.248). All 13 PCGs initiate translation using ATN codons and TAA served as complete stop codons for eight PCGs, which as incomplete stop codon "T-" for cox1, nad1, nad5-6 and "TA-" for cox2. Regarding other features, all 22 tRNAs could be folded into typical cloverleaf secondary structures. The control region is 1,717 bp and contains a long tandem repeat sequence of a 165 bp unit repeated six times. Similar sequence with variable number of tandemly repeated units from intra-genus CRs is a distinct characteristic of plant bug mt genomes. Phylogenetic relationships of 15 bugs were eventually analyzed based on Maximum likelihood (ML) and Bayesian inference (BI) methods using 17 mt genome sequences. In the phylogenetic trees, species from a same genus or subfamily are clustered into a branch with high supporting values.and the result suggest that Deraeocorinae is more closely related to Mirinae than Bryocorinae. Finally, this study revealed that mutation of tRNA anticodon is a useful phylogenetic marker that could be of significance for studies of evolutionary patterns.

Comparative mitochondrial genome analysis of three leafhopper species of the genus Abrus Dai & Zhang (Hemiptera: Cicadellidae: Deltocephalinae) from China with phylogenetic implication

BACKGROUND

The phylogenetic position and classification of Athysanini are poorly defined, as it includes a large group of polyphyletic genera that have historically been assigned to it mainly because they still exhibit the most typical deltocephaline genitalic and external body characters but lack the distinctive characteristics that other tribes possess. The bamboo-feeding leafhopper genus Abrus belong to the tribe Athysanini of subfamily Deltocephalinae, which currently comprises 19 valid described species, and are limited to the Oriental and Palaearctic regions in China. Although the taxonomy of Abrus are well updated, the references on comparative mitogenomic analyses of Abrus species are only known for a single species. In this study, we sequenced and analyzed the complete mitochondrial genomes (mitogenomes) of Abrus daozhenensis Chen, Yang & Li, 2012 (16,391bp) and A. yunshanensis Chen, Yang & Li, 2012 (15,768bp) (Athysanini), and compared with published mitogenome sequence of A. expansivus Xing & Li, 2014 (15,904bp).

RESULTS

These Abrus species shared highly conserved mitogenomes with similar gene order to that of the putative ancestral insect with 37 typical genes and a non-coding A + T-rich region. The nucleotide composition of these genomes is highly biased toward A + T nucleotides (76.2%, 76.3%, and 74.7%), AT-skews (0.091 to 0.095, and 0.095), negative GC-skews (- 0.138, - 0.161, and - 0.138), and codon usage. All 22 tRNA genes had typical cloverleaf secondary structures, except for trnS1 (AGN) which lacks the dihydrouridine arm, and distinctively trnG in the mitogenome of A. expansivus lacks the TψC arm. Phylogenetic analyses based on 13 PCGs, 2 rRNA genes, and 22 tRNA genes consistently recovered the monophyletic Opsiini, Penthimiini, Selenocephalini, Scaphoideini, and Athysanini (except Watanabella graminea, previously sequenced species as Chlorotettix nigromaculatus) based on limited available mitogenome sequence data of 37 species.

CONCLUSION

At present, Abrus belongs to the tribe Athysanini based on both morphological and molecular datasets, which is strongly supported in present phylogenetic analyses in both BI and ML methods using the six concatenated datasets: amino acid sequences and nucleotides from different combinations of protein-coding genes (PCGs), ribosomal RNA (rRNAs), and transfer RNA (tRNAs). Phylogenetic trees reconstructed herein based on the BI and ML analyses consistently recovered monophylitic Athysanini, except Watanabella graminea (Athysanini) in Opsiini with high support values.

Characterization of Complete Mitochondrial Genome and Phylogenetic Analysis of a Nocturnal Wasps-Provespa barthelemyi (Hymenoptera: Vespidae)

Genus Provespa contains nocturnal wasps mainly found in the southeastern region of Asia. There are no complete genome resources available of this genus, which hinders the study of its phylogenetic evolution and the origin of nocturnal behavior in the Vespidae family. Through high-throughput sequencing, we obtained the mitochondrial genome of Provespa barthelemyi (Hymenoptera: Vespidae), which is 17,721 base pairs in length and contains 13 protein-coding genes (PCGs), 22 tRNAs, and two rRNAs. We identified four gene rearrangement events of P. barthelemyi that frequently occur in the Vespidae family. We used Maximum Likelihood (ML) methodologies to construct a phylogenetic tree based on the sequenced mitochondrial genome and the available data of reported species belonging to Vespinae. Our findings confirmed the monophyly of Vespinae. Our study reports the first complete mitochondrial genome of Provespa and compares its characteristics with other mitochondrial genomes in the family Vespidae. This research should shed light on the phylogenetic relationships and ecological characteristics of the Vespidae family.

Investigating Monophyly of Typhlocybini Based on Complete Mitochondrial Genomes with Characterization and Comparative Analysis of 19 Species (Hemiptera: Cicadellidae: Typhlocybinae)

Tribes of the leafhopper subfamily Typhlocybinae have traditionally been defined based on differences in forewing and hindwing venation. Except for Typhlocybini (sensu lato), the classification of tribes is relatively stable. The monophyly of Typhlocybini needs to be examined, and the relationships among genera within Typhlocybini have not been resolved. Few mitogenome sequences representative of major lineages of Typhlocybini have been available to facilitate a comprehensive phylogenetic analysis of the tribe. In this study, the complete mitogenomes of 19 species of Typhlocybini were sequenced. The gene arrangements of the 19 new mitogenomes are consistent with ancestral insect mitogenomes. Phylogenetic analyses by both maximum-likelihood and Bayesian methods of 67 species of Typhlocybinae suggest that Zyginellini is paraphyletic with respect to Typhlocybini. The phylogenetic relationships within Typhlocybini are discussed, and the major results show that the Farynala and Linnavuoriana complexes previously recognized based on morphological characters correspond to monophyletic lineages.

The complete mitochondrial genome and gene rearrangements in a gall wasp species, Dryocosmus liui (Hymenoptera: Cynipoidea: Cynipidae)

Mitochondrial genomes (mitogenomes) have been widely used in comparative and evolutionary genomics, molecular evolution, phylogenetics, and population genetics, but very limited information is available for the family Cynipidae. In this report, we describe the mitogenome of Dryocosmus liui Pang, Su et Zhu, providing the first complete mitogenomic data for a cynipid gall wasp species. The mitogenome of D. liui is 16,819 bp in length, and contains the typical set of 37 genes. Two control regions were detected, with the second being a perfect inverted repeat of the major portion of the first. Gene rearrangements were found in transfer RNA (tRNA) genes, protein-coding genes (PCGs) and ribosomal RNA (rRNA) genes, compared with the putative ancestral mitogenome. Similar to two other Cynipidae species with mitogenome data available, D. liui has a novel tRNA gene cluster trnL1-trnI-trnL2-trnW-trnM-trnQ between nad1 and nad2. Phylogenetic analysis based on sequences of PCGs and rRNA genes with D. liui included obtained topologies identical to previous studies supporting the a relationship of (Cynipoidea , (Platygastroidea, Proctotrupoidea)) within the monophyletic Proctotrupomorpha and (Cynipidae, Figitidae), Ibaliidae) within the Cynipoidea.

The First Complete Mitochondrial Genome of the Genus Pachycondyla (Formicidae, Ponerinae) and Insights into the Phylogeny of Ants

Ants are the standout group among eusocial insects in terms of their exceptional species richness and ecological dominance. The phylogenetic relationships among the group remain elusive. Mitochondrial genome sequences, as a kind of molecular marker, have been widely utilized in the phylogenetic analysis of insects. However, the number of ant mitogenomes published is still very limited. In this study, we utilized next-generation sequencing to determine the complete mitogenome of Pachycondyla annamita (Formicidae, Ponerinae). This is the first mitogenome from the genus Pachycondyla. Two gene rearrangements were identified in the mitogenome, the transposition of trnQ and trnM and the transposition of trnV and rrnS. The secondary structures of tRNAs were predicted. The tRNA genes trnR and trnS1 lacked the dihydrouridine (DHU) arm, and the trnE lacked the TΨC (T) arm. Phylogenetic analyses of the mitochondrial protein-coding genes under maximum likelihood (ML) and Bayesian inference (BI) criteria resulted in conflicting hypotheses. BI analysis using amino acid data with the site-heterogeneous mixture model produced a tree topology congruent with previous studies. The Formicidae was subdivided into two main clades, namely the "poneroid" clade and the "formicoid" clade. A sister group relationship between Myrmicinae and Formicinae was recovered within the "formicoid" clade.

Molecular Phylogeny and Evolution of the Tuerkayana (Decapoda: Brachyura: Gecarcinidae) Genus Based on Whole Mitochondrial Genome Sequences

Tuerkayana is of particular interest because it has been separated, in recent years, from Cardisoma and Discoplax but studies of its taxonomic status, especially from a whole mitochondrial genome perspective, have been lacking. In this study, the mitogenomes of four species (Tuerkayana magnum, Tuerkayana rotundum, Tuerkayana hirtipes, and Tuerkayana celeste) of Tuerkayana are sequenced and contrasted with other species in Brachyura for the first time. The phylogenetic tree of Brachyura, which includes 206 crab species (189 species of Brachyuran and 17 Anomura species) with a complete mitogenome, was constructed to evaluate the phylogenetic position of Tuerkayana and Gecarcinidae within Brachyuran, and explore the monophyly of Gecarcinidae. Furthermore, two single gene trees based on cox1 and 16SrRNA separately within interspecies of Gecarcinidae were reconstructed, providing molecular evidence for Tuerkayana and further clarifying the division of genera in Gecarcinidae. Based on the mitogenome dataset of 206 crabs, the branch-site model was utilized to explore selective pressure in individual codons with CodeML. The strong selective pressure shown in nad6 indicates that it may have played a significant role in the evolution of Gecarcinidae.

First Description of the Nuclear and Mitochondrial Genomes and Associated Host Preference of Trichopoda pennipes, a Parasitoid of Nezara viridula

Trichopoda pennipes is a tachinid parasitoid of several significant heteropteran agricultural pests, including the southern green stink bug, Nezara viridula, and leaf-footed bug, Leptoglossus phyllopus. To be used successfully as a biological control agent, the fly must selectively parasitize the target host species. Differences in the host preference of T. pennipes were assessed by assembling the nuclear and mitochondrial genomes of 38 flies reared from field-collected N. viridula and L. phyllopus. High-quality de novo draft genomes of T. pennipes were assembled using long-read sequencing. The assembly totaled 672 MB distributed among 561 contigs, having an N50 of 11.9 MB and a GC of 31.7%, with the longest contig at 28 MB. The genome was assessed for completeness using BUSCO in the Insecta dataset, resulting in a score of 99.4%, and 97.4% of the genes were single copy-loci. The mitochondrial genomes of the 38 T. pennipes flies were sequenced and compared to identify possible host-determined sibling species. The assembled circular genomes ranged from 15,345 bp to 16,390 bp and encode 22 tRNAs, two rRNAs, and 13 protein-coding genes (PCGs). There were no differences in the architecture of these genomes. Phylogenetic analyses using sequence information from 13 PCGs and the two rRNAs individually or as a combined dataset resolved the parasitoids into two distinct lineages: T. pennipes that parasitized both N. viridula and L. phyllopus, and others that parasitized only L. phyllopus.

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.

Phylogenomics including the newly sequenced mitogenomes of two moths (Noctuoidea, Erebidae) reveals Ischyja manlia (incertae sedis) as a member of subfamily Erebinae

We sequenced the mitogenomes of two Erebid species, namely Ischyja manlia (Cramer, 1776) and Rusicada privata (Walker, 1865) to analyse the phylogenetic relationship and to establish the taxonomic position of incertae sedis members of the family Erebidae. The two circular genomes of I. manlia and R. privata were 15,879 bp and 15,563 bp long, respectively. The gene order was identical, containing 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and an A + T-rich region. The nucleotide compositions of the A + T-rich region of both mitogenomes were similar: 80.65% for R. privata, and 81.09% for I. manlia. The AT skew and GC skew were slightly positive in I. manlia and negative in R. privata. In I. manlia and R. privata, except for cox1 which started with CGA and TTG codons, all the other 12 PCGs started with ATN codon. The A + T-rich regions of I. manlia and R. privata were 433 and 476 bp long, respectively, and contained common characteristics of Noctuoidea moths. At present, Ischyja is treated as Erebinae incertae sedis. However, phylogenetic analysis conducted in the present study reveals that the genus Ischyja is most likely to be a member of the subfamily Erebinae.

The Complete Mitochondrial Genome of Box Tree Moth Cydalima perspectalis and Insights into Phylogenetics in Pyraloidea

To resolve and reconstruct phylogenetic relationships within Pyraloidea based on molecular data, the mitochondrial genome (mitogenome) was widely applied to understand phylogenetic relations at different taxonomic levels. In this research, a complete mitogenome of Cydalima perspectalis was recorded, and the phylogenetic position of C. perspectalis was inferred based on the sequence in combination with other available sequence data. According to the research, the circular mitochondrial genome is 15,180 bp in length. It contains 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), 13 typical protein-coding genes (PCGs), and a non-coding control region. The arrangement of a gene of the C. perspectalis mitogenome is not the same as the putative ancestral arthropod mitogenome. All of the PCGs are initiated by ATN codons, except for the cytochrome c oxidase subunit 1 (cox1) gene, which is undertaken by CGA. Five genes have incomplete stop codons that contain only ‘T’. All tRNA genes display a typical clover–leaf structure of mitochondrial tRNA, except for trnS1 (AGN). The control region contained an ‘ATAGG(A)’-like motif followed by a poly-T stretch. Based on the mitochondrial data, phylogenetic analysis within Pyraloidea was carried out using Bayesian inference (BI) and maximum likelihood (ML) analyses. Phylogenetic analysis showed that C. perspectalis is more closely related to Pygospila tyres within Spilomelinae than those of Crambidae and Pyraloidea.

Potato (Solanum tuberosum) as a New Host for Pentastiridius leporinus (Hemiptera: Cixiidae) and Candidatus Arsenophonus Phytopathogenicus

Pentastiridius leporinus is a planthopper (Hemiptera: Cixiidae) that vectors two phloem-restricted bacterial pathogens to sugar beet (Beta vulgaris (L.)): the γ-proteobacterium Candidatus Arsenophonus phytopathogenicus and the stolbur phytoplasma Candidatus Phytoplasma solani. These bacteria cause an economically important disease known as syndrome basses richesses (SBR), characterized by yellowing, deformed leaves and low beet yields. Having observed potato fields in Germany infested with cixiid planthoppers and showing signs of leaf yellowing, we used morphological criteria and COI and COII as molecular markers, to identify the planthoppers (adults and nymphs) primarily as P. leporinus. We analyzed planthoppers, potato tubers, and sugar beet roots and detected both pathogens in all sample types, confirming that P. leporinus adults and nymphs can transmit the bacteria. This is the first time that P. leporinus has been shown to transmit Arsenophonus to potato plants. We also found that two generations of P. leporinus were produced in the warm summer of 2022, which will probably increase the pest population size (and thus the prevalence of SBR) in 2023. We conclude that P. leporinus has expanded its host range to potato, and can now utilize both host plants during its developmental cycle, a finding that will facilitate the development of more efficient control strategies.

Phylogenetic Implications of Mitogenomic Sequences and Gene Rearrangements of Scale Insects (Hemiptera, Coccoidea)

Coccoidea (scale insects) are important plant parasites with high diversity of species. However, the phylogenetic relationship within Coccoidea has not been fully determined. In this study, we sequenced mitogenomes of six species belonging to five coccoid families. With the addition of three previously published mitogenomes, a total of 12 coccoid species were adopted for the phylogenetic reconstruction based on the maximum likelihood and Bayesian inference. The monophyly of Coccoidea was recovered and Aclerdidae and Coccidae were recovered as the sister group, successively sister to Cerococcidae, Kerriidae, and Eriococcidae. In addition, there were gene rearrangements occurring in all mitogenomes of coccoid species studied here. The novel gene rearrangement ND6-trnP and trnI-ND2-trnY supported the monophyly of Coccoidea and the sister relationship of Aclerdidae and Coccidae. This implies that data from the mitogenome can provide new insight for clarifying the deeper level of phylogenetic relationship within Coccoidea.

A review of Eupholidoptera (Orthoptera, Tettigoniidae) from Crete, Gavdos, Gavdopoula, and Andikithira

Being nocturnal, hiding in prickly bushes and shrubs during the day, Eupholidoptera species in Crete and its neighbouring islands are easily overlooked, and until now our knowledge about their distribution was based on some thirty sightings across 11 species. In this paper results are presented of a study of Eupholidoptera specimens collected between 1987 and 2020 by hand-catches and pitfall and fermenting traps on the Greek islands of Crete, Gavdos, Gavdopoula, and Andikithira. Diagnostic features of all known species are presented and illustrated with stacked images. An updated key to all species is provided. Eupholidoptera francisae Tilmans & Odé, sp. nov. from Andikithira and southwestern Crete and Eupholidoptera marietheresae Willemse & Kotitsa, sp. nov. from Mt. Dikti are described. Female E. cretica, E. gemellata, and E. mariannae are described, and the female of E. astyla is redescribed. Bioacoustics for E. francisae Tilmans & Odé, sp. nov., E. giuliae, and E. jacquelinae are presented for the first time. Eupholidoptera smyrnensis is reported for the first time from Crete. A substantial amount of new distribution data for Eupholidoptera species on Crete is presented. The current distribution pattern and first analyses of phylogeny based on molecular data of Eupholidoptera species on Crete are discussed in relation to paleogeographical events.

Improvement in the Identification Technology for Asian Spongy Moth, Lymantria dispar Linnaeus, 1758 (Lepidoptera: Erebidae) Based on SS-COI

Lymantria dispar (Linnaeus, 1758), which is commonly known as spongy moth, with two subspecies, is found in Asia: Lymantria dispar asiatica and Lymantria dispar japonica, collectively referred to as the Asian spongy moth (ASM). The subspecies Lymantria dispar dispar occurs in Europe and is commonly known as the European spongy moth (ESM). The ASM is on the quarantine list of many countries because it induces greater economic losses than the ESM. Accurate identification is essential to prevent the invasion of ASM into new areas. Although several techniques for identifying ASMs have been developed, the recent discovery of complex patterns of genetic variation among ASMs in China as well as new subspecies in some areas has necessitated the development of new, improved identification techniques, as previously developed techniques are unable to accurately identify ASMs from all regions in China. Here, we demonstrate the efficacy of an improved technique for the identification of the ASM using ASM-specific primers, which were designed based on cytochrome oxidase I sequences from samples obtained from all sites where ASMs have been documented to occur in China. We show that these primers are effective for identifying a single ASM at all life stages and from all ASM populations in China, and the minimum detectable concentration of genomic DNA was 30 pg. The inclusion of other Lymantria samples in our analysis confirmed the high specificity of the primers. Our improved technique allows the spread of ASMs to be monitored in real time and will help mitigate the spread of ASMs to other areas.

The complete mitochondrial genome of Aeschrocoristuberculatus and A.ceylonicus (Hemiptera, Pentatomidae) and its phylogenetic implications

Aeschrocoristuberculatus and A.ceylonicus (Hemiptera, Pentatomidae, Pentatominae) are mainly distributed in southern China, India, Myanmar, and Sri Lanka. Both species are also common agricultural pests. However, only the morphology of the genus Aeschrocoris has previously been studied, and molecular data have been lacking. In this study, the whole mitochondrial genomes of A.tuberculatus and A.ceylonicus are and annotated. The lengths of the complete mitochondrial genomes of the two species are 16,134 bp and 16,142 bp, respectively, and both contain 37 typical genes, including 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and a control region. The mitochondrial genome structure, gene order, nucleotide composition, and codon usage of A.tuberculatus and A.ceylonicus are consistent with those of typical Pentatomidae. Most PCGs of both species use ATN as the start codon, except atp8, nad1, and cox1, which use TTG as the start codon. cox1, cox2, and atp6 use a single T, and nad1 use TAG as the stop codon; the remaining PCGs have TAA as the stop codon. The A+T contents of the two species are 73.86% and 74.08%, respectively. All tRNAs have a typical cloverleaf structure, with the exception of trnS1, which lacks a dihydrouridine arm. The phylogenetic tree is reconstructed using the maximum-likelihood method based on the newly obtained mitochondrial genome sequences and 87 existing mitochondrial genomes of Pentatomoidea from the NCBI database and two species of Lygaeoidea as outgroups. The phylogenetic trees strongly support the following relationships: (Urostylididae + ((Acanthosomatidae + ((Cydnidae + (Dinidoridae + Tessaratomidae)) + (Scutelleridae + Plataspidae))) + Pentatomidae). This study enriches the mitochondrial genome database of Pentatomoidea and provides a reference for further phylogenetic studies.

Comparative analysis and phylogeny of mitochondrial genomes of Pentatomidae (Hemiptera: Pentatomoidea)

The Phyllocephalini is a group of herbivorous insects in Pentatomidae, which lack distinctive morphological characteristics and systematic studies. Up to now, there are only two complete mitochondrial genomes of Phyllocephalini have been reported. In this study, we sequenced and analyzed the complete mitochondrial genomes of three Phyllocephalini species, Gonopsis coccinea, Gonopsimorpha nigrosignata, and Chalcopis glandulosus, which were 16,534, 16,531, and 16,534 bp in length, respectively. The mitochondrial genomes contained 37 genes, including 13 protein-coding genes, two rRNA genes, 22 tRNA genes, and a control region. The gene arrangement was consistent with that of the putative ancestral insect, with no rearrangement. The cox1 gene of Pentatomidae showed the lowest evolutionary rate among the protein-coding genes, the mean genetic distance of species, genera, and subfamilies of Pentatomidae increased hierarchically based on cox1 gene. The 16S rRNA of Pentatomidae was more conserved than 12S rRNA in sequence and secondary structure. All tRNAs could be folded into a typical cloverleaf structure except trnS1. The stem region was more conserved than the loop region in the secondary structure of tRNAs within Pentatomidae. Gonopsis coccinea and Gonopsimorpha nigrosignata had one type of tandem repetition unit in the control region, while C. glandulosus had two types. The heterogeneity analysis of Pentatomidae showed that Phyllocephalinae was the most heterogeneous. Phylogenetic trees based on the newly obtain mitochondrial genomes along with other 50 mitochondrial genomes of Pentatomidae using Bayesian Inference and Maximum Likelihood strongly supported the following three relationships: (((Anaxilaus + (Plautia + Glaucias)) + (Nezara + Palomena)) + (Eysarcorini + Carpocorini)), (Hoplistoderini + (Menidini + Asopinae)), and ((Sephelini + Halyini) + (Caystrini + (Cappaeini + (Placosternum + Phyllocephalini)))). The relationships within Phyllocephalini were (Chalcopis + (Dalsira + (Gonopsimorpha + Gonopsis))). Our results provide valuable molecular data for further phylogenetic analyses of Pentatomidae.

Population genetic structure of Neoschongastia gallinarum in South China based on mitochondrial DNA markers

The genetic diversity and differentiation of four geographic populations of Neoschongastia gallinarum were evaluated using concatenated mitochondrial gene sequences (pCOI, pCOII, and pND5). Based on the results, the N. gallinarum populations had high genetic diversity and strong ecological adaptability. Genetic differentiation among paired populations calculated using concatenated mitochondrial gene sequences revealed that geographic isolation resulted in genetic differentiation among the populations of N. gallinarum, and gene flow between populations associated with human trade activities. Systematic development and molecular variance based on haplotypes revealed that genetic variation existed in different haplotypes; however, no clear rule related to geographic region was found. Further, genetic variation was mainly derived from individuals within the population. A neutral test based on concatenated mitochondrial gene sequences and nucleotide pair differences revealed that N. gallinarum did not experience an obvious population expansion in recent historical periods. Accordingly, the population size was relatively stable.

Characterization of the mitochondrial genomes of Bradysia hygida, Phytosciara flavipes and Trichosia splendens (Diptera: Sciaridae) and novel insights on the control region of sciarid mitogenomes

Sciarids, also called "fungus gnats" are small, almost entirely dark-coloured insects. Sciarid larvae feed on different substrates and can infest agricultural crops and mushroom nurseries, causing economic losses. Of the 2174 Diptera mitogenome sequences currently available in GenBank, only eight are from the Sciaridae family, none of which are complete circular molecules. Here we describe the mitogenome sequences of three sciarid species: Phytosciara flavipes, Trichosia splendens and Bradysia hygida and provide novel insights on the control region of sciarid mitogenomes. The assembled mitogenomes range from 16,062 bp in P. flavipes to 17,095 bp in B. hygida. All 13 protein coding genes, 22 tRNAs and 2 rRNAs characteristic of insect mitogenomes were identified, but the sequence of the control region could not be determined. Experimental results suggest that the B. hygida control region is about 21 kb long resulting in a 37 kb long mitogenome which constitutes the largest insect mitochondrial genome described so far. Phylogenetic analysis using all Bibionomorpha mitogenome sequences available in GenBank strongly supports the Sciaridae monophyly and led to the identification of species and subfamily specific gene rearrangements. Our study extends the knowledge of this large and diverse insect family that includes agricultural pest species.

Two complete mitochondrial genomes in Scolopendra and a comparative analysis of tRNA rearrangements in centipedes

BACKGROUND

Centipedes are one of the oldest terrestrial arthropods belonging to the sub phylum Myriapoda. With the expansion of our understanding of the application of the two centipedes Scolopendra morsitans and Scolopendra hainanum, belonging to the order Scolopendromorpha, an exhaustive classification was required. Although consensus has been reached on the phylogeny of Chilopoda based on morphological traits, recent analyses based on molecular data exhibited differences in results.

METHODS AND RESULTS

The mitochondrial genome sequences of S. morsitans and S. hainanum were obtained by next-generation sequencing. S. morsitans contains 13 PCGs, two rRNAs, 11 tRNAs, and one CR. whereas S. hainanum contains 12 PCGs, of which ATP8 remains unpredicted, two rRNAs, 14 tRNAs, and one CR. An obvious tRNA rearrangement was found in the genus Scolopendra. S. morsitans exhibited a loss of trnW, trnC, trnI, trnK, trnD, trnA, trnN, trnQ, trnF, trnT, trnS, trnL, and trnV, and a repeat of trnR and trnL. S. hainanum exhibited a loss of trnQ, trnC, trnW, trnI, trnD, trnQ, trnP, and trnV. Phylogenetic analyses of centipedes based on 12 PCGs supported the sister relationship between the orders Geophilomorpha and Lithobiomorpha and a close relationship between Scolopendra dehaani and S. hainanum.

CONCLUSIONS

The new mitogenomes determined in this study provide new genomic resources for gene rearrangements and contribute to the understanding of the evolution of gene rearrangement in Chilopoda.

Novel Mitochondrial Gene Rearrangement and Intergenic Regions Exist in the Mitochondrial Genomes from Four Newly Established Families of Praying Mantises (Insecta: Mantodea)

Simple Summary

Mantodea is regarded as an excellent material to study the gene rearrangements and large non-coding regions (LNCRs) in mitochondrial genomes. Meanwhile, as a result of the convergent evolution and parallelism, the gene rearrangements and LNCRs are specific to some taxonomic groups within Mantodea, which play an important role in phylogenetic relationship research. Nine mitochondrial genomes (mitogenomes) from four newly established families of praying mantises are obtained and annotated. Eight types of gene rearrangements, including four novel types of gene rearrangements in Mantodea, are detected, which can be explained by the tandem replication-random loss (TDRL) model. Moreover, one conserved motif between trnI-trnQ is detected in Toxoderidae. This study shed light on the formation mechanisms of these gene rearrangements and LNCRs in four newly established families of praying mantises.

Abstract

Long non-coding regions (NCRs) and gene rearrangements are commonly seen in mitochondrial genomes of Mantodea and are primarily focused on three regions: CR-I-Q-M-ND2, COX2-K-D-ATP8, and ND3-A-R-N-S-E-F-ND5. In this study, eight complete and one nearly complete mitochondrial genomes of praying mantises were acquired for the purpose of discussing mitochondrial gene rearrangements and phylogenetic relationships within Mantodea, primarily in the newly established families Haaniidae and Gonypetidae. Except for Heterochaeta sp. JZ-2017, novel mitochondrial gene arrangements were detected in Cheddikulama straminea, Sinomiopteryx graham, Pseudovates chlorophaea, Spilomantis occipitalis. Of note is the fact that one type of novel arrangement was detected for the first time in the Cyt b-S2-ND1 region. This could be reliably explained by the tandem replication-random loss (TDRL) model. The long NCR between trnT and trnP was generally found in Iridopteryginae and was similar to the ND4L or ND6 gene. Combined with gene rearrangements and intergenic regions, the monophyly of Haaniidae was supported, whereas the paraphyly of Gonypetidae was recovered. Furthermore, several synapomorphies unique to some clades were detected that conserved block sequences between trnI and trnQ and gaps between trnT and trnP in Toxoderidae and Iridopteryginae, respectively.

A LAMP Assay for the Detection of Thecodiplosis japonensis, an Alien Gall Midge Species Pest of Pine Trees

Simple Summary

Thecodiplosis japonensis is considered the most harmful pest to pines in South Korea. T. japonensis is a native species of Japan. Recently, T. japonensis was discovered in China and has caused serious damage to local pine trees. Due to the small size and little morphological difference with its related species, it is difficult to accurately identify T. japonensis by morphological methods. Accurate and efficient molecular identification methods are urgently needed to detect this invasive gall midge pest, yet there was no molecular identification method for T. japonensis. In this study, we developed a LAMP assay to detect T. japonensis based on the COI gene sequence. The LAMP assay could detect as little as 300 fg of gDNA. Using colorimetric amplification and a crude gDNA extraction method, the total procedure could be processed in 75 min. The method established in the study can be easily used in both laboratory and field conditions, enabling rapid molecular identification of T. japonensis.

Abstract

Pine needle gall midge (T. japonensis), native to Japan, has become a serious invasive pest in South Korea and, more recently in 2006, in China. It was first discovered in Qingdao, Shandong Province, and has caused serious damage to local Pinus thunbergii. The insect’s small size makes morphological-based identification difficult; therefore, molecular detection techniques are urgently needed for monitoring and preventing its further spread. At present, there is no simple and accurate field molecular identification tool. To solve this problem, a LAMP-based molecular diagnosis technology of T. japonensis was developed. Four LAMP primers were designed to specifically amplify T. japonensis DNA. Positive LAMP reactions usually produce amplification in one hour. The optimal incubation conditions for LAMP detection were determined with 4 LAMP primers for 60 min at 61 °C. The LAMP detection range of gDNA concentrations is wide, with a minimum detectable gDNA concentration of 300 fg. A non-destructive DNA-releasing procedure, HotSHOT “HS6”, which could extract “crude DNA” for LAMP assay in 10 min, was used for larval and adult samples. Therefore, we established a LAMP-based rapid molecular identification method that can be applied in the monitoring and management of T. japonensis.

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.

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

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.

Universal Mitochondrial Multi-Locus Sequence Analysis (mtMLSA) to Characterise Populations of Unanticipated Plant Pest Biosecurity Detections

Simple Summary

Agricultural and environmental sustainability requires effective biosecurity responses that prevent the establishment or spread of exotic insect pests. Understanding where new detections may have come from or if recurrent detections are connected contributes to this. Suitable population genetic markers use relatively rapidly evolving gene regions which render the PCR method species-specific at best. Because resource limitations mean these are pre-emptively developed for the highest risk species, populations of other exotic pests are unable to be characterised at the time. Here we have developed a generic method that is useful across species within the same taxonomic Order, including where there is little or no prior knowledge of their gene sequences. Markers are formed by concomitant sequencing of four gene regions. Sequence concatenation was shown to retrieve higher resolution signatures than standard DNA barcoding. The method is encouragingly universal, as illustrated across species in ten fly and 11 moth superfamilies. Although as-yet untested in a biosecurity situation, this relatively low-tech, off-the-shelf method makes a proactive contribution to the toolbox of quarantine agencies at the time of detection without the need for impromptu species-specific research and development.

Abstract

Biosecurity responses to post-border exotic pest detections are more effective with knowledge of where the species may have originated from or if recurrent detections are connected. Population genetic markers for this are typically species-specific and not available in advance for any but the highest risk species, leaving other less anticipated species difficult to assess at the time. Here, new degenerate PCR primer sets are designed for within the Lepidoptera and Diptera for the 3′ COI, ND3, ND6, and 3′ plus 5′ 16S gene regions. These are shown to be universal at the ordinal level amongst species of 14 and 15 families across 10 and 11 dipteran and lepidopteran superfamilies, respectively. Sequencing the ND3 amplicons as an example of all the loci confirmed detection of population-level variation. This supported finding multiple population haplotypes from the publicly available sequences. Concatenation of the sequences also confirmed that higher population resolution is achieved than for the individual genes. Although as-yet untested in a biosecurity situation, this method is a relatively simple, off-the-shelf means to characterise populations. This makes a proactive contribution to the toolbox of quarantine agencies at the time of detection without the need for unprepared species-specific research and development.

Macrostructural Evolution of the Mitogenome of Butterflies (Lepidoptera, Papilionoidea)

Simple Summary

Papilionoidea is a superfamily of Lepidoptera encompassing about 19,000 species. In the present work, we study the evolution of the structure of the mitogenome of these lepidopterans. The mechanisms generating the eight arrangements known for Papilionoidea were investigated analysing the movements of different mitochondrial genes. Five newly sequenced/assembled mitogenomes were included in our analysis involving more than 600 genomes. We provide new findings that help to understand the evolution of the gene orders MIQGO, IMQGO, 2S1GO, ES1GO and S1NGO in different butterflies. We demonstrate that the evolution of the 2S1GO in Lycaenidae followed a complicated pathway with multiple events of duplication and loss of trnS1 and changes in anticodon. We describe two new gene orders 2FFGO and 4QGO for Ampittia subvittatus (Hesperiidae) and Bhutanitis thaidina (Papilionidae).

Abstract

The mitogenome of the species belonging to the Papilionodea (Lepidoptera) is a double stranded circular molecule containing the 37 genes shared by Metazoa. Eight mitochondrial gene orders are known in the Papilionoidea. MIQGO is the plesiomorphic gene order for this superfamily, while other mitochondrial arrangements have a very limited distribution. 2S1GO gene order is an exception and is present in several Lycaenidae and one species of Hesperiidae. We studied the macrostructural changes generating the gene orders of butterflies by analysing a large data set (611 taxa) containing 5 new mitochondrial sequences/assemblies and 87 de novo annotated mitogenomes. Our analysis supports a possible origin of the intergenic spacer trnQ-nad2, characterising MIQGO, from trnM. We showed that the homoplasious gene order IMQGO, shared by butterflies, species of ants, beetles and aphids, evolved through different transformational pathways. We identify a complicated evolutionary scenario for 2S1GO in Lycaenidae, characterised by multiple events of duplication/loss and change in anticodon of trnS1. We show that the gene orders ES1GO and S1NGO originated through a tandem duplication random loss mechanism. We describe two novel gene orders. Ampittia subvittatus (Hesperiidae) exhibits the gene order 2FFGO, characterised by two copies of trnF, one located in the canonical position and a second placed in the opposite strand between trnR and trnN. Bhutanitis thaidina (Papilionidae) exhibits the gene order 4QGO, characterised by the quadruplication of trnQ.

Characterization of the complete mitogenome of brown-spotted locust, Cyrtacanthacris tatarica (Orthoptera: Acrididae: Cyrtacanthacridinae)

The complete mitochondrial genome of Cyrtacanthacris tatarica was firstly sequenced and analyzed. The circular mitogenome was 15,679 bp long, showing a bias of AT rich on the majority strand (42.34% of A, 29.99% of T, 11.19% of G, and 16.18% of C). It consisted of the typical 37 genes (13 protein-coding genes, two ribosomal RNA genes, and 22 transfer tRNA genes) and one longest non-coding region (called as control region). All PCGs used standard ATN initiation codons, and most PCGs were terminated with complete codons (TAA/TAG), apart from cox1 and nad5. Phylogenetic analyses based on the concatenated nucleotide sequences of PCGs supported that Cyrtacanthacridinae was monophyletic, and the sister group relationship between C. tatarica and Schistocerca gregaria gregaria was determined. Our results may provide molecular information for the genetic evolution and taxonomy of the acridid species.