Maintenance of an endemic lineage of brown trout (Salmo trutta) within the Duero river basin

The Iberian Peninsula contains diverse populations of freshwater fish, with major river basins comprising differentiated biogeographic units. The Duero River flows through the North‐Western Iberian Peninsula and is one of the most important rivers within the Iberian glacial refuge. Brown trout (Salmo trutta) populate this whole basin, and studies using both allozyme and microsatellite loci have detected a geographically sorted distribution of genetic variation in this species. In this work, sequences of the mitochondrial control region obtained from 299 brown trout from the Duero River were compared with other Iberian and European datasets. Two differentiated haplotype groups were detected inside the Duero River basin. One of them was related to the Atlantic (AT) lineage that is present in Northern European populations, whereas the other comprised an unique group that was restricted to the inner region of the basin. The amount of divergence of this Duero group from the other brown trout populations studied is consistent with a new trout lineage (Duero, DU) that is endemic to this river basin and that diverged from other Atlantic populations during the Pleistocene. The distribution of the DU and AT quaternary lineages in the Duero River was consistent with the ichthyological pattern described in the basin that originated during the Miocene–Pliocene. Evidence of selective processes that favour the haplotypes of the DU lineage may explain this discrepancy.     

Low genetic differentiation among reed warbler Acrocephalus scirpaceus populations across Europe

Migratory birds generally have higher dispersal propensity than resident species and are thus expected to show less genetic differentiation. On the other hand, specific migration patterns may promote genetic structure, such as in situations where migratory divides impede random mixing of individuals. Here we investigated population genetic structure and gene flow patterns in a polytypic passerine, the reed warbler Acrocephalus scirpaceus which shows a migratory divide in central Europe. Using ten polymorphic microsatellite loci and extensive sampling we found low but significant overall genetic differentiation (F_ST=0.013, G’_ST=0.078, D=0.063). Hierarchical F‐statistics and barrier analyses showed low but significant genetic differentiation of Iberian populations, and also slight genetic differences across the migratory divide and between subspecies (A. s. scirpaceus and A. s. fuscus). Three individual‐based Bayesian methods, however, inferred a single genetic unit. Our study thus found low levels of genetic differentiation among reed warbler populations but this genetic differentiation was not pronounced enough to detect a clear population structure using the microsatellite data and no prior information on geographic location of the sampled individuals. This result indicates high levels of gene flow and suggests a possibly recent divergence of European populations after a rapid range expansion. Further studies are necessary to assess divergence times and to reveal the evolutionary history of the reed warbler populations.     

Unravelling plant diversification: Intraspecific genetic differentiation in hybridizing Anacyclus species in the western Mediterranean Basin

Premise

The interfertile species Anacyclus clavatus, A. homogamos, and A. valentinus represent a plant complex coexisting in large anthropic areas of the western Mediterranean Basin with phenotypically mixed populations exhibiting a great floral variation. The goal of this study was to estimate the genetic identity of each species, to infer the role of hybridization in the observed phenotypic diversity, and to explore the effect of climate on the geographic distribution of species and genetic clusters.

Methods

We used eight nuclear microsatellites to genotype 585 individuals from 31 populations of three Anacyclus species for population genetic analyses by using clustering algorithms based on Bayesian models and ordination methods. In addition, we used ecological niche models and niche overlap analyses for both the species and genetic clusters. We used an expanded data set, including 721 individuals from 129 populations for ecological niche models of the genetic clusters.

Results

We found a clear correspondence between species and genetic clusters, except for A. clavatus that included up to three genetic clusters. We detected individuals with admixed genetic ancestry in A. clavatus and in mixed populations. Ecological niche models predicted similar distributions for species and genetic clusters. For the two specific genetic clusters of A. clavatus, ecological niche models predicted remarkably different areas.

Conclusions

Gene flow between Anacyclus species likely explains phenotypic diversity in contact areas. In addition, we suggest that introgression could be involved in the origin of one of the two A. clavatus genetic clusters, which also showed ecological differentiation.     

Species identification of Chilean Artemia populations based on mitochondrial DNA RFLP analysis

Aim Two species of the brine shrimp, namely Artemia franciscana Kellogg and A. persimilis Piccinelli and Prosdocimi, inhabit Chile. Most studies so far have shown that A. franciscana is the most widely distributed species in Chile, with A. persimilis present only in Chilean Patagonia. In general, there is good agreement between morphological and genetic comparisons of Chilean populations with respect to species discrimination. However, a number of results indicate an overlap with some populations tending to diverge from A. franciscana and/or resembling A. persimilis. Prior to the mid 90's the use of DNA markers in Artemia was rather limited, despite their successful application in numerous other species. In this study, we investigate whether the conclusions drawn from traditional comparative tools are congruent with the pattern of genetic divergence depicted by DNA analysis at the mitochondrial level. Location Eight sites in Chile and two reference samples of A. franciscana and A. persimilis from San Francisco Bay (USA) and Buenos Aires (Argentina), respectively. Methods Restriction fragment length polymorphism (RFLP) analysis of a 535 bp segment of the mitochondrial 16S rRNA gene with nine restriction enzymes in 240 individuals. Results No haplotype was shared between the two species. Five restriction enzymes produced species‐specific patterns, enabling the unambiguous assignment of populations to species. Very high (100%) bootstrap values supported the clustering of haplotypes in two groups corresponding to the two species. The two species were clearly differentiated with average sequence divergence of 12.3%. High genetic differentiation was also found among con‐specific populations of A. franciscana with an F_ST estimate of 91%. Main conclusions The mitochondrial DNA (mtDNA) results of this study show a broadly similar pattern to those of previous allozyme and nuclear DNA analyses, with the two New World species appearing as highly divergent. The presence of A. persimilis in southern Chile (Chilean Patagonia) was confirmed. Hence, a species previously regarded as geographically restricted mainly to Argentina, appears to have expanded its range. Populations of A. franciscana appear highly structured with a level of inter‐population genetic differentiation much higher for mtDNA than previously reported with allozymes. Clustering of these populations does not follow a clear geographic pattern. The identification of population‐specific genetic markers for A. persimilis and A. franciscana will help to tackle further aspects of the speciation patterns of these species.     

Differences in allozyme patterns between Diaphanosoma brachyurum and Diaphanosoma mongolianum,as revealed in Central European populations

Recent taxonomic revisions in the genus Diaphanosoma have made the uncertainty in species discrimination with morphological characters obvious. Therefore species characterization on a genetic basis seems to be required.Here we examined the genetic structure of two Diaphanosoma mongolianum populations and three Diaphanosoma brachyurum populations in Central Europe by allozyme electrophoresis. A genetic differentiation between both species was evident. 5 out of 9 tested loci carried diagnostic alleles. Both species differed in their habitat choice: D. mongolianum was adapted to higher trophic levels than D. brachyurum. Co-occurrence was observed in a eutrophic lake. Populations sampled from their preferred habitats showed high genotype diversities.Most genotype frequency distributions were in agreement with Hardy-Weinberg expectations and temporal (between years) changes in genotypic compositions were not common. This indicated annually occurring bouts of sexual reproduction, similar to intermittent Daphnia populations. Genetic differentiation among lakes suggested restricted gene flow.     

Phylogeography of island canary (Serinus canaria) populations

Island canaries (Serinus canaria) are characterised as a species living exclusively on North Atlantic islands, mainly on the Azores, Madeira and Canary Islands. Although they are very common in their habitats, their behaviour and breeding system has only recently been studied systematically. To advance the understanding of their ecology and to see if the rather isolated archipelagos are already promoting a genetic differentiation, we investigated their phylogeographic relationship as revealed by mtDNA sequences of the cytochrome b gene and investigated whether this measure corresponds to morphological characteristics within the islands. Genetic distances were very low throughout the distribution range of the species. Although the variation of genetic distances within the population of Pico (Azores) was larger than that on Madeira and Canary Islands, the genetic distances between island populations were very low throughout which prevented a clear phylogeographic differentiation. Moreover, morphological measurements did not reveal a consistent pattern to reliably separate the populations, although the measures of beak length and body weight revealed a clear island-specific differentiation. These data lead to the assumption that the colonisation of the Atlantic islands by the canaries occurred very recently, while there is no persisting gene flow between the populations.     

Differentiation with drift: a spatio-temporal genetic analysis of Galápagos mockingbird populations (Mimus spp.)

Small and isolated island populations provide ideal systems to study the effects of limited population size, genetic drift and gene flow on genetic diversity. We assessed genetic diversity within and differentiation among 19 mockingbird populations on 15 Galápagos islands, covering all four endemic species, using 16 microsatellite loci. We tested for signs of drift and gene flow, and used historic specimens to assess genetic change over the last century and to estimate effective population sizes. Within-population genetic diversity and effective population sizes varied substantially among island populations and correlated strongly with island size, suggesting that island size serves as a good predictor for effective population size. Genetic differentiation among populations was pronounced and increased with geographical distance. A century of genetic drift did not change genetic diversity on an archipelago-wide scale, but genetic drift led to loss of genetic diversity in small populations, especially in one of the two remaining populations of the endangered Floreana mockingbird. Unlike in other Galápagos bird species such as the Darwin''s finches, gene flow among mockingbird populations was low. The clear pattern of genetically distinct populations reflects the effects of genetic drift and suggests that Galápagos mockingbirds are evolving in relative isolation.     

Gene flow and spatio-temporal genetic variation among sympatric populations of Tetranychus kanzawai (Acari: Tetranychidae) occurring on different host plants, as estimated by microsatellite gene diversity

We investigated spatio-temporal genetic variation in allele frequency and estimated gene flow among sympatric populations of Tetranychus kanzawai on different host plants by the use of microsatellite markers. In the analysis of spatial genetic variation, no isolation by distance was detected among the populations. Gene flow between populations on Hydrangea macrophylla and those on other host plants was relatively restricted, whereas the populations on Akebia quinata and Clerodendrum trichotomum were almost panmictic. Our study on temporal genetic variation showed (1) that population differentiation was slightly reduced during the period from April to May owing to frequent gene flow among populations; and (2) that population differentiation was greatly enhanced from May to October because of bottleneck effects. Genetic differentiation among T. kanzawai populations was caused by the effect of host plants rather than by the effect of geographic distance among populations, suggesting possibility of sympatric host race formation in this species.This revised version was published online in May 2005 with a corrected cover date.     

Geographically structured genetic variation in the Medicago lupulina–Ensifer mutualism

Gene flow between genetically differentiated populations can maintain variation in species interactions, especially when population structure is congruent between interacting species. However, large‐scale empirical comparisons of the population structure of interacting species are rare, particularly in positive interspecific interactions (mutualisms). One agriculturally and ecologically important mutualism is the partnership between legume plants and rhizobia. Through characterizing and comparing the population genomic structure of the legume Medicago lupulina and two rhizobial species (Ensifer medicae and E. meliloti), we explored the spatial scale of population differentiation between interacting partners in their introduced range in North America. We found high proportions of E. meliloti in southeastern populations and high proportions of E. medicae in northwestern populations. Medicago lupulina and the Ensifer genus showed similar patterns of spatial genetic structure (isolation by distance). However, we detected no evidence of isolation by distance or population structure within either species of bacteria. Genome‐wide nucleotide diversity within each of the two Ensifer species was low, suggesting limited introduction of strains, founder events, or severe bottlenecks. Our results suggest that there is potential for geographically structured coevolution between M. lupulina and the Ensifer genus, but not between M. lupulina and either Ensifer species.     

Are populations of mayflies living in adjacent fish and fishless streams genetically differentiated?

1. Conspecific populations living in habitats with different risks of predation often show phenotypic variation in defensive traits. Traits of two species of mayflies (Baetidae: Baetis bicaudatus and Baetis sp. nov.) differ between populations living in fish and fishless streams in a high altitude drainage basin in western Colorado, U.S.A. We tested for genetic differentiation between mayfly populations in these two habitat types, assuming that lack of genetic differentiation would be consistent with the hypothesis that those traits are phenotypically plastic. 2. Previous work has shown that larvae of both species behave differently and undergo different developmental pathways in adjacent fish and fishless streams. These phenotypic differences in behaviour and development have been induced experimentally, suggesting that populations from fishless streams have the genetic capability to respond to fish. 3. During summer 2001 we collected Baetis larvae from several fish and fishless streams, and from fish and fishless sections of the same streams. We used allozymes and a fragment of the cytochrome oxidase subunit 1 mitochondrial gene to examine genetic variation of Baetis individuals within and among streams. 4. Results showed that genetic variation exists among populations of the same species of Baetis from different streams, but none of that variation was associated with the presence or absence of fish. These data confirm that populations of Baetis living in fish and fishless streams are not genetically distinct, and are consistent with the hypothesis that traits associated with environments of different risk are phenotypically plastic.     

Genetic differentiation among populations ofArabis serrata (Brassicaceae) along its geographic distribution

Levels and distribution of genetic variation were investigated in ten populations of the perennialArabis serrata distributed along a latitudinal gradient throughout Japan. Populations of this endemic species occupy predominantly three types of habitats: limestone and derived soils, volcanic and disturbed sites. Previous studies showed that plants ofA. serrata are differentiated in morphological and ecological traits under both natural conditions and common garden experiments suggesting genetic differentiation among populations. To assess the degree of genetic differentiation under different habitats ofA. serrata populations, we analyzed the isozyme genetic structure. Ten populations, located in mountains of central and northern Honshu and Hokkaido, were analyzed by starch gel enzyme electrophoresis. Fourteen loci of eight enzymes were resolved and six loci were monomorphic for all the populations. Populations sampled maintain low levels of genetic variation (P=0.16;H=0.05;A=1.16) compared to that maintained by other outcrossing seed plants. In some cases, few or no heterozygous individuals were detected, and consequently, mean observed heterozygosity was zero or near zero. Six (29%) of the fixation indices,F, estimated deviated significantly from Hardy Weinberg genotypic expectations indicating a deficiency of heterozygotes in most of the cases. The mean genetic identity (Nei'sI) between population pairs was 0.852 and indicates a moderate level of genetic differentiation among populations.Arabis serrata has most of its genetic variation distributed among rather than within populations. The among-population component of the total genetic diversity (G _ST mean value) was 0.416, indicating genetic differentiation between populations. There groups of populations were recognized in an unrooted tree generated by the Neighbor-Joining method. These results suggest groups of populations differentiated regionally. Estimates of interpopulational gene flow (Nm) were very variable (range 0.049–3.718) with a meanNm=1.203 for all populations.     

Morphological variation in a secondary contact between divergent lineages of brown trout (Salmo trutta) from the Iberian Peninsula

The aim of this study was to analyze the morphological variation of brown trout (Salmo trutta) in the Duero basin, an Atlantic river basin in the Iberian Peninsula, where a spatial segregation of two divergent lineages was previously reported, based on isozyme, microsatellite and mtDNA data. In these studies, two divergent pure regions (Pisuerga and Lower-course) and several hybrid populations between them were identified. Morphological variation was evaluated in 11 populations representative of the genetic differentiation previously observed in the Duero basin, using multivariate analysis on 12 morphometric and 4 meristic traits. A large differentiation between populations was observed (interpopulation component of variance: 41.8%), similar to that previously detected with allozymes and microsatellites. Morphometric differentiation was also reflected by the high classification success of pure and hybrid individuals to their respective populations, using multivariate discriminant functions (94.1% and 79.0%, respectively). All multivariate and clustering analyses performed demonstrated a strong differentiation between the pure regions. The hybrid populations, though showing large differentiation among them, evidenced an intermediate position between the pure samples. Head and body shape traits were the most discriminant among the morphometric characters, while pectoral rays and gillrakers were the most discriminant among the meristic traits. These results confirmed the high divergence of the brown trout from the Duero basin and suggest some traits on which selection could be acting to explain the spatial segregation observed.     

High Ag-NOR-site variation associated to a secondary contact in brown trout from the Iberian Peninsula

The analysis of nucleolar organizer regions (NORs) using silver (Ag-) staining and in situ hybridization (ISH) in brown trout (Salmo trutta) from various river basins in the Iberian Peninsula revealed high variation in the number and location of NORs. A total of 17 different Ag-NOR sites were revealed in 10 different chromosome pairs. Three different Ag-NOR patterns clustered by river basins and strongly associated to the internal transcribed spacer 1 (ITS1) variation were detected. The main variability in NOR-sites was found in a secondary contact between two divergent lineages of brown trouts at Duero basin. Our results confirmed the abrupt break in the spatial distribution of genetic variation of brown trout populations previously reported at Duero basin. We hypothesize that NOR-site variation might be a consequence of hybridization between divergent lineages of brown trouts and that NORs could play a major role in the maintenance of a hybrid zone in Duero basin via post-zygotic isolation mechanisms.     

Genetic differentiation among morphological variants of Acacia saligna (Mimosaceae)

This study investigated the pattern of genetic variation in Acacia saligna (Labill.) H.L. Wendl. to facilitate its development as a crop species for dryland salinity management. A. saligna is a morphologically variable species and four main variants are recognized. The genetic structure within A. saligna was investigated in populations across the geographic range of the species using nuclear restriction fragment length polymorphism loci. The analysis identified considerable genetic variation within A. saligna that was genetically structured into three groups. Two of the three groups corresponded to variants recognized in the field study; the third group encompassed the other two variants which, though morphologically different, were not genetically differentiated. The level of genetic differentiation between the groups suggests they may represent different taxa and a taxonomic revision of the species may be required. Identification of different taxa within A. saligna will have implications for the utilization and domestication of the species, as the taxa will need to be evaluated separately to determine their suitability for agroforestry. The high genetic variation between and within groups suggests there is a large genetic base available for breeding improved cultivars.     

Genetic Structure and Gene Flow in Elymus glaucus (blue wildrye): Implications for Native Grassland Restoration

Interest in using native grass species for restoration is increasing, yet little is known about the ecology and genetics of native grass populations or the spatial scales over which seed can be transferred and successfully grown. The purpose of this study was to investigate the genetic structure within and among populations of Elymus glaucus in order to make some preliminary recommendations for the transfer and use of this species in revegetation and restoration projects. Twenty populations from California, Oregon, and Washington were analyzed for allozyme genotype at 20 loci, and patterns of variation within and among populations were determined. Allozyme variation at the species level was high, with 80% of the loci polymorphic and an average expected heterozygosity (an index of genetic diversity) of 0.194. All but two of the populations showed some level of polymorphism. A high degree of population differentiation was found, with 54.9% of the variation at allozyme loci partitioned among populations (F_st= 0.549). A lesser degree of genetic differentiation among closely spaced subpopulations within one of the populations was also demonstrated (F_st= 0.124). Self-pollination and the patchy natural distribution of the species both likely contribute to the low level of gene flow (N_m= 0.205) that was estimated. Zones developed for the transfer of seed of commercial conifer species may be inappropriate for transfer of E. glaucus germplasm because conifer species are characterized by high levels of gene flow. Limited gene flow in E. glaucus can facilitate the divergence of populations over relatively small spatial scales. This genetic differentiation can be due to random genetic drift, localized selective pressures, or both. In order to minimize the chances of planting poorly adapted germplasm, seed of E. glaucus may need to be collected in close proximity to the proposed restoration site.     

Genetic diversity and structuring across the range of a widely distributed ladybird: focus on rear‐edge populations phenotypically divergent

Population genetics and phenotypic structures are often predicted to vary along the geographic range of a species. This phenomenon would be accentuated for species with large range areas, with discontinuities and marginal populations. We herein compare the genetic patterns of central populations of Coccinella septempunctata L. with those of two phenotypically differentiated populations considered as rear‐edge populations and subspecies based on phenotype (Algeria and Japan). According to the central‐marginal model and expected characteristics of rear‐edge populations, we hypothesize that these rear‐edge populations have (1) a reduced genetic diversity, resulting from their relative isolation over long periods of time, (2) a higher population genetic differentiation, explained by low contemporary gene flow levels, and (3) a relationship between genetic diversity characteristics and phenotypes, due to historical isolation and/or local adaptation. Based on genotyping of 28 populations for 18 microsatellite markers, several levels of regional genetic diversity and differentiation are observed between and within populations, according to their localization: low within‐population genetic diversity and higher genetic differentiation of rear‐edge populations. The genetic structuring clearly dissociates the Algerian and Eastern Asia populations from the others. Geographical patterns of genetic diversity and differentiation support the hypothesis of the central‐marginal model. The pattern observed is in agreement with the phenotypic structure across species range. A clear genetic break between populations of Algeria, the Eastern Asia, and the remaining populations is a dominant feature of the data. Differential local adaptations, absence of gene flow between marginal and central populations, and/or incapacity to mate after colonization, have contributed to their distinct genotypic and phenotypic characteristics.     

On the evolution of heavy-metal tolerant populations in Armeria maritima: evidence from allozyme variation and reproductive barriers

The process of ecological differentiation leading to the evolution of heavy-metal tolerant populations in Armeria maritima was studied by comparing population genetic structure and pattern of gene flow between populations growing on heavy-metal contaminated against non-contaminated sites using allozyme markers. In addition the evolution of reproductive isolation among populations was studied by measuring pollen fertility in interpopulational hybrids. The allozyme data suggested that in A. maritima multiple independent evolutionary origins of heavy-metal tolerant populations have occurred in the absence of strong genetic bottlenecks. The pattern of gene flow among populations was consistent with the model of isolation by distance with considerable gene flow between neighbor populations, and no reduction of gene flow between tolerant and non-tolerant populations. Hence it appears that substantial gene flow has not hampered genetic differentiation, probably because of the high selection pressure for heavy-metal tolerance. The pattern of reproductive isolation among populations suggests that evolution of heavy-metal tolerant populations has not triggered the development of reproductive barriers against non-tolerant populations. However, partial reproductive isolation has occurred under geographic separation.     

Population structure affected by excess gene flow in self-pollinating Elymus nutans and E. burchan-buddae (Triticeae: Poaceae)

Seed-mediated gene flow can considerably affect population genetic structure of strictly self-pollinating species, but little is known on the extent and nature of such gene flow among pastoral plant populations. Molecular fingerprints provide a powerful tool to address the relevant issues. Genetic structure of 22 populations of two self-pollinating pasture species, Elymus nutans and E. burchan-buddae, collected from various altitudes of the Qinghai–Tibetan Plateau was studied using fluorescence-based amplified fragment length polymorphism technique. Analysis of molecular variance revealed 42.97% and 37.63% among-population variation for the two Elymus species, respectively, indicating that the majority of the total variation presented within populations. This result contradicts the common genetic variation pattern for a selfing plant species: lower genetic variation within populations. Further analysis suggested higher level of gene flow among populations within the same region than among different regions across the sampled area for the two Elymus species. STRUCTURE analyses of the Elymus populations indicated an evident admixture genetic structure, particularly among neighboring populations from the same region, supporting the hypothesis of considerable seed dispersal among populations. The excess within-region gene flow of E. nutans and E. burchan-buddae might be caused by grazing animals that promote seed dispersal when moved around the pastoral lands during foraging. The among-population gene flow promulgated by grazing animals may promote the maintenance of genetic diversity in the pasture species, particularly in small and fragmented populations within a given region.     

Do different rates of gene flow underlie variation in phenotypic and phenological clines in a montane grasshopper community?

Species responses to environmental change are likely to depend on existing genetic and phenotypic variation, as well as evolutionary potential. A key challenge is to determine whether gene flow might facilitate or impede genomic divergence among populations responding to environmental change, and if emergent phenotypic variation is dependent on gene flow rates. A general expectation is that patterns of genetic differentiation in a set of codistributed species reflect differences in dispersal ability. In less dispersive species, we predict greater genetic divergence and reduced gene flow. This could lead to covariation in life‐history traits due to local adaptation, although plasticity or drift could mirror these patterns. We compare genome‐wide patterns of genetic structure in four phenotypically variable grasshopper species along a steep elevation gradient near Boulder, Colorado, and test the hypothesis that genomic differentiation is greater in short‐winged grasshopper species, and statistically associated with variation in growth, reproductive, and physiological traits along this gradient. In addition, we estimate rates of gene flow under competing demographic models, as well as potential gene flow through surveys of phenological overlap among populations within a species. All species exhibit genetic structure along the elevation gradient and limited gene flow. The most pronounced genetic divergence appears in short‐winged (less dispersive) species, which also exhibit less phenological overlap among populations. A high‐elevation population of the most widespread species, Melanoplus sanguinipes, appears to be a sink population derived from low elevation populations. While dispersal ability has a clear connection to the genetic structure in different species, genetic distance does not predict growth, reproductive, or physiological trait variation in any species, requiring further investigation to clearly link phenotypic divergence to local adaptation.     

Genome‐scale sampling suggests cryptic epigenetic structuring and insular divergence in Canada lynx

Determining the molecular signatures of adaptive differentiation is a fundamental component of evolutionary biology. A key challenge is to identify such signatures in wild organisms, particularly between populations of highly mobile species that undergo substantial gene flow. The Canada lynx (Lynx canadensis) is one species where mainland populations appear largely undifferentiated at traditional genetic markers, despite inhabiting diverse environments and displaying phenotypic variation. Here, we used high‐throughput sequencing to investigate both neutral genetic structure and epigenetic differentiation across the distributional range of Canada lynx. Newfoundland lynx were identified as the most differentiated population at neutral genetic markers, with demographic modelling suggesting that divergence from the mainland occurred at the end of the last glaciation (20–33 KYA). In contrast, epigenetic structure revealed hidden levels of differentiation across the range coincident with environmental determinants including winter conditions, particularly in the peripheral Newfoundland and Alaskan populations. Several biological pathways related to morphology were differentially methylated between populations, suggesting that epigenetic modifications might explain morphological differences seen between geographically peripheral populations. Our results indicate that epigenetic modifications, specifically DNA methylation, are powerful markers to investigate population differentiation in wild and non‐model systems.