Population connectivity in voles (Microtus sp.) as a gauge for tall grass prairie restoration in midwestern North America.

Ecological restoration can promote biodiversity conservation in anthropogenically fragmented habitats, but effectiveness of these management efforts need to be statistically validated to determine 'success.' One such approach is to gauge the extent of recolonization as a measure of landscape permeability and, in turn, population connectivity. In this context, we estimated dispersal and population connectivity in prairie vole (Microtus ochrogaster; N = 231) and meadow vole (M. pennsylvanicus; N = 83) within five tall-grass prairie restoration sites embedded within the agricultural matrix of midwestern North America. We predicted that vole dispersal would be constrained by the extent of agricultural land surrounding restored habitat patches, spatially isolating vole populations and resulting in significant genetic structure. We first employed genetic assignment tests based on 15 microsatellite DNA loci to validate field-derived species-designations, then tested reclassified samples with multivariate and Bayesian clustering to assay for spatial and temporal genetic structure. Population connectivity was further evaluated by calculating pairwise FST, then potential demographic effects explored by computing migration rates, effective population size (Ne), and average relatedness (r). Genetic species assignments reclassified 25% of initial field identifications (N = 11 M. ochrogaster; N = 67 M. pennsylvanicus). In M. ochrogaster population connectivity was high across the study area, reflected in little to no spatial or temporal genetic structure. In M. pennsylvanicus genetic structure was detected, but relatedness estimates identified it as kin-clustering instead, underscoring social behavior among populations rather than spatial isolation as the cause. Estimates of Ne and r were stable across years, reflecting high dispersal and demographic resilience. Combined, these metrics suggest the agricultural matrix is highly permeable for voles and does not impede dispersal. High connectivity observed confirms that the restored landscape is productive and permeable for specific management targets such as voles and also demonstrates population genetic assays as a tool to statistically evaluate effectiveness of conservation initiatives.

Timing outweighs magnitude of rainfall in shaping population dynamics of a small mammal species in steppe grassland.

Climate change-induced shifts in species phenology differ widely across trophic levels, which may lead to consumer-resource mismatches with cascading population and ecosystem consequences. Here, we examined the effects of different rainfall patterns (i.e., timing and amount) on the phenological asynchrony of population of a generalist herbivore and their food sources in semiarid steppe grassland in Inner Mongolia. We conducted a 10-y (2010 to 2019) rainfall manipulation experiment in 12 0.48-ha field enclosures and found that moderate rainfall increases during the early rather than late growing season advanced the timing of peak reproduction and drove marked increases in population size through increasing the biomass of preferred plant species. By contrast, greatly increased rainfall produced no further increases in vole population growth due to the potential negative effect of the flooding of burrows. The increases in vole population size were more coupled with increased reproduction of overwintered voles and increased body mass of young-of-year than with better survival. Our results provide experimental evidence for the fitness consequences of phenological mismatches at the population level and highlight the importance of rainfall timing on the population dynamics of small herbivores in the steppe grassland environment.

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.

Whole mitochondrial genome of long-clawed mole vole (Prometheomys schaposchnikowi) from Turkey, with its phylogenetic relationships.

The mitogenome of Prometheomys schaposchnikowi was characterized for the first time as a circular DNA molecule (16.284 bp), containing 37 coding and 2 non-coding regions. In the mitogenome, ND6 and 8 tRNA genes were encoded on the light chain, while 12 PCGs, 14 tRNAs, 2 rRNAs, D-loop and OL were encoded on the heavy chain. The most common initiation codon in PCGs was ATG. As in many mammals, incomplete stop codons in P. schaposchnikowi were in the COX3, ND1 and ND4. Phylogenetic relationships were revealed using Bayesian method and the 13 PCGs. Seven genera (Arvicola, Dicrostonyx, Lasiopodomys, Myodes, Ondatra, Proedromys and Prometheomys) formed a monophyletic group, while Eothenomys, Microtus and Neodon were paraphyletic. P. schaposchnikowi constituted the most basal group within Arvicolinae. Divergence time estimation suggested that P. schaposchnikowi diversified during the Miocene (16.28 Mya). Further molecular studies are needed to test the distinctiveness and diversity of the genus Prometheomys.

Complex Structure of Lasiopodomys mandarinus vinogradovi Sex Chromosomes, Sex Determination, and Intraspecific Autosomal Polymorphism.

The mandarin vole, Lasiopodomys mandarinus, is one of the most intriguing species among mammals with non-XX/XY sex chromosome system. It combines polymorphism in diploid chromosome numbers, variation in the morphology of autosomes, heteromorphism of X chromosomes, and several sex chromosome systems the origin of which remains unexplained. Here we elucidate the sex determination system in Lasiopodomys mandarinus vinogradovi using extensive karyotyping, crossbreeding experiments, molecular cytogenetic methods, and single chromosome DNA sequencing. Among 205 karyotyped voles, one male and three female combinations of sex chromosomes were revealed. The chromosome segregation pattern and karyomorph-related reproductive performances suggested an aberrant sex determination with almost half of the females carrying neo-X/neo-Y combination. The comparative chromosome painting strongly supported this proposition and revealed the mandarin vole sex chromosome systems originated due to at least two de novo autosomal translocations onto the ancestral X chromosome. The polymorphism in autosome 2 was not related to sex chromosome variability and was proved to result from pericentric inversions. Sequencing of microdissection derived of sex chromosomes allowed the determination of the coordinates for syntenic regions but did not reveal any Y-specific sequences. Several possible sex determination mechanisms as well as interpopulation karyological differences are discussed.

Highly divergent lineage of narrow-headed vole from the Late Pleistocene Europe.

During the Late Pleistocene, narrow-headed voles (Lasiopodomys gregalis) inhabited Eurasia's vast territories, frequently becoming the dominant small mammal species among steppe-tundra communities. We investigated the relationship between this species' European and Asiatic populations by sequencing the mtDNA genomes of two extant specimens from Russia and 10 individuals from five Central European sites, dated to the post-LGM period. Phylogenetic analyses based on a large portion of mtDNA genomes highly supported the positioning of L. gregalis within Arvicolinae. The phylogeny based on mtDNA cytochrome b sequences revealed a deep divergence of European narrow-headed voles from Asiatic ones and their sister position against the extant L. gregalis and L. raddei. The divergence of the European lineage was estimated to a minimum 230 thousand years ago. This suggest, contrary to the current biogeographic hypotheses, that during the interglacial periods narrow-headed vole did not retreat from Europe but survived the unfavourable conditions within the refugial areas. Based on this result, we propose to establish a cryptic species status for the Late Pleistocene European narrow-headed vole and to name this taxon Lasiopodomys anglicus.

Geographical distribution and hosts of the cestode Paranoplocephala omphalodes (Hermann, 1783) Lühe, 1910 in Russia and adjacent territories.

Paranoplocephala omphalodes is a widespread parasite of voles. Low morphological variability within the genus Paranoplocephala has led to erroneous identification of P. omphalodes a wide range of definitive hosts. The use of molecular methods in the earlier investigations has confirmed that P. omphalodes parasitizes four vole species in Europe. We studied the distribution of P. omphalodes in Russia and Kazakhstan using molecular tools. The study of 3248 individuals of 20 arvicoline species confirmed a wide distribution of P. omphalodes. Cestodes of this species were found in Microtus arvalis, M. levis, M. agrestis, Arvicola amphibius, and also in Chionomys gud. Analysis of the mitochondrial gene cox1 variability revealed a low haplotype diversity in P. omphalodes in Eurasia.

Subpopulation augmentation among habitat patches as a tool to manage an endangered Mojave Desert wetlands-dependent rodent during anthropogenic restricted water climate regimes.

Intensive management may be necessary to protect some highly vulnerable endangered species, particularly those dependent on water availability regimes that might be disrupted by ongoing climate change. The Amargosa vole (Microtus californicus scirpensis) is an increasingly imperiled rodent constrained to rare wetland habitat in the Mojave Desert. In 2014 and 2016, we trapped and radio-collared 30 voles, 24 were translocated and six remained at donor and recipient marshes as resident control voles. Soft-release was performed followed by remote camera and radio-telemetry monitoring. Although comparative metrics were not statistically significant, the mean maximum known distance moved (MDM) was longer for translocated (82.3 ± 14.6 m) vs. resident-control voles (74.9 ± 17.5 m) and for female (98.4 ± 19.9 m) vs. male (57.8 ± 9.1 m) voles. The mean area occupied (AO) tended to be greater in female (0.15 ± 0.04 ha) vs. male (0.12 ± 0.03 ha) voles, and control voles (0.15 ± 0.05 ha) compared with translocated voles (0.13 ± 0.03 ha). The mean minimum known time alive (MTA) was 38.2 ± 19.4 days for resident-control voles and 47.0 ± 10.6 days for translocated voles. Female survival (55.7 ± 14.3 days) exceeded that of males (31.5 ± 9.4 days) regardless of study group. Activity in bulrush/rushes mix and bulrush vegetation types was strongly and significantly overrepresented compared with salt grass and rushes alone, and habitat selection did not differ between resident and translocated voles. Our results provide ecological and methodological insights for future translocations as part of a strategy of promoting long-term survival of an extremely endangered small mammal in a wild desert environment.

Chromosome Translocations as a Driver of Diversification in Mole Voles Ellobius (Rodentia, Mammalia).

The involvement of chromosome changes in the initial steps of speciation is controversial. Here we examine diversification trends within the mole voles Ellobius, a group of subterranean rodents. The first description of their chromosome variability was published almost 40 years ago. Studying the G-band structure of chromosomes in numerous individuals revealed subsequent homologous, step-by-step, Robertsonian translocations, which changed diploid numbers from 54 to 30. Here we used a molecular cytogenetic strategy which demonstrates that chromosomal translocations are not always homologous; consequently, karyotypes with the same diploid number can carry different combinations of metacentrics. We further showed that at least three chromosomal forms with 2n = 34 and distinct metacentrics inhabit the Pamir-Alay mountains. Each of these forms independently hybridized with E. tancrei, 2n = 54, forming separate hybrid zones. The chromosomal variations correlate slightly with geographic barriers. Additionally, we confirmed that the emergence of partial or monobrachial homology appeared to be a strong barrier for hybridization in nature, in contradistinction to experiments which we reported earlier. We discuss the possibility of whole arm reciprocal translocations for mole voles. Our findings suggest that chromosomal translocations lead to diversification and speciation.

Further assessment of the Genus Neodon and the description of a new species from Nepal.

Recent molecular systematic studies of arvicoline voles of the genera Neodon, Lasiopodomys, Phaiomys, and Microtus from Central Asia suggest the inclusion of Phaiomys leucurus, Microtus clarkei, and Lasiopodomys fuscus into Neodon and moving Neodon juldaschi into Microtus (Blanfordimys). In addition, three new species of Neodon (N. linzhiensis, N. medogensis, and N. nyalamensis) have recently been described from Tibet. Analyses of concatenated mitochondrial (Cytb, COI) and nuclear (Ghr, Rbp3) genes recovered Neodon as a well-supported monophyletic clade including all the recently described and relocated species. Kimura-2-parameter distance between Neodon from western Nepal compared to N. sikimensis (K2P = 13.1) and N. irene (K2P = 13.4) was equivalent to genetic distances observed between recognized species of this genus. The specimens sampled from western Nepal were recovered sister to N. sikimensis in the concatenated analysis. However, analyses conducted exclusively with mitochondrial loci did not support this relationship. The occlusal patterns of the first lower (m1) and third upper (M3) molars were simpler in specimens from western Nepal in comparison to N. sikimensis from eastern Nepal and India. Twelve craniodental characters and four external field measurements were examined from specimens of N. sikimensis from eastern Nepal and India, N. irene, and Neodon from western Nepal. Neodon from western Nepal were significantly different from N. sikimensis from eastern Nepal and India in ten out of 16 characters measured and from N. irene for all characters except ear height. Specimens from western Nepal were smaller in size than N. sikimensis from Eastern Nepal and India and larger than N. irene. Together the results of the molecular and morphological analyses indicate that Neodon from western Nepal are distinct under the phylogenetic, genetic and morpho species concepts.

Molecular evidence of Tula virus in Microtus obscurus in the region of Yili, Xinjiang, China.

BACKGROUND: Hantaviruses are important zoonotic pathogens, and they pose a profound risk to public health. So far, there has been no evidence showing that Tula virus (TULV), one species of hantavirus, is endemic in China. In this study, we captured rodents and found that the Tula virus had infected voles in Yili region, Xinjiang, China. METHODS: Rodents were captured by flooding their burrows in mountain pasture areas in Narati, Xinyuan County, Xinjiang, China. Hantavirus L gene fragments were amplified by nest RT-PCR using genus-specific primers. Positive samples were further identified by sequencing of RT-PCR products of S gene fragment for species identification. To identify the species of captured small mammals, the rodents' cytochrome b (Cytb) was amplified by PCR and sequenced. Phylogenetic analysis was used to show the clustering and evolution relationship of the viral nucleic acids. RESULTS: Here, 31 out of 198 voles captured (16%) were infected with TULV. Host sequencing analysis showed these voles were Microtus obscurus (M. obscurs). Alignment and phylogenetic analysis of the exon region (1191 bp) of the hantavirus S gene confirmed that all of the detected amplicons were TULV, which was similar to one strain of TULV identified in Kazakhstan. CONCLUSION: This is the first identification of Tula virus in China, and we found that M. obscurus acts as a natural reservoir for carrying the virus. Although the infection rate in the local human population remains unknown, the high prevalence of TULV in the small mammals in the region constitutes a risk that this putative pathogen may spread to the local population.

Spatial incongruence in the species richness and functional diversity of cricetid rodents.

Biodiversity is multidimensional and different mechanisms can influence different dimensions. The spatial distribution of these dimensions can help in conservation decisions through the location of complementary areas with high diversity. We analyzed congruence in spatial patterns of species richness and functional diversity of cricetid rodents in the state of Oaxaca, southern Mexico, at different scales, and environmental variables related. Potential distribution models were produced for 49 species of cricetids in Maxent and superimposed to obtain potential communities in cells of 25, 50,100, 200 and 400 km2. We estimated species richness (SR) and functional diversity (SES.FD) eliminating the species richness effect through null models. The patterns and spatial congruence of species richness and functional diversity are described. The relationships between the environmental variables (elevation, temperature, precipitation, net primary productivity and potential evapotranspiration) and the SR and SES.FD were explored using Generalized Linear Models (GLMs) and Generalized Additive Models (GAMs). The highest species richness was found in mountainous ecosystems while the highest functional diversity was in tropical forests, revealing a spatial incongruence among these components of biodiversity (r = -0.14, p = 0.42; Pearson correlation). The locations of the cells of low congruence varied according to spatial resolution. In univariate models, elevation was the variable that best explained species richness (R2 = 0.77). No single variable explained the functional diversity; however, the models that included multiple environmental variables partially explained both the high and low functional diversity. The different patterns suggest that different historic, ecological and environmental processes could be responsible for the community structure of cricetid rodents in Oaxaca. These results indicate that one great challenge to be met to achieve more effective planning for biological conservation is to integrate knowledge regarding the spatial distribution of different dimensions of biodiversity.

Secondary contact between diverged host lineages entails ecological speciation in a European hantavirus.

The diversity of viruses probably exceeds biodiversity of eukaryotes, but little is known about the origin and emergence of novel virus species. Experimentation and disease outbreak investigations have allowed the characterization of rapid molecular virus adaptation. However, the processes leading to the establishment of functionally distinct virus taxa in nature remain obscure. Here, we demonstrate that incipient speciation in a natural host species has generated distinct ecological niches leading to adaptive isolation in an RNA virus. We found a very strong association between the distributions of two major phylogenetic clades in Tula orthohantavirus (TULV) and the rodent host lineages in a natural hybrid zone of the European common vole (Microtus arvalis). The spatial transition between the virus clades in replicated geographic clines is at least eight times narrower than between the hybridizing host lineages. This suggests a strong barrier for effective virus transmission despite frequent dispersal and gene flow among local host populations, and translates to a complete turnover of the adaptive background of TULV within a few hundred meters in the open, unobstructed landscape. Genetic differences between TULV clades are homogenously distributed in the genomes and mostly synonymous (93.1%), except for a cluster of nonsynonymous changes in the 5' region of the viral envelope glycoprotein gene, potentially involved in host-driven isolation. Evolutionary relationships between TULV clades indicate an emergence of these viruses through rapid differential adaptation to the previously diverged host lineages that resulted in levels of ecological isolation exceeding the progress of speciation in their vertebrate hosts.

Diet shift in bank voles induced by competition from grey-sided voles?

Grey-sided voles (Myodes rufocanus) and bank voles (Myodes glareolus) co-exist in boreal forests in northern Scandinavia. Previous studies suggest that the 2 species interact interspecifically, the grey-sided vole being the dominant species. We tested the hypothesis that bank voles shift their diet due to competition with the dominant grey-sided vole by studying stable isotope ratios in both species. Muscle samples were taken from voles in patches of old forest occupied by only bank voles and patches of old forest occupied by both grey-sided voles and bank voles. We found that: (i) stable isotope ratios of bank voles differed in areas with and without grey-sided voles; and that (ii) the stable isotope ratios of bank voles were more similar to those of grey-sided voles in areas where grey-sided voles were absent. Our data suggests that grey-sided voles forced bank voles to change their diet due to interspecific competition.

Taxonomic Interpretation of Allopatric Mammalian Forms on the Example of Two Karyoforms of Microtus (Terricola) subterraneus (Rodentia, Arvicolinae) from Eastern Europe.

New European pine vole records from the Novgorod, Kaluga, Voronezh, and Belgorod oblasts were studied by sequencing of the mtDNA cytb gene (1143 bp) and by karyotyping (routine staining and G-banding techniques). The results enabled us to summarize chromosome variability of this species throughout Eastern Europe. In the sample studied, two geographically replacing chromosomal forms have been identified: northern, 2n = 54 (Novgorod and Kaluga oblasts), and southern, 2n = 52 (Voronezh and, presumably, Belgorod oblasts). Our data make the boundaries of these two karyoforms in Eastern Europe more precise and testify to intraspecific level of their taxonomic differentiation.

Molecular systematics and phylogeography of the endemic Osgood's deermouse Osgoodomys banderanus (Rodentia: Cricetidae) in the lowlands of western Mexico.

Osgoodomys banderanus is a recognized and endemic rodent species of western Mexico, an area known for its high biodiversity and number of endemisms. Phylogeographical relationships within this taxon were analyzed based on mitochondrial (ND3, tRNA-Arginine, ND4L and partial ND4) and nuclear (GHR) nucleotide sequences. We obtained a total of 112 samples from 22 localities, covering the complete distribution of the species. Phylogenetic analyses using Maximum Likelihood and Bayesian inference confirmed that Osgoodomys is a monophyletic group. In addition, phylogenetic and phylogeographic analyses detected three major clades, which do not coincide with the recognized subspecies of O. banderanus. The genetic lineages detected are the western clade (Nayarit, Jalisco and northern Colima), the central clade (Colima, Michoacán, and northern Guerrero) and the eastern clade (central and southern Guerrero). Genetic distances among clades (5-9%) and nucleotide substitutions (30-88) among haplogroups were high, especially in the southern group. Mountain ranges such as the Transmexican Volcanic Belt and the Sierra Madre del Sur, as well as the Balsas River act as geographical barriers for Osgoodomys. Our results suggest the presence of three independent species, which need to be characterized morphologically to confirm our hypothesis.

Effects of the Mitochondrial and Nuclear Genomes on Nonshivering Thermogenesis in a Wild Derived Rodent.

A key adaptation of mammals to their environment is their ability to maintain a constant high body temperature, even at rest, under a wide range of ambient temperatures. In cold climates, this is achieved by an adaptive production of endogenous heat, known as nonshivering thermogenesis (NST), in the brown adipose tissue (BAT). This organ, unique to mammals, contains a very high density of mitochondria, and BAT correct functioning relies on the correct functioning of its mitochondria. Mitochondria enclose proteins encoded both in the maternally inherited mitochondrial genome and in the biparentally inherited nuclear genome, and one overlooked hypothesis is that both genomes and their interaction may shape NST. By housing under standardized conditions wild-derived common voles (Microtus arvalis) from two distinct evolutionary lineages (Western [W] and Central [C]), we show that W voles had greater NST than C voles. By introgressing those two lineages over at least nine generations, we then experimentally tested the influence of the nuclear and mitochondrial genomes on NST and related phenotypic traits. We found that between-lineage variation in NST and BAT size were significantly influenced by the mitochondrial and nuclear genomes, respectively, with the W mitochondrial genotype being associated with higher NST and the W nuclear genotype with a larger BAT. There were significant mito-nuclear interactions on whole animal body weight and resting metabolic rate (RMR). Hybrid voles were lighter and had higher RMR. Overall, our findings turn new light on the influence of the mitochondrial and nuclear genomes on thermogenesis and building adaptation to the environment in mammals.

Lasiopodomys fuscus as an important intermediate host for Echinococcus multilocularis: isolation and phylogenetic identification of the parasite.

BACKGROUND: Echinococcus multilocularis causes alveolar echinococcosis (AE) and is widely prevalent in Qinghai Province, China, where a number of different species have been identified as hosts. However, limited information is available on the Qinghai vole (Lasiopodomys fuscus), which is hyper endemic to Qinghai Province and may represent a potential intermediate host of E. multilocularis. Thus, L. fuscus could contribute to the endemicity of AE in the area. METHODS: Fifty Qinghai voles were captured from Jigzhi County in Qinghai Province for the clinical identification of E. multilocularis infection via anatomical examination. Hydatid fluid was collected from vesicles of the livers in suspected voles and subjected to a microscopic examination and PCR assay based on the barcoding gene of cox 1. PCR-amplified segments were sequenced for a phylogenetic analysis. E. multilocularis-infected Qinghai voles were morphologically identified and subjected to a phylogenetic analysis to confirm their identities. RESULTS: Seventeen of the 50 Qinghai voles had E. multilocularis-infection-like vesicles in their livers. Eleven out of the 17 Qinghai voles presented E. multilocularis infection, which was detected by PCR and sequencing. The phylogenetic analysis showed that all 11 positive samples belonged to the E. multilocularis Asian genotype. A morphological identification and phylogenetic analysis of the E. multilocularis-infected Qinghai voles confirmed that all captured animals were L. fuscus. CONCLUSIONS: L. fuscus can be infected with E. multilocularis and plays a potential role in the life cycle and epidemiology of E. multilocularis in the Qinghai-Tibetan Plateau of China.

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

Sequential rapid radiations pose some of the greatest difficulties in phylogenetics, especially when analysing only a small number of genetic markers. Given that most of the speciation events occur in quick succession at various points in time, this creates particular challenges in determining phylogenetic relationships, i.e. branching order and divergence times. With the development of high throughput sequencing, thousands of markers can now readily be used to tackle these issues. Microtus is a speciose genus currently composed of 65 species that evolved over the last 2 million years. Although it is a well-studied group, there is still phylogenetic uncertainty at various divergence levels. Building upon previous studies that generally used small numbers of mitochondrial and/or nuclear loci, in this genomic-scale study we used both mitochondrial and nuclear data to study the rapid radiation within Microtus, using partial mitogenomes and genotyping-by-sequencing (GBS) on seven species representing five Microtus subgenera and the main biogeographic ranges where this group occurs. Both types of genome (mitochondrial and nuclear) generated similar tree topologies, with a basal split of the Nearctic (M. ochrogaster) and Holarctic (M. oeconomus) species, and then a subdivision of the five Palearctic species into two subgroups. These data support the occurrence of two European radiations, one North American radiation, and a later expansion of M. oeconomus from Asia to both Europe and North America. We further resolved the positioning of M. cabrerae as sister group of M. agrestis and refute the claim that M. cabrerae should be elevated to its own genus (Iberomys). Finally, the data support ongoing speciation events, especially within M. agrestis, with high levels of genetic divergence between the three Evolutionarily Significant Units (ESUs) previously identified. Similar high levels of divergence were also found among ESUs within M. oeconomus and M. arvalis.

Impact of tree priors in species delimitation and phylogenetics of the genus Oligoryzomys (Rodentia: Cricetidae).

The use of genetic data and tree-based algorithms to delimit evolutionary lineages is becoming an important practice in taxonomic identification, especially in morphologically cryptic groups. The effects of different phylogenetic and/or coalescent models in the analyses of species delimitation, however, are not clear. In this paper, we assess the impact of different evolutionary priors in phylogenetic estimation, species delimitation, and molecular dating of the genus Oligoryzomys (Mammalia: Rodentia), a group with complex taxonomy and morphological cryptic species. Phylogenetic and coalescent analyses included 20 of the 24 recognized species of the genus, comprising of 416 Cytochrome b sequences, 26 Cytochrome c oxidase I sequences, and 27 Beta-Fibrinogen Intron 7 sequences. For species delimitation, we employed the General Mixed Yule Coalescent (GMYC) and Bayesian Poisson tree processes (bPTP) analyses, and contrasted 4 genealogical and phylogenetic models: Pure-birth (Yule), Constant Population Size Coalescent, Multiple Species Coalescent, and a mixed Yule-Coalescent model. GMYC analyses of trees from different genealogical models resulted in similar species delimitation and phylogenetic relationships, with incongruence restricted to areas of poor nodal support. bPTP results, however, significantly differed from GMYC for 5 taxa. Oligoryzomys early diversification was estimated to have occurred in the Early Pleistocene, between 0.7 and 2.6 MYA. The mixed Yule-Coalescent model, however, recovered younger dating estimates for Oligoryzomys diversification, and for the threshold for the speciation-coalescent horizon in GMYC. Eight of the 20 included Oligoryzomys species were identified as having two or more independent evolutionary units, indicating that current taxonomy of Oligoryzomys is still unsettled.