Hydrocarbon divergence and reproductive isolation in Timemastick insects
© Schwander et al.; licensee BioMed Central Ltd. 2013
Received: 20 November 2012
Accepted: 12 July 2013
Published: 16 July 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
Sampling locations and host plant information for interspecific analyzes, as well as for intraspecific comparisons in T. cristinae
Hydrocarbon variation and between-species mating propensities
Number of mating trials per species combination
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].
Summary of the nine character change models, and the fit of each model to Timema hydrocarbon profiles (represented by DF1 and DF2) and individual hydrocarbon components
Model for trait change
Conclusion if best fit
9Me27 + 11Me27
Time predicts the amount of change occurred (consistent with neutral divergence).
Trait not involved in speciation
The amount of change depends on the number of speciation events occurred (number of nodes).
Consistent with speciational change2
Trait values can change at any rate between speciation events
Trait not involved in speciation
Closely related species share trait values for a short time and then diverge very rapidly
Trait not involved in speciation
Trait values change very rapidly, with similar rates in different lineages
Trait values change very rapidly, with different rates in different lineages
At each speciation event, one daughter species retains the ancestral trait value, the trait in the other daughter species changes, with the amount of change dependent on time
Consistent with speciational change2, level of divergence between species also affected by the time separating them
At each speciation event, one daughter species retains the ancestral trait value, the trait in the other daughter species changes, the amount of change between speciation events is always the same
Consistent with speciational change2
At each speciation event, one daughter species retains the ancestral trait value, the trait in the other daughter species changes, the amount of change between speciation events varies freely
Consistent with speciational change2, level of divergence between species also affected by lineage-specific processes
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.
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.
- Andersson M: Sexual selection. 1994, Princeton: Princeton University PressGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Shaw KL: Interspecific genetics of mate recognition: Inheritance of female acoustic preference in Hawaiian crickets. Evolution. 2000, 54: 1303-1312.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Boake CRB: Sexual selection and speciation in Hawaiian Drosophila. Behav Genet. 2005, 35: 297-303. 10.1007/s10519-005-3221-4.PubMedView ArticleGoogle Scholar
- Eberhard WG: Sexual Selection and Animal Genitalia. 1985, Cambridge, MA: Harvard University PressView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Turelli M, Barton NH, Coyne JA: Theory and speciation. Trends Ecol Evol. 2001, 16: 330-343. 10.1016/S0169-5347(01)02177-2.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Coyne JA, Orr HA: Speciation. 2004, Sunderland, MA: Sinauer AssociatesGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Maan ME, Seehausen O: Ecology, sexual selection and speciation. Ecol Lett. 2011, 14: 591-602. 10.1111/j.1461-0248.2011.01606.x.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Butlin RK: Reinforcement of premating isolation. Speciation and Its Consequences. Edited by: Otte D, Endler JA. 1989, Sunderland, MA: Sinauer, 158-179.Google Scholar
- Lande R: Natural-selection and random genetic drift in phenotypic evolution. Evolution. 1976, 30: 314-334. 10.2307/2407703.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Albert AYK, Schluter D: Reproductive character displacement of male stickleback mate preference: Reinforcement or direct selection?. Evolution. 2004, 58: 1099-1107.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Gavrilets S: Rapid evolution of reproductive barriers driven by sexual conflict. Nature. 2000, 403: 886-889. 10.1038/35002564.PubMedView ArticleGoogle Scholar
- McPeek MA, Gavrilets S: The evolution of female mating preferences: Differentiation from species with promiscuous males can promote speciation. Evolution. 2006, 60: 1967-1980.PubMedView ArticleGoogle Scholar
- Arnqvist G, Rowe L: Antagonistic coevolution between the sexes in a group of insects. Nature. 2002, 415: 787-789. 10.1038/415787a.PubMedView ArticleGoogle Scholar
- McPeek MA: The macroevolutionary consequences of ecological differences among species. Palaeontology. 2007, 50: 111-129. 10.1111/j.1475-4983.2006.00616.x.View ArticleGoogle Scholar
- Van der Meer RK, Breed MD, Espelie KE, Winston M: Pheromone communication in social insects: ants, wasps, bees and termites. 1998, Boulder, CO: WestviewGoogle Scholar
- Ferveur JF: Cuticular hydrocarbons: Their evolution and roles in Drosophila pheromonal communication. Behav Genet. 2005, 35: 279-295. 10.1007/s10519-005-3220-5.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Vickery VR: Revision of Timema Scudder (Phasmatoptera, Timematodea) including 3 new species. Can Entomol. 1993, 125: 657-692. 10.4039/Ent125657-4.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Nosil P, Crespi BJ, Sandoval CP: Host-plant adaptation drives the parallel evolution of reproductive isolation. Nature. 2002, 417: 440-443. 10.1038/417440a.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Aitchison J: The statistical analysis of compositional data. 1986, London -New York: Chapman and Hall, 12View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.Google Scholar
- Tregenza T, Pritchard VL, Butlin RK: Patterns of trait divergence between populations of the meadow grasshopper, Chorthippus parallelus. Evolution. 2000, 54: 574-585.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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-PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Sandoval CP, Nosil P: Counteracting selective regimes and host preference evolution in ecotypes of two species of walking-sticks. Evolution. 2005, 59: 2405-2413.PubMedView ArticleGoogle Scholar
- Nosil P, Sandoval CP: Ecological niche dimensionality and the evolutionary diversification of stick insects. PLoS One. 2008, 4: e1907-View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Hey J: Using phylogenetic trees to study speciation and extinction. Evolution. 1992, 46: 627-640. 10.2307/2409633.View ArticleGoogle Scholar
- Harvey PH, Purvis A: Comparative methods for explaining adaptations. Nature. 1991, 351: 619-624. 10.1038/351619a0.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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: MassonGoogle Scholar
- 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)Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.Google Scholar
- Maddison WP, Maddison DR: Mesquite: a modular system for evolutionary analysis. 2011, Version 2.01 http://mesquiteproject.orgGoogle Scholar
- Burnham K, Anderson D: Model selection and multimodel inference: a practical information-theoretic approach. 2002, New York: SpringerGoogle Scholar
- Swofford DL: Phylogenetic analysis using parsimony and other methods 4.0b10. 2000, Sunderland, MA: Sinauer Associates, IncGoogle Scholar
- Posada D: jModelTest: Phylogenetic Model Averaging. Mol Biol Evol. 2008, 25: 1253-1256. 10.1093/molbev/msn083.PubMedView ArticleGoogle Scholar
- Felsenstein J: PHYLIP (Phylogeny Inference Package) version 3.5c. 1993, Seattle: Department of Genetics, University of WashingtonGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.