Proceedings of the Royal Society B: Biological Sciences
You have accessResearch article

Multi-modal signal evolution in birds: re-examining a standard proxy for sexual selection

Christopher R. Cooney

Christopher R. Cooney

Department of Zoology, University of Oxford, South Parks Road, Oxford OX1 3PS, UK

Department of Animal and Plant Sciences, University of Sheffield, Western Bank, Sheffield S10 2TN, UK

[email protected]

Google Scholar

Find this author on PubMed

,
Hannah E. A. MacGregor

Hannah E. A. MacGregor

Department of Zoology, University of Oxford, South Parks Road, Oxford OX1 3PS, UK

School of Biological Sciences, University of Bristol, Bristol BS8 1TQ, UK

Google Scholar

Find this author on PubMed

,
Nathalie Seddon

Nathalie Seddon

Department of Zoology, University of Oxford, South Parks Road, Oxford OX1 3PS, UK

Google Scholar

Find this author on PubMed

and
Joseph A. Tobias

Joseph A. Tobias

Department of Zoology, University of Oxford, South Parks Road, Oxford OX1 3PS, UK

Department of Life Sciences, Imperial College London, Silwood Park, Buckhurst Road, Ascot, Berkshire SL5 7PY, UK

Google Scholar

Find this author on PubMed

    Abstract

    Sexual selection is proposed to be an important driver of speciation and phenotypic diversification in animal systems. However, previous phylogenetic tests have produced conflicting results, perhaps because they have focused on a single signalling modality (visual ornaments), whereas sexual selection may act on alternative signalling modalities (e.g. acoustic ornaments). Here, we compile phenotypic data from 259 avian sister species pairs to assess the relationship between visible plumage dichromatism—a standard index of sexual selection in birds—and macroevolutionary divergence in the other major avian signalling modality: song. We find evidence for a strong negative relationship between the degree of plumage dichromatism and divergence in song traits, which remains significant even when accounting for other key factors, including habitat type, ecological divergence and interspecific interactions. This negative relationship is opposite to the pattern expected by a straightforward interpretation of the sexual selection–diversification hypothesis, whereby higher levels of dichromatism indicating strong sexual selection should be related to greater levels of mating signal divergence regardless of signalling modality. Our findings imply a ‘trade-off’ between the elaboration of visual ornaments and the diversification of acoustic mating signals, and suggest that the effects of sexual selection on diversification can only be determined by considering multiple alternative signalling modalities.

    1. Introduction

    Previous studies have provided theoretical and empirical evidence that sexual selection can stimulate the rapid divergence of traits involved in mate choice and species recognition [1,2], supporting the longstanding view that sexual selection is an important driver of speciation and lineage diversification [36]. However, direct support for this hypothesis is relatively weak and inconsistent among taxa [7], with a series of studies finding no evidence of significant correlations between sexual selection and either species richness or speciation rate when studied across birds [811], mammals, butterflies and spiders [12], and certain fish taxa [13]. Although these observations suggest that sexual selection has limited effects on diversification at macroevolutionary scales, an alternative possibility is that standard comparative analyses are simply ineffective, because they rely on crude phenotypic proxies to estimate variation in sexual selection across species.

    To quantify the intensity of sexual selection, most existing large-scale studies in birds have used visible sex-differences in plumage coloration (e.g. [2,8,9,10,14,15]). This metric—usually termed ‘plumage dichromatism’—has become a standard proxy for sexual selection because it is easily measured and positively correlated with other indices of sexual selection such as testes size, the degree of polygyny and the frequency of extra-pair paternity [1619]. Nonetheless, the extent to which plumage dichromatism provides an accurate and consistent estimate of the overall intensity of sexual selection across all lineages remains uncertain, not least because it focuses on a single sexual signalling modality, whereas many taxa engage in multimodal signalling [5]. Indeed, if the intensity of sexual selection targeted at one signalling modality (e.g. visual signals) trades off or is negatively correlated with the intensity of sexual selection targeted at another (e.g. acoustic signals) [3,20], then such interactions could lead to a breakdown in the relationship between the underlying intensity of sexual selection across species and the visual traits used as proxies for sexual selection, therefore obscuring the true relationship between sexual selection and diversification [2,8].

    Progress in resolving this question has been slow because previous studies investigating the macroevolutionary consequences of sexual selection have generally focused exclusively on visual signalling traits [2,21], leaving open the possibility that comparisons across different sexual signalling modalities may reveal contrasting patterns. Furthermore, most studies have failed to address the role of other important selection pressures potentially shaping the evolution of signal phenotypes, such as habitat differences, ecological divergence and interspecific interactions [22,23], and have typically focused on geographically, taxonomically and/or ecologically restricted datasets, rather than sampling more broadly across major clades.

    Here, we address these issues by compiling data for a global sample of 259 avian species pairs from 33 passerine families to test the relationship between visible plumage dichromatism—used as a standard proxy for sexual selection in birds and other animals [7]—and macroevolutionary divergence in the other major avian signalling modality: song. We focus on birds because they offer an unequivocal example of multimodal sexual signalling in which both traits—avian plumage coloration (a visual ornament) and song (an acoustic ornament)—are known to function in inter- and intra-sexual selection in many avian taxa [5,2427]. In addition, the availability of complementary species-level data on avian morphological traits, ecology, biogeography and phylogeny allows us to assess the importance of plumage dichromatism in relation to a suite of key variables known to influence patterns of signal evolution.

    Our analyses can be divided into three stages. First, we use published song recordings to estimate the extent of song divergence within-species pairs. Second, we assess the relationship between sexual dichromatism and degree of song divergence across pairs. Third, we use multiple regression combined with model averaging techniques to assess the relative association between dichromatism and song divergence in relation to other factors. If sexual selection has reinforcing or independent effects on traits from different signalling modalities [28], we expect the relationship between plumage dichromatism and song divergence to be positive, or non-significant, respectively. Conversely, if the effects of sexual selection on traits in different signalling modalities are negatively correlated, we expect a negative relationship between plumage dichromatism and song divergence across species pairs.

    2. Methods

    (a) Species sampling and phylogenetic framework

    We used published molecular phylogenies to select a sample of passerine species pairs for which high-quality song recordings were available [2,8]. Each pair consisted of sister species (i.e. pairs of lineages that represent each other's closest relative). We note that a few of our study pairs contain species that are not true sisters, both because of incomplete sampling in published phylogenies, and because we included some near-sisters in which one member of the pair belonged to a sister clade (or both species from a polytomy). This approach is based on the assumption that comparisons between near-sisters are informative about phenotypic divergence during recent evolutionary history [2]. Overall, our sample contained 518 species from 259 species pairs (including 243 sister species and 13 near-sisters) widely distributed across the passerine radiation (electronic supplementary material, figure S1). For full details, see the electronic supplementary material. To provide a phylogenetic framework for our analyses, we sampled 1000 molecular-only trees from www.birdtree.org [29], which were pruned to include only the species included in our dataset. We then used TreeAnnotator [30] to generate a maximum clade credibility tree, which was then pruned so that each pair was represented by a single tip.

    (b) Song divergence

    To quantify the extent of song divergence within species pairs, we downloaded songs for all species from the Macaulay Library of Natural Sounds (www.macaulaylibrary.org) and the online database Xeno Canto (www.xeno-canto.org). We digitized sound files in Raven Pro v. 1.4 using standard settings, then measured seven key temporal and spectral traits that together capture important interspecific differences in overall signal structure (for full details, see [31,32]): (i) maximum frequency (kHz), (ii) minimum frequency (kHz), (iii) peak frequency (kHz; frequency in the signal with the greatest amplitude), (iv) bandwidth (kHz; maximum frequency minus minimum frequency), (v) signal duration (s), (vi) number of notes and (vii) pace (number of notes s−1). For each species, at least three high-quality recordings were measured (mean 4.8 recordings per species), providing a total sample of 2476 songs. To reduce the dimensionality of the dataset, we conducted a principal components (PC) analysis on the covariance matrix of individual (log-transformed) song measurements. The first three PCs from this analysis accounted for over 83% of the variance in the original acoustic dataset, with each PC capturing a distinct component of overall signal structure (electronic supplementary material, table S1). Specifically, PC1 (41%) loaded heavily with variables related to song pitch, PC2 (24%) loaded heavily with variables related to song length, and PC3 (18%) primarily loaded with song pace. We therefore interpreted these PCs as axes of variation in song pitch (PC1), length (PC2) and pace (PC3), respectively, with variation in overall song structure captured by position in this three-dimensional acoustic space (electronic supplementary material, figure S2). Using these PCs, we estimated within-pair song disparity for all pairs as the Euclidean distance between species' mean PC scores in terms of overall song structure (PC1–3), and in terms of song pitch, length and pace separately (figure 1). To assess the sensitivity of our results, we also generated an alternative version of our dataset in which within-species song disparity estimates were corrected for observed levels of intraspecific variation (see the electronic supplementary material, appendix S1 for full details).

    Figure 1.

    Figure 1. A phylogenetic tree of passerine species pairs (n = 259), showing within-pair disparity in overall song (PC1–3; innermost, black), pitch (PC2; red), length (PC3; blue) and pace (PC3; outermost, green). Size of points corresponds to relative within-pair song disparity. A version including species names is available in the electronic supplementary material.

    (c) Sexual dichromatism

    To quantify the degree of sexual dichromatism within pairs, we used published [8] species-level scores of dichromatism estimated by eye from handbook illustrations (see electronic supplementary material, appendix S1 for justification). Briefly, we used standard methodology [17,33] to score the difference in plumage coloration between the sexes over five body regions (head, nape-rump-back, throat-belly, tail and wings) for each species in our sample. Each region was scored separately using three scores: 0, no difference between the sexes; 1, difference between the sexes only in shade or intensity of colour; 2, difference in colour or pattern between the sexes. The dichromatism scores for all five body regions were then summed to give species-specific scores of plumage dichromatism on a scale from 0 (monochromatic) to 10 (maximum dichromatism).

    (d) Additional predictors of song divergence

    To explore the role of other factors known to influence estimates of phenotypic (particularly song) divergence in birds, we collected data for a suite of additional explanatory variables including divergence time [34], life history and allometric effects [35], migration status [36], habitat [37,38], breeding latitude and insularity [39,40], interspecific interactions [41], niche divergence [4244] and song learning [45]. Because the key habitat attribute linked to song evolution in birds is vegetation density [32] we used a score of forest dependency (i.e. degree of association with densely forested habitat). See the electronic supplementary material, appendix S1 for full details of methods, data and data sources.

    (e) Statistical analyses

    (i) Estimating phylogenetic signal of song divergence

    We used a generalized least-squares approach to test the phylogenetic signal of song divergence in our dataset. This approach, implemented in the R package caper [46], estimates a maximum-likelihood (ML) value for phylogenetic signal (λ) [47], which typically varies between zero (trait variance is independent of phylogeny) and one (trait variance follows a Brownian motion model of evolution). In the context of analysing song divergence, a value of λ = 0 indicates that extent of song divergence within pairs is random with respect to phylogeny, whereas a value of λ = 1 implies that closely related pairs have more similar levels of song disparity than would be expected by chance. We found that ML values of λ were zero for all four measures of song divergence, with values of λ = 1 (i.e. a Brownian motion model of evolution) significantly rejected in all cases (electronic supplementary material, table S2). Results were qualitatively similar for an alternative dataset corrected for observed levels of intraspecific variation (electronic supplementary material, table S2), indicating that variation in the extent of within-pair song divergence in our dataset is unrelated to phylogeny. This allowed us to use non-phylogenetic regression techniques with more flexible error structures than are currently possible in a statistical phylogenetic comparative framework, which was necessary for our dataset (see below).

    (ii) Testing the relationship between predictors and extent of song divergence

    To model the observed variation in estimates of within-pair song divergence, accounting for the right-skewed distribution of disparity estimates (electronic supplementary material, figure S3), we used generalized linear models (GLMs) with a gamma error distribution and log link. Using this approach, we (i) examined the relationship between song disparity and degree of plumage dichromatism, (ii) tested for an interaction effect between dichromatism and habitat type (forest dependency), and (iii) assessed the combined influence of all predictor variables on the extent of song disparity using single and multipredictor regression and Akaike information criterion-based model averaging [48] corrected for small sample sizes (AICc).

    To perform model averaging, following [49,50], we fitted models encompassing all possible additive combinations of our predictor variables (see above), including a null (intercept-only) model, calculating the AICc score of each model. We then calculated the relative importance (RI) for each predictor variable as the sum of relative Akaike weights for models in which they appear. RI values scale from 0 to 1, where a variable with a score of 0 is associated with very low Akaike weights (i.e. low importance) and 1 is consistently associated with high weights (i.e. high importance). We also calculated model-averaged estimates of regression parameters and standard error values, calculated as the sum of the parameter estimates for each model including that predictor, multiplied by the relative Akaike weight of each of those models. To give further insight into the relative importance of predictor variables, we also identified the variables included in the top-ranked (i.e. best-fitting) model in each case. We used this procedure to assess the effect of predictors on response variables, including overall song disparity (PC1–3), as well as separate estimates of disparity in pitch (PC1), length (PC2) and pace (PC3) separately. For multipredictor models, we restricted the dataset to include only those species pairs for which complete data for all predictors were available (246 of 259 pairs) and pseudo-R2 values for GLMs were estimated using the method of [51]. Pair age, generation length, body mass disparity and beak disparity were ln-transformed prior to analysis and models were inspected to ensure they complied with modelling assumptions (e.g. normality of residuals). We also checked for issues related to collinearity among predictors, which we found were unlikely to affect our results (see electronic supplementary material, appendix S1 for details). To improve the interpretability of regression coefficients, predictor variables were centered and standardized prior to model fitting [52]. All analyses were conducted in R version 3.3.1 and model averaging was performed using the R package MuMIn [53].

    3. Results

    (a) Relationship between dichromatism and extent of song divergence

    Our models revealed that plumage dichromatism was significantly negatively correlated with overall song divergence between species (electronic supplementary material, table S3). Species pairs with a greater degree of plumage dichromatism tended to have less divergent songs than more monochromatic species pairs (figure 2). Analysing patterns of divergence in each song trait separately revealed that the overall effect of dichromatism was primarily driven by significant negative relationships with divergence in song pitch (PC1) and length (PC2), with more marginal effects on song pace (PC3) (electronic supplementary material, table S3). Furthermore, including an interaction effect with forest dependency in these models revealed no significant statistical support for the hypothesis that the relationship between dichromatism and song divergence is mediated by variation in habitat type across taxa (electronic supplementary material, table S3). Rerunning analyses correcting for observed levels of intraspecific variation produced highly similar results (electronic supplementary material, table S3).

    Figure 2.

    Figure 2. Scatterplot showing the relationship between total within-pair song disparity (PC1–3) and plumage dichromatism across 259 species pairs of passerine birds. Regression line (with prediction intervals, shaded) indicates the best-fitting relationship between the two variables. (Online version in colour.)

    (b) Additional predictors of song divergence

    Single predictor regression models focused on our additional predictors of song divergence identified several variables that were individually correlated with variation in song disparity across pairs (electronic supplementary material, figures S4–S7). In terms of total song disparity (PC1–3), the strongest individual predictor was pair age (electronic supplementary material, table S4). Furthermore, variation in overall song disparity was also significantly correlated with disparity in beak morphology, with more marginal effects detected for several other variables, including forest dependency and mass disparity (electronic supplementary material, table S4). We also detected additional significant correlations between individual predictors and estimates of disparity in specific components of song structure (electronic supplementary material, table S5).

    We then assessed the relative importance of all predictors using AICc model averaging techniques (figure 3). The best-supported predictor of total song disparity (PC1–3) was pair age (RI = 0.99), which exhibited a strong and highly significantly positive relationship with disparity (electronic supplementary material, table S6). However, even after accounting for this relationship, the negative effect of dichromatism remained strong (RI = 0.98) (figure 3). The AICc best model for total song disparity accounted for 17% of the total variation, and retained these two variables plus mean pair body mass, forest dependency and within-pair beak disparity mass as positive effects, and confamilial sympatry as a negative effect (figure 3; electronic supplementary material, table S6).

    Figure 3.

    Figure 3. Model-averaged coefficient estimates from multipredictor GLMs predicting variation in within-pair song disparity (PC1–3) among passerine species pairs (n = 259). Points indicate the standardized effect sizes for each of the (scaled) predictor variables and lines indicate 95% confidence intervals. Sizes of points represent the relative importance (RI) of each of the predictor variables, where a value of RI = 0 indicates low importance and a value of RI = 1 indicates high importance. Predictors included in the AICc top model are coloured (blue, negative effect; red, positive effect), with significant (p < 0.05) model-averaged coefficients shown in darker colours.

    Analysing relationships in each component song trait separately revealed that predictor variables had contrasting effects depending on the axis considered (electronic supplementary material, figure S8 and table S7). On the one hand, divergence in song pitch (PC1) was best predicted by significant effects of dichromatism (RI = 0.97), pair age (RI = 0.96), body mass (RI = 0.93), confamilial sympatry (RI = 0.92) and forest dependency (RI = 0.89), whereas divergence in song length (PC2) was best explained by significant effects of only pair age (RI = 0.92) and dichromatism (RI = 0.91). By contrast, the only significant predictor of divergence in song pace (PC3) was a positive effect of within-pair beak disparity (RI = 0.82). Overall, AICc top models for these variables accounted for 22%, 12% and 6% of the total variation in disparity in song pitch, length and pace, respectively (electronic supplementary material, table S8). Rerunning models accounting for intraspecific variation produced qualitatively similar results (electronic supplementary material, tables S9–S11).

    4. Discussion

    Our analyses reveal that the degree of sexual dichromatism is negatively related to the extent of divergence in song structure among closely related bird species, a pattern that remained strong after accounting for a suite of potentially correlated or confounding variables, as well as for intraspecific trait variation. This contrasts with the findings of previous studies reporting evidence for significant positive correlations between indicators of sexual selection and signal divergence in birds, supporting the view that sexual selection can drive parallel divergence across multiple signals [2,21]. However, these studies assessed patterns of divergence in visual signalling traits (i.e. plumage coloration), using proxies for the intensity of sexual selection derived from the same signalling modality (e.g. dichromatism). By contrast, we have focused across major avian signalling modalities, finding the opposite relationship: that dichromatism (a visual signal) is negatively associated with divergence in song (an acoustic signal). Our results are therefore incompatible with a straightforward interpretation of the sexual selection–diversification hypothesis, whereby higher levels of dichromatism indicating strong sexual selection should be related to greater levels of mating signal divergence regardless of signalling modality. Instead, our findings are consistent with the alternative view that negative interactions between alternative signalling modalities play an important role in shaping macroevolutionary patterns of signal evolution in birds.

    One intuitive explanation for the negative correlation between plumage dichromatism and song divergence is that it reflects an underlying link between sexual selection and acoustic signal divergence in species that do not rely on visual signals. This makes sense because single-species studies have demonstrated an important role for female choice and/or male–male competition in shaping many aspects of avian acoustic signal design [24], and many avian taxa with drab or monochromatic plumage are known to possess highly elaborate acoustic signals which often provide the best means of differentiating among lineages (e.g. Old World leaf warblers; Phylloscopidae) [54]. Thus, increased sexual selection on acoustic traits relative to visual traits in monochromatic taxa provides a plausible explanation for a negative relationship between plumage dichromatism and song disparity at broad macroevolutionary scales.

    A key challenge facing this interpretation is to explain why, within species, selection would favour signals from one rather than multiple signalling modalities, thus generating negative relationships across modalities at a macroevolutionary scale. It is possible that the relative costs and benefits of signalling via a given sensory modality are shaped by the prevailing environmental conditions [38], such that ecological differences among species should play a role in determining the relative prominence of one signal type over another [27]. In line with this idea, bird species inhabiting dense habitats such as reedbeds, thickets and the understorey of forests often have more elaborate songs than visual signals. However, our models including forest dependency as an interaction term provided no support for the idea that the relationship between plumage dichromatism and song divergence is mediated by broad-scale habitat differences among taxa.

    An alternative explanation is that our findings reflect the signature of evolutionary trade-offs between alternative signalling modalities. Under a resource- or cost-based trade-off scenario—such as that envisaged by Darwin [3] and later termed the ‘transfer hypothesis’ [20]—constraints on sexual selection within species make it costly for males to signal in (or females to choose between) multiple signalling modalities [55,56], generating the potential for interspecific trade-offs in ornament elaboration (and diversification) between alternative signalling modalities [57]. This explanation relies on the assumption that investment in one signalling modality constrains investment in another, which is plausible given that avian plumage and song traits may both be costly to produce [58]. However, the energetic costs of signal production may be relatively low [59], and potentially offset by differences in how such signals are produced and displayed [27]. A different trade-off scenario is suggested by the concept of ‘redundancy’ among alternative signal types. Under a redundancy-based model, the spread of an attractive signal in one modality leads simultaneously to increased selection for detecting the novel signal and a weakening of selection for elaborate signals in alternative modalities, which occurs not because of costs associated with producing or maintaining multiple sexual signals, but because sexual selection on the latter trait is weak or non-existent, due to redundancy [60]. Such redundancy-based trade-offs can theoretically occur in the absence of any habitat differences among taxa, or resource limitation underlying the production of signalling traits. Thus, whether selection favours one signal type (e.g. song) over another (e.g. plumage) largely depends on which signal type evolved first, which may largely be due to historical contingencies [6163].

    (a) Contributory factors

    In addition to variation in the strength and targets of sexual selection, our results support a role for several other factors in shaping patterns of acoustic signal divergence in birds. We found strong evidence for a positive relationship between species pair age and degree of song disparity, in line with previous studies (e.g. [2,39,64,65]), as well as the general consensus that patterns of phenotypic divergence are primarily dictated by the time available for trait differences to evolve [34,66]. Body mass also emerged as a significant predictor of song divergence, in line with previous studies indicating positive relationships between body mass and patterns of signal evolution in birds [2,19]. Furthermore, we found support for links between song divergence and both habitat and the degree of confamilial sympatry. First, we found evidence for increased pitch disparity in species pairs with higher levels of forest dependency, consistent with the idea of stronger (divergent) selection on acoustic traits in taxa signalling in densely vegetated habitats [37,38]. Second, we found that pairs that co-occurred with a greater proportion of confamilial species had lower levels of song divergence than those with lower levels of overlap. This accords with the view that interactions among related species can constrain phenotypic divergence [67], in part because acoustic communities appear to ‘partition’ finite aspects of acoustic signalling space [31,6870]. Finally, we found that divergence in song pace was significantly positively correlated with disparity in species' beak morphology. Previous studies have found evidence of correlated evolution of morphology and vocal signal structure in particular clades (e.g. Darwin's finches, Neotropical woodcreepers) [43,71,72], presumably because biophysical constraints on song production generate correlated evolution between songs and beaks. Our results in relation to beak morphology support this view, and imply that this effect holds across passerines more generally. Nonetheless, even when we accounted for these significant effects in statistical models, the strong negative association between song divergence and dichromatism was retained.

    (b) Implications for comparative studies

    Our finding that dichromatism is negatively related to song divergence across a broad sample of avian species pairs has important implications for studies testing macroevolutionary hypotheses related to sexual selection. Most importantly, it implies that plumage dichromatism provides a relatively ineffective proxy for the intensity of sexual selection in taxa primarily using non-visual signals. This potential limitation of dichromatism has previously been proposed [2,8] with reference to bird species such as the common nightingale (Luscinia megarhynchos), common whitethroat (Sylvia communis) and sedge warbler (Acrocephalus schoenobaenus), passerine species with largely monomorphic plumage coloration, elaborate song traits and strong sexual selection [7375]. Our results provide broad-scale empirical support for this view, and indicate that dichromatism will underestimate sexual selection in these taxa, potentially being negatively related to the intensity of sexual selection in samples dominated by non-visual signallers. Thus, the underlying effect of sexual selection may often be obscured in comparative studies based solely on dichromatism, perhaps helping to explain the weak or non-existent correlations between dichromatism and speciation rates in birds and other taxa with multimodal signalling [711].

    (c) Conclusion

    Taken together, our findings are consistent with the view that sexual selection plays a major role in shaping sexual signal evolution, in conjunction with ecological factors [23,76,77]. However, whereas most previous studies have focused on a single signalling modality, we found evidence of a negative relationship between visual and acoustic signalling in birds, supporting the more general view that negative interactions between signalling modalities can explain general patterns of signal evolution [3,20]. Not only do these results suggest that such ‘trade-offs’ are important in shaping phenotypic diversity, they also indicate that phylogenetic tests based on phenotypic metrics for the intensity of sexual selection will underestimate the association between sexual selection and diversification. We conclude that the rigour and accuracy of any comparative analysis testing the effects of sexual selection will be improved by considering phenotypic proxies for sexual selection that span all relevant signalling modalities, be they visual, acoustic or olfactory.

    Data accessibility

    The full dataset has been uploaded to Dryad Data Repository and is available via the link Dryad Digital Repository: http://dx.doi.org/10.5061/dryad.b4p43t7 [78].

    Authors' contributions

    C.R.C., J.A.T. and N.S. developed the conceptual framework. C.R.C., J.A.T. and H.E.A.M. collected data. C.R.C. conducted the analyses and wrote the manuscript, with input from all authors.

    Competing interests

    We declare we have no competing interests.

    Funding

    Collection of trait data was supported by Natural Environment Research Grant (NE/I028068/1) to J.A.T. and N.S. This work was further supported by the European Research Council (grant no. 615709, Project ‘ToLERates’).

    Acknowledgements

    We thank Angela Chira, Emma Hughes, Joseph Llanos and Gavin Thomas for helpful discussion and three anonymous reviewers for constructive comments on the manuscript. We thank Mark Adams and Hein van Grouw at the Natural History Museum (Tring) for logistical assistance and access to specimens.

    Footnotes

    Electronic supplementary material is available online at https://dx.doi.org/10.6084/m9.figshare.c.4244924.

    Published by the Royal Society. All rights reserved.

    References