Speciation in the Rana chensinensis species complex and its relationship to the uplift of the Qinghai-Tibetan Plateau

Speciation remains a fundamental issue in biology. Herein, we report an investigation into speciation in the Rana chensinensis species complex using DNA sequence data from one mitochondrial and five nuclear genes. A phylogenetic analysis of the data revealed four major clades in the complex, and each of them was found to likely represent a species, including one cryptic species. Ecological niche models were generated from 19 climatic variables for three of the four major clades, which were represented by widespread sampling, including R. chensinensis, Rana kukunoris and the potential cryptic species. Each clade is associated with a unique ecological unit, and this indicates that ecological divergence probably drove speciation. Ecological divergence is likely related to the late Cenozoic orogenesis of the Qinghai–Tibetan Plateau. In addition, gene flow between species was detected but only in peripheral portions of the ranges of the four major clades, thus likely had little influence on the speciation processes. Discordances between mitochondrial and nuclear genes were also found; the nominal species, R. chensinensis, contains multiple maternal clades, suggesting potential mitochondrial introgression between R. chensinensis and R. kukunoris.     

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.      

Genomic signatures of divergent selection and speciation patterns in a ‘natural experiment’, the young parallel radiations of Nicaraguan crater lake cichlid fishes

Divergent selection is the main driving force in sympatric ecological speciation and may also play a strong role in divergence between allopatric populations. Characterizing the genome‐wide impact of divergent selection often constitutes a first step in unravelling the genetic bases underlying adaptation and ecological speciation. The Midas cichlid fish (Amphilophus citrinellus) species complex in Nicaragua is a powerful system for studying evolutionary processes. Independent colonizations of isolated young crater lakes by Midas cichlid populations from the older and great lakes of Nicaragua resulted in the repeated evolution of adaptive radiations by intralacustrine sympatric speciation. In this study we performed genome scans on two repeated radiations of crater lake species and their great lake source populations (1030 polymorphic AFLPs, n ～ 30 individuals per species). We detected regions under divergent selection (0.3% in the crater lake Xiloá flock and 1.7% in the older crater lake Apoyo radiation) that might be responsible for the sympatric diversifications. We find no evidence that the same genomic regions have been involved in the repeated evolution of parallel adaptations across crater lake flocks. However, there is some genetic parallelism apparent (seven out of 51 crater lake to great lake outlier loci are shared; 13.7%) that is associated with the allopatric divergence of both crater lake flocks. Interestingly, our results suggest that the number of outlier loci involved in sympatric and allopatric divergence increases over time. A phylogeny based on the AFLP data clearly supports the monophyly of both crater lake species flocks and indicates a parallel branching order with a primary split along the limnetic‐benthic axis in both radiations.     

Mutation‐order divergence by sexual selection: diversification of sexual signals in similar environments as a first step in speciation

The origin of species remains a central question, and recent research focuses on the role of ecological differences in promoting speciation. Ecological differences create opportunities for divergent selection (i.e. ‘ecological’ speciation), a Darwinian hypothesis that hardly requires justification. In contrast, ‘mutation‐order’ speciation proposes that, instead of adapting to different environments, populations find different ways to adapt to similar environments, implying that speciation does not require ecological differences. This distinction is critical as it provides an alternative hypothesis to the prevailing view that ecological differences drive speciation. Speciation by sexual selection lies at the centre of debates about the importance of ecological differences in promoting speciation; here, we present verbal and mathematical models of mutation‐order divergence by sexual selection. We develop three general cases and provide a two‐locus population genetic model for each. Results indicate that alternative secondary sexual traits can fix in populations that initially experience similar natural and sexual selection and that divergent traits and preferences can remain stable in the face of low gene flow. This stable divergence can facilitate subsequent divergence that completes or reinforces speciation. We argue that a mutation‐order process could explain widespread diversity in secondary sexual traits among closely related, allopatric species.     

Ecological divergence combined with ancient allopatry in lizard populations from a small volcanic island

Population divergence and speciation are often explained by geographical isolation, but may also be possible under high gene flow due to strong ecology‐related differences in selection pressures. This study combines coalescent analyses of genetic data (11 microsatellite loci and 1 Kbp of mtDNA) and ecological modelling to examine the relative contributions of isolation and ecology to incipient speciation in the scincid lizard Chalcides sexlineatus within the volcanic island of Gran Canaria. Bayesian multispecies coalescent dating of within‐island genetic divergence of northern and southern populations showed correspondence with the timing of volcanic activity in the north of the island 1.5–3.0 Ma ago. Coalescent estimates of demographic changes reveal historical size increases in northern populations, consistent with expansions from a volcanic refuge. Nevertheless, ecological divergence is also supported. First, the two morphs showed non‐equivalence of ecological niches and species distribution modelling associated the northern morph with mesic habitat types and the southern morph with xeric habitat types. It seems likely that the colour morphs are associated with different antipredator strategies in the different habitats. Second, coalescent estimation of gene copy migration (based on microsatellites and mtDNA) suggest high rates from northern to southern morphs demonstrating the strength of ecology‐mediated selection pressures that maintain the divergent southern morph. Together, these findings underline the complexity of the speciation process by providing evidence for the combined effects of ecological divergence and ancient divergence in allopatry.     

Host‐driven divergence in the parasitic plant Orobanche minor Sm. (Orobanchaceae)

Many parasitic angiosperms have a broad host range and are therefore considered to be host generalists. Orobanche minor is a nonphotosynthetic root parasite that attacks a range of hosts from taxonomically disparate families. In the present study, we show that O. minor sensu lato may comprise distinct, genetically divergent races isolated by the different ecologies of their hosts. Using a three‐pronged approach, we tested the hypothesis that intraspecific taxa O. minor var. minor and O. minor ssp. maritima parasitizing either clover (Trifolium pratense) or sea carrot (Daucus carota ssp. gummifer), respectively, are in allopatric isolation. Morphometric analysis revealed evidence of divergence but this was insufficient to define discrete, host‐specific taxa. Intersimple sequence repeat (ISSR) marker‐based data provided stronger evidence of divergence, suggesting that populations were isolated from gene flow. Phylogenetic analysis, using sequence‐characterized amplified region (SCAR) markers derived from ISSR loci, provided strong evidence for divergence by clearly differentiating sea carrot‐specific clades and mixed‐host clades. Low levels of intrapopulation SCAR marker sequence variation and floral morphology suggest that populations on different hosts are probably selfing and inbreeding. Morphologically cryptic Orobanche taxa may therefore be isolated from gene flow by host ecology. Together, these data suggest that host specificity may be an important driver of allopatric speciation in parasitic plants.     

THE GENETIC ARCHITECTURE OF REPRODUCTIVE ISOLATION DURING SPECIATION‐WITH‐GENE‐FLOW IN LAKE WHITEFISH SPECIES PAIRS ASSESSED BY RAD SEQUENCING

During speciation‐with‐gene‐flow, effective migration varies across the genome as a function of several factors, including proximity of selected loci, recombination rate, strength of selection, and number of selected loci. Genome scans may provide better empirical understanding of the genome‐wide patterns of genetic differentiation, especially if the variance due to the previously mentioned factors is partitioned. In North American lake whitefish (Coregonus clupeaformis), glacial lineages that diverged in allopatry about 60,000 years ago and came into contact 12,000 years ago have independently evolved in several lakes into two sympatric species pairs (a normal benthic and a dwarf limnetic). Variable degrees of reproductive isolation between species pairs across lakes offer a continuum of genetic and phenotypic divergence associated with adaptation to distinct ecological niches. To disentangle the complex array of genetically based barriers that locally reduce the effective migration rate between whitefish species pairs, we compared genome‐wide patterns of divergence across five lakes distributed along this divergence continuum. Using restriction site associated DNA (RAD) sequencing, we combined genetic mapping and population genetics approaches to identify genomic regions resistant to introgression and derive empirical measures of the barrier strength as a function of recombination distance. We found that the size of the genomic islands of differentiation was influenced by the joint effects of linkage disequilibrium maintained by selection on many loci, the strength of ecological niche divergence, as well as demographic characteristics unique to each lake. Partial parallelism in divergent genomic regions likely reflected the combined effects of polygenic adaptation from standing variation and independent changes in the genetic architecture of postzygotic isolation. This study illustrates how integrating genetic mapping and population genomics of multiple sympatric species pairs provide a window on the speciation‐with‐gene‐flow mechanism.     

Adaptive divergence between lake and stream populations of an East African cichlid fish

Divergent natural selection acting in different habitats may build up barriers to gene flow and initiate speciation. This speciation continuum can range from weak or no divergence to strong genetic differentiation between populations. Here, we focus on the early phases of adaptive divergence in the East African cichlid fish Astatotilapia burtoni, which occurs in both Lake Tanganyika (LT) and inflowing rivers. We first assessed the population structure and morphological differences in A. burtoni from southern LT. We then focused on four lake–stream systems and quantified body shape, ecologically relevant traits (gill raker and lower pharyngeal jaw) as well as stomach contents. Our study revealed the presence of several divergent lake–stream populations that rest at different stages of the speciation continuum, but show the same morphological and ecological trajectories along the lake–stream gradient. Lake fish have higher bodies, a more superior mouth position, longer gill rakers and more slender pharyngeal jaws, and they show a plant/algae and zooplankton‐biased diet, whereas stream fish feed more on snails, insects and plant seeds. A test for reproductive isolation between closely related lake and stream populations did not detect population‐assortative mating. Analyses of F1 offspring reared under common garden conditions indicate that the detected differences in body shape and gill raker length do not constitute pure plastic responses to different environmental conditions, but also have a genetic basis. Taken together, the A. burtoni lake–stream system constitutes a new model to study the factors that enhance and constrain progress towards speciation in cichlid fishes.     

Assessing models of speciation under different biogeographic scenarios; an empirical study using multi‐locus and RNA‐seq analyses

Evolutionary biology often seeks to decipher the drivers of speciation, and much debate persists over the relative importance of isolation and gene flow in the formation of new species. Genetic studies of closely related species can assess if gene flow was present during speciation, because signatures of past introgression often persist in the genome. We test hypotheses on which mechanisms of speciation drove diversity among three distinct lineages of desert tortoise in the genus Gopherus. These lineages offer a powerful system to study speciation, because different biogeographic patterns (physical vs. ecological segregation) are observed at opposing ends of their distributions. We use 82 samples collected from 38 sites, representing the entire species' distribution and generate sequence data for mtDNA and four nuclear loci. A multilocus phylogenetic analysis in *BEAST estimates the species tree. RNA‐seq data yield 20,126 synonymous variants from 7665 contigs from two individuals of each of the three lineages. Analyses of these data using the demographic inference package ?a?i serve to test the null hypothesis of no gene flow during divergence. The best‐fit demographic model for the three taxa is concordant with the *BEAST species tree, and the ?a?i analysis does not indicate gene flow among any of the three lineages during their divergence. These analyses suggest that divergence among the lineages occurred in the absence of gene flow and in this scenario the genetic signature of ecological isolation (parapatric model) cannot be differentiated from geographic isolation (allopatric model).     

Gene expression remodelling and immune response during adaptive divergence in an African cichlid fish

Variation in gene expression contributes to ecological speciation by facilitating population persistence in novel environments. Likewise, immune responses can be of relevance in speciation driven by adaptation to different environments. Previous studies examining gene expression differences between recently diverged ecotypes have often relied on only one pair of populations, targeted the expression of only a subset of genes or used wild‐caught individuals. Here, we investigated the contribution of habitat‐specific parasites and symbionts and the underlying immunological abilities of ecotype hosts to adaptive divergence in lake–river population pairs of the cichlid fish Astatotilapia burtoni. To shed light on the role of phenotypic plasticity in adaptive divergence, we compared parasite and microbiota communities, immune response, and gene expression patterns of fish from natural habitats and a lake‐like pond set‐up. In all investigated population pairs, lake fish were more heavily parasitized than river fish, in terms of both parasite taxon composition and infection abundance. The innate immune response in the wild was higher in lake than in river populations and was elevated in a river population exposed to lake parasites in the pond set‐up. Environmental differences between lake and river habitat and their distinct parasite communities have shaped differential gene expression, involving genes functioning in osmoregulation and immune response. Most changes in gene expression between lake and river samples in the wild and in the pond set‐up were based on a plastic response. Finally, gene expression and bacterial communities of wild‐caught individuals and individuals acclimatized to lake‐like pond conditions showed shifts underlying adaptive phenotypic plasticity.     

Physiographic and climatic events in the Chihuahuan Desert lead to the speciation and distinct demographic patterns of two sister Sceloporus lizards

The biogeographic history of the Chihuahuan Desert is known to be complex, and there is evidence of the effects of physiographic and climatic events in species diversification and demographic population changes in many taxa. Here, using DNA sequence data, we studied the influence of the physiographic and climatic events that occurred in the Chihuahuan Desert during the Pliocene–Pleistocene transition on the speciation and evolutionary history of the sister lizard species Sceloporus cyanostictus and S. gadsdeni. First, based on mtDNA and nDNA sequences, we estimated the divergence times of the sister species. Then, based on mtDNA sequences, we investigated the demographic history of both species within a phylogeographic framework. The divergence time was inferred to be 1.48 Mya, date that corresponds to the existence of a large lake in the Mapimian subprovince, between the current geographic locations of S. cyanostictus and S. gadsdeni. This lake could have acted as a barrier, leading to the speciation of both species. For the demographic history of the two species, we identified two distinct patterns: the population expansion of S. gadsdeni within the Last Glacial Maximum and the potential population decline of S. cyanostictus. Our results can be used as a guide for the study of other aspects that could be critical to developing conservation actions that ensure the survival of not only S. gadsdeni and S. cyanostictus, but also other co‐occurring lizard species.     

Genomewide patterns of variation in genetic diversity are shared among populations,species and higher‐order taxa

Genomewide screens of genetic variation within and between populations can reveal signatures of selection implicated in adaptation and speciation. Genomic regions with low genetic diversity and elevated differentiation reflective of locally reduced effective population sizes (N_e) are candidates for barrier loci contributing to population divergence. Yet, such candidate genomic regions need not arise as a result of selection promoting adaptation or advancing reproductive isolation. Linked selection unrelated to lineage‐specific adaptation or population divergence can generate comparable signatures. It is challenging to distinguish between these processes, particularly when diverging populations share ancestral genetic variation. In this study, we took a comparative approach using population assemblages from distant clades assessing genomic parallelism of variation in N_e. Utilizing population‐level polymorphism data from 444 resequenced genomes of three avian clades spanning 50 million years of evolution, we tested whether population genetic summary statistics reflecting genomewide variation in N_e would covary among populations within clades, and importantly, also among clades where lineage sorting has been completed. All statistics including population‐scaled recombination rate (ρ), nucleotide diversity (π) and measures of genetic differentiation between populations (F_ST, PBS, d_xy) were significantly correlated across all phylogenetic distances. Moreover, genomic regions with elevated levels of genetic differentiation were associated with inferred pericentromeric and subtelomeric regions. The phylogenetic stability of diversity landscapes and stable association with genomic features support a role of linked selection not necessarily associated with adaptation and speciation in shaping patterns of genomewide heterogeneity in genetic diversity.     

Phylogeographic diversification of antelope squirrels (Ammospermophilus) across North American deserts

We investigated the biogeographic history of antelope squirrels, genus Ammospermophilus, which are widely distributed across the deserts and other arid lands of western North America. We combined range‐wide sampling of all currently recognized species of Ammospermophilus with a multilocus data set to infer phylogenetic relationships. We then estimated divergence times within identified clades of Ammospermophilus using fossil‐calibrated and rate‐calibrated molecular clocks. Lastly, we explored generalized distributional changes of Ammospermophilus since the last glacial maximum using species distribution models, and assessed responses to Quaternary climate change by generating demographic parameter estimates for the three wide‐ranging clades of A. leucurus. From our phylogenetic estimates we inferred strong phylogeographic structure within Ammospermophilus and the presence of three well‐supported major clades. Initial patterns of historical divergence were coincident with dynamic alterations in the landscape of western North America, and the formation of regional deserts during the Late Miocene and Pliocene. Species distribution models and demographic parameter estimates revealed patterns of recent population expansion in response to glacial retreat. When combined with evidence from co‐distributed taxa, the historical biogeography of Ammospermophilus provides additional insight into the mechanisms that impacted diversification of arid‐adapted taxa across the arid lands of western North America. We propose species recognition of populations of the southern Baja California peninsula to best represent our current understanding of evolutionary relationships among genetic units of Ammospermophilus. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 109 , 949–967.     

Differential divergence in autosomes and sex chromosomes is associated with intra‐island diversification at a very small spatial scale in a songbird lineage

Recently diverged taxa showing marked phenotypic and ecological diversity provide optimal systems to understand the genetic processes underlying speciation. We used genome‐wide markers to investigate the diversification of the Reunion grey white‐eye (Zosterops borbonicus) on the small volcanic island of Reunion (Mascarene archipelago), where this species complex exhibits four geographical forms that are parapatrically distributed across the island and differ strikingly in plumage colour. One form restricted to the highlands is separated by a steep ecological gradient from three distinct lowland forms which meet at narrow hybrid zones that are not associated with environmental variables. Analyses of genomic variation based on single nucleotide polymorphism data from genotyping‐by‐sequencing and pooled RAD‐seq approaches show that signatures of selection associated with elevation can be found at multiple regions across the genome, whereas most loci associated with the lowland forms are located on the Z sex chromosome. We identified TYRP1, a Z‐linked colour gene, as a likely candidate locus underlying colour variation among lowland forms. Tests of demographic models revealed that highland and lowland forms diverged in the presence of gene flow, and divergence has progressed as gene flow was restricted by selection at loci across the genome. This system holds promise for investigating how adaptation and reproductive isolation shape the genomic landscape of divergence at multiple stages of the speciation process.     

Calls,colours, shape,and genes: a multi‐trait approach to the study of geographic variation in the Amazonian frog Allobates femoralis

Evolutionary divergence in behavioural traits related to mating may represent the initial stage of speciation. Direct selective forces are usually invoked to explain divergence in mate‐recognition traits, often neglecting a role for neutral processes or concomitant differentiation in ecological traits. We adopted a multi‐trait approach to obtain a deeper understanding of the mechanisms behind allopatric divergence in the Amazonian frog, Allobates femoralis. We tested the null hypothesis that geographic distance between populations correlates with genetic and phenotypic divergence, and compared divergence between mate‐recognition (acoustic) and ecological (coloration, body‐shape) traits. We quantified geographic variation in 39 phenotypic traits and a mitochondrial DNA marker among 125 individuals representing eight populations. Geographic variation in acoustic traits was pronounced and tracked the spatial genetic variation, which appeared to be neutral. Thus, the evolution of acoustic traits tracked the shared history of the populations, which is unexpected for pan‐Amazonian taxa or for mate‐recognition traits. Divergence in coloration appeared uncorrelated with genetic distance, and might be partly attributed to local selective pressures, and perhaps to Batesian mimicry. Divergence in body‐shape traits was low. The results obtained depict a complex evolutionary scenario and emphasize the importance of considering multiple traits when disentangling the forces behind allopatric divergence. ©2009 The Linnean Society of London, Biological Journal of the Linnean Society, 2009, 98 , 826–838.     

Eco‐morphological differentiation in Lake Magadi tilapia,an extremophile cichlid fish living in hot,alkaline and hypersaline lakes in East Africa

Ecological diversification through divergent selection is thought to be a major force during the process of adaptive radiations. However, the large sizes and complexity of most radiations such as those of the cichlids in the African Great Lakes make it impossible to infer the exact evolutionary history of any population divergence event. The genus Alcolapia, a small cichlid lineage endemic to Lakes Magadi and Natron in East Africa, exhibits phenotypes similar to some of those found in cichlids of the radiations of the African Great Lakes. The simplicity within Alcolapia makes it an excellent model system to investigate ecological diversification and speciation. We used an integrated approach including population genomics based on RAD‐seq data, geometric morphometrics and stable isotope analyses to investigate the eco‐morphological diversification of tilapia in Lake Magadi and its satellite lake Little Magadi. Additionally, we reconstructed the demographic history of the species using coalescent simulations based on the joint site frequency spectrum. The population in Little Magadi has a characteristically upturned mouth—possibly an adaptation to feeding on prey from the water surface. Eco‐morphological differences between populations within Lake Magadi are more subtle, but are consistent with known ecological differences between its lagoons such as high concentrations of nitrogen attributable to extensive guano deposits in Rest of Magadi relative to Fish Springs Lagoon. All populations diverged simultaneously only about 1100 generations ago. Differences in levels of gene flow between populations and the effective population sizes have likely resulted in the inferred heterogeneous patterns of genome‐wide differentiation.     

Divergence along a steep ecological gradient in lake whitefish (Coregonus sp.)

To understand mechanisms structuring diversity in young adaptive radiations, quantitative and unbiased information about genetic and phenotypic diversity is much needed. Here, we present the first in‐depth investigation of whitefish diversity in a Swiss lake, with continuous spawning habitat sampling in both time and space. Our results show a clear cline like pattern in genetics and morphology of populations sampled along an ecological depth gradient in Lake Neuchâtel. Divergent natural selection appears to be involved in shaping this cline given that trait specific P_ST‐values are significantly higher than F_ST‐values when comparing populations caught at different depths. These differences also tend to increase with increasing differences in depth, indicating adaptive divergence along a depth gradient, which persists despite considerable gene flow between adjacent demes. It however remains unclear, whether the observed pattern is a result of currently stable selection‐gene flow balance, incipient speciation, or reverse speciation due to anthropogenic habitat alteration causing two formerly divergent species to collapse into a single gene pool.     

THE EVOLUTIONARY HISTORY OF DARWIN'S FINCHES: SPECIATION,GENE FLOW,AND INTROGRESSION IN A FRAGMENTED LANDSCAPE

Many classic examples of adaptive radiations take place within fragmented systems such as islands or mountains, but the roles of mosaic landscapes and variable gene flow in facilitating species diversification is poorly understood. Here we combine phylogenetic and landscape genetic approaches to understand diversification in Darwin's finches, a model adaptive radiation. We combined sequence data from 14 nuclear introns, mitochondrial markers, and microsatellite variation from 51 populations of all 15 recognized species. Phylogenetic species‐trees recovered seven major finch clades: ground, tree, vegetarian, Cocos Island, grey and green warbler finches, and a distinct clade of sharp‐beaked ground finches (Geospiza cf. difficilis) basal to all ground and tree finches. The ground and tree finch clades lack species‐level phylogenetic structure. Interisland gene flow and interspecies introgression vary geographically in predictable ways. First, several species exhibit concordant patterns of population divergence across the channel separating the Galápagos platform islands from the separate volcanic province of northern islands. Second, peripheral islands have more admixed populations while central islands maintain more distinct species boundaries. This landscape perspective highlights a likely role for isolation of peripheral populations in initial divergence, and demonstrates that peripheral populations may maintain genetic diversity through outbreeding during the initial stages of speciation.     

Here be dragons: a phylogenetic and biogeographical study of the Smaug warreni species complex (Squamata: Cordylidae) in southern Africa

Taxonomy of the Smaug warreni species complex remains contentious despite known morphological differences and geographical separation of the various taxa. This study uses an 11‐gene dataset to recover phylogenetic relationships between the seven nominal members of the S. warreni complex. Eight well‐supported clades were returned, with S. warreni barbertonensis found to be paraphyletic. A time‐calibrated analysis of molecular data indicates that all eight clades in the S. warreni complex separated in the late Miocene, much earlier than the date suggested by the existing hypothesis of vicariance through the ingression of Kalahari sands. Ecological niche modelling indicates that although all clades are allopatric, a slight decrease in temperature could potentially render them sympatric, supporting an hypothesis of range expansion through climatic change. © 2014 The Linnean Society of London