Next Article in Journal
Spatio-Temporal Distribution Patterns and Determinant Factors of Wintering Hooded Cranes (Grus monacha) Population
Next Article in Special Issue
A New Enigmatic Teleost Fish from the Mid-Cretaceous of Lebanon
Previous Article in Journal
DNA Barcoding of Lepidoptera Species from the Maltese Islands: New and Additional Records, with an Insight into Endemic Diversity
Previous Article in Special Issue
Phylogenetic Reconstruction of the Rainforest Lineage Fontainea Heckel (Euphorbiaceae) Based on Chloroplast DNA Sequences and Reduced-Representation SNP Markers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of Two New Apodemus Mitogenomes (Rodentia: Muridae) and Mitochondrial Phylogeny of Muridae

1
School of Resources and Environmental Engineering, Anhui University, Hefei 230601, China
2
Anhui Province Key Laboratory of Wetland Ecosystem Protection and Restoration, Anhui University, Hefei 230601, China
*
Author to whom correspondence should be addressed.
Diversity 2022, 14(12), 1089; https://doi.org/10.3390/d14121089
Submission received: 16 October 2022 / Revised: 26 November 2022 / Accepted: 6 December 2022 / Published: 9 December 2022
(This article belongs to the Collection Feature Papers in Phylogeny and Evolution)

Abstract

:
Apodemus is the most common small rodent species in the Palearctic realm and an ideal species for biogeographical research and understanding environmental changes. Elucidating phylogenetic relationships will help us better understand species adaptation and genetic evolution. Due to its stable structure, maternal inheritance, and rapid evolution, the mitogenome has become a hot spot for taxonomic and evolutionary studies. In this research, we determined the mitochondrial genome of Apodemus agrarius ningpoensis and Apodemus draco draco and studied the phylogeny of Muridae using ML and BI trees based on all known complete mitogenomes. The mitochondrial genome of Apodemus agrarius ningpoensis was 16,262 bp, whereas that of Apodemus draco draco was 16,222 bp, and both encoded 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. Analysis of base composition showed a clear A-T preference. All tRNAs except tRNASer and tRNALys formed a typical trilobal structure. All protein-coding genes contained T- and TAA as stop codons. Phylogeny analysis revealed two main branches in the Muridae family. Apodemus agrarius ningpoensis formed sister species with Apodemus chevrieri, whereas Apodemus draco draco with Apodemus latronum. Our findings provide theoretical basis for future studies focusing on the mitogenome evolution of Apodemus.

1. Introduction

With the continuous improvement of molecular biology technology in recent years, the construction of molecular phylogenetic trees has become an indispensable component of systematics and taxonomy [1]. An increasing number of species differentiation analyses rely on molecular phylogenetic analyses. Therefore, establishing reliable phylogenetic relationships is a prerequisite for accurately distinguishing between species and revealing the mechanisms of species formation.
Many researchers have attempted to determine the evolutionary relationships of species using different methods. Early traditional taxonomies that relied on morphology severely underestimated species richness [2]. Later researchers began to use protein electrophoresis, DNA hybridization and other molecular biology tools to explore the relationship between biological systems [3]. In recent decades, mitochondrial gene sequences have become a hotspot of taxonomic and evolutionary research of different animal groups [4,5,6]. Mitochondria have their own unique genetic system and thus have become a widely used molecular marker in the study of phylogeny of species because of their matrilineal inheritance and high mutation rate [7,8,9,10,11,12,13]. Mitochondria provide energy for the daily cellular activities, and their genes change under the influence of the environment in which the species lives [14,15]. Physiological changes and evolutionary pressures on species are directly reflected in their mitochondrial genome sequence [14]. However, single molecular markers of mitochondrial gene sequences often have limited phylogenetic information [6,16,17,18]. Increasing use of two, three, or more mitochondrial markers has been found in phylogenetic studies [19,20]. Researchers can gain more information from the whole mitotic genome sequences owing to their multiple structural genomic characteristics in different groups [21], which are reflected in functionally dependent differences associated with structural differences or phylogenetically distant organisms [22]. The mitochondrial genome contains information on gene rearrangements, tRNAs structure, genetic system that are lacking in single gene fragments [16], hence improving the precision of the study of evolutionary relationships among related species [6,10,23,24]. The ongoing refinement of the composition and structural characteristics of the mitochondrial genome of species help us understand the evolutionary history of species. This will better promote the study of species taxonomy [16,25].
The genus Apodemus of the Muridae family species is abundant and widely distributed in the Palearctic [26]. Members of this genus are highly adaptable in almost all habitats in the mountains, and are the most common small rodent in the field [4,27,28]. However, these rodents are the most important hosts of zoonotic diseases and parasites [12,29] and have been associated with multiple hazards to humans, such as disease transmission and crop destruction [30]. Their wide distribution and unique adaptations make members of the genus Apodemus ideal candidates for understanding the relationship between environmental change, biogeography, and rodent species evolution [27]. Approximately 7.17–9.95 Mya ago, Apodemus underwent two stages of radiation evolution to produce the agrarius and draco subgroups [28]. Moreover, detailed descriptions of different genetic lineages are lacking, which is not conducive to identifying these species based on morphology alone [12].
Apodemus agrarius ningpoensis and Apodemus draco draco are the two common Apodemus subspecies in the northern margin of the Oriental and southern margins of the Palearctic in eastern China. However, their molecular characteristics are limited.
Phylogenetic relationships of Apodemus agrarius ningpoensis and Apodemus draco draco are difficult to determine owing to the paucity of mitochondrial genomic data. Typically, phylogenetic errors in the taxonomy of Apodemus have been attributed to recessive cytochrome b pseudogenes [31]. In this study, we sequenced the mitogenomes to explore their phylogenetic position and the composition and evolution of the mitochondrial genomes of Apodemus agrarius ningpoensis and Apodemus draco draco. Their genome sequences, nucleotide composition, genome size, codon usage, gene overlap, and tRNA secondary structure were analyzed. The mitochondrial genome data generated herein provide a reference for future research focusing on the evolution of Apodemus mitochondrial genome. In addition, our study will lay a foundation for the identification and classification of species, better promote the study of rodent taxonomy.

2. Materials and Methods

2.1. Experimental Samples

Specimens of Apodemus agrarius ningpoensis and Apodemus draco draco were collected from two counties in Dabie Mountains, Yue-xi and Huo-shan, in eastern China. Samples were collected using the clip-night method, with fresh peanut kernel as bait. Combining with habitat characteristics, 100 standard iron plate clips were placed in each sample plot. The clips were placed at dusk and collected next morning. The captured rodents were first weighed and numbered, and then the conventional data such as head–body length, tail length, ear length, and hind foot length were measured, and the sex was identified. After morphological identification, BLAST analysis in NCBI (http://ncbi.nlm.gov/nuccore?db=nucleotide, accessed on 3 February 2022) using DNA barcoding sequences verified their classification. Specimens were stored in 95% ethanol and placed in a −80 °C refrigerator.

2.2. DNA Extraction and Sequencing

The SPIN easy DNA Kit for Tissue and Bacteria (with lysing matrix) (MP Biomedicals, Hefei, China) was used to obtain the required DNA from muscle tissue. High-quality DNA samples were library constructed and sequenced by the Genepioneer Biotechnology (Nanjing, China) on the Illumina HiSeq2500 platform [32].

2.3. Assembly and Annotation

The raw sequences in FASTA format were extracted from the FASTQ format obtained by sequencing using Python. The mitochondrial genome was assembled using the SPAdes 3.15.3 software [33] by referring to the KY851939.1 sequence in NCBI. Using the MITOS web server (http://mitos.bioinf.unileipzig.de/index.py, accessed on 10 May 2022) [34] to annotate the assembled sequence. Annotation results were compared with those of related species, the final annotation results were obtained after manual correction [35]. tRNA structures were predicted using the Mitos web server and tRNAscan-SE 2.0 [36,37]. The mitochondrial genome of the two Apodemus obtained in this experiment was mapped using OGDRAW (https://chlorobox.mpimp-golm.mpg.de/OGDraw, accessed on 15 May 2022).

2.4. Sequence Analysis

The base composition of the mitogenome of the two Apodemus species was calculated by MEGA X [38]. AT-skew = (A − T)/(A + T) and GC skew = (G − C)/(G + C) [39] was used to calculate the relative numbers of AT skew and GC skew, respectively. The relative synonymous codon usage (RSCU) of the PCGs of the two Apodemus species was obtained using PhyloSuite v1.2.2 [40]. After using the PhyloSuite v1.2.2 to obtain the coding sequence (CDS) and remove the stop codon from the sequence based on the codon matching pattern of the MAFFT algorithm and G-INS-i (accurate) strategy, the specific value of nonsynonymous to synonymous substitutions Ka/Ks values of the two Apodemus subspecies was calculated using the DnaSP software [41].

2.5. Phylogenetic Analysis

There are 21 genera in the Muridae family that are distributed in China [42]. To determine the phylogenetic position of Apodemus agrarius ningpoensis and Apodemus draco draco, we constructed a phylogenetic tree using the whole mitogenome sequences from 63 species of 14 genera in the Muridae family provided by NCBI and this study. Ratufa bicolor and Pteromys volans of the Sciuridae family were selected for outgroups [43] (Table 1). Mitochondrial genome sequences of 67 species were aligned using the MAFFT plug-in in PhyloSuitev1.2.2 based on the codon pattern and G-INS-i strategy for comparison. The aligned data were trimmed using trimA [41]. Bayesian inference (BI) phylogenetic analysis was based on the GTR model using MrBayes v3.2.2 [44,45]. The Markov chain Monte Carlo (MCMC) for BI analysis was started with a random tree, run 4 Markov chains, 1,000,000 generations, taking a sample every 1000 generations with a burn-in of 25% trees. The RAxMLv.7.3.0 [46] was used to construct the maximum likelihood (ML) tree, with GTRGAMMA as the default model; 1000 bootstrap repeat tests were used to verify branch reliability. The ChiPlot (https://www.chiplot.online/, accessed on 12 June 2022) was used to visualize and beautify the phylogenetic trees.

3. Results

3.1. Features of Mitogenomes

The mitochondrial genome length of Apodemus agrarius ningpoensis was 16,262 bp, and Apodemus draco draco was 16,222 bp. Both mitochondrial genomes were composed of 37 genes (13 protein-coding genes (PCGs), 22 tRNAs, and 2 rRNAs) (Figure 1 and Figure 2). Of these 37 genes, the light chain encoded 9 genes ND6, tRNAGlu, tRNAPro, tRNAGln, tRNAAla, tRNAAsn, tRNACys, tRNATyr, and tRNASer, whereas the rest were detected in the heavy chain.
The overlapping and gapping regions in the mitogenomes of Apodemus agrarius ningpoensis and Apodemus draco draco, ranged from 1 bp to 7 bp. ATP8 and ATP6 (Table 2) had the longest overlap area. The completed mitochondrial genomes had a similar composition, sequences, and localization to other vertebrates. Interestingly, both mitochondrial genome sequences were biased toward A and T. The AT content of Apodemus agrarius ningpoensis was 63.74%, and that of Apodemus draco draco was 64.12%. The AT-skew value was greater than 0, whereas the GC skew value was less than 0, indicating that the base composition of Apodemus agrarius ningpoensis and Apodemus draco draco showed a strong A-bias and C-bias (Table 3 and Table 4).

3.2. Protein-Coding Genes and Codon Usage

The full length of the protein-coding genes of Apodemus agrarius ningpoensis and Apodemus draco draco was 11,406 bp and 11,202 bp, respectively. In the Apodemus agrarius ningpoensis mitogenome, the start codon of all PCGs was ATN, ATG was the most common start codon, identified in 9 PCGs (COX1, COX2, ATP8, ATP6, COX3, ND4L, ND4, ND6, and CYTB), whereas ATA was the start codon in the remaining PCGs. We observed that the ND1, ND2, COX1, COX3, ND4, and CYTB genes had the stop codon T-, whereas in the remaining seven PCGs had TAA. Interestingly, we detected that the third codon of PCGs tended to select base A. In the case of Apodemus draco draco, ND1, ND2, ND3, and ND5 started with an ATT codon, unlike the rest PCGs that used ATG as the start codon. Similar to Apodemus agrarius ningpoensis, ND1, ND2, COX1, COX3, ND4, and CYTB had the stop codon T-, whereas the rest had TAA.
The codon usage of PCGs was evaluated based on RSCU values (Table 5 and Table 6, Figure 3 and Figure 4), CUA of Leu, CGA of Arg, and UCA of Ser were the most frequent codons among all 13 PCGs in Apodemus agrarius ningpoensis. In contrast, CGA of Arg, CCA of Pro, and UCA of Ser were the most commonly used codons in Apodemus draco draco PCGs.
We calculated the Ka/Ks values for 13 PCGs from 65 Muridae species (Figure 5). All Ka/Ks values were smaller than 1, indicating genes subjected to purification selection. In particular, the mean and median Ka/Ks values were the highest for ATP8, suggesting that ATP8 had the highest evolution speed, the lowest for COX1, indicating the slowest evolution.

3.3. tRNAs and rRNAs

As expected, both Apodemus agrarius ningpoensis and Apodemus draco draco had 22 tRNA genes in their mitochondrial genomes, of which 9 were encoded on the negative chain, and the remaining 13 were encoded on the positive chain (Figure 6 and Figure 7). The 12S rRNA and 16S rRNA genes were located between the tRNAPhe and tRNALeu genes and segregated by the tRNAVal gene. In addition, tRNASer lacked the entire dihydrouracil (DHU) arm, tRNALys had a DHU arm lacking the D-loop to form a stable secondary structure, while the rest could form a cloverleaf structure. We further identified 36 mismatches in the Apodemus agrarius ningpoensis mitogenome, including 28 G-U, forming weak hydrogen bonds in tRNAs, as well as atypical A-A2, U-U1, A-G1, U-U2, C-U1, and A-C1. Similarly, there were 33 mismatches in the Apodemus draco draco mitogenome, including G-U26, U-U1, A-A2, A-G1, C-A2, and C-U1.

3.4. Phylogenetic Analysis

Phylogenetic trees were constructed based on ML and BI methods and yielded a consistent topology (Figure 8 and Figure 9). Most branches of the BI and ML trees received high node support. We identified two major clades in the Muridae family. One clade consisted of the genera Rhombomys, Meriones, Brachiones, and Hapalomys, whereas the others were composed of the genera Micromys, Chiropodomys, Apodemus, Mus, Maxomys, Leopoldamys, Niviventer, Berylmys, Bandicota, and Rattus. Interestingly, the genera Hapalomys, Micromys, Chiropodomys, Apodemus, Mus, Maxomys, Leopoldamys, Niviventer, Berylmys, Bandicota, and Rattus exhibited monophyly. However, the genus Meriones did not form a monophyletic group, and was first grouped with the genus Brachiones, followed by the genus Rhombomys g, and finally with the genus Meriones. The other branch of the relationship was (Apodemus + Mus) + (Maxomys + ((Leopoldamys + Niviventer) + (Berylmys + (Bandicota + Rattus))). Several members of the genus Apodemus formed an evident monophyletic group, which was divided into two small branches, among which Apodemus agrarius ningpoensis and Apodemus chevrieri were first clustered into a clade. Subsequently, it was grouped with Apodemus peninsulae, constituting a monophyletic group. The second group was composed of Apodemus flavicollis and Apodemus sylvaticus. The monophyletic lineage of the genus Apodemus was well supported by the sisters of the genus Mus; moreover, the first genus to be grouped with the genus Apodemus was the genus Mus.

4. Discussion

The length of the mitochondrial genome analyzed in this study was consistent with that of mammals, which is generally approximately 16,000 bp long [47]. The mitochondrial genome structure of Apodemus agrarius ningpoensis and Apodemus draco draco is a double-stranded closed loop, containing a noncoding control region sequence and encoding 37 genes. The distribution of mtDNA genes was not homogeneous. Most mitochondrial genes were encoded on the heavy chain, except for the ND6 gene and eight tRNA genes, which are encoded on the light chain. The apparent A-T preference in base composition was consistent with the basic characteristics of the Muridae mitochondrial genome, which is related to a base composition bias associated with the shedding of purine or pyrimidine bases from heavy chains during replication [48]. This high AT content and mutation rate of the mitochondrial genome have also been suggested to accelerate the evolutionary process in rodents [49,50]. Most animal mitochondrial genomes usually have obvious strand-specific biases, as expressed by AT and GC skews. The mitogenomes of the two species examined in this research showed positive AT, and negative GC skews, suggesting that the base composition of the whole mitogenomes of Apodemus agrarius ningpoensis and Apodemus draco draco is biased toward A and C. Synonymous codon frequency is generally biased owing to natural selection or mutational reasons, with different pressures on gene evolution. RSCU calculations showed that the protein-coding gene region of Apodemus agrarius ningpoensis and Apodemus draco draco had 3800 and 3799 codons, respectively. The codon CUA, CUU, UUA, and UCA determining Leu and Ser with RSCU values greater than 1, occupied a higher percentage in these mitogenomes, thus explaining the negative GC skew exhibited by the PCGs. The RSCU values indicate that there are significantly more bases in the third position of the codon for A and T than for G and C. This may be attributed to the AT preference for selective avoidance of termination codon formation and amino acid loss [51].
The calculation of Ka/Ks values revealed that ATP8 had the highest mean and median Ka/Ks values, implying that ATP8 had the fastest evolution among all PCGs, and is thus likely to be under relaxed selective pressure, accumulating more mutations [52]. In contrast, COX1 had the lowest mean and median Ka/Ks values, indicating it was under greater evolutionary pressure. All 13 PCGs had Ka/Ks values smaller than 1, indicating that sites underwent purification selection. It determines the development of the mitogenome [53].
Among the predicted 22 tRNA structures, we observed that tRNASer lacked the DHU arm and the DHU loop of tRNALys failed to form a typical cloverleaf structure. Most mammalian mitochondrial genomes lack DHU arms, which might be related to the compensatory mechanism between tRNA arms [22,54,55].
The presence of gene overlap and spacing in both species and gene overlap regions implied that Muridae species need to continuously improve DNA utilization from a supply-side perspective to meet evolutionary demands, which is extremely beneficial to the evolutionary process.
Both sequences had multiple G-U mismatches, which might be mainly attributed to the low free energy of the bond, which makes it stable and neutral [56]. This usually occurs in invertebrates and is rectified by post-transcriptional processing of RNA [57].
Currently, studies on the evolution and taxonomy of the Chinese Apodemus populations have concentrated on mitochondrial single gene sequences such as D-loop region and cytb. However, the mitochondrial genome can provide a better understanding of phylogenetic relationships [6,24]. The phylogenetic trees of the gene sequences of 63 rodent species in 14 genera of the Muridae family constructed using the BI and ML methods in this study were largely consistent, with bootstrap values of evolutionary branches >70% being considered indicative of strong relationships. This indicated a close phylogenetic relationship between species. The phylogenetic analysis of Muridae in this study was consistent with a system based on COB, IRBP, and GHR genes [43,58], albeit with the use of different datasets. We identified two major clades in the Muridae family. One clade consisted of the genera Rhombomys, Meriones, Brachiones, and Hapalomys, whereas the others were composed of the genera Micromys, Chiropodomys, Apodemus, Mus, Maxomys, Leopoldamys, Niviventer, Berylmys, Bandicota, and Rattus. Except for Rhombomys opimus, which originated from the same ancestor with the genera Brachiones and Meriones, the majority of Muridae species formed their own branches on a genus basis, which was consistent with previous studies [43]. Brachiones przewalskii has been reported to be nested in a diverse taxon corresponding to Rhombomys opimus and separated from Meriones [22]. In addition, Meriones tamariscinus formation predated the separation of Brachiones from Rhombomys [59], with Meriones tamariscinus being closely related to Rhombomys but more distantly related to any other Meriones species. The results of the study by Ito et al. supported our findings [60]. These morphological features might be associated with surrounding habitats. Ecological specialization allows for molecular and phenotypic evolution independent of each other [22,61]. Interestingly, the topology of the phylogenetic tree changed when other mitochondrial genomes were added to the study. The topology of the phylogenetic tree also changed when different species were selected for tree building.
Parallel analyses in the nuclear genome and mitogenome may probably influence phylogenetic studies in the future. The cytogenetic characteristics of the analyzed species are interesting [62]. Apodemus agrarius and Apodemus draco belong to the Apodemus group and Draco group, respectively, based on morphological characteristics. The karyotype of Apodemus agrarius had 3–5 biarmed metacentric autosome pairs (2n = 48, FNa = 54), whereas that of Apodemus draco had one small metacentric pair (2n = 48, FNa = 48). The karyotype of Apodemus chevrieri, a sister species of Apodemus agrarius, possesses four small metacentric pairs (2n = 48, FNa = 54), which was consistent with the standard karyotype of Apodemus agrarius. The sister species of Apodemus draco, Apodemus latronum, has a similar karyotype to Apodemus draco. However, they have different chromosome histotypes, and central heterochromatin is present in many distal center pairs. Two DNA phylogenetic species groups were distinct. Chromosomal rearrangements at the cytological level have contributed to the formation of post-mating reproductive isolation.

5. Conclusions

We reported and analyzed the mitochondrial genomes of Apodemus agrarius ningpoensis and Apodemus draco draco for the first time. The basic features of the mitogenomes of Apodemus agrarius ningpoensis and Apodemus draco draco, such as base content, gene arrangement, codon usage, and structures of tRNA, were highly conserved, also showing a relative abundance of AT content and a relative lack of GC content, which were consistent with the basic features of the Muridae mitochondrial genome. The new mitochondrial genome data contributed to shedding light on the phylogenetic relationship of the family Muridae and also revealed the phylogenetic position of Apodemus agrarius ningpoensis and Apodemus draco draco. The phylogenetic relationships among the various groups are consistent with related studies. It will provide genetic data for further evolutionary study of Apodemus.
In the future, developing primers using the mitochondrial genome sequence can help identify the most abundant loci in different rodent polymorphic types and combine the data with the construction of DNA barcode candidate markers to distinguish rodent populations. The data can be used for comparative analysis of mitochondrial genomes of different species, to understand disease information related to the source of the samples, and to provide a reference basis for the conservation and utilization of germplasm resources of the genus Apodemus. However, because of their wide distribution, representative Apodemus species exhibit particularly high levels of morphological and genetic variations. Expanding the sample size of the geographical population of Apodemus will help to better understand the genetic differentiation, dispersal, and phylogenetic relationships of Apodemus, as well as the influence of environmental mechanisms on its evolution.

Author Contributions

Conceptualization, L.Z.; data curation, D.W.; methodology, D.W.; software, D.W.; formal analysis, D.W.; visualization, D.W.; validation, D.W.; investigation, D.W., J.X., Q.X. and L.M.; funding acquisition, L.Z.; project administration, L.Z.; resource, L.Z.; supervision, L.Z.; writing—original draft preparation, D.W.; writing—review and editing, D.W. and L.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by grants for the National Science and Technology Fundamental Resources Investigation Program of China (2019FY101804).

Institutional Review Board Statement

The animal study protocol was approved by the Institutional Animal Care and Use Committee of Anhui Zoological Society (2022102001) on 20 October 2022.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tokumasa, H. An introduction to molecular phylogenetic analysis. Rev. Agric. Sci. 2016, 4, 36–45. [Google Scholar] [CrossRef]
  2. Kryštufek, B.; Lužnik, M.; Buzan, E.V. Mitochondrial cytochrome b sequences resolve the taxonomy of field mice (Apodemus) in the western Balkan refugium. Acta Theriol. 2012, 57, 1–7. [Google Scholar] [CrossRef]
  3. Atopkin, D.M.; Bogdanov, A.S.; Chelomina, G.N. Genetic variation and differentiation in striped field mouse Apodemus agrarius inferred from RAPD-PCR analysis. Genetika 2007, 43, 804–817. [Google Scholar] [CrossRef]
  4. Yue, H.; Fan, Z.; Liu, S.; Liu, Y.; Song, Z.; Zhang, X. A mitogenome of the chevrier’s field mouse (Apodemus Chevrieri) and genetic variations inferred from the cytochrome b gene. DNA Cell Biol. 2012, 31, 460–469. [Google Scholar] [CrossRef] [PubMed]
  5. Jing, M.; Yu, H.; Wu, S.; Wang, W.; Zheng, X. Phylogenetic relationships in genus Niviventer (Rodentia: Muridae) in China inferred from complete mitochondrial cytochrome b gene. Mol. Phylogenet. Evol. 2007, 44, 521–529. [Google Scholar] [CrossRef] [PubMed]
  6. Liu, S.; Li, F.; Zhou, J.; Lv, J.; Tan, Z.; Zhang, Y.; Ge, X. The phylogeny of the Anderson’s white-bellied rat (Niviventer andersoni) based on complete mitochondrial genomes. Ecol. Evol. 2022, 12, e8663. [Google Scholar] [CrossRef]
  7. Galtier, N.; Nabholz, B.; Glemin, S.; Hurst, G.D.D. Mitochondrial DNA as a marker of molecular diversity: A reappraisal. Mol. Ecol. 2009, 18, 4541–4550. [Google Scholar] [CrossRef]
  8. An, J.; Fan, C.; Fu, Z.; Zhang, H.; Yang, P. Analyses of the complete mitochondrial genome of Paraconiothyrium sp. and gene rearrangement diversity in the Pleosporales. Diversity 2022, 14, 601. [Google Scholar] [CrossRef]
  9. Wang, X.; Pei, J.; Bao, P.; Cao, M.; Guo, S.; Song, R.; Song, W.; Liang, C.; Yan, P.; Guo, X. Mitogenomic diversity and phylogeny analysis of yak (Bos grunniens). BMC Genom. 2021, 22, 325. [Google Scholar] [CrossRef]
  10. Sullivan, K.A.; Platt, R.N.; Bradley, R.D.; Ray, D.A. Whole mitochondrial genomes provide increased resolution and indicate paraphyly in deer mice. BMC Zool. 2017, 2, 11. [Google Scholar] [CrossRef]
  11. Kim, H.R.; Park, Y.C. The complete mitochondrial genome of the striped field mouse, Apodemus agrarius (Rodentia, Murinae) from Korea. Mitochondrial DNA 2012, 23, 145–147. [Google Scholar] [CrossRef] [PubMed]
  12. Ge, D.; Feijó, A.; Cheng, J.; Lu, L.; Liu, R.; Abramov, A.V.; Xia, L.; Wen, Z.; Zhang, W.; Shi, L.; et al. Evolutionary history of field mice (Murinae: Apodemus), with emphasis on morphological variation among species in China and description of a new species. Zool. J. Linn. Soc. Lond. 2019, 187, 518–534. [Google Scholar] [CrossRef]
  13. Moshtaghi, S.; Darvish, J.; Mirshamsi, O.; Mahmoudi, A. Cryptic species diversity in the genus Allactaga (Rodentia: Dipodidae) at the edge of its distribution range. Folia Zool. 2016, 65, 142–147. [Google Scholar] [CrossRef]
  14. Plese, B.; Kenny, N.J.; Rossi, M.E.; Cárdenas, P.; Schuster, A.; Taboada, S.; Koutsouveli, V.; Koutsouveli, V. Mitochondrial evolution in the Demospongiae (Porifera): Phylogeny, divergence time, and genome biology. Mol. Phylogenet. Evol. 2020, 155, 107011. [Google Scholar] [CrossRef] [PubMed]
  15. Dowling, D.; Friberg, U.; Lindell, J. Evolutionary implications of non-neutral mitochondrial genetic variation. Trends Ecol. Evol. 2008, 23, 546–554. [Google Scholar] [CrossRef]
  16. Ye, W.; Wang, J.; Zhao, X.; Liu, H.; Zhu, S. Mitochondrial genomes of two Lycosa spiders (Araneae, Lycosidae): Genome description and phylogenetic implications. Diversity 2022, 14, 538. [Google Scholar] [CrossRef]
  17. Wu, B.; Zhou, P.; He, J.; Cui, B.; Chen, F.; Sun, L. The taxonomic research of rodents with related target genes. Chin. J. Vector Biol. Control 2016, 27, 520–524. [Google Scholar] [CrossRef]
  18. Steppan, S.J.; Schenk, J.J. Muroid rodent phylogenetics: 900-species tree reveals increasing diversification rates. PLoS ONE 2017, 12, e183070. [Google Scholar] [CrossRef] [Green Version]
  19. Bellinvia, E. A phylogenetic study of the genus Apodemus by sequencing the mitochondrial DNA control region. J. Zool. Syst. Evol. Res. 2004, 42, 289–297. [Google Scholar] [CrossRef]
  20. Yalkovskaya, L.; Sibiryakov, P.; Borodin, A. Phylogeography of the striped field mouse (Apodemus agrarius Pallas, 1771) in light of new data from central part of Northern Eurasia. PLoS ONE 2022, 17, e276466. [Google Scholar] [CrossRef]
  21. Moritz, C.; Dowling, T.E.; Brown, W.M. Evolution of animal mitochondrial DNA: Relevance for population biology and systematics. Annu. Rev. Ecol. Syst. 1987, 18, 269–292. [Google Scholar] [CrossRef]
  22. Ding, L.; Luo, G.; Zhou, Q.; Sun, Y.; Liao, J. Comparative mitogenome analysis of gerbils and the mitogenome phylogeny of Gerbillinae (Rodentia: Muridae). Biochem. Genet. 2022, 60, 2226–2249. [Google Scholar] [CrossRef] [PubMed]
  23. Yue, H.; Yan, C.; Tu, F.; Yang, C.; Ma, W.; Fan, Z.; Song, Z.; Owens, J.; Liu, S.; Zhang, X. Two novel mitogenomes of Dipodidae species and phylogeny of Rodentia inferred from the complete mitogenomes. Biochem. Syst. Ecol. 2015, 60, 123–130. [Google Scholar] [CrossRef]
  24. İbiş, O. Whole mitochondrial genome sequence and phylogenetic relationships of Williams’s jerboa (Scarturus williamsi) from Turk. PeerJ 2020, 8, e9569. [Google Scholar] [CrossRef]
  25. Chen, L.; Lin, Y.; Xiao, Q.; Lin, Y.; Du, Y.; Lin, C.; Ward-Fear, G.; Hu, C.; Qu, Y.; Li, H. Characterization of the complete mitochondrial genome of the many-lined sun skink (Eutropis multifasciata) and comparison with other Scincomorpha species. Genomics 2021, 113, 2526–2536. [Google Scholar] [CrossRef]
  26. Oh, D.; Kim, T.; Chang, M.; Han, S.; Oh, H.; Kim, S. Mitochondrial genome sequences of the striped field mice Apodemus agrarius coreae and Apodemus agrarius chejuensis. Mitochondrial DNA 2011, 22, 165–167. [Google Scholar] [CrossRef]
  27. Wei, H.; Jia, Q.; Li, F.; Liu, Y.; Chen, S.; Yong, B. Sequencing and analysis of complete mitochondrial genome of Apodemus draco (Rodentia: Arvicolinae). Mitochondrial DNA Part A 2016, 27, 2740–2741. [Google Scholar] [CrossRef]
  28. Liu, X.; Wei, F.; Li, M.; Jiang, X.; Feng, Z.; Hu, J. Molecular phylogeny and taxonomy of wood mice (genus Apodemus Kaup, 1829) based on complete mtDNA cytochrome b sequences, with emphasis on Chinese species. Mol. Phylogenet. Evol. 2004, 33, 1–15. [Google Scholar] [CrossRef]
  29. Liu, S.; He, K.; Chen, S.; Jin, W.; Murphy, R.W.; Tang, M.; Liao, R.; Li, F. How many species of Apodemus and Rattus occur in China? A survey based on mitochondrial cyt b and morphological analyses. Zool. Res. 2018, 39, 309–320. [Google Scholar] [CrossRef]
  30. Xiong, Y.; You, F.; Chen, X.; Chen, Y.; Wen, Y.; Chen, Q. Detection and phylogenetic analysis of porcine bocaviruses carried by murine rodents and house shrews in China. Transbound. Emerg. Dis. 2018, 66, 259–267. [Google Scholar] [CrossRef]
  31. Dubey, S.; Michaux, J.; Brünner, H.; Hutterer, R.; Vogel, P. False phylogenies on wood mice due to cryptic cytochrome-b pseudogene. Mol. Phylogenet. Evol. 2009, 50, 633–641. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef] [Green Version]
  34. 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. Phylogenet. Evol. 2013, 69, 313–319. [Google Scholar] [CrossRef] [PubMed]
  35. Bernt, M.; Braband, A.; Schierwater, B.; Stadler, P.F. Genetic aspects of mitochondrial genome evolution. Mol. Phylogenet. Evol. 2013, 69, 328–338. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Lowe, T.M.; Chan, P.P. tRNAscan-SE On-line: Integrating search and context for analysis of transfer RNA genes. Nucleic Acids Res. 2016, 44, W54–W57. [Google Scholar] [CrossRef]
  37. Chan, P.P.; Lowe, T.M. tRNAscan-SE: Searching for tRNA genes in genomic sequences. Methods Mol. Biol. 2019, 1962, 1–14. [Google Scholar] [CrossRef]
  38. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  39. Perna, N.T.; Kocher, T.D. Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes. J. Mol. Evol. 1995, 3, 353–358. [Google Scholar] [CrossRef]
  40. 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. 2019, 20, 348–355. [Google Scholar] [CrossRef]
  41. Librado, P.; Rozas, J. Dnasp v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009, 25, 1451–1452. [Google Scholar] [CrossRef] [Green Version]
  42. Wei, F.; Yang, Q.; Wu, Y.; Jiang, X.; Liu, S.; Li, B.; Yang, G.; Li, M.; Zhou, J.; Li, S.; et al. Catalogue of mammals in China (2021). Acta Theriol. Sin. 2021, 41, 487–501. [Google Scholar] [CrossRef]
  43. Song, S.; Yong, H.; Lim, P.; Tan, J. Mitochondrial genome of Rattus tiomanicus (Rodentia: Muridae) and molecular phylogeny of Murinae. Sains Malays. 2021, 50, 953–965. [Google Scholar] [CrossRef]
  44. Ronquist, F.; Huelsenbeck, J.P. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 2003, 19, 1572–1574. [Google Scholar] [CrossRef] [Green Version]
  45. 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] [Green Version]
  46. Stamatakis, A. Raxml version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef] [Green Version]
  47. Balakirev, A.E.; Abramov, A.V.; Rozhnov, V.V. Revision of the genus Leopoldamys (Rodentia, Muridae) as inferred from morphological and molecular data, with a special emphasis on the species composition in continental Indochina. Zootaxa 2013, 3640, 521–549. [Google Scholar] [CrossRef] [Green Version]
  48. Reyes, A.; Gissi, C.; Pesole, G.; Saccone, C. Asymmetrical directional mutation pressure in the mitochondrial genome of mammals. Mol. Biol. Evol. 1998, 15, 957–966. [Google Scholar] [CrossRef] [Green Version]
  49. Triant, D.A.; DeWoody, J.A. Accelerated molecular evolution in Microtus (Rodentia) as assessed via complete mitochondrial genome sequences. Genetica 2006, 128, 95–108. [Google Scholar] [CrossRef]
  50. Yakovchuk, P.; Protozanova, E.; Frank-Kamenetski, M.D. Base-stacking and base-pairing contributions into thermal stability of the DNA double helix. Nucleic Acids Res. 2006, 34, 564–574. [Google Scholar] [CrossRef]
  51. Charneski, C.A.; Honti, F.; Bryant, J.M.; Hurst, L.D.; Feil, E.J. Atypical AT skew in firmicute genomes results from selection and not from mutation. PLoS Genet. 2011, 7, e1002283. [Google Scholar] [CrossRef] [Green Version]
  52. Hurst, L.D. The Ka/Ks ratio: Diagnosing the form of sequence evolution. Trends Genet. 2002, 18, 486. [Google Scholar] [CrossRef]
  53. Wei, H.; Li, F.; Wang, X.; Wang, Q.; Chen, G.; Zong, H.; Chen, S. The characterization of complete mitochondrial genome and phylogenetic relationship within Rattus genus (Rodentia: Muridae). Biochem. Syst. Ecol. 2017, 71, 179–186. [Google Scholar] [CrossRef]
  54. Jiang, X.; Gao, J.; Ni, L.; Hu, J.; Li, K.; Sun, F.; Xie, J.; Bo, X.; Gao, C.; Xiao, J.; et al. The complete mitochondrial genome of Microtus fortis calamorum (Arvicolinae, Rodentia) and its phylogenetic analysis. Gene 2012, 498, 288–295. [Google Scholar] [CrossRef]
  55. Wolstenholme, D.R. Animal mitochondria DNA: Structure and evolution. Int. Rev. Cytol. 1992, 141, 173–216. [Google Scholar] [CrossRef]
  56. Liu, H.; Chen, S.; Chen, Q.; Pu, D.; Chen, Z.; Liu, Y.; Liu, X. The first mitochondrial genomes of the family Haplodiplatyidae (Insecta: Dermaptera) reveal intraspecific variation and extensive gene rearrangement. Biology 2022, 11, 807. [Google Scholar] [CrossRef]
  57. Lavrov, D.V.; Brown, W.M.; Boore, J.L. A novel type of RNA editing occurs in the mitochondrial tRNAs of the centipede Lithobius forficatus. Proc. Natl. Acad. Sci. USA 2000, 97, 13738–13742. [Google Scholar] [CrossRef] [Green Version]
  58. Lecompte, E.; Aplin, K.; Denys, C.; Catzeflis, F.; Chades, M.; Chevret, P. Phylogeny and biogeography of African Murinae based on mitochondrial and nuclear gene sequences, with a new tribal classification of the subfamily. BMC Evol. Biol. 2008, 1, 199. [Google Scholar] [CrossRef] [Green Version]
  59. Alhajeri, B.H.; Hunt, O.J.; Steppan, S.J. Molecular systematics of gerbils and deomyines (Rodentia: Gerbillinae, Deomyinae) and a test of desert adaptation in the tympanic bulla. J. Zool. Syst. Evol. Res. 2015, 53, 312–330. [Google Scholar] [CrossRef]
  60. Ito, M.; Jiang, W.; Sato, J.J.; Zhen, Q.; Jiao, W.; Goto, K.; Sato, H.; Ishiwata, K.; Oku, Y.; Chai, J.; et al. Molecular phylogeny of the subfamily Gerbillinae (Muridae, Rodentia) with emphasis on species living in the Xinjiang-Uygur Autonomous Region of China and based on the mitochondrial cytochrome b and cytochrome c oxidase subunit II genes. Zool. Sci. 2010, 27, 269–278. [Google Scholar] [CrossRef]
  61. Renaud, S.; Chevret, P.; Michaux, J. Morphological vs. molecular evolution: Ecology and phylogeny both shape the mandible of rodents. Zool. Scr. 2007, 36, 525–535. [Google Scholar] [CrossRef]
  62. Motokawa, M.; Wu, Y.; Harada, M.; Shintaku, Y.; Jiang, X.L.; Li, Y.C. Karyotypes of field mice of the genus Apodemus (Mammalia: Rodentia) from China. Zool. Res. 2018, 39, 348–355. [Google Scholar] [CrossRef]
Figure 1. Annotated mitogenome of Apodemus agrarius ningpoensis.
Figure 1. Annotated mitogenome of Apodemus agrarius ningpoensis.
Diversity 14 01089 g001
Figure 2. Annotated mitogenome of Apodemus draco draco.
Figure 2. Annotated mitogenome of Apodemus draco draco.
Diversity 14 01089 g002
Figure 3. RSCU values of Apodemus agrarius ningpoensis.
Figure 3. RSCU values of Apodemus agrarius ningpoensis.
Diversity 14 01089 g003
Figure 4. The RSCU values of Apodemus draco draco.
Figure 4. The RSCU values of Apodemus draco draco.
Diversity 14 01089 g004
Figure 5. The Ka/Ks values for the 13 PCGs of 65 Muridae mitogenomes.
Figure 5. The Ka/Ks values for the 13 PCGs of 65 Muridae mitogenomes.
Diversity 14 01089 g005
Figure 6. Structures of tRNA of Apodemus agrarius ningpoensis.
Figure 6. Structures of tRNA of Apodemus agrarius ningpoensis.
Diversity 14 01089 g006
Figure 7. Structures of tRNA of Apodemus draco draco.
Figure 7. Structures of tRNA of Apodemus draco draco.
Diversity 14 01089 g007
Figure 8. Phylogenetic tree inferred from mitochondrial genome sequences using ML. Bootstrap support values are indicated on branches. Different genera are distinguished by different colors.
Figure 8. Phylogenetic tree inferred from mitochondrial genome sequences using ML. Bootstrap support values are indicated on branches. Different genera are distinguished by different colors.
Diversity 14 01089 g008
Figure 9. Phylogenetic tree inferred from mitochondrial genome sequences using BI. Bootstrap support values are indicated on branches. Different genera are distinguished by different colors.
Figure 9. Phylogenetic tree inferred from mitochondrial genome sequences using BI. Bootstrap support values are indicated on branches. Different genera are distinguished by different colors.
Diversity 14 01089 g009
Table 1. Mitogenomes used in the analyses.
Table 1. Mitogenomes used in the analyses.
GenusSpeciesAccession No.
ApodemusApodemus chevrieriNC_017599.1
ApodemusApodemus latronumNC_019585.1
ApodemusApodemus flavicollisMN122902.1
ApodemusApodemus peninsulaeNC_016060.1
ApodemusApodemus dracoNC_019584.1
ApodemusApodemus chejuensisNC_016662.1
ApodemusApodemus agrariusNC_016428.1
ApodemusApodemus sylvaticusNC_049122.1
RattusRattus andamanensisNC_046686.1
RattusRattus tiomanicusNC_029888.1
RattusRattus hoogerwerfiNC_049040.1
RattusRattus baluensisNC_035621.1
RattusRattus sordidusNC_014871.1
RattusRattus lutreolusNC_014858.1
RattusRattus villosissimusNC_014864.1
RattusRattus tunneyiNC_014861.1
RattusRattus fuscipesNC_014867.1
RattusRattus leucopusNC_014855.1
RattusRattus niobeNC_023347.1
RattusRattus rattusNC_012374.1
RattusRattus praetorNC_012461.1
RattusRattus norvegicusNC_001665.2
RattusRattus tanezumiNC_011638.1
RattusRattus exulansNC_012389.1
RattusRattus nitidusNC_040919.1
HapalomysHapalomys delacouriMZ159976.1
NiviventerNiviventer cremoriventerNC_035822.1
NiviventerNiviventer lotipesNC_065402.1
NiviventerNiviventer sacerMZ935252.1
NiviventerNiviventer excelsiorNC_019617.1
NiviventerNiviventer confucianusNC_023960.1
NiviventerNiviventer fulvescensNC_028715.1
NiviventerNiviventer andersoniNC_060500.1
LeopoldamysLeopoldamys sabanusNC_035819.1
LeopoldamysLeopoldamys edwardsiNC_025670.1
MusMus musculus musculusNC_010339.1
MusMus cervicolorNC_025269.1
MusMus cookiiNC_025270.1
MusMus fragilicaudaNC_025287.1
MusMus terricolorNC_010650.1
MusMus musculus helgolandicusKP877610.1
MusMus musculus castaneusNC_012387.1
MusMus musculus domesticusNC_006914.1
MusMus musculus molossinusNC_006915.1
MusMus spretusNC_025952.1
MusMus famulusNC_030342.1
MusMus pahariNC_036680.1
MusMus caroliNC_025268.1
BerylmysBerylmys berdmoreiNC_036730.1
BrachionesBrachiones przewalskiiKT834972.1
MaxomysMaxomys ochraceiventerNC_056988.1
MaxomysMaxomys suriferNC_036732.1
MaxomysMaxomys whiteheadiNC_049119.1
BandicotaBandicota bengalensisNC_057104.1
BandicotaBandicota indicaKT029807.1
MicromysMicromys erythrotisNC_060316.1
MicromysMicromys minutusNC_027932.1
ChiropodomysChiropodomys gliroidesNC_049121.1
MerionesMeriones unguiculatusKF425526.1
MerionesMeriones tamariscinusNC_034314.1
MerionesMeriones meridianusNC_027684.1
MerionesMeriones libycusNC_027683.1
RhombomysRhombomys opimusMK359635.1
PteromysPteromys volansNC_019612.1
RatufaRatufa bicolorNC_023780.1
Table 2. Mitochondrial composition of Apodemus agrarius ningpoensis (left) and Apodemus draco draco (right).
Table 2. Mitochondrial composition of Apodemus agrarius ningpoensis (left) and Apodemus draco draco (right).
GeneDirectionPosition (bp)Size (bp)Start CodonStop CodonIntergenic Nucleotides
tRNAPhe+1–67\1–6867\68 0\0
12S rRNA +68–1022\69–1023955\955 0\0
tRNAVal+1023–1091\1024–109169\68 0\0
16S rRNA+1092–2665\1092–26601547\1569 0\0
tRNALeu+2664–2738\2659–273375\75 −2\−2
ND1+2736–3693\2731–3688958\958ATA\ATAT-\T-−3\−3
tRNAIle+3694–3762\3689–375769\69 0\0
tRNAGln3760–3830\3755–382571\71 −3\−3
tRNAMet+3835–3903\3830–389869\69 4\4
ND2+3094–4939\3899–49341036\1036ATA\ATAT-\T-0\0
tRNATrp+4940–5007\4935–500268\68 0\0
tRNAAla5009–5077\5004–507269\69 1\1
tRNAAsn5079–5149\5074–514471\71 1\1
tRNACys5181–5248\5176–524168\66 31\31
tRNATyr5249–5315\5242–530867\67 0\0
COX1+5317–6859\5310–68521543\1543ATG\ATGT-\T-1\1
tRNASer6859–6927\6852–692069\69 −1\−1
tRNAAsp+6931–6998\6924–699168\68 3\3
COX2+7000–7683\6993–7676684\684ATG\ATGTAA\TAA1\1
tRNALys+7687–7750\7680–774364\64 3\3
ATP8+7752–7955\7745–7948204\204ATG\ATGTAA\TAA1\1
ATP6+7913–8593\7906–8586681\681ATG\ATGTAA\TAA−43\−43
COX3+8593–9376\8586–9369784\784ATG\ATGT-\T-−1\−1
tRNAGly+9377–9444\9370–943668\67 0\0
ND3+9445–9792\9437–9784348\348ATA\ATTTAA\TAA0\0
tRNAArg+9794–9861\9786–985368\68 1\1
ND4L+9864–10160\9856–10152297\297ATG\ATGTAA\TAA2\2
ND4+10154–11531\10146–115231378\1378ATG\ATGT-\T-−7\−7
tRNAHis+11532–11599\11524–1159068\67 0\0
tRNASer+11600–11658\11591–1164959\59 0\0
tRNALeu+11658–11728\11649–1171971\71 −1\−1
ND5+11720–13549\11711–135401830\1830ATA\ATATAA\TAA−9\−9
ND613536–14054\13527–14042519\516ATG\ATGTAA\TAA−14\−14
tRNAGlu14055–14123\14043–1411169\69 0\0
CYTB+14127–15270\14115–152581144\1144ATG\ATGT-\T-3\3
tRNAThr+15271–15339\15259–1532669\68 0\0
tRNAPro15340–15406\15327–1539367\67 0\0
D-loop+15407–16262\15394–16222856\829 0\0
Table 3. Base content of the mitochondrial genome of Apodemus agrarius ningpoensis.
Table 3. Base content of the mitochondrial genome of Apodemus agrarius ningpoensis.
Apodemus agrarius
ningpoensis
Size (bp)A%T%G%C% A + T%G + C%AT SkewGC Skew
Mitogenome16,26234.0129.7412.4723.7963.7436.260.067−0.312
PCGs11,40632.0531.4512.124.463.536.50.01−0.337
tRNAs150334.1330.8118.4316.6364.9435.060.0510.051
rRNAs252937.9226.4916.7718.8264.4135.590.177−0.058
D-loop85632.7131.6611.9223.7164.3735.630.016−0.331
Table 4. Base content of the mitochondrial genome of Apodemus draco draco.
Table 4. Base content of the mitochondrial genome of Apodemus draco draco.
Apodemus draco dracoSize (bp)A%T%G%C%A + T%G + C%AT SkewGC Skew
Mitogenome16,22235.1528.9712.0823.864.1235.880.096−0.327
PCGs11,40333.3330.6411.6824.3463.9736.030.042−0.352
tRNAs149834.3131.5118.2915.8965.8234.180.0430.07
rRNAs252439.2625.7916.3218.6265.0634.940.207−0.066
D-loop82931.9729.6713.1525.2161.6438.360.037−0.314
Table 5. The codon usages of Apodemus agrarius ningpoensis.
Table 5. The codon usages of Apodemus agrarius ningpoensis.
CodonNo.RSCUCodonNo.RSCUCodonNo.RSCU
UAA(*)74AAA(K)981.98CGU(R)100.62
GCA(A)931.62AAG(K)10.02AGC(S)360.72
GCC(A)661.15CUA(L)2582.58AGU(S)200.4
GCG(A)40.07CUC(L)800.8UCA(S)1382.77
GCU(A)671.17CUG(L)190.19UCC(S)541.08
UGC(C)201.38CUU(L)1091.09UCG(S)20.04
UGU(C)90.62UUA(L)1231.23UCU(S)490.98
GAC(D)371.04UUG(L)100.1ACA(T)1612.13
GAU(D)340.96AUA(M)1981.7ACC(T)600.79
GAA(E)781.68AUG(M)350.3ACG(T)40.05
GAG(E)150.32AAC(N)921.1ACU(T)781.03
UUC(F)970.79AAU(N)760.9GUA(V)852.01
UUU(F)1481.21CCA(P)1252.51GUC(V)250.59
GGA(G)1011.92CCC(P)280.56GUG(V)90.21
GGC(G)400.76CCG(P)30.06GUU(V)501.18
GGG(G)260.5CCU(P)430.86UGA(W)931.81
GGU(G)430.82CAA(Q)771.88UGG(W)100.19
CAC(H)531.1CAG(Q)50.12UAC(Y)600.92
CAU(H)430.9CGA(R)452.81UAU(Y)701.08
AUC(I)1460.79CGC(R)70.44
AUU(I)2251.21CGG(R)20.12
Asterisk “*” indicates that it is a stop codon and has no corresponding amino acid.
Table 6. The codon usages of Apodemus draco draco.
Table 6. The codon usages of Apodemus draco draco.
CodonNo.RSCUCodonNo.RSCUCodonNo.RSCU
UAA(*)74AUU(I)2161.16CGU(R)70.44
GCA(A)891.6AAG(K)30.06CGA(R)462.88
GCC(A)761.36AAA(K)991.94AGU(S)100.2
GCG(A)10.02CUC(L)660.67UCA(S)1633.21
GCU(A)571.02CUG(L)80.08UCC(S)370.73
UGC(C)161.1CUU(L)1041.06UCG(S)50.1
UGU(C)130.9UUA(L)1341.37UCU(S)460.9
GAC(D)481.33UUG(L)130.13AGC(S)440.87
GAU(D)240.67CUA(L)2642.69ACC(T)710.92
GAA(E)841.79AUG(M)280.23ACG(T)40.25
GAG(E)100.21AUA(M)2111.77ACU(T)610.79
UUC(F)1160.95AAU(N)720.88ACA(T)1732.24
GGA(G)1182.26CCC(P)160.32GUG(V)100.24
GGC(G)270.52CCU(P)360.73GUU(V)380.91
GGG(G)190.36CCA(P)1452.94GUA(V)1012.42
GGU(G)450.86CAG(Q)10.02UGA(W)981.9
CAC(H)601.21CAA(Q)811.98UGG(W)50.1
CAU(H)390.79CGC(R)100.62UAC(Y)590.91
AUC(I)1570.84CGG(R)10.06UAU(Y)701.09
Asterisk “*” indicates that it is a stop codon and has no corresponding amino acid.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wu, D.; Zhou, L.; Xue, J.; Xia, Q.; Meng, L. Characterization of Two New Apodemus Mitogenomes (Rodentia: Muridae) and Mitochondrial Phylogeny of Muridae. Diversity 2022, 14, 1089. https://doi.org/10.3390/d14121089

AMA Style

Wu D, Zhou L, Xue J, Xia Q, Meng L. Characterization of Two New Apodemus Mitogenomes (Rodentia: Muridae) and Mitochondrial Phylogeny of Muridae. Diversity. 2022; 14(12):1089. https://doi.org/10.3390/d14121089

Chicago/Turabian Style

Wu, Di, Lizhi Zhou, Jiezhong Xue, Qiliang Xia, and Lei Meng. 2022. "Characterization of Two New Apodemus Mitogenomes (Rodentia: Muridae) and Mitochondrial Phylogeny of Muridae" Diversity 14, no. 12: 1089. https://doi.org/10.3390/d14121089

APA Style

Wu, D., Zhou, L., Xue, J., Xia, Q., & Meng, L. (2022). Characterization of Two New Apodemus Mitogenomes (Rodentia: Muridae) and Mitochondrial Phylogeny of Muridae. Diversity, 14(12), 1089. https://doi.org/10.3390/d14121089

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