COMBINING GENETIC AND ECOLOGICAL DATA TO ESTIMATE SEA TURTLE ORIGINS

Many species of sea turtles spend part of their life cycle gathered in large feeding aggregations that combine individuals from widely separated rookery populations. Biologists have applied methods of mixed-stock analysis to mitochondrial DNA samples from rookeries and mixed populations to estimate the contributions of different rookeries to the mixed stock. These methods are limited by the amount of genetic overlap between rookeries and fail to incorporate ecological covariates such as rookery size and location within major ocean currents that are strongly suspected to affect rookery contributions. A new hierarchical Bayesian model for rookery contributions incorporates these covariates (and potentially others) to draw stronger conclusions from existing data. Applying the model to various simple scenarios shows that, in some cases, it can accurately estimate turtle origins even when turtles come from rookeries with high degrees of genetic overlap. Applying the model to more complex simulations shows that it performs well in a wide range of scenarios. Applying the model to existing data on green turtles ( Chelonia mydas) narrows confidence intervals but does not change point estimates significantly. Applying it to loggerhead turtles (Caretta caretta) strengthens the dominance of the large rookery in south Florida, and brings estimates from a small data set on sea turtle strandings into line with those from rookery data. Used appropriately, hierarchical Bayesian methods offer great potential for introducing multiple levels of variation and ecological covariates into ecological models.

[1]  David R. Anderson,et al.  Bayesian Methods in Cosmology: Model selection and multi-model inference , 2009 .

[2]  W. Link,et al.  HIERARCHICAL MODELING OF POPULATION STABILITY AND SPECIES GROUP ATTRIBUTES FROM SURVEY DATA , 2002 .

[3]  J. Pella,et al.  Bayesian methods for analysis of stock mixtures from genetic characters , 2001 .

[4]  C. McCulloch,et al.  Generalized Linear Mixed Models , 2005 .

[5]  Shizhong Xu,et al.  Constrained Least Squares Estimation of Mixed Population Stock Composition from mtDNA Haplotype Frequency Data , 1994 .

[6]  John Aitchison,et al.  The Statistical Analysis of Compositional Data , 1986 .

[7]  Sylvia Richardson,et al.  Markov Chain Monte Carlo in Practice , 1997 .

[8]  M. Miyamoto,et al.  Genetic composition of a green turtle (Chelonia mydas) feeding ground population: evidence for multiple origins , 1998 .

[9]  Bradley P. Carlin,et al.  BAYES AND EMPIRICAL BAYES METHODS FOR DATA ANALYSIS , 1996, Stat. Comput..

[10]  K. Bjorndal,et al.  TRANSATLANTIC DEVELOPMENTAL MIGRATIONS OF LOGGERHEAD SEA TURTLES DEMONSTRATED BY mtDNA SEQUENCE ANALYSIS , 1998 .

[11]  P. Guttorp,et al.  Statistical Interpretation of Species Composition , 2001 .

[12]  J. J..  Marine Biology , 1914, Nature.

[13]  T. Beacham,et al.  Estimating Stock Composition in Mixed Stock Fisheries Using Morphometric, Meristic, and Electrophoretic Characteristics , 1984 .

[14]  Bradley P. Carlin,et al.  Bayesian measures of model complexity and fit , 2002 .

[15]  V. Carey,et al.  Mixed-Effects Models in S and S-Plus , 2001 .

[16]  N. Ryman,et al.  Population Genetics and Fishery Management , 1987 .

[17]  J. Aitchison.  On criteria for measures of compositional difference , 1992 .

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

[19]  Brian Dennis,et al.  Discussion: Should Ecologists Become Bayesians? , 1996 .

[20]  A. R. Gallant,et al.  Noise and Nonlinearity in Measles Epidemics: Combining Mechanistic and Statistical Approaches to Population Modeling , 1998, The American Naturalist.

[21]  M. Miyamoto,et al.  Origin of Hawksbill Turtles in a Caribbean Feeding Area as Indicated by Genetic Markers , 1996 .

[22]  H. Blankenship,et al.  Genetic Diversity Patterns of Chum Salmon in the Pacific Northwest , 1994 .

[23]  C. Walters,et al.  Are age-structured models appropriate for catch-effort data? , 1985 .

[24]  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.

[25]  B. Bolker,et al.  SEA TURTLE STOCK ESTIMATION USING GENETIC MARKERS: ACCOUNTING FOR SAMPLING ERROR OF RARE GENOTYPES , 2003 .