Rearrangement and evolution of mitochondrial genomes in Thysanoptera (Insecta)

Population Genetics, Demographic History, and Potential Distributions of the New Important Pests Monolepta signata (Coleoptera: Chrysomelidae) on Corn in China

Monolepta signata are polyphagous pest widely distributed in China, and the damage as well as economic losses it caused were increasing in recent years. Knowledge of species diversity, population structure and habitat suitability could enhance the efforts of pest control. Here, we sampled the populations of M. signata in almost all of China's major corn-producing regions. A total of 568 sequences were obtained from each gene. There were 48, 29, and 30 haplotypes of COI, ITS2 and EF-1α, respectively. The genetic distance between the HuangHuaiHai population and other populations was the largest. There were 61.90%, 71.43% and 61.90% of Nm values smaller than 1 in COI, ITS2 and EF-1α, respectively, which indicated that gene flow between most populations was weak. The degree of differentiation in most populations of M. signata was relatively high. The population of M. signata has also experienced rapid expansion. Population history dynamic analysis showed that the effective population size of M. signata remained relatively stable before 0.075 Ma. There was a slow contraction trend from 0.075 to 0.010 Ma. It has been rapidly and continuously expanding since 0.010 Ma. Among the investigated geographical populations, the "yellow-spot type" was only present in the populations of southern and southwestern regions, while the "two-spot type" and "four-spot type" were widely distributed in all other geographical populations. Predictions of the potential distribution areas of M. signata indicated that the northeast and north China regions will remain being the high suitability areas of M. signata in the future. Our results will not only facilitate studies on the phylogeography of M. signata but also benefit the effective monitoring and management of this agricultural pest.

Comparative Mitogenomes and Phylogenetic Analyses of Coccinellidae (Coleoptera: Coccinelloidea)

Coccinellidae (ladybird beetles) comprises around 6900 described species with a worldwide distribution and exhibits a broad trophic diversity. Complete mitochondrial genomes (mitogenomes) are valuable resources in many research fields, such as genomics, population genetics, molecular evolution, and phylogenetics. Here we sequenced and report the complete mitogenome of Calvia chinensis, Micraspis discolor, Harmonia eucharis, and Oenopia kirbyi. By comparing with the 36 complete mitogenomes published in GenBank, we found that the long noncoding region (LNCR) between trnI and trnQ is present in the mitogenome of Chilocorini and Coccinellini, and the size of LNCR is positively correlated with their mitogenome size. The variable number tandem repeat (VNTR) was detected in the LNCR of Calvia chinensis and Oenopia kirbyi, indicating that the LNCR may be associated with the transcriptional regulation of the mitogenome. Heterogeneity in the base composition was encountered among the mitogenomes in Coccinellidae, especially in Noviini and some species of Epilachnini and Coccinellini, which may lead to unstable phylogenetic topologies. Phylogenetic relationships have been reconstructed by maximum likelihood and Bayesian inferences based on two mitogenomic datasets, PCG_rRNA (all 13 PCGs and two rRNAs) and PCG12_rRNA (all 13 PCGs with the third codon position excluded and two rRNAs). Our results are close to the subfamily and tribe classification system reported in previous studies and suggest the maximum likelihood analysis based on the PCG12_rRNA dataset is more sensitive in avoiding the false grouping of unrelated taxa with similar base composition in the reconstruction of the phylogeny.

Unravelling the complete mitochondrial genomes of Thrips tabaci Lindeman and Thrips parvispinus Karny (Thysanoptera: Thripidae) and their phylogenetic implications

Onion (Allium cepa Linnaeus) is an important vegetable crop valued for its nutritional properties and economics worldwide. Onion cultivation faces serious threats from pests and diseases, particularly onion thrips (Thrips tabaci), which cause substantial yield losses. Recently, Black thrips (Thrips parvispinus), an invasive key pest of chili, have been reported to cause severe damage in onion crop and is likely to devastate the onion cultivation in near future. Therefore, this study was conducted to address the knowledge gap concerning the genetic basis and evolutionary history of T. tabaci and T. parvispinus through sequencing of their mitochondrial genomes. T. tabaci and T. parvispinus were collected from different locations in Maharashtra, India, and reared in the laboratory. The mitochondrial genomes of T. tabaci and T. parvispinus were sequenced to a length of 15,277 and 15,285 bp, respectively. Both genomes exhibited similar gene organization with regard to thirteen protein-coding genes and two rRNA genes. T. tabaci contained 19 tRNA genes whereas T. parvispinus contained 18 tRNA genes. The evolutionary positions of T. tabaci and T. parvispinus within the Thysanoptera order were elucidated through phylogenetic analysis of the mitogenomes of 15 thrips species. These findings provide crucial insights into the genetic makeup and evolutionary dynamics of both the thrips species, thereby aiding the development of novel and sustainable pest management strategies to mitigate their impacts on crops in the changing climate scenario.

Insect Mitochondrial Genomics: A Decade of Progress

The past decade has seen the availability of insect genomic data explode, with mitochondrial (mt) genome data seeing the greatest growth. The widespread adoption of next-generation sequencing has solved many earlier methodological limitations, allowing the routine sequencing of whole mt genomes, including from degraded or museum specimens and in parallel to nuclear genomic projects. The diversity of available taxa now allows finer-scale comparisons between mt and nuclear phylogenomic analyses; high levels of congruence have been found for most orders, with some significant exceptions (e.g., Odonata, Mantodea, Diptera). The evolution of mt gene rearrangements and their association with haplodiploidy have been tested with expanded taxonomic sampling, and earlier proposed trends have been largely supported. Multiple model systems have been developed based on findings unique to insects, including mt genome fragmentation (lice and relatives) and control region duplication (thrips), allowing testing of hypothesized evolutionary drivers of these aberrant genomic phenomena. Finally, emerging research topics consider the contributions of mt genomes to insect speciation and habitat adaption, with very broad potential impacts. Integration between insect mt genomic research and other fields within entomology continues to be our field's greatest opportunity and challenge.

Phylogenetic and Comparative Genomics Study of Cephalopina titillator Based on Mitochondrial Genomes

Camel bot fly (Cephalopina titillator) larvae cause myiasis in domesticated and wild camels, resulting in significant economic losses to the camel industry and posing a serious global public health concern. To date, only one mitochondrial genome (mitogenome) of C. titillator isolated from the Alxa Bactrian camel has been reported. Herein, C. titillator was isolated from the Junggar Bactrian camel to assemble a complete circular mitogenome with a length of 16,552 bp encoding 13 protein-coding genes, 22 tRNA genes, and two rRNA genes. The mitogenome showed a high A + T content (73.31%), positive AT-skew (0.12), and negative GC-skew (-0.34) base composition patterns. All protein-coding genes (PCGs) employed ATG, ATA, ATT, GTG, or TCG as the start codons and TAA, TAG, or single T as the stop codons. Similar to other parasites in the Oestridae subfamily, the mitogenome was structurally conserved, with genes retaining the same order and direction as those in the ancestral insect mitogenome. The phylogenetic analysis clustered this species with the Oestrinae, showing that the subfamily did not exhibit monophyly. C. titillator isolated from the Junggar Bactrian camel was found to be a sister lineage to that isolated from the Alxa Bactrian camel. Despite the lack of data on the mitogenome of C. titillator isolated from dromedaries in the Middle East, phylogenetic analysis of C. titillator isolated from Xinjiang revealed one distinct lineage of the Xinjiang camel nasal bot fly. In conclusion, this study reports the complete mitogenome of Xinjiang C. titillator for the first time, providing valuable data for future studies on the phylogenetic relationships in this subfamily.

The Complete Mitochondrial Genome of Chilo infuscatellus (Lepidoptera: Pyralidae), and Related Phylogenetic Analysis

The Chilo infuscatellus (Lepidoptera: Pyralidae) is a significant pest of sugarcane in China. The genome-level characteristics of this pest are important genetic resources for identification, phylogenetic analysis, and even management. In the present study, the complete mitogenome of C. infuscatellus was sequenced and characterized. The assembled mitochondrial genome is 15,252 bp in length and includes 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and an A + T-rich region. Except for the CGA codon for the cox1 gene, the PCGs are initiated with ATN codons (ATG, ATT, and ATA). These PCGs are terminated with TAA or an incomplete termination codon of a single T. Except for the loss of the "DHU" arm for trnS1, the tRNA genes were folded into the typical cloverleaf structure. The A + T-rich region has a high AT content of 96.19% and contains the motifs "ATAGA" and "ATTTA", as well as a 19 bp poly-T stretch and microsatellite regions. The C. infuscatellus mitogenome exhibits a conserved gene order among lepidopteran insects, with a rearrangement of the trnM gene compared to the ancestral insect gene order. Phylogenetic analysis based on the 13 PCGs using Bayesian inference (BI) and maximum likelihood (ML) methods confirmed the monophyly of Pyralidae and Crambidae within Pyraloidea. The relationships between subfamilies in Pyralidae can be described as (Galleriinae + (Phycitinae + (Pyralinae + Epipaschiinae))). The "PS clade" and "non-PS clade" were formed within the family Crambidae. These findings provide valuable genetic resources for the identification, phylogenetic analysis, and management of sugarcane borers, contributing significantly to our understanding of the phylogeny of Pyraloidea insects and their evolution.

Mitogenomics of the zoonotic parasite Echinostoma miyagawai and insights into the evolution of tandem repeat regions within the mitochondrial non-coding control region

Echinostoma miyagawai is a cosmopolitan parasite within the Echinostomatidae and is a cause of human echinostomiasis. Species within the family have been a challenge to disentangle with E. miyagawai being synonyms of several other Echinostoma species. However, complete mitochondrial genomes have been shown to be vital in distinguishing echinostomatid species, but detailed comparisons of not only gene content but also structural features have been limited. Using long range sequencing techniques, the complete mitochondrial genome of E. miyagawai was sequenced and compared to other members of Echinostomatidae. In total 12 protein coding genes, 2 ribosomal RNA genes and 22 transfer RNA genes were identified, as was an extensive noncoding control region (CR), consisting of 2 types of multiple tandem repeat units. Phylogenetic analyses of complete mitochondrial genomes corresponded to previous studies on single mitochondrial genes and nuclear ribosomal nuclear markers confirmed E. miyagawai to be within in the ‘Echinostoma revolutum’ group. The tandem repeat units found in the CR contained promoter sequences containing domains typical of initiation sites for replication and transcription as well as several palindromic regions which were shared between echinostomatid species. The study illustrates not only the utility complete mitogenomes in disentangling the relationship between these parasite species, but also provides some insight into the potential adaptations and other evolutionary processes that may govern the divergence of mitochondrial genomes for the first time in echinostomatids.

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.

Rapid evolution of mitochondrion-related genes in haplodiploid arthropods

BACKGROUND

Mitochondrial genes and nuclear genes cooperate closely to maintain the functions of mitochondria, especially in the oxidative phosphorylation (OXPHOS) pathway. However, mitochondrial genes among arthropod lineages have dramatic evolutionary rate differences. Haplodiploid arthropods often show fast-evolving mitochondrial genes. One hypothesis predicts that the small effective population size of haplodiploid species could enhance the effect of genetic drift leading to higher substitution rates in mitochondrial and nuclear genes. Alternatively, positive selection or compensatory changes in nuclear OXPHOS genes could lead to the fast-evolving mitochondrial genes. However, due to the limited number of arthropod genomes, the rates of evolution for nuclear genes in haplodiploid species, besides hymenopterans, are largely unknown. To test these hypotheses, we used data from 76 arthropod genomes, including 5 independently evolved haplodiploid lineages, to estimate the evolutionary rates and patterns of gene family turnover of mitochondrial and nuclear genes.

RESULTS

We show that five haplodiploid lineages tested here have fast-evolving mitochondrial genes and fast-evolving nuclear genes related to mitochondrial functions, while nuclear genes not related to mitochondrion showed no significant evolutionary rate differences. Among hymenopterans, bees and ants show faster rates of molecular evolution in mitochondrial genes and mitochondrion-related nuclear genes than sawflies and wasps. With genome data, we also find gene family expansions and contractions in mitochondrion-related genes of bees and ants.

CONCLUSIONS

Our results reject the small population size hypothesis in haplodiploid species. A combination of positive selection and compensatory changes could lead to the observed patterns in haplodiploid species. The elevated evolutionary rates in OXPHOS complex 2 genes of bees and ants suggest a unique evolutionary history of social hymenopterans.

Mitochondrial Genome Characteristics and Phylogenetic Analysis of Fulmekiola serrata (Kobus) (Thysanoptera: Thripidae)

Sugarcane thrips, Fulmekiola serrata (Kobus) (Thysanoptera: Thripidae), is a common foliar pest that infests sugarcane and is found throughout tropical and subtropical countries. In this study, we obtained and analyzed the complete mitochondrial genome of F. serrata for the first time and explored the phylogenetic relationships of the higher-order elements of Thysanoptera members at the mitochondrial level. The complete mitochondrial genome of F. serrata is 16,596 bp in length and includes 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and 1 noncoding control region. A+T accounted for 75% of the total bases in the mitochondrial genome of F. serrata, revealing an obvious AT bias. Among the 13 PCGs, except for nad5, which had a start codon of TTG, the remaining genes had ATNs typical of insects (ATA, ATT, ATC, and ATG); nad1, nad2, nad3, and atp8 had incomplete termination codons of TA or T. The remaining nine PCGs were complete with the termination codon TAA. Of the 22 tRNA secondary structures, all were typical cloverleaf secondary structures except for trnS1, which was missing the DHU arm. Compared with the hypothetical ancestral gene arrangement of arthropods, F. serrata presented extensive gene rearrangement, with 23 translocated genes, 8 inverted genes, and 5 shuffled genes. Both maximum likelihood (ML) and Bayesian inference (BI) phylogenetic trees resulted in similar topologies: ((Thripidae + (Stenurothripidae + Aeolothripidae)) + Phlaeothripidae), with Thripidae, Aeolothripidae and Phlaeothripidae being monophyletic groups, whereas F. serrata is closely related to Thrips palmi, and the two are sister groups.

Mitogenome-Based Phylogeny with Divergence Time Estimates Revealed the Presence of Cryptic Species within Heptageniidae (Insecta, Ephemeroptera)

Heptageniidae are known for their flat heads and bodies and are divided into three subfamilies. Despite the extensive diversity within this group and considerable efforts made to understand their evolutionary history, the internal classifications and origin time of Heptageniidae remains controversial. In this study, we newly sequenced 17 complete mitogenomes of Heptageniidae to reconstruct their phylogenetic positions within this family. Because of the ambiguous time of origin, our study also estimated the divergence time within Heptageniidae based on five fossil calibration points. The results of BI and ML trees all highly supported the monophyly of Heptageniidae and three subfamilies. The phylogenetic relationship of Rhithrogeninae + (Ecdyonurinae + Heptageniinae) was also recovered. The divergence time showed that Heptageniidae originated from 164.38 Mya (95% HPD, 150.23-181.53 Mya) in the mid-Jurassic, and Rhithrogeninae originated from 95.54 Mya (95% HPD, 73.86-120.19 Mya) in the mid-Cretaceous. Ecdyonurinae and Heptageniinae began to diverge at 90.08 Mya (95% HPD, 68.81-113.16 Mya) in the middle Cretaceous. After morphological identification, analysis of the mitogenome's composition, genetic distance calculation, phylogenetic analysis, and divergence time calculation, we suggest that two different populations of Epeorus montanus collected from Aksu, Xinjiang Uygur Autonomous Region (40°16' N, 80°26' E) and Xinyuan, Xinjiang Uygur Autonomous Region (43°20' N, 83°43' E) in China are cryptic species of E. montanus, but further detailed information on their morphological characteristics is needed to fully identify them.

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.

Phylogenetic and Comparative Analysis of Cryptochironomus, Demicryptochironomus and Harnischia Inferred from Mitogenomes (Diptera: Chironomidae)

(1) Background: Mitochondrial genomes have been extensively employed as a crucial marker in numerous dipteran families for understanding phylogenetics and systematics relations, thereby playing a pivotal role in molecular biology studies. The phylogenetic relationship of the Harnischia generic complex remains contentious due to the paucity of taxonomic and molecular data. Specifically, the evolutionary relationships among Cryptochironomus, Demicryptochironomus, and Harnischia are still unclear. (2) Methods: In this study, Polypedilum and Endochironomus were used as outgroups to analyze phylogenetic relationships among Cryptochironomus, Demicryptochironomus, and Harnischia, mitogenomes of four Cryptochironomus, two Demicryptochironomus, two Harnischia, and two Cladopelma were newly sequenced. Subsequently, we conducted a thorough analysis of the nucleotide composition, sequence length, and evolutionary rate. (3) Results: All mitogenomes exhibited structural conservation, with all genes consistently arranged in the identical order as that of the ancestral mitogenome. Nucleotide composition varied significantly among different genes, and the control region displayed the highest A + T content. All protein-coding genes undergo rigorous purification selection, with the ATP8 gene exhibiting the most rapid evolutionary rate among them. Utilizing Bayesian Inference (BI) and Maximum Likelihood (ML) methods across various databases, we reconstructed the phylogenetic relationships among the genera within the Harnischia generic complex, drawing insights from an analysis of 14 mitochondrial genomes. (4) Conclusions: Our results showed that the monophyly of the genera Harnischia was well supported in all topologies; Cryptochironomus is sister to Demicryptochironomus.

Characterizing the complete mitogenome of Odontothrips phaseoli (Thysanoptera: Thripidae) and its mitochondrial phylogeny

Described originally from Heilongjiang, China, Odontothrips phaseoli is a potential pest of threatening bean plant in northern China. The complete mitochondrial genome of O. phaseoli was sequenced and assembled, with a total length of 15,540 bp. Within this genome, 37 genes have been identified: 13 PCGs, 22 tRNAs, two rRNAs, and two putative control regions. Most PCGs terminate with TAA, while four genes (atp8, nad1, nad2 and nad4) use an incomplete 'T' and nad6 employs TAG as the stop codon. Compared to the mitogenome of the ancestral insect, O. phaseoli displays significant gene rearrangement. However, it retains three conserved gene blocks in common with its related species, Megalurothrips usitatus, both of which belong to the Megalurothrips genus-group. The phylogenetic tree, constructed based on the entire mitogenome dataset of all thrips species available in NCBI, shows that the two species cluster closely together. This alignment might underscore the close link between gene arrangements and the phylogeny relationships.

Purifying selection drove the adaptation of mitochondrial genes along with correlation of gene rearrangements and evolutionary rates in two subfamilies of Whitefly (Insecta: Hemiptera)

Insect mitochondrial genomes (mitogenomes) are usually represented by a conserved gene order. Whiteflies exhibit gene rearrangement in their mitogenomes; however, understanding how nucleotide substitution rates shape gene rearrangement in whiteflies is unclear due to the limited number of mitogenomes. Additionally, the mechanisms by which selection pressure drives adaptations in mitochondrial genes in the two subfamilies of whiteflies are not yet known. Here, we analyzed 18 whitefly mitogenomes, including one newly generated mitogenome, to compare nucleotide substitution rates, selection pressure, and gene arrangements. The newly generated mitogenome is reported along with reannotation of Pealius mori and comparisons to other whitefly mitogenomes. Comparative studies on nucleotide composition of 18 whiteflies revealed the positive GC skewness, confirming the reversal of strand asymmetry. We found 11 rearranged gene orders within two subfamilies of whiteflies with 8-18 breakpoints of gene rearrangements. Members of the subfamily Aleyrodinae exhibit more complex pathways in the evolution of gene order as compared to the subfamily Aleurodicinae. Our findings also revealed that the increase or reduction of nucleotide substitution rates does not have an impact on any of the gene rearrangement scenarios depicting neutral correlation. Selection pressure analysis revealed that the mitogenomes from members of both the subfamilies Aleurodicinae and Aleyrodinae are characterized by intense purifying selection pressure.

Evolutionary and phylogenetic insights from the mitochondrial genomic analysis of Diceraeus melacanthus and D. furcatus (Hemiptera: Pentatomidae)

The mitochondrial genomes of D. melacanthus and D. furcatus were sequenced and used to investigate the phylogenetic relationships with 54 species of Pentatomidae. Their mitogenomes are 17,197 and 15,444 bp-long, respectively, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22/21 transfer RNA genes, with conserved gene arrangement. Leu, Lys, and Ser were the most common amino acids in their PCGs. PCGs evolutionary analysis indicated their mitogenomes are under purifying selection, and the most conserved genes are from the cytochrome complex, reinforcing their suitability as markers for molecular taxonomy. We identified 490 mtSSRs in 56 Pentatomidae species, with large variation and a positive correlation between mtSSR number and genome size. Three mtSSRs were identified in each Diceraeus species. Only the mtSSR in the nad6 (D. melacanthus) and nad4 (D. furcatus) appear to have application as molecular markers for species characterization. Phylogenetic analysis confirmed the monophyly of Pentatomidae. However, our analysis challenged the monophyly of Pentatominae and Podopinae. We also detected unexpected relationships among some tribes and genera, highlighting the complexity of the internal taxonomic structure of Pentatomidae. Both Diceraeus species were grouped in the same clade with the remaining Carpocorini analyzed.

Characterization of the complete mitochondrial genome of an endemic species in China, Aulocera merlina (Lepidoptera: Nymphalidae: Satyrinae) and phylogenetic analysis within Satyrinae

The mitochondrial genome (mitogenome) has been extensively used as molecular markers in determining the insect phylogenetic relationships. In order to resolve the relationships among tribes and subtribes of Satyrinae at the mitochondrial genomic level, we obtained the complete mitogenome of Aulocera merlina (Oberthür, 1890) (Lepidoptera: Nymphalidae: Satyrinae) with a size of 15,259 bp. The mitogenome consisted of 37 typical genes, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an A + T-rich region. The gene organization and arrangement were similar to those of all other known Satyrinae mitogenomes. All PCGs were initiated with the canonical codon pattern ATN, except for the cox1 gene, which used an atypical CGA codon. Nine PCGs used the complete stop codon TAA, while the remaining PCGs (cox1, cox2, nad4, and nad5) were terminated with a single T nucleotide. The canonical cloverleaf secondary structures were found in all tRNAs, except for trnS1 which lacked a dihydrouridine arm. The 448 bp A + T-rich region was located between rrnS and trnM, and it included the motif ATAGA followed by a 19-bp poly-T stretch and a microsatellite-like (TA)6 element preceded by the ATTTA motif. The phylogenetic tree, inferred using Bayesian inference and maximum likelihood methods, generated similar tree topologies, revealing well-supported monophyletic groups at the tribe level and recovering the relationship ((Satyrini + Melanitini) + ((Amathusiini + Elymniini) + Zetherini)). The close relationship between Satyrina and Melanargiina within the Satyrini was widely accepted. Additionally, Lethina, Parargina, and Mycalesina were closely related and collectively formed a sister group to Coenonymphina. Moreover, A. merlina was closely related to Oeneis buddha within the Satyrina. These findings will provide valuable information for future studies aiming to elucidate the phylogenetic relationships of Satyrinae.

The complete mitochondrial genome of Tauraco livingstonii (Musophagidae: Tauraco)

Livingstone's turaco, Tauraco livingstonii, belongs to the family Musophagidae. In this study, we obtained the complete mitochondrial genome sequence of Livingstone's turaco by high-throughput sequencing technology and constructed a phylogenetic tree. It was found that the mitochondria of this species are 19,015 bp in length and contain a total of 37 genes, comprising 13 protein-coding genes, 22 tRNA genes, and 2 rRNA genes. The base composition of the mitochondrial genome is 31.61% A, 24.22% T, 30.64% C, and 13.52% G, with a GC content of 44%. Notably, an intriguing phenomenon of mitochondrial genome rearrangements was observed, characterized by the duplication of the tRNA Glu-L-CR gene order. In addition, the results of the phylogenetic tree analysis shed light on the taxonomic position of Livingstone's turaco and supported the taxonomy of Otidimorphae. The study provides a basis for future phylogenetic and taxonomic investigations of Musophagiformes.

The First Complete Mitochondrial Genome of the Genus Litostrophus: Insights into the Rearrangement and Evolution of Mitochondrial Genomes in Diplopoda

This study presents the complete mitochondrial genome (mitogenome) of Litostrophus scaber, which is the first mitogenome of the genus Litostrophus. The mitogenome is a circular molecule with a length of 15,081 bp. The proportion of adenine and thymine (A + T) was 69.25%. The gene ND4L used TGA as the initiation codon, while the other PCGs utilized ATN (A, T, G, C) as the initiation codons. More than half of the PCGs used T as an incomplete termination codon. The transcription direction of the L. scaber mitogenome matched Spirobolus bungii, in contrast to most millipedes. Novel rearrangements were found in the L. scaber mitogenome: trnQ -trnC and trnL1- trnP underwent short-distance translocations and the gene block rrnS-rrnL-ND1 moved to a position between ND4 and ND5, resulting in the formation of a novel gene order. The phylogenetic analysis showed that L. scaber is most closely related to S. bungii, followed by Narceus magnum. These findings enhance our understanding of the rearrangement and evolution of Diplopoda mitogenomes.

Variability and evolution of gene order rearrangement in mitochondrial genomes of arthropods (except Hexapoda)

In the species-rich Phylum Arthropoda, the mitochondrial genome is relatively well conserved both in terms of number and order of genes. However, specific clades have a 'typical' gene order that differs from the putative arthropod ancestral arrangement. The aim of this work was to compare the rate of mitochondrial gene rearrangements at inter- and intra-taxonomic levels in the Arthropoda and to postulate the most parsimonious ancestral orders representing the four major arthropod lineages. For this purpose, we performed a comparative genomic analysis of arthropod mitochondrial genomes available in the NCBI database. Using a combination of bioinformatics methods that examined mitochondrial gene rearrangements in 464 species of arthropods from three subphyla (Chelicerata, Myriapoda, and Crustacea [except Hexapoda, previously analyzed]), we observed differences in the rate of rearrangement within major lineages. A higher rate of mitochondrial genome rearrangement was observed in Crustacea and Chelicerata compared to Myriapoda. Likewise, early branching clades exhibit less variability in mitochondrial genome order than late branching clades, within each subphylum. We identified 'hot regions' in the mitochondrial genome of each studied subphylum, and postulated the most likely ancestral gene order in each subphylum and taxonomic order. Our work provides new evidence on the evolutionary dynamics of mitochondrial genome gene order in arthropods and new mitochondrial genome architectures in different taxonomic divisions within each major lineage of arthropods.

Multiple independent origins of duplicated mitochondrial control regions indicate an apomorphy in the Thysanoptera (Insecta)

The mitochondrial genome (mitogenome) of thrips is characterized by the presence of control region (CR) duplication. However, the evolution pattern of duplicated CRs in thrips is still unclear. In this study, the multiple independent origins of duplicated CR indicated that the CR duplication was not an ancestral state for Thysanoptera. The macroevolutionary pattern suggested that the earliest CR duplication event occurred in the middle Cretaceous (94.85 Ma) coincided with rearrangement events forming the ancestors of Aeolothripidae, but much later than that forming the ancestors of the suborder Terebrantia. The mitogenome with duplicated CRs showed a higher rate of gene rearrangement. The sequence similarity of the CR copies and divergence time were negatively correlated, indicating age-related deterioration of mitochondrial function. No significant differences were found in the mitochondrial DNA, the P123 and P4FD between the single and multiple-CR charactered mitogenomes, which suggested that the duplicated CRs may not affect the replication process in thrip mitogenome. The mitogenomes with duplicated CRs (mean: 0.0088 subs/s/my) show a significantly increased evolutionary rate than that with a single one (mean: 0.0058 subs/s/my). However, it seems that this higher evolutionary rate did not have adaptive mechanisms in Terebrantia. We speculated that the duplicated CRs may cause a more intense production of energy by mitochondria, and an accelerated mutation and substitution rate is expected in such mitogenomes. Our study provided new insights into the presence of CR duplications and their evolution in the mitogenomes of thrips.

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.

Revisiting mitogenome evolution in Medusozoa with eight new mitochondrial genomes

Mitogenomics has improved our understanding of medusozoan phylogeny. However, sequenced medusozoan mitogenomes remain scarce, and Medusozoa phylogeny studies often analyze mitogenomic sequences without incorporating mitogenome rearrangements. To better understand medusozoan evolution, we analyzed Medusozoa mitogenome phylogeny by sequencing and assembling eight mitogenomes from three classes (Cubozoa, Hydrozoa, and Scyphozoa). We reconstructed the mitogenome phylogeny using these mitogenomes and 84 other existing cnidarian mitogenomes to study mitochondrial gene rearrangements. All reconstructed mitogenomes had 13 mitochondrial protein-coding genes and two ribosomal genes typical for Medusozoa. Non-cubozoan mitogenomes were all linear and had typical gene orders, while arrangement of genes in the fragmented Cubozoa (Morbakka sp.) mitogenome differed from other Cubozoa mitogenomes. Gene order comparisons and ancestral state reconstruction suggest minimal rearrangements within medusozoan classes except for Hydrozoa. Our findings support a staurozoan ancestral medusozoan gene order, expand the pool of available medusozoan mitogenomes, and enhance our understanding of medusozoan phylogenetic relationships.

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.

Characterization of the complete mitochondrial genome and phylogenetic analysis of bean thrips Megalurothrips usitatus (Bagnall, 1913) (Thysanoptera: thripidae)

Megalurothrips usitatus (Bagnall, 1913) (Thysanoptera: Thripidae) is a widely distributed pest in Asia that primarily affects the production of snap beans and cowpea. The complete mitochondrial genome of Megalurothrips usitatus has been sequenced and annotated in this study, which is 17,209 bp long and contains 13 protein-coding genes (PCGs), two rRNAs, and 22 tRNA genes. Most of the protein-coding genes (PCGs) start with ATG except ND4 using TTG. Meanwhile, eight PCGs stop with TAA, four PCGs have an incomplete stop codon, and the gene Cytb ends with TAG. Phylogenetic analysis showed that M. usitatus is closely related to Frankliniella intonsa and F. occidentalis, providing a basis for the study of the mitochondrial evolution of Thripinae.

Comparative analysis of the two suborders of Thysanoptera and characterization of the complete mitochondrial genome of Thrips parvispinus

Thrips parvispinus is a serious sucking pest on a number of economically important crops in the oriental region. It has gained importance recently for its drastic range extension distribution as an invasive pest. Here, the complete mitochondrial genome (15,067 bp) of Thrips parvispinus was sequenced and characterized. It possesses 37 genes and the putative noncoding region is duplicated. Comparative analyses of nucleotide diversity, skewness, codon usage bias, and selection pressure in mitochondrial protein-coding genes of the available 31 thrips mitogenomes (24 Terebrantia + 7 Tubulifera) were performed. Phylogenetic analysis showed a sister relationship of T. parvispinus to the clade (T. florum + T. hawaiiensis). Phylogenetic analyses formed the monophyly of subfamilies Phlaeothripinae and Idolothripinae within the family Phlaeothripidae (Suborder Tubulifera). Low nucleotide diversity was indicative of reversal of strand asymmetry in the Tubulifera. Neutrality analysis showed that directional mutation plays a major role in shaping codon usage bias in both suborders. Principal component analysis indicated distinct codon usage patterns in each suborder. Our data suggested weaker selection constrains on Terebrantia than in the Tubulifera. More tubuliferan mitogenomes are required to resolve previous classification hypotheses and elucidate genome evolution in these two suborders.

Phylogenetic Analysis of Two New Mitochondrial Genomes of Singapora shinshana and Seriana bacilla from the Karst Region of Southwest China

Leafhoppers have been identified as a serious threat to different plants. To explore the characteristics of mitogenomes and reveal the phylogenetic positions of two species in the Typhlocybinae, complete mitogenomes of Singapora shinshana and Seriana bacilla were sequenced and annotated for the first time with lengths of 15,402 bp and 15,383 bp, respectively. The two mitogenomes contained 13 PCGs, 22 tRNA genes and 2 rRNA genes. The genome content, gene order, nucleotide composition, codon usage and amino acid composition are similar to those of other typical mitogenomes of Typhlocybinae. All 13 PCGs started with ATN codons, except for atp8 (TTA) and nad5 (TTG). All tRNAs were folded into a typical cloverleaf secondary structure, except for tRNA-Ser1 and tRNA-Val. Moreover, phylogenetic trees were constructed and analyzed based on all the PCGs from 42 mitogenomes using maximum likelihood (ML) and Bayesian inference (BI) methods. The results supported that eleven subfamilies are all monophyletic groups, S. shinshana and S. bacilla are members of Erythroneurini, but S. shinshana and the genus Empoascanara have a very close relationship with ((((Empoascanara sipra+ Empoascanara wengangensis) + Empoascanara dwalata) + Empoascanara gracilis) + S. shinshana), and S. bacilla is closely related to the genus Mitjaevia ((Mitjaevia dworakowskae + Mitjaevia shibingensis) + S. bacilla). These results provide valuable information for future study of evolutionary relationships in Typhlocybinae.

Nine Mitochondrial Genomes of Phasmatodea with Two Novel Mitochondrial Gene Rearrangements and Phylogeny

The classification of stick and leaf insects (Order Phasmatodea) is flawed at various taxonomic ranks due to a lack of robust phylogenetic relationships and convergent morphological characteristics. In this study, we sequenced nine new mitogenomes that ranged from 15,011 bp to 17,761 bp in length. In the mitogenome of Carausis sp., we found a translocation of trnR and trnA, which can be explained by the tandem duplication/random loss (TDRL) model. In the Stheneboea repudiosa Brunner von Wattenwyl, 1907, a novel mitochondrial structure of 12S rRNA-CR1-trnI-CR2-trnQ-trnM was found for the first time in Phasmatodea. Due to the low homology of CR1 and CR2, we hypothesized that trnI was inverted through recombination and then translocated into the middle of the control region. Control region repeats were frequently detected in the newly sequenced mitogenomes. To explore phylogenetic relationships in Phasmatodea, mtPCGs from 56 Phasmatodean species (composed of 9 stick insects from this study, 31 GenBank data, and 16 data derived from transcriptome splicing) were used for Bayesian inference (BI), and maximum likelihood (ML) analyses. Both analyses supported the monophyly of Lonchodinae and Necrosciinae, but Lonchodidae was polyphyletic. Phasmatidae was monophyletic, and Clitumninae was paraphyletic. Phyllidae was located at the base of Neophasmatodea and formed a sister group with the remaining Neophasmatodea. Bacillidae and Pseudophasmatidae were recovered as a sister group. Heteroptergidae was monophyletic, and the Heteropteryginae sister to the clade (Obriminae + Dataminae) was supported by BI analysis and ML analysis.

The Complete Mitochondrial Genome of the Chinese White Wax Scale Insect, Ericerus pela Chavannes (Hemiptera: Coccidae), with Novel Gene Arrangement and Truncated tRNA Genes

The Chinese white wax scale insect, Ericerus pela Chavannes (Hemiptera: Coccidae), is one of the scale insects with great economic value and has been dispersed and reared in China for over one thousand years. Its mitochondrial genome provides essential information for the molecular identification and genetic study of this species. We assembled the complete mitochondrial genome of E. pela based on PacBio sequencing and analyzed its genomic features. The genome was 17,766 bp in length with 13 protein-coding genes, 22 tRNAs, and two rRNA genes. The analysis results showed E. pela had significant gene rearrangements involving tRNAs compared with other Coccoidea species. Furthermore, E. pela's nine tRNAs were identified to have obvious truncated structures. The phylogenetic tree compiled of the species showed a long branch of the Coccoidea lineage, which indicated the high evolutionary rate in this group. Our study revealed the mitochondrial characteristics of E. pela and enriched the mitochondrial genetic information on Coccoidea species. It also determined the occurrence of gene rearrangement for the species in this superfamily.

Mitochondrial genome rearrangements and phylogenomics of the Hymenoptera (Insecta) using an expanded taxon sample

The order Hymenoptera is one of the most species-rich insect orders, with more than 150,000 described extant species. Many hymenopteran insects have very different mitochondrial genome (mitogenome) organizations compared to the putative ancestral organization of insects. In this study, we sequenced 18 mitogenomes of representatives in the order Hymenoptera to increase taxonomic sampling. A total of 475 species were used in phylogenetic analyses, including 18 new mitogenomes and 457 existing mitogenomes. Using a site-heterogeneous model, Bayesian's inference from amino acid data yielded more resolved relationships among Hymenoptera than maximum-likelihood analysis and coalescent-based species analyses. The monophyly of Symphyta was not supported. The Xyeloidea was the earliest branching clade in the Hymenoptera. The Orussoidea was closely related to Apocrita. Within Apocrita, the Parasitoida was non-monophyletic. The monophyly of most Parasitoida superfamilies received strong support. The Proctotrupomorpha clade was supported in Bayesian's analysis. The Apoidea was monophyletic when excluding Ampulex compressa from consideration. The superfamilies Vespoidea and Chrysidoidea were found to be non-monophyletic. Comparisons of mitochondrial gene order revealed a higher frequency of gene rearrangement among lineages with a parasitoid lifestyle, particularly prominent in Chalcidoidea. The degree of gene rearrangement ranked second in specific taxa of Cynipoidea and Ichneumonoidea.

Gene rearrangements in the mitochondrial genome of ten ascaris species and phylogenetic implications for Ascaridoidea and Heterakoidea families

The Ascaridoidea family and Heterakoidea family are the most common and typical representative of large parasites. Although our understanding of these parasites' diversity has expanded by analyses of some mitochondrial genes, there is limited information on these species' evolutionary rates. Here we determined ten complete mitogenome sequences of five subfamilies of Ascaridoidea and one subfamily of Heterakoidea. The phylogenetic tree divided the Ascaridoidea into six monophyletic major clades, and the divergence time of Heterakoidea family and Ascaridoidea family can be placed during the early Carboniferous Period (300-360 Mya). The reconstruction of the ancestral state showed that the gene orders of all species in Ascaridoidea were conserved, and the Heterakoidea had obvious genome rearrangement. The conserved blocks between them were divided into five and the main types are tandem-duplication/random loss (TDRL). These results will help to better understand the gene rearrangements and evolutionary position of ascaris species.

Comparative Mitogenomic Analyses of Hydropsychidae Revealing the Novel Rearrangement of Protein-Coding Gene and tRNA (Trichoptera: Annulipalpia)

Gene rearrangement of the mitochondrial genome of insects, especially the rearrangement of protein-coding genes, has long been a hot topic for entomologists. Although mitochondrial gene rearrangement is common within Annulipalpia, protein-coding gene rearrangement is relatively rare. As the largest family in Annulipalpia, the available mitogenomes from Hydropsychidae Curtis, 1835 are scarce, and thus restrict our interpretation of the mitogenome characteristic. In this study, we obtained 19 novel mitogenomes of Hydropsychidae, of which the mitogenomes of the genus Arctopsyche are published for the first time. Coupled with published hydropsychid mitogenome, we analyzed the nucleotide composition evolutionary rates and gene rearrangements of the mitogenomes among subfamilies. As a result, we found two novel gene rearrangement patterns within Hydropsychidae, including rearrangement of protein-coding genes. Meanwhile, our results consider that the protein-coding gene arrangement of Potamyia can be interpreted by the tandem duplication/random loss (TDRL) model. In addition, the phylogenetic relationships within Hydropsychidae constructed by two strategies (Bayesian inference and maximum likelihood) strongly support the monophyly of Arctopscychinae, Diplectroninae, Hydropsychinae, and Macronematinae. Our study provides new insights into the mechanisms and patterns of mitogenome rearrangements in Hydropsychidae.

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.

Sequencing and phylogenomics of the complete mitochondrial genome of Allodiplogaster sp. (Rhabditida: Diplogasteridae): A new gene order and its phylogenetic implications

Gene order has been utilized as a phylogenetic signal for many taxa. However, its phylogenetic performance has not been evaluated in Nematoda. As there is only one nematode mitogenome available to date, in the Diplogasteridae family, we sequenced the mitogenome of Allodiplogaster sp. and constructed a phylogeny for Nematoda using this updated mitogenome dataset. We then compared this phylogeny to one constructed using gene order. The complete mitochondrial genome of Allodiplogaster sp. was 13,953 bp in size and included 22 tRNAs, two rRNAs, and 12 protein-coding genes. To assess how Allodiplogaster sp. is related to other nematode species, we used Bayesian inference and maximum likelihood algorithms to construct phylogenetic trees of the Nematoda. We found that: 1) The target species Allodiplogaster sp. is closely related to Allodiplogaster sudhausi. The topology of the mitogenome based phylogeny was nearly identical to previous phylogenies created using 18S rRNA data, except for the placement of the Strongyloididae family. 2) The maximum likelihood tree constructed using gene order was roughly consistent with the mitogenome-based tree at the family level, but not at the species level. 3) Protein-coding genes were ordered differently in Allodiplogaster sp. versus Allodiplogaster sudhausi; this represents the first report of such a reordering in the class Chromadorea in our study. Our study confirms that gene order represents useful phylogenetic information for the Nematoda: the maximum likelihood tree based on gene order provided additional support for the nematode phylogeny constructed using molecular data.

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.

Mitogenomic phylogeny of Typhlocybinae (Hemiptera: Cicadellidae) reveals homoplasy in tribal diagnostic morphological traits

The subfamily Typhlocybinae is a ubiquitous, highly diverse group of mostly tiny, delicate leafhoppers. The tribal classification has long been controversial and phylogenetic methods have only recently begun to test the phylogenetic status and relationships of tribes. To shed light on the evolution of Typhlocybinae, we performed phylogenetic analyses based on 28 newly sequenced and 19 previously sequenced mitochondrial genomes representing all currently recognized tribes. The results support the monophyly of the subfamily and its sister-group relationship to Mileewinae. The tribe Zyginellini is polyphyletic with some included genera derived independently within Typhlocybini. Ancestral character state reconstruction suggests that some morphological characters traditionally considered important for diagnosing tribes (presence/absence of ocelli, development of hind wing submarginal vein) are homoplastic. Divergence time estimates indicate that the subfamily arose during the Middle Cretaceous and that the extant tribes arose during the Late Cretaceous. Phylogenetic results support establishment of a new genus, Subtilissimia Yan & Yang gen. nov., with two new species, Subtilissimia fulva Yan & Yang sp. nov. and Subtilissimia pellicula Yan & Yang sp. nov.; but indicate that two previously recognized species of Farynala distinguished only by the direction of curvature of the processes of the aedeagus are synonyms, that is, Farynala dextra Yan & Yang, 2017 equals Farynala sinistra Yan & Yang, 2017 syn. nov. A key to tribes of Typhlocybinae is provided.

Positive Correlation of the Gene Rearrangements and Evolutionary Rates in the Mitochondrial Genomes of Thrips (Insecta: Thysanoptera)

Simple Summary

Aeolothrips, commonly known as banded thrips, is the largest genus of the family Aeolothripidae (predatory thrips). In the current study, we sequenced the mitochondrial genome (mitogenome) of the banded thrip species Aeolothrips xinjiangensis. We found a novel gene arrangement in this mitogenome that has not been reported in Thysanoptera. By comparing the gene order and rearrangement patterns, we found seven identical gene blocks and three identical rearrangement events in two mitogenomes of banded thrips. There was marked variation in the mitochondrial gene order across thrip species, with only two conserved gene blocks shared by all 14 thrips. In addition, we found a positive correlation between the degree of gene rearrangement and evolutionary rate. Our results suggested that the mitogenomes of thrips have tended to be stable since their massive rearrangement.

Abstract

Extensive gene rearrangement is characteristic in the mitogenomes of thrips (Thysanoptera), but the historical process giving rise to the contemporary gene rearrangement pattern remains unclear. To better understand the evolutionary processes of gene rearrangement in the mitogenomes of thrips, we sequenced the mitogenome of the banded thrip species Aeolothrips xinjiangensis. First, we found a novel mitochondrial gene order in this species. This mitogenome is 16,947 bp in length and encodes the typical 37 coding genes (13 protein-coding genes, 22 tRNA genes, and two rRNA genes) of insects. The gene arrangement was dramatically different from the putative ancestral mitogenome, with 26 genes being translocated, eight of which were inverted. Moreover, we found a novel, conserved gene block, trnC-trnY, which has not been previously reported in the mitogenomes of thrips. With this newly assembled mitogenome, we compared mitogenome sequences across Thysanoptera to assess the evolutionary processes giving rise to the current gene rearrangement pattern in thrips. Seven identical gene blocks were shared by two sequenced banded thrip mitogenomes, while the reversal of ND2 combined with TDRL events resulted in the different gene orders of these two species. In phylogenetic analysis, the monophyly of the suborders and families of Thysanoptera was well supported. Across the gene orders of 14 thrips, only two conserved gene blocks, ATP8-ATP6 and ND4-ND4L, could be found. Correlation analysis showed that the degree of gene rearrangement was positively correlated with the non-synonymous substitution rate in thrips. Our study suggests that the mitogenomes of thrips remain stable over long evolutionary timescales after massive rearrangement during early diversification.

Novel mitochondrial gene rearrangements pattern in the millipede Polydesmus sp. GZCS-2019 and phylogenetic analysis of the Myriapoda

The subphylum Myriapoda included four extant classes (Chilopoda, Symphyla, Diplopoda, and Pauropoda). Due to the limitation of taxon sampling, the phylogenetic relationships within Myriapoda remained contentious, especially for Diplopoda. Herein, we determined the complete mitochondrial genome of Polydesmus sp. GZCS-2019 (Myriapoda: Polydesmida) and the mitochondrial genomes are circular molecules of 15,036 bp, with all genes encoded on + strand. The A+T content is 66.1%, making the chain asymmetric, and exhibits negative AT-skew (-0.236). Several genes rearrangements were detected and we propose a new rearrangement model: "TD (N\R) L + C" based on the genome-scale duplication + (non-random/random) loss + recombination. Phylogenetic analyses demonstrated that Chilopoda and Symphyla both were monophyletic group, whereas Pauropoda was embedded in Diplopoda to form the Dignatha. Divergence time showed the first split of Myriapoda occurred between the Chilopoda and other classes (Wenlock period of Silurian). We combine phylogenetic analysis, divergence time, and gene arrangement to yield valuable insights into the evolutionary history and classification relationship of Myriapoda and these results support a monophyletic Progoneata and the relationship (Chilopoda + (Symphyla + (Diplopoda + Pauropoda))) within myriapod. Our results help to better explain the gene rearrangement events of the invertebrate mitogenome and lay the foundation for further phylogenetic study of Myriapoda.

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.

Characterization of the complete mitochondrial genome of Myrmuslateralis (Heteroptera, Rhopalidae) and its implication for phylogenetic analyses

Mitochondrial genomes (mitogenomes) are widely used in research studies on phylogenetic relationships and evolutionary history. Here, we sequenced and analyzed the mitogenome of the scentless plant bug Myrmuslateralis Hsiao, 1964 (Heteroptera, Rhopalidae). The complete 17,309 bp genome encoded 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and a control region. The mitogenome revealed a high A+T content (75.8%), a positive AT-skew (0.092), and a negative GC-skew (-0.165). All 13 PCGs were found to start with ATN codons, except for cox1, in which TTG was the start codon. The Ka/Ks ratios of 13 PCGs were all lower than 1, indicating that purifying selection evolved in these genes. All tRNAs could be folded into the typical cloverleaf secondary structure, except for trnS1 and trnV, which lack dihydrouridine arms. Phylogenetic trees were constructed and analyzed based on the PCG+rRNA from 38 mitogenomes, using maximum likelihood and Bayesian inference methods, showed that M.lateralis and Chorosomamacilentum Stål, 1858 grouped together in the tribe Chorosomatini. In addition, Coreoidea and Pyrrhocoroidea were sister groups among the superfamilies of Trichophora, and Rhopalidae was a sister group to Alydidae + Coreidae.

The complete mitochondrial genome of Aeolothrips indicus Bhatti, 1964 (Thysanoptera: Thripidae)

Here, we have generated the complete mitochondrial sequence of Aeolothrips indicus Bhatti, 1964. So far, this is the first largest mitogenome with 17,042 bp length in order Thysanoptera. It includes 13 protein-coding genes, 22 transfer RNA genes, and two ribosomal RNA genes along with three non-coding regions. AT composition of A. indicus is 72.5% (37.7% A and 34.8% T) and GC 27.5% (15.6% C and 11.9% G). The constructed phylogeny revealed the monophyly of family Aeolothripidae in the order Thysanoptera. The data would provide further insight into the evolution and phylogeny of the order Thysanoptera.

Mitochondrial Genomes from Two Specialized Subfamilies of Reduviidae (Insecta: Hemiptera) Reveal Novel Gene Rearrangements of True Bugs

Reduviidae, a hyper-diverse family, comprise 25 subfamilies with nearly 7000 species and include many natural enemies of crop pests and vectors of human disease. To date, 75 mitochondrial genomes (mitogenomes) of assassin bugs from only 11 subfamilies have been reported. The limited sampling of mitogenome at higher categories hinders a deep understanding of mitogenome evolution and reduviid phylogeny. In this study, the first mitogenomes of Holoptilinae (Ptilocnemus lemur) and Emesinae (Ischnobaenella hainana) were sequenced. Two novel gene orders were detected in the newly sequenced mitogenomes. Combined 421 heteropteran mitogenomes, we identified 21 different gene orders and six gene rearrangement units located in three gene blocks. Comparative analyses of the diversity of gene order for each unit reveal that the tRNA gene cluster trnI-trnQ-trnM is the hotspot of heteropteran gene rearrangement. Furthermore, combined analyses of the gene rearrangement richness of each unit and the whole mitogenome among heteropteran lineages confirm Reduviidae as a 'hot-spot group' of gene rearrangement in Heteroptera. The phylogenetic analyses corroborate the current view of phylogenetic relationships between basal groups of Reduviidae with high support values. Our study provides deeper insights into the evolution of mitochondrial gene arrangement in Heteroptera and the early divergence of reduviids.

Insight into the Phylogenetic Relationships among Three Subfamilies within Heptageniidae (Insecta: Ephemeroptera) along with Low-Temperature Selection Pressure Analyses Using Mitogenomes

We determined 15 complete and two nearly complete mitogenomes of Heptageniidae belonging to three subfamilies (Heptageniinae, Rhithrogeninae, and Ecdyonurinae) and six genera (Afronurus, Epeorus, Leucrocuta, Maccaffertium, Stenacron, and Stenonema). Species of Rhithrogeninae and Ecdyonurinae had the same gene rearrangement of CR-I-M-Q-M-ND2, whereas a novel gene rearrangement of CR-I-M-Q-NCR-ND2 was found in Heptageniinae. Non-coding regions (NCRs) of 25-47 bp located between trnA and trnR were observed in all mayflies of Heptageniidae, which may be a synapomorphy for Heptageniidae. Both the BI and ML phylogenetic analyses supported the monophyly of Heptageniidae and its subfamilies (Heptageniinae, Rhithrogeninae, and Ecdyonurinae). The phylogenetic results combined with gene rearrangements and NCR locations confirmed the relationship of the subfamilies as (Heptageniinae + (Rhithrogeninae + Ecdyonurinae)). To assess the effects of low-temperature stress on Heptageniidae species from Ottawa, Canada, we found 27 positive selection sites in eight protein-coding genes (PCGs) using the branch-site model. The selection pressure analyses suggested that mitochondrial PCGs underwent positive selection to meet the energy requirements under low-temperature stress.

Mitogenomes Reveal Alternative Initiation Codons and Lineage-Specific Gene Order Conservation in Echinoderms

The mitochondrial genetic code is much more varied than the standard genetic code. The invertebrate mitochondrial code, for instance, comprises six initiation codons, including five alternative start codons. However, only two initiation codons are known in the echinoderm and flatworm mitochondrial code, the canonical ATG and alternative GTG. Here, we analyzed 23 Asteroidea mitogenomes, including ten newly sequenced species and unambiguously identified at least two other start codons, ATT and ATC, both of which also initiate translation of mitochondrial genes in other invertebrates. These findings underscore the diversity of the genetic code and expand upon the suite of initiation codons among echinoderms to avoid erroneous annotations. Our analyses have also uncovered the remarkable conservation of gene order among asteroids, echinoids, and holothuroids, with only an interchange between two gene positions in asteroids over ∼500 Ma of echinoderm evolution.

Complete mitochondrial genome sequence for the Thrips hawaiiensis (Thysanoptera: Thripidae)

Thrips hawaiiensis (Morgan) (Thysanoptera: Thripidae) is a common Thysanoptera insect widely distributed in Asia and the Pacific, it damages various plants. In this study the complete mitochondrial genome of T. hawaiiensis was sequenced and characterized by using next-generation sequencing technique. The total length of the complete genome is 15,357 bp and A + T content of 77.8% (GeneBank accession No. MW582621). The T. hawaiiensis mitochondrial genome consists of 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, 22 transfer RNA genes (tRNAs) and 2 non-coding control regions (Dloop region). According to previous studies, only a few complete mitochondrial genomes from Order Thysanoptera have been reported. Thus, T. hawaiiensis complete mitochondrial genome sequence reported will provide molecular information for mitochondrial genome research on Thysanoptera.

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.

The queen conch mitogenome: intra- and interspecific mitogenomic variability in Strombidae and phylogenetic considerations within the Hypsogastropoda

Aliger gigas is an economically important and vulnerable marine species. We present a new mitogenome of A. gigas from the Mexican Caribbean and use the eight publicly available Strombidae mitogenomes to analyze intra- and interspecific variation. We present the most complete phylogenomic understanding of Hypsogastropoda to date (17 superfamilies, 39 families, 85 genera, 109 species) to revisit the phylogenetic position of the Stromboidea and evaluate divergence times throughout the phylogeny. The A. gigas mitogenome comprises 15,460 bp including 13 PCGs, 22 tRNAs, and two rRNAs. Nucleotide diversity suggested divergence between the Mexican and Colombian lineages of A. gigas. Interspecific divergence showed high differentiation among Strombidae species and demonstrated a close relationship between A. gigas and Strombus pugilis, between Lambis lambis and Harpago chiragra, and among Tridentarius dentatus/Laevistrombus canarium/Ministrombus variabilis. At the intraspecific level, the gene showing the highest differentiation is ATP8 and the lowest is NAD4L, whereas at the interspecific level the NAD genes show the highest variation and the COX genes the lowest. Phylogenomic analyses confirm that Stromboidea belongs in the non-Latrogastropoda clade and includes Xenophoridea. The phylogenomic position of other superfamilies, including those of previously uncertain affiliation, is also discussed. Finally, our data indicated that Stromboidea diverged into two principal clades in the early Cretaceous while Strombidae diversified in the Paleocene, and lineage diversification within A. gigas took place in the Pleistocene.

Highly rearranged mitochondrial genome in Falcolipeurus lice (Phthiraptera: Philopteridae) from endangered eagles

Background

Fragmented mitochondrial (mt) genomes and extensive mt gene rearrangements have been frequently reported from parasitic lice (Insecta: Phthiraptera). However, relatively little is known about the mt genomes from the family Philopteridae, the most species-rich family within the suborder Ischnocera.

Methods

Herein, we use next-generation sequencing to decode the mt genome of Falcolipeurus suturalis and compare it with the mt genome of F. quadripustulatus. Phylogenetic relationships within the family Philopteridae were inferred from the concatenated 13 protein-coding genes of the two Falcolipeurus lice and members of the family Philopteridae using Bayesian inference (BI) and maximum likelihood (ML) methods.

Results

The complete mt genome of F. suturalis is a circular, double-stranded DNA molecule 16,659bp in size that contains 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and three non-coding regions. The gene order of the F. suturalis mt genome is rearranged relative to that of F. quadripustulatus, and is radically different from both other louse species and the putative ancestral insect. Phylogenetic analyses revealed clear genetic distinctiveness between F. suturalis and F. quadripustulatus (Bayesian posterior probabilities=1.0 and bootstrapping frequencies=100), and that the genus Falcolipeurus is sister to the genus Ibidoecus (Bayesian posterior probabilities=1.0 and bootstrapping frequencies=100).

Conclusions

These datasets help to better understand gene rearrangements in lice and the phylogenetic position of Falcolipeurus and provide useful genetic markers for systematic studies of bird lice.

Graphic abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04776-5.

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.

Mitochondrial analysis of oribatid mites provides insights into their atypical tRNA annotation, genome rearrangement and evolution

Background

The mitochondrial (mt) genomes of Sarcoptiformes mites typically contain 37 genes. Although the loss of genes is rare in Sarcoptiformes mite mitogenomes, two of the six previously reported oribatid mites (Acariforms: Sarcoptiformes) are reported to have lost parts of their tRNA genes. To confirm whether the tRNA genes were indeed lost and whether the loss is universal, we re-annotated the available oribatid mite sequences and sequenced the mitogenome of Oribatula sakamorii.

Methods

The mitogenome of O. sakamorii was sequenced using an Illumina HiSeq sequencer. The mt tRNA gene was annotated using multi-software combined with a manual annotation approach. Phylogenetic analyses were performed using the maximum likelihood and Bayesian inference methods with concatenated nucleotide and amino acid sequences.

Results

The mitogenomes of O. sakamorii contained 37 genes, including 22 tRNA genes. We identified all mt tRNA genes that were reported as “lost” in Steganacarus magnus and Paraleius leontonychus and revealed certain atypical tRNA annotation errors in oribatid mite sequences. Oribatid mite mitogenomes are characterized by low rates of genetic rearrangement, with six or seven gene blocks conserved between the mitogenome of all species and that of ancestral arthropods. Considering the relative order of the major genes (protein-coding genes and rRNAs), only one or two genes were rearranged with respect to their positions in the ancestral genome. We explored the phylogenetic relationships among the available oribatid mites, and the results confirmed the systematic position of Hermannia in the Crotonioidea superfamily. This was also supported by the synapomorphic gene-derived boundaries.

Conclusions

The tRNA “lost” phenomenon is not universal in oribatid mites. Rather, highly atypical secondary structure of the inferred mt tRNA genes made them unidentifiable using a single type of tRNA search program. The use of multi-software combined with a manual annotation approach can improve the accuracy of tRNA gene annotation. In addition, we identified the precise systematic position of Hermannia and validated that Astigmata is nested in Oribatida.

Graphic Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04719-0.