USER MANUAL - Computational Ecology Laboratory
Contents
1. ale Mdispmovepar 1 Used only for male dispersal Dispersal A function answer 5 negative ovement exponential y a 10 bx and is Parameter A the parameter a ale dispmovepar 1 Used only for male dispersal Dispersal B function answer 5 negative ovement exponential y a 10 bx and is Parameter B the parameter b Male Mdispmoveno max will A threshold option in cost Dispersal travel th ntir distance units for how far an ovement landscape as an individual male offspring can Threshold inverse square disperse You can specify max probability of to consider all individuals for cost distance on mating movement You can also the surface place an integer value in front of max to consider a percent cost distance movement for mating For example 10max would consider all mating individuals that are within 10 percent of the maximum cost distance on the surface You can also just specify a specific cost distance value Female Philopatry N philopatry If this is turned on Y Philopatry turned off then females will be able to remain in their birth locations If this is turned off N then males have the chance of coming into a female territory and taking it over When Y is specified then th generation will produce equal female and male occupied locations This option is only realistic and works with uniq2. WrightFisher The answer to hav very generation start with equal sex ratios CAREFUL this parameter is not realistic for non panmictic populations and should only be used to match Wright Fisher assumptions on equal sex ratios AtBirth This will ensure CDPOP user manual 20 that each litter is equal sex ratio but not necessary ensure that the generation will result in equal sex ratio after dispersal do to unbiased dispersal parameters and stochastic dispersal i e random offspring chosen for dispersal N This option is not used and offspring sex is assigned using the Female percent parameter above Birth Mortality birthmortper c 0 0 birth mortality Percent mortality in the offspring population Age Structure ortality agemortperc 1100 Age structured mortality here with non overlapping generations This is the mortality field for each age class Each age class must be separated with a If looking at non overlapping generations then specify numbers offspring mortality which is the same as the birth mortality in this situation and adult mortality where adult mortality is set to 100 If looking at overlapping generations with a supplied age distribution class e g age classes 0 1 2 3 4 then a mortality percentage must be entered for each age class e g 5
3. M opp requency nber of Infected The number o individuals that las a las a offspring disperse from their natal and males Ys E Ju Separated by ndividuals that are infected at each NFECT with given transmission probility of A used for CDEVOLVE f a used for CDEVOLV f B used for CDEVOLVE f b used for CDEVOLV E immigrate to the FemalesMeanMate The mean number of females matings in a generation MalesMeanMate The mean number of male matings in a generation FemalesSDMate The standard deviation in the number of variance in reproductive success matings in a generation F females MalesSDMate The standard deviation in the number of male variance in reproductive success matings in a generation OpenLocations The number of habitat available locations that CDPOP user manual 30 did not get filled during dispersal e CouldNotDispers The number of individuals that could not disperse do to high cost to moving to an open locations CDPOP user manual 31 5 General issues 5 1 How to obtain CDPOP The program is freeware and can be downloaded at http cel dbs umt edu software CDPOP with information for users including manual instructions FAQ publications ongoing research and developer involvement 5 2 Debugging and troubleshooting For
4. D Th GI DispDistl movement MateDistCD Th individuals travel to mate cost distance DispDistCD Th cost distance Euclidean distance Euclidean distance disperse from their natal location separated into female and male individual lation and Lly for lation and s travel to individual offspring distance if using a linear function as a probabil This can be converted back to cost as a probabil lity lity individual offspring disperse from their natal location separated into female and male movement distance if using a linear function MateDiststd The standard deviation Euclidean distance individuals travel to mate DispDiststd The standard deviation Euclidean distance individual offspring disperse female and males from their natal location This can be converted back to cost for both MateDiststd The standard deviation cost distance probability individual s travel to mate DispDiststd The standard deviation cost distance probability individual location for both females requency requency P1 This is th P2 This is th ql This is th requency q2 This is th SubpopMigration The nun ordered subpopulation number generation used for CD vertical infection
5. epd php for a free academic and educational usage in a single downloadable installer that has everything and then some Numpy Scipy Matplotlib and 70 modules for python We highly recommend this version for easy installs 2 2 Installation 2 2 1 Install Python NumPy and SciPy Make sure that Python and NumPy are installed and available to you You can test this by typing python at a command window If python is available you 11 get the python prompt gt gt gt If it is not a recognized command it means either that python is installed but is not in your command shell s paths or that python is not installed In the first case ask an administrator to add it to your command paths If your shell locates and loads python type import numpy Similarly type import scipy If python does not complain that there are no such modules all is well CDPOP user manual 8 The following instructions assume Python NumPy and SciPy are not yet available on your computer if they are skip to section 2 2 2 First run the Python executable installer you ve chosen either from www python org ActiveState or EPD accepting defaults for the installation directory On Windows this will typically place the executables and libraries in c Python2 7 bin and the site packages package tree for user installed Python modules in c Python2 7 1ib site packages If y
6. resolution of human evolutionary trees with polymorphic micorsatellites Nature 368 455 457 Cushman S A et al 2006 Gene Flow in Complex Landscapes Testing Multiple Hypotheses with Casual Modeling The American Naturalist 168 486 499 Cushman S A and Landguth E L 2010 Spurious correlations and inferences in landscape genetics Molecular Ecology 19 3592 3602 Holderegger R and Wagner H H 2006 A brief guide to Landscape Genetics Landscape Ecology 21 793 796 Landguth E L and Cushman S A 2010 CDPOP A spatially explicit cost distance population genetics program Molecular Ecology Resources 10 156 161 Landguth E L et al 2010a Quantifying the lag time to detect barriers in landscape genetics Molecular Ecology 19 4179 4191 Landguth E L et al 2010b Relationships between migration rates and landscape resistance assessed using individual based simulations Molecular Ecology Resources 10 854 862 Legendre P and Legendre L 1998 Numerical ecology 2 English ed Elsevier Amsterdam McRae B H and Beier P 2007 Circuit theory predicts gene flow in plant and animal populations Proceedings of the National Academy of Science USA 104 19885 19890 Nei M et al 1983 Accuracy of estimated phylogenetic trees from molecular data Journal of Molecular Evolution 19 153 170 Ray N 2005 PATHMATRIX a GIS tool to compute effective distanc
7. separating each generation to read in the next cost distance matrix by Future Mating CD Matrix Futuredispcm at 2 E EDcdmatrix16 DCdmatrixl6 A nxn future cost distance matrix for mating movement where n is the number of grid values or individuals on the landscape This is a comma delimited file format Note that this file must be the same size as the initial cost distance matrices used in the Simulations If you are using multiple surfaces i e you specified multiple cdclimate generations by separating them with then you will also need to specify multiple cost distance matrices here as well by using a Future Futuredispcd EDcdmatrix16 A nxn future cost distance CDPOP user manual 26 Dispersal CD Matrix mat EDcdmatrixl6 matrix for dispersal movement where n is the number of grid values or individuals on the landscape This is a comma delimited file format Note that this file must be the same size as the initial cost distance matrices used in the Simulations If you are using multiple surfaces i e you specified multiple cdclimate generations by separating them with then you will also need to specify multiple cost distance matrices here as well by using a ED Matrix Edmatans N a Eu
8. settles location specific fitness values fo pre processing step The program wi same as CDPOP with an additional s dispersal process E Simulating dynamic landscapes The on the connectivity of populations among scientists and land managers and spatially explicit predictions g n Mendelian reproduction In irs are one male to possibly many ating process occurs when all can have a number of offspring on a uniform probability pecified mean or a constant k allele mutation rate chosen by ffspring s genotype and the sex the mother s x y location function and the sex of the and death define whether or not or immigrants versions of CDPOP modeled three flow genetic drift and mutation nt genotypes have an equal g on their alleles to f future tion was not operating CDPOP V 1 0 ously to the adaptive or fitness ight 1932 This new functionality tic analyses to explicitly flow and selection in complex The user specifies fitness e of a single diallelic locus that ree relative fitness surfaces must s AA Aa and aa from the two emented through dif F de erential as a function of the relative f where the dispersing Lac CDPOP v1 2 reads and extracts genotype and r each n x y individual in the 11 continue all other processes the tep implement selection during the potential i
9. that the offspring has That value then becomes the individual offspring mortality percentage An offspring becomes more or less fit relative to the other offspring at that generation as a function of its genotype and where it occurs on a surface These are ASCII formatted files with 6 lines of header information and values that are space delimited and represent percent mortality 0 100 See fitvals txt for an example fitness surface and its format Fitness Fitness AA Fitvals0 When CDEVOLVE Answer is 1 then AA this is the offspring viability pa selection surface for AA offspring has AA then this THe e mortality fitness surface is used Fitness Fitness Aa Fitvals0 When CDEVOLVE Answer is 1 then Aa this is the offspring viability E selection surface for Aa offspring has Aa then this tT E mortality fitness surface is used Fitness Fitness aa Fitvals0 When CDEVOLVE Answer is 1 then aa this is the offspring viability fe selection surface for aa offspring has aa then this mortality fitness surface is used TE The following are the fitness surfaces for when CDEVOLVE Answer is 2 This corresponds to 2 loci that are under selection The x y location of an offspring is matched up with the closest x y fitness value for the offspring s corresponding fitness surface defi
10. the default The cdpop format lists the genotypes with values f for each allel gt either 0 1 or 2 Output will be labeled grid generation csv F general for a general genotype output The general format will ollow Locusla Locuslb Locus2a Locus2b LocusNa LocusNb Output will be labeled generalgrid generation csv E genalex for the program GENALEX Output will be labeled genalexgrid generation csv structure for the program structure or related programs Output will be labeled structuregrid generation stru genepop for the program GENEPOP or related programs Output will be labeled genepopgrid generation txt 4 Folders will be created in your project directory label Output Files unique time stamp dos convention e g Monte Carlo runs will be uniquely labeled mcrun0 batch run will be uniquely labeled batchr folder you wil individual s genotype spatial locations status generation The following is a summary of In addition an output csv is automatical Monte Carlo run These are population bas THe you specified a grid format op e CDPOP user manual 28 ed with a 1332964297batchrunOmcruno un0 batchrunl mcrunl and each In each ll see grid generation csv files that list each F age sex
11. to the left Ve forwardbackward This is a step wise mutation in which an allele can mutate forward or backwards only i e to the left or right with equal probability forwardAbackwardBrandomN This is a special case for the 2 loci selection model The first locus under selection can only go forward A gt a and the second locus under selection can only go backward b gt B The rest of the neutral loci are random mutations Loci loci 110 The number of loci microsatellites Initialize Genotypes intgenesans random The choice for how to initialize the genotype for CDPOP user manual 22 each n x y individuals If random is entered then the genotypes get a random assignment and the population is at a maximum genetic diversity If file is entered then the genetics get drawn from the allele frequency distribution ile specified in next column allefregfilename PUE known is entered then th genotypes are directly read from a given known file This file is very similar to the initial xyfilename and example xyED16 known csv is supplied with test data files Allele Frequency File Allefreqfile name N allele frequency file not used The allele frequency distribution for each
12. user manual pdf this file gt CDPOP history txt Notes on history and version changes 3 data Example input files gt Edcdmatrix16 csv example Euclidean distance cost distance matrix file gt xyED16 csv example n x y file for individuals gt xyED16 NAs csv example n x y file with sparse individuals gt xyED16 known csv xample n x y file with known genetic data gt agedistribution csv example age distribution file CDPOP user manual 9 allelefrequency csv example allele frequency distribution LLE fitvals0 txt and fitvals50 txt example fitness landscapes ves for natural selection inputvariablesl6pnts csv run parameters corresponding to the example files 2 3 CDPOP user manual 10 Example run 2 3 1 Command line run The example run is E or 16 points representing individuals with a cost distance matrix calculated with Euclidean distance To run the following example Tee F ollow these steps Double check that the 3 directories provided in the archive are in the same directory The included file inputvariblesl pnts csv in the data directory specifies the parameters that can be changed and used in a sample CDPOP run Open inputvariblesl pnts csv in your editor of choice A spreadsheet program like Microsoft Excel allows for easy editing of the tabular values
13. 0 20 20 30 100 Note that if you do not put 100 for the final age class then the population will be allowed to go beyond your last age class specified All age classes beyond specified in file will get lumped into last mortality category Subpopul latio n Mortal ity Subpopmortpe re v0 0 0 0 This parameter is the percent mortality for a dispersing offspring into another subpopulation Each subpopulation gets separated by a and has it s own percent mortality For example if an offspring was born in the CDPOP user manual 21 second subpopulation and disperses to the fourth subpopulation it has a 60 chance of surviving there If it stays in its own subpopulation then it would have no differential mortality consequences This parameter is very similar to cdevolve offspring viability but applied to subpopulations instead of individuals spatial locations Mutation Rate muterate 0 0005 The mutation rate Mutation Model mutationtype random the KAM model The type of mutation model random This is the kth allele mutation model Ve forward This is a step wise mutation in which an allele can mutate forwards only i e to the right backward This is a step wise mutation in which an allele can mutate backwards only i e
14. ABb Fitness AABb Fitvals0 When CD this is selection sur offspri the offspri EVOLVE Answer is 2 then ng viability P ace for AABb f ng has AABb F mortality fitness s used then this urface is Fitness AaBb Fitness AaBb Fitvals0 When CD this is selection sur offspri EVOLVE the offspri Answer is 2 then ng viability F ace for AaBb He ng has AaBb Pai mortality fitness s used then this urface is Fitness aaBb Fitness_aaBb Fitvals0 When CD this is selection sur offspri the offspri EVOLVE Answer is 2 then ng viability E ace for aaBb f ng has aaBb pl mortality fitness s used then this urface is Fitness AAbb Fitness AAbb Fitvals0 When CD this is EVOLVE the offspri Answer is 2 then ng viability selection surface for AAbb If offspri ng has AAbb E mortality fitness s used then this urface is Fitness Aabb Fitness_Aabb Fitvals0 When CD this is the offspri EVOLVE Answer is 2 then ng viability F selection surface for Aabb If offspri ng has Aabb E mortality fitness s used then this urface is Fitness aabb Fitness aabb Fi
15. CDPOP user manual 1 CDPOP USER MANUAL 2012 Version 1 2 08 Last Updated 2012 05 20 Authors E L Landguth B K Hand J M Glassy S A Cushman M Jacobi T J Julian Allen Warren 1 University of Montana Division of Biological Sciences Missoula MT 59812 USA 2 Lupine Logic Inc Missoula MT 59802 USA 3 U S Forest Service Rocky Mountain Research Station 2500 S Pine Knoll Dr Flagstaff AZ 86001 USA CDPOP user manual 2 Table of Contents I gt AMECLOQUELLON hi A AA A Sie Ee eed eae and ered ers 3 Ll Changes trom CDPOP VO spa eae a tees Ses a eae eee 3 t27 What Can CDLOP do iets fe sie ie Ge OOS Seat ew RN ee ee wd Se ee ree 3 123 How does lt CDPOP WORK siii iS dares awe wea owes 4 2 Getting Started its is el el GS wld nd WAS Ae A pea A 7 2rd DEPENGENCLES 4 ernes Sik eee ORAM E E E EEE ee es 2 1 1 Baseline requirements 2 122 EYyeaon On nOn wIindows platos snes dias a Ye 2 4 4253 PYTHON On WINAOWS se sats li le 2 lt 4 0btainang NumPy and SEL Y 2 dees ewis ees Les 262 Tnstaltat ron o cht ot ie tae te tie bt ante JA hy Mi ah Ne ce MG ita 2 2 1 Installing Python NumPy and SciPy Zi Za 2 Instaldang CDPOP 6 a Ver 2 abe 2 22 23 DESCETPELON OF EDROD TIMES a Berens cee Obey E 23 HB KAMP LS PUM gis ie ences a A es SS eth a eb ta 2003 Command LENS CUA ii a Snes 3 D
16. PI E Gate A ah Gree od oA a See eS 12 Ar QUE DU ete a a a E AS ERRADA DA T a ana eee a alese d 28 5 General ISSUES eit tase is oa eee ha ad ada 31 Diy A HOW OOD Pain CDROPD sw tai ae tee Mush ine bea eee inmate doe Daya Diz Debugging and EroubleshootinBui idaho 55 3 How to cite CDPOP 5 4 DaAsClLaARMe Pie 6 wht teow wos Sate Sh aoe lS SSS DEAS wi we etre E Referentes eere Seabee A Al Ai 32 7 Acknowledgements CDPOP user manual 3 1 Introduction The goal of this user manual is to explain the technical aspects of the current release of the CDPOP program CDPOP v1 2 is a major extension of the CDPOP program Landguth and Cushman 2010 CDPOP is an individual based program that simulates the influences of landscape structure on emergence of spatial patterns in population genetic data as functions of individual based movement breeding and dispersal 1 1 Changes from CDPOP v0 7 There are major innovations in v1 2 which were not included in the previously published v0 70 We list below the new functionalities of CDPOP v1 2 e Natural selection is implemented through differential offspring viability as functions of fitness landscapes e Gene flow and natural selection can now be simulated in dynamic landscapes e A graphical user interface provides a user friendly platform that enables users to explore analyze and model the effects of life history and differential models of
17. TRIX CIRCUITSCAPE UNICOR COSTDISTANCE etc Dispersal CD dispcdmat EDcdmatrixl6 A nxn cost distance matrix Matrix an example for dispersal movement where n Filename Euclidean distance is the number of grid values matrix used for or individuals on the the dispersal landscape This is a comma movement delimited file S th xample EDcdmatrixl6 csv for formatting this file Also note that this file can be calculated from any program you choose e g PATHMATRIX CIRCUITSCAPE UNICOR COSTDISTANCE etc It can be the same file as the Mate CD Matrix Monte Carlo mcruns 5 5 replicate The repeated number of Replicates runs denoted in Simulations to be conducted for folders labeled the Monte Carlo method i e with mcrun0 the number of replicates for 1 mcrunl batch of parameters Generation T Looptime YLO o DO Simulation run time generation ime generations or year File output indexed from 0 looptime 1 For CDPOP user manual 14 example grid0 csv gridl csv grid2 csv grid3 csv and grid4 csv would be outputted for a looptime of 5 Generations nthfile choi Sequence will The choice of a specified of Saved ce save a specified simulation run time Genotypes sequence of generation to write to file Choice generations and to calculate genetic distance matrices List read nthfile list values below Sequence read nthfile seq value b
18. There will be 3 lines of information in inputvariables csv a header line and 2 lines of information corresponding to 2 separate CDPOP runs batch process See the user manual pdf that contains a breakdown for each column header and the parameters that can be changed The Input listed is for the first row in the file Make sure you save inputvariables in the same format a comma delimited file Select Yes or OK for any Excel questions about saving in this format P Start the program For example if you use python from the command line then open a terminal window and change your shell directory to the CDPOP src home directory i e gt cd C homedirectorylocation src Run the program There are a number of ways to run this program If you are using a command shell you can run the program by typing python CDPOP py C homedirectorylocation data inputvariableslopts csv output_test Note that there are 5 arguments here that must be included with spaces in between gt python starts python for example from the command line Note that other python environments may have different calls here In PyLab the IDE distributed with EPD the call is run gt CDPOP py runs CDPOP program gt C homedirectorylocation data is the directory location of lt lt the input test files You can point this directory to other project files for example We suggest not h
19. al A function answer 5 negative ovement exponential y a 10 bx and is Parameter A the parameter a Female Fdispmovepar 1 Used only for female dispersal Dispersal B function answer 5 negative ovement exponential y a 10 bx and is Parameter B the parameter b Female Fdispmovethr 5max will A threshold option in cost Dispersal esh move 75 of the distance units for how far an ovement maximum cost individual female offspring can Threshold distance on that disperse You can specify max surface to consider all individuals for mating movement You can also place an integer value in front of max to consider a percent cost distance movement for mating For example 10max would consider all mating individuals that are within 10 percent of the maximum cost distance on the surface You can also just specify a specific cost distance value Male Mdispmoveno 2 inverse This is the function answer for Dispersal square probability movement for male dispersal Movement function for male Choice cost distance 1 Linear 1 1 Threshold dispersal Cost Distance 2 Inverse Square 1 Cost Distance 2 1 CDPOP user manual 18 3 Nearest Neighbor Use threshold to specify Moore neighborhood 4 Random Mixing 5 Negative Exponential parA 10 parB Cost Distance
20. and infection tion then you will additionally see the format followed by grid generation csv F ly created for each batch and ed metrics calculated at each each calculation Year This is the generation time or year if using overlapping generations Population The total population in each generation If you specify subpopulations in the first the subpopulation sizes in order E column of your xy csv file then this field will be separated by The first value is the total population size and each additional value corresponds to ToTFemales The total number of females in each generation ToTMales The total number of males in each generation BreedFemales The total number of breeding age females in each generation BreedMales The total number of breeding age females in each generation Fema BreedEvents This is the number of breeding events for females This number times the birth rate will give the total number of births in that generation it to the next generation E Births The number of of Migrants The number of dispersers in each generation that make DisperseDeaths The number of dispersers in each generation that do not make it to the next generation spring born at that generation AdultDeaths The number of deaths of the adult population not the offspring born that year generation The se
21. answer 5 negative exponential y a 10 bx and is ating ovement Parameter B matemoveparB 0 01 the parameter a This is only used for mating movement answer 5 negative exponential y a 10 bx and is the parameter b ating ovement Threshold matemovethre sh A threshold option in cost distance units for how far an individual can search for a mate equivalent to the cost distance kernel You can specify max to consider all individuals for mating movement You can also place an integer value in front of max to consider a percent cost distance movement for mating For example 10max would consider all mating individuals that are within 10 percent of the maximum cost distance on the surface You can also just specify a specific cost distance value Female Replacement Freplace N females mate without replacement F If you want females to mate with replacement then specify NYT TE E females to mate then you want without replacement specify N Male Replacement Mreplace Y males mate with replacement If you want males to mate with then specify Y replacement If you want males to mate without replacement then specify N Multiple Paternity Answer multiple pat ernity y females can have a litter from mul
22. aving any spaces in your directory names gt inputvariablesl6pts csv is the parameter file gt output_test is the name of the directory that will be created with CDPOP output in the directory specified by the third argument above m CDPOP user manual 11 Check for successful model run completion The program will provide step by step output in the Shell window Once completed a simulation time will be printed out and folders batchrunOmcrunO batchrunOmcrunl batchrunOmcrun2 batchrunOmcrun3 batchrunOmcrun4 and batchrunlmcrun0O will be created in your CDPOP home directory to store output from the separate batch and or Monte Carlo runs These folders are located in the data folder specified in 6 The output folder will have a unique date time stamp preceding in case you want to run multiple CDPOP runs in this same directory The program will also provide a log file with program steps in your CDPOP home directory If parameters are such that population becomes extinct before specified generation time then program will end 3 Input The foll CDPOP are for CDPOP user manual 12 owing are the general input parameters and files used in See examples provided for F the GUI each fil le or parameter and the first row ormatting Th in the inputvariables csv describing file headers listed T
23. clidean distance matrix is not calculated This is the Euclidean distance matrix answer This matrix will be called EDmatrix csv and outputted in the respectiv batch and Monte Carlo folder If Y then the Euclidean distance matrix used for the n original grid locations is calculated If N is entered then this matrix is not created GD Matrix Gendmatans Dps proportion of shared alleles genetic distance matrix is produced for every grid csv file This is the genetic distance matrix answer The genetic distance matrix used for the n original grid locations for specified generation time of the simulation run is calculated E nter braycurtis for the Bray Curtis distance measure Dps for the proportion of shared alleles or Da for Nei s genetic distance F Specify N here if you do not want to calculate these matrices Note that Nei s genetic CDPOP user manual 27 distance takes the longest to calculate and may decrease your total CDPOP simulation time All are slightly expensive in computation time and we are working on speeding these up Genotype Output Format gridformat genepop This is the genotype output format option The format for the genotype output is E specified by entering cdpop This format is
24. complex landscapes on the genetic structure of populations e Demographics allows for fluctuating population sizes FE e Sex specific dispersal F e Changes of some internal software components have allowed an overall speed increase and to improve program stability e Additional movement function option negative exponential movement e Inclusion of both a mating landscape and a dispersal landscape e Vertical transmission of an infection with giving probability e mtDNA option e Output genotype option in multiple formats genalex structure genepop general e Generation can be specified to change landscape surface for both mating and dispersal e g climate change scenarios e Subpopulation differential mortality An individual can be given an increased mortality if it disperses to another subpopulation e Philopatry option added e Multiple paternity option added i e females can have litters from multiple males but birth rate is same for each female e Additional mutation models added step wise forward and backward e Age structured mortality added for overlapping generations 1 2 What can CDPOP do CDPOP s realistic representation of the spatial environment and popu lation genetic processes provide a powerful framework to investigate CDPOP user manual 4 the impact of ecological factors on the genetic structure of popula tions Th
25. ected according to the movement function choice and without regard to any gender or mating type It is important to note that this asexual reproduction is functionally the sexual paring of hermaphroditic individuals thus it is asexual in the sense that there are no distinct sexes but is sexual in the sense that individuals mate with other individuals exchanging genes in Mendelian reproduction Reproduction Age reproage 0 individuals start reproducing at age 0 The age at which individuals can start to reproduce Use with overlapping generations i e oldmortperc not set to 100 If an age distribution file is not used then age gets initialized with 0 at the CDPOP user manual 17 starting generation Female Fdispmoveno 12 inverse This is the function answer for Dispersal square movement for female dispersal Movement probability probability Choice function for female cost 1 Linear 1 1 Threshold distance movement Cost Distance 2 Inverse Square 1 Cost Distance 2 1 3 Nearest Neighbor Use threshold to specify Moore neighborhood 4 Random Mixing 5 Negative Exponential parA 10 parB Cost Distance Female Fdispmovepar 1 Used only for female dispersal Dispers
26. elow List of nthfile list 0 3 4 the list The specified simulation run Generations of generations to time to write to file and to Saved save and here calculate genetic distance cdpop would save grid0 csv grid3 csv grid4 csv and matrices These values are used if nthfile choice List These values must be separated with a vertical bar Sequence of nthfile seq 1 will save The specified simulation run Generations the sequence of time to write to file and to Saved genotypes for calculate genetic distance grid0 csv matrices This value is used if gridl csv nthfile choice Sequence grid4 csv This is the by value in the sequence For example 1 would create values starting at 0 ending at looptime 1 by every 1 generation A value of 5 would output every 5 generation etc Mating matemoveno 1 linear Movement function answer for Movement probability mating probabililty Function function of mating cost distance 1 Linear 1 1 Threshold Cost Distance 2 Inverse Square 1 Cost Distance 2 1 3 Nearest Neighbor Use threshold to specify Moore neighborhood 4 Random Mixing 5 Negative Exponential parA 10 parB Cost Distance CDPOP user manual 15 ating ovement Parameter A matemoveparA 0 0005 pal This is only used for mating movement
27. es among samples Molecular Ecology Notes 5 177 180 1 Storfer A et al 2010 Landscape genetics where are we now Molecular Ecology 19 3496 3514 Wright S 1932 The roles of mutation inbreeding crossbreeding and selection in evolution Proceedings XI International Congress of Genetics 1 356 366 7 Acknowledgements This research was supported in part by funds provided by the Rocky Mountain Research Station Forest Service U S Department of Agriculture and by the National Science Foundation grant DGE 0504628
28. he example provided is for the first line in the inputvariables csv file File Header File Header Example Description in GUI in csv XY Filename xyfilename xyED16 example The n x y grid location supplied for 16 values This is a comma individuals delimited file with 5 column headings xyED16 NAs csv gives you an Subpopulation a unique example of how to identifier for each individual specify your corresponding to a unique initial starting subpopulation This is an population size optional tracker for This example is 16 individuals that may be located possible habitat in designated subpopulations locations that can If individuals are just be filled with continuously distributed then only 11 fill with arbitrary value like initialized with 1 If these field is used individuals You then subpopulations must be in must specify NA sequence e g 1 2 3 values in the ID and sex columns XCOORD x coordinate location but still provide the xy locations YCOORD y coordinate location that are OPEN and YCOORD a Supopulation identifier field ID a string label identifier and sex an initial sex assignment use 0 1 or F M See xyED16 csv for an example xyfilename The column order is necessary and header file included See below for s
29. help with installation problems please check first for postings at our web site Otherwise please report problems including any bugs to me at erin landguth mso umt edu 5 3 How to cite CDPOP This program was developed by Erin Landguth with help from Brian Hand Joe Glassy Sam Cushman and Tyler Julian GUI development was done by Mike Jacobi The reference to cite is as follows substituting the version number Landguth EL Cushman SA Jacobi M 2010 CDPOP A spatially explicit cost distance population genetics program Molecular Ecology Resources 10 156 161 Version X Landguth EL Cushman SA Johnson NJ 2011 Simulating natural selection in landscape genetics Molecular Ecology Resources doi 10 1111 3 1755 0998 2011 03075 x 5 4 Disclaimer The software is in the public domain and the recipient may not assert any proprietary rights thereto nor represent it to anyone as other than a University of Montana produced program version 1 x CDPOP is provided as is without warranty of any kind including but not limited to the implied warranties of merchantability and fitness for a particular purpose The user assumes all responsibility for the accuracy and suitability of this program for a specific application In no event will the authors or the University be liable for any damages including lost profits lost savings or other incidental or consequen
30. is approach has already advanced knowledge of the patterns of genetic variation in spatially explicit contexts Landguth et al 2010a Landguth et al 2010b Cushman and Landguth 2010 Landguth and Balkenhol 2012 Example simulations have included e Quantifying the time to detect barriers e Correlating migration rates and landscape resistance barriers e Testing for the effects of population sample size and number of markers e Assessing relative influence of adaptive versus neutral markers in detection of population genetic differentiation 1 3 How does CDPOP work CDPOP v1 2 models genetic exchange for a given resistance surface and n x y located individuals as functions of individual based movement through mating and dispersal vital dynamics and mutation A user must specify the input parameters through a graphical user interface or input script file As the model simulates stochastic processes most applications will quantify mean and variability of genetic structure across many runs Thus a Monte Carlo option is provided for the user to choose the number of runs to simulate given a single set of input parameters In addition a user may also frequently wish to launch several runs with different parameter values Simultaneously i e sensitivity analysis This functionality is provided through batch capability The simulation program assumes consta
31. locus used to initialize the model s n individual s genotype If you want to use a frequency distribution file you must set Initialize Genes Answer to equal file and then enter in the filename in this field See allelefrequency csv example ile for formatting this file It is basically a column of allele frequencies and make sure the length of the column equals your starting loci starting alleles Alleles alleles VA The number of starting alleles per locus mtDNA Mtdna EN If Y then last locus becomes mtDNA and every offspring inherits this locus from its mother only If N then regular Mendal inheritance occurs for this last locus CDEVOLVE Answer cdevolveans This is the answer for how many loci are under selection Use N to turn off CDEVOLVE Use 1 for natural selection CDPOP user manual 23 with 1 locus Use 2 for selection with 2 loci Alleles must be 2 if Y is entered The following are the fitness surfaces for when CDEVOLVE Answer is 1 This corresponds to 1 locus that is under selection The x y location of an offspring is matched up with the closest x y fitness value for the offspring s corresponding fitness surface defined by the genotype
32. mpacts of climate change have become an area of concern Current needs include quantitative of current and potential future patterns of fragmentation under a r Opdam amp Wascher 2003 To address to input a new landscape surface at new cost distance matrices for both ange of climate change scenarios this need CDPOP v1 2 allows users a given generation time through mating and dispersal CDPOP user manual 6 Simulating demographics CDPOP v1 2 now allows for fluctuating population sizes controlled by birth rate death rate carrying capacity and ability of offspring to disperse given a specified resistance surface The program is written in Python 2 7 and provided with installation instructions for most platforms along with sample input files CDPOP v1 2 is built on a driver module plug in docking architecture that allows for ease of future modular development CDPOP v1 2 has been debugged as carefully as possible by testing all combinations of simulation options Information for users including user manual FAQ publications ongoing research developer involvement and downloads can be found at http cel dbs umt edu software CDPOP CDPOP user manual 7 2 Getting started 2 1 Dependencies 2 1 2 Baseline Requirements CDPOP requires the Python2 7 x interpreter NumPy package and SciPy package Remember that Python modules usually require particular Python interpreter
33. ned by the genotype that the offspring has That value then becomes the individual offspring mortality percentage An offspring becomes more or less fit relative to the other offspring at that generation as a function of its genotype and where it occurs on a surface These are ASCII formatted files with 6 lines of header information and values that are space delimited and represent percent mortality 0 100 See fitvals txt for an example fitness surface and its format Fitness Fitness AABB Fitvals0 When CDEVOLVE Answer is 2 then AABB this is the offspring viability selection surface for AABB f offspring has AABB then this F mortality fitness surface is used Fitness Fitness AaBB Fitvals0 When CDEVOLVE Answer is 2 then AaBB this is the offspring viability CDPOP user manual 24 selection offspri surface for AaBB f E mortality fitness s used ng has AaBB then this urface is Fitness aaBB Fitness_aaBB Fitvals0 When CD this is selection sur offspri EVOLVE the offspri Answer is 2 then ng viability F ace for aaBB f ng has aaBB F mortality fitness s used then this urface is Fitness A
34. nt population density over time Individuals are assumed to occupy a fixed grid on the landscape that is user defined by the n x y located individuals The genotype of each locus for each individual can be initialized by randomly choosing from a file containing allel frequencies for each locus or by reading in a file containing the initial multi locus genotypes of all the individuals The initial age structure of the population is specified by an input file specifying initial age frequency The sex of each initial individual is randomly assigned P There are five movement functions that define how individuals choose a mate and disperse on the landscape as a function of cost distance linear inverse square negative exponential nearest neighbor and random mixing With the nearest neighbor movement function an individual moves to the available grid location nearest its initial location Random mixing moves an individual to a grid location that is randomly chosen from the n grids in the population In linear inverse square and negative exponential movement functions individuals move a distance from their initial location based on a draw from a probability distribution inversely proportional to a linear inverse square or negative exponential function The user specifies the maximum dispersal distance in cost units an individual can travel on the landscape The probability is one at no dis
35. ou are installing it on a network on which you do not have administrative privileges you may need to ask a system administrator to install python and the NumPy and SciPy packages in their default locations Next install NumPy and SciPy using the supplied executable superpack installer or visiting http www scipy org Download This will install NumPy and SciPy in your Python site packages directory Note that if you choose EPD that you do not need to additionally install NumPy or SciPy 2 2 2 Install CDPOP Next install the CDPOP software itself by unpacking the zip archive supplied Navigate to the directory on your PC where you wish to install CDPOP and unpack the supplied zip archive file using a free archive tool like 7Zip 7z exe Pkunzip Unzip or an equivalent Seven Zip 7Z exe is highly recommended since it can handle all common formats on Windows MAC OS X and Linux On Windows it is best to setup a project specific modeling subdirectory to perform your Simulations outside of any folder that has spaces in its name like My Documents At this point you should be able to execute th supplied test inputs 2 2 3 Description of CDPOP files 3 directories will be installed in your directory Here is a description of each 1 src CDPOP source code 2 doc gt README txt a quick how to run CDPOP instructions gt CDPOP
36. parates th age class deaths for overlapping generations for th ordered ag class in that generation Note that some years may not have all age classes and not reported Alleles This is the total number of uniqu generation This value can be calcul alleles at each lated automatically for specified subpopulations if different subpopulations were F designated in the initial xy ilename If there are first value corresponds to the total and subpopulation values follow after then the alleles in the population He This is th xpected heterozygosity value at each generation This value can be calcul lated automatically for specified subpopulations if different subpopulations were designated in the initial xyfilename If there are Sle then the CDPOP user manual 29 first value corresponds to the total He in the popul subpopulation values follow after Ho This is th observed heterozygosity value at each generation This value can be calculated automatical specified subpopulations if different subpopulations were then the designated in the initial xyfilename If there are first value corresponds to the total Ho in the popul ow after subpopulation values foll Mutations The total number of mutations at each generation MateDistED Th mate
37. pecifying constant versus non constant population sizes Age agedistribut N no age The distribution that is used Distribution ion distribution file to initialize each individuals Filename is used and age is age initialized at 0 for all N no file is not used and indidivuals age initialized with 0 CDPOP user manual 13 Filename If a filename is entered then read in the file for example agedistribution would be entered for the example provided See the agedistribution csv for formatting this file and note it must be comma delimited The file includes two columns the Age class start initial age of offspring as 0 and Number in age class which is the number of individuals in each class This number must sum to the carrying capacity of your study or the total number of xy location specified above Mate CD matecdmat EDcdmatrix16 A nxn cost distance matrix Matrix an example for mating movement where n is Filename Euclidean distance the number of grid values or matrix used for individuals on the landscape the mating This is a comma delimited file movement S th xampl EDcdmatrixl6 csv for formatting this file Also note that this file can be calculated from any program you choose e g PATHMA
38. s so be sure the version ID for any external Python module or package e g NumPy or others matches the version of your Python interpreter normally v2 7 x 2 1 2 Python on Non Windows Platforms Some common computer platforms come with Python installed These include MAC OS X and most Linux distributions To determine which Python a MAC or Linux workstation has installed start a terminal console and enter python You 1ll see the version number on the top line enter Control D to exit Replacing an older Python interpreter pre v2 4 with a newer one v 2 7 x on a Linux or MAC OS X machine can be tricky so ask a System Administrator for help if you re not sure which packages depend on the current Python installed 2 1 3 Python on Windows Windows 7 XP 2000 Server does not come with Python installed so follow the instructions below to obtain and install Python ona computer running the Windows operating system Get a windows installation of the base Python installation current v 2 7 x at http www python org download releases 2 1 4 Obtaining NumPy and SciPy A good suggestions for all in one installations is through the superpack Windows installer available from the SourceForge websit http sourceforge net project Note that more complete information for NumPy is available at www scipy org where the SciPy module is also presented Another source is http www enthought com products
39. tance from the original location and goes distance Qo Reproduction is defined by the user heterosexual With hermaphroditic m but individuals mate with other ind CDPOP user manual 5 to zero at the maximum dispersal as either hermaphroditic or ating there are no distinct sexes ividuals according to the movement F unction choice exchanging genes i heterosexual reproduction mated pa or no females and the end of the m females have mated Each mated pair that is a bounded random draw based distribution a Poisson draw with s number Mendelian inheritance with the user is used to generate the o assignment is random Dispersal of offspring occurs from according to the selected movement individual The vital rates birth the population will have emigrants Simulating natural selection Past sources of genetic variation gene These versions assumed that differe probability of surviving and passin generations and thus natural selec implements natural selection analog landscape of allele frequencies Wr nables extension of landscape gene investigate the links between gene landscapes at an individual s level landscape surfaces for each genotyp is under selection For example th be specified for the three genotype alleles A and a Selection is impl survival of dispersing individuals fitness at the location on that sur individual
40. tial damages arising from the use of or the inability to use this program We strongly urge you to read the entire documentation before ever run ning CDPOP We wish to remind users that we are not in the commercial software marketing business We are scientists who recognized the need for a tool like CDPOP to assist us in our research on landscape ecolo gy issues Therefore we do not wish to spend a great deal of time consulting on trivial matters concerning the use of CDPOP However we do recognize an obligation to provide some level of information sup port Of course we welcome and encourage your criticisms and sugges tions about the program at all times We will welcome questions about how to run CDPOP or interpret the output only after you have read the entire documentation This is only fair and will eliminate many trivi al questions Finally we are always interested in learning about how F CDPOP user manual 32 others have applied CDPOP in ecological investigation and management application Therefore w ncourage you to contact us and describe your application after using CDPOP We hope that CDPOP is of great assistance in your work and we look forward to hearing about your applications Shiny 6 References Allendorf F W and Luikart G 2007 Conservation and the genetics of populations Blackwell Malden MA Bowcock A M et al 1994 High
41. tiple males If you want multiple paternity with birth rates applied the same for each female with litter then specify Y If you want multiple paternity with birth rates applied unequally for each female with litter then specify N This answer is only functionly when Freplace Y and CDPOP user manual 16 Mreplace Y or multiple paternity option The difference is subtle For example if the birth rate is set at a constant litter size offno 3 of 2 lambda 2 and multiple paternity Y then if a female mates with 2 males the female would have 2 offspring 1 from each male However if multiple paternity N then if a female mates with 2 males then that female would have 4 offspring Selfing Answer Selfans N self ing is turned off If you want to allowing selfing i e individuals mate with themselves then specify Y If you do not want to allow for selfing then specify N Reproduction Answer sexans y sexu reproducti choosen al on is Y for sexual reproduction In sexual reproduction mated pairs consider male and females with or without replacement N for asexual reproduction With asexual reproduction all n individuals mate and bear offspring with mates sel
42. tvals0 When CD this is EVOLVE the offspri Answer is 2 then ng viability F selection surface for aabb If offspri ng has aabb mortality fitness s used then this urface is CDINFECT cdinfect SN x turned off This is the infection parameter answer This tracks vertical transmission in the population A column in grid csv denotes the infection status at each CDPOP user manual 25 generation for every individual If Y then a random status infection 0 or 1 is created and initialized for each individual If N then the status 0 is created for all individuals and initialized Transmission Probability Transmission prob Ogay This is the transmission probability for if a parent has the infection the chance that the infection will be passed along to the offspring CDCLIMATE cdclimate N changing cost distance matrices turned off This is the dynamic landscape answer If Y then a new cost distance matrix will be read in at a specified generation time in the next column TE ON used then this option is not CDCLIMATE Generation Cdclimgentim e valo The generation time that the next cost distance matrix will be read in at You can specify multiple generations by
43. ue combinations of parameters For example female dispersal threshold is 0 no movement You must start off with CDPOP user manual 19 initialized equal sexes You must specify equal sex ratio AtBirth You must specify female mating without replacement You must use male movement 1 2 or 5 Offspring Choice Offno 12 Poisson distribution This is the number of offspring each mate pair can have 1 for a random draw 2 for Poisson draw and 3 for a constant number of offspring for each mother or mate pair Lambda Lmbda 5 Poisson distribution with mean lambda 5 The parameter value used with Offspring Number If offno 1 then lmbda is the max range value between 0 imbda to draw randomly from If offno 2 then lmbda is the Poisson mean for the litter size If offno 3 then lmbda is the constant litter size value Female offspring Femalepercen E 150 50 random female assignment Percent number of female born in each litter This is a random assignment from given percentage i e even if you set this to 50 some generations could have 499 female births and 501 male births for example Equal Sex Ratio for Offspring Equalsexrati oBirth N This ensures an exact equal sex ratio for the following options
Download Pdf Manuals
Related Search