Secondary contact and changes in coastal habitat availability influence the nonequilibrium population structure of a salmonid (Oncorhynchus keta)

Numerous empirical studies have reported lack of migration–drift equilibrium in wild populations. Determining the causes of nonequilibrium population structure is challenging because different evolutionary processes acting at a variety of spatiotemporal scales can produce similar patterns. Studies of contemporary populations in northern latitudes suggest that nonequilibrium population structure is probably caused by recent colonization of the region after the last Pleistocene ice age ended ~13 000 years ago. The chum salmon's (Oncorhynchus keta) range was fragmented by dramatic environmental changes during the Pleistocene. We investigated the population structure of chum salmon on the North Alaska Peninsula (NAP) and, using both empirical data and simulations, evaluated the effects of colonization timing and founder population heterogeneity on patterns of genetic differentiation. We screened 161 single nucleotide polymorphisms and found evidence of nonequilibrium population structure when the slope of the isolation‐by‐distance relationship was examined at incremental spatial scales. In addition, simulations suggested that this pattern closely matched models of recent colonization of the NAP by secondary contact. Our results agree with geological and archaeological data indicating that the NAP was a dynamic landscape that may have been more recently colonized than during the last deglaciation because of dramatic changes in coastal hydrology over the last several thousand years.     

Rapid genetic restoration of a keystone species exhibiting delayed demographic response

Genetic founder effects are often expected when animals colonize restored habitat in fragmented landscapes, but empirical data on genetic responses to restoration are limited. We examined the genetic response of banner‐tailed kangaroo rats (Dipodomys spectabilis) to landscape‐scale grassland restoration in the Chihuahuan Desert of New Mexico, USA. Dipodomys spectabilis is a grassland specialist and keystone species. At sites treated with herbicide to remove shrubs, colonization by D. spectabilis is slow and populations persist at low density for ≥10 years (≥6 generations). Persistence at low density and low gene flow may cause strong founder effects. We compared genetic structure of D. spectabilis populations between treated sites and remnant grasslands, and we examined how the genetic response to restoration depended on treatment age, area, and connectivity to source populations. Allelic richness and heterozygosity were similar between treated sites and remnant grasslands. Allelic richness at treated sites was greatest early in the restoration trajectory, and genetic divergence did not differ between recently colonized and established populations. These results indicated that founder effects during colonization of treated sites were weak or absent. Moreover, our results suggested founder effects were not mitigated by treatment area or connectivity. Dispersal is negatively density‐dependent in D. spectabilis, and we hypothesize that high gene flow may occur early in the restoration trajectory when density is low. Our study shows genetic diversity can be recovered more rapidly than demographic components of populations after habitat restoration and that founder effects are not inevitable for animals colonizing restored habitat in fragmented landscapes.     

Clonal Diversity and Fine‐scale Genetic Structure in a High Andean Treeline Population

Clonal propagation becomes more abundant with increasing altitudes as environmental conditions worsen. To date, little attention has been paid to the way in which clonal propagation affects genetic diversity and the fine‐scale spatial genetic structure (FSGS) of clonal alpine trees. An AFLP study was undertaken to quantify the clonal and genetic diversity and FSGS of the vulnerable treeline species Polylepis reticulata in Ecuador. We successfully genotyped 32 and 75 ramets within 4 m × 100 m (coarse scale) and 4 m × 4 m (fine scale) transects of one population, respectively. Higher genotypic diversity was detected at the coarse scale than at the fine scale, while lower genetic diversity was detected for P. reticulata than other Polylepis spp. at both scales. Significantly stronger FSGS was detected at the ramet level than the genet level for P. reticulata within a spatial distance of 3 m. The studied P. reticulata population showed pronounced FSGS (Sp = 0.012 at the genet level, a statistic reflecting declining pairwise kinship with distance) revealed restricted gene dispersal, which implies restricted seed dispersal for this population, assuming pollen flow is as extensive as that described for other wind‐pollinated tree species. Our results revealed that clonal diversity is a function of both sample size and the spatial scale of the sampling area. The findings highlights that clonal propagation has affected FSGS within a spatial distance of 3 m for this species.     

Genetic signatures of a range expansion in natura: when clones play leapfrog

The genetic consequences of range expansions have generally been investigated at wide geographical and temporal scales, long after the colonization event. A unique ecological system enabled us to both monitor the colonization dynamics and decipher the genetic footprints of expansion over a very short time period. Each year an epidemic of the poplar rust (Melampsora larici‐populina) expands clonally and linearly along the Durance River, in the Alps. The colonization dynamics observed in 2004 showed two phases with different genetic outcomes. Upstream, fast colonization maintained high genetic diversity. Downstream, the colonization wave progressively faltered, diversity eroded, and differentiation increased, as expected under recurrent founder events. In line with the high dispersal abilities of rust pathogens, we provide evidence for leapfrog dispersal of clones. Our results thus emphasize the importance of colonization dynamics in shaping spatial genetic structure in the face of high gene flow.     

Short‐term variations in gene flow related to cyclic density fluctuations in the common vole

In highly fluctuating populations with complex social systems, genetic patterns are likely to vary in space and time due to demographic and behavioural processes. Cyclic rodents are extreme examples of demographically instable populations that often exhibit strong social organization. In such populations, kin structure and spacing behaviour may vary with density fluctuations and impact both the composition and spatial structure of genetic diversity. In this study, we analysed the multiannual genetic structure of a cyclic rodent, Microtus arvalis, using a sample of 875 individuals trapped over three complete cycles (from 1999 to 2007) and genotyped at 10 microsatellite loci. We tested the predictions that genetic diversity and gene flow intensity vary with density fluctuations. We found evidences for both spatial scale‐dependant variations in genetic diversity and higher gene flow during high density. Moreover, investigation of sex‐specific relatedness patterns revealed that, although dispersal is biased toward males in this species, distances moved by both sexes were lengthened during high density. Altogether, these results suggest that an increase in migration with density allows to restore the local loss of genetic diversity occurring during low density. We then postulate that this change in migration results from local competition, which enhances female colonization of empty spaces and male dispersal among colonies.     

Maintaining genetic diversity and population panmixia through dispersal and not gene flow in a holocyclic heteroecious aphid species

Heteroecious holocyclic aphids exhibit both sexual and asexual reproduction and alternate among primary and secondary hosts. Most of these aphids can feed on several related hosts, and invasions to new habitats may limit the number of suitable hosts. For example, the aphid specialist Aphis glycines survives only on the primary host buckthorn (Rhamnus spp.) and the secondary host soybean (Glycine max) in North America where it is invasive. Owing to this specialization and sparse primary host distribution, host colonization events could be localized and involve founder effects, impacting genetic diversity, population structure and adaptation. We characterized changes in the genetic diversity and structure across time among A. glycines populations. Populations were sampled from secondary hosts twice in the same geographical location: once after secondary colonization (early season), and again immediately before primary host colonization (late season). We tested for evidence of founder effects and genetic isolation in early season populations, and whether or not late-season dispersal restored genetic diversity and reduced fragmentation. A total of 24 single-nucleotide polymorphisms and 6 microsatellites were used for population genetic statistics. We found significantly lower levels of genotypic diversity and more genetic isolation among early season collections, indicating secondary host colonization occurred locally and involved founder effects. Pairwise F(ST) decreased from 0.046 to 0.017 in early and late collections, respectively, and while genetic relatedness significantly decreased with geographical distance in early season collections, no spatial structure was observed in late-season collections. Thus, late-season dispersal counteracts the secondary host colonization through homogenization and increases genetic diversity before primary host colonization.     

Spatial genetic structure in a metapopulation of the land snail Cepaea nemoralis (Gastropoda: Helicidae)

Habitat fragmentation is a major force affecting demography and genetic structure of wild populations, especially in agricultural landscapes. The land snail Cepaea nemoralis (L.) was selected to investigate the impact of habitat fragmentation on the spatial genetic structure of an organism with limited dispersal ability. Genetic and morphological patterns were investigated at a local scale of a 500 m transect and a mesoscale of 4 x 4 km in a fragmented agricultural landscape while accounting for variation in the landscape using least-cost models. Analysis of microsatellite loci using expected heterozygosity (HE), pairwise genetic distance (FST/1-FST) and spatial autocorrelograms (Moran's I) as well as shell characteristics revealed spatial structuring at both scales and provided evidence for a metapopulation structure. Genetic diversity was related to morphological diversity regardless of landscape properties. This pointed to bottlenecks caused by founder effects after (re)colonization. Our study suggests that metapopulation structure depended on both landscape features and the shape of the dispersal function. A range of genetic spatial autocorrelation up to 80 m at the local scale and up to 800 m at the mesoscale indicated leptokurtic dispersal patterns. The metapopulation dynamics of C. nemoralis resulted in a patchwork of interconnected, spatially structured subpopulations. They were shaped by gene flow which was affected by landscape features, the dispersal function and an increasing role of genetic drift with distance.     

Investigating past range dynamics for a weed of cultivation,Silene vulgaris

Since the last glacial maximum (LGM), many plant and animal taxa have expanded their ranges by migration from glacial refugia. Weeds of cultivation may have followed this trend or spread globally following the expansion of agriculture or ruderal habitats associated with human‐mediated disturbance. We tested whether the range expansion of the weed Silene vulgaris across Europe fit the classical model of postglacial expansion from southern refugia, or followed known routes of the expansion of human agricultural practices. We used species distribution modeling to predict spatial patterns of postglacial expansion and contrasted these with the patterns of human agricultural expansion. A population genetic analysis using microsatellite loci was then used to test which scenario was better supported by spatial patterns of genetic diversity and structure. Genetic diversity was highest in southern Europe and declined with increasing latitude. Locations of ancestral demes from genetic cluster analysis were consistent with areas of predicted refugia. Species distribution models showed the most suitable habitat in the LGM on the southern coasts of Europe. These results support the typical postglacial northward colonization from southern refugia while refuting the east‐to‐west agricultural spread as the main mode of expansion for S. vulgaris. We know that S. vulgaris has recently colonized many regions (including North America and other continents) through human‐mediated dispersal, but there is no evidence for a direct link between the Neolithic expansion of agriculture and current patterns of genetic diversity of S. vulgaris in Europe. Therefore, the history of range expansion of S. vulgaris likely began with postglacial expansion after the LGM, followed by more recent global dispersal by humans.     

Long‐distance dispersal and genetic structure of natural populations: an assessment of the inverse isolation hypothesis in peat mosses

It is well accepted that the shape of the dispersal kernel, especially its tail, has a substantial effect on the genetic structure of species. Theory predicts that dispersal by fat‐tailed kernels reshuffles genetic material, and thus, preserves genetic diversity during colonization. Moreover, if efficient long‐distance dispersal is coupled with random colonization, an inverse isolation effect is predicted to develop in which increasing genetic diversity per colonizer is expected with increasing distance from a genetically variable source. By contrast, increasing isolation leads to decreasing genetic diversity when dispersal is via thin‐tailed kernels. Here, we use a well‐established model group for dispersal biology (peat mosses: genus Sphagnum) with a fat‐tailed dispersal kernel, and the natural laboratory of the Stockholm archipelago to study the validity of the inverse isolation hypothesis in spore‐dispersed plants in island colonization. Population genetic structure of three species (Sphagnum fallax, Sphagnum fimbriatum and Sphagnum palustre) with contrasting life histories and ploidy levels were investigated on a set of islands using microsatellites. Our data show (, amova , IBD) that dispersal of the two most abundant species can be well approximated by a random colonization model. We find that genetic diversity per colonizer on islands increases with distance from the mainland for S. fallax and S. fimbriatum. By contrast, S. palustre deviates from this pattern, owing to its restricted distribution in the region, affecting its source pool strength. Therefore, the inverse isolation effect appears to hold in natural populations of peat mosses and, likely, in other organisms with small diaspores.     

Genetic and phenotypic divergence in an island bird: isolation by distance,by colonization or by adaptation?

Discerning the relative roles of adaptive and nonadaptive processes in generating differences among populations and species, as well as how these processes interact, is a fundamental aim in biology. Both genetic and phenotypic divergence across populations can be the product of limited dispersal and gradual genetic drift across populations (isolation by distance), of colonization history and founder effects (isolation by colonization) or of adaptation to different environments preventing migration between populations (isolation by adaptation). Here, we attempt to differentiate between these processes using island populations of Berthelot's pipit (Anthus berthelotii), a passerine bird endemic to three Atlantic archipelagos. Using microsatellite markers and approximate Bayesian computation, we reveal that the northward colonization of this species ca. 8500 years ago resulted in genetic bottlenecks in the colonized archipelagos. We then show that high levels of genetic structure exist across archipelagos and that these are consistent with a pattern of isolation by colonization, but not with isolation by distance or adaptation. Finally, we show that substantial morphological divergence also exists and that this is strongly concordant with patterns of genetic structure and bottleneck history, but not with environmental differences or geographic distance. Overall, our data suggest that founder effects are responsible for both genetic and phenotypic changes across archipelagos. Our findings provide a rare example of how founder effects can persist over evolutionary timescales and suggest that they may play an important role in the early stages of speciation.     

Strong genetic differentiation due to multiple founder events during a recent range expansion of an introduced wall lizard population

Biological invasions represent ideal systems for the study of evolutionary processes associated with colonization events. It has been hypothesized that the genetic diversity is generally decreasing from the centre of the range to the margins due to multiple founder events. Invasive populations offer the opportunity to test this hypothesis at a fine spatial and temporal scale. We analysed the genetic structure of a large expanding non-native population of the Common Wall Lizard (Podarcis muralis) in Passau (Germany) using thirteen microsatellite loci. We analyzed the genetic structure and levels of admixture across a transect reflecting the expansion process and tested for a loss of genetic diversity and an increase of genetic differentiation from the centre to the invasion front. Our results demonstrate that significant genetic population structure can emerge rapidly at a small spatial scale. We found a trend for an increase in genetic differentiation and a decrease in genetic diversity from the invasion centre to the expanding range margin, suggesting that genetic drift is the major factor causing this pattern. The correlation between genetic diversity and average genetic differentiation was significant among sites. We hypothesize that the territoriality of P. muralis generates sufficient rates of noncontiguous and stratified dispersal from longer established sites to maintain significant genetic diversity at the invasion front. Simultaneously, territoriality might restrict the colonization success of migrants at established sites, so that in combination with founder events a strong differentiation arises.     

Population genetic structure and male-biased dispersal in the queenless ant Diacamma cyaneiventre

In this study we investigated the population genetic structure of the queenless ant Diacamma cyaneiventre. This species, lacking winged queens, is likely to have a restricted female dispersal. We used both mitochondrial and microsatellite markers to assess the consequence of such restricted female dispersal at three geographical scales: within a given locality (< 1 km), between localities within a given region (< 10 km) and between regions (> 36 km). Within a locality, a strong population structure was observed for mitochondrial DNA (mtDNA) whereas weak or nonexistent population genetic structure was observed for the microsatellites (around 5% of the value for mtDNA). Male gene flow was estimated to be about 20-30 times higher than female gene flow at this scale. At a larger spatial scale, very strong genetic differentiation for both markers was observed between localities - even within a single region. Female dispersal is nonexistent at these scales and male dispersal is very restricted, especially between regions. The phylogeographical structure of the mtDNA haplotypes as well as the very low genetic diversity of mtDNA within localities indicate that new sites are colonized by a single migration event from adjacent localities, followed by successive colony fissions. These patterns of genetic variability and differentiation agree with what is theoretically expected when colonization events are kin-structured and when, following colonization, dispersion is mainly performed by males.     

Population genomic analysis suggests strong influence of river network on spatial distribution of genetic variation in invasive saltcedar across the southwestern United States

Understanding the complex influences of landscape and anthropogenic elements that shape the population genetic structure of invasive species provides insight into patterns of colonization and spread. The application of landscape genomics techniques to these questions may offer detailed, previously undocumented insights into factors influencing species invasions. We investigated the spatial pattern of genetic variation and the influences of landscape factors on population similarity in an invasive riparian shrub, saltcedar (Tamarix L.) by analysing 1,997 genomewide SNP markers for 259 individuals from 25 populations collected throughout the southwestern United States. Our results revealed a broad‐scale spatial genetic differentiation of saltcedar populations between the Colorado and Rio Grande river basins and identified potential barriers to population similarity along both river systems. River pathways most strongly contributed to population similarity. In contrast, low temperature and dams likely served as barriers to population similarity. We hypothesize that large‐scale geographic patterns in genetic diversity resulted from a combination of early introductions from distinct populations, the subsequent influence of natural selection, dispersal barriers and founder effects during range expansion.     

Pursuing glacier retreat: genetic structure of a rapidly expanding Larix decidua population

One of the greatest threats to the long-term viability of migrating plant species is the loss of genetic diversity due to founder effects. Populations can expand as a response to climate change, but it is uncertain if long-lived plant species can maintain sufficient genetic diversity at the leading edge of migrating populations. This study uses an expanding Larix decidua population investigated along a chronosequence at landscape (350 ha) and local (0.8 ha) scales to test whether accelerated migration as a result of climate warming has the potential to intensify genetic erosion. Nine SSR markers revealed similar genetic diversity among eight sub-populations along the chronosequence (overall H(e) = 0.73; SE=0.04). Sub-populations were not genetically differentiated and all sampled individuals (N=730) formed one major genetic cluster indicating homogenizing gene flow despite spatial genetic structure (SGS) up to 80 m. At the local scale, individuals at the leading edge [early successional sub-population (ESSP), N =140] and a sub-population at equilibrium [late successional sub-population (LSSP), N = 290] revealed high genetic diversity in largest-sized cohorts. SGS among juveniles occurred up to 30 m in LSSP but there was no structure in ESSP. Accordingly, a maximum likelihood paternity assignment revealed local gene dispersal in LSSP (2-48 m) and intermediate-to-long distance dispersal into ESSP (115-3132 m). The findings indicate intensive mixing of the genes in this expanding population instead of founder effects and support the view that genetic diversity can be maintained in a long-lived species during rapid population expansion driven by climate warming.     

Genetic diversity of melon aphids Aphis gossypii associated with landscape features

Despite increasing evidence that landscape features strongly influence the abundance and dispersal of insect populations, landscape composition has seldom been explicitly linked to genetic structure. We conducted a genetic study of the melon aphid, Aphis gossypii, in two counties of Beijing, China during spring migration using samples from watermelon. We performed aphid genetic analysis using restriction site associated DNA sequencing (2b‐RAD) and investigated the relationship between land cover and the genetic diversity. The percentage area of land cover (cropland, vegetable, orchard, grassland, woodland) was quantified in each particular scale (ranging from 0.5 km to 3 km) and was used as a predictor variable in our generalized linear models. We found a moderate level of genetic differentiation among nine sampled populations. Geographic distance and genetic distance were not significantly associated, indicating that geographic location was not a barrier to migration. These nine populations could be clustered depending on their level of genetic diversity (high and low). The genetic diversity (Shannon’s information index) was positively correlated with grassland at the spatial scales of 1 and 2 km and negatively with orchard and vegetable at 0.5 and 1 km. Genetic diversity was best predicted by the grassland + orchard + vegetable model at a spatial scale of 1 km. Based on the method of relative weights, orchard land had the greatest relative importance, followed by grassland and vegetable land, in that order. This study contributes to our understanding of the genetic variation of aphids in agricultural landscapes.     

Fine‐scale Spatial Genetic Structure of Ten Dipterocarp Tree Species in a Bornean Rain Forest

Fine‐scale spatial genetic structure is increasingly recognized as an important factor in the studies of tropical forest trees as it influences genetic diversity of local populations. The biologic mechanisms that generate fine‐scale spatial genetic structure are not fully understood. We studied fine‐scale spatial genetic structure in ten coexisting dipterocarp tree species in a Bornean rain forest using microsatellite markers. Six of the ten species showed statistically significant fine‐scale spatial genetic structure. Fine‐scale spatial genetic structure was stronger at smaller spatial scales (≤ 100 m) than at larger spatial scales (> 100 m) for each species. Multiple regression analysis suggested that seed dispersal distance was important at the smaller spatial scale. At the larger scale (> 100 m) and over the entire sample range (0–1000 m), pollinators and spatial distribution of adult trees were more important determinants of fine‐scale spatial genetic structure. Fine‐scale spatial genetic structure was stronger in species pollinated by less mobile small beetles than in species pollinated by the more mobile giant honeybee (Apis dorsata). It was also stronger in species where adult tree distributions were more clumped. The hypothesized mechanisms underlying the negative correlation between clump size and fine‐scale spatial genetic structure were a large overlap among seed shadows and genetic drift within clumped species.     

Habitat predictors of genetic diversity for two sympatric wetland‐breeding amphibian species

Population genetic diversity is widely accepted as important to the conservation and management of wildlife. However, habitat features may differentially affect evolutionary processes that facilitate population genetic diversity among sympatric species. We measured genetic diversity for two pond‐breeding amphibian species (Dwarf salamanders, Eurycea quadridigitata; and Southern Leopard frogs, Lithobates sphenocephalus) to understand how habitat characteristics and spatial scale affect genetic diversity across a landscape. Samples were collected from wetlands on a longleaf pine reserve in Georgia. We genotyped microsatellite loci for both species to assess population structures and determine which habitat features were most closely associated with observed heterozygosity and rarefied allelic richness. Both species exhibited significant population genetic structure; however, structure in Southern Leopard frogs was driven primarily by one outlier site. Dwarf salamander allelic richness was greater at sites with less surrounding road area within 0.5 km and more wetland area within 1.0 and 2.5 km, and heterozygosity was greater at sites with more wetland area within 0.5 km. In contrast, neither measure of Southern Leopard frog genetic diversity was associated with any habitat features at any scale we evaluated. Genetic diversity in the Dwarf salamander was strongly associated with land cover variables up to 2.5 km away from breeding wetlands, and/or results suggest that minimizing roads in wetland buffers may be beneficial to the maintenance of population genetic diversity. This study suggests that patterns of genetic differentiation and genetic diversity have associations with different habitat features across different spatial scales for two syntopic pond‐breeding amphibian species.     

Rapid divergence of mussel populations despite incomplete barriers to dispersal

Striking genetic structure among marine populations at small spatial scales is becoming evident with extensive molecular studies. Such observations suggest isolation at small scales may play an important role in forming patterns of genetic diversity within species. Isolation‐by‐distance, isolation‐by‐environment and historical priority effects are umbrella terms for a suite of processes that underlie genetic structure, but their relative importance at different spatial and temporal scales remains elusive. Here, we use marine lakes in Indonesia to assess genetic structure and assess the relative roles of the processes in shaping genetic differentiation in populations of a bivalve mussel (Brachidontes sp.). Marine lakes are landlocked waterbodies of similar age (6,000–10,000 years), but with heterogeneous environments and varying degrees of connection to the sea. Using a population genomic approach (double‐digest restriction‐site‐associated DNA sequencing), we show strong genetic structuring across populations (range F_ST: 0.07–0.24) and find limited gene flow through admixture plots. At large spatial scales (>1,400 km), a clear isolation‐by‐distance pattern was detected. At smaller spatial scales (<200 km), this pattern is maintained, but accompanied by an association of genetic divergence with degree of connection. We hypothesize that (incomplete) dispersal barriers can cause initial isolation, allowing priority effects to give the numerical advantage necessary to initiate strong genetic structure. Priority effects may be strengthened by local adaptation, which the data may corroborate by showing a high correlation between mussel genotypes and temperature. Our study indicates an often‐neglected role of (evolution‐mediated) priority effects in shaping population divergence.     

Weak founder effects but significant spatial genetic imprint of recent contraction and expansion of European beech populations

Understanding the ecological and evolutionary processes occurring during species range shifts is important in the current context of global change. Here, we investigate the interplay between recent expansion, gene flow and genetic drift, and their consequences for genetic diversity and structure at landscape and local scales in European beech (Fagus sylvatica L.) On Mont Ventoux, South-Eastern France, we located beech forest refugia at the time of the most recent population minimum, ~150 years ago, and sampled 71 populations (2042 trees) in both refugia and expanding populations over an area of 15,000 ha. We inferred patterns of gene flow and genetic structure using 12 microsatellite markers. We identified six plots as originating from planting, rather than natural establishment, mostly from local genetic material. Comparing genetic diversity and structure in refugia versus recent populations did not support the existence of founder effects: heterozygosity (He = 0.667) and allelic richness (Ar = 4.298) were similar, and F_ST was low (0.031 overall). Still, significant spatial evidence of colonization was detected, with He increasing along the expansion front, while genetic differentiation from the entire pool (β_WT) decreased. Isolation by distance was found in refugia but not in recently expanding populations. Our study indicates that beech capacities for colonization and gene flow were sufficient to preserve genetic diversity despite recent forest contraction and expansion. Because beech has long distance pollen and seed dispersal, these results illustrate a ‘best case scenario’ for the maintenance of high genetic diversity and adaptive potential under climate-change-related range change.Subject terms: Molecular ecology, Forest ecology, Genetic variation     

Genetic structure of Rhododendron ferrugineum at a wide range of spatial scales

Rhododendron ferrugineum L. (Ericaceae) is a subalpine shrub found throughout the Pyrenees and Alps at elevations of 1600-2200 m. We examined relationships between genetic and geographic distance, using 115 dominant amplified fragment length polymorphism (AFLP) markers to assess genetic structure over a wide range of spatial scales. We sampled 17 sites with distances of 4 km to more than 1000 km between them. At these scales we detected no association between geographic distance and genetic distance between populations. This suggests that genetic drift and gene flow are not in equilibrium for these populations. This pattern could have resulted from recent and rapid postglacial colonization, from more recent human disturbance, or as a function of frequent and random "natural" long-distance colonization. At two of our sites we used transects (two horizontal and two vertical with respect to slope at each site) to sample at distances ranging from 10 m to more than 5000 m. At this scale we observed a positive relationship between genetic and spatial distance along two vertical transects, one at each site. We hypothesize that isolation-by-distance at this smaller scale is a function of restricted gene flow via seed dispersal.