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 first report of complete mitogenomes of two endangered species of genus Propomacrus (Coleoptera: Scarabaeidae: Euchirinae) and phylogenetic implications

To understand the mitochondrial genome structure of two endangered and long-armed scarab beetles, Propomacrus davidi and Propomacrus bimucronatus, their complete mitogenomes were sequenced for the first time in this study. The complete mitogenomes of P. davidi and P. bimucronatus were 18, 042 bp and 18, 104 bp in length, respectively. The gene orders of their mitogenomes were highly consistent with other Coleopteran species, and the typical ATN was used as the start codon in most protein coding genes. The incomplete stop codon T was used in cox1, cox2, and nad5, and TAN was used as a complete stop codon in most protein coding genes. All predicted tRNAs could form a typical cloverleaf secondary structure, except that trnS1 lacked the dihydrouridine arm. Based on the maximum likelihood and the Bayesian inference methods, phylogenetic trees of 50 species were reconstructed. The results showed that P. davidi, P. bimucronatus, Cheirotonus jansoni and Cheirotonus gestroi clustered in the same branch, and were the most closely related. The results supported that subfamily Euchirinae is a monophyletic group of Scarabaeidae, which was consistent with the morphological classification. These molecular data enriched the complete mitogenome database of Euchirinae, and improved our understanding of the phylogenetic relationship and evolutionary characteristics of these two endangered species.

Elucidating Scarab Divergence in an Evolutionary-Ecological Context through the Comprehensive Analysis of the Complete Mitogenome of Anomala

Anomala Samouelle, 1819 is one of the specious genera of Coleoptera, with over 1000 known species, and includes some of the most destructive pests of crops or forests. Morphological convergence is a common phenomenon within this genus, making the identification of closely related species very difficult. To explore the phylogenetic placement of Anomalini and provide a basis for the classification and identification of Anomala, we comparatively analyzed the complete mitogenome of three Anomala species (A. exoleta, A. perplexa diana, and A. praecoxalis). Based on all accessible mitogenome data, we performed comparative mitochondrial genomics analysis of this genus and reconstructed the phylogenetic relationships of Scarabaeidae based on two datasets (protein-coding genes and amino acids) and two methods (Bayesian approach and maximum likelihood). The phylogenetic relationships found in this study highly support that the groups of Aphodiinae, Cetoniinae, Dynastinae, Rutelinae and Scarabaeinae are monophyletic. Interestingly, the phylogenetic clustering relationship was highly consistent with the Scarabaeidae diet, indicating that the herbivorous species and dung-eating species are clustered separately. The phylogenetic tree showed that the subfamily Melolonthinae and the genus Anomala are not monophyletic, suggesting that these two groups should be further revised with more data.

The mitochondrial genome of an important edible insect species, the African palm weevil (Rhynchophorus phoenicis)

The African palm weevil (Rhynchophorus phoenicis) is a species of high economic importance in sub-Saharan Africa, both as a culturally traditional edible insect and as an agricultural pest. Here we provide a de novo assembly and annotation for the mitochondrial genome of this species from whole-genome sequence data. The mitogenome was AT-rich and 17,161bp in length, containing 13 protein-coding, 22 transfer RNA, and two ribosomal RNA genes. Phylogenetic reconstruction showed the African palm weevil to cluster within the genus Rhynchophorus and the weevil tribe Rhynchophorini. This mitogenome will be important for future genetic research into this emerging edible insect species.

The Phylogenetic Relationship of Lamiinae (Coleoptera: Cerambycidae) Using Mitochondrial Genomes

Lamiinae is the largest subfamily of the Cerambycidae (longhorn beetles), with approximately 21,863 described species. Previous phylogenetic studies of Lamiinae showed that this subfamily was monophyletic, but the relationship between the tribes of Lamiinae is still controversial. Partial molecular data and species morphological characteristics are not sufficient to resolve species phylogenetic studies perfectly. At the same time, the full mitochondrial genome contains more comprehensive genetic data. Benefiting from the development of next-generation sequencing (NGS), mitochondrial genomes can be easily acquired and used as reliable molecular markers to investigate phylogenetic relationships within Cerambycidae. Using NGS technology, we obtained 11 mitochondrial genome sequences of Lamiinae species. Based on this newly generated mitochondrial genome dataset matrix, we reconstructed the phylogeny of Lamiinae. The Bayesian Inference and Maximum Likelihood analyses strongly support the monophyly of four tribes (Lamiini, Batocerini, Mesosini, and Saperdini), whereas the tribe Acanthocinini was identified as paraphyletic. Other mitochondrial structural features were also observed: the start codon in the nad1 gene of all 11 mitochondrial genomes is TTG; 17-22 bp intergenic spacers (IGS) with a 'TACTA' motif were found between trnS2 and nad1. Moreover, two long IGS were found in Mesosa myops and Batocera sp. Tandem repeats were found in the IGS of Batocera sp.

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.

Novel Gene Rearrangements in Mitochondrial Genomes of four families of Praying Mantises (Insecta, Mantodea) and Phylogenetic Relationships of Mantodea

The mitochondrial genome (mitogenome) plays an important role in phylogenetic studies of many species. The mitogenomes of many praying mantis groups have been well-studied, but mitogenomes of special mimic praying mantises, especially Acanthopoidea and Galinthiadoidea species, are still sorely lacking in the NCBI database. The present study analyzes five mitogenomes from four species of Acanthopoidea (Angela sp., Callibia diana, Coptopteryx sp., Raptrix fusca) and one of Galinthiadoidea (Galinthias amoena) that were sequenced by the primer-walking method. Among Angela sp. and Coptopteryx sp., three gene rearrangements were detected in ND3-A-R-N-S-E-F and COX1-L2-COX2 gene regions, two of which were novel. In addition, individual tandem repeats were found in control regions of four mitogenomes (Angela sp., C. diana, Coptopteryx sp., G. amoena). For those, plausible explanations were derived from the tandem duplication-random loss (TDRL) model and the slipped-strand mispairing model. One potential motif was found in Acanthopidae that was seen as a synapomorphy. Several conserved block sequences (CBSs) were detected within Acanthopoidea that paved the way for the design of specific primers. Via BI and ML analysis, based on four datasets (PCG12, PCG12R, PCG123, PCG123R), the merged phylogenetic tree within Mantodea was reconstructed. This showed that the monophyly of Acanthopoidea was supported and that the PCG12R dataset was the most suitable for reconstructing the phylogenetic tree within Mantodea.

Mitochondrial phylogenomics reveals deep relationships of scarab beetles (Coleoptera, Scarabaeidae)

In this study, we newly sequenced the complete mitochondrial genomes (mitogenomes) of two phytophagous scarab beetles, and investigated the deep level relationships within Scarabaeidae combined with other published beetle mitogenome sequences. The complete mitogenomes of Dicronocephalus adamsi Pascoe (Cetoniinae) and Amphimallon sp. (Melolonthinae) are 15,563 bp and 17,433 bp in size, respectively. Both mitogenomes have the typical set of 37 genes (13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes) and an A+T-rich region, with the same gene arrangement found in the majority of beetles. The secondary structures for ribosomal RNA genes (rrnL and rrnS) were inferred by comparative analysis method. Results from phylogenetic analyses provide support for major lineages and current classification of Scarabaeidae. Amino acid data recovered Scarabaeidae as monophyletic. The Scarabaeidae was split into two clades. One clade contained the subfamilies Scarabaeinae and Aphodiinae. The other major clade contained the subfamilies Dynastinae, Rutelinae, Cetoniinae, Melolonthinae and Sericini. The monophyly of Scarabaeinae, Aphodiinae, Dynastinae, Cetoniinae and Sericini were strongly supported. The Scarabaeinae was the sister group of Aphodiinae. The Cetoniinae was sister to the Dynastinae + Rutelinae clade. The Melolonthinae was a non-monophyletic group. The removal of fast-evolving sites from nucleotide dataset using a pattern sorting method (OV-sorting) supported the family Scarabaeidae as a monophyletic group. At the tribe level, the Onthophagini was non-monophyletic with respect to Oniticellini. Ateuchini was sister to a large clade comprising the tribes Onthophagini, Oniticellini and Onitini. Eurysternini was a sister group of the Phanaeini + Ateuchini clade.

Complete Mitogenomes of Ticks Ixodes acutitarsus and Ixodes ovatus Parasitizing Giant Panda: Deep Insights into the Comparative Mitogenomic and Phylogenetic Relationship of Ixodidae Species

Ticks rank second in the world as vectors of disease. Tick infestation is one of the factors threatening the health and survival of giant pandas. Here, we describe the mitogenomes of Ixodes acutitarsus and Ixodes ovatus parasitizing giant pandas, and perform comparative and phylogenetic genomic analyses on the newly sequenced and other available mitogenomes of hard ticks. All six newly determined mitogenomes contain a typical gene component and share an ancient Arthropoda gene arrangement pattern. Our study suggests that I. ovatus is a species complex with high genetic divergence, indicating that different clades of I. ovatus represent distinct species. Comparative mitogenomic analyses show that the average A + T content of Ixodidae mitogenomes is 78.08%, their GC-skews are strongly negative, while AT-skews fluctuate around 0. A large number of microsatellites are detected in Ixodidae mitogenomes, and the main microsatellite motifs are mononucleotide A and trinucleotide AAT. We summarize five gene arrangement types, and identify the trnY-COX1-trnS1-COX2-trnK-ATP8-ATP6-COX3-trnG fragment is the most conserved region, whereas the region near the control region is the rearrangement hotspot in Ixodidae mitogenomes. The phylogenetic trees based on 15 genes provide a very convincing relationship (Ixodes + (Robertsicus + ((Bothriocroton + Haemaphysalis) + (Amblyomma + (Dermacentor + (Rhipicentor + (Hyalomma + Rhipicephalus))))))) with very strong supports. Remarkably, Archaeocroton sphenodonti is embedded in the Haemaphysalis clade with strong supports, resulting in paraphyly of the Haemaphysalis genus, so in-depth morphological and molecular studies are essential to determine the taxonomic status of A. sphenodonti and its closely related species. Our results provide new insights into the molecular phylogeny and evolution of hard ticks, as well as basic data for population genetics assessment and efficient surveillance and control for the giant panda-infesting ticks.

The Genetic Diversity and the Divergence Time in Extant Primitive Mayfly, Siphluriscus chinensis Ulmer, 1920 Using the Mitochondrial Genome

In this study, the mitochondrial (mt) genomes of Siphluriscus chinensis (Ephemeroptera: Siphluriscidae) were evaluated in specimens collected from two sites in China: Niutou Mountain, Zhejiang Province (S. chinensis NTS) and Leigong Mountain, Guizhou Province (S. chinensis LGS) and were successfully sequenced. The lengths of the mt genomes of S. chinensis NTS and S. chinensis LGS were 15,904 bp (ON729390) and 15,212 bp (ON729391), respectively. However, an in-depth comparison of the two mt genomes showed significant differences between the specimens collected from the two sites. A detailed analysis of the genetic distance between S. chinensis NTS and S. chinensis LGS was undertaken to further achieve an accurate delimitation of S. chinensis. The genetic distance between S. chinensis NTS and the other three species within Siphluriscidae was a high value, above 12.2%. The two mt genomes were used to reconstruct phylogenetic relationships and estimate divergence time. The results demonstrated robust differences between S. chinensis NTS and S. chinensis LGS, which revealed that a kind of cryptic species existed. Maximum likelihood (ML) and Bayesian inference (BI) analyses produced well-supported phylogenetic trees that showed evolutionary relationships between Siphluriscidae (((S. chinensis HQ875717 + S. chinensis MF352165) + S. chinensis LGS) + S. chinensis NTS). The most recent common ancestor (MRCA) of four species within Siphluriscidae began to diversify during the Neogene [11.80 million years ago (Mya); 95% highest posterior densities (HPD) = 6.17–19.28 Mya], and S. chinensis NTS was first to diverge from the branches of S. chinensis LGS. In short, based on mitochondrial genomes, our results showed that the specimens collected from Leigong Mountain, Guizhou Province (S. chinensis LGS) belonged to S. chinensis, and the specimens collected from Niutou Mountain, Zhejiang Province (S. chinensis NTS) were a cryptic species of S. chinensis.

Revisiting trends in morphology of antennal sensilla in scarabaeoid beetles

Phytophagous scarab beetles associated with angiosperms have characteristically enlarged lamellate antennae and exhibit a striking morphological variation of sensilla. In this study, we compared the morphology of antennal surface of 62 species Scarabaeoidea using SEM microscopy, particularly also in light of their evolution in association with angiosperms. We investigated the correlation of antennal sensilla morphology, i.e., their structure and distribution, with species diversity and lineage diversification rates. A high diversity of sensilla was observed but also multiple transitional forms, even on the same antennomere. We interpreted this as evidence for a high evolutionary plasticity. We recognized clear patterns of convergence and repeated evolution of certain types of placoid sensilla. One main tendency found in the phytophagous Pleurostict chafers was a shift from sensilla trichodea to placoid-like sensilla, apparently also enhanced by the increase of the lamellate antennal surface, either by size or number of the lamellae. This trend occurred not only in the Pleurosticts, but also in Glaphyridae, a second angiosperm-associated lineage of Scarabaeoidea. However, our results suggest no direct relation between species diversity or the rate of diversification and general sensilla morphology, i.e., the origin of placoid sensilla. This could be explained not only by species-poor lineages also possessing placoid sensilla but also by otherwise successful and species rich groups having sensilla trichodea (e.g., dung beetles). Results further reveal the need to refine current phylogenetic hypotheses by more comprehensive taxon sampling and to expand the molecular characterization of pheromones and odor binding proteins to better understand the role of chemical communication in scarab diversification.

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.

The Genetic Diversity in Thereuonema tuberculata (Wood, 1862) (Scutigeromorpha: Scutigeridae) and the Phylogenetic Relationship of Scutigeromorpha Using the Mitochondrial Genome

Based on morphological characteristics to make species identification, the cryptic species of the Scutigeromorpha can be greatly underestimated. The mitochondrial genome provides a desirable tool for the biological identifications and the discovery of the cryptic species. The capacity to acquire mitochondrial genome sequences has substantially improved in recent years using next-generation sequencing (NGS) technology. On the basis of the next-generation sequencing, we obtained four complete mitochondrial genomes of Thereuonema tuberculata (Wood, 1862) from Nanyang, Henan Province (NY), Nanchang, Jiangxi Province (NC), Jinan, Shandong Province (JN), and Dali, Yunnan Province (DL) in China with GenBank numbers OK513221, OL449685, ON058988 and ON058989, respectively. The lengths of the four mitochondrial genomes ranged from 14,903 to 14,909 bp. The composition and order of genes of the four mitochondrial genomes were identical to the published mitochondrial genome of Scutigera coleoptrata (Linnaeus, 1758) (Scutigeromorpha: Scutigerdae). It was the first time that the tandem repeats in the control region were detected in Scutigeromorpha. We also calculated the corrected pairwise genetic distance of four complete mitochondrial genomes of T. tuberculata, ranging from 7.7 to 15.2%. The results showed that the T.tuberculata NC belonged to the typical sample of T. tuberculata, and T. tuberculata DL was hypothesized as a cryptic species of T. tuberculata. Meanwhile, T. tuberculata NY and T. tuberculata JN were hypothesized as potential cryptic species of T. tuberculata in this study. In both BI and ML trees, the monophyly of Scutigeromorpha, Scolopendromorpha, Geophilomorpha, and Lithobiomorpha was forcefully advocated. Moreover, Scutigeromorpha was recovered as the sister clade of (Scolopendromorpha + (Lithobiomorpha + Geophilomorpha)). Four specimens of T. tuberculata were clustered into one clade, which was the sister to the clade of S. coleoptrata.

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

Simple Summary

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

Abstract

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

The complete mitochondrial genome of Leptomantella tonkinae (Hebard, 1920) (Mantodea: Leptomantellidae) and its phylogeny

The complete mitochondrial (mt) genome of Leptomantella tonkinae (Hebard, 1920) was 15,527 bp in length and contained 13 protein-coding genes, 22 transfer RNAs, two ribosomal RNAs, and one control region. The gene arrangement of mt genome of L. tonkinae was identical to the primitive mantis. The overall AT content of the mt genome was 74%. In ML and BI phylogenetic analyses, the monophyly of Leptomantellidae was robustly supported and the clade of Leptomantellidae is a sister clade to the group of (Gonypetidae+(Leptomantellidae+(Amorphoscelidae+Nanomantidae))).