TY - JOUR AB - Pollinators display a remarkable diversity of foraging strategies with flowering plants, from primarily mutualistic interactions to cheating through nectar robbery. Despite numerous studies on the effect of nectar robbing on components of plant fitness, its contribution to reproductive isolation is unclear. We experimentally tested the impact of different pollinator strategies in a natural hybrid zone between two subspecies of Antirrhinum majus with alternate flower colour guides. On either side of a steep cline in flower colour between Antirrhinum majus pseudomajus (magenta) and A. m. striatum (yellow), we quantified the behaviour of all floral visitors at different time points during the flowering season. Using long-run camera surveys, we quantify the impact of nectar robbing on the number of flowers visited per inflorescence and the flower probing time. We further experimentally tested the effect of nectar robbing on female reproductive success by manipulating the intensity of robbing. While robbing increased over time the number of legitimate visitors tended to decrease concomitantly. We found that the number of flowers pollinated on a focal inflorescence decreased with the number of prior robbing events. However, in the manipulative experiment, fruit set and fruit volume did not vary significantly between low robbing and control treatments. Our findings challenge the idea that robbers have a negative impact on plant fitness through female function. This study also adds to our understanding of the components of pollinator-mediated reproductive isolation and the maintenance of Antirrhinum hybrid zones. AU - Andalo, Christophe AU - Burrus, Monique AU - Paute, Sandrine AU - Lauzeral, Christine AU - Field, David ID - 5680 IS - 1 JF - Botany Letters SN - 23818107 TI - Prevalence of legitimate pollinators and nectar robbers and the consequences for fruit set in an Antirrhinum majus hybrid zone VL - 166 ER - TY - JOUR AB - The evolution of new species is made easier when traits under divergent ecological selection are also mating cues. Such ecological mating cues are now considered more common than previously thought, but we still know little about the genetic changes underlying their evolution or more generally about the genetic basis for assortative mating behaviors. Both tight physical linkage and the existence of large-effect preference loci will strengthen genetic associations between behavioral and ecological barriers, promoting the evolution of assortative mating. The warning patterns of Heliconius melpomene and H. cydno are under disruptive selection due to increased predation of nonmimetic hybrids and are used during mate recognition. We carried out a genome-wide quantitative trait locus (QTL) analysis of preference behaviors between these species and showed that divergent male preference has a simple genetic basis. We identify three QTLs that together explain a large proportion (approximately 60%) of the difference in preference behavior observed between the parental species. One of these QTLs is just 1.2 (0-4.8) centiMorgans (cM) from the major color pattern gene optix, and, individually, all three have a large effect on the preference phenotype. Genomic divergence between H. cydno and H. melpomene is high but broadly heterogenous, and admixture is reduced at the preference-optix color pattern locus but not the other preference QTLs. The simple genetic architecture we reveal will facilitate the evolution and maintenance of new species despite ongoing gene flow by coupling behavioral and ecological aspects of reproductive isolation. AU - Merrill, Richard M. AU - Rastas, Pasi AU - Martin, Simon H. AU - Melo Hurtado, Maria C AU - Barker, Sarah AU - Davey, John AU - Mcmillan, W. Owen AU - Jiggins, Chris D. ID - 6022 IS - 2 JF - PLoS Biology TI - Genetic dissection of assortative mating behavior VL - 17 ER - TY - GEN AU - Merrill, Richard M. AU - Rastas, Pasi AU - Martin, Simon H. AU - Melo Hurtado, Maria C AU - Barker, Sarah AU - Davey, John AU - Mcmillan, W. Owen AU - Jiggins, Chris D. ID - 9801 TI - Raw behavioral data ER - TY - JOUR AB - Both classical and recent studies suggest that chromosomal inversion polymorphisms are important in adaptation and speciation. However, biases in discovery and reporting of inversions make it difficult to assess their prevalence and biological importance. Here, we use an approach based on linkage disequilibrium among markers genotyped for samples collected across a transect between contrasting habitats to detect chromosomal rearrangements de novo. We report 17 polymorphic rearrangements in a single locality for the coastal marine snail, Littorina saxatilis. Patterns of diversity in the field and of recombination in controlled crosses provide strong evidence that at least the majority of these rearrangements are inversions. Most show clinal changes in frequency between habitats, suggestive of divergent selection, but only one appears to be fixed for different arrangements in the two habitats. Consistent with widespread evidence for balancing selection on inversion polymorphisms, we argue that a combination of heterosis and divergent selection can explain the observed patterns and should be considered in other systems spanning environmental gradients. AU - Faria, Rui AU - Chaube, Pragya AU - Morales, Hernán E. AU - Larsson, Tomas AU - Lemmon, Alan R. AU - Lemmon, Emily M. AU - Rafajlović, Marina AU - Panova, Marina AU - Ravinet, Mark AU - Johannesson, Kerstin AU - Westram, Anja M AU - Butlin, Roger K. ID - 6095 IS - 6 JF - Molecular Ecology SN - 0962-1083 TI - Multiple chromosomal rearrangements in a hybrid zone between Littorina saxatilis ecotypes VL - 28 ER - TY - JOUR AB - Great care is needed when interpreting claims about the genetic basis of human variation based on data from genome-wide association studies. AU - Barton, Nicholas H AU - Hermisson, Joachim AU - Nordborg, Magnus ID - 6230 JF - eLife TI - Why structure matters VL - 8 ER - TY - JOUR AB - One of the most striking and consistent results in speciation genomics is the heterogeneous divergence observed across the genomes of closely related species. This pattern was initially attributed to different levels of gene exchange—with divergence preserved at loci generating a barrier to gene flow but homogenized at unlinked neutral loci. Although there is evidence to support this model, it is now recognized that interpreting patterns of divergence across genomes is not so straightforward. One problem is that heterogenous divergence between populations can also be generated by other processes (e.g. recurrent selective sweeps or background selection) without any involvement of differential gene flow. Thus, integrated studies that identify which loci are likely subject to divergent selection are required to shed light on the interplay between selection and gene flow during the early phases of speciation. In this issue of Molecular Ecology, Rifkin et al. (2019) confront this challenge using a pair of sister morning glory species. They wisely design their sampling to take the geographic context of individuals into account, including geographically isolated (allopatric) and co‐occurring (sympatric) populations. This enabled them to show that individuals are phenotypically less differentiated in sympatry. They also found that the loci that resist introgression are enriched for those most differentiated in allopatry and loci that exhibit signals of divergent selection. One great strength of the study is the combination of methods from population genetics and molecular evolution, including the development of a model to simultaneously infer admixture proportions and selfing rates. AU - Field, David AU - Fraisse, Christelle ID - 6466 IS - 7 JF - Molecular ecology TI - Breaking down barriers in morning glories VL - 28 ER - TY - JOUR AB - Fitness interactions between mutations can influence a population’s evolution in many different ways. While epistatic effects are difficult to measure precisely, important information is captured by the mean and variance of log fitnesses for individuals carrying different numbers of mutations. We derive predictions for these quantities from a class of simple fitness landscapes, based on models of optimizing selection on quantitative traits. We also explore extensions to the models, including modular pleiotropy, variable effect sizes, mutational bias and maladaptation of the wild type. We illustrate our approach by reanalysing a large dataset of mutant effects in a yeast snoRNA (small nucleolar RNA). Though characterized by some large epistatic effects, these data give a good overall fit to the non-epistatic null model, suggesting that epistasis might have limited influence on the evolutionary dynamics in this system. We also show how the amount of epistasis depends on both the underlying fitness landscape and the distribution of mutations, and so is expected to vary in consistent ways between new mutations, standing variation and fixed mutations. AU - Fraisse, Christelle AU - Welch, John J. ID - 6467 IS - 4 JF - Biology Letters SN - 17449561 TI - The distribution of epistasis on simple fitness landscapes VL - 15 ER - TY - JOUR AB - The environment changes constantly at various time scales and, in order to survive, species need to keep adapting. Whether these species succeed in avoiding extinction is a major evolutionary question. Using a multilocus evolutionary model of a mutation‐limited population adapting under strong selection, we investigate the effects of the frequency of environmental fluctuations on adaptation. Our results rely on an “adaptive‐walk” approximation and use mathematical methods from evolutionary computation theory to investigate the interplay between fluctuation frequency, the similarity of environments, and the number of loci contributing to adaptation. First, we assume a linear additive fitness function, but later generalize our results to include several types of epistasis. We show that frequent environmental changes prevent populations from reaching a fitness peak, but they may also prevent the large fitness loss that occurs after a single environmental change. Thus, the population can survive, although not thrive, in a wide range of conditions. Furthermore, we show that in a frequently changing environment, the similarity of threats that a population faces affects the level of adaptation that it is able to achieve. We check and supplement our analytical results with simulations. AU - Trubenova, Barbora AU - Krejca, Martin AU - Lehre, Per Kristian AU - Kötzing, Timo ID - 6637 IS - 7 JF - Evolution TI - Surfing on the seascape: Adaptation in a changing environment VL - 73 ER - TY - JOUR AB - This paper analyzes how partial selfing in a large source population influences its ability to colonize a new habitat via the introduction of a few founder individuals. Founders experience inbreeding depression due to partially recessive deleterious alleles as well as maladaptation to the new environment due to selection on a large number of additive loci. I first introduce a simplified version of the Inbreeding History Model (Kelly, 2007) in order to characterize mutation‐selection balance in a large, partially selfing source population under selection involving multiple non‐identical loci. I then use individual‐based simulations to study the eco‐evolutionary dynamics of founders establishing in the new habitat under a model of hard selection. The study explores how selfing rate shapes establishment probabilities of founders via effects on both inbreeding depression and adaptability to the new environment, and also distinguishes the effects of selfing on the initial fitness of founders from its effects on the long‐term adaptive response of the populations they found. A high rate of (but not complete) selfing is found to aid establishment over a wide range of parameters, even in the absence of mate limitation. The sensitivity of the results to assumptions about the nature of polygenic selection are discussed. AU - Sachdeva, Himani ID - 6680 IS - 9 JF - Evolution SN - 0014-3820 TI - Effect of partial selfing and polygenic selection on establishment in a new habitat VL - 73 ER - TY - GEN AB - Evolutionary studies are often limited by missing data that are critical to understanding the history of selection. Selection experiments, which reproduce rapid evolution under controlled conditions, are excellent tools to study how genomes evolve under selection. Here we present a genomic dissection of the Longshanks selection experiment, in which mice were selectively bred over 20 generations for longer tibiae relative to body mass, resulting in 13% longer tibiae in two replicates. We synthesized evolutionary theory, genome sequences and molecular genetics to understand the selection response and found that it involved both polygenic adaptation and discrete loci of major effect, with the strongest loci tending to be selected in parallel between replicates. We show that selection may favor de-repression of bone growth through inactivating two limb enhancers of an inhibitor, Nkx3-2. Our integrative genomic analyses thus show that it is possible to connect individual base-pair changes to the overall selection response. AU - Castro, João Pl AU - Yancoskie, Michelle N. AU - Marchini, Marta AU - Belohlavy, Stefanie AU - Hiramatsu, Layla AU - Kučka, Marek AU - Beluch, William H. AU - Naumann, Ronald AU - Skuplik, Isabella AU - Cobb, John AU - Barton, Nicholas H AU - Rolian, Campbell AU - Chan, Yingguang Frank ID - 9804 TI - Data from: An integrative genomic analysis of the Longshanks selection experiment for longer limbs in mice ER - TY - GEN AB - This paper analyzes how partial selfing in a large source population influences its ability to colonize a new habitat via the introduction of a few founder individuals. Founders experience inbreeding depression due to partially recessive deleterious alleles as well as maladaptation to the new environment due to selection on a large number of additive loci. I first introduce a simplified version of the Inbreeding History Model (Kelly, 2007) in order to characterize mutation-selection balance in a large, partially selfing source population under selection involving multiple non-identical loci. I then use individual-based simulations to study the eco-evolutionary dynamics of founders establishing in the new habitat under a model of hard selection. The study explores how selfing rate shapes establishment probabilities of founders via effects on both inbreeding depression and adaptability to the new environment, and also distinguishes the effects of selfing on the initial fitness of founders from its effects on the long-term adaptive response of the populations they found. A high rate of (but not complete) selfing is found to aid establishment over a wide range of parameters, even in the absence of mate limitation. The sensitivity of the results to assumptions about the nature of polygenic selection are discussed. AU - Sachdeva, Himani ID - 9802 TI - Data from: Effect of partial selfing and polygenic selection on establishment in a new habitat ER - TY - JOUR AB - The green‐beard effect is one proposed mechanism predicted to underpin the evolu‐tion of altruistic behavior. It relies on the recognition and the selective help of altruists to each other in order to promote and sustain altruistic behavior. However, this mechanism has often been dismissed as unlikely or uncommon, as it is assumed that both the signaling trait and altruistic trait need to be encoded by the same gene or through tightly linked genes. Here, we use models of indirect genetic effects (IGEs) to find the minimum correlation between the signaling and altruistic trait required for the evolution of the latter. We show that this correlation threshold depends on the strength of the interaction (influence of the green beard on the expression of the altruistic trait), as well as the costs and benefits of the altruistic behavior. We further show that this correlation does not necessarily have to be high and support our analytical results by simulations. AU - Trubenova, Barbora AU - Hager, Reinmar ID - 6795 IS - 17 JF - Ecology and Evolution TI - Green beards in the light of indirect genetic effects VL - 9 ER - TY - JOUR AB - * Understanding the mechanisms causing phenotypic differences between females and males has long fascinated evolutionary biologists. An extensive literature exists on animal sexual dimorphism but less information is known about sex differences in plants, particularly the extent of geographical variation in sexual dimorphism and its life‐cycle dynamics. * Here, we investigated patterns of genetically based sexual dimorphism in vegetative and reproductive traits of a wind‐pollinated dioecious plant, Rumex hastatulus, across three life‐cycle stages using open‐pollinated families from 30 populations spanning the geographic range and chromosomal variation (XY and XY1Y2) of the species. * The direction and degree of sexual dimorphism was highly variable among populations and life‐cycle stages. Sex‐specific differences in reproductive function explained a significant amount of temporal change in sexual dimorphism. For several traits, geographical variation in sexual dimorphism was associated with bioclimatic parameters, likely due to the differential responses of the sexes to climate. We found no systematic differences in sexual dimorphism between chromosome races. * Sex‐specific trait differences in dioecious plants largely result from a balance between sexual and natural selection on resource allocation. Our results indicate that abiotic factors associated with geographical context also play a role in modifying sexual dimorphism during the plant life‐cycle. AU - Puixeu Sala, Gemma AU - Pickup, Melinda AU - Field, David AU - Barrett, Spencer C.H. ID - 6831 IS - 3 JF - New Phytologist TI - Variation in sexual dimorphism in a wind-pollinated plant: The influence of geographical context and life-cycle dynamics VL - 224 ER - TY - GEN AB - Understanding the mechanisms causing phenotypic differences between females and males has long fascinated evolutionary biologists. An extensive literature exists on animal sexual dimorphism but less is known about sex differences in plants, particularly the extent of geographical variation in sexual dimorphism and its life-cycle dynamics. Here, we investigate patterns of genetically-based sexual dimorphism in vegetative and reproductive traits of a wind-pollinated dioecious plant, Rumex hastatulus, across three life-cycle stages using open-pollinated families from 30 populations spanning the geographic range and chromosomal variation (XY and XY1Y2) of the species. The direction and degree of sexual dimorphism was highly variable among populations and life-cycle stages. Sex-specific differences in reproductive function explained a significant amount of temporal change in sexual dimorphism. For several traits, geographical variation in sexual dimorphism was associated with bioclimatic parameters, likely due to the differential responses of the sexes to climate. We found no systematic differences in sexual dimorphism between chromosome races. Sex-specific trait differences in dioecious plants largely result from a balance between sexual and natural selection on resource allocation. Our results indicate that abiotic factors associated with geographical context also play a role in modifying sexual dimorphism during the plant life cycle. AU - Puixeu Sala, Gemma AU - Pickup, Melinda AU - Field, David AU - Barrett, Spencer C.H. ID - 9803 TI - Data from: Variation in sexual dimorphism in a wind-pollinated plant: the influence of geographical context and life-cycle dynamics ER - TY - JOUR AB - Many traits of interest are highly heritable and genetically complex, meaning that much of the variation they exhibit arises from differences at numerous loci in the genome. Complex traits and their evolution have been studied for more than a century, but only in the last decade have genome-wide association studies (GWASs) in humans begun to reveal their genetic basis. Here, we bring these threads of research together to ask how findings from GWASs can further our understanding of the processes that give rise to heritable variation in complex traits and of the genetic basis of complex trait evolution in response to changing selection pressures (i.e., of polygenic adaptation). Conversely, we ask how evolutionary thinking helps us to interpret findings from GWASs and informs related efforts of practical importance. AU - Sella, Guy AU - Barton, Nicholas H ID - 6855 JF - Annual Review of Genomics and Human Genetics SN - 1527-8204 TI - Thinking about the evolution of complex traits in the era of genome-wide association studies VL - 20 ER - TY - JOUR AU - Barton, Nicholas H ID - 6858 IS - 2 JF - National Science Review SN - 2095-5138 TI - Is speciation driven by cycles of mixing and isolation? VL - 6 ER - TY - JOUR AB - Gene Drives are regarded as future tools with a high potential for population control. Due to their inherent ability to overcome the rules of Mendelian inheritance, gene drives (GD) may spread genes rapidly through populations of sexually reproducing organisms. A release of organisms carrying a GD would constitute a paradigm shift in the handling of genetically modified organisms because gene drive organisms (GDO) are designed to drive their transgenes into wild populations and thereby increase the number of GDOs. The rapid development in this field and its focus on wild populations demand a prospective risk assessment with a focus on exposure related aspects. Presently, it is unclear how adequate risk management could be guaranteed to limit the spread of GDs in time and space, in order to avoid potential adverse effects in socio‐ecological systems. The recent workshop on the “Evaluation of Spatial and Temporal Control of Gene Drives” hosted by the Institute of Safety/Security and Risk Sciences (ISR) in Vienna aimed at gaining some insight into the potential population dynamic behavior of GDs and appropriate measures of control. Scientists from France, Germany, England, and the USA discussed both topics in this meeting on April 4–5, 2019. This article summarizes results of the workshop. AU - Giese, B AU - Friess, J L AU - Schetelig, M F AU - Barton, Nicholas H AU - Messer, Philip AU - Debarre, Florence AU - Meimberg, H AU - Windbichler, N AU - Boete, C ID - 6857 IS - 11 JF - BioEssays TI - Gene Drives: Dynamics and regulatory matters – A report from the workshop “Evaluation of spatial and temporal control of Gene Drives”, 4 – 5 April 2019, Vienna VL - 41 ER - TY - GEN AB - Genetic incompatibilities contribute to reproductive isolation between many diverging populations, but it is still unclear to what extent they play a role if divergence happens with gene flow. In contact zones between the "Crab" and "Wave" ecotypes of the snail Littorina saxatilis divergent selection forms strong barriers to gene flow, while the role of postzygotic barriers due to selection against hybrids remains unclear. High embryo abortion rates in this species could indicate the presence of such barriers. Postzygotic barriers might include genetic incompatibilities (e.g. Dobzhansky-Muller incompatibilities) but also maladaptation, both expected to be most pronounced in contact zones. In addition, embryo abortion might reflect physiological stress on females and embryos independent of any genetic stress. We examined all embryos of >500 females sampled outside and inside contact zones of three populations in Sweden. Females' clutch size ranged from 0 to 1011 embryos (mean 130±123) and abortion rates varied between 0 and100% (mean 12%). We described female genotypes by using a hybrid index based on hundreds of SNPs differentiated between ecotypes with which we characterised female genotypes. We also calculated female SNP heterozygosity and inversion karyotype. Clutch size did not vary with female hybrid index and abortion rates were only weakly related to hybrid index in two sites but not at all in a third site. No additional variation in abortion rate was explained by female SNP heterozygosity, but increased female inversion heterozygosity added slightly to increased abortion. Our results show only weak and probably biologically insignificant postzygotic barriers contributing to ecotype divergence and the high and variable abortion rates were marginally, if at all, explained by hybrid index of females. AU - Johannesson, Kerstin AU - Zagrodzka, Zuzanna AU - Faria, Rui AU - Westram, Anja M AU - Butlin, Roger ID - 13067 TI - Data from: Is embryo abortion a postzygotic barrier to gene flow between Littorina ecotypes? ER - TY - JOUR AB - The study of parallel ecological divergence provides important clues to the operation of natural selection. Parallel divergence often occurs in heterogeneous environments with different kinds of environmental gradients in different locations, but the genomic basis underlying this process is unknown. We investigated the genomics of rapid parallel adaptation in the marine snail Littorina saxatilis in response to two independent environmental axes (crab-predation versus wave-action and low-shore versus high-shore). Using pooled whole-genome resequencing, we show that sharing of genomic regions of high differentiation between environments is generally low but increases at smaller spatial scales. We identify different shared genomic regions of divergence for each environmental axis and show that most of these regions overlap with candidate chromosomal inversions. Several inversion regions are divergent and polymorphic across many localities. We argue that chromosomal inversions could store shared variation that fuels rapid parallel adaptation to heterogeneous environments, possibly as balanced polymorphism shared by adaptive gene flow. AU - Morales, Hernán E. AU - Faria, Rui AU - Johannesson, Kerstin AU - Larsson, Tomas AU - Panova, Marina AU - Westram, Anja M AU - Butlin, Roger K. ID - 7393 IS - 12 JF - Science Advances SN - 2375-2548 TI - Genomic architecture of parallel ecological divergence: Beyond a single environmental contrast VL - 5 ER - TY - CHAP AB - We review the history of population genetics, starting with its origins a century ago from the synthesis between Mendel and Darwin's ideas, through to the recent development of sophisticated schemes of inference from sequence data, based on the coalescent. We explain the close relation between the coalescent and a diffusion process, which we illustrate by their application to understand spatial structure. We summarise the powerful methods available for analysis of multiple loci, when linkage equilibrium can be assumed, and then discuss approaches to the more challenging case, where associations between alleles require that we follow genotype, rather than allele, frequencies. Though we can hardly cover the whole of population genetics, we give an overview of the current state of the subject, and future challenges to it. AU - Barton, Nicholas H AU - Etheridge, Alison ED - Balding, David ED - Moltke, Ida ED - Marioni, John ID - 8281 SN - 9781119429142 T2 - Handbook of statistical genomics TI - Mathematical models in population genetics ER - TY - GEN AB - The spread of adaptive alleles is fundamental to evolution, and in theory, this process is well‐understood. However, only rarely can we follow this process—whether it originates from the spread of a new mutation, or by introgression from another population. In this issue of Molecular Ecology, Hanemaaijer et al. (2018) report on a 25‐year long study of the mosquitoes Anopheles gambiae (Figure 1) and Anopheles coluzzi in Mali, based on genotypes at 15 single‐nucleotide polymorphism (SNP). The species are usually reproductively isolated from each other, but in 2002 and 2006, bursts of hybridization were observed, when F1 hybrids became abundant. Alleles backcrossed from A. gambiae into A. coluzzi, but after the first event, these declined over the following years. In contrast, after 2006, an insecticide resistance allele that had established in A. gambiae spread into A. coluzzi, and rose to high frequency there, over 6 years (~75 generations). Whole genome sequences of 74 individuals showed that A. gambiae SNP from across the genome had become common in the A. coluzzi population, but that most of these were clustered in 34 genes around the resistance locus. A new set of SNP from 25 of these genes were assayed over time; over the 4 years since near‐fixation of the resistance allele; some remained common, whereas others declined. What do these patterns tell us about this introgression event? AU - Barton, Nicholas H ID - 9805 TI - Data from: The consequences of an introgression event ER - TY - THES AB - Transcription factors, by binding to specific sequences on the DNA, control the precise spatio-temporal expression of genes inside a cell. However, this specificity is limited, leading to frequent incorrect binding of transcription factors that might have deleterious consequences on the cell. By constructing a biophysical model of TF-DNA binding in the context of gene regulation, I will first explore how regulatory constraints can strongly shape the distribution of a population in sequence space. Then, by directly linking this to a picture of multiple types of transcription factors performing their functions simultaneously inside the cell, I will explore the extent of regulatory crosstalk -- incorrect binding interactions between transcription factors and binding sites that lead to erroneous regulatory states -- and understand the constraints this places on the design of regulatory systems. I will then develop a generic theoretical framework to investigate the coevolution of multiple transcription factors and multiple binding sites, in the context of a gene regulatory network that performs a certain function. As a particular tractable version of this problem, I will consider the evolution of two transcription factors when they transmit upstream signals to downstream target genes. Specifically, I will describe the evolutionary steady states and the evolutionary pathways involved, along with their timescales, of a system that initially undergoes a transcription factor duplication event. To connect this important theoretical model to the prominent biological event of transcription factor duplication giving rise to paralogous families, I will then describe a bioinformatics analysis of C2H2 Zn-finger transcription factors, a major family in humans, and focus on the patterns of evolution that paralogs have undergone in their various protein domains in the recent past. AU - Prizak, Roshan ID - 6071 SN - 2663-337X TI - Coevolution of transcription factors and their binding sites in sequence space ER - TY - JOUR AB - Plant mating systems play a key role in structuring genetic variation both within and between species. In hybrid zones, the outcomes and dynamics of hybridization are usually interpreted as the balance between gene flow and selection against hybrids. Yet, mating systems can introduce selective forces that alter these expectations; with diverse outcomes for the level and direction of gene flow depending on variation in outcrossing and whether the mating systems of the species pair are the same or divergent. We present a survey of hybridization in 133 species pairs from 41 plant families and examine how patterns of hybridization vary with mating system. We examine if hybrid zone mode, level of gene flow, asymmetries in gene flow and the frequency of reproductive isolating barriers vary in relation to mating system/s of the species pair. We combine these results with a simulation model and examples from the literature to address two general themes: (i) the two‐way interaction between introgression and the evolution of reproductive systems, and (ii) how mating system can facilitate or restrict interspecific gene flow. We conclude that examining mating system with hybridization provides unique opportunities to understand divergence and the processes underlying reproductive isolation. AU - Pickup, Melinda AU - Barton, Nicholas H AU - Brandvain, Yaniv AU - Fraisse, Christelle AU - Yakimowski, Sarah AU - Dixit, Tanmay AU - Lexer, Christian AU - Cereghetti, Eva AU - Field, David ID - 6856 IS - 3 JF - New Phytologist SN - 0028-646X TI - Mating system variation in hybrid zones: Facilitation, barriers and asymmetries to gene flow VL - 224 ER - TY - JOUR AB - Pleiotropy is the well-established idea that a single mutation affects multiple phenotypes. If a mutation has opposite effects on fitness when expressed in different contexts, then genetic conflict arises. Pleiotropic conflict is expected to reduce the efficacy of selection by limiting the fixation of beneficial mutations through adaptation, and the removal of deleterious mutations through purifying selection. Although this has been widely discussed, in particular in the context of a putative “gender load,” it has yet to be systematically quantified. In this work, we empirically estimate to which extent different pleiotropic regimes impede the efficacy of selection in Drosophila melanogaster. We use whole-genome polymorphism data from a single African population and divergence data from D. simulans to estimate the fraction of adaptive fixations (α), the rate of adaptation (ωA), and the direction of selection (DoS). After controlling for confounding covariates, we find that the different pleiotropic regimes have a relatively small, but significant, effect on selection efficacy. Specifically, our results suggest that pleiotropic sexual antagonism may restrict the efficacy of selection, but that this conflict can be resolved by limiting the expression of genes to the sex where they are beneficial. Intermediate levels of pleiotropy across tissues and life stages can also lead to maladaptation in D. melanogaster, due to inefficient purifying selection combined with low frequency of mutations that confer a selective advantage. Thus, our study highlights the need to consider the efficacy of selection in the context of antagonistic pleiotropy, and of genetic conflict in general. AU - Fraisse, Christelle AU - Puixeu Sala, Gemma AU - Vicoso, Beatriz ID - 6089 IS - 3 JF - Molecular biology and evolution SN - 0737-4038 TI - Pleiotropy modulates the efficacy of selection in drosophila melanogaster VL - 36 ER - TY - JOUR AB - Cells need to reliably sense external ligand concentrations to achieve various biological functions such as chemotaxis or signaling. The molecular recognition of ligands by surface receptors is degenerate in many systems, leading to crosstalk between ligand-receptor pairs. Crosstalk is often thought of as a deviation from optimal specific recognition, as the binding of noncognate ligands can interfere with the detection of the receptor's cognate ligand, possibly leading to a false triggering of a downstream signaling pathway. Here we quantify the optimal precision of sensing the concentrations of multiple ligands by a collection of promiscuous receptors. We demonstrate that crosstalk can improve precision in concentration sensing and discrimination tasks. To achieve superior precision, the additional information about ligand concentrations contained in short binding events of the noncognate ligand should be exploited. We present a proofreading scheme to realize an approximate estimation of multiple ligand concentrations that reaches a precision close to the derived optimal bounds. Our results help rationalize the observed ubiquity of receptor crosstalk in molecular sensing. AU - Carballo-Pacheco, Martín AU - Desponds, Jonathan AU - Gavrilchenko, Tatyana AU - Mayer, Andreas AU - Prizak, Roshan AU - Reddy, Gautam AU - Nemenman, Ilya AU - Mora, Thierry ID - 6090 IS - 2 JF - Physical Review E TI - Receptor crosstalk improves concentration sensing of multiple ligands VL - 99 ER - TY - JOUR AB - Evolutionary studies are often limited by missing data that are critical to understanding the history of selection. Selection experiments, which reproduce rapid evolution under controlled conditions, are excellent tools to study how genomes evolve under selection. Here we present a genomic dissection of the Longshanks selection experiment, in which mice were selectively bred over 20 generations for longer tibiae relative to body mass, resulting in 13% longer tibiae in two replicates. We synthesized evolutionary theory, genome sequences and molecular genetics to understand the selection response and found that it involved both polygenic adaptation and discrete loci of major effect, with the strongest loci tending to be selected in parallel between replicates. We show that selection may favor de-repression of bone growth through inactivating two limb enhancers of an inhibitor, Nkx3-2. Our integrative genomic analyses thus show that it is possible to connect individual base-pair changes to the overall selection response. AU - Castro, João Pl AU - Yancoskie, Michelle N. AU - Marchini, Marta AU - Belohlavy, Stefanie AU - Hiramatsu, Layla AU - Kučka, Marek AU - Beluch, William H. AU - Naumann, Ronald AU - Skuplik, Isabella AU - Cobb, John AU - Barton, Nicholas H AU - Rolian, Campbell AU - Chan, Yingguang Frank ID - 6713 JF - eLife TI - An integrative genomic analysis of the Longshanks selection experiment for longer limbs in mice VL - 8 ER - TY - JOUR AB - More than 100 years after Grigg’s influential analysis of species’ borders, the causes of limits to species’ ranges still represent a puzzle that has never been understood with clarity. The topic has become especially important recently as many scientists have become interested in the potential for species’ ranges to shift in response to climate change—and yet nearly all of those studies fail to recognise or incorporate evolutionary genetics in a way that relates to theoretical developments. I show that range margins can be understood based on just two measurable parameters: (i) the fitness cost of dispersal—a measure of environmental heterogeneity—and (ii) the strength of genetic drift, which reduces genetic diversity. Together, these two parameters define an ‘expansion threshold’: adaptation fails when genetic drift reduces genetic diversity below that required for adaptation to a heterogeneous environment. When the key parameters drop below this expansion threshold locally, a sharp range margin forms. When they drop below this threshold throughout the species’ range, adaptation collapses everywhere, resulting in either extinction or formation of a fragmented metapopulation. Because the effects of dispersal differ fundamentally with dimension, the second parameter—the strength of genetic drift—is qualitatively different compared to a linear habitat. In two-dimensional habitats, genetic drift becomes effectively independent of selection. It decreases with ‘neighbourhood size’—the number of individuals accessible by dispersal within one generation. Moreover, in contrast to earlier predictions, which neglected evolution of genetic variance and/or stochasticity in two dimensions, dispersal into small marginal populations aids adaptation. This is because the reduction of both genetic and demographic stochasticity has a stronger effect than the cost of dispersal through increased maladaptation. The expansion threshold thus provides a novel, theoretically justified, and testable prediction for formation of the range margin and collapse of the species’ range. AU - Polechova, Jitka ID - 315 IS - 6 JF - PLoS Biology SN - 15449173 TI - Is the sky the limit? On the expansion threshold of a species’ range VL - 16 ER - TY - GEN AB - Both classical and recent studies suggest that chromosomal inversion polymorphisms are important in adaptation and speciation. However, biases in discovery and reporting of inversions make it difficult to assess their prevalence and biological importance. Here, we use an approach based on linkage disequilibrium among markers genotyped for samples collected across a transect between contrasting habitats to detect chromosomal rearrangements de novo. We report 17 polymorphic rearrangements in a single locality for the coastal marine snail, Littorina saxatilis. Patterns of diversity in the field and of recombination in controlled crosses provide strong evidence that at least the majority of these rearrangements are inversions. Most show clinal changes in frequency between habitats, suggestive of divergent selection, but only one appears to be fixed for different arrangements in the two habitats. Consistent with widespread evidence for balancing selection on inversion polymorphisms, we argue that a combination of heterosis and divergent selection can explain the observed patterns and should be considered in other systems spanning environmental gradients. AU - Faria, Rui AU - Chaube, Pragya AU - Morales, Hernán E. AU - Larsson, Tomas AU - Lemmon, Alan R. AU - Lemmon, Emily M. AU - Rafajlović, Marina AU - Panova, Marina AU - Ravinet, Mark AU - Johannesson, Kerstin AU - Westram, Anja M AU - Butlin, Roger K. ID - 9837 TI - Data from: Multiple chromosomal rearrangements in a hybrid zone between Littorina saxatilis ecotypes ER - TY - JOUR AB - Herd immunity, a process in which resistant individuals limit the spread of a pathogen among susceptible hosts has been extensively studied in eukaryotes. Even though bacteria have evolved multiple immune systems against their phage pathogens, herd immunity in bacteria remains unexplored. Here we experimentally demonstrate that herd immunity arises during phage epidemics in structured and unstructured Escherichia coli populations consisting of differing frequencies of susceptible and resistant cells harboring CRISPR immunity. In addition, we develop a mathematical model that quantifies how herd immunity is affected by spatial population structure, bacterial growth rate, and phage replication rate. Using our model we infer a general epidemiological rule describing the relative speed of an epidemic in partially resistant spatially structured populations. Our experimental and theoretical findings indicate that herd immunity may be important in bacterial communities, allowing for stable coexistence of bacteria and their phages and the maintenance of polymorphism in bacterial immunity. AU - Payne, Pavel AU - Geyrhofer, Lukas AU - Barton, Nicholas H AU - Bollback, Jonathan P ID - 423 JF - eLife TI - CRISPR-based herd immunity can limit phage epidemics in bacterial populations VL - 7 ER - TY - GEN AB - Herd immunity, a process in which resistant individuals limit the spread of a pathogen among susceptible hosts has been extensively studied in eukaryotes. Even though bacteria have evolved multiple immune systems against their phage pathogens, herd immunity in bacteria remains unexplored. Here we experimentally demonstrate that herd immunity arises during phage epidemics in structured and unstructured Escherichia coli populations consisting of differing frequencies of susceptible and resistant cells harboring CRISPR immunity. In addition, we develop a mathematical model that quantifies how herd immunity is affected by spatial population structure, bacterial growth rate, and phage replication rate. Using our model we infer a general epidemiological rule describing the relative speed of an epidemic in partially resistant spatially structured populations. Our experimental and theoretical findings indicate that herd immunity may be important in bacterial communities, allowing for stable coexistence of bacteria and their phages and the maintenance of polymorphism in bacterial immunity. AU - Payne, Pavel AU - Geyrhofer, Lukas AU - Barton, Nicholas H AU - Bollback, Jonathan P ID - 9840 TI - Data from: CRISPR-based herd immunity limits phage epidemics in bacterial populations ER - TY - JOUR AB - Maladapted individuals can only colonise a new habitat if they can evolve a positive growth rate fast enough to avoid extinction, a process known as evolutionary rescue. We treat log fitness at low density in the new habitat as a single polygenic trait and thus use the infinitesimal model to follow the evolution of the growth rate; this assumes that the trait values of offspring of a sexual union are normally distributed around the mean of the parents’ trait values, with variance that depends only on the parents’ relatedness. The probability that a single migrant can establish depends on just two parameters: the mean and genetic variance of the trait in the source population. The chance of success becomes small if migrants come from a population with mean growth rate in the new habitat more than a few standard deviations below zero; this chance depends roughly equally on the probability that the initial founder is unusually fit, and on the subsequent increase in growth rate of its offspring as a result of selection. The loss of genetic variation during the founding event is substantial, but highly variable. With continued migration at rate M, establishment is inevitable; when migration is rare, the expected time to establishment decreases inversely with M. However, above a threshold migration rate, the population may be trapped in a ‘sink’ state, in which adaptation is held back by gene flow; above this threshold, the expected time to establishment increases exponentially with M. This threshold behaviour is captured by a deterministic approximation, which assumes a Gaussian distribution of the trait in the founder population with mean and variance evolving deterministically. By assuming a constant genetic variance, we also develop a diffusion approximation for the joint distribution of population size and trait mean, which extends to include stabilising selection and density regulation. Divergence of the population from its ancestors causes partial reproductive isolation, which we measure through the reproductive value of migrants into the newly established population. AU - Barton, Nicholas H AU - Etheridge, Alison ID - 564 IS - 7 JF - Theoretical Population Biology TI - Establishment in a new habitat by polygenic adaptation VL - 122 ER - TY - JOUR AB - In continuous populations with local migration, nearby pairs of individuals have on average more similar genotypes than geographically well separated pairs. A barrier to gene flow distorts this classical pattern of isolation by distance. Genetic similarity is decreased for sample pairs on different sides of the barrier and increased for pairs on the same side near the barrier. Here, we introduce an inference scheme that utilizes this signal to detect and estimate the strength of a linear barrier to gene flow in two-dimensions. We use a diffusion approximation to model the effects of a barrier on the geographical spread of ancestry backwards in time. This approach allows us to calculate the chance of recent coalescence and probability of identity by descent. We introduce an inference scheme that fits these theoretical results to the geographical covariance structure of bialleleic genetic markers. It can estimate the strength of the barrier as well as several demographic parameters. We investigate the power of our inference scheme to detect barriers by applying it to a wide range of simulated data. We also showcase an example application to a Antirrhinum majus (snapdragon) flower color hybrid zone, where we do not detect any signal of a strong genome wide barrier to gene flow. AU - Ringbauer, Harald AU - Kolesnikov, Alexander AU - Field, David AU - Barton, Nicholas H ID - 563 IS - 3 JF - Genetics TI - Estimating barriers to gene flow from distorted isolation-by-distance patterns VL - 208 ER - TY - JOUR AB - Self-incompatibility (SI) is a genetically based recognition system that functions to prevent self-fertilization and mating among related plants. An enduring puzzle in SI is how the high diversity observed in nature arises and is maintained. Based on the underlying recognition mechanism, SI can be classified into two main groups: self- and non-self recognition. Most work has focused on diversification within self-recognition systems despite expected differences between the two groups in the evolutionary pathways and outcomes of diversification. Here, we use a deterministic population genetic model and stochastic simulations to investigate how novel S-haplotypes evolve in a gametophytic non-self recognition (SRNase/S Locus F-box (SLF)) SI system. For this model the pathways for diversification involve either the maintenance or breakdown of SI and can vary in the order of mutations of the female (SRNase) and male (SLF) components. We show analytically that diversification can occur with high inbreeding depression and self-pollination, but this varies with evolutionary pathway and level of completeness (which determines the number of potential mating partners in the population), and in general is more likely for lower haplotype number. The conditions for diversification are broader in stochastic simulations of finite population size. However, the number of haplotypes observed under high inbreeding and moderate to high self-pollination is less than that commonly observed in nature. Diversification was observed through pathways that maintain SI as well as through self-compatible intermediates. Yet the lifespan of diversified haplotypes was sensitive to their level of completeness. By examining diversification in a non-self recognition SI system, this model extends our understanding of the evolution and maintenance of haplotype diversity observed in a self recognition system common in flowering plants. AU - Bodova, Katarina AU - Priklopil, Tadeas AU - Field, David AU - Barton, Nicholas H AU - Pickup, Melinda ID - 316 IS - 3 JF - Genetics TI - Evolutionary pathways for the generation of new self-incompatibility haplotypes in a non-self recognition system VL - 209 ER - TY - GEN AB - File S1 contains figures that clarify the following features: (i) effect of population size on the average number/frequency of SI classes, (ii) changes in the minimal completeness deficit in time for a single class, and (iii) diversification diagrams for all studied pathways, including the summary figure for k = 8. File S2 contains the code required for a stochastic simulation of the SLF system with an example. This file also includes the output in the form of figures and tables. AU - Bod'ová, Katarína AU - Priklopil, Tadeas AU - Field, David AU - Barton, Nicholas H AU - Pickup, Melinda ID - 9813 TI - Supplemental material for Bodova et al., 2018 ER - TY - JOUR AB - Escaping local optima is one of the major obstacles to function optimisation. Using the metaphor of a fitness landscape, local optima correspond to hills separated by fitness valleys that have to be overcome. We define a class of fitness valleys of tunable difficulty by considering their length, representing the Hamming path between the two optima and their depth, the drop in fitness. For this function class we present a runtime comparison between stochastic search algorithms using different search strategies. The (1+1) EA is a simple and well-studied evolutionary algorithm that has to jump across the valley to a point of higher fitness because it does not accept worsening moves (elitism). In contrast, the Metropolis algorithm and the Strong Selection Weak Mutation (SSWM) algorithm, a famous process in population genetics, are both able to cross the fitness valley by accepting worsening moves. We show that the runtime of the (1+1) EA depends critically on the length of the valley while the runtimes of the non-elitist algorithms depend crucially on the depth of the valley. Moreover, we show that both SSWM and Metropolis can also efficiently optimise a rugged function consisting of consecutive valleys. AU - Oliveto, Pietro AU - Paixao, Tiago AU - Pérez Heredia, Jorge AU - Sudholt, Dirk AU - Trubenova, Barbora ID - 723 IS - 5 JF - Algorithmica TI - How to escape local optima in black box optimisation when non elitism outperforms elitism VL - 80 ER - TY - JOUR AB - Adaptive introgression is common in nature and can be driven by selection acting on multiple, linked genes. We explore the effects of polygenic selection on introgression under the infinitesimal model with linkage. This model assumes that the introgressing block has an effectively infinite number of genes, each with an infinitesimal effect on the trait under selection. The block is assumed to introgress under directional selection within a native population that is genetically homogeneous. We use individual-based simulations and a branching process approximation to compute various statistics of the introgressing block, and explore how these depend on parameters such as the map length and initial trait value associated with the introgressing block, the genetic variability along the block, and the strength of selection. Our results show that the introgression dynamics of a block under infinitesimal selection is qualitatively different from the dynamics of neutral introgression. We also find that in the long run, surviving descendant blocks are likely to have intermediate lengths, and clarify how the length is shaped by the interplay between linkage and infinitesimal selection. Our results suggest that it may be difficult to distinguish introgression of single loci from that of genomic blocks with multiple, tightly linked and weakly selected loci. AU - Sachdeva, Himani AU - Barton, Nicholas H ID - 282 IS - 4 JF - Genetics TI - Introgression of a block of genome under infinitesimal selection VL - 209 ER - TY - JOUR AB - We study how a block of genome with a large number of weakly selected loci introgresses under directional selection into a genetically homogeneous population. We derive exact expressions for the expected rate of growth of any fragment of the introduced block during the initial phase of introgression, and show that the growth rate of a single-locus variant is largely insensitive to its own additive effect, but depends instead on the combined effect of all loci within a characteristic linkage scale. The expected growth rate of a fragment is highly correlated with its long-term introgression probability in populations of moderate size, and can hence identify variants that are likely to introgress across replicate populations. We clarify how the introgression probability of an individual variant is determined by the interplay between hitchhiking with relatively large fragments during the early phase of introgression and selection on fine-scale variation within these, which at longer times results in differential introgression probabilities for beneficial and deleterious loci within successful fragments. By simulating individuals, we also investigate how introgression probabilities at individual loci depend on the variance of fitness effects, the net fitness of the introduced block, and the size of the recipient population, and how this shapes the net advance under selection. Our work suggests that even highly replicable substitutions may be associated with a range of selective effects, which makes it challenging to fine map the causal loci that underlie polygenic adaptation. AU - Sachdeva, Himani AU - Barton, Nicholas H ID - 39 IS - 4 JF - Genetics SN - 00166731 TI - Replicability of introgression under linked, polygenic selection VL - 210 ER - TY - JOUR AB - Genomes of closely-related species or populations often display localized regions of enhanced relative sequence divergence, termed genomic islands. It has been proposed that these islands arise through selective sweeps and/or barriers to gene flow. Here, we genetically dissect a genomic island that controls flower color pattern differences between two subspecies of Antirrhinum majus, A.m.striatum and A.m.pseudomajus, and relate it to clinal variation across a natural hybrid zone. We show that selective sweeps likely raised relative divergence at two tightly-linked MYB-like transcription factors, leading to distinct flower patterns in the two subspecies. The two patterns provide alternate floral guides and create a strong barrier to gene flow where populations come into contact. This barrier affects the selected flower color genes and tightlylinked loci, but does not extend outside of this domain, allowing gene flow to lower relative divergence for the rest of the chromosome. Thus, both selective sweeps and barriers to gene flow play a role in shaping genomic islands: sweeps cause elevation in relative divergence, while heterogeneous gene flow flattens the surrounding "sea," making the island of divergence stand out. By showing how selective sweeps establish alternative adaptive phenotypes that lead to barriers to gene flow, our study sheds light on possible mechanisms leading to reproductive isolation and speciation. AU - Tavares, Hugo AU - Whitley, Annabel AU - Field, David AU - Bradley, Desmond AU - Couchman, Matthew AU - Copsey, Lucy AU - Elleouet, Joane AU - Burrus, Monique AU - Andalo, Christophe AU - Li, Miaomiao AU - Li, Qun AU - Xue, Yongbiao AU - Rebocho, Alexandra B AU - Barton, Nicholas H AU - Coen, Enrico ID - 38 IS - 43 JF - PNAS SN - 00278424 TI - Selection and gene flow shape genomic islands that control floral guides VL - 115 ER - TY - JOUR AB - Hanemaaijer et al. (Molecular Ecology, 27, 2018) describe the genetic consequences of the introgression of an insecticide resistance allele into a mosquito population. Linked alleles initially increased, but many of these later declined. It is hard to determine whether this decline was due to counter‐selection, rather than simply to chance. AU - Barton, Nicholas H ID - 40 IS - 24 JF - Molecular Ecology SN - 1365294X TI - The consequences of an introgression event VL - 27 ER - TY - JOUR AB - We re-examine the model of Kirkpatrick and Barton for the spread of an inversion into a local population. This model assumes that local selection maintains alleles at two or more loci, despite immigration of alternative alleles at these loci from another population. We show that an inversion is favored because it prevents the breakdown of linkage disequilibrium generated by migration; the selective advantage of an inversion is proportional to the amount of recombination between the loci involved, as in other cases where inversions are selected for. We derive expressions for the rate of spread of an inversion; when the loci covered by the inversion are tightly linked, these conditions deviate substantially from those proposed previously, and imply that an inversion can then have only a small advantage. AU - Charlesworth, Brian AU - Barton, Nicholas H ID - 565 IS - 1 JF - Genetics TI - The spread of an inversion with migration and selection VL - 208 ER - TY - JOUR AB - In this issue of GENETICS, a new method for detecting natural selection on polygenic traits is developed and applied to sev- eral human examples ( Racimo et al. 2018 ). By de fi nition, many loci contribute to variation in polygenic traits, and a challenge for evolutionary ge neticists has been that these traits can evolve by small, nearly undetectable shifts in allele frequencies across each of many, typically unknown, loci. Recently, a helpful remedy has arisen. Genome-wide associ- ation studies (GWAS) have been illuminating sets of loci that can be interrogated jointly for c hanges in allele frequencies. By aggregating small signal s of change across many such loci, directional natural selection is now in principle detect- able using genetic data, even for highly polygenic traits. This is an exciting arena of progress – with these methods, tests can be made for selection associated with traits, and we can now study selection in what may be its most prevalent mode. The continuing fast pace of GWAS publications suggest there will be many more polygenic tests of selection in the near future, as every new GWAS is an opportunity for an accom- panying test of polygenic selection. However, it is important to be aware of complications th at arise in interpretation, especially given that these studies may easily be misinter- preted both in and outside the evolutionary genetics commu- nity. Here, we provide context for understanding polygenic tests and urge caution regarding how these results are inter- preted and reported upon more broadly. AU - Novembre, John AU - Barton, Nicholas H ID - 430 IS - 4 JF - Genetics TI - Tread lightly interpreting polygenic tests of selection VL - 208 ER - TY - JOUR AB - We study the Fokker-Planck equation derived in the large system limit of the Markovian process describing the dynamics of quantitative traits. The Fokker-Planck equation is posed on a bounded domain and its transport and diffusion coefficients vanish on the domain's boundary. We first argue that, despite this degeneracy, the standard no-flux boundary condition is valid. We derive the weak formulation of the problem and prove the existence and uniqueness of its solutions by constructing the corresponding contraction semigroup on a suitable function space. Then, we prove that for the parameter regime with high enough mutation rate the problem exhibits a positive spectral gap, which implies exponential convergence to equilibrium.Next, we provide a simple derivation of the so-called Dynamic Maximum Entropy (DynMaxEnt) method for approximation of observables (moments) of the Fokker-Planck solution, which can be interpreted as a nonlinear Galerkin approximation. The limited applicability of the DynMaxEnt method inspires us to introduce its modified version that is valid for the whole range of admissible parameters. Finally, we present several numerical experiments to demonstrate the performance of both the original and modified DynMaxEnt methods. We observe that in the parameter regimes where both methods are valid, the modified one exhibits slightly better approximation properties compared to the original one. AU - Bodova, Katarina AU - Haskovec, Jan AU - Markowich, Peter ID - 607 JF - Physica D: Nonlinear Phenomena TI - Well posedness and maximum entropy approximation for the dynamics of quantitative traits VL - 376-377 ER - TY - THES AB - This thesis is concerned with the inference of current population structure based on geo-referenced genetic data. The underlying idea is that population structure affects its spatial genetic structure. Therefore, genotype information can be utilized to estimate important demographic parameters such as migration rates. These indirect estimates of population structure have become very attractive, as genotype data is now widely available. However, there also has been much concern about these approaches. Importantly, genetic structure can be influenced by many complex patterns, which often cannot be disentangled. Moreover, many methods merely fit heuristic patterns of genetic structure, and do not build upon population genetics theory. Here, I describe two novel inference methods that address these shortcomings. In Chapter 2, I introduce an inference scheme based on a new type of signal, identity by descent (IBD) blocks. Recently, it has become feasible to detect such long blocks of genome shared between pairs of samples. These blocks are direct traces of recent coalescence events. As such, they contain ample signal for inferring recent demography. I examine sharing of IBD blocks in two-dimensional populations with local migration. Using a diffusion approximation, I derive formulas for an isolation by distance pattern of long IBD blocks and show that sharing of long IBD blocks approaches rapid exponential decay for growing sample distance. I describe an inference scheme based on these results. It can robustly estimate the dispersal rate and population density, which is demonstrated on simulated data. I also show an application to estimate mean migration and the rate of recent population growth within Eastern Europe. Chapter 3 is about a novel method to estimate barriers to gene flow in a two dimensional population. This inference scheme utilizes geographically localized allele frequency fluctuations - a classical isolation by distance signal. The strength of these local fluctuations increases on average next to a barrier, and there is less correlation across it. I again use a framework of diffusion of ancestral lineages to model this effect, and provide an efficient numerical implementation to fit the results to geo-referenced biallelic SNP data. This inference scheme is able to robustly estimate strong barriers to gene flow, as tests on simulated data confirm. AU - Ringbauer, Harald ID - 200 SN - 2663-337X TI - Inferring recent demography from spatial genetic structure ER - TY - JOUR AB - Genome-scale diversity data are increasingly available in a variety of biological systems, and can be used to reconstruct the past evolutionary history of species divergence. However, extracting the full demographic information from these data is not trivial, and requires inferential methods that account for the diversity of coalescent histories throughout the genome. Here, we evaluate the potential and limitations of one such approach. We reexamine a well-known system of mussel sister species, using the joint site frequency spectrum (jSFS) of synonymousmutations computed either fromexome capture or RNA-seq, in an Approximate Bayesian Computation (ABC) framework. We first assess the best sampling strategy (number of: individuals, loci, and bins in the jSFS), and show that model selection is robust to variation in the number of individuals and loci. In contrast, different binning choices when summarizing the jSFS, strongly affect the results: including classes of low and high frequency shared polymorphisms can more effectively reveal recent migration events. We then take advantage of the flexibility of ABC to compare more realistic models of speciation, including variation in migration rates through time (i.e., periodic connectivity) and across genes (i.e., genome-wide heterogeneity in migration rates). We show that these models were consistently selected as the most probable, suggesting that mussels have experienced a complex history of gene flow during divergence and that the species boundary is semi-permeable. Our work provides a comprehensive evaluation of ABC demographic inference in mussels based on the coding jSFS, and supplies guidelines for employing different sequencing techniques and sampling strategies. We emphasize, perhaps surprisingly, that inferences are less limited by the volume of data, than by the way in which they are analyzed. AU - Fraisse, Christelle AU - Roux, Camille AU - Gagnaire, Pierre AU - Romiguier, Jonathan AU - Faivre, Nicolas AU - Welch, John AU - Bierne, Nicolas ID - 139 IS - 7 JF - PeerJ TI - The divergence history of European blue mussel species reconstructed from Approximate Bayesian Computation: The effects of sequencing techniques and sampling strategies VL - 2018 ER - TY - JOUR AB - Secondary contact is the reestablishment of gene flow between sister populations that have diverged. For instance, at the end of the Quaternary glaciations in Europe, secondary contact occurred during the northward expansion of the populations which had found refugia in the southern peninsulas. With the advent of multi-locus markers, secondary contact can be investigated using various molecular signatures including gradients of allele frequency, admixture clines, and local increase of genetic differentiation. We use coalescent simulations to investigate if molecular data provide enough information to distinguish between secondary contact following range expansion and an alternative evolutionary scenario consisting of a barrier to gene flow in an isolation-by-distance model. We find that an excess of linkage disequilibrium and of genetic diversity at the suture zone is a unique signature of secondary contact. We also find that the directionality index ψ, which was proposed to study range expansion, is informative to distinguish between the two hypotheses. However, although evidence for secondary contact is usually conveyed by statistics related to admixture coefficients, we find that they can be confounded by isolation-by-distance. We recommend to account for the spatial repartition of individuals when investigating secondary contact in order to better reflect the complex spatio-temporal evolution of populations and species. AU - Bertl, Johanna AU - Ringbauer, Harald AU - Blum, Michaël ID - 33 IS - 10 JF - PeerJ TI - Can secondary contact following range expansion be distinguished from barriers to gene flow? VL - 2018 ER - TY - JOUR AB - Pedigree and sibship reconstruction are important methods in quantifying relationships and fitness of individuals in natural populations. Current methods employ a Markov chain-based algorithm to explore plausible possible pedigrees iteratively. This provides accurate results, but is time-consuming. Here, we develop a method to infer sibship and paternity relationships from half-sibling arrays of known maternity using hierarchical clustering. Given 50 or more unlinked SNP markers and empirically derived error rates, the method performs as well as the widely used package Colony, but is faster by two orders of magnitude. Using simulations, we show that the method performs well across contrasting mating scenarios, even when samples are large. We then apply the method to open-pollinated arrays of the snapdragon Antirrhinum majus and find evidence for a high degree of multiple mating. Although we focus on diploid SNP data, the method does not depend on marker type and as such has broad applications in nonmodel systems. AU - Ellis, Thomas AU - Field, David AU - Barton, Nicholas H ID - 286 IS - 5 JF - Molecular Ecology Resources TI - Efficient inference of paternity and sibship inference given known maternity via hierarchical clustering VL - 18 ER - TY - CONF AB - There has been renewed interest in modelling the behaviour of evolutionary algorithms by more traditional mathematical objects, such as ordinary differential equations or Markov chains. The advantage is that the analysis becomes greatly facilitated due to the existence of well established methods. However, this typically comes at the cost of disregarding information about the process. Here, we introduce the use of stochastic differential equations (SDEs) for the study of EAs. SDEs can produce simple analytical results for the dynamics of stochastic processes, unlike Markov chains which can produce rigorous but unwieldy expressions about the dynamics. On the other hand, unlike ordinary differential equations (ODEs), they do not discard information about the stochasticity of the process. We show that these are especially suitable for the analysis of fixed budget scenarios and present analogs of the additive and multiplicative drift theorems for SDEs. We exemplify the use of these methods for two model algorithms ((1+1) EA and RLS) on two canonical problems(OneMax and LeadingOnes). AU - Paixao, Tiago AU - Pérez Heredia, Jorge ID - 1112 SN - 978-145034651-1 T2 - Proceedings of the 14th ACM/SIGEVO Conference on Foundations of Genetic Algorithms TI - An application of stochastic differential equations to evolutionary algorithms ER - TY - JOUR AB - Variation in genotypes may be responsible for differences in dispersal rates, directional biases, and growth rates of individuals. These traits may favor certain genotypes and enhance their spatiotemporal spreading into areas occupied by the less advantageous genotypes. We study how these factors influence the speed of spreading in the case of two competing genotypes under the assumption that spatial variation of the total population is small compared to the spatial variation of the frequencies of the genotypes in the population. In that case, the dynamics of the frequency of one of the genotypes is approximately described by a generalized Fisher–Kolmogorov–Petrovskii–Piskunov (F–KPP) equation. This generalized F–KPP equation with (nonlinear) frequency-dependent diffusion and advection terms admits traveling wave solutions that characterize the invasion of the dominant genotype. Our existence results generalize the classical theory for traveling waves for the F–KPP with constant coefficients. Moreover, in the particular case of the quadratic (monostable) nonlinear growth–decay rate in the generalized F–KPP we study in detail the influence of the variance in diffusion and mean displacement rates of the two genotypes on the minimal wave propagation speed. AU - Kollár, Richard AU - Novak, Sebastian ID - 1191 IS - 3 JF - Bulletin of Mathematical Biology TI - Existence of traveling waves for the generalized F–KPP equation VL - 79 ER - TY - JOUR AB - Most phenotypes are determined by molecular systems composed of specifically interacting molecules. However, unlike for individual components, little is known about the distributions of mutational effects of molecular systems as a whole. We ask how the distribution of mutational effects of a transcriptional regulatory system differs from the distributions of its components, by first independently, and then simultaneously, mutating a transcription factor and the associated promoter it represses. We find that the system distribution exhibits increased phenotypic variation compared to individual component distributions - an effect arising from intermolecular epistasis between the transcription factor and its DNA-binding site. In large part, this epistasis can be qualitatively attributed to the structure of the transcriptional regulatory system and could therefore be a common feature in prokaryotes. Counter-intuitively, intermolecular epistasis can alleviate the constraints of individual components, thereby increasing phenotypic variation that selection could act on and facilitating adaptive evolution. AU - Lagator, Mato AU - Sarikas, Srdjan AU - Acar, Hande AU - Bollback, Jonathan P AU - Guet, Calin C ID - 570 JF - eLife SN - 2050084X TI - Regulatory network structure determines patterns of intermolecular epistasis VL - 6 ER - TY - JOUR AB - Small RNAs (sRNAs) regulate genes in plants and animals. Here, we show that population-wide differences in color patterns in snapdragon flowers are caused by an inverted duplication that generates sRNAs. The complexity and size of the transcripts indicate that the duplication represents an intermediate on the pathway to microRNA evolution. The sRNAs repress a pigment biosynthesis gene, creating a yellow highlight at the site of pollinator entry. The inverted duplication exhibits steep clines in allele frequency in a natural hybrid zone, showing that the allele is under selection. Thus, regulatory interactions of evolutionarily recent sRNAs can be acted upon by selection and contribute to the evolution of phenotypic diversity. AU - Bradley, Desmond AU - Xu, Ping AU - Mohorianu, Irina AU - Whibley, Annabel AU - Field, David AU - Tavares, Hugo AU - Couchman, Matthew AU - Copsey, Lucy AU - Carpenter, Rosemary AU - Li, Miaomiao AU - Li, Qun AU - Xue, Yongbiao AU - Dalmay, Tamas AU - Coen, Enrico ID - 611 IS - 6365 JF - Science SN - 00368075 TI - Evolution of flower color pattern through selection on regulatory small RNAs VL - 358 ER -