- Research Article
- Open Access
Whole mitochondrial genomes provide increased resolution and indicate paraphyly in deer mice
BMC Zoology volume 2, Article number: 11 (2017)
Recent phylogenies of deer mice, genus Peromyscus, have relied heavily on mitochondrial markers. These markers provided resolution at and below the level of species groups, but relationships among species groups and Peromyscus affiliated genera have received little support. Here, we present the mitochondrial genomes of 14 rodents and infer the phylogeny of Peromyscus and related taxa.
Our analyses support results from previous molecular phylogenies, but also yield support for several previously unsupported nodes throughout the Peromyscus tree. Our results also confirm several instances of paraphyly within the clade and suggest additional taxonomic work will be required to clarify some relationships.
Our findings greatly enhance our understanding of the evolution of Peromyscus providing support for previously unsupported relationships. However, the results also highlight the need to address paraphyly that may exist in this clade.
Mitochondrial loci have been the most popular phylogenetic markers in animals for over three decades. Their ease of amplification [1, 2], uniparental inheritance, lack of recombination in mammals [2,3,4], differential selection among genes , and general synteny  combine to make mitochondrial phylogenetic markers an excellent choice for many study systems. Recent increases in DNA sequencing capabilities have provided an opportunity to sequence full mitochondrial genomes (mitogenomes) rather than partial, one, or a few mitochondrial genes. As a result, previous phylogenies that relied on a portion of the mitochondrion are now undergoing re-analysis in an effort to expand support for or re-affirm previous conclusions .
The genus Peromyscus, commonly referred to as North American deer mice, encompasses more than 70 species that diverged within the last 6–10 million years . Species including P. maniculatus  and P. leucopus  are among the most common mammals in North America and have been studied extensively for over 100 years . Despite extensive study, new species  and subspecies  are being described on a regular basis. The large number of species, both described and undescribed, as well as substantial cryptic variation, has yielded numerous distinct phylogenetic hypotheses.
Previous classifications of Peromyscus recognized seven morphologically distinct subgenera (Habromys, Haplomylomys, Isthmomys, Megadontomys, Osgoodomys, Peromyscus, and Podomys) [13, 14]. Many of those subgenera (Habromys, Isthmomys, Megadontomys, Osgoodomys, and Podomys) have at times been elevated to generic status . Some recent molecular phylogenies, however, show paraphyly . Isthmomys is sister to the genus Peromyscus, representing a distinct genus, as is currently accepted, but Habromys, Megadontomys, Osgoodomys, Neotomodon, and Podomys, distinct genera, were found within Peromyscus, rendering the genus paraphyletic [7, 16, 17], though Neotomodon alstoni  previously had been Peromyscus alstoni [14, 19].
Recent attempts to resolve the history of this clade by Bradley et al. , Miller and Engstrom  and Platt et al.  were based on various combinations of nuclear and mitochondrial genes. All three studies arrived at similar topologies, but varied in their levels of support at some nodes. Species-level relationships were well-supported but the relationships among some genera and species groups were not. Given that approaches that use cytochrome b (cytb) alone were somewhat successful in resolving Peromycus, we reasoned that the use of whole mitogenomes, a more data-rich approach, may be an avenue to more robust results at the generic and species group level, thus providing additional clarity. Here, we analyze whole mitogenomes from Habromys, Isthmomys, Neotomodon, Peromyscus, Podomys, and three outgroups, Sigmodon, Neotoma, and Reithrodontomys, as identified from previous molecular phylogenies [7, 16, 17].
Sampling, DNA preparation and sequencing
Taxa selected for this study were identified based on the findings of Platt et al. . Our objective was to reduce taxon sampling for some groups to isolate specific clades. To do so, we chose representatives from selected subclades to serve as proxies. Specifically, one individual here is a proxy for a larger clade (generally represented as a species group). Members of Neotoma, Reithrodontomys, and Isthmomys, which have determined to be closely related clades by several morphological and molecular analyses [7, 15, 16] were included. Species identifications, museum catalog numbers and are included in Table 1.
Whole genomic DNA was isolated using a Qiagen DNeasy Blood and Tissue kit (Qiagen, Valencia, CA) via the manufacturer’s protocol and fragmented to ~ 400 bp. Sizes were verified on an Agilient 2100 Bioanalyzer. Illumina compatible sequencing libraries were prepared from the fragmented DNA using the NEB/KAPA library preparation. Each sample was tagged with a unique index and pooled in equal proportions, after which, the pooled libraries were sequenced on single run of an Illumina MiSeq (2 × 250 bp). Raw data are available in GenBank under BioProject ID PRJNA308567.
Mitochondrial genome assembly and annotation
Raw sequence reads were filtered and processed using Trimmomatic version 0.35 . Specifically, we clipped Illumina adapters, disregarded read ends whose Phred scores fell below 20, and utilized a four-base sliding window to trim reads once the average quality fell below 25. We assembled mitogenomes through a custom Bash script (https://github.com/KevinAMSullivan/Mitochondrial_Genomes/tree/MIRA) The script utilized two major programs, MITObim  and MIRA  to map the filtered reads to a reference genome before assembling the mitogenomes. We selected Akodon montensis , a Sigmodontine rodent and the most closely related organism with a fully sequenced mitogenome, as the reference (GenBank accession number KF769456).
MITOS  was used to annotate the mitochondrial genes for each genome. Putative genes were submitted to BLASTn to confirm sequence length. When gaps were noted, we manually checked for frameshifts but none were observed. Acceptable results were those whose top hits were to a Peromyscus or closely related species and that covered nearly 100% of the putative gene. The putative gene was shortened or lengthened to match the length of the BLAST hit if a different size. All sequences were deposited in GenBank under the accession numbers KY707299–707312. Genome coverage was estimated using bedtools genomecov  in combination with mapping of the processed reads to each MITObim assembly with BWA via a custom script (https://github.com/KevinAMSullivan/Mitochondrial_Genomes/tree/MIRA).
To test the reliability of mitogenome assemblies, we used two cytb genes per taxon from GenBank as controls (Additional file 1: Table S1), and worked under the assumption that our reference assembled cytb gene and those from GenBank should form monophyletic clades. We chose full length, non-identical entries from the selected species, aligned them with cytb from our assemblies, and used RAxML  to estimate the best tree from 1000 maximum likelihood (ML) searches. Support for each node was generated with bipartition frequencies from 10,000 bootstrap replicates. The GTR+GAMMA+I model of nucleotide substitution was used for both ML analyses.
Phylogenetic analysis of whole mitochondrial genomes
We concatenated the mitochondrial protein coding sequences and aligned them with Muscle . That alignment is available as Additional file 2. Bayesian and ML phylogenetic analyses were accomplished using MrBayes  and RAxML , respectively. For the Bayesian analysis, four independent runs on five chains were implemented for 3,000,000 generations. We evaluated stability of the final tree by continuing to an average standard deviation of split frequency less than 0.01. The data were also partitioned by codon. No substitution model was specified as a prior to let the program search across all possible models to determine the most appropriate model of evolution. We then used the selected model, GTR+GAMMA, for our ML analysis, for which we implemented 1000 bootstrap replicates. In both analyses, the A. montensis and S. hispidus  mitogenomes were specified as outgroups when drawing the trees as Neotominae and Sigmodontinae are sister subfamilies whose relationships have been supported by previous analyses [30, 31]. Our unrooted tree supports this relationship. The implementation scripts for both analyses are available on github (https://github.com/KevinAMSullivan/MIRA-MitoBim/tree/Phylogenies).
We sequenced and assembled whole mitogenomes of 14 rodents in the subfamily Neotominae and one member of Sigmodontinae, Sigmodon hispidus. Taxon sampling included seven genera and six of the 13 species groups within Peromyscus . Genome assembly data, including accession numbers, read coverage, and mean read totals across all taxa, are listed in Table 1.
Each assembled mitogenome comprises 22 tRNA, 2 rRNA, and 13 protein coding regions, totaling 37 genes. Although most genes are transcribed from the positive (heavy) strand, several genes, mainly tRNAs but also including nad6, are transcribed from the negative (light strand). Additionally, the 3′ end of atp6 overlaps with the 5′ end of atp8, a common characteristic in mitogenomes. Noncoding regions typical of mitogenomes are present, including a control region downstream of cytb. As these mitogenomes were not closed assemblies, such variation in mitogenome size is thought to be due to the highly variable and often heteroplasmic control region.
Validation of MITObim assembled and MITOS-annotated cytb took three forms. First, reconstructed genes clustered with loci amplified and sequenced using more traditional methods (Additional file 3: Figure S1), suggesting that the remaining protein coding sequences were accurately assembled and annotated for each taxon. Second, differences that were present are within acceptable ranges for such rapidly evolving loci. Despite forming monophyletic clades for each species, some of the cytb genes differed from those in GenBank. For example, cytb from our P. crinitus  mitogenome exhibited 43 and 44 differences, respectively, when compared to their GenBank counterparts, i.e. ~ 96% identity. Divergences in the range may be expected for such a rapidly evolving locus within a species . Third, we have high confidence in our base calls. For example, average coverage of our cytb loci is 142× (data not shown), and the average percentage of those reads indicating any given nucleotide is 99.3%. Given such high support, as well as the fact that P. crinitus is the top BLAST hit to the gene (96% identify, 0.0 E value), we have confidence in the veracity of our cytb sequences over alternate hypotheses for such sequence dissimilarity. Similar results were observed for cytb sequences in all taxa (Additional file 4: Table S2). The single exception was for Neotoma mexicana. The best match (100%) for our reconstructed cytb sequence was to N. isthmica. However, N. isthmica is closely related to N. mexicana , and the second best hit was to N. mexicana (99%), confirming its validity as an outgroup .
ML and Bayesian phylogenies inferred from mitochondrial protein coding genes recovered identical topologies (Fig. 1). This topology is similar to recent molecular phylogenies. First, S. hispidus, N. mexicana, I. pirrensis , and R. mexicanus  are positioned outside of Peromyscus. Unlike previous phylogenies however, high posterior probabilities are found at these nodes. In fact, every node, save that linking P. crinitus and P. polionotus , has high support. Second, our analyses indicated a paraphyletic Peromyscus. Finally, Isthmomys and Reithrodontomys comprise a clade sister to what is currently considered to Peromyscus and its affiliated genera. This pairing is of note, as both genera have shown conflicting relationships depending on the marker used [7, 16, 17, 38].
Our phylogeny using whole mitogenomes largely mirrors that of Platt et al. , save the quartet of H. ixtlani , P. aztecus , P. pectoralis , and P. attwateri . Platt et al. found P. attwateri and P. aztecus to form a clade, with P. pectoralis as sister to this grouping and H. ixtlani being the most closely related species to a clade containing all three. The whole mitogenome tree is in partial agreement with Bradley et al.  with regard to these same four species. Our close pairing of P. attwateri and P. pectoralis matches, but their phylogeny places P. aztecus as a sister group to a clade that encompasses H. ixtlani, P. attwateri, and P. pectoralis.
Regardless of the differences or similarities, our tree (Fig. 1) provides support for previously unsupported nodes within the Peromyscus phylogeny. Previous trees were largely resolved, but lacked significant support at key nodes, mainly those in middle regions of the tree that identify relationships among species groups [7, 16, 17]. Platt et al.  provided substantial support (posterior probabilities ≥ 95% and bootstrap values ≥ 85) at nodes grouping the crinitus, eremicus, and californicus species groups; the mexicanus, megalops, and melanophrys species groups; and the aztecus and boylii species groups. However, support in the Bayesian mitogenome phylogeny has posterior probabilities (> 0.95) at every node and bootstrap support (> 90) in the ML phylogeny at all nodes save one.
Two branches of our inferred phylogeny had not yet been reinforced by any molecular analysis. The first concerns H. ixtlani, P. aztecus, P. attwateri, and P. pectoralis. Although some analyses had suggested a close relationship among these taxa, no node in this clade had received substantial support in any previous molecular phylogeny. Our analysis is the first to provide such support for these relationships.
The second branch of interest suggests a close relationship between Isthmomys and Reithrodontomys. Although no phylogenetic analysis ever suggested the genus Reithrodontomys should be nested within Peromyscus, Isthmomys had previously been subsumed as a subgenus . The pairing is counterintuitive given their morphological differences. Although both occupy the same geographic region, there are obvious incongruities in size. Isthmomys pirrensis can reach well over 100 g (averaging 140 g in one study) and 300 mm [42, 43], whereas rodents in Reithrodontomys are much smaller. R. mexicanus ranges from 167 to 190 mm in length , and a collection of different surveys gave a range from 7.9–9.5 g . Although their pairing has been previously suggested [7, 16, 17], our analysis is the first to provide high support values.
One major weakness of our study is its lack of complete taxon sampling. This makes comprehensive phylogenetic reconstruction uncertain due to missing critical (i.e. representing additional ingroups) taxa. The addition of excluded taxa, especially Osgoodomys and Megadontomys, should provide a more comprehensive view of phylogenetic relationships within the genus. That being said, the high support values we recover suggest that the relationships at the deeper nodes are valid.
The Peromyscus phylogeny has been studied repeatedly and revised often, using data from several sources . This latest attempt, using next-generation sequencing, continues to suggest paraphyly. Indeed, every molecular phylogeny has been unequivocal in suggesting paraphyly. Given this observation, addressing the status of genera subsumed within Peromyscus should be seen as a priority within rodent taxonomy.
Aside from acknowledging paraphyly, future work on Peromyscus phylogenetics should focus on further clarifying the true phylogeny via more taxa and markers. The inclusion of additional mitogenomes will elucidate relationships within the genus to provide increased detail. However, despite their advantages, mitochondrial loci are imperfect markers. Mitochondria genomes are susceptible to incomplete lineage sorting (ILS) and introgression . Hybrids, especially asymmetric hybrids, exemplify a similar problem in that they may all have the same mitochondrial DNA and species boundaries in those cases would not be reflected [45,46,47]. Supplementing the mitochondrial phylogeny with large numbers of nuclear markers such as ultraconserved elements (UCEs) [48, 49] or retrotransposons [50, 51] is an important next step in understanding phylogenetic relationships within Peromyscus.
Here, we present analyses of whole mitochondrial genomes, for 14 species, ten of which are Peromyscus or close relatives. The data yield phylogenies with significant support at previously unsupported nodes, particularly the pairing of Isthmomys and Reithrodontomys, but also suggest paraphyly within the genus that could be resolved by elevating monophyletic groups to genera or subsuming currently recognized genera as subgenera. Our analyses provide evidence that additional data will help clarify the evolutionary history of this genus.
- Cytb :
Falkenberg M, Larsson NG, Gustafsson CM. DNA replication and transcription in mammalian mitochondria. Annu Rev Biochem. 2007;76:679–99.
Galtier N, Nabholz B, Glemin S, Hurst GD. Mitochondrial DNA as a marker of molecular diversity: a reappraisal. Mol Ecol. 2009;18(22):4541–50.
Birky CW Jr. Uniparental inheritance of mitochondrial and chloroplast genes: mechanisms and evolution. Proc Natl Acad Sci U S A. 1995;92(25):11331–8.
Ballard JW, Whitlock MC. The incomplete natural history of mitochondria. Mol Ecol. 2004;13(4):729–44.
Boore JL. Animal mitochondrial genomes. Nucleic Acids Res. 1999;27(8):1767–80.
Williams ST, Foster PG, Littlewood DTJ. The complete mitochondrial genome of a turbinid vetigastropod from MiSeq Illumina sequencing of genomic DNA and steps towards a resolved gastropod phylogeny. Gene. 2014;533(1):38–47.
Platt RN, Amman AM, Keith MS, Thompson CW, Bradley RD. What is Peromyscus? Evidence from nuclear and mitochondrial DNA sequences suggests the need for a new classification. J Mammal. 2015;96(4):708–19.
Wagner A. Diagnosen neuer Arten Brasiliscer Handfluger. Arch Naturgesch. 1843;9:366.
Rafinesque CS. Further discoveries in natural history made during a journey through the western region of the United States. American Monthly Magazine and Critical Review. 1818;3:445–7.
Osgood WH. Revision of the mice of the American genus Peromyscus, vol. 28 Washington DC: North American Fauna; 1909.
Bradley RD, Carroll DS, Haynie ML, Martinez RM, Hamilton MJ, Kilpatrick CW. A new species of Peromyscus from western Mexico. J Mammal. 2004;85(6):1184–93.
Bradley RD, Schmidly DJ, Amman BR, Platt RN 2nd, Neumann KM, Huynh HM, Muniz-Martinez R, Lopez-Gonzalez C, Ordonez-Garza N. Molecular and morphologic data reveal multiple species in Peromyscus pectoralis. J Mammal. 2015;96(2):446–59.
Hooper ETMG. Notes on classification of the rodent genus Peromyscus. Occasional Papers of the Museum of Zoology, University of Michigan. 1964;635:1–13.
Wilson D, Reeder DM. Mammal species of the world, vol. vol. 2. 3rd ed. Baltimore: Johns Hopkins University Press; 2005.
Carleton MD. Phylogenetic relationships of neotomine–peromyscine rodents (Muroidea) and a reappraisal of the dichotomy within new world Cricetinae. Miscellaneous Publications of the Museum of Zoology, University of Michigan. 1980;157:1–146.
Bradley RD, Durish ND, Rogers DS, Miller JR, Engstrom MD, Kilpatrick CW. Toward a molecular phylogeny for Peromyscus: evidence from mitochondrial cytochrome-b sequences. J Mammal. 2007;88(5):1146–59.
Miller JR, Engstrom MD. The relationships of major lineages within peromyscine rodents: a molecular phylogenetic hypothesis and systematic reappraisal. J Mammal. 2008;89(5):1279–95.
Merriam CH. A new genus (Neotomodon) and three new species of murine rodents from the mountains of the southern Mexico. P Biol Soc Wash. 1898;12:127–9.
Williams SL, Ramirez-Palido J, Baker RJ. Peromyscus alstoni. Mammalian Species. 1985;242:1–4.
Bolger AM, Lohse, M., & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics: btu170; 2014.
Hahn C, Bachmann L, Chevreux B. Reconstructing mitochondrial genomes directly from genomic next-generation sequencing reads-a baiting and iterative mapping approach. Nucleic Acids Res. 2013;41(13):1–9.
Chevreux B, Wetter T, Suhai S. Genome sequence assembly using trace signals and additional sequence information. Hannover: Computer Science and Biology: Proceedings of the German Conference on Bioinformatics; 1999. p. 45–6.
Thomas O. New forms of Akodon and Phyllotis, and a new genus for “Akodon” teguina. Ann Mag Nat Hist. 1913;11:404–9.
Bernt M, Donath A, Juhling F, Externbrink F, Florentz C, Fritzsch G, Putz J, Middendorf M, Stadler PF. MITOS: improved de novo metazoan mitochondrial genome annotation. Mol Phylogenet Evol. 2013;69(2):313–9.
A.R Q. BEDTools: the Swiss-Army tool for genome feature analysis. Curr Protoc Bioinformatics. 2014;47(11):1–11.
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30(9):1312–3.
Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004;5:113.
Huelsenbeck JP, Ronquist F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001;17(8):754–5.
Say T, Ord G. Description of a new species of Mammalia, whereon a genus is proposed to be founded. J Acad Natl Sci Phila. 1825;4:352–6.
D'Elia G. Phylogenetics of sigmodontinae (Rodentia, Muroidea, Cricetidae), with special reference to the akodont group, and with additional comments on historical biogeography. Cladistics. 2003;19(4):307–23.
Teta P, Canon C, Patterson BD, Pardinas UFJ. Phylogeny of the tribe Abrotrichini (Cricetidae, Sigmodontinae): integrating morphological and molecular evidence into a new classification. Cladistics. 2017;33(2):153–82.
Merriam CH. Mammals of Idaho. Washington DC: North American fauna; 1891. p. 31–88.
Baker RJ, Bradley RD. Speciation in mammals and the genetic species concept. J Mammal. 2006;87(4):643–62.
Baird SF. Characteristics of some new species of north American Mammalia. Proc Acad Natl Sci Phila. 1855;7:333–7.
Ordonez-Garza N, Thompson CW, Unkefer MK, Edwards CW, Owen JG, Bradley RD. Systematics of the Neotoma mexicana species group (Mammalia: Rodentia: Cricetidae) in Mesoamerica: new molecular evidence on the status and relationships of N. Ferruginea tomes, 1862. P Biol Soc Wash. 2014;127(3):518–32.
Goldman E. New mammals from eastern Panama. Smithsonian Miscellaneous Collections. 1912;60(2):1–18.
Saussure MH. Note Sur Quelques Mamiferes de Mexique. Revue et Magasin de Zoologie. 1860;2(13):103–9.
FB Stangl RB. Evolutionary relationships in Peromyscus: congruence in chromosomal, genic, and classical data sets. J Mammal. 1984;65(4):643–54.
Goodwin GG. A new species and a new subspecies of Peromyscus from Oaxaca. Mexico American Museum Novitates. 1964;2183:1–8.
Osgood WH. Thirty new mice of the genus Peromyscus from Mexico and Guatemala. P Biol Soc Wash. 1904;17:55–77.
Allen JA. Descriptions of new north American mammals. Bull Am Mus Nat Hist. 1895;7:327–40.
Carleton MD, G. L. Kirkland, and J. N. Layne. Advances in the study of Peromyscus (Rodentia); 1989.
Hill RW. Metabolism, thermal conductance, and body-temperature in one of largest species of Peromyscus, Peromyscus pirrensis. J Therm Biol. 1976;1(2):109–12.
Tang QY, Liu SQ, Yu D, Liu HZ, Danley PD. Mitochondrial capture and incomplete lineage sorting in the diversification of balitorine loaches (Cypriniformes, Balitoridae) revealed by mitochondrial and nuclear genes. Zool Scr. 2012;41(3):233–47.
Sarver BA, Demboski JR, Good JM, Forshee N, Hunter SS, Sullivan J. Comparative phylogenomic assessment of mitochondrial introgression among several species of chipmunks (Tamias). Genome Biol Evol. 2017;9(1):7–19.
Good JM, Hird S, Reid N, Demboski JR, Steppan SJ, Martin-Nims TR, Sullivan J. Ancient hybridization and mitochondrial capture between two species of chipmunks. Mol Ecol. 2008;17(5):1313–27.
Li G, Davis BW, Eizirik E, Murphy WJ. Phylogenomic evidence for ancient hybridization in the genomes of living cats (Felidae). Genome Res. 2016;26(1):1–11.
Faircloth BC, McCormack JE, Crawford NG, Harvey MG, Brumfield RT, Glenn TC. Ultraconserved elements anchor thousands of genetic markers spanning multiple evolutionary timescales. Syst Biol. 2012;61(5):717–26.
Crawford NG, Faircloth BC, McCormack JE, Brumfield RT, Winker K, Glenn TC. More than 1000 ultraconserved elements provide evidence that turtles are the sister group of archosaurs. Biol Lett. 2012;8(5):783–6.
Ray DA, Xing J, Salem AH, Batzer MA. SINEs of a nearly perfect character. Syst Biol. 2006;55(6):928–35.
Platt RN 2nd, Zhang Y, Witherspoon DJ, Xing J, Suh A, Keith MS, Jorde LB, Stevens RD, Ray DA. Targeted capture of phylogenetically informative Ves SINE insertions in genus Myotis. Genome Biol Evol. 2015;7(6):1664–75.
The Natural Sciences Research Laboratory at the Museum of Texas Tech University kindly provided tissues via the Genetic Resources Collection. Additional thanks goes to RTL Genomics for sequencing efforts.
This work was supported by the National Science Foundation, DEB-1355176. Additional support was provided by College of Arts and Sciences at Texas Tech University.
Availability of data and materials
The raw reads have been uploaded to NCBI’s SRA via accession numbers SAMN04393360–04393374. A nexus file used for phylogenetic analysis is available through Dryad, doi:10.5061/dryad.c2j05.
Ethics approval and consent to participate
No ethical approval was required from an institutional or national ethics review board.
Consent for publication
The authors declare they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Title of data: Cytochrome b GenBank information. Taxa names and accession numbers for cytochrome b genes downloaded from GenBank in Additional file 3: Figure S1. (XLS 28 kb)
Whole mitochondrial exomes. Description: Alignment of the protein coding genes of each assembled mitogenome along with Akodon montensis, downloaded from GenBank. (TXT 257 kb)
Cytochrome b phylogenetic validation. ML tree of the assembled cytochrome b sequences compared with cytochrome b sequences downloaded from GenBank. These results showed if the sequences from our data correlated with known cytb sequences. (PNG 60 kb)
v2. Cytochrome b blast validation. Best hits to GenBank entries for cytochrome b genes from assembled mitochondrial genomes. For N. mexicana, the best hit was to a closely related species, with the second best hit to N. mexicana. (XLS 38 kb)