Dispersal and gene flow in the rare, parasitic Large Blue butterfly Maculinea arion

Dispersal is crucial for gene flow and often determines the long‐term stability of meta‐populations, particularly in rare species with specialized life cycles. Such species are often foci of conservation efforts because they suffer disproportionally from degradation and fragmentation of their habitat. However, detailed knowledge of effective gene flow through dispersal is often missing, so that conservation strategies have to be based on mark–recapture observations that are suspected to be poor predictors of long‐distance dispersal. These constraints have been especially severe in the study of butterfly populations, where microsatellite markers have been difficult to develop. We used eight microsatellite markers to analyse genetic population structure of the Large Blue butterfly Maculinea arion in Sweden. During recent decades, this species has become an icon of insect conservation after massive decline throughout Europe and extinction in Britain followed by reintroduction of a seed population from the Swedish island of Öland. We find that populations are highly structured genetically, but that gene flow occurs over distances 15 times longer than the maximum distance recorded from mark–recapture studies, which can only be explained by maximum dispersal distances at least twice as large as previously accepted. However, we also find evidence that gaps between sites with suitable habitat exceeding ～20 km induce genetic erosion that can be detected from bottleneck analyses. Although further work is needed, our results suggest that M. arion can maintain fully functional metapopulations when they consist of optimal habitat patches that are no further apart than ～10 km.     

Limitations to the inference of gene flow at regional geographic scales–an example from the Pieris napi group (Lepidoptera: Pieridae) in Europe

We used hierarchical and pairwise F-statistics to describe genetic differentiation and infer gene flow (M) on local and regional scales within and among parapatric European butterfly taxa in the Pieris napi (L.) group. Within-population allozyme variability is consistently high, and local effective population sizes are inferred to be in the thousands of individuals. The pairwise analysis yields an average neighbourhood area of radius 3.5 km. Among populations within most regions, differentiation is low and M > 2 effective individuals population^-1generation^-1. Pairwise comparisons within the brilannica group show a disjunction indicating that it is out of equilibrium, perhaps as a result of secondary contact between highland and lowland groups. Comparison between meridionalis groups on mainland Italy and Corsica yields M > 12; this is surely loo high and lack of equilibrium resulting from initial colonization is suspected. The hierarchical analysis indicates that 23 ≤0020M≤ 88 among the taxa napi, bryoniae and meridionalis that meet in hybrid zones; no effective gene flow barrier exists among them. This high estimate could also result from recent primary contact, but such a genetic barrier should produce the ‘edge effects' seen in population genetic simulations, and no evidence of this was found among geographically close samples of napi and bryoniae populations from Switzerland. Studies of gene flow among geographic regions are greatly limited by the equilibrium assumption, though studies of local differentiation are much less so. Population studies of gene flow on local scales at regional boundaries provide limited means of testing the equilibrium assumption, and both regional and local analyses provide testable predictions about local population structure. When the equilibrium assumption is not upheld, local patterns at regional boundaries can provide historical information about primary vs. secondary contact.     

Strong gene flow and lack of stable population structure in the face of rapid adaptation to local temperature in a spring-spawning salmonid, the European grayling (Thymallus thymallus)

Gene flow has the potential to both constrain and facilitate adaptation to local environmental conditions. The early stages of population divergence can be unstable because of fluctuating levels of gene flow. Investigating temporal variation in gene flow during the initial stages of population divergence can therefore provide insights to the role of gene flow in adaptive evolution. Since the recent colonization of Lake Lesjaskogsvatnet in Norway by European grayling (Thymallus thymallus), local populations have been established in over 20 tributaries. Multiple founder events appear to have resulted in reduced neutral variation. Nevertheless, there is evidence for local adaptation in early life-history traits to different temperature regimes. In this study, microsatellite data from almost a decade of sampling were assessed to infer population structuring and its temporal stability. Several alternative analyses indicated that spatial variation explained 2-3 times more of the divergence in the system than temporal variation. Over all samples and years, there was a significant correlation between genetic and geographic distance. However, decomposed pairwise regression analysis revealed differing patterns of genetic structure among local populations and indicated that migration outweighs genetic drift in the majority of populations. In addition, isolation by distance was observable in only three of the six years, and signals of population bottlenecks were observed in the majority of samples. Combined, the results suggest that habitat-specific adaptation in this system has preceded the development of consistent population substructuring in the face of high levels of gene flow from divergent environments.     

Land clearing reduces gene flow in the granite outcrop‐dwelling lizard,Ctenophorus ornatus

An important question for the conservation of species dwelling in fragmented habitats is whether changes to the intervening landscape create a barrier to gene flow. Here, we make use of the spatial distribution of the granite outcrop‐dwelling lizard, Ctenophorus ornatus, to compare inferred levels of gene flow between outcrops in a nature reserve with that between outcrops in the adjacent agricultural land. Genetic variation, relatedness and subdivision were compared within groups of individuals from different outcrops similar in size and distance apart at each site. In the agricultural land, we found significantly lower genetic variation within outcrops and greater genetic differentiation between outcrops than in the reserve. Further, the rate at which genetic divergence between outcrops increased over geographical distance was significantly greater in the agricultural land than in the reserve. We also found that individuals were more closely related within outcrops but more distantly related between outcrops in the cleared land. These effects occur over a small spatial scale with an average distance between outcrops of less than five kilometres. Thus, even though land clearing around the outcrops leaves outcrop size unchanged, it restricts gene flow, reducing genetic variation and increasing population structure, with potentially negative consequences for the long‐term persistence of the lizards on these outcrops.     

Comparing direct vs. indirect estimates of gene flow within a population of a scattered tree species

The comparison between historical estimates of gene flow, using variance in allelic frequencies, and contemporary estimates of gene flow, using parentage assignment, is expected to provide insights into ecological and evolutionary processes at work within and among populations. Genetic variation at six microsatellite loci was used to quantify genetic structure in the insect-pollinated, animal-dispersed, low-density tree Sorbus torminalis L. Crantz, and to derive historical estimates of gene flow. The neighbourhood size and root-mean-squared dispersal distance inferred from seedling genotypes ( N _b  = 70 individuals, σ _e  = 417 m) were similar to those inferred from adult genotypes ( N _b  = 114 individuals, σ _e  = 472 m). We also used parentage analyses and a neighbourhood model approach after an evaluation of the statistical properties of this method on simulated data. From our data, we estimated even contributions of seed- and pollen-mediated dispersal to the genetic composition of established seedlings, with both fat-tailed pollen and seed dispersal kernels, and slightly higher mean distance of pollen dispersal (248 m) as compared to seed dispersal (135 m). The resulting contemporary estimate of gene dispersal distance (σ _c  = 211 m) was ∼twofold smaller than the historical estimates. Besides different assumptions and statistical nuances of both approaches, this discrepancy is likely to reflect a recent restriction in the scale of gene flow which requires manager's attention in a context of increasing forest fragmentation.     

Apparent widespread gene flow in the predominantly flightless planthopper Tumidagena minuta

1. To determine whether dispersal biology can predict the pattern of population‐genetic variation among insect populations accurately, allozyme variation was assayed for populations of a saltmarsh planthopper, Tumidagena minuta, in which > 99% of the adults are flightless. 2. The pattern of genetic isolation by distance in T. minuta was compared with that in other insects, to determine whether it was similar to isolation by distance in other sedentary insects. 3. In contrast to predictions, the pattern of isolation by distance in T. minuta was most similar to that seen in the most mobile insects in a recent review of population‐genetic variation in insects. Furthermore, population‐genetic subdivision over a spatial scale of > 400 km was weak. 4. Possible causes of the apparent contradiction between dispersal biology and population‐genetic structure in this species are discussed. The results for T. minuta highlight the fact that although mobility is generally correlated with gene flow in insects, studies of population‐genetic variation must be combined with direct studies of dispersal to understand fully the degree to which populations exchange individuals.     

A simulation-based evaluation of methods for inferring linear barriers to gene flow

Different analytical techniques used on the same data set may lead to different conclusions about the existence and strength of genetic structure. Therefore, reliable interpretation of the results from different methods depends on the efficacy and reliability of different statistical methods. In this paper, we evaluated the performance of multiple analytical methods to detect the presence of a linear barrier dividing populations. We were specifically interested in determining if simulation conditions, such as dispersal ability and genetic equilibrium, affect the power of different analytical methods for detecting barriers. We evaluated two boundary detection methods (Monmonier's algorithm and WOMBLING), two spatial Bayesian clustering methods (TESS and GENELAND), an aspatial clustering approach (STRUCTURE), and two recently developed, non-Bayesian clustering methods [PSMIX and discriminant analysis of principal components (DAPC)]. We found that clustering methods had higher success rates than boundary detection methods and also detected the barrier more quickly. All methods detected the barrier more quickly when dispersal was long distance in comparison to short-distance dispersal scenarios. Bayesian clustering methods performed best overall, both in terms of highest success rates and lowest time to barrier detection, with GENELAND showing the highest power. None of the methods suggested a continuous linear barrier when the data were generated under an isolation-by-distance (IBD) model. However, the clustering methods had higher potential for leading to incorrect barrier inferences under IBD unless strict criteria for successful barrier detection were implemented. Based on our findings and those of previous simulation studies, we discuss the utility of different methods for detecting linear barriers to gene flow.     

Isolation by distance and vicariance drive genetic structure of a coral reef fish in the Pacific Ocean

We studied the genetic diversity of a coral reef fish species to investigate the origin of the differentiation. A total of 727 Acanthurus triostegus collected from 15 locations throughout the Pacific were analyzed for 20 polymorphic loci. The genetic structure showed limited internal disequilibrium within each population; 3.7% of the loci showed significant Hardy-Weinberg disequilibrium, mostly associated with Adh*, and we subsequently removed this locus from further analysis of geographic pattern. The genetic structure of A. triostegus throughout the tropical Pacific Ocean revealed a strong geographic pattern. Overall, there was significant population differentiation (multilocus F(ST) = 0.199), which was geographically structured according to bootstraps of neighbor-joining analysis on Nei's unbiased genetic distances and AMOVA analysis. The genetic structure revealed five geographic groups in the Pacific Ocean: western Pacific (Guam, Philippines, Palau, and Great Barrier Reef); central Pacific (Solomons, New Caledonia, and Fiji); and three groups made up of the eastern populations, namely Hawaiian Archipelago (north), Marquesas (equatorial), and southern French Polynesia (south) that incorporates Clipperton Island located in the northeastern Pacific. In addition, heterozygosity values were found to be geographically structured with higher values grouped within Polynesian and Clipperton populations, which exhibited lower population size. Finally, the genetic differentiation (F(ST)) was significantly correlated with geographic distance when populations from the Hawaiian and Marquesas archipelagos were separated from all the other locations. These results show that patterns of differentiation vary within the same species according to the spatial scale, with one group probably issued from vicariance, whereas the other followed a pattern of isolation by distance. The geographic pattern for A. triostegus emphasizes the diversity of the evolutionary processes that lead to the present genetic structure with some being more influential in certain areas or according to a particular spatial scale.     

Urban landscape genetics: canopy cover predicts gene flow between white-footed mouse (Peromyscus leucopus) populations in New York City

In this study, I examine the influence of urban canopy cover on gene flow between 15 white-footed mouse (Peromyscus leucopus) populations in New York City parklands. Parks in the urban core are often highly fragmented, leading to rapid genetic differentiation of relatively nonvagile species. However, a diverse array of 'green' spaces may provide dispersal corridors through 'grey' urban infrastructure. I identify urban landscape features that promote genetic connectivity in an urban environment and compare the success of two different landscape connectivity approaches at explaining gene flow. Gene flow was associated with 'effective distances' between populations that were calculated based on per cent tree canopy cover using two different approaches: (i) isolation by effective distance (IED) that calculates the single best pathway to minimize passage through high-resistance (i.e. low canopy cover) areas, and (ii) isolation by resistance (IBR), an implementation of circuit theory that identifies all low-resistance paths through the landscape. IBR, but not IED, models were significantly associated with three measures of gene flow (Nm from F(ST) , BayesAss+ and Migrate-n) after factoring out the influence of isolation by distance using partial Mantel tests. Predicted corridors for gene flow between city parks were largely narrow, linear parklands or vegetated spaces that are not managed for wildlife, such as cemeteries and roadway medians. These results have implications for understanding the impacts of urbanization trends on native wildlife, as well as for urban reforestation efforts that aim to improve urban ecosystem processes.     

Microsatellite analysis of the phylogeography, Pleistocene history and secondary contact hypotheses for the killifish, Fundulus heteroclitus

The mummichog, Fundulus heteroclitus, exhibits extensive latitudinal clinal variation in a number of physiological and biochemical traits, coupled with phylogeographical patterns at mitochondrial and nuclear DNA loci that suggest a complicated history of spatially variable selection and secondary intergradation. This species continues to serve as a model for understanding local and regional adaptation to variable environments. Resolving the influences of historical processes on the distribution of genetic variation within and among extant populations of F. heteroclitus is crucial to a better understanding of how populations evolve in the context of contemporary environments. In this study, we analysed geographical patterns of genetic variation at eight microsatellite loci among 15 populations of F. heteroclitus distributed throughout the North American range of the species from Nova Scotia to Georgia. Genetic variation in Northern populations was lower than in Southern populations and was strongly correlated with latitude throughout the species range. The most common Northern alleles at all eight loci exhibited concordant latitudinal clinal patterns, and the existence of an abrupt transition zone in allele frequencies between Northern and Southern populations was similar to that observed for mitochondrial DNA and allozyme loci. A significant pattern of isolation by distance was observed both within and between northern and southern regions. This pattern was unexpected, particularly for northern populations, given the recent colonization history of post-Pleistocene habitats, and was inconsistent with either a recent northward population expansion or a geographically restricted northern Pleistocene refugium. The data provided no evidence for recent population bottlenecks, and estimates of historical effective population sizes suggest that post-Pleistocene populations have been large throughout the species distribution. These results suggest that F. heteroclitus was broadly distributed throughout most of its current range during the last glacial event and that the abrupt transition in allele frequencies that separate Northern and Southern populations may reflect regional disequilibrium conditions associated with the post-Pleistocene colonization history of habitats in that region.     

Microgeographic genetic structure and gene flow in Hibiscus moscheutos (Malvaceae) populations

Microgeographic genetic variation in populations of a wetland macrophyte, Hibiscus moscheutos L. (Malvaceae), was investigated using allozyme polymorphism. The species is a self-compatible insect-pollinated perennial, and seeds are water dispersed (hydrochory). Six hundred plants were analyzed from eight brackish and two freshwater populations within the Rhode River watershed/estuarine system. The genetic structure of the populations was assessed by fixation indices and spatial autocorrelation analyses. The degree of genetic differentiation among sites and gene flow between all paired combinations of sites (M ) was analyzed using three hypothetical gene flow models. Fixation indices indicated almost complete panmixia within populations, and spatial autocorrelations showed that genotypes were randomly distributed within sites, most likely the result of water dispersal of seeds. Allele frequencies were significantly different among sites, and estimated FST indicated moderate genetic differentiation (_ = 0.062). Genetic differences between populations were mostly explained by a gene flow model that accounted for the location of populations relative to the tidal stream. The importance of hydrochory in affecting spatial genetic structure was thus suggested both within and among H. moscheutos populations.     

Morphological and genetic differentiation in anadromous smelt Osmerus mordax (Mitchill): disentangling the effects of geography and morphology on gene flow

Morphological analyses were combined with genetic analyses at nine microsatellite loci to examine the determinants of gene flow at 21 spawning locations of rainbow smelt Osmerus mordax along the east coast of Canada. Associations between morphology, geography and gene flow were examined using a computational geometric approach and partial Mantel tests. Significant barriers to gene flow and discontinuities in morphology were observed between Newfoundland and mainland Canada, as well as within Newfoundland samples. On regional scales, contrasting patterns were present with restricted gene flow between Newfoundland populations ( F _ST= c . 0·11) and high gene flow between mainland populations ( F _ST= c . 0·017). Within Newfoundland populations, geographic distance was significantly associated with gene flow ( r = 0·85, P < 0·001) contrasting mainland samples where gene flow was most associated with phenotypic divergence ( r = 0·33, P < 0·001). At large spatial scales, weak ( r = 0·19, P = 0·02) associations between gene flow and geographic distance were observed, and moderate associations were also observed between gene flow and morphology ( r = 0·28, P < 0·001). The presence of significant genetic isolation by distance in Newfoundland samples and the clear discontinuity associated with the Cabot Strait suggest geography may be the primary determinant of gene flow. Interestingly, the association between genetic and morphological divergence within mainland samples and overall, supports the hypothesis that gene flow may be moderated by morphological divergence at larger spatial scales even in high gene flow environments.     

Small-scale spatial genetic structure in the Central African rainforest tree species Aucoumea klaineana: a stepwise approach to infer the impact of limited gene dispersal, population history and habitat fragmentation

Under the isolation-by-distance model, the strength of spatial genetic structure (SGS) depends on seed and pollen dispersal and genetic drift, which in turn depends on local demographic structure. SGS can also be influenced by historical events such as admixture of differentiated gene pools. We analysed the fine-scale SGS in six populations of a pioneer tree species endemic to Central Africa, Aucoumea klaineana. To infer the impacts of limited gene dispersal, population history and habitat fragmentation on isolation by distance, we followed a stepwise approach consisting of a Bayesian clustering method to detect differentiated gene pools followed by the analysis of kinship-distance curves. Interestingly, despite considerable variation in density, the five populations situated under continuous forest cover displayed very similar extent of SGS. This is likely due to an increase in dispersal distance with decreased tree density. Admixture between two gene pools was detected in one of these five populations creating a distinctive pattern of SGS. In the last population sampled in open habitat, the genetic diversity was in the same range as in the other populations despite a recent habitat fragmentation. This result may due to the increase of gene dispersal compensating the effect of the disturbance as suggested by the reduced extent of SGS estimated in this population. Thus, in A. klaineana, the balance between drift and dispersal may facilitate the maintenance of genetic diversity. Finally, from the strength of the SGS and population density, an indirect estimate of gene dispersal distances was obtained for one site: the quadratic mean parent-offspring distance, sigma(g), ranged between 210 m and 570 m.     

Population structure and genetic variability of pearl oyster Pinctada mazatlanica along Pacific coasts from Mexico to Panama

Populations of the Calafia pearl oyster Pinctada mazatlanica ofthe American Pacific coasts have been considered endangeredbecause of overfishing and/or alteration to coastal areas. Weassessed genetic variability and the pattern of populationstructure among 9 samples collected from Mexico to Panama, usingmtDNA RFLP analysis of two genes: 12S rRNA and subunit one ofCytochrome oxydase (COI). Haplotype diversity varied from 0.000to 0.856. The Panama population appeared to be monomorphic, whilethe other samples exhibited a level of haplotypic variabilitysimilar to those reported in the literature for the same kind ofanalysis on other bivalves species. A test for the impact ofdemographic history on genetic diversity was applied on thesequence data, and the results were congruent with a recentdecline of population sizes. Genetic differentiation was shown tofollow a scheme of isolation by distance, with low levels ofdifferentiation at the scales of ten to one hundred kilometres,whereas stronger and significant genetic structure was detectedat a larger scale. Three significantly distinct groups could thenbe defined, which correspond to Northern Mexico, Southern Mexico,and Panama.     

Genetic consequences of natal dispersal in the colonial lesser kestrel

Dispersal is a life-history trait that plays a fundamental role in population dynamics, influencing evolution, species distribution, and the genetics and structure of populations. In spite of the fact that dispersal has been hypothesized to be an efficient behavioural mechanism to avoid inbreeding, the expected relationship between dispersal and mate relatedness still remains controversial. Here, we examine the genetic consequences of natal dispersal, namely the higher chance of obtaining genetically less similar mates as a result of moving from natal to breeding sites, in a lesser kestrel (Falco naumanni) population. Relatedness between individuals tended to decrease with distance between their breeding colonies, indicating that the study population follows an 'isolation-by-distance' pattern of spatial genetic structure. Such a fine-scale genetic structure generates a scenario in which individuals can potentially increase the chance of obtaining genetically less similar mates by dispersing over larger distances from their natal colony. Using dispersal information and genotypic data, we showed that mate relatedness decreased with natal dispersal distance, an effect that remained significant both while including and excluding philopatric individuals from the data set. These results, together with the well known detrimental consequences of reduced genetic diversity in the study population, suggest that dispersal may have evolved, at least in part, to avoid the negative fitness consequences of mating with genetically similar individuals.     

Montane barriers to gene flow of Melanoplus bivittatus (Orthoptera: Acrididae)*

To determine the effectiveness of montane barriers to gene flow of the two-striped grasshopper, Melanoplus bivittatus, collections were made at four sites at the corners of a square location, each bisected by one of five Colorado mountains (total: five locations, 18 sites). The premise of this design is that samples from opposite sides of an effective barrier should have more genotypic differences than samples from the same side. Genotypic differences were measured for each grasshopper as 18 morphological measurements and six enzyme loci. Dissimilarities in morphology or enzyme loci were expressed as significant estimates of Mahalanobis squared distance (respectively, D² or G²) between two sample sites. The presence of foothills had no influence on the number of significant G² (age and sex classes combined) or significant D² (age and sex classes separated). For samples near the foothills, adult females had more significant D² than adult males. For samples near the two higher mountains, adult males had more significant D² for site samples from opposite sides. A combination of gene flow and selection is the most plausible explanation for most of these results.