Genetica
DOI 10.1007/s10709-006-0013-6
ORIGINAL PAPER
Natural selection and divergence in mate preference during
speciation
Patrik Nosil Æ Bernard J. Crespi Æ Regine Gries Æ
Gerhard Gries
Received: 10 December 2005 / Accepted: 25 April 2006
Springer Science+Business Media B.V. 2006
Abstract Sexual isolation can evolve due to natural
selection against hybrids (reinforcement). However,
many different forms of hybrid dysfunction, and
selective processes that do not involve hybrids, can
contribute to the evolution of sexual isolation. Here we
review how different selective processes affect the
evolution of sexual isolation, describe approaches for
distinguishing among them, and assess how they contribute to variation in sexual isolation among populations of Timema cristinae stick-insects. Pairs of
allopatric populations of T. cristinae living on different
host-plant species exhibit greater sexual isolation than
those on the same host, indicating that some sexual
isolation has evolved due to host adaptation. Sexual
isolation is strongest in regions where populations on
different hosts are in geographic contact, a pattern of
reproductive character displacement that is indicative
of reinforcement. Ecological costs to hybridization do
occur but traits under ecological selection (predation)
do not co-vary strongly with the probability of
between-population mating such that selection on
ecological traits is not predicted to produce a strong
correlated evolutionary response in mate preference.
Moreover, F1 hybrid egg inviability is lacking and the
factors contributing to reproductive character displacement require further study. Finally, we show that
sexual isolation involves, at least in part, olfactory
communication. Our results illustrate how understanding
P. Nosil (&) Æ B. J. Crespi Æ R. Gries Æ G. Gries
Department of Biological Sciences, Simon Fraser
University, V5A 1S6 Burnaby, BC, Canada
e-mail: pnosila@sfu.ca
of the evolution of sexual isolation can be enhanced by
isolating the roles of diverse ecological and evolutionary processes.
Keywords Speciation Æ Hybrid fitness Æ Mate choice Æ
Reproductive isolation Æ Pheromones
Introduction
Speciation involves the evolution of reproductive isolation between diverging populations. Understanding
speciation thus requires determining which reproductive barriers initially reduced gene flow between populations and the evolutionary forces producing them
(Mayr 1947, 1963; Coyne and Orr 2004). Barriers to
gene exchange can occur before or after mating, and
different forms of reproductive isolation are not necessarily independent. For example, selection against
hybrids (postmating isolation) can drive the evolution
of increased premating isolation (i.e. reinforcement;
Dobzhansky 1937; Servedio and Noor 2003). The
evolution of premating isolation caused by divergent
mating signals and preferences (sexual isolation hereafter) appears to be an important component of speciation in many taxa (Coyne and Orr 2004). Many
selective processes can affect the evolution of sexual
isolation, but their relative contributions are poorly
understood. In this paper, we describe selective processes which can promote the evolution of sexual isolation, present methods for distinguishing among them,
and apply the methods to explain variation in sexual
isolation among populations of walking-stick insects.
Reproductive isolation can evolve simply as a
by-product of populations adapting to different
123
Genetica
ecological environments (Funk 1998; Schluter 2000;
Jiggins et al. 2001, 2004; Podos 2001; Etges 2002;
Richmond and Reeder 2002; Rundle and Nosil 2005;
Funk et al. 2006; Vines and Schluter 2006). Divergent
natural selection acts on ecologically important traits,
resulting in population divergence in ecological traits.
If these traits, or ones genetically correlated with
them, incidentally affect mate choice, then sexual
isolation evolves as a by-product of local adaptation.
This is the classic ‘by-product’ model of allopatric
speciation (Muller 1942), but it applies to any geographical scenario. Notably, by-product speciation
invokes selection on ecological traits, and the resulting selection on mating traits can be either direct or
indirect (Rundle and Nosil 2005). If the loci affecting
selected and mating traits are the same (i.e. due to
pleiotropy; Kirkpatrick and Ryan 1991), then direct
selection occurs. (Schemske and Bradshaw 1999;
Bradshaw and Schemske 2003) If the loci are physically different, then selection acts indirectly on mate
preference alleles through their genetic association
(i.e. linkage disequilibria) with alleles at other loci,
which are under selection (Barton and Turelli 1991;
Kirkpatrick 1996; Kirkpatrick and Barton 1997;
Kirkpatrick and Servedio 1999; Servedio 2001, 2004).
Such indirect selection acting through imperfect
genetic associations is thought to be less effective at
driving speciation than is direct selection (Kirkpatrick
and Ravigne 2002).
Sexual isolation can also evolve due to direct
selection on actual mate preferences, rather than
selection on ecological traits per se (Servedio 2001).
Selection on a preference is direct when the preference
allele affects fitness independent of the genetic background in which it is found (Kirkpatrick and Ryan
1991; Kirkpatrick 1996; Kirkpatrick and Ravigné
2002). Simple examples are where preferences for
detectable signals accrue high fitness or where preference results in greater parental investment from mating partners. If divergent environments differ in their
signal transmission properties, direct and habitat-specific selection may be imposed on sensory systems and
preferences (i.e. ‘sensory drive’—Morton 1975; Ryan
et al. 1990; Endler 1992; Boughman 2002; Slabbekoorn
and Smith 2002; Seddon 2005; Patten et al. 2004; Fuller
et al. 2005 for review). Preference evolution arising
from such habitat-specific and direct selection on
preferences can be thought of as a form of local
adaptation. Thus both habitat-specific selection on
preferences and the by-product model predict greater
sexual isolation between ecologically divergent pairs of
populations than between ecologically similar pairs of
similar age.
123
Selection against hybrids can also result in the
evolution of sexual isolation (‘reinforcement’ Dobzhansky
1937). During reinforcement, selection acts directly
on genes causing low hybrid fitness, and mating
preferences evolve via their genetic association with
genes causing low hybrid fitness (i.e. via indirect
selection). When populations or species hybridize,
both selection against hybrids and the genetic association (disequilibria) between loci may be large,
such that even indirect selection acting thru imperfect genetic associations is relatively strong (Kirkpatrick
and Ravigné 2002). Initial work on reinforcement
focused on intrinsic postmating isolation (hybrid
inviability or sterility) resulting from between-locus
genetic incompatibilities (Bateson 1909; Dobzhansky
1936, 1937; Servedio and Noor 2003; Gavrilets 2004).
However, selection against hybrids can also occur
without intrinsic postmating isolation. Hybrids may
exhibit a poor ecological fit to the niches of parental
species and such ‘extrinsic hybrid inviability’ can
drive reinforcement (Fisher 1930; Hatfield and
Schluter 1999; Kirkpatrick 2001; Orr and Turelli
2001; Rundle and Whitlock 2001; Rundle 2002;
Servedio 2004). Likewise, sexual selection against
hybrids can contribute to reinforcement (Phelan and
Baker 1987; Noor 1997; Vamosi and Schluter 1999;
Naisbit et al. 2001). Finally, costs to hybridization
that do not involve hybrids can affect preference
evolution. For example, hybridization can reduce the
survival or fertility of females themselves, favoring
the evolution of sexual isolation (i.e. these are the
postmating, prezygotic incompatibilities of Howard
and Gregory 1993; Servedio 2001).
These considerations indicate that numerous selective processes can contribute to the evolution of sexual
isolation, and the relative importance of different processes is a major outstanding question in speciation research. We do not exclude a role for processes that were
not mentioned, such as genetic drift, sexual selection, or
sexual conflict (Ryan and Rand 1993; Panhuis et al.
2001; Shaw and Parsons 2002; Gavrilets 2004). Rather,
we focus on local adaptation and selection against
hybridization because they (1) appear common in nature (Coyne and Orr 2004), (2) are the topic of recent yet
generally untested theory, and (3) apply to our study
system, Timema cristinae walking sticks (drift and
environment-independent sexual selection have been
examined, but evidence suggests they do not contribute
strongly; Table 1). Speciation is often a continuous
process, with different ‘stages’ corresponding to the
magnitude of reproductive isolation that has evolved.
Ecology and geography dictate which processes act, and
genetics influences the response to selection (Fig. 1).
Shown is whether an evolutionary process occurs at all, and then if it contributes to the evolution of sexual isolation
Selection occurs, but trait reducing hybrid survival does not covary strongly with preference
Never been measured, but likely given divergent mate preferences/mating traits of parental forms
Possible
Unknown
Yes
Likely
Does not occur, so cannot contribute
Not examined
No
Unknown
No
Unknown
Selection against hybrids
Intrinsic F1 egg inviability
F1 sterility, F2 and backcross
breakdown
Natural selection
Sexual selection
Selection occurs, but male trait reducing female survival does not covary strongly with preference
Not certain if selection occurs
Possible
Maybe
Postmating, prezygotic incompatibilities
Female survival
Yes
Female fecundity
Maybe
Different populations which use the same host exhibit little or no sexual isolation
Unknown
Environment-independent
sexual selection
No (at least not
in isolation)
No (at least not
in isolation)
Genetic drift
Possible,
evidence lacking
Yes
Direct selection/sensory drive
Possible
Allopatric populations using different hosts show greater sexual isolation than those using the
same host
Trait causing postmating incompatibilities does not covary strongly with mate preference, leaving
a potential role for direct selection
Levels of sexual isolation are uncorrelated with divergence in neutral markers
Yes
Local adaptation
Yes
Explanation
Contributes?
Occurs?
Evolutionary process
Table 1 Summary of the evolutionary processes involved in the evolution of sexual isolation among populations of T. cristinae walking-stick insects
Genetica
Fig. 1 A scenario for how speciation may progress from its
earliest to latest stages (magnitude of reproductive isolation is
depicted by the wavy, black line). The scenario shown here
applies to the evolution of sexual isolation, and lists the
potential processes involved. Populations are initially allopatric,
but secondary contact can occur at any point and results in
additional processes which can drive divergence. In the
particular case shown here, secondary contact occurs roughly
midway through speciation—depicted by the dashed vertical
line. When indirect selection occurs, the genetic covariance
between traits affecting mate preference and traits under
selection is important for evolution. See the Introduction for
details and Table 1 for processes that have been most important
for the evolution of sexual isolation among populations of
T. cristinae walking-sticks. RCD refers to reproductive character displacement. The figure is modified from Rundle and Nosil
(2005)
Distinguishing among selective processes involved in
speciation
A first step towards distinguishing among processes is
to compare sexual isolation among populations with
varying degrees of ecological divergence (a proxy for
local adaptation) and geographic potential for reinforcement (Funk et al. 2002). Further steps involve
elucidating the costs to hybridization potentially
driving reinforcement, quantifying the association
between selected traits and mating preferences (which
affects the strength of indirect selection transmitted to
mate preferences), and determining the traits upon
which sexual isolation is based. We describe these
steps in more detail below, and apply them to populations of Timema walking-sticks. Some of our inferences are based upon indirect evidence but we
nonetheless dissect the role of multiple selective
processes in the evolution of sexual isolation in the
context of an explicit body of theory (Table 1 for
summary).
123
Genetica
Timema walking-stick study system
Timema walking-sticks are wingless, phytophagous
insects inhabiting the chaparral of Southwestern North
America (Crespi and Sandoval 2000). Individuals feed
and mate exclusively on the hosts upon which they rest.
We focus on T. cristinae, a species feeding upon two
different host-plant species (Ceanothus spinosus and
Adenostoma fasciculatum). We define a ‘population’ of
walking-sticks as all the individuals collected within a
homogenous patch of a single host species (as in Nosil
et al. 2002, 2003). Thus ‘parapatric’ insect populations
are in contact with an adjacent population using the
alternative host, whereas ‘allopatric’ populations are
separated from all other populations adapted to the
alternative host by distances >50 times the 12 m pergeneration gene flow distance (Sandoval 1993). Sample
sites with both hosts were chosen such that each population had only one adjacent population on the
alternate host. For simplicity, we use the term
‘hybridization’ to refer to interbreeding between populations on different hosts, but do not imply that the
host forms have achieved species status.
In this study we present much new data, but also
re-analyze and re-evaluate some previously published
data in order to synthesize the collective findings. The
data novel to this paper include: (1) analyses of pheromones and behavioral responses to them, (2) genetic
crosses examining hybrid egg inviability, (3) field collections examining the fitness of different mating-pair
types in nature, and (4) rearing experiment testing for a
heritable basis to population divergence in color-pattern. By contrast, all mate preference data stem from
previously published mate-choice experiments, and the
examination of covariance between color-pattern and
mate preference in the current study involves a novel
analysis of these data (Nosil et al. 2002, 2003).
Collectively, the different types of data provide insight
into the divergence of mating preferences during the
early stages of speciation. We first review previous
evidence for the role of host-adaptation and reinforcement in the evolution of sexual isolation. We then
present new results on the costs to hybridization
potentially driving reinforcement and the traits
affecting sexual isolation.
Host-adaptation and the evolution of sexual isolation
If local adaptation drives divergence, then sexual
isolation is predicted to be greater between ecologically divergent pairs of populations than between
ecologically similar pairs of similar age. Several studies
123
have provided evidence for such a pattern (Funk 1998;
Rundle et al. 2000; Cruz et al. 2004; McKinnon et al.
2004; Funk et al. 2006; Vines et al. 2006). In T. cristinae, a role for host-adaptation is implicated by the
observation that sexual isolation is stronger between
geographically separated pairs of populations using
different host-plant species than between geographically separated pairs of populations using the same
host (Nosil et al. 2002). Conversely, sexual isolation is
uncorrelated with divergence in mitochondrial and
nuclear DNA (COI and ITS-2, respectively) between
these populations. Thus divergence in host plant use,
rather than neutral differentiation via genetic drift, is
the predictor of sexual isolation. The weak sexual
isolation between populations using the same host also
indicates that environment-independent forms of sexual selection do not strongly affect sexual isolation.
Because selection against hybridization cannot occur
when populations are geographically separated, a role
for local adaptation is inferred. We note that sexual
isolation likely has a strong heritable genetic basis
because it is unaffected by rearing environment (i.e.
Ceanothus versus Adenostoma, Nosil et al. 2003).
Reinforcement and reproductive character
displacement
Reinforcement predicts reproductive character
displacement: increased sexual isolation in geographic
regions where hybridization occurs (sympatry/parapatry) relative to regions where it does not (allopatry; but
see Lemmon et al. 2004). Several cases of reproductive
character displacement have been documented
(Wasserman and Koepfer 1977; Zouros and d’Entremont. 1980; Noor 1995; Saetre et al. 1997; Rundle and
Schluter 1998; Higgie et al. 2000; Hobel and Gerhardt
2003; Albert and Schluter 2004; Hoskin et al. 2005), yet
few studies have distinguished among alternative
hypotheses to reinforcement, or the role of different
types of low hybrid fitness (Butlin 1995; Servedio and
Noor 2003 for review). This is not a trivial task because
different processes are not mutually exclusive and
often generate overlapping predictions (Day 2000;
Servedio and Noor 2003). For example, novel signals
might be favored in a new environment to reduce
overlap with signals in the ancestral environment
(Wasserman and Koepfer 1977). In this scenario,
‘competition’ among signals occurs only in parapatry
such that direct selection on the preference (rather
than reinforcing selection) would generate reproductive character displacement (Servedio 2001; Boughman
2002; Slabbekoorn and Smith 2002).
Genetica
In addition to the effects of host adaptation, there is
evidence for reinforcement in T. cristinae. Female
mating discrimination against males from the alternative host is stronger when populations from different
hosts are in geographic contact than when they are
fully allopatric (Nosil et al. 2003). This pattern of
reproductive character displacement could have
evolved in response to maladaptive hybridization
between the host-forms because (1) gene flow and
interbreeding occurs, (2) alternative explanations such
as ecological character displacement, differential
fusion of populations, and population ancestry were
examined, but are unsupported, and (3) indirect evidence suggests that hybrids will suffer low fitness due
to high rates of visual predation. However, costs to
hybridization have not been systematically evaluated
in previous studies. Moreover, the association between
selected traits and mating preferences affects the
potential for reinforcement, but has also not been
examined. In this paper, we present preliminary data
on these two factors, and discuss their potential role in
driving divergence in parapatric populations.
Costs to hybridization potentially contributing
to reinforcement
A clear step beyond simply documenting reproductive
character displacement is to ascertain which forms of
selection drove preference evolution. Multiple forms of
postmating isolation are measured, as well as costs to
hybridization for individuals themselves. If individual
or hybrid fitness is not reduced by a particular mechanism, a critical role for this process is unlikely. Here
we estimate two costs to hybridization, intrinsic hybrid
egg inviability and reduced survival of hybridizing
females due to predation. We also evaluate the
potential for extrinsic reductions in hybrid fitness.
(1) Intrinsic hybrid egg inviability (hatching
success)
We tested for intrinsic F1 hybrid egg inviability using
within-population and between-population crosses
(Table 2 for population pairings; n = 607 crosses).
Nymph T. cristinae were field-collected from 28 populations in the Santa Ynez Mountains, California in
spring 2003 and 2004 using sweep nets. They were
reared in glass jars (20C) with 10–15 individuals per
jar at the University of California at Santa Barbara. All
individuals used in the crosses were sexually immature
instars that were reared to sexual maturity on Ceanothus cuttings. Individuals from different populations,
and males and females, were kept separate during
rearing. Within two days of attaining sexual maturity, a
single virgin male and a single virgin female were
housed together in a Petri dish and observed until they
mated. Then they were fed Ceanothus cuttings every
second day until the female died (females lay eggs
singly).
Hatching success was quantified from egg shell
(n = 30, 958 eggs) characteristics. Hatching nymphs
Table 2 Hatching-success was not reduced for between-population versus within-population crosses
Population pair
Allopatric same hosts
P · PR
PE · WCC
LOG · BT
Parapatric different hosts
HVA · HVC
MA · MC
HA · HC
OUTA · OUTC
R12A · R12C
MBOXC · MBOXA
OGC · OGA
VPA · VPAC
Allopatric different hosts
LA · VPC
R6C · R23A
SC · LRN
Within-population
Between-population
n
F-ratio
p-Value
32 (20)
50 (33)
41 (55)
34 (29)
65 (25)
42 (32)
56
49
11
0.38
2.47
0.11
0.54
0.12
0.74
48
23
41
44
40
45
54
54
54
49
57
53
44
44
46
25
(29)
(30)
(18)
(30)
(35)
(24)
(28)
(34)
72
30
14
51
55
38
28
44
0.09
10.26
0.98
1.90
0.11
0.01
0.03
7.36
0.77
0.004
0.34
0.17
0.74
0.92
0.88
0.01
49 (34)
47 (30)
65 (26)
50
66
43
1.11
0.49
0.40
0.30
0.49
0.53
(34)
(26)
(28)
(32)
(40)
(37)
(42)
(33)
65 (35)
40 (41)
53 (39)
The table shows mean % (s.d.) hatching success for between-population and within-population crosses for 14 different pairs of
T. cristinae walking-stick populations. Also shown are the number of crosses for each population pair and the test statistic for the male
population by female population interaction term (this F-ratio stems from a separate ANOVA analysis for each population pair such
that the degrees of freedom are F1,n where n refers to the total sample size)
123
Genetica
left the egg shell fully intact except for a small opening
at one end. Thus the proportion of hatched eggs within
a brood can be calculated by counting the number of
eggs with or without a small opening at one end.
Hatching success measured this way was highly correlated with hatching success measured by monitoring
broods daily for newly hatched nymphs (r = 0.99,
p < 0.001, n = 12 broods). The repeatability of
hatching success estimates based on egg case characteristics was also estimated by recounting some broods
in the same year and in a subsequent year. Statistical
analyses were conducted on egg number from the year
of initial count, in case some eggs were lost or damaged. This does not affect our results because hatching
success was highly repeatable both within years
(r = 0.97, p < 0.001, n = 28) and between years
(r = 0.95, p < 0.001, n = 107).
In T. cristinae, most viable eggs overwinter in diapause and hatch the year after mating occurred. Thus
our analyses focus on hatching success following one
over-wintering. However, a small portion of eggs diapause over one additional winter. There is no reason to
suspect that variation in the duration of diapause
would differ for within-population versus betweenpopulation crosses, and thus hatching success following
one over-wintering likely provides an unbiased measure of egg viability. However, to ensure that our data
are not affected by variable diapause, we estimated
hatching success for 104 broods after a 2-year (2003–
2005), instead of 1-year (2003–2004), diapause. These
data are also presented.
Two types of ANOVA were used to analyze
hatching success. The first analysis pools data among
different population pairs and examines whether
hatching success varied among the four main types of
crosses (within-population, between populations using
the same host, between parapatric populations using
different hosts, and between allopatric populations
using different hosts). Year (2003 or 2004) was also
included as a factor. Mean hatching success was similar
for all types of crosses (within population
crosses = 49%, s.d. = 31; between populations using
the same host = 40%, s.d. = 28; between parapatric
populations using different hosts = 45%, s.d. = 35;
between allopatric population using different
hosts = 51%, s.d. = 39). Thus there was no evidence
for reduced hatching success of F1 individuals from
between-population crosses (main effects of cross-type
F3,607 = 1.22, p = 0.30, main effects of year
F1,607 = 0.07, p = 0.80, year · cross-type interaction
F3,607 = 1.10, p = 0.35; main effects of cross type in an
analysis excluding the interaction term F3,607 = 0.54,
p = 0.65). Congruent results were obtained when
123
hatching success was calculated after two, rather than
one, over-wintering (within population crosses = 70%,
s.d. = 21, n = 56; between parapatric populations using
different hosts = 72%, s.d. = 23, n = 48; main effects of
cross-type F1,103 = 0.36, p = 0.55).
The analyses above pool the results among population pairs and are general and useful because they yield
very large sample sizes for analysis and do not assume
independence of population pairs. However, analyses
on pooled data might obscure trends in individual
population pairs and do not treat population pairs as
the unit of replication. The second set of analyses
alleviate these problems in two ways. First, we conducted a separate ANOVA for each of the 14 population pairs. The model for each individual population
pair included male population, female population, and
the interaction between male and female population.
We are interested primarily in the interaction term
which tests whether hatching success is dependent on
population of origin for the two sexes (e.g. the
same population or not). Examination of individual
pairs of populations confirmed the absence of F1
hybrid egg inviability (Table 2). The interaction between male population and female population was
significant (p = 0.01) for only one of the 14 population
pairs, and this difference is insignificant following
Bonferroni correction. Second, we examined trends
among population pairs (i.e. using population pairs as
the unit of replication). There were no directional
trends evident such that only four of 14 pairs showed
reduced hatching success in between-population versus
within-population crosses (p = 0.82, Binomial test).
Thus strong, F1 hybrid egg inviability was not detected
in our study, and is therefore unlikely to strongly
contribute to reinforcement in this system.
Intrinsic inviability might be more pronounced in F2
or backcrossed individuals, which were not examined.
Furthermore, hybrid sterility generally evolves prior to
inviability (Coyne and Orr 2004) and was also not
examined. However, even if intrinsic inviability occurs
in advanced hybrids, or if hybrid sterility occurs, it is
likely incomplete because there is evidence for ongoing
gene flow between the host forms (Nosil et al. 2003).
Thus, at the very least, our data indicate that F1 hybrids do not show intrinsic inviability as represented by
hatching success. This is the level of hybrid breakdown
examined in most comparative studies of the evolution
of reproductive barriers (reviewed in Coyne and Orr
2004; see also Mendelson 2003). Whether there has
been sufficient time for hybrid egg inviability to evolve
cannot be answered without studying more divergent
taxa in the genus Timema (Crespi and Sandoval 2000;
Sandoval and Nosil 2005). We note that there has been
Genetica
sufficient time for other forms of reproductive isolation
to evolve in this species and that allopatric pairs of
populations of T. cristinae exhibit substantial divergence in mitochondrial and nuclear DNA, indicative of
reasonably ancient divergence times (Nosil et al. 2002).
(2) Ecological costs to hybridization I—survival of
females
Low survival of females that mate with males from the
alternative host could also select for avoidance of
between-host mating. The survival of females could be
affected by the color-pattern of the males riding on their
back. T. cristinae exhibits two main color-pattern morphs. Both morphs occur on both host species, but relative frequencies have diverged between hosts such that
the unstriped morph is more common on Ceanothus
(mean frequency = 81%), whereas the striped morph
is more common on Adenostoma (mean frequency = 72%). Thus populations on different hosts
show consistent, but not fixed, differences in morph
frequency. Population divergence has occurred via differential visual predation: the unstriped colorpattern confers high survival on Ceanothus but low
survival on Adenostoma, and vice versa for the striped
pattern (Sandoval 1994a, b; Nosil 2004; Nosil et al. 2005).
Population divergence in color-pattern morph
frequencies can result in a cost to between-population
mating, particularly for females. In T. cristinae, males
ride on the back of the female during and following
mating. This period where a male rides on a female’s
back lasts at least several hours in the field (Nosil
unpublished) and several days under laboratory conditions (Nosil et al. 2003). Most individuals within a
population are the locally cryptic color-pattern morph.
Conversely, many immigrant males from the adjacent
population on the alternative host are the locally lesscryptic morph. Thus females that mate with and carry
less-cryptic immigrant males on their back are likely to
suffer higher rates of visual predation than females
who mate with cryptic males from their own population. This reduction in survival of hybridizing females
could favor the evolution of increased mating
discrimination against males from the alternative host.
The scenario outlined above requires reduced
survival of females mating with less-cryptic males.
Studies of individuals (rather than mating pairs)
suggest that this will occur (Sandoval 1994a, b; Nosil
2004), but the fitness of mating pairs themselves has
not been examined. We tested for selection against
less-cryptic mating pair types in natural populations.
If selection against a pair type occurs, its relative frequency should decline through time (i.e. between
sequential sample periods; Endler 1986). In our analyses, the sequential sample periods are successive
months. We have used these procedures previously to
test for selection against the rare, intermediate morph
(Nosil et al. 2003). In 2004, we collected mating pairs at
one site that contained a population on each host
(Refugio 12). We sampled adult insects in April and
May 2004 (which represent periods before and after
selection, respectively), recording the color-morph
combination of each mating pair (total n = 174 mating
pairs). For each host separately, we assessed whether
the frequency of pair-types differed between months
(i.e. before and after selection) using v2-tests. We also
tested whether the change in pair types was dependent
on host species, using a loglinear analysis. This analysis
used a partial v2 value to assess the significance of the
three-way interaction between pair-type, month and
host (Norusis 1993).
We detected evidence for natural selection against
less-cryptic mating pairs in natural populations
(Fig. 2). Thus the relative frequency of the four different male/female pair types changed between months
(i.e. before versus after selection) and, as expected, the
nature of these changes was significantly dependent on
host-plant species (pair-type · month · host interaction, v2 = 8.94, d.f. = 3, p < 0.05, loglinear). On both
hosts, mating pairs where both sexes were the cryptic
morph increased after selection (thus striped/striped
increased on Adenostoma and unstriped/unstriped
increased on Ceanothus). Conversely, pairs where
both sexes were the less-cryptic morph decreased
after selection. However, differences before and
after selection were statistically significant only for
Ceanothus (v2 = 8.78, d.f. = 3, p < 0.05, n = 97; for
Adenostoma, v2 = 5.50, d.f. = 3, p = 0.14, n = 77).
Selection against less-cryptic mating pairs likely
occurs and represents a cost to hybridization potentially
contributing to reproductive character displacement.
Although our method was the only way to assess
selection using undisturbed, natural populations, it
potentially confounds differential survival with differential dispersal (Endler 1986). We note, however, that
the results of manipulative experiments have demonstrated selection against less-cryptic individuals
(Sandoval 1994a; Nosil 2004), increasing the likelihood
that our current results represent selection against lesscryptic mating pairs. The reduced survival of hybridizing females is more likely indirect than direct selection, because the cost to females depends on the traits
(color-pattern in this case) males and females carry
(Servedio 2001, p. 1913 for explicit discussion of this
issue). Thus, genetic covariance between color-pattern
and mate preference is required for mate preference
123
Genetica
Fig. 2 Analyses of natural selection using T. cristinae mating
pairs in the wild show selection against the locally less-cryptic
mating pairs; on both hosts, mating pairs with a cryptic male and
female morph increased after selection (i.e. between successive
sample months). Thus ‘‘striped/striped’’ pairs increased on
Adenostoma and ‘‘unstriped/unstriped’’ pairs increased on
Ceanothus. Conversely, pairs with a less cryptic male and female
morph decreased. UU = both sexes unstriped, US = unstriped
male on striped female, SU = striped male on unstriped female,
SS = both sexes striped
evolution due to direct selection on color-pattern. This
covariance is assessed in a subsequent section.
(3) Ecological costs to hybridization II—survival
of hybrids
Even if hybridizing females survive and produce
offspring, ecologically based natural selection against
hybrids could drive reinforcement. We were unable to
assess directly the fitness of hybrids in the wild.
Instead, we inferred the possibility of ecological
selection against hybrids using estimates of selection
on color-pattern and inferences about the color-pattern
of hybrids. As noted, populations on different hosts
exhibit divergent color-pattern morph frequencies such
that each morph is more common on the host upon
which it has higher survival. If population divergence in
color-pattern has a strong genetic basis, then it is likely
that hybrids will suffer reduced fitness due to exhibiting maladaptive color-patterns. Under additive inheritance of color-pattern, hybrids will exhibit low survival
because intermediate morphs (which bear a faint stripe
and occur at low frequency in the wild) have lower
survival on both hosts than does the locally cryptic
morph (Nosil et al. 2003; Nosil 2004). Under dominance, hybrids will still exhibit lower survival on
average, because hybrid broods have a higher tendency
to be ‘mixed’ such that they exhibit both color-pattern
123
morphs (rather than only the locally cryptic one;
Sandoval 1993; Nosil et al. 2003). It is unclear which
scenario applies because dominance in within-population genetic crosses is incomplete and between-population crosses are lacking. Because of this ambiguity,
we cannot yet provide a quantitative estimate of hybrid
fitness. Rather, we demonstrate that population
divergence in color-pattern has a strong heritable basis
and note that population divergence in color-pattern is
predicted to result in low hybrid fitness under both
modes of inheritance. To the extent that hybrids are
intermediate, or to which hybrid broods contain both
color-pattern morphs, the fitness of hybrid offspring
will be reduced relative to offspring derived from
within-population mating.
We tested for a heritable basis to population divergence in color-pattern using a reciprocal rearing
experiment. For nine populations, we raised field-collected first instar nymphs (n = 1859; 4–6 weeks of
rearing) to sexually mature adults on their native or
alternative host (Fig. 3 for population-specific sample
sizes). We used logistic regression analyses to test
whether color-pattern (proportion of striped individuals) in these populations was affected by genotype,
rearing environment (host plant reared on) or a
genotype by environment interaction (assessing
significance using likelihood ratio tests (LR). We conducted two analyses, one using population of origin as
the genotype term and one using host of origin as the
Fig. 3 Reciprocal rearing experiment demonstrates a heritable
basis to population divergence in color-pattern. Shown is the
proportion of striped adults (±95% C.I) observed when fieldcaught first instar T. cristinae were reared to adults during
4–6 weeks on Ceanothus (C) or Adenostoma (A). Populations
adapted to Adenostoma are depicted on the left, populations
adapted to Ceanothus are on the right. Genotypic effects
(population and host of origin) were large and highly statistically
significant (see Results for statistics). Values below the x-axis
refer to number of individuals
Genetica
genotype term. We report the results from a full model
that included both factors and the interaction between
them as well as the results from a reduced regression
model derived using backward elimination (initial
model included both factors and the interaction term
but then removed all terms for which the significance
of –2 log LR was >0.10).
The reciprocal-rearing experiment revealed that
population divergence in color-pattern has a strong
heritable basis (Fig. 3). A strong effect of genotype
occurred when host of origin was used as the genotype
term (full model, host of origin –2LR = 124.71,
p < 0.001; host reared upon –2LR = 17.05, p < 0.001;
interaction –2LR = 12.82, p < 0.001) and when population of origin was used as the genotype term (full
model, population of origin –2LR = 10.15, p < 0.01;
host reared upon –2LR = 49.12, p < 0.001; interaction
term –2LR = 2.32, p = 0.13; reduced model, population of origin –2LR = 62.91, p < 0.001; host reared
upon –2LR = 56.27, p < 0.001; interaction term
removed). Although rearing effects were also detected,
genotypic effects were consistent and large in magnitude. Heritable population divergence in color-pattern
suggests that hybrid offspring will be less-cryptic such
that selection against them could contribute to reinforcement. Evidence for reduced survival of hybridizing females themselves was presented above. Thus
ecological costs to hybridization likely occur. The
potential for these costs to drive the evolution of
divergent mate preferences depends on the association
between ecologically selected traits and mating
preferences.
The association between selected traits and mate
preferences
Direct selection on a trait generates indirect selection
on correlated traits, and potentially a correlated evolutionary response. Another step in understanding
preference evolution is to examine the magnitude of
the association between traits under direct selection
(e.g. ecological selection) and mate preference.
Evolution of mate preferences by indirect selection
depends on the magnitude of this genetic association. If
the association is weak, then direct selection on mate
preferences themselves may be a more likely cause of
preference evolution (Kirkpatrick and Barton 1997).
In T. cristinae, for direct selection on color-pattern
to cause the evolution of between-population mate
preference, a genetic association must exist between
color-pattern and the probability of between-population mating. We have previously shown that the
probability of between-population mating is independent from whether both sexes are the same or different
morphs (Nosil et al. 2002). These analyses test for
assortative mating by color-pattern, but do not
explicitly test for a genetic association between colorpattern and mate preference per se. A genetic association between these traits would be indicated by a
difference between color-pattern morphs within each
sex in the probability of between-population mating
(i.e. a difference in whether each morph accepts or
rejects individuals from the alternate host). We would
expect that the more-cryptic morph is less-likely to
hybridize in between-population mating trials. These
predictions have not been tested previously, and are
examined here using the mating trials from Nosil et al.
(2002, 2003).
The previous studies examined mate choice using
no-choice mating trials. A single male and a single female were housed in a Petri dish and whether copulation occurred within a 1-hour period was recorded.
We employed logistic regression to test whether the
probability of copulation in these trials was dependent
on color-pattern morph. Our prediction concerns the
probability that different color-pattern morphs within
populations will hybridize with individuals from other
populations. Thus, our analyses are restricted to the
between-population mating trials and we conducted
separate analyses for males and for females within each
of the 12 populations from Nosil et al. (2002, 2003). We
detected little or no evidence for a phenotypic association between color-pattern and mate preference. Thus
color-pattern morphs did not tend to differ in the
probability of between-population copulation when
either males or females were considered (Table 3).
Significant differences between color-pattern morphs
were detected in only two of 24 cases, one was in the
direction opposite to that predicted (i.e. cryptic morph
more likely to hybridize), and neither was significant
following Bonferroni correction (i.e. p = 0.024, 0.032
before correction for 12 comparisons within each sex).
Because the association between color-pattern and
mate preference is weak, selection on color-pattern (in
hybrids or hybridizing parental forms) is not predicted
to cause a strong correlated evolutionary response in
mate preference. Indirect selection on mate preference, generated by direct selection on color, will be
weak. Thus how strongly the observed ecological costs
to hybridization contribute to reproductive character
displacement remains unclear. The observed lack of
covariation between color-pattern and between-population mate preference is surprising because migration
between populations should generate non-random
association between alleles at different loci (i.e. linkage
123
Genetica
Table 3 The proportion of
between-population mating
trials resulting in copulation
for different color-pattern
morphs of T. cristinae. The
probability of copulation did
not tend to differ between
color-pattern morphs. PC,
VPC, PRC and LA are
allopatric populations,
whereas the others are
parapatric. Also shown is the
likelihood ratio from
likelihood ratio tests (all
d.f. = 1)
Population
Unstriped
mean (s.d.)
Striped
mean (s.d.)
Unstriped n
Striped n
LR
p-Value
Males
PC
HVC
VPC
PRC
MBOXC
OGC
HVA
VPA
MA
LA
HA
OUTA
0.32
0.26
0.30
0.29
0.27
0.18
0.49
0.33
0.25
0.23
0.35
0.29
(0.47)
(0.44)
(0.46)
(0.45)
(0.45)
(0.39)
(0.50)
(0.47)
(0.44)
(0.43)
(0.48)
(0.45)
0.29
0.36
0.13
0.00
0.29
0.17
0.35
0.29
0.26
0.26
0.36
0.34
(0.46)
(0.48)
(0.35)
(0.00)
(0.46)
(0.38)
(0.48)
(0.46)
(0.44)
(0.44)
(0.48)
(0.48)
156
54
223
168
26
78
98
148
48
52
48
153
58
59
15
1
174
86
145
66
121
162
192
73
0.15
1.24
2.19
0.67
0.06
0.01
4.61
0.40
0.01
0.17
0.01
0.69
0.70
0.27
0.14
0.41
0.80
0.93
0.03
0.53
0.93
0.68
0.95
0.41
Females
PC
HVC
VPC
PRC
MBOXC
OGC
HVA
VPA
MA
LA
HA
OUTA
0.30
0.19
0.33
0.36
0.15
0.30
0.25
0.33
0.00
0.46
0.34
0.36
(0.46)
(0.40)
(0.47)
(0.48)
(0.36)
(0.46)
(0.44)
(0.47)
(0.00)
(0.50)
(0.48)
(0.48)
0.35
0.22
0.31
0.50
0.19
0.20
0.31
0.30
0.23
0.38
0.32
0.31
(0.48)
(0.42)
(0.47)
(0.71)
(0.40)
(0.40)
(0.46)
(0.46)
(0.42)
(0.49)
(0.47)
(0.47)
263
27
282
225
98
69
92
147
11
69
38
145
54
46
29
2
84
97
173
57
48
128
135
89
0.63
0.11
0.06
0.16
0.45
2.57
0.95
0.15
5.09
1.46
0.08
0.48
0.43
0.74
0.80
0.69
0.50
0.11
0.33
0.70
0.02
0.23
0.78
0.49
disequilibrium), even when they are physically
unlinked (Nei and Li 1973; Kirkpatrick et al. 2002).
The current study examined the phenotypic covariance
between color-pattern and mate preference, whereas
the most relevant parameter for evolutionary response
is genetic covariance. The degree of correspondence
between phenotypic and genetic covariance remains a
topic of debate (Cheverud 1988; Willis et al. 1991;
Arnold 1992; Schluter 1996). If the phenotypic
covariance is not representative of the underlying
genetic covariance, an undetected association at the
genetic level could occur and have driven reinforcement. Moreover, recombination erodes genetic
covariance caused by linkage disequilibrium, potentially accounting for the lack of covariance observed.
Although further work is required, examining the
association between color-pattern and mate preference
was useful because it prevented us from concluding
prematurely that ecological costs to hybridization
(which appear strong) clearly drove reproductive
character displacement.
Given that the available evidence suggests that
indirect selection on mate preference will be weak, a
role for direct selection on mate preference in the
reproductive character displacement of mate preferences becomes highly plausible (Kirkpatrick and
Barton 1997). Further studies in Timema could focus
123
on direct selection on preference and signals, particularly now that candidate signals (pheromones) are
known.
Traits affecting sexual isolation
A different approach to understanding sexual isolation,
which is complementary to studies of geographic variation, is to identify the traits affecting mate choice
(Nagel and Schluter 1998; Jiggins et al. 2001). Determining such traits allows subsequent studies of how
selection acts on them, and thus how selection affects
the evolution of sexual isolation (Hobel and Gerhardt
2003; Boughman et al. 2005). Surprisingly, there are
relatively few systems where both geographic variation
and the key traits affecting sexual isolation have been
studied (reviewed in Coyne and Orr 2004; but see
Higgie et al. 2000; Cruz et al. 2004; McKinnon et al.
2004; Boughman et al. 2005 for some exceptions).
Previous studies did not identify the traits affecting
sexual isolation in T. cristinae. For example, betweenpopulation mating probability is independent of
differences between the sexes in color-pattern, body
size, and body shape (Nosil et al. 2002; Nosil and
Crespi 2004). Here we present evidence from pheromone analyses and behavioral experiments that sexual
Genetica
isolation between populations stems, at least in part,
from divergence in pheromones and responses to them.
Successful copulation in T. cristinae involves at least
two distinct stages. First, the male must approach the
female and attempt to mount her. Second, the female
must allow the male to mount and copulate with her, as
males cannot force copulation. Here we test whether
pheromonal communication affects the initial attraction of males to females, and thus has a potential role
in sexual isolation.
We focus on two population pairs comprised of
three populations. Thus a population adapted to
Adenostoma (HVA) is compared to a population
feeding on the same host plant (OUTA) and to a
population feeding on Ceanothus (PR). The pair of
populations using different hosts exhibits stronger
sexual isolation than the pair using the same host
(Nosil et al. 2002, 2003 for details). For example, using
an index of sexual isolation (IPSI) that ranges from –1
to +1 (with zero and one indicating random mating and
complete assortative mating, respectively; RolanAlavarez and Cabarello 2000), IPSI = 0.20 for
HVA · PR versus IPSI = 0.11 for HVA · OUTA.
Averaged across numerous pairs of populations, mean
IPSI = 0.24 and 0.08 for pairs on different versus the
same host, respectively.
Methods for analyzing pheromones and behavioral
responses
Our examination of pheromonal communication contained two main components; analytical procedures
aimed at pheromone identification and behavioral
experiments. To obtain pheromone, 16–24 male or
female HVA or PR were placed in separate Pyrex glass
chambers (5 · 10 cm) maintained at a photoperiod of
14L:10D and a temperature of 23–25C. A water
aspirator was used to draw humidified, charcoal-filtered
air at a rate of 80 ml/min through the chamber and
through a glass column (6 · 30 mm) filled with Porapak Q (50–80 mesh, Waters Associated Inc. Milford,
MA, USA). After 40 h, absorbed volatiles were
desorbed with 2 ml of redistilled pentane.
Pheromone extract was analyzed by coupled
gas chromatographic-electroantennographic detection
(GC-EAD) and GC–mass spectrometry (MS). Aliquots (3 insect-hour equivalents) of Porapak Q extract
were analyzed by GC-EAD (Arn et al. 1975; Gries
et al. 2002), employing a Hewlett-Packard (HP) 5890
gas chromatograph fitted with a GC column
(30 m · 0.25 or 0.32 mm ID) coated with DB-5,
DB-17, DB-23, or DB-210 (J&W Scientific, Folsom,
California, USA). The use of four GC columns with
different retention characteristics allowed us to calculate the retention indices (RI) (Van den Dool and
Kratz 1963) and intercolumn IR differentials of the
antennal stimulatory compounds. This analytical
method helped reveal the chemical functionality (e.g.,
secondary alcohol) of these compounds, and contributed
to their identification (Nosil et al. unpublished data).
For GC-EAD recordings, an antenna was carefully
removed from an insect’s head, and its base placed into
the opening of a glass capillary electrode (1.0 · 0.58 ·
100 mm) (A-M Systems, Inc., Carlsborg, Washington,
USA) filled with saline solution (Staddon and Everton
1980). The distal end of the antenna with its tip
removed by spring micro-scissors (Fine Science Tools
Inc., North Vancouver, British Columbia, Canada),
was inserted into the recording capillary electrode.
Coupled GC–MS analyses of pheromone extract
employed a Varian Saturn 2000 Ion Trap GC–MS
fitted with the DB-5 or DB-17 column.
GC-EAD analyses of female T. cristinae effluvia
revealed several volatiles that elicited strong responses
from male antennae (Fig. 4). The volatile and antennal
response profiles were very similar for HVA and PR
populations. GC–MS analyses confirmed that the same
components were EAD-active in HVA and PR populations. The mass spectrum and retention indices (RI)
(Van den Dool and Kratz 1963) of the most EADactive component (DB-5: RI = 1258; DB-17:
RI = 1472; DB-23: RI = 1909; DB-210: RI = 1476) in
both populations were strongly indicative of a secondary
alcohol.
We then performed no-choice experiments using
static-air olfactometers (Vet 1983, Takacs and Gries
2001) to assess whether pheromonal cues contribute to
mate discrimination. These experiments provide evidence for olfactory-based discrimination independent
of the molecular structure of the pheromone component(s). The olfactometers’ three chambers (each
10 cm diameter · 3.5 cm high) were linearly interconnected by a glass tube (each 1 cm diameter · 2.5 cm
long). A no-choice design was used because it allowed
a direct comparison to the no-choice experiments used
to quantify sexual isolation (Nosil et al. 2002, 2003),
and it avoids confounding mating propensity with mate
preference (Rolan-Alavarez and Cabarello 2000). Thus
for each trial a male was tested against a single female
that was from either the males own population or from
the alternative population. The experiment tests whether a male exhibits no response (stays in the central
release chamber), responds to the test stimulus (by
approaching the female), or responds to the control
stimulus (approaches the empty chamber).
123
Genetica
FID:
Dectector response (mV)
EAD:
5
FID:
EAD:
5
HVA volatile extract
~ 5mV
antenna
6
7
8
9
10
11
12
13
14
15
11
12
13
14
15
PR volatile extract
~ 1.5mV
antenna
6
7
8
9
10
Retention time (min)
Fig. 4 Flame ionization detector (FID) and electroantennographic detector (EAD: male Timema cristinae antenna)
responses to effluvium volatiles from female HVA (top) and
female PR (bottom). Chromatography: DB-5 GC column
(0.32 mm ID · 30 m); temperature program: 50C (1 min),
10C per min to 200C, then 25C per min to 280C. Most
antennally active compound marked by an arrow
(Experiment 1: 192 trials for HVA·PR; Experiment 2:
126 trials for HVA·OUTA).
We performed three statistical analyses for each
experiment, predicting larger between-population discrimination in the different host experiment (Experiment 1). The first logistic regression analysis uses data
from all the trials for a particular experiment. It tests
whether a male’s approach towards the female (binary
dependent variable—approached occurred or not) was
dependent on the male population, female population
or the interaction between the male and female population. We are interested in the interaction term,
which tests whether the response of a male from a
particular population is dependent on the population
of female. We also included male color morph, male
age upon capture (sexually immature versus mature),
female color morph and female age upon capture as
factors in the analysis. Models not including these
factors yielded similar results and are not shown. Significance was assessed using likelihood-ratio tests from
full and reduced models. The second analysis was the
same as the first, with one modification; we excluded all
the trials where a male remained in the central release
chamber for the entire trial. This exclusion allowed us
to test whether the results of the first analyses are
confirmed when only active males are considered. The
third analysis examined only those trials where a
female was chosen, and tested which factors predict the
time period until the female is chosen. The time until a
female was chosen is continuously-distributed such that
all the terms described in the other two analyses were
included, but analyzed using ANOVA.
Evidence for pheromone-based discrimination
The experiment proceeded in four steps. First,
olfactometers were randomly assigned to locations on
the laboratory bench. Second, females were randomly
assigned to individual olfactometers and to one of the
two side chambers. Females were confined within a
glass tube covered with wire mesh to prevent physical
contact of the bioassayed male with a female as the test
stimulus. The control chamber contained an empty
glass tube of the same type. Third, a male was released
in the central chamber. Fourth, all activity of the male
was recorded for two hours. Specifically, we recorded
when the male moved, and which chamber it entered.
We analyze the first movement of each male because
multiple movements occurred very rarely (5% of trials). Between trials run every other day, olfactometers
were rinsed with 70% ethanol, thoroughly washed with
soap, and left to dry 24 h (i.e. experiments were run
every second day). A total of 318 trials were conducted
123
We detected evidence for olfactory-based behavioral
discrimination for the population pair using different
hosts, but not for the pair using the same host (Table 4
for statistics and Fig. 5 for data). Thus for the
HVA · PR comparison, the probability that a female
was approached was significantly greater for withinpopulation trials than for between-population trials.
This pattern occurred both when all trials were considered, and when only trials where a male entered a
side chamber were considered (p < 0.05 in both cases,
see Fig. 5A–C and Table 4). Additionally, when only
the trials where a female was approached were examined, the time until a female was approached was
significantly lower for within-population versus
between-population trials, as expected if there was
between-population discrimination based upon olfactory cues. In contrast to the results with the pair using
Genetica
Table 4 Statistical analyses showing that responses of male T.
cristinae to olfactory cues are dependent on an interaction
between male population and female population (i.e. male
population · female population interaction term), but only for
the population pair using different hosts (i.e. Experiment #1)
Test
n
Model
Test-statistic
p-Value
Experiment #1—Populations on different hosts
1. Was a female chosen?
192
Full
Reduced
Full
Reduced
Full
Reduced
5.01
7.29
8.14
9.07
4.34
4.77
0.024
0.007
0.004
0.003
0.051
0.019
Full
Reduced
Full
Reduced
Full
Reduced
0.20
Removed
0.55
Removed
1.65
Removed
0.65
>0.10
0.46
>0.10
0.29
>0.10
2. If the male left the center, was a female chosen?
54
3. When a female was chosen, how long did it take?
30
Experiment #2—Populations on the same host
1. Was a female chosen?
126
2. If the male left the center, was a female chosen?
26
3. When a female was chosen, how long did it take?
15
Tests 1 and 2 use logistic regression analyses (test-statistic is -2LR) and test 3 uses ANOVA (test statistic is F-ratio). The ‘full model’
included all the factors examined and the interaction term, and ‘reduced model’ removed all terms for which p > 0.10 (Materials and
methods for details)
different hosts, there was little or no evidence for
olfactory-based discrimination in the trials using the
population pair on the same hosts. For the
HVA · OUTA comparison, significant differences for
within-population versus between-population trials
were not detected in any case (all p > 0.25). In fact,
trends were sometimes in the direction opposite to that
expected if there was between-population discrimination. For example, the time taken for a male to
approach a female tended to be greater when the
female was from the same population (Fig. 5D–F and
Table 4).
The results suggest that sexual isolation between
populations of T. cristinae that use different host plants
involves, at least in part, pheromonal communication.
Behavioral experiments revealed evidence for discrimination based upon pheromonal cues, and several
candidate pheromone components in the insects’
effluvium elicited antennal responses. Subtle, rather
than large, differences in the insects’ volatile effluvia
and antennal response pattern were detected. One
hypothesis is that the host forms of T. cristinae differ in
some qualitative aspect of pheromones, such as enantiomer composition, and that discrimination is based
on enantiospecific responses. The most antennally
active candidate pheromone component indeed can
exist as two enantiomers (Nosil et al. unpublished
data) which may or may not be so for the as yet
unidentified candidate pheromone components.
Enantiospecific pheromone production and differential
attraction to specific blends of enantiomers has been
reported in other insects (Millar et al. 1990, 1991;
Miller et al. 1996; Gries et al. 1999, 2003). Such enantiospecific responses have been shown to result in
assortative mating (Teale et al. 1994) and partial barriers to gene flow between pheromone races (Cognato
et al. 1999). Further work could examine whether the
enantiomer ratio of major pheromone components
differs between populations. We also do not know how
strongly pheromone production in T. cristinae is affected by environmental (i.e. host-plant reared upon)
versus genetic factors. However, sexual isolation itself
is unaffected by rearing environment (Nosil et al.
2003), and pheromonal discrimination in our behavioral experiments was the same for field-caught adults
and for nymphs reared to maturity in the lab in a
common environment (i.e. on Ceanothus). Other
studies of insect pheromones have detected a strong
genetic basis to their composition (Seybold et al. 1995;
Hager and Teale 1996).
We examined male responses to female pheromones, because sexual isolation between allopatric
pairs of populations involves preference divergence of
both sexes, and the male approach to females is an
important component of sexual isolation (Nosil et al.
2002). It is not yet possible to determine conclusively
whether divergence in pheromonal communication is
linked to adaptation to different host plants, because
only two population pairs have been examined. However this level of replication is highly typical of studies
of the traits underlying sexual isolation (Ryan et al.
1990; Nagel and Schluter 1998; Cruz et al. 2001; Jiggins
et al. 2001; Fordyce et al. 2002; Vines and Schluter
2006). Two lines of evidence are highly suggestive.
123
Genetica
Fig. 5 Olfactometer
experiments showing that
male discrimination based on
olfactory cues occurs only for
the T. cristinae population
pair using different hosts.
Panels A–C show the results
from a population pair that
uses different hosts
(HVA · PR), where
discrimination based on
olfactory cues was detected
(all p < 0.05). Panels D–F
show the results from a
population pair that uses the
same host (HVA · OUTA),
where discrimination based
on olfactory cues was not
detected (all p > 0.15)
First, pheromonal discrimination was detected
between the population pair on different hosts, but not
between the pair on the same host. Second, male
mating preferences are not reinforced in parapatry,
such that discrimination by males based upon pheromonal cues is more likely to involve local adaptation
than reinforcement (Nosil et al. 2003).
Sexual isolation due to local adaptation versus
reinforcement need not involve divergence in the same
traits. In T. cristinae, it appears that preferences of
both sexes diverge due to local adaptation, with only
female preferences reinforced in parapatry. Overall
variation in sexual isolation could thus involve multiple
traits, as has been observed in Drosophila. Sexual
isolation between allopatric populations of D. pseud-
123
oobscura and D. persimilis involves divergence in
courtship song (Williams et al. 2001), whereas genomic
regions affecting reinforced mating discrimination
contain genes involved in pheromonal communication
(Ortiz-Barrientos et al. 2004; see also Boughman et al.
2005 for composite mating traits in stickleback). Future
studies of T. cristinae could examine the possibility of
composite mating traits, particularly because courtship
behavior occurs but has not been examined.
Finally, we note that few studies of sexual isolation
have examined mating preferences and also determined the traits (e.g. color-patterns, body size, pheromones, behavior; but see Boughman et al. 2005) upon
which mate preferences are based. In the few cases
of reinforcement where both have been examined,
Genetica
preferences often exhibit greater displacement than
the mating traits themselves, suggesting that selection
acts mostly on receiver selectivity (Hobel and Gerhardt 2003). In cases of local adaptation, interactions
between natural and sexual selection are often involved in the divergence of traits and preferences
(Naisbit et al. 2001; Leal and Fleishman 2004; Boughman et al. 2005). It would be of interest to further
quantify geographic variation in pheromones and other
potential mating traits among populations of T. cristinae, and to compare this variation to geographic variation in preferences. Determining that phermonal
communication is involved in sexual isolation is a
major step because it will allow such studies.
General discussion
Our approach to examining the contribution of different selective processes to the evolution of sexual
isolation may be applied to other systems, and to other
forms of preference (e.g. habitat preferences, temporal
isolation). Many studies have examined the effects of
one or two processes on the evolution of sexual isolation (e.g. herbivorous insects—Funk 1998; Drosophila—Etges 2002; flowering plants—Bradshaw and
Schemske 2003; stickleback fishes—Rundle et al. 2000;
McKinnon et al. 2004; butterflies—Jiggins et al. 2001;
Fordyce et al. 2002; intertidal snails—Rolan-Alvarez
et al. 1997, Cruz et al. 2004). Outstanding questions
concern the relative importance of distinct evolutionary forces such as ecological selection, reinforcing
selection and gene flow, and the traits involved.
Ascertaining these factors will require detailed studies
of individual systems. For example, in the cricket frog
(Acris crepitans) environmental selection for transmission efficiency must be integral for call divergence
because alternative hypotheses such as reinforcement
and pleiotropic effects of body size divergence are ruled
out (Ryan et al. 1990; Ryan and Wilczynski 1991).
In cases of reproductive character displacement,
very few systems have examined multiple costs to
hybridization or the potential role of gene flow
(Servedio and Noor 2003). In T. cristinae, the effects of
reinforcing selection are greatest under intermediate
migration rates, where encounters between populations are common enough to promote reinforcement,
but low enough to prevent gene flow from eroding
divergence (Nosil et al. 2003; see also Noor 1995).
Additionally, although T. cristinae females are selected
to be more discriminating against males from an
adjacent population, this selection has indirectly
resulted in increased mating discrimination against
foreign males from multiple other populations that use
either host (Nosil et al. 2003). Such ‘universal’ effects
of reinforcement may be due to females recognizing
and preferring males from their own population based
on a ‘population-specific’ trait instead of a ‘host-specific’ or ‘species-specific’ trait (Kelly and Noor 1996;
Higgie et al. 2000; see also Zouros and D’Entremont
1980; Hoskin et al. 2005). If such ‘universal’ effects are
common, then reinforcement could contribute to speciation between ecologically similar pairs of populations,
between populations that are geographically separated
from one another, and between conspecific populations. Notably, preference for population-specific traits
makes reinforcement theoretically more likely in an
island-continent scenario (Servedio 2000), perhaps
explaining why reinforcement occurs in T. cristinae in
the face of highly asymmetric gene flow.
Although the current study focused on premating
isolation, it also presents some data on postmating
isolation. We detected no evidence for reduced F1
hatching success in between-population versus withinpopulation crosses. Our data do not allow us to disentangle the effects of differential fertilization (e.g. gametic isolation, Palumbi 1998; Swanson and Vacquier
2002) from differential mortality of embryos. However,
we detected no differences in hatching success, suggesting
that neither process occurs. An alternative explanation,
which is somewhat contrived, is that fertilization success and embryonic death work in opposite directions in
different types of crosses such that they exactly cancel
each other out (yielding equal hatching success in all
cross types). More likely, neither differential fertilization nor differential embryonic mortality occurs. A lack
of intrinsic F1 hybrid inviability between sister species
has certainly been observed in nature (Schluter 2000 for
review; Johannesson et al. 2000; Saldamando et al.
2005), but examples where it does occur are also common (Coyne and Orr 1989, 2004).
Most work on speciation via natural selection
(‘ecological speciation’) has not examined gametic
isolation and intrinsic inviability (but see Lu and Bernatchez 1998), perhaps because these barriers can
evolve via any form of selection or by genetic drift
(Rundle and Nosil 2005 for review). Previous work on
other systems has implicated a role for selection in
both the evolution of gametic isolation and intrinsic
reductions in hybrid fitness. Reproductive proteins involved in fertilization often evolve rapidly via selection
(Vacquier et al. 1997; Swanson and Vacquier 2002).
Likewise, the three genes causing intrinsic postzygotic
isolation which have been identified exhibit a history of
evolution via positive selection (Hmr, Barbash et al.
2004; Nup96, Presgraves et al. 2003; OdsH, Ting et al.
123
Genetica
1998; Wu and Ting 2004 for review; Shuker et al. 2005
for a potential counterexample). However, the causes
of selection (e.g. ecological or not) cannot be determined from these data alone. Gene flow and population subdivision are also predicted to affect the
evolution of genetic incompatibilities (Orr and Orr
1996; Church and Taylor 2002; Gavrilets 2004). The
current study examined scenarios where ecology and
gene flow varied, but the overall lack of inviability
precludes a strong test of the effects of either factor.
Further studies are clearly warranted, because the
processes driving the evolution of genetic incompatibilities remain poorly understood.
Barriers to gene flow that evolve before reproductive isolation is complete can provide particular insight
into speciation, and thus the temporal order of evolution of different forms of isolation is important. A few
comparative studies indicate that sexual isolation can
evolve before intrinsic postmating isolation but the
generality of this finding is unknown (Coyne and Orr
1989, 1997; Mendelson 2003; Ramsey et al. 2003). Every form of reproductive isolation examined previously
in T. cristinae was detected (Table 1). Populations on
different hosts are partially reproductively isolated by
divergent host preferences (Nosil et al. 2006a, b), low
survival of between-host migrants (Nosil 2004), sexual
isolation (Nosil et al. 2002, 2003), and postmating,
prezygotic incompatibilities (Nosil and Crespi 2006).
The current study focused on the same populations
examined in previous work such that in these populations premating isolation has evolved earlier than
intrinsic hybrid egg inviability.
Speciation is often a continuous process whereby
populations diverge from randomly mating units to
reproductively isolated species. Theory clearly demonstrates that the evolution of sexual isolation can involve many selective processes, and our study shows
how their contributions might be examined. The results
show that multiple processes can act simultaneously,
with a central role for ecology and geography. The
host-associated forms of T. cristinae are unlikely to
have achieved species status, as indicated by only
partial barriers to gene flow at the premating level and
weak mtDNA differentiation between adjacent populations on different hosts due to ongoing gene flow
(Nosil et al. 2003). The host forms represent either an
ongoing speciation event or population divergence that
has reached equilibrium such that we have examined
the early stages of the speciation process. Studies of
more divergent taxa are required to test how our
findings apply to the latter stages of speciation (Sandoval and Nosil 2005). The degree to which the same
traits and processes are involved at different stages of
123
speciation is poorly understood (Claridge and Morgan
1993; Ryan and Rand 1993; Mallet et al. 1998; Boake
et al. 1997; Bordenstein et al. 2000; Jiggins et al. 2004).
Studies examining both traits and preferences, and
multiple selective processes at different stages of
divergence, will likely provide the most complete picture of the entire speciation process.
Acknowledgements We thank T.E. Reimchen, M. Kirkpatrick,
A. Mooers, B. Counterman, D. Schluter, C.P. Sandoval,
J. Endler, W. Etges and three anonymous reviewers for discussion and comments on earlier versions of the paper. D. Mclaren,
B. Mickelson, M. Vankoeveringe, T. Luchin, M. Javid, T.
Leversage, H. Brackett, D. Wick, J. Smith, C. Paul, M. Zheng,
and E. Lopez provided field and technical assistance. J. Endler
(UCSB) kindly provided lab space. Financial support was provided by the Natural Sciences and Engineering Research Council
(NSERC–Canada) through a PGS B scholarship to P.N., an
NSERC grant to B.J.C., an NSERC-IRC to G.G. with Phero
Tech Inc., S.C. Johnson Canada, and Global Forest as industrial
sponsors.
References
Albert A, Schluter D (2004) Reproductive character displacement of male stickleback mate preference: reinforcement or
direct selection? Evolution 58:1099–1107
Arn H, Städler E, Rauscher S (1975) The electroantennographic
detector—a selective and sensitive tool in the gas chromatographic analysis of insect pheromones. Z Naturforsch
30c:722–725
Arnold SJ (1992) Constraints on phenotypic evolution. Am Nat
140:S85–S107
Barbash DA, Awadalla P, Tarone AM (2004) Functional
divergence caused by ancient positive selection of a Drosophila hybrid incompatibility locus. PLOS Biol 2:839–848
Barton NH, Turelli M (1991) Natural and sexual selection on
many loci. Genetics 127:229–225
Bateson W (1909) Heredity and variation in modern lights. In:
Seward AC, (ed) Darwin and modern science. Cambridge
University Press, Cambridge, pp 85–101
Boake C, DeAngelis MP, Andreadis DK (1997) Is sexual
selection and species recognition a continuum? Mating
behavior of the stalk-eyed fly Drosophila heteroneura. Proc
Natl Acad Sci USA 94:12442–12445
Bordenstein SR, Drapeau MR, Werren JH (2000) Intraspecific
variation in sexual isolation in the jewel wasp Nasonia.
Evolution 54:567–573
Boughman JW (2002) How sensory drive can promote speciation. Trends Ecol Evol 10:1–7
Boughman JW, Rundle HD, Schluter D (2005) Parallel evolution of sexual isolation in sticklebacks. Evolution 59:361–373
Bradshaw HD, Schemske DW (2003) Allele substitution at a
flower colour locus produces a pollinator shift in monkeyflowers. Nature 426:176–178
Butlin RK (1995) Reinforcement: an idea evolving. Trends Ecol
Evol 10:432–433
Cheverud JM (1988) A comparison of genetic and phenotypic
correlations. Evolution 42:958–968
Church SA, Taylor DR (2002) The evolution of reproductive
isolation in spatially structured populations. Evolution
56:1859–1862
Genetica
Claridge MF, Morgan JC (1993) Geographical variation in
acoustic signals of the planthopper, Nilaparvata bakeri
(Muir), in Asia: species recognition and sexual selection.
Biol J Linn Soc 48:267–281
Cognato AI, Seybold SJ, F.Sperling FAH (1999) Incomplete
barriers to mitochondrial gene flow between pheromone
races of the North American pine engraver, Ips pini
(Say) (Coleoptera : Scolytidae). Proc R Soc Lond B
266:1843–1850
Coyne JA, Orr HA (1989) Patterns of speciation in Drosophila.
Evolution 43:362–381
Coyne JA, Orr HA (1997) Patterns of speciation in Drosophila
revisited. Evolution 51:295–303
Coyne JA, Orr HA (2004) Speciation. Sinauer, Sunderland, MA
Crespi BJ, Sandoval CP (2000) Phylogenetic evidence for the
evolution of ecological specialization in Timema walkingsticks. J Evol Biol 13:249–262
Cruz R, Rolan-Alvarez E, Garcia C (2001) Sexual selection on
phenotypic traits in a hybrid zone of Littorina saxatilis
(Olivi). J Evol Biol 14:773–785
Cruz R, Carballo M, Conde-Padin P, Rolan-Alvarez E, (2004)
Testing alternative models for sexual isolation in natural
populations of Littorina saxatilis: indirect support for
by-product ecological speciation? J Evol Biol 17:288–293
Day T (2000) Sexual selection and the evolution of costly female
preferences: spatial effects. Evolution 54:715–730
Dobzhansky T (1936) Studies on hybrid sterility. II. Localization
of sterility factors in Drosophila pseudoobscura hybrids.
Genetics 121:113–125
Dobzhansky T (1937) Genetics and the origin of species.
Columbia University Press, New York
Endler JA (1986) Natural selection in the wild. Princeton
University Press, Princeton, New Jersey
Endler JA (1992) Signals, signal conditions, and the direction of
evolution. Am Nat 139:125–153
Etges WJ (2002) Divergence in mate choice systems: does
evolution play by the rules? Genetica 116:151–166
Fisher RA (1930) The genetical theory of natural selection.
Claredon press, Oxford
Fordyce JA, Nice CC, Forister ML, Shapiro AM (2002) The
significance of wing pattern diversity in the Lycaenidae:
mate discrimination by two recently diverged species. J Evol
Biol 15:871–879
Fuller RC, Houle D, Travis J (2005) Sensory bias as an explanation for the evolution of mate preferences. Am Nat
166:437–446
Funk DJ (1998) Isolating a role for natural selection in speciation: host adaptation and sexual isolation in Neochlamisus
bebbianae leaf beetles. Evolution 52:1744–1759
Funk DJ, Filchak KE, Feder JL (2002) Herbivorous insects:
model systems for the comparative study of speciation
ecology. Genetica 116:251–267
Funk DJ, Nosil P, Etges W (2006) Ecological divergence
exhibits consistently positive associations with reproductive isolation across disparate taxa. Proc Natl Acad Sci
USA 103:3209–3213
Gavrilets S (2004) Fitness landscapes and the origin of species.
Princeton Univ. Press, Princeton, NJ
Gries G, Gries R, Schaefer PW, Gotoh T, Higashiura Y (1999)
Sex pheromone components of pink gypsy moth, Lymantria
mathura. Naturwissenschaften 86:235–238
Gries R, Khaskin G, Gries G, Bennett RG, King GG, Morewood P, Slessor KN, Moorewood WD (2002) (Z,Z)-4,
7-Tridecadien-(S)-2-yl acetate: sex pheromone of Douglasfir cone gall midge, Contarinia oregonensis. J Chem Ecol
28:2283–2297
Gries R, Khaskin G, Khaskin E, Foltz K, Schaefer PW, Gries G
(2003) Enantiomers of (Z,Z)-6,9-heneicosadien-11-ol: sex
pheromone components of Orgyia detrita. J Chem Ecol
29:2201–2212
Hager BJ, Teale SA 1996. The genetic control of pheromone
production and response in the pine engraver beetle Ips
pini. Heredity 77:100–107
Hatfield T, Schluter D (1999) Ecological speciation in sticklebacks: environment dependent hybrid fitness. Evolution
53:866–873
Higgie M, Chenoweth S, Blows MW (2000) Natural selection
and the reinforcement of mate recognition. Science
290:519–521
Hobel G, Gerhardt HC (2003) Reproductive character displacement in the acoustic communication system of green
tree frogs (Hyla cinerea). Evolution 57:894–904
Hoskin CJ, Higgie M, McDonald KR, Moritz C (2005) Reinforcement drives rapid allopatric speciation. Nature
437:1353–1356
Howard DJ, Gregory PG (1993) Post-insemination signaling
systems and reinforcement. Phil Trans R Soc Lond B
340:231–236
Jiggins CD, Naisbit RE, Coe RL, Mallet J (2001) Reproductive
isolation caused by colour pattern mimicry. Nature 411:302–
305
Jiggins CD, Estrada C, Rodrigues A (2004) Mimicry and the
evolution of premating isolation in Heliconius melpomene
Linnaeus. J Evol Biol 17:680–691
Johannesson K, Larsson A, Cruz R, Garcia C, Rolan-Alvarez E,
(2000) Hybrid fitness seems not to be an explanation for the
partial reproductive isolation between ecotypes of Galician
Littorina saxatilis. J Moll Stud 66:149–156
Kelly JK, Noor MAF (1996) Speciation by reinforcement: a
model derived from studies of Drosophila. Genetics
143:1485–1497
Kirkpatrick M (1996) Good genes and direct selection in evolution of mating preferences. Evolution 50:2125–2140
Kirkpatrick M (2001) Reinforcement during ecological speciation. Proc R Soc London B 268:1259–1263
Kirkpatrick M, Ryan MJ (1991) The evolution of mating preferences and the paradox of the lek. Nature 350:33–38
Kirkpatrick M, Barton NH (1997) The strength of indirect
selection on female mating preferences. Proc Natl Acad Sci
USA 94:1282–1286
Kirkpatrick M, Servedio MR (1999) The reinforcement of
mating preferences on an island. Genetics 151:865–884
Kirkpatrick M, Ravigné V 2002. Speciation by natural and sexual
selection: models and experiments. Am Nat 159:S22-S35
Kirkpatrick M, Johnson T, Barton N (2002) General models of
multilocus evolution. Genetics 161:1727–1750
Leal M, Fleishman LJ (2004) Differences in visual signal design
and detectability between allopatric populations of Anolis
lizards. Am Nat 163:26–39
Lemmon AR, Smadja C, Kirkpatrick M (2004) Reproductive
character displacement is not the only possible outcome of
reinforcement. J Evol Biol 17:177–183
Lu G, Bernatchez L (1998) Experimental evidence for reduced
hybrid viability between dwarf and normal ecotypes of lake
whitefish (Coregonus clupeaformis Mitchell). Proc R Soc
Lond B 265:1025–1030
Mallet J, McMillan WO, Jiggins CD (1998) Mimicry and warning
color at the boundary between races and species In: Howard
DJ, Berlocher S (eds) Endless forms: species and speciation.
Oxford University Press, Oxford, pp 390–403
Mayr E (1947) Ecological factors in speciation. Evolution
1:263–288
123
Genetica
Mayr E (1963) Animal species and evolution. Harvard University Press, Cambridge
McKinnon JS, Mori S, Blackman BK, Kingsley DL, Jamieson L,
Chou J, Schluter D (2004) Evidence for ecology’s role in
speciation. Nature 429:294–298
Mendelson TC (2003) Sexual isolation evolves faster than hybrid
inviability in a diverse and sexually dimorphic genus of fish
(Percidae: Etheostoma). Evolution 57:317–327
Millar JG, Giblin M, Barton D, Morrison A, Underhill EW
(1990) Synthesis and field testing of enantiomers of 6Z,9Zcis-epoxydienes as sex attractants for geometrid moths.
Interaction of enantiomers and regioisomers J Chem Ecol
16:2317–2339
Millar JG, Giblin M, Barton D, Underhill EW (1991) Chiral
lepidopteran sex attractants: blends of optically active C20
and C21 diene epoxides as sex attractants for geometrid and
noctuid moths (Lepidoptera). Environ Entomol 20:450–457
Miller DR, Borden JH, Slessor KN (1996) Enantiospecific pheromone production and response profiles for populations of
pine engraver, Ips pini (Say) (Coleoptera: Scolytidae), in
British Columbia. J Chem Ecol 22:2157–2172
Morton ES (1975) Ecological sources of selection on avian
sounds. Am Nat 109:17–34
Muller HJ (1942) Isolating mechanisms, evolution and temperature. Biol Symp 6:71–125
Nagel L, Schluter D (1998) Body size, natural selection, and
speciation in sticklebacks. Evolution 52:209–218
Naisbit RE, Jiggins CD, Mallet J (2001) Disruptive sexual
selection against hybrids contributes to speciation between
Helioconius cydno and Helioconius melpomene. Proc R Soc
Lond B 268:1849–1854
Nei N, Li W-H, (1973) Linkage disequilibrium in subdivided
populations. Genetics 75:213–219
Noor MAF, (1995) Speciation driven by natural selection in
Drosophila. Nature 375:674–675
Noor MAF, (1997) Genetics of sexual isolation and courtship
dysfunction in male hybrids of Drosophila pseudoobscura
and Drosphila persimilis. Evolution 51:809–815
Norusis MN (1993) SPSS for Windows. Advanced Statistics.
SPSS Inc., Michigan
Nosil P (2004) Reproductive isolation caused by visual predation
on migrants between divergent environments. Proc R Soc
Lond B 271:1521–1528
Nosil P, Crespi BJ (2004) Does gene flow constrain adaptive
divergence or vice versa? A test using ecomorphology and
sexual isolation in Timema cristinae walking-sticks. Evolution
58:102–112
Nosil P, Crespi BJ, Sandoval CP (2002) Host-plant adaptation
drives the parallel evolution of reproductive isolation. Nature 417:441–443
Nosil P, Crespi BJ, Sandoval CP (2003) Reproductive isolation
driven by the combined effects of ecological adaptation and
reinforcement. Proc R Soc Lond B 270:1911–1918
Nosil P, Vines TH, Funk DJ (2005) Perspective: reproductive
isolation caused by natural selection against immigrants
from divergent habitats. Evolution 59:705–719
Nosil P, Crespi BJ, Sandoval CP, 2006a. The evolution of
host preference in allopatric versus parapatric populations
of T. cristinae walking-sticks. J Evol Biol 19:929–942
Nosil P, Crespi BJ, Sandoval CP, Kirkpatrick M, 2006b. Migration and the genetic covariance between habitat preference
and performance. Am Nat 167:E66–E78
Nosil P, Crespi BJ (2006) Cryptic reproductive isolation driven
by ecological divergence. Proc R Soc Lond B 273:991–997
123
Orr HA, Turelli M (2001) The evolution of postzygotic isolation: accumulating Dobzhansky-Muller incompatibilities.
Evolution 55:1085–1094
Orr HA, Orr LH (1996) Waiting for speciation: the effect of
population subdivision on the time to speciation. Evolution
50:1742–1749
Ortiz-Barrientos D., Counterman BA, Noor MAF (2004) The
genetics of speciation by reinforcement. PLOS Biol
2:2256–2263
Panhuis TM, Butlin R, Zuk M, Tregenza T (2001) Sexual
selection and speciation. Trends Ecol Evol 16:364–371
Patten MA, Rotenberry JT, Zuk M (2004) Habitat selection,
acoustic adaptation, and the evolution of reproductive isolation. Evolution 58:2144–2155
Palumbi SR (1998) Species formation and the evolution of
gamete recognition loci. In: Howard DJ, Berlocher SH (eds)
Endless forms: Species and Speciation. Oxford University
Press, New York, pp 271–278
Phelan PL, Baker TC (1987) Evolution of male pheromones in
moths: reproductive isolation through sexual selection?
Science 235:205–207
Podos J (2001) Correlated evolution of morphology and vocal
signal structure in Darwin’s finches. Nature 409:185–188
Presgraves DC, Balagopalan L, Abmayr SM, Orr HA (2003)
Adaptive evolution drives divergence of a hybrid inviability
gene between two species of Drosophila. Nature 423:715–719
Ramsey J, Bradshaw HD, Schemske DW (2003) Components
of reproductive isolation between the monkeyflowers
Mimulus lewissii and M. cardinalis (Phrymaceae). Evolution
57: 1520–1534
Richmond JQ, Reeder TW (2002) Evidence for parallel
ecological speciation in scincid lizards of the Eumeces skiltonianus species group (Squamata: Scincidae). Evolution
56:1498–1513
Rolan-Alvarez E, Johannesson K, Erlandsson J (1997) The
maintenance of a cline in the marine snail Littorina saxatilis:
the role of home site advantage and hybrid fitness. Evolution 51:1838–1847
Rolan-Alvarez E, Caballero A (2000) Estimating sexual selection and sexual isolation effects from mating frequencies.
Evolution 54:30–36
Rundle HD (2002) A test of ecologically dependent postmating
isolation between sympatric sticklebacks. Evolution 56:322–
329
Rundle HD, Schluter D (1998) Reinforcement of stickleback
mating preferences: sympatry breeds contempt. Evolution
52:200–208
Rundle HD, Nagel L, Boughman JW, Schluter D (2000) Natural
selection and parallel speciation in sympatric sticklebacks.
Science 287:306–308
Rundle HD, Whitlock M (2001) A genetic interpretation of
ecologically dependent isolation. Evolution 55:198–201
Rundle HD, Nosil P (2005) Ecological speciation. Ecol Letts
8:336–352
Ryan MJ, Cocroft RB, Wilczynski W (1990) The role of environmental selection in intraspecific divergence of mate
recognition signals in the cricket frog, Acris crepitans.
Evolution 44:1869–1872
Ryan MJ, Wilczynski W (1991) Evolution of intraspecific
variation in the advertisement call of a cricket frog (Acris
crepitans, Hylidae). Biol J Linn Soc 44:249–271
Ryan MJ, Rand AS (1993) Species recognition and sexual
selection as a unitary problem in animal communication.
Evolution 47:647–657
Genetica
Saetre GP, Moum T, Bures S, Kral M, Adamjan M, Moreno J.
1997. A sexually selected character displacement in flycatchers reinforces premating isolation. Nature 387:589–592
Saldamando CI, Tatsuta H, Butlin RK (2005) Hybrids between
Chorthippus brunneus and C-jacobsi (Orthoptera: Acrididae) do not show endogenous postzygotic isolation. Biol J
Linn Soc 84:195–203
Sandoval CP (1993) Geographic, ecological and behavioral
factors affecting spatial variation in color or morph
frequency in the walking-stick Timema cristinae. PhD
Thesis. University of California, Santa Barbara
Sandoval CP, 1994a. The effects of relative geographic scales of
gene flow and selection on morph frequencies in the walking
stick Timema cristinae. Evolution 48:1866–1879
Sandoval CP, 1994b. Differential visual predation on morphs of
Timema cristinae (Phasmatodeae: Timemidae) and its consequences for host range. Biol J Linn Soc 52:341–356
Sandoval CP, Nosil P (2005) Counteracting selective regimes and
host preference evolution in ecotypes of two species of
walking-sticks. Evolution 59:2405–2413
Schemske DW, Bradshaw HD (1999) Pollinator preference and
the evolution of floral traits in monkeyflowers (Mimulus).
Proc Natl Acad Sci USA 96:11910–11915
Schluter D (1996) Adaptive radiation along genetic lines of least
resistance. Evolution 50:1766–1774
Schluter D (2000) The ecology of adaptive radiation. Oxford
University Press, Oxford
Seddon N (2005) Ecological adaptation and species recognition
drives vocal evolution in neotropical suboscine birds.
Evolution 59:200–215
Servedio MR (2000) Reinforcement and the genetics of nonrandom mating. Evolution 54:21–29
Servedio M (2001) Beyond reinforcement: the evolution of
premating isolation by direct selection on preferences
and postmating, prezygotic incompatibilities. Evolution 55:
1909–1920
Servedio M (2004) The evolution of premating isolation: local
adaptation and natural and sexual selection against hybrids.
Evolution 58:913–924
Servedio M, Noor M (2003) The role of reinforcement in
speciation: theory and data. Annu R Ecol Syst 34:339–364
Seybold SJ, Quilici DR, Tillman JA, Vanderwel D, Wood DL,
Blomquist GJ (1995) De novo biosynthesis of the aggregation
pheromone components ipsenol and ipsdienol by the pine bark
beetles Ips paraconfusus Lanier and Ips pini (Say) (Coleoptera: Scolytidae). Proc Natl Acad Sci USA 92:8393–8397
Shaw KL, Parsons YM (2002) Divergence of mate recognition
behavior and its consequences for genetic architectures of
speciation. Am Nat 159:S61–S75
Shuker DM, Underwood K, King TM, Butlin RK (2005) Patterns
of male sterility in a grasshopper hybrid zone imply accu-
mulation of hybrid incompatibilities without selection. Proc
R Soc Lond B 272:2491–2497
Slabbekoorn H, Smith TB (2002) Habitat-dependent song
divergence in the little greenbul: an analysis of selection
pressures on acoustic signals. Evolution 56:1849–1858
Staddon BW, Everton IJ (1980) Haemolymph of milkweed bug
Oncopeltus fasciatus (Heteroptera: Lygaeidae): inorganic
constituents and amino acids. Comp Biochem Physiol
65:371–374
Swanson WJ, Vacquier VD (2002) The rapid evolution of
reproductive proteins. Nat Rev Genet 3:137–144
Takacs S, Gries G (2001) Communication ecology of webbing
cloth moth: evidence for male produced aggregation signal(s). Can Entomol 133:725–727
Teale SA, Hager BJ, Webster FX (1994) Pheromone-based
assortative mating in a bark beetle. Anim Behav 48:569–578
Ting CT, Tsaur SC, Wu ML, Wu CI (1998) A rapidly evolving
homeobox at the site of a hybrid sterility gene. Science
282:1501–1504
Vacquier VD, Swanson WJ, Lee YH (1997) Positive Darwinian
selection on two homologous fertilization proteins: what is
the selective pressure driving their divergence? J Mol Evol
44:S15–S22
Vamosi SM, Schluter D (1999) Sexual selection against hybrids
between sympatric stickleback species: evidence from a field
experiment. Evolution 53:874–879
Van den Dool H, Kratz PD (1963) A generalization of the
retention index system including linear temperature programmed gas-liquid partition chromatography. J Chromatogr 2:463–471
Vet LEM, (1983) Host-habitat location through olfactory cues by
Leptopilina clavipes (Hartig) (Hym.: Eucoilidae), a parasitoid of fungivorous Drosophila: the influence of conditioning. Netherlands J Zool 33:225–248
Vines TH, Schluter D (2006) Strong assortative mating between
allopatric sticklebacks as a by-product of adaptation to
different environments. Proc R Soc Lond B 273:911–916
Wasserman M, Koepfer HR (1977) Character displacement
for sexual isolation between Drosophila mojavensis and
Drosophila arizonensis. Evolution 31:812–823
Williams MA, Blouin AG, Noor MAF (2001) Courtship songs of
Drosophila pseudoobscura and D. persimilis. II. Genetics of
species differences. Heredity 86:68–77
Willis JH, Coyne JA, Kirkpatrick M (1991) Can one predict
the evolution of quantitative characters without genetics?
Evolution 45:441–444
Wu CI, Ting CT (2004) Genes and speciation. Nat Rev Genet
5:114–122
Zouros E, d’Entremont CJ, (1980) Sexual isolation among populations of Drosophila mojavensis: response to pressure
from a related species. Evolution 34:421–430
123