Variogram analysis of the spatial genetic structure of continuous populations using multilocus microsatellite data

Assessing fine-scale pondscape connectivity with amphibian eyes: An integrative approach using genomic and capture-mark-recapture data

In the face of habitat loss, preserving functional connectivity is essential to maintain genetic diversity and the demographic dynamics required for the viability of biotic communities. This requires knowledge of the dispersal behaviour of target species, which can be modelled as kernels, or probability density functions of dispersal distances at increasing geographic distances. We present an integrative approach to investigate the relationships between genetic connectivity and demographic parameters in organisms with low vagility focusing on five syntopic pond-breeding amphibians. We genotyped 1056 individuals of two anuran and three urodele species (1732-3913 SNPs per species) from populations located in a landscape comprising 64 ponds to characterize fine-scale genetic structure in a comparative framework, and combined these genetic data with information obtained in a previous 2-year capture-mark-recapture (CMR) study. Specifically, we contrasted graphs reconstructed from genomic data with connectivity graphs based on dispersal kernels and demographic information obtained from CMR data from previous studies, and assessed the effects of population size, population density, geographical distances, inverse movement probabilities and the presence of habitat patches potentially functioning as stepping stones on genetic differentiation. Our results show a significant effect of local population sizes on patterns of genetic differentiation at small spatial scales. In addition, movement records and cluster-derived kernels provide robust inferences on most likely dispersal paths that are consistent with genomic inferences on genetic connectivity. The integration of genetic and CMR data holds great potential for understanding genetic connectivity at spatial scales relevant to individual organisms, with applications for the implementation of management actions at the landscape level.

Effects of Forest Fragmentation on Connectivity and Genetic Diversity in an Endemic and an Invasive Rodent in Northwestern Madagascar

Habitat loss and fragmentation are of concern to conservation biologists worldwide. However, not all organisms are affected equally by these processes; thus, it is important to study the effects of living in fragmented habitats on species that differ in lifestyle and habitat requirements. In this study, we examined the dispersal and connectivity patterns of rodents, one endemic (Eliurus myoxinus) and one invasive (Rattus rattus), in two landscapes containing forest fragments and adjacent continuous forest patches in northwestern Madagascar. We generated genetic (RADseq) data for 66 E. myoxinus and 81 R. rattus individuals to evaluate differences in genetic diversity as well as inbreeding and connectivity in two landscapes. We found higher levels of inbreeding and lower levels of genetic diversity in E. myoxinus compared with R. rattus. We observed related dyads both within and between habitat patches and positive spatial autocorrelation at lower distance classes for both species, with a stronger pattern of spatial autocorrelation in R. rattus. Across each site, we identified contrasting migration rates for each species, but these did not correspond to habitat-matrix dichotomies. The relatively low genetic diversity in the endemic E. myoxinus suggests ecological constraints that require further investigation.

Evidence for high gene flow, nonrandom mating, and genetic bottlenecks of Ganoderma boninense infecting oil palm (Elaeis guineensis Jacq.) plantations in Malaysia and Indonesia

Ganoderma boninense, the causal agent of basal stem rot (BSR) disease, has been recognized as a major economic threat to commercial plantings of oil palm (Elaeis guineensis Jacq.) in Southeast Asia, which supplies 86% of the world's palm oil. High genetic diversity and gene flow among regional populations of 417 G. boninense isolates collected from Sabah, Sarawak, and Peninsular Malaysia (Malaysia) and Sumatra (Indonesia) were demonstrated using 16 microsatellite loci. Three genetic clusters and different admixed populations of G. boninense across regions were detected, and they appeared to follow the spread of the fungus from the oldest (Peninsular Malaysia and Sumatra) to younger generations of oil palm plantings (Sabah and Sarawak). Low spatial genetic differentiation of G. boninense (F_ST = 0.05) among the sampling regions revealed geographically nonrestricted gene dispersal, but isolation by distance was still evident. Analysis of molecular variance (AMOVA) confirmed the little to no genetic differentiation among the pathogen populations and the three genetic clusters defined by STRUCTURE and minimum spanning network. Despite G. boninense being highly outcrossing and spread by sexual spores, linkage disequilibrium was detected in 7 of the 14 populations. Linkage disequilibrium indicated that the reproduction of the fungus was not entirely by random mating and genetic drift could be an important structuring factor. Furthermore, evidence of population bottleneck was indicated in the oldest oil palm plantations as detected in genetic clusters 2 and 3, which consisted mainly of Peninsular Malaysia and Sumatra isolates. The population bottleneck or founding event could have arisen from either new planting or replanting after the removal of large number of palm hosts. The present study also demonstrated that migration and nonrandom mating of G. boninense could be important for survival and adaptation to new palm hosts.

The effect of sampling density and study area size on landscape genetics inferences for the Mississippi slimy salamander (Plethodon mississippi)

In landscape genetics, it is largely unknown how choices regarding sampling density and study area size impact inferences upon which habitat features impede vs. facilitate gene flow. While it is recommended that sampling locations be spaced no further apart than the average individual's dispersal distance, for low-mobility species, this could lead to a challenging number of sampling locations, or an unrepresentative study area. We assessed the effects of sampling density and study area size on landscape genetic inferences for a dispersal-limited amphibian, Plethodon mississippi, via analysis of nested datasets. Microsatellite-based genetic distances among individuals were divided into three datasets representing sparse sampling across a large study area, dense sampling across a small study area, or sparse sampling across the same small study area. These datasets were a proxy for gene flow (i.e., the response variable) in maximum-likelihood population effects models that assessed the nature and strength of their relationship with each of five land-use classes (i.e., potential predictor variables). Comparisons of outcomes were based on the rank order of effect, sign of effect (i.e., gene flow resistance vs. facilitation), spatial scale of effect, and functional relationship with gene flow. The best-fit model for each dataset had the same sign of effect for hardwood forests, manmade structures, and pine forests, indicating the impacts of these land-use classes on dispersal and gene flow in P. mississippi are robust to sampling scheme. Contrasting sampling densities led to a different inferred functional relationship between agricultural areas and gene flow. Study area size appeared to influence the scale of effect of manmade structures and the sign of effect of pine forests. Our findings provided evidence for an influence of sampling density, study area size, and sampling effort upon inferences. Accordingly, we recommend iterative subsampling of empirical datasets and continued investigation into the sensitivities of landscape genetic analyses using simulations.

Deep divergence between island populations in lichenized fungi

Macaronesia is characterized by a high degree of endemism and represents a noteworthy system to study the evolutionary history of populations and species. Here, we compare the population-genetic structure in three lichen-forming fungi, the widespread Lobaria pulmonaria and two Macaronesian endemics, L. immixta and L. macaronesica, based on microsatellites. We utilize population genetic approaches to explore population subdivision and evolutionary history of these taxa on the Canary Islands, Madeira, Azores, and the western Iberian Peninsula. A common feature in all species was the deep divergence between populations on the Azores, a pattern expected by the large geographic distance among islands. For both endemic species, there was a major split between archipelagos. In contrast, in the widespread L. pulmonaria, divergent individuals were distributed across multiple archipelagos, suggesting a complex evolutionary history involving repeated migration between islands and mainland.

A novel quantitative framework for riverscape genetics

Riverscape genetics, which applies concepts in landscape genetics to riverine ecosystems, lack appropriate quantitative methods that address the spatial autocorrelation structure of linear stream networks and account for bidirectional geneflow. To address these challenges, we present a general framework for the design and analysis of riverscape genetic studies. Our framework starts with the estimation of pairwise genetic distance at sample sites and the development of a spatially structured ecological network (SSEN) on which riverscape covariates are measured. We then introduce the novel bidirectional geneflow in riverscapes (BGR) model that uses principles of isolation-by-resistance to quantify the effects of environmental covariates on genetic connectivity, with spatial covariance defined using simultaneous autoregressive models on the SSEN and the generalized Wishart distribution to model pairwise distance matrices arising through a random walk model of geneflow. We highlight the utility of this framework in an analysis of riverscape genetics for brook trout (Salvelinus fontinalis) in north central Pennsylvania, USA. Using the fixation index (F_ST ) as the measure of genetic distance, we estimated the effects of 12 riverscape covariates on geneflow by evaluating the relative support of eight competing BGR models. We then compared the performance of the top-ranked BGR model to results obtained from comparable analyses using multiple regression on distance matrices (MRM) and the program STRUCTURE. We found that the BGR model had more power to detect covariate effects, particularly for variables that were only partial barriers to geneflow and/or uncommon in the riverscape, making it more informative for assessing patterns of population connectivity and identifying threats to species conservation. This case study highlights the utility of our modeling framework over other quantitative methods in riverscape genetics, particularly the ability to rigorously test hypotheses about factors that influence geneflow and probabilistically estimate the effect of riverscape covariates, including stream flow direction. This framework is flexible across taxa and riverine networks, is easily executable, and provides intuitive results that can be used to investigate the likely outcomes of current and future management scenarios.

Landscape genetics reveals unique and shared effects of urbanization for two sympatric pool-breeding amphibians

Metapopulation-structured species can be negatively affected when landscape fragmentation impairs connectivity. We investigated the effects of urbanization on genetic diversity and gene flow for two sympatric amphibian species, spotted salamanders (Ambystoma maculatum) and wood frogs (Lithobates sylvaticus), across a large (>35,000 km2) landscape in Maine, USA, containing numerous natural and anthropogenic gradients. Isolation-by-distance (IBD) patterns differed between the species. Spotted salamanders showed a linear and relatively high variance relationship between genetic and geographic distances (r = .057, p < .001), whereas wood frogs exhibited a strongly nonlinear and lower variance relationship (r = 0.429, p < .001). Scale dependence analysis of IBD found gene flow has its most predictable influence (strongest IBD correlations) at distances up to 9 km for spotted salamanders and up to 6 km for wood frogs. Estimated effective migration surfaces revealed contrasting patterns of high and low genetic diversity and gene flow between the two species. Population isolation, quantified as the mean IBD residuals for each population, was associated with local urbanization and less genetic diversity in both species. The influence of geographic proximity and urbanization on population connectivity was further supported by distance-based redundancy analysis and multiple matrix regression with randomization. Resistance surface modeling found interpopulation connectivity to be influenced by developed land cover, light roads, interstates, and topography for both species, plus secondary roads and rivers for wood frogs. Our results highlight the influence of anthropogenic landscape features within the context of natural features and broad spatial genetic patterns, in turn supporting the premise that while urbanization significantly restricts interpopulation connectivity for wood frogs and spotted salamanders, specific landscape elements have unique effects on these two sympatric species.

Multi-scale and multi-site resampling of a study area in spatial genetics: implications for flying insect species

The use of multiple sampling areas in landscape genetic analysis has been recognized as a useful way of generalizing the patterns of environmental effects on organism gene flow. It reduces the variability in inference which can be substantially affected by the scale of the study area and its geographic location. However, empirical landscape genetic studies rarely consider multiple sampling areas due to the sampling effort required. In this study, we explored the effects of environmental features on the gene flow of a flying long-horned beetle (Monochamus galloprovincialis) using a landscape genetics approach. To account for the unknown scale of gene flow and the multiple local confounding effects of evolutionary history and landscape changes on inference, we developed a way of resampling study areas on multiple scales and in multiple locations (sliding windows) in a single large-scale sampling design. Landscape analyses were conducted in 3*104 study areas ranging in scale from 220 to 1,000 km and spread over 132 locations on the Iberian Peninsula. The resampling approach made it possible to identify the features affecting the gene flow of this species but also showed high variability in inference among the scales and the locations tested, independent of the variation in environmental features. This method provides an opportunity to explore the effects of environmental features on organism gene flow on the whole and reach conclusions about general landscape effects on their dispersal, while limiting the sampling effort to a reasonable level.

Propagule size is not a good predictor for regional population subdivision or fine-scale spatial structure in lichenized fungi

Propagule size has important consequences on the genetic structure of wind-dispersed species, as species with small propagules have higher capability of long-distance dispersal. Here, we studied reproductive modes and compared local and regional population structures in three Macaronesian lichenized fungi differing in propagule size. First, we quantified size distribution of propagules in each species. Second, genotype simulations based on microsatellite data were used to infer the reproductive mode. Third, using spatial analysis and population genetic approaches, we quantified the local and regional scale genetic structures of the fungal species. The three species differed in size distributions of propagules. The majority of populations exhibited clonal reproductive mode. Identical reproductive modes occurred often across species in the same sites, implying a possible relationship between reproductive mode and local site conditions. Contrary to expectation, at the local scale, the species exhibited similar patterns of spatial autocorrelation in genotypes. However, in agreement with the expectation based on propagule size, the species with highest frequency of small vegetative propagules (L. pulmonaria) exhibited lowest regional genetic differentiation. Nevertheless, altogether, our results show that propagule size is not a good predictor of population subdivision in lichenized fungi, neither at local nor regional spatial scale.

The genetic structure of a Venturia inaequalis population in a heterogeneous host population composed of different Malus species

BACKGROUND

Adaptation, which induces differentiation between populations in relation to environmental conditions, can initiate divergence. The balance between gene flow and selection determines the maintenance of such a structure in sympatry. Studying these two antagonistic forces in plant pathogens is made possible because of the high ability of pathogens to disperse and of the strong selective pressures exerted by their hosts. In this article, we analysed the genetic structure of the population of the apple scab fungus, Venturia inaequalis, in a heterogeneous environment composed of various Malus species. Inferences were drawn from microsatellite and AFLP data obtained from 114 strains sampled in a single orchard on nine different Malus species to determine the forces that shape the genetic structure of the pathogen.

RESULTS

Using clustering methods, we first identified two specialist subpopulations: (i) a virulent subpopulation sampled on Malus trees carrying the Rvi6 resistance gene; and (ii) a subpopulation infecting only Malus trees that did not carry this resistance gene. A genome scan of loci on these two subpopulations did not detect any locus under selection. Additionally, we did not detect any other particular substructure linked to different hosts. However, an isolation-by-distance (IBD) pattern at the orchard scale revealed free gene flow within each subpopulation.

CONCLUSIONS

Our work shows a rare example of a very strong effect of a resistance gene on pathogen populations. Despite the high diversity of Malus hosts, the presence of Rvi6 seems sufficient to explain the observed genetic structure. Moreover, detection of an IBD pattern at the orchard scale revealed a very low average dispersal distance that is particularly significant for epidemiologists and landscape managers for the design of scab control strategies.

Hitchhiking with forests: population genetics of the epiphytic lichen Lobaria pulmonaria in primeval and managed forests in southeastern Europe

Availability of suitable trees is a primary determinant of range contractions and expansions of epiphytic species. However, switches between carrier tree species may blur co-phylogeographic patterns. We identified glacial refugia in southeastern Europe for the tree-colonizing lichen Lobaria pulmonaria, studied the importance of primeval forest reserves for the conservation of genetically diverse populations and analyzed differences in spatial genetic structure between primeval and managed forests with fungus-specific microsatellite markers. Populations belonged to either of two genepools or were admixed. Gene diversity was higher in primeval than in managed forests. At small distances up to 170 m, genotype diversity was lower in managed compared with primeval forests. We found significant associations between groups of tree species and two L. pulmonaria genepools, which may indicate "hitchhiking" of L. pulmonaria on forest communities during postglacial migration. Genepool B of L. pulmonaria was associated with European Beech (Fagus sylvatica) and we can hypothesize that genepool B survived the last glaciation associated within the refuge of European Beech on the Coastal and Central Dinarides. The allelic richness of genepool A was highest in the Alps, which is the evidence for a northern refuge of L. pulmonaria. Vicariant altitudinal distributions of the two genepools suggest intraspecific ecological differentiation.

Vertical and horizontal photobiont transmission within populations of a lichen symbiosis

Lichens are widespread symbioses and play important roles in many terrestrial ecosystems. The genetic structure of lichens is the result of the association between fungal and algal populations constituting the lichen thallus. Using eight fungus- and seven alga-specific highly variable microsatellite markers on within-population spatial genetic data from 62 replicate populations across Europe, North America, Asia and Africa, we investigated the contributions of vertical and horizontal transmission of the photobiont to the genetic structure of the epiphytic lichen Lobaria pulmonaria. Based on pairwise comparisons of multilocus genotypes defined separately for the mycobiont and for the photobiont, we inferred the transmission mode of the photobiont and the relative contribution of somatic mutation and recombination. After constraining the analysis of one symbiont to pairs of individuals with genetically identical symbiotic partners, we found that 77% of fungal and 70% of algal pairs were represented by clones. Thus, the predominant dispersal mode was by means of symbiotic vegetative propagules (vertical transmission), which dispersed fungal and algal clones co-dependently over a short distance, thus shaping the spatial genetic structure up to distances of 20m. Evidence for somatic mutation generating genetic diversity was found in both symbionts, accounting for 30% of pairwise comparisons in the alga and 15% in the fungus. While the alga did not show statistically significant evidence of recombination, recombination accounted for 7.7% of fungal pairs with identical algae. This implies that, even in a mostly vegetatively reproducing species, horizontal transmission plays a role in shaping the symbiotic association, as shown in many coral and other symbioses in nature.

Congruent genetic structure in the lichen-forming fungus Lobaria pulmonaria and its green-algal photobiont

The extent of codispersal of symbionts is one of the key factors shaping genetic structures of symbiotic organisms. Concordant patterns of genetic structure are expected in vertically transmitted symbioses, whereas horizontal transmission generally uncouples genetic structures unless the partners are coadapted. Here, we compared the genetic structures of mutualists, the lichen-forming fungus Lobaria pulmonaria and its primary green-algal photobiont, Dictyochloropsis reticulata. We performed analysis of molecular variance and variogram analysis to compare genetic structures between symbiosis partners. We simulated the expected number of multilocus-genotype recurrences to reveal whether the distribution of multilocus genotypes of either species was concordant with panmixia. Simulations and tests of linkage disequilibrium provided compelling evidence for the codispersal of mutualists. To test whether genotype associations between symbionts were consistent with randomness, as expected under horizontal transmission, we simulated the recurrence of fungal-algal multilocus genotype associations expected by chance. Our data showed nonrandom associations of fungal and algal genotypes. Either vertical transmission or horizontal transmission coupled with coadaptation between symbiont genotypes may have created these nonrandom associations. This study is among the first to show codispersal and highly congruent genetic structures in the partners of a lichen mutualism.

Quantifying the lag time to detect barriers in landscape genetics

Understanding how spatial genetic patterns respond to landscape change is crucial for advancing the emerging field of landscape genetics. We quantified the number of generations for new landscape barrier signatures to become detectable and for old signatures to disappear after barrier removal. We used spatially explicit, individual-based simulations to examine the ability of an individual-based statistic [Mantel's r using the proportion of shared alleles' statistic (Dps)] and population-based statistic (FST ) to detect barriers. We simulated a range of movement strategies including nearest neighbour dispersal, long-distance dispersal and panmixia. The lag time for the signal of a new barrier to become established is short using Mantel's r (1-15 generations). FST required approximately 200 generations to reach 50% of its equilibrium maximum, although G'ST performed much like Mantel's r. In strong contrast, FST and Mantel's r perform similarly following the removal of a barrier formerly dividing a population. Also, given neighbour mating and very short-distance dispersal strategies, historical discontinuities from more than 100 generations ago might still be detectable with either method. This suggests that historical events and landscapes could have long-term effects that confound inferences about the impacts of current landscape features on gene flow for species with very little long-distance dispersal. Nonetheless, populations of organisms with relatively large dispersal distances will lose the signal of a former barrier within less than 15 generations, suggesting that individual-based landscape genetic approaches can improve our ability to measure effects of existing landscape features on genetic structure and connectivity.

Identity and genetic structure of the photobiont of the epiphytic lichen Ramalina menziesii on three oak species in southern California

Lichens, a classic example of an obligate symbiosis between fungi and photobionts (which could be algae or cyanobacteria), are abundant in many terrestrial ecosystems. The genetic structure of the photobiont population found in association with a lichen-forming fungal species could be affected by fungal reproductive mode and by the spatial extent of gene flow in the photobiont. Using DNA sequences from one nuclear ribosomal and two chloroplast loci, we analyzed the genetic structure of the photobiont associated with the fungus Ramalina menziesii at an oak woodland study site in southern California. We had previously shown that the fungus exhibited no genetic structure among four local sites or three phorophyte species. Our goals were to identify the photobiont species and assess its genetic structure. We found that R. menziesii was highly specific in its photobiont choice and associated with one alga, Trebouxia decolorans. In contrast to the fungal population, we found significant differentiation among the algae sampled on three oak species and little genetic structure among the sites for two of the three algal loci. We hypothesize that R. menziesii is locally adapted to the phorophyte species through habitat specialization in the algal partner of the symbiosis.

Statistical methods in spatial genetics

The joint analysis of spatial and genetic data is rapidly becoming the norm in population genetics. More and more studies explicitly describe and quantify the spatial organization of genetic variation and try to relate it to underlying ecological processes. As it has become increasingly difficult to keep abreast with the latest methodological developments, we review the statistical toolbox available to analyse population genetic data in a spatially explicit framework. We mostly focus on statistical concepts but also discuss practical aspects of the analytical methods, highlighting not only the potential of various approaches but also methodological pitfalls.

A review of techniques for spatial modeling in geographical, conservation and landscape genetics

Most evolutionary processes occur in a spatial context and several spatial analysis techniques have been employed in an exploratory context. However, the existence of autocorrelation can also perturb significance tests when data is analyzed using standard correlation and regression techniques on modeling genetic data as a function of explanatory variables. In this case, more complex models incorporating the effects of autocorrelation must be used. Here we review those models and compared their relative performances in a simple simulation, in which spatial patterns in allele frequencies were generated by a balance between random variation within populations and spatially-structured gene flow. Notwithstanding the somewhat idiosyncratic behavior of the techniques evaluated, it is clear that spatial autocorrelation affects Type I errors and that standard linear regression does not provide minimum variance estimators. Due to its flexibility, we stress that principal coordinate of neighbor matrices (PCNM) and related eigenvector mapping techniques seem to be the best approaches to spatial regression. In general, we hope that our review of commonly used spatial regression techniques in biology and ecology may aid population geneticists towards providing better explanations for population structures dealing with more complex regression problems throughout geographic space.

Genetic diversity and spatial genetic structure within a population of an aromatic shrub, Lippia origanoides (Verbenaceae), in the Chicamocha Canyon, northeastern Colombia

The geographical scale of genetic structure in a continuous population is highly dependent on its breeding system and dispersion capabilities, and this knowledge is important for the study of population dynamics as well as for conservation purposes. In the present study, spatial autocorrelation statistics and intersimple sequence repeat (ISSR) markers were used to describe the genetic structure of a natural population of a prominent aromatic plant, Lippia origanoides, native to the Chicamocha Canyon in northeastern Colombia. For this purpose, individuals were sampled from two localities within the Chicamocha Canyon, where the species is abundant and continuously distributed. Cluster (principal coordinate analysis (PCO) and unweighted pair group method using arithmetic averages (UPGMA)), analysis of molecular variance (AMOVA) and Bayesian analyses revealed a low level of genetic differentiation among the two localities, suggesting that they belong to a single population. Genetic diversity levels in this population, described as the percentage of polymorphic loci (P=86.21%) and quantified using Shannon's diversity index (I=0.453) and the average panmictic heterozygosity (HB=0.484), were shown to be comparable to or higher than that in other plant species with allogamous breeding systems and to other related Verbenaceae species. Fine-scale autocorrelation analyses showed a pattern consistent with the classical model of isolation by distance with moderate but significant levels of local spatial structure. Our results suggest that sampling individuals at distances greater than approximately 1.2 km may result in the collection of different genotypes, which could help preserve the levels of genetic diversity in a propagation programme. The causes of this spatial pattern are currently unknown and could be influenced by many contemporary factors such as restricted seed dispersal and/or short-distance pollen movement, among others.

A computer program to simulate multilocus genotype data with spatially autocorrelated allele frequencies

Many models for inference of population genetic parameters are based on the assumption that the data set at hand consists of groups displaying within-group Hardy-Weinberg equilibrium at individual loci and linkage equilibrium between loci. This assumption is commonly violated by the presence of within-group spatial structure arising from nonrandom mating of individuals due to isolation by distance (IBD). This paper proposes a model and simulation method implemented in a computer program to flexibly simulate data displaying such patterns. The program permits displaying of smooth spatial variations of allele frequencies due to IBD and more abrupt variations due to presence of strong barriers to gene flow. It is useful in assessing performance of various statistical inference methods and in designing spatial sampling schemes. This is shown by a simulation study aimed at assessing the extent to which IBD patterns affect accuracy of cluster inferences performed in models assuming panmixia. The program is also used to study the effects of spatial sampling scheme (e.g. sampling individuals in clumps or uniformly across the spatial domain). The accuracy of such inferences is assessed in terms of number of inferred populations, assignment of individuals to populations and location of borders between populations. The effect of spatial sampling was weak while the effect of IBD may be substantial, leading to the inference of spurious populations, especially when IBD was strong with respect to the size of the sampling domain. The model and program are new and have been embedded in the R package Geneland, for user convenience and compliance with existing data formats.

Non-parametric smoothing of multivariate genetic distances in the analysis of spatial population structure at fine scale

Species dispersal studies provide valuable information in biological research. Restricted dispersal may give rise to a non-random distribution of genotypes in space. Detection of spatial genetic structure may therefore provide valuable insight into dispersal. Spatial structure has been treated via autocorrelation analysis with several univariate statistics for which results could dependent on sampling designs. New geostatistical approaches (variogram-based analysis) have been proposed to overcome this problem. However, modelling parametric variograms could be difficult in practice. We introduce a non-parametric variogram-based method for autocorrelation analysis between DNA samples that have been genotyped by means of multilocus-multiallele molecular markers. The method addresses two important aspects of fine-scale spatial genetic analyses: the identification of a non-random distribution of genotypes in space, and the estimation of the magnitude of any non-random structure. The method uses a plot of the squared Euclidean genetic distances vs. spatial distances between pairs of DNA-samples as empirical variogram. The underlying spatial trend in the plot is fitted by a non-parametric smoothing (LOESS, Local Regression). Finally, the predicted LOESS values are explained by segmented regressions (SR) to obtain classical spatial values such as the extent of autocorrelation. For illustration we use multivariate and single-locus genetic distances calculated from a microsatellite data set for which autocorrelation was previously reported. The LOESS/SR method produced a good fit providing similar value of published autocorrelation for this data. The fit by LOESS/SR was simpler to obtain than the parametric analysis since initial parameter values are not required during the trend estimation process. The LOESS/SR method offers a new alternative for spatial analysis.

Local genetic structure in a North American epiphytic lichen, Ramalina menziesii (Ramalinaceae)

Epiphytic lichens possess unique life history traits that can have conflicting effects on genetic structure: symbiotic mutualism between a fungus with its algal or cyanobacterial photobiont, association with a host plant, and ability to reproduce sexually and asexually. Our study species, Ramalina menziesii, has small ascospores that can facilitate long-distance gene movement, and it is capable of clonal reproduction. The goals of this study are to test whether different haplotypes were differentially distributed across host plant species, to look for evidence of asexual vs. sexual reproduction, and to assess the local genetic structure of the population. We sampled individuals from multiple trees of three oak species in four lichen subpopulations within a savanna ecosystem. Using DNA sequence data from four fungal nuclear loci, we found no tendency for host specialization. Alleles were randomly distributed across subpopulations. The frequency of multilocus genotypes was consistent with a randomly mating population. Sexual reproduction involving relichenization appeared to be the predominant mode of reproduction of R. menziesii at this study site. We found no significant local genetic structure suggesting widespread gene flow at the local scale. The genetic structure of this lichen is comparable to that of widely distributed epiphytic plants.

Finding the relevant scale: clonality and genetic structure in a marine invertebrate (Crambe crambe, Porifera)

Important changes in genetic relatedness may occur at extremely small scales in benthic invertebrates, providing key information about structuring processes in populations of these organisms. We performed a small-scale study of the population structure of the sponge Crambe crambe, in which 177 individuals from the same rocky wall (interindividual distances from 0 to 7 m) were genotyped using six microsatellite markers. 101 sponges had unique genotypes and the remaining 76 individuals formed 24 groups of sponges sharing genotypes (clones). Mean intraclone distances were found to be c. 20 cm. Spatial autocorrelation analyses showed a drastic decrease in genetic relatedness over the first 100 cm of distance. If the contribution of clonality to this pattern was eliminated, the trend was attenuated, but remained a marked one and was still significant within the first distance classes (30-40 cm). Estimated mean dispersal distances per generation were c. 35 cm, and neighbourhood sizes were estimated at c. 33 sponges. Genetic similarities with sponges of the same locality, or from other Mediterranean localities, were within the same range as those found in sponges 2-7 m apart. It is concluded that asexual reproduction plays an important role in structuring populations in this species. However, over and above the effects of clonality, a strong fine-scale genetic structure was present at distances in the range of tens of centimetres, probably as a result of short dispersal of larvae. This fine-scale genetic structure may be common in invertebrates with lecitotrophic larvae.

Landscape-level gene flow in Lobaria pulmonaria, an epiphytic lichen

Epiphytes are strongly affected by the population dynamics of their host trees. Owing to the spatio-temporal dynamics of host tree populations, substantial dispersal rates--corresponding to high levels of gene flow--are needed for populations to persist in a landscape. However, several epiphytic lichens have been suggested to be dispersal-limited, which leads to the expectation of low gene flow at the landscape scale. Here, we study landscape-level genetic structure and gene flow of a putatively dispersal-limited epiphytic lichen, Lobaria pulmonaria. The genetic structure of L. pulmonaria was quantified at three hierarchical levels, based on 923 thalli collected from 41 plots situated within a pasture-woodland landscape and genotyped at six fungal microsatellite loci. We found significant isolation by distance, and significant genetic differentiation both among sampling plots and among trees. Landscape configuration, i.e. the effect of a large open area separating two forested regions, did not leave a traceable pattern in genetic structure, as assessed with partial Mantel tests and analysis of molecular variance. Gene pools were spatially intermingled in the pasture-woodland landscape, as determined by Bayesian analysis of population structure. Evidence for local gene flow was found in a disturbed area that was mainly colonized from nearby sources. Our analyses indicated high rates of gene flow of L. pulmonaria among forest patches, which may reflect the historical connectedness of the landscape through gene movement. These results support the conclusion that dispersal in L. pulmonaria is rather effective, but not spatially unrestricted.

Putting the ‘landscape’ in landscape genetics

Landscape genetics has emerged as a new research area that integrates population genetics, landscape ecology and spatial statistics. Researchers in this field can combine the high resolution of genetic markers with spatial data and a variety of statistical methods to evaluate the role that landscape variables play in shaping genetic diversity and population structure. While interest in this research area is growing rapidly, our ability to fully utilize landscape data, test explicit hypotheses and truly integrate these diverse disciplines has lagged behind. Part of the current challenge in the development of the field of landscape genetics is bridging the communication and knowledge gap between these highly specific and technical disciplines. The goal of this review is to help bridge this gap by exposing geneticists to terminology, sampling methods and analysis techniques widely used in landscape ecology and spatial statistics but rarely addressed in the genetics literature. We offer a definition for the term "landscape genetics", provide an overview of the landscape genetics literature, give guidelines for appropriate sampling design and useful analysis techniques, and discuss future directions in the field. We hope, this review will stimulate increased dialog and enhance interdisciplinary collaborations advancing this exciting new field.

Effect of disturbances on the genetic diversity of an old-forest associated lichen

Lichens associated with old forest are commonly assumed to be negatively affected by tree logging or natural disturbances. However, in this study performed in a spruce-dominated sylvopastoral landscape in the Swiss Jura Mountains, we found that genetic diversity of the epiphytic old-forest lichen Lobaria pulmonaria depends on the type of disturbance. We collected 923 thalli from 41 sampling plots of 1 ha corresponding to the categories stand-replacing disturbance (burnt), intensive logging (logged) and uneven-aged forestry (uneven-aged), and analysed the thalli at six mycobiont-specific microsatellite loci. We found evidence for multiple independent immigrations into demes located in burnt and logged areas. Using spatial autocorrelation methods, the spatial scale of the genetic structure caused by the clonal and recombinant component of genetic variation was determined. Spatial autocorrelation of genotype diversity was strong at short distances up to 50 m in logged demes, up to 100 m in uneven-aged demes, with the strongest autocorrelation up to 150 m for burnt demes. The spatial autocorrelation was predominantly attributed to clonal dispersal of vegetative propagules. After accounting for the clonal component, we did not find significant spatial autocorrelation in gene diversity. This pattern may indicate low dispersal ranges of clonal propagules, but random dispersal of sexual ascospores. Genetic diversity was highest in logged demes, and lowest in burnt demes. Our results suggest that genetic diversity of epiphytic lichen demes may not necessarily be impacted by stand-level disturbances for extended time periods.