Genetic analysis reveals human-mediated long-distance dispersal among war cemeteries in Trifolium micranthum

Many grassland plant species have limited capacity to disperse their seeds beyond local boundaries by natural means. Meanwhile, various forms of human transportation are observed to provide long-distance dispersal. However, the contribution of human-mediated dispersal to the spatial dynamics of established (meta)populations on the regional scale has only scarcely been addressed. Trifolium micranthum is a very small legume that is rare in western Belgium. It is found (i) inland, as a lawn weed on war cemeteries that are managed by the Commonwealth War Graves Commission, and (ii) in coastal grasslands that are managed for conservation. These stations lie over 30?km apart, though interestingly, a number of satellite war cemeteries are found within the vicinity of the coastal reserves, and these also harbor T. micranthum. We have sampled plants from these populations to study their affinities through genetic analysis (AFLP). Results indicated that the coastal cemetery populations were closely related to the inland cemeteries and not to the coastal reserves. This most likely is the result of exchange of seeds with composted lawn material among regions, whereas mowing machinery and livestock may effectuate frequent dispersal within regions. Assignment tests nonetheless indicated limited genetic admixture to have taken place at the coastal cemeteries. Thus, human management is responsible for the realization of gene flow among these remote gene pools. Human-mediated long-distance dispersal is likely to play a role in the local or regional dynamics of many other native plants, including more common and well-dispersing species.     

Population genetic structure and long-distance dispersal among seabird populations: implications for colony persistence

Dramatic local population decline brought about by anthropogenic-driven change is an increasingly common threat to biodiversity. Seabird life history traits make them particularly vulnerable to such change; therefore, understanding population connectivity and dispersal dynamics is vital for successful management. Our study used a 357-base pair mitochondrial control region locus sequenced for 103 individuals and 18 nuclear microsatellite loci genotyped for 245 individuals to investigate population structure in the Atlantic and Pacific populations of the pelagic seabird, Leach's storm-petrel Oceanodroma leucorhoa leucorhoa. This species is under intense predation pressure at one regionally important colony on St Kilda, Scotland, where a disparity between population decline and predation rates hints at immigration from other large colonies. AMOVA, F(ST), Φ(ST) and Bayesian cluster analyses revealed no genetic structure among Atlantic colonies (Global Φ(ST) = -0.02 P > 0.05, Global F(ST) = 0.003, P > 0.05, STRUCTURE K = 1), consistent with either contemporary gene flow or strong historical association within the ocean basin. The Pacific and Atlantic populations are genetically distinct (Global Φ(ST) = 0.32 P < 0.0001, Global F(ST) = 0.04, P < 0.0001, STRUCTURE K = 2), but evidence for interocean exchange was found with individual exclusion/assignment and population coalescent analyses. These findings highlight the importance of conserving multiple colonies at a number of different sites and suggest that management of this seabird may be best viewed at an oceanic scale. Moreover, our study provides an illustration of how long-distance movement may ameliorate the potentially deleterious impacts of localized environmental change, although direct measures of dispersal are still required to better understand this process.     

Diverse spore rains and limited local exchange shape fern genetic diversity in a recently created habitat colonized by long-distance dispersal

Background and Aims

Populations established by long-distance colonization are expected to show low levels of genetic variation per population, but strong genetic differentiation among populations. Whether isolated populations indeed show this genetic signature of isolation depends on the amount and diversity of diaspores arriving by long-distance dispersal, and time since colonization. For ferns, however, reliable estimates of long-distance dispersal rates remain largely unknown, and previous studies on fern population genetics often sampled older or non-isolated populations. Young populations in recent, disjunct habitats form a useful study system to improve our understanding of the genetic impact of long-distance dispersal.

Methods

Microsatellite markers were used to analyse the amount and distribution of genetic diversity in young populations of four widespread calcicole ferns (Asplenium scolopendrium, diploid; Asplenium trichomanes subsp. quadrivalens, tetraploid; Polystichum setiferum, diploid; and Polystichum aculeatum, tetraploid), which are rare in The Netherlands but established multiple populations in a forest (the Kuinderbos) on recently reclaimed Dutch polder land following long-distance dispersal. Reference samples from populations throughout Europe were used to assess how much of the existing variation was already present in the Kuinderbos.

Key Results

A large part of the Dutch and European genetic diversity in all four species was already found in the Kuinderbos. This diversity was strongly partitioned among populations. Most populations showed low genetic variation and high inbreeding coefficients, and were assigned to single, unique gene pools in cluster analyses. Evidence for interpopulational gene flow was low, except for the most abundant species.

Conclusions

The results show that all four species, diploids as well as polyploids, were capable of frequent long-distance colonization via single-spore establishment. This indicates that even isolated habitats receive dense and diverse spore rains, including genotypes capable of self-fertilization. Limited gene flow may conserve the genetic signature of multiple long-distance colonization events for several decades.     

An investigation into effects of long-distance seed dispersal on organelle population genetic structure and colonization rate: a model analysis

A simulation-based modelling approach is used to examine the effects of stratified seed dispersal (representing the distribution of the majority of dispersal around the maternal parent and also rare long-distance dispersal) on the genetic structure of maternally inherited genomes and the colonization rate of expanding plant populations. The model is parameterized to approximate postglacial oak colonization in the UK, but is relevant to plant populations that exhibit stratified seed dispersal. The modelling approach considers the colonization of individual plants over a large area (three 500 km x 10 km rolled transects are used to approximate a 500 km x 300 km area). Our approach shows how the interaction of plant population dynamics with stratified dispersal can result in a spatially patchy haplotype structure. We show that while both colonization speeds and the resulting genetic structure are influenced by the characteristics of the dispersal kernel, they are robust to changes in the periodicity of long-distance events, provided the average number of long-distance dispersal events remains constant. We also consider the effects of additional physical and environmental mechanisms on plant colonization. Results show significant changes in genetic structure when the initial colonization of different haplotypes is staggered over time and when a barrier to colonization is introduced. Environmental influences on survivorship and fecundity affect both the genetic structure and the speed of colonization. The importance of these mechanisms in relation to the postglacial spread and genetic structure of oak in the UK is discussed.     

High sexual reproduction and limited contemporary dispersal in the ectomycorrhizal fungus Tricholoma scalpturatum: new insights from population genetics and spatial autocorrelation analysis

Dispersal and establishment are fundamental processes influencing the response of species to environmental changes, and the long-term persistence of populations. A previous study on the symbiotic ectomycorrhizal fungus Tricholoma scalpturatum revealed strong genetic differentiations between populations in Western Europe, suggesting restricted dispersal for this wind-dispersed cosmopolitan fungus. Two distinct genetic groups (genetic groups 1 and 2), co-occurring in some locations, were also identified and could correspond to cryptic species. In the present work, we examine the reproductive strategy and dispersal biology of the two T. scalpturatum's genetic groups. Variable molecular markers (intersimple sequence repeats and intergenic spacer 2-restriction fragment length polymorphisms) and spatial autocorrelation analyses were used to examine fine-scale patterns (< 140 m) of genetic structure, in an effort to determine the physical scale at which genetic structure exists. A total of 473 fruit bodies were mapped and collected over 3 years from two plots located in the south of France, including 219 and 254 samples from group 1 and group 2, respectively. High genetic diversity and the presence of numerous small genets were observed in both groups. Autocorrelation analyses revealed significant positive spatial genetic structures of genets at close distances (up to few metres for both groups). Mantel tests confirmed this isolation-by-distance pattern. These results clearly demonstrate high sexual reproduction and spatial structuring of genets at very small geographical scales in this wind-dispersed ectomycorrhizal fungal species, a pattern consistent with restricted contemporary dispersal of spores.     

Low genetic variability, female-biased dispersal and high movement rates in an urban population of Eurasian badgers Meles meles

1. Urban and rural populations of animals can differ in their behaviour, both in order to meet their ecological requirements and due to the constraints imposed by different environments. The study of urban populations can therefore offer useful insights into the behavioural flexibility of a species as a whole, as well as indicating how the species in question adapts to a specifically urban environment. 2. The genetic structure of a population can provide information about social structure and movement patterns that is difficult to obtain by other means. Using non-invasively collected hair samples, we estimated the population size of Eurasian badgers Meles meles in the city of Brighton, England, and calculated population-specific parameters of genetic variability and sex-specific rates of outbreeding and dispersal. 3. Population density was high in the context of badger densities reported throughout their range. This was due to a high density of social groups rather than large numbers of individuals per group. 4. The allelic richness of the population was low compared with other British populations. However, the rate of extra-group paternity and the relatively frequent (mainly temporary) intergroup movements suggest that, on a local scale, the population was outbred. Although members of both sexes visited other groups, there was a trend for more females to make intergroup movements. 5. The results reveal that urban badgers can achieve high densities and suggest that while some population parameters are similar between urban and rural populations, the frequency of intergroup movements is higher among urban badgers. In a wider context, these results demonstrate the ability of non-invasive genetic sampling to provide information about the population density, social structure and behaviour of urban wildlife.     

Pollen dispersal and genetic variation in an early‐successional forest herb in a peri‐urban forest

In order to conserve forest plant species under the particularly high constraints that represent urban surroundings, it is necessary to identify the key factors for population persistence. This study examined within‐ and between‐population pollen dispersal using fluorescent dye as pollen analogue, and genetic variation and structure using 15 allozyme loci in Centaurium erythraea, an insect‐pollinated, early‐successional forest biennial herb occurring in a peri‐urban forest (Brussels urban zone, Belgium). Dye dispersal showed an exponential decay distribution, with most dye transfers occurring at short distances (<15 m), and only a few long‐distance events (up to 743 m). Flowers of C. erythraea are mainly visited by Syrphids (Diptera) and small bees, which are usually considered as short‐distance pollen dispersers, and occasionally by bumblebees, which are usually longer‐distance pollen dispersers. Small and large dye source populations differed in dye deposition patterns. The populations showed low genetic diversity, high inbreeding coefficients (F_IS) and high genetic differentiation (F_ST), suggesting restricted gene flow, which can be expected for an early‐successional biennial species with a predominantly selfing breeding system and fluctuating population sizes. The positive relationship between recruitment rate and allelic richness and expected heterozygosity, and the absence of significant correlations between genetic variation and population size suggest seedling recruitment from the seed bank, contributing to maintain genetic diversity. Long‐distance dye dispersal events indicate pollinator movements along urban forest path and road verges. These landscape elements might therefore have a potential conservation value by contributing to connectivity of early‐successional species populations located in patchy open habitats.     

QTL and quantitative genetic analysis of beak morphology reveals patterns of standing genetic variation in an Estrildid finch

The intra- and interspecific diversity of avian beak morphologies is one of the most compelling examples for the power of natural selection acting on a morphological trait. The development and diversification of the beak have also become a textbook example for evolutionary developmental biology, and variation in expression levels of several genes is known to causally affect beak shape. However, until now, no genomic polymorphisms have been identified, which are related to beak morphology in birds. QTL mapping does reveal the location of causal polymorphisms, albeit with poor spatial resolution. Here, we estimate heritability and genetic correlations for beak length, depth and width and perform a QTL linkage analysis for these traits based on 1404 informative single-nucleotide polymorphisms genotyped in a four-generation pedigree of 992 captive zebra finches (Taeniopygia guttata). Beak size, relative to body size, was sexually dimorphic (larger in males). Heritability estimates ranged from 0.47 for beak length to 0.74 for beak width. QTL mapping revealed four to five regions of significant or suggestive genome-wide linkage for each of the three beak dimensions (nine different regions in total). Eight out of 11 genes known to influence beak morphology are located in these nine peak regions. Five QTL do not cover known candidates demonstrating that yet unknown genes or regulatory elements may influence beak morphology in the zebra finch.     

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.     

Spatial genetic structure of the mountain pine beetle (Dendroctonus ponderosae) outbreak in western Canada: historical patterns and contemporary dispersal

Environmental change has a wide range of ecological consequences, including species extinction and range expansion. Many studies have shown that insect species respond rapidly to climatic change. A mountain pine beetle epidemic of record size in North America has led to unprecedented mortality of lodgepole pine, and a significant range expansion to the northeast of its historic range. Our goal was to determine the spatial genetic variation found among outbreak population from which genetic structure, and dispersal patterns may be inferred. Beetles from 49 sampling locations throughout the outbreak area in western Canada were analysed at 13 microsatellite loci. We found significant north-south population structure as evidenced by: (i) Bayesian-based analyses, (ii) north-south genetic relationships and diversity gradients; and (iii) a lack of isolation-by-distance in the northernmost cluster. The north-south structure is proposed to have arisen from the processes of postglacial colonization as well as recent climate-driven changes in population dynamics. Our data support the hypothesis of multiple sources of origin for the outbreak and point to the need for population specific information to improve our understanding and management of outbreaks. The recent range expansion across the Rocky Mountains into the jack/lodgepole hybrid and pure jack pine zones of northern Alberta is consistent with a northern British Columbia origin. We detected no loss of genetic variability in these populations, indicating that the evolutionary potential of mountain pine beetle to adapt has not been reduced by founder events. This study illustrates a rapid range-wide response to the removal of climatic constraints, and the potential for range expansion of a regional population.     

From global to local genetic structuring in the red gorgonian Paramuricea clavata: the interplay between oceanographic conditions and limited larval dispersal

Defining the scale of connectivity among marine populations and identifying the barriers to gene flow are tasks of fundamental importance for understanding the genetic structure of populations and for the design of marine reserves. Here, we investigated the population genetic structure at three spatial scales of the red gorgonian Paramuricea clavata (Cnidaria, Octocorallia), a key species dwelling in the coralligenous assemblages of the Mediterranean Sea. Colonies of P. clavata were collected from 39 locations across the Mediterranean Sea from Morocco to Turkey and analysed using microsatellite loci. Within three regions (Medes, Marseille and North Corsica), sampling was obtained from multiple locations and at different depths. Three different approaches (measures of genetic differentiation, Bayesian clustering and spatially explicit maximum‐difference algorithm) were used to determine the pattern of genetic structure. We identified genetic breaks in the spatial distribution of genetic diversity, which were concordant with oceanographic conditions in the Mediterranean Sea. We revealed a high level of genetic differentiation among populations and a pattern of isolation by distance across the studied area and within the three regions, underlining short effective larval dispersal in this species. We observed genetic differentiation among populations in the same locality dwelling at different depths, which may be explained by local oceanographic conditions and which may allow a process of local adaptation of the populations to their environment. We discuss the implications of our results for the conservation of the species, which is exposed to various threats.     

A long-term genetic survey of an ungulate population reveals balancing selection acting on MHC through spatial and temporal fluctuations in selection

We explored a 13-year genetic survey of the major histocompatibility complex (MHC) and neutral loci of the Soay sheep population of St Kilda to test the existence and causes of balancing selection at the MHC. The sheep population experiences demographic fluctuations, partly driven by the nematode Teladorsagia circumcincta. The spatial differentiation detected at the MHC was comparable to that at neutral loci between 1988 and 1996, but significantly lower between 1996 and 2000. The rate of temporal genetic differentiation was higher at the MHC, but within the Eastern heft only. These comparisons of spatial and temporal divergence at MHC and non-MHC loci provide strong evidence of balancing selection at the MHC, acting through spatial and temporal heterogeneity in selection pressure. This heterogeneity could be due to fluctuations in the selection imposed by parasites, either directly, because the prevalence in T. circumcincta varies in space and time, or indirectly, because the fitness costs of infection may vary with resource availability.     

The NATURA 2000 database as a tool in the analysis of habitat selection at large scales: factors affecting the occurrence of pine and stone martens in Southern Europe

Most studies on habitat selection among animals are conducted at local scales, whereas reliable determination of species requirements at larger spatial scales can be problematic. We used data available for NATURA 2000 sites to determine the habitat requirements of two relatively widespread and common species—pine marten Martes martes and stone marten Martes foina—in Southern Europe. Using presence-absence data, we applied statistical models at two spatial scales. At the local scale (within the dispersal distances of the species), sites occupied by martens were compared with unoccupied sites using buffers of different sizes, whereas at the regional scale, unoccupied sites were selected randomly. To adjust for spatial autocorrelation of data, penalised quasi-likelihood approximations were used. Both species of martens demonstrated preferences for areas with higher proportions of forest cover and steeper terrain. At the local, but not at the regional, scale, pine martens occurred at lower elevations, whilst stone martens occurred at higher elevations. We found that climatic variables (mean temperature, precipitation) had no significant effect on the presence of the martens. The results of our analyses generally confirmed findings of previous studies on habitat selection of both marten species in Southern Europe. This demonstrates the utility of data collected for NATURA 2000 sites for use in various analyses such as conservation planning and evaluating the impact of climate change on the distribution of animal species.     

Estimating population boundaries using regional and local-scale spatial genetic structure: an example in Eucalyptus globulus

Eucalyptus globulus Labill is a foundation tree species over its disjunct distribution in southeastern Australia. The quality of its pulp makes it the most important hardwood species in the world. The importance of E. globulus prompted the establishment of common gardens from seed collected across its geographic range. This enabled us to study the genetic structure of the species, its population boundaries, and gene flow using 444 trees from different open-pollinated families that were genotyped at 16 microsatellite loci. A Bayesian clustering method was used to resolve five genetically distinct groups across the geographical range. These groups were identified as regions, which varied in diameter from 38 to 294?km and contain 4 to 16 putative populations. For two of these regional groups, we used spatial autocorrelation analysis based on assignment of trees to their natural stands to examine gene flow within each region. Consistent significant local-scale spatial structure occurred in both regions. Pairs of individuals within a region showed significant genetic similarity that extended beyond 40?km, suggesting distant movement of pollen. This suggests that breeding populations in E. globulus are much bigger than traditionally accepted in eucalypts. Our results are important for the management of genetic diversity and breeding populations in E. globulus. Similar studies of a variety of eucalypts pollinated by insects and birds will determine whether the local-scale genetic structure of E. globulus is unusual.     

Landscape level analysis of mountain pine beetle in British Columbia, Canada: spatiotemporal development and spatial synchrony within the present outbreak

Eruptive herbivores can exert profound landscape level influences. For example, the ongoing mountain pine beetle outbreak in British Columbia, Canada, has resulted in mortality of mature lodgepole pine over >7 million ha. Analysis of the spatio‐temporal pattern of spread can lend insights into the processes initiating and/or sustaining such phenomena. We present a landscape level analysis of the development of the current outbreak. Aerial survey assessments of tree mortality, projected onto discrete 12×12 km cells, were used as a proxy for insect population density. We examined whether the outbreak potentially originated from an epicenter and spread, or whether multiple localized populations erupted simultaneously at spatially disjunct locations. An aspatial cluster analysis of time series from 1990 to 2003 revealed four distinct time series patterns. Each time series demonstrated a general progression of increasing mountain pine beetle populations. Plotting the geographical locations of each temporal pattern revealed that the outbreak occurred first in an area of west‐central British Columbia, and then in an area to the east. The plot further revealed many localized infestations erupted in geographically disjunct areas, especially in the southern portion of the province. Autologistic regression analyses indicated a significant, positive association between areas where the outbreak first occurred and conservation lands. For example, the delineated area of west‐central British Columbia is comprised of three conservation parks and adjacent working forest. We further examined how population synchrony declines with distance at different population levels. Examination of the spatial dependence of temporal synchrony in population fluctuations during early, incipient years (i.e. 1990–1996) suggested that outbreaking mountain pine beetle populations are largely independent at scales >200 km during non‐epidemic periods. However, during epidemic years (i.e. 1999–2003), populations were clearly synchronous across the entire province, even at distances of up to 900 km. The epicentral pattern of population development can be used to identify and prioritize adjacent landscape units for both reactive and proactive management strategies intended to minimize mountain pine beetle impacts.     

Interactions among phosphorus, pH and Eh in reforested mangroves, Vietnam: a three-dimensional spatial analysis

Sediment reduction is frequently proposed to increase available phosphorus (P) but several studies found also decreases. Another important factor for P liberation is the pH. We investigated the relative importance of Eh and pH on P cycling in reforested mangroves. Sediment P compounds, pH and Eh were analysed over depth along five transects of two areas in the Saigon River Delta and compared with leaf P levels. A three-dimensional spatial approach was used to investigate pH and Eh effects on P compound distribution at different sediment depth and locate layers of predominant P uptake. Along an inundation gradient, submergence durations of 254 to 2 days per year caused a large Eh gradient within the top 20 cm, whereas Eh response was small within 20–45 cm depth. At individual layers, a correlation between Eh and Al/Fe–P was only found in the upper depth interval (0–20 cm). No significant effect of Eh or Al/Fe–P on Morgan-P (available P) was detected. Minor effects on P composition changes by the Eh were caused through generally strongly reduced sediment at deeper layers (>20 cm). In contrast, pH variations produced extreme differences in both, Ca–P and available P content at layers of apparently predominant P uptake (>20 cm). As available P was correlated with Ca–P (p < 0.001), leaf P (p < 0.001) and pH (non linear correlation) it is likely that the pH sensitive Ca–P fraction is a more effective source for mangrove tree growth than Al/Fe–P. The predominant pH effect on P uptake within these reforested mangroves differs from a proposed reduction-governed P cycling in wetlands.     

Microsatellite variation in Central Africa: an analysis of intrapopulational and interpopulational genetic diversity

As a part of a research project on molecular variation in Central Africa, we have analyzed 10 microsatellites (CD4, CSFO, D3S1358, D18S51, D21S11, F13A1, FES, TH01, TPOX, and VWA) in the Bamileke and Ewondo from Cameroon and the Sanga and Mbenzele Pygmies from the Central African Republic (a total of 390 chromosomes). A statistically significant trend towards heterozygote deficiency was detected in the Mbenzele Pygmies. This was established through the use of powerful exact tests for the Hardy-Weinberg equilibrium. A certain degree of isolation and a small effective size may explain this finding. However, the lack of any substantial reduction in allelic diversity in the Mbenzele does not support the possibility that this group has a smaller effective size in evolutionary terms. A possible explanation based on ethnographic studies suggests that the gene flow from non-Pygmies to Pygmies could have been interrupted only in relatively recent times. The analysis of association between genotypes at pairs of independent loci indicates that the level of subheterogeneity is markedly lower in the Bamileke than in other sampled populations. This may be explained by the combined effect of larger population size, more rigid respect of clanic exogamy, and higher matrimonial mobility of the Bamileke. Finally, we have analyzed interpopulational relationships among our sampled populations and other Central African populations. The results are consistent with a previous study of protein loci (Spedini et al. 1999), which suggests the recent history of the Bamileke and Ewondo has led them to aquire a substantial genetic similarity. Furthermore, the Mbenzele Pygmies diverge from Biaka Pygmies, despite their common origin and geographical proximity. This is probably due to the differentiating effect of genetic drift, which is enhanced by the small effective size of Pygmy populations.