Migration and the Genetic Covariance between Habitat Preference and Performance

Evolution repeats itself in replicate long-term studies in the wild

The extent to which evolution is repeatable remains debated. Here, we study changes over time in the frequency of cryptic color-pattern morphs in 10 replicate long-term field studies of a stick insect, each spanning at least a decade (across 30 years of total data). We find predictable "up-and-down" fluctuations in stripe frequency in all populations, representing repeatable evolutionary dynamics based on standing genetic variation. A field experiment demonstrates that these fluctuations involve negative frequency-dependent natural selection (NFDS). These fluctuations rely on demographic and selective variability that pushes populations away from equilibrium, such that they can reliably move back toward it via NFDS. Last, we show that the origin of new cryptic forms is associated with multiple structural genomic variants such that which mutations arise affects evolution at larger temporal scales. Thus, evolution from existing variation is predictable and repeatable, but mutation adds complexity even for traits evolving deterministically under natural selection.

Divergent dynamics of sexual and habitat isolation at the transition between stick insect populations and species

Speciation is often viewed as a continuum along which populations diverge until they become reproductively-isolated species. However, such divergence may be heterogeneous, proceeding in fits and bursts, rather than being uniform and gradual. We show in Timema stick insects that one component of reproductive isolation evolves non-uniformly across this continuum, whereas another does not. Specifically, we use thousands of host-preference and mating trials to study habitat and sexual isolation among 42 pairs of taxa spanning a range of genomic differentiation and divergence time. We find that habitat isolation is uncoupled from genomic differentiation within species, but accumulates linearly with it between species. In contrast, sexual isolation accumulates linearly across the speciation continuum, and thus exhibits similar dynamics to morphological traits not implicated in reproductive isolation. The results show different evolutionary dynamics for different components of reproductive isolation and highlight a special relevance for species status in the process of speciation.

The short-term, genome-wide effects of indirect selection deserve study: A response to Charlesworth and Jensen (2022)

We recently published a paper quantifying the genome-wide consequences of natural selection, including the effects of indirect selection due to the correlation of genetic regions (neutral or selected) with directly selected regions (Gompert et al., 2022). In their critique of our paper, Charlesworth and Jensen (2022) make two main points: (i) indirect selection is equivalent to hitchhiking and thus well documented (i.e., our results are not novel) and (ii) that we do not demonstrate the source of linkage disequilibrium (LD) between SNPs and the Mel-Stripe locus in the Timema cristinae experiment we analyse. As we discuss in detail below, neither of these are substantial criticisms of our work.

Pleiotropy or linkage? Their relative contributions to the genetic correlation of quantitative traits and detection by multitrait GWA studies

Genetic correlations between traits may cause correlated responses to selection. Previous models described the conditions under which genetic correlations are expected to be maintained. Selection, mutation, and migration are all proposed to affect genetic correlations, regardless of whether the underlying genetic architecture consists of pleiotropic or tightly linked loci affecting the traits. Here, we investigate the conditions under which pleiotropy and linkage have different effects on the genetic correlations between traits by explicitly modeling multiple genetic architectures to look at the effects of selection strength, degree of correlational selection, mutation rate, mutational variance, recombination rate, and migration rate. We show that at mutation-selection(-migration) balance, mutation rates differentially affect the equilibrium levels of genetic correlation when architectures are composed of pairs of physically linked loci compared to architectures of pleiotropic loci. Even when there is perfect linkage (no recombination within pairs of linked loci), a lower genetic correlation is maintained than with pleiotropy, with a lower mutation rate leading to a larger decrease. These results imply that the detection of causal loci in multitrait association studies will be affected by the type of underlying architectures, whereby pleiotropic variants are more likely to be underlying multiple detected associations. We also confirm that tighter linkage between nonpleiotropic causal loci maintains higher genetic correlations at the traits and leads to a greater proportion of false positives in association analyses.

Large-scale mutation in the evolution of a gene complex for cryptic coloration

The types of mutations affecting adaptation in the wild are only beginning to be understood. In particular, whether structural changes shape adaptation by suppressing recombination or by creating new mutations is unresolved. Here, we show that multiple linked but recombining loci underlie cryptic color morphs of Timema chumash stick insects. In a related species, these loci are found in a region of suppressed recombination, forming a supergene. However, in seven species of Timema, we found that a megabase-size "supermutation" has deleted color loci in green morphs. Moreover, we found that balancing selection likely contributes more to maintaining this mutation than does introgression. Our results show how suppressed recombination and large-scale mutation can help to package gene complexes into discrete units of diversity such as morphs, ecotypes, or species.

Adaptive zones shape the magnitude of premating reproductive isolation in Timema stick insects

Simpson's fossil-record inspired model of 'adaptive zones' proposes that evolution is dominated by small fluctuations within adaptive zones, occasionally punctuated by larger shifts between zones. This model can help explain why the process of population divergence often results in weak or moderate reproductive isolation (RI), rather than strong RI and distinct species. Applied to the speciation process, the adaptive zones hypothesis makes two inter-related predictions: (i) large shifts between zones are relatively rare, (ii) when large shifts do occur they generate stronger RI than shifts within zones. Here, we use ecological, phylogenetic and behavioural data to test these predictions in Timema stick insects. We show that host use in Timema is dominated by moderate shifts within the systematic divisions of flowering plants and conifers, with only a few extreme shifts between these divisions. However, when extreme shifts occur, they generate greater RI than do more moderate shifts. Our results support the adaptive zones model, and suggest that the net contribution of ecological shifts to diversification is dependent on both their magnitude and frequency. We discuss the generality of our findings in the light of emerging evidence from diverse taxa that the evolution of RI is not always the only factor determining the origin of species diversity. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Index of host habitat preference explored by movement-based simulations and trap captures

Animal species likely have different strengths of host habitat preference (HHP) that might be characterized by a standardized index ranging from 0 (no preference) to 1 (maximum preference). We hypothesized that in some species, HHP may result from individuals dispersing out of the host habitat having a probability of turning back at the boundary, or after entering host habitat by reducing speed or increasing size of turning angles. Computer simulations of individuals moving between various sized patches of host and nonhost habitat were conducted based on these three behaviours hypothesized to affect HHP. In the rebounding model, simulations resulted in equilibria of animal numbers inside and outside of host habitat that depend on sizes of these areas, initial number and the rebounding probability. Curvilinear regression of simulation results suggested an equation that predicted numbers in the host habitat and was solved for rebounding probability. A modified equation that sampled population densities (e.g., insect pheromone trap catches) inside and outside host habitat areas gave the rebounding probability, an index of HHP, without requiring the sizes of the areas. Simulations with traps and moving animals verified that the modified equation could predict the index correctly. The modified equation also estimates an index of HHP from sampled densities due to speed reductions and a combination of this and rebounding. Changes in angular turning size upon entering host habitat, however, did not affect habitat preference. Using pheromone trap captures, we found that the lesser date moth Batrachedra amydraula has a HHP for date Phoenix dactylifera plantations of 0.96. Host habitat preference indexes also were calculated from sampled insect densities reported in the literature. The new index of HHP is useful to characterize habitat patches of many organisms and aid understanding of animal spatial distributions and speciation processes. In addition, the index can be applied in studies of invasive species, trap crops of pest insects and conservation management.

Microgeographic morphological variation across larval wood frog populations associated with environment despite gene flow

Gene flow has historically been thought to constrain local adaptation; yet, recent research suggests that populations can diverge despite exchanging genes. Here I use a common garden experiment to assess the combined effects of gene flow and natural selection on morphological variation of 16 wood frog (Rana sylvatica) populations, a species known to experience divergent selection pressures in open- and closed-canopy ponds across relatively small geographic scales. Wood frog tadpoles from different ponds showed significant morphological variation associated with canopy type with a trade-off between tail length and body depth consistent with previous research. In contrast, neutral genetic differentiation of nine microsatellite loci as measured by Jost's D was not associated with canopy type, indicating no pattern of isolation by environment. Genetic structure analyses indicated some substructure across the 16 ponds (K = 4); however, three out of four assigned clusters included both open- and closed-canopy ponds. Together, these results suggest that morphological divergence among these wood frog populations is occurring despite gene flow and that selection within these environments is strong. Furthermore, morphological variation among ponds differed across two sampling periods during larval development, demonstrating the importance of evaluating phenotypic divergence over multiple time periods and at a time relevant to the processes being studied.

Intra- and intersexual differences in parasite resistance and female fitness tolerance in a polymorphic insect

To understand host-parasite interactions, it is necessary to quantify variation and covariation in defence traits. We quantified parasite resistance and fitness tolerance of a polymorphic damselfly (Ischnura elegans), an insect with three discrete female colour morphs but with monomorphic males. We quantified sex and morph differences in parasite resistance (prevalence and intensity of water mite infections) and morph-specific fitness tolerance in the females in natural populations for over a decade. There was no evidence for higher parasite susceptibility in males as a cost of sexual selection, whereas differences in defence mechanisms between female morphs are consistent with correlational selection operating on combinations of parasite resistance and tolerance. We suggest that tolerance differences between female morphs interact with frequency-dependent sexual conflict, which maintains the polymorphism locally. Host-parasite interactions can therefore shape intra- and intersexual phenotypic divergence and interfere with sexual selection and sexual conflict.

Density-dependent selection closes an eco-evolutionary feedback loop in the stick insect Timema cristinae

Empirical demonstrations of feedbacks between ecology and evolution are rare. Here, we used a field experiment to test the hypothesis that avian predators impose density-dependent selection (DDS) on Timema cristinae stick insects. We transplanted wild-caught T. cristinae to wild bushes at 50 : 50 cryptic : conspicuous morph ratio and manipulated density by transplanting either 24 or 48 individuals. The frequency of the conspicuous morph was reduced by 73% in the low-density treatment, but only by 50% in the high-density treatment, supporting a hypothesis of negative DDS. Coupled with previous studies on T. cristinae, which demonstrate that maladaptive gene flow reduces population density, we support an eco-evolutionary feedback loop in this system. Furthermore, our results support the hypothesis that predator satiation is the mechanism driving DDS. We found no effects of T. cristinae density on the abundance or species richness of other arthropods. Eco-evolutionary feedbacks, driven by processes like DDS, can have implications for adaptive divergence and speciation.

The magnitude of local adaptation under genotype-dependent dispersal

Dispersal moves individuals from patches where their immediate ancestors were successful to sites where their genotypes are untested. As a result, dispersal generally reduces fitness, a phenomenon known as "migration load." The strength of migration load depends on the pattern of dispersal and can be dramatically lessened or reversed when individuals move preferentially toward patches conferring higher fitness. Evolutionary ecologists have long modeled nonrandom dispersal, focusing primarily on its effects on population density over space, the maintenance of genetic variation, and reproductive isolation. Here, we build upon previous work by calculating how the extent of local adaptation and the migration load are affected when individuals differ in their dispersal rate in a genotype-dependent manner that alters their match to their environment. Examining a one-locus, two-patch model, we show that local adaptation occurs through a combination of natural selection and adaptive dispersal. For a substantial portion of parameter space, adaptive dispersal can be the predominant force generating local adaptation. Furthermore, genetic load may be largely averted with adaptive dispersal whenever individuals move before selective deaths occur. Thus, to understand the mechanisms driving local adaptation, biologists must account for the extent and nature of nonrandom, genotype-dependent dispersal, and the potential for adaptation via spatial sorting of genotypes.

Sexual dimorphism dominates divergent host plant use in stick insect trophic morphology

BACKGROUND

Clear examples of ecological speciation exist, often involving divergence in trophic morphology. However, substantial variation also exists in how far the ecological speciation process proceeds, potentially linked to the number of ecological axes, traits, or genes subject to divergent selection. In addition, recent studies highlight how differentiation might occur between the sexes, rather than between populations. We examine variation in trophic morphology in two host-plant ecotypes of walking-stick insects (Timema cristinae), known to have diverged in morphological traits related to crypsis and predator avoidance, and to have reached an intermediate point in the ecological speciation process. Here we test how host plant use, sex, and rearing environment affect variation in trophic morphology in this species using traditional multivariate, novel kernel density based and Bayesian morphometric analyses.

RESULTS

Contrary to expectations, we find limited host-associated divergence in mandible shape. Instead, the main predictor of shape variation is sex, with secondary roles of population of origin and rearing environment.

CONCLUSION

Our results show that trophic morphology does not strongly contribute to host-adapted ecotype divergence in T. cristinae and that traits can respond to complex selection regimes by diverging along different intraspecific lines, thereby impeding progress toward speciation.

Establishment of new mutations under divergence and genome hitchhiking

Theoretical models addressing genome-wide patterns of divergence during speciation are needed to help us understand the evolutionary processes generating empirical patterns. Here, we examine a critical issue concerning speciation-with-gene flow: to what degree does physical linkage (r < 0.5) of new mutations to already diverged genes aid the build-up of genomic islands of differentiation? We used simulation and analytical approaches to partition the probability of establishment for a new divergently selected mutation when the mutation (i) is the first to arise in an undifferentiated genome (the direct effect of selection), (ii) arises unlinked to any selected loci (r = 0.5), but within a genome that has some already diverged genes (the effect of genome-wide reductions in gene flow for facilitating divergence, which we term 'genome hitchhiking'), and (iii) arises in physical linkage to a diverged locus (divergence hitchhiking). We find that the strength of selection acting directly on a new mutation is generally the most important predictor for establishment, with divergence and genomic hitchhiking having smaller effects. We outline the specific conditions under which divergence and genome hitchhiking can aid mutation establishment. The results generate predictions about genome divergence at different points in the speciation process and avenues for further work.

Biogeographical survey identifies consistent alternative physiological optima and a minor role for environmental drivers in maintaining a polymorphism

The contribution of adaptive mechanisms in maintaining genetic polymorphisms is still debated in many systems. To understand the contribution of selective factors in maintaining polymorphism, we investigated large-scale (>1000 km) geographic variation in morph frequencies and fitness-related physiological traits in the damselfly Nehalennia irene. As fitness-related physiological traits, we investigated investment in immune function (phenoloxidase activity), energy storage and fecundity (abdomen protein and lipid content), and flight muscles (thorax protein content). In the first part of the study, our aim was to identify selective agents maintaining the large-scale spatial variation in morph frequencies. Morph frequencies varied considerably among populations, but, in contrast to expectation, in a geographically unstructured way. Furthermore, frequencies co-varied only weakly with the numerous investigated ecological parameters. This suggests that spatial frequency patterns are driven by stochastic processes, or alternatively, are consequence of highly variable and currently unidentified ecological conditions. In line with this, the investigated ecological parameters did not affect the fitness-related physiological traits differently in both morphs. In the second part of the study, we aimed at identifying trade-offs between fitness-related physiological traits that may contribute to the local maintenance of both colour morphs by defining alternative phenotypic optima, and test the spatial consistency of such trade-off patterns. The female morph with higher levels of phenoloxidase activity had a lower thorax protein content, and vice versa, suggesting a trade-off between investments in immune function and in flight muscles. This physiological trade-off was consistent across the geographical scale studied and supports widespread correlational selection, possibly driven by male harassment, favouring alternative trait combinations in both female morphs.

A conceptual framework for the evolution of ecological specialisation

Ecology Letters (2011) 14: 841-851 ABSTRACT: Ecological specialisation concerns all species and underlies many major ecological and evolutionary patterns. Yet its status as a unifying concept is not always appreciated because of its similarity to concepts of the niche, the many levels of biological phenomena to which it applies, and the complexity of the mechanisms influencing it. The evolution of specialisation requires the coupling of constraints on adaptive evolution with covariation of genotype and environmental performance. This covariation itself depends upon organismal properties such as dispersal behaviour and life history and complexity in the environment stemming from factors such as species interactions and spatio-temporal heterogeneity in resources. Here, we develop a view on specialisation that integrates across the range of biological phenomena with the goal of developing a more predictive conceptual framework that specifically accounts for the importance of biotic complexity and coevolutionary events.

Genetic constraints for thermal coadaptation in Drosophila subobscura

Background

Behaviour has been traditionally viewed as a driver of subsequent evolution because behavioural adjustments expose organisms to novel environments, which may result in a correlated evolution on other traits. In Drosophila subobscura, thermal preference and heat tolerance are linked to chromosomal inversion polymorphisms that show parallel latitudinal clines worldwide, such that "cold-climate" ("warm-climate") chromosome arrangements collectively favour a coherent response to colder (warmer) settings as flies carrying them prefer colder (warmer) conditions and have lower (higher) knock out temperatures. Yet, it is not clear whether a genetic correlation between thermal preference and heat tolerance can partially underlie such response.

Results

We have analyzed the genetic basis of thermal preference and heat tolerance using isochromosomal lines in D. subobscura. Chromosome arrangements on the O chromosome were known to have a biometrical effect on thermal preference in a laboratory temperature gradient, and also harbour several genes involved in the heat shock response; in particular, the genes Hsp68 and Hsp70. Our results corroborate that arrangements on chromosome O affect adult thermal preference in a laboratory temperature gradient, with cold-climate O_st carriers displaying a lower thermal preference than their warm-climate O₃₊₄ and O₃₊₄₊₈ counterparts. However, these chromosome arrangements did not have any effect on adult heat tolerance and, hence, we putatively discard a genetic covariance between both traits arising from linkage disequilibrium between genes affecting thermal preference and candidate genes for heat shock resistance. Nonetheless, a possible association of juvenile thermal preference and heat resistance warrants further analysis.

Conclusions

Thermal preference and heat tolerance in the isochromosomal lines of D. subobscura appear to be genetically independent, which might potentially prevent a coherent response of behaviour and physiology (i.e., coadaptation) to thermal selection. If this pattern is general to all chromosomes, then any correlation between thermal preference and heat resistance across latitudinal gradients would likely reflect a pattern of correlated selection rather than genetic correlation.

Ecological niche dimensionality and the evolutionary diversification of stick insects

The degree of phenotypic divergence and reproductive isolation between taxon pairs can vary quantitatively, and often increases as evolutionary divergence proceeds through various stages, from polymorphism to population differentiation, ecotype and race formation, speciation, and post-speciational divergence. Although divergent natural selection promotes divergence, it does not always result in strong differentiation. For example, divergent selection can fail to complete speciation, and distinct species pairs sometimes collapse (‘speciation in reverse’). Widely-discussed explanations for this variability concern genetic architecture, and the geographic arrangement of populations. A less-explored possibility is that the degree of phenotypic and reproductive divergence between taxon pairs is positively related to the number of ecological niche dimensions (i.e., traits) subject to divergent selection. Some data supporting this idea stem from laboratory experimental evolution studies using Drosophila, but tests from nature are lacking. Here we report results from manipulative field experiments in natural populations of herbivorous Timema stick insects that are consistent with this ‘niche dimensionality’ hypothesis. In such insects, divergent selection between host plants might occur for cryptic colouration (camouflage to evade visual predation), physiology (to detoxify plant chemicals), or both of these niche dimensions. We show that divergent selection on the single niche dimension of cryptic colouration can result in ecotype formation and intermediate levels of phenotypic and reproductive divergence between populations feeding on different hosts. However, greater divergence between a species pair involved divergent selection on both niche dimensions. Although further replication of the trends reported here is required, the results suggest that dimensionality of selection may complement genetic and geographic explanations for the degree of diversification in nature.

Heterogeneous genomic differentiation between walking-stick ecotypes: "isolation by adaptation" and multiple roles for divergent selection

Genetic differentiation can be highly variable across the genome. For example, loci under divergent selection and those tightly linked to them may exhibit elevated differentiation compared to neutral regions. These represent "outlier loci" whose differentiation exceeds neutral expectations. Adaptive divergence can also increase genome-wide differentiation by promoting general barriers to neutral gene flow, thereby facilitating genomic divergence via genetic drift. This latter process can yield a positive correlation between adaptive phenotypic divergence and neutral genetic differentiation (described here as "isolation-by-adaptation"). Here, we examine both these processes by combining an AFLP genome scan of two host plant ecotypes of Timema cristinae walking-sticks with existing data on adaptive phenotypic divergence and ecological speciation in these insects. We found that about 8% of loci are outliers in multiple population comparisons. Replicated comparisons between population-pairs using the same versus different host species revealed that 1-2% of loci are subject to host-related selection specifically. Locus-specific analyses revealed that up to 10% of putatively neutral (nonoutlier) AFLP loci exhibit significant isolation-by-adaptation. Our results suggest that selection may affect differentiation directly, via linkage, or by facilitating genetic drift. They thus illustrate the varied and sometimes nonintuitive contributions of selection to heterogeneous genomic differentiation.

NATURAL SELECTION IN POPULATIONS SUBJECT TO A MIGRATION LOAD

Migration tends to oppose the effects of divergent natural selection among populations. Numerous theoretical and empirical studies have demonstrated that this migration-selection balance constrains genetic divergence among populations. In contrast, relatively few studies have examined immigration's effects on fitness and natural selection within recipient populations. By constraining local adaptation, migration can lead to reduced fitness, known as a "migration load," which in turn causes persistent natural selection. We develop a simple two-island model of migration-selection balance that, although very general, also reflects the natural history of Timema cristinae walking-stick insects that inhabit two host plant species that favor different cryptic color patterns. We derive theoretical predictions about how migration rates affect the level of maladaptation within populations (measured as the frequency of less-cryptic color-pattern morphs), which in turn determines the selection differential (the within-generation morph frequency change). Using data on color morph frequencies from 25 natural populations, we confirm previous results showing that maladaptation is higher in populations receiving more immigrants. We then present novel evidence that this increased maladaptation leads to larger selection differentials, consistent with our model. Our results provide comparative evidence that immigration elevates the variance in fitness, which in turn leads to larger selection differentials, consistent with Fisher's Theorem of Natural Selection. However, we also find evidence that recurrent adult migration between parapatric populations may tend to obscure the effects of selection.

Effects of migration on the genetic covariance matrix

In 1996, Schluter showed that the direction of morphological divergence of closely related species is biased toward the line of least genetic resistance, represented by g(max), the leading eigenvector of the matrix of genetic variance-covariance (the G-matrix). G is used to predict the direction of evolutionary change in natural populations. However, this usage requires that G is sufficiently constant over time to have enough predictive significance. Here, we explore the alternative explanation that G can evolve due to gene flow to conform to the direction of divergence between incipient species. We use computer simulations in a mainland-island migration model with stabilizing selection on two quantitative traits. We show that a high level of gene flow from a mainland population is required to significantly affect the orientation of the G-matrix in an island population. The changes caused by the introgression of the mainland alleles into the island population affect all aspects of the shape of G (size, eccentricity, and orientation) and lead to the alignment of g(max) with the line of divergence between the two populations' phenotypic optima. Those changes decrease with increased correlation in mutational effects and with a correlated selection. Our results suggest that high migration rates, such as those often seen at the intraspecific level, will substantially affect the shape and orientation of G, whereas low migration (e.g., at the interspecific level) is unlikely to substantially affect the evolution of G.

Divergent Host Plant Adaptation and Reproductive Isolation between Ecotypes of Timema cristinae Walking Sticks

Theoretical and empirical studies have demonstrated that divergent natural selection can promote the evolution of reproductive isolation. Three unresolved questions concern the types of reproductive barriers involved, the role of geography, and the factors determining the extent of progress toward complete speciation. Here I synthesize studies of Timema cristinae host plant ecotypes to address these issues. The approach is to compare the magnitude of multiple reproductive barriers among different ecological and geographic scenarios, where pairs of populations within each scenario are the unit of replication. Application of this approach to T. cristinae revealed that divergent host adaptation can promote the evolution of diverse reproductive barriers, including those that are not inherently ecological. Gene flow in parapatry tended to constrain divergence, with the notable exception of the reinforcement of sexual isolation. Thus, geography affected progress toward speciation but did not influence all reproductive barriers in the same way. Studies of any single pair of taxa often capture only certain stages of the speciation process. For example, reproductive isolation between T. cristinae ecotypes is incomplete, and so only the stages before the completion of speciation have been examined. Studies of more divergent taxa within the genus are required to determine the factors that complete speciation.

Natural selection and divergence in mate preference during speciation

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 of the evolution of sexual isolation can be enhanced by isolating the roles of diverse ecological and evolutionary processes.