A FORTRAN program to simulate the evolution of genetic variability in a small population

This paper presents a FORTRAN-77 program that performs Monte Carlo simulation of the evolution of genetic structure in a small population under selection. The aim is to study the possibility of foreseeing a response plateau in a theoretical population, depending on population size and management, and to apply this to small populations actually selected, to predict a possible exhaustion of genetic variability. A set of subroutines describing the different steps in a selection cycle (birth, expression of phenotypic value, genetic evaluation, selection, reproduction, death) is available and the user can choose the sequence of subroutines, the characteristics of individuals submitted to each step, and also build more personal subroutines if necessary. The program is based on the generation of exact genotypes and their transmission from parents to offspring, through simulation of meiosis and pairing of gametes. Parameters concerning the genome, the initial structure of the population and its management are required. The genetic mean and variance of the population for each new cycle of selection are given as outputs. Examples of applications are given and discussed.