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.     

A likelihood-based approach to estimating and testing for isolation by distance

Simple regression of genetic similarities between pairs of populations on their corresponding geographic distances is frequently used to detect the presence of isolation by distance (IBD). However, these pairwise values are obviously not independent and there is no parametric procedure for estimating and testing for the IBD intercepts and slopes based on standard regression theory. Nonparametric tests, such as the Mantel test, and resampling techniques, such as bootstrapping, have been exploited with limited success. Here, I describe a likelihood-based analysis to allow for simultaneously detecting patterns of correlated residuals and estimating and testing for the presence of IBD. It is shown, through the analysis of two molecular datasets in pine species, that different covariance structures of the residuals exist. More over, the likelihood ratio tests under these covariance structures are less sensitive to the presence of IBD than the Mantel test and the simple regression analysis but more sensitive than the bootstrap and jackknife samples over independent populations or population pairs. Because the likelihood analysis directly models and accounts for nonindependence of residuals, it should legitimately detect the presence of IBD, thereby allowing for accurate inferences about evolutionary and demographic processes influencing the extent and patterns of IBD.     

Spatial genetic structuring in a widespread wetland plant on a plateau: Effects of elevation-driven geographic isolation and environmental heterogeneity

1. Highlands are ideal research areas for improving our understanding of the influence of ecological factors on the diversity and spatial patterns of natural species. Elevation-driven physical and environmental isolation greatly affect the evolution of plants. The mechanisms and essential drivers underlying these processes may differ among research scales, habitats and landscapes. Wetlands are important elements of the Qinghai–Tibetan Plateau, which is the highest plateau in the world, and these habitats harbour high aquatic organismal diversity. However, how the environments shape the genetic variation and structure of hydrophilous plants is poorly understood.
2. Using microsatellite markers and a chloroplast fragment, we quantified the genetic diversity and spatial genetic pattern of Stuckenia filiformis, one of the most widespread aquatic plants on the plateau. The relative contributions of geography, climate and local conditions to intra- and interpopulation variation were estimated. The results showed that intrapopulation genetic variation of the plant is moderate to high and not constrained by high-altitude environments. Topographical isolation mainly contributes to the genetic structure of S. filiformis, as inferred by simple sequence repeats and chloroplast DNA data. Significant effects of environmental variables on the spatial genetic patterns of this freshwater species were also suggested by landscape genetic analysis.
3. Infrequent long-distance dispersal, sexual recruitment and annual growth are probably important for the maintenance and distribution of this variation. Our findings imply a combined effect of geography and elevation-driven environmental heterogeneity on the evolution of aquatic organisms in highlands.

     

Isolation by environmental distance in mobile marine species: molecular ecology of franciscana dolphins at their southern range

The assessment of population structure is a valuable tool for studying the ecology of endangered species and drafting conservation strategies. As we enhance our understanding about the structuring of natural populations, it becomes important that we also understand the processes behind these patterns. However, there are few rigorous assessments of the influence of environmental factors on genetic patterns in mobile marine species. Given their dispersal capabilities and localized habitat preferences, coastal cetaceans are adequate study species for evaluating environmental effects on marine population structure. The franciscana dolphin, a rare coastal cetacean endemic to the Western South Atlantic, was studied to examine these issues. We analysed genetic data from the mitochondrial DNA and 12 microsatellite markers for 275 franciscana samples utilizing frequency‐based, maximum‐likelihood and Bayesian algorithms to assess population structure and migration patterns. This information was combined with 10 years of remote sensing environmental data (chlorophyll concentration, water turbidity and surface temperature). Our analyses show the occurrence of genetically isolated populations within Argentina, in areas that are environmentally distinct. Combined evidence of genetic and environmental structure suggests that isolation by distance and a process here termed isolation by environmental distance can explain the observed correlations. Our approach elucidated important ecological and conservation aspects of franciscana dolphins, and has the potential to increase our understanding of ecological processes influencing genetic patterns in other marine species.     

Genetic differentiation between individuals

I describe a method of analysis of genetic differentiation which is suitable for the comparison of genetic and demographic estimates of the ‘neighborhood size’– more precisely the product of density and second moment of dispersal distance σ²– in continuous populations sampled at the smallest scale. This method is based on results of models of isolation by distance common to a wide variety of dispersal distances. The performance of this method is tested by simulation for some highly leptokurtic dispersal distributions, and it is applied to a previous study of a kangaroo‐rat (Dipodomys spectabilis) population. In this case, genetic and demographic estimates are within a factor of two from each other. Thus, in line with some previous examples, this study shows that a better agreement may be attained than is usually recognized between genetic and demographic estimates.     

Drivers of population genetic differentiation in the wild: isolation by dispersal limitation,isolation by adaptation and isolation by colonization

Empirical population genetic studies have been dominated by a neutralist view, according to which gene flow and drift are the main forces driving population genetic structure in nature. The neutralist view in essence describes a process of isolation by dispersal limitation (IBDL) that generally leads to a pattern of isolation by distance (IBD). Recently, however, conceptual frameworks have been put forward that view local genetic adaptation as an important driver of population genetic structure. Isolation by adaptation (IBA) and monopolization (M) posit that gene flow among natural populations is reduced as a consequence of local genetic adaptation. IBA stresses that effective gene flow is reduced among habitats that show dissimilar ecological characteristics, leading to a pattern of isolation by environment. In monopolization, local genetic adaptation of initial colonizing genotypes results in a reduction in gene flow that fosters the persistence of founder effects. Here, we relate these different processes driving landscape genetic structure to patterns of IBD and isolation by environment (IBE). We propose a method to detect whether IBDL, IBA and M shape genetic differentiation in natural landscapes by studying patterns of variation at neutral and non‐neutral markers as well as at ecologically relevant traits. Finally, we reinterpret a representative number of studies from the recent literature by associating patterns to processes and identify patterns associated with local genetic adaptation to be as common as IBDL in structuring regional genetic variation of populations in the wild. Our results point to the importance of quantifying environmental gradients and incorporating ecology in the analysis of population genetics.     

Spatial correlations at different spatial scales are themselves highly correlated in isolation by distance processes

Although many properties of spatial autocorrelation statistics are well characterized, virtually nothing is known about possible correlations among values at different spatial scales, which ultimately would influence how inferences about spatial genetics are made at multiple spatial scales. This article reports the results of stochastic space-time simulations of isolation by distance processes, having a very wide range of amounts of dispersal for plants or animals, and analyses of the correlations among Moran's I-statistics for different mutually exclusive distance classes. In general, the stochastic correlations are extremely large (>0.90); however, the correlations bear a complex relationship with level of dispersal, spatial scale and spatial lag between distance classes. The correlations are so large that any existing or conceived statistical method that employs more than one distance class (or spatial scale) should not ignore them. This result also suggests that gains in statistical power via increasing sample size are limited, and that increasing numbers of assayed loci generally should be preferred. To the extent that sampling error for real data sets can be treated as white noise, it should be possible to account for stochastic correlations in formulating more precise statistical methods. Further, while the current results are for isolation by distance processes, they provide some guidance for some more complex stochastic space-time processes of landscape genetics. Moreover, the results hold for several popular measures other than Moran's I. In addition, in the results, the signal to noise ratios strongly decreased with distance, which also has several implications for optimal statistical methods using correlations at multiple spatial scales.     

The names of Spain: a study of the isonymy structure of Spain

In order to estimate the isonymy structure of Spain, we studied surname distribution in 283 Spanish towns based on 3.625 million telephone users selected from 6.328 million users, downloaded from a commercial CD-ROM which contains all 13 million users in the country. Since in Spain the surname is made by the paternal and the maternal surname, it was possible to classify surnames according to parental origin. Two matrices of isonymy distances, one for paternal and one for maternal surnames, were constructed and tested for correlation with geographic distance. For the whole of Spain, Euclidean distance was significantly but weakly correlated with geographic distance both for paternal and maternal surnames, with r = 0.205 +/- 0.013 and r = 0.263 +/- 0.012, respectively. Two dendrograms of the 283 sampled towns were built from the two matrices of Euclidean distance. They are largely colinear. Four main clusters identified by the dendrograms are correlated with geography. Given the surname structure of Spain, we were able to calculate from isonymy and for each town 1). total or expressed inbreeding, 2). random or expected inbreeding, and 3). local inbreeding. Total inbreeding, F(IT), was highest in the North Atlantic regions and lowest along the Mediterranean Coast. The lowest levels were found in Andalusia, Catalunyia, Valencia, and Navarra. Random inbreeding, F(ST), had a similar geographical pattern. Local inbreeding, F(IS), was relatively uniform in the whole of Spain. In towns, random inbreeding dominates over local inbreeding. From the analysis, it emerges that the northwestern area of Spain is the most inbred.     

Male and female contributions to behavioral isolation in darters as a function of genetic distance and color distance

Determining which reproductive isolating barriers arise first between geographically isolated lineages is critical to understanding allopatric speciation. We examined behavioral isolation among four recently diverged allopatric species in the orangethroat darter clade (Etheostoma: Ceasia). We also examined behavioral isolation between each Ceasia species and the sympatric rainbow darter Etheostoma caeruleum. We asked (1) is behavioral isolation present between allopatric Ceasia species, and how does this compare to behavioral isolation with E. caeruleum, (2) does male color distance and/or genetic distance predict behavioral isolation between species, and (3) what are the relative contributions of female choice, male choice, and male competition to behavioral isolation? We found that behavioral isolation, genetic differentiation, and male color pattern differentiation were present between allopatric Ceasia species. Males, but not females, discerned between conspecific and heterospecific mates. Males also directed more aggression toward conspecific rival males. The high levels of behavioral isolation among Ceasia species showed no obvious pattern with genetic distance or male color distance. However, when the E. caeruleum was included in the analysis, an association between male aggression and male color distance was apparent. We discuss the possibility that reinforcement between Ceasia and E. caeruleum is driving behavioral isolation among allopatric Ceasia species.     

Hierarchical components of genetic variation at a species boundary: population structure in two sympatric varieties of Lupinus microcarpus (Leguminosae)

Lupinus microcarpus is a self-compatible annual plant that forms a species complex of morphologically variable but indeterminate varieties. In order to examine the hypothesis that varieties of L. microcarpus comprise genetically differentiated and reproductively isolated species, populations of L. microcarpus var. horizontalis and var. densiflorus were sampled from an area of sympatry in central California and genotyped using six microsatellite loci. Bayesian clustering divided the total sample into two groups corresponding to the named varieties with extremely low levels of inferred coancestry. Similarly, maximum likelihood and distance methods for genetic assignment placed individuals in two nonoverlapping groups. Evidence for isolation by distance (IBD) within each variety was found at shorter distance classes, but varieties remained differentiated in sympatry. Furthermore, coalescent estimates of divergence time indicate separation within the past 950-5050 generations, with minimal gene flow after divergence. A four-level hierarchical analysis of molecular variance (amova) found significant levels of genetic differentiation among varieties (theta(P) = 0.292), populations within varieties (theta(S) = 0.449), subpopulations within populations (theta(SS) = 0.623), and individuals within subpopulations (f = 0.421); but the greatest degree of differentiation was at the subpopulation level. Although it is sometimes assumed that the magnitude of genetic differences (e.g. F(ST)) should be greater between species than among populations or subpopulations of the same species, shared ancestral polymorphism may lead to relatively low levels of differentiation at the species level, even as the stochastic effects of genetic drift generate higher levels of differentiation at lower hierarchical levels. These results suggest that L. microcarpus var. horizontalis and var. densiflorus are recently diverged yet reproductively isolated species, with high levels of inbreeding resulting from the combined effects of limited gene flow, demographic bottlenecks, and partial selfing in finite, geographically structured populations.     

Surnames in Chile: a study of the population of Chile through isonymy

In Chile, the Hispanic dual surname system is used. To describe the isonymic structure of this country, the distribution of 16,277,255 surnames of 8,178,209 persons was studied in the 15 regions, the 54 provinces, and the 346 communes of the nation. The number of different surnames found was 72,667. Effective surname number (Fisher's α) for the entire country was 309.0, the average for regions was 240.8 ± 17.6, for provinces 209.2 ± 8.9, and for communes 178.7 ± 4.7. These values display a variation of inbreeding between administrative levels in the Chilean population, which can be attributed to the 'Prefecture effect' of Nei and Imaizumi. Matrices of isonymic distances between units within administrative levels were tested for correlation with geographic distance. The correlations were highest for provinces (r = 0.630 ± 0.019 for Euclidean distance) and lowest for communes (r = 0.366 ± 0.009 for Lasker's). The geographical distribution of the first three-dimensions of the Euclidean distance matrix suggests that population diffusion may have taken place from the north of the country toward the center and south. The prevalence of European plus European-Amerindian (95.4%) over Amerindian ethnicity (4.6%, CIA World Factbook) is compatible with diffusion of Caucasian groups over a low-density area populated by indigenous groups. The significant excess of maternal over paternal indigenous surnames indicates some asymmetric mating between nonAmerindian and Amerindian Chileans. The available studies of Y-markers and mt-markers are in agreement with this asymmetry. In the present work, we investigate the Chilean population with the aim of detecting its structure through the study of isonymy (Crow and Mange,1965) in the three administrative levels of the nation, namely 15 regions, 54 provinces, and 346 communes.     

Isolation by environment

The interactions between organisms and their environments can shape distributions of spatial genetic variation, resulting in patterns of isolation by environment (IBE) in which genetic and environmental distances are positively correlated, independent of geographic distance. IBE represents one of the most important patterns that results from the ways in which landscape heterogeneity influences gene flow and population connectivity, but it has only recently been examined in studies of ecological and landscape genetics. Nevertheless, the study of IBE presents valuable opportunities to investigate how spatial heterogeneity in ecological processes, agents of selection and environmental variables contributes to genetic divergence in nature. New and increasingly sophisticated studies of IBE in natural systems are poised to make significant contributions to our understanding of the role of ecology in genetic divergence and of modes of differentiation both within and between species. Here, we describe the underlying ecological processes that can generate patterns of IBE, examine its implications for a wide variety of disciplines and outline several areas of future research that can answer pressing questions about the ecological basis of genetic diversity.     

Surnames in Argentina: a population study through isonymy

In order to study the isonymic structure of Argentina, the surname distributions of 22.6 million electors registered for the year 2001 were analyzed in the 24 districts (distritos) and 541 municipalities (municipios) of the country. The number of different surnames found was 414,441. Matrices of isonymic distances between districts were constructed and tested for correlation with the geographic distance between the capital towns of the districts. We found that, for the whole of Argentina, Euclidean distance was correlated with the log of geographic distance (r=0.480+/- 0.067). A dendrogram of the 24 regions was built from the matrix of Euclidean distances, using the UPGMA method. The clusters identified by the dendrogram are coincident with conterminous geographical regions of the country. Random inbreeding calculated from isonymy, F(ST), was highest in La Rioja, Corrientes, and Santiago del Estero. It was lowest in the area of Buenos Aires and in the north-central region of Santa Fé. Average Fisher's alpha for municipalities was 358; for districts, it was 422; and for Argentina as a unit, it was 602. The geographical distribution of alpha in 541 municipalities, high in the east and lower in the west of the country, is compatible with the settlement in the 20th century of subsequent waves of immigrants moving from the North Atlantic coast toward the foot of the Andes and toward the south. The present structure of Argentina indicates that migration dominates over drift.     

Genome‐wide epigenetic isolation by environment in a widespread Anolis lizard

Epigenetic changes can provide a pathway for organisms to respond to local environmental conditions by influencing gene expression. However, we still know little about the spatial distribution of epigenetic variation in natural systems, how it relates to the distribution of genetic variation and the environmental structure of the landscape, and the processes that generate and maintain it. Studies examining spatial patterns of genetic and epigenetic variation can provide valuable insights into how ecological and population processes contribute to epigenetic divergence across heterogeneous landscapes. Here, we perform a comparative analysis of spatial genetic and epigenetic variation based on 8,459 single nucleotide polymorphisms (SNPs) and 8,580 single methylation variants (SMVs) from eight populations of the Puerto Rican crested anole, Anolis cristatellus, an abundant lizard in the adaptive radiations of anoles on the Greater Antilles that occupies a diverse range of habitats. Using generalized dissimilarity modelling and multiple matrix regression, we found that genome‐wide epigenetic differentiation is strongly correlated with environmental divergence, even after controlling for the underlying genetic structure. We also detected significant associations between key environmental variables and 96 SMVs, including 42 located in promoter regions or gene bodies. Our results suggest an environmental basis for population‐level epigenetic differentiation in this system and contribute to better understanding how environmental gradients structure epigenetic variation in nature.     

Influence of mutational and sampling factors on the estimation of demographic parameters in a "continuous" population under isolation by distance

In numerous species, individual dispersal is restricted in space so that "continuous" populations evolve under isolation by distance. A method based on individual genotypes assuming a lattice population model was recently developed to estimate the product Dsigma2, where D is the population density and sigma2 is the average squared parent-offspring distance. We evaluated the influence on this method of both mutation rate and mutation model, with a particular reference to microsatellite markers, as well as that of the spatial scale of sampling. Moreover, we developed and tested a nonparametric bootstrap procedure allowing the construction of confidence intervals for the estimation of Dsigma2. These two objectives prompted us to develop a computer simulation algorithm based on the coalescent theory giving individual genotypes for a continuous population under isolation by distance. Our results show that the characteristics of mutational processes at microsatellite loci, namely the allele size homoplasy generated by stepwise mutations, constraints on allele size, and change of slippage rate with repeat number, have little influence on the estimation of Dsigma2. In contrast, a high genetic diversity (approximately 0.7-0.8), as is commonly observed for microsatellite markers, substantially increases the precision of the estimation. However, very high levels of genetic diversity (>0.85) were found to bias the estimation. We also show that statistics taking into account allele size differences give unreliable estimations (i.e., high variance of Dsigma2 estimation) even under a strict stepwise mutation model. Finally, although we show that this method is reasonably robust with respect to the sampling scale, sampling individuals at a local geographical scale gives more precise estimations of Dsigma2.     

IBDSim: a computer program to simulate genotypic data under isolation by distance

IBDSim is a package for the simulation of genotypic data under isolation by distance. It is based on a backward 'generation by generation' coalescent algorithm allowing the consideration of various isolation by distance models with discrete subpopulations as well as continuous populations. Many dispersal distributions can be considered as well as heterogeneities in space and time of the demographic parameters. Typical applications of our program include (i) the study of the effect of various sampling, mutational and demographic factors on the pattern of genetic variation; and (ii) the production of test data sets to assess the influence of these factors on inferential methods available to analyse genotypic data.     

Island isolation and habitat heterogeneity correlate with DNA variation in a marine snail (Littorina saxatilis)

The null assumption of molecular variation is that most of it is neutral to natural selection. This is in contrast to variation in morphological traits that we generally assume is maintained by selection, and therefore often by selection coupled to environmental heterogeneity in time and space. Examples of molecular variation that vary over habitat-shifts, particularly in allozymes, show that the relative impact of non-neutral variation as compared to neutral variation might be substantial in some systems. To assess the importance of habitat-generated variation in relation to variation generated by random processes in nuclear DNA markers at small spatial scales, we compared the effects of island isolation and habitat heterogeneity on genetic substructuring in a rocky shore snail ( Littorina saxatilis ). This species has a restricted migration among islands owing to the lack of free-floating larvae. Earlier studies show that allozymes vary extensively as a consequence of isolation by water barriers among islands, but also as a consequence of divergent selection among different microhabitats within islands. In the DNA markers we observed genetic differentiation owing to island isolation at three of nine loci. In addition, variation at three loci correlated with habitat type, but the correlation for two of the loci was weak. Overall, isolation contributed slightly more to the genetic variation among populations than did habitat-related factors but the difference was small. It is concluded that both island isolation, which interrupts gene flow, and a heterogeneous habitat cause genetic substructuring at the DNA level in L. saxatilis in the studied area, and thus in this species we need to be somewhat concerned about habitat heterogeneity also at DNA loci.  © 2004 The Linnean Society of London, Biological Journal of the Linnean Society , 2004, 82 , 377–384.     

Partitioning the effects of isolation by distance,environment, and physical barriers on genomic divergence between parapatric threespine stickleback

Genetic divergence between populations is shaped by a combination of drift, migration, and selection, yielding patterns of isolation‐by‐distance (IBD) and isolation‐by‐environment (IBE). Unfortunately, IBD and IBE may be confounded when comparing divergence across habitat boundaries. For instance, parapatric lake and stream threespine stickleback (Gasterosteus aculeatus) may have diverged due to selection against migrants (IBE), or mere spatial separation (IBD). To quantitatively partition the strength of IBE and IBD, we used recently developed population genetic software (BEDASSLE) to analyze partial genomic data from three lake‐stream clines on Vancouver Island. We find support for IBD within each of three outlet streams (unlike prior studies of lake‐stream stickleback). In addition, we find evidence for IBE (controlling for geographic distance): the genetic effect of habitat is equivalent to geographic separation of ～1.9 km of IBD. Remarkably, of our three lake‐stream pairs, IBE is strongest where migration between habitats is easiest. Such microgeographic genetic divergence would require exceptionally strong divergent selection, which multiple experiments have failed to detect. Instead, we propose that nonrandom dispersal (e.g., habitat choice) contributes to IBE. Supporting this conclusion, we show that the few migrants between habitats are a nonrandom subset of the phenotype distribution of the source population.