Quantification of rearrangements and evolution of mitochondrial gene order of Acari (Chelicerata: Arachnida)

Invertebrate mitogenomes are generally fixed with formal 37 genes: 13 PCGs encoded subunits of OXPHOS, 2 ribosomal RNA (rRNA) functional in the translation of these PCGs and 22 transfer RNA (tRNA) genes. The order of these genes varies greatly among organisms and named rearrangement. Rearrangement patterns of mitochondrial genomes may shed light on mutation processes and evolutionary relationships of organisms. Mitochondrial gene organization is highly variable among Acari, so rearrangement is a very common mitogenomic pattern in this group. In this study, 258 unique Acari (Acariformes + Parasitiformes) mitogenomes were downloaded from NCBI and studied about rearrangement patterns. Sixty-seven mitotypes were determined among Acari and the most rearranged genes were trnL1 and nad2. Following that, trnI, trnS1, trnN, trnE, trnT, and trnP genes are remarkably mobile (RF > 95%). Conversely, atp6, cox3, trnG, and cytb genes also appears to be quite stable (RF < 20%). Within Acari, mean distance calculations are varied from 1.210 in atp8 to 0.155 in rrnS. Contrary to expectations, among Acari mobile tRNA genes appear to be conserved in sequences, whereas PCGs have higher distance values and seem to be mutated. Consistently, tRNA genes seem saturated, but some PCGs (atp6, cox genes, cytb, nad1, and nad6) are not saturated. These values do not correlate with each other (p > 0.005). This discrepancy may indicate that the genes were rearranged after mutation load; consistent with this, DAMBE saturation values are also not correlated with RF values. Parasitiformes mitogenomes are more mobile than Acariformes mitogenomes and may be under the effect of selective sweeping.

qGO: a novel method for quantifying the diversity of mitochondrial genome organization

Quantifying the features of mitochondrial genome structural variation is crucial for understanding its contribution to complexity. Accurate quantification and interpretation of organizational diversity can help uncover biological evolutionary laws and patterns. The current qMGR approach accumulates the changes in two adjacent genes to calculate the rearrangement frequency RF of each single gene and the rearrangement score RS for specific taxa in the mitogenomes of a given taxonomic group. However, it may introduce bias, as it assigns scores to adjacent genes rather than to rearranged genes. To overcome this limitation, we propose a novel statistical method called qGO to quantify the diversity of gene organization. The qGO method, which is based on the homology of gene order, provides a more accurate representation of genome organizational diversity by partitioning gene strings and individually assigning weights to genes spanning different regions. Additionally, a comprehensive approach is employed for distance computation, generating an extensive matrix of rearrangement distances. Through experiments on more than 5500 vertebrate mitochondrial genomes, we demonstrated that the qGO method outperforms existing methods in terms of accuracy and interpretability. This method improves the comparability of genomes and allows a more accurate comparison of the diversity of mitochondrial genome organization across taxa. These findings have significant implications for unraveling genome evolution, exploring genome function, and investigating the process of molecular evolution.

Comparative mitogenomic analysis and phylogeny of Veneridae with doubly uniparental inheritance

Doubly uniparental inheritance (DUI) is an atypical animal mtDNA inheritance system, reported so far only in bivalve species, in which two mitochondrial lineages exist: one transmitted through the egg (F-type) and the other through the sperm (M-type). Although numerous species exhibit this unusual organelle inheritance, it is primarily documented in marine and freshwater mussels. The distribution, function and molecular evolutionary implications of DUI in the family Veneridae, however, remain unclear. Here, we investigated 17 species of Veneridae, compared mitochondrial genomes of DUI species and reconstructed their phylogenetic framework. Different sex-linked mitochondrial genomes have been identified in the male gonads and adductor muscles of 7 venerids, indicating the presence of DUI in these species. Analysis of the unassigned regions (URs) of the mitochondrial genome in DUI species revealed that 13 out of 44 URs contained repetitive sequences, with nine being long unassigned regions (LURs). All LURs were capable of forming secondary structures, and most of them exhibited patterns of significant sequence similarity to elements known to have specific functions in the control regions of sea urchins and mammals. The F/M phylogeny showed that DUI venerids exhibit both taxon-specific patterns and gender-specific patterns, with Gafrarium dispar experiencing masculinization events.

The first complete mitochondrial genome of the genus Laelaps with novel gene arrangement reveals extensive rearrangement and phylogenetics in the superfamily Dermanyssoidea

We collected 56 specimens of Laelaps chini from the endemic Hengduan Mountain rat species (Eothenomys miletus) and obtained the first complete mitochondrial genome of L. chini by next-generation sequencing (NGS). The L. chini mitogenome is 16,507 bp in size and contains 37 genes and a control region of 2380 bp in length. The L. chini mitogenome has a high AT content and a compact arrangement with four overlapping regions ranging from 1 to 2 bp and 16 spacer regions ranging from 1 to 48 bp. We analyzed 13 protein-coding genes of L. chini mitogenome  and found that protein-coding genes in the L. chini mitogenome preferred codons ending in A/U and codon usage pattern was mainly influenced by natural selection. Cox1 has the slowest evolution rate and cox3 has the fastest evolution rate. We combined the mitochondrial genome of eight species of gamasid mites in the superfamily Dermanyssoidea from Genbank and the L. chini mitochondrial genome to analyze its rearrangement patterns and breakpoint numbers. We found that the L. chini mitogenome showed a novel arrangement pattern and nine species of gamasid mites in the superfamily Dermanyssoidea, which have been sequenced complete mitochondrial genomes to date, all showed different degrees of rearrangement. Laelaps chini, Echinolaelaps echidninus and Echinolaelaps fukinenensis were closely related species based on genetic distance and phylogenetic analyses. Notably they are clustered with Varroa destructor of the family Varroidae, suggesting that the family Varroidae is more closely related to the family Laelapidae, but more data are needed to test whether Varroa can be classified under the family Laelapidae. The L. chini mitogenome is the first complete mitochondrial genome for the genus Laelaps, and contributes to further exploration of the mitochondrial gene rearrangements and phylogeny for the superfamily Dermanyssoidea.

Full-Length Transcriptome Profiling of the Complete Mitochondrial Genome of Sericothrips houjii (Thysanoptera: Thripidae: Sericothripinae) Featuring Extensive Gene Rearrangement and Duplicated Control Regions

The mitochondrial genome (mitogenome) of Thysanoptera has extensive gene rearrangement, and some species have repeatable control regions. To investigate the characteristics of the gene expression, transcription and post-transcriptional processes in such extensively gene-rearranged mitogenomes, we sequenced the mitogenome and mitochondrial transcriptome of Sericothrips houjii to analyze. The mitogenome was 14,965 bp in length and included two CRs contains 140 bp repeats between COIII-trnN (CR1) and trnT-trnP (CR2). Unlike the putative ancestral arrangement of insects, S. houjii exhibited only six conserved gene blocks encompassing 14 genes (trnL2-COII, trnD-trnK, ND2-trnW, ATP8-ATP6, ND5-trnH-ND4-ND4L and trnV-lrRNA). A quantitative transcription map showed the gene with the highest relative expression in the mitogenome was ND4-ND4L. Based on analyses of polycistronic transcripts, non-coding RNAs (ncRNAs) and antisense transcripts, we proposed a transcriptional model of this mitogenome. Both CRs contained the transcription initiation sites (TISs) and transcription termination sites (TTSs) of both strands, and an additional TIS for the majority strand (J-strand) was found within antisense lrRNA. The post-transcriptional cleavage processes followed the "tRNA punctuation" model. After the cleavage of transfer RNAs (tRNAs), COI and ND3 matured as bicistronic mRNA COI/ND3 due to the translocation of intervening tRNAs, and the 3' untranslated region (UTR) remained in the mRNAs for COII, COIII, CYTB and ND5. Additionally, isoform RNAs of ND2, srRNA and lrRNA were identified. In summary, the relative mitochondrial gene expression levels, transcriptional model and post-transcriptional cleavage process of S. houjii are notably different from those insects with typical mitochondrial gene arrangements. In addition, the phylogenetic tree of Thripidae including S. houjii was reconstructed. Our study provides insights into the phylogenetic status of Sericothripinae and the transcriptional and post-transcriptional regulation processes of extensively gene-rearranged insect mitogenomes.

Complete mitochondrial genome and phylogenetic analysis of Mancinella alouina

BACKGROUND

The Muricidae family in the Class Gastropoda comprises numerous species with a vast range of morphological features and a worldwide presence. The phylogeny of the Muricidae has been analyzed in previous studies; however, the evolutionary relationships among the main branches of the Muricidae remain unknown.

METHODS AND RESULTS

In the present study, the mitochondrial genome of Mancinella alouina was sequenced. The mitochondrial genome was found to be 16,671 bp in length and made up of 37 genes (13 protein-coding genes, 22 transfer RNA and 2 ribosomal RNA genes). The genome has an A-T-rich region (66.5% A + T content) and all of the PCGs use the ATN start codon and the TAG or TAA stop codons. The mitochondrial gene arrangement of Mancinella alouina is similar to that of other Muricidae, except for Ocinebrellus inornatus and Ceratostoma burnetti. On the basis of a flexible molecular clock model, time-calibrated phylogenetic results indicate that the genus Mancinella diverged roughly 18.09 Mya, and that the family Muricidae emerged in the Late Cretaceous.

CONCLUSIONS

This study reveals the structural and sequence information features of the mitochondrial genome of Mancinella alouina. This study provides evidence for the relationships within the family Muricidae at the molecular level, and infer the divergence time. The results of phylogenetic analyses strongly support the current classification.

The mitochondrial genome of Bottapotamon fukienense (Brachiura: Potamidae) is fragmented into two chromosomes

BACKGROUND

China is the hotspot of global freshwater crab diversity, but their wild populations are facing severe pressures associated with anthropogenic factors, necessitating the need to map their taxonomic and genetic diversity and design conservation policies.

RESULTS

Herein, we sequenced the mitochondrial genome of a Chinese freshwater crab species Bottapotamon fukienense, and found that it is fragmented into two chromosomes. We confirmed that fragmentation was not limited to a single specimen or population. Chromosome 1 comprised 15,111 base pairs (bp) and there were 26 genes and one pseudogene (pseudo-nad1) encoded on it. Chromosome 2 comprised 8,173 bp and there were 12 genes and two pseudogenes (pseudo-trnL2 and pseudo-rrnL) encoded on it. Combined, they comprise the largest mitogenome (23,284 bp) among the Potamidae. Bottapotamon was the only genus in the Potamidae dataset exhibiting rearrangements of protein-coding genes. Bottapotamon fukienense exhibited average rates of sequence evolution in the dataset and did not differ in selection pressures from the remaining Potamidae.

CONCLUSIONS

This is the first experimentally confirmed fragmentation of a mitogenome in crustaceans. While the mitogenome of B. fukienense exhibited multiple signs of elevated mitogenomic architecture evolution rates, including the exceptionally large size, duplicated genes, pseudogenisation, rearrangements of protein-coding genes, and fragmentation, there is no evidence that this is matched by elevated sequence evolutionary rates or changes in selection pressures.

Comparative Analysis of the Mitochondrial Genomes of Three Species of Yangiella (Hemiptera: Aradidae) and the Phylogenetic Implications of Aradidae

The mitochondrial genomes of three species of Yangiella were sequenced, annotated, and analyzed. The genome length of the three species of the genus is 15,070-15,202 bp, with a typical gene number, including a control region, 2 ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and 13 protein-coding genes (PCGs). It was found that the mitochondrial genome of Yangiella had AT bias. Except for the lack of a DHU arm of the trnS1 gene, the other tRNAs had a typical cloverleaf structure, and the codon usage preferences of the three species exhibited high similarity. In addition, tRNA gene rearrangements were observed among the three subfamilies of Aradidae (Mezirinae, Calisiinae, Aradinae), and it was found that codon usage preferences appeared to be less affected by base mutation and more by natural selection. The Pi and Ka/Ks values indicated that cox1 was the most conserved gene in the mitochondrial genome of Aradidae, while atp8 and nad6 were rapidly evolved genes. Substitution saturation level analysis showed that the nucleic acid sequence of mitochondrial protein-coding genes in Aradidae did not reach saturation, suggesting the rationality of the phylogenetic analysis data. Bayesian and maximum likelihood methods were used to analyze the phylogeny of 16 species of Hemiptera insects, which supported the monophyly of Aneurinae, Carventinae, and Mezirinae, as well as the monophyly of Yangiella. Based on fossils and previous studies, the differentiation time was inferred, indicating that Yangiella diverged about 57 million years ago.

De novo assembly of the complete mitochondrial genome of Mycetophylax simplex Emery, 1888 through organelle targeting revels no substantial expansion of gene spacers, but rather some slightly shorter genes

Mitochondria play a key role in cell biology and have their own genome, residing in a highly oxidative environment that induces faster changes than the nuclear genome. Because of this, mitochondrial markers have been exploited to reconstruct phylogenetic and phylogeographic relationships in studies of adaptation and molecular evolution. In this study, we determined the complete mitogenome of the fungus-farming ant Mycetophylax simplex (Hymenoptera, Formicidae) and conducted a comparative analysis among 29 myrmicine ant mitogenomes. Mycetophylax simplex is an endemic ant that inhabits sand dunes along the southern Atlantic coast. Specifically, the species occur in the ecosystem known as "restinga", within the Atlantic Forest biome. Due to habitat degradation, land use and decline of restinga habitats, the species is considered locally extinct in extremely urban beaches and is listed as vulnerable on the Brazilian Red List (ICMBio). We employed a mitochondrion-targeting approach to obtain the complete mitogenome through high-throughput DNA sequencing technology. This method allowed us to determine the mitogenome with high performance, coverage and low cost. The circular mitogenome has a length of 16,367 base pairs enclosing 37 genes (13 protein-coding genes, 22 tRNAs and 2 rRNAs) along with one control region (CR). All the protein-coding genes begin with a typical ATN codon and end with the canonical stop codons. All tRNAs formed the fully paired acceptor stems and fold into the typical cloverleaf-shaped secondary structures. The gene order is consistent with the shared Myrmicinae structure, and the A + T content of the majority strand is 81.51%. Long intergenic spacers were not found but some gene are slightly shorter. The phylogenetic relationships based on concatenated nucleotide and amino acid sequences of the 13 protein-coding genes, using Maximum Likelihood and Bayesian Inference methods, indicated that mitogenome sequences were useful in resolving higher-level relationship within Formicidae.

Mitogenomic architecture and evolution of the soil ciliates Colpoda

Colpoda are cosmopolitan unicellular eukaryotes primarily inhabiting soil and benefiting plant growth, but they remain one of the least understood taxa in genetics and genomics within the realm of ciliated protozoa. Here, we investigate the architecture of de novo assembled mitogenomes of six Colpoda species, using long-read sequencing and involving 36 newly isolated natural strains in total. The mitogenome sizes span from 43 to 63 kbp and typically contain 28-33 protein-coding genes. They possess a linear structure with variable telomeres and central repeats, with one Colpoda elliotti strain isolated from Tibet harboring the longest telomeres among all studied ciliates. Phylogenomic analyses reveal that Colpoda species started to diverge more than 326 million years ago, eventually evolving into two distinct groups. Collinearity analyses also reveal significant genomic divergences and a lack of long collinear blocks. One of the most notable features is the exceptionally high level of gene rearrangements between mitochondrial genomes of different Colpoda species, dominated by gene loss events. Population-level mitogenomic analysis on natural strains also demonstrates high sequence divergence, regardless of geographic distance, but the gene order remains highly conserved within species, offering a new species identification criterion for Colpoda species. Furthermore, we identified underlying heteroplasmic sites in the majority of strains of three Colpoda species, albeit without a discernible recombination signal to account for this heteroplasmy. This comprehensive study systematically unveils the mitogenomic structure and evolution of these ancient and ecologically significant Colpoda ciliates, thus laying the groundwork for a deeper understanding of the evolution of unicellular eukaryotes.IMPORTANCEColpoda, one of the most widespread ciliated protozoa in soil, are poorly understood in regard to their genetics and evolution. Our research revealed extreme mitochondrial gene rearrangements dominated by gene loss events, potentially leading to the streamlining of Colpoda mitogenomes. Surprisingly, while interspecific rearrangements abound, our population-level mitogenomic study revealed a conserved gene order within species, offering a potential new identification criterion. Phylogenomic analysis traced their lineage over 326 million years, revealing two distinct groups. Substantial genomic divergence might be associated with the lack of extended collinear blocks and relaxed purifying selection. This study systematically reveals Colpoda ciliate mitogenome structures and evolution, providing insights into the survival and evolution of these vital soil microorganisms.

The First Two Complete Mitochondrial Genomes for the Subfamily Meligethinae (Coleoptera: Nitidulidae) and Implications for the Higher Phylogeny of Nitidulidae

The phylogenetic status of the family Nitidulidae and its sister group relationship remain controversial. Also, the status of the subfamily Meligethinae is not fully understood, and previous studies have been mainly based on morphology, molecular fragments, and biological habits, rather than the analysis of the complete mitochondrial genome. Up to now, there has been no complete mitochondrial genome report of Meligethinae. In this study, the complete mitochondrial genomes of Meligethinus tschungseni and Brassicogethes affinis (both from China) were provided, and they were compared with the existing complete mitochondrial genomes of Nitidulidae. The phylogenetic analysis among 20 species of Coleoptera was reconstructed via PhyloBayes analysis and Maximum likelihood (ML) analysis, respectively. The results showed that the full lengths of Meligethinus tschungseni and Brassicogethes affinis were 15,783 bp and 16,622 bp, and the AT contents were 77% and 76.7%, respectively. Each complete mitochondrial genome contains 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and a control region (A + T-rich region). All the PCGs begin with the standard start codon ATN (ATA, ATT, ATG, ATC). All the PCGs terminate with a complete terminal codon, TAA or TAG, except cox1, cox2, nad4, and nad5, which terminate with a single T. Furthermore, all the tRNAs have a typical clover-leaf secondary structure except trnS1, whose DHU arm is missing in both species. The two newly sequenced species have different numbers and lengths of tandem repeat regions in their control regions. Based on the genetic distance and Ka/Ks analysis, nad6 showed a higher variability and faster evolutionary rate. Based on the available complete mitochondrial genomes, the results showed that the four subfamilies (Nitidulinae, Meligethinae, Carpophilinae, Epuraeinae) of Nitidulidae formed a monophyletic group and further supported the sister group relationship of Nitidulidae + Kateretidae. In addition, the taxonomic status of Meligethinae and the sister group relationship between Meligethinae and Nitidulinae (the latter as currently circumscribed) were also preliminarily explored.

Relaxed purifying selection pressure drives accelerated and dynamic gene rearrangements in thrips (Insecta: Thysanoptera) mitochondrial genomes

Insect mitochondrial genomes (mitogenome) generally present a typical gene order, which is considered as the ancestral arrangement. All sequenced mitogenomes in the Thysanoptera display high levels of gene rearrangement. Due to limited number of thrips mitogenomes sequenced, how gene rearrangement may be shaped by evolution remain unclear. Here, we analyzed 33 thrips mitogenomes, including 14 newly sequenced. These mitogenomes were diverse in organization, nucleotides substitution and gene arrangements. We found 28 highly rearranged gene orders with the breakpoints of gene rearrangements from 25 to 33. Reconstruction of the ancestors mitochondrial gene arrangements states indicated that Tubulifera have more complex pathways than Terebrantia in the gene order evolution. Molecular calibration estimated that divergence of two suborders occurred in the middle Triassic while the radiation of thrips was associated with the arose and flourish of angiosperm. Our evolutionary hypothesis testing suggests that relaxation of selection pressure enabled the early phase of Thysanoptera evolution, followed by a stronger selective pressure fixed diversification. Our analyses found gene inversion increases the nonsynonymous substitution rates and provide an evolutionary hypothesis driving the diverse gene orders.

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.

Characterization and phylogenetic analysis of the complete mitochondrial genome of Rhabdias kafunata (Rhabditida: Rhabdiasidae)

Rhabdias kafunata (Rhabditida: Rhabdiasidae) is a parasitic nematode that significantly affects bufonids. To better understand the genome-level characteristics of related species, Illumina sequencing was used to identify mitochondrial genes and analyze their basic characteristics and gene arrangements. The mitogenome of R. kafunata is 14,068 bp in length and contains 36 genes, including 12 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes and one noncoding region (NCR). The nucleotide composition is highly biased toward A + T, accounting for 75.5% of the entire mitochondrial genome. The cox1 sequence is relatively conserved in Ka/Ks analyses and can be used as a gene fragment for species identification. While 8 of the 12 PCGs use the typical ATN initiation codon, nad1-2, nad4, and cox3 utilize a TTG initiation codon. Most stop codons end with the standard TAA or TAG, except for cytb, which ends with an incomplete TA. Additionally, trnM, trnK, and trnI have the typical clover-leaf secondary structure, while the remaining tRNAs lack the DHU arm or TΨC arm. Phylogenetic analysis indicates that R. kafunata belongs to the Rhabditidae family and is closely related to Litoditis marina and Caenorhabditis angaria among sequenced lepidopteran mitochondrial genomes.

Variants in the mitochondrial genome sequence of Oryctes rhinoceros (Coleoptera: Scarabaeidae) infected with Oryctes rhinoceros nudivirus in oil palm and coconut plantations

The CRB (coconut rhinoceros beetle) haplotype was classified into CRB-S and CRB-G, based on the presence of single nucleotide polymorphisms (SNPs) in the mitochondrial cox1 gene. Mitochondrial genomes (mitogenomes) are the most widely used genetic resources for molecular evolution, phylogenetics, and population genetics in relation to insects. This study presents the mitogenome CRB-G and CRB-S which were collected in Johor, Malaysia. The mitogenome of CRB-G collected from oil palm plantations in 2020 and 2021, and wild coconut palms in 2021 was 15,315 bp, 15,475 bp, and 17,275 bp, respectively. The CRB-S was discovered in coconut and oil palms in 2021, and its mitogenome was 15,484 bp and 17,142 bp, respectively. All the mitogenomes have 37 genes with more than 99% nucleotide sequence homology, except the CRB-G haplotype collected from oil palm in 2021 with 89.24% nucleotide sequence homology. The mitogenome of Johor CRBs was variable in the natural population due to its elevated mutation rate. Substitutions and indels in cox1, cox2, nad2 and atp6 genes were able to distinguish the Johor CRBs into two haplotypes. The mitogenome data generated in the present study may provide baseline information to study the infection and relationship between the two haplotypes of Johor CRB and OrNV in the field. This study is the first report on the mitogenomes of mixed haplotypes of CRB in the field.

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.

Complete mitochondrial genomes of the slugs Deroceraslaeve (Agriolimacidae) and Ambigolimaxvalentianus (Limacidae) provide insights into the phylogeny of Stylommatophora (Mollusca, Gastropoda)

In this study, we sequenced two complete mitogenomes from Deroceraslaeve and Ambigolimaxvalentianus. The mitogenome of Ambigolimaxvalentianus represented the first such data from the family Limacidae. The lengths of the mitogenomes of Deroceraslaeve and Ambigolimaxvalentianus were 14,773 bp and 15,195 bp, respectively. The entire set of 37 mitochondrial genes were identified for both mitogenomes. Compared with the mitogenome of Achatinafulica, the trnP_trnA tRNA cluster was rearranged in both Deroceraslaeve and Ambigolimaxvalentianus. The secondary structures of tRNA and rRNA genes for the two species were predicted. Phylogenetic analyses based on amino acid sequences supported (1) monophyly of Stylommatophora, (2) division of Stylommatophora into the 'achatinoid' clade (i.e., the suborder Achatinina) and the 'non-achatinoid' clade (i.e., the suborder Helicina), (3) placement of the Orthurethra in the 'non-achatinoid' clade, and (4) monophyly of each of the superfamilies Helicoidea, Urocoptoidea, Succineoidea, Arionoidea, Pupilloidea and Limacoidea. The exemplars of Helicidae, Philomycidae and Achatinellidae displayed many more mitochondrial gene rearrangements than other species of Stylommatophora.

Evolution of mitogenomic gene order in Orthoptera

Mitochondrial gene order has contributed to the elucidation of evolutionary relationships in several animal groups. It generally has found its application as a phylogenetic marker for deep nodes. Yet, in Orthoptera limited research has been performed on the gene order, although the group represents one of the oldest insect orders. We performed a comprehensive study on mitochondrial genome rearrangements (MTRs) within Orthoptera in the context of mitogenomic sequence-based phylogeny. We used 280 published mitogenome sequences from 256 species, including three outgroup species, to reconstruct a molecular phylogeny. Using a heuristic approach, we assigned MTR scenarios to the edges of the phylogenetic tree and reconstructed ancestral gene orders to identify possible synapomorphies in Orthoptera. We found all types of MTRs in our dataset: inversions, transpositions, inverse transpositions, and tandem-duplication/random loss events (TDRL). Most of the suggested MTRs were in single and unrelated species. Out of five MTRs which were unique in subgroups of Orthoptera, we suggest four of them to be synapomorphies; those were in the infraorder Acrididea, in the tribe Holochlorini, in the subfamily Pseudophyllinae, and in the two families Phalangopsidae and Gryllidae or their common ancestor (leading to the relationship ((Phalangopsidae + Gryllidae) + Trigonidiidae)). However, similar MTRs have been found in distant insect lineages. Our findings suggest convergent evolution of specific mitochondrial gene orders in several species, deviant from the evolution of the mitogenome DNA sequence. As most MTRs were detected at terminal nodes, a phylogenetic inference of deeper nodes based on MTRs is not supported. Hence, the marker does not seem to aid resolving the phylogeny of Orthoptera, but adds further evidence for the complex evolution of the whole group, especially at the genetic and genomic levels. The results indicate a high demand for more research on patterns and underlying mechanisms of MTR events in Orthoptera.

New Dielis species and structural dichotomy of the mitochondrial cox2 gene in Scoliidae wasps

Some mitochondrial protein-coding genes of protists and land plants have split over the course of evolution into complementary genes whose products can form heteromeric complexes that likely substitute for the undivided proteins. One of these genes, cox2, has also been found to have split in animals, specifically in Scoliidae wasps (Hymenoptera: Apocrita) of the genus Dielis (Campsomerini), while maintaining the conventional structure in related Scolia (Scoliini). Here, a hitherto unrecognized Nearctic species of Dielis, D. tejensis, is described based on its phenotype and mtDNA. The mitogenome of D. tejensis sp. nov. differs from that of the sympatric sibling species Dielis plumipes fossulana by the reduced size of the cox2-dividing insert, which, however, still constitutes the fifth part of the mtDNA; an enlarged nad2-trnW intergenic region; the presence of two trnK^ttt paralogues; and other features. Both species of Dielis have a unique insertion of a threonine in COXIIA, predicted to be involved in COXIIA-COXIIB docking, and substitutions of two hydrophobic residues with redox-active cysteines around the Cu_A centre in COXIIB. Importantly, the analysis of mtDNA from another Campsomerini genus, Megacampsomeris, shows that its cox2 gene is also split. The presented data highlight evolutionary processes taking place in hymenopteran mitogenomes that do not fall within the mainstream of animal mitochondrion evolution.

Comparative analysis of the mitochondrial genome of Dermacentor steini from different regions in China

Ticks are a group of blood-sucking ectoparasites that play an important role in human health and livestock production development as vectors of zoonotic diseases. The phylogenetic tree of single genes cannot accurately reflect the true kinship between species. Based on the complete mitochondrial genome analysis one can help to elucidate the phylogenetic relationships among species. In this study, the complete mitochondrial genome of Dermacentor steini (isolate Longyan) was sequenced and compared with the mitochondrial genes of 3 other Chinese isolates (Nanchang, Jinhua and Yingtan). In Dermacentor steini 4 isolates had identical or similar mitochondrial genome lengths and an overall variation of 0.76% between sequences. All nucleotide compositions showed a distinct AT preference. The most common initiation and stop codons were ATG and TAA, respectively. Fewer base mismatches were found in the tRNA gene of D. steini (isolate Longyan), and the vicinity of the control region and tRNA gene was a hot rearrangement region of the genus Dermacentor. Maximum likelihood trees and Bayesian trees indicate that D. steini is most closely related to Dermacentor auratus. The results enrich the mitochondrial genomic data of species in the genus Dermacentor and provide novel insights for further studies on the phylogeographic classification and molecular evolution of ticks.

The complete mitochondrial genome and novel gene arrangement in Nesodiprion zhejiangensis Zhou & Xiao (Hymenoptera: Diprionidae)

The complete mitochondrial genome (mitogenome) of the sawfly, Nesodiprion zhejiangensis Zhou & Xiao, was sequenced, assembled, and deposited in GenBank (Accession Number: OM501121). The 15,660 bp N. zhejiangensis mitogenome encodes for 2 ribosomal RNAs (rrnL and rrnS), 22 transfer RNAs (tRNAs), 13 protein-coding genes (PCGs), and an AT-rich region of 450 bp in length. The nucleotide composition is biased toward adenine and thymine (A + T = 81.8%). Each PCG is initiated by an ATN codon, except for cox2, which starts with a TTG. Of 13 PCGs, 9 have a TAA termination codon, while the remainder terminate with a TAG or a single T. All tRNAs have the classic cloverleaf structure, except for the dihydrouridine (DHU) arm of tRNA^val, which forms a simple loop. There are 49 helices belonging to 6 domains in rrnL and 30 helices belonging to 4 domains in rrnS. In comparison to the ancestral architecture, N. zhejiangensis has the most rearranged mitogenome in Symphyta, in which rearrangement events of local inversion and transposition are identified in three gene clusters. Specifically, the main hotspot of gene rearrangement occurred between rrnS and trnY, and rearranged from rrnS-(AT-rich region)-I-Q-M-nd2-W-C-Y to rrnS-Q-W-C-nd2-I-M-(AT-rich region)-Y, involving a local inversion event of a large gene cluster and transposition events of some tRNAs. Transposition of trnA and trnR (rearranged from A-R to R-A) was observed at the nd3-nd5 gene junction while shuffling of trnP and trnT (rearranged from T-P to P-T) occurred at the nd4l-nd6 gene junction. While illegitimate inter-mtDNA recombination might explain the opposite orientations of transcription between rrnS and trnY, transposition events of tRNA in some gene blocks can be accounted for by the tandem duplication/random loss (TDRL) model. Our phylogenetic analysis suggests that N. zhejiangensis is closely related to congeneric species N. biremis and N. japonicus, which together form a sister lineage with the European pine sawfly, Neodiprion sertifer.

A novel, conserved and possibly functional motif "WHWGHTW" in mitochondrial transcription across Bilateria

The animal mitogenomes which undergone a reductive evolution has an obvious loss of coding capacity compared to their known closest relatives, but it has not yet been fully investigated why and how the intergenic regions do not encode protein and have no known functions, are stably maintained, replicated, and transmitted by the genome. These relatively small intergenic regions may not be under neutral evolution and they may have functional and/or regulatory roles that have yet to be identified. Here, the distribution pattern, sequence content and location of a novel sequence motif of 'WWWGHTW' were bioinformatically investigated and characterised by constructing a sampling mitogenome dataset of 1889 species from 14 phyla representing the clade of Bilateria. This motif is reverse complementary of the previously described DmTTF binding sequence and found in the nd4L- (X) -trnT gene cluster. This cluster commonly exhibits a strand displacement region and an intergenic region among the bilaterian superphylums, particularly in Ecdysozoa. This motif may be accepted as a substrate providing binding sites for the specific interaction with transcription factors because of (i) its reverse complementarity of previously described DmTTF binding sequence, and (ii) the possession of G and T nucleotides in the fourth and sixth positions, (iii) the bias on T and G nucleotides instead of C and A in the degenerated positions. This suggestion is also supported by the presence of a strand displacement region in the nd4L- (X) -trnT gene cluster, particularly in Ecdysozoa consisting of the most rearranged mitogenomes among the bilaterian superphylums.

Complete Mitogenome of Oreolalax Omeimontis Reveals Phylogenetic Status and Novel Gene Arrangement of Archaeobatrachia

Species of the genus Oreolalax displayed crucial morphological characteristics of vertebrates transitioning from aquatic to terrestrial habitats; thus, they can be regarded as a representative vertebrate genus for this landing phenomenon. But the present phylogenetic status of Oreolalax omeimontis has been controversial with morphological and molecular approaches, and specific gene rearrangements were discovered in all six published Oreolalax mitogenomes, which are rarely observed in Archaeobatrachia. Therefore, this study determined the complete mitogenome of O. omeimontis with the aim of identifying its precise phylogenetic position and novel gene arrangement in Archaeobatrachia. Phylogenetic analysis with Bayesian inference and maximum likelihood indicates O. omeimontis is a sister group to O. lichuanensis, which is consistent with previous phylogenetic analysis based on morphological characteristics, but contrasts with other studies using multiple gene fragments. Moreover, although the duplication of trnM occurred in all seven Oreolalax species, the translocation of trnQ and trnM occurred differently in O. omeimontis to the other six, and this unique rearrangement would happen after the speciation of O. omeimontis. In general, this study sheds new light on the phylogenetic relationships and gene rearrangements of Archaeobatrachia.

Evolution of Gene Arrangements in the Mitogenomes of Ensifera and Characterization of the Complete Mitogenome of Schizodactylus jimo

Gene arrangement (relative location of genes) is another evolutionary marker of the mitogenome that can provide extensive information on the evolutionary mechanism. To explore the evolution of gene arrangements in the mitogenome of diversified Ensifera, we sequenced the mitogenome of the unique dune cricket species found in China and used it for phylogenetic analysis, in combination with 84 known Ensiferan mitogenomes. The mitogenome of Schizodactylus jimo is a 16,428-bp circular molecule that contains 37 genes. We identified eight types of gene arrangement in the 85 ensiferan mitogenomes. The gene location changes (i.e., gene translocation and duplication) were in three gene blocks: I-Q-M-ND2, rrnl-rns-V, and ND3-A-R-N-S-E-F. From the phylogenetic tree, we found that Schizodactylus jimo and most other species share a typical and ancient gene arrangement type (Type I), while Grylloidea has two types (Types II and III), and the other five types are rare and scattered in the phylogenetic tree. We deduced that the tandem replication–random loss model is the evolutionary mechanism of gene arrangements in Ensifera. Selection pressure analysis revealed that purifying selection dominated the evolution of the ensiferan mitochondrial genome. This study suggests that most gene rearrangements in the ensiferan mitogenome are rare accidental events.

Novel gene rearrangement pattern in mitochondrial genome of Ooencyrtusplautus Huang & Noyes, 1994: new gene order in Encyrtidae (Hymenoptera, Chalcidoidea)

Studies of mitochondrial genomes have a wide range of applications in phylogeny, population genetics, and evolutionary biology. In this study, we sequenced and analyzed the mitochondrial genome of Ooencyrtusplautus Huang & Noyes, 1994 (Hymenoptera, Encyrtidae). The nearly complete mitogenome of O.plautus was 15,730 bp in size, including 13 PCGs (protein-coding genes), 22 tRNAs, 2 rRNAs, and a nearly complete control region. The nucleotide composition was significantly biased toward adenine and thymine, with an A + T content of 84.6%. We used the reference sequence of Chouioiacunea and calculated the Ka/Ks ratio for each set of PCGs. The highest value of the Ka/Ks ratio within 13 PCGs was found in nad2 with 1.1, suggesting that they were subjected to positive selection. This phenomenon was first discovered in Encyrtidae. Compared with other encyrtid mitogenomes, a translocation of trnW was found in O.plautus, which was the first of its kind to be reported in Encyrtidae. Comparing with ancestral arrangement pattern, wasps reflect extensive gene rearrangements. Although these insects have a high frequency of gene rearrangement, species from the same family and genus tend to have similar gene sequences. As the number of sequenced mitochondrial genomes in Chalcidoidea increases, we summarize some of the rules of gene rearrangement in Chalcidoidea, that is four gene clusters with frequent gene rearrangements. Ten mitogenomes were included to reconstruct the phylogenetic trees of Encyrtidae based on both 13 PCGs (nucleotides of protein coding genes) and AA matrix (amino acids of protein coding genes) using the maximum likelihood and Bayesian inference methods. The phylogenetic tree reconstructed by Bayesian inference based on AA data set showed that Aenasiusarizonensis and Metaphycuseriococci formed a clade representing Tetracneminae. The remaining six species formed a monophyletic clade representing Encyrtinae. In Encyrtinae, Encyrtus forms a monophyletic clade as a sister group to the clade formed by O.plautus and Diaphorencyrtusaligarhensis. Encyrtussasakii and Encyrtusrhodooccisiae were most closely related species in this monophyletic clade. In addition, gene rearrangements can provide a valuable information for molecular phylogenetic reconstruction. These results enhance our understanding of phylogenetic relationships among Encyrtidae.

A Rearrangement of the Mitochondrial Genes of Centipedes (Arthropoda, Myriapoda) with a Phylogenetic Analysis

Due to the limitations of taxon sampling and differences in results from the available data, the phylogenetic relationships of the Myriapoda remain contentious. Therefore, we try to reconstruct and analyze the phylogenetic relationships within the Myriapoda by examining mitochondrial genomes (the mitogenome). In this study, typical circular mitogenomes of Mecistocephalus marmoratus and Scolopendra subspinipes were sequenced by Sanger sequencing; they were 15,279 bp and 14,637 bp in length, respectively, and a control region and 37 typical mitochondrial genes were annotated in the sequences. The results showed that all 13 PCGs started with ATN codons and ended with TAR codons or a single T; what is interesting is that the gene orders of M. marmoratus have been extensively rearranged compared with most Myriapoda. Thus, we propose a simple duplication/loss model to explain the extensively rearranged genes of M. marmoratus, hoping to provide insights into mitogenome rearrangement events in Myriapoda. In addition, our mitogenomic phylogenetic analyses showed that the main myriapod groups are monophyletic and supported the combination of the Pauropoda and Diplopoda to form the Dignatha. Within the Chilopoda, we suggest that Scutigeromorpha is a sister group to the Lithobiomorpha, Geophilomorpha, and Scolopendromorpha. We also identified a close relationship between the Lithobiomorpha and Geophilomorpha. The results also indicate that the mitogenome can be used as an effective mechanism to understand the phylogenetic relationships within Myriapoda.

Phylogenomic Analyses of the Tenthredinoidea Support the Familial Rank of Athaliidae (Insecta, Tenthredinoidea)

The systematic status of the genus Athalia and related genera is a perennial controversy in sawfly taxonomy. Several authors have hypothesized that the placement of Athalia within the Tenthredinidae is artificial, but no studies have focused on this topic. If the hypothesis that Athalia does not belong to Tenthredinidae can be supported, the taxonomic framework of Tenthredinoidea needs revision. We present a comprehensive phylogenomic study of Tenthredinoidae, focusing on the positions of Athalia and related genera by sampling 80 representatives mainly of the Tenthredinoidea, including Heptamelinae and Blasticotomidae. Our phylogenetic reconstructions based on nuclear genes and mitochondrial (mt) sequences support Athalia and related genera as a distinct clade sister to Tenthredinidae + (Cimbicidae + Diprionidae). A comparison of symphytan mitochondrial genomes reveals an innovative gene rearrangement pattern in Athaliidae, in which Dentathalia demonstrates a more ancestral pattern than Athalia and Hypsathalia. The lineage specificity of mt rRNA secondary structures also provides sufficient support to consider Athaliidae as a separate family. In summary, the phylogeny and genomic structural changes unanimously support the taxonomic treatment of Athaliidae as a family and the re-establishment of Dentathalia as a valid genus.

The complete mitochondrial genome of Syrista parreyssii (Spinola, 1843) (Hymenoptera: Cephidae) and its phylogenetic analyses

The complete mitochondrial genome of Syrista parreyssii (Spinola, 1843) was described. The circular genome is 18,666 bp with an A + T content of 82.60%. It contains 37 genes and a 1921 bp control region. The CR-trnI (+)-trnQ (–)-trnM (+) cluster is rearranged as trnM (+)-CR-trnQ (–)-trnI (+) cluster. Phylogenetic analysis demonstrates that European Syrista and Asian Neosyrista were not sister groups. Neosyrista is a valid genus and should be reestablished. Moreover, a preliminary study based on COI showed there are at least three valid Syrista species within the European and Mediterranean regions. Whether the known Syrista parreyssii (Spinola, 1843) is a complex or there are more cryptic species needs further study.

Insights into the Evolution of Aphid Mitogenome Features from New Data and Comparative Analysis

The complete mitochondrial genomes and their rearrangement patterns can provide useful information for inferring evolutionary history of organisms. Aphids are one of the insect groups with some unique mitogenome features. In this study, to examine whether some features in aphid mitogenomes are independent species-specific evolutionary events or clade-specific events at certain taxonomic levels, we sequenced three new aphid mitogenomes (Hormaphidinae: Ceratovacuna keduensis, Pseudoregma panicola; Lachninae: Nippolachnus piri) and compared them with all known aphid mitogenomes. The three mitogenomes are 16,059–17,033 bp in length, with a set of 37 typical mitochondrial genes, a non-coding control region and a tandem repeat region. The gene orders of them are all highly rearranged. Within the subfamily Hormaphidinae, the presence of repeat region and mitogenome rearrangement in Cerataphidini species but not in the other two tribes indicate that these may be Cerataphidini-specific features. The same gene rearrangement pattern in the two Lachninae species, N. piri (Tuberolachnini) and Stomaphis sinisalicis (Stomaphidini), supports that this feature should be at least derived from the common ancestor of two tribes. Overall, our data and analyses provide new insights into the evolutionary patterns of gene rearrangement and repeat region in aphid mitogenomes, and further corroborate the potential role of gene rearrangement in elucidating the evolutionary history of different insect lineages.

The complete mitochondrial genome of Platygaster robiniae (Hymenoptera: Platygastridae): A novel tRNA secondary structure, gene rearrangements and phylogenetic implications

Platygaster robiniae is economically important as a highly specific parasitoid of the invasive pest Obolodiplosis robiniae which was introduced into the Euro-Asia region in the last decade. Despite being a critical and specific parasitoid of the invasive pest O. robiniae and its use as an effective biocontrol agent, the absence of sequence information from P. robiniae have limited its genetic applications for pest management in forests. Mitochondrial (mt) genomes generally contain abundant nucleotide information and thus are helpful for understanding species history. Here, we sequenced the complete mt genome of P. robiniae using next generation sequencing, and annotated 13 protein-coding, 22 tRNA, and 2 rRNA genes and a 702 bp noncoding region. Comparative analysis indicated that this mt genome has a normal A + T content and codons use, however possessed both the expected and unique rearrangements. Ten tRNAs at four gene blocks COII-ATP8, COIII-ND3, ND3-ND5 and the A + T-rich region-ND2 were rearranged, including gene shuffles, transpositions and inversions. Notably, two genes tRNA^Ser(UCN) and tRNA^Leu(CUN) had undergone long-range inversions, which is the first record of this rearrangement type in the superfamily Platygastroidea. The D-loops of both tRNA^Ile and tRNA^Leu(CUN) were absent from the tRNA secondary structure, which has not been reported from hymenopteran previously. Phylogenetic analysis based with the maximum likelihood and Bayesian methods showed that P. robiniae grouped with other species of Platygastridae, and that the superfamily Platygastridea is sister to the other Proctotrupomorpha superfamilies. Our tree strongly supports the monophyly of the five superfamilies of Proctotrupomorpha. This study discovered some unique characters of P. robiniae, and contributes to our understanding of genome rearrangements in the order Hymenoptera.

•

Assemble the complete mitochondria genome of the obligate parasitoid (Platygaster robiniae) of Obolodiplosis robiniae.

•

Discover some unique tRNA secondary structures.

•

Discover some unique rearrangements of Platygaster robiniae and Platygastroidea.

•

The abundance rearrangements in the order Hymenoptera.

Complete Nucleotide Sequence of the Mitogenome of Tapinoma ibericum (Hymenoptera: Formicidae: Dolichoderinae), Gene Organization and Phylogenetics Implications for the Dolichoderinae Subfamily

The ant Tapinoma ibericum Santschi, 1925 is native to the Iberian Peninsula. This species, as well as other species from the Tapinoma nigerrimum complex, could form supercolonies that make these species potentially invasive and could give rise to pests. Recently a mature colony from this species has been found in the Isle of Wight (United Kingdom). Mitogenomes have been used to study the taxonomy, biogeography and genetics of species, improving the development of strategies against pest invasion. However, the number of available mitogenomes from the subfamily Dolichoderinae is still scarce and only two of these mitogenomes belong to Tapinoma species. Herein, the complete mitogenome of T. ibericum is presented in order to increase the molecular information of the genus. The T. ibericum mitogenome, retrieved by Next-Generation Sequencing data, is 15,715 bp in length. It contains the typical set of 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNAs and the A + T-rich control region. Comparisons of the T. ibericum mitogenome with other dolichoderine mitogenomes revealed the existence of four gene rearrangements in relation with the ancestral insect mitogenome. One of these rearrangements, involving the tRNA-Ile, tRNA-Gln and tRNA-Met genes, was found in most of the analyzed ant mitogenomes. Probably this rearrangement was an ancestral or plesiomorphic character in Formicidae. Interestingly, another rearrangement that affects to tRNA-Trp, tRNA-Cys and tRNA-Tyr genes was found only in Tapinoma species. This change could be a synapomorphic character for the genus Tapinoma, and could be used as a phylogenetic marker. Additionally, a phylogenetic analysis was performed using the protein-coding gene sequences from available Dolichoderinae mitogenomes, as well as mitogenomes from representative species from other Formicidae subfamilies. Results support the monophyletic nature of the genus Tapinoma placing it within the same clade as the rest of Dolichoderinae species.

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.

Exploring mitogenome evolution in Branchiopoda (Crustacea) lineages reveals gene order rearrangements in Cladocera

The class Branchiopoda, whose origin dates back to Cambrian, includes ~ 1200 species which mainly occupy freshwater habitats. The phylogeny and systematics of the class have been debated for long time, until recent phylogenomic analyses allowed to better clarify the relationships among major clades. Based on these data, the clade Anostraca (fairy and brine shrimps) is sister to all other branchiopods, and the Notostraca (tadpole shrimps) results as sister group to Diplostraca, which includes Laevicaudata + Spinicaudata (clam shrimps) and Cladoceromorpha (water fleas + Cyclestherida). In the present analysis, thanks to an increased taxon sampling, a complex picture emerges. Most of the analyzed mitogenomes show the Pancrustacea gene order while in several other taxa they are found rearranged. These rearrangements, though, occur unevenly among taxa, most of them being found in Cladocera, and their taxonomic distribution does not agree with the phylogeny. Our data also seems to suggest the possibility of potentially homoplastic, alternative gene order within Daphniidae.

Mitochondrial Genomes Yield Insights into the Basal Lineages of Ichneumonid Wasps (Hymenoptera: Ichneumonidae)

We obtained four mitochondrial genomes of Odontocolon sp., Xorides funiuensis, Euceros kiushuensis and Euceros serricornis, which represent the first two representatives from subfamily Eucerotinae and Xoridinae (Ichneumonidae), respectively. All of the 4 newly sequenced mitochondrial genomes contain 13 protein-coding genes (PCGs) and most 24 RNA genes. Furthermore, they all have novel tRNA rearrangement patterns comparing with published mitogenomes in Ichneumonidae. For the tRNA cluster trnI-trnQ-trnM, X. funiuensis is shuffled as trnM-trnI-trnQ with trnQ inversed, while Odontocolon sp. with a remote translocation of trnK, shuffling as trnI-trnM-trnQ. E. kiushuensis and E. serricornis shared the same cluster trnQ-trnY-trnW-trnC. Finally, we reconstructed the phylogenetic relationships among the sequenced subfamilies of Ichneumonidae based on nucleotides and amino acids sequences of 13 PCGs in mitochondrial genomes, and the results of both the maximum likelihood and Bayesian inference analyses highly support that Eucerotinae is the basal ichneumonid lineage rather than Xoridinae.

Comparative Mitochondrial Genomics of 104 Darwin Wasps (Hymenoptera: Ichneumonidae) and Its Implication for Phylogeny

Simple Summary

Nearly a hundred mitochondrial genomes of ichneumonid wasps are newly reported. Comparative mitogenomics of 104 mitochondrial genomes representing 33 subfamilies of Ichneumonidae, as well as its implications for phylogeny, were studied. We found that the mitochondrial genomes of ichneumonid wasps were highly conserved in their base composition and had low evolutionary rates, but were diverse in gene order. There are 38 types of gene rearrangement events in 104 ichneumonid mitochondrial genomes, of which 30 novel rearrangement types (R3-6, R8-R10, R12-R15, R17-R18, R20-R35 and R38) and a hot spot rearrangement around R1, with a shuffled tRNA cluster trnW-trnY-trnC and trnI-trnQ-trnM, were detected. The relationships among these subfamilies are firstly discussed based on mitochondrial genomes at a large scale. We suggest five subfamily groupings of Ichneumonidae: Brachycyrtiformes, Ichneumoniformes, Ophioniformes, Pimpliformes and Xoridiformes. Two formerly unplaced subfamilies, Eucerotinae and Microleptinae, were placed in Brachycyrtiformes and Ichneumoniformes, respectively.

Abstract

Ichneumonidae is one of the largest families of insects with a mega-diversity of specialized morphological and biological characteristics. We newly sequenced 92 mitochondrial genomes of ichneumonid wasps and found that they have a conserved base composition and a lower evolutionary rate than that of other families of parasitic Hymenoptera. There are 38 types of gene order in the ichneumonid mitochondrial genome, with 30 novel types identified in 104 ichneumonids. We also found that the rearrangement events occur more frequently in Ophioniformes than in Ichneumoniformes and Pimpliformes. Furthermore, the higher Ophioniformes and their relative lineages shared the transposition of trnL2 to trnI-trnQ-trnM tRNA cluster. We confirmed five higher-level groupings of Ichneumonidae: Brachycyrtiformes, Ichneumoniformes, Ophioniformes, Pimpliformes and Xoridiformes. Two formerly unplaced subfamilies, Eucerotinae and Microleptinae, were placed in Brachycyrtiformes and Ichneumoniformes, respectively. The results will improve our understanding of the diversity and evolution of Ichneumonidae.

Novel gene rearrangement in the mitochondrial genome of Anastatus fulloi (Hymenoptera Chalcidoidea) and phylogenetic implications for Chalcidoidea

The genus Anastatus comprises a large group of parasitoids, including several biological control agents in agricultural and forest systems. The taxonomy and phylogeny of these species remain controversial. In this study, the mitogenome of A. fulloi Sheng and Wang was sequenced and characterized. The nearly full-length mitogenome of A. fulloi was 15,692 bp, compromising 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes and a control region (CR). The total A + T contents were 83.83%, 82.18%, 87.58%, 87.27%, and 82.13% in the whole mitogenome, 13 PCGs, 22 tRNA genes, 2 rRNA genes, and CR, respectively. The mitogenome presented negative AT skews and positive GC skews, except for the CR. Most PCGs were encoded on the heavy strand, started with ATN codons, and ended with TAA codons. Among the 3736 amino acid-encoding codons, TTA (Leu1), CGA (Arg), TCA (Ser2), and TCT (Ser2) were predominant. Most tRNAs had cloverleaf secondary structures, except trnS1, with the absence of a dihydrouridine (DHU) arm. Compared with mitogenomes of the ancestral insect and another parasitoid within Eupelmidae, large-scale rearrangements were found in the mitogenome of A. fulloi, especially inversions and inverse transpositions of tRNA genes. The gene arrangements of parasitoid mitogenomes within Chalcidoidea were variable. A novel gene arrangement was presented in the mitogenome of A. fulloi. Phylogenetic analyses based on the 13 protein-coding genes of 20 parasitoids indicated that the phylogenetic relationship of 6 superfamilies could be presented as Mymaridae + (Eupelmidae + (Encyrtidae + (Trichogrammatidae + (Pteromalidae + Eulophidae)))). This study presents the first mitogenome of the Anastatus genus and offers insights into the identification, taxonomy, and phylogeny of these parasitoids.

Substantial rearrangements, single nucleotide frameshift deletion and low diversity in mitogenome of Wolbachia-infected strepsipteran endoparasitoid in comparison to its tephritid hosts

Insect mitogenome organisation is highly conserved, yet, some insects, especially with parasitic life cycles, have rearranged mitogenomes. Furthermore, intraspecific mitochondrial diversity can be reduced by fitness-affecting bacterial endosymbionts like Wolbachia due to their maternal coinheritance with mitochondria. We have sequenced mitogenomes of the Wolbachia-infected endoparasitoid Dipterophagus daci (Strepsiptera: Halictophagidae) and four of its 22 known tephritid fruit fly host species using total genomic extracts of parasitised flies collected across > 700 km in Australia. This halictophagid mitogenome revealed extensive rearrangements relative to the four fly mitogenomes which exhibited the ancestral insect mitogenome pattern. Compared to the only four available other strepsipteran mitogenomes, the D. daci mitogenome had additional transpositions of one rRNA and two tRNA genes, and a single nucleotide frameshift deletion in nad5 requiring translational frameshifting or, alternatively, resulting in a large protein truncation. Dipterophagus daci displays an almost completely endoparasitic life cycle when compared to Strepsiptera that have maintained the ancestral state of free-living adults. Our results support the hypothesis that the transition to extreme endoparasitism evolved together with increased levels of mitogenome changes. Furthermore, intraspecific mitogenome diversity was substantially smaller in D. daci than the parasitised flies suggesting Wolbachia reduced mitochondrial diversity because of a role in D. daci fitness.

The Mitochondrial Genomes of 18 New Pleurosticti (Coleoptera: Scarabaeidae) Exhibit a Novel trnQ-NCR-trnI-trnM Gene Rearrangement and Clarify Phylogenetic Relationships of Subfamilies within Scarabaeidae

The availability of next-generation sequencing (NGS) in recent years has facilitated a revolution in the availability of mitochondrial (mt) genome sequences. The mt genome is a powerful tool for comparative studies and resolving the phylogenetic relationships among insect lineages. The mt genomes of phytophagous scarabs of the subfamilies Cetoniinae and Dynastinae were under-represented in GenBank. Previous research found that the subfamily Rutelinae was recovered as a paraphyletic group because the few representatives of the subfamily Dynastinae clustered into Rutelinae, but the subfamily position of Dynastinae was still unclear. In the present study, we sequenced 18 mt genomes from Dynastinae and Cetoniinae using next-generation sequencing (NGS) to re-assess the phylogenetic relationships within Scarabaeidae. All sequenced mt genomes contained 37 sets of genes (13 protein-coding genes, 22 tRNAs, and two ribosomal RNAs), with one long control region, but the gene order was not the same between Cetoniinae and Dynastinae species. All mt genomes of Dynastinae species showed the same gene rearrangement of trnQ-NCR-trnI-trnM, whereas all mt genomes of Cetoniinae species showed the ancestral insect gene order of trnI-trnQ-trnM. Phylogenetic analyses (IQ-tree and MrBayes) were conducted using 13 protein-coding genes based on nucleotide and amino acid datasets. In the ML and BI trees, we recovered the monophyly of Rutelinae, Cetoniinae, Dynastinae, and Sericinae, and the non-monophyly of Melolonthinae. Cetoniinae was shown to be a sister clade to (Dynastinae + Rutelinae).

Complete Mitochondrial Genomes of Metcalfa pruinosa and Salurnis marginella (Hemiptera: Flatidae): Genomic Comparison and Phylogenetic Inference in Fulgoroidea

The complete mitochondrial genomes (mitogenomes) of two DNA barcode-defined haplotypes of Metcalfa pruinosa and one of Salurnis marginella (Hemiptera: Flatidae) were sequenced and compared to those of other Fulgoroidea species. Furthermore, the mitogenome sequences were used to reconstruct phylogenetic relationships among fulgoroid families. The three mitogenomes, including that of the available species of Flatidae, commonly possessed distinctive structures in the 1702-1836 bp A+T-rich region, such as two repeat regions at each end and a large centered nonrepeat region. All members of the superfamily Fulgoroidea, including the Flatidae, consistently possessed a motiflike sequence (TAGTA) at the ND1 and trnS₂ junction. The phylogenetic analyses consistently recovered the familial relationships of (((((Ricaniidae + Issidae) + Flatidae) + Fulgoridae) + Achilidae) + Derbidae) in the amino acid-based analysis, with the placement of Cixiidae and Delphacidae as the earliest-derived lineages of fulgoroid families, whereas the monophyly of Delphacidae was not congruent between tree-constructing algorithms.

Novel tRNA gene rearrangements in the mitochondrial genomes of praying mantises (Mantodea: Mantidae): Translocation, duplication and pseudogenization

Gene rearrangements have been found in several mitochondrial genomes of Mantodea, located in the gene blocks CR-I-Q-M-ND2, COX1-K-D-ATP8 and ND3-A-R-N-S-E-F-ND5. We have sequenced one mitogenome of Amelidae (Yersinia mexicana) and six mitogenomes of Mantidae to discuss the mitochondrial gene rearrangement and the phylogenetic relationship within Mantidae. These mitogenomes showed rearrangements of tRNA genes except for Asiadodis yunnanensis and Hierodula zhangi. These novel gene rearrangements of Mantidae were primarily concentrated in the region of CR-I-Q-M-ND2, including gene translocation, duplication and pseudogenization. For the occurrences of these rearrangements, the tandem duplication-random loss (TDRL) model and slipped-strand mispairing model were suitable to explain. Large non-coding regions (LNCRs) located in the region of CR-I-Q-M-ND2 were detected in most Mantidae species, whereas some LNCRs had high similarity to the control region (CR). Both BI and ML phylogenetic analyses supported the monophyly of Mantidae and the paraphyly of Mantinae. The phylogenetic results with the gene order and the location of NCRs acted as forceful evidence that specific gene rearrangements and special LNCRs may be synapomorphies for several groups of mantises.

Mitochondrial Phylogenomics of Tenthredinidae (Hymenoptera: Tenthredinoidea) Supports the Monophyly of Megabelesesinae as a Subfamily

Tenthredinidae represents one of the external feeders of the most diverse superfamily, Tenthredinoidea, with diverse host plant utilization. In this study, four complete mitochondrial genomes (mitogenomes), those of Cladiucha punctata, Cladiucha magnoliae, Megabeleses magnoliae, and Megabeleses liriodendrovorax, are newly sequenced and comparatively analyzed with previously reported tenthredinid mitogenomes. The close investigation of mitogenomes and the phylogeny of Tenthredinidae leads us to the following conclusions: The subfamilial relationships and phylogenetic placements within Tenthredinidae are mostly found to be similar to the previously suggested phylogenies. However, the present phylogeny supports the monophyly of Megabelesesinae as a subfamily, with the sister-group placement of Cladiucha and Megabeleses outside of Allantinae. The occurrence of the same type of tRNA rearrangements (MQI and ANS1ERF) in the mitogenomes of Megabelesesinae species and the presence of apomorphic morphological characters also provide robust evidence for this new subfamily. The divergence and diversification times of the subfamilies appear to be directly related to colonization of the flowering plants following the Early Cretaceous. The origin time and diversification patterns of Megabelesesinae were also well matched with the divergence times of their host plants from Magnoliaceae.

How are the mitochondrial genomes reorganized in Hexapoda? Differential evolution and the first report of convergences within Hexapoda

The evolution of the Hexapoda mitochondrial genome has been the focus of several genetic and evolutionary studies over the last decades. However, they have concentrated on certain taxonomic orders of economic or health importance. The recent increase of mitochondrial genomes sequencing of diverse taxonomic orders generates an important opportunity to clarify the evolution of this group of organisms. However, there is no comparative study that investigates the evolution of the Hexapoda mitochondrial genome. In order to verify the level of rearrangement and the mitochondrial genome evolution, we performed a comparative genomic analysis of the Hexapoda mitochondrial genome available in the NCBI database. Using a combination of bioinformatics methods to carefully examine the mitochondrial gene rearrangements in 1198 Hexapoda species belonging to 32 taxonomic orders, we determined that there is a great variation in the rate of rearrangement by gene and by taxonomic order. A higher rate of genetic reassortment is observed in Phthiraptera, Thysanoptera, Protura, and Hymenoptera; compared to other taxonomic orders. Twenty-four events of convergence in the genetic order between different taxonomic orders were determined, most of them not previously reported; which proves the great evolutionary dynamics within Hexapoda.

Comparative Mitogenomics in Hyalella (Amphipoda: Crustacea)

We present the sequencing and comparative analysis of 17 mitochondrial genomes of Nearctic and Neotropical amphipods of the genus Hyalella, most from the Andean Altiplano. The mitogenomes obtained comprised the usual 37 gene-set of the metazoan mitochondrial genome showing a gene rearrangement (a reverse transposition and a reversal) between the North and South American Hyalella mitogenomes. Hyalella mitochondrial genomes show the typical AT-richness and strong nucleotide bias among codon sites and strands of pancrustaceans. Protein-coding sequences are biased towards AT-rich codons, with a preference for leucine and serine amino acids. Numerous base changes (539) were found in tRNA stems, with 103 classified as fully compensatory, 253 hemi-compensatory and the remaining base mismatches and indels. Most compensatory Watson–Crick switches were AU -> GC linked in the same haplotype, whereas most hemi-compensatory changes resulted in wobble GU and a few AC pairs. These results suggest a pairing fitness increase in tRNAs after crossing low fitness valleys. Branch-site level models detected positive selection for several amino acid positions in up to eight mitochondrial genes, with atp6 and nad5 as the genes displaying more sites under selection.

Comparative mitogenomics and phylogenetics of the stinging wasps (Hymenoptera: Aculeata)

The stinging wasps (Hymenoptera: Aculeata) include diverse groups such as vespid wasps, ants and bees. Phylogenetic relationships among major lineages of stinging wasps have been inferred from molecular and morphological data. However, the genomic features of the mitochondrial genomes and their phylogenetic utility remain to be explored. In this study, we determined 23 mitochondrial genomes from the Aculeata. Four Mutillidae species showed relatively low A + T content compared to other species of the Aculeata (69.7%-77.4%). Eleven out of 44 species, mainly from the Chrysididae and the Pompilidae, showed reversals of GC skews. Gene rearrangements occurred across the species. Patterns of tRNA rearrangement were conserved in some groups, including the Chrysididae, Bethylidae, Pompilidae, Scolioidea and Vespoidea. Rearrangement of protein-coding genes were found in 12 out of 44 species of the Aculeata, including all four species from the Chrysididae, both species from the Bethylidae, one species from the Dryinidae, all three Scolioidea species and two Apoidea species. Phylogenetic inference showed a long branch in species with unusual genomic features, such as in the Mutillidae and Bethylidae. By excluding these species, we found paraphyly of the Chrysidoidea and a sister group relationship between the Formicoidea and Vespoidea. These results improve our understanding of the evolution of mitochondrial genomes in the Aculeata and, in general, the evolution across this subclade.

Four new mitochondrial genomes for the basal bee family Melittidae (Hymenoptera: Apoidea)

We present four new complete mitochondrial genomes for Dasypoda hirtipes, Melitta schultzei, Capicola nanula and Samba griseonigra belonging to the basally branching bee family Melittidae covering four genera in three tribes (Melittini, Hesperaspini, Dasypodaini) and two subfamilies (Melittinae, Dasypodainae). The mitogenomes vary between 15,884 and 20,324 bp in length and consist of the typical set of 13 protein-coding genes, 22 tRNAs, two rRNAs and the control region. These new mitogenomes raise the number of available mitochondrial genomes for the family Melittidae to five and will help to shed light on the phylogenetic relationships within Melittidae and their position within the Anthophila.

The complete mitochondrial genome sequence of Oryctes rhinoceros (Coleoptera: Scarabaeidae) based on long-read nanopore sequencing

BACKGROUND

The coconut rhinoceros beetle (CRB, Oryctes rhinoceros) is a severe and invasive pest of coconut and other palms throughout Asia and the Pacific. The biocontrol agent, Oryctes rhinoceros nudivirus (OrNV), has successfully suppressed O. rhinoceros populations for decades but new CRB invasions started appearing after 2007. A single-SNP variant within the mitochondrial cox1 gene is used to distinguish the recently-invading CRB-G lineage from other haplotypes, but the lack of mitogenome sequence for this species hinders further development of a molecular toolset for biosecurity and management programmes against CRB. Here we report the complete circular sequence and annotation for CRB mitogenome, generated to support such efforts.

METHODS

Sequencing data were generated using long-read Nanopore technology from genomic DNA isolated from a CRB-G female. The mitogenome was assembled with Flye v.2.5, using the short-read Illumina sequences to remove homopolymers with Pilon, and annotated with MITOS. Independently-generated transcriptome data were used to assess the O. rhinoceros mitogenome annotation and transcription. The aligned sequences of 13 protein-coding genes (PCGs) (with degenerate third codon position) from O. rhinoceros, 13 other Scarabaeidae taxa and two outgroup taxa were used for the phylogenetic reconstruction with the Maximum likelihood (ML) approach in IQ-TREE and Bayesian (BI) approach in MrBayes.

RESULTS

The complete circular mitogenome of O. rhinoceros is 20,898 bp in length, with a gene content canonical for insects (13 PCGs, two rRNA genes, and 22 tRNA genes), as well as one structural variation (rearrangement of trnQ and trnI) and a long control region (6,204 bp). Transcription was detected across all 37 genes, and interestingly, within three domains in the control region. ML and BI phylogenies had the same topology, correctly grouping O. rhinoceros with one other Dynastinae taxon, and recovering the previously reported relationship among lineages in the Scarabaeidae. In silico PCR-RFLP analysis recovered the correct fragment set that is diagnostic for the CRB-G haplogroup. These results validate the high-quality of the O. rhinoceros mitogenome sequence and annotation.

Comparative Mitogenomic Analysis of Two Cuckoo Bees (Apoidea: Anthophila: Megachilidae) with Phylogenetic Implications

Bees (Hymenoptera, Apoidea and Anthophila) are distributed worldwide and considered the primary pollinators of angiosperm. Megachilidae is one of the largest families of Anthophila. In this study, two complete mitogenomes of cuckoo bees in Megachilidae, namely Coelioxys fenestrata and Euaspis polynesia, were amplified and sequenced, with a length of 17,004 bp (C. fenestrata) and 17,682 bp (E. polynesia). The obtained results show that 37 mitogenomic genes and one putative control region were conserved within Hymenoptera. Truncated stop codon T was found in the cox3 gene of E. polynesia. The secondary structure of small (rrnS) and large (rrnL) rRNA subunits contained three domains (28 helices) and five domains (44 helices) conserved within Hymenoptera, respectively. Compared with ancestral gene order, gene rearrangement events included local inversion and gene shuffling. In order to reveal the phylogenetic position of cuckoo bees, we performed phylogenetic analysis. The results supported that all families of Anthophila were monophyletic, the tribe-level relationship of Megachilidae was Osmiini + (Anthidiini + Megachilini) and Coelioxys fenestrata was clustered to the Megachile genus, which was more closely related to Megachile sculpturalis and Megachile strupigera than Euaspis polynesia.