The Journal of Heredity 2001:92(5)
© 2001 The American Genetic Association 92:447-448
Computer Note |
SGSSpatial Genetic Software: A Computer Program for Analysis of Spatial Genetic and Phenotypic Structures of Individuals and Populations
From SILVOLAB Guyane, INRA Station de Recherches Forestières, Campus agronomique, BP 709, 97387Kourou cedex, France D.O.M. Guyane (French Guiana) (Degen), and INRA Station de Recherches Forestières, Laboratoire de Génétique des Arbres Forestiers, BP 45, F-33611 Gazinet cedex, France (Petit and Kremer).
Address correspondence to B. Degen at the address above or e-mail: degen.b{at}kourou.cirad.fr.
| Abstract |
|---|
|
|
|---|
We have developed a program called Spatial Genetic Software (SGS), which provides a user-friendly Windows tool to analyze both local and broad scale genetic and phenotypic structure. It can deal with nearly any type of genetic data, codominant (allozyme, PCR-RFLP, microsatellite) or dominant (RAPD, AFLP) markers, or biparentally (nuclear) or uniparentally (cpDNA and mtDNA) inherited markers. Data based on any of these markers can be analyzed, either as individual genotypes within a single population (local scale) or as allele or haplotype frequencies from different populations (broad scale). We also include a simple approach to analysis of spatial structure for continuous quantitative traits. The program implements various parameters to analyze spatial genetic and phenotypic structure: Moran's index, Geary's index, number of alleles in common, and approaches using genetic distances and FST values. The statistical significance of all measures is verified by the use of a permutation test. The results are assessed by graphics that can be integrated, via the clipboard, to other Windows programs. The details of the computations are given in a table and can be stored as ASCII files.
| Introduction |
|---|
|
|
|---|
Over the last 20 years, following the seminal articles of Sokal and Oden (1978a, b), the analysis of spatial genetic structure has become an important issue (e.g., Epperson 1992). The majority of the studies, based on allozymes, represent either broad-scale analyses of allele frequencies (e.g., Jones et al. 1980; Sokal and Menozzi 1982) or complete inventories within a population of sessile organisms, such as plants (e.g., Bacilieri et al. 1994; Chung et al. 1998; Epperson and Allard 1989; Leonardi et al. 1996; Merzeau et al. 1994; Waser 1987). Studies using DNA markers are becoming more numerous, both on a cross-population scale (e.g., Dumolin-Lapègue et al. 1997; Petit et al. 1993, 1997) and within a single population (e.g., Streiff et al. 1998; Wagner et al. 1991).
Spatial statistical approaches differ from other statistical approaches that are common in population genetic analyses by emphasizing the interaction of a population or an individual with its neighbors. Local spatial structure is generated after a few generations as a consequence of fine-scale genetic processes, such as limited seed and pollen flow and local selection pressures. Spatial analysis also complements the more traditional methods aimed at detecting processes occurring at larger geographic scales, such as migration and colonization. The analysis of spatial structure enables an estimation of those forcing processes (e.g., Hardy and Vekemans 1999).
Different parameters have been used to quantify spatial genetic structure. Moran's index and Geary's index (Cliff and Ord 1973; Sokal and Oden 1978a) are among the most frequently used measures, but others have used the "number of alleles in common" to measure genetic divergence (Boshier et al. 1995). FST and GST statistics have also been used to quantify genetic structure (Dumolin-Lapègue et al. 1997; Streiff et al. 1998). More recently, multilocus measures of spatial autocorrelation, based on genetic distances, have been introduced (Cassens et al. 2000; Degen and Scholz 1998; Smouse and Peakall 1999; Vendramin et al. 1999).
Other population genetics programs offer complementary but generally limited features for the analysis of spatial genetic structure. Arlequin (Schneider et al. 2000) has a Mantel test procedure, which computes the correlation between genetic and other distance matrices, whereas GENEPOP (Raymond and Rousset 1995) includes a routine for regression of FST or FST/(1 - FST) onto geographic distance, and also performs a Mantel test.
A versatile software for analysis of spatial genetic and phenotypic structure is missing. Spatial Genetic Software (SGS), with its broad set of features for analysis, fills this gap which has been identified repeatedly (e.g., Smouse and Peakall 1999).
| Program Description |
|---|
|
|
|---|
Spatial Distance Classes
All calculated statistics are computed for pairs of data points belonging to a series of spatial distance classes (sq). The spatial distance between two data points is computed as the Euclidean ground distance. The dimension of all distance classes is equal. To check the influence of distance class definition on the results, it is possible to repeat the calculations with different scales for distance classes. The plots of the calculated statistic as a function of interpoint distance are called correlograms for Moran's index and Geary's index or distograms for genetic distance-based measures.
Permutation Testing
A Monte Carlo permutation procedure is applied to test significant deviations from a spatially random distribution of each calculated measure (Manly 1997). Each permutation consists of a random shuffling of genetic or phenotypic data over the spatial coordinates of the sampled points. For each of the spatial distance classes, observed values are compared with a null distribution, obtained from N Monte Carlo trials. Then a user-defined
% confidence interval for the parameters is constructed, by ordering the permuted estimates (e.g., Bacilieri et al. 1994; Streiff et al. 1998).
Analysis of the Spatial Structure
The program works with various types of genetic data, codominant (allozyme, PCR-RFLP, microsatellite) or dominant (RAPD, AFLP) markers, or biparentally (nuclear) or uniparentally (cpDNA and mtDNA) inherited markers. Data based on any of these markers can be analyzed, either as individual genotypes within a single population (local scale) or as allele or haplotype frequencies from different populations (broad scale). We also include a simple approach to analysis of spatial structure for continuous quantitative traits. The program implements various parameters to available measures. An overview is given in Table 1. The user will find details about formulas and the calculations in the manual (see below).
|
| Results |
|---|
|
|
|---|
The software offers various possibilities for viewing and saving the results: distograms or correlograms, and their confidence intervals, based on the permutation test, graphs showing the spatial distribution of data points, or diagrams to view the frequencies of pairs in all spatial distance classes. All graphs of SGS can be imported, via the clipboard, to other Windows applications. Moreover, one can generate tables and files with the results of the spatial analysis and results of the permutation test.
Availability
SGS is written in Visual Basic version 5.0 (Professional Edition) and has been compiled as 32-bit versions for the Microsoft Windows (Windows 95/98/00 and Windows NT) operating system. The program and the user's manual are available on our Internet homepage: http://kourou.cirad.fr/genetique/software.html.
| Acknowledgments |
|---|
We are grateful to Hans-Rolf Gregorius, Giovanni Vendramin, and Birgit Ziegenhagen for valuable discussions and suggestions about the program. We would also like to thank Henri Caron, Cyril Dutech, Brigitte Demesure, and Sylvie Oddou for critical testing of the program and for helpful suggestions. We also wish to thank Peter Smouse and two anonymous reviewers for their very helpful comments and improvements to the manuscript.
| Footnotes |
|---|
Corresponding Editor: Sudhir Kumar
Received April 23, 2001
Accepted June 30, 2001
| References |
|---|
|
|
|---|
-
Bacilieri R, Labbe T, and Kremer A, 1994. Intraspecific genetic structure in a mixed population of Quercus petraea (Matt.) Liebl. and Q. robur L. Heredity 73:130141.
Boshier DH, Chase MR, and Bawa KS, 1995. Population genetics of Cordia alliodora (Boraginaceae), a neotropical tree. 3. Gene flow, neighbourhood and population substructure. Am J Bot 82:484490.
Cassens I, Tiedemann R, Suchentrunk F, and Hartl GB, 2000. Mitochondrial DNA variation in the European otter (Lutra lutra) and the use of spatial autocorrelation analysis in conservation. J Hered 91:3135.
Chung MY, Chung GM, Chung MG, and Epperson B, 1998. Spatial genetic structure in populations of Cymbidium goeringii (Orchidaceae). Genes Genet Syst 73:281285.
Cliff AD and Ord JK, 1973. Spatial autocorrelation. London: Pion Limited.
Degen B and Scholz F, 1998. Spatial genetic differentiation among populations of European beech (Fagus sylvatica L.) in western Germany as identified by geostatistical analysis. For Genet 5:191199.
Deichsel G and Trampisch HJ, 1985. Clusteranalyse und Diskriminanzanalyse. Stuttgart: Gustav Fischer Verlag.
Dumolin-Lapègue S, Demesure B, Fineschi S, Corre VL, and Petit RJ, 1997. Phylogeographic structure of white oaks throughout the European continent. Genetics 146:14751487.[Abstract]
Epperson BK, 1992. Spatial structure of genetic variation within populations of forest trees. New For 6:257278.
Epperson BK and Allard RW, 1989. Spatial autocorrelation analysis of the distribution of genotypes within population of lodgepole pine. Genetics 121:369377.
Gregorius HR, 1978. The concept of genetic diversity and its formal relationship to heterozygosity and genetic distance. Math Biosci 41:253271.[Web of Science]
Hamrick JL, Murawski DA, and Nason JD, 1993. The influence of seed dispersal mechanisms on the genetic structure of tropical tree populations. Vegetation 107:281297.
Hardy OJ and Vekemans X, 1999. Isolation by distance in a continuous population: reconciliation between spatial autocorrelation analysis and population genetics models. Heredity 83:145154.
Jones JS, Selander RK, and Schnell GD, 1980. Patterns of morphological and molecular polymorphisms in the land snail Cepaea nemoralis. Biol J Linn Soc 14:359387.
Leonardi S, Raddi S, and Borghetti M, 1996. Spatial autocorrelation of allozyme traits in a Norway spruce (Picea abies) population. Can J For Res 26:6371.
Manly BFJ, 1997. Randomization, bootstrap and Monte Carlo methods in biology. London: Chapman & Hall.
Merzeau D, Comps B, Thiebaut B, Cuguen J, and Letouzey J, 1994. Genetic structure of natural stands of Fagus sylvatica L. (beech). Heredity 72:269277.
Nei M, 1972. Genetic distance between populations. Am Nat 106:283292.[Web of Science]
Petit RJ, Kremer A, and Wagner DB, 1993. Geographic structure of chloroplast DNA polymorphisms in European oaks. Theor Appl Genet 87:122128.
Petit RJ, Pineau E, Demesure B, Bacilieri R, Ducousso A, and Kremer A, 1997. Chloroplast DNA footprints of postglacial recolonization by oaks. Proc Natl Acad Sci USA 94:999610001.
Pons O and Petit RJ, 1995. Estimation, variance and optimal sampling of gene diversity. I. Haploid locus. Theor Appl Genet 90:462470.
Raymond M and Rousset F, 1995. GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Hered 86:248249.
Schneider S, Roessli D, and Exoffier L, 2000. Arlequin ver. 2.000: a software for population genetics data analysis. Geneva, Switzerland: University of Geneva.
Smouse PE and Peakall R, 1999. Spatial autocorrelation analysis of individual multiallele and multilocus genetic structure. Heredity 82:561573.
Sokal R and Menozzi IP, 1982. Spatial autocorrelation of HLA frequencies in Europe supports demic diffusion of early farmers. Am Nat 119:117.
Sokal RR and Oden NL, 1978a. Spatial autocorrelation in biology. 1. Methodology. Biol J Linn Soc 10:199228.
Sokal RR and Oden NL, 1978b. Spatial autocorrelation in biology. 2. Some biological implications and four applications of evolutionary and ecological interest. Biol J Linn Soc 10:229249.
Sokal RR and Wartenberg DE, 1983. A test of spatial autocorrelation analysis using an isolation-by-distance model. Genetics 105:219237.
Streiff R, Labbe T, Bacilieri R, Steinkellner H, Glossl J, and Kremer A, 1998. Within-population genetic structure in Quercus robur L. and Quercus petraea (Matt.) Liebl. assessed with isozymes and microsatellites. Mol Ecol 7:317328.
Surles SE, Arnold J, Schnabel A, Hamrick JL, and Bongarten BC, 1990. Genetic relatedness in open-pollinated families of two leguminous tree species, Robinia pseudoacacia L. and Gleditsia triacanthos L. Theor Appl Genet 80:4956.
Vendramin GG, Degen B, Petit RJ, Anzidei M, Madaghiele A, and Ziegenhagen B, 1999. High level of variation at Abies alba chloroplast microsatellite loci in Europe. Mol Ecol 8:11171126.[Medline]
Wagner DB, Sun ZX, Govindaraju DR, and Dancik BP, 1991. Spatial patterns of chloroplast DNA and cone morphology variation within populations of a Pinus banksianaPinus contorta sympatric region. Am Nat 138:15670.
Waser NM, 1987. Spatial genetic heterogeneity in a population of the montane perennial plant Delphinium nelsonii. Heredity 58:249256.
This article has been cited by other articles:
![]() |
V. C. R. Azevedo, M. Kanashiro, A. Y. Ciampi, and D. Grattapaglia Genetic Structure and Mating System of Manilkara huberi (Ducke) A. Chev., a Heavily Logged Amazonian Timber Species J. Hered., November 1, 2007; 98(7): 646 - 654. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Reisch, S. Schurm, and P. Poschlod Spatial Genetic Structure and Clonal Diversity in an Alpine Population of Salix herbacea (Salicaceae) Ann. Bot., April 1, 2007; 99(4): 647 - 651. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Mantovani, L. P. C. Morellato, and M. S. dos Reis Internal Genetic Structure and Outcrossing Rate in a Natural Population of Araucaria angustifolia (Bert.) O. Kuntze J. Hered., September 1, 2006; 97(5): 466 - 472. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Y. CHUNG, Y. SUH, J. LOPEZ-PUJOL, J. D. NASON, and M. G. CHUNG Clonal and Fine-scale Genetic Structure in Populations of a Restricted Korean Endemic, Hosta jonesii (Liliaceae) and the Implications for Conservation Ann. Bot., August 1, 2005; 96(2): 279 - 288. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Takahashi PSAwinD Version 1.1.1: A Program for Calculating Spatial Indices J. Hered., May 1, 2003; 94(3): 267 - 270. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

