|
Copyright © 2005, Genetics Society of America Exploring the Evolution of Wolbachia Compatibility Types A Simulation Approach 1Corresponding author: Biology Department, University College London, 4 Stephenson Way, London, NW1 2HE, United Kingdom. E-mail: s.charlat/at/ucl.ac.uk Communicating editor: H. G. Spencer Received March 14, 2003; Accepted March 7, 2005. ![]() | |||||||||
Abstract Wolbachia-induced cytoplasmic incompatibility (CI) is observed when males bearing the bacterium mate with uninfected females or with females bearing a different Wolbachia variant; in such crosses, paternal chromosomes are lost at the first embryonic mitosis, most often resulting in developmental arrest. The molecular basis of CI is currently unknown, but it is useful to distinguish conceptually the male and female sides of this phenomenon: in males, Wolbachia must do something, before it is shed from maturing sperm, that will disrupt paternal chromosomes functionality [this is usually termed “the modification (mod) function”]; in females, Wolbachia must somehow restore embryonic viability, through what is usually called “the rescue (resc) function.” The occurrence of CI in crosses between males and females bearing different Wolbachia variants demonstrates that the mod and resc functions interact in a specific manner: different mod resc pairs make different compatibility types. We are interested in the evolutionary process allowing the diversification of compatibility types. In an earlier model, based on the main assumption that the mod and resc functions can mutate independently, we have shown that compatibility types can evolve through a two-step process, the first involving drift on mod variations and the second involving selection on resc variations. This previous study has highlighted the need for simulation-based models that would include the effects of nondeterministic evolutionary forces. This study is based on a simulation program fulfilling this condition, allowing us to follow the evolution of compatibility types under mutation, drift, and selection. Most importantly, simulations suggest that in the frame of our model, the evolution of compatibility types is likely to be a gradual process, with new compatibility types remaining partially compatible with ancestral ones. | |||||||||
MATERNALLY inherited elements are subject to sex-dependent selective pressures: their fitness is increased if females produce more females or better surviving females, regardless of possible detrimental effects to males (Cosmides and Tooby 1981; Frank and Hurst 1996). The endocellular bacterium Wolbachia illustrates nicely the possible outcomes of such selection, having evolved a variety of “sex manipulation strategies” that can be interpreted within this theoretical frame (reviewed in O'Neill et al. 1997; Stouthamer et al. 1999). Cytoplasmic incompatibility (CI) is one of them, probably the most common (reviewed in Hoffmann and Turelli 1997; Charlat et al. 2001a; Bourtzis et al. 2003). In embryos resulting from crosses between males that bear a CI Wolbachia and females that do not, paternal chromosomes are lost at the first mitosis (Callaini et al. 1996, 1997; Lassy and Karr 1996; Tram and Sullivan 2002), resulting in death in most cases (that is, in diploid species), but less often in male development (in some haplo-diploids). In contrast, if the female bears the bacterium, paternal chromosomes are not lost. Infected females thus produce on average more females than do uninfected ones, allowing infected cytoplasmic lines to invade uninfected populations. Infected males suffer a fertility deficit if uninfected females remain in the populations, since some proportion of their mating will be partially or fully sterile. But Wolbachia is not affected, as it is transmitted by females only. The underlying mechanism remains to be elucidated. In males, Wolbachia must somehow affect the paternal nucleus before it is shed from maturing sperm, resulting in paternal chromosome loss after fertilization. In females, the bacterium must somehow prevent this loss and thereby rescue the embryo. This conceptual distinction between the male and female sides of CI was formalized by Werren (1997) through the modification/rescue (mod/resc) terminology. Interestingly, crosses between infected males and infected females can also be incompatible, if the two partners bear different bacteria. The mod and resc functions thus seem to interact specifically: different Wolbachia can harbor different mod resc pairs, that is, different compatibility types. We are interested in the process that allows compatibility types to evolve. In an earlier study (Charlat et al. 2001b), we showed that compatibility types are not constrained by stabilizing selection, if mod and resc are determined by different genes, which, we think, is a reasonable assumption (Poinsot et al. 2003; see also Kose and Karr 1995; Callaini et al. 1997; Tram and Sullivan 2002). We suggested that compatibility types could change through a two-step process: the first involving drift on mod variations and the second involving selection on resc variations. This work highlighted the need for simulation-based models that would incorporate the effects of nondeterministic evolutionary forces. The present article is based on a simulation program developed in an attempt to fulfill this condition, allowing us to follow the evolution of compatibility types under mutation, drift, and selection. | |||||||||
THE MODEL What defines a Wolbachia variant: Attempts have been made in earlier literature to translate the mod resc general formalization into more concrete models (reviewed in Poinsot et al. 2003): namely, the “slow-motion” model (Callaini et al. 1997; Tram and Sullivan 2002), the “titration-restitution” model (Kose and Karr 1995), and the “lock-and-key” model, more or less explicitly proposed in several articles (Breeuwer and Werren 1990; Hurst 1991; Werren 1997; Poinsot and Merçot 1999). In this later model mod and resc are seen as a lock and a key (that is, interacting physically and specifically with each other). In fertilized embryos, the lock, fixed on paternal material, would come into contact with the key, produced by Wolbachia in the egg. Depending on the conformation of the lock and the key, compatibility would range from 0 (total incompatibility) to 1 (total compatibility). In our opinion, this model is currently the most parsimonious and satisfactory (for a more detailed discussion of this issue see Poinsot et al. 2003). This view was implicitly the basis of our earlier theoretical work (Charlat et al. 2001b). The symbolism used here refers to the lock-and-key model more explicitly.A Wolbachia variant is defined here by five parameters (summarized in Table 1). Two parameters define the compatibility type: the lock and the key. In practice, lock and key are modeled as sequences of S sites, with n possible states for each site (states 1, 2, … , n). In a cross between a male bearing Wolbachia i and a female bearing Wolbachia j, a compatibility score (Cij) is calculated as the proportion of matching sites in the lock/key alignment. For example, with S = 10, we can have locki = 2222211111 and keyj = 2222222222, which gives the compatibility score Cij = 0.5. Lock and key are comparable to the modC and rescC parameters in our earlier study (Charlat et al. 2001b), but with the notable difference that the current model allows partial compatibility if S > 1. MI (mod intensity) is the efficiency of the mod function; it corresponds to what is often referred to as “CI level,” which can be measured in crosses between infected males and uninfected females. Biologically, MI can be seen as the proportion of infected males' sperm that is actually affected by CI. TE (transmission efficiency) is the proportion of infected eggs among those laid by infected females. Finally, FE (fitness effect) is the fecundity of infected females relative to that of uninfected females.
The evolution of the TE and FE parameters is not in the focus of this study, as their evolutionary trajectories have been well described by the analytical approach (Turelli 1994): selection on Wolbachia always acts to increase infected females fecundity and maternal transmission efficiency. In other words, starting from any initial condition, these parameters will rapidly reach their maximum values under mutation and selection. However, FE and TE are still relevant to the model, as we can investigate the effect of setting upper limits <1 for these two parameters. For example, the model allows us to analyze the evolution of the lock and key parameters in populations where TE cannot exceed 0.9, that is, in populations where uninfected individuals can persist. Similarly, the evolution of the MI parameter in panmictic populations has been previously worked out (Prout 1994; Turelli 1994): with random mating, MI is not under selection and evolves through drift only. However, keeping MI constant would impede the “realism” of our analysis. Consequently, the model includes the evolution of MI, although we do not focus on this aspect. From one generation to the next: Ne males and Ne females actually reproduce at every generation; Ne is thus the effective population size of cytoplasmic genes. Ne females and Ne males are randomly chosen as parents on the basis of their frequencies at generation x. Sampling errors allow frequencies to drift. For every parental pair, progeny is determined as illustrated in Figure 1, allowing natural selection to act. The Wolbachia variants present at generation x + 1 are then submitted to mutation. Thus, generations x and x + 1 are separated by a round of drift, selection, and mutation. The five above-listed parameters are allowed to mutate independently from each other. Mutation rate per generation is denoted Mu. For the lock and key parameters, Mu is multiplied by S so that Mu gives the mutation rate per generation per site while Mu × S gives the mutation rate per generation per lock (or key) sequence. Mutations affecting the lock or the key sequences result in changing one of the S sites from its state to one of the n − 1 other states. Mutations affecting the MI, TE, or FE parameters result in adding +0.05 or −0.05 with equal probability to their initial value; if the initial value is maximum, then the mutation results in adding +0.00 or −0.05 with equal probability; conversely, if the initial value is minimum, then the mutation results in adding +0.00 or +0.05 with equal probability.
Assumptions: We are aware of the following assumptions in the model: (1) unbiased sex ratio; (2) nonoverlapping generations; (3) no population structure; (4) random mating; (5) no multiple infections [a given individual host is homogeneous with regard to Wolbachia infections; when a mutation gives rise to a new variant, its host is infected by this clone only (see Dobson 2004 for an alternative model relaxing this assumption)]; (6) no recombination between Wolbachia variants; (7) no gene duplication (one bacterium bears only one mod and one resc function); and (8) no variations of host effects. | |||||||||
VALIDATING THE MODEL To validate the model, we tested whether it could retrieve earlier results, derived from the analytical approach. We first verified the basic prediction that CI allows Wolbachia to invade uninfected host populations, and more efficiently so if MI is high (Figure 2). We then investigated the combined effects of MI, TE, and FE on invasion dynamics. Caspari and Watson (1959), Fine (1978), and Hoffmann et al. (1990) showed that if TE and/or FE are <1, Wolbachia does not invade unless it first reaches a threshold frequency depending on MI, TE, and FE. Above that point, infection frequency increases toward a stable infection frequency, which is not fixation if MI and TE are <1. As illustrated in Figure 3, the simulation and analytical approaches provide congruent predictions.
| |||||||||
ILLUSTRATING RELEVANT PROCESSES WITHOUT RANDOM MUTATION Before considering “realistic evolution,” where an initial population can freely change under mutation, selection, and drift, we present here the results of simulations conducted without random mutation: several variants with specific properties are initially introduced and their frequencies followed over 1000 generations. These particular cases will allow the reader to understand which sequence of events can lead to which population state. With these processes in mind, the evolution of populations is analyzed more realistically in the next section. Infection loss: Analytical models have revealed that elevated values of the MI, TE, and FE parameters facilitate the stable maintenance of Wolbachia in host populations (Hoffmann et al. 1990). They further showed that mutations decreasing TE or FE are always selected against (Turelli 1994), so that the long-term evolution of these two parameters should stabilize the presence of Wolbachia. In contrast, mutations reducing MI are not selected against in panmictic populations (Prout 1994; Turelli 1994), so that MI is supposed to evolve under drift only (unless it is linked to other traits through pleiotropic effects). Increasing MI will stabilize the infection, but decreasing MI will have the opposite effect. Figure 4A illustrates how a random decrease of the average MI in the population can lead to infection loss: when Wolbachia variants with low MI get too frequent, the overall infection frequency (the stable equilibrium predicted by analytical models) decreases, while the threshold infection frequency (below which Wolbachia is lost deterministically) increases. Eventually, the population can get out of the conditions under which infection is maintained.
Less explicit in earlier analyses is the fact that variations affecting the lock parameter can also lead to infection loss. As detailed elsewhere (Charlat et al. 2001b), mutations affecting the lock sequence are not subject to selection, although they can give rise to self-incompatible, or “suicidal” Wolbachia. Figure 4B illustrates how this can lead to infection loss: as new lock variants (totally or partially self-incompatible) get too frequent by drift, the fitness gain provided by CI to infected females is lowered, which eventually leads to infection loss. In summary, neutral variations of MI make CI less deleterious to uninfected females, while neutral variations of the lock make CI more deleterious to infected females. The final effect is the same in the two cases: net benefit to infected cytoplasmic lines is reduced. Fixation of new compatibility types: As detailed elsewhere (Charlat et al. 2001b), random variations of lock can create the conditions for new compatibility types to invade populations. For the purpose of this section, let us define two different lock sequences, lockA (1111111111) and lockB (2222222222), and two different key sequences, keyA (1111111111) and keyB (2222222222). Consider a population harboring two Wolbachia variants: lockA/keyA (AA) and lockB/keyA (BA). This is a neutral polymorphism because the two variants harbor the same key sequence (keyA). The relative proportion of the two variants thus changes through drift only. A third variant [lockB/keyB (BB)] is introduced in the population. This new variant gets more frequent if the overall frequency of lockB exceeds that of lockA, that is, if f(BA) + f(BB) > f(AA). As f(BB) increases, lockB becomes more frequent, so that the fitness of BB increases. Eventually, BB will get fixed, so that the compatibility type will have evolved from lockA/keyA to lockB/keyB. Figure 5 shows how such a process can be visualized with our model.
Balanced suicidal polymorphism: Consider a population including two variants: lockA/keyA (AA) and lockB/keyA (BA). This is a neutral polymorphism: the relative proportion of AA and BA changes through drift only (Charlat et al. 2001b). Consider a new variant (lockA/keyB or AB) arising by mutation of the key in the AA variant. AB gets more frequent if the overall frequency of lockB exceeds that of lockA, that is, if f(BA) > f(AA) + f(AB). However, as AB becomes more frequent, lockB becomes less frequent, so that selection for AB is reduced. This leads to a situation where the three variants have the same fitness, when the frequency of lockA equals that of lockB.When such an equilibrium is reached, selection does not favor any Wolbachia variant in particular. However, simulation runs such as that presented in Figure 6 show that not all variants are maintained at stable frequencies. The frequency of BA appears to be stable, but f(AA) and f(AB) vary randomly and symmetrically, leading to the loss of either AB (Figure 6A) or AA (Figure 6B). If AB is lost, the population goes back to the initial neutral polymorphism with AA and BA. In contrast, if AA is lost, the population reaches a stable polymorphism, with BA and AB at equal frequencies.
Notably, the initial conditions (the relative frequencies of BA and AA when AB is introduced) strongly affect the respective likelihoods of the outcomes described in Figure 6, A and B. The higher the initial f(BA)/f(AA) ratio is, the higher is the frequency reached by the mutant AB under selection (in the very first generations following its introduction) and, in turn, the more likely is the outcome illustrated in Figure 6B, which depends directly on the f(AA)/f(AB) ratio at the beginning of the “drift of AA and AB” stage (Figure 6, A and B, left). A full understanding of the frequency variations during the drift of AA and AB stage is not straightforward (and not necessary for the following sections). The main point here is to show that a balanced suicidal polymorphism (illustrated in Figure 6B) can occur. The sequence of events leading to such balanced polymorphism involves a complex interaction between selection and drift: AA and AB are globally neutral relative to each other (although locally selection takes place) while the frequency of BA is locally and globally stabilized by selection. Figure 7 provides a more detailed explanation to the interested reader.
Neutral suicidal polymorphism and population extinction: Neutral variations of the lock function can greatly reduce host population mean fitness: as a lockB/keyA type gets frequent by drift, many crosses in the population are incompatible. Eventually, nothing opposes the fixation by drift of such a suicidal Wolbachia. If MI < 1, some proportion of the eggs can still survive. However, if MI = 1, and if lockB and keyA are totally incompatible, no viable progeny is produced, so that the host population simply goes extinct. | |||||||||
EVOLUTION UNDER MUTATION, DRIFT, AND SELECTION Typology on population states: Having illustrated several possible sequences of events using specific combinations of lock/key pairs without random mutation, we now analyze how these different processes globally affect the evolution of populations, by allowing all parameters to mutate randomly. For the purpose of this analysis, let us define a typology on “population states,” with six discrete states distinguished:
This typology being defined, the evolution of populations can be depicted as transitions between population states over generations. Starting from an initial condition, different simulations can lead to very different patterns, because of random mutation and drift. To describe general tendencies, we can compile the results of a sufficient number of simulation runs, which will allow us to estimate the probability of the different population states over generations. In the following sections, we use this approach to explore the evolution of compatibility types under various conditions. Preliminary analyses (not shown) have been performed using high mutation rates over a short number of generations, to moderate calculation time. These have allowed us to identify four important factors: (1) the length of the lock and key sequences (parameter S); (2) the mutation rate (parameter Mu); (3) the population size (parameter Ne); and (4) the upper limits for TE, FE, and MI (parameters TEmax, FEmax, and MImax) that will condition the possible maintenance of uninfected individuals in the populations. The effect of these factors is now examined in detail. The effect of the lock/key structure: We defined lock and key as sequences of S sites. In a cross between a male bearing Wolbachia i and a female bearing Wolbachia j, the compatibility score (Cij) is calculated as the proportion of matching sites in the lock/key alignment. To investigate the effect of varying S on the evolution of compatibility types, we followed 500 populations over 100,000 generations under two different conditions: S = 1 (Figure 8A) and S = 10 (Figure 8B). In these simulations, Ne = 103 and Mu = 10−6 (see Figure 8 legend for details on other parameters).
When S = 1 (Figure 8A), we observe that the probability of the initial configuration (state 3: predominance of the initial type) decreases very slowly over generations. Indeed, after 100,000 generations, predominance of the initial type is still observed in 90% of the simulations. Furthermore, we observe that among the remaining 10%, most populations have gone extinct due to fixation of a self-incompatible bacterium (state 2: population extinction). Thus, it appears that under these conditions, new compatibility types do not evolve. As illustrated in Figure 8B, things are clearly different when S = 10. Indeed, after 100,000 generations, <40% of the populations are still in the initial configuration, and the remaining 60% are in either state 4 (predominance of a new compatibility type) or state 6 (neutral polymorphism or predominance of partially suicidal Wolbachia). To apprehend the rationale behind these effects of S variations, one must distinguish two aspects of the differences observed between Figure 8A and 8B: (1) the slow vs. rapid decrease of state 3 (the initial configuration) and (2) the replacement of state 3 by state 2 (population extinction) vs. states 4 and 6. The explanation of difference 1 is the following: the mutation rate of the lock and key sequences is defined here as 10−6 per site; in other words, the overall mutation rate of the lock and key sequences is lower when the sequence is short (Figure 8A) than when it is long (Figure 8B); the initial configuration is lost faster in Figure 8A because the overall mutation rates of the lock and key sequences are higher. The explanation of difference 2 (population extinction vs. other states) is less straightforward. In Figure 8A, where S = 1, compatibility between a given lock and a given key can be only 0 or 1. In other words, if the intensity of mod is maximum (that is, if MI = 1), which is initially the case in these simulations, not a single viable egg is produced in incompatible crosses. Thus, if a suicidal type (with lock = 2 and key = 1) is fixed by drift, the host population goes extinct. In contrast, in Figure 8B, where S = 10, mutations of the lock and key sequences lead only to partial incompatibility. In other words, a suicidal type (partially suicidal) can reach fixation without leading the host population to extinction. In summary, the fact that populations go extinct in Figure 8A suggests that mutated lock sequences often reach fixation before a compatible key occurs by mutation. This readily leads the host population to extinction if S = 1, that is, if mutations of the lock sequence result in complete incompatibility. We investigated in a similar way the effect of varying the n parameter (the number of possible states at each site in the lock and key sequences). The results are not presented graphically, as only minor quantitative effects were observed. Let us note that simply increasing the value of n tends to slow down the evolution of new compatibility types. This can be illustrated using the following example. Consider a population where a suicidal type (e.g., with lock = 2 and key = 1) is frequent. If n = 2, mutations of the key sequence will give rise to the appropriate key (that is key = 2) with probability 1. In contrast, if n = 10, the appropriate mutation will occur with probability 1/9 only. In other words, increasing the value of n reduces the probability of favorable mutations to occur in the key sequence, resulting in slower evolution of compatibility types. In all simulations presented here, n is arbitrarily set at 10, a value that appeared to be a reasonable compromise between realism and computation time constraints. The effect of mutation rates: The results discussed so far were obtained with Mu = 10−6. To investigate the effects of mutation rates on the evolution of compatibility types, we repeated similar simulations to those presented in Figure 8, A and B (that is, with two different S values), under a 10 times higher mutation rate (Mu = 10−5). The results are presented in Figure 8, C and D. By comparing Figure 8A and 8B with 8C and 8D, respectively, one can assess the effect of varying mutation rates for two different S values. In these simulations, Ne = 103 and Mu = 10−5 (see Figure 8 legend for details on other parameters).Three notable differences can be seen between Figure 8A and 8C: (1) in Figure 8C, <40% of populations are still in state 3 after 100,000 generations, as compared to 90% in Figure 8A; (2) in Figure 8C, >30% of the populations have gone extinct (state 2) after 100,000 generations, as compared to 10% in Figure 8A; and finally, (3) in Figure 8C, >20% of the populations harbor a new lock/key pair after 100,000 generations, as compared to 1% in Figure 8A. Difference 1 (that is, state 3 is less frequent in Figure 8C than in Figure 8A) illustrates that the initial lock/key type is lost faster when mutation rate is higher, as expected. It is notable that the rate of decrease of state 3 in Figure 8C is exactly the same as that observed in Figure 8B. This is consistent with the interpretation we gave when comparing Figure 8A with 8B: state 3 is lost faster when S is larger simply because the overall mutation rate of the lock and key sequences is increased. Difference 2 (that is, state 2 is more frequent in Figure 8C than in Figure 8A) suggests that extinction risk increases with mutation rate. This can be understood by remembering that extinction results from the fixation of a suicidal type (e.g., with lock = 2 and key = 1). Mutations giving rise to such suicidal types being neutral, their rate of fixation depends only on the mutation rate (Kimura 1983). Difference 3 (that is, state 4 is more frequent in Figure 8C than in Figure 8A) suggests that the evolution of compatibility types is facilitated by increased mutation rates, which might seem straightforward. A subtle aspect deserves, however, to be discussed: in Figure 8C, the ratio state 4/state 2 (new type/extinction) is much higher than that in Figure 8A. To understand why, one must remember that spreading of new type (e.g., with lock = 2 and key = 2) can occur only in populations where a suicidal type (e.g., with lock = 2 and key = 1) is sufficiently frequent (see Figure 5). When the mutation rate is increased (as in Figure 8C), mutations giving rise to this new type are more likely to occur before the suicidal type has reached fixation (that is, before populations have gone extinct). Consider now Figure 8, B and D. Three notable differences can be seen: (1) in Figure 8D, only a tiny remnant of populations are still in state 3 after 50,000 generations, as compared to 60% in Figure 8B; (2) in Figure 8D, the frequency of state 4 reaches a plateau of 60% after 50,000 generations, while no plateau was reached in Figure 8B after 100,000 generations; and finally, (3) in Figure 8D, a small but significant and stable proportion of populations is in state 5 (balanced suicidal polymorphism), while this state was not observed in Figure 8B. Difference 1 (that is, the complete loss of state 3 in Figure 8D) suggests that state 3 (predominance of the initial compatibility type) is unstable in the long term: mutation rate in Figure 8D is high enough for its inescapable loss to be observed after 50,000 generations only. Difference 2 (that is, state 4 reaches a plateau in Figure 8D) suggests that even with high mutation rates, the proportion of populations harboring a new compatibility type can never reach 1. Difference 3 (that is, state 5 is more frequent in Figure 8D) suggests that high mutation rates facilitate the occurrence of balanced suicidal polymorphism. Overall, it is interesting to note that the conditions used in Figure 8D (Mu = 10−5 and S = 10) are such that 50,000 generations seem to be sufficient for the equilibrium distribution of population states to be observed. At equilibrium, populations are in state 4 (predominance of new type), state 5 (balanced suicidal polymorphism), or state 6 (neutral polymorphism or predominance of partially suicidal Wolbachia). This equilibrium is clearly dynamic: populations themselves are not stable, but the probability of transitions between states 4, 5, and 6 is stable. It is likely that a similar equilibrium would be observed with a mutation rate of 10−6, or even less, if populations were followed over a sufficiently large number of generations. The effect of population size: The results discussed so far were obtained with Ne = 103. To investigate the effects of population size on the evolution of compatibility types, we repeated similar simulations to those presented in Figure 8, A–D (that is, four different combinations of S and Mu values), in smaller populations (Ne = 102) and larger populations (Ne = 104).The results obtained with Ne = 102 are not presented graphically, as only minor quantitative effects were observed. Overall, reducing the population size from 103 to 102 has very little effect on the evolution of Wolbachia compatibility types. The conclusions drawn with Ne = 103 are retrieved: most importantly, the evolution of compatibility type occurs only when S is sufficiently large to prevent population extinction. More interesting are the simulations performed with Ne = 104 (presented in Figure 8, E–H). Most importantly, we note that when S = 1, increasing population size reduces risks of population extinction. In Figure 8, E and G, state 2 (population extinction) is observed in only a small fraction of all simulations, contrasting with that observed in Figure 8, A and C. To understand this result, one must recall (i) that population extinction can occur when S = 1 due to neutral fixation of a mutated lock (i.e., with lock = 2 and key = 1) and (ii) that the time separating the occurrence of a neutral allele and its fixation by drift is 4Ne on average (Kimura 1983). In other words, in large populations, the two mutation steps necessary for a new compatibility type to evolve often occur before extinction of the host population. This result is important with regard to the parameter space allowing the evolution of compatibility types: with large Ne, the evolution of compatibility types can occur even if S = 1. However, one should not conclude that large populations are protected from extinction in the long term: extinction being an absorbing state, it remains the ultimate fate of all populations if S = 1. The effect of imperfect maternal transmission and fitness costs: The evolutionary forces acting on maternal transmission rates and fitness effects to the host have been thoroughly worked out by analytical models (Turelli 1994): increased transmission rates and decreased cost to the host are always selected for. In other words, whatever the initial situation in a population, TE and FE quickly reach their maximum possible value (TEmax and FEmax). In all simulations presented so far, we assumed that maternal transmission can be perfect and that infection is not necessarily costly to the host (TEmax = 1 and FEmax = 1), so that Wolbachia variants with TE < 1 or FE < 1 were observed only at very low frequency (mutation/selection balance). In this section, we assume that transmission efficiency cannot exceed 90% and that infection reduces the host fitness by at least 10% (TEmax = FEmax = 0.9). Similarly, we assume that the intensity of the mod function cannot exceed 90% (MImax = 0.9). Such upper limits allow us to investigate the consequences of uninfected individuals persisting in the long term. With these upper limits for TE, FE and MI, we repeated simulations similar to those presented in Figure 8: we followed the evolution of populations under two S values (S = 1 or 10), two mutation rates (Mu = 10−6 or 10−5), and two different population sizes (Ne = 103 or 104). The results are presented in Figure 9.
To simplify the analysis, note first that the right columns of Figures 8 and 9 (where S = 10) are strikingly similar. In other words, setting upper limits for MI, TE, and FE does not affect the evolution of compatibility types when S is sufficiently large. We therefore focus on comparing the left columns of Figures 8 and 9 (where S = 1). Let us consider Figure 9A (Mu = 10−6, S = 1, Ne = 103, MImax = TEmax = FEmax = 0.9) and compare it with Figure 8A (same values for Mu, S, and Ne, but MImax = TEmax = FEmax = 1). In these two figures, the initial population state (state 3: predominance of the initial type) decreases at the same rate. In both cases, after 100,000 generations, predominance of the initial type is still observed in 90% of the simulations. In Figure 8A, state 3 is replaced only by state 2 (population extinction), contrasting with Figure 9A, where state 3 is replaced by state 1 (infection loss). The interpretation is the following: in Figure 8A, suicidal variants (e.g., with lock = 2 and key = 1) can get fixed by drift, which causes population extinction; conversely, in Figure 9A (where MImax = TEmax = FEmax = 0.9), the unstable equilibrium frequency (below which infection is deterministically lost) is increased by the presence of suicidal Wolbachia variants, so that the infection is always lost before the host population goes extinct. Thus, our model suggests that setting upper limits for MI, TE, and FE changes population extinction to infection loss. A similar conclusion is drawn from the comparison between Figures 8C and 9C (with Mu = 10−5): in Figure 8C, state 3 is mainly replaced by state 2 (population extinction), while it is replaced by state 1 (infection loss) in Figure 9C. Shifting to Figure 9, E and G (compared with Figure 8, E and G), we can assess the effect of MImax, TEmax, and FEmax in large populations (Ne = 104). Most importantly, we note that risks of infection loss are only partially suppressed in large populations. In others words, even large populations are prone to infection loss when S = 1. | |||||||||
SUMMARY AND CONCLUDING REMARKS We have used here a simulation approach to investigate the evolution of compatibility types under various conditions. Being aware that the flurry of conditions and results is difficult to keep track of, we summarize here our main findings, considering in turn the effect of each parameter (listed in Table 1). n is the number of possible states at each site of the lock and key sequences; it conditions the potential diversity of compatibility types. We found that the value of n has only minor effects on the evolution of compatibility types. We note, however, that increasing the value of n tends to slow down the evolution of new compatibility types, as it reduces the probability of favorable mutations in the key sequence. S is the length of the lock and key sequences; it conditions the potential for partial compatibility. If S is small (i.e., S = 1), any mutation in the lock or key sequence leads to complete incompatibility. On the contrary, if S is large (i.e., S = 10), mutation can give rise to different, but partially compatible lock and key. We found that the value of S has profound implications on the evolution of compatibility types: under small S values, population extinction or infection loss is likely to occur before new compatibility types get fixed. On the contrary, if S is large, populations are not prone to extinction or infection loss, allowing the evolution of compatibility types. Mu is the mutation rate used in the model. We found that higher mutation rates accelerate the evolution of compatibility types, as expected. We also observed that high mutation rates increase risks of population extinction or infection loss when S is small. Ne is the population size; it conditions the intensity of genetic drift and the mean time from mutation to fixation of neutral alleles. Most importantly, we observed that increasing population size tends to reduce risks of population extinction and, to a lesser extent, of infection loss. However, it must be kept in mind that population extinction and infection loss remain the only ultimate fates of all populations when S is small, since these are stable and absorbing states. MImax, TEmax, and FEmax are the upper limits of MI, TE, and FE, respectively; these condition the potential for uninfected individuals to persist in infected populations. We observed that setting these limits below their absolute maximum (that is, <1) prevents populations from going extinct when S is small. However, this does not facilitate the evolution of compatibility types, since risks of population extinction are replaced by risks of infection loss. Overall, we conclude that MImax, TEmax, and FEmax do not affect the conditions under which the evolution of compatibility types is plausible. In summary, our analysis points to S as the major parameter: the evolution of compatibility types most readily occurs when the lock and key sequences are sufficiently long, that is, when the lock/key pairs change gradually rather than suddenly. Interestingly, partial compatibility between distinct Wolbachia strains has been observed in Drosophila: the wRi variant from Drosophila simulans is partially compatible with the wMel variant from D. melanogaster (Poinsot et al. 1998). Similarly, the wCer2 variant from Rhagoletis cerasi is partially compatible with wMel and wRi (Charlat et al. 2004; Riegler et al. 2004). Such empirical data support a model with S > 1. In an alternative model, Dobson (2004) assumed that novel Wolbachia variants arising by mutation can coexist within a host individual together with the original infection. This analysis revealed interesting alternative routes for the evolution of compatibility types. In particular, it suggested that new types might evolve from an initial change in the resc function (mutation in the key sequence). One should note, however, that this remains plausible only if double infection is maintained over a sufficient number of generations. In other words, segregation from double to single infection (due to imperfect vertical transmission) has the potential to impede this particular process. Incorporating these aspects in the present model would allow one to assess the likelihood of this evolutionary trajectory. In this article, we explored the evolution of compatibility under a wide variety of conditions. This analysis relies on heavy assumptions regarding the molecular mechanism of cytoplasmic incompatibility. Progress in this domain will accelerate now that the Wolbachia genome has been fully sequenced and analyzed (Wu et al. 2004). In turn, this should enrich the models and shed light on the evolutionary processes underlying the diversity of Wolbachia compatibility types. | |||||||||
Acknowledgments We are deeply grateful to Thomas Pornin and Matthew Collette for their contribution to debugging. We also thank Frank Jiggins and two anonymous referees for critical reading and constructive comments. | |||||||||
References
| |||||||||