Effective population sizes and migration rates in fragmented populations of an endangered insect (Coenagrion mercuriale: Odonata)

1. Effective population sizes (N(e)) and migration rates (m) are critical evolutionary parameters that impact on population survival and determine the relative influence of selection and genetic drift. While the parameter m is well-studied in animal populations, N(e) remains challenging to measure and consequently is only rarely estimated, particularly in insect taxa. 2. We used demographic and genetic methods to estimate N(e) and m in a fragmented population of the endangered damselfly Coenagrion mercuriale to better understand the contrast between genetic and field estimates of these parameters and also to identify the spatial scale over which populations may become locally adapted. 3. We found a contrast between demographic- and genetic-based estimates of these parameters, with the former apparently providing overestimates of N(e), owing to substantial underestimation of the variance in reproductive success, and the latter overestimating m, because spatial genetic structure is weak. 4. The overall N(e) of sites within the population network at Beaulieu Heath, the largest C. mercuriale site in the UK, was estimated to vary between approximately 60 and 2700. 5. While N(e) was not correlated with either the total numbers of adults (N) or the area of habitat, this parameter was always less than N, because of substantial variance in reproductive success. The ratio N(e)/N varied between 0.006 and 0.42 and was generally larger in smaller populations, possibly representing some 'genetic compensation'. 6. From a simple genetic model and these data on N(e) and m, it seems that populations of C. mercuriale have the potential to respond to localized spatial variation in selection and this would need to be considered for future genetic management of this endangered species.     

Population structure and the impact of regional and local habitat isolation upon levels of genetic diversity of the endangered damselfly Coenagrion mercuriale (Odonata: Zygoptera)

1. Coenagrion mercuriale is one of Europe's most threatened damselflies. There is concern for the long‐term persistence of many of its U.K. colonies because adult lifetime movement is limited, making isolated populations susceptible to extinction. 2. Using 14 microsatellite loci we characterised levels of genetic diversity, evidence for a recent decline and the spatial genetic structure for C. mercuriale population in Wales, U.K. 3. Spatial isolation is not an absolute predictor of low genetic diversity at either local or regional scales. 4. One population inhabiting a remote, edge of range site is genetically impoverished with levels of variability (at microsatellite loci) among the lowest reported for any insect species. 5. Agricultural land and high ground are physical barriers to dispersal by adults. 6. Consistent with work from elsewhere, movement by mature C. mercuriale in Pembrokeshire is sufficient to prevent significant genetic differentiation throughout a habitat matrix of some 3–4 km if the suitable habitat sites are <2 km apart and lack barriers to movement. Even within a good habitat matrix, however, genetic isolation by distance develops within 10 km.     

Molecular and ecological evidence for small-scale isolation by distance in an endangered damselfly, Coenagrion mercuriale

Coenagrion mercuriale (Charpentier) (Odonata: Zygoptera) is one of Europe's most threatened damselflies and is listed in the European Habitats directive. We combined an intensive mark-release-recapture (MRR) study with a microsatellite-based genetic analysis for C. mercuriale from the Itchen Valley, UK, as part of an effort to understand the dispersal characteristics of this protected species. MRR data indicate that adult damselflies are highly sedentary, with only a low frequency of interpatch movement that is predominantly to neighbouring sites. This restricted dispersal leads to significant genetic differentiation throughout most of the Itchen Valley, except between areas of continuous habitat, and isolation by distance (IBD), even though the core populations are separated by less than 10 km. An urban area separating some sites had a strong effect on the spatial genetic structure. Average pairwise relatedness between individual damselflies is positive at short distances, reflecting fine-scale genetic clustering and IBD both within- and between-habitat patches. Damselflies from a fragmented habitat have higher average kinship than those from a large continuous population, probably because of poorer dispersal and localized breeding in the former. Although indirect estimates of gene flow must be interpreted with caution, it is encouraging that our results indicate that the spatial pattern of genetic variation matches closely with that expected from direct observations of movement. These data are further discussed with respect to possible barriers to dispersal within the study site and the ecology and conservation of C. mercuriale. To our knowledge, this is the first report of fine-scale genetic structuring in any zygopteran species.     

Geostatistical genetic analysis for inferring the dispersal pattern of a partially clonal species: example of the chestnut blight fungus

Spatial genetic analyses can be used to infer dispersal processes in natural populations. For partially clonal species with alternating sexual and asexual reproduction, the repetition of genotypes must be taken into account in analyses. The methods currently employed to evaluate the relevance of the spatial scale used for the estimation of gene flow are not suitable for these species. We investigated recently developed methods for taking into account repeated genotypes and for determining whether the sampling scale is large enough to capture all the spatial genetic structure existing within a population. We applied these methods to a fungal plant pathogen species, Cryphonectria parasitica, which has caused the death of many American and European chestnut trees since its introduction from Asia at the beginning of the 20th century. These methods were found to be useful for unravelling the effects of clonality and historical gene flow on the spatial genetic structure, and indicated that dispersal processes have probably occurred over a larger spatial scale than previously assumed.     

"Neighbourhood" size, dispersal and density estimates in the prickly forest skink (Gnypetoscincus queenslandiae) using individual genetic and demographic methods

Dispersal, or the amount of dispersion between an individual's birthplace and that of its offspring, is of great importance in population biology, behavioural ecology and conservation, however, obtaining direct estimates from field data on natural populations can be problematic. The prickly forest skink, Gnypetoscincus queenslandiae, is a rainforest endemic skink from the wet tropics of Australia. Because of its log-dwelling habits and lack of definite nesting sites, a demographic estimate of dispersal distance is difficult to obtain. Neighbourhood size, defined as 4piDsigma2 (where D is the population density and sigma2 the mean axial squared parent-offspring dispersal rate), dispersal and density were estimated directly and indirectly for this species using mark-recapture and microsatellite data, respectively, on lizards captured at a local geographical scale of 3 ha. Mark-recapture data gave a dispersal rate of 843 m2/generation (assuming a generation time of 6.5 years), a time-scaled density of 13 635 individuals * generation/km2 and, hence, a neighbourhood size of 144 individuals. A genetic method based on the multilocus (10 loci) microsatellite genotypes of individuals and their geographical location indicated that there is a significant isolation by distance pattern, and gave a neighbourhood size of 69 individuals, with a 95% confidence interval between 48 and 184. This translates into a dispersal rate of 404 m2/generation when using the mark-recapture density estimation, or an estimate of time-scaled population density of 6520 individuals * generation/km2 when using the mark-recapture dispersal rate estimate. The relationship between the two categories of neighbourhood size, dispersal and density estimates and reasons for any disparities are discussed.     

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.     

Plant dispersal, neighbourhood size and isolation by distance

A theoretical relationship between isolation by distance or spatial genetic structure (SGS) and seed and pollen dispersal is tested using extensive spatial-temporal simulations. Although for animals Wright's neighbourhood size N(e) = 4pisigma(2)(t) has been ascertained also, where sigma(2)(t) is the axial variance of distances between parents and offspring, and it was recently confirmed that N(e) = 4pi(sigma(2)(f) + sigma(2)(m))/2 when dispersal of females and males differ, the situation for plants had not been established. This article shows that for a very wide range of conditions, neighbourhood size defined by Crawford's formula N(e) = 4pi(sigma(2)(s) + sigma(2)(p)/2) fully determines SGS, even when dispersal variances of seed (sigma(2)(s)) and pollen sigma(2)(p)) differ strongly. Further, self-fertilization with rate s acts as zero-distance pollen dispersal, and N(e) = 4pi[sigma(2)(s) + sigma(2)(p)(1 - s)/2] fully determines SGS, for most cases where there are both likely parameter values and substantial SGS. Moreover, for most cases, there is a loglinear relationship, I(1) = 0.587 - 0.117 ln(N(e)), between SGS, as measured by I(1), Moran's coefficient for adjacent individuals, and N(e). However, there are several biologically significant exceptions, namely for very low or large N(e), SGS exceeds the loglinear values. There are also important exceptions to Crawford's formula. First, plants with low seed dispersal, high outcross pollen dispersal and high selfing rate show larger SGS than predicted. Second, in plants with very low (near zero) seed dispersal, selfing decreases SGS, opposite expectations. Finally, in some cases seed dispersal is more critical than pollen dispersal, in a manner inconsistent with Crawford's formula.     

Contrasting demographic and genetic estimates of dispersal in the endangered Coahuilan box turtle: a contemporary approach to conservation

The evolutionary viability of an endangered species depends upon gene flow among subpopulations and the degree of habitat patch connectivity. Contrasting population connectivity over ecological and evolutionary timescales may provide novel insight into what maintains genetic diversity within threatened species. We employed this integrative approach to evaluating dispersal in the critically endangered Coahuilan box turtle (Terrapene coahuila) that inhabits isolated wetlands in the desert‐spring ecosystem of Cuatro Ciénegas, Mexico. Recent wetland habitat loss has altered the spatial distribution and connectivity of habitat patches; and we therefore predicted that T. coahuila would exhibit limited movement relative to estimates of historic gene flow. To evaluate contemporary dispersal patterns, we employed mark–recapture techniques at both local (wetland complex) and regional (intercomplex) spatial scales. Gene flow estimates were obtained by surveying genetic variation at nine microsatellite loci in seven subpopulations located across the species’ geographical range. The mark–recapture results at the local spatial scale reveal frequent movement among wetlands that was unaffected by interwetland distance. At the regional spatial scale, dispersal events were relatively less frequent between wetland complexes. The complementary analysis of population genetic substructure indicates strong historic gene flow (global F_ST = 0.01). However, a relationship of genetic isolation by distance across the geographical range suggests that dispersal limitation exists at the regional scale. Our approach of contrasting direct and indirect estimates of dispersal at multiple spatial scales in T. coahuila conveys a sustainable evolutionary trajectory of the species pending preservation of threatened wetland habitats and a range‐wide network of corridors.     

Strong spatial genetic structure in peripheral but not core populations of Sitka spruce [Picea sitchensis (Bong.) Carr.

We examined spatial genetic structure within eight populations of Sitka spruce classified as core or peripheral based on ecological niche, and continuous or disjunct based on species distribution. In each population, 200 trees were spatially mapped and genotyped for eight cDNA-based sequence tagged site (STS) codominant markers. Spatial autocorrelation was assessed by estimating p(ij), the average co-ancestry coefficient, between individuals within distance intervals. The distribution of alleles and genotypes within core populations was almost random, with nonsignificant co-ancestry values among trees as close as 50 m in core populations. In contrast, the distribution of alleles and genotypes within peripheral populations revealed an aggregation of similar multilocus genotypes, with co-ancestry values greater than 0.20 among trees up to 50 m apart and significant, positive values between trees up to 500 m. The relatively high density of reproductive adults in core populations may lead to highly overlapping seed shadows that limit development of spatial genetic structure. However, in peripheral populations with a lower density of adults, the distribution of alleles and genotypes was highly structured, likely due to offspring establishment near maternal trees and subsequent biparental inbreeding, as well as more recent population establishment at the leading edge of post-Pleistocene range expansion. Conserving genetic diversity in peripheral populations may require larger reserves for in situ conservation than required in core populations. These data on spatial genetic structure can be used to provide guidance for sampling strategies for both ex situ conservation and research collections.     

Comparative analysis of the within-population genetic structure in wild cherry (Prunus avium L.) at the self-incompatibility locus and nuclear microsatellites

Gametophytic self-incompatibility (SI) systems in plants exhibit high polymorphism at the SI controlling S-locus because individuals with rare alleles have a higher probability to successfully pollinate other plants than individuals with more frequent alleles. This process, referred to as frequency-dependent selection, is expected to shape number, frequency distribution, and spatial distribution of self-incompatibility alleles in natural populations. We investigated the genetic diversity and the spatial genetic structure within a Prunus avium population at two contrasting gene loci: nuclear microsatellites and the S-locus. The S-locus revealed a higher diversity (15 alleles) than the eight microsatellites (4-12 alleles). Although the frequency distribution of S-alleles differed significantly from the expected equal distribution, the S-locus showed a higher evenness than the microsatellites (Shannon's evenness index for the S-locus: E = 0.91; for the microsatellites: E = 0.48-0.83). Also, highly significant deviations from neutrality were found for the S-locus whereas only minor deviations were found for two of eight microsatellites. A comparison of the frequency distribution of S-alleles in three age-cohorts revealed no significant differences, suggesting that different levels of selection acting on the S-locus or on S-linked sites might also affect the distribution and dynamics of S-alleles. Autocorrelation analysis revealed a weak but significant spatial genetic structure for the multilocus average of the microsatellites and for the S-locus, but could not ascertain differences in the extent of spatial genetic structure between these locus types. An indirect estimate of gene dispersal, which was obtained to explain this spatial genetic pattern, indicated high levels of gene dispersal within our population (sigma(g) = 106 m). This high gene dispersal, which may be partly due to the self-incompatibility system itself, aids the effective gene flow of the microsatellites, thereby decreasing the contrast between the neutral microsatellites and the S-locus.     

Spatial population genetic structure and limited dispersal in a Rocky Mountain alpine stream insect

Using the mitochondrial cytochrome oxidase I (COI) gene, we assessed the phylogeographic structure of Prosimulium neomacropyga, a black fly (Simuliidae) whose distribution in the US Southern Rockies ecoregion is limited to alpine tundra streams. Given high habitat specificity, lack of hydrological connection between streams, and a terrestrial environment restrictive to insect flight, we hypothesized limited gene flow. A spatially nested sampling design showed that grouping populations according to high-elevation 'islands' of alpine tundra (which typically include headwater streams of > 1 watershed) explained a significant proportion of genetic variation while grouping streams according to major watershed (across islands) did not. Nested clade analysis and isolation-by-distance (IBD) relationships further implicated limited ongoing gene flow within but not among the isolated alpine islands. IBD was strong among five streams within an individual island using each of four alternative models of pairwise landscape connectivity for flying insects. Results of all landscape models were positively correlated, suggesting that straight-line distance is an acceptable surrogate for presumably more biologically meaningful connectivity measures in this system. IBD was significantly weaker across the entire study area, comprised of three separate islands. Overall, population structure was significant with F(ST) = 0.38, suggesting limited dispersal across a small spatial extent.     

Genetic diversity and genetic structure of populations of Ranunculus japonicus Thunb. (Ranunculaceae)

Abstract In order to clarify the genetic diversity and population structure of Ranunculus japonicus , allozymic analysis was conducted on 60 populations in southwestern Japan. Considerable genetic variati ons were detected among the populations of R. japonicus . The genetic diversities within species ( H _es = 0.215) and within populations ( H _ep = 0.172) were slightly higher than those of other perennial herbs with widespread distribution and outcrossing plants. Significantly higher values of fixation index were detected in some populations, which might have arisen from restricted mating partners. The majority of genetic variation (approx. 80%) resided within a population and a moderate level of genetic differentiation ( G _ST = 0.203) was observed among populations. The F _ST value (0.203) suggests the existence of a substantial population structure in this species. The highly significant correlation between geographic distance and F _ST values indicates that isolation by distance has played an important role in the construction of the genetic structure of this species.     

Characteristics of sex-biased dispersal and gene flow in coastal river otters: implications for natural recolonization of extirpated populations

River otters (Lontra canadensis) were extirpated from much of their historic distribution because of exposure to pollution and urbanization, resulting in expansive reintroduction programmes that continue today for this and other species of otters worldwide. Bioaccumulation of toxins negatively affects fecundity among mustelids, but high vagility and different dispersal distances between genders may permit otter populations to recover from extirpation caused by localized environmental pollution. Without understanding the influence of factors such as social structure and sex-biased dispersal on genetic variation and gene flow among populations, effects of local extirpation and the potential for natural recolonization (i.e. the need for translocations) cannot be assessed. We studied gene flow among seven study areas for river otters (n = 110 otters) inhabiting marine environments in Prince William Sound, Alaska, USA. Using nine DNA microsatellite markers and assignment tests, we calculated immigration rates and dispersal distances and tested for isolation by distance. In addition, we radiotracked 55 individuals in three areas to determine characteristics of dispersal. Gender differences in sociality and spatial relationships resulted in different dispersal distances. Male river otters had greater gene flow among close populations (within 16-30 km) mostly via breeding dispersal, but both genders exhibited an equal, low probability of natal dispersal; and some females dispersed 60-90 km. These data, obtained in a coastal environment without anthropogenic barriers to dispersal (e.g. habitat fragmentation or urbanization), may serve as baseline data for predicting dispersal under optimal conditions. Our data may indicate that natural recolonization of coastal river otters following local extirpation could be a slow process because of low dispersal among females, and recolonization may be substantially delayed unless viable populations occurred nearby. Because of significant isolation by distance for male otters and low gene flow for females, translocations should be undertaken with caution to help preserve genetic diversity in this species.     

Inverted patterns of genetic diversity in continental and island populations of the heather Erica scoparia s.l.

Aim Using the heather Erica scoparia s.l. as a model, this paper aims to test theoretical predictions that island populations are genetically less diverse than continental ones and to determine the extent to which island and continental populations are connected by pollen‐ and seed‐mediated gene flow. Location Macaronesia, Mediterranean, Atlantic fringe of Europe. Methods Patterns of genetic diversity are described based on variation at two chloroplast DNA (cpDNA) loci and one nuclear DNA (nDNA) locus for 109 accessions across the entire distribution range of the species. Global patterns of genetic differentiation were investigated using principal coordinates analysis. Genetic differentiation between island and continental areas, estimations of pollen‐ and seed‐mediated gene flow, and the presence of phylogeographical signal were assessed by means of F_st /N_ST (continental scale) and F_ij/N_ij (local scale). Extant and past distribution ranges of the species were inferred from niche modelling using layers describing present and Last Glacial Maximum (LGM) macroclimatic conditions. Results The Azores exhibited a significantly higher genetic diversity than the continent. The lowest levels of genetic differentiation were observed between the Azores and the western Mediterranean, and the diversity observed in the Azores resulted from at least two colonization waves. Within the Azores, kinship coefficients showed a significant and much steeper decrease with geographical distance in the cpDNA than in the nDNA. The distribution predicted by LGM models was markedly different from the current potential distribution, particularly in western Europe, where no suitable areas were predicted by LGM models, and along the Atlantic coast of the African continent, where LGM models predicted highly suitable climatic conditions. Main conclusions The higher diversity observed in Azorean than in continental populations is inconsistent with MacArthur and Wilson’s equilibrium model and derived theoretical population genetic expectations. This inverted pattern may be the result of extinction on the continent coupled with multiple island colonization events and subsequent allopatric diversification and lineage hybridization in the Azores. The results highlight the role of allopatric diversification in explaining diversification on islands and suggest that this process has played a much more significant role in shaping Azorean biodiversity than previously thought.     

Maintenance of genetic and morphological identity in two sibling Syrrhopodon species (Calymperaceae,Bryopsida) despite extensive introgression

Bryophytes are a group of land plants in which the role of hybridization has long been challenged. Using genotyping by sequencing to circumvent the lack of molecular variation at selected loci previously used for phylogeny and morphology, we determine the level of genetic and morphological divergence and reproductive isolation between the sibling Syrrhopodon annotinus and S. simmondsii (Calymperaceae, Bryopsida) that occur in sympatry but in different habitats in lowland Amazonian rainforests. A clear morphological differentiation and a low (0.06), but significant Fst derived from the analysis of 183 single nucleotide polymorphisms were observed between the two species. Conspecific pairs of individuals consistently exhibited higher average kinship coefficients along a gradient of geographic isolation than interspecific pairs. The weak, but significant genetic divergence observed is consistent with growing evidence that ecological specialization can lead to genetic differentiation among bryophyte species. Nevertheless, the spatial genetic structures of the two species were significantly correlated, as evidenced by the significant slope of the Mantel test based on kinship coefficients between pairs of interspecific individuals and the geographic distance separating them. Interspecific pairs of individuals are thus more closely related when they are geographically closer, suggesting that isolation‐by‐distance is stronger than the interspecific reproductive barrier and pointing to interspecific gene flow. We conclude that interspecific introgression, whose role has long been questioned in bryophytes, may take place even in species wherein sporophyte production is scarce due to dioicy, raising the question as to what mechanisms maintain differentiation despite weak reproductive isolation.     

Contrasting patterns of mitochondrial and microsatellite population structure in fragmented populations of greater prairie-chickens

Greater prairie-chickens (Tympanuchus cupido pinnatus) were once found throughout the tallgrass prairie of midwestern North America but over the last century these prairies have been lost or fragmented by human land use. As a consequence, many current populations of prairie-chickens have become isolated and small. This fragmentation of populations is expected to lead to reductions in genetic variation as a result of random genetic drift and a decrease in gene flow. As expected, we found that genetic variation at both microsatellite DNA and mitochondrial DNA (mtDNA) markers was reduced in smaller populations, particularly in Wisconsin. There was relatively little range-wide geographical structure (FST) when we examined mtDNA haplotypes but there was a significant positive relationship between genetic (FST) and geographical distance (isolation by distance). In contrast, microsatellite DNA loci revealed significant geographical structure (FST) and a weak effect of isolation by distance throughout the range. These patterns were much stronger when populations with reduced levels of genetic variability (Wisconsin) were removed from the analyses. This suggests that the effects of genetic drift were stronger than gene flow at microsatellite loci, whereas these forces were in range-wide equilibrium at mtDNA markers. These differences between the two molecular markers may be explained by a larger effective population size (Ne) for mtDNA, which is expected in species such as prairie-chickens that have female-biased dispersal and high levels of polygyny. Our results suggest that historic populations of prairie-chickens were once interconnected by gene flow but current populations are now isolated. Thus, maintaining gene flow may be important for the long-term persistence of prairie-chicken populations.     

Population genetic structure and direct observations reveal sex‐reversed patterns of dispersal in a cooperative bird

Sex‐biased dispersal is pervasive and has diverse evolutionary implications, but the fundamental drivers of dispersal sex biases remain unresolved. This is due in part to limited diversity within taxonomic groups in the direction of dispersal sex biases, which leaves hypothesis testing critically dependent upon identifying rare reversals of taxonomic norms. Here, we use a combination of observational and genetic data to demonstrate a rare reversal of the avian sex bias in dispersal in the cooperatively breeding white‐browed sparrow weaver (Plocepasser mahali). Direct observations revealed that (i) natal philopatry was rare, with both sexes typically dispersing locally to breed, and (ii), unusually for birds, males bred at significantly greater distances from their natal group than females. Population genetic analyses confirmed these patterns, as (i) corrected Assignment index (AIc), F_ST tests and isolation‐by‐distance metrics were all indicative of longer dispersal distances among males than females, and (ii) spatial autocorrelation analysis indicated stronger within‐group genetic structure among females than males. Examining the spatial scale of extra‐group mating highlighted that the resulting ‘sperm dispersal’ could have acted in concert with individual dispersal to generate these genetic patterns, but gamete dispersal alone cannot account entirely for the sex differences in genetic structure observed. That leading hypotheses for the evolution of dispersal sex biases cannot readily account for these sex‐reversed patterns of dispersal in white‐browed sparrow weavers highlights the continued need for attention to alternative explanations for this enigmatic phenomenon. We highlight the potential importance of sex differences in the distances over which dispersal opportunities can be detected.