CHROMOSOMAL REARRANGEMENTS AND THE GENETICS OF REPRODUCTIVE BARRIERS IN MIMULUS (MONKEY FLOWERS)

Chromosomal rearrangements may directly cause hybrid sterility and can facilitate speciation by preserving local adaptation in the face of gene flow. We used comparative linkage mapping with shared gene‐based markers to identify potential chromosomal rearrangements between the sister monkeyflowers Mimulus lewisii and Mimulus cardinalis, which are textbook examples of ecological speciation. We then remapped quantitative trait loci (QTLs) for floral traits and flowering time (premating isolation) and hybrid sterility (postzygotic isolation). We identified three major regions of recombination suppression in the M. lewisii × M. cardinalis hybrid map compared to a relatively collinear Mimulus parishii × M. lewisii map, consistent with a reciprocal translocation and two inversions specific to M. cardinalis. These inferences were supported by targeted intraspecific mapping, which also implied a M. lewisii‐specific reciprocal translocation causing chromosomal pseudo‐linkage in both hybrid mapping populations. Floral QTLs mapped in this study, along with previously mapped adaptive QTLs, were clustered in putatively rearranged regions. All QTLs for male sterility, including two underdominant loci, mapped to regions of recombination suppression. We argue that chromosomal rearrangements may have played an important role in generating and consolidating barriers to gene flow as natural selection drove the dramatic ecological and morphological divergence of these species.     

Chromosomal rearrangements directly cause underdominant F1 pollen sterility in Mimulus lewisii–Mimulus cardinalis hybrids

Chromosomal rearrangements can contribute to the evolution of postzygotic reproductive isolation directly, by disrupting meiosis in F₁ hybrids, or indirectly, by suppressing recombination among genic incompatibilities. Because direct effects of rearrangements on fertility imply fitness costs during their spread, understanding the mechanism of F₁ hybrid sterility is integral to reconstructing the role(s) of rearrangements in speciation. In hybrids between monkeyflowers Mimulus cardinalis and Mimulus lewisii, rearrangements contain all quantitative trait loci (QTLs) for both premating barriers and pollen sterility, suggesting that they may have facilitated speciation in this model system. We used artificial chromosome doubling and comparative mapping to test whether heterozygous rearrangements directly cause underdominant male sterility in M. lewisii–M. cardinalis hybrids. Consistent with a direct chromosomal basis for hybrid sterility, synthetic tetraploid F₁s showed highly restored fertility (83.4% pollen fertility) relative to diploids F₁s (36.0%). Additional mapping with Mimulus parishii–M. cardinalis and M. parishii–M. lewisii hybrids demonstrated that underdominant male sterility is caused by one M. lewisii specific and one M. cardinalis specific reciprocal translocation, but that inversions had no direct effects on fertility. We discuss the importance of translocations as causes of reproductive isolation, and consider models for how underdominant rearrangements spread and fix despite intrinsic fitness costs.     

Growth and leaf physiology of monkeyflowers with different altitude ranges

Every species is limited both geographically and ecologically to a subset of available habitats, yet for many species the causes of distribution limits are unknown. Temperature is thought to be one of the primary determinants of species distributions along latitudinal and altitudinal gradients. This study examined leaf physiology and plant performance under contrasting temperature regimes of sister species of monkeyflower, Mimulus cardinalis and Mimulus lewisii (Phrymaceae), that differ in altitude distribution to test the hypothesis that temperature is the primary determinant of differences in fitness versus altitude. Each species attained greatest aboveground biomass, net photosynthetic rate, and effective quantum yield of photosystem II when grown under temperatures characteristic of the altitudinal range center. Although both species exhibited greater stem length, stomatal conductance, and intercellular CO₂ concentration in hot than in cold temperatures, these traits showed much greater reductions under cold temperature for M. cardinalis than for M. lewisii. Survival of M. lewisii was also sensitive to temperature, showing a striking decrease in hot temperatures. Within each temperature regime, the species native to that temperature displayed greatest growth and leaf physiological capacity. Populations from the altitude range center and range margin of each species were used to examine population differentiation, but central and marginal populations did not differ in most growth or leaf physiological responses to temperature. This study provides evidence that M. cardinalis and M. lewisii differ in survival, growth, and leaf physiology under temperature regimes characterizing their contrasting low and high altitude range centers, and suggests that the species’ altitude range limits may arise, in part, due to metabolic limitations on growth that ultimately decrease survival and limit reproduction.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Speciational history of North American Haemorhous finches (Aves: Fringillidae) inferred from multilocus data

We investigated species relationships and timing of speciation in North American Haemorhous finches by using a mitochondrial phylogeographic approach combined with a multilocus species tree reconstruction. Haemorhous purpureus and H. cassinii were strongly supported as sister taxa, and H. mexicanus was sister to H. purpureus + H. cassinii. Our divergence times indicated that diversification within Haemorhous occurred progressively from the Late Miocene into the Pleistocene. Our inferred pattern of speciation demonstrates the complexity of the origins of North American birds, and provides additional evidence that a single cause for speciation in closely related North American birds, such as Late Pleistocene glacial–interglacial cycles, is unlikely.     

Molecular phylogenetics at the population/species interface in cave spiders of the southern Appalachians (Araneae:Nesticidae:Nesticus)

This paper focuses on the relationship between population genetic structure and speciation mechanisms in a monophyletic species group of Appalachian cave spiders (Nesticus). Using mtDNA sequence data gathered from 256 individuals, I analyzed patterns of genetic variation within and between populations for three pairs of closely related sister species. Each sister-pair comparison involves taxa with differing distributional and ecological attributes; if these ecological attributes are reflected in basic demographic differences, then speciation might proceed differently across these sister taxa comparisons. Both frequency-based and gene tree analyses reveal that the genetic structure of the Nesticus species studied is characterized by similar and essentially complete population subdivision, regardless of differences in general ecology. These findings contrast with results of prior genetic studies of cave-dwelling arthropods that have typically revealed variation in population structure corresponding to differences in general ecology. Species fragmentation through both extrinsic and intrinsic evolutionary forces has resulted in discrete, perhaps independent, populations within morphologically defined species. Large sequence divergence values observed between populations suggest that this independence may extend well into the past. These patterns of mtDNA genealogical structure and divergence imply that species as morphological lineages are currently more inclusive than basal evolutionary or phylogenetic units, a suggestion that has important implications for the study of speciation mechanisms.      

Speciation slowing down in widespread and long-living tree taxa: insights from the tropical timber tree genus Milicia (Moraceae)

The long generation time and large effective size of widespread forest tree species can result in slow evolutionary rate and incomplete lineage sorting, complicating species delimitation. We addressed this issue with the African timber tree genus Milicia that comprises two morphologically similar and often confounded species: M. excelsa, widespread from West to East Africa, and M. regia, endemic to West Africa. We combined information from nuclear microsatellites (nSSRs), nuclear and plastid DNA sequences, and morphological systematics to identify significant evolutionary units and infer their evolutionary and biogeographical history. We detected five geographically coherent genetic clusters using nSSRs and three levels of genetic differentiation. First, one West African cluster matched perfectly with the morphospecies M. regia that formed a monophyletic clade at both DNA sequences. Second, a West African M. excelsa cluster formed a monophyletic group at plastid DNA and was more related to M. regia than to Central African M. excelsa, but shared many haplotypes with the latter at nuclear DNA. Third, three Central African clusters appeared little differentiated and shared most of their haplotypes. Although gene tree paraphyly could suggest a single species in Milicia following the phylogenetic species concept, the existence of mutual haplotypic exclusivity and nonadmixed genetic clusters in the contact area of the two taxa indicate strong reproductive isolation and, thus, two species following the biological species concept. Molecular dating of the first divergence events showed that speciation in Milicia is ancient (Tertiary), indicating that long-living tree taxa exhibiting genetic speciation may remain similar morphologically.     

Comparative linkage maps suggest that fission,not polyploidy,underlies near-doubling of chromosome number within monkeyflowers (Mimulus; Phrymaceae)

Changes in chromosome number and structure are important contributors to adaptation, speciation and macroevolution. In flowering plants, polyploidy and subsequent reductions in chromosome number by fusion are major sources of chromosomal evolution, but chromosome number increase by fission has been relatively unexplored. Here, we use comparative linkage mapping with gene-based markers to reconstruct chromosomal synteny within the model flowering plant genus Mimulus (monkeyflowers). Two sections of the genus with haploid numbers ⩾14 have been inferred to be relatively recent polyploids because they are phylogenetically nested within numerous taxa with low base numbers (n=8–10). We combined multiple data sets to build integrated genetic maps of the M. guttatus species complex (section Simiolus, n=14) and the M. lewisii group (section Erythranthe; n=8), and then aligned the two integrated maps using >100 shared markers. We observed strong segmental synteny between M. lewisii and M. guttatus maps, with essentially 1-to-1 correspondence across each of 16 chromosomal blocks. Assuming that the M. lewisii (and widespread) base number of 8 is ancestral, reconstruction of 14 M. guttatus chromosomes requires at least eight fission events (likely shared by Simiolus and sister section Paradanthus (n=16)), plus two fusion events. This apparent burst of fission in the yellow monkeyflower lineages raises new questions about mechanisms and consequences of chromosomal fission in plants. Our comparative maps also provide insight into the origins of a chromosome exhibiting centromere-associated female meiotic drive and create a framework for transferring M. guttatus genome resources across the entire genus.     

The Genetic Basis of a Rare Flower Color Polymorphism in Mimulus lewisii Provides Insight into the Repeatability of Evolution

A long-standing question in evolutionary biology asks whether the genetic changes contributing to phenotypic evolution are predictable. Here, we identify a genetic change associated with segregating variation in flower color within a population of Mimulus lewisii. To determine whether these types of changes are predictable, we combined this information with data from other species to investigate whether the spectrum of mutations affecting flower color transitions differs based on the evolutionary time-scale since divergence. We used classic genetic techniques, along with gene expression and population genetic approaches, to identify the putative, loss-of-function mutation that generates rare, white flowers instead of the common, pink color in M. lewisii. We found that a frameshift mutation in an anthocyanin pathway gene is responsible for the white-flowered polymorphism found in this population of M. lewisii. Comparison of our results with data from other species reveals a broader spectrum of flower color mutations segregating within populations relative to those that fix between populations. These results suggest that the genetic basis of fixed differences in flower color may be predictable, but that for segregating variation is not.     

Direct Gamete Sequencing Reveals No Evidence for Segregation Distortion in House Mouse Hybrids

Understanding the molecular basis of species formation is an important goal in evolutionary genetics, and Dobzhansky-Muller incompatibilities are thought to be a common source of postzygotic reproductive isolation between closely related lineages. However, the evolutionary forces that lead to the accumulation of such incompatibilities between diverging taxa are poorly understood. Segregation distorters are believed to be an important source of Dobzhansky-Muller incompatibilities between hybridizing species of Drosophila as well as hybridizing crop plants, but it remains unclear if these selfish genetic elements contribute to reproductive isolation in other taxa. Here, we collected viable sperm from first-generation hybrid male progeny of Mus musculus castaneus and M. m. domesticus, two subspecies of rodent in the earliest stages of speciation. We then genotyped millions of single nucleotide polymorphisms in these gamete pools and tested for a skew in the frequency of parental alleles across the genome. We show that segregation distorters are not measurable contributors to observed infertility in these hybrid males, despite sufficient statistical power to detect even weak segregation distortion with our novel method. Thus, reduced hybrid male fertility in crosses between these nascent species is attributable to other evolutionary forces.     

Genomics of speciation and introgression in Princess cichlid fishes from Lake Tanganyika

How variation in the genome translates into biological diversity and new species originate has endured as the mystery of mysteries in evolutionary biology. African cichlid fishes are prime model systems to address speciation‐related questions for their remarkable taxonomic and phenotypic diversity, and the possible role of gene flow in this process. Here, we capitalize on genome sequencing and phylogenomic analyses to address the relative impacts of incomplete lineage sorting, introgression and hybrid speciation in the Neolamprologus savoryi‐complex (the ‘Princess cichlids’) from Lake Tanganyika. We present a time‐calibrated species tree based on whole‐genome sequences and provide strong evidence for incomplete lineage sorting in the early phases of diversification and multiple introgression events affecting different stages. Importantly, we find that the Neolamprologus chromosomes show centre‐to‐periphery biases in nucleotide diversity, sequence divergence, GC content, incomplete lineage sorting and rates of introgression, which are likely modulated by recombination density and linked selection. The detection of heterogeneous genomic landscapes has strong implications on the genomic mechanisms involved in speciation. Collinear chromosomal regions can be protected from gene flow and harbour incompatibility genes if they reside in lowly recombining regions, and coupling can evolve between nonphysically linked genomic regions (chromosome centres in particular). Simultaneously, higher recombination towards chromosome peripheries makes these more dynamic, evolvable regions where adaptation polymorphisms have a fertile ground. Hence, differences in genome architecture could explain the levels of taxonomic and phenotypic diversity seen in taxa with collinear genomes and might have contributed to the spectacular cichlid diversity observed today.     

‘Lassoing’ a phylogenetic tree I: basic properties, shellings, and covers

A classical result, fundamental to evolutionary biology, states that an edge-weighted tree T with leaf set X, positive edge weights, and no vertices of degree 2 can be uniquely reconstructed from the leaf-to-leaf distances between any two elements of X. In biology, X corresponds to a set of taxa (e.g. extant species), the tree T describes their phylogenetic relationships, the edges correspond to earlier species evolving for a time until splitting in two or more species by some speciation/bifurcation event, and their length corresponds to the genetic change accumulating over that time in such a species. In this paper, we investigate which subsets of ${\binom{X}{2}}$ suffice to determine (??lasso??) the tree T from the leaf-to-leaf distances induced by that tree. The question is particularly topical since reliable estimates of genetic distance??even (if not in particular) by modern mass-sequencing methods??are, in general, available only for certain combinations of taxa.     

Temperature and photoperiodic control of developmental responses in climatic races of Mimulus

Two strongly differentiated climatic races of the Mimulus cardinalis-lewisiicomplex were grown at a variety of temperatures (3–27°C)and photoperiods (8 and 16 hr) under controlled environmentalconditions. M. cardinalis (the lowland race, 400 m) and M. lewisii(the sub-alpine race, 3200 m) were found to differ in theirphysiological responses to the varied environments in severalsignificant ways: 1) At 27°C (16 hrphotoperiod), M. lewisiisustained 100% mortality in contrast to the substantial growthand flowering of M. cardinalis under these conditions; 2) In8 hr photoperiods at all temperatures, there was little growthand no flowering in M. lewisii whereas there was considerablegrowth at all temperatures, and flowering at 23 and 27°Cin M. cardinalis; 3) At low temperatures (7–15°C),16 hr photoperiods, flowering occurred a week or two earlierin M. lewisii than in M. cardinalis. The lowland race has asignificantly wider temperature and photoperiodic tolerancethan has the sub-alpine race. Applications of gibberellic acidto rosette Mimulus plants under non-inductive conditions (15°C,8 hr photoperiod) promoted vigorous stem elongation withoutflowering. The application of steroids, other hormones and metaboliteshad no observable effects. ¹Present address: Faculty of Botany, Ohio State University,Columbus, Ohio, U.S.A. (Received March 16, 1970; )     

The geographical pattern of speciation and floral diversification in the neotropics: the tribe sinningieae (gesneriaceae) as a case study

The geographical pattern of speciation and the relationship between floral variation and species ranges were investigated in the tribe Sinningieae (Gesneriaceae), which is found mainly in the Atlantic forests of Brazil. Geographical distribution data recorded on a grid system of 0.5 x 0.5 degree intervals and a near-complete species-level phylogenetic tree of Sinningieae inferred from a simultaneous analysis of seven DNA regions were used to address the role of geographical isolation in speciation. Geographical range overlaps between sister lineages were measured across all nodes in the phylogenetic tree and analyzed in relation to relative ages estimated from branch lengths. Although there are several cases of species sympatry in Sinningieae, patterns of sympatry between sister taxa support the predominance of allopatric speciation. The pattern of sympatry between sister taxa is consistent with range shifts following allopatric speciation, except in one clade, in which the overlapping distribution of recent sister species indicates speciation within a restricted geographical area and involving changes in pollinators and habitats. The relationship between floral divergence and regional sympatry was also examined by analyzing floral contrasts, phenological overlap, and the degree of sympatry between sister clades. Morphological contrast between flowers is not increased in sympatry and phenological divergence is more apparent between allopatric clades than between sympatric clades. Therefore, our results failed to indicate a tendency for sympatric taxa to minimize morphological and phenological overlap (geographic exclusion and/or character displacement hypotheses). Instead, they point toward adaptation in phenology to local conditions and buildup of sympatries at random with respect to flower morphology. Additional studies at a lower geographical scale are needed to identify truely coexisting species and the components of their reproductive isolation.     

Nuclear Species-Diagnostic SNP Markers Mined from 454 Amplicon Sequencing Reveal Admixture Genomic Structure of Modern Citrus Varieties

Most cultivated Citrus species originated from interspecific hybridisation between four ancestral taxa (C. reticulata, C. maxima, C. medica, and C. micrantha) with limited further interspecific recombination due to vegetative propagation. This evolution resulted in admixture genomes with frequent interspecific heterozygosity. Moreover, a major part of the phenotypic diversity of edible citrus results from the initial differentiation between these taxa. Deciphering the phylogenomic structure of citrus germplasm is therefore essential for an efficient utilization of citrus biodiversity in breeding schemes. The objective of this work was to develop a set of species-diagnostic single nucleotide polymorphism (SNP) markers for the four Citrus ancestral taxa covering the nine chromosomes, and to use these markers to infer the phylogenomic structure of secondary species and modern cultivars. Species-diagnostic SNPs were mined from 454 amplicon sequencing of 57 gene fragments from 26 genotypes of the four basic taxa. Of the 1,053 SNPs mined from 28,507 kb sequence, 273 were found to be highly diagnostic for a single basic taxon. Species-diagnostic SNP markers (105) were used to analyse the admixture structure of varieties and rootstocks. This revealed C. maxima introgressions in most of the old and in all recent selections of mandarins, and suggested that C. reticulata × C. maxima reticulation and introgression processes were important in edible mandarin domestication. The large range of phylogenomic constitutions between C. reticulata and C. maxima revealed in mandarins, tangelos, tangors, sweet oranges, sour oranges, grapefruits, and orangelos is favourable for genetic association studies based on phylogenomic structures of the germplasm. Inferred admixture structures were in agreement with previous hypotheses regarding the origin of several secondary species and also revealed the probable origin of several acid citrus varieties. The developed species-diagnostic SNP marker set will be useful for systematic estimation of admixture structure of citrus germplasm and for diverse genetic studies.     

Whole genome sequencing reveals the genomic diversity,taxonomic classification,and evolutionary relationships of the genus Nocardia

Nocardia is a complex and diverse genus of aerobic actinomycetes that cause complex clinical presentations, which are difficult to diagnose due to being misunderstood. To date, the genetic diversity, evolution, and taxonomic structure of the genus Nocardia are still unclear. In this study, we investigated the pan-genome of 86 Nocardia type strains to clarify their genetic diversity. Our study revealed an open pan-genome for Nocardia containing 265,836 gene families, with about 99.7% of the pan-genome being variable. Horizontal gene transfer appears to have been an important evolutionary driver of genetic diversity shaping the Nocardia genome and may have caused historical taxonomic confusion from other taxa (primarily Rhodococcus, Skermania, Aldersonia, and Mycobacterium). Based on single-copy gene families, we established a high-accuracy phylogenomic approach for Nocardia using 229 genome sequences. Furthermore, we found 28 potentially new species and reclassified 16 strains. Finally, by comparing the topology between a phylogenomic tree and 384 phylogenetic trees (from 384 single-copy genes from the core genome), we identified a novel locus for inferring the phylogeny of this genus. The dapb1 gene, which encodes dipeptidyl aminopeptidase BI, was far superior to commonly used markers for Nocardia and yielded a topology almost identical to that of genome-based phylogeny. In conclusion, the present study provides insights into the genetic diversity, contributes a robust framework for the taxonomic classification, and elucidates the evolutionary relationships of Nocardia. This framework should facilitate the development of rapid tests for the species identification of highly variable species and has given new insight into the behavior of this genus.     

Phylogenomics and multigene phylogenies decipher two new cryptic marine algae from California,Gelidium gabrielsonii and G. kathyanniae (Gelidiales,Rhodophyta)

Molecular surveys are leading to the discovery of many new cryptic species of marine algae. This is particularly true for red algal intertidal species, which exhibit a high degree of morphological convergence. DNA sequencing of recent collections of Gelidium along the coast of California, USA, identified two morphologically similar entities that differed in DNA sequence from existing species. To characterize the two new species of Gelidium and to determine their evolutionary relationships to other known taxa, phylogenomic, multigene analyses, and morphological observations were performed. Three complete mitogenomes and five plastid genomes were deciphered, including those from the new species candidates and the type materials of two closely related congeners. The mitogenomes contained 45 genes and had similar lengths (24,963–24,964 bp). The plastid genomes contained 232 genes and were roughly similar in size (175,499–177,099 bp). The organellar genomes showed a high level of gene synteny. The two Gelidium species are diminutive, turf‐forming, and superficially resemble several long established species from the Pacific Ocean. The phylogenomic analysis, multigene phylogeny, and morphological evidence confirms the recognition and naming of two new species, describe herein as G. gabrielsonii and G. kathyanniae. On the basis of the monophyly of G. coulteri, G. gabrielsonii, G. galapagense, and G. kathyanniae, we suggest that this lineage likely evolved in California. Organellar genomes provide a powerful tool for discovering cryptic intertidal species and they continue to improve our understanding of the evolutionary biology of red algae and the systematics of the Gelidiales.     

Hybridization masks speciation in the evolutionary history of the Galápagos marine iguana

The effects of the direct interaction between hybridization and speciation—two major contrasting evolutionary processes—are poorly understood. We present here the evolutionary history of the Galápagos marine iguana (Amblyrhynchus cristatus) and reveal a case of incipient within-island speciation, which is paralleled by between-island hybridization. In-depth genome-wide analyses suggest that Amblyrhynchus diverged from its sister group, the Galápagos land iguanas, around 4.5 million years ago (Ma), but divergence among extant populations is exceedingly young (less than 50 000 years). Despite Amblyrhynchus appearing as a single long-branch species phylogenetically, we find strong population structure between islands, and one case of incipient speciation of sister lineages within the same island—ostensibly initiated by volcanic events. Hybridization between both lineages is exceedingly rare, yet frequent hybridization with migrants from nearby islands is evident. The contemporary snapshot provided by highly variable markers indicates that speciation events may have occurred throughout the evolutionary history of marine iguanas, though these events are not visible in the deeper phylogenetic trees. We hypothesize that the observed interplay of speciation and hybridization might be a mechanism by which local adaptations, generated by incipient speciation, can be absorbed into a common gene pool, thereby enhancing the evolutionary potential of the species as a whole.     

Reconciling molecules and morphology: molecular systematics and biogeography of Neotropical blennies (Acanthemblemaria)

Neotropical reef fish communities are species-poor compared to those of the Indo-West Pacific. An exception to that pattern is the blenny clade Chaenopsidae, one of only three rocky and coral reef fish families largely endemic to the Neotropics. Within the chaenopsids, the genus Acanthemblemaria is the most species-rich and is characterized by elaborate spinous processes on the skull. Here we construct a species tree using five nuclear markers and compare the results to those from Bayesian and parsimony phylogenetic analyses of 60 morphological characters. The sequence-based species tree conflicted with the morphological phylogenies for Acanthemblemaria, primarily due to the convergence of a suite of characters describing the distribution of spines on the head. However, we were able to resolve some of these conflicts by performing phylogenetic analyses on suites of characters not associated with head spines. By using the species tree as a guide, we used a quantitative method to identify suites of correlated morphological characters that, together, produce the distinctive skull phenotypes found in these fishes. A time calibrated phylogeny with nearly complete taxon sampling provided divergence time estimates that recovered a mid-Miocene origin for the genus, with a temporally and geographically complex pattern of speciation both before and after the closure of the Isthmus of Panama. Some sister taxa are broadly sympatric, but many occur in allopatry. The ability to infer the geography of speciation in Acanthemblemaria is complicated by extinctions, incomplete knowledge of their present geographic ranges and by wide-spread taxa that likely represent cryptic species complexes.     

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.