The taming of an impossible child: a standardized all-in approach to the phylogeny of Hymenoptera using public database sequences

A supermatrix phylogeny of the world's bees (Hymenoptera: Anthophila)

The increasing availability of large molecular phylogenies has provided new opportunities to study the evolution of species traits, their origins and diversification, and biogeography; yet there are limited attempts to synthesise existing phylogenetic information for major insect groups. Bees (Hymenoptera: Anthophila) are a large group of insect pollinators that have a worldwide distribution, and a wide variation in ecology, morphology, and life-history traits, including sociality. For these reasons, as well as their major economic importance as pollinators, numerous molecular phylogenetic studies of family and genus-level relationships have been published, providing an opportunity to assemble a bee 'tree-of-life'. We used publicly available genetic sequence data, including phylogenomic data, reconciled to a taxonomic database, to produce a concatenated supermatrix phylogeny for the Anthophila comprising 4,586 bee species, representing 23% of species and 82% of genera. At family, subfamily, and tribe levels, support for expected relationships was robust, but between and within some genera relationships remain uncertain. Within families, sampling of genera ranged from 67 to 100% but species coverage was lower (17-41%). Our phylogeny mostly reproduces the relationships found in recent phylogenomic studies with a few exceptions. We provide a summary of these differences and the current state of molecular data available and its gaps. We discuss the advantages and limitations of this bee supermatrix phylogeny (available online at beetreeoflife.org), which may enable new insights into long standing questions about evolutionary drivers in bees, and potentially insects more generally.

Key innovations and the diversification of Hymenoptera

The order Hymenoptera (wasps, ants, sawflies, and bees) represents one of the most diverse animal lineages, but whether specific key innovations have contributed to its diversification is still unknown. We assembled the largest time-calibrated phylogeny of Hymenoptera to date and investigated the origin and possible correlation of particular morphological and behavioral innovations with diversification in the order: the wasp waist of Apocrita; the stinger of Aculeata; parasitoidism, a specialized form of carnivory; and secondary phytophagy, a reversal to plant-feeding. Here, we show that parasitoidism has been the dominant strategy since the Late Triassic in Hymenoptera, but was not an immediate driver of diversification. Instead, transitions to secondary phytophagy (from parasitoidism) had a major influence on diversification rate in Hymenoptera. Support for the stinger and the wasp waist as key innovations remains equivocal, but these traits may have laid the anatomical and behavioral foundations for adaptations more directly associated with diversification.

The first mitogenomes of the superfamily Pamphilioidea (Hymenoptera: Symphyta): Mitogenome architecture and phylogenetic inference

The Pamphilioidea represents a small superfamily of the phytophagous suborder Symphyta (Hymenoptera). Here, nearly complete mitochondrial genomes (mitogenomes) of three pamphilioid species: Chinolyda flagellicornis (Pamphiliidae), Megalodontes spiraeae and M. cephalotes (Megalodontesidae) were newly sequenced using next generation sequencing and comparatively analysed with the previously reported symphytan mitogenomes. A positive AT skew (0.013) and a negative GC skew (-0.194) were found in pamphilioid mitogenome, and a deviation from strand asymmetry was also observed in the PCGs encoded on both strands. Several gene rearrangement events were observed in four tRNA gene clusters (WCY, IQM, ARNS1EF and TP clusters), which have not been reported from symphytan mitogenomes to date. As the most parsimonious explanation, compared with the inferred insect ancestral mitogenome architecture, the occurrence of gene rearrangements in pamphilioid mitogenomes requires totally five evolutionary steps, including four transpositions and one inversion. The predicted secondary structures of tRNAs, rrnS and rrnL genes are mostly consistent with reported hymenopteran species. Phylogenetic analyses recovered the monophyly of superfamily Pamphilioidea and indicated the relationship Tenthredinoidea + (Pamphilioidea + (Cephoidea + (Orussoidea + Apocrita))) with strong nodal supports.

Molecular Phylogeny, Diversity and Zoogeography of Net-Winged Beetles (Coleoptera: Lycidae)

We synthesize the evidence from molecular phylogenetics, extant distribution, and plate tectonics to present an insight in ancestral areas, dispersal routes and the effectiveness of geographic barriers for net-winged beetle tribes (Coleoptera: Lycidae). Samples from all zoogeographical realms were assembled and phylogenetic relationships for ~550 species and 25 tribes were inferred using nuclear rRNA and mtDNA markers. The analyses revealed well-supported clades at the rank of tribes as they have been defined using morphology, but a low support for relationships among them. Most tribes started their diversification in Southeast and East Asia or are endemic to this region. Slipinskiini and Dexorini are Afrotropical endemics and Calopterini, Eurrhacini, Thonalmini, and Leptolycini remained isolated in South America and the Caribbean after their separation from northern continents. Lycini, Calochromini, and Erotini support relationships between the Nearctic and eastern Palearctic faunas; Calochromini colonized the Afrotropical realm from East Asia and Metriorrhynchini Afrotropical and Oriental realms from the drifting Indian subcontinent. Most tribes occur in the Oriental and Sino-Japanese realms, the highest alpha-taxonomic diversity was identified in Malesian tropical rainforests. The turn-over at zoogeographical boundaries is discussed when only short distance over-sea colonization events were inferred. The lycid phylogeny shows that poor dispersers can be used for reconstruction of dispersal and vicariance history over a long time-span, but the current data are insufficient for reconstruction of the early phase of their diversification.

A new megaspilid wasp from Eocene Baltic amber (Hymenoptera: Ceraphronoidea), with notes on two non-ceraphronoid families: Radiophronidae and Stigmaphronidae

Ceraphronoids are some of the most commonly collected hymenopterans, yet they remain rare in the fossil record. Conostigmus talamasi Mikó and Trietsch, sp. nov. from Baltic amber represents an intermediate form between the type genus, Megaspilus, and one of the most species-rich megaspilid genera, Conostigmus. We describe the new species using 3D data collected with synchrotron-based micro-CT equipment. This non-invasive technique allows for quick data collection in unusually high resolution, revealing morphological traits that are otherwise obscured by the amber. In describing this new species, we revise the diagnostic characters for Ceraphronoidea and discuss possible reasons why minute wasps with a pterostigma are often misidentified as ceraphronoids. Based on the lack of ceraphronoid characteristics, we remove Dendrocerus dubitatus Brues, 1937, Stigmaphronidae, and Radiophronidae from Ceraphronoidea and consider them as incertae sedis. We also provide some guidance for their future classification.

Retention of Core Meiotic Genes Across Diverse Hymenoptera

The cellular mechanisms of meiosis are critical for proper gamete formation in sexual organisms. Functional studies in model organisms have identified genes essential for meiosis, yet the extent to which this core meiotic machinery is conserved across non-model systems is not fully understood. Moreover, it is unclear whether deviation from canonical modes of sexual reproduction is accompanied by modifications in the genetic components involved in meiosis. We used a robust approach to identify and catalogue meiosis genes in Hymenoptera, an insect order typically characterized by haplodiploid reproduction. Using newly available genome data, we searched for 43 genes involved in meiosis in 18 diverse hymenopterans. Seven of eight genes with roles specific to meiosis were found across a majority of surveyed species, suggesting the preservation of core meiotic machinery in haplodiploid hymenopterans. Phylogenomic analyses of the inventory of meiosis genes and the identification of shared gene duplications and losses provided support for the grouping of species within Proctotrupomorpha, Ichneumonomorpha, and Aculeata clades, along with a paraphyletic Symphyta. The conservation of meiosis genes across Hymenoptera provides a framework for studying transitions between reproductive modes in this insect group.

Construction of a Species-Level Tree of Life for the Insects and Utility in Taxonomic Profiling

Although comprehensive phylogenies have proven an invaluable tool in ecology and evolution, their construction is made increasingly challenging both by the scale and structure of publically available sequences. The distinct partition between gene-rich (genomic) and species-rich (DNA barcode) data is a feature of data that has been largely overlooked, yet presents a key obstacle to scaling supermatrix analysis. I present a phyloinformatics framework for draft construction of a species-level phylogeny of insects (Class Insecta). Matrix-building requires separately optimized pipelines for nuclear transcriptomic, mitochondrial genomic, and species-rich markers, whereas tree-building requires hierarchical inference in order to capture species-breadth while retaining deep-level resolution. The phylogeny of insects contains 49,358 species, 13,865 genera, 760 families. Deep-level splits largely reflected previous findings for sections of the tree that are data rich or unambiguous, such as inter-ordinal Endopterygota and Dictyoptera, the recently evolved and relatively homogeneous Lepidoptera, Hymenoptera, Brachycera (Diptera), and Cucujiformia (Coleoptera). However, analysis of bias, matrix construction and gene-tree variation suggests confidence in some relationships (such as in Polyneoptera) is less than has been indicated by the matrix bootstrap method. To assess the utility of the insect tree as a tool in query profiling several tree-based taxonomic assignment methods are compared. Using test data sets with existing taxonomic annotations, a tendency is observed for greater accuracy of species-level assignments where using a fixed comprehensive tree of life in contrast to methods generating smaller de novo reference trees. Described herein is a solution to the discrepancy in the way data are fit into supermatrices. The resulting tree facilitates wider studies of insect diversification and application of advanced descriptions of diversity in community studies, among other presumed applications. [Data integration; data mining; insects; phylogenomics; phyloinformatics; tree of life.].

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.

A phylogenetic study of the tribe Antirrhineae: Genome duplications and long-distance dispersals from the Old World to the New World

PREMISE OF THE STUDY

Antirrhineae is a large tribe within Plantaginaceae. Mostly concentrated in the Mediterranean Basin, the tribe members are present both in the Old World and the New World. Current Antirrhineae phylogenies have different views on taxonomic relationships, and they lack homogeneity in terms of geographic distribution and ploidy levels. This study aims to investigate the changes in the chromosome numbers along with dispersal routes as definitive characters identifying clades.

METHODS

With the use of multiple DNA regions and taxon sampling enriched with de novo sequences, we provide an extensive phylogeny for Antirrhineae. The reconstructed phylogeny was then used to investigate changes in ploidy levels and dispersal patterns in the tribe using ChromEvol and RASP, respectively.

KEY RESULTS

Antirrhineae is a monophyletic group with six highly supported clades. ChromEvol analysis suggests the ancestral haploid chromosome number for the tribe is six, and that the tribe has experienced several duplications and gain events. The Mediterranean Basin was estimated to be the origin for the tribe with four long-distance dispersals from the Old World to the New World, three of which were associated with genome duplications.

CONCLUSIONS

On an updated Antirrhineae phylogeny, we showed that the three out of four dispersals from the Old World to the New World were coupled with changes in ploidy levels. The observed patterns suggest that increases in ploidy levels may facilitate dispersing into new environments.

Consequences of Common Topological Rearrangements for Partition Trees in Phylogenomic Inference

In phylogenomic analysis the collection of trees with identical score (maximum likelihood or parsimony score) may hamper tree search algorithms. Such collections are coined phylogenetic terraces. For sparse supermatrices with a lot of missing data, the number of terraces and the number of trees on the terraces can be very large. If terraces are not taken into account, a lot of computation time might be unnecessarily spent to evaluate many trees that in fact have identical score. To save computation time during the tree search, it is worthwhile to quickly identify such cases. The score of a species tree is the sum of scores for all the so-called induced partition trees. Therefore, if the topological rearrangement applied to a species tree does not change the induced partition trees, the score of these partition trees is unchanged. Here, we provide the conditions under which the three most widely used topological rearrangements (nearest neighbor interchange, subtree pruning and regrafting, and tree bisection and reconnection) change the topologies of induced partition trees. During the tree search, these conditions allow us to quickly identify whether we can save computation time on the evaluation of newly encountered trees. We also introduce the concept of partial terraces and demonstrate that they occur more frequently than the original "full" terrace. Hence, partial terrace is the more important factor of timesaving compared to full terrace. Therefore, taking into account the above conditions and the partial terrace concept will help to speed up the tree search in phylogenomic inference.

Evolution of Cuticular Hydrocarbons in the Hymenoptera: a Meta-Analysis

Chemical communication is the oldest form of communication, spreading across all forms of life. In insects, cuticular hydrocarbons (CHC) function as chemical cues for the recognition of mates, species, and nest-mates in social insects. Although much is known about the function of individual hydrocarbons and their biosynthesis, a phylogenetic overview is lacking. Here, we review the CHC profiles of 241 species of Hymenoptera, one of the largest and most important insect orders, which includes the Symphyta (sawflies), the polyphyletic Parasitica (parasitoid wasps), and the Aculeata (wasps, bees, and ants). We investigated whether these taxonomic groups differed in the presence and absence of CHC classes and whether the sociality of a species (solitarily vs. social) had an effect on CHC profile complexity. We found that the main CHC classes (i.e., n-alkanes, alkenes, and methylalkanes) were all present early in the evolutionary history of the Hymenoptera, as evidenced by their presence in ancient Symphyta and primitive Parasitica wasps. Throughout all groups within the Hymenoptera, the more complex a CHC the fewer species that produce it, which may reflect the Occam’s razor principle that insects’ only biosynthesize the most simple compound that fulfil its needs. Surprisingly, there was no difference in the complexity of CHC profiles between social and solitary species, with some of the most complex CHC profiles belonging to the Parasitica. This profile complexity has been maintained in the ants, but some specialization in biosynthetic pathways has led to a simplification of profiles in the aculeate wasps and bees. The absence of CHC classes in some taxa or species may be due to gene silencing or down-regulation rather than gene loss, as demonstrated by sister species having highly divergent CHC profiles, and cannot be predicted by their phylogenetic history. The presence of highly complex CHC profiles prior to the vast radiation of the social Hymenoptera indicates a ‘spring-loaded’ system where the diversity of CHC needed for the complex communication systems of social insects were already present for natural selection to act upon, rather than having evolved independently. This diversity may have aided the multiple independent evolution of sociality within the Aculeata.

Target enrichment of ultraconserved elements from arthropods provides a genomic perspective on relationships among Hymenoptera

Gaining a genomic perspective on phylogeny requires the collection of data from many putatively independent loci across the genome. Among insects, an increasingly common approach to collecting this class of data involves transcriptome sequencing, because few insects have high-quality genome sequences available; assembling new genomes remains a limiting factor; the transcribed portion of the genome is a reasonable, reduced subset of the genome to target; and the data collected from transcribed portions of the genome are similar in composition to the types of data with which biologists have traditionally worked (e.g. exons). However, molecular techniques requiring RNA as a template, including transcriptome sequencing, are limited to using very high-quality source materials, which are often unavailable from a large proportion of biologically important insect samples. Recent research suggests that DNA-based target enrichment of conserved genomic elements offers another path to collecting phylogenomic data across insect taxa, provided that conserved elements are present in and can be collected from insect genomes. Here, we identify a large set (n = 1510) of ultraconserved elements (UCEs) shared among the insect order Hymenoptera. We used in silico analyses to show that these loci accurately reconstruct relationships among genome-enabled hymenoptera, and we designed a set of RNA baits (n = 2749) for enriching these loci that researchers can use with DNA templates extracted from a variety of sources. We used our UCE bait set to enrich an average of 721 UCE loci from 30 hymenopteran taxa, and we used these UCE loci to reconstruct phylogenetic relationships spanning very old (≥220 Ma) to very young (≤1 Ma) divergences among hymenopteran lineages. In contrast to a recent study addressing hymenopteran phylogeny using transcriptome data, we found ants to be sister to all remaining aculeate lineages with complete support, although this result could be explained by factors such as taxon sampling. We discuss this approach and our results in the context of elucidating the evolutionary history of one of the most diverse and speciose animal orders.

STBase: one million species trees for comparative biology

Comprehensively sampled phylogenetic trees provide the most compelling foundations for strong inferences in comparative evolutionary biology. Mismatches are common, however, between the taxa for which comparative data are available and the taxa sampled by published phylogenetic analyses. Moreover, many published phylogenies are gene trees, which cannot always be adapted immediately for species level comparisons because of discordance, gene duplication, and other confounding biological processes. A new database, STBase, lets comparative biologists quickly retrieve species level phylogenetic hypotheses in response to a query list of species names. The database consists of 1 million single- and multi-locus data sets, each with a confidence set of 1000 putative species trees, computed from GenBank sequence data for 413,000 eukaryotic taxa. Two bodies of theoretical work are leveraged to aid in the assembly of multi-locus concatenated data sets for species tree construction. First, multiply labeled gene trees are pruned to conflict-free singly-labeled species-level trees that can be combined between loci. Second, impacts of missing data in multi-locus data sets are ameliorated by assembling only decisive data sets. Data sets overlapping with the user's query are ranked using a scheme that depends on user-provided weights for tree quality and for taxonomic overlap of the tree with the query. Retrieval times are independent of the size of the database, typically a few seconds. Tree quality is assessed by a real-time evaluation of bootstrap support on just the overlapping subtree. Associated sequence alignments, tree files and metadata can be downloaded for subsequent analysis. STBase provides a tool for comparative biologists interested in exploiting the most relevant sequence data available for the taxa of interest. It may also serve as a prototype for future species tree oriented databases and as a resource for assembly of larger species phylogenies from precomputed trees.

Higher-level phylogeny of the Hymenoptera inferred from mitochondrial genomes

Higher-level hymenopteran relationships remain unresolved in both morphological and molecular analyses. In this study, we present the most comprehensive analyses of hymenopteran relationships based on 48 mitochondrial (mt) genomes. One complete and two nearly complete mt genomes representing three hymenopteran superfamilies were newly sequenced. We assessed the influence of inclusion/exclusion of 3rd codon positions, alignment approaches, partition schemes and phylogenetic approaches on topology and nodal support within the Hymenoptera. The results showed that the topologies were sensitive to the variation of dataset and analytical approach. However, some robust and highly supported relationships were recovered: the Ichneumonomorpha was monophyletic; the Trigonalyoidea+Megalyroidea and the Diaprioidea+Chalcidoidea were consistently recovered; the Cynipoidea was generally recovered as the sister group to the Diaprioidea+Chalcidoidea. In addition, the monophyletic Aculeata and Proctotrupomorpha were recovered in some analyses. Several gene rearrangements were detected in each of the three newly sequenced mt genomes. Specifically, the Ibalia leucospoides mt genome harbors a large inversion of a gene block from trnE to trnS2. Inverted, duplicated A+T rich regions were detected in the Ibalia leucospoides mt genome, which probably played an important role during the formation of the large gene block inversion via recombination.

Is it an ant or a butterfly? Convergent evolution in the mitochondrial gene order of Hymenoptera and Lepidoptera

Insect mitochondrial genomes (mtDNA) are usually double helical and circular molecules containing 37 genes that are encoded on both strands. The arrangement of the genes is not constant for all species, and produces distinct gene orders (GOs) that have proven to be diagnostic in defining clades at different taxonomic levels. In general, it is believed that distinct taxa have a very low chance of sharing identically arranged GOs. However, examples of identical, homoplastic local rearrangements occurring in distinct taxa do exist. In this study, we sequenced the complete mtDNAs of the ants Formica fusca and Myrmica scabrinodis (Formicidae, Hymenoptera) and compared their GOs with those of other Insecta. The GO of F. fusca was found to be identical to the GO of Dytrisia (the largest clade of Lepidoptera). This finding is the first documented case of an identical GO shared by distinct groups of Insecta, and it is the oldest known event of GO convergent evolution in animals. Both Hymenoptera and Lepidoptera acquired this GO early in their evolution. Using a phylogenetic approach combined with new bioinformatic tools, the chronological order of the evolutionary events that produced the diversity of the hymenopteran GOs was determined. Additionally, new local homoplastic rearrangements shared by distinct groups of insects were identified. Our study showed that local and global homoplasies affecting the insect GOs are more widespread than previously thought. Homoplastic GOs can still be useful for characterizing the various clades, provided that they are appropriately considered in a phylogenetic and taxonomic context.

Evolutionary dynamics of the mitochondrial genome in the evaniomorpha (hymenoptera)—a group with an intermediate rate of gene rearrangement

We determined the complete mitochondrial (mt) genomes of three evaniomorph species, Ceraphron sp. (Ceraphronoidea), Gasteruption sp. (Evanioidea), and Orthogonalys pulchella (Trigonalyoidea) as well as the nearly complete mt genome from another evaniomorph species, Megalyra sp. (Megalyroidea). Each of them possesses dramatic gene rearrangements, including protein-coding or rRNA genes. Gene inversions were identified in all of these mt genomes; for example, the two rRNA genes have inverted and moved into the nad2-cox1 junction in the Megalyra sp. mt genome. In addition, we found two copies of a 10-bp complementary repeat at the beginning of rrnS and at the end of trnL(2) in the Gasteruption sp. mt genome, consistent with recombination as the possible mechanism for gene inversion and long-range movement. Although each of the genomes contains a number of repeats of varying size, there was no consistent association of the size or number of repeats with the extent or type of gene rearrangement. The breakpoint distance analysis showed the Evaniomorpha has an intermediate rate of gene rearrangement. Sequence-based phylogenetic analyses of 13 protein-coding and 2 rRNA genes in 22 hymenopteran taxa recovered a paraphyletic Evaniomorpha with the Aculeata nested within it. Within the Evaniomorpha, our analyses confirmed the Trigonalyoidea + Megalyroidea as the sister group to the Aculeata and recovered a novel clade, Ceraphronoidea + Evanioidea. In contrast to previous hymenopteran phylogenetic studies, the internal relationships of the Evaniomorpha were highly supported and robust to the variation of alignment approach and phylogenetic inference approach.

A protocol for species delineation of public DNA databases, applied to the Insecta

Public DNA databases are composed of data from many different taxa, although the taxonomic annotation on sequences is not always complete, which impedes the utilization of mined data for species-level applications. There is much ongoing work on species identification and delineation based on the molecular data itself, although applying species clustering to whole databases requires consolidation of results from numerous undefined gene regions, and introduces significant obstacles in data organization and computational load. In the current paper, we demonstrate an approach for species delineation of a sequence database. All DNA sequences for the insects were obtained and processed. After filtration of duplicated data, delineation of the database into species or molecular operational taxonomic units (MOTUs) followed a three-step process in which (i) the genetic loci L are partitioned, (ii) the species S are delineated within each locus, then (iii) species units are matched across loci to form the matrix L × S, a set of global (multilocus) species units. Partitioning the database into a set of homologous gene fragments was achieved by Markov clustering using edge weights calculated from the amount of overlap between pairs of sequences, then delineation of species units and assignment of species names were performed for the set of genes necessary to capture most of the species diversity. The complexity of computing pairwise similarities for species clustering was substantial at the cytochrome oxidase subunit I locus in particular, but made feasible through the development of software that performs pairwise alignments within the taxonomic framework, while accounting for the different ranks at which sequences are labeled with taxonomic information. Over 24 different homologs, the unidentified sequences numbered approximately 194,000, containing 41,525 species IDs (98.7% of all found in the insect database), and were grouped into 59,173 single-locus MOTUs by hierarchical clustering under parameters optimized independently for each locus. Species units from different loci were matched using a multipartite matching algorithm to form multilocus species units with minimal incongruence between loci. After matching, the insect database as represented by these 24 loci was found to be composed of 78,091 species units in total. 38,574 of these units contained only species labeled data, 34,891 contained only unlabeled data, leaving 4,626 units composed both of labeled and unlabeled sequences. In addition to giving estimates of species diversity of sequence repositories, the protocol developed here will facilitate species-level applications of modern-day sequence data sets. In particular, the L × S matrix represents a post-taxonomic framework that can be used for species-level organization of metagenomic data, and incorporation of these methods into phylogenetic pipelines will yield matrices more representative of species diversity.

Building a phylogenomic pipeline for the eukaryotic tree of life - addressing deep phylogenies with genome-scale data

Background Understanding the evolutionary relationships of all eukaryotes on Earth remains a paramount goal of modern biology, yet analyzing homologous sequences across 1.8 billion years of eukaryotic evolution is challenging. Many existing tools for identifying gene orthologs are inadequate when working with heterogeneous rates of evolution and endosymbiotic/lateral gene transfer. Moreover, genomic-scale sequencing, which was once the domain of large sequencing centers, has advanced to the point where small laboratories can now generate the data needed for phylogenomic studies. This has opened the door for increased taxonomic sampling as individual research groups have the ability to conduct genome-scale projects on their favorite non-model organism. Results Here we present some of the tools developed, and insights gained, as we created a pipeline that combines data-mining from public databases and our own transcriptome data to study the eukaryotic tree of life. The first steps of a phylogenomic pipeline involve choosing taxa and loci, and making decisions about how to handle alleles, paralogs and non-overlapping sequences. Next, orthologs are aligned for analyses including gene tree reconstruction and concatenation for supermatrix approaches. To build our pipeline, we created scripts written in Python that integrate third-party tools with custom methods. As a test case, we present the placement of five amoebae on the eukaryotic tree of life based on analyses of transcriptome data. Our scripts available on GitHUb and may be used as-is for automated analyses of large scale phylogenomics, or adapted for use in other types of studies. Conclusion Analyses on the scale of all eukaryotes present challenges not necessarily found in studies of more closely related organisms. Our approach will be of relevance to others for whom existing third-party tools fail to fully answer desired phylogenetic questions.

Conserved class of queen pheromones stops social insect workers from reproducing

A major evolutionary transition to eusociality with reproductive division of labor between queens and workers has arisen independently at least 10 times in the ants, bees, and wasps. Pheromones produced by queens are thought to play a key role in regulating this complex social system, but their evolutionary history remains unknown. Here, we identify the first sterility-inducing queen pheromones in a wasp, bumblebee, and desert ant and synthesize existing data on compounds that characterize female fecundity in 64 species of social insects. Our results show that queen pheromones are strikingly conserved across at least three independent origins of eusociality, with wasps, ants, and some bees all appearing to use nonvolatile, saturated hydrocarbons to advertise fecundity and/or suppress worker reproduction. These results suggest that queen pheromones evolved from conserved signals of solitary ancestors.

The hymenopteran tree of life: evidence from protein-coding genes and objectively aligned ribosomal data

Previous molecular analyses of higher hymenopteran relationships have largely been based on subjectively aligned ribosomal sequences (18S and 28S). Here, we reanalyze the 18S and 28S data (unaligned about 4.4 kb) using an objective and a semi-objective alignment approach, based on MAFFT and BAli-Phy, respectively. Furthermore, we present the first analyses of a substantial protein-coding data set (4.6 kb from one mitochondrial and four nuclear genes). Our results indicate that previous studies may have suffered from inflated support values due to subjective alignment of the ribosomal sequences, but apparently not from significant biases. The protein data provide independent confirmation of several earlier results, including the monophyly of non-xyelid hymenopterans, Pamphilioidea + Unicalcarida, Unicalcarida, Vespina, Apocrita, Proctotrupomorpha and core Proctotrupomorpha. The protein data confirm that Aculeata are nested within a paraphyletic Evaniomorpha, but cast doubt on the monophyly of Evanioidea. Combining the available morphological, ribosomal and protein-coding data, we examine the total-evidence signal as well as congruence and conflict among the three data sources. Despite an emerging consensus on many higher-level hymenopteran relationships, several problems remain unresolved or contentious, including rooting of the hymenopteran tree, relationships of the woodwasps, placement of Stephanoidea and Ceraphronoidea, and the sister group of Aculeata.

Transcriptome analyses of primitively eusocial wasps reveal novel insights into the evolution of sociality and the origin of alternative phenotypes

BACKGROUND

Understanding how alternative phenotypes arise from the same genome is a major challenge in modern biology. Eusociality in insects requires the evolution of two alternative phenotypes - workers, who sacrifice personal reproduction, and queens, who realize that reproduction. Extensive work on honeybees and ants has revealed the molecular basis of derived queen and worker phenotypes in highly eusocial lineages, but we lack equivalent deep-level analyses of wasps and of primitively eusocial species, the latter of which can reveal how phenotypic decoupling first occurs in the early stages of eusocial evolution.

RESULTS

We sequenced 20 Gbp of transcriptomes derived from brains of different behavioral castes of the primitively eusocial tropical paper wasp Polistes canadensis. Surprisingly, 75% of the 2,442 genes differentially expressed between phenotypes were novel, having no significant homology with described sequences. Moreover, 90% of these novel genes were significantly upregulated in workers relative to queens. Differential expression of novel genes in the early stages of sociality may be important in facilitating the evolution of worker behavioral complexity in eusocial evolution. We also found surprisingly low correlation in the identity and direction of expression of differentially expressed genes across similar phenotypes in different social lineages, supporting the idea that social evolution in different lineages requires substantial de novo rewiring of molecular pathways.

CONCLUSIONS

These genomic resources for aculeate wasps and first transcriptome-wide insights into the origin of castes bring us closer to a more general understanding of eusocial evolution and how phenotypic diversity arises from the same genome.

Oligonucleotide primers for targeted amplification of single-copy nuclear genes in apocritan Hymenoptera

Background

Published nucleotide sequence data from the mega-diverse insect order Hymenoptera (sawflies, bees, wasps, and ants) are taxonomically scattered and still inadequate for reconstructing a well-supported phylogenetic tree for the order. The analysis of comprehensive multiple gene data sets obtained via targeted PCR could provide a cost-effective solution to this problem. However, oligonucleotide primers for PCR amplification of nuclear genes across a wide range of hymenopteran species are still scarce.

Findings

Here we present a suite of degenerate oligonucleotide primer pairs for PCR amplification of 154 single-copy nuclear protein-coding genes from Hymenoptera. These primers were inferred from genome sequence data from nine Hymenoptera (seven species of ants, the honeybee, and the parasitoid wasp Nasonia vitripennis). We empirically tested a randomly chosen subset of these primer pairs for amplifying target genes from six Hymenoptera, representing the families Chrysididae, Crabronidae, Gasteruptiidae, Leucospidae, Pompilidae, and Stephanidae. Based on our results, we estimate that these primers are suitable for studying a large number of nuclear genes across a wide range of apocritan Hymenoptera (i.e., all hymenopterans with a wasp-waist) and of aculeate Hymenoptera in particular (i.e., apocritan wasps with stingers).

Conclusions

The amplified nucleotide sequences are (a) with high probability from single-copy genes, (b) easily generated at low financial costs, especially when compared to phylogenomic approaches, (c) easily sequenced by means of an additionally provided set of sequencing primers, and (d) suitable to address a wide range of phylogenetic questions and to aid rapid species identification via barcoding, as many amplicons contain both exonic and fast-evolving intronic nucleotides.

Impact of the partitioning scheme on divergence times inferred from Mammalian genomic data sets

Data partitioning has long been regarded as an important parameter for phylogenetic inference. The division of heterogeneous multigene data sets into partitions with similar substitution patterns is known to increase the performance of probabilistic phylogenetic methods. However, the effect of the partitioning scheme on divergence time estimates has generally been ignored. To investigate the impact of data partitioning on the estimation of divergence times, we have constructed two genomic data sets. The first one with 15 nuclear genes comprising 50,928 bp were selected from the OrthoMam database; the second set was composed of complete mitochondrial genomes. We studied two partitioning schemes: concatenated supermatrices and partitioned gene analysis. We have also measured the impact of taxonomic sampling on the estimates. After drawing divergence time inferences using the uncorrelated relaxed clock in BEAST, we have compared the age estimates between the partitioning schemes. Our results show that, in general, both schemes resulted in similar chronological estimates, however the concatenated data sets were more efficient than the partitioned ones in attaining suitable effective sample sizes.

Accelerated evolution of mitochondrial but not nuclear genomes of Hymenoptera: new evidence from crabronid wasps

Mitochondrial genes in animals are especially useful as molecular markers for the reconstruction of phylogenies among closely related taxa, due to the generally high substitution rates. Several insect orders, notably Hymenoptera and Phthiraptera, show exceptionally high rates of mitochondrial molecular evolution, which has been attributed to the parasitic lifestyle of current or ancestral members of these taxa. Parasitism has been hypothesized to entail frequent population bottlenecks that increase rates of molecular evolution by reducing the efficiency of purifying selection. This effect should result in elevated substitution rates of both nuclear and mitochondrial genes, but to date no extensive comparative study has tested this hypothesis in insects. Here we report the mitochondrial genome of a crabronid wasp, the European beewolf (Philanthus triangulum, Hymenoptera, Crabronidae), and we use it to compare evolutionary rates among the four largest holometabolous insect orders (Coleoptera, Diptera, Hymenoptera, Lepidoptera) based on phylogenies reconstructed with whole mitochondrial genomes as well as four single-copy nuclear genes (18S rRNA, arginine kinase, wingless, phosphoenolpyruvate carboxykinase). The mt-genome of P. triangulum is 16,029 bp in size with a mean A+T content of 83.6%, and it encodes the 37 genes typically found in arthropod mt genomes (13 protein-coding, 22 tRNA, and two rRNA genes). Five translocations of tRNA genes were discovered relative to the putative ancestral genome arrangement in insects, and the unusual start codon TTG was predicted for cox2. Phylogenetic analyses revealed significantly longer branches leading to the apocritan Hymenoptera as well as the Orussoidea, to a lesser extent the Cephoidea, and, possibly, the Tenthredinoidea than any of the other holometabolous insect orders for all mitochondrial but none of the four nuclear genes tested. Thus, our results suggest that the ancestral parasitic lifestyle of Apocrita is unlikely to be the major cause for the elevated substitution rates observed in hymenopteran mitochondrial genomes.

Inferring the Tree of Life: chopping a phylogenomic problem down to size?

The combination of molecular sequence data and bioinformatics has revolutionized phylogenetic inference over the past decade, vastly increasing the scope of the evolutionary trees that we are able to infer. A recent paper in BMC Biology describing a new phylogenomic pipeline to help automate the inference of evolutionary trees from public sequence databases provides another important tool in our efforts to derive the Tree of Life.

See research article: http://www.biomedcentral.com/1741-7007/9/55