The Complete Mitochondrial Genomes of Aelia sibirica and A. fieberi (Hemiptera, Pentatomidae), and Phylogenetic Implications

: Species of genus Aelia are important pests of wheat crops in arid areas. In this study, the mitogenomes of A. sibirica and A. fieberi were sequenced using high-throughput sequencing technology. The mitochondrial genome characteristics of both Aelia species were compared and analyzed, and the phylogenetic relationships of Pentatomidae were constructed based on protein-coding genes. In addition, the taxonomic status of the genus Aelia was confirmed. The results showed that the total length of the mitogenome sequences of A. sibirica and A. fieberi were 15,372 bp and 15,450 bp, respectively, including 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and a control region. By comparing the mitochondrial genome structure, base composition, codon usage, RNA secondary structure, and other characteristics, it was found that the mitochondrial genome characteristics of the two species were similar. Phylogenetic analysis showed that Phyllocephali-nae and Asopinae both formed monophyletic groups, but the relationship between Podopinae and Pentatominae was not resolved. Within the subfamily Pentatominae, (Nezarini + Antestiini), (Aeliini + Carpocorini), and (Strachiini + Pentatoma ) formed stable clades. Aelia sibirica and A. fieberi were found to be a stable sibling pair, and the clade was closely related to Dolycoris baccarum .


Introduction
The mitochondrial genome is characterized by its compact size, genetic stability, rapid evolutionary rate, and maternal inheritance, and has emerged as one of the most commonly used molecular markers in evolutionary studies [1][2][3][4][5][6].In hemipteran insects, the mitochondrial genome typically consists of a circular double-stranded DNA ranging from 14 to 20 kb and encompassing 13 protein-coding genes (PCGs), 22 tRNA genes, 2 rRNA genes, and a control region [7][8][9].The known structure of hemipteran mitochondrial genomes exhibits a high degree of conservation, with a tight arrangement of genes and only sporadic occurrences of tRNA gene rearrangements observed among certain species [10].The overall gene order remains largely consistent with that of ancestral insects [11,12].The advent of high-throughput sequencing technologies has allowed the mitochondrial genome to be widely used in molecular evolutionary studies, phylogenetic analyses, and investigations pertaining to biogeography within hemipteran taxa [13][14][15].
The insect family Pentatomidae comprises a diverse range of species, and is the largest family in the superfamily Pentatomoidea.It is widely distributed across the world, with nearly 5000 species belonging to 896 genera and ten subfamilies [16].The classification of higher taxonomic ranks within Pentatomidae has been the source of considerable disagreement among researchers, despite the increasing number of studies and abundant data.Based on morphological and anatomical characteristics, Pentatomidae has been successively divided into 11, 9, 8, and 10 subfamilies [17][18][19][20].Originally, seven subfamilies were identified for Pentatomidae species in China; however, with the development of The adult specimens of A. sibirica used in this experiment were collected on 14 August 2019, from the Sanjiangyuan National Nature Reserve in Qinghai Province.The adult specimens of A. fieberi were collected on 10 August 2019 from Shanfanao in Henan Province.Specimens were preserved in 100% ethanol and stored at −20 • C until used for DNA extraction.The leg muscle tissues of A. sibirica and A. fieberi were extracted, and genomic DNA extraction was performed following the instructions of the TIANamp Genomic DNA Kit (Tiangen biotechnology, Beijing, China).Voucher specimens were deposited at the Animal Systematics Laboratory, Department of Biology, Xinzhou Normal University, Xinzhou, China.The quality of the extracted DNA was assessed using a NanoDrop 2000 nucleic acid analyzer (Thermo Fisher Scientific, Waltham, MA, USA).The mitogenome was sequenced using the whole-genome shotgun method on an Illumina Miseq platform (Personalbio, Shanghai, China).

Genome Annotations and Sequence Analyses
The sequencing data was assembled using A5-miseq v. 20150522 [31] and SPAdes v. 3.9 [32] software.The mitochondrial genome sequences were annotated by Geneious 10.1.3[33].The annotation of protein-coding genes was conducted by Open Reading Frame Finder on the NCBI website (http://www.ncbi.nlm.nih.gov/gorf/gorf.html,accessed on 13 December 2021), as well as identifying start and stop codons.The tRNA genes were automatically annotated by MITOS (http://mitos.bioinf.uni-leipzig.de/index.py/,accessed on 13 December 2021) [34] with the invertebrate mitochondrial code.The annotation of rRNA genes was based on published reference sequence (Eurydema dominulus (Scopoli, 1763)) and combined with their putative secondary structures.The locations of the control region were identified by the boundary of neighboring genes.The newly sequenced and annotated mitogenomes were submitted to GenBank with accession numbers NC_067883 (A. sibirica) and ON059969 (A. fieberi).The mitochondrial genome maps were generated by GCView Server [35].The base composition and relative synonymous codon usage (RSCU) were performed using MEGA 7.0 [36].Additionally, AT skew was calculated using the formula (A − T)/(A + T), while GC skew was calculated using (G − C)/(G + C) [37].The number of tandem repeats of control region was investigated with Tandem Repeats Finder (http://tandem.bu.edu/trf/trf.html,accessed on 20 March 2022) [38].

Phylogenetic Analysis
In order to investigate the phylogenetic relationships among different taxonomic levels within Pentatomidae, and the phylogenetic relationship between the genus Aelia and its closely related genera, this study selected 27 representative species from the family Pentatomidae to construct phylogenetic trees (Table 1).The selected 27 taxa covered 4 subfamilies of Pentatomidae and 13 tribes of Pentatominae which existed in the NCBI database.In addition, Sastragala esakii (Pentatomoidea: Acanthosomatidae) was selected as the outgroup.The maximum likelihood method (ML) and Bayesian inference (BI) were used to construct phylogenetic trees based on 13 protein-coding genes.Multiple alignments of protein-coding gene sequences were performed using MEGA 7.0, and the aligned gene sequence matrices were concatenated using the SequenceMatrix 1.8 [39].The BI tree was generated using MrBayes v. 3.2 [40], with the GTR + I + G model selected and run for 10,000,000 generations.The program was stopped when the convergence diagnostic value was below 0.01.The ML tree was constructed using RAxML v. 7.0.3[41], and branch support was assessed using 1000 bootstrap replicates.

Genomic Features
The mitochondrial genomes of two Aelia species were obtained using high-throughput sequencing.The complete mitochondrial genome sequence of A. sibirica was 15,372 bp (Figure 1A), whereas that of A. fieberi was 15,450 bp (Figure 1B).Both genomes consisted of 13 PCGs, 22 tRNA, and two rRNA genes and a control region (Table 2).This gene arrangement was consistent with the typical order observed in Heteroptera insects without any specific rearrangements.The N-strand encoded fourteen genes, including four proteincoding, nad1, nad4, nad4L, and nad5; eight tRNA, trnQ, trnC, trnY, trnF, trnH, trnP, trnL1, and trnV; and the 12S rRNA and 16S rRNA.The remaining genes were encoded on the J-strand.

Genomic Features
The mitochondrial genomes of two Aelia species were obtained using high-throughput sequencing.The complete mitochondrial genome sequence of A. sibirica was 15,372 bp (Figure 1A), whereas that of A. fieberi was 15,450 bp (Figure 1B).Both genomes consisted of 13 PCGs, 22 tRNA, and two rRNA genes and a control region (Table 2).This gene arrangement was consistent with the typical order observed in Heteroptera insects without any specific rearrangements.The N-strand encoded fourteen genes, including four protein-coding, nad1, nad4, nad4L, and nad5; eight tRNA, trnQ, trnC, trnY, trnF, trnH, trnP, trnL1, and trnV; and the 12S rRNA and 16S rRNA.The remaining genes were encoded on the J-strand.
The mitochondrial genome of A. sibirica contained ten overlapping regions, ranging from 1 to 8 bp, and seventeen intergenic spacer regions, ranging from 1 to 27 bp.In contrast, the mitochondrial genome of A. fieberi had eight overlapping regions, ranging from 1 to 8 bp, and nineteen intergenic spacer regions, ranging from 1 to 24 bp.The mitochondrial genome of A. sibirica contained ten overlapping regions, ranging from 1 to 8 bp, and seventeen intergenic spacer regions, ranging from 1 to 27 bp.In contrast, the mitochondrial genome of A. fieberi had eight overlapping regions, ranging from 1 to 8 bp, and nineteen intergenic spacer regions, ranging from 1 to 24 bp.

Nucleotide Composition and Codon Usage
In the complete mitogenomes, nucleotide composition analysis revealed that both A. sibirica and A. fieberi had significantly higher AT than GC content.The nucleotide composition in the mitochondrial genomes of both species were similar, and the overall base composition was identical (Table 3).The AT content in the control region was slightly lower than that in the complete mitochondrial genome, PCGs, rRNAs, and tRNAs with values of 70.5% and 69.2%, respectively.Among the PCGs, nad2 had the highest AT content in A. sibirica (79.8%), and nad6 had the highest AT content in A. fieberi (80.2%).The gene with the lowest AT content was cox1 (67.6% and 67.0%, respectively).The AT-skew values of the complete mitochondrial genome, 22 tRNAs, and the control region were positive in both Aelia species, whereas the AT-skew values for the two rRNAs and thirteen PCGs were negative.Among the PCGs, AT-skew values were negative for nad1, nad4, nad4L, nad5, and cytb in both Aelia species; nad6 also had a negative AT-skew value in A. fieberi.The GC-skew values for the complete mitochondrial genome, thirteen PCGs, and the control region were negative, whereas the GC-skew values for the two rRNAs and twenty-two tRNAs were positive.Among the PCGs, the AT-skew values were negative, and the GC-skew values were positive for nad1, nad4, nad4L, nad5, and cytb in both Aelia species.

PCGs
The lengths of the 13 PCGs in A. sibirica and A. fieberi were 11,033 and 11,082 bp, respectively, accounting for 71.77% and 71.73% of the complete mitochondrial genome sequence, respectively, with the exception of cox1, nad1, and atp8, which started with the TTG, TTG/GTG, and TTG codons, respectively; the genes in both species started with the ATN codon.In A. sibirica, the numbers of genes starting with the ATC, ATG, ATT, and ATA codons were one, five, two, and two, respectively; whereas in A. fieberi, the numbers of genes starting with the ATC, ATG, ATT, and ATA codons were one, four, two, and three, respectively.With the exception of nad3 and nad5, the PCGs in both species had identical stop codons.The preferred stop codon was TAA; however, some genes ended with TAG or an incomplete T.

Ribosomal and Transfer RNA Genes
The mitogenomes of A. sibirica and A. fieberi contained 22 tRNA genes with total lengths of 1458 and 1453 bp, respectively.The length ranges of the tRNA genes in A. sibirica were from 62 bp (trnC) to 71 bp (trnK), whereas in A. fieberi, they ranged from 61 bp (trnR) to 70 bp (trnK).A comparison of the tRNA secondary structures found minimal differences between the two Aelia species.Apart from trnS1 and trnV, which lacked a dihydrouridine arm, the remaining 20 tRNA genes exhibited a typical cloverleaf structure (Figure 3).The number of nonclassical pairings of A. sibirica and A. fieberi was 18 and 26, respectively, with multiple occurrences of G-U non-canonical pairing.
The lengths of 12S rRNA and 16S rRNA genes in A. sibirica were 796 and 1282 bp, respectively, whereas those in A. fieberi were 795 and 1272 bp, respectively.The 16S rRNA gene was located between trnL1 and trnV, whereas the 12S rRNA gene was located between trnV and the control region.The secondary structures of the rRNA in A. sibirica were compared with those in A. fieberi and identified five domains in the 16S rRNA, whereas the 12S rRNA contained three domains (Figures 4 and 5).The 12S rRNA secondary structure in both Aelia species showed minimal changes in the structural domains and nucleotide positions, particularly in the loop region.In contrast, domains IV and VI in the secondary structure of 16S rRNA were more conserved than domains I, II, and V, with variable sites concentrated in various loop regions.

Ribosomal and Transfer RNA Genes
The mitogenomes of A. and A. fieberi contained 22 tRNA genes with total lengths of 1458 and 1453 bp, respectively.The length ranges of the tRNA genes in A. sibirica were from 62 bp (trnC) to 71 bp (trnK), whereas in A. fieberi, they ranged from 61 bp (trnR) to 70 bp (trnK).A comparison of the tRNA secondary structures found minimal differences between the two Aelia species.Apart from trnS1 and trnV, which lacked a dihydrouridine arm, the remaining 20 tRNA genes exhibited a typical cloverleaf structure (Figure 3).The number of nonclassical pairings of A. sibirica and A. fieberi was 18 and 26, respectively, with multiple occurrences of G-U non-canonical pairing.The lengths of 12S rRNA and 16S rRNA genes in A. sibirica were 796 and 1282 bp, respectively, whereas those in A. fieberi were 795 and 1272 bp, respectively.The 16S rRNA gene was located between trnL1 and trnV, whereas the 12S rRNA gene was located between trnV and the control region.The secondary structures of the rRNA in A. sibirica were compared with those in A. fieberi and identified five domains in the 16S rRNA, whereas the 12S rRNA contained three domains (Figures 4 and 5).The 12S rRNA secondary structure in both Aelia species showed minimal changes in the structural domains and nucleotide positions, particularly in the loop region.In contrast, domains IV and VI in the secondary structure of 16S rRNA were more conserved than domains I, II, and V, with variable sites concentrated in various loop regions.

Control Region
Complete sequences of the control region in A. sibirica and A. fieberi were obtained, resulting in lengths of 717 bp and 768 bp, respectively.The control region was located between the 12S rRNA gene and trnI.The control region of A. sibirica had an AT content of 70.5%, whereas the control region of A. fieberi had an AT content of 69.2%, both of which were lower than the overall AT content of the complete mitogenomes.The repeat unit of A. fieberi was 54 bp, with five copies located between 462 and 731 bp; the repeat unit of A. sibirica was 55 bp, with four copies located between 487 and 706 bp.

Phylogenetic Analyses
Phylogenetic trees of the Pentatomidae were constructed based on 13 PCGs using both the ML and BI methods (Figure 6).The topological structures of the phylogenetic trees were consistent, and most internal nodes had high posterior values.The results showed that the subfamily Phyllocephalinae was the earliest diverging lineage and formed the basal clade of the Pentatomidae.Phyllocephalinae were a sister group to Placosternum, currently classified in Pentatominae.Both the subfamilies Phyllocephalinae and Asopinae were found to be monophyletic and showed clear relationships within the subfamilies, supported by high branch support values.In the subfamily Podopinae, Scotinophara lurida, Graphosoma rubrolineata, and Deroploa parva did not cluster into a monophyletic group, but formed sister-group relationships with Pentatominae species.Within the subfamily Pentatominae, some tribes formed stable sibling clades, such as (Aeliini + Carpocorini), (Nezarini + Antestiini), and (Strachiini + Pentatoma).

Control Region
Complete sequences of the control region in A. sibirica and A. fieberi were obtained, resulting in lengths of 717 bp and 768 bp, respectively.The control region was located between the 12S rRNA gene and trnI.The control region of A. sibirica had an AT content of 70.5%, whereas the control region of A. fieberi had an AT content of 69.2%, both of which were lower than the overall AT content of the complete mitogenomes.The repeat unit of A. fieberi was 54 bp, with five copies located between 462 and 731 bp; the repeat unit of A. sibirica was 55 bp, with four copies located between 487 and 706 bp.

Phylogenetic Analyses
Phylogenetic trees of the Pentatomidae were constructed based on 13 PCGs using both the ML and BI methods (Figure 6).The topological structures of the phylogenetic trees were consistent, and most internal nodes had high posterior values.The results showed that the subfamily Phyllocephalinae was the earliest diverging lineage and formed the basal clade of the Pentatomidae.Phyllocephalinae were a sister group to Placosternum, currently classified in Pentatominae.Both the subfamilies Phyllocephalinae and Asopinae were found to be monophyletic and showed clear relationships within the subfamilies, supported by high branch support values.In the subfamily Podopinae, Scotinophara lurida, Graphosoma rubrolineata, and Deroploa parva did not cluster into a monophyletic group, but formed sister-group relationships with Pentatominae species.Within the subfamily Pentatominae, some tribes formed stable sibling clades, such as (Aeliini + Carpocorini), (Nezarini + Antestiini), and (Strachiini + Pentatoma).

Discussion
The sequence lengths of insect mitochondrial genomes are relatively stable and primarily influenced by the length of the control region [42,43].In this study, the complete mitochondrial genome sequences of two species belonging to the genus Aelia were obtained, with lengths of 15,372 and 15,450 bp.The lengths of the control regions were 717 and 768 bp, showing only minor differences.The 37 genes were tightly arranged, and their genomic structure was consistent with that of most Pentatomidae insects, without gene

Discussion
The sequence lengths of insect mitochondrial genomes are relatively stable and primarily influenced by the length of the control region [42,43].In this study, the complete mitochondrial genome sequences of two species belonging to the genus Aelia were obtained, with lengths of 15,372 and 15,450 bp.The lengths of the control regions were 717 and 768 bp, showing only minor differences.The 37 genes were tightly arranged, and their genomic structure was consistent with that of most Pentatomidae insects, without gene rearrangements or losses [14,[44][45][46].The longest overlapping region was located between trnC and trnY, and the longest intergenic region was located between nad1 and trnL1, which is a common feature of the mitochondrial genomes of Pentatomidae [47].The A + T content of both species of Aelia was significantly higher than the G + C content, indicating a clear AT bias.For most genes, the usage of T exceeded that of A bases, whereas the usage of G exceeded that of C bases.The compositional bias of bases has important implications for studying the mechanisms of mitochondrial genome replication, transcription, and phylogenetic relationships [48].
In addition to TTG and GTG, other start codons are preferred by Pentatomidae [49,50], and ATN is the most commonly used start codon in Aelia species.In the present study, cox1, nad1, and atp8 used TTG, TTG/GTG, and TTG start codons, respectively.In total, 11 of the 13 PCGs had identical stop codons, and the most frequently used stop codon was TAA, followed by TAG and incomplete T. PCGs ending with a T codon are completed with a TAA stop codon during transcription through polyadenylation [51].The AT bias in the mitochondrial genome was also reflected in the frequency of synonymous codon usage by PCGs.Codons with higher frequencies often ended with an A or U, and their relative synonymous codon usage values were greater than one.
The lengths and nucleotide compositions of RNA genes showed high similarity between A. sibirica and A. fieberi.The secondary structures of rRNA genes play a crucial role in the functionality of PCGs.Most nucleotide position variations existed in the ring region, whereas the stem regions tended to be more conserved.In addition, the 12S rRNA gene was more conserved than the 16S rRNA gene.G-U pairing was observed multiple times in the secondary structure of tRNA, and such mismatches have been shown to play an important role in maintaining secondary structure stability [52].
In the phylogenetic analysis, there were relatively few species from Asopinae and Phyllocephalinae, and they exhibited distinct differences from other subfamilies within the family Pentatomidae, making their classification less controversial and more stable.The phylogenetic relationships obtained by different methods supported the monophyly of Asopinae and Phyllocephalinae, which was consistent with the results based on mitochondrial RNA genes and aligned with current mainstream views [20,53].
Species in the subfamily Pentatominae are abundant, and their morphological variations are great.Owing to the lack of stable and reliable identification characters for phylogenetic reconstructions, different researchers have had different opinions on the division of tribes and the monophyly of Pentatominae, and the relationships between the various groups within Pentatominae remain uncertain [24,54].Based on morphological studies, Yang divided the Chinese Pentatominae species into seven tribes [55].Whereas Rider recorded 42 tribes about the worldwide Pentatominae species [18].In the 7-tribe classification system, the genus Dolycoris was classified under the tribe Dolycorini, whereas in the 42-tribe classification system, it was placed in the tribe Carpocorini, together with Carpocoris, Palomena, and Rubiconia.In this study, Carpocorini and Aeliini were found to be sister groups, this is congruent with Roca-Cusachs et al. [20].However, Lian et al. support the finding that Eysarcorini and Carpocorini are sister groups based on mitochondrial genomes [56].The external morphological features of Placosternum urus are similar to those of the tribe Pentatomini; however, shorter or ear-shaped scent gland openings are observed.The results of this study indicated that Placosternum urus forms a sister group relationship with the subfamily Phyllocephalinae, whereas Pentatoma semiannulata forms a sister group relationship with Eurydema dominulus (belonging to Strachiini).Graphosoma rubrolineata and Scotinophara lurida, both of which belong to the subfamily Podopinae, have a complex relationship with Pentatominae.Some Podopinae species were previously assigned to Scutelleridae or Pentatominae, and were even separately promoted to family level; however, there is no widely accepted view of Podopopinae based on a comprehensive analytical approach.Yang previously classified Graphosoma and Eysarcoris in the tribe Graphosomini of the subfamily Pentatominae [55].In addition, Zhao et al. and Li et al. also support a close relationship between the genus Graphosoma and species of Pentatominae, which is corroborated by the results of this study [57,58].In this study, the representative species of Podopinae mixed in the Pentatominae species, the result combined with pre-vious studies indicate that the classified position of Graphosoma rubrolineata could not be determined [26,56].Owing to the limitations of the sample size in this study, further investigation of the monophyly of the subfamily Pentatominae and the validity of taxonomic units within the subfamily Podopinae require additional mitochondrial sequence information.

Conclusions
Currently, there are few phylogenetic studies within Pentatomidae, and the use of molecular data is limited.This study determined the mitochondrial genome sequences of A. sibirica and A. fieberi.The sequence features, base composition, codon usage, and RNA secondary structures were found to be highly similar between the two species, with no significant differences observed within the genus.The comprehensive analysis of the Aelia mitochondrial genomes performed in this study contributes new information which can be used in molecular evolutionary studies and phylogenetic analyses.More mitogenomes or other molecular data are needed to conduct further exploration of the taxonomic status and phylogenetic history of Pentatomidae species.

Figure 2 .
Figure 2. The relative synonymous codon usage (RSCU) in the mitogenomes of A. sibirica and A. fieberi.

Figure 2 .
Figure 2. The relative synonymous codon usage (RSCU) in the mitogenomes of A. sibirica and A. fieberi.

Figure 3 .
Figure 3. Secondary structures of tRNAs of A. sibirica and A. fieberi (A.sibirica as the template).Figure 3. Secondary structures of tRNAs of A. sibirica and A. fieberi (A.sibirica as the template).

Figure 3 .
Figure 3. Secondary structures of tRNAs of A. sibirica and A. fieberi (A.sibirica as the template).Figure 3. Secondary structures of tRNAs of A. sibirica and A. fieberi (A.sibirica as the template).

Figure 4 .
Figure 4. Secondary structure of 12S rRNA gene of A. sibirica and A. fieberi (A.sibirica as the template).

Figure 4 .
Figure 4. Secondary structure of 12S rRNA gene of A. sibirica and A. fieberi (A.sibirica as the template).

Figure 5 .
Figure 5. Secondary structure of 16S rRNA gene of A. sibirica and A. fieberi (A.sibirica as the template).Figure 5. Secondary structure of 16S rRNA gene of A. sibirica and A. fieberi (A.sibirica as the template).

Figure 5 .
Figure 5. Secondary structure of 16S rRNA gene of A. sibirica and A. fieberi (A.sibirica as the template).Figure 5. Secondary structure of 16S rRNA gene of A. sibirica and A. fieberi (A.sibirica as the template).

Figure 6 .
Figure 6.ML and BI phylogenetic trees of Pentatomidae based on the protein-coding genes.Numbers above each node indicate Bayesian posterior probabilities values and ML bootstrap values.

Figure 6 .
Figure 6.ML and BI phylogenetic trees of Pentatomidae based on the protein-coding genes.Numbers above each node indicate Bayesian posterior probabilities values and ML bootstrap values.
Author Contributions: Conceptualization, W.Z. and Q.Z.; methodology, W.Z. and D.L.; software, D.L.; validation, D.L., Y.W. and W.Z.; formal analysis, S.F. and H.Z.; investigation, S.F. and Q.Z.; resources, W.Z.; data curation, D.L.; writing-original draft preparation, D.L. and W.Z.; writing-review and editing, H.Z. and S.F.; visualization, D.L. and Y.W.; supervision, W.Z. and Q.Z.; project administration, W.Z.; funding acquisition, W.Z. and H.Z. All authors have read and agreed to the published version of the manuscript.Funding: This research was funded by Fundamental Research Program of Shanxi Province, grant number 20210302124140 and 202103021224331, and funded by Scientific and Technological Innovation Programs of Higher Education Institutions in Shanxi, grant number 2023L299, and funded by the Fund for Shanxi "1331 Project" Improving the quality and efficiency Plan Project subtopic, grant number 2023WTS06.Institutional Review Board Statement: Not applicable.

Table 1 .
List of sequences used to reconstruct the phylogenetic tree.

Table 2 .
Organization of the mitochondrial genomes of A. sibirica and A. fieberi.

Table 2 .
Organization of the mitochondrial genomes of A. sibirica and A. fieberi.

Table 3 .
Nucleotide composition and skewness of mitochondrial genome of A. sibirica and A. fieberi.