Next Article in Journal
Emerald Ash Borer Infestation-Induced Elevated Negative Correlations and Core Genera Shift in the Endophyte Community of Fraxinus bungeana
Previous Article in Journal
Ultrastructure and Spectral Characteristics of the Compound Eye of Asiophrida xanthospilota (Baly, 1881) (Coleoptera, Chrysomelidae)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Insects
Volume 15
Issue 7
10.3390/insects15070533
insects-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Liangpeng Ji
,
Zhancheng Jia
and
Xiaoshuan Bai
*
College of Life Sciences and Technology, Inner Mongolia Normal University, Hohhot 010022, China
*
Author to whom correspondence should be addressed.
Insects 2024, 15(7), 533; https://doi.org/10.3390/insects15070533
Submission received: 21 June 2024 / Revised: 9 July 2024 / Accepted: 12 July 2024 / Published: 14 July 2024
(This article belongs to the Section Insect Systematics, Phylogeny and Evolution)
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

Aradidae is a large family in Hemiptera, which feeds on mycelium, and its phylogenetic status is not fully understood. Mitochondrial genome sequences can be used to study species identification and phylogeny, and provide valuable molecular markers for further genetic research. In this paper, we sequenced three species of the genus Yangiella for the first time, assembled three mitochondrial genomes, and compared the general characteristics of the mitochondrial genomes of the three species. It was found that the base composition and mitochondrial genome structure of the three species were highly similar. The phylogeny of Yangiella was also discussed. Based on the phylogenetic analyses results of two data matrices, the status of Yangiella was discussed and the monophyly of this genus was verified. In addition, we infer that the genus Yangiella diverged about 57 million years ago.

Abstract

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

1. Introduction

Aradidae is a larger family in Hemiptera, with 8 subfamilies and about 2000 species [1]. Mezirinae is the largest subfamily and is found in all geographical regions of the world [2]. Yangiella is a genus of Mezirinae [3]. This genus of insects usually lives under the bark of fallen trees, and is sometimes found in gaps on the bark surface [4]. They live in groups and feed on hyphae, mainly distributed in Southwestern China and parts of Southeast Asia [5]. Due to its high concealment and lesser contact with humans, past research on this genus has mainly focused on traditional morphological classification, and molecular biology research has not been carried out, resulting in its classification status still being unclear.
Mitochondrial genomes have the characteristics of maternal inheritance, fast evolution, easy separation and purification, and are widely used in the study of phylogenetic relationship of insects [6,7] and molecular evolution [8,9]. The typical insect mitochondrial genome is a double-stranded, circular, closed DNA molecular structure with a length of 15–18 kb, containing 37 genes (composed of 13 protein-coding genes and 24 RNA genes) [10], and a non-coding control region containing the origin of replication [11]. In the past, most of the studies on the molecular system of Aradidae were based on partial gene fragments, and Heisset et al. analyzed the evolution of three subfamilies based on cox1 fragments [12,13]. Marchal et al. used cox1, 16S rRNA, and 28S rRNA gene fragments to analyze the phylogeny of 79 species in eight subfamilies. They believed that Prosympiestinae and Chinamyersiinae were paraphyletic, and the remaining six subfamilies were monophyletic [14]. However, compared with the complete mitochondrial genome, a single gene can only reflect part of the information and provide limited information. The complete mitochondrial genome sequence contains more abundant information [15], which is a powerful evidence form for phylogenetic research [16].
In this study, the mitochondrial genomes of three Yangiella species were sequenced and annotated for the first time. The basic characteristics of the mitochondrial genome, codon usage bias, and tRNA gene secondary structure were described and analyzed. The internal relationship of Aradidae was elucidated by analyzing the evolutionary rate and genome rearrangement of each species of Aradidae. In addition, based on 13 protein-coding genes and 2 RNA genes, the phylogenetic trees of 16 species (including 14 Aradidae insects and 2 outgroups) were reconstructed by Bayesian inference (BI) and maximum likelihood (ML) methods, and the divergence time of Aradidae was estimated based on phylogenetic trees and fossil data.

2. Materials and Methods

2.1. Specimen Collection, Extraction, and Sequencing

Samples of Yangiella mimetica (Hsiao, 1964 [3]), Yangiella montana (Zhang, Bai, Heiss, and Cai, 2010 [5]) and Yangiella sp. were collected in China (Table 1). The specimens were immersed in 95% ethanol during collection, and then transferred of −20 °C for long-term preservation in the Inner Mongolia Normal Entomological Museum (Hohhot, Inner Mongolia). The total DNA was extracted from the leg and chest muscle tissues of insects. Using the kit of Tiangen Biochemical Technology Co., Ltd. (Beijing, China), DNA was extracted according to the instructions in the instructions. The processed samples were sent to Berry Genomics Co., Ltd. (Beijing, China) for sequencing (DNA concentration ≥ 20 ng/μL, total ≥ 500 ng). Second-generation sequencing technology was used, and the sequencing mode was Novaseq 6000-S4-150 PE.

2.2. Gene Assembly, Annotation, and Analysis

Raw read data with a double-end sequencing length of 150 bp were obtained, and Fastp v0.23.0 was used to perform quality control on the raw data [17], and the reads below the average quality value of the reads were removed. The reads with N content less than 5 were screened to obtain clean data with a size of 5G. SPAdes v3.15.5 (https://github.com/ablab/spades (accessed on 15 March 2024)) was used for assembly [18]. Based on the assembly process, the contig sequence was obtained by iterative method to construct the local blast database. Using the cox1 sequence of Neuroroctenus yunnanensis (Hsiao, 1964), which belongs to Mezirinae, as bait, the BLAST function was used to screen contigs of the target mitochondrial genome.
The mitochondrial genome sequence was imported into the online tool MITOS (http://mitos.bioinf.uni-leipzig.de/index.py (accessed on 22 March 2024)), and set to invertebrate mitochondrial coding mode for automated annotation [19]. The annotation results were imported into Geneous Prime v11 [20], and the gene boundaries were manually corrected by comparing them with related species. The location of the tRNA gene was reconfirmed using the tRNAscan-SE online tool (Cover cutoff = 5) [21], and Adobe Illustrator CS 2024 software was used for drawing. The CGView Server online server (https://cgview.ca/ (accessed on 10 April 2024)) was used to visualize the mitochondrial map [22]. The gene rearrangement of Aradidae was studied by comparing with the mitochondrial gene arrangement of Macroscytus subaeneus (Dallas, 1851) from the Cydnidae. Finally, the newly sequenced mitochondrial genomes were submitted to the GenBank database.
The nucleotide composition of the gene was calculated using MEGA11 [23]. AT-skew and GC-skew are calculated by AT-skew = (A − T)/(A + T) and GC-skew = (G − C)/(G + C) [24], respectively. DnaSP v6.12.03 [25] was used to calculate the non-synonymous substitution rate (Ka), synonymous substitution rate (Ks), and nucleotide diversity (Pi) of protein-coding genes. Relative synonymous codon usage (RSCU) and effective number of codon (ENC) were analyzed by CodonWv 1.4.4 software.
Neutrality plot analysis was performed on Aradidae species. GC12 (the mean value of G + C content at codons 1 and 2 of protein-coding genes) was used as the ordinate, and GC3 (the G + C content at codon 3 of protein-coding genes) was used as the abscissa. The difference in base content at different sites of protein-coding genes was judged by a slope to determine whether it demonstrated neutral evolution [26]. ENC-plot analysis was performed, and the scatter plot was drawn using Microsoft Excel 2019 with GC3 as the abscissa and ENC as the ordinate [27]. The expected value was calculated according to the formula ENC = 2 + GC3 + (29/(GC32 + (1 − GC3)2)) [28], and the standard curve was drawn. Finally, DAMBE v7.3.32 software was used to evaluate the replacement saturation of sequence data [29].

2.3. Phylogenetic Analysis

We retrieved the complete mitochondrial genomes of 14 insects from NCBI. Combined with the newly determined complete mitochondrial genomes of three species of the genus Yangiella in this study, an initial dataset containing 16 species (14 species of Aradidae, Macroscytus subaeneus, Urochela quadrinotata (Reuter 1881) as outgroup, see Table S1) was complete.
The initial dataset was imported into Phylosuite v1.2.3 software [30] to extract the nucleotide sequences of protein-coding genes (PCGs) and rRNAs. Subsequently, multiple sequence alignments of each gene of the PCG and rRNA sequences were performed using the G-INS-i strategy in MAFFT v7.464 software [31]. As MAFFT v7.464 does not consider the codon structure of PCGs, it may introduce alignment errors; therefore, it was necessary to optimize the aligned PCG sequences using MACSE [32]. Then, GBlocks v0.91 b was used to prune PCGs to remove sites with alignment errors or multiple substitutions, so as to remove phylogenetic noise and retain phylogenetic signals [33]. rRNA sequences were trimmed using trimAl v1.4 software [34]. Finally, FASconCAT-G v1.04 was used to concatenate the sequences [35]. Two datasets for phylogeny were obtained: (1) a data matrix of 13 protein-coding genes and 2 rRNA genes (PCGsRNA); (2) a data matrix (PCGs12RNA) constructed by the first and second codons of 13 protein-coding genes and 2 rRNA genes.
The overall heterogeneity of the two data matrices was evaluated using Aligroove v1.06 [36]. Then, ModelFinder [37] was used to evaluate the optimal partitioning strategy and evolutionary model of the PCGsRNA and PCGs12RNA matrices, respectively. Based on the results of ModelFinder (Tables S2–S5), a Bayesian inference (BI) phylogenetic tree was constructed using Exabayes v1.5.1 software [38]. Outgroups were manually selected, MCMC generations were set to 10,000,000, sampling frequency was 1000, four MCMC chains were run, and the Burnin Fraction value was 0.25. When the average standard deviation of the split frequencies (ASDSF) value was less than 0.01, BI ran converge and the BI phylogenetic tree was obtained. Then, IQ-tree v2 software [39] was used to construct the maximum likelihood (ML) phylogenetic tree, manually select the outgroup, and calculate the Ultra-fast Bootstrap Value as the support rate for each node [40]. The online website TVBOT (https://www.chiplot.online/tvbot.html (accessed on 10 May 2024)) was used to visualize the phylogenetic trees of BI and ML [41].

2.4. Divergence Time Estimation

For these two matrices, this study used BEAST v1.10.4 to calculate the divergence time of each genus of Aradidae [42]. The partitioning scheme was set according to the results of ModelFinder. The divergence time of the previous study was used as the calibration point [43,44] (Table S6). The prior distribution selected was the normal distribution and the relaxed molecular clock model. The total generation value of the MCMC chain was set to 20 million generations, and sampling was performed every 1000 generations. The running results were viewed using Tracer v1.7.2 to test whether the parameters converged (ESS > 200) [45], TreeAnnotator v1.10.4 was used to obtain the MCC tree (the maximum clade credibility tree), and the proportion of discarded samples was set to the top 25%. Finally, FigTree v1.4.4 was used to view the MCC tree.

3. Results

3.1. General Characteristics of the Mitochondrial Genome of the Genus Yangiella

In this study, we sequenced the complete mitochondrial genomes of three species of the genus Yangiella (Figure 1), with lengths of 15,192 bp for Y. mimetica, 15,205 bp for Y. montana, and 15,070 bp for Yangiella sp., respectively. The difference in mitochondrial genome length is mainly due to changes in the length of the control region and rRNA. The mitochondrial genomes of the three species of Yangiella contain 37 genes and a non-coding control region (CR). Among them, there are 23 J-strand genes, including 14 tRNA genes and 9 PCGs. There are 14 genes on the N-strand, including 4 PCGs, 8 tRNA genes, and 2 rRNA genes (Tables S7–S9). Compared with the mitochondrial gene arrangement of Cydnidae, tRNA gene rearrangement was found in Aradidae, forming three-gene arrangement. Among them, Aradacanthia heissi (Bai, Zhang, and Cai) and Aradus compar (Kiritshenko, 1913) each all a gene arrangement, and Aneurinae, Carventinae, and Mezirinae all have the same gene order (Figure 2).
There are 12–14 gene overlap regions in the complete mitochondrial genomes of the three species, ranging from 1 to 8 bp in length. The longest overlap region is located between trnW and trnC in the three species, which is 8 bp. There are also between seven and eight intergenic spacers (IGS) between each gene, with a length of 1–50 bp. The largest intergenic spacers of each species are located between trnQ and trnI, with a length of 44–50 bp. In addition, among the 13 protein-coding genes of each Yangiella species, except for cox1, cox2, and nad1, which used TTG as the starting codon, the remaining PCGs used ATN (N = A, T, C, G) as the starting codon. The starting codons of nad3, nad4, nad4L, and nad6 were different among different species. Ten PCGs of the three species were terminated with TAA or TAG, and cox2, nad4, and nad5 were terminated with T residues (Tables S7–S9).
The mitochondrial genome of Yangiella had AT bias, and the AT content was between 69.1 and 69.7%, showing positive AT and negative GC bias. The AT and GC bias in protein-coding genes were both negative, while the genes on the J-chain had positive AT and negative GC tilt, and the N-chain was the opposite (Table S10).
The mitochondrial genome of Yangiella contains 22 tRNAs, with a length of 61–70 bp, of which trnA is the shortest, with a length of 61 bp in all three species. Y. mimetica has the longest trnW (70 bp). The longest trnQ in Y. montana was 69 bp; the trnQ, trnK, and trnD genes of the new species were all 69 bp. The trnS1 in each species lacks the DHU arm (Dihydrouridine arm) and cannot form a typical cloverleaf secondary structure. Other tRNAs can fold to form a typical clover structure (Figure 3, Figures S1–S3). A total of 102 base mismatches occurred in the tRNA of the three species, including 75 GU mismatches, 10 UU mismatches, 8 AC mismatches, 5 UC mismatches, 3 AA mismatches, and 1 AG mismatch. GU, UU, UC, and AC mismatches were found in all three species, and the number of GU mismatches was the highest in all three species, while AG mismatches were only found in trnC of Y. montana.
The Yangiella mitochondrial genome contains two rRNA genes; rrnL is located between trnL1 and trnV, and rrnS is located between trnV and control region. The length of rrnL is between 1258 and 1384 bp, and the length of rrnS is between 743 and 748 bp, both of which have obvious AT bias (Table S10).

3.2. Codon Preference Analysis

Among the 13 PCGs of Y. mimetica and Y. montana, UUA, UCA, CGA, UCU, and ACA were used more frequently, while UUA, UCA, ACA, AGA, and CGA were used more frequently in Yangiella sp. However, the synonymous codon usage frequency of UUA was the highest in the three species, while the relative synonymous usage frequency of AGG was the lowest in the three species. The most commonly used amino acids were Ile, followed by Met, Phe, and Leu2 (Figure 4). In general, the codon usage preferences of the three species of Yangiella are similar.

3.3. Neutrality Plot Analysis and ENC-Plot Analysis

The results of neutral mapping analysis (Figure 5A) showed that the GC12 values of the mitochondrial gene of Aradidae were between 0.309 and 0.391, and the GC3 values were between 0.148 and 0.341. The GC12 and GC3 were significantly correlated (R2 = 0.9131, p < 0.01), indicating that the mutation pressure affected all codon sites. In addition, all the mitochondrial genes of Aradidae were located below the ENC standard curve (Figure 5B).

3.4. Base Substitution Saturation Analysis

The transversion rate of protein-coding genes in Aradidae is often greater than the transition rate, which is more in line with the ideal state of substitution. At the saturation level, when the total genetic distance increases, the slope of the transversion saturation curve gradually decreases, but there is no obvious platform, so it does not reach saturation; the slope of the transition saturation curve is much larger than zero, and does not reach saturation (Figure 6A). In theory, the third site of the protein-coding gene is subjected to the least selection pressure and the substitution occurs more frequently. The substitution saturation of the three coding sites was calculated, respectively. The first and second sites had almost positive slope lines, which were completely unsaturated (Figure 6B,C). The slope of the transition trend line at the third site is much larger than zero, and it also does not reach the saturation state, but the slope of the transversion trend line gradually becomes smaller, showing a slight saturation state (Figure 6D). The trend line of substitution saturation of rRNA gene is almost a standard positive slope line, which is completely unsaturated (Figure 6E).

3.5. Nucleotide Diversity (Pi) and Nonsynonymous (Ka)/Synonymous (Ks) Mutation Rate Ratios

In this study, the Pi values of 13 PCGs in the mitochondrial genome of Aradidae were calculated (Figure 7). The results showed that the overall Pi changes of 13 PCGs in Aradidae were not significant. Among them, atp8, nad6, and nad4L had higher Pi values than other coding genes, and the value of cox1 was the smallest, indicating that the gene had the least variability. The Ka/Ks ratio was used to evaluate the evolutionary rate of 13 PCGs in Aradidae species (Figure 7). The results showed that the Ka/Ks ratio of all genes was less than 1, indicating that these genes were in a state of purification selection. However, atp8, nad2, nad41, nad5, and nad6 showed higher Ka/Ks ratios, and they may have higher evolutionary rates than other genes.

3.6. Heterogeneity Analysis

Heterogeneity in nucleotide divergence was evaluated via pairwise comparisons in a multiple sequence alignment. The results showed that the heterogeneity of the two matrices was low and could be used to construct phylogenetic trees (Figure 8). On the whole, PCGs12RNA has darker blue and lower heterogeneity. The heterogeneity rate of the third codon was high, and the removal of the third codon matrix can reduce the degree of heterogeneity, which can further explain that PCG base heterogeneity is mainly concentrated in the third codon.

3.7. Phylogenetic Analyses

In this study, we conducted a phylogenetic analysis that included 14 species of Aradidae, as well as outgroup species from Pentatomidae. We used maximum likelihood (ML) and Bayesian inference (BI) methods, and used two data matrices (PCGsRNA and PCGs12RNA) to generate four phylogenetic trees. The Bayesian inference (BI) and maximum likelihood (ML) trees of the two data matrices show exactly the same topology, and most nodes are strongly supported. The results showed that except for Calisiinae and Aradinae, where only one species was not discussed, the remaining subfamilies were Aneurinae (BP = 90; PP = 1), Carventinae (BP = 90; PP = 10), and Mezirinae (BP = 100; PP = 1), which are monophyletic subfamilies. Among them, Calisiinae is located at the base of the phylogenetic tree, and Mezirinae is a highly evolved subfamily, but the ML tree support of individual points in its internal branch is weak (BP = 16–19), and it is stronger in the BI tree (PP = 0.80–0.81). In addition, the three Yangiella sequenced in this study were monophyletic (BP = 100; PP = 1), and it is a sister group to Chinolyda (Figure 9).

3.8. Divergence Time Estimation

We estimated the divergence time, and there was no significant difference between the results of the two data matrices (Table S11). Based on the results of the two data matrices, it is inferred that the most recent common ancestor of Aradidae appeared in the Late Jurassic (168 MYA). Aradinae was isolated from Carventinae, Mezirinae, and Aneurinae in the early Early Cretaceous (146 MYA). In the late Early Cretaceous (128 MYA), Carventinae was separated from Aneurinae and Mezirinae. Finally, Yangiella and Neuroctenus formed and differentiated in the early Eocene (57 MYA) (Figure 10).

4. Discussion

Gene rearrangement is considered to be an important molecular marker in revealing insect evolution [46]. The arrangement and direction of the mitochondrial genes of the three species of Yangiella were consistent with those of Aneurinae, Carventinae, and Mezirinae, but the trnC and trnY genes were shifted compared with Aradus comper in Aradinae [43]. Compared with Aradacanthia heissi [47] of Calisiinae, trnC and trnW genes were reversed. The trnQ and trnI genes were rearranged in Aradidae compared to Macroscytus subaeneus [48], which is also a Hemiptera (Figure 2). In addition, it was also found that the gene rearrangements of Aradidae mainly occurred between tRNA genes and were located between the control region and cox1. Previous studies have noted that genes around the replication start point (such as CR) are more likely to be replicated to form a “hot spot” region, making the rearrangement of gene order more likely [49]. Moreover, the same gene arrangement (trnQ–trnI) was found in the mitochondrial genomes of all Aradidae, which was different from the gene order of most other hemipteran insects (trnI–trnQ) [50]. In previous studies, this rearrangement was considered to be an ancestral feature of the Aradidae [43].
The mitochondrial genome of Yangiella insects shows high AT bias, which is widely present in the mitochondrial genome of insects. This may be because guanine (G) and cytosine (C) nucleotides require more energy for biosynthesis than adenine (A) and thymine (T) nucleotides, and are not easily available in cells [51]. Therefore, the unequal use of four bases usually leads to high A + T content [8,52]. At the same time, the two strands (heavy chain and light chain) of the mitochondrial DNA of the genus insects are not uniform, and there is no reverse chain asymmetry. The AT skew is positive and the GC skew is negative, which is consistent with the characteristics of most insect genomes [53]. This phenomenon may be related to replication and transcription mechanisms [54].
The trnS1 genes of Yangiella species all lack the (DHU) arm and form a ring. Such a phenomenon exists in most insects [8]. At the same time, six mismatched base pairs were found, of which the largest number of GU mismatches can form hydrogen bonds in RNA, but the geometric shape of the base causes the GU pair to be weaker than the hydrogen bond of the AU pair [55,56]. GU base pairing also plays an important role in biological processes [57]. In addition, non-classical UC mismatches, UU mismatches, AC mismatches, AA mismatches, and AG mismatches were also found. These mismatches may be caused by error correction mismatches during RNA editing, but the mismatch phenomenon has little effect on the corresponding function of tRNA genes [58]. If the helical regions of RNA molecules are used to construct phylogeny in subsequent studies, it is necessary to understand the way these parts of the sequence evolve. The genetic relationship between species can be judged by codon usage patterns, that is, similar codon preference between genomes means a closer genetic relationship [59,60]. In this study, the RSCU values of the mitochondrial genomes of three species of the genus Yangiella were analyzed. The analysis revealed a high degree of similarity in codon preference among the three species, further confirming that they belong to the same genus.
The ENC value of Aradidae is low, indicating that its codon preference is strong [61]. All the genes in the results of ENC mapping analysis were located under the standard curve, indicating that the mitochondrial gene codons of Aradidae insects were greatly affected by natural selection. Neutrality plot analysis showed that the correlation between GC12 and GC3 in the mitochondrial genome of Aradidae was extremely significant, indicating that the codon usage preference of Aradidae was affected by base mutation. Combined with the regression coefficient, it can be inferred that natural selection and base mutation have a great influence, but the influence of natural selection is greater.
The Pi and Ka/Ks values of cox1 are the lowest, indicating that it is the most conserved gene in the mitochondrial genome [62], while atp8 and nad6 are rapidly evolving genes. In addition, the Ka/Ks < 1 of each gene in Aradidae indicates that the evolution of Aradidae is subject to purification selection [63,64], and the stability of gene function needs to be maintained by eliminating harmful mutations [65]. This is consistent with the results of the study on the evolutionary rate of Hemiptera insects in 2015 [66]. However, in the study of scale insects, nine PCGs had a higher nonsynonymous mutation rate (Ka/Ks > 1), and the dominant evolution seemed to be positive selection [50]. This difference may be due to the adaptive evolution of scale insects with host plants, while the flat bugs mostly live under the bark of rotten wood, without specific host plants, and the living environment is rarely changed.
The results of substitution saturation level analysis showed that the nucleic acid sequence of mitochondrial protein-coding genes in Aradidae did not reach saturation, which provided data rationality support for subsequent phylogenetic analysis. Heterogeneity analysis of the two data matrices showed that the compositional heterogeneity of Aradidae was low and would not lead to abnormal topological clustering. In addition, the phylogenetic results of this study are basically consistent with the traditional morphological classification and previous molecular studies [13]. Both datasets strongly support the monophyly of each subfamily and Yangiella of Aradidae. Through the analysis of divergence time, it was found that the differentiation time of the main lineages of Aradidae was consistent with its fossil record. For example, Aradinae branched with Carventinae, Mezirinae, and Aneurinae in the Early Cretaceous, which is consistent with the age of the oldest Aradinae fossil, Aradus nicholasi (Popov, 1989) (ca. 125–113 MYA) [67]. The branch of Mezirinae in the Late Cretaceous is consistent with the time of the Mezirinae fossil Myanmezira longicornis (Heiss and Poinar Jr, 2011) (100–90 MYA) in Myanmar [68]. The results of this study can provide a scientific basis for the evolution of Aradidae insects, and more comprehensive phylogenetic and more accurate molecular clocks need further study.

5. Conclusions

In this study, we report for the first time the complete mitochondrial genomes of three species in the genus Yangiella. It was found that the mitochondria of the genus have similar structural characteristics and nucleotide composition. The unique gene arrangement of Aradidae was also found, which further proved the monophyly of this family. All the phylogenetic trees have the same topological structure, which supports the monophyly of Aneurinae, Carventinae, and Mezirinae, and the monophyly of Yangiella. Unfortunately, the number of species used in this paper is limited and was not analyzed in conjunction with morphological characters. Therefore, further and more comprehensive studies are needed, as well as exploration of species divergence times.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/insects15070533/s1, Figure S1: Secondary structure prediction of the tRNA gene of Yangiella mimetica; Figure S2: Secondary structure prediction of the tRNA gene of Yangiella montana; Figure S3: Secondary structure prediction of the tRNA gene of Yangiella sp.; Table S1: Taxonomic information and GenBank accession numbers of mitogenomes used in this study. Sequences obtained in this study are marked with *; Table S2: Optimal partitioning strategy and evolution model of PCGs12RNA BI; Table S3: Optimal partitioning strategy and evolution model of PCGs12RNA IQ; Table S4: Optimal partitioning strategy and evolution model of PCGsRNA BI; Table S5: Optimal partitioning strategy and evolution model of PCGsRNA IQ; Table S6: Calibration points and references; Table S7: Gene order and basic characteristics of the Yangiella mimetica mt genome; Table S8: Gene order and basic characteristics of the Yangiella mimetica mt genome; Table S9: Gene order and basic characteristics of the Yangiella sp mt genome; Table S10. Base composition of Y. mimetica, Y. montana, and Yangiella sp. mitogenomes; Table S11: Divergence time, 95% confidence interval in parentheses.

Author Contributions

Z.J. participated in the conception and design of the paper and completed the extraction of mitochondrial genome. L.J. was responsible for the experiments, data analysis, and thesis writing. X.B. provided specimens and completed species identification, revised the manuscript, and finally determined the publication of the version. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the National Natural Science Foundation of China (No. 32060121), with funding awarded to the corresponding author, X.B.

Data Availability Statement

The genomic data in this study are available under the NCBI accession numbers PP708566, PP545373, and PP708567.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Grebennikov, V.V.; Heiss, E. DNA barcoding of flat bugs (Hemiptera: Aradidae) with phylogenetic implications. Arthropod Syst. Phylog. 2014, 72, 213–219. [Google Scholar] [CrossRef]
  2. Kormilev, N.A.; Froeschner, R.C. Flat bugs of the world. A synonymic list (Heteroptera: Aradidae). Entomography 1987, 5, 1–246. [Google Scholar]
  3. Hsiao, T.Y. Results of the zoologico-botanical expedition to Southwest China 1955–1957. Acta Entomol. Sinica 1964, 13, 587–605. [Google Scholar]
  4. Heikkala, O.; Martikainen, P.; Kouki, J. Prescribed burning is an effective and quick method to conserve rare pyrophilous forest-dwelling flat bugs. Insect Conserv. Divers. 2017, 10, 32–41. [Google Scholar] [CrossRef]
  5. Zhang, W.; Bai, X.; Heiss, E.; Cai, W. Notes on Yangiella Hsiao (Hemiptera: Aradidae: Mezirinae), with description of a new species from China. Zootaxa 2010, 2530, 29–38. [Google Scholar] [CrossRef]
  6. Liu, Y.Q.; Li, H.; Song, F.; Zhao, Y.S.; Wilson, J.J.; Cai, W.Z. Higher-level phylogeny and evolutionary history of Pentatomomorpha (Hemiptera: Heteroptera) inferred from mitochondrial genome sequences. Syst. Entomol. 2019, 44, 810–819. [Google Scholar] [CrossRef]
  7. Li, Q.; Yang, L.X.; Xiang, D.B.; Wan, Y.; Wu, Q.; Huang, W.L.; Zhao, G. The complete mitochondrial genomes of two model ectomycorrhizal fungi (Laccaria): Features, intron dynamics and phylogenetic implications. Int. J. Biol. Macromol. 2020, 145, 974–984. [Google Scholar] [CrossRef]
  8. Cameron, S.L. Insect mitochondrial genomics: Implications for evolution and phylogeny. Annu. Rev. Entomol. 2014, 59, 95–117. [Google Scholar] [CrossRef] [PubMed]
  9. Song, S.N.; Tang, P.; Wei, S.J.; Chen, X.X. Comparative and phylogenetic analysis of the mitochondrial genomes in basal hymenopterans. Sci. Rep. 2016, 6, 20972. [Google Scholar] [CrossRef]
  10. Wolstenholme, D.R. Animal mitochondrial DNA: Structure and evolution. Int. Rev. Cytol. 1992, 141, 173–216. [Google Scholar]
  11. Zhang, D.X.; Hewitt, G.M. Insect mitochondrial control region: A review of its structure, evolution and usefulness in evolutionary studies. Biochem. Syst. Ecol. 1997, 25, 99–120. [Google Scholar] [CrossRef]
  12. Heiss, E.; Grebennikov, V. Monophyly, review, six new species and DNA barcode of micropterous Afromontane Afropictinus (Heteroptera: Aradidae). Zootaxa 2016, 4208, 149–160. [Google Scholar] [CrossRef] [PubMed]
  13. Grebennikov, V.V.; Heiss, E. Survey and DNA barcoding of flat bugs (Hemiptera: Aradidae) in the Tanzanian Forest Archipelago reveal a phylogeographically structured fauna largely unknown at the species level. Eur. J. Entomol. 2018, 115, 512–523. [Google Scholar] [CrossRef]
  14. Marchal, L.; Guilbert, E. Cladistic analysis of Aradidae (Insecta, Heteroptera) based on morphological and molecular characters. Zool. Scr. 2016, 45, 273–285. [Google Scholar] [CrossRef]
  15. Yan, Z.T.; Fan, Z.H.; He, S.L.; Wang, X.Q.; Chen, B.; Luo, S.T. Mitogenomes of Eight Nymphalidae Butterfly Species and Reconstructed Phylogeny of Nymphalidae (Nymphalidae: Lepidoptera). Genes 2023, 14, 1018. [Google Scholar] [CrossRef] [PubMed]
  16. Fenn, J.D.; Song, H.; Cameron, S.L.; Whiting, M.F. A preliminary mitochondrial genome phylogeny of Orthoptera (Insecta) and approaches to maximizing phylogenetic signal found within mitochondrial genome data. Mol. Phylogenet. Evol. 2008, 49, 59–68. [Google Scholar] [CrossRef] [PubMed]
  17. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, 884–890. [Google Scholar] [CrossRef] [PubMed]
  18. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Son, P.; Prjibelski, A.D.; et al. SPAdes: A New Genome Assembly Algorithm and Its Applications to Single-Cell Sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef]
  19. Donath, A.; Juehling, F.; Al-Arab, M.; Bernhart, S.H.; Reinhardt, F.; Stadler, P.F.; Middendorf, M.; Bernt, M. Improved annotation of protein-coding genes boundaries in metazoan mitochondrial genomes. Nucleic Acids Res. 2019, 47, 10543–10552. [Google Scholar] [CrossRef]
  20. Kearse, M.; Moir, R.; Wilson, A.; Stones-Havas, S.; Cheung, M.; Sturrock, S.; Buxton, S.; Cooper, A.; Markowitz, S.; Duran, C.; et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 2012, 28, 1647–1649. [Google Scholar] [CrossRef]
  21. Chan, P.P.; Lowe, T.M. tRNAscan-SE: Searching for tRNA Genes in Genomic Sequences. Methods Mol. Biol. 2019, 1962, 1–14. [Google Scholar] [PubMed]
  22. Grant, J.R.; Stothard, P. The CGView Server: A comparative genomics tool for circular genomes. Nucleic Acids Res. 2008, 36, W181–W184. [Google Scholar] [CrossRef] [PubMed]
  23. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef] [PubMed]
  24. Perna, N.T.; Kocher, T.D. Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes. J. Mol. Evol. 1995, 41, 353–358. [Google Scholar] [CrossRef] [PubMed]
  25. Librado, P.; Rozas, J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009, 25, 1451–1452. [Google Scholar] [CrossRef] [PubMed]
  26. Sueoka, N. Directional mutation pressure and neutral molecular evolution. Proc. Natl. Acad. Sci. USA 1988, 85, 2653–2657. [Google Scholar] [CrossRef] [PubMed]
  27. Xiao, Y.; Huang, H.X.; Chen, Y.X.; Zheng, S.; Chen, J.L.; Zou, Z.J.; Mehmood, N.; Ullah, I.; Liao, X.; Wang, J.Y. Insight on genetic features prevalent in five Ipomoea species using comparative codon pattern analysis reveals differences in major codons and reduced GC content at the 5? end of CDS. Biochem. Biophys. Res. Commun. 2023, 657, 92–99. [Google Scholar] [CrossRef] [PubMed]
  28. Novembre, J.A. Accounting for background nucleotide composition when measuring codon usage bias. Mol. Biol. Evol. 2002, 19, 1390–1394. [Google Scholar] [CrossRef]
  29. Xia, X. DAMBE7: New and Improved Tools for Data Analysis in Molecular Biology and Evolution. Mol. Biol. Evol. 2018, 35, 1550–1552. [Google Scholar] [CrossRef]
  30. Xiang, C.Y.; Gao, F.; Jakovlić, I.; Lei, H.P.; Hu, Y.; Zhang, H.; Zou, H.; Wang, G.T.; Zhang, D. Using PhyloSuite for molecular phylogeny and tree-based analyses. iMeta 2023, 2, e87. [Google Scholar] [CrossRef]
  31. Nakamura, T.; Yamada, K.D.; Tomii, K.; Katoh, K. Parallelization of MAFFT for large-scale multiple sequence alignments. Bioinformatics 2018, 34, 2490–2492. [Google Scholar] [CrossRef] [PubMed]
  32. Ranwez, V.; Douzery, E.J.P.; Cambon, C.; Chantret, N.; Delsuc, F. MACSE v2: Toolkit for the Alignment of Coding Sequences Accounting for Frameshifts and Stop Codons. Mol. Biol. Evol. 2018, 35, 2582–2584. [Google Scholar] [CrossRef] [PubMed]
  33. Castresana, J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol. Biol. Evol. 2000, 17, 540–552. [Google Scholar] [CrossRef] [PubMed]
  34. Capella-Gutiérrez, S.; Silla-Martínez, J.M.; Gabaldón, T. trimAl: A tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 2009, 25, 1972–1973. [Google Scholar] [CrossRef] [PubMed]
  35. Kueck, P.; Longo, G.C. FASconCAT-G: Extensive functions for multiple sequence alignment preparations concerning phylogenetic studies. Front. Zool. 2014, 11, 81. [Google Scholar] [CrossRef] [PubMed]
  36. Kueck, P.; Meid, S.A.; Gross, C.; Waegele, J.W.; Misof, B. AliGROOVE—Visualization of heterogeneous sequence divergence within multiple sequence alignments and detection of inflated branch support. BMC Bioinf. 2014, 15, 294. [Google Scholar] [CrossRef] [PubMed]
  37. Kalyaanamoorthy, S.; Bui Quang, M.; Wong, T.K.F.; von Haeseler, A.; Jermiin, L.S. ModelFinder: Fast model selection for accurate phylogenetic estimates. Nat. Methods 2017, 14, 586–588. [Google Scholar] [CrossRef] [PubMed]
  38. Aberer, A.J.; Kobert, K.; Stamatakis, A. ExaBayes: Massively Parallel Bayesian Tree Inference for the Whole-Genome Era. Mol. Biol. Evol. 2014, 31, 2553–2556. [Google Scholar] [CrossRef] [PubMed]
  39. Minh, B.Q.; Schmidt, H.A.; Chernomor, O.; Schrempf, D.; Woodhams, M.D.; von Haeseler, A.; Lanfear, R. IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. Mol. Biol. Evol. 2020, 37, 1530–1534. [Google Scholar] [CrossRef]
  40. Hoang, D.T.; Chernomor, O.; von Haeseler, A.; Minh, B.Q.; Vinh, L.S. UFBoot2: Improving the Ultrafast Bootstrap Approximation. Mol. Biol. Evol. 2018, 35, 518–522. [Google Scholar] [CrossRef]
  41. Xie, J.; Chen, Y.; Cai, G.; Cai, R.; Hu, Z.; Wang, H. Tree Visualization By One Table (tvBOT): A web application for visualizing, modifying and annotating phylogenetic trees. Nucleic Acids Res. 2023, 51, W587–W592. [Google Scholar] [CrossRef] [PubMed]
  42. Drummond, A.J.; Rambaut, A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol. Biol. 2007, 7, 214. [Google Scholar] [CrossRef] [PubMed]
  43. Song, F.; Li, H.; Shao, R.; Shi, A.; Bai, X.; Zheng, X.; Heiss, E.; Cai, W. Rearrangement of mitochondrial tRNA genes in flat bugs (Hemiptera: Aradidae). Sci. Rep. 2016, 6, 25725. [Google Scholar] [CrossRef] [PubMed]
  44. Lin, Q.B. Late Triassic insect fauna from Toksun, Xinjiang. Acta Pal. Sin. 1992, 31, 313–335. [Google Scholar]
  45. Rambaut, A.; Drummond, A.J.; Xie, D.; Baele, G.; Suchard, M.A. Posterior Summarization in Bayesian Phylogenetics Using Tracer 1.7. Syst. Biol. 2018, 67, 901–904. [Google Scholar] [CrossRef]
  46. Dowton, M.; Castro, L.R.; Austin, A.D. Mitochondrial gene rearrangements as phylogenetic characters in the invertebrates: The examination of genome ‘morphology’. Invertebr. Syst. 2002, 16, 345–356. [Google Scholar] [CrossRef]
  47. Shi, A.; Li, H.; Bai, X.; Dai, X.; Chang, J.; Guilbert, E.; Cai, W. The complete mitochondrial genome of the flat bug Aradacanthia heissi (Hemiptera: Aradidae). Zootaxa 2012, 3238, 23–38. [Google Scholar] [CrossRef]
  48. Xu, S.W.; Wu, Y.F.; Liu, Y.Q.; Zhao, P.; Chen, Z.; Song, F.; Li, H.; Cai, W.Z. Comparative Mitogenomics and Phylogenetic Analyses of Pentatomoidea (Hemiptera: Heteroptera). Genes 2021, 12, 1306. [Google Scholar] [CrossRef]
  49. Dowton, M.; Cameron, S.L.; Dowavic, J.I.; Austin, A.D.; Whiting, M.F. Characterization of 67 Mitochondrial tRNA Gene Rearrangements in the Hymenoptera Suggests That Mitochondrial tRNA Gene Position Is Selectively Neutral. Mol. Biol. Evol. 2009, 26, 1607–1617. [Google Scholar] [CrossRef]
  50. Lu, C.C.; Huang, X.L.; Deng, J. Mitochondrial genomes of soft scales (Hemiptera: Coccidae): Features, structures and significance. BMC Genom. 2023, 24, 37. [Google Scholar] [CrossRef]
  51. Hamilton, W.L.; Claessens, A.; Otto, T.D.; Kekre, M.; Fairhurst, R.M.; Rayner, J.C.; Kwiatkowski, D. Extreme mutation bias and high AT content in Plasmodium falciparum. Nucleic Acids Res. 2017, 45, 1889–1901. [Google Scholar] [PubMed]
  52. Simon, C.; Frati, F.; Beckenbach, A.; Crespi, B.; Liu, H.; Flook, P. Evolution, Weighting, and Phylogenetic Utility of Mitochondrial Gene Sequences and a Compilation of Conserved Polymerase Chain Reaction Primers. Ann. Entomol. Soc. Am. 1994, 87, 651–701. [Google Scholar] [CrossRef]
  53. Wei, S.J.; Shi, M.; Chen, X.X.; Sharkey, M.J.; van Achterberg, C.; Ye, G.Y.; He, J.H. New Views on Strand Asymmetry in Insect Mitochondrial Genomes. PLoS ONE 2010, 5, e12708. [Google Scholar] [CrossRef] [PubMed]
  54. Hassanin, A.; Léger, N.; Deutsch, J. Evidence for multiple reversals of asymmetric mutational constraints during the evolution of the mitochondrial genome of Metazoa, and consequences for phylogenetic inferences. Syst. Biol. 2005, 54, 277–298. [Google Scholar] [CrossRef] [PubMed]
  55. Rousset, F.; Pelandakis, M.; Solignac, M. Evolution of compensatory substitutions through G.U intermediate state in Drosophila rRNA. Proc. Natl. Acad. Sci. USA 1991, 88, 10032–10036. [Google Scholar] [CrossRef] [PubMed]
  56. Golden, M.; Murrell, B.; Martin, D.; Pybus, O.G.; Hein, J. Evolutionary Analyses of Base-Pairing Interactions in DNA and RNA Secondary Structures. Mol. Biol. Evol. 2020, 37, 576–592. [Google Scholar] [CrossRef] [PubMed]
  57. Varani, G.; McClain, W.H. The G × U wobble base pair. A fundamental building block of RNA structure crucial to RNA function in diverse biological systems. EMBO Rep. 2000, 1, 18–23. [Google Scholar] [CrossRef] [PubMed]
  58. Yokobori, S.; Paabo, S. Transfer RNA editing in land snail mitochondria. Proc. Natl. Acad. Sci. USA 1995, 92, 10432–10435. [Google Scholar] [CrossRef] [PubMed]
  59. Chakraborty, S.; Yengkhom, S.; Uddin, A. Analysis of codon usage bias of chloroplast genes in Oryza species Codon usage of chloroplast genes in Oryza species. Planta 2020, 252, 67. [Google Scholar] [CrossRef] [PubMed]
  60. Bu, Y.F.; Wu, B.X.; Sun, N.; Man, Y.; Jing, Y.P. Codon usage bias predicts the functional MYB10 gene in Populus. J. Plant Physiol. 2021, 265, 153491. [Google Scholar] [CrossRef]
  61. Chakraborty, S.; Nag, D.; Mazumder, T.H.; Uddin, A. Codon usage pattern and prediction of gene expression level in Bungarus species. Gene 2017, 604, 48–60. [Google Scholar] [CrossRef] [PubMed]
  62. Demari-Silva, B.; Foster, P.G.; de Oliveira, T.M.P.; Bergo, E.S.; Sanabani, S.S.; Pessôa, R.; Sallum, M.A.M. Mitochondrial genomes and comparative analyses of Culex camposi, Culex coronator, Culex usquatus and Culex usquatissimus (Diptera: Culicidae), members of the coronator group. BMC Genom. 2015, 16, 831. [Google Scholar] [CrossRef] [PubMed]
  63. Hurst, L.D. The Ka/Ks ratio:: Diagnosing the form of sequence evolution. Trends Genet. 2002, 18, 486–487. [Google Scholar] [CrossRef] [PubMed]
  64. Meiklejohn, C.D.; Montooth, K.L.; Rand, D.M. Positive and negative selection on the mitochondrial genome. Trends Genet. 2007, 23, 259–263. [Google Scholar] [CrossRef] [PubMed]
  65. Chang, H.H.; Qiu, Z.Y.; Yuan, H.; Wang, X.Y.; Li, X.J.; Sun, H.M.; Guo, X.Q.; Lu, Y.C.; Feng, X.L.; Majid, M.; et al. Evolutionary rates of and selective constraints on the mitochondrial genomes of Orthoptera insects with different wing types. Mol. Phylogenet. Evol. 2020, 145, 106734. [Google Scholar] [CrossRef] [PubMed]
  66. Wang, Y.; Chen, J.; Jiang, L.Y.; Qiao, G.X. Hemipteran Mitochondrial Genomes: Features, Structures and Implications for Phylogeny. Int. J. Mol. Sci. 2015, 16, 12382–12404. [Google Scholar] [CrossRef]
  67. Popov, Y.A. New fossil Hemiptera (Heteroptera and Coleorrhyncha) from the Mesozoic of Mongolia. Neues Jahrb. Für Geol. Und Paläontologie. Monatshefte 1989, 3, 166–181. [Google Scholar] [CrossRef]
  68. Heiss, E.; Grimaldi, D.A. New Aradidae in Mesozoic Burmese Amber (Hemiptera, Heteroptera). Ann. Naturhist. Mus. Wien Serie A 2012, 114, 307–316. [Google Scholar]
Insects 15 00533 g001
Figure 1. Mitochondrial genome maps of three species of Yangiella. The chain is marked with an arrow indicating the direction of gene transcription. Gene lengths correspond to nucleotide lengths in the diagram. The outermost layer is the J-chain, and the second layer is the N-chain.
Figure 1. Mitochondrial genome maps of three species of Yangiella. The chain is marked with an arrow indicating the direction of gene transcription. Gene lengths correspond to nucleotide lengths in the diagram. The outermost layer is the J-chain, and the second layer is the N-chain.
Insects 15 00533 g001
Insects 15 00533 g002
Figure 2. Mitochondrial gene rearrangements in flat bugs. Except for the underlined genes, the other genes were transcribed from left to right. IGS is the intergenic region. (A) Macroscytus subaeneus. (B) Aradacanthia heissi. (C) Aradus compar. (D) Aneurus similis, Aneurus sublobatus, Libiocoris heissi, Taiwanaptera montana, Arbanatus sp., Brachyrhynchus hsiaoi, Brachyrhynchus triangulus, Mezira sp., Neuroctenus yunnanensis, Yangiella sp., Yangiella mimetica, Yangiella montana.
Figure 2. Mitochondrial gene rearrangements in flat bugs. Except for the underlined genes, the other genes were transcribed from left to right. IGS is the intergenic region. (A) Macroscytus subaeneus. (B) Aradacanthia heissi. (C) Aradus compar. (D) Aneurus similis, Aneurus sublobatus, Libiocoris heissi, Taiwanaptera montana, Arbanatus sp., Brachyrhynchus hsiaoi, Brachyrhynchus triangulus, Mezira sp., Neuroctenus yunnanensis, Yangiella sp., Yangiella mimetica, Yangiella montana.
Insects 15 00533 g002
Insects 15 00533 g003
Figure 3. Secondary structure prediction of the tRNA gene of Yangiella sp. The Watson–Crick pair is represented by a straight line, and the mismatch is represented by a dot. All of the Yangiella sp. mitogenome tRNAs except for trnS1 follow the canonical conserved three-loop cloverleaf structure consisting of the amino-acyl arm (lilac), dihydrouridine arm (DHU, green), pseudouridine arm (TΨU loop, blue), and the anticodon loop (modena). All tRNAs also contain a variable region (red).
Figure 3. Secondary structure prediction of the tRNA gene of Yangiella sp. The Watson–Crick pair is represented by a straight line, and the mismatch is represented by a dot. All of the Yangiella sp. mitogenome tRNAs except for trnS1 follow the canonical conserved three-loop cloverleaf structure consisting of the amino-acyl arm (lilac), dihydrouridine arm (DHU, green), pseudouridine arm (TΨU loop, blue), and the anticodon loop (modena). All tRNAs also contain a variable region (red).
Insects 15 00533 g003
Insects 15 00533 g004
Figure 4. Relative synonymous codon usage (RSCU) in the mitogenomes of the three species. The numbers above the bar graph indicate the frequency of amino acids. The number of codons per amino acid varies from 2 to 4. The RSCU values are color-coded based on the codons below the amino acid labels.
Figure 4. Relative synonymous codon usage (RSCU) in the mitogenomes of the three species. The numbers above the bar graph indicate the frequency of amino acids. The number of codons per amino acid varies from 2 to 4. The RSCU values are color-coded based on the codons below the amino acid labels.
Insects 15 00533 g004
Insects 15 00533 g005
Figure 5. Codon preference analysis of 14 species in Aradidae. (A) The G + C content trend line and R value of codon 1 and 2 (GC12) and codon 3 (GC3) of mitochondrial protein-coding gene in Aradidae. (B) The scatter plot of the correlation between ENC value and GC3. The purple curve represents the expected functional relationship between ENC and GC3 under no selection pressure but only under abrupt pressure.
Figure 5. Codon preference analysis of 14 species in Aradidae. (A) The G + C content trend line and R value of codon 1 and 2 (GC12) and codon 3 (GC3) of mitochondrial protein-coding gene in Aradidae. (B) The scatter plot of the correlation between ENC value and GC3. The purple curve represents the expected functional relationship between ENC and GC3 under no selection pressure but only under abrupt pressure.
Insects 15 00533 g005
Insects 15 00533 g006
Figure 6. Substitution saturation plots of PCGs and RNA in the mt genomes of 14 species of Aradidae. “S” represents the transition rate (dark red), and “V” represents the transversion rate (dark green). (A) All sites of protein-coding genes. (B) The first site of protein-coding genes. (C) Protein-coding gene second site. (D) Protein-coding gene third site. (E) rRNA gene.
Figure 6. Substitution saturation plots of PCGs and RNA in the mt genomes of 14 species of Aradidae. “S” represents the transition rate (dark red), and “V” represents the transversion rate (dark green). (A) All sites of protein-coding genes. (B) The first site of protein-coding genes. (C) Protein-coding gene second site. (D) Protein-coding gene third site. (E) rRNA gene.
Insects 15 00533 g006
Insects 15 00533 g007
Figure 7. Nucleotide diversity and Ka/Ks values of PCGs in mitochondrial genomes of 14 species of Aradidae.
Figure 7. Nucleotide diversity and Ka/Ks values of PCGs in mitochondrial genomes of 14 species of Aradidae.
Insects 15 00533 g007
Insects 15 00533 g008
Figure 8. Heterogeneity of the sequence composition of the mitochondrial genomes in different datasets. The pairwise Aliscore values are indicated by colored squares. Darker colors indicate full random similarity, and lighter colors indicate the opposite. (A) PCGsRNA matrix. (B) PCGs12RNA matrix.
Figure 8. Heterogeneity of the sequence composition of the mitochondrial genomes in different datasets. The pairwise Aliscore values are indicated by colored squares. Darker colors indicate full random similarity, and lighter colors indicate the opposite. (A) PCGsRNA matrix. (B) PCGs12RNA matrix.
Insects 15 00533 g008
Insects 15 00533 g009
Figure 9. Phylogenetic tree of 16 species of Hemiptera inferred from the PCGsRNA and PCGs12RNA data matrix using maximum likelihood (ML) and Bayesian inference (BI). The ML tree has the same topology as the BI tree, and their support values are reported above and below the nodes, respectively. The species marked as * is the research object of this study. (A) Phylogenetic tree based on PCGsRNA data matrix. (B) Phylogenetic tree based on PCGs12RNA data matrix.
Figure 9. Phylogenetic tree of 16 species of Hemiptera inferred from the PCGsRNA and PCGs12RNA data matrix using maximum likelihood (ML) and Bayesian inference (BI). The ML tree has the same topology as the BI tree, and their support values are reported above and below the nodes, respectively. The species marked as * is the research object of this study. (A) Phylogenetic tree based on PCGsRNA data matrix. (B) Phylogenetic tree based on PCGs12RNA data matrix.
Insects 15 00533 g009
Insects 15 00533 g010
Figure 10. Chronogram showing Aradidae phylogeny and divergence time. The consensus tree shows the results of the divergence time estimation of the PCGsRNA dataset using three calibration points (Table S6). The species marked as * is the research object of this study. A geological time scale is shown at the bottom.
Figure 10. Chronogram showing Aradidae phylogeny and divergence time. The consensus tree shows the results of the divergence time estimation of the PCGsRNA dataset using three calibration points (Table S6). The species marked as * is the research object of this study. A geological time scale is shown at the bottom.
Insects 15 00533 g010
Table 1. Voucher information of the specimens used for mitochondrial genome sequencing.
Table 1. Voucher information of the specimens used for mitochondrial genome sequencing.
SpecimensDate of CollectionCollection SiteLongitude (E)Latitude (N)GB Numbers
Y. mimetica16 January 2011Jinping, Yunnan103.238322.9067PP545373
Y. montana31 October 2019Yingjiang Tongbiguan, Yunnan97.936424.7031PP708566
Yangiella sp.27 May 2023Shuangbai Dutian Forest farm, Yunnan101.455624.5447PP708567
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ji, L.; Jia, Z.; Bai, X. Comparative Analysis of the Mitochondrial Genomes of Three Species of Yangiella (Hemiptera: Aradidae) and the Phylogenetic Implications of Aradidae. Insects 2024, 15, 533. https://doi.org/10.3390/insects15070533

AMA Style

Ji L, Jia Z, Bai X. Comparative Analysis of the Mitochondrial Genomes of Three Species of Yangiella (Hemiptera: Aradidae) and the Phylogenetic Implications of Aradidae. Insects. 2024; 15(7):533. https://doi.org/10.3390/insects15070533

Chicago/Turabian Style

Ji, Liangpeng, Zhancheng Jia, and Xiaoshuan Bai. 2024. "Comparative Analysis of the Mitochondrial Genomes of Three Species of Yangiella (Hemiptera: Aradidae) and the Phylogenetic Implications of Aradidae" Insects 15, no. 7: 533. https://doi.org/10.3390/insects15070533

APA Style

Ji, L., Jia, Z., & Bai, X. (2024). Comparative Analysis of the Mitochondrial Genomes of Three Species of Yangiella (Hemiptera: Aradidae) and the Phylogenetic Implications of Aradidae. Insects, 15(7), 533. https://doi.org/10.3390/insects15070533

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop