Phenotypic disparity in Iberian short-horned grasshoppers (Acrididae): the role of ecology and phylogeny
© The Author(s). 2017
Received: 30 January 2017
Accepted: 19 April 2017
Published: 4 May 2017
The combination of model-based comparative techniques, disparity analyses and ecomorphological correlations constitutes a powerful method to gain insight into the evolutionary mechanisms that shape morphological variation and speciation processes. In this study, we used a time-calibrated phylogeny of 70 Iberian species of short-horned grasshoppers (Acrididae) to test for patterns of morphological disparity in relation to their ecology and phylogenetic history. Specifically, we examined the role of substrate type and level of ecological specialization in driving different aspects of morphological evolution (locomotory traits, chemosensitive organs and cranial morphology) in this recent radiation.
We found a bimodal distribution of locomotory attributes corresponding to the two main substrate type guilds (plant vs. ground); plant-perching species tend to exhibit larger wings and thicker femora than those that remain on the ground. This suggests that life form (i.e., substrate type) is an important driving force in the evolution of morphological traits in short-horned grasshoppers, irrespective of ancestry. Substrate type and ecological specialization had no significant influence on head shape, a trait that showed a strong phylogenetic conservatism. Finally, we also found a marginal significant association between the length of antennae and the level of ecological specialization, suggesting that the development of sensory organs may be favored in specialist species.
Our results provide evidence that even in taxonomic groups showing limited morphological and ecological disparity, natural selection seems to play a more important role than genetic drift in driving the speciation process. Overall, this study suggests that morphostatic radiations should not necessarily be considered as “non-adaptive” and that the speciation process can bind both adaptive divergence mechanisms and neutral speciation processes related with allopatric and/or reproductive isolation.
Adaptive radiations, groups that have rapidly diversified from a common ancestor to exploit a wide suite of ecological niches, have intrigued evolutionary biologists for over a century [1, 2]. Considered as important biodiversity engines, these bursts of speciation are thought to arise from ecological opportunity in the form of vacant ecological niches that become available due to the colonization of new environments (spatial dispersal) or the acquisition of adaptive innovations that allow the access to novel niche dimensions (ecological dispersal) [3, 4]. Spatial and/or ecological dispersal can be driven by modifications of existing environments via climatic changes [4, 5]. When lineages first enter these new adaptive zones, morphological evolution should initially be rapid; as niche spaces become increasingly saturated, the rate of morphological evolution would be expected to slow down [6, 7]. A core prediction based on the above scenario (the so-called early-burst model) is the existence of early rapid diversification followed by a slowdown in net diversification over time [8–10]. Thereby, lineage and morphological diversification are frequently positively correlated [11, 12], but the opposite pattern (i.e., a negative relationship) has also been reported [13, 14]. Nevertheless, recent studies have shown that processes underlying phenotypic disparity and those generating species diversity can be uncoupled, suggesting that ecological opportunity is not the only diversification force [15–17]. For example, speciation by simple geographic isolation can generate a pattern of declining speciation trough time without the intervention of niche-filling processes . On the other hand, speciation bursts can occur in multiple pulses linked to environmental changes (e.g., pulses of orogenic uplift or glacial retreat), thus eroding the diversity-dependent signature of diversification. This may explain why a large number of studies performed on radiating clades have failed to detect early-bursts of phenotypic evolution, despite some of them represent the most classic examples of adaptive radiation, including Darwin’s finches from the Galapagos Islands and Anolis lizards of the Caribbean Islands (reviewed in ).
Recently, it has been recognized that there are many different types of evolutionary radiations and not all adaptive radiations conform to the early-burst model [19, 20]. For example, under a repeated radiation scenario, it is expected that subclades tend to resemble each other in morphological disparity and, contrary to that predicted under the early-burst paradigm, morphological evolution should not necessarily be initially rapid. This kind of adaptive radiation (“iterative radiation”) may arise in systems dominated by constraints, such as when evolution within clades is driven by repeated adaptation to similar environments [21–23]. This model of radiation exemplifies that episodes of ecological opportunity can be recurring over the evolutionary history of a lineage . Beyond the epithet “adaptive”, non-adaptive radiations (or morphostatic radiations) constitute another form of radiation that remains largely unstudied [25–27]. Non-adaptive radiations arise through processes that are unrelated to niche exploitation and the resulting species usually exhibit low morphological disparity and allopatric distributions . For example, non-adaptive radiations driven by sexual selection result in new species that are ecologically similar to their ancestors such are the case of Hawaiian Laupala crickets . The boundaries of adaptive and non-adaptive radiations are sometimes diffuse and it has been suggested that both processes are extremes along a continuum . Thus, some lineages may exhibit features consistent with these two concepts [31, 32].
In the present study, we used a time-calibrated phylogeny of 70 species of Iberian short-horned grasshoppers to test for patterns of morphological disparity in relation to their ecology and phylogenetic history. Specifically, we aimed at testing whether the major morphological changes occurred early in the acridid’s diversification history or if random-walk patterns (“morphological drift”) or other processes (e.g., evolutionary constraints) have driven the diversification of phenotypic attributes in this group. We also examined the potential role of substrate use and the level of ecological specialization (i.e., niche breadth of a given species) in shaping morphological attributes (locomotory traits, chemosensitive organs and cranial morphology) in this invertebrate radiation. If substrate type and ecological specialization imposes a selective pressure across taxa, it would result in the same solution -phenotypic optima- multiple times, leading to convergence among species even between distantly-related clades.
The Iberian Peninsula is renowned for its high level of biodiversity and number of endemic species, and also as one of the main refugial areas in Europe during the Pleistocene Ice Ages [42, 43]. Most of species included in this study are endemic to the Iberian Peninsula or have a distribution restricted to Iberia, France and North Africa. Our dataset included 70 taxa of short-horned grasshoppers belonging to four different subfamilies; slant-faced grasshoppers (subfamily Gomphocerinae, 43 spp.), band-winged grasshoppers (subfamily Oedipodinae, 17 spp.), spur-throated grasshoppers (Catantopinae, 5 spp.) and 5 spp. belonging to different subfamilies (Calliptaminae ×3 spp., Dericorythinae and Eyprepocnemidinae) (Additional file 1: Table S1). Our dataset accounts for around three quarters of all extant species of acridid grasshoppers that have been recorded in the Iberian Peninsula . Iberian acridids do not descend from a single common ancestor and, thus, they do not constitute a local radiation in this region, but a paraphyletic group . However, it should be noted that several recent studies have pointed out that analyses of adaptive radiation do not require to include all descendants of a particular common ancestor, as such requisite is overly restrictive (see e.g. ). In fact, no theory of adaptive radiation predicts monophyly .
We sequenced four mitochondrial gene fragments: cytochrome c oxidase subunit 1 (COI), NADH dehydrogenase subunit 5 (ND5), 12S rRNA (12S) and a fragment containing parts of 16S rRNA (16S). For some taxa we failed to obtain reliable sequences, so we complemented our dataset with sequences obtained from previous studies [45, 46]. Sequences were aligned in MAFFT online version 7  and concatenated using Sequencematrix 1.7.8. . We calculated the best-fit models of nucleotide substitution for each of the four genes using jModelTest 0.1.1 . The TIM2 + I + Γ substitution model was selected for 12S, GTR + I+ Γ for 16S, TrN + I+ Γ model for ND5 and TPM3uf + I+ Γ was selected for COI. More details about these procedures are given in .
Phylogenetic inference was carried out under both maximum likelihood (ML) and Bayesian frameworks. Maximum likelihood analyses were conducted with two search replicates and 1000 bootstrap replicates using GARLI version 2.0 . Bayesian inference analyses were conducted in BEAST 1.8.0  in order to estimate a time-calibrated phylogeny. We used an uncorrelated lognormal relaxed-clock model and applied a Yule process as tree prior. Two calibration points were used in order to calibrate the tree based on absolute times. We employed as a first calibration point the split between Gomphocerinae and Oedipodinae, estimated to have occurred ~100 Mya ago. This estimate is based on dated ancient cockroach fossils . As a second calibration point, we used the divergence between Sphingonotus azurescens (mainland species) and S. guanchus (endemic to La Gomera Island, Canary Islands), whose divergence is estimated to have occurred around 3.5 Mya . Sphingonotus guanchus was only included in BEAST analyses for calibration purposes. Note that dating should be treated with caution as our calibrations are based on previous estimates. We ran analyses for a total of 100 × 105 generations, with a sampling frequency of 1000 generations. We ensured that replicated analyses converged (effective sample size values >200) using Tracer 1.4.1. Tree and log files (9500 trees after a 5% burn-in) of the two runs were combined with LogCombiner 1.4.7 , and the maximum clade credibility (MCC) tree was compiled with TreeAnnotator 1.4.7. The obtained tree topologies were rather similar and consistent with previous studies . See  for more details.
i) Locomotory morphology
Our dataset included 316 preserved specimens representing 3 to 5 male adults for 70 species of short-horned grasshoppers. We gathered a dataset of five morphological attributes of known correlation with locomotor performance and ecology . Morphological measurements included: (1) structural (i.e., head + thorax) body length, (2) tibia length, (3) femur width, (4) femur length, and (5) forewing (tegmina) length. All measurements were taken by the same observer (VGN) using a ZEISS stereomicroscope (SteREO Discovery V.8; Carl Zeiss Microscopy GmbH, Germany). Because disparity in body size and level of sexual size dimorphism has been addressed in a previous study  and the former trait (body size) is strongly correlated with other morphological variables, we size-corrected the remaining four variables using structural body length as size measurement. The relative length and relative width of the femur are inversely correlated (see Additional file 1: Figure S1), so we calculated a width/length femur ratio. As a result, our final dataset was composed of three variables related with locomotory morphology: tibia length, femur width/length ratio and tegmina length. Tegmina are modified leathery forewings present in some insects like grasshoppers. The major role of this structure is that of protecting the hindwings (and therefore, hind- and forewing length are highly correlated; p < 0.001, R 2 = 0.98). These also have an aerodynamic function and thus can be considered an informative trait about the dispersal potential of a species . On the other hand, leg length is thought to determine the jumping performance of insects like grasshoppers, leafhoppers and froghoppers [57–59]. Thus, relative femur length can provide information on jumping ability in insects that use a catapult-like mechanism to jump. We performed a phylogenetic principal component analysis (pPCA) on these three variables using the R package phytools . The phylogenetic principal component analysis yielded a single principal component (hereafter pPClm) that explained 53.2% of the variance of our three morphological variables (loadings; tibia length: −0.518, width/length femur ratio: −0.737, tegmina length: −0.886). The positive extreme of pPClm represents species with short wings and tibiae and thicker femora (e.g., Pezotettix giornae, Podisma carpetana), and the negative extreme of this axis represents species with large wings and tibiae and more stylized femora (e.g., Chorthippus jucundus, Stethophyma grossum). We also implemented a non-phylogenetic principal component analysis (PCA) as a recent study has suggested that using phylogenetic PC axes as trait data could bias results . However, we found similar results by using both approaches (correlation pPCA vs. PCA; p = 0.028, r = 0.33).
ii) Antenna length
In conjunction with variables directly related to locomotor performance we also measured antenna length in a similar way as described above. Antenna length was corrected for body size, so we refer to relative antenna length hereafter. Members of the family Acrididae show relatively short and stout antennae in comparison with crickets and katydids. However, the length of antennae largely differs among taxa and these differences may be linked to the diet and degree of ecological specialization of each species . Acridid grasshoppers recognize their host plant by sensorial stimuli, which are perceived through chemoreceptors located on the labrum and antennae . Chemosensory sensilla on antennae express two classes of proteins involved in the recognition of odors and taste, and thus the antennae are thought to play an important role in host orientation and food selection . In this vein, we would expect a greater development of this structure in taxa with a narrower ecological niche (i.e., specialist species) and species preferring to perch on plants, which may display higher sensitivity to plant odour perception .
iii) Head shape
We aimed to examine the evolutionary relationship between grasshopper morphology (locomotory morphology, relative antenna length, and head shape) and substrate type. For this purpose, we categorized each species on the basis of substrate preference (ground vs. plant) from the literature (e.g. [68, 69]) and our own personal observations. In addition, we examined the association between our morphological variables and the level of ecological specialization (i.e., the niche breadth of a given species). Ecological specialization refers to the ability of a species to exploit a range of resources and its capacity to use each one as result of evolutionary trade-offs (the “jack-of-all-trades is master of none” hypothesis; [70, 71]). Here, we used the ‘Paired Difference Index’ (PDI) as an estimator of ecological specialization [72, 73]. PDI values were calculated from a species-habitat matrix in which we rated the level of association of each species (from complete generalist, 0, to complete specialist, 3) with the nine most common habitats in which the studied species can be found. The PDI is a robust specialization index which takes into account not only the number of resources used by a species, but also the strength of the association between the species and its resources . Scores of species-habitat association were obtained directly from the literature and our own personal observations. PDI values were computed using the R package bipartite .
Phylogenetic comparative analyses
We assessed the phylogenetic signal of our focal variables (locomotory morphology, relative antenna length, and head shape) in order to test if these traits tracked the evolutionary history of the group and were not randomly distributed across taxa. We computed Pagel’s lambda (λ; ) for the first two principal components of head shape (PChs and PChs2), and for the two remaining variables, locomotory morphology (pPClm) and relative antenna length, in phytools . Pagel’s λ can take a value from 0, which means that phylogeny has no impact on the distribution of the trait and the values can be treated as independent, to a value of 1, which means that phylogeny fully predicts the distribution of the trait. We then mapped the evolution of these continuous traits on our phylogeny by using the function ‘contMap’ in phytools, which estimates states at internal nodes using ML and interpolates the states along each edge following .
We examined the association between the two ecological factors (substrate type and degree of ecological specialization) and our morphological variables (relative antenna length, PChs and pPClm). First, in order to test if substrate type (“ground-perching” vs. “plant-perching”; categorical variable) has influence on these variables, we performed phylogenetic ANOVAs using the function ‘phy.anova’ (10,000 simulations) in the R package geiger . Because substrate type largely varies among subfamilies (most gomphocerine species are plant-perching species whereas most of the remaining species stay on the ground) we performed a second phylogenetic ANOVA only including Gomphocerinae species (n = 43). Thereby, we assessed the influence of substrate type on morphology while controlling for the non-independence of the data because of shared ancestry. Secondly, we analyzed the relationship between our morphological variables and the level of ecological specialization (PDI index; continuous variable) using phylogenetic generalized least squares (PGLSλ). PGLS regression analyses were performed using the caper package  and graphically visualized by means of phylogenetically independent contrasts (PICs) computed using the PDAP:PDTREE module in Mesquite v.3.04 .
In order to visualize the relationship between phylogeny and taxon distribution in the morphological space, we built a ‘phylomorphospace’, which projects the branches of a phylogenetic tree into a 2D morphospace defining trait variation among species . When subclades occupy limited morphospace it means that most disparity is accounted for by early divergence (early-burst of trait evolution), while if subclades occupy large regions of the morphospace (i.e., stronger overlap) it means that morphological diversity (‘disparity’) is largely explained by recent divergence. We generated a phylomorphospace by plotting head shape against locomotion morphology (PChs vs. pPClm) using phytools . We also generated a second morphospace from the two first principal components for head shape (PChs vs. PChs2).
We investigated how morphological disparity has accumulated over the group’ evolutionary history. We evaluated the rate of morphological evolution in relation to lineage diversification by means of disparity-through-time (DTT) analyses , as implemented in the geiger package . We compared observed relative disparity with average phenotypic disparity simulated under Brownian Motion (1000 simulations) in order to test whether disparity was accumulated during the early or recent history of the group. When the rate of morphological evolution is constant, as expected under BM, the DTT curve is expected to decline linearly toward zero through time. Under an early-burst model, disparity will decline sharply and much earlier, while if evolution within subclades is fast the observed DTT curve will fall above the Brownian profile. In order to quantify the magnitude of the difference between the null disparity profile computed under BM, and the observed disparity profile from our data, we computed the so-called morphological disparity index (MDI) statistic using the geiger package . A positive MDI value is indicative of greater than expected subclade disparity (i.e., disparity distributed primarily within subclades) whereas negative MDI values indicate greater than expected subclade disparity (i.e., disparity distributed primarily among subclades). A negative MDI is characteristic of adaptive radiations, as rapidly diversifying taxa are expected to evolve distinct morphologies in response to new adaptive zones and slow once niches are filled. DTT analyses were performed twice; by including all taxa and only including species of the subfamily Gomphocerinae. We then ran a node height test , which tests for accelerations or decelerations in trait evolution, by comparing the independent contrasts for a trait with the respective node height (i.e., relative age). A significant negative relationship between absolute contrast value and node age supports the hypothesis of adaptive radiations as it would imply that species are dividing niche space more finely through time, consistent with a niche-filling model. Node height tests were performed using the R package caper .
Tempo and mode of morphological evolution
Evolutionary models are conceived to infer the different possible processes shaping phenotypic evolution and provide a method of testing different predictions regarding the tempo and mode of evolution. We assessed the fit of four evolutionary models to our phenotypic variables using the functions ‘fitContinuous’ and ‘OUwie’ implemented in the R packages geiger  and OUwie , respectively. We fitted three different single-rate models; Brownian-Motion (i.e., diffusive drift), Early-Burst (i.e., exponential declining of evolutionary rates) and single-optimum Ornstein-Uhlenbeck (i.e., bounded evolution around a single phenotypic optimum), and a multi-peak Ornstein-Uhlenbeck model. Brownian-motion (BM) operates under the assumption that trait evolution proceeds a random walk wherein trait variance across lineages accumulates proportional to time (single-rate model; ). The Early-Burst (EB) model predicts rapid morphological disparity early in the radiation, followed by a slowdown in the diversification rate as ecological niches are filled over time . Support for the single-optimum OU model (OU1) would imply that there is a single phenotypic optima (θ) for all taxa. This model is often associated with a process of stabilizing selection in which variation of phenotypic traits revolves around one or more stationary peaks . Lastly, we assessed the fit of a multivariate OU model (OUVM) with separate morphological optima and separate random walk variances (σ2) for each substrate guild (1: plant, 2: ground), and one global parameter (α), which determines the strength of selection towards those optima. For the multi-peak OU model we determined the possible ancestral substrate regimes along the internal branches of the phylogenetic tree using the Stochastic Mutational Mapping on Phylogenies (SIMMAP) tool in phytools, and sampled 500 character histories in order to incorporate evolutionary uncertainty. We compared the fit of the BM, EB, OU1 and OUMV models using the Akaike information criterion corrected for small sample size (AICc), which can be employed to compare models that differ in the number of parameters and therefore have non-comparable likelihoods.
Phylogenetic signal and ancestral state reconstruction
All morphological variables except PChs2 showed a significant phylogenetic signal suggesting conservative evolution for these traits (relative antenna length, λ = 0.837, p < 0.001; pPClm, λ = 0.575, p = 0.003; PChs, λ = 0.883, p < 0.001; PChs2, λ ~ 0, p ~ 1). Visual simulation of trait evolution using ancestral character reconstruction confirmed this pattern; species in the same subfamily tend to share morphological traits (see e.g. Additional file 1: Figures S2-S3). However, within Oedipodinae, there is a greater morphological variability than in the Gomphocerinae group.
Tempo and mode of morphological evolution
Comparison of four evolutionary model fits, random-walk variances (σ2) and primary trait optima (in parentheses) for the three morphological variables: head shape (PChs), locomotory morphology (pPClm) and relative antenna length
σ2 plant (θplant)
σ2 ground (θground)
(a) All taxa
(b) Gomphocerinae subset
Like some examples of non-adaptive radiations (Plethodon: , Rattus: ), acridid grasshoppers tend to retain a rather conserved body plan with little overt ecomorphological specialization among taxa (Fig. 6). However, we observed a greater morphological disparity within Oedipodinae in comparison with that observed for the Gomphocerinae, suggesting that the level of morphological resemblance among species differ between the two subfamilies. Gomphocerinae (slant-faced grasshoppers) constitute a very recent radiation, which is suspected to be mostly the result of divergent evolution of isolated populations induced by climatic oscillations during the Pleistocene . Some gomphocerine species spread out from southern refugia and expanded their ranges northwards during interglacials, whereas many others (montane species) became restricted to high altitudes (“sky-islands”) due to their limited dispersal capacity [35, 86]. Hence, many Iberian species probably arose in association with these mountain ranges (Sistema Central, Betic ranges, Picos de Europa; ). At a smaller scale, the evolution of elaborate acoustic signals as premating reproductive isolation mechanisms (via sexual selection) may have allowed gomphocerine grasshoppers to diversify in regions in which ecologically similar and phylogenetically related species coexist [46, 86, 87]. Thus, it is likely that both divergence in allopatry [42, 86] and prezygotic isolation mechanisms [46, 88, 89] have played a crucial role in promoting speciation in this subfamily. However, our results should be interpreted cautiously because the present study assesses the dynamic of morphological evolution among a small fraction of the Acrididae’ overall diversity. Thus, our findings may reflect local processes (i.e. factors that have affected which species can occur together in the same regional theatre) rather than global mechanisms that operate across the entire clade.
Patterns of morphological disparity through time
Our results do not support the existence of an early-burst of phenotypic diversification, adding to the number of study-cases showing no support for this evolutionary model [6, 7, 90, 91]. Obtained MDI values indicate that substantial amount of variance is clustered within subclades, pointing out to the existence of low phylogenetic niche conservatism. That is, there was no tendency of closely related species to be more similar to each other in ecomorphological traits than they are to more distant relatives as expected if some processes had constrained niche divergence among phylogenetically close species [92, 93]. This is in contrast to that predicted for taxa that experience early-bursts of diversification, where the partitioning of morphological disparity through time (disparification sensu ) is expected to be low [10, 82, 95]. Overall, DTT plots showed a general tendency for relative disparity in phenotypic attributes to decrease over time, but with two main pulses of increases in disparification during the Cenozoic and the Pleistocene. The observed pattern is not consistent with that predicted under an iterative radiation scenario wherein it should be expected the appearance of repeated diversification peaks . The fact that acridids partitioned morphological disparity within rather than among clades suggests that acridid lineages did not evolve along distinct morphological trajectories through time. Thus, it is unlikely that acridid lineages explored different adaptive zones, which makes sense taking into account the limited morphological variability they exhibit.
Evolutionary models of morphological divergence
Overall, results from the multivariate model fitting analysis summarized in Table 1 favored an Ornstein-Uhlenbeck model, indicating bounded evolution around a single (OU1) or two phenotypic optima (OUVM). Specifically, in relation to head shape, we found that a single model does not uniquely explain the evolution of this trait for the entire dataset; BM and OUMV were similarly informative. However, when models of evolution were tested exclusively on the Gomphocerinae subset, BM was the most informative model. It suggests that head shape evolution displays an idiosyncratic component and that, irrespective of substrate on which they rest, Gomphocerinae species have a more conical head shape in comparison with band-winged grasshoppers due to their phylogenetic legacy. On the other hand, our results indicate that locomotory morphology in Acrididae is inconsistent with a BM process, and has not evolved under a constant rate over time. A multivariate-peak OU (OUMV) model provided the best fit in both cases (when considering all taxa and after including only Gomphocerinae), which means that the evolution of this trait revolves around adaptive peaks.
In this context, it is possible that biomechanical constraints, which underlie much of locomotory morphology variation, impose restrictions on morphological diversification [96–98]. Specifically, our results suggest that stabilizing selection pulls the locomotory morphology towards two convergent adaptive optima (θplant, θground) during the group’s evolutionary history. A key difference between the two substrate type regimes in terms of locomotory attributes relies on the different degree of development of femora (which has influence on jumping performance) and wings (which determines flight capability) that each guild of species exhibit (see more below). Interestingly, we observed a negative relationship between relative forewing length and relative femur width/length ratio in Gomphocerinae (see Additional file 1: Figure S5) pointing out to the existence of a trade-off between flying- and jumping-based dispersal investments in this group. With regard to antenna length, a single-peak OU model was the model that provided the best fit to the observed data, suggesting that antennae evolution has remained constrained by directional selection.
We investigated the role of substrate type and level of ecological specialization in driving morphological evolution in our study group. The variation in structural characteristics within the habitat (substrate type, size and incline) has direct effects on animal locomotion and can lead to the evolution of morphological and behavioral adaptations (e.g., anti-predator strategies) . For example, stick insects (order Phasmatodea) of New Caledonia and New Guinea (“land lobsters”) are ground-dwelling species and show a distinct ecomorph; they are flightless and exhibit a stocky body form and thick hind legs. This is in contrast with the majority of stick insects, which are solitary canopy-dwellers and show an instantly recognizable morphology (stick-like bodies with large and thin legs) . In particular, substrate type may be especially important for small animals that jump using a catapult-like system such as grasshoppers and leafhoppers . Here, we found that plant-perching grasshoppers tend to exhibit a slender morphology in comparison with those species that prefer to stay on the ground, which show a more compressed body (shorter and thicker femora). This result is consistent with that reported by  in a study with leafhoppers. These authors suggested that insect species that jump from plants have longer legs because it reduces the amount of energy lost to bending the leaf, whereas species that jump from a stiffer substrate can have shorter-legs . Hence, in species living on plants, natural selection would have favored long hind legs because these allow short preparation times when an emergency jump is necessary . The existence of moderate differences in morphology between ground-dwelling and plant-perching species was confirmed when restricting our analyses to the Gomphocerinae subset. Despite of members of the subfamily Gomphocerinae exhibit low variability in morphology (and as a result of this phenotypic conservatism some species can be almost only distinguished by song attributes), gomphocerine species perching on plants tend to present a less rotund morphology -larger and thinner legs- in comparison with gomphocerine species that prefer to stay on the ground, which show squat bodies and short limbs. Our results suggest that there are two possible ways to achieve superior jumping capabilities: to evolve longer limbs or to evolve more muscular limbs (as stronger legs will produce more acceleration) . This finding is in agreement with that reported by  in a study with anole species, which also seem to face trade-offs that prevent them from simultaneously optimizing different skills of jumping ability (larger vs. thicker legs).
Head shape differed significantly between ground- and plant-species across taxa. Plant-perching species show a sharper head than ground-species, which may be favored by natural selection as this head shape facilitates camouflage against potential predators by making it difficult to see the mimic against the surroundings (i.e., background matching) . However, this difference between both categories become non-significant within the Gomphocerinae, which implies that differences in head shape between ground- and plant-species have to do with its different phylogenetic history per se rather than with the possible existence of a substrate type effect. In band-winged grasshoppers (subfamily Oedipodinae), the orientation of the face is usually nearly vertical (straight head shape) whereas species of the subfamily Gomphocerinae (also known as slant-faced grasshoppers) show a sharper profile, that is, a more conical head.
Lastly, we also found a marginally significant association between the relative length of antennae and the level of ecological specialization. Larger antennae can bear a higher number of chemoreceptors, which would be favorable for grasshoppers in the search for shelter or food . Hence, in those species that can thrive only in a limited range of habitats or have a limited diet (specialist species), the development of sensory organs may be favored, whereas in species with a broader ecological niche breadth (generalist species) selection for larger antennae may be weaker. Thus, the length of antennae in short-horned grasshoppers could be employed as a proxy for the degree of ecological specialization in terms of habitat requirements of a given species.
Short-horned grasshoppers exploit a variety of habitats, but lack the extreme specializations (“ecomorphs”) observed in other taxa (see Fig. 1). Accordingly, at a first glance, it seems that acridid grasshoppers exhibit a conservative and “all-purpose” morphology allowing them to perform equally well under different environments, which is contrary to one of the premises that are frequently used to catalogue a radiation as adaptive [19, 105]. However, in this study we provide evidence that natural selection could lead to a convergent body plan in very distant clades in this recent radiation . Specifically, we found substrate type-related macroevolutionary variation in locomotory attributes and head morphology. The observed pattern was consistent with an evolutionary scenario with two adaptive peaks (selective agents) related to life-form (i.e., substrate type; plant vs. ground). Thus, these findings suggest that, although this radiation does not meet the conditions to be classified as adaptive sensu stricto (at least in comparison with textbook examples of rapid adaptive radiations such as Madagascan vangas and tetragnathid spiders on the Hawaiian Islands), it provides evidence that natural selection can act on a very small scale, and therefore it is not always detectable. So, even in a priori morphostatic radiations, adaptive processes can lead to subtle phenotypic variation. Hence, our results support the notion that ecological factors play a key role in evolutionary processes by acting as selective mechanisms and by imposing constraints that shape morphological traits . In sum, this study highlights that natural selection, and not genetic drift, constitutes the main factor driving the disparification process in this recent radiation despite of the lack of apparent phenotypic and ecological variability. This calls into the question the use of term “non-adaptive” as synonymous with “morphostatic” since the total absence of an adaptive component seems to be unlikely in any radiation.
We thank A. Hidalgo, C. Cáliz and M.J. González-Serna for their help with field work and genetic analyses. D. Chobanov and one anonymous reviewer provided useful and valuable comments on an earlier draft of this manuscript. We acknowledge support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).
VGN was supported by a “Juan de la Cierva” fellowship (FPDI-2013-16,828). VN was supported by a FPI PhD studentship (BES-2012-053741). JO was supported by a “Ramón y Cajal” research fellowship (RYC-2013-12,501). This work received financial support from the Spanish Ministry of Economy and Competitiveness (CGL2014-54671-P).
Availability of data and materials
The concatenated matrix and the MCC tree are available in LabArchives (DOI 10.6070/H4FF3QFM and 10.670/H49S1P35, respectively). All other data supporting the results of this article are included within the article and its Additional files.
VGN had the original idea for the study. VGN and JO designed the research. JO, PJC and VN collected the samples. JO supervised the genetic analyses. VGN carried out the analyses, interpreted the results and wrote the manuscript with inputs from JO, PJC and VN. All authors read and approved the final version.
The authors declare that they have no competing interests.
Consent for publication
This study did not require ethical approval. Specimens were collected under license from Spanish regional authorities.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Givnish TJ, Systma KJ. Molecular evolution and adaptive radiation. Cambridge: Cambridge University Press; 1997.Google Scholar
- Grant PR, Grant BR. How and why species multiply: the radiation of Darwin’s finches. Princeton University Press: Princeton; 2008.Google Scholar
- Simpson GG. The major features of evolution. New York: Columbia University Press; 1953.Google Scholar
- Schluter D. The ecology of adaptive radiation. Oxford: Oxford University Press; 2000.Google Scholar
- Losos JB. Adaptive radiation, ecological opportunity, and evolutionary determinism. Am Nat. 2010;175:623–39.PubMedView ArticleGoogle Scholar
- Harmon LJ, Losos JB, Davies J, Gillespie RG, Gittleman JL, Jennings WB, Kozak K, McPeek MA, Moreno-Roark F, Near TJ, Purvis A, Ricklefs RE, Schluter D, Schulte JA II, Seehausen O, Sidlauskas B, Torres-Carvajal O, Weir JT, Mooers AØ. Early bursts of body size and shape evolution are rare in comparative data. Evolution. 2010;64:2385–96.PubMedGoogle Scholar
- Mahler DL, Revell LJ, Glor RE, Losos JB. Ecological opportunity and the rate of morphological evolution in the diversification of greater Antillean anoles. Evolution. 2010;64:2731–45.PubMedView ArticleGoogle Scholar
- Rabosky DL, Lovette IJ. Explosive evolutionary radiations: decreasing speciation or increasing extinction through time? Evolution. 2008;62:1866–75.PubMedView ArticleGoogle Scholar
- Moen D, Morlon H. Why does diversification slow down? Trends Ecol Evol. 2014;29:190–7.PubMedView ArticleGoogle Scholar
- Colombo M, Damerau M, Hanel R, Salzburger W, Matschiner M. Diversity and disparity through time in the adaptive radiation of Antarctic notothenioid fishes. J Evol Biol. 2015;28:376–94.PubMedPubMed CentralView ArticleGoogle Scholar
- Price SA, Wainwright PC, Bellwood DR, Kazancioglu E, Collar DC, Near TJ. Functional innovations and morphological diversification in parrotfish. Evolution. 2010;64:3057–68.PubMedGoogle Scholar
- Rabosky DL, Santini F, Eastman J, Smith SA, Sidlauskas B, Chang J, Alfaro ME. Rates of speciation and morphological evolution are correlated across the largest vertebrate radiation. Nat Commun. 2013;4:1958.PubMedView ArticleGoogle Scholar
- Claramunt S, Derryberry EP, Brumfield RT, Remsen JV Jr. Ecological opportunity and diversification in a continental radiation of birds: climbing adaptations and cladogenesis in the Furnariidae. Am Nat. 2012;179:649–66.PubMedView ArticleGoogle Scholar
- Thacker CE. Species and shape diversification are inversely correlated among gobies and cardinalfishes (Teleostei: Gobiiformes). Org Divers Evol. 2014;14:419–36.View ArticleGoogle Scholar
- Derryberry EP, Claramunt S, Derryberry G, Chesser RT, Cracraft J, Aleixo A, Pérez-Emán J, Remsen JV Jr, Brumfield RT. Lineage diversification and morphological evolution in a large-scale continental radiation: the Neotropical ovenbirds and Woodcreepers (Aves: Furnariidae). Evolution. 2011;65:2973–86.PubMedView ArticleGoogle Scholar
- Hipsley CA, Miles DB, Müller J, Mu J. Morphological disparity opposes latitudinal diversity gradient in lacertid lizards. Biol Lett. 2014;10:20140101.PubMedPubMed CentralView ArticleGoogle Scholar
- Alhajeri BH, Schenk JJ, Steppan SJ. Ecomorphological diversification following continental colonization in muroid rodents (Rodentia: Muroidea). Biol J Linn Soc. 2016;117:463–81.View ArticleGoogle Scholar
- Pigot AL, Phillimore AB, Owens I, Orme CDL. The shape and temporal dynamics of phylogenetic trees arising from geographic speciation. Syst Biol. 2010;59:660–73.PubMedView ArticleGoogle Scholar
- Glor RE. Phylogenetic insights on adaptive radiation. Annu Rev Ecol Evol Syst. 2010;41:251–70.View ArticleGoogle Scholar
- Simões M, Breitkreuz L, Alvarado M, Baca S, Cooper JC, Heins L, Herzog K, Lieberman BS. The evolving theory of evolutionary radiations. Trends Ecol Evol. 2016;31:27–34.PubMedView ArticleGoogle Scholar
- Seehausen O. African cichlid fish: a model system in adaptive radiation research. Proc R Soc Lond B Biol Sci. 2006;273:1987–98.View ArticleGoogle Scholar
- Losos JB. Convergence, adaptation, and constraint. Evolution. 2011;65:1827–40.PubMedView ArticleGoogle Scholar
- Frédérich B, Sorenson L, Santini F, Slater GJ, Alfaro ME. Iterative ecological radiation and convergence during the evolutionary history of damselfishes (Pomacentridae). Am Nat. 2013;181:94–113.PubMedView ArticleGoogle Scholar
- Pincheira-Donoso D, Harvey LP, Ruta M. What defines an adaptive radiation? Macroevolutionary diversification dynamics of an exceptionally species-rich continental lizard radiation. BMC Evol Biol. 2015;15:153.PubMedPubMed CentralView ArticleGoogle Scholar
- Cameron RAD, Cook LM, Hallows JD. Land snails on Porto Santo: adaptive and non-adaptive radiation. Philos Trans R Soc Lond Ser B Biol Sci. 1996;351:309–27.View ArticleGoogle Scholar
- Wilke T, Benke M, Brändle M, Albrecht C, Bichain JM. The neglected side of the coin: non-adaptive radiations in spring snails (Bythinella spp.) In: Glaubrecht M, editor. Evolution in action. Berlin: Springer; 2010. p. 551–78.View ArticleGoogle Scholar
- Wellenreuther M, Sánchez-Guillén RA. Nonadaptive radiation in damselflies. Evol Appl. 2015;9:103–18.PubMedPubMed CentralView ArticleGoogle Scholar
- Gittenberger E. Radiation and adaptation, evolutionary biology and semantics. Org Divers Evol. 2004;4:135–6.View ArticleGoogle Scholar
- Mendelson TC, Shaw KL. Sexual behaviour: rapid speciation in an arthropod. Nature. 2005;433:375–6.PubMedView ArticleGoogle Scholar
- Olson ME, Arroyo-Santos A. Thinking in continua: beyond the “adaptive radiation” metaphor. BioEssays. 2009;31:1337–46.PubMedView ArticleGoogle Scholar
- Gillespie R. Community assembly through adaptive radiation in Hawaiian spiders. Science. 2004;303:356–9.PubMedView ArticleGoogle Scholar
- Rowe KC, Aplin KP, Baverstock PR, Moritz C. Recent and rapid speciation with limited morphological disparity in the genus Rattus. Syst Biol. 2011;60:188–203.PubMedView ArticleGoogle Scholar
- Soulebeau A, Aubriot X, Gaudeul M, Rouhan G, Hennequin S, Haevermans T, Dubuisson JY, Jabbour F. The hypothesis of adaptive radiation in evolutionary biology: hard facts about a hazy concept. Org Divers Evol. 2015;15:747–61.View ArticleGoogle Scholar
- Song H, Amédégnato C, Cigliano MM, Desutter-Grandcolas L, Heads SW, Huang Y, Otte D, Whiting MF. 300 million years of diversification: elucidating the patterns of orthopteran evolution based on comprehensive taxon and gene sampling. Cladistics. 2015;31:621–51.View ArticleGoogle Scholar
- Mayer F, Berger D, Gottsberger B, Schulze W. Non-ecological radiations in acoustically communicating grasshoppers? In: Glaubrecht M, editor. Evolution in action. Case studies in adaptive radiation, speciation and the origin of biodiversity. Berlin: Springer Verlag; 2010. p. 451–64.Google Scholar
- Kozak KH, Weisrock DW, Larson A. Rapid lineage accumulation in a non-adaptive radiation: phylogenetic analysis of diversification rates in eastern north American woodland salamanders (Plethodontidae: Plethodon). Proc R Soc Lond B Biol Sci. 2006;273:539–46.View ArticleGoogle Scholar
- Blankers T, Adams DC, Wiens JJ. Ecological radiation with limited morphological diversification in salamanders. J Evol Biol. 2012;25:634–46.PubMedView ArticleGoogle Scholar
- Moen DS, Irschick DJ, Wiens JJ. Evolutionary conservatism and convergence both lead to striking similarity in ecology, morphology and performance across continents in frogs. Proc R Soc Lond B Biol Sci. 2013;280:20132156.View ArticleGoogle Scholar
- Vidal-García M, Byrne PG, Roberts JD, Keogh JS. The role of phylogeny and ecology in shaping morphology in 21 genera and 127 species of Australo-Papuan myobatrachid frogs. J Evol Biol. 2014;27:181–92.PubMedView ArticleGoogle Scholar
- Vidal-García M, Keogh JS. Convergent evolution across the Australian continent: ecotype diversification drives morphological convergence in two distantly related clades of Australian frogs. J Evol Biol. 2015;28:2136–51.PubMedView ArticleGoogle Scholar
- Edwards S, Tolley KA, Vanhooydonck B, Measey GJ, Herrel A. Is dietary niche breadth linked to morphology and performance in Sandveld lizards Nucras (Sauria: Lacertidae)? Biol J Linn Soc. 2013;110:674–88.View ArticleGoogle Scholar
- Hewitt GM. Some genetic consequences of ice ages, and their role in divergence and speciation. Biol J Linn Soc. 1996;58:247–76.View ArticleGoogle Scholar
- Gómez A, Lunt DH. Refugia within refugia: patterns of phylogeographic concordance in the Iberian Peninsula. In: Weiss S, Ferrand N, editors. Phylogeography of southern European refugia. Dordrecht: Springer; 2007. p. 155–88.View ArticleGoogle Scholar
- Presa JJ, García MD, Clemente ME. Catalogue of Orthoptera Caelifera from the Iberian Peninsula and Balearic Islands (Orthoptera: Caelifera). J Orthoptera Res. 2007;16:175–9.View ArticleGoogle Scholar
- Contreras D, Chapco W. Molecular phylogenetic evidence for multiple dispersal events in gomphocerines grasshoppers. J Orthoptera Res. 2006;15:91–8.View ArticleGoogle Scholar
- Nattier R, Robillard T, Amedegnato C, Couloux A, Cruaud C, Desutter-Grandcolas L. Phylogeny and evolution of acoustic communication in Gomphocerinae (Orthoptera: Caelifera: Acrididae). Zool Scr. 2011;40:479–97.View ArticleGoogle Scholar
- Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.PubMedPubMed CentralView ArticleGoogle Scholar
- Vaidya G, Lohman DJ, Meier R. SequenceMatrix: concatenation software for the fast assembly of multi-gene datasets with character set and codon information. Cladistics. 2011;27:171–80.View ArticleGoogle Scholar
- Posada D. jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008;25:1253–6.PubMedView ArticleGoogle Scholar
- García-Navas V, Noguerales V. Cordero PJ, Ortego J. Ecological drivers of body size evolution and sexual size dimorphism in short-horned grasshoppers (Orthoptera: Acrididae). bioRxiv/2017/119560. doi:https://doi.org/10.1101/119560.
- Bazinet AL, Zwickl DJ, Cummings MP. A gateway for phylogenetic analysis powered by grid computing featuring GARLI 2.0. Syst Biol. 2014;63:812–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.PubMedPubMed CentralView ArticleGoogle Scholar
- Fries M, Chapco W, Contreras D. A molecular phylogenetic analysis of the Oedipodinae and their intercontinental relationships. J Orthoptera Res. 2007;16:115–25.View ArticleGoogle Scholar
- Husemann M, Depperman J, Hochkirch A. Multiple independent colonization of the Canary Islands by the winged grasshopper genus Sphingonotus Fieber, 1852. Mol Phylogenet Evol. 2014;81:174–81.PubMedView ArticleGoogle Scholar
- Chapman RF, Joern A. Biology of grasshoppers. New York: Wiley and Sons; 1990.Google Scholar
- Heidinger IMM, Hein S, Bonte D. Patch connectivity and sand dynamics affect dispersal-related morphology of the blue-winged grasshopper Oedipoda caerulescens in coastal grey dunes. Insect Conserv Diver. 2010;3:205–12.Google Scholar
- Burrows M. Froghopper insects leap to new heights. Nature. 2003;424:509.PubMedView ArticleGoogle Scholar
- Burrows M, Sutton GP. The effect of leg length on jumping performance of short- and long-legged leafhopper insects. J Exp Biol. 2008;211:1317–25.PubMedView ArticleGoogle Scholar
- Bennet-Clark HC. Jumping in Orthoptera. In: Chapman RF, Joern A, editors. Biology of grasshoppers. New York: John Wiley and Sons; 1990. p. 173–203.Google Scholar
- Revell LJ. Two new graphical methods for mapping trait evolution on phylogenies. Methods Ecol Evol. 2013;4:754–9.View ArticleGoogle Scholar
- Uyeda JC, Caetano DS, Pennell MW. Comparative analysis of principal components can be misleading. Syst Biol. 2015;64:677–89.PubMedView ArticleGoogle Scholar
- Uvarov BP. Grasshoppers and locusts. A handbook of general Acridology vol 2. London: Centre for Overseas Pest Research; 1977.Google Scholar
- Yu Y, Zhang S, Zhang L, Zhao X. Developmental expression of odorant-binding proteins and chemosensory proteins in the embryos of Locusta migratoria. Arch Insect Biochem Physiol. 2009;71:105–15.PubMedView ArticleGoogle Scholar
- Chen HH, Zhao YX, Kang L. Antennal sensilla of grasshoppers (Orthoptera: Acrididae) in relation to food preferences and habits. J Biosci. 2003;28:743–52.PubMedView ArticleGoogle Scholar
- Barros FC, Herrel A, Kohldorf T. Head shape evolution in Gymnophthalmidae: does habitat use constrain the evolution of cranial design in fossorial lizards? J Evol Biol. 2011;24:2423–33.PubMedView ArticleGoogle Scholar
- Adams DC, Otarola-Castillo E. Geomorph: an R package for the collection and analysis of geometric morphometric shape data. Methods Ecol Evol. 2013;4:393–9.View ArticleGoogle Scholar
- Rohlf FJ, Slice DE. Extensions of the Procrustes method for the optimal superimposition of landmarks. Syst Zool. 1990;39:40–59.View ArticleGoogle Scholar
- Harz K. Orthopteren Europas/the Orthoptera of Europe Vol. II. The Hague: Dr. W. Junk; 1975.View ArticleGoogle Scholar
- Llucià-Pomares D. Revisión de los ortópteros (Insecta: Orthoptera) de Cataluña (España). Monografías SEA. 2002;7:1–426.Google Scholar
- Futuyma DJ, Moreno G. The evolution of ecological specialization. Annu Rev Ecol Evol Syst. 1998;19:207–33.View ArticleGoogle Scholar
- Devictor V, Clavel J, Julliard R, Lavergne S, Mouillot D, Thuiller W, Venail P, Villéger S, Mouquet N. Defining and measuring ecological specialization. J Appl Ecol. 2010;47:15–25.View ArticleGoogle Scholar
- Poisot T, Bever JD, Nemri A, Thrall PH, Hochberg ME. A conceptual framework for the evolution of ecological specialisation. Ecol Lett. 2011;14:841–51.PubMedPubMed CentralView ArticleGoogle Scholar
- Poisot T, Thrall PH, Hochberg ME. Trophic network structure emerges through antagonistic coevolution in temporally varying environments. Proc R Soc Lond B Biol Sci. 2012;279:299–308.View ArticleGoogle Scholar
- Dormann CF, Fruend J, Gruber B. Bipartite: Visualising Bipartite Networks and Calculating Some (Ecological) Indices; 2016. https://cran.rproject.org/web/packages/bipartite/index.html.
- Pagel MD. Inferring the historical patterns of biological evolution. Nature. 1999;401:877–84.PubMedView ArticleGoogle Scholar
- Revell LJ. Phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol Evol. 2011;3:217–23.View ArticleGoogle Scholar
- Felsenstein J. Phylogenies and the comparative method. Am Nat. 1985;125:1–15.View ArticleGoogle Scholar
- Harmon LJ, Weir J, Brock C, Glor RE, Challenger W. GEIGER: investigating evolutionary radiations. Bioinformatics. 2008;24:129–31.PubMedView ArticleGoogle Scholar
- Orme D. CAPER: Comparative Analyses of Phylogenetics and Evolution in R. ver. 0.52; 2013. http://cran.r-project.org/web/packages/caper.
- Midford PE, Garland T Jr., Maddison WP. PDAP Package of Mesquite. Version 1.07; 2005. http://mesquiteproject.org/pdap_mesquite/.
- Sidlauskas B. Continuous and arrested morphological diversification in sister clades of characiform fishes: a phylomorphospace approach. Evolution. 2008;62:3135–56.PubMedView ArticleGoogle Scholar
- Harmon LJ, Schulte JA II, Larson A, Losos JB. Tempo and mode of evolutionary radiation in iguanian lizards. Science. 2003;301:961–4.PubMedView ArticleGoogle Scholar
- Freckleton RP, Harvey PH. Detecting non-Brownian trait evolution in adaptive radiations. PLoS Biol. 2003;4:2104–11.Google Scholar
- Beaulieu JM, Jhwueng D-C, Boettiger C, O’Meara BC. Modeling stabilizing selection: expanding the Ornstein-Uhlenbeck model of adaptive evolution. Evolution. 2012;66:2369–83.PubMedView ArticleGoogle Scholar
- Butler MA, King AA. Phylogenetic comparative analysis: a modelling approach for adaptive evolution. Am Nat. 2004;164:683–95.View ArticleGoogle Scholar
- Berger D, Chobanov DP, Mayer F. Interglacial refugia and range shifts of the alpine grasshopper Stenobothrus cotticus (Orthoptera: Acrididae: Gomphocerinae). Org Divers Evol. 2010;10:123–33.View ArticleGoogle Scholar
- Gangwere SK, Morales-Agacino F. The biogeography of Iberian orthopteroids. Miscelánea Zoológica. 1990;2:9–75.Google Scholar
- Bridle JR, Saldamando CI, Koning W, Butlin RK. Assortative preferences and discrimination by females against hybrid male song in the grasshoppers Chorthippus brunneus and Chorthippus jacobsi (Orthoptera: Acrididae). J Evol Biol. 2006;19:1248–56.PubMedView ArticleGoogle Scholar
- Butlin RK, Beaumont M, Hewitt GM. Selection for assortative mating between parapatric subspecies of grasshopper. Anim Behav. 1992;43:1045–7.View ArticleGoogle Scholar
- Hopkins MJ, Smith AB. Dynamic evolutionary change in post-Paleozoic echinoids and the importance of scale when interpreting changes in rates of evolution. Proc Natl Acad Sci U S A. 2015;112:3758–63.PubMedPubMed CentralView ArticleGoogle Scholar
- Ingram T, Harmon LJ, Shurin JB. When should we expect early bursts of trait evolution in comparative data? Predictions from an evolutionary food web model. J Evol Biol. 2012;25:1902–10.PubMedView ArticleGoogle Scholar
- Wiens JJ, Ackerly DD, Allen AP, Anacker BL, Buckley LB, Cornell HV, Damschen EI, Davies TJ, Grytnes JA, Harrison SP, Hawkins BA, Holt RD, Mccain CM, Stephens PR. Niche conservatism as an emerging principle in ecology and conservation biology. Ecol Lett. 2010;13:1310–24.PubMedView ArticleGoogle Scholar
- Münkemüller T, Boucher F, Thuiller W, Lavergne S. Phylogenetic niche conservatism - common pitfalls and ways forward. Funct Ecol. 2015;29:627–39.PubMedPubMed CentralView ArticleGoogle Scholar
- Evans MEK, Smith SA, Flynn R, Donoghue MJ. Climate, niche evolution, and diversification of the “bird-cage” evening primroses (Oenothera, sections Anogra and Kleinia). Am Nat. 2009;173:225–40.PubMedView ArticleGoogle Scholar
- López-Fernández H, Arbour JH, Winemiller KO, Honeycutt RL. Testing for ancient adaptive radiations in Neotropical cichlid fishes. Evolution. 2013;67:1321–37.PubMedGoogle Scholar
- Toro E, Herrel A, Irschick D. The evolution of jumping performance in Caribbean Anolis lizards: solutions to biomechanical trade-offs. Am Nat. 2004;163:844–56.PubMedView ArticleGoogle Scholar
- Calsbeek R, Irschick DJ. The quick and the dead: correlational selection on morphology, performance, and habitat use in island lizards. Evolution. 2007;61:2493–503.PubMedView ArticleGoogle Scholar
- Irschick DJ, Meyers JJ, Husak JF, Le Galliard JF. How does selection operate on whole-organism functional performance capacities? A review and synthesis. Evol Ecol Res. 2008;10:177–96.Google Scholar
- Biewener A. Animal Locomotion. Oxford: Oxford University Press; 2003.Google Scholar
- Buckley TR, Attanayake D, Bradler S. Extreme convergence in stick insect evolution: phylogenetic placement of the Lord Howe Island tree lobster. Proc R Soc Lond B Biol Sci. 2009;276:1055–62.View ArticleGoogle Scholar
- Gilman CA, Irschick DJ. Foils of flexion: the effects of perch compliance on lizard locomotion and perch choice in the wild. Funct Ecol. 2013;27:374–81.View ArticleGoogle Scholar
- Burrows M. Anatomy of the hind legs and actions of their muscles during jumping in leafhopper insects. J Exp Biol. 2007;210:3590–600.PubMedView ArticleGoogle Scholar
- Endler JA. A predator's view of animal color patterns. Evol Biol. 1978;11:319–64.Google Scholar
- Chapman RF. Contact chemoreception in feeding by phytophagous insects. Annu Rev Entomol. 2003;48:455–84.PubMedView ArticleGoogle Scholar
- Losos JB, Mahler DL. Adaptive radiation: the interaction of ecological opportunity, adaptation, and speciation. In: Bell MA, Futuyma DJ, Eanes WF, Levinon JS, editors. Evolution since Darwin: the first 150 years. Sunderland: Sinauer Assoc; 2010.Google Scholar
- Blom MPK, Horner P, Moritz C. Convergence across a continent: adaptive diversification in a recent radiation of Australian lizards. Proc R Soc Lond B Biol Sci. 2016;283:20160181.View ArticleGoogle Scholar
- Friedman ST, Price SA, Hoey AS, Wainwright PC. Ecomorphological convergence in planktivorous surgeonfishes. J Evol Biol. 2016;29:965–78.PubMedView ArticleGoogle Scholar