GENE DISPERSAL AND SPATIAL GENETIC STRUCTURE

Spatial autocorrelation statistics have been studied in theoretical population genetic models and widely used in experimental studies of spatial structure in many plant and animal populations. However, the statistical properties of spatial autocorrelation statistics have remained uncharacterized. Little is known about how values of spatial autocorrelation statistics in population samples depend on the level of dispersal and scheme of sampling. In this paper, we characterize the statistical properties of join-count spatial autocorrelation statistics for population genetic surveys under various conditions of dispersal and sampling. The results indicate generally high statistical power. These results can provide a method to estimate gene dispersal based on standing spatial patterns of genetic variation observed within populations.     

Spatial and population genetic structure of microsatellites in white pine

We evaluated the population genetic structure of seven microsatellite loci for old growth and second growth populations of eastern white pine (Pinus strobus). From each population, located within Hartwick Pines State Park, Grayling, Michigan, USA, 120-122 contiguous trees were sampled for genetic analysis. Within each population, genetic diversity was high and inbreeding low. When comparing these populations, there is a significant, but small (less than 1%), genetic divergence between populations. Spatial distance between populations or timber harvest at the second growth site were reasonable explanations for the observed minor differences in allele frequencies between populations. Spatial autocorrelation analysis suggested that, for the old growth population, weak positive structuring at 15 m fits the isolation by distance model for a neighbourhood size of about 100 individuals. In comparison, genotypes were randomly distributed in the second growth population. Thus, logging may have decreased spatial structuring at the second growth site, suggesting that management practices may be used to alter natural spatial patterns. In addition, the amount of autocorrelation in the old growth population appears to be lower for some of the microsatellites, suggesting higher numbers of rare alleles and that higher mutation rates may have directly affected spatial statistics by reducing structure.     

High levels of allozyme variation within populations and low allozyme divergence within and among species of Hemerocallis (Liliaceae)

Thirty populations from five species of Hemerocallis in Korea were analyzed by starch gel electrophoresis to measure genetic diversity and to determine genetic population structure and the amount of genetic divergence within and between species at 12 isozyme loci. In addition, Moran's I spatial autocorrelation statistics were used to examine the spatial distribution of allozyme polymorphisms in populations of H. thunbergii and H. hakuunensis. Populations of five Korean species maintain high levels of genetic variation and little differentiation among populations and species. Mean expected heterozygosities range from 0.165 in H. hongdoensis, an island endemic, to 0.265 in H. taeanensis, and a total of 81 alleles across the 12 loci were detected in the five species. G(ST) values for each of the five species were low, ranging from 0.051 in H. taeanensis to 0.078 in H. hakuunensis. Mean intraspecific Nei's genetic identities (I) between populations of the five species were all above 0.97. However, a considerable level of heterozygote deficiencies within populations was detected, ranging from 0.242 to 0.411 measured as F(IS) statistics. This deficiencies may be due to inbreeding, limited pollen and seed dispersal, or from the pooling of subpopulations that differ in allele frequencies. A small spatial scale population substructuring (<12 m) was found in H. thunbergii and H. hakuunensis. A group of populations from each of the five previously designated Hemerocallis species (based on their morphology, ecology, and phenology) agrees with our allozyme data, though pairwise comparisons among species had high I values (from 0.862, H. middendorffii vs. H. hongdoensis, to 0.969, H. thunbergii vs. H. taeanensis). This is attributed to the presence of the same high-frequency alleles in different species at seven loci. In addition, no "diagnostic allele" that appears in all populations of one species, but is absent in other species, was detected at the 12 isozyme loci. These all suggest that species of Hemerocallis in Korea may have recently derived from an ancestor or progenitor harboring high levels of genetic diversity.     

Estimating dispersal from short distance spatial autocorrelation

A series of theoretical studies has formed a strong connection between spatial statistics observed in populations and summary measures of the amount of dispersal. Synthesized, these developments allow dispersal to be indirectly estimated from standing spatial patterns of genetic variation under a range of conditions broad enough to be likely met in most populations of either plants or animals. The spatial correlations at the shortest distances are particularly robust to range of conditions and have disproportionately high statistical power. This review integrates theoretical results in a way that maximizes robustness and flexibility in the use of short distance autocorrelation to estimate Wright's neighborhood size, or the total variance in dispersal distances. Empirical guidelines are developed that are meant to be as practical and broad as possible. The guidelines focus on Moran's I-statistics for diploid genotypes converted to allele frequencies, but are also extended to or compared with several other approaches.     

On methods of spatial analysis for genotyped individuals

Spatial autocorrelation methods have commonly been applied to individual-based spatial genetic studies, although their properties and the relations among the statistics have not been carefully examined. This paper first introduces a reformulation of widely used spatial statistics using point processes. When Moran's I statistics are applied to allele frequencies within an individual, the frequencies are no longer continuous variables but have only three discrete values and specific interpretations of Moran's I statistics and the number of alleles in common (NAC) can be expressed as the weighted sum of join-count statistics. The distributions of minor genotypes are amplified in Moran's I depending on the allele frequency in the population, while NAC uses a constant weighting system. Under the point process framework, spatial analysis can be conducted on the common theoretical base, from individual locations to genetic distributions of different levels, (for example, genotype and allele). The methodology is demonstrated by application to field data for molecular ecological studies of Fagus crenata population dynamics.     

Multilocus estimation of genetic structure within populations

Spatial structure of genetic variation within populations is well measured by statistics based on the distribution of pairs of individual genotypes, and various such statistics have been widely used in experimental studies. However, the problem of uncharacterized correlations among statistics for different alleles has limited the applications of multiallelic, multilocus summary measures, since these had unknown sampling distributions. Usually multiple alleles and/or multiple loci are required in order to precisely measure spatial structures, and to provide precise indirect estimates of the amount of dispersal in samples of reasonable size. This article examines the correlations among pair-wise statistics, including Moran I-statistics and various measures of conditional kinship, for different alleles of a locus. First the correlations are mathematically derived for random spatial distributions, which allow averages over alleles and loci to be used as more powerful yet exact test statistics for the null hypothesis. Then extensive computer simulations are conducted to examine the correlations among values for different alleles under isolation by distance processes. For loci with more than three alleles, the results show that the correlations are remarkably and perhaps surprisingly small, establishing the principle that then alleles behave as nearly independent realizations of space-time stochastic processes. The results also show that the correlations are largely robust with respect to the degree of spatial structure, and they can be used in a straightforward manner to form confidence intervals for averages. The results allow a precise connection between observations in experimental studies and levels of dispersal in theoretical models.     

Spatial genetic structure in wild cherry (Prunus avium L.): II. Effect of density and clonal propagation on spatial genetic structure based on simulation studies

We conducted simulations to disentangle the effect of density and clonal propagation on spatial genetic structure and genetic diversity parameters in Prunus avium. In a previous paper, we observed stronger family structure in populations exhibiting high density and high clonal propagation, whereas one low-density and low clonal propagation population showed negative spatial autocorrelation. We tried to understand these results by simulating 200 years of growth, mating and dispersal with the model Eco-Gene, using two levels of density and three levels of clonal propagation in one high-density population (Spargründe) and one low-density population (Chorin). We used allele frequencies from eight microsatellite loci to generate the populations used for simulations. In order to detect effects of the initial structure on the results, we also ran the simulations starting from the real data. We observed positive effect of clonal propagation on the strength of SGS, while high densities exhibited a negative impact. Genetic diversity was maintained at high densities, while pollen dispersal was shorter. Heterozygosity increased at higher clonal propagation rates, although genetic diversity (accumulated number of genotypes) was lower. We discuss the results according to mating processes and potential management scenario.     

THE GENETIC STRUCTURE OF LARGE-LAKE DAPHNIA POPULATIONS

Cyclical parthenogenesis allows study of the genetic and evolutionary characteristics of groups exhibiting both asexual and sexual reproduction. The cladoceran genus Daphnia contains species which vary with respect to the relative incidence of sexual reproduction; pond species tend to undergo sexual reproduction more regularly than species found in large lakes. Previous genetic studies have focused on pond populations, generating expectations about large-lake populations that have not been fully met by recent studies. The present study of the Palearctic species Daphnia galeata further examines the genetic structure of large-lake populations. Nine local populations, from lakes in northern Germany, are examined for genetic variation at seven enzyme loci. Populations exhibit similar allelic arrays and often similar allele frequencies at the five polymorphic loci; values of Nei's genetic distance (D) ranged from 0.002 to 0.239, with a mean of 0.084. F_ST values range from 0.012 to 0.257, and spatial autocorrelation coefficients range from -0.533 to 0.551, for the eight alleles analyzed. With few exceptions, within-population genotypic frequencies were in Hardy-Weinberg equilibrium. There was, however, significant heterogeneity in genotypic frequencies among populations. The number of coexisting clonal groups, as determined by three locus genotypes, is high within populations. Clonal groups are widely distributed among localities. The amount of genetic divergence observed among these large-lake populations is smaller than that previously observed among pond populations and suggests that different processes may be important in determining the genetic structure and subsequent phenotypic divergence of lake versus pond populations.     

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

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

DISPERSAL AND GENETIC STRUCTURE IN KANGAROO RATS

We used spatial autocorrelation of allele frequencies to examine local structure in a population of bannertailed kangaroo rats for which Wright's isolation-by-distance model seems applicable, and for which we can estimate neighborhood size based on 10 years of data on demography and dispersal. The uniform dispersion and strong philopatric tendencies of this species provide a test case for the idea that restricted dispersal can lead to local genetic structure in small mammals. Whether we considered such complications as nonnormal dispersal distances, variation in lifetime reproductive success, fluctuating population density, and adult as well as juvenile dispersal, our estimate of effective population size was fewer than 15 animals. Nevertheless, data from four polymorphic allozyme loci analyzed over a range of separations between 50 m (approximately one home range diameter) and 1,000 m detected no evidence for spatial clustering of alleles. One resolution of this apparent paradox is that “gamete dispersal,” caused by the movements of males away from their residences during the breeding season, may be a significant (and unmeasured) component of gene dispersal. Our analyses also demonstrate that a decline in population density may actually increase neighborhood size. A more general implication is that even extremely philopatric mammals have effective population sizes large enough to prevent the development of local genetic structure.     

Sampling for Microsatellite-Based Population Genetic Studies: 25 to 30 Individuals per Population Is Enough to Accurately Estimate Allele Frequencies

One of the most common questions asked before starting a new population genetic study using microsatellite allele frequencies is “how many individuals do I need to sample from each population?” This question has previously been answered by addressing how many individuals are needed to detect all of the alleles present in a population (i.e. rarefaction based analyses). However, we argue that obtaining accurate allele frequencies and accurate estimates of diversity are much more important than detecting all of the alleles, given that very rare alleles (i.e. new mutations) are not very informative for assessing genetic diversity within a population or genetic structure among populations. Here we present a comparison of allele frequencies, expected heterozygosities and genetic distances between real and simulated populations by randomly subsampling 5–100 individuals from four empirical microsatellite genotype datasets (Formica lugubris, Sciurus vulgaris, Thalassarche melanophris, and Himantopus novaezelandia) to create 100 replicate datasets at each sample size. Despite differences in taxon (two birds, one mammal, one insect), population size, number of loci and polymorphism across loci, the degree of differences between simulated and empirical dataset allele frequencies, expected heterozygosities and pairwise F_ST values were almost identical among the four datasets at each sample size. Variability in allele frequency and expected heterozygosity among replicates decreased with increasing sample size, but these decreases were minimal above sample sizes of 25 to 30. Therefore, there appears to be little benefit in sampling more than 25 to 30 individuals per population for population genetic studies based on microsatellite allele frequencies.     

Fine-scale genetic structure within plots of Polygala reinii (Polygalaceae) having an ant-dispersal seed

Fine-scale genetic structure within a population was analyzed for the myrmecochorous forest perennial Polygala reinii (Polygalaceae) using allozyme loci. In the analysis, two sampling plots were established to cover the isolated patchy distribution within the study population. Size and spatial structure were also examined for the plots to assess their interaction with the genetic structuring. Allozyme analysis based on 13 putative loci encoding 10 enzyme systems showed high genetic variation and low values of fixation indices at the two plots. Spatial autocorrelation analysis based on the multilocus coancestry coefficient (f _ij ) revealed significant genetic structuring in both plots, suggesting limited gene-, especially seed-, dispersal within the population. The spatial structure within the plots, assessed by O-ring statistics, was characterized by the occurrence of spatial clustering of individuals within a few meters. In particular, the range of the spatial clustering roughly corresponded to that of the genetic structuring. While the size structure did not significantly differ between the plots, these results indicate that the fine-scale genetic structure reflects the formation of spatial clustering of related individuals within the population, partly due to the restricted ant-mediated seed dispersal in P. reinii.     

Spatial genetic structure of two sympatric neotropical palms with contrasting life histories

The spatial genetic structure within sympatric populations of two neotropical dioecious palm species with contrasting life histories was characterized to evaluate the influence of life history traits on the extent of genetic isolation by distance. Chamaedorea tepejilote is a common wind-pollinated arboreal understory palm. Chamaedorea elatior is an uncommon climbing subcanopy palm with entomophilous pollination syndrome. A total of 59 allozyme alleles for C. tepejilote and 53 alleles for C. elatior was analyzed using both unweighted (Iu) and weighted (Iw) Moran's I spatial autocorrelation statistics. The spatial genetic structure detected within these populations is consistent with those reported for highly dispersed plants. A significance test for differences between mean Moran's I-coefficients revealed less spatial genetic structure within the C. tepejilote population than that in the C. elatior population. Adjacent individuals of C. elatior exhibited significant spatial genetic autocorrelation (Iu=0.039, Iw=0.034), indicating a Wright's neighborhood size of about 100 individuals. For C. tepejilote, nonrandom genetic distribution among nearest neighbors was detected, even from small spatial autocorrelation values (Iu=0.008, Iw=0.009), consistent with a neighborhood size of about 300 individuals. For both species, seed dispersal, mortality among life cycle stages, overlapping generations, and contrasting traits of mating and reproduction influence the standing spatial genetic structure within populations.     

Genetic structure of age classes in Camellia japonica (Theaceae)

Camellia japonica L. (Theaceae), an insect- and bird-pollinated, broad-leaved evergreen tree, is widely distributed in Japan and the southern Korean peninsula. The species has a relatively even age distribution within populations, which may influence the spatial genetic structure of different age classes relative to species with typical L-shaped age distributions. To determine whether the internal spatial genetic structure found in seedlings and young individuals carries over into adults, we used allozyme loci, F-statistics, spatial autocorrelation statistics (Moran's I), and coancestry measures to examine changes in genetic structure among seven age classes in a population (60-m x 100-m area) in southern Korea. In seedlings, weak but significant positive values of Moran's I-statistics and coancestry measures were found for distances less than 14 m, which is consistent with a mechanism of limited seed dispersal combined with overlapping seed shadows. This spatial structure, however, dissipates in older age classes, and in adults genetic variation has an essentially random spatial distribution. Morisita's index of dispersion of individuals in each age class showed that seedlings and juveniles are more highly clustered than are older individuals. These results suggest that self-thinning changes the spatial relationships of individuals, and thus genotypes. A multilocus estimate of FST (0.008) shows a small but statistically significant difference in allele frequencies among age classes. In summary, intrapopulation genetic structure within and among age classes of C. japonica was significant but weak. Despite presumably limited seed dispersal, weak spatial genetic structure in juveniles suggests overlapping seed shadows followed by self-thinning during recruitment. The present study also demonstrates that studies of spatial genetic structure focusing on limited numbers of generations may not be sufficient to reveal the entire picture of genetic structure in populations with overlapping generations.     

Genetic neighbourhood of clone structures in eelgrass meadows quantified by spatial autocorrelation of microsatellite markers

Limited dispersal distances in plant populations frequently cause local genetic structure, which can be quantified by spatial autocorrelation. In clonal plants, three levels of spatial organization can contribute to positive autocorrelation; namely, the neighbourhood of (a) ramets, (b) clone fragments and (c) entire clones. Here we use data from an exhaustive sampling scheme on a clonal plant to measure the contribution of the neighbourhoods of each distinct clonal structure to total spatial autocorrelation. Four plots (256 grid points each) within dense meadows of the marine clonal plant Zostera marina (eelgrass) were sampled for clone structure with nine microsatellite markers ( approximately 80 alleles). We found significant coancestry (f(ij)), at all three levels of spatial organization, even when not allowing for joins between samples of identical genets. In addition, absolute values of f(ij) and the maximum distance with significant positive f(ij) decreased with the progressive exclusion of joins between alike genotypes. The neighbourhood of this clonal plant thus consists of three levels of organization, which are reflected in different kinship structures. Each of these kinship structures may affect the level of biparental inbreeding and the physical distance between flowering shoots and their outcrossing neighbourhood. These results also emphasize the notion that spatial autocorrelation crucially depends on the scale and intensity of sampling.     

用分子标记揭示植物随机大居群中亚居群的遗传结构——茅栗自然居群空间遗传结构的SSR分析

采用微卫星标记对茅栗(Castanea sequinii)随机大居群以及其中各亚居群的遗传结构进行了空间自相关分析,以探讨植物自然居群内遗传变异的分布特征及其形成机制。通过9对微卫星引物所产生的119个多态位点,测定了大别山区域内茅栗居群以及各亚居群的空间自相关系数Moran's I值。结果表明:大别山分布的野生茅栗为一个缺乏空间结构的随机大居群,茅栗亚居群之间频繁的花粉流削弱了地理隔离导致的遗传漂变或分化作用在维系居群随机遗传结构中具有的重要作用。但是,在接近亚居群大小的地域范围内(0.228 km)具有一定的空间结构,即小地域尺度中的亚居群存在着空间遗传结构。取样的3个亚居群在小格局范围内都存在一定的空间结构,遗传变异基本上呈非随机分布,在短距离内(61 m)3个亚居群一致表现出不同程度的显著正相关,而随着距离的增加,Moran's I值虽然在不同亚居群间存在一定差异变化,但是总体而言趋向预期值,即不存在空间结构,说明其遗传变异在亚居群内只是在短距离内形成一定的空间结构。研究认为有限的种子散播以及微生境选择等因素可能是产生这种小格局的遗传结构的主要原因。上述研究结果有助于进一步了解植物随机大居群的进化历史和生态过程,同时也为栗属植物中国特有种的保育策略提供了科学依据。     

Methods for estimating gene frequencies and detecting selection in bacterial populations

Recent breakthroughs in molecular technology, most significantly the polymerase chain reaction (PCR) and in situ hybridization, have allowed the detection of genetic variation in bacterial communities without prior cultivation. These methods often produce data in the form of the presence or absence of alleles or genotypes, however, rather than counts of alleles. Using relative allele frequencies from presence-absence data as estimates of population allele frequencies tends to underestimate the frequencies of common alleles and overestimate those of rare ones, potentially biasing the results of a test of neutrality in favor of balancing selection. In this study, a maximum-likelihood estimator (MLE) of bacterial allele frequencies designed for use with presence-absence data is derived using an explicit stochastic model of the host infection (or bacterial sampling) process. The performance of the MLE is evaluated using computer simulation and a method is presented for evaluating the fit of estimated allele frequencies to the neutral infinite alleles model (IAM). The methods are applied to estimate allele frequencies at two outer surface protein loci (ospA and ospC) of the Lyme disease spirochete, Borrelia burgdorferi, infecting local populations of deer ticks (Ixodes scapularis) and to test the fit to a neutral IAM.     

Fine-scale genetic structure of grape phylloxera from the roots and leaves of Vitis

Patterns of variation at microsatellite loci suggest that root populations of the pest grape phylloxera (Daktulosphaira vitifoliae) are largely parthenogenetic in Australian vineyards. To investigate reproduction in leaf galling phylloxera and the association between these individuals and phylloxera on roots, we examined in detail genetic variation in phylloxera from a vineyard block. Some genotypes found on leaf galls within this block were not present on roots, whereas others spanned both zones. There was no evidence that genotypes on roots were the product of sexual reproduction in leaf galls. mtDNA variation was not associated with the location of the phylloxera clones. The spatial distribution of genotypes within a root population was further investigated by intensively sampling phylloxera from another vineyard block. Join-count spatial autocorrelation statistics were used to explore fine-scale spatial structure. Clones were nonrandomly distributed within the block and there was evidence that the distribution of clones followed rows. These findings suggest firstly that there is limited dispersal of root and leaf feeding phylloxera, and secondly that factors, other than vine host, are likely to be important and contribute to clonal structure within populations.     

THE EVOLVING GENETIC HISTORY OF A POPULATION OF LATHYRUS SYLVESTRIS: EVIDENCE FROM TEMPORAL AND SPATIAL GENETIC STRUCTURE

We analyze patterns of genetic microdifferentiation within a natural population of Lathyrus sylvestris, a perennial herb with both sexual reproduction and clonal growth. In a population from the northern foothills of the Pyrénées in southwestern France, a combined demographic and genetic investigation enabled the study not only of spatial genetic structure of the population, but also of the history of the population's spatial genetic structure over time. Excavation of all individuals allowed identification of clonemates. Age of each individual was determined by counting annual growth rings in the taproot, a method tested with individuals of known age planted in experimental gardens. Each individual was mapped, and genotypes of all individuals were determined using allozyme markers for a number of polymorphic loci. Distribution patterns and spatial genetic structure, both for all individuals and for different age classes, were analyzed using spatial autocorrelation statistics (Geary's Index, Moran's Index). Patterns of gene flow within the population were also studied using F-statistics and tests for random associations of alleles. Because age, allele frequencies, and location were known for each individual, it was possible to study how spatial genetic structure changed over time. Results from all these diverse approaches are consistent with one another, and clearly show the following: (1) founder effects, with the study transect being first colonized by individuals at either end of the transect that were homozygous for different alleles at one marker locus; (2) a difference in spatial distribution of individuals originated from sexual reproduction (seedlings) and from clonal growth (connected individuals); (3) restricted gene flow, due to inbreeding among related, clumped individuals; and (4) increase in heterozygote deficit within the youngest cohort of individuals. The results indicate that genetic differentiation in time was much less marked than differentiation in space. Nevertheless, the results revealed that the studied population is experiencing demographic and genetic variation in time, suggesting that it is not at equilibrium. On the one hand, spatial structuring is becoming less marked due to the recombination of founder genotypes; on the other hand, as establishment of new individuals increases, a new spatial structure emerges due to mating between relatives.