Complete Mitochondrial Genome Characterization of Schrankia costaestrigalis (Insecta: Erebidae: Hypenodinae) and Its Phylogenetic Implication

The pinion-streaked snout Schrankia costaestrigalis is a new potato pest that has recently been recorded in China. In this study, we analyzed the complete mitochondrial genome of S. costaestrigalis. The results revealed the mitogenome (GenBank: OQ181231) to occur as a circular DNA molecule of 16,376 bp with 51.001% AT content, including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and 1 control region. Notably, the PCGs exhibited typical ATN (Met) start codons, including cox1, which deviated from the usual CGA start codon observed in other lepidopteran mitogenomes, and followed the conventional TAN stop codons. The 22 tRNA genes demonstrated the ability to form a cloverleaf structure, with the exception of trnS1-NCU, which lacked the DHU arm present in other Erebidae mitogenomes. Additionally, conserved motifs like “ATAGA + poly-T (19 bp) stretch” and five microsatellite-like elements (TA) were identified in the AT-rich region. The phylogenetic trees revealed that the Hypenodinae subfamily forms an independent lineage closely related to Erebinae and Catocalinae. The comprehensive mitogenome of S. costaestrigalis will greatly enhance future studies focused on the molecular classification and phylogenetic understanding of the Hypenodinae subfamily within the larger family Erebidae.


Introduction
The pinion-streaked snout Schrankia costaestrigalis (Stephens, 1834), a species of Lepidoptera in the family Erebidae, subfamily Hypenodinae, belongs to the genus Schrankia Hübner (1825) (Figure 1).It is distributed in Europe, Central and North Asia, North Africa, and Australia [1].This pest was previously defined as a general insect that feeds on various herbaceous and woody plants [2].However, the pest broke out in the potato's main producing area of Yulin City, Guangxi, China in February 2017, affecting an area of 287.34 ha.This was the first report of S. costaestrigalis in China [3].The only other report of this pest-damaging crop was on the broad bean Vicia faba in Japan [4].The damage rate to potato tubers generally ranged between 50% and 80%, with values reaching 100% in some fields.This resulted in a yield loss of 5 million kg and caused a significant impact on the winter potato industry in Guangxi.Furthermore, there has been a rise in potato production in China since 2015, making it the fourth largest crop in the country after rice, wheat, and corn [5].Thus, the occurrence of this pest poses a potential threat to the commercial production of potatoes in China and food security in general.
family.By utilizing a combination of next-generation sequencing (NGS) and Sanger sequencing techniques, we successfully obtained the full sequence.Additionally, we constructed phylogenetic trees using 24 mitogenomes, including 22 mitogenomes of the Erebidae family and 2 outgroup mitogenomes.This analysis enhances our understanding of the phylogenetic placement of S. costaestrigalis within the broader context of the Erebidae family.The findings of this study are valuable for reconstructing the phylogenetic relationships among species within the Erebidae family.

Adult Larva
Genes 2023, 14, x FOR PEER REVIEW 3 of 15 The feeding lava Potato damage

Animal Materials and DNA Extraction
The individuals of S. costaestrigalis were collected from Renhou Village (22.639°N, 110.056°E), Renhou Town, Yuzhou District, Yulin City, Guangxi Zhuang Autonomous Region, China, on 20 March 2022.They were continuously raised in the lab.To extract genomic DNA (gDNA) from the pupae, the Qiagen DNeasy Blood and Tissue Extraction kit (Qiagen, Germantown, MD, USA) were employed.The purified gDNA was detected using a NanoPhotometer ® spectrophotometer (Implen, Los Angeles, CA, USA), while the concentration was determined using a Qubit ® 2.0 fluorometer (Life Technologies, Carlsbad, CA, USA).

Illumina and Sanger Sequencing
Sequencing libraries for the quality-checked gDNA were generated using a TrueLib DNA Library Rapid Prep Kit for Illumina sequencing (Illumina, Inc., San Diego, CA, USA).The libraries were subjected to size distribution analysis using an Agilent 2100 bioanalyzer (Agilent Technologies, Inc., Santa Clara, CA, USA), followed by real-time PCR quantitative testing.The successfully generated libraries were sequenced using an Illumina NovaSeq 6000 platform (Illumina, Inc., San Diego, CA, USA).A total of 150-bp paired-end reads with a 300-bp insert library were generated.Three pairs of primers designed to match generally conserved regions of the published target mtDNA were used to amplify S. costaestrigalis is not typically regarded as an agricultural pest, and research on the pest is lacking [3,4].In particular, the pest's origin and development in China remain unclear.Mitochondria, which are cellular organelles, play a crucial role in cellular energy metabolism and respiration.Studying the mitochondrial genome of a rare agricultural pest can help us to better understand its population genetic structure and genetic diversity, and thus allow us to infer its origin and distribution.Moreover, by comparing the mitochondrial genome of this pest with those of other related species, scientists can also clarify its classification and evolutionary relationships, providing a scientific basis for pest control [6].The interrelationships between major subfamilies and tribes within the Erebidae family were analyzed utilizing molecular data from nuclear and mitochondrial genes.Clade Hypenodinae was lifted to a subfamily status, which made the genera Schrankia Hübner and Luceria Walker be associated more closely [7].And by comparing the mitochondrial genome of Orthaga olivacea Warre with that of other lepidopteran insects, it was confirmed that O. Olivacea belongs to the Pyralidae family, which provides a reference for the pest's control [8].
Insect mitochondrial genomes, known as mitogenomes, are circular, double-stranded DNA molecules that typically range in length from 15 to 19 kb.They comprise a total of 37 genes, encompassing various types including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and a control region (CR) [9].Despite the usefulness of next-generation sequencing (NGS) in mitogenome assembly and phylogenetic analysis of insects [9], it can be challenging to bridge the gap between contigs, particularly in regions with high A+T content.To address this issue, PCR amplification and Sanger sequencing can be employed to successfully fill in the gaps in the mitogenome sequence.
This research focuses on the sequencing and characterization of the mitogenome of S. costaestrigalis, representing the first complete mitogenome within the Hypenodinae subfamily.By utilizing a combination of next-generation sequencing (NGS) and Sanger sequencing techniques, we successfully obtained the full sequence.Additionally, we constructed phylogenetic trees using 24 mitogenomes, including 22 mitogenomes of the Erebidae family and 2 outgroup mitogenomes.This analysis enhances our understanding of the phylogenetic placement of S. costaestrigalis within the broader context of the Erebidae family.The findings of this study are valuable for reconstructing the phylogenetic relationships among species within the Erebidae family.

Animal Materials and DNA Extraction
The individuals of S. costaestrigalis were collected from Renhou Village (22.639 • N, 110.056 • E), Renhou Town, Yuzhou District, Yulin City, Guangxi Zhuang Autonomous Region, China, on 20 March 2022.They were continuously raised in the lab.To extract genomic DNA (gDNA) from the pupae, the Qiagen DNeasy Blood and Tissue Extraction kit (Qiagen, Germantown, MD, USA) were employed.The purified gDNA was detected using a NanoPhotometer ® spectrophotometer (Implen, Los Angeles, CA, USA), while the concentration was determined using a Qubit ® 2.0 fluorometer (Life Technologies, Carlsbad, CA, USA).

Illumina and Sanger Sequencing
Sequencing libraries for the quality-checked gDNA were generated using a Tru-eLib DNA Library Rapid Prep Kit for Illumina sequencing (Illumina, Inc., San Diego, CA, USA).The libraries were subjected to size distribution analysis using an Agilent 2100 bioanalyzer (Agilent Technologies, Inc., Santa Clara, CA, USA), followed by realtime PCR quantitative testing.The successfully generated libraries were sequenced using an Illumina NovaSeq 6000 platform (Illumina, Inc., San Diego, CA, USA).A total of 150-bp paired-end reads with a 300-bp insert library were generated.Three pairs of primers designed to match generally conserved regions of the published target mtDNA were used to amplify short fragments from nad3-nad5, rrnL and the control region (CR) (Supplementary Table S1).The PCR products were cloned into pMD18-T vectors (Takara, Kyoto, Japan) and subsequently sequenced, or they were sequenced directly by the dideoxy nucleotide procedure, using an ABI 3730 automatic sequencer (Applied Biosystems, Foster City, CA, USA) [10] (Supplementary File S1).

Raw Reads Cleaning and Mitogenome Assembly
The high-quality clean reads were filtered from the raw data.This was achieved by employing the fastp tool (version 0.23.2) [11], which is a widely used software for read quality filtering (source code available at: https://github.com/OpenGene/fastp,accessed on 9 January 2023).The quality control (QC) criteria applied to ensure the reliability of the raw reads were as stated: (1) Trimming adapter sequences longer than six bases; (2) Removal of reads with >0 unidentified nucleotides (N); (3) Removal of reads with >20% bases with Phred quality < Q20; (4) Removal of reads with <15 bases.

Phylogenetic Inference
In this study, we performed a phylogenetic analysis using the mitogenomic sequences of 22 species belonging to the Erebidae family.Two additional species were used as outgroup references.(Table 1).Lepisma saccharina and Corydidarum magnifica were selected as the outgroups.Nucleotide sequences of 13 PCGs from mitogenoms were used to construct the phylogenetic relationships within S. costaestrigalis utilizing PhyloSuite version 1.2.2 [24], including MAFFT [25], ModelFinder [26,27], and MrBayes version 3.2.7 [28].Alignment was performed using MAFFT version 7 with the default settings.The nucleotide matrix was used for the phylogenetic analysis with two methods: Bayesian inference (BI) in MrBayes 3.2.7 and the maximum likelihood (ML) method with MEGA v11.0.13 [21].The best-fit edge-unlinked partition model for BI was selected using PhyloSuite version 1.2.2 [24] with the Bayesian Information Criterion (BIC), employing ModelFinder [26,27] (Supplementary Table S2).A BI analysis was conducted for each matrix (two parallel runs, 20,000,000 generations).The initial 25% of the generated trees were discarded as burn-in, and the average standard deviation of split frequencies for the remaining trees was 0.002572 (<0.01), indicating convergence.The ML tree was constructed using MEGA v11.0.13 [21], and the best-fit model, as determined by BIC scores, was the General Time Reversible (GTR) model coupled with a discrete γ distribution (+G) consisting of 5 rate categories.For the ML analysis, 1000 bootstrap resampling replicates were performed to assess node support values.The generated phylogenetic trees were visualized using the Interactive Tree Of Life (iTOL) tool (source: https://itol.embl.de/)(accessed on 9 January 2023) [29].From a 300 bp insert library, a total of approximately 5.64 Gb raw reads were generated.Subsequently, the fastp software [9] was employed to obtain approximately 5.37 Gb of high-quality clean reads.The Q20 (percentage of bases with quality value ≥ 20), Q30 (percentage of bases with quality value ≥ 30), and G+C content values of the clean reads were 97.82%, 92.84%, and 34.84%, respectively (Table 2).Using NOVOPlasty, the high-quality cleaned short reads, which accounted for 0.20% of the total reads from the mitogenome, allowed for the near-complete assembly of the S. costaestrigalis (OQ181231.1)mitogenome.This assembly achieved 100% coverage of the mitogenome with a high average-read depth of 675 times.By reducing repetitive sequences, the assembly process was able to generate a comprehensive and accurate representation of the S. costaestrigalis mitogenome [12].This contains the mitochondrial sequence results from four possible combinations, ranging in length from 16,094 to 16,101 bp.Although the circular mitochondrial genome sequence was successfully assembled using GetOrganelle [14], the length of the sequence was 15,550 bp, and filling the gap in CR with high A+T content between contigs proved to be difficult.We sequenced the three regions using sanger sequencing, including nad3-nad5, rrnL, and CR, and manually assembled a complete mitogenome which consists of traditional circular DNA molecules.The mitogenome exhibited the longest length of 16,376 bp in the family Erebidae (Table 1).The mitogenome of S. costaestrigalis contains 39.57% A, 41.68% T, 7.32% G, and 11.43% C, showing an obvious AT bias with a 81.25% A+T content, which was slightly lower compared to Dysgonia stuposa [38].Specifically, the major strand of the S. costaestrigalis mitogenome exhibited an AT-skew of −0.026 and a GC-skew of −0.219.These values indicated a compositional bias on the major strand, with a slight excess of T nucleotides over A nucleotides, and a strong excess of C nucleotides over G nucleotides.The AT and GC bias were similar to other mitogenomes, in the family Erebidae, such as Dysgonia stuposa [38] and Hydrillodes repugnalis [40].
The S. costaestrigalis mitogenome comprised 13 PCGs, 1 CR, and 22 tRNA genes, and 2 rRNA genes (Figure 1).The arrangement and orientation of genes in the S. costaestrigalis mitogenome were consistent with other mitogenomes in the family Erebidae.On the majority strand (J-strand), 23 genes, including 9 PCGs and 14 tRNAs, were encoded.On the minority strand (N-strand), there were four PCGs, eight tRNAs, and two rRNAs (Figure 2 and Table 3).Within these genes, there were 12 instances of overlap, totaling 37 bp, with the longest overlapping region found between trnW-UCA and trnC-GCA (Table 3).Additionally, there were 17 intergenic spacer regions, spanning a total of 113 bp, with the longest spacer occurring between trnS1-GCU and trnE-UUC (Table 3).
regions using sanger sequencing, including nad3-nad5, rrnL, and CR, and manually assembled a complete mitogenome which consists of traditional circular DNA molecules.The mitogenome exhibited the longest length of 16,376 bp in the family Erebidae (Table 1).The mitogenome of S. costaestrigalis contains 39.57% A, 41.68% T, 7.32% G, and 11.43% C, showing an obvious AT bias with a 81.25% A+T content, which was slightly lower compared to Dysgonia stuposa [38].Specifically, the major strand of the S. costaestrigalis mitogenome exhibited an AT-skew of −0.026 and a GC-skew of −0.219.These values indicated a compositional bias on the major strand, with a slight excess of T nucleotides over A nucleotides, and a strong excess of C nucleotides over G nucleotides.The AT and GC bias were similar to other mitogenomes, in the family Erebidae, such as Dysgonia stuposa [38] and Hydrillodes repugnalis [40].
The S. costaestrigalis mitogenome comprised 13 PCGs, 1 CR, and 22 tRNA genes, and 2 rRNA genes (Figure 1).The arrangement and orientation of genes in the S. costaestrigalis mitogenome were consistent with other mitogenomes in the family Erebidae.On the majority strand (J-strand), 23 genes, including 9 PCGs and 14 tRNAs, were encoded.On the minority strand (N-strand), there were four PCGs, eight tRNAs, and two rRNAs (Figure 2 and Table 3).Within these genes, there were 12 instances of overlap, totaling 37 bp, with the longest overlapping region found between trnW-UCA and trnC-GCA (Table 3).Additionally, there were 17 intergenic spacer regions, spanning a total of 113 bp, with the longest spacer occurring between trnS1-GCU and trnE-UUC (Table 3).Abbreviations: J, J-strand (the majority strand); N, N-strand (the minority strand).
The analysis of relative synonymous codon usage (RSCU) was performed on the 13 protein-coding genes, which consisted of a total of 3711 codons excluding the start and stop codons.The RSCU analysis unveiled the presence of codon usage bias within the S. costaestrigalis mitogenome (refer to Figure 3 and Supplementary Table S3).Furthermore, when examining the frequency of amino acids (as shown in Supplementary Table S3), it was observed that Leu was the most abundant (558 occurrences), followed by Ile (455) and Phe (365).Looking specifically at the codon usage counts, it was found that UUA (482) was the predominant codon for Leu, AUU (435) followed Ile, and UUU (343) followed Phe.While most codons were present in the S. costaestrigalis mitogenome, there were a few exceptions.Specifically, the codons CUG for Leu, CCG for Pro, and AGG for Ser1 were absent, indicating a lack of GC-rich synonymous codons with G at the third codon position [38].Furthermore, the thirteen protein-coding genes displayed a preference for A and T nucleotides, as outlined in Supplementary Table S4.
The analysis of relative synonymous codon usage (RSCU) was performed on the 13 protein-coding genes, which consisted of a total of 3711 codons excluding the start and stop codons.The RSCU analysis unveiled the presence of codon usage bias within the S. costaestrigalis mitogenome (refer to Figure 3 and Supplementary Table S3).Furthermore, when examining the frequency of amino acids (as shown in Supplementary Table S3), it was observed that Leu was the most abundant (558 occurrences), followed by Ile (455) and Phe (365).Looking specifically at the codon usage counts, it was found that UUA (482) was the predominant codon for Leu, AUU (435) followed Ile, and UUU (343) followed Phe.While most codons were present in the S. costaestrigalis mitogenome, there were a few exceptions.Specifically, the codons CUG for Leu, CCG for Pro, and AGG for Ser1 were absent, indicating a lack of GC-rich synonymous codons with G at the third codon position [38].Furthermore, the thirteen protein-coding genes displayed a preference for A and T nucleotides, as outlined in Supplementary Table S4.

Transfer and Ribosomal RNA Genes
The 22 tRNA genes in the mitochondrial genome of S. costaestrigalis were distributed among the protein-coding genes (PCGs).Among these tRNAs, 14 were encoded on the majority strand, while 8 were encoded on the minority strand (Figure 2 and Table 3).Their lengths varied between 64 bp (trnT-UGU) and 73 bp (trnQ-UUG and trnK-CUU) (Table 3 and Supplementary Table S4), and they all possessed the characteristic cloverleaf secondary structure (depicted in Figure 4).The secondary structures of most tRNAs closely reflected similarity with the one observed in Dysgonia stuposa [38] and Hydrillodes repugnalis [40] mitogenomes.However, trnS1-GCU in the S. costaestrigalis mitogenome exhibited a different structure, featuring a dihydrouridine (DHU) stem (Figure 4) of 3 bp.The diversity in tRNA secondary structures reflected the evolutionary variations among species.The length of the anticodon tRNA stems ranged from 4 bp (trnI-GAU, trnK-CUU, trnH-GUG, and trnS2-UGA) to 9 bp (trnS1-GCU) (Figure 4).The DHU stem length varied from 2 bp (trnF-GAA and trnL1-UAG) to 6 bp (trnP-UGG) (Figure 5), with most falling between 3 and 4 bp.The TΨC stem length varied between 3 bp (trnN-GUU) to 8 bp (trnL1-UAG) (Figure 4), with the majority being between 4 and 5 bp.Two types of mismatched base pairs were identified in tRNAs, namely A-G base pairs and non-canonical G-U base pairs (Figure 5).The DHU stem of trnP-UGG and the TΨC stem of trnC-GCA contained A-G base pairs (Figure 5).The anticodon stems of trnC-GCA, trnS1-GCU, trnF-GAA, and trnT-UGU; the amino acid acceptor stems of trnW-UCA, trnC-GCA, trnA-UGC, and trnL1-UAG; the DHU stem of trnL2-UAA; and G-U base pairs were present in the TΨC stem of trnS1-GCU (Figure 5).The rrnL and rrnS were 1400 and 781 bp long, respectively, and were encoded on the N strand.The rrnL was located between trnL1 and trnV, while rrnS resided between trnV and CR.This is similar to that of the Dysgonia stuposa and Hydrillodes repugnalis mitogenomes [38,40].The AU-skews of the rrnL and rrnS were 0.062 and 0.063, respectively, and corresponding GC-skews were determined as 0.406 and 0.345, respectively (Supplementary Table S4).The rrnL and rrnS showed a markedly high A+T content (83.64 and 85.53%, respectively) bias (Supplementary Table S4), similar to that of the Laelia suffusa mitogenome [44].

Control Region
The CR, also denoted as the AT-rich region or D-loop, had a length of 1421 base pairs and exhibited a high AT content of 91.77%.It was situated between the rrnS and trnM-CAU genes (Figure 2 and Supplementary Table S3).The AT content of the CR in the mitogenome of S. costaestrigalis was marginally greater compared to that of Leucoma salicis [45].Analysis of the CR revealed compositional biases on the majority strand.Specifically, the AT-skew was calculated to be −0.051,indicating a slight excess of T nucleotides over A nucleotides.Similarly, the GC-skew was determined to be −0.197,indicating a significant excess of C nucleotides over G nucleotides.These findings suggest a notable compositional bias in favor of certain nucleotides on the major strand of the CR.The conserved structures that connected the motif "ATAGA + poly-T (19 bp) stretch" [31,36,44,45,49] were located right flanking of the rrnS gene within the CR (Figure 5 and Supplementary File S2).TR1, TR3, TR6, TR8, and TR10 were varied and typical microsatellite-like elements (TA) (Supplementary Table S5).

The Construction of Phylogenetic Trees
The phylogenetic trees were constructed using nucleotide sequences of the 13 PCGs from 24 mitogenomes (as depicted in Figures 6 and 7).Although the BI tree demonstrated considerably higher support values compared to the ML tree in the same dataset, the trees exhibited significantly low support values in the same branches.In the two trees, nine subfamilies were clustered as ((((Arctiinae + (Herminiinae + Aganainae)) + ((Erebinae + Catocalinae) + Hypenodinae) + Calpinae) + Hypeninae) + Lymantriinae) (Figures 6 and 7), which mirrored earlier phylogenetic investigations [50].At the subfamily level, the subfamily Hypenodinae is an independent lineage and is closely related to the subfamilies Erebinae and Catocalinae based on the topology of the trees (Figures 6 and 7).As is known, the species of the subfamily Hypenodinae are small-sized moths and distributed only in the world [51].They were placed as a tribe before, but had a kinship with the subfamily Erebinae and Catocalinae [7].These results will contribute to establishing a molecular framework for the classification and phylogeny within the Erebidae family, specifically focusing on the Hypenodinae subfamily.

Conclusions
The pinion-streaked snout S. costaestrigalis is distributed in Europe, Central and North Asia, North Africa, and Australia [1], yet information on its mitogenome molecular

Conclusions
The pinion-streaked snout S. costaestrigalis is distributed in Europe, Central and North Asia, North Africa, and Australia [1], yet information on its mitogenome molecular

Conclusions
The pinion-streaked snout S. costaestrigalis is distributed in Europe, Central and North Asia, North Africa, and Australia [1], yet information on its mitogenome molecular phy-logenetic is lacking.Therefore, in the current study, we assembled the S. costaestrigalis mitogenome.This was the first complete mitogenome in the subfamily Hypenodinae and is observed to have similar structural characteristics and nucleotide composition, compared to other previously reported mitogenomes of the family Erebidae.We identified 13 PCGs, 1 CR, and 22 tRNA genes, and 2 rRNA genes in our assembled mitogenome.All of the 13 protein-coding genes were initiated with a common start codon, typically ATN (Met), including cox1, which was initiated with the CGA start codon in most of the lepidopteran mitogenomes.These genes were terminated using the standard TAN stop codons.The 22 tRNAs exhibit a characteristic cloverleaf structure, typical for mitochondrial genomes.Notably, trnS1-NCU lacks the DHU arm, distinguishing it from other mitogenomes within the family Erebidae.Additionally, the CR presents a conserved motif "ATAGA + poly-T (19 bp) stretch".Phylogenetic analysis has revealed that the subfamily Hypenodinae forms an independent lineage, closely related to the subfamilies Erebinae and Catocalinae.Given the diverse nature of the Erebidae family and the current limitations in mitogenome data, a more comprehensive understanding of the phylogeny within the family Erebidae would necessitate the inclusion of additional mitogenomes.

Figure 1 .
Figure 1.The images of S. costaestrigales and its damage.

Figure 1 .
Figure 1.The images of S. costaestrigales and its damage.

Figure 2 .
Figure 2. Mitogenome pattern map of S. costaestrigalis.Arrows indicated the orientation of gene transcription.

Genes 2023 , 15 Figure 4 .
Figure 4. Secondary structure of 22 tRNAs from the S. costaestrigalis mitogenome.Mismatched base pairs were visually highlighted in yellow, while matched base pairs were shown in green.Bases within loops were represented in blue, while the nucleotide outlines of the anticodon were indicated in red.

Figure 4 .
Figure 4. Secondary structure of 22 tRNAs from the S. costaestrigalis mitogenome.Mismatched base pairs were visually highlighted in yellow, while matched base pairs were shown in green.Bases within loops were represented in blue, while the nucleotide outlines of the anticodon were indicated in red.

Figure 4 .
Figure 4. Secondary structure of 22 tRNAs from the S. costaestrigalis mitogenome.Mismatched base pairs were visually highlighted in yellow, while matched base pairs were shown in green.Bases within loops were represented in blue, while the nucleotide outlines of the anticodon were indicated in red.

Figure 5 .
Figure 5. Alignment of initiation site for the control region of four species mitogenomes.The Yellow indicates the conserved motif ATAGA.Figure 5. Alignment of initiation site for the control region of four species mitogenomes.The Yellow indicates the conserved motif ATAGA.

Figure 5 .
Figure 5. Alignment of initiation site for the control region of four species mitogenomes.The Yellow indicates the conserved motif ATAGA.Figure 5. Alignment of initiation site for the control region of four species mitogenomes.The Yellow indicates the conserved motif ATAGA.

Figure 6 .
Figure 6.The phylogenetic tree was inferred from the nucleotide sequences of 13 mitogenome PCGs using BI methods.Red-colored node indicators represent the posterior probability values on the Bayesian inference (BI) phylogenetic tree.Red star indicates the newly determined S. costaestrigalis.Lepisma saccharina and Corydidarum magnifica are outgroups.

Figure 7 .
Figure 7.The phylogenetic tree was inferred from the nucleotide sequences of 13 mitogenome PCGs using ML methods.The nodes on the tree are represented with red-colored ML bootstrap support values.Red star indicates the newly determined S. costaestrigalis, Lepisma saccharina and Corydidarum magnifica as outgroups.

Figure 6 .
Figure 6.The phylogenetic tree was inferred from the nucleotide sequences of 13 mitogenome PCGs using BI methods.Red-colored node indicators represent the posterior probability values on the Bayesian inference (BI) phylogenetic tree.Red star indicates the newly determined S. costaestrigalis.Lepisma saccharina and Corydidarum magnifica are outgroups.

Genes 2023 , 15 Figure 6 .
Figure 6.The phylogenetic tree was inferred from the nucleotide sequences of 13 mitogenome PCGs using BI methods.Red-colored node indicators represent the posterior probability values on the Bayesian inference (BI) phylogenetic tree.Red star indicates the newly determined S. costaestrigalis.Lepisma saccharina and Corydidarum magnifica are outgroups.

Figure 7 .
Figure 7.The phylogenetic tree was inferred from the nucleotide sequences of 13 mitogenome PCGs using ML methods.The nodes on the tree are represented with red-colored ML bootstrap support values.Red star indicates the newly determined S. costaestrigalis, Lepisma saccharina and Corydidarum magnifica as outgroups.

Figure 7 .
Figure 7.The phylogenetic tree was inferred from the nucleotide sequences of 13 mitogenome PCGs using ML methods.The nodes on the tree are represented with red-colored ML bootstrap support values.Red star indicates the newly determined S. costaestrigalis, Lepisma saccharina and Corydidarum magnifica as outgroups.

Table 1 .
The sequences of 24 mitogenomes employed for the construction of phylogenetic trees.

Table 2 .
Summary of sequencing reads for S. costaestrigalis.