Cryptic divergent lineages of Pultenaea pauciflora M.B. Scott (Fabaceae: Mirbelieae) exhibit different evolutionary history

Genetic structure among disjunct population groups of Pultenaea pauciflora was assessed to determine the evolutionary history of this species as a basis for conservation management strategies. Analysis of individuals from all extant populations using 1737 amplified length polymorphism markers revealed two highly divergent genetic entities with strong geographical structuring. Populations located at Narrogin and Brookton clustered together in Bayesian assignment analysis with every individual optimally placed in a single cluster with complete membership. Genetic differentiation between populations in these two areas was very low. Populations at Boddington were highly divergent from those located at Narrogin and Brookton. All individuals from Boddington populations were optimally placed into a second cluster with complete membership. Populations located at Boddington maintain lower levels of allelic diversity, yet greater levels of mean heterozygosity than populations located at Narrogin and Brookton. The degree of genetic differentiation and different patterns of genetic diversity strongly suggest historical divergence and separate evolutionary influences on the two lineages that occur in different ecological habitat. These Evolutionary Significant Units are likely to represent two cryptic sister taxa in the extant populations currently recognized as P. pauciflora, and the reassessment of taxonomic and conservation status of both lineages is required. © 2013 State of Western Australia. Biological Journal of the Linnean Society © 2013 The Linnean Society of London     

Identifying evolutionarily significant units for conservation of the endangered Malus sieversii using genome‐wide RADseq data

Malus sieversii, a wild progenitor of the domesticated apple, is an endangered species and is assigned second conservation priority by the China Plant Red Data Book. It is urgent to carry out in situ conservation of this species, but previous studies have not identified evolutionarily significant units (ESUs) for conservation management. In this study, we investigated the genetic diversity and relationships of six M. sieversii populations from China using integrated analysis of microsatellite (nSSR) data, genome‐wide SNPs and previous results in order to propose a reasonable conservation management. The results showed that levels of genetic diversity were inconsistently reflected by our nSSR and previous studies, suggesting that indices of genetic diversity are not effective to identify priority conservation areas for M. sieversii. Based on the selection criteria of ESUs for endangered species conservation, ESUs should reflect lineage divergence, geographical separation and different adaptive variation. Our phylogenetic tree based on genome‐wide SNPs yielded a clear relationship of divergent lineages among M. sieversii populations, leading to new different from those of previous studies. Three independent lineages, including the pairs of populations Huocheng‐Yining, Gongliu‐Xinyuan and Tuoli‐Emin, were identified. The geographic distances between populations among the different phylogenetic lineages were much greater than those within the same phylogenetic lineage. A cluster analysis on environmental variables showed that the three independent lineages inhabit different environmental conditions, suggesting that they may have adapted to different environments. Based on the results, we propose that three independent ESUs should be recognized as conservation units for M. sieversii in China.     

Ecological and genetic evidence for cryptic ecotypes in a rare sexually deceptive orchid,Drakaea elastica

Species with specialized ecological interactions present significant conservation challenges. In plants that attract pollinators with pollinator‐specific chemical signals, geographical variation in pollinator species may indicate the presence of cryptic plant taxa. We investigated this phenomenon in the rare sexually deceptive orchid Drakaea elastica using a molecular phylogenetic analysis to resolve pollinator species boundaries, pollinator choice experiments and a population genetic study of the orchid. Pollinator choice experiments demonstrated the existence of two ecotypes within D. elastica, each attracting their own related but phylogenetically distinct pollinator species. Despite the presence of ecotypes, population genetic differentiation was low across populations at six microsatellite loci (F_ST = 0.026). However, Bayesian STRUCTURE analysis revealed two genetic clusters, broadly congruent with the ecotype distributions. These ecotypes may represent adaptation to regional variation in pollinator availability and perhaps the early stages of speciation, with pronounced morphological and genetic differences yet to evolve. Resolution of the taxonomic status of the D. elastica ecotypes is required as this has implications for conservation efforts and allocation of management funding. Furthermore, any reintroduction programmes must incorporate knowledge of ecotype distribution and pollinator availability to ensure reproductive success in restored populations. © 2014 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 177 , 124–140.     

Estimating divergence time for two evolutionarily significant units of a protected fish species

The evolutionarily significant unit concept provides a powerful tool for conserving biodiversity below the species level, but temporal criteria are often used explicitly or implicitly in the operational definitions of evolutionarily significant units (ESUs). Such temporal considerations have important implications for recently diverged taxa, as is the case with the White sands pupfish (Cyprinodon tularosa). This species consists of two native populations previously designated as the Malpais Spring and Salt Creek ESUs based on allele frequency differences at nuclear markers and their ecologically divergent habitats; despite a lack of reciprocal monophyly. Isolation of these two ESUs presumably occurred during the mid-Holocene, but an alternative hypothesis is that the populations were isolated due to changes in surface hydrology associated with overgrazing in the late 19th and early 20th centuries. We assayed 13 microsatellite loci and applied an Approximate Bayesian Computation analysis to estimate time of divergence between the two populations. Our reference table consisted of 1,000,000 simulated data sets, and we used three different models, each having different combinations of summary statistics. Estimates of median divergence time varied from approximately 6,500–11,000 generations (3,250–11,000 years). These findings support the hypothesis that Malpais Spring and Salt Creek having been isolated for a least a few millennia, and together with previously documented adaptive divergence, argues for continued management as separate conservation units. We consider the temporal constraints for defining evolutionary significance as it relates to recently diverged populations occupying ecologically divergent habitats.     

Fish conservation in the land of steppe and sky: Evolutionarily significant units of threatened salmonid species in Mongolia mirror major river basins

Mongolia's salmonids are suffering extensive population declines; thus, more comprehensive fisheries management and conservation strategies are required. To assist with their development, a better understanding of the genetic structure and diversity of these threatened species would allow a more targeted approach for preserving genetic variation and ultimately improve long‐term species recoveries. It is hypothesized that the unfragmented river basins that have persisted across Mongolia provide unobstructed connectivity for resident salmonid species. Thus, genetic structure is expected to be primarily segregated between major river basins. We tested this hypothesis by investigating the population structure for three salmonid genera (Hucho, Brachymystax and Thymallus) using different genetic markers to identify evolutionarily significant units (ESUs) and priority rivers to focus conservation efforts. Fish were assigned to separate ESUs when the combined evidence of mitochondrial and nuclear data indicated genetic isolation. Hucho taimen exhibited a dichotomous population structure forming two ESUs, with five priority rivers. Within the Brachymystax genus, there were three B. lenokESUs and one B. tumensisESU, along with six priority rivers. While B. tumensiswas confirmed to display divergent mtDNA haplotypes, haplotype sharing between these two congeneric species was also identified. For T. baicalensis,only a single ESU was assigned, with five priority rivers identified plus Lake Hovsgol. Additionally, we confirmed that T. nigrescens from Lake Hovsgol is a synonym of T. baicalensis. Across all species, the most prominent pattern was strong differentiation among major river basins with low differentiation and weak patterns of isolation by distance within river basins, which corroborated our hypothesis of high within‐basin connectivity across Mongolia. This new genetic information provides authorities the opportunity to distribute resources for management between ESUs while assigning additional protection for the more genetically valuable salmonid rivers so that the greatest adaptive potential within each species can be preserved.     

Microsatellite genetic variation in the Chinese endemic Eucommia ulmoides (Eucommiaceae): implications for conservation

Eucommia ulmoides, a traditional Chinese medicinal plant, is endangered as a consequence of long‐term and widespread harvest in the late 20th century. It has been widely cultivated as a source of herbal medicine and for use in the organic chemical industry in China. In this study, eight microsatellite markers were applied to investigate genetic diversity in E. ulmoides. Three hundred individuals from one semi‐wild population and nine cultivated populations across its main production area were collected. A high level of genetic diversity at population levels (H_E = 0.716) was observed. The highly outcrossed mating system, high longevity of E. ulmoides and seed admixture may be responsible for high genetic variation within populations. A genetic bottleneck was observed in one population. Populations were only slightly differentiated from one another (F_ST = 0.063); this was also supported by AMOVA, which revealed that 94.05% of the total variation resided within populations. This is probably attributable to long‐distance gene flow mediated by the exchange of seeds by local farmers. Implications of these results for the conservation of genetic resources of E. ulmoides are discussed. © 2013 The Linnean Society of London, Botanical Journal of the Linnean Society, 2013, 173 , 775–785.     

Resurrecting an extinct salmon evolutionarily significant unit: archived scales, historical DNA and implications for restoration

Archival scales from 603 sockeye salmon (Oncorhynchus nerka), sampled from May to July 1924 in the lower Columbia River, were analysed for genetic variability at 12 microsatellite loci and compared to 17 present‐day O. nerka populations—exhibiting either anadromous (sockeye salmon) or nonanadromous (kokanee) life histories—from throughout the Columbia River Basin, including areas upstream of impassable dams built subsequent to 1924. Statistical analyses identified four major genetic assemblages of sockeye salmon in the 1924 samples. Two of these putative historical groupings were found to be genetically similar to extant evolutionarily significant units (ESUs) in the Okanogan and Wenatchee Rivers (pairwise F_ST = 0.004 and 0.002, respectively), and assignment tests were able to allocate 77% of the fish in these two historical groupings to the contemporary Okanogan River and Lake Wenatchee ESUs. A third historical genetic grouping was most closely aligned with contemporary sockeye salmon in Redfish Lake, Idaho, although the association was less robust (pairwise F_ST = 0.060). However, a fourth genetic grouping did not appear to be related to any contemporary sockeye salmon or kokanee population, assigned poorly to the O. nerka baseline, and had distinctive early return migration timing, suggesting that this group represents a historical ESU originating in headwater lakes in British Columbia that was probably extirpated sometime after 1924. The lack of a contemporary O. nerka population possessing the genetic legacy of this extinct ESU indicates that efforts to reestablish early‐migrating sockeye salmon to the headwater lakes region of the Columbia River will be difficult.     

Selection and localised genetic structure in the threatened Manauense Harlequin Frog (Bufonidae: Atelopus manauensis)

Knowledge of genome-wide variation and the processes influencing gene flow are critical to managing threatened species. Here, we characterise genetic diversity and the environmental features associated with connectivity for a narrowly distributed and threatened Amazonian frog, Atelopus manauensis. We sampled 94 individuals throughout the upper, middle and lower courses from each of the five major rivers covering the species’ known geographic range and genotyped each individual at 3859 single nucleotide polymorphisms (SNPs). Genetic variation was significantly subdivided into six groups, each mostly containing individuals sampled within the same major river system. The genetic distances among these groups increased with geographic distance, and open forests and the extremities of the altitude gradient were associated with less genetic connectivity. Using F_ST outlier approaches and environmental association analyses, we identified SNPs indicative of localised adaptation, with 28 SNPs significantly associated with forest biomass and altitude. Evidence of divergent selection among the six genetic clusters suggests the presence of six Evolutionary Significant Units (ESUs). Overall, the ESUs were characterized by low contemporary effective population sizes (N_E?<?100) suggesting that genetic variation will be lost by random genetic drift. We demonstrate surprisingly high levels of divergence across the limited distribution of A. manauensis and suggest that each of the six adaptively divergent lineages be considered in conservation planning.

     

Molecular and morphometric data demonstrate the presence of cryptic lineages among freshwater crabs (Decapoda: Potamonautidae: Potamonautes) from the Drakensberg Mountains, South Africa

The evolutionary importance of cryptic taxa is well documented, yet few studies have examined them with a view to conservation. In the present study, the significance of cryptic speciation in freshwater crabs is examined. Allozyme and 16S rRNA sequence data were used to explore the degree of population differentiation between mountain stream populations of two distinct freshwater crab species. Marked patterns of differentiation were evident among populations; those in close geographical proximity were characterized by the fixation of alternate alleles at certain loci, indicating that currently there is no gene flow among populations. Both allozyme and sequence data provide evidence for the recognition of at least five distinct evolutionary lineages with pronounced levels of genetic differentiation. Morphometric data failed to detect any distinct geographically intermediate population groupings. Our findings indicate the presence of five phylogeographical units, all worthy of conservation, three of which are evolutionarily significant.  © 2003 The Linnean Society of London, Biological Journal of the Linnean Society , 2003, 78 , 129–147.     

A comparative analysis of population structure of a forest tree, Eucalyptus globulus (Myrtaceae), using microsatellite markers and quantitative traits

Eucalyptus globulus (Myrtaceae) is a forest tree native to southeastern Australia, but is grown globally for pulpwood and timber. Eight microsatellite loci were used to determine the degree of selectively neutral differentiation between native populations of the geographic races of E. globulus that are used in a national breeding programme. Spatial differentiation was detected among 340 samples from across the species range (F _ST=0.09±0.02). Analysis of molecular variance showed that there was significant variation between the races, and an unweighted pair group method with arithmetic mean analysis of Nei’s genetic distance between races showed that geographically proximal races tended to be more closely related than geographically distant races. This contrasted markedly with analyses based on quantitative genetic data, where some races appeared to be highly divergent from their geographically closest neighbours. Comparison of racial differentiation based on quantitative (Q _ST) and molecular (F _ST) data suggested that at least five of the quantitative traits used for defining races of E. globulus have been influenced by natural selection, resulting in cases of both phenotypic divergence of parapatric races and phenotypic convergence of allopatric races. We conclude that selectively neutral molecular markers are more useful than quantitative genetic data for identifying the evolutionary affinities and lineages within E. globulus. However, both sources of information should be used in defining evolutionarily important units for conservation. The population structure observed in E. globulus has important consequences for future association studies and may also affect breeding strategies if significant genome co-adaptation has occurred.     

How genetically variable are Neottia ovata (Orchidaceae) populations in northeast Poland?

Using 32 allozyme loci, we examined the genetic diversity of ten populations of Neottia ovata differing in size and located in different regions of northeast Poland. Within‐population genetic variation was low (P_pl = 9.4–31.2%, A = 1.09–1.34 and H_E = 0.044–0.128) relative to taxa with similar life histories. In the majority of N. ovata populations (excluding those in Biebrza), we noted high genotypic diversity (0.49?1.00) and high proportions of unique genotypes (86–100%). Neottia ovata shows a strong pattern of genetic differentiation among populations in northeast Poland, reflected in the high overall F_ST value (0.409). There was a weak, but significant, relationship between genetic and geographical distance (r = 0.09, P < 0.05). We discuss the breeding system, small population sizes and population subdivision as the most important factors affecting the genetic diversity of this species. © 2012 The Linnean Society of London, Botanical Journal of the Linnean Society, 2012, 170 , 40–49     

Identification of evolutionarily significant units in the Cuban endemic damselfly <Emphasis Type="Italic">Hypolestes trinitatis</Emphasis> (Odonata: Hypolestidae)

Species classification may not reflect the underlying/cryptic genetic diversity which should otherwise be conserved as it represents the potential of populations to evolve and adapt. The identification of evolutionarily significant units (ESUs) allows cryptic genetic diversity to be taken into account when designating conservation priorities. Here, we used mitochondrial and nuclear DNA sequences integrated with ecological niche models (ENM) to identify ESUs in Hypolestes trinitatis, a threatened Cuban endemic damselfly species. We found that this species comprises two distinct genetic groups in Central and Eastern Cuba respectively, which are also geographically isolated, as shown by ENM. Therefore, we propose these groups to be considered as different ESUs. According with their extent of occurrence, number of locations and inferred decline of habitat extent and quality, Central and Eastern ESUs qualify as Endangered [EN B1b(iii)] and Vulnerable [VU B1b(iii)], respectively.     

Population structure and cryptic speciation in bonnethead sharks Sphyrna tiburo in the south‐eastern U.S.A. and Caribbean

Population structure and lineage diversification within a small, non‐dispersive hammerhead shark species, the bonnethead shark Sphyrna tiburo, was assessed. Sphyrna tiburo is currently described as one continuously distributed species along the Atlantic continental margins of North, Central and South America, but recent genetic analysis of an insular population (Trinidad) suggests the possibility of cryptic speciation. To address this issue S. tiburo were sampled at six sites along c. 6200 km of continuous, continental coastline and from one island location (Grand Bahama) across a discontinuity in their distribution (the Straits of Florida), in order to test if they constitute a single lineage over this distribution. A total of 1030 bp of the mitochondrial control region (CR) was obtained for 239 S. tiburo, revealing 73 distinct haplotypes, high nucleotide diversity (0·01089) and a pair of highly divergent lineages estimated to have separated 3·61–5·62 million years ago. Mitochondrial cytochrome oxidase I and nuclear internal transcribed spacer loci show the same pattern. Divergence is similar within S. tiburo to that observed between established elasmobranch sister species, providing further evidence of cryptic speciation. A global AMOVA based on CR confirms that genetic diversity is primarily partitioned among populations (Φ_ST = 0·828, P < 0·001) because the divergent lineages are almost perfectly segregated between Belize and North America–The Bahamas. An AMOVA consisting solely of the North American and Bahamian samples is also significantly different from zero (Φ_ST = 0·088, P < 0·001) and pairwise F_ST is significantly different between all sites. These findings suggest that S. tiburo comprises a species complex and supports previous research indicating fine population structure, which has implications for fisheries management and biodiversity conservation.     

Molecular taxonomy and population structure of the rough‐toothed dolphin Steno bredanensis (Cetartiodactyla: Delphinidae)

Several delphinid species have shown genetic population structure, both between and within ocean basins. We investigated genetic differentiation in the rough‐toothed dolphin, Steno bredanensis, using mitochondrial control region sequences from several localities worldwide (N = 112). Preliminary analyses indicated high levels of genetic differentiation between the Atlantic and Pacific/Indian Oceans, which were further investigated using complete cytochrome b sequences and mitogenomes. Phylogenetic analyses were inconclusive about the existence of cryptic speciation in the genus Steno. Notwithstanding this result, analysis of molecular variance and Φ‐statistics analyses revealed strong population differentiation not only between the Atlantic and Pacific, but also within the Atlantic, where three populations were detected: Caribbean, southeastern Brazil, and southern Brazil. We propose that these populations be considered management units for conservation purposes. Our results provide the first perspective on the worldwide genetic differentiation of S. bredanensis. © 2015 The Linnean Society of London     

Microsatellite variation in Bavarian populations of European grayling (Thymallus thymallus): Implications for conservation

European grayling populations in Bavaria have shown steady declines during the last 10–20years. In order to provide guidelines for conservation strategies and future management programs, we investigated the genetic structure of 15 grayling populations originating from three major Central European drainages (the Danube, the Elbe and the Rhine/Main) using 20 microsatellite loci. Genetic divergence between the three drainage systems was substantial as illustrated by highly significant heterogeneity of genotype frequencies, high number of drainage-specific private alleles, high between-drainage F _ST values, high assignment success of individuals to their drainage of origin and the high bootstrap support for the genetic distance based drainage-specific population clusters. In agreement with earlier studies, microsatellites revealed relatively low levels of intrapopulational genetic diversity in comparison to the overall level of variation across populations. Maximum likelihood methods using the coalescent approach revealed that the proportion of common ancestors was generally high in native populations and that the estimates of N _e were correlated with the genetic diversity parameters in all drainages. The number of effective immigrants per generation (N _e m) was less than one for all pairwise comparisons of populations within the drainages, indicating restricted interpopulational gene flow. Based on these findings we recommend a drainage and sub-drainage specific conservation of grayling populations in order to preserve their overall genetic diversity and integrity. For large-scale stocking actions to supplement declining or to restore extinct populations, creation of separate broodstocks for major conservation units (ESUs and MUs) is warranted. This revised version was published online in July 2006 with corrections to the Cover Date.     

Genetic structure of Korean populations of bumblebees Bombus ignitus (Hymenoptera: Apidae) as revealed by microsatellite markers

Bumblebee Bombus ignitus, which is indigenous to Korea, Japan, and China, has been recognized as a valuable pollinator for both crops and wild plants. Bombus ignitus has now become commercially important as a pollinator because of its use in the agricultural industry, particularly for greenhouse pollination. For long‐term management and effective conservation of B. ignitus, an understanding of the genetic structure of its naturally occurring populations is practically important. In this study, the genetic structure among the five populations of B. ignitus in South Korea was assessed using nine microsatellite loci. These markers showed high levels of genetic variability, with the number of alleles ranging from 6 to 22 (mean = 13.4) and the frequency of the most common allele ranging from 0.11 to 0.66. Only the Sabuk (SB) population showed a genetic signature of a recent bottleneck, which was further supported by the lowest level of allelic richness (AR) (mean AR = 3.944). Genetic differentiation was highly significant among all pairs of populations (P < 0.001) across the nine microsatellite markers, suggesting a lack of gene flow among those populations. Interestingly, F_ST (and R_ST and D_est) values were always greater for the Taebaek population than for the four remaining populations. The phylogenetic analysis showed evidence supporting our hypothesis that the Taebaek population is genetically more divergent than the other populations. Overall, our results suggest that the Korean populations of B. ignitus might have undergone geographic isolation and have become highly separated spatially from one another, thereby having a limited range of migration among their local habitats.     

Conservation genetics of Blanding’s turtle and its application in the identification of evolutionarily significant units

Blanding’s turtle is a North American freshwater turtle whose main range occurs south of the Great Lakes; disjunct populations occur east of the Appalachian Mountains from New York to Nova Scotia. The species is listed as threatened or endangered in most of its range. We employed five variable microsatellites to examine samples of 300 individuals in 12 populations. Estimates of F _ST based on pairwise comparisons of populations ranged from 0.000 to 0.465. Phylogenetic analysis of these F _ST values reveals that the Appalachian Mountains and the Hudson River appear to present major barriers to gene flow in Blanding’s turtle. The extent of fine-scale genetic structure previously reported in the Nova Scotian populations was not found in other parts of the species’ range. We recommend that populations separated by the Appalachian Mountains as well as the highly disjunct Nova Scotian populations of Blanding’s turtle be recognized as evolutionarily significant units.     

Cryptic speciation and phylogeographic relationships in the elephant ear sponge Ianthella basta (Porifera,Ianthellidae) from northern Australia

Morphological delineation of sponge species is hindered by the narrow range of fixed diagnostic characters and our limited knowledge of how much phenotypic plasticity the sponge body plan assumes in response to environmental conditions. Here, we make use of the partial mitochondrial cytochrome c oxidase subunit I (COI) gene and the second internal transcribed spacer (ITS2) region to assess the taxonomic validity of colour morphotypes observed in the elephant ear sponge Ianthella basta (Pallas, 1776) across its distribution range in northern Australia, and explore levels and patterns of genetic diversity among populations of the species collected from both sides of the Torres Strait. Molecular phylogenetic analyses revealed congruent topologies consistent with three evolutionarily significant units (ESUs) that were independent of the morphology of the sponge. ESU I includes previously morphologically and genetically delineated western Pacific specimens of I. basta (Guam), and probably corresponds to the type specimen originally described from Indonesia. ESU I occurs in almost all sampling sites across northern Australia, suggesting considerable levels of connectivity among reefs throughout the Torres Strait. ESUs II and III are each exclusively associated with a geographic region of high sponge species richness separated by Torres Strait, and probably represent the result of historical population fragmentation. © 2012 The Linnean Society of London, Zoological Journal of the Linnean Society, 2012, 166 , 225–235.     

Mitochondrial phylogeography and population history of finless porpoises in Sino‐Japanese waters

Phylogeographical patterns of intraspecific variation can provide insights into the population‐level processes responsible for speciation and yield information useful for conservation purposes. In the present study, three hundred and forty‐five base pairs of the mitochondrial DNA control region were sequenced to analyse the genetic diversity, population structure and history, and phylogeography of finless porpoises (Neophocaena phocaenoides) in Chinese and Japanese waters of the North Pacific. Nucleotide and haplotype diversities were 0.44% and 0.79 ± 0.01 for Chinese and Japanese waters, respectively, but varied significantly among populations. Analysis of molecular variance showed a high level of genetic structure between populations (Φ_ST = 0.61, P < 0.001; F_ST = 0.52, P < 0.001). Eleven of 18 haplotypes were restricted to a single population, common haplotypes were found in two to four populations, but no haplotype was found throughout Sino‐Japanese waters, suggesting multiple colonization events followed by limited gene flow. The inferred age of demographic expansion was from the end stage of the last ice age to the Holocene. No obvious phylogeographical pattern was revealed, including between saline and fresh water populations. A low level of genetic diversity for each population and high among‐population differentiation in haplotype frequency were revealed, which suggest a role for random genetic drift, recent demographic bottlenecks, and reduced or limited gene flow in these populations. Some conservation considerations, with special reference to the unique Yangtze population, are discussed. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 193–204.     

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.