Patterns of differentiation among wild rabbit populations Oryctolagus cuniculus L. in arid and semiarid ecosystems of north-eastern Australia

Feral rabbit populations in Australia have generally been managed using localized control procedures. While these procedures may result in local extinctions, persistence of populations will depend on the probability of recolonization. Genetic markers developed using temperature gradient gel electrophoresis (TGGE) combined with heteroduplex analysis (HA) of mitochondrial DNA (mtDNA) were used to characterize the degree of subdivision and extent of gene flow within and among rabbit populations distributed over large distances (up to 1000 km) in southern Queensland (QLD) and north-west New South Wales (NSW), Australia. TGGE analyses revealed significant heterogeneity in mtDNA control region haplotype frequencies. From heterogeneity χ² tests, it was evident that the differentiation observed was largely attributable to five sites which were located in the semiarid eastern region, whereas haplotype frequencies were homogeneous throughout the arid western region. These results suggest that there are independent population systems within the study area. The extent of gene flow among local populations within each system is related to the spatial configuration of acceptable habitat patches and the persistence of the populations is determined by the probability of recolonization following local extinction. These data suggest that to provide better overall control of rabbit populations, different management strategies may be necessary in arid and semiarid ecosystems. In arid south-west QLD and north-west NSW, where extensive gene flow occurs over large distances, rabbit populations should be managed at a regional level. In semiarid eastern QLD, where gene flow is restricted and populations are more isolated, localized control procedures may provide effective short-term relief. These results indicate that in nonequilibrium systems with patchy distribution of individuals, the interpretation of migration rate from estimates of gene flow obtained using existing genetic models must include an understanding of the spatial and temporal scales over which population processes operate.     

Local genetic structure in red grouse (Lagopus lagopus scoticus): evidence from microsatellite DNA markers

Allelic variation at seven hypervariable tri- and tetranucleotide microsatellite loci was used to determine levels of population differentiation between 14 populations of red grouse ( Lagopus lagopus scoticus ) in northeast Scotland, UK. Despite the potential for long-distance dispersal in grouse, and a semicontinuous habitat, significant population divergence was observed (mean R _ST = 0.153; P < 0.01) and an isolation-by-distance effect detected (Mantel test: P < 0.001). Examination of the spatial trend in principal component scores derived from allele frequencies among populations highlighted a barrier to gene flow that was confounding a simple isolation-by-distance effect. This barrier corresponded to an area of unsuitable habitat for grouse associated with a river system that bisected the study area. Mean genetic relatedness was higher for males than for females in all but one of the study populations, suggesting that the territorial behaviour and natal philopatry displayed by cocks have a manifold effect in generating the observed spatial genetic structure. Lower female relatedness values suggest a higher level of female-mediated gene flow, which is sufficient to prevent the loss of genetic variation from within populations and the onset of inbreeding effects. The potential consequences of local subdivision for red grouse populations are discussed.     

A further example of natural selection on phenotypes of the land snail Arianta arbustorum (L.)

The land snail Arianta arbustorum is polymorphic for colour and banding of the shell. This paper reports the results of a survey of natural populations of the species in the valley of the River Dove in Derbyshire and Staffordshire. There is some evidence of climatic selection in non-woodland populations such as has been observed elsewhere in Britain. The woodland sites can be divided into several classes depending upon the tree species present. The morph-frequencies are significantly different between class, but homogeneous within. It is argued that the greenness of the ground cover depends upon the trees growing overhead, and that the morph-frequency is related to this. It is suggested that visually discriminating predators may be partially responsible for this correlation. SUMMARY Arianta arbustorum (L.) is a terrestrial land snail which is polymorphic for colour and banding. Brown shell colour is dominant to yellow, and banded is dominant to unhanded. The two loci are closely linked. Communities of larger helicid snails were sampled along the valley of the River Dove in Derbyshire and Staffordshire. The samples were scored for species composition, and the proportions of the various morphs of A. arbustorum were determined. Outside the woodland areas, A. arbustorum shows slight evidence of micro-climatic selection similar to that found elsewhere in Britain. There is little evidence of differences in the ecological genetics of A. arbustorum between woodland and non-woodland habitats despite fairly heavy predation. It is argued that predation may be lower outside the woods allowing climatic effects to become detectable. The woodland samples are heterogeneous with respect to morph frequency. The heterogeneity can be partitioned between three classes of woodland which are internally homogeneous. It is argued that birch woodland may be lighter in background than hazel, with the elm-ash-beech woods being darker than either. The proportion of unhanded yellow A. arbustorum is highest in the first of the classes, and lowest in the last. It is argued that this result could be due to the actions of visually discriminating predators upon these populations of A. arbustorum.     

EVIDENCE FOR SYMPATRIC GENETIC DIVERGENCE OF ANADROMOUS AND NONANADROMOUS MORPHS OF SOCKEYE SALMON (ONCORHYNCHUS NERKA)

Anadromous and nonanadromous morphs of the Pacific salmon Oncorhynchus nerka spawn in close physical proximity in tributaries to Takla Lake, British Columbia, yet differ in morphology, gill raker number, allozyme allele frequencies, and reproductive traits. Both morphs are semelparous typically maturing at age four, the anadromous morph (sockeye) at fork lengths of 38–65 cm and the nonanadromous morph (kokanee) at 17–22 cm. When reared together, pure and hybrid morphs also exhibited different growth rates and maturity schedules. Collectively, these large differences between the morphs confirm that sockeye and kokanee exist as reproductively isolated populations. Average gene flow (m) was estimated to be 0.1–0.8% between morphs, 1.7–3.7% among tributaries for kokanee, and 0.3–5.6% among tributaries for sockeye. We conclude that divergence has occurred in sympatry and examine potential isolating mechanisms.     

Isolation by resistance across a complex coral reef seascape

A detailed understanding of the genetic structure of populations and an accurate interpretation of processes driving contemporary patterns of gene flow are fundamental to successful spatial conservation management. The field of seascape genetics seeks to incorporate environmental variables and processes into analyses of population genetic data to improve our understanding of forces driving genetic divergence in the marine environment. Information about barriers to gene flow (such as ocean currents) is used to define a resistance surface to predict the spatial genetic structure of populations and explain deviations from the widely applied isolation-by-distance model. The majority of seascape approaches to date have been applied to linear coastal systems or at large spatial scales (more than 250 km), with very few applied to complex systems at regional spatial scales (less than 100 km). Here, we apply a seascape genetics approach to a peripheral population of the broadcast-spawning coral Acropora spicifera across the Houtman Abrolhos Islands, a high-latitude complex coral reef system off the central coast of Western Australia. We coupled population genetic data from a panel of microsatellite DNA markers with a biophysical dispersal model to test whether oceanographic processes could explain patterns of genetic divergence. We identified significant variation in allele frequencies over distances of less than 10 km, with significant differentiation occurring between adjacent sites but not between the most geographically distant ones. Recruitment probabilities between sites based on simulated larval dispersal were projected into a measure of resistance to connectivity that was significantly correlated with patterns of genetic divergence, demonstrating that patterns of spatial genetic structure are a function of restrictions to gene flow imposed by oceanographic currents. This study advances our understanding of the role of larval dispersal on the fine-scale genetic structure of coral populations across a complex island system and applies a methodological framework that can be tailored to suit a variety of marine organisms with a range of life-history characteristics.     

Genetic differentiation between sun and shade habitats in populations of Lindera benzoin L.

Differences in selection patterns among habitats can alter the distribution of genetic diversity even when this is estimated with neutral markers. For plants, light is an essential resource that can influence both abiotic and biotic components of habitat. We examined genetic differentiation between sun and shade habitats in Lindera benzoin L. (Spicebush), a perennial understory shrub. Genetic diversity of 127 plants from sun and shade habitats in two populations of L. benzoin was determined using 12 polymorphic microsatellite markers. We analyzed patterns of genetic diversity using analysis of molecular variance (AMOVA), and we assessed correlation between genetic and geographic distance using Mantel tests. We found (1) low levels of differentiation among populations (F _ST = 0.028), (2) little evidence of genetic structure within populations due to isolation-by-distance, and (3) some evidence of habitat-based genetic differentiation. Specifically, the AMOVA showed a small (0.5%) but significant portion of overall variation could be explained by differences between habitats. The overall low levels of differentiation we saw were likely a result of extensive gene flow in this dioecious, bird-dispersed species.     

GENETIC AND MORPHOLOGICAL DIVERGENCE AMONG MORPHOTYPES OF THE ISOPOD EXCIROLANA ON THE TWO SIDES OF THE ISTHMUS OF PANAMA

Excirolana braziliensis is a dioecious marine isopod that lives in the high intertidal zone of sandy beaches on both sides of Central and South America. It possesses no larval stage and has only limited means of adult dispersal. Indirect estimates of gene flow have indicated that populations from each beach exchange less than one propagule per generation. Multivariate morphometrics have discovered three morphs of this species in Panama, two of them closely related and found on opposite sides of Central America (“C morph” in the Caribbean and “C′ morph” in the eastern Pacific), the third found predominantly in the eastern Pacific (“P morph”). Though the P and C′ morphs are seldom found on the same beach, they have overlapping latitudinal ranges in the eastern Pacific. A related species, Excirolana chamensis, has been described from the Pacific coast of Panama. Each beach contains populations that remain morphologically and genetically stable, but a single drastic change in both isozymes and morphology has been documented. We studied isozymes and multivariate morphology of 10 populations of E. braziliensis and of one population of E. chamensis. Our objective was to assess the degree of genetic and morphological variation, the correlation of divergence on these two levels of integration, the phylogenetic relationships between morphs, and the possible contributions of low vagility, low gene flow, and occasional extinction and recolonization to the genetic structuring of populations. Genetic distance between the P morph, on one hand, and the other two morphotypes of E. braziliensis, on the other, was as high as the distance between E. braziliensis and E. chamensis. Several lines of evidence agree that E. chamensis and the P morph had diverged from other morphs of E. braziliensis before the rise of the Panama Isthmus separated the C and C′ forms, and that the P morph constitutes a different species. A high degree of genetic differentiation also exists between populations of the same morph. On the isozyme level, every population can be differentiated from every other on the basis of at least one diagnostically different locus, regardless of geographical distance or morphological affiliation. Morphological and genetic distances between populations are highly correlated. However, despite the high degree of local variation, evolution of E. braziliensis as a whole has not been particularly rapid; divergence between the C and C′ morphs isolated for 3 million yr by the Isthmus of Panama is not high by the standard of within-morph differentiation or by comparison with other organisms similarly separated. Alleles that are common in one population may be absent from another of the same morph, yet they appear in a different morph in a separate ocean. The high degree of local differentiation, the exclusive occupation of a beach by one genotype with rare arrival of foreign individuals that cannot interbreed freely with the residents, the genetic stability of populations with infrequent complete replacement by another genetic population, and the sharing by morphs of polymorphisms that are not shared by local populations, all suggest a mode of evolution concentrated in rare episodes of extinction and recolonization, possibly coupled with exceptional events of gene flow that help preserve ancestral variability in both oceans.     

Genetic Diversity and Spatial Genetic Structure of the Grassland Perennial Saxifraga granulata along Two River Systems

Due to changes in land use, the natural habitats of an increasing number of plant species have become more and more fragmented. In landscapes that consist of patches of suitable habitat, the frequency and extent of long-distance seed dispersal can be expected to be an important factor determining local genetic diversity and regional population structure of the remaining populations. In plant species that are restricted to riparian habitats, rivers can be expected to have a strong impact on the dynamics and spatial genetic structure of populations as they may enable long-distance seed dispersal and thus maintain gene flow between fragmented populations. In this study, we used polymorphic microsatellite markers to investigate the genetic diversity and the spatial genetic structure of 28 populations of Saxifraga granulata along two rivers in central Belgium. We hypothesized that rivers might be essential for gene flow among increasingly isolated populations of this species. Genetic diversity was high (H_S = 0.68), which to a certain extent can be explained by the octoploid nature of S. granulata in the study area. Populations along the Dijle and Demer rivers were also highly differentiated (G” _ST = 0.269 and 0.164 and D _EST = 0.190 and 0.124, respectively) and showed significant isolation-by-distance, indicating moderate levels of gene flow primarily between populations that are geographically close to each other. Along the river Demer population genetic diversity was higher upstream than downstream, suggesting that seed dispersal via the water was not the primary mode of dispersal. Overall, these results indicate that despite increasing fragmentation populations along both rivers were highly genetically diverse. The high ploidy level and longevity of S. granulata have most likely buffered negative effects of fragmentation on genetic diversity and the spatial genetic structure of populations in riparian grasslands.     

Pollinator visitation in populations of tristylous Eichhornia paniculata in northeastern Brazil

Summary The frequencies of floral morphs in populations of tristylous Eichhornia paniculata often deviate from the theoretical expectation of equality. This variation is associated with the breakdown of tristyly and the evolution of self-fertilization. Differences in morph frequencies could result from selection pressures due to variable levels of insect visitation to populations and contrasting foraging behavior among the floral morphs. We estimated pollinator densities in 16 populations and quantified visitation sequences to morphs in five populations of E. paniculata in northeastern Brazil. Foraging behavior among floral morphs was measured as the frequency of visits to morphs relative to their frequency in the population (preference) and number of flights between inflorescences of the same versus different morphs (constancy). Pollinator density (number/m²/minute) was not correlated with population size, plant density or morph diversity. Pollinator densities varied most among populations of less than 200 plants. Whether pollinators discriminated among the morphs, depended on whether they primarily collected nectar or pollen. In four populations, nectar-feeding bees (Ancyloscelis and Florilegus spp.) and butterflies showed no consistent preference or constancy among the morphs. In contrast, pollen-collecting bees (Trigona sp.) visited a lower proportion of longstyled inflorescences than expected and tended to visit more mid-and short-styled inflorescences in succession, once they were encountered. Pollinator constancy for morphs did not result from differences in inflorescence production or spatial patchiness among the morphs. Although non-random pollinator visitation to morphs in heterostylous populations could potentially affect mating and hence morph frequencies, the observed visitation patterns in this study do not provide evidence that pollinators play a major role in influencing floral morph frequencies.     

Genetic variability in white-tailed deer (Odocoileus virginianus) and its relationship to environmental parameters and herd origin (Cervidae)

Allozyme variation was examined in 1571 white-tailed deer (Odocoileus virginianus) from 29 localities in Tennessee by starch gel electrophoresis. For 11 polymorphic loci, sex-related, age-related and temporal differences were minimal. However, significant spatial hererogeneity was evident in genotypes (contingency table results), allele frequencies (F _ST=0.057) and heterozygosity. Heterozygosity ranged from 16.9% to 26.8% with a mean of 22.9%. The spatial pattern of allele frequencies determined from Rogers' coefficients of genetic similarity indicated associations based on geographic proximity and stocking history. In hierarchial analyses, physiographic regions accounted for more of the total gene diversity than herd origin groups (populations of similar origin) but less than individual populations. For five loci, physiographic regions accounted for more of the gene diversity than populations, suggesting a selection role in the observed genetic variability. Bivariate and canonical correlation analyses revealed significant associations between environmental and genetic variables. Temperature variables and allele frequencies for three loci (alcohol dehydrogenase, alpha-glycerophosphate dehydrogenase, sorbitol dehydrogenase) had the prominent roles in the multivariate association between environmental and genetic variables. Herd origin, gene flow and selection appear to be involved in the gene diversity in deer from Tennesee.     

Parasites and parallel divergence of the number of individual MHC alleles between sympatric three‐spined stickleback Gasterosteus aculeatus morphs in Iceland

Two pairs of sympatric three‐spined stickleback Gasterosteus aculeatus morphs and two single morph populations inhabiting mud and lava or rocky benthic habitats in four Icelandic lakes were screened for parasites and genotyped for MHC class IIB diversity. Parasitic infection differed consistently between G. aculeatus from different benthic habitats. Gasterosteus aculeatus from the lava or rocky habitats were more heavily infected in all lakes. A parallel pattern was also found in individual MHC allelic variation with lava G. aculeatus morphs exhibiting lower levels of variation than the mud morphs. Evidence for selective divergence in MHC allele number is ambiguous but supported by two findings in addition to the parallel pattern observed. MHC allele diversity was not consistent with diversity reported at neutral markers (microsatellites) and in Þingvallavatn the most common number of alleles in each morph was associated with lower infection levels. In the Þingvallavatn lava morph, lower infection levels by the two most common parasites, Schistocephalus solidus and Diplostomum baeri, were associated with different MHC allele numbers.     

Habitat differentiation vs. isolation-by-distance: the genetic population structure of Elymus athericus in European salt marshes

We investigated genetic differentiation among populations of the clonal grass Elymus athericus, a common salt-marsh species occurring along the Wadden Sea coast of Europe. While E. athericus traditionally occurs in the high salt marsh, it recently also invaded lower parts of the marsh. In one of the first analyses of the genetic population structure in salt-marsh species, we were interested in population differentiation through isolation-by-distance, and among strongly divergent habitats (low and high marsh) in this wind- and water-dispersed species. High and low marsh habitats were sampled at six sites throughout the Wadden Sea. Based on reciprocal transplantation experiments conducted earlier revealing lower survival of foreign genotypes we predicted reduced gene flow among habitats. Accordingly, an analysis with polymorphic cross-species microsatellite primers revealed significant genetic differentiation between high and low marsh habitats already on a very small scale (< 100 m), while isolation-by-distance was present only on larger scales (60-443 km). In an analysis of molecular variance we found that 14% of the genetic variance could be explained by the differentiation between habitats, as compared to only 8.9% to geographical (isolation-by-distance) effects among six sites 2.5-443 km distant from each other. This suggests that markedly different selection regimes between these habitats, in particular intraspecific competition and herbivory, result in habitat adaptation and restricted gene flow over distances as small as 80 m. Hence, the genetic population structure of plant species can only be understood when considering geographical and selection-mediated restrictions to gene flow simultaneously.     

Pattern‐oriented modelling of population genetic structure

Although several statistical approaches can be used to describe patterns of genetic variation and infer stochastic differentiation, selective responses, or interruptions of gene flow due to physical or environmental barriers, it is worthwhile to note that similar processes, controlled by several parameters in theoretical models, frequently give rise to similar patterns. Here, we develop a Pattern‐Oriented Modelling (POM) approach that allows us to determine how a complex set of parameters potentially driving empirical genetic differentiation among populations generate alternative scenarios that can be fitted to observed data. We generated 10 000 random combinations of parameters related to population size, gene flow and response to gradients (both driven by dispersal and selection) in a spatially explicit model, and analysed simulated patterns with F_ST statistics and mean correlograms using Moran's I spatial autocorrelation coefficients. These statistics were compared with observed patterns for a tree species endemic to the Brazilian Cerrado. For a best match with observed F_ST (equal to 0.182), the important parameters driving simulated scenario are mainly related to population structure, including low population size with closed populations (low N_m), strong distance decay of gene flow, in addition to a strong effect of the initial variance of allele frequencies. These scenarios present a low autocorrelation of allele frequencies. Best matching of correlograms, on the other hand, appears in simulations with a large population size, high N_m and low population differentiation and F_ST (as well as more gene flow). Thus, targeting the two statistics (correlograms and F_ST) shows that best matches with empirical data with two distinct sets of parameters in the simulations, because observed patterns involve both a relatively high F_ST and significant autocorrelation. This conflict can be resolved by assuming that initial variance in allele frequencies can be interpreted as reflecting deep‐time historical variation and evolutionary dynamics of allele frequencies, creating a relatively high level of population differentiation, whereas current patterns in gene flow creates spatial autocorrelation. This make sense in terms of the previous knowledge on population differentiation in D. alata, especially if patterns are explained by a combination of isolation‐by‐distance and allelic surfing due to range expansion after the last glacial maximum. This reveals the potential for more complex applications of POM in population genetics. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 1152–1161.     

Spatial patterns of human gene frequencies in Europe

The aims of this study of spatial patterns of human gene frequencies in Europe are twofold. One is to present new methodology developed for the analysis of such data. The other is to report on the diversity of spatial patterns observed in Europe and their interpretation as evidence of population processes. Spatial variation in 59 allele and haplotype frequencies (26 genetic systems) for polymorphisms in blood antigens, enzymes, and proteins is analyzed for an aggregate of 3,384 localities, using homogeneity tests, one-dimensional and directional spatial correlograms, and SYMAP interpolated surfaces. The data matrices are reduced to reveal the principal patterns by clustering techniques. The findings of this study can be summarized as follows: 1) There is significant heterogeneity in allele frequencies among the localities for all but one genetic system. 2) There are significant spatial patterns for most allele frequencies. 3) There is a substantial minority of clinal patterns in these populations. Clinal trends are found more frequently in HLA alleles than for other variables. North-south and northwest-southwest gradients predominate. 4) There is a strong decline in overall genetic similarity with geographic distance for most variables. 5) There are few, if any, appreciable correlations in pairs of allele frequencies over the continent, and there is little interesting correlation structure in the resulting correlation matrix. 6) Few spatial correlograms are markedly similar to each other, yet they form well-defined clusters. Spatial variation patterns, therefore, differ among allele frequencies. Patterns of human gene frequencies in modern Europe are diverse and complex. No single model suffices for interpretation of the observed genetic structure. Some clinal patterns reported here support the Neolithic demic-expansion hypothesis, others suggest latitudinal selection. Most of the clinal patterns are in HLA alleles, but there is also evidence from ABO for east-west migration diffusion. The majority of patterns are patchy, consistent with hypotheses of isolation by distance or of settlement of genetically differing, subsequently expanding ethnic groups. While undoubtedly there has been an ongoing stochastic process of differentiation consistent with the isolation-by-distance model, this has not obscured the directional patterns caused by migration (demic diffusion), and has perhaps only reinforced the contribution from settlement of ethnic units to patterns of genetic variation. However, the impact of the latter is most difficult to discern and requires further methodological developments.     

Spatial structure of genetic variation and primary succession in the pioneer tree species Antirhea borbonica on La Réunion

In habitats where colonization and extinction are recurrent, the distribution of gene frequencies among patches of suitable habitat may reflect the age structure of different populations. In this study, we quantify population genetic structure for a pioneer tree species, Antirhea borbonica, in a chrono-sequence of primary succession on the lava flows of the Piton de La Fournaise volcano (La Réunion). Using microsatellite loci and amplified fragment length polymorphism (AFLP) markers, we quantified genetic variation within and among populations for early- and late-succession populations in a landscape where extinction and recolonization are recurrent (the ‘Grand Brûlé’) and for late-succession populations in a more stable landscape. This study produced three main results. First, we detected no evidence that founder events increase genetic differentiation among colonizing populations; F_ST values among early- and among late-succession populations were similar. Second, we found no evidence for isolation by distance; genetic distance was not correlated with spatial distance within and among populations. Third, F_IS values are consistently high in all populations, despite the fact that A. borbonica populations are functionally close to dioecy and thus expected to have an outcrossing mating system. Multiple colonization events from different sources may limit differentiation among young populations and spatial isolation may enhance differentiation among late-succession populations. Ecological processes acting during colonization may create the conditions for spatial aggregation within pioneer populations, and thus contribute to the high F_IS values.     

Gene Flow and the Geographical Distribution of a Molecular Polymorphism in DROSOPHILA PSEUDOOBSCURA

This paper discusses the relation between the geographical distribution of an enzyme polymorphism and population structure in Drosophila pseudoobscura. California populations of this species living in very different montane and lowland habitats separated by several kilometers are similar to each other in the frequency of an esterase allele. Previous estimates suggest that gene flow is too limited to account for this homogeneity of genetic structure, so that it must reflect some balancing force of natural selection. We show, however, that dispersal over unfavorable habitats is much greater than earlier supposed. Isolated populations of D. pseudoobscura separated by 15 km from other populations are subject to large amounts of immigration. This is shown by changes in the seasonal abundance of this species and in the annual pattern of lethal alleles in such populations. The genetic structure of an experimentally perturbed isolated population in an oasis returned to normal within a single year, suggesting that such populations are ephemeral and that the oasis is subject to annual recolonization by distant migrants. Direct assessment of marked flies shows that they can move at least 10 km in 24 hours over a desert. Such extensive gene flow may help explain the distribution of the esterase allele, and is relevant to the high level of molecular polymorphism and its general lack of geographic differentiation throughout the range of D. pseudoobscura.     

THE EFFECTS OF THE RELATIVE GEOGRAPHIC SCALES OF GENE FLOW AND SELECTION ON MORPH FREQUENCIES IN THE WALKING-STICK TIMEMA CRISTINAE

Gene frequencies in large populations are determined by a balance between selection and gene flow between neighborhoods of different selection regimes. This balance is affected by the area of the patches of a given selection regime relative to the gene-flow distance. If patches are small relative to gene-flow distance, similarity in the total area occupied by different patch types is a crucial condition for the stability of polymorphisms. However, if patches are larger than the gene-flow distance, then the relative area of different patch types is less important because of reduced gene flow resulting from isolation by distance. Two morphs (striped and unstriped) of the walking-stick Timema cristinae were each strongly associated with patches of distinct species of food plants on which they are most cryptic. The frequency of a morph was high on the plant on which it is most cryptic when either: (1) the area occupied by the food plant (patch) was very large; (2) the patch was completely isolated from other patches; or (3) the patch was larger than adjacent patches. Results (1) and (2) are consistent with isolation-by-distance models, and result (3) is consistent with Levene's multiple-niche polymorphism model.     

Lake Erie water snakes revisited: Morph- and age-specific variation in relative crypsis

Summary Populations of water snakes (Nerodia sipedon insularum) on islands in western Lake Erie are variable in colour pattern, consisting of unbanded, intermediate, and banded morphs. In contrast, mainland populations (N. s. sipedon) consist solely of banded morphs. Previous investigators hypothesized that natural selection favoured unbanded morphs on exposed island shorelines and banded morphs in overgrown mainland habitats and that gene flow from mainland populations was responsible for the persistence of banded morphs on islands. To clarify the potential role of natural selection, I quantified relative crypsis among morphs and age classes of water snakes by comparing the size of patches making up their colour patterns with the size of patches in island and mainland backgrounds. This analysis reveals that if unbanded morphs are more cryptic than intermediate and banded morphs on islands, it is only in the young-of-the-year age class. For older snakes on islands and for all snakes on the mainland, unbanded morphs are consistently less cryptic than intermediate and banded morphs. Given these results, the net direction of selection in island populations should depend on the intensity of predation on different age classes of snakes. Overall, selection may favour unbanded morphs (e.g. if predation occurs primarily on young-of-the-year), intermediate and banded morphs (e.g. if predation occurs primarily on older snakes), or be weak or absent (e.g. certain combinations of predation on young-of-the-year and older snakes). Using estimates of relative crypsis to guide reanalysis of morph frequency data, I find support for the hypothesis that unbanded morphs are favoured by natural selection in island populations.     

Rapid Homogenization of Multiple Sources: Genetic Structure of a Recolonizing Population of Fishers

Abstract: Fishers (Martes pennanti) were extirpated from much of southern Ontario, Canada, prior to the 1950s. We hypothesised that the recent recolonization of this area originated from an expansion of the population in Algonquin Provincial Park, which historically served as a refuge for fishers. To test this hypothesis, we created a sampling lattice to encompass Algonquin and the surrounding area, and we collected contemporaneous DNA samples. We sampled fishers from each of 35 sites and genotyped them at 16 microsatellite loci. Using a Bayesian assignment approach, with no a priori geographic information, we inferred 5 discrete genetic populations and used genetic population assignment as a means to cluster sites together. We concluded that the Algonquin Park fisher population has not been a substantial source for recolonization and expansion, which has instead occurred from a number of remnant populations within Ontario, Quebec, and most recently from the Adirondacks in New York, USA. The genetic structure among sampling sites across the entire area revealed a pattern of isolation-by-distance (IBD). However, an examination of the distribution of genetic structure (F_ST/1^- F_ST) at different distances showed higher rates of gene flow than predicted under a strict IBD model at small distances (40 km) within clusters and at larger distances up to 100 km among clusters. This pattern of genetic structure suggests increased migration and gene flow among expanding reproductive fronts.     

Geographic and genetic boundaries of brown bear (Ursus arctos) population in the Caucasus

The taxonomic status of brown bears in the Caucasus remains unclear. Several morphs or subspecies have been identified from the morphological (craniological) data, but the status of each of these subspecies has never been verified by molecular genetic methods. We analysed mitochondrial DNA sequences (control region) to reveal phylogenetic relationships and infer divergence time between brown bear subpopulations in the Caucasus. We estimated migration and gene flow from both mitochondrial DNA and microsatellite allele frequencies, and identified possible barriers to gene flow among the subpopulations. Our suggestion is that all Caucasian bears belong to the nominal subspecies of Ursus arctos. Our results revealed two genetically and geographically distinct maternal haplogroups: one from the Lesser Caucasus and the other one from the Greater Caucasus. The genetic divergence between these haplogroups dates as far back as the beginning of human colonization of the Caucasus. Our analysis of the least‐cost distances between the subpopulations suggests humans as a major barrier to gene flow. The low genetic differentiation inferred from microsatellite allele frequencies indicates that gene flow between the two populations in the Caucasus is maintained through the movements of male brown bears. The Likhi Ridge that connects the Greater and Lesser Caucasus mountains is the most likely corridor for this migration.