Posted on Authorea 7 Jul 2021 — The copyright holder is the author/funder. All rights reserved. No reuse without permission. — https://doi.org/10.22541/au.162566526.66282591/v1 — This a preprint and has not been peer reviewed. Data may be preliminary. . N. with infectious together, common clustered the most that initially showed the the 1981), gene comprised among variety (cytb) (Musser, and a also b together studies inhabit cytochrome are previous mitochondrial They They on to the [6]. based valleys. tails According divisions: on dry primary tail-tip to [7]. two All the forests humans and damp in craniums [5]. from flat agents Islands slender, ranging Sunda long, habitats Great the of by the to murid China other and Himalayas during China the the southward Southeast expand of to regions potential Niviventer altitude its low suggesting [3,4]. Provinces, the colder Guizhou to turned and extended Chongqing climate species in this Quaternary late that the showed records Fossil [2]. of species congeneric Rodentia. Rat White-bellied Anderson’s Introduction 1. genus the in evolution mitogenomic of understanding our improve and that clearer tree of Niviventer gene reconsideration Our further (cytb) PCGs topologies. to b well-resolved PCGs, helpful cytochrome identical, be 13 one nearly all from produce of hypothesis D-loop, the sequences the reject RNA nucleotide for ribosomal the except two of that sequence (PCGs), analyses support mitogenome genes results Phylogenetic entire protein-coding the 13 (D-Loop). and genes, region ND6 RNA control minus transfer non-coding 22 one of and consists genes, genome of mitochondrial genome bp mitochondrial complete 16,291 first The the sequenced have we study, present the Rat, White-bellied Anderson’s Abstract 2021 7, July 4 3 2 1 Ge Xingyi Liu Shujing history evolutionary its into insight ( Rat White-bellied andersoni Anderson’s the of phylogeny The hns cdm fIseto n urnieIsiueo nmlQuarantine of Institute Quarantine and Prevention Inspection and of Control Academy Diseases Chinese Endemic of Institute Yunnan University Dali University Hunan . .andersoni N. otis1 eonzdseis ihaohr6 r eonzda yoys curn from occurring synonyms, as recognized are 65 another with species, recognized 17 contains 1 1 Dvso.Hwvr igegn eune r oeie iie nterpoiino useful of provision their in limited sometimes are sequences gene single However, -Division. iiFu Lili , ae ncmlt iohnra eoe provides genomes mitochondrial complete on based ) .andersoni N. Niviventer 2 saseisedmct hn,adhstelretbd-iewe oprdwith compared when body-size largest the has and China, to endemic species a is iu Zhou Jihua , iietrandersoni Niviventer .andersoni N. .andersoni N. lsee with clustered iietrandersoni Niviventer , .confucianus N. 1.Te iei aiu id ffrs nbt olnsadmountains and lowlands both in forest of kinds various in live They [1]. 3 .excelsior N. ihuLv Jizhou , Dvso 6.Meanwhile, [6]. -Division Dvso n the and -Division Toa,11)(, 1911) (Thomas, stenx lss eaieadfre h itrgopof, group sister the formed and relative closest next the as n omasse ru with group sister a form and eog ogenus to belongs 1 4 hnyn Tan Zhongyang , .confucianus N. .niviventer N. .andersoni N. Niviventer .fulvescens N. .andersoni N. Niviventer Niviventer ssse to sister is Dvso 8.Pyoeei trees Phylogenetic [8]. -Division 1 .confucianus N. and uziZhang Yunzhi , Niviventer pce r itnuse from distinguished are species sa pce nei oCia In China. to endemic species an is ) sn etgnrto sequencing. next-generation using .fulvescens N. aiyMrde n order and Muridae, family , and .excelsior N. Niviventer .cremoriventer N. a eaae into separated was n hystatistically they and , u eerhmay research Our . eeclustered were 2 and , were Posted on Authorea 7 Jul 2021 — The copyright holder is the author/funder. All rights reserved. No reuse without permission. — https://doi.org/10.22541/au.162566526.66282591/v1 — This a preprint and has not been peer reviewed. Data may be preliminary. i.1 h ioeoeognzto of organization used mitogenome and the obtained 1, were reads Fig. clean quality of high mitogenome 1,578,672 of of total sequence a the reads, assemble raw to the filtering quality After iterations. 1,000 X with MEGA test organization using bootstrap Genome performed a 3.1 also via was discussion evaluated rat were and of tree Results species [24]. ML 13 3. the model in of substitution PCGs values nucleotide 13 support best Aikake the The the corrected of [25]. determine using analysis evaluated to (ML) before were 2.1.10 likelihood settings evolution jModelTest Maximum default of in with Models (AICc) ClustalX Criteria [22,23]. using Information aligned 7.2) (Version were DAMBE sequences nucleotide by The concatenation 1). (Table comparing GenBank from performed skew was CodonW GC analysis using and Phylogenetic calculated T) were + values (A (RSCU) of analysis (A-T)/ usage Phylogenetic map = codon 2.3 genome skew synonymous mitochondrial AT Relative Compo- circular Search The [18]. formulas: [19]. tRNAscan-SE source [20]. the C) the 1.4.2 code + to genetic using (G according mitochondrial identified (G-C)/ calculated vertebrate = were was the using analysis genes mode, skew were tRNA mitochon- search sition (PCGs) The default complete genes the a coding under code. into Server BLASTx ( Protein genetic assembled using Finder [17]. novo mitochondrial ORF database 2019.1.1) de vertebrate NCBI sequence (Version and the package protein picked using software non-redundant were annotated Geneious 0.39) NCBI (Version reads using Trimmomatic to instrument genome Mitochondrial using aligned HiSeq-PE150 drial filtered were [16]. on and purified reads DIAMOND sequenced trimmed The cleaned and by were library reads [14]. The raw sequencing described The previously the [15]. as construct China). (rPCR) to Beijing, PCR used (TIANGEN, random were to products subjected kit rPCR was extraction DNA mitochondrial mitochondrial using The muscle annotation and the assembly from sequencing, extracted Mitogenome laboratory was 2.2 the DNA to -80degC. transported -80degC. at then at directly. stored storage, tubes stored and short cryopreservation for (Solarbio) and the nitrogen in ice liquid kept in dry were put muscle in organ immediately and for were kidney dissected samples lung, and the spleen, sacrificed All liver, were heart, individuals The These collection. 2018. August in H=1875.43m), 102deg10’15.7296”, of extraction Individuals DNA genomic and genetic collection useful Sample information. 2.1 provide mitogenome findings methods complete Our and with Materials species species. 2. other rat to the other phylogeny comparisons evolution several of phylogenetic the of with features for of rate relationships data the understanding fast phylogenetic provided our its and has limited compare recombination has study data of The genetic level of . low lack inheritance, The animal maternal diverse of [10-13]. size, of rodents) analyses small in phylogenetic genus its and (particularly the taxonomic to within for due species used 4 groups, been of have [21,24]. genomes mitogenomes species mitochondrial complete related Complete the closely resolution between now, single improved relationships to to provide evolutionary Compared Up can the [9]. into sequences scales investigations genome time for mitochondrial and sensitivity pressures complete and evolutionary sequences, different gene under evolves mitochondrial gene each since data, .andersoni N. .andersoni N. .andersoni N. .andersoni N. .andersoni N. ntepeetsuy esqecdtecmlt iohnra eoeof genome mitochondrial complete the sequenced we study, present the In . edpstdi nNB ihGnakacsinnme W314 The MW030174. number accession GenBank with NCBI in it deposited we eecletdfo uegCut,Yna rvne hn (24 China province, Yunnan County, Lufeng from collected were iohnra eoe fe bann h opeemtcodilgenome mitochondrial complete the obtaining After genome. mitochondrial a iclrDAmlcl hc a 621b nlnt.A hw in shown As length. in bp 16,291 was which molecule DNA circular a was .andersoni N. https://www.ncbi.nlm.nih.gov/orffinder/ .andersoni N. .andersoni N. 2 .andersoni N. a iia ota fms l te oet [26]. rodents other all most of that to similar was iohnra eoeadhsalwdu to us allowed has and genome mitochondrial n 2ohrrtmtgnmsdownloaded mitogenomes rat other 12 and a rw sn GRW131[21]. 1.3.1 OGDRAW using drawn was Niviventer eeaalbei GenBank. in available were n LSxwt the with BLASTx and ) .andersoni N. ° 57’45.774”; Posted on Authorea 7 Jul 2021 — The copyright holder is the author/funder. All rights reserved. No reuse without permission. — https://doi.org/10.22541/au.162566526.66282591/v1 — This a preprint and has not been peer reviewed. Data may be preliminary. rmtepeiu td sn nytect ee xetfrtepyoeei oiinof position those phylogenetic that with the sequence showed for congruent studies complete except previous generally the gene, in were on trees cytb results gene based the cremoriventer identical.Our cytb only constructed Single were using trees preform genome study . to likelihood previous mitochondrial maximum used (except this the the genome in were the from mitochondrial rate of complete species of mutation PCGs the rapid 13 topologies on the and based The of constructed 4). because was excluded tree D-loop). was likelihood (Fig. region maximum genome D-loop The The region. mitochondrial 1). complete (Table analysis the phylogenetic on of relationships based phylogenetic When reconstructed the same. investigate the further almost To were analyses phylogenetic on two species, based the genus tree rat of within and phylogenetic results other method heterogeneity another The high with likelihood constructed 3B). its maximum compared we to (Fig. Thus, by due ND6 analysis tree excluding [29]. phylogenetic PCGs phylogenetic performance during phylogenetic a excluded be poor established should consistently we sequences gene species, set ND6 rat 3). which that (Table 13 replications Gs of 1,000 and PCGs As with than 13 rRNA on the (0.204) AT-skew Based in its present and were 63.43%, Cs was analysis and rRNA Phylogenetic As of 3.4 content more on A+T that identified The showed were many (-0.099) respectively. subunits, GC-skew bp, of ribosomal 957 and large mitogenomes and and the 1,567 small are in the lengths of common encoding strand is srRNA) L- arm (lrRNA, the genes DHU had rRNA the tRNAs that tRNA two of eight arm the The genes, Loss dihydroxyuridine All tRNA a 2). lacked shape. these (Table H-strand which Among ring [28]. the trns1, a bp. on except to 75 encoded structure, were simplified to cloverleaf 14 bp been typical remaining 59 a the from and exhibited length L-strand genes in the ranging on encoded biased, were A+T may highly This be to 2). used of (Fig. frequently mitogenome genes acids The most RNA amino Transfer three and PCG the codons. PCGs. RNA frequent were 3,805 of Ribosomal most (M) AT-skew has 3.3 positive AUA the region and were and gene GC-skew (I) threonine protein-coding negative AUU and the the (L), explain isoleucine that CUA shows Leucine, also analyses, values (RSCU) which codons. RSCU usage 4, the codon Table synonymous to relative in According The displayed 2). are (Table CAT PCGs with of codons terminated Initiation which ND2 mitogenome. for of the except mitogenome of 70.1% the for of accounted of which PCGs mitogenome bp, in 11,420 PCGs was all PCGs the of 13 in the the of T whereas length Total than G than frequently (PCGs) [27]. common more genes genome more occurs Protein-coding mitochondrial is 3.2 A C of that that transcription biases 3). suggesting indicating and (Table composition (-0.347), (0.092) replication genome respectively, base value mitochondrial the positive 38.4%, Such skew and in was GC 61.6% role skewness 3). negative of (Table AT vital a percentage AT a showed GC towards play also and biased AT to It as is the reported makes on genome been which mitochondrial encoded have T, were for the 30.0% tRNAs that and 8 indicating G and for ND6 12.1% while C, of (H)-strand, composition base Heavy total the The (Table on rRNAs encoded 2 were and tRNAs (L)-strand. genes 22 Light of PCGs, 13 Most including identified, 2). were genes mitochondrial typical Thirty-seven Niviventer hnto than .andersoni N. .andersoni N. .andersoni N. . .andersoni N. .andersoni N. .andersoni N. n eelctdbtentRNA between located were and , otie h yia 2tN ee hogottegnm n appeared and genome the throughout genes tRNA 22 typical the contained u musculus Mus .andersoni N. and .excelsior N. emntdwt opee(A)o rnae T tpcodons, stop (T) truncated or (TAA) complete with terminated iohnra eoewsetmtdt e3.%frA 58 for 25.8% A, for 33.7% be to estimated was genome mitochondrial eetpclAN xetfrN1 hc tre ihGG All GTG. with started which ND1, for except ATN, typical were a hlgntclycoe to closer phylogenetically was sotru Fg A.Sm eerhr aesuggested have researchers Some 3A). (Fig. outgroup as 3 63,1.Orrslssgetthat suggest results Our [6,30,31]. .andersoni N. .confucianus N. Phe n tRNA and h hlgntcrltosiswere relationships phylogenetic the , a lsrto closer was Leu .excelsior N. h rN n srRNA and lrRNA The . .fulvescens N. .andersoni N. .confucianus N. n clustered and .andersoni N. and and N. Posted on Authorea 7 Jul 2021 — The copyright holder is the author/funder. All rights reserved. No reuse without permission. — https://doi.org/10.22541/au.162566526.66282591/v1 — This a preprint and has not been peer reviewed. Data may be preliminary. i:Mrde nCiaifre rmcmlt iohnra yohoebgn.Mlclrphylogenetics Molecular gene. b (Roden- cytochrome Niviventer mitochondrial genus doi:10.1016/j.ympev.2007.04.003 in complete (2):521-529. relationships from Phylogenetic 44 (2007) inferred evolution X China and Reference Zheng in W, Geographic Wang Muridae) and SH, Taxonomic Wu tia: HT, A Yu M, World. pp, Jing the 2142 6. of Press, University Species Hopkins (1):134-135 Mammal Johns 111 (2005) Anthropologie DM ed). archaeology. (3rd Reeder pleistocene DE, Chinese Wilson in 5. research 43:333–359 Current Perspect (2007) Asian Panxian JJ at China. Bahain abundance Taxonomic South 4. (2004) in H Cave Weiwen Pleistocene H, Yamei middle S, rat a Miller-Antonio endemic Dadong, and LA, variables Schepartz An D, morphometric (Thomas, Bekken (2017) cranial andersoni) with 3. Q (Niviventer verified Yang Niviventer China 2196:48-58 Q, Anderson’s in Zootaxa the subspecies Liu new characteristic. of Y, two variation pelage of Du Geographic Muridae) X, (2009) (Rodentia: J Lv 1911) Yang Z, China S, doi:10.1038/srep46127 Wen of Li centre 7:46127. Y, biodiversity 2. reports Chang terrestrial Scientific the L, in Quaternary. Xia changes late environmental J, the moderate through of Cheng evidence L, is Lu complex mi- species D, complete Ge The 1. sequences: DNA of availability References the regarding supplied of was togenome information following DLDXLL2017007. The number: approval ethics animal Institute Yunnan the accessibility the with by Data approval Prevention ethics and the under Control were Disease study this Endemic in of experiments the all and collection Sample participate interests. to of consent conflict and no talents approval high-level is of Ethics there project that training declare health “Research Yunnan authors China National the The the of and and Program (81660558) (L-2017027). 2017YFF0210302), Development China SN. of and Foundation (project/subject Research Science infrastructure” Key Natural quality National national the of by application funded and jointly was work This interest of declaration and genus Acknowledgements the within evolution understand rodents. murid better other that to showed helpful be mitogenome should complete the on cremoriventer based D-loop, analysis without Phylogenetic confucianus were mitogenome features characteristics. structural complete morphological genome the to mitochondrial similar and most The ND6 the genus excluding Muridae. in family PCGs species the the of other to species of related similar genome closely mitochondrial with complete it the complete sequenced and have We tree population gene same mitochondrial genomes the Conclusion mitochondrial single for complete The 4. trees using gene tree [9]. other phylogenetic genes from that warranted. differ the suggesting was may conflicting, of were population selection tree a subjective mitogenome for the tree on gene depending one that known been cremoriventer N. excelsior N. a h lss eainhpto relationship closest the had lsee oehr hnwith then together, clustered iietrandersoni Niviventer speiul ugse.Teaalblt fcmlt iohnra eoeof genome mitochondrial complete of availability The suggested. previously as ic ahgn vle ne ieeteouinr rsueadtm cl,i has it scale, time and pressure evolutionary different under evolves gene each Since . .excelsior N. ossetwt rvossuisbsdo igect eesqecsand sequences gene cytb single on based studies previous with consistent , Niviventer sdpstdi eBn fNB ne ceso ubrMW030174. number accession under NCBI of GenBank in deposited is .confucianus N. .andersoni N. ntepyoeei nlsso eune fte1 PCGs, 13 the of sequences of analysis phylogenetic the In . 4 and n omdasse ru of group sister a formed and , .andersoni N. .excelsior N. Niviventer, ahrthan rather , o h rttm n compared and time first the for .andersoni N. swl sisrltosi to relationship its as well as .fulvescens N. .fulvescens N. a consistently was .andersoni N. and and N. N. Posted on Authorea 7 Jul 2021 — The copyright holder is the author/funder. All rights reserved. No reuse without permission. — https://doi.org/10.22541/au.162566526.66282591/v1 — This a preprint and has not been peer reviewed. Data may be preliminary. M imA oe ,Topo D isnT,HgisD 20)CutlWadCutlXvrin2.0. version X Clustal and W Clustal (2007) Wallace DG F, Valentin Higgins H, TJ, doi:10.1093/bioinformatics/btm404 McWilliam Gibson PA, JD, (21):2947-2948. McGettigan Thompson 23 R, R, Bioinformatics Chenna Evolution. Lopez NP, A, Molecular Brown Wilm and G, IM, Blackshields Phylogenetics, MA, Genomics, Larkin Microbial doi:10.1093/jhered/esx033 23. for (4):431-437. (W1):W59-W64. Tools 108 47 New heredity research of DAMBE6: Journal acids (2017) The Nucleic X Xia expanded genomes. of 22. 1.3.1: organellar study version the of (OGDRAW) in OrganellarGenomeDRAW visualization progress (2019) graphical and doi:10.1093/nar/gkz238 R the method Bock for analysis P, Lehwark [The toolkit doi:10.1360/yc-007-0420 S, (2007) (4):420-426. FC Greiner 29 He 21. Hereditas YP, = Zhu chuan DM, animal Yi Ren of bias]. sites SF, degenerate Wu codon fold XM, four Wu 41:353–358 at 20. Evolution composition Molecular nucleotide of of Patterns Journal Methods (1995) genomes. Sequences. TD mitochondrial Genomic Kocher NT, in Perna Genes tRNA 19. and for doi:10.1007/978-1-4939-9173-0 integrated Searching 1962:1-14. an tRNAscan-SE: biology (2019) molecular Basic: TM 28 in Lowe Geneious Bioinformatics PP, data. (2012) Chan Markowitz sequence A A, 18. of Drummond Cooper analysis S, P, and Buxton organization Meintjes doi:10.1093/bioinformatics/bts199 S, the B, Sturrock for (12):1647-1649. Ashton M, platform software Cheung T, Nature desktop S, Thierer DIAMOND. Stones-Havas extendable using A, C, Wilson alignment Duran R, protein S, Moir sensitive M, and data. Kearse Fast sequence 17. (2015) Illumina DH doi:10.1038/nmeth.3176 for Huson (1):59-60. trimmer C, 12 flexible Xie methods a B, Buchfink Trimmomatic: prevalence, range, (2014) 16. Host B (2010) Usadel doi:10.1093/bioinformatics/btu170 Z Shi M, (15):2114-2120. LF, 30 doi:10.1128/JVI.02497-09 Lohse Wang Bioinformatics J, (8):3889-3897. AM, Yuan 84 Or- (9):2896-2913. Bolger Y, virology mt-Genome Zhang of 8 Journal Y, 15. Zhu It: bats. evolution P, in Zhou Know adenoviruses and H, of Not diversity Zhang biology genetic X, Do and Ge Genome We Y, Li as 14. DNA Lineages. Mitochondrial Nonbilaterian Animal in and (2016) rules doi:10.1093/gbe/evw195 Evolution W DNA: and Pett mitochondrial DV, animal ganization Lavrov of inheritance doi:10.1186/s40709-017-0060-4 and 13. 24:2. research Evolution biological (2017) of E Journal Biochem Zouros exceptions. ED, Muridae). Ladoukakis complete (Rodentia: of 12. characterization genus The (2017) doi:10.1016/j.bse.2017.01.012 SD within Chen 71:179-186. H, relationship (Cephalobo- Ecol Zong National varius phylogenetic Syst GY, the Chen Acrobeloides and Q, of Wang of genome XM, genome Wang Proceedings mitochondrial FJ, mitochondrial doi:10.7717/peerj.9108 Li The HX, 8:e9108. kingdom. Wei (2020) PeerJ Fungi 11. JK (Nematoda). Park the Tylenchina HJ, for of doi:10.1073/pnas.1711939114 Kil non-monophyly Life Y, confirms (35):9391-9396. of Lee morpha) 114 Tree T, America Kim of genome States A 10. United (2017) the of SH Sciences Bul- Kim of . J, Academy Indo-Malayan the 168((3):236–256 Choi of and History systematics Sulawesi, Natural on fromCeylon, 9. of Notes species Museum 105. and American genera No. the new of of Expeditions. of transmission letin Archbold descriptions and the and of emergence rodents, Results the murid (1981) on GG biodiversity Mitchell Musser KE, of doi:10.1038/nature09575 Jones 8. Impacts A, (7324):647-652. Jolles (2010) P, 468 RS Hudson Nature RD, Ostfeld Holt diseases. T, CD, infectious Harvell Bogich A, SS, Dobson P, Myers Daszak CE, LK, Belden F, Keesing 7. 5 1 Posted on Authorea 7 Jul 2021 — The copyright holder is the author/funder. All rights reserved. No reuse without permission. — https://doi.org/10.22541/au.162566526.66282591/v1 — This a preprint and has not been peer reviewed. Data may be preliminary. al opeemtcodilgnmsue o hlgntcaayi nti study this in analysis phylogenetic for used genomes mitochondrial Complete 1 Table genes tein-coding 4 Table 3 Table andersoni Niviventer 2 Table 1 Table D-loop. the values. minus RSCU genome all 4. of sum Fig. the represents above species. column rat the 13 of of height the and 3. acid, amino of each mitogenome encoding the in of (RSCU) phylogeny, China 1. Molecular in (2018) Fig. species Q Yang rat AV, wild Abramov doi:10.1093/jmammal/gyy117 common J, Legends (6):1350-1374. of Cheng Figure 99 Muridae) , complex Mammalogy Wen species (Rodentia, of hiXin a Journal sacer Z of Y, ). confucianus revision Du (Rodentia, Niviventer 959:137-159. systematic L, (2020) and Xia ZooKeys diversity, Li, Y morphological Lu evidence. Li D, morphological X, Ge and Rong 31. karyotyping, J, molecular, Wang on H, based Li species doi:10.3897/zookeys.959.53426 Y, distinct a Li a phylogenetics: Y, and is molecular phylogenetics Li teleostean Molecular in 30. information criterion. mitogenomic doi:10.1006/mpev.2000.0839 optimality of (3):437-455. maximum-parsimony Use 17 the evolution (2000) under M exploration Nishida of tree-based M, (8):1767-1780. strand review Miya International on views 27 evolution. 29. New and (2010) research structure JH DNA: doi:10.1016/s0074-7696(08)62066-5 He mitochondrial acids 141:173-216. GY, doi:10.1371/journal.pone.0012708 Animal cytology Ye (1992) Nucleic (9):e12708. C, DR 5 Achterberg Wolstenholme one Evolutionary van 28. PloS MJ, Molecular Sharkey genomes. genomes. mitochondrial XX, X: insect Chen mitochondrial in M, MEGA (6):1547-1549. asymmetry Shi Animal SJ, (2018) 35 Wei (1999) 27. evolution K JL Tamura and C, biology Boore doi:10.1093/nar/27.8.1767 Knyaz Molecular M, 26. Li Platforms. and G, heuristics Computing new Stecher across doi:10.1093/molbev/msy096 models, S, more Analysis 2: Kumar Genetics jModelTest (2012) doi:10.1038/nmeth.2109 D (8):772. 25. Posada 9 R, methods Doallo Nature GL, computing. Taboada parallel D, Darriba 24. h aiu ieiodaaye fpyoeei eainhpbsdo A 3PG n B 2PCGs 12 (B) and PCGs 13 (A) on based relationship phylogenetic of analyses likelihood maximum The iohnra eoempof map genome Mitochondrial eu pce omnnm GenBank name Common confucianus Niviventer Niviventer Species edwardsi Leopoldamys Genus u u musculus Mus Mus uloiecmoiinadA-Csens fthe of skewness AT-GC and composition Nucleotide of genome study mitochondrial this the in of analysis Characteristics phylogenetic for used genomes mitochondrial Complete h aiu ieiodaaye fpyoeei eainhp ae ncmlt mitochondrial complete on based relationships phylogenetic of analyses likelihood maximum The eaiesnnmu oo sg n oo ubr in numbers codon and usage codon synonymous Relative iietrexcelsior Niviventer cremoriventer Niviventer sabanus Leopoldamys iietrandersoni Niviventer iietrandersoni Niviventer hns ht-ele a NC rat white-bellied Chinese NC rat giant long-tailed Edwards’s os os NC NC NC MN964122.1 rat white-bellied large rat tree dark-tailed mouse house rat giant long-tailed 6 h o eo h a hr ersnsalcodons all represents chart bar the below box The . . i.2. Fig. iietrandersoni Niviventer iietrandersoni Niviventer h eaiesnnmu oo usage codon synonymous relative The mitogenome iohnra pro- mitochondrial 019617.1 035822.1 023960.1 005089.1 025670.1 Fig. Posted on Authorea 7 Jul 2021 — The copyright holder is the author/funder. All rights reserved. No reuse without permission. — https://doi.org/10.22541/au.162566526.66282591/v1 — This a preprint and has not been peer reviewed. Data may be preliminary. -op14312189H64.2% 62.7% H 68.7% L 69.6% H 58.3% L 61.3% 66.2% 62.7% 61.3% 73.5% H L H H 80.9% H 62.0% H TAA 63.3% 64.7% TAA H 62.9% H TAA H 57.4% H ATG 889 ATG H 66.2% 16291 62.7% ATA 15403 64.2% T– H 67 82.4% TAA 1159 15402 H 67 59.4% TAA H 15283 59.5% 519 15336 15335 D-loop H 14125 1830 14050 69 15269 H 54.5% ATG tRNA-Pro 13554 13532 14119 ATG 50.0% T– L tRNA-Thr H 11725 14051 TAA ATC 66.2% 59.7% CYTB 69.6% L TAA 71 tRNA-Glu 1378 L 59 63.6% 11724 ND6 ATG 68 11527 297 H L TAA 11654 11654 ND5 ATG 11595 10150 L 10156 348 11596 52.2% tRNA-Leu1 ATG 11528 H 62.0% 68 9860 tRNA-Ser1 9788 63.7% 72.5% tRNA-His TAA 9857 784 ATG 9441 H 68 ND4 681 9790 L 9372 H 57.6% ND4L H 9440 204 54.7% 63.9% 8589 tRNA-Arg 8589 62.6% 9373 7951 58.8% 67.6% ATG 7909 ND3 684 65 H 7748 tRNA-Gly H H CAT 7679 COX3 7747 H H H 68 ATP6 6996 1545 7683 69 ATP8 6994 6857 tRNA-Lys T– ATC 6923 6927 66 COX2 5313 68 6855 tRNA-Asp 5311 71 tRNA-Ser2 5245 5246 1036 69 GTG COX1 5143 5178 66 tRNA-Tyr 4935 5071 5073 tRNA-Cys 5001 3900 5003 69 tRNA-Asn 955 4936 71 tRNA-Ala 69 3899 tRNA-Trp 3690 3827 3759 3831 A+T ND2 2736 3757 1567 Strand tRNA-Met 3691 75 Codon 957 tRNA-Gln 2660 Stop 2735 68 68 tRNA-lle Codon 1025 Start 1094 2661 ND1 1093 Length(bp) 68 tRNA-Leu2 1026 69 Stop l-rRNA tRNA-Val Start 1 s-rRNA tRNA-Phe Gene atsRtu andamanensis Rattus Rattus eu pce omnnm GenBank name Common Species Genus atstiomanicus Rattus tanezumi Rattus norvegicus Rattus baluensis Rattus iietrfulvescens Niviventer aaa edrtMN126562.1 NC NC NC NC rat field Malayan rat house Oriental rat Norway rat forest Indochinese hsntwieblidrtNC rat white-bellied Chestnut 7 011638.1 001665.2 035621.1 046686.1 028715.1 Posted on Authorea 7 Jul 2021 — The copyright holder is the author/funder. All rights reserved. No reuse without permission. — https://doi.org/10.22541/au.162566526.66282591/v1 — This a preprint and has not been peer reviewed. Data may be preliminary. A()5 .2AUS 906 U()1 .4UCF 3 1.12 1.72 0.88 132 0.28 UUC(F) 104 81 1.34 UUU(F) GAA(E) 0.24 13 1.67 0.66 GAG(E) 51 1.76 0.80 0.33 GAC(D) 16 25 25 1.26 0.17 mitochondrial 116 UUG(L) GAU(D) UGC(C) 1.50 using 91 0.13 0.85 reconstructed 5 UUA(L) were CUU(L) 76 19 0.67 Muridae UGU(C) 2.19 0.14 of 1.11 91 0.00 GCA(A) 1.33 97 lineages CUG(L) 8 0.38 major 1.86 GCC(A) 251 among genomes. 0.00 CUC(L) GCG(A) 19 of relationships 0.92 annotation 7 0.06 and 67 Phylogenetic CUA(L) sequence 4.00 38 genome 0 1.08 AGU(S) mitochondrial AAG(K) 0.88 complete GCU(A) 0.54 94 0.00 6 first AGC(S) The 170 AGG(*) 0 2.69 AAA(K) 0.92 CGU(R) 200 0.12 1 1.46 8 AUC(I) Highlights AGA(*) 1.08 14 26 1.79 CGG(R) AUU(I) 1.88 0 0.15 UAA(*) 43 CGC(R) CAU(H) 59 2.18 1.12 71 UAG(*) 5 CGA(R) 1.85 69 1.04 94 UAU(Y) 0.42 1.75 81 CAC(H) 0.08 CAG(Q) 109 GGA(G) UAC(Y) 8 0.67 59 1.01 CAA(Q) RSCU 98 UGG(W) CCA(P) 52 0.70 GGC(G) 0.37 22 0.69 71 Count UGA(W) 4 0.86 CCC(P) GGG(G) 35 41 0.08 GUA(V) 1.31 CCG(P) GGU(G) Codon 35 15 GUC(V) 54 1.11 1.72 RSCU 35 CCU(P) GUG(V) 102 0.76 AAU(N) 0.28 Count 6 GUU(V) 198 AAC(N) 86 2.05 ACG(T) 0.87 AUA(M) Codon 59 32 ACC(T) 0.16 RSCU 159 AUG(M) ACU(T) 1.14 Count 55 ACA(T) 1.83 10 UCU(S) Codon 72 UCG(S) RSCU 115 UCC(S) Count UCA(S) Codon in numbers codon genes and coding usage codon synonymous Relative 4 Table the of skewness AT-GC and composition Nucleotide of 3 genome Table mitochondrial the of Characteristics 2 Table RA 543.91.82.42.96.3024-0.099 0.057 -0.365 -0.370 0.204 0.067 63.43 0.068 -0.347 64.91 20.09 0.030 64.23 16.54 25.24 56.25 24.41 30.29 16.48 25.01 29.92 18.55 38.19 27.29 11.36 0.092 34.62 11.51 34.31 61.62 28.96 25.85 2524 27.97 1499 889 12309 12.53 33.65 16291 region Control rRNAs tRNAs PCGs mitogenome andersoni Niviventer ie(p + Tsens Cskewness GC skewness AT A+T C T G A (bp) Size 8 iietrandersoni Niviventer iietrandersoni Niviventer iietrandersoni Niviventer iietrandersoni Niviventer mitogenome iohnra protein- mitochondrial . Posted on Authorea 7 Jul 2021 — The copyright holder is the author/funder. All rights reserved. No reuse without permission. — https://doi.org/10.22541/au.162566526.66282591/v1 — This a preprint and has not been peer reviewed. Data may be preliminary. figures/Figure1/Figure1-eps-converted-to.pdf 9 Posted on Authorea 7 Jul 2021 — The copyright holder is the author/funder. All rights reserved. No reuse without permission. — https://doi.org/10.22541/au.162566526.66282591/v1 — This a preprint and has not been peer reviewed. Data may be preliminary. figures/Figure2/Figure2-eps-converted-to.pdf 10 Posted on Authorea 7 Jul 2021 — The copyright holder is the author/funder. All rights reserved. No reuse without permission. — https://doi.org/10.22541/au.162566526.66282591/v1 — This a preprint and has not been peer reviewed. Data may be preliminary. figures/Figure3/Figure3-eps-converted-to.pdf 11 Posted on Authorea 7 Jul 2021 — The copyright holder is the author/funder. All rights reserved. No reuse without permission. — https://doi.org/10.22541/au.162566526.66282591/v1 — This a preprint and has not been peer reviewed. Data may be preliminary. figures/Figure4/Figure4-eps-converted-to.pdf 12