TY - JOUR AB - We analyse how migration from a large mainland influences genetic load and population numbers on an island, in a scenario where fitness-affecting variants are unconditionally deleterious, and where numbers decline with increasing load. Our analysis shows that migration can have qualitatively different effects, depending on the total mutation target and fitness effects of deleterious variants. In particular, we find that populations exhibit a genetic Allee effect across a wide range of parameter combinations, when variants are partially recessive, cycling between low-load (large-population) and high-load (sink) states. Increased migration reduces load in the sink state (by increasing heterozygosity) but further inflates load in the large-population state (by hindering purging). We identify various critical parameter thresholds at which one or other stable state collapses, and discuss how these thresholds are influenced by the genetic versus demographic effects of migration. Our analysis is based on a ‘semi-deterministic’ analysis, which accounts for genetic drift but neglects demographic stochasticity. We also compare against simulations which account for both demographic stochasticity and drift. Our results clarify the importance of gene flow as a key determinant of extinction risk in peripheral populations, even in the absence of ecological gradients. This article is part of the theme issue ‘Species’ ranges in the face of changing environments (part I)’. AU - Sachdeva, Himani AU - Olusanya, Oluwafunmilola O AU - Barton, Nicholas H ID - 10658 IS - 1846 JF - Philosophical Transactions of the Royal Society B SN - 0962-8436 TI - Genetic load and extinction in peripheral populations: The roles of migration, drift and demographic stochasticity VL - 377 ER - TY - JOUR AB - Many studies have quantified the distribution of heterozygosity and relatedness in natural populations, but few have examined the demographic processes driving these patterns. In this study, we take a novel approach by studying how population structure affects both pairwise identity and the distribution of heterozygosity in a natural population of the self-incompatible plant Antirrhinum majus. Excess variance in heterozygosity between individuals is due to identity disequilibrium, which reflects the variance in inbreeding between individuals; it is measured by the statistic g2. We calculated g2 together with FST and pairwise relatedness (Fij) using 91 SNPs in 22,353 individuals collected over 11 years. We find that pairwise Fij declines rapidly over short spatial scales, and the excess variance in heterozygosity between individuals reflects significant variation in inbreeding. Additionally, we detect an excess of individuals with around half the average heterozygosity, indicating either selfing or matings between close relatives. We use 2 types of simulation to ask whether variation in heterozygosity is consistent with fine-scale spatial population structure. First, by simulating offspring using parents drawn from a range of spatial scales, we show that the known pollen dispersal kernel explains g2. Second, we simulate a 1,000-generation pedigree using the known dispersal and spatial distribution and find that the resulting g2 is consistent with that observed from the field data. In contrast, a simulated population with uniform density underestimates g2, indicating that heterogeneous density promotes identity disequilibrium. Our study shows that heterogeneous density and leptokurtic dispersal can together explain the distribution of heterozygosity. AU - Surendranadh, Parvathy AU - Arathoon, Louise S AU - Baskett, Carina AU - Field, David AU - Pickup, Melinda AU - Barton, Nicholas H ID - 11411 IS - 3 JF - Genetics TI - Effects of fine-scale population structure on the distribution of heterozygosity in a long-term study of Antirrhinum majus VL - 221 ER - TY - JOUR AB - Selection accumulates information in the genome—it guides stochastically evolving populations toward states (genotype frequencies) that would be unlikely under neutrality. This can be quantified as the Kullback–Leibler (KL) divergence between the actual distribution of genotype frequencies and the corresponding neutral distribution. First, we show that this population-level information sets an upper bound on the information at the level of genotype and phenotype, limiting how precisely they can be specified by selection. Next, we study how the accumulation and maintenance of information is limited by the cost of selection, measured as the genetic load or the relative fitness variance, both of which we connect to the control-theoretic KL cost of control. The information accumulation rate is upper bounded by the population size times the cost of selection. This bound is very general, and applies across models (Wright–Fisher, Moran, diffusion) and to arbitrary forms of selection, mutation, and recombination. Finally, the cost of maintaining information depends on how it is encoded: Specifying a single allele out of two is expensive, but one bit encoded among many weakly specified loci (as in a polygenic trait) is cheap. AU - Hledik, Michal AU - Barton, Nicholas H AU - Tkačik, Gašper ID - 12081 IS - 36 JF - Proceedings of the National Academy of Sciences SN - 0027-8424 TI - Accumulation and maintenance of information in evolution VL - 119 ER - TY - THES AB - In evolve and resequence experiments, a population is sequenced, subjected to selection and then sequenced again, so that genetic changes before and after selection can be observed at the genetic level. Here, I use these studies to better understand the genetic basis of complex traits - traits which depend on more than a few genes. In the first chapter, I discuss the first evolve and resequence experiment, in which a population of mice, the so-called "Longshanks" mice, were selected for tibia length while their body mass was kept constant. The full pedigree is known. We observed a selection response on all chromosomes and used the infinitesimal model with linkage, a model which assumes an infinite number of genes with infinitesimally small effect sizes, as a null model. Results implied a very polygenic basis with a few loci of major effect standing out and changing in parallel. There was large variability between the different chromosomes in this study, probably due to LD. In chapter two, I go on to discuss the impact of LD, on the variability in an allele-frequency based summary statistic, giving an equation based on the initial allele frequencies, average pairwise LD, and the first four moments of the haplotype block copy number distribution. I describe this distribution by referring back to the founder generation. I then demonstrate how to infer selection via a maximum likelihood scheme on the example of a single locus and discuss how to extend this to more realistic scenarios. In chapter three, I discuss the second evolve and resequence experiment, in which a small population of Drosophila melanogaster was selected for increased pupal case size over 6 generations. The experiment was highly replicated with 27 lines selected within family and a known pedigree. We observed a phenotypic selection response of over one standard deviation. I describe the patterns in allele frequency data, including allele frequency changes and patterns of heterozygosity, and give ideas for future work. AU - Belohlavy, Stefanie ID - 11388 SN - 978-3-99078-018-3 TI - The genetic basis of complex traits studied via analysis of evolve and resequence experiments ER - TY - JOUR AB - Realistic models of biological processes typically involve interacting components on multiple scales, driven by changing environment and inherent stochasticity. Such models are often analytically and numerically intractable. We revisit a dynamic maximum entropy method that combines a static maximum entropy with a quasi-stationary approximation. This allows us to reduce stochastic non-equilibrium dynamics expressed by the Fokker-Planck equation to a simpler low-dimensional deterministic dynamics, without the need to track microscopic details. Although the method has been previously applied to a few (rather complicated) applications in population genetics, our main goal here is to explain and to better understand how the method works. We demonstrate the usefulness of the method for two widely studied stochastic problems, highlighting its accuracy in capturing important macroscopic quantities even in rapidly changing non-stationary conditions. For the Ornstein-Uhlenbeck process, the method recovers the exact dynamics whilst for a stochastic island model with migration from other habitats, the approximation retains high macroscopic accuracy under a wide range of scenarios in a dynamic environment. AU - Bod'ová, Katarína AU - Szep, Eniko AU - Barton, Nicholas H ID - 10535 IS - 12 JF - PLoS Computational Biology SN - 1553-734X TI - Dynamic maximum entropy provides accurate approximation of structured population dynamics VL - 17 ER - TY - JOUR AB - The Mytilus complex of marine mussel species forms a mosaic of hybrid zones, found across temperate regions of the globe. This allows us to study ‘replicated’ instances of secondary contact between closely related species. Previous work on this complex has shown that local introgression is both widespread and highly heterogeneous, and has identified SNPs that are outliers of differentiation between lineages. Here, we developed an ancestry‐informative panel of such SNPs. We then compared their frequencies in newly sampled populations, including samples from within the hybrid zones, and parental populations at different distances from the contact. Results show that close to the hybrid zones, some outlier loci are near to fixation for the heterospecific allele, suggesting enhanced local introgression, or the local sweep of a shared ancestral allele. Conversely, genomic cline analyses, treating local parental populations as the reference, reveal a globally high concordance among loci, albeit with a few signals of asymmetric introgression. Enhanced local introgression at specific loci is consistent with the early transfer of adaptive variants after contact, possibly including asymmetric bi‐stable variants (Dobzhansky‐Muller incompatibilities), or haplotypes loaded with fewer deleterious mutations. Having escaped one barrier, however, these variants can be trapped or delayed at the next barrier, confining the introgression locally. These results shed light on the decay of species barriers during phases of contact. AU - Simon, Alexis AU - Fraisse, Christelle AU - El Ayari, Tahani AU - Liautard‐Haag, Cathy AU - Strelkov, Petr AU - Welch, John J AU - Bierne, Nicolas ID - 8708 IS - 1 JF - Journal of Evolutionary Biology SN - 1010061X TI - How do species barriers decay? Concordance and local introgression in mosaic hybrid zones of mussels VL - 34 ER - TY - JOUR AB - Montane cloud forests are areas of high endemism, and are one of the more vulnerable terrestrial ecosystems to climate change. Thus, understanding how they both contribute to the generation of biodiversity, and will respond to ongoing climate change, are important and related challenges. The widely accepted model for montane cloud forest dynamics involves upslope forcing of their range limits with global climate warming. However, limited climate data provides some support for an alternative model, where range limits are forced downslope with climate warming. Testing between these two models is challenging, due to the inherent limitations of climate and pollen records. We overcome this with an alternative source of historical information, testing between competing model predictions using genomic data and demographic analyses for a species of beetle tightly associated to an oceanic island cloud forest. Results unequivocally support the alternative model: populations that were isolated at higher elevation peaks during the Last Glacial Maximum are now in contact and hybridizing at lower elevations. Our results suggest that genomic data are a rich source of information to further understand how montane cloud forest biodiversity originates, and how it is likely to be impacted by ongoing climate change. AU - Salces-Castellano, Antonia AU - Stankowski, Sean AU - Arribas, Paula AU - Patino, Jairo AU - Karger, Dirk N. AU - Butlin, Roger AU - Emerson, Brent C. ID - 8743 IS - 2 JF - Evolution SN - 0014-3820 TI - Long-term cloud forest response to climate warming revealed by insect speciation history VL - 75 ER - TY - JOUR AB - Domestication is a human‐induced selection process that imprints the genomes of domesticated populations over a short evolutionary time scale and that occurs in a given demographic context. Reconstructing historical gene flow, effective population size changes and their timing is therefore of fundamental interest to understand how plant demography and human selection jointly shape genomic divergence during domestication. Yet, the comparison under a single statistical framework of independent domestication histories across different crop species has been little evaluated so far. Thus, it is unclear whether domestication leads to convergent demographic changes that similarly affect crop genomes. To address this question, we used existing and new transcriptome data on three crop species of Solanaceae (eggplant, pepper and tomato), together with their close wild relatives. We fitted twelve demographic models of increasing complexity on the unfolded joint allele frequency spectrum for each wild/crop pair, and we found evidence for both shared and species‐specific demographic processes between species. A convergent history of domestication with gene flow was inferred for all three species, along with evidence of strong reduction in the effective population size during the cultivation stage of tomato and pepper. The absence of any reduction in size of the crop in eggplant stands out from the classical view of the domestication process; as does the existence of a “protracted period” of management before cultivation. Our results also suggest divergent management strategies of modern cultivars among species as their current demography substantially differs. Finally, the timing of domestication is species‐specific and supported by the few historical records available. AU - Arnoux, Stéphanie AU - Fraisse, Christelle AU - Sauvage, Christopher ID - 8928 IS - 2 JF - Journal of Evolutionary Biology SN - 1010061X TI - Genomic inference of complex domestication histories in three Solanaceae species VL - 34 ER - TY - JOUR AB - Marine environments are inhabited by a broad representation of the tree of life, yet our understanding of speciation in marine ecosystems is extremely limited compared with terrestrial and freshwater environments. Developing a more comprehensive picture of speciation in marine environments requires that we 'dive under the surface' by studying a wider range of taxa and ecosystems is necessary for a more comprehensive picture of speciation. Although studying marine evolutionary processes is often challenging, recent technological advances in different fields, from maritime engineering to genomics, are making it increasingly possible to study speciation of marine life forms across diverse ecosystems and taxa. Motivated by recent research in the field, including the 14 contributions in this issue, we highlight and discuss six axes of research that we think will deepen our understanding of speciation in the marine realm: (a) study a broader range of marine environments and organisms; (b) identify the reproductive barriers driving speciation between marine taxa; (c) understand the role of different genomic architectures underlying reproductive isolation; (d) infer the evolutionary history of divergence using model‐based approaches; (e) study patterns of hybridization and introgression between marine taxa; and (f) implement highly interdisciplinary, collaborative research programmes. In outlining these goals, we hope to inspire researchers to continue filling this critical knowledge gap surrounding the origins of marine biodiversity. AU - Faria, Rui AU - Johannesson, Kerstin AU - Stankowski, Sean ID - 9100 IS - 1 JF - Journal of Evolutionary Biology SN - 1010061X TI - Speciation in marine environments: Diving under the surface VL - 34 ER - TY - JOUR AB - Interspecific crossing experiments have shown that sex chromosomes play a major role in reproductive isolation between many pairs of species. However, their ability to act as reproductive barriers, which hamper interspecific genetic exchange, has rarely been evaluated quantitatively compared to Autosomes. This genome-wide limitation of gene flow is essential for understanding the complete separation of species, and thus speciation. Here, we develop a mainland-island model of secondary contact between hybridizing species of an XY (or ZW) sexual system. We obtain theoretical predictions for the frequency of introgressed alleles, and the strength of the barrier to neutral gene flow for the two types of chromosomes carrying multiple interspecific barrier loci. Theoretical predictions are obtained for scenarios where introgressed alleles are rare. We show that the same analytical expressions apply for sex chromosomes and autosomes, but with different sex-averaged effective parameters. The specific features of sex chromosomes (hemizygosity and absence of recombination in the heterogametic sex) lead to reduced levels of introgression on the X (or Z) compared to autosomes. This effect can be enhanced by certain types of sex-biased forces, but it remains overall small (except when alleles causing incompatibilities are recessive). We discuss these predictions in the light of empirical data comprising model-based tests of introgression and cline surveys in various biological systems. AU - Fraisse, Christelle AU - Sachdeva, Himani ID - 9168 IS - 2 JF - Genetics SN - 1943-2631 TI - The rates of introgression and barriers to genetic exchange between hybridizing species: Sex chromosomes vs autosomes VL - 217 ER -