Research article | Open | Published:
Hydrocarbon divergence and reproductive isolation in Timemastick insects
BMC Evolutionary Biologyvolume 13, Article number: 151 (2013)
Individuals commonly prefer certain trait values over others when choosing their mates. If such preferences diverge between populations, they can generate behavioral reproductive isolation and thereby contribute to speciation. Reproductive isolation in insects often involves chemical communication, and cuticular hydrocarbons, in particular, serve as mate recognition signals in many species. We combined data on female cuticular hydrocarbons, interspecific mating propensity, and phylogenetics to evaluate the role of cuticular hydrocarbons in diversification of Timema walking-sticks.
Hydrocarbon profiles differed substantially among the nine analyzed species, as well as between partially reproductively-isolated T. cristinae populations adapted to different host plants. In no-choice trials, mating was more likely between species with similar than divergent hydrocarbon profiles, even after correcting for genetic divergences. The macroevolution of hydrocarbon profiles, along a Timema species phylogeny, fits best with a punctuated model of phenotypic change concentrated around speciation events, consistent with change driven by selection during the evolution of reproductive isolation.
Altogether, our data indicate that cuticular hydrocarbon profiles vary among Timema species and populations, and that most evolutionary change in hydrocarbon profiles occurs in association with speciation events. Similarities in hydrocarbon profiles between species are correlated with interspecific mating propensities, suggesting a role for cuticular hydrocarbon profiles in mate choice and speciation in the genus Timema.
Individuals in natural populations commonly prefer certain trait values over others when choosing their mates [1, 2]. Traits that are used for mate selection include vocal signals (e.g., [3, 4]), pheromones (e.g., ), behavioral repertoires (e.g., ), and morphological compatibility (e.g., [7, 8]). Theory predicts that divergent preferences for certain trait values can generate behavioral reproductive isolation at both the intra- and interspecific levels [9–11]. Consistent with this prediction, mating preferences have been observed to vary among populations and closely related species in nature (e.g., [1, 12, 13]), and in multiple taxa evidence suggests that the resulting behavioral isolation has been involved in speciation .
Species-specific mating preferences can arise via different mechanisms, and one of the critical components to understanding the process of speciation is determining which factors promote divergence in traits used for mate selection. For example, if mating with individuals from other populations is associated with costs, selection for population-level recognition can drive the evolution of mating signals (e.g., [1, 8, 14]). The resulting behavioral isolation can promote speciation and may lead directly to increased levels of prezygotic isolation, as via reinforcement [11, 15]. Alternatively, differences in mate selection traits between species can arise via neutral processes [16, 17], as a by-product of divergent ecological adaptation [13, 18], or they can evolve in parallel with reproductive isolation, for example as a consequence of sexual selection. Under these scenarios, the same forces of evolution that affect signal properties for intra-specific mate choice would also affect reproductive isolation between species.
The different mechanisms that can cause species-specific mating preferences are expected to generate distinct patterns of evolutionary trait divergence . If mating preferences contribute to the evolution and maintenance of reproductive isolation among species, mating signals are expected to change rapidly during the speciation process. Little change might occur between speciation events, because stabilizing selection for reliable species discrimination should decrease variance in mating signals [20, 21]. Stabilizing selection for species discrimination should also favor large signal differences between species relative to intraspecific variation [22, 23]. By contrast, under sexual selection or drift, mating signals may be expected to change continuously over evolutionary time, without accelerated change during periods of species formation [1, 24]. Thus, distinguishing between gradual vs. speciation-associated change of mating signals in a taxon can provide insights on the selective forces underlying phenotypic change [8, 19, 25].
Here, we investigate whether cuticular hydrocarbons may provide a signal for mate recognition in Timema stick insects, and whether hydrocarbon profiles have diverged gradually between species or if most change is associated with speciation events. Chemical communication among insects is extremely widespread (e.g., [26–28]), and cuticular hydrocarbons, in particular, serve as mate recognition signals in many species (reviewed in ). Analyses of courtship behavior suggested that chemical signals also underlie species recognition and premating isolation in Timema stick insects , a genus that comprises 21 described species (16 of them sexual, the others asexual) of plant-feeding insects, distributed primarily in California . In Timema, sexual isolation has been shown to represent an important reproductive barrier separating ecologically isolated populations within, as well as between, species [30, 32, 33].
We analyzed cuticular hydrocarbon components of nine closely-related sexual species of Timema (T. bartmani, T. boharti, T. californicum, T. chumash, T. cristinae, T. knulli, T. petita, T. podura, and T. poppensis) to characterize differences among species and to test whether interspecific matingsare more likely between species with similar than distinct hydrocarbon profiles. Using the phylogeny of these species, we also evaluated the tempo and mode of hydrocarbon profile evolution to develop insights into the processes underlying divergence between species for this trait. Finally, in order to test for possible links between species recognition, intraspecific mate choice, and the early stages of the speciation process in Timema, we tested for divergence in hydrocarbon components between populations within one of the nine species, T. cristinae. This species was chosen because there is extensive evidence for sexual isolation between populations occurring on Ceanothus versus Adenostoma host plants (e.g., [33–35]) but the proximate mechanisms underlying mate discrimination have remained speculative. This analysis allows us to assess whether the same traits that are used to distinguish conspecific from heterospecific individuals may also be used in mate choice within species.
Cuticular hydrocarbon variation between species and populations
For each of the nine species, we quantified the amount of five cuticular hydrocarbon components (3-, 5-, 7-, 9/11- and 13-methylheptacosane, henceforth 3Me27, 5Me27, 7Me27, 9/11Me27, and 13Me27, respectively) of three to 11 virgin females. Because our approach required analyzing nine species in parallel, we focused on only one sex for simplicity, and we used females because a previous study indicated that mate discrimination, at least in its early stages, appears to mainly depend on males. Specifically, males are more likely to pair with, and court, females of their own than of other species . For eight of the nine species we used females from two geographically distant locations, to account for intraspecific variation when evaluating species differences (Table 1). Except for T. bartmani and T. petita, which are known to occur only on a single host plant (Abies concolor and Ceanothus spp, respectively), the two populations of each species were collected on different hosts (Table 1). For the ninth species, T. boharti, we were able to collect individuals from only one location; across the nine species, we therefore included individuals from 17 locations in total. Sixteen of these 17 locations comprised only the focal species, with no other Timema species present on the same or different host plants. For one of the populations (ED, Table 1), the focal species T. chumash overlapped with a second species, T. podura. Because past work implicates reinforcement of mating preferences between co-occurring populations of T. cristinae , we specifically chose isolated populations to avoid potentially increased levels of discrimination at locations where species meet.
In a first step, we tested whether Timema species differed in the composition of their cuticular hydrocarbons, as expected if these chemicals contribute to sexual isolation. Accordingly, a non-parametric MANOVA revealed significant hydrocarbon profile differences among Timema species (populations nested within species; species effect: F8,8 = 10.4, p <0.0001, populations within species: F8,59 = 3.8, p <0.0001; Figure 1). To simplify the graphical representation of hydrocarbon profile differences between species, we combined the hydrocarbon components using discriminant function analysis, after log-contrast transformation [14, 36]. Log contrasts were calculated by dividing the value for each hydrocarbon by the component 13Me27, and then taking the log of these new variables, resulting in four log-contrast transformed values [log-contrasts for 3Me27, 5Me27, 7Me27 and 9/11Me27, referred to as LC3 (Log Contrast 3), LC5, LC7 and LC9-11, respectively] for every individual. Results using other components as the divisor were qualitatively similar. The two first discriminant functions (DF1, DF2) captured 50.1% and 30.4% of the total variance between species, respectively (loadings: DF1= −2.95 × LC3 + 0.03 × LC5 + 0.01 × LC7 + 5.8 × LC9-11; DF2= −5.60 × LC3 - 0.34 × LC5 + 0.63 × LC7 + 0.10 × LC9-11). DF1 and DF2 also encompass sufficient variance to allow for the visual separation of species (Figure 1). While the between-population differences in hydrocarbon profiles were relatively small compared to the between-species differences (78.0% of the variation between species vs. 7.5% between populations within species), T. chumash stood out with individuals from the two populations being characterized by very distinct hydrocarbon profiles (Figure 1, MANOVA: Wilks’ λ = 0.02, F1,7 = 65.4, p = 0.0007). This difference in T. chumash is consistent with previous data on notable between-population divergence in chromosome numbers , and mtDNA (Crespi, unpubl. data). It is unlikely to represent a profile shift in areas of overlap with other species as the hydrocarbon profiles of T. chumash individuals in areas of overlap with T. podura (population ED) appeared more similar to T. podura than the profile of individuals from the isolated T. chumash population (population HW2; Figure 2). Significant hydrocarbon profile differences of much smaller effect size were also found between the two populations of the species T. cristinae (Wilks’ λ = 0.09, F1,8 = 13.0, p = 0.007) and T. poppensis (Wilks’ λ = 0.18, F1,8 = 5.6, p = 0.04) but not for any of the remaining species (T. bartmani, T. californicum, T. knulli, T. petita, or T. podura, all p > 0.27).
Hydrocarbon variation and between-species mating propensities
Next, we tested whether hydrocarbon profile differences may be involved in interspecific mating decisions by pairing individual males and females and recording whether mating occurred. These no-choice mating trials (> 1400 in total, Table 2) revealed large-scale differences among species combinations in their propensity for interspecific matings, ranging from 0% to almost 95%. Some of this variation appeared to stem from variation among species in male “eagerness” for mating. Males of some species courted any female shortly after encountering her, whereas males from other species were reluctant to court females of even their own species. This pattern is reflected by the significant differences among species in mating proportions measured for intraspecific no-choice trials (χ2= 55.3, df = 8, p < 0.0001, range 0.42 to 1.00) as well as by the fact that the interspecific mating proportions measured for the reciprocal species combinations were only moderately correlated (Pearson's product–moment correlation: rho= 0.45; t = 2.9, df = 33, p = 0.007). We therefore analyzed the interspecific mating propensities in two ways (1) using the average values of the reciprocal combinations for each species pair, and (2) using two separate points for each species pair and using male species as the first explanatory factor in the statistical analyses (based on behavioral evidence for a strong male effect on pairing and mating ).
Both types of approaches are consistent with the idea that cuticular hydrocarbon profiles function as signals for interspecific mate discrimination in Timema. The mean interspecific mating proportions were significantly negatively correlated with hydrocarbon profile differences between species (as measured by species distances in the multivariate hydrocarbon component space; Mantel’s r: -0.43, p= 0.001). Interspecific mating propensities were also negatively correlated with mitochondrial genetic divergences between species (Mantel’s r: -0.53, p= 0.002). The correlation between hydrocarbon profile differences and mating propensity was not, however, explained by genetic divergences, as these two variables were not significantly correlated with each other (Mantel’s r: 0.22, p= 0.09), and hydrocarbon profile differences were still negatively correlated with interspecific mating propensities after correcting for genetic divergences between species (partial Mantel: r = −0.39, p= 0.002).
We obtained qualitatively the same results when using each male–female species combination as a separate data point and using male species as an explanatory variable (to correct for male “eagerness”). Thus, interspecific mating propensity was significantly affected by cuticular hydrocarbon profile divergence, even after correction for genetic divergence (glm with quasibinomial error distribution; effect of genetic divergence: t= −5.5, p< 0.0001, hydrocarbon divergence: t= −2.9, p= 0.006). However, the significance values from these analyses must be considered with caution due to the non-independence of the different pairwise species divergences.
To develop insights into possible contributions of individual hydrocarbon components to the association between interspecific mating propensity and hydrocarbon profile divergences, we also tested for correlations between mating proportions and differences between species in their levels of each individual hydrocarbon (i.e., the five non log-transformed hydrocarbon components). Differences for 3Me27, 7Me27, and 13Me27 were significantly negatively correlated with interspecific mating propensity (Mantel’s r: -0.30 to −0.51, all p<0.02), whereas the two other components were not or only marginally non-significantly correlated with interspecific mating propensity (5Me27:r= 0.13, p = 0.77; 9Me27 + 11Me27: r= −0.24, p = 0.06). These patterns indicate that putative mating cues would likely be derived from global or combinatorial profiles rather than individual components, and/or that different components are used as signals by different species.
Macroevolutionary change in CHCs
Finally, we wanted to infer whether cuticular hydrocarbon profile differences between species develop gradually over time, as expected under constant intraspecific selection or neutral divergence, or whether there is accelerated change during speciation events. To this end, we first built a maximum-likelihood (ML) phylogeny describing the relationships between the nine Timema species used for hydrocarbon analysis (Figure 2). The branch lengths in this phylogeny can be interpreted as proportional to time, as the ML tree with a global clock constraint did not differ significantly from the best (unconstrained) ML tree (likelihood ratio test: p = 0.15). Maximum parsimony analyses and Bayesian inferences yielded the same tree topology (Figure 2), which is fully compatible with previous phylogenies of the genus [38, 39].
We then used this phylogeny to investigate the evolutionary tempo and mode of hydrocarbon profile divergence, represented by DF1 and DF2, as well as of the individual hydrocarbon components. In particular, we inferred whether hydrocarbon profiles tend to diverge continuously between species, or whether the amount of profile changeis correlated with the numberof speciation events . To this end, we evaluated the fit of nine different diffusion-based maximum-likelihood models, each representing a specific evolutionary scenario of character change. These scenarios include a 'neutral divergence' model whereby character change occurs on all branches in the phylogeny and is proportional to time, as well as different types of 'speciation' models where character change coincides with species splits (nodes) in the phylogeny (see Table 3 for a summary of models and interpretations and Methods and Additional file 1 for details). Akaike Information Criterion (AIC) values are calculated to determine which of the nine models of evolutionary change best characterizes the tempo and mode of diversification of hydrocarbon profiles in Timema.
The evolution of global hydrocarbon profiles (described by DF1 and DF2) were better explained by change occurring proportionally to speciation events (nodes in the phylogeny) than by neutral models in which change occurs proportionally to time (branch lengths). Specifically, the evolution of DF1 on the Timema phylogeny was best explained by punctuated models, in which at each node, one daughter branch retains the ancestral character value, and the other daughter branch changes(AIC differences between punctuated and other models>2; Table 3). For DF2, models with equal branch lengths provided a better fit to the data than any other model, but it was not possible to distinguish between nonphylogenetic and punctuated models (AIC differences <2; Table 3). The evolutionary patterns returned by the analyses for individual components varied widely. Support for speciational change was found for one of the five analyzed hydrocarbon components, 7Me27, where a phylogenetic model with equal branch lengths provided the best fit to the data (Table 3). Phenotypic change was therefore dependent on the number of speciation events but independent of time since divergence. By contrast, we found that evolutionary change in 3Me27 was best described by a pure-phylogenetic model, with change proportional to branch lengths, consistent with the significant correlation of LC3 differences between species and their genetic divergences. The same pattern was also revealed for 13Me27, as expected given the strong correlation between 13Me27 and 3Me27 (Pearson's product–moment correlation: rho= −0.95; t = −26.2, df = 77, p < 0.0001). The evolutionary change of the two remaining components (5Me27 and 9Me27 + 11Me27) was best described by non-phylogenetic models, which provide little information on processes driving character changes (Table 3). These diverse patterns for individual components could indicate that different species use different components as cues, or, more likely, that combinations of several components are used for species discrimination. The latter hypothesis is notably supported by DF1, representing a profile combination, and being best described by punctuational (i.e., speciational) models of change.
Altogether, our data thus indicate that cuticular hydrocarbon profiles vary among Timema species, and that most evolutionary change in global hydrocarbon profiles occurs in association with speciation events. Individual hydrocarbon components appear to change according to different scenarios, with some following patterns expected under neutral evolution and others with the amount of change more strongly correlated with speciation events than with divergence time. Similarities in hydrocarbon profiles between species are correlated with inter-specific mating propensities which in combination with the punctuational change of DF1 and DF2 on the Timema phylogeny supports the idea that hydrocarbon profiles may function as inter-specific mating signals in this group.
Cuticular hydrocarbon variation among T. cristinae populations
To test whether cuticular hydrocarbon profiles are also correlated with mating decisions at the intraspecific level, we combined an additional set of four populations of the species T. cristinae with the two T. cristinae populations included in the interspecific analyzes. This species was chosen because there is extensive evidence for host-plant associated mating preferences, and partial reproductive isolation, between populations occurring on Ceanothus host plants and populations occurring on Adenostoma host plants (e.g., [33–35]), but the proximate mechanisms underlying such mate discrimination have remained speculative.
A MANOVA revealed significant cuticular hydrocarbon profile differences between populations on different host plants but also differences between populations within hosts (populations nested in hosts; host effect: Wilks’ λ = 0.63, approx. F1,4 = 8.2, p = 0.046; populations within hosts: λ = 0.26, approx. F4,75 = 5.9, p = 0.0004). This multivariate difference appeared to be mainly due to variation in hydrocarbon component 13Me27 (82% variation between hosts, as opposed to 14% variation between populations within hosts, p = 0.048). None of the remaining hydrocarbon components displayed significantly more variation between populations on different hostplants as compared to populations from the same host-plant.
Discussion and conclusion
We have used a combination of data from female cuticular hydrocarbon profiles, interspecific mating trials, and phylogenetics to evaluate the potential role of cuticular hydrocarbons in mating and diversification of Timema walking sticks. Our primary findings are that (1) hydrocarbon profiles differ substantially between species, and between T. cristinae populations adapted to different host plants, such that species-specific and population-specific sets of hydrocarbons may serve as signals for mate choice between species and populations; (2) in no-choice trials, mating was significantly more likely between pairs of species with similar hydrocarbon profiles than between species with divergent profiles; and (3) the macroevolution of major components of hydrocarbon profiles (DF1) fits best with a model of punctuational phenotypic change, as expected under change driven by selection processes. Taken together, these convergent lines of evidence strongly suggest a role for cuticular hydrocarbon variation in mate choice and speciation in the genus Timema, although elucidating the exact nature of this role will require additional studies. Notably, the experimental manipulation of hydrocarbon profiles are required to demonstrate that these profiles are indeed used as mating cues in Timema, and will allow direct assessment of the importance of hydrocarbons for mating preferences at the individual level and at different degrees of reproductive isolation.
Among other insects, cuticular hydrocarbon variation has been implicated in a range of developmental, ecological and behavioral contexts, including ecological adaptation to the local abiotic or biotic environment , effects from temperature and nutrition during juvenile stages [41, 42], interactions with host plants [43, 44], and mate choice and sexual selection within species [45, 46]. Evidence for cuticular hydrocarbons mediating the evolution of reproductive isolation has thus far centered mainly on associations between variation in hydrocarbon profiles and variation in propensities to mate . For example, in Chorthippus parallelus grasshoppers, male hydrocarbon variation, but not other potential mate-choice signals, was significantly linked with the degree of assortative mating among individuals from 12 populations . Similarly, in Chrysochus beetles , and some Drosophila vinegar flies (e. g., [50, 51]), mating between individuals from closely-related species or strains is affected in part by variation in cuticular hydrocarbon profiles. In Timema, assuming individuals indeed detect aspects of hydrocarbon profiles, decreased propensity to accept mates with lower hydrocarbon profile similarity would limit the likelihood of potentially costly interspecific matings. Therefore, cuticular hydrocarbon profiles and preferences may act as effective isolation mechanisms among divergent and potentially post-zygotically incompatible populations and species. Alternatively, reproductive isolation, and cuticular hydrocarbon profiles, may evolve as byproducts of species divergence. The precise roles of hydrocarbon variation in processes underlying the evolution of reproductive isolation, and the timing of evolutionary changes in hydrocarbon profiles during the process of speciation, remain to be elucidated, and require integration of data on selective pressures, mechanisms, and underlying genetic changes from different temporal and spatial scales.
In this study, we have analyzed cuticular hydrocarbon variation across nine species, and multiple populations, at varying levels of phylogenetic and genetic divergence. Inference of the tempo and mode of hydrocarbon profile evolution provides evidence generally consistent with the hypothesis that change in this set of characters does not occur gradually during phyletic evolution, but relatively rapidly, in proportion to the number of speciation events. This pattern of speciational change has also been found for the macro evolution of Timema courtship behavior, but not for the macroevolution of male genitalic morphology, which has diversified in proportion to phyletic branch lengths, apparently under forces of sexual selection or sexual conflict . Among other insect species, previous studies of macroevolutionary change in hydrocarbon profiles have reported rapid evolution in conjunction with species specificity and a notable degree of phylogenetic conservatism [50, 53, 54], patterns that are largely concordant with those described here for Timema, but have not been partitioned into speciational versus phyletic concentrations of change. If changes in cuticular hydrocarbon profile drive speciation, then directional selection due to mate choice with profile-based criteria should characterize processes of divergence or reinforcement, while stabilizing selection should predominate at other times. Alternatively, hydrocarbon profiles may diverge between populations primarily due to differences in ecological variables such as temperature, humidity, habitat preference, or host-plant use, promoting the evolution of pre- and post-zygotic isolation under circumstances where such profiles also function in intraspecific or interspecific mate choices [29, 40].
In Timema, natural selection for cryptic coloration patterns that match those of their host plants represents a major force in phenotypic divergence among populations and species [34, 55–57]. Our intraspecific data showing hydrocarbon profile differences between T. cristinae from their two host plants, Adenostoma and Ceanothus, which are distributed in large- and small-scale mosaic patterns across their chaparral habitat, indicate that hydrocarbon-based chemical signals could be used as indicators of hostplant and other environmental adaptations  of a prospective mating partner. In our interspecific analyses of hydrocarbon profiles, we also found significant differences between the two populations sampled for T. poppensis and T. chumash (Figure 1). In these cases, the study of additional populations from each host plant may indicate whether these profile divergences are also correlated to specific hostplant adaptations as in T. cristinae.
The degree to which cuticular hydrocarbon profiles are genetically-based, compared to being acquired from the host plant or other features of the local environment, remains to be investigated in Timema; both heritable and environmental effects have been well established among other insects (e. g., [29, 40, 59]). Recent studies have identified sets of genes, mainly desaturases, that underlie variation in hydrocarbon profiles among Drosophila. Findings of rapid evolution and positive selection of such genes [60, 61] suggest that hydrocarbon evolution may exert strong, direct effects on phenotypic changes that mediate speciation. In Timema, the rapid changes in CHC profiles associated with speciation events and their correlation with reproductive isolation, provide further support for the idea that such traits can play an integral role in population divergence and speciation.
Insect collection and mating trials
We collected data on intra- and interspecific mating propensity, hydrocarbon profile variation and genetic divergences for nine of the 16 described sexual Timema species (T. bartmani, T. boharti, T. californicum, T. chumash, T. cristinae, T. knulli, T. petita, T. podura, and T. poppensis). These species include all but one sexual species from three subclades in the Timema genus, mainly distributed in central California. The six remaining sexual species not included in this study, from southernmost California, Arizona, and Mexico, are basal to these subdivisions [38, 62]. Such a near-complete coverage of species in a clade is necessary for valid tests of speciational change hypotheses, as missing species introduce errors into estimations of the number of speciation events separating a given species pair . Another factor affecting the inference of speciation events are species extinctions, which are assumed to be either negligible  or to have occurred at random on the phylogeny [19, 64].
Insects were collected using sweep nets between March and May of 2008 to 2011. Only individuals collected as juvenileswere used for analyses, and males and females were housed separately on their original host plant and raised to adults. This protocol ensured that all individuals used in experiments were virgins and of similar age. Different Timema species occur over latitudinally and altitudinally spread locations, such that the average developmental stage (number of molts to maturity) of individuals at a given date may vary greatly among populations. For the mating trials we needed individuals of each species that reached maturity at approximately the same time. For population combinations where developmental stages were very different, we therefore maintained juvenile males and females of the more advanced population in the refrigerator at 7°C for up to one week to slow down their development. Even though this treatment is unlikely to influence our results given only early juvenile stages were concerned, we nevertheless used these individuals only for mating trials, not for other experiments.
No-choice mating trials were conducted by introducing one male and one female into a 6-cm Petri dish. As in previous studies (e.g., [30, 33]), we recorded after 1 h whether the pair was copulating. Mating proportions for each male-species by female-species combination were then estimated as the proportion of copulating pairs out of the total number of trials conducted for that combination. Across the nine species (81 pairwise combinations) we conducted >1400 no-choice trials (Table 2).
Characterization of cuticular hydrocarbons
For each location, cuticular hydrocarbons of three to six adult virgin females (not used in the mating trials) were analyzed. Live individuals were anaesthetized by freezing for 1h, and submerged singly in 1 ml of HPLC-grade hexane for 5 min to extract the cuticular hydrocarbons from their body surface. The supernatant of each sample was then withdrawn, concentrated to circa 100 μl, and (E9)-octadecen-1-yl acetate (2 μl of a 100 ng/μl solution) was added as an internal standard (IS) for quantitative analyzes. The total amount of each target cuticular hydrocarbon (see below) was determined by multiplying the area count of the respective chromatographic peak (see below) with the 200 ng of the IS and by dividing the product by the area count of the IS.
Samples were analyzed with a Varian 3800 gas chromatograph (GC) coupled to a Varian Saturn Ion Trap mass spectrometer (MS) equipped with a DB-5 MS column (50 m ×0.25 mm i.d. for analyzes in 2009 and 2010; 30 m ×0.25 mm i.d. for analyzes in in 2011), using the following temperature program: 50°C for 2 min, then 20°C per min to 280°C (in 2009 and 2010) or to 240°C (in 2011). The final temperature of 280°C or 240°C was held for 46 min or 28 min, respectively. The injector temperature was 300°C. The mass spectrometer was set to scan for fragment ions between m/z 41 to m/z 500.
In 2009, analyses of body surface extracts of T. cristinae, T. knulli and T. poppensis revealed quantitative and qualitative differences in cuticular hydrocarbons between species, especially components eluting between heptacosane and octacosane. These components were selected as potential indicators of species-specific profiles. The components were hypothesized to be 3-, 5-, 7-, 9-, 11- and 13-methylheptacosane (henceforth 3Me27, 5Me27, 7Me27, 9Me27, 11Me27 and 13Me27, respectively) based on diagnostic fragment ions of their mass spectra . The presence of these fragment ions was confirmed by re-analyzes of the samples on a Hewlett Packard GCD Quadrupole mass spectrometer fitted with a DB-5 column (30 m × 0.25 mm i.d.), scanning for fragment ions between m/z 41 and m/z 425. To confirm the structural assignment for each of these components, previously known 3Me27 , 5Me27 [67, 68], and 7Me27, 9Me27, 11Me27 and 13Me27 [68, 69] were synthesized by coupling methylketones with phosphorus ylides, and by hydrogenating the resulting olefins. Identical mass spectra and retention times of each of the insect-produced and corresponding synthetic components confirmed all structural assignments. Because in a large number of samples it was not possible to completely separate 9Me27 and 11Me27, we pooled them for further analyses.
We analyzed proportional rather than absolute abundances of components, to remove species differences stemming from body size variation (up to threefold among Timema females) and to reduce experimental error [14, 70]. Proportional cuticular hydrocarbon components were calculated by dividing the amount of each component in a given sample by the sum of all components in that sample. These hydrocarbon proportions were then transformed using log-contrasts [14, 36] to remove the non-independence among analyzed variables. Log contrasts were calculated by dividing the value for each hydrocarbon by the component 13Me27, and then taking the log of these new variables, resulting in four log-contrast transformed values [log-contrasts for 3Me27, 5Me27, 7Me27 and 9/11Me27, referred to as LC3 (Log Contrast 3), LC5, LC7 and LC9-11, respectively] for every individual. Results using other components as the divisor were qualitatively similar.
To test for species and population differences in hydrocarbon profiles we performed a multivariate analysis of variance (MANOVA; ), with populations nested within species. To test for significance of the main factor (species) against the nested term, we used the nested.npmanova command in the BiodiversityR package  which evaluates the statistical significance of the F-ratios by permutation. To test for the effect of species, populations were randomized between species. To test for the effect of populations, individuals were randomized among populations. Given these permutation schemes, the indicated p-values correspond to the proportion of randomized datasets producing a larger or equal F-ratio than the original dataset. We calculated p-values for 1000 randomizations. Populations within each individual species were compared using parametric MANOVAS (based on Wilk’s test statistic).
The expression of hydrocarbon profiles as well as other signal traits, are known to be affected by environmental influences as well individual condition (e.g., [43, 73, 74]). While our sample of two populations per species, from different host plants in most cases, allows us to assess species differences while including intra-specific variation, the sampling scheme is not designed to investigate the influence of environmental effects. Thus, future work will be required to distinguish the relative contribution of plastic effects, especially those linked to host plant rearing environment, versus genetic influences on cuticular hydrocarbon profiles (see for example ).
To test whether interspecific matings were more likely between species with similar than divergent hydrocarbon profiles, we estimated profile divergences between species using the multivariate Euclidian distance with the species median values for each individual hydrocarbon component as a reference. In other words, each hydrocarbon component defines an axis in a multi-dimensional space, in which each species is represented by a cloud of points – the Euclidian distance between two species is the distance between the two cloud centers (given by the median value on each axis) in the multi-dimensional space. We used (partial) Mantel tests implemented in the R package vegan1.17-4  to test whether this distance is correlated with interspecific mating propensity (proportion of between species trials resulting in successful copulation) and genetic distance.
We next wanted to develop insights into the processes that shape hydrocarbon profile evolution among sexual Timema species. To this end, we evaluated the fit of different diffusion-based maximum-likelihood models, each representing a specific evolutionary scenario. In particular, we wanted to infer whether hydrocarbon profiles tend to diverge continuously between species, or whether profile changes occur in proportion to speciation events. These evolutionary scenarios can be evaluated by applying different evolutionary models of character change in a phylogenetic framework . We used the program CoMET  within the Mesquite system  to infer which of nine different character change models best describes the diversification of hydrocarbon profiles in Timema. Specifically, the nine models implemented in CoMET represent all possible combinations of three different phylogenetic models with three different tempo-of-change models, in a 3 × 3 matrix  (Additional file 1). The phylogenetic models describe the pattern of character change as pure-phylogenetic, non-phylogenetic or punctuational. Under the pure-phylogenetic model, character change occurs along all branches in the phylogeny so that the level of character divergence is correlated with the level of phylogenetic divergence. The non-phylogenetic model assumes a star phylogeny, which means that character divergence occurs independently of phylogenetic divergence. The punctuational model assumes that at each bifurcation in the tree, one daughter branch retains the ancestral character value, and the other daughter branch changes. The tempo-of-change models evaluate three different rates (distance, equal and free) at which the character can change on each branch. In the distance model, change is proportional to genetic distance (i.e., branch length). The equal model assumes that all branches in the phylogeny have the same lengths, and phenotypic change is therefore independent of time since divergence. In the free model, the length of each branch is proportional to the amount of character change occurred, rather than to genetic divergence. Under this model, the characters change at different rates in different lineages, and character change is thus not proportional to time since divergence. Overall, CoMET thus estimates the fit of nine distinct models of character change to the Timema phylogeny (Additional file 1), whereby Akaike Information Criterion (AIC) values are calculated to determine the model of evolutionary change that best characterizes the tempo and mode of diversification of hydrocarbon profiles in Timema. Note that while second order information criterion (AICc) values are usually preferred over AIC values for model selection (because under AICc, the relative penalty for model complexity increases for small datasets; ), AICc cannot be applied for the free rates tempo-of-change models, given the number of parameters estimated in these models (see Additional file 1 for details). This issue has, however, no influence on our conclusions as even with the less stringent penalizing for model complexity, we find the best fit for the simple models (in which only one parameter is estimated; see results).
Several previous studies have focused on species relationships in Timema[38, 39]. However, because the model evaluation with CoMET requires character values for each tip in the phylogeny, we built a new phylogeny using the nine species for which we determined hydrocarbon profiles and interspecific mating propensities. This phylogeny was based on concatenated mitochondrial COI and nuclear Hsp70 sequences (total of 1879 bp) that we generated for a previous study (see Additional file 2 for GenBank accession numbers). Maximum likelihood phylogenetic analyzes (with heuristic tree searches) were carried out using PAUP*4.0b10  on the freely available Bioportal (http://www.bioportal.uio.no). We used the optimal model of sequence evolution (GTR+I+G) identified with the Akaike information criterion (AIC) as implemented in jMODELTEST 0.1.1 . Branch support was evaluated by a maximum-likelihood bootstrapping analysis using Seqboot (500 replicates), DNAml, and Consense within the Phylip 3.68 package . For completeness, we also inferred Bayesian posterior probabilities with MrBayes3.1.2 , with Markov chains run for 106 generations, 104 generations burn-in values, and trees sampled every 100 generations.
To test whether hydrocarbon profile distances are also correlated with intraspecific mating propensities, we used four additional T. cristinae populations. Timema cristinae populations are morphologically adapted to different host plants (Adenostoma and Ceanothus; ), and there is extensive evidence for weak to moderate levels of sexual isolation between populations on different hosts [32–34]. However, sexual isolation is not based on morphology , and thus further work on the causes of sexual isolation is warranted. For logistic reasons, it was not possible to use the same methods as described for analyses of interspecific hydrocarbon profiles. Instead of using lab-raised virgin females, we directly collected 30 adult females from two Ceanothus populations and 41 from two Adentostoma populations. Although females thus represent a pool of different ages and mating status, this protocol is more likely to introduce error than systematic bias in terms of population divergence in hydrocarbon profiles. Cuticular hydrocarbons of each individual were directly extracted in the field and stored on ice until they were analyzed as described above. Because it was not possible to objectively separate 11Me27 and 13Me27 in these field samples, the two components were pooled; all remaining analyzes were conducted as described for interspecific comparisons.
Availability of supporting data
Data archived in the DRYAD repository under doi:10.5061/dryad.98f8c.
Andersson M: Sexual selection. 1994, Princeton: Princeton University Press
Arnqvist G, Nilsson T: The evolution of polyandry: multiple mating and female fitness in insects. Anim Behav. 2000, 60: 145-164. 10.1006/anbe.2000.1446.
Shaw KL: Interspecific genetics of mate recognition: Inheritance of female acoustic preference in Hawaiian crickets. Evolution. 2000, 54: 1303-1312.
Gerhardt HC: Acoustic spectral preferences in two cryptic species of grey treefrogs: implications for mate choice and sensory mechanisms. Anim Behav. 2005, 70: 39-48. 10.1016/j.anbehav.2004.09.021.
Blows MW: Interaction between natural and sexual selection during the evolution of mate recognition. Proc R Soc Lond B. 2002, 269: 1113-1118. 10.1098/rspb.2002.2002.
Boake CRB: Sexual selection and speciation in Hawaiian Drosophila. Behav Genet. 2005, 35: 297-303. 10.1007/s10519-005-3221-4.
Eberhard WG: Sexual Selection and Animal Genitalia. 1985, Cambridge, MA: Harvard University Press
McPeek MA, Shen L, Torrey JZ, Farid H: The tempo and mode of three-dimensional morphological evolution in male reproductive structures. Am Nat. 2008, 171: E158-E178. 10.1086/587076.
Turelli M, Barton NH, Coyne JA: Theory and speciation. Trends Ecol Evol. 2001, 16: 330-343. 10.1016/S0169-5347(01)02177-2.
Boake CRB: Flying apart: Mating behavior and speciation. BioScience. 2000, 50: 501-508. 10.1641/0006-3568(2000)050[0501:FAMBAS]2.0.CO;2.
Coyne JA, Orr HA: Speciation. 2004, Sunderland, MA: Sinauer Associates
Edward DA, Chapman T: The evolution and significance of male mate choice. Trends Ecol Evol. 2011, 26: 647-654. 10.1016/j.tree.2011.07.012.
Maan ME, Seehausen O: Ecology, sexual selection and speciation. Ecol Lett. 2011, 14: 591-602. 10.1111/j.1461-0248.2011.01606.x.
Blows MW, Allan RA: Levels of mate recognition within and between two Drosophila species and their hybrids. Am Nat. 1998, 152: 826-837. 10.1086/286211.
Butlin RK: Reinforcement of premating isolation. Speciation and Its Consequences. Edited by: Otte D, Endler JA. 1989, Sunderland, MA: Sinauer, 158-179.
Lande R: Natural-selection and random genetic drift in phenotypic evolution. Evolution. 1976, 30: 314-334. 10.2307/2407703.
Lynch M: The rate of morphological evolution in mammals from the standpoint of the neutral expectation. Am Nat. 1990, 136: 727-741. 10.1086/285128.
Albert AYK, Schluter D: Reproductive character displacement of male stickleback mate preference: Reinforcement or direct selection?. Evolution. 2004, 58: 1099-1107.
Mooers AO, Vamosi SM, Schluter D: Using phylogenies to test macroevolutionary hypotheses of trait evolution in Cranes (Gruinae). Am Nat. 1999, 154: 249-259. 10.1086/303226.
Butlin RK, Hewitt GM, Webb SF: Sexual selection for intermediate optimum in Chorthippus brunneus (Orthoptera, Acrididae). Anim Behav. 1985, 33: 1281-1292. 10.1016/S0003-3472(85)80188-3.
Barlow GW, Siri P: Does sexual selection account for the conspicuous head dimorphism in the Midas cichlid?. Anim Behav. 1997, 53: 573-584. 10.1006/anbe.1996.0332.
Gavrilets S: Rapid evolution of reproductive barriers driven by sexual conflict. Nature. 2000, 403: 886-889. 10.1038/35002564.
McPeek MA, Gavrilets S: The evolution of female mating preferences: Differentiation from species with promiscuous males can promote speciation. Evolution. 2006, 60: 1967-1980.
Arnqvist G, Rowe L: Antagonistic coevolution between the sexes in a group of insects. Nature. 2002, 415: 787-789. 10.1038/415787a.
McPeek MA: The macroevolutionary consequences of ecological differences among species. Palaeontology. 2007, 50: 111-129. 10.1111/j.1475-4983.2006.00616.x.
Van der Meer RK, Breed MD, Espelie KE, Winston M: Pheromone communication in social insects: ants, wasps, bees and termites. 1998, Boulder, CO: Westview
Ferveur JF: Cuticular hydrocarbons: Their evolution and roles in Drosophila pheromonal communication. Behav Genet. 2005, 35: 279-295. 10.1007/s10519-005-3220-5.
Nojima S, Schal C, Webster FX, Santangelo RG, Roelofs WL: Identification of the sex pheromone of the German cockroach, Blattella germanica. Science. 2005, 307: 1104-1106. 10.1126/science.1107163.
Howard RW, Blomquist GJ: Ecological, behavioral, andbiochemical aspects of insect hydrocarbons. Annu Rev Entomol. 2005, 50: 371-393. 10.1146/annurev.ento.50.071803.130359.
Arbuthnott D, Crespi BJ: Courtship and mate discrimination within and between species of Timema walking-sticks. Anim Behav. 2009, 78: 53-59. 10.1016/j.anbehav.2009.02.028.
Vickery VR: Revision of Timema Scudder (Phasmatoptera, Timematodea) including 3 new species. Can Entomol. 1993, 125: 657-692. 10.4039/Ent125657-4.
Nosil P: Reproductive isolation caused by visual predation on migrants between divergent environments. Proc R Soc Lond B. 2004, 271: 1521-1528. 10.1098/rspb.2004.2751.
Nosil P, Crespi BJ, Sandoval CP: Host-plant adaptation drives the parallel evolution of reproductive isolation. Nature. 2002, 417: 440-443. 10.1038/417440a.
Nosil P: Divergent host plant adaptation and reproductive isolation between ecotypes of Timema cristinae walking sticks. Am Nat. 2007, 169: 151-162. 10.1086/510634.
Nosil P, Crespi BJ, Sandoval CP: Reproductive isolation driven by the combined effects of ecological adaptation and reinforcement. Proc R Soc Lond B. 2003, 270: 1911-1918. 10.1098/rspb.2003.2457.
Aitchison J: The statistical analysis of compositional data. 1986, London -New York: Chapman and Hall, 12
Schwander T, Crespi BJ: Multiple direct transitions from sexual reproduction to apomictic parthenogenesis in Timema stick insects. Evolution. 2009, 63: 84-103. 10.1111/j.1558-5646.2008.00524.x.
Law JH, Crespi BJ: The evolution of geographic parthenogenesis in Timema walking-sticks. Mol Ecol. 2002, 11: 1471-1489. 10.1046/j.1365-294X.2002.01547.x.
Schwander T, Henry L, Crespi BJ: Molecular evidence for ancient asexuality in Timema stick insects. Curr Biol. 2011, 21: 1129-1134. 10.1016/j.cub.2011.05.026.
Buckley SH, Tregenza T, Butlin RK: Transitions in cuticular composition across a hybrid zone: historical accident or environmental adaptation?. Biol J Linn Soc. 2003, 78: 193-201. 10.1046/j.1095-8312.2003.00147.x.
Toolson EC, Markow TA, Jackson LL, Howard RW: Epicuticular hydrocarbon composition of wild and laboratory-reared Drosophila mojavensis Patterson and Crow (Diptera Drosophilidae). Ann Entomol Soc Am. 1990, 83: 1165-1176.
Etges WJ, Veenstra CL, Jackson LL: Premating isolation is determined by larval rearing substrates in cactophilic Drosophila mojavensis. VII. Effects of larval dietary fatty acids on adult epicuticular hydrocarbons. J Chem Ecol. 2006, 32: 2629-2646. 10.1007/s10886-006-9187-8.
Etges WJ, De Oliveira CC, Ritchie MG, Noor MAF: Genetics of incipient speciation in Drosophila mojavensis: II. Host plants and mating status influence cuticular hydrocarbon QTL expression and G x E interactions. Evolution. 2009, 63: 1712-1730. 10.1111/j.1558-5646.2009.00661.x.
Grace T, Wisely SM, Brown SJ, Dowell FE, Joern A: Divergent host plant adaptation drives the evolution of sexual isolation in the grasshopper Hesperotettix viridis (Orthoptera: Acrididae) in the absence of reinforcement. Biol J Linn Soc. 2010, 100: 866-878. 10.1111/j.1095-8312.2010.01458.x.
Thomas ML, Simmons LW: Sexual selection on cuticular hydrocarbons in the Australian field cricket, Teleogryllus oceanicus. BMC Evol Biol. 2009, 9: 162-10.1186/1471-2148-9-162.
Everaerts C, Farine JP, Cobb M, Ferveur JF: Drosophila cuticular hydrocarbons revisited: mating status alters cuticular profiles. PLoS One. 2010, 5: e9607-10.1371/journal.pone.0009607.
Wicker-Thomas C: Evolution of insect pheromones and their role in reproductive isolation and speciation. Ann Soc Entomol Fr. 2011, 47: 55-62. 10.1080/00379271.2011.10697696.
Tregenza T, Pritchard VL, Butlin RK: Patterns of trait divergence between populations of the meadow grasshopper, Chorthippus parallelus. Evolution. 2000, 54: 574-585.
Peterson MA, Dobler S, Larson EL, Juarez D, Schlarbaum T, Monsen KJ, Francke W: Profiles of cuticular hydrocarbons mediate male mate choice and sexual isolation between hybridising Chrysochus (Coleoptera : Chrysomelidae). Chemoecology. 2007, 17: 87-96. 10.1007/s00049-007-0366-z.
Alves H, Rouault JD, Kondoh Y, Nakano Y, Yamamoto D, Kim YK, Jallon JM: Evolution of cuticular hydrocarbons of Hawaiian Drosophilidae. Behav Genet. 2010, 40: 694-705. 10.1007/s10519-010-9364-y.
Grillet M, Everaerts C, Houot B, Ritchie MG, Cobb M, Ferveur JF: Incipient speciation in Drosophila melanogaster involves chemical signals. Sci Rep. 2012, 2: 224-
Arbuthnott D, Elliot MG, McPeek MA, Crespi BJ: Divergent patterns of diversification in courtship and genitalic characters of Timema walking-sticks. J Evol Biol. 2010, 23: 1399-1411. 10.1111/j.1420-9101.2010.02000.x.
Mullen SP, Mendelson TC, Schal C, Shaw KL: Rapid evolution of cuticular hydrocarbons in a species radiation of acoustically diverse Hawaiian crickets (Gryllidae : Trigonidiinae : Laupala). Evolution. 2007, 61: 223-231. 10.1111/j.1558-5646.2007.00019.x.
De Oliveira CC, Manfrin MH, Sene FD, Jackson LL, Etges WJ: Variations on a theme: diversification of cuticular hydrocarbons in a clade of cactophilic Drosophila. BMC Evol Biol. 2011, 11: 179-10.1186/1471-2148-11-179.
Sandoval CP, Nosil P: Counteracting selective regimes and host preference evolution in ecotypes of two species of walking-sticks. Evolution. 2005, 59: 2405-2413.
Nosil P, Sandoval CP: Ecological niche dimensionality and the evolutionary diversification of stick insects. PLoS One. 2008, 4: e1907-
Sandoval CP, Crespi BJ: Adaptive evolution of cryptic coloration: the shape of host plants and dorsal stripes in Timema walking-sticks. Biol J Linn Soc. 2008, 94: 1-5. 10.1111/j.1095-8312.2007.00941.x.
Nosil P, Gompert Z, Farkas TE, Comeault AA, Feder JL, Buerkle CA, Parchman TL: Genomic consequences of multiple speciation processes in a stick insect. Proc R Soc Lond B. 2012, 279: 5058-5065. 10.1098/rspb.2012.0813.
Gleason JM, James RA, Wicker-Thomas C, Ritchie MG: Identification of quantitative trait loci function through analysis of multiple cuticular hydrocarbons differing between Drosophila simulans and Drosophila sechellia females. Heredity. 2009, 103: 416-424. 10.1038/hdy.2009.79.
Shirangi TR, Dufour HD, Williams TM, Carroll SB: Rapid evolution of sex pheromone-producing enzyme expression in Drosophila. PLoS Biol. 2009, 7: e1000168-10.1371/journal.pbio.1000168.
Keays MC, Barker D, Wicker-Thomas C, Ritchie MG: Signatures of selection and sex-specific expression variation of a novel duplicate during the evolution of the Drosophila desaturase gene family. Mol Ecol. 2011, 20: 3617-3630.
Sandoval CP, Carmean DA, Crespi BJ: Molecular phylogenetics of sexual and parthenogenetic Timema walking-sticks. Proc R Soc Lond B. 1998, 265: 589-595. 10.1098/rspb.1998.0335.
Hey J: Using phylogenetic trees to study speciation and extinction. Evolution. 1992, 46: 627-640. 10.2307/2409633.
Harvey PH, Purvis A: Comparative methods for explaining adaptations. Nature. 1991, 351: 619-624. 10.1038/351619a0.
Krkosova Ž, Kubinec R, Soják L, Amann A: Temperature-programmed gas chromatography linear retention indices of all C-4-C-30 monomethylalkanes on methylsilicone OV-1 stationary phase - Contribution towards a better understanding of volatile organic compounds in exhaled breath. J Chromatogr A. 2008, 1179: 59-68. 10.1016/j.chroma.2007.10.081.
Holman L, Jørgensen CG, Nielsen J, D'Ettorre P: Identification of an ant queen pheromone regulating worker sterility. Proc R Soc Lond B. 2010, 277: 3793-3800. 10.1098/rspb.2010.0984.
Lafrance D, Shani A, Margalit J: Biological activity of synthetic hydrocarbon mixtures of cuticular components of the female housefly (Musca domestica L.). J Chem Ecol. 1989, 15: 1475-1490. 10.1007/BF01012377.
Hemptinne JL, Lognay G, Dixon AFG: Mate recognition in the two-spot ladybird beetle, Adalia bipunctata: role of chemical and behavioural cues. J Insect Physiol. 1998, 44: 1163-1171. 10.1016/S0022-1910(98)00081-X.
Stransky K, Streibl M, Sorm F: Über Naturwachse IV.: Über einen neuen Typ verzweigter Paraffine aus dem Wachs der Honigbiene (Apis mellifera L.). Collect Czechoslov Chem Commun. 1966, 31: 4694-4702. 10.1135/cccc19664694.
Rundle HD, Chenoweth SF, Doughty P, Blows MW: Divergent selection and the evolution of signal traits and mating preferences. PLoS Biol. 2005, 3: 1988-1995.
Legendre L, Legendre P: La structure des données écologiques, vol. 2, 2 edn. Paris et les Presses de l'Université du. 1984, Québec: Masson
Kindt R, Coe R: Tree diversity analysis. A manual and software for common statistical methods for ecological and biodiversity studies. 2005, Nairobi: World Agroforestry Centre (ICRAF)
Danielson-Francois AM, Kelly JK, Greenfield MD: Genotype x environment interaction for male attractiveness in an acoustic moth: evidence for plasticity and canalization. J Evol Biol. 2006, 19: 532-542. 10.1111/j.1420-9101.2005.01006.x.
Rodriguez-Munoz R, Bretman A, Hadfield JD, Tregenza T: Sexual selection in the cricket Gryllus bimaculatus: no good genes?. Genetica. 2008, 134: 129-136. 10.1007/s10709-008-9250-1.
Vegan : Community Ecology Package. R package version 1.17-4 http://CRAN.R-project.org/package=vegan
Lee CH, Blay S, Mooers AO, Singh A, Oakley TH: CoMET: A Mesquite package for comparing models of continuous character evolution on phylogenies. Evol Bioinforma. 2006, 2: 183-186.
Maddison WP, Maddison DR: Mesquite: a modular system for evolutionary analysis. 2011, Version 2.01 http://mesquiteproject.org
Burnham K, Anderson D: Model selection and multimodel inference: a practical information-theoretic approach. 2002, New York: Springer
Swofford DL: Phylogenetic analysis using parsimony and other methods 4.0b10. 2000, Sunderland, MA: Sinauer Associates, Inc
Posada D: jModelTest: Phylogenetic Model Averaging. Mol Biol Evol. 2008, 25: 1253-1256. 10.1093/molbev/msn083.
Felsenstein J: PHYLIP (Phylogeny Inference Package) version 3.5c. 1993, Seattle: Department of Genetics, University of Washington
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.
Sandoval CP: Differential visual predation on morphs of Timema cristinae (Phasmatodeae, Timemidae) and its consequences for host-range. Biol J Linn Soc. 1994, 52: 341-356. 10.1111/j.1095-8312.1994.tb00996.x.
Nosil P, Crespi BJ: Does gene flow constrain adaptive divergence or vice versa? A test using ecomorphology and sexual isolation in Timema cristinae walking-sticks. Evolution. 2004, 58: 102-112.
We thank Grigori Khaskin for methylheptacosane syntheses and Aaron Comeault and Timothy Farkasfor help with collecting and extracting cuticular hydrocarbons for T. cristinae. This project was supported by funding from the NWO to TS, NSERC to BJC and NSERC-IRC to GG, with Contech Enterprises Inc. and Global Forest Sciences as sponsors.
The authors declare that they have no competing interests.
TS and BJC planned the study, RG and GG performed hydrocarbon experiments, DA and TS conducted interspecific mating tests, TS analyzed the data, TS and BJC wrote the manuscript with input from all authors. All authors read and approved the final manuscript.