The origins of species richness in the Hymenoptera: insights from a family-level supertree
© Davis et al; licensee BioMed Central Ltd. 2010
Received: 24 September 2009
Accepted: 27 April 2010
Published: 27 April 2010
The order Hymenoptera (bees, ants, wasps, sawflies) contains about eight percent of all described species, but no analytical studies have addressed the origins of this richness at family-level or above. To investigate which major subtaxa experienced significant shifts in diversification, we assembled a family-level phylogeny of the Hymenoptera using supertree methods. We used sister-group species-richness comparisons to infer the phylogenetic position of shifts in diversification.
The supertrees most supported by the underlying input trees are produced using matrix representation with compatibility (MRC) (from an all-in and a compartmentalised analysis). Whilst relationships at the tips of the tree tend to be well supported, those along the backbone of the tree (e.g. between Parasitica superfamilies) are generally not. Ten significant shifts in diversification (six positive and four negative) are found common to both MRC supertrees. The Apocrita (wasps, ants, bees) experienced a positive shift at their origin accounting for approximately 4,000 species. Within Apocrita other positive shifts include the Vespoidea (vespoid wasps/ants containing 24,000 spp.), Anthophila + Sphecidae (bees/thread-waisted wasps; 22,000 spp.), Bethylidae + Chrysididae (bethylid/cuckoo wasps; 5,200 spp.), Dryinidae (dryinid wasps; 1,100 spp.), and Proctotrupidae (proctotrupid wasps; 310 spp.). Four relatively species-poor families (Stenotritidae, Anaxyelidae, Blasticotomidae, Xyelidae) have undergone negative shifts. There are some two-way shifts in diversification where sister taxa have undergone shifts in opposite directions.
Our results suggest that numerous phylogenetically distinctive radiations contribute to the richness of large clades. They also suggest that evolutionary events restricting the subsequent richness of large clades are common. Problematic phylogenetic issues in the Hymenoptera are identified, relating especially to superfamily validity (e.g. "Proctotrupoidea", "Mymarommatoidea"), and deeper apocritan relationships. Our results should stimulate new functional studies on the causes of the diversification shifts we have identified. Possible drivers highlighted for specific adaptive radiations include key anatomical innovations, the exploitation of rich host groups, and associations with angiosperms. Low richness may have evolved as a result of geographical isolation, specialised ecological niches, and habitat loss or competition.
One of the greatest challenges in evolutionary biology is to explain heterogeneity in species richness amongst taxa, and in particular why a few notable taxa comprise the majority of species [1–4]. With over half of all described species, the insects pose perhaps the most obvious target group for biologists attempting to tackle this problem . In this paper we address the phylogenetic location of shifts in diversification within one of the largest insect orders, the Hymenoptera (bees, ants, wasps and sawflies), containing some eight percent of all described species.
Phylogenies are useful tools for understanding the evolution of species richness. Since they specify shared common ancestry and absolute or relative taxon age they allow appropriate comparisons to be made amongst taxa, [6–8]. Taxon age in turn is important because for a given positive net rate of cladogenesis, species richness will increase over time. Thus, the species richness of a taxon can only be identified as anomalous if its absolute or relative age is also known. The cladistic and molecular revolutions, which have advanced phylogenetic information, have also therefore stimulated the development of statistical techniques that can best use the available phylogenetic information for macroevolutionary inference [1–5].
One of the most useful pieces of macroevolutionary information that can be extracted from a phylogeny is the identity of clades that are different, relative to others, in their rates of speciation and/or extinction. Once the identity of these exceptional clades is known, hypotheses about underlying causes may be tested , for example relating to adaptive radiations  or key innovations  although this may not always be straightforward . Within the insects, some studies have attempted to do this at level of order or family [13–15], but within orders macroevolutionary studies have generally focussed on a small subset of taxa [16–20], which places obvious constraints on the explanatory potential of the study. A notable exception is the study of Hunt et al.  using a phylogeny of nearly 2,000 beetle species to estimate shifts in diversification across the order. Consistent with a comparable study across the angiosperm families , they detected numerous, both positive and negative, shifts in diversification.
There are four insect orders with over 100,000 described species, of which the Hymenoptera is one [23–25]. Little work has addressed the evolutionary origins of this diversity. Below family-level, ant diversification has been attributed to co-evolution with angiosperms , fig wasp diversification has been attributed to co-speciation with their host plants , nematine sawfly diversification has been linked to the exploitation of a larger host plant range , and the limited diversification of bumble bees has been attributed to biogeographical constraints , but the influence of such factors at higher taxonomic levels within the order have not been investigated. Studies addressing the role of parasitism in insect diversification more generally have also included hymenopteran data. For example, Mitter et al.  show increased diversification in phytophagous groups of insect (including the hymenopterous sawflies), but the same authors  found no evidence to suggest carnivorous parasitism (including the hymenopterous Parasitica) itself, enhanced diversification.
Clearly, if studies of Hymenoptera macroevolution are to be conducted at large taxonomic scales, appropriate phylogenies will be needed. Few large-scale hymenopteran phylogenies have been produced, and those that have are based solely on morphological data [30–34]. The most recent large-scale morphological phylogeny was constructed in 1999 . Since 1999, and Ronquist's influential "State of the Art" hymenopteran phylogenetics symposium  a wealth of Hymenoptera phylogenies have been published. Perhaps most significant is the increased contribution of molecular work, in its infancy in 1999, alongside morphological work. Most molecular work however, has been fairly focussed, concentrating on certain subsets of hymenopteran families from the two traditional suborders Symphyta (sawflies) [36, 37], or Apocrita (bees, ants and wasps) [38, 39], or at even lower taxonomic levels such as superfamily [40, 41].
Here we take available molecular, morphological and palaeontological information to construct a family-level supertree of the Hymenoptera including extinct taxa, presenting, in terms of the scope of its evidence base, the largest phylogenetic analysis of Hymenoptera to date. The supertree is used to examine the points at which Hymenoptera have experienced significant shifts in diversification. We then explore possible reasons for these shifts in the light of previous hypotheses, and their ability to explain Hymenoptera species richness.
For consistency with other analyses by us on other orders  the taxonomic nomenclature is taken from a family list generated from Gordh & Headrick  for living families. There are slight differences between this list and other authoritative lists for Hymenoptera (see additional file 1): the list generally follows Gauld & Bolton , although departs significantly in splitting Apidae (sensu lato) into several families, following Goulet & Huber . Gordh & Headrick  also include taxa that were recognised prior to Gauld & Bolton  but have since been assigned subfamily status by most researchers (e.g. Loboscelidiidae). By including such taxa in our input tree search, we allowed maximal opportunity for data coverage where phylogenetic studies distinguish them. Where possible we also included extinct families (despite their exclusion from diversification analyses below). A benefit of including extinct taxa is that they can help resolve phylogenetic problems , and may aid future attempts to date the tree. We used the nomenclature of Ross & Jarzembowski  and the regularly updated EDNA Fossil Insect Database http://edna.palass-hosting.org/search.php for extinct families. Our diversification analyses were only conducted with extant taxa (see below). For a full list of taxa and alternative taxonomies see additional file 1.
Input trees were sought online via Web of Science and Google Scholar using the following keywords: hymenoptera*, apocrita*, aculeata*, parasitica*, symphyta*, phylogen*, cladistic, cladogram. (* where any word begins with those letters). Further source studies were collected from the reference sections of these papers. Potential input trees from the collected literature were extracted manually and compiled into a dataset of hymenopteran phylogenies. Hennig  is often cited as a milestone in hexapod phylogenetics (synthesizing his initial cladistic views), and only studies post-1969 were considered. New papers on hymenopteran phylogeny are regularly published, so a cut-off date of December 2008 was used to facilitate future updates of our work.
A total of 101 input trees were obtained from the literature, which contain information on hymenopteran family relationships. However, not all of these could be retained for the final analysis as there is data non-independence (i.e. where different studies overlap in the primary data they use) between them. This is a major challenge with supertree analyses when they are used to summarise published phylogenies (but not when supertrees are used as phylogenomic tools ) as it is often impossible to fully remove data non-independence without unnecessary loss of other information. Following, and slightly modifying, the guidelines of Bininda-Emonds et al.  we attempt to minimise data non-independence as far as possible. Choosing which input trees to keep is done as follows:
Where data overlap is under 50% we retain both trees
Where data overlap is over 50% but taxon overlap is under 50% we keep both trees
Where data and taxon overlap is 50% or over, we first take the author's preferred hypothesis, because supertrees are designed to synthesise current hypotheses.
Where no tree is explicitly preferred, we take the most comprehensive tree, because supertrees are designed to be comprehensive.
Where both trees are equally comprehensive we take the most recent tree, under the notion that inferences converge on the truth with time.
Where both trees are equally recent, we use a tree that has support measures provided over one that does not, as a more reliable phylogeny.
Where trees cannot be chosen based on the above criteria we use all available trees and downweight them in the final analysis accordingly.
Of the 101 trees in the dataset, 24 were removed based on the above protocol. Of the families listed in Gauld & Bolton , only four were absent from input trees: the Rotoitidae (Chalcidoidea), the Eucoilidae (Cynipoidea), the Peredeniidae (Proctotrupoidea) and the Austroniidae (Proctotrupoidea). For a list of input trees and data non-independence see additional file 2. Input tree files are available on request from RBD.
Within these trees, valid taxa (see taxonomic nomenclature) were sometimes portrayed either as paraphyletic or polyphyletic, for example in molecular data sets where a single family could be represented by several species. In these instances polyphyletic families were removed entirely as their placement is uncertain. Paraphyletic families were condensed into a single branch. This issue could be dealt with in other ways. For example weighted sets of input trees representing alternative placements of such taxa could be used. Our approach is more pragmatic than this. Alternatively, such taxa could be entirely removed from supertree analysis, but this could result in unnecessary loss of information. The effects of other approaches to dealing with this issue are unknown.
Many previous supertree studies [22, 55–64] have used the popular matrix representation with parsimony (MRP) method [65, 66] only, though some more recent studies use multiple methods [67, 68]. We implement three supertree methods; two matrix-based and one distance-based, to avoid over-reliance on a single method, and to ensure we obtain as optimal a phylogeny as possible. In addition to MRP, we use matrix representation with compatibility (MRC (also known as split-fit)) , and the average consensus method . Although the average consensus approach can use branch length information, here we do not as many trees do not have such information available. We avoided using strict supertree methods as they only identify relationships common to all input trees . For software and settings see additional file 3.
Analyses of the complete data were successfully conducted using the average consensus and MRC methods. However, analysis of the full data set could not be conducted using MRP. There is a high amount of conflict in the dataset, which makes the analysis more computationally challenging. Furthermore poorly represented taxa (i.e. with much missing data in the MRP matrix) cloud the analysis. We note that different supertree methods are implemented in different software and the number of optimal trees returned will reflect not only the data but also the efficiency of tree searching. To enable analysis using MRP, in addition to MRC and the average consensus the dataset was split into three widely accepted subgroups within Hymenoptera: Symphyta, Parasitica and Aculeata. Since Symphyta and Parasitica are likely not monophyletic (Figure 1) the most widely represented member of the relevant nested clade, here Vespidae, was included in these analyses to test for paraphyly and to provide a reference point in order to graft the three parts of the hymenopteran tree together. Including more nested aculeate families has no effect on the outcome of the analyses. Separate supertree analyses were run for these three groups of hymenopteran taxa. The following summarises the analyses undertaken using each technique:
Standard MRP: Separate Symphyta, Parasitica and Aculeata supertrees produced and grafted together. This constrains the monophyly of Aculeata and Apocrita. Due to large amounts of conflict in the data set, Parasitica families that were underrepresented (i.e. < 10 pseudocharacters in the matrix) or only in input tree polytomies had to be removed to ensure the analysis could be completed. Fully bifurcating trees were required for our diversification analysis. Extended majority rule trees display all relationships in 50% or more of the most parsimonious trees (MPTs) recovered in supertree analysis. Where less than 50% of these MPTs agree on one relationship, the relationship with the largest proportion of sub-50% agreement between MPTs is displayed producing a fully bifurcating tree. In the extended majority rule setting symphytan and aculeate taxa in polytomies were additionally removed. Their inclusion cannot be justified as the extended majority rule method will arbitrarily place them in bifurcating splits without any precise information on this.
MRC: Separate Symphyta, Parasitica and Aculeata supertrees were produced. Problematic Parasitica did not need removing, but could be removed for direct comparison with standard MRP results. Analyses were also run using the all-inclusive dataset and the supertree compared to the re-grafted one to test the effects of constraining monophyly. In the extended majority rule setting symphytan and aculeate taxa in polytomies were additionally removed.
Average consensus: Separate Symphyta, Parasitica and Aculeata supertrees were produced. Underrepresented Parasitica did not need removing, but all taxa only in polytomies were removed as the average consensus method produces only one fully bifurcating tree. The reasons behind the removal of these taxa are the same as if an extended majority rule tree was to be constructed (see above). Analyses were also run using the all-inclusive dataset.
We use the V index of Wilkinson et al.  to measure support. This considers the number of input trees in agreement and in conflict with relationships in the supertree on a scale running from -1 (all conflicting) to +1 (all supporting). For a given relationship in a supertree, each input tree can either support or conflict with this, where the taxa involved in the relationship are represented. Each input tree can therefore contribute a score of +1 (indicating support) or -1 (indicating conflict) to the score for a relationship in a supertree. All +1s and -1s are added and a mean calculated to provide scores for individual relationships. The mean of these nodal scores can also be calculated to provide an overall V score for a supertree allowing comparison between supertrees produced using different methods. V+ is a more liberal score, which considers permitting relationships in input trees (e.g. polytomies) as support (+1), and is also used here. V scores were calculated using the software stsupport (obtained in January 2007 from http://webspace.qmul.ac.uk/jacotton/software/stsupport.html#Download). For clarity of presentation we use branch thickness to visualise support for relationships in figures and provide a table of V and V+ scores in additional file 4. An alternative supertree-specific support measure, rQS [73, 74] was not used as this has been found to be unreliable [75, 76]. Output files are available on request from RBD.
We compared sister taxon species richness to assess where significant shifts in diversification have occurred during the evolution of the Hymenoptera, as used previously in studies of insect species richness [13, 14]. This analysis does not require information on dates of branching events, which we have not attempted to infer. If sister taxa radiate at equal, but not necessarily constant, rates over time, including extinction  all possible sizes of splits of N species between these two clades are equally likely . The two-tailed probability of an equal or greater magnitude of split under the null hypothesis of equal net diversification rates is given by 2 [N small /(N small + N large - 1)].
Species richness for most families was taken from Goulet & Huber . Most of the values are estimates of described species richness. However, the number of described species is not accurately known for the two largest families, Braconidae and Ichneumonidae. Estimates vary markedly from source to source, and both are known to be markedly underdescribed. We therefore also follow Goulet & Huber  in using the widely cited estimates of 60,000 spp. for Ichneumonidae  and 40,000 spp. for Braconidae , which although far larger than the number of described species are still possibly underestimates of true species richness. Alternative lower estimates have no effect on the results below. If estimates were not given in Goulet & Huber , they were taken from Parker , except for Stenotritidae  and Pterygophoridae . Lastly, an outgroup to all Hymenoptera was required in order to test for a significant difference in species richness at the basal node of the tree. Whether we take the panorpidan orders (Diptera, Lepidoptera, Mecoptera, Siphonaptera, Strepsiptera and Trichoptera) [76, 85, 86] or all other Holometabola [87–89] as the sister group to Hymenoptera, has no effect on our overall findings. We take a combined described species estimate from Grimaldi & Engel  for this purpose.
Shifts observed in supertrees from all-in and compartmentalised MRC analyses.
Group in which shift occurs. Species numbers in brackets (overall/those involved in shift (all-in analysis))
Direction of shift
Sister group. Species number in brackets (overall/those involved in shift (all-in analysis except *))
Other Tenthredinoidea (with or without Xyelotomidae)
Bethylidae + Chrysididae
Sclerogibbidae + Dryinidae + Embolemidae
Vespoidea (excl. Sierolomorphidae)
Apoidea* (21,842/177) or Sierolomorphidae (10/10)
Anthophila + Sphecidae
Two-Way Shifts (Both Trees)
Agaonidae + Mymaridae + Signiphoridae + Mymarommatidae + Aphelinidae + Trichogrammatidae
Aphelinidae + Trichogrammatidae
Anthophoridae + Apidae + Megachilidae
The Vespoidea (excluding Sierolomorphidae) and the Anthophila + Sphecidae (forming the majority of Apoidea) are also recovered as significant positive shifts in diversification with possibly sizeable effects on species richness (Table 1). Other shifts are confined to smaller clades comprising only one or two families: upshifts in diversification in Bethylidae + Chrysididae, in Dryinidae, and in Proctotrupidae, and downshifts in diversification in Stenotritidae, Anaxyelidae, Blasticotomidae and Xyelidae. Only one of these shifts is linked to a phylogenetic topology not recovered in the standard MRP tree: in the MRP tree Sierolomorphidae nests within Vespoidea as sister to Formicidae + Falsiformicidae (i.e. ants), and not at the base of this superfamily.
In addition to these positive and negative shifts in diversification are six examples of a significant difference in species richness between sister taxa, but not a significant difference with respect to their outgroup. Two of these involve Chalcidoidea families, two involve Apoidea families, and one each involves cynipoid and vespoid families (Figure 4). Each tree also produces a small number of shifts in diversification not found in the other tree (see additional file 6); these are generally away from the backbone of the tree with the exception of a possible positive shift of all Hymenoptera to the exclusion of Xyelidae in the compartmentalised MRC tree. Many of these are poorly supported groups or the relationship with their sister group is poorly supported (see V scores in additional file 4).
We have produced a family-level phylogeny of Hymenoptera that is the most comprehensive to date in terms of its evidence base. We include almost all living and many fossil families. As a result we have also identified several likely shifts in diversification. Below we first discuss the implications of our trees for the future of Hymenoptera phylogenetics. We then discuss the implications of the identified shifts in diversification.
This supertree approach has yielded findings both supporting and questioning previous hypotheses across the hymenopteran tree. Several relationships receive full support (V = 1) and appear in both MRC trees, indicating that the phylogenetic evidence, as used in our analysis, is in full agreement: Argidae + Pergidae, Cephoidea, Orussoidea (Symphyta); Aphelinidae + Trichogrammatidae, Ceraphronoidea, Megaspilidae + Ceraphronidae, Evaniidae + Cretevaniidae (Parasitica); Anthophila, Megachilidae + Apidae + Anthophoridae (Aculeata). In addition to these, two sets of relationships receive a V score of over 0.8: Argidae + Pergidae + Pterygophoridae (Symphyta), and Eupelmidae + Encyrtidae (Parasitica). Furthermore, symphytan superfamily relationships are consistently recovered with a mostly well-supported branching order in supertree analyses of (Xyelidae, (Tenthredinoidea, (Megalodontoidea, (Cephoidea, (Siricoidea, (Xiphydriidae, (Orussoidea + Apocrita))))))), with Xiphydriidae outside of any existing superfamily. The latter has often remained unclassified at a superfamily level . Here we find no evidence to place it in an existing superfamily either, and its position, is in agreement with more recent symphytan studies [37, 42] Full V scores are provided in additional file 4. Hence, our analyses suggest that the above sets of relationships can be regarded as well known and supported.
The most concerning lack of resolution regards apocritan superfamily relationships, and is highlighted by two MRC and one MRP trees with alternative topologies (Figures 2, 4, additional files 7, 8). However, the shifts in diversification we identified earlier are stable to these differences (see below). No group is strongly supported as the sister group to Ichneumonoidea. Tentatively, the Stephanidae might be closely related as in both MRC and MRP analyses they are either singularly or in a more inclusive clade as sister to Aculeata yet this is poorly supported (V = -0.728 or -0.642 respectively). Furthermore, the validity of some traditionally recognised superfamilies is uncertain. Not all extended majority rule analyses show Vespoidea to be monophyletic. Vespoid monophyly has previously been questioned , though based on features of Pompilidae (here shown to group in Vespoidea) and not Sierolomorphidae, which may sit outside Vespoidea.
Arguably our most interesting finding lies in the phylogeny of Chalcidoidea and Proctotrupoidea. Chalcidoidea are usually well accepted as a superfamily [24, 90] but here Mymarommatidae (the single family in Mymarommatoidea) always nest within them with good support in agreement with Krogmann & Vilhelmsen . The status of Mymarommatoidea is dubious though. It has previously been included within Chalcidoidea when Mymarommatidae nests within them, and as a separate superfamily otherwise. Its taxonomy seems phylogeny-dependent, and results here suggest Mymarommatidae be included in the Chalcidoidea . Despite this, the Chalcidoidea are not always recovered as monophyletic (position of Tanaostigmatidae - Figure 4), but evidence in input trees for this is lacking (V = 0).
The validity of Proctotrupoidea has previously been questioned, for example in some input studies [39, 41] and through agreement in supertree analyses there may be a case for revising their higher taxonomy. It is of interest to see which families consistently group together as a basis for revising the taxonomy of this group, though we do not here suggest any specific revisions. With strong support (V = 0.6) Proctotrupidae and Vanhorniidae are recovered as sister groups in agreement with Castro & Dowton  and Dowton & Austin  with Heloridae basal to this in all analyses. Diapriidae and Monomachidae are recovered in all analyses with reasonable support (V = 0.179), while Roproniidae and Pelecinidae are recovered in all analyses but with poor support (V = -0.234). This relationship is also in agreement with Castro & Dowton  and Dowton & Austin . These problematic relationships and arguably a revision of "proctotrupoid" taxonomy should be priorities for future research.
Diversification Shift Analyses
Because we have conducted the first quantitative tests of shifts in diversification at higher taxonomic levels in Hymenoptera, most of the inferred shifts are entirely novel hypotheses. Two important resulting questions are: how robust are these inferences, and what might have caused them? Whilst our analyses can directly inform us about the former, the phylogenetic positions of shifts constrain their possible causes but are not themselves hard evidence about causes: any difference between the two clades concerned might be responsible. A particular problem is that some ecological differences between clades noted below, such as the taxonomic extent of hosts and geographical distribution, might have caused shifts in diversification or might equally well be the result of them. Nonetheless, such differences help us to raise hypotheses for further tests.
The deepest consistent shift in diversification we infer, although not having the greatest effect on species richness (Table 1; for species richness calculations see additional file 9), occurs at the origin of Apocrita. The sister pairing of Orussidae (plus Paroryssidae) and Apocrita receives poor support like most of the nodes along the supertree backbone, but this relationship does appear in several studies where it is often strongly supported [33, 37, 42, 43, 92]. Where alternative sister-pairings exist based on both morphological  and molecular data [36, 93] the Apocrita are still shown to be sister to another single symphytan family and our findings are not affected.
Orussidae are the only carnivorous parasitic symphytans . Apocritan taxa are united by the move to a more refined carnivorous parasitic lifestyle, which probably originated out of an orussid-style life cycle. Therefore it is possible that carnivorous parasitism is a prerequisite for other key innovations, which are discussed below (i.e. a case of contingent radiation ). Apocritan larvae have a closed hindgut [44, 94, 95]. This might function to increase assimilation of nutrients or avoid fouling the host , and avoid detection by ovipositing orussids, which use faeces to detect hosts . In the adult, the constriction between abdominal segments 1 and 2, (i.e. the "wasp waist") is a defining feature of apocritan taxa [44, 49, 94]. This might function to facilitate oviposition or improve flight . Furthermore, Sharkey  suggests that the variety of ovipositor morphologies has had a key role in facilitating the diversity of apocritan species and may represent an adaptive radiation towards using a range of host species.
We observe a positive shift in diversification the Vespoidea (Table 1) to the exclusion of Sierolomorphidae in both MRC trees despite different proposed sister groups (either Sierolomorphidae (10 species) or Apoidea). The shift detected could have led to the origin of over 40,000 species had negative shifts in diversification further up the tree not occurred (where these are detected depends on the phylogeny used). Though a sister relationship with Apoidea gains greater support, the size of the potential shift in diversification is greater when Sierolomorphidae are the sister group (p = 0.0004 in contrast to 0.0146). The other Vespoidea differ from Sierolomorphidae in a few anatomical features, though they lack the constriction observed in other Vespoidea (except Pompilidae) between the second and third abdominal segments [48, 49]. This extra vespoid constriction may act to enhance any dexterity conferred by the apocritan "wasp waist" in an example of contingent radiation , and could enable effective use of the sting in defence . The specialised articulation between the second and third abdominal segments is also absent from Apoidea.
The Anthophila (bees) and Sphecidae, together containing over 21,500 species, form a well-supported clade (V = 0.467) and together represent a positive shift in diversification. The taxonomy used here  includes sphecid subfamilies which have family status elsewhere (e.g. Crabronidae, Pemphredonidae, Astatidae). Even when the species numbers, which these families contribute to our Sphecidae, are removed (a ten-fold reduction), this positive shift is still observed. In a similar, but less drastic, manner to that observed in Vespoidea, the potential number of species we might have expected from this shift could be greater (another potential 2,000 species) had negative shifts in diversification further up the tree not occurred. Their sister group Heterogynaidae (well supported, V = 0.467) constitutes just five species and is poorly known. While heterogynaid females are brachypterous (i.e. have reduced, non-functional wings) [49, 98] those of the Sphecidae and the bees are not [48, 49]. Reduced wings may result in reduced dispersal. Heterogynaidae species are found only from Madagascar to the Mediterranean , whereas bees and sphecids have a wider distribution . Furthermore, adult sphecids feed at flowers as do all bees, but whilst larvae of bees are also phytophagous, the larvae of sphecids retain a carnivorous habit . A move towards an association with angiosperms has been cited as a reason for insect diversification in similar studies on other groups [17, 19, 20].
The Bethylidae and Chrysididae are a fairly well-supported sister pairing (V = 0.111) and combined have experienced a significant positive shift in diversification compared to their sister group of Sclerogibbidae, Dryinidae and Embolemidae (Table 1). The monophyly of these five families is not well-supported (V = -0.091, but V+ = 0.333), but these taxa are sometimes in polytomies [32, 34], and Bethylidae and Chrysididae appear in trees that do not include the other families and are consequently shown to be sister to alternative families [100, 101]. Bethylids are parasitoids of Lepidoptera and Coleoptera, while chrysidids are parasitoids primarily of other Hymenoptera and Lepidoptera. These host groups are themselves extremely speciose, compared to the host groups of Sclerogibbidae (Embioptera hosts), Dryinidae and Embolemidae (both Homoptera hosts) [48, 49].
A significant positive shift is observed in the Dryinidae compared to their sister family Embolemidae (Table 1), a sister pairing that receives strong support (V = 0.667). Most dryinids possess an enlarged claw on the female protarsus, which could confer increase raptorial ability, consequently enhancing capacity for parasitisation [49, 102]. The host groups of Dryinidae are two large superfamilies of Hemiptera (Cicadelloidea and Fulgoroidea) , making their potential ecological niche broad.
The Proctotrupidae (310 species) are a positive radiation relative to their sister group Vanhorniidae (5 species), a sister pairing with good support. The Vanhorniidae are parasites on just one particular coleopteran family - Eucnemidae (1100 species ), while the Proctotrupidae parasitise members of the Diptera and Coleoptera  including two of the largest coleopteran families Carabidae (20,000 species) and Staphylinidae (32,000 species - ).
The only significant negative shift occurring within the Apocrita concerns the family Stenotritidae (Table 1). This family is sister to Colletidae (V = 0.5). While Colletidae (comprising 2,000 species) has a wide distribution through the southern hemisphere (with genera Colletes and Hylaeus wider still ), Stenotritidae is found solely in Australia .
Three symphytan families have likely undergone significant negative shifts in diversification. The Anaxyelidae comprises just one modern species (Syntexis libocedrii) surviving from a previously larger Mesozoic/Cenozoic group . It is sister to Siricidae (V = -0.270 but V+ = 0.316). Syntexis libocedrii has a very specialised lifestyle, where the female oviposits only into recently burnt or still burning cedar wood specifically . The Blasticotomidae, comprising 9 species, also occupies a specialist niche. Blasticotomidae larvae exhibit a phloem-feeding strategy known as "phloemining", which is unique amongst modern families. This habit may have been more widespread previously within tenthredinoid families . They are strongly supported as sister to the other Tenthredinoidea (either including or excluding the extinct family Xyelotomidae), which comprises over 7,500 species (V = 0.556 or 0.569 respectively).
The Xyelidae may also represent a significant negative shift in diversification away from the other Hymenoptera (over 190,000 species). How large this shift is inferred to be depends on the phylogeny studied, and therefore the shifts in diversification detected and adjusted for further up the tree. They are a relict family containing now just 46 species dating back to the Triassic. Carpenter  reports 36 fossil genera compared to the current five. Other relict groups are often inferred to have lost original niche space over time  or to have been hindered in further diversification by the radiation of more derived taxa . In the compartmentalised MRC analysis, a positive shift is also detected in the remaining Hymenoptera. If this inference is true, it could be in part due to competitive replacement of Xyelidae. It is notable that there are no positive shifts in diversification detected amongst symphytan taxa, in contrast to the Apocrita.
In addition to the above shifts in diversification, which are all in one direction, there are six two-way shifts identified common to both MRC trees, where both sister taxa have shifted in opposite directions leading to a significant difference in species richness (Table 1). All occur at tipward nodes within Apocrita. Only one relationship of these six is dubious (Ormyridae v Agaonidae + Mymaridae + Signiphoridae + Mymarommatidae + Aphelinidae + Trichogrammatidae (V = -0.667)). Within this group, there is a two-way shift between Mymarommatidae and Trichogrammatidae + Aphelinidae. Mymarommatidae are globally distributed, but little is known of their biology . Trichogrammatidae and Aphelinidae are egg parasites of taxonomically diverse insect groups including Lepidoptera and some Hemiptera. Similarly, Mutillidae use speciose host groups, while their sister group Sapygidae have a specialist cleptoparasitic lifestyle. Another positive shift is that of the gall-forming Cynipidae, compared to the negative shift in their sister group Ibaliidae, which are parasites of the species-poor Siricoidea . The other two two-way shifts in diversification both occur in Apoidea. A negative shift in Oxaeidae, a geographically restricted group  occurs opposite a positive shift in Adrenidae, which have a wider distribution. Lastly, a negative shift in Melittidae occurs against a positive shift in Anthophoridae + Apidae + Megachilidae. Melittidae are oligolectic, feeding on a narrow range of host plants , and the other families comprise the long-tongued bees, which may represent an adaptive radiation to feed on flowers with more complex structures .
Some further shifts in diversification were inferred in the all-in analysis (Figure 4) but not the compartmentalized one, hence are more tentative. The Ichneumonoidea make up 100,000 species in a possible positive shift. This superfamily displays a number of innovations relating to parasitism (e.g. ovipositor steering mechanisms, polydnaviruses etc ). In the compartmentalised tree the Ichneumonoidea fall within a larger clade representing a positive shift. Other positive shifts in diversification are in large, but poorly supported groups with no obvious connection. Negative shifts in diversification may have occurred in Pterygophoridae and Monomachidae, which have restricted geography [49, 114], and also in the poorly known Rhopalosomatidae and Bradynobaenidae.
This study provides the most comprehensive phylogenetic framework of Hymenoptera to date with which to test evolutionary hypotheses. Such large-scale phylogenies must be regarded as work in progress since there are inevitably many topological uncertainties. With the current evidence we identify many shifts in diversification robustly, which have shaped the evolution of the group, from suborder, through superfamily, to family-level clades. The location of other shifts in diversification is hard to place given current phylogenetic uncertainty. Our findings support other large phylogenetic studies suggesting that the origin of species richness in large clades (such as beetles and angiosperms, and insects as a whole) are due to multiple smaller adaptive radiations with possibly diverse independent causes [13–15, 21]. Thus, there has been much heterogeneity in net rate of change in species richness within Hymenoptera, without which their richness would certainly have been considerably lower, but also might in places have been higher. Notable differences between the clades involved often include morphological innovations in one of the groups, but also differences in host groups and geographic distribution, and these suggest causal hypotheses for future microevolutionary studies on speciation and extinction processes. Our tree provides a benchmark against which future large scale hymenopteran phylogenies can be evaluated.
RBD was supported by BBSRC PhD studentship BB/D527026/1.
- Brooks DR, McLennan DA: The Nature of Diversity. 2002, Chicago: University of Chicago PressGoogle Scholar
- Schluter D: The Ecology of Adaptive Radiation. 2000, Oxford: Oxford University PressGoogle Scholar
- Simpson GG: The Major Features of Evolution. 1953, New York: Columbia University PressGoogle Scholar
- Stanley SM: Macroevolution: pattern and process. 1979, London: John Hopkins Paperbacks;, 2Google Scholar
- Mayhew PJ: Why are there so many insects? Perspectives from fossils and phylogenies. Biological Reviews. 2007, 82: 425-454. 10.1111/j.1469-185X.2007.00018.x.PubMedGoogle Scholar
- Mooers AØ, Heard SJ: Evolutionary process from phylogenetic tree shape. Quarterly Review of Biology. 1997, 72: 31-54. 10.1086/419657.Google Scholar
- Nee S: Inferring speciation rates from phylogenies. Evolution. 2001, 55: 661-668. 10.1554/0014-3820(2001)055[0661:ISRFP]2.0.CO;2.PubMedGoogle Scholar
- Purvis A: Using interspecies phylogenies to test macroevolutionary hypotheses. New Uses for New Phylogenies. Edited by: Harvey PH, Brown AJL, Smith JM, Nee S. 1996, Oxford: Oxford University Press, 153-168.Google Scholar
- Nee S, Barraclough TG, Harvey PH: Temporal changes in biodiversity: detecting patterns and identifying causes. Biodiversity: a biology of numbers and difference. Edited by: Gaston KJ. 1996, Oxford: Oxford University Press, 230-252.Google Scholar
- Gavrilets S, Losos JB: Adaptive radiation: contrasting theory with data. Science. 2009, 323: 732-737. 10.1126/science.1157966.PubMedGoogle Scholar
- Hunter JP: Key innovations and the ecology of macroevolution. Trends in Ecology and Evolution. 1998, 13: 31-36. 10.1016/S0169-5347(97)01273-1.PubMedGoogle Scholar
- de Queiroz A: Contingent predictability in evolution: key traits and diversification. Systematic Biology. 2002, 51: 917-929. 10.1080/10635150290102627.PubMedGoogle Scholar
- Davis RB, Baldauf SL, Mayhew PJ: Eusociality and the success of the termites: insights from a supertree of dictyopteran families. Journal of Evolutionary Biology. 2009, 22: 1750-1761. 10.1111/j.1420-9101.2009.01789.x.PubMedGoogle Scholar
- Mayhew PJ: Shifts in hexapod diversification and what Haldane could have said. Proceedings of the Royal Society of London B. 2002, 269: 969-974. 10.1098/rspb.2002.1957.Google Scholar
- Mayhew PJ: A tale of two analyses: estimating the consequences of shifts in hexapod diversification. Biological Journal of the Linnean Society. 2003, 80: 23-36. 10.1046/j.1095-8312.2003.00217.x.Google Scholar
- Barraclough TG, Vogler AP: Recent diversification rates in North American tiger beetles estimated from a dated mtDNA phylogenetic tree. Molecular Biology and Evolution. 2002, 19: 1706-1716.PubMedGoogle Scholar
- Farrell BD: 'Inordinate fondness' explained: why are there so many beetles?. Science. 1998, 25: 196-198.Google Scholar
- Hines HM: Historical biogeography, divergence times, and diversification patterns of bumble bees (Hymenoptera: Apidae: Bombus). Systematic Biology. 2008, 57: 58-75. 10.1080/10635150801898912.PubMedGoogle Scholar
- Moreau CS, Bell CD, Vila R, Archibald SB, Pierce NP: Phylogeny of the ants: diversification in the age of angiosperms. Science. 2006, 312: 101-104. 10.1126/science.1124891.PubMedGoogle Scholar
- Wheat CW, Vogel H, Wittstock U, Braby MF, Underwood D, Mitchell-Olds T: The genetic basis of a plant-insect coevolutionary key innovation. Proceedings of the National Academy of Sciences of the United States of America. 2007, 104: 20427-20431. 10.1073/pnas.0706229104.PubMed CentralPubMedGoogle Scholar
- Hunt T, Bergsten J, Levkanicova Z, Papadopoulou A, St. John O, Wild R, Hammond PM, Ahrens D, Balke M, Caterino MS, et al: A comprehensive phylogeny of beetles reveals the evolutionary origins of a superradiation. Science. 2007, 318: 1913-1916. 10.1126/science.1146954.PubMedGoogle Scholar
- Davies TJ, Barraclough TG, Chase MW, Soltis PS, Savolainen V: Darwin's abominable mystery: insights from a supertree of angiosperms. Proceedings of the National Academy of Sciences of the United States of America. 2004, 101: 1904-1909. 10.1073/pnas.0308127100.PubMed CentralPubMedGoogle Scholar
- Gaston KJ: The magnitude of global insect species richness. Conservation Biology. 1991, 5: 283-296. 10.1111/j.1523-1739.1991.tb00140.x.Google Scholar
- Grimaldi DA, Engel MS: Evolution of the Insects. 2005, Cambridge: Cambridge University PressGoogle Scholar
- Gullan PJ, Cranston PS: The Insects: An Outline of Entomology. 2005, Oxford: Blackwell PublishingGoogle Scholar
- Jackson AP: Cophylogeny of the Ficus microcosm. Biological Reviews. 2004, 79: 751-768. 10.1017/S1464793104006463.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Mitter C, Farrell B, Wiegmann B: The phylogentic study of adaptive zones: has phytophagy promoted insect diversification?. The American Naturalist. 1988, 132: 107-128. 10.1086/284840.Google Scholar
- Wiegmann BM, Mitter C, Farrell B: Diversification of carnivorous parasitic insects: extraordinary radiation or specialized dead end?. The American Naturalist. 1993, 142: 737-754. 10.1086/285570.Google Scholar
- Königsmann E: Das phylogenetische System der Hymenoptera Teil 2: Symphyta. Deutsche Entomologische Zeitschrift Neue Folge. 1977, 24: 1-40.Google Scholar
- Königsmann E: Das phylogenetische Suystem der Hymenoptera Teil3: Terebrantes (Unterordnung Apocrita). Deutsche Entomologische Zeitschrift Neue Folge. 1978, 25: 1-55.Google Scholar
- Königsmann E: Das phylogenetische System der Hymenoptera Teil 4: Aculeata (Unterordnung Apocrita). Deutsche Entomologische Zeitschrift Neue Folge. 1978, 25: 365-435.Google Scholar
- Rasnitsyn AP: An outline of the evolution of the hymenopterous insects (Order Vespida). Oriental Insects. 1988, 22: 115-145.Google Scholar
- Ronquist F, Rasnitsyn AP, Roy A, Eriksson K, Lindgren M: Phylogeny of the Hymenoptera: a cladisitic reanalysis of Rasnitsyn's (1988) data. Zoologica Scripta. 1999, 28: 13-50. 10.1046/j.1463-6409.1999.00023.x.Google Scholar
- Ronquist F: Phylogeny of the Hymenoptera (Insecta): the state of the art. Zoologica Scripta. 1999, 28: 3-11. 10.1046/j.1463-6409.1999.00019.x.Google Scholar
- Dowton M, Austin AD: Evidence for AT-transversion bias in wasp (Hymenoptera: Symphyta) mitochondrial genes and its implications for the origin of parasitism. Journal of Molecular Evolution. 1997, 44: 398-405. 10.1007/PL00006159.PubMedGoogle Scholar
- Schulmeister S, Wheeler WC, Carpenter JM: Simultaneous analysis of the basal lineages of Hymenoptera (Insecta) using sensitivity analysis. Cladistics. 2002, 18: 455-484.Google Scholar
- Dowton M, Austin AD: Molecular phylogeny of the insect order Hymenoptera: apocritan relationships. Proceedings of the National Academy of Sciences of the United States of America. 1994, 91: 9911-9915. 10.1073/pnas.91.21.9911.PubMed CentralPubMedGoogle Scholar
- Dowton M, Austin AD: Simultaneous analysis of 16S, 28S, COI and morphology in the Hymenoptera: Apocrita - evolutionary transitions among parasitic wasps. Biological Journal of the Linnean Society. 2001, 74: 87-111. 10.1006/bijl.2001.0577.Google Scholar
- Campbell B, Heraty J, Rasplus J-Y, Chan K, Steffen-Campbell J, Babcock C: Molecular systematics of the Chalcidoidea using 28S-D2 rDNA. Hymenoptera: Evolution, Biodiversity and Biological Control. Edited by: Austin AD, Dowton M. 2000, Collingwood, Australia: CSIRO Publishing, 59-71.Google Scholar
- Dowton M, Austin AD, Dillon N, Bartowsky E: Molecular phylogeny of the apocritan wasps: the Proctotrupomorpha and Evaniomorpha. Systematic Entomology. 1997, 22: 245-255. 10.1046/j.1365-3113.1997.d01-42.x.Google Scholar
- Schulmeister S: Review of morphological evidence on the phylogeny of basal Hymenoptera (Insecta) with a discussion of the ordering of characters. Biological Journal of the Linnean Society. 2003, 79: 209-243. 10.1046/j.1095-8312.2003.00232.x.Google Scholar
- Vilhelmsen L: Phylogeny and classification of the extant basal lineages of Hymenoptera. Zoological Journal of the Linnean Society. 2001, 131: 393-442. 10.1111/j.1096-3642.2001.tb01320.x.Google Scholar
- Sharkey MJ: Phylogeny and classification of Hymenoptera. Zootaxa. 2007, 1668: 521-548.Google Scholar
- Gibson GAP: Sister-group relationships of the Platygastroidea and Chalcidoidea (Hymenoptera) - an alternate hypothesis to Rasnitsyn (1988). Zoologica Scripta. 1999, 28: 125-138. 10.1046/j.1463-6409.1999.00015.x.Google Scholar
- Krogmann L, Vilhelmsen L: Phylogenetic implications of the mesosomal skeleton in Chalcidoidea (Hymenoptera, Apocrita) - tree searches in a jungle of homoplasy. Invertebrate Systematics. 2006, 20: 615-674. 10.1071/IS06012.Google Scholar
- Gordh G, Headrick DH: A Dictionary of Entomology. 2000, Wallingford: CABI PublishingGoogle Scholar
- Gauld I, Bolton B: The Hymenoptera. 1988, Oxford: Oxford University PressGoogle Scholar
- Goulet H, Huber JT: Hymenoptera of the World. 1993, Ottawa: Research Branch, Agriculture CanadaGoogle Scholar
- Béthoux O: Gaps and nodes between fossil and extant insects. Systematic Entomology. 2009, 34: 599-604. 10.1111/j.1365-3113.2009.00484.x.Google Scholar
- Ross AJ, Jarzembowski EA: Arthropoda (Hexapoda: Insecta). The Fossil Record 2. Edited by: Benton MJ. 1993, London: Chapman & Hall, 363-426.Google Scholar
- Hennig W: Die Stammesgeschichte der Insekten. 1969, Frankfurt am Main: KramerGoogle Scholar
- Pisani D, Cotton JA, McInerney JO: Supertrees disentangle the chimeric origin of eukaryotic genomes. Molecular Biology and Evolution. 2007, 24: 1752-1760. 10.1093/molbev/msm095.PubMedGoogle Scholar
- Bininda-Emonds ORP, Jones KE, Price SA, Cardillo M, Grenyer R, Purvis A: Garbage in, garbage out: data issues in supertree construction. Phylogenetic Supertrees: Combining Information to Reveal the Tree of Life. Edited by: Bininda-Emonds ORP. 2004, Dordrecht: Kluwer Academic Publishers, 267-280.Google Scholar
- Beck RMD, Bininda-Emonds ORP, Cardillo M, Liu F-GR, Purvis A: A higher-level MRP supertree of placental mammals. BMC Evolutionary Biology. 2006, 6: 93-10.1186/1471-2148-6-93.PubMed CentralPubMedGoogle Scholar
- Bininda-Emonds ORP, Gittleman JL, Purvis A: Building large trees by combining phylogenetic information: a complete taxonomy of the extant Carnivora (Mammalia). Biological Reviews. 1999, 74: 143-175. 10.1017/S0006323199005307.PubMedGoogle Scholar
- Cavalcanti MJ: A phylogenetic supertree of the hammerhead sharks (Carcharhiniformes, Sphyrnidae). Zoological Studies. 2007, 46: 6-11.Google Scholar
- Eo SH, Bininda-Emonds ORP, Carroll JP: A phylogenetic supertree of the fowls (Galloanserae, Aves). Zoologica Scripta. 2009, 38: 465-481. 10.1111/j.1463-6409.2008.00382.x.Google Scholar
- Grotkopp E, Rejmánek M, Sanderson MJ, Rost TL: Evolution of genome sizes in pines (Pinus) and its life-history correlates: supertree analyses. Evolution. 2004, 58: 1705-1729.PubMedGoogle Scholar
- Jones KE, Purvis A, MacLarnon A, Bininda-Emonds ORP, Simmons NB: A phylogenetic supertree of the bats (Mammalia, Chiroptera). Biological Reviews. 2002, 77: 223-259. 10.1017/S1464793101005899.PubMedGoogle Scholar
- Pisani D, Yates AM, Langer MC, Benton MJ: A genus-level supertree of the Dinosauria. Proceedings of the Royal Society of London B. 2002, 269: 915-921. 10.1098/rspb.2001.1942.Google Scholar
- Purvis A: A composite estimate of primate phylogeny. Philosophical Transactions of the Royal Society London B. 1995, 348: 405-421. 10.1098/rstb.1995.0078.Google Scholar
- Ruta M: A species-level supertree for stylophoran echinoderms. Acta Palaeontologica Polonica. 2003, 48: 559-568.Google Scholar
- Salamin N, Hodkinson TR, Savolainen V: Building supertrees: an empirical assessment using the grass family (Poaceae). Systematic Biology. 2002, 51: 136-150. 10.1080/106351502753475916.PubMedGoogle Scholar
- Baum BR: Combining trees as a way of combining datasets for phylogenetic inference, and the desirability of combining gene trees. Taxon. 1992, 41: 3-10. 10.2307/1222480.Google Scholar
- Ragan MA: Phylogenetic inference based on matrix representation of trees. Molecular Phylogenetics and Evolution. 1992, 1: 53-58. 10.1016/1055-7903(92)90035-F.PubMedGoogle Scholar
- Lloyd GT, Davis KE, Pisani D, Tarver JE, Ruta M, Sakamoto M, Hone DWE, Jennings R, Benton MJ: Dinosaurs and the Cretaceous terrestrial revolution. Proceedings of the Royal Society of London B. 2008, 275: 2483-2490. 10.1098/rspb.2008.0715.Google Scholar
- Ruta M, Pisani D, Lloyd GT, Benton MJ: A supertree of Temnospondyli: cladogenetic patterns in the most species-rich group of early tetrapods. Proceedings of the Royal Society of London B. 2007Google Scholar
- Ross HA, Rodrigo AJ: An assessment of matrix representation with compatibility in supertree construction. Phylogenetic Supertrees: Combining Information to Reveal the Tree of Life. Edited by: Bininda-Emonds ORP. 2004, Dordrecht: Kluwer Academic Publishers, 35-63.Google Scholar
- Lapointe F-J, Cucumel G: The average consensus procedure: combination of weighted trees containing identical or overlapping sets of taxa. Systematic Biology. 1997, 46: 306-312.Google Scholar
- Wilkinson M, Thorley JL, Pisani D, Lapointe F-J, McInerney JO: Some desiderata for liberal supertrees. Phylogenetic Supertrees: Combining Information to Reveal the Tree of Life. Edited by: Bininda-Emonds ORP. 2004, Dordrecht: Kluwer Academic Publishers, 227-246.Google Scholar
- Wilkinson M, Pisani D, Cotton JA, Corfe I: Measuring support and finding unsupported relationships in supertrees. Systematic Biology. 2005, 54: 823-831. 10.1080/10635150590950362.PubMedGoogle Scholar
- Bininda-Emonds ORP: Novel versus unsupported clades: assessing the qualitative support for clades in MRP supertrees. Systematic Biology. 2003, 52: 839-848.PubMedGoogle Scholar
- Price SA, Bininda-Emonds ORP, Gittleman JL: A complete phylogeny of the whales, dolphins and even-toed hoofed mammals (Cetartiodactyla). Biological Reviews. 2005, 80: 445-473. 10.1017/S1464793105006743.PubMedGoogle Scholar
- Baker WJ, Savolainen V, Asmussen-Lange CB, Chase MW, Dransfield J, Forest F, Harley MM, Uhl NW, Wilkinson M: Complete generic-level phylogenetic analyses of palms (Arecaceae) with comparisons of supertree and supermatrix approaches. Systematic Biology. 2009, 58: 240-256. 10.1093/sysbio/syp021.PubMedGoogle Scholar
- Davis RB, Baldauf SL, Mayhew PJ: Many hexapod groups originated earlier and withstood extinction better than previously realized: inferences from supertrees. Proceedings of the Royal Society of London B.Google Scholar
- Nee S, May RM, Harvey PH: The reconstructed evolutionary process. Philosophical Transactions of the Royal Society London B. 1994, 344: 305-311. 10.1098/rstb.1994.0068.Google Scholar
- Farris JS: Expected asymmetry of evolutionary rates. Systematic Zoology. 1976, 25: 196-198. 10.2307/2412748.Google Scholar
- Moore B, Chan KMA, Donoghue MJ: Detecting diversification rate variation in supertrees. Phylogenetic Supertrees: Combining Information to Reveal the Tree of Life. Edited by: Bininda-Emonds ORP. 2004, Dordrecht: Kluwer Academic Publishers, 487-533.Google Scholar
- Townes H: The genera of Ichneumonidae 1. Memoirs of the American Entomological Institute. 1969, 11: 1-300.Google Scholar
- van Achterberg C: Essay on the phylogeny of Braconidae. Entomologisk Tidskrift. 1984, 105: 41-58.Google Scholar
- Parker SP: Synopsis and Classification of Living Organisms. 1982, New York: McGraw-HillGoogle Scholar
- Naumann ID: Hymenoptera. The Insects of Australia. Edited by: Naumann ID. 1991, Carlton: Melbourne University Press, 2: 916-1000.Google Scholar
- Schmidt S, Smith DR: An annotated systematic world catalogue of the Pergidae (Hymenoptera). Contributions of the American Entomological Institute. 2006, 34 (3):Google Scholar
- Kristensen NP: Phylogeny of endopterygote insects, the most successful lineage of living organisms. European Journal of Entomology. 1999, 96: 237-253.Google Scholar
- Wheeler WC, Whiting M, Wheeler QD, Carpenter JM: The phylogeny of the extant hexapod orders. Cladistics. 2001, 17: 113-169. 10.1111/j.1096-0031.2001.tb00115.x.Google Scholar
- Ayme-Southgate AJ, Southgate RJ, Philipp RA, Sotka EE, Kramp C: The myofibrillar protein, projectin, is highly conserved across insect evolution except for its PEVK domain. Journal of Molecular Evolution. 2008, 67: 653-669. 10.1007/s00239-008-9177-2.PubMed CentralPubMedGoogle Scholar
- Savard J, Tautz D, Richards S, Weinstock GM, Gibbs RA, Werren JH, Tettelin H, Lercher MJ: Phylogenomic analysis reveals bees and wasps (Hymenoptera) at the base of the radiation of holometabolous insects. Genome Research. 2006, 16: 1334-1338. 10.1101/gr.5204306.PubMed CentralPubMedGoogle Scholar
- Wiegmann BM, Trautwein MD, Kim J-W, Cassel BK, Bertone MA, Winterton SL, Yeates DK: Single-copy nuclear genes resolve the phylogeny of the holometabolous insects. BMC Biology. 2009, 7: 34-10.1186/1741-7007-7-34.PubMed CentralPubMedGoogle Scholar
- Gibson GAP, Heraty JM, Woolley JB: Phylogenetics and classification of Chalcidoidea and Mymarommatoidea - a review of current concepts (Hymenoptera, Apocrita). Zoologica Scripta. 1999, 28: 87-124. 10.1046/j.1463-6409.1999.00016.x.Google Scholar
- Castro LR, Dowton M: Molecular analyses of the Apocrita (Insecta: Hymenoptera) suggest that the Chalcidoidea are sister to the diaprioid complex. Invertebrate Systematics. 2006, 20: 603-614. 10.1071/IS06002.Google Scholar
- Schulmeister S: Simultaneous analysis of basal Hymenoptera (Insecta): introducing robust-choice sensitivity analysis. Biological Journal of the Linnean Society. 2003, 79: 245-275. 10.1046/j.1095-8312.2003.00233.x.Google Scholar
- Dowton M, Austin AD: Increased genetic diversity in mitochondrial genes is correlated with the evolution of parasitism in the Hymenoptera. Journal of Molecular Evolution. 1995, 41: 958-965. 10.1007/BF00173176.PubMedGoogle Scholar
- Hennig W: Insect Phylogeny. 1981, Chichester: John Wiley & SonsGoogle Scholar
- Wharton RA, Vilhelmsen L, Gibson GAP: Characterizing basal apocritans (Hymenoptera: Apocrita). Proceedings of the Russian Entomological Society St Petersburg. 2004, 75: 17-23.Google Scholar
- Powell JA, Turner WJ: Observations on oviposition behavior and host selection in Orussus occidentalis (Hymenoptera: Siricoidea). Journal of the Kansas Entomological Society. 1975, 48: 299-307.Google Scholar
- Ward PS: Ants. Current Biology. 2006, 16: R152-R155. 10.1016/j.cub.2006.02.054.PubMedGoogle Scholar
- Day MC: The enigmatic genus Heterogyna Nagy (Hymenoptera: Sphecidae; Heterogyninae). Systematic Entomology. 1984, 9: 293-307. 10.1111/j.1365-3113.1984.tb00054.x.Google Scholar
- Arévalo HA, Frank JH: Nectar sources for Larra Bicolor (Hymenoptera: Sphecidae), a parasitoid of Scapteriscus mole crickets (Orthoptera: Gryllotalpidae), in Northern Florida. Florida Entomologist. 2005, 88: 146-151. 10.1653/0015-4040(2005)088[0146:NSFLBH]2.0.CO;2.Google Scholar
- Carpenter JM, Wheeler WC: Towards simultaneous analysis of morphological and molecular data in Hymenoptera. Zoologica Scripta. 1999, 28: 251-260. 10.1046/j.1463-6409.1999.00009.x.Google Scholar
- Whiting MF: Phylogeny of the holometabolous insect orders: molecular evidence. Zoologica Scripta. 2002, 31: 3-15. 10.1046/j.0300-3256.2001.00093.x.Google Scholar
- Olmi M: The Dryinidae and Embolemidae (Hymenoptera: Chrysidoidea) of Fennoscandia and Denmark: Chrysidoidea) of Fennoscandia and Denmark. 1994, Leiden: Brill Academic PublishersGoogle Scholar
- Almeida EAB, Danforth BN: Phylogeny of colletid bees (Hymenoptera: Colletidae) inferred from four nuclear genes. Molecular Phylogenetics and Evolution. 2009, 50: 290-309. 10.1016/j.ympev.2008.09.028.PubMedGoogle Scholar
- Zhang H, Rasnitsyn AP: Two new anaxyelid sawflies (Insecta, Hymenoptera, Siricoidea) from the Yixian Formation of western Liaoning, China. Cretaceous Research. 2006, 27: 279-284. 10.1016/j.cretres.2005.11.001.Google Scholar
- Ortega-Blanco J, Rasnitsyn AP, Delclòs X: First record of anaxyelid woodwasps (Hymenoptera: Anaxyelidae) in Lower Cretaceous Spanish amber. Zootaxa. 2008, 1937: 39-50.Google Scholar
- Shcherbakov DE: Fern sawfly larvae Blasticotoma filiceti Klug 1834 (Hymenoptera: Blasticotomidae) are visited by ants: a new kind of trophobiosis. Russian Entomological Journal. 2006, 15: 67-72.Google Scholar
- Carpenter FM: Arthropoda: Superclass Hexapoda. 1992, Boulder: The Geological Society of AmericaGoogle Scholar
- Signor PW: Patterns of diversification. Palaeobiology, A Synthesis. Edited by: Briggs DEG. 1990, Crowther PR: Blackwell, 130-135.Google Scholar
- Connell JH: Diversity and the coevolution of competitors or the ghost of competition past. Oikos. 1980, 35: 131-138. 10.2307/3544421.Google Scholar
- Hurd PD, Linsley EG: The bee family Oxaeidae with a revision of the North American species (Hymenoptera: Apoidea). Smithsonian Contributions to Zoology. 1976, 220: 1-75.Google Scholar
- Michez D, Patiny S, Rasmont P, Timmermann K, Vereecken NJ: Phylogeny and host-plant evolution in Melittidae s.l. (Hymenoptera: Apoidea). Apidologie. 2008, 39: 146-162. 10.1051/apido:2007048.Google Scholar
- Winston ML: The Biology of the Honey Bee. 1991, Cambridge, MA: Harvard University PressGoogle Scholar
- Quicke DLJ: Parasitic Wasps. 1997, London: Chapman & HallGoogle Scholar
- Oelrichs PB, MacLeod JK, Seawright AA, Grace PB: Isolation and identification of the toxic peptides from Lophyrotoma zonalis (Pergidae) sawfly larvae. Toxicon. 2001, 39: 1933-1936. 10.1016/S0041-0101(01)00144-1.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.