GENETIC ANALYSIS OF COMPLEX DEMOGRAPHIC SCENARIOS: SPATIALLY EXPANDING POPULATIONS OF THE CANE TOAD, BUFO MARINUS

Abstract Inferring the spatial expansion dynamics of invading species from molecular data is notoriously difficult due to the complexity of the processes involved. For these demographic scenarios, genetic data obtained from highly variable markers may be profitably combined with specific sampling schemes and information from other sources using a Bayesian approach. The geographic range of the introduced toad Bufo marinus is still expanding in eastern and northern Australia, in each case from isolates established around 1960. A large amount of demographic and historical information is available on both expansion areas. In each area, samples were collected along a transect representing populations of different ages and genotyped at 10 microsatellite loci. Five demographic models of expansion, differing in the dispersal pattern for migrants and founders and in the number of founders, were considered. Because the demographic history is complex, we used an approximate Bayesian method, based on a rejection‐regression algorithm, to formally test the relative likelihoods of the five models of expansion and to infer demographic parameters. A stepwise migration‐foundation model with founder events was statistically better supported than other four models in both expansion areas. Posterior distributions supported different dynamics of expansion in the studied areas. Populations in the eastern expansion area have a lower stable effective population size and have been founded by a smaller number of individuals than those in the northern expansion area. Once demographically stabilized, populations exchange a substantial number of effective migrants per generation in both expansion areas, and such exchanges are larger in northern than in eastern Australia. The effective number of migrants appears to be considerably lower than that of founders in both expansion areas. We found our inferences to be relatively robust to various assumptions on marker, demographic, and historical features. The method presented here is the only robust, model‐based method available so far, which allows inferring complex population dynamics over a short time scale. It also provides the basis for investigating the interplay between population dynamics, drift, and selection in invasive species.

[1]  Rousset.  Genetic differentiation between individuals , 1999 .

[2]  J. Tufto.  Effects of Releasing Maladapted Individuals: A Demographic‐Evolutionary Model , 2001, The American Naturalist.

[3]  M. Beaumont.  Detecting population expansion and decline using microsatellites. , 1999, Genetics.

[4]  H. Mooney,et al.  The evolutionary impact of invasive species , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[5]  T. C. Beebee,et al.  Population genetics of a successful invader: the marsh frog Rana ridibunda in Britain , 2003, Molecular ecology.

[6]  P. Donnelly,et al.  Inferring coalescence times from DNA sequence data. , 1997, Genetics.

[7]  David J. Balding,et al.  Inferences from DNA data: population histories, evolutionary processes and forensic match probabilities , 2003 .

[8]  W. Freeland,et al.  The rate of range expansion by Bufo marinus in Northern Australia, 1980-84 , 1985 .

[9]  P. Chesson,et al.  Community ecology theory as a framework for biological invasions , 2002 .

[10]  A. Estoup,et al.  Phylogeography and population structure of an ecotonal marsupial, Bettongia tropica, determined using mtDNA and microsatellites , 2000, Molecular ecology.

[11]  L Lacey Knowles,et al.  Statistical phylogeography. , 2002, Molecular ecology.

[12]  M. Nordborg,et al.  Coalescent Theory , 2019, Handbook of Statistical Genomics.

[13]  Templeton,et al.  Nested clade analyses of phylogeographic data: testing hypotheses about gene flow and population history , 1998, Molecular ecology.

[14]  D. Reznick,et al.  The population ecology of contemporary adaptations: what empirical studies reveal about the conditions that promote adaptive evolution , 2004, Genetica.

[15]  M. Nei.  Molecular Evolutionary Genetics , 1987 .

[16]  François Rousset,et al.  GENEPOP (version 1.2): population genetic software for exact tests and ecumenicism , 1995 .

[17]  B. Rannala,et al.  Detecting immigration by using multilocus genotypes. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[18]  M. Kirkpatrick,et al.  Evolution of a Species' Range , 1997, The American Naturalist.

[19]  M. Stephens,et al.  Inference Under the Coalescent , 2004 .

[20]  S. Easteal.  The history of introductions of Bufo marinus (Amphibia: Anura); a natural experiment in evolution , 1981 .

[21]  P. Grant.  Evolution on islands , 1998 .

[22]  S. Easteal,et al.  Continuing Geographical Spread of Bufo marinus in Australia: Range Expansion between 1974 and 1980 , 1985 .

[23]  C. Lee.  Evolutionary genetics of invasive species , 2002 .

[24]  T. Case,et al.  POPULATION GENETICS AND COLONY STRUCTURE OF THE ARGENTINE ANT (LINEPITHEMA HUMILE) IN ITS NATIVE AND INTRODUCED RANGES , 2001, Evolution; international journal of organic evolution.

[25]  Peter Beerli,et al.  Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[26]  A. Estoup,et al.  Bayesian inferences on the recent island colonization history by the bird Zosterops lateralis lateralis , 2003, Molecular ecology.

[27]  J. Felsenstein,et al.  Maximum-likelihood estimation of migration rates and effective population numbers in two populations using a coalescent approach. , 1999, Genetics.

[28]  Populations of Cane Toad, Bufo-Marinus, in Relation to Time Since Colonization , 1986 .

[29]  J. Cheverud,et al.  EPISTASIS AS A SOURCE OF INCREASED ADDITIVE GENETIC VARIANCE AT POPULATION BOTTLENECKS , 1996, Evolution; international journal of organic evolution.

[30]  R Gomulkiewicz,et al.  The effects of density dependence and immigration on local adaptation and niche evolution in a black-hole sink environment. , 1999, Theoretical population biology.

[31]  R. Griffiths,et al.  Inference from gene trees in a subdivided population. , 2000, Theoretical population biology.

[32]  Ross Ihaka,et al.  Gentleman R: R: A language for data analysis and graphics , 1996 .

[33]  Jon A Yamato,et al.  Maximum likelihood estimation of population growth rates based on the coalescent. , 1998, Genetics.

[34]  Wen-Hsiung Li,et al.  Coalescing into the 21st century: An overview and prospects of coalescent theory. , 1999, Theoretical population biology.

[35]  R. Alford,et al.  Nomadic movement in tropical toads , 2002 .

[36]  S. Wright,et al.  Evolution in Mendelian Populations. , 1931, Genetics.

[37]  POPULATION GENETICS AND COLONY STRUCTURE OF THE ARGENTINE ANT (LINEPITHEMA HUMILE) IN ITS NATIVE AND INTRODUCED RANGES , 2001, Evolution; international journal of organic evolution.

[38]  D. Balding,et al.  Handbook of statistical genetics , 2004 .

[39]  S. Easteal.  THE ECOLOGICAL GENETICS OF INTRODUCED POPULATIONS OF THE GIANT TOAD, BUFO MARINUS. III. GEOGRAPHICAL PATTERNS OF VARIATION , 1985, Evolution; international journal of organic evolution.

[40]  W. Boughton,et al.  Expansion of the Range of the Introduced Toad Bufo marinus in Australia from 1935 to 1974 , 1981 .

[41]  J. Cornuet,et al.  Inferring population history from microsatellite and enzyme data in serially introduced cane toads, Bufo marinus. , 2001, Genetics.

[42]  Thomas Lenormand,et al.  Gene flow and the limits to natural selection , 2002 .

[43]  M. Feldman,et al.  Population growth of human Y chromosomes: a study of Y chromosome microsatellites. , 1999, Molecular biology and evolution.

[44]  R. Freckleton,et al.  The Ecological Detective: Confronting Models with Data , 1999 .

[45]  T. Burke,et al.  Effective number of breeding adults in Bufo bufo estimated from age‐specific variation at minisatellite loci , 1997, Molecular ecology.

[46]  B. Weir,et al.  ESTIMATING F‐STATISTICS FOR THE ANALYSIS OF POPULATION STRUCTURE , 1984, Evolution; international journal of organic evolution.

[47]  F. Rousset,et al.  Absence of evidence for isolation by distance in an expanding cane toad (Bufo marinus) population: an individual‐based analysis of microsatellite genotypes , 2000, Molecular ecology.

[48]  D. Balding,et al.  Genealogical inference from microsatellite data. , 1998, Genetics.

[49]  E. Beurden.  An Isolated and Expanding Population of the Introduced Toad Bufo marinus in New South Wales , 1980 .

[50]  A. Estoup,et al.  Polymerase chain reaction primers for polymorphic microsatellite loci in the invasive toad species Bufo marinus , 2000, Molecular ecology.

[51]  R. Sutherst,et al.  The Potential Geographical Distribution of the Cane Toad, Bufo marinus L. in Australia , 1996 .

[52]  P Donnelly,et al.  Coalescents and genealogical structure under neutrality. , 1995, Annual review of genetics.

[53]  François Rousset,et al.  Influence of mutational and sampling factors on the estimation of demographic parameters in a "continuous" population under isolation by distance. , 2003, Molecular biology and evolution.

[54]  S. Easteal,et al.  The ecological genetics of introduced populations of the giant toad, Bufo marinus (Amphibia: Anura): dispersal and neighbourhood size , 1986 .

[55]  P. Donnelly,et al.  Inference in molecular population genetics , 2000 .

[56]  R. Nielsen,et al.  Distinguishing migration from isolation: a Markov chain Monte Carlo approach. , 2001, Genetics.

[57]  Range expansion of the introduced cane toad Bufo marinus in New South Wales , 1991 .

[58]  H. Chiang,et al.  The biology of invasions. , 1977 .

[59]  C. Loader.  Local Likelihood Density Estimation , 1996 .

[60]  M. Beaumont,et al.  Estimation of admixture proportions: a likelihood-based approach using Markov chain Monte Carlo. , 2001, Genetics.

[61]  Louis Legendre,et al.  The Importance of Being Digital , 1963 .

[62]  M. Crawley,et al.  Exotic plant invasions and the enemy release hypothesis , 2002 .

[63]  G. Zug,et al.  The marine toad, Bufo marinus : a natural history resumé of native populations , 1979 .

[64]  D. Balding,et al.  Approximate Bayesian computation in population genetics. , 2002, Genetics.

[65]  J M Cornuet,et al.  Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. , 1996, Genetics.

[66]  Montgomery Slatkin,et al.  ISOLATION BY DISTANCE IN EQUILIBRIUM AND NON‐EQUILIBRIUM POPULATIONS , 1993, Evolution; international journal of organic evolution.