Dispersal ability and habitat requirements determine landscape‐level genetic patterns in desert aquatic insects

Species occupying the same geographic range can exhibit remarkably different population structures across the landscape, ranging from highly diversified to panmictic. Given limitations on collecting population‐level data for large numbers of species, ecologists seek to identify proximate organismal traits—such as dispersal ability, habitat preference and life history—that are strong predictors of realized population structure. We examined how dispersal ability and habitat structure affect the regional balance of gene flow and genetic drift within three aquatic insects that represent the range of dispersal abilities and habitat requirements observed in desert stream insect communities. For each species, we tested for linear relationships between genetic distances and geographic distances using Euclidean and landscape‐based metrics of resistance. We found that the moderate‐disperser Mesocapnia arizonensis (Plecoptera: Capniidae) has a strong isolation‐by‐distance pattern, suggesting migration–drift equilibrium. By contrast, population structure in the flightless Abedus herberti (Hemiptera: Belostomatidae) is influenced by genetic drift, while gene flow is the dominant force in the strong‐flying Boreonectes aequinoctialis (Coleoptera: Dytiscidae). The best‐fitting landscape model for M. arizonensis was based on Euclidean distance. Analyses also identified a strong spatial scale‐dependence, where landscape genetic methods only performed well for species that were intermediate in dispersal ability. Our results highlight the fact that when either gene flow or genetic drift dominates in shaping population structure, no detectable relationship between genetic and geographic distances is expected at certain spatial scales. This study provides insight into how gene flow and drift interact at the regional scale for these insects as well as the organisms that share similar habitats and dispersal abilities.     

ESTIMATION OF GENETIC NEIGHBORHOOD PARAMETERS FROM POLLEN AND SEED DISPERSAL IN THE MARINE ANGIOSPERM ZOSTERA MARINA L.

The relative importance of random genetic drift and local adaptation in causing population substructuring in plant species remains an important empirical question. Here I estimate the effective size of the genetic neighborhood, N_b, as a means of evaluating the likely role of genetic drift in creating genetic differentiation within a population of a marine plant, Zostera marina L. (eelgrass). Calculations of effective neighborhood size are based on field estimates of pollen and seed-dispersal distributions, an electrophoretic estimate of the mating system using open-pollinated progeny arrays, and determination of the effective density of reproductive individuals in the population. Neighborhood area calculated from the parent-offspring dispersal variances was equal to N_a = 524 m²; variance in the seed-dispersal distribution contributes more than twice as much as variance in pollen dispersal to N_a. Including an outcrossing rate slightly different from random, estimated neighborhood size for Z. marina is N_b = 6255. This estimate is one of the largest reported for plants or animals and indicates that genetic drift in small neighborhoods is highly unlikely to cause genetic substructuring in the study population. High gene-flow levels provided by the marine environment appear to prevent genetic isolation by distance among eelgrass patches, but the importance of drift through founder events in this population characterized by high patch turnover cannot be discounted and is the subject of ongoing study.     

Endangered species in small habitat patches can possess high genetic diversity: the case of the Tana River red colobus and mangabey

We used mtDNA sequence data from the Tana River red colobus and mangabey to determine how their population genetic structure was influenced by dispersal and habitat fragmentation. The colobus and mangabey are critically endangered primates endemic to gallery forests in eastern Kenya. The forests are a Pliocene–Pleistocene refugium that has recently undergone significant habitat loss and fragmentation due to human activities. We expected both primates to exhibit low levels of genetic diversity due to elevated genetic drift in their small populations, and to show a strong correspondence between genetic and geographic distance due to disruption of gene flow between forests by habitat fragmentation. Additionally, because mangabey females are philopatric, we expected their mtDNA variation to be homogeneous within forest patches but to be heterogeneous between patches. In contrast, colobus have a female-biased dispersal and so we expected their mtDNA variation to be homogeneous within and between forest patches. We found high levels of haplotype and nucleotide diversity as well as high levels of sequence divergence between haplotype groups in both species. The red colobus had significantly higher genetic variation than the mangabey did. Most of the genetic variation in both primates was found within forest fragments. Although both species showed strong inter-forest patch genetic structure we found no correspondence between genetic and geographic distances for the two primates. We attributed the high genetic diversity to recent high effective population size, and high sequence divergence and strong genetic structures to long-term habitat changes in the landscape.     

Population structure in the catfish Trichogenes longipinnis: drift offset by asymmetrical migration in a tiny geographic range

Based on population genetic theory and empirical studies of small populations, we expect that species with very small ranges (narrow endemics) will exhibit reduced genetic diversity, increasing their susceptibility to the negative effects of genetic homogeneity. Although this pattern of reduced diversity applies to most narrow endemics, conservation biologists have yet to identify a general pattern for the degree of spatial population genetic structure expected in species with very small ranges. In part, this is because the degree of population structure within narrow endemics will be highly variable depending on the equilibrium between the homogenizing effects of dispersal and the diversifying effects of drift and local selection in small populations, thus precluding general predictions about the relative importance of small range, small population sizes, and habitat patchiness for maintaining genetic diversity in narrowly-distributed species. We document a striking example of high population structure in the tiny geographic range of a stream-dwelling catfish, Trichogenes longipinnis , endemic to the Atlantic Forest of Brazil. The maintenance of this diversity results from a combination of asymmetrical and limited dispersal, and drift in small populations. Our results highlight the need to understand population structure, and not only overall genetic diversity, of narrowly-distributed species for their conservation planning.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 97 , 259–274.     

Diffusion approximations for one-locus multi-allele kin selection,mutation and random drift in group-structured populations: a unifying approach to selection models in population genetics

Diffusion approximations are ascertained from a two-time-scale argument in the case of a group-structured diploid population with scaled viability parameters depending on the individual genotype and the group type at a single multi-allelic locus under recurrent mutation, and applied to the case of random pairwise interactions within groups. The main step consists in proving global and uniform convergence of the distribution of the group types in an infinite population in the absence of selection and mutation, using a coalescent approach. An inclusive fitness formulation with coefficient of relatedness between a focal individual J affecting the reproductive success of an individual I, defined as the expected fraction of genes in I that are identical by descent to one or more genes in J in a neutral infinite population, given that J is allozygous or autozygous, yields the correct selection drift functions. These are analogous to the selection drift functions obtained with pure viability selection in a population with inbreeding. They give the changes of the allele frequencies in an infinite population without mutation that correspond to the replicator equation with fitness matrix expressed as a linear combination of a symmetric matrix for allozygous individuals and a rank-one matrix for autozygous individuals. In the case of no inbreeding, the mean inclusive fitness is a strict Lyapunov function with respect to this deterministic dynamics. Connections are made between dispersal with exact replacement (proportional dispersal), uniform dispersal, and local extinction and recolonization. The timing of dispersal (before or after selection, before or after mating) is shown to have an effect on group competition and the effective population size. In memory of Sam Karlin.     

Spatial genetic structuring in a vagile species, the European wood mouse

We examined the genetic structure of natural populations of the European wood mouse Apodemus sylvaticus at the microgeographic (<3 km) and macrogeographic (>30 km) scales. Ecological and behavioural studies indicate that this species exhibits considerable dispersal relative to its home-range size. Thus, there is potential for high gene flow over larger geographic areas. As levels of population genetic structure are related to gene flow, we hypothesized that population genetic structuring at the microgeographic level should be negligible, increasing only with geographic distance. To test this, four sites were sampled within a microgeographic scale with two additional samples at the macrogeographic level. Individuals ( n =415) were screened and analysed for seven polymorphic microsatellite loci. Contrary to our hypothesis, significant levels of population structuring were detected at both scales. Comparing genetic differentiation with geographic distance suggests increasing genetic isolation with distance. However, this distance effect was non-significant being confounded by surprisingly high levels of differentiation among microgeographic samples. We attribute this pattern of genetic differentiation to the effect of habitat fragmentation, splitting large populations into components with small effective population sizes resulting in enhanced genetic drift. Our results indicate that it is incorrect to assume genetic homogeneity among populations even where there is no evidence of physical barriers and dispersal can occur freely. In the case of A. sylvaticus , it is not clear whether dispersal does not occur across habitat barriers or behavioural dispersal occurs without consequent gene flow.     

Unpredictable selection in a structured population leads to local genetic differentiation in evolved reaction norms

Unpredictability during development of the optimum phenotype under future selection leads to a compromise reaction norm with a slope that is shallower than the slope of the optimum reaction norm. Unpredictability of selection can lead to an evolved curved reaction norm when genetic variation for curvature is available even if the optimum reaction norm is linear. This requires asymmetry in the frequency distribution of the habitats of selection; at small population size, stochasticity in the number of individuals per selection habitat is sufficient to generate such asymmetry. Unpredictability of selection in structured populations leads to local genetic differentiation of reaction norms. The mean habitat of a subpopulation is defined as the subpopulation's focal habitat. The evolved mean reaction norm of each subpopulation is anchored at the optimum genotypic value in its focal habitat. Linear reaction norms are parallel if the conditional distribution of adults around the focal habitats is the same for each subpopulation. Adult migration and absence of zygote dispersal represents the ultimate structured population, each habitat playing the role of focal habitat. Absence of zygote dispersal requires that the flow of individuals through the habitats is used instead of the habitats’ frequencies in the prediction of the evolved reaction norm. Adult migration in absence of zygote dispersal leads to an evolved pattern of locally differentiated reaction norms with optimum genotypic value anchored in the focal habitat and, for linear reaction norms, parallel slopes.     

Leptokurtic pollen-flow,non-leptokurtic gene-flow in a wind-pollinated herb,Plantago lanceolata L.

Summary The purpose of this study was to simultaneously measure pollen dispersal distance and actual pollen-mediated gene-flow distance in a wind-pollinated herb, Plantago lanceolata. The pollen dispersal distribution, measured as pollen deposition in a wind tunnel, is leptokurtic, as expected from previous studies of wind-pollinated plants. Gene-flow, measured as seeds produced on rows of male-sterile inflorescences in the wind tunnel, is non-leptokurtic, peaking at an intermediate distance. The difference between the two distributions results from the tendency of the pollen grains to cluster. These pollen clusters are the units of gene dispersal, with clusters of intermediate and large size contributing disproportionately to gene-flow. Since many wind-pollinated species show pollen clustering (see text), the common assumption for wind-pollinated plants that gene-flow is leptokurtic requires re-examination. Gene-flow was also measured in an artifical outdoor population of male-steriles, containing a single pollen source plant in the center of the array. The gene flow distribution is significantly platykurtic, and has the same general properties outdoors, where wind speed and turbulence are uncontrolled, as it does in the wind tunnel. I estimated genetic neighborhood size based on my measure of gene-flow in the outdoor population. The estimate shows that populations of Plantago lanceolata will vary in effective number from a few tens of plants to more than five hundred plants, depending on the density of the population in question. Thus, the measured pollen-mediated gene-flow distribution and population density will interact to produce effective population sizes ranging from those in which there is no random genetic drift to those in which random genetic drift plays an important role in determining gene frequencies within and among populations. Despite the platykurtosis in the distribution, pollen-mediated gene dispersal distances are still quite limited, and considerable within and among-population genetic differentiation is to be expected in this species.     

Individual variation in dispersal associated with phenotype influences fine-scale genetic structure in weasels

In general, landscape genetic studies have ignored the potential role that the phenotype of individuals plays in determining fine-scale genetic structure in species. This potential over-simplification ignores an important component that dispersal is both condition- and phenotype-dependent. In order to investigate the relationship between potential dispersal, habitat selection and phenotype, we examined the spatial ecology, body mass and fine-scale genetic structure of weasels (Mustela nivalis) in Bia?owie?a Forest in Poland. Our study population is characterized by an almost three-fold phenotypic variation in adult body mass and weasels were segregated in certain habitats according to size. We detected significant genetic structuring associated with habitat within the studied area and analyses of radio-tracking and re-capture data showed that the maximum extent of movement was achieved by weasels of medium body size, whereas the smallest and largest individuals exhibited higher site fidelity. With the unrestricted movement of the medium-sized individuals across optimal habitat, genetic admixture does occur. However, the presence of a barrier leads to unidirectional gene flow, with larger individuals outcompeting smaller individuals and therefore maintaining the genetic break in the study area. This highlights the importance of considering both intrinsic (phenotype) and extrinsic (environmental) factors in understanding dispersal patterns and ultimately, gene flow in complex landscapes.     

Body Size and Geographic Range Do Not Explain Long Term Variation in Fish Populations: A Bayesian Phylogenetic Approach to Testing Assembly Processes in Stream Fish Assemblages

We combine evolutionary biology and community ecology to test whether two species traits, body size and geographic range, explain long term variation in local scale freshwater stream fish assemblages. Body size and geographic range are expected to influence several aspects of fish ecology, via relationships with niche breadth, dispersal, and abundance. These traits are expected to scale inversely with niche breadth or current abundance, and to scale directly with dispersal potential. However, their utility to explain long term temporal patterns in local scale abundance is not known. Comparative methods employing an existing molecular phylogeny were used to incorporate evolutionary relatedness in a test for covariation of body size and geographic range with long term (1983 – 2010) local scale population variation of fishes in West Fork White River (Indiana, USA). The Bayesian model incorporating phylogenetic uncertainty and correlated predictors indicated that neither body size nor geographic range explained significant variation in population fluctuations over a 28 year period. Phylogenetic signal data indicated that body size and geographic range were less similar among taxa than expected if trait evolution followed a purely random walk. We interpret this as evidence that local scale population variation may be influenced less by species-level traits such as body size or geographic range, and instead may be influenced more strongly by a taxon’s local scale habitat and biotic assemblages.     

Microhabitat contributes to microgeographic divergence in threespine stickleback

Since the New Synthesis, most migration-selection balance theory has predicted that there should be negligible differentiation over small spatial scales (relative to dispersal), because gene flow should erode any effect of divergent selection. Nevertheless, there are classic examples of microgeographic divergence, which theory suggests can arise under specific conditions: exceptionally strong selection, phenotypic plasticity in philopatric individuals, or nonrandom dispersal. Here, we present evidence of microgeographic morphological variation within lake and stream populations of threespine stickleback (Gasterosteus aculeatus). It seems reasonable to assume that a given lake or stream population of fish is well-mixed. However, we found this assumption to be untenable. We examined trap-to-trap variation in 34 morphological traits measured on stickleback from 16 lakes and 16 streams. Most traits varied appreciably among traps within populations. Both between-trap distance and microhabitat characteristics such as depth and substrate explained some of the within-population morphological variance. Microhabitat was also associated with genotype at particular loci but there was no genetic isolation by distance, implying that heritable habitat preferences may contribute to microgeographic variation. Our study adds to growing evidence that microgeographic divergence can occur across small spatial scales within individuals’ daily dispersal neighborhood where gene flow is expected to be strong.     

Geographical Variation in a Quantitative Character

A model for the evolution of the local averages of a quantitative character under migration, selection, and random genetic drift in a subdivided population is formulated and investigated. Generations are discrete and nonoverlapping; the monoecious, diploid population mates at random in each deme. All three evolutionary forces are weak, but the migration pattern and the local population numbers are otherwise arbitrary. The character is determined by purely additive gene action and a stochastically independent environment; its distribution is Gaussian with a constant variance; and it is under Gaussian stabilizing selection with the same parameters in every deme. Linkage disequilibrium is neglected. Most of the results concern the covariances of the local averages. For a finite number of demes, explicit formulas are derived for (i) the asymptotic rate and pattern of convergence to equilibrium, (ii) the variance of a suitably weighted average of the local averages, and (iii) the equilibrium covariances when selection and random drift are much weaker than migration. Essentially complete analyses of equilibrium and convergence are presented for random outbreeding and site homing, the Levene and island models, the circular habitat and the unbounded linear stepping-stone model in the diffusion approximation, and the exact unbounded stepping-stone model in one and two dimensions.     

The role of environment and core‐margin effects on range‐wide phenotypic variation in a montane grasshopper

The integration of genetic information with ecological and phenotypic data constitutes an effective approach to gain insight into the mechanisms determining interpopulation variability and the evolutionary processes underlying local adaptation and incipient speciation. Here, we use the Pyrenean Morales grasshopper (Chorthippus saulcyi moralesi) as study system to (i) analyse the relative role of genetic drift and selection in range‐wide patterns of phenotypic differentiation and (ii) identify the potential selective agents (environment, elevation) responsible for variation. We also test the hypothesis that (iii) the development of dispersal‐related traits is associated with different parameters related to population persistence/turnover, including habitat suitability stability over the last 120 000 years, distance to the species distribution core and population genetic variability. Our results indicate that selection shaped phenotypic differentiation across all the studied morphological traits (body size, forewing length and shape). Subsequent analyses revealed that among‐population differentiation in forewing length was significantly explained by a temperature gradient, suggesting an adaptive response to thermoregulation or flight performance under contrasting temperature regimes. We found support for our hypothesis predicting a positive association between the distance to the species distribution core and the development of dispersal‐related morphology, which suggests an increased dispersal capability in populations located at range edges that, in turn, exhibit lower levels of genetic variability. Overall, our results indicate that range‐wide patterns of phenotypic variation are partially explained by adaptation in response to local environmental conditions and differences in habitat persistence between core and peripheral populations.     

Evidence of spatial genetic structure in a California bunchgrass population

We investigated the scale of genetic variation of purple needlegrass (Nassella pulchra), a species commonly used in California for grassland restoration. Common garden and field data revealed evidence of genetic differentiation between two intermixed microhabitats characterized by differences in soil depth and community composition. We assessed the genetic variation within a single population using randomly amplified polymorphic DNA (RAPD) data collected from clusters of five individuals in 40 locations. We found no evidence for genetic structure at the whole population level. At smaller spatial scales, however, we found strong evidence that genetic subdivision of the population occurs at the level of the maternal neighborhood. We suggest that the interaction between widespread pollen dispersal and restricted seed dispersal may be the primary factor generating these results; panmictic pollen dispersal will make detection of genetic patterning difficult at larger spatial scales while limited seed dispersal will generate local genetic structure. As a result, the detection of population genetic structure will depend on the spatial scale of analysis. Local selection gradients related to topography and soil depth are also likely to play a role in structuring local genetic variation. Since N. pulchra is widely used in California in grassland and woodland habitat restoration, we suggest that, as a general rule, care should be exercised in transferring germplasm for the purposes of conservation when little is known about the within-population genetic subdivision of a plant species. Received: 23 December 1996 / Accepted: 20 May 1997     

Partitioning the effects of isolation by distance,environment, and physical barriers on genomic divergence between parapatric threespine stickleback

Genetic divergence between populations is shaped by a combination of drift, migration, and selection, yielding patterns of isolation‐by‐distance (IBD) and isolation‐by‐environment (IBE). Unfortunately, IBD and IBE may be confounded when comparing divergence across habitat boundaries. For instance, parapatric lake and stream threespine stickleback (Gasterosteus aculeatus) may have diverged due to selection against migrants (IBE), or mere spatial separation (IBD). To quantitatively partition the strength of IBE and IBD, we used recently developed population genetic software (BEDASSLE) to analyze partial genomic data from three lake‐stream clines on Vancouver Island. We find support for IBD within each of three outlet streams (unlike prior studies of lake‐stream stickleback). In addition, we find evidence for IBE (controlling for geographic distance): the genetic effect of habitat is equivalent to geographic separation of ～1.9 km of IBD. Remarkably, of our three lake‐stream pairs, IBE is strongest where migration between habitats is easiest. Such microgeographic genetic divergence would require exceptionally strong divergent selection, which multiple experiments have failed to detect. Instead, we propose that nonrandom dispersal (e.g., habitat choice) contributes to IBE. Supporting this conclusion, we show that the few migrants between habitats are a nonrandom subset of the phenotype distribution of the source population.     

Microgeographic population structure of green swordail fish: genetic differentiation despite abundant migration

Swordtails (Xiphophorus; Poeciliidae) have figured prominently in research on fish mating behaviours, sexual selection, and carcinogenesis, but their population structures and dispersal patterns have been relatively neglected. Using nine microsatellite loci, we estimated genetic differentiation in Xiphophorus helleri within and between adjacent streams in Belize. The genetic data were complemented by a tagging study of movement within one stream. In the absence of physical dispersal barriers (waterfalls), population structure followed an isolation by distance (IBD) pattern. Genetic differentiation (F_ST up to 0.07) was significant between and within creeks, despite high dispersal in the latter as judged by the tagging data. Such heterogeneity apparently was a result of genetic drift in local demes, due to small population sizes and highly skewed paternity. The IBD pattern was interrupted by waterfalls, boosting F_ST above 0.30 between adjacent samples across these barriers. Overall, our results are helpful in understanding the interplay of evolutionary forces and population dynamics in a small fish living in a changeable habitat.     

Population biology of Avena

Summary Pollen and seed dispersal patterns were analyzed in both natural and experimental populations of Avena barbata. Localized estimates of gene flow rates and plant densities gave estimates of neighborhood size in the range of 40 to 400 plants; the estimates of mean rate and distance of gene flow seemed to vary widely due to variable wind direction, rodent activity, microsite heterogeneity, etc. The relative sizes of neighborhoods in several populations were correlated with the patchy distribution of different genotypes (scored for lemma color and leaf sheath hairiness) within short distances, but patch sizes had a wide range among different sites. Highly localized gene flow patterns seemed to account for the observed pattern of highly patchy variation even when the dispersal curves for both pollen and seed were platykurtic in many cases. Measures of the stability of patches in terms of their size, dispersion in space and genetic structure in time are needed in order to sort out the relative roles of founder effects, random drift (due to small neighborhood size), and highly localized selection. However, our observations suggest that many variables and stochastic processes are involved in such studies so as to allow only weak inference about the underlying role of natural selection, drift and factors of population regulatien.     

Genetic drift widens the expected cline but narrows the expected cline width

Random genetic drift shifts clines in space, alters their width, and distorts their shape. Such random fluctuations complicate inferences from cline width and position. Notably, the effect of genetic drift on the expected shape of the cline is opposite to the naive (but quite common) misinterpretation of classic results on the expected cline. While random drift on average broadens the overall cline in expected allele frequency, it narrows the width of any particular cline. The opposing effects arise because locally, drift drives alleles to fixation--but fluctuations in position widen the expected cline. The effect of genetic drift can be predicted from standardized variance in allele frequencies, averaged across the habitat: . A cline maintained by spatially varying selection (step change) is expected to be narrower by a factor of √1- relative to the cline in the absence of drift. The expected cline is broader by the inverse of this factor. In a tension zone maintained by underdominance, the expected cline width is narrower by about 1- relative to the width in the absence of drift. Individual clines can differ substantially from the expectation, and we give quantitative predictions for the variance in cline position and width. The predictions apply to clines in almost one-dimensional circumstances such as hybrid zones in rivers, deep valleys, or along a coast line and give a guide to what patterns to expect in two dimensions.     

Coyotes demonstrate how habitat specialization by individuals of a generalist species can diversify populations in a heterogeneous ecoregion

The tendency for individuals to disperse into habitat similar to their natal habitat has been observed in a wide range of species, although its population genetic consequences have received little study. Such behavior could lead to discrete habitat-specific population subdivisions even in the absence of physical dispersal barriers or habitat gaps. Previous studies of coyotes have supported this hypothesis in a small region of California, but its evolutionary significance ultimately depends on the extent and magnitude of habitat-specific subdivision. Here, we investigated these questions using autosomal, Y chromosome, and mitochondrial markers and >2,000 coyotes from a broad region, including 2 adjacent ecoregions with contrasting levels of habitat heterogeneity--the California Floristic Province (CFP) (heterogeneous landscape) and the Desert-Prairie ecoregion (DPE) (homogeneous landscape). Consistent with predictions, we found a close correspondence between population genetic structure and habitat subdivisions throughout the CFP and virtual panmixia over the larger DPE. Conversely, although genetic diversity was similar in these 2 ecoregions overall, it was lower within sites of the CFP, as would be the expected consequence of greater genetic drift within subregions. The magnitude of habitat-specific genetic subdivisions (i.e., genetic distance) in the CFP varied considerably, indicating complexity (e.g., asymmetric gene flow or extinction/recolonization), but, in general, was higher than that due to geographic distance or recent human-related barriers. Because habitat-specific structure can enhance a species' adaptive potential and resilience to changing environments, these findings suggest the CFP may constitute an evolutionarily important portion of the range for coyotes and sympatric species exhibiting habitat-specific population structure.