Incipient ring speciation revealed by a migratory divide

Ever since Ernst Mayr (1942) called ring species the ‘perfect demonstration of speciation’, they have attracted much interest from researchers examining how two species evolve from one. In a ring species, two sympatric and reproductively isolated forms are connected by a long chain of intermediate populations that encircle a geographic barrier. Ring species have the potential to demonstrate that speciation can occur without complete geographic isolation, in contrast to the classic model of allopatric speciation. They also allow researchers to examine the causes of reproductive isolation in the contact zone and to use spatial variation to infer the steps by which speciation occurs. According to the classical definition, a ring species must have (i) gradual variation through a chain of populations connecting two divergent and sympatric forms, and (ii) complete or nearly complete reproductive isolation between the terminal forms. But evolutionary biologists now recognize that the process of speciation might often occur with some periods of geographic contact and hybridization between diverging forms; during these phases, even partial reproductive isolation can limit gene flow and permit further divergence to occur. In this issue Bensch et al. (2009) make an exciting and important contribution by extending the ring species concept to a case in which the divergence is much younger and not yet advanced to full reproductive isolation. Their study of geographic variation in willow warblers (Phylloscopus trochilus; Fig. 1 ) provides a beautiful example of gradual variation through a ring of populations connecting two forms that are partially reproductively isolated where they meet, possibly due to divergent migratory behaviours of the terminal forms.

Figure 1 Open in figure viewer PowerPoint A male willow warbler resembling the southeastern‐migrating form (Phylloscopus trochilus acredula), on its breeding territory in central Sweden. (Photo: Anders Hedenström).     

A tree of tree frogs around the Black Sea

Speciation, the process by which one species evolves into two or more, is a major focus of ongoing debate, particularly regarding the geographic context in which it occurs. Geographic models of speciation tend to fall into discrete categories, typically referred to as allopatric, parapatric and sympatric speciation, according to whether two groups evolve reproductive isolation while geographically isolated, differentiated but connected by gene flow, or completely co‐occurring. Yet molecular studies indicate that full development of reproductive isolation can take very long compared with the timescale at which climatic oscillations occur, such that the geographic context of differentiating forms might change often during the long process to full species. Studies of genetic relationships across the ranges of organisms with low‐dispersal distances have the potential to reveal these complex histories. In a particularly elegant example in this issue, Dufresnes et al. ( 2016 ) use genetic variation and ecological niche modelling to show that a ring of populations of the eastern tree frog (Hyla orientalis) surrounding the Black Sea had a complex history of geographic differentiation. Alternating phases of geographic fragmentation and phases of gene flow between neighbouring populations have produced a pattern of gradual genetic change connecting the western, southern and eastern sides of the ring, with the northwestern and northeastern forms being most differentiated. In the north, a population in Crimea appears to have been produced through mixture of the two extreme forms. The overall genetic relationships are reminiscent of those found in ring species, which have been used as prime demonstrations of the process of speciation. The difference, however, is that the terminal forms appear to have mixed rather than be reproductively isolated, although more research is needed to infer whether there might be some reproductive isolation on the northern side of the ring.     

Repeated and predictable patterns of ecotypic differentiation during a biological invasion: lake–stream divergence in parapatric Swiss stickleback

The relative importance of ecological selection and geographical isolation in promoting and constraining genetic and phenotypic differentiation among populations is not always obvious. Interacting with divergent selection, restricted opportunity for gene flow may in some cases be as much a cause as a consequence of adaptation, with the latter being a hallmark of ecological speciation. Ecological speciation is well studied in parts of the native range of the three‐spined stickleback. Here, we study this process in a recently invaded part of its range. Switzerland was colonized within the past 140 years from at least three different colonization events involving different stickleback lineages. They now occupy diverse habitats, ranging from small streams to the pelagic zone of large lakes. We use replicated systems of parapatric lake and stream populations, some of which trace their origins to different invasive lineages, to ask (i) whether phenotypic divergence occurred among populations inhabiting distinct habitats, (ii) whether trajectories of phenotypic divergence follow predictable parallel patterns and (iii) whether gene flow constrains divergent adaptation or vice versa. We find consistent phenotypic divergence between populations occupying distinct habitats. This involves parallel evolution in several traits with known ecological relevance in independent evolutionary lineages. Adaptive divergence supersedes homogenizing gene flow even at a small spatial scale. We find evidence that adaptive phenotypic divergence places constraints on gene flow over and above that imposed by geographical distance, signalling the early onset of ecological speciation.     

Genomics of Rapid Incipient Speciation in Sympatric Threespine Stickleback

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this.     

Song variation in an avian ring species

Abstract.— Divergence of mating signals can occur rapidly and be of prime importance in causing reproductive isolation and speciation. A ring species, in which two reproductively isolated taxa are connected by a chain of intergrading populations, provides a rare opportunity to use spatial variation to reconstruct the history of divergence. I use geographic variation in the song of a likely ring species, the greenish warbler ( Phylloscopus trochiloides ) to reconstruct the microevolutionary steps that occurred during divergence of a trait that is often important in speciation in birds. Populations of a western Siberian ( P. t. viridanus ) and an eastern Siberian ( P. t. plumbeitarsus ) form of the greenish warbler meet, but do not interbreed in central Siberia; these forms are connected by a chain of interbreeding populations extending in a ring to the south around the treeless Tibetan Plateau. I show that: (1) song structure differs greatly between the two Siberian forms, which share the same habitat; (2) song structure changes gradually around the ring; (3) singing behavior is relatively simple in the Himalayas, but becomes increasingly complex to the north, both to the west and east of the Tibetan Plateau; and (4) song varies along independent axes of complexity in the western and eastern south-north clines. By comparing geographic variation in singing behavior and ecological variables, I distinguish among possible causes of song divergence, including selection based on the acoustic environment, stochastic effects of sexual selection, and selection for species recognition. I suggest that parallel south-to-north ecological gradients have caused a greater intensity of sexual selection on song in northern populations and that the stochastic effects of sexual selection have led to divergence in song structure.     

Speciation along a depth gradient in a marine adaptive radiation

Oceans are home to much of the world''s biodiversity, but we know little about the processes driving speciation in marine ecosystems with few geographical barriers to gene flow. Ecological speciation resulting from divergent natural selection between ecological niches can occur in the face of gene flow. Sister species in the young and ecologically diverse rockfish genus Sebastes coexist in the northeast Pacific, implying that speciation may not require geographical isolation. Here, I use a novel phylogenetic comparative analysis to show that rockfish speciation is instead associated with divergence in habitat depth and depth-associated morphology, consistent with models of parapatric speciation. Using the same analysis, I find no support for alternative hypotheses that speciation involves divergence in diet or life history, or that speciation involves geographic isolation by latitude. These findings support the hypothesis that rockfishes undergo ecological speciation on an environmental gradient.     

Population differentiation and restricted gene flow in Spanish crossbills: not isolation-by-distance but isolation-by-ecology

Divergent selection stemming from environmental variation may induce local adaptation and ecological speciation whereas gene flow might have a homogenizing effect. Gene flow among populations using different environments can be reduced by geographical distance (isolation-by-distance) or by divergent selection stemming from resource use (isolation-by-ecology). We tested for and encountered phenotypic and genetic divergence among Spanish crossbills utilizing different species of co-occurring pine trees as their food resource. Morphological, vocal and mtDNA divergence were not correlated with geographical distance, but they were correlated with differences in resource use. Resource diversity has now been found to repeatedly predict crossbill diversity. However, when resource use is not 100% differentiated, additional characters (morphological, vocal, genetic) must be used to uncover and validate hidden population structure. In general, this confirms that ecology drives adaptive divergence and limits neutral gene flow as the first steps towards ecological speciation, unprevented by a high potential for gene flow.     

Evaluation of genetic isolation within an island flora reveals unusually widespread local adaptation and supports sympatric speciation

It is now recognized that speciation can proceed even when divergent natural selection is opposed by gene flow. Understanding the extent to which environmental gradients and geographical distance can limit gene flow within species can shed light on the relative roles of selection and dispersal limitation during the early stages of population divergence and speciation. On the remote Lord Howe Island (Australia), ecological speciation with gene flow is thought to have taken place in several plant genera. The aim of this study was to establish the contributions of isolation by environment (IBE) and isolation by community (IBC) to the genetic structure of 19 plant species, from a number of distantly related families, which have been subjected to similar environmental pressures over comparable time scales. We applied an individual-based, multivariate, model averaging approach to quantify IBE and IBC, while controlling for isolation by distance (IBD). Our analyses demonstrated that all species experienced some degree of ecologically driven isolation, whereas only 12 of 19 species were subjected to IBD. The prevalence of IBE within these plant species indicates that divergent selection in plants frequently produces local adaptation and supports hypotheses that ecological divergence can drive speciation in sympatry.     

Ernst Mayr and the integration of geographic and ecological factors in speciation

Mayr's best recognized scientific contributions include the biological species concept and the theory of geographic speciation. In the latter, reproductive isolation evolves as an incidental by‐product of genetic divergence between allopatric populations. Mayr noted that divergent natural selection could accelerate speciation, but also argued that gene flow so strongly retards divergence that, even with selection, non‐allopatric speciation is unlikely. However, current theory and data demonstrate that substantial divergence, and even speciation, in the face of gene flow is possible. Here, I attempt to connect some opposing views about speciation by integrating Mayr's ideas about the roles of ecology and geography in speciation with current data and theory. My central premise is that the speciation process (i.e. divergence) is often continuous, and that the opposing processes of selection and gene flow interact to determine the degree of divergence (i.e. the degree of progress towards the completion of speciation). I first establish that, in the absence of gene flow, divergent selection often promotes speciation. I then discuss how population differentiation in the face of gene flow is common when divergent selection occurs. However, such population differentiation does not always lead to the evolution of discontinuities, strong reproductive isolation, and thus speciation per se. I therefore explore the genetic and ecological circumstances that facilitate speciation in the face of gene flow. For example, particular genetic architectures or ecological niches may tip the balance between selection and gene flow strongly in favour of selection. The circumstances allowing selection to overcome gene flow to the extent that a discontinuity develops, and how often these circumstances occur, are major remaining questions in speciation research. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 26–46.     

Model‐based demographic inference of introgression history in European whitefish species pairs'

Parallel phenotypic differentiation is generally attributed to parallel adaptive divergence as an evolutionary response to similar environmental contrasts. Such parallelism may actually originate from several evolutionary scenarios ranging from repeated parallel divergence caused by divergent selection to a unique divergence event followed by gene flow. Reconstructing the evolutionary history underlying parallel phenotypic differentiation is thus fundamental to understand the relative contribution of demography and selection on genomic divergence during speciation. In this study, we investigate the divergence history of replicate European whitefish (Coregonus lavaretus), limnetic and benthic species pairs from two lakes in Norway and two lakes in Switzerland. Demographic models accounting for semi‐permeability and linked selection were fitted to the unfolded joint allele frequency spectrum built from genome‐wide SNPs and compared to each other in each species pair. We found strong support for a model of asymmetrical post‐glacial secondary contact between glacial lineages in all four lakes. Moreover, our results suggest that heterogeneous genomic differentiation has been shaped by the joint action of linked selection accelerating lineage sorting during allopatry, and heterogeneous migration eroding divergence at different rates along the genome following secondary contact. Our analyses reveal how the interplay between demography, selection and historical contingency has influenced the levels of diversity observed in previous whitefish phylogeographic studies. This study thus provides new insights into the historical demographic and selective processes that shaped the divergence associated with ecological speciation in European whitefish.     

Asymmetric introgression between the M and S forms of the malaria vector, Anopheles gambiae, maintains divergence despite extensive hybridization

The suggestion that genetic divergence can arise and/or be maintained in the face of gene flow has been contentious since first proposed. This controversy and a rarity of good examples have limited our understanding of this process. Partially reproductively isolated taxa have been highlighted as offering unique opportunities for identifying the mechanisms underlying divergence with gene flow. The African malaria vector, Anopheles gambiae s.s., is widely regarded as consisting of two sympatric forms, thought by many to represent incipient species, the M and S molecular forms. However, there has been much debate about the extent of reproductive isolation between M and S, with one view positing that divergence may have arisen and is being maintained in the presence of gene flow, and the other proposing a more advanced speciation process with little realized gene flow because of low hybrid fitness. These hypotheses have been difficult to address because hybrids are typically rare (<1%). Here, we assess samples from an area of high hybridization and demonstrate that hybrids are fit and responsible for extensive introgression. Nonetheless, we show that strong divergent selection at a subset of loci combined with highly asymmetric introgression has enabled M and S to remain genetically differentiated despite extensive gene flow. We propose that the extent of reproductive isolation between M and S varies across West Africa resulting in a 'geographic mosaic of reproductive isolation'; a finding which adds further complexity to our understanding of divergence in this taxon and which has considerable implications for transgenic control strategies.     

Parallel habitat specialization within the wolf spider genus Hogna from the Galápagos

Within most island archipelagos, such as the Galápagos, similar ecological gradients are found on geographically isolated islands. Species radiations in response to these ecological gradients may follow different scenarios being (i) a single habitat specialization event followed by secondary colonization of each ecotype on the different islands or (ii) repeated and parallel habitat specialization on each island separately. This latter scenario has been considered less likely as gene flow might hamper such ecotypic differentiation. At least for the Galápagos, the extent to which this process is involved in species radiations remains yet poorly understood. Within the wolf spider genus Hogna, seven species are described that can be divided into three different ecotypes based on general morphology and habitat preference i.e. species that inhabit the pampa vegetation in the highlands, species that occur in coastal dry habitats and one generalist species. Comparison of the species phylogeny based on one mitochondrial (COI) and one nuclear (28S) gene fragment convincingly demonstrates that ‘pampa’ and ‘coastal dry’ species evolved in parallel on the islands Santa Cruz and San Cristóbal. Despite the observation that allozymes analysis indicated that each species forms a distinct genetic cluster, phylogenetic divergence within these species complexes was very low and paraphyletic and most likely due to hybridization rather than incomplete lineage sorting, as demonstrated for the Santa Cruz species complex. This suggests that within‐island speciation occurred under low levels of gene flow. Species phylogeny in general did not follow the progression of island emergence as a molecular clock analysis suggested that island endemic species may have diverged after as well as before the emergence of the islands. This represents the first clear example of parallel and within‐island speciation because of habitat specialization on the Galápagos and that such divergence most likely occurred under historic gene flow.     

Radiating on oceanic islands: patterns and processes of speciation in the land snail genus Theba (Risso 1826)

Island radiations have played a major role in shaping our current understanding of allopatric, sympatric and parapatric speciation. However, the fact that species divergence correlates with island size emphasizes the importance of geographic isolation (allopatry) in speciation. Based on molecular and morphological data, we investigated the diversification of the land snail genus Theba on the two Canary Islands of Lanzarote and Fuerteventura. Due to the geological history of both islands, this study system provides ideal conditions to investigate the interplay of biogeography, dispersal ability and differentiation in generating species diversity. Our analyses demonstrated extensive cryptic diversification of Theba on these islands, probably driven mainly by non-adaptive allopatric differentiation and secondary gene flow. In a few cases, we observed a complete absence of gene flow among sympatrically distributed forms suggesting an advanced stage of speciation. On the Jandía peninsula genome scans suggested genotype-environment associations and potentially adaptive diversification of two closely related Theba species to different ecological environments. We found support for the idea that genetic differentiation was enhanced by divergent selection in different environments. The diversification of Theba on both islands is therefore best explained by a mixture of non-adaptive and adaptive speciation, promoted by ecological and geomorphological factors.     

Phenotypic divergence and reproductive isolation between sympatric forms of Japanese threespine sticklebacks

The threespine stickleback ( Gasterosteus aculeatus ) species complex is well suited for identifying the types of phenotypic divergence and isolating barriers that contribute to reproductive isolation at early stages of speciation. In the present study, we characterize the patterns of genetic and phenotypic divergence as well as the types of isolating barriers that are present between two sympatric pairs of threespine sticklebacks in Hokkaido, Japan. One sympatric pair consists of an anadromous and a resident freshwater form and shows divergence in body size between the forms, despite the lack of genetic differentiation between them. The second sympatric pair consists of two anadromous forms, which originated before the last glacial period and are currently reproductively isolated. These two anadromous forms have diverged in many morphological traits as well as in their reproductive behaviours. Both sexual isolation and hybrid male sterility contribute to reproductive isolation between the anadromous species pair. We discuss the shared and unique aspects of phenotypic divergence and reproductive isolation in the Japanese sympatric pairs compared with postglacial stickleback species pairs. Further studies of these divergent species pairs will provide a deeper understanding of the mechanisms of speciation in sticklebacks.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 671–685.     

Theoretical models of the influence of genomic architecture on the dynamics of speciation

A long‐standing problem in evolutionary biology has been determining whether and how gradual, incremental changes at the gene level can account for rapid speciation and bursts of adaptive radiation. Using genome‐scale computer simulations, we extend previous theory showing how gradual adaptive change can generate nonlinear population transitions, resulting in the rapid formation of new, reproductively isolated species. We show that these transitions occur via a mechanism rooted in a basic property of biological heredity: the organization of genes in genomes. Genomic organization of genes facilitates two processes: (i) the build‐up of statistical associations among large numbers of genes and (ii) the action of divergent selection on persistent combinations of alleles. When a population has accumulated a critical amount of standing, divergently selected variation, the combination of these two processes allows many mutations of small effect to act synergistically and precipitously split one population into two discontinuous, reproductively isolated groups. Periods of allopatry, chromosomal linkage among loci, and large‐effect alleles can facilitate this process under some conditions, but are not required for it. Our results complement and extend existing theory on alternative stable states during population divergence, distinct phases of speciation and the rapid emergence of multilocus barriers to gene flow. The results are thus a step towards aligning population genomic theory with modern empirical studies.     

Genomic regions with a history of divergent selection affect fitness of hybrids between two butterfly species

Speciation is the process by which reproductively isolated lineages arise, and is one of the fundamental means by which the diversity of life increases. Whereas numerous studies have documented an association between ecological divergence and reproductive isolation, relatively little is known about the role of natural selection in genome divergence during the process of speciation. Here, we use genome-wide DNA sequences and Bayesian models to test the hypothesis that loci under divergent selection between two butterfly species (Lycaeides idas and L. melissa) also affect fitness in an admixed population. Locus-specific measures of genetic differentiation between L. idas and L. melissa and genomic introgression in hybrids varied across the genome. The most differentiated genetic regions were characterized by elevated L. idas ancestry in the admixed population, which occurs in L. idas-like habitat, consistent with the hypothesis that local adaptation contributes to speciation. Moreover, locus-specific measures of genetic differentiation (a metric of divergent selection) were positively associated with extreme genomic introgression (a metric of hybrid fitness). Interestingly, concordance of differentiation and introgression was only partial. We discuss multiple, complementary explanations for this partial concordance.     

Divergence with gene flow and fine-scale phylogeographical structure in the wedge-billed woodcreeper, Glyphorynchus spirurus, a Neotropical rainforest bird

Determining the relative roles of vicariance and selection in restricting gene flow between populations is of central importance to the evolutionary process of population divergence and speciation. Here we use molecular and morphological data to contrast the effect of isolation (by mountains and geographical distance) with that of ecological factors (altitudinal gradients) in promoting differentiation in the wedge-billed woodcreeper, Glyphorynchus spirurus , a tropical forest bird, in Ecuador. Tarsus length and beak size increased relative to body size with altitude on both sides of the Andes, and were correlated with the amount of moss on tree trunks, suggesting the role of selection in driving adaptive divergence. In contrast, molecular data revealed a considerable degree of admixture along these altitudinal gradients, suggesting that adaptive divergence in morphological traits has occurred in the presence of gene flow. As suggested by mitochondrial DNA sequence data, the Andes act as a barrier to gene flow between ancient subspecific lineages. Genome-wide amplified fragment length polymorphism markers reflected more recent patterns of gene flow and revealed fine-scale patterns of population differentiation that were not detectable with mitochondrial DNA, including the differentiation of isolated coastal populations west of the Andes. Our results support the predominant role of geographical isolation in driving genetic differentiation in G. spirurus , yet suggest the role of selection in driving parallel morphological divergence along ecological gradients.     

Genome‐wide tests for introgression between cactophilic Drosophila implicate a role of inversions during speciation

Models of speciation‐with‐gene‐flow have shown that the reduction in recombination between alternative chromosome arrangements can facilitate the fixation of locally adaptive genes in the face of gene flow and contribute to speciation. However, it has proven frustratingly difficult to show empirically that inversions have reduced gene flow and arose during or shortly after the onset of species divergence rather than represent ancestral polymorphisms. Here, we present an analysis of whole genome data from a pair of cactophilic fruit flies, Drosophila mojavensis and D. arizonae, which are reproductively isolated in the wild and differ by several large inversions on three chromosomes. We found an increase in divergence at rearranged compared to colinear chromosomes. Using the density of divergent sites in short sequence blocks we fit a series of explicit models of species divergence in which gene flow is restricted to an initial period after divergence and may differ between colinear and rearranged parts of the genome. These analyses show that D. mojavensis and D. arizonae have experienced postdivergence gene flow that ceased around 270 KY ago and was significantly reduced in chromosomes with fixed inversions. Moreover, we show that these inversions most likely originated around the time of species divergence which is compatible with theoretical models that posit a role of inversions in speciation with gene flow.     

Body size evolution simultaneously creates and collapses species boundaries in a clade of scincid lizards

Speciation is generally viewed as an irreversible process, although habitat alterations can erase reproductive barriers if divergence between ecologically differentiated species is recent. Reversed speciation might also occur if geographical contact is established between species that have evolved the same reproductive isolating barrier in parallel. Here, we demonstrate a loss of intrinsic reproductive isolation in a clade of scincid lizards as a result of parallel body size evolution, which has allowed for gene flow where large-bodied lineages are in secondary contact. An mtDNA phylogeny confirms the monophyly of the Plestiodon skiltonianus species complex, but rejects that of two size-differentiated ecomorphs. Mate compatibility experiments show that the high degree of body size divergence imposes a strong reproductive barrier between the two morphs; however, the strength of the barrier is greatly diminished between parallel-evolved forms. Since two large-bodied lineages are in geographical contact in the Sierra Nevada Mountains of California, we were also able to test for postzygotic isolation under natural conditions. Analyses of amplified fragment length polymorphisms show that extensive gene exchange is occurring across the contact zone, resulting in an overall pattern consistent with isolation by distance. These results provide evidence of reversed speciation between clades that diverged from a common ancestor more than 12Myr ago.     

Glacial cycles drive rapid divergence of cryptic field vole species

Understanding the factors that contribute to the generation of reproductively isolated forms is a fundamental goal of evolutionary biology. Cryptic species are an especially interesting challenge to study in this context since they lack obvious morphological differentiation that provides clues to adaptive divergence that may drive reproductive isolation. Geographical isolation in refugial areas during glacial cycling is known to be important for generating genetically divergent populations, but its role in the origination of new species is still not fully understood and likely to be situation dependent. We combine analysis of 35,434 single‐nucleotide polymorphisms (SNPs) with environmental niche modeling (ENM) to investigate genomic and ecological divergence in three cryptic species formerly classified as the field vole (Microtus agrestis). The SNPs demonstrate high genomic divergence (pairwise F_ST values of 0.45–0.72) and little evidence of gene flow among the three field vole cryptic species, and we argue that genetic drift may have been a particularly important mechanism for divergence in the group. The ENM reveals three areas as potential glacial refugia for the cryptic species and differing climatic niches, although with spatial overlap between species pairs. This evidence underscores the role that glacial cycling has in promoting genetic differentiation and reproductive isolation by subdivision into disjunct distributions at glacial maxima in areas relatively close to ice sheets. Future investigation of the intrinsic barriers to gene flow between the field vole cryptic species is required to fully assess the mechanisms that contribute to reproductive isolation. In addition, the Portuguese field vole (M. rozianus) shows a high inbreeding coefficient and a restricted climatic niche, and warrants investigation into its conservation status.