@inbook{3277, abstract = {The problem of the origin of metazoa is becoming more urgent in the context of astrobiology. By now it is clear that clues to the understanding of this crucial transition in the evolution of life can arise in a fourth pathway besides the three possibilities in the quest for simplicity outlined by Bonner in his classical book. In other words, solar system exploration seems to be one way in the long-term to elucidate the simplicity of evolutionary development. We place these ideas in the context of different inheritance systems, namely the genotypic and phenotypic replicators with limited or unlimited heredity, and ask which of these can support multicellular development, and to which degree of complexity. However, the quest for evidence on the evolution of biotas from planets around other stars does not seem to be feasible with present technology with direct visualization of living organisms on exoplanets. But this may be attempted on the Galilean moons of Jupiter where there is a possibility of detecting reliable biomarkers in the next decade with the Europa Jupiter System Mission, in view of recent progress by landing micropenetrators on planetary, or satellite surfaces. Mars is a second possibility in the inner Solar System, in spite of the multiple difficulties faced by the fleet of past, present and future missions. We discuss a series of preliminary ideas for elucidating the origin of metazoan analogues with available instrumentation in potential payloads of feasible space missions to the Galilean moons.}, author = {de Vladar, Harold and Chela Flores, Julian}, booktitle = {Life on Earth and other planetary bodies}, pages = {387 -- 405}, publisher = {Springer}, title = {{Can the evolution of multicellularity be anticipated in the exploration of the solar system?}}, doi = {10.1007/978-94-007-4966-5_22}, volume = {24}, year = {2012}, } @article{498, abstract = {Understanding patterns and correlates of local adaptation in heterogeneous landscapes can provide important information in the selection of appropriate seed sources for restoration. We assessed the extent of local adaptation of fitness components in 12 population pairs of the perennial herb Rutidosis leptorrhynchoides (Asteraceae) and examined whether spatial scale (0.7-600 km), environmental distance, quantitative (QST) and neutral (FST) genetic differentiation, and size of the local and foreign populations could predict patterns of adaptive differentiation. Local adaptation varied among populations and fitness components. Including all population pairs, local adaptation was observed for seedling survival, but not for biomass, while foreign genotype advantage was observed for reproduction (number of inflorescences). Among population pairs, local adaptation increased with QST and local population size for biomass. QST was associated with environmental distance, suggesting ecological selection for phenotypic divergence. However, low FST and variation in population structure in small populations demonstrates the interaction of gene flow and drift in constraining local adaptation in R. leptorrhynchoides. Our study indicates that for species in heterogeneous landscapes, collecting seed from large populations from similar environments to candidate sites is likely to provide the most appropriate seed sources for restoration.}, author = {Pickup, Melinda and Field, David and Rowell, David and Young, Andrew}, journal = {Evolutionary Applications}, number = {8}, pages = {913 -- 924}, publisher = {Wiley-Blackwell}, title = {{Predicting local adaptation in fragmented plant populations: Implications for restoration genetics}}, doi = {10.1111/j.1752-4571.2012.00284.x}, volume = {5}, year = {2012}, } @misc{13075, abstract = {Little is known about the stability of trophic relationships in complex natural communities over evolutionary timescales. Here, we use sequence data from 18 nuclear loci to reconstruct and compare the intraspecific histories of major Pleistocene refugial populations in the Middle East, the Balkans and Iberia in a guild of four Chalcid parasitoids (Cecidostiba fungosa, C. semifascia, Hobbya stenonota and Mesopolobus amaenus) all attacking Cynipid oak galls. We develop a likelihood method to numerically estimate models of divergence between three populations from multilocus data. We investigate the power of this framework on simulated data, and - using triplet alignments of intronic loci - quantify the support for all possible divergence relationships between refugial populations in the four parasitoids. Although an East to West order of population divergence has highest support in all but one species, we cannot rule out alternative population tree topologies. Comparing the estimated times of population splits between species, we find that one species, M. amaenus, has a significantly older history than the rest of the guild and must have arrived in central Europe at least one glacial cycle prior to other guild members. This suggests that although all four species may share a common origin in the East, they expanded westwards into Europe at different times.}, author = {Lohse, Konrad and Barton, Nicholas H and Stone, Graham and Melika, George}, publisher = {Dryad}, title = {{Data from: A likelihood-based comparison of population histories in a parasitoid guild}}, doi = {10.5061/DRYAD.0G0FS}, year = {2012}, } @misc{9758, abstract = {We propose a two-step procedure for estimating multiple migration rates in an approximate Bayesian computation (ABC) framework, accounting for global nuisance parameters. The approach is not limited to migration, but generally of interest for inference problems with multiple parameters and a modular structure (e.g. independent sets of demes or loci). We condition on a known, but complex demographic model of a spatially subdivided population, motivated by the reintroduction of Alpine ibex (Capra ibex) into Switzerland. In the first step, the global parameters ancestral mutation rate and male mating skew have been estimated for the whole population in Aeschbacher et al. (Genetics 2012; 192: 1027). In the second step, we estimate in this study the migration rates independently for clusters of demes putatively connected by migration. For large clusters (many migration rates), ABC faces the problem of too many summary statistics. We therefore assess by simulation if estimation per pair of demes is a valid alternative. We find that the trade-off between reduced dimensionality for the pairwise estimation on the one hand and lower accuracy due to the assumption of pairwise independence on the other depends on the number of migration rates to be inferred: the accuracy of the pairwise approach increases with the number of parameters, relative to the joint estimation approach. To distinguish between low and zero migration, we perform ABC-type model comparison between a model with migration and one without. Applying the approach to microsatellite data from Alpine ibex, we find no evidence for substantial gene flow via migration, except for one pair of demes in one direction.}, author = {Aeschbacher, Simon and Futschik, Andreas and Beaumont, Mark}, publisher = {Dryad}, title = {{Data from: Approximate Bayesian computation for modular inference problems with many parameters: the example of migration rates}}, doi = {10.5061/dryad.274b1}, year = {2012}, } @article{3290, abstract = {Analysis of genomic data requires an efficient way to calculate likelihoods across very large numbers of loci. We describe a general method for finding the distribution of genealogies: we allow migration between demes, splitting of demes [as in the isolation-with-migration (IM) model], and recombination between linked loci. These processes are described by a set of linear recursions for the generating function of branch lengths. Under the infinite-sites model, the probability of any configuration of mutations can be found by differentiating this generating function. Such calculations are feasible for small numbers of sampled genomes: as an example, we show how the generating function can be derived explicitly for three genes under the two-deme IM model. This derivation is done automatically, using Mathematica. Given data from a large number of unlinked and nonrecombining blocks of sequence, these results can be used to find maximum-likelihood estimates of model parameters by tabulating the probabilities of all relevant mutational configurations and then multiplying across loci. The feasibility of the method is demonstrated by applying it to simulated data and to a data set previously analyzed by Wang and Hey (2010) consisting of 26,141 loci sampled from Drosophila simulans and D. melanogaster. Our results suggest that such likelihood calculations are scalable to genomic data as long as the numbers of sampled individuals and mutations per sequence block are small.}, author = {Lohse, Konrad and Harrison, Richard and Barton, Nicholas H}, journal = {Genetics}, number = {3}, pages = {977 -- 987}, publisher = {Genetics Society of America}, title = {{A general method for calculating likelihoods under the coalescent process}}, doi = {10.1534/genetics.111.129569}, volume = {189}, year = {2011}, } @misc{9762, abstract = {Defining population structure and genetic diversity levels is of the utmost importance for developing efficient conservation strategies. Overfishing has caused mean annual catches of the European spiny lobster (Palinurus elephas) to decrease alarmingly along its distribution area. In this context, there is a need for comprehensive studies to evaluate the genetic health of the exploited populations. The present work is based on a set of 10 nuclear markers amplified in 331 individuals from 10 different localities covering most of P. elephas distribution area. Samples from Atlantic and Mediterranean basins showed small but significant differences, indicating that P. elephas populations do not behave as a single panmictic unit but form two partially-overlapping groups. Despite intense overfishing, our dataset did not recover a recent bottleneck signal, and showed a large and stable historical effective size instead. This result could be accounted for by specific life history traits (reproduction and longevity) and the limitations of molecular markers in covering very recent timescales for non temporal samples. Our study emphasizes the necessity of integrating information on effective population sizes and life history parameters when evaluating population connectivity levels from genetic data.}, author = {Palero, Ferran and Abello, Pere and Macpherson, Enrique and Beaumont, Mark and Pascual, Marta}, publisher = {IST Austria}, title = {{Data from: Effect of oceanographic barriers and overfishing on the population genetic structure of the European spiny lobster (Palinurus elephas)}}, doi = {10.5061/dryad.299h8}, year = {2011}, } @article{3372, abstract = {Nowak et al.1 argue that inclusive fitness theory has been of little value in explaining the natural world, and that it has led to negligible progress in explaining the evolution of eusociality. However, we believe that their arguments are based upon a misunderstanding of evolutionary theory and a misrepresentation of the empirical literature. We will focus our comments on three general issues.}, author = {Abbot, Patrick and Abe, Jun and Alcock, John and Alizon, Samuel and Alpedrinha, Joao and Andersson, Malte and Andre, Jean and Van Baalen, Minus and Balloux, Francois and Balshine, Sigal and Barton, Nicholas H and Beukeboom, Leo and Biernaskie, Jay and Bilde, Trine and Borgia, Gerald and Breed, Michael and Brown, Sam and Bshary, Redouan and Buckling, Angus and Burley, Nancy and Burton Chellew, Max and Cant, Michael and Chapuisat, Michel and Charnov, Eric and Clutton Brock, Tim and Cockburn, Andrew and Cole, Blaine and Colegrave, Nick and Cosmides, Leda and Couzin, Iain and Coyne, Jerry and Creel, Scott and Crespi, Bernard and Curry, Robert and Dall, Sasha and Day, Troy and Dickinson, Janis and Dugatkin, Lee and El Mouden, Claire and Emlen, Stephen and Evans, Jay and Ferriere, Regis and Field, Jeremy and Foitzik, Susanne and Foster, Kevin and Foster, William and Fox, Charles and Gadau, Juergen and Gandon, Sylvain and Gardner, Andy and Gardner, Michael and Getty, Thomas and Goodisman, Michael and Grafen, Alan and Grosberg, Rick and Grozinger, Christina and Gouyon, Pierre and Gwynne, Darryl and Harvey, Paul and Hatchwell, Ben and Heinze, Jürgen and Helantera, Heikki and Helms, Ken and Hill, Kim and Jiricny, Natalie and Johnstone, Rufus and Kacelnik, Alex and Kiers, E Toby and Kokko, Hanna and Komdeur, Jan and Korb, Judith and Kronauer, Daniel and Kümmerli, Rolf and Lehmann, Laurent and Linksvayer, Timothy and Lion, Sébastien and Lyon, Bruce and Marshall, James and Mcelreath, Richard and Michalakis, Yannis and Michod, Richard and Mock, Douglas and Monnin, Thibaud and Montgomerie, Robert and Moore, Allen and Mueller, Ulrich and Noë, Ronald and Okasha, Samir and Pamilo, Pekka and Parker, Geoff and Pedersen, Jes and Pen, Ido and Pfennig, David and Queller, David and Rankin, Daniel and Reece, Sarah and Reeve, Hudson and Reuter, Max and Roberts, Gilbert and Robson, Simon and Roze, Denis and Rousset, Francois and Rueppell, Olav and Sachs, Joel and Santorelli, Lorenzo and Schmid Hempel, Paul and Schwarz, Michael and Scott Phillips, Tom and Shellmann Sherman, Janet and Sherman, Paul and Shuker, David and Smith, Jeff and Spagna, Joseph and Strassmann, Beverly and Suarez, Andrew and Sundström, Liselotte and Taborsky, Michael and Taylor, Peter and Thompson, Graham and Tooby, John and Tsutsui, Neil and Tsuji, Kazuki and Turillazzi, Stefano and Úbeda, Francisco and Vargo, Edward and Voelkl, Bernard and Wenseleers, Tom and West, Stuart and West Eberhard, Mary and Westneat, David and Wiernasz, Diane and Wild, Geoff and Wrangham, Richard and Young, Andrew and Zeh, David and Zeh, Jeanne and Zink, Andrew}, journal = {Nature}, number = {7339}, pages = {E1 -- E4}, publisher = {Nature Publishing Group}, title = {{Inclusive fitness theory and eusociality}}, doi = {10.1038/nature09831}, volume = {471}, year = {2011}, } @article{3375, abstract = {By exploiting an analogy between population genetics and statistical mechanics, we study the evolution of a polygenic trait under stabilizing selection, mutation and genetic drift. This requires us to track only four macroscopic variables, instead of the distribution of all the allele frequencies that influence the trait. These macroscopic variables are the expectations of: the trait mean and its square, the genetic variance, and of a measure of heterozygosity, and are derived from a generating function that is in turn derived by maximizing an entropy measure. These four macroscopics are enough to accurately describe the dynamics of the trait mean and of its genetic variance (and in principle of any other quantity). Unlike previous approaches that were based on an infinite series of moments or cumulants, which had to be truncated arbitrarily, our calculations provide a well-defined approximation procedure. We apply the framework to abrupt and gradual changes in the optimum, as well as to changes in the strength of stabilizing selection. Our approximations are surprisingly accurate, even for systems with as few as five loci. We find that when the effects of drift are included, the expected genetic variance is hardly altered by directional selection, even though it fluctuates in any particular instance. We also find hysteresis, showing that even after averaging over the microscopic variables, the macroscopic trajectories retain a memory of the underlying genetic states.}, author = {de Vladar, Harold and Barton, Nicholas H}, journal = {Journal of the Royal Society Interface}, number = {58}, pages = {720 -- 739}, publisher = {The Royal Society}, title = {{The statistical mechanics of a polygenic character under stabilizing selection mutation and drift}}, doi = {10.1098/rsif.2010.0438}, volume = {8}, year = {2011}, } @article{3380, abstract = {Linkage between markers and genes that affect a phenotype of interest may be determined by examining differences in marker allele frequency in the extreme progeny of a cross between two inbred lines. This strategy is usually employed when pooling is used to reduce genotyping costs. When the cross progeny are asexual, the extreme progeny may be selected by multiple generations of asexual reproduction and selection. We analyse this method of measuring phenotype in asexual progeny and examine the changes in marker allele frequency due to selection over many generations. Stochasticity in marker frequency in the selected population arises due to the finite initial population size. We derive the distribution of marker frequency as a result of selection at a single major locus, and show that in order to avoid spurious changes in marker allele frequency in the selected population, the initial population size should be in the low to mid hundreds.}, author = {Logeswaran, Sayanthan and Barton, Nicholas H}, journal = {Genetical Research}, number = {3}, pages = {221 -- 232}, publisher = {Cambridge University Press}, title = {{Mapping Mendelian traits in asexual progeny using changes in marker allele frequency}}, doi = {10.1017/S0016672311000115}, volume = {93}, year = {2011}, } @article{3390, abstract = {What determines the genetic contribution that an individual makes to future generations? With biparental reproduction, each individual leaves a 'pedigree' of descendants, determined by the biparental relationships in the population. The pedigree of an individual constrains the lines of descent of each of its genes. An individual's reproductive value is the expected number of copies of each of its genes that is passed on to distant generations conditional on its pedigree. For the simplest model of biparental reproduction analogous to the Wright-Fisher model, an individual's reproductive value is determined within ~10 generations, independent of population size. Partial selfing and subdivision do not greatly slow this convergence. Our central result is that the probability that a gene will survive is proportional to the reproductive value of the individual that carries it, and that conditional on survival, after a few tens of generations, the distribution of the number of surviving copies is the same for all individuals, whatever their reproductive value. These results can be generalized to the joint distribution of surviving blocks of ancestral genome. Selection on unlinked loci in the genetic background may greatly increase the variance in reproductive value, but the above results nevertheless still hold. The almost linear relationship between survival probability and reproductive value also holds for weakly favored alleles. Thus, the influence of the complex pedigree of descendants on an individual's genetic contribution to the population can be summarized through a single number: its reproductive value.}, author = {Barton, Nicholas H and Etheridge, Alison}, journal = {Genetics}, number = {4}, pages = {953 -- 973}, publisher = {Genetics Society of America}, title = {{The relation between reproductive value and genetic contribution}}, doi = {10.1534/genetics.111.127555}, volume = {188}, year = {2011}, } @article{3391, abstract = {Evolutionary biology shares many concepts with statistical physics: both deal with populations, whether of molecules or organisms, and both seek to simplify evolution in very many dimensions. Often, methodologies have undergone parallel and independent development, as with stochastic methods in population genetics. Here, we discuss aspects of population genetics that have embraced methods from physics: non-equilibrium statistical mechanics, travelling waves and Monte-Carlo methods, among others, have been used to study polygenic evolution, rates of adaptation and range expansions. These applications indicate that evolutionary biology can further benefit from interactions with other areas of statistical physics; for example, by following the distribution of paths taken by a population through time}, author = {de Vladar, Harold and Barton, Nicholas H}, journal = {Trends in Ecology and Evolution}, number = {8}, pages = {424 -- 432}, publisher = {Cell Press}, title = {{The contribution of statistical physics to evolutionary biology}}, doi = {10.1016/j.tree.2011.04.002}, volume = {26}, year = {2011}, } @article{3393, abstract = {Unlike unconditionally advantageous “Fisherian” variants that tend to spread throughout a species range once introduced anywhere, “bistable” variants, such as chromosome translocations, have two alternative stable frequencies, absence and (near) fixation. Analogous to populations with Allee effects, bistable variants tend to increase locally only once they become sufficiently common, and their spread depends on their rate of increase averaged over all frequencies. Several proposed manipulations of insect populations, such as using Wolbachia or “engineered underdominance” to suppress vector-borne diseases, produce bistable rather than Fisherian dynamics. We synthesize and extend theoretical analyses concerning three features of their spatial behavior: rate of spread, conditions to initiate spread from a localized introduction, and wave stopping caused by variation in population densities or dispersal rates. Unlike Fisherian variants, bistable variants tend to spread spatially only for particular parameter combinations and initial conditions. Wave initiation requires introduction over an extended region, while subsequent spatial spread is slower than for Fisherian waves and can easily be halted by local spatial inhomogeneities. We present several new results, including robust sufficient conditions to initiate (and stop) spread, using a one-parameter cubic approximation applicable to several models. The results have both basic and applied implications.}, author = {Barton, Nicholas H and Turelli, Michael}, issn = {1537-5323}, journal = {American Naturalist}, number = {3}, pages = {E48 -- E75}, publisher = {The University of Chicago Press}, title = {{Spatial waves of advance with bistable dynamics: Cytoplasmic and genetic analogues of Allee effects}}, doi = {10.1086/661246}, volume = {178}, year = {2011}, } @article{3394, abstract = {Random genetic drift shifts clines in space, alters their width, and distorts their shape. Such random fluctuations complicate inferences from cline width and position. Notably, the effect of genetic drift on the expected shape of the cline is opposite to the naive (but quite common) misinterpretation of classic results on the expected cline. While random drift on average broadens the overall cline in expected allele frequency, it narrows the width of any particular cline. The opposing effects arise because locally, drift drives alleles to fixation—but fluctuations in position widen the expected cline. The effect of genetic drift can be predicted from standardized variance in allele frequencies, averaged across the habitat: 〈F〉. A cline maintained by spatially varying selection (step change) is expected to be narrower by a factor of relative to the cline in the absence of drift. The expected cline is broader by the inverse of this factor. In a tension zone maintained by underdominance, the expected cline width is narrower by about 1 – 〈F〉relative to the width in the absence of drift. Individual clines can differ substantially from the expectation, and we give quantitative predictions for the variance in cline position and width. The predictions apply to clines in almost one-dimensional circumstances such as hybrid zones in rivers, deep valleys, or along a coast line and give a guide to what patterns to expect in two dimensions.}, author = {Polechova, Jitka and Barton, Nicholas H}, journal = {Genetics}, number = {1}, pages = {227 -- 235}, publisher = {Genetics Society of America}, title = {{Genetic drift widens the expected cline but narrows the expected cline width}}, doi = {10.1534/genetics.111.129817}, volume = {189}, year = {2011}, } @article{3395, abstract = {Defining population structure and genetic diversity levels is of the utmost importance for developing efficient conservation strategies. Overfishing has caused mean annual catches of the European spiny lobster (Palinurus elephas) to decrease alarmingly along its distribution area. In this context, there is a need for comprehensive studies aiming to evaluate the genetic health of the exploited populations. The present study is based on a set of ten nuclear markers amplified in 331 individuals from ten different localities covering most of P. elephas distribution area. Samples from Atlantic and Mediterranean basins showed small but significant differences, indicating that P. elephas populations do not behave as a single panmictic unit but form two partially-overlapping groups. Despite intense overfishing, our dataset did not recover a recent bottleneck signal, and instead showed a large and stable historical effective size. This result could be accounted for by specific life-history traits (reproduction and longevity) and the limitations of molecular markers in covering recent timescales for nontemporal samples. The findings of the present study emphasize the need to integrate information on effective population sizes and life-history parameters when evaluating population connectivity levels from genetic data.}, author = {Palero, Ferran and Abello, Pere and Macpherson, Enrique and Beaumont, Mark and Pascual, Marta}, journal = {Biological Journal of the Linnean Society}, number = {2}, pages = {407 -- 418}, publisher = {Wiley-Blackwell}, title = {{Effect of oceanographic barriers and overfishing on the population genetic structure of the European spiny lobster Palinurus elephas}}, doi = {10.1111/j.1095-8312.2011.01728.x}, volume = {104}, year = {2011}, } @article{3778, author = {Barton, Nicholas H}, journal = {Heredity}, number = {2}, pages = {205 -- 206}, publisher = {Nature Publishing Group}, title = {{Estimating linkage disequilibria}}, doi = {10.1038/hdy.2010.67}, volume = {106}, year = {2011}, } @article{3784, abstract = {Advanced stages of Scyllarus phyllosoma larvae were collected by demersal trawling during fishery research surveys in the western Mediterranean Sea in 2003–2005. Nucleotide sequence analysis of the mitochondrial 16S rDNA gene allowed the final-stage phyllosoma of Scyllarus arctus to be identified among these larvae. Its morphology is described and illustrated. This constitutes the second complete description of a Scyllaridae phyllosoma with its specific identity being validated by molecular techniques (the first was S. pygmaeus). These results also solved a long lasting taxonomic anomaly of several species assigned to the ancient genus Phyllosoma Leach, 1814. Detailed examination indicated that the final-stage phyllosoma of S. arctus shows closer affinities with the American scyllarid Scyllarus depressus or with the Australian Scyllarus sp. b (sensu Phillips et al., 1981) than to its sympatric species S. pygmaeus.}, author = {Palero, Ferran and Guerao, Guillermo and Clark, Paul and Abello, Pere}, journal = {Journal of the Marine Biological Association of the United Kingdom}, number = {2}, pages = {485 -- 492}, publisher = {Cambridge University Press}, title = {{Scyllarus arctus (Crustacea: Decapoda: Scyllaridae) final stage phyllosoma identified by DNA analysis, with morphological description}}, doi = {10.1017/S0025315410000287}, volume = {91}, year = {2011}, } @article{3303, abstract = {Biological traits result in part from interactions between different genetic loci. This can lead to sign epistasis, in which a beneficial adaptation involves a combination of individually deleterious or neutral mutations; in this case, a population must cross a “fitness valley” to adapt. Recombination can assist this process by combining mutations from different individuals or retard it by breaking up the adaptive combination. Here, we analyze the simplest fitness valley, in which an adaptation requires one mutation at each of two loci to provide a fitness benefit. We present a theoretical analysis of the effect of recombination on the valley-crossing process across the full spectrum of possible parameter regimes. We find that low recombination rates can speed up valley crossing relative to the asexual case, while higher recombination rates slow down valley crossing, with the transition between the two regimes occurring when the recombination rate between the loci is approximately equal to the selective advantage provided by the adaptation. In large populations, if the recombination rate is high and selection against single mutants is substantial, the time to cross the valley grows exponentially with population size, effectively meaning that the population cannot acquire the adaptation. Recombination at the optimal (low) rate can reduce the valley-crossing time by up to several orders of magnitude relative to that in an asexual population. }, author = {Weissman, Daniel and Feldman, Marcus and Fisher, Daniel}, journal = {Genetics}, number = {4}, pages = {1389 -- 1410}, publisher = {Genetics Society of America}, title = {{The rate of fitness-valley crossing in sexual populations}}, doi = {10.1534/genetics.110.123240}, volume = {186}, year = {2010}, } @article{474, abstract = {Classical models of gene flow fail in three ways: they cannot explain large-scale patterns; they predict much more genetic diversity than is observed; and they assume that loosely linked genetic loci evolve independently. We propose a new model that deals with these problems. Extinction events kill some fraction of individuals in a region. These are replaced by offspring from a small number of parents, drawn from the preexisting population. This model of evolution forwards in time corresponds to a backwards model, in which ancestral lineages jump to a new location if they are hit by an event, and may coalesce with other lineages that are hit by the same event. We derive an expression for the identity in allelic state, and show that, over scales much larger than the largest event, this converges to the classical value derived by Wright and Malécot. However, rare events that cover large areas cause low genetic diversity, large-scale patterns, and correlations in ancestry between unlinked loci.}, author = {Barton, Nicholas H and Kelleher, Jerome and Etheridge, Alison}, journal = {Evolution}, number = {9}, pages = {2701 -- 2715}, publisher = {Wiley-Blackwell}, title = {{A new model for extinction and recolonization in two dimensions: Quantifying phylogeography}}, doi = {10.1111/j.1558-5646.2010.01019.x}, volume = {64}, year = {2010}, } @misc{9764, author = {Rosas, Ulises and Barton, Nicholas H and Copsey, Lucy and Barbier De Reuille, Pierre and Coen, Enrico}, publisher = {Public Library of Science}, title = {{Heterosis and the drift load}}, doi = {10.1371/journal.pbio.1000429.s003}, year = {2010}, } @article{3604, abstract = {We investigated temporal changes in hybridization and introgression between native red deer (Cervus elaphus) and invasive Japanese sika (Cervus nippon) on the Kintyre Peninsula, Scotland, over 15 years, through analysis of 1513 samples of deer at 20 microsatellite loci and a mtDNA marker. We found no evidence that either the proportion of recent hybrids, or the levels of introgression had changed over the study period. Nevertheless, in one population where the two species have been in contact since ∼1970, 44% of individuals sampled during the study were hybrids. This suggests that hybridization between these species can proceed fairly rapidly. By analysing the number of alleles that have introgressed from polymorphic red deer into the genetically homogenous sika population, we reconstructed the haplotypes of red deer alleles introduced by backcrossing. Five separate hybridization events could account for all the recently hybridized sika-like individuals found across a large section of the Peninsula. Although we demonstrate that low rates of F1 hybridization can lead to substantial introgression, the progress of hybridization and introgression appears to be unpredictable over the short timescales.}, author = {Senn, Helen and Goodman, Simon and Swanson, Graeme and Barton, Nicholas H and Pemberton, Josephine}, journal = {Molecular Ecology}, number = {5}, pages = {910 -- 924}, publisher = {Wiley-Blackwell}, title = {{Investigating temporal changes in hybridisation and introgression between invasive sika (Cervus nippon) and native red deer (Cervus elaphus) on the Kintyre Peninsula, Scotland}}, doi = {10.1111/j.1365-294X.2009.04497.x}, volume = {19}, year = {2010}, }