The Journal of Heredity 2001:92(3)
© 2001 The American Genetic Association 92:301-302
Computer Note |
EASYPOP (Version 1.7): A Computer Program for Population Genetics Simulations
From the Institute of Zoology, University of Bern, CH3032, Hinterkappelen-Bern, Switzerland.
Address correspondence to Francois Balloux at the address above or e-mail: francois. balloux{at}esh.unibe.ch
 |
Introduction |
|---|
 Top Introduction Program DescriptionInput and OutputHow to get the...References |
|---|
If there are many programs for population genetics data analysis, less effort has been devoted to simulation software in this field. However, there are many situations where simulated population genetics datasets would be very useful. Simulations allow exploring situations too complex to be solved analytically. Indeed, outside very simple population genetics models, it is often extremely difficult to obtain the expectations for quantities such as gene correlation or effective size. Simulations also permit us to validate inferences drawn from empirical studies. They can be run with different estimated parameters and compared to the original dataset. It is also possible to build any null hypotheses against which inferences drawn from real data can be tested. Another nice feature of simulated datasets is the possibility of testing new software devoted to population genetics analysis under nontrivial conditions. In addition to their value in research, simulated data can be very useful for educational purposes. Various datasets can be generated and provided to students to train their skills in population genetics data analysis. All these applications prompted me to develop the software EASYPOP, an individual-based model that simulates neutral loci datasets under a very broad range of conditions.
|
Program Description |
|---|
|
TopIntroductionProgram DescriptionInput and OutputHow to get the...References |
|---|
EASYPOP can simulate haploid, diploid, or haplodiploid organisms. For haploid organisms, the proportion of clonal and sexual reproduction (recombination) can be set, whereas for diploid organisms with one sex, the proportion of selfing can also be selected. For diploid organisms with two sexes, there is the choice between hermaphrodites or sexuals. For hermaphrodites, the proportion of selfing can be set, whereas for sexuals complex breeding structures can be simulated (e.g., monogamy with a given proportion of extra-pair matings). Numbers of individuals (males and females) can be set independently for each population and dispersal is sex specific. Various migration models are implemented, such as the two-dimensional stepping stone or hierarchical island model. In addition, there is an isolation by distance option that works with the coordinates of the populations in any number of dimensions. The migration model and migration rates can be modified after a certain number of generations. This later option allows us to explore the genetic consequences of a modification of the migration pattern or to impose a bottleneck on a population. There are also several mutation models implemented, which are oriented toward the simulation of microsatellite loci. Mutation rates, mutation models, and the maximal number of allelic states can be set individually for each locus. Genotypes are real multilocus; that is, genotypes are not independent replicates over loci.
EASYPOP can handle very large simulations on standard personal computers. Most parameters are limited only by the memory of the computer. However, the number of possible allelic states and the number of populations are limited to 999 and 10,000, respectively. It is possible to run up to 99 independent replicates. The program works with probabilities (e.g., migration, mutation), and thus allows fitting to analytical variance expectations.
EASYPOP has been debugged as carefully as possible by testing simulation results against analytical expectations in all cases where they could be computed. Some of these tests can be found in Balloux et al. (2000).
|
Input and Output |
|---|
|
TopIntroductionProgram DescriptionInput and OutputHow to get the...References |
|---|
The input consists only of the parameters of the simulation chosen by the user. The outputs of the program are the multilocus genotypes after a given number of generations. Output files fitting FSTAT (Goudet 1995), GENEPOP (Raymond and Rousset 1995) or ARLEQUIN (Schneider et al. 1997) are generated. This allows for direct analysis of the results using the mentioned population genetics packages. An additional output file giving the number of alleles, Nei and Chesser's (1983) unbiased heterozygosities, and F statistics for each 10th generation is provided. This file can be used to follow these quantities over time.
|
How to get the Program |
|---|
|
TopIntroductionProgram DescriptionInput and OutputHow to get the...References |
|---|
EASYPOP is available in Mac PPC and Windows 95/98/NT4 versions. The program can be downloaded free of charge with its user manual at the following website: http://www.unil.ch/izea/softwares/easypop.html. After each substantial development of EASYPOP, people who register the program will be contacted about updates. The C source code can be obtained by e-mailing: francois.balloux;caesh.unibe.ch.
|
Acknowledgments |
|---|
It is a pleasure to thank Jean Lehman who introduced me to programming. I am also deeply indebted to Max Reuter who showed me how to program in C, and helped me numerous times to find out some vicious bugs. I also thank Jerôme Goudet, Nicolas Perrin, Eric Petit, Thierry De Meeûs, Pierre Fontanillas, Laurent Lehman, Patrick Basset, and Renaud Kaeuffer for the time they invested to test this program and for comments and suggestions.
|
Footnotes |
|---|
Corresponding Editor: Sudhir Kumar
Received June 14, 2000
Accepted October 31, 2000
|
References |
|---|
|
TopIntroductionProgram DescriptionInput and OutputHow to get the...References |
|---|
-
Balloux F, Brünner H, Lugon-Moulin N, Hausser J, and Goudet J, 2000. Microsatellites can be misleading: an empirical and simulation study. Evolution 54:1414–1422.[ISI][Medline]
Goudet J, 1995. FSTAT: A computer program to calculate F-statistics. J Hered 86:485–486.
Nei M and Chesser RK, 1983. Estimation of fixation indices and gene diversities. Ann Hum Genet 47:253–259.[ISI][Medline]
Raymond M and Rousset F, 1995. GENEPOP: population genetics software for exact tests and ecumenicism. J Hered 86:248–249.
Schneider S, Kueffer JM, Roessli D, and Excoffier L, 1997. ARLEQUIN (version 1.1): a software for population genetic data analysis. Geneva, Switzerland: Genetics and Biometry Laboratory, University of Geneva.
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  B. Padhukasahasram, P. Marjoram, J. D. Wall, C. D. Bustamante, and M. Nordborg Exploring Population Genetic Models With Recombination Using Efficient Forward-Time Simulations Genetics, April 1, 2008; 178(4): 2417 - 2427. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 T. Antao, A. Beja-Pereira, and G. Luikart MODELER4SIMCOAL2: A user-friendly, extensible modeler of demography and linked loci for coalescent simulations Bioinformatics, July 15, 2007; 23(14): 1848 - 1850. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
F. Guillaume and J. Rougemont Nemo: an evolutionary and population genetics programming framework Bioinformatics, October 15, 2006; 22(20): 2556 - 2557. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
M. Foll and O. Gaggiotti Identifying the Environmental Factors That Determine the Genetic Structure of Populations Genetics, October 1, 2006; 174(2): 875 - 891. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
A. Vasemagi The Adaptive Hypothesis of Clinal Variation Revisited: Single-Locus Clines as a Result of Spatially Restricted Gene Flow Genetics, August 1, 2006; 173(4): 2411 - 2414. [Full Text] [PDF]
|
|
|
|
|
|
J. Goudet and L. Buchi The Effects of Dominance, Regular Inbreeding and Sampling Design on QST, an Estimator of Population Differentiation for Quantitative Traits Genetics, February 1, 2006; 172(2): 1337 - 1347. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
B. Peng and M. Kimmel simuPOP: a forward-time population genetics simulation environment Bioinformatics, September 15, 2005; 21(18): 3686 - 3687. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
G. Guillot, A. Estoup, F. Mortier, and J. F. Cosson A Spatial Statistical Model for Landscape Genetics Genetics, July 1, 2005; 170(3): 1261 - 1280. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
F. Balloux, L. Lehmann, and T. de Meeus The Population Genetics of Clonal and Partially Clonal Diploids Genetics, August 1, 2003; 164(4): 1635 - 1644. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
O. J. Hardy, N. Charbonnel, H. Freville, and M. Heuertz Microsatellite Allele Sizes: A Simple Test to Assess Their Significance on Genetic Differentiation Genetics, April 1, 2003; 163(4): 1467 - 1482. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 M. C. Fisher, B. Rannala, V. Chaturvedi, and J. W. Taylor Disease surveillance in recombining pathogens: Multilocus genotypes identify sources of human Coccidioides infections PNAS, June 25, 2002; 99(13): 9067 - 9071. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||