Download supplementary materials

1
FEDER ET AL. ONLINE SUPPLEMENT
2
ADAPTIVE CHROMOSOMAL DIVERGENCE DRIVEN BY MIXED GEOGRAPHIC
3
MODE OF EVOLUTION
4
5
6
ANALYTICAL APPROXIMATION
In addition to computer simulations, we also examined an analytical approach to estimate
7
the probabilities of establishment for new inversions suggested to us by S. Yeaman (pers.
8
comm.). The approach involved splicing the results of Kirkpatrick and Barton (2006; Eq. 3) for
9
approximating the rate of increase in frequency of an inversion,
10
11
λ = 1+[2r/((2r + n-1)ms)](n-1)m
12
13
into Kimura’s diffusion equation for the probability of fixation of a new mutation (see Crow and
14
Kimura 1970; Eq. 8.8.3.13) to estimate the probability of establishment of an inversion (Pr[fix]),
15
as done with single locus models by Yeaman and Otto (2011). Substituting Kirkpatrick and
16
Barton’s (λ – 1) for Kimura’s s yields:
17
18
Pr[ fix] 
1  e 4 N(  1)p
1  e 4 N (  1)
19
20  where, m is the migration rate between populations, s is the selection coefficient for each
21
individual locus, n is the number of loci, N is the population size, r the recombination rate
22
between adjacent loci, and p the starting frequency of the inversion (k/2N).
23
Yeaman and Otto (2011) found that this splicing approach provided an accurate estimate
24
of the establishment probability of a single mutation. We found, however, that this approach did
25
not perform especially well when applied to predict the probability of establishment of an
26
inversion. In general, the estimated probability of establishment for the analytical formula was
27
on the scale of an order of magnitude (or more) higher than that estimated for the mixed mode
28
simulations (see Figs. S1 and S2 for results for an inversion containing four loci under divergent
29
selection). This was true whether or not we included a deleterious meiotic effect for single
30
recombination events in inversion heterokaryotypes (data not shown). In addition, regardless of
31
the rate of gene flow, the analytical approximation was not sensitive to differences in the level of
32
divergent selection affecting loci (see dashed lines in Figs. S1 A-C). In contrast, the interplay
33
between migration and selection strongly influenced establishment probabilities in the mixed
34
mode simulations, as would be intuitively expected (Figs. S1 A-C). When selection is high
35
relative to migration rate, reduced recombination is not as strongly favored (locally favored
36
alleles are at high frequency for all loci), and the probability of establishment for an inversion
37
drops (see right hand side of solid line curves for mixed mode simulations in Figs. S1 A-C). The
38
same is true when selection is weak relative to migration, but in this case the cause is that gene
39
flow tends to swamp local adaptation (see left hand side of solid line curves for mixed mode
40
simulations in Figs. S1 B, C). The analytical and mixed mode simulations results did appear to
41
converge, however, with weak selection acting on loci (s = 0.001) and low migration rate (m =
42
0.001; Fig. S1 A). Moreover, the probabilities of inversion establishment derived from the
43
analytical approximation equation and mixed mode simulations were similarly affected by
44
recombination rate (Fig. S2) under conditions of relatively high migrations rate (m = 0.1) and
45
moderate selection (s = 0.1). However, despite the curves for the analytical approximations and
46
the simulations being similar in shape, the magnitudes of the difference in the probability of
47
establishment for an inversion were still over an order of magnitude higher for the analytical
48
approximation (Fig. S2).
49
We suspect that two factors compromise the analytical approximations. First, the
50
conditions most relevant to the mixed mode model, where migration rates (m) and selection
51
coefficients (s) are high and s is not >> m, negate a number of simplifying analytical assumptions
52
of Kirkpatrick and Barton’s estimate of λ and of this application of Kimura’s diffusion
53
approximation. Second, the selective advantage of reduced recombination afforded by an
54
inversion changes through time in relation to the genetic composition of populations 1 and 2 for
55
locally favored alleles. This is particularly true for the mixed model, where populations are not at
56
selection-migration equilibrium for locally favored alleles at the time of secondary contact and
57
the rate and degree to which introgressed genes accumulate between populations are prime
58
factors influencing the changing selective advantage of the inversion. As the analytical approach
59
assumes a fixed s value, this could have compromised its effectiveness. In contrast, our
60
simulations allowed for fluctuating selection favoring the inversion. Future analytical attempts to
61
estimate the establishment probabilities for new inversions might concentrate on branching
62
approximations in which probabilities of transitions are dependent on the states of populations.
63
64
65
PRESTANDING INVERSIONS IN BOTH POPULATIONS
In the simulations reported in the main text, we considered standing inversion
66
polymorphism to be present at low frequency (k = 1 to 200) in only one of the two populations.
67
Prior to secondary contact, however, it is possible that both populations 1 and 2 contain standing
68
inversion variation for a given genomic region. This will generally increase the probability of
69
establishment of an inversion polymorphism (usually by a factor for k approximately equal to the
70
sum of the number of inversion copies in the two populations combined prior to contact) (Fig.
71
S3). But standing variation in the two populations can complicate the dynamics of the process, as
72
it will usually require one or the other of the inversions (usually the one at lower initial
73
frequency) to be essentially lost while the other is retained. The issues of standing and partially
74
overlapping inversions in both populations prior to secondary contact are topics warranting
75
further investigation.
76
77
78
DELETERIOUS MEIOTIC EFFECTS
The simulation runs described in the main text considered heterokaryotypes to have a
79
selective disadvantage of 10-5 due to meiotic irregularities associated with single exchange
80
events. Varying the level of this selective disadvantage from 10-3 to 10-7 did not greatly affect our
81
results, especially given that divergent selection pressures (s) between populations were over
82
several orders of magnitude higher. However, this does not mean that meiotic problems in
83
heterokaryotypes are unimportant for the dynamics of chromosomal evolution, even when they
84
contribute only slight underdominance to fitness. This is because when migration rates are low
85
relative to divergent selection (m < s), negative frequency dependent selection resulting from
86
meiotic irregularities can still impede the establishment of an inversion polymorphism. Selection
87
favoring reduced recombination is not as strong under this condition, so slight deleterious effects
88
in heterokaryotypes can gain in significance for impeding the establishment of an inversion. In
89
addition, negative frequency dependent selection against the rarer arrangement can also
90
contribute to the fixation of chromosomal arrangements if populations experience a period of
91
allopatry following secondary contact and introgression.
92
93
94
SELECTION ACTING AFTER MATING
Divergent selection acting after mating is generally less effective than selection occurring
95
immediately after migration in maintaining genetic differentiation between populations (Fry
96
2003; Nosil et al. 2005). This is because when selection occurs after mating, migrant alleles are
97
not selected against until after they have a chance to introgress into the gene pool of the alternate
98
population in the form of F1 hybrids. But this is not the case when selection occurs prior to
99
mating. Here, migrant genes are selected against before they occur in F1 hybrids.
100
In the simulations conducted in the main text, we considered divergent selection to occur
101
after migration and before mating. However, we also modeled the consequences of selection
102
acting after mating for its effects on inversion establishment. Selection acting after mating tended
103
to increase the probability of inversion establishment under the mixed mode and sympatric
104
origins models for low (m = 0.001) and modest (m = 0.01) levels of migration (Fig. S4). This
105
was true because the increased rate of effective introgression for low and modest migration rates
106
resulting from selection following mating increased the selective advantage of reduced
107
recombination associated with the inversion. For high migration rate (m = 0.1), there was little
108
effect of when selection occurred on the probabilities of inversion establishment (Fig. S4), as
109
effective introgression was similar between populations whether selection occurred prior to or
110
after mating.
111
GENE FLUX
112
In the stochastic simulations reported in the main text, we did not allow for gene flux
113
between inverted and standard arrangements (i.e., there was no double recombination or gene
114
conversion in heterokaryotypes). In nature, gene flux does occur between inverted and collinear
115
genomic regions. It is not uncommon to observe genetic exchange on the order of 10-6 to 10-9,
116
and sometimes much higher, in heterokaryotypes (Navarro et al. 1997; Schaeffer and Anderson
117
2005). Allowing for gene flux did not greatly affect the probabilities of retention of an inversion
118
in our stochastic simulations (data not shown, spreadsheets of full results available upon
119
request). This is because when an inversion is lost, it generally is lost in the first few generations
120
after it occurs as a new mutation, especially for the sympatric origins model. At this time, the
121
inversion is present at extremely low frequency in rare heterokaryotypes. Thus, under the
122
sympatric origins model, if a new inversion failed to capture all locally favored alleles across
123
loci, it is unlikely to obtain them through gene flux with the standard arrangement before being
124
selectively lost. In contrast, inversions containing all favorable alleles will experience only a
125
very slight drain of positively selected alleles and influx of deleterious alleles due to gene flux
126
during the critical stages of establishment under the sympatric origins and mixed modes models.
127
Gene flux can be an important factor under certain circumstances, however, in facilitating the
128
eventual fixation of alternative arrangements between populations when they have become
129
established. Once an inversion polymorphism is established, gene flux can help sort new, locally
130
favored mutations differentially into inverted vs. standard arrangements when they arise in the
131
wrong genetic background. But gene flux can also impede the evolution of intrinsic postmating
132
isolation caused by negative incompatibilities between universally favored alleles (Navarro and
133
Barton 2003). Future work could examine these possibilities more thoroughly.
134
135
136
POSITIVE EPISTASIS
In the stochastic simulations in the main text, we considered loci to independently and
137
multiplicatively affect fitness. We also examined how positive epistasis fitness interactions
138
between a pair of alleles at two different loci 1 and 2 captured within a rearrangement containing
139
a total of four loci influenced its establishment. We present results in Figures 5 A and B for the
140
sympatric and mixed mode models, respectively, with a recombination rate of r = 0.1 between
141
the four loci and a baseline level of divergent selection of s = 0.1 per locus. Positive epistasis
142
was introduced by adding an epistasis term of e = 0.1 to the net fitness of locus 1 and 2
143
genotypes in population 1 that contained the locally favored alleles a at both loci and subtracted
144
a term e = 0.1 when the locally unfavored allele A was present at both loci. The reverse was true
145
for populations 2, where the epistasis term e = 0.1 was subtracted or added to allele a and allele
146
A containing genotypes at both locus 1 and 2. We then multiplied the fitnesses of locus 1 and 2
147
genotypes by that for loci 3 and 4 to get the overall four locus fitness for each genotype.
148
When positive epistatic fitness interactions between loci were considered, the
149
probabilities of inversion establishment were affected in both models (Figs. S5A, B). For the
150
sympatric origin, positive epistasis can help relax some of the constraint that a new inversion
151
must capture all of the locally adapted alleles across all loci to establish as a polymorphism.
152
However, even so, the inversion must still capture all of the favorable alleles having large effects
153
on fitness in order to establish. Strong positive epitasis can also influence the dynamics of
154
inversion establishment by changing gene frequencies at selection-migration equilibrium under
155
the sympatric origins model. In this case, positive epistasis generally enhances genetic
156
divergence between populations prior to the inversion arising, elevating the probabilities that a
157
new rearrangement will capture locally favored alleles. However, it lessens the subsequent
158
strength of selection favoring reduced recombination, thereby inhibiting the establishment of the
159
new inversion. The latter argument also holds for the mixed mode model. The results showed
160
that positive fitness interactions (e = 0.1) between two of the four loci reduced the probability of
161
establishment of an inversion for low (m= 0.001) to modest (m = 0.01) migration rates. For high
162
migration rate (m = 0.1), positive epistasis had minimal effects under the mixed mode model, but
163
increased the probability of establishment under the sympatric origins model.
164
165
166
NEGATIVE EPISTASIS
To examine the consequences of intrinsic postmating isolation, we performed an analysis
167
of the mixed mode and sympatric models in which one locally adapted a allele at locus 1 that
168
was favored in population 1 negatively interacted with an A allele at a second locus 2 that was
169
favored in population 2. The two negatively interacting loci were considered to reside within a
170
rearrangement that contained a total of four loci, with a recombination rate of r = 0.1 between the
171
four loci, a migration rate m = 0.1 between populations, and a baseline level of divergent
172
selection of s = 0.1 per locus. The fitnesses of all genotypes that contained an a allele at locus 1
173
and an A allele at locus 2 were set equal to 1-e (where e was varied from 0.001 to 0.95)
174
regardless of the other alleles present at locus 1 and locus 2. The two locus fitnesses for locus 1
175
and locus 2 genotypes were then multiplied by that for loci 3 and 4 to get the overall four locus
176
fitness values.
177
Computer simulations of the mixed mode and sympatric origins models indicated that
178
negative epistasis of this type generally made it harder for an inversion polymorphism to become
179
established (Fig. S5C). Indeed, under the sympatric model incompatible alleles were effectively
180
selectively eliminated, making it extremely improbable that a new inversion captured them. As a
181
consequence, there was little or no chance for a new inversion to establish under the sympatric
182
origins model that contained negative epistatically interacting loci. These results suggest that
183
when intrinsic postmating isolation is associated with an inversion it may often evolve after the
184
establishment of the rearrangement in sympatry (Navarro and Barton 2003) or later during a
185
period of secondary allopatry after inversion fixation (Kirkpatrick and Barton 2006). The effects
186
of epistasis warrant further analysis.
187
188
Supplementary Figure Legends
189
Figure S1. Comparisons of the probabilities of establishment between analytical approximation
190
(dashed lines) and computer simulations for the mixed mode model (solid lines) for varying
191
levels of divergent selection (s) acting on loci within the inversion under conditions of relatively
192
(A) low migration rate (m = 0.001), (B) modest migration rate (m = 0.01), (C) and high migration
193
rate (m = 0.1) between equal-sized populations.. Shown on a log scale are the probabilities of
194
establishment for an inversion in population 1 starting from a single copy (k1 = 1) when the
195
rearrangement contained four loci, with a recombination rate of r = 0.1 between loci.
196
197
Figure S2. Comparisons of the probabilities of establishment between analytical approximation
198
(dashed lines) and computer simulations for the mixed mode model (solid lines) for varying
199
levels of recombination (r) between loci within the inversion. Shown on a log scale are the
200
probabilities of establishment for an inversion in population 1 starting from a single copy (k1 =
201
1) when the rearrangement contained four loci, divergent selection of s = 0.1 per locus, and
202
migration rate m = 0.1 between equal-sized populations.
203
204
Figure S3. The effects of prestanding rearrangements present in both populations 1 and 2 on the
205
establishment of an inversion following secondary contact under the mixed mode model. Shown
206
on a log scale are the probabilities of establishment for an inversion in either population 1 or 2
207
estimated from 100,000 stochastic simulation runs with an initial copy number of k1 in
208
population 1 and k2 in population 2. The rearrangement contained four loci, with a
209
recombination rate of r = 0.1 between loci, divergent selection of s = 0.1 per locus, and migration
210
rate m = 0.1 between equal-sized populations.
211
212
Figure S4. Comparison of the effects of selection acting before vs. after mating (sel/mate vs.
213
mate/sel) on the establishment of an inversion under the mixed mode (solid lines) and sympatric
214
origins (stippled lines) models. Shown on a log scale are the probabilities of establishment for
215
an inversion in population 1 derived from 100,000 stochastic simulation runs between equal-
216
sized populations. The rearrangement contained four loci, with a recombination rate of r = 0.1
217
between loci and divergent selection of s = 0.1 per locus.
218
219
Figure S5. The effect of fitness interactions between a pair of loci on the establishment of an
220
inversion under (A) the sympatric origins model with positive epistasis, (B) the mixed mode
221
model with positive epistasis, and (C) the mixed mode model with negative epistasis. Shown on
222
a log scale are the probabilities of establishment for an inversion in population 1 derived from
223
100,000 stochastic simulations with m = 0.001, 0.01 and 0.1 for A) and B) and m = 0.1 for C)
224
between equal-sized populations . The rearrangement contained four loci, with a recombination
225
rate of r = 0.1 between loci, divergent selection of s = 0.1 per locus. For (A) and (B), there was
226
positive epistasis of e = 0.1 between two of the four loci. For C), varying levels of negative
227
epistasis of (e) were considered to act between two of the four loci to reduce fitness to 1-e for the
228
two loci, regardless of habitat. For the sympatric origins model, selection acting against negative
229
epistatic alleles quickly eliminates them from populations, making it extremely improbable that a
230
new inversion will capture them and become established.