Gene-Net

a program to simulate dispersal and gene flow in synthetic river networks


page updated: 05/01/2008


Gene-Net  v 1.0
 
 

Introduction

Model principles

References

Installation and Downloads

Input files

User manual

Contacts

 

Contributors by alphabetic order: E. d'Arlhac, A. Chaput-Bardy, C. Fleurant, C. Lemaire, J.M. Richer, J. Secondi.

Introduction

Gene-Net software allows to test hypotheses on genetic structure of freshwater organisms related to river network geometry. Users start to select synthetic river network and initial population parameters. Configuration of the network is characterized by the length and the number of branches. Each individual of a population located on the network is characterised by its sex, a life cycle, its genotype and a probability to move taking into account the population density. Then, the program simulates dispersal and breeding along the river network. Genetic data resulting from the simulation are available for any located population at any time in ARLEQUIN format.
 
 


Model principles

Network generation module:
Demographic simulations are performed along a synthetic network discretized into evenly spaced nodes. Each node is connected to one or two upper nodes (upstream) and one lower node (downstream) to allow migration in both directions. Nodes and river network characteristics are computed from a fractal river network generation module. This module needs four parameters that are enough to describe a fractal river network: two Horton’s ratios RB and RL, the maximum order of the river network n and a scale factor. The scale factor is used to multiply the space-step (distance between two nodes) and thus generate real size networks. The precise position of a stream segment within the whole river network can be identified by classifying the network. The classification system is as follows (Figure 1):
• headwaters are first order stream segments;
• when two stream segments within the same order w = i merge, the stream segment resulting from this confluence is within order  w = i+1;
• when two stream segments within different orders, w = i and  w = j merge, the stream segment resulting from this confluence is within order  w = max(i,j).
This classification can be put into mathematical equation and sum-up by equation (1).
        (1)

where d is the Kronecker index :


Figure 1 : Stralher Classification.

Strahler's classification can organize the different segments of a stream network into a hierarchy. Consequently, the stream outlet will have the highest index value, corresponding to the river network order. This classification puts forward general geometric laws. Among them, Horton's laws (1945) describe the way stream networks are organized (Figure 1). These laws express the so-called bifurcation ratio RB and length ratio RL, also known as Horton's ratios. A great number of experimental studies on stream networks (e.g.  La Barbera, Rosso, 1990; Tarboton et al., 1990; Rosso et al., 1991) revealed that these ratios are rather stable and fluctuate between 3 and 5 for RB and between 1.5 and 3.5 for RL. The calculation of the fractal dimension of a stream network varies according to the observation scale; the stream network is therefore a multifractal object (Rodriguez-Iturbe, Rinaldo, 1997). Horton's laws also make it possible to work out the RB and RL (equation 2).
      (2)
Where  is the average value of the morphometric lengths of w order and   is the number of morphometric lengths of w order.
Then the fractal dimension D can be computed by (equation 3).
        (3)

Population module:
Populations located on the river network are characterized by the life cycle of the study species. Life cycles are designed for freshwater organisms restricted to watercourses, with or without larval stages. Equation (4) was used for a damselfly species, Calopteryx splendens (Harris, 1782). Life cycle of odonates is composed by egg, larval, immature adult (young) and breeder stages. Each stage is characterized by a survival rate which allows calculation of the population growth rate l.
        (4)
Where Segg, Slarvae, Syoung and Sbreeder are respectively the survival rates of eggs, larvae, young adults and breeders. f, the fertility rate represents the number of eggs laid by a female in average.
A population consists in individuals whose sex and alleles are randomly drawn. Allele frequencies follow the Dirichlet model. Furthermore, Hardy Weinberg equilibrium is assumed in each population. The implementation of the life cycle is applied sequentially by the reproduction of individuals within a population and survival rates for different stages. When reproduction occurs, females are fertilized by at least one male. An egg results in the random sampling of one allele from each parent. One or more populations are located on selected places on the network. For a given population, the displacement of individuals takes place according to a migration rate parameter. Displacement from a node to another is governed by a uniform probability density function which allows upstream and downstream migrations. This probability density function gives the individual ability to move along some distance. This travel distance is assigned to an individual at birth and remains constant throughout the life cycle. The spatial migration process leads to an increase in the size of populations that were originally under the carrying capacity threshold. Boundary nodes (boundaries of the river catchment) are linked to upper nil nodes. For these nodes special conditions are applied. These are Dirichlet’s boundary conditions where the number of individuals remains constant when the density threshold is reached.

Simulation module:
The simulated demographic and genetic information are stored into output text files for later use with any visualization data software. Information available are  and  values and the regression slope between geographic and genetic distances, isolation by distance (Rousset, 1997) of any population, at any time and at any location. The figure 2 shows an isolation by distance pattern obtained from the species Calopteryx splendens in a river network whose parameters are RB = 3 and RL = 2 and n = 3 after 100 generations. Gene-Net outputs files containing sampled genetic diversity in the ARLEQUIN format (Excoffier et al., 2005).


Figure 2 : Differentiation among aquatic insects. F-Statistics are plotted against distance. River networks parameters are RB= 3 and RL = 2 and n = 3.

References

Campbell Grant EH, Lowe WH, Fagan WF (2007) Living in the branches: population dynamics and ecological processes in dendritic networks. Ecology Letters 10, 165-175.
Excoffier L, Laval G, S. S (2005) Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evolutionary Bioinformatics Online 1, 47-50.
Gibbs HL, Gibbs KL, Siebenmann M, Collins L (1998) Genetic differentiation among populations of the rare mayfly Siphlonisca aerodromia Needham. J. N. Am. Benthol. Soc. 17, 461-474.
La Barbera P, Rosso R (1990) On the fractal dimension of stream networks. Water Resources Research 26, 2245-2248.
Lowe WH, Likens GE, Power ME (2006) Linking scales in stream ecology. BioScience 56, 591-597.
Neuenschwanger S (2006) AQUASPLATCHE: a program to simulate genetic diversity in populations living in linear habitats. Molecular Ecology Notes 6, 583-585.
Power ME, Dietrich WE (2002) Food webs in river networks. Ecological Research 17, 451-471.
Rodriguez-Iturbe I, Rinaldo A (1997) Fractal river basins Cambridge Unviversity Press.
Rosso R, Bacchi B, La Barbera P (1991) Fractal relation of mainstream length to catchment area in river networks. Water Resources Research 27, 381-387.
Rousset F (1997) Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics 145, 1219-1228.
Schmidt SK, Hughes JM, Bunn SE (1995) Gene flow among conspecific populations of Baetis spp. (Ephemeroptera): adult flight and larval drift. J. N. Am. Benthol. Soc. 14, 147-157.
Tarboton D, Bras R, Rodriguez-Iturbe I (1990) Comment on "on the fractal dimension of stream networks" by La Barbera and Rosso. Water Resources Research 26, 2243-2244.
 
 

Installation and Downloads

Two Windows versions of Gene-Net are available, a GUI (Graphical User Interface) version and a consol version (see below).

Installation in GUI version:

    1. To run Gene-Net in GUI version you need a Java Runtime Machine (JRM) installed on your PC. Most of the time JRM is installed by default. If you do not have a JRM, download the JRM on the site http://www.java.com/fr/download/.

    2. Then download the archive genenet.zip which contains two applications: the graphical user interface (genenet.jar) and the computing program (genenet.exe). Unzip the archive into a directory (D:/ is for instance).

    3. Double click genenet.jar to launch the GUI (the program genenet.exe is called by the GUI).

    4. It is necessary to configure the GUI. Go to "Settings". In "Gene-Net directory" choose the directory where genenet.jar and genenet.exe are located (if you have unziped the archive in D:/, then genenet.jar and genenet.exe are in D:/genenet/. Then in "Gene-Net binary" selecte the program genenet.exe (if you have unziped the archive in D:/, then genenet.exe is in D:/genenet/

    5. Click "Start". The Settings window opens. With the two top buttons "Clear all parameters" and "Default Parameters" you can respectively delete all parameters value and give default values. Click "Start" again to run the simulation.
 

Installation in Consol version:

    1. Download the archive genenetconsol.zipwhich contains two files: the computing program (genenet.exe) and the inputfile of the parameters (parameters.input). Unzip the archive into a directory (D:/ for instance).

    2. In a script consol (e.g. DOS) type the command : genenet.exe parameters.input
 
 

Input files

Here are some input files (parameters.inputs) which could be used in consol version to test several scenarios. These input files had been tested on a PC intel(R) 4, CPU 3.20 GHz and 2.00 Go RAM:
 
 
  Description of the scenario
Download the input file
CPU time
Download Results (Graphs)
INSECT 
- River network : RB=3, RL=2, n=3, scale=1
- Initial population : 200 individuals
- Population: damselfly species, Calopteryx splendens
- Simulation: 100 generations
parameters.inputs
3 s
<Fis> vs distance to outlet
<Fis>,<Fit>,<Fst> vs generations
"Genetic distance" vs distance
INSECT
- River network : RB=4, RL=2, n=3, scale=1
- Initial population : 200 individuals
- Population: damselfly species, Calopteryx splendens
- Simulation: 100 generations
parameters.inputs
3 s
<Fis> vs distance to outlet
<Fis>,<Fit>,<Fst> vs generations
"Genetic distance" vs distance
FISH
- River network : RB=3, RL=2, n=3, scale=10
- Initial population : 5000 individuals
- Population: Brown trout, Salmo trutta fario
- Simulation: 100 generations
parameters.inputs
 31 s
 <Fis> vs distance to outlet
<Fis>,<Fit>,<Fst> vs generations
"Genetic distance" vs distance
FISH
- River network : RB=4, RL=2, n=3, scale=10
- Initial population : 5000 individuals
- Population: Brown trout, Salmo trutta fario
- Simulation: 100 generations
parameters.inputs
 42 s
 <Fis> vs distance to outlet
<Fis>,<Fit>,<Fst> vs generations
"Genetic distance" vs distance
MAMMAL
- River network : RB=3, RL=2, n=3, scale=10
- Initial population : 5 individuals
- Population: european mink, Mustela lutreola
- Simulation: 100 generations
parameters.inputs
 8 s
 <Fis> vs distance to outlet
<Fis>,<Fit>,<Fst> vs generations
"Genetic distance" vs distance
MAMMAL
- River network : RB=4, RL=2, n=3, scale=10
- Initial population : 5 individuals
- Population: european mink, Mustela lutreola
- Simulation: 100 generations
parameters.inputs
18 s
 <Fis> vs distance to outlet
<Fis>,<Fit>,<Fst> vs generations
"Genetic distance" vs distance

 
 

User manual

Download here.
 
 

Contacts

Audrey Chaput-Bardi
Cyril Fleurant
Christophe Lemaire
Jean Secondi