Evolutionary History of the Hymenoptera.

Hymenoptera (sawflies, wasps, ants, and bees) are one of four mega-diverse insect orders, comprising more than 153,000 described and possibly up to one million undescribed extant species [1, 2]. As parasitoids, predators, and pollinators, Hymenoptera play a fundamental role in virtually all terrestrial ecosystems and are of substantial economic importance [1, 3]. To understand the diversification and key evolutionary transitions of Hymenoptera, most notably from phytophagy to parasitoidism and predation (and vice versa) and from solitary to eusocial life, we inferred the phylogeny and divergence times of all major lineages of Hymenoptera by analyzing 3,256 protein-coding genes in 173 insect species. Our analyses suggest that extant Hymenoptera started to diversify around 281 million years ago (mya). The primarily ectophytophagous sawflies are found to be monophyletic. The species-rich lineages of parasitoid wasps constitute a monophyletic group as well. The little-known, species-poor Trigonaloidea are identified as the sister group of the stinging wasps (Aculeata). Finally, we located the evolutionary root of bees within the apoid wasp family "Crabronidae." Our results reveal that the extant sawfly diversity is largely the result of a previously unrecognized major radiation of phytophagous Hymenoptera that did not lead to wood-dwelling and parasitoidism. They also confirm that all primarily parasitoid wasps are descendants of a single endophytic parasitoid ancestor that lived around 247 mya. Our findings provide the basis for a natural classification of Hymenoptera and allow for future comparative analyses of Hymenoptera, including their genomes, morphology, venoms, and parasitoid and eusocial life styles.

The Trichoptera barcode initiative: a strategy for generating a species-level Tree of Life.

DNA barcoding was intended as a means to provide species-level identifications through associating DNA sequences from unknown specimens to those from curated reference specimens. Although barcodes were not designed for phylogenetics, they can be beneficial to the completion of the Tree of Life. The barcode database for Trichoptera is relatively comprehensive, with data from every family, approximately two-thirds of the genera, and one-third of the described species. Most Trichoptera, as with most of life's species, have never been subjected to any formal phylogenetic analysis. Here, we present a phylogeny with over 16 000 unique haplotypes as a working hypothesis that can be updated as our estimates improve. We suggest a strategy of implementing constrained tree searches, which allow larger datasets to dictate the backbone phylogeny, while the barcode data fill out the tips of the tree. We also discuss how this phylogeny could be used to focus taxonomic attention on ambiguous species boundaries and hidden biodiversity. We suggest that systematists continue to differentiate between 'Barcode Index Numbers' (BINs) and 'species' that have been formally described. Each has utility, but they are not synonyms. We highlight examples of integrative taxonomy, using both barcodes and morphology for species description.This article is part of the themed issue 'From DNA barcodes to biomes'.

Progress, pitfalls and parallel universes: a history of insect phylogenetics.

The phylogeny of insects has been both extensively studied and vigorously debated for over a century. A relatively accurate deep phylogeny had been produced by 1904. It was not substantially improved in topology until recently when phylogenomics settled many long-standing controversies. Intervening advances came instead through methodological improvement. Early molecular phylogenetic studies (1985-2005), dominated by a few genes, provided datasets that were too small to resolve controversial phylogenetic problems. Adding to the lack of consensus, this period was characterized by a polarization of philosophies, with individuals belonging to either parsimony or maximum-likelihood camps; each largely ignoring the insights of the other. The result was an unfortunate detour in which the few perceived phylogenetic revolutions published by both sides of the philosophical divide were probably erroneous. The size of datasets has been growing exponentially since the mid-1980s accompanied by a wave of confidence that all relationships will soon be known. However, large datasets create new challenges, and a large number of genes does not guarantee reliable results. If history is a guide, then the quality of conclusions will be determined by an improved understanding of both molecular and morphological evolution, and not simply the number of genes analysed.

Mitochondrial capture enriches mito-DNA 100 fold, enabling PCR-free mitogenomics biodiversity analysis.

Biodiversity analyses based on next-generation sequencing (NGS) platforms have developed by leaps and bounds in recent years. A PCR-free strategy, which can alleviate taxonomic bias, was considered as a promising approach to delivering reliable species compositions of targeted environments. The major impediment of such a method is the lack of appropriate mitochondrial DNA enrichment ways. Because mitochondrial genomes (mitogenomes) make up only a small proportion of total DNA, PCR-free methods will inevitably result in a huge excess of data (>99%). Furthermore, the massive volume of sequence data is highly demanding on computing resources. Here, we present a mitogenome enrichment pipeline via a gene capture chip that was designed by virtue of the mitogenome sequences of the 1000 Insect Transcriptome Evolution project (1KITE, www.1kite.org). A mock sample containing 49 species was used to evaluate the efficiency of the mitogenome capture method. We demonstrate that the proportion of mitochondrial DNA can be increased by approximately 100-fold (from the original 0.47% to 42.52%). Variation in phylogenetic distances of target taxa to the probe set could in principle result in bias in abundance. However, the frequencies of input taxa were largely maintained after capture (R(2) = 0.81). We suggest that our mitogenome capture approach coupled with PCR-free shotgun sequencing could provide ecological researchers an efficient NGS method to deliver reliable biodiversity assessment.

The evolutionary history of holometabolous insects inferred from transcriptome-based phylogeny and comprehensive morphological data.

BACKGROUND: Despite considerable progress in systematics, a comprehensive scenario of the evolution of phenotypic characters in the mega-diverse Holometabola based on a solid phylogenetic hypothesis was still missing. We addressed this issue by de novo sequencing transcriptome libraries of representatives of all orders of holometabolan insects (13 species in total) and by using a previously published extensive morphological dataset. We tested competing phylogenetic hypotheses by analyzing various specifically designed sets of amino acid sequence data, using maximum likelihood (ML) based tree inference and Four-cluster Likelihood Mapping (FcLM). By maximum parsimony-based mapping of the morphological data on the phylogenetic relationships we traced evolutionary transformations at the phenotypic level and reconstructed the groundplan of Holometabola and of selected subgroups. RESULTS: In our analysis of the amino acid sequence data of 1,343 single-copy orthologous genes, Hymenoptera are placed as sister group to all remaining holometabolan orders, i.e., to a clade Aparaglossata, comprising two monophyletic subunits Mecopterida (Amphiesmenoptera + Antliophora) and Neuropteroidea (Neuropterida + Coleopterida). The monophyly of Coleopterida (Coleoptera and Strepsiptera) remains ambiguous in the analyses of the transcriptome data, but appears likely based on the morphological data. Highly supported relationships within Neuropterida and Antliophora are Raphidioptera + (Neuroptera + monophyletic Megaloptera), and Diptera + (Siphonaptera + Mecoptera). ML tree inference and FcLM yielded largely congruent results. However, FcLM, which was applied here for the first time to large phylogenomic supermatrices, displayed additional signal in the datasets that was not identified in the ML trees. CONCLUSIONS: Our phylogenetic results imply that an orthognathous larva belongs to the groundplan of Holometabola, with compound eyes and well-developed thoracic legs, externally feeding on plants or fungi. Ancestral larvae of Aparaglossata were prognathous, equipped with single larval eyes (stemmata), and possibly agile and predacious. Ancestral holometabolan adults likely resembled in their morphology the groundplan of adult neopteran insects. Within Aparaglossata, the adult's flight apparatus and ovipositor underwent strong modifications. We show that the combination of well-resolved phylogenies obtained by phylogenomic analyses and well-documented extensive morphological datasets is an appropriate basis for reconstructing complex morphological transformations and for the inference of evolutionary histories.

Advances in insect phylogeny at the dawn of the postgenomic era.

Most species on Earth are insects and thus, understanding their evolutionary relationships is key to understanding the evolution of life. Insect relationships are increasingly well supported, due largely to technological advances in molecular sequencing and phylogenetic computational analysis. In this postgenomic era, insect systematics will be furthered best by integrative methods aimed at hypothesis corroboration from molecular, morphological, and paleontological evidence. This review of the current consensus of insect relationships provides a foundation for comparative study and offers a framework to evaluate incoming genomic evidence. Notable recent phylogenetic successes include the resolution of Holometabola, including the identification of the enigmatic Strepsiptera as a beetle relative and the early divergence of Hymenoptera; the recognition of hexapods as a crustacean lineage within Pancrustacea; and the elucidation of Dictyoptera orders, with termites placed as social cockroaches. Regions of the tree that require further investigation include the earliest winged insects (Palaeoptera) and Polyneoptera (orthopteroid lineages).

Potential pitfalls of modelling ribosomal RNA data in phylogenetic tree reconstruction: evidence from case studies in the Metazoa.

BACKGROUND: Failure to account for covariation patterns in helical regions of ribosomal RNA (rRNA) genes has the potential to misdirect the estimation of the phylogenetic signal of the data. Furthermore, the extremes of length variation among taxa, combined with regional substitution rate variation can mislead the alignment of rRNA sequences and thus distort subsequent tree reconstructions. However, recent developments in phylogenetic methodology now allow a comprehensive integration of secondary structures in alignment and tree reconstruction analyses based on rRNA sequences, which has been shown to correct some of these problems. Here, we explore the potentials of RNA substitution models and the interactions of specific model setups with the inherent pattern of covariation in rRNA stems and substitution rate variation among loop regions. RESULTS: We found an explicit impact of RNA substitution models on tree reconstruction analyses. The application of specific RNA models in tree reconstructions is hampered by interaction between the appropriate modelling of covarying sites in stem regions, and excessive homoplasy in some loop regions. RNA models often failed to recover reasonable trees when single-stranded regions are excessively homoplastic, because these regions contribute a greater proportion of the data when covarying sites are essentially downweighted. In this context, the RNA6A model outperformed all other models, including the more parametrized RNA7 and RNA16 models. CONCLUSIONS: Our results depict a trade-off between increased accuracy in estimation of interdependencies in helical regions with the risk of magnifying positions lacking phylogenetic signal. We can therefore conclude that caution is warranted when applying rRNA covariation models, and suggest that loop regions be independently screened for phylogenetic signal, and eliminated when they are indistinguishable from random noise. In addition to covariation and homoplasy, other factors, like non-stationarity of substitution rates and base compositional heterogeneity, can disrupt the signal of ribosomal RNA data. All these factors dictate sophisticated estimation of evolutionary pattern in rRNA data, just as other molecular data require similarly complicated (but different) corrections.

Molecular dating and biogeography of fig-pollinating wasps.

Figs and fig-pollinating wasps are obligate mutualists that have coevolved for over 60 million years. But when and where did pollinating fig wasps (Agaonidae) originate? Some studies suggest that agaonids arose in the Late Cretaceous and the current distribution of fig-wasp faunas can be explained by the break-up of the Gondwanan landmass. However, recent molecular-dating studies suggest divergence time estimates that are inconsistent with the Gondwanan vicariance hypothesis and imply that long distance oceanic dispersal could have been an important process for explaining the current distribution of both figs and fig wasps. Here, we use a combination of phylogenetic and biogeographical data to infer the age, the major period of diversification, and the geographic origin of pollinating fig wasps. Age estimates ranged widely depending on the molecular-dating method used and even when using the same method but with slightly different constraints, making it difficult to assess with certainty a Gondwanan origin of agaonids. The reconstruction of ancestral areas suggests that the most recent common ancestor of all extant fig-pollinating wasps was most likely Asian, although a southern Gondwana origin cannot be rejected. Our analysis also suggests that dispersal has played a more important role in the development of the fig-wasp biota than previously assumed.

Ancient rapid radiations of insects: challenges for phylogenetic analysis.

Phylogenies of major groups of insects based on both morphological and molecular data have sometimes been contentious, often lacking the data to distinguish between alternative views of relationships. This paucity of data is often due to real biological and historical causes, such as shortness of time spans between divergences for evolution to occur and long time spans after divergences for subsequent evolutionary changes to obscure the earlier ones. Another reason for difficulty in resolving some of the relationships using molecular data is the limited spectrum of genes so far developed for phylogeny estimation. For this latter issue, there is cause for current optimism owing to rapid increases in our knowledge of comparative genomics. At least some historical patterns of divergence may, however, continue to defy our attempts to completely reconstruct them with confidence, at least using current strategies.

Why weight?

Whether phylogenetic data should be differentially or equally weighted is currently debated. Further, if differential weighting is to be explored, there is no consensus among investigators as to which weighting scheme is most appropriate. Mitochondrial genome data offer a powerful tool in assessment of differential weighting schemes because taxa can be selected from which a highly corroborated phylogeny is available (so that accuracy can be assessed), and it can be assumed that different data partitions share the same history (so that gene-sorting issues are not so problematic). Using mitochondrial data from 17 mammalian genomes, we evaluated the most commonly used weighting schemes, such as successive weighting, transversion weighting, codon-based weighting, and amino acid coding, and compared them to more complex weighting schemes including a 6-parameter weighting, pseudoreplicate reweighting, and tri-level weighting. We found that the most commonly used weighting schemes perform the worst with these data. Some of the more complex schemes perform well, however, none of them is consistently superior. These results support ones biases; if one has a predilection to avoid differential weighting, these data support equally weighted parsimony and maximum likelihood. Others might be encouraged by these results to try weighting as a form of data exploration.

Site specific rates of mitochondrial genomes and the phylogeny of eutheria.

BACKGROUND: Traditionally, most studies employing data from whole mitochondrial genomes to diagnose relationships among the major lineages of mammals have attempted to exclude regions that potentially complicate phylogenetic analysis. Components generally excluded are 3rd codon positions of protein-encoding genes, the control region, rRNAs, tRNAs, and the ND6 gene (encoded on the opposite strand). We present an approach that includes all the data, with the exception of the control region. This approach is based on a site-specific rate model that accommodates excessive homoplasy and that utilizes secondary structure as a reference for proper alignment of rRNAs and tRNAs. RESULTS: Mitochondrial genomic data for 78 eutherian mammals, 8 metatherians, and 3 monotremes were analyzed with a Bayesian analysis and our site specific rate model. The resultant phylogeny revealed strong support for most nodes and was highly congruent with more recent phylogenies based on nuclear DNA sequences. In addition, many of the conflicting relationships observed by earlier mitochondrial-based analyses were resolved without need for the exclusion of large subsets of the data. CONCLUSION: Rather than exclusion of data to minimize presumed noise associated with non-protein encoding genes in the mitochondrial genome, our results indicate that selection of an appropriate model that accommodates rate heterogeneity across data partitions and proper treatment of RNA genes can result in a mitochondrial genome-based phylogeny of eutherian mammals that is reasonably congruent with recent phylogenies derived from nuclear genes.

Phylogenetic position of the ant genus Acropyga Roger (Hymenoptera: Formicidae) and the evolution of trophophoresy.

Trophophoresy is exhibited in two ant genera: Acropyga (Formicinae), in which all 37 species are thought to be trophophoretic, and Tetraponera (Pseudomyrmecinae), in which it has been observed in only one species, T. binghami. This study analyses a dataset comprised of both morphological and molecular (D2 region of 28S rRNA and EF1-alpha) data. Evidence is presented in favor of Acropyga being monophyletic, hence trophophoresy has evolved only once within the Formicinae and twice within the ants overall. The data further suggests that Acropyga belongs within a clade containing Anoplolepis, Aphomomyrmex, and Petalomyrmex. Aphomomyrmex and Petalomyrmex were found to be the sister group to Acropyga. The results indicate that the Lasiini and Plagiolepidini are not monophyletic and are in need of reexamination. Given the extant pantropical distribution of Acropyga it is speculated that Acropyga maybe of Gondwanan origin and that trophobiosis was the first form of agriculture to evolve in the ants.

Aligned 18S and insect phylogeny.

The nuclear small subunit rRNA (18S) has played a dominant role in the estimation of relationships among insect orders from molecular data. In previous studies, 18S sequences have been aligned by unadjusted automated approaches (computer alignments that are not manually readjusted), most recently with direct optimization (simultaneous alignment and tree building using a program called "POY"). Parsimony has been the principal optimality criterion. Given the problems associated with the alignment of rRNA, and the recent availability of the doublet model for the analysis of covarying sites using Bayesian MCMC analysis, a different approach is called for in the analysis of these data. In this paper, nucleotide sequence data from the 18S small subunit rRNA gene of insects are aligned manually with reference to secondary structure, and analyzed under Bayesian phylogenetic methods with both GTR+I+G and doublet models in MrBayes. A credible phylogeny of Insecta is recovered that is independent of the morphological data and (unlike many other analyses of 18S in insects) not contradictory to traditional ideas of insect ordinal relationships based on morphology. Hexapoda, including Collembola, are monophyletic. Paraneoptera are the sister taxon to a monophyletic Holometabola but weakly supported. Ephemeroptera are supported as the sister taxon of Neoptera, and this result is interpreted with respect to the evolution of direct sperm transfer and the evolution of flight. Many other relationships are well-supported but several taxa remain problematic, e.g., there is virtually no support for relationships among orthopteroid orders. A website is made available that provides aligned 18S data in formats that include structural symbols and Nexus formats.

Phylogeny and host-plant association in the leaf beetle genus Trirhabda LeConte (Coleoptera: Chrysomelidae).

The leaf beetle genus Trirhabda contains 26 described species from the United States and Canada, feeding on host plants from the families Asteraceae and Hydrophyllaceae. In this study, we present a phylogeny for the genus that was reconstructed from mitochondrial COI and 12S rRNA fragments, nuclear ITS2 rRNA, and morphological characters. Both parsimony and mixed-model Bayesian likelihood analyses were performed. Under both methods, the mitochondrial and nuclear partitions support the same backbone phylogeny, as do the combined data. The utility of the molecular data is contrasted with the low phylogenetic signal among morphological characters. The phylogeny was used to trace the evolution of the host-plant association in Trirhabda. The recovered phylogeny shows that although the host-plant association is phylogenetically conservative, Trirhabda experienced one shift to a distantly related host-plant family, 6 shifts between host-plant tribes, and 6 between genera within tribes. The phylogeny reveals that Trirhabda were plesiomorphically adapted to tolerate complex secondary compounds of its host plants and this adaptation is retained in Trirhabda species, as evidenced by multiple shifts from chemically simpler host plants back to the more complex host plants.

Convergent evolution of cucurbitacin feeding in spatially isolated rootworm taxa (Coleoptera: Chrysomelidae; Galerucinae, Luperini).

Historically, chemical ecologists assumed that cucurbitacin feeding and sequestration in rootworm leaf beetles is a remnant of an ancient association between the Luperini (Coleoptera: Chrysomelidae; Galerucinae) and Cucurbitaceae (ancestral host hypothesis). Under this premise, rootworms that do not develop on cucurbits but undergo pharmacophagous forays for cucurbitacins are thought to do so to supplement novel host diets that lack these bitter compounds. The ancestral host hypothesis is supported from studies of pyrrolizidine alkaloid pharmacophagy in Lepidoptera but has not been subjected to phylogenetic analysis within the Luperini. New evidence that this feeding behavior is better correlated with an adult affinity for pollen than with larval host offers the possibility that Old and New World rootworm species with an affinity for cucurbitacins converged on this behavior through apomorphic taste receptor modifications (loose receptor hypothesis). Here we test the monophyly of cucurbitacin feeding within the Luperini by using nuclear and mitochondrial sequence data to infer phylogenetic relationships among 49 taxa representing tribes of the Galerucinae and subtribes of the Luperini. The resulting phylogenetic hypothesis is mostly concordant with existing tribal and subtribal delineations within the Subfamily Galerucinae sensu stricto (Galerucinae not including the flea beetles). The establishment of ancestry among the subtribes of the Luperini refutes the monophyly of cucurbitacin feeding and cucurbit specialization, with the New World Diabroticina being paraphyletic to the Old World Aulacophorina and cosmopolitan Luperina. These data unambiguously support the convergent evolution of cucurbitacin feeding in rootworms and are inconsistent with the ancestral host hypothesis.

Phylogeny of the photosynthetic euglenophytes inferred from the nuclear SSU and partial LSU rDNA.

Previous studies using the nuclear SSU rDNA have indicated that the photosynthetic euglenoids are a monophyletic group; however, some of the genera within the photosynthetic lineage are not monophyletic. To test these results further, evolutionary relationships among the photosynthetic genera were investigated by obtaining partial LSU nuclear rDNA sequences. Taxa from each of the external clades of the SSU rDNA-based phylogeny were chosen to create a combined dataset and to compare the individual LSU and SSU rDNA datasets. Conserved areas of the aligned sequences for both the LSU and SSU rDNA were used to generate parsimony, log-det, maximum-likelihood and Bayesian trees. The SSU and LSU rDNA consistently generated the same seven terminal clades; however, the relationship among those clades varied depending on the type of analysis and the dataset used. The combined dataset generated a more robust phylogeny, but the relationships among clades still varied. The addition of the LSU rDNA dataset to the euglenophyte phylogeny supports the view that the genera Euglena, Lepocinclis and Phacus are not monophyletic and substantiates the existence of several well-supported clades. A secondary structural model for the D2 region of the LSU rDNA was proposed on the basis of compensatory base changes found in the alignment.

18S ribosomal RNA and tetrapod phylogeny.

Previous phylogenetic analyses of tetrapod 18S ribosomal RNA (rRNA) sequences support the grouping of birds with mammals, whereas other molecular data, and morphological and paleontological data favor the grouping of birds with crocodiles. The 18S rRNA gene has consequently been considered odd, serving as "definitive evidence of different genes providing significantly different estimates of phylogeny in higher organisms" (p. 156; Huelsenbeck et al., 1996, Trends Ecol. Evol. 11:152-158). Our research indicates that the previous discrepancy of phylogenetic results between the 18S rRNA gene and other genes is caused mainly by (1) the misalignment of the sequences, (2) the inappropriate use of the frequency parameters, and (3) poor sequence quality. When the sequences are aligned with the aide of the secondary structure of the 18S rRNA molecule and when the frequency parameters are estimated either from all sites or from the variable domains where substitutions have occurred, the 18S rRNA sequences no longer support the grouping of the avian species with the mammalian species.

MITOCHONDRIAL DNA SEQUENCE-BASED PHYLOGENY AND THE EVOLUTION OF VIVIPARITY IN THE SCELOPORUS SCALARIS GROUP (REPTILIA, SQUAMATA).

The lizard genus Sceloporus contains both oviparous and viviparous species. The scalaris complex is the only monophyletic group within the genus that includes both reproductive modes, thus it is particularly well suited for studies of the evolution of viviparity. Approximately 874 nucleotides of mtDNA sequence data, collected from 38 specimens, comprising 25 populations of all five recognized species within the group, were used in a phylogenetic analysis of the origin of viviparity. Viviparity appears to have evolved twice in this group: once in S. goldmani, included in a clade formed by a northern group consisting of S. scalaris, S. chaneyi, and S. goldmani, and one more time in S. bicanthalis, included in the southern group formed by S. bicanthalis and S. aeneus. An oviparous population of S. bicanthalis nested within that viviparous clade, indicates that reversal from viviparity to oviparity may be possible. Degree of sequence divergence among several S. bicanthalis individuals pertaining to a population in which both parity modes occur, was no larger between oviparous and viviparous lizards than among viviparous lizards. This suggests that this population is a single species, and it may represent a transition from oviparity to viviparity or vice-versa.