MetaPopGen: an r package to simulate population genetics in large size metapopulations

Population genetics simulation models are useful tools to study the effects of demography and environmental factors on genetic variation and genetic differentiation. They allow for studying species and populations with complex life histories, spatial distribution and many other complicating factors that make analytical treatment impracticable. Most simulation models are individual-based: this poses a limitation to simulation of very large populations because of the limits in computer memory and long computation times. To overcome these limitations, we propose an intermediate approach that allows modelling of very complex demographic scenarios, which would be intractable with analytical models, and removes the limitations imposed by large population size, which affect individual-based simulation models. We implement this approach in a software package for the r environment, MetaPopGen. The innovative concept of this approach with respect to the other population genetic simulators is that it focuses on genotype numbers rather than on individuals. Genotype numbers are iterated through time by using random number generators for appropriate probabilistic distributions to reproduce the stochasticity inherent to Mendelian segregation, survival, dispersal and reproduction. Features included in the model are age structure, monoecious and dioecious (or separate sexes) life cycles, mutation, dispersal and selection. The model simulates only one locus at a time. All demographic parameters can be genotype-, sex-, age-, deme- and time-dependent. MetaPopGen is therefore indicated to study large populations and very complex demographic scenarios. We illustrate the capabilities of MetaPopGen by applying it to the case of a marine fish metapopulation in the Mediterranean Sea.