The Microtus voles: Resolving the phylogeny of one of the most speciose mammalian genera using genomics

Low Genetic Variability of the Tundra Vole in Lithuania

The distribution and spread of the tundra vole (Alexandromys oeconomus) in Lithuania have been documented over the last 70 years, but the genetic diversity of the species has not been studied. In this study, we examined A. oeconomus trapped in three sites in northern and western Lithuania using mtDNA sequence analysis of the cytb and control region. The western and northern sites are separated by anthropogenic landscape barriers. The western site is subject to regular spring flooding. Phylogenetic analyses of the studied individuals placed them in the Central European phylogroup, suggesting that Lithuanian A. oeconomus originated from northeastern Poland. In Lithuania, the genetic diversity of A. oeconomus at both mtDNA loci was relatively low (Hd < 0.6, π < 0.002) compared to that found in other European samples (Hd = 0.833-0.958; π = 0.00402-0.01552). Individuals analyzed in Lithuania were genetically different from samples collected in Poland and Northern Europe (ΦST > 0.15, p < 0.05). The genetic divergence between the western and northern samples of A. oeconomus in Lithuania, together with the low genetic variability among the voles studied, provides new insights into the phylogeography of the species and the influence of barriers on the colonization of the country.

A new perspective on the taxonomy and systematics of Arvicolinae (Gray, 1821) and a new time-calibrated phylogeny for the clade

Background

Arvicoline rodents are one of the most speciose and rapidly evolving mammalian lineages. Fossil arvicolines are also among the most common vertebrate fossils found in sites of Pliocene and Pleistocene age in Eurasia and North America. However, there is no taxonomically robust, well-supported, time-calibrated phylogeny for the group.

Methods

Here we present well-supported hypotheses of arvicoline rodent systematics using maximum likelihood and Bayesian inference of DNA sequences of two mitochondrial genes and three nuclear genes representing 146 (82% coverage) species and 100% of currently recognized arvicoline genera. We elucidate well-supported major clades, reviewed the relationships and taxonomy of many species and genera, and critically compared our resulting molecular phylogenetic hypotheses to previously published hypotheses. We also used five fossil calibrations to generate a time-calibrated phylogeny of Arvicolinae that permitted some reconciliation between paleontological and neontological data.

Results

Our results are largely congruent with previous molecular phylogenies, but we increased the support in many regions of the arvicoline tree that were previously poorly-sampled. Our sampling resulted in a better understanding of relationships within Clethrionomyini, the early-diverging position and close relationship of true lemmings (Lemmus and Myopus) and bog lemmings (Synaptomys), and provided support for recent taxonomic changes within Microtini. Our results indicate an origin of ∼6.4 Ma for crown arvicoline rodents. These results have major implications (e.g., diversification rates, paleobiogeography) for our confidence in the fossil record of arvicolines and their utility as biochronological tools in Eurasia and North America during the Quaternary.

Novel frontier in wildlife monitoring: Identification of small rodent species from fecal pellets using near-infrared reflectance spectroscopy (NIRS)

Small rodents are prevalent and functionally important across the world's biomes, making their monitoring salient for ecosystem management, conservation, forestry, and agriculture. There is a growing need for cost-effective and noninvasive methods for large-scale, intensive sampling. Fecal pellet counts readily provide relative abundance indices, and given suitable analytical methods, feces could also allow for the determination of multiple ecological and physiological variables, including community composition. In this context, we developed calibration models for rodent taxonomic determination using fecal near-infrared reflectance spectroscopy (fNIRS). Our results demonstrate fNIRS as an accurate and robust method for predicting genus and species identity of five coexisting subarctic microtine rodent species. We show that sample exposure to weathering increases the method's accuracy, indicating its suitability for samples collected from the field. Diet was not a major determinant of species prediction accuracy in our samples, as diet exhibited large variation and overlap between species. fNIRS could also be applied across regions, as calibration models including samples from two regions provided a good prediction accuracy for both regions. We show fNIRS as a fast and cost-efficient high-throughput method for rodent taxonomic determination, with the potential for cross-regional calibrations and the use on field-collected samples. Importantly, appeal lies in the versatility of fNIRS. In addition to rodent population censuses, fNIRS can provide information on demography, fecal nutrients, stress hormones, and even disease. Given the development of such calibration models, fNIRS analytics could complement novel genetic methods and greatly support ecosystem- and interaction-based approaches to monitoring.

Complete mitogenomes of Turkish tree squirrels, Sciurus anomalus and S. vulgaris, (Sciuridae: Rodentia: Mammalia) and their phylogenetic status within the tribe Sciurini

The genus Sciurus, a member of the family Sciuridae, is widely distributed in the Holarctic region. To better understand mitogenomic characteristics and to reveal internal phylogenetic relationships of the genus, 20 complete mitogenomes of Turkish tree squirrels were successfully sequenced for the first time, including 19 for S. anomalus (from 16,505 bp to 16,510 bp) and one for S. vulgaris (16,511 bp). The mitogenomes of two species were AT-biased. All tRNAs for two species displayed a typical clover-leaf structure, except for tRNA^Ser(AGY). The tRNA Serine1 (S1)-GCT structure lacked the dihydrouridine (DHU) loop and stem. Based on mitogenomic dataset for phylogeny of Sciurinae, phylogenetic analyses (Bayesian Inference and Maximum Likelihood) did not support monophyly of Sciurus and proposed that S. anomalus, the most basal taxa in the Sciurini tribe, had at least five mitogenome lineages, which were also supported by network analysis. The dissimilarities among the five linegaes of S. anomalus ranged from 0.0042 (0.42%) to 0.0062 (0.62%) using K2P sequence pairwise distances. In addition to this mitogenomic analysis result, phylogenetic analyses using the CYTB + D-loop dataset proposed the existence of at least nine lineages for S. anomalus, which was different than those of the previous studies. The current study proposed that the use of mitogenomic data for reconstructing the phylogeny of Turkey' Sciurus holds an important value for revealing evolutionary relationships.

A mitochondrial genome phylogeny of voles and lemmings (Rodentia: Arvicolinae): Evolutionary and taxonomic implications

Arvicolinae is one of the most impressive placental radiations with over 150 extant and numerous extinct species that emerged since the Miocene in the Northern Hemisphere. The phylogeny of Arvicolinae has been studied intensively for several decades using morphological and genetic methods. Here, we sequenced 30 new mitochondrial genomes to better understand the evolutionary relationships among the major tribes and genera within the subfamily. The phylogenetic and molecular dating analyses based on 11,391 bp concatenated alignment of protein-coding mitochondrial genes confirmed the monophyly of the subfamily. While Bayesian analysis provided a high resolution across the entire tree, Maximum Likelihood tree reconstruction showed weak support for the ordering of divergence and interrelationships of tribal level taxa within the most ancient radiation. Both the interrelationships among tribes Lagurini, Ellobiusini and Arvicolini, comprising the largest radiation and the position of the genus Dinaromys within it also remained unresolved. For the first time complex relationships between genus level taxa within the species-rich tribe Arvicolini received full resolution. Particularly Lemmiscus was robustly placed as sister to the snow voles Chionomys in the tribe Arvicolini in contrast with a long-held belief of its affinity with Lagurini. Molecular dating of the origin of Arvicolinae and early divergences obtained from the mitogenome data were consistent with fossil records. The mtDNA estimates for putative ancestors of the most genera within Arvicolini appeared to be much older than it was previously proposed in paleontological studies.

The Tien Shan vole (Microtus ilaeus; Rodentia: Cricetidae) as a new species in the Late Pleistocene of Europe

Grey voles (subgenus Microtus) represent a complex of at least seven closely related and partly cryptic species. The range of these species extends from the Atlantic to the Altai Mountains, but most of them occur east of the Black Sea. Using ancient DNA analyses of the Late Pleistocene specimens, we identified a new mtDNA lineage of grey voles in Europe. Phylogenetic analysis of mitochondrial DNA cytochrome b sequences from 23 voles from three caves, namely, Emine-Bair-Khosar (Crimea, Ukraine), Cave 16 (Bulgaria), and Bacho Kiro (Bulgaria), showed that 14 specimens form a previously unrecognized lineage, sister to the Tien Shan vole. The average sequence divergence of this lineage and the extant Tien Shan vole was 4.8%, which is similar to the divergence of grey vole forms, which are considered distinct species or being on the verge of speciation; M. arvalis and M. obscurus or M. mystacinus and M. rossiaemeridionalis. We estimated the time to the most recent common ancestor of the grey voles to be 0.66 Ma, which is over twice the recent estimates, while the divergence of the extant Tien Shan vole and the new lineage to be 0.29 Ma. Our discovery suggests that grey voles may have been more diversified in the past and that their ranges may have differed substantially from current ones. It also underlines the utility of ancient DNA to decipher the evolutionary history of voles.

Genomic Resources for the North American Water Vole (Microtus richardsoni) and the Montane Vole (Microtus montanus)

Voles of the genus Microtus are important research organisms, yet genomic resources are lacking. Such resources would benefit future studies of immunology, phylogeography, cryptic diversity, and more. We sequenced and assembled nuclear genomes from two subspecies of water vole (Microtus richardsoni) and from the montane vole (Microtus montanus). The water vole genomes were sequenced with Illumina and 10× Chromium plus Illumina sequencing, resulting in assemblies with ∼1600,000 and ∼30,000 scaffolds, respectively. The montane vole was also assembled into ∼13,000 scaffolds using Illumina sequencing. Mitochondrial genome assemblies were also performed for both species. Structural and functional annotation for the best water vole nuclear genome resulted in ∼24,500 annotated genes, with 83% of these having functional annotations. Assembly quality statistics for our nuclear assemblies fall within the range of genomes previously published in the genus Microtus, making the water vole and montane vole genomes useful additions to currently available genomic resources.

Mitochondrial Genome Evolution, Genetic Diversity, and Population Structure in British Water Voles (Arvicola amphibius)

The European water vole (Arvicola amphibius) is a rodent within the subfamily Arvicolinae. In Britain, water voles have declined rapidly during the last century, making them a conservation priority. The relationship of Arvicola to other genera within Arvicolinae remains debated. Additionally, the impact that captive breeding programs in Britain are having on the genetic diversity of water voles is unknown. We use available mitochondrial genomes to construct the phylogeny of species within Arvicolinae, followed by sequencing the mitochondrial DNA control region of 17 individuals from a captive population of water voles in Britain to assess their genetic diversity and population structure. Our study first provides an updated phylogenetic tree of Arvicolinae using the mitochondrial genome of 31 species. Second, our results show considerable genetic diversity in the captive population of water voles, when compared with natural populations in Britain. We confirm the grouping of British water voles into two clades, with all captive individuals found in the English/Welsh clade. Moreover, captive water voles clustered closely with populations in the South East and East of England. The mitochondrial genome provides a useful marker to study the phylogenetics of this rodent clade and in addition, our study provides support for the breeding program at Wildwood Trust and provides a framework for future conservation genetics studies in this species.

Complete Mitochondrial Genome of Three Species of the Genus Microtus (Arvicolinae, Rodentia)

The 65 species of the genus Microtus have unusual sex-related genetic features and a high rate of karyotype variation. However, only nine complete mitogenomes for these species are currently available. We describe the complete mitogenome sequences of three Microtus, which vary in length from 16,295 bp to 16,331 bp, contain 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes and a control region. The length of the 13 PCGs and the coded proteins is the same in all three species, and the start and stop codons are conserved. The non-coding regions include the L-strand origin of replication, with the same sequence of 35 bp, and the control region, which varies between 896 bp and 930 bp in length. The control region includes three domains (Domains I, II and III) with extended termination-associated sequences (ETAS-1 and ETAS-2) in Domain I. Domain II and Domain III include five (CSB-B, C, D, E and F) and three (CSB-1, CSB-2, and CSB-3) conserved sequence blocks, respectively. Phylogenetic reconstructions using the mitochondrial genomes of all the available Microtus species and one representative species from another genus of the Arvicolinae subfamily reproduced the established phylogenetic relationships for all the Arvicolinae genera that were analyzed.

Complete mitochondrial genomes confirm the generic placement of the plateau vole, Neodon fuscus

The plateau vole, Neodon fuscus is endemic to China and is distributed mainly in Qinghai Province. It is of public health interest, as it is, a potential reservoir of Toxoplasma gondii and the intermediate host of Echinococcus multilocularis. However, genetic data of this species are lacking, and its name and taxonomy are still a controversy. In the present study, we determined the nucleotide sequence of the entire mitochondrial (mt) genome of N. fuscus and analyzed its evolutionary relationship. The mitogenome was 16328 bp in length and contained 13 protein-coding genes, 22 genes for transfer RNAs (tRNA), two ribosomal RNA genes and two major noncoding regions (O_L region and D-loop region). Most genes were located on the heavy strand. All tRNA genes had typical cloverleaf structures except for tRNA^{Ser (GCU)}. The mt genome of N. fuscus was rich in A+T (58.45%). Maximum likelihood (ML) and Bayesian methods yielded phylogenetic trees from 33 mt genomes of Arvicolinae, in which N. fuscus formed a sister group with Neodon irene and Neodon sikimensis to the exclusion of species of Microtus and other members of the Arvicolinae. Further phylogenetic analyses (ML only) based on the cytb gene sequences also demonstrated that N. fuscus had a close relationship with N. irene. The complete mitochondrial genome was successfully assembled and annotated, providing the necessary information for the phylogenetic analyses. Although the name Lasiopodomys fuscus was used in the book ‘Wilson & Reeder’s Mammal Species of the World’, we have confirmed here that its appropriate name is N. fuscus through an analysis of the evolutionary relationships.

Multiple intrasyntenic rearrangements and rapid speciation in voles

Remarkably stable genomic chromosome elements (evolutionary conserved segments or syntenies) are the basis of large-scale chromosome architecture in vertebrate species. However, these syntenic elements harbour evolutionary important changes through intrachromosomal rearrangements such as inversions and centromere repositioning. Here, using FISH with a set of 20 region-specific probes on a wide array of 28 species, we analyzed evolution of three conserved syntenic regions of the Arvicolinae ancestral karyotype. Inside these syntenies we uncovered multiple, previously cryptic intrachromosomal rearrangements. Although in each of the three conserved blocks we found inversions and centromere repositions, the blocks experienced different types of rearrangements. In two syntenies centromere repositioning predominated, while in the third region, paracentric inversions were more frequent, whereas pericentric inversions were not detected. We found that some of the intrachromosomal rearrangements, mainly paracentric inversions, were synapomorphic for whole arvicoline genera or tribes: genera Alexandromys and Microtus, tribes Ellobini and Myodini. We hypothesize that intrachromosomal rearrangements within conserved syntenic blocks are a major evolutionary force modulating genome architecture in species-rich and rapidly-evolving rodent taxa. Inversions and centromere repositioning may impact speciation and provide a potential link between genome evolution, speciation, and biogeography.