Next Article in Journal
Sensitivity of Vanessa cardui to Temperature Variations: A Cost-Effective Experiment for Environmental Education
Previous Article in Journal
Seasonal Variability in the Prevalence of DWV Strains in Individual Colonies of European Honeybees in Hawaii
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Complete Mitochondrial Genome of the Eggplant Fruit and Shoot Borer, Leucinodes orbonalis Guenée (Lepidoptera: Crambidae), and Comparison with Other Pyraloid Moths

by
Joshua B. Despabiladeras
and
Ma. Anita M. Bautista
*
Functional Genomics Laboratory, National Institute of Molecular Biology and Biotechnology, College of Science, University of the Philippines-Diliman, Quezon City 1101, Philippines
*
Author to whom correspondence should be addressed.
Submission received: 18 February 2024 / Revised: 17 March 2024 / Accepted: 18 March 2024 / Published: 25 March 2024
(This article belongs to the Section Insect Molecular Biology and Genomics)

Abstract

:

Simple Summary

Eggplant is an important agricultural produce in the Philippines, but consistent production is hampered by frequent infestations of the eggplant fruit and shoot borer (EFSB). The management and monitoring of EFSB infestation requires full understanding of its biology, including its evolution and adaptations. The use of genomic resources provides a starting point for understanding the biology of EFSB. However, these resources are lacking for EFSB. This study aimed to sequence and characterize the mitochondrial genome (mitogenome) of the EFSB and use the mitogenome in reconstructing the evolutionary relationships of the EFSB with respect to other members of the economically important superfamily, Pyraloidea. The annotated mitogenome of the EFSB contains the typical elements expected for a lepidopteran mitochondrion, such as gene content, gene syntenies, and nucleotide composition. Phylogenetic analyses using 72 mitogenomes from other moths robustly placed the EFSB as basal to its supposed subfamily, Spilomelinae, despite recovering the expected subfamily relationships of earlier works. The mitogenome resource for the EFSB generated here will add to the scarce genomic resources for the EFSB and be useful for future phylogenetic reconstruction of Pyraloidea and Spilomelinae.

Abstract

The eggplant fruit and shoot borer (EFSB) (Leucinodes orbonalis Guenée) is a devastating lepidopteran pest of eggplant (Solanum melongena L.) in the Philippines. Management of an insect pest like the EFSB requires an understanding of its biology, evolution, and adaptations. Genomic resources provide a starting point for understanding EFSB biology, as the resources can be used for phylogenetics and population structure studies. To date, genomic resources are scarce for EFSB; thus, this study generated its complete mitochondrial genome (mitogenome). The circular mitogenome is 15,244 bp-long. It contains 37 genes, namely 13 protein-coding, 22 tRNA, and 2 rRNA genes, and has conserved noncoding regions, motifs, and gene syntenies characteristic of lepidopteran mitogenomes. Some protein-coding genes start and end with non-canonical codons. The tRNA genes exhibit a conserved cloverleaf structure, with the exception in trnS1. Partitioned phylogenetic analysis using 72 pyraloids generated highly supported maximum likelihood and Bayesian inference trees revealing expected basal splits between Crambidae and Pyralidae, and Spilomelinae and Pyraustinae. Spilomelinae was recovered to be paraphyletic, with the EFSB robustly placed before the split of Spilomelinae and Pyraustinae. Overall, the EFSB mitogenome resource will be useful for delineations within Spilomelinae and population structure analysis.

1. Introduction

The mitochondrion is a fundamental, energy-generating organelle responsible for oxidative phosphorylation [1,2]. Mitochondria are believed to have originated from the endosymbiosis of an ancestral Archaeal prokaryote with an alphaproteobacterium [3]. The energy-generating capacity of the mitochondria is fundamental to life that almost all extant eukaryotes contain fully functional or derivative forms of mitochondria [4], with the only exception to date being the flagellated excavates, Monocercomonoides [5]. Although most of the mitochondrial proteins are encoded by the nuclear genome, most mitochondria preserve a highly reduced genome encoding 13 protein-coding genes important for the electron transport chain and respiratory function. Mitochondrial genomes also contain a self-sufficient translation machinery consisting of a complete set of 22 tRNAs and 2 rRNAs [6]. Some mitochondrial genes (e.g., cytochrome oxidase I (COI) subunit) are already widely established for use in phylogenetics [7]. Animal mitogenomes have conserved gene content and maternal inheritance, properties that enable extensive use not only in phylogenetics but also in population genetics [8], diagnostics and species delineation [9], and molecular evolution [10,11]. Moreover, the use of multi-gene datasets in phylogenetics has additional advantages such as reducing sampling error, making it the modern standard in inferring species relationships between organisms [12,13].
The eggplant fruit and shoot borer (EFSB), Leucinodes orbonalis Guenée, is a devastating crambid pest of the eggplant, Solanum melongena L. The eggplant is one of the top agricultural products of the Philippines, but consistent production is hampered by frequent EFSB infestations [14]. EFSB is a member of the species-rich lepidopteran superfamily Pyraloidea, which contains at least 15,500 genera [15]. The superfamily is ecologically diverse and economically significant, having numerous major crop pests [16]. Within Pyraloidea, the EFSB is a member of the largest subfamily, Spilomelinae [17]. Obtaining reference mitogenomes for this superfamily can aid in pest monitoring and management schemes since mitochondrial genes have been used for population structure determination [18], and the development of possible control strategies such as the Trojan female [19]. Additionally, as previously mentioned, mitogenomes are useful in establishing evolutionary relationships between organisms. Reliable phylogenies are an important tool that could help understand, for instance, how host-plant ranges evolve within the family [16]. Systematics on Spilomelinae has been historically difficult due to morphological heterogeneity, but significant progress has been made with the help of mitochondrial and nuclear markers [17,20]. The use of whole mitogenomes could help resolve the systematics of Spilomelinae. However, as of 13 March 2024, only 23 whole mitochondrial sequences of Spilomelinae are on the RefSeq database. Thus, adding more mitogenome sequences and representatives could improve resolution and taxonomic delineation in this economically important subfamily.
The continued development of high-throughput sequencing (HTS) technologies has increased accessibility for most laboratories, leading to significant developments in phylogenetics, population genetics, barcoding, and biodiversity studies. Philippine EFSB research could also benefit from genomic resources produced by HTS, which can significantly impact the implementation of pest management. Genetic resources (mitochondrial COI) and morphological traits have been used by Sagarbarria et al. (2018) to study the population structure of EFSB in the Philippines [18]. The authors concluded that a single, widespread haplotype predominates Philippine EFSB populations. However, the authors are also aware of the limitations of using COI as the sole genetic marker for population genetic analysis and have recommended the use of other markers. The analysis of whole mitogenomes is attractive since they represent a compromise between single markers and whole genome analysis. The conserved properties of animal mitogenomes, the advantage of using multi-gene datasets, and the continued ease of acquisition would provide a promising strategy for species identification and population genetics using whole mitogenomes [21]. Whole mitogenomes were successfully used to determine the population structure of the Chinese seabass (Lateolabrax maculatus) in China [22] and Mansonia mosquitoes in Brazil [23]. Here, we present a report on the fully annotated sequence of the EFSB mitochondrial genome from short Illumina sequencing reads and phylogenetic reconstruction of Pyraloidea using whole mitogenomes. The EFSB mitogenome presented here will contribute to the genetic resources of EFSB and be useful in future population studies of EFSB in the Philippines.

2. Materials and Methods

2.1. Sample Collection and Genomic DNA Extraction

Eggplants obtained from farmers from different parts of Luzon, Philippines (i.e., Tarlac, Ilocos Norte, Laguna) were individually dissected and probed for EFSB infestations. The Institute of Plant Breeding (IPB) at the College of Agriculture and Food Science, University of the Philippines Los Baños, also provided samples. EFSB samples were stored at −80 °C until extraction.
Extraction of EFSB genomic DNA was performed using the Qiagen MagAttract HMW DNA kit (Qiagen Inc., Mckinley West Park, Taguig, Philippines). The protocol described by the manufacturer with some modifications recommended by Kingan et al. (2019) was followed: proteinase K and buffer ATL were pre-mixed, and the resulting solution was the one used for micropestle homogenization; proteinase K digestion was shortened to 2 h; all mixing steps were performed via gentle flicking; RNase A incubation time was increased to 5 min; wide-bore tips were used to transfer all solutions containing nucleic acids; and all centrifugation steps were performed at 4 °C [24]. Extracts were run on a 0.7% agarose gel for 45 min at 70 V and visualized with GelRed (Biotium Inc., Fremont, CA, USA) staining. Extracts were sent to the Philippine Genome Center for Pulse Field Gel Electrophoresis (PFGE) analysis.

2.2. Library Preparation and Sequencing

Suitable extracts determined by PFGE were chosen as the starting material for Illumina library preparation using the TruSeq DNA PCR-Free library prep kit (Illumina Inc., Biopolis Way, Singapore), with the IDT for Illumina TruSeq DNA UD indexes plate (Illumina Inc., Biopolis Way, Singapore) serving as the index adapters. Libraries were constructed, and quality checked using Agilent TapeStation (Agilent Technologies Inc., Yishun Ave. 7, Singapore). Libraries that were determined to have a 550 bp insert size were then sequenced on an Illumina NovaSeq 6000 system (Illumina, Biopolis Way, Singapore) using a 2 × 150 bp run for 300 cycles at the Philippine Genome Center.

2.3. Pre-Processing, Read Filtering, Mitogenome Assembly, and Annotation

Raw sequences were pre-processed using fastp 0.23.2 [25] to remove adapter and polyG sequences and assessed using FastQC 0.11.9 [26]. The processed sequences were then mapped to reference mitochondrial sequences using Bowtie2 2.3.5.1 [27]. Mitochondrial references were obtained from the NCBI organelle genome database and chosen based on the following criteria: close phylogenetic relationship with EFSB, tagged as a complete reference by NCBI, and must be reported or part of a published material. Following said criteria, the following mitochondrial genomes were chosen: the Glyphodes quadrimaculalis mitogenome, NC_022699.1 [28]; Maruca testulalis mitogenome, NC_024283.1 [29]; and Omiodes indicata mitogenome, NC_039177.1 [30]. Mapped sequences were then extracted and used as the starting input for assembly using metaSPAdes 3.15.4 [31] packaged within MitoZ 3.4 [32], with k-mer lengths of 21, 33, 55, and 77. Downstream annotation of the mitochondrial genes was performed using the annotation module of MitoZ. An additional annotation was performed using MITOS 2.0.8 [33] to confirm the annotations performed by MitoZ, determine the position of the control region, and determine the structures of the tRNAs and rRNAs. Secondary structures were verified using tRNAscan-SE 2.0.9 [34], and subsequently drawn using the StructureEditor 1.0 packaged within the RNAstructure 6.4 software package [35]. The final annotations were also manually checked to determine if the gene boundaries were consistent. Coverage analysis was performed by mapping the reads to the assembled mitogenome using bwa 0.7.17-r1188 [36] and assessing the coverage using samtools 1.16 [37] and mosdepth 0.3.2 [38]. The final mitogenome of the EFSB has been deposited and is available in GenBank under accession number PP493058.

2.4. Nucleotide Sequence Composition Analysis

The nucleotide composition, codon usage, and relative synonymous codon usage (RSCU) of the EFSB mitogenome were calculated using Phylosuite 1.2.3 [39]. GC skew analysis using 100 bp windows and a 20 bp step size was performed using SkewIT [40]. GC content in 20 bp sliding windows and sequencing depth per position was calculated within MitoZ. A Circos plot summarizing gene content and organization, GC content, GC skew, and coverage across all positions was created using Circos v0.69-8 [41].

2.5. Multiple Sequence Alignment and Phylogenetic Analysis

Reference sequences of representative members of Pyraloidea were obtained from the NCBI with the following search parameters: “Pyraloidea[ORGN] AND (mitochondrion[TITL] or mitochondria[TITL]) AND 10,000:20,000[SLEN]”. The downloaded Genbank files were imported to Phylosuite [39] for sequence cleaning and standardization. Only sequences flagged as “NC” were retained for analysis. The filtered sequences were then manually checked for ambiguous bases (R, Y, N, etc.), and sequences containing such bases were removed. A total of 72 working mitogenome sequences from Pyraloidea were retained after filtering and standardization. The following subfamilies were represented: Spilomelinae, Pyraustinae, Odontinae, Acentropinae, Schoenobiinae, Glaphyriinae, and Crambinae of Crambidae; and Gallerinae, Phycitinae, Pyralinae, and Epipaschiinae of Pyralidae. The full list of organisms used for the alignment is shown in Table 1. Each of the 37 genes was extracted from the reference sequences and aligned separately using MAFFT 7.505 [42], and uninformative sites were removed using Gblocks 0.91b [43]. The cleaned alignments were then concatenated to form a linear sequence of 13,953 bp for each organism using Phylosuite [39]. Partitioned maximum likelihood and Bayesian inference phylogenetic analyses were performed using IQ-tree 2.1.2 [44,45] and MrBayes 3.2.7 [46]. For maximum likelihood phylogenetic analysis, each of the 37 genes forms 1 partition initially. Automatic model selection was performed using ModelFinder in IQ-tree, and gene partitions were automatically merged to avoid over-parametrization [47]. Branch support was assessed using the SH-aLRT test with 10,000 replicates. For Bayesian inference, Metropolis-Coupled, Markov chain Monte Carlo (MC3) sampling was performed to estimate the tree with the highest posterior probability. The parallel version of MrBayes [46] was used alongside the BEAGLE library [48] to speed up calculations. In brief, ModelFinder was used to find the best-of-fit models, with the -mset parameter to restrict the search to models utilized by MrBayes. Partitions were then merged, and the merged partitions and their respective prior rates, gamma shape parameter, and proportion of invariant sites were then used as inputs to MrBayes. A total of three independent MCMC runs were performed, with each run containing 1 cold and 5 hot chains. The MCMC runs were run for 10,000,000 generations, with sampling every 1000 generations. Convergence was assessed when the effective sample size (ESS) for the parameters was greater than 200, and the potential scale reduction factor (PSRF) approached 1. The first 25% of the trees were discarded as burn-in, and the BI consensus trees were built using a strict consensus rule (Allcompat). The final consensus trees for both ML and BI were then drawn and annotated using ggtree [49]. Topological concordance was assessed by plotting both trees side by side and comparing the tips.

3. Results and Discussion

3.1. EFSB Mitogenome Structure and Organization

A total of 975,248 sequences were obtained after Bowtie2 mapping, and FastQC analysis revealed that the sequences are of high quality (Figure S1), with most of the sequences having a Phred score of 30 or above. The cleaned sequences were then used for mitogenome assembly using metaSPAdes and subsequent annotation via MitoZ and MITOS. The sequencing data obtained here yielded a minimum coverage depth value of 2215× and a maximum coverage depth value of 19,989×, yielding a mean coverage depth across all positions of 9571×. Cumulative coverage distribution obtained using mosdepth shows that 100% of the EFSB mitogenome has been covered at least 2215 times (Figure S2). The assembled mitochondrial genome has a total length of 15,244 bp and contains the 37 genes and control region expected for lepidopteran mitochondria (Figure 1).
The EFSB mitogenome contains 13 protein-coding genes, 22 tRNA genes, and 2 rRNA genes. Of the 37 genes, 23 genes (14 tRNAs and 9 PCGs) are found in the majority strand, while 14 genes (8 tRNAs, 4 PCGs, and 2 rRNAs) are found in the minority strand. The obtained mitogenome shows the typical circular double-stranded DNA structure and exhibits the conserved lepidopteran synteny of trnM-trnI-trnQ-nad2 similarly to other crambid mitogenomes [50,51]. This synteny is characteristic of lepidopterans and differs from the ancestral insect mitogenome synteny of trnI-trnQ-trnM-nad2 [52]. The EFSB mitogenome is highly AT-rich, having an AT% of 80.9%. Lepidopterans consistently have AT-rich mitogenomes with negative GC skew values (Table 1). The EFSB mitogenome is consistent with this observation, having a slightly negative AT skew of −0.01 and a negative GC-skew of −0.17. The EFSB mitogenome contains a control region with high AT content and is situated between the s-rna and trnM genes (Figure 1A,B). The control regions are also characterized by the switch in sign in the GC-skew diagrams (Figure 1C). The control region contains conserved motifs necessary for mitochondrial replication (see below) [53].
Table 1. Reference mitochondrial genomes chosen for composition comparison and phylogenetic analysis.
Table 1. Reference mitochondrial genomes chosen for composition comparison and phylogenetic analysis.
OrganismRefSeq IDTaxonomyAT%AT-SkewGC%GC-SkewReferences
Leucinodes orbonalisTHIS STUDYCrambidae, Spilomelinae80.9−0.01319.1−0.172THIS STUDY
Acrobasis inoueiNC_061244.1Pyralidae, Phycitinae80.3−0.02119.7−0.219[54]
Aedes albopictusNC_006817.1Diptera79.60.00820.5−0.181Direct submission
Aglossa dimidiataNC_058009.1Pyralidae, Pyralinae 79.1−0.04320.9−0.226Direct submission
Amyelois transitellaNC_028443.1Pyralidae, Phycitinae79.6−0.04820.4−0.237Direct submission
Botyodes diniasalisNC_073002.1Crambidae, Spilomelinae80.8−0.02119.1−0.182Direct submission
Botyodes principalisNC_061248.1Crambidae, Spilomelinae80.7−0.01419.3−0.189[54]
Cataclysta lemnataNC_050323.1Crambidae, Acentropinae79.50.00420.6−0.214Direct submission
Cathayia obliquellaNC_053657.1Pyralidae, Gallerinae80.6−0.03719.4−0.225[55]
Chilo sacchariphagusNC_029716.1Crambidae, Crambinae81−0.01219.1−0.246Direct submission
Chilo suppressalisNC_015612.1Crambidae, Crambinae80.60.00819.3−0.235[56]
Cnaphalocrocis
medinalis
NC_015985.1Crambidae, Pyraustinae82−0.01518−0.175[56]
Cnaphalocrocis
patnalis
NC_060868.1Crambidae, Pyraustinae81.8−0.02218.2−0.165Direct submission
Conogethes
punctiferalis
NC_021389.1Crambidae, Spilomelinae80.6−0.02519.4−0.207[57]
OrganismRefSeq IDTaxonomyAT%AT-SkewGC%GC-SkewReferences
Crambus perlellusNC_061606.1Crambidae, Crambinae81.3−0.00818.7−0.197Direct submission
Culex
quinquefasciatus
NC_014574.1Diptera780.00722−0.173[58]
Cydalima
perspectalis
NC_042150.1Crambidae, Spilomelinae80.9−0.01619.1−0.193[59]
Dausara
latiterminalis
NC_056799.1Crambidae, Odontinae 80.5−0.00319.5−0.201[60]
Diatraea saccharalisNC_013274.1Crambidae, Crambinae80.10.02120−0.258[61]
Dioryctria rubellaNC_061242.1Pyralidae, Phycitinae79.8−0.02320.1−0.233[54]
Drosophila
melanogaster
NC_024511.2Diptera82.20.01617.8−0.147Direct submission
Dusungwua
basinigra
NC_061240.1Pyralidae, Phycitinae80−0.01420−0.213[54]
Elophila
interruptalis
NC_021756.1Crambidae, Acentropinae80.3−0.01119.7−0.229[50]
Elophila turbataNC_068592.1Crambidae, Acentropinae81.2−0.00518.9−0.225Direct submission
Endotricha consociaNC_037501.1Pyralidae, Pyralinae79.7−0.03920.2−0.226[62]
Endotricha
kuznetzovi
NC_061642.1Pyralidae, Pyralinae80.7−0.03319.2−0.206Direct submission
Ephestia elutellaNC_039716.1Pyralidae, Phycitinae80.7−0.04319.4−0.217[63]
Ephestia kuehniellaNC_022476.1Pyralidae, Phycitinae79.7−0.04920.2−0.234[64]
Evergestis extimalisNC_071781.1Crambidae, Glaphyriinae80.7−0.01819.3−0.168Direct submission
Evergestis junctalisNC_030509.1Crambidae, Glaphyriinae81−0.01519−0.168Direct submission
Galleria mellonellaNC_028532.1Pyralidae, Gallerinae80.4−0.03919.6−0.237Direct submission
Glyphodes pyloalisNC_025933.1Crambidae, Spilomelinae80.7−0.01619.3−0.194Direct submission
Glyphodes
quadrimaculalis
NC_022699.1Crambidae, Spilomelinae80.8−0.00719.2−0.192[28]
Heortia vitessoidesNC_056800.1Crambidae, Odontinae80.6−0.01219.4−0.172[60]
Hypsopygia reginaNC_030508.1Pyralidae, Pyralinae78.7−0.03721.3−0.228Direct submission
Lamoria adaptellaNC_062173.1Pyralidae, Gallerinae80.1−0.01219.9−0.24Direct submission
Lista haraldusalisNC_024535.1Pyralidae, Epipaschiinae81.5−0.00718.5−0.171[65]
Loxostege sticticalisNC_027174.1Crambidae, Pyraustinae80.80.00219.2−0.191Direct submission
Maruca testulalisNC_024283.1Crambidae, Spilomelinae80.8−0.00519.2−0.171[29]
Maruca vitrataNC_024099.1Crambidae, Spilomelinae80.7−0.00219.3−0.172Direct submission
Meroptera pravellaNC_035242.1Pyralidae, Phycitinae80.5−0.01919.3−0.199[66]
Nagiella inferiorNC_040973.1Crambidae, Spilomelinae81.50.00918.5−0.22Direct submission
Nomophila noctuellaNC_025764.1Crambidae, Spilomelinae81.40.00218.6−0.176[67]
Omiodes indicataNC_039177.1Crambidae, Spilomelinae81.6−0.01218.4−0.162[30]
Omphisa
fuscidentalis
NC_066444.1Crambidae, Spilomelinae790.01321−0.274Direct submission
Orthaga euadrusalisNC_061246.1Pyralidae, Epipaschiinae80.2−0.0219.8−0.199[54]
OrganismRefSeq IDTaxonomyAT%AT-SkewGC%GC-SkewReferences
Orthaga olivaceaNC_046504.1Pyralidae, Epipaschiinae79−0.04321−0.215Direct submission
Orthopygia
glaucinalis
NC_047304.1Pyralidae, Pyralinae79.2−0.04420.8−0.223[68]
Orybina regalisNC_061247.1Pyralidae, Pyralinae81−0.01619−0.205[54]
Ostrinia furnacalisNC_056248.1Crambidae, Pyraustinae80.90.03119.1−0.196[69]
Ostrinia kasmiricaNC_059846.1Crambidae, Pyraustinae810.03119−0.192Direct submission
Ostrinia nubilalisNC_054270.1Crambidae, Pyraustinae80.50.03319.4−0.195[70]
Ostrinia scapulalisNC_048887.1Crambidae, Pyraustinae810.0319.1−0.196[51]
Ostrinia zealisNC_048888.1Crambidae, Pyraustinae80.90.03119.1−0.193[51]
Palpita hypohomaliaNC_039632.1Crambidae, Spilomelinae81−0.00118.9−0.196Direct submission
Paracymoriza
distinctalis
NC_023471.1Crambidae, Acentropinae82.2−0.00217.7−0.155[71]
Paracymoriza
prodigalis
NC_020094.1Crambidae, Acentropinae81.50.00218.4−0.183[72]
Paralipsa gularisNC_054356.1Pyralidae, Gallerinae79.5−0.01420.5−0.239Direct submission
Parapediasia
teterrellus
NC_068594.1Crambidae, Crambinae80.5−0.00319.5−0.231Direct submission
Parapoynx crisonalisNC_031151.1Crambinade, Acentropinae820.01718−0.153Direct submission
Perula sp. NC_066226.1Pyralidae, Pyralinae81−0.03719−0.213Direct submission
Plodia interpunctellaNC_027961.1Pyralidae, Phycitinae80.1−0.0519.9−0.233Direct submission
Polythlipta liquidalisNC_073109.1Crambidae, Spilomelinae81−0.00519−0.21Direct submission
Prophantis adustaNC_067853.1Crambidae, Spilomelinae81.50.00118.5−0.196Direct submission
Pseudargyria
interruptella
NC_029751.1Crambidae, Crambinae79.4−0.01120.6−0.216Direct submission
Pseudonoorda
nigropunctalis
NC_056801.1Crambidae, Odontinae81−0.00319−0.201[60]
Pycnarmon
lactiferalis
NC_033540.1Crambidae Spilomelinae81.7−0.00418.3−0.173[73]
Pygospila tyresNC_066087.1Crambidae, Spilomelinae81.3−0.00818.7−0.158Direct submission
Pyrausta despicataNC_046050.1Crambidae, Pyraustinae80.90.00919−0.204Direct submission
Scirpophaga
incertulas
NC_031329.1Crambidae, Schoenobiinae77.20.02922.9−0.32Direct submission
Sinomphisa plagialisNC_061243.1Crambidae, Spilomelinae80.6−0.00819.4−0.216[54]
Sitochroa verticalisNC_062118.1Crambidae, Pyraustinae80.60.00519.5−0.203Direct submission
Syllepte taiwanalisNC_061245.1Crambidae, Pyraustinae81.7−0.00918.3−0.182[54]
Tyspanodes hypsalisNC_025569.1Crambidae, Spilomelinae81.4−0.01718.6−0.175[74]
Tyspanodes striataNC_030510.1Crambidae, Spilomelinae81.3−0.01818.7−0.177Direct submission
The EFSB mitogenome is highly compact, containing mostly genes and the control region with few and small intergenic spacers. The longest spacer is the 69 bp region between trnQ and ND2. The trnQ-ND2 spacer is considered to be a feature of lepidopteran mitogenomes [75] and is thought to have arisen from partial duplication of the ND2 gene [76]. Another spacer in the EFSB mitogenome is the one found between trnS2 and nad1. This spacer is conserved in all insects and contains the motif “ATACTAA” believed to be the recognition site for the mitochondrial transcription termination (mTERM) protein [30,75] (Figure 2A). In addition, overlapping sequences were also found for the EFSB mitogenome. In particular, the 7 bp “ATGATAA” overlap found at the interface of ATP6 and ATP8 is another conserved feature of lepidopteran mitogenomes [30,75] (Figure 2B). For the control region, four regions are of interest that are also conserved features of lepidopterans. The “ATTTA” motif followed by a poly-AT stretch of length 12, and the 13 bp poly-T stretch has also been documented in other lepidopterans [30,75,77]. The “ATAG” box near the 5’end of the s-rna is believed to be the origin of light strand replication for the mitogenome [53]. Interestingly, for other lepidopteran mitogenomes, the “ATAG” box precedes the poly-T stretch [78,79], but the motif is “TTAG” in EFSB (Figure 2C).

3.2. Protein-Coding Genes

The EFSB mitogenome contains 13 protein-coding genes with a total concatenated length of 11,191 bp, comprising 73.41% of the total mitogenome length. All the PCGs are initiated with ATK codons except for nad2 and cox1, which are initiated by TTG codons (Table 2). TTG codons are established initiation codons in invertebrate mitochondria [80], with the TTG codons having relatively high frequencies in the nad2 and cox1 genes of invertebrates [81]. With the PCGs having ATK start codons, ATG (7/13) is slightly preferred over ATT (4/13) codons (Table 2). On the other hand, some of the PCGs of the EFSB mitogenome have incomplete stop codons. The genes COX1, ND5, and CYTB have non-canonical T stop codons, while ATP6 and ND4 have non-canonical TA (Table 2). Interestingly, both MITOS and MitoZ have annotated canonical TAA stop codons for all the PCGs in the EFSB mitogenomes. However, significant overlap of the annotated stop codons with downstream genes has been observed for some PCGs. Thus, manual checking of the annotated PCG gene boundaries was performed to determine if the gene boundaries were consistent with downstream genes. The presence of incomplete stop codons has been documented for invertebrate mitogenomes and are processed into full TAA stop codons post transcriptionally [82,83,84]. This observation is supported by the analysis of Donath et al. (2019), where the authors show that invertebrate COX1, ND5, CYTB, ATP6, and ND4 have significant frequencies of incomplete stop codons [81]. After manual checking, the gene boundaries of the aforementioned genes were adjusted to end in incomplete stop codons. Total amino acids encoded shows that the EFSB mitochondrion codes for a higher proportion of hydrophobic amino acids, with leucine (L), isoleucine (I), and phenylalanine (F) representing three of the top four amino acids used (Figure 3B). This is expected since most of the proteins encoded by mitogenome are components of membrane-bound, hydrophobic complexes important for energy generation [85]. The relative synonymous codon usage is shown in Figure 3A. RSCU analysis indicates that AT rich codons are preferentially used compared to GC rich codons. For example, the EFSB mitochondrion prefers to utilize the codons UUA, AAU, CAA, and AUU to code for the amino acids leucine (Leu), asparagine (Asn), glutamine (Gln), and isoleucine (Ile), respectively, compared to their counterpart codons. The high AT content of the EFSB mitogenome (Table 1) may be correlated to the bias of using AT-rich codons. Overall, the patterns in codon usage are similar to those observed in other members of Pyraloidea [51,86].

3.3. Ribosomal and Transfer RNA Genes

The tRNA genes have a total concatenated length of 1479 bp, contributing 9.7% of the total EFSB mitogenome while the rRNA genes have a total length of 2167 bp, contributing 14.21% of the mitogenome. There are 22 tRNA genes, with serine and leucine having two isoacceptors, while the other 18 amino acids have 1. All the tRNAs have been predicted and verified by MITOS [33] and tRNAScanSE [34] except for trnS2, which yielded no homologous search using tRNAScanSE. This may be due to the four unique unpaired uracils in the anticodon loop (Figure 4). However, such a case of four unpaired uracils was also demonstrated in the trnS2 of the hemp borer, Grapholita delineana [87], the trnS2 of Plodia interpunctella [86], and the trnS2 of five skippers [53]. The length of the tRNAs ranges from 64 bp (trnA and trnR) to 72 bp (trnL1) (Table 2). The secondary structure of tRNAs is generally conserved, forming the typical three-loop, one-stem cloverleaf structures consisting of an amino acid acceptor stem, a dihydrouridine loop (D loop), a pseudouridine loop (ΨU loop), and an anticodon loop. tRNAs also contain a variable region between the anticodon loop and the ΨU loop. All EFSB mitochondrial tRNAs exhibited this typical cloverleaf, except for trnS1, where the D loop fails to form a canonical stem–loop structure (Figure 4). The noncanonical structure of the D-armless trnS1 is a common and conserved feature of metazoan mitogenomes [88]. All the tRNAs have a structurally conserved 7 bp amino acid acceptor stem, except for trnA which contains a 6 bp stem. The tRNAs trnA and trnL2 also contain a unique unpaired uracil in the acceptor stem (Figure 4). Similarly, all the tRNAs contain a conserved 5 bp stem, and a 7 bp loop in the anticodon loop, except for trnL2, which contains a 4 bp stem, and a 9 bp loop; trnK and trnS1 have a 4 bp stem with a 7 bp loop; and trnS2 contains a unique anticodon loop with four unpaired uracils (Figure 4). Noncanonical base pairings also occur in the EFSB mitochondrial tRNAs, with G-U pairs comprising all such pairings (Figure 4). Despite the conservation of both the acceptor stem and anticodon loop, the D and ΨU loops show diversity in the number of stem pairings and the number of nucleotides contained in the hairpin loops. Variability in the D and ΨU may indicate slight deviations from the canonical D/T loop interactions important for aminoacylation [89], but are still functional. The EFSB mitogenome contains two ribosomal RNA genes, l-rna and s-rna. The l-rna and s-rna genes are 1358 bp- and 809 bp-long, respectively. The length of the EFSB mitogenome rRNA genes is similar to that of other pyraloid moths [30].

3.4. Phylogenetic Relationships

A whole mitogenome phylogenetic analysis was used to determine the placement of EFSB within Pyraloidea. The organisms, taxonomy, and reference sequence ID of the mitochondrial genomes used for the analysis are shown in Table 1. For maximum likelihood (ML) and Bayesian inference (BI) analyses, ModelFinder was first used to determine the best-of-fit models for the dataset [47]. For the ML analysis, the initial 37 partitions were reduced to 9 partitions, as summarized in Table S1. ModelFinder found that the general time reversible model (GTR) is the most suitable model for all partitions except partition 3, where the TPM2 model is the most suitable. The GTR model is the model where all six substitution rates and all four base frequencies are unequal. The TPM2 model, on the other hand, assumes that the base frequencies are equal, and the substitution rates have the following relationship: AC=AT, AG=CT, and CG=GT [45]. A similar case was observed when applying ModelFinder for the BI analysis. The initial 37 partitions were also reduced to 9, and ModelFinder found that the GTR model is the most suitable for all partitions for BI. In addition, the corresponding parameters for the best model for each partition were used as the priors for the BI run. All the parameters used as priors are summarized in Table S2. Upon assessing the individual MCMC runs, all three runs converged, as evidenced by the similarly shaped violin plots of the total branch length parameter (Figure S3A). The trace plot for the total branch length across all generations also showed a fairly constant value with small fluctuations around the mean. Lastly, the diagnostic parameters average effective sample size (ESS) and the potential scale reduction factor (PSRF) for the total branch length are 6742.56 and 1.000, respectively (Table S3), indicating that the three MCMC runs converged.
The side-by-side comparison of the obtained phylogenetic trees is shown in Figure 5, while the consensus ML and BI trees are shown in Figure S4. The ML and BI trees are mostly concordant with each other. The general subfamily relationships within the Pyraloidea agree for both trees, with Pyraloidea having the basal split to Crambidae and Pyralidae. The general relationships across Pyralidae also agree with both trees, having recovered the (Gallerinae, (Phycitinae, (Pyralinae, Epipaschiinae))) relationships. Similarly, the ((Spilomelinae, Pyraustinae), ((Crambinae, (Schoenobiinae, Acentropinae)), (Odontinae, Glaphyriinae))) relationships in Crambidae were recovered for both trees (Figure 5). The only disagreements between the ML and BI trees were the order of emergence between Polythlipta liquidalis, Sinomphisa plagialis, and Conogethes punctiferalis. For the ML tree, the clade containing S. plagialis and C. punctiferalis is basal to P. liquidalis, but the BI tree shows that C. punctiferalis is basal to the clade containing P. liquidalis and S. plagialis (Figure 5). Nevertheless, on all other tips, the two trees agree.
The root node for Pyraloidea is highly supported (100% bootstrap support (BS); 1 posterior probability (PP)), indicating the monophyly and shared ancestry for this superfamily. The expected basal split and monophyly of both Crambidae and Pyralidae were recovered, as shown in Figure 5 and Figure S4. The split is highly supported (95.1% BS; 0.999 PP), concordant with the earlier works of Regier et al. (2012) [16] and Léger et al. (2020) [17]. The relationships between the four subfamilies of Pyralidae show high support values. A split in the clade consisting of Gallerinae and the clade consisting of Phycitinae, Pyralinae, and Epipaschiinae is moderately supported in the ML tree (77.8% BS) but well supported in the BI tree (1 PP). Phycitinae was recovered to be a sister clade to the clade containing Pyralinae and Epipaschiinae, and this is also highly supported (100% BS; 1 PP) and consistent with earlier works [17,54,62]. Interestingly, the analysis failed to resolve the subfamily relationships between Pyralinae and Epipaschiinae, exemplified by Orthaga olivacea, which is classified under Epipaschiinae, showing closer relationships to Aglossa dimidiata than Lista haraldusalis, which are Pyralinae and Epipaschiinae, respectively. Similarly, Perula sp., which is Pyralinae, clusters closer to Lista haraldusalis than to other Pyralinae. The relationships between the subfamilies of Crambidae also show relatively good bootstrap support values. Crambidae was recovered to be monophyletic and showed the split between the clade consisting of Pyraustinae and Spilomelinae (PS clade) and the clade containing all other subfamilies (non-PS clade). This expected split is highly supported (95.1% BS; 0.999 PP). The non-PS clade shows a split between the Odontinae-Glaphyriinae (OG) clade and the Crambinae, Acentropinae, Midilinae, Musotiminae, Schoenobiinae, and Scopariinae (CAMMSS) clade, and this split is well-supported (92.8% BS; 0.999 BI). The sister relationships between CAMMSS and OG, and Odontinae and Glaphyriinae, are consistent with the findings of Qi et al. (2021) [60]. In the CAMMSS, Schoenobiinae was recovered to be a sister to Acentropinae. This result is consistent with earlier mitogenomic works [54,90], but different from others [60,62], which places Schoenobiinae as a sister to Crambinae. The subfamilies represented in the non-PS clade were recovered to be monophyletic, and this is well supported for in both trees (Glaphyriinae: 100% BS; 1 PP, Odontinae: 100% BS; 1 PP, Crambinae: 100% BS, 1 PP, Acentropinae: 99.9% BS; 1 PP) (Figure 5).
The general subfamily relationships in Pyraloidea obtained in Figure 5 and Figure S4 are consistent with earlier works [16,17]. Interestingly, while Pyraustinae was recovered to be monophyletic (100% BS; 1 PP), Spilomelinae was recovered to be paraphyletic. The EFSB is traditionally classified under Spilomelinae, but both ML and BI trees show that EFSB emerged before the split of Pyraustinae and Spilomelinae (Figure 5). The placement of EFSB within the PS clade is highly supported (100% BS; 1 PP). Spilomelinae without EFSB is also highly supported (100% BS; 1 PP), indicating that EFSB may have emerged before the PS split. Such a case was also documented for the pyralid moth Orybina regalis. Traditionally placed within Pyralinae, Orybina regalis was recovered to be basal to Gallerinae in the analysis by Liu et al. (2021) [54]. The analysis conducted here also recovered O. regalis to be basal to Gallerinae in both ML and BI trees (Figure 5 and Figure S4). Additionally, while the subfamily relationships recovered in this study are consistent with earlier works, the lower-level relationships vary. There are some species which are robustly placed in other subfamilies. For example, Omphisa fuscidentalis is classified under Spilomelinae but is robustly placed in Pyraustinae in both ML and BI trees (Figure 5 and Figure S4). A similar case was observed for Syllepte taiwanalis and both Cnaphalocrocis species. Both are classified under Pyraustinae but were robustly placed in Spilomelinae. It is interesting to note, however, that such placement for the said taxa was also observed in earlier works [54,90]. The earlier work by Tang and Du (2023) shows consistent subfamily relationships with those obtained in Figure 5 but differ greatly in the lower-level relationships within Spilomelinae [90]. In Figure 5, the tribe Spilomelini is a sister to the clade containing Pycnarmon lactiferalis and Syllepte taiwanalis, and their collective clade is a sister to the tribe Margaroniini, which in Figure 5 is the clade containing the genus Glyphodes and Maruca. This relationship is concordant with Tang and Du (2023) [90]. However, the relationships within Margaroniini are different. For example, Palpita hypohomalia has close relationships with the genus Glyphodes in Figure 5, but the same species has a closer relationship with Cydalima perspectalis in Tang and Du (2023) [90]. Additionally, while the sister relationships between the tribes Trichaeini (Prophantis adusta) and Nomophilini (Nomophila noctuella) can be observed, Nagiella inferior clusters close to these two, which is different from the relationship observed in Tang and Du (2023) [90].
Based on the above analyses, the phylogenetic trees obtained in this work are largely consistent with earlier works. Both ML and BI trees robustly place Leucinodes orbonalis in the PS clade but are basal to the split of Pyraustinae and Spilomelinae. However, there is still the possibility of long-branch attraction. The incorporation of more mitochondrial sequences from Pyraloidea and the incorporation of nuclear genes may be able to resolve the lower-level relationships with Spilomelinae and produce a more conclusive phylogenetic analysis for this economically important subfamily.

4. Summary and Conclusions

The whole mitochondrial genome of the eggplant fruit and shoot borer was assembled and annotated in this study, producing a complete, circular 15,244 bp-long mitogenome that contains the 37 genes expected for a pyraloid mitochondrion. The EFSB mitogenome contains conserved lepidopteran mitogenome features such as the control region motifs, and the trnM-trnI-trnQ-nad2 synteny. Partitioned phylogenetic analysis using the EFSB sequence assembled in this study recovered the expected basal split of the families Crambidae and Pyralidae and the split between the subfamilies Spilomelinae and Pyraustinae. However, Spilomelinae was recovered to be paraphyletic, as indicated by the robust placing of EFSB before the split of Spilomelinae and Pyraustinae in both ML and BI trees. The addition of more mitochondrial sequences from Pyraloidea and the incorporation of nuclear genes should yield a more resolved phylogeny that details the genus-level relationships within the economically important Spilomelinae. Overall, the mitogenome produced here contributes to the scarce genetic resources for the EFSB and will be a great help in future phylogenetic reconstruction of Spilomelinae and in future population studies and studies that will impact EFSB management.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/2075-4450/15/4/220/s1, Figure S1. Sequence quality of the EFSB mitogenome reads; Figure S2. Cumulative coverage distribution for the EFSB mitogenome; Figure S3. Estimation of the total branch length parameter across 3 MCMC runs; Figure S4. Partitioned phylogenetic analysis of 72 members of the superfamily Pyraloidea and 3 dipterans using maximum likelihood and Bayesian inference; Table S1. Final gene partitions for the Maximum Likelihood Phylogenetic analysis; Table S2. Final gene partitions for the Bayesian Inference Phylogenetic analysis; Table S3. 95% Credibility interval of the sampled parameters during the MCMC run.

Author Contributions

Conceptualization, J.B.D. and M.A.M.B.; methodology, J.B.D. and M.A.M.B.; formal analysis, J.B.D.; investigation, J.B.D.; resources, M.A.M.B.; writing—original draft preparation, J.B.D. and M.A.M.B.; writing—review and editing, J.B.D. and M.A.M.B.; visualization, J.B.D.; supervision, M.A.M.B.; funding acquisition, M.A.M.B. All authors have read and agreed to the published version of the manuscript.

Funding

The funding for this work was provided by the Department of Science and Technology Philippine Council for Agriculture, Aquatic and Natural Resources Research and Development (DOST-PCAARRD) to University of the Philippines Los Baños and University of the Philippines Diliman Additional funding for the bioinformatics work was also obtained from the Accelerated Science and Technology Human Resource Development Program (ASTHRDP) of the Department of Science and Technology Science Education Institute (DOST-SEI). Funding was also provided by the National Institute of Molecular Biology and Biotechnology, University of the Philippines Diliman (NIMBB-UPD) in-house funds to sustain further analysis.

Data Availability Statement

The sequence of the final mitogenome assembly has been deposited into GenBank under accession number PP493058. The final mitogenome assembly, intermediate files, and scripts used for the analysis of the mitogenome are deposited and openly available in Zenodo at 10.5281/zenodo.10653202.

Acknowledgments

The authors would also like to thank the DNA Sequencing Core Facility (DSCF) and Core Facility for Bioinformatics (CFB) of the Philippine Genome Center (PGC) for their invaluable help in sequencing, and preliminary bioinformatics analysis for this work, respectively. Lastly, the authors would also like to thank the Computing and Archiving Research Environment (COARE) of the Department of Science and Technology Advanced Science and Technology Institute (DOST-ASTI) for the HPC services that facilitated the phylogenetic analyses covered in this work.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of the data; in the writing of the manuscript; or in the decision to publish the results.

Correction Statement

This article has been republished with a minor correction to the readability of table. This change does not affect the scientific content of the article.

References

  1. Martin Embley, T. Multiple Secondary Origins of the Anaerobic Lifestyle in Eukaryotes. Phil. Trans. R. Soc. B 2006, 361, 1055–1067. [Google Scholar] [CrossRef]
  2. Muñoz-Gómez, S.A.; Slamovits, C.H.; Dacks, J.B.; Baier, K.A.; Spencer, K.D.; Wideman, J.G. Ancient Homology of the Mitochondrial Contact Site and Cristae Organizing System Points to an Endosymbiotic Origin of Mitochondrial Cristae. Curr. Biol. 2015, 25, 1489–1495. [Google Scholar] [CrossRef] [PubMed]
  3. Emelyanov, V.V. Mitochondrial Connection to the Origin of the Eukaryotic Cell. Eur. J. Biochem. 2003, 270, 1599–1618. [Google Scholar] [CrossRef] [PubMed]
  4. Gray, M.W. Mitochondrial Evolution. Cold Spring Harb. Perspect. Biol. 2012, 4, a011403. [Google Scholar] [CrossRef]
  5. Karnkowska, A.; Vacek, V.; Zubáčová, Z.; Treitli, S.C.; Petrželková, R.; Eme, L.; Novák, L.; Žárský, V.; Barlow, L.D.; Herman, E.K.; et al. A Eukaryote without a Mitochondrial Organelle. Curr. Biol. 2016, 26, 1274–1284. [Google Scholar] [CrossRef]
  6. Chan, D.C. Mitochondria: Dynamic Organelles in Disease, Aging, and Development. Cell 2006, 125, 1241–1252. [Google Scholar] [CrossRef] [PubMed]
  7. Hebert, P.D.N.; Ratnasingham, S.; De Waard, J.R. Barcoding Animal Life: Cytochrome c Oxidase Subunit 1 Divergences among Closely Related Species. Proc. R. Soc. Lond. B 2003, 270, S96–S99. [Google Scholar] [CrossRef]
  8. Ma, C.; Yang, P.; Jiang, F.; Chapuis, M.; Shali, Y.; Sword, G.A.; Kang, L. Mitochondrial Genomes Reveal the Global Phylogeography and Dispersal Routes of the Migratory Locust. Mol. Ecol. 2012, 21, 4344–4358. [Google Scholar] [CrossRef]
  9. Nelson, L.A.; Lambkin, C.L.; Batterham, P.; Wallman, J.F.; Dowton, M.; Whiting, M.F.; Yeates, D.K.; Cameron, S.L. Beyond Barcoding: A Mitochondrial Genomics Approach to Molecular Phylogenetics and Diagnostics of Blowflies (Diptera: Calliphoridae). Gene 2012, 511, 131–142. [Google Scholar] [CrossRef]
  10. Castro, L.R.; Austin, A.D.; Dowton, M. Contrasting Rates of Mitochondrial Molecular Evolution in Parasitic Diptera and Hymenoptera. Mol. Biol. Evol. 2002, 19, 1100–1113. [Google Scholar] [CrossRef]
  11. Shao, R. Rates of Gene Rearrangement and Nucleotide Substitution Are Correlated in the Mitochondrial Genomes of Insects. Mol. Biol. Evol. 2003, 20, 1612–1619. [Google Scholar] [CrossRef]
  12. Som, A. Causes, Consequences and Solutions of Phylogenetic Incongruence. Brief. Bioinform. 2015, 16, 536–548. [Google Scholar] [CrossRef]
  13. Philippe, H.; Brinkmann, H.; Lavrov, D.V.; Littlewood, D.T.J.; Manuel, M.; Wörheide, G.; Baurain, D. Resolving Difficult Phylogenetic Questions: Why More Sequences Are Not Enough. PLoS Biol. 2011, 9, e1000602. [Google Scholar] [CrossRef] [PubMed]
  14. Navasero, M.V.; Candano, R.N.; Hautea, D.M.; Hautea, R.A.; Shotkoski, F.A.; Shelton, A.M. Assessing Potential Impact of Bt Eggplants on Non-Target Arthropods in the Philippines. PLoS ONE 2016, 11, e0165190. [Google Scholar] [CrossRef] [PubMed]
  15. van Nieukerken, E.J.; Kaila, L.; Kitching, I.; Kristensen, N.P.; Lees, D.C.; Minet, J.; Zwick, A. Order Lepidoptera Linnaeus, 1758. In Animal Biodiversity: An Outline of Higher-Level Classification and Survey of Taxonomic Richness; Zootaxa: Auckland, New Zealand, 2011; Volume 3148. [Google Scholar]
  16. Regier, J.C.; Mitter, C.; Solis, M.A.; Hayden, J.E.; Landry, B.; Nuss, M.; Simonsen, T.J.; Yen, S.; Zwick, A.; Cummings, M.P. A Molecular Phylogeny for the Pyraloid Moths (Lepidoptera: Pyraloidea) and Its Implications for Higher-level Classification. Syst. Entomol. 2012, 37, 635–656. [Google Scholar] [CrossRef]
  17. Léger, T.; Mally, R.; Neinhuis, C.; Nuss, M. Refining the Phylogeny of Crambidae with Complete Sampling of Subfamilies (Lepidoptera, Pyraloidea). Zool. Scr. 2021, 50, 84–99. [Google Scholar] [CrossRef]
  18. Sagarbarria, M.G.S.; Marasigan, J.C.M.; Taylo, L.D.; Hautea, D.M. Identity and Genetic Structure of Eggplant Fruit and Shoot Borer, Leucinodes Orbonalis Guenée (Lepidoptera:Crambidae) Populations in the Philippines Inferred from Morphological Traits and COI Sequence Data. J. Asia-Pac. Entomol. 2018, 21, 1009–1019. [Google Scholar] [CrossRef]
  19. Dowling, D.K.; Tompkins, D.M.; Gemmell, N.J. The T Rojan F Emale T Echnique for Pest Control: A Candidate Mitochondrial Mutation Confers Low Male Fertility across Diverse Nuclear Backgrounds in Drosophila melanogaster. Evol. Appl. 2015, 8, 871–880. [Google Scholar] [CrossRef]
  20. Mally, R.; Hayden, J.E.; Neinhuis, C.; Jordal, B.H.; Nuss, M. The Phylogenetic Systematics of Spilomelinae and Pyraustinae (Lepidoptera: Pyraloidea: Crambidae) Inferred from DNA and Morphology. Arthropod Syst. Phylogeny 2019, 77, 141–204. [Google Scholar] [CrossRef]
  21. Yi, M.-R.; Hsu, K.-C.; Gu, S.; He, X.-B.; Luo, Z.-S.; Lin, H.-D.; Yan, Y.-R. Complete Mitogenomes of Four Trichiurus Species: A Taxonomic Review of the T. Lepturus Species Complex. ZK 2022, 1084, 1–26. [Google Scholar] [CrossRef]
  22. Gong, J.; Chen, B.; Li, B.; Zhou, Z.; Shi, Y.; Ke, Q.; Zhang, D.; Xu, P. Genetic Analysis of Whole Mitochondrial Genome of Lateolabrax Maculatus (Perciformes: Moronidae) Indicates the Presence of Two Populations along the Chinese Coast. Zoologia 2020, 37, e49046. [Google Scholar] [CrossRef]
  23. Alonso, D.P.; Alvarez, M.V.N.; Amorim, J.A.; De Sá, I.L.R.; De Carvalho, D.P.; Ribeiro, K.A.N.; Ribolla, P.E.M.; Sallum, M.A.M. Mansonia Spp. Population Genetics Based on Mitochondrion Whole-Genome Sequencing alongside the Madeira River near Porto Velho, Rondonia, Brazil. Infect. Genet. Evol. 2022, 103, 105341. [Google Scholar] [CrossRef] [PubMed]
  24. Kingan, S.; Heaton, H.; Cudini, J.; Lambert, C.; Baybayan, P.; Galvin, B.; Durbin, R.; Korlach, J.; Lawniczak, M. A High-Quality De Novo Genome Assembly from a Single Mosquito Using PacBio Sequencing. Genes 2019, 10, 62. [Google Scholar] [CrossRef] [PubMed]
  25. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. Fastp: An Ultra-Fast All-in-One FASTQ Preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef]
  26. Andrews, S. FastQC A Quality Control Tool for High Throughput Sequence Data. Available online: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 8 May 2022).
  27. Langmead, B.; Salzberg, S.L. Fast Gapped-Read Alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef]
  28. Park, J.S.; Kim, M.J.; Ahn, S.-J.; Kim, I. Complete Mitochondrial Genome of the Grass Moth Glyphodes quadrimaculalis (Lepidoptera: Crambidae). Mitochondrial DNA 2015, 26, 247–249. [Google Scholar] [CrossRef] [PubMed]
  29. Zou, Y.; Ma, W.; Zhang, L.; He, S.; Zhang, X.; Tao, Z. The Complete Mitochondrial Genome of the Bean Pod Borer, Maruca Testulalis (Lepidoptera: Crambidae: Spilomelinae). Mitochondrial DNA 2016, 27, 740–741. [Google Scholar] [CrossRef]
  30. Yang, M.; Song, L.; Mao, J.; Shi, Y.; Wu, C.; Zhang, Y.; Huang, L.; Peng, W.; Liu, X. Complete Mitochondrial Genome of the Soybean Leaffolder, Omiodes Indicata (Lepidoptera: Pyraloidea: Crambidae), and Phylogenetic Analysis for Pyraloidea. Int. J. Biol. Macromol. 2018, 115, 53–60. [Google Scholar] [CrossRef]
  31. Nurk, S.; Meleshko, D.; Korobeynikov, A.; Pevzner, P.A. metaSPAdes: A New Versatile Metagenomic Assembler. Genome Res. 2017, 27, 824–834. [Google Scholar] [CrossRef]
  32. Meng, G.; Li, Y.; Yang, C.; Liu, S. MitoZ: A Toolkit for Animal Mitochondrial Genome Assembly, Annotation and Visualization. Nucleic Acids Res. 2019, 47, e63. [Google Scholar] [CrossRef]
  33. Bernt, M.; Donath, A.; Jühling, F.; Externbrink, F.; Florentz, C.; Fritzsch, G.; Pütz, J.; Middendorf, M.; Stadler, P.F. MITOS: Improved de Novo Metazoan Mitochondrial Genome Annotation. Mol. Phylogenetics Evol. 2013, 69, 313–319. [Google Scholar] [CrossRef]
  34. Chan, P.P.; Lin, B.Y.; Mak, A.J.; Lowe, T.M. tRNAscan-SE 2.0: Improved Detection and Functional Classification of Transfer RNA Genes. Nucleic Acids Res. 2021, 49, 9077–9096. [Google Scholar] [CrossRef]
  35. Reuter, J.S.; Mathews, D.H. RNAstructure: Software for RNA Secondary Structure Prediction and Analysis. BMC Bioinform. 2010, 11, 129. [Google Scholar] [CrossRef]
  36. Li, H. Aligning Sequence Reads, Clone Sequences and Assembly Contigs with BWA-MEM. arXiv 2013, arXiv:1303.3997. [Google Scholar] [CrossRef]
  37. Danecek, P.; Bonfield, J.K.; Liddle, J.; Marshall, J.; Ohan, V.; Pollard, M.O.; Whitwham, A.; Keane, T.; McCarthy, S.A.; Davies, R.M.; et al. Twelve Years of SAMtools and BCFtools. GigaScience 2021, 10, giab008. [Google Scholar] [CrossRef]
  38. Pedersen, B.S.; Quinlan, A.R. Mosdepth: Quick Coverage Calculation for Genomes and Exomes. Bioinformatics 2018, 34, 867–868. [Google Scholar] [CrossRef] [PubMed]
  39. Zhang, D.; Gao, F.; Jakovlić, I.; Zou, H.; Zhang, J.; Li, W.X.; Wang, G.T. PhyloSuite: An Integrated and Scalable Desktop Platform for Streamlined Molecular Sequence Data Management and Evolutionary Phylogenetics Studies. Mol. Ecol. Resour. 2020, 20, 348–355. [Google Scholar] [CrossRef] [PubMed]
  40. Lu, J.; Salzberg, S.L. SkewIT: The Skew Index Test for Large-Scale GC Skew Analysis of Bacterial Genomes. PLoS Comput. Biol. 2020, 16, e1008439. [Google Scholar] [CrossRef]
  41. Krzywinski, M.; Schein, J.; Birol, İ.; Connors, J.; Gascoyne, R.; Horsman, D.; Jones, S.J.; Marra, M.A. Circos: An Information Aesthetic for Comparative Genomics. Genome Res. 2009, 19, 1639–1645. [Google Scholar] [CrossRef]
  42. Katoh, K.; Standley, D.M. MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef]
  43. 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]
  44. Chernomor, O.; Von Haeseler, A.; Minh, B.Q. Terrace Aware Data Structure for Phylogenomic Inference from Supermatrices. Syst. Biol. 2016, 65, 997–1008. [Google Scholar] [CrossRef] [PubMed]
  45. 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] [PubMed]
  46. Ronquist, F.; Teslenko, M.; Van Der Mark, P.; Ayres, D.L.; Darling, A.; Höhna, S.; Larget, B.; Liu, L.; Suchard, M.A.; Huelsenbeck, J.P. MrBayes 3.2: Efficient Bayesian Phylogenetic Inference and Model Choice Across a Large Model Space. Syst. Biol. 2012, 61, 539–542. [Google Scholar] [CrossRef] [PubMed]
  47. Kalyaanamoorthy, S.; Minh, B.Q.; Wong, T.K.F.; Von Haeseler, A.; Jermiin, L.S. ModelFinder: Fast Model Selection for Accurate Phylogenetic Estimates. Nat. Methods 2017, 14, 587–589. [Google Scholar] [CrossRef] [PubMed]
  48. Ayres, D.L.; Darling, A.; Zwickl, D.J.; Beerli, P.; Holder, M.T.; Lewis, P.O.; Huelsenbeck, J.P.; Ronquist, F.; Swofford, D.L.; Cummings, M.P.; et al. BEAGLE: An Application Programming Interface and High-Performance Computing Library for Statistical Phylogenetics. Syst. Biol. 2012, 61, 170–173. [Google Scholar] [CrossRef] [PubMed]
  49. Yu, G.; Smith, D.K.; Zhu, H.; Guan, Y.; Lam, T.T. GGTREE: An R Package for Visualization and Annotation of Phylogenetic Trees with Their Covariates and Other Associated Data. Methods Ecol. Evol. 2017, 8, 28–36. [Google Scholar] [CrossRef]
  50. Park, J.S.; Kim, M.J.; Kim, S.-S.; Kim, I. Complete Mitochondrial Genome of an Aquatic Moth, Elophila interruptalis (Lepidoptera: Crambidae). Mitochondrial DNA 2014, 25, 275–277. [Google Scholar] [CrossRef]
  51. Zhou, N.; Dong, Y.; Qiao, P.; Yang, Z. Complete Mitogenomic Structure and Phylogenetic Implications of the Genus Ostrinia (Lepidoptera: Crambidae). Insects 2020, 11, 232. [Google Scholar] [CrossRef]
  52. Liu, D.; Basso, A.; Babbucci, M.; Patarnello, T.; Negrisolo, E. Macrostructural Evolution of the Mitogenome of Butterflies (Lepidoptera, Papilionoidea). Insects 2022, 13, 358. [Google Scholar] [CrossRef]
  53. Kim, M.J.; Wang, A.R.; Park, J.S.; Kim, I. Complete Mitochondrial Genomes of Five Skippers (Lepidoptera: Hesperiidae) and Phylogenetic Reconstruction of Lepidoptera. Gene 2014, 549, 97–112. [Google Scholar] [CrossRef]
  54. Liu, X.; Qi, M.; Xu, H.; Wu, Z.; Hu, L.; Yang, M.; Li, H. Nine Mitochondrial Genomes of the Pyraloidea and Their Phylogenetic Implications (Lepidoptera). Insects 2021, 12, 1039. [Google Scholar] [CrossRef]
  55. Roh, S.J.; Jeon, J.H.; Kim, D.-S.; Byun, B.-K. The Complete Mitochondrial Genome of Unique Snout Moth, Cathayia Obliquella (Pyralidae: Galleriinae) with Its Phylogenetic Implications and a Redescription of External Morphology. J. Asia-Pac. Biodivers. 2020, 13, 613–624. [Google Scholar] [CrossRef]
  56. Chai, H.-N.; Du, Y.-Z.; Zhai, B.-P. Characterization of the Complete Mitochondrial Genomes of Cnaphalocrocis medinalis and Chilo suppressalis (Lepidoptera: Pyralidae). Int. J. Biol. Sci. 2012, 8, 561–579. [Google Scholar] [CrossRef]
  57. Wu, Q.-L.; Gong, Y.-J.; Shi, B.; Gu, Y.; Wei, S.-J. The Complete Mitochondrial Genome of the Yellow Peach Moth Dichocrocis punctiferalis (Lepidoptera: Pyralidae). Mitochondrial DNA 2013, 24, 105–107. [Google Scholar] [CrossRef] [PubMed]
  58. Behura, S.K.; Lobo, N.F.; Haas, B.; deBruyn, B.; Lovin, D.D.; Shumway, M.F.; Puiu, D.; Romero-Severson, J.; Nene, V.; Severson, D.W. Complete Sequences of Mitochondria Genomes of Aedes Aegypti and Culex Quinquefasciatus and Comparative Analysis of Mitochondrial DNA Fragments Inserted in the Nuclear Genomes. Insect Biochem. Mol. Biol. 2011, 41, 770–777. [Google Scholar] [CrossRef] [PubMed]
  59. Que, S.; Yu, A.; Tu, Y.; Xiong, C.; Liu, X. Complete Mitochondrial Genome and Phylogenetic Analysis of Diaphania perspectalis. Mitochondrial DNA Part. B 2019, 4, 933–934. [Google Scholar] [CrossRef]
  60. Qi, M.; Zhao, H.; Yu, F.; Zhang, A.; Li, H. The First Mitogenomes of the Subfamily Odontiinae (Lepidoptera, Crambidae) and Phylogenetic Analysis of Pyraloidea. Insects 2021, 12, 486. [Google Scholar] [CrossRef] [PubMed]
  61. Li, W.; Zhang, X.; Fan, Z.; Yue, B.; Huang, F.; King, E.; Ran, J. Structural Characteristics and Phylogenetic Analysis of the Mitochondrial Genome of the Sugarcane Borer, Diatraea saccharalis (Lepidoptera: Crambidae). DNA Cell Biol. 2011, 30, 3–8. [Google Scholar] [CrossRef] [PubMed]
  62. Zhu, W.; Yan, J.; Song, J.; You, P. The First Mitochondrial Genomes for Pyralinae (Pyralidae) and Glaphyriinae (Crambidae), with Phylogenetic Implications of Pyraloidea. PLoS ONE 2018, 13, e0194672. [Google Scholar] [CrossRef] [PubMed]
  63. Liu, Q.; Jiang, X.; Hou, X.; Yang, H.; Chen, W. The Mitochondrial Genome of Ephestia elutella (Insecta: Lepidoptera: Pyralidae). Mitochondrial DNA Part. B 2018, 3, 189–190. [Google Scholar] [CrossRef] [PubMed]
  64. Traut, W.; Vogel, H.; Glöckner, G.; Hartmann, E.; Heckel, D.G. High-Throughput Sequencing of a Single Chromosome: A Moth W Chromosome. Chromosome Res. 2013, 21, 491–505. [Google Scholar] [CrossRef] [PubMed]
  65. Ye, F.; Yu, H.-L.; Li, P.-F.; You, P. The Complete Mitochondrial Genome of Lista haraldusalis (Lepidoptera: Pyralidae). Mitochondrial DNA 2015, 26, 853–854. [Google Scholar] [CrossRef] [PubMed]
  66. Living Prairie Mitogenomics Consortium the Complete Mitochondrial Genome of the Lesser Aspen Webworm Moth Meroptera pravella (Insecta: Lepidoptera: Pyralidae). Mitochondrial DNA Part. B 2017, 2, 344–346. [CrossRef] [PubMed]
  67. Tang, M.; Tan, M.; Meng, G.; Yang, S.; Su, X.; Liu, S.; Song, W.; Li, Y.; Wu, Q.; Zhang, A.; et al. Multiplex Sequencing of Pooled Mitochondrial Genomes—a Crucial Step toward Biodiversity Analysis Using Mito-Metagenomics. Nucleic Acids Res. 2014, 42, e166. [Google Scholar] [CrossRef] [PubMed]
  68. Mao, X.; Li, H.-X.; Yang, M.-F. The Complete Mitochondrial Genome Sequences of Two Insect-Tea Producers in Pyralidae (Lepidoptera) from South China: Pyralis farinalis and Orthopygia glaucinalis. Mitochondrial DNA Part. B 2019, 4, 3850–3851. [Google Scholar] [CrossRef]
  69. Li, C.; Li, L.; Ren, Y.; Lu, Z.; Song, Y.; Liu, L.; Lv, S.; Yu, Y.; Men, X. Characterization of the Complete Mitochondrial Genome of Asia Corn Borer, Ostrinia furnacalis (Lepidoptera: Crambidae). Mitochondrial DNA Part. B 2020, 5, 936–937. [Google Scholar] [CrossRef]
  70. Fisher, K.E.; Bradbury, S.P.; Coates, B.S. Prediction of Mitochondrial Genome-Wide Variation through Sequencing of Mitochondrion-Enriched Extracts. Sci. Rep. 2020, 10, 19123. [Google Scholar] [CrossRef]
  71. Ye, F.; You, P. The Complete Mitochondrial Genome of Paracymoriza distinctalis (Lepidoptera: Crambidae). Mitochondrial DNA 2016, 27, 28–29. [Google Scholar] [CrossRef]
  72. Ye, F.; Shi, Y.; Xing, L.; Yu, H.; You, P. The Complete Mitochondrial Genome of Paracymoriza prodigalis (Leech, 1889) (Lepidoptera), with a Preliminary Phylogenetic Analysis of Pyraloidea. Aquat. Insects 2013, 35, 71–88. [Google Scholar] [CrossRef]
  73. Chen, S.; Li, F.-H.; Lan, X.-E.; You, P. The Complete Mitochondrial Genome of Pycnarmon lactiferalis (Lepidoptera: Crambidae). Mitochondrial DNA Part. B 2016, 1, 638–639. [Google Scholar] [CrossRef]
  74. Wang, J.; Li, P.; You, P. The Complete Mitochondrial Genome of Tyspanodes hypsalis (Lepidoptera: Crambidae). Mitochondrial DNA 2014, 27, 1821–1822. [Google Scholar] [CrossRef]
  75. Wu, Y.-P.; Zhao, J.-L.; Su, T.-J.; Li, J.; Yu, F.; Chesters, D.; Fan, R.-J.; Chen, M.-C.; Wu, C.-S.; Zhu, C.-D. The Complete Mitochondrial Genome of Leucoptera malifoliella Costa (Lepidoptera: Lyonetiidae). DNA Cell Biol. 2012, 31, 1508–1522. [Google Scholar] [CrossRef]
  76. Kim, M.I.; Baek, J.Y.; Kim, M.J.; Jeong, H.C.; Kim, K.-G.; Bae, C.H.; Han, Y.S.; Jin, B.R.; Kim, I. Complete Nucleotide Sequence and Organization of the Mitogenome of the Red-Spotted Apollo Butterfly, Parnassius bremeri (Lepidoptera: Papilionidae) and Comparison with Other Lepidopteran Insects. Mol. Cells 2009, 28, 347–364. [Google Scholar] [CrossRef]
  77. Ramírez-Ríos, V.; Alvarez, J.C.; Villanueva-Mejia, D. Mitochondrial Genomes of Lepidopteran Insects Considered Crop Pests. In Lepidoptera; Perveen, F.K., Ed.; InTech: London, UK, 2017; ISBN 978-953-51-3659-0. [Google Scholar]
  78. Yang, Y.; Wu, Z.; Xu, H.; Zheng, X.; Tian, J.; Lu, Z. Sequence and Analysis of mtDNA Control Region in Cnaphalacrocis medinalis Guenée and Marasmia patnalis Bradley (Lepidoptera: Pyralidae). Mitochondrial DNA Part. A 2016, 27, 3991–3994. [Google Scholar] [CrossRef] [PubMed]
  79. He, S.-L.; Zou, Y.; Zhang, L.-F.; Ma, W.-Q.; Zhang, X.-Y.; Yue, B.-S. The Complete Mitochondrial Genome of the Beet Webworm, Spoladea Recurvalis (Lepidoptera: Crambidae) and Its Phylogenetic Implications. PLoS ONE 2015, 10, e0129355. [Google Scholar] [CrossRef] [PubMed]
  80. Quek, Z.B.R.; Chang, J.J.M.; Ip, Y.C.A.; Chan, Y.K.S.; Huang, D. Mitogenomes Reveal Alternative Initiation Codons and Lineage-Specific Gene Order Conservation in Echinoderms. Mol. Biol. Evol. 2021, 38, 981–985. [Google Scholar] [CrossRef] [PubMed]
  81. Donath, A.; Jühling, 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] [PubMed]
  82. Ojala, D.; Montoya, J.; Attardi, G. tRNA Punctuation Model of RNA Processing in Human Mitochondria. Nature 1981, 290, 470–474. [Google Scholar] [CrossRef] [PubMed]
  83. Stewart, J.B.; Beckenbach, A.T. Characterization of Mature Mitochondrial Transcripts in Drosophila, and the Implications for the tRNA Punctuation Model in Arthropods. Gene 2009, 445, 49–57. [Google Scholar] [CrossRef]
  84. Taanman, J.-W. The Mitochondrial Genome: Structure, Transcription, Translation and Replication. Biochim. Et. Biophys. Acta (BBA) - Bioenerg. 1999, 1410, 103–123. [Google Scholar] [CrossRef]
  85. Johnston, I.G.; Williams, B.P. Evolutionary Inference across Eukaryotes Identifies Specific Pressures Favoring Mitochondrial Gene Retention. Cell Syst. 2016, 2, 101–111. [Google Scholar] [CrossRef] [PubMed]
  86. Liu, Q.-N.; Chai, X.-Y.; Bian, D.-D.; Zhou, C.-L.; Tang, B.-P. The Complete Mitochondrial Genome of Plodia interpunctella (Lepidoptera: Pyralidae) and Comparison with Other Pyraloidea Insects. Genome 2016, 59, 37–49. [Google Scholar] [CrossRef] [PubMed]
  87. Song, L.; Shi, Y.; Zhang, H.; Wang, Z.; Liu, X.; Yang, M. Complete Mitochondrial Genome of the Hemp Borer, Grapholita Delineana (Lepidoptera: Tortricidae): Gene Variability and Phylogeny among Grapholita. J. Asia-Pac. Entomol. 2021, 24, 250–258. [Google Scholar] [CrossRef]
  88. Watanabe, Y.; Suematsu, T.; Ohtsuki, T. Losing the Stem-Loop Structure from Metazoan Mitochondrial tRNAs and Co-Evolution of Interacting Factors. Front. Genet. 2014, 5, 109. [Google Scholar] [CrossRef]
  89. Salinas-Giegé, T.; Giegé, R.; Giegé, P. tRNA Biology in Mitochondria. IJMS 2015, 16, 4518–4559. [Google Scholar] [CrossRef]
  90. Tang, C.; Du, X. Complete Mitochondrial Genomes of Two Moths in the Tribe Trichaeini (Lepidoptera: Crambidae) and Their Phylogenetic Implications. Ecol. Evol. 2023, 13, e10188. [Google Scholar] [CrossRef]
Figure 1. Circos plot detailing the assembled EFSB mitogenome and various statistics. (A) Gene organization of the EFSB mitogenome. Yellow, purple, pink, and grey boxes indicate protein-coding genes, tRNA genes, rRNA genes, and the control region, respectively. Boxes oriented inward are genes transcribed in the majority strand, while the boxes oriented outward are genes transcribed in the minority strand. (B) The GC content histogram was calculated in 20 bp sliding windows. Black and red lines correspond to 50% and 25% GC content, respectively. (C) GC-skew was calculated using a window size of 100 bp and a step size of 20 bp. Black line corresponds to 0 GC-skew value and histogram bins colored yellow and red have negative and positive GC skew values, respectively. (D) Coverage distribution across all positions. The EFSB mitogenome has a mean coverage of 9571×. Positions having coverage values above and below the mean coverage value are colored purple and green, respectively. The highlighted sector in dotted lines indicates the control region of the EFSB mitogenome characterized by low GC content and a switch in the sign of the GC skew values.
Figure 1. Circos plot detailing the assembled EFSB mitogenome and various statistics. (A) Gene organization of the EFSB mitogenome. Yellow, purple, pink, and grey boxes indicate protein-coding genes, tRNA genes, rRNA genes, and the control region, respectively. Boxes oriented inward are genes transcribed in the majority strand, while the boxes oriented outward are genes transcribed in the minority strand. (B) The GC content histogram was calculated in 20 bp sliding windows. Black and red lines correspond to 50% and 25% GC content, respectively. (C) GC-skew was calculated using a window size of 100 bp and a step size of 20 bp. Black line corresponds to 0 GC-skew value and histogram bins colored yellow and red have negative and positive GC skew values, respectively. (D) Coverage distribution across all positions. The EFSB mitogenome has a mean coverage of 9571×. Positions having coverage values above and below the mean coverage value are colored purple and green, respectively. The highlighted sector in dotted lines indicates the control region of the EFSB mitogenome characterized by low GC content and a switch in the sign of the GC skew values.
Insects 15 00220 g001
Figure 2. Intergenic spacers, overlapping regions, and conserved motifs in the EFSB mitogenome. (A) Intergenic spacers found between trnQ-ND2 and trnS2-ND1. Sequences in black underlines are the spacers (B) Gene overlap regions between ATP8-ATP6. Sequences highlighted in grey are the overlapping regions. (C) Control region motifs found in the EFSB mitogenome. Sequences with the red line indicate conserved motifs, and the arrow direction shows the transcription direction of the genes, which are colored for reference.
Figure 2. Intergenic spacers, overlapping regions, and conserved motifs in the EFSB mitogenome. (A) Intergenic spacers found between trnQ-ND2 and trnS2-ND1. Sequences in black underlines are the spacers (B) Gene overlap regions between ATP8-ATP6. Sequences highlighted in grey are the overlapping regions. (C) Control region motifs found in the EFSB mitogenome. Sequences with the red line indicate conserved motifs, and the arrow direction shows the transcription direction of the genes, which are colored for reference.
Insects 15 00220 g002
Figure 3. Relative synonymous codon usage and total amino acids encoded by the EFSB mitogenome. (A) Relative synonymous codon usage (RSCU) is a measure of the deviation of each codon from the assumption that all codons are used equally. The bar graph shows the relative distribution of codon usage across all amino acid codons. The number of codons for each amino acid ranges from two to eight. The RSCU values are color-coded based on the codons below the amino acid labels. (B) Amino acids encoded by the EFSB mitogenome. Each amino acid is colored differently, and the labels in the circles correspond to the one-letter amino acid abbreviations. The number of times a given amino acid was coded is proportional to the area of the circles and the size of the text labels.
Figure 3. Relative synonymous codon usage and total amino acids encoded by the EFSB mitogenome. (A) Relative synonymous codon usage (RSCU) is a measure of the deviation of each codon from the assumption that all codons are used equally. The bar graph shows the relative distribution of codon usage across all amino acid codons. The number of codons for each amino acid ranges from two to eight. The RSCU values are color-coded based on the codons below the amino acid labels. (B) Amino acids encoded by the EFSB mitogenome. Each amino acid is colored differently, and the labels in the circles correspond to the one-letter amino acid abbreviations. The number of times a given amino acid was coded is proportional to the area of the circles and the size of the text labels.
Insects 15 00220 g003
Figure 4. Structures of the EFSB mitogenome tRNAs. All of the EFSB mitogenome tRNAs except for trnS1 follow the canonical conserved 3-loop cloverleaf structure consisting of the acceptor stem (dark green), dihydrouridine loop (D loop, orange), the pseudouridine arm (ΨU loop, purple), and the anticodon loop (light green). All tRNAs also contain a variable region (pink). In the case of trnS1, the D loop fails to form a stem–loop structure. The anticodon for each tRNA is indicated in yellow, while all non-canonical base pairs are indicated by a red bond. All tRNA structures of the EFSB mitogenome were predicted using MITOS, validated using tRNAScan-SE, and drawn using StructureEditor.
Figure 4. Structures of the EFSB mitogenome tRNAs. All of the EFSB mitogenome tRNAs except for trnS1 follow the canonical conserved 3-loop cloverleaf structure consisting of the acceptor stem (dark green), dihydrouridine loop (D loop, orange), the pseudouridine arm (ΨU loop, purple), and the anticodon loop (light green). All tRNAs also contain a variable region (pink). In the case of trnS1, the D loop fails to form a stem–loop structure. The anticodon for each tRNA is indicated in yellow, while all non-canonical base pairs are indicated by a red bond. All tRNA structures of the EFSB mitogenome were predicted using MITOS, validated using tRNAScan-SE, and drawn using StructureEditor.
Insects 15 00220 g004
Figure 5. Topological concordance between the partitioned maximum likelihood and Bayesian inference phylogenetic trees. The maximum likelihood (ML) and Bayesian inference (BI) phylogenetic trees are shown on the left- and right-hand side of the labels, respectively. The labeling scheme follows the tip order from the ML tree as indicated by the dotted lines. The subfamilies are coded according to color. The grey lines connecting the tips of both trees indicate tips that agree for both trees, while the red lines indicate the tips with a different order and topology between the trees. SH_aLRT and posterior probability support values are shown for the ML and BI trees, respectively. The subfamilies included in the analysis are Spilomelinae, Pyraustinae, Acentropinae, Schoenobiinae, Crambinae, Odontinae, and Glaphyriinae of Crambidae; and Gallerinae, Phycitinae, Pyralinae, and Epipaschiinae of Pyralidae. The outgroups chosen for the analysis consists of the dipterans Drosophila melanogaster, Culex quinquefasciatus, and Aedes albopictus. The subfamilies are highlighted according to color. The EFSB is highlighted in bold for emphasis. The PS (Pyraustinae-Spilomelinae), CAMMSS (Crambinae, Acentropinae, Midilinae, Musotiminae, Schoenobiinae, and Scopariinae), and OG (Odontinae-Glaphyriinae) clades are the clades first defined by Regier et al. (2012) [16].
Figure 5. Topological concordance between the partitioned maximum likelihood and Bayesian inference phylogenetic trees. The maximum likelihood (ML) and Bayesian inference (BI) phylogenetic trees are shown on the left- and right-hand side of the labels, respectively. The labeling scheme follows the tip order from the ML tree as indicated by the dotted lines. The subfamilies are coded according to color. The grey lines connecting the tips of both trees indicate tips that agree for both trees, while the red lines indicate the tips with a different order and topology between the trees. SH_aLRT and posterior probability support values are shown for the ML and BI trees, respectively. The subfamilies included in the analysis are Spilomelinae, Pyraustinae, Acentropinae, Schoenobiinae, Crambinae, Odontinae, and Glaphyriinae of Crambidae; and Gallerinae, Phycitinae, Pyralinae, and Epipaschiinae of Pyralidae. The outgroups chosen for the analysis consists of the dipterans Drosophila melanogaster, Culex quinquefasciatus, and Aedes albopictus. The subfamilies are highlighted according to color. The EFSB is highlighted in bold for emphasis. The PS (Pyraustinae-Spilomelinae), CAMMSS (Crambinae, Acentropinae, Midilinae, Musotiminae, Schoenobiinae, and Scopariinae), and OG (Odontinae-Glaphyriinae) clades are the clades first defined by Regier et al. (2012) [16].
Insects 15 00220 g005
Table 2. Summary of the genes for the Leucinodes orbonalis Guenée mitogenome.
Table 2. Summary of the genes for the Leucinodes orbonalis Guenée mitogenome.
GeneLocationStrandStart CodonStop CodonLength (bp)Anticodon
trnM1–69Majority--69CAU
trnI70–134Majority--65GAU
trnQ132–200Minority--69UUG
ND2269–1270MajorityTTGTAA1002-
trnW1285–1351Majority--67UCA
trnC1344–1409Minority--66GCA
trnY1413–1479Minority--67GUA
COX11493–3026MajorityTTGT1534-
trnL23027–3093Majority--67UAA
COX23094–3774MajorityATGTAA681-
trnK3779–3849Majority--71CUU
trnD3850–3918Majority--69GUC
ATP83919–4086MajorityATTTAA168-
ATP64080–4759MajorityATGTA680-
COX34760–5548MajorityATGTAA789-
trnG5551–5616Majority--66UCC
ND35617–5970MajorityATTTAA354-
trnA5975–6038Majority--64UGC
trnR6039–6102Majority--64UCG
trnN6102–6168Majority--67GUU
trnS16171–6236Majority--66GCU
trnE6240–6306Majority--67UUC
trnF6326–6392Minority--67GAA
ND56393–8124MinorityATTT1732-
trnH8125–8194Minority--70GUG
ND48195–9534MinorityATGTA1340-
ND4L9544–9834MinorityATGTAA291-
trnT9840–9905Majority--66UGU
trnP9906–9971Minority--66UGG
ND69974–10507MajorityATTTAA534-
CYTB10511–11657MajorityATGT1147-
trnS211658–11726Majority--69UGA
ND111743–12681MinorityATGTAA939-
trnL112683–12754Minority--72UAG
l-rna12764–14121Minority--1358-
trnV14110–14174Minority--65UAC
s-rna14175–14983Minority--809-
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

Despabiladeras, J.B.; Bautista, M.A.M. Complete Mitochondrial Genome of the Eggplant Fruit and Shoot Borer, Leucinodes orbonalis Guenée (Lepidoptera: Crambidae), and Comparison with Other Pyraloid Moths. Insects 2024, 15, 220. https://0-doi-org.brum.beds.ac.uk/10.3390/insects15040220

AMA Style

Despabiladeras JB, Bautista MAM. Complete Mitochondrial Genome of the Eggplant Fruit and Shoot Borer, Leucinodes orbonalis Guenée (Lepidoptera: Crambidae), and Comparison with Other Pyraloid Moths. Insects. 2024; 15(4):220. https://0-doi-org.brum.beds.ac.uk/10.3390/insects15040220

Chicago/Turabian Style

Despabiladeras, Joshua B., and Ma. Anita M. Bautista. 2024. "Complete Mitochondrial Genome of the Eggplant Fruit and Shoot Borer, Leucinodes orbonalis Guenée (Lepidoptera: Crambidae), and Comparison with Other Pyraloid Moths" Insects 15, no. 4: 220. https://0-doi-org.brum.beds.ac.uk/10.3390/insects15040220

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