Increased genetic divergence between two closely related fir species in areas of range overlap

Because of introgressive hybridization, closely related species can be more similar to each other in areas of range overlap (parapatry or sympatry) than in areas where they are geographically isolated from each other (allopatry). Here, we report the reverse situation based on nuclear genetic divergence between two fir species, Abies chensiensis and Abies fargesii, in China, at sites where they are parapatric relative to where they are allopatric. We examined genetic divergence across 126 amplified fragment length polymorphism (AFLP) markers in a set of 172 individuals sampled from both allopatric and parapatric populations of the two species. Our analyses demonstrated that AFLP divergence was much greater between the species when comparisons were made between parapatric populations than between allopatric populations. We suggest that selection in parapatry may have largely contributed to this increased divergence.     

The role of climate for the range limits of parapatric European land salamanders

Abrupt range limits of parapatric species may serve as a model system to understand the factors that determine species’ range borders. Theory suggests that parapatric range limits can be caused by abiotic conditions along environmental gradients, biotic interactions or a combination of both. Geographic ranges of the parapatric salamanders, Salamandra salamandra and S. atra, meet in small contact zones in the European Alps and to date, the cause of parapatry and the restricted range of S. atra remain elusive. We combine multivariate approaches and climatic data analysis to explore niche differentiation among the two salamanders with respect to the available climatic environment at their contact zones. Our purpose is to evaluate whether climatic conditions explain the species’ sharp range limits or if biotic interactions may play a role for range delimitation. Analyses were carried out in three contact zones in Switzerland to assess possible geographic variation. Our results indicate that both species occur at localities with different climatic conditions as well as the presence of a strong climatic gradient across the species’ range limits. Although the species’ climatic niches differ moderately (with a wider niche breadth for S. atra), interspecific niche overlap is found. Comparisons among the contact zones confirm geographic variation in the species’ climatic niches as well as in the conditions within the geographically available space. Our results suggest that the change in climatic conditions along the recognized gradient represents a determining factor for species’ range limits within contact zones. However, our analyses of geographic variation in climatic conditions reveal that both salamander species can occur in a much wider range of conditions than observed within contact zones. This finding and the interspecific climatic niche overlap within each contact zone provides indirect evidence that biotic interactions (likely competition) between the two species may also determine their range limits.     

Applications of ecological niche modeling for species delimitation: a review and empirical evaluation using day geckos (Phelsuma) from Madagascar

Although the systematic utility of ecological niche modeling is generally well known (e.g., concerning the recognition and discovery of areas of endemism for biogeographic analyses), there has been little discussion of applications concerning species delimitation, and to date, no empirical evaluation has been conducted. However, ecological niche modeling can provide compelling evidence for allopatry between populations, and can also detect divergent ecological niches between candidate species. Here we present results for two taxonomically problematic groups of Phelsuma day geckos from Madagascar, where we integrate ecological niche modeling with mitochondrial DNA and morphological data to evaluate species limits. Despite relatively modest levels of genetic and morphological divergence, for both species groups we find divergent ecological niches between closely related species and parapatric ecological niche models. Niche models based on the new species limits provide a better fit to the known distribution than models based upon the combined (lumped) species limits. Based on these results, we elevate three subspecies of Phelsuma madagascariensis to species rank and describe a new species of Phelsuma from the P. dubia species group. Our phylogeny continues to support a major endemic radiation of Phelsuma in Madagascar, with dispersals to Pemba Island and the Mascarene Islands. We conclude that ecological niche modeling offers great potential for species delimitation, especially for taxonomic groups exhibiting low vagility and localized endemism and for groups with more poorly known distributions. In particular, niche modeling should be especially sensitive for detecting recent parapatric speciation driven by ecological divergence, when the environmental gradients driving speciation are represented within the ecological niche models.     

Along the speciation continuum in sticklebacks

Speciation can be viewed as a continuum, potentially divisible into several states: (1) continuous variation within panmictic populations, (2) partially discontinuous variation with minor reproductive isolation, (3) strongly discontinuous variation with strong but reversible reproductive isolation and (4) complete and irreversible reproductive isolation. Research on sticklebacks (Gasterosteidae) reveals factors that influence progress back and forth along this continuum, as well as transitions between the states. Most populations exist in state 1, even though some of these show evidence of disruptive selection and positive assortative mating. Transitions to state 2 seem to usually involve strong divergent selection coupled with at least a bit of geographic separation, such as parapatry (e.g. lake and stream pairs and mud and lava pairs) or allopatry (e.g. different lakes). Transitions to state 3 can occur when allopatric or parapatric populations that evolved under strong divergent selection come into secondary contact (most obviously the sympatric benthic and limnetic pairs), but might also occur between populations that remained in parapatry or allopatry. Transitions to state 4 might be decoupled from these selective processes, because the known situations of complete, or nearly complete, reproductive isolation (Japan Sea and Pacific Ocean pair and the recognized gasterosteid species) are always associated with chromosomal rearrangements and environment‐independent genetic incompatibilities. Research on sticklebacks has thus revealed complex and shifting interactions between selection, adaptation, mutation and geography during the course of speciation.     

The niches of nuthatches affect their lineage evolution differently across latitude

Ecological niche evolution can promote or hinder the differentiation of taxa and determine their distribution. Niche‐mediated evolution may differ among climatic regimes, and thus, species that occur across a wide latitudinal range offer a chance to test these heterogeneous evolutionary processes. In this study, we examine (a) how many lineages have evolved across the continent‐wide range of the Eurasian nuthatch (Sitta europaea), (b) whether the lineages’ niches are significantly divergent or conserved and (c) how their niche evolution explains their geographic distribution. Phylogenetic reconstruction and ecological niche models (ENMs) showed that the Eurasian nuthatch contained six parapatric lineages that diverged within 2 Myr and did not share identical climatic niches. However, the niche discrepancy between these distinct lineages was relatively conserved compared with the environmental differences between their ranges and thus was unlikely to drive lineage divergence. The ENMs of southern lineages tended to cross‐predict with their neighbouring lineages whereas those of northern lineages generally matched with their abutting ranges. The coalescence‐based analyses revealed more stable populations for the southern lineages than the northern ones during the last glaciation cycle. In contrast to the overlapping ENMs, the smaller parapatric distribution suggests that the southern lineages might have experienced competitive exclusion to prevent them from becoming sympatric. On the other hand, the northern lineages have expanded their ranges and their current abutting distribution might have resulted from lineages adapting to different climatic conditions in allopatry. This study suggests that niche evolution may affect lineage distribution in different ways across latitude.     

NATURAL HYBRIDIZATION BETWEEN IRIS FULVA AND IRIS HEXAGONA: PATTERN OF RIBOSOMAL DNA VARIATION

Genetic variation associated with the natural hybridization of Iris fulva and I. hexagona was investigated to test for the occurrence of introgression. These species have been viewed as a classic example of the process of introgressive hybridization (Anderson, 1949). However, more recent studies have concluded that there has not been an exchange of genetic material between I. fulva and I. hexagona (Randolph et al., 1967). Our analysis has involved the examination of both allopatric and parapatric populations of I. fulva and I. hexagona with reference to diagnostic ribosomal DNA markers. The pattern of variation in the parapatric population indicates the presence of the repeated backcrossing necessary to the process of introgressive hybridization. Indeed, in the region of parapatry, we suggest that localized introgression of ribosomal sequences has occurred into both I. fulva and I. hexagona. Significantly, we have also detected the presence of the diagnostic ribosomal markers from each species in allopatric populations of the alternate species. Our findings suggest that not only is introgressive hybridization presently occurring in parapatry between I. fulva and I. hexagona, but that past hybridization between these species has resulted in introgression into areas of allopatry.     

Competition drives trait evolution and character displacement between Mimulus species along an environmental gradient

Closely related species may evolve to coexist stably in sympatry through niche differentiation driven by in situ competition, a process termed character displacement. Alternatively, past evolution in allopatry may have already sufficiently reduced niche overlap to permit establishment in sympatry, a process called ecological sorting. The relative importance of each process to niche differentiation is contentious even though they are not mutually exclusive and are both mediated via multivariate trait evolution. We explore how competition has impacted niche differentiation in two monkeyflowers, Mimulus alsinoides and M. guttatus, which often co‐occur. Through field observations, common gardens, and competition experiments, we demonstrate that M. alsinoides is restricted to marginal habitats in sympatry and that the impacts of character displacement on niche differentiation are complex. Competition with M. guttatus alters selection gradients and has favored taller M. alsinoides with earlier seasonal flowering at low elevation and floral shape divergence at high elevation. However, no trait exhibits the pattern typically associated with character displacement, higher divergence between species in sympatry than allopatry. Thus, although character displacement was unlikely the process driving initial divergence along niche axes necessary for coexistence, we conclude that competition in sympatry has likely driven trait evolution along additional niche axes.     

Lack of congruence of genetic and niche divergence in Podarcis hispanicus complex

Niche divergence among closely related lineages can be informative on the ecological and evolutionary processes involved in differentiation, particularly in the case of cryptic species complexes. Here we compared phylogenetic relationships and niche similarity between pairs of lineages included in the Podarcis hispanicus complex to examine patterns of niche divergence and its role in the organization of current diversity patterns, as allopatric, parapatric, and sympatric lineages occur in the Western Mediterranean Basin. First, we used ecological niche models to characterize the realized climatic niche of each Podarcis hispanicus complex lineage based on topographic and climatic variables, to identify important variables, and to test for niche conservatism or divergence between pairs of lineages. Variables related to precipitation generally exhibited the highest contribution to niche models, highlighting the importance of rainfall levels in shaping distributions of Podarcis wall lizards. We found that most forms have significant differences in realized climatic niches that do not follow the pattern of mitochondrial divergence. These results lend support to the hypothesis that genetic divergence across Podarcis hispanicus complex most likely occurred in allopatric conditions, mostly with significant niche divergence. Competition after secondary contact is also suggested by the common occurrence of niche overlap between lineages that exhibit strictly parapatric distribution. The almost continuous distribution of Podarcis lizards in the study area appears to be a result of a combination of complementary suitable niches and competition, which seem two important mechanisms limiting geographic distributions and restricting the existence of extensive contact zones.     

Ecological correlates of the distribution limits of two poeciliid species along a salinity gradient

Identifying the environmental factors responsible for the formation of a species' distribution limit is challenging because organisms interact in complex ways with their environments. However, the use of statistical niche models in combination with the analysis of phenotypic variation along environmental gradients can help to reduce such complexity and identify a subset of candidate factors. In the present study, we used such approaches to describe and identify factors responsible for the parapatric distribution of two closely‐related livebearer fish species along a salinity gradient in the lowlands of Trinidad, West Indies. The downstream distribution limits of Poecilia reticulata were strongly correlated with the brackish–freshwater interface. We did not observe significant phenotypic variation in life‐history traits for this species when comparing marginal with more central populations, suggesting that abrupt changes in conditions at the brackish–freshwater interface limit its distribution. By contrast, Poecilia picta was present across a wide range of salinities, although it gradually disappeared from upstream freshwater localities. In addition, P. picta populations exhibited an increase in offspring size in localities where they coexist with P. reticulata, suggesting a role for interspecific competition. The parapatric distribution of these two species, suggests that P. reticulata distributions are limited by an abiotic factor (salinity), whereas P. picta is limited by a biotic factor (interspecific competition). Similar parapatric patterns have been previously described in other systems, suggesting they might be a common pattern in nature. © 2013 The Linnean Society of London     

ASYMMETRIC RESPONSE PATTERNS TO SUBSPECIES‐SPECIFIC SONG DIFFERENCES IN ALLOPATRY AND PARAPATRY IN THE GRAY‐BREASTED WOOD‐WREN

Song divergence between closely related taxa may play a critical role in the evolutionary processes of speciation and hybridization. We explored song variation between two Ecuadorian subspecies of the gray‐breasted wood‐wren (Henicorhina leucophrys) and tested the impact of song divergence on response behaviors. Songs were significantly different between the two subspecies, even between two parapatric populations 10 km apart. Playback experiments revealed an asymmetric response pattern to these divergent subspecies specific songs; one subspecies responded more to songs of its own subspecies than to the other subspecies’ songs, whereas the second responded equally strongly to songs of both subspecies. While song parameters revealed a mixed pattern of divergence between allopatric and parapatric populations, the majority of spectral characteristics showed increased divergence in parapatry, suggestive of character displacement. This increased song divergence in parapatry appeared to affect behavioral responses to playback as discriminating responses were most prominent in parapatry and against parapatric songs. The clear behavioral impact of subspecies‐specific song differences supports a potential role for song as an acoustic barrier to gene flow. The asymmetric nature of the responses suggests that song divergence could affect the direction of gene flow and the position of the subspecies‐specific transition.     

Diversification in North American arid lands: Niche conservatism, divergence and expansion of habitat explain speciation in the genus Ephedra

A lineage of 12 arid land shrubby species in the gymnosperm genus Ephedra (Gnetales) from North America is used to evaluate the influence of climate on speciation. With a long evolutionary history, and a well documented fossil record this lineage is an ideal model for understanding the process of speciation under a niche conservatism scenario. Using seven DNA molecular markers, Bayesian inference is carried out to uncover sister species and to estimate time of divergence of the lineages. Ecological niche models are generated for four parapatric and sympatric sister species and two analyses of niche evolution are performed, one based on ecological niche models and another using raw data and multivariate analysis. As previous analyses suggest, the diversification of North America Ephedra species may be the result of a recent secondary radiation. Both parapatric and sympatric species diverged mostly in a scenario of climatic niche conservatism. However, we also found strong evidence for niche divergence for one of the sister species pairs (E. californica-E. trifurca). Moreover, the multivariate analysis found environmental differences for some variables between sister species. The estimated divergence time of three pairs of sister species distributed in southwestern North America (E. cutleri-E. aspera, E. californica-E. trifurca and E. torreyana-E. viridis) is inferred to have occurred in the Late Miocene to Pliocene and for the sister species pair E. antisyphilitica-E. coryi distributed in the southern United States and northeastern Mexico, it was inferred from the Pliocene to Pleistocene. The orogenetic and climatic changes documented for these regions related to expansion of arid lands, may have contributed to the diversification in North American Ephedra, rather than adaptations to new climatic conditions.     

Evidence for spatial clines and mixed geographic modes of speciation for North American cherry‐infesting Rhagoletis (Diptera: Tephritidae) flies

An important criterion for understanding speciation is the geographic context of population divergence. Three major modes of allopatric, parapatric, and sympatric speciation define the extent of spatial overlap and gene flow between diverging populations. However, mixed modes of speciation are also possible, whereby populations experience periods of allopatry, parapatry, and/or sympatry at different times as they diverge. Here, we report clinal patterns of variation for 21 nuclear‐encoded microsatellites and a wing spot phenotype for cherry‐infesting Rhagoletis (Diptera: Tephritidae) across North America consistent with these flies having initially diverged in parapatry followed by a period of allopatric differentiation in the early Holocene. However, mitochondrial DNA (mtDNA) displays a different pattern; cherry flies at the ends of the clines in the eastern USA and Pacific Northwest share identical haplotypes, while centrally located populations in the southwestern USA and Mexico possess a different haplotype. We hypothesize that the mitochondrial difference could be due to lineage sorting but more likely reflects a selective sweep of a favorable mtDNA variant or the spread of an endosymbiont. The estimated divergence time for mtDNA suggests possible past allopatry, secondary contact, and subsequent isolation between USA and Mexican fly populations initiated before the Wisconsin glaciation. Thus, the current genetics of cherry flies may involve different mixed modes of divergence occurring in different portions of the fly''s range. We discuss the need for additional DNA sequencing and quantification of prezygotic and postzygotic reproductive isolation to verify the multiple mixed‐mode hypothesis for cherry flies and draw parallels from other systems to assess the generality that speciation may commonly involve complex biogeographies of varying combinations of allopatric, parapatric, and sympatric divergence.     

Evolutionary and demographic history of the Californian scrub white oak species complex: an integrative approach

Understanding the factors promoting species formation is a major task in evolutionary research. Here, we employ an integrative approach to study the evolutionary history of the Californian scrub white oak species complex (genus Quercus). To infer the relative importance of geographical isolation and ecological divergence in driving the speciation process, we (i) analysed inter‐ and intraspecific patterns of genetic differentiation and employed an approximate Bayesian computation (ABC) framework to evaluate different plausible scenarios of species divergence. In a second step, we (ii) linked the inferred divergence pathways with current and past species distribution models (SDMs) and (iii) tested for niche differentiation and phylogenetic niche conservatism across taxa. ABC analyses showed that the most plausible scenario is the one considering the divergence of two main lineages followed by a more recent pulse of speciation. Genotypic data in conjunction with SDMs and niche differentiation analyses support that different factors (geography vs. environment) and modes of speciation (parapatry, allopatry and maybe sympatry) have played a role in the divergence process within this complex. We found no significant relationship between genetic differentiation and niche overlap, which probably reflects niche lability and/or that multiple factors, have contributed to speciation. Our study shows that different mechanisms can drive divergence even among closely related taxa representing early stages of species formation and exemplifies the importance of adopting integrative approaches to get a better understanding of the speciation process.     

Divergence maintained by climatic selection despite recurrent gene flow: a case study of Castanopsis carlesii (Fagaceae)

Local adaptation to different environments has the potential to maintain divergence between populations despite recurrent gene flow and is an important driver for generating biological diversity. In this study, we investigate the role of adaptation in the maintenance of two parapatric varieties of a forest tree. We used sequence variation of chloroplastic DNA and restriction site‐associated DNA to investigate the genetic structure of two varieties of Castanopsis carlesii in subtropical China and relate it to climatic variation. We used niche reconstruction methods to investigate niche differentiation between the two varieties and to estimate the past distribution of this species. A deep divergence was observed between the two varieties, but evidence of introgression and genetic admixture was detected in two phenotypically and geographically intermediate populations. Niche reconstruction suggests that the distribution of the two varieties was disjunct during periods of global cooling and that the two varieties occupy significantly different niches. The genetic structure was mainly driven by environmental factors, and 13 outlier loci under divergent selection were correlated with climatic variation. These results suggest that the two varieties evolved in allopatry and came back into secondary contact after the last glacial maximum and that they are an evolutionary example of divergence maintained by climatic selection despite recurrent gene flow.     

Can tooth differentiation help to understand species coexistence? The case of wood mice in China

Five wood mice Apodemus species occur across China, in allopatry but also in sympatry up to cases of syntopy. They all share a similar external appearance, similar habitats of grasslands and forests and a generalist feeding behaviour. This overall similarity raises questions about the mechanisms insuring competition avoidance and allowing the coexistence of the species. In this context, a morphometric analysis of two characters related to feeding (mandible and molar) addressed the following issues: (1) Were the species actually different in size and/or shape of these characters, supporting their role in resource partitioning? (2) Did this pattern of phenotypic divergence match the neutral genetic differentiation, suggesting that differentiation might have occurred in a former phase of allopatry as a result of stochastic processes? (3) Did the species provide evidence of character displacement when occurring in sympatry, supporting an ongoing role of competition in the interspecific divergence? Results evidenced first that different traits, here mandibles and molars, provided discrepant pictures of the evolution of the Apodemus group in China. Mandible shape appeared as prone to vary in response to local conditions, blurring any phylogenetic or ecological pattern, whereas molar shape evolution appeared to be primarily driven by the degree of genetic differentiation. Molar size and shape segregated the different species in the morphospace, suggesting that these features may be involved in a resource partitioning between Apodemus species. The morphological segregation of the species, likely achieved by processes of differentiation in isolation promoted by the complex landscape of China, could contribute to competition avoidance and hence explain why no evidence was found of character displacement.     

Morphological differentiation correlates with ecological but not with genetic divergence in a Gehyra gecko

Body size affects life history, the ecological niche of an organism and its interactions with other organisms. Resultantly, marked differences in body size between related organisms are often an indication of a species boundary. This is particularly evident in the Gehyra variegata species complex of geckos, which displays differential body sizes between genetically divergent species, but high levels of intraspecific morphological conservatism. We report on a Gehyra population that displays extraordinary body size differentiation in comparison with other G. variegata species. We used morphological and environmental data to show this population is phenotypically and ecologically distinct from its parapatric congener Gehyra lazelli and that morphology and ecology are significantly correlated. Contrastingly, mtDNA analysis indicates paraphyly between the two groups, and allele frequencies at six microsatellite loci show no population structure concordant with morpho‐/ecotype. These results suggest either ecological speciation or environmentally induced phenotypic polymorphism, in an otherwise morphologically conservative group.     

Foraging niche dynamics and overlap in a guild of passerine birds in a south Swedish coniferous woodland

Summary Three foraging niche dimensions were measured at four seasons in six passerine bird species in a pine-dominated woodland. Total niche size increased in the order Certhia familiaris. In the first two species niche sizewas relatively constant at all seasons in spite of fluctuations in separatedimensions suggesting intraspecific complementarity between dimensionbreadths. Diffuse niche overlap was much lower in summer than at otherseasons in all species. Generally, this was due to movements of exploitationcurves along resource axes away from zones of interspecific overlap andto enlargement of total guild niche size, rather than to a reduction in totalniche size of the species. In some species a positive correlation was indicatedbetween breadth of niche dimension and diffuse overlap in the same dimension.Total niche size was correlated with total diffuse overlap up to acertain overlap value. Interspecific niche size complementarity was stronglysuggested by total guild niche being much less variable between seasons thanthe separate niche size of each species. The narrow niche of P. major, usuallyregarded a generalist, may be due to it being primarily adapted to broadleavedrather than coniferous woodland.     

Morphological and niche divergence of pinyon pines

The environmental variables that define a species ecological niche should be associated with the evolutionary patterns present in the adaptations that resulted from living in these conditions. Thus, when comparing across species, we can expect to find an association between phylogenetically independent phenotypic characters and ecological niche evolution. Few studies have evaluated how organismal phenotypes might mirror patterns of niche evolution if these phenotypes reflect adaptations. Doing so could contribute on the understanding of the origin and maintenance of phenotypic diversity observed in nature. Here, we show the pattern of niche evolution of the pinyon pine lineage (Pinus subsection Cembroides); then, we suggest morphological adaptations possibly related to niche divergence, and finally, we test for correlation between ecological niche and morphology. We demonstrate that niche divergence is the general pattern within the clade and that it is positively correlated with adaptation.     

Divergence and gene flow in the globally distributed blue‐winged ducks

The ability to disperse over long distances can result in a high propensity for colonizing new geographic regions, including uninhabited continents, and lead to lineage diversification via allopatric speciation. However, high vagility can also result in gene flow between otherwise allopatric populations, and in some cases, parapatric or divergence‐with‐gene‐flow models might be more applicable to widely distributed lineages. Here, we use five nuclear introns and the mitochondrial control region along with Bayesian models of isolation with migration to examine divergence, gene flow, and phylogenetic relationships within a cosmopolitan lineage comprising six species, the blue‐winged ducks (genus Anas), which inhabit all continents except Antarctica. We found two primary sub‐lineages, the globally‐distributed shoveler group and the New World blue‐winged/cinnamon teal group. The blue‐winged/cinnamon sub‐lineage is composed of sister taxa from North America and South America, and taxa with parapatric distributions are characterized by low to moderate levels of gene flow. In contrast, our data support strict allopatry for most comparisons within the shovelers. However, we found evidence of gene flow from the migratory, Holarctic northern shoveler A. clypeata and the more sedentary, African Cape shoveler A. smithii into the Australasian shoveler A. rhynchotis, although we could not reject strict allopatry. Given the diverse mechanisms of speciation within this complex, the shovelers and blue‐winged/cinnamon teals can serve as an effective model system for examining how the genome diverges under different evolutionary processes and how genetic variation is partitioned among highly dispersive taxa.     

Deep mitochondrial introgression and hybridization among ecologically divergent vole species

The completion of speciation is typically difficult to ascertain in rapidly diverging taxa but the amount of hybridization and gene flow in sympatry or parapatry contains important information about the level of reproductive isolation achieved. Here, we examined the progress in speciation between the Mediterranean (Microtus duodecimcostatus) and the Lusitanian pine vole (M. lusitanicus), which are part of the most rapid radiation of species known in mammals. These two Iberian pine voles are classified as separate species because of differences in morphology and ecology, but relatively many ambiguous individuals can be found in sympatric conditions. Our phylogenetic analyses of rangewide data from the mitochondrial cytochrome b gene (mtDNA) demonstrated high levels of diversity and a basal separation in two parapatric lineages. However, mtDNA affiliation was at odds with morphological classification or geographical distribution of the taxa. In contrast, statistical analyses of microsatellites (nucDNA) showed two clear genetic clusters in allopatry and sympatry generally matching morphological classification. This cytonuclear discordance over a large geographic area suggests historical introgression of mtDNA from M. duodecimcostatus to M. lusitanicus. There was statistical evidence for at least two recent hybrids in the sympatry zone but gene flow is apparently low given clear‐cut differences in nucDNA. Our results indicate a relatively advanced speciation process in these Iberian pine voles without fully established reproductive isolation. This situation enables use of combined population genomic and experimental approaches for the separation of patterns and mechanisms in the ongoing explosive diversification of these and other Arvicoline rodents in the future.