- Research article
- Open Access
How common is ecological speciation in plant-feeding insects? A 'Higher' Nematinae perspective
BMC Evolutionary Biologyvolume 10, Article number: 266 (2010)
Ecological speciation is a process in which a transiently resource-polymorphic species divides into two specialized sister lineages as a result of divergent selection pressures caused by the use of multiple niches or environments. Ecology-based speciation has been studied intensively in plant-feeding insects, in which both sympatric and allopatric shifts onto novel host plants could speed up diversification. However, while numerous examples of species pairs likely to have originated by resource shifts have been found, the overall importance of ecological speciation in relation to other, non-ecological speciation modes remains unknown. Here, we apply phylogenetic information on sawflies belonging to the 'Higher' Nematinae (Hymenoptera: Tenthredinidae) to infer the frequency of niche shifts in relation to speciation events.
Phylogenetic trees reconstructed on the basis of DNA sequence data show that the diversification of higher nematines has involved frequent shifts in larval feeding habits and in the use of plant taxa. However, the inferred number of resource shifts is considerably lower than the number of past speciation events, indicating that the majority of divergences have occurred by non-ecological allopatric speciation; based on a time-corrected analysis of sister species, we estimate that a maximum of c. 20% of lineage splits have been triggered by a change in resource use. In addition, we find that postspeciational changes in geographic distributions have led to broad sympatry in many species having identical host-plant ranges.
Our analysis indicates that the importance of niche shifts for the diversification of herbivorous insects is at present implicitly and explicitly overestimated. In the case of the Higher Nematinae, employing a time correction for sister-species comparisons lowered the proportion of apparent ecology-based speciation events from c. 50-60% to around 20%, but such corrections are still lacking in other herbivore groups. The observed convergent but asynchronous shifting among dominant northern plant taxa in many higher-nematine clades, in combination with the broad overlaps in the geographic distributions of numerous nematine species occupying near-identical niches, indicates that host-plant shifts and herbivore community assembly are largely unconstrained by direct or indirect competition among species. More phylogeny-based studies on connections between niche diversification and speciation are needed across many insect taxa, especially in groups that exhibit few host shifts in relation to speciation.
Ecological speciation is a process in which a shift in resource or habitat use within an ancestral species triggers the formation of two new sister species, each adapted to exploit different niches [1, 2]. The speciation process is thought to involve an initial period of resource polymorphism, during which the parent lineage utilizes multiple environments or niches [3, 4]; subsequently, tradeoffs in the efficiency by which individuals can use different resources lead to disruptive selection and, eventually, to lineage splitting [5, 6]. Ecological speciation has recently been focus of intensive research, and it is becoming increasingly clear that niche-based selection may underlie or at least speed up the diversification of, for example, many bird [7, 8], lizard , fish , and invertebrate [11, 12] groups.
Many of the best examples of species pairs that may have formed as a result of ecological speciation come from plant-feeding insects [2, 13, 14]. Insect herbivores are typically highly specialized in their host-plant preferences, and exhibit elaborate physiological , behavioural , and morphological  adaptations for utilizing their respective host plants. However, although only a small fraction of available plants constitute a suitable food source for most insect species, contrasts of herbivore versus plant phylogenies have in many cases revealed drastic discrepancies between the phylogenetic trees [18–20]. In essence, this means that host-plant associations are evolvable and change occasionally during the evolutionary history of insect lineages [21–23]. Occasional colonizations of novel hosts can theoretically cause ecological speciation, which could provide an explanation for the enormous species diversity of plant-feeding insects on the Earth [24–26].
A considerable proportion of evolutionary insect-plant research has been devoted to the possibility that ecological speciation in plant-feeding insects occurs in sympatry, so that both sister lineages are formed within a continuous geographical area [4, 13, 27]. However, a more likely scenario is that ecological speciation is initiated in allopatric or partially allopatric (para- or peripatric) settings; in these cases, increasing specialization onto different hosts in different parts of the geographical range of an insect species is thought to reduce the probability of hybridization if the populations later come into contact again [2, 28]. Hence, the main difference to 'ordinary' allopatric speciation is that the cause for reproductive isolation lies in the disparate ecology of the incipient species, rather than in gradual accumulation of genetic incompatibilities between geographically isolated populations [29, 30].
While numerous putative cases of ecology-based speciation in plant-feeding insects are known, we still lack an understanding of the actual frequency or importance of ecological shifts in the formation of new species [22, 27, 31]. Investigations of ecologically divergent species pairs can provide insights into traits or circumstances that enhance the likelihood of niche-based divergence, but studying the frequency of ecological speciation requires use of a broader, phylogeny-based approach [22, 30, 32, 33]. The usefulness of phylogenies stems from the fact that different speciational processes should produce very different distributions of niches on the phylogenetic trees of insect herbivores: if speciation is mainly allopatric and non-ecological, particular host taxa should be closely clustered on the tips of the insect phylogeny (Fig. 1a). Conversely, if speciation is mainly ecology-based, closely related insects should tend to have different hosts, meaning that plants would occur in an intermixed fashion along the tips of the insect phylogeny (Fig. 1b). In the first case, the number of inferred host-plant shifts should be distinctly lower than the number of speciation events that are required to produce the extant herbivore species, whereas, in the latter case, the number of inferred host shifts should be close to the number of lineage splits.
Here, we apply a molecular phylogenetic analysis of sawflies belonging to the so-called 'Higher' Nematinae (Hymenoptera: Tenthredinidae) to estimate the relative importance of ecological versus non-ecological speciation in plant-feeding insects. This group of over 700 species comprises most of the taxonomic and ecological diversity found within the tenthredinid subfamily Nematinae, which prompted Ross  to speculate: "The higher group of genera must have evolved some highly beneficial biological characteristics, because they are at present the most abundant boreal sawfly group in number of species and probably also in population." Indeed, higher nematines are ubiquitous in most habitats across the Northern Hemisphere, and their larvae feed on a wide variety of northern plant taxa [35–37]. Their larval feeding habits are equally diverse: in addition to 'normal' external folivores, the group includes berry miners, flower and catkin feeders, leaf folders, and various gall-inducing species (Fig. 2) . Combined with the high species number, such broad diversity in species-specific resource use presents many possibilities for studying the tempo and mode of speciation.
Taxon sampling, amplification, and sequencing
Our study builds on a previous phylogenetic analysis of the whole subfamily Nematinae  by adding 78 new species and sequences of a third gene (Cytochrome b) to the published dataset. The current taxon sample includes 127 exemplars of 125 Higher Nematinae species, meaning that nearly all higher-nematine species groups and main ecological niches (host-plant taxa and larval habits) are represented [35, 36]. Multiple representatives were included for all large genera and species groups (Additional file 1). Trees were rooted by including three non-nematine tenthredinids and ten species belonging to the nematine tribes Hoplocampini, Stauronematini, Pseudodineurini, Caulocampini, Susanini, Dineurini, and Cladiini as outgroups in the analyses. These small 'Lower' Nematinae groups form a paraphyletic grade with respect to the ingroup .
Sequence data were collected from two mitochondrial genes (Cytochrome oxidase I [CoI]: 810 bp; Cytochrome b [Cytb]: 718 bp) and from two exons (501 bp + 276 bp = 777 bp) of the F2 copy of the nuclear Elongation factor-1α (EF-1α) gene following previously-described protocols [39, 40]. The concatenated data matrix consists of 2305 bp of sequence data for 140 species. Sequences are missing for three, nine, and six species for CoI, Cytb, and EF-1α, respectively, but every included species has full-length sequences from at least two genes. New sequences have been submitted to GenBank under accession numbers HM237366-HM237589, and the Nexus-formatted data matrix, together with resultant phylogenetic trees, is available as Additional file 2.
Modeltest 3.5  was implemented in conjunction with PAUP* 4.0b10  to identify the least complex substitution model for use in Bayesian phylogenetic analyses in MrBayes 3.1.2 . Hierarchical likelihood ratio tests indicated a GTR+I+Γ4 model as optimal for each of the three genes. A separate, unlinked substitution model was allowed for each gene in a three-partition analysis. A single run employing default priors was run for eight million generations with eight incrementally heated (t = 0.1) chains; tree sampling was done from the current cold chain every 100th generation, and the first 10,001 trees recovered prior to reaching stationarity were discarded as a burnin. The consensus tree showing all compatible groupings (Fig. 3) was calculated on the basis of the remaining 70,000 trees. A corresponding maximum-likelihood (ML) analysis was performed using RAxML 7.0.4 . This analysis employed a separate GTR+I+Γ4 model for each gene, but branch lengths were estimated jointly for the whole data (Additional file 2). Clade support was estimated on the basis of 500 bootstrap replicates of the data matrix (Fig. 3).
BEAST 1.4.8  was used to estimate the relative ages of various nematine groups based on a Bayesian relaxed molecular clock method. The topologically unconstrained analysis allowed a separate GTR+I+ Γ4 model of substitution for each gene and employed an uncorrelated relaxed lognormal clock model for rate variation among branches, a Yule prior on speciation, and default priors for other parameters except for the mean of branch rates (ucld.mean), which was fixed to 1. Three independent runs with automatic tuning of operators were run for 80 million generations, and parameters and trees were sampled every 1,000 generations (the XML file is available as Additional file 3). After inspection of adequate convergence of runs and effective sample sizes of the parameters in Tracer 1.4.1 , the tree files were combined in LogCombiner 1.4.8 (part of the BEAST package). The first 40,000 trees from each file were discarded as a burnin, and the tree file was subsequently thinned by resampling trees every 3,000 generations; the maximum clade credibility (MCC) tree showing mean branch lengths (Fig. 4) is based on the 40,001 post-stationarity trees that remained after thinning.
To reconstruct ancestral host-plant families and feeding habits, these traits were treated as unordered multistate characters and maximum-parsimony optimized on the phylogenetic trees using Mesquite 2.6 . Oligo- and polyphagous taxa were coded with all used host families.
To estimate the number of ecological shifts that have occurred during the radiation of the Higher Nematinae, we first identified all distinct ecological niches (feeding habit × host plant(s)) found in the ingroup species included in the phylogenetic analysis, and coded each niche with a separate state within a single character (outgroup states were coded as unknown). Because the aim was to calculate numbers of changes, the typical number of steps between two different states was 1. However, we also created 'generalist' states for species that utilize multiple plant taxa and then used the step-matrix option in Mesquite to define the cost to these states, from the plant taxa that are included within the generalist host range, as being zero. By doing so, we essentially assumed that a clear overlap in the host ranges of different species implies that they have not speciated ecologically (theoretical models of resource-based speciation typically assume distinct, non-overlapping niches as the cause of divergent selection [2, 30]). Phylogenetic uncertainty in the estimate was taken into account by recording the numbers of steps in the niche character across the 70,000 post-burnin trees that were sampled by MrBayes during the phylogenetic analysis .
As a separate estimate of the proportion of lineage splits accompanied by a shift in resource use, we indentified all terminal sister-species pairs across the MCC tree (Fig. 4), and then separated these 35 pairs into those in which both species have identical or overlapping niches, and into those in which the species have different niches. Thereafter, we performed a logistic regression in SPSS for Windows 17.0 (SPSS, Inc., 233 S. Wacker Drive, Chicago, IL 60606-6307, USA) to test whether the probability that sister species have a different niche depends on the time elapsed since their most recent common ancestor (= relative node height in the MCC tree).
Proportions of higher nematine species feeding on different plant genera (Fig. 5) were extracted from Lacourt's  list of host-plant affiliations of sawflies of the Western Palearctic region. Only species with known hosts were included, and proportions were calculated separately for the tribe Pristiphorini and for the Nematini+Mesoneurini clade (see Figs. 3, 4 and 5). Oligo- and polyphagous species were counted as an additional species for each plant genus on which they feed (for example, the oligophagous Craesus latipes (Villaret) was treated as one species on Alnus and another on Betula).
The Bayesian and ML analyses of the sequence data produced relatively well-supported trees showing that the Higher Nematinae constitutes a monophyletic clade within the subfamily Nematinae (Figs. 3 and 4; the ML tree is included in Additional file 2). The three topologies are largely congruent, with discrepancies mainly evident in tree regions that are weakly supported. Major divisions within the ingroup correspond closely with the traditional tribes Nematini, Pristiphorini, and Mesoneurini (Figs. 3 and 4). Conflict with traditional classifications is mostly evident in that the largest nematine genera (Nematus, Pristiphora, and Amauronematus) come out as para- and polyphyletic (see ). Posterior probabilities and bootstrap proportions of many groupings within Nematini are low, but because these uncertainties concern mainly relationships among strongly supported middle-level clades, they have only minor importance for the conclusions below.
The trees document a pattern of frequent faunal exchange across the Holarctic region, because European, North American, and Asian exemplar species (see Additional file 1) are intermixed throughout the trees. Nearctic species are found scattered among European ones in all major higher-nematine genera (Nematus, Amauronematus, and Pristiphora), but also in smaller groups such as Eitelius, Pontopristia, and Pikonema (Figs. 3 and 4). Our analysis also confirms Smith's  hypothesis of a close relationship between the North American Nematus erythrogaster -group and the Holarctic genus Craesus, as well as his earlier  suggestion of a connection between the exclusively Nearctic genus Neopareophora and European Nepionema.
Evolutionary dynamics of resource use
Diversification of the Higher Nematinae has been a dynamic process in which host-plant associations and larval lifestyles change continually, although feeding-habit changes are distinctly rarer than shifts in host-plant use (Figs. 3 and 4). Various forms of internal feeding, such as gall induction, leaf folding, catkin feeding, and berry mining, have evolved repeatedly from ancestors whose larvae were external feeders on leaves or on needles (Fig. 4). Higher-nematine larvae are currently found on plants belonging to over 16 families, but most species are concentrated on plants in Salicaceae, Betulaceae, Rosaceae, Ericaceae, and Pinaceae (Fig. 5), partly because colonization and recolonization events have occurred repeatedly among these plant taxa (Fig. 3). Interestingly, the distribution of species across utilized plant families differs between the clade formed by Pristiphorini and its sister clade formed by Nematini+Mesoneurini (Fig. 5; χ2= 155.93, df = 15, P < 0.0001), the clearest discrepancy being the high proportion of Salix-feeding species within the most species-rich tribe Nematini.
Speciation and niche shifts
The existence of the 125 sampled ingroup species demands 124 past speciation events, but explaining the current distribution of species-level niches on the Bayesian consensus tree (Fig. 3) and on the ML tree requires only 68 shifts in feeding habits and/or host taxa. This estimate is robust against phylogenetic uncertainty, because when niches are similarly maximum-parsimony optimized on the 70,000 Bayesian post-burnin trees, 67-69 shifts (mean = 68.13) are needed. Even if shifts to generalist host-use states are treated as a true niche shift (= 1 step from other states), the Bayesian consensus tree is only 75 steps long (ML tree = 74 steps), and all trees in the Bayesian tree sample require between 73 and 75 changes (mean = 74.73 steps).
When only sister-taxon pairs on the MCC tree are considered, species in 19 out of 35 pairs (54.3%) have non-overlapping host ranges and/or a distinct difference in their larval feeding modes. However, the logistic regression (Fig. 6) shows that the probability of sister species having different niches is strongly affected by the time since their most recent common ancestor: Pdifference = 1/(1+e-(-1.29+15.12*split age)), the constant (P = 0.045) and the effect of split age (P = 0.013) being statistically significant.
Research on ecological speciation has traditionally focussed on sister-species pairs or on small groups of closely related lineages that differ in their resource use, and which could therefore have originated by niche shifts. While such studies have convincingly shown that ecology-based diversification is possible in highly disparate taxa and under many plausible scenarios (e.g., [8, 51, 52]), the overall frequency of ecological speciation remains unknown. Recently, broader phylogenetic approaches have provided insights into the relative importance of alternative speciation modes [53–55]. For example, comparative analyses employing age-range corrections have shown that closely related species tend to have less overlap in their geographical ranges than do distantly related species, which indicates that speciation rarely occurs in sympatric settings [53, 56–58]. However, the finding that speciation is largely allopatric does not exclude the possibility that the build-up of reproductive isolation between incipient species has an ecological basis: as mentioned above, ecological divergence can, and in fact is more likely to, occur in complete or partial allopatry [27, 28, 30]. Therefore, phylogenetic studies on the frequency of niche shifts in relation to the number of past speciation events are more likely to produce a correct view of the prevalence of ecological speciation [22, 31, 33, 56].
The traditional paradigm in plant-herbivore research is that host shifts are a major factor promoting species divergence [23, 59, 60]. However, there are two good reasons for suspecting that the importance of niche shifts for speciation is overestimated: First, there has been a huge--most likely disproportionate--interest in the intriguing possibility that host-associated speciation in insect herbivores could occur in sympatry, i.e., without geographical isolation [4, 13, 27, 58]. Second, it seems probable--and perfectly logical--that insect groups that are chosen for phylogenetic studies on host-plant shifts are selected preferentially from taxa in which species are known a priori to be relatively specialized and to exhibit clear interspecific variation in host-plant use (e.g., [19, 40, 61]).
In the case of the Higher Nematinae, our phylogenetic analysis reveals a pattern of frequent niche shifts both in terms of larval lifestyles and host-plant use (Figs. 3 and 4). Despite this, optimizing larval niches on the Bayesian tree sample shows that at most about 60% of lineage splits could have been caused by ecological factors. This value should be considered as an upper limit, because our ecologically overdispersed taxon-sampling scheme will raise the relative number of niche shifts, and the sampling design should also override the tendency of maximum parsimony to underestimate the frequency of changes in fast-evolving traits . When only sister species are considered, the percentage of pairs having divergent niches is 54. This raw value is intriguingly close to Winkler & Mitter's  recent 'fifty-fifty' estimate of the proportions of ecological vs. nonecological speciation, which was based on a broad literature survey of sister species of herbivorous insects. However, our logistic regression (Fig. 6) shows that immediately after speciation only an estimated 21.6% of higher-nematine sister-species pairs would have non-overlapping niches. The discrepancy between the methods is most likely explained by postspeciational host shifts, which can inflate the apparent frequency of ecological speciation in uncorrected sister-species comparisons. Although denser taxon sampling will be needed for a more exact estimate of the prevalence of ecology-based diversification within the Higher Nematinae, the marked drop in the inferred proportion in the age-adjusted analysis suggests that phylogenetic time corrections would be useful also in surveys of other insect herbivore taxa.
It is possible, however, that the frequency of ecological speciation varies among clades [22, 31]. In particular, it appears that the extreme diversity (400-500 species [63, 64]) of gall-inducing nematines in the subtribe Euurina (Figs. 3 and 4) has been spurred by host-shifting among Salix species. Like other gall-inducing insects , Euurina gallers are very host-specific compared to willow-associated nematines having external-feeding larvae, which tend to utilize multiple host species [36, 66], and which therefore probably may have radiated mainly allopatrically. Higher-nematine subgenera and species-groups feeding externally on other plant taxa are likewise often dominated by species having identical or broadly overlapping host-plant ranges [36, 38, 66], and niche shifting seems to be particularly infrequent in relation to speciation in groups associated with Picea, Larix, Vaccinium, and plants within Betulaceae (Fig. 3). It remains to be studied whether such non-ecological radiations can be used to estimate the proportion of non-ecological speciation in related groups in which species differ in their host use, because some proportion of host switches also in these groups undoubtedly have occurred well before or after speciation events (cf. ).
While earlier hypotheses on insect diversification emphasized coevolution of plant defenses and herbivore counterdefenses as a major driver of insect diversification [59, 67, 68], recent studies applying dated phylogenies have uncovered a possible role of long-term climatic conditions in determining rates of speciation and extinction in various herbivore groups [69–71]. Throughout the Earth's history, climatic changes have lead to major shifts in plant communities and global vegetation patterns, with direct negative or positive consequences for associated herbivores . In particular, systems experiencing repeated cycles of range contractions, expansions, and faunal mixing can constitute 'speciation machines' that lead to escalating diversification across multiple trophic levels [73–75]. A Cenozoic high-latitude speciation machine could explain why higher nematines seem to have an inordinate fondness for willows: nearly a third of Western Palearctic species in the tribe Pristiphorini feed on Salix species, and in the Nematini the percentage is as high as 68 (Fig. 5; this general pattern holds also in North America [35, 37]). The largest willow-associated radiations have occurred within the aforementioned gall-inducing subtribe Euurina, and in the predominantly willow-feeding Amauronematus+Pontopristia clade, which includes at least 112 species [36, 76]. These temporally overlapping radiations began c. 30 million years ago (Mya), assuming that the most recent common ancestor of the Higher Nematinae lived about 70 Mya (Fig. 4; see ). Interestingly, this would place the onset of these radiations at the time of strong cooling of the global climate, which began in the early Oligocene c. 35 Mya, and which was followed by alternating periods of cold ice ages and warmer interglacials . Such climatic oscillations, with resultant long-distance migrations of whole ecosystems, could have promoted the diversification of willows that are concentrated in relatively cool habitats and that currently comprise over 400 species [78, 79]. The conditions that generated diversity in Salix would simultaneously have acted also on the insect groups that depend on them and, like in higher nematines, willows currently support a considerable proportion of species also in, for example, northern butterflies and moths , phytophagous beetles [80, 81], and leafhoppers .
The role of competition in directing the historical assembly and present structure of herbivore communities has been debated for decades [82–84]. If competition was a force directing host switching, shifts would tend to occur towards un- or underused plant taxa, meaning that, over time, herbivore host-plant associations would become overdispersed with regard to plant phylogeny. By contrast, higher-nematine host use is strongly underdispersed, shifts having occurred repeatedly and in many directions among a handful of dominant northern plant families, while a large proportion of the Holarctic flora apparently has been effectively ignored for tens of millions of years. This shifting pattern conforms to the 'resource island model' [23, 85, 86] of herbivore diversification, in which phylogenetically biased colonizations and back-colonizations among plant taxa, in combination with abundance-dependent extinction, lead to accumulation of herbivore species on common plants that have many relatives [80, 87, 88]. Competition could still operate more subtly, if recruitment follows a 'macroevolutionary ideal free distribution' (cf. ), so that the number of herbivore species that can be supported depends on the commonness (or overall biomass) of a given plant . However, the convergent, asynchronous, and undoubtedly ongoing colonizations of many plant taxa by various higher-nematine groups (Figs. 3 and 4) indicates that ecological pre-emption of host taxa does not occur, and that northern insect-plant communities are still unsaturated and could therefore soak up even more herbivore species in the future. The broad overlaps in the geographical distributions [35, 36] of many closely related, ecologically near-identical higher-nematine species--that necessarily must have diverged in allopatry and then brought to sympatry by postspeciational range shifts--provides further support for the view that interspecific competition, either via direct resource competition or via indirect competition caused by shared natural enemies, is of minor importance in structuring herbivore communities [73, 82, 84, 91].
Our phylogeny-based analysis of the Higher Nematinae strongly indicates that the importance of niche shifts for speciation in plant-feeding insects is at present explicitly and implicitly overestimated. In particular, applying a time correction for sister-group comparisons lowered the proportion of apparent ecology-based speciation events from roughly 50% to around 20%. The vast majority of lineage splits in higher nematines therefore seem to have occurred non-ecologically in allopatry, and this may well be true also for most other plant-feeding insects. Reconciling this result with the finding of Janz et al.  that species richness in nymphalid butterfly clades correlates positively with collective host ranges requires further work; we propose that the correlation follows from reduced extinction probabilities in ecologically versatile groups, rather than from increased ecological speciation within them.
Evolutionary dynamics observed within the Higher Nematinae favour a largely non-interactive, non-equilibrium view of community assembly in northern plant-herbivore networks: geographical shifts across the whole Northern Hemisphere have been commonplace in many higher-nematine groups, and the frequent co-occurrence of related species utilizing seemingly identical niches indicates that distributional changes occur largely unimpeded by direct or indirect competitive interactions. More detailed surveys of local communities are, however, necessary in order to exclude the possibility that competitive repulsion occurs on smaller geographical scales, which could lead to a mosaic pattern of patch occupancy by ecologically equivalent relatives .
The Higher Nematinae comprises well over 70% of the species in the subfamily Nematinae, but the main part of higher-nematine diversity lies within the tribe Nematini, in which a strikingly high proportion of species use willows as hosts. This suggests that the success of higher nematines was caused, not by the evolution of superior biological characteristics as suggested by Ross , but by a fortuitous association with willows at a time of a cyclically cooling global climate. Reliably dated molecular-phylogenetic analyses of Salix and Salix-associated herbivores are desperately needed to test our hypothesis that the diversification of willow-based food webs was accelerated during the latter half of the Cenozoic Era. The genus Salix has thus far proven to be an extremely challenging target for such studies [94, 95], but even comparative analyses across herbivore taxa would surely provide interesting insights into the evolutionary history of Holarctic plant-herbivore communities.
Schluter D: The Ecology of Adaptive Radiation. 2000, Oxford, Oxford University Press
Rundle HD, Nosil P: Ecological speciation. Ecol Lett. 2005, 8: 336-352. 10.1111/j.1461-0248.2004.00715.x.
Smith TB, Skúlason S: Evolutionary significance of resource polymorphisms in fishes, amphibians, and birds. Annu Rev Ecol Syst. 1996, 27: 111-133. 10.1146/annurev.ecolsys.27.1.111.
Drès M, Mallet J: Host races in plant-feeding insects and their importance in sympatric speciation. Phil Trans R Soc Lond B. 2002, 357: 471-492. 10.1098/rstb.2002.1059.
Räsänen K, Hendry AP: Disentangling interactions between adaptive divergence and gene flow when ecology drives diversification. Ecol Lett. 2008, 11: 624-636. 10.1111/j.1461-0248.2008.01176.x.
Nosil P, Harmon LJ, Seehausen O: Ecological explanations for (incomplete) speciation. Trends Ecol Evol. 2009, 24: 145-156. 10.1016/j.tree.2008.10.011.
Ryan PG, Bloomer P, Moloney CL, Grant TJ, Delport W: Ecological speciation in South Atlantic island finches. Science. 2007, 315: 1420-1423. 10.1126/science.1138829.
Irestedt M, Fjeldså J, Dalén L, Ericson PGP: Convergent evolution, habitat shifts and variable diversification rates in the ovenbird-woodcreeper family (Furnariidae). BMC Evol Biol. 2009, 9: 268-10.1186/1471-2148-9-268.
Thorpe RS, Surget-Groba Y, Johansson H: The relative importance of ecology and geographic isolation for speciation in anoles. Phil Trans R Soc B. 2008, 363: 3071-3081. 10.1098/rstb.2008.0077.
Schluter D, Conte GL: Genetics and ecological speciation. Proc Natl Acad Sci USA. 2009, 106: 9955-9962. 10.1073/pnas.0901264106.
Parent CE, Crespi BJ: Sequential colonization and diversification of Galápagos endemic land snail genus Bulimulus (Gastropoda: Stylommatophora). Evolution. 2006, 60: 2311-2328.
Nyman T, Bokma F, Kopelke J-P: Reciprocal diversification in a complex plant-herbivore-parasitoid food web. BMC Biol. 2007, 5: 49-10.1186/1741-7007-5-49.
Via S: Sympatric speciation in animals: the ugly duckling grows up. Trends Ecol Evol. 2001, 16: 381-390. 10.1016/S0169-5347(01)02188-7.
Peccoud J, Ollivier A, Plantegenest M, Simon J-C: A continuum of genetic divergence from sympatric host races to species in the pea aphid complex. Proc Natl Acad Sci USA. 2009, 106: 7495-7500. 10.1073/pnas.0811117106.
Després L, David J-P, Gallet C: The evolutionary ecology of insect resistance to plant chemicals. Trends Ecol Evol. 2007, 22: 298-307. 10.1016/j.tree.2007.02.010.
Karban R, Agrawal AA: Herbivore offense. Annu Rev Ecol Syst. 2002, 33: 641-664. 10.1146/annurev.ecolsys.33.010802.150443.
Toju H: Natural selection drives the fine-scale divergence of a coevolutionary arms race involving a long-mouthed weevil and its obligate host plant. BMC Evol Biol. 2009, 9: 273-10.1186/1471-2148-9-273.
Percy DM, Page RDM, Cronk QCB: Plant-insect interactions: double-dating associated insect and plant lineages reveals asynchronous radiations. Syst Biol. 2004, 53: 120-127. 10.1080/10635150490264996.
Lopez-Vaamonde C, Wikström N, Labandeira C, Godfray HCJ, Goodman SJ, Cook JM: Fossil-calibrated molecular phylogenies reveal that leaf-mining moths radiated several million years after their host plants. J Evol Biol. 2006, 19: 1314-1326. 10.1111/j.1420-9101.2005.01070.x.
Hunt T, Bergsten J, Levkanicova Z, Papadopoulou A, St. John O, Wild R, Hammond PM, Ahrens D, Balke M, Caterino MS, Gómez-Zurita H, Ribera I, Barraclough TG, Bocakova M, Bocak L, Vogler AP: A comprehensive phylogeny of beetles reveals the evolutionary origins of a superradiation. Science. 2007, 318: 1913-1916. 10.1126/science.1146954.
Janz N, Nylin S: The oscillation hypothesis of host plant-range and speciation. Specialization, Speciation, and Radiation: the Evolutionary Biology of Herbivorous Insects. Edited by: Tilmon KJ. 2008, Berkeley, University of California Press, 203-215.
Winkler IS, Mitter C: The phylogenetic dimension of insect-plant interactions: a review of recent evidence. Specialization, Speciation, and Radiation: the Evolutionary Biology of Herbivorous Insects. Edited by: Tilmon KJ. 2008, Berkeley, University of California Press, 240-263.
Nyman T: To speciate, or not to speciate? Resource heterogeneity, the subjectivity of similarity, and the macroevolutionary consequences of niche-width shifts in plant-feeding insects. Biol Rev. 2010, 85: 393-411. 10.1111/j.1469-185X.2009.00109.x.
Mitter C, Farrell B, Wiegmann B: The phylogenetic study of adaptive zones: has phytophagy promoted insect diversification?. Am Nat. 1988, 132: 107-128. 10.1086/284840.
Dyer LA, Singer MS, Lill JT, Stireman JO, Gentry GL, Marquis RJ, Ricklefs RE, Greeney HF, Wagner DL, Morais HC, Diniz IR, Kursar TA, Coley PD: Host specificity of Lepidoptera in tropical and temperate forests. Nature. 2007, 448: 696-700. 10.1038/nature05884.
Futuyma DJ, Agrawal AA: Macroevolution and the biological diversity of plants and herbivores. Proc Natl Acad. Sci USA. 2009, 106: 18054-18061. 10.1073/pnas.0904106106.
Futuyma DJ: Sympatric speciation: norm or exception?. Specialization, Speciation, and Radiation: the Evolutionary Biology of Herbivorous Insects. Edited by: Tilmon KJ. 2008, Berkeley, University of California Press, 136-148.
Funk DJ: Isolating a role for natural selection in speciation: host adaptation and sexual isolation in Neochlamisus bebbianae leaf beetles. Evolution. 1998, 52: 1744-1759. 10.2307/2411347.
Funk DJ, Nosil P, Etges WJ: Ecological divergence exhibits consistently positive associations with reproductive isolation across disparate taxa. Proc Natl Acad Sci USA. 2006, 103: 3209-3213. 10.1073/pnas.0508653103.
Rundell RJ, Price TD: Adaptive radiation, nonadaptive radiation, ecological speciation and nonecological speciation. Trends Ecol Evol. 2009, 24: 394-399. 10.1016/j.tree.2009.02.007.
Ross HH: An uncertainty principle in ecological evolution. Univ Arkansas Mus Occ Pap. 1972, 4: 133-157.
Barraclough TG, Hogan JE, Vogler AP: Testing whether ecological factors promote cladogenesis in a group of tiger beetles (Coleoptera: Cicindelidae). Proc R Soc Lond B. 1999, 266: 1061-1067. 10.1098/rspb.1999.0744.
Peterson AT, Soberón J, Sánchez-Cordero V: Conservatism of ecological niches in evolutionary time. Science. 1999, 285: 1265-1267. 10.1126/science.285.5431.1265.
Ross HH: A generic classification of the Nearctic sawflies. Ill Biol Monogr. 1937, 34: 1-173.
Smith DR: Suborder Symphyta. Catalog of Hymenoptera in America North of Mexico. Edited by: Krombein KV, Hurd PD Jr, Smith DR, Burks BD. 1979, Washington, Smithsonian Institution Press, 3-137.
Lacourt J: Répertoire des Tenthredinidae ouest-paléarctiques. Mém SEF. 1999, 3: 1-432.
Kouki J: Latitudinal gradients in species richness in northern areas: some exceptional patterns. Ecol Bull. 1999, 47: 30-37.
Nyman T, Farrell BD, Zinovjev AG, Vikberg V: Larval habits, host-plant associations, and speciation in nematine sawflies (Hymenoptera: Tenthredinidae). Evolution. 2006, 60: 1622-1637.
Nyman T, Zinovjev AG, Vikberg V, Farrell BD: Molecular phylogeny of the sawfly subfamily Nematinae (Hymenoptera: Tenthredinidae). Syst Entomol. 2006, 31: 569-583. 10.1111/j.1365-3113.2006.00336.x.
Nyman T, Widmer A, Roininen H: Evolution of gall morphology and host-plant relationships in willow-feeding sawflies (Hymenoptera: Tenthredinidae). Evolution. 2000, 54: 526-533.
Posada D, Crandall KA: Modeltest: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.
Swofford DL: PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods), version 40b10. 2002, Sunderland, Sinauer
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.
Stamatakis A, Hoover P, Rougemont J: A rapid bootstrap algorithm for the RAxML web-servers. Syst Biol. 2008, 75: 758-771. 10.1080/10635150802429642.
Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.
Rambaut A, Drummond AJ: Tracer, version 1.4.1. 2007, [http://beast.bio.ed.ac.uk/Tracer]
Maddison WP, Maddison DR: Mesquite: a modular system for evolutionary analysis, version 2.6. 2009, [http://mesquiteproject.org]
Huelsenbeck JP, Rannala B, Masly JP: Accommodating phylogenetic uncertainty in phylogenetic studies. Science. 2000, 288: 2349-2350. 10.1126/science.288.5475.2349.
Smith DR: The abbotii and erythrogaster groups of Nematus Panzer (Hymenoptera: Tenthredinidae) in North America. Proc Entomol Soc Wash. 2008, 110: 647-667.
Smith DR: Nepionema, a nematine sawfly genus new to North America, and an unusual new species of Nematus (Hymenoptera: Tenthredinidae). Proc Entomol Soc Wash. 1994, 96: 133-138.
Savolainen V, Anstett M-C, Lexer C, Hutton I, Clarkson JJ, Norup MV, Powell MP, Springate D, Salamin N, Baker WJ: Sympatric speciation in palms on an oceanic island. Nature. 2006, 441: 210-213. 10.1038/nature04566.
Kempf F, Boulinier T, De Meeûs T, Arnathau C, McCoy KD: Recent evolution of host-associated divergence in the seabird tick Ixodes uriae. Mol Ecol. 2009, 18: 4450-4462. 10.1111/j.1365-294X.2009.04356.x.
Barraclough TG, Vogler AP: Detecting the geographical pattern of speciation from species-level phylogenies. Am Nat. 2000, 155: 419-434. 10.1086/303332.
Barraclough TG, Nee S: Phylogenetics and speciation. Trends Ecol Evol. 2001, 16: 391-399. 10.1016/S0169-5347(01)02161-9.
Wood TE, Takebayashi N, Barker MS, Mayrose I, Greenspoon PB, Riesenberg LH: The frequency of polyploid speciation in vascular plants. Proc Natl Acad Sci USA. 2009, 106: 13875-13879. 10.1073/pnas.0811575106.
Kozak KH, Wiens JJ: Does niche conservatism promote speciation? A case study in North American salamanders. Evolution. 2006, 60: 2604-2621.
Perret M, Chautems A, Spichiger R, Barraclough TG, Savolainen V: The geographical pattern of speciation and floral diversification in the Neotropics: the tribe Sinningieae (Gesneriaceae) as a case study. Evolution. 2007, 61: 1641-1660. 10.1111/j.1558-5646.2007.00136.x.
Linnen CR, Farrell BD: A test of the sympatric host race formation hypothesis in Neodiprion(Hymenoptera: Diprionidae). Proc R Soc Lond B. 2010,
Ehrlich PR, Raven PH: Butterflies and plants: a study in coevolution. Evolution. 1964, 18: 586-608. 10.2307/2406212.
Diehl SR, Bush GL: An evolutionary and applied perspective of insect biotypes. Annu Rev Entomol. 1984, 29: 471-504. 10.1146/annurev.en.29.010184.002351.
Ronquist F, Liljeblad J: Evolution of the gall wasp-host plant association. Evolution. 2001, 55: 2503-2522.
Schluter D, Price T, Mooers AØ, Ludwig D: Likelihood of ancestor states in adaptive radiation. Evolution. 1997, 51: 1699-1711. 10.2307/2410994.
Kopelke J-P: Gall-forming Nematinae, their willow hosts (Salix spp.) and biological strategies (Insecta, Hymenoptera, Symphyta, Tenthredinidae, Nematinae: Euura, Phyllocolpa, Pontania). Senckenb Biol. 2003, 82: 163-189.
Roininen H, Nyman T, Zinovjev AG: Biology, ecology, and evolution of gall-inducing sawflies (Hymenoptera: Tenthredinidae and Xyelidae). Biology, Ecology, and Evolution of Gall-inducing Arthropods. Edited by: Raman A, Schaefer CW, Withers TM. 2005, Enfield, Science Publishers Inc, 467-494.
Price PW: Adaptive radiation of gall-inducing insects. Basic Appl Ecol. 2005, 6: 413-421. 10.1016/j.baae.2005.07.002.
Taeger A, Altenhofer E, Blank SM, Jansen E, Kraus M, Pschorn-Walcher H, Ritzau C: Kommentare zur Biologie, Verbreitung und Gefährdung der Pflanzenwespen Deutschlands (Hymenoptera, Symphyta). Pflanzenwespen Deutschlands (Hymenoptera: Symphyta). Edited by: Taeger A, Blank SM. 1998, Keltern, Goecke & Evers, 49-135.
Farrell BD, Mitter C, Futuyma DJ: Diversification at the insect-plant interface. Insights from phylogenetics. BioScience. 1992, 42: 34-42. 10.2307/1311626.
Becerra JX, Noge K, Venable DL: Macroevolutionary chemical escalation in an ancient plant-herbivore arms race. Proc Natl Acad Sci USA. 2009, 106: 18062-18066. 10.1073/pnas.0904456106.
McKenna DD, Farrell BD: Tropical forests are both evolutionary cradles and museums of leaf beetle diversity. Proc Natl Acad Sci USA. 2006, 103: 10947-10951. 10.1073/pnas.0602712103.
McLeish MJ, Chapman TW, Schwarz MP: Host-driven diversification of gall-inducing Acacia thrips and the aridification of Australia. BMC Biol. 2007, 5: 3-10.1186/1741-7007-5-3.
Winkler IS, Mitter C, Scheffer SJ: Repeated climate-linked host shifts have promoted diversification in a temperate clade of leaf-mining flies. Proc Natl Acad Sci USA. 2009, 106: 18103-18108. 10.1073/pnas.0904852106.
Wilf P, Labandeira CC, Johnson KR, Coley PD, Cutter AD: Insect herbivory, plant defense, and early Cenozoic climate change. Proc Natl Acad Sci USA. 2001, 98: 6221-6226. 10.1073/pnas.111069498.
Ross HH: The origin of species diversity in ecological communities. Taxon. 1972, 21: 253-259. 10.2307/1218192.
Esselstyn JA, Timm RM, Brown RM: Do geological or climatic processes drive speciation in dynamic archipelagos? The tempo and mode of diversification in Southeast Asian shrews. Evolution. 2009, 63: 2595-2610. 10.1111/j.1558-5646.2009.00743.x.
Procheş Ş, Forest F, Veldtman R, Chown SL, Cowling RM, Johnson SD, Richardson DM, Savolainen V: Dissecting the plant-insect diversity relationship in the Cape. Mol Phylogenet Evol. 2009, 51: 94-99. 10.1016/j.ympev.2008.05.040.
Taeger A, Blank SM: ECatSym - Electronic World Catalog of Symphyta (Insecta, Hymenoptera). Program version 3.9, data version 34. 2008, Müncheberg, Digital Entomological Information
Zachos JC, Dickens GR, Zeebe RE: An early Cenozoic perspective on greenhouse warming and carbon-cycle dynamics. Nature. 2008, 451: 279-283. 10.1038/nature06588.
Argus GW: Infrageneric classification of Salix (Salicaceae) in the New World. Syst Bot Monogr. 1997, 52: 1-121. 10.2307/25096638.
Skvortsov AK: Willows of Russia and adjacent countries. Taxonomical and geographical revision. Univ Joensuu Fac Math Nat Sci Rep Ser. 1999, 39: 1-307.
Kennedy CEJ, Southwood TRE: The number of species of insects associated with British trees: a re-analysis. J Anim Ecol. 1984, 53: 455-478. 10.2307/4528.
Mardulyn P, Milinkovitch MC, Pasteels JM: Phylogenetic analyses of DNA and allozyme data suggest that Gonioctena leaf beetles (Coleoptera; Chrysomelidae) experienced convergent evolution in their history of host-plant family shifts. Syst Biol. 1997, 46: 722-747. 10.1093/sysbio/46.4.722.
Lawton JH, Strong DR: Community patterns and competition in folivorous insects. Am Nat. 1981, 118: 317-338. 10.1086/283826.
Denno RF, McClure MS, Ott JR: Interspecific interactions in phytophagous insects: competition reexamined and resurrected. Annu Rev Entomol. 1995, 40: 297-331. 10.1146/annurev.en.40.010195.001501.
Kaplan I, Denno RF: Interspecific interactions in phytophagous insects revisited: a quantitative assessment of competition theory. Ecol Lett. 2007, 10: 977-994. 10.1111/j.1461-0248.2007.01093.x.
Janzen DH: Host plants as islands in evolutionary and contemporary time. Am Nat. 1968, 102: 592-595. 10.1086/282574.
McKenna DD, Sequeira AS, Marvaldi AE, Farrell BD: Temporal lags and overlap in the diversification of weevils and flowering plants. Proc Natl Acad Sci USA. 2009, 106: 7083-7088. 10.1073/pnas.0810618106.
Kelly CK, Southwood TRE: Species richness and resource availability: a phylogenetic analysis of insects associated with trees. Proc Natl Acad Sci USA. 1999, 96: 8013-8016. 10.1073/pnas.96.14.8013.
Lewinsohn TM, Novotny V, Basset Y: Insects on plants: diversity of herbivore assemblages revisited. Annu Rev Ecol Evol Syst. 2005, 36: 597-620. 10.1146/annurev.ecolsys.36.091704.175520.
Fretwell SD, Lucas HL: On territorial behavior and other factors influencing habitat distribution in birds. I. Theoretical development. Acta Biotheor. 1970, 19: 16-36. 10.1007/BF01601953.
Rabosky DL: Ecological limits and diversification rate: alternative paradigms to explain the variation in species richness among clades and regions. Ecol Lett. 2009, 12: 735-743. 10.1111/j.1461-0248.2009.01333.x.
Tack AJM, Ovaskainen O, Harrison PJ, Roslin T: Competition as a structuring force in leaf miner communities. Oikos. 2009, 118: 809-818. 10.1111/j.1600-0706.2008.17397.x.
Janz N, Nylin S, Wahlberg N: Diversity begets diversity: host expansions and the diversification of plant-feeding insects. BMC Evol Biol. 2006, 6: 4-10.1186/1471-2148-6-4.
Emerson BC, Gillespie RG: Phylogenetic analysis of community assembly and structure over space and time. Trends Ecol Evol. 2008, 23: 619-630. 10.1016/j.tree.2008.07.005.
Azuma T, Kajita T, Yokoyama J, Ohashi H: Phylogenetic relationships of Salix (Salicaceae) based on rbcL sequence data. Am J Bot. 2000, 87: 67-75. 10.2307/2656686.
Chen J-H, Sun H, Wen J, Yang Y-P: Molecular phylogeny of Salix L. (Salicaceae) inferred from three chloroplast datasets and its systematic implications. Taxon. 2010, 59: 29-37.
We especially wish to thank colleagues who provided samples of higher-nematine species, many of which would otherwise have been impossible to obtain: Lauri Kapari, Heikki Roininen, Alexey Zinovjev, Ewald Altenhofer, Mikk Heidemaa, Marko Prous, Akira Yamagami, Matti Viitasaari, Herbert R. Jacobson, Valerie Caron, Urs Schaffner, Jens Rydell, Jeffrey Joy, and B. DeJonge. We also thank two anonymous referees for their constructive criticisms, which helped to improve the manuscript. The Centre of Scientific Computing in Helsinki and the Cyberinfrastructure for Phylogenetic Research (CIPRES) project in San Diego allocated computing resources for phylogenetic analyses on their servers. This research was initiated with support from the SYNTHESYS Project http://www.synthesys.info which is financed by the European Community Research Infrastructure Action under the FP6 "Structuring the European Research Area Programme" (project BE-TAF-1462), and the main part of the funding was provided by the Academy of Finland (project 124695 for TN).
The study was conceived and designed by TN and JLB. TN performed laboratory and data analyses, prepared figures, and wrote the manuscript. JLB planned taxon sampling and study setup, and assisted in writing. VV and DRS provided taxonomic and ecological background information and identified specimens used in the analyses. All authors read and approved the manuscript.