bioRxiv preprint doi: https://doi.org/10.1101/2020.11.20.391151; this version posted March 13, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
1 Molecular signatures of resource competition: clonal interference 2 favors ecological diversification and can lead to incipient speciation
3 Massimo Amicone+, * and Isabel Gordo+, *
4 +Instituto Gulbenkian de Ciência (IGC)
5 *corresponding authors: [email protected], [email protected] 6
7 Short running title 8 Clonal interactions drive ecotypes formation. 9 Author contributions 10 MA and IG designed the study, all authors wrote the manuscript and provided final 11 approval for publication.
12 Conflict of interest 13 The authors declare no conflict of interest.
14 Acknowledgements 15 The authors acknowledge C Bank and the members of the Evolutionary Dynamics and 16 Evolutionary Biology labs of the IGC for their assistance throughout the development of 17 this work; PR Campos, T Paixão, S Miller, G Sgarlata and R Ramiro for their comments 18 on the manuscript. This work was supported by Portuguese Science Foundation (FCT) 19 Grant PTDC/BIA-EVL/31528/2017; Deutsche Forschungs-gemeinschaft (DFG) Grant 20 SFB 1310; and a cooperation agreement between University of Cologne and Gulbenkian 21 Institute. MA was supported by FCT Grant PD/BD/138735/2018. 22 23 Data archival location will be made public upon acceptance.
24
1 bioRxiv preprint doi: https://doi.org/10.1101/2020.11.20.391151; this version posted March 13, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
25 Abstract
26 Microbial ecosystems harbor an astonishing diversity that can persist for long times. To 27 understand how such diversity is structured and maintained, ecological and 28 evolutionary processes need to be integrated at similar timescales. Here, we study a 29 model of resource competition that allows for evolution via de novo mutation, and focus 30 on rapidly adapting asexual populations with large mutational inputs, as typical of many 31 bacteria species. We characterize the adaptation and diversification of an initially 32 maladapted population and show how the eco-evolutionary dynamics are shaped by the 33 interaction between simultaneously emerging lineages – clonal interference. We find 34 that in large populations, more intense clonal interference fosters diversification under 35 sympatry, increasing the probability that phenotypically and genetically distinct 36 clusters stably coexist. In smaller populations, the accumulation of deleterious and 37 compensatory mutations can push further the diversification process and kick-start 38 speciation. Our findings have implications beyond microbial populations, providing 39 novel insights about the interplay between ecology and evolution in clonal populations.
40 Keywords: Eco-evolutionary dynamics, clonal interference, resource competition, 41 competitive exclusion, diversification, speciation. 42
2 bioRxiv preprint doi: https://doi.org/10.1101/2020.11.20.391151; this version posted March 13, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
43 Introduction
44 Understanding the mechanisms behind the evolution of biodiversity and the formation 45 of communities remains a difficult challenge. One must integrate ecology and evolution 46 over similar timescales, as taken together, they can give rise to phenomena that could 47 not be explained by either alone (Schoener, 2011). The competitive exclusion principle 48 (first stated by Hardin, 1960) theoretically binds the number of species by the number 49 of limiting resources. This principle generated an apparent contradiction between 50 theoretical expectations and observations, often referred to as the “paradox of the 51 plankton” (Hutchinson, 1961). In fact, ecosystems can be replete of diversity even in 52 limiting environments, both in nature (Hutchinson, 1961; Tilman, 1982; Huston, 1994) 53 and in more controlled laboratory conditions (Maharjan et al., 2006; Gresham et al., 54 2008; Kinnersley, Holben and Rosenzweig, 2009; Herron and Doebeli, 2013; Good et al., 55 2017). Different theoretical approaches have been adopted to resolve such controversy. 56 One approach is to assume an existing diversity and identify mechanisms that can 57 maintain it. Following this, several ecological properties were proposed to maintain 58 diversity, including heterogeneity in space (Abrams, 1988) and time (Litchman and 59 Klausmeier, 2001), trade-offs on the species’ traits (Posfai, Taillefumier and Wingreen, 60 2017), or gene regulation (Pacciani-Mori et al., 2020). Here we explore a classical 61 model of competition for resources, where extensive diversity can be maintained by a 62 metabolic trade-off (Posfai, Taillefumier and Wingreen, 2017), and ask a different 63 question: in an initially monomorphic population, what diversity can be generated and 64 maintained if the species’ traits continuously evolve?
65 In eco-evolutionary frameworks, mutations generate new genetic variants whose fate 66 depends on the state of the ecosystem and, in turn, their increase in frequency can alter 67 the populations. A common outcome of such eco-evolutionary feedbacks is that 68 evolution limits diversity by reducing the effectiveness of coexistence mechanisms 69 (Edwards et al., 2018). The diversity that would be possible by ecological principles 70 alone, is reduced by selection of the fittest and competitive exclusion. Several studies 71 have produced novel understanding on the evolution of diversity (Dieckmann and 72 Doebeli, 1999; Shoresh, Hegreness and Kishony, 2008; Doebeli, 2011; Kremer and 73 Klausmeier, 2017), but the majority rely on the strong-selection-weak-mutation 74 assumption, i.e. on a timescale separation between ecological and evolutionary
3 bioRxiv preprint doi: https://doi.org/10.1101/2020.11.20.391151; this version posted March 13, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
75 processes (but see Farahpour et al., 2018). The emergence of mutations is assumed to 76 be much slower than the ecological dynamics, thus, before a new lineage arises, the 77 population has already reached ecological equilibrium. While allowing for analytical 78 tractability the weak mutation assumption comes at a cost: it neglects the overlap 79 between multiple evolving lineages – clonal interference. Clonal interference has been 80 extensively observed in microbial communities in vitro and in vivo (Desai, Fisher and 81 Murray, 2007; Barroso-Batista et al., 2014) and can occur under different regimes of 82 intensity: a weak regime where a few lineages compete for fixation (Philip J. Gerrish & 83 Richard E. Lenski, 1998; Billiard and Smadi, 2020) or one where many different 84 haplotypes segregate (Good et al., 2012). Population genetics models incorporating 85 clonal interference have generated predictions for the adaptation rate, fixation 86 probabilities and genetic diversity in a population (Gerrish and Lenski, 1998; Park and 87 Krug, 2007; Good et al., 2012; de Sousa et al., 2016), yet they typically ignore ecological 88 interactions (but see Good, Martis and Hallatschek, 2018).
89 When the population size (N) and/or the mutation rate to new beneficial mutations (Ub)
90 are not small (NUb >>1), multiple lineages can increase in frequency simultaneously, 91 ecologically interact with each other and evolve in non-trivial ways. Although these 92 processes are inevitably intertwined in real ecosystems (Lawrence et al., 2012; Barroso- 93 Batista et al., 2014, 2020; Garud et al., 2019), theoretical work is still needed to 94 investigate how they act in chorus.
95 Good and colleagues (Good, Martis and Hallatschek, 2018) have recently developed a 96 theoretical work that incorporates ecological and evolutionary mechanisms via a 97 combination of frequency-dependent and directional selection. Their eco-evolutionary 98 model is able to reproduce empirical patterns of co-existence and fixation of new 99 mutations in experimentally evolved clonal populations (Good et al., 2017). It also 100 shows how diversification depends on the ratio between the rates of strategy mutations 101 and unconditionally beneficial mutations. However, in order to derive analytical 102 expression for the eco-evolutionary dynamics, the authors have focused on the weak
103 mutation limit NUb <<1 and only briefly investigated clonal interference.
104 Here, we study a similar model to that of (Posfai, Taillefumier and Wingreen, 2017; 105 Good, Martis and Hallatschek, 2018) but assume a trade-off that only affects well 106 adapted genotypes (fitness-dependent trade-off) and conduct a more systematic
4 bioRxiv preprint doi: https://doi.org/10.1101/2020.11.20.391151; this version posted March 13, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
107 simulation study of the different mutation regimes, including extensive clonal 108 interference, in large and small populations. We follow an initially isogenic population 109 throughout time and characterize the patterns of adaptation at both phenotypic and 110 genetic levels, by common statistics used to analyze molecular evolution data. We 111 focused only on mutations that affect the ability of consuming resources (instead of 112 differentiating between strategy and unconditionally beneficial mutations) and we do 113 not impose restrictive assumptions on mutation rates and on timescales, as common in 114 other models (Geritz et al., 1998; Shoresh, Hegreness and Kishony, 2008; Good, Martis 115 and Hallatschek, 2018). Albeit at the cost of analytical tractability, our approach is to 116 focus on the phenomena that emerge from the stochastic simulations where many more 117 lineages compete, compared to previous studies (Billiard and Smadi, 2020).
118 We find that: high levels of intra-specific variation can be generated and maintained via 119 a balance between selection and mutation; functionally distinct clusters of genotypes – 120 ecotypes - can emerge and stably coexist; and the interaction between large mutational 121 inputs and the energetic trade-off can lead to incipient speciation. Taken together, our 122 results describe how clonal populations can give rise to extensive diversity and 123 establish a first form of community, even in simple and constant environments.
124
125 Model and Methods
126 Eco-Evolutionary model
127 We model the dynamics of a single clonal lineage evolving to consume a set of different 128 substitutable resources, constantly replenished in a well-mixed environment (Fig.1A). 129 Individuals mutate at a rate U (per-genome, per-generation rate of non-lethal 130 mutations) and the fate of the emerging mutations depends on their phenotypic effects, 131 on the resource concentration, on the other individuals present in the environment and 132 on drift.
133 The underlying dynamics are based on the MacArthur’s consumer resource model (Mac 134 Arthur, 1969), recently formalized to explain high levels of diversity in the presence of a 135 metabolic trade-off (Posfai, Taillefumier and Wingreen, 2017) and further extended to 136 study adaptation (Good, Martis and Hallatschek, 2018). Briefly, let M be the number of
5 bioRxiv preprint doi: https://doi.org/10.1101/2020.11.20.391151; this version posted March 13, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
137 types present at time t with densities (#cells/V) ni, (i=1…M) and R the number of
138 substitutable resources with input concentrations rj, (j=1…R). The expected density 139 dynamics of each type are: